Bound electron states in a charged chain within the Dirac description

Alexander Eremko; Larissa Brizhik; Vadim Loktev

REVIEW 3 major objections 3 minor 15 references

Bound electron states in a charged chain within the Dirac description

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

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

Pith's one-line read Paper derives exact analytical three-dimensional Dirac bound states for an electron in the Coulomb field of a chain of positive ions.

desk verdict A careful solution of the Dirac equation for an effective 2D Coulomb potential is presented as an exact solution for a charged chain, but the central exactness claim is unsupported because the coupling constant is free and the true potential is logarithmic. read the letter →

arxiv 2411.14212 v1 pith:RS6V4JLB submitted 2024-11-21 cond-mat.mes-hall

Alexander Eremko , Larissa Brizhik , Vadim Loktev This is my paper · ORCID

classification cond-mat.mes-hall

keywords Diracequationspinorinvariantboundelectronstatelow-dimensionalsystemchargedchainspin-orbitcouplingLaguerrepolynomialsCoulombpotential

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 aims to provide the first exact analytical description of three-dimensional bound electron states in the Coulomb field of an infinite chain of positively charged ions, solved within the Dirac equation rather than the nonrelativistic Schrödinger equation. Using the chain's symmetries and a newly found spinor invariant, the authors obtain orthonormal bispinors whose radial parts are Laguerre polynomials and whose energies are $E_{k,n,M}=E_k\,\Delta_{n,M}$ with $\Delta_{n,M}$ given by Eq. (32). The point of the exercise is that spin and longitudinal motion become coupled automatically, so spin-orbit interaction is a consequence of the Dirac structure instead of an added term. A sympathetic reader would care because these states are the natural zero-order building blocks for a relativistic theory of self-trapped Davydov solitons on molecular chains.

What carries the argument

The central object is the newly found spinor invariant (45), an operator that commutes with the Dirac Hamiltonian for the charged chain together with the longitudinal momentum $\hat p_z$, the total angular momentum $\hat J_z$, and the spin-polarization operator $\hat S_z$. It supplies the fourth quantum label and accounts for the accidental spin degeneracy of the energy spectrum, just as the Johnson–Lippmann invariant explains degeneracies in relativistic hydrogen. The solution mechanism is the invariance algebra of the Dirac equation: after imposing the joint eigenstates of these operators, the radial equations reduce to the hypergeometric equation, and termination of the hypergeometric series gives Laguerre polynomials and the quantization condition (28). The effective potential $V_{\mathrm{eff}}=C_{\mathrm{eff}}e^2Z_v/r_\perp$ is the input that makes the reduction exact, with $\gamma_M=\sqrt{M^2-C_{\mathrm{eff}}^2Z_v^2\alpha^2}$ controlling the radial falloff $r_\perp^{\gamma_M-1/2}$ of the bound states.

What would settle it

A numerical solution of the radial Dirac equation with the true potential $V_0(r_\perp)=-(2e^2Z_v/a)\ln v(r_\perp)$ from Eq. (6) would settle the claim: if the bound-state energies obtained numerically cannot be matched by $E_k\Delta_{n,M}$ for any fixed choice of $C_{\mathrm{eff}}$, then the effective-potential replacement is not valid and the paper's central claim fails.

Watch

Extended reading notes

Core claim

The paper claims to obtain, for the first time, exact analytical expressions for three-dimensional bound electron states in the Coulomb field of a chain of positively charged ions, using the Dirac equation. Concretely, the eigen-bispinors (40)–(44) have the form of a plane wave $e^{i(kz+M\varphi)}$ along the chain times $e^{-\xi/2}r_\perp^{\gamma_M-1/2}$ times combinations of normalized Laguerre polynomials, with $\xi=2\kappa_\perp r_\perp$, $\gamma_M=\sqrt{M^2-C_{\mathrm{eff}}^2Z_v^2\alpha^2}$, and $\kappa_\perp$ given by Eq. (33). The bound-state energy is $E=E_k\Delta_{n,M}$, Eqs. (31)–(32), where $E_k=\sqrt{m^2c^4+c^2\hbar^2k^2}$ and $\Delta_{n,M}$ depends on the radial quantum number $n$, the half-integer angular momentum $M$, and the fine structure constant $\alpha$. The derivation works with an effective potential $V_{\mathrm{eff}}=C_{\mathrm{eff}}e^2Z_v/r_\perp$ replacing the true chain potential for $r_\perp>a$. The paper also constructs a new spinor invariant, Eq. (45), which explains why states with the same $n,k,M$ are degenerate in the spin label $\sigma$, in analogy with the Johnson–Lippmann invariant of the relativistic Kepler problem. In this description spin and longitudinal propagation are entangled automatically, so no hand-written spin-orbit term is needed.

