Quantitative Disentanglement of Terahertz Spin and Orbital Pumping in 3d Ferromagnetic Heterostructures

Resumen Técnico: Desentrelazamiento Cuantitativo del Bombeo de Espín y Orbital de Terahertz en Heteroestructuras Ferromagnéticas de 3d

Planteamiento del Problema
El bombeo de espín y orbital —la inyección de momento angular de espín y orbital desde un ferromagneto (FM) excitado hacia una capa no magnética (NM) adyacente— son fundamentales para el transporte de momento angular en heteroestructuras magnéticas. Si bien la excitación óptica de femtosegundos ha extendido estos fenómenos al régimen ultrafast, detectables mediante la espectroscopia de emisión de terahertz de espín (STES), una limitación crítica persiste: las corrientes de espín y orbital se generan simultáneamente y contribuyen de forma aditiva a la misma señal de emisión de terahertz. Esta superposición impide su separación cuantitativa. Aunque la teoría microscópica predice que la eficiencia del bombeo orbital aumenta de Fe a Ni a través de la serie 3d, no está claro si estas predicciones se mantienen bajo excitación ultrafast o cómo se comparan cuantitativamente las contribuciones relativas de las corrientes de espín y orbital. Los intentos previos de aislar los canales orbitales utilizando convertidores de espín a carga débiles (p. ej., Cu, Ti) no lograron determinar las amplitudes relativas, polaridades o longitudes de propagación dentro de un único marco de trabajo.

Metodología
Para superar el desafío de la generación simultánea y la detección aditiva, los autores desarrollaron una estrategia de desentrelazamiento cuantitativo que combina STES con un control preciso del espesor de muestras en cuña. El diseño experimental se basa en dos propiedades complementarias de las corrientes de espín y orbital:

1. Polaridades Opuestas: Los investigadores seleccionaron capas NM (Ta, W y Nb) donde los ángulos de Hall de espín ($\gamma_{SH}^{NM}$) y de orbital ($\gamma_{OH}^{NM}$) poseen signos opuestos. En consecuencia, las emisiones de terahertz inducidas por la corriente de espín y las inducidas por la corriente de orbital exhiben polaridades opuestas.
2. Longitudes de Decaimiento Distintas: Las corrientes de espín y orbital decaen sobre diferentes longitudes características dentro de la capa NM. Mediante el uso de muestras en forma de cuña con gradientes de espesor continuos (0.75–2.5 nm/mm), los autores realizaron mediciones de alta precisión dependientes del espesor para separar las contribuciones basadas en sus perfiles de decaimiento espacial.

El estudio investigó sistemáticamente heteroestructuras de Fe, Co y Ni (fijadas en 5 nm de espesor) recubiertas con capas en cuña de Ta, W o Nb. Las formas de onda de terahertz se registraron mediante muestreo electro-óptico de espacio libre. Para extraer cuantitativamente las contribuciones individuales, los autores emplearon un Modelo de Emisión de Terahertz de Múltiples Componentes (MC-TEM). Este modelo trata la emisión total como una suma coherente de tres componentes: la emisión del FM puro, la conversión de espín a carga mediada por el efecto Hall de espín inverso (ISHE) y la conversión de orbital a carga mediada por el efecto Hall de orbital inverso (IOHE). Los efectos de Rashba-Edelstein interfaciales se excluyeron mediante experimentos de control.

Resultos Clave

