REVIEW 3 major objections 5 minor 52 references
Sturmian basis set for the Dirac equation with finite nuclear size: Application to polarizability, Zeeman and hyperfine splitting, and vacuum polarization
T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read A modified Sturmian basis with a free scale parameter reproduces the Dirac spectrum and QED properties of hydrogen-like ions with finite nuclear size, including the all-order vacuum polarization density.
desk verdict A practical Sturmian basis with solid benchmarks on standard observables, but the headline vacuum-polarization result needs a convergence study and an independent reference overlay before I'd trust it. 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 object is the modified relativistic Coulomb-Sturmian basis: $\pi_L(r)=S_{n,l_L}(2\lambda r)$ and $\pi_S(r)=S_{n,l_S}(2\lambda r)$, where $S_{nl}$ is the Coulomb-Sturmian radial function (a Laguerre polynomial times $r^{l+1}e^{-\lambda r}$) and $l_{L,S}=|\kappa\pm 1/2|-1/2$. Because these functions reproduce the Dirac radial asymptotics at the origin, they can represent wave functions inside the finite-size nucleus. Rayleigh-Ritz reduction of the radial Dirac equation produces a symmetric generalized eigenvalue problem whose matrix elements are evaluated analytically from Laguerre integrals; the free parameter $\lambda$ lets the basis concentrate at the nuclear scale ($\lambda\sim 1/r_n$ for vacuum polarization) without changing the matrix structure.
What would settle it
Compute the $|\kappa|\ge 4$ partial-wave contributions to the Wichmann-Kroll vacuum polarization density with increasing basis size $n$: if the basis is complete in the nuclear interior these contributions must converge to the known Green's-function results, whereas a $\lambda$-plateau that shifts with $n$, or a failure to converge at short distances, would show the interior is not fully resolved.
Extended reading notes
Core claim
The central claim is that Coulomb-Sturmian functions can be adapted to a finite nuclear charge distribution by replacing the single orbital index $l$ with two indices, $l_L$ and $l_S$, for the large and small components, and by treating the Sturmian scale $\lambda$ as a free parameter instead of tying it to the ground-state energy. The resulting basis has the correct $r\to 0$ asymptotics of the Dirac radial functions for a non-singular potential, and so is designed to resolve the nuclear interior. With this basis and the Rayleigh-Ritz method, the authors reproduce the finite-nuclear-size corrections to binding energies, static dipole polarizability, first-order hyperfine splitting, the bound-electron $g$ factor, and nuclear magnetic shielding in hydrogen-like ions, at the level of dual-kinetic-balance B-spline calculations and analytical formulas. For the Wichmann-Kroll vacuum polarization density, the basis performs on par with, and for the $|\kappa|=2$ partial wave slightly better than, a machine-precision Gaussian basis, while avoiding the Gaussian basis's growing linear dependence.
Load-bearing premise
The load-bearing premise is that the modified Sturmian functions, chosen only to have the right shape near $r=0$ and at infinity, still form a complete enough set to represent Dirac wave functions inside a finite-size nucleus; the paper checks this numerically at each step but does not prove it.
Editorial extensions
If this is right
- Finite-nuclear-size corrections to hydrogenic binding energies, polarizabilities, hyperfine structure, and $g$ factors can be obtained in one basis, without kinetic-balance construction, matching DKB B-spline and analytical benchmarks.
- The Wichmann-Kroll vacuum polarization density, previously a challenge for finite-basis methods, is accessible with this basis at both short and large distances from the nucleus.
- The absence of linear dependence lets the basis be enlarged to improve accuracy without the precision collapse that limits Gaussian-basis calculations.
- The construction depends only on spherical symmetry and the correct $r\to 0$ asymptotics, so the same machinery extends in principle to arbitrary binding potentials.
- For each quantity the free-$\lambda$ plateau supplies an internal convergence check and an uncertainty estimate, as shown for the hyperfine and shielding factors.
Reading between the lines
- A natural extension, not pursued in the paper, is to apply the same $\lambda$-tuning to one- and two-loop QED diagrams such as electron self-energy, where the short-distance part of the electron propagator is decisive and linear independence of the basis could remove a practical obstacle.
- The $\lambda$-plateau behavior could be developed into an automatic uncertainty procedure: scan $\lambda$, locate the plateau or extremum, and use its width as a numerical error band for any computed matrix element.
- If completeness in the nuclear interior is confirmed by stronger tests or a proof, the same recipe of fixing large- and small-component asymptotics separately may transfer to relativistic molecular calculations, where exponential-type bases and linear independence are both attractive.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a relativistic Coulomb-Sturmian-like basis set for solving the radial Dirac equation with a finite-size nucleus. The basis functions are obtained from ordinary Coulomb Sturmians by replacing the orbital index with the asymptotic exponents l_L and l_S appropriate for the large and small components in the finite-nucleus interior. The authors apply the basis in a Rayleigh-Ritz procedure to compute binding energies, static dipole polarizability, first-order hyperfine splitting, the bound-electron g factor, and the nuclear magnetic shielding factor, comparing with analytical formulas and DKB B-spline results. As the main new application, they compute the all-order (Wichmann-Kroll) vacuum polarization charge density for hydrogen-like uranium and compare it with their earlier Gaussian-basis results. The paper argues that the Sturmian basis avoids linear-dependence problems and is well suited for finite-nucleus problems, with the scale parameter λ treated as a free parameter.
Significance. If the vacuum polarization results are fully validated, the proposed basis set would be a useful tool for relativistic atomic-structure and QED calculations, offering an alternative to B-splines and Gaussian bases while being linearly independent and providing analytic matrix elements. The benchmark comparisons for energies, polarizability, hyperfine splitting, and g factor are credible and agree with independent results to the stated precision. The resolution-of-identity test and the absence of spurious states are positive indicators. However, the strongest new claim—that the basis reproduces the all-order vacuum polarization density—is not yet backed by the same level of convergence evidence as the other observables; the paper lacks an n-convergence study, error bars, and an overlay of an independent reference for that quantity. The overall approach is promising, but the central VP claim needs strengthening before the paper can be accepted.
major comments (3)
- [Sec. IV.E, Figs. 3–8] The vacuum polarization density, which is the paper's headline application, is presented without a convergence study in the basis size n. The text states that the quality of results reaches a plateau as n increases, but no table or figure analogous to Table VI documents this for ρ(r). Since the VP density is dominated by the nuclear interior where the modified basis differs most from the exact Sturm-Liouville functions, a plateau-like insensitivity to λ could mask systematic truncation error. The authors should provide a quantitative n- and λ-convergence analysis for the VP density, with either tabulated values at selected radii or a figure showing the dependence on n.
- [Sec. IV.E, Figs. 3–8] The text claims that the VP curves 'can be compared' with the Green's-function results of Ref. [33] and the finite-basis results of Ref. [30], but no such overlay is shown. The only comparison presented is with the authors' own Gaussian-basis results from Ref. [31]. Adding an overlay of an independent calculation (e.g., the Green's function integration from Ref. [33]) would directly test the VP claim and address the concern that both basis-set calculations might share a common systematic error. Without this, the statement that the CS basis provides 'slightly better' results is not well quantified.
- [Sec. IV.D, Figs. 1–2 and Table VI] The selection of the free parameter λ is done manually for each observable and each ion, based on identifying a plateau or a maximum in the λ-dependence. The paper states that the plateau width implies an uncertainty, but this is a heuristic criterion. For the VP calculation, only two λ values (λ=94 and λ=54) are shown, without demonstrating that they lie in a converged plateau for the VP density. The authors should either provide a more formal convergence criterion (e.g., monitoring the plateau as a function of n and λ) or at least present a convergence table for the VP quantity analogous to Table VI, so that the reader can judge whether the chosen λ values are justified.
minor comments (5)
- [Abstract and throughout] The manuscript would benefit from a careful proofreading pass; for example, 'The tin hydrogen-like ion' (Sec. IV.A) should be 'The hydrogen-like tin ion', and 'Strumian' appears in the conclusion instead of 'Sturmian'.
- [Eq. (52)] In the sentence following Eq. (52), 'where where for brevity' contains a duplicated 'where'.
- [Eq. (41)] In the definition of the shell-nucleus potential, the notation 'r > r n' has an unwanted space; it should be 'r > r_n'.
- [Sec. IV.E] The parameter choice for the VP calculation is stated as 'λ ~ 1/r_n', but the two values actually used (λ=94 and λ=54) are not derived from this estimate for Z=92 with r_n = 5.8507 fm; a brief explanation of how these values were chosen would improve reproducibility.
- [Sec. III] The text surrounding Eq. (45) refers to 'B-orthogonality' and the footnote reference to Ref. [31], but the definition of the B matrix is not given in the present paper; a self-contained definition would be helpful.
Circularity Check
No significant circularity: the basis is benchmarked against independent DKB B-splines, analytical formulas, and Green's-function results; λ is a convergence parameter, not fitted to the reported observables.
full rationale
The derivation chain is self-contained. The new ingredient is the modified Sturmian basis (22), whose large/small-component indices lL,S are set by the finite-nucleus asymptotic relations (18)-(20); this is a basis construction, not a result derived from the target quantities. The free scale λ is selected by convergence-plateau scans (Figs. 1-2, Table VI) and is not optimized to reproduce any of the tabulated observables. The central tests are external: binding energies are compared with DKB B-splines and analytic point-nucleus energies (Table I); static polarizability with Szmytkowski's analytic formulas (Table II); hyperfine A(Zα) with Shabaev's formula and DKB (Table III); g factor with analytic Eqs. (62)-(63) and DKB (Table IV); shielding S(Zα) with Eq. (65) and Moskovkin et al. (Table V). The VP density is computed from the spectral formula (66)-(72) with the same regularization used in earlier finite-basis work; the paper cites Mohr et al.'s independent Green's-function method [33] as the comparison reference. The only self-references ([31] for Gaussian-basis VP curves and a B-orthogonality footnote) are non-load-bearing numerical benchmarks, not inputs to the derivation. The manuscript's own admission that exact finite-nucleus Sturmians would require Whittaker functions (Sec. II.C) is a stated design choice, not circularity. A verification gap remains: the VP figures do not overlay [33] and no separate n-convergence table is given for the VP density, but this is an evidentiary/completeness issue, not a case where a prediction reduces to its input by construction.
Assumptions & free parameters
free parameters (3)
- lambda (Sturmian scale parameter) =
varies per case; e.g., 2.74 for 1s energy (Table I), 94 for VP density (Fig. 3), 0.2 to 20 for HFS (Table III)
- delta (VP subtraction step) =
10^-6
- kappa cutoff =
|kappa| <= 3
assumptions (5)
- standard math Coulomb Sturmian functions are complete for the point-Coulomb radial Schrodinger problem (and their relativistic analogues for the point-nucleus Dirac problem).
- domain assumption The finite-nucleus Dirac wavefunction has zero-distance asymptotics P ~ r^{l_L+1}, Q ~ r^{l_S+1} as given by Eqs. (19)-(20), following Grant [19, Sec. 5.4].
- domain assumption The generalized Rayleigh-Ritz eigenvalue problem Hv = epsilon Cv for the Dirac Hamiltonian yields the physical spectrum without variational collapse or spurious states.
- ad hoc to paper A flat plateau or a maximum in the computed value as a function of lambda marks the converged basis-set limit.
- ad hoc to paper Equation (72), with C-symmetrization (71) and delta=10^-6, isolates the (Zalpha)^{>=3} Wichmann-Kroll vacuum polarization density.
Cite this review
Pith. "Pith review of Sturmian basis set for the Dirac equation with finite nuclear size: Application to polarizability, Zeeman and hyperfine splitting, and vacuum polarization." pith.science (2026). https://pith.science/paper/LUSHAHB3
@misc{pith2026250611988,
author = {Pith},
title = {Pith review of: Sturmian basis set for the Dirac equation with finite nuclear size: Application to polarizability, Zeeman and hyperfine splitting, and vacuum polarization},
year = {2026},
howpublished = {\url{https://pith.science/paper/LUSHAHB3}},
note = {Machine review of arXiv:2506.11988}
}
abstract
We investigate the application of the Sturmian basis set in relativistic atomic structure calculations. We propose a simple implementation of this approach and demonstrate its ability to provide various quantities for hydrogen-like ions, including binding energies, static dipole polarizability, $g$ factor, hyperfine splitting, and nuclear magnetic shielding. Finally, we calculate the all-order (Wichmann-Kroll) vacuum polarization charge density, which was a challenge for the finite-basis-set approach until recently. Comparison of the obtained results with the previously published numerical and analytical calculations is presented. All calculations are performed with the finite size of the nucleus and can in principle be extended to arbitrary binding potentials.
Figures
Figures from the paper (3 more)
Reference graph
Works this paper leans on
-
[33]
G. W. F. Drake, ed., Springer handbook of atomic, molecular and optical physics , Springer Handbooks of Atomic, Molecular, and Optical Physics (Springer Sci- ence+Business Media Inc., New York, NY, 2006)
work page 2006
-
[30]
M. F. Herbst, J. E. Avery, and A. Dreuw, Phys. Rev. A 99, 012512 (2019)
work page 2019
- [31]
-
[1]
W. R. Johnson and J. Sapirstein, Phys. Rev. Lett. 57, 1126 (1986)
1986
-
[2]
J. Sapirstein and W. R. Johnson, Journal of Physics B: Atomic, Molecular and Optical Physics 29, 5213 (1996)
work page 1996
-
[3]
V. M. Shabaev, I. I. Tupitsyn, V. A. Yerokhin, G. Plu- nien, and G. Soff, Phys. Rev. Lett. 93, 130405 (2004)
work page 2004
-
[4]
I. P. Grant, Journal of Physics B: Atomic, Molecular and Optical Physics 42, 055002 (2009)
work page 2009
-
[5]
A. V. Volotka, D. A. Glazov, G. Plunien, and V. M. Shabaev, Annalen der Physik 525, 636 (2013)
2013
Show all 52 references
-
[6]
Shabaev, in Comprehensive Computational Chemistry (First Edition) , edited by M
V. Shabaev, in Comprehensive Computational Chemistry (First Edition) , edited by M. Y´ a˜ nez and R. J. Boyd (El- sevier, Oxford, 2024) first edition ed., pp. 94–128
2024
-
[7]
V. A. Yerokhin and V. M. Shabaev, Journal of Physical and Chemical Reference Data 44, 033103 (2015)
2015
-
[8]
V. A. Yerokhin, K. Pachucki, and V. Patk´ oˇ s, Annalen der Physik 531, 1800324 (2019)
2019
-
[9]
D. A. Glazov, A. V. Volotka, V. M. Shabaev, I. I. Tupit- syn, and G. Plunien, Physical Review A 81, 062112 (2010)
2010
-
[10]
V. M. Shabaev, D. A. Glazov, G. Plunien, and A. V. Volotka, Journal of Physical and Chemical Reference Data 44, 031205 (2015)
2015
-
[11]
V. A. Yerokhin, K. Pachucki, M. Puchalski, C. H. Keitel, and Z. Harman, Physical Review A 102, 022815 (2020)
2020
-
[12]
D. A. Glazov, D. V. Zinenko, V. A. Agababaev, A. D. Moshkin, E. V. Tryapitsyna, A. M. Volchkova, and A. V. Volotka, Atoms 11 (2023), 10.3390/atoms11090119
2023 doi
-
[13]
QED effects in quadratic Zeeman splitting in highly charged hydrogen-like ions,
V. A. Agababaev, E. A. Prokhorchuk, D. A. Glazov, A. V. Malyshev, V. M. Shabaev, and A. V. Volotka, “QED effects in quadratic Zeeman splitting in highly charged hydrogen-like ions,” (2025), arXiv:2505.08929 [physics.atom-ph]
2025
-
[14]
V. A. Yerokhin, A. N. Artemyev, V. M. Shabaev, and G. Plunien, Phys. Rev. A 72, 052510 (2005)
2005
-
[15]
A. V. Volotka, D. A. Glazov, G. Plunien, V. M. Shabaev, and I. I. Tupitsyn, The European Physical Journal D - Atomic, Molecular, Optical and Plasma Physics 38, 293 (2006)
2006
-
[16]
M. G. Kozlov, V. A. Yerokhin, M. Y. Kaygorodov, and E. V. Tryapitsyna, Phys. Rev. A 110, 062805 (2024)
2024
-
[17]
Ferenc, M
D. Ferenc, M. Salman, and T. Saue, Phys. Rev. A 111, L040802 (2025)
2025
-
[18]
Quiney, I
H. Quiney, I. Grant, and S. Wilson, Physica Scripta 36, 460 (2006)
2006
-
[19]
Grant, Relativistic Quantum Theory of Atoms and Molecules: Theory and Computation , Springer Series on Atomic, Optical, and Plasma Physics (Springer New York, 2007)
I. Grant, Relativistic Quantum Theory of Atoms and Molecules: Theory and Computation , Springer Series on Atomic, Optical, and Plasma Physics (Springer New York, 2007)
2007
-
[20]
Rotenberg, Annals of Physics 19, 262 (1962)
M. Rotenberg, Annals of Physics 19, 262 (1962)
1962
-
[21]
Rotenberg, Advances in atomic and molecular physics 6, 233 (1970)
M. Rotenberg, Advances in atomic and molecular physics 6, 233 (1970)
1970
-
[22]
Truskova, USSR Computational Mathematics and Mathematical Physics 24, 49 (1984)
N. Truskova, USSR Computational Mathematics and Mathematical Physics 24, 49 (1984)
1984
-
[23]
A. I. Sherstyuk, Theoretical and Mathematical Physics 56, 813 (1983)
1983
-
[24]
Sherstyuk, Optics and Spectroscopy 87, p
A. Sherstyuk, Optics and Spectroscopy 87, p. 696–702 (1999)
1999
-
[25]
A. A. Krylovetsky, N. L. Manakov, and S. I. Marmo, Journal of Experimental and Theoretical Physics 92, 37 12 (2001)
2001
-
[26]
left” and “right
below. The Coulomb Sturmians are defined as the solution to a Sturm-Liouville problem: − d2 dr2 + l(l + 1) r2 − 2E0 + αn,lV (r) Snl(r) = 0 (1) (E0 and αn,l are parameters) with boundary conditions Snl(r)|r=0 = Snl(r)|r→∞ = 0. (2) The solution to this problem for the Coulomb po...
1994
-
[27]
I. P. Grant and H. M. Quiney, Phys. Rev. A 62, 022508 (2000)
2000
-
[28]
Szmytkowski, Journal of Physics B: Atomic, Molecular and Optical Physics 30, 825 (1997)
R. Szmytkowski, Journal of Physics B: Atomic, Molecular and Optical Physics 30, 825 (1997)
1997
-
[29]
208Pb nuclear charge ra- dius revisited: closing the fine-structure-anomaly gap,
Z. Sun, K. Beyer, Z. A. Mandrykina, I. A. Valuev, C. H. Keitel, and N. S. Oreshkina, “ 208Pb nuclear charge ra- dius revisited: closing the fine-structure-anomaly gap,” (2025), arXiv:2504.19977 [physics.atom-ph]
2025
-
[32]
V. K. Ivanov, S. S. Baturin, D. A. Glazov, and A. V. Volotka, Phys. Rev. A 110, 032815 (2024)
2024
-
[34]
P. J. Mohr, G. Plunien, and G. Soff, Physics Reports 293, 227 (1998)
1998
-
[35]
Hannesson and S
T. Hannesson and S. M. Blinder, International Journal of Quantum Chemistry 15, 7 (1979)
1979
-
[36]
G. W. F. Drake and S. P. Goldman, Phys. Rev. A 23, 2093 (1981)
1981
-
[37]
Prudnikov, I
A. Prudnikov, I. Brychkov, and O. Marichev, Inte- grals and Series, Vol. 2: Special Functions (Gordon and Breach Science Publishers, New York, 1986)
1986
-
[38]
Kemble, The Fundamental Principles of Quantum Me- chanics: With Elementary Applications , Dover Books on Science (Dover, 1958)
E. Kemble, The Fundamental Principles of Quantum Me- chanics: With Elementary Applications , Dover Books on Science (Dover, 1958)
1958
-
[39]
The mpmath development team, mpmath: a Python library for arbitrary-precision floating-point arithmetic (version 1.3.0) (2023), http://mpmath.org/
2023
-
[40]
CO- DATA Recommended Values of the Fundamental Physi- cal Constants: 2022,
P. Mohr, D. Newell, B. Taylor, and E. Tiesinga, “CO- DATA Recommended Values of the Fundamental Physi- cal Constants: 2022,” (2024), arXiv:2409.03787 [hep-ph]
2024 arXiv
-
[41]
V. B. Berestetskii, E. M. Lifshitz, and L. P. Pitaevskii, Quantum electrodynamics, Course of Theoret- ical Physics, Vol. 4 (Pergamon Press, Oxford, UK, 1982)
1982
-
[42]
Szmytkowski, Journal of Mathematical Chemistry 42, 397 (2007)
R. Szmytkowski, Journal of Mathematical Chemistry 42, 397 (2007)
2007
-
[43]
V. M. Shabaev, Journal of Physics B: Atomic, Molecular and Optical Physics 27, 5825 (1994)
1994
-
[44]
Glazov and V
D. Glazov and V. Shabaev, Physics Letters A 297, 408 (2002)
2002
-
[45]
S. G. Karshenboim, R. N. Lee, and A. I. Milstein, Phys. Rev. A 72, 042101 (2005)
2005
-
[46]
D. L. Moskovkin, N. S. Oreshkina, V. M. Shabaev, T. Beier, G. Plunien, W. Quint, and G. Soff, Phys. Rev. A 70, 032105 (2004)
2004
-
[47]
V. A. Yerokhin, K. Pachucki, Z. Harman, and C. H. Keitel, Physical Review Letters 107, 043004 (2011)
2011
-
[48]
V. A. Yerokhin, K. Pachucki, Z. Harman, and C. H. Keitel, Phys. Rev. A 85, 022512 (2012)
2012
-
[49]
Nu- clear magnetic shielding in highly charged ions,
A. M. Volchkova, D. A. Glazov, and V. M. Shabaev, “Nu- clear magnetic shielding in highly charged ions,” (2021), arXiv:2108.12311 [physics.atom-ph]
2021 arXiv
-
[50]
G. A. Rinker and L. Wilets, Phys. Rev. A 12, 748 (1975)
1975
-
[51]
E. H. Wichmann and N. M. Kroll, Phys. Rev. 101, 843 (1956)
1956
-
[52]
Soff and P
G. Soff and P. J. Mohr, Phys. Rev. A 38, 5066 (1988)
1988
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.