Asymptotic Momentum of Dirac Particles in One Space Dimension

Abigail Perryman; A. Shadi Tahvildar-Zadeh; Kabir Narayanan

REVIEW 3 major objections 5 minor 1 cited by

For a free Dirac particle with Gaussian wave-packet data, momentum k survives while energy becomes either +E or −E depending on the particle's initial position, sending it at constant speed to the left or right.

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 · deepseek-v4-flash

2026-08-03 14:06 UTC pith:L3Z2BH5T

load-bearing objection Real SPA theorem, but the advertised long-time result for fixed parameters isn't what the proof shows. the 3 major comments →

arxiv 2512.21423 v2 pith:L3Z2BH5T submitted 2025-12-24 math-ph math.APmath.MPquant-ph

Asymptotic Momentum of Dirac Particles in One Space Dimension

Kabir Narayanan , Abigail Perryman , A. Shadi Tahvildar-Zadeh This is my paper

classification math-ph math.APmath.MPquant-ph MSC 35Q4181Q0581P05 PACS 03.65.Pm03.65.Sq

keywords Dirac equationBohmian mechanicsstationary phase approximationGaussian wave packetasymptotic momentumnegative energyZitterbewegungclassical limit

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 proves that a massive spin-½ particle in one dimension, starting as a Gaussian wave packet with non-zero mean momentum k and positive mean energy E=√(m²+k²), asymptotically separates into two non-overlapping packets: one with energy E moving in the direction of k, and one with energy −E moving oppositely. Because of this separation, a Bohmian trajectory guided by the wave function eventually straightens into uniform motion: which direction and which energy sign depend on the initial position, so the particle's asymptotic momentum equals k for almost all initial data. The result matters because it grounds the common physics practice of treating free particles as plane waves with definite energy-momentum (as in Compton scattering) in a rigorous wave-packet picture, and because it shows how the 'negative energy' part of a Dirac wave function acts on trajectories.

Core claim

The central claim is Theorem 2.1: for Gaussian initial data of the form (37), the Dirac evolution approaches, as t→∞, the superposition of two wave packets with the same momentum k and opposite energies ±E that travel in opposite directions. Consequently the worldline of a typical particle — one with initial position outside a small exceptional set — eventually becomes a straight line with constant velocity, and the sign of the asymptotic Bohmian energy is fixed by which side of a critical initial position the particle starts on. The paper establishes this by computing a stationary-phase approximation to the exact integral representation of the solution and by constructing explicit barrier c

What carries the argument

The proof rests on (i) an exact integral representation of the Dirac solution (14) with Bessel kernels; (ii) a two-dimensional stationary-phase approximation for the oscillatory integrals, valid in the classical limit ω≫1, whose leading terms (52)–(53) are exactly the two counter-propagating Gaussian packets; (iii) the rigorous error bound (54) for the remainder; and (iv) barrier hyperbolas s±(t)=σ²C±/(v₀t) in the guiding ODE that confine trajectories so that each one is eventually guided by only one packet. The key identity is the cancellation in Lemma 3.1 that removes the free transport term, making ψ+ an oscillatory integral amenable to SPA.

Load-bearing premise

The load-bearing premise is that the very large Compton-frequency (classical) limit and the very large time limit can be interchanged — that is, for fixed k and m the two packets separate and their overlap vanishes as t→∞, even though the only controlled error bound is for ω→∞ on a bounded macroscopic time interval.

What would settle it

Compute the exact solution (14) for fixed parameters, say m=k=1, and at a sequence of growing times T measure the overlap integral ∫|ψ−(T,s)||ψ+(T,s)|ds or the set of positions where the velocity field v=J/ρ differs from ±v₀. The theorem predicts this overlap and the exceptional set both tend to zero in probability as T→∞; if either saturates at a positive value at fixed ω, the claim as worded fails.

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

If this is right

