Problemas com Física Clássica

Por volta do início do século 20, a física clássica, baseada na mecânica newtoniana e nas equações de eletricidade e magnetismo de Maxwell, descrevia a natureza como a conhecíamos. A Mecânica Estatística também foi uma disciplina bem desenvolvida que descreve sistemas com um grande número de graus de liberdade. Por volta dessa época, Einstein introduziu a Relatividade Especial, que era compatível com as equações de Maxwell, mas mudou nossa compreensão do espaço-tempo e da mecânica modificada.

Muitas coisas permaneceram sem explicação. Embora o elétron como constituinte dos átomos tenha sido encontrado, a estrutura atômica era rica e bastante misteriosa. Havia problemas com a física clássica, incluindo a radiação do corpo negro, o efeito fotoelétrico, a teoria atômica básica, a dispersão de Compton e, eventualmente, a difração de todos os tipos de partículas. Plank formulou a hipótese de que a energia EM sempre foi emitida em quanta

$\ begin {displaymath} E = h \ nu = \ hbar \ omega \ end {displaymath}$

para resolver o problema do corpo negro. Muito mais tarde, deBroglie derivou o comprimento de onda das partículas.

$\ begin {displaymath} \ lambda = {h \ over p} \ end {displaymath}$

No final das contas, os problemas levaram ao desenvolvimento da Mecânica Quântica, em que todas as partículas têm comportamento tanto de onda quanto de partícula.

Pensamentos sobre experimentos de difração

A difração de fótons, elétrons e nêutrons foi observada e usada para estudar a estrutura cristalina.

Para entender a entrada experimental de uma forma simplificada, consideramos alguns experimentos mentais sobre a difração de fótons, elétrons e balas através de duas fendas. Por exemplo, os fótons, que compõem todas as ondas eletromagnéticas, mostram um padrão de difração exatamente como previsto pela teoria de ondas eletromagnéticas, mas sempre detectamos um número inteiro de fótons com a relação de Plank,

$\ bgroup \ color {black} $ E = h \ nu $ \ egroup$

satisfazendo entre a freqüência das ondas e a energia das partículas.

Elétrons, nêutrons e tudo o mais se comportam exatamente da mesma maneira, exibindo difração em forma de onda, mas detectando e satisfazendo um número inteiro de partículas. Esta fórmula de comprimento de onda deBroglie relaciona a propriedade da onda à propriedade da partícula.

$\ bgroup \ color {black} $ \ lambda = {h \ over p} $ \ egroup$ $\ bgroup \ color {black} $ \ lambda $ \ egroup$ $\ bgroup \ color {black} $ p $ \ egroup$

Amplitudes de probabilidade

Na Mecânica Quântica, entendemos essa dualidade onda-partícula usando amplitudes de probabilidade (complexas) que satisfazem uma equação de onda.

$\ begin {displaymath} \ bgroup \ color {black} \ psi (\ vec {x}, t) = e ^ {i (\ vec {k} \ cdot \ vec {x} - \ omega t)} = e ^ {i (\ vec {p} \ cdot \ vec {x} -Et) / \ hbar} \ egroup \ end {displaymath}$

A probabilidade de encontrar uma partícula em uma posição $\ bgroup \ color {black} $ \ vec {x} $ \ egroup$ em algum momento $\ bgroup \ color {black} $ t $ \ egroup$ é o quadrado absoluto da amplitude da probabilidade.

$\ bgroup \ color {black} $ \ psi (\ vec {x}, t) $ \ egroup$

$\ begin {displaymath} \ bgroup \ color {black} P (\ vec {x}, t) = \ left \ vert \ psi (\ vec {x}, t) \ right \ vert ^ 2 \ egroup \ end {displaymath }$

Para calcular a probabilidade de encontrar um elétron em nosso detector de experimento mental adicionamos, à amplitude da probabilidade de chegar ao detector pela fenda 1, a amplitude de chegar ao detector pela fenda 2 e obter o quadrado absoluto.

$\ begin {displaymath} \ bgroup \ color {preto} P _ {\ mathrm {detector}} = \ left \ vert \ psi_1 + \ psi_2 \ right \ vert ^ 2 \ egroup \ end {displaymath}$

A Mecânica Quântica muda completamente nossa visão do mundo. Em vez de um mundo determinístico, agora temos apenas probabilidades . Não podemos nem mesmo medir a posição e o momento de uma partícula (com precisão) ao mesmo tempo. A mecânica quântica exigirá que usemos a matemática de operadores, transformadas de Fourier, espaços vetoriais e muito mais.

Pacotes de onda e incerteza

A amplitude de probabilidade para uma partícula livre com momento $\ bgroup \ color {black} $ \ vec {p} $ \ egroup$ e energia é a função de onda complexa

$\ bgroup \ color {black} $ E = \ sqrt {(pc) ^ 2 + (mc ^ 2) ^ 2} $ \ egroup$

$\ begin {displaymath} \ bgroup \ color {black} \ psi _ {\ mathrm {partícula livre}} (\ vec {x}, t) = e ^ {i (\ vec {p} \ cdot \ vec {x} - Et) / \ hbar}. \ Egroup \ end {displaymath}$

Observe que em $\ bgroup \ color {black} $ \ vert \ psi \ vert ^ 2 = 1 $ \ egroup$ todos os lugares isso não representa uma partícula localizada. Na verdade, reconhecemos na propriedade de onda que, para termos exatamente uma frequência, uma onda deve ser espalhada no espaço.

Podemos construir pacotes de ondas localizados que representam partículas individuais somando essas funções de onda de partículas livres (com alguns coeficientes).

$\ begin {displaymath} \ bgroup \ color {black} \ psi (x, t) = {1 \ over \ sqrt {2 \ pi \ hbar}} \ int \ limits _ {- \ infty} ^ {+ \ infty} \ phi (p) e ^ {i (px-Et) / \ hbar} dp \ egroup \ end {displaymath}$

(Mudamos para uma dimensão para simplificar.) Da mesma forma, podemos calcular o coeficiente para cada momento

$\ begin {displaymath} \ bgroup \ color {black} \ phi (p) = {1 \ over \ sqrt {2 \ pi \ hbar}} \ int \ limits _ {- \ infty} ^ \ infty \ psi (x) e ^ {- ipx / \ hbar} dx. \ egroup \ end {displaymath}$

Esses coeficientes $\ bgroup \ color {black} $ \ phi (p) $ \ egroup$ são na verdade a função de estado da partícula no espaço de momento. Podemos descrever o estado de uma partícula em espaço de posição com $\ bgroup \ color {black} $ \ psi (x) $ \ egroup$ ou em espaço de momento com $\ bgroup \ color {black} $ \ phi (p) $ \ egroup$ .

Podemos usar $\ bgroup \ color {black} $ \ phi (p) $ \ egroup$ para calcular a função de distribuição de probabilidade para o momento.

$\ begin {displaymath} \ bgroup \ color {black} P (p) = \ left \ vert \ phi (p) \ right \ vert ^ 2 \ egroup \ end {displaymath}$

Mostraremos que pacotes de ondas como esses se comportam corretamente no limite clássico, justificando a escolha que fizemos.

$\ bgroup \ color {black} $ \ psi _ {\ mathrm free particle} (\ vec {x}, t) $ \ egroup$

O Princípio da Incerteza de Heisenberg é uma propriedade das ondas que podemos deduzir de nosso estudo de pacotes de ondas localizados.

$\ begin {displaymath} \ bgroup \ color {black} \ Delta p \ Delta x \ geq {\ hbar \ over 2} \ egroup \ end {displaymath}$

Isso mostra que, devido à natureza ondulatória das partículas, não podemos localizar uma partícula em um pequeno volume sem aumentar sua energia. Por exemplo, podemos estimar a energia do estado fundamental (e o tamanho) de um átomo de hidrogênio muito bem a partir do princípio da incerteza.

O próximo passo na construção da Mecânica Quântica é determinar como uma função de onda se desenvolve com o tempo - particularmente útil se um potencial for aplicado. A equação diferencial que as funções de onda devem satisfazer é chamada de Equação de Schrödinger.

Operadores

A Equação de Schrödinger vem diretamente de nossa compreensão dos pacotes de ondas. Para passar de pacotes de ondas a uma equação diferencial, usamos o novo conceito de operadores (lineares).

Determinamos os operadores de momento e energia requerendo que, quando um operador de alguma variável $\ bgroup \ color {black} $ v $ \ egroup$ venha a atuar em nossa função de onda simples, obtenhamos $\ bgroup \ color {black} $ v $ \ egroup$ vezes a mesma função de onda.

$\ begin {displaymath} \ bgroup \ color {preto} p_x ^ {(op)} = {\ hbar \ over i} {\ partial \ over \ partial x} \ egroup \ end {displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} p_x ^ {(op)} e ^ {i (\ vec {p} \ cdot \ vec {x} -Et) / \ h ... ... ec { x} -Et) / \ hbar} = p_xe ^ {i (\ vec {p} \ cdot \ vec {x} -Et) / \ hbar} \ egroup \ end {displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} E ^ {(op)} = i \ hbar {\ partial \ over \ partial t} \ egroup \ end {displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} E ^ {(op)} e ^ {i (\ vec {p} \ cdot \ vec {x} -Et) / \ hba ... ... \ vec {x} -Et) / \ hbar} = Ee ^ {i (\ vec {p} \ cdot \ vec {x} -Et) / \ hbar} \ egroup \ end {displaymath}$

Expectativas de valores

Podemos usar operadores para nos ajudar a calcular o valor esperado de uma variável física. Se uma partícula está no estado $\ bgroup \ color {black} $ \ psi (x) $ \ egroup$ , a maneira normal de calcular o valor esperado de $\ bgroup \ color {black} $ f (x) $ \ egroup$ é

