{"id":"8c806b00-a207-453e-9932-dee6011eb803","arxiv_id":"2607.04840","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"In CH3NH3PbBr3, femtosecond TA is dominated by photoinduced Coulomb screening; on picosecond scales, lattice vibrations reshape X-ray TA while lattice expansion blueshifts optical TA.","lead":"A first-principles method combining nonequilibrium Bethe-Salpeter theory with ab initio molecular dynamics separates electronic from lattice effects in perovskite pump-probe spectra. It shows which mechanisms dominate X-ray and optical transient absorption on femtosecond versus picosecond timescales in CH3NH3PbBr3.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.5","headline":"Linear recombination of separately computed electronic and lattice TA components with free weights is the softest link in the quantitative mechanism assignment.","rationale":"The Reader correctly isolates the free-weight linear recombination and the Fermi-Dirac assumption as the weakest points. Those are not cosmetic: the strongest claim is a quantitative ranking of four microscopic channels across two timescales and two spectral windows. That ranking is obtained only after the separately computed spectra are scaled and added with free coefficients (and an ad-hoc 0.2 factor). Because the paper already flags the occupation limitation for the optical case and never tests the neglected electron–lattice cross terms, the concern is load-bearing for the quantitative part of the claim while leaving the qualitative femtosecond screening dominance intact. Hence the verdict remains CONDITIONAL; no stronger or weaker adjustment is warranted. The proposed joint AIMD+BSE calculation is a direct, feasible check that would either vindicate or overturn the additivity premise.","tokens_in":12389,"tokens_out":537,"duration_ms":4958,"concrete_test":"Recompute one joint nonequilibrium BSE spectrum on AIMD snapshots that already carry the experimental lattice expansion and the same photoexcited occupations used for the electronic-only curves (ne = 3e18 cm-3). If the joint XTA/optical TA lineshape cannot be recovered (within ~10–15 % of peak amplitudes) by the published free-weight linear combination of the four separate components, the additivity assumption fails and the mechanism ranking must be revised.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim assigns distinct dominant mechanisms (screening vs Pauli; vibrations vs expansion) by energy window and delay. That assignment rests on treating the four contributions as additive and then fitting free linear weights to experiment (Fig. 4b: 0.53/0.04/0.43; Fig. 7b: 0.65/0.30/0.05), plus an unexplained 0.2 scale factor on the vibrational XTA curve in Fig. 4a. The paper never demonstrates that the cross terms (carrier-modified screening on a vibrating/expanded lattice, or non-Fermi-Dirac occupations at 100 ps) remain small. The authors themselves note that the Fermi-Dirac cDFT occupation underestimates the optical negative feature near 2.35 eV (§4.4). If the omitted cross terms or the occupation error are comparable to the fitted weights, the ranking of mechanisms (especially “vibrations dominate XTA, expansion dominates optical TA”) is not uniquely fixed by the calculation.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript develops a first-principles workflow that combines a nonequilibrium Bethe–Salpeter equation (NE-BSE) with ab initio molecular dynamics (AIMD) to separate electronic and lattice contributions to pump–probe transient absorption (TA) of CH3NH3PbBr3. Electronic effects are obtained by inserting photoexcited occupations (from RT-TDDFT at ~2 fs or constrained DFT/Fermi–Dirac at ~100 ps) either into the BSE transition coefficients (Pauli blocking) or into the screened Coulomb interaction (photoinduced screening). Lattice effects are obtained from expanded cells and from spectral averages over AIMD snapshots. Equilibrium X-ray (Br K-edge) and optical spectra match experiment. On the femtosecond scale both XTA and optical TA are reported to be dominated by screening-induced blueshifts of excitonic resonances, with Pauli blocking negligible (XTA) or secondary (optical). On the picosecond scale the authors conclude that XTA is governed by lattice vibrations plus residual screening, whereas optical TA is governed by lattice expansion plus Pauli blocking, supported by spectral-moment analysis and by linear fits of the calculated components to experimental TA lineshapes.","tokens_in":12639,"tokens_out":1305,"duration_ms":9671,"significance":"If the mechanism ranking holds, the work supplies a concrete, transferable protocol for assigning TA features in soft, strongly electron–phonon-coupled materials where electronic and lattice responses are entangled. The selective insertion of occupations into transition coefficients versus W cleanly separates Pauli blocking from screening within a single many-body framework, and the AIMD averaging plus spectral-moment analysis gives a transparent microscopic account of vibrational redistribution. Equilibrium spectra agree well with experiment, and the data are deposited in NOMAD, which strengthens reproducibility. The result is of clear interest to the ultrafast spectroscopy and perovskite communities and goes beyond purely electronic NE-BSE treatments by placing lattice vibrations on equal footing.","major_comments":[{"comment":"The quantitative ranking of mechanisms on the 100 ps scale (XTA: vibrations + screening; optical TA: expansion + Pauli) rests on treating the four contributions as additive and recombining them with free linear weights (Fig. 4b: 0.53/0.04/0.43; Fig. 7b: 0.65/0.30/0.05), plus an unexplained 0.2 scale factor on the vibrational XTA curve in Fig. 4a. The manuscript does not demonstrate that cross terms (screening evaluated on vibrating/expanded lattices, or non-Fermi–Dirac occupations) remain small relative to these weights. Because the authors themselves note that the Fermi–Dirac cDFT occupation underestimates the optical negative feature near ~2.35 eV (§4.4), the fitted ranking is not uniquely fixed by the calculation. A controlled test—e.g., NE-BSE on a few AIMD snapshots with photoexcited occupations, or a sensitivity analysis of the fit coefficients under occupation variations—is needed","section":null},{"comment":"The vibrational XTA component is scaled by 0.2 in Fig. 4a before the linear combination in Fig. 4b. No physical or numerical justification is given for this factor (normalization convention, absolute intensity mismatch, or otherwise). Without it the vibrational weight relative to the electronic component cannot be interpreted, which directly affects the statement that lattice vibrations are essential for the XTA pre-edge and overall lineshape (§4.2, Conclusions).","section":null}],"minor_comments":[{"comment":"Eq. (4.1) is written as a display equation but is never numbered; later text refers to “Eq. 4.1” and “Eq. S4” inconsistently. Number all main-text equations.","section":null},{"comment":"The excitation densities used for femtosecond spectra (3.0×10^20 and 1.2×10^21 cm^-3) are two orders of magnitude above the picosecond experimental density (3.0×10^18 cm^-3). A short statement on why the high-density regime is still representative for the screening mechanism would help the reader.","section":null},{"comment":"Table 1 and Table 2 report spectral moments with different energy windows (13.46–13.49 keV vs 1–6 eV). Explicitly state the integration limits in the table captions and confirm that the relative changes are robust to modest window variations.","section":null},{"comment":"Figure 1 workflow labels “Pauli” and “Screening” are clear, but the main text sometimes uses “photoinduced Coulomb screening” and sometimes only “screening”; a single consistent term would improve readability.","section":null},{"comment":"The Supporting Information is cited for formalism (S1–S2) and spectral moments (S4) but is not available in the review package; ensure SI equations for the nonequilibrium BSE (especially the selective insertion into transition coefficients vs W) are complete and self-contained.","section":null}],"recommendation":"major_revision","confidential_remarks":"The methodological core (NE-BSE + selective occupation insertion + AIMD) is solid and builds cleanly on the authors’ prior work (refs 28–29). The load-bearing weakness is the free linear recombination used to assign mechanism dominance at 100 ps; if the authors can either justify the 0.2 scale factor and show that cross terms are small, or reframe the picosecond conclusions as qualitative rather than quantitative rankings, the paper would be suitable for a high-quality materials/condensed-matter journal. Scope and novelty are appropriate."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The useful takeaway is that this group has extended their nonequilibrium BSE (already in refs 28–29) by averaging AIMD snapshots, then applied the full four-way split—screening vs Pauli, expansion vs vibrations—to both X-ray and optical TA of MAPbBr3 at 2 fs and 100 ps. That is new, and the mechanism ranking they get is concrete and useful: fs spectra are screening-dominated blueshifts; ps XTA is vibrations plus residual screening; ps optical TA is expansion plus Pauli.\n\nWhat they do well is the selective insertion of occupations. Putting photoexcited carriers only into the transition coefficients isolates Pauli blocking; putting them only into the screened Coulomb interaction isolates screening. Equilibrium spectra match experiment, spectral moments quantify the vibrational redistribution, and the data are deposited in NOMAD. The fs results look especially clean: Pauli is negligible in XTA and only modest in optical, while screening produces the observed nonlinear blueshift. That part does not rest on free parameters.\n\nThe soft spot the stress-test flags is real but secondary. On the ps scale they recombine the four separately computed components with free linear weights (0.53/0.04/0.43 for XTA; 0.65/0.30/0.05 for optical) and an unexplained 0.2 scale on the vibrational XTA curve. They also note themselves that the Fermi–Dirac cDFT occupation underestimates the optical negative feature near 2.35 eV. Cross terms (screening on a vibrating lattice, non-thermal occupations) are not computed. So the precise numerical ranking of “vibrations dominate XTA, expansion dominates optical” is not uniquely fixed. Still, the qualitative lineshapes already point the same way before the fit, and the fs conclusions stand without it.\n\nThis is for people who actually interpret pump–probe spectra of soft hybrid perovskites or who build first-principles TA methods. It deserves a serious referee; the free-weight issue is fixable by reporting unweighted residuals or a joint calculation. I would cite the mechanism map and the AIMD+NE-BSE workflow.","headline":"Solid first-principles decomposition of MAPbBr3 TA that cleanly separates screening, Pauli, expansion, and vibrations; the free linear weights are a real but secondary soft spot that does not erase the qualitative mechanism map.","tokens_in":13253,"tokens_out":531,"would_cite":true,"duration_ms":4598,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"First-principles method separates electronic from lattice effects in perovskite transient absorption, showing screening rules femtoseconds while vibrations and expansion rule picoseconds.","keywords":["transient absorption spectrum","electronic and lattice contributions","nonequilibrium Bethe-Salpeter equation","ab initio molecular dynamics","metal halide perovskites","photoinduced Coulomb screening","Pauli blocking","CH3NH3PbBr3"],"falsifier":"A simultaneous X-ray and optical transient-absorption measurement on the same CH3NH3PbBr3 crystal at controlled carrier density and lattice temperature that cannot be reproduced by any linear combination of the four calculated components would falsify the claimed separation.","tokens_in":13261,"feed_emoji":"⚡","tokens_out":664,"duration_ms":4976,"temperature":0.7,"pith_summary":"Soft metal-halide perovskites couple electrons and lattice so tightly that pump-probe spectra mix carrier and thermal effects, making experimental lineshapes hard to assign. This paper builds a first-principles workflow that folds photoexcited carrier distributions into a nonequilibrium Bethe-Salpeter equation and averages spectra over ab-initio molecular-dynamics snapshots, so electronic and lattice pieces can be switched on and off separately. Applied to CH3NH3PbBr3, the calculation shows that at two femtoseconds both X-ray and optical transient absorption are almost pure electronic: photoinduced Coulomb screening weakens the electron-hole attraction and blueshifts excitonic peaks, while Pauli blocking is negligible. At one hundred picoseconds the lattice matters, but differently in each window: vibrations redistribute spectral weight at the Br K-edge, whereas lattice expansion blueshifts the optical edge; residual electronic pieces are still screening (X-ray) or Pauli blocking (optical). The result gives a concrete map of which microscopic mechanism dominates which energy window and delay, so experimentalists can stop treating every derivative-like feature as pure carrier dynamics.","feed_headline":"Screening, not blocking, blueshifts perovskite spectra in femtoseconds","feed_subtitle":"First-principles split shows lattice vibrations and expansion take over only on the picosecond scale","key_machinery":"Nonequilibrium Bethe-Salpeter equation fed by RT-TDDFT or constrained-DFT carrier occupations, combined with AIMD snapshot averaging, that isolates Pauli blocking, photoinduced screening, lattice expansion and vibrations by selective inclusion.","core_discovery":"On the femtosecond scale both X-ray and optical transient absorption of CH3NH3PbBr3 are dominated by photoinduced Coulomb screening that blueshifts excitonic resonances (Pauli blocking is negligible); on the picosecond scale X-ray transient absorption is governed by lattice vibrations plus screening, whereas optical transient absorption is dominated by lattice expansion and Pauli blocking.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Screening not Pauli blocking blueshifts perovskite TA in femtoseconds","Photoinduced Coulomb screening dominates fs perovskite exciton shifts","Electronic screening blueshifts CH3NH3PbBr3 spectra before lattice effects","Lattice vibrations and expansion control perovskite TA only at picoseconds","First-principles splits electronic vs lattice origins of perovskite transient absorption"],"cache_read_input_tokens":128,"weakest_assumption_plain":"Electronic and lattice pieces are treated as independent and can be linearly recombined with free weights to match experiment, while picosecond carriers are assumed to follow a simple Fermi-Dirac distribution.","fun_headline_variants_meta":{"raw":{"variants":["Screening not Pauli blocking blueshifts perovskite TA in femtoseconds","Photoinduced Coulomb screening dominates fs perovskite exciton shifts","Electronic screening blueshifts CH3NH3PbBr3 spectra before lattice effects","Lattice vibrations and expansion control perovskite TA only at picoseconds","First-principles splits electronic vs lattice origins of perovskite transient absorption"]},"model":"grok-4.5","effort":"low","cost_usd":0.004698,"raw_usage":{"total_tokens":1357,"prompt_tokens":765,"num_sources_used":0,"completion_tokens":92,"cost_in_usd_ticks":46980000,"prompt_tokens_details":{"text_tokens":765,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":500,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":765,"tokens_out":92,"duration_ms":4226,"temperature":1.0,"reasoning_tokens":500,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-11T12:43:09.039979+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"A simultaneous X-ray and optical transient-absorption measurement on the same CH3NH3PbBr3 crystal at controlled carrier density and lattice temperature that cannot be reproduced by any linear combination of the four calculated components would falsify the claimed separation.","supporting_citations":[],"review_version":1}