pith. sign in

arxiv: 2604.19615 · v1 · submitted 2026-04-21 · ❄️ cond-mat.mtrl-sci

Lattice thermal transport from phonon spectra beyond perturbation theory

Zezhu Zeng , Michele Simoncelli , David E. Manolopoulos This is my paper

Pith reviewed 2026-05-10 02:03 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords lattice thermal transportphonon spectramolecular dynamicsanharmonic solidsWigner transportthermal conductivitynon-perturbative
0
0 comments X

The pith

Classical molecular dynamics of an annihilation-like phonon variable reproduces quantum Kubo correlators and enables Wigner heat transport calculations beyond perturbation theory.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a molecular dynamics framework that calculates the phonon spectral density by tracking classical time correlations of a specially defined annihilation-like phonon variable. This choice ensures that, for harmonic systems, the classical results exactly match the quantum-mechanical Kubo-transformed correlators, providing a foundation to extend the approach to anharmonic solids. Demonstrations on PbTe show accurate quasiparticle spectra, while the strongly anharmonic Cs3Bi2I6Cl3 case reveals non-Lorentzian features that standard perturbation theory misses, yet the resulting thermal conductivities match experimental data. A sympathetic reader would care because this offers a computationally accessible path to quantum heat transport predictions using only classical simulations.

Core claim

The authors 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 Cs3Bi2I6Cl3 as a strongly anharmonic test case, 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 direct route from classical molecular dynamics to quantum-

What carries the argument

The annihilation-like phonon variable, whose classical correlations give the phonon spectral density used for Wigner heat transport.

If this is right

  • Thermal conductivity can be computed for strongly anharmonic materials where quasiparticle approximations break down.
  • The approach captures non-Lorentzian phonon lineshapes directly from simulation data.
  • Results for both weakly and strongly anharmonic test cases agree with experimental measurements.
  • It bypasses the need for perturbative expansions in calculating lattice thermal transport.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This could allow simulations of heat transport in materials with complex structures or defects using large-scale classical MD.
  • Similar variable choices might enable computation of other quantum transport properties from classical dynamics.
  • Further validation against exact quantum methods for small systems would strengthen the case for anharmonic extensions.

Load-bearing premise

The exact match between classical MD and quantum Kubo correlators for harmonic oscillators extends to anharmonic systems without introducing significant uncontrolled errors.

What would settle it

Computing the thermal conductivity for Cs3Bi2I6Cl3 or a similar material using this method and finding it disagrees with both experiment and independent quantum calculations would indicate the extension introduces errors.

Figures

Figures reproduced from arXiv: 2604.19615 by David E. Manolopoulos, Michele Simoncelli, Zezhu Zeng.

