expectation

Expectations for SHO Hamiltonian, and virial theorem.

October 15, 2015 phy1520 , , , , ,

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

Question: Expectations for SHO Hamiltonian, and virial theorem. ([1] pr. 2.3)

(a)

For a 1D SHO, compute \(
\bra{m} x \ket{n},
\bra{m} x^2 \ket{n},
\bra{m} p \ket{n},
\bra{m} p^2 \ket{n} \) and \( \bra{m} \symmetric{x}{p} \ket{n} \).

(b)

Verify the virial theorem is satisfied for energy eigenstates.

Answer

(a)

Using

\begin{equation}\label{eqn:shoExpectations:20}
\begin{aligned}
x &= \frac{x_0}{\sqrt{2}} \lr{ a + a^\dagger } \\
p &= \frac{i\Hbar}{x_0 \sqrt{2}} \lr{ a^\dagger – a} \\
a(t) &= a(0) e^{-i \omega t} \\
a(0) \ket{n} &= \sqrt{n} \ket{n-1} \\
a^\dagger(0) \ket{n} &= \sqrt{n+1} \ket{n+1} \\
x_0^2 &= \frac{\Hbar}{\omega m},
\end{aligned}
\end{equation}

we have

\begin{equation}\label{eqn:shoExpectations:40}
\begin{aligned}
\bra{m} x \ket{n}
&=
\frac{x_0}{\sqrt{2}} \bra{m} \lr{ a + a^\dagger } \ket{n} \\
&=
\frac{x_0}{\sqrt{2}} \bra{m}
\lr{
e^{-i \omega t} \sqrt{n} \ket{n-1}
+
e^{i \omega t} \sqrt{n+1} \ket{n+1}
} \\
&=
\frac{x_0}{\sqrt{2}} \lr{
\delta_{m, n-1} e^{-i \omega t} \sqrt{n}
+
\delta_{m, n+1} e^{i \omega t} \sqrt{n+1}
},
\end{aligned}
\end{equation}

\begin{equation}\label{eqn:shoExpectations:60}
\begin{aligned}
\bra{m} x^2 \ket{n}
&=
\frac{x_0^2}{2} \bra{m} \lr{ a + a^\dagger }^2 \ket{n} \\
&=
\frac{x_0^2}{2}
\lr{
e^{i \omega t} \sqrt{m} \bra{m-1}
+
e^{-i \omega t} \sqrt{m+1} \bra{m+1}
}
\lr{
e^{-i \omega t} \sqrt{n} \ket{n-1}
+
e^{i \omega t} \sqrt{n+1} \ket{n+1}
} \\
&=
\frac{x_0^2}{2}
\lr{
\delta_{m+1,n+1} \sqrt{(m+1)(n+1)}
+\delta_{m+1,n-1} \sqrt{(m+1)n} e^{-2 i \omega t}
+\delta_{m-1,n+1} \sqrt{m(n+1)} e^{2 i \omega t}
+\delta_{m-1,n-1} \sqrt{m n}
},
\end{aligned}
\end{equation}

\begin{equation}\label{eqn:shoExpectations:80}
\begin{aligned}
\bra{m} p \ket{n}
&=
\frac{i\Hbar}{\sqrt{2} x_0} \bra{m} \lr{ a^\dagger – a} \ket{n} \\
&=
\frac{i\Hbar}{\sqrt{2} x_0} \bra{m} \lr{
e^{i \omega t} \sqrt{n+1} \ket{n+1}

e^{-i \omega t} \sqrt{n} \ket{n-1}
} \\
&=
\frac{i\Hbar}{\sqrt{2} x_0} \lr{
\delta_{m,n+1} e^{i \omega t} \sqrt{n+1}

\delta_{m,n-1} e^{-i \omega t} \sqrt{n}
},
\end{aligned}
\end{equation}

