Hamiltonian

Quantum SHO ladder operators as a diagonal change of basis for the Heisenberg EOMs

August 19, 2015 phy1520 , , , , , , , , , , , , , , , , , , ,

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

Many authors pull the definitions of the raising and lowering (or ladder) operators out of their butt with no attempt at motivation. This is pointed out nicely in [1] by Eli along with one justification based on factoring the Hamiltonian.

In [2] is a small exception to the usual presentation. In that text, these operators are defined as usual with no motivation. However, after the utility of these operators has been shown, the raising and lowering operators show up in a context that does provide that missing motivation as a side effect.
It doesn’t look like the author was trying to provide a motivation, but it can be interpreted that way.

When seeking the time evolution of Heisenberg-picture position and momentum operators, we will see that those solutions can be trivially expressed using the raising and lowering operators. No special tools nor black magic is required to find the structure of these operators. Unfortunately, we must first switch to both the Heisenberg picture representation of the position and momentum operators, and also employ the Heisenberg equations of motion. Neither of these last two fit into standard narrative of most introductory quantum mechanics treatments. We will also see that these raising and lowering “operators” could also be introduced in classical mechanics, provided we were attempting to solve the SHO system using the Hamiltonian equations of motion.

I’ll outline this route to finding the structure of the ladder operators below. Because these are encountered trying to solve the time evolution problem, I’ll first show a simpler way to solve that problem. Because that simpler method depends a bit on lucky observation and is somewhat unstructured, I’ll then outline a more structured procedure that leads to the ladder operators directly, also providing the solution to the time evolution problem as a side effect.

The starting point is the Heisenberg equations of motion. For a time independent Hamiltonian \( H \), and a Heisenberg operator \( A^{(H)} \), those equations are

\begin{equation}\label{eqn:harmonicOscDiagonalize:20}
\ddt{A^{(H)}} = \inv{i \Hbar} \antisymmetric{A^{(H)}}{H}.
\end{equation}

Here the Heisenberg operator \( A^{(H)} \) is related to the Schrodinger operator \( A^{(S)} \) by

\begin{equation}\label{eqn:harmonicOscDiagonalize:60}
A^{(H)} = U^\dagger A^{(S)} U,
\end{equation}

where \( U \) is the time evolution operator. For this discussion, we need only know that \( U \) commutes with \( H \), and do not need to know the specific structure of that operator. In particular, the Heisenberg equations of motion take the form

\begin{equation}\label{eqn:harmonicOscDiagonalize:80}
\begin{aligned}
\ddt{A^{(H)}}
&= \inv{i \Hbar}
\antisymmetric{A^{(H)}}{H} \\
&= \inv{i \Hbar}
\antisymmetric{U^\dagger A^{(S)} U}{H} \\
&= \inv{i \Hbar}
\lr{
U^\dagger A^{(S)} U H
– H U^\dagger A^{(S)} U
} \\
&= \inv{i \Hbar}
\lr{
U^\dagger A^{(S)} H U
– U^\dagger H A^{(S)} U
} \\
&= \inv{i \Hbar} U^\dagger \antisymmetric{A^{(S)}}{H} U.
\end{aligned}
\end{equation}

The Hamiltonian for the harmonic oscillator, with Schrodinger-picture position and momentum operators \( x, p \) is

\begin{equation}\label{eqn:harmonicOscDiagonalize:40}
H = \frac{p^2}{2m} + \inv{2} m \omega^2 x^2,
\end{equation}

so the equations of motions are

\begin{equation}\label{eqn:harmonicOscDiagonalize:100}
\begin{aligned}
\ddt{x^{(H)}}
&= \inv{i \Hbar} U^\dagger \antisymmetric{x}{H} U \\
&= \inv{i \Hbar} U^\dagger \antisymmetric{x}{\frac{p^2}{2m}} U \\
&= \inv{2 m i \Hbar} U^\dagger \lr{ i \Hbar \PD{p}{p^2} } U \\
&= \inv{m } U^\dagger p U \\
&= \inv{m } p^{(H)},
\end{aligned}
\end{equation}