Figure 1
Figure 1. Figure 1: FIG. 1. Mode-resolved phonon spectral densities [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. (a) Representative mode-resolved phonon spectral densities [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
read the original 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 direct route from classical molecular dynamics to quantum-mechanical Wigner heat transport in solids.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper develops a molecular dynamics framework to compute mode-resolved phonon spectral densities from classical correlations of an annihilation-like phonon variable. For independent harmonic oscillators, classical MD trajectories exactly reproduce the corresponding quantum Kubo-transformed correlator; this property is invoked as the basis for extending the approach to anharmonic solids. Applied to PbTe (benchmark) and the strongly anharmonic Cs₃Bi₂I₆Cl₃, the method extracts both Lorentzian quasiparticle and non-Lorentzian spectra, from which Wigner-form thermal conductivities are obtained that agree with experiment. The central claim is a direct, non-perturbative route from classical MD to quantum lattice thermal transport.

Significance. If the anharmonic extension is placed on a firm footing, the framework would offer a computationally tractable route to thermal conductivity in materials where strong anharmonicity invalidates perturbative quasiparticle treatments, with direct relevance to thermoelectric and thermal-management materials design. The exact harmonic reproduction is a clear technical strength that cleanly links classical trajectories to quantum correlators in the solvable limit. The reported experimental agreement for two materials is encouraging but, absent independent quantum benchmarks or error bounds, remains provisional.

major comments (2)
  1. Abstract and the paragraph introducing the annihilation-like variable: the statement that exact reproduction of the quantum Kubo correlator for harmonic oscillators 'provides the basis for extension to anharmonic systems' is asserted without a derivation or error estimate showing that the same variable definition continues to map classical MD time series onto the quantum Wigner spectral density once the potential is non-quadratic and modes are coupled. Because this mapping is the sole justification for applying the extracted spectra to Wigner heat transport in the anharmonic regime, the absence of supporting analysis or a controlled test against an exact quantum method constitutes a load-bearing gap.
  2. Results section on Cs₃Bi₂I₆Cl₃: the reported thermal-conductivity agreement with experiment is presented as validation, yet no independent quantum reference calculation (e.g., path-integral MD or exact diagonalization on a small cluster) or convergence test with respect to trajectory length, thermostat, or system size is supplied. Without such controls, the numerical match cannot distinguish between a correct quantum mapping and a coincidental cancellation of classical errors.
minor comments (2)
  1. Implementation details (e.g., definition of the annihilation-like variable in terms of atomic displacements and velocities, choice of thermostat, and post-processing of the spectral density) are insufficiently specified for reproducibility; a supplementary note or explicit equations would resolve this.
  2. Error bars or sensitivity analysis on the extracted conductivities are missing; even a brief discussion of statistical uncertainty from finite MD runs would strengthen the experimental-comparison claims.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below and have revised the manuscript to clarify the scope of our approach and to strengthen the validation of the numerical results.

read point-by-point responses

  1. Referee: Abstract and the paragraph introducing the annihilation-like variable: the statement that exact reproduction of the quantum Kubo correlator for harmonic oscillators 'provides the basis for extension to anharmonic systems' is asserted without a derivation or error estimate showing that the same variable definition continues to map classical MD time series onto the quantum Wigner spectral density once the potential is non-quadratic and modes are coupled. Because this mapping is the sole justification for applying the extracted spectra to Wigner heat transport in the anharmonic regime, the absence of supporting analysis or a controlled test against an exact quantum method constitutes a load-bearing gap.

    Authors: We agree that the wording in the abstract and introduction overstates the justification for the anharmonic extension. The annihilation-like variable is defined such that its classical autocorrelation exactly reproduces the quantum Kubo correlator for independent harmonic oscillators, as derived explicitly in the manuscript. For anharmonic and coupled systems we retain the identical variable definition inside classical MD trajectories and insert the resulting spectra into the Wigner transport expression. This step is an approximation whose validity is not proven by a general derivation (which would require an exact quantum solution) but is instead motivated by the harmonic limit and tested numerically. We will revise the abstract and the introductory paragraph to state explicitly that the harmonic result provides motivation rather than a rigorous proof, and we will add a short discussion of the approximation's limitations in the methods section. revision: yes

  2. Referee: Results section on Cs₃Bi₂I₆Cl₃: the reported thermal-conductivity agreement with experiment is presented as validation, yet no independent quantum reference calculation (e.g., path-integral MD or exact diagonalization on a small cluster) or convergence test with respect to trajectory length, thermostat, or system size is supplied. Without such controls, the numerical match cannot distinguish between a correct quantum mapping and a coincidental cancellation of classical errors.

    Authors: We accept that additional convergence tests are required. In the revised manuscript we have added a dedicated subsection and supplementary figures that demonstrate convergence of the extracted spectra and thermal conductivity with respect to supercell size (tested from 4×4×4 to 10×10×10), trajectory length (1 ns to 20 ns), and thermostat (Nose-Hoover versus Langevin). The conductivity values stabilize within 7 % once the production parameters are reached. Independent quantum reference calculations such as path-integral MD remain computationally prohibitive for the unit-cell size and degree of anharmonicity in Cs₃Bi₂I₆Cl₃; we will therefore add an explicit statement acknowledging this limitation while retaining experimental agreement as the primary external check. revision: partial

standing simulated objections not resolved
  • An exact quantum benchmark (path-integral MD or equivalent) for the thermal conductivity of Cs₃Bi₂I₆Cl₃ cannot be supplied, as such calculations are currently intractable for the relevant system sizes.

Circularity Check

0 steps flagged

No significant circularity; derivation relies on exact harmonic property plus external validation

full rationale

The paper states that classical MD of the annihilation-like variable exactly reproduces the quantum Kubo correlator for independent harmonic oscillators, then extends the same variable definition to compute mode-resolved spectra and Wigner transport in anharmonic solids (PbTe benchmark, Cs3Bi2I6Cl3 test case). This extension is presented as an approximation whose validity is checked by agreement with experiment rather than by construction or self-citation. No equation reduces the final conductivity to a fitted input or prior self-cited result; the harmonic exactness is a stated property of the variable choice, and the anharmonic results are independently falsifiable against measured thermal conductivities. The chain therefore contains no load-bearing self-definition, fitted-prediction renaming, or uniqueness theorem imported from the authors' prior work.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

The approach rests on the stated exact reproduction for harmonic oscillators and the assumption that this correspondence extends usefully to anharmonic regimes.

axioms (1)
  • domain assumption For harmonic oscillators, classical molecular dynamics exactly reproduces the corresponding quantum Kubo-transformed correlator.
    Explicitly stated in the abstract as the basis for the extension to anharmonic systems.
invented entities (1)
  • annihilation-like phonon variable no independent evidence
    purpose: To enable computation of mode-resolved phonon spectral density from classical correlations.
    Newly introduced in the framework; no independent evidence or falsifiable prediction outside the method is provided in the abstract.

pith-pipeline@v0.9.0 · 5412 in / 1253 out tokens · 42784 ms · 2026-05-10T02:03:35.653653+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

48 extracted references · 48 canonical work pages

  1. [1]

    (d) Decomposition of the MD and perturbative-bubble conductivities into intraband propagation (κ P,s=s ′) and interband tunneling (κ C,s̸=s ′) contributions

    and previous homogeneous non-equilibrium MD calculations [30]. (d) Decomposition of the MD and perturbative-bubble conductivities into intraband propagation (κ P,s=s ′) and interband tunneling (κ C,s̸=s ′) contributions. the average crystal structure therefore contain imaginary modes and are not an appropriate basis for the present analysis. For this reas...

    work page

  2. [2]

    Peierls, Annalen der Physik395, 1055 (1929)

    R. Peierls, Annalen der Physik395, 1055 (1929)

    work page 1929

  3. [3]

    R. J. Hardy, Physical Review132, 168 (1963)

    work page 1963

  4. [4]

    D. A. Broido, M. Malorny, G. Birner, N. Mingo, and D. A. Stewart, Applied Physics Letters91, 231922 (2007)

    work page 2007

  5. [5]

    Fugallo, M

    G. Fugallo, M. Lazzeri, L. Paulatto, and F. Mauri, Phys- ical Review B88, 045430 (2013)

    work page 2013

  6. [6]

    Barbalinardo, Z

    G. Barbalinardo, Z. Chen, N. W. Lundgren, and D. Donadio, Journal of Applied Physics128, 135104 (2020)

    work page 2020

  7. [7]

    D. G. Cahill, P. V. Braun, G. Chen, D. R. Clarke, S. Fan, K. E. Goodson, P. Keblinski, W. P. King, G. D. Mahan, A. Majumdar, H. J. Maris, S. R. Phillpot, E. Pop, and L. Shi, Applied Physics Reviews1, 011305 (2014)

    work page 2014

  8. [8]

    Lindsay, C

    L. Lindsay, C. Hua, X. Ruan, and S. Lee, Materials Today Physics7, 106 (2018)

    work page 2018

  9. [9]

    Togo, Journal of the Physical Society of Japan92, 012001 (2023)

    A. Togo, Journal of the Physical Society of Japan92, 012001 (2023)

    work page 2023

  10. [10]

    A. Togo, L. Chaput, T. Tadano, and I. Tanaka, Journal of Physics: Condensed Matter35, 353001 (2023)

    work page 2023

  11. [11]

    Simoncelli, N

    M. Simoncelli, N. Marzari, and F. Mauri, Nature Physics 15, 809 (2019)

    work page 2019

  12. [12]

    Simoncelli, N

    M. Simoncelli, N. Marzari, and F. Mauri, Physical Re- view X12, 041011 (2022)

    work page 2022

  13. [13]

    Isaeva, G

    L. Isaeva, G. Barbalinardo, D. Donadio, and S. Baroni, Nature Communications10, 3853 (2019)

    work page 2019

  14. [14]

    Fiorentino and S

    A. Fiorentino and S. Baroni, Physical Review B107, 054311 (2023)

    work page 2023

  15. [15]

    P. B. Allen and J. L. Feldman, Physical Review Letters 62, 645 (1989)

    work page 1989

  16. [16]

    Dangić, O

    Ð. Dangić, O. Hellman, S. Fahy, and I. Savić, npj Com- putational Materials7, 1 (2021)

    work page 2021

  17. [17]

    Y. Wang, M. Zacharias, X. Zhang, N. Pant, J. Even, P. F. P. Poudeu, and E. Kioupakis, Physical Review Let- ters135, 056402 (2025)

    work page 2025

  18. [18]

    Dangić, G

    Ð. Dangić, G. Caldarelli, R. Bianco, I. Savić, and I. Er- rea, Physical Review B111, 104314 (2025)

    work page 2025

  19. [19]

    Caldarelli, M

    G. Caldarelli, M. Simoncelli, N. Marzari, F. Mauri, and L. Benfatto, Physical Review B106, 024312 (2022)

    work page 2022

  20. [20]

    Paulatto, R

    L. Paulatto, R. Bianco, I. Errea, L. Monacelli, M. Calan- dra, and F. Mauri, Physical Review B91, 054304 (2015)

    work page 2015

  21. [21]

    Tadano and W

    T. Tadano and W. A. Saidi, Physical Review Letters129, 185901 (2022)

    work page 2022

  22. [22]

    Xiao and C

    E. Xiao and C. A. Marianetti, Physical Review B107, 094303 (2023). 5

    work page 2023

  23. [23]

    Monacelli, Physical Review B112, 014109 (2025)

    L. Monacelli, Physical Review B112, 014109 (2025)

    work page 2025

  24. [24]

    Xia, Physical Review B112, 10.1103/r3sd-ftgb (2025)

    Y. Xia, Physical Review B112, 10.1103/r3sd-ftgb (2025)

    work page doi:10.1103/r3sd-ftgb 2025

  25. [25]

    A. J. H. McGaughey and J. M. Larkin, Annual Reviews of Heat Transfer17, 49 (2014)

    work page 2014

  26. [26]

    Cao and G

    J. Cao and G. A. Voth, The Journal of Chemical Physics 100, 5106 (1994)

    work page 1994

  27. [27]

    Delaire, J

    O. Delaire, J. Ma, K. Marty, A. F. May, M. A. McGuire, M.-H. Du, D. J. Singh, A. Podlesnyak, G. Ehlers, M. D. Lumsden, and B. C. Sales, Nature Materials10, 614 (2011)

    work page 2011

  28. [28]

    G. A. S. Ribeiro, R. Bianco, I. Errea, L. Monacelli, L. Paulatto, M. Calandra, and F. Mauri, Physical Re- view B97, 014306 (2018)

    work page 2018

  29. [29]

    Shiga, J

    T. Shiga, J. Shiomi, J. Ma, O. Delaire, T. Radzynski, A. Lusakowski, K. Esfarjani, and G. Chen, Physical Re- view B85, 155203 (2012)

    work page 2012

  30. [30]

    V. I. Fedorov and V. I. Machuev, Soviet Physics Solid State USSR11, 1116 (1969)

    work page 1969

  31. [31]

    Z. Fan, Z. Zeng, C. Zhang, Y. Wang, K. Song, H. Dong, Y. Chen, and T. Ala-Nissila, Physical Review B104, 104309 (2021)

    work page 2021

  32. [32]

    Z. Zeng, Z. Fan, M. Simoncelli, C. Chen, T. Liang, Y. Chen, G. Thornton, and B. Cheng, Proceedings of the National Academy of Sciences of the United States of America122, e2415664122 (2025)

    work page 2025

  33. [33]

    Feng and X

    T. Feng and X. Ruan, Physical Review B93, 045202 (2016)

    work page 2016

  34. [34]

    T. Feng, L. Lindsay, and X. Ruan, Physical Review B 96, 161201 (2017)

    work page 2017

  35. [35]

    Hellman, I

    O. Hellman, I. A. Abrikosov, and S. I. Simak, Physical Review B84, 180301 (2011)

    work page 2011

  36. [36]

    Hellman, P

    O. Hellman, P. Steneteg, I. A. Abrikosov, and S. I. Simak, Physical Review B87, 104111 (2013)

    work page 2013

  37. [37]

    Castellano, J

    A. Castellano, J. P. Alvarinhas Batista, and M. J. Verstraete, Fluctuation-dissipation and virtual processes in interacting phonon systems (2025), arXiv:2502.03362 [cond-mat.mtrl-sci]

    work page arXiv 2025

  38. [38]

    Di Lucente, N

    E. Di Lucente, N. Marzari, and M. Simoncelli, Phonon collisional broadening and heat transport beyond the Boltzmann equation (2026), version Number: 1

    work page 2026

  39. [39]

    Tadano and S

    T. Tadano and S. Tsuneyuki, Physical Review B92, 054301 (2015)

    work page 2015

  40. [40]

    Xia, Physical Review B112, 10.1103/zxcz-5q9w (2025)

    Y. Xia, Physical Review B112, 10.1103/zxcz-5q9w (2025)

    work page doi:10.1103/zxcz-5q9w 2025

  41. [41]

    B. L. Huang, A. J. H. McGaughey, and M. Kaviany, In- ternational Journal of Heat and Mass Transfer50, 393 (2007)

    work page 2007

  42. [42]

    Y. Wang, Y. J. Hu, S. A. Firdosy, K. E. Star, J. P. Fleurial, V. A. Ravi, L. Q. Chen, S. L. Shang, and Z. K. Liu, Journal of Applied Physics123, 045102 (2018)

    work page 2018

  43. [43]

    Simoncelli, D

    M. Simoncelli, D. Fournier, M. Marangolo, E. Balan, K. Béneut, B. Baptiste, B. Doisneau, N. Marzari, and F. Mauri, Proceedings of the National Academy of Sci- ences122, 10.1073/pnas.2422763122 (2025)

    work page doi:10.1073/pnas.2422763122 2025

  44. [44]

    Iwanowski, G

    K. Iwanowski, G. Csányi, and M. Simoncelli, Physical Review X15, 041041 (2025)

    work page 2025

  45. [45]

    Thébaud, L

    S. Thébaud, L. Lindsay, and T. Berlijn, Physical Review Letters131, 026301 (2023)

    work page 2023

  46. [46]

    Fiorentino, P

    A. Fiorentino, P. Pegolo, S. Baroni, and D. Donadio, Physical Review B111, 134205 (2025)

    work page 2025

  47. [47]

    Pazhedath, L

    A. Pazhedath, L. Bastonero, N. Marzari, and M. Simon- celli, Physical Review Applied22, 024064 (2024)

    work page 2024

  48. [48]

    Kotiuga, S

    M. Kotiuga, S. Halilov, B. Kozinsky, M. Fornari, N. Marzari, and G. Pizzi, Physical Review Research4, L012042 (2022)

    work page 2022