Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Problem set 5

These are my solutions to the problem sets assigned to me during the graduate-level quantum mechanics course.

Problem 1

The operator D(z)=exp⁡[za†−z∗a]D\left(z\right)=\exp\left[za^{\dagger}-z^{*}a\right], for a harmonic oscillator or a single field mode, is a displacement operator in phase space. This is apparent from the relations

D(q0,p0)QD−1(q0,p0)=Q−q0ID(q0,p0)PD−1(q0,p0)=P−p0I\begin{align*} D\left(q_{0},p_{0}\right)QD^{-1}\left(q_{0},p_{0}\right) & =Q-q_{0}\mathbb{I}\\ D\left(q_{0},p_{0}\right)PD^{-1}\left(q_{0},p_{0}\right) & =P-p_{0}\mathbb{I} \end{align*}

where D(z)=exp⁡[i(p0Q−q0P)/ℏ]D\left(z\right)=\exp\left[i\left(p_{0}Q-q_{0}P\right)/\hbar\right]. Determine the composition law for the product of two successive displacements in phase space, D(z1)D(z1)D\left(z_{1}\right)D\left(z_{1}\right). Why is the result not simply equal to D(z1+z2)D\left(z_{1}+z_{2}\right)?

Temos o operador D(z)=exp⁡(za†−z∗a)D\left(z\right)=\exp\left(za^{\dagger}-z^{*}a\right) e queremos calcular D(z1)D(z2)D\left(z_{1}\right)D\left(z_{2}\right). Como temos [a,a†]=I\left[a,a^{\dagger}\right]=\mathbb{I} e de acordo com a fórmula de Zasenhaus (uma fórmula de uma popular série de fórmulas conhecidas pelo nome genérico de fórmulas Baker-Campbell-Hausdorff - BCH): Lima (2013):

et(A+B)=etAetBe−t22![A,B]et33!(2[B,[A,B]]+[A,[A,B]])…e^{t\left(A+B\right)}=e^{tA}e^{tB}e^{-\frac{t^{2}}{2!}\left[A,B\right]}e^{\frac{t^{3}}{3!}\left(2\left[B,\left[A,B\right]\right]+\left[A,\left[A,B\right]\right]\right)}\dots

Logo:

D(z1+z2)=exp⁡((z1+z2)a†−(z1+z2)∗a)=exp⁡((z1+z2)a†−(z1∗+z2∗)a)=exp⁡([z1a†−z1∗a]+[z2a†−z2∗a])\begin{align*} D\left(z_{1}+z_{2}\right) & =\exp\left(\left(z_{1}+z_{2}\right)a^{\dagger}-\left(z_{1}+z_{2}\right)^{*}a\right)\\ & =\exp\left(\left(z_{1}+z_{2}\right)a^{\dagger}-\left(z^{*}_{1}+z^{*}_{2}\right)a\right)\\ & =\exp\left(\left[z_{1}a^{\dagger}-z^{*}_{1}a\right]+\left[z_{2}a^{\dagger}-z^{*}_{2}a\right]\right) \end{align*}

Sendo A=z1a†−z1∗aA=z_{1}a^{\dagger}-z^{*}_{1}a e B=z2a†−z2∗aB=z_{2}a^{\dagger}-z^{*}_{2}a , calculando o comutador:

[z1a†−z1∗a,z2a†−z2∗a]=(z1a†−z1∗a)(z2a†−z2∗a)−(z2a†−z2∗a)(z1a†−z1∗a)[A,B]=(z2z1a†a†−z2∗z1a†a−z2z1∗aa†+z1∗az2∗a)−(z1z2a†a†−z1∗z2a†a−z2∗az1a†+z2∗az1∗a)=−z2∗z1a†a−z2z1∗aa†+z1∗z2a†a+z2∗az1a†=(z1∗z2−z2∗z1)a†a+(z2∗z1−z2z1∗)aa†=(z2∗z1−z2z1∗)aa†−(z2∗z1−z1∗z2)a†a=(z2∗z1−z2z1∗)(aa†−a†a)=(z2∗z1−z2z1∗)[a,a†]=z2∗z1−z2z1∗\begin{align*} \left[z_{1}a^{\dagger}-z^{*}_{1}a,z_{2}a^{\dagger}-z^{*}_{2}a\right] & =\left(z_{1}a^{\dagger}-z^{*}_{1}a\right)\left(z_{2}a^{\dagger}-z^{*}_{2}a\right)-\left(z_{2}a^{\dagger}-z^{*}_{2}a\right)\left(z_{1}a^{\dagger}-z^{*}_{1}a\right)\\ \left[A,B\right] & =\left(z_{2}z_{1}a^{\dagger}a^{\dagger}-z^{*}_{2}z_{1}a^{\dagger}a-z_{2}z^{*}_{1}aa^{\dagger}+z^{*}_{1}az^{*}_{2}a\right)-\left(z_{1}z_{2}a^{\dagger}a^{\dagger}-z^{*}_{1}z_{2}a^{\dagger}a-z^{*}_{2}az_{1}a^{\dagger}+z^{*}_{2}az^{*}_{1}a\right)\\ & =-z^{*}_{2}z_{1}a^{\dagger}a-z_{2}z^{*}_{1}aa^{\dagger}+z^{*}_{1}z_{2}a^{\dagger}a+z^{*}_{2}az_{1}a^{\dagger}\\ & =\left(z^{*}_{1}z_{2}-z^{*}_{2}z_{1}\right)a^{\dagger}a+\left(z^{*}_{2}z_{1}-z_{2}z^{*}_{1}\right)aa^{\dagger}\\ & =\left(z^{*}_{2}z_{1}-z_{2}z^{*}_{1}\right)aa^{\dagger}-\left(z^{*}_{2}z_{1}-z^{*}_{1}z_{2}\right)a^{\dagger}a\\ & =\left(z^{*}_{2}z_{1}-z_{2}z^{*}_{1}\right)\left(aa^{\dagger}-a^{\dagger}a\right)\\ & =\left(z^{*}_{2}z_{1}-z_{2}z^{*}_{1}\right)\left[a,a^{\dagger}\right]\\ & =z^{*}_{2}z_{1}-z_{2}z^{*}_{1} \end{align*}

Então [B,[A,B]]=[A,[A,B]]=0\left[B,\left[A,B\right]\right]=\left[A,\left[A,B\right]\right]=0, de forma que ficamos apenas com:

D(z1+z2)=e(z1a†−z1∗a)a†e(z2a†−z2∗a)ae−12(z2∗z1−z2z1∗)\begin{align*} D\left(z_{1}+z_{2}\right) & =e^{\left(z_{1}a^{\dagger}-z^{*}_{1}a\right)a^{\dagger}}e^{\left(z_{2}a^{\dagger}-z^{*}_{2}a\right)a}e^{-\frac{1}{2}\left(z^{*}_{2}z_{1}-z_{2}z^{*}_{1}\right)} \end{align*}

E pela definiçãoD(z)=exp⁡(za†−z∗a)D\left(z\right)=\exp\left(za^{\dagger}-z^{*}a\right), logo:

D(z1+z2)=e12(z2z1∗−z2∗z1)D(z1)D(z2)D\left(z_{1}+z_{2}\right)=e^{\frac{1}{2}\left(z_{2}z^{*}_{1}-z^{*}_{2}z_{1}\right)}D\left(z_{1}\right)D\left(z_{2}\right)

Abrindo os números imaginários zi=ui+iviz_{i}=u_{i}+iv_{i}, então:

z2z1∗−z2∗z1=(u2+iv2)(u1−iv1)−(u2−iv2)(u1+iv1)=(u2u1−iu2v1+iv2u1+v2v1)−(u2u1+iu2v1−iv2u1+v1v2)=2i(v2u1−u2v1)\begin{align*} z_{2}z^{*}_{1}-z^{*}_{2}z_{1} & =\left(u_{2}+iv_{2}\right)\left(u_{1}-iv_{1}\right)-\left(u_{2}-iv_{2}\right)\left(u_{1}+iv_{1}\right)\\ & =\left(u_{2}u_{1}-iu_{2}v_{1}+iv_{2}u_{1}+v_{2}v_{1}\right)-\left(u_{2}u_{1}+iu_{2}v_{1}-iv_{2}u_{1}+v_{1}v_{2}\right)\\ & =2i\left(v_{2}u_{1}-u_{2}v_{1}\right) \end{align*}