$\ begin {displaymath} \ bgroup \ color {black} \ langle f (x) \ rangle = \ int \ limits _ {- \ inft ... ... = \ int \ limits _ {- \ infty} ^ \ infty \ psi ^ * (x) \ psi (x) f (x) dx. \ egroup \ end {displaymath}$

Se valor esperado da variável que desejamos calcular não for o de uma função simples de $\ bgroup \ color {black} $ x $ \ egroup$ (como $\ bgroup \ color {black} $ p $ \ egroup$ ) , deixe o operador $\ bgroup \ color {black} $ \ psi (x) $ \ egroup$ atuar

$\ begin {displaymath} \ bgroup \ color {black} \ langle p \ rangle = \ int \ limits _ {- \ infty} ^ \ infty \ psi ^ * (x) p ^ {(op)} \ psi (x) dx . \ egroup \ end {displaymath}$

Temos uma notação abreviada para o valor esperado de uma variável $\ bgroup \ color {black} $ v $ \ egroup$ no estado $\ bgroup \ color {black} $ \ psi $ \ egroup$ que é bastante útil.

$\ begin {displaymath} \ bgroup \ color {black} \ langle \ psi \ vert v \ vert \ psi \ rangle \ equiv \ int \ limits _ {- \ infty} ^ \ infty \ psi ^ * (x) v ^ {( op)} \ psi (x) dx. \ egroup \ end {displaymath}$

Estendemos a notação de apenas valores de expectativa para

$\ begin {displaymath} \ bgroup \ color {black} \ langle \ psi \ vert v \ vert \ phi \ rangle \ equiv \ int \ limits _ {- \ infty} ^ \ infty \ psi ^ * (x) v ^ {( op)} \ phi (x) dx \ egroup \ end {displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} \ langle \ psi \ vert \ phi \ rangle \ equiv \ int \ limits _ {- \ infty} ^ \ infty \ psi ^ * (x) \ phi (x) dx \ egroup \ end {displaymath}$

Usamos muito essa notação abreviada de Dirac Bra-Ket .

Comutadores

Operadores (ou variáveis na mecânica quântica) não necessariamente comutam. Podemos calcular o comutador de duas variáveis, por exemplo

$\ begin {displaymath} \ bgroup \ color {black} [p, x] \ equiv px-xp = {\ hbar \ over i}. \ egroup \ end {displaymath}$

Mais tarde, aprenderemos a derivar a relação de incerteza para duas variáveis de seu comutador. Também usaremos comutadores para resolver vários problemas importantes.

A Equação de Schrödinger

As funções de onda devem satisfazer a equação de Schrödinger, que na verdade é uma equação de onda.

$\ begin {displaymath} \ bgroup \ color {black} {- \ hbar ^ 2 \ over 2m} \ nabla ^ 2 \ psi (\ vec {x}, t ... ... vec {x}, t) = i \ hbar {\ partial \ psi (\ vec {x}, t) \ over \ partial t} \ egroup \ end {displaymath}$

Nós o usaremos para resolver muitos problemas neste curso. Em termos de operadores, isso pode ser escrito como

$\ begin {displaymath} \ bgroup \ color {black} H \ psi (\ vec {x}, t) = E \ psi (\ vec {x}, t) \ egroup \ end {displaymath}$

onde, descartando o rótulo $\ bgroup \ color {black} $ (op) $ \ egroup$ ,

$\ bgroup \ color {black} $ H = {p ^ 2 \ over 2m} + V (\ vec {x}) $ \ egroup$

é o operador hamiltoniano. Portanto, a Equação de Schrödinger é, em certo sentido, simplesmente a afirmação (nos operadores) de que a energia cinética mais a energia potencial é igual à energia total.

Autofunções, autovalores e espaços vetoriais

Para qualquer problema físico, as soluções da equação de Schrödinger separadas (entre tempo e espaço)

$\ bgroup \ color {black} $ \ psi (x, t) = u (x) T (t), $ \ egroup$

formam um conjunto extremamente importante. Se assumirmos que a equação se separa, obtemos as duas equações (em uma dimensão para simplificar)

$\ begin {displaymath} \ bgroup \ color {black} i \ hbar {\ partial T (t) \ over \ partial t} = E \; T (t) \ egroup \ end {displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} Hu (x) = E \; u (x) \ egroup \ end {displaymath}$

A segunda equação é chamada de equação de Schrödinger independente do tempo. Para estados ligados, existem apenas soluções para essa equação para algum conjunto quantizado de energias

$\ begin {displaymath} \ bgroup \ color {black} Hu_i (x) = E_iu_i (x). \ egroup \ end {displaymath}$

Para estados que não são limitados, é permitida uma faixa contínua de energia.

A equação de Schrödinger independente do tempo é um exemplo de uma equação de autovalor.

$\ begin {displaymath} \ bgroup \ color {black} H \ psi_i (\ vec {x}) = E_i \ psi_i (\ vec {x}) \ egroup \ end {displaymath}$

Se operarmos $\ bgroup \ color {black} $ \ psi_i $ \ egroup$ com $\ bgroup \ color {black} $ H $ \ egroup$ , obteremos a mesma função $\ bgroup \ color {black} $ \ psi_i $ \ egroup$ vezes alguma constante. Nesse caso, $\ bgroup \ color {black} $ \ psi_i $ \ egroup$ seria chamada de função autônoma e $\ bgroup \ color {black} $ E_i $ \ egroup$ seria chamada de autovalor. Normalmente, há um número infinito de soluções, indicado pelo índice $\ bgroup \ color {black} $ i $ \ egroup$ aqui.

Os operadores para variáveis físicas devem ter autovalores reais. Eles são chamados de operadores hermitianos. Podemos mostrar que as autofunções dos operadores Hermitianos são ortogonais (e podem ser normalizadas).

$\ begin {displaymath} \ bgroup \ color {preto} \ langle \ psi_i \ vert \ psi_j \ rangle = \ delta_ {ij} \ egroup \ end {displaymath}$

(No caso de autofunções com o mesmo autovalor, chamadas autofunções degeneradas, podemos escolher combinações lineares ortogonais entre si.) Vamos assumir que as autofunções também formam um conjunto completo de modo que qualquer função de onda possa ser expandida nelas

$\ begin {displaymath} \ bgroup \ color {preto} \ phi (\ vec {x}) = \ sum \ limits_i \ alpha_i \ psi_i (\ vec {x}) \ egroup \ end {displaymath}$

onde $\ bgroup \ color {black} $ \ alpha_i $ \ egroup$ são coeficientes que podem ser facilmente calculados (devido à ortonormalidade) por

$\ begin {displaymath} \ bgroup \ color {black} \ alpha_i = \ langle \ psi_i \ vert \ phi \ rangle. \ egroup \ end {displaymath}$

Portanto, agora temos outra maneira de representar um estado (além de espaço de posição e espaço de momento). Podemos representar um estado fornecendo os coeficientes em soma acima. (Observe que é apenas uma autofunção do operador momentum e é apenas uma autofunção do operador de posição (no espaço p), portanto, eles também representam uma expansão do estado em termos de autofunções.)

$\ bgroup \ color {black} $ \ psi_p (x) = e ^ {i (px-Et) / \ hbar} $ \ egroup$

$\ bgroup \ color {black} $ \ phi_x (p) = e ^ {- i (px-Et) / \ hbar} $ \ egroup$

Como a $\ bgroup \ color {black} $ \ psi_i $ \ egroup$ forma é um conjunto ortonormal completo, eles podem ser considerados como vetores unitários de um espaço vetorial. A função de onda arbitrária $\ bgroup \ color {black} $ \ phi $ \ egroup$ seria um vetor naquele espaço e poderia ser representada por seus coeficientes.

$\ begin {displaymath} \ bgroup \ color {black} \ phi = \ pmatrix {\ alpha_1 \ cr \ alpha_2 \ cr \ alpha_3 \ cr ... \ cr} \ egroup \ end {displaymath}$

O bra-ket pode ser pensado como um produto escalar entre o vetor arbitrário e um dos vetores unitários. Podemos usar a expansão em termos de estados próprios de energia para calcular muitas coisas. Em particular, uma vez que o tempo de desenvolvimento dos estados próprios de energia é muito simples,

$\ bgroup \ color {black} $ \ langle \ phi \ vert \ psi_i \ rangle $ \ egroup$ $\ bgroup \ color {black} $ \ phi $ \ egroup$

$\ begin {displaymath} \ bgroup \ color {black} \ psi (\ vec {x}, t) = \ psi (\ vec {x}) e ^ {- iE_it / \ hbar} \ egroup \ end {displaymath}$

podemos usar esses estados próprios para acompanhar o desenvolvimento de um estado arbitrário no tempo $\ bgroup \ color {black} $ \ phi $ \ egroup$

$\ begin {displaymath} \ bgroup \ color {black} \ phi (t) = \ pmatrix {\ alpha_1e ^ {- iE_1t / \ hbar} ... ... ^ {- iE_2t / \ hbar} \ cr \ alpha_3e ^ {-iE_3t / \ hbar} \ cr ... \ cr} \ egroup \ end {displaymath}$

simplesmente calculando os coeficientes $\ bgroup \ color {black} $ \ alpha_i $ \ egroup$ em $\ bgroup \ color {black} $ t = 0 $ \ egroup$ .

Podemos definir o conjugado hermitiano $\ bgroup \ color {black} $ O ^ \ dagger $ \ egroup$ do operador $\ bgroup \ color {black} $ O $ \ egroup$ por

$\ begin {displaymath} \ bgroup \ color {black} \ langle \ psi \ vert O \ vert \ psi \ rangle = \ lang ... ... vert O \ psi \ rangle = \ langle O ^ \ dagger \ psi \ vert \ psi \ rangle. \ egroup \ end {displaymath}$

Os operadores hermitianos $\ bgroup \ color {black} $ H $ \ egroup$ têm essa propriedade em que $\ bgroup \ color {black} $ H ^ \ dagger = H $ \ egroup$ .

Uma partícula em uma caixa

Como uma ilustração concreta dessas idéias, estudamos a partícula em uma caixa (em uma dimensão). Esta é apenas uma partícula (de massa $\ bgroup \ color {black} $ m $ \ egroup$ ) que é livre para se mover dentro das paredes de uma caixa $\ bgroup \ color {black} $ 0 <x <a $ \ egroup$ , mas que não pode penetrar nas paredes. Representamos isso por um potencial que é zero dentro da caixa e infinito fora. Resolvemos a equação de Schrödinger dentro da caixa e percebemos que a probabilidade de a partícula estar fora da caixa e, portanto, a função de onda lá, deve ser zero. Uma vez que não há potencial dentro, a equação de Schrödinger é

$\ begin {displaymath} \ bgroup \ color {black} Hu_n (x) = - {\ hbar ^ 2 \ over 2m} {d ^ 2u_n (x) \ over dx ^ 2} = E_n u_n (x) \ egroup \ end {displaymath}$

onde previmos que haverão muitas soluções indexadas por $\ bgroup \ color {black} $ n $ \ egroup$ . Conhecemos quatro funções (apenas 2 linearmente independentes) que têm um segundo derivado que é uma constante vezes a mesma função:

$\ bgroup \ color {black} $ u (x) = e ^ {ikx} $ \ egroup$ ,

$\ bgroup \ color {black} $ u (x) = e ^ {- ikx} $ \ egroup$ ,

$\ bgroup \ color {black} $ u (x) = \ sin (kx) $ \ egroup$ e

$\ bgroup \ color {black} $ u (x) = \ cos (kx) $ \ egroup$ .

A função de onda deve ser contínua, portanto necessitamos das condições de contorno

$\ begin {displaymath} \ bgroup \ color {black} u (0) = u (a) = 0. \ egroup \ end {displaymath}$

A função seno é sempre zero em $\ bgroup \ color {black} $ x = 0 $ \ egroup$ e nenhuma das outras é. Para tornar a função seno zero em $\ bgroup \ color {black} $ x = a $ \ egroup$ precisamos $\ bgroup \ color {black} $ ka = n \ pi $ \ egroup$ ou $\ bgroup \ color {black} $ k = {n \ pi \ over a} $ \ egroup$ . Assim, as autofunções de energia são dadas por

$\ begin {displaymath} \ bgroup \ color {preto} u_n (x) = C \ sin \ left ({n \ pi x \ over a} \ right) \ egroup \ end {displaymath}$

onde permitimos a constante geral $\ bgroup \ color {black} $ C $ \ egroup$ porque satisfaz à equação diferencial. Voltando a $\ bgroup \ color {black} $ \ sin \ left ({n \ pi x \ over a} \ right) $ \ egroup$ na equação de Schrödinger, descobrimos que

$\ begin {displaymath} \ bgroup \ color {black} E_n = {n ^ 2 \ pi ^ 2 \ hbar ^ 2 \ over 2ma ^ 2}. \ egroup \ end {displaymath}$

Somente energias quantizadas são permitidas quando resolvemos este problema de estado limitado. Temos uma tarefa restante. Os estados próprios devem ser normalizados para representar uma partícula.

$\ begin {displaymath} \ bgroup \ color {black} \ langle u_n \ vert u_n \ rangle = \ int \ limits_0 ... ... \ left ({n \ pi x \ over a} \ right) dx = \ vert C \ vert ^ 2 {a \ over 2} \ egroup \ end {displaymath}$

Portanto, a função de onda será normalizada se escolhermos

$\ bgroup \ color {black} $ C = \ sqrt {2 \ over a} $ \ egroup$

$\ begin {displaymath} \ bgroup \ color {black} u_n (x) = \ sqrt {2 \ over a} \ sin \ left ({n \ pi x \ over a} \ right) \ egroup \ end {displaymath}$

Sempre podemos multiplicar por qualquer número complexo de magnitude 1, mas isso não muda a física. Este exemplo mostra muitos dos recursos que veremos para outros problemas de estado vinculado. A única diferença é que, por causa de uma mudança infinita no potencial nas paredes da caixa, não precisamos manter a primeira derivada da função de onda contínua. Em todos os outros problemas, teremos que prestar mais atenção a isso.

Potenciais constantes por partes em uma dimensão

Agora estudamos a física de vários potenciais simples em uma dimensão . Primeiro, uma série de potenciais constantes por partes para os quais a equação de Schrödinger é

$\ begin {displaymath} \ bgroup \ color {black} {- \ hbar ^ 2 \ over 2m} {d ^ 2u (x) \ over dx ^ 2} + Vu (x) = Eu (x) \ egroup \ end { displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} {d ^ 2u (x) \ over dx ^ 2} + {2m \ over \ hbar ^ 2} (EV) u (x) = 0 \ egroup \ end {displaymath }$

e a solução geral , para $\ bgroup \ color {black} $ E> V $ \ egroup$ , pode ser escrita como

$\ begin {displaymath} \ bgroup \ color {black} u (x) = Ae ^ {ikx} + Be ^ {- ikx} \ egroup \ end {displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} u (x) = A \ sin (kx) + B \ cos (kx) \ egroup \ end {displaymath}$

com

$\ bgroup \ color {black} $ k = \ sqrt {2m (EV) \ over \ hbar ^ 2} $ \ egroup$

Também precisaremos de soluções para as regiões classicamente proibidas, onde a energia total é menor do que a energia potencial.

$\ bgroup \ color {black} $ E <V $ \ egroup$

$\ begin {displaymath} \ bgroup \ color {black} u (x) = Ae ^ {\ kappa x} + Be ^ {- \ kappa x} \ egroup \ end {displaymath}$

com

$\ bgroup \ color {black} $ \ kappa = \ sqrt {2m (VE) \ over \ hbar ^ 2} $ \ egroup$

(Ambos $\ bgroup \ color {black} $ k $ \ egroup$ e $\ bgroup \ color {black} $ \ kappa $ \ egroup$ são números reais positivos.) Os problemas de espalhamento 1D costumam ser análogos aos problemas em que a luz é refletida ou transmitida na superfície do vidro.

Em primeiro lugar, calcula-se a probabilidade da partícula de energia $\ bgroup \ color {black} $ E $ \ egroup$ refletir-se por um passo de potencial de altura

$\ bgroup \ color {black} $ V_0 $ \ egroup$ : $\ bgroup \ color {black} $ P_R = \ left ({\ sqrt {E} - \ sqrt {E-V_0} \ over \ sqrt {E} + \ sqrt {E-V_0}} \ right) ^ 2 $ \ egroup$ .

Também usamos este exemplo para entender a probabilidade atual.

$\ bgroup \ color {black} $ j = {\ hbar \ over 2im} [u ^ * {du \ over dx} - {du ^ * \ over dx} u] $ \ egroup$

Em segundo lugar, investigamos o poço de potencial quadrado do poço de potencial quadrado ( $\ bgroup \ color {black} $ V (x) = - V_0 $ \ egroup$ para $\ bgroup \ color {black} $ - a <x <a $ \ egroup$ e $\ bgroup \ color {black} $ V (x) = 0 $ \ egroup$ em outros lugares), para o caso em que a partícula não está ligada $\ bgroup \ color {black} $ E> 0 $ \ egroup$ . Assumindo um feixe de partículas incidente da esquerda, precisamos combinar as soluções nas três regiões nos limites em $\ bgroup \ color {black} $ x = \ pm um $ \ egroup$ . Depois de alguma aritmética difícil, as probabilidades de serem transmitidas ou refletidas são calculadas. Verifica-se que a probabilidade de ser transmitido vai para 1 para algumas energias particulares.

$\ begin {displaymath} \ bgroup \ color {black} E = -V_0 + {n ^ 2 \ pi ^ 2 \ hbar ^ 2 \ over 8ma ^ 2} \ egroup \ end {displaymath}$

Este tipo de comportamento é exibido pelo espalhamento de elétrons dos átomos. Em algumas energias, a probabilidade de espalhamento vai para zero.

Terceiro, estudamos a barreira de potencial quadrada ( $\ bgroup \ color {black} $ V (x) = + V_0 $ \ egroup$ para $\ bgroup \ color {black} $ - a <x <a $ \ egroup$ e $\ bgroup \ color {black} $ V (x) = 0 $ \ egroup$ em outro lugar), para o caso em que $\ bgroup \ color {black} $ E <V_0 $ \ egroup$ . Classicamente a probabilidade de ser transmitida seria zero, uma vez que a partícula é excluída energeticamente de estar dentro da barreira. O cálculo quântico dá a probabilidade de ser transmitido através da barreira como

$\ begin {displaymath} \ bgroup \ color {black} \ vert T \ vert ^ 2 = {(2k \ kappa) ^ 2 \ over (k ^ 2 + \ k ... ... htarrow ({4k \ kappa \ sobre k ^ 2 + \ kappa ^ 2}) ^ 2 e ^ {- 4 \ kappa a} \ egroup \ end {displaymath}$

onde

$\ bgroup \ color {black} $ k = \ sqrt {2mE \ over \ hbar ^ 2} $ \ egroup$ $\ bgroup \ color {black} $ \ kappa = \ sqrt {2m (V_0-E) \ over \ hbar ^ 2} $ \ egroup$

O estudo desta expressão mostra que a probabilidade de ser transmitido diminui à medida que a barreira fica mais alta ou mais larga. No entanto, a penetração da barreira é um fenômeno quântico importante.

Também estudamos bem o quadrado para o caso de estado ligado em que $\ bgroup \ color {black} $ E <0 $ \ egroup$ . Aqui, precisamos resolver uma equação transcendental para determinar as energias do estado vinculado. O número de estados vinculados aumenta com a profundidade e a largura do poço, mas sempre há pelo menos um estado vinculado.

O oscilador harmônico em uma dimensão

Em seguida, resolvemos os autoestados de energia do potencial do oscilador harmônico , onde eliminamos a constante da mola usando a frequência clássica do oscilador . Os valores próprios de energia são

$\ bgroup \ color {black} $ V (x) = {1 \ over 2} kx ^ 2 = {1 \ over 2} m \ omega ^ 2x ^ 2 $ \ egroup$ $\ bgroup \ color {black} $ k $ \ egroup$

$\ bgroup \ color {black} $ \ omega = \ sqrt {k \ over m} $ \ egroup$

$\ begin {displaymath} \ bgroup \ color {black} E_n = \ left (n + {1 \ over 2} \ right) \ hbar \ omega. \ egroup \ end {displaymath}$

Os autoestados de energia acabam sendo um polinômio (em $\ bgroup \ color {black} $ x $ \ egroup$ ) de $\ bgroup \ color {black} $ n $ \ egroup$ tempos de grau . Portanto, o estado fundamental, devidamente normalizado, é de apenas $\ bgroup \ color {black} $ e ^ {- m \ omega x ^ 2 / \ hbar} $ \ egroup$

$\ begin {displaymath} \ bgroup \ color {black} u_0 (x) = \ left ({m \ omega \ over \ pi \ hbar} \ right) ^ {1 \ over 4} e ^ {- m \ omega x ^ 2 / \ hbar}. \ Egroup \ end {displaymath}$

Posteriormente, retornaremos o oscilador harmônico para resolver o problema pelos métodos do operador.

Potenciais de função delta em uma dimensão

O potencial de função delta é muito útil para fazer modelos simples de moléculas e sólidos. Primeiro, resolvemos o problema com uma função delta atrativa

$\ bgroup \ color {black} $ V (x) = - aV_0 \ delta (x) $ \ egroup$

Como o estado ligado tem energia negativa, as soluções normalizáveis são $\ bgroup \ color {black} $ Ce ^ {\ kappa x} $ \ egroup$ para $\ bgroup \ color {black} $ x <0 $ \ egroup$ e $\ bgroup \ color {black} $ Ce ^ {- \ kappa x} $ \ egroup$ para $\ bgroup \ color {black} $ x> 0 $ \ egroup$ . Tornar $\ bgroup \ color {black} $ u (x) $ \ egroup$ contínuo e sua primeira derivada ter uma descontinuidade calculada a partir da equação de Schrödinger em $\ bgroup \ color {black} $ x = 0 $ \ egroup$ , nos dá exatamente um estado ligado com

$\ begin {displaymath} \ bgroup \ color {black} E = - {ma ^ 2V_0 ^ 2 \ over 2 \ hbar ^ 2}. \ egroup \ end {displaymath}$

Em seguida, usamos duas funções delta para modelar uma molécula. Resolvendo esse problema combinando funções de onda nas fronteiras em $\ bgroup \ color {black} $ \ pm d $ \ egroup$ , encontramos novamente equações transcendentais para duas energias de estado vinculadas. A energia do estado fundamental é mais negativa do que para uma função delta, indicando que a molécula seria ligada. Uma olhada na função de onda mostra que as duas funções delta podem diminuir a energia cinética em comparação com o estado de uma função delta, reduzindo a curvatura da função de onda. O estado excitado tem mais curvatura do que o estado atômico, então não esperaríamos ligações moleculares nesse estado.

$\ bgroup \ color {black} $ V (x) = - aV_0 \ delta (x + d) -aV_0 \ delta (xd) $ \ egroup$

Nosso potencial 1D final é o modelo de um sólido.

$\ begin {displaymath} \ bgroup \ color {black} V (x) = - aV_0 \ sum \ limits_ {n = - \ infty} ^ \ infty \ delta (x-na) \ egroup \ end {displaymath}$

Isso tem uma matriz periódica infinita de funções delta, então isso pode ser aplicável a um cristal. A solução para isso é um pouco complicada, mas se resume a

$\ begin {displaymath} \ bgroup \ color {black} \ cos (\ phi) = \ cos (ka) + {2maV_0 \ over \ hbar ^ 2k} \ sin (ka). \ egroup \ end {displaymath}$

Como o lado direito da equação pode ser maior que 1,0 (ou menor que -1), existem regiões para as quais não há soluções. Também existem faixas de energias com soluções. Essas faixas de energia são vistas em cristais (como Si).

$\ bgroup \ color {black} $ E = {\ hbar ^ 2k ^ 2 \ over 2m} $ \ egroup$

Solução de oscilador harmônico com operadores

Podemos resolver o problema do oscilador harmônico usando métodos de operadores. Escrevemos o hamiltoniano em termos do operador

$\ begin {displaymath} \ bgroup \ color {black} A \ equiv \ left (\ sqrt {m \ omega \ over2 \ hbar} x + i {p \ over \ sqrt {2m \ hbar \ omega}} \ right) \ egroup \ end {displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} H = {p ^ 2 \ over 2m} + {1 \ over 2} m \ omega ^ 2x ^ 2 = \ hbar \ omega (A ^ \ dagger A + {1 \ sobre 2}) \ egroup \ end {displaymath}$

Calculamos os comutadores

$\ begin {displaymath} \ bgroup \ color {black} [A, A ^ \ dagger] = {i \ over 2 \ hbar} (- [x, p] + [p, x]) = 1 \ egroup \ end { displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} [H, A] = \ hbar \ omega [A ^ \ dagger A, A] = \ hbar \ omega [A ^ \ dagger, A] A = - \ hbar \ omega A \ egroup \ end {displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} [H, A ^ \ dagger] = \ hbar \ omega [A ^ \ dagger A, A ^ \ ... ... \ hbar \ omega A ^ \ dagger [ A, A ^ \ dagger] = \ hbar \ omega A ^ \ dagger \ egroup \ end {displaymath}$

Se aplicamos o comutador $\ bgroup \ color {black} $ [H, A] $ \ egroup$ à autofunção $\ bgroup \ color {black} $ u_n $ \ egroup$ , obtemos $\ bgroup \ color {black} $ [H, A] u_n = - \ hbar \ omega Au_n $ \ egroup$ o qual rearranja a equação de autovalor

$\ begin {displaymath} \ bgroup \ color {black} H (Au_n) = (E_n- \ hbar \ omega) (Au_n). \ egroup \ end {displaymath}$

Isso diz que $\ bgroup \ color {black} $ (Au_n) $ \ egroup$ é uma autofunção de $\ bgroup \ color {black} $ H $ \ egroup$ com autovalor $\ bgroup \ color {black} $ (E_n- \ hbar \ omega) $ \ egroup$ , portanto, reduz a energia em $\ bgroup \ color {black} $ \ hbar \ omega $ \ egroup$ . Uma vez que a energia deve ser positiva para este hamiltoniano, a redução deve parar em algum lugar, no estado fundamental, onde teremos

$\ begin {displaymath} \ bgroup \ color {black} Au_0 = 0. \ egroup \ end {displaymath}$

Isso nos permite calcular a energia do estado fundamental como

$\ begin {displaymath} \ bgroup \ color {black} Hu_0 = \ hbar \ omega (A ^ \ dagger A + {1 \ over 2}) u_0 = {1 \ over 2} \ hbar \ omega u_0 \ egroup \ end {displaymath }$

mostrando que a energia do estado fundamental é $\ bgroup \ color {black} $ {1 \ over 2} \ hbar \ omega $ \ egroup$ . Da mesma forma, $\ bgroup \ color {black} $ A ^ \ dagger $ \ egroup$ aumenta a energia em $\ bgroup \ color {black} $ \ hbar \ omega $ \ egroup$ . Podemos subir e descer a escada de energia usando $\ bgroup \ color {black} $ A ^ \ dagger $ \ egroup$ e $\ bgroup \ color {black} $ A $ \ egroup$ , sempre em etapas de $\ bgroup \ color {black} $ \ hbar \ omega $ \ egroup$ . Os valores próprios de energia são, portanto,

$\ begin {displaymath} \ bgroup \ color {black} E_n = \ left (n + {1 \ over 2} \ right) \ hbar \ omega. \ egroup \ end {displaymath}$

Um pouco mais de cálculos mostra que

$\ begin {displaymath} \ bgroup \ color {black} Au_n = \ sqrt {n} u_ {n-1} \ egroup \ end {displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} A ^ \ punhal u_n = \ sqrt {n + 1} u_ {n + 1}. \ egroup \ end {displaymath}$

Essas fórmulas são úteis para todos os tipos de cálculos dentro do importante sistema oscilador harmônico. Ambos $\ bgroup \ color {black} $ p $ \ egroup$ e $\ bgroup \ color {black} $ x $ \ egroup$ podem ser escritos em termos de $\ bgroup \ color {black} $ A $ \ egroup$ e $\ bgroup \ color {black} $ A ^ \ dagger $ \ egroup$ .

$\ begin {displaymath} \ bgroup \ color {black} x = \ sqrt {\ hbar \ over 2m \ omega} (A + A ^ \ dagger) \ egroup \ end {displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} p = -i \ sqrt {m \ hbar \ omega \ over 2} (AA ^ \ dagger) \ egroup \ end {displaymath}$

Mais diversão com os operadores

Encontramos o operador de desenvolvimento de tempo resolvendo a equação .

$\ bgroup \ color {black} $ i \ hbar {\ partial \ psi \ over \ partial t} = H \ psi $ \ egroup$

$\ begin {displaymath} \ bgroup \ color {black} \ psi (t) = e ^ {- iHt / \ hbar} \ psi (t = 0) \ egroup \ end {displaymath}$

Isso implica que $\ bgroup \ color {black} $ e ^ {- iHt / \ hbar} $ \ egroup$ é o operador de desenvolvimento de tempo. Em alguns casos, podemos calcular o operador real da série de potências para o exponencial.

$\ begin {displaymath} \ bgroup \ color {black} e ^ {- iHt / \ hbar} = \ sum \ limits_ {n = 0} ^ \ infty {(- iHt / \ hbar) ^ n \ over n!} \ egroup \ end {displaymath}$

Temos trabalhado no que é chamado de quadro de Schrödinger, no qual as funções de onda (ou estados) se desenvolvem com o tempo. Existe a imagem alternativa de Heisenberg em que os operadores se desenvolvem com o tempo enquanto os estados não mudam.

Por exemplo, se quisermos calcular o valor esperado do operador $\ bgroup \ color {black} $ B $ \ egroup$ como uma função do tempo na imagem de Schrödinger usual, obtemos

$\ begin {displaymath} \ bgroup \ color {black} \ langle \ psi (t) \ vert B \ vert \ psi (t) \ rangle = ... ... 0) \ vert e ^ {iHt / \ hbar} Seja ^ {- iHt / \ hbar} \ vert \ psi (0) \ rangle. \ Egroup \ end {displaymath}$

No quadro de Heisenberg, o operador

$\ bgroup \ color {black} $ B (t) = e ^ {iHt / \ hbar} Be ^ {- iHt / \ hbar} $ \ egroup$

Usamos métodos de operação para calcular a relação de incerteza entre as variáveis não comutáveis

$\ begin {displaymath} \ bgroup \ color {black} (\ Delta A) (\ Delta B) \ geq {i \ over 2} \ langle [A, B] \ rangle \ egroup \ end {displaymath}$

o que dá o resultado que deduzimos dos pacotes de ondas para $\ bgroup \ color {black} $ p $ \ egroup$ e $\ bgroup \ color {black} $ x $ \ egroup$ .

Novamente, usamos métodos de operador para calcular a derivada de tempo de um valor esperado.

$\ begin {displaymath} \ bgroup \ color {black} {d \ over dt} \ langle \ psi \ vert A \ vert \ psi \ ra ... ... rtial A \ over \ partial t} \ right \ vert \ psi \ right \ rangle_ \ psi \ egroup \ end {displaymath}$

(A maioria dos operadores que usamos não tem dependência de tempo explícita, então o segundo termo geralmente é zero.) Isso novamente mostra a importância do operador hamiltoniano para o desenvolvimento do tempo. Podemos usar isso para mostrar que na mecânica quântica os valores esperados $\ bgroup \ color {black} $ p $ \ egroup$ e $\ bgroup \ color {black} $ x $ \ egroup$ se comportam como esperaríamos da mecânica newtoniana (teorema de Ehrenfest) .

$\ begin {displaymath} \ bgroup \ color {black} {d \ langle x \ rangle \ over dt} = {i \ over \ hbar} ... ... 2 \ over 2m}, x] \ rangle = \ left \ langle {p \ over m} \ right \ rangle \ egroup \ end {displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} {d \ langle p \ rangle \ over dt} = {i \ over \ hbar} ... ... \ right \ rangle = - \ left \ langle {dV ( x) \ over dx} \ right \ rangle \ egroup \ end {displaymath}$

Qualquer operador $\ bgroup \ color {black} $ A $ \ egroup$ que comuta com o hamiltoniano tem um valor de expectativa independente do tempo . As autofunções de energia também podem ser autofunções (simultâneas) do operador de comutação $\ bgroup \ color {black} $ A $ \ egroup$ . Geralmente é uma simetria do $\ bgroup \ color {black} $ H $ \ egroup$ que leva a um operador de comutação e, portanto, uma constante adicional do movimento.

Duas partículas em 3 dimensões

Até agora, trabalhamos com estados de apenas uma partícula em uma dimensão. A extensão para duas partículas diferentes e para três dimensões é direta. As coordenadas e os momentos de diferentes partículas e das dimensões adicionais comutam entre si como poderíamos esperar da física clássica. As únicas coisas que não comutam são uma coordenada com seu impulso, por exemplo.

$\ begin {displaymath} \ bgroup \ color {black} [p _ {(2) z}, z _ {(2)}] = {\ hbar \ over i} \ egroup \ end {displaymath}$

enquanto

$\ begin {displaymath} \ bgroup \ color {black} [p _ {(1) x}, x _ {(2)}] = [p _ {(2) z}, y _ {(2)}] = 0. \ egroup \ end {displaymath}$

Podemos escrever estados para duas partículas que não estão correlacionadas, como $\ bgroup \ color {black} $ u_0 (\ vec {x} _ {(1)}) u_3 (\ vec {x} _ {(2)}) $ \ egroup$ , ou podemos escrever estados em que as partículas estão correlacionadas. O hamiltoniano para duas partículas em 3 dimensões simplesmente se torna

