{"id":"7ea63821-db79-4a8a-b1ce-b74b7c98362a","arxiv_id":"2608.05067","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"First systematic ab initio study of 57 alkali-metal and alkali-alkaline-earth-metal diatomic anions, predicting ground-state curves, dipole moments, polarizabilities, and select excited states.","lead":"This paper presents high-level quantum-chemistry calculations for 57 molecular anions made of alkali and alkaline-earth atoms, giving their binding energies, bond lengths, dipole moments, and polarizabilities. The data fill a gap for cold-molecule research and include predicted crossings with neutral molecules that could mediate resonant electron attachment.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Accuracy claim rests on a validation set that excludes the heaviest species; neglect of spin–orbit coupling for Fr- and Ra-containing anions is the least secure link in the extrapolation to 'few-percent' accuracy for all 57 systems.","rationale":"The reader's weakest-assumption statement identifies the same load-bearing point: the composite method is validated on only KRb− and RbSr− and is assumed uniform across all 57 anions, with spin–orbit coupling neglected for Fr and Ra. I agree, and I narrow the concern to the specific extrapolation step that is least supported. The existing experimental benchmarks for alkali dimer anions are a real strength; they cover nine species and give a 3.4% mean deviation, which supports the method for Li–Cs systems. However, none of those species contains Fr or Ra, so the empirical anchor does not constrain the heavy end. The atomic benchmarks in Table I show the expected degradation for Fr/Ra, making it plausible that molecular well depths could shift by more than the quoted few percent once spin–orbit and higher-order relativistic effects are included. The paper's central contribution is a uniform reference dataset; if the heavy-subset uncertainties are underestimated, the headline claim of first accurate predictions is only partially supported. This does not invalidate the paper; it argues for the same conditional acceptance the reader recommended, perhaps with the condition that a heavy-element benchmark or explicit uncertainty per species be added. The missing Supplemental Material is a reproducibility concern, but it is secondary because the numerical PECs will presumably be provided by the publisher; the accuracy extrapolation is the deeper issue.","tokens_in":25975,"tokens_out":7538,"duration_ms":88779,"concrete_test":"Repeat the Section III E convergence protocol for at least one Fr-containing and one Ra-containing anion, e.g. FrBe− and FrRa− or Fr2−, using a spin–orbit-coupled electronic-structure treatment such as two-component CCSD(T) with spin–orbit pseudopotentials or MRCI+Q with the Breit–Pauli operator, and compare Re, De, and ωe against the scalar CCSD(T)+ΔT values. If the spin–orbit-induced shift in De exceeds the claimed few-percent band (e.g., more than roughly 100 cm−1 for a ~3000 cm−1 well), the uniform few-percent accuracy assertion fails for the heavy subset, and the manuscript would need either explicit spin–orbit corrections or per-species error bars before the Fr/Ra predictions can be called reference quality.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III E validates the composite CCSD(T)+ΔT scheme of Eq. (1) on only two species, KRb− and RbSr−, and all other 55 PECs inherit the statement that well depths are accurate to within a few percent. The extrapolation is least secure at the heavy end: Fr and Ra are described with scalar-relativistic ECPs and no spin–orbit coupling, and the atomic benchmarks in Table I already show larger deviations for these elements (e.g., Fr ionization potential off by 421 cm−1 and Fr S–P excitation energy by 300 cm−1, even after spin–orbit averaging). The ground states are Σ states, so first-order spin–orbit effects vanish, but second-order spin–orbit coupling to low-lying Π states grows steeply with Z and is entirely absent from the error budget. Critically, the nine alkali-dimer anions with experimental well depths that support the reported 3.4% mean absolute deviation contain no Fr or Ra species, so that empirical check does not bound the error for the 12 Fr/Ra-containing anions. If spin–orbit shifts a heavy-anion potential by more than a few percent, the 'first accurate predictions' label and the reference-quality claim for those species are not yet supported. A per-species uncertainty estimate, or at least one heavy-element benchmark, is missing.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports a systematic computational study of 57 diatomic molecular anions formed from alkali-metal and alkali-metal–alkaline-earth-metal atoms. Using a composite CCSD(T)+ΔT scheme with large Gaussian basis sets and small-core relativistic energy-consistent pseudopotentials, the authors compute ground-state potential energy curves, spectroscopic constants, permanent dipole moments, and static electric dipole polarizabilities for all species. They also investigate selected excited states with MRCISD and EOM-EA-CCSD, including dipole-bound states for six polar molecules, and predict crossings between neutral ground-state and excited anionic-state curves. The results are benchmarked against atomic properties and experimental well depths for nine alkali-dimer anions (mean absolute deviation 3.4%), and basis-set and correlation-convergence tests are performed for KRb− and RbSr−.","tokens_in":26195,"tokens_out":8947,"duration_ms":83460,"significance":"If the accuracy claims hold, this is a valuable reference dataset: it provides the first high-level predictions for most of the alkali-metal–alkaline-earth-metal anions, which are largely unexplored experimentally, and it identifies several strongly polar anions with large dipole moments and potential dipole-bound states. The methodological validation is a notable strength: no parameters are fitted to the target molecular data, the atomic benchmarks are external, and the basis-set and correlation convergence is demonstrated for representative species. The main weakness is that the global few-percent accuracy claim is extrapolated to the heaviest (Fr- and Ra-containing) species without a dedicated benchmark or spin–orbit treatment, and the crossing prediction rests on unvalidated MRCISD excited-state curves.","major_comments":[{"comment":"The statement in Section IV that the calculated well depths are accurate to within a few percent is not supported for the 12 anions containing Fr or Ra. The only molecular validations are the two convergence tests on KRb− and RbSr− in Section III E and the empirical comparison with nine alkali-dimer anions in Table II, none of which contains Fr or Ra. The atomic benchmarks in Table I show the largest deviations for Fr (ionization potential off by 421 cm⁻¹, S–P excitation energy by 300 cm⁻¹) and Ra (IP off by 365 cm⁻¹), and the method uses scalar-relativistic ECPs with no spin–orbit coupling. Since second-order spin–orbit effects on Σ states grow with Z, an error of several percent in the well depths of Fr/Ra-containing anions cannot be excluded. I recommend adding a benchmark for at least one heavy system (e.g., CsFr− or FrBa− with spin–orbit treatment) or, failing that, restricting the few-percent claim to species without Fr/Ra and assigning a larger uncertainty to the heavy end.","section":"§III E and §IV (surface accuracy claim)"},{"comment":"The predicted crossings between the neutral ground state and the excited A²Σ+ anionic state, which motivate the resonant-electron-attachment scenario in the abstract and conclusions, are based on MRCISD curves for only Li₂⁻, RbCs⁻, and NaCs⁻. The authors appropriately label the continuum-embedded portions as approximate diabatic continuations, but they do not quantify how the crossing position depends on the level of theory or on the uncertainty of the atomic asymptotes. Given that the atomic excitation energies in Table I carry errors of up to 300 cm⁻¹, the crossing could shift significantly or disappear for some species. A sensitivity analysis (e.g., shifting the curves by the atomic EA/excitation-energy uncertainties, or comparing the A-state with EOM-EA-CCSD) should be included before the crossing prediction is presented as a robust finding.","section":"§III B, Fig. 1"},{"comment":"Although the abstract states that the authors assess convergence and uncertainties of their results, the uncertainty analysis is limited to two representative species and a global few-percent statement. Tables II and III present D_e, R_e, ω_e, dipole moments, and polarizabilities for 57 species without any per-species uncertainty estimate. The spread among the CI variants in Fig. 7 and the basis-set differences in Fig. 6 suggest that the error could vary substantially across the periodic table (e.g., the ΔT correction ranges from 2.4% to 3.4% between the two classes, and the CBS extrapolation residual differs between KRb− and RbSr−). I recommend reporting at least the estimated absolute uncertainty in D_e and R_e for each species, or providing a clear scaling relationship based on the atomic benchmarks.","section":"§III D–E and Tables II–III"}],"minor_comments":[{"comment":"Reference [100] appears to be malformed: the bibliographic entry is replaced by a list of bond-function exponents. Please restore the proper citation (apparently to Tao and Pan, J. Chem. Phys. 97, 4989 (1992)) and move the exponents to the text or Supplemental Material.","section":"Reference [100]"},{"comment":"The text says that the electronic adiabatic electron affinities use the experimental electron affinities of the constituent atoms, but for Fr the Table I value is theoretical (Ref. [125]). Please clarify that for species with no experimental atomic EA, the calculated atomic EA is used.","section":"§III D, electron affinity formula"},{"comment":"The 3.4% mean absolute deviation is a central validation number, but Table II lists two experimental well depths for Rb₂⁻ and Cs₂⁻ (from Refs. [76] and [77]) without specifying which is used in the nine-species statistics. Please state explicitly which experimental values enter the mean deviation and error budget.","section":"Table II and §III D1"},{"comment":"The basis-set label 'apVQZ' is used without definition in the figure caption or text; please define it (e.g., as aug-cc-pVQZ) and clarify that the diffuse functions are the even-tempered sets described in Section II.","section":"Figure 3 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the main dataset is potentially valuable. The two main concerns are the unsupported few-percent accuracy claim for Fr/Ra-containing species and the lack of a sensitivity analysis for the excited-state crossings. The reliance on the authors' own previous paper (Ref. [96]) for the neutral curves in the electron-affinity formula is a notable dependency; it is not circular for the PECs, but the EAs quoted here inherit any errors in that earlier work. I would encourage the editor to request the additional benchmarks or a clearly scoped accuracy statement before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a competent computational survey that delivers what it promises: ground-state PECs, spectroscopic constants, dipole moments, and polarizabilities for 57 molecular anions, most of which had no prior accurate data. The composite CCSD(T)+ΔT scheme is standard but well executed, and the validation is genuinely useful—nine experimental alkali-dimer anion well depths are reproduced with a 3.4% mean absolute deviation, and the basis-set/correlation convergence tests on KRb− and RbSr− are thorough. No parameters are fitted to the target molecules, so the circularity burden is low.\n\nThe soft spots are real but not disqualifying. The accuracy claim is extrapolated from two validation systems to all 57 species, and the heaviest cases (Fr and Ra) are exactly where the evidence is thinnest: no experimental anion well depth includes Fr or Ra, spin–orbit coupling is neglected beyond the scalar ECPs, and the atomic benchmarks in Table I already show the largest deviations for Fr. Because the ground states are Σ states, first-order spin–orbit effects vanish, but second-order coupling to low-lying Π states can matter at this level of precision and it is absent from the error budget. That does not undermine the dataset for the lighter species, but it does mean the \"few-percent\" label should be qualified for the 12 Fr/Ra-containing anions. A per-species uncertainty estimate, or at least one heavy-element benchmark, would close this cleanly.\n\nThe excited-state and dipole-bound-state sections are more exploratory. The predicted neutral–anion crossings are clearly flagged as suggestive rather than quantitative, and the EOM-EA-CCSD binding energies for six dipole-bound states are well characterized with basis convergence. These are nice additions, not the core of the paper. One practical gripe: the arXiv version points to Supplemental Material that is not included, so the actual numerical PECs are not yet available to the reader.\n\nThis paper deserves a serious referee. The central ground-state dataset is a useful reference for cold-molecule spectroscopy, electron-attachment experiments, and future theory. I would accept it after minor revision, with the request to either provide a heavy-element benchmark or soften the accuracy claim for those species, and to make the supplemental data available with the preprint.","headline":"A solid, comprehensive dataset paper for 57 alkali/alkaline-earth diatomic anions, with a real but localized soft spot: the few-percent accuracy claim is extrapolated to the heavy Fr/Ra species without a spin-orbit or heavy-element benchmark.","tokens_in":665,"tokens_out":792,"would_cite":true,"duration_ms":27152,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["31.15.A-","31.15.Ar","33.15.-e"],"model":"deepseek-v4-flash","headline":"A uniform computational survey maps the ground states of all 57 alkali-metal and alkali-alkaline-earth diatomic anions.","keywords":["molecular anions","alkali-metal dimers","alkaline-earth atoms","coupled cluster","potential energy curves","dipole-bound states","electron affinity","ultracold collisions"],"falsifier":"Photoelectron spectroscopy of an unstudied heteronuclear anion such as LiBa$^-$ or CsRa$^-$ would measure the adiabatic electron affinity and vibrational spacings; if the measured binding energy differed from the predicted value by more than roughly 200 cm$^{-1}$ for LiBa$^-$, the uniform accuracy claim would fail for that class. Alternatively, a microwave or Rydberg-spectroscopy measurement of the dipole-bound state of NaCs$^-$ could test the predicted binding energy directly.","tokens_in":25764,"feed_emoji":"⚛️","tokens_out":4768,"duration_ms":42546,"temperature":0.7,"pith_summary":"This paper sets out to provide a uniform, high-accuracy computational description of all 57 diatomic anions formed from an alkali-metal atom plus either a second alkali-metal atom or an alkaline-earth-metal atom. For most of these species no reliable potential energy curves or spectroscopic constants existed before, so the paper aims to supply reference data where experiment and theory are missing. It computes ground-state potential curves, well depths, vibrational and rotational constants, permanent dipole moments, and static polarizabilities, along with excited valence and dipole-bound states for selected species. If the calculations are right, these molecules become accessible to planned cold-collision and Rydberg-electron experiments, and photoelectron spectra of the anions can be assigned against a consistent theoretical grid.","feed_headline":"First accurate maps of 57 alkali-metal molecular anions","feed_subtitle":"Predicted potentials, dipoles, and polarizabilities give ultracold and Rydberg experiments a roadmap to these species.","key_machinery":"The load-bearing object is the composite interaction energy $V_{\\mathrm{int}}(R) = V^{\\mathrm{apwCV5Z+bf}}_{\\mathrm{CCSD(T)}}(R) + \\delta V^{\\mathrm{apwCVTZ}}_{\\mathrm{CCSDT}}(R)$, where the second term adds the iterative-triples correction computed as CCSDT minus CCSD(T) in a smaller basis. Around this the paper wraps counterpoise-corrected supermolecule energies, a core-valence quintuple-zeta basis with bond functions, Stuttgart small-core pseudopotentials, and a custom even-tempered diffuse Gaussian sequence at the bond midpoint for dipole-bound states. EOM-EA-CCSD supplies the dipole-bound electron binding energies, MRCISD supplies the excited-state curves, and Dunham fits turn the potentials into spectroscopic constants.","core_discovery":"The central claim is that a composite CCSD(T)+ΔT method, with large core-valence basis sets and small-core relativistic pseudopotentials, yields ground-state potential energy curves accurate to a few percent for all 21 alkali-metal diatomic anions (X$^{2}\\Sigma^{+}$) and all 36 alkali-metal--alkaline-earth-metal anions (X$^{1}\\Sigma^{+}$). The paper further claims that the excited A$^{2}\\Sigma^{+}$ states of the alkali-metal anions cross the neutral ground-state curves, producing temporary anionic states embedded in the electron-detachment continuum, and that six polar molecules support dipole-bound states whose binding energies grow near-linearly with the neutral dipole moment. These crossings and dipole-bound states are presented as a mechanism for resonant electron attachment in ultracold mixtures of ground-state molecules and Rydberg atoms.","pith_inferences":["If the crossing picture is right, the same temporary-anion mechanism should operate for other polar neutrals with sufficiently large excited-state dipole moments, extending the Rydberg-attachment idea beyond the alkali-metal set.","The near-linear dipole-bound binding energy versus dipole moment correlation, demonstrated here for six related molecules, could be tested as a scaling rule for other polar molecules and would let experimentalists estimate binding from a measured dipole alone.","The absence of spin-orbit coupling for Fr and Ra means the true heavy-anion curves may split into multiple components; resolving this would require a two-component treatment the present data cannot distinguish.","The predicted data could underpin a search for laser-coolable molecular anions: any anion with a bound excited state of opposite parity within the alkali set would be a candidate for optical cycling."],"forward_implications":["Photoelectron spectra of the nine experimentally studied anions can be re-examined against a single consistent set of curves, and the remaining 48 species now have predicted well depths and vibrational spacings to test.","The predicted neutral--anion curve crossings identify internuclear distances where electron attachment should be resonantly enhanced, giving Rydberg--molecule experiments specific targets.","Dipole-bound binding energies for NaCs$^-$ and LiCs$^-$ are large enough relative to the rotational constants that many rotational levels should lie below the detachment threshold, a regime open to high-resolution spectroscopy.","The reported dipole moments and polarizabilities provide the input needed to estimate trap-induced Stark shifts and long-range interactions in hybrid ion-atom systems.","The uniform treatment makes possible systematic trends---well depth falling with increasing atomic size, Ba-containing anions deepest---that can guide which species to try to form and cool."],"supporting_citations":[{"why":"Supplies the neutral-molecule potential energy curves, well depths, and dipole moments used to define anionic adiabatic electron affinities and dipole-bound thresholds.","marker":"[96]"},{"why":"Li$_2^-$ photoelectron experiment providing the principal experimental benchmark for the alkali-metal dimer anions.","marker":"[75]"},{"why":"Photoelectron measurements for Na$_2^-$, NaK$^-$, K$_2^-$, KRb$^-$, KCs$^-$, Rb$_2^-$, RbCs$^-$, Cs$_2^-$ against which well depths and electron affinities are validated.","marker":"[77]"},{"why":"Earlier theoretical Rb$_2^-$ and Cs$_2^-$ curves used for comparison in the alkali-metal series.","marker":"[90]"},{"why":"Prior calculations for LiBe$^-$, NaBe$^-$, LiMg$^-$, NaMg$^-$ used to benchmark the alkaline-earth series.","marker":"[91]"},{"why":"Previous KCa$^-$ calculation compared against the present result.","marker":"[92]"},{"why":"EOM-EA-CCSD method used to compute dipole-bound electron binding energies.","marker":"[108]"},{"why":"Experimental observation of dipole-bound states in KI$^-$ motivating and contextualizing the dipole-bound-state predictions.","marker":"[52]"}],"fun_headline_variants":["57 alkali-metal anion states mapped to few-percent accuracy","Anion-neutral crossings may enable Rydberg electron capture","Dipole-bound states link polar anions to ultracold experiments","First accurate potentials for 57 diatomic anions","All 57 alkali-metal diatomic anion curves computed"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire dataset inherits the accuracy of a single composite method whose few-percent error bar is inferred from two representative anions (KRb$^-$ and RbSr$^-$) and from atomic benchmarks, with spin-orbit coupling neglected for the heaviest elements, so individual heavy-anion well depths could lie outside the stated uncertainty.","fun_headline_variants_meta":{"raw":{"variants":["57 alkali-metal anion states mapped to few-percent accuracy","Anion-neutral crossings may enable Rydberg electron capture","Dipole-bound states link polar anions to ultracold experiments","First accurate potentials for 57 diatomic anions","All 57 alkali-metal diatomic anion curves computed"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00057,"raw_usage":{"total_tokens":2702,"prompt_tokens":959,"completion_tokens":1743,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":575,"completion_tokens_details":{"reasoning_tokens":1678}},"tokens_in":575,"tokens_out":1743,"duration_ms":11445,"temperature":1.0,"reasoning_tokens":1678,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T10:03:30.006393+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Photoelectron spectroscopy of an unstudied heteronuclear anion such as LiBa$^-$ or CsRa$^-$ would measure the adiabatic electron affinity and vibrational spacings; if the measured binding energy differed from the predicted value by more than roughly 200 cm$^{-1}$ for LiBa$^-$, the uniform accuracy claim would fail for that class. Alternatively, a microwave or Rydberg-spectroscopy measurement of the dipole-bound state of NaCs$^-$ could test the predicted binding energy directly.","supporting_citations":[{"cited_title":"Gronowski, A","cited_arxiv_id":null,"evidence_quote":"Supplies the neutral-molecule potential energy curves, well depths, and dipole moments used to define anionic adiabatic electron affinities and dipole-bound thresholds."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Li$_2^-$ photoelectron experiment providing the principal experimental benchmark for the alkali-metal dimer anions."},{"cited_title":"Dobrzyniecki and M","cited_arxiv_id":null,"evidence_quote":"Photoelectron measurements for Na$_2^-$, NaK$^-$, K$_2^-$, KRb$^-$, KCs$^-$, Rb$_2^-$, RbCs$^-$, Cs$_2^-$ against which well depths and electron affinities are validated."},{"cited_title":"Nasiri and M","cited_arxiv_id":null,"evidence_quote":"Earlier theoretical Rb$_2^-$ and Cs$_2^-$ curves used for comparison in the alkali-metal series."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Prior calculations for LiBe$^-$, NaBe$^-$, LiMg$^-$, NaMg$^-$ used to benchmark the alkaline-earth series."},{"cited_title":"Partridge, D","cited_arxiv_id":null,"evidence_quote":"Previous KCa$^-$ calculation compared against the present result."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"EOM-EA-CCSD method used to compute dipole-bound electron binding energies."}],"review_version":1}