Anatomy of Spin--Orbit Torques in Monolayer Fe$_3$GeTe$_2$ and Fe$_3$GaTe$_2$: Insights from atomistic and momentum-space decompositions

Resumo Técnico: Anatomia dos Torques de Spin–Órbita em Monocamadas de Fe$_3$GeTe$_2$ e Fe$_3$GaTe$_2$

Definição do Problema
O controle elétrico da magnetização via torques de spin–órbita (SOT) é um objetivo central na espintrônica. Embora os ímãs de van der Waals (vdW) como o Fe$_3$GeTe$_2$ (FGT) e o Fe$_3$GaTe$_2$ (FGaT) ofereçam oportunidades únicas devido à sua baixa simetria cristalina e densidade de portadores sintonizável por porta (gate), as origens microscópicas de seus torques induzidos por corrente permanecem incompletamente compreendidas. Especificamente, a dependência dos harmônicos de torque em relação ao preenchimento de banda, o papel de estados eletrônicos específicos próximos ao nível de Fermi ($\varepsilon_F$) e a existência de torques "escondidos" decorrentes de graus de liberdade de sublatice em monocamadas não centrosimétricas requerem investigação sistemática. Embora o FGT e o FGaT compartilhem a mesma estrutura cristalina $D_{3h}$ e caráter dominante de Fe 3d, eles diferem por um elétron de valência (Ga vs. Ge), tornando efetivamente o FGaT um análogo dopado com lacunas (holes) do FGT. O artigo aborda como essa modificação estequiométrica altera o cenário de SOT sem a necessidade de controle eletrostático por porta.

Metodologia
Os autores empregam uma abordagem sistemática de primeira ordem combinando a teoria do funcional da densidade (DFT) com a teoria de resposta linear:

• Estrutura Eletrônica: Os cálculos foram realizados utilizando o Quantum ESPRESSO com pseudopotenciais de Vanderbilt otimizados (ONCV) totalmente relativísticos. O acoplamento spin-órbita (SOC) foi tratado de forma autoconsistente dentro de um arcabouço não colinear. O funcional Perdew–Wang LDA foi utilizado para os cálculos de resposta para garantir momentos magnéticos precisos, consistente com estudos anteriores sobre a monocamada de FGT.
• Interpolação de Wannier: Para permitir uma amostragem densa da zona de Brillouin, funções de Wannier adaptadas à simetria (SAWFs) foram construídas usando o WannierBerri. Crucialmente, o SOC foi incluído desde o início nos estados de Bloch, garantindo que o modelo interpolado herde o emaranhamento spin-orbital ab initio em vez de adicionar o SOC como um termo perturbativo. O conjunto de base incluiu 48 funções spinor (projetadas sobre orbitais Fe 3d, X p e Te p) abrangendo o manifesto $p$–$d$ relevante.
• Formalismo de Torkance: O torque de spin–órbita (SOT) foi computado usando o formalismo de resposta linear de Kubo–Bastin. O tensor de torkance ($t_{ab}$), definido como torque por unidade de campo elétrico, foi decomposto em componentes pares e ímpares em relação à reversão de tempo (TR) (mar de Fermi e superfície de Fermi).
• Rotação da Magnetização: Para mapear a completa dependência angular do torque na esfera de magnetização, um esquema de "rotação de troca rígida" foi empregado. O Hamiltoniano ab initio foi computado uma vez para a magnetização ao longo de $\hat{z}$, então dividido em partes TR-pares e TR-ímpares. O campo de troca TR-ímpar foi rotacionado via transformações de spinor SU(2) para simular direções arbitrárias de magnetização ($\hat{m}$), enquanto a parte TR-par (termos orbitais e de SOC) permaneceu fixa. Esta abordagem foi testada contra a DFT não colinear totalmente autoconsistente, mostrando concordância dentro de $\sim 5\%$ para a maioria das orientações.
• Decomposição: Os torques resultantes foram analisados via expansão de harmônicos esféricos vetoriais, decomposição resolvida por átomo (sublatices Fe$_I$, Fe$_{II}$, Te, X) e mapeamento de densidade no espaço de momento.

Resultos Principais

