Entropy when density operator has zero eigenvalues

September 20, 2015 phy1520 , , , , ,

[Click here for a PDF of this post with nicer formatting]

In the class notes and the text [1] the Von Neumann entropy is defined as

\begin{equation}\label{eqn:densityMatrixEntropy:20}
S = -\textrm{Tr} \rho \ln \rho.
\end{equation}

In one of our problems I had trouble evaluating this, having calculated a density operator matrix representation

\begin{equation}\label{eqn:densityMatrixEntropy:40}
\rho = E \wedge E^{-1},
\end{equation}

where

\begin{equation}\label{eqn:densityMatrixEntropy:60}
E = \inv{\sqrt{2}}
\begin{bmatrix}
1 & 1 \\
1 & -1
\end{bmatrix},
\end{equation}

and
\begin{equation}\label{eqn:densityMatrixEntropy:100}
\wedge =
\begin{bmatrix}
1 & 0 \\
0 & 0
\end{bmatrix}.
\end{equation}

The usual method of evaluating a function of a matrix is to assume the function has a power series representation, and that a similarity transformation of the form \( A = E \wedge E^{-1} \) is possible, so that

\begin{equation}\label{eqn:densityMatrixEntropy:80}
f(A) = E f(\wedge) E^{-1},
\end{equation}

however, when attempting to do this with the matrix of \ref{eqn:densityMatrixEntropy:40} leads to an undesirable result

\begin{equation}\label{eqn:densityMatrixEntropy:120}
\ln \rho =
\inv{2}
\begin{bmatrix}
1 & 1 \\
1 & -1
\end{bmatrix}
\begin{bmatrix}
\ln 1 & 0 \\
0 & \ln 0
\end{bmatrix}
\begin{bmatrix}
1 & 1 \\
1 & -1
\end{bmatrix}.
\end{equation}

The \( \ln 0 \) makes the evaluation of this matrix logarithm rather unpleasant. To give meaning to the entropy expression, we have to do two things, the first is treating the trace operation as a higher precedence than the logarithms that it contains. That is

\begin{equation}\label{eqn:densityMatrixEntropy:140}
\begin{aligned}
-\textrm{Tr} ( \rho \ln \rho )
&=
-\textrm{Tr} ( E \wedge E^{-1} E \ln \wedge E^{-1} ) \\
&=
-\textrm{Tr} ( E \wedge \ln \wedge E^{-1} ) \\
&=
-\textrm{Tr} ( E^{-1} E \wedge \ln \wedge ) \\
&=
-\textrm{Tr} ( \wedge \ln \wedge ) \\
&=
– \sum_k \wedge_{kk} \ln \wedge_{kk}.
\end{aligned}
\end{equation}

Now the matrix of the logarithm need not be evaluated, but we still need to give meaning to \( \wedge_{kk} \ln \wedge_{kk} \) for zero diagonal entries. This can be done by considering a limiting scenerio

\begin{equation}\label{eqn:densityMatrixEntropy:160}
\begin{aligned}
-\lim_{a \rightarrow 0} a \ln a
&=
-\lim_{x \rightarrow \infty} e^{-x} \ln e^{-x} \\
&=
\lim_{x \rightarrow \infty} x e^{-x} \\
&=
0.
\end{aligned}
\end{equation}

The entropy can now be expressed in the unambiguous form, summing over all the non-zero eigenvalues of the density operator

\begin{equation}\label{eqn:densityMatrixEntropy:180}
\boxed{
S = – \sum_{ \wedge_{kk} \ne 0} \wedge_{kk} \ln \wedge_{kk}.
}
\end{equation}

References

[1] Jun John Sakurai and Jim J Napolitano. Modern quantum mechanics. Pearson Higher Ed, 2014.

PHY1520H Graduate Quantum Mechanics. Lecture 1: Lighting review. Taught by Prof. Arun Paramekanti

September 17, 2015 phy1520 , , , , , , , , , , , , , ,

[Click here for a PDF of this post with nicer formatting]

