REVIEW 3 major objections 4 minor 1 cited by
Finding the ultra-narrow $^3\!P_2 \rightarrow \, ^3\!P_0$ electric quadrupole transition in Ni$^{12+}$ ion for an optical clock
T0 review · 3 major / 4 minor · reviewed 2026-08-08 · deepseek-v4-flash
Pith's one-line read The energy of the ultra-narrow 3P2 → 3P0 electric quadrupole transition in Ni12+ is predicted at 20081(10) cm^-1 and measured at 20078.984(10) cm^-1, an agreement of 2 cm^-1.
desk verdict First optical excitation of the Ni12+ clock transition, with a 10 cm^-1 ab initio prediction that landed 2 cm^-1 from the measured value; a strong letter with one modest caveat about the full-CI residual estimate. 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 tool is a pure CI calculation for all 16 electrons of the Ni12+ ground configuration, broken into additive contributions from excitations of the outer six electrons, inner ten electrons, extra reference configurations, and an estimated residual full-CI term. The residual is estimated through a linear relation between a configuration's relative weight and its energy contribution, extrapolated from computed configurations. A hybrid CI+all-order method with a coupled-cluster effective Hamiltonian provides an independent cross-check, and a neural-network-based configuration selection reduced memory and time needs.
What would settle it
If another independent computation of the 3P0 energy, or a future more precise measurement, deviates from 20081 $cm^{-1}$ by more than 10 $cm^{-1}$, the claimed uncertainty and the extrapolation of the weight-energy relation would be falsified. A direct check would be an independent CI calculation with a different code or basis that does not use the same linear extrapolation.
Extended reading notes
Core claim
The central claim is that a fully converged 16-electron configuration-interaction calculation can predict the energy of the forbidden 3P2 → 3P0 clock transition in Ni12+ with an uncertainty of only 10 $cm^{-1}$, sufficient to locate the transition with rapid laser scans. The prediction, 20081(10) $cm^{-1}$, differs from the measured value 20078.984(10) $cm^{-1}$ by 2 $cm^{-1}$. The authors argue this validates their method for other complex atomic systems and enables a Ni12+ optical clock.
Load-bearing premise
The estimate of the residual full-CI contribution assumes that the linear weight–energy relation found for computed configurations holds for the far more numerous configurations that cannot be computed directly, with the 3P0 correction being small (-3 $cm^{-1}$) but not independently verified.
Editorial extensions
If this is right
- The Ni12+ clock transition is now accessible for quantum-logic spectroscopy and the development of a high-precision optical clock.
- The demonstrated accuracy of the 16-electron CI method supports its application to other highly charged ions where clock transitions have not yet been located.
- The separation of CI into additive contributions and the weight-energy estimate extends the practical size of CI calculations beyond the current limits of exact diagonalization.
Reading between the lines
- The neural-network configuration-selection algorithm, developed for this calculation, could be applied to other large CI problems where the full expansion is intractable.
- The measured 2 cm^-1 discrepancy suggests the remaining uncertainty may be dominated by the estimated full-CI residual or QED corrections; a systematic study of similar transitions could refine these estimates.
- The same fast-scanning search strategy, combined with precise predictions, could be used to find other forbidden transitions in HCIs that were previously considered too weak to search for.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper predicts the energy of the strongly forbidden 3s2 3p4 3P2 → 3s2 3p4 3P0 electric quadrupole transition in Ni12+ as 20081(10) cm−1, using two independent many-body approaches: a hybrid CI+all-order method and a pure 16-electron CI calculation. The prediction was then used to search for and find the transition experimentally at 20078.984(10) cm−1, only 2 cm−1 from the theoretical value. The authors claim that this level of agreement, with a stated 0.05% uncertainty for a 16-electron system, is unprecedented and qualifies their method for future atomic-structure calculations.
Significance. If the uncertainty estimate holds, this is a landmark result: it demonstrates that a full 16-electron configuration-interaction calculation can predict a state energy to about 0.05% accuracy, and the experimental confirmation in a few hours is a powerful validation. The paper also contributes an open-source parallel CI package, a neural-network-based configuration selection tool, and a two-method cross-check that agree to 5 cm−1 for the target level. The main weakness is the unquantified extrapolation used to estimate the residual full-CI contribution, which is central to the claimed 10 cm−1 uncertainty budget.
major comments (3)
- [Theory: pure CI method (Table II, column 'Estimate full CI')] The estimate of the contribution of all omitted configurations, labeled 'Estimate full CI', is obtained from a linear relation between configuration weight and energy contribution, but no uncertainty is assigned to this extrapolation. The 10 cm−1 uncertainty budget appears to rely on the difference between the two methods and the agreement for the 3P1 level, and does not explicitly include the error in this extrapolation. The 1S0 level in Table II shows a discrepancy of −66 cm−1 from experiment, with an 'Estimate full CI' correction of only −24 cm−1, demonstrating that the method's accuracy is strongly level-dependent. Therefore, the small −3 cm−1 correction for 3P0 does not by itself guarantee that the extrapolation error is small. Please provide a quantitative bound for the extrapolation error, for example by testing the linear relation on additional levels or by varying the slope within a plausible range, or explicitly state that the 10 cm−1 budget excludes this source and justify why it is negligible.
- [Theory: pure CI method, paragraph on 'Estimate full CI'] The description 'a relative weight of 0.01% leading to a correction of about −11 cm−1' is ambiguous: it is unclear whether this refers to the contribution of a single configuration or to the integrated contribution of all configurations with weights below a threshold. The manuscript does not provide the actual fitting parameters, the range of weights over which the linear relation was validated, or a plot of the data supporting the 'linear dependence (similar for all four levels)' claim. Without this information, the 'Estimate full CI' column in Table II is not reproducible and its reliability cannot be independently assessed.
- [Conclusion] The conclusion states that the work demonstrates 'full convergence of a 16-electron CI computation', but this overstates the status of the calculation: the largest part of the configuration space is not computed directly but estimated via the extrapolated linear relation. A more precise wording would be 'convergence with respect to explicitly included configurations and an estimated residual contribution', which would be consistent with the actual procedure described in the text.
minor comments (4)
- [Throughout] Several formulas are typeset incorrectly in the text as provided, e.g., '3s23p4' instead of '3s^2 3p^4' and 'Ni12+' instead of 'Ni$^{12+}$'. These should be corrected in the final version.
- [References] Reference [12] is not a standard citation but a footnote explaining the treatment of k and l partial waves. It should be renumbered as a footnote or integrated into the main text.
- [Table I and Table II] The column groups such as '17 spd fg', 'l >6', 'Extra conf.', and 'Estimate full CI' are not fully defined in the captions. Please add a sentence explaining the meaning of the grouped columns, e.g., that '17 spd fg' refers to orbitals with n up to 17 and partial waves s, p, d, f, g, and that the contributions are incremental from the previous column.
- [Figure 1] The caption 'Illustration of the basis set upscale' is vague. The figure should describe what the blocks represent (e.g., principal quantum number ranges and partial waves) and how it relates to the convergence procedure.
Circularity Check
No significant circularity: the 3P2–3P0 energy is computed ab initio and the experimental value is used only as a posteriori validation, not as an input.
full rationale
The paper's central prediction (20081(10) cm^-1 for the 3P2–3P0 E2 transition) is obtained from two independent many-body treatments, CI+all-order and 16-electron CI. No parameter is fitted to the 3P0 energy or to the measured clock transition. The experimental 3P1, 1D2, and 1S0 levels are used as benchmarks in Tables I and II, not as constraints in the calculation. The self-citations to the CI package pCI [13,14], the CI+all-order formalism [10], the recurrent basis [11,18], and the neural-network selection [19] are method citations; they do not smuggle the predicted transition energy into the calculation. The 'Estimate full CI' procedure in the pure CI section is an empirical extrapolation of a weight–energy relation tested on computed configurations; it is explicitly acknowledged as an estimate ('impossible to compute directly'), is applied uniformly to the low-lying levels, and contributes only -3 cm^-1 to the 3P0 energy. Its reliability concerns the uncertainty budget, but it is not circular because the correction is not constructed from the target 3P0 energy. The agreement with the subsequently measured value (20078.984(10) cm^-1) is an external check, not an input. Thus no step in the derivation reduces by construction to its own input.
Assumptions & free parameters
free parameters (1)
- Full-CI weight-to-energy conversion factor =
-11 cm^-1 per 0.01% relative weight
assumptions (4)
- domain assumption The CI model space with single and double excitations from the chosen reference configurations is sufficient to converge the low-lying state energies to a few cm^-1.
- ad hoc to paper The relationship between configuration relative weight and energy contribution is linear and can be extrapolated to all omitted configurations.
- domain assumption QED corrections computed in Ref. [15] are accurate to 10% or better for these levels.
- domain assumption The frozen-core approximation in the CI+all-order method, with core-valence correlations treated by coupled-cluster, is valid for this system.
Cite this review
Pith. "Pith review of Finding the ultra-narrow $^3\!P_2 \rightarrow \, ^3\!P_0$ electric quadrupole transition in Ni$^{12+}$ ion for an optical clock." pith.science (2026). https://pith.science/paper/PO242C34
@misc{pith2026250205386,
author = {Pith},
title = {Pith review of: Finding the ultra-narrow $^3\!P_2 \rightarrow \, ^3\!P_0$ electric quadrupole transition in Ni$^12+$ ion for an optical clock},
year = {2026},
howpublished = {\url{https://pith.science/paper/PO242C34}},
note = {Machine review of arXiv:2502.05386}
}
abstract
The Ni$^{12+}$ ion features an electronic transition with a natural width of only 8 mHz, allowing for a highly stable optical clock. We predict that the energy of this strongly forbidden $3s^2 3p^4\, ^3\!P_2 \rightarrow 3s^2 3p^4 \, ^3\!P_0$ electric quadrupole transition is 20081(10) cm$^{-1}$. For this, we use both a hybrid approach combining configuration interaction (CI) with coupled-cluster (CC) method and a pure CI calculation for the complete 16-electron system, ensuring convergence. The resulting very small theoretical uncertainty of only 0.05\% allowed us to find the transition experimentally in a few hours, yielding an energy of 20078.984(10) cm$^{-1}$. This level of agreement for a 16-electron system is unprecedented and qualifies our method for future calculations of many other complex atomic systems. While paving the way for a high-precision optical clock based on Ni$^{12+}$, our theory and code development will also enable better predictions for other highly charged ions and other complex atomic systems.
Figures
Forward citations
Cited by 1 Pith paper
-
Coulomb Crystallization of Highly Charged Ni^12+ Ions in a Linear Paul Trap
Ni12+ ions were sympathetically cooled with laser-cooled Be+ in a room-temperature Paul trap to form a two-species Coulomb crystal at the 100 mK level.
Reference graph
Works this paper leans on
-
[1]
M. G. Kozlov, M. S. Safronova, J. R. Crespo L´ opez- Urrutia, and P. O. Schmidt, Rev. Mod. Phys. 90, 045005 (2018)
work page 2018
-
[2]
M. S. Safronova, D. Budker, D. DeMille, D. F. J. Kimball, A. Derevianko, and C. W. Clark, Rev. Mod. Phys. 90, 025008 (2018)
2018
- [3]
-
[4]
S. A. King, L. J. Spieß, P. Micke, A. Wilzewski, T. Leopold, E. Benkler, R. Lange, N. Huntemann, A. Surzhykov, V. A. Yerokhin, J. R. Crespo L´ opez- Urrutia, and P. O. Schmidt, Nature 611, 43 (2022)
work page 2022
-
[5]
A. D. Ludlow, M. M. Boyd, J. Ye, E. Peik, and P. O. Schmidt, Rev. Mod. Phys. 87, 637 (2015)
2015
- [6]
-
[7]
S. Chen, Z. Zhou, J. Li, T. Zhang, C. Li, T. Shi, Y. Huang, K. Gao, and H. Guan, Phys. Rev. Res. 6, 013030 (2024)
work page 2024
-
[8]
S. Chen, L. J. Spieß, A. Wilzewski, M. Wehrheim, K. Di- etze, I. Vybornyi, K. Hammerer, J. R. C. L´ opez-Urrutia, and P. O. Schmidt, Phys. Rev. Appl. 22, 054059 (2024)
work page 2024
Show all 27 references
-
[9]
Kramida, Yu
A. Kramida, Yu. Ralchenko, J. Reader, and and NIST 6 ASD Team, NIST Atomic Spectra Database (ver. 5.9), [Online]. Available: https://physics.nist.gov/asd [Mon Jun 20 2022]. National Institute of Standards and Technology, Gaithersburg, MD. (2022)
2022
-
[10]
M. S. Safronova, M. G. Kozlov, W. R. Johnson, and D. Jiang, Phys. Rev. A 80, 012516 (2009)
2009
-
[11]
M. G. Kozlov, S. G. Porsev, and V. V. Flambaum, J. Phys. B 29, 689 (1996)
1996
-
[12]
We have verified with h and i partial waves that, as expected, the total contri- bution of a partial wave does not depend on the choice of the basis
We find the calculations of k and l partial waves to be intractable with such a large basis. We have verified with h and i partial waves that, as expected, the total contri- bution of a partial wave does not depend on the choice of the basis. We use far more compact basis with...
-
[13]
Cheung, M
C. Cheung, M. S. Safronova, and S. G. Porsev, Symmetry 13 (2021)
2021
-
[14]
Cheung, M
C. Cheung, M. G. Kozlov, S. G. Porsev, M. S. Safronova, I. I. Tupitsyn, and A. I. Bondarev, Computer Physics Communications 308, 109463 (2025)
2025
-
[15]
I. I. Tupitsyn, M. G. Kozlov, M. S. Safronova, V. M. Shabaev, and V. A. Dzuba, Phys. Rev. Lett. 117, 253001 (2016)
2016
-
[16]
M. G. Kozlov, M. S. Safronova, S. G. Porsev, and I. I. Tupitsyn, Phys. Rev. A 94, 032512 (2016)
2016
-
[17]
C. Shah, M. Togawa, M. Botz, J. Danisch, J. J. Goes, S. Bernitt, M. Maxton, K. K¨ obnick, J. Buck, J. Selt- mann, M. Hoesch, M. F. Gu, F. S. Porter, T. Pfeifer, M. A. Leutenegger, C. Cheung, M. S. Safronova, and J. R. Crespo L´ opez-Urrutia, Astrophys. J. 969, 52 (2024)
2024
-
[18]
M. G. Kozlov, S. G. Porsev, M. S. Safronova, and I. I. Tupitsyn, Comp. Phys. Comm. 195, 199 (2015)
2015
-
[19]
Bilous, C
P. Bilous, C. Cheung, and M. Safronova, Phys. Rev. A 110, 042818 (2024)
2024
-
[20]
Micke, S
P. Micke, S. K¨ uhn, L. Buchauer, J. R. Harries, T. M. B¨ ucking, K. Blaum, A. Cieluch, A. Egl, D. Hollain, S. Kraemer, T. Pfeifer, P. O. Schmidt, R. X. Sch¨ ussler, C. Schweiger, T. St¨ ohlker, S. Sturm, R. N. Wolf, S. Bernitt, and J. R. Crespo L´ opez-Urrutia, Review of Scie...
2018
-
[21]
Leopold, S
T. Leopold, S. A. King, P. Micke, A. Bautista-Salvador, J. C. Heip, C. Ospelkaus, J. R. Crespo L´ opez-Urrutia, and P. O. Schmidt, Review of Scientific Instruments 90, 073201 (2019)
2019
-
[22]
P. O. Schmidt, T. Rosenband, C. Langer, W. M. Itano, J. C. Bergquist, and D. J. Wineland, Science 309, 749 (2005)
2005
-
[23]
Micke, T
P. Micke, T. Leopold, S. A. King, E. Benkler, L. J. Spieß, L. Schm¨ oger, M. Schwarz, J. R. Crespo L´ opez-Urrutia, and P. O. Schmidt, Nature 578, 60 (2020)
2020
-
[24]
F. Wolf, Y. Wan, J. C. Heip, F. Gebert, C. Shi, and P. O. Schmidt, Nature 530, 457 (2016)
2016
-
[25]
In total, we covered in six hours 100 GHz within the uncertainty range until we found the clock transition
allowed us to detect the excitation from the ground state with nearly 100% fidelity. In total, we covered in six hours 100 GHz within the uncertainty range until we found the clock transition. Figure 2 shows an example of an off-resonant clock-laser scan leading to high motion...
2020
-
[27]
D. B. Hume, T. Rosenband, and D. J. Wineland, Phys. Rev. Lett. 99, 120502 (2007)
2007
-
[35]
17 spd fg
The Breit and Coulomb interactions were included on the same footing for constructing the basis set, which is large enough for the coupled-cluster computations to accurately generate the effective Hamiltonian. We apply several approximations of increasing accu- racy to constru...
2000
Reviewed August 8, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.