Influence of thermostatting on nonequilibrium molecular dynamics simulations of heat conduction in solids
Top Cited Papers
- 16 December 2019
- journal article
- research article
- Published by AIP Publishing in The Journal of Chemical Physics
- Vol. 151 (23), 234105
- https://doi.org/10.1063/1.5132543
Abstract
Nonequilibrium molecular dynamics (NEMD) has been extensively used to study thermal transport at various length scales in many materials. In this method, two local thermostats at different temperatures are used to generate a nonequilibrium steady state with a constant heat flux. Conventionally, the thermal conductivity of a finite system is calculated as the ratio between the heat flux and the temperature gradient extracted from the linear part of the temperature profile away from the local thermostats. Here, we show that, with a proper choice of the thermostat, the nonlinear part of the temperature profile should actually not be excluded in thermal transport calculations. We compare NEMD results against those from the atomistic Green’s function method in the ballistic regime and those from the homogeneous nonequilibrium molecular dynamics method in the ballistic-to-diffusive regime. These comparisons suggest that in all the transport regimes, one should directly calculate the thermal conductance from the temperature difference between the heat source and sink and, if needed, convert it into the thermal conductivity by multiplying it with the system length. Furthermore, we find that the Langevin thermostat outperforms the Nosé-Hoover (chain) thermostat in NEMD simulations because of its stochastic and local nature. We show that this is particularly important for studying asymmetric carbon-based nanostructures, for which the Nosé-Hoover thermostat can produce artifacts leading to unphysical thermal rectification.Keywords
Funding Information
- Jiangsu Provincial Natural Science Foundation (BK20160308)
- National Natural Science Foundation of China (11804242)
- Academy of Finland (312298)
- National Natural Science Foundation of China (11974059, 11502217, 51676121)
This publication has 75 references indexed in Scilit:
- A gentler approach to RNEMD: Nonisotropic velocity scaling for computing thermal conductivity and shear viscosityThe Journal of Chemical Physics, 2010
- Theoretical Methods for Calculating the Lattice Thermal Conductivity of MineralsReviews in Mineralogy and Geochemistry, 2010
- Quantum thermal transport in nanostructuresZeitschrift für Physik B Condensed Matter, 2008
- Canonical sampling through velocity rescalingThe Journal of Chemical Physics, 2007
- Solid-State Thermal RectifierScience, 2006
- Phonon transport in nanowires coated with an amorphous material: An atomistic Green’s function approachPhysical Review B, 2003
- Comparison of atomic-level simulation methods for computing thermal conductivityPhysical Review B, 2002
- Explicit reversible integrators for extended systems dynamicsMolecular Physics, 1996
- Lattice-dynamical calculation of the Kapitza resistance between fcc latticesPhysical Review B, 1989
- Molecular dynamics with coupling to an external bathThe Journal of Chemical Physics, 1984