pith. sign in

arxiv: 2509.05960 · v5 · pith:5VQZRVBYnew · submitted 2025-09-07 · ❄️ cond-mat.mtrl-sci

Correlation Matrix Method for Phonon Quasiparticles

Pith reviewed 2026-05-18 18:33 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords phonon anharmonicitymolecular dynamicscorrelation matricesquasiparticlestemperature-dependent phononsphonon lifetimesanharmonic effectssilicon
0
0 comments X

The pith

Anharmonic phonon modes are obtained by maximizing vibration stability when fitting atomic trajectories from molecular dynamics, with all information contained in two small correlation matrices.

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

The paper establishes a method for extracting anharmonic phonon quasiparticles directly from molecular dynamics simulations by maximizing their vibration stability during trajectory fitting. This matters because anharmonicity governs thermal transport and structural stability in real materials, yet traditional approaches often require expensive higher-order force constants. The authors prove that the full quasiparticle information resides in two compact correlation matrices S and Q that are built straight from the simulation data. An optimization scheme then delivers the temperature-dependent modes together with their frequencies and lifetimes. Tests on silicon and cubic CaSiO3 reproduce known phonon softening and other temperature trends.

Core claim

Anharmonic phonon modes can be obtained by maximizing their vibration stability during fitting the atomic trajectory. All information about these quasiparticles is contained in two small correlation matrices S and Q, which can be constructed directly from molecular dynamics simulations. Based on these matrices, an optimization scheme efficiently determines temperature-dependent phonon modes along with their frequencies and lifetimes.

What carries the argument

Two correlation matrices S and Q constructed from molecular dynamics trajectories, which encode all quasiparticle information and enable optimization by maximizing vibration stability.

If this is right

  • Temperature-dependent phonon frequencies and lifetimes can be computed from standard MD runs without explicit higher-order force constants.
  • The method reproduces the known phonon softening in cubic CaSiO3 as temperature rises.
  • Phonon quasiparticle properties become accessible for materials where full anharmonic expansions are computationally prohibitive.
  • The same matrix construction and optimization framework can be applied to other quasiparticles such as electrons, holes, and magnons.

Where Pith is reading between the lines

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

  • The approach may lower the cost of predicting thermal conductivity in strongly anharmonic solids by bypassing direct calculation of cubic and quartic terms.
  • It could be combined with existing MD packages to enable routine temperature-dependent phonon analysis during high-throughput material screening.
  • Connections to other quasiparticle renormalization schemes in condensed-matter theory might allow cross-validation of lifetimes across different formalisms.

Load-bearing premise

Maximizing vibration stability during fitting of the atomic trajectory produces the physically correct anharmonic phonon modes and lifetimes.

What would settle it

Running the optimization on silicon and obtaining phonon frequencies or lifetimes that differ from those measured by inelastic neutron scattering or computed by established anharmonic perturbation theory would falsify the claim that the maximized-stability modes are the correct quasiparticles.

Figures

Figures reproduced from arXiv: 2509.05960 by Fawei Zheng, Wenjing Li, Yong Lu.