\begin{equation}\label{eqn:shoExpectations:100}
\begin{aligned}
\bra{m} p^2 \ket{n}
&=
\frac{\Hbar^2}{2 x_0^2} \ket{m} \lr{ a – a^\dagger } \lr{ a^\dagger – a}
\ket{n} \\
&=
\frac{\Hbar^2}{2 x_0^2}
\lr{
-e^{-i \omega t} \sqrt{m+1} \bra{m+1}
+
e^{i \omega t} \sqrt{m} \bra{m-1}
}
\lr{
e^{i \omega t} \sqrt{n+1} \ket{n+1}

e^{-i \omega t} \sqrt{n} \ket{n-1}
} \\
&=
\frac{\Hbar^2}{2 x_0^2}
\lr{
\delta_{m+1,n+1} \sqrt{(m+1)(n+1)}
+\delta_{m+1,n-1} \sqrt{(m+1)n} e^{-2 i \omega t}
+\delta_{m-1,n+1} \sqrt{m(n+1)} e^{2 i \omega t}
+\delta_{m-1,n-1} \sqrt{m n}
}.
\end{aligned}
\end{equation}

For the anticommutator \( \symmetric{x}{p} \), we have

\begin{equation}\label{eqn:shoExpectations:120}
\begin{aligned}
\symmetric{x}{p}
&=
\frac{i\Hbar}{2}
\lr{
\lr{ a e^{-i \omega t} + a^\dagger e^{i \omega t} } \lr{ a^\dagger e^{i \omega t} – a e^{-i \omega t} }

\lr{ a^\dagger e^{i \omega t} – a e^{-i \omega t} }
\lr{ a e^{-i \omega t} + a^\dagger e^{i \omega t} }
} \\
&=
\frac{i\Hbar}{2}
\lr{
– a^2 e^{- 2 i \omega t}
+ (a^\dagger)^2 e^{ 2 i \omega t}
+ a a^\dagger
– a^\dagger a
+ a^2 e^{- 2 i \omega t}
– (a^\dagger)^2 e^{ 2 i \omega t}
– a^\dagger a
+ a a^\dagger
} \\
&=
i\Hbar
\lr{
a a^\dagger – a^\dagger a
},
\end{aligned}
\end{equation}

so

\begin{equation}\label{eqn:shoExpectations:140}
\begin{aligned}
\bra{m} \symmetric{x}{p} \ket{n}
&=
i\Hbar
\bra{m}
\lr{
a a^\dagger – a^\dagger a
}
\ket{n} \\
&=
i\Hbar
\bra{m}
\lr{
\sqrt{(n+1)^2}\ket{n}
-\sqrt{n^2}\ket{n}
} \\
&=
i\Hbar
\bra{m}
\lr{
2 n + 1
}
\ket{n}.
\end{aligned}
\end{equation}

(b)

For the SHO, the virial theorem requires \( \expectation{p^2/m} = \expectation{m \omega x^2} \). That momentum expectation with respect to the eigenstate \( \ket{n} \) is

\begin{equation}\label{eqn:shoExpectations:160}
\begin{aligned}
\expectation{p^2/m}
&=
\frac{\Hbar^2}{2 x_0^2 m}
\lr{
\sqrt{(n+1)(n+1)}
+
\sqrt{n n}
} \\
&=
\frac{\Hbar^2 m \omega}{2 \Hbar m} \lr{ 2 n + 1 } \\
&=
\Hbar \omega \lr{ n + \inv{2} }.
\end{aligned}
\end{equation}

For the position expectation we’ve got

\begin{equation}\label{eqn:shoExpectations:180}
\begin{aligned}
\expectation{m \omega x^2}
&=
\frac{m \omega^2 x_0^2}{2}
\lr{
\sqrt{(n+1)(n+1)}
+ \sqrt{n n}
} \\
&=
\frac{m \omega^2 \Hbar}{2 m \omega}
\lr{
\sqrt{(n+1)(n+1)}
+ \sqrt{n n}
} \\
&=
\frac{\omega \Hbar}{2 }
\lr{ 2 n + 1 } \\
&=
\omega \Hbar
\lr{ n + \inv{2} }.
\end{aligned}
\end{equation}

