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

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.

arXiv comment19 Pages, 8 Figures

References in corpus (12)

• Optimized Tersoff and Brenner empirical potential parameters for lattice dynamics and phonon thermal transport in carbon nanotubes and graphene
• Force and heat current formulas for many-body potentials in molecular dynamics simulation with applications to thermal conductivity calculations
• Efficient molecular dynamics simulations with many-body potentials on graphics processing units
• Homogeneous nonequilibrium molecular dynamics method for heat transport and spectral decomposition with many-body potentials
• Frequency-dependent phonon mean free path in carbon nanotubes from non-equilibrium molecular dynamics
• Thermal conductivity decomposition in two-dimensional materials: Application to graphene
• Divergence of the Thermal Conductivity in Uniaxially Strained Graphene
• Quantitatively Analyzing Phonon Spectral Contribution of Thermal Conductivity Based on Non-Equilibrium Molecular Dynamics Simulation I: From Space Fourier Transform
• Quantitatively Analyzing Phonon Spectral Contribution of Thermal Conductivity Based on Non-Equilibrium Molecular Dynamics Simulation II: From Time Fourier Transform
• Thermal Transport in MoS${}_2$ from Molecular Dynamics using Different Empirical Potentials
• Reduced Thermal Conductivity of Supported and Encased Monolayer and Bilayer MoS${}_2$
• A minimal Tersoff potential for diamond silicon with improved descriptions of elastic and phonon transport properties

Cited by in corpus (19)

• Neuroevolution machine learning potentials: Combining high accuracy and low cost in atomistic simulations and application to heat transport
• GPUMD: A package for constructing accurate machine-learned potentials and performing highly efficient atomistic simulations
• Molecular dynamics simulations of heat transport using machine-learned potentials: A mini review and tutorial on GPUMD with neuroevolution potentials
• Quantum-corrected thickness-dependent thermal conductivity in amorphous silicon predicted by machine-learning molecular dynamics simulations
• Spin-Dependent Graph Neural Network Potential for Magnetic Materials
• Variable thermal transport in black, blue, and violet phosphorene from extensive atomistic simulations with a neuroevolution potential
• Combining the D3 dispersion correction with the neuroevolution machine-learned potential
• Dissimilar thermal transport properties in $\kappa$-Ga${}_2$O${}_3$ and $\beta$-Ga${}_2$O${}_3$ revealed by machine-learning homogeneous nonequilibrium molecular dynamics simulations
• Heat transport across graphene/hexagonal-BN tilted grain boundaries from phase-field crystal model and molecular dynamics simulations
• Phonon coherence and minimum thermal conductivity in disordered superlattice
• Nontrivial impact of interlayer coupling on thermal conductivity: opposing trends in in-plane and out-of-plane phonons
• Copper delocalization leads to ultralow thermal conductivity in chalcohalide CuBiSeCl2
• Ultrahigh interfacial thermal conductance for cooling gallium oxide electronics using cubic boron arsenide
• Atomic-Scale Probing of Heterointerface Phonon Bridges in Nitride Semiconductor
• Role of octahedral tilting induced acoustic softening on limiting thermal transport in SrSnO3
• Uniaxial stress tuning of interfacial thermal conductance in cubic BAs/4H-SiC heterostructures
• Substrate-Dependence of Monolayer MoS${}_2$ Thermal Conductivity and Thermal Boundary Conductance
• LLM-driven discovery for carbon allotropes with bond-network entropy
• Interface phonon modes governing the ideal limit of thermal transport across diamond/cubic boron nitride interfaces