Equação de Dirac: mudanças entre as edições
m (→Referências) |
Sem resumo de edição |
||
(2 revisões intermediárias pelo mesmo usuário não estão sendo mostradas) | |||
Linha 169: | Linha 169: | ||
=Discretização= | =Discretização= | ||
A equação de Dirac 1D pode ser escrita, na forma matricial, como: | |||
<center> | |||
<math> | |||
i \partial_t \mathbf{\Phi} = [-i\sigma_1\partial_x + \sigma_3] \mathbf{\Phi} + [V(t,x)I_2] \mathbf{\Phi} | |||
</math> | |||
</center> | |||
Onde <math>\mathbf{\Phi} = (\phi_1, \phi_4)^T</math> e <math>I_2</math> é matriz identidade de dimensão 2.\\ | |||
Para discretizar uma equação diferencial parcial como essa, é necessário discretizar o espaço e o tempo. Convenciona-se <math>\Delta t</math> como um passo finito de tempo e <math>h</math> como um passo finito no espaço, de tal forma que <math>x_j = x_0 + jh, t_n = t_0 + n\Delta t</math>, onde <math>j,n</math> são números inteiros. Define-se a notação <math>\mathbf{\Phi}(t_n, x_j) = \mathbf{\Phi}_j ^n</math> e também <math>V(t_n,x_n) = V^n _j</math>. Discretiza-se as derivadas parciais explicitamente com uma expansão em série de taylor da própria função: | |||
<center> | |||
<math> | |||
\mathbf{\Phi^{n+1} _j} = \mathbf{\Phi^n _j} + \partial_t \mathbf{\Phi^n _j} \Delta t + \mathcal{O}(\Delta t ^2) | |||
</math> | |||
</center> | |||
Considerando uma derivada discretizada <math>\delta_t \approx \partial_t</math> e truncando na primeira ordem: | |||
<center> | |||
<math> | |||
\delta_t\mathbf{\Phi^n _j} = \frac{\mathbf{\Phi^{n+1} _j} - \mathbf{\Phi^n _j}}{\Delta t} | |||
</math> | |||
</center> | |||
O processo é completamente análogo para a derivada espacial, porém para facilitar a aplicação do método mantém-se o espaço centrado em <math>x_j</math>, em outras palavras faz-se uma expansão em torno de <math>x_{j-1}</math>, obtendo: | |||
<center> | |||
<math> | |||
\delta_x\mathbf{\Phi^n _j} = \frac{\mathbf{\Phi^n _{j+1}} - \mathbf{\Phi^n _{j-1}}}{2h} | |||
</math> | |||
</center> | |||
Com isso, obtém-se uma equação para um método explícito no tempo da equação de Dirac 1D. | |||
<center> | |||
<math> | |||
i \delta_t \mathbf{\Phi^n _j} = [-i\sigma_1\delta_x + \sigma_3] \mathbf{\Phi^n _j} + [V^n _jI_2] \mathbf{\Phi^n _j} | |||
</math> | |||
</center> | |||
<center> | |||
<math> | |||
i \frac{\mathbf{\Phi^{n+1} _j} - \mathbf{\Phi^n _j}}{\Delta t} = [-i\sigma_1\delta_x + \sigma_3] \mathbf{\Phi^n _j} + [V^n _jI_2] \mathbf{\Phi^n _j} | |||
</math> | |||
</center> | |||
<center> | |||
<math> | |||
i \mathbf{\Phi^{n+1} _j} = \mathbf{\Phi^n _j} + \Delta t[-i\sigma_1\delta_x + \sigma_3] \mathbf{\Phi^n _j} + \Delta t[V^n _jI_2] \mathbf{\Phi^n _j} | |||
</math> | |||
</center> | |||
Pode-se também desenvolver um método implícito no tempo fazendo a expansão de <math>\mathbf{\Phi^{n-1} _j}</math> em torno de <math>t_n</math>, obtendo: | |||
<center> | |||
<math> | |||
\delta_t\mathbf{\Phi^n _j} = \frac{\mathbf{\Phi^{n} _j} - \mathbf{\Phi^{n-1} _j}}{\Delta t} | |||
</math> | |||
</center> | |||
Ao aplicar esta aproximação na equação discretizada basta dar um passo a frente em todos os elementos, obtendo um método implícito no tempo, já que há dependência com <math>t_{n+1}</math>. | |||
<center> | |||
<math> | |||
i \mathbf{\Phi^{n+1} _j} = \mathbf{\Phi^n _j} + \Delta t[-i\sigma_1\delta_x + \sigma_3] \mathbf{\Phi^{n+1} _j} + \Delta t[V^{n+1} _jI_2] \mathbf{\Phi^{n+1} _j} | |||
</math> | |||
</center> | |||
=Método de Crank-Nicholson= | |||
O método de Crank-Nicholson(CN) consiste em uma média entre um método explícito e outro implícito no espaço. Utilizará-se a notação <math>\mathbf{\Phi^{n+1/2} _j}</math> para representar justamente a média entre ambos os métodos, ou seja: | |||
<center> | |||
<math> | |||
\mathbf{\Phi^{n+1/2} _j} = \frac{\mathbf{\Phi^{n+1} _j} + \mathbf{\Phi^n _j}}{2} | |||
</math> | |||
</center> | |||
Define-se a notação: | |||
<center> | |||
<math> | |||
V^{n+1/2} _j\mathbf{\Phi^{n+1/2} _j} = \frac{ V^{n+1} _j\mathbf{\Phi^{n+1} _j} + V^{n} _j\mathbf{\Phi^n _j}}{2} | |||
</math> | |||
</center> | |||
Dessa maneira, enuncia-se o método CN para a equação de Dirac 1D como: | |||
<center> | |||
<math> | |||
i \delta_t \mathbf{\Phi^n _j} = [-i\sigma_1\delta_x + \sigma_3] \mathbf{\Phi^{n+1/2} _j} + [V^{n+1/2} _jI_2] \mathbf{\Phi^{n+1/2} _j} | |||
</math> | |||
</center> | |||
Onde <math>\delta_t, \delta_x </math> são as discretizações explícitas das derivadas. | |||
Para que seja possível aplicar e estudar o método é necessário passar da notação matricial para escalar: | |||
<center> | |||
<math> | |||
i \delta_t \mathbf{\Phi^n _j} = -i\sigma_1\delta_x(\frac{\mathbf{\Phi^{n+1} _j} + \mathbf{\Phi^n _j}}{2}) + \sigma_3(\frac{\mathbf{\Phi^{n+1} _j} + \mathbf{\Phi^n _j}}{2}) + I_2 (\frac{V^{n+1} _j \mathbf{\Phi^{n+1} _j} + V^n _j\mathbf{\Phi^n _j}}{2}) | |||
</math> | |||
</center> | |||
<center> | |||
<math> | |||
i \frac{\mathbf{\Phi^{n+1} _j} - \mathbf{\Phi^n _j}}{\Delta t} = -\frac{i}{2}\sigma_1[\frac{\mathbf{\Phi^{n+1} _{j+1}} - \mathbf{\Phi^{n+1} _{j-1}}}{2h} + \frac{\mathbf{\Phi^n _{j+1}} - \mathbf{\Phi^n _{j-1}}}{2h} ] + \sigma_3(\frac{\mathbf{\Phi^{n+1} _j} + \mathbf{\Phi^n _j}}{2}) + I_2 (\frac{V^{n+1} _j \mathbf{\Phi^{n+1} _j} + V^n _j\mathbf{\Phi^n _j}}{2}) | |||
</math> | |||
</center> | |||
Isolando cada tempo em um lado da igualdade: | |||
<center> | |||
<math> | |||
[I_2 + \frac{i}{2}\Delta t\sigma_3 + \frac{i}{2}\Delta t V^{n+1} _j I_2]\mathbf{\Phi^{n+1} _j} + \frac{\Delta t}{4h}\sigma_1[\mathbf{\Phi^{n+1} _{j+1}} - \mathbf{\Phi^{n+1} _{j-1}}] = </math> | |||
<math> | |||
[I_2 - \frac{i}{2}\Delta t\sigma_3 - \frac{i}{2}\Delta t V^{n} _j I_2]\mathbf{\Phi^{n} _j} - \frac{\Delta t}{4h}\sigma_1[\mathbf{\Phi^{n} _{j+1}} - \mathbf{\Phi^{n} _{j-1}}] | |||
</math> | |||
</center> | |||
Abrindo as matrizes <math>\sigma_1, \sigma_3</math> e <math>I_2</math> , e operando-as sobre o vetor <math>\mathbf{\Phi}</math> na equação tem-se: | |||
<center> | |||
<math> | |||
[\begin{bmatrix} | |||
1 & 0 \\ | |||
0 & 1 \\ | |||
\end{bmatrix} | |||
+\frac{i\Delta t}{2} | |||
\begin{bmatrix} | |||
1 & 0 \\ | |||
0 & -1 \\ | |||
\end{bmatrix} | |||
+ \frac{i}{2}\Delta t V^{n+1} _j | |||
\begin{bmatrix} | |||
1 & 0 \\ | |||
0 & 1 \\ | |||
\end{bmatrix}] | |||
\begin{bmatrix} | |||
\psi^{n+1} _{1,j} \\ | |||
\psi^{n+1} _{4,j} \\ | |||
\end{bmatrix} | |||
+\frac{\Delta t}{4h} | |||
\begin{bmatrix} | |||
0 & 1 \\ | |||
1 & 0 \\ | |||
\end{bmatrix} | |||
\begin{bmatrix} | |||
\psi^{n+1} _{1,j+1} - \psi^{n+1} _{1,j-1} \\ | |||
\psi^{n+1} _{4,j+1} - \psi^{n+1} _{4,j-1} \\ | |||
\end{bmatrix} = </math> | |||
<math> | |||
[\begin{bmatrix} | |||
1 & 0 \\ | |||
0 & 1 \\ | |||
\end{bmatrix} | |||
-\frac{i\Delta t}{2} | |||
\begin{bmatrix} | |||
1 & 0 \\ | |||
0 & -1 \\ | |||
\end{bmatrix} | |||
-\frac{i}{2}\Delta t V^{n} _j | |||
\begin{bmatrix} | |||
1 & 0 \\ | |||
0 & 1 \\ | |||
\end{bmatrix}] | |||
\begin{bmatrix} | |||
\psi^{n} _{1,j} \\ | |||
\psi^{n} _{4,j} \\ | |||
\end{bmatrix} | |||
-\frac{\Delta t}{4h} | |||
\begin{bmatrix} | |||
0 & 1 \\ | |||
1 & 0 \\ | |||
\end{bmatrix} | |||
\begin{bmatrix} | |||
\psi^{n} _{1,j+1} - \psi^{n} _{1,j-1} \\ | |||
\psi^{n} _{4,j+1} - \psi^{n} _{4,j-1} \\ | |||
\end{bmatrix} | |||
</math> | |||
</center> | |||
Pode-se realizar as operações matriciais e escrever duas equações escalares. Para facilitar notação utilizará-se <math>f^n _j = \psi^n _{1,j}</math> e <math>g^n _j = \psi^n _{4,j}</math>: | |||
<center> | |||
<math> | |||
\begin{cases} | |||
[1 + \frac{i \Delta t}{2}(V^{n+1} _j + 1)]f^{n+1} _j + \frac{\Delta t}{4h}(g^{n+1}_{j+1} - g^{n+1}_{j-1}) = [1 - \frac{i \Delta t}{2}(V^{n} _j + 1)]f^{n} _j - \frac{\Delta t}{4h}(g^{n}_{j+1} - g^{n}_{j-1}) \\ | |||
[1 + \frac{i \Delta t}{2}(V^{n+1} _j - 1)]g^{n+1} _j + \frac{\Delta t}{4h}(f^{n+1}_{j+1} - f^{n+1}_{j-1} ) = [1 - \frac{i \Delta t}{2}(V^{n} _j - 1)]g^{n} _j - \frac{\Delta t}{4h}(f^{n}_{j+1} - f^{n}_{j-1}) \\ | |||
\end{cases} | |||
</math> | |||
</center> | |||
Tem-se então um número <math>n</math> de equações onde <math>n</math> é o número de elementos do espaço discretizado. Portanto o primeiro termo das duas equações gera uma matriz diagonal pois multiplica os termos espaciais dependentes de <math>x_j</math>, já o segundo termo gera uma matriz tridiagonal com diagonal principal nula. Nota-se que os primeiros termos dos dois lados da igualdade são um o conjugado do outro, define-se, portanto, <math>\alpha^n = 1 + \frac{i \Delta t}{2}(V^{n+1} _j + 1)</math> e <math>\beta = 1 + \frac{i \Delta t}{2}(V^{n+1} _j - 1)</math>. | |||
<center> | |||
<math> | |||
\begin{cases} | |||
\alpha^{n+1}f^{n+1} _j + \dfrac{\Delta t}{4h}(g^{n+1}_{j+1} - g^{n+1}_{j-1}) = | |||
\alpha^{n^*}f^{n} _j - \dfrac{\Delta t}{4h}(g^{n}_{j+1} - g^{n}_{j-1} ) \\ | |||
\beta^{n+1}g^{n+1} _j + \dfrac{\Delta t}{4h}(f^{n+1}_{j+1} - f^{n+1}_{j-1} ) = | |||
\beta^{n^*}g^{n} _j - \dfrac{\Delta t}{4h}(f^{n}_{j+1} - f^{n}_{j-1}) \\ | |||
\end{cases} | |||
</math> | |||
</center> | |||
Considerando que o potencial V é só função da posição, escreve-se o método como: | |||
<center> | |||
<math> | |||
\begin{cases} | |||
Af^{n+1} + Bg^{n+1} = A^*f^n - Bg^n \\ | |||
Cg^{n+1} + Bf^{n+1} = C^*g^n - Bf^n \\ | |||
\end{cases} | |||
</math> | |||
</center> | |||
Onde: | |||
<center> | |||
<math> | |||
A = \begin{bmatrix} | |||
\alpha & 0 & 0 & \cdots & 0\\ | |||
0 & \alpha & 0 & \cdots & 0 \\ | |||
0 & 0 & \alpha & \cdots & 0 \\ | |||
\vdots & \vdots & \vdots & \ddots & \vdots \\ | |||
0 & \cdots & \cdots & \cdots & \alpha \\ | |||
\end{bmatrix}; | |||
B = \begin{bmatrix} | |||
0 & \frac{\Delta t}{4h} & 0 & \cdots & 0\\ | |||
-\frac{\Delta t}{4h} & 0 & \frac{\Delta t}{4h} & \cdots & 0 \\ | |||
0 & -\frac{\Delta t}{4h} & 0 & \cdots & 0 \\ | |||
\vdots & \vdots & \ddots & \ddots & \vdots \\ | |||
0 & \cdots & \cdots & -\frac{\Delta t}{4h} & 0 \\ | |||
\end{bmatrix}; | |||
C = \begin{bmatrix} | |||
\beta & 0 & 0 & \cdots & 0\\ | |||
0 & \beta & 0 & \cdots & 0 \\ | |||
0 & 0 & \beta & \cdots & 0 \\ | |||
\vdots & \vdots & \vdots & \ddots & \vdots \\ | |||
0 & \cdots & \cdots & \cdots & \beta \\ | |||
\end{bmatrix}; | |||
</math> | |||
</center> | |||
Por fim, pode-se escrever o método resolvendo o sistema: | |||
<center> | |||
<math> | |||
\begin{cases} | |||
f^{n+1} = F^{-1}D f^n -F^{-1}E g^n \\ | |||
g^{n+1} = J^{-1}G g^n - J^{-1}H f^n \\ | |||
\end{cases} | |||
</math> | |||
</center> | |||
Onde <math>D = (B^{-1}A^* + C^{-1}B)$, $E = (I - C^{-1}C^*)$, $F = (B^{-1}A - C^{-1}B)$, $G = (A^{-1}A^* + I)$, $H = (A^{-1}B + B^{-1}C^*)$ e $J= (A^{-1}B - B^{-1}C)</math>. | |||
Com isso, condições iniciais e de contorno bem estabelecidas, já é possível aplicar o método, dado que todas essas matrizes dependem somente dos parâmetros do sistema. | |||
=Estabilidade Crank-Nicholson= | |||
Utilizará-se o método de Von-Neumann para analisar a estabilidade do método para a equação de Dirac unidimensional, para supõe-se que a função $\mathbf{\Phi^{n} _j}$ pode ser dada pela série de fourier: | |||
<center> | |||
<math> | |||
\mathbf{\Phi^{n} _j} = \sum_{k=0}^{\inf} \mathbf{A}^{n}e^{ikqjh} | |||
</math> | |||
</center> | |||
Devido à indepêndencia linear de cada termo do somatório, ao substituir na equação do método haverá uma equação para cada ente do somatório. Se o módulo da razão $\frac{\mathbf{A}^{n+1}}{\mathbf{A}^{n}} \le 1$ então é pode-se dizer que o método estável, já que dessa forma garante-se uma não divergência. | |||
Aplica-se um termo geral da série de índice $k$ no método CN para a equação de dirac 1D. | |||
<center> | |||
<math> | |||
[I_2 + \frac{i}{2}\Delta t\sigma_3 + \frac{i}{2}\Delta t V^{n+1} _j I_2]\mathbf{A}^{n+1}e^{ikqjh} + \frac{\Delta t}{4h}\sigma_1\mathbf{A}^{n+1}[e^{ikq(j+1)h} - e^{ikq(j-1)h}] =</math> | |||
<math> | |||
[I_2 - \frac{i}{2}\Delta t\sigma_3 - \frac{i}{2}\Delta t V^{n} _jI_2]\mathbf{A}^{n}e^{ikqjh} - \frac{\Delta t}{4h}\sigma_1\mathbf{A}^{n}[e^{ikq(j+1)h} - e^{ikq(j-1)h}] | |||
</math> | |||
</center> | |||
Divide-se tudo por <math>e^{ikqjh}</math> e isola-se <math>\mathbf{A}</math>: | |||
<center> | |||
<math> | |||
[I_2 + \frac{i}{2}\Delta t\sigma_3 + \frac{i}{2}\Delta t V^{n+1} _j I_2 + \frac{\Delta t}{4h}\sigma_1(e^{ikqh} - e^{-ikqh}) ]\mathbf{A}^{n+1}=</math> | |||
<math> | |||
[I_2 - \frac{i}{2}\Delta t\sigma_3 - \frac{i}{2}\Delta t V^{n} _jI_2 - \frac{\Delta t}{4h}\sigma_1(e^{ikqh} - e^{-ikqh})]\mathbf{A}^{n} | |||
</math> | |||
<math> | |||
[I_2 + \frac{i}{2}\Delta t\sigma_3 + \frac{i}{2}\Delta t V^{n+1} _j I_2 + \frac{\Delta t}{4h}\sigma_1(\frac{\sin (kqh)}{2i}) ]\mathbf{A}^{n+1}= | |||
</math> | |||
<math> | |||
[I_2 - \frac{i}{2}\Delta t\sigma_3 - \frac{i}{2}\Delta t V^{n} _jI_2 - \frac{\Delta t}{4h}\sigma_1(\frac{\sin (kqh)}{2i})]\mathbf{A}^{n} | |||
</math> | |||
<math> | |||
[I_2 + \frac{i}{2}\Delta t\sigma_3 + \frac{i}{2}\Delta t V^{n+1} _j I_2 - \frac{i\Delta t}{4h}\sigma_1(\frac{\sin (kqh)}{2}) ]\mathbf{A}^{n+1}= | |||
</math> | |||
<math> | |||
[I_2 - \frac{i}{2}\Delta t\sigma_3 - \frac{i}{2}\Delta t V^{n} _jI_2 + \frac{i\Delta t}{4h}\sigma_1(\frac{\sin (kqh)}{2})]\mathbf{A}^{n} | |||
</math> | |||
</center> | |||
Nota-se que os termos que multiplicam o fator <math>\mathbf{A}</math> são o conjugado um do outro, define-se <math>z = I_2 + \frac{i}{2}\Delta t\sigma_3 + \frac{i}{2}\Delta t V^{n+1} _j I_2 - \frac{i\Delta t}{4h}\sigma_1(\frac{\sin (kqh)}{2})</math>, dessa maneira: | |||
<center> | |||
<math> | |||
\frac{\mathbf{A}^{n+1}}{\mathbf{A}^{n}} = \frac{z^*}{z} | |||
</math> | |||
<math> | |||
|\frac{\mathbf{A}^{n+1}}{\mathbf{A}^{n}}| = |\frac{z^*}{z}| = \frac{|z^*|}{|z|} = 1 | |||
</math> | |||
</center> | |||
Onde <math>z</math> é sempre diferente de zero, dado que a parte real é dada por uma matriz identidade constante. | |||
Mostra-se, portanto, que a razão entre os coeficientes da série de fourier nunca divergem, ou seja, o método é incondicionalmente estável. | |||
Edição das 20h27min de 3 de maio de 2024
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 , 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:
onde, como anteriormente, os autovalores de 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:
Assim, o Hamiltoniano é modificado para representar a medida da energia relativística total.
Diferentemente da equação de Schrödinger, aqui não representa apenas uma função de onda, mas sim um conjunto de quatro delas. Usando a notação
,
as componentes de representam as funções de onda associadas ao elétron e ao pósitron: () representa a função de onda do elétron com spin up (down), e () representa a função de onda do pósitron com spin up (down). O objeto é chamado de spinor.
Dedução da equação de Dirac em duas dimensões
Consideraremos neste trabalho a equação de Dirac em duas dimensões, e . 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 completo
Consideremos uma partícula sob ação de um potencial (onde ), que afeta a energia potencial da partícula, e de um potencial "escalar" , que afeta a massa de repouso da mesma. Seguindo uma das propostas possíveis para o Hamiltoniano, temos
onde ; e são matrizes 4x4 adimensionais e é o vetor momento linear da partícula.
Pode-se mostrar que e devem satisfazer certos vínculos, limitando as escolhas possíveis para essas matrizes. A escolha mais simples e usualmente adotada consiste em tomar
Sendo , podemos escrever o produto escalar como
Portanto, em notação matricial o Hamiltoniano pode ser escrito como
Unidades naturais e redução para duas dimensões
A fim de simplificar o formalismo, adotamos as chamadas "unidades naturais", onde . Note que isso equivale a reescalar as quantidades físicas do problema por um fator adequado. Ao fazer , também assumimos que a partícula está no limite relativístico.
Além disso, queremos estudar o problema em duas dimensões. Observamos que ; logo, . Portanto, temos o Hamiltoniano simplificado
Forma explícita final
Retornando ao problema original, queremos resolver
Novamente utilizando a notação matricial, obtemos
Realizando a multiplicação matricial, pode-se ver que se obtém dois sistemas acoplados: com e com . Escolhendo o sistema de com :
Simplificando e isolando a derivada temporal, obtemos por fim
Discretização
A equação de Dirac 1D pode ser escrita, na forma matricial, como:
Onde e é matriz identidade de dimensão 2.\\ Para discretizar uma equação diferencial parcial como essa, é necessário discretizar o espaço e o tempo. Convenciona-se como um passo finito de tempo e como um passo finito no espaço, de tal forma que , onde são números inteiros. Define-se a notação e também . Discretiza-se as derivadas parciais explicitamente com uma expansão em série de taylor da própria função:
Considerando uma derivada discretizada e truncando na primeira ordem:
O processo é completamente análogo para a derivada espacial, porém para facilitar a aplicação do método mantém-se o espaço centrado em , em outras palavras faz-se uma expansão em torno de , obtendo:
Com isso, obtém-se uma equação para um método explícito no tempo da equação de Dirac 1D.
Pode-se também desenvolver um método implícito no tempo fazendo a expansão de em torno de , obtendo:
Ao aplicar esta aproximação na equação discretizada basta dar um passo a frente em todos os elementos, obtendo um método implícito no tempo, já que há dependência com .
Método de Crank-Nicholson
O método de Crank-Nicholson(CN) consiste em uma média entre um método explícito e outro implícito no espaço. Utilizará-se a notação para representar justamente a média entre ambos os métodos, ou seja:
Define-se a notação:
Dessa maneira, enuncia-se o método CN para a equação de Dirac 1D como:
Onde são as discretizações explícitas das derivadas. Para que seja possível aplicar e estudar o método é necessário passar da notação matricial para escalar:
Isolando cada tempo em um lado da igualdade:
Abrindo as matrizes e , e operando-as sobre o vetor na equação tem-se:
Pode-se realizar as operações matriciais e escrever duas equações escalares. Para facilitar notação utilizará-se e :
Falhou ao verificar gramática (Erro de conversão. Servidor ("https://wikimedia.org/api/rest_") reportou: "Cannot get mml. TeX parse error: Bracket argument to \\ must be a dimension"): {\displaystyle {\begin{cases}[1+{\frac {i\Delta t}{2}}(V_{j}^{n+1}+1)]f_{j}^{n+1}+{\frac {\Delta t}{4h}}(g_{j+1}^{n+1}-g_{j-1}^{n+1})=[1-{\frac {i\Delta t}{2}}(V_{j}^{n}+1)]f_{j}^{n}-{\frac {\Delta t}{4h}}(g_{j+1}^{n}-g_{j-1}^{n})\\[1+{\frac {i\Delta t}{2}}(V_{j}^{n+1}-1)]g_{j}^{n+1}+{\frac {\Delta t}{4h}}(f_{j+1}^{n+1}-f_{j-1}^{n+1})=[1-{\frac {i\Delta t}{2}}(V_{j}^{n}-1)]g_{j}^{n}-{\frac {\Delta t}{4h}}(f_{j+1}^{n}-f_{j-1}^{n})\\\end{cases}}}
Tem-se então um número de equações onde é o número de elementos do espaço discretizado. Portanto o primeiro termo das duas equações gera uma matriz diagonal pois multiplica os termos espaciais dependentes de , já o segundo termo gera uma matriz tridiagonal com diagonal principal nula. Nota-se que os primeiros termos dos dois lados da igualdade são um o conjugado do outro, define-se, portanto, e .
Considerando que o potencial V é só função da posição, escreve-se o método como:
Onde:
Por fim, pode-se escrever o método resolvendo o sistema:
Onde .
Com isso, condições iniciais e de contorno bem estabelecidas, já é possível aplicar o método, dado que todas essas matrizes dependem somente dos parâmetros do sistema.
Estabilidade Crank-Nicholson
Utilizará-se o método de Von-Neumann para analisar a estabilidade do método para a equação de Dirac unidimensional, para supõe-se que a função $\mathbf{\Phi^{n} _j}$ pode ser dada pela série de fourier:
Devido à indepêndencia linear de cada termo do somatório, ao substituir na equação do método haverá uma equação para cada ente do somatório. Se o módulo da razão $\frac{\mathbf{A}^{n+1}}{\mathbf{A}^{n}} \le 1$ então é pode-se dizer que o método estável, já que dessa forma garante-se uma não divergência.
Aplica-se um termo geral da série de índice $k$ no método CN para a equação de dirac 1D.
Divide-se tudo por e isola-se :
Nota-se que os termos que multiplicam o fator são o conjugado um do outro, define-se , dessa maneira:
Onde é sempre diferente de zero, dado que a parte real é dada por uma matriz identidade constante. Mostra-se, portanto, que a razão entre os coeficientes da série de fourier nunca divergem, ou seja, o método é incondicionalmente estável.
Referências
- The quantum theory of the electron. Proceedings of the Royal Society of London A, v. 117, n. 778, p. 610–624, fev. 1928.
- SAKURAI, J. J. Mecânica quântica moderna. 2. ed. Porto Alegre: Bookman, 2012.
- 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.
- SOFF, G. et al. Solution of the Dirac Equation for Scalar Potentials and its Implications in Atomic Physics. Zeitschrift für Naturforschung A, v. 28, n. 9, p. 1389–1396, 1 set. 1973.
- THALLER, B. The Dirac equation. Berlin: Springer, 2010.