The Critical Role of Itinerant Contributions to Orbital Angular Momentum Relaxation and Dynamics

Resumen Técnico: El Papel Crítico de las Contribuciones Itinerantes en la Relajación y Dinámica del Momento Angular Orbital

Planteamiento del Problema
El momento angular orbital (OAM) es reconocido como un grado de libertad prometedor para el transporte de baja disipación y el control de la magnetización. Sin embargo, existe una discrepancia significativa entre las mediciones experimentales y las predicciones teóricas con respecto a los mecanismos de relajación del OAM. Experimentalmente, las longitudes de difusión orbital varían ampliamente, desde ~70–80 nm en algunos sistemas hasta apenas 2–9 nm en otros. Teóricamente, las simulaciones de transporte cuántico de primeros principios basadas en la aproximación centrada en átomos (ACA) predicen longitudes de relajación orbital de menos de 1 nm, mientras que los cálculos de transporte no local en redes multi-orbitales sugieren respuestas que persisten durante cientos o miles de constantes de red.

Los autores postulan que este conflicto no surge de diferencias materiales o fallos metodológicos, sino de la definición subyacente del OAM. Los enfoques existentes basados en ACA tratan el OAM como una cantidad puramente local proyectada sobre orbitales atómicos, omitiendo las contribuciones itinerantes asociadas con la estructura geométrica de los estados de Bloch (contribuciones de la fase de Berry). Por el contrario, algunas simulaciones de transporte capturan implícitamente estos efectos itinerantes al inferir el transporte a partir del carácter orbital en lugar de una evaluación directa del OAM. El artículo pretende resolver esto desarrollando un marco que trate el OAM local e itinerante en igualdad de condiciones, incorporando simultáneamente la dispersión cuántica y la dinámica de relajación.

Metodología
El estudio emplea dos enfoques complementarios de primeros principios para desentrañar los canales de relajación de espín y orbital en MoS$_2$ monocapa (fuerte acoplamiento espín-órbita, SOC) y siliceno (SOC débil):

1. Dinámica de la Matriz de Densidad de Lindblad de Primeros Principios (FPDMD):

   • Este método calcula el OAM utilizando un formalismo de fase de Berry derivado de las derivadas en el espacio $k$ de los estados de Bloch, consistente con la teoría moderna de la polarización. Esta formulación incluye intrínsecamente las contribuciones itinerantes.
   • La dinámica está gobernada por la ecuación de movimiento de Lindblad para la matriz de densidad de un electrón, incluyendo explícitamente la dispersión electrón-fonón y el acoplamiento espín-órbita.
   • Los observables (espín y OAM) se calculan como la traza del operador con la matriz de densidad en evolución temporal, y las vidas medias se extraen ajustando el decaimiento a funciones exponenciales.

2. Transporte Cuántico de Red de Tight-Binding (TB) en Espacio Real:

   • Este enfoque utiliza un método de escala lineal basado en la expansión de polinomios de Chebyshev de las funciones espectrales.
   • Crucialmente, este método permite la separación explícita de los componentes del OAM: el componente centrado en el átomo ($L^{\text{ACA}}$) y el componente itinerante completo que incluye las contribuciones de la fase de Berry ($L^{\text{Berry}}$).
   • Los Hamiltonianos de TB están parametrizados mediante datos de primeros principios (por ejemplo, vía PAOFLOW para MoS$_2$ y estructuras de bandas de DFT para el siliceno). Se introduce el desorden (Anderson) para coincidir con las tasas de dispersión observadas en FPDMD.

Resultados Clave

