{"id":"b560f80e-e349-48a8-b4e6-9265846f8630","arxiv_id":"2607.07636","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":5,"one_line_summary":"Precision spectroscopy of high-n Rydberg-Stark states yields new experimental values for the hyperfine coupling constant b_F of D₂⁺(v⁺=1, N⁺=0), the spin-rotation constant c_e of H₂⁺(v⁺=1, N⁺=2), and the fundamental vibrational interval of ortho-D₂⁺.","lead":"The paper measures fine and hyperfine structure constants of molecular hydrogen ions (H₂⁺ and D₂⁺) by exploiting long-lived high-n Rydberg-Stark states in weak electric fields. A smart generalist might read it because it validates a new spectroscopic method applicable to any molecular ion and provides benchmark data for tests of fundamental physics.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"Systematic uncertainty from d-series quantum defects is the load-bearing assumption; scatter in c_e fits at n=43 suggests the uncertainty budget may be stretched but not broken, given agreement with independent theory.","rationale":"The reader's verdict of ACCEPT is appropriate. The central claims are supported by three independent cross-checks against ab initio theory, all of which agree within uncertainties. The methodology is transparent, with individual fit results and error budgets provided in the Supplemental Material. The load-bearing concern — systematic uncertainty from d-series quantum defects — is real and acknowledged by the authors, but the agreement with independent theoretical calculations provides strong evidence that the systematic budget is not significantly underestimated. The scatter in individual measurements (particularly c_e at n=43) and the high RMS values of some fits are worth noting, but they do not undermine the central claims given the overall consistency of the results. The concrete test I propose (excluding n=43 H₂ data and performing a jackknife analysis) would provide additional confidence in the robustness of the uncertainty estimates, but the current evidence is sufficient for an ACCEPT verdict. The paper represents a genuine advance: first experimental determination of b_F for D₂⁺(v⁺=1, N⁺=0) and c_e for H₂⁺(v⁺=1, N⁺=2), with precision competitive with or exceeding prior work.","tokens_in":14304,"tokens_out":3418,"duration_ms":325376,"concrete_test":"Recompute the final c_e value for H₂⁺(v⁺=1, N⁺=2) using only the n=34 and n=37 data (excluding the n=43 measurements, which have 7× larger systematic uncertainties and the largest scatter). If the resulting c_e shifts by more than the stated 0.11 MHz uncertainty, the n=43 data are either biasing the result or the systematic uncertainty at n=43 is underestimated. Additionally, perform a jackknife analysis (leave-one-out) across all individual measurements for both b_F and c_e to verify that no single measurement dominates the weighted mean beyond what its stated uncertainty would warrant.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The reader correctly identifies the quantum defects of the d series as the load-bearing systematic. The main text states that black error bar extensions 'reflect systematic uncertainties in the quantum defects of the d series, which are the most strongly coupled to the f series and to the linear Stark manifolds,' but does not describe how these systematics are computed — that methodology resides in the companion paper (Ref. [49], arXiv:2602.17511). This matters because the systematic uncertainties vary dramatically with n: for c_e in H₂, they range from 50 kHz (n=37) to 350 kHz (n=43), a factor of 7. The n=43 H₂ measurements themselves show large scatter (c_e = 40.313(140)(350) and 39.491(386)(350) MHz), with the first value lying ~2σ from the final mean of 39.62(11). For b_F in D₂, the n=34 measurements scatter from 139.842 to 140.059 MHz (range ~217 kHz) against per-point total uncertainties of ~90–145 kHz. The weighted RMS values for several individual fits exceed 2 (up to 4.570 for one D₂ fit and 4.095 for one H₂ fit), indicating that the model does not always reproduce the data at the level of the statistical errors. The final values do agree with independent ab initio theory (b_F: 139.84(5) vs 139.837(10); c_e: 39.62(11) vs 39.5716; vibrational interval: 47,279,980.8(1.9) vs 47,279,981.5898(12) MHz), which provides a strong independent cross-check and suggests the systematic budget is not grossly underestimated. However, the agreement with theory cannot fully substitute for an internally validated systematic uncertainty estimate, because the same MQDT framework underlies both the extraction and the theoretical comparison in some respects. The concern is therefore not that the results are wrong, but that the stated uncertainties rely on a systematic estimation procedure that is not independently verified within this paper.","agreement_with_reader":"agree"},"referee_report":{"model":"glm-5.2","summary":"This manuscript reports precision spectroscopy of high-n molecular Rydberg-Stark states of para-H2 and ortho-D2 in weak electric fields, using the zero-quantum-defect method with frequency-comb calibration. By analyzing the Stark spectra with a combination of multichannel quantum-defect theory (MQDT) and matrix diagonalization, the authors extract new experimental values for the Fermi-contact hyperfine coupling constant b_F = 139.84(5) MHz of D2+(v+=1, N+=0), the spin-rotation coupling constant c_e = 39.62(11) MHz of H2+(v+=1, N+=2), and the fundamental vibrational interval of ortho-D2+ (47,279,980.8(1.9) MHz). All three values agree with independent ab initio theoretical calculations within their stated uncertainties. The paper also presents a physical interpretation of the strikingly different spectral patterns in H2 versus D2, attributing them to the role of the anisotropic charge-quadrupole interaction in the N+=2 core of H2 and its absence in the rotationless N+=0 core of D2.","tokens_in":14729,"tokens_out":1458,"duration_ms":203935,"significance":"The manuscript reports the first experimental determination of the hyperfine coupling constant b_F of D2+(v+=1, N+=0) and the spin-rotation constant c_e of H2+(v+=1, N+=2), both achieving sub-MHz precision. The agreement with independent theory (b_F: 139.84(5) vs 139.837(10) MHz; c_e: 39.62(11) vs 39.5716 MHz; vibrational interval: 47,279,980.8(1.9) vs 47,279,981.5898(12) MHz) provides a strong cross-check on the systematic uncertainty budget. The methodology is general and applicable to other molecular cations, and the experimental infrastructure (frequency-comb calibration, Doppler compensation, detailed error budgets in Tables SIII/SIV) is rigorous. The approach of using Rydberg-Stark states as a probe of ion-core fine and hyperfine structure is a valuable contribution to molecular ion metrology.","major_comments":[{"comment":"The systematic uncertainties from the d-series quantum defects are load-bearing for the extracted values of b_F and c_e, but the methodology for computing these systematics is not described in this manuscript and is deferred to the companion paper (Ref. [49], arXiv:2602.17511). The main text (page 5) states that the black error bar extensions 'reflect systematic uncertainties in the quantum defects of the d series,' but provides no further detail on how they are evaluated. Given that these systematics vary dramatically with n (e.g., from 50 kHz at n=37 to 350 kHz at n=43 for c_e in H2, per Table SII), a brief summary of the methodology here would strengthen the paper's self-containedness and allow the reader to assess the robustness of the uncertainty budget without consulting the companion paper.","section":null},{"comment":"Several individual fits exhibit weighted RMS values well above 1 (up to 4.570 for one D2 fit at n=34, F_z=697.9 mV/cm, and 4.095 for one H2 fit at n=43, F_z=897.7 mV/cm, per Tables SI and SII). This indicates that the model does not always reproduce the data at the level of the statistical errors. The paper does not discuss these elevated RMS values or their potential implications for the final extracted parameters. A brief comment on whether these reflect underestimated statistical errors, residual model deficiencies, or specific spectral features not captured by the MQDT treatment would help the reader evaluate the reliability of the fits.","section":null},{"comment":"The n=43 H2 measurements show notable scatter in c_e: the two fitted values are 40.313(140)(350) and 39.491(386)(350) MHz (Table SII), with the first value lying approximately 2 sigma from the final mean of 39.62(11) MHz. Similarly, the n=34 D2 measurements of b_F range from 139.842 to 140.059 MHz (Table SI), a spread of ~217 kHz against per-point total uncertainties of ~90-145 kHz. While the final weighted means agree with theory, the paper would benefit from a brief discussion of the source of this scatter and why the scatter does not undermine the final quoted values.","section":null}],"minor_comments":[{"comment":"Page 2, first paragraph: 'stuctures' should be 'structures'.","section":null},{"comment":"Page 3, line describing the D2 ion core: the text states 'ortho-D2 (I=0,2)' but later refers to 'I=0' and 'I=2' components separately. The notation could be clarified for readers unfamiliar with the nuclear-spin symmetry species.","section":null},{"comment":"Figure 2 caption: the line widths of 5 MHz are mentioned for the calculated spectra, but the experimental line widths are not stated. Including the experimental line width (or a range) would aid comparison.","section":null},{"comment":"Table I: the theoretical value of c_e is cited as 39.5716 MHz from Ref. [53] (Korobov, Hilico, and Karr, 2006). It would be helpful to note whether more recent theoretical calculations exist, given the 20-year gap.","section":null},{"comment":"The Supplemental Material reference [56] contains a placeholder URL '[url]' that should be replaced with the actual link upon publication.","section":null},{"comment":"Page 5, last paragraph: the GK(2,2)-GK(0,2) interval of 42,639,437.7(9) MHz is mentioned as derived from a 'separate measurement' but no reference or further detail is given. A brief citation or description would be appropriate.","section":null}],"recommendation":"minor_revision","confidential_remarks":"The central experimental results are sound and the agreement with independent theory is compelling. The main concern is the self-containedness of the systematic uncertainty methodology, which is deferred to the companion paper. Since the companion paper (Ref. [49]) appears to be submitted concurrently and is not yet published, the authors should ensure that the present manuscript can stand on its own with respect to the key systematic uncertainty evaluation. The elevated RMS values and scatter in individual fits are not disqualifying given the overall agreement with theory, but they warrant brief discussion. I recommend minor revision to address these presentation and discussion points."},"author_rebuttal":{"model":"glm-5.2","summary":"We thank the referee for the careful reading of our manuscript and for the constructive comments. We agree with all three major comments and will revise the manuscript accordingly.","responses":[{"response":"We agree that a brief summary of how the d-series quantum-defect systematic uncertainties are evaluated would improve the self-containedness of the manuscript. In the revised manuscript, we will add a concise paragraph (approximately 4-5 sentences) on page 5 summarizing the methodology: the systematic uncertainties arise from the limited precision of the d-series quantum defects, which are the most strongly coupled to the f series and the linear Stark manifolds. These uncertainties are propagated through the MQDT calculations by varying the d-series quantum defects within their experimentally determined confidence intervals and recording the resulting shifts in the fitted values of b_F and c_e. The resulting systematic uncertainty scales with n because the sensitivity of the Stark manifold positions to the d-series quantum defects increases with n, as reflected in the n-dependent error bar extensions in Fig. 4. We will also add an explicit cross-reference to Section [X] of Ref. [49] for readers seeking the full derivation. We note that the methodology and its results are fully documented in the companion paper; the addition here is purely for the reader's convenience.","revision_made":"yes","referee_comment":"The systematic uncertainties from the d-series quantum defects are load-bearing for the extracted values of b_F and c_e, but the methodology for computing these systematics is not described in this manuscript and is deferred to the companion paper (Ref. [49]). A brief summary of the methodology here would strengthen the paper's self-containedness."},{"response":"The referee is correct that several individual fits have weighted RMS values significantly above 1, and we agree that the manuscript should comment on this. In the revised manuscript, we will add a brief discussion following the presentation of the fit results. The elevated RMS values arise from a combination of two factors: (1) residual model deficiencies, particularly in the treatment of line intensities (as already noted in the manuscript, saturation and lifetime broadening affect the low-frequency side of the spectra), and (2) the fact that the statistical uncertainties assigned to individual line positions are dominated by the frequency-comb calibration precision and line-center determination, which in some cases are smaller than the residual discrepancies between the MQDT model and the data. Importantly, the fitted parameters (b_F, c_e, ionization energies) are determined from the positions of many lines simultaneously, and the scatter of the individual fit results around the weighted mean (Fig. 4) is consistent with the error bars that include both statistical and systematic contributions. The final quoted uncertainties are derived from the scatter of the individual results, not from the individual fit chi-squares, so the elevated RMS values do not lead to underestimated final uncertainties. We will make this explicit in the revised text.","revision_made":"yes","referee_comment":"Several individual fits exhibit weighted RMS values well above 1 (up to 4.570 for one D2 fit at n=34, F_z=697.9 mV/cm, and 4.095 for one H2 fit at n=43, F_z=897.7 mV/cm). The paper does not discuss these elevated RMS values or their potential implications for the final extracted parameters."},{"response":"We agree that the scatter in the individual results, particularly at n=43 for H2 and n=34 for D2, warrants discussion. In the revised manuscript, we will add a comment addressing this point. The scatter at n=43 in H2 is primarily attributable to the increased systematic uncertainty from the d-series quantum defects at high n (350 kHz, as shown in Table SII), which reflects the growing sensitivity of the Stark manifold structure to the short-range quantum-defect parameters. The two n=43 measurements were taken at different field strengths (995.3 and 897.7 mV/cm), and the different Stark states probed at these fields have different sensitivities to the quantum-defect parameters, leading to the observed scatter. For the n=34 D2 measurements, the spread of ~217 kHz is larger than the per-point statistical uncertainties but is comparable to the systematic uncertainties (55 kHz per point) when combined with the field-dependent variations in the MQDT model accuracy. Crucially, the final values are determined as weighted means over all n values and field strengths (Fig. 4), and the quoted uncertainties are derived from the scatter of these individual results following the procedure described in Ref. [58]. The agreement of the final means with independent ab initio theory provides an additional cross-check. We will add this discussion in the context of the presentation of Fig. 4.","revision_made":"yes","referee_comment":"The n=43 H2 measurements show notable scatter in c_e, and the n=34 D2 measurements of b_F show a spread of ~217 kHz against per-point total uncertainties of ~90-145 kHz. The paper would benefit from a brief discussion of the source of this scatter and why it does not undermine the final quoted values."}],"tokens_in":14736,"tokens_out":1054,"duration_ms":124238,"standing_objections":[]},"desk_editor":{"model":"glm-5.2","letter":"This paper reports the first experimental determinations of the hyperfine coupling constant b_F = 139.84(5) MHz in D₂⁺(v⁺=1, N⁺=0) and the spin-rotation constant c_e = 39.62(11) MHz in H₂⁺(v⁺=1, N⁺=2), plus a vibrational interval in ortho-D₂⁺. These are genuine new results — nobody has measured these before — and the agreement with independent ab initio theory is excellent (b_F: 139.84(5) vs 139.837(10); c_e: 39.62(11) vs 39.5716; vibrational interval: 47,279,980.8(1.9) vs 47,279,981.5898(12) MHz). The experimental side is solid: frequency-comb calibration, Doppler compensation, magnetic shielding, and transparent error budgets in the supplement. The method — using Stark-ℓ-mixing to make autoionizing Rydberg states long-lived enough for precision spectroscopy — is clever and generalizable. The physical picture explaining why H₂ (N⁺=2, quadrupole coupling) gives irregular spectra while D₂ (N⁺=0, clean hyperfine triplets) does not is well argued and convincing. The companion paper (Ref. [49]) provides the MQDT + matrix diagonalization framework, and the calculated spectra genuinely reproduce the measured ones, including intensities. That said, the systematic uncertainty budget has a soft spot worth probing. The dominant systematic comes from d-series quantum defects, and the methodology for estimating it lives in the companion paper rather than here. The n=43 H₂ measurements show real scatter — one point at 40.313(140)(350) MHz sits about 2σ from the final mean — and a couple of individual fits have weighted RMS values above 4 (one D₂ fit at 4.570, one H₂ fit at 4.095). The paper handles this by quoting final uncertainties that encompass the scatter, which is reasonable, but the model doesn't always reproduce data at the statistical error level. The unexplained I-dependent autoionization dynamics in D₂ is an open question the authors flag honestly; it doesn't affect the extracted constants but suggests the MQDT framework has limits. The stress-test concern about systematic estimation not being independently verified within this paper is valid but proportionate: the theory agreement is a strong cross-check, and the final uncertainties look defensible. This is a specialized result aimed at precision spectroscopy and few-body QED practitioners. It deserves a serious referee who can check the MQDT methodology in the companion paper and assess whether the systematic treatment is as robust as claimed. I'd accept for peer review.","headline":"New experimental values for b_F and c_e in molecular hydrogen ions, with systematic uncertainty budget that deserves scrutiny but holds up","tokens_in":15378,"tokens_out":660,"would_cite":true,"duration_ms":171455,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"glm-5.2","headline":"Rydberg-Stark spectroscopy measures molecular ion hyperfine structure","keywords":[],"falsifier":"If the fitted b_F and c_e values disagreed with independent theoretical calculations or with future direct measurements of the same ion-core constants beyond the stated uncertainties, the Rydberg-Stark spectroscopic extraction method would be called into question.","tokens_in":14467,"feed_emoji":"⚛️","tokens_out":776,"duration_ms":135293,"temperature":0.7,"pith_summary":"The paper demonstrates that high-n molecular Rydberg-Stark states, when studied via precision spectroscopy in weak electric fields, serve as a direct window into the fine and hyperfine structure of molecular ions that are otherwise extremely difficult to access. The core mechanism is that an external electric field mixes the orbital states of the distant Rydberg electron, decoupling it from the ion core and making the normally short-lived autoionizing states long-lived and sharp enough for precision measurement. The spectrum of these Stark states then mirrors the internal structure of the ion core. The authors apply this method to two contrasting systems — para-H₂ with a rotating (N⁺=2) ion core and ortho-D₂ with a non-rotating (N⁺=0) ion core — and show that the strikingly different spectral patterns arise from the presence or absence of anisotropic charge-quadrupole interactions between the rotating ion core and the Rydberg electron. By fitting the spectra with a combined multichannel quantum-defect theory and matrix diagonalization approach, they extract new experimental values for the Fermi-contact hyperfine coupling constant b_F of D₂⁺(v⁺=1, N⁺=0) [139.84(5) MHz], the spin-rotation coupling constant c_e of H₂⁺(v⁺=1, N⁺=2) [39.62(11) MHz], and the fundamental vibrational interval of ortho-D₂⁺ [47,279,980.8(1.9) MHz], all in agreement with ab initio theory.","feed_headline":"Rydberg-Stark spectra reveal molecular ion hyperfine structure at MHz precision","feed_subtitle":"Electric-field-decoupled Rydberg states act as a probe of H₂⁺ and D₂⁺ fine and hyperfine couplings, yielding new experimental values that ab","key_machinery":"High-n Rydberg-Stark states in weak electric fields; multichannel quantum-defect theory (MQDT) combined with matrix diagonalization; the distinction between anisotropic charge-quadrupole interactions (present for N⁺≥1 rotating cores) and isotropic Fermi-contact hyperfine interactions (dominant for N⁺=0 non-rotating cores).","core_discovery":"The central discovery is that the fine and hyperfine structure of molecular hydrogen ions can be precisely determined from the spectra of high-n Rydberg-Stark states, because the Stark effect makes these states long-lived and their spectral structure directly encodes the ion-core interactions. The method works for both rotating cores (where anisotropic charge-quadrupole coupling complicates the spectrum) and non-rotating cores (where the Rydberg electron motion is nearly separable from the core hyperfine dynamics), and yields coupling constants competitive with or exceeding the precision of prior experimental determinations.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Stark-decoupled Rydberg states probe molecular ion hyperfine structure","Rydberg-Stark spectra yield MHz-precision hyperfine constants for H₂⁺ and D₂⁺","Long-lived Rydberg-Stark states encode fine and hyperfine structure of H₂⁺ and D₂⁺","Electric-field-decoupled Rydberg states reveal ion-core couplings at MHz precision","Rydberg-Stark spectroscopy measures D₂⁺ vibrational interval and hyperfine constants"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The extraction of coupling constants relies on the quantum defects and MQDT parameters used in the matrix diagonalization model being complete and accurate, particularly for the d-series quantum defects that couple most strongly to the f-series and the linear Stark manifolds; if these parameters carry unaccounted errors, the fitted b_F and c_e values would shift.","fun_headline_variants_meta":{"raw":{"variants":["Stark-decoupled Rydberg states probe molecular ion hyperfine structure","Rydberg-Stark spectra yield MHz-precision hyperfine constants for H₂⁺ and D₂⁺","Long-lived Rydberg-Stark states encode fine and hyperfine structure of H₂⁺ and D₂⁺","Electric-field-decoupled Rydberg states reveal ion-core couplings at MHz precision","Rydberg-Stark spectroscopy measures D₂⁺ vibrational interval and hyperfine constants"]},"model":"glm-5.2","effort":"high","cost_usd":0.0,"raw_usage":{"total_tokens":958,"prompt_tokens":851,"completion_tokens":107,"prompt_tokens_details":null},"tokens_in":851,"tokens_out":107,"duration_ms":45390,"temperature":1.0,"reasoning_tokens":null,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-09T04:13:56.031192+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"If the fitted b_F and c_e values disagreed with independent theoretical calculations or with future direct measurements of the same ion-core constants beyond the stated uncertainties, the Rydberg-Stark spectroscopic extraction method would be called into question.","supporting_citations":[],"review_version":1}