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 6

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

Problem 1

Note that if one associates with the coordinates xμx^{\mu} of a space-time point a matrix

X=(x0+x3x1−ix2x1+ix2x0−x3)=xμσμX=\left(\begin{array}{cc} x^{0}+x^{3} & x^{1}-ix^{2}\\ x^{1}+ix^{2} & x^{0}-x^{3} \end{array}\right)=x^{\mu}\sigma_{\mu}

then det⁡X=D(X)=gμνxμxν\det X=D\left(X\right)=g_{\mu\nu}x^{\mu}x^{\nu} has the value D(X)=gμνxμxνD\left(X\right)=g_{\mu\nu}x^{\mu}x^{\nu} .

a) If the matrix

W=(αβγδ)W=\left(\begin{array}{cc} \alpha & \beta\\ \gamma & \delta \end{array}\right)

with complex elements α\alpha, β\beta, γ\gamma, δ\delta has determinant 1 and

W∗=(α‾γ‾β‾δ‾)W^{*}=\left(\begin{array}{cc} \overline{\alpha} & \overline{\gamma}\\ \overline{\beta} & \overline{\delta} \end{array}\right)

is its Hermitian adjoint, show that

X′=WXW∗X'=WXW^{*}

defines a proper Lorentz transformation x′μ=Λνμxνx^{'\mu}=\Lambda^{\mu}_{\nu}x^{\nu} , (i.e., show that det⁡Λ=+1,Λ00≥1\det\Lambda=+1,\Lambda^{0}_{0}\ge1, and that D(X)=D(X′)D\left(X\right)=D\left(X'\right))

Temos:

X=(x0+x3x1−ix2x1+ix2x0−x3)=xμσμX=\left(\begin{array}{cc} x^{0}+x^{3} & x^{1}-ix^{2}\\ x^{1}+ix^{2} & x^{0}-x^{3} \end{array}\right)=x^{\mu}\sigma_{\mu}

e det⁡X=D(X)=gμνxμxν\det X=D\left(X\right)=g_{\mu\nu}x^{\mu}x^{\nu}. Então se temos uma matrix com determinante 1:

W=(αβγδ)W=\left(\begin{array}{cc} \alpha & \beta\\ \gamma & \delta \end{array}\right)

E sua adjunta hermitiana é:

W∗=(α‾γ‾β‾δ‾)W^{*}=\left(\begin{array}{cc} \overline{\alpha} & \overline{\gamma}\\ \overline{\beta} & \overline{\delta} \end{array}\right)

Então queremos mostrar que:

X′=WXW∗X'=WXW^{*}

determine uma transormação de Lorentz x′μ=Λνμxνx^{'\mu}=\Lambda^{\mu}_{\nu}x^{\nu}. Isto é, mostrar que det⁡Λ=+1\det\Lambda=+1, Λ00≥1\Lambda^{0}_{0}\ge1 e D(X)=D(X′)D\left(X\right)=D\left(X'\right).

Um vetor de espaço-tempo contravariante de 4 dimensões pode ser dado por xμ=(ct,x)=(x0,x1,x2,x3)x^{\mu}=\left(ct,\boldsymbol{x}\right)=\left(x^{0},x^{1},x^{2},x^{3}\right), lembrando que o vetor contravariante inverte o sinal das componentes epaciais quando comparado com o covariante. Isto é, temo:

xμ=gμνxνx_{\mu}=g_{\mu\nu}x^{\nu}

Onde o tensor métrico é responsável por subir ou baixar os índices:

gμν=(10000−10000−10000−1)g_{\mu\nu}=\left(\begin{array}{cccc} 1 & 0 & 0 & 0\\ 0 & -1 & 0 & 0\\ 0 & 0 & -1 & 0\\ 0 & 0 & 0 & -1 \end{array}\right)

O produto interno então é dado entre o vetor covariante e contravariante:

xμxμ=xμxμx_{\mu}x^{\mu}=x^{\mu}x_{\mu}

Dessa forma temos que D(X)=gμνxμxν=xνxν=xνxνD\left(X\right)=g_{\mu\nu}x^{\mu}x^{\nu}=x_{\nu}x^{\nu}=x^{\nu}x_{\nu}. Definindo uma transformação linear do vetor:

