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

Resumen Técnico: Estudio de Primeros Principios sobre el Apilamiento en Masa, la Imagen $J_{\text{eff}}$, el Hamiltoniano Magnético, los Factores $g$ y las Distorsiones Estructurales de $\alpha$-RuCl$_3$

Planteamiento del Problema
El $\alpha$-RuCl$_3$ es un candidato principal para la realización del estado de líquido de espín de Kitaev, aunque sus propiedades magnéticas fundamentales siguen siendo objeto de debate. Los problemas clave no resueltos incluyen:

1. Estructura en Masa: Aunque experimentos recientes sugieren que la fase de baja temperatura pertenece al grupo espacial $R\bar{3}$ (con una transición a $C2/m$ a $\sim$150 K), no existía un estudio teórico previo que confirmara la estructura $R\bar{3}$ como la fase de baja temperatura energéticamente favorable.
2. Estructura Electrónica: La validez de la imagen $J_{\text{eff}} = 1/2$ para el $\alpha$-RuCl$_3$ requiere verificación, particularmente respecto a la orientación del eje de cuantización del momento angular relativo al vector de Néel.
3. Hamiltoniano Magnético: Estudios previos reportan discrepancias significativas en los parámetros de intercambio magnético (por ejemplo, el signo del intercambio de Kitaev de primer vecino $K_1$) y a menudo omiten las interacciones de segundo vecino.
4. Anisotropía del Factor $g$: Los informes experimentales sobre la anisotropía del factor $g$ ($g_{XY}$ vs. $g_Z$) son contradictorios, variando desde anisotropías pequeñas hasta grandes. Además, los métodos teóricos para calcular estos factores (específicamente la proyección de orbitales atómicos) pueden ser insuficientes.
5. Distorsiones Estructurales: El impacto de distorsiones estructurales específicas en el octaedro de RuCl$_6$ (distorsión trigonal vs. giro relativo de los triángulos de Cl) sobre la anisotropía magnética no ha sido aislado sistemáticamente.

Metodología
Los autores emplearon la Teoría de la Funcional de la Densidad con Restricciones (cDFT) dentro del marco de Quantum ESPRESSO utilizando la Aproximación de Densidad Local (LDA) y métodos LDA+U rotacionalmente invariantes.

• Relajación Estructural: Los cálculos se realizaron tanto para las configuraciones de apilamiento $C2/m$ como $R\bar{3}$ bajo condiciones de relajación de celda fija, de celda variable y experimental para determinar la estabilidad energética.
• Análisis de $J_{\text{eff}}$: La estructura electrónica se analizó proyectando las funciones de onda sobre estados de base $J_{\text{eff}}$. Crucialmente, el eje de cuantización del momento angular se alineó paralelamente al vector de Néel para probar la imagen de $J_{\text{eff}}$.
• Construcción del Hamiltoniano Magnético: Los autores ajustaron las energías totales de diversas configuraciones magnéticas (zigzag, Néel, stripy, ferromagnética y una configuración específica de 1 por 4) a tres Hamiltonianos de modelo:
  1. Modelo convencional de Kitaev-Heisenberg (KH).
  2. Modelo de hielo de espín (SI).
  3. Un nuevo modelo Anisotrópico $J_1$-$J_2$-$J_3$ (A-$J_1J_2J_3$) que incluye interacciones anisotrópicas de segundo y tercer vecino.
• Cálculo del Factor $g$: Se compararon dos métodos:
  1. Proyección de Orbitales Atómicos: Proyectando los momentos sobre el manifold de Hubbard.
  2. Inclinación del Momento con Interpolación de Wannier: Inclinando el momento $J_{\text{eff}}$ por $1^\circ$ y computando el cambio en el momento magnético total mediante interpolación de Wannier de equivalencia de traslación (WannierBerri).
• Análisis de Distorsión: Se aislaron los efectos de dos distorsiones específicas de un octaedro perfecto: la distorsión trigonal (cambiando el ángulo Ru-Cl-Z) y el giro relativo (rotando los triángulos de Cl superior e inferior en el plano XY).

Contribuciones Clave y Resultados