and
\begin{equation}\label{eqn:harmonicOscDiagonalize:120}
\begin{aligned}
\ddt{p^{(H)}}
&= \inv{i \Hbar} U^\dagger \antisymmetric{p}{H} U \\
&= \inv{i \Hbar} U^\dagger \antisymmetric{p}{\inv{2} m \omega^2 x^2 } U \\
&= \frac{m \omega^2}{2 i \Hbar} U^\dagger \lr{ -i \Hbar \PD{x}{x^2} } U \\
&= -m \omega^2 U^\dagger x U \\
&= -m \omega^2 x^{(H)}.
\end{aligned}
\end{equation}

In the Heisenberg picture the equations of motion are precisely those of classical Hamiltonian mechanics, except that we are dealing with operators instead of scalars

\begin{equation}\label{eqn:harmonicOscDiagonalize:140}
\begin{aligned}
\ddt{p^{(H)}} &= -m \omega^2 x^{(H)} \\
\ddt{x^{(H)}} &= \inv{m } p^{(H)}.
\end{aligned}
\end{equation}

In the text the ladder operators are used to simplify the solution of these coupled equations, since they can decouple them. That’s not really required since we can solve them directly in matrix form with little work

\begin{equation}\label{eqn:harmonicOscDiagonalize:160}
\ddt{}
\begin{bmatrix}
p^{(H)} \\
x^{(H)}
\end{bmatrix}
=
\begin{bmatrix}
0 & -m \omega^2 \\
\inv{m} & 0
\end{bmatrix}
\begin{bmatrix}
p^{(H)} \\
x^{(H)}
\end{bmatrix},
\end{equation}

or, with length scaled variables

\begin{equation}\label{eqn:harmonicOscDiagonalize:180}
\begin{aligned}
\ddt{}
\begin{bmatrix}
\frac{p^{(H)}}{m \omega} \\
x^{(H)}
\end{bmatrix}
&=
\begin{bmatrix}
0 & -\omega \\
\omega & 0
\end{bmatrix}
\begin{bmatrix}
\frac{p^{(H)}}{m \omega} \\
x^{(H)}
\end{bmatrix} \\
&=
-i \omega
\begin{bmatrix} 0 & -i \\ i & 0 \\ \end{bmatrix}
\begin{bmatrix}
\frac{p^{(H)}}{m \omega} \\
x^{(H)}
\end{bmatrix} \\
&=
-i \omega
\sigma_y
\begin{bmatrix}
\frac{p^{(H)}}{m \omega} \\
x^{(H)}
\end{bmatrix}.
\end{aligned}
\end{equation}

Writing \( y = \begin{bmatrix} \frac{p^{(H)}}{m \omega} \\ x^{(H)} \end{bmatrix} \), the solution can then be written immediately as

\begin{equation}\label{eqn:harmonicOscDiagonalize:200}
\begin{aligned}
y(t)
&=
\exp\lr{ -i \omega \sigma_y t } y(0) \\
&=
\lr{ \cos \lr{ \omega t } I – i \sigma_y \sin\lr{ \omega t } } y(0) \\
&=
\begin{bmatrix}
\cos\lr{ \omega t } & \sin\lr{ \omega t } \\
-\sin\lr{ \omega t } & \cos\lr{ \omega t }
\end{bmatrix}
y(0),
\end{aligned}
\end{equation}

or

\begin{equation}\label{eqn:harmonicOscDiagonalize:220}
\begin{aligned}
\frac{p^{(H)}(t)}{m \omega} &= \cos\lr{ \omega t } \frac{p^{(H)}(0)}{m \omega} + \sin\lr{ \omega t } x^{(H)}(0) \\
x^{(H)}(t) &= -\sin\lr{ \omega t } \frac{p^{(H)}(0)}{m \omega} + \cos\lr{ \omega t } x^{(H)}(0).
\end{aligned}
\end{equation}