x′μ=Λνμxνx^{'\mu}=\Lambda^{\mu}_{\nu}x^{\nu}

Temo então o produto escalar

x′μgμνx′ν=xβΛβμgαβΛαμxα=xμgαβxν\begin{align*} x^{'\mu}g_{\mu\nu}x^{'\nu} & =x^{\beta}\Lambda^{\mu}_{\beta}g_{\alpha\beta}\Lambda^{\mu}_{\alpha}x^{\alpha}\\ & =x^{\mu}g_{\alpha\beta}x^{\nu} \end{align*}

Se demandamos que seja preservado o produto interno esalar de Minkowski:

ΛαμgμνΛβμ=gαβ\Lambda^{\mu}_{\alpha}g_{\mu\nu}\Lambda^{\mu}_{\beta}=g_{\alpha\beta}

Então agora voltando às nossas matrizes, podemos ver que, utilizando o Wolfram Mathematica:

Det{[}\{\{x0~+~x3,~x1~-~I~x2\},~\{x1~+~I~x2,~x0~-~x3\}\}{]}
D(X)=xνxν=xνgμνxμ=(x0)2−(x1)2−(x2)2−(x3)2\begin{align*} D\left(X\right) & =x_{\nu}x^{\nu}\\ & =x_{\nu}g_{\mu\nu}x_{\mu}\\ & =\left(x^{0}\right)^{2}-\left(x^{1}\right)^{2}-\left(x^{2}\right)^{2}-\left(x^{3}\right)^{2} \end{align*}

E:

X~:=~\{\{x0~+~x3,~x1~-~I~x2\},~\{x1~+~I~x2,~x0~-~x3\}\};

W~:=~\{\{a,~b\},~\{g,~d\}\};

WC~:=~\{\{Conjugate{[}a{]},~Conjugate{[}g{]}\},\{Conjugate{[}b{]},~Conjugate{[}d{]}\}\};

Simplify{[}Det{[}W~.~X~.~WC{]}{]}
D(X′)=xν′x′νD(X′)=xν′gμνxμ′D(WXW∗)=xβΛβμgμνΛαμxα=(αδ−βγ)[(x0)2−(x1)2−(x2)2−(x3)2](α‾δ‾−β‾γ‾)\begin{align*} D\left(X'\right) & =x^{'}_{\nu}x^{'\nu}\\ D\left(X'\right) & =x'_{\nu}g_{\mu\nu}x'_{\mu}\\ D\left(WXW^{*}\right) & =x^{\beta}\Lambda^{\mu}_{\beta}g_{\mu\nu}\Lambda^{\mu}_{\alpha}x^{\alpha}\\ & =\left(\alpha\delta-\beta\gamma\right)\left[\left(x^{0}\right)^{2}-\left(x^{1}\right)^{2}-\left(x^{2}\right)^{2}-\left(x^{3}\right)^{2}\right]\left(\overline{\alpha}\overline{\delta}-\overline{\beta}\overline{\gamma}\right) \end{align*}

Temos então, como queríamos demonstrar, que:

D(X′)=(x0)2−(x1)2−(x2)2−(x3)2=D(X)D\left(X'\right)=\left(x^{0}\right)^{2}-\left(x^{1}\right)^{2}-\left(x^{2}\right)^{2}-\left(x^{3}\right)^{2}=D\left(X\right)

Tendo então X\mathbb{X} representando uma coleção de matrizes hermitianas 2×22\times2:

X={W,∈Mat2×2(C):X=X†}\mathbb{X}=\left\{ W,\in Mat_{2\times2}\left(\mathbb{C}\right):X=X^{\dagger}\right\}

Isso pode ser considerado uma forma alternativa de R1,3\mathbb{R}^{1,3} através do mapa que nos foi dado:

R1,3→X,(x0,x1,x2,x3)↦X=(x0+x3x1−ix2x1+ix2x0−x3)\mathbb{R}^{1,3}\rightarrow\mathbb{X},\quad\left(x^{0},x^{1},x^{2},x^{3}\right)\mapsto X=\left(\begin{array}{cc} x^{0}+x^{3} & x^{1}-ix^{2}\\ x^{1}+ix^{2} & x^{0}-x^{3} \end{array}\right)

Lembrando que det⁡X=xνxν\det X=x_{\nu}x^{\nu}, e então a base padrão para R1,3\mathbb{R}^{1,3} pode ser escrita para X\mathbb{X} usando as matrizes de pauli. Com essa base, então um ponto (x0,x1,x2,x3)\left(x^{0},x^{1},x^{2},x^{3}\right) no espaço-tempo de Minkowsku corresponde a uma matriz 2×22\times2 , X=xμσμX=x^{\mu}\sigma_{\mu}. Tendo o produto interno de Frobenius dado por:

⟨A,B⟩=Tr(A†B)\left\langle A,B\right\rangle =Tr\left(A^{\dagger}B\right)

E a norma de Frobenius ∥A∥=⟨A,A⟩\left\Vert A\right\Vert =\sqrt{\left\langle A,A\right\rangle } Então é fácil checar que a base é ortonormal:

⟨σμ,σν⟩=δμν\left\langle \sigma_{\mu},\sigma_{\nu}\right\rangle =\delta_{\mu\nu}

Desde que definimos o produto interno como: ⟨A,B⟩=12Tr(A†B)\left\langle A,B\right\rangle =\frac{1}{2}Tr\left(A^{\dagger}B\right)Lima (2013). Temos então uma forma simples de escrever o mapa inverso de X\mathbb{X} para R1,3\mathbb{R}^{1,3}.

X↦(⟨X,σ0⟩,⟨X,σ1⟩,⟨X,σ2⟩,⟨X,σ3⟩)X\mapsto\left(\left\langle X,\sigma_{0}\right\rangle ,\left\langle X,\sigma_{1}\right\rangle ,\left\langle X,\sigma_{2}\right\rangle ,\left\langle X,\sigma_{3}\right\rangle \right)

Calculando para X=xμσμX=x^{\mu}\sigma_{\mu}, podemos notar facilmente que:

X↦(x0,x1,x2,x3)X\mapsto\left(x^{0},x^{1},x^{2},x^{3}\right)

Então temos que:

X′↦(x′0,x′1,x′2,x′3)↦(Λν0xν,Λν1xν,Λν2xν,Λν3xν)↦(⟨X′,σ0⟩,⟨X′,σ1⟩,⟨X′,σ2⟩,⟨X′,σ3⟩)↦(⟨WXW∗,σ0⟩,⟨WXW∗,σ1⟩,⟨WXW∗,σ2⟩,⟨WXW∗,σ3⟩)\begin{align*} X' & \mapsto\left(x^{'0},x^{'1},x^{'2},x^{'3}\right)\\ & \mapsto\left(\Lambda^{0}_{\nu}x^{\nu},\Lambda^{1}_{\nu}x^{\nu},\Lambda^{2}_{\nu}x^{\nu},\Lambda^{3}_{\nu}x^{\nu}\right)\\ & \mapsto\left(\left\langle X',\sigma_{0}\right\rangle ,\left\langle X',\sigma_{1}\right\rangle ,\left\langle X',\sigma_{2}\right\rangle ,\left\langle X',\sigma_{3}\right\rangle \right)\\ & \mapsto\left(\left\langle WXW^{*},\sigma_{0}\right\rangle ,\left\langle WXW^{*},\sigma_{1}\right\rangle ,\left\langle WXW^{*},\sigma_{2}\right\rangle ,\left\langle WXW^{*},\sigma_{3}\right\rangle \right) \end{align*}

Abrindo então o primero termo:

Λν0xν=⟨WXW∗,σ0⟩=⟨Wσ0W∗x0+Wσ1W∗x1+Wσ2W∗x2+Wσ3W∗x3,σ0⟩=⟨Wσ0W∗x0+Wσ1W∗x1+Wσ2W∗x2+Wσ3W∗x3,σ0⟩=⟨Wσ0W∗,σ0⟩x0+⟨Wσ1W∗,σ0⟩x1+⟨Wσ2W∗,σ0⟩x2+⟨Wσ3W∗,σ0⟩x3=⟨WσνW∗,σ0⟩xν\begin{align*} \Lambda^{0}_{\nu}x^{\nu} & =\left\langle WXW^{*},\sigma_{0}\right\rangle \\ & =\left\langle W\sigma_{0}W^{*}x^{0}+W\sigma_{1}W^{*}x^{1}+W\sigma_{2}W^{*}x^{2}+W\sigma_{3}W^{*}x^{3},\sigma_{0}\right\rangle \\ & =\left\langle W\sigma_{0}W^{*}x^{0}+W\sigma_{1}W^{*}x^{1}+W\sigma_{2}W^{*}x^{2}+W\sigma_{3}W^{*}x^{3},\sigma_{0}\right\rangle \\ & =\left\langle W\sigma_{0}W^{*},\sigma_{0}\right\rangle x^{0}+\left\langle W\sigma_{1}W^{*},\sigma_{0}\right\rangle x^{1}+\left\langle W\sigma_{2}W^{*},\sigma_{0}\right\rangle x^{2}+\left\langle W\sigma_{3}W^{*},\sigma_{0}\right\rangle x^{3}\\ & =\left\langle W\sigma_{\nu}W^{*},\sigma_{0}\right\rangle x^{\nu} \end{align*}

Temos então que:

Λν0=⟨WσνW∗,σ0⟩\Lambda^{0}_{\nu}=\left\langle W\sigma_{\nu}W^{*},\sigma_{0}\right\rangle

Ou ainda mais geral:

Λνμ=⟨WσνW∗,σμ⟩\Lambda^{\mu}_{\nu}=\left\langle W\sigma_{\nu}W^{*},\sigma_{\mu}\right\rangle

Logo:

Λ=(Λ00Λ10Λ20Λ30Λ01Λ11Λ21Λ31Λ02Λ12Λ22Λ32Λ03Λ13Λ23Λ33)=(⟨Wσ0W∗,σ0⟩⟨Wσ1W∗,σ0⟩⟨Wσ2W∗,σ0⟩⟨Wσ3W∗,σ0⟩⟨Wσ0W∗,σ1⟩⟨Wσ1W∗,σ1⟩⟨Wσ2W∗,σ1⟩⟨Wσ3W∗,σ1⟩⟨Wσ0W∗,σ2⟩⟨Wσ1W∗,σ2⟩⟨Wσ2W∗,σ2⟩⟨Wσ3W∗,σ2⟩⟨Wσ0W∗,σ3⟩⟨Wσ1W∗,σ3⟩⟨Wσ2W∗,σ3⟩⟨Wσ3W∗,σ3⟩)\Lambda=\left(\begin{array}{cccc} \Lambda^{0}_{0} & \Lambda^{0}_{1} & \Lambda^{0}_{2} & \Lambda^{0}_{3}\\ \Lambda^{1}_{0} & \Lambda^{1}_{1} & \Lambda^{1}_{2} & \Lambda^{1}_{3}\\ \Lambda^{2}_{0} & \Lambda^{2}_{1} & \Lambda^{2}_{2} & \Lambda^{2}_{3}\\ \Lambda^{3}_{0} & \Lambda^{3}_{1} & \Lambda^{3}_{2} & \Lambda^{3}_{3} \end{array}\right)=\left(\begin{array}{cccc} \left\langle W\sigma_{0}W^{*},\sigma_{0}\right\rangle & \left\langle W\sigma_{1}W^{*},\sigma_{0}\right\rangle & \left\langle W\sigma_{2}W^{*},\sigma_{0}\right\rangle & \left\langle W\sigma_{3}W^{*},\sigma_{0}\right\rangle \\ \left\langle W\sigma_{0}W^{*},\sigma_{1}\right\rangle & \left\langle W\sigma_{1}W^{*},\sigma_{1}\right\rangle & \left\langle W\sigma_{2}W^{*},\sigma_{1}\right\rangle & \left\langle W\sigma_{3}W^{*},\sigma_{1}\right\rangle \\ \left\langle W\sigma_{0}W^{*},\sigma_{2}\right\rangle & \left\langle W\sigma_{1}W^{*},\sigma_{2}\right\rangle & \left\langle W\sigma_{2}W^{*},\sigma_{2}\right\rangle & \left\langle W\sigma_{3}W^{*},\sigma_{2}\right\rangle \\ \left\langle W\sigma_{0}W^{*},\sigma_{3}\right\rangle & \left\langle W\sigma_{1}W^{*},\sigma_{3}\right\rangle & \left\langle W\sigma_{2}W^{*},\sigma_{3}\right\rangle & \left\langle W\sigma_{3}W^{*},\sigma_{3}\right\rangle \end{array}\right)

Ou como Tr(AB)=Tr(BA)Tr\left(AB\right)=Tr\left(BA\right)

Λ=(⟨σ0,Wσ0W∗⟩⟨σ0,Wσ1W∗⟩⟨σ0,Wσ2W∗⟩⟨σ0,Wσ3W∗⟩⟨σ1,Wσ0W∗⟩⟨σ1,Wσ1W∗⟩⟨σ1,Wσ2W∗⟩⟨σ1,Wσ3W∗⟩⟨σ2,Wσ0W∗⟩⟨σ2,Wσ1W∗⟩⟨σ2,Wσ2W∗⟩⟨σ2,Wσ3W∗⟩⟨σ3,Wσ0W∗⟩⟨σ3,Wσ1W∗⟩⟨σ3,Wσ2W∗⟩⟨σ3,Wσ3W∗⟩)\Lambda=\left(\begin{array}{cccc} \left\langle \sigma_{0},W\sigma_{0}W^{*}\right\rangle & \left\langle \sigma_{0},W\sigma_{1}W^{*}\right\rangle & \left\langle \sigma_{0},W\sigma_{2}W^{*}\right\rangle & \left\langle \sigma_{0},W\sigma_{3}W^{*}\right\rangle \\ \left\langle \sigma_{1},W\sigma_{0}W^{*}\right\rangle & \left\langle \sigma_{1},W\sigma_{1}W^{*}\right\rangle & \left\langle \sigma_{1},W\sigma_{2}W^{*}\right\rangle & \left\langle \sigma_{1},W\sigma_{3}W^{*}\right\rangle \\ \left\langle \sigma_{2},W\sigma_{0}W^{*}\right\rangle & \left\langle \sigma_{2},W\sigma_{1}W^{*}\right\rangle & \left\langle \sigma_{2},W\sigma_{2}W^{*}\right\rangle & \left\langle \sigma_{2},W\sigma_{3}W^{*}\right\rangle \\ \left\langle \sigma_{3},W\sigma_{0}W^{*}\right\rangle & \left\langle \sigma_{3},W\sigma_{1}W^{*}\right\rangle & \left\langle \sigma_{3},W\sigma_{2}W^{*}\right\rangle & \left\langle \sigma_{3},W\sigma_{3}W^{*}\right\rangle \end{array}\right)

Usando o Wolfram Mathematica, podemos calcular o determinante:

O0~:=~\{\{1,~0\},~~\{0,~1\}\};~

O1~:=~\{\{0,~1\},~~\{1,~0\}\};~

O2~:=~\{\{0,~-I\},~\{I,~0\}\};~

O3~:=~\{\{1,~0\},~\{0,~-1\}\};

Simplify{[}Det{[}\{

\{Tr{[}O0.W.O0.WC{]}/2,Tr{[}O0.W.O1.WC{]}/2,~Tr{[}O0.W.O2.WC{]}/2,Tr{[}O0.W.O3.WC{]}/2\},

\{Tr{[}O1.W.O0.WC{]}/2,Tr{[}O1.W.O1.WC{]}/2,~Tr{[}O1.W.O2.WC{]}/2,Tr{[}O1.W.O3.WC{]}/2\},~

\{Tr{[}O2.W.O0.WC{]}/2,Tr{[}O2.W.O1.WC{]}/2,~Tr{[}O2.W.O2.WC{]}/2,Tr{[}O2.W.O3.WC{]}/2\},~

\{Tr{[}O3.W.O0.WC{]}/2,Tr{[}O3.W.O1.WC{]}/2,~Tr{[}O3.W.O2.WC{]}/2,Tr{[}O3.W.O3.WC{]}/2\}\}{]}{]}
det⁡(Λ)=(αδ−βγ)2(α‾δ‾−β‾γ‾)2=+1\det\left(\Lambda\right)=\left(\alpha\delta-\beta\gamma\right)^{2}\left(\overline{\alpha}\overline{\delta}-\overline{\beta}\overline{\gamma}\right)^{2}=+1

Como queríamoss demonstrar. Agora calculando o primeiro termo:

Tr{[}O0.W.O0.WC{]}/2
Λ00=,Wσ0W∗=12Tr(σ0Wσ0W∗)=12(∣α∣2+∣δ∣2+∣β∣2+∣γ∣2)=12(∣α∣2+∣δ∣2+∣β∣2+∣γ∣2)\begin{align*} \Lambda^{0}_{0} & =,W\sigma_{0}W^{*}\\ & =\frac{1}{2}Tr\left(\sigma_{0}W\sigma_{0}W^{*}\right)\\ & =\frac{1}{2}\left(\left|\alpha\right|^{2}+\left|\delta\right|^{2}+\left|\beta\right|^{2}+\left|\gamma\right|^{2}\right)\\ & =\frac{1}{2}\left(\left|\alpha\right|^{2}+\left|\delta\right|^{2}+\left|\beta\right|^{2}+\left|\gamma\right|^{2}\right) \end{align*}

Como vimos em aula 32, as constantes Λβα\Lambda^{\alpha}_{\beta} é estrita as condições:

ΛΓαΛγβηαβ=ηΓγΛ0αΛ0βηαβ=η00(Λ00)2+(Λ01)2+(Λ02)2+(Λ03)=1(Λ00)2=1+∑i=13(Λ0i)2\begin{align*} \Lambda^{\alpha}_{\varGamma} & \Lambda^{\beta}_{\gamma}\eta_{\alpha\beta}=\eta_{\varGamma\gamma}\\ \Lambda^{\alpha}_{0} & \Lambda^{\beta}_{0}\eta_{\alpha\beta}=\eta_{00}\\ \left(\Lambda^{0}_{0}\right)^{2}+\left(\Lambda^{1}_{0}\right)^{2}+\left(\Lambda^{2}_{0}\right)^{2}+\left(\Lambda^{3}_{0}\right) & =1\\ \left(\Lambda^{0}_{0}\right)^{2}= & 1+\sum^{3}_{i=1}\left(\Lambda^{i}_{0}\right)^{2} \end{align*}

Onde usamos γ=Γ=0\gamma=\Gamma=0 e ηαβ\eta_{\alpha\beta} é a “métrica e minkowski”. Calulando então o somatório:

FullSimplify{[}Expand{[}(Tr{[}O1.W.O0.WC{]}/2)\textasciicircum 2+(Tr{[}O2.W.O0.WC{]}/2)\textasciicircum 2+(Tr{[}O3.W.O0.WC{]}/2)\textasciicircum 2{]}{]}
∑i=13(Λ0i)2=14((∣α∣4+∣β∣4+∣δ∣4+∣γ∣4+2(α‾(−δ‾(αδ−βγ−βγ)+α∣β∣2+α∣γ∣2)+β‾(γ‾(αδ+αδ−βγ)+β∣δ∣2)+∣δ∣2∣γ∣2)))=14((∣α∣4+∣β∣4+∣δ∣4+∣γ∣4+2(α‾δ‾βγ−(δ‾α‾−γ‾β‾)+αδγ‾β‾+∣α∣2∣γ∣2+∣α∣2∣β∣2+∣β∣2∣δ∣2+∣δ∣2∣γ∣2)))=14((∣α∣4+∣β∣4+∣δ∣4+∣γ∣4+2(α‾δ‾βγ+αδγ‾β‾−1+∣α∣2∣γ∣2+∣α∣2∣β∣2+∣β∣2∣δ∣2+∣δ∣2∣γ∣2)))=14((∣α∣4+∣β∣4+∣δ∣4+∣γ∣4+2(∣α∣2∣δ∣2+∣β∣2∣γ∣2−2+∣α∣2∣γ∣2+∣α∣2∣β∣2+∣β∣2∣δ∣2+∣δ∣2∣γ∣2)))=−1+14(∣α∣4+2∣α∣2∣δ∣2+2∣α∣2∣β∣2+2∣α∣2∣γ∣2+∣δ∣4+2∣β∣2∣δ∣2+2∣δ∣2∣γ∣2+∣β∣4+2∣β∣2∣γ∣2+∣γ∣4)=−1+[12(∣α∣2+∣δ∣2+∣β∣2+∣γ∣2)]2\begin{align*} \sum^{3}_{i=1}\left(\Lambda^{i}_{0}\right)^{2} & =\frac{1}{4}\left(\left(\left|\alpha\right|^{4}+\left|\beta\right|^{4}+\left|\delta\right|^{4}+\left|\gamma\right|^{4}+2\left(\overline{\alpha}\left(-\overline{\delta}(\alpha\delta-\beta\gamma-\beta\gamma)+\alpha\left|\beta\right|^{2}+\alpha\left|\gamma\right|^{2}\right)+\overline{\beta}\left(\overline{\gamma}(\alpha\delta+\alpha\delta-\beta\gamma)+\beta\left|\delta\right|^{2}\right)+\left|\delta\right|^{2}\left|\gamma\right|^{2}\right)\right)\right)\\ & =\frac{1}{4}\left(\left(\left|\alpha\right|^{4}+\left|\beta\right|^{4}+\left|\delta\right|^{4}+\left|\gamma\right|^{4}+2\left(\overline{\alpha}\overline{\delta}\beta\gamma-\left(\overline{\delta}\overline{\alpha}-\overline{\gamma}\overline{\beta}\right)+\alpha\delta\overline{\gamma}\overline{\beta}+\left|\alpha\right|^{2}\left|\gamma\right|^{2}+\left|\alpha\right|^{2}\left|\beta\right|^{2}+\left|\beta\right|^{2}\left|\delta\right|^{2}+\left|\delta\right|^{2}\left|\gamma\right|^{2}\right)\right)\right)\\ & =\frac{1}{4}\left(\left(\left|\alpha\right|^{4}+\left|\beta\right|^{4}+\left|\delta\right|^{4}+\left|\gamma\right|^{4}+2\left(\overline{\alpha}\overline{\delta}\beta\gamma+\alpha\delta\overline{\gamma}\overline{\beta}-1+\left|\alpha\right|^{2}\left|\gamma\right|^{2}+\left|\alpha\right|^{2}\left|\beta\right|^{2}+\left|\beta\right|^{2}\left|\delta\right|^{2}+\left|\delta\right|^{2}\left|\gamma\right|^{2}\right)\right)\right)\\ & =\frac{1}{4}\left(\left(\left|\alpha\right|^{4}+\left|\beta\right|^{4}+\left|\delta\right|^{4}+\left|\gamma\right|^{4}+2\left(\left|\alpha\right|^{2}\left|\delta\right|^{2}+\left|\beta\right|^{2}\left|\gamma\right|^{2}-2+\left|\alpha\right|^{2}\left|\gamma\right|^{2}+\left|\alpha\right|^{2}\left|\beta\right|^{2}+\left|\beta\right|^{2}\left|\delta\right|^{2}+\left|\delta\right|^{2}\left|\gamma\right|^{2}\right)\right)\right)\\ & =-1+\frac{1}{4}\left(\left|\alpha\right|^{4}+2\left|\alpha\right|^{2}\left|\delta\right|^{2}+2\left|\alpha\right|^{2}\left|\beta\right|^{2}+2\left|\alpha\right|^{2}\left|\gamma\right|^{2}+\left|\delta\right|^{4}+2\left|\beta\right|^{2}\left|\delta\right|^{2}+2\left|\delta\right|^{2}\left|\gamma\right|^{2}+\left|\beta\right|^{4}+2\left|\beta\right|^{2}\left|\gamma\right|^{2}+\left|\gamma\right|^{4}\right)\\ & =-1+\left[\frac{1}{2}\left(\left|\alpha\right|^{2}+\left|\delta\right|^{2}+\left|\beta\right|^{2}+\left|\gamma\right|^{2}\right)\right]^{2} \end{align*}

Onde usamos o fato de que −1=(αδ−βγ)(γ‾β‾−α‾δ‾)=α‾δ‾βγ−∣α∣2∣δ∣2−∣β∣2∣γ∣2+αδγ‾β‾-1=\left(\alpha\delta-\beta\gamma\right)\left(\overline{\gamma}\overline{\beta}-\overline{\alpha}\overline{\delta}\right)=\overline{\alpha}\overline{\delta}\beta\gamma-\left|\alpha\right|^{2}\left|\delta\right|^{2}-\left|\beta\right|^{2}\left|\gamma\right|^{2}+\alpha\delta\overline{\gamma}\overline{\beta}. Temos então a seguinte desigualdade:

∑i=13(Λ0i)2≥0⟶[12(∣α∣2+∣δ∣2+∣β∣2+∣γ∣2)]2≥1\sum^{3}_{i=1}\left(\Lambda^{i}_{0}\right)^{2}\geq0\longrightarrow\left[\frac{1}{2}\left(\left|\alpha\right|^{2}+\left|\delta\right|^{2}+\left|\beta\right|^{2}+\left|\gamma\right|^{2}\right)\right]^{2}\geq1

Então:

∣12(∣α∣2+∣δ∣2+∣β∣2+∣γ∣2)∣≥1\left|\frac{1}{2}\left(\left|\alpha\right|^{2}+\left|\delta\right|^{2}+\left|\beta\right|^{2}+\left|\gamma\right|^{2}\right)\right|\geq1

Uma vez que se ∣x∣>1(∣x∣<1),\left|x\right|>1\left(\left|x\right|<1\right),então x2>1(x2<1)x^{2}>1\left(x^{2}<1\right). E lembrando que

\[
\Lambda^{0}_{0}=\frac{1}{2}\left(\left|\alpha\right|^{2}+\left|\delta\right|^{2}+\left|\beta\right|^{2}+\left|\gamma\right|^{2}\right)
\]

Concluímos então que:

∣Λ00∣≥1\left|\Lambda^{0}_{0}\right|\geq1

Problem 2

Operate on (iγμ∂μ−m)ψ=0\left(i\gamma^{\mu}\partial_{\mu}-m\right)\psi=0 with γμ∂μ\gamma^{\mu}\partial_{\mu} and show that each of the four components ψi\psi_{i} satisfies the Klein-Gordon equation (□2+m2)ψi=0\left(\square^{2}+m^{2}\right)\psi_{i}=0.

Primeiro temos que o oprador d’Alembertiano é dado por:

□2=ηαβ∂∂xβ∂∂xα=∇2−∂2∂t2\square^{2}=\eta^{\alpha\beta}\frac{\partial}{\partial x^{\beta}}\frac{\partial}{\partial x^{\alpha}}=\nabla^{2}-\frac{\partial^{2}}{\partial t^{2}}

Sendo a função onda ψ\psi uma matriz:

(iγμ∂μ−m)ψ=0\left(i\gamma^{\mu}\partial_{\mu}-m\right)\psi=0

Vamos operar γν∂ν\gamma^{\nu}\partial_{\nu}. Mas primeiro temos que:

∂μ=∂∂xμ=(∂∂t,∇)\partial_{\mu}=\frac{\partial}{\partial x^{\mu}}=\left(\frac{\partial}{\partial t},\boldsymbol{\nabla}\right)

E como temos

p→−iℏ∇,E=iℏ∂∂t,pμ=(E,−p)\boldsymbol{p}\rightarrow-i\hbar\nabla,\quad E=i\hbar\frac{\partial}{\partial t},\qquad p_{\mu}=\left(E,-\boldsymbol{p}\right)

Então vamos escrever γνi∂ν=γνpν/ℏ\gamma^{\nu}i\partial_{\nu}=\gamma^{\nu}p_{\nu}/\hbar, sendo ℏ=1\hbar=1. Temos então:

(iγμ∂μ−m)ψ=0(γμpμ−m)ψ=0γμpμψ=−mψ\begin{align*} \left(i\gamma^{\mu}\partial_{\mu}-m\right)\psi & =0\\ \left(\gamma^{\mu}p_{\mu}-m\right)\psi & =0\\ \gamma^{\mu}p_{\mu}\psi & =-m\psi \end{align*}

E ainda sendo γμpμ=p\gamma^{\mu}p_{\mu}=\cancel{p}:

pψ=−mψ\cancel{p}\psi=-m\psi

Operando então: −ip-i\cancel{p}:

ippψ=−ipmψ\begin{align*} i\cancel{p}\cancel{p} & \psi=-i\cancel{p}m\psi \end{align*}

E uma vez que na aula 37, na equação 159 vimos que γμγη+γηγμ=I2gμη\gamma^{\mu}\gamma^{\eta}+\gamma^{\eta}\gamma^{\mu}=\mathbb{I}2g^{\mu\eta}:

pp=γμpμγνpν=γμγνpμpν=12(γμγν+γνγμ)pμpν=I122gμηpμpν=Ipμpν=Ip2\begin{align*} \cancel{p}\cancel{p} & =\gamma^{\mu}p_{\mu}\gamma^{\nu}p_{\nu}\\ & =\gamma^{\mu}\gamma^{\nu}p_{\mu}p_{\nu}\\ & =\frac{1}{2}\left(\gamma^{\mu}\gamma^{\nu}+\gamma^{\nu}\gamma^{\mu}\right)p_{\mu}p_{\nu}\\ & =\mathbb{I}\frac{1}{2}2g^{\mu\eta}p_{\mu}p_{\nu}\\ & =\mathbb{I}p^{\mu}p_{\nu}\\ & =\mathbb{I}p^{2} \end{align*}

Lembando que [pi,pj]=0\left[p_{i},p_{j}\right]=0, o que nos permite escrever que γμγνpμpν=γνγμpνpμ=γνγμpμpν\gamma^{\mu}\gamma^{\nu}p_{\mu}p_{\nu}=\gamma^{\nu}\gamma^{\mu}p_{\nu}p_{\mu}=\gamma^{\nu}\gamma^{\mu}p_{\mu}p_{\nu}.Também usamos o fato de que as matrizes γν\gamma^{\nu} são constantes. E usando o próprio resultado obtido em 42, então:

p2ψ=m2ψpμpμψ=m2ψ(i∂μ)(i∂μ)ψ=m2ψ−(∂2∂t2−∇2)ψ=m2ψ−□2ψ=m2ψ\begin{align*} p^{2}\psi & =m^{2}\psi\\ p^{\mu}p_{\mu}\psi & =m^{2}\psi\\ \left(i\partial^{\mu}\right)\left(i\partial_{\mu}\right)\psi & =m^{2}\psi\\ -\left(\frac{\partial^{2}}{\partial t^{2}}-\nabla^{2}\right)\psi & =m^{2}\psi\\ -\square^{2}\psi & =m^{2}\psi \end{align*}

Ou então, como queríamos mostrar:

(□2+m2)ψ=0\left(\square^{2}+m^{2}\right)\psi=0

Ou para cada uma das componentes de ψ\psi:

(□2+m2)ψi=0\left(\square^{2}+m^{2}\right)\psi_{i}=0

Problem 3

Show that the helicity operator

h^=∑^⋅p^≡(σ⋅p^00σ⋅p^)\widehat{h}=\widehat{\sum}\cdot\widehat{\boldsymbol{p}}\equiv\left(\begin{array}{cc} \boldsymbol{\sigma}\cdot\widehat{\boldsymbol{p}} & 0\\ 0 & \boldsymbol{\sigma}\cdot\widehat{\boldsymbol{p}} \end{array}\right)

commutes with the Dirac Hamiltonian H^=α⋅p+βm\widehat{H}=\alpha\cdot\boldsymbol{p}+\beta m. Here p^\widehat{\boldsymbol{p}} is the unit vector pointing in the direction of the momentum, p/∣p∣\boldsymbol{p}/\left|\boldsymbol{p}\right|. The “spin” component in the direction of motion, 12σ⋅p^\frac{1}{2}\boldsymbol{\sigma}\cdot\widehat{\boldsymbol{p}}, , is therefore a “good” quantum number.

O operador de helicidade é dado por:

h^=∑^⋅p^≡(σ⋅p^00σ⋅p^)\widehat{h}=\widehat{\sum}\cdot\widehat{\boldsymbol{p}}\equiv\left(\begin{array}{cc} \boldsymbol{\sigma}\cdot\widehat{\boldsymbol{p}} & 0\\ 0 & \boldsymbol{\sigma}\cdot\widehat{\boldsymbol{p}} \end{array}\right)

Onde p^=p/∣p∣\widehat{\boldsymbol{p}}=\boldsymbol{p}/\left|\boldsymbol{p}\right|. Queremos mostrar que h^\widehat{h} comuat com o Hamiltoniano de Dirac H^=α⋅p+βm\widehat{H}=\alpha\cdot\boldsymbol{p}+\beta m.

Temos da aula 35 e equações 101 e 102 que:

αi=(0σiσi0),β=(100−1)\alpha_{i}=\left(\begin{array}{cc} 0 & \sigma_{i}\\ \sigma_{i} & 0 \end{array}\right),\qquad\beta=\left(\begin{array}{cc} \boldsymbol{1} & 0\\ 0 & -\boldsymbol{1} \end{array}\right)

Então:

H^=α⋅p+βm=(0σ⋅pσ⋅p0)+(m100−m1)=(mσ⋅pσ⋅pm)\begin{align*} \widehat{H} & =\alpha\cdot\boldsymbol{p}+\beta m=\left(\begin{array}{cc} 0 & \boldsymbol{\sigma}\cdot\boldsymbol{p}\\ \boldsymbol{\sigma}\cdot\boldsymbol{p} & 0 \end{array}\right)+\left(\begin{array}{cc} m\boldsymbol{1} & 0\\ 0 & -m\boldsymbol{1} \end{array}\right)=\left(\begin{array}{cc} \boldsymbol{m} & \boldsymbol{\sigma}\cdot\boldsymbol{p}\\ \boldsymbol{\sigma}\cdot\boldsymbol{p} & \boldsymbol{m} \end{array}\right) \end{align*}

Então podemos verificar explicitamente que:

[H^,h^]=(mσ⋅pσ⋅pm)(σ⋅p^00σ⋅p^)−(σ⋅p^00σ⋅p^)(mσ⋅pσ⋅pm)=(m(σ⋅p^)(σ⋅p)(σ⋅p^)(σ⋅p)(σ⋅p^)m(σ⋅p^))−((σ⋅p^)m(σ⋅p^)(σ⋅p)(σ⋅p^)(σ⋅p)(σ⋅p^)m)\begin{align*} \left[\widehat{H},\widehat{h}\right] & =\left(\begin{array}{cc} \boldsymbol{m} & \boldsymbol{\sigma}\cdot\boldsymbol{p}\\ \boldsymbol{\sigma}\cdot\boldsymbol{p} & \boldsymbol{m} \end{array}\right)\left(\begin{array}{cc} \boldsymbol{\sigma}\cdot\widehat{\boldsymbol{p}} & 0\\ 0 & \boldsymbol{\sigma}\cdot\widehat{\boldsymbol{p}} \end{array}\right)-\left(\begin{array}{cc} \boldsymbol{\sigma}\cdot\widehat{\boldsymbol{p}} & 0\\ 0 & \boldsymbol{\sigma}\cdot\widehat{\boldsymbol{p}} \end{array}\right)\left(\begin{array}{cc} \boldsymbol{m} & \boldsymbol{\sigma}\cdot\boldsymbol{p}\\ \boldsymbol{\sigma}\cdot\boldsymbol{p} & \boldsymbol{m} \end{array}\right)\\ & =\left(\begin{array}{cc} \boldsymbol{m}\left(\boldsymbol{\sigma}\cdot\widehat{\boldsymbol{p}}\right) & \left(\boldsymbol{\sigma}\cdot\boldsymbol{p}\right)\left(\boldsymbol{\sigma}\cdot\widehat{\boldsymbol{p}}\right)\\ \left(\boldsymbol{\sigma}\cdot\boldsymbol{p}\right)\left(\boldsymbol{\sigma}\cdot\widehat{\boldsymbol{p}}\right) & \boldsymbol{m}\left(\boldsymbol{\sigma}\cdot\widehat{\boldsymbol{p}}\right) \end{array}\right)-\left(\begin{array}{cc} \left(\boldsymbol{\sigma}\cdot\widehat{\boldsymbol{p}}\right)\boldsymbol{m} & \left(\boldsymbol{\sigma}\cdot\widehat{\boldsymbol{p}}\right)\left(\boldsymbol{\sigma}\cdot\boldsymbol{p}\right)\\ \left(\boldsymbol{\sigma}\cdot\widehat{\boldsymbol{p}}\right)\left(\boldsymbol{\sigma}\cdot\boldsymbol{p}\right) & \left(\boldsymbol{\sigma}\cdot\widehat{\boldsymbol{p}}\right)\boldsymbol{m} \end{array}\right) \end{align*}

E como queríamos demonstrar:

[H^,h^]=0\left[\widehat{H},\widehat{h}\right]=0

Problem 4

We have already seen that the equation Hu=(α⋅p+βm)u=EuHu=\left(\boldsymbol{\alpha}\cdot\boldsymbol{p}+\beta m\right)u=Eu reduces, for p≠0\boldsymbol{p}\neq0, to

σ⋅puB=(Ec−mc)uAσ⋅puA=(Ec+mc)uB\begin{align*} \boldsymbol{\sigma}\cdot\boldsymbol{p}u_{B} & =\left(\frac{E}{c}-mc\right)u_{A}\\ \boldsymbol{\sigma}\cdot\boldsymbol{p}u_{A} & =\left(\frac{E}{c}+mc\right)u_{B} \end{align*}

where uu has been divided into two two-components spinors, uAu_{A} and uBu_{B}.

a) For a nonrelativistic electron of velocity vv, use (1) to show that uAu_{A} is larger than uBu_{B} by a factor of the order v/cv/c. In nonrelativistic problems, ψA\psi_{A} and ψB\psi_{B} are referred to as the “large” and “small” components of the electron wavefunction ψ\psi.

b) In the nonrelativistic limit, show that the Dirac equation for a electron (charge −e-e) in an electromagnetic field Aμ=(A0,A)A^{\mu}=\left(A^{0},\boldsymbol{A}\right) reduces to the Schrodinger-Pauli equation

