{"id":"44e405b2-02c6-4fab-84b0-d470bfd318e3","arxiv_id":"2606.15853","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Ultrafast low-energy electron diffuse scattering shows that photoexcited Td-WTe2 first populates phonons along the W-chain direction via electron–phonon coupling before anharmonic scattering spreads them across the Brillouin zone.","lead":"By bouncing ultrafast low-energy electrons off the surface of the semimetal WTe2, researchers watched how laser-deposited energy first excites atomic vibrations along the tungsten-chain direction and then spreads through the lattice. The measurement maps, in momentum space, a hierarchical phonon cascade that helps separate electron–phonon from phonon–phonon relaxation.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Missing forward model of diffuse intensity leaves open the possibility that the measured Γ–X anisotropy is a phonon structure-factor artifact, not an anisotropic phonon population.","rationale":"The reader's weakest_assumption correctly identifies the key gap: the measured diffuse-intensity anisotropy is interpreted as phonon-population anisotropy without a forward model that accounts for the ULEED structure factor. This is indeed the single most load-bearing concern. If the structure-factor weighting along Γ–X is intrinsically stronger—which is plausible given the soft, out-of-plane-polarized low-energy branches along the W-chain direction—then the central microscopic attribution to anisotropic el–ph coupling would be significantly weakened, and the 'hierarchical cascade' could become largely a detector-response artifact. The paper has all the ingredients to perform the missing calculation (DFT phonons, polarization vectors, experimental geometry), but it stops at color-coding the dispersion by out-of-plane contribution. Thus the verdict should remain CONDITIONAL: the interpretation is plausible and internally consistent, but it hinges on a check that has not been shown. The proposed test—comparing the measured early-time anisotropy to a computed thermal TDS map, or normalizing by the late-time thermal map—would settle the question with existing data and code.","tokens_in":8939,"tokens_out":2958,"duration_ms":37536,"concrete_test":"Compute the one-phonon thermal diffuse scattering (TDS) map from the DFT phonons (Methods A) for an isotropic phonon population n_j(q)=[exp(ħω_j(q)/k_B T)-1]^{-1} at T=30 K, using the same ULEED weighting: I(q) ∝ Σ_j |e_j(q)·Q|^2 / ω_j(q) × (n_j(q)+1/2), with Q along the surface normal. Compare the resulting Γ–X/Γ–Y diffuse-intensity ratio to the measured early-time (0–5 ps) map. Alternatively, normalize the measured early-time diffuse map by the measured late-time (>100 ps) quasi-thermal map; if the Γ–X enhancement persists after this normalization, the anisotropy reflects population dynamics, whereas if it disappears, the phonon structure factor alone explains the observed anisotropy.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that anisotropic electron–phonon coupling preferentially populates finite-momentum phonons along Γ–X—rests on interpreting the measured diffuse-background anisotropy as population anisotropy. But the ULEEDS diffuse intensity is not a direct phonon population: it is weighted by the phonon polarization projection onto the scattering vector (here, predominantly out-of-plane), by the inverse-frequency structure factor, and by detector geometry. The paper explicitly acknowledges these weights (Sec. II, Fig. 4b) but never computes the expected diffuse intensity for an isotropic phonon population. In Td-WTe2, the low-energy branches along Γ–X and Γ–Y differ substantially: the soft ZA/TA modes along the W-chain direction are lower in frequency and have larger out-of-plane displacement components than modes along Γ–Y. If that branch-dependent structure-factor weighting is stronger along Γ–X, then even a perfectly isotropic phonon population would produce a Γ–X-enhanced diffuse map. Consequently, the assignment of the fast 3.5 ps rise along Γ–X to anisotropic el–ph coupling is conditional on a cancellation that is asserted, not demonstrated. The same issue clouds the delayed Γ–Y and zone-center rises: what is called 'redistribution' could partly be the momentum-dependent structure factor of a slowly thermalizing population. This is a missing control, not an internal contradiction; it is directly testable and the necessary DFT inputs are already in the paper.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses ultrafast low-energy electron diffraction and diffuse scattering (ULEEDS) on Td-WTe2 to follow momentum-resolved phonon population dynamics after 1030 nm optical excitation. The main experimental results are: (i) a biexponential increase in the Debye-Waller-derived mean-squared displacement; (ii) a prompt rise of diffuse intensity at the Γ-X zone boundary with a ~3.5 ps time constant; (iii) a delayed rise at Γ-Y (~10 ps); and (iv) a slow accumulation near the zone center (>100 ps). The authors interpret this sequence as a hierarchical relaxation: hot electrons preferentially emit finite-momentum phonons along the W-chain (Γ-X) direction via anisotropic electron-phonon coupling, followed by anharmonic phonon-phonon scattering that redistributes energy across the Brillouin zone and finally into low-frequency zone-center modes. The central claim is that this is a momentum-resolved cascade rather than instantaneous thermalization.","tokens_in":9300,"tokens_out":7184,"duration_ms":76162,"significance":"If established, the result would be a valuable demonstration of momentum-resolved, surface-sensitive probing of phonon thermalization in an anisotropic semimetal, complementing previous work on coherent zone-center phonons. The paper introduces no new theory but leverages a relatively new ULEEDS diffuse-scattering capability. The separation of fast and slow MSD components and the fluence dependence of the slow time constant are useful observations. However, the decisive step—interpreting the diffuse-intensity anisotropy as phonon-population anisotropy—is not backed by a forward calculation of the diffuse scattering. The necessary DFT phonon data are available in the paper (Fig. 4b), making the missing control straightforward. The significance of the paper thus hinges on this quantitative step.","major_comments":[{"comment":"The diffuse ULEED intensity is not a direct phonon population. The paper explicitly states it is structure-factor-weighted and sensitive to out-of-plane displacements and inverse-frequency weighting (Sec. II, paragraphs after Fig. 2; Fig. 4b). Yet the analysis treats the measured Γ-X vs Γ-Y diffuse enhancement as a direct population anisotropy. Because the low-energy acoustic modes along Γ-X in Td-WTe2 are lower in frequency and have different polarization content than those along Γ-Y, an isotropic phonon population could by itself produce a Γ-X anisotropic diffuse pattern through the weighting factors. The authors should compute the expected one-phonon diffuse intensity for an isotropic/thermal phonon distribution using their DFT phonon frequencies and polarization vectors, and compare with the measured maps. This is a missing control, not an internal contradiction, and the necessary in","section":"§II (introductory paragraphs) and Methods C; Fig. 4b"},{"comment":"Even with the structure-factor control in place, the attribution of the fast Γ-X rise specifically to electron-phonon coupling is qualitative. The supporting evidence cited is the elongated Fermi surface and band extrema along Γ-X, and the weak fluence dependence of the fast time constant. These are necessary but not sufficient: an anisotropic phonon density of states or direction-dependent anharmonic decay could also produce a momentum-selective early signal. A quantitative estimate (e.g., q-resolved el-ph matrix elements or an EPC-weighted phonon emission rate from DFT) is needed to back the claim that the Γ-X build-up is dominated by el-ph scattering. Alternatively, a control experiment varying the excitation photon energy or polarization could strengthen the assignment.","section":"§II, 'Along Γ−X... identifies anisotropic el–ph coupling'"}],"minor_comments":[{"comment":"The text reads 'where τ_fast and τ_fast are...'—the second symbol should be τ_slow. Please fix this typo.","section":"Eq. (2)"},{"comment":"The abstract states a '30–100 ps timescale' for the final accumulation, while the text says the zone-center intensity 'saturates only after more than 100 ps.' Please make the reported times consistent.","section":"Abstract vs. §II"},{"comment":"The sentence 'Assuming 1−R for the absorption, 70 % of the beam is absorbed for the applied angle of incidence and p-polarization' is unclear. Specify how the incident-angle projection and the (1−R) factor combine to give the absorbed fluence.","section":"Methods B"},{"comment":"The ROIs used for diffuse analysis are not defined in detail. Please specify the size and position of the integration windows (e.g., relative to the zone boundaries and the Γ point), and describe the background-subtraction procedure more explicitly. Also report fit amplitudes and uncertainties for the single-exponential diffuse rises.","section":"Methods C / Fig. 4"},{"comment":"The statement that the 3.5 ps rise 'is further used as an upper boundary for a Gaussian instrument response function' needs clarification: was the IRF width fixed to 3.5 ps or fitted with that as a bound? How sensitive are the reported MSD time constants to this assumption?","section":"Methods C"},{"comment":"There is a typo: 'This assginment' should be 'This assignment.'","section":"§II"}],"recommendation":"major_revision","confidential_remarks":"The main missing control—a forward model of the diffuse intensity—is feasible and directly testable with the DFT data already in the paper. If the authors can provide this calculation and, if possible, a quantitative el-ph coupling estimate, the central claim would be substantially strengthened. I would not reject because the experimental observations appear solid and the proposed interpretation is plausible, but the current manuscript does not yet establish the anisotropic el-ph assignment over structure-factor artifacts."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a solid experimental paper, and the measurement is genuinely new and worth knowing about. The authors use ULEEDS to track the diffuse background in Td-WTe2 and see a clean sequence: fast rise along Gamma-X, slower along Gamma-Y, and late accumulation near the zone center. That momentum-resolved cascade is not in the earlier WTe2 literature, and the fluence dependence of the MSD time constants supports a two-step picture. The distortion correction, fitting, and error treatment look careful, and the paper is honest about the surface sensitivity and the absence of the coherent interlayer mode.\n\nThe main soft spot is exactly what the stress-test flags: the diffuse intensity is structure-factor-weighted, and the paper never computes what an isotropic phonon population would look like in their geometry. The low-frequency branches along Gamma-X are softer and have larger out-of-plane components, so a q-independent population could easily produce a Gamma-X enhancement by itself. If that weighting explains the 3.5 ps rise, then the el-ph attribution is weaker. This is a missing control, not an internal contradiction. The temporal ordering (Gamma-X before Gamma-Y) is suggestive, and the Fermi-surface phase-space argument is plausible, but it is not a calculation. A forward model using the DFT inputs already in the paper would settle it.\n\nThere's also a minor circularity: they use the 3.5 ps Gamma-X rise as an upper bound for the instrument response in the MSD fit. It doesn't break the analysis because the MSD and diffuse signals are independent in principle, but it should be flagged. No raw data or code was available for independent checking, which is common but limits verification.\n\nOverall, this deserves peer review. It's a new result in a well-studied material, the technique is on solid ground, and the authors should be asked to add the isotropic-population forward model—or at least discuss the structure-factor weighting quantitatively. If the weighting doesn't create the anisotropy, this is a strong paper. If it can't be ruled out, the conclusion should be softened to momentum-dependent detection weighting rather than anisotropic el-ph coupling.","headline":"New momentum-resolved phonon cascade data in WTe2; the el-ph anisotropy claim needs a forward model before I'd sign off.","tokens_in":9794,"tokens_out":2133,"would_cite":true,"duration_ms":23768,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In photoexcited WTe2, energy first flows into phonons along the tungsten-chain axis before spreading across the surface Brillouin zone.","keywords":["T_d-WTe2","ultrafast low-energy electron diffraction","diffuse scattering","phonon dynamics","electron-phonon coupling","phonon-phonon scattering","anisotropic thermalization","layered semimetal"],"falsifier":"Compute the ULEED diffuse-intensity map for a momentum-independent (isotropic) phonon population, using the DFT phonon dispersion, polarization vectors, and the inverse-frequency/out-of-plane structure-factor weighting of the backscattering geometry. If that map already shows more diffuse intensity along Gamma-X than Gamma-Y, then the observed anisotropy would not establish anisotropic electron–phonon coupling.","tokens_in":8872,"feed_emoji":"🔬","tokens_out":10540,"duration_ms":99530,"temperature":0.7,"pith_summary":"The paper sets out to show that energy relaxation in the semimetal T_d-WTe2 after an optical pulse is a momentum-resolved cascade, not a single jump to thermal equilibrium. Using ultrafast low-energy electron diffraction and diffuse scattering, the authors track phonon populations in distinct momentum regions of the surface Brillouin zone. They find a rapid build-up of diffuse intensity along the tungsten-chain (Gamma-X) direction within a few picoseconds, which they attribute to anisotropic electron–phonon coupling during electronic cooling. That is followed by a slower redistribution perpendicular to the chains (Gamma-Y) and a gradual accumulation near the zone centre over 30–100 ps, assigned to anharmonic phonon–phonon scattering. The separation of these steps, supported by the fluence dependence of the slow time constant, identifies a hierarchical pathway in which energy is deposited into selected phonon modes before being shared with the broader lattice bath.","feed_headline":"Phonon energy in WTe2 flows along tungsten chains first","feed_subtitle":"Ultrafast diffraction tracks the cascade: energy first enters tungsten-chain phonons, then spreads across the Brillouin zone.","key_machinery":"Ultrafast low-energy electron diffraction with diffuse-scattering analysis (ULEED) is the technique carrying the argument. Low-energy (90 eV) electrons backscattered from the surface are inelastically scattered by phonons, and the transient diffuse background at a given in-plane momentum transfer reflects the structure-factor-weighted population of phonons with that momentum. Because the scattering vector is nearly perpendicular to the surface, the measurement is preferentially sensitive to phonons with out-of-plane displacement components, and low-frequency modes dominate through an inverse-frequency weighting. This momentum resolution is what allows the authors to separate the Gamma-X, Gam","core_discovery":"After optical excitation, T_d-WTe2 does not thermalize as a uniform lattice bath. The diffuse electron-scattering signal shows that energy is first deposited into finite-momentum phonons along the tungsten-chain (Gamma-X) direction on a few-picosecond time scale, consistent with anisotropic electron–phonon coupling originating from the elongated electron and hole pockets aligned with that axis. Diffuse intensity along the orthogonal Gamma-Y direction rises on a roughly 10 ps time scale, and intensity near the zone centre only saturates after more than 100 ps; the slow component's rate increases with pump fluence. Combined with a biexponential increase of the Debye-Waller mean-squared displac","pith_inferences":["If the interpretation is correct, the technique becomes a direct probe of which phonon momenta are populated first in any low-symmetry or topological material, simply by reading the shape of the diffuse pattern at early times.","The instrument resolution (about 3 ps) means the reported 3.5 ps rise is an upper bound; if a future instrument with sub-picosecond resolution shows the anisotropy is already fully developed within the first pulse, the electron–phonon stage would be even more sharply momentum-selective.","The surface-specific backscattering geometry leaves open how the cascade proceeds in the bulk; comparing these results with bulk-sensitive diffuse probes on the same material could separate surface and interior relaxation channels."],"forward_implications":["If the relaxation is hierarchical, the transient phonon distribution in T_d-WTe2 remains nonthermal for tens of picoseconds, so properties that depend on the phonon spectrum (electrical transport, band renormalization) evolve directionally in time and are not captured by a single temperature.","The early Gamma-X phonon population is a momentum-resolved fingerprint of the anisotropic electron–phonon coupling, which can be used to benchmark first-principles calculations of scattering rates in low-symmetry materials.","The inverse scaling of the slow rise time with fluence supports anharmonic phonon–phonon scattering as the redistribution mechanism, implying that the thermalization time can be tuned by excitation density.","Because the backscattering geometry suppresses interlayer sliding modes and highlights out-of-plane acoustic modes, the measurement gives a surface-specific view that supplements and must be reconciled with bulk-sensitive diffuse x-ray or transmission electron diffraction."],"fun_headline_variants":["WTe2 phonons flow along tungsten chains before spreading","Ultrafast diffraction tracks phonon cascade in WTe2","Energy in WTe2 hits chain-axis phonons first","Phonon energy hierarchy in WTe2: chains first, then lattice","How phonons cool in WTe2: anisotropic cascade tracked"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The measured anisotropy in diffuse intensity is read as anisotropy in the phonon population, which requires that the probe's own sensitivity—the inverse-frequency and out-of-plane polarization weighting—does not intrinsically favor the Gamma-X direction when the phonon population is uniform.","fun_headline_variants_meta":{"raw":{"variants":["WTe2 phonons flow along tungsten chains before spreading","Ultrafast diffraction tracks phonon cascade in WTe2","Energy in WTe2 hits chain-axis phonons first","Phonon energy hierarchy in WTe2: chains first, then lattice","How phonons cool in WTe2: anisotropic cascade tracked"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000555,"raw_usage":{"total_tokens":2470,"prompt_tokens":728,"completion_tokens":1742,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":472,"completion_tokens_details":{"reasoning_tokens":1655}},"tokens_in":472,"tokens_out":1742,"duration_ms":12699,"temperature":1.0,"reasoning_tokens":1655,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T11:14:55.616964+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the ULEED diffuse-intensity map for a momentum-independent (isotropic) phonon population, using the DFT phonon dispersion, polarization vectors, and the inverse-frequency/out-of-plane structure-factor weighting of the backscattering geometry. If that map already shows more diffuse intensity along Gamma-X than Gamma-Y, then the observed anisotropy would not establish anisotropic electron–phonon coupling.","supporting_citations":[],"review_version":1}