Load-bearing premise

The calculation becomes exact only after replacing the actual chain potential with a fitted $1/r$ potential whose constant $C_{\mathrm{eff}}$ is not determined from the chain itself; if the fitted constant cannot be fixed from the actual potential, the exactness claim collapses.

Editorial extensions

If this is right

Bound states of an electron on a charged chain form a discrete family labelled by longitudinal momentum $k$, half-integer angular momentum $M$, radial quantum number $n$, and spin $\sigma$, with energy $E_k\Delta_{n,M}$.
The transverse decay length $\kappa_\perp^{-1}$ depends on the longitudinal energy through $E_k^2-E^2$, so the confinement of the wavefunction and the electron's motion along the chain are coupled.
Spin projection and longitudinal motion enter the bispinor through $\sigma E_k$ and $c\hbar k$ in relations (15)–(16), so the Dirac equation alone generates spin-orbit-type coupling.
The spin degeneracy of $E_{k,n,M}$ is exact and is explained by the new invariant (45), giving a symmetry-based explanation rather than a numerical coincidence.
These wavefunctions are the stated starting point for adding the chain's periodicity and lattice deformation, i.e. for constructing a relativistic Davydov soliton.

Reading between the lines

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

If $C_{\mathrm{eff}}$ is fixed by calibrating one bound-state energy against a numerical or experimental value, the remaining states become quantitative predictions for electrons on charged nanowires or molecular chains; the paper does not perform this calibration.
The invariant (45) is likely to survive as an approximate or exact symmetry when weak periodicity is added, which would make the transverse eigenfunctions a useful basis for band-structure and polaron calculations.
Because $\Delta_{n,M}$ depends on $Z_v\alpha$ and on $M^2$, the spectrum should exhibit observable relativistic fine-structure splittings of transverse levels; measuring the level ordering could discriminate this Dirac description from a nonrelativistic one.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the letter

The algebra in this paper is mostly sound, but the central claim is not. The explicit bound states (40)-(44) are exact for the Dirac equation with the 2D Coulomb potential V_eff = C_eff e^2 Z_v / r_perp, where C_eff is a free parameter. They are not exact states for the actual chain potential V0(r_perp) of Eq. (6), which is logarithmic at small r and only asymptotically 1/r. C_eff is never fixed by V0, so the spectrum (31)-(32) is a function of an undetermined coupling, not a prediction for the charged chain. The stress-test note is right on this point.

What is genuinely useful: the reduction to 1D radial equations using the constants of motion is clean, the Laguerre solution is explicit and includes normalization, and the two spin states are worked out symmetrically. If you work on Dirac particles in a transverse 2D Coulomb-like field with conserved axial momentum, this is a convenient reference for wavefunctions. The claimed new spinor invariant (45) is intriguing but simply asserted; no proof or verification is offered, so it currently reads as a conjecture.

The main problem is overstatement. The abstract and conclusions claim "exact analytical expressions" for the "realistic Coulomb potential (6)", but that is not accurate. The effective potential is introduced with an undetermined C_eff, and the true V0 is not 1/r. Also, the solved equation is essentially the 2D Dirac-Coulomb problem, which has known solutions; the "for the first time" claim is too strong without citing that literature. The unproven invariant is another gap, but the potential replacement alone is enough to break the exactness claim.

Who is this for? Someone working on Dirac equations in low-dimensional systems who wants explicit cylinder-symmetric spinor states in an effective potential. It is not a paper to rely on for the actual chain potential. It deserves referee time because the derivation is careful and the invariant, if proven, could be a real result. As written it should be rejected, but the authors could fix it by either determining C_eff from V0 (e.g., by a matching procedure) or explicitly presenting this as an effective-model solution, and by proving or properly citing the invariant.

Recommendation: send to peer review, but expect the exactness claims to be substantially revised.

Referee Report

3 major / 3 minor

