{"id":"dde72116-e52b-4c88-8bd7-219240dee068","arxiv_id":"2502.08521","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Cation substitution tunes exciton-phonon coupling in 2D perovskites, with stronger polaronic coupling correlating with reduced thermal dephasing.","lead":"Researchers measured how strongly excitons couple to lattice vibrations in three nearly identical 2D perovskite materials and found that swapping one atom on the organic spacer molecule tunes that coupling. The material with the strongest coupling also showed the least thermal decoherence, supporting the idea that lattice dressing can shield excitons from scattering.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Correlation trend rests on fitted Huang-Rhys parameters with no reported uncertainties; the M1 two-PES fit is underdetermined.","rationale":"The reader's weakest assumption identifies the same vulnerable link: the fitting model and the extracted S values. I read the paper in good faith and find it a well-conceived spectroscopy study with genuinely suggestive trends; the risk is not internal inconsistency but underdetermination of the extracted parameter. The absence of error bars is decisive because the headline correlation is monotonic across only three materials, and the smallest reported differences appear to be between PEA and Cl-PEA. The two-PES extension in Eq. 4 makes the M1 fits even harder to constrain, since the ad hoc weighting absorbs intensity information from the absorption spectrum. The 2DES dephasing comparison uses only two materials and leans on a prior published value, so it should be labeled preliminary rather than supporting a general protective mechanism. A concrete bootstrap or repeat-measurement uncertainty analysis would settle whether the trend survives. Therefore I maintain the reader's CONDITIONAL verdict rather than moving to ACCEPT or REJECT.","tokens_in":17897,"tokens_out":1135,"duration_ms":11600,"concrete_test":"Perform a global nonlinear least-squares fit of Eqs. 2-4 to the vertical spectral cuts, bootstrap-resampling the noise from the Fourier-domain residuals or from repeated independent measurements, and report S +/- uncertainty for both M1 and M2. If the uncertainty intervals for F-PEA, PEA, and Cl-PEA overlap, the three-point trend in Figs. 4(g)-(h) is not statistically supported.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is a three-point correlation between octahedral distortion (bond angle variance) and exciton-lattice displacement extracted from RISRS lineshape fits. The load-bearing link is the assumption that the displaced-harmonic-oscillator model in Eqs. 2-4 returns reliable Huang-Rhys parameters. The paper reports no error bars, no fit residuals, and no convergence or uniqueness analysis for the S values, and the lowest-energy M1 mode explicitly 'deviates from the two-peak structure' and requires an ad hoc two-PES expansion (Eq. 4) whose relative weights are taken from absorption intensities rather than fit freely. If anharmonicity or mode mixing shifts the fitted S values, a three-point trend with no uncertainties could reorder. The 2DES polaronic-protection comparison is separately a two-material comparison that relies on a self-cited prior value for PEA, so the inverse correlation between S and alpha_LO is not established within this manuscript.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports a combined crystallographic, resonant impulsive stimulated Raman scattering (RISRS), and two-dimensional electronic spectroscopy (2DES) study of (PEA)2PbI4 and its 4-fluoro and 4-chloro substituted derivatives. From RISRS beat maps the authors extract Huang-Rhys parameters S by fitting probe-energy-dependent modulation lineshapes with a displaced harmonic oscillator model (Eqs. 2-3), extended for the lowest-energy mode with a weighted sum of two excited-state potential energy surfaces (Eq. 4). They report a three-point trend in which the exciton-lattice displacement correlates with octahedral bond angle variance, with F-PEA largest and Cl-PEA smallest, and they interpret this as evidence that organic cation substitution tunes polaronic coupling. In addition, temperature-dependent 2DES linewidths for F-PEA are fitted to a single-phonon activation model (Eq. 7), giving a smaller effective exciton-phonon scattering parameter than previously reported for PEA, which the authors interpret as support for polaronic protection against thermal dephasing.","tokens_in":18185,"tokens_out":4630,"duration_ms":50513,"significance":"If the central claims are quantitatively reliable, the paper would establish a practical design rule for 2D metal-halide perovskites: organic cation substitution can tune the degree of exciton-lattice dressing, and stronger dressing may suppress thermal dephasing. The study combines structural analysis, time-domain vibrational spectroscopy, and nonlinear optical linewidth measurements in a coherent way, and the use of chemically similar cations that preserve the average electronic structure is a sound design. The authors are also transparent about the functional forms used in the fits. However, the central claims currently rest on fitted Huang-Rhys parameters without reported uncertainties and on a two-material comparison for the dephasing claim, so the quantitative support is not yet commensurate with the strength of the conclusions.","major_comments":[{"comment":"The central correlation between lattice displacement and bond angle variance is based on Huang-Rhys parameters extracted from nonlinear least-squares fits of the RISRS lineshapes, yet the manuscript reports no uncertainties on the fitted S values and no goodness-of-fit, residual, or covariance diagnostics. Since the ordering F-PEA > PEA > Cl-PEA is the main experimental result, a three-point trend without error bars is not statistically supported. I request that the authors report confidence intervals for each S (e.g., from fit covariance or a bootstrap over the beating-map noise), and that they state explicitly which parameters were fixed and which were free. The text says that γ and ωeg can be approximated from linear spectra and therefore leave S as the 'only true fitting parameter'; the sensitivity of S to the chosen fixed values should also be quantified.","section":"§II, Eqs. (2)-(4), Figs. 4(g)-(h)"},{"comment":"For the lowest-energy mode M1, the manuscript states that the lineshape 'deviates from the two-peak structure' and requires a weighted sum of two excited-state potential energy surfaces, with weights A and B fixed from absorption intensities. This is a strong modelling assumption. There is no test of whether adding a third PES, allowing the weights to vary, or including anharmonicity or excited-state absorption would change the extracted S values. The M1 Huang-Rhys parameters, including those for XB, are therefore model-dependent in a way that is not quantified. I request a sensitivity analysis: at minimum, vary the Eq. (4) weights over a reasonable range and report how S changes, and provide the fit residuals for the M1 lineshapes.","section":"§II, Eq. (4) and Fig. 4"},{"comment":"The polaronic-protection claim rests on a two-material comparison of the effective exciton-phonon scattering parameter αLO: αLO = 1.58 meV for F-PEA measured here versus αLO = 33 meV for PEA from Ref. [86]. Two points cannot establish an inverse correlation between S and αLO, and the comparison is not self-contained because the reference value comes from a different study and no uncertainty is reported for either αLO. Furthermore, the manuscript does not specify which Huang-Rhys parameter (XA/M2, XB/M1, or another) is used in the comparison, while the 2DES linewidths are reported for exciton B (Fig. S4). The authors should either add the Cl-PEA dephasing data to complete the three-point trend, or explicitly frame the inverse correlation as a tentative two-point observation, and they should state which S value is being compared.","section":"§III, Eq. (7), Fig. S4"},{"comment":"The methods for the dephasing analysis need clarification: the main text says 'we fit γ to Eq. 3', but the temperature-dependent linewidth model is Eq. (7), whereas Eq. (3) is the auxiliary function in the FCS model. This equation reference error makes the 2DES analysis difficult to follow. Additionally, Eq. (7) assumes a single thermally populated phonon mode; the manuscript should justify this assumption for (F-PEA)2PbI4 or discuss how a multi-mode model would affect αLO.","section":"Supplement, Eq. (7) and §III"}],"minor_comments":[{"comment":"There are numerous typographical errors that should be corrected in revision, including 'dipplacment' (Discussion), 'octehedra' (Introduction), 'it's' for 'its' (Abstract), 'Frank-Condon' for 'Franck-Condon', 'steady-sate' for 'steady-state', and 'Samll Structures' in Ref. [64]. These do not affect the scientific content but should be fixed.","section":"Throughout"},{"comment":"The manuscript describes a blue shift of the lowest-energy phonon from F-PEA to Cl-PEA, but Table I gives M1 = 3.8 meV for F-PEA and 4.0 meV for Cl-PEA. The shift is very small; please comment on whether this difference is within the spectral resolution and whether it is statistically significant.","section":"Fig. 4 and Table I"},{"comment":"The text states that the primary exciton peak and continuum edge appear at the same energy across all samples, but Fig. 1(d) appears to show subtle differences. Please clarify whether 'same energy' refers to within the experimental resolution and quantify any shifts.","section":"§II, Fig. 1(d)"},{"comment":"The Franck-Condon analysis of the XB resonance in panels (f)-(h) is used to explain the quenched/enhanced XB intensity, but the connection between the fitted S values and the observed absorption spectra is only qualitative. It would be helpful to overlay the calculated vibronic progressions on the measured absorption spectra to show the level of agreement.","section":"§III, Fig. 5"}],"recommendation":"major_revision","confidential_remarks":"The manuscript addresses an interesting and timely topic, and the experimental dataset is potentially valuable. My main concern is that the signature quantitative claims—the three-point correlation with bond angle variance and the two-material inverse correlation supporting polaronic protection—are not yet supported by the reported statistical and model-validation evidence. The authors should also make sure that all equation references and parameter definitions are consistent between the main text and the Supplement. I believe a major revision that adds uncertainty analysis, sensitivity tests for Eq. (4), and a clearer presentation of the dephasing comparison is within reach."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new piece is the first Huang-Rhys extraction for F-PEA and Cl-PEA and a new 2DES dephasing measurement for F-PEA. The paper uses an established FCS model to connect RISRS lineshapes to lattice displacement, and the data look carefully taken. But the headline correlation across three materials has no error bars on the fitted S values, and the dephasing comparison is effectively two materials, so the claims run ahead of the evidence.\n\nWhat is good: the measurements are demanding—RISRS on three derivatives at 15 K, plus 2DES—and the analysis follows the Turner/Arpin displaced-harmonic-oscillator framework rather than inventing a new one. The authors are honest that the lowest-energy mode deviates from the two-peak structure and that their crystallography does not explain the phonon blueshift. Extracting S for two new compounds and correlating with bond-angle variance is a useful step, and the 2DES data for F-PEA are new.\n\nSoft spots: no uncertainties on S. Equations 2–4 have several adjustable parameters; the two-PES expansion for M1 uses weights fixed by absorption rather than fit freely, and there is no uniqueness or residual analysis. With only three materials, the correlation is a tantalizing trend, not a demonstrated relation. The polaronic-protection claim rests on a two-point comparison between the new alpha_LO for F-PEA and a prior self-cited value for PEA; the alpha values come from different datasets and possibly different exciton features. That is not enough to establish an inverse correlation. These are not fatal flaws—the core idea is plausible and worth testing—but the manuscript should present uncertainties and ideally more cation variants before claiming a design rule.\n\nThis paper is for experimentalists working on 2D perovskites and exciton-phonon coupling, and for theorists who want to model polaron formation. I would send it to peer review because the new measurements deserve scrutiny, but I would push for major revision: error bars on S, a robustness check of the M1 fit, and a more cautious framing of the dephasing comparison. I would not yet cite the quantitative values as established.","headline":"Useful new Huang-Rhys data for two halogenated 2D perovskites, but the central trends rest on a three-point correlation without error bars and a two-material dephasing comparison, so the design-rule claims are premature.","tokens_in":18755,"tokens_out":3045,"would_cite":false,"duration_ms":31349,"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":"The organic cation in 2D perovskites tunes how strongly excitons are dressed by the lattice.","keywords":["2D hybrid perovskites","exciton polarons","Huang-Rhys parameter","RISRS","octahedral distortion","bond angle variance","thermal dephasing","2D electronic spectroscopy"],"falsifier":"Measure the Huang-Rhys parameter by an independent route—for example, the temperature dependence of the phonon-sideband intensities in photoluminescence, or non-resonant Raman cross-sections—on the same three crystals. If the ordering of $S$ is not F-PEA > PEA > Cl-PEA, or if the two-PES expansion gives materially different displacements when the relative $X_A$/$X_B$ weights are varied away from the absorption intensities, the claimed correlation between octahedral distortion and polaronic coupling collapses.","tokens_in":17683,"feed_emoji":"🧪","tokens_out":6371,"duration_ms":58570,"temperature":0.7,"pith_summary":"The paper aims to show that the organic cation in a 2D hybrid perovskite is a tunable handle on exciton polaron strength. In the model system (PEA)2PbI4 and its para-fluorinated and para-chlorinated derivatives, the authors use resonant impulsive stimulated Raman scattering (RISRS) to measure how far the lead-iodide lattice displaces when an exciton forms, quantified by the Huang-Rhys parameter $S = \\frac{1}{2}\\Delta^2$. They find a direct correlation between this displacement and static octahedral distortion (bond angle variance), with F-PEA showing the largest displacement and Cl-PEA the smallest, even though the average lattice and electronic structure stay nearly unchanged. They also report that F-PEA, the most strongly polaronic compound, shows the weakest thermal dephasing in 2D electronic spectroscopy, which they read as support for the idea that a self-induced lattice deformation cloud protects excitons from scattering. If correct, the work turns organic-cation substitution into a practical design lever for fine structure and coherence in 2D metal-halide perovskites.","feed_headline":"Strongest polaron dressing gives least dephasing in 2D perovskites","feed_subtitle":"Halogenating the organic cation tunes lattice flexibility, letting engineers suppress exciton scattering.","key_machinery":"The engine of the analysis is the displaced harmonic oscillator model of femtosecond coherence spectra (FCS), applied to probe-energy-resolved cuts of the RISRS beating maps. In this model the ground and excited exciton states are harmonic surfaces with the same frequency, offset by a displacement $\\Delta$, and the Huang-Rhys parameter $S = \\frac{1}{2}\\Delta^2$ measures how strongly the lattice dresses the exciton. Equations (2)–(3) express the Fourier-domain modulation amplitude as a weighted sum of Lorentzians broadened by a dephasing parameter $\\gamma$, with fit parameters $\\{S, \\gamma, \\omega_{eg}\\}$; because $\\gamma$ and the exciton energy are constrained by linear absorption, $S$ is effectively the only free parameter. For the lowest-energy mode, whose line shape deviates from the two-peak structure, the model is extended in Eq. (4) to a weighted sum of two excited-state potential energy surfaces ($X_A$ and $X_B$), with weights fixed by the absorption spectrum. The structural input is the bond angle variance $\\sigma^2$ computed from single-crystal XRD, used as the measure of octahedral distortion that is then correlated with the extracted displacements.","core_discovery":"The central claim is that exciton–phonon coupling in 2D metal-halide perovskites is set by the static flexibility of the inorganic octahedral framework, and that this flexibility can be tuned by halogen substitution on the organic cation without disturbing the excitonic landscape. Specifically, the authors claim that the equilibrium displacement $\\Delta$ of the excited-state potential energy surfaces along the two dominant lead-iodide phonon modes, extracted from RISRS spectra through the Huang-Rhys parameter $S = \\frac{1}{2}\\Delta^2$, follows the bond-angle-variance ordering F-PEA > PEA > Cl-PEA. They further claim that this ordering is mirrored inversely in thermal dephasing: the compound with the largest lattice displacement, F-PEA, has the smallest exciton–phonon scattering parameter $\\alpha_{\\mathrm{LO}}$ from temperature-dependent 2DES linewidths, consistent with polaronic protection. Taken together, the paper asserts that lattice engineering through the organic cation provides a pathway to control the polaronic character of excitons and thereby the optical fine structure and many-body scattering in these materials.","pith_inferences":["If the bond-angle-variance correlation is causal rather than coincidental, then a cheap static XRD measurement could serve as a screening metric for polaronic strength in newly synthesized 2D halide perovskites, which the paper does not explicitly claim.","The inverse relation between Huang-Rhys displacement and the 2DES scattering parameter suggests a trade-off worth testing: the same lattice flexibility that maximizes dressing may also soften the phonons that scatter excitons, potentially unifying the two parameters through a single structural coordinate.","A natural extension would be to vary the halogen on the inorganic site as well (for example, bromide or chloride octahedra) while holding the cation fixed; the framework here predicts the polaron displacement should track the resulting octahedral distortion, a test the present data cannot perform.","The authors compare dephasing only for PEA and F-PEA; completing the trend with Cl-PEA would test whether polaronic protection scales monotonically with $S$."],"forward_implications":["Organic-cation halogenation becomes a design knob for the exciton fine structure: in F-PEA the higher-energy $X_B$ resonance is quenched, while in Cl-PEA it is enhanced, matching the Franck-Condon weights built from the measured $S$ values.","Stronger polaronic dressing is associated with weaker thermal dephasing, so samples with greater octahedral distortion should show narrower homogeneous linewidths at elevated temperatures.","The same RISRS plus FCS protocol can be used to rank polaronic coupling in other 2D halide perovskites without relying on hard-to-deconvolute Raman sidebands.","Device-oriented consequences follow: tailoring the cation may suppress many-body scattering pathways relevant for polariton lasing and transport applications."],"supporting_citations":[{"why":"Establishes that RISRS phonon coherences reveal the polaronic character of excitons in (PEA)2PbI4 and identifies the low-energy modes as lead-iodide octahedral vibrations; the present RISRS data are compared directly to it.","marker":"[61]"},{"why":"Supplies the femtosecond-coherence-spectra lineshape model (displaced harmonic oscillator with Huang-Rhys parameter) that Eqs. (2)-(3) implement.","marker":"[117]"},{"why":"Reports the 2DES measurement of thermal dephasing and the exciton-phonon scattering parameter for (PEA)2PbI4 that the F-PEA comparison extends.","marker":"[86]"},{"why":"Articulates the exciton-polaron framework and the polaronic-protection hypothesis that the dephasing results are claimed to support.","marker":"[59]"},{"why":"Provides the single-crystal X-ray diffraction data from which the octahedral distortion parameters for the three compounds are estimated.","marker":"[101]"},{"why":"Links octahedral distortion to electronic structure and lattice flexibility, motivating the use of bond angle variance as a structural predictor.","marker":"[75]"}],"fun_headline_variants":["Lattice engineering tunes polaron protection in 2D perovskites","Polaron strength inverts dephasing in 2D perovskites","Halogen substitution tunes exciton-phonon coupling","Strongest polaron dressing gives least dephasing"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the displaced-harmonic-oscillator lineshape fits return a unique and faithful measure of the exciton-lattice displacement from the measured RISRS spectra, with no independent check on anharmonicity, mode mixing, or fit uniqueness.","fun_headline_variants_meta":{"raw":{"variants":["Lattice engineering tunes polaron protection in 2D perovskites","Polaron strength inverts dephasing in 2D perovskites","Halogen substitution tunes exciton-phonon coupling","Strongest polaron dressing gives least dephasing"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000825,"raw_usage":{"total_tokens":3695,"prompt_tokens":1123,"completion_tokens":2572,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":739,"completion_tokens_details":{"reasoning_tokens":2501}},"tokens_in":739,"tokens_out":2572,"duration_ms":19546,"temperature":1.0,"reasoning_tokens":2501,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T04:47:17.494539+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the Huang-Rhys parameter by an independent route—for example, the temperature dependence of the phonon-sideband intensities in photoluminescence, or non-resonant Raman cross-sections—on the same three crystals. If the ordering of $S$ is not F-PEA > PEA > Cl-PEA, or if the two-PES expansion gives materially different displacements when the relative $X_A$/$X_B$ weights are varied away from the absorption intensities, the claimed correlation between octahedral distortion and polaronic coupling collapses.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes that RISRS phonon coherences reveal the polaronic character of excitons in (PEA)2PbI4 and identifies the low-energy modes as lead-iodide octahedral vibrations; the present RISRS data are compared directly to it."},{"cited_title":"Principles of nonlinear optical spectroscopy (1995)","cited_arxiv_id":null,"evidence_quote":"Supplies the femtosecond-coherence-spectra lineshape model (displaced harmonic oscillator with Huang-Rhys parameter) that Eqs. (2)-(3) implement."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the 2DES measurement of thermal dephasing and the exciton-phonon scattering parameter for (PEA)2PbI4 that the F-PEA comparison extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Articulates the exciton-polaron framework and the polaronic-protection hypothesis that the dephasing results are claimed to support."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the single-crystal X-ray diffraction data from which the octahedral distortion parameters for the three compounds are estimated."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Links octahedral distortion to electronic structure and lattice flexibility, motivating the use of bond angle variance as a structural predictor."}],"review_version":1}