Preprint

Article

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

Altmetrics

Downloads

129

Views

Comments

A peer-reviewed article of this preprint also exists.

Asher Yahalom^*

This version is not peer-reviewed

Submitted:

21 December 2022

Posted:

12 January 2023

You are already at the latest version

Alerts

Abstract

Electrons moving at slow speeds much lower that the speed of light are described by a wave function which is a solution of Pauli’s equation. This is a low velocity limit of the relativistic Dirac equation. Here we compare two approaches, one which is the more conservative Copenhagen’s interpretation denying a trajectory of the electron but allowing a trajectory to the electron expectation value through Ehrenfest theorem. The said expectation value is of course calculated using a solution of Pauli’s equation. A less orthodox approach is championed by Bohm, and attributes a velocity field to the electron also derived from the Pauli wave function. It is thus interesting to compare the trajectory followed by the electron according to Bohm and its expectation value according to Ehrenfest. Both similarities and differences will be considered.

Keywords:

Subject: Physical Sciences - Quantum Science and Technology

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

\vec{x} (t)

, mass m and charge e interacting with a given electromagnetic vector potential

\vec{A} (\vec{x}, t)

and scalar potential

φ (\vec{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{matrix} A & = & \int_{t 1}^{t 2} L d t, L = L_{0} + L_{i} \\ L_{0} & \equiv & \frac{1}{2} m v^{2}, L_{i} \equiv e (\vec{A} \cdot \vec{v} - φ), \vec{v} \equiv \frac{d \vec{x}}{d t} \equiv \dot{\vec{x}}, v = | \vec{v} | . \end{matrix}

(1)

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

δ L_{0} = m \dot{\vec{x}} \cdot δ \dot{\vec{x}} = \frac{d (m \vec{v} \cdot δ \vec{x})}{d t} - m \dot{\vec{v}} \cdot δ \vec{x}

(2)

\begin{matrix} δ L_{i} & = & e (δ \vec{A} \cdot \vec{v} + \vec{A} \cdot δ \dot{\vec{x}} - δ φ) \\ = & e (\partial_{k} \vec{A} \cdot \vec{v} δ x_{k} + \frac{d (\vec{A} \cdot δ \vec{x})}{d t} - δ \vec{x} \cdot \frac{d \vec{A}}{d t} - \vec{\nabla} φ \cdot δ \vec{x}), \end{matrix}

(3)

in the above

\partial_{k} \equiv \frac{\partial}{\partial x_{k}}

and

\vec{\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

\vec{A}

as:

\frac{d \vec{A} (\vec{x} (t), t)}{d t} = \partial_{t} \vec{A} + v_{l} \partial_{l} \vec{A}, \partial_{t} \equiv \frac{\partial}{\partial t} .

(4)

Thus, the variation

δ L_{i}

can be written in the following form:

δ L_{i} = \frac{d (e \vec{A} \cdot δ \vec{x})}{d t} + e [(\partial_{k} A_{l} - \partial_{l} A_{k}) v_{l} - \partial_{t} A_{k} - \partial_{k} φ] δ x_{k},

(5)

Defining the electric and magnetic fields in the standard way:

\vec{B} \equiv \vec{\nabla} \times \vec{A}, \vec{E} \equiv - \partial_{t} \vec{A} - \vec{\nabla} φ,

(6)

it follows that:

ϵ_{k l n} B_{n} = \partial_{k} A_{l} - \partial_{l} A_{k}, E_{k} = - \partial_{t} A_{k} - \partial_{k} φ,

(7)

in which

ϵ_{k l n}

is the three index antisymmetric tensor. Thus, we may write

δ L_{i}

as:

δ L_{i} = \frac{d (e \vec{A} \cdot δ \vec{x})}{d t} + e [ϵ_{k l n} v_{l} B_{n} + E_{k}] δ x_{k} = \frac{d (e \vec{A} \cdot δ \vec{x})}{d t} + e [\vec{v} \times \vec{B} + \vec{E}] \cdot δ \vec{x} .

(8)

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

{\vec{F}}_{L} \equiv e [\vec{v} \times \vec{B} + \vec{E}]

(9)

to write:

δ L_{i} = \frac{d (e \vec{A} \cdot δ \vec{x})}{d t} + {\vec{F}}_{L} \cdot δ \vec{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:

δ L = δ L_{0} + δ L_{i} = \frac{d ((m \vec{v} + e \vec{A}) \cdot δ \vec{x})}{d t} + (- m \dot{\vec{v}} + {\vec{F}}_{L}) \cdot δ \vec{x} .

(11)

Thus, the variation of the action is:

δ A = \int_{t 1}^{t 2} δ L d t = {(m \vec{v} + e \vec{A}) \cdot δ \vec{x}|}_{t 1}^{t 2} - \int_{t 1}^{t 2} (m \dot{\vec{v}} - {\vec{F}}_{L}) \cdot δ \vec{x} d t .

(12)

Since the classical trajectory is such that the variation of the action on it vanishes for a small modification of the trajectory

δ \vec{x}

that vanishes at

t 1

and

t 2

but is otherwise arbitrary it follows that:

m \dot{\vec{v}} = {\vec{F}}_{L} = e [\vec{v} \times \vec{B} + \vec{E}] \Rightarrow \dot{\vec{v}} = \frac{e}{m} [\vec{v} \times \vec{B} + \vec{E}] .

(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{\vec{v}} = k [\vec{v} \times \vec{B} + \vec{E}] .

(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:

{\vec{A}}^{'} = \vec{A} + \vec{\nabla} Λ, φ^{'} = φ - \partial_{t} Λ .

(15)

we obtain the same fields:

{\vec{B}}^{'} = \vec{\nabla} \times {\vec{A}}^{'} = \vec{\nabla} \times \vec{A} = \vec{B}, {\vec{E}}^{'} = - \partial_{t} {\vec{A}}^{'} - \vec{\nabla} φ^{'} = - \partial_{t} \vec{A} - \vec{\nabla} φ = \vec{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 ℏ \dot{ψ} = {\hat{H}}_{S} ψ, {\hat{H}}_{S} = - \frac{1}{2 m} {(ℏ \vec{\nabla} - i e \vec{A})}^{2} + e φ

(17)

in the above

i = \sqrt{- 1}

and

ψ

is the complex wave function.

\dot{ψ} = \frac{\partial ψ}{\partial t}

is the partial time derivative of the wave function.

ℏ = \frac{h}{2 π}

is Planck’s constant divided by

2 π

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

ϕ

ψ = a e^{i ϕ} .

(18)

We define the velocity field:

{\vec{v}}_{S} = \frac{ℏ}{m} \vec{\nabla} ϕ - \frac{e}{m} \vec{A}

(19)

and the mass density is defined as:

\hat{ρ} = m a^{2} .

(20)

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

\frac{\partial \hat{ρ}}{\partial t} + \vec{\nabla} \cdot (\hat{ρ} {\vec{v}}_{S}) = 0

(21)

Hence

{\vec{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

ρ = e a^{2}

. The equation for the phase

ϕ

derived from equation (17) is as follows:

ℏ \frac{\partial ϕ}{\partial t} + \frac{1}{2 m} {(ℏ \vec{\nabla} ϕ - e \vec{A})}^{2} + e φ = \frac{ℏ^{2} \nabla^{2} a}{2 m a}

(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 {\vec{v}}_{S}}{d t} = \frac{\partial {\vec{v}}_{S}}{\partial t} + ({\vec{v}}_{S} \cdot \vec{\nabla}) {\vec{v}}_{S} = - \vec{\nabla} \frac{Q}{m} + k (\vec{E} + \vec{v} \times \vec{B})

(23)

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

Q = - \frac{ℏ^{2}}{2 m} \frac{{\vec{\nabla}}^{2} \sqrt{\hat{ρ}}}{\sqrt{\hat{ρ}}} .

(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

ℏ \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:

{\vec{F}}_{Q} = - \vec{\nabla} Q ≃ \frac{ℏ^{2}}{2 m L_{R}^{3}}, L_{R} ≃ \frac{R}{| \vec{\nabla} R |}

(25)

in which

L_{R}

is the typical length of the amplitudes gradient. Thus:

| F_{Q} | < < | F_{L} | \Rightarrow L_{R} > > L_{R c} = {(\frac{ℏ^{2}}{2 m F_{L}})}^{\frac{1}{3}} .

(26)

in which

{\vec{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 ℏ \dot{ψ} = {\hat{H}}_{P} ψ, {\hat{H}}_{P} = - \frac{ℏ^{2}}{2 m} {[\vec{\nabla} - \frac{i e}{ℏ} \vec{A}]}^{2} + μ \vec{B} \cdot \vec{σ} + e φ = {\hat{H}}_{S} I + μ \vec{B} \cdot \vec{σ}

(27)

ψ

here is a two dimensional complex column vector (also denoted as spinor),

{\hat{H}}_{P}

is a two dimensional hermitian operator matrix,

μ

is the magnetic moment of the particle, and I is a two dimensional unit matrix.

\vec{σ}

is a vector of two dimensional Pauli matrices which can be represented as follows:

σ_{1} = (\begin{matrix} 0 & 1 \\ 1 & 0 \end{matrix}), σ_{2} = (\begin{matrix} 0 & - i \\ i & 0 \end{matrix}), σ_{3} = (\begin{matrix} 1 & 0 \\ 0 & - 1 \end{matrix}) .

(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

ψ

satisfying equation (27) must also satisfy a continuity equation of the form:

\frac{\partial ρ_{p}}{\partial t} + \vec{\nabla} \cdot \vec{j} = 0 .

(29)

In the above:

ρ_{p} = ψ^{†} ψ, \vec{j} = \frac{ℏ}{2 m i} [ψ^{†} \vec{\nabla} ψ - (\vec{\nabla} ψ^{†}) ψ] - k \vec{A} ρ_{p} .

(30)

The symbol

ψ^{†}

represents a row spinor (the transpose) whose components are equal to the complex conjugate of the column spinor

ψ

. Comparing the standard continuity equation to equation (29) suggests the definition of a velocity field as follows [3]:

\vec{v} = \frac{\vec{j}}{ρ_{p}} = \frac{ℏ}{2 m i ρ_{p}} [ψ^{†} \vec{\nabla} ψ - (\vec{\nabla} ψ^{†}) ψ] - k \vec{A} .

(31)

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

ψ = R e^{i \frac{χ}{2}} (\begin{matrix} cos (\frac{θ}{2}) e^{i \frac{ϕ}{2}} \\ i sin (\frac{θ}{2}) e^{- i \frac{ϕ}{2}} \end{matrix}) \equiv (\begin{matrix} ψ_{↑} \\ ψ_{↓} \end{matrix}) .

(32)

In terms of this representation the density is given as:

R^{2} = ψ^{†} ψ = ρ_{p} \Rightarrow R = \sqrt{ρ_{p}} .

(33)

The mass density is given as:

\hat{ρ} = m ψ^{†} ψ = m R^{2} = m ρ_{p} .

(34)

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

a_{↑} = |ψ_{↑}| = R |cos \frac{θ}{2}|, a_{↓} = |ψ_{↓}| = R |sin \frac{θ}{2}|

(35)

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

< \frac{ℏ}{2} \vec{σ} > = \frac{ℏ}{2} \int ψ^{†} \vec{σ} ψ d^{3} x = \frac{ℏ}{2} \int (\frac{ψ^{†} \vec{σ} ψ}{ρ_{p}}) ρ_{p} d^{3} x

(36)

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

\hat{s} \equiv \frac{ψ^{†} \vec{σ} ψ}{ρ_{p}} = (sin θ sin ϕ, sin θ cos ϕ, cos θ), | \hat{s} | = \sqrt{\hat{s} \cdot \hat{s}} = 1 .

(37)

This gives an easy physical interpretation to the variables

θ, ϕ

as angles which describe the projection of the spin density on the axes.

θ

is the elevation angle of the spin density vector and

ϕ

is the azimuthal angle of the same. The velocity field can now be calculated by inserting

ψ

given in equation (32) into equation (31):

\vec{v} = \frac{ℏ}{2 m} (\vec{\nabla} χ + cos θ \vec{\nabla} ϕ) - k \vec{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 \vec{v}}{d t} = - \vec{\nabla} (\frac{Q}{m}) - {(\frac{ℏ}{2 m})}^{2} \frac{1}{ρ_{p}} \partial_{k} (ρ_{p} \vec{\nabla} {\hat{s}}_{j} \partial_{k} {\hat{s}}_{j}) + k (\vec{E} + \vec{v} \times \vec{B}) - \frac{μ}{m} (\vec{\nabla} B_{j}) {\hat{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:

{\vec{F}}_{Q S} \equiv - \frac{ℏ^{2}}{4 m} \frac{1}{ρ_{p}} \partial_{k} (ρ_{p} \vec{\nabla} {\hat{s}}_{j} \partial_{k} {\hat{s}}_{j}) = - \frac{ℏ^{2}}{4 m} [\partial_{k} (\vec{\nabla} {\hat{s}}_{j} \partial_{k} {\hat{s}}_{j}) + \frac{\partial_{k} ρ_{p}}{ρ_{p}} \vec{\nabla} {\hat{s}}_{j} \partial_{k} {\hat{s}}_{j}]

(40)

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

{\vec{F}}_{g r a d B S} \equiv - μ (\vec{\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

{\vec{F}}_{Q S}

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}} {| \vec{\nabla} {\hat{s}}_{i} |}^{- 1}

(42)

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

F_{Q S} \approx \frac{ℏ^{2}}{4 m} [\frac{1}{L_{s}^{3}} + \frac{1}{L_{s}^{2} L_{R}}] = \frac{ℏ^{2}}{4 m L_{s}^{2}} [\frac{1}{L_{s}} + \frac{1}{L_{R}}]

(43)

this suggested the definition of the hybrid typical length:

L_{s R} = {[\frac{1}{L_{s}} + \frac{1}{L_{R}}]}^{- 1} = \{\begin{matrix} L_{s} & L_{s} ≪ L_{R} \\ L_{R} & L_{R} ≪ L_{s} \end{matrix} .

(44)

In terms of this typical length we may write:

F_{Q S} \approx \frac{ℏ^{2}}{4 m L_{s}^{2} L_{s R}}

(45)

Thus the conditions for a classical trajectory become:

F_{Q S} ≪ F_{L} \Rightarrow L_{s}^{2} L_{s R} ≫ \frac{ℏ^{2}}{4 m F} \Rightarrow L_{s} ≫ \{\begin{matrix} {(\frac{ℏ^{2}}{4 m F_{L}})}^{\frac{1}{3}} & L_{s} ≪ L_{R} \\ {(\frac{ℏ^{2}}{4 m F_{L} L_{R}})}^{\frac{1}{2}} & L_{R} ≪ L_{s} \end{matrix} .

(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 \hat{s}}{d t} = \frac{2 μ}{ℏ} {\vec{B}}_{e f f} \times \hat{s}, {\vec{B}}_{e f f} = \vec{B} - \frac{ℏ^{2}}{4 μ m R^{2}} \partial_{i} (ρ \partial_{i} \hat{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 ψ^{†} A_{o} ψ

(48)

the following equality holds:

\frac{d < A_{o} >}{d t} = < \frac{\partial A_{o}}{\partial t} > + \frac{1}{i ℏ} < [A_{o}, \hat{H}] >, [A_{o}, \hat{H}] \equiv A_{o} \hat{H} - \hat{H} A_{o} .

(49)

The position and velocity operators defined as [31]:

{\vec{x}}_{o} \equiv \vec{x}, {\vec{v}}_{o} \equiv \frac{1}{m} ({\vec{p}}_{o} - e \vec{A}) = \frac{1}{m} (- i ℏ \vec{\nabla} - e \vec{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{matrix} \frac{d < {\vec{x}}_{o} >}{d t} & = & \frac{1}{i ℏ} < [\vec{x}, {\hat{H}}_{S}] > = < {\vec{v}}_{o} > \\ \frac{d < {\vec{v}}_{o} >}{d t} & = & < \frac{\partial {\vec{v}}_{o}}{\partial t} > + \frac{1}{i ℏ} < [{\vec{v}}_{o}, {\hat{H}}_{S}] > \\ = & \frac{k}{2} < {\vec{v}}_{o} \times \vec{B} - \vec{B} \times {\vec{v}}_{o} > + k < \vec{E} > \end{matrix}

(51)

Thus the electron position expectation value satisfies the equation:

\frac{d^{2} < {\vec{x}}_{o} >}{d t^{2}} = \frac{k}{2} < {\vec{v}}_{o} \times \vec{B} - \vec{B} \times {\vec{v}}_{o} > + k < \vec{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

\vec{B} = 0

. In this case:

\frac{d^{2} < \vec{x} >}{d t^{2}} = k < \vec{E} (\vec{x}, t) > \neq k \vec{E} (< \vec{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

\vec{B} \neq 0

which only takes a conventional "Lorentz force" form for a constant magnetic field

\vec{B}

[31]:

\frac{d^{2} < \vec{x} >}{d t^{2}} = k (< {\vec{v}}_{o} > \times \vec{B} + < \vec{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{matrix} \frac{d < {\vec{x}}_{o} >}{d t} = \frac{1}{i ℏ} < [\vec{x}, {\hat{H}}_{P}] > = \frac{1}{i ℏ} < [\vec{x}, {\hat{H}}_{S}] > = < {\vec{v}}_{o} > \\ \frac{d < {\vec{v}}_{o} >}{d t} = < \frac{\partial {\vec{v}}_{o}}{\partial t} > + \frac{1}{i ℏ} < [{\vec{v}}_{o}, {\hat{H}}_{P}] > \\ = < \frac{\partial {\vec{v}}_{o}}{\partial t} > + \frac{1}{i ℏ} < [{\vec{v}}_{o}, {\hat{H}}_{S}] > + \frac{1}{i ℏ} < [{\vec{v}}_{o}, μ \vec{B} \cdot \vec{σ}] > \\ = \frac{k}{2} < {\vec{v}}_{o} \times \vec{B} - \vec{B} \times {\vec{v}}_{o} > + k < \vec{E} > + \frac{μ}{i ℏ} < [{\vec{v}}_{o}, B_{i}] σ_{i} >, \end{matrix}

(55)

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

[{\vec{v}}_{o}, B_{i}] = \frac{1}{m} [- i ℏ \vec{\nabla} - e \vec{A}, B_{i}] = - \frac{i ℏ}{m} [\vec{\nabla}, B_{i}]

(56)

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

ψ

[{\vec{v}}_{o}, B_{i}] ψ = - \frac{i ℏ}{m} [\vec{\nabla}, B_{i}] ψ = - \frac{i ℏ}{m} (\vec{\nabla} (B_{i} ψ) - B_{i} \vec{\nabla} ψ) = - \frac{i ℏ}{m} (\vec{\nabla} B_{i}) ψ

(57)

Or in operator jargon:

[{\vec{v}}_{o}, B_{i}] = - \frac{i ℏ}{m} (\vec{\nabla} B_{i})

(58)

It thus follows that:

\frac{μ}{i ℏ} < [{\vec{v}}_{o}, B_{i}] σ_{i} > = - \frac{μ}{m} < (\vec{\nabla} B_{i}) σ_{i} > = - \frac{μ}{m} < (\vec{\nabla} B_{i}) {\hat{s}}_{i} >

(59)

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

\frac{d^{2} < {\vec{x}}_{o} >}{d t^{2}} = \frac{d < {\vec{v}}_{o} >}{d t} = \frac{k}{2} < {\vec{v}}_{o} \times \vec{B} - \vec{B} \times {\vec{v}}_{o} > + k < \vec{E} > - \frac{μ}{m} < {\hat{s}}_{i} \vec{\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_{g r a d B S z} = - μ (\partial_{z} B_{j}) {\hat{s}}_{j}

(61)

hence depending on the values

{\hat{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 = < {\hat{H}}_{P} > = < {\hat{H}}_{S} > + μ < \vec{B} \cdot \vec{σ} > .

(62)

If the direction of the field

\vec{B}

is defined as a z direction it follows that:

E = < {\hat{H}}_{S} > + μ < B_{z} σ_{z} > = < {\hat{H}}_{S} > + μ \int d^{3} x B_{z} (a_{↑}^{2} - a_{↓}^{2}) .

(63)

So particles with a definite spin direction (up or down) may have an upper or lower energy depending on the value of

μ

. For an electron

μ

is the Bohr magneton:

μ = μ_{B} = \frac{| e | ℏ}{2 m} .

(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 = < {\hat{H}}_{S} > - \vec{μ} \cdot \vec{B}, \vec{μ} = - μ \int d^{3} x (a_{↑}^{2} - a_{↓}^{2}) \hat{z} = μ \int d^{3} x (a_{↓}^{2} - a_{↑}^{2}) \hat{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:

\vec{M} = ρ_{P} \vec{μ}

(66)

Satisfies the Bloch phenomenological equations:

\begin{matrix} \frac{d M_{x}}{d t} & = & γ_{g y r o} {(\vec{M} \times \vec{B})}_{x} - \frac{M_{x}}{T_{2}} \\ \frac{d M_{y}}{d t} & = & γ_{g y r o} {(\vec{M} \times \vec{B})}_{y} - \frac{M_{y}}{T_{2}} \\ \frac{d M_{z}}{d t} & = & γ_{g y r o} {(\vec{M} \times \vec{B})}_{z} - \frac{M_{z} - M_{0}}{T_{1}} \end{matrix}

(67)

in the above

T_{1}

and

T_{2}

are typical relaxation times and

γ

is a gyromagnetic ratio which a for an electron takes the value:

γ_{g y r o} = \frac{2 μ_{B}}{ℏ}

(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.

References

Kant, I. (1781). Critik der reinen Vernunft.
D. Bohm, Quantum Theory (Prentice Hall, New York, 1966) section 12.6.
P.R. Holland The Quantum Theory of Motion (Cambridge University Press, Cambridge, 1993).
D. Durr & S. Teufel Bohmian Mechanics: The Physics and Mathematics of Quantum Theory (Springer-Verlag, Berlin Heidelberg, 2009).
E. Madelung, Z. Phys., 40 322 (1926).
R. Englman and A. Yahalom "Complex States of Simple Molecular Systems" a chapter of the volume "The Role of Degenerate States in Chemistry" edited by M. Baer and G. Billing in Adv. Chem. Phys. Vol. 124 (John Wiley & Sons 2002). [Los-Alamos Archives physics/0406149]. [CrossRef]
W. Pauli (1927) Zur Quantenmechanik des magnetischen Elektrons Zeitschrift für Physik (43) 601-623.
R. Englman, A.Yahalom and M. Baer, J. Chem. Phys.109 6550 (1998).
R. Englman, A. Yahalom and M. Baer, Phys. Lett. A 251 223 (1999).
R. Englman and A. Yahalom, Phys. Rev. A 60 1802 (1999).
R. Englman, A.Yahalom and M. Baer, Eur. Phys. J. D 8 1 (2000).
R. Englman and A.Yahalom, Phys. Rev. B 61 2716 (2000).
R. Englman and A.Yahalom, Found. Phys. Lett. 13 329 (2000).
R. Englman and A.Yahalom, The Jahn Teller Effect: A Permanent Presence in the Frontiers of Science in M.D. Kaplan and G. Zimmerman (editors), Proceedings of the NATO Advanced Research Workshop, Boston, Sep. 2000 (Kluwer, Dordrecht, 2001).
M. Baer and R. Englman, Chem. Phys. Lett. 335 85 (2001).
A. Mebel, M. Baer, R. Englman and A. Yahalom, J.Chem. Phys. 115 3673 (2001).
R. Englman & A. Yahalom, "Signed Phases and Fields Associated with Degeneracies" Acta Phys. et Chim., 34-35, 283 (2002). [Los-Alamos Archives - quant-ph/0406194]. [CrossRef]
R. Englman, A. Yahalom and M. Baer,"Hierarchical Construction of Finite Diabatic Sets, By Mathieu Functions", Int. J. Q. Chemistry, 90, 266-272 (2002). [Los-Alamos Archives -physics/0406126]. [CrossRef]
R. Englman, A. Yahalom, M. Baer and A.M. Mebel "Some Experimental. and Calculational Consequences of Phases in Molecules with Multiple Conical Intersections" International Journal of Quantum Chemistry, 92, 135-151 (2003). [CrossRef]
R. Englman & A. Yahalom, "Phase Evolution in a Multi-Component System", Physical Review A, 67, 5, 054103-054106 (2003). [Los-Alamos Archives -quant-ph/0406195].
R. Englman & A. Yahalom, "Generalized "Quasi-classical" Ground State of an Interacting Doublet" Physical Review B, 69, 22, 224302 (2004). [Los-Alamos Archives - cond-mat/0406725].
Ehrenfest, P. (1927). "Bemerkung über die angenäherte Gültigkeit der klassischen Mechanik innerhalb der Quantenmechanik". Zeitschrift für Physik. 45 (7–8): 455–457. [CrossRef]
A. Yahalom "The Fluid Dynamics of Spin". Molecular Physics, Published online: 13 Apr 2018.http://dx.doi.org/10.1080/00268976.2018.1457808 (arXiv:1802.09331v1 [physics.flu-dyn]). [CrossRef]
E. Schrödinger, Ann. d. Phys. 81 109 (1926). English translation appears in E. Schrödinger, Collected Papers in Wave Mechanics (Blackie and Sons, London, 1928) p. 102.
R. A. Fisher Phil. Trans. R. Soc. London 222, 309.
L. Mandel and E. Wolf, Optical Coherence and Quantum Optics (University Press, Cambridge, 1995) section 3.1.
A. Yahalom "The Fluid Dynamics of Spin - a Fisher Information Perspective" arXiv:1802.09331v2 [cond-mat.] 6 Jul 2018. Proceedings of the Seventeenth Israeli - Russian Bi-National Workshop 2018 "The optimization of composition, structure and properties of metals, oxides, composites, nano and amorphous materials".
B. R. Frieden Science from Fisher Information: A Unification (Cambridge University Press, Cambridge, 2004).
Asher Yahalom "The Fluid Dynamics of Spin - a Fisher Information Perspective and Comoving Scalars" Chaotic Modeling and Simulation (CMSIM) 1: 17-30, 2020.
Yahalom, A. Fisher Information Perspective of Pauli’s Electron. Entropy 2022, 24, 1721. [CrossRef]
David J. Griffiths Introduction to Quantum Mechanics, Cambridge University Press; 3rd edition (August 16, 2018).

Figure 1. A schematics of Stern-Gerlach experiment. Neutral particles travelling through an inhomogeneous magnetic field, and being deflected up or down depending on their spin; (1) particle source, (2) beam of particles, (3) inhomogeneous magnetic field, (4) classically expected result (neglecting the quantum spin force), (5) observed result.

Figure 2. A Stern-Gerlach simulation (Holland [3] figure 9.13). Notice how the spin components align in opposite directions in each beam.

Figure 3. The evolution of magnetization towards relaxation, the tip of the magnetization vector is described by the orange line.

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.