Negative and Zero Linear Compressibility in MCN (M = Ag, Au, Cu): A First-Principles Study

Resumen Técnico: Compresibilidad Lineal Negativa y Cero en MCN (M = Ag, Au, Cu)

Planteamiento del Problema
La compresibilidad lineal negativa (NLC) y la compresibilidad lineal cero (ZLC) son fenómenos mecánicos poco comunes donde un eje cristalográfico se expande o permanece constante, respectivamente, bajo presión hidrostática. Aunque se conocen en materiales de estructura de armazón (framework), óxidos inorgánicos y cristales moleculares, lograr una NLC/ZLC fuerte y estable en un amplio rango de presión sigue siendo un desafío debido a las inestabilidades mecánicas o efectos débiles en los materiales candidatos. Los mecanismos existentes, como el de bisagra tipo "wine-rack" (estante de vino), la inclinación de poliedros y el deslizamiento de capas, están bien documentados pero son limitados. La familia de los cianuros metálicos (MCN, donde M = Ag, Au, Cu) ha mostrado un comportamiento anómalo, incluyendo expansión térmica negativa y una NLC gigante en fases específicas (por ejemplo, CuCN $C222_1$), sin embargo, el potencial para ZLC o NLC en las fases $R\bar{3}m$ y $P6mm$ de AgCN, AuCN y CuCN sigue sin explorarse suficientemente.

Metodología
El estudio emplea cálculos de primeros principios basados en la Teoría del Funcional de la Densidad (DFT) utilizando el paquete VASP (Vienna Ab-initio Simulation Package).

• Estructuras Investigadas: Los autores analizaron seis estructuras: fases reportadas experimentalmente (AgCN $R\bar{3}m$, AuCN $P6mm$, CuCN $R\bar{3}m$) y tres fases no reportadas previamente (AgCN $P6mm$, CuCN $P6mm$, AuCN $R\bar{3}m$). Las fases no reportadas fueron generadas mediante la sustitución de sitios metálicos en estructuras conocidas.
• Funcionales: Se probaron seis funcionales de intercambio y correlación para determinar su fiabilidad: PBE, PBE-D3, PBEsol, PBEsol-D3, r2SCAN y r2SCAN-D3. Basándose en las comparaciones de los parámetros de red con los datos experimentales, se seleccionaron PBE-D3, PBEsol y r2SCAN para el análisis detallado de las propiedades.
• Cálculos: Se realizaron relajaciones estructurales completas. Las constantes elásticas ($C_{ij}$) se computaron mediante métodos de diferencia finita para derivar las compresibilidades lineales ($\beta$), los módulos de Young ($E$) y los módulos de compresibilidad ($K$). Se aplicaron presiones hidrostáticas de hasta 13 GPa para evaluar la estabilidad y la evolución de la respuesta. Las curvas de dispersión de fonones se calcularon para investigar la rigidez de los enlaces y las vibraciones atómicas.

Contribuciones Clave y Resultados

1. Predicción de Nuevas Fases: El estudio predice la existencia de fases energéticamente competitas no reportadas experimentalmente: la fase $P6mm$ para AgCN y CuCN, y la fase $R\bar{3}m$ para AuCN. Las diferencias de energía entre estas fases propuestas y sus contrapartes experimentales son pequeñas (a menudo < 50 meV/u.f.), lo que sugiere que son candidatos viables.
2. Extrema Anisotropía Elástica y Compresibilidad Anómala: Las seis fases de MCN investigadas exhiben una extrema anisotropía elástica, donde la constante elástica a lo largo de la dirección $c$ ($C_{33}$) es un orden de magnitud mayor que en el plano $ab$($C_{11}$).
   • Predicción de ZLC/NLC: Esta anisotropía conduce a una compresibilidad lineal cercana a cero o negativa a lo largo del eje $c$. Específicamente, el AgCN en la fase $R\bar{3}m$ exhibe NLC (expansión del eje $c$) bajo presión. Otras fases muestran un comportamiento cercano a ZLC.
   • Estabilidad de Presión: Estas respuestas anómalas persisten en un amplio rango de presión (hasta 13 GPa) sin inestabilidad estructural.
3. Identificación de un Mecanismo Novedoso: El artículo identifica un mecanismo distinto de los motivos tradicionales de "wine-rack" o de bisagra. Las estructuras de MCN forman una geometría de "bosque de bambú" compuesta por varas (rods) de MCN unidimensionales y rígidas.
   • Varas Rígidas: Las varas están mantenidas por enlaces extremadamente fuertes entre átomos ligeros (M-C-N), evidenciado por modos de fonones de alta frecuencia (hasta 68 THz).
   • Empaquetamiento Disperso: Las varas interactúan mediante fuerzas muy débiles en el plano perpendicular, evidenciado por modos de fonones de muy baja frecuencia (~1.2 THz).
   • Respuesta: Bajo presión hidrostática, la reducción de volumen es acomodada casi enteramente por el movimiento de las varas para acercarse entre sí (aumentando la densidad de empaquetamiento en el plano $ab$), en lugar de comprimir las propias varas. En casos específicos (AgCN $R\bar{3}m$), este reordenamiento provoca que las varas rígidas se expandan ligeramente a lo largo de su longitud (eje $c$) para mantener los ángulos de enlace, lo que resulta en NLC.

Significancia
Los autores afirman que este trabajo expande la familia relativamente pequeña de materiales conocidos por exhibir ZLC o NLC. Al identificar el mecanismo de "bosque de bambú" —donde el empaquetamiento disperso de cadenas 1D rígidas permite una compresibilidad anómala— el estudio proporciona un nuevo origen microscópico para estos fenómenos, distinto de los mecanismos reportados anteriormente. Los hallazgos ofrecen una comprensión más profunda de la relación entre la anisotropía estructural y la respuesta mecánica, guiando potencialmente el diseño de materiales para componentes mecánicos de precisión y sensores de presión, aunque el artículo se centra principalmente en la comprensión fundamental más que en aplicaciones de dispositivos específicos. El estudio también destaca la utilidad de funcionales de DFT específicos (r2SCAN y PBE-D3) para modelar con precisión estos complejos sistemas de cianuros.