Efeito Hall Anômalo Quântico em Monocamadas de Favô Ti/Zr/Hf-Bi/Sb: Previsão de Primeiros Princípios

Resumo

Avanços em pesquisa de materiais bidimensionais continuamente revelam oportunidades para novas fases topológicas. Embora o efeito Hall anômalo quântico (QAHE) tenha sido recentemente demonstrado em isolantes topológicos dopados magneticamente, ele permanece confinado a temperaturas criogênicas. Usando teoria do funcional da densidade, prevemos que monocamadas de favô compostas de metais de transição Ti, Zr ou Hf combinados com pnictogênios Bi ou Sb (MPn, M = Ti, Zr, Hf; Pn = Bi, Sb) hospedam QAHE intrínseco. Entre os sistemas estudados, TiBi planar e HfBi planar emergem como isolantes QAH com lacunas de banda de 15 meV e 7 meV, respectivamente, sob modesta tensão de tração. Revelamos como variação de constantes de rede e distâncias de flambagem impulsiona transições de fase topológica, observamos inversões de banda induzidas por acoplamento spin-órbita, e confirmamos modos de borda quirais em nanofitas zigzag via modelos de ligação forte derivados de Wannier. Estes achados expandem o catálogo de materiais QAH realistas e sugerem plataformas viáveis para dispositivos espintrônicos de baixa energia.

Introdução

Isolantes topológicos bidimensionais (2D TIs) exibem estados de borda sem lacuna, polarizados por spin, enquanto mantêm um bulk isolante, uma marca do efeito Hall quântico de spin (QSH). Quando a simetria de reversão temporal (TRS) é quebrada — seja por magnetismo intrínseco ou perturbações externas — canais de borda helicoidais são convertidos em quirais, dando origem ao efeito Hall anômalo quântico (QAHE). O transporte de carga sem dissipação, livre de campo magnético, do QAHE o torna atraente para espintrônica de ultra-baixa potência, mas realizações experimentais têm sido limitadas a isolantes topológicos magnéticos diluídos. Trabalhos teóricos indicam que introduzir ferromagnetismo e forte acoplamento spin-órbita (SOC) em um material QSH pode abrir uma lacuna de banda topológica e gerar número de Chern não-nulo. Enquanto muitos sistemas QSH — como monocamadas do grupo IV (Sn) e grupo V (Bi, Sb) — têm sido extensivamente estudados, o potencial de favôs baseados em metais de transição para QAHE permanece amplamente inexplorado. Aqui, preenchemos esta lacuna explorando favôs MPn (M = Ti, Zr, Hf; Pn = Bi, Sb), baseando-se em previsões recentes de que favôs III-V podem hospedar fases QSH e que dopagem com metal de transição pode induzir magnetismo. Nosso trabalho emprega cálculos de primeiros princípios para avaliar se substituição de Ti, Zr ou Hf em favôs de Bi/Sb pode estabilizar uma fase QAH intrínseca.

Resultados e Discussões

A substituição de metade dos átomos de Bi em uma rede de favô pristina por um metal de transição (Ti, Zr ou Hf) rende monocamadas MPn que podem adotar geometrias de baixo flambagem ou planas. A Figura 1a mostra a vista superior da rede M-Pn, enquanto Figuras 1b e 1c ilustram vistas laterais das configurações flambadas e planas, respectivamente. A correspondente primeira zona de Brillouin (BZ) com pontos de alta simetria é mostrada na Figura 1d.

a Estrutura cristalina de favô M-Sb/Bi. b, c Vistas laterais de estruturas flambadas e planas, respectivamente. d Primeira zona de Brillouin (BZ) com pontos de alta simetria

Avaliamos estabilidade estrutural e efeitos de tensão variando a constante de rede e relaxando posições atômicas para ambas as fases flambada e planar. Curvas de energia total para TiBi, ZrBi e HfBi são plotadas na Figura 2a-c. Estruturas MPn de baixo flambagem consistentemente possuem energias menores que suas contrapartes planas. No entanto, o QAHE aparece predominantemente no regime planar tensionado, com apenas uma estreita janela de constante de rede suportando uma fase QAH em HfBi flambado (ver Fig. 2c).

Diagrama de fase de a TiBi, b ZrBi e c HfBi mostrando energia total em diferentes constantes de rede. O diagrama é dividido em várias regiões rotuladas como QAH (fase Hall anômalo quântico), I (isolante) e SM (semi-metal). Círculos azuis e triângulos vermelhos representam casos flambados e planares, respectivamente

A Tabela 1 lista as constantes de rede de equilíbrio, momentos magnéticos, lacunas de banda e números de Chern calculados para todos os sistemas MPn. Números de Chern inteiros não-nulos identificam isolantes QAH: TiBi planar (C = +1, Δ = 15 meV) e HfBi planar (C = –1, Δ = 7 meV) são os candidatos mais promissores. A transição de fase em TiBi pode ser impulsionada pela redução da distância de flambagem (Fig. 3), enquanto tensão de tração induz transição em HfBi (Fig. 4).

