Home On solutions of the Dirac equation for 1D hydrogenic atoms or ions
Article Open Access

On solutions of the Dirac equation for 1D hydrogenic atoms or ions

  • Eugene Oks EMAIL logo
Published/Copyright: July 22, 2025

Abstract

Studies of 1D hydrogenic atoms or ions reveal interesting physics and could be relevant to the following experimental situations: (1) hydrogenic atoms in very strong magnetic fields where the electron is effectively confined in a finite-radius cylinder surrounding a field line; (2) an electron constrained to the surface of a nanotube in the field of an ion; (3) excitons in research of semiconductors, superconductivity, and polymers; (4) electron at the surface of liquid helium. In the present theoretical paper, we revisited the 1D Dirac–Coulomb problem. We employed a model where the Coulomb potential is truncated by a constant value at some relatively small but finite distance R from the origin, thus allowing for the finite size of the nucleus of charge Z. We derived analytically the wave functions in the exterior and interior regions. We showed that for nonzero energies, the wave functions in the exterior and interior regions and their derivatives cannot be matched at the boundary, so that the corresponding mathematical solutions are not physically acceptable. For zero energy, we determined the value of R(Z), at which the interior and exterior wave functions and their derivatives can be matched at the boundary. We demonstrated that for the electronic hydrogenic atoms or ions, the matching value of R(Z) is greater than the corresponding nuclear sizes by at least two orders of magnitude. However, for the muonic hydrogenic atoms or ions, the matching value of R(Z) turned out to be approximately equal to the nuclear size R 0(Z) (within the accuracy of the expressions for R 0(Z)) for the range of the nuclear charges 14 ≤ Z ≤ 25 (i.e., from Si to Mn), the best match of R and R 0 being for Z = 18 (Ar) and Z = 19 (K).

1 Introduction

There is a very large number of studies of 1D hydrogen atoms (many of the studies are referred in this work). From the theoretical viewpoint, they are fascinating for several reasons. First, they reveal interesting physics. Second, on some occasions, they contradict each other, causing disputes (the disputes being the driving force in search for the scientific truth). For example, one disputed issue is whether it should be considered an impenetrable barrier at the origin, while another disputed issue is even whether there are bound states (specific references are given here).

Third, the total number of the relevant theoretical publications is ∼102. The overwhelming majority of them are dealing with the Schrödinger equation: therefore, here we refer only to some works published in the last 30 years (out of about 70 years of such studies) – namely, studies by Golovinskii and Preobrazhenskii [1,2], Loudon [3], Palma and Raff [4], Nieto [5], Gordeyev and Chhajlany [6], Fischer et al. [7], Xianxi et al. [8], and references therein. Andrews [9,10] insisted that an impenetrable barrier should be considered at the origin, while some other authors disagreed. Of course, all the results of these studies were non-relativistic since they were obtained based on the Schrödinger equation.

There are also relevant relativistic theoretical studies. For such studies based on the Klein–Fock–Gordon equation, we refer to the studies by Barton [11], Spector and Lee [12], and references therein. For the relevant studies based on the Dirac equation, we refer to the studies by Audinet and Toulouse [13], Gusynin et al. [14], Ho and Khalilov [15], Benvegnu [16], Moss [17], Lapidus [18], and references therein.

As for the disputes in the studies based on the Dirac equation, Ho and Khalilov [15] set up an impenetrable barrier at the origin, while Audinet and Toulouse [13], Benvegnu [16], and Moss [17] did not. Last but not least, Moss [17] claimed that there are no bound states, while Ho and Khalilov [15] clearly showed the presence of the bound states.

From the experimental viewpoint, the above theoretical results could be useful for some of the following studies:

  1. Hydrogenic atoms in very strong magnetic fields where the electron is effectively confined in a finite-radius cylinder surrounding a field line (e.g., book by Ruder et al. [19]);

  2. An electron constrained to the surface of a nanotube in the field of an ion (e.g., study by Banyai et al. [20] and references therein);

  3. Excitons in research of semiconductors (Elliot and Loudon [21,22]), superconductivity (Ginzburg [23]), Wigner crystal (Carr [24], Wigner [25]), and polymers (Abe and Su [26], Abe [27], and Heeger et al. [28]);

  4. Electron at the surface of liquid helium – e.g., studies by Kawakami et al. [29], Golovinskii and Preobrazhenskii [2], Nieto [5], Cheng et al. [30], Bridges and McGill [31], Grimes et al. [32], and references therein.

We note that among the above 1D theoretical studies, some of them incorporated the effects of the finite nuclear size and some of them did not. In the 3D situation, there are a number of studies incorporating the effects of the finite nuclear size. Examples are the studies by Shabaev [33], Franzke and Gertsmann [34], Lee and Milstein [35], Salman and Saue [36], Pachucki et al. [37], Oks [38,39,40,41].

Since the present work focuses on the 1D situation, the above 3D theoretical studies of the effects of the finite nuclear size are only tangentially related to our work.

In the present theoretical study, we revisit the 1D Dirac–Coulomb problem. We start from the Dirac equation in the form utilized in the illuminating study by Ho and Khalilov [15] (hereafter, HK study). However, HK employed the model where the potential is the Coulomb one down to the origin with the impenetrable wall at the origin, the finite potential being defined only for the positive values of the coordinate. In distinction, we use a more realistic model where the Coulomb potential is truncated by a constant value at some relatively small but finite distance R from the origin (thus allowing for the nuclear size) and is defined both for the positive and negative values of the coordinate with no impenetrable wall at the origin. Besides, while HK considered only a 1D hydrogen atom, we consider a hydrogenic atom or ion of any nuclear charge Z. To the best of our knowledge, there were no studies of the solutions of the 1D Dirac equation for the potential that we use.