For Gaussian positive-energy initial data, the asymptotic Bohmian momentum equals the conserved quantum expectation of the momentum operator, k, for typical trajectories.
Negative-energy packets transport opposite to their momentum, so a particle starting on that side recedes with velocity −v₀ while carrying momentum +k.
The supports of the two packets become almost disjoint at large times, so interference and zitterbewegung cease after the packets separate.
The result justifies treating the incoming electron as a plane wave with dispersion E=√(m²+k²) in scattering calculations such as Compton's.
By the same stationary-phase formula applied to the general superposition (46), the phenomena extend to any mixture of the two energy eigenstates with fixed k.

Where Pith is reading between the lines

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

The proof is carried out in the classical limit ω≫1 on a bounded macroscopic time interval; a fully rigorous reading of 't→∞' at fixed k and m would require a separate uniformity argument in t that the stated error bound (54) does not supply.
Because the paper's numerical experiments use initial spinors not constrained by the energy-eigenfunction condition, the result likely extends to arbitrary Gaussian two-component data, with a critical initial position separating the two asymptotic behaviours.
The deterministic sign flip at the bifurcation point suggests a hidden-variable threshold: the asymptotic energy sign is fixed by initial position, so 'measurement' of energy here corresponds to a pre-existing property of the trajectory.
The massless limit is left open; since Bohmian momentum p=m cotΘ secΩ diverges as m→0, the plane-wave justification for photons would need a separate limiting argument, exactly the issue the paper flags for future work.

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

Real SPA theorem, but the advertised long-time result for fixed parameters isn't what the proof shows.

read the letter

Colleague,

The paper does something real: it gives a rigorous stationary-phase approximation for the 1D free Dirac equation with Gaussian initial data and proves an explicit error bound. That is the genuinely new piece. The split of a positive-energy packet into two packets with the same momentum but opposite energies and velocities is known qualitatively, and the authors cite the right Bohmian literature. The barrier argument for the trajectories is also a nice touch. The numerics are illustrative and consistent.

The soft spot is exactly where the stress-test lands. Theorem 2.1 is stated for fixed k and m, as t→∞. The proof, however, lives in the classical limit of Section 3: macroscopic length L, ω = L/λ_CWE ≫ 1, and initial momentum k=ωp0, mass m=ω. The SPA error bound (54) is A|p0|^5 e^{Bt/σ}/√ω, which is small for large ω on a bounded macroscopic time interval [|v0|T/2,T]. It does not become small as t→∞ with ω fixed; the exponential factor grows. The trajectory argument in Section 3.3 inherits that: it shows that for each T and sufficiently large ω, typical trajectories reach constant velocity on [t0,T], but that is not the same as an asymptotic statement for a fixed electron over infinite time.

So the central claim as worded overclaims. That is a real issue, but it is fixable. The authors should restate the main theorem as a classical-limit theorem: for large ω and times up to a fixed macroscopic T, the solution and the Bohmian trajectories are approximated by the two-packet SPA with controlled error. That is what the proof actually supports. The missing numerical code and the unspecified constants A and B are minor by comparison.

The paper deserves a serious referee. The technical machinery is sound, and the gap between statement and proof is clear and addressable rather than a hidden circularity. I would send it to review with the expectation of a major revision. For a reading group, it is a good case study in matching a theorem statement to the actual limit being taken.

Best

Referee Report

3 major / 5 minor

Summary. The paper studies Bohmian trajectories for a single massive Dirac particle in one space dimension, with Gaussian wave-packet initial data. It claims (Theorem 2.1) that for fixed momentum k and mass m, as t→∞ the solution splits into two packets with the same momentum k and energies ±E = ±√(m²+k²), traveling in opposite directions, and that typical Bohmian trajectories consequently become straight lines with asymptotic speed ±v0 = ±k/E, the sign determined by the initial position. The proof is based on a stationary-phase approximation (SPA) of the exact integral representation of the Dirac solution. A rigorous SPA error bound is derived using Chako's method, and the trajectory statement is approached through barrier curves for the SPA velocity field and a typicality/equivariance argument. Numerical experiments in Appendix B illustrate the predicted behavior.

Significance. If the claims were fully established, the paper would provide a rigorous connection between Bohmian and operator notions of energy and momentum for free Dirac particles, and would offer some justification for using plane-wave dispersion relations in simplified scattering models. The paper contains genuine technical contributions: the SPA computation is explicit (critical points (0,y±), Hessians, phases, stationary-point contributions), and the error analysis via Chako's method is careful, including bookkeeping of the dependence on p0, σ, and the exponential-in-t factor. The numerical experiments are a useful sanity check. However, as detailed below, the main theorem as stated is stronger than what the proof actually establishes, and the trajectory argument has a significant gap.

major comments (3)

[Theorem 2.1 vs. §3, Eq. (54)] Theorem 2.1 is stated for fixed k, m and asserts behavior as t→∞. The proof, however, is carried out in the classical limit of §3: m=ω, k=ωp0, and Theorem 3.2 proves the SPA only for t∈[|v0|T/2,T] with the error bound |R±|≤A|p0|^5 e^{Bt/σ}/√ω. This error is made small by taking ω→∞ with T fixed; the factor e^{Bt/σ} grows in t and cannot control t→∞ at fixed physical parameters. Thus the proof establishes a high-mass/high-momentum (classical-limit) statement on bounded macroscopic time, not the fixed-parameter long-time claim of Theorem 2.1 and the abstract. Please either prove the t→∞ statement by different estimates or restate the theorem/abstract as a classical-limit result.
[§3.3, proof of Theorem 2.1] The trajectory conclusion is not derived from the estimates. The barriers s±(t)=σ²C±/(v0t) tend to 0 as t grows, whereas the packet centers are at ±v0t. Being to the right of H+ (or of Q+) only gives s·t ≥ σ²C+; it does not imply that |φ+|/|φ−| is small or that the SPA velocity v(U) is close to v0. The sentence that “it will thus be (almost) entirely guided by φ− and will have a constant velocity v0 by the time it reaches t=T” is asserted, not proved. Lemma 3.4 controls |q−q̃| once such velocity closeness is known, but it does not establish that the trajectory has entered the region where one Gaussian dominates. A quantitative estimate showing that q(t) reaches s(t)≫σ²/(v0t) (or an equivalent condition) and stays there is needed.
[Eqs. (28), (35), (43), (50)] The connection between the initial data in Theorem 2.1 and the initial data actually analyzed is not explicit. Theorem 2.1 uses the spinor (k/E,m/E) in (37), while Theorem 3.2 analyzes (0,1)^T in (50), and §3.1 claims a WLOG reduction using the basis (35),(36). However, (35) appears to write components as (cosΘ,sinΘ) rather than the half-angle components (cos(Θ/2),sin(Θ/2)) of the eigenfunctions in (28). The two are not proportional for finite k. Since the SPA formulas in Corollary 3.2.1 are expressed in half-angle variables, the paper needs to show explicitly how (37) is represented in that basis and how Theorem 2.1 follows from Theorem 3.2/Corollary 3.2.1. As written, the hypotheses of the theorem and the proof setup are not visibly matched.

minor comments (5)

[Abstract and Introduction] The abstract says the wave function “becomes locally a plane wave,” while the main theorem states a superposition of two packets. Please align the wording. Also “at large times” should be qualified as “in the classical limit, on bounded macroscopic time intervals” if the theorem is restated accordingly.
[§3.1, Eq. (45)] The scale-separation assumption 1/(ωp0)≪σ≪1 is introduced as “we may also assume,” but Theorem 3.2 is stated for arbitrary p0,σ,T. Please clarify whether (45) is actually used in the proof or is only motivation.
[§3.3, Eq. (65)-(69)] The notation C+ and C− is introduced in (69); it would help to state explicitly that C+ can vanish at Θ0=π/4 and that the hyperbola degenerates, as done in the text, but the proof of Theorem 2.1 should handle this case consistently (no separate treatment is given).
[Appendix B] The numerical figures are informative but the text does not describe the numerical integration method or the discretization parameters. A brief statement of the numerical scheme and error control would strengthen the appendix.
[General] There are minor typographical issues (e.g., “deBroglie,” “the above” repeated, occasional missing articles). They do not affect the mathematics.

Circularity Check

0 steps flagged

No circularity: the asymptotic packet decomposition is obtained by stationary-phase evaluation of the exact integral representation, with the error bound estimating the remainder; no fitted parameter is renamed a prediction. The only self-citations are motivational. A fixed-parameter t→∞ wording overreaches the ω→∞, bounded-T proof, but that is a correctness gap, not circularity.

full rationale

The asymptotic formulas (52)-(53) are not inserted as an ansatz: they are derived from the exact solution (47)-(48) by SPA, with the phases e^{iω(p0s∓E0t)} arising from the two critical points y±=p0±√(1+p0²) of the phase function ϕ(x,y) in (79). The result is therefore a rigorous consequence of the free Dirac evolution, not an equivalence with the input by construction. The error bound (54), |R±|≤A|p0|^5 e^{Bt/σ}/√ω, is an estimate for the difference between the exact and approximate solutions; it is not a fitted value and does not presuppose the theorem. The initial data (37) contains k and E only as parameters; the theorem shows that typical trajectories asymptotically move with the corresponding v0=±k/E, which is nontrivial because the Bohmian velocity field is time-dependent and the wave packet splits. The only self-citation, [7] (Kiessling, Lienert, and one of the present authors), is used in the introduction as motivation and for numerical context; it is not load-bearing for Theorem 2.1 or Theorem 3.2. The proof is self-contained, relying on the exact solution (14)/(41)-(42), Lemma 3.1, and Chako's external SPA estimates. One limitation is flagged explicitly by the paper itself: "The proof of the above theorem relies on the Stationary Phase Approximation (SPA) method ... Thus to apply SPA we need to study the limit ω→∞" (Section 2, after Theorem 2.1). Consequently, Theorem 3.2 and the trajectory estimate Lemma 3.4 are proven only on t∈[|v0|T/2,T] for ω>Ω(p0,T,σ), not for fixed m,k as t→∞. This is an overclaim/inference gap about exchanging large time with the classical limit, not a circularity: the asymptotic momentum is still derived from the dynamics and is not assumed equal to the input. Thus the circularity score is 0.

Axiom & Free-Parameter Ledger

2 free parameters · 6 axioms · 0 invented entities

No new entities are introduced. The only free parameters are physical initial data (mean momentum and width), not fitted constants. The main interpretive load is Bohmian mechanics and the classical-limit scaling.

free parameters (2)

k (p0 after rescaling)
Initial mean momentum; the asymptotic result states it is preserved. It is an input to the problem, not fitted to the conclusion.
sigma (Gaussian width)
Initial position spread; constrained to satisfy 1/(ωp0)≪σ≪1 in the classical-limit proof.

axioms (6)

domain assumption Free Dirac equation with gamma matrix representation (10)-(11) governs a spin-half particle in 1D
Starting physical model; standard relativistic quantum mechanics.
domain assumption Bohmian guiding equation dQ/dt = (|ψ−|²−|ψ+|²)/(|ψ−|²+|ψ+|²) and typicality/quantum equilibrium (Definition 3.3)
Interpretive framework of the paper; without it, 'trajectories' have no meaning.
domain assumption Bohmian local momentum and energy definitions p=m cotΘ secΩ, E=m cscΘ secΩ (22)
These define the asymptotic momentum/energy claims; from Vigier [9].
standard math Stationary phase approximation error estimates of Chako [1] apply to oscillatory integrals with parameters (s,t,p0,σ), with phase Hessian nondegenerate and amplitude analytic in a neighborhood of critical points
Basis for the claimed rigorous error bound.
domain assumption Separation of scales 1/(ωp0)≪σ≪1 and k=ωp0 with ω=L/λ_CWE≫1
The proof is in the classical limit; not satisfied for fixed k as ω→∞.
standard math Uniqueness, existence, and finite speed of propagation for the Dirac equation and the ODE trajectories
Used in barrier/typicality argument.

pith-pipeline@v1.3.0-alltime-deepseek · 19603 in / 16339 out tokens · 155204 ms · 2026-08-03T14:06:29.386266+00:00 · methodology

0 comments

read the original abstract

We analyze the trajectories of a massive particle in one space dimension whose motion is guided by a spin-half wave function that evolves according to the free Dirac equation, with its initial wave function being a Gaussian wave packet with a nonzero expected value of momentum $k$. We prove that at large times, the wave function is approximately equal to the superposition of two wave packets traveling in opposite directions, which results in trajectories with approximately constant asymptotic momentum $k$ and asymptotic energy $\pm c^2\sqrt{m^2+k^2}$, with $m$ the rest mass of the particle and $c$ the speed of light. The sign of the asymptotic energy is determined by the initial position of the particle. Particles with negative energy will have an asymptotic velocity that is in the opposite direction of their momentum. The proof uses the stationary phase approximation method, for which we establish a rigorous error bound.

Figures

Figures reproduced from arXiv: 2512.21423 by Abigail Perryman, A. Shadi Tahvildar-Zadeh, Kabir Narayanan.

**Figure 1.** Figure 1: The curves y = y0(x) and y = B±(x) for Θ0 = π 8 (on the left) and Θ0 = 3π 8 (on the right). Hence the velocity, and therefore the momentum, of the particle becomes a constant, as advertised. To find the B± curves, we first find y0(x) such that F(x, y0(x)) = 0. Let X := e xy > 0. Then F(x, y) = cos Θ0 X2 + η 2 (X2 + 2 tan Θ0 cos(aωx)ηX − η 2 ) so that if we define y0(x) = 1 x ln tan Θ0 2 tan Θ0 cos aωx + p … view at source ↗

**Figure 2.** Figure 2: Barrier construction in the case Θ0 < π/4 stationary phase approximation v (U) are sufficiently close, so that the trajectories guided by them are uniformly close for t ∈ [t0, T]. Lemma 3.4. Given ϵ0 > 0, p0 ̸= 0, T > 0, and σ > 0 there exists Ω0 = Ω0(ϵ0, p0, T, σ) > 0 such that if ω > Ω0, for the trajectories s = q(t) and s = ˜q(t) defined in the above, we have |q(t) − q˜(t)| < ϵ0 ∀t ∈ [t1, T]. (70) Proof… view at source ↗

**Figure 3.** Figure 3: 50 electron trajectories. Each row demonstrates the effect of varying one of the initial [PITH_FULL_IMAGE:figures/full_fig_p025_3.png] view at source ↗