This shows that the virial theorem holds for the SHO Hamiltonian for eigenstates.

References

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

1D SHO linear superposition that maximizes expectation

October 7, 2015 phy1520 , , , , , ,

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

Question: 1D SHO linear superposition that maximizes expectation ([1] pr. 2.17)

For a 1D SHO

(a)

Construct a linear combination of \( \ket{0}, \ket{1} \) that maximizes \( \expectation{x} \) without using wave functions.

(b)

How does this state evolve with time?

(c)

Evaluate \( \expectation{x} \) using the Schrodinger picture.

(d)

Evaluate \( \expectation{x} \) using the Heisenberg picture.

(e)

Evaluate \( \expectation{(\Delta x)^2} \).

Answer

(a)

Forming

\begin{equation}\label{eqn:shoSuperposition:20}
\ket{\psi} = \frac{\ket{0} + \sigma \ket{1}}{\sqrt{1 + \Abs{\sigma}^2}}
\end{equation}

the position expectation is

\begin{equation}\label{eqn:shoSuperposition:40}
\bra{\psi} x \ket{\psi}
=
\inv{1 + \Abs{\sigma}^2} \lr{ \bra{0} + \sigma^\conj \bra{1} } \frac{x_0}{\sqrt{2}} \lr{ a^\dagger + a } \lr{ \ket{0} + \sigma \ket{1} }.
\end{equation}

Evaluating the action of the operators on the kets, we’ve got

\begin{equation}\label{eqn:shoSuperposition:60}
\lr{ a^\dagger + a } \lr{ \ket{0} + \sigma \ket{1} }
=
\ket{1} + \sqrt{2} \sigma \ket{2} + \sigma \ket{0}.
\end{equation}

The \( \ket{2} \) term is killed by the bras, leaving

\begin{equation}\label{eqn:shoSuperposition:80}
\begin{aligned}
\expectation{x}
&=
\inv{1 + \Abs{\sigma}^2} \frac{x_0}{\sqrt{2}} \lr{ \sigma + \sigma^\conj} \\
&=
\frac{\sqrt{2} x_0 \textrm{Re} \sigma}{1 + \Abs{\sigma}^2}.
\end{aligned}
\end{equation}

Any imaginary component in \( \sigma \) will reduce the expectation, so we are constrained to picking a real value.

The derivative of

\begin{equation}\label{eqn:shoSuperposition:100}
f(\sigma) = \frac{\sigma}{1 + \sigma^2},
\end{equation}

is

\begin{equation}\label{eqn:shoSuperposition:120}
f'(\sigma) = \frac{1 – \sigma^2}{(1 + \sigma^2)^2}.
\end{equation}

That has zeros at \( \sigma = \pm 1 \). The second derivative is

\begin{equation}\label{eqn:shoSuperposition:140}
f”(\sigma) = \frac{-2 \sigma (3 – \sigma^2)}{(1 + \sigma^2)^3}.
\end{equation}

That will be negative (maximum for the extreme value) at \( \sigma = 1 \), so the linear superposition of these first two energy eigenkets that maximizes the position expectation is

\begin{equation}\label{eqn:shoSuperposition:160}
\psi = \inv{\sqrt{2}}\lr{ \ket{0} + \ket{1} }.
\end{equation}

That maximized position expectation is

\begin{equation}\label{eqn:shoSuperposition:180}
\expectation{x}
=
\frac{x_0}{\sqrt{2}}.
\end{equation}

(b)

The time evolution is given by

\begin{equation}\label{eqn:shoSuperposition:200}
\begin{aligned}
\ket{\Psi(t)}
&= e^{-iH t/\Hbar} \inv{\sqrt{2}}\lr{ \ket{0} + \ket{1} } \\
&= \inv{\sqrt{2}}\lr{ e^{-i(0+ \ifrac{1}{2})\Hbar \omega t/\Hbar} \ket{0} +
e^{-i(1+ \ifrac{1}{2})\Hbar \omega t/\Hbar} \ket{1} } \\
&= \inv{\sqrt{2}}\lr{ e^{-i \omega t/2} \ket{0} + e^{-3 i \omega t/2} \ket{1} }.
\end{aligned}
\end{equation}