We obtain analytically the wave functions in the exterior and interior regions and then deduce the value of R(Z), at which the interior and exterior wave functions and their derivatives can be matched at the boundary. We show that for the electronic hydrogenic atoms or ions, the matching value of R(Z) exceeds the corresponding nuclear sizes by at least two orders of magnitude. However, for the muonic hydrogenic atoms or ions, the matching value of R(Z) is approximately equal to the nuclear size R 0(Z) (within the accuracy of the formula for R 0(Z)) for the range of the nuclear charges 14 ≤ Z ≤ 25 (i.e., from Si to Mn).

2 New results

2.1 General calculations

In the HK study [15], the solutions of the 1D Dirac equation have been sought on the positive half-line x ≥ 0 for a negative fermion in the field described by the following model potential energy:

(1) V HK ( x ) = q / x , ( x > 0 ) , V HK ( 0 ) = .

In distinction, we use the potential

(2) V ( x ) = q / x , ( x R ) , V ( x ) = q / R , ( R x 0 ) ,

where R = const > 0.

The fermion can be either electron or muon – so we use the natural units ħ = c = m = 1, where m is either the electron mass or the muon mass, respectively. For the quantity q, we take

(3) q = Z α ,

where Z is the nuclear charge and α is the fine structure constant. (Thus, we extend the study to hydrogen-like ions, while HK considered only hydrogen atoms.) Then, the Dirac equation for the two components ψ 1(x) and ψ 2(x) of the wave function can be written similar to Eqs. (5) and (6) from the HK study, but only for xR:

(4) d ψ 1 / d x + ( 1 q / x ) ψ 1 = E ψ 2 ,

(5) d ψ 2 / d x + ( 1 q / x ) ψ 2 = E ψ 1 ,

where E is the eigenenergy. Below we attempt to find physically acceptable analytical solutions of Eqs. (4) and (5) in the range of x ≥ 0. After finding such solution, we will properly extend it to the range of x < 0.

On expressing ψ 1 via ψ 2 from Eq. (5) and substituting in Eq. (4), we obtain

(6) x 2 d 2 ψ 2 / d x 2 + [ ( E 2 1 ) x 2 + 2 q x q ( q + 1 ) ] ψ 2 = 0 .

There is a reference book titled “Handbook of Exact Solutions for Ordinary Differential Equations” by Polyanin and Zaitsev [42], which is one of the most comprehensive handbooks on the subject. As one of the numerous examples, the authors of the handbook [42] consider the following equation (example 2.1.2.110):

(7) x 2 d 2 y / d x 2 + ( a x 2 + b x + c ) y = 0 .

They note that the first step to the exact solution is the substitution y(x) = x γ u(x), where γ is the solution of the equation γ 2γ + c = 0, which in our case is

(8) γ 2 γ q ( q + 1 ) = 0 .

The two solutions of Eq. (8) are

(9) γ + = 1 + q , γ = q .

In the case of Eq. (6), the substitution

(10) ψ 2 ( x ) = x γ u ( x ) .

yields

(11) x d 2 u / d x 2 + 2 γ d u / d x + [ ( E 2 1 ) x + 2 q ] u = 0 .

Here and below, we temporarily omitted the subscript of the parameter γ for brevity. The exact solutions of the equations of the type of Eq. (11) are given in the handbook [42] in example 2.1.2.103

(12) u ± ( x ) = const exp [ ± ( 1 E 2 ) 1 / 2 x ] J [ γ ± q / ( 1 E 2 ) 1 / 2 , 2 γ , ± 2 ( 1 E 2 ) 1 / 2 ( x ) ] ,

where J(A, B, w) is either the Kummer function Φ(A, B, w) or the Tricomi function U(A, B, w).

First, we explore the case where there is the Kummer function in Eq. (12) and the minus sign in the argument of the exponential. The intermediate result is

(13) ψ 2 ( x ) = const x γ exp [ ( 1 E 2 ) 1 / 2 x ] Φ [ γ q / ( 1 E 2 ) 1 / 2 , 2 γ , 2 ( 1 E 2 ) 1 / 2 x ] .

The final step in the fulfillment of the boundary condition ψ 2(∞) = 0 is to ensure that the Kummer function in Eq. (13) reduces to a polynomial of the finite degree. This is achieved by the requirement that

(14) γ q / ( 1 E 2 ) 1 / 2 = k , k = 0 , 1 , 2 ,

On solving Eq. (14) for the energy E, we obtain (for the positive set of E)

(15) E n = { 1 q 2 / ( γ + k ) 2 } 1 / 2 , k = 0 , 1 , 2 ,

For γ = γ + = 1 + q, Eq. (15) for E n coincides with the corresponding expression from the HK study (i.e., with their Eq. (8)), which in our notations (i.e., with m = 1) is

(16) E n , HK = [ 1 q 2 / ( q + n ) 2 ] 1 / 2 , n = k + 1 ,

where the subscript HK refers to the result from the HK study.

Now, we go back to the general form of ψ 2(x) from Eq. (12). It turned out that the choice of the Tricomi function coupled with the choice of the minus sign in the argument of the exponential, as well as the choice of the Kummer function coupled with the choice of the plus sign in the argument of the exponential, yielded the results for ψ 2(x) – for the energies E n from Eq. (15) – that are identical to the above case (presented in detail) where the choice was of the Kummer function coupled with the choice of the minus sign in the argument of the exponential. As for the last remaining choice of the Tricomi function coupled with the choice of the plus sign in the argument of the exponential, it yielded a (complex-valued) function that does not satisfy the condition ψ 2(∞) = 0 for any energy.

Next we will solve the Dirac equation in the range Rx ≥ 0, where the potential energy is V(x) = – q/R = const. In this range, the Dirac equation has the form

(17) d ψ 1 / d x + ( 1 q / R ) ψ 1 = E ψ 2 ,

(18) d ψ 2 / d x + ( 1 q / R ) ψ 2 = E ψ 1 .

At x = 0, the solutions of these equations should be limited, but do not have to vanish – because we did not require V(0) = ∞. Below we show that all solutions are limited at x = 0. So, the crucial test of whether the solutions are physically acceptable will be to try to match them with the corresponding exterior solutions and their logarithmic derivatives at x = R. As x approaches R from the outside, for the exterior functions we obtain (Eq. (13))

