Dynamical Hubbard approach to correlated materials: the case of transition-metal monoxides

Resumen Técnico: Enfoque de Hubbard Dinámico para Materiales Correlacionados

Planteamiento del Problema
Los monóxidos de metales de transición (MnO, FeO, CoO, NiO) sirven como sistemas prototípicos para el estudio de correlaciones electrónicas complejas, categorizados como aislantes de Mott y de transferencia de carga. Mientras que la teoría del funcional de la densidad de Kohn-Sham estándar con funcionales semi-locales falla al predecir la naturaleza aislante de CoO y FeO y subestima las brechas de banda en MnO y NiO, los métodos más avanzados enfrentan compromisos significativos.

• Enfoques de Campo Medio Estático: Métodos como DFT+U y funcionales híbridos pueden recuperar los estados aislantes y las magnitudes de la brecha, pero son inherentemente estáticos. No logran capturar las autoenergías dependientes de la frecuencia necesarias para describir la renormalización de bandas, la transferencia de peso espectral y los tiempos de vida de las partículas; firmas clave de las correlaciones electrónicas.
• Teoría de Perturbación de Muchos Cuerpos (MBPT): Los enfoques como $G_0W_0$ y $GW_0$ a menudo tienen dificultades cuando parten de estados fundamentales metálicos de DFT. Aunque extensiones como QSG$\hat{W}$ mejoran la descripción de la brecha, frecuentemente descartan la dependencia completa de la frecuencia debido al costo computacional, lo que limita su capacidad para describir la transferencia de peso espectral resuelta en $k$.
• DFT+DMFT: La Teoría de Campo Medio Dinámico (DFT+DMFT) modela con éxito las correlaciones dinámicas locales y los momentos magnéticos fluctuantes mediante el mapeo del subespacio correlacionado a un modelo de impureza de Anderson. Sin embargo, esto conlleva altos costos computacionales y complejidad numérica, ya que las soluciones exactas escalan exponencialmente con el tamaño del sistema. Además, resultados recientes de DFT+DMFT han predicho brechas anchas para los cuatro monóxidos y no han logrado reproducir ciertas características espectrales observadas en la Espectroscopía de Fotoemisión con Resolución Angular (ARPES), como bandas dispersivas específicas en NiO.

Metodología
Los autores proponen un funcional de "Hubbard Dinámico" (dynH), una generalización de su formulación introducida previamente para manejar sistemas magnéticos de múltiples sitios y polarizados en el espín. El enfoque está diseñado para ser numéricamente eficiente y físicamente transparente, capturando simultáneamente las correlaciones dinámicas.

