Prefácio

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} $ \ egroupem 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 $ \ egrouptodos 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) $ \ egroupsã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) $ \ egroupou em espaço de momento com \ bgroup \ color {black} $ \ phi (p) $ \ egroup.

Podemos usar \ bgroup \ color {black} $ \ phi (p) $ \ egrouppara 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 $ \ egroupvenha a atuar em nossa função de onda simples, obtenhamos \ bgroup \ color {black} $ v $ \ egroupvezes 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) $ \ egroupatuar

\ 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 $ \ egroupque é 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}

e

\ 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 $ \ egroupcom \ bgroup \ color {black} $ H $ \ egroup, obteremos a mesma função \ bgroup \ color {black} $ \ psi_i $ \ egroupvezes alguma constante. Nesse caso, \ bgroup \ color {black} $ \ psi_i $ \ egroupseria chamada de função autônoma e \ bgroup \ color {black} $ E_i $ \ egroupseria chamada de autovalor. Normalmente, há um número infinito de soluções, indicado pelo índice \ bgroup \ color {black} $ i $ \ egroupaqui.

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 $ \ egroupsã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 $ \ egroupforma é 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 $ \ egroupseria 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 $ \ egroupem \ bgroup \ color {black} $ t = 0 $ \ egroup.

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

\ 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 $ \ egrouptê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 $ \ egroupe nenhuma das outras é. Para tornar a função seno zero em \ bgroup \ color {black} $ x = a $ \ egroupprecisamos \ bgroup \ color {black} $ ka = n \ pi $ \ egroupou \ 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 $ \ egroupporque satisfaz à equação diferencial. Voltando a \ bgroup \ color {black} $ \ sin \ left ({n \ pi x \ over a} \ right) $ \ egroupna 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}

ou

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

ou

\ 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\ bgroup \ color {black} $ \ kappa $ \ egroupsã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 $ \ egrouprefletir-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 $ \ egrouppara \ bgroup \ color {black} $ - a <x <a $ \ egroupe \ bgroup \ color {black} $ V (x) = 0 $ \ egroupem 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 $ \ egrouppara \ bgroup \ color {black} $ - a <x <a $ \ egroupe \ bgroup \ color {black} $ V (x) = 0 $ \ egroupem 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} $ \ egrouppara\ bgroup \ color {black} $ x <0 $ \ egroup e \ bgroup \ color {black} $ Ce ^ {- \ kappa x} $ \ egrouppara \ bgroup \ color {black} $ x> 0 $ \ egroup. Tornar \ bgroup \ color {black} $ u (x) $ \ egroupcontí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 $ \ egroupo 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 $ \ egroupcom 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 $ \ egroupaumenta 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}

e

\ 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\ bgroup \ color {black} $ x $ \ egrouppodem ser escritos em termos de \ bgroup \ color {black} $ A $ \ egroup\ 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 $ \ egroupcomo 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\ 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\ 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 $ \ egroupque 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 $ \ egroupque 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} $ \ egrouppartí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 $ \ egroupo 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\ 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}

equação de Schrödinger, portanto, separa em uma parte angular (o \ bgroup \ color {black} $ L ^ 2 $ \ egrouptermo) 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) $ \ egroupserá 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) $ \ egroupque 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\ bgroup \ color {black} $ m $ \ egroupinteiros 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$ z $

\ 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 $ \ egrouppotenciais 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 $ \ egroupem geral, em seguida, encontrando o comportamento de\ bgroup \ color {black} $ r $ \ egrouppequeno , 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}$\egroupe 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 $ \ egroupe 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 $ \ egroupestados para \ bgroup \ color {black} $ \ ell = 1 $ \ egroupserã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 $ \ egroupsã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 $ \ egroupmas é 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 $ \ egroupessas 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 $ \ egroupque 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 $ \ egroupsã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 $ \ egroupe da identidade.

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

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} $ \ egroupe 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 $ \ egrouppode 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 $ \ egroupcampo, mudaremos o padrão de difração , embora \ bgroup \ color {black} $ B = 0 $ \ egroupao 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\ bgroup \ color {black} $ \ vec {L_2} $ \ egrouppara 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\ 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 $ \ egroupestados, obtenho estados de momento angular total com \ bgroup \ color {black} $ j = 0,1,2,3 $ \ egroupe 4. Existem 25 estados de produto, pois cada \ bgroup \ color {black} $ \ ell = 2 $ \ egroupestado tem 5 diferentes \ bgroup \ color {black} $ m $ \ egroups 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 $ \ egroupexiste 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 $ \ egroupcomparação com as energias de hidrogênio que calculamos.

  1. A correção relativística da energia cinética do elétron.
  2. A correção Spin-Orbit.
  3. 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$ g $\ 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\ bgroup \ color {black} $ m_s $ \ egroupe calcular os efeitos da estrutura fina, \ bgroup \ color {black} $ H_1 + H_2 $ \ egroupcomo 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

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 $ \ egroupcomo 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} $ \ egroupestá 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 $ \ egroupestados. 

\ 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\ bgroup \ color {black} $ m_f $ \ egroupque 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} $ \ egrouppartí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 $ \ egroupcm.

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.

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 $ \ egroupestado 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.

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} $ \ egroupcomportamento 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 $ \ egroupestados 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 $ \ egroupestados 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

  1. Acople os elétrons de valência (ou lacunas) para dar o spin total máximo .
  2. 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).
  3. 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

$ \ 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 $ \ egroupsistema 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) $ \ egrouppotencial 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

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 $ \ egrouptermo 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} $ \ egrouppor \ bgroup \ color {black} $ \ sqrt {N + 1} $ \ egroupno 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 $ \ egroupcampo 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 $ \ egroupvai 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 $ \ egrouprepresenta 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 $ \ egrouppara obter a \ bgroup \ color {black} $ - $ \ egroupno 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 $ \ egrouppode 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} $ \ egroupUm 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} $ \ egroupe 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 $ \ egroupcomo 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\ 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 $ \ egroupe 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 $ \ egroupno 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\ bgroup \ color {black} $ N + 1 $ \ egrouppara 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 $ \ egrouptermo 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 $ \ egrouptermo para ir de um fóton de vetor de onda \ bgroup \ color {black} $ \ vec {k} $ \ egroupe um estado atômico \ bgroup \ color {black} $ i $ \ egroupa um fóton espalhado de vetor de onda \ bgroup \ color {black} $ \ vec {k '} $ \ egroupe 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} $ \ egroupoperador.

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 $ \ egroupeV. 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 $ \ egroupe estados intermediários \ bgroup \ color {black} $ \ psi_j $ \ egrouppodem 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}

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 $ m_ec ^ 2 $, 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} $$ 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} $ \ egroupcom 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 *}

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 $ \ egroupnegativo, 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 $ \ egroupsinais 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 $ \ egrouppara \ bgroup \ color {black} $ + e $ \ egroupe 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

equação de Dirac é invariante em rotações em torno do \ bgroup \ color {black} $ k $ \ egroupeixo 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 $ \ egroupcomuta com, \ bgroup \ color {black} $ S $ \ egrouppois 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 $ \ egroupmatrizes 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 $ \ egroupmatrizes acabam repetindo as mesmas quantidades porque o quadrado de qualquer \ bgroup \ color {black} $ \ gamma $ \ egroupmatriz é 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 $ \ egrouptrajeto 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 $ \ egrouptermo 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 $ \ egroupao 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 $ \ egroupe os quatro componentes de \ bgroup \ color {black} $ \ bar {\ psi} $ \ egroup. A equação de Euler-Lagrange usando os \ bgroup \ color {black} $ \ bar {\ psi} $ \ egroupcampos independentes é simples, pois não há derivada de \ bgroup \ color {black} $ \ bar {\ psi} $ \ egroupno 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 $ \ egroupem 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 *}

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\ bgroup \ color {black} $ b_r $ \ egroupsã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 *}
 

Feynman Lectures, Volume III, Capítulo I.

Gasciorawicz realmente não faz isso.

Capítulos 5 de Rohlf

Griffiths não cobre isso.

Difração de Duas Fendas

As ondas de água exibirão um padrão de interferência difrativa em um experimento de 2 fendas, conforme diagramado abaixo. O diagrama mostra as cristas das ondas de água em algum momento. A jusante das fendas, veremos a interferência construtiva onde as ondas das fendas estão em fase e a interferência destrutiva onde estão 180 graus defasados, por exemplo, onde uma crista de uma fenda encontra uma depressão da outra. O gráfico marcado \ bgroup \ color {black} $ I_ {12} $ \ egroupmostra o padrão de interferência no local do absorvedor. A difração é um fenômeno de onda simples .

 

\ epsfig {file = figs / 2slit_water.eps, height = 3in}

difração da luz era um fenômeno bem conhecido no final do século 19 e era bem explicado na teoria eletromagnética clássica, uma vez que a luz era considerada uma onda EM. O padrão de difração de duas fendas estreitas é particularmente fácil de entender em termos de ondas. A configuração é mostrada no diagrama abaixo.

 

\ epsfig {file = figs / 2slit_diag.eps, height = 3in}

As ondas EM de comprimento de onda \ bgroup \ color {black} $ \ lambda $ \ egroupsão emitidas por uma única fonte de luz, como um laser. Eles viajam para duas fendas estreitas, (para simplificar) equidistantes da fonte e a uma distância uma da \ bgroup \ color {black} $ d $ \ egroupoutra. A luz viaja das fendas para uma tela de detecção. Um padrão de difração pode ser visto na tela de detecção, como o mostrado na imagem abaixo.

 

\ epsfig {file = figs / 2slit_phot.eps, height = 3in}

