Equação de Klein-Gordon

INTRODUÇÃO

A equação de Klein-Gordon é uma das equações fundamentais na teoria quântica relativística. Ela descreve partículas escalares (partículas sem spin, como os mésons, em seu modelo básico) e é uma extensão relativística da equação de Schrödinger, incorporando a relação de energia relativística de Einstein $E=p^2{c}^2+m^2{c}^4$. A equação é nomeada em homenagem a Oskar Klein e Walter Gordon, que a formularam independentemente. De maneira geral, a equação pode ser escrita como:

$(◻+\frac{m^2{c}^2}{{\hslash }^2})\psi (x,t)=0$

onde $(◻=\frac{{\partial }^2}{c^2\partial t^2}-{\mathrm{\nabla }}^2)$ é chamado operador de d'Alambert.

Abrindo a equação, é obtido:

$\frac{{\partial }^2\psi }{\partial t^2}=c^2{\mathrm{\nabla }}^2\psi -\frac{m^2{c}^4}{{\hslash }^2}\psi$

$\frac{{\partial }^2\psi }{\partial t^2}=c^2\frac{{\partial }^2\psi }{\partial x^2}-\frac{m^2{c}^4}{{\hslash }^2}\psi$ (em uma dimensão)

MÉTODO DAS DIFERENÇAS FINITAS

O método das diferenças finitas é uma técnica numérica amplamente utilizada para resolver EDPs. Ele envolve a discretização das variáveis contínuas (geralmente no tempo ou no espaço), transformando as equações diferenciais em sistemas algébricos que podem ser resolvidos numericamente. Os primeiros passos para utilizar o método é fazer a discretização no tempo e no espaço. Para uma equação no tempo você discretiza o tempo em intervalos $\mathrm{\Delta }t$ criando uma sequência de pontos $t_n=n\mathrm{\Delta }t$. Para uma equação no espaço você discretiza o espaço em intervalos $\mathrm{\Delta }x$ criando uma sequência de pontos $x_i=i\mathrm{\Delta }x$. Depois de discretizar o espaço e o tempo, as derivadas contínuas são aproximadas por diferenças finitas. Isso envolve substituir as derivadas por aproximações baseadas nos valores de uma função nos pontos discretos:

$\frac{\partial u}{\partial t}\approx \frac{u_i^{n+1}-u_i^n}{\mathrm{\Delta }t}$ e $\frac{{\partial }^{2}u}{\partial t^2}\approx \frac{u_i^{n+1}-2u_i^n+u_i^{n-1}}{(\mathrm{\Delta }t{)}^2}$ para o tempo.

$\frac{\partial u}{\partial x}\approx \frac{u_{i+1}^n-u_i^n}{\mathrm{\Delta }x}e\frac{{\partial }^{2}u}{\partial x^2}\approx \frac{u_{i+1}^n-2u_i^n+u_{i-1}^n}{(\mathrm{\Delta }x{)}^2}$ para o espaço.

Na equação de Klein-Gordon, escrevemos desta o método das diferenças finitas:

$\frac{{\partial }^2\psi (x,t)}{\partial t^2}\approx \frac{\psi (x,t+\mathrm{\Delta }t)-2\psi (x,t)+\psi (x,t-\mathrm{\Delta }t)}{(\mathrm{\Delta }t{)}^2}$

$\frac{{\partial }^2\psi (x,t)}{\partial x^2}\approx \frac{\psi (x+\mathrm{\Delta }x,t)-2\psi (x,t)+\psi (x-\mathrm{\Delta }x,t)}{(\mathrm{\Delta }x{)}^2}$

ou seja:

$\frac{\psi (x,t+\mathrm{\Delta }t)-2\psi (x,t)+\psi (x,t-\mathrm{\Delta }t)}{(\mathrm{\Delta }t{)}^2}=c^2\frac{\psi (x+\mathrm{\Delta }x,t)-2\psi (x,t)+\psi (x,t-\mathrm{\Delta }t)}{(\mathrm{\Delta }x{)}^2}-\frac{m^2{c}^4}{{\hslash }^2}\psi$

isso nos leva a equação final:

$\psi (x,t+\mathrm{\Delta }t)-2\psi (x,t)+\psi (x,t-\mathrm{\Delta }t)=\frac{c^2\mathrm{\Delta }{t}^2}{\mathrm{\Delta }{x}^2}\psi (x+\mathrm{\Delta }x,t)-2\psi (x,t)+\psi (x,t-\mathrm{\Delta },t)-\frac{m^2{c}^4\mathrm{\Delta }{t}^2}{{\hslash }^2}\psi$

chamarei $\alpha =\frac{c\mathrm{\Delta }t}{\mathrm{\Delta }x}$ e $\beta =\frac{m{c}^2\mathrm{\Delta }t}{\hslash }$