(12m(P+eA)2+e2mσ⋅B−eA0)ψA=ENRψA\left(\frac{1}{2m}\left(\boldsymbol{P}+e\boldsymbol{A}\right)^{2}+\frac{e}{2m}\boldsymbol{\sigma}\cdot\boldsymbol{B}-eA^{0}\right)\psi_{A}=E_{NR}\psi_{A}

where the magnetic field B=∇×A\boldsymbol{B}=\nabla\times\boldsymbol{A} and ENR=E−mE_{NR}=E-m. Assume ∣eA0∣≪m\left|eA^{0}\right|\ll m.

Hint: Make the substitution Pμ→Pμ+eAμP^{\mu}\rightarrow P^{\mu}+eA^{\mu} in eqs. (1) written in terms of ψA,B\psi_{A,B}, eliminate ψB\psi_{B}, and use P×A+A×P=i∇×A\boldsymbol{P}\times\boldsymbol{A}+\boldsymbol{A}\times\boldsymbol{P}=i\nabla\times\boldsymbol{A} where P=−i∇\boldsymbol{P}=-i\nabla.

Não usando unidades naturais pois queremos a razão v/cv/c, temos então:

σ⋅puB=(Ec−mc)uAσ⋅puA=(Ec+mc)uB\begin{align} \boldsymbol{\sigma}\cdot\boldsymbol{p}u_{B} & =\left(\frac{E}{c}-mc\right)u_{A}\\ \boldsymbol{\sigma}\cdot\boldsymbol{p}u_{A} & =\left(\frac{E}{c}+mc\right)u_{B} \end{align}