1. Diferenças na Estrutura Eletrônica: Apesar de funções de trabalho semelhantes ($\sim 28$ meV de diferença), a dopagem por lacunas no FGaT altera significativamente o preenchimento de banda próximo aos pontos $K$ e $K'$. A densidade de estados em $\varepsilon_F$ no FGaT é reduzida por um fator de três em comparação ao FGT, e a polarização de spin no nível de Fermi inverte de majoritária (FGT) para minoritária (FGaT). O desdobramento (splitting) de valência-Zeeman permanece robusto em $\sim 40$ meV em ambos os materiais, ditado pelo ambiente de Te.
2. Dependência Angular dos Torques:
   • Torques TR-Ímpares (Superfície de Fermi): Estes dominam a resposta total em ambos os materiais. A dependência angular é notavelmente semelhante entre FGT e FGaT, com magnitudes de pico no FGaT apenas $\sim 15\%$ menores. Isso sugere que o canal TR-ímpar é amplamente transferível apesar das distintas configurações de superfície de Fermi.
   • Torques TR-Pares (Mar de Fermi): Estes exibem diferenças marcantes. O FGT apresenta uma contribuição de quarto harmônico pronunciada para magnetização no plano, ao passo que este componente é fortemente suprimido no FGaT. Os autores atribuem essa supressão à redução dos bolsos de Fermi nos pontos $K$ no FGaT, o que leva a uma compensação mais completa no espaço de momento das contribuições de harmônicos superiores.
3. Análise por Átomo e Torques Escondidos:
   • A análise revela que os dois átomos de Fe$_I$ (relacionados pela simetria de espelho horizontal $\sigma_h$) contribuem de forma diferente para o torque quando a magnetização se inclina da direção fora do plano.
   • Torques Escondidos: Embora a simetria global $D_{3h}$ proíba um acúmulo de spin induzido por corrente uniforme, a simetria de sublatice permite campos spin-órbita "escondidos". No FGaT, essas contribuições escondidas nas sublatices Fe$_I$ podem exceder os componentes permitidos pela simetria em quase um fator de 30. Embora se cancelem exatamente na monocamada ideal, o artigo as identifica como um recurso potencial para engenharia de torque se a simetria de sublatice for quebrada.
   • A densidade de torque TR-ímpar está confinada a contornos nítidos próximos aos pontos $K$, enquanto a densidade TR-par está espalhada pela zona, sendo intensificada em anticrossings de bandas.

Significância e Alegações
O artigo afirma fornecer uma compreensão microscópica dos torques induzidos por corrente em ferromagnetos 2D ao isolar os efeitos do preenchimento de banda da simetria cristalina. Ao comparar FGT e FGaT, os autores demonstram que:

• Estequiometria como um Controle (Knob): Mudanças qualitativas nos harmônicos de torque (especificamente a supressão do quarto harmônico no canal TR-par) podem ser alcançadas intrinsecamente através da modificação estequiométrica (Ge $\to$ Ga) em vez de via controle eletrostático externo.
• Graus de Liberdade Escondidos: O estudo destaca que monocamadas com quebra de inversão de simetria hospedam torques escondidos decorrentes de graus de liberdade específicos de sítio (sublatices) que não são capturados pelas visões convencionais de acúmulo de spin induzido por corrente. Esses torques escondidos, embora se cancelem na resposta líquida da camada ideal, sugerem que a engenharia de defeitos ou a quebra de simetria em heteroestruturas vdW poderia ativar componentes de torque grandes e, de outra forma, inacessíveis.
• Rigor Metodológico: O trabalho estabelece um pipeline computacional robusto para calcular SOTs em sistemas de baixa simetria, enfatizando a necessidade de tratamento de SOC autoconsistente e funções de Wannier adaptadas à simetria para capturar corretamente a interação entre campos de troca e acoplamento spin-órbita.

Os autores concluem que seus resultados oferecem orientação para a "engenharia de defeitos e de van der Waals de materiais magnéticos em camadas", sugerindo que manipular a simetria de sublatice pode ser uma rota viável para aumentar ou moldar os torques de spin-órbita em magnetos 2D.