How to Improve The Accuracy of Equilibrium Molecular Dynamics For Computation of Thermal Conductivity?

Abstract

Equilibrium molecular dynamics (EMD) simulations through Green-Kubo formula (GKF) have been widely used in the study of thermal conductivity of various materials. However, there exist controversial simulation results which have huge discrepancies with experimental ones in literatures. In this paper, we demonstrate that the fluctuation in calculated thermal conductivity is due to the uncertainty in determination of the truncation time, which is related to the ensemble and size dependent phonon relaxation time. We thus propose a new scheme in the direct integration of heat current autocorrelation function (HCACF) and a nonzero correction in the double-exponential-fitting of HCACF to describe correctly the contribution to thermal conductivity from low frequency phonons. By using crystalline Silicon (Si) and Germanium (Ge) as examples, we demonstrate that our method can give rise to the values of thermal conductivity in an excellent agreement with experimental ones.

21 pages, 3 figures

References in corpus (3)

• Calculation of nanowire thermal conductivity using complete phonon dispersion relations
• Carbon Nanocone: A Promising Thermal Rectifier
• Thermal rectification in carbon nanotube intramolecular junctions: Molecular dynamics calculations

Cited by in corpus (9)

• Extreme Low Thermal Conductivity in Nanoscale 3D Si Phononic Crystal with Spherical Pores
• Thermal transport in nanostructures
• Ab initio Green-Kubo Approach for the Thermal Conductivity of Solids
• Phonon Coherent Resonance and Its Effect on Thermal Transport In Core-Shell Nanowires
• Machine Learning for Predicting Thermal Transport Properties of Solids
• Accurate thermal conductivities from optimally short molecular dynamics simulations
• Ab initio Green-Kubo simulations of heat transport in solids: Method and implementation
• Noncontact friction: Role of phonon damping and its nonuniversality
• Anomalous Temperature Induced Transition and Convergence of Thermal Conductivity in Germanene Monolayer