This solution depends on being lucky enough to recognize that the matrix has a Pauli matrix as a factor (which squares to unity, and allows the exponential to be evaluated easily.)

If we hadn’t been that observant, then the first tool we’d have used instead would have been to diagonalize the matrix. For such diagonalization, it’s natural to work in completely dimensionless variables. Such a non-dimensionalisation can be had by defining

\begin{equation}\label{eqn:harmonicOscDiagonalize:240}
x_0 = \sqrt{\frac{\Hbar}{m \omega}},
\end{equation}

and dividing the working (operator) variables through by those values. Let \( z = \inv{x_0} y \), and \( \tau = \omega t \) so that the equations of motion are

\begin{equation}\label{eqn:harmonicOscDiagonalize:260}
\frac{dz}{d\tau}
=
\begin{bmatrix}
0 & -1 \\
1 & 0
\end{bmatrix}
z.
\end{equation}

This matrix can be diagonalized as

\begin{equation}\label{eqn:harmonicOscDiagonalize:280}
A
=
\begin{bmatrix}
0 & -1 \\
1 & 0
\end{bmatrix}
=
V
\begin{bmatrix}
i & 0 \\
0 & -i
\end{bmatrix}
V^{-1},
\end{equation}

where

\begin{equation}\label{eqn:harmonicOscDiagonalize:300}
V =
\inv{\sqrt{2}}
\begin{bmatrix}
i & -i \\
1 & 1
\end{bmatrix}.
\end{equation}

The equations of motion can now be written

\begin{equation}\label{eqn:harmonicOscDiagonalize:320}
\frac{d}{d\tau} \lr{ V^{-1} z } =
\begin{bmatrix}
i & 0 \\
0 & -i
\end{bmatrix}
\lr{ V^{-1} z }.
\end{equation}

This final change of variables \( V^{-1} z \) decouples the system as desired. Expanding that gives

\begin{equation}\label{eqn:harmonicOscDiagonalize:340}
\begin{aligned}
V^{-1} z
&=
\inv{\sqrt{2}}
\begin{bmatrix}
-i & 1 \\
i & 1
\end{bmatrix}
\begin{bmatrix}
\frac{p^{(H)}}{x_0 m \omega} \\
\frac{x^{(H)}}{x_0}
\end{bmatrix} \\
&=
\inv{\sqrt{2} x_0}
\begin{bmatrix}
-i \frac{p^{(H)}}{m \omega} + x^{(H)} \\
i \frac{p^{(H)}}{m \omega} + x^{(H)}
\end{bmatrix} \\
&=
\begin{bmatrix}
a^\dagger \\
a
\end{bmatrix},
\end{aligned}
\end{equation}

where
\begin{equation}\label{eqn:harmonicOscDiagonalize:n}
\begin{aligned}
a^\dagger &= \sqrt{\frac{m \omega}{2 \Hbar}} \lr{ -i \frac{p^{(H)}}{m \omega} + x^{(H)} } \\
a &= \sqrt{\frac{m \omega}{2 \Hbar}} \lr{ i \frac{p^{(H)}}{m \omega} + x^{(H)} }.
\end{aligned}
\end{equation}

Lo and behold, we have the standard form of the raising and lowering operators, and can write the system equations as

\begin{equation}\label{eqn:harmonicOscDiagonalize:360}
\begin{aligned}
\ddt{a^\dagger} &= i \omega a^\dagger \\
\ddt{a} &= -i \omega a.
\end{aligned}
\end{equation}

It is actually a bit fluky that this matched exactly, since we could have chosen eigenvectors that differ by constant phase factors, like

\begin{equation}\label{eqn:harmonicOscDiagonalize:380}
V = \inv{\sqrt{2}}
\begin{bmatrix}
i e^{i\phi} & -i e^{i \psi} \\
1 e^{i\phi} & e^{i \psi}
\end{bmatrix},
\end{equation}

so