(c)

The position expectation in the Schrodinger representation is

\begin{equation}\label{eqn:shoSuperposition:220}
\begin{aligned}
\expectation{x(t)}
&=
\inv{2}
\lr{ e^{i \omega t/2} \bra{0} + e^{3 i \omega t/2} \bra{1} } \frac{x_0}{\sqrt{2}} \lr{ a^\dagger + a }
\lr{ e^{-i \omega t/2} \ket{0} + e^{-3 i \omega t/2} \ket{1} } \\
&=
\frac{x_0}{2\sqrt{2}}
\lr{ e^{i \omega t/2} \bra{0} + e^{3 i \omega t/2} \bra{1} }
\lr{ e^{-i \omega t/2} \ket{1} + e^{-3 i \omega t/2} \sqrt{2} \ket{2} + e^{-3 i \omega t/2} \ket{0} } \\
&=
\frac{x_0}{\sqrt{2}} \cos(\omega t).
\end{aligned}
\end{equation}

(d)

\begin{equation}\label{eqn:shoSuperposition:240}
\begin{aligned}
\expectation{x(t)}
&=
\inv{2}
\lr{ \bra{0} + \bra{1} } \frac{x_0}{\sqrt{2}}
\lr{ a^\dagger e^{i\omega t} + a e^{-i \omega t} }
\lr{ \ket{0} + \ket{1} } \\
&=
\frac{x_0}{2 \sqrt{2}}
\lr{ \bra{0} + \bra{1} }
\lr{ e^{i\omega t} \ket{1} + \sqrt{2} e^{i\omega t} \ket{2} + e^{-i \omega t} \ket{0} } \\
&=
\frac{x_0}{\sqrt{2}} \cos(\omega t),
\end{aligned}
\end{equation}

matching the calculation using the Schrodinger picture.

(e)

Let’s use the Heisenberg picture for the uncertainty calculation. Using the calculation above we have

\begin{equation}\label{eqn:shoSuperposition:260}
\begin{aligned}
\expectation{x^2}
&=
\inv{2} \frac{x_0^2}{2}
\lr{ e^{-i\omega t} \bra{1} + \sqrt{2} e^{-i\omega t} \bra{2} + e^{i \omega t} \bra{0} }
\lr{ e^{i\omega t} \ket{1} + \sqrt{2} e^{i\omega t} \ket{2} + e^{-i \omega t} \ket{0} } \\
&=
\frac{x_0^2}{4} \lr{ 1 + 2 + 1} \\
&=
x_0^2.
\end{aligned}
\end{equation}

The uncertainty is
\begin{equation}\label{eqn:shoSuperposition:280}
\begin{aligned}
\expectation{(\Delta x)^2}
&=
\expectation{x^2} – \expectation{x}^2 \\
&=
x_0^2 – \frac{x_0^2}{2} \cos^2(\omega t) \\
&=
\frac{x_0^2}{2} \lr{ 2 – \cos^2(\omega t) } \\
&=
\frac{x_0^2}{2} \lr{ 1 + \sin^2(\omega t) }
\end{aligned}
\end{equation}

References

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

SHO translation operator expectation

September 2, 2015 phy1520 , , , ,

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

Question: SHO translation operator expectation ([1] pr. 2.12)

Using the Heisenberg picture evaluate the expectation of the position operator \( \expectation{x} \) with respect to the initial time state

\begin{equation}\label{eqn:translationExpectation:20}
\ket{\alpha, 0} = e^{-i p_0 a/\Hbar} \ket{0},
\end{equation}

where \( p_0 \) is the initial time position operator, and \( a \) is a constant with dimensions of position.

