{"id":"72552b5a-823e-452a-b37e-2838fcc45651","arxiv_id":"2411.16137","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"New 5D quantum scattering rates for C2- vibrational quenching by para- and ortho-H2 are computed, but agreement with the 20 K experiment requires a rotationally restrictive, acknowledged nonphysical channel selection.","lead":"This paper calculates how hydrogen molecules collide with and cool down carbon dimer anions, using a new five-dimensional interaction surface. The results are meant to guide experiments on laser cooling molecular anions, but the computed cooling rates only match the single experimental measurement after selecting a specific subset of collision outcomes.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Agreement with the measured 20 K quenching rate rests on a post hoc j1'=0 filter that the authors themselves call nonphysical; the fully summed 5D rates are one to two orders of magnitude above experiment.","rationale":"The reader's weakest_assumption correctly identifies the most load-bearing concern: the j1'=0 final-state restriction is the only thing that brings the computed 5D rates close to the single experimental point, and the paper itself acknowledges this restriction is nonphysical and that the experiment did not observe final rotational states. This concern directly undermines the abstract's headline claim of agreement with experiment, because the measured rate is an inclusive vibrational quenching rate to which all final rotational channels contribute. The paper is nevertheless a substantial and transparent computational contribution: it provides a new 5D PES, detailed CC scattering calculations, convergence checks, and a clear statement of the limitation. The appropriate verdict is therefore CONDITIONAL, as the reader concluded; my stress-test does not change that verdict. I find no additional load-bearing concern that would move the paper to ACCEPT or REJECT. The main open question is whether some as-yet-unidentified physical mechanism (e.g., rapid secondary rotational relaxation before detection, or a detection bias toward cold rotors) could justify the j1'=0 projection; until such a mechanism is demonstrated, the agreement with experiment remains conditional.","tokens_in":27370,"tokens_out":3703,"duration_ms":42513,"concrete_test":"Recompute the comparison made in Figure 15 using the inclusive summed rate coefficient defined by Eq. (16) in Section IV.C: thermally average over the initial j1 states (e.g., j1=0,2,4 at 20 K) and sum over all final j1' and j2' states, exactly as done for the full red/purple curves, and compare that summed value at 20 K with the experimental point from Ref. [52]. If the ratio k_full(20 K)/k_j1'=0(20 K) is greater than 10, then the reported agreement is not robust and the central claim should be restated as a conditional match under an unjustified channel restriction. A complementary decisive check is to determine from the experimental detection scheme used in Ref. [52] whether it can distinguish final rotational states of C2-; if it cannot, the inclusive comparison is the only valid one.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's only point of agreement with experiment is the dashed red curve in Figure 15, constructed in Section IV.C by restricting the vibrationally quenched C2- to the final rotational state j1'=0. The authors explicitly state in Section IV.B that the experiments of Ref. [52] 'were not observing which rotational states will be populated after the anion's decay to its ν = 0 level,' and in Section V they call the j1'=0 restriction 'a nonphysical option.' The measured vibrational quenching rate is therefore an inclusive observable: it should include all inelastic ν=1→0 transitions into any final rotational state of C2- and any concurrent rotational change of H2. Under that inclusive comparison, the fully summed 5D rate coefficients shown in Figure 15 are about one to two orders of magnitude larger than the 20 K experimental point. The reported closeness of 5.0e-13 to 4.0e-13 cm3 molecule-1 s-1 is thus obtained by selecting the one final channel that happens to match the measurement, with no physical or experimental justification for why the detection would be blind to final states with j1'>0. This is a selection-bias concern about the central claim of agreement, not an internal inconsistency of the scattering calculation; the quantum dynamics, PES construction, and convergence checks are otherwise standard and transparently reported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs a new five-dimensional ab initio potential energy surface for C2-(X2Σg+) + H2, treating H2 as a rigid rotor and C2- as a vibrating-rotating diatomic, and solves the quantum close-coupling scattering problem to obtain state-resolved ro-vibrational quenching cross sections and rate coefficients for both para- and ortho-H2. The authors compare the resulting vibrational de-excitation rates with the single experimental datum at 20 K from their earlier work. They report that the full 5D rates, summed over final anion rotational states, are one to two orders of magnitude larger than the measured value, and that agreement with experiment is recovered only when the final C2- is restricted to its ground rotational state j1'=0, an option they later label as nonphysical.","tokens_in":27603,"tokens_out":6768,"duration_ms":67023,"significance":"The new 5D PES and the state-to-state close-coupling results are a substantial step beyond the earlier 3D treatment and provide a useful resource for buffer-gas cooling and possibly astrochemical modeling of C2-. The paper is transparent about the numerical setup, reports convergence checks, and makes the fitting data available in the Supplementary Information. The state-resolved rate distributions in Figures 14-16 are potentially valuable for planning state-specific experiments. However, the central claim of quantitative agreement with the measured 20 K quenching rate is not supported, because it rests on a post hoc restriction of the final anion rotational state that the authors themselves call nonphysical and for which the experimental detection scheme provides no justification.","major_comments":[{"comment":"The only point of agreement with experiment is the dashed red curve in Figure 15, constructed by summing vibrationally quenched final states with j1'=0 only. The authors state in Section IV.B that the experiments of Ref. [52] were not observing which rotational states are populated after the anion decays to its nu=0 level, and in Section V they call the j1'=0 restriction 'a nonphysical option.' The measured vibrational quenching rate is therefore an inclusive observable over all final C2- rotational states. When the fully summed 5D rate coefficients are used, as shown by the solid curves in Figure 15, they exceed the measured 20 K point by one to two orders of magnitude. The closeness of the reported 5.0e-13 cm3 molecule-1 s-1 to the experimental 4.0e-13 cm3 molecule-1 s-1 is thus a consequence of selecting the one final channel that matches the measurement. This is a load-bearing issue for the central claim of the paper, and it needs to be resolved either by providing an experimental or physical justification for why the detection is blind to j1'>0 final states, or by removing the agreement claim and reporting the computed state-resolved and summed rates without claiming validation by the experiment.","section":"Section IV.C, Figure 15, Section V"},{"comment":"The manuscript states that 'the vibrationally inelastic cross sections remained reasonably consistent and thus our R range is sufficiently large to obtain cross sections which are to the correct order of magnitude.' This convergence criterion is insufficient for the quantitative comparison made in Section IV.C, where the claimed agreement is a factor of about 1.25 (5.0e-13 versus 4.0e-13 cm3 molecule-1 s-1). The authors should provide quantitative convergence tests for the specific nu=1, j1=0 to nu'=0, j1'=0 channel, including convergence with respect to the number of propagator steps, R_max, and the size of the angular basis, rather than an order-of-magnitude consistency statement.","section":"Section III, vibrational convergence discussion"},{"comment":"The equation defining k_{nu->nu'}(T) appears malformed: it begins with '1P' and contains unbalanced parentheses and unclear placement of the degeneracy and Boltzmann factors. Since this equation underlies the construction of the thermally averaged and state-restricted curves in Figure 15, the authors must rewrite it clearly so that the averaging procedure is unambiguous and reproducible.","section":"Section IV.C, unnumbered equation after Eq. (16)"}],"minor_comments":[{"comment":"The maximum radial distance for the ab initio points is given as 25.0 Å in Section II but as 41 Å in Section III; please clarify which value is correct and how this affects the long-range switching at R0 = 20.5 Å.","section":"Section II.A and Section III"},{"comment":"In the text following Eq. (5), references to 'eq.(2)' and 'eq.(4)' appear to be misnumbered; the potential expansion is Eq. (5) and the bi-spherical harmonics are defined in Eq. (6). Please correct the cross-references.","section":"Section II.B"},{"comment":"The caption of Figure 15 and the discussion in Section IV.C describe the red and purple curves in opposite ways: the caption attributes the j2'=2 para-H2 excitation to the purple curve and the anion-only rotational channels to the red curve, while the main text says the red curve involves concurrent para-H2 rotational excitation and the purple curve involves only anion rotational states. These statements need to be reconciled.","section":"Figure 15 caption and main text"},{"comment":"The experimental point is shown without an uncertainty estimate; because the claimed agreement is quantitative, the authors should report the experimental error bar and state explicitly whether the agreement is within that uncertainty.","section":"Section IV.C, Figure 15"},{"comment":"The numerical values of the fully summed 5D rate coefficients and of the j1'=0-restricted rate at 20 K are only given in the text for the restricted case; a small table listing the summed rates for para-, ortho-, and the 3:1 weighted average at 20 K would make the comparison transparent.","section":"Section IV.C"}],"recommendation":"major_revision","confidential_remarks":"The new PES and the state-resolved scattering results are useful, and the manuscript is generally careful in its technical execution. However, the headline comparison with experiment is built on a channel restriction that the authors themselves describe as nonphysical, so the stated agreement cannot be taken as validation. The revision path is clear: either justify the j1'=0 restriction with a concrete detection mechanism, or reframe the paper to present the state-resolved predictions and the discrepancy with the inclusive experimental value as an open problem rather than as agreement. I therefore recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real contribution here is the 5D interaction potential and the state-to-state close-coupling rates. Treating H2 as a rigid rotor and C2- as a vibrating rotor is a clear step beyond the earlier 3D treatment, and the paper is transparent about the PES construction: 76k ab initio points, an ANN fit with documented residuals, a long-range switching function, and checks against the old 3D surface. The scattering calculations look carefully converged, with basis sets, partial waves, and energy grids all discussed. The para/ortho-H2 comparison is a useful addition, and the data appear to be available in the SI. This is a solid piece of computational cold-collision physics.\n\nThe soft spot is exactly where the reader's report lands. The only point of agreement with the measured 20 K rate is the dashed red curve in Figure 15, which restricts the final C2- to j1'=0. The authors explicitly say in Section IV.B that the experiment did not observe final rotational states, and in Section V they call the j1'=0 restriction a nonphysical option. The measured rate is inclusive, so the full summed 5D rates are the relevant comparison, and those sit one to two orders of magnitude above the experimental point. Choosing the one channel that happens to match, and then putting that agreement in the abstract, overstates what the theory actually explains. The authors deserve credit for flagging the issue, but they do not resolve it, and the framing of the paper's main claim should change.\n\nMinor points: the ANN fit has some large maximum deviations at high energies, though probably not consequential for the dynamics studied; the angular constraints for large r1 are a bit ad hoc but documented. The citation pattern leans heavily on the authors' own prior work, which is acceptable here since this is a direct extension.\n\nOverall: the computational machinery and the new PES are worth having, and the paper is honest about its main difficulty. But the headline comparison to experiment is fragile and should be reframed as a state-selected possibility rather than a confirmed match. A serious referee should see it; I would send it to peer review with a request for major revision.","headline":"A genuinely new 5D PES and state-resolved quenching rates for C2- + H2, but the paper's only agreement with the 20 K experiment comes from a final-state filter the authors themselves call nonphysical.","tokens_in":28212,"tokens_out":1462,"would_cite":true,"duration_ms":16889,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A new five-dimensional potential surface for C2- + H2 changes the predicted quenching rates by one to two orders of magnitude and brings theory into agreement with experiment only when the final anion is in its lowest rotational state.","keywords":["C2- anion","vibrational quenching","H2 buffer gas","5D potential energy surface","close-coupling scattering","rate coefficients","cold ion trap","laser cooling"],"falsifier":"Measure the final rotational state distribution of C2-($\\nu=0$) produced by H2 quenching in the trap, for example by state-selective photodetachment or photoelectron imaging after the quenching step. If a substantial population appears in $j_1'>0$ states, especially around $j_1'=6$, then the $j_1'=0$-only comparison is not the experimental channel and the claimed agreement would collapse; if the population is overwhelmingly $j_1'=0$, the agreement stands.","tokens_in":2000,"feed_emoji":"🧊","tokens_out":7218,"duration_ms":109976,"temperature":0.7,"pith_summary":"This paper tries to establish that the collisional vibrational quenching of the molecular anion C2- by H2 must be treated on a new five-dimensional (5D) interaction potential, not the earlier three-dimensional one, and that the added dimensions change the predicted rate coefficients by one to two orders of magnitude. The authors build the 5D ab initio surface, run close-coupling scattering calculations for the $\\nu=1\\to 0$ transition, and find that para-H2 cools more efficiently than ortho-H2. The full rotationally summed rate at 20 K overshoots the only measured value by a factor of up to a hundred, but when they restrict the sum to collisions that leave C2- in its ground rotational state $j_1'=0$, the weighted average rate of about $5.0\\times 10^{-13}$ cm$^3$ molecule$^{-1}$ s$^{-1}$ comes close to the experimental $4.0\\times 10^{-13}$. The authors label the $j_1'=0$-only option as nonphysical, so the central claim is that the 5D dynamics are right and the experiment is likely seeing only a subset of final rotational states.","feed_headline":"New 5D model matches measured C2- quenching rate at 20 K","feed_subtitle":"The full 5D rates overshoot by 10 to 100 times; agreement appears only when the anion ends in j=0.","key_machinery":"The load-bearing object is a new five-dimensional ab initio potential energy surface for C2-($X\\,^2\\Sigma_g^+$) + H2, built from 76,474 RCCSD(T)-F12 points and fitted by a (50,50,50) artificial neural network, with H2 frozen as a rigid rotor and C2- treated as a rotating-vibrating diatomic. The dynamics are run with close-coupled scattering using the vibrationally averaged coupling matrix elements $V_{\\nu\\nu'}(R)$ expanded in bispherical harmonics; the off-diagonal $V_{01}$ term is what drives the $\\nu=1\\to 0$ transition. What this machinery does is turn the question of how fast H2 quenches C2- vibration into a state-resolved statement about which final rotational states of the anion and of H2 are populated, which is exactly the information needed to compare with and interpret the trap experiment.","core_discovery":"On the paper's own terms, the central discovery is that realistic 5D dynamics reverses the earlier comparison with experiment: with H2 treated as a rigid rotor and the C2- bond allowed to vibrate, every rotationally summed quenching rate coefficient becomes larger than the earlier 3D results and larger than the measured 20 K point, often by one to two orders of magnitude. The purely vibrational channel that leaves C2- in $j_1'=0$ is the slowest of all the quenching paths, yet it is the one whose Boltzmann-averaged, ortho/para-weighted rate, about $5.0\\times 10^{-13}$ cm$^3$ molecule$^{-1}$ s$^{-1}$ at 20 K, sits close to the experimental $4.0\\times 10^{-13}$. The paper therefore claims that the state-resolved rates are trustworthy enough to identify which final rotational states are being populated, and that the experimental rate likely corresponds to a subset of final states rather than to the full rotationally summed rate.","pith_inferences":["If the final-state restriction turns out to be wrong and the experiment sums over all $j_1'$, then the paper's own numbers imply the measured 20 K rate should be one to two orders of magnitude larger than observed; that would point to a missing loss or detection channel in the experiment rather than to a failure of the surface.","A direct test would be rotationally resolved detection of C2-($\\nu=0$) after the quenching pulse; the computed state-to-state distributions predict specific peak final states, such as $j_1'=6$ for cold initial states, which a photoelectron or action spectrum could look for.","The same 5D machinery, with C2- vibrational states extended beyond $\\nu=1$, could supply the repumping-cycle rates needed to design closed laser-cooling schemes, since the bottleneck is knowing how fast each excited vibrational level is quenched by the buffer gas."],"forward_implications":["The full 5D rate coefficients for C2-($\\nu=1$) quenching by H2 are one to two orders of magnitude larger than the earlier 3D values, so any future cooling model that uses the 3D rates will underestimate the vibrational quenching efficiency.","para-H2 is a more efficient buffer-gas partner than ortho-H2 for cooling C2- vibration, whether or not the H2 rotor is excited during the collision.","The state-resolved calculations identify which final rotational states of the anion are populated after quenching; the distribution peaks at $j_1'=6$ for initial $j_1=0$ and shifts to higher $j_1'$ for hotter initial states.","If the experimental signal really comes from anions that end in $j_1'=0$, the computed 20 K rate of about $5.0\\times 10^{-13}$ cm$^3$ molecule$^{-1}$ s$^{-1}$ matches the measured value, supporting the 5D surface and the close-coupling dynamics."],"supporting_citations":[{"why":"Supplies the earlier 3D potential, the measured 20 K rate coefficient, and the comparison the 5D results are meant to supersede.","marker":"[52]"},{"why":"Supplies the coupled-channel scattering program used to propagate the 5D dynamics.","marker":"[63]"},{"why":"Sets out the coupled-channel scattering theory for diatom-diatom collisions on which the dynamics are based.","marker":"[64]"},{"why":"Provides the earlier scattering code whose 3D results are reproduced to validate the reduced-dimension limit of the new surface.","marker":"[65]"},{"why":"Provides the earlier He, Ne, Ar rovibrational quenching rates used to emphasize how much faster H2 quenches the anion.","marker":"[51]"},{"why":"Supplies the pseudo-singlet treatment of the doublet anion that the present scattering calculations adopt.","marker":"[49]"},{"why":"Identifies C2- as a laser-cooling candidate, giving the practical reason the quenching rates matter.","marker":"[5]"},{"why":"Shows the vibrational levels of C2- live for seconds, so collisional quenching is the only efficient removal path.","marker":"[24]"}],"fun_headline_variants":["State-specific 5D rates match C2- experiment","C2- quenching: only j=0 channel explains 20 K data","5D model shows total C2- rate overshoots 100x","Para-H2 more efficient than ortho for C2- cooling","New 5D PES resolves C2- laser cooling rate puzzle"],"cache_read_input_tokens":30208,"weakest_assumption_plain":"The whole agreement with experiment rests on assuming that the measured quenching rate counts only collisions that leave the C2- anion in its lowest rotational state $j_1'=0$, with no rotational excitation; the paper itself calls this choice nonphysical, and if the experiment actually collects anions in all final rotational states the computed rates overshoot the data by one to two orders of magnitude.","fun_headline_variants_meta":{"raw":{"variants":["State-specific 5D rates match C2- experiment","C2- quenching: only j=0 channel explains 20 K data","5D model shows total C2- rate overshoots 100x","Para-H2 more efficient than ortho for C2- cooling","New 5D PES resolves C2- laser cooling rate puzzle"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000455,"raw_usage":{"total_tokens":2361,"prompt_tokens":1099,"completion_tokens":1262,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":715,"completion_tokens_details":{"reasoning_tokens":1168}},"tokens_in":715,"tokens_out":1262,"duration_ms":10529,"temperature":1.0,"reasoning_tokens":1168,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:30:28.728668+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the final rotational state distribution of C2-($\\nu=0$) produced by H2 quenching in the trap, for example by state-selective photodetachment or photoelectron imaging after the quenching step. If a substantial population appears in $j_1'>0$ states, especially around $j_1'=6$, then the $j_1'=0$-only comparison is not the experimental channel and the claimed agreement would collapse; if the population is overwhelmingly $j_1'=0$, the agreement stands.","supporting_citations":[{"cited_title":"Fa¨ y and H","cited_arxiv_id":null,"evidence_quote":"Supplies the coupled-channel scattering program used to propagate the 5D dynamics."},{"cited_title":"Civi˘ s, Y","cited_arxiv_id":null,"evidence_quote":"Sets out the coupled-channel scattering theory for diatom-diatom collisions on which the dynamics are based."},{"cited_title":"L´ opez-Dur´ ann, E","cited_arxiv_id":null,"evidence_quote":"Provides the earlier scattering code whose 3D results are reproduced to validate the reduced-dimension limit of the new surface."},{"cited_title":"Rosmus and H.-J","cited_arxiv_id":null,"evidence_quote":"Provides the earlier He, Ne, Ar rovibrational quenching rates used to emphasize how much faster H2 quenches the anion."},{"cited_title":"Zeitz, S","cited_arxiv_id":null,"evidence_quote":"Supplies the pseudo-singlet treatment of the doublet anion that the present scattering calculations adopt."},{"cited_title":"Yzombard, M","cited_arxiv_id":null,"evidence_quote":"Identifies C2- as a laser-cooling candidate, giving the practical reason the quenching rates matter."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows the vibrational levels of C2- live for seconds, so collisional quenching is the only efficient removal path."}],"review_version":1}