{"id":"546f1758-059c-4913-b84b-d72eebe73f10","arxiv_id":"2411.09375","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Vibrationally resolved electron-D2 resonant cross sections, including transitions between excited electronic states, are computed with the local complex potential method and compared to experiment.","lead":"This paper computes vibrationally resolved cross sections for electron collisions with deuterium molecules, covering vibrational excitation, electron attachment, and dissociation through three temporary negative-ion states. The data are aimed at fusion divertor and negative-ion-source models, where deuterium-specific collision rates are currently sparse.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 15% uncertainty claim is extrapolated to high-v and excited-state channels with no benchmark; the dataset's fusion-applicability rests on exactly these unvalidated channels.","rationale":"The paper is a credible computational application of an established method. The v=0 comparisons in Figure 5 (elastic, VE1, DA1/DA2, DE1/DE2) provide real validation, and the H2/D2 isotope comparison in Figure 6 is a useful sanity check, though it covers only a few channels. I do not fully share the reader's choice of weakest assumption as the primary concern: the neutral X, b, and B potential curves and the D2− widths are electronic quantities that are, to first order, isotopically invariant, so using H2-based molecular data with the D2 reduced mass captures most of the isotope physics. The more load-bearing premise is the uniform uncertainty claim. Every experimental and literature comparison in Section III.A is for the ground electronic state and v=0; the channels that make the dataset novel, namely high-v ground-state transitions and all B-state-initiated processes, are never benchmarked, yet the paper asserts the same 15% uncertainty for them and presents them as ready for fusion plasma applications. This is a correctness risk rather than a demonstrated inconsistency: the LCP equations are standard, the v=0 comparisons are plausible, and no internal contradiction is visible, but the central claim that the full dataset carries validated uncertainty is stronger than the evidence. A single targeted comparison of a few unvalidated channels with an independent nonlocal or R-matrix calculation would settle whether the 15% bound actually holds across the full dataset. I therefore keep the reader's CONDITIONAL verdict unchanged.","tokens_in":14963,"tokens_out":12455,"duration_ms":136898,"concrete_test":"Run an independent nonlocal-resonance calculation (e.g., Horáček et al. 2006) for D2 and compare VE1(10→11), DA1(v=10), VE3(0→1), and DE5(v=0) against the present LCP results. If any magnitude, peak position, or near-threshold structure shifts by more than ~15%, the uniform uncertainty claim is refuted and the dataset should be flagged as validated only for the v=0 ground-state channels.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III.A validates the calculation exclusively against v=0 ground-state channels: elastic, VE1(0→1, 0→2, 0→3), DA1/DA2 v=0, and DE1/DE2 v=0. No experimental or independent theoretical benchmark is shown for any v≥1 initial level or for any B-state-initiated process (VE3, DA3, DE4–DE6), although these unvalidated channels are the unique content of the dataset and dominate the rate behavior in Figures 3 and 4. The statement in Section III.A that 'we can assume for the presented cross sections an uncertainty of the same order of magnitude' therefore rests on an extrapolation across vibrational quantum number, electronic state, and process type. The local-complex-potential approximation and the treatment of the two 2Σ_g^+ resonances (B− and C−) as independent, incoherently summed local resonances are least controlled in exactly this unvalidated region: near-threshold dissociation from high v and the narrow C− Rydberg resonance. If local-width errors, same-symmetry resonance interference, or continuum normalization shift these channels, the claimed 15% uncertainty and the statement that the dataset is directly usable for fusion divertor and negative-ion-source modeling fail, even though the v=0 comparisons look good.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents local-complex-potential (LCP) calculations of vibrationally resolved resonant electron-collision cross sections for D2: vibrational excitation (VE), dissociative attachment (DA), and dissociative excitation (DE), including transitions among the ground X 1Σg+ state and the first two excited states b 3Σu+ and B 1Σu+, proceeding through the X-, B-, and C- resonances of D2-. The authors compare selected v=0 ground-state channels with experimental and theoretical data, report agreement of roughly 15%, discuss the H2/D2 isotopologue effect, and make the full dataset available through the IAEA HCDB and LXCat databases.","tokens_in":15195,"tokens_out":3952,"duration_ms":41324,"significance":"If the claimed completeness and accuracy hold, the paper would supply a useful and previously missing dataset for fusion divertor and negative-ion-source modeling, since it extends the LCP method to vibrational transitions between electronic excited states for the first time. The availability of the results in HCDB and LXCat is a concrete strength for downstream plasma-kinetics use. However, the central accuracy assertion rests on an extrapolation from a small set of v=0 ground-state validation cases to the full dataset, and the input resonance data are taken from H2 calculations without isotope-specific recomputation, so the significance hinges on these unvalidated assumptions.","major_comments":[{"comment":"The resonance potentials and autoionization widths for D2- are all taken from H2-based calculations (refs. [49,50] for X- and B-, [51] for C-); only the nuclear reduced mass is changed to that of D2. Since V_s^r in Eq. (5) is proportional to the square root of the partial width, any isotope dependence of the widths propagates directly into every cross section defined by Eqs. (1)-(3). The paper should either provide evidence that these widths are isotope-independent for the present purposes or estimate the sensitivity of the final cross sections to plausible isotopic changes in the widths. At minimum, this transfer of input data should be stated as an explicit uncertainty source rather than only as a modeling choice.","section":"Section II, Eqs. (4)-(6) and Figure 1"},{"comment":"All validation comparisons are for v=0 initial states of the X state: elastic scattering, VE1(0→1, 0→2, 0→3), DA1/DA2, and DE1/DE2. No comparison or independent benchmark is shown for any v≥1 initial level or for any B-state-initiated process (VE3, DA3, DE4-DE6), although these unvalidated channels are the novel content of the dataset and dominate the rate behavior shown in Figures 3 and 4. The statement that 'we can assume for the presented cross sections an uncertainty of the same order of magnitude' is therefore an extrapolation across vibrational quantum number, electronic state, and process type. The authors should provide at least spot checks for a few high-v or excited-state channels, or alternatively reformulate the uncertainty claim as applying only to the validated v=0 channels.","section":"Section III.A, Figure 5"},{"comment":"Equations (1)-(3) sum incoherently over the three resonances r. Two of these, B- and C-, have the same 2Σg+ symmetry, and the paper treats them as independent local resonances. For overlapping resonances of the same symmetry, interference effects can be significant, particularly for the narrow C- Rydberg resonance that produces the sharp structures near 11-12 eV in Figure 2. The incoherent-sum approximation should be justified or tested, for example by comparing with a two-resonance nonlocal model or by examining whether the local-width treatment is valid for the Rydberg resonance. Without such a test, the accuracy of the high-energy structures in the VE and DE channels remains uncontrolled.","section":"Section II, Eqs. (1)-(3)"},{"comment":"The paper presents the dataset as 'full sets of cross sections' for fusion plasma applications, but the calculations include only resonant channels. For elastic scattering and DE, the paper itself states that non-resonant contributions dominate at higher energies and that the LCP results underestimate the cross sections there. For VE through the B/C- states, the non-resonant direct excitation contribution is also not quantified. The completeness claim should be qualified in the abstract and conclusions by specifying the energy range in which the resonant cross sections are intended to be used and by stating that the dataset must be supplemented with non-resonant cross sections outside that range.","section":"Section III.A and Table I"}],"minor_comments":[{"comment":"The text refers to 'Frank-Condon factors'; the correct spelling is 'Franck-Condon factors'.","section":"Section III (text near Figure 2)"},{"comment":"There are several typographical errors in the closing paragraph: 'improuve' should be 'improve', and 'intent to' should be 'intend to'.","section":"Section IV, Conclusions"},{"comment":"The sentence 'integration over internuclear distances carrier out over the interval R ∈ [0.4, 15] a.u.' should read 'is carried out'.","section":"Section II, paragraph after Eq. (6)"},{"comment":"The references cited for the VE1 comparisons are listed as 'Phelps et al.' in one place and 'Buckman and Phelps' in another; the names should be used consistently throughout the text and figure caption.","section":"Section III.A, Figure 5 discussion"},{"comment":"The sentence 'the thresholds for DA and DE processes for H2 decrease compared with those for D2' is grammatically confusing; it should be rephrased to clarify that the threshold electron energies are lower for H2 than for D2 at the same vibrational quantum number.","section":"Section III.B"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the central LCP calculation is not circular; the v=0 comparisons are encouraging. The main risk is that the 15% uncertainty is presented as global when only v=0 ground-state channels are validated, and the input widths are H2-based. These are fixable with additional sensitivity tests or more cautious wording, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a competent application of the established LCP method to D2, producing a genuinely new vibrationally resolved dataset including transitions involving the b and B excited states. The v=0 ground-state validation against Golden, Buckman-Phelps, Schulz, Rapp, Krishnakumar, and the AN calculations looks good, and the data are deposited in IAEA and LXCat, which makes them immediately usable. The central claim—that these are the first LCP results for vibrational transitions and dissociation between different electronic states of D2—checks out.\n\nWhat is new: the full VE1/VE2/VE3 and DA1-DA3, DE1-DE6 vibrationally resolved sets for D2, especially B-state initiated channels. That fills a real gap for divertor and negative-ion-source modeling, where people currently use H2 data or older D2 estimates. The isotopologue comparison with H2, including the factor ~200 suppression of DA1 at v=0, is a useful cross-check.\n\nSoft spots, in proportion: the 15% uncertainty claim is extrapolated. The validation section is entirely v=0 ground-state channels; there is no benchmark for v≥1 initial levels or any B-state channel (VE3, DA3, DE4-DE6), which are precisely the unique content. The statement that \"we can assume ... same order of magnitude\" is not backed by propagated errors or an independent test. Also, the input widths and resonance potentials are taken from prior H2-based calculations without recomputation; that is physically defensible at Born-Oppenheimer level, but it should be stated as an assumption with an estimated sensitivity, especially for the narrow C− Rydberg resonance. The omission of non-resonant background is clearly disclosed and is appropriate for the low-energy resonant channels, but users need to remember the dataset is not total cross sections.\n\nThe stress-test concern about extrapolation is fair and should be addressed by the authors—ideally with a few high-v benchmarks or a sensitivity scan on the widths—but it does not undermine the basic calculation. The equations are standard, the code lineage is established, and the comparisons that exist are encouraging. I would send this to a serious referee, and ask them to focus on the uncertainty claim and the excited-state channels.\n\nWho this is for: plasma modelers and fusion data compilers. It deserves peer review; with a modest revision it would become the standard D2 reference dataset.","headline":"Useful new D2 resonant cross-section dataset with honest v=0 validation; the 15% uncertainty statement is broader than the evidence, but the paper is solid enough for serious review.","tokens_in":15706,"tokens_out":1947,"would_cite":true,"duration_ms":19027,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["34.80.Ht","34.80.Gs"],"model":"deepseek-v4-flash","headline":"Electron impacts on deuterium molecules are computed across three electronic states, giving vibrationally resolved cross sections for vibrational excitation, dissociative attachment, and dissociative excitation via the local complex…","keywords":["electron-D2 collisions","vibrational excitation","dissociative electron attachment","dissociative excitation","local complex potential","cross sections","deuterium plasma","isotope effect"],"falsifier":"Measure the absolute D− yield from D2(X 1Σ+g, v=0) for electron energies from 4 to 15 eV with enough resolution to separate the predicted 10 eV structure from the 14 eV peak, and compare against the DA1+DA2 sums in this paper; a mismatch beyond the claimed ~15% in either peak would show the adopted resonance widths need revision. A second check would be a crossed-beam measurement of the D+D yield from the DE2 channel just above threshold, which would test the paper's claim that dissociation through the b 3Σ+u state dominates.","tokens_in":14767,"feed_emoji":"⚛️","tokens_out":10169,"duration_ms":97013,"temperature":0.7,"pith_summary":"The paper aims to close a data gap: electron-impact cross sections specific to deuterium, rather than scaled hydrogen data, for the resonant processes that control molecular chemistry in fusion divertor plasmas and negative-ion sources. It computes vibrationally resolved cross sections for resonant vibrational excitation, dissociative electron attachment, and dissociative excitation of D2, covering the ground state X 1Σ+g and the first two excited states b 3Σ+u and B 1Σ+u, with electron capture through the three lowest resonances of D−2. The calculation extends the local complex potential method to transitions between different electronic states, not only vibrational levels of the ground state, and it reports agreement with available experiments to about 15%. A sympathetic reading is that the paper establishes that this method can now produce a complete, vibrationally resolved deuterium data set across electronic states, which plasma modelers can adopt directly.","feed_headline":"Deuterium electron-collision data now cover all key breakup channels","feed_subtitle":"Vibration-resolved excitation, attachment, and dissociation data now feed fusion plasma and ion-source models.","key_machinery":"The load-bearing object is the local complex potential (LCP) nuclear equation for the D−2 resonances, $\\left[-\\frac{\\hbar^2}{2\\mu}\\frac{d^2}{dR^2}+V^-_r(R)-\\frac{i}{2}\\Gamma_r(R)-E\\right]\\xi^r_{s,v}(R)=-V^s_r(R)\\chi^s_v(R)$, with $\\xi^r_{s,v}$ the resonance wavefunction and $\\chi^s_v$ the initial molecular vibrational state. The coupling $|V^s_r|^2=\\hbar\\Gamma^s_r/(2\\pi k)$ ties each resonance to each electronic state of D2, and the total width $\\Gamma_r$ is the sum of partial widths $\\Gamma^s_r$. This set of equations is what turns the input potentials and widths into all of the paper's vibrationally resolved cross sections, and it is what makes transitions between electronic states (not just vibrational levels) computable for the first time in this method.","core_discovery":"The central claim is that the local complex potential formalism describes resonant electron–D2 collisions well enough to generate vibrationally resolved cross sections for all three main processes across three electronic states. For each resonance r, the nuclear wavefunction ξr_s,v solves a driven Schrödinger equation with a complex potential V−r(R) − iΓr(R)/2, and the discrete-to-continuum coupling V_s^r is fixed by the partial width Γ_s^r via |V_s^r|^2 = ℏ Γ_s^r/(2πk). All VE, DA, and DE cross sections are squares of overlaps of this wavefunction with vibrational or continuum states of the neutral molecule. The results show that the X 2Σ+u resonance controls low-energy vibrational excitation and attachment, the B 2Σ+g and C 2Σ+g resonances feed dissociation through the b 3Σ+u and B 1Σ+u states, and the C 2Σ+g Rydberg resonance produces the narrow structures seen near 11 eV and the 14 eV attachment peak. The paper validates the set against elastic, vibrational excitation, dissociative attachment, and dissociative excitation measurements, finding agreement around 15%.","pith_inferences":["Because the C 2Σ+g Rydberg resonance drives narrow structures near 11 eV and the 14 eV attachment peak, the same machinery can be applied to T2 and DT, for which comparable data do not yet exist, using credible widths as the main input.","The factor of about 200 isotope suppression in DA1 at v=0 means deuterium divertor models should not rescale hydrogen cross sections by reduced mass; the vibrational dynamics changes the near-threshold resonance shapes, not just energy scales.","A useful next test is to compare Maxwellian rate coefficients computed from these cross sections against measured plasma decay or D− production in low-pressure deuterium discharges, which would validate the data under conditions closer to divertor and ion-source operation."],"forward_implications":["For low vibrational levels of the ground state, mono- and bi-quantic vibrational excitation dominates below about 1 eV, while dissociative excitation through the b 3Σ+u state becomes the leading molecular breakup channel at higher energies.","For medium vibrational levels such as v=10, DA1+DA2 overtakes dissociation and remains the dominant destruction path across the full electron-temperature range shown.","Transitions starting from the B 1Σ+u excited state are dominated by DE5 through the b 3Σ+u final state, while the X←B coupling is weak; VE2 from the ground state to B is inefficient except for narrow spikes near 11 eV.","The full vibrationally resolved cross section set is what plasma codes need to predict the vibrational distribution, dissociation degree, and negative-ion production in detached divertor and negative-ion-source conditions.","The validated agreement with experiment to about 15% supports using these data as a quantitative replacement for the previously available deuterium resonant cross sections."],"supporting_citations":[{"why":"Defines the local complex potential method used to compute all cross sections.","marker":"[35–38]"},{"why":"Supplies the complex potentials for the low-lying X 2Σ+u and B 2Σ+g resonances of D−2.","marker":"[49, 50]"},{"why":"Provides the R-matrix data for the C 2Σ+g Rydberg resonance.","marker":"[51]"},{"why":"Supplies elastic-scattering measurements used to validate the low-energy resonant peak.","marker":"[52]"},{"why":"Provides vibrational-excitation measurements used to validate the VE1 transitions.","marker":"[54]"},{"why":"Provides dissociative-attachment measurements that validate the low-energy DA1 channel.","marker":"[56]"},{"why":"Provides dissociative-attachment measurements used to validate the higher-energy DA structures.","marker":"[57]"},{"why":"Provides isotope-dependent dissociative-attachment measurements used for the DA comparison and the isotope discussion.","marker":"[13]"},{"why":"Supplies adiabatic-nuclei dissociative-excitation cross sections used for the DE comparison.","marker":"[58]"},{"why":"Supplies adiabatic-nuclei dissociative-excitation data used for the DE comparison.","marker":"[59]"}],"fun_headline_variants":["D2 electron breakup: all channels now resolved","All D2 breakup channels from electron impact now available","D2 electron collisions: all breakup routes mapped","Fusion-ready D2 electron cross sections: all breakup channels","Vibrational data for all D2 electron-breakup channels"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the resonance potentials and autoionization widths of D−2, taken from earlier hydrogen-based calculations, are accurate for deuterium; if those widths are wrong, every cross section in the paper inherits the error.","fun_headline_variants_meta":{"raw":{"variants":["D2 electron breakup: all channels now resolved","All D2 breakup channels from electron impact now available","D2 electron collisions: all breakup routes mapped","Fusion-ready D2 electron cross sections: all breakup channels","Vibrational data for all D2 electron-breakup channels"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00131,"raw_usage":{"total_tokens":5320,"prompt_tokens":908,"completion_tokens":4412,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":524,"completion_tokens_details":{"reasoning_tokens":4333}},"tokens_in":524,"tokens_out":4412,"duration_ms":32447,"temperature":1.0,"reasoning_tokens":4333,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T20:42:27.477700+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the absolute D− yield from D2(X 1Σ+g, v=0) for electron energies from 4 to 15 eV with enough resolution to separate the predicted 10 eV structure from the 14 eV peak, and compare against the DA1+DA2 sums in this paper; a mismatch beyond the claimed ~15% in either peak would show the adopted resonance widths need revision. A second check would be a crossed-beam measurement of the D+D yield from the DE2 channel just above threshold, which would test the paper's claim that dissociation through the b 3Σ+u state dominates.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the R-matrix data for the C 2Σ+g Rydberg resonance."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies elastic-scattering measurements used to validate the low-energy resonant peak."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides vibrational-excitation measurements used to validate the VE1 transitions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides dissociative-attachment measurements that validate the low-energy DA1 channel."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides dissociative-attachment measurements used to validate the higher-energy DA structures."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies adiabatic-nuclei dissociative-excitation cross sections used for the DE comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies adiabatic-nuclei dissociative-excitation data used for the DE comparison."}],"review_version":1}