Summary. The paper solves the three-dimensional Dirac equation for an electron in the field of an infinite chain of positively charged ions, retaining only the zero-order Fourier component of the chain potential, which is the logarithmic potential V0(r_perp) of Eq. (6). It then replaces this potential by an effective 1/r_perp potential with an arbitrary parameter C_eff, reduces the Dirac equation to two radial systems, and obtains explicit Laguerre-polynomial eigenfunctions (40)-(44) and energy eigenvalues (31)-(32). The authors claim these are the exact analytical bound states of the Dirac equation with the realistic chain potential (6), and attribute an accidental degeneracy to a new spinor invariant (45). The central exactness claim is not supported because the solution is constructed for the fitted 1/r_perp potential, not for V0, and C_eff is never determined from V0.

Significance. If the central claim were correct, the paper would provide a useful relativistic description of bound states on a charged chain with naturally emerging spin-orbit coupling, and the new invariant would be of interest. The algebraic reduction of the Dirac equation in cylindrical coordinates is carried out in detail, and the resulting eigenfunctions for the effective 1/r_perp model are explicit and normalized. However, the undetermined parameter C_eff enters the spectrum directly, so the results do not constitute predictions for the actual chain potential. The claimed invariant (45) is not proved and is not used in the derivation. Thus the significance of the paper as stated is not established; its concrete value is limited to a solvable model with a free effective charge.

major comments (3)

[§3.2, after Eq. (26); Conclusions] The central claim of exactness for the chain potential is unsupported. The actual potential V0(r_perp) in Eq. (6) is logarithmic, V0 = -2 e^2 Z_v/a ln( a/(2 r_perp) + sqrt(1 + a^2/(4 r_perp^2)) ), and is not proportional to 1/r_perp. The derivation replaces it with V_eff = C_eff e^2 Z_v/r_perp, introducing a free parameter C_eff that is never fixed from V0. All subsequent results, including the quantization condition (28), the spectrum (31)-(32), the damping parameter (33), and the wavefunctions (40)-(44), are solutions of that effective model. Because the energies depend explicitly on C_eff, the spectrum is a one-parameter family of fitted curves, not a prediction. The sentence in the Conclusions claiming 'the exact analytical solution of the DE with the realistic Coulomb potential (6)' therefore contradicts the actual derivation.
[§3.2, Eq. (45)] The new spinor invariant (45) is presented without a proof that it commutes with the Dirac Hamiltonian. It is not used to obtain the eigenfunctions, yet the abstract and the Conclusions credit the solution to this invariant, and the claimed accidental degeneracy rests on it. The paper should either provide a direct commutator calculation or state explicitly for which Hamiltonian (V0 or V_eff) the invariant is valid. As it stands, the existence of this invariant is an unverified assertion.
[§3.2, Eq. (6) and the passage introducing C_eff] The replacement of V0 by V_eff is not justified uniformly in r_perp. For small r_perp, V0 behaves as -2 e^2 Z_v/a ln(a/(2 r_perp)), whereas V_eff diverges as -C_eff e^2 Z_v/r_perp. These are qualitatively different; no single constant C_eff can reproduce the logarithmic behavior. Since bound-state wavefunctions and normalization are sensitive to the small-r region, the undetermined C_eff cannot absorb this discrepancy. This reinforces that the obtained states are not approximate or exact bound states of the chain potential (6), even in a limiting sense.

minor comments (3)

[Eqs. (40)-(44) and Eq. (43)] The symbol b appears in the polynomials P_n,± and Q_n,± and in the normalization constant (43) but is never defined. From the relation between constants in Eq. (35) and the definition of δ in Eq. (36), b presumably equals δ, but the paper should state this explicitly.
[Eq. (44) and surrounding text] The argument of the polynomials is written as 2κ⊥ρ, but ρ was already defined in Eq. (25) as the ratio sqrt((E_k-E)/(E_k+E)). The intended argument is the radial coordinate 2κ⊥r_perp; the notation should be changed to avoid ambiguity.
[Eq. (28) and Eq. (30)] The quantization condition is stated for integer n, but the allowed range of n is not discussed in relation to the half-integral values of M and the requirement that γ_M be real (M^2 > C_eff^2 Z_v^2 α^2). A short discussion of parameter ranges would improve clarity.

Circularity Check

1 steps flagged · score 6.0 of 10

The claimed exact solution is exact only for the fitted 1/r potential V_eff, not for the actual chain potential V0 of Eq. (6), because C_eff is never fixed from V0.

fitted input called prediction [Section 3.2, after Eq. (26), before Eq. (27); and Eqs. (31)-(32)]
"At r⊥ >a, the potential can be expanded in the series V0(r⊥ ) = −2e2Zv/a ln(...) ≈ −e2Zv/r⊥(1 − a2/24r2⊥ + ...), and, to obtain the solution for the electrons bound by the chain, we can use the eﬀective potential Vef f =Cef fe2Zv/r ⊥ with the parameter Cef f."