O centro do padrão de difração ocorre no local da tela equidistante de cada fenda onde as ondas das duas fendas estão em fase (porque elas viajaram exatamente a mesma distância) e os campos se somam, então as ondas interferem construtivamente e há um máximo de intensidade. A alguma distância deste centro do padrão de difração, haverá interferência destrutiva entre as ondas das duas fendas e a intensidade será zero. Isso ocorrerá quando a distância percorrida por duas ondas for diferente em \ bgroup \ color {black} $ \ lambda / 2 $ \ egroup, então as ondas estão 180 graus fora de fase e os campos das duas fendas se cancelam .

Podemos calcular esta localização observando o diagrama acima. Assumimos que a distância até a tela é muito maior do que \ bgroup \ color {black} $ d $ \ egroup. Para luz detectada em um ângulo \ bgroup \ color {black} $ \ theta $ \ egroup, a distância extra percorrida desde a fenda 1 é justa \ bgroup \ color {black} $ d \ sin \ theta $ \ egroup. Portanto, o ângulo do primeiro mínimo (ou nulo) pode ser encontrado na equação . \ bgroup \ color {black} $ d \ sin \ theta = {\ lambda \ over 2} $ \ egroup

De maneira mais geral, obteremos um máximo se os caminhos das fendas diferirem em um número inteiro de comprimentos de onda e obteremos um valor nulo quando os caminhos diferirem em meio número inteiro de comprimentos de onda . . \ bgroup \ color {black} $ d \ sin \ theta = {n \ lambda} $ \ egroup \ bgroup \ color {black} $ d \ sin \ theta_ {null} = {(n + 1) \ lambda \ over 2} $ \ egroup

Embora seja muito difícil porque os elétrons estão carregados, a difração de elétrons de 2 fendas também foi observada.

\ epsfig {file = figs / 2slit_e.eps, height = 4in}

Assim, todos os tipos de partículas parecem difratar, indicando que há algum tipo de onda envolvida. Vamos agora continuar com alguns experimentos mentais sobre difração para ilustrar a física que a Mecânica Quântica precisava corresponder.Difração de fenda única

Existem muitos outros exemplos de difração. A imagem abaixo mostra a difração de uma única fenda onde ondas de diferentes partes da fenda interferem umas nas outras. Para obter o campo em algum ponto da tela de detecção, deve-se integrar sobre a fenda.

 

Exemplo: derivar a localização dos nós no padrão de difração de uma fenda de largura $ a $. Agora tente calcular a distribuição de intensidade para difração de fenda única. *

Exemplo: derivar a localização dos nós no padrão de difração de duas fendas estreitas a uma distância uma da $ d $outra. Agora tente calcular a distribuição de intensidade. *

Difração de fenda única

Existem muitos outros exemplos de difração. A imagem abaixo mostra a difração de uma única fenda onde ondas de diferentes partes da fenda interferem umas nas outras. Para obter o campo em algum ponto da tela de detecção, deve-se integrar sobre a fenda.

Exemplo: derivar a localização dos nós no padrão de difração de uma fenda de largura $ a $. Agora tente calcular a distribuição de intensidade para difração de fenda única.

Difração de Cristais

As ondas de elétrons foram demonstradas pela primeira vez medindo a difração de cristais . Davison e Germer observaram a difração de elétrons de um cristal de níquel em 1928. Eles variaram a energia do elétron para medir o comprimento de onda do elétron, concordando bem com a expectativa de deBroglie.

Primeiro, vemos a difração de um único cristal .

\ epsfig {file = figs / single_crystal_e_diff.eps, height = 4in}

A difração de material policristalino dá anéis concêntricos em vez de pontos.

\ epsfig {file = figs / polycrystaline.eps, height = 4in}

A difração de cristais é uma ferramenta poderosa. Primeiro, vemos a difração de raios-X de um único cristal de sódio que tem localizações periódicas dos átomos.

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

Outro exemplo de difração é mostrado abaixo. Aqui os nêutrons difratam de um único cristal de sódio .

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

O comprimento de onda DeBroglie

A transformação de Lorentz foi postulada para ondas eletromagnéticas antes de Einstein desenvolver a relatividade especial. As ondas EM eram inteiramente consistentes com a Relatividade. Por exemplo, a fase de uma onda EM em algum ponto é a mesma do ponto transformado de Lorentz.

DeBroglie aplicou este requisito de invariância de Lorentz na fase das ondas de matéria para determinar qual deve ser o comprimento de onda. É fácil para nós derivar o comprimento de onda usando 4 vetores. Posição e tempo formam um vetor 4.

 

\ begin {displaymath} \ bgroup \ color {black} x ^ \ mu = (ct, \ vec {x}) \ egroup \ end {displaymath}

 

Energia e impulso formam outro.

 

\ begin {displaymath} \ bgroup \ color {black} p ^ \ mu = (E, \ vec {p} c) \ egroup \ end {displaymath}

 

Lembre-se de que os vetores de Lorentz devem ser transformados, mas os escalares de Lorentz são automaticamente invariantes sob transformações. Por exemplo, o escalar formado pontuando o momento 4 em si mesmo é

 

\ begin {displaymath} \ bgroup \ color {black} p ^ \ mu p_ \ mu = -E ^ 2 + p ^ 2c ^ 2 = -m ^ 2c ^ 4. \ egroup \ end {displaymath}

 

A massa de uma partícula é claramente invariante de Lorentz.

Para calcular o comprimento de onda de nossas ondas de matéria, vamos usar o escalar

 

\ begin {displaymath} \ bgroup \ color {preto} p ^ \ mu x_ \ mu / c = -Et + \ vec {p} \ cdot \ vec {x} \ egroup \ end {displaymath}

 

Já é fácil ver que a fase de uma onda como

 

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

 

é invariante de Lorentz.

Para ler o comprimento de onda, vamos escolher a parte da expressão que corresponde a . Vemos isso e, portanto, o comprimento de onda DeBroglie é. \ bgroup \ color {black} $ e ^ {2 \ pi i (x / \ lambda- \ nu t)} $ \ egroup \ bgroup \ color {black} $ {2 \ pi \ over \ lambda} = {p \ over \ hbar} $ \ egroup

\ bgroup \ color {black} $ \ displaystyle \ lambda = {2 \ pi \ hbar \ over p} = {h \ over p} $ \ egroup
DeBroglie derivou esse resultado em 1923.

Também podemos ler a frequência \ bgroup \ color {black} $ \ nu $ \ egroup.

 

\ begin {displaymath} \ bgroup \ color {black} \ nu = {E \ over 2 \ pi \ hbar} = E / h \ egroup \ end {displaymath}

 

Em certo sentido, essa foi a entrada para o nosso cálculo.

O comprimento de onda deBroglie será nosso principal input de física para o desenvolvimento da Mecânica Quântica. Não é que este trabalho tenha sido o mais significativo, mas este comprimento de onda resume a maior parte do que aconteceu antes de 1923.

Computando comprimentos de onda DeBroglie

Costumamos citar a energia de uma partícula em termos de sua energia cinética em Volts, eV do elétron (ou Volts do Milhão de elétrons, MeV). A razão para isso é que as partículas são geralmente aceleradas para alguma energia por um campo elétrico. Se eu deixar um elétron (ou próton ...) ser acelerado por uma diferença de potencial de 100 Volts, ele terá uma energia cinética de 100eV.

Todo o problema de calcular um comprimento de onda deBroglie é converter a energia cinética em momento. Se você sempre quer estar correto sem nenhuma necessidade de pensar, use a fórmula relativisticamente correta para a energia cinética

 

\ begin {displaymath} \ bgroup \ color {black} T = \ sqrt {(mc ^ 2) ^ 2 + p ^ 2c ^ 2} -mc ^ 2 \ egroup \ end {displaymath}

 

 

e resolvê-lo para \ bgroup \ color {black} $ pc $ \ egroup,

 

\ begin {displaymath} \ bgroup \ color {black} pc = \ sqrt {(T + mc ^ 2) ^ 2- (mc ^ 2) ^ 2} \ egroup \ end {displaymath}

 

 

em seguida, use esta fórmula útil para obter a resposta.

 

\ begin {displaymath} \ bgroup \ color {black} \ lambda = {h \ over p} = {2 \ pi \ hbar c \ over pc} \ egroup \ end {displaymath}

 

 

Lembro-me disso permitindo-me manter todo o cálculo em eV. Eu também conheço as massas das partículas. \ bgroup \ color {black} $ \ hbar c = 1973 \ mathrm {eV \ AA} $ \ egroup

 

\ begin {displaymath} \ bgroup \ color {black} m_ec ^ 2 = 0,51 \ mathrm {MeV} \ egroup \ end {displaymath}

 

 

 

 

\ begin {displaymath} \ bgroup \ color {black} m_pc ^ 2 = 938.3 \ mathrm {MeV} \ egroup \ end {displaymath}

(Se \ bgroup \ color {black} $ T «<mc ^ 2 $ \ egroup, certifique-se de que a precisão de sua calculadora é suficiente ou use o método não relativístico abaixo.)

Se você sabe que a partícula é super-relativista , de modo que \ bgroup \ color {black} $ T »mc ^ 2 $ \ egroup, em seguida, basta usar \ bgroup \ color {black} $ pc = T $ \ egroupea vida é fácil.

Se você sabe que a partícula é altamente não-relativista , \ bgroup \ color {black} $ T «mc ^ 2 $ \ egroup, então você pode usar doação . \ bgroup \ color {black} $ T = {p ^ 2 \ over 2m} = {(pc) ^ 2 \ over 2mc ^ 2} $ \ egroup \ bgroup \ color {black} $ pc = \ sqrt {2mc ^ 2T} $ \ egroup