1. Confirmación del Apilamiento $R\bar{3}$:
   El estudio demuestra que la estructura $R\bar{3}$ es de menor energía que la estructura $C2/m$ por aproximadamente 5–6 meV/Ru a través de diferentes valores de Hubbard $U$ y tipos de relajación. Esto proporciona el primer apoyo teórico para la fase $R\bar{3}$ observada experimentalmente a baja temperatura.

2. Imagen Refinada de $J_{\text{eff}}$:
   Se muestra que los estados del mínimo de la banda de conducción son casi puramente de carácter $J_{\text{eff}} = 1/2$ y $m_{\text{eff}} = -1/2$, pero solo cuando el eje de cuantización del momento angular se establece paralelo al vector de Néel. Esta perspectiva, previamente omitida, es esencial para identificar correctamente el carácter orbital de los estados. La magnitud del momento $J_{\text{eff}}$ resultó ser robusta a través de las configuraciones, aunque ligeramente menor al $1/2$ ideal debido a la hibridación y las distorsiones.

3. Hamiltoniano Magnético y Parámetros de Intercambio:

   • Los modelos convencionales KH y de hielo de espín fallan en reproducir con precisión las energías totales de DFT, particularmente para la configuración Néel.
   • El modelo propuesto A-$J_1J_2J_3$ (y una variante simplificada A-$J_1J_2$-I-$J_3$) reproduce con precisión las energías de DFT.
   • Las interacciones de segundo vecino son críticas. Los componentes anisotrópicos del intercambio de segundo vecino ($K_2, \Gamma_2, \Gamma'_2$) son no despreciables y, en magnitud, mayores que los términos de intercambio anisotrópicos de primer vecino ($\Gamma'_1$) frecuentemente considerados en modelos previos. La inclusión de estos términos redujo el error cuadrático medio (RMS) en la energía total en un factor de 3.6 en comparación con el modelo convencional.

4. Anisotropía del Factor $g$ y Método de Cálculo:

   • El método de inclinación del momento combinado con la interpolación de Wannier produce factores $g$ de $g_{XY} \approx 2.22$ y $g_Z \approx 1.84$. Estos valores indican una anisotropía pequeña, alineándose con reportes experimentales recientes (por ejemplo, Loidl et al.) y contradiciendo afirmaciones anteriores de una gran anisotropía.
   • El método de proyección de orbitales atómicos subestima significativamente la magnitud del factor $g$ ($g_{XY} \approx 1.54$), perdiendo una parte no despreciable del momento magnético total. Los autores concluyen que el cálculo preciso del magnetismo orbital requiere el enfoque de interpolación de Wannier.

5. Papel de las Distorsiones Estructurales:
   Al aislar las distorsiones, el estudio encuentra que el giro relativo de los triángulos de Cl (dentro del plano XY) es el factor dominante que determina la anisotropía magnética y la dirección correcta del momento del estado fundamental ($\theta \sim 60^\circ$). La distorsión trigonal por sí sola no logra reproducir el estado fundamental experimental. Este giro de la distorsión, presente en el material en masa sin perturbación externa, ha sido pasado por alto en modelos teóricos previos.

Significancia
El artículo afirma proporcionar una comprensión de primeros principios exhaustiva de $\alpha$-RuCl$_3$ al resolver contradicciones de larga data en sus descripciones estructurales y magnéticas. Específicamente:

• Valida teóricamente la estructura de baja temperatura $R\bar{3}$.
• Establece la necesidad de alinear el eje de cuantización con el vector de Néel para interpretar correctamente la imagen de $J_{\text{eff}}$.
• Demuestra que el modelado magnético preciso requiere la inclusión de interacciones de intercambio anisotrópicas de segundo vecino.
• Corrige la metodología para calcular los factores $g$, mostrando que la proyección de orbitales atómicos es insuficiente y que la verdadera anisotropía es pequeña.
• Identifica el giro relativo de los triángulos de Cl como un parámetro estructural crítico, previamente omitido, que gobierna la anisotropía magnética, sugiriendo una nueva vía para la ingeniería de la magnetismo en materiales de Kitaev.

Los autores concluyen que su trabajo ofrece una perspectiva diferente sobre el magnetismo del $\alpha$-RuCl$_3$ y sugiere que los futuros modelos teóricos y experimentos de alta presión deben incorporar estas distorsiones estructurales y términos de intercambio específicos para una interpretación precisa.