De orma que:

D(z1+z2)=ei(v2u1−u2v1)D(z1)D(z2)D\left(z_{1}+z_{2}\right)=e^{i\left(v_{2}u_{1}-u_{2}v_{1}\right)}D\left(z_{1}\right)D\left(z_{2}\right)

Logo dois doslocamentos na sequência é quase a mesma coisa que fazer um único deslocamento, a única diferença é uma fase. Outra forma de resolver esta questão é usando a fórmula:

et(A+B)=etAetBe−t22![A,B]et33!(2[B,[A,B]]+[A,[A,B]])…e^{t\left(A+B\right)}=e^{tA}e^{tB}e^{-\frac{t^{2}}{2!}\left[A,B\right]}e^{\frac{t^{3}}{3!}\left(2\left[B,\left[A,B\right]\right]+\left[A,\left[A,B\right]\right]\right)}\dots

Então:

D(z)=exp⁡([za†]+[−z∗a])=eza†e−z∗ae−12[za†,−z∗a]e13!(2[−z∗a,[za†,−z∗a]]+[za†,[za†,−z∗a]])…\begin{align*} D\left(z\right) & =\exp\left(\left[za^{\dagger}\right]+\left[-z^{*}a\right]\right)\\ & =e^{za^{\dagger}}e^{-z^{*}a}e^{-\frac{1}{2}\left[za^{\dagger},-z^{*}a\right]}e^{\frac{1}{3!}\left(2\left[-z^{*}a,\left[za^{\dagger},-z^{*}a\right]\right]+\left[za^{\dagger},\left[za^{\dagger},-z^{*}a\right]\right]\right)}\dots \end{align*}

E como:

[za†,−z∗a]=za†(−z∗a)−(−z∗a)za†=∣z∣2[a,a†]=∣z∣2\begin{align*} \left[za^{\dagger},-z^{*}a\right] & =za^{\dagger}\left(-z^{*}a\right)-\left(-z^{*}a\right)za^{\dagger}=\left|z\right|^{2}\left[a,a^{\dagger}\right]=\left|z\right|^{2} \end{align*}

Ficamos com:

D(z)=eza†e−z∗ae−12∣z∣2e13!(2[−z∗a,∣z∣2]+[za†,∣z∣2])…=eza†e−z∗ae−∣z∣2/2\begin{align*} D\left(z\right) & =e^{za^{\dagger}}e^{-z^{*}a}e^{-\frac{1}{2}\left|z\right|^{2}}e^{\frac{1}{3!}\left(2\left[-z^{*}a,\left|z\right|^{2}\right]+\left[za^{\dagger},\left|z\right|^{2}\right]\right)}\dots\\ & =e^{za^{\dagger}}e^{-z^{*}a}e^{-\left|z\right|^{2}/2} \end{align*}

Ainda podemos fazer:

D(z)=exp⁡([−z∗a]+[za†])=e−z∗aeza†e−12[−z∗a,za†]e13(2[za†,[−z∗a,za†]]+[−z∗a,[−z∗a,za†]])…=e−z∗aeza†e12∣z∣2et33!(2[za†,−∣z∣2]+[−z∗a,−∣z∣2])=e−z∗aeza†e∣z∣2/2\begin{align*} D\left(z\right) & =\exp\left(\left[-z^{*}a\right]+\left[za^{\dagger}\right]\right)\\ & =e^{-z^{*}a}e^{za^{\dagger}}e^{-\frac{1}{2}\left[-z^{*}a,za^{\dagger}\right]}e^{\frac{1}{3}\left(2\left[za^{\dagger},\left[-z^{*}a,za^{\dagger}\right]\right]+\left[-z^{*}a,\left[-z^{*}a,za^{\dagger}\right]\right]\right)}\dots\\ & =e^{-z^{*}a}e^{za^{\dagger}}e^{\frac{1}{2}\left|z\right|^{2}}e^{\frac{t^{3}}{3!}\left(2\left[za^{\dagger},-\left|z\right|^{2}\right]+\left[-z^{*}a,-\left|z\right|^{2}\right]\right)}\\ & =e^{-z^{*}a}e^{za^{\dagger}}e^{\left|z\right|^{2}/2} \end{align*}

Onde usamos o fato de que [A,B]=AB−BA=−(BA−AB)=−[B,A]\left[A,B\right]=AB-BA=-\left(BA-AB\right)=-\left[B,A\right]. Temos então que:

D(z)=eza†e−z∗ae−∣z∣2/2=e−z∗aeza†e∣z∣2/2D\left(z\right)=e^{za^{\dagger}}e^{-z^{*}a}e^{-\left|z\right|^{2}/2}=e^{-z^{*}a}e^{za^{\dagger}}e^{\left|z\right|^{2}/2}

E então:

D(z1+z2)=exp⁡([(z1+z2)a†]+[−(z1+z2)∗a])=exp⁡([(z1+z2)a†]+[−(z1∗+z2∗)a])=e(z1+z2)a†e−(z1∗+z2∗)ae−12[(z1+z2)a†,−(z1∗+z2∗)a]et33!(2[−(z1∗+z2∗)a,[(z1+z2)a†,−(z1∗+z2∗)a]]+[(z1+z2)a†,[(z1+z2)a†,−(z1∗+z2∗)a]])…\begin{align*} D\left(z_{1}+z_{2}\right) & =\exp\left(\left[\left(z_{1}+z_{2}\right)a^{\dagger}\right]+\left[-\left(z_{1}+z_{2}\right)^{*}a\right]\right)\\ & =\exp\left(\left[\left(z_{1}+z_{2}\right)a^{\dagger}\right]+\left[-\left(z^{*}_{1}+z^{*}_{2}\right)a\right]\right)\\ & =e^{\left(z_{1}+z_{2}\right)a^{\dagger}}e^{-\left(z^{*}_{1}+z^{*}_{2}\right)a}e^{-\frac{1}{2}\left[\left(z_{1}+z_{2}\right)a^{\dagger},-\left(z^{*}_{1}+z^{*}_{2}\right)a\right]}e^{\frac{t^{3}}{3!}\left(2\left[-\left(z^{*}_{1}+z^{*}_{2}\right)a,\left[\left(z_{1}+z_{2}\right)a^{\dagger},-\left(z^{*}_{1}+z^{*}_{2}\right)a\right]\right]+\left[\left(z_{1}+z_{2}\right)a^{\dagger},\left[\left(z_{1}+z_{2}\right)a^{\dagger},-\left(z^{*}_{1}+z^{*}_{2}\right)a\right]\right]\right)\dots} \end{align*}

E então:

[(z1+z2)a†,−(z1∗+z2∗)a]=−(z1+z2)a†(z1∗+z2∗)a−(−(z1∗+z2∗)a(z1+z2)a†)=−(z1+z2)(z1∗+z2∗)a†a+(z1∗+z2∗)(z1+z2)aa†=(z1+z2)(z1∗+z2∗)(aa†−a†a)=(z1+z2)(z1∗+z2∗)[a,a†]=(z1+z2)(z1∗+z2∗)\begin{align*} \left[\left(z_{1}+z_{2}\right)a^{\dagger},-\left(z^{*}_{1}+z^{*}_{2}\right)a\right] & =-\left(z_{1}+z_{2}\right)a^{\dagger}\left(z^{*}_{1}+z^{*}_{2}\right)a-\left(-\left(z^{*}_{1}+z^{*}_{2}\right)a\left(z_{1}+z_{2}\right)a^{\dagger}\right)\\ & =-\left(z_{1}+z_{2}\right)\left(z^{*}_{1}+z^{*}_{2}\right)a^{\dagger}a+\left(z^{*}_{1}+z^{*}_{2}\right)\left(z_{1}+z_{2}\right)aa^{\dagger}\\ & =\left(z_{1}+z_{2}\right)\left(z^{*}_{1}+z^{*}_{2}\right)\left(aa^{\dagger}-a^{\dagger}a\right)\\ & =\left(z_{1}+z_{2}\right)\left(z^{*}_{1}+z^{*}_{2}\right)\left[a,a^{\dagger}\right]\\ & =\left(z_{1}+z_{2}\right)\left(z^{*}_{1}+z^{*}_{2}\right) \end{align*}

