Pith. sign in

REVIEW 2 major objections 2 minor 44 references

The Dirac equation admits exact propagating normalizable positive-energy wave-packet solutions in the axisymmetric potential V = -v0/ρ.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

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

T0 review · grok-4.3

2026-06-27 00:02 UTC pith:NIJSVMUF

load-bearing objection The paper claims exact Dirac wave-packet solutions in the V=-v0/ρ potential via a Helmholtz mapping, but the key cancellation step remains unverified from the abstract alone. the 2 major comments →

arxiv 2606.18470 v1 pith:NIJSVMUF submitted 2026-06-16 quant-ph hep-thmath-phmath.MP

Exact propagating Dirac wave packets in an attractive Coulomb-like potential

Siddhant Das This is my paper
classification quant-ph hep-thmath-phmath.MP
keywords Dirac equationwave packetsexact solutionsaxisymmetric potentialCoulomb-like potentialrelativistic quantum mechanicsspin decoupling
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper constructs exact positive-energy normalizable wave-packet solutions of the Dirac equation for the attractive axisymmetric potential V = -v0/ρ. These are presented as the first such solutions known for the Dirac equation in any external potential. One family uses only elementary functions, with longitudinal profiles that recover the free Schrödinger Hermite-Gauss packets in the nonrelativistic limit. The solutions display probability density independent of spin orientation and complete freezing of time evolution when the coupling reaches v0 = ħc/2. Additional solutions arise from a direct mapping of 2D Helmholtz equation solutions.

Core claim

We construct exact, positive-energy, normalizable wave-packet solutions of the Dirac equation in the axisymmetric potential V=-v0/ρ -- to our knowledge, the first such solutions in any external potential. Remarkably, one family comprises only elementary functions whose longitudinal profiles reproduce the free-Schrödinger Hermite--Gauss wave packets in the nonrelativistic limit. All packets share two striking features: (i) a probability density that is pointwise decoupled from spin orientation -- despite the inherent spin-orbit coupling of the Dirac equation -- and (ii) a complete freezing of their time evolution at the critical coupling v0→ħc/2. We also present a simple scheme that maps solu

What carries the argument

An ansatz that maps solutions of the 2D Helmholtz equation onto exact solutions of the Dirac equation for the potential V=-v0/ρ.

Load-bearing premise

The axisymmetric potential V=-v0/ρ permits an ansatz derived from the 2D Helmholtz equation to satisfy the Dirac equation exactly.

What would settle it

Substitution of the constructed wave-packet forms into the time-dependent Dirac equation for V=-v0/ρ yields a nonzero residual.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Probability density remains pointwise independent of spin orientation for all solutions.
  • Time evolution of every packet freezes completely when v0 approaches ħc/2.
  • One explicit family reduces exactly to free Schrödinger Hermite-Gauss packets in the nonrelativistic limit.
  • The Helmholtz mapping generates additional families of exact Dirac packets beyond the elementary ones.

Where Pith is reading between the lines

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

  • These exact solutions could serve as benchmarks for numerical relativistic quantum dynamics in 2D or cylindrical geometries.
  • The spin-decoupled density may simplify modeling of transport or interference in systems with strong spin-orbit effects.
  • The critical-coupling freezing suggests a regime where relativistic packets behave as stationary states despite nonzero energy.
  • The mapping technique might be tested on other axisymmetric potentials that reduce to Helmholtz in some limit.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 2 minor

Summary. The manuscript constructs exact positive-energy normalizable wave-packet solutions of the Dirac equation for the axisymmetric potential V = -v0/ρ. One family uses elementary functions whose longitudinal profiles recover free-Schrödinger Hermite-Gauss packets nonrelativistically; a second family is obtained by mapping solutions of the 2D Helmholtz equation to Dirac spinors. All solutions are reported to exhibit spin-independent probability density and complete freezing of time evolution at the critical value v0 = ħc/2.

Significance. If the central derivations are correct, the results would constitute a notable advance by supplying the first exact propagating normalizable solutions of the Dirac equation in a nontrivial external potential, together with an explicit mapping procedure that generates further solutions from the Helmholtz equation. The reported decoupling of probability density from spin and the exact freezing at critical coupling are potentially useful for understanding relativistic spin-orbit effects.