Figure 1
Figure 1. Figure 1: FIG. 1. Optimization process and ACF for silicon phonons. (a) [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Phonon dispersions for silicon at (a) 500 K and (b) 1000 K. [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Phonon dispersions and modes for CaSiO [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
read the original abstract

Phonon anharmonicity is ubiquitous in real materials and is crucial for understanding thermal properties and phase stability. In this work, we show that anharmonic phonon modes can be obtained by maximizing their vibration stability during fitting the atomic trajectory. We prove that all information about these quasiparticles is contained in two small correlation matrices $\mathcal{S}$ and $\mathcal{Q}$, which can be constructed directly from molecular dynamics simulations. Based on these matrices, we proposed an optimization scheme, which allows us to efficiently determine temperature-dependent phonon modes along with their frequencies and lifetimes. We verified this method by applying it to silicon and cubic CaSiO$_3$, where it successfully captured their temperature-dependent phonon behaviors and the well-known phonon softening in cubic CaSiO$_3$. This theory provides a convenient tool for investigating phonon quasiparticles and can be extended to study other quasiparticles, such as electrons, holes, and magnons.

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 introduces a correlation matrix method for phonon quasiparticles. It claims that all information about anharmonic phonon modes is contained in two small correlation matrices S and Q constructed directly from molecular dynamics trajectories. An optimization scheme that maximizes vibration stability is proposed to extract temperature-dependent modes together with their frequencies and lifetimes. The method is applied to silicon and cubic CaSiO3, where it reproduces known temperature trends and the phonon softening in CaSiO3.

Significance. If validated, the approach would offer an efficient route to anharmonic phonons that avoids repeated full dynamical-matrix calculations, using only small matrices from standard MD runs. The direct construction of S and Q and the successful qualitative reproduction of established behaviors in two materials are positive features. Significance remains provisional until the optimization is shown to recover modes that match independent definitions such as spectral-function peaks.

major comments (2)
  1. [Theory section on the optimization scheme] Theory section on the optimization scheme: the central assertion that maximizing vibration stability recovers the physically correct quasiparticle modes lacks a derivation establishing that this maximum coincides with the poles of the retarded phonon Green's function or the peaks of the spectral function obtained from the same MD data. This link is load-bearing for the claim that S and Q contain all quasiparticle information.
  2. [Results section, application to Si and CaSiO3] Results section, application to Si and CaSiO3: the fitting procedure uses the identical MD trajectory to build both S/Q and to perform the stability maximization, yet no cross-validation (e.g., comparison of extracted lifetimes to those from velocity-autocorrelation Fourier transforms on held-out segments or independent simulations) is reported. This leaves open whether the modes are physical or artifacts of the chosen metric.
minor comments (2)
  1. [Notation] Notation: provide explicit matrix-element definitions for S and Q immediately after their introduction, including the precise time-window and averaging conventions used in the MD construction.
  2. [Figures] Figures: add direct overlays of literature or independent-calculation data points on the temperature-dependent frequency and lifetime plots for quantitative assessment.

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 that will be incorporated to strengthen the theoretical foundation and validation of the method.

read point-by-point responses
  1. Referee: [Theory section on the optimization scheme] Theory section on the optimization scheme: the central assertion that maximizing vibration stability recovers the physically correct quasiparticle modes lacks a derivation establishing that this maximum coincides with the poles of the retarded phonon Green's function or the peaks of the spectral function obtained from the same MD data. This link is load-bearing for the claim that S and Q contain all quasiparticle information.

    Authors: We agree that an explicit derivation connecting the stability maximization to the poles of the retarded Green's function would strengthen the presentation. The manuscript already proves that the correlation matrices S and Q constructed from MD trajectories contain all information needed to describe the anharmonic quasiparticles. The optimization is formulated to identify the modes that maximize vibration stability within this reduced description. In the revised manuscript we will add a dedicated subsection deriving that the stationary point of the stability functional coincides with the quasiparticle poles obtained from the same S and Q matrices, thereby establishing the link to the spectral function peaks. revision: yes

  2. Referee: [Results section, application to Si and CaSiO3] Results section, application to Si and CaSiO3: the fitting procedure uses the identical MD trajectory to build both S/Q and to perform the stability maximization, yet no cross-validation (e.g., comparison of extracted lifetimes to those from velocity-autocorrelation Fourier transforms on held-out segments or independent simulations) is reported. This leaves open whether the modes are physical or artifacts of the chosen metric.

    Authors: The referee correctly notes that the reported results use the full trajectory for both matrix construction and optimization. Validation in the current manuscript rests on the reproduction of established temperature trends for silicon and the well-documented phonon softening in cubic CaSiO3. To address the concern about possible metric-specific artifacts, we will add a cross-validation analysis in the revised Results section: the MD trajectory will be partitioned into independent segments, S and Q will be built from one segment, the optimization performed, and the resulting lifetimes compared against those obtained from velocity-autocorrelation Fourier transforms on the held-out segments. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper constructs S and Q directly from MD trajectories as the complete information carriers for quasiparticles, then introduces a separate optimization that maximizes vibration stability to extract modes, frequencies, and lifetimes. This procedure is presented as a new scheme and is validated against known temperature-dependent behaviors in silicon and CaSiO3. No quoted step reduces the extracted modes or the stability maximum to a tautological redefinition of the input trajectory, a fitted parameter renamed as a prediction, or a self-citation chain. The central claim therefore retains independent mathematical and computational content beyond the raw MD data.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The method rests on the domain assumption that phonon quasiparticles are well-defined objects whose properties are fully captured by the two correlation matrices and recovered by stability maximization; no free parameters or new entities are explicitly introduced in the abstract.

axioms (1)
  • domain assumption Anharmonic phonon modes can be obtained by maximizing their vibration stability during fitting the atomic trajectory.
    This premise underpins the optimization scheme and is required for the claim that the procedure yields physically meaningful quasiparticles.

pith-pipeline@v0.9.0 · 5686 in / 1223 out tokens · 36244 ms · 2026-05-18T18:33:15.689636+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.

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

62 extracted references · 62 canonical work pages

  1. [1]

    Zein, On density functional calculations of crystal elastic modula and phonon spectra, Fizika Tverdogo Tela26, 3028 (1984)

    N. Zein, On density functional calculations of crystal elastic modula and phonon spectra, Fizika Tverdogo Tela26, 3028 (1984)

  2. [2]

    Baroni, P

    S. Baroni, P. Giannozzi, and A. Testa, Green’s-function approach to linear response in solids, Phys. Rev. Lett.58, 1861 (1987)

  3. [3]

    Gonze, Adiabatic density-functional perturbation theory, Phys

    X. Gonze, Adiabatic density-functional perturbation theory, Phys. Rev. A52, 1096 (1995)

  4. [4]

    Baroni, S

    S. Baroni, S. De Gironcoli, A. Dal Corso, and P. Giannozzi, Phonons and related crystal properties from density-functional perturbation theory, Rev. Mod. Phys.73, 515 (2001)

  5. [5]

    Parlinski, Z

    K. Parlinski, Z. Li, and Y. Kawazoe, First-principles determina- tion of the soft mode in cubic ZrO 2, Phys. Rev. Lett.78, 4063 (1997)

  6. [6]

    Alf `e, PHON: A program to calculate phonons using the small displacement method, Comput

    D. Alf `e, PHON: A program to calculate phonons using the small displacement method, Comput. Phys. Commun.180, 2622 (2009)

  7. [7]

    Wang, L.-Q

    Y. Wang, L.-Q. Chen, and Z.-K. Liu, YPHON: A package for calculating phonons of polar materials, Comput. Phys. Com- mun.185, 2950 (2014)

  8. [8]

    Togo and I

    A. Togo and I. Tanaka, First principles phonon calculations in materials science, Scr. Mater.108, 1 (2015)

  9. [9]

    Tadano, Y

    T. Tadano, Y. Gohda, and S. Tsuneyuki, Anharmonic force con- stants extracted from first-principles molecular dynamics: ap- plications to heat transfer simulations, J. Phys.: Condens. Matter 26, 225402 (2014)

  10. [10]

    W. Li, J. Carrete, N. A. Katcho, and N. Mingo, ShengBTE: A solver of the Boltzmann transport equation for phonons, Com- put. Phys. Commun.185, 1747 (2014)

  11. [11]

    Chernatynskiy and S

    A. Chernatynskiy and S. R. Phillpot, Phonon transport simulator (PhonTS), Comput. Phys. Commun.192, 196 (2015)

  12. [12]

    A. Togo, L. Chaput, and I. Tanaka, Distributions of phonon lifetimes in Brillouin zones, Phys. Rev. B91, 094306 (2015)

  13. [13]

    Feng and X

    T. Feng and X. Ruan, Quantum mechanical prediction of four- phonon scattering rates and reduced thermal conductivity of solids, Phys. Rev. B93, 045202 (2016)

  14. [14]

    Carrete, B

    J. Carrete, B. Vermeersch, A. Katre, A. van Roekeghem, T. Wang, G. K. Madsen, and N. Mingo, almaBTE: A solver of the space–time dependent Boltzmann transport equation for phonons in structured materials, Comput. Phys. Commun220, 351 (2017)

  15. [15]

    Eriksson, E

    F. Eriksson, E. Fransson, and P. Erhart, The hiphive package for the extraction of high-order force constants by machine learning, Adv. Theory Simul.2, 1800184 (2019)

  16. [16]

    F. Zhou, W. Nielson, Y. Xia, and V. Ozolin ¸ˇs, Compressive sens- ing lattice dynamics. I. General formalism, Phys. Rev. B100, 184308 (2019)

  17. [17]

    Karch, F

    K. Karch, F. Bechstedt, P. Pavone, and D. Strauch, Pressure- dependent properties of SiC polytypes, Phys. Rev. B53, 13400 (1996)

  18. [18]

    D. A. Broido, M. Malorny, G. Birner, N. Mingo, and D. A. Stewart, Intrinsic lattice thermal conductivity of semiconductors from first principles, Appl. Phys. Lett.91(2007)

  19. [19]

    Lindsay, D

    L. Lindsay, D. A. Broido, and T. L. Reinecke, Thermal conduc- tivity and large isotope effect in GaN from first principles, Phys. Rev. Lett.109, 095901 (2012)

  20. [20]

    L. Xie, J. H. Feng, R. Li, and J. Q. He, First-principles study of anharmonic lattice dynamics in low thermal conductivity AgCrSe2: evidence for a large resonant four-phonon scattering, Phys. Rev. Lett.125, 245901 (2020)

  21. [21]

    Y. Xia, V. I. Hegde, K. Pal, X. Hua, D. Gaines, S. Patel, J. He, M. Aykol, and C. Wolverton, High-throughput study of lattice thermal conductivity in binary rocksalt and zinc blende com- pounds including higher-order anharmonicity, Phys. Rev. X10, 041029 (2020)

  22. [22]

    Debernardi, Phonon linewidth in III-V semiconductors from density-functional perturbation theory, Phys

    A. Debernardi, Phonon linewidth in III-V semiconductors from density-functional perturbation theory, Phys. Rev. B57, 12847 (1998)

  23. [23]

    N. K. Ravichandran and D. Broido, Phonon-phonon interactions in strongly bonded solids: Selection rules and higher-order pro- cesses, Phys. Rev. X10, 021063 (2020)

  24. [24]

    Togo, First-principles phonon calculations with phonopy and phono3py, J

    A. Togo, First-principles phonon calculations with phonopy and phono3py, J. Phys. Soc. Jpn.92, 012001 (2023)

  25. [25]

    Feng and X

    T. Feng and X. Ruan, Four-phonon scattering reduces intrinsic thermal conductivity of graphene and the contributions from flexural phonons, Phys. Rev. B97, 045202 (2018)

  26. [26]

    Tristant, A

    D. Tristant, A. Cupo, X. Ling, and V. Meunier, Phonon anhar- monicity in few-layer black phosphorus, ACS nano13, 10456 (2019)

  27. [27]

    G. Sun, J. Ma, C. Liu, Z. Xiang, D. Xu, T.-H. Liu, and X. Luo, Four-phonon and normal scattering in 2D hexagonal structures, Int. J. Heat Mass Transfer.215, 124475 (2023)

  28. [28]

    Danayat, Z

    S. Danayat, Z. T. Mona, A. S. Nayal, R. S. Annam, and J. Garg, Role of four-phonon processes in thermal conductivity of two-dimensional materials and thermal-transport enhance- ment arising from interconnected nanofiller networks in poly- mer/nanofiller composites, Nanoscale16, 13755 (2024)

  29. [29]

    Errea, M

    I. Errea, M. Calandra, and F. Mauri, Anharmonic free energies and phonon dispersions from the stochastic self-consistent har- monic approximation: Application to platinum and palladium hydrides, Phys. Rev. B89, 064302 (2014). 6

  30. [30]

    Errea, B

    I. Errea, B. Rousseau, and A. Bergara, Isotope effect in the superconducting high-pressure simple cubic phase of calcium from first principles, J. Appl. Phys.111(2012)

  31. [31]

    Monacelli, R

    L. Monacelli, R. Bianco, M. Cherubini, M. Calandra, I. Errea, and F. Mauri, The stochastic self-consistent harmonic approxi- mation: calculating vibrational properties of materials with full quantum and anharmonic effects, J. Phys.: Conden. Matter33, 363001 (2021)

  32. [32]

    T. R. Koehler, Theory of the self-consistent harmonic approx- imation with application to solid neon, Phys. Rev. Lett.17, 89 (1966)

  33. [33]

    Souvatzis, O

    P. Souvatzis, O. Eriksson, M. Katsnelson, and S. Rudin, En- tropy driven stabilization of energetically unstable crystal struc- tures explained from first principles theory, Phys. Rev. Lett.100, 095901 (2008)

  34. [34]

    Errea, M

    I. Errea, M. Calandra, and F. Mauri, Anharmonic free energies and phonon dispersions from the stochastic self-consistent har- monic approximation: Application to platinum and palladium hydrides, Phys. Rev. B89, 064302 (2014)

  35. [35]

    Tadano and S

    T. Tadano and S. Tsuneyuki, Self-consistent phonon calcula- tions of lattice dynamical properties in cubic SrTiO3 with first- principles anharmonic force constants, Phys. Rev. B92, 054301 (2015)

  36. [36]

    Simoncelli, N

    M. Simoncelli, N. Marzari, and F. Mauri, Unified theory of thermal transport in crystals and glasses, Nat. Phys.15, 809 (2019)

  37. [37]

    Simoncelli, N

    M. Simoncelli, N. Marzari, and F. Mauri, Wigner formulation of thermal transport in solids, Phys. Rev. X12, 041011 (2022)

  38. [38]

    van Roekeghem, J

    A. van Roekeghem, J. Carrete, and N. Mingo, Quantum self- consistent ab-initio lattice dynamics, Comput. Phys. Commun. 263, 107945 (2021)

  39. [39]

    L. T. Kong, Phonon dispersion measured directly from molec- ular dynamics simulations, Comput. Phys. Commun.182, 2201 (2011)

  40. [40]

    Hellman, P

    O. Hellman, P. Steneteg, I. A. Abrikosov, and S. I. Simak, Tem- perature dependent effective potential method for accurate free energy calculations of solids, Phys. Rev. B87, 104111 (2013)

  41. [41]

    J.-J. Zhou, O. Hellman, and M. Bernardi, Electron-phonon scat- tering in the presence of soft modes and electron mobility in SrTiO3 perovskite from first principles, Phys. Rev. Lett.121, 226603 (2018)

  42. [42]

    Esfarjani and H

    K. Esfarjani and H. T. Stokes, Method to extract anharmonic force constants from first principles calculations, Phys. Rev. B 77, 144112 (2008)

  43. [43]

    Esfarjani, H

    K. Esfarjani, H. Stokes, S. Nayeb Sadeghi, Y. Liang, B. Timalsina, H. Meng, J. Shiomi, B. Liao, and R. Sun, ALAT- DYN: A set of Anharmonic LATtice DYNamics codes to com- pute thermodynamic and thermal transport properties of crys- talline solids, Comput. Phys. Commun.312, 109575 (2025)

  44. [44]

    de Koker, Thermal conductivity of MgO periclase from equi- librium first principles molecular dynamics, Phys

    N. de Koker, Thermal conductivity of MgO periclase from equi- librium first principles molecular dynamics, Phys. Rev. Lett. 103, 125902 (2009)

  45. [45]

    Zhang, T

    D.-B. Zhang, T. Sun, and R. M. Wentzcovitch, Phonon quasi- particles and anharmonic free energy in complex systems, Phys. Rev. Lett.112, 058501 (2014)

  46. [46]

    Sun, D.-B

    T. Sun, D.-B. Zhang, and R. M. Wentzcovitch, Dynamic stabi- lization of cubic CaSiO 3 perovskite at high temperatures and pressures from ab initio molecular dynamics, Phys. Rev. B89, 094109 (2014)

  47. [47]

    Kresse and J

    G. Kresse and J. Furthm¨ uller, Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set, Phys. Rev. B54, 11169 (1996)

  48. [48]

    C. Wang, C. T. Chan, and K. Ho, Tight-binding molecular- dynamics study of phonon anharmonic effects in silicon and diamond, Phys. Rev. B42, 11276 (1990)

  49. [49]

    D. Kim, H. L. Smith, J. L. Niedziela, C. W. Li, D. L. Abernathy, and B. Fultz, Phonon anharmonicity in silicon from 100 to 1500 K, Phys. Rev. B91, 014307 (2015)

  50. [50]

    D. S. Kim, O. Hellman, N. Shulumba, C. N. Saunders, J. Y. Y. Lin, H. L. Smith, J. E. Herriman, J. L. Niedziela, D. L. Aber- nathy, C. W. Li, and B. Fultz, Temperature-dependent phonon lifetimes and thermal conductivity of silicon by inelastic neutron scattering and ab initio calculations, Phys. Rev. B102, 174311 (2020)

  51. [51]

    Zhang, D.-B

    Z. Zhang, D.-B. Zhang, K. Onga, A. Hasegawa, K. Ohta, K. Hi- rose, and R. M. Wentzcovitch, Thermal conductivity of CaSiO3 perovskite at lower mantle conditions, Phys. Rev. B104, 184101 (2021)

  52. [52]

    Y. Shin, E. Di Lucente, N. Marzari, and L. Monacelli, The thermodynamics of CaSiO 3 in Earth’s lower mantle, arXiv:2411.18489 (2024)

  53. [53]

    W. F. McDonough and S.-S. Sun, The composition of the Earth, Chem. Geol.120, 223 (1995)

  54. [54]

    J. J. Valencia-Cardona, G. Shukla, Z. Wu, C. Houser, D. A. Yuen, and R. M. Wentzcovitch, Influence of the iron spin crossover in ferropericlase on the lower mantle geotherm, Geo- phys. Res. Lett.44, 4863 (2017)

  55. [55]

    Komabayashi, K

    T. Komabayashi, K. Hirose, N. Sata, Y. Ohishi, and L. S. Dubrovinsky, Phase transition in CaSiO 3 perovskite, Earth Planet. Sci. Lett.260, 564 (2007)

  56. [56]

    Kurashina, K

    T. Kurashina, K. Hirose, S. Ono, N. Sata, and Y. Ohishi, Phase transition in Al-bearing CaSiO 3 perovskite: implications for seismic discontinuities in the lower mantle, Phys. Earth Planet. In.145, 67 (2004)

  57. [57]

    S. Ono, Y. Ohishi, and K. Mibe, Phase transition of Ca- perovskite and stability of Al-bearing Mg-perovskite in the lower mantle, Am. Mineral.89, 1480 (2004)

  58. [58]

    Thomson, W

    A. Thomson, W. Crichton, J. Brodholt, I. Wood, N. Siersch, J. Muir, D. Dobson, and S. A. Hunt, Seismic velocities of CaSiO3 perovskite can explain LLSVPs in Earth’s lower mantle, Nature 572, 643 (2019)

  59. [59]

    Tschauner, S

    O. Tschauner, S. Huang, S. Yang, M. Humayun, W. Liu, S. N. Gilbert Corder, H. A. Bechtel, J. Tischler, and G. R. Rossman, Discovery of davemaoite, CaSiO3-perovskite, as a mineral from the lower mantle, Science374, 891 (2021)

  60. [60]

    Immoor, L

    J. Immoor, L. Miyagi, H.-P. Liermann, S. Speziale, K. Schulze, J. Buchen, A. Kurnosov, and H. Marquardt, Weak cubic CaSiO3 perovskite in the Earth’s mantle, Nature603, 276 (2022)

  61. [61]

    Parlinski, Z

    K. Parlinski, Z. Q. Li, and Y. Kawazoe, First-principles deter- mination of the soft mode in cubic ZrO 2, Phys. Rev. Lett.78, 4063 (1997)

  62. [62]

    Y. Wang, J. J. Wang, W. Y. Wang, Z. G. Mei, S. L. Shang, L. Q. Chen, and Z. K. Liu, A mixed-space approach to first-principles calculations of phonon frequencies for polar materials, J. Phys.: Condens. Matter22, 202201 (2010)