We propose a novel technique of molecular dynamics simulation to evaluate the relaxation time of phonons in solids for investigation of solid heat conductivity. The basic idea is to observe relaxation behavior of the power spectrum of atomic velocities after energetically stimulating modes in a specific frequency region. The transient entropy S(t) is defined with the power spectrum based on non-equilibrium statistical mechanics to quantitatively evaluate the relaxation speed. In this paper, two example systems are shown: a Lennard-Jones model crystal and a silicon crystal. For both systems, we found that the observed S(t) is well fitted to a single exponential function, from which we can obtain a frequency-dependent relaxation time.
- Heat Transfer Division
Relaxation of Phonons in Classical MD Simulation
- Views Icon Views
- Share Icon Share
- Search Site
Matsumoto, M, Kunisawa, T, & Xiao, P. "Relaxation of Phonons in Classical MD Simulation." Proceedings of the ASME/JSME 2007 Thermal Engineering Heat Transfer Summer Conference collocated with the ASME 2007 InterPACK Conference. ASME/JSME 2007 Thermal Engineering Heat Transfer Summer Conference, Volume 1. Vancouver, British Columbia, Canada. July 8–12, 2007. pp. 77-80. ASME. https://doi.org/10.1115/HT2007-32659
Download citation file: