^{1}

^{2}

In this paper, we examine quantum systems with relativistic dynamics. We show that for a successful description of these systems, the application of Galilei invariant nonrelativistic Hamiltonian is necessary. To modify this Hamiltonian to relativistic dynamics, we require precise relativistic kinetic energy operators instead of nonrelativistic ones for every internal (Jacobi) coordinate. Finally, we introduce and investigate the Schrödinger equation with relativistic dynamics for two-particle systems with harmonic oscillator and Coulomb potentials.

The definition of relativistic quantum mechanical momentum operator, introduced in [

The core of SRT are two postulates [

- the mass of particle m (sometimes called as rest mass) is constant in all reference frames;

- the relativistic momentum of particle, moving in reference frame with velocity v, is p = γ m v , where γ = ( 1 − v 2 / c 2 ) − 1 / 2 is the Lorentz factor;

- the total energy of the particle equals E = γ m c 2 , and the kinetic energy is T = ( γ − 1 ) m c 2 ;

- the general relation between radius vector r and time t of the event as seen in the reference frame S and the radius vector r ′ and time t ′ of the same event in the reference frame S ′ , which is moving with uniform velocity relative to S, is defined by Lorentz transformation, which implies invariance of the space-time interval, defined as the scalar product of these four-vectors:

s 2 = ( c t ) 2 − r 2 = ( c t ′ ) 2 − ( r ′ ) 2 ; (1)

- the mass invariance with respect to transformations from one to another inertial reference frame of the above type also defines invariance of the four- vector, composed of relativistic energy and relativistic momentum:

m 2 c 4 = E 2 − ( p c ) 2 = ( E ′ ) 2 − ( p ′ c ) 2 . (2)

The main problem is―how and to what extent can these statements be applied in the theory of quantum systems?

The base for modifications is the many-particle Schrödinger equation

i ℏ ∂ ∂ t Ψ ( ρ 1 , ρ 2 , ⋯ , ρ N , t ) = H ^ ( ρ 1 , ρ 2 , ⋯ , ρ N ) Ψ ( ρ 1 , ρ 2 , ⋯ , ρ N , t ) , (3)

written in a laboratory reference system. Here ρ j marks the set of one-particle variables (spatial radius vector r j and internal degrees of freedom τ j like mass, charge, spin, isospin and so on for every particle). The Hamiltonian is independent of time, so the formal solution of (3) is

Ψ ( ρ 1 , ρ 2 , ⋯ , ρ N , t ) = exp ( − i ℏ H ^ t ) Φ ( ρ 1 , ρ 2 , ⋯ , ρ N ) . (4)

The stationary states of (3) correspond to the functions Φ , defined as the eigenfunctions of the Hamiltonian:

H ^ Φ = E Φ . (5)

The probability density is therefore independent of time. The quantum system is stable with respect to breakup into composing particles if the corresponding eigenvalue E is negative (−E equals the binding energy of the system). We will consider here the application of relativistic dynamics to quantum systems defined this way.

The introduced presentation gives the simplest way to ensure symmetry of the Hamiltonian with respect to permutations of one-particle degrees of freedom, because quantum systems or their subsystems are composed of identical particles. Moreover, these particles are fermions, so the eigenfunction of the Hamiltonian has to be antisymmetric with respect to permutations of the variables of identical fermions. It is especially simple to project the wave function, written in single-particle variables, to the corresponding subspace of Hilbert space. Various methods for antisymmetrization of wave functions, dependent on single-particle coordinates, are thoroughly developed from application of Slater determinants or Racah coefficients of fractional parentage [

The relativistic dynamics introduction for Schrödinger Equation (3), if even possible, is an immensely difficult problem. The equation is written in a laboratory reference system. Due to the necessary equivalence of all inertial reference systems, the equation has to be invariant with respect to time and space variables for the Lorentz transformation of every particle, keeping invariant the length of the four vector, defined by Equation (1). Even for one free particle such an equation (known as the Dirac equation) is possible only after defining the quantum mechanical momentum operator as corresponding to the relativistic momentum. However, this operator does not have any dependence on velocity, necessary for relativistic momentum, so the corresponding equation can not be considered completely relativistic.

Moreover, even in the nonrelativistic approximation the wave functions defined in this way are not acceptable for a precise description of the quantum system. These functions do not satisfy Galilei invariance and, in many cases, contains uncontrolled center of mass excitations. For the best illustration of this statement, let us consider a two-particle system. The potential energy operator is translationally invariant, i.e. independent of the center of mass coordinate, hence the only problem is in the kinetic energy operator. The corresponding wave function written in one-particle spatial variables is not under application, because the elimination of the center of mass kinetic energy operator in a two-particle Hamiltonian is an easy task. Elimination of the center of mass kinetic energy operator is not complicated for a system of N particles either. However, after this the Hamiltonian of the quantum system appears as a function of 3 ( N − 1 ) internal spatial variables and N sets of individual degrees of freedom of particles and antisymmetrization of eigenfunctions of this operator becomes a very complex problem.

Therefore, the solution of the nonrelativistic Schrödinger equation for a quantum system with internal Hamiltonian is a real challenge, but only these wave functions can give an adequate description of the quantum system. The experimental investigation of a quantum system is based on observing the behavior of a system after excitation by an external agent. Every act of this kind excites the center of mass of the system and, at the same time, the internal degrees of freedom. Unfortunately, control of the center of mass movement after a collision is very problematic. Among the attempts to minimize the energy loss due to center of mass excitation most notable are the Mösbauer effect and colliding beams accelerators. To avoid gamma ray energy loss due to nuclear recoil, Mösbauer used nuclei bound in a solid. On the other hand, in an experiment involving head-on colliding beams with identical particles and equal energies, the kinetic energy of the center of mass is zero before, and hence after, the collision. Either way, as the exciting particle leaves the range of the interaction force, the center of mass velocity is constant, so this point can serve as the origin of an inertial reference frame. All particles, composing the system, move in the same inertial frame, hence any transformations of radius vector and time coordinate or energy and momentum of given particle between different inertial systems are not necessary. Both the nonrelativistic and modified relativistic (introduced in [

The transformation of the initial Schrödinger Equation (3), present in a laboratory system, to the aforementioned form, starts with the nonrelativistic Hamiltonian transformation to center of mass and internal variables. The canonical choice is the Jacobian coordinates, defined as in [

∑ j = 1 N 1 m j Δ j = ∑ β = 0 N − 1 1 μ β Δ ( β ) . (6)

Here β is the subscript of the Jacobi variable ξ β , the corresponding Laplasian is Δ ( β ) , and reduced mass equals μ β . The Jacobi coordinate ξ 0 and mass μ 0 are defined as the system’s center of mass radius vector and mass:

ξ 0 = 1 μ 0 ∑ j = 1 N m j r j , μ 0 = ∑ j = 1 N m j . (7)

The Hamiltonian appears with separated center of mass and internal variables:

H ^ ( ρ 1 , ρ 2 , ⋯ , ρ N ) = H ^ c m ( ξ 0 ) + H ^ i n t ( ξ 1 , ξ 2 , ⋯ , ξ N − 1 , τ 1 , τ 2 , ... , τ N ) . (8)

Furthermore, to transform to a relativistic form, we need to take into account the rest energies of the particles in the Hamiltonian:

∑ j = 1 N m j c 2 = μ 0 c 2 . (9)

Now the quantum system’s Schrödinger Equation (3) takes the following form:

( i ℏ ∂ ∂ t − H ^ c m − μ 0 c 2 − H ^ int ) Ψ ( ξ 0 , ξ 1 , ⋯ , ξ N − 1 , τ 1 , τ 2 , ⋯ , τ N , t ) = 0 (10)

The variables of the complete Hamiltonian are separated, so the simplest solutions of this equation are products of functions, dependent on different variables:

Ψ ( ξ 0 , ξ 1 , ⋯ , ξ N − 1 , τ 1 , τ 2 , ⋯ , τ N , t ) = Ψ c m ( ξ 0 , t ) Ψ int ( ξ 1 , ξ 2 , ⋯ , ξ N − 1 , τ 1 , τ 2 , ⋯ , τ N ) , (11)

where the introduced functions are eigenfunctions of the following operators, corresponding to the same eigenvalue E:

( i ℏ ∂ ∂ t − H ^ c m − μ 0 c 2 − E ) Ψ c m ( ξ 0 , t ) = 0 (12)

and

[ H ^ i n t − E ] Ψ i n t ( ξ 1 , ξ 2 , ⋯ , ξ N − 1 , τ 1 , τ 2 , ⋯ , τ N ) = 0. (13)

The eigenvalue E, by definition of the internal Hamiltonian, equals the binding energy of system with opposite sign. Thus μ 0 c 2 + E is the reduced, slightly smaller, rest energy of the bound system μ ¯ 0 c 2 . Necessarily, the kinetic energy operator has to be expressed with the reduced rest energy of system, i.e. as

H ^ c m = − ( ℏ c ) 2 μ ¯ 0 c 2 Δ ( 0 ) . (14)

We go on to further modify the center of mass Schrödinger Equation (12), taking into account the relativistic dynamics, by applying the technique introduced in [

It is necessary to modify the kinetic energy operators for all Jacobi variables for the introduction of relativistic dynamics into the internal Hamiltonian. The natural subject for the application of these ideas is the two-particle system.

The internal Hamiltonian of a two-particle system, whose interaction is defined by a local potential, dependent on the radius vector r = r 2 − r 1 (here r 2 and r 1 are radius vectors, defining locations of particles in laboratory system) and on internal degrees of freedom of particles τ 1 and τ 2 , is

H ^ 0 ( r , τ 1 , τ 2 ) = T ^ ( 0 ) ( r , μ ) + V ^ ( r , τ 1 , τ 2 ) . (15)

Here T ^ ( 0 ) is the nonrelativistic kinetic energy operator

T ^ ( 0 ) ( r , μ ) = − ℏ 2 2 μ Δ = − ( ℏ c ) 2 2 μ c 2 [ 1 r ∂ 2 ∂ r 2 r − L ^ 2 ( θ φ ) r 2 ] , (16)

with the reduced mass of system μ = m 1 m 2 / ( m 1 + m 2 ) and the orbital mo- mentum operator L ^ 2 ( θ φ ) . The eigenvalues of the stationary Schrödinger equation with this Hamiltonian due to symmetries of potential are identified at least by such quantum numbers as parity π and total angular momentum J . The simpler is the potential, the larger is the set of exact quantum numbers, characterizing the eigenfunction of this Hamiltonian, i.e. the state of quantum system.

The Hamiltonian with relativistic dynamics

H ^ ( r , τ 1 , τ 2 ) = T ^ ( r , μ ) + V ^ ( r , τ 1 , τ 2 ) , (17)

is defined with the relativistic kinetic energy operator

T ^ ( r , μ ) = μ c 2 [ ( 1 − 2 T ^ ( 0 ) ( r , μ ) μ c 2 ) − 1 / 2 − 1 ] . (18)

This Hamiltonian is translationally invariant, describing the relative movement of particles in an inertial reference frame, whose origin is at the center of mass of system. The eigenfunctions of the Hamiltonian (17), as in the nonrelativistic case, can be identified by the same set of exact quantum numbers, because the operators H ^ 0 ( r τ 1 τ 2 ) and H ^ ( r , τ 1 , τ 2 ) share the common set of commuting operators.

However, a new problem appears. In general, the eigenfunctions of the nonrelativistic Hamiltonian and the Hamiltonian with relativistic dynamics are different due to presence of potential. For a free particle they coincide, but here the possibility to define a relation of eigenvalues of these Hamiltonians is not as simple. Let us restrict the investigation by applying wave functions, obtained solving the nonrelativistic Schrödinger equation. The nonrelativistic energy expectation value equals

E ( 0 ) = 〈 T ^ ( 0 ) ( r , μ ) 〉 + 〈 V ^ ( r , τ 1 , τ 2 ) 〉 (19)

The corrected value for this wave function is

E = 〈 T ^ ( r , μ ) 〉 + 〈 V ^ ( r , τ 1 , τ 2 ) 〉 . (20)

For a relativistic kinetic energy expectation value evaluation we will apply three different approximations.

The first approximation is based on a series expansion of the relativistic kinetic energy operator:

T ^ ( r , μ ) = T ^ ( 0 ) ( r , μ ) + 3 ( T ^ ( 0 ) ( r , μ ) ) 2 2 μ c 2 + ⋯ , (21)

where we can successfully evaluate the expectation value of ( T ^ ( 0 ) ( r , μ ) ) 2 and therefore the kinetic energy with first order correction

T ( 1 ) = 〈 T ^ ( 0 ) 〉 + 3 2 μ c 2 〈 ( T ^ ( 0 ) ) 2 〉 (22)

can be successfully evaluated. This result is precise enough for quantum systems with dynamics not far from nonrelativistic.

The second approximation is

T ( 2 ) = μ c 2 [ ( 1 − 2 〈 T ^ ( 0 ) 〉 μ c 2 ) − 1 / 2 − 1 ] . (23)

With this approximation, we ensure, according to the main idea of SRT, the expectation value of the nonrelativistic kinetic energy to be less than μ c 2 / 2 , same as for a free particle. It is the really original result, consistent with the basic postulate of SRT stating that the relative velocity of particles cannot be larger than speed of light at any interaction. However, this result is only exact in the case when the eigenfunctions of the kinetic energy operator are applied for the calculation of the expectation value. For a quantum system this approximation at a low nonrelativistic kinetic energy gives a lower bound for relativistic kinetic energy expectation value due to a relation, valid in an arbitrary basis:

〈 ( T ^ ( 0 ) ) 2 〉 ≥ 〈 T ^ ( 0 ) 〉 2 . (24)

We obtain the third evaluation of relativistic kinetic energy in a given state of a quantum system by applying spectral representation of the corresponding operator. Both kinetic energy operators have the same set of eigenfunctions. The internal Schrödinger equation for the nonrelativistic operator is

T ^ ( 0 ) ( r ) ψ k ( r ) = T k ( 0 ) ψ k ( r ) . (25)

We can rewrite it as

Δ ψ k ( r ) = − k 2 ψ k ( r ) , (26)

where

k 2 = 2 μ c 2 T k ( 0 ) ( ℏ c ) 2 (27)

and

ψ k ( r ) = 1 ( 2 π ) 3 / 2 e i ( k ⋅ r ) . (28)

The spectral representation of this operator becomes

T ^ ( 0 ) ( r , r ′ ) = ∫ d k ψ k ( r ) T k ( 0 ) ψ k + ( r ′ ) , (29)

where the kernel of this operator

T k ( 0 ) = ( ℏ c ) 2 2 μ c 2 k 2 . (30)

Due to equivalence of eigenfunctions, the spectral representation of the relativistic kinetic energy operator is

T ^ ( r , r ′ ) = ∫ d k ψ k ( r ) T k ψ k + ( r ′ ) . (31)

The kernel of this operator is defined by Equation (18):

T k = μ c 2 [ ( 1 − ( ℏ c μ c 2 ) 2 k 2 ) − 1 / 2 − 1 ] . (32)

The expectation value of the nonrelativistic kinetic energy in a given state of the internal two-particle Hamiltonian, whose eigenfunction is Ψ ( r ) , equals

T ( 0 ) = ∫ d r ∫ d r ′ Ψ + ( r ) T ^ ( 0 ) ( r , r ′ ) Ψ ( r ′ ) , (33)

while the expectation value of the relativistic kinetic energy is

T ( 3 ) = ∫ d r ∫ d r ′ Ψ + ( r ) T ^ ( r , r ′ ) Ψ ( r ′ ) . (34)

From operator definitions given by Equations ((29), (31) and (28)), alternative definitions of expectation values follow:

T ( 0 ) = ∫ d k Ψ + ( k ) T k ( 0 ) Ψ ( k ) (35)

and

T ( 3 ) = ∫ d k Ψ + ( k ) T k Ψ ( k ) , (36)

where Ψ ( k ) is Fourier image of wave function Ψ ( r ) :

Ψ ( k ) = ∫ d r Ψ ( r ) ψ k + ( r ) . (37)

One of the best known problems of nonrelativistic quantum mechanics, having successful applications, is a system of two particles with a harmonic oscillator potential. The internal Schrödinger equation for this system is

( − ( ℏ c ) 2 2 μ c 2 Δ + 1 2 μ ω 2 r 2 ) Ψ = E ( 0 ) Ψ , (38)

where μ is the reduced mass of the interacting particles. The equation has square integrable solutions

Ψ n l m ( r ) = N n l x l e x p ( − x 2 / 2 ) L n l + 1 / 2 ( x 2 ) Y l m ( ϑ φ ) (39)

only for eigenvalues E n l ( 0 ) = ℏ ω ( 2 n + l + 3 / 2 ) . Here the argument of the radial function is x = r / b , with b = ℏ / ( μ ω ) . n = 0 , 1 , 2 , ⋯ equals the number of nodes of the radial function and the normalization constant is

N n l 2 = 2 n ! b 3 Γ ( n + l + 3 / 2 ) . (40)

To evaluate the expectation value of the corresponding relativistic kinetic energy operator 〈 T ^ 〉 n l , we require the expectation value of the nonrelativistic kinetic energy. According to Virial theorem for a harmonic oscillator potential

〈 T ^ ( 0 ) 〉 n l = E n l ( 0 ) / 2 = 〈 V ^ 〉 n l , (41)

hence

E n l ( 0 ) μ c 2 = β ( 2 n + l + 3 / 2 ) . (42)

Here and in the following expressions, for the sake of simplicity and comparison with different approximations, we use a dimension-free parameter

β = ℏ ω μ c 2 . (43)

The first order correction of the nonrelativistic value, defined in Equation (22), simplifies due to Virial theorem:

〈 ( T ^ ( 0 ) ) 2 〉 n l = 〈 V ^ 2 〉 n l = ( ℏ ω 2 ) 2 〈 x 4 〉 n l .

Thus, the value with the first order correction taken into account equals

E n l ( 1 ) μ c 2 = β ( 2 n + l + 3 / 2 ) + 3 8 β 2 [ ( 2 n + l + 3 / 2 ) 2 + ( 2 n + l + 3 / 2 ) + 2 n ( n + l + 1 / 2 ) ] . (44)

The second approximation, defined by Equation (23), gives:

E n l ( 2 ) μ c 2 = [ 1 − β ( 2 n + l + 3 / 2 ) ] − 1 / 2 − 1 + β ( 2 n + l + 3 / 2 ) / 2 . (45)

In a nonrelativistic approximation ( β ≪ 1 ) , this expression leads to E n l ( 2 ) = E n l ( 0 ) . The new result of relativistic kinematics is the existence of a restricted set of eigenvalues, satisfying the condition

2 n + l + 3 / 2 < 1 β and β < 2 / 3 . (46)

Finally, for the expectation value of the third approximation we need the Fourier image of the eigenfunction, defined by Equation (39). This problem is not complicated, since the transformation operator between the Hamiltonians of the harmonic oscillator in coordinates and momenta representations is well-known. These Hamiltonians, written in dimensionless variables, are correspondingly

H ^ ( r ) = ℏ ω 2 ( − Δ x + x 2 ) , (47)

H ^ ( p ) = ℏ ω 2 ( − Δ y + y 2 ) . (48)

Here x = r / b and y = p b / ℏ . The mentioned above transformation operator is

1 ( 2 π ) 3 / 2 e i ( x ⋅ y ) . (49)

Therefore, the eigenfunctions of these operators are connected by this transformation. The Fourier image of the original function is

Ψ n l m ( p ) = N ¯ n l y l e x p ( − y 2 / 2 ) L n l + 1 / 2 ( y 2 ) Y l m ( ϑ ¯ φ ¯ ) , (50)

where y = p / d , d = ℏ / b = ℏ ω μ and y = ( y , ϑ ¯ , φ ¯ ) . The normalization constant equals

N ¯ n l 2 = 2 n ! d 3 Γ ( n + l + 3 / 2 ) . (51)

This wave function with the kernel, defined in Equation (30) produces the expectation value of the nonrelativistic kinetic energy, equal E n l ( 0 ) / 2 . According to the definition we express in Equation (36), the third approximation for relativistic kinetic energy equals

T n l ( 3 ) μ c 2 = n ! Γ ( n + l + 3 / 2 ) ∫ 0 1 / β d z ⋅ z l + 1 / 2 [ ( 1 − β z ) − 1 / 2 − 1 ] e − z [ L n l + 1 / 2 ( z ) ] 2 , (52)

and thus the approximation for total energy is

E n l ( 3 ) μ c 2 = T n l ( 3 ) μ c 2 + β ( 2 n + l + 3 / 2 ) / 2 . (53)

The system of a charged particle moving in a Coulomb field has many interesting applications. We consider a particle with charge − e (electron or muon) moving in a Coulomb well, created by an atomic nucleus, whose mass equals M and charge is Z e . The nonrelativistic Schrödinger equation for this system is

( ( ℏ c ) 2 2 μ c 2 Δ + α Z ℏ c r + E n ( 0 ) ) ψ n l m ( r θ φ ) = 0. (54)

Here μ is the reduced mass of light particle and atomic nucleus, α is the fine structure constant. The square integrable solutions of this equation:

ψ n l m ( r θ φ ) = N n l x l + 1 exp ( − x / 2 ) L n − l − 1 2 l + 1 ( x ) Y l m ( ϑ φ ) , (55)

where x = 2 r / n a , n = l + 1 , l + 2 , ⋯ , a = ℏ c / α Z μ c 2 ,

N n l 2 = ( n − l − 1 ) ! a n 2 ( n + l ) ! , (56)

exists only for eigenvalues

E n ( 0 ) = − μ c 2 2 ( α Z n ) 2 . (57)

Virial theorem for Coulomb potential gives

〈 V 〉 n l = 2 E n , 〈 T ^ ( 0 ) 〉 n l = − E n . (58)

The first order correction, defined by Equation (22), is

E n l ( 1 ) μ c 2 = − 1 2 ( α Z n ) 2 + 3 2 ( α Z n ) 4 ( n l + 1 / 2 − 3 4 ) . (59)

The correction, defined by Equation (23), equals

E n l ( 2 ) μ c 2 = [ 1 − ( α Z n ) 2 ] − 1 / 2 − ( α Z n ) 2 − 1. (60)

Following from this equation, the condition of bound states existence in a Coulomb well looks like

n > α Z . (61)

For existing stable atomic nuclei α Z < 1 , so relativistic kinematics do not restrict the number of bound states in a Coulomb well. On the other hand, significant relativistic corrections can arise for Hydrogen-like ions of heavy atoms.

Finally, for the evaluation of the correction, defined by Equation (36) we require the Fourier image of wave function ψ n l m ( r ) . Corresponding expre- ssions are complex enough, so we will only present results, obtained for the ground state function

ψ 100 ( r ) = e − r / a π a 3 . (62)

The Fourier image of this function is

ψ 100 ( k ) = ( 2 a ) 3 / 2 π [ 1 + ( a k ) 2 ] 2 . (63)

This function with kernel (30) gives a nonrelativistic kinetic energy, equal − E 10 ( 0 ) = μ c 2 ( α Z ) 2 / 2 . With the relativistic kernel (32), we find the corrected energy value to be

E 10 ( 3 ) μ c 2 = 32 π ∫ 0 1 α Z x 2 ( 1 + x 2 ) 4 [ ( 1 − ( α Z x ) 2 ) − 1 / 2 − 1 ] − ( α Z ) 2 . (64)

The introduced method of relativistic corrections of bound state energies of two- particle systems with harmonic oscillator potential presents the eigenvalues as functions of one dimensionless parameter β < 2 / 3 .

β | |||||
---|---|---|---|---|---|

0.001 0.01 0.1 0.2 0.3 0.4 0.5 0.6 0.666 | 0.0015 0.0150 0.1500 0.3000 0.4500 0.6000 0.7500 0.9000 0.9990 | 0.00150 0.01514 0.16406 0.35625 0.57656 0.82500 1.10156 1.40625 1.62275 | 0.00150 0.01509 0.15965 0.34523 0.57340 0.88114 1.37500 2.61228 31.1223 | 0.00150 0.01515 0.17199 0.42617 0.66631 0.82711 0.92913 0.99864 1.03511 | 0.00150 0.01495 0.14531 0.28125 0.40781 0.52500 0.63281 0.73125 0.79108 |

corresponding to relativistic momentum, i.e. when the kinetic energy expansion

T = m 2 c 4 + p 2 c 2 − m c 2 = p 2 2 m − p 4 8 m 3 c 2 + ⋯ (65)

is approximated as

T ¯ = m 2 c 4 + p 0 2 c 2 − m c 2 = p 0 2 2 m − p 0 4 8 m 3 c 2 + ⋯ (66)

From this energy expression, found in numerous textbooks and basic for Klein-Gordon and Dirac equations, it follows that

E 2 = γ 2 m 2 c 4 = m 2 c 4 + p 0 2 c 2 = ( 1 + ( p 0 m c ) 2 ) m 2 c 4 , (67)

therefore it is not the completely relativistic expression for energy, but rather a crude approximation for the Lorentz factor

γ = 1 1 − v 2 / c 2 (68)

by 1 + v 2 / c 2 . Therefore, the accuracy of this presentation does not exceed ( v 2 / c 2 ) 2 . Moreover, this approximation results in the absence of relativistic momentum and energy dependence of v in the limit v → c . The corresponding expectation value for harmonic oscillator potential, given in the last column of

E n l ( − ) μ c 2 = β ( 2 n + l + 3 / 2 ) − 1 8 β 2 [ ( 2 n + l + 3 / 2 ) 2 + ( 2 n + l + 3 / 2 ) + 2 n ( n + l + 1 / 2 ) ] . (69)

The problem is to understand how sensitive to relativistic corrections are the energies of different quantum systems, performing vibrational movement. The first one among them is nuclear motion in diatomic molecules. The parameter β takes a larger value at a small reduced mass of nuclei, hence the best system for the evaluation is the Hydrogen molecule H_{2}. Corresponding reduced mass equals half the proton mass. The molecular internal energy is approximately the sum of the electromagnetic energy, the vibrational energy ( v + 1 / 2 ) ℏ ω and the rotational energy [

For the ground state ( n = 1 , l = 0 ) of Hydrogen-like ions, corresponding results are present in

1 | 13.6057 | 13.6030 | 13.6051 | 13.6030 | 13.6066 | 13.5984 |

2 | 54.4228 | 54.3793 | 54.4141 | 54.3787 | 54.4373 | 54.4178 |

3 | 122.4512 | 122.2313 | 122.4072 | 122.2266 | 122.5246 | 122.4544 |

4 | 217.6911 | 216.9955 | 217.5519 | 216.9773 | 217.9229 | 217.7186 |

5 | 340.1423 | 338.4442 | 339.8023 | 338.3912 | 340.7084 | 340.2260 |

6 | 489.8049 | 486.2838 | 489.0996 | 486.1588 | 490.9787 | 489.9932 |

40 | 21,769.11 | 14,813.70 | 20,271.28 | 15,924.90 | 24,087.58 | 22,236.68 |

60 | 48,980.49 | 13,768.75 | 40,584.86 | 29,189.06 | 60,717.74 | 51,515.58 |

70 | 66,667.89 | 1,433.78 | 49,936.25 | 38,807.72 | 88,412.60 | 71,574.80 |

ionization potentials are expressed in electronvolts, equal to the energies with opposite sign. Experimental data, taken from [

An interesting conclusion follows from these results about the Hydrogen atom ionization potential. As usual, the starting point is the nonrelativistic value J 10 ( 0 ) = 13. 6057 eV with the traditional first order correction, predicting relativistic kinetic energy smaller than the nonrelativistic one, i.e. J 10 ( − ) = 13. 6066 eV . Comparing with the experimental data, this result is worse than zero approximation. The possibility of correction in approximation, proportional ( α Z ) 4 is found―it is the so called Darvin term [

Our investigation has shown that successful relativistic dynamics application needs Galilei invariant formalism for the description of a quantum system. Galilei invariance, basic for relativity, appears as a bridge between relativity and quantum mechanics for relativistic quantum systems theory. Internal Hamiltonian of the quantum system corresponds to the nonrelativistic dynamics of particles, performing mutual movement, when all particles are present in the same inertial system. No time coordinate and Lorentz transformations, describing relations between movements in different inertial systems, are necessary for these particles. The only modification of such a Hamiltonian, taking into account relativistic dynamics, is the introduction of the precise relativistic kinetic energy operator for every internal spatial variable instead of the corresponding nonrelativistic one. These operators, defined in [

The introduced internal relativistic Hamiltonian in limit of small (in comparison with μ c 2 / 2 ) kinetic energies of particles appears as the nonrelativistic internal Hamiltonian. Both Hamiltonians for a system of identical particles are symmetrical with respect to particles’ spatial coordinates and intrinsic degrees of freedom permutations. However, the antisymmetrization of the eigenfunctions of these Hamiltonians is problematic, because permutations of particles coordinates produce linear transformations of Jacobi variables. Traditional methods of this operation are not applicable. It is necessary to apply original methods, whose historical overview and the simplest known solution of the problem is present in [

Our investigation has shown that known potentials of particles interactions are weak enough to produce the quantum system with significant relativity. As a result, the relativity is soft, so it is acceptable to apply the introduced approximate methods for correction of nonrelativistic energies of quantum systems bound states.

The really interesting quantum system, whose solution needs more sophisticated relativistic effects consideration, is the nucleon. For weakly bound states (such as proton and neutron) to exist in a deep confining well, significant kinetic energy is necessary, a few times larger than the relativistic limit. Therefore, for the solution of this problem, application of relativistic theory is very promising.

Kamuntavičius, G.P. and Kamuntavičius, G. (2017) Relativistic Dynamics of a Quantum System. Journal of Applied Mathematics and Physics, 5, 1478-1492. https://doi.org/10.4236/jamp.2017.57121