Pith. sign in

REVIEW 2 major objections 4 minor 37 references

Coherent phonon modulations in TR-ARPES yield state-resolved electron-phonon couplings and unify two generation mechanisms.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · grok-4.5

2026-07-15 03:22 UTC pith:MFTLPMPR

load-bearing objection Solid first-principles TD-aGW + classical e-ph simulation that cleanly extracts state-resolved g from TR-ARPES modulations and unifies ISRS/DECP under one driving term; classical/q=0 limits are stated, not hidden. the 2 major comments →

arxiv 2607.12783 v1 pith:MFTLPMPR submitted 2026-07-14 cond-mat.mtrl-sci

Photogeneration and signatures of coherent phonons in time-resolved photoemission spectroscopy: First-principles time-dependent adiabatic GW approach

classification cond-mat.mtrl-sci
keywords coherent phononsTR-ARPESelectron-phonon couplingtime-dependent GWmonolayer MoS2ISRSDECPRaman selection rule
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

This paper shows that a first-principles time-dependent adiabatic GW calculation that includes electron-phonon coupling can both generate coherent phonons with a light pulse and predict how those phonons shift and modulate the TR-ARPES spectrum of monolayer MoS2. The amplitude of the intensity oscillation at each band and crystal momentum is linearly proportional to the phonon displacement amplitude, and the slope of that line recovers the microscopic electron-phonon matrix element of that state. The same driving force that appears in the phonon equation of motion produces the classic signatures of impulsive stimulated Raman scattering below the optical gap and of displacive excitation near and above the gap; the two regimes are therefore continuous rather than mutually exclusive. Because the driving term shares the same optical matrix product that appears in the Raman tensor, only Raman-active modes are generated, yet the frequency dependence of the two spectra is not identical. The practical payoff is a route to extract state-resolved couplings from existing pump-probe ARPES data and to design coherent-phonon control of band structure.

Core claim

State-resolved electron-phonon coupling strengths can be read off from the linear relation between the amplitude of coherent-phonon-induced TR-ARPES intensity modulations and the phonon displacement amplitude; the same density-matrix response that drives the phonon also produces both the impulsive-stimulated-Raman and displacive-excitation signatures depending on pump frequency, with no sharp transition between them.

What carries the argument

The coupled equations of motion for the single-particle density matrix (in the adiabatic GW Hamiltonian that includes electron-electron, electron-hole and classical electron-phonon terms) and for the zone-center phonon displacement expectation value, whose driving force is the instantaneous expectation value of the electron-phonon matrix elements.

Load-bearing premise

Treating phonons classically, restricting them to zone-center modes, and omitting phonon-phonon and higher-order electron-phonon couplings is assumed enough to capture both generation mechanisms and the TR-ARPES signatures.

What would settle it

Measure the pump-frequency dependence of the phase and the time-averaged displacement of the A1' coherent phonon in monolayer MoS2 TR-ARPES; if the phase never shifts by approximately π/2 and the average displacement never becomes finite near the A-exciton energy, or if the extracted dAamp/dQamp fails to track first-principles g_nn for multiple valence bands, the central claim is wrong.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • TR-ARPES intensity oscillations at known phonon frequencies can be inverted to obtain band- and momentum-resolved deformation potentials without additional spectroscopy.
  • Coherent-phonon amplitudes can be used as a quantitative knob to shift selected bands or to open/close gaps on the femtosecond scale.
  • The continuous crossover between ISRS and DECP implies that intermediate pump frequencies produce mixed sine-cosine phases rather than a binary switch.
  • Only Raman-active modes appear, but their generation spectrum differs from the Raman intensity spectrum below the optical gap and above the C-exciton peak.

Where Pith is reading between the lines

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

  • The same linear-extraction protocol should work for other layered semiconductors once the probe energy resolution is sufficient to isolate non-degenerate bands.
  • Including dephasing rates of a few tens of meV will collapse the time-averaged displacement, making experimental phase analysis more ambiguous and favoring amplitude-based rather than phase-based mechanism assignment.
  • Finite-q or multi-phonon processes, if they dominate in thicker samples or different materials, would invalidate the direct mapping from dA/dQ to g_nn and would require an extended theory.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 4 minor

