Equação de Dirac

Grupo: André Luis Della Valentina, Lucas dos Santos Assmann, Vinícius Bayne Müller

Introdução

A equação de Dirac descreve uma partícula relativística de spin ${\frac {1}{2}}$, como o elétron, com estrutura análoga a da equação de Schrödinger. Ela surgiu inicialmente como tentativa de explicar discrepâncias entre experimentos e teoria para o espectro do átomo de hidrogênio, possibilitando correções para o cálculo da energia do elétron em diferentes níveis (as chamadas correções de estrutura fina). Além disso, por meio dela foi possível prever a existência de antimatéria: descrevendo o elétron, ela também descreve o pósitron.

A fim de compatibilizar a Mecânica Quântica com a Relatividade Especial, a equação diferencial parcial é de primeira ordem tanto no tempo quanto no espaço (diferentemente da equação de Schrödinger, que é de segunda ordem no espaço). A equação de Dirac pode ser escrita de diversas formas; aqui, apresentamo-la explicitamente como um sistema de EDPs acopladas, mais conveniente para os propósitos do trabalho.

Assim como a equação de Schrödinger, a construção da equação de Dirac vem a partir do operador Hamiltoniano, que descreve a evolução temporal do estado quântico em questão:

$i\hbar {\frac {\partial }{\partial t}}|\psi (t)\rangle =H|\psi (t)\rangle$

onde, como anteriormente, os autovalores de $H$ correspondem aos valores possíveis de energia que o sistema pode assumir.

A mudança com relação à Mecânica Quântica não relativística acontece quando consideramos a energia relativística da partícula:

$E^{2}=p^{2}c^{2}+m^{2}c^{4}$

Assim, o Hamiltoniano é modificado para representar a medida da energia relativística total.

Dedução da equação de Dirac em duas dimensões

Consideraremos neste trabalho a equação de Dirac em duas dimensões, $x$ e $y$. A escolha dessas coordenadas se dá pela conveniência do acoplamento das EDPs: nesse caso, as quatro equações acopladas passam a ser acopladas duas a duas, facilitando o estudo do sistema.

Construção do Hamiltoniano

Consideremos uma partícula sob ação de um potencial $V({\boldsymbol {x}};t)$ (onde ${\boldsymbol {x}}=(x,y,z)^{T}$), que afeta a energia potencial da partícula, e de um potencial "escalar" $V_{sc}({\boldsymbol {x}};t)$, que afeta a massa de repouso da mesma. Seguindo uma das propostas possíveis para o Hamiltoniano, temos

$H=c{\boldsymbol {\alpha }}\cdot {\boldsymbol {p}}+\beta (mc^{2}+V_{sc})+VI_{4}$

onde ${\boldsymbol {\alpha }}=\alpha _{x}{\hat {i}}+\alpha _{y}{\hat {j}}+\alpha _{z}{\hat {k}}$; $\alpha _{i}$ e $\beta$ são matrizes 4x4 adimensionais e ${\boldsymbol {p}}$ é o vetor momento linear da partícula.

Pode-se mostrar que ${\boldsymbol {\alpha }}$ e $\beta$ devem satisfazer certos vínculos, limitando as escolhas possíveis para essas matrizes. A escolha mais simples e usualmente adotada consiste em tomar

$\beta ={\begin{pmatrix}I_{2}&0\\0&-I_{2}\end{pmatrix}}={\begin{pmatrix}1&0&0&0\\0&1&0&0\\0&0&-1&0\\0&0&0&-1\end{pmatrix}}$

$\alpha _{x}={\begin{pmatrix}0&\sigma _{x}\\\sigma _{x}&0\end{pmatrix}}={\begin{pmatrix}0&0&0&1\\0&0&1&0\\0&1&0&0\\1&0&0&0\end{pmatrix}}$

$\alpha _{y}={\begin{pmatrix}0&\sigma _{y}\\\sigma _{y}&0\end{pmatrix}}={\begin{pmatrix}0&0&0&-i\\0&0&i&0\\0&-i&0&0\\i&0&0&0\end{pmatrix}}$

$\alpha _{z}={\begin{pmatrix}0&\sigma _{z}\\\sigma _{z}&0\end{pmatrix}}={\begin{pmatrix}0&0&1&0\\0&0&0&-1\\1&0&0&0\\0&-1&0&0\end{pmatrix}}$

Sendo ${\boldsymbol {p}}=-i\hbar \nabla$, podemos escrever o produto escalar ${\boldsymbol {\alpha }}\cdot {\boldsymbol {p}}$ como

${\boldsymbol {\alpha }}\cdot {\boldsymbol {p}}=-i\hbar \left(\alpha _{x}{\frac {\partial }{\partial x}}+\alpha _{y}{\frac {\partial }{\partial y}}+\alpha _{z}{\frac {\partial }{\partial z}}\right)$

Referências

1. The quantum theory of the electron. Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character, v. 117, n. 778, p. 610–624, fev. 1928.
2. SAKURAI, J. J. Mecânica quântica moderna. 2. ed. Porto Alegre: Bookman, 2012.
3. BAO, W. et al. Numerical Methods and Comparison for the Dirac Equation in the Nonrelativistic Limit Regime. Journal of Scientific Computing, v. 71, n. 3, p. 1094–1134, jun. 2017.