(19) ψ 2,ext R γ , ( 1 / ψ 2 , ext ) d ψ 2,ext / d x γ / R .

If E ≠ 0, then (Eq. (5))

(20) ψ 1,ext ( γ / E ) R γ 1 , ( 1 / ψ 1,ext ) d ψ 1,ext / d x ( γ 1 ) / R .

For the interior solutions, in the case of E ≠ 0, from Eq. (18), we obtain

(21) ψ 1 = ( 1 / E ) [ d ψ 2 / d x + ( 1 q / R ) ψ 2 ] ,

and on substituting Eq. (21) in Eq. (17), we obtain

(22) d 2 ψ 2 / d x 2 + [ E 2 ( 1 q / R ) 2 ] ψ 2 = 0 .

Below we explore three subcases of the case E ≠ 0.

2.2 Subcase of E 2 < (1 – q/R)2

If E 2 < (1 – q/R)2, then the solutions of the Eq. (20) are

(23) ψ 2,int ( x ) = exp { ± [ ( 1 q / R ) 2 E 2 ] 1 / 2 x } .

Then, for the logarithmic derivative at x = R, we obtain

(24) ( 1 / ψ 2,int ) d ψ 2,int / d x = ± [ ( 1 q / R ) 2 E 2 ] 1 / 2 .

By equating the logarithmic derivatives of ψ 2(x) on both sides of x = R, we obtain

(25) ± [ ( 1 q / R ) 2 E 2 ] 1 / 2 = γ ± / R .

On substituting ψ 2(x) from Eq. (23) in Eq. (22), we obtain

(26) ψ 1,int ( x ) = ( 1 / E ) { 1 q / R ± ( 1 ) [ ( 1 q / R ) 2 E 2 ] 1 / 2 } exp { ± [ ( 1 q / R ) 2 E 2 ] 1 / 2 x } .

Then, for the logarithmic derivative at x = R, we obtain

(27) ( 1 / ψ 1,int ) d ψ 1,int / d x = ± [ ( 1 q / R ) 2 E 2 ] 1 / 2 .

By equating the logarithmic derivatives of ψ 1(x) on both sides of x = R, we get

(28) ± [ ( 1 q / R ) 2 E 2 ] 1 / 2 = ( γ ± 1 ) / R .

Obviously, Eqs. (25) and (28) are incompatible with each other. Thus, in the subcase of E 2 < (1 – q/R)2], there are no physically acceptable solutions.

2.3 Subcase of E 2 > (1 – q/R)2

Now we proceed to the subcase where E 2 > (1 – q/R)2. We denote

(29) p = [ E 2 ( 1 q / R ) 2 ] 1 / 2 .

Then, from Eq. (22) follows that either

(30) ψ 2,int ( x ) = sin ( p x ) ,

or

(31) ψ 2,int ( x ) = cos ( p x ) .

For the sub-subcase where ψ 2,int(x) = sin(px), by equating the corresponding logarithmic derivatives at x = R, we obtain

(32) p / tan ( p R ) = γ / R ,

so that

(33) tan ( p R ) = p R / γ .

Here we temporarily omitted the subscript ± at γ for brevity. For ψ 1,int(x) in this sub-subcase we obtain (Eq. (21))

(34) ψ 1,int ( x ) = ( 1 / E ) [ ( 1 q / R ) sin ( p x ) p cos ( p x ) ] .

The logarithmic derivative calculated at x = R is

(35) ( 1 / ψ 1,int ) d ψ 1,int / d x = p [ p tan ( p R ) + 1 q / R ] / [ ( 1 q / R ) tan ( p R ) p ] .

By equating it to the corresponding logarithmic derivative in the exterior, we obtain (Eq. (20))

(36) p [ p tan ( p R ) + 1 q / R ] / [ ( 1 q / R ) tan ( p R ) p ] = ( γ 1 ) / R .

On substituting Eq. (33) in Eq. (36), after some elementary calculations, we obtain the following equation with respect to unknown R:

(37) R 2 + R / p 2 ( R ) + γ ( γ 1 ) / p 2 ( R ) = 0 .

In view of Eq. (9), for γ + = 1 + q, the quantity γ(γ – 1) = q(1 + q); for γ = – q, the result for the quantity γ(γ – 1) is the same: γ(γ – 1) = q(1 + q). So, Eq. (37) becomes

(38) R 2 + R / p 2 ( R ) + q ( 1 + q ) / p 2 ( R ) = 0 .

Since R > 0, all three terms in the left side of Eq. (38) are positive, so that Eq. (38) does not have solutions for R > 0. Consequently, the sub-subcase where ψ 2,int = sin(px), does not lead to a physically acceptable solution.

Now, we analyze the sub-subcase where ψ 2,int(x) = cos(px). By equating the corresponding logarithmic derivatives at x = R, we obtain

(39) tan ( p R ) = γ / ( p R ) .

For ψ 1,int(x) in this subcase, we obtain

(40) Ψ 1,int ( x ) = ( 1 / E ) [ p sin ( p x ) + ( 1 q / R ) cos ( p x ) ] .

The logarithmic derivative calculated at x = R is

(41) ( 1 / ψ 1,int ) d ψ 1,int / d x = p [ p ( 1 q / R ) tan ( p R ) ] / [ p tan ( p R ) + 1 q / R ] .

By equating it to the corresponding logarithmic derivative in the exterior, we obtain

(42) p [ p ( 1 q / R ) tan ( p R ) ] / [ p tan ( p R ) + 1 q / R ] = ( γ 1 ) / R .

On substituting Eq. (39) in Eq. (42), after denoting y = 1/R and performing some elementary calculations, we obtain the following equation with respect to unknown y:

(43) [ γ ( γ 1 ) q ] y 2 + y + p 2 = 0 .

We found above that both for γ + and for γ , we have γ(γ – 1) = q(1 + q), then Eq. (43) becomes

