Charge-Transfer Electronic Structure of NiX$_2$ (X = S, Se)

Resumen Técnico: Estructura Electrónica de Transferencia de Carga de NiX$_2$ (X = S, Se)

Planteamiento del Problema
La naturaleza electrónica de los dicalcogenuros de níquel (NiS$_2$ y NiSe$_2$) sigue siendo objeto de debate respecto a la interacción entre las correlaciones electrónicas, el enlace covalente y el papel específico de los dímeros de calcógeno. Mientras que el NiS$_2$ se clasifica tradicionalmente como un aislante de transferencia de carga (CT) con una brecha pequeña (~0.2 eV) y el NiSe$_2$ como un metal, el mecanismo que impulsa la transición metal-aislante (MIT) en la solución sólida NiS$_{2-x}$Se$_x$ no está totalmente resuelto. Estudios previos de la Teoría del Funcional de la Densidad combinada con la Teoría de Campo Medio Dinámica (DFT+DMFT) situaron los estados del dímero antibondante de azufre-selenio por debajo de la banda de Hubbard superior (UHB) del Ni 3d. Sin embargo, datos espectroscópicos experimentales recientes sugieren discrepancias con estos modelos, insinuando una energía de transferencia de carga significativamente menor o incluso negativa y un ordenamiento diferente de los estados electrónicos. La posición relativa precisa de los estados Ni 3d y los estados $p$ del calcógeno, gobernada por la corrección de doble conteo ($\mu_{dc}$), requiere una reevaluación para reconciliar la teoría con los espectros experimentales de la banda de valencia, de niveles de núcleo y de absorción.

Metodología
Los autores emplearon un marco de trabajo DFT+DMFT para investigar las estructuras electrónicas de NiS$_2$ y NiSe$_2$.

• Configuración Computacional: Los cálculos de DFT se realizaron utilizando el paquete WIEN2k dentro de la Aproximación de Densidad Local (LDA) para las estructuras cristalinas experimentales. Se construyó un modelo de enlace fuerte (TB) utilizando wien2wannier y wannier90, incorporando orbitales Ni 3d, y calcógeno $p_\pi$ y $p_\sigma$.
• Parámetros de Interacción: Las interacciones de la capa Ni 3d fueron tratadas con un Hubbard $U = 5.0$ eV y una acoplamiento de Hund $J = 1.0$ eV. El problema de impurezas se resolvió utilizando el resolvedor de Monte Carlo de tiempo continuo de expansión de hibridación (CT-QMC) a 580 K.
• Estrategia de Corrección de Doble Conteo: A diferencia de estudios previos que utilizaron esquemas fijos, los autores trataron la corrección de doble conteo ($\mu_{dc}$) como un parámetro ajustable. Optimizaron $\mu_{dc}$ variando sistemáticamente el parámetro de diferencia de energía $\Delta_{dp}$ (definido como la diferencia de energía entre los orbitales Ni 3d y calcógeno $p$) para reproducir los datos experimentales de fotoemisión de la banda de valencia (PES).
• Validación Experimental: Los espectros teóricos se compararon con:
  • PES de rayos X duros (HX-PES, $h\nu = 6.5$ keV) y PES de rayos X blandos (SX-PES, $h\nu = 1.2$ keV) medidos en SPring-8 y Photon Factory.
  • Espectroscopía de fotoelectrones de rayos X (XPS) de nivel de núcleo Ni 2p y Espectroscopía de Absorción de Rayos X (XAS) en los bordes Ni $L_{2,3}$.
• Variaciones del Modelo: El estudio comparó "modelos completos" (que incluyen orbitales de calcógeno $p_\sigma$) con modelos reducidos "Ni-$p_\pi$" (excluyendo $p_\sigma$) para aislar el papel de los estados antibondantes del dímero.

Resultados Clave

