Projeto Computacional de Materiais com Banda Plana para Aplicações em Espintrônica

Resumo

A mecânica quântica estabelece que a integral de salto entre orbitais locais torna as bandas de energia dispersivas. No entanto, em casos especiais, surgem bandas sem dispersão devido à interferência quântica — as chamadas bandas planas. Diversos modelos teóricos preveem propriedades físicas fascinantes para esses sistemas, mas nenhum composto real com banda plana foi encontrado em 25 anos de pesquisa intensa.

Neste trabalho, demonstramos por cálculos de primeiros princípios que alguns óxidos de pirocloro possuem quase-banda plana logo abaixo do nível de Fermi. Suas bandas de valência são bem descritas por um modelo de ligação forte (tight-binding) na rede de pirocloro com integral de salto isotrópica entre vizinhos mais próximos. Esse modelo pertence à classe de Mielke, cujo estado fundamental é ferromagnético com dopagem adequada de portadores e interação de Coulomb intra-atômica repulsiva.

Cálculos de banda com polarização de spin para o sistema dopado com lacunas confirmam o estado fundamental ferromagnético em certa faixa de dopagem. Notavelmente, esses compostos não contêm elementos magnéticos (metais de transição ou terras-raras), o que representa um avanço conceitual significativo.

Introdução

As propriedades eletrônicas e magnéticas dos materiais são amplamente determinadas pela sua dispersão de energia. Condutividade elevada requer grande dispersão de banda; já o magnetismo surge quando a dispersão é pequena. Geralmente, a dispersão depende do caráter da função de onda atômica, por isso a maioria dos materiais magnéticos contém elementos magnéticos.

Se conseguirmos sintetizar um material magnético sem elementos magnéticos, o impacto será imensurável. Propomos aqui candidatos a ferromagnetos sem elementos magnéticos, usando cálculos de primeiros princípios. A largura de banda de orbitais de átomos não magnéticos costuma ser pequena, mas em certos casos torna-se extremamente reduzida — a banda plana. Se o nível de Fermi estiver exatamente sobre essa banda, um estado fundamental ferromagnético torna-se possível.

Modelo de Ligação Forte e Rede de Pirocloro

Considere o Hamiltoniano simples de ligação forte:

H₀ = ε Σ c_i⁺c_i + Σ t_ijc_i⁺c_j

onde c_i é o operador de aniquilação no sítio i, ε a energia no sítio, e t_ij a integral de salto (finita e isotrópica = t apenas para vizinhos mais próximos). Em redes simples (ex.: quadrada), a dispersão é E(k) = ε + 2t(cosk_x + cosk_y), com largura de banda W = 8t ∝ t. Como t depende da sobreposição de funções de onda atômicas, bandas formadas por orbitais s ou p são largas, desfavorecendo o magnetismo.

Contudo, em redes específicas como a rede de pirocloro (sublattice A da estrutura A₂B₂O₇), aparecem bandas duplamente degeneradas e sem dispersão. Isso pode ser provado matematicamente. Outras redes geradoras de bandas planas incluem a rede quadriculada 2D e a rede de kagome 2D.

Interessantemente, prova-se que se essa banda plana estiver meio-preenchida, o sistema possui estado fundamental ferromagnético único para qualquer valor positivo da interação de Coulomb intra-atômica U. Essas redes são conhecidas como "redes com frustração geométrica" no contexto de spins localizados. De fato, óxidos de pirocloro R₂Ti₂O₇ (R = terra-rara) exibem propriedades magnéticas exóticas: líquido de spin quântico, gelo de spin, monopolo magnético.

Condições para Banda Plana em Óxidos de Pirocloro

Para que a banda plana apareça no topo da banda de valência em A₂B₂O₇, três condições são necessárias:

1. O íon B deve ser inerte (camada fechada: nd⁰ ou np⁰);
2. A integral de salto deve ser isotrópica, exigindo que o nível de Fermi esteja no orbital s do sítio A;
3. A banda plana deve se formar no sublattice A (rede de pirocloro).

Isso leva a três famílias de candidatos:

• (a) A¹⁺₂B⁶⁺₂O₇ → A = Tl¹⁺; B = Mo⁶⁺, W⁶⁺
• (b) A²⁺₂B⁵⁺₂O₇ → A = Sn²⁺, Pb²⁺; B = Nb⁵⁺, Ta⁵⁺
• (c) A³⁺₂B⁴⁺₂O₇ → A = Bi³⁺; B = Ti⁴⁺, Sn⁴⁺

