Abstract:Nano-size confinement induces many intriguing non-Fourier heat conduction phenomena, such as nonlinear temperature gradients, temperature jumps near the contacts, and size-dependent thermal conductivity. Over the past decades, these effects have been studied and interpreted by nonequilibrium molecular dynamics (NEMD) and phonon Boltzmann transport equation (BTE) simulations separately, but no theory that unifies these two methods has ever been established. In this work, we unify these methods using a quantitative mode-level comparison and demonstrate that they are equivalent for various thermostats. We show that different thermostats result in different non-Fourier thermal transport characteristics due to the different mode-level phonon excitations inside the thermostats, which explains the different size-dependent thermal conductivities calculated using different reservoirs, even though they give the same bulk thermal conductivity. Specifically, the Langevin thermostat behaves like a thermalizing boundary in phonon BTE and provides mode-level thermal-equilibrium phonon outlets, while the Nose-Hoover chain thermostat and velocity rescaling method behave like biased reservoirs, which provide a spatially uniform heat generation and mode-level nonequilibrium phonon outlets. These findings explain why different experimental measurement methods can yield different size-dependent thermal conductivity. They also indicate that the thermal conductivity of materials can be tuned for various applications by specifically designing thermostats. The unification of NEMD and phonon BTE will largely facilitate the study of thermal transport in complex systems in the future by, e.g., replacing computationally unaffordable first-principles NEMD simulations with computationally less expensive spectral BTE simulations.

Effects and correction of angular momentum non-conservation in RNEMD for calculating thermal conductivity

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

Interpretation of Apparent Thermal Conductivity in Finite Systems from Equilibrium Molecular Dynamics Simulations

Nonequilibrium Molecular Dynamics Calculation Of The Thermal Conductivity Based On An Improved Relaxation Scheme

Equivalence of the Equilibrium and the Nonequilibrium Molecular Dynamics Methods for Thermal Conductivity Calculations: from Bulk to Nanowire Silicon

A Unified Phonon Interpretation for the Non-Fourier Heat Conduction by Non-equilibrium Molecular Dynamics Simulations

Capturing the Nuclear Quantum Effects in Molecular Dynamics for Lattice Thermal Conductivity Calculations: Using Ice As Example

Unification of Nonequilibrium Molecular Dynamics and the Mode-Resolved Phonon Boltzmann Equation for Thermal Transport Simulations

How to Improve the Accuracy of Equilibrium Molecular Dynamics for Computation of Thermal Conductivity?

Modulating Thermal Conduction by the Axial Strain

A Uniform Source-And-Sink Scheme For Calculating Thermal Conductivity By Nonequilibrium Molecular Dynamics

Thermal conductivity prediction by atomistic simulation methods: Recent advances and detailed comparison

KINETIC ENERGY EFFECT OF THERMAL MASS ON MICROSCALE HEAT CONDUCTION

Theory for Thermal Wave Nonreciprocity: Angular Momentum Bias

Some Critical Issues in the Characterization of Nanoscale Thermal Conductivity by Molecular Dynamics Analysis

MD Simulations of Nanoparticle Thermal Conductivity and Nanoscale Thermal Contact Resistance

A Novel Energy Conversion Based Method for Velocity Correction in Molecular Dynamics Simulations

Heat-bath approach to anomalous thermal transport: effects of inelastic scattering

Revisit Many-body Interaction Heat Current and Thermal Conductivity Calculation in Moment Tensor Potential/LAMMPS Interface

Thermal Conductivity Calculation using Homogeneous Non-equilibrium Molecular Dynamics Simulation with Allegro

An enhanced version of the heat exchange algorithm with excellent energy conservation properties