Comparing the Schr\"odinger and Dirac Descriptions of an Electron in a Uniform Magnetic Field

TL;DR

This paper compares non-relativistic and relativistic quantum descriptions of an electron in a uniform magnetic field, highlighting differences and relativistic effects using Schrödinger and Dirac equations.

Contribution

It provides a detailed comparison of Schrödinger and Dirac models for an electron in a magnetic field, emphasizing relativistic modifications.

Findings

Relativistic effects alter the electron's behavior in a magnetic field.

Differences between Schrödinger and Dirac descriptions are systematically analyzed.

Relativistic phenomena impact the quantum state evolution of the electron.

Abstract

In this article we present a detailed description of an electron in a uniform magnetic field evolving under the Schr\"odinger equation using ladder operators. Based on this analysis, we describe the same physical system using the Dirac equation, known from relativistic quantum mechanics. The main differences between these two quantum mechanical approaches are discussed and we observe specifically how the relativistic phenomena modify the description of this particular quantum system.

Tables

Table 1. Table 1: Radial eigenfunction common to the two-dimensional harmonic oscillator and the observable L ^ z subscript ^ 𝐿 𝑧 \hat{L}_{z} , for the first six values of n 𝑛 n . Negative angular momentum projection m l subscript 𝑚 𝑙 m_{l} entails negative azimuthal phase.

$n$	$m_l$	$F_{n,m_l}(\rho ,\phi )$
0	0	$F_{0,0}=\frac{\beta }{\sqrt{\pi }}{e}^{-{\beta }^2{\rho }^2/2}$
1	$1$	$F_{1,1}=\frac{\beta }{\sqrt{\pi }}\left(\beta \rho \right)e^{-{\beta }^2{\rho }^2/2}{e}^{i\phi }$
2	$2$	$F_{2,2}=\frac{\beta }{\sqrt{2\pi }}{\left(\beta \rho \right)}^2{e}^{-{\beta }^2{\rho }^2/2}{e}^{i2\phi }$
	0	$F_{2,0}=\frac{\beta }{\sqrt{\pi }}\left[{\left(\beta \rho \right)}^2-1\right]e^{-{\beta }^2{\rho }^2/2}$
3	$3$	$F_{3,3}=\frac{\beta }{\sqrt{6\pi }}{\left(\beta \rho \right)}^3{e}^{-{\beta }^2{\rho }^2/2}{e}^{i3\phi }$
	$1$	$F_{3,1}=\frac{\beta }{\sqrt{2\pi }}\left[{\left(\beta \rho \right)}^3-2\left(\beta \rho \right)\right]e^{-{\beta }^2{\rho }^2/2}{e}^{i\phi }$
4	$4$	$F_{4,4}=\frac{\beta }{2\sqrt{6\pi }}{\left(\beta \rho \right)}^4{e}^{-{\beta }^2{\rho }^2/2}{e}^{i4\phi }$
	$2$	$F_{4,2}=\frac{\beta }{\sqrt{6\pi }}\left[{\left(\beta \rho \right)}^4-3{\left(\beta \rho \right)}^2\right]e^{-{\beta }^2{\rho }^2/2}{e}^{i2\phi }$
	0	$F_{4,0}=\frac{\beta }{\sqrt{4\pi }}\left[{\left(\beta \rho \right)}^4-4{\left(\beta \rho \right)}^2+2\right]e^{-{\beta }^2{\rho }^2/2}$
5	$5$	$F_{5,5}=\frac{\beta }{2\sqrt{30\pi }}{\left(\beta \rho \right)}^5{e}^{-{\beta }^2{\rho }^2/2}{e}^{i5\phi }$
	$3$	$F_{5,3}=\frac{\beta }{2\sqrt{6\pi }}\left[{\left(\beta \rho \right)}^5-4{\left(\beta \rho \right)}^3\right]e^{-{\beta }^2{\rho }^2/2}{e}^{i3\phi }$
	$1$	$F_{5,1}=\frac{\beta }{\sqrt{12\pi }}\left[{\left(\beta \rho \right)}^5-6{\left(\beta \rho \right)}^3+6\left(\beta \rho \right)\right]e^{-{\beta }^2{\rho }^2/2}{e}^{i\phi }$

Equations

$$\hat{H}=\frac{e\hbar B}{2m_{0}}\left(\begin{array}[]{lr}1&0\\ 0&-1\end{array}\right).$$

$$\vec{A}(\vec{r})=\frac{1}{2}\big{(}\vec{B}\times\vec{r}\big{)}.$$

$$\begin{split}\hat{H}&=\frac{1}{2m_{0}}\big{(}\hat{p}+e\vec{A}\big{)}^{2}+\frac{e}{m_{0}}\vec{B}\!\cdot\!\vec{S},\\ &=\frac{-\hbar^{2}}{2m_{0}}\!\nabla^{2}\!+\!\frac{e^{2}A^{2}}{2m_{0}}\!-\!\frac{ie\hbar}{2m_{0}}\big{(}\nabla\!\cdot\!\vec{A}\!+\!\vec{A}\!\cdot\!\nabla\big{)}\!+\!\frac{e}{m_{0}}\vec{B}\!\cdot\!\vec{S}.\end{split}$$