$\ begin {displaymath} \ bgroup \ color {black} H = {- \ hbar ^ 2 \ over 2m _ {(1)}} \ left ({\ parti ... ... artial z _ {(2)} ^ 2 } \ right) + V (\ vec {x} _ {(1)}, \ vec {x} _ {(2)}) \ egroup \ end {displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} H = {- \ hbar ^ 2 \ over 2m _ {(1)}} \ nabla ^ 2 _ {(1)} + ... ... 2m _ {(2) }} \ nabla ^ 2 _ {(1)} + V (\ vec {x} _ {(1)}, \ vec {x} _ {(2)}) \ egroup \ end {displaymath}$

Se duas partículas interagem entre si, sem potencial externo,

$\ begin {displaymath} \ bgroup \ color {black} H = {- \ hbar ^ 2 \ over 2m _ {(1)}} \ nabla ^ 2 _ {(1)} + ... ... 2m _ {(2) }} \ nabla ^ 2 _ {(1)} + V (\ vec {x} _ {(1)} - \ vec {x} _ {(2)}) \ egroup \ end {displaymath}$

o hamiltoniano tem uma simetria translacional e permanece invariante sob a tradução . Podemos mostrar que esta simetria translacional implica na conservação do momento total . Da mesma forma, mostraremos que a simetria rotacional implica na conservação do momento angular e que a simetria do tempo implica na conservação da energia. $\ bgroup \ color {black} $ \ vec {x} \ rightarrow \ vec {x} + \ vec {a} $ \ egroup$

Para duas partículas interagindo através de um potencial que depende apenas da diferença nas coordenadas,

$\ begin {displaymath} \ bgroup \ color {black} H = {\ vec {p} _1 ^ 2 \ over 2m} + {\ vec {p} _2 ^ 2 \ over 2m} + V (\ vec {r} _1 - \ vec {r} _2) \ egroup \ end {displaymath}$

podemos fazer a transformação usual para o centro de massa feita na mecânica clássica

$\ begin {displaymath} \ bgroup \ color {preto} \ vec {r} \ equiv \ vec {r} _1- \ vec {r} _2 \ egroup \ end {displaymath}$

$\ begin {displaymath} \ bgroup \ color {preto} \ vec {R} \ equiv {m_1 \ vec {r} _1 + m_2 \ vec {r} _2 \ over m_1 + m_2} \ egroup \ end {displaymath}$

e reduzir o problema para o CM se movendo como uma partícula livre

$\ begin {displaymath} \ bgroup \ color {black} M = m_1 + m_2 \ egroup \ end {displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} H = {- \ hbar ^ 2 \ over 2M} \ vec {\ nabla} _R ^ 2 \ egroup \ end {displaymath}$

mais um problema potencial em 3 dimensões com a massa reduzida usual.

$\ begin {displaymath} \ bgroup \ color {black} {1 \ over \ mu} = {1 \ over m_1} + {1 \ over m_2} \ egroup \ end {displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} H = - {\ hbar ^ 2 \ over 2 \ mu} \ vec {\ nabla} _r ^ 2 + V (\ vec {r}) \ egroup \ end {displaymath }$

Portanto, agora temos um problema 3D para resolver (3 variáveis em vez de 6).

Partículas Idênticas

Partículas idênticas nos apresentam outra simetria na natureza. Os elétrons, por exemplo, são indistinguíveis uns dos outros, então devemos ter uma simetria do hamiltoniano sob o intercâmbio de qualquer par de elétrons. Vamos chamar a operadora que intercambia elétron-1 e elétron-2 $\ bgroup \ color {black} $ P_ {12} $ \ egroup$

$\ begin {displaymath} \ bgroup \ color {black} [H, P_ {12}] = 0 \ egroup \ end {displaymath}$

Portanto, podemos fazer com que nossos estados próprios de energia também sejam estados próprios de $\ bgroup \ color {black} $ P_ {12} $ \ egroup$ . É fácil ver (operando em um estado próprio duas vezes com $\ bgroup \ color {black} $ P_ {12} $ \ egroup$ ) que os valores próprios possíveis são $\ bgroup \ color {black} $ \ pm 1 $ \ egroup$ . É uma lei da física que $\ bgroup \ color {black} $ {1 \ over 2} $ \ egroup$ partículas de spin chamadas férmions (como elétrons) sempre são antissimétricas em intercâmbio , enquanto partículas com spin inteiro chamadas de bósons (como fótons) sempre são simétricas em intercâmbio . A anti-simetria sob intercâmbio leva ao princípio de exclusão de Pauli de que dois elétrons (por exemplo) não podem estar no mesmo estado.

Alguns problemas 3D separáveis em coordenadas cartesianas

Começamos nosso estudo da Mecânica Quântica em 3 dimensões com alguns casos simples de problemas que podem ser separados em coordenadas cartesianas. Isso é possível quando o hamiltoniano pode ser escrito

$\ begin {displaymath} \ bgroup \ color {black} H = H_x + H_y + H_z. \ egroup \ end {displaymath}$

Um bom exemplo de separação de variável em coordenadas cartesianas é o oscilador harmônico 3D

$\ begin {displaymath} \ bgroup \ color {black} V (r) = {1 \ over 2} m \ omega ^ 2r ^ 2 \ egroup \ end {displaymath}$

que tem energias que dependem de três números quânticos.

$\ begin {displaymath} \ bgroup \ color {black} E_ {n_xn_yn_z} = \ left (n_x + n_y + n_z + {3 \ over 2} \ right) \ hbar \ omega \ egroup \ end {displaymath}$

Ele realmente se comporta como 3 osciladores harmônicos unidimensionais independentes.

Outro problema que separa é a partícula em uma caixa 3D . Novamente, as energias dependem de três números quânticos

$\ begin {displaymath} \ bgroup \ color {black} E_ {n_xn_yn_z} = {\ pi ^ 2 \ hbar ^ 2 \ over 2mL ^ 2} \ left (n_x ^ 2 + n_y ^ 2 + n_z ^ 2 \ right) \ egroup \ end {displaymath}$

para uma caixa cúbica de lado $\ bgroup \ color {black} $ L $ \ egroup$ . Investigamos o efeito do princípio de exclusão de Pauli preenchendo nossa caixa 3D com férmions idênticos, que devem estar todos em estados diferentes. Podemos usar isso para modelar anãs brancas ou estrelas de nêutrons.

Na física clássica, são necessárias três coordenadas para fornecer a localização de uma partícula em 3D. Na mecânica quântica, estamos descobrindo que são necessários três números quânticos para rotular e autoestar de energia (sem incluir o spin).

Momentum Angular

Para o comum problema de potenciais centrais $\ bgroup \ color {black} $ V (r) $ \ egroup$ , usamos o óbvio simetria rotacional ao descobrir que o momento angular , , operadores comutar com , $\ bgroup \ color {black} $ \ vec {L} = \ vec {x} \ times \ vec {p} $ \ egroup$ $\ bgroup \ color {black} $ H $ \ egroup$

$\ begin {displaymath} \ bgroup \ color {black} [H, L_z] = [H, L_x] = [H, L_y] = 0 \ egroup \ end {displaymath}$

mas eles não se deslocam entre si.

$\ begin {displaymath} \ bgroup \ color {black} [L_x, L_y] \ neq 0 \ egroup \ end {displaymath}$

Queremos encontrar dois operadores mutuamente comutantes que comutam com $\ bgroup \ color {black} $ H $ \ egroup$ , então voltamos para $\ bgroup \ color {black} $ L ^ 2 = L_x ^ 2 + L_y ^ 2 + L_z ^ 2 $ \ egroup$ o qual comuta com cada componente de $\ bgroup \ color {black} $ L $ \ egroup$ .

$\ begin {displaymath} \ bgroup \ color {black} [L ^ 2, L_z] = 0 \ egroup \ end {displaymath}$

Escolhemos nossos dois operadores para serem $\ bgroup \ color {black} $ L ^ 2 $ \ egroup$ e $\ bgroup \ color {black} $ L_z $ \ egroup$ .

Alguns cálculos revelam que podemos escrever

$\ begin {displaymath} \ bgroup \ color {black} p ^ 2 = {1 \ over r ^ 2} \ left (L ^ 2 + (\ vec {r} \ cdot \ vec {p}) ^ 2-i \ hbar \ vec {r} \ cdot \ vec {p} \ right). \ egroup \ end {displaymath}$

Com isso, a parte da energia cinética de nossa equação só terá derivadas ao $\ bgroup \ color {black} $ r $ \ egroup$ assumir que temos autoestados de $\ bgroup \ color {black} $ L ^ 2 $ \ egroup$ .

$\ begin {displaymath} \ bgroup \ color {black} {- \ hbar ^ 2 \ over 2 \ mu} \ left [{1 \ over r ^ 2} \ le ... ... r ^ 2} \ right] u_E (\ vec {r}) + V (r) u_E (\ vec {r}) = Eu_E (\ vec {r}) \ egroup \ end {displaymath}$

A equação de Schrödinger, portanto, separa em uma parte angular (o $\ bgroup \ color {black} $ L ^ 2 $ \ egroup$ termo) e uma parte radial (o resto). Com esta separação, esperamos (antecipando um pouco a solução angular)

$\ begin {displaymath} \ bgroup \ color {black} u_E (\ vec {r}) = R_ {E \ ell} (r) Y _ {\ ell m} (\ theta, \ phi) \ egroup \ end {displaymath}$

será uma solução. O termo $\ bgroup \ color {black} $ Y _ {\ ell m} (\ theta, \ phi) $ \ egroup$ será de autofunções $\ bgroup \ color {black} $ L ^ 2 $ \ egroup$

$\ begin {displaymath} \ bgroup \ color {black} L ^ 2 Y _ {\ ell m} (\ theta, \ phi) = \ ell (\ ell + 1) \ hbar ^ 2 Y _ {\ ell m} (\ theta , \ phi) \ egroup \ end {displaymath}$

então a equação radial se torna

$\ begin {displaymath} \ bgroup \ color {black} {- \ hbar ^ 2 \ over 2 \ mu} \ left [{1 \ over r ^ 2} \ le ... ... r ^ 2} \ right] R_ {E \ ell} (r) + V (r) R_ {E \ ell} (r) = ER_ {E \ ell} (r) \ egroup \ end {displaymath}$

Devemos voltar a esta equação para cada $\ bgroup \ color {black} $ V (r) $ \ egroup$ que desejamos resolver.

Nós resolver a parte angular do problema, em geral, usando operadores momento angular. Descobrimos que o momento angular é quantizado .

$\ begin {displaymath} \ bgroup \ color {black} L_zY _ {\ ell m} (\ theta, \ phi) = m \ hbar Y _ {\ ell m} (\ theta, \ phi) \ egroup \ end {displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} L ^ 2 Y _ {\ ell m} (\ theta, \ phi) = \ ell (\ ell + 1) \ hbar ^ 2 Y _ {\ ell m} (\ theta , \ phi) \ egroup \ end {displaymath}$

com $\ bgroup \ color {black} $ \ ell $ \ egroup$ e $\ bgroup \ color {black} $ m $ \ egroup$ inteiros que satisfaçam a condição . Os operadores que aumentam e diminuem o componente do momento angular são $\ bgroup \ color {black} $ - \ ell \ leq m \ leq \ ell $ \ egroup$

$\ begin {displaymath} \ bgroup \ color {black} L_ \ pm = L_x \ pm iL_y \ egroup \ end {displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} L_ \ pm Y _ {\ ell m} = \ hbar \ sqrt {\ ell (\ ell + 1) -m (m \ pm 1)} Y _ {\ ell (m \ pm 1)} \ egroup \ end {displaymath}$

Derivamos a forma funcional dos Harmônicos esféricos usando a forma diferencial dos operadores de momento angular. $\ bgroup \ color {black} $ Y _ {\ ell m} (\ theta, \ phi) $ \ egroup$

Soluções para a equação radial para potenciais constantes

Soluções para a equação radial em potencial constante são importantes, pois são as soluções para grandes $\ bgroup \ color {black} $ r $ \ egroup$ potenciais de faixa limitada. Eles são, portanto, usados em problemas de dispersão como os estados de entrada e saída. As soluções são as funções esféricas de Bessel e as funções esféricas de Neumann.

