{"id":"d0a8434e-ac4e-4f27-b8f7-ac97c7dcbce1","arxiv_id":"2505.00532","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"In OCS, low-energy electron attachment proceeds through shape resonances followed by three mode-specific vibronic Feshbach resonances, which explains the observed kinetic energy bands of S- fragments and the shift between 5.5 and 6.0 eV electron beams.","lead":"This paper combines velocity-map imaging experiments with coupled-cluster and wavepacket calculations to show that low-energy electrons attaching to OCS molecules pass through shape resonances and then three mode-specific vibronic Feshbach resonances before dissociating into S- fragments. The work points to a general mechanism by which tuning electron energy can control fragmentation pathways in polyatomic molecules.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Bound-state EA-EOMCCSD curves along a 1D stretch cut cannot by themselves establish the three vibronic Feshbach resonances; no resonance-width or scattering calculation is presented.","rationale":"The reader's weakest assumption identifies exactly the load-bearing gap: the theoretical mechanism is inferred from bound-state EA-EOMCCSD potential curves along a one-dimensional coordinate, while the states of interest are metastable anions in the continuum. I agree and add two concrete aggravating details. First, the supplement itself states that the bending mode contributes substantially to high-energy S− formation, so the 1D stretch-only PECs cannot fully determine mode-specific vibronic resonances. Second, the claimed MCTDH-CAP validation is not accompanied by any Hamiltonian or computational parameters, making the main quantitative support uncheckable. The experimental dataset and angular-distribution analysis are valuable and appear carefully executed, but the causal chain from shape resonance through avoided crossings to three mode-specific Feshbach resonances and three KE bands is not uniquely established. A 2D wavepacket calculation with absorbing boundaries and a check of resonance widths would directly test whether the observed bands require the proposed mechanism. Since the paper's experimental core can stand while the mechanistic interpretation needs verification, the verdict remains CONDITIONAL rather than moving to ACCEPT or REJECT.","tokens_in":13314,"tokens_out":5006,"duration_ms":57445,"concrete_test":"Run a 2D (stretch+bend) MCTDH wavepacket propagation on the EA-EOMCCSD surfaces with an absorbing boundary, starting from the Franck-Condon region at 5.5 and 6.0 eV, and compare the computed S− kinetic-energy spectrum with Fig. 1(d). The three-band structure and the 5.5-to-6.0 eV intensity shift must appear only when the avoided-crossing couplings are active; additionally, compute CAP-based resonance widths along the stretch coordinate to confirm the resonances are narrow relative to vibrational spacings.","verdict_should_be":"UNCHANGED","load_bearing_attack":"To support the central claim, the three kinetic-energy (KE) bands of S− (Fig. 1d) must be caused by three vibronic Feshbach resonances formed by non-adiabatic resonant tunneling. The evidence offered is the avoided crossings among EA-EOMCCSD potential curves (Fig. 2a) and a derived three-band decomposition of vibrational levels (Fig. 2b,c). This is the least secure link. EA-EOMCCSD is a bound-state method; temporary anion states above the neutral threshold are not true bound states, and no resonance widths, lifetimes, or scattering cross sections are computed. Whether the avoided crossings correspond to narrow resonances capable of redistributing vibrational energy before dissociation is therefore not established. The problem is aggravated by dimensionality: the PECs are 1D cuts along CO-S stretch, whereas the supplement (Sec. II) explicitly states that the OCS bending mode contributes substantially to high-energy S− formation, and Sec. I concedes the axial-recoil model fails for highly excited bending. A single stretch coordinate cannot determine mode-specific vibronic resonances. Finally, the only quantitative theory-experiment comparison claimed to validate the dynamics, the MCTDH quantum flux in Fig. 1(a), is presented without Hamiltonian, basis, CAP parameters, or initial-state details, so it cannot be independently checked. If the avoided crossings are artifacts of the bound-state/basis treatment or of the 1D cut, the assignment of the three observed bands to specific vibronic Feshbach resonances collapses.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a kinematically complete study of dissociative electron attachment (DEA) in OCS using velocity map imaging at electron beam energies of 5.5 and 6.0 eV. The experiment resolves three kinetic energy bands of S− fragments and shows that the population shifts toward the highest band as the beam energy increases. The authors combine these measurements with EA-EOMCCSD potential energy curves along the CO–S stretch coordinate and an MCTDH wavepacket flux calculation to propose that, beyond the Franck–Condon region, non-adiabatic resonant tunneling through avoided crossings dynamically generates three mode-specific vibronic Feshbach resonances, which produce the three observed kinetic energy bands. The paper also presents angular distribution fits and a vibrational analysis of the transient anion with revival times.","tokens_in":13596,"tokens_out":4195,"duration_ms":44895,"significance":"If substantiated, the proposed mechanism would be an important advance in understanding resonance-mediated dissociation in polyatomic molecules, connecting shape resonances, vibronic Feshbach resonances, and mode-specific fragmentation. The experimental methodology is careful and current, particularly the use of conical wedge slicing to reduce artifacts in the momentum reconstruction, and the dataset is a valuable contribution to DEA studies of OCS. The theoretical tools (EA-EOMCCSD and MCTDH) are state of the art, and the paper explicitly reports the revival and super-revival times of the anion. However, the central mechanistic claim is presently supported only by qualitative comparisons, and the calculations do not yet establish that the three vibrational bands in the theory correspond to the three fragment kinetic energy bands.","major_comments":[{"comment":"The central claim of three vibronic Feshbach resonances is inferred from avoided crossings among EA-EOMCCSD potential energy curves, but EA-EOMCCSD is a bound-state method and the relevant anion states lie in the continuum. No complex absorbing potential, complex scaling, or scattering calculation is presented, so the resonance positions and, more importantly, the resonance widths and lifetimes are not established. Without widths, the assertion that the avoided crossings correspond to narrow Feshbach resonances capable of redistributing vibrational energy before dissociation remains speculative.","section":"Fig. 2(a) and main text discussion of PECs"},{"comment":"The three-band decomposition shown in Fig. 2(b,c) is obtained by partitioning the vibrational eigenvalues of the transient anion (listed in Supplementary Table I), not from a quantum dynamics calculation of the fragmentation products. The only dynamics shown, the MCTDH quantum flux in Fig. 1(a), is a total excitation function and is not decomposed into kinetic energy bands or final-state modes. Therefore the connection between the theoretical three vibrational bands and the experimentally observed three S− kinetic energy bands is not quantitatively demonstrated; the assignment is a qualitative hypothesis.","section":"Fig. 2(b,c) and Supplementary Sec. IV"},{"comment":"The theoretical model is strictly one-dimensional along the CO–S stretch, yet the paper claims mode-specific Feshbach resonances and mode-selective fragmentation. The supplement itself states that the bending mode contributes substantially to high-energy S− formation and that the axial recoil approximation fails for highly excited bending. A single stretch coordinate cannot determine three mode-specific resonances, and the manuscript provides no justification for why the stretch coordinate alone captures the relevant non-adiabatic dynamics.","section":"Supplementary Sec. II and main text p. 3"},{"comment":"The quantitative theory–experiment comparison in Fig. 1(a) is not independently checkable because the MCTDH calculation is described only as 'Normalized Quantum Flux' with no Hamiltonian, basis set, number of degrees of freedom, potential energy surfaces, CAP parameters, initial wavepacket, or propagation time. Since the flux is normalized and the experimental yield is also normalized, the agreement may reflect overall line-shape similarity rather than a parameter-free prediction. This weakens the statement that the theoretical results probe the experimental observations quantitatively.","section":"Fig. 1(a) and MCTDH description"}],"minor_comments":[{"comment":"The acronym 'MCDTH' in the caption should be corrected to 'MCTDH'.","section":"Fig. 1(a) caption"},{"comment":"The sentence 'This results rationales with our experimental observation' is ungrammatical; it should read 'This result rationalizes our experimental observation.'","section":"Last paragraph before references"},{"comment":"The text says 'FIG. 1(a) shows a minimum at 6.5 eV in the 2Σ symmetric excited state of anionic OCS,' but Fig. 1(a) shows the ion yield and quantum flux; the potential energy curve minimum is shown in Fig. 2(a). Please correct the cross-reference.","section":"Main text discussion of Fig. 1(a)"},{"comment":"The caption states '5.5 eV beam energy (a, b, c) and 5.5 eV beam energy (d, e, f)'; the second energy should presumably be 6.0 eV, consistent with the main text.","section":"Supplementary Fig. S1 caption"},{"comment":"The phrase 'Our theoretical results probe most of the experimental observations quantitatively and qualitatively' is vague; 'reproduce' or 'capture' would be more specific, and the extent of the quantitative agreement should be stated explicitly.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"The experimental data are solid and likely worth publishing, but the theoretical interpretation currently overreaches the calculations presented. A revision that adds resonance characterization (e.g., CAP-EA-EOMCCSD or a scattering calculation), a multi-dimensional treatment including bending, and a state- or energy-resolved wavepacket analysis would substantially strengthen the paper. Without these, the title and abstract claim of 'observation' of specific Feshbach resonances is not fully supported."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: the experiment is careful and the beam-energy effect is genuinely new, but the paper's central mechanistic claim—that three vibronic Feshbach resonances, formed by resonant tunneling, produce the three kinetic-energy bands—is plausible rather than demonstrated. The gap between the observed bands and the computed avoided crossings is the weak spot.\n\nWhat's actually good: the VMI measurements use conical wedge slicing, which avoids the low-momentum artifacts of the parallel-slicing method, and the calibration against the well-characterized O−/O2 resonance is standard. The documented shift of S− population from lower to higher kinetic-energy bands between 5.5 and 6.0 eV is a solid, new experimental observation, and the symmetry assignment from angular distributions is a reasonable addition. The EA-EOMCCSD curves are state-of-the-art for the bound-state problem, and the paper does not oversell the electronic-structure calculation itself.\n\nSoft spots: the main problem is that the three Feshbach resonances are inferred from avoided crossings along a one-dimensional CO–S stretch cut, using a bound-state method for states that actually lie in the continuum. No resonance widths or scattering cross sections are computed, so the assignment of the three KE bands to specific resonances is an interpretive leap. The problem is compounded by the authors' own supplementary material, which says the bending mode contributes substantially to high-energy S−; a single stretch coordinate cannot support a mode-specific vibrational assignment. The MCTDH flux in Fig. 1(a) is shown without enough detail (Hamiltonian, basis, CAP) to be checked independently, and the revival-time estimates from a cubic fit to vibrational levels are peripheral overreach. These are real limitations, but they are fixable and the experimental observations survive.\n\nWho it's for: specialists in dissociative electron attachment and electron-molecule scattering, and anyone interested in resonance-mediated fragmentation. The paper is worth reading for the experiment, and the theoretical hypothesis is worth having on record.\n\nRecommendation: this deserves a serious referee, not a desk rejection. A referee should push for either a scattering calculation (even a simple local-complex-potential one) that produces resonance widths and positions, or a clear statement that the Feshbach assignment is tentative. With that addition, the paper would be much stronger.","headline":"Solid experiment, plausible but under-supported mechanism—deserves referee time with a push for scattering calculations.","tokens_in":14138,"tokens_out":3415,"would_cite":true,"duration_ms":33886,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["34.80.Ht"],"model":"deepseek-v4-flash","headline":"Dissociative electron attachment in carbonyl sulfide is shown to run from shape resonances through non-adiabatic resonant tunneling at avoided crossings into three mode-specific vibronic Feshbach resonances, which fragment into three…","keywords":["dissociative electron attachment","OCS","vibronic Feshbach resonance","shape resonance","resonant tunneling","velocity map imaging","EA-EOMCCSD","mode-specific fragmentation"],"falsifier":"A scattering-level calculation (for example, R-matrix or complex-absorbing-potential electron-molecule scattering) that locates the temporary anion resonance poles along the same coordinates would settle the claim: if no Feshbach resonance appears near the three avoided crossings of the EA-EOMCCSD curves, or if the calculated resonance widths are too broad to yield the sharp kinetic-energy bands, the assignment fails.","tokens_in":13091,"feed_emoji":"⚛️","tokens_out":11206,"duration_ms":97902,"temperature":0.7,"pith_summary":"This paper reports kinematically complete measurements and matching theoretical calculations of dissociative electron attachment (DEA) in linear OCS. The authors argue that the simple picture of a shape resonance decaying near the Franck-Condon region is incomplete: beyond that region, non-adiabatic resonant tunneling through avoided crossings dynamically generates three mode-specific vibronic Feshbach resonances, each of which fragments into a distinct kinetic-energy band of S−. The interpretation explains both the three-band kinetic energy distribution of S− and why raising the electron beam energy from 5.5 to 6.0 eV shifts the S− population into higher-energy bending vibrational states. If correct, the work establishes a concrete polyatomic benchmark for how resonance-mediated dissociation channels are selected by electron energy and molecular geometry.","feed_headline":"Tunneling behind three S- bands in OCS electron capture","feed_subtitle":"Coupled-cluster curves show how shape resonances tunnel into vibronic Feshbach states as beam energy rises.","key_machinery":"The central objects are potential energy curves for OCS and OCS− calculated with the equation-of-motion coupled-cluster singles and doubles method for electron attachment (EA-EOMCCSD). Near the Franck-Condon region these curves show shape resonances; at larger CO-S separation they show avoided crossings, which the paper treats as the tunneling barriers through which the transient anion passes into three mode-specific vibronic Feshbach resonances. The companion Multi-Configurational Time-Dependent Hartree (MCTDH) wavepacket propagation supplies a normalized quantum flux that is compared with the experimental ion yield, and a partial-wave angular-distribution formula assigns the resonance symmetries.","core_discovery":"Using velocity-map imaging with conical time-gated wedge slicing, the paper measures S− fragments from DEA of OCS at beam energies around 5.5 and 6.0 eV and observes three characteristic kinetic energy bands, with the third band gaining population as the beam energy increases. Potential energy curves computed with EA-EOMCCSD show that in the Franck-Condon region the electron attaches through shape resonances ($^{2}\\Delta$, $^{2}\\Pi$, and four $^{2}\\Sigma$ states) that decay into broad vibronic features. Beyond the Franck-Condon region, the curves exhibit avoided crossings, and the paper identifies resonant tunneling through these crossings as the mechanism that forms three distinct vibronic mode-specific Feshbach resonances. These resonances redistribute vibrational energy among modes before the OCS− dissociates, producing the three kinetic energy bands and the energy-dependent shift of S− population between 5.5 and 6.0 eV.","pith_inferences":["If the avoided-crossing assignment is right, the same three-band kinetic-energy pattern could serve as a mode-specific fingerprint for vibronic Feshbach resonances in other linear triatomics, such as CO2 or CS2.","A direct test would be to compute true resonance widths and positions with a scattering method; if the bound-state EA-EOMCCSD curves happen to give the right crossings, they become a cheap screening tool for polyatomic DEA systems.","The predicted ~8 ps revival implies that a coherent vibrational wavepacket in the transient anion might be observable with pump-probe techniques, connecting DEA to quantum-control studies.","Measuring angular distributions continuously across 5.5-6.0 eV could separate the $\\Sigma+\\Pi$ and $\\Sigma+\\Delta$ contributions and test whether the two pathways follow different avoided crossings."],"forward_implications":["The three kinetic-energy bands of S− are direct experimental signatures of three mode-specific vibronic Feshbach resonances formed by resonant tunneling through avoided crossings.","Raising the beam energy from 5.5 to 6.0 eV activates dense high-energy bending states by vibronic intensity borrowing and overlapping resonances, explaining the abrupt increase in high-momentum S− fragments.","The linear-to-bent geometry change during dissociation shows up in angular distributions consistent with $\\Sigma \\to \\Sigma+\\Pi$ or $\\Sigma \\to \\Sigma+\\Delta$ transitions, linking symmetry to the fragmentation channel.","The strong anharmonicity of the transient OCS− anion predicts a vibrational revival near 8.2 ps, a time scale that could be probed in ultrafast experiments.","MCTDH quantum-flux calculations reproduce the experimental ion-yield energy dependence, providing a benchmark for resonance-mediated dissociation in polyatomic molecules."],"supporting_citations":[{"why":"Supplies the EA-EOMCCSD electronic-structure method that generates the transient anion potential energy curves used to identify avoided crossings.","marker":"[22]"},{"why":"Supplies the MCTDH wavepacket propagation method used to compute the normalized quantum flux and population dynamics.","marker":"[25]"},{"why":"Provides the MCTDH formalism used for time-resolved nuclear dynamics of the transient anion.","marker":"[26]"},{"why":"Supplies the partial-wave angular-distribution formula used to fit the S− angular distributions and assign resonance symmetries.","marker":"[28]"},{"why":"Provides the O−/O2 6.5 eV resonance used to calibrate electron beam energy and fragment kinetic energy.","marker":"[33]"},{"why":"Reports the earlier S−/OCS ion-yield resonances at 5.0 and 6.5 eV that this study reproduces and extends.","marker":"[35]"},{"why":"Supplies the conical wedge-slice imaging technique and earlier S−/OCS angular-distribution analysis that the present experiment builds on.","marker":"[15]"}],"fun_headline_variants":["Tunneling from shape to Feshbach states yields three S- bands from OCS","Resonant tunneling into vibronic Feshbach states drives OCS's three S- bands","Shape resonance tunneling explains OCS's mode-specific triple S- fragmentation","Tunneling through avoided crossings splits OCS's S- into three bands"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that potential energy curves from a bound-state coupled-cluster calculation, taken along a single CO-S stretch coordinate, faithfully represent the metastable OCS− states that lie in the continuum and the avoided crossings that supposedly create the three Feshbach resonances, even though the paper's own supplement says the bending mode contributes substantially to high-energy S− formation.","fun_headline_variants_meta":{"raw":{"variants":["Tunneling from shape to Feshbach states yields three S- bands from OCS","Resonant tunneling into vibronic Feshbach states drives OCS's three S- bands","Shape resonance tunneling explains OCS's mode-specific triple S- fragmentation","Tunneling through avoided crossings splits OCS's S- into three bands"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000859,"raw_usage":{"total_tokens":3744,"prompt_tokens":978,"completion_tokens":2766,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":594,"completion_tokens_details":{"reasoning_tokens":2678}},"tokens_in":594,"tokens_out":2766,"duration_ms":19602,"temperature":1.0,"reasoning_tokens":2678,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:39:03.523504+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A scattering-level calculation (for example, R-matrix or complex-absorbing-potential electron-molecule scattering) that locates the temporary anion resonance poles along the same coordinates would settle the claim: if no Feshbach resonance appears near the three avoided crossings of the EA-EOMCCSD curves, or if the calculated resonance widths are too broad to yield the sharp kinetic-energy bands, the assignment fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the EA-EOMCCSD electronic-structure method that generates the transient anion potential energy curves used to identify avoided crossings."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the MCTDH wavepacket propagation method used to compute the normalized quantum flux and population dynamics."},{"cited_title":"Meyer, U","cited_arxiv_id":null,"evidence_quote":"Provides the MCTDH formalism used for time-resolved nuclear dynamics of the transient anion."},{"cited_title":"Nandi, V","cited_arxiv_id":null,"evidence_quote":"Provides the O−/O2 6.5 eV resonance used to calibrate electron beam energy and fragment kinetic energy."},{"cited_title":"Iga and S","cited_arxiv_id":null,"evidence_quote":"Reports the earlier S−/OCS ion-yield resonances at 5.0 and 6.5 eV that this study reproduces and extends."}],"review_version":1}