Pauli’s Electron in Ehrenfest and Bohm Theories, a Comparative Study

1. Introduction

Quantum mechanics is usually interpreted by the Copenhagen approach. This approach objects to the physical reality of the quantum wave function and declares it to be epistemological (a tool for estimating probability of measurements) in accordance with the Kantian [1] depiction of reality, and its denial of the human ability to grasp any thing in its reality (ontology). However, we also see the development of another approach of prominent scholars that think about quantum mechanics differently. This school believes in the ontological existence of the wave function. According to this approach the wave function is an element of reality much like an electromagnetic field. This was supported by Einstein and Bohm [2,3,4] has resulted in different understandings of quantum mechanics among them the fluid realization championed by Madelung [5,6] which stated that the modulus square of the wave function is a fluid density and the phase is a potential of the velocity field of the fluid.

A non relativistic quantum equation for a spinor was first introduced by Wolfgang Pauli in 1927 [7], this was motivated by the need to explain the Stern-Gerlach experiments. Later it was shown that the Pauli equation is a low velocity limit of the relativistic Dirac equation. This equation is based on a two dimensional operator matrix Hamiltonian. Two dimensional operator matrix Hamiltonians are common in the literature [8,9,10,11,12,13,14,15,16,17,18,19,20,21] and describe many types of quantum systems. A Bohmian analysis of the Pauli equation was given by Holland and others [3,4], however, the analogy of the Pauli theory to fluid dynamics and the notion of spin vorticity were not considered. In [23] spin fluid dynamics was introduces for a single electron with a spin. One thus must contemplate where do those internal energies originate? The answer to this question seems to come from measurement theory [25,27]. Fisher information is a basic notion of measurement theory, and is a figure of merit of a measurement quality of any quantity. It was shown [27] that this notion is the internal energy of a spin less electron (up to a proportionality constant) and can be used to partially interpret the internal energy of an electron with spin. An attempt to derive most physical theories from Fisher information is due to Frieden [28]. It was suggested [29] that there exist a velocity field such that the Fisher information will given a complete explanation for the spin fluid internal energy. It was also suggested that one may define comoving scalar fields as in ideal fluid mechanics, however, this was only demonstrated implicitly but not explicitly. A common feature of previous work on the fluid & Fisher information interpretation of quantum mechanics, is the negligence of electromagnetic interaction thus setting the vector potential to zero, this was recently corrected in [30].

Ehrenfest [22] published his paper in 1927 as well with the title: "Remark on the approximate validity of classical mechanics within quantum mechanics". Using this approach we can accept the orthodox Copenhagen’s interpretation denying a trajectory of the electron but at the same time accept the existence of a trajectory of the electron’s position vector expectation value through Ehrenfest theorem. The Ehrenfest approach is thus independent of interpretation, and can be applied according to both the Copenhagen and Bohm schools. However, only in the Bohm approach one may compare the trajectory of the electron to that of its expectation value.

We will begin this paper by reminding the reader the basic equation describing the motion of a classical electron. This will be followed by a discussion of Schrödinger equation with a non trivial vector potential and its interpretation in terms of Bohmian equation of motion with a quantum force correction. Then we introduce Pauli’s equation with a vector potential and interpret it in terms a Bohmian equation of motion with a quantum force correction which is different from the Schrödinger case . Finally we derive an equation for Pauli’s electron position vector expectation value using Ehrenfest theorem and compare the result to the results obtained in Bohm’s approach, similarities and differences will arise, a concluding section will follow discussing the Stern-Gerlach experiment.

2. A Classical Charged Particle

Consider a classical particle with the coordinates $\overrightarrow{x}(t)$, mass m and charge e interacting with a given electromagnetic vector potential $\overrightarrow{A}(\overrightarrow{x},t)$ and scalar potential $\phi (\overrightarrow{x},t)$. We will not be interested in the effects of the particle on the field and thus consider the field as “external”. The action of said particle is:

$$\begin{array}{ccc}\hfill \mathcal{A}& =& {\int }_{t1}^{t2}Ldt,L=L_0+L_i\hfill \\ \hfill L_0& \equiv & \frac{1}{2}m{v}^2,L_i\equiv e(\overrightarrow{A}·\overrightarrow{v}-\phi ),\overrightarrow{v}\equiv \frac{d\overrightarrow{x}}{dt}\equiv \dot{\overrightarrow{x}},v=\left|\overrightarrow{v}\right|.\hfill \end{array}$$

(1)

The variation of the two parts of the Lagrangian are given by:

$$\delta L_0=m\dot{\overrightarrow{x}}·\delta \dot{\overrightarrow{x}}=\frac{d(m\overrightarrow{v}·\delta \overrightarrow{x})}{dt}-m\dot{\overrightarrow{v}}·\delta \overrightarrow{x}$$

(2)

$$\begin{array}{ccc}\hfill \delta L_i& =& e\left(\delta \overrightarrow{A}·\overrightarrow{v}+\overrightarrow{A}·\delta \dot{\overrightarrow{x}}-\delta \phi \right)\hfill \\ & =& e\left({\partial }_k\overrightarrow{A}·\overrightarrow{v}\delta x_k+\frac{d(\overrightarrow{A}·\delta \overrightarrow{x})}{dt}-\delta \overrightarrow{x}·\frac{d\overrightarrow{A}}{dt}-\overrightarrow{\nabla }\phi ·\delta \overrightarrow{x}\right),\hfill \end{array}$$

(3)

in the above ${\partial }_k\equiv \frac{\partial }{\partial x_k}$ and $\overrightarrow{\nabla }\equiv (\frac{\partial }{\partial x},\frac{\partial }{\partial y},\frac{\partial }{\partial z})\equiv (\frac{\partial }{\partial x_1},\frac{\partial }{\partial x_2},\frac{\partial }{\partial x_3})$. We use the Einstein summation convention in which a Latin index (say $k,l$) takes one of the values $k,l\in [1,2,3]$. We may write the total time derivative of $\overrightarrow{A}$ as:

$$\frac{d\overrightarrow{A}(\overrightarrow{x}(t),t)}{dt}={\partial }_t\overrightarrow{A}+v_l{\partial }_l\overrightarrow{A},{\partial }_t\equiv \frac{\partial }{\partial t}.$$

(4)

Thus, the variation $\delta L_i$ can be written in the following form:

$$\delta L_i=\frac{d(e\overrightarrow{A}·\delta \overrightarrow{x})}{dt}+e\left[\left({\partial }_k{A}_l-{\partial }_l{A}_k\right)v_l-{\partial }_t{A}_k-{\partial }_k\phi \right]\delta x_k,$$

(5)

Defining the electric and magnetic fields in the standard way:

$$\overrightarrow{B}\equiv \overrightarrow{\nabla }\times \overrightarrow{A},\overrightarrow{E}\equiv -{\partial }_t\overrightarrow{A}-\overrightarrow{\nabla }\phi ,$$

(6)

it follows that:

$${ϵ}_{kln}{B}_n={\partial }_k{A}_l-{\partial }_l{A}_k,E_k=-{\partial }_t{A}_k-{\partial }_k\phi ,$$

(7)

in which ${ϵ}_{kln}$ is the three index antisymmetric tensor. Thus, we may write $\delta L_i$ as:

$$\delta L_i=\frac{d(e\overrightarrow{A}·\delta \overrightarrow{x})}{dt}+e\left[{ϵ}_{kln}{v}_l{B}_n+E_k\right]\delta x_k=\frac{d(e\overrightarrow{A}·\delta \overrightarrow{x})}{dt}+e\left[\overrightarrow{v}\times \overrightarrow{B}+\overrightarrow{E}\right]·\delta \overrightarrow{x}.$$

(8)

We use the standard definition of the Lorentz force (MKS units):

$${\overrightarrow{F}}_L\equiv e\left[\overrightarrow{v}\times \overrightarrow{B}+\overrightarrow{E}\right]$$

(9)

to write:

$$\delta L_i=\frac{d(e\overrightarrow{A}·\delta \overrightarrow{x})}{dt}+{\overrightarrow{F}}_L·\delta \overrightarrow{x}.$$

(10)

Combining the variation of $L_i$ given in equation (10) and the variation of $L_0$ given in equation (2), it follows from equation (1) that the variation of L is:

$$\delta L=\delta L_0+\delta L_i=\frac{d\left((m\overrightarrow{v}+e\overrightarrow{A})·\delta \overrightarrow{x}\right)}{dt}+(-m\dot{\overrightarrow{v}}+{\overrightarrow{F}}_L)·\delta \overrightarrow{x}.$$

(11)

Thus, the variation of the action is:

$$\delta A={\int }_{t1}^{t2}\delta Ldt={\left.(m\overrightarrow{v}+e\overrightarrow{A})·\delta \overrightarrow{x}\right|}_{t1}^{t2}-{\int }_{t1}^{t2}(m\dot{\overrightarrow{v}}-{\overrightarrow{F}}_L)·\delta \overrightarrow{x}dt.$$

(12)

Since the classical trajectory is such that the variation of the action on it vanishes for a small modification of the trajectory $\delta \overrightarrow{x}$ that vanishes at $t1$ and $t2$ but is otherwise arbitrary it follows that:

$$m\dot{\overrightarrow{v}}={\overrightarrow{F}}_L=e\left[\overrightarrow{v}\times \overrightarrow{B}+\overrightarrow{E}\right]\Rightarrow \dot{\overrightarrow{v}}=\frac{e}{m}\left[\overrightarrow{v}\times \overrightarrow{B}+\overrightarrow{E}\right].$$

(13)

Thus, the dynamics of a classical particle in a given electric and magnetic field is described by a single number, the ratio between its charge and mass:

$$k\equiv \frac{e}{m}\Rightarrow \dot{\overrightarrow{v}}=k\left[\overrightarrow{v}\times \overrightarrow{B}+\overrightarrow{E}\right].$$

(14)

The reader is reminded that the connection between the electromagnetic potentials and the fields is not unique. Indeed performing a gauge transformation to obtain a new set of potentials:

$${\overrightarrow{A}}^{\prime }=\overrightarrow{A}+\overrightarrow{\nabla }\Lambda ,{\phi }^{\prime }=\phi -{\partial }_t\Lambda .$$

(15)

we obtain the same fields:

$${\overrightarrow{B}}^{\prime }=\overrightarrow{\nabla }\times {\overrightarrow{A}}^{\prime }=\overrightarrow{\nabla }\times \overrightarrow{A}=\overrightarrow{B},{\overrightarrow{E}}^{\prime }=-{\partial }_t{\overrightarrow{A}}^{\prime }-\overrightarrow{\nabla }{\phi }^{\prime }=-{\partial }_t\overrightarrow{A}-\overrightarrow{\nabla }\phi =\overrightarrow{E}.$$

(16)

3. Schrödinger’s Theory

Quantum mechanics according to the Copenhagen interpretation has lost faith in our ability to predict precisely the whereabouts of even a single particle. What the theory does predict precisely is the evolution in time of a quantity denoted "the quantum wave function", which is related to a quantum particle whereabouts in a statistical manner. This evolution is described by an equation suggested by Schrödinger [24]:

$$i\hslash \dot{\psi }={\widehat{H}}_S\psi ,{\widehat{H}}_S=-\frac{1}{2m}{\left(\hslash \overrightarrow{\nabla }-ie\overrightarrow{A}\right)}^2+e\phi$$

(17)

in the above $i=\sqrt{-1}$ and $\psi$ is the complex wave function. $\dot{\psi }=\frac{\partial \psi }{\partial t}$ is the partial time derivative of the wave function. $\hslash =\frac{h}{2\pi }$ is Planck’s constant divided by $2\pi$ and m is the particles mass. However, this presentation of quantum mechanics is rather abstract and does not give any physical picture regarding the meaning of the quantities involved. Thus we write the quantum wave function using its modulus a and phase $\varphi$:

$$\psi =a{e}^{i\varphi }.$$

(18)

We define the velocity field:

$${\overrightarrow{v}}_S=\frac{\hslash }{m}\overrightarrow{\nabla }\varphi -\frac{e}{m}\overrightarrow{A}$$

(19)

and the mass density is defined as:

$$\widehat{\rho }=m{a}^2.$$

(20)

It is easy to show from equation (17) that the continuity equation is satisfied:

$$\frac{\partial \widehat{\rho }}{\partial t}+\overrightarrow{\nabla }·(\widehat{\rho }{\overrightarrow{v}}_S)=0$$

(21)

Hence ${\overrightarrow{v}}_S$ field is the velocity associated with mass conservation. However, it is also the mass associate with probability $a^2$ (by Born’s interpretational postulate) and charge density $\rho =e{a}^2$. The equation for the phase $\varphi$ derived from equation (17) is as follows:

$$\hslash \frac{\partial \varphi }{\partial t}+\frac{1}{2m}{\left(\hslash \overrightarrow{\nabla }\varphi -e\overrightarrow{A}\right)}^2+e\phi =\frac{{\hslash }^2{\nabla }^{2}a}{2ma}$$

(22)

In term of the velocity defined in equation (19) one obtains the following equation of motion (see Madelung [5] and Holland [3]):

$$\frac{d{\overrightarrow{v}}_S}{dt}=\frac{\partial {\overrightarrow{v}}_S}{\partial t}+({\overrightarrow{v}}_S·\overrightarrow{\nabla }){\overrightarrow{v}}_S=-\overrightarrow{\nabla }\frac{Q}{m}+k(\overrightarrow{E}+\overrightarrow{v}\times \overrightarrow{B})$$

