pith. sign in

arxiv: 2604.11063 · v1 · submitted 2026-04-13 · 🧮 math.NA · cs.NA

An Adaptive Log-Laguerre Spectral Method for the Radial Dirac Equation: Resolving Asymptotic Decay and Core Singularities in Atomic Calculations

Sheng Chen , Sihong Shao , Shuai Wu This is my paper

Pith reviewed 2026-05-10 16:01 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords radial Dirac equationspectral methodLaguerre functionslog-orthogonal functionsatomic calculationsnumerical analysishigh-precisionrelativistic quantum chemistry
0
0 comments X

The pith

An adaptive spectral method using optimized Laguerre scaling and log-orthogonal functions solves the radial Dirac equation to 10^{-10} relative accuracy across atomic potentials.

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

The paper shows that standard basis approaches cannot simultaneously handle non-polynomial singularities at the origin and state-dependent multi-scale decays at infinity when solving the radial Dirac equation. It introduces an adaptive generalized Laguerre method that dynamically tunes scaling factors for the semi-infinite domain together with log-orthogonal functions that approximate r^s singularities without prior knowledge of the exponent s. An inverse operator formulation removes spurious states. This combination restores exponential convergence and reaches relative accuracies of 10^{-10} in Hartree units for Coulomb, finite-nucleus, and screened potentials, supplying a clean numerical kernel for relativistic atomic calculations.

Core claim

The adaptive generalized Laguerre spectral method with dynamically optimized scaling factors, combined with Log-Orthogonal Functions and an inverse operator formulation, resolves both the asymptotic decay on [0, infinity) and the core r^s singularities, restores exponential convergence, and consistently achieves relative accuracies of 10^{-10} without domain truncation or basis pollution for Coulomb, finite-nucleus, and screened potentials.

What carries the argument

Adaptive generalized Laguerre spectral method with dynamically optimized scaling factors, augmented by Log-Orthogonal Functions that intrinsically approximate non-polynomial r^s singularities.

If this is right

  • Supplies a non-polluting computational kernel for high-precision atomic structure calculations.
  • Enables reliable generation of pseudopotentials and all-electron data for relativistic quantum chemistry.
  • Extends without modification to Coulomb, finite-nucleus, and screened potentials.
  • Eliminates the need for arbitrary domain truncation while preserving spectral accuracy.

Where Pith is reading between the lines

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

  • The same singularity-handling technique could be tested on time-dependent or multi-electron Dirac problems.
  • Replacing traditional bases with this adaptive pair in existing atomic codes might raise benchmark precision.
  • The inverse operator step may generalize to other singular eigenvalue problems in quantum mechanics.

Load-bearing premise

That adaptive optimization of Laguerre scaling factors together with log-orthogonal functions can capture arbitrary state-dependent decays and unknown power singularities for any potential without introducing pollution or requiring advance knowledge of the singularity exponent.

What would settle it

A test case on a screened potential with an extreme multi-scale decay or an unapproximated singularity exponent where the computed eigenvalues deviate from reference values by more than 10^{-10} or spurious states appear.

Figures

Figures reproduced from arXiv: 2604.11063 by Sheng Chen, Shuai Wu, Sihong Shao.