Transição de fase após variar distância flambada. a Diagrama de fase de TiBi em a=4,6 Å. A seta mostra o caminho da transição. b-f Transição da estrutura de banda conforme distância de flambagem (δ) foi reduzida de 0,44 a 0,4 Å. A transição ocorre em δ=0,41 Å

Transição de fase após variar constante de rede. a Diagrama de fase de HfBi flambado. A seta mostra o caminho da transição. b-h Transição da estrutura de banda conforme constante de rede foi aumentada de 4,7 a 5,1 Å

Estruturas de banda eletrônica para geometrias de equilíbrio são apresentadas na Figura 5. Na geometria planar, TiBi exibe clara inversão de banda e lacuna de 15 meV, enquanto HfBi planar mostra lacuna de 7 meV com C = –1. HfBi flambado, no entanto, é um semi-metal com alto número de Chern (C = –3). Outras configurações, como ZrBi flambado e TiSb, permanecem isolantes triviais.

Estruturas de banda eletrônica de M-Pn (M=Ti, Zr e Hf; Pn=Sb e Bi) em suas constantes de rede de equilíbrio para a planar e b casos flambados. Círculos vermelhos e azuis indicam contribuições +s_z e −s_z, respectivamente

Para dissecar o papel de SOC e magnetismo, focamos em TiBi planar (a = 4,76 Å). A Figura 6 mostra que o cálculo não-magnético rende estado metálico, enquanto ordenação ferromagnética (momento magnético ≈ 1,05 µB por célula, dominado por Ti) produz espectro semi-metálico. Inclusão de SOC abre lacuna de 15 meV, confirmando que inversão de banda é impulsionada por SOC e responsável pelo QAHE.

Estruturas de banda eletrônica de filme TiBi planar em a=4,76 Å para cálculos não-magnéticos (a) sem SOC e (c) com SOC, bem como cálculos ferromagnéticos (b) sem SOC e (d) com SOC. Círculos vermelhos e azuis indicam contribuições +s_z e −s_z, respectivamente, para (c) não-magnético (d) cálculos ferromagnéticos com SOC

Estados de borda são impressões digitais essenciais de uma fase topológica. Usando Hamiltonianos de ligação forte derivados de Wannier, construímos uma nanofita zigzag de HfBi de 127 Å de largura (Fig. 7). A estrutura de banda revela modos de borda quirais cruzando o nível de Fermi um número ímpar de vezes, igual a |C| = 1, confirmando a natureza QAH do HfBi planar.

Estrutura de banda ao longo da borda de nanofita zigzag de HfBi flambado com a=4,9 Å e largura de 127 Å. Círculos azuis (vermelhos) indicam contribuição das bordas esquerda (direita). Bandas de bulk são denotadas pela região preenchida laranja

Cálculos de fônon mostram pequenas frequências imaginárias para monocamadas livres, sugerindo que um substrato pode ser necessário para realização experimental. Além disso, dentro de uma supercélula 2×2, a configuração ferromagnética permanece energeticamente favorável comparada a ordem antiferromagnética em casos flambados, enquanto sistemas planares exibem quase-degenerescência entre estados FM e AFM.

Conclusões

Nosso estudo de primeiros princípios demonstra que substituição de Ti, Zr ou Hf em favôs de Bi ou Sb pode induzir QAHE intrínseco. Enquanto geometrias flambadas são energeticamente preferidas, modesta tensão de tração ou redução da distância de flambagem impulsiona transição para configuração planar que hospeda número de Chern não-nulo. TiBi planar e HfBi planar emergem como robustos isolantes QAH com lacunas de banda de 15 meV e 7 meV, respectivamente. Estes materiais adicionam ao crescente rol de candidatos QAH bidimensionais e são plataformas promissoras para dispositivos espintrônicos de baixa potência.

Métodos/Experimental

Todos os cálculos foram realizados dentro da estrutura da teoria do funcional da densidade usando aproximação de gradiente generalizada (GGA) e potenciais de onda aumentada projetada (PAW) como implementado no VASP 5.3. O corte de onda plana foi definido para 350 eV, e estruturas foram relaxadas até forças residuais abaixo de 5×10⁻³ eV/Å. Convergência do loop de auto-consistência eletrônica foi atingida a 10⁻⁶ eV. Espaçamento de vácuo de 20 Å ao longo do eixo z preveniu interações intercamada, e uma grade Monkhorst-Pack 24×24×1 centrada em Gamma amostrou a zona de Brillouin 2D. Funções de Wannier maximamente localizadas foram geradas com WANNIER90, permitindo construção de modelos de ligação forte e cálculo de estados de borda. Invariantes topológicos foram obtidos via pacote Z2Pack, que rastreia centros de carga de Wannier híbridos para render o número de Chern.