$\ begin {displaymath} \ bgroup \ color {black} j_ \ ell (\ rho) = (- \ rho) ^ \ ell \ left ({1 \ over \ rh ... ... r \ rho} \ rightarrow { \ sin (\ rho - {\ ell \ pi \ over 2}) \ over \ rho} \ egroup \ end {displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} n_ \ ell (\ rho) = - (- \ rho) ^ \ ell \ left ({1 \ over \ r ... ... \ rho} \ rightarrow { - \ cos (\ rho - {\ ell \ pi \ over 2}) \ over \ rho} \ egroup \ end {displaymath}$

onde $\ bgroup \ color {black} $ \ rho = kr $ \ egroup$ . A combinação linear dessas que cai adequadamente em geral $\ bgroup \ color {black} $ r $ \ egroup$ é chamada de função de Hankel do primeiro tipo.

$\ begin {displaymath} \ bgroup \ color {black} h ^ {(1)} _ \ ell (\ rho) = j_ \ ell (\ rho) + i n_ \ ell (... ...} \ rightarrow - {i \ over \ rho} e ^ {i (\ rho - {\ ell \ pi \ over 2})} \ egroup \ end {displaymath}$

Usamos essas soluções para fazer uma análise de onda parcial de espalhamento , resolver para estados vinculados de um poço de potencial esférico , resolver para estados vinculados de um poço esférico infinito (uma "caixa'' esférica) e resolver para espalhamento de um potencial esférico.

Hidrogênio

A equação radial do hidrogênio (potencial de Coulomb) é resolvida encontrando o comportamento de $\ bgroup \ color {black} $ r $ \ egroup$ em geral, em seguida, encontrando o comportamento de $\ bgroup \ color {black} $ r $ \ egroup$ pequeno , em seguida usando uma solução de séries de potências para obter

$\ begin {displaymath} \ bgroup \ color {black} R (\ rho) = \ rho ^ \ ell \ sum \ limits_ {k = 0} ^ \ infty a_k \ rho ^ ke ^ {- \ rho / 2} \ egroup \ end {displaymath}$

com

$\ bgroup \ color {black} $ \ rho = \ sqrt {-8 \ mu E \ over \ hbar ^ 2} r $ \ egroup$

Para manter a função de onda normalizável, a série de potências deve terminar, dando-nos nossa condição de autovalor de energia.

$\ begin {displaymath} \ bgroup \ color {black} E_n = - {Z ^ 2 \ alpha ^ 2mc ^ 2 \ over 2n ^ 2} \ egroup \ end {displaymath}$

Aqui $\ bgroup \ color {black} $ n $ \ egroup$ é chamado de número quântico principal e é dado por

$\ begin {displaymath} \ bgroup \ color {black} n = n_r + \ ell + 1 \ egroup \ end {displaymath}$

onde $\ bgroup \ color {black} $ n_r $ \ egroup$ é o número de nós na função de onda radial. É uma característica estranha do hidrogênio que uma excitação radial e uma excitação angular tenham a mesma energia.

Portanto, um autoestado

$\ bgroup \ color {black} $ \ psi_ {n \ ell m} (\ vec {x}) = R_ {n \ ell} (r) Y _ {\ ell m} (\ theta, \ phi) $ \ egroup$

de energia de hidrogênio é descrito por três números quânticos inteiros com requisitos de

$\ bgroup \ color {black} $ n \ geq 1 $ \ egroup$ , $\ bgroup \ color {black} $ \ ell <n $ \ egroup$

e também um inteiro, e

$\ bgroup \ color {black} $ - l \ leq m \ leq \ ell $ \ egroup$ .

O estado fundamental do Hidrogênio é $\bgroup\color{black}$\psi_{100}$\egroup$ e tem energia de -13,6 eV. Calculamos vários dos autoestados de menor energia.

O diagrama abaixo mostra os estados de limite de energia mais baixos do hidrogênio e seus decaimentos típicos.

$\ epsfig {arquivo = figs / H_spect_simp.eps, altura = 4in}$

Solução do problema 3D do OH em coordenadas esféricas

Como exemplo de outro problema com simetria esférica, resolvemos o problema do oscilador harmônico simétrico 3D . Já resolvemos esse problema em coordenadas cartesianas. Agora usamos coordenadas esféricas e autofunções de momento angular.

As energias próprias são

$\ begin {displaymath} \ bgroup \ color {black} E = \ left (2n_r + \ ell + {3 \ over 2} \ right) \ hbar \ omega \ egroup \ end {displaymath}$

onde $\ bgroup \ color {black} $ n_r $ \ egroup$ é o número de nós na função de onda radial e $\ bgroup \ color {black} $ \ ell $ \ egroup$ é o número quântico do momento angular total. Isso dá exatamente o mesmo conjunto de energias próprias que obtivemos na solução cartesiana, mas os estados próprios são agora estados de momento angular total definido e componente z do momento angular.

Representação Matricial de Operadores e Estados

Podemos definir os componentes de um vetor de estado $\ bgroup \ color {black} $ \ psi $ \ egroup$ como as projeções do estado em um conjunto ortonormal completo de estados, como as autofunções de um operador Hermitiano.

$\ begin {eqnarray *} \ psi_i & \ equiv & \ langle u_i \ vert \ psi \ rangle \\ \ vert \ psi \ rangle & = & \ sum \ limits_i \ psi_i \ vert u_i \ rangle \ end {eqnarray *}$

Da mesma forma, podemos definir o elemento de matriz de um operador em termos de um par desses estados de base ortonormais

$\ begin {displaymath} \ bgroup \ color {black} O_ {ij} \ equiv \ langle u_i \ vert O \ vert u_j \ rangle. \ egroup \ end {displaymath}$

Com essas definições, problemas de Mecânica Quântica podem ser resolvidos usando os operadores e estados de representação de matrizes. Um operador agindo em um estado é uma matriz vezes um vetor.

$\ begin {displaymath} \ bgroup \ color {black} \ left (\ matrix {(O \ psi) _1 \ cr (O \ psi) _2 \ cr ... ... matrix {\ psi_1 \ cr \ psi_2 \ cr ... \ cr \ psi_j \ cr ...} \ right) \ egroup \ end {displaymath}$

O produto dos operadores é o produto das matrizes. Operadores que não comutam são representados por matrizes que não comutam.

Um estudo de operadores $\ bgroup \ color {black} $ \ ell = 1 $ \ egroup$ e autofunções

O conjunto de estados com o mesmo momento angular total e os operadores de momento angular que atuam sobre eles são frequentemente representados por vetores e matrizes. Por exemplo, os diferentes $\ bgroup \ color {black} $ m $ \ egroup$ estados para $\ bgroup \ color {black} $ \ ell = 1 $ \ egroup$ serão representados por um vetor de 3 componentes e os operadores de momento angular representados por matrizes 3X3.

Existem razões práticas e teóricas pelas quais esse conjunto de estados é separado dos estados com diferentes números quânticos de momento angular total.

Os estados são frequentemente (quase) degenerados e, portanto, devem ser tratados como um grupo por razões práticas.

Além disso, uma rotação dos eixos coordenados não mudará o número quântico total do momento angular, portanto, o operador de rotação trabalha dentro deste grupo de estados.

Escrevemos nossos 3 vetores de componentes da seguinte maneira.

$\ begin {displaymath} \ bgroup \ color {black} \ psi = \ left (\ matrix {\ psi_ + \ cr \ psi_0 \ cr \ psi _-} \ right) \ egroup \ end {displaymath}$

As matrizes que representam os operadores de momento angular para $\ bgroup \ color {black} $ \ ell = 1 $ \ egroup$ são as seguintes

$\ begin {eqnarray *} L_x = {\ hbar \ over \ sqrt {2}} \ left (\ matrix {0 & 1 & 0 \ cr 1 & 0 & 1 \ cr0 & 1 & 0} \ ... ... d \ qquad L_z = \ hbar \ left (\ matriz {1 & 0 & 0 \ cr 0 & 0 & 0 \ cr 0 & 0 & -1} \ right) \ end {eqnarray *}$

As mesmas matrizes também representam o spin 1, $\ bgroup \ color {black} $ s = 1 $ \ egroup$ mas é claro que atuariam em um espaço vetorial diferente.

Os operadores de rotação (operadores de simetria) são dados por

$\ begin {eqnarray *} R_z (\ theta_z) = e ^ {i \ theta_zL_z / \ hbar} \ qquad \ qquad & R_x (\ theta _... ... a_xL_x / \ hbar} \ qquad \ qquad & R_y (\ theta_y) = e ^ {i \ theta_yL_y / \ hbar} \ end {eqnarray *}$

para a forma diferencial ou a forma matricial dos operadores.

Para $\ bgroup \ color {black} $ \ ell = 1 $ \ egroup$ essas são matrizes 3X3 (unitárias). Nós as usamos quando precisamos redefinir a direção de nossos eixos coordenados. As rotações dos estados do momento angular não são iguais às rotações dos vetores no espaço 3. Os componentes dos vetores representam quantidades diferentes e, portanto, se transformam de maneira bastante diferente. Os "vetores'' que estamos usando para o momento angular na verdade deveriam ser chamados de espinores quando nos referimos às suas propriedades sob rotações e impulsos de Lorentz.

Spin 1/2 e outros 2 sistemas de estado

A álgebra do momento angular definida pelas relações de comutação entre os operadores requer que o número quântico do momento angular total seja um inteiro ou meio inteiro. A possibilidade de meio inteiro não era útil para o momento angular orbital porque não havia função harmônica esférica correspondente (valor único) para representar a amplitude para uma partícula estar em alguma posição.

A possibilidade de meio inteiro é usada para representar o momento angular interno de algumas partículas. O caso mais simples e importante é girar pela metade. Existem apenas dois estados possíveis com diferentes componentes z de spin: spin para cima $\ bgroup \ color {black} $ \ pmatrix {1 \ cr 0} $ \ egroup$ , com componente z do momento angular $\ bgroup \ color {black} $ + {\ hbar \ over 2} $ \ egroup$ , e spin para baixo $\ bgroup \ color {black} $ \ pmatrix {0 \ cr 1} $ \ egroup$ , com $\ bgroup \ color {black} $ - {\ hbar \ over 2} $ \ egroup$ . Os operadores de rotação correspondentes são

$\ begin {eqnarray *} S_x = {\ hbar \ over 2} \ left (\ matrix {0 & 1 \ cr 1 & 0 \ cr} \ right) \ qqu ... ... qquad S_z = {\ hbar \ over 2 } \ left (\ matrix {1 & 0 \ cr 0 & -1 \ cr} \ right) \ end {eqnarray *}$

que satisfazem as relações de comutação usuais das quais derivamos as propriedades dos operadores de momento angular.

É comum definir as Matrizes de Pauli $\ bgroup \ color {black} $ \ sigma_i $ \ egroup$ que possuem as seguintes propriedades.

$\ begin {eqnarray *} S_i & \ equiv & {\ hbar \ over 2} \ sigma_i. \\ \ vec {S} & = & {\ hbar \ ov ... ... ma_y + \ sigma_y \ sigma_z = 0 \\ \ {\ sigma_i, \ sigma_j \} & = & 2 \ delta_ {ij} \ end { eqnarray *}$

As duas últimas linhas afirmam que as matrizes de Pauli são anti-comutação. As matrizes $\ bgroup \ color {black} $ \ sigma $ \ egroup$ são as matrizes Hermitianas e sem rastros de dimensão 2. Qualquer matriz 2 por 2 pode ser escrita como uma combinação linear das matrizes $\ bgroup \ color {black} $ \ sigma $ \ egroup$ e da identidade.

Mecânica Quântica em um Campo Eletromagnético

O hamiltoniano clássico para uma partícula em um campo eletromagnético é

$\ begin {displaymath} \ bgroup \ color {black} H = {1 \ over 2m} \ left (\ vec {p} + {e \ over c} \ vec {A} \ right) ^ 2-e \ phi \ egroup \ end {displaymath}$

onde $\ bgroup \ color {black} $ e $ \ egroup$ é definido como um número positivo. Este hamiltoniano fornece a lei de força de Lorentz correta. Observe que o operador de momento agora incluirá o momento no campo, não apenas o momento da partícula. Como este hamiltoniano é escrito, $\ bgroup \ color {black} $ \ vec {p} $ \ egroup$ é a variável conjugada a $\ bgroup \ color {black} $ \ vec {r} $ \ egroup$ e está relacionada à velocidade de

$\ bgroup \ color {black} $ \ vec {p} = m \ vec {v} - {e \ over c} \ vec {A} $ \ egroup$ .

Na Mecânica Quântica, o operador momentum é substituído da mesma maneira para incluir os efeitos dos campos magnéticos e, eventualmente, da radiação.

$\ begin {displaymath} \ bgroup \ color {preto} \ vec {p} \ rightarrow \ vec {p} + {e \ over c} \ vec {A} \ egroup \ end {displaymath}$

Partindo do hamiltoniano acima, derivamos o hamiltoniano para uma partícula em um campo magnético constante .

$\ begin {displaymath} \ bgroup \ color {black} {- \ hbar ^ 2 \ over 2m} \ nabla ^ 2 \ psi + {e \ over 2m ... ... (r ^ 2B ^ 2 - (\ vec { r} \ cdot \ vec {B}) ^ 2 \ right) \ psi = (E + e \ phi) \ psi \ egroup \ end {displaymath}$

Isso tem o efeito familiar de um momento magnético paralelo ao vetor do momento angular, além de alguns termos adicionais que são muito pequenos para átomos em campos realizáveis em laboratório.

Então, para átomos, o termo adicional dominante é

$\ begin {displaymath} \ bgroup \ color {black} H_B = {e \ over 2mc} \ vec {B} \ cdot \ vec {L} = - \ vec {\ mu} \ cdot \ vec {B}, \ egroup \ end {displaymath}$

onde

$\ bgroup \ color {black} $ \ vec {\ mu} = - {e \ over 2mc} \ vec {L} $ \ egroup$ .

Este é, efetivamente, o momento magnético devido ao momento angular orbital do elétron.

Os outros termos podem ser importantes se um estado está espalhado por uma região muito maior do que um átomo.

Trabalhamos o exemplo de um plasma em um campo magnético constante. Uma partícula carregada no plasma tem o seguinte espectro de energia

$\ begin {displaymath} \ bgroup \ color {black} E_n = {eB \ hbar \ over m_ec} \ left (n + {1 \ over 2} \ right) + {\ hbar ^ 2k ^ 2 \ over 2m_e}. \ egroup \ end {displaymath}$

que depende de 2 números quânticos.

$\ bgroup \ color {black} $ \ hbar k $ \ egroup$ é o momento conservado ao longo da direção do campo que pode assumir qualquer valor.

$\ bgroup \ color {black} $ n $ \ egroup$ é um número inteiro que lida com o estado em x e y.

Este problema pode ser simplificado usando-se alguns operadores de simetria diferentes. Nós o trabalhamos de duas maneiras diferentes: em uma se reduz à equação radial para o átomo de Hidrogênio; no outro, reduz-se à equação do oscilador harmônico, mostrando que esses dois problemas que podemos resolver são de alguma forma equivalentes.

Simetria de fase local em mecânica quântica e a simetria de calibre

Existe uma simetria na física que podemos chamar de Simetria de Fase Local na mecânica quântica. Nessa simetria, mudamos a fase da função de onda (do elétron) em uma quantidade diferente em todo o espaço-tempo. Para compensar essa mudança, precisamos também fazer uma transformação de calibre dos potenciais eletromagnéticos. Todos eles devem ir juntos assim.

$\ begin {eqnarray *} \ psi (\ vec {r}, t) & \ rightarrow & e ^ {- i {e \ over \ hbar c} f (\ vec {r}, ... ... ghtarrow & \ phi + {1 \ over c} {\ partial f (\ vec {r}, t) \ over \ partial t} \ end {eqnarray *}$

A simetria de fase local requer que o eletromagnetismo exista e tenha uma simetria de calibre para que possamos manter a equação de Schrödinger invariante sob esta transformação de fase.

Exploramos a simetria do calibre em EM para mostrar que, em regiões livres de campo , a função $\ bgroup \ color {black} $ f $ \ egroup$ pode ser simplesmente igual a uma integral de linha do potencial vetorial (se escolhermos o calibre correto).

$\ begin {displaymath} \ bgroup \ color {black} f (\ vec {r}) = \ int \ limits _ {\ vec {r} _0} ^ {\ vec {r}} d \ vec {r} \ cdot \ vec {A}. \ egroup \ end {displaymath}$

Usamos isso para mostrar que o fluxo magnético encerrado por um supercondutor é quantizado.

Também mostramos que os campos magnéticos podem ser usados para alterar os efeitos de interferência na mecânica quântica. O Efeito Aharanov Böhm nos traz de volta ao experimento de difração de duas fendas, mas adiciona campos magnéticos.

$\ epsfig {arquivo = figs / aharanov.eps, altura = 3,5 pol.}$

Os feixes de elétrons viajam por duas fendas em regiões livres de campo, mas temos a capacidade de variar um campo magnético delimitado pelo caminho dos elétrons. Na tela, as amplitudes das duas fendas interferem . Vamos começar com e em todos os lugares. Quando mudamos o campo, as funções de onda devem mudar. $\ bgroup \ color {black} $ \ psi = \ psi_1 + \ psi_2 $ \ egroup$ $\ bgroup \ color {black} $ B = 0 $ \ egroup$ $\ bgroup \ color {black} $ A = 0 $ \ egroup$ $\ bgroup \ color {black} $ B $ \ egroup$

$\ begin {eqnarray *} \ psi_1 & \ rightarrow & \ psi_1 e ^ {- i {e \ over \ hbar c} \ int \ limits ... ... 2 \ right) e ^ {- i {e \ over \ hbar c} \ int \ limits_ {2} d \ vec {r} \ cdot \ vec {A}} \ end {eqnarray *}$

A fase relativa das duas fendas depende do fluxo entre as fendas. Variando o $\ bgroup \ color {black} $ B $ \ egroup$ campo, mudaremos o padrão de difração , embora $\ bgroup \ color {black} $ B = 0 $ \ egroup$ ao longo de todo o caminho dos elétrons.

Adição de Momentum Angular

Freqüentemente, é necessário adicionar o momento angular de duas (ou mais) fontes para obter estados de momento angular total definido. Por exemplo, na ausência de campos externos, os estados próprios de energia do Hidrogênio (incluindo todos os efeitos de estrutura fina) também são estados próprios do momento angular total . Isso quase tem que ser verdade se houver simetria esférica para o problema.

Como exemplo, vamos supor que estamos adicionando o momento angular orbital de dois elétrons $\ bgroup \ color {black} $ \ vec {L_1} $ \ egroup$ e $\ bgroup \ color {black} $ \ vec {L_2} $ \ egroup$ para obter um momento angular total $\ bgroup \ color {black} $ \ vec {J} $ \ egroup$ . Mostraremos que o número quântico do momento angular total assume todos os valores na faixa

$\ begin {displaymath} \ bgroup \ color {preto} \ vert \ ell_1- \ ell_2 \ vert \ leq j \ leq \ ell_1 + \ ell_2. \ egroup \ end {displaymath}$

Podemos entender isso qualitativamente no modelo vetorial ilustrado abaixo. Estamos adicionando dois vetores quânticos.

$\ epsfig {arquivo = figs / vector.eps, altura = 2in}$

O comprimento do vetor resultante está em algum lugar entre a diferença de suas magnitudes e a soma de suas magnitudes, uma vez que não sabemos para qual direção os vetores estão apontando.

Os estados de momento angular total definido com números quânticos $\ bgroup \ color {black} $ j $ \ egroup$ e $\ bgroup \ color {black} $ m $ \ egroup$ , podem ser escritos em termos de produtos dos estados individuais (como o elétron 1 está neste estado E o elétron 2 está nesse estado). A expansão geral é chamada de série Clebsch-Gordan:

$\ begin {displaymath} \ bgroup \ color {black} \ psi_ {jm} = \ sum \ limits_ {m_1m_2} \ langle \ el ... ... 2 \ vert jm \ ell_1 \ ell_2 \ rangle Y _ {\ ell_1m_1} Y _ {\ ell_2m_2} \ egroup \ end {displaymath}$

ou em termos dos vetores ket

$\ begin {displaymath} \ bgroup \ color {black} \ vert jm \ ell_1 \ ell_2 \ rangle = \ sum \ limits _ {... ... rt jm \ ell_1 \ ell_2 \ rangle \ vert \ ell_1m_1 \ ell_2m_2 \ rangle \ egroup \ end {displaymath}$

Os coeficientes de Clebsch-Gordan são tabulados, embora possamos calcular muitos deles nós mesmos.

Ao combinar estados de partículas idênticas, o estado de momento angular total mais alto,, $\ bgroup \ color {black} $ s = s_1 + s_2 $ \ egroup$ será sempre simétrico sob intercâmbio. A simetria sob o intercâmbio irá alternar conforme $\ bgroup \ color {black} $ j $ \ egroup$ é reduzida.

O número total de estados é sempre preservado. Por exemplo, se eu adicionar dois $\ bgroup \ color {black} $ \ ell = 2 $ \ egroup$ estados, obtenho estados de momento angular total com $\ bgroup \ color {black} $ j = 0,1,2,3 $ \ egroup$ e 4. Existem 25 estados de produto, pois cada $\ bgroup \ color {black} $ \ ell = 2 $ \ egroup$ estado tem 5 diferentes $\ bgroup \ color {black} $ m $ \ egroup$ s possíveis . Compare isso com a soma do número de estados que acabamos de listar.

$\ begin {displaymath} \ bgroup \ color {black} 5 \ otimes 5 = 9_S \ oplus 7_A \ oplus 5_S \ oplus 3_A \ oplus 1_S \ egroup \ end {displaymath}$

onde os números são o número de estados no multipleto.

Usaremos adição de momento angular para:

Adicione o momento angular orbital ao momento angular de rotação de um elétron em um átomo ; $\ vec {J} = \ vec {L} + \ vec {S}$
Some os momentos angulares orbitais de dois elétrons em um átomo ; $\ vec {L} = \ vec {L_1} + \ vec {L_2}$
Some os spins de duas partículas ; $\ vec {S} = \ vec {S_1} + \ vec {S_2}$
Adicione o spin nuclear ao momento angular atômico total ; $\ vec {F} = \ vec {J} + \ vec {I}$
Some os momentos angulares totais de dois elétrons ; $\ vec {J} = \ vec {J_1} + \ vec {J_2}$
Some o momento angular orbital total ao momento angular de spin total para uma coleção de elétrons em um átomo ; $\ vec {J} = \ vec {L} + \ vec {S}$
Escreva o produto dos harmônicos esféricos em termos de uma soma de harmônicos esféricos.

Teoria da Perturbação Independente do Tempo

Suponha que já resolvemos um problema de autovalor de energia e agora precisamos incluir um termo adicional no hamiltoniano. Podemos usar a teoria de perturbação independente do tempo para calcular correções para os autovalores e estados próprios da energia. Se a equação de Schrödinger para o problema completo for

$\ begin {displaymath} \ bgroup \ color {black} (H_0 + H_1) \ psi_n = E_n \ psi_n \ egroup \ end {displaymath}$

e já resolvemos o problema dos autovalores $\ bgroup \ color {black} $ H_0 $ \ egroup$ , pois podemos usar uma série de perturbações, para expandir nossos autovalores de energia e autossados em potências da pequena perturbação.

$\ begin {eqnarray *} & E_n = E_n ^ {(0)} + E_n ^ {(1)} + E_n ^ {(2)} + ... \\ & \ psi_n = N \ left (\ p .. . ... n} c_ {nk} \ phi_k \ right) \\ & c_ {nk} = c_ {nk} ^ {(1)} + c_ {nk} ^ {(2)} + ... \ end {eqnarray *}$

onde o sobrescrito (0), (1), (2) são os termos zero, de primeira e de segunda ordem na série. $\ bgroup \ color {black} $ N $ \ egroup$ existe para manter a função de onda normalizada, mas não terá um papel importante em nossos resultados.

Resolvendo a equação de Schrödinger em cada ordem da série de perturbações, calculamos as correções para as energias e autofunções. Fornecemos apenas os primeiros termos acima.

$\ begin {eqnarray *} & E_n ^ {(1)} = \ langle \ phi_n \ vert H_1 \ vert \ phi_n \ rangle \\ & c _ {... ..._ k \ vert H_1 \ vert \ phi_n \ rangle \ vert ^ 2 \ over E_n ^ {(0)} - E_k ^ {(0)}} \ end {eqnarray *}$

Um problema surge no caso de estados degenerados ou quase degenerados. O denominador de energia na última equação acima é pequeno e a série não converge. Para lidar com este caso, precisamos rediagonalizar o hamiltoniano completo no subespaço de estados quase degenerados.

$\ begin {displaymath} \ bgroup \ color {black} \ sum \ limits_ {i \ in {\ cal N}} \ langle \ phi_n ^ {(j)} \ vert H \ vert \ phi_n ^ {(i)} \ rangle \ alpha_i = E_n \ alpha_j. \ egroup \ end {displaymath}$

Este é apenas o problema de autovalor padrão para o hamiltoniano completo no subespaço de estados (quase) degenerados.

Usaremos a teoria de perturbação independente do tempo para calcular estruturas finas e correções hiperfinas para energias de hidrogênio, bem como para muitos outros cálculos. A teoria de perturbação de estado degenerado será usada para o Efeito Stark e para a divisão hiperfina em Hidrogênio.

A Estrutura Fina do Hidrogênio

Resolvemos o problema de um elétron não relativístico e sem spin em um potencial coulomb exatamente. Os átomos de hidrogênio reais têm várias pequenas correções para essa solução simples. Se dissermos que o spin do elétron é um efeito relativístico, todas elas podem ser chamadas de correções relativísticas, que estão fora de ordem em $\ bgroup \ color {black} $ \ alpha ^ 2 $ \ egroup$ comparação com as energias de hidrogênio que calculamos.

A correção relativística da energia cinética do elétron.
A correção Spin-Orbit.
A correção do `` Termo de Darwin '' para os estados da equação de Dirac.

Calculando esses efeitos de estrutura fina separadamente e somando-os, descobrimos que obtemos um bom cancelamento que resulta em uma fórmula simples.

$\ begin {displaymath} \ bgroup \ color {black} E_ {nlm} = E ^ {(0)} _ n + {{E ^ {(0)} _ n} ^ 2 \ over {2mc ^ 2}} \ left [ 3 - {4n \ over {j + {1 \ over 2}}} \ right] \ egroup \ end {displaymath}$

A correção depende apenas do número quântico angular total e não depende $\ bgroup \ color {black} $ \ ell $ \ egroup$ dos estados de diferentes momentos angulares totais divididos em energia, mas ainda há uma boa dose de degenerescência. Faz sentido, para um problema com simetria esférica, que os estados de momento angular total definido sejam os autoestados de energia e dos quais o resultado dependa $\ bgroup \ color {black} $ j $ \ egroup$ .

Também calculamos o efeito Zeeman no qual um campo magnético externo é aplicado ao hidrogênio. O campo externo é muito importante, pois quebra a simetria esférica e divide os estados degenerados, permitindo entender o Hidrogênio por meio da espectroscopia.

A correção devido a um campo magnético fraco é encontrada para ser

$\ begin {displaymath} \ bgroup \ color {black} \ Delta E = \ left <\ psi_ {n \ ell jm_j} \ left \ ve ... ... r B \ over {2mc}} m_j \ left (1 \ pm {1 \ over {2 \ ell +1}} \ right) \ egroup \ end {displaymath}$

O fator é conhecido como Fator Lande porque o estado se divide como se tivesse essa razão giromagnética. Sabemos que é de fato uma combinação dos fatores orbitais e de spin g em um estado de definido . Assumimos que o efeito do campo é pequeno em comparação com as correções de estrutura fina. Podemos escrever a energia total em um campo magnético fraco. $\ bgroup \ color {black} $ \ left (1 \ pm {1 \ over {2 \ ell +1}} \ right) $ \ egroup$ $\ bgroup \ color {black} $ j $ \ egroup$

$\ begin {displaymath} \ bgroup \ color {black} E_ {njm_j \ ell s} = - {1 \ over 2} \ alpha ^ 2mc ^ 2 \ ... ... er j + {1 \ over 2}} - {3 \ over 4n} \ right] \ right) + g_L \ mu_BBm_j \ egroup \ end {displaymath}$

Assim, em um campo fraco, a degeneração é completamente quebrada para os estados $\ psi_ {njm_j \ ell s}$ . Todos os estados podem ser detectados espectroscopicamente.

No limite forte campo poderíamos usar estados de definitiva $\ bgroup \ color {black} $ m_ \ ell $ \ egroup$ e $\ bgroup \ color {black} $ m_s $ \ egroup$ e calcular os efeitos da estrutura fina, $\ bgroup \ color {black} $ H_1 + H_2 $ \ egroup$ como uma perturbação. Em um campo de força intermediária, da ordem de 500 Gauss, a combinação da Hamiltoniana da estrutura fina do Hidrogênio e o termo para o campo B deve ser diagonalizado no conjunto de estados com o mesmo número quântico principal $\ bgroup \ color {black} $ n $ \ egroup$ .

Estrutura Hiperfina

A interação entre o spin do núcleo e o momento angular do elétron causa uma divisão adicional (hiperfina) dos estados atômicos. É chamado de hiperfino porque também é ordenado $\ bgroup \ color {black} $ \ alpha ^ 2 $ \ egroup$ como as correções de estrutura fina, mas é menor por um fator de aproximadamente por causa da dependência de massa do momento magnético de rotação para as partículas. $\ bgroup \ color {black} $ {m_e \ over m_p} $ \ egroup$

O momento magnético do núcleo é

$\ begin {displaymath} \ bgroup \ color {black} \ vec \ mu_N = {Zeg_N \ over {2M_Nc}} \ vec I \ egroup \ end {displaymath}$

onde $\ bgroup \ color {black} $ \ vec {I} $ \ egroup$ está o vetor de spin nuclear . Como o núcleo, o próton e o nêutron têm estrutura interna, a razão giromagnética nuclear não é apenas 2. Para o próton, é . $\ bgroup \ color {black} $ g_p \ approx 5,56 $ \ egroup$

Calculamos a contribuição hiperfina para o hamiltoniano para $\ bgroup \ color {black} $ \ ell = 0 $ \ egroup$ estados.

$\ begin {displaymath} \ bgroup \ color {black} H_ {hf} = \ left <{e \ over {mc}} \ vec S \ cdot \ ve ... ...) g_N {1 \ over {n ^ 3}} {\ vec {S} \ cdot \ vec {I} \ over {\ hbar ^ 2}} \ egroup \ end {displaymath}$

Agora, assim como no caso da interação spin-órbita, vamos definir o momento angular total $\ bgroup \ color {black} $ \ vec {L} \ cdot \ vec {S} $ \ egroup$

$\ begin {displaymath} \ bgroup \ color {preto} \ vec {F} = \ vec {S} + \ vec {I}. \ egroup \ end {displaymath}$

Está nos estados de definido $\ bgroup \ color {black} $ f $ \ egroup$ e $\ bgroup \ color {black} $ m_f $ \ egroup$ que a perturbação hiperfina será diagonal. Em essência, estamos fazendo a teoria da perturbação de estado degenerada. Poderíamos diagonalizar a matriz 4 por 4 para a perturbação para resolver o problema ou podemos usar o que sabemos para escolher os estados corretos para começar. Novamente, como a interação da órbita do spin, os estados de momento angular total serão os estados corretos porque podemos escrever a perturbação em termos de números quânticos desses estados.

$\ begin {displaymath} \ bgroup \ color {black} \ vec {S} \ cdot \ vec {I} = {1 \ over 2} \ left (F ^ 2 -... ... er 2} \ hbar ^ 2 \ left (f (f + 1) - {3 \ over 4} - {3 \ over 4} \ right) \ egroup \ end {displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} \ Delta E = {2 \ over 3} (Z \ alpha) ^ 4 \ left (... ... mc ^ 2) g_N {1 \ over {n ^ 3}} \ left (f (f + 1) - {3 \ over 2} \ right). \ egroup \ end {displaymath}$

Para o estado fundamental do hidrogênio, estamos apenas adicionando duas $\ bgroup \ color {black} $ {1 \ over 2} $ \ egroup$ partículas de spin , então os valores possíveis são $\ bgroup \ color {black} $ f = 0,1 $ \ egroup$ . A transição entre os dois estados dá origem a ondas EM com $\ bgroup \ color {black} $ \ lambda = 21 $ \ egroup$ cm.

Vamos trabalhar o efeito de um campo B externo nos estados hiperfinos do Hidrogênio, tanto no campo forte quanto na aproximação do campo fraco. Também resolvemos o problema sem uma aproximação de intensidade de campo. O resultado da intensidade do campo intermediário sempre aplicável é que os quatro estados têm energias que dependem da intensidade do campo B. Dois dos estados próprios de energia se misturam de uma forma que também depende de B. As quatro energias são

$\ begin {eqnarray *} E & = & E_ {n00} + {{\ cal A} \ hbar ^ 2 \ over 4} \ pm \ mu_BB \\ E & = & E_ {n00} -... ... eft ({{ \ cal A} \ hbar ^ 2 \ over 2} \ right) ^ 2 + \ left (\ mu_BB \ right) ^ 2}. \ end {eqnarray *}$

O átomo de hélio

O hamiltoniano para o hélio tem os mesmos termos do hidrogênio, mas tem uma grande perturbação devido à repulsão entre os dois elétrons.

$\ begin {displaymath} \ bgroup \ color {black} H = {p ^ 2_1 \ over {2m}} + {p ^ 2_2 \ over {2m}} ... ... {r_2}} + {e ^ 2 \ over {\ vert \ vec {r} _1 - \ vec {r} _2 \ vert}} \ egroup \ end {displaymath}$

Observe que a perturbação devido à repulsão entre os dois elétrons é aproximadamente do mesmo tamanho que o resto do hamiltoniano, então a teoria de perturbação de primeira ordem é improvável que seja precisa.

O estado fundamental do Hélio tem dois elétrons no nível 1s . Como o estado espacial é simétrico, a parte do spin do estado deve ser anti-simétrica $\ bgroup \ color {black} $ s = 0 $ \ egroup$ (como sempre é para cascas fechadas). Para nossos autoestados de energia de ordem zero, usaremos estados de produto de funções de onda de hidrogênio

$\ begin {displaymath} \ bgroup \ color {black} u (\ vec {r} _1, \ vec {r} _2) = \ phi_ {n_1 \ ell_1 m_1} (\ vec {r} _1) \ phi_ {n_2 \ ell_2 m_2} (\ vec {r} _2) \ egroup \ end {displaymath}$

e ignorar a perturbação. A energia para dois elétrons no estado (1s) para $\ bgroup \ color {black} $ Z = 2 $ \ egroup$ é então eV. $\ bgroup \ color {black} $ 4 \ alpha ^ 2mc ^ 2 = 108,8 $ \ egroup$

Podemos estimar a energia do estado fundamental na teoria de perturbação de primeira ordem , usando o termo de repulsão de elétrons como uma perturbação (muito grande). Isso não é muito preciso.

Podemos melhorar a estimativa da energia do estado fundamental usando o princípio variacional . O principal problema com nossa estimativa da teoria de perturbação é que não estamos contabilizando as mudanças na função de onda dos elétrons devido à triagem . Podemos fazer isso em alguma aproximação razoável, reduzindo a carga do núcleo na função de onda (não no hamiltoniano). Com o parâmetro $\ bgroup \ color {black} $ Z ^ * $ \ egroup$ , obtemos uma melhor estimativa da energia.

Cálculo	Energia	$\ bgroup \ color {black} $ Z_ {wfn} $ \ egroup$
$\ bgroup \ color {black} $ 0 ^ {th} $ \ egroup$ Pedido	-108,8	2
$\ bgroup \ color {black} $ 1 ^ {st} $ \ egroup$ Teoria de perturbação de ordem	-74,8	2
$\ bgroup \ color {black} $ 1 ^ {st} $ \ egroup$ Ordem Variacional	-77,38	$\ bgroup \ color {black} $ {27 \ over 16} $ \ egroup$
Real	-78.975

Observe que o cálculo variacional ainda usa a teoria de perturbação de primeira ordem. Ele apenas adiciona um parâmetro variável à função de onda que usamos para minimizar a energia. Isso só funciona para o estado fundamental e para outros estados especiais.

Existe apenas um $\ bgroup \ color {black} $ (1s) ^ 2 $ \ egroup$ estado permitido e é o estado fundamental. Para estados excitados , os estados espaciais são (geralmente) diferentes, então eles podem ser simétricos ou anti-simétricos (sob o intercâmbio dos dois elétrons). Acontece que o estado anti-simétrico mantém os elétrons mais separados, então a repulsão é menor e a energia é menor. Se o estado espacial for antissimétrico, então o estado de spin é simétrico, s = 1. Portanto, os estados tripletos são geralmente significativamente mais baixos em energia do que os estados singuletos de spin correspondentes. Esta parece ser uma forte interação dependente do spin, mas na verdade é apenas o efeito da repulsão entre os elétrons, tendo um grande efeito dependendo da simetria do estado espacial e, portanto, da simetria do estado de spin.

O primeiro estado de saída tem o conteúdo do estado hidrogenico de (1s) (2s) e tem s = 1. Calculamos a energia deste estado.

Aprenderemos mais tarde que as transições eletromagnéticas que mudam o spin são fortemente suprimidas, fazendo com que o tripleto de spin (ortohelium) e os estados de spin singlet (parahelium) tenham cadeias de decaimento quase separadas.

Física Atômica

O hamiltoniano para um átomo com elétrons e prótons Z tem muitos termos que representam a repulsão entre cada par de elétrons.

$\ begin {displaymath} \ bgroup \ color {black} \ left [\ sum \ limits ^ Z_ {i = 1} \ left ({p ^ 2_i \ o ... ... \ vert \ vec {r_i} - \ vec {r_j} \ right \ vert}} \ right] \ psi = E \ psi. \ egroup \ end {displaymath}$

Vimos que a repulsão coulombiana entre elétrons é uma correção muito grande no Hélio e que o problema dos três corpos na mecânica quântica só é resolvido por aproximação.

A física das conchas fechadas e do momento angular nos permite entender até mesmo os átomos mais complexos. Quando temos elétrons suficientes para preencher uma camada, digamos 1s ou 2p, a distribuição de elétrons resultante é esfericamente simétrica porque

$\ begin {displaymath} \ bgroup \ color {black} \ sum ^ \ ell_ {m = - \ ell} \ left \ vert Y _ {\ ell m} \ ... ... a, \ phi \ right) \ right \ vert ^ 2 = {2 \ ell + 1 \ over {4 \ pi}}. \ egroup \ end {displaymath}$

Com todos os estados preenchidos e as fases relativas determinadas pela antissimetria exigida por Pauli, os números quânticos da casca fechada são determinados. Existe apenas um estado possível que representa uma casca fechada e os números quânticos são

$\ begin {eqnarray *} & s = 0 \\ & \ ell = 0 \\ & j = 0 \ end {eqnarray *}$

A cápsula fechada protege a carga nuclear. Por causa da triagem , o potencial não tem mais um $\ bgroup \ color {black} $ {1 \ over r} $ \ egroup$ comportamento puro . Os elétrons que estão distantes do núcleo percebem menos a carga nuclear e aumentam sua energia. Vemos que as camadas atômicas são preenchidas na ordem 1s, 2s, 2p, 3s, 3p, 4s, 3d, 4p, 5s, 4d, 5p, 6s, 4f, 5d, 6p. O efeito da triagem aumentando a energia dos $\ bgroup \ color {black} $ \ ell $ \ egroup$ estados superiores é claro. Não é de admirar que a tabela periódica não seja totalmente periódica.

Um conjunto de diretrizes, conhecido como regras de Hund, nos ajuda a determinar os números quânticos para os estados fundamentais dos átomos. As conchas hidrogenadas se enchem dando $\ bgroup \ color {black} $ j = 0 $ \ egroup$ estados bem definidos para as conchas fechadas. À medida que adicionamos elétrons de valência , seguimos as regras de Hund para determinar o estado fundamental. Conseguimos uma grande simplificação tratando cascas quase fechadas como cascas fechadas mais buracos de rotação com carga positiva . Por exemplo, se um átomo tem dois elétrons aquém de uma camada fechada, nós o tratamos como uma camada fechada mais dois buracos positivos.) $\ bgroup \ color {black} $ {1 \ over 2} $ \ egroup$

