REVIEW 3 major objections 5 minor 42 references
Phonon angular momentum induced by Terahertz electric field
T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read A terahertz electric field can create a steady phonon angular momentum of order $\hbar$ per unit cell in time-reversal-symmetric polar crystals, through off-diagonal elements of the phonon angular momentum operator.
desk verdict A clean mechanism for THz-driven phonon AM in TRS crystals, but the headline magnitude needs an audit trail. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The mechanism is the second-order nonlinear response coefficient $\Pi^\alpha_{\nu\lambda}(\omega,-\omega)$ of Eq. (6), written as sums over distinct optical modes $\sigma\ne\sigma'$ of products of the field-coupling constants $B^\nu_\sigma$, the off-diagonal phonon angular momentum matrix elements $J^\alpha_{\sigma\sigma'}$, and resonance denominators containing the phonon lifetime $\tau$. Because $J^\alpha_{\sigma\sigma'}$ is pure imaginary under time-reversal symmetry, the sharp resonance sits in the imaginary part of $\Pi$, and for a field $\mathbf E(t)=\hat{\mathbf x}\cos(\omega_\Gamma t)+\hat{\mathbf y}\cos(\omega_\Gamma t+\varphi)$ the dc angular momentum reduces to $2\pi^2\,\mathrm{Im}\,\Pi^z_{xy}(\omega_\Gamma,-\omega_\Gamma)\sin\varphi$. This object converts an elliptically polarized THz field into a steady atomic rotation, with only $\Gamma$-point infrared-active optical modes participating because of the dipole coupling in Eq. (4).
What would settle it
Measure the orbital magnetization of wurtzite GaN driven by circularly polarized THz pulses at the 16.6 THz TO resonance with $E_0=10^5$ V/m: the paper predicts $M_z$ about $8.8\times10^{-3}\,\mu_B$/nm$^3$ with sign set by helicity and a $\sin\varphi$ phase dependence. Absence of that helicity-odd magnetization at this scale, or a resonance at a different frequency, would falsify the central claim.
Extended reading notes
Core claim
The paper's central claim is that a dc phonon angular momentum $J^\alpha_{\rm ph}(0)=\sum_{\nu\lambda}\int \Pi^\alpha_{\nu\lambda}(\omega,-\omega)E_\nu(\omega)E_\lambda(-\omega)\,d\omega$ emerges in time-reversal-symmetric polar crystals under alternating THz driving, with the response coefficient $\Pi^\alpha_{\nu\lambda}$ built from off-diagonal phonon angular momentum matrix elements $J^\alpha_{\sigma\sigma'}$ between distinct $\Gamma$-point optical modes. In such crystals the diagonal elements vanish by symmetry, so the effect comes entirely from one-phonon interband transitions, and it peaks when the driving frequency matches the optical phonon energy. The resonance appears in the imaginary part of $\Pi$, and the induced angular momentum varies as $\sin\varphi$ with the phase difference between orthogonal field components, reversing with the handedness of circularly polarized light. For wurtzite GaN, only the infrared-active $E_1$(TO) modes resonate, and at $|E_0|=10^5$ V/m the maximum phonon angular momentum is on the order of $\hbar$ per unit cell, corresponding to an orbital magnetization of about $8.8\times10^{-3}\,\mu_B$/nm$^3$.
Load-bearing premise
The prediction relies on the assumption that the THz-driven phonon system is well described by the free-phonon Hamiltonian plus a single constant lifetime $\tau$, with the field coupling only to $\Gamma$-point infrared modes; if anharmonicity, phonon-phonon interactions, or thermal occupation factors renormalize the modes, the resonant magnitude of order $\hbar$ per unit cell will change.
Editorial extensions
If this is right
- In any of the 32 point-group symmetric crystals, time-reversal symmetry no longer blocks macroscopic phonon angular momentum: a THz field can create a nonzero dc value through interband transitions.
- Resonant circularly polarized THz driving of wurtzite GaN at about 16.6 THz yields a phonon angular momentum of order $\hbar$ per unit cell at $10^5$ V/m, detectable as an orbital magnetization near $10^{-2}\,\mu_B$/nm$^3$.
- The generated angular momentum is proportional to $\sin\varphi$, so left- and right-circular polarization produce opposite angular momentum and the magnitude can be tuned continuously by the phase difference.
- Since the effect scales quadratically with field strength, stronger THz sources should push the magnetization well above noise floors without changing the mechanism.
Reading between the lines
- This implies the same formula extends to any polar crystal with infrared-active modes; the main material requirement is a pair of near-degenerate optical modes with nonzero off-diagonal phonon angular momentum, so the effect may be tunable by alloying or strain.
- Because the response is quadratic in the field and controlled by phase, this offers a THz-rate, helicity-sensitive knob for phononic angular momentum, akin to optical orientation but for lattice rotation.
- The constant-lifetime approximation likely underestimates temperature dependence: including anharmonic phonon lifetimes would make the resonance amplitude and line shape temperature dependent in a way that could be checked by varying temperature in the same experiment.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a quantum nonlinear response theory for phonon angular momentum (AM) induced by a terahertz electric field in polar insulators. Starting from a free-phonon Hamiltonian and a dipole-approximation coupling to the field, the authors write the dc phonon AM as J^α_ph(0) = ∫ Π^α_νλ(ω,-ω) E_ν(ω) E_λ(-ω) dω, with the response coefficient Π given by Eq. (6) in terms of field-phonon couplings B and off-diagonal phonon AM matrix elements J^α_σσ'. In time-reversal-symmetric systems the diagonal AM elements vanish, so the effect is attributed to interband phonon transitions. The paper predicts a resonant feature in the imaginary part of Π and a sin(φ) phase dependence of the induced AM, and it demonstrates these in a two-dimensional hexagonal toy model and in first-principles wurtzite GaN, where it claims a maximum AM of order ℏ per unit cell at E0 = 10^5 V/m, measurable through the induced orbital magnetization.
Significance. If the quantitative claims are correct, this would be an important advance: it offers a general mechanism for generating phonon AM in ordinary polar crystals with time-reversal symmetry, thereby addressing a long-recognized symmetry bottleneck. The explicit formula Eq. (6), the predicted imaginary-part resonance, and the phase-difference control are crisp, testable predictions. The use of first-principles phonon eigenvectors and Born effective charges for GaN, with a clear contrast to chiral-phonon selection-rule scenarios, is a strength. However, the headline magnitude for GaN is not currently reproducible from the text, and the central derivation is deferred to missing supplemental material, so the significance is conditional on those points being resolved.
major comments (3)
- [Phonon model in dipole approximation, Eq. (6)] The central formula, Eq. (6), is stated without derivation in the main text; the derivation is relegated to 'Supplemental Material I and II', and the posted manuscript contains no supplemental material. The symmetry analysis and the first-principles details are likewise deferred to Supplemental Materials III, IV, and V. Because Eq. (6) is the basis for all subsequent predictions, the manuscript as posted is not self-contained and cannot be independently validated. I request that the derivation, or at least a complete outline with the essential intermediate steps and the domain of validity of the constant-relaxation-time approximation, be included in the paper or in a posted supplement.
- [Wurtzite GaN, Fig. 3] The headline quantitative claim—AM on the order of ℏ per unit cell at |E0|=10^5 V/m—is not supported by the information provided and appears inconsistent with the paper's own formulas by more than an order of magnitude. Reducing Eq. (6) to the resonant E1(TO) doublet (modes 7 and 8) and combining with Eq. (9) at resonance, with the stated parameters Z*=2.83e, M = m_Ga m_N/(m_Ga+m_N), ω=16.6 THz, and τ=33.95 ps, gives a per-cell AM J^z ≈ 4π (Z*e)^2 (ℏ/(2Mω)) τ² E0² / ℏ ≈ 0.07 ℏ, roughly fourteen times smaller than the claimed order ℏ. The manuscript does not display the actual numerical evaluation of Eq. (6), the phonon eigenvector normalization, the computed matrix elements J^z_78, or the field convention used in Fig. 3(b). Please provide the complete evaluation and reconcile the discrepancy; if the full multiband sum or a different eigenvector basis changes the result, that enhancement must be shown explicitly.
- [Introduction / Discussion] The statement that the phenomenon is 'applicable to all 32-point group symmetric systems' is too broad as presented. According to Eq. (6), the magnitude is governed by the denominator (ω_{σ'}-ω_σ - i/(2τ)), so the effect is largest when two IR-active modes are degenerate or nearly degenerate, as in the E1 doublet of GaN. In point groups where all IR modes are non-degenerate, the same tensor may be symmetry-allowed but the response is suppressed by the mode-splitting denominator. The paper should either present the missing symmetry analysis from Supplemental Material III or qualify the generality claim so that it accurately reflects this sensitivity to mode degeneracy.
minor comments (5)
- [Eq. (6)] The symbol ω′_σ in the second denominator is undefined; it should presumably be ω_σ or ω_{σ′}. Please correct this and define all frequencies explicitly.
- [Eq. (7)] The notation 'ω−ωΓ+ 1/(4τ²)/(ω−ωΓ)' is ambiguous; add parentheses to clarify that this means (ω−ωΓ) + (1/(4τ²))/(ω−ωΓ), and in general use consistent bracketing in the displayed equations.
- [Fig. 2(d-f)] The vertical axes of the nonlinear coefficient plots are not labeled with units; please specify the units of Π and the normalization (e.g., per primitive cell or per unit volume) used in the calculation.
- [Abstract / Eq. (9)] The phrase 'nonlinear polarized light' is unclear; describe the field as two orthogonal linearly polarized components with a controllable phase difference, since the nonlinearity refers to the response order, not to the light itself.
- [Wurtzite GaN / Eq. (6)] The temperature at which τ=33.95 ps applies is not stated, although the Discussion emphasizes the temperature dependence; please give the temperature or otherwise specify the lifetime's range of validity.
Circularity Check
No circular derivation: the response formula is a standard second-order perturbation result; the only self-citations are structural operator inputs, not fitted predictions.
full rationale
The derivation chain is: Hamiltonian (1)-(4) -> nonlinear response formula (5)-(6) -> off-diagonal phonon AM matrix elements -> 2D hexagonal model and wurtzite GaN evaluation. Equation (6) is not defined in terms of the predicted phonon AM; it combines independently available inputs (Born effective charges, phonon eigenvectors and frequencies, and a literature relaxation time). The reported GaN magnitude is a forward calculation from these inputs, not a parameter fitted to the target observable. The paper cites Ref. [1], which shares an author, for the second-quantized form of the phonon AM operator; this is a structural operator definition whose stated assumptions (harmonic phonons, TRS) do not include the target nonlinear response result, so the citation is not load-bearing circularity. Ref. [16], also by overlapping authors, is cited for the physical importance of off-diagonal phonon AM contributions, but the derivation does not reduce to that prior work. The phase dependence in Eq. (9) and the imaginary-part resonance are explicit consequences of the derived coefficient and the symmetry properties of J, not imported conclusions. The quantitative claim of roughly hbar per unit cell at 1e5 V/m may warrant a reproducibility check (a simple two-mode evaluation of Eq. (6) with the stated parameters suggests a smaller magnitude), but this is a numerical/correctness concern rather than evidence that the prediction is equivalent to an input. No circular step is identifiable from the manuscript text.
Assumptions & free parameters
free parameters (3)
- Phenomenological phonon lifetime tau =
33.95 ps for GaN (Ref. [35])
- Driving field amplitude E0 =
10^5 V/m
- Off-diagonal phonon AM matrix elements J^alpha_sigma sigma' =
Computed from phonon eigenvectors and masses
assumptions (4)
- domain assumption Kubo-style nonlinear response theory for bosonic phonons with a classical electric field is valid for computing the dc phonon AM (Supplemental Materials I-II).
- domain assumption Dipole approximation: the spatially uniform THz field couples only to q=0 optical phonons through Born effective charges (Eqs. (3)-(4)).
- domain assumption Constant relaxation time approximation with a single phenomenological lifetime tau for all phonon modes.
- domain assumption In time-reversal-symmetric systems, diagonal matrix elements of the phonon AM operator vanish and off-diagonal elements are pure imaginary, so only interband transitions contribute.
Cite this review
Pith. "Pith review of Phonon angular momentum induced by Terahertz electric field." pith.science (2026). https://pith.science/paper/EE3ZGXTL
@misc{pith2026250605715,
author = {Pith},
title = {Pith review of: Phonon angular momentum induced by Terahertz electric field},
year = {2026},
howpublished = {\url{https://pith.science/paper/EE3ZGXTL}},
note = {Machine review of arXiv:2506.05715}
}
read the original abstract
Despite the growing interest in phonon angular momentum (AM) in recent years, current studies remain limited to a few materials due to the constraints imposed by time reversal symmetry on macroscopic phonon AM. In this work, we theoretically investigate the generation of total phonon AM through alternating terahertz electric fields in polarized materials. In contrast to previous studies on phonon AM, here the off-diagonal elements of the phonon AM operator play an essential role. According to our formula, the large AM is generated when the energy of incident electric fields matches the frequency of optical phonons at {\Gamma} point. Furthermore, a specific resonance on the imaginary part of the response coefficient, as well as periodic regulation of the phonon AM by the phase difference of the driving field, is observed. In polar material GaN, the oscillation maximum is observed as \hbar per unit cell which can be experimentally measured through orbital magnetization induced by phonon AM. Our work offers a promising approach to generate observable phonon AM in a wider range of materials, advancing both the understanding of phonon fundamental physics and potential applications in phononic devices.
Figures
Reference graph
Works this paper leans on
- [1]
- [2]
- [3]
-
[4]
D. M. Juraschek and N. A. Spaldin, Phys. Rev. Mater3, 064405 (2019)
work page 2019
- [5]
-
[6]
J. J. Nakane and H. Kohno, Phys. Rev. B 97, 174403 (2018)
work page 2018
-
[7]
Z. Li, T. Wang, C. Jin, Z. Lu, Z. Lian, Y. Meng, M. Blei, M. Gao, T. Taniguchi, K. Watanabe,et al., ACS nano 13, 14107 (2019)
work page 2019
-
[8]
G. Grissonnanche, S. Th ´eriault, A. Gourgout, M.-E. Boulanger, E. Lefranc ¸ois, A. Ataei, F. Lalibert ´e, M. Dion, J.-S. Zhou, S. Pyon, et al., Nat. Phys. 16, 1108 (2020)
work page 2020
Show all 42 references
-
[9]
S. R. Tauchert, M. Volkov, D. Ehberger, D. Kazenwadel, M. Evers, H. Lange, A. Donges, A. Book, W. Kreuzpaintner, U. Nowak, et al., Nature 602, 73 (2022)
2022
-
[10]
Vittmann, J
C. Vittmann, J. Lim, D. Tamascelli, S. F. Huelga, and M. B. Plenio, J. Phys. Chem. Lett. 14, 340 (2023)
2023
-
[11]
K. Kim, E. Vetter, L. Yan, C. Yang, Z. Wang, R. Sun, Y. Yang, A. H. Comstock, X. Li, J. Zhou, et al. , Nat. Mater. 22, 322 6 (2023)
2023
-
[12]
Y. Gao, Y. Pan, J. Zhou, and L. Zhang, Phys. Rev. B108, 064510 (2023)
2023
-
[13]
Komiyama and S
H. Komiyama and S. Murakami, Phys. Rev. B 103, 214302 (2021)
2021
-
[14]
K. C. Lee, M. R. Sprague, B. J. Sussman, J. Nunn, N. K. Lang- ford, X.-M. Jin, T. Champion, P. Michelberger, K. F. Reim, D. England, et al., Science 334, 1253 (2011)
2011
-
[15]
Hamada, E
M. Hamada, E. Minamitani, M. Hirayama, and S. Murakami, Phys. Rev. Lett121, 175301 (2018)
2018
-
[16]
Zhong, H
J. Zhong, H. Sun, Y. Pan, Z. Wang, X. Xu, L. Zhang, and J. Zhou, Phys. Rev. B107, 125147 (2023)
2023
-
[17]
Rostami, F
H. Rostami, F. Guinea, and E. Cappelluti, Phys. Rev. B 105, 195431 (2022)
2022
-
[18]
Pan and F
Y. Pan and F. Caruso, NPJ 2D Mater.Appl. 8, 42 (2024)
2024
-
[19]
Basini, M
M. Basini, M. Pancaldi, B. Wehinger, M. Udina, V. Unikanda- nunni, T. Tadano, M. Hoffmann, A. Balatsky, and S. Bonetti, Nature , 1 (2024)
2024
-
[20]
J. Luo, T. Lin, J. Zhang, X. Chen, E. R. Blackert, R. Xu, B. I. Yakobson, and H. Zhu, Science382, 698 (2023)
2023
-
[21]
Kuznetsov and C
A. Kuznetsov and C. Stanton, Phys. Rev. B 51, 7555 (1995)
1995
-
[22]
Fu and M
Z. Fu and M. Yamaguchi, Sci. Rep. 6, 38264 (2016)
2016
-
[23]
Ishizuka and M
H. Ishizuka and M. Sato, Phys. Rev. Lett 129, 107201 (2022)
2022
-
[24]
Ishizuka and M
H. Ishizuka and M. Sato, Phys. Rev. B 110, L020303 (2024)
2024
-
[25]
D. M. Juraschek, M. Fechner, A. V. Balatsky, and N. A. Spaldin, Phys. Rev. Mater.1, 014401 (2017)
2017
-
[26]
R. M. Geilhufe, V. Juri ˇci´c, S. Bonetti, J.-X. Zhu, and A. V. Balatsky, Phys. Rev. Res.3, L022011 (2021)
2021
-
[27]
R. M. Geilhufe, Phys. Rev. Res. 4, L012004 (2022)
2022
-
[28]
R. M. Geilhufe and W. Hergert, Phys. Rev. B 107, L020406 (2023)
2023
-
[29]
Gonze and C
X. Gonze and C. Lee, Phys. Rev. B 55, 10355 (1997)
1997
-
[30]
Kraut and R
W. Kraut and R. von Baltz, Phys. Rev. B 19, 1548 (1979)
1979
-
[31]
Huang, Y.-H
Y.-S. Huang, Y.-H. Chan, and G.-Y. Guo, Phys. Rev. B 108, 075413 (2023)
2023
-
[32]
Ishito, H
K. Ishito, H. Mao, Y. Kousaka, Y. Togawa, S. Iwasaki, T. Zhang, S. Murakami, J.-i. Kishine, and T. Satoh, Nat. Phys. 19, 35 (2023)
2023
-
[33]
Ishito, H
K. Ishito, H. Mao, K. Kobayashi, Y. Kousaka, Y. Togawa, H. Kusunose, J.-i. Kishine, and T. Satoh, Chirality 35, 338 (2023)
2023
-
[34]
W. J. Choi, S. H. Lee, B. C. Park, and N. A. Kotov, J. Am. Chem. Soc. 144, 22789 (2022)
2022
-
[35]
Choi and S.-H
Y.-J. Choi and S.-H. Jhi, Phys. Rev. B 106, 094311 (2022)
2022
-
[36]
Togo and I
A. Togo and I. Tanaka, SCRIPTA MATER 108, 1 (2015)
2015
-
[37]
Kresse and J
G. Kresse and J. Furthm¨ uller, Phys. Rev. B54, 11169 (1996)
1996
-
[38]
Kroumova, M
E. Kroumova, M. Aroyo, J. Perez-Mato, A. Kirov, C. Capillas, S. Ivantchev, and H. Wondratschek, Phase Transitions 76, 155 (2003)
2003
-
[39]
Carlos, J
W. Carlos, J. Freitas Jr, M. A. Khan, D. Olson, and J. Kuznia, Phys. Rev. B48, 17878 (1993)
1993
-
[40]
C. Chen, D. Zhai, C. Xiao, and W. Yao, Phys. Rev. Res. 6, L012059 (2024)
2024
-
[41]
Z. Zhou, C. Min, H. Ma, Y. Zhang, X. Xie, H. Zhan, and X. Yuan, Opt. Express 30, 13416 (2022)
2022
-
[42]
X. Feng, J. Cao, Z.-F. Zhang, L. K. Ang, S. Lai, H. Jiang, C. Xiao, and S. A. Yang, arXiv preprint arXiv:2409.09669 (2024)
2024 arXiv
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.