Accurate thermal conductivities from optimally short molecular dynamics simulations

Abstract

The evaluation of transport coefficients in extended systems, such as thermal conductivity or shear viscosity, is known to require impractically long simulations, thus calling for a paradigm shift that would allow to deploy state-of-the-art quantum simulation methods. We introduce a new method to compute these coefficients from optimally short molecular dynamics simulations, based on the Green-Kubo theory of linear response and the cepstral analysis of time series. Information from the full sample power spectrum of the relevant current for a single and relatively short trajectory is leveraged to evaluate and optimally reduce the noise affecting its zero-frequency value, whose expectation is proportional to the corresponding conductivity. Our method is unbiased and consistent, in that both the resulting bias and statistical error can be made arbitrarily small in the long-time limit. A simple data-analysis protocol is proposed and validated with the calculation of thermal conductivities in the paradigmatic cases of elemental and molecular fluids (liquid Ar and H${}_2$O) and of crystalline and glassy solids (MgO and a-SiO${}_2$). We find that simulation times of one to a few hundred picoseconds are sufficient in these systems to achieve an accuracy of the order of 10% on the estimated thermal conductivities.

12 pages, 5 figures; published version

References in corpus (2)

• Ab initio Green-Kubo Approach for the Thermal Conductivity of Solids
• Gauge invariance of thermal transport coefficients

Cited by in corpus (27)

• Modeling heat transport in crystals and glasses from a unified lattice-dynamical approach
• Heat transport in liquid water from first-principles and deep-neural-network simulations
• Computational prediction of lattice thermal conductivity -- a comparison of molecular dynamics and Boltzmann transport approaches
• Viscosity in water from first-principles and deep-neural-network simulations
• Topological quantisation and gauge invariance of charge transport in liquid insulators
• Theory and numerical simulation of heat transport in multi-component systems
• Computing the heat conductivity of fluids from density fluctuations
• Temperature- and vacancy-concentration-dependence of heat transport in Li${}_3$ClO from multi-method numerical simulations
• Hydrodynamic finite-size scaling of the thermal conductivity in glasses
• Unsupervised landmark analysis for jump detection in molecular dynamics simulations
• SporTran: a code to estimate transport coefficients from the cepstral analysis of (multivariate) current time series
• Ab initio Green-Kubo simulations of heat transport in solids: Method and implementation
• Invariance principles in the theory and computation of transport coefficients
• Heat transport in insulators from ab initio Green-Kubo theory
• Mode Localization and Suppressed Heat Transport in Amorphous Alloys
• Sampling mobility profiles of confined fluids with equilibrium molecular dynamics simulations
• Oxidation states, Thouless' pumps, and non-trivial ionic transport in non-stoichiometric electrolytes
• Reduced variance analysis of molecular dynamics simulations by linear combination of estimators
• Thermal conductivity of Li${}_3$PS${}_4$ solid electrolytes with ab initio accuracy
• Transport properties in liquids from first principles: the case of liquid water and liquid Argon
• Self-interaction and transport of solvated electrons in molten salts
• Heat conductivity from energy-density fluctuations
• Seebeck coefficient of liquid water from equilibrium molecular dynamics
• STable AutoCorrelation Integral Estimator (STACIE): Robust and accurate transport properties from molecular dynamics simulations
• Investigating finite-size effects in molecular dynamics simulations of ion diffusion, heat transport, and thermal motion in superionic materials
• Reliable Viscosity Calculation from High-Pressure Equilibrium Molecular Dynamics: Case Study of 2,2,4-Trimethylhexane
• Topology, oxidation states, and charge transport in ionic conductors