Dessa forma análogo ao caso anterior temos:

[−(z1∗+z2∗)a,[(z1+z2)a†,−(z1∗+z2∗)a]]=[−(z1∗+z2∗)a,(z1+z2)(z1∗+z2∗)]=0\left[-\left(z^{*}_{1}+z^{*}_{2}\right)a,\left[\left(z_{1}+z_{2}\right)a^{\dagger},-\left(z^{*}_{1}+z^{*}_{2}\right)a\right]\right]=\left[-\left(z^{*}_{1}+z^{*}_{2}\right)a,\left(z_{1}+z_{2}\right)\left(z^{*}_{1}+z^{*}_{2}\right)\right]=0

E então:

D(z1+z2)=e(z1+z2)a†e−(z1∗+z2∗)ae−12(z1+z2)(z1∗+z2∗)=e(z1+z2)a†e−(z1∗+z2∗)ae−12(∣z1∣2+∣z2∣2+z1z2∗+z2z1∗)=e−12(z1z2∗+z2z1∗)e(z1+z2)a†e−(z1∗+z2∗)ae−12(∣z1∣2+∣z2∣2)\begin{align*} D\left(z_{1}+z_{2}\right) & =e^{\left(z_{1}+z_{2}\right)a^{\dagger}}e^{-\left(z^{*}_{1}+z^{*}_{2}\right)a}e^{-\frac{1}{2}\left(z_{1}+z_{2}\right)\left(z^{*}_{1}+z^{*}_{2}\right)}\\ & =e^{\left(z_{1}+z_{2}\right)a^{\dagger}}e^{-\left(z^{*}_{1}+z^{*}_{2}\right)a}e^{-\frac{1}{2}\left(\left|z_{1}\right|^{2}+\left|z_{2}\right|^{2}+z_{1}z^{*}_{2}+z_{2}z^{*}_{1}\right)}\\ & =e^{-\frac{1}{2}\left(z_{1}z^{*}_{2}+z_{2}z^{*}_{1}\right)}e^{\left(z_{1}+z_{2}\right)a^{\dagger}}e^{-\left(z^{*}_{1}+z^{*}_{2}\right)a}e^{-\frac{1}{2}\left(\left|z_{1}\right|^{2}+\left|z_{2}\right|^{2}\right)} \end{align*}

E abrindo novamente:

e(z1+z2)a=ez1aez2ae−12[z1a,z2a]⋯=ez1aez2a\begin{align*} e^{\left(z_{1}+z_{2}\right)a} & =e^{z_{1}a}e^{z_{2}a}e^{-\frac{1}{2}\left[z_{1}a,z_{2}a\right]}\dots=e^{z_{1}a}e^{z_{2}a} \end{align*}

Então:

D(z1+z2)=e−12(z1z2∗+z2z1∗)e(z1+z2)a†e−(z1∗+z2∗)ae−12(∣z1∣2+∣z2∣2)=e−12(z1z2∗+z2z1∗)ez1a†[ez2a†e−z1∗a]e−z2∗ae−12(∣z1∣2+∣z2∣2)\begin{align*} D\left(z_{1}+z_{2}\right) & =e^{-\frac{1}{2}\left(z_{1}z^{*}_{2}+z_{2}z^{*}_{1}\right)}e^{\left(z_{1}+z_{2}\right)a^{\dagger}}e^{-\left(z^{*}_{1}+z^{*}_{2}\right)a}e^{-\frac{1}{2}\left(\left|z_{1}\right|^{2}+\left|z_{2}\right|^{2}\right)}\\ & =e^{-\frac{1}{2}\left(z_{1}z^{*}_{2}+z_{2}z^{*}_{1}\right)}e^{z_{1}a^{\dagger}}\left[e^{z_{2}a^{\dagger}}e^{-z^{*}_{1}a}\right]e^{-z^{*}_{2}a}e^{-\frac{1}{2}\left(\left|z_{1}\right|^{2}+\left|z_{2}\right|^{2}\right)} \end{align*}

Queremos agora inverter ez2a†e−z1∗a→e−z1∗aez2a†e^{z_{2}a^{\dagger}}e^{-z^{*}_{1}a}\rightarrow e^{-z^{*}_{1}a}e^{z_{2}a^{\dagger}}. Usando apenas e(A+B)=eAeBe−12[A,B]e^{\left(A+B\right)}=e^{A}e^{B}e^{-\frac{1}{2}\left[A,B\right]} porque jábemos que os outros termo zeram, então:

e(z2a†−z1∗a)=ez2a†e−z1∗ae−12[z2a†,−z1∗a]=ez2a†e−z1∗ae−z1∗z22e(−z1∗a+z2a†)=e−z1∗aez2a†e−12[−z1∗a,z2a†]=e−z1∗aez2a†ez1∗z22\begin{align*} e^{\left(z_{2}a^{\dagger}-z^{*}_{1}a\right)} & =e^{z_{2}a^{\dagger}}e^{-z^{*}_{1}a}e^{-\frac{1}{2}\left[z_{2}a^{\dagger},-z^{*}_{1}a\right]}=e^{z_{2}a^{\dagger}}e^{-z^{*}_{1}a}e^{-\frac{z^{*}_{1}z_{2}}{2}}\\ e^{\left(-z^{*}_{1}a+z_{2}a^{\dagger}\right)} & =e^{-z^{*}_{1}a}e^{z_{2}a^{\dagger}}e^{-\frac{1}{2}\left[-z^{*}_{1}a,z_{2}a^{\dagger}\right]}=e^{-z^{*}_{1}a}e^{z_{2}a^{\dagger}}e^{\frac{z^{*}_{1}z_{2}}{2}} \end{align*}

Logo:

e−z1∗aez2a†ez1∗z22=ez2a†e−z1∗ae−z1∗z22e−z1∗aez2a†ez1∗z2=ez2a†e−z1∗a\begin{align*} e^{-z^{*}_{1}a}e^{z_{2}a^{\dagger}}e^{\frac{z^{*}_{1}z_{2}}{2}} & =e^{z_{2}a^{\dagger}}e^{-z^{*}_{1}a}e^{-\frac{z^{*}_{1}z_{2}}{2}}\\ e^{-z^{*}_{1}a}e^{z_{2}a^{\dagger}}e^{z^{*}_{1}z_{2}} & =e^{z_{2}a^{\dagger}}e^{-z^{*}_{1}a} \end{align*}

Substituindo então:

D(z1+z2)=e−12(z1z2∗+z2z1∗)ez1a†e−z1∗aez2a†ez1∗z2e−z2∗ae−12(∣z1∣2+∣z2∣2)=e−12(z2z1∗−z1z2∗)ez1a†e−z1∗aez2a†e−z2∗ae−12(∣z1∣2+∣z2∣2)\begin{align*} D\left(z_{1}+z_{2}\right) & =e^{-\frac{1}{2}\left(z_{1}z^{*}_{2}+z_{2}z^{*}_{1}\right)}e^{z_{1}a^{\dagger}}e^{-z^{*}_{1}a}e^{z_{2}a^{\dagger}}e^{z^{*}_{1}z_{2}}e^{-z^{*}_{2}a}e^{-\frac{1}{2}\left(\left|z_{1}\right|^{2}+\left|z_{2}\right|^{2}\right)}\\ & =e^{-\frac{1}{2}\left(z_{2}z^{*}_{1}-z_{1}z^{*}_{2}\right)}e^{z_{1}a^{\dagger}}e^{-z^{*}_{1}a}e^{z_{2}a^{\dagger}}e^{-z^{*}_{2}a}e^{-\frac{1}{2}\left(\left|z_{1}\right|^{2}+\left|z_{2}\right|^{2}\right)} \end{align*}