\begin{equation}\label{eqn:harmonicOscDiagonalize:341}
\begin{aligned}
V^{-1} z
&=
\frac{e^{-i(\phi + \psi)}}{\sqrt{2}}
\begin{bmatrix}
-i e^{i\psi} & e^{i \psi} \\
i e^{i\phi} & e^{i \phi}
\end{bmatrix}
\begin{bmatrix}
\frac{p^{(H)}}{x_0 m \omega} \\
\frac{x^{(H)}}{x_0}
\end{bmatrix} \\
&=
\inv{\sqrt{2} x_0}
\begin{bmatrix}
-i e^{i\phi} \frac{p^{(H)}}{m \omega} + e^{i\phi} x^{(H)} \\
i e^{i\psi} \frac{p^{(H)}}{m \omega} + e^{i\psi} x^{(H)}
\end{bmatrix} \\
&=
\begin{bmatrix}
e^{i\phi} a^\dagger \\
e^{i\psi} a
\end{bmatrix}.
\end{aligned}
\end{equation}

To make the resulting pairs of operators Hermitian conjugates, we’d want to constrain those constant phase factors by setting \( \phi = -\psi \). If we were only interested in solving the time evolution problem no such additional constraints are required.

The raising and lowering operators are seen to naturally occur when seeking the solution of the Heisenberg equations of motion. This is found using the standard technique of non-dimensionalisation and then seeking a change of basis that diagonalizes the system matrix. Because the Heisenberg equations of motion are identical to the classical Hamiltonian equations of motion in this case, what we call the raising and lowering operators in quantum mechanics could also be utilized in the classical simple harmonic oscillator problem. However, in a classical context we wouldn’t have a justification to call this more than a change of basis.

References

[1] Eli Lansey. The Quantum Harmonic Oscillator Ladder Operators, 2009. URL http://behindtheguesses.blogspot.ca/2009/03/quantum-harmonic-oscillator-ladder.html. [Online; accessed 18-August-2015].

[2] Jun John Sakurai and Jim J Napolitano. Modern quantum mechanics, chapter {Time Development of the Oscillator}. Pearson Higher Ed, 2014.

Heisenberg picture position commutator

August 14, 2015 phy1520 , , , ,

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

Question: Heisenberg picture position commutator ([1] pr. 2.5)

Evaluate

\begin{equation}\label{eqn:positionCommutator:20}
\antisymmetric{x(t)}{x(0)},
\end{equation}

for a Heisenberg picture operator \( x(t) \) for a free particle.

Answer

The free particle Hamiltonian is

\begin{equation}\label{eqn:positionCommutator:40}
H = \frac{p^2}{2m},
\end{equation}

so the time evolution operator is

\begin{equation}\label{eqn:positionCommutator:60}
U(t) = e^{-i p^2 t/(2 m \Hbar)}.
\end{equation}

The Heisenberg picture position operator is

\begin{equation}\label{eqn:positionCommutator:80}
\begin{aligned}
x^\textrm{H}
&= U^\dagger x U \\
&= e^{i p^2 t/(2 m \Hbar)} x e^{-i p^2 t/(2 m \Hbar)} \\
&= \sum_{k = 0}^\infty \inv{k!} \lr{ \frac{i p^2 t}{2 m \Hbar} }^k
x
e^{-i p^2 t/(2 m \Hbar)} \\
&= \sum_{k = 0}^\infty \inv{k!} \lr{ \frac{i t}{2 m \Hbar} }^k p^{2k} x
e^{-i p^2 t/(2 m \Hbar)} \\
&=
\sum_{k = 0}^\infty \inv{k!} \lr{ \frac{i t}{2 m \Hbar} }^k \lr{ \antisymmetric{p^{2k}}{x} + x p^{2k} }
e^{-i p^2 t/(2 m \Hbar)} \\
&= x +
\sum_{k = 0}^\infty \inv{k!} \lr{ \frac{i t}{2 m \Hbar} }^k \antisymmetric{p^{2k}}{x}
e^{-i p^2 t/(2 m \Hbar)} \\
&= x +
\sum_{k = 0}^\infty \inv{k!} \lr{ \frac{i t}{2 m \Hbar} }^k \lr{ -i \Hbar \PD{p}{p^{2k}} }
e^{-i p^2 t/(2 m \Hbar)} \\
&= x +
\sum_{k = 0}^\infty \inv{k!} \lr{ \frac{i t}{2 m \Hbar} }^k \lr{ -i \Hbar 2 k p^{2 k -1} }
e^{-i p^2 t/(2 m \Hbar)} \\
&= x +
-2 i \Hbar p \frac{i t}{2 m \Hbar} \sum_{k = 1}^\infty \inv{(k-1)!} \lr{ \frac{i t}{2 m \Hbar} }^{k-1} p^{2(k – 1)}
e^{-i p^2 t/(2 m \Hbar)} \\
&= x + t \frac{p}{m}.
\end{aligned}
\end{equation}

