pith. sign in

arxiv: 2605.20096 · v1 · pith:VZZP7M6Wnew · submitted 2026-05-19 · ❄️ cond-mat.mtrl-sci

Lattice thermal conductivity decomposition: Peierls vs. non-Peierls contributions

Andrey Pereverzev This is my paper

Pith reviewed 2026-05-20 03:30 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords lattice thermal conductivityGreen-Kubo methodPeierls heat currentphonon transportalpha-quartzsolid argonrelaxation time approximationnon-Peierls contributions
0
0 comments X

The pith

Quadratic and Peierls heat currents produce nearly identical lattice thermal conductivities in crystals.

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

This paper compares the Green-Kubo lattice thermal conductivity from the full classical heat current against its quadratic component and the standard Peierls heat current. Calculations on solid argon, a mass-alternating argon model, and alpha-quartz show only slight differences between the quadratic and Peierls results. In alpha-quartz the optical phonon contribution exceeds the acoustic one. The relaxation time approximation underestimates conductivity across all three systems. The work therefore tests whether non-Peierls terms matter for heat transport in these solids.

Core claim

The Green-Kubo lattice thermal conductivity computed using the full classical heat current of a crystalline solid is compared with results obtained from the quadratic component of the heat current and from the commonly used Peierls heat current. For all materials considered, the thermal conductivities calculated using the quadratic and Peierls heat currents differ only slightly. In the case of α-quartz, the optical phonon contribution to the thermal conductivity is found to exceed that of the acoustic modes. The relaxation time approximation systematically underestimates the thermal conductivity in all three systems.

What carries the argument

Decomposition of the classical heat current into quadratic and Peierls components inside the Green-Kubo formula.

If this is right

  • The Peierls heat current remains a reliable approximation for thermal conductivity in the studied crystals.
  • Optical phonon modes can carry more heat than acoustic modes in materials such as alpha-quartz.
  • The relaxation time approximation yields systematically lower values than the full Green-Kubo approach.

Where Pith is reading between the lines

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

  • The close agreement suggests that higher-order heat-current terms can often be neglected in similar insulators.
  • The optical-mode dominance in quartz points to a need to recheck mode-resolved contributions in other non-cubic crystals.
  • If the pattern holds, simpler Peierls-based models could be used for rapid screening of thermal transport in related materials.

Load-bearing premise

The decomposition of the classical heat current into quadratic and Peierls components is assumed to be physically meaningful and sufficient to isolate non-Peierls contributions without significant cross terms or higher-order effects.

What would settle it

A calculation or measurement that finds a large numerical difference between the quadratic and Peierls thermal conductivities in solid argon or alpha-quartz.

Figures

Figures reproduced from arXiv: 2605.20096 by Andrey Pereverzev.

