Nome do Projeto
Estudo Analítico e Numérico de Estruturas de Bandas e Espalhamento Quântico em Potenciais Periódicos Unidimensionais
Ênfase
Pesquisa
Data inicial - Data final
11/08/2026 - 10/08/2030
Unidade de Origem
Coordenador Atual
Área CNPq
Ciências Exatas e da Terra
Resumo
A análise de potenciais periódicos na equação de Schrödinger permite compreender a estrutura de bandas de energia e as lacunas proibidas em sistemas quânticos. Embora o modelo tradicional de Kronig-Penney obtenha soluções analíticas exatas ao recorrer ao potencial Delta de Dirac, a descrição de arranjos mais complexos exige métodos numéricos aplicáveis a perfis de potencial arbitrários. Este projeto visa desenvolver um modelo computacional em Python, empregando o Método das Diferenças Finitas (MDF) com condições de contorno de Bloch, para resolver a equação de Schrödinger unidimensional em redes periódicas. O trabalho envolve o mapeamento da primeira zona de Brillouin e o cálculo de autovalores e autovetores, com foco na modelagem de poços de potencial generalizados e no estudo de seu comportamento assintótico à medida que tendem ao limite do potencial Delta de Dirac. Para garantir a precisão dos algoritmos, os resultados numéricos serão validados por testes de convergência e confrontados diretamente com as equações transcendentais analíticas tradicionais.
Objetivo Geral
Desenvolver um modelo computacional em Script ou linguagem de programação para resolver a equação de Schrodinger unidimensional em potenciais periódicos utilizando o Método das Diferenças Finitas (MDF), analisando a evolução das estruturas de bandas de energia e lacunas proibidas a medida que os poços de potencial contínuos tendem assintoticamente ao limite do potencial Delta de Dirac.
Justificativa
A análise do comportamento quântico de partículas sob potenciais periódicos é a base para entender a estrutura de bandas e as propriedades de transporte na matéria condensada [1, 2]. O modelo de Kronig-Penney com barreiras retangulares e potenciais Delta de Dirac, tornou-se o ponto de partida clássico para descrever lacunas de energia (band gaps) e relações de dispersão E(k) em redes cristalinas [1-3].
Apesar do modelo tradicional de Kronig-Penney seja conceitualmente simples, a sua importância reside em sua capacidade de servir como ponto de partida para a descrição de sistemas físicos reais mais complexos. Estudo recentes da estrutura de bandas em redes 1D expandiu-se além dos livros-texto, encontrando aplicação direta em sistemas reais de semicondutores, superredes (como junções GaAs/AlGaAs) e, mais recentemente, em simulações quânticas com átomos frios e cadeias híbridas de átomos e íons em armadilhas ópticas [4-6].
Este é um tema de interesse atual demonstrada pela intensa produção científica nos últimos anos [7-16]. A busca por descrições mais realistas motivou a extensão das equações transcendentais de energia para poços de potenciais assimétricas e acoplados, permitindo calcular analiticamente grandezas de transporte como velocidade de grupo e massa efetiva de banda [7, 8]. Do ponto de vista computacional, metodologias baseados na matriz de transferência [9], no método das diferenças finitas (FDM) [10] e na mecânica matricial em base de ondas planas [2] viabilizaram a resolução numérica precisa de estruturas de bandas em células unitárias complexas, bem como a avaliação do impacto de impurezas e desordem na rede [11].
A generalização de potenciais periódicos 1D facilitou a descrição de fenômenos quânticos emergentes. No âmbito dos sistemas 1D sob confinamento harmônico transversal, a consideração de comprimentos de espalhamento dependentes da energia permitiu descrever a estrutura de bandas de experimentos reais em guias de onda atômicos e simuladores quânticos átomo-íon [4,6]. Da mesma forma, a investigação de redes periódicas atrativas revelou a existência de bandas de energia negativa, estados ligados altamente degenerados e assimetrias marcantes entre as massas efetivas de elétrons e buracos [12]. Outro avanço envolve a introdução de dimensões virtuais para mapear a estrutura de bandas em invariantes topológicos [13], o que expõe estados de borda quiralmente protegidos, modos planos e espectros fractais [11]. Além disso, a aplicação do potencial de Kronig-Penney a gases quânticos em limites de forte correlação desvelou a abertura de lacunas nas excitações de Lieb e o aparecimento de modos planos topológicos em sua função espectral [14]. Mais recentemente, investigações sobre a equação de Schrödinger fracionária para potenciais periódicos 1D demonstraram que a estrutura de bandas em mídias quânticas periódicas pode ser sintonizada e até invertida mediante alterações no índice de Lévy do operador quântico [15].
Do ponto de vista pedagógico, a integração de simulações numéricas com demostrações analíticas didáticas reforça o valor educativo do modelo, preenchendo a lacuna entre os conceitos básicos de mecânica quântica e a pesquisa de ponta em física da matéria condensada.
Diante desse cenário, o projeto justifica-se pela oportunidade de integrar o desenvolvimento analítico e computacional de potenciais periódicos 1D, investigando transporte quântico, topologia e sistemas correlacionados, além de disponibilizar ferramentas numéricas de código aberto para a comunidade.
REFERÊNCIAS BIBLIOGRÁFICAS
[1] KRONIG, R. de L.; PENNEY, W. G. Quantum mechanics of electrons in crystal lattices. Proceedings of the Royal Society of London. Series A, v. 130, n. 814, p. 499–513, 1931.
[2] PAVELICH, R. L.; MARSIGLIO, F. The Kronig-Penney model extended to arbitrary potentials via numerical matrix mechanics. American Journal of Physics, v. 83, n. 9, p. 773–781, 2015.
[3] BELLONI, M.; ROBINETT, R. W. The infinite well and Dirac delta function potentials as pedagogical, mathematical and physical models in quantum mechanics. Physics Reports, v. 540, n. 2, p. 25–122, 2014.
[4] NEGRETTI, A. et al. Generalized Kronig-Penney model for ultracold atomic quantum systems. Physical Review B, v. 90, n. 15, p. 155426, 2014.
[5] SROCZYŃSKA, M.; WASAK, T.; IDZIASZEK, Z. Analytically solvable quasi-one-dimensional Kronig-Penney model. Physical Review A, v. 111, n. 2, p. 023313, 2025.
[6] KESHAVARZ, E.; MAZHARIMOUSAVI, S. H. One-dimensional quantum scattering from multiple Dirac δ-potentials: a Python-based solution. Physica Scripta, v. 99, n. 7, p. 075229, 2024.
[7] MENDOZA-VILLA, F.; RAMOS-GUIVAR, J. A.; ESPINOZA-BERNARDO, R. M. Generalized One-Dimensional Periodic Potential Wells Tending to the Dirac Delta Potential. Physics, v. 6, n. 1, p. 75–93, 2024.
[8] MENDOZA-VILLA, F. et al. Towards the Analytical Generalization of the Transcendental Energy Equation, Group Velocity, and Effective Mass in One-Dimensional Periodic Potential Wells with a Computational Application to Common Coupled Potentials. Applied Sciences, v. 14, n. 10, p. 3987, 2024.
[9] FIGUEROA, J.; GONZALEZ, I.; SALINAS-ARIZMENDI, D. A novel transfer matrix framework for multiple Dirac delta potentials. Physics Letters A, v. 555, p. 130785, 2025.
[10] LI, L.; HU, J.; ZHANG, Y. One-dimension periodic potentials in Schrödinger equation solved by the finite difference method. European Journal of Physics, v. 46, n. 1, p. 015402, 2025.
[11] RESHODKO, I. et al. Topological states in the Kronig–Penney model with arbitrary scattering potentials. New Journal of Physics, v. 21, n. 1, p. 013010, 2019.
[12] JOHNSTON, D. C. Attractive Kronig-Penney Band Structures and Wave Functions. arXiv preprint, arXiv:1905.12084 [cond-mat.mtrl-sci], 2019.
[13] ZHENG, Y.; YANG, S.-J. Topological bands in one-dimensional periodic potentials. Physica B: Condensed Matter, v. 454, p. 93–97, 2014.
[14] HE, W.-B. et al. Lieb excitations and topological flat mode of the spectral function of a Tonks-Girardeau gas in a Kronig-Penney potential. Physical Review A, v. 111, n. 1, p. 013312, 2025.
[15] GUYETTE, B. R.; LEWIS, J. M.; CARR, L. D. Tunable Bands in One-Dimensional Fractional Quantum Media. Physical Review E, v. 114, n. 1, p. 014106, 2026.
Apesar do modelo tradicional de Kronig-Penney seja conceitualmente simples, a sua importância reside em sua capacidade de servir como ponto de partida para a descrição de sistemas físicos reais mais complexos. Estudo recentes da estrutura de bandas em redes 1D expandiu-se além dos livros-texto, encontrando aplicação direta em sistemas reais de semicondutores, superredes (como junções GaAs/AlGaAs) e, mais recentemente, em simulações quânticas com átomos frios e cadeias híbridas de átomos e íons em armadilhas ópticas [4-6].
Este é um tema de interesse atual demonstrada pela intensa produção científica nos últimos anos [7-16]. A busca por descrições mais realistas motivou a extensão das equações transcendentais de energia para poços de potenciais assimétricas e acoplados, permitindo calcular analiticamente grandezas de transporte como velocidade de grupo e massa efetiva de banda [7, 8]. Do ponto de vista computacional, metodologias baseados na matriz de transferência [9], no método das diferenças finitas (FDM) [10] e na mecânica matricial em base de ondas planas [2] viabilizaram a resolução numérica precisa de estruturas de bandas em células unitárias complexas, bem como a avaliação do impacto de impurezas e desordem na rede [11].
A generalização de potenciais periódicos 1D facilitou a descrição de fenômenos quânticos emergentes. No âmbito dos sistemas 1D sob confinamento harmônico transversal, a consideração de comprimentos de espalhamento dependentes da energia permitiu descrever a estrutura de bandas de experimentos reais em guias de onda atômicos e simuladores quânticos átomo-íon [4,6]. Da mesma forma, a investigação de redes periódicas atrativas revelou a existência de bandas de energia negativa, estados ligados altamente degenerados e assimetrias marcantes entre as massas efetivas de elétrons e buracos [12]. Outro avanço envolve a introdução de dimensões virtuais para mapear a estrutura de bandas em invariantes topológicos [13], o que expõe estados de borda quiralmente protegidos, modos planos e espectros fractais [11]. Além disso, a aplicação do potencial de Kronig-Penney a gases quânticos em limites de forte correlação desvelou a abertura de lacunas nas excitações de Lieb e o aparecimento de modos planos topológicos em sua função espectral [14]. Mais recentemente, investigações sobre a equação de Schrödinger fracionária para potenciais periódicos 1D demonstraram que a estrutura de bandas em mídias quânticas periódicas pode ser sintonizada e até invertida mediante alterações no índice de Lévy do operador quântico [15].
Do ponto de vista pedagógico, a integração de simulações numéricas com demostrações analíticas didáticas reforça o valor educativo do modelo, preenchendo a lacuna entre os conceitos básicos de mecânica quântica e a pesquisa de ponta em física da matéria condensada.
Diante desse cenário, o projeto justifica-se pela oportunidade de integrar o desenvolvimento analítico e computacional de potenciais periódicos 1D, investigando transporte quântico, topologia e sistemas correlacionados, além de disponibilizar ferramentas numéricas de código aberto para a comunidade.
REFERÊNCIAS BIBLIOGRÁFICAS
[1] KRONIG, R. de L.; PENNEY, W. G. Quantum mechanics of electrons in crystal lattices. Proceedings of the Royal Society of London. Series A, v. 130, n. 814, p. 499–513, 1931.
[2] PAVELICH, R. L.; MARSIGLIO, F. The Kronig-Penney model extended to arbitrary potentials via numerical matrix mechanics. American Journal of Physics, v. 83, n. 9, p. 773–781, 2015.
[3] BELLONI, M.; ROBINETT, R. W. The infinite well and Dirac delta function potentials as pedagogical, mathematical and physical models in quantum mechanics. Physics Reports, v. 540, n. 2, p. 25–122, 2014.
[4] NEGRETTI, A. et al. Generalized Kronig-Penney model for ultracold atomic quantum systems. Physical Review B, v. 90, n. 15, p. 155426, 2014.
[5] SROCZYŃSKA, M.; WASAK, T.; IDZIASZEK, Z. Analytically solvable quasi-one-dimensional Kronig-Penney model. Physical Review A, v. 111, n. 2, p. 023313, 2025.
[6] KESHAVARZ, E.; MAZHARIMOUSAVI, S. H. One-dimensional quantum scattering from multiple Dirac δ-potentials: a Python-based solution. Physica Scripta, v. 99, n. 7, p. 075229, 2024.
[7] MENDOZA-VILLA, F.; RAMOS-GUIVAR, J. A.; ESPINOZA-BERNARDO, R. M. Generalized One-Dimensional Periodic Potential Wells Tending to the Dirac Delta Potential. Physics, v. 6, n. 1, p. 75–93, 2024.
[8] MENDOZA-VILLA, F. et al. Towards the Analytical Generalization of the Transcendental Energy Equation, Group Velocity, and Effective Mass in One-Dimensional Periodic Potential Wells with a Computational Application to Common Coupled Potentials. Applied Sciences, v. 14, n. 10, p. 3987, 2024.
[9] FIGUEROA, J.; GONZALEZ, I.; SALINAS-ARIZMENDI, D. A novel transfer matrix framework for multiple Dirac delta potentials. Physics Letters A, v. 555, p. 130785, 2025.
[10] LI, L.; HU, J.; ZHANG, Y. One-dimension periodic potentials in Schrödinger equation solved by the finite difference method. European Journal of Physics, v. 46, n. 1, p. 015402, 2025.
[11] RESHODKO, I. et al. Topological states in the Kronig–Penney model with arbitrary scattering potentials. New Journal of Physics, v. 21, n. 1, p. 013010, 2019.
[12] JOHNSTON, D. C. Attractive Kronig-Penney Band Structures and Wave Functions. arXiv preprint, arXiv:1905.12084 [cond-mat.mtrl-sci], 2019.
[13] ZHENG, Y.; YANG, S.-J. Topological bands in one-dimensional periodic potentials. Physica B: Condensed Matter, v. 454, p. 93–97, 2014.
[14] HE, W.-B. et al. Lieb excitations and topological flat mode of the spectral function of a Tonks-Girardeau gas in a Kronig-Penney potential. Physical Review A, v. 111, n. 1, p. 013312, 2025.
[15] GUYETTE, B. R.; LEWIS, J. M.; CARR, L. D. Tunable Bands in One-Dimensional Fractional Quantum Media. Physical Review E, v. 114, n. 1, p. 014106, 2026.
Metodologia
A execução do projeto está organizada em quatro etapas integradas:
Etapa 1: Discretização do Sistema e Condições de Bloch
O domínio espacial da célula unitária será dividido em uma grade uniforme de N pontos. A derivada segunda da função de onda na equação de Schrödinger independente do tempo será aproximada pelo operador de diferença finita central de segunda ordem. Para aplicar as condições de contorno de Bloch, os pontos extremos da malha serão acoplados, convertendo a equação diferencial em um problema de autovalores associado a uma matriz Hamiltoniana Hermitiana e complexa.
Etapa 2: Implementação das Rotinas em Python
A montagem e a solução das matrizes serão feitas em Python, utilizando os pacotes NumPy e SciPy. O algoritmo construirá a matriz Hamiltoniana com base no número de pontos N da malha, e a extração dos autovalores (energias) e autovetores (funções de onda) será realizada por algoritmos de diagonalização dedicados a matrizes Hermitianas.
Etapa 3: Modelação de Potenciais e Limite do Delta de Dirac
Nesta etapa, serão codificados diferentes perfis de potenciais periódicos. Para estudar a aproximação em direção ao potencial Delta de Dirac, a largura da barreira será reduzida progressivamente enquanto sua altura cresce, mantendo constante a área integrada do potencial. O cálculo da estrutura de bandas será realizado efetuando-se a varredura do vetor de onda k ao longo de toda a primeira zona de Brillouin.
Etapa 4: Análise de Convergência e Validação dos Resultados
A precisão dos resultados obtidos por diferenças finitas será avaliada de duas formas. Primeiro, se fará o refinamento da malha (variando N) para testar a estabilidade e a convergência das energias. Em segundo lugar, os autovalores numéricos serão confrontados diretamente com as soluções exatas obtidas pelas equações transcendentais do modelo de Kronig-Penney clássico, quantificando os erros residuais do método.
Etapa 1: Discretização do Sistema e Condições de Bloch
O domínio espacial da célula unitária será dividido em uma grade uniforme de N pontos. A derivada segunda da função de onda na equação de Schrödinger independente do tempo será aproximada pelo operador de diferença finita central de segunda ordem. Para aplicar as condições de contorno de Bloch, os pontos extremos da malha serão acoplados, convertendo a equação diferencial em um problema de autovalores associado a uma matriz Hamiltoniana Hermitiana e complexa.
Etapa 2: Implementação das Rotinas em Python
A montagem e a solução das matrizes serão feitas em Python, utilizando os pacotes NumPy e SciPy. O algoritmo construirá a matriz Hamiltoniana com base no número de pontos N da malha, e a extração dos autovalores (energias) e autovetores (funções de onda) será realizada por algoritmos de diagonalização dedicados a matrizes Hermitianas.
Etapa 3: Modelação de Potenciais e Limite do Delta de Dirac
Nesta etapa, serão codificados diferentes perfis de potenciais periódicos. Para estudar a aproximação em direção ao potencial Delta de Dirac, a largura da barreira será reduzida progressivamente enquanto sua altura cresce, mantendo constante a área integrada do potencial. O cálculo da estrutura de bandas será realizado efetuando-se a varredura do vetor de onda k ao longo de toda a primeira zona de Brillouin.
Etapa 4: Análise de Convergência e Validação dos Resultados
A precisão dos resultados obtidos por diferenças finitas será avaliada de duas formas. Primeiro, se fará o refinamento da malha (variando N) para testar a estabilidade e a convergência das energias. Em segundo lugar, os autovalores numéricos serão confrontados diretamente com as soluções exatas obtidas pelas equações transcendentais do modelo de Kronig-Penney clássico, quantificando os erros residuais do método.
Indicadores, Metas e Resultados
Indicadores de Avaliação:
• Número de scripts e pacotes computacionais em Python desenvolvidos, testados e validados.
• Número de artigos submetidos e aceitos para publicação em periódicos científicos qualificados pelo sistema Qualis/Capes.
• Número de resumos publicados em anais de congressos e apresentações realizadas pela equipe.
• Número de orientações de estudantes de iniciação científica e/ou pôs graduação.
Metas Cientificas:
• Concluir a programação do solver quântico e sua validação matemática frente ao modelo analítico de Kronig-Penney ate o final do segundo ano do projeto.
• Publicar ao menos dois artigos em revistas científicas indexadas entre o terceiro e o quarto ano do projeto.
• Levar pelo menos três trabalhos para congressos nacionais ou internacionais nas áreas de física ou modelagem computacional.
Resultados Esperados:
• Espera-se obter um algoritmo em Python, estável e documentado, capaz de resolver a equação de Schrodinger para potenciais periódicos arbitrários por Diferenças Finitas.
• Espera-se gerar um conjunto de dados gráficos e numéricos que demonstrem o comportamento assintótico das bandas de energia e funções de onda no limite de transição para o potencial Delta de Dirac.
• O projeto servira como base para a inserção e treinamento de estudantes em técnicas de computação científica, programação em Python aplicada a física e métodos numéricos avançados.
• Número de scripts e pacotes computacionais em Python desenvolvidos, testados e validados.
• Número de artigos submetidos e aceitos para publicação em periódicos científicos qualificados pelo sistema Qualis/Capes.
• Número de resumos publicados em anais de congressos e apresentações realizadas pela equipe.
• Número de orientações de estudantes de iniciação científica e/ou pôs graduação.
Metas Cientificas:
• Concluir a programação do solver quântico e sua validação matemática frente ao modelo analítico de Kronig-Penney ate o final do segundo ano do projeto.
• Publicar ao menos dois artigos em revistas científicas indexadas entre o terceiro e o quarto ano do projeto.
• Levar pelo menos três trabalhos para congressos nacionais ou internacionais nas áreas de física ou modelagem computacional.
Resultados Esperados:
• Espera-se obter um algoritmo em Python, estável e documentado, capaz de resolver a equação de Schrodinger para potenciais periódicos arbitrários por Diferenças Finitas.
• Espera-se gerar um conjunto de dados gráficos e numéricos que demonstrem o comportamento assintótico das bandas de energia e funções de onda no limite de transição para o potencial Delta de Dirac.
• O projeto servira como base para a inserção e treinamento de estudantes em técnicas de computação científica, programação em Python aplicada a física e métodos numéricos avançados.
Equipe do Projeto
| Nome | CH Semanal | Data inicial | Data final |
|---|---|---|---|
| JAVIER ANTONIO GOMEZ ROMERO | 14 | ||
| JULIA OLIVEIRA NOBRE | |||
| LARA FERNANDA SANT ANA RIBEIRO |