Resumo
Nesse trabalho introduzimos os operadores de Schrödinger de forma ingênua e despretensiosa, mas munida de certa precisão e didática. São mencionados os tipos de potenciais mais estudados, os periódicos e os aleatórios e, via simulação numérica, comparados os autovalores da matriz truncada de um operador munido de um potencial periódico e com este mesmo potencial perturbado por uma variável aleatória contínua com distribuição uniforme. Avaliamos qual a magnitude necessária da variável aleatória para produzir diferença significativa entre os conjuntos de autovalores.
Palavras-chave:
Operadores de Schrödinger; potencial periódico; potencial perturbado
Abstract
In this paper we introduce Schrödinger's operators in a naive and unpretentious way, but with some precision and didactics. The most studied periodic and random potential types are mentioned and, by numerical simulation, compare the eigenvalues of an operator's truncated matrix with a periodic potential and the same potential disturbed by an evenly distributed continuous random variable. We evaluated the magnitude of the random variable required to produce a significant difference between the eigenvalue sets.
Keywords:
Schrödinger operators; periodic potential; disturbed potential
1. Introdução
A física de matéria condensada de um elétron é uma simplificação que tem tido um sucesso notável na descrição de certas propriedades dos cristais. As noções básicas de materiais isolantes, condutores e semi-condutores podem ser explicadas pelo espectro de energia de um único elétron que se move sob a influência de um arranjo periódico de átomos e do princípio de exclusão de Pauli [1][1] M. Massimi, Pauli's Exclusion Principle: The Origin and Validation of a Scientific Principle (Cambridge University Press, Cambridge, 2005).. O último determina a ocupação dos estados por elétrons não interagentes. No modelo de um elétron, a interação entre os elétrons é negligenciada e as propriedades espectrais e de transporte do material são descritas pelo modelo de uma partícula do operador Hamiltoniano de Schrödinger.
De fato, em condutores e semi-condutores, elétrons livres têm que partilhar o espaço total disponível no interior do material - e por isso seus níveis de energia se empilham criando uma estrutura de bandas a partir de cada nível de energia atômico. Em bons condutores (metais) os elétrons estão tão fortemente degenerados que eles acabam por não contribuir de forma significativa para a capacidade térmica do metal. Muitas propriedades mecânicas, elétricas, magnéticas, ópticas e químicas dos sólidos são consequências diretas da repulsão de Pauli entre elétrons livres ou semi-livres. Para mais detalhes, veja [1][1] M. Massimi, Pauli's Exclusion Principle: The Origin and Validation of a Scientific Principle (Cambridge University Press, Cambridge, 2005). e [2][2] D.J. Griffiths, Introduction to Quantum Mechanics (Prentice Hall, New Jersey, 2004), 2ª ed..
Descrevendo a evolução temporal do elétron, a equação de Erwin Schrödinger dependente do tempo é dada por
onde t é a variável temporal, x é a variável espacial e é a constante de Planck dividida por . O operador Hamiltoniano H será, por simplicidade tomando um sistema de unidades adequado para o qual e , escrito como
Ainda, é possível reescrever a equação (1) usando o operador H como sendo para o qual o problema de existência e unicidade das soluções da equação (1) com o estado inicial dado pode ser reduzido ao problema do estabelecimento da auto-adjunticidade para o operador linear H, ou seja, se o operador for hermitiano está garantida a existência e unicidade de soluções com estado inicial .
A equação de Schrödinger (1) tem raízes no modelo ondulatório da mecânica quântica, começando da ideia que uma partícula livre se comportaria como uma onda plana e poderia então ser descrita por uma função de onda da forma
para todo e para todo , a qual se propaga com velocidade constante na direção . O A da equação é um número complexo não nulo que representa a constante de normalização. Usando algumas das leis fundamentais da teoria quântica, a saber, a equação de Einstein
e a relação de Broglie para o comprimento de onda
onde
chega-se a equação (1) quando se quer determinar a evolução temporal de começando de algum estado inicial . Com efeito, a energia mecânica total E de uma partícula de massa m em termos da energia potencial , sendo p o momentum da partícula, é
Usando a equação (6), temos
ou seja, , e pela equação (8), a equação (7) fica
Nos termos da equação (9), é particularmente difícil de escrever uma equação de onda para uma partícula por conta da energia potencial, por isso a ignoraremos por enquanto. Se a partícula tem massa nula, por exemplo um fóton, então a equação de onda eletromagnética seria apropriada, isto é, a equação de movimento é
onde é a amplitude dependente do tempo e do espaço (campo elétrico ou magnético) da onda, e c é a velocidade da luz. A solução da equação (10) é a onda plana
pois o movimento de onda sempre é descrito como função de , onde é a velocidade da onda, e . Para ondas eletromagnéticas, , a velocidade da luz.
Aplicaremos a equação (10) ao caso de uma partícula material (massa não nula), mas, por simplicidade, continuemos com , caso denominado de "partícula livre". Substituindo a amplitude pela função de onda , considerando e assumindo uma solução tipo onda plana análoga a equação (11), inserimos esta solução na equação (10), e após dividi-la por de ambos os lados da equação, obtemos
o que é inconsistente com a equação (9) pois, com , e não . Para obter apenas a primeira potência de , devemos derivar apenas uma vez com respeito ao tempo. Isto resulta na modificação da equação (10), alterando por e por , obtendo
onde K é uma constante. Inserindo a equação (11) na equação (13) e isolando K,
que nos dá
Mas, da equação (9), com , temos que implica em , assim
e portanto
Multiplicando ambos os lados por temos
Esta equação da onda é aplicável somente a uma partícula livre, isto é, uma partícula para a qual o comprimento de onda de de Broglie é constante por toda a parte. O comprimento de onda de de Broglie é constante pois a energia total é presumidamente constante, então a energia cinética, e portanto o momentum, são constantes. No caso de uma energia potencial constante, ainda é simples, pois o comprimento de onda de de Broglie é constante, e é um caso particular de constante. Portanto podemos, sem restrições, escrever
que é consistente com a equação (9) com no lugar de . Em mecânica quântica, quantidades observáveis como a energia, momentum e posição são representadas por operadores. Não é um grande salto trocar em (21) por , para então obter a forma unidimensional da equação de Schrödinger dependente do tempo apresentada em (1),
Para mais detalhes, veja [3][3] C.E. Burkhardt e J.J. Leventhal, Foundations of Quantum Physics (Springer, New York, 2008)..
Por simplicidade vamos considerar o exemplo em que dado por , numa caixa rígida satisfazendo se ou . Suponhamos também . Por continuidade, temos , logo e , de onde , ou seja, de onde , com e
Escrevendo como uma série de Fourier, temos
com
Como , temos
Como é uma função senoidal, sua segunda derivada também é uma função senoidal e satisfaz
de onde
Fazendo em (7), por (4) e (9), temos
logo
e como supomos , na verdade temos , onde K é a energia cinética, de onde, usando (30), a equação (28) se torna
Para um potencial não nulo, temos que , logo
nos dá
chamada de equação de Schrödinger independente do tempo unidimensional que é obtida resolvendo a equação (1) por meio da separação de variáveis, como em [4][4] A.M. Hinz, em: Spectral Theory of Schrödinger Operators, editado por R. del Rio e C. Villegas-Blas (American Mathematical Society, Rhode Island, 2004)..
Para o caso de maior dimensão, digamos d, o operador diferencial de segunda ordem dá lugar ao operador Laplaciano
e a equação de Schrödinger independente do tempo é escrita como
Escrevendo o operador Hamiltoniano como , então a equação (35) é de fato uma equação de autovalores (e autovetores) para H, que pode ser escrita de forma resumida simplesmente como
1.1. Resolvendo a equação de Schrödinger
Utilizando o método de separação de variáveis, podemos escrever
para e . Então, para com , temos
de onde e . Além disso, segue
e
Avaliando em (2), temos
e para , tal que , f satisfaz a equação diferencial ordinária
onde
cujas soluções gerais são do tipo , para . Estas soluções gerais caracterizam completamente a função com relação ao tempo.
Como , ao substituirmos obtemos
de onde vem
para todo , ou seja, a função u satisfaz a equação de Schrödinger independente do tempo.
Seja . Dizemos que f é limitada se existe algum tal que , para todo . Então, a evolução temporal dada por f é limitada se, e só se, o autovalor é real, uma vez que
e só é limitado se . Essa é uma situação fisicamente relevante. Para um maior aprofundamento, consulte [4][4] A.M. Hinz, em: Spectral Theory of Schrödinger Operators, editado por R. del Rio e C. Villegas-Blas (American Mathematical Society, Rhode Island, 2004).. Vale mencionar que os autovalores são de fato reais, uma vez que o operador é auto-adjunto. Com efeito, considerando , onde H é um operador Hamiltoniano autoadjunto (Hermitiano), logo , temos
Por definição, temos
então
assim , e portanto . O número E é o autovalor de H associado a autofunção .
2. Operadores de Schrödinger
Discretos
Os operadores de Schrödinger discretos são definidos num espaço de sequências, a saber, o espaço de Hilbert
cuja norma é denotada por
Quando uma partícula se move no reticulado , supondo que , a probabilidade de encontrarmos tal partícula no ponto é dada por , onde é uma sequência indexada por . É conveniente equipar com a norma
Esta norma reflete a estrutura de reticulado de . Dois vértices n e m de são conectados por uma aresta, se eles forem vizinhos próximos, ou seja, se . Para a norma dá o comprimento do menor caminho, não necessariamente único, entre n e m. O operador , a energia cinética, é o análogo discreto de , a saber
Este operador é chamado de Laplaciano discreto para o reticulado . Sua representação matricial é dada por
ou seja,
para o caso unidimensional. Por simplicidade, consideramos , assim o operador de Schrödinger discreto unidimensional é dado por na matriz (62):
A equação de Schrödinger independente do tempo, , toma a forma matricial (63):
O operador é representado pela equação de recorrência, para e , dada por
cuja equação de autovalores e autoestados é dada por
3. Potenciais
Considere o potencial periódico em uma dimensão satisfazendo , para , onde a é a constante da rede discreta e n é um inteiro. A equação de Schrödinger para autovalores de energia para um elétron é dada por
Claramente, a função de onda também satisfaz a mesma equação. Assim, se diferencia de somente por uma fase independente de x, ou seja,
onde k é um número real. A equação (67) representa, em uma de suas formas, o Teorema de Bloch[5][5] F. Bloch, Zeitschrift für Physik 52, 555 (1929)., a qual gera importantes consequências na teoria de condução de elétrons. Considerando que e satisfazem a mesma equação, elas são linearmente dependentes, ou seja,
Usando o mesmo argumento outra vez, obtemos
Como a função de onda não se anula em geral, obtemos , ou seja,
conhecida como equação funcional exponencial de Cauchy.
Lema 3.1A equação funcional exponencial de Cauchytem por solução geral, ondeé uma função aditiva, ou.
Demonstração. De fato, é solução da equação (69), para todo . Suponha que para algum . Logo,
para todo , então se a função c se anula num ponto, ela é identicamente nula. Note que satisfaz a equação (69), logo esta admite solução não trivial. Além disso, a função é estritamente positiva. Com efeito, tome por exemplo na equação (69), logo teremos
Tomando o logaritmo natural em ambos os membros de (69), obtemos , e chamando , verificamos que
ou seja, a função A é uma função aditiva. Além disso, é uma solução geral para a equação (69), o que demonstra o resultado. ▪
Em particular, é uma solução da equação (69), onde é uma constante. Como a função na equação (67) deve ser finita para todos os intervalos, a constante deve ser um número imaginário puro , onde k é um número real, logo
Combinando as equações (68) com (73), obtemos
o que mostra o Teorema 3.2 seguinte:
Teorema 3.2 (de Bloch) As soluções da equação de Schrödinguer
para um potencial periódico satisfazem
onde r é o período da função de onda.
Para mais detalhes, veja [6][6] S. Fujita e K. Ito, Quantum Theory of Conducting Matter (Springer, New York, 2007). e [7][7] C. Kittel, Introduction to Solid State Physics (John Wiley and Sons, New Jersey, 1996), 7ª ed..
Nos cristais, os átomos ou os núcleos são distribuídos em uma rede periódica (digamos a rede por simplicidade) de uma maneira completamente regular. Vamos assumir que uma partícula (elétron) no ponto é estimulada por um potencial da forma devido a um átomo (ou íon, ou núcleo) localizado no ponto . Aqui, a constante q (carga ou constante de acoplamento, em termos físicos) poderia ser absorvida na função f. No entanto, uma vez que vamos variar essa quantidade de átomo para átomo mais tarde, é útil escrever o potencial da forma acima. Então, em um cristal regular, nossa partícula está exposta a um potencial total
A função f é chamada de potencial de um único sítio para distinguir do potencial total V. O potencial V em (77) é periódico com respeito ao reticulado , ou seja, para todo e . A teoria matemática dos operadores de Schrödinger com potenciais periódicos está bem desenvolvida e baseia-se em uma análise minuciosa das propriedades de simetria dos operadores periódicos. Por exemplo, sabe-se que tais operadores apresentam, em geral, um espectro com estrutura de banda, isto é, para ,
Este espectro também é conhecido como sendo espectro absolutamente contínuo.%, todavia, a maioria dos sólidos não constitui um cristal ideal.
Nas últimas décadas, os operadores de Schrödinger aleatórios, que se originaram na física da matéria condensada, foram estudados intensamente e de forma muito produtiva. A teoria está na interface de uma série de campos matemáticos: a teoria dos operadores, as equações diferenciais parciais, a teoria das probabilidades, em particular o estudo de processos estocásticos, o de passeios aleatórios e o movimento browniano em um ambiente aleatório. Os operadores de Schrödinger aleatórios modelam sólidos desordenados. Os sólidos aparecem na natureza de várias formas, sendo que em alguns casos são (quase) totalmente ordenados.
A maioria dos sólidos não constitui um cristal ideal. As posições dos átomos podem se desviar das posições de rede ideais de maneira não regular devido a imperfeições no processo de cristalização, ou as posições dos átomos podem estar completamente desordenadas, como é o caso em materiais amorfos ou vítreos. O sólido também pode ser uma mistura de vários materiais, como é o caso, por exemplo, de ligas ou semicondutores dopados (com impurezas intencionalmente introduzidas com propósito de modificar suas propriedades elétricas). Em todos esses casos, parece razoável considerar o potencial como uma quantidade aleatória. Por exemplo, se o material for puro, mas as posições dos átomos se desviam das posições da rede ideais aleatoriamente, podemos considerar um potencial aleatório da forma
Aqui, são variáveis aleatórias que descrevem o desvio do i-ésimo átomo da posição i da rede. Pode-se, por exemplo, assumir que as variáveis aleatórias são independentes e identicamente distribuídas. O sub-índice do potencial V aparece para deixar claro que depende dos parâmetros aleatórios (desconhecidos). Para modelar um material amorfo como vidro ou borracha, assumimos que os átomos do material estão localizados em pontos completamente aleatórios no espaço. Esse potencial aleatório pode ser formalmente escrito como
Para mais detalhes, veja [8][8] W. Kirsch, in: Random Schrödinger operators (Société Mathématique de France, Paris, 2008), v. 25, p. 1..
4. Metodologia
Realizamos simulações no software Mathematica, para um caso bastante simplificado, encontrando o espectro (conjunto de autovalores, na verdade) do operador periódico cujo potencial foi dado pela tal função. Acoplamos ao potencial uma variável aleatória a fim de tornar o potencial periódico aleatório e executamos simulações computacionais para comparar estatisticamente os resultados do espectro do operador de Schrödinger com cada um dos potenciais. Primeiramente simulamos usando uma matriz truncada do operador com potencial periódico, que usamos para comparar com todos os outros. Após encontrarmos diferenças significativas entre o espectro do operador periódico e do perturbado, realizamos 20 simulações similares para efeito de amostragem.
5. Resultados
Conforme [8][8] W. Kirsch, in: Random Schrödinger operators (Société Mathématique de France, Paris, 2008), v. 25, p. 1., em sistemas unidimensinais, um distúrbio arbitrariamente pequeno em um operador com potencial periódico, tornando-o levemente desordenado, mudará o espectro total de absolutamente contínuo (composto de bandas) para puramente pontual e, portanto, de condutor para isolante.
Vamos considerar um "caminho" unidimensional na rede bidimensional. Para investigar nosso caso, consideramos o potencial
onde o denominador , que é o dobro da variância, foi calculado de modo a ajustar o período da função, , e o fator 10 que multiplica a série serve simplesmente para facilitar a visualização dos autovalores na simulação. De fato, a função V está bem definida para todo , um vez que
e a integral da distribuição Gaussiana converge. O 40 compensa o valor no qual a série é superior a integral numa vizinhança de .
Para fins de simulação, truncamos a matriz (bi-infinita) (Equação (62)) do operador em uma matriz de com 201 linhas e 201 colunas, sendo que o valor 201 foi escolhido pois tomamos e . Assim, o potencial fica dado por
e o período da função foi dividido em 20 partes, de modo a facilitar a visualização. O primeiro bloco da matriz, com 20 linhas e 20 colunas (com valores arredondados) é dado pela matriz (83):
Na Figura 1, temos o conjunto dos autovalores no plano complexo. Embora com 201 valores, o espectro está longe de ser normalmente distribuído, uma vez que o teste de normalidade de Shapiro-Wilk nos deu um p-valor de , o que é bem menor que os 0,05 costumeiros, para um intervalo de 95% de confiança. Para avaliarmos o espectro com "impurezas" induzidas pela aleatoriedade, considere as variáveis aleatórias e com e distribuição uniforme, e o vetor aleatório . Uma vez que estamos lidando com o truncamento da matriz do operador, não poderemos perceber diferenças significativas simplesmente comparando o retrato de cada espectro. Com propósito de comparar se há diferenças significativas entre os dois espectros, utilizamos um teste de hipóteses. Para o potencial com perturbação ,
O espectro do operador com potencial perturbado (com fator ) no plano complexo é apresentado na Figura 2. O espectro não é normalmente distribuído, com p-valor de no teste de Shapiro-Wilk.
Utilizando o teste de hipóteses para dados não-paramétricos de Mann-Whitney, obtemos um p-valor médio de 0,073495, que é pouco maior que 0,05 para um intervalo de 95% de confiança, quando comparado com o espectro do operador periódico, logo os espectros não são significativamente diferentes. Realizamos 20 simulações para este caso, as quais são apresentadas na Tabela 1.
Finalmente, para o potencial com perturbação ,
O espectro do operador com potencial perturbado (com fator ) no plano complexo é apresentado na Figura 3. O espectro não é normalmente distribuído, com p-valor de no teste de Shapiro-Wilk. Utilizando o teste de hipóteses para dados não-paramétricos de Mann-Whitney, obtemos um p-valor médio de 0,03579, que é pouco menor que 0,05 para um intervalo de 95% de confiança, quando comparado com o espectro do operador periódico, logo os espectros, nesse caso, são significativamente diferentes. Realizamos 20 simulações para este caso, as quais são apresentadas na Tabela 2.
6. Considerações Finais
Fomos capazes de perceber diferenças significativas entre o espectro do operador com potencial periódico e do perturbado, todavia precisamos de uma perturbação da ordem de , com uniformemente distribuída. Para os valores de perturbações menores testados, não podemos distinguir estatisticamente os dois espectros. Para mais informações sobre programação em Mathematica, veja [9][9] P.R. Wellin, R.J. Gaylord e S.N. Kamin, An Introduction to Programming with Mathematica (Cambridge University Press, New York, 2005), 3ª ed..
Agradecimentos
Os autores agradecem ao professor Leonardo de Lima pela revisão minuciosa de parte deste trabalho e pelo incentivo; e ao professor Tobias Heimfarth pelas conversas sobre a implementação computacional.
Referências
-
[1]M. Massimi, Pauli's Exclusion Principle: The Origin and Validation of a Scientific Principle (Cambridge University Press, Cambridge, 2005).
-
[2]D.J. Griffiths, Introduction to Quantum Mechanics (Prentice Hall, New Jersey, 2004), 2ª ed.
-
[3]C.E. Burkhardt e J.J. Leventhal, Foundations of Quantum Physics (Springer, New York, 2008).
-
[4]A.M. Hinz, em: Spectral Theory of Schrödinger Operators, editado por R. del Rio e C. Villegas-Blas (American Mathematical Society, Rhode Island, 2004).
-
[5]F. Bloch, Zeitschrift für Physik 52, 555 (1929).
-
[6]S. Fujita e K. Ito, Quantum Theory of Conducting Matter (Springer, New York, 2007).
-
[7]C. Kittel, Introduction to Solid State Physics (John Wiley and Sons, New Jersey, 1996), 7ª ed.
-
[8]W. Kirsch, in: Random Schrödinger operators (Société Mathématique de France, Paris, 2008), v. 25, p. 1.
-
[9]P.R. Wellin, R.J. Gaylord e S.N. Kamin, An Introduction to Programming with Mathematica (Cambridge University Press, New York, 2005), 3ª ed.
Datas de Publicação
-
Publicação nesta coleção
16 Mar 2020 -
Data do Fascículo
2020
Histórico
-
Recebido
23 Set 2019 -
Revisado
18 Jan 2020 -
Aceito
25 Jan 2020