Ou manipulando:

uA=σ⋅p(Ec−mc)uBuB=σ⋅p(Ec+mc)uA\begin{align} u_{A} & =\frac{\boldsymbol{\sigma}\cdot\boldsymbol{p}}{\left(\frac{E}{c}-mc\right)}u_{B}\\ u_{B} & =\frac{\boldsymbol{\sigma}\cdot\boldsymbol{p}}{\left(\frac{E}{c}+mc\right)}u_{A} \end{align}

Pegando então a segunda equação, e calculando o módulo de uBu_{B}:

∣uB∣=1(Ec+mc)∣(σ⋅p)uA∣\left|u_{B}\right|=\frac{1}{\left(\frac{E}{c}+mc\right)}\left|\left(\boldsymbol{\sigma}\cdot\boldsymbol{p}\right)u_{A}\right|

E abrindo:

σ⋅p=pkσk→(σ⋅p)⋅(σ⋅p)=pkσkpjσj=pkpj(σkσj)=pjpj=∣p∣2\boldsymbol{\sigma}\cdot\boldsymbol{p}=p_{k}\sigma_{k}\rightarrow\left(\boldsymbol{\sigma}\cdot\boldsymbol{p}\right)\cdot\left(\boldsymbol{\sigma}\cdot\boldsymbol{p}\right)=p_{k}\sigma_{k}p_{j}\sigma_{j}=p_{k}p_{j}\left(\sigma_{k}\sigma_{j}\right)=p_{j}p_{j}=\left|p\right|^{2}

