Real-space first-principles approach to orbitronic phenomena in metallic multilayers

Resumen Técnico: Enfoque de primeros principios en el espacio real para fenómenos orbitrónicos en multicapas metálicas

Planteamiento del problema
El campo de la orbitrónica busca explotar el grado de libertad orbital de los electrones para el procesamiento de información y aplicaciones espintrónicas. Si bien se han logrado avances significativos en la comprensión de los efectos orbitrónicos en materiales cuánticos de volumen y sistemas 2D, los metales de transición (MT) siguen siendo una plataforma central tanto para la investigación fundamental como para las aplicaciones tecnológicas debido a sus complejas estructuras de banda $d$ y su fuerte hibridación. Sin embargo, los enfoques teóricos estándar suelen depender de cálculos de la teoría del funcional de la densidad (DFT) en el espacio de momentos, seguidos de la construcción de Hamiltonianos de red de enlaces mediante funciones de Wannier o proyecciones de orbitales pseudo-atómicos. Estos métodos enfrentan desafíos cuando se aplican a heteroestructuras complejas, sistemas desordenados o interfaces donde la simetría de traslación se rompe. Específicamente, existe la necesidad de un marco escalable de primeros principios que pueda calcular directamente el transporte y la acumulación de Hall orbital (de espín) en el espacio real, acomodando naturalmente los efectos de tamaño finito, la rugosidad de la interfaz y las imperfecciones estructurales sin aproximaciones intermedias.

Metodología
Los autores desarrollan un método de primeros principios en el espacio real basado en el método de Orbitales Muffin-Tin Lineal en el Espacio Real dentro de la Aproximación de la Esfera Atómica (RS-LMTO-ASA). Este enfoque retiene el Hamiltoniano de DFT completo en el espacio real, evitando la necesidad de esquemas de proyección intermedios. La metodología consta de dos componentes principales:

1. Cálculo de la estructura electrónica: Se utiliza el método RS-LMTO-ASA para obtener estructuras electrónicas autoconsistentes y Hamiltonianos en el espacio real. Para resolver el problema de autovalores y calcular la Densidad Local de Estados (LDOS) de manera eficiente, los autores emplean un esquema de recursión basado en polinomios de Chebyshev, el cual ofrece una estabilidad numérica superior en comparación con el algoritmo de Lanczos tradicional, particularmente para sistemas con contribuciones sustanciales de orbitales $s$ y $p$. La LDOS se calcula mediante una expansión de Chebyshev de las funciones de Green, utilizando el Método de Polinomios de Kernel (KPM) con un kernel de Jackson para asegurar la convergencia y obtener funciones espectrales suaves.
2. Transporte de respuesta lineal: Las cantidades de transporte, tales como las conductividades de Hall orbital y de espín y las acumulaciones fuera del equilibrio, se calculan directamente en el espacio real utilizando el formalismo de Kubo-Bastin. Los autores emplean una expansión de polinomios de Chebyshev de las funciones de Green para evaluar la fórmula de Kubo-Bastin. Esto permite el cálculo sistemático de coeficientes de respuesta lineal (conductividades y susceptibilidades magnetoeléctricas) sin depender de evaluaciones estocásticas de la traza. En su lugar, los observables se computan localmente en sitios representativos, explotando la simetría de traslación en las direcciones periódicas mientras se incorporan plenamente las contribuciones de todos los demás sitios mediante la recursión de Chebyshev.

Contribuciones clave

• Desarrollo de un marco escalable: El artículo presenta una metodología que escala linealmente con el número de átomos no equivalentes, permitiendo la simulación de heteroestructuras y multicapas a gran escala con geometrías complejas.
• Transporte directo en el espacio real: Al integrar la expansión de Chebyshev con el Hamiltoniano RS-LMTO-ASA, el método calcula las propiedades de transporte directamente en el espacio real, eliminando la necesidad de proyecciones de Wannier y permitiendo la inclusión natural de condiciones de contorno abiertas y desorden.
• Análisis resuelto por capas: El enfoque proporciona acceso a resolución atómica de los perfiles espaciales de la acumulación de espín y de orbital a través de interfaces y dentro de placas finitas, distinguiendo entre contribuciones de volumen y de superficie.

Resultados
El método fue aplicado a una serie de heteroestructuras basadas en metales de transición y sistemas elementales:

• Metales de transición de volumen: Los cálculos para fcc Pt, fcc Ti, fcc Cu, fcc Ni y bcc W reprodujeron las tendencias establecidas en la conductividad de Hall orbital (OHC) y la conductividad de Hall de espín (SHC). Los resultados confirmaron que los orbitales $d$ dominan la respuesta cerca del nivel de Fermi en metales de transición pesados como Pt y W, mostrando también un buen acuerdo con estudios previos en el espacio de momentos (por ejemplo, Salemi et al., Go et al.).
• Ni ferromagnético: En el Ni ferromagnético, el método capturó con éxito la ruptura de la simetría de inversión temporal, revelando valores finitos para ambos componentes (pares e impares bajo inversión temporal) del tensor de conductividad (por ejemplo, $\sigma^S_{xz}$ y $\sigma^S_{zx}$), consistentes con las predicciones teóricas de los efectos de Hall de espín y orbital anómalos.
• Acumulaciones en multicapas: En láminas de 18 capas de Cu, Ti, Pt, W y Ni, el estudio demostró que ocurre una acumulación sustancial de orbital y de espín cerca de las superficies, decayendo hacia el centro de la lámina. Notablemente, en el Tungsteno (W), el signo de la acumulación orbital superficial difiere del signo de la OHC de volumen. Esta discrepancia se atribuyó a diferencias significativas entre la densidad de estados (DOS) de la superficie y del volumen cerca del nivel de Fermi, resaltando el papel crítico de la estructura electrónica de la superficie.
• Heteroestructuras: Las simulaciones de bicapas consistentes en dos capas ferromagnéticas (Fe, Co, Ni) adyacentes a nueve capas de un metal de transición no magnético (Ti o W) revelaron que las mayores acumulaciones ocurren cerca de la interfaz. Las heteroestructuras basadas en Fe exhibieron la acumulación orbital más fuerte, particularmente en los sistemas Fe/Ti. El estudio también señaló que, si bien la acumulación es un indicador del torque, el torque real está gobernado por la divergencia del flujo de momento angular, una distinción que los autores señalan requiere cálculos de flujo futuros.

Significancia
El artículo afirma que la metodología desarrollada en el espacio real proporciona un marco escalable y flexible para simulaciones realistas de fenómenos de transporte orbital. Su principal significancia radica en su capacidad para modelar sistemas donde la inhomogeneidad espacial y los efectos de tamaño finito son cruciales, tales como interfaces y multicapas desordenadas, los cuales son difíciles de tratar con los enfoques estándar en el espacio de momentos. Los autores demuestran que puede surgir una acumulación sustancial de orbital (espín) incluso en sistemas centrosimétricos, impulsada por asimetrías en la estructura de bandas y dispersión interfacial. Además, este trabajo establece que la relación entre las conductividades de volumen y las acumulaciones de superficie no es trivial y está fuertemente influenciada por la estructura electrónica de la superficie. La metodología se presenta como una herramienta poderosa para investigar fenómenos de la orbitrónica en heteroestructuras complejas, ofreciendo acceso directo a las respuestas locales y allanando el camino para el estudio de sistemas con simetría reducida y desorden, relevantes para la próxima generación de dispositivos espintrónicos y orbitrónicos.