Então o produto vai ser dado por:

D(z1)D(z2)=ez1a†e−z1∗ae−12∣z1∣2ez2a†e−z2∗ae−12∣z2∣2=ez1a†e−z1∗aez2a†e−z2∗ae−12(∣z1∣2+∣z2∣2)=e12(z1z2∗+z2z1∗)D(z1+z2)\begin{align*} D\left(z_{1}\right)D\left(z_{2}\right) & =e^{z_{1}a^{\dagger}}e^{-z^{*}_{1}a}e^{-\frac{1}{2}\left|z_{1}\right|^{2}}e^{z_{2}a^{\dagger}}e^{-z^{*}_{2}a}e^{-\frac{1}{2}\left|z_{2}\right|^{2}}\\ & =e^{z_{1}a^{\dagger}}e^{-z^{*}_{1}a}e^{z_{2}a^{\dagger}}e^{-z^{*}_{2}a}e^{-\frac{1}{2}\left(\left|z_{1}\right|^{2}+\left|z_{2}\right|^{2}\right)}\\ & =e^{\frac{1}{2}\left(z_{1}z^{*}_{2}+z_{2}z^{*}_{1}\right)}D\left(z_{1}+z_{2}\right) \end{align*}

Problem 2

If the initial state vector of a harmonic oscillator is the coherent state, ∣Ψ(0)⟩=∣z⟩\left|\Psi\left(0\right)\right\rangle =\left|z\right\rangle, show that the state remains coherent, and the time evolution of ∣Ψ(t)⟩\left|\Psi\left(t\right)\right\rangle corresponds to a classical orbit in (q0,p0)\left(q_{0},p_{0}\right) space.

A demonstração que um estado coerente permanece coerente, e que sua evolução no tempo corresponde a uma órbita clássica no espaço (q0,p0)\left(q_{0},p_{0}\right) foi demonstrado na lista 4, mas repetirei o cálculo a seguir. Sendo então ∣α⟩\left|\alpha\right\rangle um estado coerente, é um auto estado do operador a^\widehat{a}, de forma que temos a^∣α⟩=α∣α⟩\widehat{a}\left|\alpha\right\rangle =\alpha\left|\alpha\right\rangle. O operador a^\widehat{a} por sua vez é dado por:

a^=mω2ℏx^+i2ℏmωp^\widehat{a}=\sqrt{\frac{m\omega}{2\hbar}}\widehat{x}+\frac{i}{\sqrt{2\hbar m\omega}}\widehat{p}

Então, se α\alpha é um autovalor complexo, e por normalização temos ⟨α∣a^∣α⟩=α\left\langle \alpha\left|\widehat{a}\right|\alpha\right\rangle =\alpha, então os valores quânticos médios da posição e do momento são respectivamente x0=⟨α∣x^∣α⟩x_{0}=\left\langle \alpha\left|\widehat{x}\right|\alpha\right\rangle e p0=⟨α∣p^∣α⟩p_{0}=\left\langle \alpha\left|\widehat{p}\right|\alpha\right\rangle onde:

α=⟨α∣a^∣α⟩=⟨α∣mω2ℏx^+i2ℏmωp^∣α⟩=mω2ℏ⟨α∣x^∣α⟩+i2ℏmω⟨α∣p^∣α⟩=mω2ℏx0+i2ℏmωp0\begin{align} \alpha & =\left\langle \alpha\left|\widehat{a}\right|\alpha\right\rangle \\ & =\left\langle \alpha\left|\sqrt{\frac{m\omega}{2\hbar}}\widehat{x}+\frac{i}{\sqrt{2\hbar m\omega}}\widehat{p}\right|\alpha\right\rangle \\ & =\sqrt{\frac{m\omega}{2\hbar}}\left\langle \alpha\left|\widehat{x}\right|\alpha\right\rangle +\frac{i}{\sqrt{2\hbar m\omega}}\left\langle \alpha\left|\widehat{p}\right|\alpha\right\rangle \\ & =\sqrt{\frac{m\omega}{2\hbar}}x_{0}+\frac{i}{\sqrt{2\hbar m\omega}}p_{0} \end{align}

Queremos reescrever U^∣α⟩\widehat{U}\left|\alpha\right\rangle, onde U^=U^(t−t0)=exp⁡[−iH(t−t0)ℏ]\widehat{U}=\widehat{U}\left(t-t_{0}\right)=\exp\left[\frac{-iH\left(t-t_{0}\right)}{\hbar}\right]. Calculando então U^†a^U^\widehat{U}^{\dagger}\widehat{a}\widehat{U}:

U^†a^U^=exp⁡[i(t−t0)ℏH^]a^exp⁡[−i(t−t0)ℏH^]=eλA^B^e−λA^\begin{align*} \widehat{U}^{\dagger}\widehat{a}\widehat{U} & =\exp\left[\frac{i\left(t-t_{0}\right)}{\hbar}\widehat{H}\right]\widehat{a}\exp\left[-\frac{i\left(t-t_{0}\right)}{\hbar}\widehat{H}\right]\\ & =e^{\lambda\widehat{A}}\widehat{B}e^{-\lambda\widehat{A}} \end{align*}

Onde A^=H^\widehat{A}=\widehat{H}, B^=a^\widehat{B}=\widehat{a} e λ=i(t−t0)ℏ\lambda=\frac{i\left(t-t_{0}\right)}{\hbar}. Esse produto é dado por outra das fórmulas Baker-Hausdorff, conforme equação:

eλA^B^e−λA^=B^+λ[A^,B^]+λ22![A^,[A^,B^]]+λ33![A^,[A^,[A^,B^]]]+…e^{\lambda\widehat{A}}\widehat{B}e^{-\lambda\widehat{A}}=\widehat{B}+\lambda\left[\widehat{A},\widehat{B}\right]+\frac{\lambda^{2}}{2!}\left[\widehat{A},\left[\widehat{A},\widehat{B}\right]\right]+\frac{\lambda^{3}}{3!}\left[\widehat{A},\left[\widehat{A},\left[\widehat{A},\widehat{B}\right]\right]\right]+\dots

E como:

a^†a^=(mω2ℏx^−i2ℏmωp^)(mω2ℏx^+i2ℏmωp^)=(mω2ℏx^2+i2ℏx^p^−i2ℏp^x^+12ℏmωp^2)=mω2ℏx^2+i2ℏ(x^p^−p^x^)+12ℏmωp^2=mω2ℏx^2+i2ℏ[x^,p^]+12ℏmωp^2\begin{align*} \widehat{a}^{\dagger}\widehat{a} & =\left(\sqrt{\frac{m\omega}{2\hbar}}\widehat{x}-\frac{i}{\sqrt{2\hbar m\omega}}\widehat{p}\right)\left(\sqrt{\frac{m\omega}{2\hbar}}\widehat{x}+\frac{i}{\sqrt{2\hbar m\omega}}\widehat{p}\right)\\ & =\left(\frac{m\omega}{2\hbar}\widehat{x}^{2}+\frac{i}{2\hbar}\widehat{x}\widehat{p}-\frac{i}{2\hbar}\widehat{p}\widehat{x}+\frac{1}{2\hbar m\omega}\widehat{p}^{2}\right)\\ & =\frac{m\omega}{2\hbar}\widehat{x}^{2}+\frac{i}{2\hbar}\left(\widehat{x}\widehat{p}-\widehat{p}\widehat{x}\right)+\frac{1}{2\hbar m\omega}\widehat{p}^{2}\\ & =\frac{m\omega}{2\hbar}\widehat{x}^{2}+\frac{i}{2\hbar}\left[\widehat{x},\widehat{p}\right]+\frac{1}{2\hbar m\omega}\widehat{p}^{2} \end{align*}

Lembrando que [x^,p^]=iℏ\left[\widehat{x},\widehat{p}\right]=i\hbar:

a^†a^=mω2ℏx^2+12ℏmωp^2−12=1ℏω(12mp^2+mω22x^2−ℏω2)\begin{align*} \widehat{a}^{\dagger}\widehat{a} & =\frac{m\omega}{2\hbar}\widehat{x}^{2}+\frac{1}{2\hbar m\omega}\widehat{p}^{2}-\frac{1}{2}\\ & =\frac{1}{\hbar\omega}\left(\frac{1}{2m}\widehat{p}^{2}+\frac{m\omega^{2}}{2}\widehat{x}^{2}-\frac{\hbar\omega}{2}\right) \end{align*}

O Hamiltoniano para um oscilador harmônico pode ser escrito como:

H^=12mp^2+12mω2x^2=ℏω(a^†a^+12)\begin{align*} \widehat{H} & =\frac{1}{2m}\widehat{p}^{2}+\frac{1}{2}m\omega^{2}\widehat{x}^{2}\\ & =\hbar\omega\left(\widehat{a}^{\dagger}\widehat{a}+\frac{1}{2}\right) \end{align*}

Definindo o operador número N^=a^†a^\widehat{N}=\widehat{a}^{\dagger}\widehat{a}, temos que:

H^=ℏω(N^+12)\widehat{H}=\hbar\omega\left(\widehat{N}+\frac{1}{2}\right)

E como e [N^,a^]=−a^\left[\widehat{N},\widehat{a}\right]=-\widehat{a}, então, por simplicidade ignorando ωℏ2\frac{\omega\hbar}{2}:

eλA^B^e−λA^=a^+λℏω[N^,a^]+(λℏω)22![N^,[N^,a^]]+(λℏω)33![N^,[N^,[N^,a^]]]+…=a^+(−λℏω)a^−(λℏω)22![N^,a^]+(−λℏω)33![N^,[N^,a^]]+…=a^+(−λℏω)a^+(λℏω)22!a^−(−λℏω)33![N^,a^]+…=a^+(−λℏω)a^+(λℏω)22!a^+(−λℏω)33!a^+…=a^[1+(−λℏω)+(λℏω)22!+(−λℏω)33!+… ]=a^e−λℏω\begin{align*} e^{\lambda\widehat{A}}\widehat{B}e^{-\lambda\widehat{A}} & =\widehat{a}+\lambda\hbar\omega\left[\widehat{N},\widehat{a}\right]+\frac{\left(\lambda\hbar\omega\right)^{2}}{2!}\left[\widehat{N},\left[\widehat{N},\widehat{a}\right]\right]+\frac{\left(\lambda\hbar\omega\right)^{3}}{3!}\left[\widehat{N},\left[\widehat{N},\left[\widehat{N},\widehat{a}\right]\right]\right]+\dots\\ & =\widehat{a}+\left(-\lambda\hbar\omega\right)\widehat{a}-\frac{\left(\lambda\hbar\omega\right)^{2}}{2!}\left[\widehat{N},\widehat{a}\right]+\frac{\left(-\lambda\hbar\omega\right)^{3}}{3!}\left[\widehat{N},\left[\widehat{N},\widehat{a}\right]\right]+\dots\\ & =\widehat{a}+\left(-\lambda\hbar\omega\right)\widehat{a}+\frac{\left(\lambda\hbar\omega\right)^{2}}{2!}\widehat{a}-\frac{\left(-\lambda\hbar\omega\right)^{3}}{3!}\left[\widehat{N},\widehat{a}\right]+\dots\\ & =\widehat{a}+\left(-\lambda\hbar\omega\right)\widehat{a}+\frac{\left(\lambda\hbar\omega\right)^{2}}{2!}\widehat{a}+\frac{\left(-\lambda\hbar\omega\right)^{3}}{3!}\widehat{a}+\dots\\ & =\widehat{a}\left[1+\left(-\lambda\hbar\omega\right)+\frac{\left(\lambda\hbar\omega\right)^{2}}{2!}+\frac{\left(-\lambda\hbar\omega\right)^{3}}{3!}+\dots\right]\\ & =\widehat{a}e^{-\lambda\hbar\omega} \end{align*}

E como λ=i(t−t0)ℏ\lambda=\frac{i\left(t-t_{0}\right)}{\hbar}, temos então que:

U^†a^U^=a^e−iω(t−t0)\widehat{U}^{\dagger}\widehat{a}\widehat{U}=\widehat{a}e^{-i\omega\left(t-t_{0}\right)}

Aplicando o operador U^\widehat{U}:

a^=U^†a^U^eiω(t−t0)U^a^=U^U^†a^U^eiω(t−t0)=a^U^eiω(t−t0)\begin{align*} \widehat{a} & =\widehat{U}^{\dagger}\widehat{a}\widehat{U}e^{i\omega\left(t-t_{0}\right)}\\ \widehat{U}\widehat{a} & =\widehat{U}\widehat{U}^{\dagger}\widehat{a}\widehat{U}e^{i\omega\left(t-t_{0}\right)}\\ & =\widehat{a}\widehat{U}e^{i\omega\left(t-t_{0}\right)} \end{align*}

Mmanipulando a^∣α⟩=α∣α⟩\widehat{a}\left|\alpha\right\rangle =\alpha\left|\alpha\right\rangle, temos:

a^∣α⟩=α∣α⟩U^a^∣α⟩=αU^∣α⟩a^eiω(t−t0)U^∣α⟩=αU^∣α⟩a^(U^∣α⟩)=(αeiω(t−t0))(U^∣α⟩)a^∣α(t)⟩=α(t)∣α(t)⟩\begin{align*} \widehat{a}\left|\alpha\right\rangle & =\alpha\left|\alpha\right\rangle \\ \widehat{U}\widehat{a}\left|\alpha\right\rangle & =\alpha\widehat{U}\left|\alpha\right\rangle \\ \widehat{a}e^{i\omega\left(t-t_{0}\right)}\widehat{U}\left|\alpha\right\rangle & =\alpha\widehat{U}\left|\alpha\right\rangle \\ \widehat{a}\left(\widehat{U}\left|\alpha\right\rangle \right) & =\left(\alpha e^{i\omega\left(t-t_{0}\right)}\right)\left(\widehat{U}\left|\alpha\right\rangle \right)\\ \widehat{a}\left|\alpha\left(t\right)\right\rangle & =\alpha\left(t\right)\left|\alpha\left(t\right)\right\rangle \end{align*}

Isto significa que o estado coerente permanece um estado coerente durante a evolução do tempo com um autovalor dependente do tempo, que é a primeira coisa demonstração que estávamos buscando. Agora fazendo:

⟨α(t)∣a^∣α(t)⟩=⟨α∣U^†a^U^∣α⟩=⟨α∣a^e−iω(t−t0)∣α⟩=e−iω(t−t0)⟨α∣a^∣α⟩=e−iω(t−t0)α\begin{align} \left\langle \alpha\left(t\right)\left|\widehat{a}\right|\alpha\left(t\right)\right\rangle & =\left\langle \alpha\left|\widehat{U}^{\dagger}\widehat{a}\widehat{U}\right|\alpha\right\rangle \\ & =\left\langle \alpha\left|\widehat{a}e^{-i\omega\left(t-t_{0}\right)}\right|\alpha\right\rangle \\ & =e^{-i\omega\left(t-t_{0}\right)}\left\langle \alpha\left|\widehat{a}\right|\alpha\right\rangle \\ & =e^{-i\omega\left(t-t_{0}\right)}\alpha \end{align}

E usando a definição de a^\widehat{a} temos:

⟨α(t)∣a^∣α(t)⟩=⟨α(t)∣mω2ℏx^+i2ℏmωp^∣α(t)⟩=mω2ℏ⟨α(t)∣x^∣α(t)⟩+i2ℏmω⟨α(t)∣p^∣α(t)⟩=mω2ℏx(t)+i2ℏmωp(t)\begin{align} \left\langle \alpha\left(t\right)\left|\widehat{a}\right|\alpha\left(t\right)\right\rangle & =\left\langle \alpha\left(t\right)\left|\sqrt{\frac{m\omega}{2\hbar}}\widehat{x}+\frac{i}{\sqrt{2\hbar m\omega}}\widehat{p}\right|\alpha\left(t\right)\right\rangle \\ & =\sqrt{\frac{m\omega}{2\hbar}}\left\langle \alpha\left(t\right)\left|\widehat{x}\right|\alpha\left(t\right)\right\rangle +\frac{i}{\sqrt{2\hbar m\omega}}\left\langle \alpha\left(t\right)\left|\widehat{p}\right|\alpha\left(t\right)\right\rangle \\ & =\sqrt{\frac{m\omega}{2\hbar}}x\left(t\right)+\frac{i}{\sqrt{2\hbar m\omega}}p\left(t\right) \end{align}