This has the structure of a classical free particle \( x(t) = x + v t \), but in this case \( x,p \) are operators.

The evolved position commutator is
\begin{equation}\label{eqn:positionCommutator:100}
\begin{aligned}
\antisymmetric{x(t)}{x(0)}
&=
\antisymmetric{x + t p/m}{x} \\
&=
\frac{t}{m} \antisymmetric{p}{x} \\
&=
-i \Hbar \frac{t}{m}.
\end{aligned}
\end{equation}

Compare this to the classical Poisson bracket
\begin{equation}\label{eqn:positionCommutator:120}
\antisymmetric{x(t)}{x(0)}_{\textrm{classical}}
=
\PD{x}{}\lr{x + p t/m} \PD{p}{x} – \PD{p}{}\lr{x + p t/m} \PD{x}{x}
=
– \frac{t}{m}.
\end{equation}

This has the expected relation \( \antisymmetric{x(t)}{x(0)} = i \Hbar \antisymmetric{x(t)}{x(0)}_{\textrm{classical}} \).

References

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

Heisenberg picture spin precession

August 13, 2015 phy1520 , , , , ,

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

Question: Heisenberg picture spin precession ([1] pr. 2.1)

For the spin Hamiltonian

\begin{equation}\label{eqn:heisenbergSpinPrecession:20}
H = -\frac{e B}{m c} S_z = \omega S_z,
\end{equation}

express and solve the Heisenberg equations of motion for \( S_x(t), S_y(t) \), and \( S_z(t) \).

Answer

The equations of motion are of the form

\begin{equation}\label{eqn:heisenbergSpinPrecession:40}
\begin{aligned}
\frac{dS_i^\textrm{H}}{dt}
&= \inv{i \Hbar} \antisymmetric{S_i^\textrm{H}}{H} \\
&= \inv{i \Hbar} \antisymmetric{U^\dagger S_i U}{H} \\
&= \inv{i \Hbar} \lr{U^\dagger S_i U H – H U^\dagger S_i U } \\
&= \inv{i \Hbar} U^\dagger \lr{ S_i H – H S_i } U \\
&= \frac{\omega}{i \Hbar} U^\dagger \antisymmetric{ S_i}{S_z } U.
\end{aligned}
\end{equation}

These are

\begin{equation}\label{eqn:heisenbergSpinPrecession:60}
\begin{aligned}
\frac{dS_x^\textrm{H}}{dt} &= -\omega U^\dagger S_y U \\
\frac{dS_y^\textrm{H}}{dt} &= \omega U^\dagger S_x U \\
\frac{dS_z^\textrm{H}}{dt} &= 0.
\end{aligned}
\end{equation}

To completely specify these equations, we need to expand \( U(t) \), which is

\begin{equation}\label{eqn:heisenbergSpinPrecession:80}
\begin{aligned}
U(t)
&= e^{-i H t /\Hbar} \\
&= e^{-i \omega S_z t /\Hbar} \\
&= e^{-i \omega \sigma_z t /2} \\
&= \cos\lr{ \omega t/2 } -i \sigma_z \sin\lr{ \omega t/2 } \\
&=
\begin{bmatrix}
\cos\lr{ \omega t/2 } -i \sin\lr{ \omega t/2 } & 0 \\
0 & \cos\lr{ \omega t/2 } + i \sin\lr{ \omega t/2 }
\end{bmatrix} \\
&=
\begin{bmatrix}
e^{-i\omega t/2} & 0 \\
0 & e^{i\omega t/2}
\end{bmatrix}.
\end{aligned}
\end{equation}

