arxiv: 2511.15346 · v1 · submitted 2025-11-19 · ⚛️ physics.atom-ph · nucl-ex· nucl-th· physics.chem-ph

Comprehensive Assessment of Th³⁺ Properties for Nuclear Clock and Fundamental Physics Applications

A. Chakraborty , B. K. Sahoo This is my paper

Pith reviewed 2026-05-17 20:53 UTC · model grok-4.3

classification ⚛️ physics.atom-ph nucl-exnucl-thphysics.chem-ph

keywords thorium-3+ ionnuclear clockisotope shiftsnuclear charge radiihyperfine structurerelativistic coupled-clusteratomic polarizabilitiesfundamental physics

0 comments p. Extension

The pith

Relativistic calculations give accurate nuclear charge radii and moments for thorium ions to support nuclear clocks.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper seeks to compute many atomic properties of the Th^{3+} ion with a high-level quantum method so that these values can help develop a nuclear clock using the 229Th isotope and allow tests of fundamental physics. A reader would care because better atomic data reduces timing errors in the clock and sharpens searches for changes in physical constants or new forces. The authors combine their computed isotope shift factors with measured data to obtain precise differences in nuclear size between 232Th and 229Th, and between the ground and isomeric states of 229Th. They also extract the nuclear magnetic dipole and electric quadrupole moments for both states of 229Th by fitting hyperfine constants to their results. The calculations further supply electric dipole polarizabilities and hyperfine-induced quadrupole moments needed to estimate systematic errors in a working nuclear clock.

Core claim

By using the relativistic coupled-cluster method with singles, doubles, and triples excitations, the authors calculate a broad set of atomic properties for the Th^{3+} ion. Combining the resulting isotope shift parameters with experimental data yields highly accurate differential nuclear charge radii for 232,229Th and 229m,229Th. Nuclear magnetic dipole and electric quadrupole moments are obtained for both the ground and isomeric states of 229Th. Electric dipole polarizabilities and hyperfine-induced quadrupole moments are evaluated to assess clock systematics. The work finds unexpectedly large contributions from higher-order relativistic effects and excitations involving higher angular-mom-

What carries the argument

Relativistic coupled-cluster framework with singles, doubles, and triples excitations, which generates accurate electronic wave functions and properties for the heavy thorium ion by systematically including electron correlations and relativistic corrections.

If this is right

Differential nuclear charge radii for 232,229Th and 229m,229Th become available at the stated high accuracy.
Nuclear magnetic dipole and electric quadrupole moments are fixed for the ground and isomeric states of 229Th.
Systematic uncertainties in a 229Th^{3+}-based nuclear clock can be evaluated using the computed polarizabilities and hyperfine-induced moments.
Atomic energy predictions for these ions require inclusion of higher-order relativistic effects and excitations to high angular momentum orbitals.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The derived nuclear properties could tighten limits on possible time variation of fundamental constants when thorium clocks are compared with other frequency standards.
The same computational approach could be applied to other heavy ions to identify additional candidates for nuclear clocks.
The noted importance of higher angular momentum orbitals indicates that future calculations should employ even larger orbital sets to reduce remaining uncertainty.
Linking these atomic results to nuclear structure models may improve predictions of the 229Th isomer energy and lifetime.

Load-bearing premise

The coupled-cluster calculations, when merged with experimental hyperfine constants, produce nuclear radii and moments whose uncertainties are smaller than those from direct measurement alone.

What would settle it

An independent experimental measurement of the differential nuclear charge radius between 229Th and 232Th that lies outside the uncertainty interval reported from the combined theory-experiment approach.

Figures

Figures reproduced from arXiv: 2511.15346 by A. Chakraborty, B. K. Sahoo.

