Spectral decomposition of thermal conductivity: Comparing velocity decomposition methods in homogeneous molecular dynamics simulations
Loading...
Access rights
openAccess
publishedVersion
URL
Journal Title
Journal ISSN
Volume Title
A1 Alkuperäisartikkeli tieteellisessä aikakauslehdessä
This publication is imported from Aalto University research portal.
View publication in the Research portal (opens in new window)
View/Open full text file from the Research portal (opens in new window)
Other link related to publication (opens in new window)
View publication in the Research portal (opens in new window)
View/Open full text file from the Research portal (opens in new window)
Other link related to publication (opens in new window)
Date
2021-05-17
Department
Major/Subject
Mcode
Degree programme
Language
en
Pages
11
Series
Physical Review B, Volume 103, issue 20
Abstract
The design of 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 (MA) 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 MA implementation in the Graphics Processing Units Molecular Dynamics code on a single graphics processing unit is over 1000 times faster than the existing implementation in the Large-scale Atomic/Molecular Massively Parallel Simulator code on one central processing unit.Description
Funding Information: Some of the computing for this project was performed on the Sherlock cluster at Stanford University. We would like to thank Stanford University and the Stanford Research Computing Center (SRCC) for providing computational resources and support that contributed to these results. This paper was also partially supported by ASCENT, one of the six centers in JUMP, a Semiconductor Research Corporation (SRC) program sponsored by DARPA. A.J.G. also acknowledges support from the Achievement Rewards for College Scientists (ARCS) Northern California Chapter. T.A-N. has been supported in part by the Academy of Finland through its QTF Center of Excellence (Project No. 312298). Z.F. has been supported by the National Natural Science Foundation of China (NSFC; No. 11974059). Publisher Copyright: © 2021 American Physical Society.
Keywords
Other note
Citation
Gabourie, A J, Fan, Z, Ala-Nissila, T & Pop, E 2021, ' Spectral decomposition of thermal conductivity : Comparing velocity decomposition methods in homogeneous molecular dynamics simulations ', Physical Review B, vol. 103, no. 20, 205421 . https://doi.org/10.1103/PhysRevB.103.205421