Figure 1
Figure 1. Figure 1: Comparison of relativistic large-component radial wavefunctions [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Comparison of relativistic large-component radial wavefunctions [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Schematic of the bipartite domain decomposition on the semi-unbounded interval. The physical [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Example: Coulomb Potential for Uranium ion ( [PITH_FULL_IMAGE:figures/full_fig_p011_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Example: Coulomb Potential. Maximum relative error of ALLSM for different [PITH_FULL_IMAGE:figures/full_fig_p012_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Example: Gaussian Nuclear Model. Left: Comparison with established benchmarks from Andrae et al. [23]; the slight error increase at high basis sizes occurs as the ALLSM precision exceeds the literature’s limits. Right: Internal consistency test against an ultra-high-fidelity reference (N1 = 120, N2 = 60), confirming strict convergence. that of the Coulomb potential solved via ALSM. This strong similarity i… view at source ↗
Figure 7
Figure 7. Figure 7: Convergence behavior of the numerical error for different nuclear potentials using ALSM. [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Example: Yukawa Potential. Left: Maximum relative error of the ground-state energy for Case 1 (V0 = 0.1c, λ = 0.01c) and Case 2 (V0 = 0.4c, λ = 0.07c). The results demonstrate that ALLSM maintains high precision across different screening regimes. Right: Evolution of the adaptive scaling parameter as a function of the degrees of freedom N2. The scaling factor effectively converges from a sub-optimal starti… view at source ↗
Figure 9
Figure 9. Figure 9: Example: Yukawa Potential. Left: Comparative evolution of the ground-state energy Egs and the optimal scaling parameters, β and β 2 , as a function of the screening parameter λ. The synchronized decaying trend suggests that the adaptive scaling factor effectively tracks the energy shift as the system approaches the continuum. Right: Radial wavefunctions for various λ values approaching λcrit. As the bindin… view at source ↗
Figure 10
Figure 10. Figure 10: The number of converged physical eigenvalues [PITH_FULL_IMAGE:figures/full_fig_p017_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: The number of eigenvalues with degree of freedom for fixed ratio. [PITH_FULL_IMAGE:figures/full_fig_p017_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Topological stabilization mechanism of the Dirac spectrum via operator reformulations. [PITH_FULL_IMAGE:figures/full_fig_p019_12.png] view at source ↗
read the original abstract

The high-precision solution of the radial Dirac equation is fundamental to relativistic quantum chemistry, essential for reliable pseudopotential generation and all-electron electronic structure methods. However, standard basis-set approaches struggle to simultaneously capture two distinct physical regimes: the non-polynomial singularities at the origin and the state-dependent, multi-scale asymptotic decay of wavefunctions on semi-infinite domains. In this work, we propose a high-precision adaptive spectral-element framework designed to rigorously resolve these spatial challenges. To capture the diverse exponential decay behavior on $[0, \infty)$ without arbitrary domain truncation, an adaptive generalized Laguerre spectral method is introduced, dynamically optimizing the basis scaling factors. Concurrently, near-origin non-polynomial {$r^s$} singularities are resolved utilizing Log-Orthogonal Functions, a basis that intrinsically approximates complex singular behaviors without requiring prior knowledge of the exact analytical exponent {$s$}. Furthermore, the framework incorporates an inverse operator formulation to guarantee spectral purity and eliminate spurious states. Validated across diverse physical regimes, including Coulomb, finite-nucleus, and screened potentials, the proposed method restores exponential convergence and consistently achieves relative accuracies of $10^{-10}$ {in Hartree atomic units or electron volts}. This work provides a robust, non-pollution computational kernel for atomic structure calculations, establishing a numerical standard for generating high-precision atomic data in complex molecular simulations.

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.

Referee Report

2 major / 0 minor

Summary. The manuscript proposes an adaptive log-Laguerre spectral method for the radial Dirac equation. It introduces dynamically optimized generalized Laguerre basis functions to capture state-dependent multi-scale exponential decay on [0, ∞) without domain truncation, combined with Log-Orthogonal Functions to resolve non-polynomial r^s singularities at the origin without requiring prior knowledge of the exponent s. An inverse-operator formulation is employed to guarantee spectral purity and remove spurious modes. The method is validated on Coulomb, finite-nucleus, and screened potentials, with claims of restored exponential convergence and consistent relative accuracies of 10^{-10} in Hartree atomic units or electron volts.

Significance. If the convergence and accuracy claims hold under rigorous verification, the work would represent a useful contribution to high-precision numerical methods in relativistic quantum chemistry. The adaptive handling of both origin singularities and asymptotic decay in a single framework, without explicit singularity exponent input or artificial truncation, could streamline all-electron atomic calculations and pseudopotential generation. The emphasis on spectral purity via the inverse formulation addresses a known practical issue in Dirac spectral methods.

major comments (2)
  1. Abstract: The headline claims of exponential convergence and 10^{-10} relative accuracy across multiple potentials rest entirely on assertions without any accompanying error analysis, convergence plots, implementation details, or tabulated numerical results. This is load-bearing because the central contribution is the restoration of spectral accuracy for arbitrary potentials; absent supporting evidence, the performance assertions cannot be evaluated.
  2. Abstract: No explicit definition, weight function, truncation rule, or construction of the Log-Orthogonal Functions is supplied, nor is the algorithm for dynamically optimizing the generalized Laguerre scaling factors described. Without these, it is impossible to assess whether the two adaptations can jointly resolve unknown r^s singularities and state-dependent decay without introducing pollution or requiring implicit tuning, as asserted.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful review and constructive feedback. The comments focus on the need for clearer linkage between the abstract claims and the supporting material in the manuscript. We address each point below and outline the revisions we will implement to improve clarity and accessibility of the evidence.

read point-by-point responses

  1. Referee: Abstract: The headline claims of exponential convergence and 10^{-10} relative accuracy across multiple potentials rest entirely on assertions without any accompanying error analysis, convergence plots, implementation details, or tabulated numerical results. This is load-bearing because the central contribution is the restoration of spectral accuracy for arbitrary potentials; absent supporting evidence, the performance assertions cannot be evaluated.

    Authors: The manuscript provides the requested supporting evidence in the body of the paper. Section 4 contains convergence studies with plots (Figures 4--7) that demonstrate restored exponential convergence rates for the Coulomb, finite-nucleus, and screened potentials. Section 5 presents tabulated relative errors reaching 10^{-10} (or better) in Hartree units, together with a detailed error analysis. Implementation details, including basis construction and solver parameters, appear in Section 3. We will revise the abstract to add a short clause referencing these sections and the nature of the numerical verification, thereby directing readers to the evidence while preserving the abstract's brevity. revision: yes

  2. Referee: Abstract: No explicit definition, weight function, truncation rule, or construction of the Log-Orthogonal Functions is supplied, nor is the algorithm for dynamically optimizing the generalized Laguerre scaling factors described. Without these, it is impossible to assess whether the two adaptations can jointly resolve unknown r^s singularities and state-dependent decay without introducing pollution or requiring implicit tuning, as asserted.

    Authors: The abstract is intentionally concise and therefore omits the technical specifications. These are supplied explicitly in the main text: Section 2.1 defines the Log-Orthogonal Functions, states the weight function, truncation rule, and construction procedure that approximates r^s singularities without prior knowledge of s; Section 2.3 and Algorithm 1 describe the dynamic optimization of the generalized Laguerre scaling parameters, including the safeguards against mode pollution and the absence of implicit tuning. We will revise the abstract to include a brief pointer to Section 2 for these constructions, enabling immediate location of the full specifications. revision: yes

Circularity Check

0 steps flagged

No circularity: method is an explicit numerical construction with external validation

full rationale

The paper proposes an adaptive generalized Laguerre spectral method with dynamic scaling optimization and Log-Orthogonal Functions for the radial Dirac equation, plus an inverse-operator formulation. These are presented as algorithmic choices (scaling-factor optimization, basis construction without prior s knowledge, spurious-mode elimination) whose performance is then validated numerically on Coulomb, finite-nucleus, and screened potentials to 10^{-10} accuracy. No equation or central claim reduces by construction to a fitted parameter, self-definition, or self-citation chain; the derivation chain consists of standard spectral-method steps whose correctness is checked against independent physical benchmarks rather than tautologically assumed.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 1 invented entities

The approach rests on standard spectral convergence theory for Laguerre polynomials and the assumption that log-orthogonal functions can represent the required singular behavior; no free parameters are explicitly fitted in the abstract description, and no new physical entities are postulated.

axioms (2)
  • standard math Spectral methods based on generalized Laguerre polynomials converge exponentially for functions with appropriate decay at infinity.
    Invoked implicitly when claiming restoration of exponential convergence on [0, infinity).
  • domain assumption Log-orthogonal functions can approximate non-polynomial r^s singularities at the origin without prior knowledge of the exponent s.
    Central to the near-origin treatment described in the abstract.
invented entities (1)
  • Log-Orthogonal Functions no independent evidence
    purpose: To intrinsically approximate complex singular behaviors at the origin without requiring the exact analytical exponent s.
    Introduced as the basis for resolving core singularities.

pith-pipeline@v0.9.0 · 5552 in / 1427 out tokens · 37679 ms · 2026-05-10T16:01:08.777717+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

3 extracted references · 3 canonical work pages

  1. [1]

    Purgatorio—a new implementation of the Inferno algorithm.Journal of Quantitative Spectroscopy and Radiative Transfer2006,99, 658–679

    (1) Wilson, B.; Sonnad, V.; Sterne, P.; Isaacs, W. Purgatorio—a new implementation of the Inferno algorithm.Journal of Quantitative Spectroscopy and Radiative Transfer2006,99, 658–679. (2) Hohenberg, P.; Kohn, W. Inhomogeneous electron gas.Physical review1964,136, B864. (3) Jönsson, P.; He, X.; Fischer, C. F.; Grant, I. The grasp2K relativistic atomic str...

    work page 1982

  2. [2]

    Stabilized finite element method for the radial Dirac equation.Journal of Computational Physics2013,236, 426–442

    (16) Almanasreh, H.; Salomonson, S.; Svanstedt, N. Stabilized finite element method for the radial Dirac equation.Journal of Computational Physics2013,236, 426–442. (17) Ainsworth, M.; Senior, B. Aspects of an adaptive hp-finite element method: Adaptive strategy, con- forming approximation and efficient solvers.Computer Methods in Applied Mechanics and En...

    work page 1997

  3. [3]

    Finite nuclear charge density distributions in electronic structure calculations for atoms and molecules.Physics Reports2000,336, 413–525

    (23) Andrae, D. Finite nuclear charge density distributions in electronic structure calculations for atoms and molecules.Physics Reports2000,336, 413–525. (24) Roy, A. K. Accurate ro-vibrational spectroscopy of diatomic molecules in a Morse oscillator potential. Results in physics2013,3, 103–108. (25) Chen, S.; Shen, J.; Zhang, Z.; Zhou, Z. A spectrally a...

    work page 2002