The actual chain potential V0 from Eq. (6) is logarithmic, not a 1/r potential, and the constant C_eff is introduced as an adjustable parameter without being fixed by V0. Every subsequent result — the quantization condition (28), the energy E = E_k Δ_{n,M} with Δ_{n,M} depending on C_eff in Eq. (32), and the eigenfunctions (40)-(44) — is a solution of the Dirac equation with V_eff, not with V0. Thus the paper's predicate of 'exact analytical expressions for bound electron states in the Coulomb field of the chain' reduces to a one-parameter fit: the spectrum is a function of the fitted input C_eff rather than a prediction forced by the chain potential.

full rationale

The central derivation is self-contained algebraically, but its physical claim is not. The paper replaces the logarithmic potential V0 of Eq. (6) by the effective 1/r potential V_eff = C_eff e^2 Z_v/r_perp and never determines C_eff from the actual potential; the final spectrum and states carry C_eff as a free parameter. Consequently the 'exact analytical solution with the realistic Coulomb potential (6)' asserted in the Conclusions is not a prediction about the charged chain: it is the exact solution for a different, fitted potential. The self-citations [10,11] and the new invariant (45) are not the load-bearing part of this problem; the load-bearing defect is that the adjustable C_eff enters the energy and normalization, so the claimed result is parameterized by its own input. This warrants a partial-circularity score of 6 rather than a higher score, because the algebra after the substitution is legitimate and the paper does not rely on a self-citation chain.

Assumptions & free parameters 1 free parameters · 4 assumptions · 0 invented entities

The central claim rests on a free parameter C_eff and on the unproven commutation of the invariant (45). No new physical entities are introduced.

free parameters (1)

C_eff
Effective potential strength parameter introduced to replace the true chain potential with a 1/r Coulomb form. The energy spectrum and wavefunctions depend on it, but it is not determined from first principles.

assumptions (4)

domain assumption Zero-order Fourier component V0 represents the electron-chain interaction.
The paper uses V0 (Eq. (6)) as the potential, neglecting higher harmonics without justifying their smallness for bound states near the chain.
ad hoc to paper Effective Coulomb potential V_eff = C_eff e^2 Z_v / r_perp approximates V0.
Introduced with undetermined constant C_eff; the solution depends on this choice and no prescription for C_eff is given.
ad hoc to paper The operator in Eq. (45) commutes with the Dirac Hamiltonian.
Stated without proof; used to explain degeneracy. The commutation is not demonstrated in the paper.
standard math Standard Dirac equation and matrix algebra.
Background used throughout, including the form of Dirac matrices and the separation ansatz.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Bound electron states in a charged chain within the Dirac description." pith.science (2026). https://pith.science/paper/RS6V4JLB

@misc{pith2026241114212,
  author       = {Pith},
  title        = {Pith review of: Bound electron states in a charged chain within the Dirac description},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/RS6V4JLB}},
  note         = {Machine review of arXiv:2411.14212}
}

read the original abstract

For the first time the exact analytical expressions for the three-dimensional bound electron states in the Coulomb field of the chain consisting of positively charged ions, are obtained within the Dirac description, using the new spinor invariant found for this problem. It is demonstrated that within such approach the coupling between electron spin and its one-dimensional propagation along the chain naturally arise, without any need to include artificially into the equations the so-called spin-orbit interaction.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

15 extracted references · 13 canonical work pages

[1]

Dauxois, M

T. Dauxois, M. Peyrard. Physics of Solitons. Cambridge, Cambridge University Press, 2006

work page 2006
[2]

Kosevich, B.A

A.M. Kosevich, B.A. Ivanov, and A.S. Kovalev. Magnetic s olitons. Phys. Rep. 194, No 3/4, 117-238 (1990)

work page 1990
[3]

Davydov, N.I

A.S. Davydov, N.I. Kislukha. Solitary excitons in one-d imensional molecular chains. Phys. Stat. Sol. B (1973) 5 465. https://doi.org/10.1002/pssb.2220590212

doi:10.1002/pssb.2220590212 1973
[4]

A.S. Davydov. Solitons in Molecular Systems. Dordrecht : Reidel; 1991

work page 1991
[5]

A.C. Scott. Davydov’s soliton. Phys. Rep. (1992) 217 1

work page 1992
[6]

Kosevich, A.S

