REVIEW 1 major objections 4 minor 21 references
Para- and diamagnetic contributions to magnetic shielding constants of relativistic hydrogenlike atoms in some low-lying discrete energy eigenstates
T0 review · 1 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read This paper presents tabulated values of the relative paramagnetic and diamagnetic contributions to the magnetic shielding constant of relativistic hydrogenlike atoms for all nuclear charges Z from 1 to 137.
desk verdict Table 8's 3p3/2 (mu=+/-1/2) columns are corrupted for Z=70-137: the sum sigma_d+sigma_p exceeds 1, the sigma_p column is copied from Table 4, and 68 rows must be recomputed before the paper's central data can be trusted. 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 Gordon decomposition of the Dirac current, which separates the magnetic interaction into a convection (diamagnetic-like) and a spin (paramagnetic-like) piece, giving $\sigma = \sigma_d + \sigma_p$. The paper works with the explicit closed forms (1.2)–(1.6), expressed through the radial quantum number $n$, the Dirac quantum number $\kappa$, the magnetic quantum number $\mu$, and $\gamma_\kappa = \sqrt{\kappa^2 - (\alpha Z)^2}$. These formulas, imported from the author's earlier derivation that used a Sturmian expansion of the Dirac–Coulomb Green function, are what the tables evaluate.
What would settle it
Recompute the ratios by an independent method, such as solving the Dirac equation with a magnetic perturbation numerically for a few representative cases (for example $1s_{1/2}$ with $Z = 80$ and $2p_{1/2}$ with $Z = 118$), and compare against the tabulated values at the stated precision. Any disagreement beyond rounding would falsify the formulas or the evaluation.
Extended reading notes
Core claim
The paper's central claim is that for any discrete energy eigenstate of a relativistic hydrogenlike atom, the magnetic shielding constant splits into a diamagnetic part $\sigma_d$ and a paramagnetic part $\sigma_p$ whose relative sizes are fixed by the three quantum numbers $(n, \kappa, \mu)$ through the closed expressions in Eqs. (1.2)–(1.6). The tabulated ratios $\sigma_d/\sigma$ and $\sigma_p/\sigma$ are the numerical content of those expressions. The data exhibit regular trends: for states with $\kappa < 0$ at maximal $\mu$, $\sigma_d$ dominates; for $\kappa > 0$, $\sigma_p$ dominates; for $s$-states there is a crossover $Z_c$ (67 for $1s$, 79 for $2s$, 86 for $3s$) after which the paramagnetic term is larger. For $p_{1/2}$ states, the paramagnetic share has a minimum at a high $Z_c$ (99 for $2p_{1/2}$, 105 for $3p_{1/2}$). The same formulas are used to compile total shielding constants for three CODATA values of the fine-structure constant.
Load-bearing premise
The correctness of every table rests on the imported closed-form expressions for $\sigma_d$ and $\sigma_p$ being valid for all discrete states and for all $Z$ up to 137; the paper does not re-derive them here.
Editorial extensions
If this is right
- The tables give benchmarks for approximate relativistic quantum-chemical calculations of magnetic shielding in heavy ions.
- The crossover charges $Z_c$ provide a simple diagnostic for when paramagnetic shielding overtakes diamagnetic shielding in $s$-states.
- Table 17 quantifies how sensitive the total shielding constant is to the fine-structure constant across CODATA 2014, 2018, and 2024 values.
- Because the formulas hold for arbitrary discrete states, the same evaluation can be extended to higher excited states and to arbitrary $\mu$ values without further derivation.
Reading between the lines
- The one-electron results are the natural zero-order reference for many-electron heavy atoms, so these ratios could be used to isolate electron-correlation and finite-nuclear-size effects on shielding.
- The observed $Z_c$ pattern across $n$ suggests a nearly linear relation that an empirical formula could capture; one could test $Z_c(n)$ against a larger set of $n$ values.
- For states with $|\kappa|=1$, the restriction $Z < \alpha^{-1}\sqrt{3/2} \approx 118.67$ means the tables stop before the Dirac critical charge; extending to the supercritical regime would require a different boundary condition.
- The sharp alpha-dependence seen near high $Z$ (e.g., $1s$ at $Z=118$) implies these ratios could in principle be used to constrain $\alpha$ if measured in highly charged ions.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper tabulates relative diamagnetic and paramagnetic contributions to the magnetic shielding constant of relativistic hydrogenlike atoms. The author evaluates closed-form expressions from her earlier work (Eqs. (1.2)-(1.6)) for the ground state and the first two excited-state manifolds, with Z from 1 to 137 and alpha^-1 = 137.035999177, and presents the ratios sigma_d/sigma and sigma_p/sigma in Tables 1-14. It also compares absolute sigma_d and sigma_p values with Pyper and Zhang (Tables 15-16) and lists total sigma for selected ions under three CODATA values of alpha (Table 17).
Significance. The paper fills a gap by providing a systematic tabulation of the separation of sigma into para- and diamagnetic parts for many states and Z values, using parameter-free analytical formulas and agreeing with published values where comparisons exist. The formulas are stated explicitly, and no parameters are fitted; the comparisons in Tables 15-16 provide independent numerical checks for a subset of cases. However, because the contribution of the paper is exclusively numerical data, internal consistency of every table is essential. The error in Table 8 described below undermines a substantial block of the data until corrected.
major comments (1)
- [Table 8 (3p3/2, mu = +/- 1/2), rows Z = 70-137] These rows violate the identity sigma_d + sigma_p = sigma stated in Eq. (1.1). For example, at Z = 70 the entries sigma_d/sigma = 1.056872965322 and sigma_p/sigma = 1.138700966968 sum to 2.195573932290, not 1; at Z = 137 they sum to 3.543724464168. Every valid row must satisfy sigma_d/sigma + sigma_p/sigma = 1, and all other tables in the paper do. The sigma_p values in this block coincide with the 2p3/2 (mu = +/- 1/2) values in Table 4 (e.g., 1.138700966968 at Z = 70), rather than with a 3p3/2 evaluation. Since the abstract claims Tables 1-14 cover 1 <= Z <= 137, the 68 affected rows must be recomputed from Eqs. (1.2)-(1.4); this is a load-bearing data error, not a typographical issue.
minor comments (4)
- [Abstract and Introduction] The text contains several typographical errors: "diamagentic" in the abstract, "analitycal" in the Introduction, and "theese" in the Introduction; these should be corrected.
- [Tables 4, 8, 13, 14] Several numerical entries contain stray spacing in the printed mantissas (for example, "1.13870096696 8(0)" in Table 8 and "2.010634177 934(0)" in Table 4); the final typeset version should ensure each number is printed as one continuous token.
- [Section 2] The author states that the results were obtained with the formulas from Eqs. (1.1)-(1.6), but the paper does not identify which numerical precision or rounding convention was used in the tables; a brief statement of the estimated numerical accuracy would help readers judge the reliability of the last digits.
- [References [18] and [20]] References [18] and [20] are given as web addresses without full bibliographic details; they should be completed with the standard citation information for the CODATA reports.
Circularity Check
No significant circularity: the tables are direct evaluations of explicitly stated formulas, with external comparisons, not fits renamed as predictions.
full rationale
The paper's deliverable is numerical tables of sigma_d/sigma and sigma_p/sigma for hydrogenlike ions. The formulas evaluated (Eqs. 1.2-1.6) are written out in full in the paper; they are not fitted to the tabulated data. The only free parameters are the CODATA fine-structure constant and the scanned nuclear charge Z. Eq. (1.1) defines sigma = sigma_d + sigma_p and is used as a consistency check, not as an input whose output is then relabeled. The author's earlier paper [8] is cited as the source of the closed-form expressions, but since the expressions themselves are stated and the numerical evaluation is transparent, the self-citation is not load-bearing. Tables 15 and 16 compare against independent results of Pyper and Zhang for a subset of states, providing an external benchmark. A reader could in principle reproduce every entry from Eqs. (1.2)-(1.6) without accepting Ref. [8]; therefore no step reduces to its own inputs. A data-integrity issue such as Table 8 rows failing the identity sigma_d/sigma + sigma_p/sigma = 1 would be a correctness defect, not a circularity.
Assumptions & free parameters
assumptions (4)
- ad hoc to paper Correctness of the previously derived formulas in Eqs. (1.2)-(1.6)
- domain assumption Validity of the Sturmian expansion of the generalized Dirac-Coulomb Green function
- domain assumption Gordon decomposition into diamagnetic and paramagnetic currents
- domain assumption Pointlike, spinless, motionless nucleus of charge Ze
Cite this review
Pith. "Pith review of Para- and diamagnetic contributions to magnetic shielding constants of relativistic hydrogenlike atoms in some low-lying discrete energy eigenstates." pith.science (2026). https://pith.science/paper/APSTZETL
@misc{pith2026241201364,
author = {Pith},
title = {Pith review of: Para- and diamagnetic contributions to magnetic shielding constants of relativistic hydrogenlike atoms in some low-lying discrete energy eigenstates},
year = {2026},
howpublished = {\url{https://pith.science/paper/APSTZETL}},
note = {Machine review of arXiv:2412.01364}
}
abstract
We present tabulated data for numerical calculations of relative para- and diamagentic contributions to the magnetic shielding constant ($\sigma$) of the Dirac one-electron atoms with a pointlike, spinless and motionless nuclei of charge $Ze$. Exploiting the analytical formulas for the diamagnetic ($\sigma_{d}$) and paramagnetic ($\sigma_{p}$) terms of $\sigma$, valid for an arbitrary discrete energy state, recently derived by us with the aid of the Gordon decomposition technique, we have found the numerical values of $\sigma_{d}/\sigma$ and $\sigma_{p}/\sigma$ for the ground state and for the first two sets of excited states (i.e.: $2s_{1/2}$, $2p_{1/2}$, $2p_{3/2}$, $3s_{1/2}$, $3p_{1/2}$, $3p_{3/2}$, $3d_{3/2}$, and $3d_{5/2}$) of the relativistic hydrogenic ions with the nuclear charge number from the range $1 \leqslant Z \leqslant 137$. The comparisons of our results with those reported by other authors for some atomic states are also presented. We also compile here the numerical values of the total magnetic shielding constants for the ground state $1s_{1/2}$ and for each state belonging to the first set of excited states of selected hydrogenlike ions, obtained with the use of three different values of the fine-structure constant, i.e.: $\alpha^{-1}=137.035 \: 999 \: 139$ (from CODATA 2014), $\alpha^{-1}=137.035 \: 999 \: 084$ (from CODATA 2018) and $\alpha^{-1}=137.035 \: 999 \: 177$ (from CODATA 2024).
Reference graph
Works this paper leans on
- [8]
-
[1]
R. Szmytkowski, J. Phys. B 30 (1997) 825 [erratum: J. Phys. B 3 0 (1997) 2747]
work page 1997
-
[2]
R. Szmytkowski and P. Stefa´ nska, e-print arXiv:1102.1811
-
[3]
P. Stefa´ nska and R. Szmytkowski, Int. J. Quantum Chem. 112 (2012) 1363
work page 2012
- [4]
- [5]
- [6]
- [7]
Show all 21 references
-
[9]
Stefa´ nska, At
P. Stefa´ nska, At. Data Nucl. Data Tables 108 (2016) 193
2016
-
[10]
Stefa´ nska, At
P. Stefa´ nska, At. Data Nucl. Data Tables 113 (2017) 316
2017
-
[11]
Stefa´ nska, At
P. Stefa´ nska, At. Data Nucl. Data Tables 120 (2018) 352
2018
-
[12]
Stefa´ nska, Atoms 12 (2024) 54
P. Stefa´ nska, Atoms 12 (2024) 54
2024
-
[13]
Stefa´ nska, At
P. Stefa´ nska, At. Data Nucl. Data Tables 135-136 (2020) 10 1360
2020
-
[14]
J. H. Van Vleck, The theory of electric and magnetic susceptibilit ies, Oxford University Press, London, 1932
1932
-
[15]
Gordon, Z
W. Gordon, Z. Phys. 50 (1928) 630
1928
-
[16]
N. C. Pyper, Mol. Phys. 97 (1999) 381
1999
-
[17]
N. C. Pyper and Z. C. Zhang, Mol. Phys. 97 (1999) 391
1999
-
[18]
P. J. Mohr, E. Tiesinga, D. B. Newell and B. N. Taylor, Codata In ternation- ally Recommended 2022 Values of the Fundamental Physical Consta nts. Available at http://physics.nist.gov/constants
2022
-
[19]
E. R. Cohen and B. N. Taylor, J. Phys. Chem. Ref. Data 17 (198 8) 1795
-
[20]
Tiesinga, P
E. Tiesinga, P. J. Mohr, D. B. Newell and B. N. Taylor, The 2018 C ODATA Recommended Values of the Fundamental Physical Constants, 20 20. Available at http://physics.nist.gov/constants
2018
-
[21]
P. J. Mohr, D. B. Newell and B. N. Taylor, Rev. Mod. Phys. 88 (2 016) 035009. 21
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.