First-principles study of bulk stacking, $J_{\rm eff}$ picture, magnetic Hamiltonian, $g$ factors, and structural distortions of $α$-RuCl$_3$

Resumo Técnico: Estudo de Primeiros Princípios sobre Empilhamento de Volume, Imagem $J_{\text{eff}}$, Hamiltoniana Magnética, Fatores $g$ e Distorções Estruturais da $\alpha$-RuCl$_3$

Definição do Problema
A $\alpha$-RuCl$_3$ é uma candidata primária para a realização do estado de líquido de spin de Kitaev, embora suas propriedades magnéticas fundamentais permaneçam objeto de debate. As principais questões não resolvidas incluem:

1. Estrutura de Volume: Embora estudos recentes sugiram que a fase de baixa temperatura pertence ao grupo de espaço $R\bar{3}$ (com uma transição para $C2/m$ a $\sim$150 K), não havia estudo teórico prévio confirmando a estrutura $R\bar{3}$ como a fase energeticamente favorável de baixa temperatura.
2. Estrutura Eletrônica: A validade da imagem $J_{\text{eff}} = 1/2$ para a $\alpha$-RuCl$_3$ requer verificação, particularmente em relação à orientação do eixo de quantização do momento angular em relação ao vetor Néel.
3. Hamiltoniana Magnética: Estudos anteriores relatam discrepâncias significativas nos parâmetros de troca magnética (por exemplo, o sinal da troca de Kitaev de vizinho mais próximo $K_1$) e frequentemente negligenciam interações de segundo vizinho mais próximo.
4. Anisotropia do Fator $g$: Relatos experimentais sobre a anisotropia do fator $g$ ($g_{XY}$ vs. $g_Z$) são contraditórios, variando de anisotropias pequenas a grandes. Além disso, os métodos teóricos para calcular esses fatores (especificamente a projeção de orbitais atômicos) podem ser insuficientes.
5. Distorções Estruturais: O impacto de distorções estruturais específicas no octaedro de RuCl$_6$ (distorção trigonal vs. torção relativa dos triângulos de Cl) na anisotropia magnética não foi sistematicamente isolado.

Metodologia
Os autores empregaram DFT restrita (cDFT) dentro do framework Quantum ESPRESSO usando a Aproximação de Densidade Local (LDA) e métodos LDA+U rotacionalmente invariantes.

• Relaxação Estrutural: Cálculos foram realizados para ambas as configurações de empilhamento $C2/m$ e $R\bar{3}$ sob condições experimentais, de célula fixa e de relaxação de célula variável para determinar a estabilidade energética.
• Análise $J_{\text{eff}}$: A estrutura eletrônica foi analisada projetando as funções de onda sobre estados de base $J_{\text{eff}}$. Crucialmente, o eixo de quantização do momento angular foi alinhado paralelamente ao vetor Néel para testar a imagem $J_{\text{eff}}$.
• Construção da Hamiltoniana Magnética: Os autores ajustaram as energias totais de várias configurações magnéticas (zigzag, Néel, stripy, ferromagnética e uma configuração específica de 1 por 4) a três Hamiltonianas de modelo:
  1. Modelo convencional de Kitaev-Heisenberg (KH).
  2. Modelo de gelo de spin (SI).
  3. Um novo modelo Anisotrópico $J_1$-$J_2$-$J_3$ (A-$J_1J_2J_3$) incluindo interações anisotrópicas de segundo e terceiro vizinho mais próximo.
• Cálculo do Fator $g$: Dois métodos foram comparados:
  1. Projeção de Orbital Atômico: Projetando momentos sobre o manifold de Hubbard.
  2. Inclinação de Momento com Interpolação de Wannier: Inclinando o momento $J_{\text{eff}}$ em $1^\circ$ e computando a mudança no momento magnético total usando interpolação de Wannier transladacionalmente equivariante (WannierBerri).
• Análise de Distorção: Os efeitos de duas distorções específicas de um octaedro perfeito foram isolados: distorção trigonal (mudando o ângulo Ru-Cl-Z) e torção relativa (rotacionando os triângulos de Cl superiores e inferiores no plano XY).

Principais Contribuições e Resultados

