Dirac equation

Alternate Dirac equation representation

November 27, 2015 phy1520 No comments , ,

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

Given an alternate representation of the Dirac equation

\begin{equation}\label{eqn:diracAlternate:20}
H =
\begin{bmatrix}
m c^2 + V_0 & c \hat{p} \\
c \hat{p} & – m c^2 + V_0
\end{bmatrix},
\end{equation}

calculate the constant momentum solutions, the Heisenberg velocity operator \( \hat{v} \), and find the form of the probability density current.

Plane wave solutions

The action of the Hamiltonian on

\begin{equation}\label{eqn:diracAlternate:40}
\psi =
e^{i k x – i E t/\Hbar}
\begin{bmatrix}
\psi_1 \\
\psi_2
\end{bmatrix}
\end{equation}

is
\begin{equation}\label{eqn:diracAlternate:60}
\begin{aligned}
H \psi
&=
\begin{bmatrix}
m c^2 + V_0 & c (-i \Hbar) i k \\
c (-i \Hbar) i k & – m c^2 + V_0
\end{bmatrix}
\begin{bmatrix}
\psi_1 \\
\psi_2
\end{bmatrix}
e^{i k x – i E t/\Hbar} \\
&=
\begin{bmatrix}
m c^2 + V_0 & c \Hbar k \\
c \Hbar k & – m c^2 + V_0
\end{bmatrix}
\psi.
\end{aligned}
\end{equation}

Writing

\begin{equation}\label{eqn:diracAlternate:80}
H_k
=
\begin{bmatrix}
m c^2 + V_0 & c \Hbar k \\
c \Hbar k & – m c^2 + V_0
\end{bmatrix}
\end{equation}

the characteristic equation is

\begin{equation}\label{eqn:diracAlternate:100}
0
=
(m c^2 + V_0 – \lambda)
(-m c^2 + V_0 – \lambda)
– (c \Hbar k)^2
=
\lr{ (\lambda – V_0)^2 – (m c^2)^2 } – (c \Hbar k)^2,
\end{equation}

so

\begin{equation}\label{eqn:diracAlternate:120}
\lambda = V_0 \pm \epsilon,
\end{equation}

where
\begin{equation}\label{eqn:diracAlternate:140}
\epsilon^2 = (m c^2)^2 + (c \Hbar k)^2.
\end{equation}

We’ve got

\begin{equation}\label{eqn:diracAlternate:160}
\begin{aligned}
H – ( V_0 + \epsilon )
&=
\begin{bmatrix}
m c^2 – \epsilon & c \Hbar k \\
c \Hbar k & – m c^2 – \epsilon
\end{bmatrix} \\
H – ( V_0 – \epsilon )
&=
\begin{bmatrix}
m c^2 + \epsilon & c \Hbar k \\
c \Hbar k & – m c^2 + \epsilon
\end{bmatrix},
\end{aligned}
\end{equation}

so the eigenkets are

\begin{equation}\label{eqn:diracAlternate:180}
\begin{aligned}
\ket{V_0+\epsilon}
&\propto
\begin{bmatrix}
-c \Hbar k \\
m c^2 – \epsilon
\end{bmatrix} \\
\ket{V_0-\epsilon}
&\propto
\begin{bmatrix}
-c \Hbar k \\
m c^2 + \epsilon
\end{bmatrix}.
\end{aligned}
\end{equation}

Up to an arbitrary phase for each, these are

\begin{equation}\label{eqn:diracAlternate:200}
\begin{aligned}
\ket{V_0 + \epsilon}
&=
\inv{\sqrt{ 2 \epsilon ( \epsilon – m c^2) }}
\begin{bmatrix}
c \Hbar k \\
\epsilon -m c^2
\end{bmatrix} \\
\ket{V_0 – \epsilon}
&=
\inv{\sqrt{ 2 \epsilon ( \epsilon + m c^2) }}
\begin{bmatrix}
-c \Hbar k \\
\epsilon + m c^2
\end{bmatrix} \\
\end{aligned}
\end{equation}

We can now write

\begin{equation}\label{eqn:diracAlternate:220}
H_k =
E
\begin{bmatrix}
V_0 + \epsilon & 0 \\
0 & V_0 – \epsilon
\end{bmatrix}
E^{-1},
\end{equation}

