Accurately predicting the thermal conductivity of quantum solids
Researchers have developed a numerical method that enables the thermal behaviour of solids at very low temperatures to be calculated without approximations, when quantum effects influence heat diffusion
References:
Vladislav Efremkin, Stefano Mossa, Jean-Louis Barrat, and Markus Holzmann, Non-perturbative computation of thermal conductivity based on path integral Monte Carlo methods, Physical Review Letters - Published 25 August, 2026.
DOI: https://doi.org/10.1103/2f5l-scbv
Open access: arXiv
Heat transfer is a fundamental property of materials that occurs whenever there is a temperature difference, and is described by Fourier’s well-known law, which involves thermal conductivity – a coefficient that reflects a material’s greater or lesser ability to diffuse energy under the influence of a given temperature gradient. Understanding the microscopic factors that contribute to higher or lower thermal conductivity is therefore crucial for applications ranging from energy technologies to electronics and the design of high-performance insulation materials.
This research was carried out in the following CNRS laboratories:
In a recent study, researchers have proposed a new computational method for studying thermal conduction in insulating or semiconducting solid systems, based on path integral Monte Carlo (PIMC) simulations. This complex method calculates the dynamics of the nuclei whilst taking precise account of all quantum effects, which are present in all solids but are particularly pronounced at low temperatures and in crystals containing light atoms. It utilises the formal relationship which, in quantum physics, enables the calculation—via Monte Carlo simulations at equilibrium—of the time-dependent correlation functions, which are obtained in the imaginary time domain and must then be analytically extended into the real time domain to obtain the transport coefficients. Unlike many recently developed methods, this approach requires neither simplifying assumptions nor perturbative approximations, and it takes into account all the degrees of freedom involved in heat transport. The researchers validate the method by studying crystals of noble gases, argon and neon, for which they obtain results consistent with experimental measurements. The results are obtained by constraining the spectral function associated with the heat flux to take the form of a physically relevant model, constructed from the effective phonon frequencies. They show that the qualitative change in behaviour at low temperatures is associated with a transport characteristic time that diverges as T→ 0, differing from phonons, whose lifetime remains finite in this limit.
This new method, which provides a robust framework for studying heat transport beyond the capabilities of classical molecular dynamics, without resorting to perturbative or semi-classical approximations, should lead to a better understanding of the microscopic characteristics of matter that determine heat transport properties, particularly in low-temperature regimes where quantum effects become significant. These results have been published in the Physical Review Letters.