Acople os elétrons de valência (ou lacunas) para dar o spin total máximo .
Agora escolha o estado de máximo $\ ell$ (sujeito ao princípio de Pauli. O princípio de Pauli em vez da regra, muitas vezes determina tudo aqui).
Se a casca estiver mais da metade cheia, escolha o estado de momento angular total mais alto, $j = \ ell + s$ caso contrário, escolha o mais baixo $j = \ vert \ ell-s \ vert$ .

Moléculas

Podemos estudar moléculas simples para entender os fenômenos físicos das moléculas em geral. A molécula mais simples com a qual podemos trabalhar é o íon $\ mathrm {H} _2 ^ +$ . Possui dois núcleos (A e B) compartilhando um elétron (1).

$\ begin {displaymath} \ bgroup \ color {black} H_0 = {p ^ 2_e \ over {2m}} - {e ^ 2 \ over {r_ {1A}}} - {e ^ 2 \ over {r_ {1B} }} + {e ^ 2 \ over {R_ {AB}}} \ egroup \ end {displaymath}$

$\ bgroup \ color {black} $ R_ {AB} $ \ egroup$ é a distância entre os dois núcleos. Calculamos a energia do estado fundamental usando os estados do hidrogênio como base.

A função de onda de energia mais baixa pode ser considerada uma combinação linear (anti) simétrica de um elétron no estado fundamental próximo ao núcleo A e o estado fundamental próximo ao núcleo B

$\ begin {displaymath} \ bgroup \ color {black} \ psi_ \ pm \ left (\ vec {r}, \ vec {R} \ right) = C_ \ pm (R) \ left [\ psi_A \ pm \ psi_B \ direita] \ egroup \ end {displaymath}$

onde é gs ao redor do núcleo A. e não são ortogonais; há sobreposição. O estado simétrico (ligação) tem uma grande probabilidade de o elétron ser encontrado entre os núcleos . O estado antissimétrico (anticoagulante) tem uma pequena probabilidade e, portanto, uma energia muito maior. Lembre-se de que essa simetria é a da função de onda de um elétron ao redor dos dois núcleos. $\ bgroup \ color {black} $ \ psi_A = \ sqrt {1 \ over {\ pi a ^ 3_0}} e ^ {- r _ {\ scriptscriptstyle 1A} / a_0} $ \ egroup$ $\ bgroup \ color {black} $ \ psi_A $ \ egroup$ $\ bgroup \ color {black} $ \ psi_B $ \ egroup$

A $\ mathrm {H} _2$ molécula também é simples e sua energia pode ser calculada com a ajuda do cálculo anterior. O estado simétrico do espaço será o estado fundamental.

$\ begin {displaymath} \ bgroup \ color {black} \ left <\ psi \ vert H \ vert \ psi \ right> = 2E_ {H _... ... si \ left \ vert {e ^ 2 \ over {r_ {12}}} \ right \ vert \ psi \ right> \ egroup \ end {displaymath}$

A molécula pode vibrar no potencial criado quando o elétron compartilhado une os átomos, dando origem a um espectro de energia do oscilador harmônico.

As moléculas podem girar como corpos rígidos clássicos sujeitos à restrição de que o momento angular é quantizado em unidades de $\ bgroup \ color {black} $ \ hbar $ \ egroup$ .

$\ begin {displaymath} \ bgroup \ color {black} E_ {rot} = {1 \ over 2} {L ^ 2 \ over I} = {\ ell (\ ell ... ... c ^ 2 \ over 2 } \ approx {m \ over M} E \ approx {1 \ over {1000}} eV \ egroup \ end {displaymath}$

Teoria da Perturbação Dependente do Tempo

Usamos a teoria de perturbação independente do tempo para encontrar as mudanças de energia dos estados e para encontrar a mudança nos autoestados de energia na presença de uma pequena perturbação. Agora consideramos o caso de uma perturbação ${\ cal V}$ que é dependente do tempo . Tal perturbação pode causar transições entre estados próprios de energia . Vamos calcular a taxa dessas transições.

Derivamos uma equação para a taxa de variação da amplitude no estado próprio de $n ^ {th}$ energia .

$\ begin {eqnarray *} i \ hbar {\ partial c_n (t) \ over \ partial t} & = & \ sum \ limits_k {\ cal V} _ {nk} (t) c_k (t) e ^ {i ( E_n-E_k) t / \ hbar} \\ \ end {eqnarray *}$

Supondo que no $\ bgroup \ color {black} $ t = 0 $ \ egroup$ sistema quântico comece em algum estado inicial $\ phi_i$ , derivamos a amplitude para estar em um estado final $\ phi_n$ .

$\ begin {displaymath} \ bgroup \ color {black} c_n (t) = {1 \ over i \ hbar} \ int \ limits_0 ^ te ^ {i \ omega_ {ni} t '} {\ cal V} _ {ni } (t ') dt' \ egroup \ end {displaymath}$

Um caso importante de potencial dependente do tempo é uma perturbação oscilante senoidal pura (harmônica) . Podemos compensar qualquer dependência de tempo a partir de uma combinação linear de ondas seno e cosseno. Com alguns cálculos, derivamos a taxa de transição em um potencial harmônico de frequência $\ bgroup \ color {black} $ \ omega $ \ egroup$ .

$\ begin {displaymath} \ bgroup \ color {black} \ Gamma_ {i \ rightarrow n} \ equiv {dP_n \ over ... ... = {2 \ pi V_ {ni} ^ 2 \ over \ hbar} \ delta (E_n-E_i + \ hbar \ omega) \ egroup \ end {displaymath}$

Ele contém uma função delta de conservação de energia. A função delta pode parecer estranha. A taxa de transição seria zero se a energia não fosse conservada e infinita se a energia fosse exatamente conservada. Podemos entender isso se houver uma função de distribuição do $\ bgroup \ color {black} $ P (\ omega) $ \ egroup$ potencial de perturbação ou se houver um contínuo de estados finais que precisamos integrar. Em qualquer caso, a função delta nos ajuda a fazer a integral de maneira simples.

Radiação em átomos

A interação de átomos com ondas eletromagnéticas pode ser calculada usando a teoria de perturbação dependente do tempo. O problema atômico é resolvido na ausência de ondas EM, então os termos do vetor potencial no hamiltoniano podem ser tratados como uma perturbação.

$\ begin {displaymath} \ bgroup \ color {black} H = {1 \ over 2m} \ left (\ vec {p} + {e \ over c} \ vec {A} \ right) ^ 2 + V (r) . \ egroup \ end {displaymath}$

Em um medidor em que , a perturbação é $\ bgroup \ color {black} $ \ vec {\ nabla} \ cdot \ vec {A} = 0 $ \ egroup$

$\ begin {displaymath} \ bgroup \ color {black} {\ cal V} = {e \ over mc} \ vec {A} \ cdot \ vec {p} + {e ^ 2 \ over 2mc ^ 2} A ^ 2 . \ egroup \ end {displaymath}$

Para a maioria dos decaimentos atômicos, o $\ bgroup \ color {black} $ A ^ 2 $ \ egroup$ termo pode ser desprezado, pois é muito menor do que o termo. Tanto a decadência dos estados atômicos excitados com a emissão da radiação quanto a excitação dos átomos com a absorção da radiação podem ser calculadas. $\ bgroup \ color {black} $ \ vec {A} \ cdot \ vec {p} $ \ egroup$

Um campo EM arbitrário pode ser analisado por Fourier para fornecer uma soma de componentes de frequência definida. Considere o potencial vetorial para um desses componentes ,. A energia no campo é . Se o campo for quantizado (como mostraremos mais tarde) com fótons de energia , podemos escrever a intensidade do campo em termos do número de fótons .

$\ bgroup \ color {black} $ \ vec {A} (\ vec {r}, t) \ equiv 2 \ vec {A} _0 \ cos (\ vec {k} \ cdot \ vec {r} - \ omega t ) $ \ egroup$

$\ bgroup \ color {black} $ Energia = {\ omega ^ 2 \ over 2 \ pi c ^ 2} V \ vert A_0 \ vert ^ 2 $ \ egroup$ $\ bgroup \ color {black} $ E = \ hbar \ omega $ \ egroup$ $\ bgroup \ color {black} $ N $ \ egroup$

$\ begin {eqnarray *} \ vec {A} (\ vec {r}, t) & = & \ left [{2 \ pi \ hbar c ^ 2 N \ over \ omega V} \ righ ... ... ot \ vec {r} - \ omega t)} + e ^ {- i (\ vec {k} \ cdot \ vec {r} - \ omega t)} \ right) \ end {eqnarray *}$

A direção do campo é dada pelo vetor de polarização unitária . O termo cosseno foi dividido em exponenciais positivas e negativas. Na teoria das perturbações dependentes do tempo, o exponencial positivo corresponde à absorção de um fóton e excitação do átomo e o exponencial negativo corresponde à emissão de um fóton e decaimento do átomo para um estado de energia inferior. $\ bgroup \ color {black} $ \ hat {\ epsilon} $ \ egroup$

Pense no campo EM como um oscilador harmônico em cada frequência, o exponencial negativo corresponde a um operador de aumento para o campo e o exponencial positivo a um operador de redução. Em analogia ao oscilador harmônico 1D quântico, substituímos $\ bgroup \ color {black} $ \ sqrt {N} $ \ egroup$ por $\ bgroup \ color {black} $ \ sqrt {N + 1} $ \ egroup$ no caso do operador ascendente.

$\ begin {eqnarray *} \ vec {A} (\ vec {r}, t) & = & \ left [{2 \ pi \ hbar c ^ 2 \ over \ omega V} \ right] ... ... \ omega t)} + \ sqrt {N + 1} e ^ {- i (\ vec {k} \ cdot \ vec {r} - \ omega t)} \ right) \ end {eqnarray *}$

Com esta mudança, que posteriormente será justificada com a quantização do campo, ocorre uma perturbação mesmo sem campo aplicado ( $\ bgroup \ color {black} $ N = 0 $ \ egroup$ )

$\ begin {displaymath} \ bgroup \ color {black} {\ cal V} _ {N = 0} = V_ {N = 0} e ^ {i \ omega t} = {e \ o ... ... vec {k} \ cdot \ vec {r} - \ omega t)} \ hat {\ epsilon} \ cdot \ vec {p} \ egroup \ end {displaymath}$

que pode causar decaimento de estados atômicos.

Conectando este $\ bgroup \ color {black} $ N = 0 $ \ egroup$ campo nas equações de perturbação dependentes do tempo de primeira ordem, a taxa de decaimento para um estado atômico pode ser calculada.

$\ begin {eqnarray *} \ Gamma_ {i \ rightarrow n} & = & {(2 \ pi) ^ 2e ^ 2 \ over m ^ 2 \ omega V} \ vert ... ... vec {p} \ vert \ phi_i \ rangle \ vert ^ 2 \: \ delta (E_n-E_i + \ hbar \ omega) \\ \ end {eqnarray *}$

O quadrado absoluto da integral de tempo da teoria de perturbação produz a função delta de conservação de energia.

Para obter a taxa de decaimento total, devemos somar os estados finais permitidos. Podemos supor que o átomo permanece em repouso como uma boa aproximação, mas os estados finais do fóton devem ser cuidadosamente considerados. Aplicando condições de contorno periódicas em um volume cúbico $\ bgroup \ color {black} $ V $ \ egroup$ , a integral sobre os estados finais pode ser feita conforme indicado abaixo.

$\ begin {eqnarray *} k_xL = 2 \ pi n_x & dn_x = {L \ over 2 \ pi} dk_x \\ k_yL = 2 \ pi n_y & dn _... ... i) ^ 3} d ^ 3k \\ & \ Gamma_ {tot} = \ int \ Gamma_ {i \ rightarrow n} d ^ 3n \\ \ end {eqnarray *}$

Com esta integral de espaço de fase feita com a ajuda da função delta, a fórmula geral para a taxa de decaimento é

$\ begin {displaymath} \ bgroup \ color {black} \ Gamma_ {tot} = {e ^ 2 (E_i-E_n) \ over 2 \ pi \ hba ... ... on} ^ {(\ lambda)} \ cdot \ vec {p_e} \ vert \ phi_i \ rangle \ vert ^ 2. \ egroup \ end {displaymath}$

Essa taxa de decaimento ainda contém a integral sobre as direções dos fótons e uma soma sobre a polarização do estado final.

