High-pressure electride superconductor Li5N for multifunctional applications: A theoretical insight into the physical properties

Resumen Técnico: Electrido superconductor de alta presión Li5N para aplicaciones multifuncionales

Planteamiento del problema
Los electridos, sólidos iónicos donde los electrones actúan como aniones confinados en cavidades intersticiales, exhiben propiedades físicas únicas como alta conductividad y actividad catalítica. Mientras que los electridos orgánicos suelen ser térmicamente inestables y sensibles a factores ambientales, los electridos inorgánicos ofrecen una mayor estabilidad. Sin embargo, los electridos tridimensionales (3D) establecidos experimentalmente son escasos, y sus propiedades bajo condiciones extremas no se comprenden completamente. Específicamente, aunque se predijo previamente que el compuesto Li5N sería un electrido 3D estable a altas presiones (80–100 GPa), su comportamiento más allá de los 100 GPa, particularmente respecto a su naturaleza de electrido, superconductividad y propiedades físicas integrales (mecánicas, ópticas, térmicas), permaneció en gran medida inexplorado. Además, la coexistencia de estados de electrido y superconductividad representa una frontera significativa en la física de la materia condensada; no obstante, los electridos superconductores de alta temperatura de transición (alta Tc) son poco comunes en condiciones ambientales.

Metodología
Este estudio emplea cálculos de primeros principios basados en la Teoría de la Funcional de la Densidad (DFT) para investigar las propiedades físicas del Li5N hexagonal (grupo espacial P6/mmm) bajo presiones hidrostáticas altas que oscilan entre 150 y 350 GPa. Las simulaciones se realizaron utilizando el paquete CASTEP con el funcional de intercambio-correlación PBEsol (Perdew-Burke-Ernzerhof) con corrección de sólidos. Los parámetros computacionales clave incluyeron un corte de base de ondas planas de 500 eV y una rejilla k de Monkhorst-Pack de 14×14×9.

La metodología comprendió:

• Análisis estructural y termodinámico: Cálculo de los parámetros de red de equilibrio, energía de formación ($E_f$) y energía cohesiva ($E_c$) para evaluar la estabilidad.
• Estabilidad dinámica: Se computaron curvas de dispersión de fonones y densidades de estados (PHDOS) utilizando la teoría de perturbación de la densidad funcional (DFPT) y el método de supercelda de desplazamientos finitos para identificar los rangos de estabilidad dinámica.
• Propiedades mecánicas: Las constantes de rigidez elástica ($C_{ij}$) se derivaron mediante el método de esfuerzo-deformación. Los módulos policristalinos (Bulk $B$, Cizalla $G$, Young $E$), la razón de Poisson ($\nu$), la razón de Pugh ($G/B$) y la dureza de Vickers ($H_V$) se estimaron utilizando la aproximación de Voigt–Reuss–Hill (VRH). La anisotropía elástica se analizó mediante factores de anisotropía de cizalladura y construcciones de superficie 3D.
• Propiedades térmicas: Las temperaturas de Debye ($\Theta_D$), velocidades de sonido, conductividades térmicas mínimas ($k_m$), conductividades térmicas de red ($k_l$) y temperaturas de fusión ($T_m$) se calcularon utilizando los modelos de Anderson, Cahill-Pohl, Clarke y Slack.
• Propiedades electrónicas y ópticas: Se realizaron cálculos de estructuras de bandas electrónicas, densidad de estados totales y parciales (TDOS/PDOS) y análisis de población de Mulliken para determinar la naturaleza del enlace y el carácter metálico. Las propiedades ópticas, incluyendo funciones dieléctricas, índices de refracción y reflectividad, se modelaron incorporando un término de Drude para dar cuenta de las transiciones intrabanda metálicas.
• Superconductividad: La constante de acoplamiento electrón-fonón ($\lambda_{ep}$) y el pseudopotencial de Coulomb ($\mu^*$) se estimaron utilizando la ecuación de McMillan y la fórmula de Bennemann-Garland, utilizando las temperaturas críticas ($T_c$) predichas teóricamente de la literatura previa.