• Formalismo: El método extiende la formulación de DFT+U rotacionalmente invariante de Dudarev mediante la introducción de un potencial de interacción apantallado, local y dependiente de la frecuencia. La autoenergía se deriva de la derivada funcional del funcional de Hubbard dinámico con respecto a la función de Green local:
  $$\Sigma^{\sigma, I}_{\text{dynH}}(\omega) = -\int d\omega' U_I(\omega') G^{\sigma}_I(\omega + \omega') + \frac{U^{\infty}_I}{2}$$
  Aquí, $U_I(\omega)$ representa el promedio local de la interacción de Coulomb apantallada (calculada mediante la polarización de spin-DFT de RPA), y $U^{\infty}_I$ es el potencial de Coulomb desnudo sin apantallar.
• Implementación:
  • La autoenergía local se traslada hacia arriba (upfolded) al espacio $k$ para resolver la ecuación de Dyson ($G^{-1} = G^{-1}_0 - \Sigma$).
  • La solución utiliza el Método de Inversión Algorítmica sobre Suma de Polos (AIM-SOP). Esta técnica representa tanto la interacción apantallada local como la función de Green como sumas sobre polos de primer orden.
  • Crucialmente, AIM-SOP opera directamente en el eje de frecuencia real, evitando el uso de frecuencias de Matsubara imaginarias y la posterior continuación analítica requerida en los solucionadores de impurezas estándar de DMFT. Esto mejora la precisión del cálculo.
  • Los cálculos del estado fundamental emplean el funcional r2SCAN, el cual predice correctamente las fases aislantes para los cuatro compuestos de manera equitativa.
  • La interacción apantallada efectiva se calcula en el primer orden de la teoría de perturbaciones, con un límite de frecuencia cero de $U(0) = 5.3$ eV.

Resultados Clave
El marco se aplicó a los estados fundamentales de MnO, FeO, CoO y NiO, comparando los resultados con datos experimentales de Fotoemisión (PES) e Fotoemisión Inversa (IPES), así como con cálculos de vanguardia de DFT+DMFT y QSG$\hat{W}$.

• MnO: La Densidad de Estados (DOS) ocupada se alinea con los resultados de DFT+DMFT, mostrando una ligera contracción en el ancho de banda general en comparación con el experimento. La excitación LUMO a 5 eV coincide en posición, ancho e intensidad con los datos de IPES. Los satélites de alta energía siguen siendo esquivos, una limitación compartida con DFT+DMFT y QSG$\hat{W}$.
• FeO: El método predice una brecha pequeña de aproximadamente 2 eV, consistente con los datos de IPES. Esto contrasta con DFT+DMFT y QSG$\hat{W}$, que predicen brechas anchas para los cuatro monóxidos. La disminución de la intensidad después de las primeras excitaciones está bien descrita.
• CoO: La DOS se compara bien con los datos de DFT+DMFT y PES. Una diferencia notable es la distribución de los estados vacíos: DFT+dynH muestra pesos comparables para los estados de simetría $t_{2g}$ y $e_g$, mientras que DFT+DMFT transfiere el peso principalmente al estado de mayor energía.
• NiO: El enfoque logra una excelente concordancia con los datos experimentales de ARPES y los espectros PES/IPES.
  • El ancho de banda general (6 eV) es consistente con ARPES.
  • El método reproduce con éxito el doblete de Zhang-Rice y las bandas planas de Ni(3d)/O(2p).
  • Hallazgo Crítico: A diferencia de DFT+DMFT, que predice un peso espectral evanescente para bandas específicas, DFT+dynH reproduce la banda dispersiva en -8 eV y la banda plana antiferromagnética (AFM) alrededor de -4 eV. Estas características están presentes en los datos experimentales de ARPES, pero anteriormente estaban ausentes en las soluciones de impureza de DMFT. Los autores atribuyen la ligera subestimación de las bandas AFM al punto de partida r2SCAN, la falta de auticonsistencia o la naturaleza aproximada del término de doble conteo.

Significancia y Reivindicaciones
El artículo sostiene que el enfoque DFT+dynH ofrece una alternativa "simple y transparente" a los métodos computacionalmente costosos como DFT+DMFT. Sus principales contribuciones son:

1. Generalización: Extensión exitosa del funcional de Hubbard dinámico a sistemas magnéticos de múltiples sitios y polarizados en el espín.
2. Precisión vs. Costo: Lograr resultados comparables o superiores a las teorías de vanguardia (DFT+DMFT, QSG$\hat{W}$) para propiedades espectrales, incluyendo la renormalización de bandas y la transferencia de peso espectral, pero con un costo computacional y complejidad numérica significativamente menores.
3. Resolución de Discrepancias: Corrección de fallos específicos de teorías previas, como la predicción de una brecha pequeña para FeO y la recuperación de características espectrales en NiO que estaban ausentes en las soluciones de DMFT.
4. Perspectiva Física: Proporcionar acceso a las firmas fundamentales de las correlaciones electrónicas (tiempos de vida de las partículas, satélites, transferencia de peso espectral) sin la escala prohibitiva de los solucionadores de impurezas exactos.

Los autores posicionan a DFT+dynH como una herramienta prometedora para simular materiales correlacionados realistas y complejos, particularmente aquellos con ordenamientos orbitales de simetría rota e hibridaciones, donde los métodos que emplean soluciones exactas pueden resultar prohibitivos.