Bayesian Methods for the Investigation of Temperature-Dependence in Conductivity
Pith reviewed 2026-05-25 07:55 UTC · model grok-4.3
The pith
Bayesian methods supply a single framework for estimating parameters in temperature-dependent conductivity models, selecting among them, and extrapolating with uncertainty.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Bayesian methods offer a coherent framework that addresses challenges in parameter estimation, model selection, and extrapolation with uncertainty propagation for temperature-dependent transport data, illustrated with examples from molecular dynamics simulations of superionic materials.
What carries the argument
Bayesian inference applied to empirical functional forms such as the Arrhenius equation, enabling joint treatment of parameter posteriors, model probabilities, and predictive distributions.
If this is right
- Fitted activation energies and prefactors come with full posterior distributions rather than point estimates.
- Model choice between candidate functional forms is decided by explicit posterior odds that incorporate data quality.
- Extrapolated conductivities at unmeasured temperatures include credible intervals that reflect both parameter and model uncertainty.
- The same workflow applies directly to diffusion coefficients and other temperature-dependent transport quantities.
Where Pith is reading between the lines
- The framework could be extended to joint analysis of conductivity and diffusion data from the same simulation trajectories.
- It would allow systematic comparison of multiple empirical models across a library of materials rather than one at a time.
- Adoption would change how simulation papers report extrapolated performance at device operating temperatures.
Load-bearing premise
Empirical models like the Arrhenius equation are suitable functional forms whose adequacy can be assessed from the available simulation data.
What would settle it
A dataset generated from a known non-Arrhenius mechanism where the Bayesian procedure still assigns high probability to the Arrhenius model would falsify the claim that the framework reliably assesses model adequacy.
Figures
read the original abstract
Temperature-dependent transport data, including diffusion coefficients and ionic conductivities, are routinely analysed by fitting empirical models such as the Arrhenius equation. These fitted models yield parameters such as the activation energy, and can be used to extrapolate to temperatures outside the measured range. Researchers frequently face challenges in this analysis: quantifying the uncertainty of fitted parameters, assessing whether the data quality is sufficient to support a particular empirical model, and using these models to predict behaviour at temperatures outside the measured range. Bayesian methods offer a coherent framework that addresses all of these challenges. This tutorial introduces the use of Bayesian methods for analysing temperature-dependent transport data, covering parameter estimation, model selection, and extrapolation with uncertainty propagation, with illustrative examples from molecular dynamics simulations of superionic materials.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript is a tutorial on applying Bayesian inference to temperature-dependent transport data (diffusion coefficients and ionic conductivities) obtained from molecular dynamics simulations of superionic materials. It demonstrates parameter estimation for empirical models such as the Arrhenius equation, model selection to assess data quality and model adequacy, and extrapolation to unmeasured temperatures with full uncertainty propagation via posterior predictive distributions.
Significance. If the tutorial examples are reproducible and correctly illustrate the claimed advantages, the work provides a practical, coherent framework that directly addresses three recurring challenges in materials-science transport analysis: parameter uncertainty, model adequacy testing, and reliable extrapolation. The explicit use of marginal likelihoods or information criteria for model selection and posterior predictive checks for extrapolation constitutes a clear pedagogical contribution that can be adopted without new theoretical machinery.
minor comments (3)
- The abstract states that the tutorial covers 'illustrative examples from molecular dynamics simulations' but does not indicate whether the underlying conductivity data, prior specifications, or Stan/JAGS/PyMC code are provided as supplementary material or a public repository; this is essential for a methods tutorial.
- In the section describing the Arrhenius model, the functional form and the definition of the likelihood (including any assumed error model on the conductivity values) should be written explicitly with equation numbers so that readers can reproduce the posterior exactly.
- The discussion of model selection would benefit from a short table comparing the marginal likelihoods (or WAIC/LOO values) across the candidate models for at least one of the MD datasets, rather than only qualitative statements.
Simulated Author's Rebuttal
We thank the referee for their positive summary of the manuscript as a tutorial on Bayesian methods for temperature-dependent transport data and for recommending minor revision. No specific major comments appear in the report.
Circularity Check
No significant circularity
full rationale
The paper is a tutorial demonstrating standard Bayesian inference (posteriors for parameters, marginal likelihoods for model selection, posterior predictive distributions for extrapolation) applied to Arrhenius-style fits on external MD-derived conductivity data. No equations, derivations, or predictions are shown that reduce by construction to quantities defined or fitted within the paper itself. The central claim follows directly from the general properties of Bayesian methods and requires no domain-specific self-referential steps, self-citation chains, or ansatzes smuggled via prior work by the same authors. The illustrative examples serve only to demonstrate application and do not form load-bearing premises.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
Ohno, S.et al.Materials design of ionic conductors for solid state batteries.Prog. in Energy2, 022001 (2020)
work page 2020
-
[2]
R.et al.Neutron studies of Na-ion battery materials.J
Shah, A. R.et al.Neutron studies of Na-ion battery materials.J. Phys.: Mater.4, 042008 (2021)
work page 2021
-
[3]
Energy Mater.5, 14119–14126 (2022)
Kumar, S.et al.Li-Ion Diffusion Correlations in LiAlGeO4: Quasielastic Neutron Scattering and Ab Ini- tio Simulation.ACS Appl. Energy Mater.5, 14119–14126 (2022)
work page 2022
-
[4]
Eames, C.et al.Ionic transport in hybrid lead iodide perovskite solar cells.Nat. Commun.6, 7497 (2015)
work page 2015
-
[5]
Forrester, F. N., Quirk, J. A., Famprikis, T. & Daw- son, J. A. Disentangling Cation and Anion Dynamics in Li3PS4 Solid Electrolytes.Chem. Mater.34, 10561– 10571 (2022)
work page 2022
-
[6]
Krenzer, G.et al.Nature of the Superionic Phase Tran- sition of Lithium Nitride from Machine Learning Force Fields.Chem. Mater.35, 6133–6140 (2023)
work page 2023
-
[7]
Usler, A. L., Kemp, D., Bonkowski, A. & De Souza, R. A. A general expression for the statistical error in a diffusion coefficient obtained from a solid-state ¡scp¿molecular-dynamics¡/scp¿ simulation.J. Comput. Chem.44, 1347–1359 (2023)
work page 2023
-
[8]
Pushing Forward Simulation Techniques of Ion Transport in Ion Conductors for Energy Materials
Canepa, P. Pushing Forward Simulation Techniques of Ion Transport in Ion Conductors for Energy Materials. ACS Mater. Au3, 75–82 (2022)
work page 2022
- [9]
-
[10]
Das Temperaturabhaengigkeitsgesetz der Viskositaet von Fluessigkeiten.Phys
Vogel, H. Das Temperaturabhaengigkeitsgesetz der Viskositaet von Fluessigkeiten.Phys. Z.22, 645 (1921)
work page 1921
-
[11]
Fulcher, G. S. Analysis of Recent Measurements of the Viscosity of Glasses.J. Am. Ceram. Soc.8, 339–355 (1925)
work page 1925
-
[12]
Tammann, G. & Hesse, W. Die Abh¨ angigkeit der Viscosit¨ at von der Temperatur bie unterk¨ uhlten Fl¨ ussigkeiten.Z. Anorg. Allg. Chem.156, 245–257 (1926)
work page 1926
-
[13]
de Souza, V. K. & Wales, D. J. Super-Arrhenius Diffusion in an Undercooled Binary Lennard-Jones Liquid Results from a Quantifiable Correlation Effect.Phys. Rev. Lett. 96, 057802 (2006)
work page 2006
-
[14]
Yeandel, S. R., Scanlon, D. O. & Goddard, P. Enhanced Li-ion dynamics in trivalently doped lithium phosphi- dosilicate Li 2SiP2: a candidate material as a solid Li electrolyte.J. Mater. Chem. A7, 3953–3961 (2019)
work page 2019
-
[15]
A.et al.Rotational stacking faults in the ionic conductor Li 3ScCl6.Chem
Goldmann, B. A.et al.Rotational stacking faults in the ionic conductor Li 3ScCl6.Chem. Mater.(2025)
work page 2025
-
[16]
McCluskey, A. R., Squires, A. G., Dunn, J., Coles, S. W. & Morgan, B. J. kinisi: Bayesian analysis of mass transport from molecular dynamics simulations.J. Open Source Softw.9, 5984 (2024)
work page 2024
-
[17]
A new look at the statistical model identifi- cation.IEEE Trans
Akaike, H. A new look at the statistical model identifi- cation.IEEE Trans. Autom. Contr.19, 716–723 (1974)
work page 1974
-
[18]
Estimating the Dimension of a Model.Ann
Schwarz, G. Estimating the Dimension of a Model.Ann. Stat.6, 461–464 (1978)
work page 1978
-
[19]
McCluskey, A. R. Is There Still a Place for Linearization in the Chemistry Curriculum?J. Chem. Educ.100, 4174–4176 (2023)
work page 2023
-
[20]
Mayer, J., Khairy, K. & Howard, J. Drawing an elephant with four complex parameters.Am. J. Phys.78, 648–649 (2010)
work page 2010
-
[21]
Kass, R. E. & Raftery, A. E. Bayes Factors.J. Am. Stat. Assoc.90, 773–795 (1995)
work page 1995
-
[22]
Sivia, D. S. & Skilling, J.Data Analysis: A Bayesian Tutorial(Oxford University Press, Oxford, UK, 2006), 2nd edn
work page 2006
-
[23]
Mean-squared displacements exhibit ballistic behaviour at short times and diffusive behaviour at longer times [26, 27]. Because diffusion coefficients are estimated by linear regression of MSD against time, only data from the dif- fusive regime, where the MSD is approximately linear in time, is used for fitting
-
[24]
This extrapolated distribution captures parameter uncer- tainty but not estimated measurement noise. For extrap- olation far from the measured range, parameter uncer- tainty typically dominates, so this distinction is rarely important in practice
-
[25]
McCluskey, A. R., Coles, S. W. & Morgan, B. J. Accu- rate Estimation of Diffusion Coefficients and their Un- 9 certainties from Computer Simulation.J. Chem. Theory Comput.21, 79–87 (2024)
work page 2024
-
[26]
He, X., Zhu, Y., Epstein, A. & Mo, Y. Statistical vari- ances of diffusional properties from ab initio molecular dynamics simulations.npj Comput. Mater.4, 18 (2018)
work page 2018
-
[27]
Bonkowski, A. & De Souza, R. A. From setup to analy- sis: A compact guide to performing Molecular Dynamics simulations of ion transport in solids.Solid State Ionics 429, 116967 (2025)
work page 2025
-
[28]
McCluskey, A. R., Coles, S. W. & Morgan, B. J. model-arrhenius-0.0.1. https://github.com/scams- research/model-arrhenius (2025)
work page 2025
-
[29]
Open Source Softw.5, 2373 (2020)
Marin-Lafl` eche, A.et al.MetalWalls: A classical molec- ular dynamics software dedicated to the simulation of electrochemical systems.J. Open Source Softw.5, 2373 (2020)
work page 2020
-
[30]
Burbano, M., Carlier, D., Boucher, F., Morgan, B. J. & Salanne, M. Sparse Cyclic Excitations Explain the Low Ionic Conductivity of Stoichiometric Li7La3Xr2O12. Phys. Rev. Lett.116, 135901 (2016)
work page 2016
-
[31]
A unified formulation of the constant temper- ature molecular dynamics methods.J
Nos´ e, S. A unified formulation of the constant temper- ature molecular dynamics methods.J. Chem. Phys.81, 511–519 (1984)
work page 1984
-
[32]
Hoover, W. G. Canonical dynamics: Equilibrium phase- space distributions.Phys. Rev. A31, 1695–1697 (1985)
work page 1985
-
[33]
Martyna, G. J., Klein, M. L. & Tuckerman, M. Nos´ e–Hoover chains: The canonical ensemble via contin- uous dynamics.J. Chem. Phys.97, 2635–2643 (1992)
work page 1992
-
[34]
Coles, S. W. LLZO Simulation Trajectories. https://doi.org/10.5281/zenodo.17702491 (2025)
-
[35]
Michaud-Agrawal, N., Denning, E. J., Woolf, T. B. & Beckstein, O. MDAnalysis: A toolkit for the analysis of molecular dynamics simulations.J. Comput. Chem.32, 2319–2327 (2011)
work page 2011
-
[36]
InProceedings of the 15th Python in Science Conference, 98–105 (2016)
Gowers, R.et al.MDAnalysis: A Python Package for the Rapid Analysis of Molecular Dynamics Simulations. InProceedings of the 15th Python in Science Conference, 98–105 (2016)
work page 2016
-
[37]
Wang, J., Ding, J., Delaire, O. & Arya, G. AgCrSe2 Simulation Trajectories. https://doi.org/10.5281/zenodo.17910336 (2025)
-
[38]
Skilling, J. Nested Sampling. InAIP Conference Pro- ceedings, vol. 735, 395–405 (AIP, 2004)
work page 2004
-
[39]
Bayesian Methods for the Investigation of Temperature-Dependence in Conductivity
Speagle, J. S. dynesty: a dynamic nested sampling package for estimating bayesian posteriors and evidences. Mon. Not. R. Astron. Soc.493, 3132–3158 (2020). Supplemental Material for “Bayesian Methods for the Investigation of Temperature-Dependence in Conductivity” Andrew R. McCluskey, 1, 2,∗ Samuel W. Coles, 3, 4 and Benjamin J. Morgan 5, 6,† 1Centre for ...
work page 2020
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.