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

$\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 },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,-\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,-\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 }$

$\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,-\mathrm{\Delta },t)-{\beta }^2\psi$