**Figure 2.** Figure 2: FIG. 2. Comparison of (a) [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

read the original abstract

By employing singles, doubles, and triples excitations within the relativistic coupled-cluster framework, we perform comprehensive calculations of a wide range of atomic properties for the Th$^{3+}$ ion. These properties are essential for advancing nuclear clock technology and probing fundamental physics. Combining our isotope shift parameters with experimental data, we estimate highly accurate values of the differential nuclear charge radii for $^{232,229}$Th and $^{229m,229}$Th. Additionally, we determine the nuclear magnetic dipole and electric quadrupole moments for both the ground and isomeric states of $^{229}$Th by combining measured hyperfine structure constants with our theoretical calculations. Our precise evaluations of electric dipole polarizabilities and hyperfine-induced quadrupole moments are critical for assessing systematic uncertainties in $^{229}$Th$^{3+}$-based nuclear clock. Notably, we observe unexpectedly significant contributions from higher-order relativistic effects and excitations involving orbitals with higher angular momentum, which markedly influence the energies of the ground state and its fine-structure partner. These results highlight the substantial challenges in achieving highly accurate predictions for these properties.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

Th^{3+} RCCSDT calculations are standard but the extracted nuclear radii and moments rest on unshown convergence and benchmarks, so the accuracy claims need checking.

read the letter

The main point is that this paper applies relativistic coupled-cluster with singles, doubles, and triples to Th^{3+} and pulls nuclear charge radii and moments out by dividing experimental isotope shifts and hyperfine constants by the computed electronic factors. They also flag that higher-order relativistic corrections and high-angular-momentum excitations turn out larger than expected for the ground-state energies. That combination is the practical deliverable for the nuclear-clock effort. The work does a reasonable job of covering a wide list of properties—polarizabilities, hyperfine-induced moments, and the isotope-shift parameters—that directly feed into systematic budgets for a ^{229}Th clock. The method is established for open-shell ions, and the authors are upfront that the system is difficult. The soft spot is the accuracy story. The abstract and the stress-test note give no basis-set extrapolation numbers, no quadruples test, and no reproduction of known hyperfine constants in lighter homologues. Without those, it is not clear that the theoretical factors carry smaller uncertainty than the final nuclear values they produce. If the full manuscript has those checks in tables or appendices, the claim strengthens; if not, the propagated error on the radii and moments stays unquantified. This paper is for the small group working on ^{229}Th nuclear clocks and related fundamental-physics tests. A reader who needs order-of-magnitude estimates for planning can take the numbers as a starting point. Anyone who wants to publish a new limit on constant variation or dark-matter coupling using these values should first verify the error budget themselves. I would send it to peer review. The calculations are non-trivial, the application is timely, and referees can ask for the missing convergence data in one round.

Referee Report

2 major / 2 minor

Summary. The manuscript reports relativistic coupled-cluster calculations (RCCSDT) of atomic properties for Th^{3+}, including energies, hyperfine structure constants, isotope-shift parameters, electric dipole polarizabilities, and hyperfine-induced quadrupole moments. These theoretical electronic factors are combined with experimental hyperfine and isotope-shift data to extract differential nuclear charge radii for ^{232,229}Th and ^{229m,229}Th as well as nuclear magnetic dipole and electric quadrupole moments for the ground and isomeric states of ^{229}Th. The work emphasizes large higher-order relativistic and high-angular-momentum contributions to the low-lying states and discusses implications for ^{229}Th-based nuclear clocks.

Significance. If the computed electronic factors can be shown to carry uncertainties smaller than the target nuclear precision, the extracted radii and moments would constitute useful benchmarks for nuclear-structure models and systematic-error budgets in nuclear-clock proposals. The explicit recognition of substantial higher-order relativistic and high-l effects is a constructive observation for the community working on heavy open-shell ions.

major comments (2)

[Abstract and § on isotope-shift calculations] Abstract and results section on isotope shifts: the claim that the extracted differential charge radii are 'highly accurate' with uncertainties smaller than experiment presupposes that the theoretical isotope-shift parameters (field-shift and mass-shift factors) are known to sub-percent level, yet no basis-set extrapolation, no test of quadruple excitations, and no reproduction of known hyperfine constants in lighter homologues (e.g., Ac^{2+} or Pa^{4+}) are presented to quantify the residual theoretical error.
[Hyperfine structure and nuclear moments] Section on hyperfine constants and nuclear moments: the semi-empirical extraction of μ and Q for ^{229}Th ground and isomer states inherits the full theoretical uncertainty of the computed electronic hyperfine parameters; without an explicit error budget or comparison against independent all-order or QED-corrected methods, the assertion that the final nuclear moments are more accurate than existing literature values cannot be verified.

minor comments (2)

[Notation] Notation for the differential mean-square charge radii should be standardized (e.g., consistent use of δ⟨r²⟩ versus Δ⟨r²⟩) and defined at first appearance.
[Energy results] Table of computed energies would benefit from an additional column listing the size of the triple-excitation contribution to allow readers to judge the importance of the T3 terms.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the constructive and detailed report. We address the two major comments point by point below. Where the comments identify gaps in uncertainty quantification, we have revised the manuscript to incorporate additional analysis and moderated language; we note one area where full independent benchmarks remain outside present computational reach.

read point-by-point responses

Referee: Abstract and results section on isotope shifts: the claim that the extracted differential charge radii are 'highly accurate' with uncertainties smaller than experiment presupposes that the theoretical isotope-shift parameters (field-shift and mass-shift factors) are known to sub-percent level, yet no basis-set extrapolation, no test of quadruple excitations, and no reproduction of known hyperfine constants in lighter homologues (e.g., Ac^{2+} or Pa^{4+}) are presented to quantify the residual theoretical error.

Authors: We accept that the original wording overstated the precision without a dedicated uncertainty section. In the revised manuscript we have added an explicit discussion of basis-set convergence (showing changes when the basis is extended by one cardinal number) and a perturbative estimate of quadruple-excitation contributions to the field- and mass-shift factors. We have also inserted a short paragraph referencing our earlier RCCSDT benchmarks on Ac^{2+} and Pa^{4+} that used the same code and active-space strategy. The abstract has been revised to replace 'highly accurate' with 'with quantified theoretical uncertainties smaller than the experimental isotope-shift precision'. revision: yes
Referee: Section on hyperfine constants and nuclear moments: the semi-empirical extraction of μ and Q for ^{229}Th ground and isomer states inherits the full theoretical uncertainty of the computed electronic hyperfine parameters; without an explicit error budget or comparison against independent all-order or QED-corrected methods, the assertion that the final nuclear moments are more accurate than existing literature values cannot be verified.

Authors: We agree that an explicit error budget is required. The revised version now contains a dedicated paragraph that propagates the estimated 1–2 % theoretical uncertainty in the hyperfine electronic factors (obtained from basis-set and active-space variations) into the extracted nuclear moments, and compares the final values with the most recent literature determinations. Full all-order or QED-corrected calculations for this open-shell ion are not yet available in the literature; we therefore qualify the claim of improved accuracy by stating that the new moments are consistent with prior values within the combined experimental and theoretical uncertainties. revision: partial

standing simulated objections not resolved

Direct numerical comparison of our hyperfine constants against independent QED-corrected all-order calculations, which do not currently exist for Th^{3+}.

Circularity Check

0 steps flagged

No significant circularity; derivation uses independent ab initio computation to extract nuclear quantities from experiment

full rationale

The paper computes electronic factors (isotope shift parameters, hyperfine constants, polarizabilities) via relativistic coupled-cluster with singles, doubles and triples excitations applied to the Th^{3+} ion. These factors are then combined with separate experimental isotope shifts and hyperfine constants to extract differential nuclear charge radii and nuclear moments. This is a standard semi-empirical extraction in which the theoretical electronic quantities are obtained from first-principles many-body calculations without being fitted to the nuclear data under study. No equation reduces an output to an input by construction, no self-citation supplies a uniqueness theorem or ansatz, and no fitted parameter is relabeled as a prediction. The derivation chain therefore remains self-contained against external benchmarks and does not exhibit circularity.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central results rest on the validity of the relativistic coupled-cluster ansatz for an open-shell actinide ion and on the assumption that experimental hyperfine constants can be inverted using the computed electronic factors without significant higher-order nuclear effects.

axioms (1)

domain assumption Relativistic coupled-cluster method with singles, doubles, and triples excitations captures the dominant correlation and relativistic effects for Th³⁺ valence properties
Invoked throughout the abstract as the computational framework

pith-pipeline@v0.9.0 · 5495 in / 1389 out tokens · 22917 ms · 2026-05-17T20:53:05.153632+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

By employing singles, doubles, and triples excitations within the relativistic coupled-cluster framework, we perform comprehensive calculations of a wide range of atomic properties for the Th^{3+} ion... Combining our isotope shift parameters with experimental data, we estimate highly accurate values of the differential nuclear charge radii
IndisputableMonolith/Foundation/DimensionForcing.lean alexander_duality_circle_linking unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Notably, we observe unexpectedly significant contributions from higher-order relativistic effects and excitations involving orbitals with higher angular momentum

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Towards better nuclear charge radii
nucl-ex 2026-04 unverdicted novelty 2.0

An effort is described to produce more precise and transparent recommended values for nuclear charge radii through integrated experimental and theoretical approaches.

Reference graph

Works this paper leans on

33 extracted references · 33 canonical work pages · cited by 1 Pith paper

[1]

Kraemer, J

S. Kraemer, J. Moens, M. Athanasakis-Kaklamanakis, S. Bara, K. Beeks, P. Chhetri, K. Chrysalidis, A. Claessens, T. E. Cocolios, J. G. M. Correia et al., Nature (London) 617, 706 (2023)

work page 2023
[2]

Peik and C

E. Peik and C. Tamm, Europhys. Lett.61, 181 (2003)

work page 2003
[3]

E. Peik, T. Schumm, M. S. Safronova, A. P´ alffy, J. Weit- enberg, and P. G. Thirolf, Quantum Sci. Technol.6, 034002 (2021)

work page 2021
[4]

C. J. Campbell, A. G. Radnaev, A. Kuzmich, V. A. Dzuba, V. V. Flambaum, and A. Derevianko, Phys. Rev. Lett.108, 120802 (2012)

work page 2012
[5]

V. V. Flambaum, Phys. Rev. Lett.97, 092502 (2006)

work page 2006
[6]

M. S. Safronova, D. Budker, D. DeMille, D. F. J. Kimball, A. Derevianko, and C. W. Clark, Rev. Mod. Phys.90, 025008 (2018)

work page 2018
[7]

C. J. Campbell, A. G. Radnaev, and A. Kuzmich, Phys. Rev. Lett.106, 223001 (2011)

work page 2011
[8]

Zitzer, J

G. Zitzer, J. Tiedau, Ch. E. D¨ ullmann, M. V. Okhapkin, and E. Peik, Phys. Rev. A111, L050802 (2025)

work page 2025
[9]

M. S. Safronova, U. I. Safronova, A. G. Radnaev, C. J. Campbell, and A. Kuzmich, Phys. Rev. A88, 060501(R) (2013)

work page 2013
[10]

J. C. Berengut, V. A. Dzuba, V. V. Flambaum, and S. G. Porsev, Phys. Rev. Lett.102, 210801 (2009)

work page 2009
[11]

V. A. Dzuba and V. V. Flambaum, Phys. Rev. Lett.131, 263002 (2023)

work page 2023
[12]

China Phys

Ran Si, Chaofan Shi, Nan Xue, Xiangjin Kong, Chongyang Chen, Bingsheng Tu, Yu-Gang Ma, Sci. China Phys. Mech. Astron.68, 272011 (2025)

work page 2025
[13]

Safronova, S

M.S. Safronova, S. G. Porsev, M.G. Kozlov, J. Thielking, M.V. Okhapkin, P. Glowacki, D.M. Meier, and E. Peik, Phys. Rev. Lett.121, 213001 (2018)

work page 2018
[14]

W. H. King,Isotope shifts in atomic spectra, Springer Science & Business Media, (2013)

work page 2013
[15]

B. K. Sahoo, S. Blundell, A. V. Oleynichenko, R. F. G. Ruiz, L. V. Skripnikov, and B. Ohayon, J. Phys. B58, 042001 (2025)

work page 2025
[16]

Chakraborty and B

A. Chakraborty and B. K. Sahoo, Phys. Rev. A111, 042807 (2025)

work page 2025
[17]

Chakraborty and B

A. Chakraborty and B. K. Sahoo, Phys. Rev. A111, 062812 (2025)

work page 2025
[18]

Chakraborty, B

Vaibhav Katyal, A. Chakraborty, B. K. Sahoo, Ben Ohayon, Chien-Yeah Seng, Mikhail Gorchtein, and John Behr, Phys. Rev. A111, 042813 (2025)

work page 2025
[19]

S. G. Porsev, M. S. Safronova, and M. G. Kozlov, Phys. Rev. Lett.127, 253001 (2021)

work page 2021
[20]

Kramida, Yu

A. Kramida, Yu. Ralchenko, J. Reader, and NIST ASD Team, NIST Atomic Spectra Database, version 5.8, https://physics.nist.gov/asd (2020)

work page 2020
[21]

Blaise and J

J. Blaise and J. Wyart,Selected con- stants, energy levels, and atomic spectra of actinides, https://www.lac.universite-paris- saclay.fr/Data/Database/

work page
[22]

Thielking, M

J. Thielking, M. V. Okhapkin, P. Glowacki, D. M. Meier, L. von der Wense, B. Seiferle, C. E. Dullmann, P. G. Thirolf, and E. Peik, Nature556, 321 (2018)

work page 2018
[23]

Zhang, Tian Ooi, Jacob S

C. Zhang, Tian Ooi, Jacob S. Higgins, Jack F. Doyle, Lars von der Wense, Kjeld Beeks, Adrian Leitner, Georgy A. Kazakov, Peng Li, Peter G. Thirolf, Thorsten Schumm, and Jun Ye, Nature633, 63 (2024)

work page 2024
[24]

R. R. Chasman, I. Ahmad, A. M. Friedman, and J. R. Erskine, Rev. Mod. Phys.49, 833 (1977)

work page 1977
[25]

Minkov and A

N. Minkov and A. P´ alffy, Phys. Rev. Lett.122, 162502 (2019)

work page 2019
[26]

Gerstenkorn, P

S. Gerstenkorn, P. Luc, J. Verges, D. W. Englekemeir, J. E. Gindler, and F. S. Tomkins, J. Phys. (Paris) 35, 483 (1974)

work page 1974
[27]

C. E. Bemis, F. K. McGowan, J. L. C. Ford, Jr., W. T. Milner, R. L. Robinson, P. H. Stelson, G. A. Leander, and C. W. Reich, Phys. Scr.38, 657 (1988)

work page 1988
[28]

Li, H.-X

F.-C. Li, H.-X. Qiao, Y.-B. Tang, and T.-Y. Shi, Phys. Rev. A104, 062808 (2021)

work page 2021
[29]

Schwartz, Phys

C. Schwartz, Phys. Rev.97, 380 (1955)

work page 1955
[30]

B. K. Sahoo, Phys. Rev. A111, L060801 (2025)

work page 2025
[31]

Derevianko, Phys

A. Derevianko, Phys. Rev. A93, 012503 (2016)

work page 2016
[32]

Shi-Cheng Yu, Wen-Ting Gan, Xia Hua, Xin Tong and Cheng-Bin Li, Phys. Rev. A109, 063115 (2024)

work page 2024
[33]

M. S. Safronova and U. I. Safronova, Phys. Rev.A87, 062509 (2013)

work page 2013