Pith. sign in

REVIEW 3 major objections 5 minor 94 references

Consistency between the Green-Kubo formula and Lorentz model for predicting the infrared dielectric function of polar materials

T0 review · 3 major / 5 minor · reviewed 2026-08-16 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The Green-Kubo formula and the Lorentz model give unified predictions of infrared dielectric functions when the Lorentz model is parameterized with a frequency-dependent phonon self-energy from molecular dynamics and electronic…

desk verdict A useful, honest benchmark paper showing that a Lorentz model with an MD-derived frequency-dependent phonon self-energy reproduces the Green-Kubo dielectric function, plus a real ε∞ correction for rigid-ion models; the main caveat is that both sides come from the same classical trajectories. read the letter →

arxiv 2504.17464 v1 pith:NGZRHLB3 submitted 2025-04-24 cond-mat.mes-hall

Wei-Zhe Yuan , Yangyu Guo , Hong-Liang Yi This is my paper · ORCID
classification cond-mat.mes-hall
keywords infrareddielectricfunctionGreen-KuboformulaLorentzmodelphononself-energymulti-phononabsorptionmoleculardynamicspolarmaterialsBorneffectivecharge
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper argues that two standard routes to the infrared dielectric function of polar crystals—the Green-Kubo formula computed from equilibrium molecular dynamics and the phenomenological Lorentz model—can be made quantitatively consistent. It shows that the usual Lorentz model with a constant linewidth misses the multi-phonon absorption that the Green-Kubo correlation function contains, and that replacing the linewidth by a frequency-dependent phonon self-energy removes the discrepancy. It also identifies how electronic polarization must enter each approach: a multiplicative $\varepsilon_\infty$ factor for rigid-ion models, and Born effective charges for machine-learned potentials. Demonstrations on MgO and LiH support the claim, with classical molecular dynamics accurate at high temperature and perturbation theory better at cryogenic temperature. If the reconciliation is right, molecular dynamics becomes a dependable predictor of temperature-dependent infrared spectra, including multiphonon features.

What carries the argument

The central object is the frequency-dependent phonon self-energy $\Pi(\omega)$ of the infrared-active transverse-optical phonon, which generalizes the constant linewidth of the textbook Lorentz model. It is extracted from equilibrium molecular dynamics by projecting atomic velocities onto phonon modes, constructing the retarded single-phonon Green's function through the Kubo transform, and inverting Dyson's equation; the imaginary part gives the frequency-dependent damping and the real part the anharmonic frequency shift. Inserting $\Pi(\omega)$ into the generalized Lorentz model (Eq. 8) is what makes the dielectric function carry multi-phonon absorption. The second piece of machinery is the electronic-polarization correction: a derivation from the standard oscillator equations shows that rigid-ion models must multiply the ionic dielectric function by $\varepsilon_\infty$ rather than add $\varepsilon_\infty-1$, while machine-learned potentials supply the electronic response implicitly through Born effective charges when the dipole moment is evaluated.

What would settle it

Repeat the two-route comparison on a third strongly anharmonic polar crystal at a temperature where classical statistics should hold—say 300 K for a heavier-ion oxide—using the same machine-learned potential; if the Lorentz model fed with the MD-derived self-energy does not match the Green-Kubo dielectric function wherever multi-phonon features appear, the claimed reconciliation is not general.

Watch

Extended reading notes

Core claim

The central claim is that the Green-Kubo formula and the Lorentz model are the same physics expressed twice, once through dipole-moment fluctuations and once through damped oscillator response, and that they agree quantitatively when the Lorentz model is fed the full frequency-dependent phonon self-energy $\Pi(\omega)=\Delta(\omega)-i\Gamma(\omega)$ extracted from MD velocity-correlation functions via the Kubo transform and Dyson's equation. With this generalized damping, the Lorentz model reproduces the multi-phonon absorption structure of the Green-Kubo dielectric function for MgO and LiH, which a constant linewidth $\tau^{-1}$ cannot do. The paper further claims that electronic polarization is not a simple additive background: for a rigid-ion model the correct correction is multiplicative, $\varepsilon=\varepsilon_\infty\varepsilon_{\mathrm{ion}}$, derived from the standard ionic-polarization equations, while a machine-learned potential already encodes the electron-polarization response in its atomic dynamics and needs only Born effective charges to compute the dipole moment. Classical molecular dynamics with a machine-learned potential achieves quantitative agreement with experiment at 295 K and 950 K for MgO but fails at 8 K, where perturbation theory and path-integral methods perform better.

Load-bearing premise

The argument holds only if the thermal motion of atoms in a classical molecular dynamics simulation reproduces the same anharmonic vibrations that cause real infrared absorption, so that the damping pulled out of simulated velocity correlations faithfully represents the crystal's actual multiphonon processes.

Editorial extensions

