Solução da Equação de Schrödinger para o átomo de Hidrogênio

Diagrama das coordenadas esféricas.

Em mecânica quântica, a solução da Equação de Schrödinger para o átomo de Hidrogênio é o processo de resolução da equação diferencial parcial formulada por Erwin Schrödinger em 1925 para o caso particular de duas partículas de cargas elétricas de mesmo módulo e sinal oposto (um elétron e um próton), em que a função energia potencial a que o elétron, a menor das duas partículas, está sujeito é da forma:

$V(r)=-\frac{k{e}^2}{r}$

Em que $r$ é a distância entre as partículas. Isto é, corresponde à interação, provocada pela força elétrica, entre próton e elétron característica de um átomo de Hidrogênio. Devido às consequências dramáticas se comparadas ao resultado do modelo clássico para essa situação, tal resolução é de extrema importância para a teoria atômica e, portanto, para a química, já que introduz os conceitos de orbital e números quânticos, a serem generalizados para outros elementos químicos, para além do Hidrogênio.

Para resolvê-la, a equação de Schrödinger independente do tempo é escrita em coordenadas esféricas, de modo que a função de onda $\psi$ seja tal que $\psi =\psi (r,\theta ,\varphi )$. Dessa forma a equação, genericamente escrita como

$-\frac{{\hslash }^2}{2m_e}{\mathrm{\nabla }}^2\psi +V\psi =E\psi$

será escrita da seguinte maneira:

$-\frac{{\hslash }^2}{2m_e}(\frac{1}{r}\frac{{\partial }^2}{\partial r^2}(r\psi )+\frac{1}{r^2\mathrm{sen}\theta }\frac{\partial }{\partial \theta }(\mathrm{sen}\theta \frac{\partial \psi }{\partial \theta })+\frac{1}{r^2{\mathrm{sen}}^2\theta }\frac{{\partial }^2\psi }{\partial {\varphi }^2})-\frac{k{e}^2}{r}\psi =E\psi$

Em que $m_e$ é a massa do elétron e $E$ é uma energia bem definida, isto é, representa os autovalores possíveis da energia, que, por sua vez, representam estados estacionários.

Esse problema é um caso particular do problema de força central em três dimensões espaciais em mecânica quântica, caracterizado por uma função energia potencial da forma $V=V(r)$, podendo este potencial ter outras formas que não somente proporcionais ao inverso do quadrado da distância.

Solução por separação de variáveis

Predefinição:Artigo principal Primeiramente, torna-se prático isolar o operador laplaciano do lado esquerdo da equação e retornar a escrever a energia potencial $-\frac{k{e}^2}{r}$ como $V(r)$, por compacticidade. Logo, a equação fica:

$\frac{1}{r}\frac{{\partial }^2}{\partial r^2}(r\psi )+\frac{1}{r^2\mathrm{sen}\theta }\frac{\partial }{\partial \theta }(\mathrm{sen}\theta \frac{\partial \psi }{\partial \theta })+\frac{1}{r^2{\mathrm{sen}}^2\theta }\frac{{\partial }^2\psi }{\partial {\varphi }^2}=\frac{2m_e}{{\hslash }^2}(V(r)-E)\psi$

Para encontrar as funções $\psi =\psi (r,\theta ,\varphi )$ que são soluções da equação, utiliza-se o método da separação de variáveis, baseando-se no teorema que afirma que toda solução de uma equação diferencial parcial linear pode ser escrita como uma combinação linear de um conjunto (talvez infinito) de soluções separáveis, isto é, que possuem a seguinte forma:

$\psi (r,\theta ,\varphi )=R(r)\mathrm{\Theta }(\theta )\mathrm{\Phi }(\varphi )$

Substituindo na equação de Schrödinger, percebe-se que a derivação só afetará uma das funções que compõem $\psi$, ou seja, já que não há derivadas cruzadas, somente uma dentre as três funções $R(r)$, $\mathrm{\Theta }(\theta )$ e $\mathrm{\Phi }(\varphi )$ serão afetadas para cada termo. Além disso, as derivadas parciais tornar-se-ão derivadas ordinárias, já que cada uma dessas funções depende somente da variável em relação a qual a função está sendo derivada, resultando em:

$\frac{\mathrm{\Theta }\mathrm{\Phi }}{r}\frac{d^2}{d{r}^2}(rR)+\frac{R\mathrm{\Phi }}{r^2\mathrm{sen}\theta }\frac{d}{d\theta }(\mathrm{sen}\theta \frac{d\mathrm{\Theta }}{d\theta })+\frac{R\mathrm{\Theta }}{r^2{\mathrm{sen}}^2\theta }\frac{d^2\mathrm{\Phi }}{d{\varphi }^2}=\frac{2m_e}{{\hslash }^2}(V(r)-E)R\mathrm{\Theta }\mathrm{\Phi }$

Multiplicando ambos os lados da equação por $\frac{r^2{\mathrm{sen}}^2\theta }{R\mathrm{\Theta }\mathrm{\Phi }}$, obtém-se:

$\frac{r{\mathrm{sen}}^2\theta }{R}\frac{d^2}{d{r}^2}(rR)+\frac{\mathrm{sen}\theta }{\mathrm{\Theta }}\frac{d}{d\theta }(\mathrm{sen}\theta \frac{d\mathrm{\Theta }}{d\theta })+\frac{1}{\mathrm{\Phi }}\frac{d^2\mathrm{\Phi }}{d{\varphi }^2}=\frac{2m_e}{{\hslash }^2}(V(r)-E)r^2{\mathrm{sen}}^2\theta$

Ao atentar-se, percebe-se que o termo que envolve o terceiro termo aditivo do lado esquerdo da equação é a única dependência em $\varphi$ da equação inteira. Isolando-o, obtém-se:

$\frac{1}{\mathrm{\Phi }}\frac{d^2\mathrm{\Phi }}{d{\varphi }^2}=\frac{2m_e}{{\hslash }^2}(V(r)-E)r^2{\mathrm{sen}}^2\theta -\frac{r{\mathrm{sen}}^2\theta }{R}\frac{d^2}{d{r}^2}(rR)-\frac{\mathrm{sen}\theta }{\mathrm{\Theta }}\frac{d}{d\theta }(\mathrm{sen}\theta \frac{d\mathrm{\Theta }}{d\theta })$

Como o lado esquerdo, depende somente de $\varphi$ e como o lado direito depende somente de $r$ e $\theta$, podemos escrever que o lado esquerdo é uma função $f=f(\varphi )$ e o direito uma função $g=g(r,\theta )$. Logo, vale a seguinte igualdade:

$f(\varphi )=g(r,\theta )$

Contudo, como as três coordenadas esféricas são variáveis independentes, a única maneira da igualdade entre as funções ser verdadeira é se ambas forem iguais a uma constante. Tal constante é tradicionalmente escrita como sendo $-m^2$, em que $m$ seja, a princípio, um número qualquer. Assim, as equações foram separadas, tornando-se duas equações diferenciais ordinárias distintas: uma dependente da variável $\varphi$ e outra dependente das variáveis $r$ e $\theta$ . Elas encontram-se abaixo:

Equação diferencial ordinária em $\varphi$ $\frac{1}{\mathrm{\Phi }}\frac{d^2\mathrm{\Phi }}{d{\varphi }^2}=-m^2$

Equação diferencial ordinária em $r$ e $\theta$ $-\frac{2m_e}{{\hslash }^2}(V(r)-E)r^2{\mathrm{sen}}^2\theta +\frac{r{\mathrm{sen}}^2\theta }{R}\frac{d^2}{d{r}^2}(rR)+\frac{\mathrm{sen}\theta }{\mathrm{\Theta }}\frac{d}{d\theta }(\mathrm{sen}\theta \frac{d\mathrm{\Theta }}{d\theta })=m^2$

Nota-se que a segunda equação ainda possui variáveis acopladas e, portanto, deve-se separá-las de maneira similar à anterior. Para isso, divide-se a segunda equação por $se{n}^2\theta$ e se coloca o terceiro termo do lado esquerdo para o lado direito da equação, obtendo a seguinte expressão:

$-\frac{2m_e}{{\hslash }^2}(V(r)-E)r^2+\frac{r}{R}\frac{d^2}{d{r}^2}(rR)+\frac{1}{\mathrm{\Theta }sen\theta }\frac{d}{d\theta }(\mathrm{sen}\theta \frac{d\mathrm{\Theta }}{d\theta })=\frac{m^2}{se{n}^2\theta }⟹-\frac{2m_e}{{\hslash }^2}(V(r)-E)r^2+\frac{r}{R}\frac{d^2}{d{r}^2}(rR)=\frac{m^2}{se{n}^2\theta }-\frac{1}{\mathrm{\Theta }sen\theta }\frac{d}{d\theta }(\mathrm{sen}\theta \frac{d\mathrm{\Theta }}{d\theta })$