• Evolución de la Polaridad y Comportamiento de Umbral:

  • Fe|Ta: Exhibió una polaridad negativa a lo largo de todo el rango de espesor, pero mostró un régimen distintivo de "tipo umbral" por debajo de $d_{Ta} \approx 0.5$ nm, donde el crecimiento de la amplitud se vio suprimido. Esto se atribuyó a la casi cancelación de las contribuciones negativas de espín y positivas de orbital.
  • Ni|Ta: Mostró una reversión de polaridad única de negativa a positiva en $d_{Ta} \approx 0.3$ nm, indicando dominancia orbital en espesores delgados.
  • Co|Ta: Mostró una doble reversión de polaridad (negativa $\to$ positiva $\to$ negativa), revelando una competencia compleja donde las contribuciones orbitales dominan en capas delgadas y las contribuciones de espín dominan en capas más gruesas.
  • FM|W: Fe|W y Co|W mantuvieron una polaridad negativa en todo el rango de espesor, indicando dominancia de espín. Ni|W exhibió dos reversiones de polaridad: una primera reversión en $d_W \approx 0.2$ nm (dominancia orbital) y una segunda en $d_W \approx 2.6$ nm (reemergencia de la dominancia de espín).
  • FM|Nb: Fe|Nb mantuvo una polaridad negativa. Ni|Nb exhibió una reversión de polaridad única de negativa a positiva cerca de $d_{Nb} \approx 0.6$ nm. Co|Nb mostró dos reversiones de polaridad (negativa $\to$ positiva $\to$ negativa), ocurriendo la segunda en $d_{Nb} \approx 5.4$ nm.
  • Se observaron regímenes de tipo umbral (crecimiento de amplitud suprimido en capas NM delgadas) en Fe|Ta, Co|W y Fe|Nb, derivados de la casi cancelación de las señales de espín y orbital.

• Descomposición Cuantitativa:

  • El ajuste MC-TEM reprodujo con éxito las evoluciones de amplitud y forma de onda para todos los sistemas.
  • Las razones de amplitud de corriente orbital a espín ($\tilde{A}_L / \tilde{A}_S$) extraídas aumentaron sistemáticamente de Fe a Co a Ni.
  • Para Fe y Co, la contribución de la corriente orbital osciló entre unos pocos puntos porcentuales y aproximadamente el 20% de la contribución de la corriente de espín.
  • Para Ni, la contribución orbital alcanzó varias decenas de puntos porcentuales, con estimaciones de hasta ~70% en Ni|Nb, superando significativamente las predicciones teóricas actuales para la serie 3d.
  • Incluso el Fe, tradicionalmente visto como una fuente dominante de espín, resultó generar una corriente orbital no despreciable, la cual es esencial para interpretar el comportamiento de umbral en regímenes de NM delgada.

• Longitudes de Difusión:

  • Las longitudes de difusión orbital ($\lambda_L$) extraídas fueron consistentemente más cortas que las longitudes de difusión de espín ($\lambda_S$) correspondientes en todos los materiales NM.
  • Notablemente, $\lambda_L$ fue marcadamente más larga en Nb (~1.5 nm) en comparación con los metales 5d más pesados Ta y W (0.3–0.8 nm), lo que sugiere que la propagación orbital depende de la fuerza del acoplamiento espín-órbita (SOC) de la capa NM.

Significancia y Reivindicaciones
El artículo establece un marco cuantitativo para investigar el bombeo de espín y orbital ultrafast en heteroestructuras magnéticas. Al resolver las contribuciones coexistentes de las corrientes de espín y orbital dentro de una única plataforma STES, este trabajo proporciona una base para evaluar las teorías microscópicas de la generación y el transporte de momento angular orbital (OAM).

Los autores destacan que, si bien sus resultados confirman el aumento progresivo del bombeo orbital de Fe a Ni, la magnitud de la contribución orbital es sustancialmente mayor de lo predicho por las teorías existentes. Atribuyen esta discrepancia a posibles diferencias entre la dinámica coherente y adiabática (base teórica) y la dinámica fuera de equilibrio ultrafast (realidad experimental), o a variaciones en las eficiencias de transmisión interfacial. Además, la observación de que las longitudes de difusión orbital varían con el SOC de la capa NM desafía las simulaciones recientes que sugieren la independencia del SOC.

El estudio concluye que la STES, combinada con el análisis de muestras en cuña y el modelado de múltiples componentes, sirve como una sonda cuantitativa robusta para el acoplamiento de los bombeos de espín y orbital. Este enfoque abre nuevas vías para la ingeniería de la interconversión espín-orbital ultrafast en materiales y dispositivos orbitrónicos, aunque los autores mantienen la modestia respecto a la cuantificación precisa de las densidades de corriente debido a las incertidumbres en los valores teóricos de las conductividades de Hall de espín y orbital.