Platform for zero-field isolated skyrmions: 4$d$/Co atomic bilayers on Re(0001)

Resumen Técnico: Plataforma para Escarmiones Aislados de Campo Cero en Bilas 4d/Co sobre Re(0001)

Planteamiento del Problema
Los escarmiones magnéticos, texturas de espín con protección topológica, son candidatos prometedores para dispositivos espintrónicos de próxima generación, computación neuromórfica y tecnologías cuánticas debido a su tamaño a escala nanométrica y su estabilidad. Aunque se han realizado escarmiones aislados experimentalmente en diversas películas ultradelgadas y multicapas de metales de transición, estos suelen requerir un campo magnético externo para su estabilización. La estabilización sin campo es altamente deseable para la integración de dispositivos con el fin de eliminar los componentes generadores de campo y reducir el consumo de energía. Hasta la fecha, dentro de la familia de películas ultradelgadas de metales de transición, solo se han observado los sistemas de la bilayer atómica Rh/Co sobre Ir(111) capaces de albergar escarmiones aislados de menos de 10 nm a campo cero, estabilizados por frustración de intercambio. Predicciones teóricas sugieren que las bilayers 4d/Co sobre la superficie de Re(0001) también podrían albergar tales estados, pero estudios previos dependieron de modelos de espín efectivos que omitieron interacciones de espín de orden superior (HOI, por sus siglas en inglés) y carecieron de simulaciones de espín atomísticas para confirmar la estabilidad.

Metodología
Este estudio emplea un enfoque combinado de teoría del funcional de la densidad (DFT) de primeros principios y simulaciones de espín atomísticas para investigar tres sistemas específicos: bilayers atómicas de fcc-Rh/Co, hcp-Pd/Co y fcc-Ru/Co crecidas sobre la superficie de Re(0001).

1. Cálculos de DFT: Las estructuras electrónica y magnética se calcularon utilizando el código fleur basado en el método de ondas planas aumentadas lineales de potencial completo (FP-LAPW). El estudio modeló las bilayers 4d/Co sobre nueve capas de Re(0001).

   • Dispersión de Energía: Se calculó la dispersión de energía de las espirales de espín homogéneas en función del vector de propagación $q$ a lo largo de direcciones de alta simetría ($\Gamma K M$ y $\Gamma M$), tanto con como sin acoplamiento espín-órbita (SOC).
   • Extracción de Interacciones: Las constantes de intercambio de pares ($J_{ij}$) y de interacción Dzyaloshinskii-Moriya (DMI, $D_{ij}$) se extrajeron mapeando las energías de DFT a Hamiltonianos de Heisenberg y DMI.
   • Interacciones de Orden Superior (HOI): Para ir más allá de los modelos estándar, el estudio calculó tres constantes de HOI: bicuadrática ($B_1$), de cuatro espines de tres sitios ($Y_1$) y de cuatro espines de cuatro sitios ($K_1$). Estas se derivaron de las diferencias de energía entre estados multi-Q específicos (uudd y 3Q) y los estados correspondientes de espiral de espín de un solo Q.
   • Anisotropía Magnetocristalina (MAE): Se calculó de forma autoconsistente para determinar el eje de magnetización fácil.

2. Simulaciones de Espín Atomístico: Se construyó un modelo de espín atomístico extendido incorporando todos los parámetros derivados de DFT (intercambio de Heisenberg, DMI, MAE y HOI).

   • Relajación: Las estructuras de espín se relajaron utilizando la ecuación de Landau-Lifshitz en una red hexagonal de $150 \times 150$ para determinar el diagrama de fases magnéticas a temperatura cero.
   • Análisis de Estabilidad: Se utilizó el método de Banda Elástica de Nudged Geodésica (GNEB) para calcular la ruta de mínima energía (MEP) para el colapso de los escarmiones aislados hacia el estado fundamental ferromagnético (FM), determinando las barreras de energía de aniquilación y creación.

Contribuciones Clave y Resultados

• Identificación de Plataformas de Escarmiones de Campo Cero: Las simulaciones del modelo de espín extendido revelan que los escarmiones aislados emergen espontáneamente sobre el fondo ferromagnético a campo cero para fcc-Rh/Co/Re(0001) y hcp-Pd/Co/Re(0001). En contraste, no se pudieron estabilizar escarmiones aislados en fcc-Ru/Co/Re(0001) ni a campo cero ni a campos finitos debido a una fuerza de DMI insuficiente.
• Dimensiones de los Escarmiones: El radio de los escarmiones aislados de campo cero es aproximadamente de 6 nm para fcc-Rh/Co/Re(0001) y de 12 nm para hcp-Pd/Co/Re(0001). Estos radios disminuyen a medida que aumenta el campo magnético externo.
• Papel de las Interacciones de Orden Superior: El estudio destaca que, más allá de las interacciones de intercambio de vecinos más cercanos y DMI, se induce frustración en estas películas. Crucialmente, la inclusión de HOI (específicamente la interacción de cuatro espines de tres sitios) es necesaria para una descripción precisa del sistema, aunque el modelo de espín efectivo (que omite las HOI) predijo correctamente el estado fundamental.
• Barreras de Energía y Estabilidad: Los cálculos de la teoría del estado de transición muestran que los escarmiones aislados en los sistemas estables están protegidos por barreras de energía sustanciales de aproximadamente 150 meV (barrera de aniquilación) a campo cero.
  • Composición de la Barrera: El DMI proporciona la contribución positiva dominante a la barrera de energía, con una contribución positiva menor de la interacción de cuatro espines de tres sitios. Los términos de intercambio de pares y MAE actúan para reducir la barrera (favoreciendo el estado FM). La interacción de cuatro espines de cuatro sitios proporciona una contribución insignificante.
  • Observabilidad: La magnitud de estas barreras sugiere que estos escarmiones podrían observarse en experimentos de baja temperatura.
• Propiedades de los Materiales:
  • Todas las películas exhiben un eje de magnetización fácil fuera del plano.
  • El DMI en fcc-Rh/Co y hcp-Pd/Co es impulsado principalmente por el sustrato de Re, mientras que en fcc-Ru/Co, las capas de Ru y Co contribuyen de manera diferente, resultando en un DMI total significativamente menor.
  • Los momentos magnéticos de las capas de Co y 4d muestran una variación significativa con los vectores de onda de la espiral de espín, una característica que no es capturada por modelos de Heisenberg simples con momentos fijos, aunque los momentos permanecen bastante constantes cerca del estado FM.

Significancia y Reivindicaciones
Los autores proponen las bilayers 4d/Co sobre la superficie de Re(0001) como una nueva y prometedora plataforma para la realización de escarmiones aislados a escala nanométrica y de campo cero. El trabajo demuestra que combinaciones específicas de metales de transición 4d (Rh, Pd) con Co sobre Re pueden estabilizar estas texturas topológicas sin campos externos, impulsadas por una compleja interacción de frustración de intercambio, DMI e interacciones de espín de orden superior.

El artículo afirma que las barreras de energía calculadas son lo suficientemente altas como para permitir la observación experimental mediante microscopía de efecto túnel de espín polarizado (SP-STM) a bajas temperaturas. Además, los autores señalan que, dado que el sustrato de Re se vuelve superconductor por debajo de 1.7 K, estos sistemas podrían servir como sistemas híbridos magnet-superconductor de modelo, potencialmente relevantes para la superconductividad topológica emergente. El estudio enfatiza que, si bien el modelo de espín efectivo predijo la existencia de estos estados, el modelo extendido que incluye HOI es esencial para evaluar con precisión la estabilidad y las propiedades de los escarmiones.