Resultos Clave

• Estabilidad: El Li5N es termodinámicamente estable dentro del rango de 150–350 GPa, evidenciado por energías de formación negativas y energías cohesivas positivas. La estabilidad dinámica se confirma entre 100 GPa y 382 GPa; por debajo de 100 GPa y por encima de 382 GPa, el sistema exhibe frecuencias de fonones imaginarias, lo que indica inestabilidad.
• Comportamiento mecánico: El compuesto es mecánicamente estable y exhibe un comportamiento dúctil a través del rango de presión estudiado, respaldado por razones de Poisson (0.29–0.34), razones de Pugh (<0.57) y presiones de Cauchy positivas. El material es anisotrópico, con constantes elásticas que aumentan significativamente con la presión. La dureza de Vickers aumenta de 26.62 GPa a 150 GPa a 34.48 GPa a 350 GPa, acercándose al umbral de los materiales superduros.
• Estructura electrónica: El Li5N retiene su carácter metálico a lo largo del rango de 150–350 GPa, sin apertura de brecha de banda (band gap) en el nivel de Fermi. La densidad de estados en el nivel de Fermi ($N(E_F)$) disminuye de 0.207 estados/eV/f.u. a 150 GPa a 0.077 estados/eV/f.u. a 350 GPa debido al ensanchamiento de bandas inducido por la presión. El análisis de Mulliken y el mapeo de densidad de carga revelan una transición donde el enlace iónico se debilita y el carácter covalente se fortalece conforme aumenta la presión.
• Propiedades térmicas y acústicas: Las velocidades de sonido y las temperaturas de Debye aumentan con la presión, alcanzando ~2134 K a 350 GPa, lo que sugiere una alta conductividad térmica. La temperatura de fusión es extremadamente alta, variando entre ~5156 K y ~9457 K. La conductividad térmica de red aumenta con la presión pero disminuye con la temperatura.
• Respuesta óptica: El material muestra características ópticas metálicas típicas, incluyendo una función dieléctrica real negativa a bajas energías (comportamiento tipo Drude), alta reflectividad estática (0.97) y una fuerte absorción en la región ultravioleta (10$^5$ cm$^{-1}$). La respuesta óptica es anisotrópica, diferenciándose entre las direcciones de polarización <100> y <001>.
• Superconductividad: El estudio indica que, si bien el Li5N posee características superconductoras, la fuerza de acoplamiento electrón-fonón ($\lambda_{ep}$) y la fuerza de interacción efectiva disminuyen significativamente con el aumento de la presión (de 0.569 a 150 GPa a 0.118 a 350 GPa). Esta supresión se atribuye a la reducción de $N(E_F)$, lo que sugiere un efecto perjudicial sobre la temperatura de transición superconductora a presiones más altas.

Significancia y Reivindicaciones
Los autores afirman que este estudio proporciona una visión teórica integral de las propiedades multifuncionales del electrido Li5N bajo alta presión, llenando un vacío dejado por estudios previos que se centraron principalmente en la estabilidad estructural hasta los 100 GPa. El trabajo establece al Li5N como un electrido dinámica y mecánicamente estable en el rango de 150–350 GPa, caracterizado por ductilidad, alta dureza y conductividad metálica.

El artículo destaca la variación sistemática de las propiedades estructurales, mecánicas, optoelectrónicas y superconductoras como una función de la presión, presentando una estrategia computacional efectiva para ajustar estas propiedades. Aunque el estudio señala que el Li5N exhibe superconductividad, concluye modestamente que el aumento de la presión suprime el acoplamiento electrón-fonón, lo que probablemente reduce la temperatura de transición superconductora. Los autores posicionan sus resultados como una guía valiosa para la futura investigación experimental, particularmente en espectroscopia y en la síntesis de fases de electridos de alta presión, señalando que los datos experimentales para muchos de estos parámetros específicos (por ejemplo, frecuencias de fonones, conductividad térmica de red) no están disponibles actualmente.