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.

MOLECULAR DYNAMICS STUDY ON THERMAL CONDUCTIVITY AND DISCUSSION ON SOME RELATED TOPICS

Application of Molecular Dynamics Simulation in Thermal Problems

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

Exploring Thermal Interface Materials for Electronics Via Molecular Dynamics Simulations: A Review

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

MOLECULAR SIMULATION OF SOME THERMOPHYSICAL PROPERTIES AND PHENOMENA

MD Simulations of Nanoparticle Thermal Conductivity and Nanoscale Thermal Contact Resistance

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

Achievements and Prospects of Molecular Dynamics Simulations in Thermofluid Sciences

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

Analysis of the mechanism of enhanced heat transfer by nanofluids

Molecular dynamics simulation and study on some related topics of phase change heat transfer

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

Non-equilibrium Molecular Dynamics Study of Nanoscale Thermal Contact Resistance

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

Molecular dynamics simulations of self-diffusion coefficient and thermal conductivity of methane at low and moderate densities

Molecular dynamics study on thermal conductivity of diamond nanoparticles

Molecular Dynamics Method for Supercritical CO2 Heat Transfer: A Review

Molecular Dynamics Studies on the Thermal Conductivity of Single-Walled Carbon Nanotubes

Molecular Dynamics Study of Solid Thin-Film Thermal Conductivity

MOLECULAR DYNAMICS SIMULATION TO THE THERMAL CONDUCTIVITY OF COPPER THIN-FILM