Quantum Mechanics Term 1 Problems Class 4

Dr A. Donos
(Michaelmas 2021)

Consider the 1-D Simple Harmonic Oscillator with frequency w, If a^ is the annihilation operator
a) Construct the eigenstates |z⟩
b) Normalise it
c) Find its time evolution
d) Describe the time evolution of the corresponding wavefunction in position space.

Hint: Baker-Campell-Haussdorff formula

eX^⁢eY^=eY^+[X^,Y^]+12⁢[X^,[X^,Y^]]+13!⁢[X^,[X^,[X^,Y^]]]+…⁢eX^ (1)

Other facts of the 1-D SHO

H^ =ℏ⁢w⁢(a^+⁢a^+12)
|n⟩ =(a^+)nn!⁢|0⟩
a^ =12⁢m⁢w⁢ℏ⁢(m⁢w⁢x^+i⁢p^)

a) We need to solve a^⁢|z⟩=z⁢|z⟩. An important thing to note is that since a^ is not hermitian, its eigenvalues does not have to be real. We also choose to write |z⟩ in terms of the eigenbasis of the SHO.

|z⟩ =∑n=0∞cn⁢(z)⁢|n⟩
z⁢|z⟩ =∑n=0∞z⁢cn⁢(z)⁢|n⟩
a^⁢|z⟩ =∑n=0∞cn⁢(z)⁢a^⁢|n⟩
=∑n=0∞cn⁢(z)n!⁢a^⁢(a^+)n⁢|0⟩
=∑n=0∞cn⁢(z)n!⁢((a^+)n⁢a^+[a^,(a^+)n])⁢|0⟩
=∑n=0∞cn⁢(z)n!⁢n⁢(a^+)n-1⁢|0⟩
=∑n=1∞n⁢cn⁢(z)⁢|n-1⟩
=∑n=0∞n+1⁢cn+1⁢(z)⁢|n⟩
⇒z⁢cn⁢(z) =n+1⁢cn+1⁢(z)
⇒cn⁢(z) =c⁢znn!
⇒|z⟩ =c⁢∑n=0∞znn!⁢(a^+)nn!⁢|0⟩
=c⁢∑n=0∞znn!⁢(a^+)n⁢|0⟩
=c⁢ez⁢a^+⁢|0⟩

b) We need ⟨z|z⟩=1.

⟨z| =c*⁢ez*⁢a^⁢⟨0|
⟨z|z⟩ =|c|2⁢⟨0|ez*⁢a^⁢ez⁢a^+|0⟩

We also note the following fact

ez*⁢a^⁢|0⟩ =(I^+z^*⁢a^+…)⁢|0⟩=|0⟩
⇒⟨0|⁢ez⁢a^+ =⟨0|

So we use the hint with X^=z*⁢a^,Y^=z⁢a^+.

[X^,Y^] =|z|2⁢[a^,a^+]=|z|2
⇒[X^,[X^,Y^]] =[z*⁢a^,|z|2]=0
ez*⁢a^⁢ez⁢a^+ =ez⁢a^++|z|2⁢ez*⁢a^
=e|z|2⁢ez⁢a^+⁢ez*⁢a^
⇒⟨z|z⟩ =|c|2⁢e|z|2⁢⟨0|ez⁢a^+⁢ez*⁢a^|0⟩
=|c|2⁢e|z|2
⇒c =e-|z|22

c) If |ψ(t=0)⟩=|z⟩,

|ψ⁢(t)⟩ =e-i⁢tℏ⁢H^⁢|z⟩
=c⁢(z)⁢e-i⁢tℏ⁢H^⁢ez⁢a^+⁢|0⟩

We know again from a 1-D SHO,

H^⁢|0⟩ =ℏ⁢w2⁢|0⟩
⇒e-i⁢tℏ⁢H^⁢|0⟩ =e-i⁢w⁢t2⁢|0⟩

So we again use the hint to commute the exponentials. This time X^=-i⁢tℏ⁢H^,Y^=z⁢a^+.

[X^,Y^] =-i⁢tℏ⁢z⁢ℏ⁢w⁢[a^+⁢a^,a^+]
=-i⁢t⁢w⁢z⁢(a^+⁢[a^,a^+]+[a^+,a^+]⁢a^)
=-i⁢t⁢w⁢z⁢a^+
] =-i⁢t⁢w.-i⁢t⁢w⁢z⁢[a^+⁢a^,a^+]
=(-i⁢t⁢w)2⁢z⁢a^+

We make a hypothesis that [X^,[X^,…,[X^,Y^]]⁢…]=(-i⁢t⁢w)n⁢z⁢a^+. We prove this using induction.