Então, por exemplo, calcule o comprimento de onda de um elétron de 100 eV. Isso é não relativístico, pois 100 eV «510000 eV. Então eV ou 10000 eV. \ bgroup \ color {black} $ pc = \ sqrt {10 ^ 6 \ times 100} $ \ egroup

\ begin {displaymath} \ bgroup \ color {black} \ lambda = {2 \ pi 1973 \ over pc} = {12000 \ over 10000} = 1,2 \ mathrm {\ AA} \ egroup \ end {displaymath}

Dualidade de partícula de onda (experimentos de pensamento)

Richard Feynman (Prêmio Nobel de Eletrodinâmica Quântica ...) apresenta vários experimentos mentais em suas Lectures on Physics, terceiro volume.

 

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

Para nosso primeiro experimento mental, consideraremos duas difração de luz de silte . Suponha que, em vez de usar a tela, eu use um fotodetector sensível . Medo a intensidade da luz em função da posição perto da tela e encontro o mesmo conjunto de máximos e mínimos que fiz usando meus olhos.

 

\ epsfig {file = figs / 2slit_phot.eps, height = 3.5in}

Agora vamos diminuir a intensidade da fonte de luz. Para intensidade muito baixa, acho que meu detector coleta um fóton por vez . Nunca coleta meio fóton. (Com o detector certo, eu poderia verificar novamente que, para cada fóton, o efeito fotoelétrico é visto e tal \ bgroup \ color {black} $ E = h \ nu $ \ egroup.) Então as ondas que estão difratando são de alguma forma formadas por fótons. Com uma intensidade suficientemente baixa , posso garantir que apenas um fóton está presente no aparelho por vez. Posso operar meu detector e coletar dados por um longo tempo, somando o número de fótons detectados em função da posição. O que vou receber pela distribuição? Eu obtenho exatamente a mesma distribuição de antes com máximos e mínimos. Não importa quão baixa seja a intensidade, (1 partícula / minuto!)ainda vemos difração . Nunca detectamos uma fração de um elétron ou um fóton, apenas números inteiros .

Como um único fóton interfere em si mesmo? Deve de alguma forma passar por ambas as fendas.

Vamos voltar para a difração de elétrons por um minuto. Em nosso experimento mental, novamente temos duas fendas.

 

\ epsfig {file = figs / 2slit_electron.eps, height = 3.2in}

Usamos nosso detector para medir o padrão de difração semelhante ao dos fótons rotulados \ bgroup \ color {black} $ P_ {12} $ \ egroupabaixo. Se cobrirmos a fenda 2 , obteremos a distribuição rotulada \ bgroup \ color {black} $ P_1 $ \ egroupe se / emphtextcovermos a fenda 1 obtermos a distribuição de intensidade rotulada \ bgroup \ color {black} $ P_2 $ \ egroup. Poderíamos ter feito o mesmo com nossos fótons acima.

Agora tentaremos ver por qual fenda cada elétron passa. Colocamos uma fonte de luz brilhante perto das fendas e detectamos a luz refletida no elétron para que possamos ver por qual fenda ele passa.

 

\ epsfig {file = figs / 2slit_electron_see.eps, height = 4in}

Que distribuição vemos agora? Na verdade, veremos \ bgroup \ color {black} $ P_1 + P_2 $ \ egroupse podemos dizer por qual fenda o elétron passou. Nossa observação do elétron ao passar pela fenda mudou a distribuição de intensidade resultante . Se desligarmos a luz, voltamos a medir \ bgroup \ color {black} $ P_ {12} $ \ egroup.

Você pode explicar por que a luz faz com que o padrão de difração desapareça?
É a mera observação? A luz muda a fase do elétron?

Existem muitos exemplos de um observador alterando o resultado de um experimento Quantum. Na verdade, é sustentado que quando um estado é observado, sua função de onda colapsa no estado `` visto ''. Nesse caso, bastava acender a luz. Não precisamos olhar.

Finalmente, faremos um experimento de difração de duas fendas com balas . Devemos fazer fendas grandes o suficiente para as balas passarem.

 

\ epsfig {file = figs / 2slit_bullets.eps, height = 3.0in}

Não importa a distância entre as fendas que escolhemos, nunca observamos difração para as balas. Eles sempre fornecem o \ bgroup \ color {black} $ P_1 + P_2 $ \ egrouppadrão (provavelmente diferente para os marcadores). Você pode explicar isso?
O comprimento de onda da bala é muito menor do que o tamanho real da bala.

Por que não vimos difração para as balas. As balas são objetos macroscópicos com grandes momentos em comparação com partículas individuais. O momento elevado implica um pequeno comprimento de onda deBroglie. Portanto, para ver a difração, devemos tornar a distância entre as fendas muito menor do que fizemos para os fótons (talvez \ bgroup \ color {black} $ 10 ^ {20} $ \ egroupvezes menor). Mas as balas também são grandes. Eles só passam por grandes fendas, que devem estar mais afastadas do que as fendas usadas para os fótons. Para balas, o comprimento de onda é minúsculo em comparação com qualquer fenda pela qual uma bala passaria . Portanto, nenhuma difração de bala é possível.

Feynman Lectures on Physics, Vol. III Capítulo 1

Exemplos

Distribuição de intensidade para difração de duas fendas *

Derive a localização dos nós no padrão de difração a partir de duas fendas estreitas \ bgroup \ color {black} $ d $ \ egroupseparadas por uma distância . Agora tente calcular a distribuição de intensidade.

Este é um exercício de laboratório .

 

Distribuição de intensidade para difração de fenda única *

Derive a localização dos nós no padrão de difração de uma fenda de largura \ bgroup \ color {black} $ a $ \ egroup. Agora tente calcular a distribuição de intensidade para difração de fenda única.

Este é um exercício de laboratório .

Problemas de teste de amostra

  1. Qual é o comprimento de onda deBroglie de um elétron com 13,6 eV de energia cinética? Qual é o comprimento de onda deBroglie de um elétron com 10 MeV de energia cinética?
    A resposta
    13.6 eV é muito menor, então isso não é relativístico. $ mc ^ 2 = 0,511 MeV $

     

    \ begin {eqnarray *} {p ^ 2 \ over 2m} & = & 13.6 \\ {p ^ 2c ^ 2 \ over 2mc ^ 2} & = & 13.6 \\ pc & = & \ ... ... er \ sqrt {2 (0,511 \ vezes 10 ^ 6 eV) (13,6 eV)}} \ aproximadamente 3,33 \ AA \\ \ end {eqnarray *}

     


     

    10 MeV é muito maior do que $ mc ^ 2 $um elétron, então é super-relativístico e podemos usar $ E = pc $.

     

    \ begin {displaymath} \ lambda = {2 \ pi \ hbar c \ over pc} = {2 \ pi (197,3) MeV F) \ over 10 MeV} \ aproximadamente 120 F \ end {displaymath}

     

  2. Qual é o comprimento de onda DeBroglie para cada uma das seguintes partículas? As energias fornecidas são as energias cinéticas.
    a)
    um elétron de 10 eV
    b)
    um elétron de 1 MeV
    c)
    um próton de 10 MeV
  3. Um experimento de difração de elétrons de 2 fendas é configurado como (não) mostrado abaixo. A distribuição de intensidade de elétrons observada é plotada na figura. Agora, uma fonte de luz intensa é introduzida perto das duas fendas. Com esta luz, pode-se "ver" por qual fenda cada elétron passa. Faça um gráfico qualitativo da nova distribuição de intensidade do elétron de cada fenda e das 2 fendas combinadas. Qual é a condição do comprimento de onda da luz para que esse efeito ocorra?
  4. Qual é o comprimento de onda DeBroglie para cada uma das seguintes partículas? As energias fornecidas são as energias cinéticas.
    a)
    um elétron de 1 eV
    b)
    um $ 10 ^ 4 $próton MeV
  5. Qual KE deve ter um átomo de hidrogênio para que seu comprimento de onda DeBroglie seja menor que o tamanho do átomo? (Fatores de 2 não são importantes.)
  6. Calcule o comprimento de onda DeBroglie para (a) um próton com energia cinética de 10 MeV, (b) Um elétron com energia cinética de 10 MeV e (c) uma bola de chumbo de 1 grama movendo-se com uma velocidade de 10 cm / s (um erg é um grama cm $ ^ 2 $/ seg $ ^ 2 $). Certifique-se de levar em conta a relatividade quando necessário. 

 

No final do século XIX, as leis da física baseavam-se na mecânica e na lei da gravitação de Newton, nas equações de Maxwell que descreviam a eletricidade e o magnetismo e na mecânica estatística que descrevia o estado de uma grande coleção de matéria. Essas leis da física descreveram a natureza muito bem na maioria das condições, no entanto, algumas medições do final do século 19 e início do século 20 não puderam ser compreendidas. Os problemas com a física clássica levaram ao desenvolvimento da Mecânica Quântica e da Relatividade Especial.

