Equivalence of the equilibrium and the nonequilibrium molecular dynamics methods for thermal conductivity calculations: From bulk to nanowire silicon

Abstract

Molecular dynamics simulations play an important role in studying heat transport in complex materials. The lattice thermal conductivity can be computed either using the Green-Kubo formula in equilibrium MD (EMD) simulations or using Fourier's law in nonequilibrium MD (NEMD) simulations. These two methods have not been systematically compared for materials with different dimensions and inconsistencies between them have been occasionally reported in the literature. Here we give an in-depth comparison of them in terms of heat transport in three allotropes of Si: three dimensional bulk silicon, two-dimensional silicene, and quasi-one-dimensional silicon nanowire. By multiplying the correlation time in the Green-Kubo formula with an appropriate effective group velocity, we can express the running thermal conductivity in the EMD method as a function of an effective length and directly compare it with the length-dependent thermal conductivity in the NEMD method. We find that the two methods quantitatively agree with each other for all the systems studied, firmly establishing their equivalence in computing thermal conductivity.

arXiv comment8 pages, 7 figures

References in corpus (11)

• Thermal conduction in classical low-dimensional lattices
• Evidence for Dirac Fermions in a honeycomb lattice based on silicon
• Molecular dynamics calculation of the thermal conductivity of vitreous silica
• Force and heat current formulas for many-body potentials in molecular dynamics simulation with applications to thermal conductivity calculations
• Efficient molecular dynamics simulations with many-body potentials on graphics processing units
• Thermal conductivity decomposition in two-dimensional materials: Application to graphene
• An enhanced version of the heat exchange algorithm with excellent energy conservation properties
• Divergence of the Thermal Conductivity in Uniaxially Strained Graphene
• Length dependence of thermal conductivity by approach-to-equilibrium molecular dynamics
• Bimodal grain-size scaling of thermal transport in polycrystalline graphene from large-scale molecular dynamics simulations
• Size Effects in Thermal Conduction by Phonons

Cited by in corpus (16)

• Neuroevolution machine learning potentials: Combining high accuracy and low cost in atomistic simulations and application to heat transport
• GPUMD: A package for constructing accurate machine-learned potentials and performing highly efficient atomistic simulations
• Homogeneous nonequilibrium molecular dynamics method for heat transport and spectral decomposition with many-body potentials
• Influence of Boundaries and Thermostatting on Nonequilibrium Molecular Dynamics Simulations of Heat Conduction in Solids
• Molecular dynamics simulations of heat transport using machine-learned potentials: A mini review and tutorial on GPUMD with neuroevolution potentials
• Thermal Transport in MoS${}_2$ from Molecular Dynamics using Different Empirical Potentials
• Heat transport in pristine and polycrystalline single-layer hexagonal boron nitride
• Twist-Dependent Anisotropic Thermal Conductivity in Homogeneous MoS${}_2$ Stacks
• Thermal transport properties of single-layer black phosphorous from extensive molecular dynamics simulations
• Thermal conductivity reduction in carbon nanotube by fullerene encapsulation: A molecular dynamics study
• Solving Optimization Problems through Fully Convolutional Networks: an Application to the Travelling Salesman Problem
• A minimal Tersoff potential for diamond silicon with improved descriptions of elastic and phonon transport properties
• Interpretation of apparent thermal conductivity in finite systems from equilibrium molecular dynamics simulations
• Effects of transverse geometry on the thermal conductivity of Si and Ge nanowires
• Thermal conductivity of MgO in giant planetary interior conditions predicted by deep potential
• Exactly equivalent thermal conductivity in finite systems from equilibrium and nonequilibrium molecular dynamics simulations