Novamente, temos o lado esquerdo da equação dependente somente de $r$ e o lado direito dependente somente de $\theta$, podendo escrever que o lado esquerdo é uma função $F=F(r)$ e o direito uma função $G=G(\theta )$. Logo,é válida a seguinte igualdade:

$F(r)=G(\theta )$

E, mais uma vez, a única maneira de isso se verificar é se ambas forem iguais a uma constante. Essa constante, por motivos que ficarão claros mais adiante na resolução, é convenientemente escrita como $l(l+1)$, em que $l$ é um número qualquer, a princípio. Assim, obtemos duas novas equações:

$\frac{m^2}{se{n}^2\theta }-\frac{1}{\mathrm{\Theta }sen\theta }\frac{d}{d\theta }(\mathrm{sen}\theta \frac{d\mathrm{\Theta }}{d\theta })=l(l+1)$

$-\frac{2m_e}{{\hslash }^2}(V(r)-E)r^2+\frac{r}{R}\frac{d^2}{d{r}^2}(rR)=l(l+1)$

Na primeira equação, passando os termos da esquerda para a direita e multiplicando-a pela função incógnita $\mathrm{\Theta }$, obtemos:

Equação diferencial ordinária em $\theta$ $\frac{1}{sen\theta }\frac{d}{d\theta }(\mathrm{sen}\theta \frac{d\mathrm{\Theta }}{d\theta })+(l(l+1)-\frac{m^2}{se{n}^2\theta })\mathrm{\Theta }=0$

Na segunda equação, denominada equação radial, primeiramente se isola o segundo termo do lado esquerdo e depois se multiplica a equação por $\frac{R}{r}$, obtendo:

$\frac{d^2}{d{r}^2}(rR)=[\frac{l(l+1)}{r}+\frac{2m_e}{{\hslash }^2}(V(r)-E)r]R$

Por fim, fatorando o termo $\frac{2m_{e}r}{{\hslash }^2}$ do lado direito, conclui-se com a seguinte equação:

Equação diferencial ordinária em $r$ $\frac{d^2}{d{r}^2}(rR)=\frac{2m_e}{{\hslash }^2}[\frac{l(l+1){\hslash }^2}{2m_e{r}^2}+V(r)-E](rR)$

Solução da equação em $\varphi$

A equação mais simples e a primeira a ser obtida é a equação da coordenada angular $\varphi$ que, de fato, é idêntica à equação do oscilador harmônico simples.

Equação diferencial em $\varphi$ $\frac{1}{\mathrm{\Phi }}\frac{d^2\mathrm{\Phi }}{d{\varphi }^2}=-m^2$

É também possível escrevê-la, usando a notação ${\mathrm{\Phi }}^{″}=\frac{d^2\mathrm{\Phi }}{d{\varphi }^2}$, como:

${\mathrm{\Phi }}^{″}+m^2\mathrm{\Phi }=0$

Há duas maneiras principais de como escrever a solução geral dessa equação: como combinação linear de senos e cossenos ou como combinação linear de exponenciais imaginárias. Isto é:

$\mathrm{\Phi }(\varphi )=A\mathrm{sen}(m\varphi )+B\cos (m\varphi )$

$\mathrm{\Phi }(\varphi )=A{e}^{im\varphi }+B{e}^{-im\varphi }$

É conveniente escolher a segunda forma, admitindo a possibilidade que $m$ possa ser tanto positivo quanto negativo, e omitindo a constante $A$, devido à propriedade de normalização da função de onda. Assim, a solução fica:

$\mathrm{\Phi }(\varphi )=e^{im\varphi }$, $m\ge 0$ ou $m<0$.

Quantização de $m$

Usando a propriedade cíclica da função $\mathrm{\Phi }$, isto é, que $\mathrm{\Phi }(\varphi +2\pi )=\mathrm{\Phi }(\varphi )$, obtemos o seguinte:

$\mathrm{\Phi }(\varphi +2\pi )=\mathrm{\Phi }(\varphi )⟹e^{im(\varphi +2\pi )}=e^{im\varphi }⟹e^{im\varphi }.e^{im(2\pi )}=e^{im\varphi }⟹e^{im(2\pi )}=1$

Para que essa última exponencial imaginária seja igual a 1, $m$ tem de ser um número inteiro. Esse número é denominado o número quântico magnético.

$m=0,1,-1,2,-2,...$