REVIEW 3 major objections 2 minor
Precision calculation of the bound-electron $g$ factor in molecular hydrogen ions
T0 review · 3 major / 2 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read The bound-electron g factor in molecular hydrogen ions is computed to a relative accuracy of 4–5 parts in 10^11, an improvement of three orders of magnitude.
desk verdict Abstract promises a three-order precision jump for molecular bound-electron g-factors; the hybrid method is plausible, but the matching between NRQED and the Dirac calculation needs a close look. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing machinery is the two-center Dirac equation, solved with a minmax finite-element method; it supplies all relativistic corrections of order $(Z\alpha)^4$ and higher. The nonrelativistic QED expansion supplies the lower-order and radiative corrections up to order $\alpha^5$. The paper's accuracy claim depends on this hybrid scheme being both complete through order $\alpha^5$ and free of double counting at the matching point between the two parts.
What would settle it
Compare the computed scalar $g$ factor for a specific rovibrational state (e.g., the ground state of $\mathrm{H}_2^+$) with an independent fully relativistic treatment at the claimed $10^{-11}$ level; if the two disagree beyond $5 \times 10^{-11}$, the hybrid scheme has an error. A future Penning-trap measurement of the $g$ factor of a single trapped $\mathrm{H}_2^+$ ion at comparable precision would settle the claim empirically.
Extended reading notes
Core claim
The paper claims that the bound-electron $g$ factor in $\mathrm{H}_2^+$ and $mathrm{HD}^+$ can be computed to a relative accuracy of $4\text{–}5 \times 10^{-11}$ for the scalar component by splitting the calculation into two parts. Contributions through order $\alpha^5$ are handled in a nonrelativistic QED framework, except for relativistic corrections of order $(Z\alpha)^4$ and above, which are obtained from a minmax finite-element solution of the two-center Dirac equation. The result improves on earlier calculations by more than three orders of magnitude and is delivered for a wide range of rovibrational states, making it directly relevant to Penning-trap experiments with single molecular
Load-bearing premise
The central premise is that combining the nonrelativistic QED expansion with the Dirac-solver results covers every contribution through order $\alpha^5$ exactly once, with no gap or double counting between the two methods.
Editorial extensions
If this is right
- At $4\text{–}5 \times 10^{-11}$ accuracy, the scalar $g$ factor can distinguish closely spaced rovibrational states in Penning-trap experiments with single $\mathrm{H}_2^+$ or $\mathrm{HD}^+$ ions, making internal-state identification reliable.
- The results open a new route to precision QED tests, since the two-center Coulomb field of a molecular ion introduces bound-state effects absent in single-electron atoms.
- The hybrid approach demonstrates that combining a nonrelativistic QED expansion with a high-precision Dirac solver works for molecular systems, implying the same strategy can be applied to other small molecules.
- The improvement by three orders of magnitude shifts the limiting uncertainty in molecular $g$-factor calculations to the individual higher-order terms, inviting further refinement of each contribution.
Reading between the lines
- If the claimed accuracy is independently confirmed, the molecular-ion $g$ factor becomes a sensitive probe of how bound-state QED scales with nuclear charge distribution, because the two-center potential makes the $(Z\alpha)^4$ terms considerably harder than in single-electron atoms.
- The same hybrid method could be extended to hyperfine structure or the rotational $g$ factor, turning the scalar $g$ factor into one component of a complete precision Zeeman model for molecular ions.
- A natural cross-check, not described in the paper, is to compute the $g$ factor at the matching boundary between the two methods in both ways and verify that the difference is smaller than the claimed $5 \times 10^{-11}$ error bar.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports a calculation of the bound-electron g factor for rovibrational states of H2+ and HD+. The authors state that relativistic and QED corrections through α^5 are included, using a nonrelativistic QED framework for lower-order terms and a minmax finite-element solution of the two-center Dirac equation for relativistic corrections of order (Zα)^4 and above. They claim a relative accuracy of 4–5 × 10^-11 for the scalar g-factor component, more than three orders of magnitude better than previous calculations, with applications to Penning-trap spectroscopy and QED tests. This review is based solely on the abstract, as the full text is not available.
Significance. If the stated accuracy is correct, this would represent a substantial advance for precision molecular spectroscopy and for tests of bound-state QED in molecular systems. The claimed improvement by more than three orders of magnitude over previous calculations would be notable. However, because the full derivation, numerical convergence study, and comparison with existing data are not available in the abstract, the significance cannot currently be assessed beyond the assertion. The approach of combining NRQED with a high-precision two-center Dirac solver is plausible and potentially powerful, but the central claim depends on a careful matching procedure that is not described.
major comments (3)
- [Abstract (matching procedure)] The central claim of 4–5×10^-11 accuracy depends on the hybrid scheme combining NRQED with the two-center Dirac solution. The abstract does not specify how the Dirac-based 'relativistic g factor' is matched to the NRQED expansion. A full Dirac expectation value contains terms of order (Zα)^2 and (Zα)^3 that are already part of the NRQED calculation; unless these are explicitly subtracted, they are double counted. Conversely, if a subtraction is performed, the remaining all-order terms may contaminate the α^4 coefficient. No matching condition or counterterm subtraction is described, so the claimed order-α^5 completeness is not supported.
- [Abstract (error budget and validation)] The stated relative accuracy of 4–5×10^-11 is asserted without an uncertainty budget. There is no breakdown of numerical convergence (basis size, finite-element mesh, extrapolation), omitted higher-order terms beyond α^5, recoil corrections, finite-nuclear-size effects, or numerical precision of the Dirac solver. Without such a budget, the accuracy claim cannot be verified. The abstract also provides no comparison with known one-electron g factors (e.g., hydrogen-like ions) or with previous molecular-ion calculations, which is needed to establish the claimed improvement by three orders of magnitude.
- [Abstract (completeness of α^5 terms)] The abstract says 'relativistic and QED corrections of orders up to α^5 are taken into account' but does not list which α^5 terms are included. In particular, radiative corrections of order α^5, such as self-energy and vacuum polarization contributions, must be added to the Dirac-based relativistic contribution without double counting any α^5 pieces that may already appear in the expansion of the Dirac expectation value. The hybrid procedure requires a well-defined matching of the QED expansion to the all-order Dirac result; the abstract provides no evidence that such matching is done consistently.
minor comments (2)
- [Abstract (typo)] 'rovibraional' should be 'rovibrational'.
- [Abstract (reference to prior work)] The phrase 'improvement by more than three orders of magnitude over previous calculations' lacks a citation or specific baseline. The reader cannot identify which previous calculation is referenced or how the comparison is made.
Circularity Check
No significant circularity detected in the abstract-only claim.
full rationale
The abstract describes a first-principles calculation combining nonrelativistic QED and a finite-element solution of the two-center Dirac equation. No fitted parameters, no data-driven inputs, and no self-citation chains are visible. The potential concern about double counting between NRQED and Dirac contributions, or about the completeness of the alpha^5 expansion, is a correctness/matching question, not a circularity question. The claimed improvement over previous calculations is a comparison against external results, not a prediction forced by construction. Because the available evidence is limited to the abstract, there is no exhibited reduction of the output to the input, so the circularity score is 0.
Assumptions & free parameters
assumptions (3)
- domain assumption The QED expansion in powers of alpha is valid for the bound-electron g factor in molecular hydrogen ions up to order alpha^5.
- domain assumption The two-center Dirac equation solution using the minmax finite element method accurately captures relativistic corrections of order (Z alpha)^4 and above.
- domain assumption The separation into NRQED for lower orders and a Dirac-based relativistic g factor does not double count or omit contributions through order alpha^5.
Cite this review
Pith. "Pith review of Precision calculation of the bound-electron $g$ factor in molecular hydrogen ions." pith.science (2026). https://pith.science/paper/6MVT3OA5
@misc{pith2026250804242,
author = {Pith},
title = {Pith review of: Precision calculation of the bound-electron $g$ factor in molecular hydrogen ions},
year = {2026},
howpublished = {\url{https://pith.science/paper/6MVT3OA5}},
note = {Machine review of arXiv:2508.04242}
}
abstract
We calculate the bound-electron $g$ factor for a wide range of rovibrational states of the molecular hydrogen ions H$_2^+$ and HD$^+$. Relativistic and QED corrections of orders up to $\alpha^5$ are taken into account. All contributions are calculated in a nonrelativistic QED framework, except for relativistic corrections of order $(Z\alpha)^4$ and above, which are obtained by calculating the relativistic $g$ factor using a precise minmax finite element solution of the two-center Dirac equation. A relative accuracy of $4-5 \times 10^{-11}$ is achieved for the scalar $g$ factor component, which represents an improvement by more than three orders of magnitude over previous calculations. These results are useful for internal state identification and rovibraional spectroscopy of single molecular hydrogen ions in Penning traps, and open a new avenue towards precision tests of QED.
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.