Alguns dos problemas que levaram ao desenvolvimento da Mecânica Quântica estão listados aqui.

  • Radiação de corpo negro: a física clássica previu que objetos quentes irradiariam instantaneamente todo o seu calor em ondas eletromagnéticas. O cálculo, que foi baseado nas equações de Maxwell e na Mecânica Estatística, mostrou que a taxa de radiação foi para o infinito quando o comprimento de onda EM foi para zero, "A Catástrofe Ultravioleta". Plank resolveu o problema postulando que a energia EM foi emitida em quanta com $ E = h \ nu $.
  • O efeito fotoelétrico: quando a luz foi usada para tirar os elétrons dos sólidos, os resultados foram completamente diferentes do que se esperava das equações de Maxwell. As medições eram fáceis de explicar (para Einstein) se a luz fosse composta de partículas com as energias postuladas por Plank.
  • Átomos: Depois que Rutherford descobriu que a carga positiva nos átomos estava concentrada em um núcleo muito pequeno, a física clássica previu que os elétrons atômicos orbitando o núcleo irradiariam sua energia e entrariam em espiral no núcleo. Isso claramente não aconteceu. A energia irradiada pelos átomos também saiu em quantidades quantizadas em contradição com as previsões da física clássica. O átomo de Bohr postulou uma regra de quantização de momento angular, $ L = n \ hbar $pois $ n = 1,2,3 ... $, isso deu o resultado certo para o hidrogênio, mas acabou se revelando errado, uma vez que o estado fundamental do hidrogênio tem momento angular zero. Foi necessário um conhecimento completo da Mecânica Quântica para explicar os espectros de energia atômica.
  • Espalhamento Compton: Quando a luz foi espalhada pelos elétrons, ela se comportou como uma partícula, mas muda o comprimento de onda no espalhamento; mais evidências para a natureza das partículas da luz e o postulado de Plank.
  • Ondas e partículas: em experimentos de difração, a luz mostrou se comportar como uma onda, enquanto em experimentos como o efeito fotoelétrico, a luz se comportou como uma partícula. Experimentos de difração mais difíceis mostraram que os elétrons (assim como as outras partículas) também se comportavam como uma onda, mas só podemos detectar um número inteiro de elétrons (ou fótons).

A Mecânica Quântica incorpora uma dualidade onda-partícula e explica todos os fenômenos acima. Ao fazer isso, a Mecânica Quântica muda nossa compreensão da natureza de maneiras fundamentais. Enquanto as leis clássicas da física são determinísticas, QM é probabilística. Podemos apenas prever a probabilidade de uma partícula ser encontrada em alguma região do espaço.

Ondas eletromagnéticas como a luz são feitas de partículas que chamamos de fótons. Einstein, com base na fórmula de Plank, formulou a hipótese de que as partículas de luz tinham energia proporcional à sua frequência.

\ begin {displaymath} \ bgroup \ color {black} E = h \ nu \ egroup \ end {displaymath}

A nova ideia da Mecânica Quântica é que a probabilidade de cada partícula (em função da posição e do tempo) é igual ao quadrado de uma função de amplitude de probabilidade e que essas amplitudes de probabilidade obedecem a uma equação de onda. É muito parecido com o caso do eletromagnetismo, onde a densidade de energia é como o quadrado do campo e, portanto, a densidade de probabilidade do fóton é como o quadrado do campo, embora o campo seja feito de ondas. Portanto, as amplitudes de probabilidade são como os campos que conhecemos do eletromagnetismo em muitos aspectos.

DeBroglie assumiu \ bgroup \ color {black} $ E = h \ nu $ \ egrouppara fótons e outras partículas e usou a invariância de Lorentz (da relatividade especial) para derivar o comprimento de onda de partículas como elétrons.

\ begin {displaymath} \ bgroup \ color {black} \ lambda = {h \ over p} \ egroup \ end {displaymath}

O resto da mecânica ondulatória foi construída em torno dessas idéias, dando uma imagem completa que poderia explicar as medições acima e poderia ser testada com uma precisão muito alta, particularmente no átomo de hidrogênio. Passaremos vários capítulos explorando essas idéias.

Exemplo: Suponha que o fóton seja uma partícula com o comprimento de onda deBroglie padrão. Use a cinemática para derivar o comprimento de onda do fóton espalhado como uma função do ângulo para a dispersão de Compton. *

Gasiorowicz Capítulo 1

Capítulos Rohlf 3,6

Griffiths realmente não cobre isso.

Cohen-Tannoudji et al. Capítulo

Radiação de corpo negro *

Um corpo negro é aquele que absorve toda a radiação EM (luz ...) que o atinge. Para ficar em equilíbrio térmico, ele deve emitir radiação na mesma taxa em que a absorve, de modo que um corpo negro também irradie bem. (Fogões são pretos.)

A radiação de um objeto quente é familiar para nós. Objetos em torno da temperatura ambiente irradiam principalmente no infravermelho, como pode ser visto no gráfico abaixo.

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

Se aquecermos um objeto até cerca de 1.500 graus, começaremos a ver um brilho vermelho opaco e diremos que o objeto está incandescente. Se aquecermos algo até cerca de 5000 graus, perto da temperatura da superfície do Sol, ela irradia bem por todo o espectro visível e dizemos que é branco quente.

Ao considerar as placas em equilíbrio térmico, pode-se mostrar que a potência emissiva sobre o coeficiente de absorção deve ser a mesma em função do comprimento de onda, mesmo para placas de materiais diferentes.

\ begin {displaymath} \ bgroup \ color {black} {E_1 (\ lambda, T) \ over A_1 (\ lambda)} = {E_2 (\ lambda, T) \ over A_2 (\ lambda)} \ egroup \ end { displaymath}

 

Se houvesse diferenças, poderia haver um fluxo líquido de energia de uma placa para a outra, violando a condição de equilíbrio.

\ epsfig {file = figs / blackbody.eps, height = 4in}

Um corpo negro é aquele que absorve toda a radiação incidente sobre ele.

\ begin {displaymath} \ bgroup \ color {black} A_ {BB} = 1 \ egroup \ end {displaymath}

 

Assim, o poder Emissivo do corpo negro,, \ bgroup \ color {black} $ E (\ nu, T) $ \ egroupé universal e pode ser derivado dos primeiros princípios.

Um bom exemplo de corpo negro é uma cavidade com um pequeno orifício. Qualquer luz incidente sobre o orifício entra na cavidade e essencialmente nunca é refletida, uma vez que teria que sofrer um grande número de reflexos nas paredes da cavidade. Se tornarmos as paredes absorventes (talvez pintando-as de preto), a cavidade formará um corpo preto perfeito.

\ epsfig {file = figs / cavity.eps, height = 2in}

Há uma relação simples entre a densidade de energia em uma cavidade \ bgroup \ color {black} $ u (\ nu, T) $ \ egroupe o poder emissivo de corpo negro de um corpo negro que simplesmente vem de uma análise de quanta radiação, viajando na velocidade da luz, fluirá de um buraco no cavidade em um segundo.

\ begin {displaymath} \ bgroup \ color {black} E (\ nu, T) = {c \ over 4} u (\ nu, T) \ egroup \ end {displaymath}

 

A única parte que requer um pouco de reflexão é o 4 na equação acima.

Rayleigh e Jeans calculado t ele densidade de energia (em ondas EM) no interior de uma cavidade e, por conseguinte, o espectro de emissão de um corpo negro. Seu cálculo foi baseado na teoria EM simples e equipartição. Não apenas não concordou com os dados; dizia que toda a energia seria irradiada instantaneamente em radiação EM de alta frequência. Isso foi chamado de catástrofe ultravioleta .

\ begin {displaymath} \ bgroup \ color {black} u (\ nu, T) = {8 \ pi \ nu ^ 2 \ over c ^ 3} kT \ egroup \ end {displaymath}

 

\ epsfig {arquivo = figs / Rayleigh.eps, altura = 3in}

Plank encontrou uma fórmula que se ajusta bem aos dados em comprimentos de onda longos e curtos.

\ begin {displaymath} \ bgroup \ color {black} u (\ nu, T) = {8 \ pi \ nu ^ 2 \ over c ^ 3} {h \ nu \ over e ^ {h \ nu / kT} - 1} \ egroup \ end {displaymath}

 

Sua fórmula se ajustou aos dados tão bem que ele tentou encontrar uma maneira de derivá-los. Em poucos meses, ele foi capaz de fazer isso, postulando que a energia era emitida em quanta com \ bgroup \ color {black} $ E = h \ nu $ \ egroup. Embora haja um grande número de modos de cavidade em alta frequência, a probabilidade de emitir tais quanta de alta energia desaparece exponencialmente de acordo com a distribuição de Boltzmann. Assim, Plank suprimiu a radiação de alta frequência no cálculo e concordou com o experimento. Observe que a fórmula de corpo negro de Plank é a mesma no limite que \ bgroup \ color {black} $ h \ nu «kT $ \ egroup vai para zero em geral, \ bgroup \ color {black} $ \ nu $ \ egroupenquanto a fórmula de Rayleigh vai para o infinito.

\ epsfig {arquivo = figs / Plank.eps, altura = 3in}

É interessante notar que as ondas EM clássicas sugariam toda a energia térmica da matéria, tornando o universo um lugar muito frio para nós. A figura abaixo compara os dois cálculos com alguns dados em \ bgroup \ color {black} $ T = 1600 $ \ egroupgraus. (Também é surpreendente que o início da revolução quântica tenha vindo da radiação do corpo negro.)

\ epsfig {arquivo = figs / PlankBB1600.eps, altura = 3in}

Portanto, a potência emissiva por unidade de área é

\ begin {displaymath} \ bgroup \ color {black} E (\ nu, T) = {2 \ pi \ nu ^ 2 \ over c ^ 2} {h \ nu \ over e ^ {h \ nu / kT} - 1} \ egroup \ end {displaymath}

 

Podemos integrar isso sobre frequência para obter a potência total irradiada por unidade de área.

