Abstract:The design of new applications, especially those based on heterogeneous integration, must rely on detailed knowledge of material properties, such as thermal conductivity (TC). To this end, multiple methods have been developed to study TC as a function of vibrational frequency. Here, we compare three spectral TC methods based on velocity decomposition in homogenous molecular dynamics simulations: Green-Kubo modal analysis (GKMA), the spectral heat current (SHC) method, and a method we propose called homogeneous nonequilibrium modal analysis (HNEMA). First, we derive a convenient per-atom virial expression for systems described by general many-body potentials, enabling compact representations of the heat current, each velocity decomposition method, and other related quantities. Next, we evaluate each method by calculating the spectral TC for carbon nanotubes, graphene, and silicon. We show that each method qualitatively agrees except at optical phonon frequencies, where a combination of mismatched eigenvectors and a large density of states produces artificial TC peaks for modal analysis methods. Our calculations also show that the HNEMA and SHC methods converge much faster than the GKMA method, with the SHC method being the most computationally efficient. Finally, we demonstrate that our single-GPU modal analysis implementation in GPUMD (Graphics Processing Units Molecular Dynamics) is over 1000 times faster than the existing LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) implementation on one CPU.

Homogeneous Nonequilibrium Molecular Dynamics Evaluation of Thermal Conductivity in 2D Yukawa Liquids

Kinetic Studies of Exchange-Correlation Effect on the Collective Excitations of Warm Dense Plasmas.

Time correlation functions and transport coefficients of two-dimensional Yukawa liquids

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

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

Spectral Decomposition of Thermal Conductivity: Comparing Velocity Decomposition Methods in Homogeneous Molecular Dynamics Simulations

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

Shear viscosity and shear thinning in two-dimensional Yukawa liquids

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

Heat Conduction In Homogeneous And Heterogeneous Billiard Systems

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

Vibrational model of heat transfer in strongly coupled Yukawa fluids (dusty plasma liquids)

Heat conduction in 2D strongly-coupled dusty plasmas

Thermal conduction in two-dimensional complex plasma layers

Isomorph invariance and thermodynamics of repulsive dense bi-Yukawa one-component plasmas

Structure of inhomogeneous attractive and repulsive hard-core yukawa fluid: grand canonical Monte Carlo simulation and density functional theory study.

Longitudinal viscosity of 2D Yukawa liquids

Thermal motions in complex liquids: the 2D Lennard-Jones liquid

Homogeneous nonequilibrium molecular dynamics method for heat transport and spectral decomposition with many-body potentials

Rapid heating and cooling in two-dimensional Yukawa systems