The equations of motion can now be written out in full. To do so seems a bit silly since we also know that \( S_x^\textrm{H} = U^\dagger S_x U, S_y^\textrm{H} U^\dagger S_x U \). However, if that is temporarily forgotten, we can show that the Heisenberg equations of motion can be solved for these too.

\begin{equation}\label{eqn:heisenbergSpinPrecession:100}
\begin{aligned}
U^\dagger S_x U
&=
\frac{\Hbar}{2}
\begin{bmatrix}
e^{i\omega t/2} & 0 \\
0 & e^{-i\omega t/2}
\end{bmatrix}
\begin{bmatrix} 0 & 1 \\ 1 & 0 \\ \end{bmatrix}
\begin{bmatrix}
e^{-i\omega t/2} & 0 \\
0 & e^{i\omega t/2}
\end{bmatrix} \\
&=
\frac{\Hbar}{2}
\begin{bmatrix}
0 & e^{i\omega t/2} \\
e^{-i\omega t/2} & 0
\end{bmatrix}
\begin{bmatrix}
e^{-i\omega t/2} & 0 \\
0 & e^{i\omega t/2}
\end{bmatrix} \\
&=
\frac{\Hbar}{2}
\begin{bmatrix}
0 & e^{i\omega t} \\
e^{-i\omega t} & 0
\end{bmatrix},
\end{aligned}
\end{equation}

and
\begin{equation}\label{eqn:heisenbergSpinPrecession:120}
\begin{aligned}
U^\dagger S_y U
&=
\frac{\Hbar}{2}
\begin{bmatrix}
e^{i\omega t/2} & 0 \\
0 & e^{-i\omega t/2}
\end{bmatrix}
\begin{bmatrix} 0 & -i \\ i & 0 \\ \end{bmatrix}
\begin{bmatrix}
e^{-i\omega t/2} & 0 \\
0 & e^{i\omega t/2}
\end{bmatrix} \\
&=
\frac{i\Hbar}{2}
\begin{bmatrix}
0 & -e^{i\omega t/2} \\
e^{-i\omega t/2} & 0
\end{bmatrix}
\begin{bmatrix}
e^{-i\omega t/2} & 0 \\
0 & e^{i\omega t/2}
\end{bmatrix} \\
&=
\frac{i \Hbar}{2}
\begin{bmatrix}
0 & -e^{i\omega t} \\
e^{-i\omega t} & 0
\end{bmatrix}.
\end{aligned}
\end{equation}

The equations of motion are now fully specified

\begin{equation}\label{eqn:heisenbergSpinPrecession:140}
\begin{aligned}
\frac{dS_x^\textrm{H}}{dt} &=
-\frac{i \Hbar \omega}{2}
\begin{bmatrix}
0 & -e^{i\omega t} \\
e^{-i\omega t} & 0
\end{bmatrix} \\
\frac{dS_y^\textrm{H}}{dt} &=
\frac{\Hbar \omega}{2}
\begin{bmatrix}
0 & e^{i\omega t} \\
e^{-i\omega t} & 0
\end{bmatrix} \\
\frac{dS_z^\textrm{H}}{dt} &= 0.
\end{aligned}
\end{equation}

Integration gives

\begin{equation}\label{eqn:heisenbergSpinPrecession:160}
\begin{aligned}
S_x^\textrm{H} &=
\frac{\Hbar}{2}
\begin{bmatrix}
0 & e^{i\omega t} \\
e^{-i\omega t} & 0
\end{bmatrix} + C \\
S_y^\textrm{H} &=
\frac{\Hbar}{2}
\begin{bmatrix}
0 & -i e^{i\omega t} \\
i e^{-i\omega t} & 0
\end{bmatrix} + C \\
S_z^\textrm{H} &= C.
\end{aligned}
\end{equation}