Summary. The manuscript develops a first-principles time-dependent adiabatic GW (TD-aGW) framework that couples the electronic density matrix to classical zone-center phonon displacements via electron-phonon matrix elements. Applied to monolayer MoS2, the method simultaneously generates coherent phonons under a short optical pump and computes the resulting time-resolved ARPES spectral function. The central claims are that (i) the amplitude of coherent-phonon-induced modulations in the TR-ARPES difference spectrum is linearly proportional to the phonon displacement amplitude, with the slope equal (up to a known factor) to the state-resolved e-ph coupling g_nn^ u(k,q=0), and (ii) both the impulsive stimulated Raman scattering (ISRS) and displacive excitation of coherent phonons (DECP) mechanisms emerge from the same driving term in the phonon equation of motion, with a continuous crossover as the pump frequency approaches the optical gap. Supporting evidence includes intensity scaling of amplitudes, a π/2 phase jump across the gap, and direct comparison of fitted dA_amp/dQ_amp against independently computed DFPT e-ph matrix elements for the A1' mode.

Significance. If the extraction procedure holds under realistic experimental conditions, the work supplies a concrete, state- and momentum-resolved route to e-ph couplings that optical pump-probe methods cannot access. The unification of ISRS and DECP within a single first-principles density-matrix response, together with an explicit comparison of the coherent-phonon amplitude spectrum against the Raman intensity, clarifies long-standing phenomenological distinctions and selection-rule coincidences. The calculations are parameter-light once the electronic structure and e-ph matrix elements are fixed, and the linear relations demonstrated in Fig. 2c-d constitute a falsifiable prediction that can be tested by combined TR-ARPES and structural probes. These features make the paper a useful methodological and conceptual advance for the ultrafast spectroscopy community.

major comments (2)
  1. The classical phonon approximation and restriction to q=0 modes (Eqs. 3-4 and surrounding text) are stated explicitly, yet they remain load-bearing for both the generation mechanism assignment and the extraction of g. The manuscript should quantify, at least for one representative pump frequency, the expected magnitude of finite-q scattering or quantum phonon fluctuations that would invalidate the linear dA_amp/dQ_amp relation, or else clearly demarcate the regime of validity in the abstract and conclusions.
  2. Dissipation is omitted from the main simulations (Sec. III and Fig. 1-4) and only briefly explored in the SM. Because dephasing/depopulation on the 100-fs scale converts the finite mean displacement Q_bar into a decaying oscillation about zero and renders phase extraction ambiguous, the claim that TR-ARPES can routinely distinguish ISRS from DECP needs a quantitative statement of the dephasing rates below which the distinction survives.
minor comments (4)
  1. Fig. 2(d) axis labels and legend use mixed notation (g%%! vs g_nn^ u); a single consistent symbol set would improve readability.
  2. The precise definition of the reference time t=237 fs used for ΔA is given only in the caption of Fig. 1; it should be restated when the fitting procedure is introduced in Sec. III.
  3. A short sentence clarifying that the off-diagonal force terms were numerically verified to cancel for the A1' mode (rather than assumed a priori) would strengthen the discussion around Fig. 3.
  4. Typographical inconsistencies appear in phonon-mode labels (A'1 vs A1') and in the units of E_max across figure panels; these should be standardized.

Circularity Check

1 steps flagged

Minor self-definitional consistency check that spectral-modulation slopes recover the input g; core generation-mechanism and selection-rule results follow independently from the EOMs.

specific steps
  1. self definitional [Sec. III, Fig. 2(c–d) and surrounding text]
    "In Fig. 2 (d) we compare the fitted slopes against the e-ph coupling matrix elements g!!$(k,q=0) of the corresponding state n for u=A1'. The perfect proportionality confirms that dAamp/dQamp is proportional to g!!$(k,q=0) and, in principle, the state-dependent e-ph coupling strength can be extracted from measuring riangle A(k,\omega,t) and Q."

    By construction H_e-ph(t)=g〈Q(t)〉 (Eq. 3) shifts the instantaneous quasiparticle energies by exactly g_nn〈Q〉. Consequently the peak positions of the spectral function A (and therefore the amplitude of the intensity-difference signal riangle A at fixed energy) must scale linearly with g Q_amp. Plotting the fitted slopes versus the input g is guaranteed to produce a straight line through the origin; the “confirmation” is tautological within the model and does not constitute an independent extraction.

full rationale

