Effective Caldirola-Kanai Model for Accelerating Twisted Dirac States in Nonuniform Axial Fields
Pith reviewed 2026-05-21 13:27 UTC · model grok-4.3
The pith
The transverse envelope of twisted Dirac particles in longitudinally varying axial fields obeys an effective Caldirola-Kanai Hamiltonian whose solutions follow from an Ermakov mapping of Landau states.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Starting from the Dirac equation and applying controlled spinless and paraxial approximations, the transverse envelope obeys an effective nonstationary Schrödinger equation governed by a Caldirola-Kanai Hamiltonian. The longitudinal energy gain or loss encoded in f(z) = [E_0 - V(z)]^2 - m^2 generates an effective gain or damping rate γ̃(z) = ∂_z f(z) / [2 f(z)] and a z-dependent oscillator frequency ω̃(z) = p_0 Ω(z) / √f(z). Exploiting the Ermakov mapping, the authors obtain a closed-form propagated twisted wave function by transforming the stationary Landau basis; the transverse evolution is controlled by a single scaling function b(z) that satisfies a generalized Ermakov-Pinney equation.
What carries the argument
The Ermakov mapping that establishes unitary equivalence between the Caldirola-Kanai system and a stationary harmonic oscillator, thereby transforming the Landau-level basis into propagating solutions via the scaling function b(z).
If this is right
- For uniform acceleration with vanishing magnetic field the propagated wave function reduces to previously known analytic expressions for accelerating beams.
- For pure solenoid focusing with negligible electric acceleration the solution recovers the standard Landau-level propagation under magnetic focusing.
- The beam radius, phase, and orbital-angular-momentum content at any longitudinal position are determined by evaluating the single scaling function b(z) for the given field profiles.
- The same mapping supplies explicit control over how longitudinal field inhomogeneities reshape the transverse envelope of any twisted state.
Where Pith is reading between the lines
- The framework could be tested by preparing twisted electron or ion beams in a linear accelerator segment with deliberately engineered longitudinal field gradients and measuring the output beam radius against the predicted b(z).
- Relaxing the spinless approximation while retaining the same mapping would yield spin-dependent corrections that might be observable as small polarization-dependent shifts in the focal position.
- The reduction to a single scaling function suggests that similar analytic control may exist for other relativistic wave equations whose transverse dynamics can be cast into time-dependent oscillator form.
Load-bearing premise
The spinless and paraxial approximations remain valid for the entire propagation region even when both the solenoid field and the accelerating electric field vary along the axis, without uncontrolled higher-order corrections from the energy-variation mapping.
What would settle it
Direct numerical solution of the full Dirac equation for a chosen pair of inhomogeneous E_z(z) and B_z(z) profiles, followed by quantitative comparison of the extracted transverse envelope to the analytic form predicted by the scaling function b(z).
Figures
read the original abstract
We study relativistic twisted (orbital-angular-momentum) states of a massive charged particle propagating through an axially symmetric, longitudinally inhomogeneous solenoid field and a co-directed accelerating or decelerating electric field. Starting from the Dirac equation and using controlled spinless and paraxial approximations, we show that the transverse envelope obeys an effective nonstationary Schr\"odinger equation governed by a Caldirola--Kanai Hamiltonian. The longitudinal energy gain or loss encoded in $f(z)=[E_0-V(z)]^2-m^2$ generates an effective gain or damping rate $\widetilde{\gamma}(z)=\partial_z f(z)/[2f(z)]$ and a $z$-dependent oscillator frequency $\widetilde{\omega}(z)=p_0\Omega(z)/\sqrt{f(z)}$. Exploiting the Ermakov mapping (unitary equivalence of Caldirola--Kanai systems), we obtain a closed-form propagated twisted wave function by transforming the stationary Landau basis. The transverse evolution is controlled by a single scaling function $b(z)$ that satisfies a generalized Ermakov--Pinney equation with coefficients determined by $E_z(z)$ and $B_z(z)$. In the limiting cases of uniform acceleration with $B_z=0$ and of solenoid focusing with negligible acceleration, our solution reduces to previously known analytic results, providing a direct bridge to established models.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript derives an effective Caldirola-Kanai Hamiltonian for the transverse envelope of relativistic twisted Dirac states propagating in axially symmetric, longitudinally inhomogeneous solenoid (B_z(z)) and accelerating electric (E_z(z)) fields. Starting from the Dirac equation, it applies spinless and paraxial approximations to obtain a nonstationary 2D Schrödinger equation whose longitudinal energy variation f(z) = [E_0 - V(z)]^2 - m^2 is mapped to an effective damping rate γ̃(z) = ∂_z f(z)/(2f(z)) and z-dependent frequency ω̃(z). The Ermakov-Pinney mapping then yields a closed-form propagated wave function expressed via a single scaling function b(z) that satisfies a generalized Ermakov-Pinney equation with coefficients fixed by E_z(z) and B_z(z). The solution reduces to known analytic results for uniform acceleration (B_z=0) and pure solenoid focusing.
Significance. If the controlled approximations hold, the work supplies a parameter-free analytic bridge between uniform-acceleration and solenoid-focusing regimes for twisted states, with explicit closed-form propagation via the Ermakov mapping and direct reduction to established limits. This is a genuine strength for beam-physics modeling where numerical propagation of Dirac spinors is costly.
major comments (1)
- [derivation of effective Hamiltonian and f(z)→γ̃(z) mapping] The central claim that the spinless/paraxial reduction remains valid under simultaneous longitudinal inhomogeneity in both E_z(z) and B_z(z) is load-bearing for the effective Caldirola-Kanai model. The manuscript asserts the approximations are “controlled” but supplies no explicit bound relating the inhomogeneity length scales (e.g., |dB_z/dz|^{-1} or |dE_z/dz|^{-1}) to the local cyclotron radius or transverse de Broglie wavelength. Without such a quantitative criterion, the neglected transverse-longitudinal coupling and non-paraxial corrections could become O(1) precisely in the regime where the Ermakov mapping is applied.
minor comments (2)
- [Ermakov mapping section] Notation for the scaling function b(z) and the generalized Ermakov-Pinney equation should be cross-referenced to the explicit coefficients involving E_z(z) and B_z(z) for immediate readability.
- [limiting cases] The abstract states that the solution reduces to “previously known analytic results” in the two limiting cases; a brief side-by-side comparison (perhaps in a table or short paragraph) would strengthen the bridge claim.
Simulated Author's Rebuttal
We thank the referee for the careful reading and for identifying the need for explicit validity criteria on the spinless and paraxial approximations. We address the major comment below and will revise the manuscript accordingly.
read point-by-point responses
-
Referee: [derivation of effective Hamiltonian and f(z)→γ̃(z) mapping] The central claim that the spinless/paraxial reduction remains valid under simultaneous longitudinal inhomogeneity in both E_z(z) and B_z(z) is load-bearing for the effective Caldirola-Kanai model. The manuscript asserts the approximations are “controlled” but supplies no explicit bound relating the inhomogeneity length scales (e.g., |dB_z/dz|^{-1} or |dE_z/dz|^{-1}) to the local cyclotron radius or transverse de Broglie wavelength. Without such a quantitative criterion, the neglected transverse-longitudinal coupling and non-paraxial corrections could become O(1) precisely in the regime where the Ermakov mapping is applied.
Authors: We agree that an explicit quantitative criterion is absent from the current manuscript and that its inclusion would strengthen the presentation. In the revised version we will add a short subsection (or appendix paragraph) that derives approximate validity conditions from the neglected terms in the Dirac equation. These conditions take the form L_B, L_E ≫ r_c and L_B, L_E ≫ λ_⊥, where L_B = |dB_z/dz|^{-1}, L_E = |dE_z/dz|^{-1}, r_c = p_⊥/(e B_z) is the local cyclotron radius, and λ_⊥ = 2π ħ / p_⊥ is the transverse de Broglie wavelength. The estimates follow from requiring that the longitudinal derivatives of the vector potential and scalar potential induce corrections smaller than the retained paraxial terms. With these bounds stated, the Ermakov mapping remains applicable inside the indicated regime. The core analytic results and limiting cases are unaffected. revision: yes
Circularity Check
No circularity: derivation starts from Dirac equation with standard approximations and known Ermakov mapping
full rationale
The paper derives the effective Caldirola-Kanai Hamiltonian directly from the Dirac equation via controlled spinless and paraxial approximations, then maps f(z) to the damping rate γ̃(z) and frequency ω̃(z) by explicit differentiation and rescaling. The closed-form solution follows from the established Ermakov-Pinney equation and unitary equivalence to the stationary Landau basis. No step reduces to a self-definition, a fitted parameter renamed as prediction, or a load-bearing self-citation; the limiting cases recover independent prior results without tautology. The derivation chain remains self-contained against the input Dirac equation and standard mathematical tools.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math Dirac equation governs the relativistic charged particle
- domain assumption Spinless and paraxial approximations are controlled and sufficient
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
the transverse envelope obeys an effective nonstationary Schrödinger equation governed by a Caldirola–Kanai Hamiltonian... Exploiting the Ermakov mapping we obtain a closed-form propagated twisted wave function
-
IndisputableMonolith/Foundation/AlphaCoordinateFixation.leanJ_uniquely_calibrated_via_higher_derivative unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
b(z) that satisfies a generalized Ermakov–Pinney equation with coefficients determined by Ez(z) and Bz(z)
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
2c ˜K2 arctanh ˜K2 q 1 + 2 ˜K1z+ ˜K 2 2 z2 −1− ˜K1z ˜K 2 2 z− ˜K1 q 1 + 2 ˜K1z+ ˜K 2 2 z2 −1 2 .(45) 7 We use the connection between hyperbolic tangent and natural logarithm in form arctanhx= 1 2 ln 1 +x 1−x ,|x|<1 and obtain in previous notation b(z) = 1√c vuut1 + (q cb2 0 −1 + sign(b ′
-
[2]
The corresponding differential equation can be analytically solved through the presented algorithm
cp0z V(z) ln " E0 +p 0 E0 −p 0 p f(z)−p 0 +V(z)p f(z)−p 0 −V(z) #)2 .(46) The same logic could be applied to the evolution in the homogeneous magnetic field and the inhomogeneous electric field. The corresponding differential equation can be analytically solved through the presented algorithm. II. HOMOGENEOUS ELECTRIC FIELD CROSS-CHECK In Ref. [25] the au...
- [3]
- [4]
- [5]
- [6]
- [7]
-
[8]
B. J. McMorran, A. Agrawal, I. M. Anderson, A. A. Herz- ing, H. J. Lezec, J. J. McClelland, and J. Unguris, Science 331, 192 (2011)
work page 2011
-
[9]
T. Schachinger, S. L¨ offler, M. St¨ oger-Pollach, and P. Schattschneider, Ultramicroscopy158, 17 (2015)
work page 2015
-
[10]
E. Mafakheri, A. H. Tavabi, P.-H. Lu, R. Balboni, F. Ven- turi, C. Menozzi, G. C. Gazzadi, S. Frabboni, A. Sit, R. E. Dunin-Borkowski, E. Karimi, and V. Grillo, Ap- plied Physics Letters110, 093113 (2017)
work page 2017
-
[11]
I. P. Ivanov, Phys. Rev. D83, 093001 (2011)
work page 2011
-
[12]
I. P. Ivanov, Progress in Particle and Nuclear Physics 127, 103987 (2022)
work page 2022
-
[13]
D. V. Karlovets, S. S. Baturin, G. Geloni, G. K. Sizykh, and V. G. Serbo, Eur. Phys. J. C83, 372 (2023)
work page 2023
-
[14]
A. D. Chaikovskaia, D. V. Karlovets, and V. G. Serbo, Phys. Rev. A109, 012222 (2024)
work page 2024
- [15]
-
[16]
N. Sheremet, A. Chaikovskaia, D. Grosman, and D. Karlovets, Phys. Rev. A111, 052810 (2025)
work page 2025
-
[17]
A. S. Dyatlov, M. A. Nozdrin, A. N. Sergeev, N. E. Sheremet, S. S. Stafeev, and D. V. Karlovets, Appl. Opt. 64, 10776 (2025)
work page 2025
-
[18]
A. Dyatlov, V. Bleko, K. Cherepanov, V. Kobets, M. Martyanov, M. Nozdrin, A. Sergeev, N. Sheremet, A. Zhemchugov, and D. Karlovets, in2024 International Conference Laser Optics (ICLO)(2024) pp. 438–438
work page 2024
-
[19]
A. S. Dyatlov, A. V. Afanasyev, V. V. Kobets, A. E. Levichev, M. V. Maksimov, D. A. Nikiforov, M. A. Noz- drin, K. Popov, K. A. Sibiryakova, K. E. Yunenko, and D. V. Karlovets, Classical and quantum beam dynam- ics simulation of the RF photoinjector test bench (2025), arXiv:2509.00732 [physics.acc-ph]
-
[20]
X. Liu, Q. Meng, Z. Yang, W. Ma, L. Lu, P. Zhang, and L. Zou, Journal of Applied Physics138, 124303 (2025)
work page 2025
-
[21]
L. Zou, P. Zhang, and A. J. Silenko, Phys. Rev. A103, L010201 (2021)
work page 2021
-
[22]
G. Guzzinati, P. Schattschneider, K. Y. Bliokh, F. Nori, and J. Verbeeck, Phys. Rev. Lett.110, 093601 (2013)
work page 2013
-
[23]
G. K. Sizykh, A. D. Chaikovskaia, D. V. Grosman, I. I. Pavlov, and D. V. Karlovets, Progress of Theoretical and Experimental Physics2024, 053A02 (2024)
work page 2024
-
[24]
G. K. Sizykh, A. D. Chaikovskaia, D. V. Grosman, I. I. Pavlov, and D. V. Karlovets, Phys. Rev. A109, L040201 (2024)
work page 2024
- [25]
-
[26]
N. V. Filina and S. S. Baturin, Phys. Rev. A113, L021302 (2026)
work page 2026
-
[27]
Q. Meng, Z. Huang, X. Liu, W. Ma, Z. Yang, L. Lu, A. J. Silenko, P. Zhang, and L. Zou, Phys. Rev. Res.7, 043213 (2025)
work page 2025
- [28]
-
[29]
Kanai, Progress of Theoretical Physics3, 440 (1948)
E. Kanai, Progress of Theoretical Physics3, 440 (1948)
work page 1948
-
[30]
F. F. L´ opez-Ruiz and J. Guerrero, Journal of Physics: Conference Series538, 012015 (2014)
work page 2014
-
[31]
N. V. Filina and S. S. Baturin, Phys. Rev. A108, 012219 (2023)
work page 2023
-
[32]
H. R. Lewis, Phys. Rev. Lett.18, 510 (1967)
work page 1967
-
[33]
H. R. Lewis, Phys. Rev. Lett.18, 636 (1967)
work page 1967
-
[34]
H. R. Lewis and W. B. Riesenfeld, J. Math. Phys.10, 1458 (1969)
work page 1969
- [35]
-
[36]
Q. Meng, X. Liu, W. Ma, Z. Yang, L. Lu, A. J. Silenko, 9 P. Zhang, and L. Zou, Phys. Rev. Res.7, 023306 (2025)
work page 2025
- [37]
- [38]
-
[39]
See Supplemental Material at [URL] for details on the (i) an analytical solution of the Ermakov–Pinney equa- tion for an arbitrary longitudinal electric field in the ab- sence of a magnetic field; (ii) the homogeneous-electric- field limit and a cross-check against Ref. [25]; and (iii) an explicit construction of the Ermakov–Pinney solution from two linea...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.