HW2

Bogliubov transformation

a.) Consider a Hamiltonian of the form Ĥ=D12(ââ+ââ)+D22(ââ+ââ)+D3\hat{H}=\frac{D_{1}}{2}\left(\hat{a}^{\dagger} \hat{a}^{\dagger}+\hat{a} \hat{a}\right)+\frac{D_{2}}{2}\left(\hat{a}^{\dagger} \hat{a}+\hat{a} \hat{a}^{\dagger}\right)+D_{3} with D2>D1>0D_{2}>D_{1}>0. Here â\hat{a} and â\hat{a}^{\dagger} are simple harmonic oscillator raising and lowering operators with [â1â]=1\left[\hat{a}_{1} \hat{a}^{\dagger}\right]=1. Determine an analytic formula for the energy levels in terms of D1,D2,D3D_{1}, D_{2}, D_{3} and integers.

Hint: Consider a new set of raising and lowering operators Â=âcoshθ+âsinhθ\hat{A}=\hat{a} \cosh \theta+\hat{a}^{\dagger} \sinh \theta and Â=âcoshθ+âsinhθ\hat{A}^{\dagger}=\hat{a}^{\dagger} \cosh \theta+\hat{a} \sinh \theta. Verify that [Â,Â]=1\left[\hat{A}, \hat{A}^{\dagger}\right]=1 and then find a way to pick a θ\theta such that Ĥ=C1ÂÂ+C2\hat{H}=C_{1} \hat{A}^{\dagger} \hat{A}+C_{2}. From this form, you should be able to read off the spectrum.

b.) Consider our squeezed state Ŝ(ξ,η)|0.\hat{S}(\xi,\eta) |0\rangle. Show that there is a linear combination Â=âcoshϕ+âsinhϕ\hat{A}=\hat{a} \cosh \phi+\hat{a}^{\dagger} \sinh \phi for some ϕ\phi such that Â\hat{A} annihilates the squeezed state: ÂŜ(ξ,η)|0=0\hat{A} \hat{S}(\xi, \eta)|0\rangle=0

This means that we can think of the squeezed vacuum as the ground state of a Hamiltonian of the form C1ÂÂ+C2C_{1} \hat{A}^{\dagger} \hat{A}+C_{2}. You need to find ϕ\phi as a function of ξ\xi and η\eta.
Note that ϕ\phi is generically complex and you can write the answer in terms of inverse functions.

Simple harmonic oscillator wavefunction in momentum space

a.) Repeat the derivation of the wavefunctions for the simple harmonic oscillator, but now in momentum space. First, verify that |p=eipx̂|p=0|p\rangle=e^{i \frac{p \hat{x}}{\hbar}} | p{=}0\rangle.

Second, define ϕn(p)=(i)np|n=(i)n1n!p=0|eiipx̂(â)n|0.\begin{aligned} \phi_{n}(p) & =(i)^{n}\langle p | n\rangle\\ & =(i)^{n} \frac{1}{\sqrt{n!}}\langle p{=}0| e^{-\frac{i}{\hbar}i p \hat{x}}\left(\hat{a}^{\dagger}\right)^{n}|0\rangle. \end{aligned} Use operator methods to find the wavefunction, which looks schematically similar to (polynomial in p) times exp (polynomial in p)

b.) When computing the wavefunction in position space, we argued we needed to convert the p̂\hat{p} operator in the exponent into an x̂\hat{x} operator, so it can annihilate against x=0|\langle x{=}0|. We did this by breaking the exp(ixp̂)\exp \left(\tfrac{i}{\hbar}x \hat{p}\right) into an exp(â\hat{a}^{\dagger}) and exp(â\hat{a}) factors, moved the exp(â\hat{a}) factor to the right where it annihilated against |0|0\rangle. Then introduced a new exp(â)(\hat{a}) factor and moved to the left, finally combining with exp(â)\exp (\hat{a}^{\dagger})to get exp(x̂)\exp (\hat{x}). We can shorten the derivation by introducing the correct exp(â)\exp (\hat{a}) factor on the right acting on |0|0\rangle (multiply by one trick) move it to the left and combine with the exp(p̂\hat{p}) factor to make an exp(x̂\hat{x}) factor. Find the correct exp(αâ)\exp (\alpha \hat{a}) factor to introduce on the right and show the steps needed to verify that

ψn(x)=1n!emw0x22x=0|(a+2mw0)n|0\psi_{n}(x)=\frac{1}{\sqrt{n!}} e^{-\frac{m w_{0} x^{2}}{2 \hbar}}\langle x{=}0|\left(a^{\dagger}+\sqrt{\frac{2 m w_{0}}{\hbar}}\right)^{n}|0\rangle en route to finding the position wavefunction.

Displacement operator

Consider an arbitrary position translation followed by a momentum translation: exp(ip0x̂)exp(ix0p̂).\exp \left(\frac{i}{\hbar} p_{0} \hat{x}\right) \exp \left(-\frac{i}{\hbar} x_{0} \hat{p}\right). Combine both operators into one single exponent and replace x̂\hat{x} and p̂\hat{p} in terms of the â\hat{a} and â\hat{a}^{\dagger} operators to rewrite the operator as D(α)eiϕ.D(\alpha) e^{i \phi}. Determine α\alpha and ϕ\phi in terns of x0x_{0} and p0p_{0}. Since the overall factor eiϕe^{i \phi} plays no role in wavefunctions, we can drop it when constructing the coherent state D(α)|0=|αD(\alpha)|0\rangle=|\alpha\rangle

How does the state change if we translate momentum first and then translate position?