Igualando então 37 e 38:

α=eiω(t−t0)(mω2ℏx(t)+i2ℏmωp(t))\alpha=e^{i\omega\left(t-t_{0}\right)}\left(\sqrt{\frac{m\omega}{2\hbar}}x\left(t\right)+\frac{i}{\sqrt{2\hbar m\omega}}p\left(t\right)\right)

Para facilitar, vamos denotar τ=t−t0\tau=t-t_{0}. Esta equação é então satisfeita por uma solução correspondente à órbita clássica no espaço (x0,p0)\left(x_{0},p_{0}\right):

x(τ)≡x0cos⁡(ωτ)+p0mωsin⁡(ωτ)p(τ)≡−mωx0sin⁡(ωτ)+p0cos⁡(ωτ)\begin{align*} x\left(\tau\right) & \equiv x_{0}\cos\left(\omega\tau\right)+\frac{p_{0}}{m\omega}\sin\left(\omega\tau\right)\\ p\left(\tau\right) & \equiv-m\omega x_{0}\sin\left(\omega\tau\right)+p_{0}\cos\left(\omega\tau\right) \end{align*}

Pois:

α=eiωτ[mω2ℏ(x0cos⁡(ωτ)+p0mωsin⁡(ωτ))+i2ℏmω(−mωx0sin⁡(ωτ)+p0cos⁡(ωτ))]=eiωτ[mω2ℏx0cos⁡(ωτ)+p02mωsin⁡(ωτ)−imω2ℏx0sin⁡(ωτ)+p02ℏmωcos⁡(ωτ)]=eiωτ[mω2ℏx0(eiωτ+e−iωτ2)+p02mω(eiωτ−e−iωτ2i)−imω2ℏx0(eiωτ−e−iωτ2i)+p02ℏmω(eiωτ+e−iωτ2)]=[mω2ℏx0(e2iωτ+12)−ip02mω(e2iωτ−12)−mω2ℏx0(e2iωτ−12)+p02ℏmω(e2iωτ+12)]=mω2ℏx0+ip02mω\begin{align*} \alpha & =e^{i\omega\tau}\left[\sqrt{\frac{m\omega}{2\hbar}}\left(x_{0}\cos\left(\omega\tau\right)+\frac{p_{0}}{m\omega}\sin\left(\omega\tau\right)\right)+\frac{i}{\sqrt{2\hbar m\omega}}\left(-m\omega x_{0}\sin\left(\omega\tau\right)+p_{0}\cos\left(\omega\tau\right)\right)\right]\\ & =e^{i\omega\tau}\left[\sqrt{\frac{m\omega}{2\hbar}}x_{0}\cos\left(\omega\tau\right)+\frac{p_{0}}{\sqrt{2 m\omega}}\sin\left(\omega\tau\right)-i\sqrt{\frac{m\omega}{2\hbar}}x_{0}\sin\left(\omega\tau\right)+\frac{p_{0}}{\sqrt{2\hbar m\omega}}\cos\left(\omega\tau\right)\right]\\ & =e^{i\omega\tau}\left[\sqrt{\frac{m\omega}{2\hbar}}x_{0}\left(\frac{e^{i\omega\tau}+e^{-i\omega\tau}}{2}\right)+\frac{p_{0}}{\sqrt{2 m\omega}}\left(\frac{e^{i\omega\tau}-e^{-i\omega\tau}}{2i}\right)-i\sqrt{\frac{m\omega}{2\hbar}}x_{0}\left(\frac{e^{i\omega\tau}-e^{-i\omega\tau}}{2i}\right)+\frac{p_{0}}{\sqrt{2\hbar m\omega}}\left(\frac{e^{i\omega\tau}+e^{-i\omega\tau}}{2}\right)\right]\\ & =\left[\sqrt{\frac{m\omega}{2\hbar}}x_{0}\left(\frac{e^{2i\omega\tau}+1}{2}\right)-\frac{ip_{0}}{\sqrt{2 m\omega}}\left(\frac{e^{2i\omega\tau}-1}{2}\right)-\sqrt{\frac{m\omega}{2\hbar}}x_{0}\left(\frac{e^{2i\omega\tau}-1}{2}\right)+\frac{p_{0}}{\sqrt{2\hbar m\omega}}\left(\frac{e^{2i\omega\tau}+1}{2}\right)\right]\\ & =\sqrt{\frac{m\omega}{2\hbar}}x_{0}+\frac{ip_{0}}{\sqrt{2 m\omega}} \end{align*}

Que é o resultado obtido inicialmente na equação 26, e assim demonstramos a segunda parte do exercício.

Problem 3

Consider a harmonic oscillator with frequency ω\omega, mass mm and Hamiltonian

H=(N+12)ℏωH=\left(N+\frac{1}{2}\right)\hbar\omega

where N≡a†aN\equiv a^{\dagger}a with N∣n⟩=n∣n⟩N\left|n\right\rangle =n\left|n\right\rangle and we recall that the lowering and raising operators satisfy

a∣n⟩=n∣n−1⟩a†∣n⟩=n+1∣n+1⟩\begin{align*} a\left|n\right\rangle & =\sqrt{n}\left|n-1\right\rangle \\ a^{\dagger}\left|n\right\rangle & =\sqrt{n+1}\left|n+1\right\rangle \end{align*}

Assuming that at the timet=0t=0 the oscillator is in a coherent state ∣n⟩\left|n\right\rangle defined by

a∣z⟩=z∣z⟩a\left|z\right\rangle =z\left|z\right\rangle

answer:

a) What is the value of ⟨z∣N∣z⟩\left\langle z\left|N\right|z\right\rangle for z=12exp⁡(iπ/4)z=\frac{1}{2}\exp\left(i\pi/4\right), assuming that ∣z⟩\left|z\right\rangle is normalized?

b) Assuming that at t=0t=0 the oscillator is in the ground state ∣0⟩\left|0\right\rangle, determine the state form

at time t=1/10st=1/10s for ω=5πs−1\omega=5\pi s^{-1}.

c) What is the value of cnc_{n} (as a function of nn and zz) so that the coherent state ∣z⟩=∑n=0+∞cn∣n⟩\left|z\right\rangle =\sum^{+\infty}_{n=0}c_{n}\left|n\right\rangle (expanded in the eigenstates ∣n⟩\left|n\right\rangle of the number operator NN) be normalized?

(Remember that ex=∑n=0+∞xn/n!e^{x}=\sum^{+\infty}_{n=0}x^{n}/n!).

d) Use the result of the previous item and calculate the numerical value of ∣⟨z′∣z⟩∣2\left|\left\langle z'|z\right\rangle \right|^{2} for z=12exp⁡(iπ/4)z=\frac{1}{2}\exp\left(i\pi/4\right)

and z′=15exp⁡(iπ/4)z'=\frac{1}{5}\exp\left(i\pi/4\right).

Temos então o seguinte hamiltoniano de um oscilador harmônico com frequência ω\omega

H=(N+12)ℏωH=\left(N+\frac{1}{2}\right)\hbar\omega

Onde N≡a†aN\equiv a^{\dagger}a, N∣n⟩=n∣n⟩N\left|n\right\rangle =n\left|n\right\rangle, e os operadores aa e a†a^{\dagger}são definidos como:

a∣n⟩=n∣n−1⟩a†∣n⟩=n+1∣n+1⟩\begin{align*} a\left|n\right\rangle & =\sqrt{n}\left|n-1\right\rangle \\ a^{\dagger}\left|n\right\rangle & =\sqrt{n+1}\left|n+1\right\rangle \end{align*}

Assumimos também que no tempo t=0t=0 o oscilador está no estado coerente ∣z⟩\left|z\right\rangle definido por:

a∣z⟩=z∣z⟩a\left|z\right\rangle =z\left|z\right\rangle

Queremos então:

a) O valor ⟨z∣N∣z⟩\left\langle z\left|N\right|z\right\rangle de z=12exp⁡(iπ/4)z=\frac{1}{2}\exp\left(i\pi/4\right) assumindo que ∣z⟩\left|z\right\rangle é normalizado. A partir da definições de NNe da definiçao do estado coerente, podemos escrever então:

⟨z∣N∣z⟩=⟨z∣a†a∣z⟩=(a∣z⟩)†a∣z⟩=(z∣z⟩)†z∣z⟩=⟨z∣z∗z∣z⟩=∣z∣2⟨z∣z⟩=∣z∣2\begin{align*} \left\langle z\left|N\right|z\right\rangle & =\left\langle z\left|a^{\dagger}a\right|z\right\rangle \\ & =\left(a\left|z\right\rangle \right)^{\dagger}a\left|z\right\rangle \\ & =\left(z\left|z\right\rangle \right)^{\dagger}z\left|z\right\rangle \\ & =\left\langle z\left|z^{*}z\right|z\right\rangle \\ & =\left|z\right|^{2}\left\langle z|z\right\rangle \\ & =\left|z\right|^{2} \end{align*}

Então sendo z=12exp⁡(iπ/4)z=\frac{1}{2}\exp\left(i\pi/4\right):

⟨z∣N∣z⟩=12exp⁡(iπ/4)12exp⁡(−iπ/4)\begin{align*} \left\langle z\left|N\right|z\right\rangle & =\frac{1}{2}\exp\left(i\pi/4\right)\frac{1}{2}\exp\left(-i\pi/4\right) \end{align*}

Logo temos:

⟨z∣N∣z⟩=14\left\langle z\left|N\right|z\right\rangle =\frac{1}{4}

b) Assumindo que em t=0t=0 o oscilador está no estado fundamental ∣0⟩\left|0\right\rangle queremos determinar o estado no tempo t=1/10t=1/10 para ω=5πs−1\omega=5\pi s^{-1}. Sendo a dinâmica dada por U^=exp⁡[−iHtℏ]\widehat{U}=\exp\left[\frac{-iHt}{\hbar}\right], temos então:

∣ψ(t)⟩=U^∣0⟩=exp⁡[−iHtℏ]∣0⟩=exp⁡[−i(N+12)ℏωtℏ]∣0⟩=exp⁡[−i12ωt]exp⁡[−iNωtℏ]∣0⟩=exp⁡[−i12ωt]exp⁡[−iωtℏ0]∣0⟩=exp⁡[−i12ωt]∣0⟩=[cos⁡(−ωt2)+sin⁡(−ωt2)]∣0⟩=[cos⁡(ωt2)−sin⁡(ωt2)]∣0⟩\begin{align*} \left|\psi\left(t\right)\right\rangle & =\widehat{U}\left|0\right\rangle \\ & =\exp\left[\frac{-iHt}{\hbar}\right]\left|0\right\rangle \\ & =\exp\left[\frac{-i\left(N+\frac{1}{2}\right)\hbar\omega t}{\hbar}\right]\left|0\right\rangle \\ & =\exp\left[-i\frac{1}{2}\omega t\right]\exp\left[\frac{-iN\omega t}{\hbar}\right]\left|0\right\rangle \\ & =\exp\left[-i\frac{1}{2}\omega t\right]\exp\left[\frac{-i\omega t}{\hbar}0\right]\left|0\right\rangle \\ & =\exp\left[-i\frac{1}{2}\omega t\right]\left|0\right\rangle \\ & =\left[\cos\left(-\frac{\omega t}{2}\right)+\sin\left(-\frac{\omega t}{2}\right)\right]\left|0\right\rangle \\ & =\left[\cos\left(\frac{\omega t}{2}\right)-\sin\left(\frac{\omega t}{2}\right)\right]\left|0\right\rangle \end{align*}

Onde usamos que F(A)∣a⟩=F(a)∣a⟩F\left(A\right)\left|a\right\rangle =F\left(a\right)\left|a\right\rangle. Substituindo os valores dados então:

∣ψ(t)⟩=[cos⁡(5π10⋅2)−isin⁡(5π10⋅2)]∣0⟩\begin{align*} \left|\psi\left(t\right)\right\rangle & =\left[\cos\left(\frac{5\pi}{10\cdot2}\right)-i\sin\left(\frac{5\pi}{10\cdot2}\right)\right]\left|0\right\rangle \end{align*}

O estado que buscamos é então:

∣ψ(t)⟩=(1−i2)∣0⟩\begin{align*} \left|\psi\left(t\right)\right\rangle & =\left(\frac{1-i}{\sqrt{2}}\right)\left|0\right\rangle \end{align*}

c) Queremos agora o valor de cnc_{n}(como uma função de nn e zz) de modo que o estado coerente expandido nos autoestados do operador NN,∣z⟩=∑n=0+∞cn∣n⟩\left|z\right\rangle =\sum^{+\infty}_{n=0}c_{n}\left|n\right\rangle seja normalizado. Obsevação: ex=∑n=0+∞xn/n!e^{x}=\sum^{+\infty}_{n=0}x^{n}/n!. Temos então:

∣z⟩=∑n=0+∞⟨n∣z⟩∣n⟩\begin{align*} \left|z\right\rangle & =\sum^{+\infty}_{n=0}\left\langle n|z\right\rangle \left|n\right\rangle \end{align*}

E partindo de a†∣n⟩=n+1∣n+1⟩a^{\dagger}\left|n\right\rangle =\sqrt{n+1}\left|n+1\right\rangle, podemos calcular:

a†∣0⟩=∣1⟩→∣1⟩=a†1∣0⟩a†∣1⟩=2∣2⟩→∣2⟩=a†2∣1⟩=(a†)22⋅1∣0⟩a†∣2⟩=3∣3⟩→∣3⟩=a†3∣2⟩=(a†)33⋅2⋅1∣0⟩⋮∣n⟩=(a†)nn!∣0⟩\begin{align*} a^{\dagger}\left|0\right\rangle & =\left|1\right\rangle & \rightarrow\left|1\right\rangle & =\frac{a^{\dagger}}{\sqrt{1}}\left|0\right\rangle \\ a^{\dagger}\left|1\right\rangle & =\sqrt{2}\left|2\right\rangle & \rightarrow\left|2\right\rangle & =\frac{a^{\dagger}}{\sqrt{2}}\left|1\right\rangle =\frac{\left(a^{\dagger}\right)^{2}}{\sqrt{2\cdot1}}\left|0\right\rangle \\ a^{\dagger}\left|2\right\rangle & =\sqrt{3}\left|3\right\rangle & \rightarrow\left|3\right\rangle & =\frac{a^{\dagger}}{\sqrt{3}}\left|2\right\rangle =\frac{\left(a^{\dagger}\right)^{3}}{\sqrt{3\cdot2\cdot1}}\left|0\right\rangle \\ & & \vdots\\ & & \left|n\right\rangle & =\frac{\left(a^{\dagger}\right)^{n}}{\sqrt{n!}}\left|0\right\rangle \end{align*}

Temos então:

∣z⟩=∑n=0+∞((a†)nn!∣0⟩)†∣z⟩∣n⟩=∑n=0+∞⟨0∣(a)nn!∣z⟩∣n⟩=∑n=0+∞1n!⟨0∣an∣z⟩∣n⟩=∑n=0+∞znn!⟨0∣z⟩∣n⟩\begin{align*} \left|z\right\rangle & =\sum^{+\infty}_{n=0}\left(\frac{\left(a^{\dagger}\right)^{n}}{\sqrt{n!}}\left|0\right\rangle \right)^{\dagger}\left|z\right\rangle \left|n\right\rangle \\ & =\sum^{+\infty}_{n=0}\left\langle 0\right|\frac{\left(a\right)^{n}}{\sqrt{n!}}\left|z\right\rangle \left|n\right\rangle \\ & =\sum^{+\infty}_{n=0}\frac{1}{\sqrt{n!}}\left\langle 0\right|a^{n}\left|z\right\rangle \left|n\right\rangle \\ & =\sum^{+\infty}_{n=0}\frac{z^{n}}{\sqrt{n!}}\left\langle 0|z\right\rangle \left|n\right\rangle \end{align*}

