English

Computation of thermal conductivity based on Path Integral Monte Carlo methods

Statistical Mechanics 2026-02-19 v1

Abstract

The calculation of thermal conductivity in insulating solids at temperatures below the Debye temperature is problematic, due to the breakdown of classical and semi-classical approaches. In this work, we present a fully quantum methodology to compute thermal conductivity based on Path Integral Monte Carlo (PIMC) simulations combined with the Green-Kubo linear response theory. The method is applied to crystalline argon modeled by a Lennard-Jones potential, a paradigmatic system where quantum effects strongly affect both thermodynamic and transport properties. From PIMC simulations, we obtain the temperature-dependent phonon frequencies, lifetimes, and specific heat. From the imaginary time correlations of the energy current, we extract the thermal transport coefficients based on a physically motivated prior. We show that the experimentally observed increase of the thermal conductivity at low temperatures cannot be explained within a standard Peierls-Boltzmann framework or quasi-harmonic approximation using phonon lifetimes alone. Instead, a distinct transport lifetime emerges from the analysis of heat-current correlations. Our results demonstrate that quantum Monte Carlo methods provide a robust, non-perturbative framework to investigate heat transport in insulating solids, beyond the limits of classical molecular dynamics and quasi-harmonic approximations.

Keywords

Cite

@article{arxiv.2602.16405,
  title  = {Computation of thermal conductivity based on Path Integral Monte Carlo methods},
  author = {Vladislav Efremkin and Stefano Mossa and Jean-Louis Barrat and Markus Holzmann},
  journal= {arXiv preprint arXiv:2602.16405},
  year   = {2026}
}

Comments

6 pages, 5 figures

R2 v1 2026-07-01T10:41:14.089Z