The paper’s load-bearing outputs—the time-dependent phonon displacements, the pump-frequency crossover between persistent-force (DECP-like) and vanishing-force (ISRS-like) regimes, the associated phase shifts, and the comparison of coherent-phonon amplitude versus Raman intensity—are obtained by direct numerical solution of the coupled density-matrix and classical-phonon equations (Eqs. 1–4) that take as input independently computed GW quasiparticle energies, optical matrix elements and DFPT e-ph matrix elements. No free parameters are adjusted to produce these features. The linear relation dA_amp/dQ_amp o g_nn shown in Fig. 2d is a direct algebraic consequence of the same linear e-ph term that is already present in the Hamiltonian (Eq. 3); it therefore functions only as an internal consistency check of the extraction protocol and does not circularly generate the mechanism or selection-rule claims. Self-citations supply the prior electronic TD-aGW propagator and are not invoked to justify the phonon results. The derivation chain is therefore essentially self-contained.

Axiom & Free-Parameter Ledger

3 free parameters · 5 axioms · 0 invented entities

The central claims rest on standard many-body approximations (GW, adiabatic COHSEX, classical phonons) plus the authors' prior TD-aGW propagator; free parameters are only the external pump settings and optional phenomenological dephasing rates used in SM. No new particles or forces are invented. The ledger therefore records domain assumptions of the electronic-structure community rather than ad-hoc constructs.

free parameters (3)
  • pump duration and envelope = 10 fs (main text)
    Fixed at 10 fs truncated half-period sine-squared (and 20 fs in SM); chosen to be short compared with phonon periods but not derived from first principles.
  • maximum field amplitude E_max = 10^{-3} a.u. (reference calculations)
    Scanned from 10^{-4} to 10^{-3} a.u. to demonstrate linearity; absolute scale is arbitrary within the linear-response regime.
  • phenomenological dephasing/depopulation rates = 0.01-0.05 eV
    Optional rates 0.01-0.05 eV introduced only in SM to illustrate dissipation; not used for main claims but affect phase extraction when present.
axioms (5)
  • domain assumption Static COHSEX approximation for the time-dependent GW self-energy (adiabatic limit)
    Invoked in Eq. (2) and the method section; standard for TD-aGW but neglects dynamical screening during the pulse.
  • domain assumption Classical treatment of phonon displacement operators (expectation-value dynamics only)
    Eqs. (3)-(4); quantum phonon fluctuations and zero-point motion are omitted by construction.
  • domain assumption Restriction to zone-center (q=0) phonon modes
    Stated after Eq. (3); justified by optical selection rules but excludes finite-q scattering channels.
  • domain assumption GW quasiparticle Hamiltonian plus Berry-connection treatment of the position operator
    Inherited from the authors' prior TD-aGW framework (Ref. 28); assumed accurate for MoS2 excitons.
  • ad hoc to paper Off-diagonal elements of the e-ph force contribute only secondary, rapidly oscillating terms that can be neglected for the A1' mode
    Numerical observation stated in the discussion of Fig. 3; not proven generally.

pith-pipeline@v1.1.0-grok45 · 16340 in / 3307 out tokens · 32762 ms · 2026-07-15T03:22:07.620429+00:00 · methodology

0 comments
read the original abstract

Coherent lattice dynamics can be observed in pump-probe time-resolved and angle-resolved photoemission spectroscopy (TR-ARPES) as a periodic modulation of intensity and energy of photoelectrons over probe time. We present an ab initio time-dependent GW approach including electron-phonon (e-ph) couplings to simulate the photogeneration of coherent phonons and their effects on the TR-ARPES of monolayer MoS2. We demonstrate that state-resolved e-ph coupling strength can be obtained from analyzing coherent phonon modulations on TR-ARPES. Features of both the impulsive stimulated Raman scattering mechanism and the displacive excitation mechanism of coherent phonon generation are identified in our simulations. We clarify their origins and the coincident selection rule of coherent phonon generation and Raman scattering intensity. This method provides supports to analyze coherent phonon dynamics and e-ph couplings in TR-ARPES and enables quantitative engineering of band structure through coherent phonons.

Figures

Figures reproduced from arXiv: 2607.12783 by Steven G. Louie, Yang-hao Chan, Zhenglu Li.