E como ∣v∣=v⋅v\left|\boldsymbol{v}\right|=\sqrt{\boldsymbol{v}\cdot\boldsymbol{v}}, então temos que ∣σ⋅p∣=pjpj=∣p∣\left|\boldsymbol{\sigma}\cdot\boldsymbol{p}\right|=\sqrt{p_{j}p_{j}}=\left|\boldsymbol{p}\right|. E devido a multiplicidade do módulo:

∣(σ⋅p)uA∣=∣p∣∣uA∣\left|\left(\boldsymbol{\sigma}\cdot\boldsymbol{p}\right)u_{A}\right|=\left|\boldsymbol{p}\right|\left|u_{A}\right|

Então sendo E≈mc2E\approx mc^{2} e p=mvp=mv

∣uB∣≈mv(mc2c+mc)∣uA∣\left|u_{B}\right|\approx\frac{mv}{\left(\frac{mc^{2}}{c}+mc\right)}\left|u_{A}\right|

Ou como queríamos demonstrar:

∣uB∣≈12vc∣uA∣\left|u_{B}\right|\approx\frac{1}{2}\frac{v}{c}\left|u_{A}\right|

Alguma observações é que encontrei na literatura a razão entre a componentes após substituir pμ→pμ+eAμp^{\mu}\rightarrow p^{\mu}+eA^{\mu}, mas optei por ir por outro caminho para ver se encontrava a mesma solução, uma vez que acabam desprezando os termos eAμeA^{\mu} durante o desenvolvimento. Independentemente da nota, eu também fiquei um pouco pensativo pois sempre substituem EE na segunda equação em 55 e não em 56. Na formulação em 57 levaria a uma indeterminação e em 56 ficamos com σ⋅puB≈0\boldsymbol{\sigma}\cdot\boldsymbol{p}u_{B}\approx0 então assumindo que σ⋅p≠0\boldsymbol{\sigma}\cdot\boldsymbol{p}\neq0 temos que uB≈0u_{B}\approx0. Eu acredito que isso decorre do fato que fizemos uma aproximação na energia e como pudemoz constatar, uBu_{B} é a componente menor por v/cv/c e como estamos em um caso não relativístico v/c≈0v/c\approx0. Mas se estou errado, gostaria de um retorno.

b)

Agora voltando a equação 99 da aula 35, temos que:

H=α⋅p+βmH=\boldsymbol{\alpha}\cdot\boldsymbol{p}+\beta m
(iγμ∂μ−Im)ψ=0\left(i\gamma^{\mu}\partial_{\mu}-\boldsymbol{I}m\right)\psi=0

Lembrando que γμ≡(β,βα)\gamma^{\mu}\equiv\left(\beta,\beta\boldsymbol{\alpha}\right) e

αi=(0σiσi0),β=(100−1)\alpha_{i}=\left(\begin{array}{cc} 0 & \sigma_{i}\\ \sigma_{i} & 0 \end{array}\right),\qquad\beta=\left(\begin{array}{cc} \boldsymbol{1} & 0\\ 0 & -\boldsymbol{1} \end{array}\right)

Reescrevendo:

(iγμpμ−Im)ψ=0\left(i\gamma^{\mu}p_{\mu}-\boldsymbol{I}m\right)\psi=0

Fazendo a substituição pμ→pu+eAμp^{\mu}\rightarrow p^{u}+eA^{\mu}, onde Aμ=(A0,A)A^{\mu}=\left(A^{0},\boldsymbol{A}\right)

(iγμ(∂μ+eAμ)−Im)ψ=0(γμi∂μ+eγμAμ−Im)ψ=0(γν(i∂ν+eAν)−Im)ψ=−γ0(i∂0+eA0)ψ\begin{align*} \left(i\gamma^{\mu}\left(\partial_{\mu}+eA_{\mu}\right)-\boldsymbol{I}m\right)\psi & =0\\ \left(\gamma^{\mu}i\partial_{\mu}+e\gamma^{\mu}A_{\mu}-\boldsymbol{I}m\right)\psi & =0\\ \left(\gamma^{\nu}\left(i\partial_{\nu}+eA_{\nu}\right)-\boldsymbol{I}m\right)\psi & =-\gamma^{0}\left(i\partial_{0}+eA_{0}\right)\psi \end{align*}

