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.

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

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

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

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

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

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

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

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

Numerical Experiment Of Thermal Conductivity In Two-Dimensional Yukawa Liquids

Exactly Equivalent Thermal Conductivity in Finite Systems from Equilibrium and Nonequilibrium Molecular Dynamics Simulations

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

Heat Current Calculation in Equilibrium Molecular Dynamics Simulations of Thermal Conductivity

Thermal Conductivity Calculation of Complex (dusty) Plasmas

Modulating Thermal Conduction by the Axial Strain

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

The Calculation of Thermal Conductivities by Three Dimensional Direct Simulation Monte Carlo Method.

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

Molecular Dynamics Simulations of Relaxation and Heat Conduction in One-Dimensional FPU Β Lattice

Thermal transport and non-newtonian behaviors of 3D complex liquids using molecular simulations

Molecular Dynamics Simulations Of Thermal Conductivity In 2d Complex Yukawa Liquids