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

Resumen Técnico: Anatomía de los Torques de Espín-Órbita en Monocapas de Fe$_3$GeTe$_2$ y Fe$_3$GaTe$_2$

Planteamiento del Problema
El control eléctrico de la magnetización mediante torques de espín-órbita (SOT) es un objetivo central en la espintrónica. Si bien los imanes de van der Waals (vdW) como el Fe$_3$GeTe$_2$ (FGT) y el Fe$_3$GaTe$_2$ (FGaT) ofrecen oportunidades únicas debido a su baja simetría cristalina y densidad de portadores sintonizable mediante puerta, los orígenes microscópicos de sus torques inducidos por corriente no se comprenden completamente. Específicamente, la dependencia de los armónicos de torque respecto al llenado de bandas, el papel de los estados electrónicos específicos cerca del nivel de Fermi ($\varepsilon_F$) y la existencia de torques "ocultos" derivados de los grados de libertad de la subred en monocapas no centrosimétricas requieren una investigación sistemática. Aunque el FGT y el FGaT comparten la misma estructura cristalina $D_{3h}$ y un carácter dominante de Fe 3d, difieren por un electrón de valencia (Ga frente a Ge), lo que convierte efectivamente al FGaT en un análogo dopado con huecos del FGT. El artículo aborda cómo esta modificación estequiométrica altera el panorama de los SOT sin la necesidad de un control electrostático por puerta.

Metodología
Los autores emplean un enfoque sistemático de primeros principios que combina la teoría del funcional de la densidad (DFT) con la teoría de respuesta lineal:

• Estructura Electrónica: Los cálculos se realizaron utilizando Quantum ESPRESSO con pseudopotenciales de Vanderbilt (ONCV) optimizados y totalmente relativistas. El acoplamiento espín-órbita (SOC) se trató de forma autoconsistente dentro de un marco no colineal. Se utilizó el funcional Perdew–Wang LDA para los cálculos de respuesta para asegurar momentos magnéticos precisos, consistentes con estudios previos sobre la monocapa de FGT.
• Interpolación de Wannier: Para permitir un muestreo denso de la zona de Brillouin, se construyeron funciones de Wannier adaptadas a la simetría (SAWF) utilizando WannierBerri. Crucialmente, el SOC se incluyó desde el principio en los estados de Bloch, asegurando que el modelo interpolado herede el entrelazamiento espín-orbital ab initio en lugar de añadir el SOC como un término perturbativo. El conjunto de base incluyó 48 funciones espinoras (proyectadas sobre orbitales Fe 3d, X p y Te p) que abarcan el manifold $p$–$d$ relevante.
• Formalismo de Torkance: El torque de espín-órbita (SOT) se calculó utilizando el formalismo de respuesta lineal de Kubo–Bastin. El tensor de torkance ($t_{ab}$), definido como el torque por unidad de campo eléctrico, se descompuso en componentes pares (mar de Fermi) e impares (superficie de Fermi) respecto a la inversión temporal (TR).
• Rotación de la Magnetización: Para mapear la dependencia angular completa del torque en la esfera de magnetización, se empleó un esquema de "rotación de intercambio rígido". Se calculó una vez el Hamiltoniano ab initio para la magnetización a lo largo de $\hat{z}$, y luego se dividió en partes TR-pares y TR-impares. El campo de intercambio TR-impar se rotó mediante transformaciones de espínores SU(2) para simular direcciones de magnetización arbitrarias ($\hat{m}$), mientras que la parte TR-par (términos orbitales y de SOC) permaneció fija. Este enfoque fue validado frente a la DFT no colineal totalmente autoconsistente, mostrando una concordancia del $\sim 5\%$ para la mayoría de las orientaciones.
• Descomposición: Los torques resultantes se analizaron mediante expansión de armónicos esféricos vectoriales, descomposición resuelta por átomo (subredes Fe$_I$, Fe$_{II}$, Te, X) y mapeo de densidad en el espacio de momentos.

Resultos Clave

