Predicting Phonon Properties From Molecular Dynamics Simulations Using the Spectral Energy Density

Author(s):  
Joseph E. Turney ◽  
John A. Thomas ◽  
Alan J. H. McGaughey ◽  
Cristina H. Amon

Using lattice dynamics theory, we derive the spectral energy density and the relation between the spectral energy density and the phonon frequencies and relaxation times. We then calculate the spectral energy density and phonon frequencies and relaxation times for a test system of Lennard-Jones argon using velocities obtained from molecular dynamics simulations. The phonon properties, which can be used to calculate thermal conductivity, are compared to predictions made using (i) anharmonic lattice dynamics calculations and (ii) a technique that performs normal mode analysis on the positions and velocities obtained from molecular dynamics simulations.

Author(s):  
Jason M. Larkin ◽  
Alexandre D. Massicotte ◽  
Joseph E. Turney ◽  
Alan J. H. McGaughey ◽  
Cristina H. Amon

To predict the thermal conductivity of a dielectric or insulating material requires the phonon frequencies and lifetimes. Techniques for predicting these quantities have been proposed based in molecular dynamics simulation and lattice dynamics calculations. Here, two expressions for the phonon spectral energy density are described and applied to two test systems: Lennard-Jones argon and Stillinger-Weber silicon. One spectral energy density expression is derived from lattice dynamics theory, while the other uses only the atomic velocities from molecular dynamics simulation. We find that while the spectral energy density using only atomic velocities can predict the phonon frequencies, it is not generally able to predict the lifetimes due to terms omitted in the derivation.


Author(s):  
Asegun S. Henry ◽  
Gang Chen

Silicon's material properties, have been studied extensively because of its technological significance in a variety of industries, including microelectronics. Yet, questions surrounding the phonon relaxation times in silicon continue to linger.1,2 Previous theoretical works3-5 have generated qualitative expressions for phonon relaxation times, however these approaches require fitting parameters that cannot be determined reliably. This paper first discusses implementation issues associated with using the Green-Kubo method in molecular dynamics (MD) simulations. We compare various techniques used in similar works and discusses several implementation issues that have arisen in the literature. We then describe an alternative procedure for analyzing the normal modes of a crystal to extract phonon relaxation times. As an example material we study bulk crystalline silicon using equilibrium MD simulations and lattice dynamics. The environment dependent interatomic potential6 is used to model the interactions and frequency dependent phonon properties are extracted from the MD simulations.


Author(s):  
J. A. Thomas ◽  
J. E. Turney ◽  
R. M. Iutzi ◽  
A. J. H. McGaughey ◽  
C. H. Amon

Thermal transport by phonons in water/carbon nanotube (CNT) composite systems is investigated using molecular dynamics (MD) simulation. We calculate the spectral energy density of empty and water-filled CNTs and use it to extract the mode-specific phonon group velocities, relaxation times, and thermal conductivities. The total thermal conductivity predicted from the spectral energy density is consistent with what we predict using a direct application of the Fourier law in an MD simulation. The number of atoms and simulation runtime required to predict the spectral energy density, however, are both at least one order-of-magnitude smaller.


2015 ◽  
Vol 7 (17) ◽  
pp. 2317-2331 ◽  
Author(s):  
Gautier Moroy ◽  
Olivier Sperandio ◽  
Shakti Rielland ◽  
Saurabh Khemka ◽  
Karen Druart ◽  
...  

ChemMedChem ◽  
2010 ◽  
Vol 5 (3) ◽  
pp. 443-454 ◽  
Author(s):  
Torsten Luksch ◽  
Andreas Blum ◽  
Nina Klee ◽  
Wibke E. Diederich ◽  
Christoph A. Sotriffer ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document