major comments (2)
  1. [Mapping scheme section] The mapping from 2D Helmholtz solutions to Dirac spinors (described after the elementary family) is load-bearing for the central claim of exact solutions. The manuscript must supply the explicit substitution of the four-component ansatz into the Dirac operator to demonstrate that every term generated by the potential (the βV contribution) and by the cylindrical derivatives cancels exactly, leaving only the Helmholtz operator acting on the profile; without this step-by-step cancellation shown, it remains possible that residual 1/ρ or derivative terms survive.
  2. [Results for the elementary family] The normalizability and positive-energy assertions for both families require explicit verification that the constructed packets are square-integrable and lie in the positive continuum; the manuscript should state the integration measure in cylindrical coordinates and confirm that the longitudinal and transverse factors separately yield finite norms for the reported parameter ranges.
minor comments (2)
  1. [Abstract] Notation for the cylindrical radius (ρ) and the potential strength (v0) should be introduced once in the introduction and used consistently; the abstract uses both without prior definition.
  2. [Nonrelativistic limit paragraph] The nonrelativistic limit statement would be strengthened by an explicit expansion of the Dirac spinor components to O(1/c) rather than a qualitative remark.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive evaluation of the work's significance and for the detailed comments. Both major points identify places where additional explicit derivations will strengthen the manuscript; we will incorporate them in the revision.

read point-by-point responses

  1. Referee: [Mapping scheme section] The mapping from 2D Helmholtz solutions to Dirac spinors (described after the elementary family) is load-bearing for the central claim of exact solutions. The manuscript must supply the explicit substitution of the four-component ansatz into the Dirac operator to demonstrate that every term generated by the potential (the βV contribution) and by the cylindrical derivatives cancels exactly, leaving only the Helmholtz operator acting on the profile; without this step-by-step cancellation shown, it remains possible that residual 1/ρ or derivative terms survive.

    Authors: We agree that the explicit substitution must be shown. In the revised manuscript we will add a dedicated paragraph (or short subsection) that inserts the four-component ansatz into the Dirac operator, carries out the term-by-term cancellation of the βV and cylindrical-derivative contributions, and isolates the 2D Helmholtz operator acting on the profile function. This will remove any ambiguity about residual terms. revision: yes

  2. Referee: [Results for the elementary family] The normalizability and positive-energy assertions for both families require explicit verification that the constructed packets are square-integrable and lie in the positive continuum; the manuscript should state the integration measure in cylindrical coordinates and confirm that the longitudinal and transverse factors separately yield finite norms for the reported parameter ranges.

    Authors: We will add an explicit verification section (or appendix) that states the cylindrical integration measure ∭ |ψ|² ρ dρ dφ dz and demonstrates that the product of the longitudinal (Hermite–Gauss) and transverse factors produces finite norms for the stated parameter ranges. We will also indicate how the positive-energy character follows from the construction (the solutions are built from the positive-energy branch of the Dirac operator). revision: yes

Circularity Check

0 steps flagged

No significant circularity; explicit mapping from Helmholtz to Dirac is presented as a derivation step

full rationale

The paper's central construction is an explicit ansatz/mapping that takes solutions of the 2D Helmholtz equation and produces candidate spinors for the Dirac equation in V=-v0/ρ. This is not self-definitional because the mapping is offered as a concrete scheme whose validity rests on algebraic cancellation of all potential and derivative terms in the Dirac operator, which is an independent verification task rather than a redefinition. No parameters are fitted to data and then relabeled as predictions, no uniqueness theorem is imported from self-citation, and the nonrelativistic limit and probability-density decoupling are presented as emergent properties rather than inputs. The derivation chain therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Review based on abstract only; no explicit free parameters, new entities, or ad-hoc axioms beyond the standard Dirac equation are mentioned. The critical value ħc/2 is a standard constant combination.

axioms (2)
  • ad hoc to paper The Dirac equation in cylindrical coordinates for axisymmetric potentials admits exact normalizable positive-energy solutions of the claimed form.
    This is the load-bearing premise enabling the construction of the wave packets.
  • ad hoc to paper Solutions of the 2D Helmholtz equation can be mapped to solutions of the Dirac equation in this potential.
    Invoked in the final sentence as the source of further exact packets.

pith-pipeline@v0.9.1-grok · 5661 in / 1391 out tokens · 48932 ms · 2026-06-27T00:02:04.399421+00:00 · methodology

0 comments
read the original abstract

