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 4

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

Problem 1

K(x′′,t;x′,t0)=m2πiℏ(t−t0)exp⁡[im(x′′−x′)22ℏ(t−t0)]K\left(x'',t;x',t_{0}\right)=\sqrt{\frac{m}{2\pi i\hbar\left(t-t_{0}\right)}}\exp\left[\frac{im\left(x''-x'\right)^{2}}{2\hbar\left(t-t_{0}\right)}\right]

Show that the propagator of a one-dimensional free particle is given by

Queremos mostrar que o propagador de uma partícula livre unidimensional é dado por:

K(x′′,t;x′,t0)=m2πiℏ(t−t0)exp⁡[im(x′′−x′)22ℏ(t−t0)]K\left(x'',t;x',t_{0}\right)=\sqrt{\frac{m}{2\pi i\hbar\left(t-t_{0}\right)}}\exp\left[\frac{im\left(x''-x'\right)^{2}}{2\hbar\left(t-t_{0}\right)}\right]

O propagador pode ser escrito elegantemente como:

K(x′′,t;x′,t0)=⟨x′′∣exp⁡[−iH(t−t0)ℏ]∣x′⟩K\left(x'',t;x',t_{0}\right)=\left\langle x''\left|\exp\left[\frac{-iH\left(t-t_{0}\right)}{\hbar}\right]\right|x'\right\rangle

Sendo então o hamiltoniano de uma partícula livre dado por H=p22mH=\frac{p^{2}}{2m} Uma vez que ∣p′⟩\left|p'\right\rangle é um autoket do operador pp e HH, logo p∣p′⟩=p′∣p′⟩p\left|p'\right\rangle =p'\left|p'\right\rangle e H∣p′⟩=p′22m∣p′⟩H\left|p'\right\rangle =\frac{p'^{2}}{2m}\left|p'\right\rangle. A autofunção do momento é somente a função transformação:

⟨x′∣p′⟩=12πℏexp⁡(ip′x′ℏ)\left\langle x'|p'\right\rangle =\sqrt{\frac{1}{2\pi\hbar}}\exp\left(\frac{ip'x'}{\hbar}\right)

E como (H∣p′⟩)†=⟨p∣H=⟨p∣p′22m\left(H\left|p'\right\rangle \right)^{\dagger}=\left\langle p\right|H=\left\langle p\right|\frac{p'^{2}}{2m} Então:

K(x′′,t;x′,t0)=⟨x′′∣exp⁡[−iH(t−t0)ℏ]∣x′⟩=∫dp′⟨x′′∣p′⟩⟨p′∣exp⁡[−iH(t−t0)ℏ]∣x′⟩=∫dp′12πℏexp⁡(ip′x′′ℏ)⟨p′∣exp⁡[−iH(t−t0)ℏ]∣x′⟩=12πℏ∫dp′exp⁡(ip′x′′ℏ)⟨p′∣x′⟩exp⁡(−ip′2(t−t0)2ℏm)=12πℏ∫dp′exp⁡(ip′(x′′−x′)ℏ)exp⁡(−ip′2(t−t0)2ℏm)=12πℏ∫dp′exp⁡(ip′(x′′−x′)ℏ−ip′2(t−t0)2ℏm)\begin{align*} K\left(x'',t;x',t_{0}\right) & =\left\langle x''\left|\exp\left[\frac{-iH\left(t-t_{0}\right)}{\hbar}\right]\right|x'\right\rangle \\ & =\int dp'\left\langle x''|p'\right\rangle \left\langle p'\left|\exp\left[\frac{-iH\left(t-t_{0}\right)}{\hbar}\right]\right|x'\right\rangle \\ & =\int dp'\sqrt{\frac{1}{2\pi\hbar}}\exp\left(\frac{ip'x''}{\hbar}\right)\left\langle p'\left|\exp\left[\frac{-iH\left(t-t_{0}\right)}{\hbar}\right]\right|x'\right\rangle \\ & =\sqrt{\frac{1}{2\pi\hbar}}\int dp'\exp\left(\frac{ip'x''}{\hbar}\right)\left\langle p'|x'\right\rangle \exp\left(\frac{-ip'^{2}\left(t-t_{0}\right)}{2\hbar m}\right)\\ & =\frac{1}{2\pi\hbar}\int dp'\exp\left(\frac{ip'\left(x''-x'\right)}{\hbar}\right)\exp\left(\frac{-ip'^{2}\left(t-t_{0}\right)}{2\hbar m}\right)\\ & =\frac{1}{2\pi\hbar}\int dp'\exp\left(\frac{ip'\left(x''-x'\right)}{\hbar}-\frac{ip'^{2}\left(t-t_{0}\right)}{2\hbar m}\right) \end{align*}

Sendo −2αβ=i(x′′−x′)ℏ-2\alpha\beta=\frac{i\left(x''-x'\right)}{\hbar}e β2=i(t′−t0)2mℏ\beta^{2}=\frac{i\left(t'-t_{0}\right)}{2m\hbar}. Completando quadrado:

K(x′′,t;x′,t0)=12πℏ∫dp′exp⁡(−[2αβp′+(βp′)2])=12πℏ∫dp′exp⁡[−(2αβ+βp′)2+α2]=12πℏe−α2∫dp′exp⁡[−(α+βp′)2]=12πℏeα2β∫e−u2du\begin{align*} K\left(x'',t;x',t_{0}\right) & =\frac{1}{2\pi\hbar}\int dp'\exp\left(-\left[2\alpha\beta p'+\left(\beta p'\right)^{2}\right]\right)\\ & =\frac{1}{2\pi\hbar}\int dp'\exp\left[-\left(2\alpha\beta+\beta p'\right)^{2}+\alpha^{2}\right]\\ & =\frac{1}{2\pi\hbar}e^{-\alpha^{2}}\int dp'\exp\left[-\left(\alpha+\beta p'\right)^{2}\right]\\ & =\frac{1}{2\pi\hbar}\frac{e^{\alpha^{2}}}{\beta}\int e^{-u^{2}}du \end{align*}

Onde foi feito u=(α+βp′)u=\left(\alpha+\beta p'\right) e dudp′=−β\frac{du}{dp'}=-\beta. Uma vez que ∫−∞+∞e−ax2dx=πa\int^{+\infty}_{-\infty}e^{-ax^{2}}dx=\sqrt{\frac{\pi}{a}}, então:

K(x′′,t;x′,t0)=12πℏeα2β\begin{align*} K\left(x'',t;x',t_{0}\right) & =\frac{1}{2\sqrt{\pi}\hbar}\frac{e^{\alpha^{2}}}{\beta} \end{align*}

Agora como:

−2αβ=i(x′′−x′)ℏα=−i(x′′−x′)2ℏβα2=(x′′−x′)222ℏ21β2=(x′′−x′)222ℏ2(2imℏ(t′−t0))=im(x′′−x′)22ℏ(t′−t0)\begin{align*} -2\alpha\beta & =\frac{i\left(x''-x'\right)}{\hbar}\\ \alpha & =-\frac{i\left(x''-x'\right)}{2\hbar\beta}\\ \alpha^{2} & =\frac{\left(x''-x'\right)^{2}}{2^{2}\hbar^{2}}\frac{1}{\beta^{2}}\\ & =\frac{\left(x''-x'\right)^{2}}{2^{2}\hbar^{2}}\left(\frac{2im\hbar}{\left(t'-t_{0}\right)}\right)\\ & =\frac{im\left(x''-x'\right)^{2}}{2\hbar\left(t'-t_{0}\right)} \end{align*}

Então:

K(x′′,t;x′,t0)=12πℏeα2β=12πℏ2mℏi(t′−t0)exp⁡[im(x′′−x′)22ℏ(t−t0)]\begin{align*} K\left(x'',t;x',t_{0}\right) & =\frac{1}{2\sqrt{\pi}\hbar}\frac{e^{\alpha^{2}}}{\beta}\\ & =\frac{1}{2\sqrt{\pi}\hbar}\sqrt{\frac{2m\hbar}{i\left(t'-t_{0}\right)}}\exp\left[\frac{im\left(x''-x'\right)^{2}}{2\hbar\left(t-t_{0}\right)}\right] \end{align*}

Então obtemos:

K(x′′,t;x′,t0)=m2πiℏ(t−t0)exp⁡[im(x′′−x′)22ℏ(t−t0)]K\left(x'',t;x',t_{0}\right)=\sqrt{\frac{m}{2\pi i\hbar\left(t-t_{0}\right)}}\exp\left[\frac{im\left(x''-x'\right)^{2}}{2\hbar\left(t-t_{0}\right)}\right]

Que é o que queríamos demonstrar desde o começo.

Problem 2

Show, using the a,a†a,a^{\dagger} operator method, that the propagator of a one-dimensional simple harmonic oscillator is given by

K(x′′,t;x′,t0)=mω2πiℏsin⁡[ωt(t−t0)]exp⁡[imω2ℏsin⁡[ω(t−t0)]((x′′2+x′2)cos⁡[ω(t−t0)]−2x′′x′)]K\left(x'',t;x',t_{0}\right)=\sqrt{\frac{m\omega}{2\pi i\hbar\sin\left[\omega t\left(t-t_{0}\right)\right]}}\exp\left[\frac{im\omega}{2\hbar\sin\left[\omega\left(t-t_{0}\right)\right]}\left(\left(x''^{2}+x'^{2}\right)\cos\left[\omega\left(t-t_{0}\right)\right]-2x''x'\right)\right]

Usando o método do operador a,a†a,a^{\dagger} queremo msotrar que o propagador unidimensional de um oscilador harmônico simples é dado por:

K(x′′,t;x′,t0)=mω2πiℏsin⁡[ωt(t−t0)]exp⁡[imω2ℏsin⁡[ω(t−t0)]((x′′2+x′2)cos⁡[ω(t−t0)]−2x′′x′)]K\left(x'',t;x',t_{0}\right)=\sqrt{\frac{m\omega}{2\pi i\hbar\sin\left[\omega t\left(t-t_{0}\right)\right]}}\exp\left[\frac{im\omega}{2\hbar\sin\left[\omega\left(t-t_{0}\right)\right]}\left(\left(x''^{2}+x'^{2}\right)\cos\left[\omega\left(t-t_{0}\right)\right]-2x''x'\right)\right]

Partindodo propagador escrito como:

K(x′′,t;x′,t0)=⟨x′′∣exp⁡[−iH^(t−t0)ℏ]∣x′⟩K\left(x'',t;x',t_{0}\right)=\left\langle x''\left|\exp\left[\frac{-i\widehat{H}\left(t-t_{0}\right)}{\hbar}\right]\right|x'\right\rangle

E sendo 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], então podemos escrever simplesmente:

K(x′′,t;x′,t0)=⟨x′′∣U^∣x′⟩K\left(x'',t;x',t_{0}\right)=\left\langle x''\left|\widehat{U}\right|x'\right\rangle

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. E o operador a^\widehat{a} é 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 ⟨α∣a^∣α⟩=α\left\langle \alpha\left|\widehat{a}\right|\alpha\right\rangle =\alpha, então os valores quântios 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}

Vamos utilizar o seguinte operador identidadeShao (2016)John Wiley & Sons, Ltd (n.d.):

12πℏ∫−∞+∞∫−∞+∞dx0dp0∣α⟩⟨α∣=I\frac{1}{2\pi\hbar}\int^{+\infty}_{-\infty}\int^{+\infty}_{-\infty}dx_{0}dp_{0}\left|\alpha\right\rangle \left\langle \alpha\right|=\mathbb{I}

Inserindo esta identidade no propagador, temos que:

K(x′′,t;x′,t0)=12πℏ∫−∞+∞∫−∞+∞dx0dp0⟨x′′∣U^∣α⟩⟨α∣x′⟩K\left(x'',t;x',t_{0}\right)=\frac{1}{2\pi\hbar}\int^{+\infty}_{-\infty}\int^{+\infty}_{-\infty}dx_{0}dp_{0}\left\langle x''\right|\widehat{U}\left|\alpha\right\rangle \left\langle \alpha|x'\right\rangle

Agora queremos reescrever U^∣α⟩\widehat{U}\left|\alpha\right\rangle. Calculando então U^†a^U^\widehat{U}^{\dagger}\widehat{a}\widehat{U}. Então:

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 pelo lema de Baker-Hausdorff, conforme equaçãoShao (2016)Maziero (2021)Sakurai & Napolitano (2011):

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*}

Então 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} Wheeler (2017), 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)}

Ou então, 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*}

Então manipulando 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. Voltando então ao propagador, podemos reescrever como:

K(x′′,t;x′,t0)=12πℏ∫−∞+∞∫−∞+∞dx0dp0⟨x′′∣α(t)⟩⟨α∣x′⟩K\left(x'',t;x',t_{0}\right)=\frac{1}{2\pi\hbar}\int^{+\infty}_{-\infty}\int^{+\infty}_{-\infty}dx_{0}dp_{0}\left\langle x''|\alpha\left(t\right)\right\rangle \left\langle \alpha|x'\right\rangle

Agora queremos manipular o termo ⟨α∣x′⟩\left\langle \alpha|x'\right\rangle. Na configuração de espaço temos:

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

Onde usamos o resultado obtido na aula 8 ⟨x∣p∣a⟩=−iℏ∂∂x⟨x∣a⟩\left\langle x\left|p\right|a\right\rangle =-i\hbar\frac{\partial}{\partial x}\left\langle x|a\right\rangle Temos então uma equação diferencial de primeira ordem:

ℏ2mω∂∂x⟨x∣α⟩=−(mω2ℏx−α)⟨x∣α⟩\begin{align*} \sqrt{\frac{\hbar}{2m\omega}}\frac{\partial}{\partial x}\left\langle x|\alpha\right\rangle & =-\left(\sqrt{\frac{m\omega}{2\hbar}}x-\alpha\right)\left\langle x|\alpha\right\rangle \end{align*}

Em que a solução é:

⟨x∣α⟩=Cexp⁡[−(mω2ℏx−α)2]\begin{align*} \left\langle x|\alpha\right\rangle & =C\exp\left[-\left(\sqrt{\frac{m\omega}{2\hbar}}x-\alpha\right)^{2}\right] \end{align*}

Confore pode ser verificado facilmente substituindo na equação diferencial. Abrindo mais temos:

⟨x∣α⟩=Cexp⁡[−(mω2ℏx2+α2−2mωℏxα)]\left\langle x|\alpha\right\rangle =C\exp\left[-\left(\frac{m\omega}{2\hbar}x^{2}+\alpha^{2}-\sqrt{\frac{2m\omega}{\hbar}}x\alpha\right)\right]

Usando o obtido na equação 16:

α2=(mω2ℏx0+i2ℏmωp0)2=mω2ℏx02−12ℏmωp02+ix0p02ℏ\begin{align*} \alpha^{2} & =\left(\sqrt{\frac{m\omega}{2\hbar}}x_{0}+\frac{i}{\sqrt{2\hbar m\omega}}p_{0}\right)^{2}\\ & =\frac{m\omega}{2\hbar}x^{2}_{0}-\frac{1}{2\hbar m\omega}p^{2}_{0}+i\frac{x_{0}p_{0}}{2\hbar} \end{align*}

Então:

⟨x∣α⟩=Cexp⁡[−mω2ℏx2+2mωℏxα−(mω2ℏx02−12ℏmωp02−ix0p02ℏ)]=Cexp⁡[−mω2ℏ(x2+x02)+2mωℏxα+12ℏmωp02−ix0p02ℏ]=Cexp⁡[−mω2ℏ(x2+x02)+2mωℏx(mω2ℏx0+i2ℏmωp0)+12ℏmωp02−ix0p02ℏ]=Cexp⁡[−mω2ℏ(x2+x02)+2mω2ℏxx0+iℏxp0+12ℏmωp02−ix0p02ℏ]==Cexp⁡[−mω2ℏ(x2−2xx0+x02)+iℏxp0+12ℏmωp02−ix0p02ℏ]=Cexp⁡[−mω2ℏ(x−x0)2+p022ℏmω+i(xp0ℏ−x0p02ℏ)]\begin{align*} \left\langle x|\alpha\right\rangle & =C\exp\left[-\frac{m\omega}{2\hbar}x^{2}+\sqrt{\frac{2m\omega}{\hbar}}x\alpha-\left(\frac{m\omega}{2\hbar}x^{2}_{0}-\frac{1}{2\hbar m\omega}p^{2}_{0}-i\frac{x_{0}p_{0}}{2\hbar}\right)\right]\\ & =C\exp\left[-\frac{m\omega}{2\hbar}\left(x^{2}+x^{2}_{0}\right)+\sqrt{\frac{2m\omega}{\hbar}}x\alpha+\frac{1}{2\hbar m\omega}p^{2}_{0}-i\frac{x_{0}p_{0}}{2\hbar}\right]\\ & =C\exp\left[-\frac{m\omega}{2\hbar}\left(x^{2}+x^{2}_{0}\right)+\sqrt{\frac{2m\omega}{\hbar}}x\left(\sqrt{\frac{m\omega}{2\hbar}}x_{0}+\frac{i}{\sqrt{2\hbar m\omega}}p_{0}\right)+\frac{1}{2\hbar m\omega}p^{2}_{0}-i\frac{x_{0}p_{0}}{2\hbar}\right]\\ & =C\exp\left[-\frac{m\omega}{2\hbar}\left(x^{2}+x^{2}_{0}\right)+\frac{2m\omega}{2\hbar}xx_{0}+\frac{i}{\hbar}xp_{0}+\frac{1}{2\hbar m\omega}p^{2}_{0}-i\frac{x_{0}p_{0}}{2\hbar}\right]=\\ & =C\exp\left[-\frac{m\omega}{2\hbar}\left(x^{2}-2xx_{0}+x^{2}_{0}\right)+\frac{i}{\hbar}xp_{0}+\frac{1}{2\hbar m\omega}p^{2}_{0}-i\frac{x_{0}p_{0}}{2\hbar}\right]\\ & =C\exp\left[-\frac{m\omega}{2\hbar}\left(x-x_{0}\right)^{2}+\frac{p^{2}_{0}}{2\hbar m\omega}+i\left(\frac{xp_{0}}{\hbar}-\frac{x_{0}p_{0}}{2\hbar}\right)\right] \end{align*}

Para a condição de normalização da função de onda, então temos que:

1=∫−∞+∞∣⟨x∣α⟩∣2dx1=∣C∣2∫−∞+∞exp⁡[−mω2ℏ(x−x0)2+p02ℏmω]dxexp⁡[−p02ℏmω]=∣C∣2∫−∞+∞exp⁡[−mω2ℏ(x−x0)2]dx\begin{align*} 1 & =\int^{+\infty}_{-\infty}\left|\left\langle x|\alpha\right\rangle \right|^{2}dx\\ 1 & =\left|C\right|^{2}\int^{+\infty}_{-\infty}\exp\left[-\frac{m\omega}{2\hbar}\left(x-x_{0}\right)^{2}+\frac{p^{2}_{0}}{\hbar m\omega}\right]dx\\ \exp\left[-\frac{p^{2}_{0}}{\hbar m\omega}\right] & =\left|C\right|^{2}\int^{+\infty}_{-\infty}\exp\left[-\frac{m\omega}{2\hbar}\left(x-x_{0}\right)^{2}\right]dx \end{align*}

Fazendo u=mω2ℏx−αu=\sqrt{\frac{m\omega}{2\hbar}}x-\alpha, então dudx=mω2ℏ\frac{du}{dx}=\sqrt{\frac{m\omega}{2\hbar}}:

mω2ℏexp⁡[−p02ℏmω]=∣C∣2∫−∞+∞e−2u2dumω2ℏexp⁡[−p02ℏmω]=∣C∣2π2(mωπℏ)1/4exp⁡[−p022ℏmω]=C\begin{align*} \sqrt{\frac{m\omega}{2\hbar}}\exp\left[-\frac{p^{2}_{0}}{\hbar m\omega}\right] & =\left|C\right|^{2}\int^{+\infty}_{-\infty}e^{-2u^{2}}du\\ \sqrt{\frac{m\omega}{2\hbar}}\exp\left[-\frac{p^{2}_{0}}{\hbar m\omega}\right] & =\left|C\right|^{2}\sqrt{\frac{\pi}{2}}\\ \left(\frac{m\omega}{\pi\hbar}\right)^{1/4}\exp\left[-\frac{p^{2}_{0}}{2\hbar m\omega}\right] & =C \end{align*}

Onde usamos a integral de gauss ∫−∞+∞e−ax2dx=πa\int^{+\infty}_{-\infty}e^{-ax^{2}}dx=\sqrt{\frac{\pi}{a}}. Então:

⟨x∣α⟩=(mωπℏ)1/4exp⁡[−p022ℏmω]exp⁡[−mω2ℏ(x−x0)2+p022ℏmω+i(xp0ℏ−x0p02ℏ)]⟨x∣α⟩=(mωπℏ)1/4exp⁡[−mω2ℏ(x−x0)2+i(xp0ℏ−x0p02ℏ)]\begin{align*} \left\langle x|\alpha\right\rangle & =\left(\frac{m\omega}{\pi\hbar}\right)^{1/4}\exp\left[-\frac{p^{2}_{0}}{2\hbar m\omega}\right]\exp\left[-\frac{m\omega}{2\hbar}\left(x-x_{0}\right)^{2}+\frac{p^{2}_{0}}{2\hbar m\omega}+i\left(\frac{xp_{0}}{\hbar}-\frac{x_{0}p_{0}}{2\hbar}\right)\right]\\ \left\langle x|\alpha\right\rangle & =\left(\frac{m\omega}{\pi\hbar}\right)^{1/4}\exp\left[-\frac{m\omega}{2\hbar}\left(x-x_{0}\right)^{2}+i\left(\frac{xp_{0}}{\hbar}-\frac{x_{0}p_{0}}{2\hbar}\right)\right] \end{align*}

Para ⟨x∣α(t)⟩\left\langle x|\alpha\left(t\right)\right\rangle temos um resultado análogo, substituindo p0p_{0} e q0q_{0} respectivamente por p(t)=⟨α(t)∣p^∣α(t)⟩p\left(t\right)=\left\langle \alpha\left(t\right)\left|\widehat{p}\right|\alpha\left(t\right)\right\rangle e q(t)=⟨α(t)∣q^∣α(t)⟩q\left(t\right)=\left\langle \alpha\left(t\right)\left|\widehat{q}\right|\alpha\left(t\right)\right\rangle. Para encontrarmos estamos funções, fazendo então:

⟨α(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*}

Então:

α=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}. Então podemos escrever:

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{2m\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{2m\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{2m\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{2m\omega}} \end{align}

Que é o resultado obtido na equação 16. Substituindo x(τ)x\left(\tau\right) e p(τ)p\left(\tau\right) em ⟨x′′∣α(t)⟩⟨α∣x′⟩=⟨x′′∣α(t)⟩⟨x′∣α⟩†\left\langle x''|\alpha\left(t\right)\right\rangle \left\langle \alpha|x'\right\rangle =\left\langle x''|\alpha\left(t\right)\right\rangle \left\langle x'|\alpha\right\rangle ^{\dagger} , Com o auxílio do matemática, e escrevendo x′=x1x'=x_{1} e x′′=x2x''=x_{2}:

Subscript[Q,11]:=((m*w)/(Pi*h))^(1/4)*Exp[-(m*w)/(2*h)*(Subscript[x,1]-Subscript[x,0])^2
+I*Subscript[x,1]*Subscript[p,0]/h-I*Subscript[x,0]*Subscript[p,0]/(2*h)]
Subscript[Q,12]:=((m*w)/(Pi*h))^(1/4)*Exp[-(m*w)/(2*h)*(Subscript[x,1]-Subscript[x,0])^2
-I*Subscript[x,1]*Subscript[p,0]/h+I*Subscript[x,0]*Subscript[p,0]/(2*h)]
Subscript[Q,21]:=((m*w)/(Pi*h))^(1/4)*Exp[-(m*w)/(2*h)*(Subscript[x,2]-Subscript[x,t])^2
+I*Subscript[x,2]*Subscript[p,t]/h-I*Subscript[x,t]*Subscript[p,t]/(2*h)]
Subscript[Q,22]:=Subscript[Q,21]/.{Subscript[x,t]->Subscript[x,0]*Cos[w*t]+
Subscript[p,0]*Sin[w*t]/(m*w),
Subscript[p,t]->-m*w*Subscript[x,0]*Sin[w*t]+Subscript[p,0]*Cos[w*t]}
Subscript[Q,12]*Subscript[Q,22]/(2*Pi*h)

Obtemos:

K(x′′,t;x′,t0)=mωℏ2π3/2ℏexp⁡(ix′′(p0cos⁡(ωτ)−mωx0sin⁡(ωτ))ℏ−mω(−p0sin⁡(ωτ)mω−x0cos⁡(ωτ)+x′′)22h)×exp⁡(−i(p0sin⁡(ωτ)mω+x0cos⁡(τω))(p0cos⁡(ωτ)−mωx0sin⁡(ωτ))2ℏ−mω(x′−x0)22ℏ+ip0x02ℏ−ip0x′ℏ)\begin{align*} K\left(x'',t;x',t_{0}\right) & =\frac{\sqrt{\frac{m\omega}{\hbar}}}{2\pi^{3/2}\hbar}\exp\left(\frac{ix''\left(p_{0}\cos(\omega\tau)-m\omega x_{0}\sin(\omega\tau)\right)}{\hbar}-\frac{m\omega\left(-\frac{p_{0}\sin(\omega\tau)}{m\omega}-x_{0}\cos(\omega\tau)+x''\right){}^{2}}{2h}\right)\\ & \times\exp\left(-\frac{i\left(\frac{p_{0}\sin(\omega\tau)}{m\omega}+x_{0}\cos(\tau\omega)\right)\left(p_{0}\cos(\omega\tau)-m\omega x_{0}\sin(\omega\tau)\right)}{2\hbar}-\frac{m\omega\left(x'-x_{0}\right){}^{2}}{2\hbar}+\frac{ip_{0}x_{0}}{2\hbar}-\frac{ip_{0}x'}{\hbar}\right) \end{align*}

Utilizando novamente o mathematica para simplificar o expoente com a função ’TrigToExp’, obtemos que:

β=p02e−2iωτ4ℏmω−p024ℏmω+mωx0x2e−iωτℏ−mωx02e−2iωτ4ℏ+mωx0x′ℏ−mωx′22ℏ−mωx′′22ℏ−3mωx024ℏ+ip0x′′e−iωτℏ−ip0x0e−2iωτ2ℏ+ip0x02ℏ−ip0x′ℏ\begin{align*} \beta & =\frac{p^{2}_{0}e^{-2i\omega\tau}}{4\hbar m\omega}-\frac{p^{2}_{0}}{4\hbar m\omega}+\frac{m\omega x_{0}x_{2}e^{-i\omega\tau}}{\hbar}-\frac{m\omega x^{2}_{0}e^{-2i\omega\tau}}{4\hbar}+\frac{m\omega x_{0}x'}{\hbar}-\frac{m\omega x'^{2}}{2\hbar}\\ & -\frac{m\omega x''^{2}}{2\hbar}-\frac{3m\omega x^{2}_{0}}{4\hbar}+\frac{ip_{0}x''e^{-i\omega\tau}}{\hbar}-\frac{ip_{0}x_{0}e^{-2i\omega\tau}}{2\hbar}+\frac{ip_{0}x_{0}}{2\hbar}-\frac{ip_{0}x'}{\hbar} \end{align*}

Abrindo os termos e multiplicando o termo que aompanha p02p^{2}_{0} por ii\frac{i}{i}:

β=−mω2ℏ(x′2+x′′2)−mω2ℏ(e−2iωτ2+12+1)x02−i2ℏmω(1−e−2iωτ2i)p02+mωℏ(x′′e−iωτ+x′)x0+iℏ(x′′e−iωτ−x′)p0−iℏ(e−2iωτ2−12)p0x0\begin{align*} \beta & =-\frac{m\omega}{2\hbar}\left(x'^{2}+x''^{2}\right)-\frac{m\omega}{2\hbar}\left(\frac{e^{-2i\omega\tau}}{2}+\frac{1}{2}+1\right)x^{2}_{0}-\frac{i}{2\hbar m\omega}\left(\frac{1-e^{-2i\omega\tau}}{2i}\right)p^{2}_{0}\\ & +\frac{m\omega}{\hbar}\left(x''e^{-i\omega\tau}+x'\right)x_{0}+\frac{i}{\hbar}\left(x''e^{-i\omega\tau}-x'\right)p_{0}-\frac{i}{\hbar}\left(\frac{e^{-2i\omega\tau}}{2}-\frac{1}{2}\right)p_{0}x_{0} \end{align*}

Agora como

cos(ωt)e−iωt=(eiωt+e−iωt2)e−iωt=(12+e−2iωt)=(12+e−2iωt)cos\left(\omega t\right)e^{-i\omega t}=\left(\frac{e^{i\omega t}+e^{-i\omega t}}{2}\right)e^{-i\omega t}=\left(\frac{1}{2}+e^{-2i\omega t}\right)=\left(\frac{1}{2}+e^{-2i\omega t}\right)
sin⁡(ωt)e−ωit=(eiωt−e−iωt2i)e−ωit=(1−e−2ωit2i)e−ωit=i(e−2ωit2−12)e−iωt\sin\left(\omega t\right)e^{-\omega it}=\left(\frac{e^{i\omega t}-e^{-i\omega t}}{2i}\right)e^{-\omega it}=\left(\frac{1-e^{-2\omega it}}{2i}\right)e^{-\omega it}=i\left(\frac{e^{-2\omega it}}{2}-\frac{1}{2}\right)e^{-i\omega t}

Ajustando os sinais e manipulando:

β=−mω2ℏ(x′2+x′′2)−[mω2ℏ(1+cos(ωτ)e−iωτ)]x02−[i2ℏmω(sin⁡(ωτ)e−iωτ)]p02−[1ℏsin⁡(ωτ)e−iωτ]p0x0−[−mωℏ(x′′e−iωτ+x′)]x0−[−iℏ(x′′e−iωτ−x′)]p0\begin{align*} \beta & =-\frac{m\omega}{2\hbar}\left(x'^{2}+x''^{2}\right)-\left[\frac{m\omega}{2\hbar}\left(1+cos\left(\omega\tau\right)e^{-i\omega\tau}\right)\right]x^{2}_{0}-\left[\frac{i}{2\hbar m\omega}\left(\sin\left(\omega\tau\right)e^{-i\omega\tau}\right)\right]p^{2}_{0}\\ & -\left[\frac{1}{\hbar}\sin\left(\omega\tau\right)e^{-i\omega\tau}\right]p_{0}x_{0}-\left[-\frac{m\omega}{\hbar}\left(x''e^{-i\omega\tau}+x'\right)\right]x_{0}-\left[-\frac{i}{\hbar}\left(x''e^{-i\omega\tau}-x'\right)\right]p_{0} \end{align*}

Aqui precisamos notar que o term o c22c_{22} possui um fator multiplicativo 12\frac{1}{2} em comparação com a equação (37) do artigo, possívelmente um erro de digitação. Além disso na equação (34) foi negligenciado o termo que acompanha a exponencial, talvez intencionlmente. E rearranjando:

β=−mωπℏ(x12+x22)−c11x02−c22p02−c1x0−c12p0x0−c1p0\begin{align*} \beta & =-\frac{m\omega}{\pi\hbar}\left(x^{2}_{1}+x^{2}_{2}\right)-c_{11}x^{2}_{0}-c_{22}p^{2}_{0}-c_{1}x_{0}-c_{12}p_{0}x_{0}-c_{1}p_{0} \end{align*}

Substituindo no propagador e integrando primeiro em p0p_{0}temos então:

K(x′′,t;x′,t0)=mωπℏexp⁡[−mω2ℏ[(x2+x′2)]]2πℏ∫−∞+∞exp⁡[−c11x02−c1x0]∫−∞+∞exp⁡[−c22p02+(−c12x0−c2)p0]dp0dx0K\left(x'',t;x',t_{0}\right)=\frac{\sqrt{\frac{m\omega}{\pi\hbar}}\exp\left[-\frac{m\omega}{2\hbar}\left[\left(x^{2}+x'^{2}\right)\right]\right]}{2\pi\hbar}\int^{+\infty}_{-\infty}\exp\left[-c_{11}x^{2}_{0}-c_{1}x_{0}\right]\int^{+\infty}_{-\infty}\exp\left[-c_{22}p^{2}_{0}+\left(-c_{12}x_{0}-c_{2}\right)p_{0}\right]dp_{0}dx_{0}

Evidentemente poderíamos completar quadrado para realizar uma integração de gauss, mas usando uma forma alternativa generalizada como ∫−∞+∞e−ax2+bx+c=πaeb24a+c\int^{+\infty}_{-\infty}e^{-ax^{2}+bx+c}=\sqrt{\frac{\pi}{a}}e^{\frac{b^{2}}{4a}+c}, tendo a=c22a=c_{22}, b=(c12x0+c2)b=\left(c_{12}x_{0}+c_{2}\right) e c=0c=0:

K(x′′,t;x′,t0)=mωπℏexp⁡[−mω2ℏ[(x2+x′2)]]2πℏπc22∫−∞+∞exp⁡[−c11x02−c1x0+(c12x0+c2)24c22]dx0=mωℏc2212πℏexp⁡[−mω2ℏ[(x2+x′2)]]∫−∞+∞exp⁡[−c11x02−c1x0+c122x02+c22+2c12c2x04c22]dx0=mωℏc2212πℏexp⁡[−mω2ℏ[(x2+x′2)]]∫−∞+∞exp⁡[−(c11−c1224c22)x02+(c12c22c22−c1)x0+(c224c22)]dx0=mωℏc2212πℏπ(c11−14c122c22)exp⁡[−mω2ℏ[(x2+x′2)]+14(c12c22c22−c1)2(c11−c1224c22)+14c22c22]=mωℏc2212πℏ114c22π(4c22c11−c122)exp⁡[−mω2ℏ[(x2+x′2)]+14(c12c2−2c22c1)214c222(4c22c11−c122)14c22+14c22c22]=mωℏ1πℏπ(4c22c11−c122)exp⁡[−mω2ℏ[(x2+x′2)]+14(c12c2−2c22c1)2(4c22c11−c122)c22+14c22c22]=mωπℏ1ℏ14c22c11−c122exp⁡[−mω2ℏ[(x2+x′2)]+14c122c22+4c222c12−4c12c2c22c1+4c22c22c11−c22c122(4c22c11−c122)c22]=mωπℏ1ℏ14c22c11−c122exp⁡[−mω2ℏ[(x2+x′2)]+c22c11−c12c2c1+c222c124c22c11−c122]\begin{align*} K\left(x'',t;x',t_{0}\right) & =\sqrt{\frac{m\omega}{\pi\hbar}}\frac{\exp\left[-\frac{m\omega}{2\hbar}\left[\left(x^{2}+x'^{2}\right)\right]\right]}{2\pi\hbar}\sqrt{\frac{\pi}{c_{22}}}\int^{+\infty}_{-\infty}\exp\left[-c_{11}x^{2}_{0}-c_{1}x_{0}+\frac{\left(c_{12}x_{0}+c_{2}\right)^{2}}{4c_{22}}\right]dx_{0}\\ & =\sqrt{\frac{m\omega}{\hbar c_{22}}}\frac{1}{2\pi\hbar}\exp\left[-\frac{m\omega}{2\hbar}\left[\left(x^{2}+x'^{2}\right)\right]\right]\int^{+\infty}_{-\infty}\exp\left[-c_{11}x^{2}_{0}-c_{1}x_{0}+\frac{c^{2}_{12}x^{2}_{0}+c^{2}_{2}+2c_{12}c_{2}x_{0}}{4c_{22}}\right]dx_{0}\\ & =\sqrt{\frac{m\omega}{\hbar c_{22}}}\frac{1}{2\pi\hbar}\exp\left[-\frac{m\omega}{2\hbar}\left[\left(x^{2}+x'^{2}\right)\right]\right]\int^{+\infty}_{-\infty}\exp\left[-\left(c_{11}-\frac{c^{2}_{12}}{4c_{22}}\right)x^{2}_{0}+\left(\frac{c_{12}c_{2}}{2c_{22}}-c_{1}\right)x_{0}+\left(\frac{c^{2}_{2}}{4c_{22}}\right)\right]dx_{0}\\ & =\sqrt{\frac{m\omega}{\hbar c_{22}}}\frac{1}{2\pi\hbar}\sqrt{\frac{\pi}{\left(c_{11}-\frac{1}{4}\frac{c^{2}_{12}}{c_{22}}\right)}}\exp\left[-\frac{m\omega}{2\hbar}\left[\left(x^{2}+x'^{2}\right)\right]+\frac{1}{4}\frac{\left(\frac{c_{12}c_{2}}{2c_{22}}-c_{1}\right)^{2}}{\left(c_{11}-\frac{c^{2}_{12}}{4c_{22}}\right)}+\frac{1}{4}\frac{c^{2}_{2}}{c_{22}}\right]\\ & =\sqrt{\frac{m\omega}{\hbar c_{22}}}\frac{1}{2\pi\hbar}\frac{1}{\sqrt{\frac{1}{4c_{22}}}}\sqrt{\frac{\pi}{\left(4c_{22}c_{11}-c^{2}_{12}\right)}}\exp\left[-\frac{m\omega}{2\hbar}\left[\left(x^{2}+x'^{2}\right)\right]+\frac{1}{4}\frac{\left(c_{12}c_{2}-2c_{22}c_{1}\right)^{2}\frac{1}{4c^{2}_{22}}}{\left(4c_{22}c_{11}-c^{2}_{12}\right)\frac{1}{4c_{22}}}+\frac{1}{4}\frac{c^{2}_{2}}{c_{22}}\right]\\ & =\sqrt{\frac{m\omega}{\hbar}}\frac{1}{\pi\hbar}\sqrt{\frac{\pi}{\left(4c_{22}c_{11}-c^{2}_{12}\right)}}\exp\left[-\frac{m\omega}{2\hbar}\left[\left(x^{2}+x'^{2}\right)\right]+\frac{1}{4}\frac{\left(c_{12}c_{2}-2c_{22}c_{1}\right)^{2}}{\left(4c_{22}c_{11}-c^{2}_{12}\right)c_{22}}+\frac{1}{4}\frac{c^{2}_{2}}{c_{22}}\right]\\ & =\sqrt{\frac{m\omega}{\pi\hbar}}\frac{1}{\hbar}\frac{1}{\sqrt{4c_{22}c_{11}-c^{2}_{12}}}\exp\left[-\frac{m\omega}{2\hbar}\left[\left(x^{2}+x'^{2}\right)\right]+\frac{1}{4}\frac{c^{2}_{12}c^{2}_{2}+4c^{2}_{22}c^{2}_{1}-4c_{12}c_{2}c_{22}c_{1}+4c^{2}_{2}c_{22}c_{11}-c^{2}_{2}c^{2}_{12}}{\left(4c_{22}c_{11}-c^{2}_{12}\right)c_{22}}\right]\\ & =\sqrt{\frac{m\omega}{\pi\hbar}}\frac{1}{\hbar}\frac{1}{\sqrt{4c_{22}c_{11}-c^{2}_{12}}}\exp\left[-\frac{m\omega}{2\hbar}\left[\left(x^{2}+x'^{2}\right)\right]+\frac{c^{2}_{2}c_{11}-c_{12}c_{2}c_{1}+c^{2}_{22}c^{2}_{1}}{4c_{22}c_{11}-c^{2}_{12}}\right] \end{align*}

Manipulando então:

4c22c11−c122=4mω2ℏ(1+cos⁡(ωτ)e−iωτ)isin⁡(ωτ)e−iωτ2ℏmω−(1ℏsin⁡(ωτ)e−iωτ)2=1ℏ2[isin⁡(ωτ)e−iωτicos⁡(ωτ)sin⁡(ωτ)e−i2ωτ−sin⁡2(ωτ)e−i2ωτ]=1ℏ2[isin⁡(ωτ)e−iωτ+(icos⁡(ωτ)−sin⁡(ωτ))sin⁡(ωτ)e−i2ωτ]=1ℏ2[sin⁡(ωτ)e−iωτ+(cos⁡(ωτ)+isin⁡(ωτ))sin⁡(ωτ)e−i2ωτ]=1ℏ2[sin⁡(ωτ)e−iωτ+e−iωtsin⁡(ωτ)e−i2ωτ]=−isin⁡(ωτ)h2[e−iωτ+e−iωτ]=2isin⁡(ωt)ℏ2e−iωτ\begin{align} \sqrt{4c_{22}c_{11}-c^{2}_{12}} & =4\frac{m\omega}{2\hbar}\left(1+\cos\left(\omega\tau\right)e^{-i\omega\tau}\right)i\frac{\sin\left(\omega\tau\right)e^{-i\omega\tau}}{2\hbar m\omega}-\left(\frac{1}{\hbar}\sin\left(\omega\tau\right)e^{-i\omega\tau}\right)^{2}\\ & =\frac{1}{\hbar^{2}}\left[i\sin\left(\omega\tau\right)e^{-i\omega\tau}i\cos\left(\omega\tau\right)\sin\left(\omega\tau\right)e^{-i2\omega\tau}-\sin^{2}\left(\omega\tau\right)e^{-i2\omega\tau}\right]\\ & =\frac{1}{\hbar^{2}}\left[i\sin\left(\omega\tau\right)e^{-i\omega\tau}+\left(i\cos\left(\omega\tau\right)-\sin\left(\omega\tau\right)\right)\sin\left(\omega\tau\right)e^{-i2\omega\tau}\right]\\ & =\frac{1}{\hbar^{2}}\left[\sin\left(\omega\tau\right)e^{-i\omega\tau}+\left(\cos\left(\omega\tau\right)+i\sin\left(\omega\tau\right)\right)\sin\left(\omega\tau\right)e^{-i2\omega\tau}\right]\\ & =\frac{1}{\hbar^{2}}\left[\sin\left(\omega\tau\right)e^{-i\omega\tau}+e^{-i\omega t}\sin\left(\omega\tau\right)e^{-i2\omega\tau}\right]\\ & =-\frac{i\sin\left(\omega\tau\right)}{h^{2}}\left[e^{-i\omega\tau}+e^{-i\omega\tau}\right]\\ & =2\frac{i\sin\left(\omega t\right)}{\hbar^{2}}e^{-i\omega\tau} \end{align}

Substituindo então

mωπℏ1ℏ14c22c11−c122=mωπℏ1ℏ12isin⁡(ωτ)h2e−iωτ=mω2πℏisin⁡(ωτ)e−ωτ\begin{align*} \sqrt{\frac{m\omega}{\pi\hbar}}\frac{1}{\hbar}\frac{1}{\sqrt{4c_{22}c_{11}-c^{2}_{12}}} & =\sqrt{\frac{m\omega}{\pi\hbar}}\frac{1}{\hbar}\frac{1}{\sqrt{2\frac{i\sin\left(\omega\tau\right)}{h^{2}}e^{-i\omega\tau}}}\\ & =\sqrt{\frac{m\omega}{2\pi\hbar i\sin\left(\omega\tau\right)e^{-\omega\tau}}} \end{align*}

Que é o resultado que estávamos buscando, o termo eωτ/2e^{\omega\tau/2} é uma fase global que pode ser ignorandoShao (2016). Podemos obter o expoente através do Mathematica:

(\!\(\*SubsuperscriptBox[\(c\),\(2\),\(2\)]\ \*SubscriptBox[\(c\),\(11\)]\)- 
Subscript[c,1]Subscript[c,2]Subscript[c,12]+\!\( \*SubsuperscriptBox[\(c\),\(1\),
\(2\)]\ \*SubscriptBox[\(c\),\(22\)]\))/(\!\( \*SubsuperscriptBox[\(c\),\(12\),
\(2\)]-\(4\ \*SubscriptBox[\(c\),\(11\)]\ \*SubscriptBox[\(c\),\(22\)]\)\))/.{Subscript[c,11]
-> m*w*(1+Cos[w*t]*Exp[-I*w*t])/(2*h), Subscript[c,22]->I*Sin[w*t]*Exp[-I*w*t]/(2*h*m*w), 
Subscript[c,12]->Sin[w*t]*Exp[-I*w*t]/(h), Subscript[c,1]->-m*w*(Subscript[x,2]*Exp[-I*w*t]
+Subscript[x,1])/(h), Subscript[c,2]->-I*(Subscript[x,2]*Exp[-I*w*t]-Subscript[x,1])/(h)}

E então usar o comando ’Simplify’ , onde obtemos:

β=−mω2ℏ[(x22+x12)]+mω(−2x′′x′eiωτ+x′2+x′′2)ℏ(−1+e2itw)=−mω2ℏ[(x22+x12)]+mω(x′2+x′′2−eiωτ2x′′x′)eiωτℏ(eiωτ−e−itw)=−mω2ℏ[(x22+x12)]+imω(x′2+x′′2−eiωτ2x′′x′)ℏeiωτ2sin⁡(ωτ)=imω2ℏsin⁡(ωτ)[isin⁡(ωτ)[(x′′2+x′2)]+e−iωτ(x′2+x′′2−eiωτ2x′′x′)]=imω2ℏsin⁡(ωτ)[isin⁡(ωτ)[(x′′2+x′2)]+e−iωτ(x′′2+x′2)−2x′′x′]=imω2ℏsin⁡(ωτ)[(eiωτ−e−iωτ)2[(x′′2+x′2)]+e−iωτ(x′′2+x′2)−2x′′x′]=imω2ℏsin⁡(ωτ)[(eiωτ+e−iωτ)2[(x′′2+x′2)]−2x′′x′]=imω2ℏsin⁡(ωτ)[cos⁡(ωτ)[(x′′2+x′2)]−2x′′x′]\begin{align*} \beta & =-\frac{m\omega}{2\hbar}\left[\left(x^{2}_{2}+x^{2}_{1}\right)\right]+\frac{m\omega\left(-2x''x'e^{i\omega\tau}+x'^{2}+x''^{2}\right)}{\hbar\left(-1+e^{2itw}\right)}\\ & =-\frac{m\omega}{2\hbar}\left[\left(x^{2}_{2}+x^{2}_{1}\right)\right]+\frac{m\omega\left(x'^{2}+x''^{2}-e^{i\omega\tau}2x''x'\right)}{e^{i\omega\tau}\hbar\left(e^{i\omega\tau}-e^{-itw}\right)}\\ & =-\frac{m\omega}{2\hbar}\left[\left(x^{2}_{2}+x^{2}_{1}\right)\right]+\frac{im\omega\left(x'^{2}+x''^{2}-e^{i\omega\tau}2x''x'\right)}{\hbar e^{i\omega\tau}2\sin\left(\omega\tau\right)}\\ & =\frac{im\omega}{2\hbar\sin\left(\omega\tau\right)}\left[i\sin\left(\omega\tau\right)\left[\left(x''^{2}+x'^{2}\right)\right]+e^{-i\omega\tau}\left(x'^{2}+x''^{2}-e^{i\omega\tau}2x''x'\right)\right]\\ & =\frac{im\omega}{2\hbar\sin\left(\omega\tau\right)}\left[i\sin\left(\omega\tau\right)\left[\left(x''^{2}+x'^{2}\right)\right]+e^{-i\omega\tau}\left(x''^{2}+x'^{2}\right)-2x''x'\right]\\ & =\frac{im\omega}{2\hbar\sin\left(\omega\tau\right)}\left[\frac{\left(e^{i\omega\tau}-e^{-i\omega\tau}\right)}{2}\left[\left(x''^{2}+x'^{2}\right)\right]+e^{-i\omega\tau}\left(x''^{2}+x'^{2}\right)-2x''x'\right]\\ & =\frac{im\omega}{2\hbar\sin\left(\omega\tau\right)}\left[\frac{\left(e^{i\omega\tau}+e^{-i\omega\tau}\right)}{2}\left[\left(x''^{2}+x'^{2}\right)\right]-2x''x'\right]\\ & =\frac{im\omega}{2\hbar\sin\left(\omega\tau\right)}\left[\cos\left(\omega\tau\right)\left[\left(x''^{2}+x'^{2}\right)\right]-2x''x'\right] \end{align*}

Lembrando que τ=t−t0\tau=t-t_{0}, então o propagador é, como queríamos demonstrar:

K(x′′,t;x′,t0)=mω2πiℏsin⁡[ωt(t−t0)]exp⁡[imω2ℏsin⁡[ω(t−t0)]((x′′2+x′2)cos⁡[ω(t−t0)]−2x′′x′)]K\left(x'',t;x',t_{0}\right)=\sqrt{\frac{m\omega}{2\pi i\hbar\sin\left[\omega t\left(t-t_{0}\right)\right]}}\exp\left[\frac{im\omega}{2\hbar\sin\left[\omega\left(t-t_{0}\right)\right]}\left(\left(x''^{2}+x'^{2}\right)\cos\left[\omega\left(t-t_{0}\right)\right]-2x''x'\right)\right]

Problem 3

The propagator can also be written as

K(x′′,t;x′,t0)=⟨x′′,t∣x′,t9⟩K\left(x'',t;x',t_{0}\right)=\left\langle x'',t|x',t_{9}\right\rangle

where ∣x′,t0⟩\left|x',t_{0}\right\rangle and ⟨x′′,t∣\left\langle x'',t\right| are to be understood as an eigenket and an eigenbra of the position operator in the Heisenberg picture. The propagator in momentum space analogous to (1)(1) is given by ⟨p′′,t∣p′,t0⟩\left\langle p'',t|p',t_{0}\right\rangle. Derive an explicit expression for ⟨p′′,t∣p′,t0⟩\left\langle p'',t|p',t_{0}\right\rangle for the free-particle case.

Análogo ao que foi feito em aula, podemos fazer então:

⟨p′′∣α,t0;t⟩=⟨p′′∣e−H(t−t0)/ℏ∣α,t0⟩=∑a′⟨p′′∣a′⟩⟨a′∣α,t0⟩e−Ea′(t−t0)/ℏ=∑a′⟨p′′∣a′⟩⟨a′∣α,t0⟩e−Ea′(t−t0)/ℏ=∫d3p′∑a′⟨p′′∣a′⟩⟨a′∣p′⟩⟨p′∣α,t0⟩e−Ea′(t−t0)/ℏ=∫d3p′[∑a′⟨p′′∣a′⟩⟨a′∣p′⟩e−Ea′(t−t0)/ℏ]⟨p′∣α,t0⟩=∫d3p′K(p′′,t;p′,t0)⟨p′∣α,t0⟩\begin{align*} \left\langle p''|\alpha,t_{0};t\right\rangle & =\left\langle p''\left|e^{-H\left(t-t_{0}\right)/\hbar}\right|\alpha,t_{0}\right\rangle \\ & =\sum_{a'}\left\langle p''|a'\right\rangle \left\langle a'|\alpha,t_{0}\right\rangle e^{-E_{a'}\left(t-t_{0}\right)/\hbar}\\ & =\sum_{a'}\left\langle p''|a'\right\rangle \left\langle a'|\alpha,t_{0}\right\rangle e^{-E_{a'}\left(t-t_{0}\right)/\hbar}\\ & =\int d^{3}p'\sum_{a'}\left\langle p''|a'\right\rangle \left\langle a'|p'\right\rangle \left\langle p'|\alpha,t_{0}\right\rangle e^{-E_{a'}\left(t-t_{0}\right)/\hbar}\\ & =\int d^{3}p'\left[\sum_{a'}\left\langle p''|a'\right\rangle \left\langle a'|p'\right\rangle e^{-E_{a'}\left(t-t_{0}\right)/\hbar}\right]\left\langle p'|\alpha,t_{0}\right\rangle \\ & =\int d^{3}p'K\left(p'',t;p',t_{0}\right)\left\langle p'|\alpha,t_{0}\right\rangle \end{align*}

Temos então o propagador:

K(p′′,t;p′,t0)=∑a′⟨p′′∣a′⟩⟨a′∣p′⟩e−Ea′(t−t0)/ℏK\left(p'',t;p',t_{0}\right)=\sum_{a'}\left\langle p''|a'\right\rangle \left\langle a'|p'\right\rangle e^{-E_{a'}\left(t-t_{0}\right)/\hbar}

Para derivar uma expressão explícita para o caso da partícula livre, podemos escrever então mais elegantemente como:

K(p′′,t;p′,t0)=⟨p′′∣exp⁡[−iH(t−t0)ℏ]∣p′⟩K\left(p'',t;p',t_{0}\right)=\left\langle p''\left|\exp\left[\frac{-iH\left(t-t_{0}\right)}{\hbar}\right]\right|p'\right\rangle

Fazendo análogo à primeira questão, onde o hamiltoniano de uma partícula livre é dado por H=p22mH=\frac{p^{2}}{2m} Uma vez que ∣p′⟩\left|p'\right\rangle é um autoket do operador pp e HH, logo p∣p′⟩=p′∣p′⟩p\left|p'\right\rangle =p'\left|p'\right\rangle e H∣p′⟩=(p′22m)∣p′⟩H\left|p'\right\rangle =\left(\frac{p'^{2}}{2m}\right)\left|p'\right\rangle. Logo

K(p′′,t;p′,t0)=⟨p′′∣exp⁡[−iH(t−t0)ℏ]∣p′⟩=∫dp′′′⟨p′′∣p′′′⟩⟨p′′′∣exp⁡[−iH(t−t0)ℏ]∣p′⟩=∫dp′′′δ(p′′−p′′′)⟨p′′′′∣exp⁡[−iH(t−t0)ℏ]∣p′⟩=∫⟨p′′′∣p′⟩δ(p′′−p′′′)dp′′′exp⁡(−ip′2(t−t0)ℏ2m)=∫δ(p′′′−p′′)δ(p′′′−p′)dp′′′exp⁡(−ip′2(t−t0)ℏ2m)\begin{align*} K\left(p'',t;p',t_{0}\right) & =\left\langle p''\left|\exp\left[\frac{-iH\left(t-t_{0}\right)}{\hbar}\right]\right|p'\right\rangle \\ & =\int dp'''\left\langle p''|p'''\right\rangle \left\langle p'''\left|\exp\left[\frac{-iH\left(t-t_{0}\right)}{\hbar}\right]\right|p'\right\rangle \\ & =\int dp'''\delta\left(p''-p'''\right)\left\langle p''''\left|\exp\left[\frac{-iH\left(t-t_{0}\right)}{\hbar}\right]\right|p'\right\rangle \\ & =\int\left\langle p'''|p'\right\rangle \delta\left(p''-p'''\right)dp'''\exp\left(\frac{-ip'^{2}\left(t-t_{0}\right)}{\hbar2m}\right)\\ & =\int\delta\left(p'''-p''\right)\delta\left(p'''-p'\right)dp'''\exp\left(\frac{-ip'^{2}\left(t-t_{0}\right)}{\hbar2m}\right) \end{align*}

Agora sendo f(p′′′)=δ(p′′′−p′′)f\left(p'''\right)=\delta\left(p'''-p''\right) e ∫f(p′′′)δ(p′′′−p′)dp′′′=f(p′)=δ(p′−p′′)\int f\left(p'''\right)\delta\left(p'''-p'\right)dp'''=f\left(p'\right)=\delta\left(p'-p''\right). Ficamos então com:

K(p′′,t;p′,t0)=δ(p′−p′′)exp⁡(−ip′2(t−t0)2mℏ)K\left(p'',t;p',t_{0}\right)=\delta\left(p'-p''\right)\exp\left(\frac{-ip'^{2}\left(t-t_{0}\right)}{2m\hbar}\right)
References
  1. Shao, J. (2016). Elementary derivation of the quantum propagator for the harmonic oscillator. American Journal of Physics, 84(10), 770–774. 10.1119/1.4960479
  2. Coherent States in Quantum Physics. (n.d.). John Wiley & Sons, Ltd.
  3. Maziero, J. (2021). Funções matriciais de matrizes. https://nbviewer.org/github/jonasmaziero/algebra_linear/blob/ead5a0acc24539eb1b7510fe656749a39a05b2be/04_functions.ipynb
  4. Sakurai, J. J., & Napolitano, J. (2011). Modern Quantum Mechanics. Addison-Wesley. https://books.google.com.br/books?id=N4I-AQAACAAJ
  5. Wheeler, G. (2017). The number operator. http://www.physics.usu.edu/Wheeler/QuantumMechanics/QM16SHOQuestions.pdf