[X^,[X^,…,[X^,Y^]]⁢…] =-i⁢tℏ⁢ℏ⁢w⁢(-i⁢t⁢w)n⁢z⁢[a^+⁢a^,a^+]
=(-i⁢t⁢w)n+1⁢z⁢a^+
⇒Y^+[X^,Y^]+12!⁢[X^,[X^,Y^]]+… =z⁢∑n=0∞(-i⁢t⁢w)nn!⁢a^+
=z⁢e-i⁢t⁢w⁢a^+
|ψ⁢(t)⟩ =c⁢(z)⁢ez⁢e-i⁢t⁢w⁢a^+⁢e-i⁢tℏ⁢H^⁢|0⟩
=c⁢(z)⁢e-i⁢t⁢w2⁢ez⁢e-i⁢t⁢w⁢a^+⁢|0⟩
=e-i⁢t⁢w2⁢|z⁢e-i⁢t⁢w⟩

Just as a reminder we are now working in the eigenbasis of z^ as z⁢e-i⁢t⁢w is no longer an integer and hence not in the eigenbasis of k^(H^=ℏwk^).

d) In equation form, the question is really asking for ψz⁢(x)=⟨x|z⟩.

a^⁢|z⟩ =z⁢|z⟩
⇒⟨x|a^|z⟩ =z⁢⟨x|z⟩
⇔12⁢m⁢ℏ⁢w⁢⟨x|m⁢w⁢x^+i⁢p^|z⟩ =z⁢ψz⁢(x)
⇔12⁢m⁢ℏ⁢w⁢(m⁢w⁢x⁢ψz⁢(x)+ℏ⁢ψz′⁢(x)) =z⁢ψz⁢(x)
⇒ℏ2⁢m⁢w⁢ψz′⁢(x) =(z-m⁢w2⁢ℏ⁢x)⁢ψz⁢(x)
⇒ℏ2⁢m⁢w⁢ψz′⁢(x)ψz⁢(x) =z-m⁢w2⁢ℏ⁢x
⇔dd⁢x⁢ln⁡(ψz⁢(x)) =2⁢m⁢wℏ⁢dd⁢x⁢(z⁢x-12⁢m⁢w2⁢ℏ⁢x2)
⇒ln⁡(ψz⁢(x)) =2⁢m⁢wℏ⁢(z⁢x-12⁢m⁢w2⁢ℏ⁢x2)+λ
⇒ψz⁢(x) =d⁢e2⁢m⁢wℏ⁢(z⁢x-12⁢m⁢w2⁢ℏ⁢x2)
=d⁢exp⁡(2⁢m⁢wℏ⁢z⁢x)⁢exp⁡(-m⁢w2⁢ℏ⁢x2)
=d⁢exp⁡(-m⁢w2⁢ℏ⁢(x2-8⁢ℏm⁢w⁢z⁢x))
=d⁢exp⁡(-m⁢w2⁢ℏ⁢((x-2⁢ℏm⁢w⁢z)2-2⁢ℏm⁢w⁢z2))
=d⁢exp⁡(z2)⁢exp⁡(-m⁢w2⁢ℏ⁢(x-2⁢ℏm⁢w⁢z)2)
|ψz⁢(x)|2 =|d|2⁢exp⁡(2⁢m⁢wℏ⁢(z+z*)⁢x)⁢exp⁡(-m⁢wℏ⁢x2)
=|d|2⁢exp⁡(2⁢2⁢m⁢wℏ⁢ℜ⁡(z)⁢x)⁢exp⁡(-m⁢wℏ⁢x2)
=|d|2⁢exp⁡(-m⁢wℏ⁢(-2⁢2⁢ℏm⁢w⁢ℜ⁡(z)⁢x+x2))
=|d|2⁢exp⁡(-m⁢wℏ⁢((x-2⁢ℏm⁢w⁢ℜ⁡(z))2-2⁢ℏm⁢w⁢(ℜ⁡(z))2))
=|d|2⁢exp⁡(2⁢(ℜ⁡(z))2)⁢exp⁡(-m⁢wℏ⁢(x-2⁢ℏm⁢w⁢ℜ⁡(z))2)

However, we must remember that our |z⟩ is really e-i⁢t⁢w2⁢|e-i⁢t⁢w⁢z0⟩. Hence,

ℜ⁡(z)=ℜ⁡(z0)⁢cos⁡(w⁢t)+ℑ⁡(z0)⁢sin⁡(w⁢t)

This shows that the time evolution corresponding to the wavefunction is a Gaussian which does not change form but oscillates with frequency w.