\ begin {displaymath} \ bgroup \ color {black} R (T) = {\ pi ^ 2 c \ over 60 (\ hbar c) ^ 3} k ^ 4 T ^ 4 = (5,67 \ vezes 10 ^ {- 8 } \ mathrm {W / m ^ 2 / \; ^ \ circ K ^ 4}) T ^ 4 \ egroup \ end {displaymath}

 

Exemplo: Qual é a temperatura na superfície solar? Use a intensidade da radiação na Terra e o pico do espectro em cerca de 500 nm para obter respostas. *

Exemplo: O fundo de micro-ondas cósmico é uma radiação de corpo negro com uma temperatura de 2,7 graus. Para qual frequência (e qual comprimento de onda) a intensidade atinge o pico? *

O efeito fotoelétrico

O efeito fotoelétrico mostra que a hipótese de Plank, usada para ajustar os dados do corpo negro, é realmente correta para a radiação EM. Einstein foi mais longe e propôs, em 1905, que a luz era composta de partículas com energia relacionada à frequência da luz \ bgroup \ color {black} $ E = h \ nu $ \ egroup. (Ele recebeu seu prêmio Nobel pelo efeito fotoelétrico, não por Relatividade Especial ou Geral.)

\ epsfig {arquivo = figs / Einstein.eps, altura = 4in}
\ epsfig {file = figs / photoelecd.eps, height = 2in}

Quando a luz atinge uma superfície polida (metal), os elétrons são ejetados.

\ epsfig {file = figs / photoelectric.eps, height = 4in}

As medições foram feitas da energia máxima do elétron versus frequência e intensidade da luz. A física clássica previu que a energia do elétron deveria aumentar com a intensidade, conforme o campo elétrico aumenta. Isso não é observado. A energia do elétron é independente da intensidade e depende linearmente da freqüência da luz, como pode ser visto na figura acima. A energia cinética dos elétrons é dada pela constante de Plank vezes a frequência da luz menos uma função de trabalho \ bgroup \ color {black} $ W $ \ egroup que depende do material.

\ begin {displaymath} \ bgroup \ color {black} {1 \ over 2} mv ^ 2 = h \ nu-W \ egroup \ end {displaymath}

 

Esta equação apenas expressa a conservação da energia \ bgroup \ color {black} $ h \ nu $ \ egroupsendo a energia do fóton e \ bgroup \ color {black} $ W $ \ egroup a energia de ligação dos elétrons no sólido. Os dados do efeito fotoelétrico apoiaram fortemente a hipótese de que a luz é composta de partículas (fótons).

O átomo de Rutherford *

A teoria clássica dos átomos afirmava que os elétrons estavam ligados a uma grande carga positiva do tamanho do átomo. Rutherford espalhou \ bgroup \ color {black} $ \ alpha $ \ egrouppartículas carregadas ( ) dos átomos para ver qual era a distribuição de carga positiva. Com uma distribuição de carga aproximadamente uniforme, suas \ bgroup \ color {black} $ \ alpha $ \ egrouppartículas de 5,5 MeV nunca deveriam ter retrodifundido porque tinham energia suficiente para superar a força coulomb de uma distribuição de carga, essencialmente arando direto pelo meio.

\ epsfig {file = figs / Rutherford_diag.eps, height = 4in}

Ele descobriu que as \ bgroup \ color {black} $ \ alpha $ \ egrouppartículas costumavam se espalhar em ângulos maiores que 90 graus. Seus dados podem ser explicados se o núcleo positivo de um átomo for muito pequeno. Para um núcleo muito pequeno, a força de Coulomb continua a aumentar à medida que se \ bgroup \ color {black} $ \ alpha $ \ egroup aproxima do núcleo, e o retroespalhamento é possível.

\ epsfig {file = figs / Rutherford.eps, height = 4in}

Exemplo: Use a eletrostática para estimar o quão pequeno um núcleo de ouro deve ser para retroespalhar uma partícula alfa de 5,5 MeV. *

Isso trouxe um novo problema. O tamanho atômico era conhecido por vários tipos de experimentos. Se os elétrons orbitam em torno do núcleo atômico , de acordo com as equações de Maxwell, eles devem irradiar energia à medida que se aceleram. Esta radiação não é observada e os estados básicos dos átomos são estáveis.

Na Mecânica Quântica, a localização do elétron ao redor do núcleo é limitada por causa da natureza ondulatória do elétron . Para o hidrogênio, onde não há problema de múltiplos corpos para tornar o cálculo desnecessariamente difícil, os níveis de energia podem ser calculados com muita precisão. O hidrogênio foi usado para testar a mecânica quântica à medida que se desenvolvia. Também usaremos bastante hidrogênio neste curso.

O espalhamento das \ bgroup \ color {black} $ \ alpha $ \ egrouppartículas de alta energia permitiu a Rutherford "ver" dentro do átomo e determinar que o núcleo atômico é muito pequeno. (Ele provavelmente destruiu todos os átomos que "viu".) A figura abaixo mostra a distribuição angular de Rutherford em seu experimento de espalhamento junto com vários usos subsequentes da mesma técnica, com partículas de energia cada vez mais altas. Vemos a descoberta de Rutherford do minúsculo núcleo , a descoberta da estrutura nuclear , a descoberta de um próton pontual dentro do núcleo , a descoberta da estrutura do próton , a descoberta de quarks dentro do próton e, finalmente, a falta de descoberta, até agora , de qualquerestrutura de quark .

\ epsfig {file = figs / Scattering.eps, height = 6in}

Para "ver'' essas coisas em distâncias cada vez menores, precisamos usar feixes de partículas com comprimento de onda cada vez menor e, portanto, com energia mais alta.

Espectro Atômico *

O hidrogênio foi, em última análise, o verdadeiro teste da teoria quântica. Medições de altíssima precisão foram feitas usando grades de difração. Estes foram bem compreendidos na QM não relativística e compreendidos ainda melhor na Teoria Quântica de Campos totalmente relativística, Eletrodinâmica Quântica.

A figura abaixo mostra os níveis de energia em hidrogênio , as transições entre os níveis de energia e o comprimento de onda da luz produzida na transição. A série Lyman cobre as transições para o estado fundamental e está além da parte visível do espectro. A série Balmer é devida a transições para o primeiro estado excitado e está no visível.

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

Na época de Plank \ bgroup \ color {black} $ E = h \ nu $ \ egroupmuitos dados existiam sobre as energias discretas nas quais os átomos irradiavam. Cada átomo tinha sua impressão digital de radiação única. Absorção em energias discretas também foi observada.

fórmula de Rydberg para as energias dos fótons emitidos pelo hidrogênio foi desenvolvida bem antes da explicação QM aparecer.

\ begin {displaymath} \ bgroup \ color {black} E = 13.6 \ mathrm {eV} \ left ({1 \ over n_1 ^ 2} - {1 \ over n_2 ^ 2} \ right) \ egroup \ end {displaymath}

Alguns dos estados nos átomos mais pesados ​​seguiram o mesmo tipo de fórmula. Melhores experimentos mostraram que as linhas espectrais eram freqüentemente divididas em um multipleto de linhas. Vamos entender essas divisões muito mais tarde no curso.

Átomos mais pesados ​​fornecem um espectro ainda mais rico, mas são muito mais difíceis de calcular. Muito boas técnicas de aproximação foram desenvolvidas. Com computadores e boa técnica, os níveis de energia de átomos mais complexos podem ser calculados. O espectro de mercúrio mostrado abaixo tem muito mais linhas do que as vistas em hidrogênio.

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

A figura abaixo mostra a parte visível do espectro para várias fontes atômicas ou moleculares. Para o Hidrogênio, vemos principalmente a série Balmer, com uma linha da Paschen-Bach. Os espectros de diferentes átomos são bastante distintos. As moléculas têm muitos outros tipos de excitações e podem produzir muitas frequências de luz.

\ epsfig {file = figs / color_spectra.eps, height = 4in}

O átomo de Bohr

Bohr postulou que os elétrons orbitavam o núcleo como planetas orbitando o sol. Ele conseguiu ajustar os dados para o Hidrogênio postulando que os elétrons orbitavam o núcleo em órbitas circulares, e que o momento angular é quantizado de forma que \ bgroup \ color {black} $ L = n \ hbar $ \ egroup, para \ bgroup \ color {black} $ n = 1,2,3 ... $ \ egroup. Isso é natural, pois \ bgroup \ color {black} $ \ hbar $ \ egrouptem unidades de momento angular. Bohr postou corretamente as energias do hidrogênio e o tamanho do átomo.

Equilíbrio de forças para órbitas circulares.

\ begin {displaymath} \ bgroup \ color {black} {mv ^ 2 \ over r} = {1 \ over 4 \ pi \ epsilon_0} {e ^ 2 \ over r ^ 2} \ egroup \ end {displaymath}

Suposição de quantização de momento angular.

\ begin {displaymath} \ bgroup \ color {black} L = mvr = n \ hbar \ egroup \ end {displaymath}

Resolva a velocidade.

\ begin {displaymath} \ bgroup \ color {black} v = {n \ hbar \ over mr} \ egroup \ end {displaymath}

Conecte a equação de força para obter a fórmula \ bgroup \ color {black} $ r $ \ egroup.

\ begin {eqnarray *} {mn ^ 2 \ hbar ^ 2r ^ 2 \ over m ^ 2r ^ 2r} & = & {e ^ 2 \ over 4 \ pi \ epsilon_0} ... ... {1 \ over r } & = & {me ^ 2 \ over 4 \ pi \ epsilon_0 \ hbar ^ 2} {1 \ over n ^ 2} \\ \ end {eqnarray *}