A.M. Kosevich, A.S. Kovalev. Autolocalization of vibra tions in one-dimensional an- harmonic chain. (In Russian). Zh. Exp. Theor. Phys., (1974) 67 1795

work page 1974
[7]

A.S. Davydov. Biology and Quantum Mechanics. Pergamon P ress, 1982. - Science - 229 pages. ISBN-10:0080263925; ISBN-13: 008026392.2

work page 1982
[8]

Brizhik, A.A

L.S. Brizhik, A.A. Eremko. Ground state diagram of a 1D el ectron-phonon system. Synth. Met., 2000, v. 109, No 1-3, 117-121

work page 2000

Show all 15 references

[9]

Brizhik, J

L.S. Brizhik, J. Luo, B.M.A.G. Piette, W.J. Zakrzewski. Long-range electron trans- port mediated by alpha-helices. Phys. Rev. E (2019) 100 062205. DOI: 10.1103/Phys- RevE.100.062205

2019 doi
[10]

Eremko, L

A.A. Eremko, L. Brizhik, V.M. Loktev. Algebra of the spi nor invariants and the relativistic hydrogen atom. Annals of Physics, 2023 , 451, 169266. https://doi.org/10.1016/j.aop.2023.169266

2023
[11]

Eremko, L.S

A.A. Eremko, L.S. Brizhik,V.M. Loktev, Spin relevant i nvariants and the general solution of the Dirac equation for the Coulomb ﬁeld. Annals o f Physics, 439 (2022) 168786 (20 pp). https://doi.org/10.1016/j.aop.2022.168 786

2022 doi
[12]

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

1950
[13]

Bound-state solution s of the Dirac equation for the Kratzer potential with pseudoscalar-Coulomb term

Altuˇ g Arda, and Ramazan Sever, "Bound-state solution s of the Dirac equation for the Kratzer potential with pseudoscalar-Coulomb term", Eu r. Phys. J. Plus (2019) 134: 29 DOI 10.1140/epjp/i2019-12421-9

2019 doi
[14]

Bound state solutions of t he Dirac equation for the trigonometric and hyperbolic Scarf-Grosche potential s using the Nikiforov-Uvarov method

W. A. Yahya and K. J. Oyewumi, "Bound state solutions of t he Dirac equation for the trigonometric and hyperbolic Scarf-Grosche potential s using the Nikiforov-Uvarov method", J. Math. Phys. 54, 013508 (2013); https://doi.org/10.1063/1.4772478

2013 doi
[15]

Analytical solution of the Dirac equation for the hyperbolic potential by the extended Nikiforov-Uvarov method

H. Karayer, "Analytical solution of the Dirac equation for the hyperbolic potential by the extended Nikiforov-Uvarov method", Eur. Phys. J. Plu s (2019) 134: 452 DOI 10.1140/epjp/i2019-12828-2. 12

2019 doi

Pith tools

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

[1] [1]

Dauxois, M

T. Dauxois, M. Peyrard. Physics of Solitons. Cambridge, Cambridge University Press, 2006

work page 2006

[2] [2]

Kosevich, B.A

A.M. Kosevich, B.A. Ivanov, and A.S. Kovalev. Magnetic s olitons. Phys. Rep. 194, No 3/4, 117-238 (1990)

work page 1990

[3] [3]

Davydov, N.I

A.S. Davydov, N.I. Kislukha. Solitary excitons in one-d imensional molecular chains. Phys. Stat. Sol. B (1973) 5 465. https://doi.org/10.1002/pssb.2220590212

doi:10.1002/pssb.2220590212 1973

[4] [4]

A.S. Davydov. Solitons in Molecular Systems. Dordrecht : Reidel; 1991

work page 1991

[5] [5]

A.C. Scott. Davydov’s soliton. Phys. Rep. (1992) 217 1

work page 1992

[6] [6]

Kosevich, A.S

A.M. Kosevich, A.S. Kovalev. Autolocalization of vibra tions in one-dimensional an- harmonic chain. (In Russian). Zh. Exp. Theor. Phys., (1974) 67 1795

work page 1974

[7] [7]

A.S. Davydov. Biology and Quantum Mechanics. Pergamon P ress, 1982. - Science - 229 pages. ISBN-10:0080263925; ISBN-13: 008026392.2

work page 1982

[8] [8]

Brizhik, A.A

L.S. Brizhik, A.A. Eremko. Ground state diagram of a 1D el ectron-phonon system. Synth. Met., 2000, v. 109, No 1-3, 117-121

work page 2000