• MoS$_2$ (Régimen de SOC Fuerte):

  • Dinámica de Múltiples Escalas Temporales: La relajación orbital exhibe dos escalas de tiempo distintas. Un decaimiento inicial rápido (0.25 ps) que corresponde a la redistribución de valle de las poblaciones polarizadas en el valle, donde los portadores se dispersan hacia valles opuestos con orbital opuesto pero con la misma polarización de espín. Sigue un decaimiento más lento (6.1 ps), donde los grados de libertad de espín y orbital están acoplados por el fuerte SOC.
  • Dominio Itinerante: Al comparar los componentes, la vida media orbital solo de ACA se calcula en < 3 fs, lo cual es tres a cuatro órdenes de magnitud más corta que el resultado completo de FPDMD. Esto indica que, en sistemas de fuerte SOC, la relajación está gobernada por el OAM itinerante, el cual la ACA omite.

• Siliceno (Régimen de SOC Débil):

  • Dinámicas Desacopladas: El espín y la dinámica orbital se desacoplan. Mientras que un campo eléctrico externo potencia drásticamente la relajación de espín (vía el desdoblamiento y la precesión inducidos por el campo), las vidas medias orbitales permanecen mayormente independientes del campo.
  • Comportamiento ACA vs. Itinerante: El componente ACA exhibe oscilaciones rápidas (frecuencia $\omega_{p_z} \approx 0.7 \times \Delta_{p_z}/\hbar$) impulsadas por el desdoblamiento del campo cristalino, decayendo en una escala de tiempo sub-picosegundo. En contraste, el OAM itinerante muestra una oscilación mínima y decae en una escala de tiempo de decenas de picosegundos.
  • Discrepancia de Vida Media: La vida media orbital de la ACA es al menos un orden de magnitud más corta que la vida media itinerante en todos los campos eléctricos. El decaimiento de la ACA es impulsado por el desfase de la precesión coherente del campo cristalino, mientras que el decaimiento itinerante está limitado por la dispersión.

• Mecanismos de Relajación:

  • Escalamiento ACA: La tasa de relajación orbital de la ACA escala aproximadamente de forma lineal con la tasa de relajación de momento ($\tau_p^{-1}$), imitando el mecanismo Elliott-Yafet (EY) pero impulsado por la interrupción de la precesión del campo cristalino.
  • Escalamiento Itinerante: La relajación itinerante está determinada por la distribución de la textura orbital en el espacio de momentos. En el siliceno, la dispersión entre valles domina, conectando bolsas de OAM opuesto. Aunque la contribución relativa de la dispersión entre valles frente a la intra-valle cambia con el campo eléctrico, la vida media total permanece casi independiente del campo.

Significación y Reivindicaciones
El artículo sostiene que la amplia gama de longitudes de difusión orbital predichas y medidas se origina en la definición misma del OAM. Al demostrar que los modelos basados en ACA producen vidas medias al menos un orden de magnitud más cortas que aquellos que incluyen contribuciones itinerantes, los autores concluyen que:

1. El OAM itinerante es la cantidad gobernante: En los sistemas estudiados (MoS$_2$ y siliceno), el componente itinerante, asociado con la estructura geométrica de los estados de Bloch, dicta la relajación y difusión orbital real.
2. Limitaciones de los Modelos Locales: Las aproximaciones centradas en átomos fallan al capturar la naturaleza de larga duración de la dinámica orbital porque omiten las contribuciones itinerantes y están dominadas por la precesión ultrafast del campo cristalino.
3. Aplicabilidad General: Los hallazgos sugieren que los estudios teóricos de la dinámica orbital y la interpretación de los datos experimentales deben incluir las contribuciones del OAM itinerante. Esto se afirma esencial para materiales 2D como MoS$_2$ y siliceno, y es probablemente aplicable a la familia más amplia de materiales 2D.
4. Dependencia del Material: En los metales, el equilibrio entre los canales atómico e itinerante se espera que dependa del carácter orbital (localizado vs. deslocalizado) y de la calidad de la muestra (desorden), lo que podría explicar la variabilidad en las longitudes de difusión experimental observadas en metales de transición.

El estudio establece una imagen microscópica unificada donde la relajación orbital no es un único mecanismo, sino una competencia entre los efectos locales del campo cristalino y las contribuciones geométricas itinerantes, siendo estas últimas las que a menudo dominan las propiedades de transporte observables.