Agora queremos apenas conectar \ bgroup \ color {black} $ v $ \ egroup\ bgroup \ color {black} $ r $ \ egroupentrar na fórmula de energia. Escrevemos o potencial de hidrogênio em termos da constante de estrutura fina . \ bgroup \ color {black} $ \ alpha_ {SI} = {1 \ over 4 \ pi \ epsilon_0} {e ^ 2 \ over \ hbar c} \ approx {1 \ over 137} $ \ egroup

\ begin {eqnarray *} V (r) = {- \ alpha \ hbar c \ over r} \\ {1 \ over r} = {\ alpha mc \ over \ hbar} {1 \ over n ^ 2} \\ \ end {eqnarray *}

Agora calculamos os níveis de energia.

\ begin {eqnarray *} E & = & {1 \ over 2} mv ^ 2 + V (r) \\ E & = & {1 \ over 2} m \ left ({n \ hbar \ over ... ... 2 \ direita) - {\ alpha ^ 2 mc ^ 2 \ over n ^ 2} = - {\ alpha ^ 2 mc ^ 2 \ over 2n ^ 2} \\ \ end {eqnarray *}

O eV constante . A fórmula de Bohr fornece o espectro de energia do hidrogênio correto .

 \ bgroup \ color {black} $ {\ alpha ^ 2 mc ^ 2 \ over 2} = {511000 \ over 2 (137) ^ 2} = 13,6 $ \ egroup

Também podemos calcular o raio do estado fundamental da órbita de Bohr.

\ begin {eqnarray *} {1 \ over r} & = & {\ alpha mc ^ 2 \ over \ hbar c} \\ r & = & {(137) (1973) \ over 511000} = 0,53 \ AA \\ \ fim de {eqnarray *} 

Isso também é sobre o raio certo. O nome desse número, raio de Bohr , persiste até hoje, apesar do fato de que a suposição de Bohr está errada.

Embora o momento angular seja quantizado em unidades de \ bgroup \ color {black} $ \ hbar $ \ egroupo estado fundamental do Hidrogênio tem momento angular zero . Isso colocaria o elétron de Bohr no núcleo.

Bohr ajustou os dados, com algum elemento de verdade, mas seu modelo está ERRADO .

Derivações e Cálculos

Black Body Radiation Formulas *

Ainda não disponível

A constante de estrutura fina e o potencial de Coulomb

Vamos agora enfrentar pela primeira vez o problema de qual conjunto de unidades usar. Textos avançados normalmente usam unidades CGS em que a energia potencial é

\ begin {displaymath} \ bgroup \ color {black} V (r) = {- e ^ 2 \ over r} \ egroup \ end {displaymath}

enquanto as unidades internacionais padrão

\ begin {displaymath} \ bgroup \ color {black} V (r) = {1 \ over 4 \ pi \ epsilon_0} {- e ^ 2 \ over r} \ egroup \ end {displaymath}

Podemos contornar o problema definindo a constante de estrutura fina adimensional \ bgroup \ color {black} $ \ alpha $ \ egroup.

\ begin {eqnarray *} \ alpha_ {SI} & = & {1 \ over 4 \ pi \ epsilon_0} {e ^ 2 \ over \ hbar c} \ appro ... ... 37} \\ \ alpha_ {CGS } & = & {e ^ 2 \ over \ hbar c} \ approx {1 \ over 137} \\ \ end {eqnarray *}

Portanto, em qualquer conjunto de unidades, o potencial de hidrogênio é

\ begin {displaymath} \ bgroup \ color {black} V (r) = {- \ alpha \ hbar c \ over r} \ egroup \ end {displaymath}

Exemplos

A Temperatura Solar

Estime a temperatura solar usando

  • a intensidade da radiação solar na terra de 1,4 quilowatts por metro quadrado. ( m, m) $ r_ {sun} = 7 \ vezes 10 ^ 8 $ $ d_ {sol} = 1,5 \ vezes 10 ^ {11} $
  • e o espectro solar que atinge o pico a cerca de 500 nm.

Primeiro, calculamos a energia irradiada por unidade de área na superfície solar.

\ begin {displaymath} \ bgroup \ color {black} R = (1400 \ mathrm {W / m ^ 2}) (4 \ pi d_ {sol} ^ 2) / (4 \ pi r_ {sol} ^ 2) = 6,4 \ times 10 ^ 7 \ mathrm {W / m ^ 2} \ egroup \ end {displaymath}

 

Comparamos isso com a expectativa em função da temperatura.

\ begin {displaymath} \ bgroup \ color {black} R (T) = (5,67 \ times 10 ^ {- 8} \ mathrm {W / m ^ 2 / \; ^ \ circ K ^ 4}) T ^ 4 \ egroup \ end {displaymath}

 

e pegue

\ begin {eqnarray *} T ^ 4 & = & {6,4 \ times 10 ^ 7 \ over 5,67 \ times 10 ^ {- 8}} \\ T & = & 5800 \; ^ \ circ K \\ \ end {eqnarray *}

 

 

Vamos supor que os \ bgroup \ color {black} $ E (\ lambda, T) $ \ egrouppicos sejam de 500 nm, como mostra um dos gráficos. Precisamos nos transformar \ bgroup \ color {black} $ E (\ nu, T) $ \ egroup. Lembre-se de funções de distribuição. \ bgroup \ color {black} $ f (\ nu) d \ nu = g (\ lambda) d \ lambda $ \ egroup

\ begin {eqnarray *} E (\ nu, T) & = & {2 \ pi \ nu ^ 2 \ over c ^ 2} {h \ nu \ over e ^ {h \ nu / kT} -1} \\ ... ... ^ {h \ nu / kT} -1} = {2 \ pi \ nu ^ 4 \ over c ^ 3} {h \ nu \ over e ^ {h \ nu / kT} -1} \\ \ end {eqnarray *}

 

 

Este pico quando

\ begin {displaymath} \ bgroup \ color {black} {\ nu ^ 5 \ over e ^ {h \ nu / kT} -1} \ egroup \ end {displaymath}

 

é máximo.

\ begin {eqnarray *} {5 \ nu ^ 4 \ over e ^ {h \ nu / kT} -1} - {\ nu ^ 5 (h / kT) e ^ {h \ nu / kT} \ over (e ^ ... ...} - 1) \ over e ^ {h \ nu / kT}} = h \ nu / kT \\ 5 (1-e ^ {- h \ nu / kT}) = h \ nu / kT \\ \ end {eqnarray *}

Vamos definir \ bgroup \ color {black} $ y = h \ nu / kT $ \ egroupe resolver a equação

\ begin {displaymath} \ bgroup \ color {black} y = 5 (1-e ^ {- y}) \ egroup \ end {displaymath}

 

Resolvi isso iterativamente começando em y = 5 e comecei a \ bgroup \ color {black} $ y = 4,97 $ \ egroupsugerir

\ begin {displaymath} \ bgroup \ color {black} h \ nu = 4,97kT \ egroup \ end {displaymath}

Se tomarmos \ bgroup \ color {black} $ \ lambda = 500 $ \ egroupnm e \ bgroup \ color {black} $ \ nu = 6 \ vezes 10 ^ {14} $ \ egroup

\ begin {displaymath} \ bgroup \ color {black} T = {h \ nu \ over 5 k} = {(6,6 \ vezes 10 ^ {- 34}) (6 ... ... 10 ^ {14}) \ over (5) (1,4 \ vezes 10 ^ {- 23})} = 6 \ vezes 10 ^ 3 = 5700 \ egroup \ end {displaymath}

 

Isso concorda bem.

Radiação de corpo negro do universo inicial *

Encontre a frequência \ bgroup \ color {black} $ \ nu $ \ egroupem que o Emissivo \ bgroup \ color {black} $ E (\ nu, T) $ \ egroupé máximo para a radiação cósmica de fundo de 2,7 graus. Encontre o comprimento \ bgroup \ color {black} $ \ lambda $ \ egroupde onda \ bgroup \ color {black} $ E (\ lambda, T) $ \ egroupmáximo.

A radiação cósmica de fundo foi produzida quando o universo era muito mais quente do que agora. A maioria dos átomos do universo era ionizada e os fótons interagiam frequentemente com os íons ou elétrons livres. À medida que o universo esfriava, formando átomos neutros, os fótons se separavam da matéria e se propagavam pelo espaço. Hoje vemos esses fótons como a radiação de fundo. Como o universo está se expandindo, a radiação mudou para uma temperatura muito mais baixa do vermelho. Observamos cerca de 2,7 graus. A radiação de fundo é muito uniforme, mas estamos começando a observar não uniformidades no \ bgroup \ color {black} $ 10 ^ {- 5} $ \ egroupnível.

Do problema anterior, podemos dizer que o pico \ bgroup \ color {black} $ \ lambda $ \ egroupocorre quando

\ begin {eqnarray *} h \ nu & = & 5kT \\ \ nu & = & 5kT / h \\ \ lambda & = & ch / (5kT) = {(3 \ vezes 1 ... ... 6,6 \ vezes 10 ^ {- 34 }) \ over (5) (1,4 \ vezes 10 ^ {- 23}) (2,7)} = 1 mm \\ \ end {eqnarray *}

Da mesma forma, o pico \ bgroup \ color {black} $ \ nu $ \ egroupocorre quando

\ begin {displaymath} \ bgroup \ color {black} \ nu = 2,8 kT / h = {(1,4 \ vezes 10 ^ {- 23}) (2,7) \ over (6,6 \ vezes 10 ^ {- 34})} = 6 \ vezes 10 ^ {10} Hz \ egroup \ end {displaymath}

Dispersão Compton *

