Lattice thermal transport from phonon spectra beyond perturbation theory

Abstract

We develop a molecular dynamics framework to compute the mode-resolved phonon spectral density from classical correlations of an annihilation-like phonon variable. For harmonic oscillators, classical molecular dynamics exactly reproduces the corresponding quantum Kubo-transformed correlator, providing the basis for extension to anharmonic systems. Using PbTe as a benchmark and Cs${}_3$Bi${}_2$I${}_6$Cl${}_3$ as a strongly anharmonic test case, we show that the method captures both quasiparticle and non-Lorentzian spectra beyond perturbative quasiparticle theory, while yielding thermal conductivity in good agreement with experiment. This framework provides a classical molecular dynamics route to mode-resolved phonon spectral densities for spectral Wigner heat transport in strongly anharmonic solids above their Debye temperatures.