Corrigendum: Heisenberg uncertainty relation for relativistic electrons (2019 New J. Phys. 21 073036)

The Heisenberg uncertainty relation is derived for relativistic electrons described by the Dirac equation. The standard nonrelativistic lower bound 3 / 2 ÿ is attained only in the limit and the wave function that reproduces this value is singular. At the other end, in the ultrarelativistic limit, the bound is the same as that found before for photons.

In the equations (19) and (23) the mixed terms containing the products of the amplitudes p f , + ( )and p f , -( ) are missing. There was also an additional misprint.
The corrected form of equation (19) where , These changes do not influence our final conclusion since all omitted terms contain the derivatives with respect to angular variables while the functions that minimize the product of uncertainties are spherically symmetric.

Introduction
Recent progress in the generation of relativistic beams of electrons calls for adequate theoretical tools, based on the Dirac equation, that would allow to account precisely for quantum properties of these electrons. The problem raised and answered in this paper is that of Heisenberg uncertainty relation (HUR) for relativistic electrons. It is shown that the relativistic HUR is markedly different from its original nonrelativistic form.
In most monographs [1][2][3][4][5][6][7][8] devoted to relativistic quantum mechanics of electrons the HURs has not been mentioned. In those few places where there is a reference to uncertainty relations in a relativistic theory, the authors [9, 10] follow the early paper by Landau and Peierls [11] giving only a qualitative analysis of limitations imposed by relativity on measurements. In this work we present a complementary approach to uncertainty relations. It could be said that we do what Kennard [12] had done for the nonrelativistic form of the uncertainty relation. He derived an exact mathematical inequality which expressed precisely the heuristic explanation (Bohr-Heisenberg microscope) of the nonrelativistic uncertainty relation. We do the same for relativistic electrons.
In nonrelativistic quantum mechanics one may rely directly on the interpretation of the modulus squared of the wave function y | ( )| r t , 2 and its Fourier transform y |˜( )| p t , 2 as the probability densities in position and momentum space. In relativistic quantum mechanics there is an ambiguity in the definition of the probability density in position space because there are various properties of relativistic particles that may be used to define their position. In particular, we may use the energy distribution, as we have done for photons [13], or the rest mass distribution.
In this work we follow Dirac [3] who postulated (p 258) that 'The square of the modulus of the wave function, summed over the four components, should give the probability per unit volume of the electron being at a certain place'. This quantity is also ( [3] , , and y a ( ) r t , are the components of the Dirac bispinor. Thus, the localization of the electron is determined by its interaction with the electromagnetic field: the electron is where its charge is. A natural definition of the probability density in momentum space is the modulus squared of the Fourier transform of y ( ) r , Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
Thus, the relation between r ( ) r t , r and r ( ) p t , p , based on the Fourier transformation, is the same as in the nonrelativistic quantum mechanics.
From the probability densities in position and in momentum space we construct the standard expressions for the dispersions Δr 2 and Δp 2 : where N 2 is the normalization constant ò ò We omitted the time dependence because the uncertainty relations are always formulated at a fixed time.
A simple textbook proof of the HUR in nonrelativistic quantum mechanics [14,15] can be extended to the relativistic theory. This extension starts from the following non-negative expression built from the Dirac bispinor: After the rearrangement, this inequality becomes a quadratic form in ζ: where we used the expression for Δp 2 in the position representation The requirement that Q(ζ) has at most one real zero (the discriminant cannot be positive) leads to the standard nonrelativistic HUR in three-dimensions Our task is to determine whether in the relativistic theory this inequality is saturated or whether it is a strict inequality. Saturation would mean that there exists a bispinor describing an electronic state which obeys the set of three equations In the nonrelativistic quantum mechanics, these equations are satisfied by a Gaussian wave function but in the relativistic case it is not so simple. Of course, a Gaussian wave function is a solution the equation (12) but such a wave function does not describe a state of an electron but it will also have a part describing the positron. We shall analyze this problem in detail in the next section with the use of a variational procedure.

Evaluation of the variational functional
The variational procedure is simplified if we put á ñ = r 0 and á ñ = p 0 in the definitions (5) and (6). In the case of Δr 2 the prescription is simple. We just choose the origin of the coordinate system at the center of the probability distribution. In the case of Δp 2 the solution is more subtle because the space of momenta is not homogeneous and we cannot choose its origin at will. We may, however, make the following replacement: y y  á ñ ( ) ( · ) ( ) r p r r exp i which does not restrict the freedom in the choice of the variational trial function. Such a replacement does not change Δr 2 but it eliminates á ñ p in (10), so that we obtain: In this way we can eliminate completely the average values of position and momentum from the formulas for the dispersions.
We will express now the product Δr 2 Δp 2 in terms of the independent degrees of freedom that can be subjected to unrestricted variations. The solution of the differential equation obtained from the variational procedure will determine the minimal value of D D r p s , of momentum variables and the spin index s are not subjected to any restrictions; they represent the independent degrees of freedom of an electron moving in free space. They form a representation of the Poincaré group (inhomogeneous Lorentz group) for massive spin 1/2 particles, as has been described by Wigner [16].
Two orthonormal bispinors a ( ) p u s , in the Weyl (chiral) representation of γ matrices [17] will be chosen in the form: The corresponding formula for Δr 2 , Double arrows in (19) denote the antisymmetrized derivative, i.e. ¶ « = ¶ - ¶ f g f g g f x x x . We now express these formulas in spherical coordinates to simplify the variational equations. and we explicitly indicated the dependence of the function f on the parameter d. Thus, our variational equation has the form of the eigenvalue problem for the radial wave equation with a modified harmonic oscillator potential. The lowest value of γ (the energy of the ground state) will be the best bound in the uncertainty relation. Should we allow for the angular dependence, the effective potential would pick up the centrifugal term l(l+1)/q 2 . This would definitely increase the ground state energy.
The eigenvalue problem (25) has no analytic solution but considering the potential V(q), as a function of the parameter d we find (see figure 1) that the smaller the value of d, the deeper the potential. Therefore, we will obtain lower values of the eigenvalue γ for smaller values of d. In the formal nonrelativistic limit (when  ¥ m or  ¥ c ) we have d→0 and the variational equation is that of the harmonic oscillator in 3D and we obtain the nonrelativistic result γ 0 =3/2. The physical limit of d→0 is achieved when either the dispersion in p tends to 0 or the dispersion in r tends to infinity. Of course, these limits cannot be reached for regular, integrable functions. Thus, in relativistic quantum mechanics of electrons the HUR has the form of the strict inequality: 3 0 2 2 In contrast to the nonrelativistic case, this inequality is never saturated. There is no square integrable function that would give the equality in (30), but one may get arbitrarily close to 3 2 . An exact form of the relativistic uncertainty relation can be obtained also in the other extreme case; in the ultra-relativistic limit, when  ¥ d . In this case, the variational equation again has an analytic solution. The wave function of the lowest energy state and the corresponding value of γ are: Quite unexpectedly, this value of g ¥ is the same as for photons [13], even though the definition of Δr 2 for photons is different; it is based on the energy density and not on the charge density.

Two explicit examples
The general results presented here are well illustrated by two analytic solutions of the Dirac equation: Ground state of the electron in the Coulomb potential and the Dirac hopfion (discovered by us recently [19]).

Hydrogen-like ion
The Dirac bispinor describing the ground state of the electron in a hydrogen-like ion has the form ( [5] p 55): We use the standard notation for the hydrogenic wave functions, hoping that the use of γ in this subsection will not be confused with its use in the uncertainty relations. For the bispinor (33) both dispersions Δr 2 and Δp 2 can be analytically evaluated and the HUR reads: parameter d in our calculations. The parameter d allows one to compare directly the values of the two dispersions Δr 2 and Δp 2 . One may tell that we are in the nonrelativistic regime, when the wave function is more squeezed in the momentum space than in the position space. The nonrelativistic limit is reached in two distinct physical situations, when D  ¥ r 2 or Δp 2 →0. The dependence of the lower bound on d in the uncertainty relation, as illustrated in figure 1 stands in contrast to the standard nonrelativistic quantum mechanics where all Gaussians saturate the uncertainty relation.