{"id":"cf1a8f5a-87e0-44aa-b114-61496dc518b7","arxiv_id":"1908.02795","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A data reduction framework for ultrafast electron diffuse scattering yields time- and wavevector-resolved phonon branch populations in graphite and a mode-projected electron-phonon coupling of 0.035 ± 0.001 eV^2 for the A'1 phonon.","lead":"Ultrafast electron diffuse scattering measurements on graphite are analyzed with a new method that separates the scattered signal into individual phonon modes and maps how each mode's population changes across the Brillouin zone after laser excitation. The method extracts mode-resolved electron-phonon coupling from a tabletop experiment, opening a path to studying superconductors, charge-density waves, and heat transport without large X-ray facilities.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (B10) evaluates the electronic DOS at 1.55 eV − ℏω_A1', but a vertical 1.55 eV excitation puts the emitting carrier near 0.775 eV; the DOS is roughly doubled, so <g²_e,A1'> could shift by ~2×, outside the quoted ±0.001 eV².","rationale":"We did not find a flaw in the pixel-level reduction; Eq. 7 is well-posed, and the DFPT polarization calculation is a standard tool. The central quantitative claim, however, depends entirely on Appendix B. The reader's weakest assumption (graphene DOS for graphite) is real, but even within the Dirac approximation the energy at which D_e is evaluated appears to be off by a factor of two: the photon creates two carriers, not one. If the check confirms this, ⟨g²_e,A1'⟩ becomes ≈0.08 eV², no longer 'in agreement' at the stated precision; if it does not (e.g., because ℏων is defined as the quasiparticle energy in Eq. B10), then the concern reduces to a labeling problem, but that still needs to be fixed. The verdict remains conditional: the population-dynamics method stands, but the absolute coupling value and its quoted error require a corrected or justified conversion and a systematic-error budget.","tokens_in":18616,"tokens_out":11793,"duration_ms":143916,"concrete_test":"Re-derive the conversion from G_e,A1' to ⟨g²_e,A1'⟩ with D_e evaluated at ℏΩ/2 − ℏω_A1' ≈ 0.605 eV (and, for sensitivity, at the graphite DOS from a full DFT band structure at that energy); if the resulting ⟨g²_e,A1'⟩ differs from 0.035 eV² by more than 50%, the headline value and its agreement claim are not supported. A complementary check is to compare the same conversion against the independently measured A1' hot-electron decay rate in a trARPES experiment with 800 nm excitation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Appendix B is the only place where the measured heating rate is converted to the headline ⟨g²_e,A1'⟩, so Eq. (B10) is load-bearing. Eq. (B10) evaluates D_e at ℏων − ℏω_j,k, with ℏων stated as the 1.55 eV pump photon energy. But an 800 nm photon in a Dirac cone creates an electron at E_k = ℏΩ/2 ≈ 0.775 eV above the Dirac point (and a hole at −0.775 eV); subsequent phonon emission leaves that electron at 0.775 − ℏω_A1' ≈ 0.605 eV. Since D_e(ε) ∝ |ε| (Eq. B9), using D_e(1.38 eV) instead of D_e(0.605 eV) inflates the DOS by ≈2.3×, and because ⟨g²⟩ is inversely proportional to D_e in Eq. (B10), the inferred coupling would be ≈0.08 eV² rather than 0.035 eV². The quoted uncertainty ±0.001 eV² is only the NLM fit covariance; it does not include this energy-argument error, the graphite-vs-graphene DOS approximation, the choice of v_F, or the provenance of Eq. (B10) from the overlapping preprint Ref. 62. The population-decomposition method (Eq. 7) is unaffected; the concern targets the absolute mode-projected coupling value and the claimed agreement with other measurements.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a generally-applicable data reduction method for ultrafast electron diffuse scattering (UEDS) that recovers time- and momentum-resolved phonon populations from measured diffuse intensities and DFPT-computed phonon polarization vectors, without free parameters in the linear inversion (Eq. 7). The method is demonstrated on graphite, where the authors extract transient populations of the in-plane phonon branches and, using a non-thermal lattice model fit, report a mode-projected electron-phonon coupling of <g^2_{e,A'_1}> = 0.035 ± 0.001 eV^2 for the A'_1 phonon at K, obtained from the fitted heating rate via the conversion in Appendix B. The paper claims agreement with trARPES measurements and simulations.","tokens_in":19020,"tokens_out":3584,"duration_ms":40618,"significance":"If the central results hold, the paper establishes UEDS as a tabletop probe of full-Brillouin-zone, mode-resolved phonon dynamics, directly analogous to trARPES for the lattice system. The population-decomposition framework is internally consistent and carefully validated: the linear inversion uses M=44 reflections with no fitted parameters, and the authors explicitly check that time-dependent Debye-Waller changes are negligible. The paper also ships open-source tooling (scikit-ued) and detailed DFPT inputs, which aids reproducibility. The headline coupling value, however, is only as solid as the conversion in Appendix B, and the concerns raised below affect that absolute value and the claimed agreement with other measurements.","major_comments":[{"comment":"The conversion from the fitted heating rate to <g^2_{e,A'_1}> evaluates the electronic DOS D_e at ℏων − ℏω_{j,k}, with ℏων stated to be the 1.55 eV pump photon energy. For an 800 nm excitation in a Dirac-cone material, the photoinjected electron sits approximately 0.775 eV above the Dirac point, not 1.55 eV. After emitting an A'_1 phonon (ℏω ≈ 0.17 eV), the relevant final electron energy is ≈ 0.605 eV. Since D_e(ε) ∝ |ε| (Eq. B9), evaluating at 1.38 eV instead of 0.605 eV overestimates the DOS by ≈ 2.3×, which—because <g^2> is inversely proportional to D_e in Eq. (B10)—underestimates the coupling by the same factor. The corrected value would be ≈ 0.08 eV^2, well outside the quoted ±0.001 eV^2. This error directly challenges the claimed agreement with other techniques. The authors must either correct the energy argument or explicitly justify why D_e should be evaluated at the full photon energy in this two-band picture.","section":"Appendix B, Eq. (B10)"},{"comment":"The conversion chain rests on two unquantified approximations: (i) the graphene 2D linear DOS, Eq. (B9), is taken to represent graphite near the K point, and (ii) Eq. (B10) is adopted from a preprint (Ref. [62]) without derivation or a stated range of validity. The reported uncertainty ±0.001 eV^2 is only the NLM fit covariance; it does not include the systematic error in the DOS approximation (including the choice of v_F) or the uncertainty in the borrowed conversion relation. Since these choices can shift the headline value substantially, the authors should derive Eq. (B10), or cite a peer-reviewed derivation, and provide a sensitivity analysis of <g^2> with respect to the DOS parameters.","section":"Appendix B, Eqs. (B8)–(B10)"}],"minor_comments":[{"comment":"The phrase \"anndmode-dependent\" should read \"and mode-dependent\".","section":"§IV C"},{"comment":"\"Similarily\" should be \"Similarly\".","section":"§IV D"},{"comment":"The symbol ℏων denotes the optical excitation energy, but ων could be confused with a phonon frequency; use a distinct symbol such as E_pump to avoid ambiguity.","section":"Appendix B, Eq. (B10)"},{"comment":"The caption of Fig. 6 states an uncertainty of ≈4×10^-5 on the population values, while Fig. 7 shows the A'_1 population change in arbitrary units; the relationship between these two error statements should be clarified.","section":"§IV E, Fig. 7"},{"comment":"The term \"kinematical scattering\" should be \"kinematic scattering\" for consistency with standard usage.","section":"§II B"}],"recommendation":"major_revision","confidential_remarks":"The population-decomposition method is innovative and likely to be well received. The main concern is the Appendix B conversion, where the energy at which D_e is evaluated appears physically incorrect for a Dirac cone; this is a fixable but load-bearing issue. The paper also leans heavily on an unpublished preprint (Ref. [62]) with acknowledged input from M. X. Na; the editor may wish to monitor the overlap between the two works. With the conversion corrected and systematic uncertainties quantified, the paper would be a strong contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth reading for the population-inversion method; the mode-projected coupling number has a shaky appendix. The core trick--solving a linear system over ~44 Brillouin zones to get branch-resolved delta-n_j,k(t) from UEDS intensity--is genuinely new and parameter-free. That part is solid. They check Debye-Waller effects, use DFPT polarization vectors, apply a non-negativity constraint, and produce clean maps of nonthermal TA and LA populations in graphite. This formalism will be reused by the community.\n\nThe weak link is Appendix B. Equation (B10) converts the fitted heating rate G_e,A'1 into <g^2> using a graphene-like DOS evaluated at 1.55 eV - hbar-omega. For a 1.55 eV pump in a Dirac cone, the emitting electron sits near 0.775 eV, so the DOS should be evaluated near 0.6 eV after phonon emission. D_e is proportional to |E|; that is about a factor 2.3 in D_e and hence a factor ~2 in g^2. The quoted 0.035 +/- 0.001 eV^2 would become ~0.08 eV^2. The error bar is only fit covariance, not model uncertainty. Since agreement with trARPES and simulations is asserted without a quantitative comparison table, this needs fixing.\n\nAlso, the conversion formula comes from an overlapping preprint (Na et al., ref. 62), and M. X. Na is thanked in the acknowledgements. Not a fatal circularity, but it should be acknowledged and, ideally, re-derived or cross-checked against independent data.\n\nMinor: the non-thermal lattice model has several approximations--heat capacities, decay channels--but those are documented. The data and code are not deposited, which is a pity for a methods paper.\n\nBottom line: the population decomposition is a real advance and deserves citation; the absolute coupling value is not yet reliable. This paper should go to peer review with a clear request for a sensitivity analysis and a re-derivation of Appendix B.","headline":"Solid new method for branch-resolved phonon populations; the headline e-ph coupling value rests on a sloppy energy argument in Appendix B that could flip the number by ~2.","tokens_in":19539,"tokens_out":3149,"would_cite":true,"duration_ms":38955,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Diffuse electron scattering yields mode-resolved phonon populations","keywords":["ultrafast electron diffuse scattering","phonon populations","electron-phonon coupling","Brillouin zone","graphite","density functional perturbation theory","nonequilibrium lattice dynamics","one-phonon structure factor"],"falsifier":"Compute the mode-projected coupling using the full three-dimensional graphite band structure rather than the graphene linear density of states, and compare the resulting $\\langle g_{e,A'_1}^2\\rangle$ with $0.035$ eV$^2$; a substantial difference would show the DOS approximation dominates the error. Separately, apply the same linear inversion to static thermal diffuse scattering at a known temperature: the recovered populations should match Bose-Einstein occupations for every branch and wavevector. Any systematic branch-dependent mismatch would falsify the polarization-vector-based structure factors.","tokens_in":18453,"feed_emoji":"⚛️","tokens_out":8433,"duration_ms":80494,"temperature":0.7,"pith_summary":"This paper establishes a quantitative, parameter-free data-reduction method that converts ultrafast electron diffuse scattering (UEDS) patterns into transient phonon populations resolved by phonon branch and wavevector across the entire Brillouin zone. The central claim is that, given intensity data from enough Brillouin zones and density-functional-perturbation-theory (DFPT) phonon polarization vectors, a linear inversion at each reduced wavevector recovers the eight in-plane mode populations; no additional free parameters enter. Demonstrated on photodoped graphite, the method tracks energy flow out of the strongly coupled $A'_1$ optical phonon at K, into mid-zone acoustic modes, and yields a mode-projected electron-phonon coupling strength $\\langle g_{e,A'_1}^2\\rangle = 0.035 \\pm 0.001$ eV$^2$ consistent with other measurements and simulations. If the method holds generally, it puts full-Brillouin-zone, mode-resolved lattice dynamics within reach of a laboratory-scale electron-scattering setup.","feed_headline":"Diffuse electron scattering yields mode-resolved phonon populations","feed_subtitle":"Applied to graphite, it fixes the A'_1 electron-phonon coupling at 0.035 ± 0.001 eV^2 from lab-scale data.","key_machinery":"The load-bearing object is the one-phonon structure factor $|F_{1j}(\\mathbf{q}, t_0)|^2$, computed from DFPT phonon polarization vectors $e_{j,s,\\mathbf{k}}$ and Debye-Waller factors. It maps each phonon mode's contribution to diffuse intensity at every scattering vector, and it varies strongly with mode and wavevector: near $\\Gamma$, longitudinal modes scatter radially and transverse modes azimuthally, and only 37% of the measured wavevectors have a single mode contributing more than half of the one-phonon signal. The paper's procedure uses this mode-dependent contrast to turn the measured diffuse intensity into a solvable linear system per reduced wavevector, then converts the fitted heating rate for the $A'_1$ mode into a mode-projected electron-phonon coupling through the relation $\\hbar/\\tau_{e,j}=2\\pi\\langle g_{ep,j}^2\\rangle D_e(\\hbar\\omega_\\nu-\\hbar\\omega_{j,\\mathbf{k}})$ with a graphene-like linear electronic density of states.","core_discovery":"The paper's central discovery is that the time-dependent diffuse scattering intensity change $\\Delta I(\\mathbf{q},t)$ can be decomposed into branch-resolved population changes $\\Delta n_{j,\\mathbf{k}}(t)$ by solving the linear system $\\mathbf{I}_\\mathbf{k}(t)=\\mathbf{F}_\\mathbf{k}\\,\\mathbf{n}_\\mathbf{k}(t)$ at each reduced wavevector $\\mathbf{k}$, where the design matrix is built from one-phonon structure factors computed with DFPT polarization vectors. In graphite, intensity data from 44 Brillouin zones are used to solve for the eight in-plane phonon modes, with non-negative solutions stable for $|\\mathbf{k}|>0.45$ Å$^{-1}$. The recovered populations show the expected fast population of the $A'_1$ phonon at K and the subsequent transfer of energy into transverse acoustic modes along $\\Gamma$–M. Fitting a non-thermal-lattice model to the $A'_1$ population transients gives coupling constants $G_{e,A'_1}=(6.8\\pm0.3)\\times10^{17}$ W m$^{-3}$ K$^{-1}$, $G_{A'_1,l}=(8.0\\pm0.5)\\times10^{17}$ W m$^{-3}$ K$^{-1}$, and $G_{e,l}=(0.0\\pm6.0)\\times10^{15}$ W m$^{-3}$ K$^{-1}$, from which the mode-projected coupling $\\langle g_{e,A'_1}^2\\rangle=0.035\\pm0.001$ eV$^2$ follows through a standard heating-rate conversion.","pith_inferences":["The quoted $\\pm0.001$ eV$^2$ is likely an underestimate of the total uncertainty: the graphene-DOS conversion and the borrowed heating-rate relation introduce systematic errors not reflected in the fit covariance, so a fair comparison with other methods should quote a wider error bar.","The method's accuracy depends on the strength of the polarization-based selection rules; in materials or geometries where multiple scattering redistributes diffuse intensity between Brillouin zones, the linear inversion will be less reliable unless a correction for secondary Bragg scattering is included.","A natural test of the method's internal consistency would be to compare population dynamics extracted at two different pump fluences: the extracted coupling constants should be fluence-independent, while the transient populations should scale with absorbed energy density.","If generalized to three-dimensional polarization information, the same data-reduction route could in principle supply momentum-resolved phonon populations in anisotropic materials relevant to charge-density-wave and superconducting phases, where mode-resolved coupling is the key unknown."],"forward_implications":["The measured $A'_1$ and acoustic-mode populations provide a quantitative, wavevector-resolved view of the energy flow from strongly coupled optical phonons to mid-zone acoustic modes, including the predicted buildup of TA population at $M$ and $\\frac{1}{3}M$.","Extraction of momentum- and mode-dependent electron-phonon and phonon-phonon coupling constants from a single UEDS experiment is possible, moving beyond the two-temperature model's assumption of internally thermalized phonons.","The same formalism carries over to ultrafast X-ray diffuse scattering with minimal changes, extending full-Brillouin-zone phonon population measurements to X-ray free-electron-laser sources.","Pre-photoexcitation thermal diffuse data can be used within the same framework to infer phonon frequencies, then the time-resolved data give population changes, so the method remains applicable when phonon frequencies are not known a priori.","With improved time resolution the approach could detect early-time populations of the $E_{2g}$ phonon in graphite and test the applicability of Kramers-Heisenberg-Dirac theory to Raman scattering in graphene/graphite."],"supporting_citations":[{"why":"Supplies the prior UEDS measurement platform and the qualitative picture of strongly coupled optical phonons that this paper makes quantitative.","marker":"[22]"},{"why":"Provides the kinematical one-phonon scattering formalism used for the diffuse intensity expression.","marker":"[36]"},{"why":"Gives the phonon-population form of the one-phonon diffuse scattering intensity used in Equation (2).","marker":"[37]"},{"why":"Computes the DFPT phonon frequencies and polarization vectors that build the one-phonon structure factors.","marker":"[29]"},{"why":"Provides the anharmonic decay probabilities used to define the effective heat capacity of the modes into which the A'_1 phonon decays.","marker":"[49]"},{"why":"Supplies the relation between the mode-projected coupling and the electron-phonon heating rate used in the final conversion.","marker":"[62]"},{"why":"One of the time- and angle-resolved photoemission measurements whose reported coupling is consistent with the extracted value.","marker":"[12]"},{"why":"Another time-resolved photoemission determination used as a comparison for the A'_1 coupling.","marker":"[14]"}],"fun_headline_variants":["Mode-resolved phonon coupling from diffuse electron scattering","Pinning A1' electron-phonon coupling in graphite via UEDS","Ultrafast diffuse scattering reveals phonon branch populations","Extracting mode-dependent couplings from phonon transients","Phonon population dynamics mapped across Brillouin zone"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conversion from the fitted heating rate to the headline coupling value assumes that graphite's electronic density of states has the same two-dimensional, linear Dirac form as graphene's, with graphene's Fermi velocity; the quoted $\\pm0.001$ eV$^2$ uncertainty covers only the fit covariance, not this electronic-structure approximation.","fun_headline_variants_meta":{"raw":{"variants":["Mode-resolved phonon coupling from diffuse electron scattering","Pinning A1' electron-phonon coupling in graphite via UEDS","Ultrafast diffuse scattering reveals phonon branch populations","Extracting mode-dependent couplings from phonon transients","Phonon population dynamics mapped across Brillouin zone"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00036,"raw_usage":{"total_tokens":2026,"prompt_tokens":1102,"completion_tokens":924,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":718,"completion_tokens_details":{"reasoning_tokens":843}},"tokens_in":718,"tokens_out":924,"duration_ms":10241,"temperature":1.0,"reasoning_tokens":843,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:34:22.212286+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the mode-projected coupling using the full three-dimensional graphite band structure rather than the graphene linear density of states, and compare the resulting $\\langle g_{e,A'_1}^2\\rangle$ with $0.035$ eV$^2$; a substantial difference would show the DOS approximation dominates the error. Separately, apply the same linear inversion to static thermal diffuse scattering at a known temperature: the recovered populations should match Bose-Einstein occupations for every branch and wavevector. Any systematic branch-dependent mismatch would falsify the polarization-vector-based structure factors.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the prior UEDS measurement platform and the qualitative picture of strongly coupled optical phonons that this paper makes quantitative."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the kinematical one-phonon scattering formalism used for the diffuse intensity expression."},{"cited_title":"Xu \\ and\\ author T","cited_arxiv_id":null,"evidence_quote":"Gives the phonon-population form of the one-phonon diffuse scattering intensity used in Equation (2)."},{"cited_title":"Giannozzi , author O","cited_arxiv_id":null,"evidence_quote":"Computes the DFPT phonon frequencies and polarization vectors that build the one-phonon structure factors."},{"cited_title":"Bonini , author M","cited_arxiv_id":null,"evidence_quote":"Provides the anharmonic decay probabilities used to define the effective heat capacity of the modes into which the A'_1 phonon decays."},{"cited_title":"Direct determination of mode-projected electron-phonon coupling in the time-domain","cited_arxiv_id":"1902.05572","evidence_quote":"Supplies the relation between the mode-projected coupling and the electron-phonon heating rate used in the final conversion."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"One of the time- and angle-resolved photoemission measurements whose reported coupling is consistent with the extracted value."},{"cited_title":"Stange , author C","cited_arxiv_id":null,"evidence_quote":"Another time-resolved photoemission determination used as a comparison for the A'_1 coupling."}],"review_version":1}