1. Diferencias en la Estructura Electrónica: A pesar de tener funciones de trabajo similares ($\sim 28$ meV de diferencia), el dopaje con huecos en el FGaT altera significamente el llenado de bandas cerca de los puntos $K$ y $K'$. La densidad de estados en $\varepsilon_F$ en el FGaT se reduce en un factor de tres en comparación con el FGT, y la polarización de espín en el nivel de Fermi se invierte de mayoritaria (FGT) a minoritaria (FGaT). El desdoblamiento de valle-Zeeman permanece robusto en $\sim 40$ meV en ambos materiales, dictado por el entorno de Te.
2. Dependencia Angular de los Torques:
   • Torques TR-impares (Superficie de Fermi): Estos dominan la respuesta total en ambos materiales. La dependencia angular es notablemente similar entre el FGT y el FGaT, con magnitudes máximas en el FGaT solo un $\sim 15\%$ menores. Esto sugiere que el canal TR-impar es ampliamente transferible a pesar de las distintas configuraciones de la superficie de Fermi.
   • Torques TR-pares (Mar de Fermi): Exhiben diferencias marcadas. El FGT muestra una contribución pronunciada del cuarto armónico para la magnetización en el plano, mientras que este componente se ve fuertemente suprimido en el FGaT. Los autores atribuyen esta supresión a la reducción de las bolsas de Fermi en el punto $K$ en el FGaT, lo que conduce a una compensación más completa en el espacio de momentos de las contribuciones de altos armónicos.
3. Torques Resueltos por Átomo y Torques Ocultos:
   • El análisis revela que los dos átomos de Fe$_I$ (relacionados por el espejo horizontal $\sigma_h$) contribuyen de manera diferente al torque cuando la magnetización se inclina desde la dirección fuera del plano.
   • Torques Ocultos: Mientras que la simetría $D_{3h}$ global prohíbe una acumulación de espín inducida por corriente uniforme, la simetría de la subred permite la existencia de campos espín-órbita "ocultos". En el FGaT, estas contribuciones ocultas en las subredes de Fe$_I$ pueden exceder los componentes permitidos por la simetría por casi un factor de 30. Aunque se cancelan exactamente en la monocapa ideal, el artículo los identifica como un recurso potencial para la ingeniería de torque si se rompe la simetría de la subred.
   • La densidad de torque TR-impar se confina a contornos agudos cerca de los puntos $K$, mientras que la densidad TR-par se distribuye por toda la zona, potenciada en las anticruzaciones de bandas.

Significancia y Reivindicaciones
El artículo afirma proporcionar una comprensión microscópica de los torques inducidos por corriente en ferromagnetos 2D al aislar los efectos del llenado de bandas de la simetría cristalina. Al comparar FGT y FGaT, los autores demuestran que:

• La Estequiometría como Palanca de Control: Los cambios cualitativos en los armónicos de torque (específicamente la supresión del cuarto armónico en el canal TR-par) pueden lograrse intrínsecamente mediante la modificación estequiométrica (Ge $\to$ Ga) en lugar de mediante el control electrostático externo.
• Grados de Libertad Ocultos: El estudio destaca que las monocapas con asimetría de inversión albergan torques ocultos derivados de grados de libertad específicos del sitio (subredes) que no son capturados por las imágenes convencionales de la acumulación de espín inducida por corriente. Estos torques ocultos, aunque se cancelan en la respuesta neta de la capa ideal, sugieren que la ingeniería de defectos o la ruptura de simetría en heteroestructuras de vdW podría activar componentes de torque grandes y de otro modo inaccesibles.
• Rigor Metodológico: El trabajo establece un flujo de trabajo computacional robusto para calcular SOTs en sistemas de baja simetría, enfatizando la necesidad de un tratamiento de SOC autoconsistente y funciones de Wannier adaptadas a la simetría para capturar correctamente la interacción entre los campos de intercambio y el acoplamiento espín-órbita.

Los autores concluyen que sus resultados ofrecen orientación para la "ingeniería de defectos y de van der Waals de materiales magnéticos en capas", sugiriendo que manipular la simetría de la subred podría ser una vía viable para mejorar o adaptar los torques de espín-órbita en imanes 2D.