If this is right

  • Molecular dynamics and the Lorentz model can now be cross-validated: either route yields the same infrared dielectric function when the Lorentz model uses an MD-derived frequency-dependent self-energy, so results from one method can be checked against the other.
  • The conventional rigid-ion correction $\varepsilon_{\mathrm{GK}}+\varepsilon_\infty-1$ underestimates the ionic infrared response; the physical correction is $\varepsilon_{\mathrm{GK}}$ multiplied by $\varepsilon_\infty$, which changes the predicted longitudinal-optical frequency and the reflectance outside the main reflection band.
  • Machine-learned potentials trained on DFT data, combined with Born effective charges, reproduce the measured infrared reflectance of MgO at room and elevated temperatures, including multi-phonon absorption.
  • Classical molecular dynamics is reliable for infrared spectra at elevated temperatures but not at cryogenic temperatures, where quantum nuclear effects and isotope scattering dominate; perturbation theory handles that regime better, and path-integral molecular dynamics partially restores the agreement.
  • The generalized Lorentz model with a frequency-dependent self-energy, not a constant linewidth, is the correct phenomenological target for atomistic predictions of temperature-dependent infrared optical properties.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The same self-energy route should extend to other infrared-active crystals and to anisotropic or two-dimensional polar materials, where multi-phonon features and frequency-dependent damping are even more prominent; a test on such a material would sharpen the generality of the reconciliation.
  • The multiplicative $\varepsilon_\infty$ correction implies that the static dielectric constant predicted by rigid-ion molecular dynamics should be multiplied by $\varepsilon_\infty$, a prediction that can be checked directly against measured $\varepsilon(0)$ across a series of alkali halides.
  • If the equivalence is exact in the classical regime, the phonon self-energy from molecular dynamics can be used as a temperature-dependent input to mesoscale radiative-heat-transfer models, giving them an atomistic grounding without running MD at every design point.
  • The observed failure of the additive correction suggests prior rigid-ion molecular-dynamics studies of near-field radiative heat transfer may have systematically underestimated infrared responses; recomputing those quantities with the multiplicative correction is a concrete test of the paper's mechanism.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The manuscript compares two routes to the infrared dielectric function of polar crystals, MgO and LiH: the Green–Kubo formula evaluated from equilibrium molecular dynamics and the Lorentz model parameterized either by a constant TO linewidth from spectral energy density analysis or by a frequency-dependent phonon self-energy extracted from modal velocity correlation functions. Using both a rigid-ion BKS potential and a NEP machine-learned potential, it shows that the simple Lorentz model misses multi-phonon absorption features, while the Lorentz model with an MD-derived self-energy reproduces the Green–Kubo spectra. It further argues that, with the rigid-ion model, the electronic polarization correction is multiplicative through the factor ε∞ (Eq. 28), whereas with the NEP the electronic contribution is captured by assigning Born effective charges when computing the dipole moment. The paper also compares with experiment and perturbation theory and includes a path-integral MD analysis of nuclear quantum effects.

Significance. If the claims hold, the paper supplies a useful practical recipe: MD-derived phonon self-energies can be injected into a Lorentz-form dielectric function to capture anharmonic and multi-phonon infrared absorption, and the ε∞ correction for rigid-ion MD clarifies an often approximate treatment. Strengths include the two-material/two-potential benchmark, the comparison against experimental reflectance and perturbation theory in Fig. 4, the force-error sensitivity test in Appendix C, and openly available NEP training data. The principal caveat is that the Green–Kubo and Lorentz-plus-self-energy curves are derived from the same MD trajectories, so their agreement is a consistency check between two reduction formulas rather than independent cross-validation; the manuscript should state this explicitly and qualify the regime of validity of the claimed unification.

major comments (3)
  1. [Section III C and Fig. 3] The Lorentz model with Π(ω) obtained from Eqs. (9)–(15) and the Green–Kubo susceptibility of Eq. (1) are both computed from the same equilibrium MD trajectories. The close agreement in Figs. 3(a)–3(f) therefore demonstrates internal consistency of two spectral reductions, not an independent cross-validation of the underlying dynamics. The Introduction's statement that 'the Green-Kubo formula reconciles the Lorentz model parameterized with phonon self-energy extracted from MD simulation' and the Abstract's 'cross-validation' wording should be tempered: the independent tests are the experimental reflectance and the perturbation-theory comparison, while the GK–Lorentz agreement tests the single-mode parametrization. I recommend adding an explicit statement of this distinction and, if feasible, an analytical demonstration of how Eq. (8) with Π from Eq. (14) reduces to the same spectral function as Eq. (1) for a single infrared-active mode.
  2. [Section II B, Eq. (9), and Section III D, Fig. 4] The extraction of Γ(ω) uses the classical Kubo transform, so the phonon self-energy inherits classical statistics: there is no zero-point motion and no isotope-disorder scattering. The authors' own Fig. 4 shows that the MD-NEP reflectance at 8 K deviates from experiment while perturbation theory is much closer, and the PIMD result is propagated only at the Green–Kubo level, not through the self-energy/Lorentz construction. Consequently, the reconciliation is demonstrated only in the regime of classical nuclear statistics and moderate anharmonicity. This regime restriction should appear near the Introduction's central claim and in the Abstract, not only in the Conclusion.
  3. [Section II D 2] For the NEP calculations, fixed DFPT Born effective charges are assigned only when computing the dipole moment, while the underlying potential has no explicit charge or polarization degrees of freedom and was trained on neutral DFT cells. The paper correctly notes that the mechanism remains to be fully elucidated, but this approximation is load-bearing for the MLP branch of the unification claim. A quantitative test would strengthen the argument, such as comparing the ε∞εion result obtained from the RIM charges with the Z*-based NEP result, or checking against density-functional perturbation theory mode effective charges. At minimum, the statement that MLP 'automatically' includes the electronic contribution should be softened to an effective description in which the fitted forces absorb polarization effects while the dipole is assigned via fixed Born effective charges.