(44) q 2 y 2 + y + p 2 = 0 .

Since R > 0, so that y > 0, then all three terms in the left side of Eq. (44) are positive, so that Eq. (38) does not have solutions for R > 0. Consequently, the sub-subcase where ψ 2,int = cos(px), does not lead to a physically acceptable solution either.

2.4 Subcase of E 2 = (1 – q/R)2

Now, we proceed to the subcase where E 2 = (1 – q/R)2]. Then, the solutions of Eq. (22) are

(45) ψ 2,int ( x ) = x β , β = 0 or β = 1 .

If β = 0, so that ψ 2,int(x) = const and dψ 2,int(x)/dx = 0, then the logarithmic derivative (1/ψ 2,int)dψ 2,int/dx = 0: it does not match the corresponding logarithmic derivative of the exterior function ψ 2,ext(x) from Eq. (19), which is equal to γ/R.

If β = 1, so that ψ 2,int(x) = x and dψ 2,int(x)/dx = 1, so that the logarithmic derivative of the function ψ 2,int(x) calculated at x = R is 1/R, it still does not match the corresponding logarithmic derivative of the exterior function ψ 2,ext(x) from Eq. (19), which is equal to γ/R.

Alternatively, we can express ψ 2 via ψ 1 from Eq. (17), substitute in Eq. (18), and obtain (for E 2 = (1 – q/R)2]): d2 ψ 1,int/dx 2 = 0. Then,

(46) ψ 1,int ( x ) = x β , β = 0 or β = 1 .

If β = 0, then (1/ψ 1,int)dψ 1,int/dx = 0, so that it does not match the corresponding logarithmic derivative of the exterior function ψ 1,ext(x) from Eq. (19), which is equal to (γ – 1)/R.

If β = 1, so that ψ 1,int(x) = x and dψ 1,int(x)/dx = 1, so that the logarithmic derivative of the function ψ 1,int(x) calculated at x = R is 1/R, it still does not match the corresponding logarithmic derivative of the exterior function ψ 1,ext(x) from Eq. (20), which is equal to (γ – 1)/R.

Thus, there are no physically acceptable solutions for E ≠ 0.

2.5 Subcase of E = 0

Finally, we study the case of E = 0. In this case, the Dirac equation simplifies to

(47) d ψ 1 / d x + [ 1 + V ( x ) ] ψ 1 = 0 ,

(48) d ψ 2 / d x + [ 1 + V ( x ) ] ψ 2 = 0 .

In HK study, the case of E = 0 was studied in the simpler model where the potential energy was V HK(x) = –q/x, (x > 0), V HK(0) = ∞. Within this simpler model, HK obtained the following result for the components of the Dirac bispinor:

(49) ψ 2 ( x ) = 0 , ψ 1 ( x ) = x q exp ( x ) .

After we normalized that solution – what was not done by HK – we found that the normalization constant in front of x q in Eq. (29) should be 2 q /[qΓ(2q)]1/2, where Γ(2q) is the gamma function.

For our more complicated model, where the potential energy is given by Eq. (2), the solution given by Eq. (49) is valid only in the exterior, i.e., at xR. In the interior region (|x| ≤ R), where V(x) = –q/R, from Eqs. (47) and (48), we find

(50) ψ 2,int ( x ) = 0 , ψ 1,int ( x ) = C exp [ ( q / R 1 ) x ] ,

where C = const will be determined below.

We emphasize that in the interior region (|x| ≤ R), there is no antisymmetric solution, but only the symmetric one. Therefore, the solution must be symmetric in the entire range of –∞ < x < ∞. In particular, the exterior solution for ψ 1,ext(x) must be |x| q exp(–|x|).

Obviously, ψ 2(x), being identically zero both in the interior and exterior regions, satisfies the matching condition at x = R.

So, we focus on ψ 1(x). In the exterior region, where ψ 1(x) has the form given in Eq. (49) with the substitution of x by |x|, its logarithmic derivative, calculated at x = R, is

(51) ( 1 / ψ 1,ext ) d ψ 1,ext / d x = q / R 1 .

The logarithmic derivative of ψ 1,int(x) given in Eq. (50), being calculated at x = R, yields

(52) ( 1 / ψ 1,int ) d ψ 1,int / d x = q / R 1 .

Thus, the logarithmic derivatives of ψ 1,int(x) and ψ 1,ext(x) match at x = R regardless of the value of R.

Now, let us require the functions ψ 1,int(x) and ψ 1,ext(x) to match at x = R:

(53) C exp [ q R ] = R q exp ( R ) .

Thus, we find the value of the constant C

(54) C = R q exp ( q ) .

The normalization integral N is as follows:

(55) N = 2 0 R d x R 2 q exp ( 2 q 2 + 2 q / R ) + 2 R d x x 2 q exp ( 2 x ) .

On calculating the integral, we obtain

(56) N = R 2 q [ exp ( 2 R ) exp ( 2 q ) ] / ( q / R 1 ) + 2 2 q Γ ( 1 + 2 q , 2 R ) ,

where Γ(1 + 2q, 2R) is the incomplete gamma function.

After equating

(57) N = 1 ,

we obtain the following relation between the value of boundary R (for which the interior and exterior wave functions and their derivatives can be matched at x = R) and q

(58) R 2 q [ exp ( 2 R ) exp ( 2 q ) ) q / R 1 + 2 2 q Γ ( 1 + 2 q , 2 R ) = 1 .

Figure 1 presents the dependence of the boundary R (in natural units) on the nuclear charge Z (by using q = ).

Figure 1 
                  Dependence of the boundary R (in natural units) on the nuclear charge Z (by using q = Zα).
Figure 1

Dependence of the boundary R (in natural units) on the nuclear charge Z (by using q = ).

Figure 2 presents the dependence of R (in fm) on the nuclear charge Z (by using q = ) – specifically for the situation where the negative fermion is electron. We remind that for electrons, the natural unit of length is 386.1 fm.