O cálculo do elemento da matriz atômica é geralmente feito na aproximação do dipolo elétrico

$\ begin {displaymath} \ bgroup \ color {black} e ^ {- i \ vec {k} \ cdot \ vec {r}} \ approx 1 \ egroup \ end {displaymath}$

o que é válido se o comprimento de onda do fóton for muito maior que o tamanho do átomo. Com a ajuda de algumas relações de comutação, a fórmula da taxa de decaimento torna-se

$\ begin {displaymath} \ bgroup \ color {black} \ Gamma_ {tot} = {\ alpha \ omega_ {in} ^ 3 \ over 2 ... ... cdot \ langle \ phi_n \ vert \ vec {r} \ vert \ phi_i \ rangle \ vert ^ 2. \ egroup \ end {displaymath}$

O elemento da matriz atômica do operador vetorial $\ bgroup \ color {black} $ \ vec {r} $ \ egroup$ é zero, a menos que certas restrições no momento angular dos estados inicial e final sejam satisfeitas. As regras de seleção para transições de dipolo elétrico (E1) são:

$\ begin {displaymath} \ bgroup \ color {black} \ Delta \ ell = \ pm 1 \ qquad \ Delta m = 0, \ pm 1 \ qquad \ Delta s = 0. \ egroup \ end {displaymath}$

Este é o resultado do teorema de Wigner-Eckart que afirma que o elemento da matriz de um operador vetorial $\ bgroup \ color {black} $ V ^ q $ \ egroup$ , onde o inteiro $\ bgroup \ color {black} $ q $ \ egroup$ vai de -1 a +1, é dado por

$\ begin {displaymath} \ bgroup \ color {black} \ langle \ alpha'j'm '\ vert V ^ q \ vert \ alpha ... ... \ langle \ alpha'j' \ vert \ vert V \ vert \ vert \ alpha j \ rangle \ egroup \ end {displaymath}$

Aqui $\ bgroup \ color {black} $ \ alpha $ \ egroup$ representa todos os (outros) números quânticos do estado, não os números quânticos do momento angular. No caso de um operador espacial simples como $\ bgroup \ color {black} $ \ vec {r} $ \ egroup$ , apenas o momento angular orbital está envolvido.

$\ begin {displaymath} \ bgroup \ color {black} \ Gamma_ {2p \ rightarrow 1s} = {2 \ alpha \ omega ... ... vert 4 \ sqrt {6} \ left ({2 \ over 3} \ direita) ^ 5a_0 \ right \ vert ^ 2 \ egroup \ end {displaymath}$

Derivamos um resultado simples para a taxa de decaimento total de um estado, somada sobre a polarização final do fóton e integrada sobre a direção do fóton

$\ begin {displaymath} \ bgroup \ color {black} \ Gamma_ {tot} = {4 \ alpha \ omega_ {in} ^ 3 \ over 3c ^ 2} \ vert \ vec {r} _ {ni} \ vert ^ 2 \ egroup \ end {displaymath}$

Isso pode ser usado para calcular facilmente as taxas de decaimento do hidrogênio, por exemplo, a taxa de decaimento 2p.

$\ begin {displaymath} \ bgroup \ color {black} \ Gamma_ {2p \ rightarrow 1s} = {4 \ alpha \ omega ... ... vert 4 \ sqrt {6} \ left ({2 \ over 3} \ direita) ^ 5a_0 \ right \ vert ^ 2 \ egroup \ end {displaymath}$

A taxa de decaimento total está relacionada à largura de energia de um estado excitado, como pode ser esperado do princípio da incerteza. A largura total na metade do máximo (FWHM) da distribuição de energia de um estado é . A distribuição de frequência segue uma distribuição Breit-Wigner. $\ bgroup \ color {black} $ \ hbar \ Gamma_ {tot} $ \ egroup$

$\ begin {displaymath} \ bgroup \ color {black} I_i (\ omega) = \ vert \ phi_i (\ omega) \ vert ^ 2 = {1 \ over (\ omega- \ omega_0) ^ 2 + {\ Gamma ^ 2 \ over 4}} \ egroup \ end {displaymath}$

Além da largura de energia inerente de um estado, outros efeitos podem influenciar as larguras medidas, incluindo alargamento de colisão, alargamento Doppler e recuo atômico.

A teoria quântica da radiação EM pode ser usada para entender muitos fenômenos, incluindo distribuições angulares de fótons, polarização de fótons, LASERs, o efeito Mössbauer, o efeito fotoelétrico, o espalhamento de luz e a absorção de raios-x.

Teoria de Campo Clássica

Uma revisão da teoria de campo clássica é útil para fundamentar nosso desenvolvimento de teorias quânticas relativísticas de campo para fótons e elétrons. Vamos trabalhar com 4 vetores como o vetor de coordenadas abaixo

$\ begin {displaymath} \ bgroup \ color {black} (x_1, x_2, x_3, x_4) = (x, y, z, ict) \ egroup \ end {displaymath}$

usando o $\ bgroup \ color {black} $ i $ \ egroup$ para obter a $\ bgroup \ color {black} $ - $ \ egroup$ no termo de tempo em um produto escalar (em vez de usar um tensor métrico).

Uma densidade Lagrangiana escalar de Lorentz será derivada para cada teoria de campo que construirmos. Do Lagrangiano podemos derivar uma equação de campo chamada equação de Euler-Lagrange.

$\ begin {displaymath} \ bgroup \ color {black} {\ partial \ over \ partial x_ \ mu} \ left ({\ par ... ... l x_ \ mu)} \ right) - {\ partial {\ cal L} \ over \ partial \ phi} = 0 \ egroup \ end {displaymath}$

O Lagrangiano para um campo escalar massivo $\ bgroup \ color {black} $ \ phi $ \ egroup$ pode ser deduzido do requisito de que seja um campo escalar

$\ begin {displaymath} \ bgroup \ color {black} {\ cal L} = - {1 \ over 2} \ left ({\ partial \ phi \ ... ... al \ phi \ over \ partial x_ \ nu } + \ mu ^ 2 \ phi ^ 2 \ right) + \ phi \ rho \ egroup \ end {displaymath}$

onde o último termo é a interação com uma fonte. A equação de Euler-Lagrange dá

$\ begin {displaymath} \ bgroup \ color {black} {\ parcial \ over \ parcial x_ \ mu} {\ parcial \ over \ parcial x_ \ mu} \ phi- \ mu ^ 2 \ phi = \ rho \ egroup \ end {displaymath}$

que é conhecida como a equação de Klein-Gordon com uma fonte e é uma equação relativística razoável para um campo escalar.

Usando as transformadas de Fourier, o campo de uma fonte pontual pode ser calculado

$\ begin {eqnarray *} \ phi (\ vec {x}) & = & {- G e ^ {- \ mu r} \ over 4 \ pi r} \\ \ end {eqnarray *}$

Este é um campo que cai muito mais rápido do que $\ bgroup \ color {black} $ {1 \ over r} $ \ egroup$ . Um campo escalar massivo cai exponencialmente e quanto maior a massa, mais rápida é a queda. Isso se encaixa muito bem na forma da força entre os núcleons, embora a força nuclear real precise de um estudo muito mais detalhado.

O Campo Eletromagnético Clássico

Para o estudo do campo de Maxwell, é mais conveniente fazer uma pequena modificação no sistema de unidades usado. Nas unidades Racionalizadas de Heaviside-Lorentz, os campos são todos reduzidos por um fator de $\ bgroup \ color {black} $ \ sqrt {4 \ pi} $ \ egroup$ e as cargas são aumentadas pelo mesmo fator. Com essa mudança, as equações de Maxwell, assim como as Lagrangianas que usamos, são simplificadas. Teria simplificado muitas coisas se Maxwell tivesse começado com este conjunto de unidades.

Como é bem conhecido pela eletricidade e magnetismo clássicos, os componentes do campo elétrico e magnético são, na verdade, elementos de um tensor de Lorentz de categoria 2.

$\ begin {displaymath} \ bgroup \ color {black} F _ {\ mu \ nu} = \ pmatrix {0 & B_z & -B_y & -i ... ... cr B_y & -B_x & 0 & -iE_z \ cr iE_x & iE_y & iE_z & 0} \ egroup \ end {displaymath}$

Este tensor de campo pode simplesmente ser escrito em termos do potencial vetorial , (que é um vetor de Lorentz).

$\ begin {eqnarray *} A_ \ mu & = & (\ vec {A}, i \ phi) \\ F _ {\ mu \ nu} & = & = & {\ parcial A_ \ nu \ over \ parcial x_ \ mu} - { \ parcial A_ \ mu \ over \ parcial x_ \ nu} \\ \ end {eqnarray *}$

Observe que $\ bgroup \ color {black} $ F _ {\ mu \ nu} $ \ egroup$ é automaticamente anti-simétrico sob o intercâmbio dos índices.

Com os campos derivados do potencial vetorial, duas das equações de Maxwell são satisfeitas automaticamente. As duas equações restantes podem ser escritas como uma equação de 4 vetores.

$\ begin {displaymath} \ bgroup \ color {black} {\ partial F _ {\ mu \ nu} \ over \ partial x_ \ nu} = {j_ \ mu \ over c} \ egroup \ end {displaymath}$

Agora queremos escolher um Lagrangiano escalar. Visto que E&M é uma teoria bem compreendida, o Lagrangiano que é conhecido por fornecer as equações corretas também é conhecido.

$\ begin {displaymath} \ bgroup \ color {black} {\ cal L} = - {1 \ over 4} F _ {\ mu \ nu} F _ {\ mu \ nu} + {1 \ over c} j_ \ mu A_ \ mu \ egroup \ end {displaymath}$

Observe que (além da velocidade da luz não ser definida como 1), o Lagrangiano não contém constantes desnecessárias neste conjunto de unidades. O último termo é um termo fonte que fornece a interação entre o campo EM e partículas carregadas. Ao trabalhar com este Lagrangiano, trataremos cada componente de $\ bgroup \ color {black} $ A $ \ egroup$ como um campo independente. Neste caso, a equação de Euler-Lagrange é a equação de Maxwell conforme escrita acima.

A densidade hamiltoniana de campo livre pode ser calculada de acordo com a prescrição padrão que produz

$\ begin {eqnarray *} {\ partial {\ cal L} \ over \ partial (\ partial A_ \ mu / \ partial x_4)} & = ... ... A_ \ mu \ over \ partial x_4} - {\ cal L} \\ & = & {1 \ over 2} (E ^ 2 + B ^ 2) \\ \ end {eqnarray *}$

se não houver termos de origem na região.

A simetria de calibre pode ser usada para colocar uma condição no potencial do vetor.

$\ begin {displaymath} \ bgroup \ color {black} {\ partial A_ \ nu \ over \ partial x_ \ nu} = 0. \ egroup \ end {displaymath}$

Isso é chamado de condição de Lorentz . Mesmo com isso satisfeito, ainda há uma liberdade de medida substancial possível. As transformações de calibre podem ser feitas conforme mostrado abaixo.

$\ begin {eqnarray *} A_ \ mu \ rightarrow A_ \ mu + {\ partial \ Lambda \ over \ partial x_ \ mu} \\ \ Box \ Lambda = 0 \\ \ end {eqnarray *}$

Quantização do campo EM

O hamiltoniano para o campo de Maxwell pode ser usado para quantizar o campo da mesma forma que a mecânica de uma onda dimensional foi quantizada. O campo de radiação pode ser mostrado como a parte transversal do campo, $\ bgroup \ color {black} $ \ vec {A} _ \ perp $ \ egroup$ enquanto as cargas estáticas dão origem a $\ bgroup \ color {black} $ A_ \ parallel $ \ egroup$ e $\ bgroup \ color {black} $ A_0 $ \ egroup$ .

Nós decompor o campo de radiação em seus componentes de Fourier

$\ begin {displaymath} \ bgroup \ color {black} \ vec {A} (\ vec {x}, t) = {1 \ over \ sqrt {V}} \ sum \ ... ... c {x} } + c_ {k, \ alpha} ^ * (t) e ^ {- i \ vec {k} \ cdot \ vec {x}} \ right) \ egroup \ end {displaymath}$

onde são vetores unitários reais, e é o coeficiente da onda com o vetor de onda e o vetor de polarização . Uma vez que o vector de onda é escolhido, os dois vectores de polarização deve ser escolhido de modo que , e formar um sistema ortogonal destro . $\ bgroup \ color {black} $ \ hat {\ epsilon} ^ {(\ alpha)} $ \ egroup$ $\ bgroup \ color {black} $ c_ {k, \ alpha} $ \ egroup$ $\ bgroup \ color {black} $ \ vec {k} $ \ egroup$ $\ bgroup \ color {black} $ \ hat {\ epsilon} ^ {(\ alpha)} $ \ egroup$ $\ bgroup \ color {black} $ \ hat {\ epsilon} ^ {(1)} $ \ egroup$ $\ bgroup \ color {black} $ \ hat {\ epsilon} ^ {(2)} $ \ egroup$ $\ bgroup \ color {black} $ \ vec {k} $ \ egroup$

Conectando a decomposição de Fourier na fórmula para a densidade hamiltoniana e usando a natureza transversal do campo de radiação, podemos calcular o hamiltoniano (densidade integrada sobre o volume).

$\ begin {eqnarray *} H & = & \ sum \ limits_ {k, \ alpha} \ left ({\ omega \ over c} \ right) ^ 2 \ left ... ... (t) c_ {k, \ alpha} ^ * (t) + c_ {k, \ alpha} ^ * (t) c_ {k, \ alpha} (t) \ right] \\ \ end {eqnarray *}$

Este hamiltoniano será usado para quantizar o campo EM. Ao calcular o hamiltoniano, foi tomado cuidado para não comutar os coeficientes de Fourier e seus conjugados.

A coordenada canônica e os momentos podem ser encontrados

$\ begin {displaymath} \ bgroup \ color {black} Q_ {k, \ alpha} = {1 \ over c} (c_ {k, \ alpha} + c_ {k, \ alpha} ^ *) \ egroup \ end { displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} P_ {k, \ alpha} = - {i \ omega \ over c} (c_ {k, \ alpha} -c_ {k, \ alpha} ^ *) \ egroup \ end {displaymath}$

para o oscilador harmônico em cada frequência. Assumimos que uma coordenada e seu momento conjugado têm o mesmo comutador que na mecânica das ondas e que as coordenadas de diferentes osciladores comutam.

$\ begin {eqnarray *} \ left [Q_ {k, \ alpha}, P_ {k ', \ alpha'} \ right] & = & i \ hbar \ delta_ {kk '} \ ... ... \ right] & = & 0 \\ \ left [P_ {k, \ alpha}, P_ {k ', \ alpha'} \ right] & = & 0 \\ \ end {eqnarray *}$

Como foi feito para o oscilador harmônico 1D, escrevemos o hamiltoniano em termos de operadores de aumento e redução que têm as mesmas relações de comutação do oscilador harmônico 1D.

$\ begin {eqnarray *} a_ {k, \ alpha} & = & {1 \ over \ sqrt {2 \ hbar \ omega}} (\ omega Q_ {k, \ alpha} ... ... \ alpha '} ^ \ dagger \ right] & = & \ delta_ {kk '} \ delta _ {\ alpha \ alpha'} \\ \ end {eqnarray *}$

Isso significa que tudo o que sabemos sobre os operadores de subida e descida se aplica aqui. As energias estão em etapas $\ bgroup \ color {black} $ \ hbar \ omega $ \ egroup$ e deve haver um estado fundamental. Os estados podem ser rotulados por um número quântico $\ bgroup \ color {black} $ n_ {k, \ alpha} $ \ egroup$ .

$\ begin {eqnarray *} H & = & \ left (a_ {k, \ alpha} ^ \ dagger a_ {k, \ alpha} + {1 \ over 2} \ right) \ ... ... r \ omega \ \ N_ {k, \ alpha} & = & a_ {k, \ alpha} ^ \ dagger a_ {k, \ alpha} \\ \ end {eqnarray *}$

Os coeficientes de Fourier agora podem ser escritos em termos dos operadores de aumento e redução para o campo.

$\ begin {eqnarray *} c_ {k, \ alpha} & = & \ sqrt {\ hbar c ^ 2 \ over 2 \ omega} a_ {k, \ alpha} \\ c ... ... {x}} + a_ {k, \ alpha} ^ \ dagger (t) e ^ {- i \ vec {k} \ cdot \ vec {x}} \ right) \\ \ end {eqnarray *}$

$\ begin {eqnarray *} H & = & {1 \ over 2} \ sum \ limits_ {k, \ alpha} \ hbar \ omega \ left [a_ {k, \ alp ... ... ts_ {k, \ alpha } \ hbar \ omega \ left (N_ {k, \ alpha} + {1 \ over 2} \ right) \\ \ end {eqnarray *}$

Os estados do campo são dados pelo número de ocupação de cada estado possível de fóton.

$\ begin {displaymath} \ bgroup \ color {black} \ vert n_ {k_1, \ alpha_1}, n_ {k_2, \ alpha_2}, ... ..., \ alpha_i}} \ over \ sqrt {n_ {k_i, \ alpha_i}!}} \ vert \ rangle \ egroup \ end {displaymath}$

Qualquer estado pode ser construído operando com operadores de criação no estado de vácuo. Qualquer estado com vários fótons será automaticamente simétrico sob o intercâmbio de pares de fótons porque os operadores comutam.

$\ begin {displaymath} \ bgroup \ color {black} a_ {k, \ alpha} ^ \ dagger a_ {k ', \ alpha'} ^ \ dag ... ... {k ', \ alpha'} ^ \ punhal a_ {k, \ alpha} ^ \ dagger \ vert \ rangle \ egroup \ end {displaymath}$

Este é essencialmente o mesmo resultado que nossa estimativa anterior de colocar um $\ bgroup \ color {black} $ n + 1 $ \ egroup$ no operador de emissão.

Agora podemos escrever o campo de radiação quantizado em termos dos operadores em $\ bgroup \ color {black} $ t = 0 $ \ egroup$ .

$\ begin {eqnarray *} A_ \ mu & = & {1 \ over \ sqrt {V}} \ sum \ limits_ {k \ alpha} \ sqrt {\ hbar c ^ 2 \ ... ... ho x_ \ rho } + a_ {k, \ alpha} ^ \ dagger (0) e ^ {- ik_ \ rho x_ \ rho} \ right) \\ \ end {eqnarray *}$

Além da aproximação do dipolo elétrico, o próximo termo na expansão é . Este termo é dividido de acordo com suas propriedades de rotação e transformação de Lorentz em termos de Quadrupolo Elétrico e Dipolo Magnético. A interação do spin do elétron com o campo magnético é da mesma ordem e deve ser incluída junto com os termos E2 e M1. $\ bgroup \ color {black} $ e ^ {i \ vec {k} \ cdot \ vec {x}} $ \ egroup$ $\ bgroup \ color {black} $ i \ vec {k} \ cdot \ vec {x} $ \ egroup$

$\ begin {displaymath} \ bgroup \ color {black} {e \ hbar \ over 2mc} (\ vec {k} \ times \ hat {\ epsilon} ^ {(\ lambda)}) \ cdot \ vec {\ sigma} \ egroup \ end {displaymath}$

O termo Quadrupolo Elétrico (E2) não muda a paridade e nos dá a regra de seleção.

$\ begin {displaymath} \ bgroup \ color {preto} \ vert \ ell_n- \ ell_i \ vert \ leq 2 \ leq \ ell_n + \ ell_i \ egroup \ end {displaymath}$

O termo dipolo magnético (M1) não altera a paridade, mas pode alterar o spin. Como é um operador vetorial (axial), ele muda o momento angular em 0, +1 ou -1 unidade.

O campo quantizado é muito útil na derivação da fórmula de radiação do corpo negro de Plank que deu início à revolução quântica. Ao equilibrar as taxas de reação proporcionais $\ bgroup \ color {black} $ N $ \ egroup$ e $\ bgroup \ color {black} $ N + 1 $ \ egroup$ para absorção e emissão em equilíbrio, a densidade de energia no campo de radiação dentro de uma cavidade é facilmente derivada.

$\ begin {eqnarray *} U (\ nu) = U (\ omega) {d \ omega \ over d \ nu} = {8 \ pi \ over c ^ 3} {h \ nu ^ 3 \ over e ^ {\ hbar \ omega / kT} -1} \\ \ end {eqnarray *}$

Dispersão de fótons

34 O campo de fótons quantizado pode ser usado para calcular a seção transversal para espalhamento de fótons. A aproximação de dipolo elétrico é usada para simplificar o elemento da matriz atômica em baixa energia, onde o comprimento de onda é longo em comparação com o tamanho atômico.

Para espalhar um fóton, o campo deve agir duas vezes, uma para aniquilar o fóton de estado inicial e outra para criar o fóton de estado final. Uma vez que o campo quantizado contém operadores de criação e aniquilação,

$\ begin {displaymath} \ bgroup \ color {black} A_ \ mu (x) = {1 \ over \ sqrt {V}} \ sum \ limits_ {k \ ... ... \ rho} + a_ {k, \ alpha} ^ \ dagger (0) e ^ {- ik_ \ rho x_ \ rho} \ right) \ egroup \ end {displaymath}$

tanto o $\ bgroup \ color {black} $ A ^ 2 $ \ egroup$ termo de primeira ordem quanto o termo de segunda ordem podem contribuir para a dispersão. Ambas as amplitudes são adequadas . $\ bgroup \ color {black} $ \ vec {A} \ cdot \ vec {p} $ \ egroup$ $\ bgroup \ color {black} $ e ^ 2 $ \ egroup$