minor comments (5)
  1. [Fig. 3 caption] The caption contains a duplicate panel label '(d)' for the MgO reflectance panel; the MgO panels should be labeled (a)–(c) and the LiH panels (d)–(f).
  2. [Eq. (6)] The summation notation '3,n' in the SED expression is unclear; the index ranges should be written explicitly, for example as a sum over the 3 Cartesian components and the n basis atoms.
  3. [Section II D] The sentence contains the typo 'simulation detals'; also, the statement that the 10×10×10 supercell 'has been tested to be sufficiently capture the long-range interaction well' should be rephrased.
  4. [Section III C] A brief statement on statistical uncertainty or convergence of the MD spectra with respect to trajectory length and supercell size would help the reader assess how robust the agreement in Fig. 3 is.
  5. [Section III D] The text says PIMD brings the 8 K MgO reflectance closer to experiment, but it should state explicitly at which temperature the PIMD curve in Fig. 4 is computed and whether PIMD data at 295 K and 950 K are also shown or only the 8 K case.

Circularity Check

1 steps flagged · score 4.0 of 10

The Lorentz/GK 'unification' is a same-trajectory identity rather than an independent cross-validation, though external experimental benchmarks keep the paper's spectral predictions independently supported.

  1. other [Section II.B (Eqs. 8-14) and Section III.C (Fig. 3)]
    "By utilizing the phonon self-energy of MgO and LiH extracted from MD, we calculate the dielectric function according to Eq. 8 and derive the corresponding reflectance, as presented in Fig. 3. The calculated results demonstrate excellent consistency with those predicted by the Green-Kubo formula."

    The Green-Kubo susceptibility in Eq. (1b) is computed from the current autocorrelation, with J(t) = sum q_b v_lb(t) from Eq. (3). The mode-projected velocities v_k,nu(t) in Eq. (11) are the same linear combinations of the atomic velocities, so for the infrared-active TO mode the current correlation is proportional to <v* v>_omega. Equation (9) defines the retarded Green's function from exactly this <v* v>_omega, Eq. (14) extracts Gamma(omega) from that Green's function, and Eq. (8) then builds the Lorentz dielectric function from Gamma(omega). Thus the Lorentz-with-self-energy spectrum is a rearrangement of the same modal velocity autocorrelation that already determines the Green-Kubo spectrum. The excellent consistency in Fig.

full rationale

No load-bearing self-citation was found; the authors' prior work appears only as background citation. The paper is otherwise honest about its limitations: the 8 K classical-MD failure is openly attributed to isotope disorder and nuclear quantum effects, and the PIMD comparison is performed only at the Green-Kubo level and is never propagated through the self-energy/Lorentz construction, so the claimed unification is not extended to the quantum regime. The one genuine circularity concern is the internal 'unification' claim: because the phonon self-energy is extracted from the same modal velocity correlation function that, through the current correlation, already determines the Green-Kubo susceptibility, the Lorentz model with MD-derived self-energy is not an independent route. For a single infrared-active TO mode, Eq. (1b) and Eqs. (9)-(14) are two representations of the same correlation function, so the Fig. 3 overlap is a consistency check rather than a cross-validation. However, the paper's practical spectral predictions are independently supported by reflectance measurements and by the comparison of the MD self-energy with experimental fits, so the circularity is partial and confined to the 'reconciliation' framing, yielding a moderate score rather than a high one.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

The central claims rest on standard linear response theory (fluctuation-dissipation, Kubo transform, Dyson equation) plus a set of modeling choices: classical statistics for MD, a 1D Born-Huang chain with local field for the ε∞ correction, and the use of DFT/DFPT-fitted potentials and charges. Each of these is either a standard domain assumption or a paper-specific modeling assumption that limits the range of validity (e.g., the 8 K failure).

free parameters (4)
  • Born effective charges Z* (MgO=1.95, LiH=1.02) = 1.95 (MgO), 1.02 (LiH)
    Chosen from DFPT, used to compute dipole moments in MLP-based Green-Kubo; affects the magnitude of infrared absorption and static dielectric constant; the paper notes DFPT overestimates LiH's effective charge.
  • Experimental high-frequency dielectric constants ε∞ (MgO=3.01, LiH=3.61) = 3.01 (MgO), 3.61 (LiH)
    Adopted instead of DFPT values (3.16 and 4.72) to match optical measurements; directly controls the ε∞ε_GK correction and reflectance outside the Reststrahlen band.
  • LiH BKS potential parameters (9 potential parameters + 2 partial charges) = charges ±0.46 e; parameters in Table II
    Newly fitted in this work to DFT phonon dispersion; the resulting RIM fails to reproduce anharmonicity and experimental reflectance, which the authors attribute to the two-body form.
  • SED Lorentzian fit parameters C, ω, Γ (Eq. 7) = not tabulated
    Phonon frequencies and linewidths used in the standard Lorentz model are obtained by fitting SED spectra with Lorentzian functions; these are model fits to MD data, not to dielectric data.
assumptions (5)
  • domain assumption Fluctuation-dissipation theorem and the Green-Kubo relation between dipole moment/current fluctuations and dielectric susceptibility (Eq. 1)
    Foundation of the GK dielectric function; standard linear response theory for classical systems.
  • domain assumption The retarded phonon Green's function computed from classical modal velocity correlations via the Kubo transform (Eq. 9) equals the phonon Green's function entering the anharmonic Lorentz model
    Connects MD correlation functions to self-energy; assumes classical statistics adequately represent quantum phonon populations at the temperatures of interest; the paper's 8 K MgO result shows this can fail.
  • ad hoc to paper Born-Huang 1D diatomic chain with local field E_eff = E + P/(3ε0) captures the electronic polarization correction for 3D cubic crystals, yielding ε = ε∞ ε_ion (Eq. 28)
    Used to derive the multiplicative ε∞ correction; neglects frequency dependence of electronic screening and 3D specifics; validated indirectly by comparison with the Lorentz model.
  • domain assumption Lyddane-Sachs-Teller relation and the standard Lorentz oscillator form (Eq. 5) apply to the anharmonic crystal with damping
    Used to rewrite Eq. 26 into Eq. 27 and to define the linewidth-parameterized Lorentz model; standard, but the linewidth-only form is the one shown to fail.
  • ad hoc to paper NEP trained on DFT energies/forces, combined with fixed Born effective charges for dipole moments, captures electronic polarization dynamics without explicit charges in the potential
    Assumption underlying the MLP-based calculations; the paper notes remaining discrepancies for LiH below TO that may indicate long-range or nonlinear dipole effects.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Consistency between the Green-Kubo formula and Lorentz model for predicting the infrared dielectric function of polar materials." pith.science (2026). https://pith.science/paper/NGZRHLB3