Figure 2 
                  Dependence of R (in fm) on the nuclear charge Z (by using q = Zα) – specifically for the situation where the negative fermion is electron.
Figure 2

Dependence of R (in fm) on the nuclear charge Z (by using q = ) – specifically for the situation where the negative fermion is electron.

Figure 3 presents the dependence of R (in fm) on the nuclear charge Z (by using q = ) – specifically for the situation where the negative fermion is muon. We remind that for muons, the natural unit of length is 1.867 fm.

Figure 3 
                  Dependence of R (in fm) on the nuclear charge Z (by using q = Zα) – specifically for the situation where the negative fermion is muon.
Figure 3

Dependence of R (in fm) on the nuclear charge Z (by using q = ) – specifically for the situation where the negative fermion is muon.

From Figures 13 (that were created using the software Mathematica) and the corresponding calculations, we see the following.

  1. As Z varies from 1 to 137, i.e., by 2 orders of magnitude, the value of the boundary R decreases only by about 12%.

  2. For the situation where the negative fermion is electron, the value of the boundary R exceeds the corresponding nuclear radii by at least 2 orders of magnitude. So, this situation seems to have only academic interest.

  3. For the situation where the negative fermion is muon, the value of the boundary R can be of the same order of magnitude as the corresponding nuclear radii. Thus, this situation seems to be more reasonable.

Below we provide more details related to the above item C. The dependence of the nuclear size R 0 on the nuclear charge is

(59) R 0 = 1.2 ( λ Z ) 1 / 3 fm , 2 λ 2.5 .

For the muonic case, the boundary R is approximately equal to the nuclear size R 0 within the accuracy of 10% (the uncertainty of R 0 from Eq. (59) being about 10%) for

(60) 14 Z 25 ,

that is, from Si to Mn. The best match of R and R 0 occurs for Z = 18 (Ar) and Z = 19 (K).

3 Conclusion

We revisited the 1D Dirac–Coulomb problem. Instead of the HK model [15] where the potential was the Coulomb one down to the origin with the impenetrable wall at the origin and where the finite potential was defined only for the positive values of the coordinate, we employed a more realistic model where the Coulomb potential is truncated by a constant value at some relatively small but finite distance R from the origin (thus allowing for the nuclear size) and is defined both for the positive and negative values of the coordinate with no impenetrable wall at the origin. Besides, while HK considered only a 1D hydrogen atom, we considered a hydrogenic atom or ion of any nuclear charge Z.

For the case of the non-zero energy E, we derived analytically the wave functions in the exterior and interior regions for all possible subcases and showed that the exterior and interior wavefunctions cannot be matched (together with their derivatives) at the boundary. In other words, we demonstrated that there are no physically acceptable solutions for any non-zero energy.

For the case of zero energy, we also derived analytically the wave functions in the exterior and interior regions. We determined the value of R(Z), at which the interior and exterior wave functions can be matched at the boundary (i.e., ensuring the continuity of the wave function and of its derivative at the boundary). We demonstrated that for the electronic hydrogenic atoms or ions, the matching value of R(Z) is greater than the corresponding nuclear sizes by at least two orders of magnitude. However, for the muonic hydrogenic atoms or ions, the matching value of R(Z) turned out to be approximately equal to the nuclear size R 0(Z) (within the accuracy of the expressions for R 0(Z)) for the range of the nuclear charges 14 ≤ Z ≤ 25 (i.e., from Si to Mn), the best match of R and R 0 being for Z = 18 (Ar) and Z = 19 (K).

We note that the presence of the zero-energy state causes the vacuum fermion number to be nonzero and to be fractional, as mentioned by HK. This happens also when the conjugation-symmetric Dirac Hamiltonian possesses an energy spectrum that is symmetric about zero, and a soliton is present, so that there is also an unpaired zero-energy state, which is self-conjugate (Jackiw [43]).

We also note that Barton [11] found a zero-energy solution while considering a truncated Coulomb potential in frames of the Klein–Gordon equation for a spinless charged particle. So, it seems that the relativistic effect of the presence of the zero-energy state is unrelated to the spin.

  1. Funding information: The author states no funding involved.

  2. Author contribution: The author has accepted responsibility for the entire content of this manuscript and approved its submission.

  3. Conflict of interest: The author states no conflict of interest.

  4. Data availability statement: All data generated or analyzed during this study are included in this published article.

References

[1] Golovinskii PA, Preobrazhenskii MA. Quasi-classical description of the Stark effect for an electron in an image-potential state. Opt Spectrosc. 2017;122:120.10.1134/S0030400X1701009XSearch in Google Scholar

[2] Golovinskii PA, Preobrazhenskii MA. Stimulated radiative transitions of an electron in an image potential. Opt Spectrosc. 2015;122:191.10.1134/S0030400X1502006XSearch in Google Scholar

[3] Loudon R One-dimensional hydrogen atom. Proc R Soc. 2016;472: 20150534.10.1098/rspa.2015.0534Search in Google Scholar

[4] Palma G, Raff U. The one-dimensional hydrogen atom revisited. Canad J Phys. 2006;84:787.10.1139/p06-072Search in Google Scholar

[5] Nieto MM. Electrons above a helium surface and one-dimensional Rydberg atom. Phys Rev A. 2005;61:034901.10.1103/PhysRevA.61.034901Search in Google Scholar

[6] Gordeyev AN, Chhajlany SC. One-dimensional hydrogen atom: a singular potential in quantum mechanics. J Phys A: Math Gen. 1997;30:6893.10.1088/0305-4470/30/19/025Search in Google Scholar

[7] Fischer W, Leschke H, Muller P. Comment on ‘On the Coulomb potential in one dimension’ by P Kurasov. J Phys A: Math Gen. 1997;30:5579.10.1088/0305-4470/30/15/037Search in Google Scholar

[8] Xianxi D, Dai J, Dai J. Orthogonality criteria for singular states and the nonexistence of stationary states with even parity for the one-dimensional hydrogen atom. Phys Rev A. 1997;55:2617.10.1103/PhysRevA.55.2617Search in Google Scholar