Onde usamos que a∣z⟩=z∣z⟩a\left|z\right\rangle =z\left|z\right\rangle e F(A)∣a⟩=F(a)∣a⟩F\left(A\right)\left|a\right\rangle =F\left(a\right)\left|a\right\rangle novamente. Então uma vez que temos cn=⟨n∣z⟩→c0=⟨0∣z⟩c_{n}=\left\langle n|z\right\rangle \rightarrow c_{0}=\left\langle 0|z\right\rangle, logo:

cn=znn!c0c_{n}=\frac{z^{n}}{\sqrt{n!}}c_{0}

Precisamos então idetnficar c0c_{0} garantindo a normalização. Dessa forma temos:

⟨z∣z⟩=(∑m=0+∞⟨m∣cm∗)(∑n=0+∞cn∣n⟩)1=∑m,n=0+∞cm∗cn⟨m∣n⟩=∑m,n=0+∞cm∗cnδmn=∑n=0+∞∣cn∣2\begin{align*} \left\langle z|z\right\rangle & =\left(\sum^{+\infty}_{m=0}\left\langle m\right|c^{*}_{m}\right)\left(\sum^{+\infty}_{n=0}c_{n}\left|n\right\rangle \right)\\ 1 & =\sum^{+\infty}_{m,n=0}c^{*}_{m}c_{n}\left\langle m|n\right\rangle \\ & =\sum^{+\infty}_{m,n=0}c^{*}_{m}c_{n}\delta_{mn}\\ & =\sum^{+\infty}_{n=0}\left|c_{n}\right|^{2} \end{align*}

Substituindo então:

1=∑n=0+∞∣znn!c0∣2=∣c0∣2∑n=0+∞(∣z∣)nn!=∣c0∣2∑n=0+∞(∣z∣2)nn!=∣c0∣2e∣z∣2\begin{align*} 1 & =\sum^{+\infty}_{n=0}\left|\frac{z^{n}}{\sqrt{n!}}c_{0}\right|^{2}\\ & =\left|c_{0}\right|^{2}\sum^{+\infty}_{n=0}\frac{\left(\left|z\right|\right)^{n}}{n!}\\ & =\left|c_{0}\right|^{2}\sum^{+\infty}_{n=0}\frac{\left(\left|z\right|^{2}\right)^{n}}{n!}\\ & =\left|c_{0}\right|^{2}e^{\left|z\right|^{2}} \end{align*}

Onde usamos que ex=∑n=0+∞xn/n!e^{x}=\sum^{+\infty}_{n=0}x^{n}/n!. Temos:

∣c0∣=e−∣z∣2/2\left|c_{0}\right|=e^{-\left|z\right|^{2}/2}

Assim o valor que buscamos para cnc_{n} é:

cn=znn!e−∣z∣2/2c_{n}=\frac{z^{n}}{\sqrt{n!}}e^{-\left|z\right|^{2}/2}

Uma vez que e−∣z∣2/2e^{-\left|z\right|^{2}/2} é sempre real e maior que 0.

d)

Usando o resultado do item anterior, vamos calcular ∣⟨z′∣z⟩∣2\left|\left\langle z'|z\right\rangle \right|^{2} onde z=12eiπ/4z=\frac{1}{2}e^{i\pi/4} e z′=14eiπ/4z'=\frac{1}{4}e^{i\pi/4} , ete é um cálculo bastante direto:

⟨z′∣z⟩=(∑n′=0+∞⟨n′∣cn′)(∑n=0+∞cn∣n⟩)=∑n,n′=0+∞cn′∗cn⟨n′n⟩=∑n,n′=0+∞((z′n′)∗n′!e−∣z′∣2/2)(znn!e−∣z∣2/2)δn′n=e−(∣z∣2+∣z′∣2)/2∑n=0+∞((z′n′)∗n!)(znn!)=e−(∣z∣2+∣z′∣2)/2∑n=0+∞((z′∗z)nn!)=e−(∣z∣2+∣z′∣2)/2∑n=0+∞((14e−iπ/412eiπ/4)nn!)=e−(∣z∣2+∣z′∣2)/2∑n=0+∞((1/8)nn!)\begin{align*} \left\langle z'|z\right\rangle & =\left(\sum^{+\infty}_{n'=0}\left\langle n'\right|c_{n'}\right)\left(\sum^{+\infty}_{n=0}c_{n}\left|n\right\rangle \right)\\ & =\sum^{+\infty}_{n,n'=0}c^{*}_{n'}c_{n}\left\langle n'n\right\rangle \\ & =\sum^{+\infty}_{n,n'=0}\left(\frac{\left(z'^{n'}\right)^{*}}{\sqrt{n'!}}e^{-\left|z'\right|^{2}/2}\right)\left(\frac{z^{n}}{\sqrt{n!}}e^{-\left|z\right|^{2}/2}\right)\delta_{n'n}\\ & =e^{-\left(\left|z\right|^{2}+\left|z'\right|^{2}\right)/2}\sum^{+\infty}_{n=0}\left(\frac{\left(z'^{n'}\right)^{*}}{\sqrt{n!}}\right)\left(\frac{z^{n}}{\sqrt{n!}}\right)\\ & =e^{-\left(\left|z\right|^{2}+\left|z'\right|^{2}\right)/2}\sum^{+\infty}_{n=0}\left(\frac{\left(z'^{*}z\right){}^{n}}{n!}\right)\\ & =e^{-\left(\left|z\right|^{2}+\left|z'\right|^{2}\right)/2}\sum^{+\infty}_{n=0}\left(\frac{\left(\frac{1}{4}e^{-i\pi/4}\frac{1}{2}e^{i\pi/4}\right){}^{n}}{n!}\right)\\ & =e^{-\left(\left|z\right|^{2}+\left|z'\right|^{2}\right)/2}\sum^{+\infty}_{n=0}\left(\frac{\left(1/8\right){}^{n}}{n!}\right) \end{align*}

Novamente usando ex=∑n=0+∞xn/n!e^{x}=\sum^{+\infty}_{n=0}x^{n}/n!, e sendo ∣z∣2=14,∣z′∣2=116\left|z\right|^{2}=\frac{1}{4},\left|z'\right|^{2}=\frac{1}{16}, ficamos com:

⟨z′∣z⟩=exp⁡[−(14+116)12]exp⁡[18]=exp⁡[−(432+132)+432]=exp⁡[−116]\begin{align*} \left\langle z'|z\right\rangle & =\exp\left[-\left(\frac{1}{4}+\frac{1}{16}\right)\frac{1}{2}\right]\exp\left[\frac{1}{8}\right]\\ & =\exp\left[-\left(\frac{4}{32}+\frac{1}{32}\right)+\frac{4}{32}\right]\\ & =\exp\left[-\frac{1}{16}\right] \end{align*}

Logo:

⟨z′∣z⟩=exp⁡[−(14+116)12]exp⁡[18]=exp⁡[−(432+132)+432]=exp⁡[−132]\begin{align*} \left\langle z'|z\right\rangle & =\exp\left[-\left(\frac{1}{4}+\frac{1}{16}\right)\frac{1}{2}\right]\exp\left[\frac{1}{8}\right]\\ & =\exp\left[-\left(\frac{4}{32}+\frac{1}{32}\right)+\frac{4}{32}\right]\\ & =\exp\left[-\frac{1}{32}\right] \end{align*}

Fazendo o quadrado temos

∣⟨z′∣z⟩∣2=exp⁡[−116]\left|\left\langle z'|z\right\rangle \right|^{2}=\exp\left[-\frac{1}{16}\right]

Ou em valores numéricos ∣⟨z′∣z⟩∣2=0.94\left|\left\langle z'|z\right\rangle \right|^{2}=0.94.

References
  1. Lima, P. C. de. (2013). Equações Diferenciais Parciais I. http://150.164.25.15/ lima/apostilas/EDP1-aulas-22-8-13.pdf