@misc{pith2026250417464,
  author       = {Pith},
  title        = {Pith review of: Consistency between the Green-Kubo formula and Lorentz model for predicting the infrared dielectric function of polar materials},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/NGZRHLB3}},
  note         = {Machine review of arXiv:2504.17464}
}
read the original abstract

Accurate prediction of infrared dielectric functions in polar materials is fundamental for thermal and photonic applications, yet it remains unexplored whether the two main methods, Green-Kubo formula and Lorentz model, can give unified predictions. In this work, we present a detailed comparison of these two approaches using MgO and LiH as prototypical cases employing both empirical rigid ion model (RIM) and machine learning potential (MLP). We demonstrate that the conventional Lorentz model fails to capture the multi-phonon absorption inherent in Green-Kubo method, which can be resolved via using the phonon self-energy as a generalization of the usual linewidth. In addition, with RIM, a correction factor is required in the ionic contribution to infrared response to account for the electronic polarization effect, which is yet captured by MLP using the Born effective charges for calculating dipole moment. The present benchmark study thus enables cross-validation of dielectric function calculations while providing mechanistic insights into the polarization dynamics.

Figures

Figures reproduced from arXiv: 2504.17464 by the authors.

Figure 1
Figure 1. FIG. 1. The partial plots of total energy, atomic forces, virial for MgO (a-c) and LiH (d-f). The insets show the RMSE of the [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. The real and imaginary part of infrared dielectric function and semi-infinite normal reflectance of MgO (a)-(c) and [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. The real and imaginary part of infrared dielectric function and semi-infinite normal reflectance of MgO and LiH [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figures from the paper (8 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Infrared semi-infinite normal reflectance of MgO [PITH_FULL_IMAGE:figures/full_fig_p009_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5. The temperature-dependent infrared dielectric func [PITH_FULL_IMAGE:figures/full_fig_p009_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. The infrared dielectric function of LiH calculated [PITH_FULL_IMAGE:figures/full_fig_p010_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7. Phonon dispersion relations of MgO from: (a) [PITH_FULL_IMAGE:figures/full_fig_p010_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8. Phonon dispersion relations of LiH from: (a) [PITH_FULL_IMAGE:figures/full_fig_p010_8.png]
Figure 9
Figure 9. Figure 9: FIG. 9. Reflectance of (a) MgO and (b) LiH at 300 K calcu [PITH_FULL_IMAGE:figures/full_fig_p011_9.png]
Figure 11
Figure 11. Figure 11: FIG. 11. Real (a) and imaginary part (b) of phonon self [PITH_FULL_IMAGE:figures/full_fig_p011_11.png]
Figure 12
Figure 12. Figure 12: FIG. 12. The infrared dielectric function of MgO at 300 K [PITH_FULL_IMAGE:figures/full_fig_p012_12.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

94 extracted references · 70 canonical work pages

  1. [1]

    Thus, Eq.(25) will be: εion (ω) = 1 + [εion (0)− 1]ω2 0 ω2 0−ω2

    In the other limit of ω→∞ , εion (∞) = 1. Thus, Eq.(25) will be: εion (ω) = 1 + [εion (0)− 1]ω2 0 ω2 0−ω2 . (26) Considering the damping due to phonon scattering and the Lyddane–Sachs–Teller (LST) relationship, we obtain the Lorentz formula without electronic degree of freedom: εion (ω) = 1 + ω2 LO−ω2 TO ω2 TO−ω2−iωτ−1. (27) Compared Eq.(27) with Eq.(5), ...

  2. [2]

    The other atomic parameters for MgO are given in Table I [61]

    Rigid ion empirical potential We use a pairwise van Beest, Kramer, and van San- ten (BKS) potential form ϕij between atoms i and j composed of a short-range Buckingham term and a long- range Coulomb interaction which can be written as [60] ϕij = qiqj rij +A exp −rij ρ − C r6 ij , (29) where rij is the distance between atom i and atom j, qi and qj are the ...

    work page arXiv

  3. [3]

    A total of 1000 structures were sampled for MgO and LiH, with different lattice expansion respect to ground state (-2% ∼ +5%) and random displacement with magnitude of 0.1 ˚A

    Machine learning potential We use the NEP approach [55, 66] to construct ac- curate MLP models of MgO and LiH. A total of 1000 structures were sampled for MgO and LiH, with different lattice expansion respect to ground state (-2% ∼ +5%) and random displacement with magnitude of 0.1 ˚A. Each structure contains 128 atoms. The Quantum Espresso package [67] w...

  4. [4]

    J. R. Howell, M. P. Meng¨ u¸ c, and K. J. Daun,Thermal Radiation: An Introduction (CRC press, 2023)

    2023

  5. [5]

    Zernik, Theory of infrared optical properties of fluo- rite crystals, Rev

    W. Zernik, Theory of infrared optical properties of fluo- rite crystals, Rev. Mod. Phys. 39, 432 (1967)

    1967

  6. [6]

    J. R. Howell, M. P. Meng¨ u¸ c, K. J. Daun, and R. Siegel, Thermal radiation heat transfer (CRC press, 2020)

    2020

  7. [7]

    M. F. Modest and S. Mazumder, Radiative heat transfer (Academic press, 2021)

    2021

  8. [8]

    Overall, the phonon spectra derived from NEP, DFT 0 10 20 30 40 50 Frequency (THz) −40 −20 0 20 40 Re ε (a) cutoff 5 ˚ A cutoff 10 ˚ A 0 10 20 30 40 50 Frequency (THz) 0 20 40 60 80 100 Im ε (b) FIG. 6. The infrared dielectric function of LiH calculated from NEP with cutoff radius of 5 ˚A (blue solid lines) and 10 ˚A (red dashed lines): (a) real part; (b) i...

Show all 94 references
  1. [9]

    Zhao, Thermal Radiation: From Macro to Nano (Cambridge University Press, 2024)

    C. Zhao, Thermal Radiation: From Macro to Nano (Cambridge University Press, 2024)

    2024

  2. [10]

    Y. Wu, J. Duan, W. Ma, Q. Ou, P. Li, P. Alonso- Gonz´ alez, J. D. Caldwell, and Q. Bao, Manipulating po- laritons at the extreme scale in van der waals materials, Nat. Rev. Phys. 4, 578 (2022)

    2022

  3. [11]

    Foteinopoulou, G

    S. Foteinopoulou, G. C. R. Devarapu, G. S. Subrama- nia, S. Krishna, and D. Wasserman, Phonon-polaritonics: enabling powerful capabilities for infrared photonics, Nanophotonics 8, 2129 (2019)

    2019

  4. [12]

    Z. Tong, X. Yang, T. Feng, H. Bao, and X. Ruan, First- principles predictions of temperature-dependent infrared dielectric function of polar materials by including four- phonon scattering and phonon frequency shift, Phys. Rev. B 101, 125416 (2020)

    2020

  5. [13]

    C.-L. Zhou, Z. Torbatian, S.-H. Yang, Y. Zhang, H.-L. Yi, M. Antezza, D. Novko, and C.-W. Qiu, Unconventional thermophotonic charge density wave, Phys. Rev. Lett. 133, 066902 (2024)

    2024

  6. [14]

    E. D. Palik, Handbook of optical constants of solids , Vol. 3 (Academic press, San Diego, 1998)

    1998

  7. [15]

    Born and K

    M. Born and K. Huang, Dynamical theory of crystal lat- tices (Oxford university press, 1996)

    1996

  8. [16]

    Z. Han, C. Lee, J. Song, H. Wang, P. Bermel, and X. Ruan, Temperature-dependent full spectrum dielec- tric function of semiconductors from first principles, Phys. Rev. B 107, L201202 (2023)

    2023

  9. [17]

    H. Zhou, J. Tiwari, and T. Feng, Understanding the flat thermal conductivity of La 2Zr2O7 at ultrahigh tempera- tures, Phys. Rev. Mater. 8, 043804 (2024)

    2024

  10. [18]

    Tiwari and T

    J. Tiwari and T. Feng, Accurate prediction of thermal conductivity of Al 2O3 at ultrahigh temperatures, Phys. Rev. B 109, 075201 (2024)

    2024

  11. [19]

    W.-Z. Yuan, Y. Guo, C.-L. Zhou, and H.-L. Yi, Radia- tive thermal diode driven by the synergy between isotope engineering and lattice thermal expansion, Phys. Rev. B 110, 165412 (2024)

    2024

  12. [20]

    Z. Guo, P. Sokalski, Z. Han, Y. Cheng, L. Shi, T. Taniguchi, K. Watanabe, and X. Ruan, First- principles prediction of zone-center optical phonon linewidths and IR spectra of hexagonal boron nitride, Appl. Phys. Lett. 125, 231106 (2024)

    2024

  13. [21]

    Simoncelli, N

    M. Simoncelli, N. Marzari, and F. Mauri, Wigner formu- lation of thermal transport in solids, Phys. Rev. X 12, 041011 (2022)

    2022

  14. [22]

    H. Bao, B. Qiu, Y. Zhang, and X. Ruan, A first- principles molecular dynamics approach for predicting optical phonon lifetimes and far-infrared reflectance of polar materials, J. Quant. Spectrosc. Radiat. Transfer 113, 1683 (2012)

    2012

  15. [23]

    A. A. Maradudin and R. F. Wallis, Lattice anharmonic- ity and optical absorption in polar crystals. ii. classical treatment in the linear approximation, Phys. Rev. 123, 777 (1961). 13

    1961

  16. [24]

    Gangemi, A

    F. Gangemi, A. Carati, L. Galgani, R. Gangemi, and A. Maiocchi, Agreement of classical Kubo theory with the infrared dispersion curves n(ω) of ionic crystals, Eu- rophys. Lett. 110, 47003 (2015)

    2015

  17. [25]

    Domingues, A

    G. Domingues, A. M. Monthe, S. Gu´ evelou, and B. Rousseau, Study by molecular dynamics of the influ- ence of temperature and pressure on the optical proper- ties of undoped 3C-SiC structures, J. Quant. Spectrosc. Radiat. Transfer 205, 220 (2018)

    2018

  18. [26]

    Chen and L.-S

    W. Chen and L.-S. Li, The study of the optical phonon frequency of 3C-SiC by molecular dynamics simulations with deep neural network potential, J. Appl. Phys. 129, 244104 (2021)

    2021

  19. [27]

    Y. Guo, M. G´ omez Viloria, R. Messina, P. Ben-Abdallah, and S. Merabia, Atomistic modeling of extreme near-field heat transport across nanogaps between two polar dielec- tric materials, Phys. Rev. B 108, 085434 (2023)

    2023

  20. [28]

    D. J. Adams and I. R. McDonald, Rigid-ion models of the interionic potential in the alkali halides, J. Phys. C: Solid State Phys. 7, 2761 (1974)

    1974

  21. [29]

    G. K. Horton and A. A. Maradudin, Dynamical proper- ties of solids , Vol. 1 (North-Holland Publishing Company, Amsterdam, 1974)

    1974

  22. [30]

    Szigeti, Polarisability and dielectric constant of ionic crystals, Trans

    B. Szigeti, Polarisability and dielectric constant of ionic crystals, Trans. Faraday Soc. 45, 155 (1949)

    1949

  23. [31]

    Szigeti, Compressibility and absorption frequency of ionic crystals, Proc

    B. Szigeti, Compressibility and absorption frequency of ionic crystals, Proc. R. Soc. (London) A 204, 51 (1950)

    1950

  24. [32]

    Shanker and S

    J. Shanker and S. Dixit, Dielectric constants and their pressure and temperature derivatives for ionic crystals, Phys. Stat. Sol. (a) 123, 17 (1991)

    1991

  25. [33]

    Sirdeshmukh, L

    D. Sirdeshmukh, L. Sirdeshmukh, and K. Subhadra, Fifty years of Szigeti’s dielectric theory—a review, Pramana 69, 491 (2007)

    2007

  26. [34]

    Carati and A

    A. Carati and A. M. Maiocchi, Replacement of the lorentz law for the shape of the spectral lines in the infrared region, J. Opt. Soc. Am. A 33, 1193 (2016)

    2016

  27. [35]

    Denham, G

    P. Denham, G. R. Field, P. L. R. Morse, G. R. Wilkinson, and W. C. Price, Optical and dielectric properties and lattice dynamics of some fluorite structure ionic crystals, Proc. R. Soc. (London) A 317, 55 (1970)

    1970

  28. [36]

    D. D. S. Meneses, J.-F. Brun, P. Echegut, and P. Simon, Contribution of semi-quantum dielectric function models to the analysis of infrared spectra, Appl. Spectrosc. 58, 969 (2004)

    2004

  29. [37]

    T. Sun, P. B. Allen, D. G. Stahnke, S. D. Jacobsen, and C. C. Homes, Infrared properties of ferropericlase Mg1−xFexO: Experiment and theory, Phys. Rev. B 77, 134303 (2008)

    2008

  30. [38]

    Fugallo, B

    G. Fugallo, B. Rousseau, and M. Lazzeri, Infrared re- flectance, transmittance, and emittance spectra of MgO from first principles, Phys. Rev. B 98, 184307 (2018)

    2018

  31. [39]

    Benshalom, G

    N. Benshalom, G. Reuveni, R. Korobko, O. Yaffe, and O. Hellman, Dielectric response of rock-salt crystals at fi- nite temperatures from first principles, Phys. Rev. Mater. 6, 033607 (2022)

    2022

  32. [40]

    Amano, T

    T. Amano, T. Yamazaki, R. Akashi, T. Tadano, and S. Tsuneyuki, Lattice dielectric properties of rutile TiO2: First-principles anharmonic self-consistent phonon study, Phys. Rev. B 107, 094305 (2023)

    2023

  33. [41]

    J. Y. Yang, W. J. Zhang, L. H. Liu, J. Qiu, K. Wang, and J. Y. Tan, Temperature-dependent infrared dielec- tric functions of MgO crystal: An ellipsometry and first- principles molecular dynamics study, J. Chem. Phys. 141, 104703 (2014)

    2014

  34. [42]

    J. A. Thomas, J. E. Turney, R. M. Iutzi, C. H. Amon, and A. J. H. McGaughey, Predicting phonon dispersion relations and lifetimes from the spectral energy density, Phys. Rev. B 81, 081411 (2010)

    2010

  35. [43]

    T. Feng, B. Qiu, and X. Ruan, Anharmonicity and neces- sity of phonon eigenvectors in the phonon normal mode analysis, J. Appl. Phys. 117, 195102 (2015)

    2015

  36. [44]

    A. A. Maradudin and A. E. Fein, Scattering of neutrons by an anharmonic crystal, Phys. Rev. 128, 2589 (1962)

    1962

  37. [45]

    A. J. C. Ladd, B. Moran, and W. G. Hoover, Lattice ther- mal conductivity: A comparison of molecular dynamics and anharmonic lattice dynamics, Phys. Rev. B 34, 5058 (1986)

    1986

  38. [46]

    Monacelli and N

    L. Monacelli and N. Marzari, Electrostatic interactions in atomistic and machine-learned potentials for polar mate- rials (2024), arXiv:2412.01642 [cond-mat.mtrl-sci]

    2024 arXiv

  39. [47]

    Gonze and C

    X. Gonze and C. Lee, Dynamical matrices, born effective charges, dielectric permittivity tensors, and interatomic force constants from density-functional perturbation the- ory, Phys. Rev. B 55, 10355 (1997)

    1997

  40. [48]

    Cowley, The lattice dynamics of an anharmonic crys- tal, Adv

    R. Cowley, The lattice dynamics of an anharmonic crys- tal, Adv. Phys. 12, 421 (1963)

    1963

  41. [49]

    G. K. Horton and A. A. Maradudin, Dynamical proper- ties of solids , Vol. 2 (North-Holland Publishing Company, 1975)

    1975

  42. [50]

    T. Sun, X. Shen, and P. B. Allen, Phonon quasiparticles and anharmonic perturbation theory tested by molecular dynamics on a model system, Phys. Rev. B 82, 224304 (2010)

    2010

  43. [51]

    Shih and T

    P. Shih and T. C. Berkelbach, Anharmonic lattice dy- namics from vibrational dynamical mean-field theory, Phys. Rev. B 106, 144307 (2022)

    2022

  44. [52]

    Jasrasaria and T

    D. Jasrasaria and T. C. Berkelbach, Nonperturbative simulation of anharmonic rattler dynamics in type-i clathrates with vibrational dynamical mean-field theory, Phys. Rev. B 110, 064312 (2024)

    2024

  45. [53]

    Lu and F

    Y. Lu and F. Zheng, Lattice thermal conductivity of ZrSe2 based on the anharmonic phonon approach and on- the-fly machine learning force fields, Phys. Rev. B 109, 014305 (2024)

    2024

  46. [54]

    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)

    1957

  47. [55]

    T. Feng, W. Yao, Z. Wang, J. Shi, C. Li, B. Cao, and X. Ruan, Spectral analysis of nonequilibrium molecu- lar dynamics: Spectral phonon temperature and local nonequilibrium in thin films and across interfaces, Phys. Rev. B 95, 195202 (2017)

    2017

  48. [56]

    I. P. Ipatova, A. A. Maradudin, and R. F. Wallis, Tem- perature dependence of the width of the fundamental lattice-vibration absorption peak in ionic crystals. ii. ap- proximate numerical results, Phys. Rev. 155, 882 (1967)

    1967

  49. [57]

    Huang and M

    K. Huang and M. Born, On the interaction between the radiation field and ionic crystals, P. Roy. Soc. A-Math. Phy. 208, 352 (1951)

    1951

  50. [58]

    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 mat...

    2022

  51. [59]

    Z. Fan, Y. Wang, P. Ying, K. Song, J. Wang, Y. Wang, Z. Zeng, K. Xu, E. Lindgren, J. M. Rahm, A. J. Gabourie, J. Liu, H. Dong, J. Wu, Y. Chen, Z. Zhong, J. Sun, P. Erhart, Y. Su, and T. Ala-Nissila, GPUMD: A package for constructing accurate machine-learned po- tentials and pe...

    2022

  52. [60]

    Eastwood, R

    J. Eastwood, R. Hockney, and D. Lawrence, P3M3DP—the three-dimensional periodic particle- particle/ particle-mesh program, Comp. Phys. Comm. 19, 215 (1980)

    1980

  53. [61]

    Puligheddu, Y

    M. Puligheddu, Y. Xia, M. Chan, and G. Galli, Com- putational prediction of lattice thermal conductivity: A comparison of molecular dynamics and boltzmann trans- port approaches, Phys. Rev. Mater. 3, 085401 (2019)

    2019

  54. [62]

    H. Zhou, S. Zhou, Z. Hua, K. Bawane, and T. Feng, Im- pact of classical statistics on thermal conductivity predic- tions of bas and diamond using machine learning molec- ular dynamics, Appl. Phys. Lett. 125, 172202 (2024)

    2024

  55. [63]

    P. Ying, W. Zhou, L. Svensson, E. Berger, E. Frans- son, F. Eriksson, K. Xu, T. Liang, J. Xu, B. Song, S. Chen, P. Erhart, and Z. Fan, Highly efficient path- integral molecular dynamics simulations with gpumd us- ing neuroevolution potentials: Case studies on thermal propertie...

    2025

  56. [64]

    B. W. H. van Beest, G. J. Kramer, and R. A. van Santen, Force fields for silicas and aluminophosphates based on ab initio calculations, Phys. Rev. Lett. 64, 1955 (1990)

    1990

  57. [65]

    Matsui, Molecular dynamics study of the structural and thermodynamic properties of MgO crystal with quantum correction, J

    M. Matsui, Molecular dynamics study of the structural and thermodynamic properties of MgO crystal with quantum correction, J. Chem. Phys. 91, 489 (1989)

    1989

  58. [66]

    J. D. Gale and A. L. Rohl, The general utility lattice program GULP, Mol. Simul. 29, 291 (2003)

    2003

  59. [67]

    J. D. Gale, Empirical potential derivation for ionic ma- terials, Phil. Mag. B 73, 3 (1996)

    1996

  60. [68]

    Brodsky and E

    M. Brodsky and E. Burstein, Infrared lattice vibrations of single crystal lithium hydride and some of its isotopic derivations, J. Phys. Chem. Solids 28, 1655 (1967)

    1967

  61. [69]

    Shanker, A

    J. Shanker, A. K. G. Lashkari, and V. P. Gupta, Effect of three-body interactions on the electronic polarizabilities and photoelastic behaviour of alkali halides, Phys. Stat. Sol. (b) 91, 263 (1979)

    1979

  62. [70]

    Z. Fan, Z. Zeng, C. Zhang, Y. Wang, K. Song, H. Dong, Y. Chen, and T. Ala-Nissila, Neuroevolution machine learning potentials: Combining high accuracy and low cost in atomistic simulations and application to heat transport, Phys. Rev. B 104, 104309 (2021)

    2021

  63. [71]

    Giannozzi, O

    P. Giannozzi, O. Baseggio, P. Bonf` a, D. Brunato, R. Car, I. Carnimeo, C. Cavazzoni, S. de Gironcoli, P. Delugas, F. Ferrari Ruffino, A. Ferretti, N. Marzari, I. Timrov, A. Urru, and S. Baroni, Quantum ESPRESSO toward the exascale, J. Chem. Phys. 152, 154105 (2020)

    2020

  64. [72]

    G. I. Csonka, J. P. Perdew, A. Ruzsinszky, P. H. T. Philipsen, S. Leb` egue, J. Paier, O. A. Vydrov, and J. G. ´Angy´ an, Assessing the performance of recent density functionals for bulk solids, Phys. Rev. B 79, 155107 (2009)

    2009

  65. [73]

    Z. Deng, C. Chen, X.-G. Li, and S. P. Ong, An elec- trostatic spectral neighbor analysis potential for lithium nitride, npj Comput. Mater. 5, 75 (2019)

    2019

  66. [74]

    https://github.com/Yuanwz00/nep_data

  67. [75]

    J. R. Jasperse, A. Kahan, J. N. Plendl, and S. S. Mitra, Temperature dependence of infrared dispersion in ionic crystals LiF and MgO, Phys. Rev. 146, 526 (1966)

    1966

  68. [76]

    Giura, L

    P. Giura, L. Paulatto, F. He, R. P. S. M. Lobo, A. Bosak, E. Calandrini, L. Paolasini, and D. Antonangeli, Multi- phonon anharmonicity of MgO, Phys. Rev. B 99, 220304 (2019)

    2019

  69. [77]

    Tangney and S

    P. Tangney and S. Scandolo, A many-body interatomic potential for ionic systems: Application to MgO, J. Chem. Phys. 119, 9673 (2003)

    2003

  70. [78]

    X. Wu, W. Zhou, H. Dong, P. Ying, Y. Wang, B. Song, Z. Fan, and S. Xiong, Correcting force error-induced underestimation of lattice thermal conductivity in ma- chine learning molecular dynamics, J. Chem. Phys. 161, 014103 (2024)

    2024

  71. [79]

    W. Zhou, N. Liang, X. Wu, S. Xiong, Z. Fan, and B. Song, Insight into the effect of force error on the thermal con- ductivity from machine-learned potentials, Mater. Today Phys. 50, 101638 (2025)

    2025

  72. [80]

    J. T. Gourley and W. A. Runciman, Multiphonon in- frared absorption spectra of MgO and CaO, J. Phys. C: Solid State Phys. 6, 583 (1973)

    1973

  73. [81]

    Zhang, H

    L. Zhang, H. Wang, M. C. Muniz, A. Z. Panagiotopou- los, R. Car, and W. E, A deep potential model with long-range electrostatic interactions, J. Chem. Phys.156, 124107 (2022)

    2022

  74. [82]

    Borik, Non-linear dipole effects in the infra-red ab- sorption of GaAs, Phys

    H. Borik, Non-linear dipole effects in the infra-red ab- sorption of GaAs, Phys. Stat. Sol. (b) 39, 145 (1970)

    1970

  75. [83]

    R. M. Pick, M. H. Cohen, and R. M. Martin, Microscopic theory of force constants in the adiabatic approximation, Phys. Rev. B 1, 910 (1970)

    1970

  76. [84]

    Salanne and P

    M. Salanne and P. A. Madden, Polarization effects in ionic solids and melts, Mol. Phys. 109, 2299 (2011)

    2011

  77. [85]

    Ghosez, J.-P

    P. Ghosez, J.-P. Michenaud, and X. Gonze, Dynamical atomic charges: The case of ABO 3 compounds, Phys. Rev. B 58, 6224 (1998)

    1998

  78. [86]

    Behler, Four generations of high-dimensional neural network potentials, Chem

    J. Behler, Four generations of high-dimensional neural network potentials, Chem. Rev. 121, 10037 (2021)

    2021

  79. [87]

    Schall, M

    M. Schall, M. Walther, and P. Uhd Jepsen, Fundamental and second-order phonon processes in CdTe and ZnTe, Phys. Rev. B 64, 094301 (2001)

    2001

  80. [88]

    Tamura, Isotope scattering of dispersive phonons in ge, Phys

    S.-i. Tamura, Isotope scattering of dispersive phonons in ge, Phys. Rev. B 27, 858 (1983)

    1983

  81. [89]

    Rossi, M

    M. Rossi, M. Ceriotti, and D. E. Manolopoulos, How to remove the spurious resonances from ring polymer molec- ular dynamics, J. Chem. Phys. 140, 234116 (2014)

    2014

  82. [90]

    T. E. Markland and M. Ceriotti, Nuclear quantum effects enter the mainstream, Nat. Rev. Chem. 2, 0109 (2018)

    2018

  83. [91]

    A. R. Oganov, M. J. Gillan, and G. D. Price, Ab ini- tio lattice dynamics and structural stability of MgO, J. Chem. Phys. 118, 10174 (2003)

    2003

  84. [92]

    G. Roma, C. M. Bertoni, and S. Baroni, The phonon spectra of LiH and LiD from density-functional pertur- bation theory, Solid State Commun. 98, 203 (1996)

    1996

  85. [93]

    Posternak, A

    M. Posternak, A. Baldereschi, H. Krakauer, and R. Resta, Non-nominal value of the dynamical effective charge in alkaline-earth oxides, Phys. Rev. B 55, R15983 (1997)

    1997

  86. [94]

    L. L. Boyer, J. A. Harrington, M. Hass, and H. B. Rosen- stock, Multiphonon absorption in ionic crystals, Phys. Rev. B 11, 1665 (1975)

    1975

Pith tools

Reviewed August 16, 2026 · model on record in the stance chip above.