1. Energía de Transferencia de Carga Optimizada: Al ajustar los espectros de la banda de valencia a los datos experimentales de PES, se determinó que el $\Delta_{dp}$ óptimo es de 1.82 eV para NiS$_2$ y 1.58 eV para NiSe$_2$. Estos valores indican una energía de transferencia de carga más pequeña de lo asumido anteriormente en los primeros estudios de DFT+DMFT, lo cual es consistente con análisis recientes de modelos de clúster.
2. Reproducción de Características Espectrales: El modelo refinado reproduce con éxito:
   • Los espectros de la banda de valencia tanto para HX-PES como para SX-PES, capturando correctamente las variaciones de la sección eficaz dependiente de los orbitales.
   • Los espectros de PES de nivel de núcleo Ni 2p y XAS de NiS$_2$, incluyendo las características satélite específicas (St1 y St2) en el componente Ni 2$p_{3/2}$, las cuales solo se capturaron cerca de los valores óptimos de $\Delta_{dp}$.
3. Estructura Electrónica y Transferencia de Carga:
   • Los parámetros optimizados revelan una transferencia de carga ligando-a-Ni sustancial, lo que resulta en una ocupación de Ni 3d de aproximadamente 8.5 (significativamente mayor que la 8.2 encontrada en NiO).
   • Esta mayor mezcla entre las configuraciones $|t_6e_2\rangle$ y $|t_6e_3p_\pi\rangle$ conduce a una reducción del momento local de Ni a un valor apantallado de 1.38 $\mu_B$, consistente con las observaciones experimentales a baja temperatura.
4. Mecanismo de la Transición Metal-Aislante (MIT):
   • NiS$_2$: El sistema es un aislante donde la UHB está separada de la banda de transferencia de carga (CT). Los estados del dímero antibondante S-S ($p_\sigma$) se encuentran por encima del nivel de Fermi ($E_F$) y actúan principalmente como espectadores pasivos en el escenario óptimo de estado Ni $d$ profundo.
   • NiSe$_2$: El desdoblamiento del dímero Se-Se es más débil, situando el estado antibondante $p_\sigma$ a una energía menor en comparación con el NiS$_2$. Aunque el subsistema Ni-$p_\pi$ por sí solo podría permanecer con brecha para ciertos parámetros, la débil hibridación con la banda antibondante Se $p_\sigma$ proporciona el "impulso" necesario para empujar al sistema hacia el estado metálico. La banda Se $p_\sigma$ cruza el $E_F$, cerrando la brecha.
   • El estudio aclara que la MIT en la serie NiS$_{2-x}$Se$_x$ es impulsada por la posición relativa y la interacción entre la banda de Hubbard superior de Ni y los estados antibondantes del dímero de calcógeno, más que por un simple control de ancho de banda.

Significancia y Reivindicaciones
El artículo afirma resolver las discrepancias entre los modelos teóricos anteriores y los datos espectroscópicos experimentales recientes mediante el establecimiento de una estructura electrónica realista a través de la optimización de la corrección de doble conteo. Los autores aseveran que su modelo proporciona una descripción unificada de la estructura electrónica de baja energía, las excitaciones de nivel de núcleo y los momentos magnéticos para ambos NiS$_2$ y NiSe$_2$.

La principal significancia radica en:

• Corrección de la Escala de Energía de Transferencia de Carga: Establecer una energía de transferencia de carga más pequeña que se alinea con los datos experimentales de nivel de núcleo y de banda de valencia, desafiando las suposiciones previas sobre la separación de energía entre la UHB y el tope de la banda de valencia.
• Clarificación del Mecanismo de la MIT: Demostrar que la transición de NiS$_2$ aislante a NiSe$_2$ metálico está gobernada por la alineación de energía específica de los estados del dímero de calcógeno antibondante respecto a la UHB de Ni 3d. El estudio destaca que la naturaleza metálica de NiSe$_2$ surge del solapamiento de la banda antibondante Se $p_\sigma$ con el nivel de Fermi, facilitado por el reducido desdoblamiento del dímero y la energía de transferencia de carga específica.
• Validación Metodológica: Probar que el ajuste de $\mu_{dc}$ contra la PES experimental es una estrategia robusta para determinar la estructura electrónica de materiales correlacionados, como lo evidencia la reproducción simultánea de los espectros de banda de valencia, de nivel de núcleo y de absorción.

Los autores concluyen que el NiS$_2$ y el NiSe$_2$ poseen fuerzas de hibridación Ni-calcógeno muy similares, con el NiS$_2$ situado cerca de la MIT en el lado aislante, mientras que la posición específica de los estados antibondantes de Se impulsa al NiSe$_2$ hacia el régimen metálico.