Resolution in the Case of the Hydrogen Atom of an Improved Dirac Equation ()
1. Introduction
The improved Dirac equation was obtained from Lochak’s theory of a leptonic magnetic monopole [1] [2] [3]. In this theory, the quantum wave has two U(1) gauge invariances. The mass term of the leptonic magnetic monopole is then able to replace the usual mass term of the Dirac equation. First, the non-linear mass term was read in the frame of Hestenes’s space-time algebra [4] [5] [6] - [11]. The improved equation may be obtained from the simplification of the Lagrangian density of the Dirac equation, where the
mass term is replaced by
. The Dirac wave of the linear equation or of the improved equation may be read as functions of space-time with values in the Clifford algebra
[12] - [34].
Since 1928, the relativistic invariance of Dirac’s theory used the previous Pauli matrices for the spin of the electron: the space-time variable
was replaced by
(1)
This is equivalent to saying that the three Pauli matrices:
(2)
form a orthogonal oriented basis in space. We shall put arrows on vectors in space, so any vector reads
(3)
The geometric algebra of space
and
are isomorphic algebras on the real field; the sum and the product of matrices are familiar in quantum physics. This matrix representation identifies complex numbers and scalar matrices in the Pauli algebra. With this identification, we write the x of (1) as
, we consider
as a basis in space-time and we use the Einstein summation convention, with superscripts and subscripts, where Latin indices are in space and Greek indices in space-time. Any element z in the Clifford algebra of space
is a sum of a real part x, a vector part
, an axial-vector part
and a pseudo-scalar part iy, and we need:
(4)
The application
is the main automorphism of
. The reverse is also the adjoint (transposed conjugate matrix), so
is the reversion. The third conjugation,
is the product of the two previous ones and we shall need:
(5)
Space-time is then made of the auto-adjoint part of the space algebra. We use:
(6)
The main reason to use the geometric algebra
is the ability to read all relativistic quantum physics in this algebra: the electron wave is a function of space and time in
:
(7)
The link between
and the complex formalism is simple only if we use the left and right Weyl spinors
and
, by letting:
(8)
The usual formalism uses a
and
matrices defined as following:
(9)
(10)
Our improved wave equation of the electron, which has the Dirac equation as linear approximation, reads [6] - [19]:
(11)
where
,
,
,
,
. This equation is invariant under any transformation D defined by an element M of the Lie group
(group of invertible elements of
):
(12)
(13)
(14)
Multiplying (11) on the left side by
, we obtain the usual form in the Pauli algebra of the improved equation:
(15)
where
is the Yvon-Takabayasi angle. In the usual formalism of complex
matrices, much more complicated than the Pauli algebra, this wave equation reads:
(16)
Our wave equation has the Dirac equation as linear approximation: if
is null or negligible, the improved Equation (15) is reduced to the Dirac equation that reads in the Pauli algebra and in the Dirac algebra:
(17)
(18)
Our equation is an improvement from the linear Dirac equation for many reasons explained in our previous papers [6] - [34]: Charge conjugation gives a positron with positive mass-energy. The gauge group comprises a local electric gauge linked to the conservative current of probability and a second gauge (only global) that is the chiral gauge of the Yvon-Takabayasi angle, linked to a second conservative current. This
gauge invariance corresponds to the first part of the
group of electro-weak interactions. Therefore, our improved equation allows us a generalization, with mass term, to all fermion waves of leptons and quarks. In this paper, we shall show how this improved, non-linear wave equation is nevertheless able to obtain the following results: the true quantum numbers, the true number of energy levels and the true energy levels in the case of the hydrogen atom. This was previously never obtained from a non-linear wave equation; therefore, everyone was thinking that the linearity of quantum mechanics was necessary, even for the linking to the non-linear relativistic gravitation.
2. Separating Variables
To solve the Dirac equation or the improved equation in the case of the hydrogen atom, two methods exist. We shall use here, not the initial method based on the non-relativistic approximation of the wave equation, but the new method invented by H. Krüger [35], separating the variables in spherical coordinates:
(19)
We use the following notation:
(20)
(21)
(22)
H. Krüger obtained the remarkable identity:
(23)
which with
(24)
gives also:
(25)
Aiming for the separation of the temporal variable
and the angular variable
from the radial variable r and the other angular variable
, we let:
(26)
where X is a function (with value in the Pauli algebra) of only r and
,
is the energy of the electron,
is an arbitrary phase (that plays no role here because the wave equations are electric gauge-invariant), and
is a real constant which will be interpreted as the magnetic quantum number (we name here this quantum number
because m is used in the mass term). We then get:
(27)
(28)
We also have:
(29)
(30)
(31)
Then if we let
(32)
we get:
(33)
Then with the form (26) of the wave, the Yvon-Takabayasi angle
depends neither on the time nor on the
angle. It depends only on r and
. Hence the separation of variables can be similarly obtained for the Dirac equation or for the improved equation. We have:
(34)
(35)
(36)
(37)
(38)
We then get:
(39)
For the hydrogen atom we have:
(40)
where
is the fine structure constant. We have:
(41)
Also the wave Equation (15) becomes:
(42)
which means:
(43)
The Dirac equation, in contrast, gives:
(44)
Now we let,
being the complex conjugate of z:
(45)
where
are functions with complex values of the real variables r and
. We get:
(46)
We then obtain the following equations:
(47)
(48)
(49)
(50)
(51)
Therefore, the improved equation is equivalent to:
(52)
Conjugating the equations containing the terms of the right columns we get the system:
(53)
Moreover we have:
(54)
then we get:
(55)
In the Equation (53) only two angular operators are present, then we let:
(56)
where A, B, C et D are functions of r while U and V are functions of
. The (53) system becomes:
(57)
Then if a
constant exists such that:
(58)
the (57) system becomes:
(59)
To get the system equivalent to the Dirac equation it is enough to suppress the
angle. This does not change the angular system (58), while in the place of (59) we get the system:
(60)
3. Kinetic Momentum Operators
We established in [8] the form that the operators of kinetic momentum have in space-time. With the Pauli algebra we have (a detailed calculation is in [13] A.3):
(61)
(62)
(63)
We indeed have also:
(64)
From (26) we have the following equivalence for
:
(65)
Then the
wave satisfying (26) is a proper vector of
and
is the magnetic quantum number. Moreover for a wave
satisfying (26), we have:
(66)
if and only if:
(67)
And (58) implies at the second order:
(68)
(69)
(70)
therefore
is a proper vector of
, with the proper value
, if and only if:
(71)
with the definition of S in (21) and with (26) we can see that the change of
into
conserves the same value for the wave if and only if
has a half-integer value. It is only in this case that the wave is correctly defined. The general results on the angular momentum operators imply then:
(72)
To solve the angular system we let, if
and with
:
(73)
while if
we let:
(74)
The angular system (58) is then equivalent [6] to the differential equation:
(75)
The change of variable:
(76)
gives then the differential equation of the Gegenbauer polynomials1:
(77)
And we get, as only integrable function:
(78)
with:
(79)
The
term is a factor of U and V, its argument may be absorbed by the
of (26), and its modulus may be transferred to the radial functions. We can then let
, which gives:
(80)
Since we have the (71) conditions on
and
, an integer n always exists such as:
(81)
This constrains the (80) series to be a finite sum, so then U and V are integrable. And since U and V have real values, we have:
(82)
4. Resolution of the Linear Radial System
We employ the following transformations:
(83)
The (60) system becomes:
(84)
And the (59) system becomes:
(85)
(86)
(87)
(88)
Since this radial system has the same asymptotic behavior as its linear approximation (84) we are left with the following as the only integrable solution:
(89)
We now study the case where
. The (88) equation is equivalent to:
(90)
Similarly, Equation (87) is equivalent to:
(91)
Next, Equation (86) is equivalent to:
(92)
Finally, Equation (85) is equivalent to:
(93)
We then arrive at three kinds of systems: index 0, index between 0 and n, and index n. For the null index the system is independent of
:
(94)
This system is the same as in the linear case. We obtain a non-null solution only if the determinant is null, then only if s satisfies:
(95)
And the system (94) is then equivalent to:
(96)
For the n index the resulting system is:
(97)
From the definition of
in (89) we obtain the cancellation of the determinants in the two previous sub-systems, and the (97) system is equivalent to:
(98)
Next the system with index between 0 and n reads:
(99)
(100)
(101)
(102)
Now we multiply (99) by
:
(103)
Next (100) implies:
(104)
And (103) becomes:
(105)
Using (97) we get:
(106)
Substituting these relations into (105) we obtain:
(107)
Since
we get:
(108)
(109)
(110)
Squaring the equation and substituting the definition of
we get:
(111)
(112)
(113)
which gives the Sommerfeld energy levels:
(114)
5. Constant Radial Polynomials
To arrive at all the results of the Dirac equation, there is one last thing to explain: why we have
different bound states with principal quantum number
. For this we must return to the particular case where the radial polynomial is reduced to a constant. We begin directly from (84), and we let:
(115)
We obtain the following set of equations from (59):
(116)
This is equivalent to the set formed by the four following systems:
(117)
(118)
(119)
(120)
The cancellation of the determinant in (117) and (118) gives again (95). The cancellation of the determinant in (117) and (118) is simply equivalent to
, which results from the definition of
. Each system (117) to (120) is then reduced to only one equation:
(121)
We then obtain:
(122)
We have a non-null solution only if:
(123)
Since
and
are positive we finally get:
(124)
This equality gives again the formula of the energy levels (114) with
. Since
comes with an absolute value, we can so well have
as
. But Darwin’s solutions, which work with real constants and not with complex ones at this stage of the calculation, forbid
to be negative and it is this interdiction which induces, for a principal quantum number
, a result of
states. This is what really happens: The change of sign of
is equivalent to changing V into −V. And if we change the sign of
and of V, then
are invariant if
, and the wave remains unchanged. Changing the sign of
brings no more supplementary solution and we can then use only the solutions with
, in the case
. And this allows us to get the true number of states.
6. Concluding Remarks
Sommerfeld’s formula (114) for energy levels does not account for the Lamb effect which is, if
, a very small shift between the energy levels with the same quantum numbers but with opposite values of
. If the (114) formula was not the same for opposite values of
, we should not be able to get four polynomial radial functions with only one condition which gives the quantification of the energy levels. The Standard Model has a precise answer using vacuum polarization. But the calculation must be revised, both to avoid divergences and to employ the improved wave equation which accounts also for weak interactions [33] [34]. Since the Lamb shift is of the same order as the hyperfine structure coming from the interaction between the magnetic moment of the proton and that of the electron, a true calculation must account for the origin of the magnetic moment of the proton from the waves of the three quarks inside, and for the true potential seen by the electron wave, with both the electric charges and the potential vector of the moving charges.
NOTES
1When we solve the Dirac equation with Darwin’s method, that means with the ad-hoc operators, we get some Legendre polynomials and spherical harmonics. Here, working with
, which is equivalent to working with the Weyl spinors
and
, we get the Gegenbauer polynomials, and it is the degree of these polynomials that gives the needed quantum number.