Todos os íons A têm configuração (5s)² ou (6s)². Focamos em três compostos: Tl₂Mo₂O₇, Sn₂Nb₂O₇ e Bi₂Ti₂O₇. Os dois últimos já foram sintetizados; o análogo Tl₂Ru₂O₇ também existe e mostra transição metal-isolante única.

Métodos

Calculamos a estrutura eletrônica dos três compostos por primeiros princípios (DFT), assumindo estrutura de pirocloro ideal A₂B₂O₆O′ (grupo espacial Fd-3m). Usamos o método FLAPW (onda plana de potencial total) com aproximação de gradiente generalizado (GGA) e o pacote WIEN2k. Parâmetros: RK_max = 7.0, malha k com ~1000 pontos na zona de Brillouin. A estrutura foi otimizada mantendo a simetria. Para Sn₂Nb₂O₇ e Bi₂Ti₂O₇, usamos parâmetros de rede experimentais; para Tl₂Mo₂O₇, otimizamos a = 10.517 Å (próximo ao Tl₂Ru₂O₇ experimental).

Resultados e Discussão

Estrutura de Bandas

A Figura 2 mostra as dispersões de banda. Em Sn₂Nb₂O₇, o topo da banda de valência (−3 a 0 eV) reproduz fielmente o modelo de ligação forte da Figura 1b — surpreendente, pois o modelo usa apenas dois parâmetros (ε e t). A banda de valência é bem descrita por orbitais "Sn-s" (na verdade, orbitais anti-ligantes Sn-s/O′-p). A principal diferença é a planicidade da banda perto de 0 eV, indicando que integrais de salto além dos vizinhos mais próximos são necessárias para ajuste preciso.

Em Tl₂Mo₂O₇, a banda de valência é quase idêntica à do Sn₂Nb₂O₇, confirmando a banda plana. Porém, a banda de condução desce e o gap colapsa: a banda Mo-d é parcialmente ocupada, sugerindo que a configuração iônica formal Tl¹⁺₂Mo⁶⁺₂O²⁻₇ não é adequada. A banda plana Tl-s se entrelaça com a banda Mo-d, como no análogo Tl₂Ru₂O₇, onde a transição metal-isolante é atribuída à banda plana Tl-s "oculta". Espera-se fenômeno semelhante no Tl₂Mo₂O₇ se sintetizado.

Em Bi₂Ti₂O₇, o topo da banda de valência é quase plano na maior parte dos eixos de simetria, embora a forma geral difira do modelo de pirocloro. Ainda assim, a quase-banda plana e o alto densidade de estados (DOS) resultante permitem esperar ferromagnetismo com dopagem de lacunas.

Estados Ferromagnéticos

Para induzir ferromagnetismo nos isolantes Sn₂Nb₂O₇ e Bi₂Ti₂O₇, introduzimos lacunas na quase-banda plana. No caso de banda plana perfeita, qualquer U > 0 causa ferromagnetismo com banda meio-preenchida. Para quase-banda plana, estudo numérico mostra que U > U_c (ordem da largura de banda W) induz ferromagnetismo. Como estimar U e U_c em compostos reais é difícil, realizamos cálculos ab initio com polarização de spin.

Simulamos dopagem com lacunas substituindo N por O′ (Sn₂Nb₂O₆N e Bi₂Ti₂O₆N), reduzindo dois elétrons por célula primitiva (um por unidade fórmula), meio-preenchendo a quase-banda plana.

A Figura 3 mostra as curvas DOS: a quase-banda plana forma pico agudo no nível de Fermi. Ambos compostos tornam-se semi-metálicos (spin totalmente polarizado em E_F). Momento magnético total M = 2.00 μ_B/célula primitiva, confirmando polarização total. A divisão de troca é ~0,3 eV (Sn₂Nb₂O₆N) e ~0,4 eV (Bi₂Ti₂O₆N) — menor que no Fe bcc (~2 eV), mas excede a largura de banda, estabilizando o estado semi-metálico.

Conclusões

Estabelecemos um princípio orientador para projetar compostos com banda plana. Seguindo-o, selecionamos três óxidos de pirocloro e investigamos sua estrutura eletrônica por primeiros princípios. Combinado com análise de ligação forte, mostramos que alguns compostos possuem quase-banda plana real. Demontramos também que a dopagem com lacunas leva a estado fundamental ferromagnético, apesar da ausência de elementos magnéticos. Esses achados são um grande passo para realizar não apenas sistemas de banda plana em compostos, mas também ferromagnetos sem elementos magnéticos.

Abreviações

DOS:

Densidade de estados

FLAPW:

Onda plana de potencial total (Full-potential augmented plane wave)

TB:

Ligação forte (Tight-binding)