Compton espalhou fótons de alta energia (essencialmente) de elétrons livres em 1923. Ele mediu o comprimento de onda dos fótons espalhados como uma função do ângulo de espalhamento. A figura abaixo mostra o estado inicial (a) e o estado final, com o fóton espalhado por um ângulo \ bgroup \ color {black} $ \ theta $ \ egroupe o elétron recuando em um ângulo \ bgroup \ color {black} $ \ phi $ \ egroup. Os fótons eram provenientes do decaimento nuclear e, portanto, eram de energia alta o suficiente para que não importasse que os elétrons estivessem realmente ligados aos átomos. Queremos derivar a fórmula para o comprimento de onda do fóton espalhado em função do ângulo .

\ epsfig {arquivo = figs / Compton.eps, altura = 2,5 pol.}

Resolvemos o problema usando apenas a conservação de energia e momento. Vamos trabalhar em unidades em que \ bgroup \ color {black} $ c = 1 $ \ egroup por enquanto. Colocaremos as \ bgroup \ color {black} $ c $ \ egroupcostas no final. Suponha que o fóton esteja inicialmente se movendo na \ bgroup \ color {black} $ z $ \ egroupdireção com energia E e que ele se espalhe no \ bgroup \ color {black} $ yz $ \ egroupplano de forma que \ bgroup \ color {black} $ p_x = 0 $ \ egroupA conservação do momento dá

\ begin {displaymath} \ bgroup \ color {black} E = E '\ cos \ theta + p_e \ cos \ phi \ egroup \ end {displaymath}

e

\ begin {displaymath} \ bgroup \ color {black} E '\ sin \ theta = p_e \ sin \ phi. \ egroup \ end {displaymath}

Conservação de energia dá

\ begin {displaymath} \ bgroup \ color {black} E + m = E '+ \ sqrt {p_e ^ 2 + m ^ 2} \ egroup \ end {displaymath}

Nosso objetivo é resolver para \ bgroup \ color {black} $ E '$ \ egroupem termos de, \ bgroup \ color {black} $ \ cos \ theta $ \ egroupentão vamos nos certificar de que eliminamos o \ bgroup \ color {black} $ \ phi $ \ egroup. Continuando da equação de energia

\ begin {displaymath} \ bgroup \ color {black} E-E '+ m = \ sqrt {p_e ^ 2 + m ^ 2} \ egroup \ end {displaymath}

quadratura e cálculo \ bgroup \ color {black} $ p_e ^ 2 $ \ egroupdos componentes