Figure 1
Figure 1. Figure 1: FIG. 1. (a) Solid argon treated as a tetragonal lattice. The tetrago [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Phonon branches of the SA (red) and SAAM (blue) models [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
read the original abstract

The Green-Kubo lattice thermal conductivity computed using the full classical heat current of a crystalline solid is compared with results obtained from the quadratic component of the heat current and from the commonly used Peierls heat current. In addition, thermal conductivity within the relaxation time approximation is evaluated. Three crystalline systems are investigated: solid argon, a model of solid argon with alternating masses, and $\alpha$-quartz. For all materials considered, the thermal conductivities calculated using the quadratic and Peierls heat currents differ only slightly. In the case of $\alpha$-quartz, the optical phonon contribution to the thermal conductivity is found to exceed that of the acoustic modes. The relaxation time approximation systematically underestimates the thermal conductivity in all three systems.

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 manuscript compares Green-Kubo lattice thermal conductivities computed from the full classical heat current against those obtained from its quadratic component and from the standard Peierls heat current. Three systems are examined: solid argon, an alternating-mass model of solid argon, and α-quartz. The quadratic and Peierls conductivities are reported to differ only slightly in all cases. In α-quartz the optical-phonon contribution is found to exceed the acoustic contribution, while the relaxation-time approximation systematically underestimates the conductivity in every system examined.

Significance. If the decomposition is shown to be free of large cross terms, the work supplies concrete numerical evidence that the Peierls form remains a reasonable proxy for the full classical heat current even in systems where anharmonic or optical-mode effects are non-negligible. The optical-versus-acoustic split in α-quartz supplies a falsifiable prediction that can be tested against other mode-resolved calculations or experiments.

major comments (2)
  1. [Results and discussion of heat-current decomposition] Green-Kubo decomposition (implied in the comparison of separate conductivities): the claim that the quadratic and Peierls conductivities “differ only slightly” does not automatically establish that the cross term 2∫⟨J_quad(t)·J_Peierls(0)⟩dt is negligible. Without an explicit plot or bound on this cross-correlation function, the physical separation into Peierls versus non-Peierls contributions remains unverified and is load-bearing for the central interpretation.
  2. [α-quartz results] α-quartz optical/acoustic split: the statement that optical modes contribute more than acoustic modes depends on a specific projection of the heat current onto phonon eigenvectors. The manuscript should state whether this projection assumes harmonic eigenvectors, includes anharmonic corrections, or employs any orthogonality approximation, because any such assumption directly affects the reported dominance.
minor comments (2)
  1. The reported conductivity values lack statistical uncertainties or error bars derived from the finite MD trajectories; these should be added so that the significance of the small observed differences can be assessed quantitatively.
  2. Notation for the quadratic heat-current component should be defined explicitly (e.g., which terms in the microscopic energy flux are retained) to avoid ambiguity with other common decompositions in the literature.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major point below and describe the revisions we will implement.

read point-by-point responses

  1. Referee: [Results and discussion of heat-current decomposition] Green-Kubo decomposition (implied in the comparison of separate conductivities): the claim that the quadratic and Peierls conductivities “differ only slightly” does not automatically establish that the cross term 2∫⟨J_quad(t)·J_Peierls(0)⟩dt is negligible. Without an explicit plot or bound on this cross-correlation function, the physical separation into Peierls versus non-Peierls contributions remains unverified and is load-bearing for the central interpretation.

    Authors: We agree that an explicit evaluation of the cross term provides stronger support for interpreting the results as a decomposition into largely independent Peierls and non-Peierls channels. In the revised manuscript we will add the integrated cross-correlation contribution (or a plot of the time-dependent cross-correlation function) for all three systems. Preliminary checks indicate this term is small relative to the diagonal contributions, consistent with the observed similarity between quadratic and Peierls conductivities, but we will present the data explicitly. revision: yes

  2. Referee: [α-quartz results] α-quartz optical/acoustic split: the statement that optical modes contribute more than acoustic modes depends on a specific projection of the heat current onto phonon eigenvectors. The manuscript should state whether this projection assumes harmonic eigenvectors, includes anharmonic corrections, or employs any orthogonality approximation, because any such assumption directly affects the reported dominance.

    Authors: The projection is performed using the eigenvectors obtained from the harmonic dynamical matrix evaluated at the relevant wave-vectors; no anharmonic corrections to the eigenvectors are included, and the standard orthogonality of the harmonic basis is assumed. We will add an explicit statement to this effect in the Methods section of the revised manuscript, together with a brief justification that this is the conventional approach for classical molecular-dynamics studies of mode-resolved transport. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected in derivation chain

full rationale

The paper computes lattice thermal conductivities for three crystalline systems via direct molecular-dynamics trajectories and the standard Green-Kubo integral applied to the full classical heat current, its quadratic component, and the Peierls heat current. These are independent numerical evaluations of correlation functions without parameter fitting to the target conductivities, without self-definitional equations that equate inputs to outputs by construction, and without load-bearing reliance on prior self-citations for uniqueness or ansatz choices. The relaxation-time-approximation comparison is likewise a separate standard evaluation. The central claims therefore remain self-contained against external simulation benchmarks rather than reducing to the paper's own inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The paper rests on standard classical statistical mechanics and the validity of the chosen interatomic potentials for argon and quartz; no new entities are postulated and no free parameters are fitted beyond those implicit in the potentials themselves.

axioms (1)
  • domain assumption Classical mechanics governs the atomic trajectories and heat current correlations in the Green-Kubo formula.
    Invoked throughout the Green-Kubo and relaxation-time calculations for all three systems.

pith-pipeline@v0.9.0 · 5650 in / 1337 out tokens · 42138 ms · 2026-05-20T03:30:53.767952+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    The third and fourth terms ... represent the quadratic part of the heat current, Jquad. Peierls heat current constitute a part of the quadratic current ... obtained by converting the stress-dependent part ... to normal mode coordinates and keeping only the terms that are diagonal in the phonon branch index

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

23 extracted references · 23 canonical work pages

  1. [1]

    with the same potential parameters. In order to achieve splitting of the phonon spectrum into acoustic and optical branches, two alternating atomic masses which differ by a factor of 10 were used, namely 7.26327 and 72.6327. The average of these two values is 39.948, which corresponds the atomic mass of argon. We refer to this crystal as the solid argon w...

    work page

  2. [2]

    R. Kubo, M. Toda, and N. Hashitsume, Statistical Physics II , 2nd ed. (Springer, Berlin, 1995) p. 155

    work page 1995

  3. [3]

    A. J. C. Ladd, B. Moran, and W. G. Hoover, Phys. Rev. B 34, 5058 (1986)

    work page 1986

  4. [4]

    J. Che, T. Cagin, W. Deng, and W. A. Goddard III, J. Chem. Phys. 113, 6888 (2000)

    work page 2000

  5. [5]

    A. J. H. McGaughey and M. Kaviany, Int. J. Heat Mass Trans. 47, 1783 (2004)

    work page 2004

  6. [6]

    A. J. H. McGaughey and M. Kaviany, Int. J. Heat Mass Trans. 47, 1799 (2004)

    work page 2004

  7. [7]

    Izvekov, P

    S. Izvekov, P . W. Chung, and B. M. Rice, Int. J. Heat Mass Trans. 54, 5623 (2011). 6

    work page 2011

  8. [9]

    Pereverzev and T

    A. Pereverzev and T. Sewell, Phys. Rev. B 97, 104308 (2018)

    work page 2018

  9. [10]

    Marcolongo, P

    A. Marcolongo, P . Umari, and S. Baroni, Nature Phys. 12, 80 (2016)

    work page 2016

  10. [11]

    Simoncelli, N

    M. Simoncelli, N. Marzari, and F. Mauri, Nat. Phys. 15, 809 (2019)

    work page 2019

  11. [12]

    Sun and P

    T. Sun and P . B. Allen, Phys. Rev. B 84, 224305 (2010)

    work page 2010

  12. [13]

    A. P . Thompson, H. M. Aktulga, R. Berger, D. S. Bolineanu, W. M. Brown, P . S. Crozier, P . J. in ’t V eld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton, Comp. Phys. Comm. 271, 108171 (2022), also see https://www.lammps.org/

    work page 2022

  13. [14]

    R. J. Hardy, Phys. Rev. 132, 168 (1963)

    work page 1963

  14. [15]

    P . B. Allen and J. L. Feldman, Phys. Rev. B 48, 12581 (1993)

    work page 1993

  15. [16]

    P . K. Schelling, S. R. Phillpot, and P . Keblinski, Phys. Rev. B 65, 144306 (2002)

    work page 2002

  16. [17]

    Marcolongo, L

    A. Marcolongo, L. Ercole, and S. Baroni, J. Chem. Theory Comput. 16, 3352 (2020)

    work page 2020

  17. [18]

    Pereverzev and T

    A. Pereverzev and T. Sewell, Crystals 10, 1123 (2020)

    work page 2020

  18. [19]

    N. W. Ashcroft and N. D. Mermin, Solid State Physics (Har- court College Publishers, Forth Worth, 1976)

    work page 1976

  19. [20]

    Touloukian, Thermophysical Properties of Matter , V ol

    Y . Touloukian, Thermophysical Properties of Matter , V ol. 2 (Plenum, New Y ork, 1970) p. 182

    work page 1970

  20. [21]

    B. W. H. van Beest, G. J. Kramer, and R. A. van Santen, Phys. Rev. Lett. 64, 1955 (1990)

    work page 1955

  21. [22]

    Pereverzev and T

    A. Pereverzev and T. Sewell, Int. J. Heat Mass Trans. 188, 122647 (2022)

    work page 2022

  22. [23]

    J. M. Ziman, Electrons and Phonons (Oxford University Press,

    work page

  23. [24]

    Pereverzev, to appear

    A. Pereverzev, to appear

    work page