[9] Andrews M. The one-dimensional hydrogen atom. Am J Phys. 1988;56:776.10.1119/1.15476Search in Google Scholar

[10] Andrews M. Singular potentials in one dimension. Am J Phys. 1976;44:1064.10.1119/1.10585Search in Google Scholar

[11] Barton GJ. On the 1D Coulomb Klein–Gordon equation. Phys A: Math Theor. 2007;40:1011.10.1088/1751-8113/40/5/010Search in Google Scholar

[12] Spector HN, Lee J. Relativistic one-dimensional hydrogen atom. Am J Phys. 1985;53:248.10.1119/1.14132Search in Google Scholar

[13] Audinet T, Toulouse J. Effective quantum electrodynamics: One-dimensional model of the relativistic hydrogen-like atom. J Chem Phys. 2023;158:244108.10.1063/5.0152956Search in Google Scholar PubMed

[14] Gusynin VP, Sobol OO, Zolotaryuk AV, Zolotaryuk Y. Bound states of a one-dimensional Dirac equation with multiple delta-potentials. Low Temp Phys. 2022;48:1022.10.1063/10.0015111Search in Google Scholar

[15] Ho C-L, Khalilov VR. Fractional fermion number in a (1 + 1)-dimensional Dirac equation with a scalar Coulomb field. Phys Rev D. 2000;63:027701.10.1103/PhysRevD.63.027701Search in Google Scholar

[16] Benvegnu S. Relativistic point interaction with Coulomb potential in one dimension. J Math Phys. 1997;38:556.10.1063/1.531877Search in Google Scholar

[17] Moss RE. The hydrogen atom in one dimension. Am J Phys. 1987;55:397.10.1119/1.15144Search in Google Scholar

[18] Lapidus IR. Relativistic one‐dimensional hydrogen atom. Am J Phys. 1983;51:1036.10.1119/1.13445Search in Google Scholar

[19] Ruder H, Wunner G, Herold H, Geyer F. Atoms in strong magnetic fields. Berlin: Springer; 1994.10.1007/978-3-642-78820-8Search in Google Scholar

[20] Banyai L, Galbraith I, Ell C, Haug H. Excitons and biexcitons in semiconductor quantum wires. Phys Rev B. 1987;36:6099.10.1103/PhysRevB.36.6099Search in Google Scholar

[21] Elliot RJ, Loudon R. Theory of the absorption edge in semiconductors in a high magnetic field. J Phys Chem Solids. 1960;15:196.10.1016/0022-3697(60)90243-2Search in Google Scholar

[22] Elliot RL, Loudon R. Theory of fine structure on the absorption edge in semiconductors. J Phys Chem Solids. 1959;8:382.10.1016/0022-3697(59)90371-3Search in Google Scholar

[23] Ginzburg VL. Once again about high-temperature superconductivity. Contemp Phys. 1992;33:15.10.1080/00107519208219137Search in Google Scholar

[24] Carr Jr WJ. Energy, specific heat, and magnetic properties of the low-density electron gas. Phys Rev. 1961;122:1437.10.1103/PhysRev.122.1437Search in Google Scholar

[25] Wigner EP. Effects of the electron interaction on the energy levels of electrons in metals. Trans Faradey Soc. 1938;34:678.10.1039/tf9383400678Search in Google Scholar

[26] Abe S, Su WP. Excitons and charge transfer states in one-dimensional semiconductors. Mol Cryst Liq Cryst. 1991;194:357.10.1080/00268949108041188Search in Google Scholar

[27] Abe S. Exciton versus interband absorption in Peierls insulators. J Phys Soc Jap. 1989;58:62.10.1143/JPSJ.58.62Search in Google Scholar

[28] Heeger AJ, Kivelson S, Schrieffer JR, Su WP. Solitons in conducting polymers. Rev Mod Phys. 1988;60:781.10.1103/RevModPhys.60.781Search in Google Scholar

[29] Kawakami E, Elarabi A, Konstantinov D. Relaxation of the excited Rydberg states of surface electrons on liquid helium. Phys Rev Lett. 2021;126:106802.10.1103/PhysRevLett.126.106802Search in Google Scholar PubMed

[30] Cheng E, Cole MW, Cohen MH. Binding of electrons to the surface of liquid helium. Phys Rev B. 1994;50:1136.10.1103/PhysRevB.50.1136Search in Google Scholar PubMed

[31] Bridges F, McGill JF. Mobility of electrons on the surface of liquid helium. Phys Rev B. 1977;15:1324.10.1103/PhysRevB.15.1324Search in Google Scholar

[32] Grimes CC, Brown TR, Burns ML, Zipfel CL. Spectroscopy of electrons in image-potential-induced surface states outside liquid helium. Phys Rev B. 1976;13:140.10.1103/PhysRevB.13.140Search in Google Scholar

[33] Shabaev VM. Finite nuclear size corrections to the energy levels of the multicharged ions. J Phys B: Mol Opt Phys. 1993;26:1103.10.1088/0953-4075/26/6/011Search in Google Scholar

[34] Franzke KL, Gertsmann U. Finite nuclear size effect on the relativistic hyperfine splittings of 2s and 2p excited states of hydrogen-like atoms. Foundations. 2024;4:513.10.3390/foundations4040034Search in Google Scholar

[35] Lee RN, Milstein AI. Finite nuclear size and vacuum polarization in heavy atoms. Phys Lett A. 1994;189:72.10.1016/0375-9601(94)90820-6Search in Google Scholar

[36] Salman M, Saue T. Calculating the many-potential vacuum polarization density of the Dirac equation in the finite-basis approximation. Phys Rev A. 2023;108:012808.10.1103/PhysRevA.108.012808Search in Google Scholar

[37] Pachucki K, Patcos V, Yerokhin VA. Recoil corrections with finite nuclear size in hydrogenic systems. Phys Rev A. 2025;111:032820.10.1103/PhysRevA.111.032820Search in Google Scholar