$$\begin{split}\frac{e^{2}A^{2}}{2m_{0}}&=\frac{e^{2}}{8m_{0}}\big{(}\vec{B}\times\vec{r}\big{)}^{2},\\ &=\frac{e^{2}}{8m_{0}}B^{2}r^{2}\sin^{2}(\theta_{Br}),\\ &=\frac{e^{2}}{8m_{0}}\big{(}B^{2}r^{2}-(\vec{B}\!\cdot\!\vec{r})^{2}\big{)},\end{split}$$

$$\nabla\!\cdot\!\big{(}\vec{A}\phi\big{)}=\big{(}\nabla\!\cdot\vec{A}\big{)}\phi+\vec{A}\!\cdot\!\nabla\phi.$$

$$\begin{split}\big{(}\nabla\!\cdot\vec{A}\big{)}\phi&=\frac{1}{2}\nabla\!\cdot\!\big{(}\vec{B}\times\vec{r}\big{)}\phi,\\ &=\frac{1}{2}\big{(}\nabla\times\vec{B}\big{)}\!\cdot\!\vec{r}\phi-\frac{1}{2}\big{(}\nabla\times\vec{r}\big{)}\!\cdot\!\vec{B}\phi,\\ &=0.\end{split}$$

$$\begin{split}\vec{A}\!\cdot\nabla\phi&=\frac{1}{2}\big{(}\vec{B}\times\vec{r}\big{)}\!\cdot\!\nabla\phi,\\ &=\frac{1}{2}\vec{B}\!\cdot\!\big{(}\vec{r}\times\nabla\big{)}\phi.\end{split}$$

$$\begin{split}\vec{A}\!\cdot\nabla\phi&=\frac{i}{2\hbar}\vec{B}\!\cdot\!\big{(}\vec{r}\times\hat{P}\big{)}\phi,\\ &=\frac{i}{2\hbar}\vec{B}\!\cdot\!\hat{L}\phi.\end{split}$$

$$\hat{H}\!=\!\frac{-\hbar^{2}}{2m_{0}}\nabla^{2}\!+\frac{e^{2}}{8m_{0}}\big{(}B^{2}r^{2}\!-\!(\vec{B}\cdot\vec{r})^{2}\big{)}+\frac{e}{2m_{0}}\vec{B}\!\cdot\!\big{(}\hat{L}+2\hat{S}\big{)}.$$

$$\begin{split}\hat{H}&=\frac{-\hbar^{2}}{2m_{0}}\nabla^{2}+\frac{e^{2}B^{2}}{8m_{0}}\big{(}x^{2}+y^{2}\big{)}+\frac{eB}{2m_{0}}\big{(}\hat{L}_{z}+2\hat{S}_{z}\big{)},\\ &=\frac{-\hbar^{2}}{2m_{0}}\frac{\partial^{2}}{\partial z^{2}}+\frac{-\hbar^{2}}{2m_{0}}\left(\frac{\partial^{2}}{\partial x^{2}}+\frac{\partial^{2}}{\partial y^{2}}\right)\\ &\hskip 42.67912pt+\frac{e^{2}B^{2}}{8m_{0}}\big{(}x^{2}+y^{2}\big{)}+\frac{eB}{2m_{0}}\big{(}\hat{L}_{z}+2\hat{S}_{z}\big{)}.\end{split}$$

$$\psi (\vec{r},\sigma )=F(x,y)e^{i{p}_{z}z/\hslash }\Gamma ,$$

$$\Gamma=\left\{\begin{array}[]{ll}{\tiny\left(\begin{array}[]{c}1\\ 0\end{array}\right)},&m_{s}=+1/2\\ {\tiny\left(\begin{array}[]{c}0\\ 1\end{array}\right)},&m_{s}=-1/2\end{array}\right.,$$

$$\underbrace{\left[\frac{-\hbar^{2}}{2m_{0}}\!\left(\!\frac{\partial^{2}}{\partial x^{2}}\!+\!\frac{\partial^{2}}{\partial y^{2}}\!\right)\!+\!\frac{1}{2}m_{0}\omega^{2}\!\big{(}\!x^{2}\!+\!y^{2}\!\big{)}\!\right]}_{\displaystyle H^{\prime}_{xy}}\!F\!=\!E^{\prime}F,$$

$$\begin{array}{rl}{E}^{\prime }& =\hslash \omega (n_x+n_y+1),\\ & =\hslash \omega (n+1),n\ge 0.\end{array}$$

$$F_{n,0},F_{n-1,1},\cdots ,F_{0,n}.$$

$$\begin{array}{rl}{\hat{a}}_x& =\frac{1}{\sqrt{2}}\left(\hat{x}\beta +i\frac{{\hat{P}}_x}{\beta \hslash }\right),{\hat{a}}_x^{†}=\frac{1}{\sqrt{2}}\left(\hat{x}\beta -i\frac{{\hat{P}}_x}{\beta \hslash }\right),\\ {\hat{a}}_y& =\frac{1}{\sqrt{2}}\left(\hat{y}\beta +i\frac{{\hat{P}}_y}{\beta \hslash }\right),{\hat{a}}_y^{†}=\frac{1}{\sqrt{2}}\left(\hat{y}\beta -i\frac{{\hat{P}}_y}{\beta \hslash }\right),\end{array}$$

$$\begin{split}\hat{x}&=\frac{1}{\beta\sqrt{2}}\big{(}\hat{a}_{x}+\hat{a}^{\dagger}_{x}\big{)},\qquad\hat{P}_{x}=-i\frac{\hbar\beta}{\sqrt{2}}\big{(}\hat{a}_{x}-\hat{a}^{\dagger}_{x}\big{)},\\ \hat{y}&=\frac{1}{\beta\sqrt{2}}\big{(}\hat{a}_{y}+\hat{a}^{\dagger}_{y}\big{)},\qquad\hat{P}_{y}=-i\frac{\hbar\beta}{\sqrt{2}}\big{(}\hat{a}_{y}-\hat{a}^{\dagger}_{y}\big{)}.\end{split}$$

$$\hat{L}_{z}=i\hbar\big{(}\hat{a}_{x}\hat{a}^{\dagger}_{y}-\hat{a}^{\dagger}_{x}\hat{a}_{y}\big{)}.$$

$$\hat{H}^{\prime}_{xy}=\hbar\omega\big{(}\hat{a}^{\dagger}_{x}\hat{a}_{x}+\hat{a}^{\dagger}_{y}\hat{a}_{y}+1\big{)}.$$

$$\begin{split}\big{[}\hat{H}^{\prime}_{xy},\hat{L}_{z}\big{]}&=\!\bigg{[}\hbar\omega\big{(}\!\hat{a}^{\dagger}_{x}\hat{a}_{x}\!+\!\hat{a}^{\dagger}_{y}\hat{a}_{y}\!+\!1\!\big{)},i\hbar\big{(}\!\hat{a}_{x}\hat{a}^{\dagger}_{y}\!-\!\hat{a}^{\dagger}_{x}\hat{a}_{y}\!\big{)}\bigg{]},\\ &=\!i\hbar^{2}\omega\big{(}\!-\!\hat{a}_{x}\hat{a}^{\dagger}_{y}\!+\!\hat{a}^{\dagger}_{x}\hat{a}_{y}\!+\!\hat{a}^{\dagger}_{y}\hat{a}_{x}\!-\!\hat{a}_{y}\hat{a}^{\dagger}_{x}),\\ &=\!0.\end{split}$$

$$\hat{a}_{R}=\frac{1}{\sqrt{2}}\big{(}\hat{a}_{x}\!-\!i\hat{a}_{y}\big{)},\qquad\hat{a}_{L}=\frac{1}{\sqrt{2}}\big{(}\hat{a}_{x}\!+\!i\hat{a}_{y}\big{)}.$$

$$\big{[}\hat{a}_{R},\hat{a}^{\dagger}_{R}\big{]}=\big{[}\hat{a}_{L},\hat{a}^{\dagger}_{L}\big{]}=1.$$

$$\begin{split}\hat{a}^{\dagger}_{R}\hat{a}_{R}=\frac{1}{2}\big{(}\hat{a}^{\dagger}_{x}\hat{a}_{x}+\hat{a}^{\dagger}_{y}\hat{a}_{y}+i(\hat{a}_{x}\hat{a}^{\dagger}_{y}-\hat{a}^{\dagger}_{x}\hat{a}_{y})\big{)},\\ \hat{a}^{\dagger}_{L}\hat{a}_{L}=\frac{1}{2}\big{(}\hat{a}^{\dagger}_{x}\hat{a}_{x}+\hat{a}^{\dagger}_{y}\hat{a}_{y}-i(\hat{a}_{x}\hat{a}^{\dagger}_{y}-\hat{a}^{\dagger}_{x}\hat{a}_{y})\big{)}.\end{split}$$

$$\hat{H}^{\prime}_{xy}=\hbar\omega\big{(}\hat{N}_{R}+\hat{N}_{L}+1\big{)},\qquad\hat{L}_{z}=\hbar\big{(}\hat{N}_{R}-\hat{N}_{L}\big{)},$$

$$F_{n_R,n_L}=\frac{1}{\sqrt{(n_R)!(n_L)!}}({\hat{a}}_R^{†}{)}^{n_R}({\hat{a}}_L^{†}{)}^{n_L}{F}_{0,0}.$$

$$n=n_R+n_L,m_l=n_R-n_L.$$

$$\begin{array}[]{ll}n_{R}=n&;\;n_{L}=0\\ n_{R}=n-1&;\;n_{L}=1,\\ \qquad\vdots&\\ n_{R}=0&;\;n_{L}=n.\\ \end{array}$$

$$\begin{split}n=0&\quad\Rightarrow\quad\big{\{}\;\;n_{R}=0\;,\;n_{L}=0\;\rightarrow\;m_{l}=0\\ n=1&\quad\Rightarrow\quad\left\{\begin{array}[]{l}n_{R}=1\;,\;n_{L}=0\;\rightarrow\;m_{l}=+1\\ n_{R}=0\;,\;n_{L}=1\;\rightarrow\;m_{l}=-1\end{array}\right..\end{split}$$

$$m_l=n,n-2,n-4,\cdots ,-n+2,-n.$$

$$F_{\leavevmode\resizebox{140.92792pt}{}{$\displaystyle n_{R}\!=\!\frac{n\!+\!m_{l}}{2},n_{L}\!=\!\frac{n\!-\!m_{l}}{2}$}}.$$

Full text

I Introduction

In this article we present the detailed calculations required for a didactic comparison between the dynamics of the Schrödinger and Dirac equations for an electron in a uniform magnetic field. This is done in order to juxtapose these two quantum mechanical descriptions of the same system as a pedagogic tool as this comparison is not found in any quantum mechanics textbook. For that reason, this analysis would be an invaluable resource to instructors, who could make use this example, either for an advanced undergraduate quantum mechanics course or for a beginning graduate course.

It is well known that one of the simplest phenomenon at the quantum scale is the interaction between an electron (which is a spin one half fermion) and an external uniform axial magnetic field Peter_Rowe ; Tannoudji ; Haugset ; Yoshioka ; Thaller ; Greenshields . Although this physical system comprises only a single particle, whether it is at rest or in motion, its description is valuable as it furnishes the basic concepts required to understand other phenomena. For specific examples see the electron vortex beams Babiker ; Chowdhury ; Bliokh and the interaction of solid state materials with magnetism, known as the integer Quantum Hall effect Landau ; Feldman_Kahn ; Pekka ; laughlin ; Schattschneider . Additionally, this physical system shows an isomorphism to quantum optics, demonstrating a narrow relation to describe the Gaussian beam profile of electromagnetic radiation, leading to orthogonal states known as Orbital Angular Momentum of light Ole . This optical phenomenon has been obtained experimentally, due to the high accuracy now achievable and its great applicability Allen of simulations in innumerable quantum systems of interest, for example, the orbital angular momentum Hall effectZhang , which is an analogue to the integer quantum Hall effect.

A brief outline of the paper is as follows, in section II we describe the interaction between spin $1/2$ and a uniform magnetic field. In section III we explain in detail the Schrödinger equation of the electron in motion under the interaction of the external magnetic field, performing the analysis by using the ladder operator method. Consequently, we can notice that the wave function of the electron exhibit a cylindrical symmetry, this feature is related with a uniform external magnetic field. Also, we explain briefly the equivalence between the electron eigenstates and the Laguerre-Gaussian modes, which are characteristic of the orbital angular momentum of light in cylindrical coordinates. We then further present a short discussion about the electron energy levels and the Landau levels. Finally, in section IV we present a relativistic analysis using the Dirac equation for the electron wave function and discuss the main differences to the Schrödinger equation description.

II Spin-magnetic field interaction

Here we consider the quantum description of an electron initially at rest in a uniform external magnetic field. We consider the electron magnetic dipole moment \vec{\mu}\!=\!-\leavevmode\resizebox{35.23018pt}{}{\left(!\displaystyle\frac{e}{m_{0}}!\right)}\vec{\hat{S}}, where ($e\!=\!1,6021\!\!\times\!\!10^{-19}$C) is the modulus of the fundamental electric charge and the vectorial spin operator $\vec{\hat{S}}\!=\!\frac{\hbar}{2}\vec{\hat{\sigma}}$, where $\vec{\hat{\sigma}}\!=\!(\hat{\sigma}_{x},\hat{\sigma}_{y},\hat{\sigma}_{z})$ is the vector of Pauli matrices. Neglecting any external force, the Hamiltonian of the electron is an interaction potential of the magnetic dipole moment with the magnetic field is $\hat{H}\!=\!-\vec{\mu}\!\cdot\!\vec{B}$. To simplify the calculations we assume the orientation of the magnetic field is along the $\hat{z}$ axis, accordingly the electrons Hamiltonian is given by,

[TABLE]

The energy levels are $E_{\pm}\!=\!\pm\hbar\omega/2$, where $\omega\!=\!\frac{eB}{m_{0}}$ is the precession frequency of the magnetic dipole momentum around the external field, known as Larmor Frequency. In addition, the eigenstate of positive energy is $|0\rangle\!\!=$

\left(\!\begin{array}[]{c}1\\ 0\end{array}\!\right)

and the negative energy is $|1\rangle\!\!=$

\left(\!\begin{array}[]{c}0\\ 1\end{array}\!\right)

. These eigenstate are related with the spin orientation when it is aligned or anti-aligned to the magnetic field, respectively.

III Schrödinger equation of an electron in magnetic field

In contrast to last section, here we make a quantum description of an electron in motion, which is in a region with a uniform external magnetic field $\vec{B}$; in the classical description, the electron follows a helical path along the magnetic field showing an axial symmetry. According to this spatial feature, we apply a gauge transformation over the magnetic vector potential to take into account this symmetry. In particular, this specific symmetric gauge transformation is well know as the Landau gauge,

[TABLE]

The linear momentum operator is $\hat{p}\!\rightarrow\!\hat{p}-q\vec{A}$, where the electron charge $q\!=\!-e$. Then, in terms of momentum, the Hamiltonian is,

[TABLE]

The second term of the Hamiltonian can be expanded,

[TABLE]

where $\theta_{Br}$ is the angle between the magnetic field and the electron vector position. The third term in Eq. (3) will be calculated using,

[TABLE]

The first term in last expression is equal to,

[TABLE]

This result shows that the Landau gauge satisfies the Coulomb gauge. Continuing the calculation of the second term in Eq. (3) in a similar way,

[TABLE]

Representing the differential operator in term of the quantum linear momentum operator it follows that,

[TABLE]

Eliminating $\phi$ from the above results allows us to write the Hamiltonian (3) as,

[TABLE]

For simplicity and without loss of generality, we are going to assume a uniform magnetic field aligned along $\hat{z}$ axis. Consequently, the Hamiltonian is can be written,

[TABLE]

According to this result, we notice clearly the contribution of various terms of different natures in the Hamiltonian. These are, a free particle along $\hat{z}$ axis, a two-dimensional harmonic oscillator in $x\!-\!y$ plane with its projection of angular momentum along $\hat{z}$ axis and, finally, the electron spin interaction with the magnetic field. We can express $\psi(\vec{r},\sigma)$, the associated electron wave function, as the product of functions in according to their nature,

[TABLE]

where the spinorial function is,

[TABLE]

and $m_{s}$ is the quantum number equal to the spin projection along the magnetic field direction. In virtue of this, we can analyze the two-dimensional harmonic movement of the electron, getting the Schrödinger equation of a quantum harmonic oscillator,

[TABLE]

where $\omega\!=\!\frac{eB}{2m_{0}}$ is the harmonic frequency. The quantum harmonic oscillator is,

[TABLE]

We note that each energy level $E^{\prime}_{n}$, is related to ($n+1$) degenerate eigenstates $F_{n_{x},n_{y}}$,

[TABLE]

In relation to this degeneracy, we look for a physical quantity which enables a better characterization of states belonging to a certain energy level. Given the spatial symmetry of the harmonic potential under a rotation around the $\hat{z}$ axis and the axial component of the orbital angular momentum of the electron $\hat{L}_{z}\!=\!\hat{x}\hat{p}_{y}\!-\!\hat{y}\hat{p}_{x}$, we can analyze the two-dimensional motion of the electron by applying the ladder operator method to each Cartesian coordinate,

[TABLE]

where the constant $\beta\!\!=\!\!\sqrt{\frac{m_{0}\omega}{\hbar}}$ is reciprocal to the natural length of the harmonic oscillator Tannoudji . Moreover, these operators satisfy the following commutation relations $[\hat{a}_{x},\hat{a}^{\dagger}_{x}]\!=\![\hat{a}_{y},\hat{a}^{\dagger}_{y}]\!=\!1$. Furthermore, we can express $\hat{x}$, $\hat{y}$, $\hat{p}_{x}$, $\hat{p}_{y}$, as function of the ladder operators. In that way, we have,

[TABLE]

The angular momentum operator is given by,

[TABLE]

Similarly, we can express the harmonic Hamiltonian as,

[TABLE]

Calculating the commutation relation between the angular momentum and the Hamiltonian yields:

[TABLE]

Thus, we can show that $\hat{L}_{z}$ is a constant of motion. Consequently, the angular momentum operator has the same set of eigenstates as $\hat{H}^{\prime}_{xy}$; in this way, we can expose an axial symmetry in the system. As such, we are going to implement the cylindrical coordinate system to describe the quantum operators corresponding to the electron position and momentum. Additionally, we will borrow from geometrical optics the notion of the Jones’ vectors for the circular polarization of light Jones_vector . One can introduce new ladder operators related to the axial symmetry $\hat{a}_{R}$ and $\hat{a}_{L}$, which are right and left circular operators respectively. These are denoted as a function of the Cartesian ladder operators as,

[TABLE]

Particularly, when these operators $\hat{a}_{R}$ and $\hat{a}_{L}$ act on the eigenstate $F_{n_{x},n_{y}}$, they generate a quantum state which is a linear combination of $F_{n_{x}-1,n_{y}}$ and $F_{n_{x},n_{y}-1}$. This results in a state with energy smaller by $\hbar\omega$. Contrastingly, $\hat{a}^{\dagger}_{R}$ and $\hat{a}^{\dagger}_{L}$ generate a quantum state with energy bigger by $\hbar\omega$.

Moreover, these operators satisfy the following commutation relations,

[TABLE]

Hence, we can express,

[TABLE]

It leads us to express $\hat{H}^{\prime}_{xy}$ and $\hat{L}_{z}$ as,

[TABLE]

where $\hat{N}_{R}\!=\!\hat{a}^{\dagger}_{R}\hat{a}_{R}$ and $\hat{N}_{L}\!=\!\hat{a}^{\dagger}_{L}\hat{a}_{L}$ are known as the Hermitian number operator. When acting on the eigenstate $F_{n_{x},n_{y}}$, these operators simply return a positive integer $n_{R}$ and $n_{L}$, respectively. The physical meaning of these numbers are related to the current energy level, as we will show later.

If we denote the ground state as $F_{0,0}$, as the state of the lowest energy. In relation to this state, we can express any other eigenstate as a consecutive action of the operators $\hat{a}^{\dagger}_{R}$ and $\hat{a}^{\dagger}_{L}$. Thus,

[TABLE]

Since $F_{n_{R},n_{L}}$ is a set of eigenstates in common to $\hat{H}^{\prime}_{xy}$ and $\hat{L}_{z}$, we can obtain the energy levels $\hbar\omega(n\!+\!1)$ and the eigenvalue $m_{l}\hbar$. In this way, we can define the main quantum number $n$ and the orbital quantum number $m_{l}$ as function of $n_{R}$ and $n_{L}$,

[TABLE]

Since $n_{R}$ and $n_{L}$ are positive integers, there are ($n\!+\!1$) degenerate eigenstate for every energy level, as shown in Eq. (15). This proves the next relation,

[TABLE]

Observing the particular case $n\!=\!0,1$ and its possible values of $m_{l}$, we have,

[TABLE]

In general, for every fixed value of $n$, there is a set of values of $m_{l}$,

[TABLE]

In virtue of this transformation of notation, we can express the eigenstate $F_{n_{R},n_{L}}$ by using the pair of quantum numbers $n,m_{l}$, which are more appropriate to characterize this quantum system due to the observable $\hat{L}_{z}$. Thus, we have,

[TABLE]

Then, we use the cylindrical coordinates,

[TABLE]

We can rewrite the Hamiltonian of the two-dimensional harmonic oscillator in this coordinate system. Also, we can transform the set of operator $\hat{a}_{R}$ and $\hat{a}_{L}$. In that way, we obtain,

[TABLE]

The radial distance and the azimuthal angle are defined as $\varrho\!=\!\sqrt{x^{2}\!+\!y^{2}}$ and $\varphi\!=\!\textnormal{arctan}\left(y/x\right)$. We can verify that,

[TABLE]

and the Cartesian differentiation will be written as,

[TABLE]

Then,

[TABLE]

Thus, the operators $\hat{a}_{R}$ and $\hat{a}^{\dagger}_{R}$ are expressed in cylindrical coordinate as,

[TABLE]

Analogously, we can obtain,

[TABLE]

In particular, the action of $\hat{a}_{R}$ or $\hat{a}_{L}$ on the ground state $F_{n_{R}=0,n_{L}=0}$ is,

[TABLE]

This result means physically the impossibility of annihilating a quantum of energy $\hbar\omega$ from the ground state. In virtue of this fact, we can solve the differential equation Eq. (38) and obtain the normalized eigenfunction,

[TABLE]

In contrast, when $\hat{a}^{\dagger}_{R}$ acting $n_{R}$-times on $F_{0,0}(\varrho,\varphi)$ yields the creation of $n_{R}$ quanta of energy $\hbar\omega$ on the ground state,

[TABLE]

Similarly, the action of $\hat{a}^{\dagger}_{L}$ $n_{L}$-times on $F_{0,0}(\varrho,\varphi)$ yields,

[TABLE]

After obtaining these results, it’s worth mentioning that every energy level $E_{n}\!=\!\hbar\omega(n\!+\!1)$ corresponds to the maximum projection of the angular momentum $m_{l}\!=\!n$ in Eq. (40), or the minimum value $m_{l}\!=\!-n$ in Eq. (41). Nonetheless, by a consecutive action of the operator $\hat{a}^{\dagger}_{L}$ $n_{L}$-times on the eigenfunction $F_{n_{R},0}$ or $\hat{a}^{\dagger}_{R}$ $n_{R}$-times on the eigenfunction $F_{0,n_{L}}$, we obtain the eigenfunction $F_{n_{R},n_{L}}(\varrho,\varphi)$ for any pair of quantum numbers $n_{R}$, $n_{L}$. Naturally, after calculating a lot of eigenfunctions by the action of the two creation operators, we start to obtain eigenfunctions which are proportional to the generalized Laguerre polynomials, denoted by $L^{(\alpha)}_{n}(x)$, multiplied by a Gaussian function which limits spatially the probability density of the electron in transversal directionTannoudji ; Yoshioka . In general, any wave function associated with the eigenstate ($n\!=\!n_{R}\!+\!n_{L}$) with angular momentum ($m_{l}\!=\!n_{R}\!-\!n_{L}$) can be written as follows,

[TABLE]

where $\beta^{2}\varrho^{2}$ is the argument of the polynomial function. The normalization constant is,

[TABLE]

The generalized Laguerre polynomials can be explicitly written using the formula,

[TABLE]

Exemplifying a particular case in Figure 1, we show the probability density of the eigenstate with null angular momentum $|F_{4,0}(\varrho,\varphi)|^{2}$. Additionally, in Table 1 we show the electron radial eigenfunction to the first six values of $n$. In relation to these eigenfunctions, we obtain more information about every electron eigenstate when we calculate its radial probability density as,

[TABLE]

Thus, in Figure 2 we show the radial probability densities of every eigenstates belonging to the first six values of $n$.

We want to briefly mention a close relationship between the radial eigenfunction of the electron and Quantum Optics. These eigenfunctions are the same as those we can obtain as a natural solutions of the wave equation under the paraxial approximation Goubau ; Andrews , for a Gaussian beam of light in cylindrical coordinates. Commonly, in quantum optics these functions are called Laguerre-Gaussian modes, which are expressed through the generalized Laguerre polynomials $L^{|l|}_{p}$, where $p\!\geq\!0$ is the radial index which is related to the number of rings in the probability density of every eigenfunction; we can notice this number in Figure 2 along every diagonal from lower left to upper right. The integer $l$ which is the azimuthal index, has the same physical meaning as the axial electron angular momentum.

This physical feature lets us establish an isomorphic behavior between the quantum mechanics description of an electron in a uniform magnetic field and the Laguerre-Gaussian modes of lightOle . Where the amplitude of the paraxial beam of light corresponds to the spatial probability density of the electron and the Gouy phase of opticsSimin_Feng assumes the meaning of time for the quantum harmonic oscillator using an unitary dynamical evolution. Recently, these Laguerre-Gaussian beams have been of considerable practical interest, particularly in the field of optical trapping, where they are applied to study the driving of micromachined elements with light as optical tweezersPaula .

We can obtain the Schrödinger equation in cylindrical coordinates by using the full definition of the electron wave function,

[TABLE]

In this way, we obtain the total energy of the electron in relation to its eigenstates,

[TABLE]

In general, the electron energy levels are $E\!\geq\!0$, being the zero of energy associated to eigenstates with a null linear momentum $p_{z}$ and the lowest possible projection of angular momentum $m_{l}\!=\!-n\!$ in relation to the main quantum number $n$, also, a spin orientation opposite to the magnetic field $m_{s}\!=\!-1/2$. Additionally, under some of the last conditions and getting free $m_{s}$, we can obtain the energy levels associated to the spin interaction with the external magnetic field as we shown in section II, but in this case we see the contribution of the zero-point energy which is related to the harmonic potential and that is the reason to have this energy shift.

[TABLE]

Furthermore, we notice that the addition of $n$ and $m_{l}$ is always an even number. Then, we can rewrite the energy expression using a new integer number,

[TABLE]

where $r$ is associated to every energy level and known as the Landau level. Also, we can see easily that eigenstates with the same spatial probability density belong to different energy levels. For instance, $F_{1,1}$ and $F_{1,-1}$ are associated to different energies, but their spatial probability density are equals. In general we have,

[TABLE]

In relation to Eq. (49), we can see that every Landau Level is highly degenerate, because the electron eigenstates with different quantum numbers $n$ and $m_{l}$ are related to the same energy level, as we see in Figure 3. This feature can be observed in the quantization of the cyclotron orbits of charged particles in magnetic field Mikhailov . Such charged particles can only occupy orbits with discrete energy values, but these levels are degenerate, where the degeneracy is associated to the number of electrons per level, which is directly proportional to the strength of the applied magnetic field.

Moreover, this physical aspect can be experimentally appreciated within solid materials, where there are electronic oscillations under the action of an external magnetic field; when we apply a differential of electric potential through the material, that leads to discrete values of electric current directly related to the Landau levels of the electrons. This phenomenon is commonly called Integer Quantum Hall Effect Jeanneret , and its useful applications have been evidenced in quantum metrology in order to acquire more information about microscopic details of semiconductors. Also, evidence of Landau levels has been obtained in the propagation of electron vortex beams along an external longitudinal magnetic field Franco .

IV Dirac equation for an electron in magnetic field

In this section we present the relativistic description via Dirac equation for an electron in a region with a uniform magnetic field $\vec{B}$. According to this quantum description, the electron eigenstate is expressed by the spinorial formalism. In general, a Dirac spinor for the electron is a column vector with four-elements, each one of these are related to the eigenfunction (42) obtained from the relevant Schrödinger equation. Also, the order of the spinor’s elements are associated with the spin orientation and the energy sign of the electron Strange ; Greiner ; Bransden , where the energy negative elements are associated to the electron’s antiparticle, commonly called positron.

After this short description, the electron Dirac equation is,

[TABLE]

where $U$ is a general electron spinor. We can write the linear momentum operator according to the canonical transformation $\hat{\vec{p}}\!\rightarrow\!\hat{\vec{p}}+e\vec{A}$, where $\vec{A}$ is the magnetic vector potential. Additionally, we have $\tilde{\alpha}$ and $\tilde{\beta}$ which are four Hermitian matrices with dimension $4\!\times\!4$, satisfying the condition $\tilde{\alpha}^{2}\!=\!\tilde{\beta}^{2}\!=\!\hat{\mathbb{1}}$,

[TABLE]

Therefore, we obtain the time-independent Dirac Hamiltonian $\hat{H}\!=\!c\tilde{\alpha}\cdot\hat{\vec{p}}+\tilde{\beta}m_{0}c^{2}$. We can express the electron’s spinor evolution through an unitary transformation,

[TABLE]

where $\phi$ and $\chi$ are two scalar functions that represent the first two elements of positive energy and the last two elements of negative energy of the electron’s spinor respectively. When we substitute the spinor (53) in Eq. (51), and using the matrix notation, we obtain,

[TABLE]

From the first and second row of the matrix we can obtain two linear equations equal to the electron total energy times the respective scalar functions defined before,

[TABLE]

When we isolate $\phi$ from (55a) and $\chi$ from (55b) we find both,

[TABLE]

Replacing $\chi$ by $\phi$ and vice versa, results in,

[TABLE]

The first term on the right is equal to

[TABLE]

Similarly, the second term on the right side of Eq. (57) is equal to

[TABLE]

Then we obtain,

[TABLE]

According to the last results, we rewrite Eq. (57), getting,

[TABLE]

It’s worth mentioning that $e\hbar\hat{\sigma}\cdot\vec{B}$ is the term associated with the interaction potential between the spin of the electron and the external magnetic field. Usually, this term is put in by hand as we, see Eq. (3). This is done in order to have a good quantum description by Schrödinger equation of an electron with spin $1/2$ in a region with an external magnetic field. On the other hand, we have that this interaction potential appears naturally by using the Dirac equation to describe the same physical system.

Continuing the analysis, we make a simplification (without loss of physical generality), assuming an external uniform magnetic field along $\hat{z}$ axis. In this way, we can express the magnetic vector potential by the Landau gauge which lets us to put in evidence the axial symmetry as we explained in last section,

[TABLE]

Replacing the vector potential in Eq. (61) and dividing this equation by $2m_{0}$, we have,

[TABLE]

The first three terms in the right hand side correspond to the Hamiltonian of the two-dimensional harmonic oscillator and the magnetic field interaction between angular momentum and spin. The fourth term corresponds to the Hamiltonian of a free particle along $\hat{z}$ axis. We can express the first two elements of the electron spinor by the product of the eigenfunctions of the two-dimensional harmonic oscillator obtained in previous section, and a plane wave function of a free particle along the direction of the external field,

[TABLE]

where $\Gamma$ is the spinorial function that we defined previously in Eq. (12). In relation to the first two spinor elements, which are related to the positive energy particle, we have the spin orientation up or down. We further have the last two spinor elements which are related with negative energy particle (i.e., an antiparticle); these are related to the spin orientation as before. Therefore, we can denote the two positive spinor elements as,

[TABLE]

When the spin operator acts on the spinor positive energy elements, we obtain the relation,

[TABLE]

Additionally, the commutation relations between the axial angular momentum and the harmonic Hamiltonian are equal to zero, we can say that $\hat{L}_{z}$ is a conserved physical quantity, and its eigenfunctions are the eigenfunctions of the two-dimensional harmonic oscillator. Thus, we get the following relation,

[TABLE]

where $m_{l}$ is the quantum number of the axial angular momentum. When we replace Eq. (64) in Eq. (63) we obtain,

[TABLE]

Rearranging the last equation,

[TABLE]

we notice that the left hand terms in brackets correspond to the Hamiltonian of the two-dimensional harmonic oscillator in Cartesian coordinates with a characteristic frequency $\omega\!=\!eB/2m_{0}$. Therefore, we can say that the right hand terms is equal to the energy of the oscillator. We accordingly have,

[TABLE]

Due to the total energy is to second power, when we isolate $E$, we obtain two signs for the total energy of the particle, being the positive sign related to the electron and the negative sign to the positron,

[TABLE]

If we have the particular case of an electron with null linear momentum, an angular momentum projection equal to $m_{l}\!=\!-n$ and a weak magnetic field in relation to the rest mass of the electron, then we can approximate the energy expression as,

[TABLE]

In relation to this result we obtain the energy levels of the spin interaction with an external magnetic field as we shown in the last two sections, also the energy shift related to the zero-point energy of the harmonic potential, and the rest mass of the electron, which is relevant by a relativistic description as we done using Dirac equation.

According to the axial symmetry exhibited by the two dimensional harmonic Hamiltonian Eq. (13), we can make a coordinate system transformation in order to rewrite the Hamiltonian in cylindrical coordinates, as we have done in previous section. We thus obtain the same radial eigenfunctions of the electron, Eq. (42) in the case where the electron state is expressed by the Dirac spinor.

Specifically, if we describe the electron using the Dirac equation, we have to consider the first two elements of the Dirac spinor as the radial eigenfunctions obtained in last section. Thus, we implement the relation (56), according to the orientation of the electron’s spin. In the same way, we can obtain the other two elements of the spinor which are related to negative energy. Then, we express the complete Dirac spinor associated with a particular quantum state of the electron.

To exemplify this process, we calculate a particular case to obtain the Dirac spinor for an electron in a specific state. We choose the ground state of the electron $n\!=\!0$, $m_{l}\!=\!0$ and its spin orientation $m_{s}\!=\!+1/2$. Then, the electron relativistic energy is,

[TABLE]

Previously, we have obtained the radial wave function of the ground state via Schrödinger equation,

[TABLE]

We can write the first two elements of the spinor,

[TABLE]

when we transform the relations (56) in cylindrical coordinates.

Firstly, we have to express every unitary vectors of Cartesian coordinate as,

[TABLE]

In this way, the magnetic vector potential is,

[TABLE]

Similarly, the nabla operator in cylindrical coordinates is,

[TABLE]

Also, the Pauli matrices vector is written as,

[TABLE]

Once we have calculated all the relevant transformations, we are able to use Eq. (56) to calculate the two negative energy elements of the spinor as follows,

[TABLE]

where $\hat{a}^{\dagger}_{R}$ and $\hat{a}_{R}$ are the ladder operators defined in Eq. (37). Therefore, the two electron spinor elements corresponding to the negative energy of the ground state can be expressed as,

[TABLE]

Therefore, we obtain the complete Dirac spinor associated to the ground state of the electron,

[TABLE]

where its normalization constant is,

[TABLE]

We would like to emphasize the relevance of the two negative elements of the Dirac spinor for an electron in two specific cases: i) When the electron has a great linear momentum and it is the dominant physical quantity, we have that the electron speed is very close to the speed of light and the third spinor element becomes relevant for this description. ii) When we have a strong confinement of the electron and it is trapped in a region less than or equal to its Compton wavelength Leary ; Velasco . Thereby, the fourth spinor element becomes relevant to the physical description. Thus, we can obtain this last case when the electron is in a region with a very strong magnetic field, yielding a high characteristic frequency $\omega$ of the harmonic potential, as consequence a very small natural length of the harmonic oscillator $\sqrt{\hbar/m_{0}\omega}$. As such, we show in Figure 4 a different radial probability density of the electron ground state via Dirac equation, in comparison to that obtained before via Schrödinger equation.

The fourth element of the spinor shows an excited eigenstate $n,m_{l}\!=\!1,1$; whose contribution is comparable to the first element. Subsequently, this element of the spinor becomes relevant for the description of the electron quantum state; and this is why we observe a more open electron radial probability density via the Dirac equation.

We advise to the reader to see a short animation animation about the radial probability density according to both theories in relation to the external magnetic field strength.

In summary, we can conclude that the quantum mechanical description via the Dirac equation requires more complicated mathematical steps, but that it also brings new concepts, such as the negative energy of the electron and a different treatment through the application of Dirac spinors; this yields a more general way to describe the quantum state of an electron when relativistic effects are taken into account which is richer than the Schrödinger description.

V Acknowledgements

The first author wish to thank to his doctoral advisor Eduardo, for all fruitful discussions and his great listening and grateful understanding. Also to Benjamin for his great collaboration and finally the Brazilian agency CAPES for its financial support.