Answer

Recall that the Heisenberg picture position operator expands to

\begin{equation}\label{eqn:translationExpectation:40}
x^{\textrm{H}}(t)
= U^\dagger x U
= x_0 \cos(\omega t) + \frac{p_0}{m \omega} \sin(\omega t),
\end{equation}

so the expectation of the position operator is
\begin{equation}\label{eqn:translationExpectation:60}
\begin{aligned}
\expectation{x}
&=
\bra{0} e^{i p_0 a/\Hbar} \lr{ x_0 \cos(\omega t) + \frac{p_0}{m \omega}
\sin(\omega t) } e^{-i p_0 a/\Hbar} \ket{0} \\
&=
\bra{0} \lr{ e^{i p_0 a/\Hbar} x_0 \cos(\omega t) e^{-i p_0 a/\Hbar} \cos(\omega t) + \frac{p_0}{m \omega} \sin(\omega t) } \ket{0}.
\end{aligned}
\end{equation}

The exponential sandwich above can be expanded using the Baker-Campbell-Hausdorff [2] formula

\begin{equation}\label{eqn:translationExpectation:80}
\begin{aligned}
e^{i p_0 a/\Hbar} x_0 e^{-i p_0 a/\Hbar}
&=
x_0
+ \frac{i a}{\Hbar} \antisymmetric{p_0}{x_0}
+ \inv{2!} \lr{\frac{i a}{\Hbar}}^2 \antisymmetric{p_0}{\antisymmetric{p_0}{x_0}}
+ \cdots \\
&=
x_0
+ \frac{i a}{\Hbar} \lr{ -i \Hbar }
+ \inv{2!} \lr{\frac{i a}{\Hbar}}^2 \antisymmetric{p_0}{-i \Hbar}
+ \cdots \\
&=
x_0 + a.
\end{aligned}
\end{equation}

The position expectation with respect to this translated state is

\begin{equation}\label{eqn:translationExpectation:100}
\begin{aligned}
\expectation{x}
&= \bra{0} \lr{ (x_0 + a)\cos(\omega t) + \frac{p_0}{m \omega} \sin(\omega t)
}\ket{0} \\
&= a \cos(\omega t).
\end{aligned}
\end{equation}

The final simplification above follows from \( \bra{n} x \ket{n} = \bra{n} p \ket{n} = 0 \).

References

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

[2] Wikipedia. Baker-campbell-hausdorff formula — wikipedia, the free encyclopedia, 2015. URL https://en.wikipedia.org/w/index.php?title=Baker\%E2\%80\%93Campbell\%E2\%80\%93Hausdorff_formula&oldid=665123858. [Online; accessed 16-August-2015].

Update to old phy356 (Quantum Mechanics I) notes.

February 12, 2015 math and physics play , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , ,

It’s been a long time since I took QM I. My notes from that class were pretty rough, but I’ve cleaned them up a bit.

The main value to these notes is that I worked a number of introductory Quantum Mechanics problems.

These were my personal lecture notes for the Fall 2010, University of Toronto Quantum mechanics I course (PHY356H1F), taught by Prof. Vatche Deyirmenjian.

The official description of this course was:

The general structure of wave mechanics; eigenfunctions and eigenvalues; operators; orbital angular momentum; spherical harmonics; central potential; separation of variables, hydrogen atom; Dirac notation; operator methods; harmonic oscillator and spin.

This document contains a few things

• My lecture notes.
Typos, if any, are probably mine(Peeter), and no claim nor attempt of spelling or grammar correctness will be made. The first four lectures had chosen not to take notes for since they followed the text very closely.
• Notes from reading of the text. This includes observations, notes on what seem like errors, and some solved problems. None of these problems have been graded. Note that my informal errata sheet for the text has been separated out from this document.
• Some assigned problems. I have corrected some the errors after receiving grading feedback, and where I have not done so I at least recorded some of the grading comments as a reference.
• Some worked problems associated with exam preparation.