[38] Oks E. High-energy tail of the linear momentum distribution in the ground state of hydrogen atoms or hydrogen-like ions. J Phys B: Mol Opt Phys. 2001;34:2235.10.1088/0953-4075/34/11/315Search in Google Scholar

[39] Oks E. Alternative kind of hydrogen atoms as a possible explanation for the latest puzzling observation of the 21 cm radio line from the early Universe. Res Astron Astrophys. 2020;20:109.10.1088/1674-4527/20/7/109Search in Google Scholar

[40] Oks E. Symmetry of the non-analytic solution of the Dirac equation inside the proton of hydrogen atoms. Symmetry. 2025;17:517.10.3390/sym17040517Search in Google Scholar

[41] Oks E. Resolution of the neutron lifetime puzzle and the conceptual design of its experimental confirmation. Nucl Phys B. 2025;1014:116879.10.1016/j.nuclphysb.2025.116879Search in Google Scholar

[42] Polyanin AD, Zaitsev VF. Handbook of exact solutions for ordinary differential equations. Boca Raton, FL: CRS Press; 1995.Search in Google Scholar

[43] Jackiw R. Fractional charge and zero modes for planar systems in a magnetic field. Phys Rev D. 1984;29:2375.10.1103/PhysRevD.29.2375Search in Google Scholar

Received: 2025-04-25
Revised: 2025-06-11
Accepted: 2025-06-12
Published Online: 2025-07-22