We construct exact, positive-energy, normalizable wave-packet solutions of the Dirac equation in the axisymmetric potential $V=-\,v_0/\rho$ -- to our knowledge, the first such solutions in any external potential. Remarkably, one family comprises only elementary functions whose longitudinal profiles reproduce the free-Schr\"odinger Hermite--Gauss wave packets in the nonrelativistic limit. All packets share two striking features: (i) a probability density that is pointwise decoupled from spin orientation -- despite the inherent spin-orbit coupling of the Dirac equation -- and (ii) a complete freezing of their time evolution at the critical coupling $v_0\to\hbar c/2$. We also present a simple scheme that maps solutions of the 2D Helmholtz equation to further exact Dirac wave packets.

Figures

Figures reproduced from arXiv: 2606.18470 by Siddhant Das.

Figure 1
Figure 1. Figure 1: FIG. 1. Time evolution of the marginal density ( [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Marginal density ( [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

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

44 extracted references · 5 canonical work pages · 2 internal anchors

  1. [1]

    Das, (2021), arXiv:2112.13180 [quant-ph]

    S. Das, (2021), arXiv:2112.13180 [quant-ph]

    work page arXiv 2021

  2. [2]

    V. G. Bagrov and D. Gitman,The Dirac Equation and its So- lutions, Studies in Mathematical Physics, Vol. 4 (de Gruyter, Berlin, 2014)

    2014

  3. [3]

    Bia lynicki-Birula and Z

    I. Bia lynicki-Birula and Z. Bia lynicka-Birula, J. Phys. A: Math. Theor.56, 285302 (2023)

    2023

  4. [4]

    A. G. Campos, K. Z. Hatsagortsyan, and C. H. Keitel, Phys. Rev. Res.6, 023040 (2024)

    2024

  5. [5]

    G. V. Shishkin and V. M. Villalba, J. Math. Phys.30, 2373 (1989)

    1989

  6. [6]

    Greiner,Relativistic Quantum Mechanics–Wave Equations, 3rd ed., Texts and Monographs in Physics (Springer-Verlag, Berlin Heidelberg, 2000)

    W. Greiner,Relativistic Quantum Mechanics–Wave Equations, 3rd ed., Texts and Monographs in Physics (Springer-Verlag, Berlin Heidelberg, 2000)

    2000

  7. [7]

    K. B. Korotchenko and Y. P. Kunashenko, Nuovo Cimento C 34, 111 (2011)

    2011

  8. [8]

    H. A. Olsen and Y. Kunashenko, Phys. Rev. A56, 527 (1997)

    1997

  9. [9]

    Das and D

    S. Das and D. D¨ urr, Sci. Rep.9, 2242 (2019). 5

    2019

  10. [10]

    S. Das, M. N ¨oth, and D. D¨ urr, Phys. Rev. A99, 052124 (2019)

    2019

  11. [11]

    J. M. Wilkes, Eur. J. Phys.41, 035402 (2020)

    2020

  12. [12]

    Das and S

    S. Goldstein, R. Tumulka, and N. Zangh `ı, Sci. Rep.14, 3775 (2024); Found. Phys.54, 1 (2024); S. Das and S. Aristarhov, (2023), arXiv:2312.01802 [quant-ph]

    work page arXiv 2024

  13. [13]

    Drezet, Symmetry16(2024), 10.3390/sym16101325

    A. Drezet, Symmetry16(2024), 10.3390/sym16101325

    work page doi:10.3390/sym16101325 2024

  14. [14]

    Poirier and R

    B. Poirier and R. Lombardini, Entropy26, 336 (2024)

    2024

  15. [15]

    Arrival-time distributions as a probe of the preferred foliation in relativistic Bohmian mechanics

    A. Amblard and A. Drezet, (2026), arXiv:2604.17507 [quant- ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  16. [16]

    Detection Time Distribution Predicted Using Absorbing Boundary Conditions and Imaginary Potentials

    A. Jozani and R. Tumulka, (2026), arXiv:2603.22044 [quant- ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  17. [17]

    The empirical consequences of the no-signaling hypothesis,

    W. Cavendish, “The empirical consequences of the no-signaling hypothesis,” (2026), forthcoming

    2026

  18. [18]

    Ananthaswamy, Sci

    A. Ananthaswamy, Sci. Am.326(2022), German translation: Stoppuhr f ¨ur die Quantenwelt, Spektrum der Wissenschaft (2022)

    2022

  19. [19]

    C. J. Meier, S¨ uddeutsche Zeitung (2025)

    2025

  20. [20]

    A. J. Bracken and G. F. Melloy, J. Phys. A: Math. Gen.27, 2197 (1994)

    1994

  21. [21]

    G. F. Melloy and A. J. Bracken, Found. Phys.28, 505 (1998)

    1998

  22. [22]

    Bia lynicki-Birula, Z

    I. Bia lynicki-Birula, Z. Bia lynicka-Birula, and S. Au- gustynowicz, J. Phys. A: Math. Theor.55, 255702 (2022)

    2022

  23. [23]

    Several possibilities exist—for instance, the MIT bag, the chiral bag, or the boundary condition(1+𝛽)Ψ =0 adopted in [1]

  24. [24]

    M¨ uller, N

    C. M¨ uller, N. Gr¨ un, and W. Scheid, Phys. Lett. A242, 245 (1998)

    1998

  25. [25]

    re- duced

    The physical constants may be restored with the following dimensionless substitutions:𝜌↦→𝜌/o,𝑧↦→𝑧/o,𝑡↦→𝑐𝑡/o, 𝑣 0 ↦→𝑣 0/(ℏ𝑐), andΨ↦→o 3/2Ψ, whereo=ℏ/(𝑚𝑐)is the “re- duced” Compton wavelength

  26. [26]

    Gallone, inAdvances in Quantum Mechanics, Springer INdAM Series, Vol

    M. Gallone, inAdvances in Quantum Mechanics, Springer INdAM Series, Vol. 18, edited by A. Michelangeli and G. Dell’ Antonio (Springer, Cham, 2017) pp. 169–185

    2017

  27. [27]

    Dolbeault, M

    J. Dolbeault, M. J. Esteban, and E. S ´er´e, J. Spectr. Theory13, 491 (2023)

    2023

  28. [28]

    The↑,↓labels are mnemonic: these states arenoteigenstates of ˆ𝑆𝑧; instead, they are eigenstates of ˆ𝐽𝑧 = ˆ𝐿 𝑧 + ˆ𝑆𝑧

  29. [29]

    Bia lynicki-Birula and Z

    I. Bia lynicki-Birula and Z. Bia lynicka-Birula, Phys. Rev. Res. 2, 038001 (2020)

    2020

  30. [30]

    Zaslow and M

    B. Zaslow and M. E. Zandler, Am. J. Phys.35, 1118 (1967)

    1967

  31. [31]

    C. C. Leary, D. Reeb, and M. G. Raymer, New J. Phys.10, 103022 (2008)

    2008

  32. [32]

    Mathematica, Version 14.3,

    W. R. Inc., “Mathematica, Version 14.3,” Champaign, IL, 2025

    2025

  33. [33]

    G. C. Hegerfeldt, Phys. Rev. D10, 3320 (1974)

    1974

  34. [34]

    Equation (10) may be viewed as a mass-deformed analog of the Cauchy–Riemann system, with the HE (18) in the role of Laplace’s equation; the map (19) correspondingly reconstructs a holomorphic function from a single harmonic seed

  35. [35]

    Bia lynicki-Birula and Z

    I. Bia lynicki-Birula and Z. Bia lynicka-Birula, Phys. Rev. Lett. 118, 114801 (2017)

    2017

  36. [36]

    M. H. Johnson and B. A. Lippmann, Phys. Rev.78, 329 (1950)

    1950

  37. [37]

    A. M. Essin and D. J. Griffiths, Am. J. Phys.74, 109 (2006)

    2006

  38. [38]

    Leviatan, Phys

    A. Leviatan, Phys. Rev. Lett.103, 042502 (2009)

    2009

  39. [39]

    Here ˆ𝑅𝑧 acts as an independent generator, unlike for (14), where it coincides with𝜕 𝑤0

  40. [40]

    I. S. Gradshteyn and I. M. Ryzhik,Table of integrals, series, and products, seventh ed. (Elsevier, New York, 2007)

    2007

  41. [41]

    Andrews, Am

    M. Andrews, Am. J. Phys.76, 1102 (2008)

    2008

  42. [42]

    Mawhin and A

    J. Mawhin and A. Ronveaux, Arch. Hist. Exact Sci.64, 429 (2010)

    2010

  43. [43]

    In cylindrical coordinates,r↦→ −ris implemented via(𝜌, 𝜙, 𝑧) ↦→ (𝜌, 𝜙+𝜋,−𝑧)

  44. [44]

    A. P. Prudnikov, Y. A. Brychkov, and O. I. Marichev,Integrals and Series. Volume 4: Direct Laplace Transforms(Gordon and Breach Science Publishers, New York, 1992). END MATTER Of further interest to specialists, we record two closed-form families that complement the elementary solution (14): a complementary-error-function family from a second weight, and ...

    1992