portanto, $\psi (x,t+\mathrm{\Delta }t)-2\psi (x,t)+\psi (x,t-\mathrm{\Delta }t)={\alpha }^2\psi (x+\mathrm{\Delta }x,t)-2\psi (x,t)+\psi (x,t-\mathrm{\Delta }t)-{\beta }^2\psi$

ou, mais usualmente: ${\psi }_i^{n+1}=2{\psi }_i^n-{\psi }_i^{n-1}+{\alpha }^2({\psi }_{i+1}^n-2{\psi }_i^n+{\psi }_{i-1}^n)-{\beta }^2\psi$

ESTABILIDADE ta errado, estou testando outras coisas

Para analisar a estabilidade do método utilizaremos os Modos de Furrier.

${\psi }_i^{n+1}=2{\psi }_i^n-{\psi }_i^{n-1}+{\alpha }^2({\psi }_{i+1}^n-2{\psi }_i^n+{\psi }_{i-1}^n)-{\beta }^2{\psi }_i^n$.

sendo ${\psi }_i^n=A^n{e}^{iki\mathrm{\Delta }x}$ o Modo de Furrier.

Substituímos ${\psi }_i^n,{\psi }_{i+1}^n$, e ${\psi }_{i-1}^n$na equação:

${\psi }_{i+1}^n=A^n{e}^{ik(i+1)\mathrm{\Delta }x}=A^n{e}^{iki\mathrm{\Delta }x}{e}^{ik\mathrm{\Delta }x}$

${\psi }_{i-1}^n=A^n{e}^{ik(i-1)\mathrm{\Delta }x}=A^n{e}^{iki\mathrm{\Delta }x}{e}^{-ik\mathrm{\Delta }x}$

Usamos a identidade $e^{ik\mathrm{\Delta }x}+e^{-ik\mathrm{\Delta }x}=2\cos (k\mathrm{\Delta }x)$:

$A^{n+1}=2A^n-A^{n-1}+{\alpha }^2[2\cos (k\mathrm{\Delta }x)-2]A^n-{\beta }^2{A}^n$

Fatoramos:

$A^{n+1}=(2-2{\alpha }^2-{\beta }^2+2{\alpha }^2\cos (k\mathrm{\Delta }x))A^n-A^{n-1}.$

A relação de recorrência é:

$A^{n+1}-\lambda A^n+A^{n-1}=0$

onde

$\lambda =2-2{\alpha }^2-{\beta }^2+2{\alpha }^2\cos (k\mathrm{\Delta }x)$.

Definimos $\mu =|\frac{A^{n+1}}{A^n}|$ como sendo o fator de amplificação. Assim, a equação fica $\mu A^n=\lambda A^n-A^n\frac{A^{n-1}}{A^n}$

Dividindo tudo por $A^n$ : $\mu =\lambda -\frac{1}{\mu }$

Portanto, a equação característica associada é:

${\mu }^2-\lambda \mu +1=0$

onde $\mu$ são as raízes que representam o fator de amplificação $\mu =|\frac{A^{n+1}}{A^n}|$.

Para que o método seja estável, as raízes $\mu$ devem satisfazer $\Vert \mu |\le 1$. Isso é garantido se o discriminante da equação característica satisfizer:

${\lambda }^2-4\le 0$.

Substituímos $\lambda$:

${(2-2{\alpha }^2-{\beta }^2+2{\alpha }^2\cos (k\mathrm{\Delta }x))}^2-4\le 0$.

O caso crítico ocorre para o maior valor de $cos(k\mathrm{\Delta }x)$, que é $cos(k\mathrm{\Delta }x)=1$, e o menor valor, $cos(k\mathrm{\Delta }x)=-1$:

$cos(k\mathrm{\Delta }x)=1$:

$\lambda =2-2{\alpha }^2-{\beta }^2+2{\alpha }^2$.

Isso simplifica para:

$\lambda =2-{\beta }^2$

Para estabilidade:

$(2-{\beta }^2{)}^2-4\le 0$.

$\beta \ge 0$.

Para $cos(k\mathrm{\Delta }x)=-1$:

$\lambda =2-2{\alpha }^2-{\beta }^2-2{\alpha }^2$.

$\lambda =2-4{\alpha }^2-{\beta }^2$.

Ou seja, para que seja estável:

$(2-4{\alpha }^2-{\beta }^2{)}^2-4\le 0$.

Após expandir:

$\alpha \ge \frac{\beta }{2}$.

A condição de estabilidade combinada é:

$\alpha \ge \frac{\beta }{2}$ e $\beta \ge 0$ .

Essas desigualdades controlam os passos de tempo e espaço, garantindo a estabilidade do método.

CRITÉRIO DE ESTABILIDADE

Para analisar a estabilidade, considera-se uma solução modal da forma: ${\varphi }_j^n=e^{i(kj\mathrm{\Delta }x-\omega n\mathrm{\Delta }t)},$

onde $k$ é o número de onda e $\omega$ é a frequência angular.