1. Confirmação do Empilhamento $R\bar{3}$:
   O estudo demonstra que a estrutura $R\bar{3}$ é de menor energia do que a estrutura $C2/m$ por aproximadamente 5–6 meV/Ru através de diferentes valores de Hubbard $U$ e tipos de relaxação. Isso fornece o primeiro suporte teórico para a fase $R\bar{3}$ observada experimentalmente em baixa temperatura.

2. Imagem $J_{\text{eff}}$ Refinada:
   Mostra-se que os estados do mínimo da banda de condução são quase puramente de caráter $J_{\text{eff}} = 1/2$ e $m_{\text{eff}} = -1/2$, mas apenas quando o eixo de quantização do momento angular é definido paralelo ao vetor Néel. Esta perspectiva, anteriormente negligenciada, é essencial para identificar corretamente o caráter orbital dos estados. A magnitude do momento $J_{\text{eff}}$ mostrou-se robusta entre as configurações, embora ligeiramente menor que o ideal de $1/2$ devido à hibridização e distorções.

3. Hamiltoniana Magnética e Parâmetros de Troca:

   • Os modelos convencionais de KH e de gelo de spin falham em reproduzir com precisão as energias totais de DFT, particularmente para a configuração Néel.
   • O proposto modelo A-$J_1J_2J_3$ (e uma variante simplificada A-$J_1J_2$-I-$J_3$) reproduz com precisão as energias de DFT.
   • Interações de segundo vizinho mais próximo são críticas. Os componentes anisotrópicos da troca de segundo vizinho mais próximo ($K_2, \Gamma_2, \Gamma'_2$) são não desprezíveis e, em magnitude, maiores do que os termos de troca de vizinho mais próximo anisotrópicos ($\Gamma'_1$) frequentemente considerados em modelos anteriores. A inclusão desses termos reduziu o erro quadrático médio (RMS) na energia total em um fator de 3,6 em comparação com o modelo convencional.

4. Anisotropia do Fator $g$ e Método de Cálculo:

   • O método de inclinação de momento combinado com interpolação de Wannier produz fatores $g$ de $g_{XY} \approx 2,22$ e $g_Z \approx 1,84$. Estes valores indicam uma pequena anisotropia, alinhando-se com relatórios experimentais recentes (ex: Loidl et al.) e contradizendo alegações anteriores de grande anisotropia.
   • O método de projeção de orbital atômico subestima significativamente a magnitude do fator $g$ ($g_{XY} \approx 1,54$), perdendo uma parte não desprezível do momento magnético total. Os autores concluem que o cálculo preciso do magnetismo orbital requer a abordagem de interpolação de Wannier.

5. Papel das Distorções Estruturais:
   Ao isolar as distorções, o estudo descobre que a torção relativa dos triângulos de Cl (dentro do plano XY) é o fator dominante que determina a anisotropia magnética e a direção correta do momento do estado fundamental ($\theta \sim 60^\circ$). A distorção trigonal sozinha falha em reproduzir o estado fundamental experimental. Esta torção estrutural, presente no material bulk sem perturbação externa, foi negligenciada em modelos teóricos anteriores.

Significância
O artigo afirma fornecer uma compreensão de primeiros princípios abrangente da $\alpha$-RuCl$_3$ ao resolver contradições de longa data em suas descrições estruturais e magnéticas. Especificamente:

• Valida a estrutura de baixa temperatura $R\bar{3}$ teoricamente.
• Estabelece a necessidade de alinhar o eixo de quantização com o vetor Néel para interpretar corretamente a imagem $J_{\text{eff}}$.
• Demonstra que o modelagem magnética precisa requer a inclusão de interações de troca de segundo vizinho mais próximo anisotrópicas.
• Corrige a metodologia para o cálculo de fatores $g$, mostrando que a projeção de orbital atômico é insuficiente e que a verdadeira anisotropia é pequena.
• Identifica a torção relativa dos triângulos de Cl como um parâmetro estrutural crítico e anteriormente negligenciado que governa a anisotropia magnética, sugerindo uma nova via para a engenharia de magnetismo em materiais de Kitaev.

Os autores concluem que seu trabalho oferece uma visão diferente sobre o magnetismo da $\alpha$-RuCl$_3$ e sugere que futuros modelos teóricos e experimentos de alta pressão devem incorporar essas distorções estruturais e termos de troca específicos para uma interpretação precisa.