\ begin {displaymath} \ bgroup \ color {black} E ^ 2 + E '^ 2 + m ^ 2-2EE' + 2mE-2mE '= (E-E' \ cos \ theta) ^ 2 + (E '\ sen \ theta) ^ 2 + m ^ 2 \ egroup \ end {displaymath}

 

e escrever os quadrados do lado direito

\ begin {displaymath} \ bgroup \ color {black} E ^ 2 + E '^ 2 + m ^ 2-2EE' + 2mE-2mE '= E ^ 2 + E' ^ 2 -2EE '\ cos \ theta + m ^ 2 \ egroup \ end {displaymath}

e removendo coisas que aparecem em ambos os lados

\ begin {displaymath} \ bgroup \ color {black} -2EE '+ 2mE-2mE' = - 2EE '\ cos \ theta \ egroup \ end {displaymath}

e agrupamento

\ begin {eqnarray *} m (E-E ') & = & EE' (1- \ cos \ theta) \\ {(E-E ') \ over EE'} & = & {(1- \ cos \. .. ... sobre m} \\ {1 \ sobre E '} - {1 \ sobre E} & = & {(1- \ cos \ theta) \ sobre m} \\ \ end {eqnarray *}

Já que em nossas unidades finas, \ bgroup \ color {black} $ \ lambda = h / p = h / E $ \ egroup

\ begin {displaymath} \ bgroup \ color {black} \ lambda '- \ lambda = {h \ over m} (1- \ cos \ theta). \ egroup \ end {displaymath}

Agora aplicamos a velocidade da luz para fazer com que as unidades tenham um comprimento.

\ begin {displaymath} \ bgroup \ color {black} \ lambda '- \ lambda = {hc \ over mc ^ 2} \ left (1- \ cos \ theta \ right) \ egroup \ end {displaymath}

Esses cálculos podem ser bastante frustrantes se você não decidir quais variáveis ​​deseja manter e quais deve eliminar de suas equações. Neste caso, eliminamos \ bgroup \ color {black} $ \ phi $ \ egroupusando a equação de energia e computação \ bgroup \ color {black} $ p_e ^ 2 $ \ egroup.

Tamanho nuclear de Rutherford

Se a carga positiva nos átomos de ouro fosse uniformemente distribuída em uma esfera ou raio de 5 Angstroms, qual é a \ bgroup \ color {black} $ \ alpha $ \ egroupenergia cinética máxima da partícula pela qual o \ bgroup \ color {black} $ \ alpha $ \ egrouppode ser espalhado de volta na direção de onde veio?

Para resolver isso, precisamos calcular o potencial no centro da distribuição de carga em relação ao potencial no infinito (que diremos que é zero). Isso nos diz diretamente a energia cinética em eV necessária para arar através da distribuição de carga.

O potencial na superfície do núcleo é onde Z é o número de prótons no átomo e R é o raio nuclear. Essa é a parte fácil. Agora precisamos integrar nosso caminho para o centro. \ bgroup \ color {black} $ {1 \ over 4 \ pi \ epsilon_0} {Ze \ over R} $ \ egroup

\ begin {displaymath} \ bgroup \ color {black} V = {1 \ over 4 \ pi \ epsilon_0} {Ze \ over R} - \ in ... ... R ^ 0 {1 \ over 4 \ pi \ epsilon_0} {r ^ 3 \ over R ^ 3} {Ze \ over r ^ 2} dr \ egroup \ end {displaymath}

O dá a fração da carga nuclear dentro de um raio . \ bgroup \ color {black} $ {r ^ 3 \ over R ^ 3} $ \ egroup\ bgroup \ color {black} $ r $ \ egroup

\ begin {eqnarray *} V = {1 \ over 4 \ pi \ epsilon_0} {Ze \ over R} - {1 \ over 4 \ pi \ epsilon_0 ... ... t ({Ze \ over R} + { Ze \ over 2R} \ right) = {3Ze \ over 8 \ pi \ epsilon_0R} \\ \ end {eqnarray *}

Então

\ begin {displaymath} \ bgroup \ color {black} V = {(3) (79) (1,6 \ times 10 ^ {- 19} \ mathrm {C}) \ ... ... {C ^ 2 / Nm ^ 2}) R} = {1,7 \ vezes 10 ^ {- 7} \ sobre R} \ mathrm {Nm / C} \ egroup \ end {displaymath}

Então, é a energia cinética em eV necessária para que uma partícula de carga \ bgroup \ color {black} $ + e $ \ egroupatravesse o centro de uma distribuição esférica de carga. A \ bgroup \ color {black} $ \ alpha $ \ egrouppartícula realmente tem carga, \ bgroup \ color {black} $ + 2e $ \ egroupentão precisamos multiplicar por 2. Para um raio nuclear de 5 Å ou metros, precisamos de cerca de 680 eV para arar através do núcleo. Para o raio nuclear real de cerca de 5 Fermis ou metros, precisamos de 68 MeV para arar. \ bgroup \ color {black} $ 5 \ times10 ^ {- 10} $ \ egroup \ bgroup \ color {black} $ 5 \ times10 ^ {- 15} $ \ egroup

Vamos comparar os números das unidades SI acima com meu método sugerido de usar a constante de estrutura fina ... Colocando a carga alfa de \ bgroup \ color {black} $ 2e $ \ egroup.

\ begin {displaymath} \ bgroup \ color {black} U = {6Ze ^ 2 \ over 8 \ pi \ epsilon_0R} = {3Z \ alpha \ ... ...) (79) (1973) \ over (137) (5)} \ mathrm {eV} = 683 \ mathrm {eV} \ egroup \ end {displaymath}

Isso é mais fácil. 

Problemas de teste de amostra

  1. Qual é o comprimento de onda máximo da radiação eletromagnética que pode ejetar elétrons de um metal com função de trabalho de 3 eV? (3 pontos)
    Resposta
    A energia mínima do fóton necessária para eliminar um elétron é 3 eV. Só precisamos converter isso para comprimento de onda.
    \ begin {eqnarray *} E & = & h \ nu = 3 eV \\ \ nu & = & {c \ over \ lambda} \\ {2 \ pi \ hbar c \ o ... ... over 3 eV} & = & {2 \ pi (1973 eV \ AA) \ acima de 3 eV} \ aprox 4000 \ AA \\ \ end {eqnarray *}
  2. Com base no eletromagnetismo clássico e na mecânica estatística, Rayleigh calculou a densidade de energia dentro de uma cavidade. Ele descobriu que, a uma temperatura T, a densidade de energia em função da frequência era
    \ begin {displaymath} u (\ nu, T) = {8 \ pi \ nu ^ 2 \ over c ^ 3} k_B T. \ end {displaymath}

    Por que isso está relacionado à radiação do corpo negro? Por que isso estava em desacordo óbvio com a observação?
  3. Qual é o comprimento de onda máximo da radiação eletromagnética que pode ejetar elétrons de um metal com função de trabalho de 2 eV?
  4. Diga simplesmente que problema com a radiação de corpo negro fez com que Plank propusesse a relação$ E = h \ nu $para a luz.
  5. A função de trabalho de um metal é 2 eV. Suponha que um feixe de luz de comprimento de onda $ \ lambda $incida sobre uma superfície polida do metal. Trace a energia máxima do elétron (em eV) dos elétrons ejetados do metal $ \ lambda $em relação aos Angstroms. Certifique-se de rotular ambos os eixos com alguns valores numéricos.

 

 

Para ondas EM , a intensidade e, portanto, a probabilidade de encontrar um fóton, é proporcional ao quadrado dos campos. Os campos obedecem à equação de onda. Os campos de duas fendas podem adicionar padrões de interferência de forma construtiva ou destrutiva. Os campos \ bgroup \ color {black} $ E $ \ egroup\ bgroup \ color {black} $ B $ \ egroupestão 90 graus defasados ​​e ambos contribuem para a intensidade.

Usaremos as mesmas idéias para elétrons , embora os detalhes do campo variem um pouco porque elétrons e fótons são tipos de partículas um tanto diferentes. Para ambas as partículas, o comprimento de onda é dado por 

\ begin {displaymath} \ bgroup \ color {black} \ lambda = {h \ over p} \ egroup \ end {displaymath}

e a frequência por

\ begin {displaymath} \ bgroup \ color {black} E = h \ nu = \ hbar \ omega. \ egroup \ end {displaymath}

Usaremos uma amplitude de probabilidade complexa \ bgroup \ color {black} $ \ psi (x, t) $ \ egroup para o elétron. As partes reais e imaginárias estão defasadas como os campos EM. A onda progressiva com momentum \ bgroup \ color {black} $ p $ \ egroupe energia, \ bgroup \ color {black} $ E $ \ egroupentão, é

 

\ begin {displaymath} \ bgroup \ color {black} \ psi (x, t) \ propto e ^ {i (kx- \ omega t)} = e ^ {i (px-Et) / \ hbar} \ egroup \ fim {displaymath}

probabilidade de encontrar um elétron é igual ao quadrado absoluto da amplitude de probabilidade complexa.

 

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

 

 

(Superaremos o problema de que essa probabilidade é 1 em todos os lugares para nossa função de onda simples.)

Acabamos de incluir a maior parte da física da Mecânica Quântica. Muito do que fizermos no restante do curso será deduzido do parágrafo acima. Nossa contribuição veio de deBroglie e Plank, com suporte de experimentos.

Vamos resumir a entrada da física novamente.

  • Partículas livres são representadas por funções de onda complexas com uma relação entre suas propriedades de partícula - energia e momento, e suas propriedades de onda - frequência e comprimento de onda dados por Plank e deBroglie.
  • O quadrado absoluto da função de onda fornece a função de distribuição de probabilidade. A mecânica quântica apenas nos diz a probabilidade.
  • Podemos fazer sobreposições de nossas funções de onda de partículas livres para criar estados que não têm momento definido. Veremos que qualquer estado pode ser feito a partir da superposição de estados de partículas livres com diferentes momentos.

Agora temos uma dualidade onda-partícula para todas as partículas; no entanto, a física agora apenas nos diz a probabilidade de ocorrerem alguns eventos quânticos. Perdemos o poder preditivo completo da física clássica.

Gasiorowicz Capítulo 1

Rohlf Capítulo 5

Griffiths 1.2, 1.3

Cohen-Tannoudji et al. Capítulo

 

Derivações e cálculos

Revisão de números complexos

Esta é uma revisão simples, mas você deve certificar-se de usar números complexos corretamente. Um dos erros mais comuns em problemas de teste é esquecer de usar o conjugado complexo ao calcular uma probabilidade.

Um número complexo \ bgroup \ color {black} $ c = a + ib $ \ egroupconsiste em uma parte real \ bgroup \ color {black} $ a $ \ egroupe uma parte imaginária \ bgroup \ color {black} $ ib $ \ egroup. (Escolhemos \ bgroup \ color {black} $ a $ \ egroup\ bgroup \ color {black} $ b $ \ egrouppara ser números reais.) \ bgroup \ color {black} $ i $ \ egroupÉ a raiz quadrada de -1.

O conjugado complexo de \ bgroup \ color {black} $ c $ \ egroupé \ bgroup \ color {black} $ c ^ * = a-ib $ \ egroup. (Basta mudar o sinal de todos os \ bgroup \ color {black} $ i $ \ egroup.)

O quadrado absoluto de um número complexo é calculado multiplicando-o por seu conjugado complexo.

 

\ begin {displaymath} \ bgroup \ color {black} \ vert c \ vert ^ 2 = c ^ * c = (a-ib) (a + ib) = a ^ 2 + iab-iab + b ^ 2 = a ^ 2 + b ^ 2 \ egroup \ end {displaymath}

Isso dá a magnitude ao quadrado do número complexo. O quadrado absoluto é sempre real.

Usaremos exponenciais complexas o tempo todo. 

\ begin {eqnarray *} e ^ {i \ theta} & = & \ cos \ theta + i \ sin \ theta \\ e ^ {- i \ theta} & = & = & \ cos \ theta-i \ sin \ theta \ \ \ end {eqnarray *}

Você pode verificar se o quadrado absoluto dessas exponenciais é sempre 1. Elas costumam ser chamadas de fator de fase.

Podemos escrever e . \ bgroup \ color {black} $ \ sin \ theta = {e ^ {i \ theta} -e ^ {- i \ theta} \ over 2i} $ \ egroup \ bgroup \ color {black} $ \ cos \ theta = {e ^ {i \ theta} + e ^ {- i \ theta} \ over 2} $ \ egroup

Como acontece com outros exponenciais, podemos multiplicá-los adicionando os expoentes.

 

\ begin {displaymath} \ bgroup \ color {black} e ^ {ikx} e ^ {- i \ omega t} = e ^ {i (kx- \ omega t)} \ egroup \ end {displaymath}

Revisão das Ondas Viajantes

Uma onda viajante normal pode ser dada por

 

\ begin {displaymath} \ bgroup \ color {black} \ cos (kx- \ omega t). \ egroup \ end {displaymath}

 

A fase da onda passa \ bgroup \ color {black} $ 2 \ pi $ \ egroupem um comprimento de onda em \ bgroup \ color {black} $ x $ \ egroup. Portanto, o comprimento de onda \ bgroup \ color {black} $ \ lambda $ \ egroupsatisfaz

 

\ begin {displaymath} \ bgroup \ color {black} k \ lambda = 2 \ pi. \ egroup \ end {displaymath}

 

Da mesma forma, a fase passa \ bgroup \ color {black} $ 2 \ pi $ \ egroupem um período \ bgroup \ color {black} $ \ tau $ \ egroupde tempo.

 

\ begin {displaymath} \ bgroup \ color {black} \ omega \ tau = 2 \ pi \ egroup \ end {displaymath}

 

\ bgroup \ color {black} $ \ omega $ \ egroupé a frequência angular. Ele muda a \ bgroup \ color {black} $ 2 \ pi $ \ egroupcada ciclo. A frequência \ bgroup \ color {black} $ \ nu $ \ egroupaumenta em 1 a cada ciclo, então

 

\ begin {displaymath} \ bgroup \ color {black} \ omega = 2 \ pi \ nu. \ egroup \ end {displaymath}

 

Não há razão para memorizar essas equações. Eles devem ser óbvios.

Vamos ver o quão rápido um dos picos da onda se move. Isso é chamado de velocidade de fase. No momento \ bgroup \ color {black} $ t = 0 $ \ egroup, há um pico em \ bgroup \ color {black} $ x = 0 $ \ egroup. Este é o pico para o qual o argumento do cosseno é 0. No momento \ bgroup \ color {black} $ t = 1 $ \ egroup, o argumento é zero quando \ bgroup \ color {black} $ kx = \ omega t $ \ egroupou em . Se calcularmos a velocidade de fase tomando , obtemos \ bgroup \ color {black} $ x = {\ omega \ over k} $ \ egroup \ bgroup \ color {black} $ {\ Delta x \ over \ Delta t} $ \ egroup

 

\ begin {displaymath} \ bgroup \ color {black} v _ {\ mathrm {phase}} = {\ omega \ over k}. \ egroup \ end {displaymath}

 

Ou seja, um dos picos dessa onda viaja com velocidade de \ bgroup \ color {black} $ \ omega \ over k $ \ egroup.

 

\ begin {displaymath} \ bgroup \ color {black} v = {\ omega \ over k} = {2 \ pi \ nu \ over {2 \ pi \ over \ lambda}} = \ nu \ lambda \ egroup \ end { displaymath}

 

Em QM não-relativista, temos \ bgroup \ color {black} $ \ hbar k = p $ \ egroup\ bgroup \ color {black} $ E = \ hbar \ omega $ \ egroupe , por isso, \ bgroup \ color {black} $ E = {p ^ 2 \ over 2m} $ \ egroup

 

\ begin {displaymath} \ bgroup \ color {black} \ omega (k) = {E \ over \ hbar} = {\ hbar ^ 2k ^ 2 \ over 2m \ hbar} = {\ hbar k ^ 2 \ over 2m} \ egroup \ end {displaymath}

 

Você deve se lembrar que um pulso se moverá na velocidade do grupo que é dada por

 

\ begin {displaymath} \ bgroup \ color {black} v_g = \ left ({d \ omega \ over dk} \ right) = {2 \ hbar k \ over 2m} = {\ hbar k \ over m} = {p \ over m}. \ egroup \ end {displaymath}

 

(A velocidade de fase para o caso não relativístico é .) \ bgroup \ color {black} $ v_p = {p \ over 2m} $ \ egroup

Problemas de teste de amostra

  1. Escreva as duas funções de onda de partícula livre (não normalizadas) para uma partícula de energia cinética E. Inclua a dependência de tempo adequada e expressões para outras constantes em termos de E.