Pith. sign in

REVIEW 3 major objections 3 minor 3 references

Nondegenerate, Lamb-shift Solution of the Dirac Hydrogen Atom

T0 review · 3 major / 3 minor · reviewed 2026-08-16 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper claims that adding a state-dependent Lamb-shift constant to the Coulomb potential produces a hydrogen-atom solution of the Dirac equation in which every energy level is distinct.

desk verdict The algebra is consistent, but the state-dependent λ means each state solves a different Hamiltonian, so the 'nondegenerate solution' is the input Lamb shift reproduced, not derived. read the letter →

arxiv 2504.18477 v1 pith:27WLOL5W submitted 2025-04-25 quant-ph physics.atom-ph

Nathan S. Roberts This is my paper · ORCID
classification quant-phphysics.atom-ph MSC 81Q0533C4581V45 PACS 31.30.jf03.65.Pm
keywords DiracequationhydrogenatomLambshiftdegeneracyremovalLaguerrepolynomialsrelativisticquantummechanicsfinestructureenergylevels
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper sets out to remove the last degeneracy in the relativistic hydrogen atom. It adds a dimensionless Lamb-shift constant $\lambda$, multiplied by $\mu c^2$, to the Coulomb potential and solves the resulting radial Dirac equations with Laguerre polynomials, obtaining the closed energy formula $E=\mu c^2\left(1/\sqrt{1+(Z\alpha/(n-|k|+\sqrt{k^2-(Z\alpha)^2}))^2}+\lambda-1\right)$. The claim is that every state now receives its own energy, with the historically degenerate $2s_{1/2}$ and $2p_{1/2}$ pair split by the classic 1058 MHz. The radial functions reduce to the Schrödinger hydrogen eigenfunctions in the nonrelativistic limit, so the paper presents the result as one route from the Dirac equation to observed atomic spectra.

What carries the argument

The central object is the dimensionless Lamb shift $\lambda$ of Eq. (14), which packages the Bethe logarithm, form-factor terms, and Uehling and vertex corrections into one state-dependent number. The carrying mechanism is the Laguerre-polynomial ansatz and the coefficient-matching procedure that turns the trial radial functions $R_\pm=N_\pm e^{-bx}x^{u/2-1}(C_\pm y+x y')$ into the standard Laguerre equation $xy''+(u+1-x)y'+gy=0$. Requiring the unwanted low-order and high-order terms to vanish fixes $u=2\sqrt{k^2-(Z\alpha)^2}$, $b=1/2$, and $g=vw-u/2$; combining $w=1/\sqrt{1+(Z\alpha/(g+u/2))^2}$ with $w=E/(\mu c^2)-\lambda$ yields the energy formula. The dimensionless coordinate $x=2r/(a_0 v)$ is what allows the solution to reduce cleanly to the Schrödinger wavefunctions in the appropriate limits.

What would settle it

A decisive check is to plug the $2s_{1/2}$ radial function from (37) into equation (16) using the $2p_{1/2}$ value of $\lambda$ from Table I; if the equation is not satisfied identically, the two states do not come from one Hamiltonian. A reader can perform this substitution with the tables and formulas given.

Watch

Extended reading notes

Core claim

On the paper's own terms, the discovery is a solution of the radial Dirac equations for a potential $V=-Z\alpha c\hbar/r+\lambda\mu c^2$, where $\lambda$ is a dimensionless Lamb shift computed state by state from quantum electrodynamics. The radial functions take the closed form $R_\pm = N_\pm \sqrt{e^{-x}x^{u-2}}\left[\pm(v-k)L_g^u(x)-(g+u)L_{g-1}^u(x)\right]$ with $g=n-|k|$, $u=2\sqrt{k^2-(Z\alpha)^2}$, $v=\sqrt{(g+u/2)^2+(Z\alpha)^2}$, and $x=2r/(a_0 v)$. Substituting $w=1/\sqrt{1+(Z\alpha/(g+u/2))^2}$ back into $w=E/(\mu c^2)-\lambda$ yields the energy formula (40), in which $n$, $k$, and $\lambda$ together single out every state. The paper also shows that the limit $\lambda\to0$ recovers the Sommerfeld fine-structure levels and that $k\to\ell$, $\alpha\to0$ recovers the Schrödinger wavefunctions.

Load-bearing premise

The load-bearing premise is that the dimensionless Lamb shift $\lambda$, which the paper computes separately for each state, can be inserted as a single constant term in one effective potential; if $\lambda$ is not the same for all states, the various 'solutions' are solutions of different equations, not one hydrogen atom.

Editorial extensions

If this is right

  • Every hydrogen bound state $n\ell_j$ has its own predicted energy, with the Dirac degeneracy of $2s_{1/2}$ and $2p_{1/2}$ replaced by the 1058 MHz Lamb split.
  • The single formula (40) contains the Lamb shift directly, so no separate post-Dirac correction is needed to compare with measured intervals.
  • The radial eigenfunctions (37)–(38) give normalized probability densities for individual states, plotted here for the low-lying levels.
  • In the limit $k\to\ell$, $\alpha\to0$, the wavefunctions coincide with Schrödinger's hydrogen functions; in the limit $\lambda\to0$, the energies coincide with the Sommerfeld fine-structure values.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Beyond the paper: because $\lambda$ is state dependent and built from measured and calculated quantum-electrodynamic corrections, the energy formula reproduces known Lamb splittings by construction; the new content is the closed-form packaging, not a fresh numerical prediction.
  • Beyond the paper: applying the same ansatz to hydrogenic ions or muonic hydrogen, where $Z\alpha$ is larger, would make the radial shifts visible and would test whether a state-by-state $\lambda$ can be assigned a single Hamiltonian.
  • Beyond the paper: transition rates computed from (37)–(38), such as $2s_{1/2}\to1s_{1/2}$, should be compared with standard quantum-electrodynamic values; any mismatch would locate the cost of moving the Lamb shift into the potential.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 3 minor

Summary. The paper claims a new solution of the Dirac equation for the hydrogen atom that removes all degeneracy by adding a dimensionless Lamb-shift term λ to the Coulomb potential. The derivation defines λ from known QED corrections in Eq. (14), inserts it as a constant shift in the potential in Eqs. (15)–(16), solves the radial Dirac equation with Laguerre polynomials, and obtains the eigenfunctions (37)–(38) and the energy formula (40). The authors show that in the limits λ → 0 and α → 0 the wavefunctions reduce to the Schrödinger results. The central claim is that Eq. (40) provides a single nondegenerate formula for all hydrogen energy levels.

Significance. If the proposal were correct, it would offer a compact formula for hydrogenic energy levels including the Lamb shift, with explicit wavefunctions. The algebraic solution from Eq. (15) to Eq. (40) appears internally consistent for a fixed λ, and the reduction to Schrödinger wavefunctions is a useful check. However, the central physical claim is not supported because λ is state-dependent and the degeneracy removal is imported from input data rather than derived from the Dirac equation. The paper also does not produce a falsifiable prediction: the 1058 MHz splitting in Table IV is exactly the input Lamb shift. The manuscript therefore does not meet the standard for a substantive advance in bound-state QED.

major comments (3)
  1. [§2, Eqs. (14)–(16)] The potential V = −Zα cℏ/r + λ μc² is inserted into the Dirac equation, but λ in Eq. (14) and Table I is state-dependent (it depends on n, ℓ, j). A single Hamiltonian requires one common λ; using a different λ for each state means that each state solves a different operator, so there is no single Dirac hydrogen atom being solved. The paper does not identify any fixed potential that produces all the levels in Table IV.
  2. [§3, Eq. (40) and Table IV] Equation (40) is the ordinary Dirac energy plus λ μc². Because λ is the input Lamb shift for each state, the 2s1/2–2p1/2 splitting in Table IV is exactly the 1058 MHz input, not a prediction of the Dirac equation. The abstract's statement that the paper 'use[s] the Lamb shift to give each atomic state a unique energy level' confirms that the degeneracy removal is imported from experiment rather than derived.
  3. [§4, final paragraph] The limit λ → 0 returning to the degenerate Dirac levels shows that the proposed solution contains no mechanism that removes degeneracy; the nondegeneracy is entirely due to the added shift. Since the wavefunctions in Eqs. (37)–(38) are independent of λ, the claimed 'new solution' is the standard Dirac solution with a state-dependent energy offset, not a structurally new solution.
minor comments (3)
  1. [Table I] The last column has formatting issues; for example the entry '66 180×10−15' is ambiguous and should be written as a standard mantissa and exponent.
  2. [Figure 1 caption] The caption admits that the claimed leftward shifts are 'too small to see clearly at this scale,' which weakens the visual evidence for the nondegenerate wavefunctions.
  3. [Notation] The notation '1s12', '2p12' etc. should be typeset as 1s_{1/2}, 2p_{1/2} to avoid confusion with multiplication or powers.

Circularity Check

3 steps flagged · score 8.0 of 10

Degeneracy removal is the input Lamb shift added to the ordinary Dirac energy, not a derived prediction.

  1. self definitional [Abstract]
    "We work in terms of dimensionless quantities and use the Lamb shift to give each atomic state in our model a unique, nondegenerate energy level."

    The paper's central claim of removing all degeneracy is stated here as being achieved by using the Lamb shift, not by a new consequence of the Dirac equation. The final energy formula (40) is the ordinary Dirac energy plus lambda mu c^2, so the uniqueness of each level is imposed by inserting the measured/known Lamb shift as an input, not derived from solving the equation.

  2. fitted input called prediction [Eqs. (14)-(16), Table IV, Section 4 (Discussion)]
    "Multiplying 𝜆 by 𝜇𝑐2, we combine it with Coulomb's law: 𝑉=−𝑍𝛼𝑐ℏ𝑟+𝜆𝜇𝑐2 (15) The effective potential 𝑉 in (15) is identical to (13); it is simply written more concisely using 𝜆 (14)."

    Lambda in (14) is assembled from the known Bethe logarithm, Uehling, vertex, and other Lamb-shift contributions, and is state-dependent. It is then inserted as a constant in the potential, and the final energy (40) differs between states only by the additive lambda mu c^2 term. Table IV's 2s1/2-2p1/2 difference is exactly [lambda_2s - lambda_2p] mu c^2, which the paper identifies as the classic 1058 MHz Lamb shift. Thus the 'prediction' is the input value returned as output.

1 more flagged steps
  1. other [Eqs. (34), (37)-(38), Section 3 (Results)]
    "We observe that 𝑤, which contains the Lamb shift 𝜆 in (17), cancels out in (34). Values like 𝑣 (36) can be expressed independently of 𝜆, so 𝜆 vanishes from our eigenfunctions."

    The radial eigenfunctions, which are the actual solution of the Dirac equation, do not depend on lambda at all. The only place lambda appears is the additive energy term in (40). This demonstrates that the removal of degeneracy is not a property of the new wavefunctions or of a new Dirac solution; it is entirely supplied by the externally inserted Lamb shift constant.

full rationale

The algebraic reduction of the Dirac equation to Laguerre-polynomial radial functions is internally consistent, and the eigenfunctions do reduce to Schrodinger wavefunctions in the stated limits. However, the paper's central claim, that it 'finally removes all degeneracy from the hydrogen atom,' is circular by construction. Lambda in (14) is explicitly defined from known Lamb-shift contributions and is state-dependent, so inserting it into the potential (15) means each state is governed by a different Hamiltonian. The radial functions (37)-(38) are independent of lambda, and the final energy (40) is simply the standard Dirac energy plus lambda mu c^2. Consequently, the 1058 MHz 2s1/2-2p1/2 splitting in Table IV is the classic Lamb shift entered at (14), not a prediction of the Dirac equation; the paper itself says the levels 'replicate the classic Lamb shift' and that lambda -> 0 restores the degenerate Sommerfeld levels. The nondegeneracy is therefore imported by definition rather than derived. Score 8 reflects that the central result reduces by construction, even though the Laguerre algebra and Schrodinger-limit checks are not themselves circular.

Assumptions & free parameters 2 free parameters · 3 assumptions · 0 invented entities

The central claim rests on a state-dependent constant lambda that is inserted into the potential and then recovered in the energy. No new entities are introduced, but the derivation relies on an ad hoc assumption that the Lamb shift can be represented as a constant potential shift, which is not self-evident and is not derived. The free parameters are the external input values for lambda and the Bethe logarithms.

free parameters (2)
  • lambda (dimensionless Lamb shift) per state = 2s1/2: 8.461e-12; 2p1/2: -1.04e-13; 1s1/2: 6.618e-11; 2p3/2: 1.01e-13 (from Table I)
    lambda is defined in Eq. (14) using known QED quantities, and the values used in Table I are taken from external calculations. It is inserted into the potential and then appears linearly in the final energy formula (40), so it is an input parameter chosen to reproduce the known Lamb shift.
  • Bethe logarithm values beta_nℓ = beta_1s = 2.984128555, beta_2s = 2.811769893, beta_2p = -0.030016708 (Eq. 9)
    These numeric values are taken from the cited literature and used in lambda. They are not derived in this paper and are treated as external inputs.
assumptions (3)
  • standard math The Dirac equation with a static Coulomb potential describes the hydrogen atom.
    This is the starting point of the paper, Eq. (6), and is a standard result in relativistic quantum mechanics.
  • ad hoc to paper The known Lamb shift can be represented as a constant added to the Coulomb potential.
    The paper inserts lambda μc² into the potential in Eq. (15). This is not a consequence of the Dirac equation; it is an ad hoc modeling choice, and the value of lambda is state-dependent, so the potential is different for each state.
  • domain assumption The numerical values of QED corrections (Bethe logarithm, form factors, etc.) in Eq. (14) are correct and applicable.
    The paper adopts the standard QED results for the Lamb shift from the cited literature, without re-deriving them. This is a domain assumption: the paper trusts the external calculations.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Nondegenerate, Lamb-shift Solution of the Dirac Hydrogen Atom." pith.science (2026). https://pith.science/paper/27WLOL5W

@misc{pith2026250418477,
  author       = {Pith},
  title        = {Pith review of: Nondegenerate, Lamb-shift Solution of the Dirac Hydrogen Atom},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/27WLOL5W}},
  note         = {Machine review of arXiv:2504.18477}
}
read the original abstract

When the Dirac equation was first published in 1928, three solutions appeared immediately within the same year, each describing the most important problem in physics at that time: the hydrogen atom. These solutions lifted some of the degeneracy from earlier atomic models, but not all of it--they still predicted the same degenerate energy levels for the 2s1/2 and 2p1/2 states, for example. In this paper, we introduce a new solution of the Dirac equation, which finally removes all degeneracy from the hydrogen atom. We work in terms of dimensionless quantities and use the Lamb shift to give each atomic state in our model a unique, nondegenerate energy level. We obtain radial eigenfunctions in terms of the Laguerre polynomials, demonstrate how they can be reduced to the Schrodinger wavefunctions by applying limits, and plot our results.

Figures

Figures reproduced from arXiv: 2504.18477 by the authors.

Figure 1
Figure 1. FIG 1. The Schrodinger ground state 1s (solid red curve) [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

3 extracted references · 3 canonical work pages

  1. [1]

    Paul A. M. Dirac, The quantum theory of the electron, Proc. R. Soc. A 117, 610 (1928). [2] Walter Gordon, The energy levels of the hydrogen atom according to Dirac’s quantum theory of the electron (in German), Z. Phys. 48, 11 (received 23 Feb 1928). [3] Charles G. Darwin, The wave equation of the electron, Proc. R. Soc. A 118, 654 (received 6 Mar 1928). [...

    work page 1928

  2. [2]

    apparent

    METHODS Equation (13) depicts the Coulomb field centered at X exchanging a single photon with an electron to change its momentum (green). Sometimes however, a second photon escapes during transit to emerge as bremsstrahlung (blue). Other times, the electron reabsorbs this second photon before it can escape (vertex function). The photon can even split temp...

    work page 1930

  3. [11]

    Feynman, Space-time approach to quantum electrodynamics, Phys

    Richard P. Feynman, Space-time approach to quantum electrodynamics, Phys. Rev. 76, 769 (1949). [12] Walter Gordon, The current of Dirac’s electron theory (in German), Z. Phys. 50, 630 (1928). [13] Julian Schwinger, On quantum-electrodynamics and the magnetic moment of the electron, Phys. Rev. 73, 416 (1947). [14] James D. Bjorken and Sidney D. Drell, 𝑅𝑒𝑙𝑎...

    work page 1949

Pith tools

Reviewed August 16, 2026 · model on record in the stance chip above.