Substituindo essa expressão na equação discreta e usando as propriedades das exponenciais complexas, temos: $\frac{e^{-i\omega \mathrm{\Delta }t}-2+e^{i\omega \mathrm{\Delta }t}}{\mathrm{\Delta }{t}^2}=\frac{c^2}{\mathrm{\Delta }{x}^2}(e^{ik\mathrm{\Delta }x}-2+e^{-ik\mathrm{\Delta }x})-\frac{m^2{c}^4}{{\hslash }^2}.$

${\varphi }_{j+1}^n=e^{ik\mathrm{\Delta }x}{\varphi }_j^n,{\varphi }_{j-1}^n=e^{-ik\mathrm{\Delta }x}{\varphi }_j^n,$ ${\varphi }_j^{n+1}=e^{-i\omega \mathrm{\Delta }t}{\varphi }_j^n,{\varphi }_j^{n-1}=e^{i\omega \mathrm{\Delta }t}{\varphi }_j^n.$

Substituindo essas expressões na equação discreta: $e^{-i\omega \mathrm{\Delta }t}{\varphi }_j^n-2{\varphi }_j^n+e^{i\omega \mathrm{\Delta }t}{\varphi }_j^n=\frac{c^2\mathrm{\Delta }{t}^2}{\mathrm{\Delta }{x}^2}(e^{ik\mathrm{\Delta }x}{\varphi }_j^n-2{\varphi }_j^n+e^{-ik\mathrm{\Delta }x}{\varphi }_j^n)-\frac{m^2{c}^4\mathrm{\Delta }{t}^2}{{\hslash }^2}{\varphi }_j^n.$

Dividindo todos os termos por ${\varphi }_j^n$ (assumindo ${\varphi }_j^n\ne 0$) e rearranjando, obtém-se: $\frac{e^{-i\omega \mathrm{\Delta }t}-2+e^{i\omega \mathrm{\Delta }t}}{\mathrm{\Delta }{t}^2}=\frac{c^2}{\mathrm{\Delta }{x}^2}(e^{ik\mathrm{\Delta }x}-2+e^{-ik\mathrm{\Delta }x})-\frac{m^2{c}^4}{{\hslash }^2}.$

Utilizando identidades trigonométricas para simplificar, a equação se torna: $\frac{-4{\sin }^2(\omega \mathrm{\Delta }t/2)}{\mathrm{\Delta }{t}^2}=\frac{-4c^2{\sin }^2(k\mathrm{\Delta }x/2)}{\mathrm{\Delta }{x}^2}-\frac{m^2{c}^4}{{\hslash }^2}.$

Multiplicando por $-1$, obtemos: $\frac{4{\sin }^2(\omega \mathrm{\Delta }t/2)}{\mathrm{\Delta }{t}^2}=\frac{4c^2{\sin }^2(k\mathrm{\Delta }x/2)}{\mathrm{\Delta }{x}^2}+\frac{m^2{c}^4}{{\hslash }^2}.$ Critério de Estabilidade

A condição de estabilidade exige que ${\sin }^2(\omega \mathrm{\Delta }t/2)$ seja real e menor ou igual a 1. Isso leva à restrição: $\frac{4c^2{\sin }^2(k\mathrm{\Delta }x/2)}{\mathrm{\Delta }{x}^2}+\frac{m^2{c}^4}{{\hslash }^2}\le \frac{4}{\mathrm{\Delta }{t}^2}.$

Logo, o passo temporal $\mathrm{\Delta }t$ deve satisfazer:

$\mathrm{\Delta }t\le 2\sqrt{\frac{\mathrm{\Delta }{x}^2}{4c^2{\sin }^2(k\mathrm{\Delta }x/2)}+\frac{{\hslash }^2}{m^2{c}^4}}.$

No limite de $m\to 0$ (ondas livres), a condição reduz-se ao critério CFL clássico: $\frac{c\mathrm{\Delta }t}{\mathrm{\Delta }x}\le 1$

C.C e C.I

Condições iniciais e condições de contorno são fundamentais para a resolução da equação, já que elas ditam o comportamento da função oa longo do tempo e ao longo do espaço, para plotar a evolução temporal, utilizarei as seguintes condilções iniciais e de contorno:

$\psi (x,0)=A{e}^{-\frac{x-x_0}{2{\sigma }^2}}$ que define um pulso gaussiano como condição inicial.

e $\frac{\partial \psi (x,0)}{\partial t}=0$ que define que, no instante de tempo t=0, a função não possui velocidade inicial, o que implica que o pulso está parado inicialmente e sua evolução se deve pela propagação de flutuações espaciais.

Nesta condição, A é a altura do pulso, $x_0$ é a posição central do pulso e $\sigma$ é a largura do pulso.

Utilizarei também as condições de contorno em que $\psi (0,t)=0$ e $\psi (L,t)=0$ o que garante que a função 'morra' nas pontas.

Utilizando estas condições iniciais e condições de contorno, foi feito um gif que mostra a evolução temporal da equação de Klein-Gordon utilizando o método das diferenças finitas: