Pressure-Driven Evolution of the Physical Properties of Layered Y2GaB and Y2InB MAB Phases: A Comparative First-Principles Study

Resumen Técnico: Evolución Impulsada por la Presión de las Propiedades Físicas de las Fases MAB de Capas de Y2GaB e Y2InB

Planteamiento del Problema
Los boruros ternarios laminados conocidos como fases MAB exhiben una combinación única de comportamiento electrónico metálico y rigidez de tipo cerámico debido a su estructura nanolaminada, donde coexisten enlaces metal-boro fuertes en el plano con interacciones intercapas más débiles. Si bien estudios previos han caracterizado compuestos MAB específicos a presión ambiente, existe una falta de comprensión unificada sobre cómo la presión hidrostática y la sustitución en el sitio A (específicamente Ga frente a In) influyen en su estabilidad estructural, anisotropía mecánica, topología electrónica y propiedades termofísicas. Además, trabajos computacionales previos sobre Y2GaB e Y2InB reportaron modos de fonones de frecuencia negativa, lo que plantea interrogantes sobre su estabilidad dinámica y la sensibilidad de estas evaluaciones a los parámetros computacionales (tratamientos de celda unitaria frente a supercelda). Este estudio aborda estas brechas mediante la comparación sistemática de la evolución dependiente de la presión de Y2GaB e Y2InB de 0 a 20 GPa para establecer ventanas de presión mecánicamente viables y aclarar las relaciones estructura-propiedad.

Metodología
El estudio emplea cálculos de primeros principios basados en la teoría del funcional de la densidad (DFT) de Kohn-Sham utilizando el código CASTEP con la aproximación de gradiente generalizado (GGA) de Perdew-Burke-Ernzerhof (PBE). Las características metodológicas clave incluyen:

• Optimización Estructural: Se optimizaron estructuras hexagonales $P6_3/mmc$ bajo presiones hidrostáticas (0, 5, 10, 15 y 20 GPa) utilizando el algoritmo BFGS.
• Estabilidad Elástica: Las constantes elásticas se derivaron mediante el método de deformación-esfuerzo finito. La estabilidad se evaluó utilizando criterios corregidos por presión para la simetría hexagonal, distinguiendo entre estabilidad de deformación homogénea y estabilidad de fonones.
• Derivación de Propiedades: Se calcularon los módulos policristalinos (promedios de Voigt-Reuss-Hill), velocidades acústicas, temperaturas de Debye y estimaciones de conductividad térmica (modelos de Clarke y Slack) para los estados elásticamente estables.
• Análisis Electrónico y Óptico: Las estructuras de bandas, las densidades de estados (DOS) y las superficies de Fermi se computaron con y sin acoplamiento espín-órbita (SOC) utilizando Quantum ESPRESSO. Las propiedades ópticas se resolvieron para polarizaciones de campo eléctrico a lo largo de [100] y [001].
• Enlace y Dinámica: Se utilizaron mapas de diferencia de densidad electrónica, análisis de carga de Mulliken/Hirshfeld y poblaciones de solapamiento de enlaces para caracterizar el enlace. Las dispersiones de fonones se calcularon mediante el método de desplazamiento finito tanto para configuraciones de celda unitaria como para superceldas expandidas para evaluar la estabilidad dinámica y la sensibilidad numérica.

Resultados Clave

• Evolución Estructural: Ambos compuestos mantienen la estructura $P6_3/mmc$ bajo presión. La compresión es anisotrópica, con el parámetro de red basal ($a$) contrayéndose más significativamente que el parámetro axial ($c$). En consecuencia, la relación $c/a$ aumenta para ambos, y el volumen de la celda unitaria disminuye monótonamente (aproximadamente un 21.8% para Y2GaB y un 20.3% para Y2InB entre 0 y 20 GPa).
• Estabilidad Elástica y Propiedades Mecánicas:
  • Ventanas de Estabilidad: Y2GaB satisface los criterios de estabilidad elástica corregidos por presión a 0, 5 y 10 GPa, pero se vuelve inestable a 15 y 20 GPa. Y2InB es elásticamente estable solo a presión ambiente (0 GPa).
  • Respuesta Mecánica: Dentro de su rango estable, Y2GaB exhibe un aumento en los módulos de bulk, cizalla y Young con la presión. Ambos compuestos se predicen como dúctiles (relación de Pugh $G/B < 0.57$). Y2GaB muestra una alta anisotropía elástica a 0 GPa ($A^U = 5.25$), la cual disminuye con la presión, mientras que Y2InB es menos anisotrópico a presión ambiente.
  • Propiedades Termofísicas: Para el caso estable de Y2GaB, la compresión aumenta las velocidades del sonido, la temperatura de Debye y las temperaturas de fusión estimadas. La conductividad térmica de red estimada por Slack a 300 K aumenta de 0.605 W m⁻¹ K⁻¹ (0 GPa) a 1.543 W m⁻¹ K⁻¹ (10 GPa). Ambos compuestos exhiben conductividades térmicas mínimas muy bajas, lo que sugiere potencial para aplicaciones de barrera térmica a bajas temperaturas, limitado por sus puntos de fusión relativamente bajos.
• Estructura Electrónica: Ambos materiales son metales multibanda. Y2GaB tiene nueve hojas de superficie de Fermi, e Y2InB tiene ocho a presión ambiente. El acoplamiento espín-órbita (SOC) modifica las dispersiones de banda locales pero no abre una brecha ni altera el número de cruces de la superficie de Fermi. No se observa una transición de Lifshitz (reconstrucción topológica de la superficie de Fermi) entre 0 y 10 GPa para Y2GaB, aunque el bolsillo central se contrae significativamente. Los estados electrónicos cerca del nivel de Fermi están dominados por los orbitales Y-4d, con contribuciones secundarias de los orbitales B-2p y los orbitales p del sitio A.
• Carácter de Enlace: El análisis revela un enlace mixto iónico-covalente entre Y y B, caracterizado por la acumulación de carga alrededor de B y el agotamiento alrededor de Y. Los elementos del sitio A (Ga/In) exhiben un enlace metálico más difuso.
• Propiedades Ópticas: Ambos compuestos muestran respuestas Drude metálicas y una fuerte anisotropía óptica entre las polarizaciones [100] y [001]. Actúan como reflectores eficientes en la región visible y absorbentes en el ultravioleta. La presión induce desplazamientos no uniformes en las características ópticas, elevando generalmente las energías características para Y2GaB.
• Dinámica de Fonones: Los cálculos de fonones de celda unitaria muestran frecuencias positivas para ambos compuestos a 0 y 10 GPa, lo que indica estabilidad dinámica dentro de este tratamiento restringido. Sin embargo, los cálculos de supercelda expandida revelan modos imaginarios, lo que resalta la sensibilidad de la estabilidad de los fonones de baja frecuencia al tamaño de la celda y al tratamiento de las constantes de fuerza. Notablemente, Y2InB permanece elásticamente inestable a 10 GPa a pesar de mostrar frecuencias positivas en el espectro de fonones de la celda unitaria, lo que subraya la distinción entre los criterios de estabilidad elástica y dinámica.

Significancia y Reivindicaciones
El artículo establece un marco comparativo para comprender cómo la sustitución del sitio A (Ga frente a In) y la presión hidrostática modulan las propiedades de las fases MAB. Sus principales contribuciones son:

1. Definición de Ventanas de Estabilidad: Identifica rangos de presión específicos donde Y2GaB es mecánicamente viable (0–10 GPa) y aclara que Y2InB es estable solo a presión ambiente bajo los criterios adoptados.
2. Desacoplamiento de Tipos de Estabilidad: El estudio demuestra explícitamente que la estabilidad elástica (resistencia a la deformación homogénea) y la estabilidad dinámica (resistencia a los desplazamientos atómicos) son distintas y deben evaluarse por separado, señalando particularmente la sensibilidad de los cálculos de fonones al tamaño de la supercelda.
3. Relaciones Estructura-Propiedad: Detalla cómo la presión mejora sistemáticamente la rigidez, el transporte térmico y la respuesta óptica en Y2GaB, al tiempo que reduce la anisotropía elástica, proporcionando una base para futuras validaciones experimentales.
4. Alcance Modesto: Los autores declaran explícitamente que no existen derivados de MABene de estos compuestos establecidos experimentalmente; el estudio se centra en las fases parentales como una base necesaria para cualquier evaluación futura de derivados de dimensiones reducidas. Los valores de conductividad térmica se presentan como estimaciones de modelos (Slack/Clarke) en lugar de predicciones de transporte precisas, reconociendo la ausencia de cálculos explícitos de dispersión electrón-fonón y anarmónica.