The integration constants are fixed by the boundary condition \( S_i^\textrm{H}(0) = S_i \), so

\begin{equation}\label{eqn:heisenbergSpinPrecession:180}
\begin{aligned}
S_x^\textrm{H} &=
\frac{\Hbar}{2}
\begin{bmatrix}
0 & e^{i\omega t} \\
e^{-i\omega t} & 0
\end{bmatrix} \\
S_y^\textrm{H} &=
\frac{i \Hbar}{2}
\begin{bmatrix}
0 & – e^{i\omega t} \\
e^{-i\omega t} & 0
\end{bmatrix} \\
S_z^\textrm{H} &= S_z.
\end{aligned}
\end{equation}

Observe that these integrated values \( S_x^\textrm{H}, S_y^\textrm{H} \) match \ref{eqn:heisenbergSpinPrecession:100}, and \ref{eqn:heisenbergSpinPrecession:120} as expected.

References

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

Dynamics of non-Hermitian Hamiltonian

August 13, 2015 phy1520 , , , , ,

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

Question: Dynamics of non-Hermitian Hamiltonian ([1] pr. 2.2)

Revisiting an earlier Hamiltonian, but assuming it was entered incorrectly as

\begin{equation}\label{eqn:dynamicsNonHermitian:20}
H = H_{11} \ket{1}\bra{1}
+ H_{22} \ket{2}\bra{2}
+ H_{12} \ket{1}\bra{2}.
\end{equation}

What principle is now violated? Illustrate your point explicitly by attempting to solve the most generaqtl time-dependent problem using an illegal Hamiltonian of this kind. You may assume that \( H_{11} = H_{22} \) for simplicity.

Answer

In matrix form this Hamiltonian is

\begin{equation}\label{eqn:dynamicsNonHermitian:40}
\begin{aligned}
H
&=
\begin{bmatrix}
\bra{1} H \ket{1} & \bra{1} H \ket{2} \\
\bra{2} H \ket{1} & \bra{2} H \ket{2} \\
\end{bmatrix} \\
&=
\begin{bmatrix}
H_{11} & H_{12} \\
0 & H_{22} \\
\end{bmatrix}.
\end{aligned}
\end{equation}

This is not a Hermitian operator. What is the physical implication of this non-Hermicity? Consider the simpler case where \( H_{11} = H_{22} \). Such a Hamiltonian has the form

\begin{equation}\label{eqn:dynamicsNonHermitian:60}
H =
\begin{bmatrix}
a & b \\
0 & a
\end{bmatrix}.
\end{equation}