Figure 1
Figure 1. Figure 1: FIG. 1. (Color online) (a) Snapshot of computed [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. (Color online) (a) Computed temporal evolution of spectrum intensity difference from the spectrum at [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. (Color online) (a) The temporal evolution of [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. (Color online) The temporal dependence of (a) the [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗

discussion (0)

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

Reference graph

Works this paper leans on

37 extracted references · 3 canonical work pages

  1. [1]

    /𝑠𝑖𝑛T𝜔/0𝑡+𝜙78V+Δ𝐴WWWW, where Δ𝐴6

    give an excitation energy for the lowest optically bright ’A’ exciton of 1.76 eV and a band gap of 2.27 eV. Fig. 1 (a) shows the simulated TR-ARPES for monolayer MoS2 with a 1.7 eV pump light (near resonant with the A exciton excitation energy) at a probe time of 237 fs after the pump pulse. We define 𝑡=0 as the time when the pump starts as shown in Eq. S...

  2. [2]

    / and 𝑄6

    From these results we conclude that the coherent phonon generated with a 1.7 eV pump frequency is due to the e-ph coupling from the excited carriers, which is in line with the DECP mechanism [18]. Next, we investigate the pump frequency dependence of the TR-ARPES difference spectrum and coherent phonon. Fig. 4 (a) and (b) show the temporal evolution of Δ𝐴...

  3. [3]

    Heinrich, H.-T

    T. Heinrich, H.-T. Chang, S. Zayko, K. Rossnagel, M. Sivis, and C. Ropers, Phys. Rev. X 13, 021033 (2023)

  4. [4]

    M. Hase, M. Katsuragawa, A. M. Constantinescu, and H. Petek, New Journal of Physics 15, 055018 (2013)

  5. [5]

    J.-Y. Shan, M. Ye, H. Chu, S. Lee, J.-G. Park, L. Balents, and D. Hsieh, Nature 600, 235 (2021)

  6. [6]

    Schlipf and F

    M. Schlipf and F. Gygi, Computer Physics Communications 196, 36 (2015)

  7. [7]

    X. Li, T. Qiu, J. Zhang, E. Baldini, J. Lu, A. M. Rappe, and K. A. Nelson, Science 364, 1079 (2019), https://www.science.org/doi/pdf/10.1126/science.aaw4913

  8. [8]

    D. N. Basov, R. D. Averitt, and D. Hsieh, Nature Materials 16, 1077 (2017)

  9. [9]

    A. S. Disa, T. F. Nova, and A. Cavalleri, Nature Physics 17, 1087 (2021)

  10. [10]

    de la Torre, D

    A. de la Torre, D. M. Kennes, M. Claassen, S. Gerber, J. W. McIver, and M. A. Sentef, Rev. Mod. Phys. 93, 041002 (2021)

  11. [11]

    Baldini, T

    E. Baldini, T. Palmieri, A. Dominguez, P. Ruello, A. Rubio, and M. Chergui, Nano Lett. 18, 5007 (2018)

  12. [12]

    T. Y. Jeong, B. M. Jin, S. H. Rhim, L. Debbichi, J. Park, Y. D. Jang, H. R. Lee, D.-H. Chae, D. Lee, Y.-H. Kim, S. Jung, and K. J. Yee, ACS Nano 10, 5560 (2016)

  13. [13]

    S. Mor, V. Gosetti, A. Molina-S´anchez, D. Sangalli, S. Achilli, V. F. Agekyan, P. Franceschini, C. Giannetti, L. Sangaletti, and S. Pagliara, Phys. Rev. Res. 3, 043175 (2021)

  14. [14]

    C. J. Sayers, A. Genco, C. Trovatello, S. D. Conte, V. O. Khaustov, J. Cervantes-Villanueva, D. Sangalli, A. Molina- Sanchez, C. Coletti, C. Gadermaier, and G. Cerullo, Nano Lett. 23, 9235 (2023)

  15. [15]

    De Giovannini, H

    U. De Giovannini, H. Huebner, S. A. Sato, and A. Rubio, Phys. Rev. Lett. 125, 136401 (2020)

  16. [16]

    Gerber, S.-L

    S. Gerber, S.-L. Yang, D. Zhu, H. Soifer, J. A. Sobota, S. Rebec, J. J. Lee, T. Jia, B. Moritz, C. Jia, A. Gauthier, Y. Li, D. Leuenberger, Y. Zhang, L. Chaix, W. Li, H. Jang, J.-S. Lee, M. Yi, G. L. Dakovski, S. Song, J. M. Glownia, S. Nelson, K. W. Kim, Y.-D. Chuang, Z. Hussain, R. G. Moore, T. P. Devereaux, W.-S. Lee, P. S. Kirchmann, and Z.-X. Shen, S...

  17. [17]

    P. Hein, S. Jauernik, H. Erk, L. Yang, Y. Qi, Y. Sun, C. Felser, and M. Bauer, Nature Communications 11, 2613 (2020)

  18. [18]

    Suzuki, Y

    T. Suzuki, Y. Shinohara, Y. Lu, M. Watanabe, J. Xu, K. L. Ishikawa, H. Takagi, M. Nohara, N. Katayama, H. Sawa, M. Fujisawa, T. Kanai, J. Itatani, T. Mizokawa, S. Shin, and K. Okazaki, Phys. Rev. B 103, L121105 (2021)

  19. [19]

    H. J. Zeiger, J. Vidal, T. K. Cheng, E. P. Ippen, G. Dresselhaus, and M. S. Dresselhaus, Phys. Rev. B 45, 768 (1992)

  20. [20]

    L. Dhar, J. A. Rogers, and K. A. Nelson, Chem. Rev. 94, 157 (1994)

  21. [21]

    T. E. Stevens, J. Kuhl, and R. Merlin, Phys. Rev. B 65, 144304 (2002)

  22. [22]

    K. G. Nakamura, Y. Shikano, and Y. Kayanuma, Phys. Rev. B 92, 144304 (2015)

  23. [23]

    Caruso and M

    F. Caruso and M. Zacharias, Phys. Rev. B 107, 054102 (2023)

  24. [24]

    G. A. Garrett, T. F. Albrecht, J. F. Whitaker, and R. Merlin, Phys. Rev. Lett. 77, 3661 (1996)

  25. [25]

    Shinohara, K

    Y. Shinohara, K. Yabana, Y. Kawashita, J.-I. Iwata, T. Otobe, and G. F. Bertsch, Phys. Rev. B 82, 155110 (2010)

  26. [26]

    Shinohara, S

    Y. Shinohara, S. A. Sato, K. Yabana, J.-I. Iwata, T. Otobe, and G. F. Bertsch, The Journal of Chemical Physics 137, 22A527 (2012), https://pubs.aip.org/aip/jcp/article-pdf/doi/10.1063/1.4739844/14005795/22a527 1 online.pdf

  27. [27]

    Trovatello, H

    C. Trovatello, H. P. C. Miranda, A. Molina-S´anchez, R. Borrego-Varillas, C. Manzoni, L. Moretti, L. Ganzer, M. Maiuri, J. Wang, D. Dumcenco, A. Kis, L. Wirtz, A. Marini, G. Soavi, A. C. Ferrari, G. Cerullo, D. Sangalli, and S. D. Conte, ACS Nano 14, 5700 (2020), pMID: 32233453, https://doi.org/10.1021/acsnano.0c00309

  28. [28]

    Y.-H. Chan, D. Y. Qiu, F. H. da Jornada, and S. G. Louie, Proceedings of the National Academy of Sciences 118, e1906938118 (2021), https://www.pnas.org/doi/pdf/10.1073/pnas.1906938118

  29. [29]

    M. S. Hybertsen and S. G. Louie, Phys. Rev. B 34, 5390 (1986)

  30. [30]

    Rohlfing and S

    M. Rohlfing and S. G. Louie, Phys. Rev. B 62, 4927 (2000)

  31. [31]

    Deslippe, G

    J. Deslippe, G. Samsonidze, D. A. Strubbe, M. Jain, M. L. Cohen, and S. G. Louie, Computer Physics Communications 183, 1269 (2012)

  32. [32]

    Y.-H. Chan, D. Y. Qiu, F. H. da Jornada, and S. G. Louie, Proceedings of the National Academy of Sciences 120, e2301957120 (2023), https://www.pnas.org/doi/pdf/10.1073/pnas.2301957120

  33. [33]

    J. E. Sipe and A. I. Shkrebtii, Phys. Rev. B 61, 5337 (2000)

  34. [34]

    Chan, Zhenglu Li, and S

    Y.-H. Chan, Zhenglu Li, and S. G. Louie, Phys. Rev. B 112, 024308 (2025)

  35. [35]

    P. Y. Yu and M. Cardona, Fundamentals of Semiconductors: Physics and Materials Properties, 4th ed. (Springer, Berlin, 2010)

  36. [36]

    Loudon, The Quantum Theory of Light (Oxford University Press, Oxford, 2000)

    R. Loudon, The Quantum Theory of Light (Oxford University Press, Oxford, 2000)

  37. [37]

    Reichardt and L

    S. Reichardt and L. Wirtz, Nonadiabatic exciton-phonon coupling in Raman spectroscopy of layered materials, Sci. Adv. 6, eabb5915 (2020) 12