REVIEW 4 major objections 4 minor 38 references
Roles of electron correlation effects for accurate determination of $g_j$ factors of low-lying states of $^{113}$Cd$^+$ and their applications to atomic clock
T0 review · 4 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read Relativistic coupled-cluster calculations determine the $^{113}$Cd$^+$ ground-state $g_j$ factor to be 2.00286(53), making the second-order Zeeman shift in the hyperfine clock transition negligible at the $10^{-16}$ level.
desk verdict Solid RCC calculation of Cd+ g_j factors that fills a real gap for the 113Cd+ clock effort, but the ground-state uncertainty rests on a heuristic 50%-of-triples rule and one Table II parenthesis looks internally inconsistent. 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 central object is the RCCSDpT method: a relativistic coupled-cluster expansion truncated at single and double excitations, with valence triple excitations added perturbatively through operators $T_3$ and $S_{3v}$. The $g_j$ factor is evaluated as the expectation value of the magnetic-moment operator, together with the QED correction operator, using the coupled-cluster wavefunction. The load-bearing pieces are the four triple-excitation terms that enter this expectation value; they give large and strongly cancelling contributions for the low-lying states. The quoted uncertainty is half of the net triple contribution, under the expectation that omitted higher-order clusters shift the result within that range.
What would settle it
Measure the $5s ^2S_{1/2}$ $g_j$ of $^{113}$Cd$^+$ by trapped-ion microwave spectroscopy; a result outside 2.00286(53) would refute the claimed accuracy and the clock conclusion. Alternatively, recompute $g_j$ with full, nonperturbative triple excitations; if the change exceeds half of the perturbative-triples correction, the uncertainty heuristic is falsified.
Extended reading notes
Core claim
The central claim is that electron correlation, not the additional relativistic and QED interactions, controls the accuracy of Cd$^+$ $g_j$ factors, and that a coupled-cluster treatment with perturbative triples determines them with uncertainties around $10^{-4}$. The paper reports $g_j = 2.00286(53)$ for the $5s ^2S_{1/2}$ ground state plus ten excited-state values, and validates the underlying wavefunctions by comparing computed electron affinities with experimental energies, finding sub-percent agreement. Because no measured $g_j$ exists for Cd$^+$, the authors argue from this energy agreement that their $g_j$ values carry comparable accuracy, and they show the ground-state value makes the second-order Zeeman shift in the $^{113}$Cd$^+$ hyperfine clock transition a negligible contributor at the $10^{-16}$ level. They also find triple-excitation contributions decisive for the lower $S$ and $P$ states and assign each uncertainty as half of the net triple contribution.
Load-bearing premise
The accuracy and clock conclusion rest on the unverified premise that omitted higher-order correlation effects are no larger than half of the computed perturbative-triples contribution, since there is no experimental $g_j$ value for Cd$^+$ to check against.
Editorial extensions
If this is right
- The ground-state $g_j = 2.00286(53)$ can enter the $^{113}$Cd$^+$ clock error budget directly, keeping the second-order Zeeman shift below $10^{-16}$ fractional frequency at $B = 10^{-7}$ T.
- The same $g_j$ enables calibration of the applied magnetic field to better than $10^{-10}$ T using the field-sensitive transitions, which keeps the Zeeman-related clock error below $10^{-17}$.
- The ten excited-state $g_j$ values give predictions for future measurements and for evaluating magnetic-field systematics in other Cd$^+$ transitions.
- The calculation identifies perturbative triple excitations as the dominant uncertainty source for the low-lying states, so accuracies beyond $10^{-16}$ would require a full treatment of triple or higher excitations.
Reading between the lines
- If an experimental $g_j$ for $^{113}$Cd$^+$ becomes available, the half-triples uncertainty heuristic can be tested directly: agreement within 0.00053 would validate the error budget, while a larger deviation would require revisiting the clock's Zeeman systematics.
- The same RCCSDpT machinery could supply $g_j$ values for other trapped-ion microwave-clock species, such as Hg$^+$ or Sr$^+$, where measured magnetic-moment data are scarce.
- Because the second-order Zeeman shift depends on $B^2$, the conclusion is sensitive to the assumed field; in situ calibration with the field-sensitive transitions could tolerate slightly larger $B$ while still meeting the $10^{-16}$ target.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports relativistic coupled-cluster (RCCSD and RCCSDpT) calculations of Landé g_j factors for eleven low-lying states of 113Cd+ (5s–7s 2S, 5p–6p 2P, 5d 2D, 4f 2F), including Breit and lower-order QED corrections. Electron affinities are computed at the RCCSD level and compared with NIST values as a wave-function benchmark. Using the ground-state value g_j = 2.00286(53), the authors estimate the second-order Zeeman shift in the |5s, F=0, m_F=0⟩ → |5s, F=1, m_F=0⟩ microwave clock transition and conclude that the 10^-16 fractional-frequency target is not limited by this shift.
Significance. The paper fills a practical gap: no experimental g_j values exist for Cd+, and a reliable g_j is needed for the ongoing 113Cd+ microwave clock program. The method is standard, and the error budget is itemized: DHF/RCCSD/RCCSDpT contributions, Breit, QED, and a table of individual triple-excitation terms. The EA comparison with NIST (0.1–0.5% agreement) provides a nontrivial external check on the orbitals and correlation treatment. If the uncertainty estimate can be made robust, the g_j table and the clock-shift analysis are useful to the atomic-clock and atomic-structure communities. A notable mitigating factor is that the clock conclusion is insensitive to the quoted ground-state uncertainty: even a doubling of δg_j would keep the estimated fractional Zeeman shift below 10^-16.
major comments (4)
- [Sec. IV, Table II] The uncertainty rule of taking half of the perturbative-triples contribution is load-bearing, but the manuscript provides no convergence evidence supporting it. For the ground state, g_D^j moves from 1.999876 (DHF) to 2.001624 (RCCSD) and then by -0.001064 from triples, so the final value is a residue of two large opposite-sign corrections. Because the triples contribution is comparable to (and for 5s larger than) the RCCSD correlation correction, the perturbation series is not in a regime where half-of-triples is self-evidently a conservative bound. I recommend adding a benchmark, such as an RCCSDT (or iterative-triples) calculation for at least the 5s and 5p states, a basis-set extrapolation study, or an application of the same pipeline to an ion with a measured g_j.
- [Table II] For 6p 2P3/2 the quoted final uncertainty (1), i.e., 0.00001, is an order of magnitude smaller than half of the quoted triples contribution of 0.000189 (half = 0.0000945). This contradicts the stated uncertainty rule and appears to be a typo; the entry should be corrected and the effect on any derived quantities checked.
- [Sec. IV, Table I] The caption and text state that the EA uncertainties are estimated from perturbative triple excitations, but no triple corrections to EAs are reported anywhere in the paper; Eq. (19) and the surrounding text describe RCCSD EAs only. If triples contributions to EAs were computed, they should be tabulated; otherwise the EA uncertainties are unsupported. This matters because the EA agreement is the principal external validation of the wave functions.
- [Sec. III, Table III] The list of triple operator combinations in the text and the columns of Table III do not match. The text identifies seven combinations (T2† O T3, S2v† O T3, S2v† O S3v, T2† O S3v, S1v† T2† O S3v, T3† O T3, S3v† O S3v), while Table III tabulates only four, and the status of the first three is not stated. Since the central 'Triples' column in Table II is the sum of these contributions, the manuscript should either tabulate all terms or state explicitly that the omitted terms were evaluated and found negligible.
minor comments (4)
- [Sec. IV, Eq. (29)] The notation ∂(Δν_Zeem(B))/∂g_j A_hf δg_j is ambiguous; the partial derivative should be written with respect to g_j alone, and the text should state clearly that the quoted 1.70×10^-14 δg_j is a fractional-frequency shift relative to A_hf.
- [Fig. 1] The axis labels and some state labels in Fig. 1 appear garbled in the submitted version; please verify that all subplot axes are legible and that the 5d 2D3/2 and 5d 2D5/2 rows are correctly distinguished.
- [Throughout] There are several minor stylistic and typographical issues, including 'systematical shift' in the abstract, 'singly ionized cadmium (Cd+) ion', and the use of 'Howbeit'; these should be corrected in a final pass.
- [References] The NIST database is cited only by URL [31]; please include the specific database version or access date so the comparison values are reproducible.
Circularity Check
No circularity: the g_j values are produced by an ab initio RCCSDpT calculation and benchmarked against external NIST electron affinities; self-citations are methodological, and the uncertainty heuristic is a limitation, not an input-equivalent reduction.
full rationale
I find no circular step in the derivation chain. The paper computes g_j factors from the relativistic coupled-cluster wave functions defined in Eqs. (16)-(23) rather than fitting them to experimental g_j data, and no experimental Cd+ g_j values exist to be used as an input. The energy benchmark in Table I is external: the authors state that electron affinities are compared with "the experimental values listed in the National Institute of Science and Technology (NIST) database" and that this comparison is used "to validate our calculations"; these are independent data sets, not the target g_j values. The self-citations to the authors' earlier works [23] and [30] supply the perturbative-triples implementation and QED effective potentials, but the numerical g_j results are obtained in this paper from the explicit expressions in Eqs. (21)-(23) and Tables II-III; the citations are methodological lineage rather than an imported result that by itself fixes the g_j values. The clock application uses the calculated g_j in the textbook Breit-Rabi expressions, Eqs. (24)-(33), which is an application of the computed value rather than a definitional source of that value. The paper's uncertainty rule, "we assign uncertainties as 50% of the perturbative excitation contributions to the final values of the g_j factors," is a heuristic and is not backed by a convergence study; this is a correctness and robustness concern, not a circularity, because the central g_j values are not constructed from those uncertainties. The manuscript itself flags the key limitation: "There are no experimental values of the g_j factors of the considered states in Cd+ available to compare with our calculations." I also note an internal inconsistency in Table II: for 6p 2P3/2 the triples contribution is 0.000189, half of which is 0.0000945, yet the quoted final uncertainty is (1), much smaller than half of the triples contribution; for the 4f states the quoted uncertainties are also smaller than half of the triples contributions. This is an inconsistency in the stated uncertainty rule, not a circular reduction, and it would be a referee concern about the uncertainty estimate rather than about circularity. No fitted parameter is renamed as a prediction, no uniqueness theorem is imported from the authors' prior work to forbid alternatives, and no known empirical pattern is merely relabeled. The central claim is therefore self-contained against external benchmarks and does not reduce to its inputs.
Assumptions & free parameters
free parameters (1)
- uncertainty fraction for triple contributions =
0.5
assumptions (5)
- domain assumption The Dirac-Coulomb Hamiltonian, augmented by Breit and effective QED potentials, is an adequate representation of 113Cd+
- domain assumption RCCSD with perturbative triple excitations captures the dominant electron correlation effects
- domain assumption The bound-state QED correction to g_j is captured by the 0.001160 spin-scaling term of Eq. (10), with nuclear recoil negligible
- ad hoc to paper Omitted higher excitations produce g_j corrections smaller than half of the perturbative triple contributions
- domain assumption Agreement of electron affinities with NIST at the 0.1-0.5 percent level implies comparable accuracy for g_j factors
Cite this review
Pith. "Pith review of Roles of electron correlation effects for accurate determination of $g_j$ factors of low-lying states of $^{113}$Cd$^+$ and their applications to atomic clock." pith.science (2026). https://pith.science/paper/O7UIG6H3
@misc{pith2026190802007,
author = {Pith},
title = {Pith review of: Roles of electron correlation effects for accurate determination of $g_j$ factors of low-lying states of $^113$Cd$^+$ and their applications to atomic clock},
year = {2026},
howpublished = {\url{https://pith.science/paper/O7UIG6H3}},
note = {Machine review of arXiv:1908.02007}
}
abstract
We investigate roles of electron correlation effects in the determination of $g_j$ factors of the $ns~^2S_{1/2}$ ($n$=5,6,7), $np~^2P_{1/2,3/2}$ ($n$=5,6), $5d~^2D_{3/2,5/2}$, and $4f~^2F_{5/2,7/2}$ states of the singly ionized cadmium (Cd$^+$) ion. Single and double excited configurations along with important valence triple excited configurations through relativistic coupled-cluster (RCC) theory are taken into account for incorporating electron correlation effects in our calculations. We find significant contributions from the triples to the lower $S$ and $P$ states for attaining high accuracy results. The contributions of Breit interaction and lower-order quantum electrodynamics effects, such as vacuum polarization and self-energy corrections, are also estimated using the RCC theory and are quoted explicitly. In addition, we present energies of the aforementioned states from our calculations and compare them with the experimental results to validate $g_j$ values. Using the $g_j$ factor of the ground state, systematical shift due to the Zeeman effect in the microwave clock frequency of the $|5s~^2S_{1/2}, F=0,m_F=0 \rangle \leftrightarrow |5s~^2S_{1/2}, F=1,m_F=0 \rangle$ transition in $^{113}$Cd$^+$ ion has been estimated.
Figures
Reference graph
Works this paper leans on
- [1]
-
[2]
P. Phoonthong, M. Mizuno, K. Kido, and N. Shiga, Appl. Phys. B 117, 673 (2014)
work page 2014
-
[3]
E. A. Burt, L. Yi, B. Tucker, R. Hamell, and R. L. Tjoelker, IEEE Trans. Ultrason., Ferroelectr., Freq. Con- trol. 63, 1013 (2016)
work page 2016
-
[4]
along with values for other variables as defined above. 7 According to this, the calibration of the B value will be affected by the uncertainty in the value of gj. Using our estimated gj value, we anticipate the uncertainty in B would be less than 10 −10 Tesla. This is sufficiently low to maintain the uncertainty in the fractional second-order Zeeman shift wi...
-
[5]
K. Miao, J. W. Zhang, X. L. Sun, S. G. Wang, A. M. Zhang, K. Liang, and L. J. Wang, Opt. Letter 40, 4249 (2015)
work page 2015
-
[6]
J. W. Zhang, Z. B. Wang, S. G. Wang, K. Miao, B. Wang, and L. J. Wang, Phys. Rev. A 86, 022523 (2012)
work page 2012
-
[7]
S. G. Wang, J. W. Zhang, K. Miao, Z. B. Wang, and L. J. Wang, Opt. Express 21, 12434 (2013)
work page 2013
-
[8]
J. W. Zhang, S. G. Wang, K. Miao, Z. B. Wang, and L. J. Wang, Appl. Phys. B 114, 183 (2014)
work page 2014
Show all 38 references
-
[9]
Y. N. Zuo, J. Z. Han, J. W. Zhang, and L. J. Wang, arXiv:1902.10907 [physics.atom-ph]. (2019)
2019 arXiv
-
[10]
C. W. White, W. M. Hughes, G. S. Hayne, and H. G. Robinson, Phys. Rev. A 7, 1178 (1973)
1973
-
[11]
Verd´ u, S
J. Verd´ u, S. Djeki´ c, S. Stahl, T. Valenzuela, M. Vogel, G. Werth, T. Beier, H. J. Kluge, and W. Quint, Phys. Rev. Lett. 92, 093002 (2004)
2004
-
[12]
Sturm, A
S. Sturm, A. Wagner, M. Kretzschmar, W. Quint, G. Werth, and K. Blaum, Phys. Rev. A 87, 030501(R) (2013)
2013
-
[13]
Sturm, F
S. Sturm, F. K¨ ohler, J. Zatorski, A. Wagner, Z. Harman, G. Werth, W. Quint, C. H. Keitel, and K. Blaum, Nature 506, 467 (2014)
2014
-
[14]
K¨ ohler, S
F. K¨ ohler, S. Sturm, A. Kracke, G. Werth, W. Quint, and K. Blaum, J. Phys. B: At. Mol. Pot. Phys. 48, 144032 (2015)
2015
-
[15]
Wagner, S
A. Wagner, S. Sturm, F. K¨ ohler, D. A. Glazov, A. V. Volotka, G. Plunien, W. Quint, G. Werth, V. M. Shabaev, and K. Blaum, Phys. Rev. Lett. 110, 033003 (2013)
2013
-
[16]
A. V. Volotka, D. A. Glazov, V. M. Shabaev, I. I. Tupit- syn, and G. Plunien, Phys. Rev. Lett 112, 253004 (2014)
2014
-
[17]
K¨ ohler, K
F. K¨ ohler, K. Blaum, M. Block, S. Chenmarev, S. Eliseev, D. A. Glazov, M. Goncharov, J. Hou, A. Kracke, D. A. Nesterenko, Y. N. Novikov, W. Quint, E. M. Ramirez, V. M. Shabaev, S. Sturm, A. V. Volotka, and G. Werth, Nature Communications 7, 10246 (2016)
2016
-
[18]
Sturm, M
S. Sturm, M. Vogel, F. K¨ ohler-Langes, W. Quint, K. Blaum, and G. Werth, Atoms 5, 4 (2017)
2017
-
[19]
Veseth, Phys
L. Veseth, Phys. Rev. A 22, 803 (1980)
1980
-
[20]
Veseth, J
L. Veseth, J. Phys. B: At. Mol. Phys. 16, 2891 (1983)
1983
-
[21]
V. A. Dzuba, V. V. Flambaum, P. G. Silvestrov, and 0. P. Sushkov, Phys. Scr. 31, 275 (1985)
1985
-
[22]
G. H. Gossel, V. A. Dzuba, and V. V. Flambaum, Phys. Rev. A 88, 034501 (2013)
2013
-
[23]
Lindroth and A
E. Lindroth and A. Ynnerman, Phys. Rev. A 47, 961 (1993)
1993
-
[24]
B. K. Sahoo and P. Kumar, Phys. Rev. A 96, 012511 (2017)
2017
-
[25]
Y. M. Yu and B. K. Sahoo, Phys. Rev. A 96, 050502(R) (2017)
2017
-
[26]
C. B. Li, Y. M. Yu, and B. K. Sahoo, Phys. Rev. A 97, 022512 (2018)
2018
-
[27]
J. J. Sakurai, Advanced Quantum Mechanics , Addison- Wesley Publishing Company, Virginia, USA, 1967
1967
-
[28]
Czarnecki, U
A. Czarnecki, U. D. Jentschura, K. Pachucki, and V. A. Yerokhin, Can. J. Phys. 84, 453 (2005)
2005
-
[29]
A. J. Akhiezer and V. B. Berestetskii, Quantum Electro- dynamics, Interscience, New York, 1965, Chap. 8, Sec. 50.2
1965
-
[30]
K. T. Cheng and W. J. Childs, Phys. Rev. A 31, 2775 8 (1985)
1985
-
[31]
Y. M. Yu and B. K. Sahoo, Phys. Rev. A 99, 022513 (2019)
2019
-
[32]
http://physics.nist.gov/PhysRefData/ASD/
-
[33]
D. A. Glazov, V. M. Shabaev, I. I. Tupitsyn, A. V. Volotka, V. A. Yerokhin, G. Plunien, and G. Soff, Phys. Rev. A 70, 062104 (2004)
2004
-
[34]
V. M. Shabaev, D. A. Glazov, G. Plunien, and A. V. Volotka, J. Phys. Chem. Ref. Data 44, 031205 (2015)
2015
-
[35]
Breit and I
G. Breit and I. I. Rabi, Phys. Rev. 38, 2082 (1931)
1931
-
[36]
Vanier and C
J. Vanier and C. Audoin, The Quantum Physics of Atomic Frequency Standard Vol. 1, (Adam Hilger, Bris- tol and Philadelphia, 1989)
1989
-
[37]
Itano, J
Wayne M. Itano, J. Res. Natl. Inst. Stand. Technol. 105, 829 (2000)
2000
-
[38]
P. W. Spence and M. N. McDermott, Phys. Lett. 42A, 273 (1972)
1972
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.