This has only one unique eigenvector ( \( (1,0) \), but we can still solve the time evolution equation

\begin{equation}\label{eqn:dynamicsNonHermitian:80}
i \Hbar \PD{t}{U} = H U,
\end{equation}

since for constant \( H \), we have

\begin{equation}\label{eqn:dynamicsNonHermitian:100}
U = e^{-i H t/\Hbar}.
\end{equation}

To exponentiate, note that we have

\begin{equation}\label{eqn:dynamicsNonHermitian:120}
{\begin{bmatrix}
a & b \\
0 & a
\end{bmatrix}}^n
=
\begin{bmatrix}
a^n & n a^{n-1} b \\
0 & a^n
\end{bmatrix}.
\end{equation}

To prove the induction, the \( n = 2 \) case follows easily

\begin{equation}\label{eqn:dynamicsNonHermitian:140}
\begin{bmatrix}
a & b \\
0 & a
\end{bmatrix}
\begin{bmatrix}
a & b \\
0 & a
\end{bmatrix}
=
\begin{bmatrix}
a^2 & 2 a b \\
0 & a^2
\end{bmatrix},
\end{equation}

as does the general case

\begin{equation}\label{eqn:dynamicsNonHermitian:160}
\begin{bmatrix}
a^n & n a^{n-1} b \\
0 & a^n
\end{bmatrix}
\begin{bmatrix}
a & b \\
0 & a
\end{bmatrix}
=
\begin{bmatrix}
a^{n+1} & (n +1 ) a^{n} b \\
0 & a^{n+1}
\end{bmatrix}.
\end{equation}

The exponential sum is thus
\begin{equation}\label{eqn:dynamicsNonHermitian:180}
e^{H \tau}
=
\begin{bmatrix}
e^{a \tau} & 0 + \frac{b \tau}{1!} + \frac{2 a b \tau^2}{2!} + \frac{3 a^2 b \tau^3}{3!} + \cdots \\
0 & e^{a \tau}
\end{bmatrix}.
\end{equation}

That sum simplifies to

\begin{equation}\label{eqn:dynamicsNonHermitian:200}
\frac{b \tau}{0!} + \frac{a b \tau^2}{1!} + \frac{a^2 b \tau^3}{2!} + \cdots \\
=
b \tau \lr{ 1 + \frac{a \tau}{1!} + \frac{(a \tau)^2}{2!} + \cdots }
=
b \tau e^{a \tau}.
\end{equation}

The exponential is thus
\begin{equation}\label{eqn:dynamicsNonHermitian:220}
e^{H \tau}
=
\begin{bmatrix}
e^{a\tau} & b \tau e^{a\tau} \\
0 & e^{a\tau}
\end{bmatrix}
=
\begin{bmatrix}
1 & b \tau \\
0 & 1
\end{bmatrix}
e^{a\tau}.
\end{equation}

In particular

\begin{equation}\label{eqn:dynamicsNonHermitian:240}
U = e^{-i H t/\Hbar} =
\begin{bmatrix}
1 & -i b t/\Hbar \\
0 & 1
\end{bmatrix}
e^{-i a t /\Hbar }.
\end{equation}

We can verify that this is a solution to \ref{eqn:dynamicsNonHermitian:80}. The left hand side is

\begin{equation}\label{eqn:dynamicsNonHermitian:260}
\begin{aligned}
i \Hbar \PD{t}{U}
&=
i \Hbar
\begin{bmatrix}
-i a/\Hbar & -i b /\Hbar + (-i b t/\Hbar)(-i a/\Hbar) \\
0 & -i a /\Hbar
\end{bmatrix}
e^{-i a t /\Hbar } \\
&=
\begin{bmatrix}
a & b – i a b t/\Hbar \\
0 & a
\end{bmatrix}
e^{-i a t /\Hbar },
\end{aligned}
\end{equation}

and for the right hand side
\begin{equation}\label{eqn:dynamicsNonHermitian:280}
\begin{aligned}
H U
&=
\begin{bmatrix}
a & b \\
0 & a
\end{bmatrix}
\begin{bmatrix}
1 & -i b t/\Hbar \\
0 & 1
\end{bmatrix}
e^{-i a t /\Hbar } \\
&=
\begin{bmatrix}
a & b – i a b t/\Hbar \\
0 & a
\end{bmatrix}
e^{-i a t /\Hbar } \\
&=
i \Hbar \PD{t}{U}.
\end{aligned}
\end{equation}

While the Schr\”{o}dinger is satisfied, we don’t have the unitary invertion physical property that is desired for the time evolution operator \( U \). Namely

\begin{equation}\label{eqn:dynamicsNonHermitian:300}
\begin{aligned}
U^\dagger U
&=
\begin{bmatrix}
1 & 0 \\
i b t/\Hbar & 1
\end{bmatrix}
e^{i a t /\Hbar }
\begin{bmatrix}
1 & -i b t/\Hbar \\
0 & 1
\end{bmatrix}
e^{-i a t /\Hbar } \\
&=
\begin{bmatrix}
1 & -i b t/\Hbar \\
i b t/\Hbar & (b t)^2/\Hbar^2
\end{bmatrix} \\
&\ne I.
\end{aligned}
\end{equation}

We required \( U^\dagger U = I \) for the time evolution operator, but don’t have that property for this non-Hermitian Hamiltonian.

References

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