Onde μ=0,1,2,3\mu=0,1,2,3 e ν=1,2,3\nu=1,2,3. Além disso temos que pμ=i∂μp_{\mu}=i\partial_{\mu}e E=i∂0E=i\partial_{0}. Então:

(βα⋅(p+eA)−Im)ψ=−(E+eA0)βψ((0σi⋅(p+eA)−σi⋅(p+eA)0)−(m00m))(uAuB)=−(E+eA0)(100−1)(uAuB)(−mσi⋅(p+eA)−σi⋅(p+eA)−m)(uAuB)=(E+eA0)(−uAuB)((E+eA0−m)σi⋅(p+eA)−σi⋅(p+eA)−(E+eA0+m))(uAuB)=(00)\begin{align*} \left(\beta\boldsymbol{\alpha}\cdot\left(\boldsymbol{p}+e\boldsymbol{A}\right)-\boldsymbol{I}m\right)\psi & =-\left(E+eA_{0}\right)\beta\psi\\ \left(\left(\begin{array}{cc} 0 & \boldsymbol{\sigma}_{i}\cdot\left(\boldsymbol{p}+e\boldsymbol{A}\right)\\ -\boldsymbol{\sigma}_{i}\cdot\left(\boldsymbol{p}+e\boldsymbol{A}\right) & 0 \end{array}\right)-\left(\begin{array}{cc} m & 0\\ 0 & m \end{array}\right)\right)\left(\begin{array}{c} u_{A}\\ u_{B} \end{array}\right) & =-\left(E+eA_{0}\right)\left(\begin{array}{cc} \boldsymbol{1} & 0\\ 0 & -\boldsymbol{1} \end{array}\right)\left(\begin{array}{c} u_{A}\\ u_{B} \end{array}\right)\\ \left(\begin{array}{cc} -m & \boldsymbol{\sigma}_{i}\cdot\left(\boldsymbol{p}+e\boldsymbol{A}\right)\\ -\boldsymbol{\sigma}_{i}\cdot\left(\boldsymbol{p}+e\boldsymbol{A}\right) & -m \end{array}\right)\left(\begin{array}{c} u_{A}\\ u_{B} \end{array}\right) & =\left(E+eA_{0}\right)\left(\begin{array}{c} -u_{A}\\ u_{B} \end{array}\right)\\ \left(\begin{array}{cc} \left(E+eA_{0}-m\right) & \boldsymbol{\sigma}_{i}\cdot\left(\boldsymbol{p}+e\boldsymbol{A}\right)\\ -\boldsymbol{\sigma}_{i}\cdot\left(\boldsymbol{p}+e\boldsymbol{A}\right) & -\left(E+eA_{0}+m\right) \end{array}\right)\left(\begin{array}{c} u_{A}\\ u_{B} \end{array}\right) & =\left(\begin{array}{c} 0\\ 0 \end{array}\right) \end{align*}

O que nos leva a um novo conjunto de equações:

σ⋅(p+eA)uB=(E+eA0−m)uAσ⋅(p+eA)uA=(E+eA0+m)uB\begin{align*} \boldsymbol{\sigma}\cdot\left(\boldsymbol{p}+e\boldsymbol{A}\right)u_{B} & =\left(E+eA^{0}-m\right)u_{A}\\ \boldsymbol{\sigma}\cdot\left(\boldsymbol{p}+e\boldsymbol{A}\right)u_{A} & =\left(E+eA^{0}+m\right)u_{B} \end{align*}

Substituindo uBu_{B} na segunda equação na primeira:

(E+eA0−m)uA=[σ⋅(p+eA)]2(E+eA0+m)uA\begin{align*} \left(E+eA^{0}-m\right)u_{A} & =\frac{\left[\boldsymbol{\sigma}\cdot\left(\boldsymbol{p}+e\boldsymbol{A}\right)\right]^{2}}{\left(E+eA^{0}+m\right)}u_{A} \end{align*}

Então reescrevendo ENR=E−mE_{NR}=E-m:

(ENR+eA0)uA=[σ⋅(p+eA)]2(ENR+eA0+2m)uA\begin{align*} \left(E_{NR}+eA^{0}\right)u_{A} & =\frac{\left[\boldsymbol{\sigma}\cdot\left(\boldsymbol{p}+e\boldsymbol{A}\right)\right]^{2}}{\left(E_{NR}+eA^{0}+2m\right)}u_{A} \end{align*}

E usando eA0≪meA^{0}\ll m:

1(ENR+eA0+2m)=12m[1(ENR2m+eA02m+1)]≈12m[1(ENR2m)+1]\begin{align*} \frac{1}{\left(E_{NR}+eA^{0}+2m\right)} & =\frac{1}{2m}\left[\frac{1}{\left(\frac{E_{NR}}{2m}+\frac{eA^{0}}{2m}+1\right)}\right]\approx\frac{1}{2m}\left[\frac{1}{\left(\frac{E_{NR}}{2m}\right)+1}\right] \end{align*}

E fazendo expansão em série de taylor com ENR2m=x\frac{E_{NR}}{2m}=x, uma vez que se E≈mE\approx m em unidades naturais, então ENR≈0E_{NR}\approx0:

1(ENR+eA0+2m)≈12m[1−ENR2m+(ENR2m)2+… ]\frac{1}{\left(E_{NR}+eA^{0}+2m\right)}\approx\frac{1}{2m}\left[1-\frac{E_{NR}}{2m}+\left(\frac{E_{NR}}{2m}\right)^{2}+\dots\right]

Pegando só o primeiro termo da expansão:

:

1(ENR+eA0+2m)≈12m\frac{1}{\left(E_{NR}+eA^{0}+2m\right)}\approx\frac{1}{2m}

E voltando:

(ENR+eA0)uA=12m[σ⋅(p+eA)]2uA\begin{align*} \left(E_{NR}+eA^{0}\right)u_{A} & =\frac{1}{2m}\left[\boldsymbol{\sigma}\cdot\left(\boldsymbol{p}+e\boldsymbol{A}\right)\right]^{2}u_{A} \end{align*}

Abrindo então o termo ao quadrado, e lembrando que pμ=−iℏ∂μp_{\mu}=-i\hbar\partial_{\mu}, σμ\sigma_{\mu} é costante e as componentes do vetor potencial A\boldsymbol{A} são escalares:

[σ⋅(p+eA)]2=[σipi+eσiAi]2=[σjpjσipi+eσjpjσiAi+eσjAjσipi+e2σjAjσiAi]=σjσi[pjpi+epjAi+eAjpi+e2AjAi]=σjσi[pjpi+e2AjAi]+σjσi[epjAi+eAjpi]=12(σjσi+σiσj)[pjpi+e2AjAi]+σjσie[pjAi+Ajpi]\begin{align*} \left[\boldsymbol{\sigma}\cdot\left(\boldsymbol{p}+e\boldsymbol{A}\right)\right]^{2} & =\left[\sigma_{i}p_{i}+e\sigma_{i}A_{i}\right]^{2}\\ & =\left[\sigma_{j}p_{j}\sigma_{i}p_{i}+e\sigma_{j}p_{j}\sigma_{i}A_{i}+e\sigma_{j}A_{j}\sigma_{i}p_{i}+e^{2}\sigma_{j}A_{j}\sigma_{i}A_{i}\right]\\ & =\sigma_{j}\sigma_{i}\left[p_{j}p_{i}+ep_{j}A_{i}+eA_{j}p_{i}+e^{2}A_{j}A_{i}\right]\\ & =\sigma_{j}\sigma_{i}\left[p_{j}p_{i}+e^{2}A_{j}A_{i}\right]+\sigma_{j}\sigma_{i}\left[ep_{j}A_{i}+eA_{j}p_{i}\right]\\ & =\frac{1}{2}\left(\sigma_{j}\sigma_{i}+\sigma_{i}\sigma_{j}\right)\left[p_{j}p_{i}+e^{2}A_{j}A_{i}\right]+\sigma_{j}\sigma_{i}e\left[p_{j}A_{i}+A_{j}p_{i}\right] \end{align*}

Como {σi,σj}=2δij\left\{ \sigma_{i},\sigma_{j}\right\} =2\delta_{ij} e σiσj=δij+iεijkσk\sigma_{i}\sigma_{j}=\delta_{ij}+i\varepsilon_{ijk}\sigma_{k}:

[σ⋅(p+eA)]2=δij[pjpi+e2AjAi]+σjσie[pjAi+Ajpi]=p2+e2A2+(δij+iεijkσk)e[pjAi+Ajpi]\begin{align*} \left[\boldsymbol{\sigma}\cdot\left(\boldsymbol{p}+e\boldsymbol{A}\right)\right]^{2} & =\delta_{ij}\left[p_{j}p_{i}+e^{2}A_{j}A_{i}\right]+\sigma_{j}\sigma_{i}e\left[p_{j}A_{i}+A_{j}p_{i}\right]\\ & =\boldsymbol{p}^{2}+e^{2}\boldsymbol{A}^{2}+\left(\delta_{ij}+i\varepsilon_{ijk}\sigma_{k}\right)e\left[p_{j}A_{i}+A_{j}p_{i}\right] \end{align*}

Abrindo então último termo e sendo p=−i∇\boldsymbol{p}=-i\nabla, lembrando que tudo opera sobre uAu_{A}:

(δij+iεijkσk)(pjAi+Ajpi)uA=[δijpjAi+δijAjpi+iεjikσkpjAi+iεjikσkAjpi]uA=(p⋅A+A⋅p)uA+εjikσk∂j(AiuA)+εjikσkAj∂i(uA)=(p⋅A+A⋅p)uA+εjikσk∂j(Ai)uA+[εjikσkAi∂j(uA)+εjikσkAj∂i(uA)]\begin{align*} \left(\delta_{ij}+i\varepsilon_{ijk}\sigma_{k}\right)\left(p_{j}A_{i}+A_{j}p_{i}\right)u_{A} & =\left[\delta_{ij}p_{j}A_{i}+\delta_{ij}A_{j}p_{i}+i\varepsilon_{jik}\sigma_{k}p_{j}A_{i}+i\varepsilon_{jik}\sigma_{k}A_{j}p_{i}\right]u_{A}\\ & =\left(\boldsymbol{p}\cdot\boldsymbol{A}+\boldsymbol{A}\cdot\boldsymbol{p}\right)u_{A}+\varepsilon_{jik}\sigma_{k}\partial_{j}\left(A_{i}u_{A}\right)+\varepsilon_{jik}\sigma_{k}A_{j}\partial_{i}\left(u_{A}\right)\\ & =\left(\boldsymbol{p}\cdot\boldsymbol{A}+\boldsymbol{A}\cdot\boldsymbol{p}\right)u_{A}+\varepsilon_{jik}\sigma_{k}\partial_{j}\left(A_{i}\right)u_{A}+\left[\varepsilon_{jik}\sigma_{k}A_{i}\partial_{j}\left(u_{A}\right)+\varepsilon_{jik}\sigma_{k}A_{j}\partial_{i}\left(u_{A}\right)\right] \end{align*}

Importante pestar atenção nos últimos dois termos. Por exemplo, mantendo 3=13=1 fixo:

εji3Ai∂j(uA)+εji3Aj∂i(uA)=ε123A2∂1(uA)+ε123A1∂2(uA)+ε213A1∂2(uA)+ε213A2∂1(uA)=(ε123+ε213)A2∂1(uA)+(ε123+ε213)A1∂2(uA)=0\begin{align*} \varepsilon_{ji3}A_{i}\partial_{j}\left(u_{A}\right)+\varepsilon_{ji3}A_{j}\partial_{i}\left(u_{A}\right) & =\varepsilon_{123}A_{2}\partial_{1}\left(u_{A}\right)+\varepsilon_{123}A_{1}\partial_{2}\left(u_{A}\right)+\varepsilon_{213}A_{1}\partial_{2}\left(u_{A}\right)+\varepsilon_{213}A_{2}\partial_{1}\left(u_{A}\right)\\ & =\left(\varepsilon_{123}+\varepsilon_{213}\right)A_{2}\partial_{1}\left(u_{A}\right)+\left(\varepsilon_{123}+\varepsilon_{213}\right)A_{1}\partial_{2}\left(u_{A}\right)\\ & =0 \end{align*}

É análogo para outras componentes. Então ficamos com:

(δij+iεijkσk)(pjAi+Ajpi)uA=(p⋅A+A⋅p)uA+εjikσk∂j(Ai)uA=[(p⋅A+A⋅p)+σ⋅(∇×A)]uA=[(p⋅A+A⋅p)+σ⋅B]\begin{align*} \left(\delta_{ij}+i\varepsilon_{ijk}\sigma_{k}\right)\left(p_{j}A_{i}+A_{j}p_{i}\right)u_{A} & =\left(\boldsymbol{p}\cdot\boldsymbol{A}+\boldsymbol{A}\cdot\boldsymbol{p}\right)u_{A}+\varepsilon_{jik}\sigma_{k}\partial_{j}\left(A_{i}\right)u_{A}\\ & =\left[\left(\boldsymbol{p}\cdot\boldsymbol{A}+\boldsymbol{A}\cdot\boldsymbol{p}\right)+\boldsymbol{\sigma}\cdot\left(\nabla\times\boldsymbol{A}\right)\right]u_{A}\\ & =\left[\left(\boldsymbol{p}\cdot\boldsymbol{A}+\boldsymbol{A}\cdot\boldsymbol{p}\right)+\boldsymbol{\sigma}\cdot\boldsymbol{B}\right] \end{align*}

E então:

[σ⋅(p+eA)]2=p2+e2A2+(δij+iεijkσk)e[pjAi+Ajpi]=p2+e2A2+ep⋅A+eA⋅p+eσ⋅B=\begin{align*} \left[\boldsymbol{\sigma}\cdot\left(\boldsymbol{p}+e\boldsymbol{A}\right)\right]^{2} & =\boldsymbol{p}^{2}+e^{2}\boldsymbol{A}^{2}+\left(\delta_{ij}+i\varepsilon_{ijk}\sigma_{k}\right)e\left[p_{j}A_{i}+A_{j}p_{i}\right]\\ & =\boldsymbol{p}^{2}+e^{2}\boldsymbol{A}^{2}+e\boldsymbol{p}\cdot\boldsymbol{A}+e\boldsymbol{A}\cdot\boldsymbol{p}+e\boldsymbol{\sigma}\cdot\boldsymbol{B}\\ & = \end{align*}

E finalmente:

(ENR+eA0)uA=12m[σ⋅(p+eA)]2uA(ENR+eA0)uA=(p+eA)2+eσ⋅B2muA\begin{align*} \left(E_{NR}+eA^{0}\right)u_{A} & =\frac{1}{2m}\left[\boldsymbol{\sigma}\cdot\left(\boldsymbol{p}+e\boldsymbol{A}\right)\right]^{2}u_{A}\\ \left(E_{NR}+eA^{0}\right)u_{A} & =\frac{\left(\boldsymbol{p}+e\boldsymbol{A}\right)^{2}+e\boldsymbol{\sigma}\cdot\boldsymbol{B}}{2m}u_{A} \end{align*}

Ou na formulação que o exercício pediu:

ENRuA=(12m(p+eA)2+e2mσ⋅B−eA0)uAE_{NR}u_{A}=\left(\frac{1}{2m}\left(\boldsymbol{p}+e\boldsymbol{A}\right)^{2}+\frac{e}{2m}\boldsymbol{\sigma}\cdot\boldsymbol{B}-eA^{0}\right)u_{A}

Problem 5

Derive the completeness relations

∑s=1,2u(s)(p)u‾(s)(p)=p+m∑s=1,2v(s)(p)v‾(s)(p)=p−m\begin{align*} \sum_{s=1,2} & u^{\left(s\right)}\left(p\right)\overline{u}^{\left(s\right)}\left(p\right)=\cancel{p}+m\\ \sum_{s=1,2} & v^{\left(s\right)}\left(p\right)\overline{v}^{\left(s\right)}\left(p\right)=\cancel{p}-m \end{align*}

Queremos desenvolver:

∑s=1,2u(s)(p)u‾(s)(p)=p+m∑s=1,2v(s)(p)v‾(s)(p)=p−m\begin{align*} \sum_{s=1,2} & u^{\left(s\right)}\left(p\right)\overline{u}^{\left(s\right)}\left(p\right)=\cancel{p}+m\\ \sum_{s=1,2} & v^{\left(s\right)}\left(p\right)\overline{v}^{\left(s\right)}\left(p\right)=\cancel{p}-m \end{align*}

Temos então da equaçao 177 da aula 37:

χ(1)=(10),χ(2)=(01)\chi^{\left(1\right)}=\left(\begin{array}{c} 1\\ 0 \end{array}\right),\qquad\chi^{\left(2\right)}=\left(\begin{array}{c} 0\\ 1 \end{array}\right)

E da 179 e 181, sendo u(s)=u(s)u^{\left(s\right)}=u^{\left(s\right)} e u(s+2)=v(s)u^{\left(s+2\right)}=v^{\left(s\right)} :

u(1)=N(χ(1)σ⋅pE+mχ(1))=N(10p3E+mip2+p1E+m)\begin{align*} u^{\left(1\right)} & =N\left(\begin{array}{c} \chi^{\left(1\right)}\\ \frac{\boldsymbol{\sigma}\cdot\boldsymbol{p}}{E+m}\chi^{\left(1\right)} \end{array}\right) & =N\left(\begin{array}{c} \begin{array}{c} 1\\ 0 \end{array}\\ \frac{p_{3}}{E+m}\\ \frac{ip_{2}+p_{1}}{E+m} \end{array}\right) \end{align*}

Pois:

(σ⋅p)χ(1)=(σ1p1+σ2p2+σ3p3)χ(1)=(p3p1−ip2ip2+p1−p3)(10)=(p3ip2+p1)\left(\boldsymbol{\sigma}\cdot\boldsymbol{p}\right)\chi^{\left(1\right)}=\left(\sigma_{1}p_{1}+\sigma_{2}p_{2}+\sigma_{3}p_{3}\right)\chi^{\left(1\right)}=\left(\begin{array}{cc} p_{3} & p_{1}-ip_{2}\\ ip_{2}+p_{1} & -p_{3} \end{array}\right)\left(\begin{array}{c} 1\\ 0 \end{array}\right)=\left(\begin{array}{c} p_{3}\\ ip_{2}+p_{1} \end{array}\right)

De maneira análoga:

u(2)=N(χ(2)σ⋅pE+mχ(2))=N(01p1−ip2E+m−p3E+m)v(1)=N(−σ⋅pE+mχ(1)χ(1))=N(−p3E+m−ip2−p1E+m10)v(2)=N(−σ⋅pE+mχ(2)χ(2))=N(ip2−p1E+mp3E+m01)\begin{align*} u^{\left(2\right)} & =N\left(\begin{array}{c} \chi^{\left(2\right)}\\ \frac{\boldsymbol{\sigma}\cdot\boldsymbol{p}}{E+m}\chi^{\left(2\right)} \end{array}\right) & =N\left(\begin{array}{c} \begin{array}{c} 0\\ 1 \end{array}\\ \frac{p_{1}-ip_{2}}{E+m}\\ \frac{-p_{3}}{E+m} \end{array}\right)\\ v^{\left(1\right)} & =N\left(\begin{array}{c} -\frac{\boldsymbol{\sigma}\cdot\boldsymbol{p}}{E+m}\chi^{\left(1\right)}\\ \chi^{\left(1\right)} \end{array}\right) & =N\left(\begin{array}{c} \begin{array}{c} \frac{-p_{3}}{E+m}\\ \frac{-ip_{2}-p_{1}}{E+m} \end{array}\\ 1\\ 0 \end{array}\right)\\ v^{\left(2\right)} & =N\left(\begin{array}{c} -\frac{\boldsymbol{\sigma}\cdot\boldsymbol{p}}{E+m}\chi^{\left(2\right)}\\ \chi^{\left(2\right)} \end{array}\right) & =N\left(\begin{array}{c} \begin{array}{c} \frac{ip_{2}-p_{1}}{E+m}\\ \frac{p_{3}}{E+m} \end{array}\\ 0\\ 1 \end{array}\right) \end{align*}

