{"id":"a8e227ca-2b99-421b-b10b-cb2010c87737","arxiv_id":"1908.10130","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"Displaced harmonic oscillators with different curvatures produce a predicted extra substructure in vibronic absorption spectra, quantified by a new closed-form Franck-Condon expression; anharmonic and dissipative simulations for H2 and stiff-stilbene illustrate its effects.","lead":"This paper shows that when the ground and excited molecular potentials have very different curvatures, absorption spectra develop an extra repeating substructure, called the s-progression, on top of the usual vibronic peaks. It derives the Franck-Condon formula behind this effect and tests how anharmonicity and a dissipative environment reshape spectra and population dynamics in models of H2 and stiff-stilbene.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified; Eq. (39) is internally consistent and the sub-structure prediction survives scrutiny, with the dissipator caveat already acknowledged by the authors.","rationale":"The reader's verdict is CONDITIONAL, based on the Markovian harmonic-ladder dissipator and the hand-tuned stiff-stilbene model. I agree that these are genuine limitations for the quantitative dissipative spectra and dynamics. However, the strongest claim singled out by the reader is the analytic Franck-Condon result of Eq. (39). I checked the derivation step by step and found no internal inconsistency: the completing-the-square and Hermite-shift manipulations are valid, the missing binomial coefficient in Eq. (B17) is restored in Eq. (B19), and the sign difference in beta only affects phases, not absolute squares. The reduction to the standard displaced harmonic result, to the zero-displacement result, and to the large-curvature Hermite limit all work. The proposed s-progression is therefore a mathematically genuine feature of differing-curvature harmonic models, not an artifact of the dissipator or of the stochastic simulation. The dissipator concern affects the Morse and stiff-stilbene applications, but the paper explicitly acknowledges this in Sec. III C. Since my review does not uncover a flaw in the central analytic claim, I do not change the reader's verdict; the conditional recommendation should remain, pending the reproducibility/validation issues already identified.","tokens_in":21589,"tokens_out":22355,"duration_ms":220049,"concrete_test":"Verify Eq. (39) by direct numerical quadrature of Eq. (38) for the Fig. 2b parameters (alpha_g = 10 alpha_e, alpha_e = 1, displacement chosen so the equal-frequency reference has D = 30). Evaluate n = 0..80 and require (i) agreement with Eq. (39) to at least 1e-8, (ii) sum_n FC(n) = 1 to numerical precision, and (iii) first s-progression maxima appearing near the n values stated in Sec. III A (approximately n = 38 and n = 40).","verdict_should_be":"UNCHANGED","load_bearing_attack":"I re-examined the central analytic claim, namely that curvature difference produces an additional s-progression in the Franck-Condon envelope, as expressed in Eq. (39). The derivation in Appendix B checks out: the Gaussian completing-the-square step (B11), the Hermite shift identity (B17), and the treatment of even/odd terms are all consistent once the binomial coefficient is restored in (B19). The sign of beta in (B14) differs from the standard identity, but this only changes the overall phase of each overlap, not the squared Franck-Condon coefficient, so Eq. (39) is unaffected. The limiting cases are also correct: equal curvature reduces Eq. (39) to e^{-D}D^n/n!, zero displacement reduces to Eq. (35), and the large-curvature limit gives the Hermite-polynomial form of Eq. (42). The physical interpretation in Fig. 3, based on a narrow ground-state wavefunction sampling individual oscillations of high-n excited wavefunctions, is consistent with the mathematics. The one substantive limitation is the Markovian harmonic-ladder Lindblad dissipator used for the Morse and stiff-stilbene simulations; this is a real concern for the quantitative dissipative predictions, but the authors themselves flag it in Sec. III C, and it does not bear on the validity of the derived Franck-Condon coefficients or the existence of the s-progression in the isolated-molecule model.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper investigates linear absorption spectra and wavepacket dynamics for molecular vibrations beyond the standard equal-curvature displaced harmonic oscillator model. In Sec. III A the authors study a harmonic model with different ground- and excited-state curvatures, derive an analytic Franck-Condon expression (Eq. (39) with the derivation in Appendix B), and identify an additional 's-progression' caused by the narrow ground-state wavefunction sampling individual oscillations of high-lying excited vibrational wavefunctions. In Sec. III B they use a Morse potential with parameters for H2 and compare harmonic versus Morse raising/lowering operators within a Markovian stochastic Schrödinger equation. In Sec. III C they construct a model stiff-stilbene potential with barriers and curvature difference, and simulate absorption spectra and cis/trans population dynamics. The central analytic derivation is internally consistent and reduces correctly to Eq. (26) in the equal-curvature limit and to Eq. (35) in the zero-displacement limit.","tokens_in":21888,"tokens_out":5920,"duration_ms":64772,"significance":"If correct, Eq. (39) provides a quantitative, closed-form account of curvature-induced vibronic substructure and suggests an experimental route to estimate curvature differences from the width of the s-progression. The paper also gives a useful demonstration that the choice of Lindblad ladder operators can matter for anharmonic systems at large Huang-Rhys factors. Strengths are the complete analytic derivation in Appendix B, the explicit limiting-case checks, and the physical interpretation in Fig. 3. The main caveats are that the stochastic simulation results lack reported statistical uncertainty, the stiff-stilbene parameters are tuned to produce targeted dynamics, and the dissipator model is Markovian with harmonic ladder operators applied to anharmonic systems.","major_comments":[{"comment":"The stochastic Schrödinger equation results are presented without any measure of statistical uncertainty. No trajectory count, no error bars, and no convergence test are reported for the absorption spectra in Figs. 4, 5, and 7 or for the population dynamics in Fig. 8. Several conclusions are drawn from small amplitude differences, such as the statements that 'the ZPL has a larger amplitude' and that 'the peaks and troughs of the spectra are increased'; without sampling error the reader cannot distinguish these features from Monte Carlo noise. Please add trajectory counts, standard errors, or stated convergence criteria for all SSE averages.","section":"Sec. III B, Sec. III C, Figs. 4-9"},{"comment":"The stiff-stilbene calculations are not fully predictive because the central dynamical parameters are chosen to produce the target behavior: omega_e is set to 'ensure the desired period of 400 fs', the moment of inertia I is then computed from Eq. (45) using this choice, and gamma is chosen to 'ensure appropriate broadening in absorption spectra and significant population trapping in the cis-S1 state at 400 fs'. The later interpretation of the 400 fs population dynamics and the cis/trans photoselectivity therefore partly reports the input assumptions. Please separate fitted quantities from predictions, for example by showing how the spectra and dynamics respond to reasonable variations of omega_e, I, and gamma.","section":"Sec. III C, Eq. (45), Fig. 1"},{"comment":"The stiff-stilbene dissipative dynamics and spectra rest on harmonic lowering operators as Lindblad operators even though the system is anharmonic. The authors correctly flag in Sec. III C that the asymmetric broadening may differ if system-specific operators are used, and for the Morse model they test this in Fig. 7; however, no analogous sensitivity test is presented for stiff-stilbene. Since the broadened spectra and trapping dynamics are central to that section, please add a comparison with an anharmonic (or otherwise different) dissipator, or justify why the harmonic-ladder choice is sufficient for the stiff-stilbene model.","section":"Sec. III C, Sec. II B (Eqs. 13-17)"}],"minor_comments":[{"comment":"The Hermite addition identity in Eq. (B17) is missing the binomial coefficient; the correct coefficient appears later in Eq. (B19), so this appears to be a typo rather than a substantive error.","section":"Appendix B, Eq. (B17)"},{"comment":"Eq. (35) is valid for even n, and the factorial notation (n/2)! makes this implicit; please state the even-n condition explicitly in the main text.","section":"Sec. III A, Eq. (35)"},{"comment":"The axes of Fig. 2 are unlabeled; please add axis labels and a clear legend for the standard harmonic model and the differing-curvature model.","section":"Sec. III A, Fig. 2"},{"comment":"The Huang-Rhys parameter D is used for the Morse oscillator without an explicit definition of how the displacement is chosen for a given D; please clarify the relation between D and the Morse potential parameters.","section":"Sec. III B"}],"recommendation":"major_revision","confidential_remarks":"The central analytic result is sound and the s-progression interpretation is a meaningful contribution, although Eq. (39) is a special case of the earlier Chang formula. The main weakness is the numerical application: missing statistical characterization of the SSE results and the hand-tuned stiff-stilbene parameters. I would encourage the editor to require the authors to report trajectory counts/error bars and to reframe the stiff-stilbene section as a model demonstration rather than a predictive calculation, or to add sensitivity analyses. No concerns about citation pattern or novelty disclosure."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper's central analytic result is the closed-form Franck-Condon expression for displaced harmonic oscillators with different curvatures, Eq. (39), and the identification of the s-progression substructure that appears when the excited state is much shallower than the ground state. That part is genuinely good. I checked the Appendix B derivation and it is complete: the completing-the-square step, the Hermite shift identity, and the even/odd bookkeeping all work, and the expression reduces to the standard displaced harmonic result for equal curvature, to the un-displaced expression for zero displacement, and to a Hermite-polynomial form for large curvature difference. The physical explanation via a narrow ground wavefunction sampling individual oscillations of high-n excited wavefunctions is convincing. Credit where due: this goes beyond the existing literature, including the Chang (2005) formula they cite, because they connect the displaced and un-displaced cases and make the substructure interpretable.\n\nThe Morse and stiff-stilbene applications are more mixed. The Morse results show that harmonic raising/lowering operators overestimate asymmetric broadening for large Huang-Rhys factors, and that the difference is small for D=1.0, which is useful. But the stochastic simulations are presented without trajectory counts or error bars, so I cannot tell whether the small differences between the red and blue curves in Fig 7 are converged or noise. That is a real weakness in the numerical part.\n\nThe stiff-stilbene section is the softest. The model PES is hand-tuned to reproduce a desired 200 fs barrierless isomerisation time and the dissipation rate is chosen to give appropriate broadening and population trapping. That is fine for a proof-of-principle, but the resulting absorption spectra and dynamics should be described as illustrative, not as quantitative predictions. The authors do acknowledge the main caveats: they note that harmonic ladder operators may change the asymmetric broadening in anharmonic systems and that a continuum of states may contribute above S1(0). The reader's weakest assumption is exactly this, and it is fair.\n\nCitation pattern looks fine, with relevant prior work including Chang, Fidler and Engel, and the experimental stiff-stilbene papers. The paper is somewhat long and the connection between Secs IIIB and IIIC is not always tight, but it is readable.\n\nWho should read this: people computing vibronic spectra of molecules with large curvature differences, and people building open-quantum-system models for photoswitches. I would send it to a serious referee. The analytic part deserves publication even if the numerical parts need more care. A referee should ask for error bars or convergence data for the SSE simulations and a clearer statement that the stiff-stilbene parameters are a model, not an ab initio surface.","headline":"A solid analytic result on Franck-Condon substructure for differing curvatures, wrapped in a numerical paper whose stiff-stilbene part is illustrative rather than quantitative.","tokens_in":22442,"tokens_out":1976,"would_cite":true,"duration_ms":19587,"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":"When ground and excited potential surfaces curve differently, absorption spectra gain an additional weak substructure—the s-progression—and a closed-form Franck-Condon expression predicts where these extra peaks appear and how they decay.","keywords":["absorption spectra","Franck-Condon coefficients","curvature difference","Morse oscillator","stochastic Schrödinger equation","vibronic progression","anharmonicity","stiff-stilbene photoswitch"],"falsifier":"Record a low-temperature, high-resolution absorption spectrum of a molecule with a strongly steeper ground potential than excited potential and a substantial displacement (stiff-stilbene is the paper's candidate); the predicted s-progression should appear as a weak, decaying set of peaks on the short-wavelength side of the main vibronic envelope, with spacings set by the excited-state vibrational frequency and maxima at energies where the excited wavefunction has a local extremum. Observing a smooth Poisson-like envelope with no such substructure, or substructure independent of the curvature ratio, would refute the mechanism; equivalently, recomputing the same spectra with a non-Markovian bath and finding that the substructure vanishes would show that the Markovian Lindblad assumption, not the curvature difference, is doing the work.","tokens_in":21348,"feed_emoji":"🔬","tokens_out":7079,"duration_ms":62647,"temperature":0.7,"pith_summary":"This paper argues that the usual displaced-harmonic-oscillator picture of absorption spectra is incomplete when the ground and excited potential energy surfaces have different curvatures. Using a stochastic Schrödinger equation to include a dissipative environment, the authors show that curvature mismatch produces an additional, decaying substructure—dubbed the s-progression—in the vibronic progression, on top of the familiar main envelope. They derive a closed-form expression for the Franck-Condon coefficients that quantifies this substructure and reduces to the standard Huang-Rhys Poisson distribution when curvatures are equal. The same ideas are applied to a Morse-oscillator model of H2 and to a model potential for the stiff-stilbene photoswitch, where curvature difference, anharmonicity, and dissipation together explain the observed narrowing, shift, and extra peaks in the absorption lineshape.","feed_headline":"Unequal potential curvatures add hidden peaks to spectra","feed_subtitle":"A new Franck-Condon formula locates the extra s-progression and shows how to read curvature ratios from absorption data.","key_machinery":"The central object is the Franck-Condon coefficient $|\\langle \\psi_g^{n=0}|\\psi_e^n\\rangle|^2$ for harmonic oscillators whose ground ($\\omega_g$) and excited ($\\omega_e$) frequencies differ, generalised from the equal-curvature Poisson expression. The load-bearing identity is the closed-form result of Eq. (39), built from an explicit Hermite-polynomial summation containing the factor $(1-\\alpha_e/\\alpha)^l$; this factor is zero in the equal-curvature limit (recovering the standard $e^{-D}D^n/n!$ distribution) and is responsible for the s-progression when curvature differs. Mechanistically, the substructure appears because the ground-state wavefunction in the steeper potential is narrow enough to overlap individual oscillations of the excited vibrational wavefunction, so the Franck-Condon amplitude is set by whether the wavefunction at the Franck-Condon centre sits at a node (vanishing overlap) or a local extremum (large overlap). The paper also uses the stochastic Schrödinger equation with Lindblad operators (Eqs. 13–17) to include the environment, comparing harmonic ladder operators $L=a$ with Morse raising and lowering operators whose $L^\\dagger L|n\\rangle = \\Gamma(n-n/\\nu)|n\\rangle$, and the dipole correlation function $C_{\\mu\\mu}(t)$, whose Fourier transform is the absorption lineshape, to connect wavepacket dynamics to spectra.","core_discovery":"In a molecule with different ground- and excited-state potential curvatures and a sizable displacement, the absorption spectrum no longer follows the smooth Gaussian (Poisson) envelope predicted by identical-curvature displaced oscillators. Instead, the main vibronic progression narrows and shifts to larger vibrational quantum numbers, and a second, decaying series of Franck-Condon peaks appears at higher energies—the s-progression. The paper shows that this substructure arises because the narrow ground-state wavefunction (steeper ground potential) samples only isolated oscillations of the excited-state vibrational wavefunctions, so the overlap alternates between constructive and destructive depending on whether the excited wavefunction has a local extremum or a node at the excitation point. This mechanism is quantified by the derived Franck-Condon expression (Eq. 39), which in the equal-curvature limit recovers the standard Huang-Rhys distribution and in the zero-displacement limit yields the even-n-only progression of the un-displaced differing-curvature model. Applied to an anharmonic Morse oscillator, the paper finds that dissipation broadens the peaks asymmetrically in a way that deviates from harmonic predictions, and that the choice of harmonic versus Morse raising and lowering operators as Lindblad dissipators matters at large displacements; applied to stiff-stilbene, the curvature-mismatch and barrier features explain the narrow, red-shifted absorption band and the s-progression near 350 nm.","pith_inferences":["The s-progression could serve as a spectroscopic ruler for the curvature ratio independent of absolute displacement, since its onset is set by where the narrow ground wavefunction first overlaps individual excited-state oscillations.","A complementary substructure should appear in fluorescence emission spectra, with the roles of ground and excited curvatures swapped, offering an independent check of the mechanism.","The substructure's visibility depends on the Markovian dissipator assumption; experiments in solvents with different spectral densities, or non-Markovian simulations, could test how robust the predicted peaks are.","Analogous curvature-induced substructures may modulate cross-peak intensities in two-dimensional electronic spectra, providing a nonlinear-spectroscopy test of the same Franck-Condon mechanism."],"forward_implications":["Absorption spectra of molecules with strongly different ground and excited curvatures should show a measurable weak substructure (the s-progression) on the short-wavelength side of the main vibronic envelope; measuring its spacing and decay can estimate the curvature ratio.","Fitting spectra with equal-curvature displaced harmonic oscillators overestimates the Huang-Rhys parameter and misassigns peak widths and shifts for molecules like stiff-stilbene.","For anharmonic potentials at large displacements, using harmonic Lindblad operators overestimates high-frequency broadening; system-specific Morse operators change the lineshape qualitatively at strong dissipation.","In the stiff-stilbene model, damping combined with excited-state barriers traps population in the cis conformation on the excited state and increases the cis yield at 400 fs, suggesting environment tuning can control photoselectivity.","The derived Franck-Condon expression generalises the textbook Huang-Rhys formula, reducing to the Poisson distribution when curvatures are equal and to the even-n-only progression when displacement vanishes."],"supporting_citations":[{"why":"Supplies the stochastic Schrödinger equation and Lindblad master equation used for all dissipative wavepacket simulations.","marker":"[22]"},{"why":"Established that curvature difference shifts the absorption peak maximum and width, the starting point this paper extends to large curvature differences and displacement.","marker":"[47]"},{"why":"Derived the general Franck-Condon factors for differing-curvature displaced oscillators; the paper's Eq. (39) is a connected, more interpretable form.","marker":"[60]"},{"why":"Provides TD-DFT energies (e.g., S1(0) = 3.5 eV) used to build the stiff-stilbene model potential.","marker":"[50]"},{"why":"Supplies the schematic excited-state barriers and the 200 fs barrierless isomerisation time used in the stiff-stilbene dynamics.","marker":"[51]"},{"why":"Supplies the H2 Morse parameters (De, omega-x, beta) used in the anharmonic-oscillator calculations.","marker":"[61]"},{"why":"Gives the Morse raising and lowering operators used to test the dependence of spectra on the Lindblad-dissipator choice.","marker":"[64]"}],"fun_headline_variants":["Curvature mismatch adds extra peaks to absorption","Unequal curvatures generate hidden spectral substructure","New Franck-Condon rule for curved potentials","Beyond identical curvatures: spectra gain s-progression","How curvature mismatch alters vibronic spectra"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The environment is treated as a Markovian (memoryless) bath whose dissipative action is captured by harmonic raising and lowering operators; if realistic condensed-phase systems require a non-Markovian bath or system-specific dissipators, the predicted s-progression, asymmetric broadening, and stiff-stilbene dynamics would change.","fun_headline_variants_meta":{"raw":{"variants":["Curvature mismatch adds extra peaks to absorption","Unequal curvatures generate hidden spectral substructure","New Franck-Condon rule for curved potentials","Beyond identical curvatures: spectra gain s-progression","How curvature mismatch alters vibronic spectra"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000603,"raw_usage":{"total_tokens":2868,"prompt_tokens":1049,"completion_tokens":1819,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":665,"completion_tokens_details":{"reasoning_tokens":1749}},"tokens_in":665,"tokens_out":1819,"duration_ms":13023,"temperature":1.0,"reasoning_tokens":1749,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:51:55.665343+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Record a low-temperature, high-resolution absorption spectrum of a molecule with a strongly steeper ground potential than excited potential and a substantial displacement (stiff-stilbene is the paper's candidate); the predicted s-progression should appear as a weak, decaying set of peaks on the short-wavelength side of the main vibronic envelope, with spacings set by the excited-state vibrational frequency and maxima at energies where the excited wavefunction has a local extremum. Observing a smooth Poisson-like envelope with no such substructure, or substructure independent of the curvature ratio, would refute the mechanism; equivalently, recomputing the same spectra with a non-Markovian bath and finding that the substructure vanishes would show that the Markovian Lindblad assumption, not the curvature difference, is doing the work.","supporting_citations":[{"cited_title":"Improta \\ and\\ author F","cited_arxiv_id":null,"evidence_quote":"Provides TD-DFT energies (e.g., S1(0) = 3.5 eV) used to build the stiff-stilbene model potential."},{"cited_title":"Quick , author F","cited_arxiv_id":null,"evidence_quote":"Supplies the schematic excited-state barriers and the 200 fs barrierless isomerisation time used in the stiff-stilbene dynamics."},{"cited_title":"Micciarelli , author R","cited_arxiv_id":null,"evidence_quote":"Supplies the H2 Morse parameters (De, omega-x, beta) used in the anharmonic-oscillator calculations."},{"cited_title":"\\ Dong , author R","cited_arxiv_id":null,"evidence_quote":"Gives the Morse raising and lowering operators used to test the dependence of spectra on the Lindblad-dissipator choice."}],"review_version":1}