Pressure-Induced Structural Phase Transitions in MBANP: A First-Principles Study of Elastic, Phonon, and Optical Properties

Resumen Técnico: Transiciones de Fase Estructural Inducidas por Presión en MBANP

Planteamiento del Problema
Los cristales orgánicos no lineales (NLO), específicamente la 2-(α-metilbencilamino)-5-nitropiridina (MBANP), son críticos para las tecnologías fotónicas y optoelectrónicas de próxima generación debido a su alta eficiencia de generación de segundo armónico (SHG) y su flexibilidad estructural. Aunque la MBANP ha sido caracterizada extensamente con respecto a su comportamiento dependiente de la temperatura y sus interacciones intermoleculares, su respuesta a la presión hidrostática permanece completamente inexplorada. Comprender esta respuesta es fundamental para dilucidar cómo las interacciones no covalentes gobiernan la estabilidad cristalina y la actividad NLO, y es prácticamente esencial para evaluar la estabilidad mecánica de estos materiales en dispositivos miniaturizados sometidos a estrés.

Metodología
Este estudio emplea un enfoque de primeros principios basado en la Teoría de la Funcional de la Densidad (DFT) utilizando el código CASTEP. Los cálculos utilizan la Aproximación de Gradiente Generalizado (GGA) con el funcional Perdew–Burke–Ernzerhof (PBE), aumentada por la corrección de van der Waals de Tkatchenko–Scheffler (TS) para modelar con precisión las interacciones intermoleculares débiles. Se utilizaron pseudopotenciales de conservación de norma y un conjunto de base de ondas planas con una energía de corte de 830 eV. El estudio investiga sistemáticamente monocristales de MBANP bajo presión hidrostática en un rango de 0 a 20 GPa. Los análisis clave incluyen:

• Evolución Estructural: Monitoreo de los parámetros de red y la conformación molecular.
• Propiedades Elásticas: Cálculo de constantes elásticas ($C_{ij}$) y derivación de módulos mecánicos (Bulk, Shear, Young, razón de Poisson) mediante el promedio de Voigt-Reuss-Hill, junto con verificaciones de los criterios de estabilidad de Born.
• Estructura Electrónica: Análisis de estructuras de bandas, densidad de estados (PDOS), densidad de carga, función de localización electrónica (ELF) y mecanismos de transferencia de carga.
• Interacciones Intermoleculares: Utilización de análisis de superficie de Hirshfeld y gráficos de huella digital 2D para cuantificar los tipos de contacto.
• Propiedades Vibracionales y Ópticas: Simulación de espectros Raman y cálculo de métricas ópticas dependientes de la frecuencia (función dieléctrica, índice de refracción, absorción, reflectividad).

Contribuciones Clave y Resultados

1. Transición de Fase Estructural a ~14 GPa:
   El estudio identifica una transición estructural isosimétrica inducida por presión cerca de los 14 GPa. Por debajo de este umbral, el cristal exhibe una compresibilidad anisotrópica pronunciada (eje $a$ > eje $c$ > eje $b$). A 14 GPa, una contracción anómala del eje $a$ es acompañada por una expansión discontinua de los ejes $b$ y $c$, y un salto en el ángulo $\beta$. Esto significa un reordenamiento molecular concertado o una distorsión de tipo cizalla, resultando en una fase de alta presión con propiedades mecánicas cuasi-isotrópicas. A pesar de estos cambios, la simetría del grupo espacial ($P2_1$) se preserva.

2. Estabilidad Mecánica y Endurecimiento:
   La MBANP permanece mecánicamente y dinámicamente estable a lo largo del rango de 0–20 GPa, confirmado por la satisfacción de los criterios de estabilidad de Born modificados y la ausencia de frecuencias imaginarias en los espectros de dispersión de fonones. El material experimenta un endurecimiento progresivo, con el Módulo de Bulk aumentando de 10.7 GPa a 105.2 GPa y el Módulo de Cizalla de 3.76 GPa a 27.2 GPa. La razón de Poisson ($\nu$) aumenta monotonicamente de 0.343 a 0.4 en el rango de 0–12 GPa antes de estabilizarse alrededor de 0.38, indicando una naturaleza dúctil. La relación $B/G$ exhibe una variación significativa, aumentando de 2.84 a 4.7 entre 0 y 12 GPa, y luego disminuyendo a 3.86 a 20 GPa, con un punto de inflexión a 14 GPa que correlaciona con la transición estructural.

3. Mejoras en la Estructura Electrónica y Transferencia de Carga:
   La presión aumenta significativamente la transferencia de carga intramolecular (CT) a lo largo del eje donante-aceptor (metilbencilamino al grupo nitro). La brecha de banda (bandgap) se estrecha de 2.226 eV a 0 GPa a 1.056 eV a 20 GPa, impulsada por el ensanchamiento de la banda de valencia (orbitales $\pi$ de la piridina) y la banda de conducción (estados $\pi^*$ del nitro). El análisis de densidad de carga revela un aumento de la depleción de electrones en el nitrógeno donante y la acumulación en el oxígeno aceptor, además de un acortamiento de la longitud del enlace de hidrógeno N–H⋯O, confirmando un efecto "push-pull" fortalecido.

4. Reorganización de las Interacciones Intermoleculares:
   El análisis de superficie de Hirshfeld cuantifica un cambio en la dominancia de las interacciones impulsado por la presión. Los contactos de van der Waals débiles e isotrópicos (H⋯H y C⋯H) disminuyen su contribución, mientras que las interacciones fuertes y direccionales aumentan. Notablemente, los contactos C⋯C (indicativos de apilamiento $\pi$–$\pi$) aumentan súbitamente, y los contactos específicos C–H⋯O se transforman en enlaces de hidrógeno genuinos (distancias < 2.6 Å) por encima de los 14 GPa. Esta reorganización se identifica como un mecanismo clave para mantener la estabilidad bajo compresión extrema.

5. Respuestas Vibracionales y Ópticas:
   Los espectros Raman exhiben un desplazamiento al azul (blue-shift) uniforme en todos los modos, atribuido a la compresión de los enlaces y al aumento de las constantes de fuerza. Una discontinuidad distinta en las intensidades de los picos y el ensanchamiento ocurre a 14 GPa, corroborando la transición de fase. Ópticamente, la constante dieléctrica estática y el índice de refracción aumentan con la presión (subiendo de ~3.56 a ~4.4 y de ~1.86 a ~2.1, respectivamente). Sin embargo, el borde de absorción y los picos principales de la función dieléctrica se desplazan hacia energías más altas (desplazamiento al azul), consistente con la conclusión de que el borde de absorción se desplaza a una energía mayor. Consecuentemente, la absorción del material en la región visible se debilita, lo que conduce a una disminución de la transmitancia (de ~87% a ~73% para una muestra de 100 nm) y un aumento de la reflectividad.

Significancia
Este trabajo proporciona la primera investigación teórica exhaustiva de la MBANP bajo alta presión, llenando un vacío crítico en la comprensión de sus relaciones estructura-propiedad. El estudio demuestra que la presión hidrostática actúa como una herramienta potente para ajustar las propiedades funcionales de los cristales orgánicos NLO sin alterar su simetría. Al elucidar los mecanismos de endurecimiento mecánico, mejora de la transferencia de carga y reorganización intermolecular, los hallazgos afirman la estabilidad de la MBANP para aplicaciones en dispositivos optoelectrónicos miniaturizados y ofrecen orientación para el diseño de cristales moleculares de próxima generación con características electromecánicas y ópticas ajustables.