{"id":"e7a53267-c556-47d7-a73c-2cb275e8d8de","arxiv_id":"2509.07612","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A cDFT/RT-TDDFT + non-equilibrium BSE method reproduces transient absorption spectra in three solids and attributes the blue shift to Coulomb screening, the red shift to lattice expansion, with Pauli blocking minor.","lead":"This paper presents a first-principles computational workflow for simulating ultrafast pump-probe X-ray spectroscopy in solids, combining constrained DFT, real-time TDDFT, and the Bethe-Salpeter equation, and applies it to WSe2, CsPbBr3, and TiO2. It matters because it attributes the transient spectral blue shift to Coulomb screening and the red shift to lattice heating, giving experimentalists a quantitative interpretation and a way to tune exciton resonances.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Thermal red shift may be driven by fitted temperature-dependent scissor shifts rather than lattice expansion","rationale":"The central novel claim is the decomposition of transient absorption into electronic (Coulomb screening, Pauli blocking) and thermal (lattice expansion) contributions. The electronic part is credible and partially supported by independent experiments and prior work (Ref. 26). The thermal part, however, rests on Eq. 2 computed with different scissor shifts at high T. Since the paper advertises 'first-principles,' any fitted parameter used only at high T weakens that claim. The scissor changes are ~0.4–1.7 eV, far larger than the displayed red shifts, so unless the high-T scissor shift is itself derived from the same lattice expansion—which the text does not explain—the computed ΔImε(ω,T) conflates two effects. The most direct, decisive test is to rerun the thermal calculation with a fixed scissor shift. This concern does not invalidate the electronic conclusions, but it makes the thermal conclusion conditional on the scissor treatment. The reader's weakest assumption identified exactly this issue, and my analysis agrees. Thus the verdict remains CONDITIONAL, unchanged from the reader, pending the clarifying calculation/report.","tokens_in":15131,"tokens_out":4520,"duration_ms":48636,"concrete_test":"Recompute the high-temperature spectra (T = 350–798 K) with the scissor shifts fixed at their 0 K values (6.9, 163.5, 106.2 eV for WSe2, CsPbBr3, TiO2), using only the expanded lattice geometries. If the red-shift peaks in Figs. 1b–3b largely disappear, the claim that lattice expansion causes the red shift is an artifact of the scissor adjustment; if they persist, the concern is resolved. The authors should also report how the high-T scissor shifts were determined.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The thermal contribution in Eq. 2 is computed using different scissor shifts at 0 K and at elevated T (SI Table I: WSe2 6.9→6.5 eV; CsPbBr3 163.5→162.0 eV; TiO2 106.2→104.5 eV). These shifts are far larger than the reported red shifts (Figs. 1b–3b) and the text does not explain how the high-temperature scissor values are obtained. If they are fitted to experiment, then ΔImε(ω,T) in Eq. 2 is not a first-principles prediction of thermal effects; the red shift could originate from the imposed scissor correction, not from homogeneous lattice expansion. This directly undermines central claim (iii) and the paper's 'first-principles' characterization.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a first-principles workflow for pump-probe transient absorption (TA) spectroscopy by combining constrained DFT (cDFT) and real-time TDDFT with a non-equilibrium Bethe-Salpeter equation (BSE) treatment. The method is implemented in the all-electron code exciting and applied to core-level TA spectra of three materials: the Se M-edge of WSe2, the Br K-edge of CsPbBr3, and the Ti K-edge of anatase TiO2. The electronic contribution is decomposed into Pauli blocking and Coulomb screening, and a thermal contribution is added by computing equilibrium-like spectra on homogeneously expanded lattices at elevated temperatures. The central claims are that (i) photoinduced Coulomb screening dominates the electronic response and produces a blue shift of core-exciton resonances, (ii) Pauli blocking is minor, and (iii) thermal lattice expansion produces a red shift. The paper also explores the dependence of TA spectra on excitation density, carrier distribution temperature, pump polarization, and pump wavelength.","tokens_in":15373,"tokens_out":4845,"duration_ms":47890,"significance":"If the central thermal claim is established, this would be a useful contribution: it provides an all-electron implementation, a transparent decomposition of electronic effects, and application to three distinct material classes with core-level sensitivity. The electronic part of the calculation is internally coherent: excited-state occupations enter the RPA screening and dipole weights, and a common scissor shift cancels in the difference spectra. The authors also provide data availability via NOMAD, which supports reproducibility. The main significance is conditional, however, because the thermal contribution is computed with different scissor shifts at elevated temperature, and the text does not explain their origin. The paper's 'first-principles' characterization is therefore not yet fully supported for the thermal component, which is a central claim rather than a peripheral detail.","major_comments":[{"comment":"The thermal contribution is defined as Imε(ω,T) − Imε_eq(ω), but SI Table I lists different scissor shifts at high temperature: WSe2 6.9→6.5 eV, CsPbBr3 163.5→162.0 eV, and TiO2 106.2→104.5 eV. The text does not state how the high-temperature scissor values are obtained. If they are fitted or adjusted to experimental high-temperature spectra, then the red shifts in Figs. 1b–3b are not a prediction from homogeneous lattice expansion; the scissor changes alone are 0.4–1.7 eV, comparable to or larger than the spectral shifts shown. The authors must either show that the high-temperature scissor is determined from first principles (e.g., from a quasiparticle calculation on the expanded lattice) or repeat the thermal calculation with the same scissor as the equilibrium reference, so that ΔImε(ω,T) isolates the lattice-expansion effect.","section":"Eq. (2) and SI Table I"},{"comment":"The equilibrium spectra are aligned to experiment by material-specific scissor shifts (SI Table I), so the reported 'good agreement' and 'excellent agreement' with experiment is partly by construction for absolute peak positions. This is not by itself fatal for difference spectra, but it weakens the claim of a 'quantitative and predictive framework' and needs to be stated explicitly. The authors should clarify which conclusions depend on the fitted absolute energy scale and which are robust to the scissor choice. Ideally, they should show that the electronic TA spectra remain unchanged when a common scissor is used for the equilibrium and non-equilibrium calculations.","section":"Methods and Figs. 1–3"},{"comment":"The thermal calculation is described only as 'homogeneous lattice expansion at this respective temperature.' The lattice parameters used for WSe2 at 700 K, CsPbBr3 at 350 K, and TiO2 at 798 K are not given, nor is the expansion recipe (linear versus volumetric, and whether the thermal expansion coefficients are experimental or computed). Without these details the thermal calculation is not reproducible, and one cannot separate the effect of lattice expansion from the changed scissor shifts. Report the expanded structures, the expansion procedure, and the scissor values derived from those structures.","section":"Method (IV) and Supplementary Section II"}],"minor_comments":[{"comment":"The text refers to Fig. 3e and Fig. 3f when discussing carrier-temperature dependence in CsPbBr3; the corresponding panels for CsPbBr3 appear to belong to Fig. 2. Correct the cross-references.","section":"CsPbBr3 section"},{"comment":"Typo: 'absprption' should be 'absorption' in the sentence beginning 'The top panels of Figs. 1a-b also show the impact of photoexcited carriers on the absprption spectra.'","section":"WSe2 section"},{"comment":"Extra punctuation in 'as illustrated in Figs. 1c-d,.'—remove the comma after the figure reference.","section":"Figs. 1c-d caption"},{"comment":"The figures label the ordinate as 'Im M' while the text uses Imε. Use a consistent notation or define the conversion explicitly.","section":"Figure notation"},{"comment":"The BSE broadening for CsPbBr3 is listed as 2.4 eV, which is very large. Please justify this value and discuss its effect on the TA amplitudes and peak positions.","section":"SI Table I"},{"comment":"Spin-orbit coupling is neglected. For WSe2, justify this approximation for the Se M4,5-edge spectra, where spin-orbit effects may be relevant.","section":"Methods"}],"recommendation":"major_revision","confidential_remarks":"To the editor: The manuscript is in scope and the electronic decomposition is valuable. The main risk is the thermal attribution: the temperature-dependent scissor shifts are unexplained and could dominate the computed red shift, which would undermine central claim (iii). I recommend major revision rather than rejection because the issue is local and fixable: the authors can recompute the thermal contribution with a consistent scissor or provide a first-principles derivation of the high-temperature scissor values. The data-availability statement and implementation in exciting are positive features."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this one. The core electronic result — that photoexcited carriers blue-shift core excitons mainly through Coulomb screening, with Pauli blocking playing a minor role — is well supported and consistent with the group's earlier ZnO work. The genuinely new piece is the workflow: cDFT and RT-TDDFT carrier occupations fed into a non-equilibrium BSE implemented in the exciting code, demonstrated for TMDC, halide perovskite, and oxide. Data are available in NOMAD. That is a useful capability, and the authors do a careful job of separating screening from Pauli blocking.\n\nThe soft spot is the thermal claim. The thermal TA is defined as the difference between the spectrum on the thermally expanded lattice and the equilibrium spectrum. But SI Table I lists different scissor shifts at high T (WSe2 6.9→6.5 eV; CsPbBr3 163.5→162.0 eV; TiO2 106.2→104.5 eV). Since the equilibrium spectrum is aligned to experiment with a fitted scissor, and the high-T spectrum uses a different fitted value, part of the computed red shift could be baked into the fit rather than arising from lattice expansion. The text does not explain how the high-T scissor values are obtained. That does not kill the electronic part, but it does undercut the paper's \"first-principles\" framing for the thermal decomposition and the specific claim that lattice expansion drives the red shift. The authors need to show that a fixed scissor with the expanded lattice still gives a red shift, or justify the T-dependent scissor from first principles.\n\nMinor points: the Lorentzian broadening and the low-T scissor shifts are free parameters, but that is normal for BSE core-spectrum calculations. The T-dependent scissor discrepancy deserves a sentence in the main text, not just the SI.\n\nOverall, the electronic machinery is credible and worth citing; the thermal attribution needs a supplementary calculation. I would send it to peer review, because the method and three-material demonstration are genuinely useful and the thermal issue is addressable.","headline":"Worth reading: the electronic transient-absorption machinery is credible and useful, but the thermal red-shift attribution rests on an unexplained change of scissor shift between 0 K and high T, so that claim needs a fixed-scissor check.","tokens_in":15873,"tokens_out":1347,"would_cite":true,"duration_ms":14014,"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":"The paper claims that photoinduced Coulomb screening, not Pauli blocking, drives the blue shift of core-exciton resonances in transient X-ray absorption spectra, while thermal lattice expansion shifts them red.","keywords":["pump-probe spectroscopy","transient absorption","Bethe-Salpeter equation","Coulomb screening","core excitons","WSe2","CsPbBr3","anatase TiO2"],"falsifier":"Recompute the picosecond thermal spectra with the scissor shift held at its 0 K value while keeping the expanded lattice; if the red shift disappears, then Eq. (2) is caused by the fitted scissor change rather than by lattice expansion.","tokens_in":15022,"feed_emoji":"⚛️","tokens_out":5416,"duration_ms":50859,"temperature":0.7,"pith_summary":"The paper tries to establish a first-principles route to simulating ultrafast pump-probe spectroscopy and to identifying the physical mechanisms behind transient X-ray absorption in solids. It computes excited-state carrier distributions with constrained DFT for thermalized carriers and real-time TDDFT for non-thermal carriers, then feeds them into a non-equilibrium Bethe-Salpeter equation for the absorption spectrum. Applied to WSe2, CsPbBr3, and anatase TiO2, the method reproduces measured transient spectra and separates the response into electronic and thermal parts. The central finding is that photoinduced Coulomb screening, not Pauli blocking, drives the blue shift of core-exciton resonances, while lattice heating shifts the spectra red.","feed_headline":"Coulomb screening, not Pauli blocking, shifts core excitons","feed_subtitle":"A first-principles route to pump-probe X-ray spectra separates electronic and thermal shifts in three solids.","key_machinery":"The load-bearing object is the non-equilibrium Bethe-Salpeter equation, the two-particle equation whose eigenvalues give exciton energies and whose macroscopic dielectric function gives absorption. The novelty is that excited-state occupations enter selectively: occupations modify the screened Coulomb interaction to produce Coulomb screening, and the dipole transition coefficients to produce Pauli blocking, so the two mechanisms can be computed separately. The carrier distributions themselves are generated by constrained DFT for thermalized picosecond populations and real-time TDDFT for femtosecond non-thermal populations, in an all-electron full-potential implementation that treats core and","core_discovery":"For three materials spanning different solid classes, the paper claims that the transient absorption spectrum at a core edge is quantitatively captured by a non-equilibrium Bethe-Salpeter calculation whose input occupations come from either thermal Fermi-Dirac constrained DFT or real-time TDDFT. In every case the dominant electronic contribution is Coulomb screening by photoexcited carriers: it weakens the electron-hole attraction, reduces exciton binding energies (in CsPbBr3 from about 360 meV to 90 meV at an excitation density of 3.2e21 cm^-3), and blue-shifts the resonances. Pauli blocking contributes negligibly on picosecond delays and only mildly on femtosecond delays. Heating the latti","pith_inferences":["If Coulomb screening is the dominant measurable response, transient resonance shifts could serve as an in-situ probe of excited-carrier density and exciton binding energy in materials where direct transport measurements are difficult.","The paper's mechanism suggests a design rule: choosing a shorter pump wavelength increases carrier polarizability and screening, so pump color can tune the transient response independently of pump fluence.","The saturation of screening with carrier temperature hints that experimental temperature scans could directly test whether carrier delocalization or carrier density is the limiting factor for resonance shifts.","Because the method separates Pauli blocking from screening, it could be extended to valence-edge or infrared transient spectra where Pauli blocking is often assumed to be central; the present results suggest it may often be secondary."],"forward_implications":["In all three studied materials, transient X-ray absorption at core edges can be predicted from first principles, and the blue shift in measured spectra can be read as a signature of photoinduced Coulomb screening rather than Pauli blocking.","The drop in exciton binding energy with excitation density provides a quantitative link between pump fluence and resonance energy shift.","Thermal red shifts at picosecond delays are attributed to lattice expansion, so the method separates electronic carrier effects from lattice heating in the same spectrum.","Pump polarization and pump wavelength allow femtosecond-scale tuning of exciton resonances, raising the possibility of controlling screening without altering the material.","Carrier delocalization at higher carrier temperatures enhances screening, but this enhancement saturates around 500 K in all three materials."],"supporting_citations":[{"why":"Supplies the non-equilibrium Bethe-Salpeter equation formalism for transient photoabsorption that this paper extends to real materials.","marker":"[17]"},{"why":"Demonstrates the combined theoretical-experimental approach for nonequilibrium optical properties in bulk silicon.","marker":"[18]"},{"why":"Provides the all-electron full-potential code in which the non-equilibrium BSE and TDDFT methods are implemented.","marker":"[25]"},{"why":"Earlier demonstration that the same machinery reproduces experimental transient absorption at a core edge in ZnO.","marker":"[26]"},{"why":"Experimental WSe2 extreme-ultraviolet transient absorption spectra used as the reference for the first material.","marker":"[27]"},{"why":"Experimental CsPbBr3 spectra used as the reference for the halide perovskite case.","marker":"[49]"},{"why":"Experimental anatase TiO2 femtosecond X-ray absorption data used as the reference for the oxide case.","marker":"[52]"},{"why":"Supplies the all-electron real-time TDDFT projection method used to obtain non-thermal carrier distributions.","marker":"[55]"}],"fun_headline_variants":["Screening, not Pauli blocking, drives exciton shifts","Coulomb screening blue-shifts excitons in three solids","Pump-probe theory: screening dominates exciton dynamics","Why excitons blue-shift: screening, not Pauli blocking","Excitons shift by screening, not blocking"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The thermal contribution is attributed entirely to homogeneous lattice expansion, yet the high-temperature spectra also use fitted scissor shifts that differ from the 0 K values; if those scissor changes drive the computed red shift, the paper's thermal mechanism is not established.","fun_headline_variants_meta":{"raw":{"variants":["Screening, not Pauli blocking, drives exciton shifts","Coulomb screening blue-shifts excitons in three solids","Pump-probe theory: screening dominates exciton dynamics","Why excitons blue-shift: screening, not Pauli blocking","Excitons shift by screening, not blocking"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000177,"raw_usage":{"total_tokens":1141,"prompt_tokens":766,"completion_tokens":375,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":510,"completion_tokens_details":{"reasoning_tokens":293}},"tokens_in":510,"tokens_out":375,"duration_ms":4864,"temperature":1.0,"reasoning_tokens":293,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T21:57:04.533770+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the picosecond thermal spectra with the scissor shift held at its 0 K value while keeping the expanded lattice; if the red shift disappears, then Eq. (2) is caused by the fitted scissor change rather than by lattice expansion.","supporting_citations":[{"cited_title":"Perfetto, D","cited_arxiv_id":null,"evidence_quote":"Supplies the non-equilibrium Bethe-Salpeter equation formalism for transient photoabsorption that this paper extends to real materials."},{"cited_title":"Sangalli, S","cited_arxiv_id":null,"evidence_quote":"Demonstrates the combined theoretical-experimental approach for nonequilibrium optical properties in bulk silicon."},{"cited_title":"Gulans, S","cited_arxiv_id":null,"evidence_quote":"Provides the all-electron full-potential code in which the non-equilibrium BSE and TDDFT methods are implemented."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier demonstration that the same machinery reproduces experimental transient absorption at a core edge in ZnO."},{"cited_title":"Carrier and Phonon Dynamics in Multilayer WSe2 captured by Extreme Ultraviolet Transient Absorption Spectroscopy","cited_arxiv_id":"2211.08731","evidence_quote":"Experimental WSe2 extreme-ultraviolet transient absorption spectra used as the reference for the first material."},{"cited_title":"Cannelli, N","cited_arxiv_id":null,"evidence_quote":"Experimental CsPbBr3 spectra used as the reference for the halide perovskite case."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Experimental anatase TiO2 femtosecond X-ray absorption data used as the reference for the oxide case."},{"cited_title":"Rodrigues Pela and C","cited_arxiv_id":null,"evidence_quote":"Supplies the all-electron real-time TDDFT projection method used to obtain non-thermal carrier distributions."}],"review_version":1}