Podemos definir no Wolfram Mathematica:

u1:=\{\{1\},\{0\},\{p3/(e+M)\},\{(I{*}p2+p1)/(e+M)\}\};~

u2:=\{\{0\},\{1\},\{(p1-I{*}p2)/(E+M)\},\{-p3/(e+M)\}\};~

v1:=\{\{-p3/(e+M)\},\{(-I{*}p2-p1)/(e+M)\},\{1\},\{0\}\};~

v2:=\{\{(I{*}p2-p1)/(e+M)\},\{p3/(e+M)\},\{0\},\{1\}\};

E agora u‾=u†γ0\overline{u}=u^{\dagger}\gamma^{0}, calculando u‾(1)\overline{u}^{\left(1\right)}:

u‾(1)=N(χ(1)(σ⋅p)†E+mχ(1))(100−1)=N(χ(1)−(σ⋅p)†E+mχ(1))=(10−p3E+mip2−p1E+m)\begin{align*} \overline{u}^{\left(1\right)} & =N\left(\begin{array}{cc} \chi^{\left(1\right)} & \frac{\left(\boldsymbol{\sigma}\cdot\boldsymbol{p}\right)^{\dagger}}{E+m}\chi^{\left(1\right)}\end{array}\right)\left(\begin{array}{cc} \boldsymbol{1} & 0\\ 0 & \boldsymbol{-1} \end{array}\right) & =N\left(\begin{array}{cc} \chi^{\left(1\right)} & -\frac{\left(\boldsymbol{\sigma}\cdot\boldsymbol{p}\right)^{\dagger}}{E+m}\chi^{\left(1\right)}\end{array}\right) & =\left(\begin{array}{cccc} 1 & 0 & \frac{-p_{3}}{E+m} & \frac{ip_{2}-p_{1}}{E+m}\end{array}\right) \end{align*}

Checando no Wolfram Mathematica:

b:=\{\{1,0,0,0\},\{0,1,0,0\},\{0,0,-1,0\},\{0,0,0,-1\}\}

lu1:=Conjugate{[}Transpose{[}u1{]}{]}.b

MatrixForm{[}v2{]}

Uma vez que está certo, declaramos o resto da mesma maneira:

u‾(2)=(u(2))†γ0,v‾(1)=(v(1))†γ0,v‾(2)=(v(2))†γ0\overline{u}^{\left(2\right)}=\left(u^{\left(2\right)}\right)^{\dagger}\gamma^{0},\quad\overline{v}^{\left(1\right)}=\left(v^{\left(1\right)}\right)^{\dagger}\gamma^{0},\quad\overline{v}^{\left(2\right)}=\left(v^{\left(2\right)}\right)^{\dagger}\gamma^{0}

O que retorna:

u‾(2)=N(01p1+ip2E+mp3E+m),v‾(1)=N(−p3E+m−p1+ip2E+m−10),v‾(1)=N(−p1−ip2E+mp3E+m−10)\overline{u}^{\left(2\right)}=N\left(\begin{array}{cccc} 0 & 1 & \frac{p_{1}+ip_{2}}{E+m} & \frac{p_{3}}{E+m}\end{array}\right),\quad\overline{v}^{\left(1\right)}=N\left(\begin{array}{cccc} \frac{-p_{3}}{E+m} & \frac{-p_{1}+ip_{2}}{E+m} & -1 & 0\end{array}\right),\quad\overline{v}^{\left(1\right)}=N\left(\begin{array}{cccc} \frac{-p_{1}-ip_{2}}{E+m} & \frac{p_{3}}{E+m} & -1 & 0\end{array}\right)

Fazendo então primeiro:

MatrixForm{[}u1.lu1+u2.lu2{]}
∑s=1,2u(s)(p)u‾(s)(p)=(E+m0−p3−p1+p20E+m−p1−ip2p3p3p1−ip2−(p12+p22+p32)E+m0p1+ip2−p30−(p12+p22+p32)E+m)\begin{align*} \sum_{s=1,2} & u^{\left(s\right)}\left(p\right)\overline{u}^{\left(s\right)}\left(p\right)=\left(\begin{array}{cccc} E+m & 0 & -p_{3} & -p_{1}+p_{2}\\ 0 & E+m & -p_{1}-ip_{2} & p_{3}\\ p_{3} & p_{1}-ip_{2} & -\frac{\left(p^{2}_{1}+p^{2}_{2}+p^{2}_{3}\right)}{E+m} & 0\\ p_{1}+ip_{2} & -p_{3} & 0 & -\frac{\left(p^{2}_{1}+p^{2}_{2}+p^{2}_{3}\right)}{E+m} \end{array}\right) \end{align*}

Uma vez que N=E+MN=\sqrt{E+M}. Além disso temos que E2=(mc2)2+(pc)2E^{2}=\left(mc^{2}\right)^{2}+\left(pc\right)^{2}, sendo c=1c=1 então p2=E2−m2p^{2}=E^{2}-m^{2}. E como:

E2−m2E+m=(E−m)(E+m)E+m=E−m\frac{E^{2}-m^{2}}{E+m}=\frac{\left(E-m\right)\left(E+m\right)}{E+m}=E-m

Temos apenas:

∑s=1,2u(s)(p)u‾(s)(p)=(E+m0−p3−p1+p20E+m−p1−ip2p3p3p1−ip2−(E−m)0p1+ip2−p30−(E−m))\sum_{s=1,2}u^{\left(s\right)}\left(p\right)\overline{u}^{\left(s\right)}\left(p\right)=\left(\begin{array}{cccc} E+m & 0 & -p_{3} & -p_{1}+p_{2}\\ 0 & E+m & -p_{1}-ip_{2} & p_{3}\\ p_{3} & p_{1}-ip_{2} & -\left(E-m\right) & 0\\ p_{1}+ip_{2} & -p_{3} & 0 & -\left(E-m\right) \end{array}\right)

Agora lembrando que:

γ0=(1000010000−10000−1)γ1=(000100100−100−1000)γ1=(000i00i00i00−i000)γ1=(00100001−10000100)\begin{align*} \gamma^{0}=\left(\begin{array}{cccc} 1 & 0 & 0 & 0\\ 0 & 1 & 0 & 0\\ 0 & 0 & -1 & 0\\ 0 & 0 & 0 & -1 \end{array}\right) & \gamma^{1}=\left(\begin{array}{cccc} 0 & 0 & 0 & 1\\ 0 & 0 & 1 & 0\\ 0 & -1 & 0 & 0\\ -1 & 0 & 0 & 0 \end{array}\right)\\ \gamma^{1}=\left(\begin{array}{cccc} 0 & 0 & 0 & i\\ 0 & 0 & i & 0\\ 0 & i & 0 & 0\\ -i & 0 & 0 & 0 \end{array}\right) & \gamma^{1}=\left(\begin{array}{cccc} 0 & 0 & 1 & 0\\ 0 & 0 & 0 & 1\\ -1 & 0 & 0 & 0\\ 0 & 1 & 0 & 0 \end{array}\right) \end{align*}

Neste ponto, uma observação é importante de ser eita. Até aqui estávamos usando pip_{i} apenas para indicar a componente ii do momento como estamos acostumados a trabalhar com vetores. Mas utilizando a notação de variante e covarianete, estavámos lidando com pμ=(E,p)p^{\mu}=\left(E,\boldsymbol{p}\right), é necessário mudar o sinal das componentes espaciais do momento pip^{i} para obtermos as componentes pi=(E,−p)p_{i}=\left(E,-\boldsymbol{p}\right) a fim de calcularmos γμpμ=p\gamma^{\mu}p_{\mu}=\cancel{p}. Desta forma podemos ver que é exatamente:

∑s=1,2u(s)(p)u‾(s)(p)=γμpμ+Im\sum_{s=1,2}u^{\left(s\right)}\left(p\right)\overline{u}^{\left(s\right)}\left(p\right)=\gamma^{\mu}p_{\mu}+\mathbb{I}m

Como queríamos demonstrar. Fazendo um processo análogo:

MatrixForm{[}v1.lv1+v2.lv2{]}

E já sabendo que (p12+p22+p32)/(E+m)=E−m\left(p^{2}_{1}+p^{2}_{2}+p^{2}_{3}\right)/(E+m)=E-m

∑s=1,2v(s)(p)v‾(s)(p)=(E−m0−p3−p1+p20E−m−p1−ip2p3p3p1−ip2−E−m0p1+ip2−p30−E−m)\begin{align*} \sum_{s=1,2} & v^{\left(s\right)}\left(p\right)\overline{v}^{\left(s\right)}\left(p\right)=\left(\begin{array}{cccc} E-m & 0 & -p_{3} & -p_{1}+p_{2}\\ 0 & E-m & -p_{1}-ip_{2} & p_{3}\\ p_{3} & p_{1}-ip_{2} & -E-m & 0\\ p_{1}+ip_{2} & -p_{3} & 0 & -E-m \end{array}\right) \end{align*}

Que de maneira análoga, nos leva a :

∑s=1,2v(s)(p)v‾(s)(p)=γμpμ−Im\sum_{s=1,2}v^{\left(s\right)}\left(p\right)\overline{v}^{\left(s\right)}\left(p\right)=\gamma^{\mu}p_{\mu}-\mathbb{I}m
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