{"id":"213e1013-3a60-4eb8-9fd5-8ae595f0b8af","arxiv_id":"2505.20823","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Converged quantum dynamics of the cis-PSB3 model reveals coherent oscillations in ground-state cis/trans populations that were missed by earlier, less converged calculations.","lead":"This paper revisits quantum dynamics simulations of a three-mode model of a retinal chromophore fragment and shows that a larger basis along the twist coordinate changes the picture. It reports coherent oscillations in ground-state populations that were missed by earlier simulations, and argues that surface-hopping trajectories actually agree well with the corrected quantum result.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Benchmark correction is credible, but the physical oscillation claim rests on the untested constant-metric kinetic operator for large-amplitude torsion.","rationale":"The reader's CONDITIONAL verdict is corroborated: the convergence evidence is strong, including multiple grid sizes, cross-code comparison with ElVibRot, MCTDH tests with up to 160 SPFs, and close agreement with TSH/TSH-EDC. I do not see an internal inconsistency in the benchmark claim that Ref.[13] used an unconverged Tors basis. The remaining load-bearing gap is the physical interpretation: the oscillations are obtained under the constant-diagonal-metric KEO, which is explicitly approximate and untested against the exact G-matrix for large-amplitude torsion. Since the paper itself traces the new dynamics to Tors kinetic energy, the metric approximation is not a peripheral modeling detail—it is the operator that defines that kinetic energy. A direct test with the exact KEO would settle whether the predicted coherence is a model artifact. This is essentially the weakest assumption the reader identified, though I emphasize the KEO component more sharply than the broader envir\\onmental and mode-omission concerns.","tokens_in":11330,"tokens_out":8129,"duration_ms":93244,"concrete_test":"Build the exact coordinate-dependent G-matrix kinetic operator for the BLA/Tors/HOOP model (e.g., with ElVibRot-TnumTana, which the authors already use) and rerun the converged grid calculation with the same 256-point Fourier basis and identical initial state. If the ~30 fs P_S0_cis steps and the ~100 fs P_S0_trans oscillation survive to within the grid-convergence tolerance, the concern is resolved; if they shift by more than the convergence error or disappear, the physical prediction is an artifact of the constant-metric KEO.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's new physical result—the ~30 fs step-like oscillations in P_S0_cis and the ~100 fs trans-to-cis back-reaction in P_S0_trans—depends on the kinetic energy operator used for the three curvilinear coordinates. Section 2 states: 'we only adopted the approximate kinetic energy operator of Ref.[13] with a constant and diagonal metric tensor.' For a dihedral torsion undergoing large-amplitude motion, which the authors invoke to explain the oscillations, the exact G-matrix is coordinate-dependent; constant 1/I_theta is an approximation inherited from Ref.[13]. The analysis in Section 3 attributes the new behavior to a larger Tors kinetic energy, but that kinetic energy is computed with this approximate operator. If the exact kinetic coupling between Tors, HOOP, and BLA redistributes energy differently, the oscillations and back-reaction could be an artifact of the constant-metric approximation rather than a robust property of cis-PSB3. The benchmark claim that Ref.[13]'s basis was unconverged is not threatened, because the same approximate operator was used there; but the title-level claim 'Quantum Dynamics Predicts Coherent Oscillatory Behavior' is. The acknowledged neglect of environment and additional modes further limits physical extrapolation, but the KEO approximation is the most immediate unquantified threat.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper revisits the quantum dynamics of the two-state, three-mode model of the cis-PSB3 photoisomerization originally studied in Ref. [13]. Using much larger primitive basis sets along the torsional coordinate (up to 256 Fourier grid points, versus 136 in Ref. [13]), as well as cross-checks with the ElVibRot code and MCTDH calculations with up to 160 SPFs, the authors find that the earlier quantum dynamics was unconverged. The converged calculations reveal two coherent oscillatory features in the ground-state populations: step-like oscillations in the cis population with a period of about 30 fs starting around 60 fs, and a trans-to-cis back-reaction with a period of about 100 fs. The authors further compare with trajectory-based methods and report that TSH and TSH-EDC are in close agreement with the converged QD results, whereas Ehrenfest and CT-MQC show only qualitative agreement. The paper highlights the challenge of converging QD for large-amplitude torsional motion.","tokens_in":11650,"tokens_out":7095,"duration_ms":68199,"significance":"If the results are correct, the paper provides an important correction to a widely used benchmark for nonadiabatic dynamics methods, demonstrating that the previous QD reference was unconverged and that surface-hopping methods perform better than previously assessed. The identification of coherent oscillations and a ground-state back-reaction is a new physical prediction of the model. The study also provides a valuable cautionary example of the large number of basis functions or SPFs needed to converge quantum dynamics for strongly anharmonic large-amplitude motion. Notable strengths are the extensive convergence testing across two independent codes (Quantics and ElVibRot), the systematic MCTDH convergence study, and the direct comparison with multiple trajectory-based methods.","major_comments":[{"comment":"The central physical claim of coherent oscillations and the trans-to-cis back-reaction relies on the kinetic energy operator, which the authors state is the approximate constant-diagonal-metric operator of Ref. [13]. For the large-amplitude torsional motion invoked to explain the oscillations, the exact G-matrix is coordinate-dependent; the constant metric is an untested approximation. Please test the sensitivity of the populations and kinetic energies to this approximation (e.g., by using a coordinate-dependent G-matrix or an exact kinetic operator for the three curvilinear coordinates), or at a minimum provide a quantitative argument for why the constant-metric operator is accurate in the region explored by the wavepacket. Without this, the title-level claim that 'Quantum Dynamics Predicts Coherent Oscillatory Behavior' may be a property of the approximate operator rather than of the model molecular system.","section":"Section 2, kinetic energy operator discussion"},{"comment":"The claim that the calculations are 'fully converged' is based on visual overlap of population curves, with no quantitative error measure. Given the non-monotonic approach to convergence shown in Fig. 3 (magenta, light-green, orange, dark-green for 94, 136, 162, 192 grid points), please provide quantitative measures such as integrated absolute differences between successive basis sizes (e.g., Nq=256 vs 384, or 192 vs 256) or norm-based wavepacket differences. Also, although the text states tests up to Nq=384 were performed, no 384-point curve is shown; please display it or report the difference from the 256-point result. This quantitative evidence is needed to support the central assertion that Ref. [13] was unconverged and that the new oscillations are converged features.","section":"Section 2, convergence tests"},{"comment":"The statement that TSH and TSH-EDC are in 'quantitative agreement' and that Ehrenfest/CT-MQC are not is based on visual inspection of the population curves and kinetic energy plots. Please provide quantitative error metrics (e.g., time-integrated absolute deviations from the QD reference for each population and for the quantum yield) to support the revised benchmark conclusions. This is particularly important because the paper's message includes a reassessment of the performance of trajectory-based methods.","section":"Section 3, comparison with quantum-classical methods"}],"minor_comments":[{"comment":"The light-green curves are labeled both as the Ref. [13] result (32-136-60 grid) and as one of the convergence series with Nq=136; this dual labeling is confusing and should be clarified in the caption or text.","section":"Section 2, Figure 3"},{"comment":"There are two typos: 'Tors ans HOOP' should be 'Tors and HOOP', and 'are non correctly reproduced' should be 'are not correctly reproduced'.","section":"Section 2, text"},{"comment":"The caption does not specify the color-to-grid mapping for the various curves; please add a legend or explicit description in the text.","section":"Section 2, Figure 2 caption"},{"comment":"The matching of the HOOP kinetic energy period with the cis population oscillation period is stated qualitatively; please quantify both periods and, if possible, perform a simple Fourier analysis to support the claim.","section":"Section 3, HOOP period analysis"},{"comment":"The cis/trans partition at Tors=±90 degrees is not tested; a brief statement about the sensitivity of the oscillation features to the dividing surface would be helpful.","section":"Section 3, cis/trans partition"},{"comment":"The MCTDH calculations use 80, 256, and 90 grid points for BLA, Tors, and HOOP, which differs from the primitive basis in Table 1 (30, 256, 60); please state the reason for the increased BLA and HOOP grids in MCTDH.","section":"Appendix A, MCTDH grid"},{"comment":"Reference [10] is cited as a preprint; please consider citing the published version if it has appeared.","section":"Introduction, reference [10]"}],"recommendation":"major_revision","confidential_remarks":"The paper is a useful benchmark correction, but the KEO approximation and the lack of quantitative convergence metrics need to be addressed. The authors should be encouraged to either test an exact kinetic operator or clearly restrict their claims to the approximate model. The manuscript fits the journal's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: the central claim holds—Ref. [13]'s quantum dynamics benchmark for the cis-PSB3 model was unconverged, and the converged calculation produces qualitatively new oscillatory S0 populations. The work is careful and the new result is genuine, not a refit.\n\nWhat is new: with 256 Fourier functions for Tors (vs 136) and checks across 94–384 grids plus MCTDH up to 160 SPFs, the S0 populations develop step-like cis oscillations with ~30 fs period starting at ~60 fs and a trans→cis back-reaction with ~100 fs period. These features were absent in Ref. [13]. The two-code cross-check (Quantics/ElVibRot) and single-set vs multi-set MCTDH agreement support the numerical conclusion. The correction to the trajectory-method assessment is also valuable: TSH and TSH-EDC track the converged QD closely, Ehrenfest/CT-MQC less so.\n\nSoft spots, in proportion. First, convergence is shown by visual overlap only; no quantitative error estimates, and no input files or raw populations are provided. For a paper whose message is 'the previous benchmark was wrong,' that is an unnecessary reproducibility gap. Second, and more substantive: the oscillations are attributed to large-amplitude Tors motion, but the KEO is the constant diagonal metric from Ref. [13]. For a torsion, the exact G-matrix is coordinate-dependent. The basis-size correction is unaffected because both runs used the same KEO, but the title-level claim that quantum dynamics 'predicts coherent oscillatory behavior' is only as strong as that approximation. The paper could have tested this, or at least flagged it quantitatively. The known model limitations (no environment, no extra modes) are acknowledged, but the KEO question is the one that matters for the physical interpretation.\n\nSerious thinker: yes—the paper is honest about what it does and does not show, and the numerics are extensive. I'd send it to review. It deserves a referee who can check the convergence claims and push for archived inputs and a KEO sensitivity test. I'd cite it if I worked on nonadiabatic benchmarks.","headline":"A credible benchmark correction with a real new oscillation signal, though the physical interpretation leans on an approximate kinetic operator.","tokens_in":12173,"tokens_out":2334,"would_cite":true,"duration_ms":24293,"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":"Converged quantum dynamics predicts two coherent oscillatory patterns in the early-time photoisomerization of a retinal chromophore model that earlier simulations missed.","keywords":["quantum dynamics","photoisomerization","cis-PSB3","retinal chromophore model","nonadiabatic dynamics","conical intersection","coherent oscillations","surface hopping"],"falsifier":"A decisive calculation is to propagate the same model with a coordinate-dependent kinetic energy operator: if the ~30 fs and ~100 fs oscillations vanish or move outside numerical error, the constant-metric approximation is their source rather than a robust dynamical feature.","tokens_in":11132,"feed_emoji":"🧪","tokens_out":10970,"duration_ms":103375,"temperature":0.7,"pith_summary":"This paper revisits the quantum dynamics of a two-electronic-state, three-vibrational-mode model of the retinal chromophore fragment cis-PSB3, a standard benchmark for trajectory-based methods in photoisomerization. The authors show that the earlier quantum-dynamics results for this model were not converged: with a much larger torsional basis (256 Fourier functions instead of 136), the ground-state populations develop two coherent oscillatory patterns that were absent before. A step-like rise in the ground-state cis population with about a 30 fs period appears after roughly 60 fs, and an $S_0$ trans-to-cis back-reaction produces a second oscillation with roughly a 100 fs period. Because the converged results closely match surface-hopping trajectories, the paper concludes that the earlier assessment of quantum-classical methods must be revised. The wider point is that even reduced-dimensionality photoisomerization models need very large bases to capture large-amplitude torsional motion.","feed_headline":"Quantum convergence reveals two oscillations in photoisomerization","feed_subtitle":"A larger torsional basis exposes ~30 fs cis steps and a ~100 fs trans-to-cis back-reaction missed by the prior benchmark.","key_machinery":"The argument is carried by a two-electronic-state, three-vibrational-mode analytical model of cis-PSB3 whose coordinates are the bond-length-alternation stretch (BLA), the reactive double-bond torsion (Tors), and the hydrogen-out-of-plane wag (HOOP). The numerical machinery is a converged primitive basis: 256 Fourier functions for Tors (up to 384 tested), 30 harmonic-oscillator functions for BLA, 60 for HOOP, with Short Iterative Lanczos propagation cross-checked against MCTDH runs using up to 160 single-particle functions on Tors in $S_0$ and 110 in $S_1$. The physical mechanism that carries the discovery is kinetic-energy redistribution: after population transfer through the $S_0/S_1$ conical intersection, the low-frequency Tors mode receives a large amount of kinetic energy, enabling the large-amplitude motion and the ground-state back-reaction, while the ~30 fs period of the $S_0$ cis steps matches the HOOP kinetic-energy oscillation.","core_discovery":"On the paper's own terms, the central discovery is that the quantum dynamics of the two-state ($S_0/S_1$), three-mode cis-PSB3 model was not converged in the earlier benchmark: enlarging the torsional Fourier basis from 136 to 256 points (with tests up to 384) changes the ground-state populations qualitatively, not just quantitatively. After convergence, the $S_0$ cis population rises in step-like increments with roughly a 30 fs period starting near 60 fs, and the $S_0$ trans population shows a maximum near 85 fs, a minimum near 125 fs, and renewed growth toward 195 fs. The latter is interpreted as a trans-to-cis back-reaction on the ground state driven by the large kinetic energy accumulated in the torsional mode, with a correlated loss in $S_0$ trans population of about 0.1 between 86 and 130 fs. The same converged dynamics puts the short-time quantum yield near 0.75 at about 75 fs and brings surface-hopping methods (TSH and TSH-EDC) into close agreement with the quantum reference, while Ehrenfest and CT-MQC agree only qualitatively.","pith_inferences":["This reading suggests that the need for very large torsional bases is likely generic for reduced-dimensionality photoisomerization models, because photoexcitation concentrates kinetic energy in the low-frequency large-amplitude mode.","A testable extension is to repeat the dynamics with a coordinate-dependent kinetic energy operator; if the oscillation periods or amplitudes change substantially, the predicted coherence is partly an artifact of the constant-metric approximation.","Should the oscillations survive coupling to additional modes and a dissipative environment, they could be observable in femtosecond pump-probe experiments as periodic modulations of ground-state recovery in retinal chromophore models.","The near-quantitative agreement between converged quantum dynamics and surface hopping raises the option of using TSH as a substitute benchmark when full quantum convergence is impractical, with a torsional convergence check standing in for the quantum basis check."],"forward_implications":["The earlier quantum-dynamics benchmark for this model is superseded; assessments of trajectory methods based on it must be revised, with TSH and TSH-EDC now matching the converged quantum reference closely.","The ground-state cis and trans populations of the model are non-monotonic, so the time-dependent quantum yield oscillates rather than rising to a plateau.","The HOOP mode actively shapes early-time relaxation: its kinetic-energy period coincides with the ~30 fs steps in the ground-state cis population.","MCTDH users need to converge the number of single-particle functions along strongly excited torsional coordinates; too few SPFs reproduce the same missing oscillations as the small primitive grid.","The model, once converged, provides a consistent reference for the photoisomerization quantum yield at short times, near 0.75 at about 75 fs."],"supporting_citations":[{"why":"Defines the two-state, three-mode diabatic potential energy surfaces and couplings for cis-PSB3 on which all dynamics is run.","marker":"[11]"},{"why":"Supplies the kinetic energy operator, the earlier unconverged quantum dynamics, and the trajectory-based results that the paper re-benchmarks.","marker":"[13]"},{"why":"Quantics code used for the converged Short Iterative Lanczos quantum-dynamics propagation.","marker":"[14]"},{"why":"ElVibRot code used as an independent check of the converged populations with matching basis sets.","marker":"[15]"},{"why":"Provide the MCTDH method whose single-particle-function convergence is tested in Appendix A.","marker":"[16,17]"},{"why":"Tully surface hopping, the trajectory method closest to the converged quantum dynamics.","marker":"[20]"},{"why":"Energy-based decoherence correction (TSH-EDC) whose results nearly coincide with TSH.","marker":"[21]"},{"why":"Defines CT-MQC, the coupled-trajectory method compared against the quantum reference.","marker":"[22]"}],"fun_headline_variants":["Converged quantum dynamics exposes cis-trans oscillations","Bigger basis flips photoisomerization picture: reveals steps","Quantum model: 30 fs steps, 100 fs back-reaction","Unconverged? No: torsional basis reveals coherent steps","Photoisomerization oscillations: cis steps, trans back-reaction"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the two-state, three-mode model, with its simplified kinetic-energy operator that treats the coordinates as independent with fixed inertial factors, faithfully represents the large-amplitude torsional and HOOP motion of cis-PSB3; if that reduced model is unrepresentative, the predicted oscillations hold only for the model.","fun_headline_variants_meta":{"raw":{"variants":["Converged quantum dynamics exposes cis-trans oscillations","Bigger basis flips photoisomerization picture: reveals steps","Quantum model: 30 fs steps, 100 fs back-reaction","Unconverged? No: torsional basis reveals coherent steps","Photoisomerization oscillations: cis steps, trans back-reaction"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000162,"raw_usage":{"total_tokens":1236,"prompt_tokens":938,"completion_tokens":298,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":554,"completion_tokens_details":{"reasoning_tokens":210}},"tokens_in":554,"tokens_out":298,"duration_ms":3396,"temperature":1.0,"reasoning_tokens":210,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:46:42.910917+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive calculation is to propagate the same model with a coordinate-dependent kinetic energy operator: if the ~30 fs and ~100 fs oscillations vanish or move outside numerical error, the constant-metric approximation is their source rather than a robust dynamical feature.","supporting_citations":[{"cited_title":"H.; Yang, X.; De Vico, L.; Olivucci, M","cited_arxiv_id":null,"evidence_quote":"Defines the two-state, three-mode diabatic potential energy surfaces and couplings for cis-PSB3 on which all dynamics is run."},{"cited_title":"Quantum and quantum-classical studies of the photoisomerization of a retinal chromophore model","cited_arxiv_id":null,"evidence_quote":"Supplies the kinetic energy operator, the earlier unconverged quantum dynamics, and the trajectory-based results that the paper re-benchmarks."},{"cited_title":"Quantics: A general purpose package for Quantum molecular dynamics simulations","cited_arxiv_id":null,"evidence_quote":"Quantics code used for the converged Short Iterative Lanczos quantum-dynamics propagation."},{"cited_title":"ElVibRot-TnumTana Quantum Dynamics Code","cited_arxiv_id":null,"evidence_quote":"ElVibRot code used as an independent check of the converged populations with matching basis sets."},{"cited_title":"Critical appraisal of the fewest switches algorithm for surface hopping","cited_arxiv_id":null,"evidence_quote":"Energy-based decoherence correction (TSH-EDC) whose results nearly coincide with TSH."},{"cited_title":"K.; Agostini, F.; Gross, E","cited_arxiv_id":null,"evidence_quote":"Defines CT-MQC, the coupled-trajectory method compared against the quantum reference."}],"review_version":1}