© 2025 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Research Articles
  2. Single-step fabrication of Ag2S/poly-2-mercaptoaniline nanoribbon photocathodes for green hydrogen generation from artificial and natural red-sea water
  3. Abundant new interaction solutions and nonlinear dynamics for the (3+1)-dimensional Hirota–Satsuma–Ito-like equation
  4. A novel gold and SiO2 material based planar 5-element high HPBW end-fire antenna array for 300 GHz applications
  5. Explicit exact solutions and bifurcation analysis for the mZK equation with truncated M-fractional derivatives utilizing two reliable methods
  6. Optical and laser damage resistance: Role of periodic cylindrical surfaces
  7. Numerical study of flow and heat transfer in the air-side metal foam partially filled channels of panel-type radiator under forced convection
  8. Water-based hybrid nanofluid flow containing CNT nanoparticles over an extending surface with velocity slips, thermal convective, and zero-mass flux conditions
  9. Dynamical wave structures for some diffusion--reaction equations with quadratic and quartic nonlinearities
  10. Solving an isotropic grey matter tumour model via a heat transfer equation
  11. Study on the penetration protection of a fiber-reinforced composite structure with CNTs/GFP clip STF/3DKevlar
  12. Influence of Hall current and acoustic pressure on nanostructured DPL thermoelastic plates under ramp heating in a double-temperature model
  13. Applications of the Belousov–Zhabotinsky reaction–diffusion system: Analytical and numerical approaches
  14. AC electroosmotic flow of Maxwell fluid in a pH-regulated parallel-plate silica nanochannel
  15. Interpreting optical effects with relativistic transformations adopting one-way synchronization to conserve simultaneity and space–time continuity
  16. Modeling and analysis of quantum communication channel in airborne platforms with boundary layer effects
  17. Theoretical and numerical investigation of a memristor system with a piecewise memductance under fractal–fractional derivatives
  18. Tuning the structure and electro-optical properties of α-Cr2O3 films by heat treatment/La doping for optoelectronic applications
  19. High-speed multi-spectral explosion temperature measurement using golden-section accelerated Pearson correlation algorithm
  20. Dynamic behavior and modulation instability of the generalized coupled fractional nonlinear Helmholtz equation with cubic–quintic term
  21. Study on the duration of laser-induced air plasma flash near thin film surface
  22. Exploring the dynamics of fractional-order nonlinear dispersive wave system through homotopy technique
  23. The mechanism of carbon monoxide fluorescence inside a femtosecond laser-induced plasma
  24. Numerical solution of a nonconstant coefficient advection diffusion equation in an irregular domain and analyses of numerical dispersion and dissipation
  25. Numerical examination of the chemically reactive MHD flow of hybrid nanofluids over a two-dimensional stretching surface with the Cattaneo–Christov model and slip conditions
  26. Impacts of sinusoidal heat flux and embraced heated rectangular cavity on natural convection within a square enclosure partially filled with porous medium and Casson-hybrid nanofluid
  27. Stability analysis of unsteady ternary nanofluid flow past a stretching/shrinking wedge
  28. Solitonic wave solutions of a Hamiltonian nonlinear atom chain model through the Hirota bilinear transformation method
  29. Bilinear form and soltion solutions for (3+1)-dimensional negative-order KdV-CBS equation
  30. Solitary chirp pulses and soliton control for variable coefficients cubic–quintic nonlinear Schrödinger equation in nonuniform management system
  31. Influence of decaying heat source and temperature-dependent thermal conductivity on photo-hydro-elasto semiconductor media
  32. Dissipative disorder optimization in the radiative thin film flow of partially ionized non-Newtonian hybrid nanofluid with second-order slip condition
  33. Bifurcation, chaotic behavior, and traveling wave solutions for the fractional (4+1)-dimensional Davey–Stewartson–Kadomtsev–Petviashvili model
  34. New investigation on soliton solutions of two nonlinear PDEs in mathematical physics with a dynamical property: Bifurcation analysis
  35. Mathematical analysis of nanoparticle type and volume fraction on heat transfer efficiency of nanofluids
  36. Creation of single-wing Lorenz-like attractors via a ten-ninths-degree term
  37. Optical soliton solutions, bifurcation analysis, chaotic behaviors of nonlinear Schrödinger equation and modulation instability in optical fiber
  38. Chaotic dynamics and some solutions for the (n + 1)-dimensional modified Zakharov–Kuznetsov equation in plasma physics
  39. Fractal formation and chaotic soliton phenomena in nonlinear conformable Heisenberg ferromagnetic spin chain equation
  40. Single-step fabrication of Mn(iv) oxide-Mn(ii) sulfide/poly-2-mercaptoaniline porous network nanocomposite for pseudo-supercapacitors and charge storage
  41. Novel constructed dynamical analytical solutions and conserved quantities of the new (2+1)-dimensional KdV model describing acoustic wave propagation
  42. Tavis–Cummings model in the presence of a deformed field and time-dependent coupling
  43. Spinning dynamics of stress-dependent viscosity of generalized Cross-nonlinear materials affected by gravitationally swirling disk
  44. Design and prediction of high optical density photovoltaic polymers using machine learning-DFT studies
  45. Robust control and preservation of quantum steering, nonlocality, and coherence in open atomic systems
  46. Coating thickness and process efficiency of reverse roll coating using a magnetized hybrid nanomaterial flow
  47. Dynamic analysis, circuit realization, and its synchronization of a new chaotic hyperjerk system
  48. Decoherence of steerability and coherence dynamics induced by nonlinear qubit–cavity interactions
  49. Finite element analysis of turbulent thermal enhancement in grooved channels with flat- and plus-shaped fins
  50. Modulational instability and associated ion-acoustic modulated envelope solitons in a quantum plasma having ion beams
  51. Statistical inference of constant-stress partially accelerated life tests under type II generalized hybrid censored data from Burr III distribution
  52. On solutions of the Dirac equation for 1D hydrogenic atoms or ions
  53. Entropy optimization for chemically reactive magnetized unsteady thin film hybrid nanofluid flow on inclined surface subject to nonlinear mixed convection and variable temperature
  54. Stability analysis, circuit simulation, and color image encryption of a novel four-dimensional hyperchaotic model with hidden and self-excited attractors
  55. A high-accuracy exponential time integration scheme for the Darcy–Forchheimer Williamson fluid flow with temperature-dependent conductivity
  56. Novel analysis of fractional regularized long-wave equation in plasma dynamics
  57. Development of a photoelectrode based on a bismuth(iii) oxyiodide/intercalated iodide-poly(1H-pyrrole) rough spherical nanocomposite for green hydrogen generation
  58. Investigation of solar radiation effects on the energy performance of the (Al2O3–CuO–Cu)/H2O ternary nanofluidic system through a convectively heated cylinder
  59. Quantum resources for a system of two atoms interacting with a deformed field in the presence of intensity-dependent coupling
  60. Studying bifurcations and chaotic dynamics in the generalized hyperelastic-rod wave equation through Hamiltonian mechanics
  61. A new numerical technique for the solution of time-fractional nonlinear Klein–Gordon equation involving Atangana–Baleanu derivative using cubic B-spline functions
  62. Interaction solutions of high-order breathers and lumps for a (3+1)-dimensional conformable fractional potential-YTSF-like model
  63. Hydraulic fracturing radioactive source tracing technology based on hydraulic fracturing tracing mechanics model
  64. Numerical solution and stability analysis of non-Newtonian hybrid nanofluid flow subject to exponential heat source/sink over a Riga sheet
  65. Numerical investigation of mixed convection and viscous dissipation in couple stress nanofluid flow: A merged Adomian decomposition method and Mohand transform
  66. Effectual quintic B-spline functions for solving the time fractional coupled Boussinesq–Burgers equation arising in shallow water waves
  67. Analysis of MHD hybrid nanofluid flow over cone and wedge with exponential and thermal heat source and activation energy
  68. Solitons and travelling waves structure for M-fractional Kairat-II equation using three explicit methods
  69. Impact of nanoparticle shapes on the heat transfer properties of Cu and CuO nanofluids flowing over a stretching surface with slip effects: A computational study
  70. Computational simulation of heat transfer and nanofluid flow for two-sided lid-driven square cavity under the influence of magnetic field
  71. Irreversibility analysis of a bioconvective two-phase nanofluid in a Maxwell (non-Newtonian) flow induced by a rotating disk with thermal radiation
  72. Hydrodynamic and sensitivity analysis of a polymeric calendering process for non-Newtonian fluids with temperature-dependent viscosity
  73. Exploring the peakon solitons molecules and solitary wave structure to the nonlinear damped Kortewege–de Vries equation through efficient technique
  74. Modeling and heat transfer analysis of magnetized hybrid micropolar blood-based nanofluid flow in Darcy–Forchheimer porous stenosis narrow arteries
  75. Activation energy and cross-diffusion effects on 3D rotating nanofluid flow in a Darcy–Forchheimer porous medium with radiation and convective heating
  76. Insights into chemical reactions occurring in generalized nanomaterials due to spinning surface with melting constraints
  77. Influence of a magnetic field on double-porosity photo-thermoelastic materials under Lord–Shulman theory
  78. Soliton-like solutions for a nonlinear doubly dispersive equation in an elastic Murnaghan's rod via Hirota's bilinear method
  79. 10.1515/phys-2025-0146
  80. 10.1515/phys-2025-0143
  81. Review Article
  82. Examination of the gamma radiation shielding properties of different clay and sand materials in the Adrar region
  83. Special Issue on Fundamental Physics from Atoms to Cosmos - Part II
  84. Possible explanation for the neutron lifetime puzzle
  85. Special Issue on Nanomaterial utilization and structural optimization - Part III
  86. Numerical investigation on fluid-thermal-electric performance of a thermoelectric-integrated helically coiled tube heat exchanger for coal mine air cooling
  87. Special Issue on Nonlinear Dynamics and Chaos in Physical Systems
  88. Analysis of the fractional relativistic isothermal gas sphere with application to neutron stars
  89. Abundant wave symmetries in the (3+1)-dimensional Chafee–Infante equation through the Hirota bilinear transformation technique
  90. Successive midpoint method for fractional differential equations with nonlocal kernels: Error analysis, stability, and applications
Downloaded on 17.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/phys-2025-0185/html
Scroll to top button