Thermal conductivity of seifertite and pyrite-type SiO₂: A comparative study
Pith reviewed 2026-05-20 16:37 UTC · model grok-4.3
The pith
Pyrite-type SiO2 conducts heat 19 percent less than seifertite, suggesting an insulating layer inside super-Earths.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using Green-Kubo molecular dynamics with two machine learning potentials trained on SCAN and PBEsol functionals, the lattice thermal conductivity of pyrite-type SiO2 is 19 percent lower than that of seifertite at the relevant pressures and temperatures. The Green-Kubo results follow a T^{-1} temperature dependence and are substantially higher than those from the phonon quasiparticle approach, which misses diffusive phonon contributions.
What carries the argument
Green-Kubo method applied to molecular dynamics trajectories generated by machine learning potentials fitted to SCAN and PBEsol density-functional calculations.
If this is right
- Heat flow across the seifertite to pyrite-type boundary is reduced, which can steepen the temperature gradient in a planetary mantle.
- The lower-conductivity pyrite-type phase could act as a thermal blanket that slows cooling of a super-Earth core.
- The T^{-1} dependence implies that the insulating effect strengthens at the higher temperatures found deeper in large planets.
- Models of super-Earth thermal evolution must incorporate this phase-dependent drop rather than assuming uniform silica conductivity.
Where Pith is reading between the lines
- If the 19 percent drop persists under the even higher pressures of the deepest mantles, it could create a stable barrier layer that affects magnetic dynamo generation.
- Similar conductivity reductions may occur at other silica phase boundaries, altering heat transport in both Earth and exoplanet interiors.
- The difference between Green-Kubo and phonon-quasiparticle results highlights the need to include anharmonic and diffusive phonon modes when modeling transport at mantle temperatures.
Load-bearing premise
The machine learning potentials reproduce ab initio accuracy for the vibrational properties and scattering rates that control heat transport in these phases.
What would settle it
An independent first-principles calculation or high-pressure experiment that measures the conductivity ratio between the two phases and finds a difference significantly smaller or larger than 19 percent.
Figures
read the original abstract
Thermal conductivity is a fundamental material property that plays a crucial role in understanding the dynamics and evolution of planetary interiors. Despite its importance, the thermal conductivity of seifertite and pyrite-type SiO$_2$ remains unknown. Here, we calculate the lattice thermal conductivities of seifertite and pyrite-type SiO$_2$ using the Green-Kubo method based on molecular dynamics (MD) simulations driven by two machine learning potentials (MLPs) constructed from the SCAN and PBEsol exchange-correlation functionals, with $\textit{ab initio}$-level accuracy. To demonstrate our methodology, we also compute thermal conductivities using the phonon quasiparticle approach for comparison. Overall, the Green-Kubo method predicts up to 119 % higher thermal conductivity with a temperature dependence close to $T^{-1}$, as it fully captures diffusion-like phonons at high temperatures that are missed by the phonon quasiparticle approach. The 19 % reduction in thermal conductivity across the phase transition from seifertite to the pyrite-type phase suggests the potential formation of a thermally insulating layer in the mantle of super-Earths.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript calculates the lattice thermal conductivities of seifertite and pyrite-type SiO₂ via the Green-Kubo method applied to MD trajectories generated by two machine-learning potentials (MLPs) trained on SCAN and PBEsol functionals. It compares these results to the phonon quasiparticle approach, reports that Green-Kubo yields up to 119% higher values with a T^{-1} temperature dependence due to capture of diffusive modes, finds a 19% reduction across the seifertite-to-pyrite transition, and infers possible formation of a thermally insulating layer in super-Earth mantles.
Significance. If the MLPs deliver the claimed ab initio accuracy for anharmonic phonon scattering and diffusive modes under the relevant high-P, high-T conditions, the work would supply useful constraints on heat transport in high-pressure silica phases relevant to planetary interiors. The explicit comparison between Green-Kubo and quasiparticle methods usefully illustrates the limitations of the latter at elevated temperatures.
major comments (2)
- [Computational Methods] Computational Methods section: The manuscript asserts 'ab initio-level accuracy' for the SCAN- and PBEsol-trained MLPs but provides no direct cross-validation of the resulting Green-Kubo thermal conductivities (or underlying force constants) against explicit DFT-based MD at the simulated pressures and temperatures; given the sensitivity of lattice thermal conductivity to anharmonic scattering, this omission leaves both the absolute values and the reported 19% reduction without demonstrated support.
- [Results] Results section: No error bars, statistical uncertainties, or convergence tests with respect to simulation cell size or trajectory length are reported for the Green-Kubo conductivities, despite the abstract quoting precise figures (119% higher, 19% reduction) that are central to the planetary implication.
minor comments (2)
- [Abstract] Abstract: The statement 'up to 119 % higher' does not specify the temperature, pressure, or phase for which the maximum difference occurs.
- [Computational Methods] The manuscript would benefit from a brief table comparing key MLP validation metrics (force RMSE, energy RMSE) against the underlying DFT data for both functionals.
Simulated Author's Rebuttal
We thank the referee for their thorough review and constructive feedback on our manuscript. We appreciate the recognition of the potential implications for planetary interiors. Below, we address each of the major comments in detail.
read point-by-point responses
-
Referee: [Computational Methods] Computational Methods section: The manuscript asserts 'ab initio-level accuracy' for the SCAN- and PBEsol-trained MLPs but provides no direct cross-validation of the resulting Green-Kubo thermal conductivities (or underlying force constants) against explicit DFT-based MD at the simulated pressures and temperatures; given the sensitivity of lattice thermal conductivity to anharmonic scattering, this omission leaves both the absolute values and the reported 19% reduction without demonstrated support.
Authors: We acknowledge the referee's concern regarding the lack of direct validation of the Green-Kubo results against DFT MD. Direct ab initio MD simulations for thermal conductivity calculations are computationally prohibitive for the large system sizes and long trajectories needed for convergence in the Green-Kubo formalism. Our MLPs were trained on extensive DFT datasets and validated for structural, energetic, and vibrational properties under high-pressure and high-temperature conditions relevant to the study. We will revise the manuscript to include additional details on the MLP validation, such as force and energy errors, and discuss how these ensure accuracy for anharmonic phonon scattering. We also note that the comparison with the phonon quasiparticle method provides indirect support for the trends observed. revision: yes
-
Referee: [Results] Results section: No error bars, statistical uncertainties, or convergence tests with respect to simulation cell size or trajectory length are reported for the Green-Kubo conductivities, despite the abstract quoting precise figures (119% higher, 19% reduction) that are central to the planetary implication.
Authors: We agree that providing statistical uncertainties and convergence tests is crucial for the reliability of the reported values. In the revised version of the manuscript, we will add error bars derived from the standard deviation across multiple independent MD runs or block averaging techniques. Additionally, we will include convergence tests demonstrating the stability of the thermal conductivity values with respect to simulation cell size and trajectory length. revision: yes
Circularity Check
No circularity: results derived from external functionals via explicit MD
full rationale
The paper computes lattice thermal conductivities of seifertite and pyrite-type SiO2 via Green-Kubo MD driven by MLPs trained on standard SCAN and PBEsol functionals, which are external ab initio methods. The reported 19% reduction and absolute values emerge from explicit simulation trajectories rather than any internal fit, self-definition, or self-citation chain that reduces the output to the input by construction. The phonon quasiparticle comparison is presented as an independent methodological check. No load-bearing step invokes a uniqueness theorem, ansatz smuggled via prior work, or renaming of known results; the derivation remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Machine learning potentials trained on SCAN and PBEsol exchange-correlation functionals achieve ab initio-level accuracy for lattice thermal conductivity in high-pressure SiO2 phases.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
we calculate the lattice thermal conductivities of seifertite and pyrite-type SiO2 using the Green-Kubo method based on molecular dynamics (MD) simulations driven by two machine learning potentials (MLPs) constructed from the SCAN and PBEsol exchange-correlation functionals, with ab initio-level accuracy
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
-
[1]
D. G. Cahill, W. K. Ford, K. E. Goodson, G. D. Ma- han, A. Majumdar, H. J. Maris, R. Merlin, and S. R. Phillpot, Nanoscale thermal transport, J. Appl. Phys.93, 793 (2003)
work page 2003
-
[2]
L. T. Elkins-Tanton, Magma oceans in the inner solar system, Annu. Rev. Earth Planet. Sci.40, 113 (2012)
work page 2012
-
[3]
B. A. Buffett, Estimates of heat flow in the deep man- tle based on the power requirements for the geodynamo, Geophys. Res. Lett.29, 7 (2002)
work page 2002
-
[4]
J. B. Naliboff and L. H. Kellogg, Dynamic effects of a step-wise increase in thermal conductivity and viscosity in the lowermost mantle, Geophys. Res. Lett.33, L12301 (2006)
work page 2006
-
[5]
K. Umemoto, R. M. Wentzcovitch, S. Wu, M. Ji, C.- Z. Wang, and K.-M. Ho, Phase transitions in mgsio 3 post-perovskite in super-earth mantles, Earth Planet. Sci. Lett.478, 40 (2017)
work page 2017
-
[6]
X. Deng, J. Deng, M. J. Walter, Y. Fei, J. Yang, and R. E. Cohen, Phase relations and melting of sio 2 under mantle conditions determined from machine learning potentials, Geophys. Res. Lett.53, e2025GL119400 (2026)
work page 2026
-
[7]
D. Park, X. Deng, and J. Deng, Melting phase rela- tion of seifertite and pyrite-type sio 2 determined by ma- chine learning potentials, Phys. Rev. B 10.1103/b44t- 2x7s (2026)
-
[8]
Y. Kang and K. Morita, Thermal conductivity of the cao–al2o3–sio2 system, ISIJ Int.46, 420 (2006)
work page 2006
-
[9]
H. Hasegawa, Y. Hoshino, T. Kasamoto, Y. Akaida, T. Kowatari, Y. Shiroki, H. Shibata, H. Ohta, and Y. Waseda, Thermal conductivity measurements of some synthetic al2o3–cao–sio2 slags by means of a front-heating and front-detection laser-flash method, Metall. Mater. Trans. B43, 1405 (2012)
work page 2012
-
[10]
P. G. Klemens, Thermal conductivity and lattice vibra- tional modes, inSolid State Physics, Vol. 7, edited by F. Seitz and D. Turnbull (Elsevier, 1958) pp. 1–98
work page 1958
-
[11]
F. M¨ uller-Plathe, A simple nonequilibrium molecular dy- namics method for calculating the thermal conductivity, J. Chem. Phys.106, 6082 (1997)
work page 1997
-
[12]
M. S. Green, Markoff random processes and the sta- tistical mechanics of time-dependent phenomena. ii. ir- reversible processes in fluids, J. Chem. Phys.22, 398 (1954)
work page 1954
-
[13]
Kubo, Statistical-mechanical theory of irreversible processes
R. Kubo, Statistical-mechanical theory of irreversible processes. i. general theory and simple applications to magnetic and conduction problems, J. Phys. Soc. Jpn. 12, 570 (1957)
work page 1957
-
[14]
X. Tang and J. Dong, Lattice thermal conductivity of mgo at conditions of earth’s interior, Proc. Natl. Acad. Sci. U.S.A.107, 4539 (2010)
work page 2010
-
[15]
H. Dekura and T. Tsuchiya, Lattice thermal conductiv- ity of mgsio3 postperovskite under the lowermost mantle conditions from ab initio anharmonic lattice dynamics, Geophys. Res. Lett.46, 12919 (2019)
work page 2019
-
[16]
Z. Zhang, D.-B. Zhang, K. Onga, A. Hasegawa, K. Ohta, K. Hirose, and R. M. Wentzcovitch, Thermal conductiv- ity of casio3 perovskite at lower mantle conditions, Phys. Rev. B104, 184101 (2021)
work page 2021
-
[17]
G. Barbalinardo, Z. Chen, N. W. Lundgren, and D. Donadio, Efficient anharmonic lattice dynamics calcu- lations of thermal transport in crystalline and disordered solids, J. Appl. Phys.128, 135104 (2020)
work page 2020
-
[18]
D. P. Sellan, E. S. Landry, J. E. Turney, A. J. H. Mc- Gaughey, and C. H. Amon, Size effects in molecular dy- namics thermal conductivity predictions, Phys. Rev. B 81, 214305 (2010)
work page 2010
-
[19]
S. Stackhouse, L. Stixrude, and B. B. Karki, Ther- mal conductivity of periclase (mgo) from first principles, Phys. Rev. Lett.104, 208501 (2010)
work page 2010
-
[20]
P. K. Schelling, S. R. Phillpot, and P. Keblinski, Com- parison of atomic-level simulation methods for computing thermal conductivity, Phys. Rev. B65, 144306 (2002)
work page 2002
-
[21]
A. Marcolongo, P. Umari, and S. Baroni, Microscopic theory and quantum simulation of atomic heat transport, Nat. Phys.12, 80 (2016)
work page 2016
-
[22]
C. Carbogno, R. Ramprasad, and M. Scheffler, Ab ini- tio green–kubo approach for the thermal conductivity of solids, Phys. Rev. Lett.118, 175901 (2017)
work page 2017
-
[23]
J. Behler and M. Parrinello, Generalized neural-network representation of high-dimensional potential-energy sur- faces, Phys. Rev. Lett.98, 146401 (2007)
work page 2007
- [24]
-
[25]
J. Deng and L. Stixrude, Thermal conductivity of sili- cate liquid determined by machine learning potentials, Geophys. Res. Lett.48, e2021GL093806 (2021)
work page 2021
-
[26]
J. P. Perdew, A. Ruzsinszky, G. I. Csonka, O. A. Vydrov, G. E. Scuseria, L. A. Constantin, X. Zhou, and K. Burke, Restoring the density-gradient expansion for exchange in solids and surfaces, Phys. Rev. Lett.100, 136406 (2008)
work page 2008
-
[27]
J. Sun, A. Ruzsinszky, and J. P. Perdew, Strongly con- strained and appropriately normed semilocal density functional, Phys. Rev. Lett.115, 036402 (2015)
work page 2015
-
[28]
Z. Zhang, D.-B. Zhang, T. Sun, and R. M. Wentzcov- itch, phq: A fortran code to compute phonon quasiparti- cle properties and dispersions, Comput. Phys. Commun. 243, 110 (2019)
work page 2019
- [29]
- [30]
-
[31]
J. Zeng, D. Zhang, A. Peng, X. Zhang, S. He, Y. Wang, X. Liu, H. Bi, Y. Li, C. Cai, C. Zhang, Y. Du, J.-X. Zhu, P. Mo, Z. Huang, Q. Zeng, S. Shi, X. Qin, Z. Yu, and H. Wang, Deepmd-kit v3: A multiple-backend frame- work for machine learning potentials, J. Chem. Theory Comput.21, 4375 (2025)
work page 2025
-
[32]
J. Deng, H. Niu, J. Hu, M. Chen, and L. Stixrude, Melt- ing of mgsio3 determined by machine learning potentials, Phys. Rev. B107, 064103 (2023)
work page 2023
-
[33]
Plimpton, Fast parallel algorithms for short-range molecular dynamics, J
S. Plimpton, Fast parallel algorithms for short-range molecular dynamics, J. Comput. Phys.117, 1 (1995)
work page 1995
-
[34]
A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolin- tineanu, W. M. Brown, P. S. Crozier, P. J. in ’t Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton, Lammps—a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales, Comput. Phys. Commun.27...
work page 2022
-
[35]
W. G. Hoover, Canonical dynamics: Equilibrium phase- space distributions, Phys. Rev. A31, 1695 (1985)
work page 1985
- [36]
-
[37]
French, Thermal conductivity of dissociating water— an ab initio study, New J
M. French, Thermal conductivity of dissociating water— an ab initio study, New J. Phys.21, 023007 (2019)
work page 2019
-
[38]
See the Supplemental Material at for details of the con- vergence tests for the methods and parameters
- [39]
-
[40]
D. Surblys, H. Matsubara, G. Kikugawa, and T. Ohara, Application of atomic stress to compute heat flux via molecular dynamics for systems with many-body inter- actions, Phys. Rev. E99, 051301 (2019)
work page 2019
-
[41]
D. Surblys, H. Matsubara, G. Kikugawa, and T. Ohara, Methodology and meaning of computing heat flux via atomic stress in systems with constraint dynamics, J. Appl. Phys.130, 215104 (2021)
work page 2021
-
[42]
H. R. Seyf, K. Gordiz, F. DeAngelis, and A. Henry, Using green–kubo modal analysis (gkma) and interface conduc- tance modal analysis (icma) to study phonon transport with molecular dynamics, J. Appl. Phys.125, 081101 (2019)
work page 2019
-
[43]
R. J. Hardy, Energy-flux operator for a lattice, Phys. Rev. 132, 168 (1963)
work page 1963
-
[44]
A. Togo and I. Tanaka, First principles phonon calcula- tions in materials science, Scr. Mater.108, 1 (2015)
work page 2015
-
[45]
P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococ- cioni, I. Dabo, A. D. Corso, S. de Gironcoli, S. Fab- ris, G. Fratesi, R. Gebauer, U. Gerstmann, C. Gougous- sis, A. Kokalj, M. Lazzeri, and R. M. Wentzcovitch, Quantum espresso: A modular and open-source software project for quantum simulations ...
work page 2009
-
[46]
G. M. Manthilake, N. de Koker, D. J. Frost, and C. A. McCammon, Lattice thermal conductivity of lower man- tle minerals and heat flux from earth’s core, Proc. Natl. Acad. Sci. U.S.A.108, 17901 (2011)
work page 2011
- [47]
-
[48]
J. L. Feldman, M. D. Kluge, P. B. Allen, and F. Wooten, Thermal conductivity and localization in glasses: Numer- ical study of a model of amorphous silicon, Phys. Rev. B 48, 12589 (1993)
work page 1993
- [49]
-
[50]
M. Simoncelli, N. Marzari, and F. Mauri, Wigner formu- lation of thermal transport in solids, Phys. Rev. X12, 041011 (2022)
work page 2022
-
[51]
H. R. Seyf and A. Henry, A method for distinguishing between propagons, diffusions, and locons, J. Appl. Phys. 120, 025101 (2016)
work page 2016
-
[52]
W. Li, J. Carrete, G. K. H. Madsen, and N. Mingo, In- fluence of the optical-acoustic phonon hybridization on phonon scattering and thermal conductivity, Phys. Rev. B93, 205203 (2016)
work page 2016
-
[53]
T. Jia, X. Liu, Y. Zhang, and S.-H. Wei, Origin of the variation in lattice thermal conductivities in pyrite-type dichalcogenides, Phys. Rev. B107, 115204 (2023)
work page 2023
-
[54]
S. Mukhopadhyay, D. S. Parker, B. C. Sales, A. A. Puret- zky, M. A. McGuire, and L. Lindsay, Two-channel model for ultralow thermal conductivity of crystalline tl 3vse4, Science360, 1455 (2018)
work page 2018
-
[55]
D. Wang, Z. Wu, and X. Deng, Thermal conductivity of fe-bearing bridgmanite and post-perovskite: Implications for the heat flux from the core, Earth Planet. Sci. Lett. 621, 118368 (2023)
work page 2023
-
[56]
Y. Peng and J. Deng, Thermal conductivity of mgsio 3– h2o system determined by machine learning potentials, Geophys. Res. Lett.51, e2023GL107245 (2024)
work page 2024
-
[57]
Y. Xu, T. J. Shankland, S. Linhardt, D. C. Rubie, F. Lan- genhorst, and K. Klasinski, Thermal diffusivity and con- ductivity of olivine, wadsleyite and ringwoodite to 20 gpa and 1373 k, Phys. Earth Planet. Inter.143–144, 321 (2004)
work page 2004
-
[58]
S. A. Hunt, A. M. Walker, R. J. McCormack, D. P. Dob- son, A. S. Wills, and L. Li, The effect of pressure on thermal diffusivity in pyroxenes, Mineral. Mag.75, 2597 (2011)
work page 2011
-
[59]
F. Deschamps and W.-P. Hsieh, Lowermost mantle ther- mal conductivity constrained from experimental data and tomographic models, Geophys. J. Int.219, S115 (2019)
work page 2019
-
[60]
D. Wang, Z. Wu, and X. Deng, Thermal conductivity of hydrous wadsleyite determined by nonequilibrium molec- ular dynamics based on machine learning, Geophys. Res. Lett.49, e2022GL100337 (2022)
work page 2022
- [61]
-
[62]
J. Wei, Z. Xia, Y. Xia, and J. He, Hierarchy of exchange- correlation functionals in computing lattice thermal con- ductivities of rocksalt and zinc-blende semiconductors, Phys. Rev. B110, 035205 (2024)
work page 2024
-
[63]
J. M. Guerrero, F. Deschamps, Y. Li, W.-P. Hsieh, and P. J. Tackley, Influence of heterogeneous thermal con- ductivity on the long-term evolution of the lower-mantle thermochemical structure: Implications for primordial reservoirs, Solid Earth14, 119 (2023)
work page 2023
-
[64]
Christensen, Effects of phase transitions on mantle convection, Annu
U. Christensen, Effects of phase transitions on mantle convection, Annu. Rev. Earth Planet. Sci.23, 65 (1995)
work page 1995
-
[65]
U. Christensen, Geodynamo models with a stable layer and heterogeneous heat flow at the top of the core, Geo- 10 phys. J. Int.215, 1338 (2018)
work page 2018
-
[66]
T. V. Hoolst, L. Noack, and A. Rivoldini, Exoplanet in- teriors and habitability, Adv. Phys. X4, 1630316 (2019). [67]https://osf.io/y9qnmwith DOI: 10.17605/OSF.IO/Y9QNM. Thermal conductivity of seifertite and pyrite-type SiO 2: A comparative study Doyoon Park ,1,∗ Yihang Peng ,2 and Jie Deng 2,† 1Department of Mechanical and Aerospace Engineering, Princet...
-
[67]
H. R. Seyf and A. Henry, A method for distinguishing between propagons, diffusions, and locons, J. Appl. Phys.120, 025101 (2016)
work page 2016
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.