O elemento de matriz do $\ bgroup \ color {black} $ A ^ 2 $ \ egroup$ termo para ir de um fóton de vetor de onda $\ bgroup \ color {black} $ \ vec {k} $ \ egroup$ e um estado atômico $\ bgroup \ color {black} $ i $ \ egroup$ a um fóton espalhado de vetor de onda $\ bgroup \ color {black} $ \ vec {k '} $ \ egroup$ e um estado atômico $\ bgroup \ color {black} $ n $ \ egroup$ é particularmente simples, pois não contém coordenadas atômicas ou momentos.

$\ begin {eqnarray *} {e ^ 2 \ over 2mc ^ 2} \ langle n; \ vec {k} '\ hat {\ epsilon} ^ {(\ alpha')} \ v ... ... \ epsilon_ \ mu ^ {(\ alpha ')} e ^ {- i (\ omega- \ omega') t} \ delta_ {ni} \\ \ end {eqnarray *}$

Os termos de segunda ordem podem mudar os estados atômicos por causa do $\ bgroup \ color {black} $ \ vec {p} $ \ egroup$ operador.

A seção transversal para espalhamento de fótons é então dada pelo

$\ begin {displaymath} \ bgroup \ color {black} {d \ sigma \ over d \ Omega} = \ left ({e ^ 2 \ over ... ... ângulo \ over \ omega_ {ji} + \ omega ' } \ right] \ right \ vert ^ 2 \ egroup \ end {displaymath}$

Fórmula de Kramers-Heisenberg . Os três termos vêm dos três diagramas de Feynman que contribuem para o espalhamento na ordem $\ bgroup \ color {black} $ e ^ 2 $ \ egroup$ .

Este resultado pode ser especializado para o caso de espalhamento elástico, com o auxílio de alguns comutadores.

$\ begin {displaymath} \ bgroup \ color {black} {d \ sigma_ {elas} \ over d \ Omega} = \ left ({e ^ ... ... rangle \ over \ omega_ {ji} + \ omega} \ right] \ right \ vert ^ 2 \ egroup \ end {displaymath}$

Lord Rayleigh calculou o espalhamento elástico de baixa energia da luz dos átomos usando o eletromagnetismo clássico. Se a energia do fóton espalhado é menor do que a energia necessária para excitar o átomo, então a seção transversal é proporcional a $\ bgroup \ color {black} $ \ omega ^ 4 $ \ egroup$ , de modo que a luz azul se espalha mais do que a luz vermelha nos gases incolores em nossa atmosfera.

Se a energia do fóton espalhado for muito maior do que a energia de ligação do átomo, $\ bgroup \ color {black} $ \ omega »1 $ \ egroup$ eV. então a seção transversal se aproxima daquela do espalhamento de um elétron livre, Thomson Scattering .

$\ begin {displaymath} \ bgroup \ color {black} {d \ sigma \ over d \ Omega} = \ left ({e ^ 2 \ over ... ... ft \ vert \ hat {\ epsilon} \ cdot \ hat {\ epsilon} '\ right \ vert ^ 2 \ egroup \ end {displaymath}$

O espalhamento é aproximadamente independente de energia e a única dependência angular é da polarização. A luz dispersa pode ser polarizada mesmo que a luz incidente não o seja.

Electron Self Energy

Mesmo no eletromagnetismo clássico, se alguém puder calcular a energia necessária para montar um elétron, o resultado é infinito, mas os elétrons existem. A correção da autoenergia quântica também é infinita, embora possa ser considerada finita se aceitarmos o fato de que nossas teorias não são válidas até energias infinitas.

A correção da autoenergia quântica tem efeitos importantes e mensuráveis. Ele causa mudanças observáveis de energia no Hidrogênio e nos ajuda a resolver o problema dos infinitos devido aos denominadores de energia de estados intermediários.

As equações diferenciais acopladas da teoria de perturbação de primeira ordem para o estado em estudo $\ bgroup \ color {black} $ \ phi_n $ \ egroup$ e estados intermediários $\ bgroup \ color {black} $ \ psi_j $ \ egroup$ podem ser resolvidas para a correção de energia própria.

$\ begin {eqnarray *} \ Delta E_n & = & \ sum \ limits _ {\ vec {k}, \ alpha} \ sum \ limits_j \ vert H _... ... {1-e ^ {i (\ omega_ {nj } - \ omega) t} \ over \ hbar (\ omega_ {nj} - \ omega)} \\ \ end {eqnarray *}$

O resultado é, em geral, complexo. A parte imaginária da correção da energia própria está diretamente relacionada à largura do estado.

$\ begin {displaymath} \ bgroup \ color {black} - {2 \ over \ hbar} \ Im (\ Delta E_n) = \ Gamma_n \ egroup \ end {displaymath}$

A dependência do tempo da função de onda para o estado $\ psi_n$ é modificada pela correção de energia própria .

$\ begin {displaymath} \ bgroup \ color {black} \ psi_n (\ vec {x}, t) = \ psi_n (\ vec {x}) e ^ {- i (E_n + \ Re (\ Delta E_n)) t / \ hbar} e ^ {- \ Gamma_nt \ over 2} \ egroup \ end {displaymath}$

Isso nos dá o comportamento de decaimento exponencial que esperamos, evitando que as seções transversais de espalhamento ressonante cheguem ao infinito .

A parte real da correção deve ser estudada para compreender as mudanças relativas de energia dos estados. É a diferença entre a energia própria do elétron ligado e a de um elétron livre na qual estamos interessados. A correção de autoenergia para uma partícula livre pode ser calculada.

$\ begin {eqnarray *} \ Delta E_ {free} & = & - {2 \ alpha E_ {cut-off} \ over 3 \ pi m ^ 2c ^ 2} p ^ 2 \\ \ end {eqnarray *}$

Nós automaticamente contabilizamos essa correção por uma mudança na massa observada do elétron. Para a definição não relativística da energia de um elétron livre, um aumento na massa diminui a energia.

$\ begin {eqnarray *} m_ {obs} & = & (1+ {4 \ alpha E_ {cut-off} \ over 3 \ pi mc ^ 2}) m_ {bare} \\ \ end {eqnarray *}$

Se cortarmos a integral em , a correção da massa é de apenas cerca de 0,3% ,

Uma vez que a massa observada do elétron já é responsável pela maior parte da correção de energia própria para um estado ligado, devemos corrigir esse efeito para evitar a contagem dupla da correção. A correção de energia própria para um estado limitado então é.

$\ begin {eqnarray *} \ Delta E_n ^ {(obs)} & = & \ Delta E_n + {2 \ alpha E_ {corte} \ over 3 \ pi m ^ 2c ^ 2} \ langle n \ vert p ^ 2 \ vert n \ rangle \\ \ end {eqnarray *}$

Em 1947, Willis E. Lamb e RC Retherford usaram técnicas de microondas para determinar a divisão entre os estados $2S_ {1 \ sobre 2}$ e $2P_ {1 \ sobre 2}$ no hidrogênio . O resultado pode ser bem explicado pela correção da energia própria, pelo menos quando a mecânica quântica relativística é usada. Nosso cálculo não relativístico fornece uma explicação qualitativa do efeito.

$\ begin {displaymath} \ bgroup \ color {black} \ Delta E_n ^ {(obs)} = {4 \ alpha ^ 5 \ over 3 \ pi ... ... {\ omega} _ {nj}} \ right ) + {11 \ over 24} - {1 \ over 5} \ right) mc ^ 2 \ egroup \ end {displaymath}$

A Equação de Dirac

Nosso objetivo é encontrar o análogo da equação de Schrödinger para partículas relativísticas de spin meio, entretanto, devemos notar que mesmo na equação de Schrödinger, a interação do campo com o spin era bastante ad hoc. Não houve explicação da razão giromagnética de 2. Pode-se incorporar o spin na equação não relativística usando o hamiltoniano de Schrödinger-Pauli que contém o produto escalar das matrizes de Pauli com o operador momentum.

$\ begin {displaymath} \ bgroup \ color {black} H = {1 \ over 2m} \ left (\ vec {\ sigma} \ cdot [\ ve ... ... {e \ over c} \ vec {A } (\ vec {r}, t)] \ right) ^ 2-e \ phi (\ vec {r}, t) \ egroup \ end {displaymath}$

Um pequeno cálculo mostra que isso dá a interação correta com o spin.

$\ begin {displaymath} \ bgroup \ color {black} H = {1 \ over 2m} [\ vec {p} + {e \ over c} \ vec {A} (... ... t) + {e \ hbar \ over 2mc} \ vec {\ sigma} \ cdot \ vec {B} (\ vec {r}, t) \ egroup \ end {displaymath}$

Este hamiltoniano atua em um spinor de dois componentes.

Podemos estender este conceito para usar a equação da energia relativística . A ideia é substituir $\ bgroup \ color {black} $ \ vec {p} $ \ egroup$ com na equação da energia relativística. $\ bgroup \ color {black} $ \ vec {\ sigma} \ cdot \ vec {p} $ \ egroup$

$\ begin {eqnarray *} \ left ({E \ over c} \ right) ^ 2-p ^ 2 = (mc) ^ 2 \\ \ left ({E \ over c} - \ vec ... ... 0} -i \ hbar \ vec {\ sigma} \ cdot \ vec {\ nabla} \ right) \ phi = (mc) ^ 2 \ phi \\ \ end {eqnarray *}$

Em vez de uma equação de segunda ordem na derivada do tempo, podemos fazer uma equação de primeira ordem, como a equação de Schrödinger, estendendo essa equação para quatro componentes.

$\ begin {eqnarray *} \ phi ^ {(L)} & = & \ phi \\ \ phi ^ {(R)} & = & {1 \ over mc} \ left (i \ hbar {\ pa ... ... al x_0} -i \ hbar \ vec {\ sigma} \ cdot \ vec {\ nabla} \ right) \ phi ^ {(L)} \\ \ end {eqnarray *}$

Agora reescrevendo em termos de e e ordenando-o como uma equação de matriz, obtemos uma equação que pode ser escrita como um produto escalar entre 4 vetores. $\ bgroup \ color {black} $ \ psi_A = \ phi ^ {(R)} + \ phi ^ {(L)} $ \ egroup$ $\ bgroup \ color {black} $ \ psi_B = \ phi ^ {(R)} - \ phi ^ {(L)} $ \ egroup$

$\ begin {eqnarray *} \ pmatrix {-i \ hbar {\ partial \ over \ partial x_0} & -i \ hbar \ vec {\ sig ... ...] = \ hbar \ left [\ gamma_ \ mu { \ partial \ over \ partial x_ \ mu} \ right] \\ \ end {eqnarray *}$

Defina as matrizes 4 por 4 $\ bgroup \ color {black} $ \ gamma_ \ mu $ \ egroup$ por.

$\ begin {eqnarray *} \ gamma_i & = & \ pmatrix {0 & -i \ sigma_i \ cr i \ sigma_i & 0 \ cr} \\ \ gamma_4 & = & \ pmatrix {1 & 0 \ cr 0 & -1 \ cr} \\ \ end {eqnarray *}$

Com esta definição, a equação relativística pode ser bastante simplificada

$\ begin {eqnarray *} \ left (\ gamma_ \ mu {\ partial \ over \ partial x_ \ mu} + {mc \ over \ hbar} \ right) \ psi = 0 \\ \ end {eqnarray *}$

onde as matrizes gama são dadas por

$\ bgroup \ color {black} $ \ gamma_1 = \ pmatrix {0 & 0 & 0 & -i \ cr 0 & 0 & -i & 0 \ cr 0 & i & 0 & 0 & 0 \ cr i & 0 & 0 & 0 \ cr} $ \ egroup$

$\ bgroup \ color {black} $ \ gamma_2 = \ pmatrix {0 & 0 & 0 & -1 \ cr 0 & 0 & 1 & 0 \ cr 0 & 1 & 0 & 0 & 0 \ cr -1 & 0 & 0 & 0 \ cr} $ \ egroup$ $\ bgroup \ color {black} $ \ gamma_3 = \ pmatrix {0 & 0 & -i & 0 \ cr 0 & 0 & 0 & i \ cr i & 0 & 0 & 0 & 0 & 0 \ cr 0 & -i & 0 & 0 \ cr} $ \ egroup$

e eles satisfazem as relações anti-comutação.

$\ begin {displaymath} \ bgroup \ color {black} \ {\ gamma_ \ mu, \ gamma_ \ nu \} = 2 \ delta _ {\ mu \ nu} \ egroup \ end {displaymath}$

Na verdade, qualquer conjunto de matrizes que satisfaça as relações anticomutação produziria resultados físicos equivalentes, no entanto, trabalharemos na representação explícita acima das matrizes gama.

Definindo , $\ bgroup \ color {black} $ \ bar {\ psi} = \ psi ^ \ dagger \ gamma_4 $ \ egroup$

$\ begin {displaymath} \ bgroup \ color {black} j_ \ mu = ic \ bar {\ psi} \ gamma_ \ mu \ psi \ egroup \ end {displaymath}$

satisfaz a equação de uma corrente conservada de 4 vetores

$\ begin {displaymath} \ bgroup \ color {black} {\ partial \ over \ partial x_ \ mu} j_ \ mu = 0 \ egroup \ end {displaymath}$

e também se transforma como um vetor 4. O quarto componente do vetor mostra que a densidade de probabilidade é . Isso indica que a normalização do estado inclui todos os quatro componentes dos espinores de Dirac. $\ bgroup \ color {black} $ \ psi ^ \ dagger \ psi $ \ egroup$

Para elétrons não relativísticos, os dois primeiros componentes do espinor de Dirac são grandes, enquanto os dois últimos são pequenos.

$\ begin {eqnarray *} \ psi = \ pmatrix {\ psi_A \ cr \ psi_B} \\ \ psi_B \ approx {c \ over 2mc ... ... p} + {e \ over c} \ vec {A} \ direita) \ psi_A \ aprox {pc \ over 2mc ^ 2} \ psi_A \\ \ end {eqnarray *}$

Usamos esse fato para escrever uma equação aproximada de dois componentes derivada da equação de Dirac no limite não relativístico.

$\ begin {eqnarray *} \ left ({p ^ 2 \ over 2m} - {Ze ^ 2 \ over 4 \ pi r} - {p ^ 4 \ over 8m ^ 3c ^ 2} + {Ze ^ .... .. ^ 2 \ over 8m ^ 2c ^ 2} \ delta ^ 3 (\ vec {r}) \ right) \ psi & = & E ^ {(NR)} \ psi \\ \ end {eqnarray *}$

Esta `` equação de Schrödinger '', derivada da equação de Dirac , concorda bem com a que usamos para entender a estrutura fina do Hidrogênio. Os primeiros dois termos são os termos de energia cinética e potencial para o hamiltoniano de hidrogênio não perturbado. O terceiro termo é a correção relativística da energia cinética . O quarto termo é a interação spin-órbita correta , incluindo o efeito Thomas Precession que não perdemos tempo para entender quando fizemos a estrutura fina NR. O quinto termo é o assim chamado termo de Darwin, que dissemos vir da equação de Dirac; e agora está.

Para uma partícula livre, cada componente do spinor de Dirac satisfaz a equação de Klein-Gordon.

$\ begin {displaymath} \ bgroup \ color {black} \ psi _ {\ vec {p}} = u _ {\ vec {p}} e ^ {i (\ vec {p} \ cdot \ vec {x} -Et) / \ hbar} \ egroup \ end {displaymath}$

Isso é consistente com a relação de energia relativística.

As quatro soluções normalizadas para uma partícula de Dirac em repouso são.

$\ begin {eqnarray *} \ psi ^ {(1)} = \ psi_ {E = + mc ^ 2, + \ hbar / 2} & = & {1 \ over \ sqrt {V}} \ pmatrix {... ... \ over \ sqrt {V}} \ pmatrix {0 \ cr 0 \ cr 0 \ cr 1 \ cr} e ^ {+ imc ^ 2t / \ hbar} \\ \ end {eqnarray *}$

O primeiro e o terceiro giram para cima, enquanto o segundo e o quarto giram para baixo. A primeira e a segunda são soluções de energia positiva, enquanto a terceira e a quarta são `` soluções de energia negativa '' , que ainda precisamos entender.

O próximo passo é encontrar as soluções com momentum definido. As quatro soluções de ondas planas para a equação de Dirac são

$\ begin {displaymath} \ bgroup \ color {black} \ psi ^ {(r)} _ {\ vec {p}} \ equiv \ sqrt {mc ^ 2 \ ov ... ... V} u ^ {( r)} _ {\ vec {p}} e ^ {i (\ vec {p} \ cdot \ vec {x} -Et) / \ hbar} \ egroup \ end {displaymath}$

onde os quatro espinores são dados por.

$\ begin {eqnarray *} u ^ {(1)} _ {\ vec {p}} = \ sqrt {E + mc ^ 2 \ over 2mc ^ 2} \ pmatrix {1 \ cr 0 \ cr {.... .._ x-ip_y) c \ over -E + mc ^ 2} \ cr {p_zc \ over -E + mc ^ 2} \ cr 0 \ cr 1 \ cr} \\ \ end {eqnarray *}$

$\ bgroup \ color {black} $ E $ \ egroup$ é positivo para as soluções 1 e 2 e negativo para as soluções 3 e 4. Os espinores são ortogonais

$\ begin {displaymath} \ bgroup \ color {black} u ^ {(r) \ dagger} _ {\ vec {p}} u ^ {(r ')} _ {\ vec {p}} = {\ vert E \ vert \ over mc ^ 2} \ delta_ {rr '} \ egroup \ end {displaymath}$

e as constantes de normalização foram definidas para que os estados sejam normalizados corretamente e os spinors sigam a convenção dada acima, com a normalização proporcional à energia.

As soluções não são em geral estados próprios de qualquer componente do spin, mas são estados próprios de helicidade , o componente do spin ao longo da direção do momento.

Observe que com o $\ bgroup \ color {black} $ E $ \ egroup$ negativo, o exponencial tem a velocidade de fase, a velocidade do grupo e o fluxo de probabilidade, todos na direção oposta do momento como o definimos. Isso claramente não faz sentido. As soluções 3 e 4 precisam ser entendidas de uma forma para a qual os operadores não relativísticos não nos prepararam. Vamos simplesmente renomear as soluções 3 e 4 de forma que $\ bgroup \ color {black} $ e ^ {i (\ vec {p} \ cdot \ vec {x} -Et) / \ hbar} $ \ egroup$

$\ begin {eqnarray *} \ vec {p} \ rightarrow - \ vec {p} \\ E \ rightarrow -E \\ \ end {eqnarray *}$

de forma que todas as energias são positivas e os momentos apontam na direção das velocidades. Isso significa que alteramos os sinais nas soluções 3 e 4 da seguinte maneira.

$\ begin {eqnarray *} \ psi ^ {(1)} _ {\ vec {p}} & = & \ sqrt {E + mc ^ 2 \ over 2EV} \ pmatrix {1 \ cr 0 \ ... .. .er E + mc ^ 2} \ cr 0 \ cr 1 \ cr} e ^ {- i (\ vec {p} \ cdot \ vec {x} -Et) / \ hbar} \\ \ end {eqnarray *}$

Temos ondas planas da forma

$\ begin {displaymath} \ bgroup \ color {black} e ^ {\ pm ip_ \ mu x_ \ mu / \ hbar} \ egroup \ end {displaymath}$

com o sinal de mais para as soluções 1 e 2 e o sinal de menos para as soluções 3 e 4. Esses $\ bgroup \ color {black} $ \ pm $ \ egroup$ sinais no exponencial não são muito surpreendentes do ponto de vista das soluções possíveis para uma equação diferencial. O problema agora é que para as soluções 3 e 4 os operadores de momentum e energia devem ter um sinal negativo adicionado a eles e a fase da função de onda em uma posição fixa se comporta de maneira oposta em função do tempo do que esperamos e de soluções 1 e 2. É como se as soluções 3 e 4 estivessem retrocedendo no tempo.

Se mudarmos a carga do elétron de $\ bgroup \ color {black} $ - e $ \ egroup$ para $\ bgroup \ color {black} $ + e $ \ egroup$ e mudarmos o sinal do expoente, a equação de Dirac permanece invariante. Assim, podemos transformar a solução do expoente negativo (retrocedendo no tempo) na solução do expoente positivo convencional se mudarmos a carga para $\ bgroup \ color {black} $ + e $ \ egroup$ . Podemos interpretar as soluções 3 e 4 como pósitrons. Faremos essa mudança com mais cuidado quando estudarmos o operador de conjugação de carga.

A equação de Dirac deve ser invariante sob impulsos de Lorentz e rotações, sendo que ambos são apenas mudanças na definição de um sistema de coordenadas inercial. Sob aumentos de Lorentz, transforma-se como um vetor 4, mas as matrizes são constantes. A equação de Dirac mostra-se invariante sob impulsos ao longo da direção se transformarmos o espinor de Dirac de acordo com $\ bgroup \ color {black} $ {\ partial \ over \ partial x_ \ mu} $ \ egroup$ $\ bgroup \ color {black} $ \ gamma_ \ mu $ \ egroup$ $\ bgroup \ color {black} $ x_i $ \ egroup$

$\ begin {eqnarray *} \ psi '& = & S_ {boost} \ psi \\ S_ {boost} & = & \ cosh {\ chi \ over 2} + i \ gamma_i \ gamma_4 \ sinh {\ chi \ over 2} \\ \ end {eqnarray *}$

com . $\ bgroup \ color {black} $ \ tanh \ chi = \ beta $ \ egroup$

A equação de Dirac é invariante em rotações em torno do $\ bgroup \ color {black} $ k $ \ egroup$ eixo se transformarmos o espinor de Dirac de acordo com

$\ begin {eqnarray *} \ psi '& = & S_ {rot} \ psi \\ S_ {rot} & = & \ cos {\ theta \ over 2} + \ gamma_i \ gamma_j \ sin {\ theta \ over 2} \ fim de {eqnarray *}$

com $\ bgroup \ color {black} $ ijk $ \ egroup$ é uma permutação cíclica.

Outra simetria relacionada à escolha do sistema de coordenadas é a paridade. Sob uma operação de inversão de paridade, a equação de Dirac permanece invariante se

$\ begin {displaymath} \ bgroup \ color {black} \ psi '= S_P \ psi = \ gamma_4 \ psi \ egroup \ end {displaymath}$

Desde então , o terceiro e o quarto componentes do spinor mudam de sinal, enquanto os dois primeiros não. Como poderíamos ter escolhido , tudo o que sabemos é que os componentes 3 e 4 têm a paridade oposta dos componentes 1 e 2 . $\ bgroup \ color {black} $ \ gamma_4 = \ pmatrix {1 & 0 & 0 & 0 \ cr 0 & 1 & 0 & 0 \ cr 0 & 0 & -1 & 0 \ cr 0 & 0 & 0 & -1 \ cr} $ \ egroup$ $\ bgroup \ color {black} $ - \ gamma_4 $ \ egroup$

De 4 por 4 matrizes, podemos derivar 16 componentes independentes de objetos covariantes. Nós definir o produto de todas as matrizes gama .

$\ begin {displaymath} \ bgroup \ color {black} \ gamma_5 = \ gamma_1 \ gamma_2 \ gamma_3 \ gamma_4 \ egroup \ end {displaymath}$

que obviamente anticomuta com todas as matrizes gama.

$\ begin {displaymath} \ bgroup \ color {black} \ {\ gamma_ \ mu, \ gamma_5 \} = 0 \ egroup \ end {displaymath}$

Para rotações e impulsos, $\ bgroup \ color {black} $ \ gamma_5 $ \ egroup$ comuta com, $\ bgroup \ color {black} $ S $ \ egroup$ pois comuta com o par de matrizes gama. Para uma inversão de paridade, é anti-mudo com $\ bgroup \ color {black} $ S_P = \ gamma_4 $ \ egroup$ .

O conjunto mais simples de covariantes que podemos fazer a partir de espinores e $\ bgroup \ color {black} $ \ gamma $ \ egroup$ matrizes de Dirac estão tabulados abaixo.

Classificação	Forma covariante	não. de componentes

Escalar	$\ bgroup \ color {black} $ \ bar {\ psi} \ psi $ \ egroup$	1
Pseudoescalar	$\ bgroup \ color {black} $ \ bar {\ psi} \ gamma_5 \ psi $ \ egroup$	1
Vetor	$\ bgroup \ color {black} $ \ bar {\ psi} \ gamma_ \ mu \ psi $ \ egroup$	4
Vetor Axial	$\ bgroup \ color {black} $ \ bar {\ psi} \ gamma_5 \ gamma_ \ mu \ psi $ \ egroup$	4
Tensor antissimétrico de Rank 2	$\ bgroup \ color {black} $ \ bar {\ psi} \ sigma _ {\ mu \ nu} \ psi $ \ egroup$	6
Total		16

Produtos de mais $\ bgroup \ color {black} $ \ gamma $ \ egroup$ matrizes acabam repetindo as mesmas quantidades porque o quadrado de qualquer $\ bgroup \ color {black} $ \ gamma $ \ egroup$ matriz é 1.

Para muitos propósitos, é útil escrever a equação de Dirac na forma tradicional $\ bgroup \ color {black} $ H \ psi = E \ psi $ \ egroup$ . Para fazer isso, devemos separar as derivadas de espaço e tempo, tornando a equação menos covariante.

$\ begin {eqnarray *} \ left (\ gamma_ \ mu {\ partial \ over \ partial x_ \ mu} + {mc \ over \ hbar} \ r ... ... mc ^ 2} \ gamma_4 \ right) \ psi = - \ hbar {\ partial \ over \ partial t} \ psi \\ \ end {eqnarray *}$