Disclaimer

Peeter’s lecture notes from class. These may be incoherent and rough.

These are notes for the UofT course PHY1520, Graduate Quantum Mechanics, taught by Prof. Paramekanti, covering [1] chap. 1 content.

Classical mechanics

We’ll be talking about one body physics for most of this course. In classical mechanics we can figure out the particle trajectories using both of \( (\Br, \Bp \), where

\begin{equation}\label{eqn:qmLecture1:20}
\begin{aligned}
\ddt{\Br} &= \inv{m} \Bp \\
\ddt{\Bp} &= \spacegrad V
\end{aligned}
\end{equation}

A two dimensional phase space as sketched in fig. 1 shows the trajectory of a point particle subject to some equations of motion

lectureOnePhaseSpaceClassicalFig1

fig. 1. One dimensional classical phase space example

Quantum mechanics

For this lecture, we’ll work with natural units, setting

\begin{equation}\label{eqn:qmLecture1:480}
\boxed{
\Hbar = 1.
}
\end{equation}

In QM we are no longer allowed to think of position and momentum, but have to start asking about state vectors \( \ket{\Psi} \).

We’ll consider the state vector with respect to some basis, for example, in a position basis, we write

\begin{equation}\label{eqn:qmLecture1:40}
\braket{ x }{\Psi } = \Psi(x),
\end{equation}

a complex numbered “wave function”, the probability amplitude for a particle in \( \ket{\Psi} \) to be in the vicinity of \( x \).

We could also consider the state in a momentum basis

\begin{equation}\label{eqn:qmLecture1:60}
\braket{ p }{\Psi } = \Psi(p),
\end{equation}

a probability amplitude with respect to momentum \( p \).

More precisely,

\begin{equation}\label{eqn:qmLecture1:80}
\Abs{\Psi(x)}^2 dx \ge 0
\end{equation}

is the probability of finding the particle in the range \( (x, x + dx ) \). To have meaning as a probability, we require

\begin{equation}\label{eqn:qmLecture1:100}
\int_{-\infty}^\infty \Abs{\Psi(x)}^2 dx = 1.
\end{equation}

The average position can be calculated using this probability density function. For example

\begin{equation}\label{eqn:qmLecture1:120}
\expectation{x} = \int_{-\infty}^\infty \Abs{\Psi(x)}^2 x dx,
\end{equation}

or
\begin{equation}\label{eqn:qmLecture1:140}
\expectation{f(x)} = \int_{-\infty}^\infty \Abs{\Psi(x)}^2 f(x) dx.
\end{equation}

Similarly, calculation of an average of a function of momentum can be expressed as

\begin{equation}\label{eqn:qmLecture1:160}
\expectation{f(p)} = \int_{-\infty}^\infty \Abs{\Psi(p)}^2 f(p) dp.
\end{equation}

Transformation from a position to momentum basis

We have a problem, if we which to compute an average in momentum space such as \( \expectation{p} \), when given a wavefunction \( \Psi(x) \).

How do we convert

\begin{equation}\label{eqn:qmLecture1:180}
\Psi(p)
\stackrel{?}{\leftrightarrow}
\Psi(x),
\end{equation}

or equivalently
\begin{equation}\label{eqn:qmLecture1:200}
\braket{p}{\Psi}
\stackrel{?}{\leftrightarrow}
\braket{x}{\Psi}.
\end{equation}

Such a conversion can be performed by virtue of an the assumption that we have a complete orthonormal basis, for which we can introduce identity operations such as

\begin{equation}\label{eqn:qmLecture1:220}
\int_{-\infty}^\infty dp \ket{p}\bra{p} = 1,
\end{equation}

or
\begin{equation}\label{eqn:qmLecture1:240}
\int_{-\infty}^\infty dx \ket{x}\bra{x} = 1
\end{equation}

Some interpretations:

  1. \( \ket{x_0} \leftrightarrow \text{sits at} x = x_0 \)
  2. \( \braket{x}{x’} \leftrightarrow \delta(x – x’) \)
  3. \( \braket{p}{p’} \leftrightarrow \delta(p – p’) \)
  4. \( \braket{x}{p’} = \frac{e^{i p x}}{\sqrt{V}} \), where \( V \) is the volume of the box containing the particle. We’ll define the appropriate normalization for an infinite box volume later.

The delta function interpretation of the braket \( \braket{p}{p’} \) justifies the identity operator, since we recover any state in the basis when operating with it. For example, in momentum space

\begin{equation}\label{eqn:qmLecture1:260}
\begin{aligned}
1 \ket{p}
&=
\lr{ \int_{-\infty}^\infty dp’
\ket{p’}\bra{p’} }
\ket{p} \\
&=
\int_{-\infty}^\infty dp’
\ket{p’}
\braket{p’}{p} \\
&=
\int_{-\infty}^\infty dp’
\ket{p’}
\delta(p – p’) \\
&=
\ket{p}.
\end{aligned}
\end{equation}

This also the determination of an integral operator representation for the delta function

\begin{equation}\label{eqn:qmLecture1:500}
\begin{aligned}
\delta(x – x’)
&=
\braket{x}{x’} \\
&=
\int dp \braket{x}{p} \braket{p}{x’} \\
&=
\inv{V} \int dp e^{i p x} e^{-i p x’},
\end{aligned}
\end{equation}

or

\begin{equation}\label{eqn:qmLecture1:520}
\delta(x – x’)
=
\inv{V} \int dp e^{i p (x- x’)}.
\end{equation}

Here we used the fact that \( \braket{p}{x} = \braket{x}{p}^\conj \).

FIXME: do we have a justification for that conjugation with what was defined here so far?

The conversion from a position basis to momentum space is now possible

\begin{equation}\label{eqn:qmLecture1:280}
\begin{aligned}
\braket{p}{\Psi}
&= \Psi(p) \\
&= \int_{-\infty}^\infty \braket{p}{x} \braket{x}{\Psi} dx \\
&= \int_{-\infty}^\infty \frac{e^{-ip x}}{\sqrt{V}} \Psi(x) dx.
\end{aligned}
\end{equation}

The momentum space to position space conversion can be written as

\begin{equation}\label{eqn:qmLecture1:300}
\Psi(x)
= \int_{-\infty}^\infty \frac{e^{ip x}}{\sqrt{V}} \Psi(p) dp.
\end{equation}

Now we can go back and figure out the an expectation

\begin{equation}\label{eqn:qmLecture1:320}
\begin{aligned}
\expectation{p}
&=
\int \Psi^\conj(p) \Psi(p) p d p \\
&=
\int dp
\lr{
\int_{-\infty}^\infty \frac{e^{ip x}}{\sqrt{V}} \Psi^\conj(x) dx
}
\lr{
\int_{-\infty}^\infty \frac{e^{-ip x’}}{\sqrt{V}} \Psi(x’) dx’
}
p \\
&=\int dp dx dx’
\Psi^\conj(x)
\inv{V} e^{ip (x-x’)} \Psi(x’) p \\
&=
\int dp dx dx’
\Psi^\conj(x)
\inv{V} \lr{ -i\PD{x}{e^{ip (x-x’)}} }\Psi(x’) \\
&=
\int dp dx
\Psi^\conj(x) \lr{ -i \PD{x}{} }
\inv{V} \int dx’ e^{ip (x-x’)} \Psi(x’) \\
&=
\int dx
\Psi^\conj(x) \lr{ -i \PD{x}{} }
\int dx’ \lr{ \inv{V} \int dp e^{ip (x-x’)} } \Psi(x’) \\
&=
\int dx
\Psi^\conj(x) \lr{ -i \PD{x}{} }
\int dx’ \delta(x – x’) \Psi(x’) \\
&=
\int dx
\Psi^\conj(x) \lr{ -i \PD{x}{} }
\Psi(x)
\end{aligned}
\end{equation}

Here we’ve essentially calculated the position space representation of the momentum operator, allowing identifications of the following form

\begin{equation}\label{eqn:qmLecture1:380}
p \leftrightarrow -i \PD{x}{}
\end{equation}
\begin{equation}\label{eqn:qmLecture1:400}
p^2 \leftrightarrow – \PDSq{x}{}.
\end{equation}

Alternate starting point.

Most of the above results followed from the claim that \( \braket{x}{p} = e^{i p x} \). Note that this position space representation of the momentum operator can also be taken as the starting point. Given that, the exponential representation of the position-momentum braket follows

\begin{equation}\label{eqn:qmLecture1:540}
\bra{x} P \ket{p}
=
-i \Hbar \PD{x}{} \braket{x}{p},
\end{equation}

but \( \bra{x} P \ket{p} = p \braket{x}{p} \), providing a differential equation for \( \braket{x}{p} \)

\begin{equation}\label{eqn:qmLecture1:560}
p \braket{x}{p} = -i \Hbar \PD{x}{} \braket{x}{p},
\end{equation}

with solution

\begin{equation}\label{eqn:qmLecture1:580}
i p x/\Hbar = \ln \braket{x}{p} + \text{const},
\end{equation}

or
\begin{equation}\label{eqn:qmLecture1:600}
\braket{x}{p} \propto e^{i p x/\Hbar}.
\end{equation}

Matrix interpretation

  1. Ket’s \( \ket{\Psi} \leftrightarrow \text{column vector} \)
  2. Bra’s \( \bra{\Psi} \leftrightarrow {(\text{row vector})}^\conj \)
  3. Operators \( \leftrightarrow \) matrices that act on vectors.

\begin{equation}\label{eqn:qmLecture1:420}
\hat{p} \ket{\Psi} \rightarrow \ket{\Psi’}
\end{equation}

Time evolution

For a state subject to the equations of motion given by the Hamiltonian operator \( \hat{H} \)

\begin{equation}\label{eqn:qmLecture1:440}
i \PD{t}{} \ket{\Psi} = \hat{H} \ket{\Psi},
\end{equation}

the time evolution is given by
\begin{equation}\label{eqn:qmLecture1:460}
\ket{\Psi(t)} = e^{-i \hat{H} t} \ket{\Psi(0)}.
\end{equation}

Incomplete information

We’ll need to introduce the concept of Density matrices. This will bring us to concepts like entanglement.

References

[1] Jun John Sakurai and Jim J Napolitano. Modern quantum mechanics. Pearson Higher Ed, 2014.

Gauge transformed probability current

September 17, 2015 math and physics play , , , ,

[Click here for a PDF of this post with nicer formatting]

Question: Gauge transformed probability current ([1] pr. 2.37 (b))

For the gauge transformed Schrodinger equation

\begin{equation}\label{eqn:gaugeTxCurrent:20}
\inv{2m} \BPi(\Bx) \cdot \BPi(\Bx) \psi(\Bx, t) + e \phi(\Bx) \psi(\Bx, t) = i \Hbar \PD{t}{}\psi(\Bx, t),
\end{equation}

where

\begin{equation}\label{eqn:gaugeTxCurrent:40}
\BPi(\Bx) = -i \Hbar \spacegrad – \frac{e}{c} \BA(\Bx),
\end{equation}

find the probability current defined by

\begin{equation}\label{eqn:gaugeTxCurrent:60}
\PD{t}{\psi} + \spacegrad \cdot \Bj.
\end{equation}

Answer

Equation \ref{eqn:gaugeTxCurrent:20} and its conjugate are

\begin{equation}\label{eqn:gaugeTxCurrent:22}
\begin{aligned}
\inv{2m} \BPi \cdot \BPi \psi + e \phi \psi &= i \Hbar \PD{t}{\psi} \\
\inv{2m} \BPi^\conj \cdot \BPi^\conj \psi^\conj + e \phi \psi^\conj &= -i \Hbar \PD{t}{\psi^\conj}
\end{aligned}
\end{equation}

which can be used immediately in a chain rule expansion of the probability time derivative

\begin{equation}\label{eqn:gaugeTxCurrent:80}
\begin{aligned}
i \Hbar \PD{t}{\rho}
&=
i \Hbar \psi^\conj \PD{t}{\psi} +
i \Hbar \psi \PD{t}{\psi^\conj} \\
&=
\psi^\conj \lr{ \inv{2m} \BPi \cdot \BPi \psi + e \phi \psi } –
\psi \lr{ \inv{2m} \BPi^\conj \cdot \BPi^\conj \psi^\conj + e \phi \psi^\conj
} \\
&=
\inv{2m} \lr{
\psi^\conj \BPi \cdot \BPi \psi
-\psi \BPi^\conj \cdot \BPi^\conj \psi^\conj
}.
\end{aligned}
\end{equation}

We have a difference of conjugates, so can get away with expanding just the first term

\begin{equation}\label{eqn:gaugeTxCurrent:100}
\begin{aligned}
\psi^\conj \BPi \cdot \BPi \psi
&=
\psi^\conj
\psi \\
&=
\psi^\conj
\lr{ -i \Hbar \spacegrad – \frac{e}{c} \BA } \cdot \lr{ -i \Hbar \spacegrad – \frac{e}{c} \BA }
\psi \\
&=
\psi^\conj
\lr{
-\Hbar^2 \spacegrad^2 + \frac{i \Hbar e}{c} \lr{ \BA \cdot \spacegrad + \spacegrad \cdot \BA }
+ \frac{e^2}{c^2} \BA^2
}
\psi.
\end{aligned}
\end{equation}

Note that in the directional derivative terms, the gradient operates on everything to its right, including \( \BA \). Also note that the last term has no imaginary component, so it will not contribute to the difference of conjugates.

This gives

\begin{equation}\label{eqn:gaugeTxCurrent:120}
\begin{aligned}
\psi^\conj \BPi \cdot \BPi \psi – \psi \BPi^\conj \cdot \BPi^\conj \psi^\conj
&=
\psi^\conj
\lr{
-\Hbar^2 \spacegrad^2 \psi + \frac{i \Hbar e}{c} \lr{ \BA \cdot \spacegrad \psi + \spacegrad \cdot (\BA \psi) }
} \\
&\quad –
\psi
\lr{
-\Hbar^2 \spacegrad^2 \psi^\conj – \frac{i \Hbar e}{c} \lr{ \BA \cdot \spacegrad \psi^\conj + \spacegrad \cdot (\BA \psi^\conj) }
} \\
&=
-\Hbar^2 \lr{ \psi^\conj \spacegrad^2 \psi – \psi \spacegrad^2 \psi^\conj } \\
&\quad +
\frac{i \Hbar e}{c}
\lr{
\psi^\conj
\BA \cdot \spacegrad \psi + \psi^\conj \spacegrad \cdot (\BA \psi)
+
\psi
\BA \cdot \spacegrad \psi^\conj + \psi \spacegrad \cdot (\BA \psi^\conj)
}
\end{aligned}
\end{equation}

The first term is recognized as a divergence

\begin{equation}\label{eqn:gaugeTxCurrent:140}
\begin{aligned}
\spacegrad \cdot \lr{ \psi^\conj \spacegrad \psi – \psi \spacegrad \psi^\conj }
&=
\psi^\conj \spacegrad \cdot \spacegrad \psi
+
\spacegrad \psi \cdot \spacegrad \psi^\conj
–
\psi \spacegrad \cdot \spacegrad \psi^\conj
–
\spacegrad \psi^\conj \cdot \spacegrad \psi \\
&= \psi^\conj \spacegrad^2 \psi – \psi \spacegrad^2 \psi^\conj.
\end{aligned}
\end{equation}

The second term can also be factored into a divergence operation

\begin{equation}\label{eqn:gaugeTxCurrent:160}
\begin{aligned}
\psi^\conj
\BA \cdot \spacegrad \psi &+ \psi^\conj \spacegrad \cdot (\BA \psi)
+
\psi
\BA \cdot \spacegrad \psi^\conj + \psi \spacegrad \cdot (\BA \psi^\conj) \\
&=
\lr{ \psi^\conj\BA \cdot \spacegrad \psi
+\psi \spacegrad \cdot (\BA \psi^\conj)
}
+\lr{
\psi \BA \cdot \spacegrad \psi^\conj
+\psi^\conj \spacegrad \cdot (\BA \psi)
} \\
&= 2 \spacegrad \cdot \lr{ \BA \psi \psi^\conj } \\
\end{aligned}
\end{equation}

Putting all the pieces back together we have

\begin{equation}\label{eqn:gaugeTxCurrent:180}
\begin{aligned}
\PD{t}{\rho}
&=
\inv{2m i \Hbar} \lr{
\psi^\conj \BPi \cdot \BPi \psi
-\psi \BPi^\conj \cdot \BPi^\conj \psi^\conj
} \\
&=
\spacegrad \cdot
\inv{2m i \Hbar} \lr{
-\Hbar^2
\lr{ \psi^\conj \spacegrad \psi – \psi \spacegrad \psi^\conj }
+ \frac{ i \Hbar e}{c} 2 \BA \psi \psi^\conj
} \\
&=
\spacegrad \cdot
\lr{
\frac{i \Hbar}{2 m} \lr{ \psi^\conj \spacegrad \psi – \psi \spacegrad \psi^\conj }
+ \frac{e}{m c} \BA \psi \psi^\conj
}.
\end{aligned}
\end{equation}

From \ref{eqn:gaugeTxCurrent:60}, the probability current must be

\begin{equation}\label{eqn:gaugeTxCurrent:200}
\Bj
=
\frac{\Hbar}{2 i m} \lr{ \psi^\conj \spacegrad \psi – \psi \spacegrad \psi^\conj }
– \frac{e}{m c} \BA \psi \psi^\conj,
\end{equation}

or
\begin{equation}\label{eqn:gaugeTxCurrent:220}
\boxed{
\Bj
=
\frac{\Hbar}{m} \textrm{Im} \lr{ \psi^\conj \spacegrad \psi }
– \frac{e}{m c} \BA \psi \psi^\conj.
}
\end{equation}

References

[1] Jun John Sakurai and Jim J Napolitano. Modern quantum mechanics. Pearson Higher Ed, 2014.

Gauge transformation of free particle Hamiltonian

September 15, 2015 math and physics play , , , , , ,


[Click here for a PDF of this post with nicer formatting]

Question:

Given a gauge transformation of the free particle Hamiltonian to

\begin{equation}\label{eqn:gaugeTx:20}
H = \inv{2 m} \BPi \cdot \BPi + e \phi,
\end{equation}

where

\begin{equation}\label{eqn:gaugeTx:40}
\BPi = \Bp – \frac{e}{c} \BA,
\end{equation}

calculate \( m d\Bx/dt \), \( \antisymmetric{\Pi_i}{\Pi_j} \), and \( m d^2\Bx/dt^2 \), where \( \Bx \) is the Heisenberg picture position operator, and the fields are functions only of position \( \phi = \phi(\Bx), \BA = \BA(\Bx) \).

Answer

The final results for these calculations are found in [1], but seem worth deriving to exercise our commutator muscles.

Heisenberg picture velocity operator

The first order of business is the Heisenberg picture velocity operator, but first note

\begin{equation}\label{eqn:gaugeTx:60}
\begin{aligned}
\BPi \cdot \BPi
&= \lr{ \Bp – \frac{e}{c} \BA} \cdot \lr{ \Bp – \frac{e}{c} \BA} \\
&= \Bp^2 – \frac{e}{c} \lr{ \BA \cdot \Bp + \Bp \cdot \BA } + \frac{e^2}{c^2} \BA^2.
\end{aligned}
\end{equation}

The time evolution of the Heisenberg picture position operator is therefore

\begin{equation}\label{eqn:gaugeTx:80}
\begin{aligned}
\ddt{\Bx}
&= \inv{i\Hbar} \antisymmetric{\Bx}{H} \\
&= \inv{i\Hbar 2 m} \antisymmetric{\Bx}{\BPi^2} \\
&= \inv{i\Hbar 2 m} \antisymmetric{\Bx}{\Bp^2 – \frac{e}{c} \lr{ \BA \cdot \Bp
+ \Bp \cdot \BA } + \frac{e^2}{c^2} \BA^2 } \\
&= \inv{i\Hbar 2 m}
\lr{
\antisymmetric{\Bx}{\Bp^2}
– \frac{e}{c} \antisymmetric{\Bx}{ \BA \cdot \Bp + \Bp \cdot \BA }
}
.
\end{aligned}
\end{equation}

For the \( \Bp^2 \) commutator we have

\begin{equation}\label{eqn:gaugeTx:100}
\antisymmetric{x_r}{\Bp^2}
=
i \Hbar \PD{p_r}{\Bp^2}
=
2 i \Hbar p_r,
\end{equation}

or
\begin{equation}\label{eqn:gaugeTx:120}
\antisymmetric{\Bx}{\Bp^2}
=
2 i \Hbar \Bp.
\end{equation}

Computing the remaining commutator, we’ve got

\begin{equation}\label{eqn:gaugeTx:140}
\begin{aligned}
\antisymmetric{x_r}{\Bp \cdot \BA + \BA \cdot \Bp}
&= x_r p_s A_s – p_s A_s x_r \\
&\quad+ x_r A_s p_s – A_s p_s x_r \\
&= \lr{ \antisymmetric{x_r}{p_s} + p_s x_r } A_s – p_s A_s x_r \\
&\quad+ x_r A_s p_s – A_s \lr{ \antisymmetric{p_s}{x_r} + x_r p_s } \\
&= \antisymmetric{x_r}{p_s} A_s + {p_s A_s x_r – p_s A_s x_r} \\
&\quad+ {x_r A_s p_s – x_r A_s p_s} + A_s \antisymmetric{x_r}{p_s} \\
&= 2 i \Hbar \delta_{r s} A_s \\
&= 2 i \Hbar A_r,
\end{aligned}
\end{equation}

so

\begin{equation}\label{eqn:gaugeTx:160}
\antisymmetric{\Bx}{\Bp \cdot \BA + \BA \cdot \Bp} = 2 i \Hbar \BA.
\end{equation}

Assembling these results gives

\begin{equation}\label{eqn:gaugeTx:180}
\boxed{
\ddt{\Bx} = \inv{m} \lr{ \Bp – \frac{e}{c} \BA } = \inv{m} \BPi,
}
\end{equation}

as asserted in the text.

Kinetic Momentum commutators

\begin{equation}\label{eqn:gaugeTx:200}
\begin{aligned}
\antisymmetric{\Pi_r}{\Pi_s}
&=
\antisymmetric{p_r – e A_r/c}{p_s – e A_s/c} \\
&=
{\antisymmetric{p_r}{p_s}}
– \frac{e}{c} \lr{ \antisymmetric{p_r}{A_s} + \antisymmetric{A_r}{p_s}}
+ \frac{e^2}{c^2} {\antisymmetric{A_r}{A_s}} \\
&=
– \frac{e}{c} \lr{ (-i\Hbar) \PD{x_r}{A_s} + (i\Hbar) \PD{x_s}{A_r} } \\
&=
– \frac{i e \Hbar}{c} \lr{ -\PD{x_r}{A_s} + \PD{x_s}{A_r} } \\
&=
– \frac{i e \Hbar}{c} \epsilon_{t s r} B_t,
\end{aligned}
\end{equation}

or
\begin{equation}\label{eqn:gaugeTx:220}
\boxed{
\antisymmetric{\Pi_r}{\Pi_s}
=
\frac{i e \Hbar}{c} \epsilon_{r s t} B_t.
}
\end{equation}

Quantum Lorentz force

For the force equation we have

\begin{equation}\label{eqn:gaugeTx:240}
\begin{aligned}
m \frac{d^2 \Bx}{dt^2}
&= \ddt{\BPi} \\
&= \inv{i \Hbar} \antisymmetric{\BPi}{H} \\
&= \inv{i \Hbar 2 m } \antisymmetric{\BPi}{\BPi^2}
+ \inv{i \Hbar } \antisymmetric{\BPi}{e \phi}.
\end{aligned}
\end{equation}

For the \( \phi \) commutator consider one component

\begin{equation}\label{eqn:gaugeTx:260}
\begin{aligned}
\antisymmetric{\Pi_r}{e \phi}
&=
e \antisymmetric{p_r – \frac{e}{c} A_r}{\phi} \\
&=
e \antisymmetric{p_r}{\phi} \\
&=
e (-i\Hbar) \PD{x_r}{\phi},
\end{aligned}
\end{equation}

or
\begin{equation}\label{eqn:gaugeTx:280}
\inv{i \Hbar} \antisymmetric{\BPi}{e \phi}
=
– e \spacegrad \phi
=
e \BE.
\end{equation}

For the \( \BPi^2 \) commutator I initially did this the hard way (it took four notebook pages, plus two for a false start.) Realizing that I didn’t use \ref{eqn:gaugeTx:220} for that expansion was the clue to doing this more expediently.

Considering a single component

\begin{equation}\label{eqn:gaugeTx:300}
\begin{aligned}
\antisymmetric{\Pi_r}{\BPi^2}
&=
\antisymmetric{\Pi_r}{\Pi_s \Pi_s} \\
&=
\Pi_r \Pi_s \Pi_s – \Pi_s \Pi_s \Pi_r \\
&=
\lr{ \antisymmetric{\Pi_r}{\Pi_s} + {\Pi_s \Pi_r} }
\Pi_s
– \Pi_s
\lr{ \antisymmetric{\Pi_s}{\Pi_r} + {\Pi_r \Pi_s} } \\
&= i \Hbar \frac{e}{c} \epsilon_{r s t}
\lr{ B_t \Pi_s + \Pi_s B_t },
\end{aligned}
\end{equation}

or

\begin{equation}\label{eqn:gaugeTx:320}
\begin{aligned}
\inv{ i \Hbar 2 m} \antisymmetric{\BPi}{\BPi^2}
&= \frac{e}{2 m c } \epsilon_{r s t} \Be_r
\lr{ B_t \Pi_s + \Pi_s B_t } \\
&= \frac{e}{ 2 m c }
\lr{
\BPi \cross \BB
– \BB \cross \BPi
}.
\end{aligned}
\end{equation}

Putting all the pieces together we’ve got the quantum equivalent of the Lorentz force equation

\begin{equation}\label{eqn:gaugeTx:340}
\boxed{
m \frac{d^2 \Bx}{dt^2} = e \BE + \frac{e}{2 c} \lr{
\frac{d\Bx}{dt} \cross \BB
– \BB \cross \frac{d\Bx}{dt}
}.
}
\end{equation}

While this looks equivalent to the classical result, all the vectors here are Heisenberg picture operators dependent on position.

References

[1] Jun John Sakurai and Jim J Napolitano. Modern quantum mechanics. Pearson Higher Ed, 2014.

Notes from reading Sakurai’s QM, prep for graduate Quantum Mechanics (PHY1520H)

September 14, 2015 phy1520 , ,

I start the lectures for grad Quantum Mechanics (PHY1520H) tomorrow, taught by Prof. Arun Paramekanti.

The curriculum has scary range of material, although some of it will be review. As prep for that, I’ve been working some problems, and have also explored a few interesting details and digressions. My notes for that preparatory work, collecting a number of individual posts from over the summer, is now posted.