- Tutorial
- Open Access
Bayesian Methods for the Investigation of Temperature Dependence in Conductivity
PRX Energy 5, 022001 – Published 21 May, 2026
DOI: https://doi.org/10.1103/nxgq-lmp6
Abstract
Temperature-dependent transport data, including diffusion coefficients and ionic conductivities, are routinely analyzed 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 behavior 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 analyzing temperature-dependent transport data, covering parameter estimation, model selection, and extrapolation with uncertainty propagation, with illustrative examples from molecular dynamics simulations of superionic materials.
Physics Subject Headings (PhySH)
Popular Summary
How ion transport changes with temperature is central to understanding energy materials such as battery electrolytes and fuel cell membranes. To characterize this dependence, researchers typically fit empirical models to data; this requires choosing from possible competing models and deriving meaningful uncertainties in the fitted parameters. These models can then be used to predict behavior at unmeasured temperatures, provided the uncertainties are properly propagated. This Tutorial introduces Bayesian methods as a coherent framework for these tasks. Rather than returning a single best-fit value for each model parameter, Bayesian analysis yields a full probability distribution. This allows researchers to rigorously quantify uncertainties in fitted parameters, objectively compare competing models, and propagate those uncertainties into predictions at new temperatures. Using examples from computer simulations of fast ion-conducting materials, the authors demonstrate the use of Bayesian methods to quantify the evidence for non-Arrhenius behavior and provide meaningful uncertainty estimates when extrapolating to unmeasured temperatures.
Article Text
Open Reports
The author and at least one reviewer agreed to include their comments with the published article as part ofAPS' Open Reports Efforts.
References (34)
- P. Canepa, Pushing forward simulation techniques of ion transport in ion conductors for energy materials, ACS Mater. Au 3, 75 (2023).
- J. Wang, J. Ding, O. Delaire, and G. Arya, Atomistic mechanisms underlying non-Arrhenius ion transport in superionic conductor , ACS Appl. Energy Mater. 4, 7157 (2021).
- The Nernst-Einstein relation between and introduces a factor, so the Arrhenius form applies to rather than .
- H. Vogel, Das Temperaturabhaengigkeitsgesetz der Viskositaet von Fluessigkeiten, Phys. Z. 22, 645 (1921).
- G. S. Fulcher, Analysis of recent measurements of the viscosity of glasses, J. Am. Ceram. Soc. 8, 339 (1925).
- G. Tammann and W. Hesse, Die Abhängigkeit der Viscosität von der Temperatur bie unterkühlten Flüssigkeiten, Z. Anorg. Allg. Chem. 156, 245 (1926).
- V. K. de Souza and D. J. Wales, Super-Arrhenius diffusion in anundercooled binary Lennard-Jones liquid results from a quantifiable correlation effect, Phys. Rev. Lett. 96, 057802 (2006).
- S. R. Yeandel, D. O. Scanlon, and P. Goddard, Enhanced -ion dynamics in trivalently doped lithium phosphidosilicate : A candidate material as a solid Li electrolyte, J. Mater. Chem. A 7, 3953 (2019).
- B. A. Goldmann, C. Rosenbach, H. A. Evans, B. Helm, B. Wankmiller, O. Maus, E. Suard, L. F. Nazar, M. R. Hansen, B. J. Morgan, M. S. Islam, and W. G. Zeier, Rotational stacking faults in the ionic conductor , Chem. Mater. 37, 9858 (2025).
- A. R. McCluskey, A. G. Squires, J. Dunn, S. W. Coles, and B. J. Morgan, Kinisi: Bayesian analysis of mass transport from molecular dynamics simulations, J. Open Source Softw. 9, 5984 (2024).
- Information criteria such as AIC and BIC address this problem [32, 33], but rely on assumptions such as normally distributed parameters that may not hold in practice.
- This follows from the principle of maximum entropy [34].
- A. R. McCluskey, Is there still a place for linearization in the chemistry curriculum? J. Chem. Educ. 100, 4174 (2023).
- J. Mayer, K. Khairy, and J. Howard, Drawing an elephant with four complex parameters, Am. J. Phys. 78, 648 (2010).
- R. E. Kass and A. E. Raftery, Bayes factors, J. Am. Stat. Assoc. 90, 773 (1995).
- D. S. Sivia and J. Skilling, Data Analysis: A Bayesian Tutorial (Oxford University Press, Oxford, UK, 2006), 2nd ed.
- This extrapolated distribution captures parameter uncertainty but not estimated measurement noise. For extrapolation far from the measured range, parameter uncertainty typically dominates, so this distinction is rarely important in practice.
- A. R. McCluskey, S. W. Coles, and B. J. Morgan, Accurate estimation of diffusion coefficients and their uncertainties from computer simulation, J. Chem. Theory Comput. 21, 79 (2025).
- A temperature-independent Haven ratio would uniformly scale the conductivity values without affecting the temperature dependence; a temperature-dependent ratio would introduce additional curvature, but the analysis procedure would be unchanged.
- N. Michaud-Agrawal, E. J. Denning, T. B. Woolf, and O. Beckstein, MDAnalysis: A toolkit for the analysis of molecular dynamics simulations, J. Comput. Chem. 32, 2319 (2011).
- R. Gowers, M. Linke, J. Barnoud, T. Reddy, M. Melo, S. Seyler, J. Domański, D. Dotson, S. Buchoux, I. Kenney, and O. Beckstein, in Proceedings of the 15th Python in Science Conference (Austin, US, 2016), pp. 98–105.
- A. R. McCluskey, S. W. Coles, and B. J. Morgan, model-arrhenius-1.0.0, 2025, https://github.com/scams-research/model-arrhenius.
- A. Marin-Laflèche, M. Haefele, L. Scalfi, A. Coretti, T. Dufils, G. Jeanmairet, S. Reed, A. Serva, R. Berthin, C. Bacon, S. Bonella, B. Rotenberg, P. Madden, and M. Salanne, MetalWalls: A classical molecular dynamics software dedicated to the simulation of electrochemical systems, J. Open Source Softw. 5, 2373 (2020).
- M. Burbano, D. Carlier, F. Boucher, B. J. Morgan, and M. Salanne, Sparse cyclic excitations explain the low ionic conductivity of stoichiometric , Phys. Rev. Lett. 116, 135901 (2016).
- S. Nosé, A unified formulation of the constant temperature molecular dynamics methods, J. Chem. Phys. 81, 511 (1984).
- W. G. Hoover, Canonical dynamics: Equilibrium phase-space distributions, Phys. Rev. A 31, 1695 (1985).
- G. J. Martyna, M. L. Klein, and M. Tuckerman, Nosé–Hoover chains: The canonical ensemble via continuous dynamics, J. Chem. Phys. 97, 2635 (1992).
- Samuel W. Coles, LLZO Simulation Trajectories, 2025, https://doi.org/10.5281/zenodo.17702491.
- J. Wang, J. Ding, O. Delaire, and G. Arya, Simulation Trajectories, 2025, https://doi.org/10.5281/zenodo.17910336.
- J. Skilling, in AIP Conference Proceedings (AIP, Garching, DE, 2004), Vol. 735, pp. 395–405.
- J. S Speagle, dynesty: A dynamic nested sampling package for estimating Bayesian posteriors and evidences, Mon. Not. R. Astron. Soc. 493, 3132 (2020).
- H. Akaike, A new look at the statistical model identification, IEEE Trans. Autom. Control 19, 716 (1974).
- G. Schwarz, Estimating the dimension of a model, Ann. Stat. 6, 461 (1978).
- E. T. Jaynes, Information theory and statistical mechanics, Phys. Rev. 106, 620 (1957).