(23)

The right hand side of the above equation contains the "quantum correction":

$$Q=-\frac{{\hslash }^2}{2m}\frac{{\overrightarrow{\nabla }}^2\sqrt{\widehat{\rho }}}{\sqrt{\widehat{\rho }}}.$$

(24)

For the meaning of this correction in terms of information theory see: [23,27,29]. These results illustrates the advantages of using the two variables, phase and modulus, to obtain equations of motion that have a substantially different form than the familiar Schrödinger equation (although having the same mathematical content) and have straightforward physical interpretations [2].

The quantum correction Q will of course disappear in the classical limit $\hslash \to 0$, but even if one intends to consider the quantum equation in its full rigor, one needs to take into account the expansion of an unconfined wave function. As Q is related to the typical gradient of the wave function amplitude it follows that as the function becomes smeared over time and the gradient becomes small the quantum correction becomes negligible. To put in quantitative terms:

$${\overrightarrow{F}}_Q=-\overrightarrow{\nabla }Q\simeq \frac{{\hslash }^2}{2m{L}_R^3},L_R\simeq \frac{R}{|\overrightarrow{\nabla }R|}$$

(25)

in which $L_R$ is the typical length of the amplitudes gradient. Thus:

$$|F_Q|<<|F_L|\Rightarrow L_R>>L_{Rc}={\left(\frac{{\hslash }^2}{2m{F}_L}\right)}^{\frac{1}{3}}.$$

(26)

in which ${\overrightarrow{F}}_L$ is the classical Lorentz force given in equation (9). For the current application in which an electron transverses many rotations of macroscopic size this terms in not important.

4. Pauli’s Theory

Schrödinger’s quantum mechanics is limited to the description of spin less particles. Indeed the need for spin became necessary as Schrödinger equation could not account for the result of the Stern Gerlach experiments, predicting a single spot instead of the two spots obtained for hydrogen atoms. Thus Pauli introduced his equation for a non-relativistic particle with spin is given by:

$$i\hslash \dot{\psi }={\widehat{H}}_P\psi ,{\widehat{H}}_P=-\frac{{\hslash }^2}{2m}{[\overrightarrow{\nabla }-\frac{ie}{\hslash }\overrightarrow{A}]}^2+\mu \overrightarrow{B}·\overrightarrow{\sigma }+e\phi ={\widehat{H}}_{S}I+\mu \overrightarrow{B}·\overrightarrow{\sigma }$$

(27)

$\psi$ here is a two dimensional complex column vector (also denoted as spinor), ${\widehat{H}}_P$ is a two dimensional hermitian operator matrix, $\mu$ is the magnetic moment of the particle, and I is a two dimensional unit matrix. $\overrightarrow{\sigma }$ is a vector of two dimensional Pauli matrices which can be represented as follows:

$${\sigma }_1=\left(\begin{array}{cc}0& 1\\ 1& 0\end{array}\right),{\sigma }_2=\left(\begin{array}{cc}0& -i\\ i& 0\end{array}\right),{\sigma }_3=\left(\begin{array}{cc}1& 0\\ 0& -1\end{array}\right).$$

(28)

The ad-hoc nature of this equation was later amended as it became clear that this is the non relativistic limit of the relativistic Dirac equation. A spinor $\psi$ satisfying equation (27) must also satisfy a continuity equation of the form:

$$\frac{\partial {\rho }_p}{\partial t}+\overrightarrow{\nabla }·\overrightarrow{j}=0.$$

(29)

In the above:

$${\rho }_p={\psi }^{†}\psi ,\overrightarrow{j}=\frac{\hslash }{2mi}[{\psi }^{†}\overrightarrow{\nabla }\psi -(\overrightarrow{\nabla }{\psi }^{†})\psi ]-k\overrightarrow{A}{\rho }_p.$$

(30)

The symbol ${\psi }^{†}$ represents a row spinor (the transpose) whose components are equal to the complex conjugate of the column spinor $\psi$. Comparing the standard continuity equation to equation (29) suggests the definition of a velocity field as follows [3]:

$$\overrightarrow{v}=\frac{\overrightarrow{j}}{{\rho }_p}=\frac{\hslash }{2mi{\rho }_p}[{\psi }^{†}\overrightarrow{\nabla }\psi -(\overrightarrow{\nabla }{\psi }^{†})\psi ]-k\overrightarrow{A}.$$

(31)

Holland [3] has suggested the following representation of the spinor:

$$\psi =R{e}^{i\frac{\chi }{2}}\left(\begin{array}{c}cos\left(\frac{\theta }{2}\right)e^{i\frac{\varphi }{2}}\\ isin\left(\frac{\theta }{2}\right)e^{-i\frac{\varphi }{2}}\end{array}\right)\equiv \left(\begin{array}{c}{\psi }_{\uparrow }\\ {\psi }_{\downarrow }\end{array}\right).$$

(32)

In terms of this representation the density is given as:

$$R^2={\psi }^{†}\psi ={\rho }_p\Rightarrow R=\sqrt{{\rho }_p}.$$

(33)

The mass density is given as:

$$\widehat{\rho }=m{\psi }^{†}\psi =m{R}^2=m{\rho }_p.$$

(34)

The probability amplitudes for spin up and spin down electrons are given by:

$$a_{\uparrow }=\left|{\psi }_{\uparrow }\right|=R\left|cos\frac{\theta }{2}\right|,a_{\downarrow }=\left|{\psi }_{\downarrow }\right|=R\left|sin\frac{\theta }{2}\right|$$

(35)

Let us now look at the expectation value of the spin:

$$<\frac{\hslash }{2}\overrightarrow{\sigma }>=\frac{\hslash }{2}\int {\psi }^{†}\overrightarrow{\sigma }\psi d^{3}x=\frac{\hslash }{2}\int \left(\frac{{\psi }^{†}\overrightarrow{\sigma }\psi }{{\rho }_p}\right){\rho }_p{d}^{3}x$$

(36)

The spin density can be calculated using the representation given in equation (32) as:

$$\widehat{s}\equiv \frac{{\psi }^{†}\overrightarrow{\sigma }\psi }{{\rho }_p}=(sin\theta sin\varphi ,sin\theta cos\varphi ,cos\theta ),\left|\widehat{s}\right|=\sqrt{\widehat{s}·\widehat{s}}=1.$$

(37)

This gives an easy physical interpretation to the variables $\theta ,\varphi$ as angles which describe the projection of the spin density on the axes. $\theta$ is the elevation angle of the spin density vector and $\varphi$ is the azimuthal angle of the same. The velocity field can now be calculated by inserting $\psi$ given in equation (32) into equation (31):

$$\overrightarrow{v}=\frac{\hslash }{2m}(\overrightarrow{\nabla }\chi +cos\theta \overrightarrow{\nabla }\varphi )-k\overrightarrow{A}.$$

(38)

We are now in a position to calculate the material derivative of the velocity and obtain the equation of motion for a particle with ([3] p. 393 equation (9.3.19)):

$$\frac{d\overrightarrow{v}}{dt}=-\overrightarrow{\nabla }(\frac{Q}{m})-{\left(\frac{\hslash }{2m}\right)}^2\frac{1}{{\rho }_p}{\partial }_k({\rho }_p\overrightarrow{\nabla }{\widehat{s}}_j{\partial }_k{\widehat{s}}_j)+k(\overrightarrow{E}+\overrightarrow{v}\times \overrightarrow{B})-\frac{\mu }{m}(\overrightarrow{\nabla }{B}_j){\widehat{s}}_j.$$

(39)

The Pauli equation of motion differs from the classical equation motion and the Schrödinger equation of motion. In addition to the Schrödinger quantum force correction we have an additional spin quantum force correction:

$${\overrightarrow{F}}_{QS}\equiv -\frac{{\hslash }^2}{4m}\frac{1}{{\rho }_p}{\partial }_k({\rho }_p\overrightarrow{\nabla }{\widehat{s}}_j{\partial }_k{\widehat{s}}_j)=-\frac{{\hslash }^2}{4m}\left[{\partial }_k(\overrightarrow{\nabla }{\widehat{s}}_j{\partial }_k{\widehat{s}}_j)+\frac{{\partial }_k{\rho }_p}{{\rho }_p}\overrightarrow{\nabla }{\widehat{s}}_j{\partial }_k{\widehat{s}}_j\right]$$

(40)

as well as a term characterizing the interaction of the spin with a gradient of the magnetic field.

$${\overrightarrow{F}}_{gradBS}\equiv -\mu (\overrightarrow{\nabla }{B}_j)s_j$$

(41)

As both the upper and lower spin components of the wave function are expanding in free space the gradients which appear in ${\overrightarrow{F}}_{QS}$ will tend to diminish for any macroscopic scale making this force negligible. To estimate the condition qualitatively we introduce the typical spin length:

$$L_s={\min }_{i\in \{1,2,3\}}{\left|\overrightarrow{\nabla }{\widehat{s}}_i\right|}^{-1}$$

(42)

Using the above definition we may estimate the spin quantum force:

$$F_{QS}\approx \frac{{\hslash }^2}{4m}[\frac{1}{L_s^3}+\frac{1}{L_s^2{L}_R}]=\frac{{\hslash }^2}{4m{L}_s^2}[\frac{1}{L_s}+\frac{1}{L_R}]$$

(43)

this suggested the definition of the hybrid typical length:

$$L_{sR}={[\frac{1}{L_s}+\frac{1}{L_R}]}^{-1}=\left\{\begin{array}{cc}{L}_s& L_s\ll L_R\\ L_R& L_R\ll L_s\end{array}\right..$$

(44)

In terms of this typical length we may write:

$$F_{QS}\approx \frac{{\hslash }^2}{4m{L}_s^2{L}_{sR}}$$

(45)

Thus the conditions for a classical trajectory become:

$$F_{QS}\ll F_L\Rightarrow L_s^2{L}_{sR}\gg \frac{{\hslash }^2}{4mF}\Rightarrow L_s\gg \left\{\begin{array}{cc}{\left(\frac{{\hslash }^2}{4m{F}_L}\right)}^{\frac{1}{3}}& L_s\ll L_R\\ {\left(\frac{{\hslash }^2}{4m{F}_L{L}_R}\right)}^{\frac{1}{2}}& L_R\ll L_s\end{array}\right..$$

(46)

Another important equation derived from equation (27) is the equation of motion for the spin orientation vector ([3] p. 392 equation (9.3.16)):

$$\frac{d\widehat{s}}{dt}=\frac{2\mu }{\hslash }{\overrightarrow{B}}_{eff}\times \widehat{s},{\overrightarrow{B}}_{eff}=\overrightarrow{B}-\frac{{\hslash }^2}{4\mu m{R}^2}{\partial }_i(\rho {\partial }_i\widehat{s})$$

(47)

The quantum correction to the magnetic field explains [3] why a spin pick up the orientation of the filed in a Stern-Gerlach experiment instead of precessing around it as a classical magnetic dipole would.

5. Ehrenfest Theorem

The above electron equations of motion are not accepted by all quantum physicists. Physicist who follow the Copenhagen school of quantum mechanics declare that a quantum electron does not have a trajectory. However, all quantum physicists agree that one can describe the trajectory of the expectation value of various operators such as position and momentum associated with the electron’s trajectory. This calculation is done through Ehrenfest Theorem [31]. The theorem states that for every quantum operator $A_o$ with expectation value:

$$<A_o>=\int d^{3}x{\psi }^{†}{A}_o\psi$$

(48)

the following equality holds:

$$\frac{d<A_o>}{dt}=<\frac{\partial A_o}{\partial t}>+\frac{1}{i\hslash }<[A_o,\widehat{H}]>,[A_o,\widehat{H}]\equiv A_o\widehat{H}-\widehat{H}{A}_o.$$

(49)

The position and velocity operators defined as [31]:

$${\overrightarrow{x}}_o\equiv \overrightarrow{x},{\overrightarrow{v}}_o\equiv \frac{1}{m}({\overrightarrow{p}}_o-e\overrightarrow{A})=\frac{1}{m}(-i\hslash \overrightarrow{\nabla }-e\overrightarrow{A})$$

(50)

For a Schrödinger’s electron (that is without spin) the following results are obtained by Griffiths [31] by inserting the above operators into equation (49):

$$\begin{array}{ccc}\hfill \frac{d<{\overrightarrow{x}}_o>}{dt}& =& \frac{1}{i\hslash }<[\overrightarrow{x},{\widehat{H}}_S]>=<{\overrightarrow{v}}_o>\hfill \\ \hfill \frac{d<{\overrightarrow{v}}_o>}{dt}& =& <\frac{\partial {\overrightarrow{v}}_o}{\partial t}>+\frac{1}{i\hslash }<[{\overrightarrow{v}}_o,{\widehat{H}}_S]>\hfill \\ & =& \frac{k}{2}<{\overrightarrow{v}}_o\times \overrightarrow{B}-\overrightarrow{B}\times {\overrightarrow{v}}_o>+k<\overrightarrow{E}>\hfill \end{array}$$

(51)

Thus the electron position expectation value satisfies the equation:

$$\frac{d^2<{\overrightarrow{x}}_o>}{d{t}^2}=\frac{k}{2}<{\overrightarrow{v}}_o\times \overrightarrow{B}-\overrightarrow{B}\times {\overrightarrow{v}}_o>+k<\overrightarrow{E}>$$

(52)

this equation resembles the classical and quantum equations of motion but also differs from them in many important aspects. First let us compare it with equation (13), let us also assume that $\overrightarrow{B}=0$. In this case:

$$\frac{d^2<\overrightarrow{x}>}{d{t}^2}=k<\overrightarrow{E}(\overrightarrow{x},t)>\ne k\overrightarrow{E}(<\overrightarrow{x}>,t),$$

(53)

thus as noted by many authors, even in this case the expectation value equations differ from the classical equation of motion except for a very restrictive class of linear electric fields. The difference is even more pronounced for the case $\overrightarrow{B}\ne 0$ which only takes a conventional "Lorentz force" form for a constant magnetic field $\overrightarrow{B}$ [31]:

$$\frac{d^2<\overrightarrow{x}>}{d{t}^2}=k(<{\overrightarrow{v}}_o>\times \overrightarrow{B}+<\overrightarrow{E}>).$$

(54)

Comparing equation (52) with the quantum motion equation (23), we see that the expectation value equation does not contain a quantum force term, which is a further justification to our assumption that this term may be neglected on macroscopic scales.

Let us now turn our attention to the more realistic Pauli electron which does posses spin, equation (51) now take the form:

$$\begin{array}{ccc}& & \frac{d<{\overrightarrow{x}}_o>}{dt}=\frac{1}{i\hslash }<[\overrightarrow{x},{\widehat{H}}_P]>=\frac{1}{i\hslash }<[\overrightarrow{x},{\widehat{H}}_S]>=<{\overrightarrow{v}}_o>\hfill \\ & & \frac{d<{\overrightarrow{v}}_o>}{dt}=<\frac{\partial {\overrightarrow{v}}_o}{\partial t}>+\frac{1}{i\hslash }<[{\overrightarrow{v}}_o,{\widehat{H}}_P]>\hfill \\ & & =<\frac{\partial {\overrightarrow{v}}_o}{\partial t}>+\frac{1}{i\hslash }<[{\overrightarrow{v}}_o,{\widehat{H}}_S]>+\frac{1}{i\hslash }<[{\overrightarrow{v}}_o,\mu \overrightarrow{B}·\overrightarrow{\sigma }]>\hfill \\ & & =\frac{k}{2}<{\overrightarrow{v}}_o\times \overrightarrow{B}-\overrightarrow{B}\times {\overrightarrow{v}}_o>+k<\overrightarrow{E}>+\frac{\mu }{i\hslash }<[{\overrightarrow{v}}_o,B_i]{\sigma }_i>,\hfill \end{array}$$

(55)

in which in the last term we use the Einstein summation convention. However:

$$[{\overrightarrow{v}}_o,B_i]=\frac{1}{m}[-i\hslash \overrightarrow{\nabla }-e\overrightarrow{A},B_i]=-\frac{i\hslash }{m}[\overrightarrow{\nabla },B_i]$$

(56)

The value of this commutation relation can be deduced by operating with the above operator on an arbitrary wave function $\psi$.

$$[{\overrightarrow{v}}_o,B_i]\psi =-\frac{i\hslash }{m}[\overrightarrow{\nabla },B_i]\psi =-\frac{i\hslash }{m}\left(\overrightarrow{\nabla }(B_i\psi )-B_i\overrightarrow{\nabla }\psi \right)=-\frac{i\hslash }{m}(\overrightarrow{\nabla }{B}_i)\psi$$

(57)

Or in operator jargon:

$$[{\overrightarrow{v}}_o,B_i]=-\frac{i\hslash }{m}(\overrightarrow{\nabla }{B}_i)$$

(58)

It thus follows that:

$$\frac{\mu }{i\hslash }<[{\overrightarrow{v}}_o,B_i]{\sigma }_i>=-\frac{\mu }{m}<(\overrightarrow{\nabla }{B}_i){\sigma }_i>=-\frac{\mu }{m}<(\overrightarrow{\nabla }{B}_i){\widehat{s}}_i>$$

(59)

in which we have used equation (37). Inserting equation (59) into equation (55):

$$\frac{d^2<{\overrightarrow{x}}_o>}{d{t}^2}=\frac{d<{\overrightarrow{v}}_o>}{dt}=\frac{k}{2}<{\overrightarrow{v}}_o\times \overrightarrow{B}-\overrightarrow{B}\times {\overrightarrow{v}}_o>+k<\overrightarrow{E}>-\frac{\mu }{m}<{\widehat{s}}_i\overrightarrow{\nabla }{B}_i>$$

(60)

comparing the above equation to equation (39) it follows that the only quantum force surviving the expectation value averaging is the one describing the effect of the magnetic field gradient on the spin vector which is in accordance with what should be expected in the macroscopic limit thus leading to the Stern-Gerlach experiment.

6. Conclusion: Spin orientation and the Stern-Gerlach experiment

We have seen how the Ehrenfest theorem approach causes "quantum forces" to disappear, those force disappear also in Bohm’s approach if one consider macroscopic scales of propagation.

The Stern-Gerlach experiment is an example of using the spin force term given in equation (41) to separate spin up and spin down particles thus obtaining from a single ray of particles two spots. If the magnetic field gradient is dominant in a single direction (say z) we may write equation (41) as:

$$F_{gradBSz}=-\mu ({\partial }_z{B}_j){\widehat{s}}_j$$

(61)

hence depending on the values ${\widehat{s}}_j$ some particle will move up and some will move down creating two spots (see Figure 1).

The Stern-Gerlach experiment is usually performed with natural particle not with charged particles like electrons, the reason for this is that generally speaking the classical Lorentz forces are much stronger than the quantum spin force and thus the two spot effect is not observed. Holland shows by simulating equation (47) that the spins in a Stern-Gerlach rotate in the direction or opposite to the direction of the magnetic field depending on the trajectory of the particle, that is to which spot it belongs (see Figure 2 (Holland [3] figure 9.13)). Notice, however, that from an energy perspective the lowest energy belong to the case in which the spin (and thus its related magnetic dipole) point at the direction of the field. The energy value is given by the expectation value of the Hamiltonian:

$$E=<{\widehat{H}}_P>=<{\widehat{H}}_S>+\mu <\overrightarrow{B}·\overrightarrow{\sigma }>.$$

(62)

If the direction of the field $\overrightarrow{B}$ is defined as a z direction it follows that:

$$E=<{\widehat{H}}_S>+\mu <B_z{\sigma }_z>=<{\widehat{H}}_S>+\mu \int d^{3}x{B}_z(a_{\uparrow }^2-a_{\downarrow }^2).$$

(63)

So particles with a definite spin direction (up or down) may have an upper or lower energy depending on the value of $\mu$. For an electron $\mu$ is the Bohr magneton:

$$\mu ={\mu }_B=\frac{\left|e\right|\hslash }{2m}.$$

(64)

and thus lower spin electrons will have a lower energy, if $B_z$ is constant we may write this term in the "classical" form using a magnetic dipole:

$$E=<{\widehat{H}}_S>-\overrightarrow{\mu }·\overrightarrow{B},\overrightarrow{\mu }=-\mu \int d^{3}x(a_{\uparrow }^2-a_{\downarrow }^2)\widehat{z}=\mu \int d^{3}x(a_{\downarrow }^2-a_{\uparrow }^2)\widehat{z}.$$

(65)

hence the magnetic dipole will point in the direction of the field for a lower energy configurations (spin down) and in the opposite direction for the higher energy configuration (spin up). As systems tend to relax to their lower energy state, one may ask why do the particles in a Stern-Gerlach experiment do not relax to the spin down configuration and instead split to beams of spin up and spin down with about the same size (see Figure 2)? The answer may be connected to the fact that in this type of experiment the electrons feel the magnetic field for only a short while and do not have enough time to relax to their minimum energy configurations.

This is not the case in NMR and MRI experiment in which the magnetic dipoles are under the influence of a strong magnetic field, for a long duration. In those cases the magnetization defined as:

$$\overrightarrow{M}={\rho }_P\overrightarrow{\mu }$$

(66)

Satisfies the Bloch phenomenological equations:

$$\begin{array}{ccc}\hfill \frac{d{M}_x}{dt}& =& {\gamma }_{gyro}{(\overrightarrow{M}\times \overrightarrow{B})}_x-\frac{M_x}{T_2}\hfill \\ \hfill \frac{d{M}_y}{dt}& =& {\gamma }_{gyro}{(\overrightarrow{M}\times \overrightarrow{B})}_y-\frac{M_y}{T_2}\hfill \\ \hfill \frac{d{M}_z}{dt}& =& {\gamma }_{gyro}{(\overrightarrow{M}\times \overrightarrow{B})}_z-\frac{M_z-M_0}{T_1}\hfill \end{array}$$

(67)

in the above $T_1$ and $T_2$ are typical relaxation times and $\gamma$ is a gyromagnetic ratio which a for an electron takes the value:

$${\gamma }_{gyro}=\frac{2{\mu }_B}{\hslash }$$

(68)

The magnetization satisfying the above equation eventually relaxes to the direction of the field which is the minimal energy configuration.

As a future direction to the current research it may be interesting to study the same problem for a fully relativistic electron, however, this will require using a Dirac equation rather than Pauli’s equation.