**Figure 4.** Figure 4: The Bloch sphere representation of a single electron trajectory over 8 seconds ( [PITH_FULL_IMAGE:figures/full_fig_p026_4.png] view at source ↗

**Figure 5.** Figure 5: Numerical vs. predicted values of Bloch variables for 50 trajectories after 6 seconds, for [PITH_FULL_IMAGE:figures/full_fig_p026_5.png] view at source ↗

discussion (0)

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Scattering Cross Section Formula Derived From Macroscopic Model of Detectors
quant-ph 2026-01 conditional novelty 7.0

The usual quantum scattering cross-section formula is derived as the actual detection distribution from imaginary-potential and repeated-measurement detector models in the far-field limit.

Reference graph

Works this paper leans on

9 extracted references · cited by 1 Pith paper

[1]

4, 372–422

Nicholas Chako,Asymptotic expansions of double and multiple integrals occurring in diffraction theory, IMA Journal of Applied Mathematics1(1965), no. 4, 372–422

1965
[2]

Compton,A quantum theory of the scattering of x-rays by light elements, Phys

Arthur H. Compton,A quantum theory of the scattering of x-rays by light elements, Phys. Rev. 21(1923), 483–502

1923
[3]

Series A, Containing Papers of a Mathematical and Physical Character117 (1928), no

Paul Adrien Maurice Dirac,The quantum theory of the electron, Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character117 (1928), no. 778, 610–624

1928
[4]

D¨ urr, S

D. D¨ urr, S. Goldstein, K. M¨ unch-Berndl, and N. Zangh ` ı,Hypersurface Bohm– Dirac models, Phys. Rev. A60(1999), 2729–2736

1999
[5]

5, 843–907

Detlef D¨ urr, Sheldon Goldstein, and Nino Zangh ` ı,Quantum equilibrium and the origin of ab- solute uncertainty, Journal of Statistical Physics67(1992), no. 5, 843–907

1992
[6]

Peter R Holland,The quantum theory of motion: an account of the de Broglie-Bohm causal interpretation of quantum mechanics, Cambridge university press, 1995

1995
[7]

Kiessling, Matthias Lienert, and A

Michael K.-H. Kiessling, Matthias Lienert, and A. Shadi Tahvildar-Zadeh,A Lorentz-covariant interacting electron–photon system in one space dimension, Letters in Mathematical Physics 110(2020), no. 12, 3153–3195

2020
[8]

thesis, LMU, 2010

Sarah R¨ omer,The classical limit of Bohmian mechanics, Ph.D. thesis, LMU, 2010

2010
[9]

Sci.(Paris)235(1952), 1107–1109

Jean-Pierre Vigier,Forces s’exer¸ cant sur les lignes de courant usuelles des particules de spin 0, 1/2, et 1 en th´ eorie de l’onde pilote, CR Acad. Sci.(Paris)235(1952), 1107–1109. 24 Figure 3: 50 electron trajectories. Each row demonstrates the effect of varying one of the initial parameters Θ0,k 0,m el, andσ. 25 Figure 4: The Bloch sphere representati...

1952

[1] [1]

4, 372–422

Nicholas Chako,Asymptotic expansions of double and multiple integrals occurring in diffraction theory, IMA Journal of Applied Mathematics1(1965), no. 4, 372–422

1965

[2] [2]

Compton,A quantum theory of the scattering of x-rays by light elements, Phys

Arthur H. Compton,A quantum theory of the scattering of x-rays by light elements, Phys. Rev. 21(1923), 483–502

1923

[3] [3]

Series A, Containing Papers of a Mathematical and Physical Character117 (1928), no

Paul Adrien Maurice Dirac,The quantum theory of the electron, Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character117 (1928), no. 778, 610–624

1928

[4] [4]

D¨ urr, S

D. D¨ urr, S. Goldstein, K. M¨ unch-Berndl, and N. Zangh ` ı,Hypersurface Bohm– Dirac models, Phys. Rev. A60(1999), 2729–2736

1999

[5] [5]

5, 843–907

Detlef D¨ urr, Sheldon Goldstein, and Nino Zangh ` ı,Quantum equilibrium and the origin of ab- solute uncertainty, Journal of Statistical Physics67(1992), no. 5, 843–907

1992

[6] [6]

Peter R Holland,The quantum theory of motion: an account of the de Broglie-Bohm causal interpretation of quantum mechanics, Cambridge university press, 1995

1995

[7] [7]

Kiessling, Matthias Lienert, and A

Michael K.-H. Kiessling, Matthias Lienert, and A. Shadi Tahvildar-Zadeh,A Lorentz-covariant interacting electron–photon system in one space dimension, Letters in Mathematical Physics 110(2020), no. 12, 3153–3195

2020

[8] [8]

thesis, LMU, 2010

Sarah R¨ omer,The classical limit of Bohmian mechanics, Ph.D. thesis, LMU, 2010

2010

[9] [9]

Sci.(Paris)235(1952), 1107–1109

Jean-Pierre Vigier,Forces s’exer¸ cant sur les lignes de courant usuelles des particules de spin 0, 1/2, et 1 en th´ eorie de l’onde pilote, CR Acad. Sci.(Paris)235(1952), 1107–1109. 24 Figure 3: 50 electron trajectories. Each row demonstrates the effect of varying one of the initial parameters Θ0,k 0,m el, andσ. 25 Figure 4: The Bloch sphere representati...

1952