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 Simulation on the Out-of Plane Thermal Conductivity of Single-Crystal Carbon Thin Films

Thermal Conductivity of A New Carbon Nanotube Analog: the Diamond Nanothread

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

Molecular dynamics study on the dependence of thermal conductivity on size and strain in GaN nanofilms

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

Molecular Dynamics Study of Solid Thin-Film Thermal Conductivity

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

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

Research on the Out-of-plane thermal conductivity of nanometer silicon film

A perspective on effective medium models of thermal conductivity in (ultra)nanocrystalline diamond films

Ultrahigh Convergent Thermal Conductivity of Carbon Nanotubes from Comprehensive Atomistic Modeling

Molecular Dynamics Simulation of thermal conductivity of diamond/epoxy resin composites

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

Size Effect Of Lattice Thermal Conductivity Across Nanoscale Thin Films By Molecular Dynamics Simulations

Thermal Conductive Performance of Deposited Amorphous Carbonmaterials by Molecular Dynamics Simulation

Thermal Conductivity Of Nanocrystalline Silicon By Direct Molecular Dynamics Simulation

MD Simulations of Nanoparticle Thermal Conductivity and Nanoscale Thermal Contact Resistance

Molecular Dynamics Prediction of Thermal Conductivity of GaN Films and Wires at Realistic Length Scales

Molecular Dynamics Simulations of Thermal Conductivity of Silicon Nanotubes

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

Thermal Conductivity of Nanoscale Thin Nickel Films