Influence of Boundaries and Thermostatting on Nonequilibrium Molecular Dynamics Simulations of Heat Conduction in Solids

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 to 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. Our findings are important to obtain correct results from molecular dynamics simulations of nanoscale heat transport as the accuracy of the interatomic potentials is rapidly improving.

13 pages, 12 figures

References in corpus (36)

• Canonical sampling through velocity-rescaling
• Thermal conduction in classical low-dimensional lattices
• Phononics: Manipulating heat flow with electronic analogs and beyond
• Optimized Tersoff and Brenner empirical potential parameters for lattice dynamics and phonon thermal transport in carbon nanotubes and graphene
• Thermal diode: Rectification of heat flux
• Length-dependent thermal conductivity in suspended single-layer graphene
• Heat Transport in low-dimensional systems
• Thermal Conductivity and Thermal Rectification in Graphene Nanoribbons: a Molecular Dynamics Study
• Controlling the energy flow in nonlinear lattices: a model for a thermal rectifier
• Quantum thermal transport in nanostructures
• Molecular dynamics calculation of the thermal conductivity of vitreous silica
• Accurate sampling using Langevin dynamics
• Role of anharmonic phonon scattering in the spectrally decomposed thermal conductance at planar interfaces
• 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
• Frequency-dependent phonon mean free path in carbon nanotubes from non-equilibrium molecular dynamics
• Homogeneous nonequilibrium molecular dynamics method for heat transport and spectral decomposition with many-body potentials
• Heat transport in low dimensions: introduction and phenomenology
• Thermal conductivity decomposition in two-dimensional materials: Application to graphene
• Effect of Grain Boundaries on Thermal Transport in Graphene
• Ab initio Green-Kubo Approach for the Thermal Conductivity of Solids
• Microscopic theory and quantum simulation of atomic heat transport
• A revisit to phonon-phonon scattering in single-layer graphene
• Molecular Dynamics Simulations of Heat Conduction in Nanostructures: Effect of Heat Bath
• Equivalence of the equilibrium and the nonequilibrium molecular dynamics methods for thermal conductivity calculations: From bulk to nanowire silicon
• Accelerated molecular dynamics force evaluation on graphics processing units for thermal conductivity calculations
• Thermal rectification in thickness-asymmetric graphene nanoribbons
• Thermal Transport in MoS${}_2$ from Molecular Dynamics using Different Empirical Potentials
• Boltzmann Transport in Nanostructures as a Friction Effect
• Non-Ergodicity of Nose-Hoover chain thermostat in computationally achievable time
• Heat transport in pristine and polycrystalline single-layer hexagonal boron nitride
• Thermal transport properties of single-layer black phosphorous from extensive molecular dynamics simulations
• Ingredients for an efficient thermal diode
• Size Effects in Thermal Conduction by Phonons
• An efficient thermal diode with ballistic spacer
• Thermal Transport Across Graphene Step Junctions

Cited by in corpus (30)

• GPUMD: A package for constructing accurate machine-learned potentials and performing highly efficient atomistic simulations
• Van der Waals stacked multilayer in-plane graphene/hexagonal boron nitride heterostructure: its interfacial thermal transport properties
• Molecular dynamics simulations of heat transport using machine-learned potentials: A mini review and tutorial on GPUMD with neuroevolution potentials
• Quantum-corrected thickness-dependent thermal conductivity in amorphous silicon predicted by machine-learning molecular dynamics simulations
• Machine learned Force-Fields for an ab-initio Quality Description of Metal-Organic Frameworks
• Anisotropic and high thermal conductivity in monolayer quasi-hexagonal fullerene: A comparative study against bulk phase fullerene
• Sub-micrometer phonon mean free paths in metal-organic frameworks revealed by machine-learning molecular dynamics simulations
• Non-Fourier heat transport in nanosystems
• Engineering thermal transport across layered graphene-MoS2 superlattices
• Ultrahigh convergent thermal conductivity of carbon nanotubes from comprehensive atomistic modeling
• Thermal conductivity reduction in carbon nanotube by fullerene encapsulation: A molecular dynamics study
• Moiré-Driven Interfacial Thermal Transport in Twisted Transition Metal Dichalcogenides
• Accurate estimation of interfacial thermal conductance between silicon and diamond enabled by a machine learning interatomic potential
• Heat transport in carbon nanotubes: Length dependence of phononic conductivity from the Boltzmann transport equation and molecular dynamics
• Thermal rectification and interface thermal resistance in hybrid pillared-graphene and graphene: A molecular dynamics approach
• Heat transport across graphene/hexagonal-BN tilted grain boundaries from phase-field crystal model and molecular dynamics simulations
• Interpretation of apparent thermal conductivity in finite systems from equilibrium molecular dynamics simulations
• Probing the ideal limit of interfacial thermal conductance in two-dimensional van der Waals heterostructures
• Effects of transverse geometry on the thermal conductivity of Si and Ge nanowires
• Simulations of Heat Transport in Single-Molecule Junctions: Investigations of the Thermal Diode Effect
• Exactly equivalent thermal conductivity in finite systems from equilibrium and nonequilibrium molecular dynamics simulations
• Ultrahigh interfacial thermal conductance for cooling gallium oxide electronics using cubic boron arsenide
• Heat Transport with a Twist
• Competing mechanisms govern the thermal rectification behavior in semi-stochastic polycrystalline graphene with graded grain-size distribution
• Anisotropic and isotropic elasticity and thermal transport in monolayer C${}_{24}$ networks from machine-learning molecular dynamics
• Non-local thermal transport modeling using the thermal distributor
• Substrate-Dependence of Monolayer MoS${}_2$ Thermal Conductivity and Thermal Boundary Conductance
• Uniaxial stress tuning of interfacial thermal conductance in cubic BAs/4H-SiC heterostructures
• Convergence of Nonequilibrium Langevin Dynamics for Planar Flows
• Effects of ballistic transport on the thermal resistance and temperature profile in nanowires