Uncertainty in coherent and squeezed states

a.) When examining the general coherent states as a function of time, we found eiĤtD(α)|0=eiω0t2|αeiω0t=eiω0t2D(αeiω0t)|0.\begin{aligned} e^{-i \hat{H} t} D(\alpha)|0\rangle & =e^{-i \frac{\omega_{0} t}{2}}\left|\alpha e^{-i \omega_{0} t}\right\rangle \\ & =e^{\frac{-i \omega_{0} t}{2}} D\left(\alpha e^{-i \omega_{0} t}\right)|0\rangle. \end{aligned} Compute the expectation value of x̂\hat{x} and p̂\hat{p} as functions of time along with (Δx)αeiω0t(\Delta x)_{\alpha e^{- i \omega_{0} t}} and (Δp)αeiω0t(\Delta p)_{\alpha e^{- i \omega_{0} t}}. Show that the system is always in a minimum uncertainty state. Explain how the uncertainty in position and momentum change with time. Express your results in terms of x0x_{0} and p0p_{0}, using the α\alpha you found in problem 3.

Recall: (ΔÔ)2ψ=ψ|Ô2|ψψ|O|ψ2(\Delta \hat{O})^{2} \psi =\langle\psi| \hat{O}^{2}|\psi\rangle -\langle\psi| O|\psi\rangle^{2}
b.) For the squeezed vacuum (not the displaced squeezed vacuum), we saw that

eiĤtŜ(ξ,η)(0)=eiω0t2Ŝ(ξe2iω0t,η)|0=|ξe2iω0t,η\begin{aligned} e^{-i \hat{H} t} \hat{S}(\xi, \eta)(0) & =e^{-i \frac{\omega_{0} t}{2}} \hat{S}\left(\xi e^{-2i \omega_{0} t}, \eta)|0\rangle \right. \\ & =\left|\xi e^{-2i \omega_{0} t} , \eta \right \rangle \end{aligned} Pick ξ=reiϕ\xi=r e^{i \phi} and η=0\eta=0 and determine the expectation value of x̂\hat{x} and p̂\hat{p} as functions of time along with (Δx)ξeiω0t,η=0(\Delta x)_{\xi e^{-i\omega_{0}t}, \eta =0} and (Δp)ξeiω0t,η=0(\Delta p)_{\xi e^{-i\omega_{0}t}, \eta =0}. Show that ΔxΔp=/2\Delta x \Delta p = \hbar/2. Explain how Δx\Delta x and Δp\Delta p vary with time.

Factorization method for the particle in a box

a.) In class, we showed the original Schrödinger factorization method for a particle in an infinite square well. Schrödinger described this as "shooting sparrows with artillery". We can proceed in another fashion.

Take the potential to be zero between L2xL2-\frac{L}{2} \leq x \leq \frac{L}{2}. Consider the lowering operator Âk=12m(p̂iktan(kx̂)).\hat{A}_{k}=\frac{1}{\sqrt{2} m}\big (\hat{p}-i \hbar k \tan (k \hat{x})\big ). Show that Ĥ=ÂkÂk+Ek\hat{H}=\hat{A}_{k}^{\dagger} \hat{A}_{k}^{\phantom{\dagger}}+E_{k}, where you need to determine EkE_{k}.
Now, consider increasing kk. Starting from k=0k = 0, we see that EkE_{k} increases until kk reaches π2\frac{\pi}{2}. This is the same solution we examined in class. But now, for excited states, instead of using the Sctrödinger factorization method again, lets just consider increasing kk further. The EkE_k continues to increase, but we will find that when tankL2=\tan k \frac{L}{2}=\infty again, we find another excited state and so on. The idea is that we increase kk until each time tankL2\tan k \frac{L}{2} diverges. This condition coincides with ψ(±L2)=0\psi\left( \pm \frac{L}{2}\right)=0.

Verify that the energies and wavefunctions are given by the well-known results for the particle in a box.

b.) Now consider a potential that is finite V(x)={V0|x|L20|x|>L2V(x)=\left\{\begin{array}{cc} -V_{0} & |x| \leqslant \frac{L}{2} \\ 0 & |x|>\frac{L}{2} \end{array}\right. for |x|L2|x| \leqslant \frac{L}{2}. Use Ak=12(p̂iktan(kx̂))A_{k}=\frac{1}{\sqrt{2}}(\hat{p}-i \hbar k \tan (k \hat{x})) and
for |x|12|x| \geqslant \frac{1}{2} use Ak=12m(p̂±iK)A_{k}=\frac{1}{\sqrt{2 m}}(\hat{p} \pm i \hbar K) (decide whether + or - for x<L2x < -\frac{L}{2} and x>L2x > \frac{L}{2}).
Let ϕ=kL2=2m(V0+E)2L2,κ=2mE2\phi=\frac{k L}{2}=\sqrt{\frac{2 m\left(V_{0}+E\right)}{\hbar^{2}}} \frac{L}{2}, \quad \kappa=\sqrt{-\frac{2 m E}{\hbar^{2}}} and ϕ0=2mV02L2\phi_{0}=\sqrt{\frac{2 m V_{0}}{\hbar^{2}}} \frac{L}{2} (recall E<0 for bound states). Use the requirement that Âk\hat{A}_{k} is continuous at x=±L2x= \pm \frac{L}{2}, to find a transcendental equation that determines a valid solution (this requirement comes from conservation of probability current). Note that any kk value that satisfies this equation yields a valid solution.

Determine the wavefunctions (unnormalized).

This yields all of the even solutions. One can also find the odd ones but I won’t ask you to.