Assim, podemos identificar o operador abaixo como o Hamiltoniano.

$\ begin {displaymath} \ bgroup \ color {black} H = ic \ gamma_4 \ gamma_jp_j + mc ^ 2 \ gamma_4 \ egroup \ end {displaymath}$

O hamiltoniano nos ajuda a identificar constantes do movimento. Se um operador comuta com $\ bgroup \ color {black} $ H $ \ egroup$ , ele representa uma quantidade conservada.

É fácil ver o $\ bgroup \ color {black} $ p_k $ \ egroup$ trajeto com o hamiltoniano para uma partícula livre, de modo que o momentum seja conservado . Os componentes do momento angular orbital não comutam com $\ bgroup \ color {black} $ H $ \ egroup$ .

$\ begin {displaymath} \ bgroup \ color {black} [H, L_z] = ic \ gamma_4 [\ gamma_jp_j, xp_y-yp_x] = \ hbar c \ gamma_4 (\ gamma_1p_y- \ gamma_2 p_x) \ egroup \ end {displaymath}$

Os componentes do spin também não comutam com $\ bgroup \ color {black} $ H $ \ egroup$ .

$\ begin {displaymath} \ bgroup \ color {black} {[H, S_z]} = \ hbar c \ gamma_4 [\ gamma_2p_x- \ gamma_1p_y] \ egroup \ end {displaymath}$

Mas, do acima exposto, os componentes do momento angular total comutam com $\ bgroup \ color {black} $ H $ \ egroup$ .

$\ begin {eqnarray *} [H, J_z] = [H, L_z] + [H, S_z] = \ hbar c \ gamma_4 (\ gamma_1p_y- \ gamma_2 p_x) + \ hbar c \ gamma_4 [\ gamma_2p_x- \ gamma_1p_y] = 0 \\ \ end {eqnarray *}$

A equação de Dirac naturalmente conserva o momento angular total, mas não as partes orbitais ou de spin dele.

Também podemos ver que a helicidade , ou giro ao longo da direção do movimento, comuta.

$\ begin {displaymath} \ bgroup \ color {black} [H, \ vec {S} \ cdot \ vec {p}] = [H, \ vec {S}] \ cdot \ vec {p} = 0 \ egroup \ fim {displaymath}$

Para qualquer cálculo, precisamos saber o termo de interação com o campo eletromagnético. Com base na interação do campo com uma corrente

$\ begin {displaymath} \ bgroup \ color {black} H_ {int} = - {1 \ over c} j_ \ mu A_ \ mu \ egroup \ end {displaymath}$

e a corrente que encontramos para a equação de Dirac, o hamiltoniano de interação é.

$\ begin {displaymath} \ bgroup \ color {black} H_ {int} = ie \ gamma_4 \ gamma_k A_k \ egroup \ end {displaymath}$

Isso é mais simples do que o caso não relativístico, sem $\ bgroup \ color {black} $ A ^ 2 $ \ egroup$ termo e apenas uma potência de $\ bgroup \ color {black} $ e $ \ egroup$ .

A equação de Dirac tem alguns fenômenos inesperados que podemos derivar. Os autovalores de velocidade para elétrons estão sempre $\ bgroup \ color {black} $ \ pm c $ \ egroup$ ao longo de qualquer direção. Portanto, os únicos valores de velocidade que podemos medir são $\ bgroup \ color {black} $ \ pm c $ \ egroup$ .

Estados localizados, expandidos em ondas planas, contêm todos os quatro componentes das soluções de ondas planas. A mistura dos componentes 1 e 2 com os componentes 3 e 4 dá origem a Zitterbewegung , a oscilação muito rápida da velocidade e posição dos elétrons.

$\ begin {eqnarray *} \ langle v_k \ rangle & = & \ sum \ limits _ {\ vec {p}} \ sum \ limits_ {r = 1} ^ 4 \ ... ... mma_k u ^ {(r ') } _ {\ vec {p}} e ^ {2i \ vert E \ vert t / \ hbar} \ right] \\ \ end {eqnarray *}$

A última soma que contém os termos cruzados entre energia negativa e positiva representa oscilações de frequência extremamente altas no valor esperado da velocidade , conhecido como Zitterbewegung. O valor esperado da posição tem oscilações rápidas semelhantes.

É possível resolver a equação de Dirac exatamente para o hidrogênio de uma forma muito semelhante à solução não relativística. Uma diferença é que está claro desde o início que o momento angular total é uma constante do movimento e é usado como um número quântico básico. Há outro número quântico conservado relacionado ao componente de spin ao longo da direção de $\ bgroup \ color {black} $ \ vec {J} $ \ egroup$ . Com esses números quânticos, a equação radial pode ser resolvida de maneira semelhante ao caso não relativístico que produz a relação de energia .

$\ begin {displaymath} \ bgroup \ color {black} E = {mc ^ 2 \ over \ sqrt {1+ {Z ^ 2 \ alpha ^ 2 \ over \ le ... ... rt {\ left (j + { 1 \ over 2} \ right) ^ 2-Z ^ 2 \ alpha ^ 2} \ right) ^ 2}}} \ egroup \ end {displaymath}$

Podemos identificar o número quântico do princípio padrão neste caso como . Esse resultado fornece a mesma resposta que nosso cálculo não relativístico para ordem, mas também está correto para ordem superior . É uma solução exata para o problema de mecânica quântica proposto, mas não inclui os efeitos da teoria de campos , como o deslocamento de Lamb e o momento magnético anômalo do elétron. $\ bgroup \ color {black} $ n = n_r + j + {1 \ over 2} $ \ egroup$ $\ bgroup \ color {black} $ \ alpha ^ 4 $ \ egroup$

Um cálculo de espalhamento Thomson mostra que mesmo o espalhamento simples de fótons de baixa energia depende da "energia negativa" ou estados de pósitron para obter uma resposta diferente de zero. Se o cálculo for feito com os dois diagramas em que um fóton é absorvido e depois emitido por um elétron (e vice-versa), o resultado é zero em baixa energia porque a interação hamiltoniana conecta o primeiro e o segundo estados de onda plana com o terceiro e o quarto em momento zero. Isso está em contradição com os cálculos clássicos e não relativísticos, bem como com a medição. Existem diagramas adicionais se considerarmos a possibilidade de que o fóton pode criar um par de pósitrons de elétrons que se aniquila com o elétron inicial emitindo um fóton (ou com os fótons inicial e final trocados). Esses dois termos fornecem a resposta certa.

A equação de Dirac é invariante sob a conjugação de carga, definida como a mudança dos estados do elétron para os estados de pósitrons carregados opostos com o mesmo momento e spin (e mudança do sinal dos campos externos). Para fazer isso, o spinor de Dirac é transformado de acordo com.

$\ begin {displaymath} \ bgroup \ color {black} \ psi '= \ gamma_2 \ psi ^ * \ egroup \ end {displaymath}$

É claro que uma segunda operação de conjugação de carga leva o estado de volta ao original $\ bgroup \ color {black} $ \ psi $ \ egroup$ . Aplicar isso às soluções de onda plana dá

$\ begin {eqnarray *} \ psi ^ {(1)} _ {\ vec {p}} = \ sqrt {mc ^ 2 \ over \ vert E \ vert V} u ^ {(1)} _ {... ... {- \ vec {p}} e ^ {i (- \ vec {p} \ cdot \ vec {x} - \ vert E \ vert t) / \ hbar} \\ \ end {eqnarray *}$

que define novos espinores de pósitrons e que são conjugados de carga de e . $\ bgroup \ color {black} $ v ^ {(1)} _ {\ vec {p}} $ \ egroup$ $\ bgroup \ color {black} $ v ^ {(2)} _ {\ vec {p}} $ \ egroup$ $\ bgroup \ color {black} $ u ^ {(1)} _ {\ vec {p}} $ \ egroup$ $\ bgroup \ color {black} $ u ^ {(2)} _ {\ vec {p}} $ \ egroup$

A Equação de Dirac

Para prosseguir em direção a uma teoria de campo para elétrons e quantização do campo de Dirac, desejamos encontrar um Lagrangiano escalar que produza a equação de Dirac. A partir do estudo das covariantes de Lorentz, sabemos que é um escalar e que podemos formar um escalar a partir do produto escalar de dois 4 vetores como no Lagrangiano abaixo. O Lagrangiano não pode depender explicitamente das coordenadas. $\ bgroup \ color {black} $ \ bar {\ psi} \ psi $ \ egroup$

$\ begin {displaymath} \ bgroup \ color {black} {\ cal L} = - c \ hbar \ bar {\ psi} \ gamma_ \ mu {\ partial \ over \ partial x_ \ mu} \ psi-mc ^ 2 \ bar {\ psi} \ psi \ egroup \ end {displaymath}$

(Também poderíamos adicionar um termo tensorial, mas não é necessário para obter a equação de Dirac.) Os campos independentes são considerados os 4 componentes de $\ bgroup \ color {black} $ \ psi $ \ egroup$ e os quatro componentes de $\ bgroup \ color {black} $ \ bar {\ psi} $ \ egroup$ . A equação de Euler-Lagrange usando os $\ bgroup \ color {black} $ \ bar {\ psi} $ \ egroup$ campos independentes é simples, pois não há derivada de $\ bgroup \ color {black} $ \ bar {\ psi} $ \ egroup$ no Lagrange.

$\ begin {eqnarray *} {\ partial \ over \ partial x_ \ mu} \ left ({\ partial {\ cal L} \ over \ part ... ... \ mu {\ partial \ over \ partial x_ \ mu } + {mc \ over \ hbar} \ right) \ psi = 0 \\ \ end {eqnarray *}$

Isso nos dá a equação de Dirac indicando que este Lagrangiano é o correto. A equação de Euler-Lagrange derivada usando os campos $\ bgroup \ color {black} $ \ psi $ \ egroup$ é a equação adjunta de Dirac ,

A densidade Hamiltoniana pode ser derivada da Lagrangiana na forma padrão e a Hamiltoniana total calculada pela integração no espaço. Observe que a densidade hamiltoniana é a mesma que a hamiltoniana derivada da equação de Dirac diretamente.

$\ begin {displaymath} \ bgroup \ color {black} H = \ int \ psi ^ \ dagger \ left (\ hbar c \ gamma_4 \ ... ... \ partial \ over \ partial x_k} + mc ^ 2 \ gamma_4 \ direita) \ psi d ^ 3x \ egroup \ end {displaymath}$

Podemos expandir $\ bgroup \ color {black} $ \ psi $ \ egroup$ em ondas planas para entender o hamiltoniano como uma soma de osciladores.

$\ begin {displaymath} \ bgroup \ color {black} \ psi (\ vec {x}, t) = \ sum \ limits _ {\ vec {p}} \ su ... ... u ^ {(r)} _ {\ vec {p}} e ^ {i (\ vec {p} \ cdot \ vec {x} -Et) / \ hbar} \ egroup \ end {displaymath}$

$\ begin {displaymath} \ bgroup \ color {black} \ psi ^ \ dagger (\ vec {x}, t) = \ sum \ limits _ {\ ve ... ... dagger} _ {\ vec {p}} e ^ {- i (\ vec {p} \ cdot \ vec {x} -Et) / \ hbar} \ egroup \ end {displaymath}$

Escrevendo o hamiltoniano em termos desses campos , a fórmula pode ser simplificada, resultando

$\ begin {eqnarray *} H & = & \ sum \ limits _ {\ vec {p}} \ sum \ limits_ {r = 1} ^ 4 E c _ {\ vec {p}, r} ^ * c _ {\ vec {p }, r}. \\ \ end {eqnarray *}$

Por analogia com o eletromagnetismo, podemos substituir os coeficientes de Fourier para as ondas planas de Dirac por operadores.

$\ begin {eqnarray *} H & = & \ sum \ limits _ {\ vec {p}} \ sum \ limits_ {r = 1} ^ 4 E b ^ {(r) ^ \ dagger ... ... {(r ) \ dagger} _ {\ vec {p}} e ^ {- i (\ vec {p} \ cdot \ vec {x} -Et) / \ hbar} \\ \ end {eqnarray *}$

A criação de uma operadores de aniquilação e Satisfazer as relações anticomutadoras. $\ bgroup \ color {black} $ b ^ {(r) ^ \ dagger} _ {\ vec {p}} $ \ egroup$ $\ bgroup \ color {black} $ b ^ {(r)} _ {\ vec {p}} $ \ egroup$

$\ begin {eqnarray *} \ {b ^ {(r)} _ {\ vec {p}}, b ^ {(r ') ^ \ dagger} _ {\ vec {p'}} \} & = & \ delta_ {rr ... ... {(r)} _ {\ vec {p}} & = & b ^ {(r) ^ \ dagger} _ {\ vec {p}} b ^ {(r)} _ {\ vec {p}} \\ \ end {eqnarray *}$

$\ bgroup \ color {black} $ N ^ {(r)} _ {\ vec {p}} $ \ egroup$ é o operador do número de ocupação. As relações anti-comutação restringem o número de ocupação a 1 ou 0 .

O campo de Dirac e o hamiltoniano podem agora ser reescritos em termos de campos de elétron e pósitron para os quais a energia é sempre positiva, substituindo o operador para aniquilar um `` estado de energia negativa '' por um operador para criar um estado de pósitron com o momento certo e rodar.

$\ begin {eqnarray *} d ^ {(1)} _ {\ vec {p}} & = & - b ^ {(4) ^ \ dagger} _ {\ vec {p}} \\ d ^ {(2 )} _ {\ vec {p}} & = & b ^ {(3) ^ \ dagger} _ {\ vec {p}} \\ \ end {eqnarray *}$

Estes anti-comutação com tudo o mais, exceto que

$\ begin {displaymath} \ bgroup \ color {black} \ {d ^ {(s)} _ {\ vec {p}}, d ^ {(s ') ^ \ dagger} _ {\ vec {p'}} \} = \ delta_ {ss '} \ delta _ {\ vec {p} \ vec {p'}} \ egroup \ end {displaymath}$

Agora reescreva os campos e o hamiltoniano.

$\ begin {eqnarray *} \ psi (\ vec {x}, t) & = & \ sum \ limits _ {\ vec {p}} \ sum \ limits_ {s = 1} ^ 2 \ sqr ... ... vec {p}} + d ^ {(s) ^ \ dagger} _ {\ vec {p}} d ^ {(s)} _ {\ vec {p}} - 1 \ right) \\ \ end {eqnarray *}$

Todas as energias desses estados são positivas .

Existe uma energia constante (infinita), semelhante mas de sinal oposto ao do campo EM quantizado, que devemos adicionar para fazer com que o estado de vácuo tenha energia zero. Observe que, se tivéssemos usado operadores de deslocamento (Bose-Einstein) em vez de anti-deslocamento, não haveria estado fundamental de energia mais baixa, de modo que essa subtração de energia não teria sido possível. As estatísticas de Fermi-Dirac são necessárias para partículas que satisfaçam a equação de Dirac .

Uma vez que os operadores que criam estados de férmions são anti-comutação , os estados de férmions devem ser anti-simétricos sob o intercâmbio. Assuma $\ bgroup \ color {black} $ b_r ^ \ dagger $ \ egroup$ e $\ bgroup \ color {black} $ b_r $ \ egroup$ são os operadores de criação e aniquilação de férmions e que eles anti-comutam.

$\ begin {displaymath} \ bgroup \ color {black} \ {b_r ^ \ dagger, b_ {r '} ^ \ dagger \} = 0 \ egroup \ end {displaymath}$

Os estados são então anti-simétricos sob intercâmbio de pares de férmions.

$\ begin {displaymath} \ bgroup \ color {black} b_r ^ \ dagger b_ {r '} ^ \ dagger \ vert \ rangle = -b_ {r'} ^ \ dagger b_ {r} ^ \ dagger \ vert \ rangle \ egroup \ end {displaymath}$

Não é difícil mostrar que o número de ocupação para estados de férmions é zero ou um .

Observe que os spinors satisfazem as seguintes equações ligeiramente diferentes .

$\ begin {eqnarray *} (i \ gamma_ \ mu p_ \ mu + mc) u ^ {(s)} _ {\ vec {p}} = 0 \\ (-i \ gamma_ \ mu p_ \ mu + mc) v ^ {(s)} _ {\ vec {p}} = 0 \\ \ end {eqnarray *}$

Sumário

Problemas com Física Clássica

Pensamentos sobre experimentos de difração

Amplitudes de probabilidade

Pacotes de onda e incerteza

Operadores

Expectativas de valores

Comutadores

A Equação de Schrödinger

Autofunções, autovalores e espaços vetoriais

Uma partícula em uma caixa

Potenciais constantes por partes em uma dimensão

O oscilador harmônico em uma dimensão

Potenciais de função delta em uma dimensão

Solução de oscilador harmônico com operadores

Mais diversão com os operadores

Duas partículas em 3 dimensões

Partículas Idênticas

Alguns problemas 3D separáveis ​​em coordenadas cartesianas

Momentum Angular

Soluções para a equação radial para potenciais constantes

Hidrogênio

Solução do problema 3D do OH em coordenadas esféricas

Representação Matricial de Operadores e Estados

Um estudo de operadores e autofunções

Spin 1/2 e outros 2 sistemas de estado

Mecânica Quântica em um Campo Eletromagnético

Simetria de fase local em mecânica quântica e a simetria de calibre

Adição de Momentum Angular

Teoria da Perturbação Independente do Tempo

A Estrutura Fina do Hidrogênio

Estrutura Hiperfina

O átomo de hélio

Física Atômica

Moléculas

Teoria da Perturbação Dependente do Tempo

Radiação em átomos

Teoria de Campo Clássica

O Campo Eletromagnético Clássico

Quantização do campo EM

Dispersão de fótons

Electron Self Energy

A Equação de Dirac

A Equação de Dirac

História

Estudos de Química

Alguns problemas 3D separáveis em coordenadas cartesianas

Um estudo de operadores $\ bgroup \ color {black} $ \ ell = 1 $ \ egroup$ e autofunções