{"id":"2dde12a8-5cbf-49e8-a483-1133bb7ab527","arxiv_id":"2502.01529","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Coherent phonon decoherence in crystals is shown to equal the imaginary part of the nonequilibrium phonon self-energy, and ab initio calculations for Bi and Sb reproduce experimental lifetime trends.","lead":"This paper builds a first-principles theory for how laser-excited atomic vibrations in crystals lose their coherence, and it tests the theory on the semimetals antimony and bismuth. The calculations reproduce the observed dependence of phonon lifetimes on temperature and laser fluence, giving a predictive tool for ultrafast materials experiments.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The fluence validation is underdetermined: the fluence-to-electronic-temperature conversion is confined to the Supplemental Material, and the q=0 intraband rate in Eq. (24) depends on a hand-set γ=15 meV; the claimed agreement is therefore not yet an independent test of the self-energy formula.","rationale":"The reader's weakest assumption—that the photoexcited semimetal can be represented by two quasi-equilibrium distributions with a fluence-to-Tel mapping confined to the Supplemental Material—is exactly the load-bearing point I would stress. The central evidence for the paper is quantitative agreement in Fig. 3, and that evidence depends on converting the experimental control variable (fluence) into the theoretical input (Tel) through an unshown prescription. If that prescription has free parameters, or if the actual electron distribution after the pump is not Fermi-Dirac, the agreement is at least in part a fit rather than a prediction. The additional hand-set Lorentzian broadening γ=15 meV in Eq. (24) is not a minor detail: at q=0 the Fan-Migdal contribution vanishes, so the reported electronic decoherence rate is essentially determined by this tail. The paper's own stated omission of photoexcited phonon populations in the phonon-phonon channel further weakens the decomposition of the total rate into electron-phonon and phonon-phonon contributions. I do not see a fatal internal inconsistency in the formal derivation; the remaining uncertainty is empirical and calibrational. Therefore I retain the conditional verdict rather than moving to acceptance or rejection. The proposed test—replacing the calibrated quasi-equilibrium distribution with the actual nonequilibrium electron distribution for the same fluences—settles the key question of whether the fluence dependence is a genuine prediction of the self-energy framework or an artifact of the auxiliary mapping.","tokens_in":15452,"tokens_out":9381,"duration_ms":92713,"concrete_test":"Recompute Fig. 3(a)-(b) using Eq. (24) with the full non-equilibrium electronic occupation generated by the actual pump excitation (for example, via the time-dependent Boltzmann equation with the reported experimental fluence and pulse parameters), instead of a Fermi-Dirac distribution at a converted electronic temperature, keeping γ = 15 meV and Γ_eff fixed. If the predicted lifetime-versus-fluence curve tracks the experimental points within their error bars, the quasi-equilibrium/Tel mapping is harmless; if it departs from the data, the fluence validation is an artifact of the calibration.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the self-energy framework reproduces the measured fluence dependence is anchored in Fig. 3(a)-(b). The experimental data are plotted against an effective electronic temperature obtained from the experimental fluence through a prescription confined to the Supplemental Material (Sec. III, text after Eq. (24)). This conversion is not shown to be parameter-free in the main text; if the mapping involves an assumed absorbed fraction, electronic heat capacity, or thermalization time, the x-axis can be shifted to force agreement with the computed τ ∝ 1/Tel trend. A second calibration enters inside Eq. (24): at the zone center the Fan-Migdal rate vanishes, and the self-consistent linewidth result Γ_SL is controlled by the Lorentzian tail δγ with γ = 15 meV, asserted as 'representative' rather than computed. The quoted Γ_ep ≈ 0.9 ps^-1 (Bi) and 0.4 ps^-1 (Sb) depend sensitively on this width; a different γ would change the balance between electron-phonon and phonon-phonon contributions and could alter the paper's conclusion about which channel dominates. The manuscript itself acknowledges a third uncontrolled piece when it neglects photoexcited phonon distributions in Γ_pp. Together these calibrated elements leave the central validation underdetermined: the formal framework may be correct, but the comparison as presented does not yet establish it.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript derives quantum kinetic equations for coherent phonons from a non-adiabatic electron-phonon Hamiltonian extended with cubic anharmonicity, and identifies the decoherence rate and frequency renormalization with the imaginary and real parts of the non-equilibrium phonon self-energy (Eqs. 13-15 and 21-23). It then evaluates these expressions for the A1g mode of Bi and Sb using first-principles electron-phonon and phonon-phonon calculations, comparing the resulting coherent-phonon lifetimes with fluence-dependent and temperature-dependent pump-probe experiments (Fig. 3). The paper concludes that both electron-phonon and phonon-phonon scattering contribute to decoherence, with the former controlling the fluence dependence and the latter controlling the temperature dependence.","tokens_in":15696,"tokens_out":9508,"duration_ms":85483,"significance":"If the formal mapping between coherent-phonon decoherence and the non-equilibrium phonon self-energy is correct, this is a valuable bridge between coherent phonon dynamics and existing equilibrium self-energy implementations. The manuscript has clear strengths: the derivation starts from an explicit many-body Hamiltonian, the final expressions reduce to known equilibrium self-energies when distributions are thermalized, the problem of the vanishing Fan-Migdal rate at q=0 is addressed through a self-consistent linewidth scheme, and the computational setup uses widely available first-principles tools. However, the numerical validation currently relies on several external or empirical inputs, so the reported agreement with experiments is not yet an independent confirmation of the formalism. The main text is not self-contained for the key validation steps, and the quantitative claims are stronger than the evidence presented.","major_comments":[{"comment":"The finite zone-center rate Γ_SL_A1g is produced by replacing the energy-conserving delta function with a Lorentzian of width γ = 15 meV, described only as 'representative of the electron linewidths.' Because the Fan-Migdal rate vanishes exactly at q=0, the quoted electron-phonon decoherence rates Γ_ep ≈ 0.9 ps⁻¹ (Bi) and 0.4 ps⁻¹ (Sb) are controlled by this uncomputed width rather than by the ab initio couplings alone. The authors should provide a sensitivity study over γ, or preferably compute the electron linewidths that enter the self-consistent scheme, before the electron-phonon channel can be regarded as quantitatively predictive.","section":"§III, Eq. (24) and Figs. 2(b),(e)"},{"comment":"The x-axis of the fluence comparison is an effective electronic temperature obtained from the experimental fluence through a prescription confined to the Supplemental Material. The main text states that the temperatures are 'chosen to reproduce the experimental conditions,' but it does not show the conversion or list its assumptions, such as absorbed fraction, electronic specific heat, or thermalization time. If any part of this mapping is adjustable, the reported τ ∝ 1/T_el agreement with experiment can be shifted along the x-axis and is not an independent test of Eq. (24). The conversion and its material-specific parameters should be presented in the main text or at least summarized with a sensitivity analysis.","section":"§III, Fig. 3(a)-(b) and text after Eq. (24)"},{"comment":"The phonon-phonon comparison introduces an empirical constant rate Γ_eff_A1g taken from Refs. [50,51], which are analyses of the same experiments whose data are plotted in the figure. Adding this offset makes the absolute agreement partly constructed rather than predicted. The authors should either compute the residual rate from an independent mechanism, clearly label it as a fit parameter, or restrict the validation claim to the temperature slope. A related uncontrolled simplification is the neglect of photoexcited phonon distributions in Γ_pp, which the manuscript acknowledges but does not quantify; this affects the separation of the fluence dependence into an electron-phonon part and a temperature-only phonon-phonon part.","section":"§III, text after Eq. (25) and Figs. 3(c)-(d)"},{"comment":"The displayed Heisenberg equation of motion is incorrect as written. For a coordinate U and conjugate momentum P, the double commutator [U,[U,H]] is proportional to [U,P] and does not yield the acceleration term; the correct form is d²U/dt² = −ℏ⁻²⟨[H,[U,H]]⟩ (equivalently +ℏ⁻²⟨[[U,H],H]⟩), which is what leads to Eq. (8). This appears to be a typographical error, but since Eq. (7) is the stated starting point of the central derivation, it should be corrected and the corresponding steps in the Supplemental Material should be checked.","section":"§II B, Eq. (7)"}],"minor_comments":[{"comment":"The phrases 'robust agreement' and 'good quantitative agreement' overstate the evidence given the empirical constants and calibrated temperatures; 'consistent with' or similar wording would be more proportionate unless the above issues are resolved.","section":"Abstract and Conclusions"},{"comment":"The caption contains a typo: 'ans Sb' should read 'and Sb'.","section":"Fig. 3 caption"},{"comment":"The manuscript depends heavily on the Supplemental Material for both the derivation leading to Eqs. (12)-(15) and the fluence-to-temperature mapping; if the SM is not included with the arXiv posting, the referees and readers cannot assess these steps. The authors should ensure the SM is available and its equations are numbered for cross-reference.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The formal part of the paper is plausible and potentially useful, and I do not see a fundamental error in the self-energy identification. The main weakness is the numerical validation: three externally calibrated ingredients (Tel, γ, Γ_eff) carry much of the agreement, and the fluence-to-temperature conversion is hidden in the SM. I would accept the paper after the authors provide the missing details, add a sensitivity analysis for γ, and temper the claims. The Eq. (7) commutator error should also be fixed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe useful new thing here is the formal bridge: the authors show, from the non-adiabatic nuclear Hamiltonian, that the decoherence rate and frequency shift of a coherent phonon are the imaginary and real parts of the nonequilibrium phonon self-energy (Eqs. 13–15 and 21–23). That is a genuinely useful organizing result, and it correctly reduces to the equilibrium self-energy in the appropriate limits. The second genuinely new piece is the self-consistent linewidth (SL) treatment, which fixes the unphysical vanishing of the Fan–Migdal rate at zone center; without it the calculation would say undamped coherent phonons at Γ, which is obviously wrong. The Bi/Sb comparison is a reasonable first application.\n\nThe soft spots are real but mostly in the validation, not the derivation. The fluence-to-electronic-temperature conversion is confined to the SM, so the x-axis of Fig. 3(a)-(b) cannot be independently checked from the main text. The Lorentzian width γ = 15 meV is asserted as representative, and the quoted Γ_ep values depend on it. The residual rate Γ_eff is added to match the low-temperature offset, and the phonon-phonon calculation neglects photoexcited phonon distributions. The paper acknowledges each of these explicitly, which is good, but together they mean the comparison is calibrated rather than a clean ab initio test.\n\nI would not call the agreement robust; I would call it suggestive. The formal framework is likely to be reused, and the zone-center fix is valuable, but the numbers do not yet close the case. That said, the paper deserves a serious referee. The right referee will ask for the SM derivation, a sensitivity study over γ and the Tel mapping, and a discussion of whether the remaining offsets undermine the channel assignment. It is a useful paper for the ultrafast community, and with revision it could be a solid reference.\n\nRecommendation: send it out, but ask for the sensitivity analysis before acceptance.","headline":"A formally grounded derivation of coherent-phonon decoherence from the nonequilibrium self-energy, with a validation that is suggestive but underdetermined by calibrated inputs.","tokens_in":16267,"tokens_out":2015,"would_cite":false,"duration_ms":18052,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["63.20.kd","63.20.kg"],"model":"deepseek-v4-flash","headline":"Phonon decoherence rates are the imaginary part of a non-equilibrium phonon self-energy.","keywords":["phonon decoherence","coherent phonons","non-equilibrium phonon self-energy","electron-phonon coupling","phonon-phonon scattering","bismuth antimony semimetals","pump-probe spectroscopy","first-principles calculations"],"falsifier":"A time- and angle-resolved photoemission experiment on bismuth at the fluence the paper maps to $T_{\\mathrm{el}}=3000$ K, taken simultaneously with a measurement of the A1g decay rate, would settle the point: if the decay rate does not track the instantaneous occupations of the electronic states near the Fermi surface, the quasi-equilibrium mapping used in Eq. (24) fails.","tokens_in":15213,"feed_emoji":"⚛️","tokens_out":7463,"duration_ms":60615,"temperature":0.7,"pith_summary":"This paper claims that the decay of coherent phonons—the oscillating lattice vibrations set in motion by an ultrafast light pulse—has a concrete microscopic origin: the decoherence rate equals the imaginary part, and the frequency shift the real part, of a non-equilibrium phonon self-energy constructed from electron-phonon and phonon-phonon interactions. The authors derive this result from quantum kinetic equations rather than inserting a phenomenological damping term, and they show that existing first-principles machinery can evaluate the relevant self-energies. Applying it to the A1g optical phonon of the semimetals bismuth and antimony, they reproduce the measured pump-fluence dependence and lattice-temperature dependence of the coherent phonon lifetime with quantitative accuracy. The practical payoff is that coherent phonon lifetimes, a central timescale for light-induced structural control, become a calculable material property rather than an adjustable parameter.","feed_headline":"Coherent phonon lifetimes come from a quantum self-energy","feed_subtitle":"First-principles calculations for Bi and Sb match pump-probe lifetimes across fluence and temperature.","key_machinery":"The central object is the non-equilibrium phonon self-energy $\\Pi_{\\mathbf q\\nu}$—a many-body correction describing how a phonon exchanges energy with electrons and with other phonons. For the electron-phonon channel it is built from the retarded density-density response function (Eq. 15); for the phonon-phonon channel it is summed over three-phonon scattering processes weighted by Bose occupation factors (Eq. 23). To make the electron-phonon self-energy usable for coherent modes at the zone center, the paper adopts the self-consistent linewidth approximation (Eq. 24), which replaces the strict energy-conserving delta with a Lorentzian of width $\\gamma=15$ meV and thereby restores the intraband transitions that the Fan-Migdal approximation misses. These objects are what convert the exact equation of motion into the damped-oscillator form whose solutions are the observed decaying coherent oscillations.","core_discovery":"Starting from a many-body Hamiltonian that avoids the Born-Oppenheimer approximation and couples the phonon displacement to the time-dependent electron density, the paper derives an exact equation of motion for the coherent displacement $U_{\\mathbf q\\nu}$. The equation has the form of a damped driven oscillator, and the damping emerges rather than being assumed: the decoherence rate is $\\Gamma_{\\mathbf q\\nu}=-\\mathrm{Im}\\,\\Pi_{\\mathbf q\\nu}^{\\mathrm{NA}}$ for electron-phonon coupling and $\\Gamma_{\\mathbf q\\nu}=-\\mathrm{Im}\\,\\Pi_{\\mathbf q\\nu}^{\\mathrm{pp}}$ for phonon-phonon coupling, with frequency renormalization given by $2\\omega_{\\mathbf q\\nu}\\mathrm{Re}\\,\\Pi_{\\mathbf q\\nu}$. A technical obstacle is that the standard Fan-Migdal self-energy for the electron-phonon channel vanishes exactly at the zone center, where coherent phonons live; the paper overcomes this with a self-consistent linewidth approximation that includes intraband transitions and yields finite decoherence rates. First-principles calculations for the A1g mode in Bi and Sb then give lifetimes that shrink with rising electronic temperature (pump fluence) and with rising lattice temperature, in agreement with pump-probe experiments; the paper concludes that electron-phonon and phonon-phonon scattering each dominate decoherence in different experimentally accessible regimes.","pith_inferences":["The effective-temperature mapping from pump fluence to a single electronic temperature is the least controlled step in the calculation; I would test it by computing the same lifetimes with a time-resolved non-equilibrium occupation and checking whether the decay rate follows the instantaneous distribution rather than a Fermi-Dirac one.","If the self-energy view is right, decoherence and equilibrium phonon linewidths are the same object only when occupations are thermal; out of equilibrium, stimulated-emission and absorption channels can differ, so a coherent phonon in a strongly pumped material could decay faster or slower than any equilibrium linewidth would suggest.","A natural extension is to multimode coherent states: because self-energy contributions are additive at linear order, the framework predicts a hierarchy of decoherence times for simultaneously excited modes, which multi-color pump-probe experiments could map.","The framework should be testable in insulators and semiconductors with intense mid-infrared pumping, where the electron-phonon channel is weak and anharmonic decay should set the lifetime; a failure of the predicted temperature scaling there would point to missing electron-hole or four-phonon terms."],"forward_implications":["Coherent phonon lifetimes in semimetals can be computed from ground-state electronic structure plus electronic and lattice temperatures, so experiments can be interpreted without a full time-dependent simulation.","The phenomenological damping terms used in earlier coherent-phonon models are replaced by specific self-energy diagrams, giving a physical meaning to each contribution to the decay.","In bismuth and antimony, electron-phonon scattering controls decoherence at high pump fluences and low lattice temperature, while phonon-phonon scattering dominates at high lattice temperature; the two channels have different signatures that experiments can separate.","The Fan-Migdal approximation alone is inadequate for coherent phonons at the Brillouin-zone center, and the self-consistent linewidth correction becomes necessary for any material where a Raman-active or A1g-like mode is studied.","The same formalism extends to other driven solids, so the timescales of light-induced phase transitions and structural switching can be predicted from first principles."],"supporting_citations":[{"why":"Supplies the non-equilibrium many-body Hamiltonian and the Green's-function formulation from which the self-energy expressions for electron-phonon coupling are derived.","marker":"[30]"},{"why":"Supplies the self-consistent linewidth scheme that restores finite phonon decoherence at the zone center by including intraband transitions.","marker":"[62]"},{"why":"Supplies the Wannier-interpolation method for electron-phonon matrix elements on which the first-principles calculations rely.","marker":"[25]"},{"why":"Provides the fluence-dependent experimental coherent-phonon lifetimes for bismuth that the calculations reproduce.","marker":"[17]"},{"why":"Provides experimental data and a first-principles treatment of coherent phonon excitation in antimony that the calculations reproduce.","marker":"[53]"},{"why":"Provides the temperature-dependent experimental lifetimes of the A1g mode in bismuth used as a comparison for the phonon-phonon channel.","marker":"[50]"},{"why":"Provides the temperature-dependent experimental lifetimes of the A1g mode in antimony used as a comparison for the phonon-phonon channel.","marker":"[51]"},{"why":"Supplies the equilibrium Green's-function electron-phonon formalism to which the non-equilibrium self-energy reduces for Fermi-Dirac occupations.","marker":"[35]"},{"why":"Supplies the equilibrium anharmonic phonon self-energy to which the phonon-phonon expression reduces at thermal occupations.","marker":"[45–49]"}],"fun_headline_variants":["Phonon decoherence traced to electron and phonon scattering","First-principles theory nails phonon lifetimes in Bi and Sb","New framework predicts how coherent phonons decay","Why phonons lose coherence: a first-principles answer","Electron or phonon scattering: which kills phonon coherence?"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a photoexcited semimetal behaves as if its electrons have one temperature set by the pump intensity while the lattice keeps its own temperature; if the real electron distribution is not of this two-temperature form, the computed fluence dependence is not a clean test of the self-energy formula.","fun_headline_variants_meta":{"raw":{"variants":["Phonon decoherence traced to electron and phonon scattering","First-principles theory nails phonon lifetimes in Bi and Sb","New framework predicts how coherent phonons decay","Why phonons lose coherence: a first-principles answer","Electron or phonon scattering: which kills phonon coherence?"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000945,"raw_usage":{"total_tokens":4088,"prompt_tokens":1049,"completion_tokens":3039,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":665,"completion_tokens_details":{"reasoning_tokens":2957}},"tokens_in":665,"tokens_out":3039,"duration_ms":20840,"temperature":1.0,"reasoning_tokens":2957,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T15:01:34.850073+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A time- and angle-resolved photoemission experiment on bismuth at the fluence the paper maps to $T_{\\mathrm{el}}=3000$ K, taken simultaneously with a measurement of the A1g decay rate, would settle the point: if the decay rate does not track the instantaneous occupations of the electronic states near the Fermi surface, the quasi-equilibrium mapping used in Eq. (24) fails.","supporting_citations":[{"cited_title":"Cappelluti, Electron-phonon effects on the raman spectrum in Mgb 2, Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the non-equilibrium many-body Hamiltonian and the Green's-function formulation from which the self-energy expressions for electron-phonon coupling are derived."},{"cited_title":"Ishioka and O","cited_arxiv_id":null,"evidence_quote":"Supplies the self-consistent linewidth scheme that restores finite phonon decoherence at the zone center by including intraband transitions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Wannier-interpolation method for electron-phonon matrix elements on which the first-principles calculations rely."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the fluence-dependent experimental coherent-phonon lifetimes for bismuth that the calculations reproduce."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides experimental data and a first-principles treatment of coherent phonon excitation in antimony that the calculations reproduce."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the temperature-dependent experimental lifetimes of the A1g mode in bismuth used as a comparison for the phonon-phonon channel."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the temperature-dependent experimental lifetimes of the A1g mode in antimony used as a comparison for the phonon-phonon channel."},{"cited_title":"Giustino, M","cited_arxiv_id":null,"evidence_quote":"Supplies the equilibrium Green's-function electron-phonon formalism to which the non-equilibrium self-energy reduces for Fermi-Dirac occupations."}],"review_version":1}