where
\begin{equation}\label{eqn:diracAlternate:240}
\begin{aligned}
E &=
\inv{\sqrt{2 \epsilon} }
\begin{bmatrix}
\frac{c \Hbar k}{ \sqrt{ \epsilon – m c^2 } } & -\frac{c \Hbar k}{ \sqrt{ \epsilon + m c^2 } } \\
\sqrt{ \epsilon – m c^2 } & \sqrt{ \epsilon + m c^2 }
\end{bmatrix}, \qquad k > 0 \\
E &=
\inv{\sqrt{2 \epsilon} }
\begin{bmatrix}
-\frac{c \Hbar k}{ \sqrt{ \epsilon – m c^2 } } & -\frac{c \Hbar k}{ \sqrt{ \epsilon + m c^2 } } \\
-\sqrt{ \epsilon – m c^2 } & \sqrt{ \epsilon + m c^2 }
\end{bmatrix}, \qquad k < 0. \end{aligned} \end{equation} Here the signs have been adjusted to ensure the transformation matrix has a unit determinant. Observe that there's redundancy in this matrix since \( \ifrac{c \Hbar \Abs{k}}{ \sqrt{ \epsilon - m c^2 } } = \sqrt{ \epsilon + m c^2 } \), and \( \ifrac{c \Hbar \Abs{k}}{ \sqrt{ \epsilon + m c^2 } } = \sqrt{ \epsilon - m c^2 } \), which allows the transformation matrix to be written in the form of a rotation matrix \begin{equation}\label{eqn:diracAlternate:260} \begin{aligned} E &= \inv{\sqrt{2 \epsilon} } \begin{bmatrix} \frac{c \Hbar k}{ \sqrt{ \epsilon - m c^2 } } & -\frac{c \Hbar k}{ \sqrt{ \epsilon + m c^2 } } \\ \frac{c \Hbar k}{ \sqrt{ \epsilon + m c^2 } } & \frac{c \Hbar k}{ \sqrt{ \epsilon - m c^2 } } \end{bmatrix}, \qquad k > 0 \\
E &=
\inv{\sqrt{2 \epsilon} }
\begin{bmatrix}
-\frac{c \Hbar k}{ \sqrt{ \epsilon – m c^2 } } & -\frac{c \Hbar k}{ \sqrt{ \epsilon + m c^2 } } \\
\frac{c \Hbar k}{ \sqrt{ \epsilon + m c^2 } } & -\frac{c \Hbar k}{ \sqrt{ \epsilon – m c^2 } }
\end{bmatrix}, \qquad k < 0 \\ \end{aligned} \end{equation} With \begin{equation}\label{eqn:diracAlternate:280} \begin{aligned} \cos\theta &= \frac{c \Hbar \Abs{k}}{ \sqrt{ 2 \epsilon( \epsilon - m c^2) } } = \frac{\sqrt{ \epsilon + m c^2} }{ \sqrt{ 2 \epsilon}}\\ \sin\theta &= \frac{c \Hbar k}{ \sqrt{ 2 \epsilon( \epsilon + m c^2) } } = \frac{\textrm{sgn}(k) \sqrt{ \epsilon - m c^2}}{ \sqrt{ 2 \epsilon } }, \end{aligned} \end{equation} the transformation matrix (and eigenkets) is \begin{equation}\label{eqn:diracAlternate:300} \boxed{ E = \begin{bmatrix} \ket{V_0 + \epsilon} & \ket{V_0 - \epsilon} \end{bmatrix} = \begin{bmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{bmatrix}. } \end{equation} Observe that \ref{eqn:diracAlternate:280} can be simplified by using double angle formulas \begin{equation}\label{eqn:diracAlternate:320} \begin{aligned} \cos(2 \theta) &= \frac{\lr{ \epsilon + m c^2} }{ 2 \epsilon } - \frac{\lr{ \epsilon - m c^2}}{ 2 \epsilon } \\ &= \frac{1}{ 2 \epsilon } \lr{ \epsilon + m c^2 - \epsilon + m c^2 } \\ &= \frac{m c^2 }{ \epsilon }, \end{aligned} \end{equation} and \begin{equation}\label{eqn:diracAlternate:340} \sin(2\theta) = 2 \frac{1}{2 \epsilon} \textrm{sgn}(k ) \sqrt{ \epsilon^2 - (m c^2)^2 } = \frac{\Hbar k c}{\epsilon}. \end{equation} This allows all the \( \theta \) dependence on \( \Hbar k c \) and \( m c^2 \) to be expressed as a ratio of momenta \begin{equation}\label{eqn:diracAlternate:360} \boxed{ \tan(2\theta) = \frac{\Hbar k}{m c}. } \end{equation}

Hyperbolic solutions

For a wave function of the form

\begin{equation}\label{eqn:diracAlternate:380}
\psi =
e^{k x – i E t/\Hbar}
\begin{bmatrix}
\psi_1 \\
\psi_2
\end{bmatrix},
\end{equation}

some of the work above can be recycled if we substitute \( k \rightarrow -i k \), which yields unnormalized eigenfunctions

\begin{equation}\label{eqn:diracAlternate:400}
\begin{aligned}
\ket{V_0+\epsilon}
&\propto
\begin{bmatrix}
i c \Hbar k \\
m c^2 – \epsilon
\end{bmatrix} \\
\ket{V_0-\epsilon}
&\propto
\begin{bmatrix}
i c \Hbar k \\
m c^2 + \epsilon
\end{bmatrix},
\end{aligned}
\end{equation}

where

\begin{equation}\label{eqn:diracAlternate:420}
\epsilon^2 = (m c^2)^2 – (c \Hbar k)^2.
\end{equation}

The squared magnitude of these wavefunctions are

\begin{equation}\label{eqn:diracAlternate:440}
\begin{aligned}
(c \Hbar k)^2 + (m c^2 \mp \epsilon)^2
&=
(c \Hbar k)^2 + (m c^2)^2 + \epsilon^2 \mp 2 m c^2 \epsilon \\
&=
(c \Hbar k)^2 + (m c^2)^2 + (m c^2)^2 \mp (c \Hbar k)^2 – 2 m c^2 \epsilon \\
&= 2 (m c^2)^2 \mp 2 m c^2 \epsilon \\
&= 2 m c^2 ( m c^2 \mp \epsilon ),
\end{aligned}
\end{equation}

so, up to a constant phase for each, the normalized kets are

\begin{equation}\label{eqn:diracAlternate:460}
\begin{aligned}
\ket{V_0+\epsilon}
&=
\inv{\sqrt{ 2 m c^2 ( m c^2 – \epsilon ) }}
\begin{bmatrix}
i c \Hbar k \\
m c^2 – \epsilon
\end{bmatrix} \\
\ket{V_0-\epsilon}
&=
\inv{\sqrt{ 2 m c^2 ( m c^2 + \epsilon ) }}
\begin{bmatrix}
i c \Hbar k \\
m c^2 + \epsilon
\end{bmatrix},
\end{aligned}
\end{equation}

After the \( k \rightarrow -i k \) substitution, \( H_k \) is not Hermitian, so these kets aren’t expected to be orthonormal, which is readily verified

\begin{equation}\label{eqn:diracAlternate:480}
\begin{aligned}
\braket{V_0+\epsilon}{V_0-\epsilon}
&=
\inv{\sqrt{ 2 m c^2 ( m c^2 – \epsilon ) }}
\inv{\sqrt{ 2 m c^2 ( m c^2 + \epsilon ) }}
\begin{bmatrix}
-i c \Hbar k &
m c^2 – \epsilon
\end{bmatrix}
\begin{bmatrix}
i c \Hbar k \\
m c^2 + \epsilon
\end{bmatrix} \\
&=
\frac{ 2 ( c \Hbar k )^2 }{2 m c^2 \sqrt{(\Hbar k c)^2} } \\
&=
\textrm{sgn}(k)
\frac{
\Hbar k }{m c } .
\end{aligned}
\end{equation}

Heisenberg velocity operator

\begin{equation}\label{eqn:diracAlternate:500}
\begin{aligned}
\hat{v}
&= \inv{i \Hbar} \antisymmetric{ \hat{x} }{ H} \\
&= \inv{i \Hbar} \antisymmetric{ \hat{x} }{ m c^2 \sigma_z + V_0 + c \hat{p} \sigma_x } \\
&= \frac{c \sigma_x}{i \Hbar} \antisymmetric{ \hat{x} }{ \hat{p} } \\
&= c \sigma_x.
\end{aligned}
\end{equation}

Probability current

Acting against a completely general wavefunction the Hamiltonian action \( H \psi \) is

\begin{equation}\label{eqn:diracAlternate:520}
\begin{aligned}
i \Hbar \PD{t}{\psi}
&= m c^2 \sigma_z \psi + V_0 \psi + c \hat{p} \sigma_x \psi \\
&= m c^2 \sigma_z \psi + V_0 \psi -i \Hbar c \sigma_x \PD{x}{\psi}.
\end{aligned}
\end{equation}

Conversely, the conjugate \( (H \psi)^\dagger \) is

\begin{equation}\label{eqn:diracAlternate:540}
-i \Hbar \PD{t}{\psi^\dagger}
= m c^2 \psi^\dagger \sigma_z + V_0 \psi^\dagger +i \Hbar c \PD{x}{\psi^\dagger} \sigma_x.
\end{equation}

These give

\begin{equation}\label{eqn:diracAlternate:560}
\begin{aligned}
i \Hbar \psi^\dagger \PD{t}{\psi}
&=
m c^2 \psi^\dagger \sigma_z \psi + V_0 \psi^\dagger \psi -i \Hbar c \psi^\dagger \sigma_x \PD{x}{\psi} \\
-i \Hbar \PD{t}{\psi^\dagger} \psi
&= m c^2 \psi^\dagger \sigma_z \psi + V_0 \psi^\dagger \psi +i \Hbar c \PD{x}{\psi^\dagger} \sigma_x \psi.
\end{aligned}
\end{equation}

Taking differences
\begin{equation}\label{eqn:diracAlternate:580}
\psi^\dagger \PD{t}{\psi} + \PD{t}{\psi^\dagger} \psi
=
– c \psi^\dagger \sigma_x \PD{x}{\psi} – c \PD{x}{\psi^\dagger} \sigma_x \psi,
\end{equation}

or

\begin{equation}\label{eqn:diracAlternate:600}
0
=
\PD{t}{}
\lr{
\psi^\dagger \psi
}
+
\PD{x}{}
\lr{
c \psi^\dagger \sigma_x \psi
}.
\end{equation}

The probability current still has the usual form \( \rho = \psi^\dagger \psi = \psi_1^\conj \psi_1 + \psi_2^\conj \psi_2 \), but the probability current with this representation of the Dirac Hamiltonian is

\begin{equation}\label{eqn:diracAlternate:620}
\begin{aligned}
j
&= c \psi^\dagger \sigma_x \psi \\
&= c
\begin{bmatrix}
\psi_1^\conj &
\psi_2^\conj
\end{bmatrix}
\begin{bmatrix}
\psi_2 \\
\psi_1
\end{bmatrix} \\
&= c \lr{ \psi_1^\conj \psi_2 + \psi_2^\conj \psi_1 }.
\end{aligned}
\end{equation}

PHY1520H Graduate Quantum Mechanics. Lecture 10: 1D Dirac scattering off potential step. Taught by Prof. Arun Paramekanti

October 20, 2015 phy1520 No comments , ,

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

Dirac scattering off a potential step

For the non-relativistic case we have

\begin{equation}\label{eqn:qmLecture10:20}
\begin{aligned}
E < V_0 &\Rightarrow T = 0, R = 1 \\ E > V_0 &\Rightarrow T > 0, R < 1.
\end{aligned}
\end{equation}

What happens for a relativistic 1D particle?

Referring to fig. 1.

fig. 1. Potential step

fig. 1. Potential step

the region I Hamiltonian is

\begin{equation}\label{eqn:qmLecture10:40}
H =
\begin{bmatrix}
\hat{p} c & m c^2 \\
m c^2 & – \hat{p} c
\end{bmatrix},
\end{equation}

for which the solution is

\begin{equation}\label{eqn:qmLecture10:60}
\Phi = e^{i k_1 x }
\begin{bmatrix}
\cos \theta_1 \\
\sin \theta_1
\end{bmatrix},
\end{equation}

where
\begin{equation}\label{eqn:qmLecture10:80}
\begin{aligned}
\cos 2 \theta_1 &= \frac{ \Hbar c k_1 }{E_{k_1}} \\
\sin 2 \theta_1 &= \frac{ m c^2 }{E_{k_1}} \\
\end{aligned}
\end{equation}

To consider the \( k_1 < 0 \) case, note that

\begin{equation}\label{eqn:qmLecture10:100}
\begin{aligned}
\cos^2 \theta_1 – \sin^2 \theta_1 &= \cos 2 \theta_1 \\
2 \sin\theta_1 \cos\theta_1 &= \sin 2 \theta_1
\end{aligned}
\end{equation}

so after flipping the signs on all the \( k_1 \) terms we find for the reflected wave

\begin{equation}\label{eqn:qmLecture10:120}
\Phi = e^{-i k_1 x}
\begin{bmatrix}
\sin\theta_1 \\
\cos\theta_1
\end{bmatrix}.
\end{equation}

FIXME: this reasoning doesn’t entirely make sense to me. Make sense of this by trying this solution as was done for the form of the incident wave solution.

The region I wave has the form

\begin{equation}\label{eqn:qmLecture10:140}
\Phi_I
=
A e^{i k_1 x}
\begin{bmatrix}
\cos\theta_1 \\
\sin\theta_1 \\
\end{bmatrix}
+
B e^{-i k_1 x}
\begin{bmatrix}
\sin\theta_1 \\
\cos\theta_1 \\
\end{bmatrix}.
\end{equation}

By the time we are done we want to have computed the reflection coefficient

\begin{equation}\label{eqn:qmLecture10:160}
R =
\frac{\Abs{B}^2}{\Abs{A}^2}.
\end{equation}

The region I energy is

\begin{equation}\label{eqn:qmLecture10:180}
E = \sqrt{ \lr{ m c^2}^2 + \lr{ \Hbar c k_1 }^2 }.
\end{equation}

We must have
\begin{equation}\label{eqn:qmLecture10:200}
E
=
\sqrt{ \lr{ m c^2}^2 + \lr{ \Hbar c k_2 }^2 } + V_0
=
\sqrt{ \lr{ m c^2}^2 + \lr{ \Hbar c k_1 }^2 },
\end{equation}

so

\begin{equation}\label{eqn:qmLecture10:220}
\begin{aligned}
\lr{ \Hbar c k_2 }^2
&=
\lr{ E – V_0 }^2 – \lr{ m c^2}^2 \\
&=
\underbrace{\lr{ E – V_0 + m c }}_{r_1}\underbrace{\lr{ E – V_0 – m c }}_{r_2}.
\end{aligned}
\end{equation}

The \( r_1 \) and \( r_2 \) branches are sketched in fig. 2.

fig. 2. Energy signs

fig. 2. Energy signs

For low energies, we have a set of potentials for which we will have propagation, despite having a potential barrier. For still higher values of the potential barrier the product \( r_1 r_2 \) will be negative, so the solutions will be decaying. Finally, for even higher energies, there will again be propagation.

The non-relativistic case is sketched in fig. 3.

fig. 3. Effects of increasing potential for non-relativistic case

fig. 3. Effects of increasing potential for non-relativistic case

For the relativistic case we must consider three different cases, sketched in fig 4, fig 5, and fig 6 respectively. For the low potential energy, a particle with positive group velocity (what we’ve called right moving) can be matched to an equal energy portion of the potential shifted parabola in region II. This is a case where we have transmission, but no antiparticle creation. There will be an energy region where the region II wave function has only a dissipative term, since there is no region of either of the region II parabolic branches available at the incident energy. When the potential is shifted still higher so that \( V_0 > E + m c^2 \), a positive group velocity in region I with a given energy can be matched to an antiparticle branch in the region II parabolic energy curve.

lecture10Fig4a

Fig 4. Low potential energy

lecture10Fig4b

fig. 5. High enough potential energy for no propagation

lecture10Fig4c

fig 6. High potential energy

 

Boundary value conditions

We want to ensure that the current across the barrier is conserved (no particles are lost), as sketched in fig. 7.

 

fig. 7. Transmitted, reflected and incident components.

fig. 7. Transmitted, reflected and incident components.

Recall that given a wave function

\begin{equation}\label{eqn:qmLecture10:240}
\Psi =
\begin{bmatrix}
\psi_1 \\
\psi_2
\end{bmatrix},
\end{equation}

the density and currents are respectively

\begin{equation}\label{eqn:qmLecture10:260}
\begin{aligned}
\rho &= \psi_1^\conj \psi_1 + \psi_2^\conj \psi_2 \\
j &= \psi_1^\conj \psi_1 – \psi_2^\conj \psi_2
\end{aligned}
\end{equation}

Matching boundary value conditions requires

  1. For both the relativistic and non-relativistic cases we must have\begin{equation}\label{eqn:qmLecture10:280}
    \Psi_{\textrm{L}} = \Psi_{\textrm{R}}, \qquad \mbox{at \( x = 0 \).}
    \end{equation}
  2. For the non-relativistic case we want
    \begin{equation}\label{eqn:qmLecture10:300}
    \int_{-\epsilon}^\epsilon -\frac{\Hbar^2}{2m} \PDSq{x}{\Psi} =
    {\int_{-\epsilon}^\epsilon \lr{ E – V(x) } \Psi(x)}.
    \end{equation}The RHS integral is zero, so

    \begin{equation}\label{eqn:qmLecture10:320}
    -\frac{\Hbar^2}{2m} \lr{ \evalbar{\PD{x}{\Psi}}{{\textrm{R}}} – \evalbar{\PD{x}{\Psi}}{{\textrm{L}}} } = 0.
    \end{equation}

    We have to match

    For the relativistic case

    \begin{equation}\label{eqn:qmLecture10:460}
    -i \Hbar \sigma_z \int_{-\epsilon}^\epsilon \PD{x}{\Psi} +
    {m c^2 \sigma_x \int_{-\epsilon}^\epsilon \psi}
    =
    {\int_{-\epsilon}^\epsilon \lr{ E – V_0 } \psi},
    \end{equation}

the second two integrals are wiped out, so

\begin{equation}\label{eqn:qmLecture10:340}
-i \Hbar c \sigma_z \lr{ \psi(\epsilon) – \psi(-\epsilon) }
=
-i \Hbar c \sigma_z \lr{ \psi_{\textrm{R}} – \psi_{\textrm{L}} }.
\end{equation}

so we must match

\begin{equation}\label{eqn:qmLecture10:360}
\sigma_z \psi_{\textrm{R}} = \sigma_z \psi_{\textrm{L}} .
\end{equation}

It appears that things are simpler, because we only have to match the wave function values at the boundary, and don’t have to match the derivatives too. However, we have a two component wave function, so there are still two tasks.

Solving the system

Let’s look for a solution for the \( E + m c^2 > V_0 \) case on the right branch, as sketched in fig. 8.

 

fig. 8. High potential region. Anti-particle transmission.

fig. 8. High potential region. Anti-particle transmission.

While the right branch in this case is left going, this might work out since that is an antiparticle. We could try both.

Try

\begin{equation}\label{eqn:qmLecture10:480}
\Psi_{II} = D e^{i k_2 x}
\begin{bmatrix}
-\sin\theta_2 \\
\cos\theta_2
\end{bmatrix}.
\end{equation}

This is justified by

\begin{equation}\label{eqn:qmLecture10:500}
+E \rightarrow
\begin{bmatrix}
\cos\theta \\
\sin\theta
\end{bmatrix},
\end{equation}

so

\begin{equation}\label{eqn:qmLecture10:520}
-E \rightarrow
\begin{bmatrix}
-\sin\theta \\
\cos\theta \\
\end{bmatrix}
\end{equation}

At \( x = 0 \) the exponentials vanish, so equating the waves at that point means

\begin{equation}\label{eqn:qmLecture10:380}
\begin{bmatrix}
\cos\theta_1 \\
\sin\theta_1 \\
\end{bmatrix}
+
\frac{B}{A}
\begin{bmatrix}
\sin\theta_1 \\
\cos\theta_1 \\
\end{bmatrix}
=
\frac{D}{A}
\begin{bmatrix}
-\sin\theta_2 \\
\cos\theta_2
\end{bmatrix}.
\end{equation}

Solving this yields

\begin{equation}\label{eqn:qmLecture10:400}
\frac{B}{A} = – \frac{\cos(\theta_1 – \theta_2)}{\sin(\theta_1 + \theta_2)}.
\end{equation}

This yields

\begin{equation}\label{eqn:qmLecture10:420}
\boxed{
R = \frac{1 + \cos( 2 \theta_1 – 2 \theta_2) }{1 – \cos( 2 \theta_1 – 2 \theta_2)}.
}
\end{equation}

As \( V_0 \rightarrow \infty \) this simplifies to

\begin{equation}\label{eqn:qmLecture10:440}
R = \frac{ E – \sqrt{ E^2 – \lr{ m c^2 }^2 } }{ E + \sqrt{ E^2 – \lr{ m c^2 }^2 } }.
\end{equation}

Filling in the details for these results part of problem set 4.

Second update of aggregate notes for phy1520, Graduate Quantum Mechanics

October 20, 2015 phy1520 No comments , , , , , , , , , , , ,

I’ve posted a second update of my aggregate notes for PHY1520H Graduate Quantum Mechanics, taught by Prof. Arun Paramekanti. In addition to what was noted previously, this contains lecture notes up to lecture 9, my ungraded solutions for the second problem set, and some additional worked practise problems.

Most of the content was posted individually in the following locations, but those original documents will not be maintained individually any further.

PHY1520H Graduate Quantum Mechanics. Lecture 9: Dirac equation (cont.). Taught by Prof. Arun Paramekanti

October 15, 2015 phy1520 No comments , , , , , , , ,

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

Where we left off

\begin{equation}\label{eqn:qmLecture9:20}
-i \Hbar \PD{t}{}
\begin{bmatrix}
\psi_1 \\
\psi_2
\end{bmatrix}
=
\begin{bmatrix}
-i \Hbar c \PD{x}{} & m c^2 \\
m c^2 & i \Hbar c \PD{x}{} \\
\end{bmatrix}.
\end{equation}

With a potential this would be

\begin{equation}\label{eqn:qmLecture9:40}
-i \Hbar \PD{t}{}
\begin{bmatrix}
\psi_1 \\
\psi_2
\end{bmatrix}
=
\begin{bmatrix}
-i \Hbar c \PD{x}{} + V(x) & m c^2 \\
m c^2 & i \Hbar c \PD{x}{} + V(x) \\
\end{bmatrix}.
\end{equation}

This means that the potential is raising the energy eigenvalue of the system.

Free Particle

Assuming a form

\begin{equation}\label{eqn:qmLecture9:60}
\begin{bmatrix}
\psi_1(x,t) \\
\psi_2(x,t)
\end{bmatrix}
=
e^{i k x}
\begin{bmatrix}
f_1(t) \\
f_2(t) \\
\end{bmatrix},
\end{equation}

and plugging back into the Dirac equation we have

\begin{equation}\label{eqn:qmLecture9:80}
-i \Hbar \PD{t}{}
\begin{bmatrix}
f_1 \\
f_2
\end{bmatrix}
=
\begin{bmatrix}
k \Hbar c & m c^2 \\
m c^2 & – \Hbar k c \\
\end{bmatrix}
\begin{bmatrix}
f_1 \\
f_2
\end{bmatrix}.
\end{equation}

We can use a diagonalizing rotation

\begin{equation}\label{eqn:qmLecture9:100}
\begin{bmatrix}
f_1 \\
f_2
\end{bmatrix}
=
\begin{bmatrix}
\cos\theta_k & -\sin\theta_k \\
\sin\theta_k & \cos\theta_k \\
\end{bmatrix}
\begin{bmatrix}
f_{+} \\
f_{-} \\
\end{bmatrix}.
\end{equation}

Plugging this in reduces the system to the form

\begin{equation}\label{eqn:qmLecture9:140}
-i \Hbar \PD{t}{}
\begin{bmatrix}
f_{+} \\
f_{-} \\
\end{bmatrix}
=
\begin{bmatrix}
E_k & 0 \\
0 & -E_k
\end{bmatrix}
\begin{bmatrix}
f_{+} \\
f_{-} \\
\end{bmatrix}.
\end{equation}

Where the rotation angle is found to be given by

\begin{equation}\label{eqn:qmLecture9:160}
\begin{aligned}
\sin(2 \theta_k) &= \frac{m c^2}{\sqrt{(\Hbar k c)^2 + m^2 c^4}} \\
\cos(2 \theta_k) &= \frac{\Hbar k c}{\sqrt{(\Hbar k c)^2 + m^2 c^4}} \\
E_k &= \sqrt{(\Hbar k c)^2 + m^2 c^4}
\end{aligned}
\end{equation}

See fig. 1 for a sketch of energy vs momentum. The asymptotes are the limiting cases when \( m c^2 \rightarrow 0 \). The \( + \) branch is what we usually associate with particles. What about the other energy states. For Fermions Dirac argued that the lower energy states could be thought of as “filled up”, using the Pauli principle to leave only the positive energy states available. This was called the “Dirac Sea”. This isn’t a good solution, and won’t work for example for Bosons.

fig. 1. Dirac equation solution space

fig. 1. Dirac equation solution space

Another way to rationalize this is to employ ideas from solid state theory. For example consider a semiconductor with a valence and conduction band as sketched in fig. 2.

fig. 2. Solid state valence and conduction band transition

fig. 2. Solid state valence and conduction band transition

A photon can excite an electron from the valence band to the conduction band, leaving all the valence band states filled except for one (a hole). For an electron we can use almost the same picture, as sketched in fig. 3.

fig. 3. Pair creation

fig. 3. Pair creation

A photon with energy \( E_k – (-E_k) \) can create a positron-electron pair from the vacuum, where the energy of the electron and positron pair is \( E_k \).

At high enough energies, we can see this pair creation occur.

Zitterbewegung

If a particle is created at a non-eigenstate such as one on the asymptotes, then oscillations between the positive and negative branches are possible as sketched in fig. 4.

fig. 4. Zitterbewegung oscillation

fig. 4. Zitterbewegung oscillation

Only “vertical” oscillations between the positive and negative locations on these branches is possible since those are the points that match the particle momentum. Examining this will be the aim of one of the problem set problems.

Probability and current density

If we define a probability density

\begin{equation}\label{eqn:qmLecture9:180}
\rho(x, t) = \Abs{\psi_1}^2 + \Abs{\psi_2}^2,
\end{equation}

does this satisfy a probability conservation relation

\begin{equation}\label{eqn:qmLecture9:200}
\PD{t}{\rho} + \PD{x}{j} = 0,
\end{equation}

where \( j \) is the probability current. Plugging in the density, we have

\begin{equation}\label{eqn:qmLecture9:220}
\PD{t}{\rho}
=
\PD{t}{\psi_1^\conj} \psi_1
+
\psi_1^\conj \PD{t}{\psi_1}
+
\PD{t}{\psi_2^\conj} \psi_2
+
\psi_2^\conj \PD{t}{\psi_2}.
\end{equation}

It turns out that the probability current has the form

\begin{equation}\label{eqn:qmLecture9:240}
j(x,t) = c \lr{ \psi_1^\conj \psi_1 + \psi_2^\conj \psi_2 }.
\end{equation}

Here the speed of light \( c \) is the slope of the line in the plots above. We can think of this current density as right movers minus the left movers. Any state that is given can be thought of as a combination of right moving and left moving states, neither of which are eigenstates of the free particle Hamiltonian.

Potential step

The next logical thing to think about, as in non-relativistic quantum mechanics, is to think about what occurs when the particle hits a potential step, as in fig. 5.

fig. 5. Reflection off a potential barrier

fig. 5. Reflection off a potential barrier

The approach is the same. We write down the wave functions for the \( V = 0 \) region (I), and the higher potential region (II).

The eigenstates are found on the solid lines above the asymptotes on the right hand movers side as sketched in fig. 6. The right and left moving designations are based on the phase velocity \( \PDi{k}{E} \) (approaching \( \pm c \) on the top-right and top-left quadrants respectively).

fig. 6. Right movers and left movers

fig. 6. Right movers and left movers

For \( k > 0 \), an eigenstate for the incident wave is

\begin{equation}\label{eqn:qmLecture9:261}
\Bpsi_{\textrm{inc}}(x) =
\begin{bmatrix}
\cos\theta_k \\
\sin\theta_k
\end{bmatrix}
e^{i k x},
\end{equation}

For the reflected wave function, we pick a function on the left moving side of the positive energy branch.

\begin{equation}\label{eqn:qmLecture9:260}
\Bpsi_{\textrm{ref}}(x) =
\begin{bmatrix}
? \\
?
\end{bmatrix}
e^{-i k x},
\end{equation}

We’ll go through this in more detail next time.

Question: Calculate the right going diagonalization

Prove (7).

Answer

To determine the relations for \( \theta_k \) we have to solve

\begin{equation}\label{eqn:qmLecture9:280}
\begin{bmatrix}
E_k & 0 \\
0 & -E_k
\end{bmatrix}
= R^{-1} H R.
\end{equation}

Working with \( \Hbar = c = 1 \) temporarily, and \( C = \cos\theta_k, S = \sin\theta_k \), that is

\begin{equation}\label{eqn:qmLecture9:300}
\begin{aligned}
\begin{bmatrix}
E_k & 0 \\
0 & -E_k
\end{bmatrix}
&=
\begin{bmatrix}
C & S \\
-S & C
\end{bmatrix}
\begin{bmatrix}
k & m \\
m & -k
\end{bmatrix}
\begin{bmatrix}
C & -S \\
S & C
\end{bmatrix} \\
&=
\begin{bmatrix}
C & S \\
-S & C
\end{bmatrix}
\begin{bmatrix}
k C + m S & -k S + m C \\
m C – k S & -m S – k C
\end{bmatrix} \\
&=
\begin{bmatrix}
k C^2 + m S C + m C S – k S^2 & -k S C + m C^2 -m S^2 – k C S \\
-k C S – m S^2 + m C^2 – k S C & k S^2 – m C S -m S C – k C^2
\end{bmatrix} \\
&=
\begin{bmatrix}
k \cos(2 \theta_k) + m \sin(2 \theta_k) & m \cos(2 \theta_k) – k \sin(2 \theta_k) \\
m \cos(2 \theta_k) – k \sin(2 \theta_k) & -k \cos(2 \theta_k) – m \sin(2 \theta_k) \\
\end{bmatrix}.
\end{aligned}
\end{equation}

This gives

\begin{equation}\label{eqn:qmLecture9:320}
\begin{aligned}
E_k
\begin{bmatrix}
1 \\
0
\end{bmatrix}
&=
\begin{bmatrix}
k \cos(2 \theta_k) + m \sin(2 \theta_k) \\
m \cos(2 \theta_k) – k \sin(2 \theta_k) \\
\end{bmatrix} \\
&=
\begin{bmatrix}
k & m \\
m & -k
\end{bmatrix}
\begin{bmatrix}
\cos(2 \theta_k) \\
\sin(2 \theta_k) \\
\end{bmatrix}.
\end{aligned}
\end{equation}

Adding back in the \(\Hbar\)’s and \(c\)’s this is

\begin{equation}\label{eqn:qmLecture9:340}
\begin{aligned}
\begin{bmatrix}
\cos(2 \theta_k) \\
\sin(2 \theta_k) \\
\end{bmatrix}
&=
\frac{E_k}{-(\Hbar k c)^2 -(m c^2)^2}
\begin{bmatrix}
– \Hbar k c & – m c^2 \\
– m c^2 & \Hbar k c
\end{bmatrix}
\begin{bmatrix}
1 \\
0
\end{bmatrix} \\
&=
\inv{E_k}
\begin{bmatrix}
\Hbar k c \\
m c^2
\end{bmatrix}.
\end{aligned}
\end{equation}

Question: Verify the Dirac current relationship.

Prove \ref{eqn:qmLecture9:240}.

Answer

The components of the Schrodinger equation are

\begin{equation}\label{eqn:qmLecture9:360}
\begin{aligned}
-i \Hbar \PD{t}{\psi_1} &= -i \Hbar c \PD{x}{\psi_1} + m c^2 \psi_2 \\
-i \Hbar \PD{t}{\psi_2} &= m c^2 \psi_1 + i \Hbar c \PD{x}{\psi_2},
\end{aligned}
\end{equation}

The conjugates of these are
\begin{equation}\label{eqn:qmLecture9:380}
\begin{aligned}
i \Hbar \PD{t}{\psi_1^\conj} &= i \Hbar c \PD{x}{\psi_1^\conj} + m c^2 \psi_2^\conj \\
i \Hbar \PD{t}{\psi_2^\conj} &= m c^2 \psi_1^\conj – i \Hbar c \PD{x}{\psi_2^\conj}.
\end{aligned}
\end{equation}

This gives
\begin{equation}\label{eqn:qmLecture9:400}
\begin{aligned}
i \Hbar \PD{t}{\rho}
&=
\lr{ i \Hbar c \PD{x}{\psi_1^\conj} + m c^2 \psi_2^\conj } \psi_1 \\
&+ \psi_1^\conj \lr{ i \Hbar c \PD{x}{\psi_1} – m c^2 \psi_2 } \\
&+ \lr{ m c^2 \psi_1^\conj – i \Hbar c \PD{x}{\psi_2^\conj} } \psi_2 \\
&+ \psi_2^\conj \lr{ -m c^2 \psi_1 – i \Hbar c \PD{x}{\psi_2} }.
\end{aligned}
\end{equation}

All the non-derivative terms cancel leaving

\begin{equation}\label{eqn:qmLecture9:420}
\inv{c} \PD{t}{\rho}
=
\PD{x}{\psi_1^\conj} \psi_1
+ \psi_1^\conj \PD{x}{\psi_1}
– \PD{x}{\psi_2^\conj} \psi_2
– \psi_2^\conj \PD{x}{\psi_2}
=
\PD{x}{}
\lr{
\psi_1^\conj \psi_1 –
\psi_2^\conj \psi_2
}.
\end{equation}

PHY1520H Graduate Quantum Mechanics. Lecture 8: Dirac equation in 1D. Taught by Prof. Arun Paramekanti

October 13, 2015 phy1520 No comments , , , ,

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

Schrodinger Derivation

Recall that a “derivation” of the Schrodinger equation can be associated with the following equivalences

\begin{equation}\label{eqn:qmLecture8:300}
E \leftrightarrow \Hbar \omega \leftrightarrow i \Hbar \PD{t}{}
\end{equation}
\begin{equation}\label{eqn:qmLecture8:320}
p \leftrightarrow \Hbar k \leftrightarrow -i \Hbar \PD{t}{}
\end{equation}

so that the classical energy relationship

\begin{equation}\label{eqn:qmLecture8:20}
E = \frac{p^2}{2m}
\end{equation}

takes the form

\begin{equation}\label{eqn:qmLecture8:40}
i \Hbar \PD{t}{} = -\frac{\Hbar^2}{2m}.
\end{equation}

How do we do this in a relativistic context where the energy momentum relationship is

\begin{equation}\label{eqn:qmLecture8:60}
E = \sqrt{ p^2 c^2 + m^2 c^4 } \approx m c^2 + \frac{p^2}{2m} + \cdots
\end{equation}

where \( m \) is the rest mass and \( c \) is the speed of light.

Attempt I

\begin{equation}\label{eqn:qmLecture8:80}
E = m c^2 + \frac{p^2}{2m} + (…) p^4 + (…) p^6 + \cdots
\end{equation}

First order in time, but infinite order in space \( \partial/\partial x \). Useless.

Attempt II

\begin{equation}\label{eqn:qmLecture8:100}
E^2 = p^2 c^2 + m^2 c^4.
\end{equation}

This gives

\begin{equation}\label{eqn:qmLecture8:120}
-\Hbar^2 \PDSq{t}{\psi} = – \Hbar^2 c^2 \PDSq{x}{\psi} + m^2 c^4 \psi.
\end{equation}

This is the Klein-Gordon equation, which is second order in time.

Attempt III

Suppose that we have the matrix

\begin{equation}\label{eqn:qmLecture8:140}
\begin{bmatrix}
p c & m c^2 \\
m c^2 & – p c
\end{bmatrix},
\end{equation}

or

\begin{equation}\label{eqn:qmLecture8:160}
\begin{bmatrix}
m c^2 & i p c \\
-i p c & – m c^2
\end{bmatrix},
\end{equation}

These both happen to have eigenvalues \( \lambda_{\pm} = \pm \sqrt{p^2 c^2} \). For those familiar with the Dirac matrices, this amounts to a choice for different representations of the gamma matrices.

Working with \ref{eqn:qmLecture8:140}, which has some nicer features than other possible representations, we seek a state

\begin{equation}\label{eqn:qmLecture8:180}
\Bpsi =
\begin{bmatrix}
\psi_1(x, t) \\
\psi_2(x, t) \\
\end{bmatrix},
\end{equation}

where we aim to write down an equation for this composite state.

\begin{equation}\label{eqn:qmLecture8:200}
i \Hbar \PD{t}{\Bpsi} = \BH \Bpsi
\end{equation}

Assuming the matrix is the Hamiltonian, multiplying that with the composite state gives

\begin{equation}\label{eqn:qmLecture8:220}
\begin{aligned}
\begin{bmatrix}
i \Hbar \PD{t}{\psi_1} \\
i \Hbar \PD{t}{\psi_1}
\end{bmatrix}
&=
\begin{bmatrix}
\hat{p} c & m c^2 \\
m c^2 & – \hat{p} c
\end{bmatrix}
\begin{bmatrix}
\psi_1(x, t) \\
\psi_2(x, t) \\
\end{bmatrix} \\
&=
\begin{bmatrix}
\hat{p} c \psi_1 + m c^2 \psi_2 \\
m c^2 \psi_1 – \hat{p} c \psi_2
\end{bmatrix}.
\end{aligned}
\end{equation}

What happens when we square this?

\begin{equation}\label{eqn:qmLecture8:240}
\begin{aligned}
\lr{ i \Hbar \PD{t}{} }^2 \Bpsi
&= \BH \BH \Bpsi \\
&=
\begin{bmatrix}
\hat{p} c & m c^2 \\
m c^2 & – \hat{p} c
\end{bmatrix}
\begin{bmatrix}
\hat{p} c & m c^2 \\
m c^2 & – \hat{p} c
\end{bmatrix}
\Bpsi \\
&=
\begin{bmatrix}
\hat{p}^2 c^2 + m^2 c^4 & 0 \\
0 & \hat{p}^2 c^2 + m^2 c^4 \\
\end{bmatrix}.
\end{aligned}
\end{equation}

That is
\begin{equation}\label{eqn:qmLecture8:260}
– \Hbar^2 \PDSq{t}{} \Bpsi
=
\lr{ \hat{p}^2 c^2 + m^2 c^4 } \mathbf{1} \Bpsi,
\end{equation}

or more exactly

\begin{equation}\label{eqn:qmLecture8:280}
– \Hbar^2 \PDSq{t}{} \psi_{1,2}
=
\lr{ \hat{p}^2 c^2 + m^2 c^4 } \psi_{1,2}.
\end{equation}

This recovers the Klein Gordon equation for each of the wave functions \( \psi_1, \psi_2 \).

Instead of squaring the operators, lets try to solve the first order equation. To do so we’ll want to diagonalize \( \BH \).

Before doing that, let’s write out the Hamiltonian in an alternate but useful form

\begin{equation}\label{eqn:qmLecture8:340}
\BH
=
\hat{p} c
\begin{bmatrix}
1 & 0 \\
0 & -1
\end{bmatrix}
+
m c^2
\begin{bmatrix}
0 & 1 \\
1 & 0
\end{bmatrix}
= \hat{p} c \hat{\sigma}_z + m c^2 \hat{\sigma}_x.
\end{equation}

We have two types of operators in the mix here. We have matrix operators that act on the wave function matrices, as well as derivative operators that act on the components of those matrices.

We have

\begin{equation}\label{eqn:qmLecture8:360}
\hat{\sigma}_z
\begin{bmatrix}
\psi_1 \\
\psi_2 \\
\end{bmatrix}
=
\begin{bmatrix}
\psi_1 \\
-\psi_2 \\
\end{bmatrix},
\end{equation}

and

\begin{equation}\label{eqn:qmLecture8:380}
\hat{\sigma}_x
\begin{bmatrix}
\psi_1 \\
\psi_2 \\
\end{bmatrix}
=
\begin{bmatrix}
\psi_2 \\
\psi_1 \\
\end{bmatrix}.
\end{equation}

Because the derivative actions of \( \hat{p} \) and the matrix operators are independent, we see that these operators commute. For example

\begin{equation}\label{eqn:qmLecture8:400}
\hat{\sigma}_z \hat{p}
\begin{bmatrix}
\psi_1 \\
\psi_2 \\
\end{bmatrix}
=
\hat{\sigma}_z
\begin{bmatrix}
-i \Hbar \PD{x}{\psi_1} \\
-i \Hbar \PD{x}{\psi_2} \\
\end{bmatrix}
=
\begin{bmatrix}
-i \Hbar \PD{x}{\psi_1} \\
i \Hbar \PD{x}{\psi_2} \\
\end{bmatrix}
=
\hat{p}
\hat{\sigma}_z
\begin{bmatrix}
\psi_1 \\
\psi_2 \\
\end{bmatrix}.
\end{equation}

Diagonalizing it

Suppose the wave function matrix has the structure

\begin{equation}\label{eqn:qmLecture8:420}
\Bpsi =
\begin{bmatrix}
f_{+} \\
f_{-} \\
\end{bmatrix}
e^{i k x}.
\end{equation}

We’ll plug this into the Schrodinger equation and see what we get.