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

Technisches Resümee: Druckinduzierte strukturelle Phasenübergänge in MBANP

Problemstellung
Organische nichtlineare optische (NLO) Kristalle, speziell 2-(α-Methylbenzylamino)-5-nitropyridin (MBANP), sind aufgrund ihrer hohen Frequenzverdopplungseffizienz (Second-Harmonic Generation, SHG) und ihrer strukturellen Flexibilität entscheidend für die nächste Generation photonischer und optoelektronischer Technologien. Während MBANP hinsichtlich seines temperaturabhängigen Verhaltens und seiner intermolekularen Wechselwirkungen bereits umfassend charakterisiert wurde, bleibt seine Reaktion auf hydrostatischen Druck bislang gänzlich unerforscht. Das Verständnis dieser Reaktion ist grundlegend für die Aufklärung darüber, wie nichtkovalente Wechselwirkungen die Kristallstabilität und die NLO-Aktivität steuern, und praktisch essenziell für die Bewertung der mechanischen Stabilität dieser Materialien in miniaturisierten Bauteilen unter mechanischer Belastung.

Methodik
Diese Studie verwendet einen First-Principles-Ansatz basierend auf der Dichtefunktionaltheorie (DFT) unter Verwendung des CASTEP-Codes. Die Berechnungen nutzen die verallgemeinerte Gradienten-Approximation (GGA) mit dem Perdew–Burke–Ernzerhof-Funktional (PBE), ergänzt durch die Tkatchenko–Scheffler (TS) Van-der-Waals-Korrektur, um schwache intermolekulare Wechselwirkungen präzise zu modellieren. Normerhaltende Pseudopotentiale und ein Plane-Wave-Basissatz mit einer Cutoff-Energie von 830 eV wurden eingesetzt. Die Studie untersucht systematisch MBANP-Einkristalle unter hydrostatischem Druck im Bereich von 0 bis 20 GPa. Zu den zentralen Analysen gehören:

• Strukturelle Entwicklung: Überwachung der Gitterparameter und der molekularen Konformation.
• Elastische Eigenschaften: Berechnung der elastischen Konstanten ($C_{ij}$) und Ableitung der mechanischen Moduli (Bulk-, Scher-, Youngs-Modul, Poisson-Zahl) mittels Voigt-Reuss-Hill-Mittelung, ergänzt durch Überprüfungen der Born-Stabilitätskriterien.
• Elektronische Struktur: Analyse der Bandstrukturen, der Zustandsdichte (PDOS), der Ladungsdichte, der Elektronenlokalisationsfunktion (ELF) und der Ladungstransfermechanismen.
• Intermolekulare Wechselwirkungen: Nutzung der Hirshfeld-Oberflächenanalyse und 2D-Fingerprint-Plots zur Quantifizierung der Kontaktarten.
• Vibrations- und optische Eigenschaften: Simulation von Raman-Spektren und Berechnung frequenzabhängiger optischer Metriken (Dielektrizitätsfunktion, Brechungsindex, Absorption, Reflexion).

Wesentliche Beiträge und Ergebnisse

1. Struktureller Phasenübergang bei ~14 GPa:
   Die Studie identifiziert einen druckinduzierten isosymmetrischen strukturellen Übergang nahe 14 GPa. Unterhalb dieser Schwelle weist der Kristall eine ausgeprägte anisotrope Kompressibilität auf ($a$-Achse > $c$-Achse > $b$-Achse). Bei 14 GPa geht eine anomale Kontraktion der $a$-Achse mit einer diskontinuierlichen Expansion der $b$- und $c$-Achsen sowie einem Sprung des $\beta$-Winkels einher. Dies deutet auf eine koordinierte molekulare Umlagerung oder eine scherähnliche Verzerrung hin, was zu einer Hochdruckphase mit quasi-isotropen mechanischen Eigenschaften führt. Trotz dieser Änderungen bleibt die Raumgruppensymmetrie ($P2_1$) erhalten.

2. Mechanische Stabilität und Versteifung:
   MBANP bleibt im Bereich von 0–20 GPa mechanisch und dynamisch stabil, was durch die Erfüllung modifizierter Born-Stabilitätskriterien und das Fehlen imaginärer Frequenzen in den Phonon-Dispersionsspektren bestätigt wird. Das Material erfährt eine progressive Versteifung, wobei das Kompressionsmodul (Bulk Modulus) von 10,7 GPa auf 105,2 GPa und der Schubmodul (Shear Modulus) von 3,76 GPa auf 27,2 GPa ansteigt. Die Poisson-Zahl ($\nu$) steigt im Bereich von 0–12 GPa monoton von 0,343 auf 0,4 an, bevor sie sich bei etwa 0,38 stabilisiert, was auf eine duktile Natur hindeutet. Das $B/G$-Verhältnis zeigt erhebliche Variationen und steigt von 2,84 auf 4,7 zwischen 0 und 12 GPa an, um dann auf 3,86 bei 20 GPa zu sinken, wobei ein Wendepunkt bei 14 GPa mit dem strukturellen Übergang korreliert.

3. Verbesserung des elektronischen Zustands und des Ladungstransfers:
   Der Druck verstärkt den intramolekularen Ladungstransfer (CT) entlang der Donor-Akzeptor-Achse (Methylbenzylamino zu der Nitrogruppe) signifikant. Die Bandlücke verengt sich von 2,226 eV bei 0 GPa auf 1,056 eV bei 20 GPa, getrieben durch die Verbreiterung des Valenzbandes (Pyridin-$\pi$-Orbitale) und des Leitungsbandes (Nitro-$\pi^*$-Zustände). Die Ladungsdichteanalyse zeigt eine erhöhte Elektronenverarmung am Donor-Stickstoff und eine Akkumulation am Akzeptor-Sauerstoff sowie eine verkürzte N–H$\cdots$O Wasserstoffbrückenbindung, was den verstärkten „Push-Pull“-Effekt bestätigt.

4. Reorganisation der intermolekularen Wechselwirkungen:
   Die Hirshfeld-Oberflächenanalyse quantifiziert eine druckgetriebene Verschiebung der Dominanz der Wechselwirkungen. Schwache, isotrope Van-der-Waals-Kontakte (H$\cdots$H und C$\cdots$H) nehmen an Beitrag ab, während starke, gerichtete Wechselwirkungen zunehmen. Insbesondere steigen die C$\cdots$C-Kontakte (Indikator für $\pi$–$\pi$-Stapelung) stark an, und spezifische C–H$\cdots$O-Kontakte transformieren sich oberhalb von 14 GPa in echte Wasserstoffbrückenbindungen (Abstände < 2,6 Å). Diese Reorganisation wird als Schlüsselmechanismus zur Aufrechterhaltung der Stabilität unter extremer Kompression identifiziert.

5. Vibrations- und optische Antworten:
   Die Raman-Spektren zeigen eine einheitliche Blauverschiebung über alle Moden hinweg, was auf die Bindungskompression und erhöhte Kraftkonstanten zurückzuführen ist. Eine deutliche Diskontinuität in den Peak-Intensitäten und eine Verbreiterung treten bei 14 GPa auf, was den Phasenübergang untermauert. Optisch steigen die statische Dielektrizitätskonstante und der Brechungsindex mit dem Druck (Anstieg von ~3,56 auf ~4,4 bzw. ~1,86 auf ~2,1). Jedoch verschieben sich die Absorptionskante und die Hauptmaxima der Dielektrizitätsfunktion zu höheren Energien (Blauverschiebung), was konsistent mit der Schlussfolgerung ist, dass die Absorptionskante zu höheren Energien wandert. Infolgedessen schwächt sich die Absorption im sichtbaren Bereich ab, was zu einer Abnahme der Transmission (von ~87 % auf ~73 % bei einer 100 nm Probe) und einer Zunahme der Reflexion führt.

Bedeutung
Diese Arbeit liefert die erste umfassende theoretische Untersuchung von MBANP unter hohem Druck und schließt eine kritische Lücke im Verständnis der Struktur-Eigenschafts-Beziehungen. Die Studie zeigt, dass hydrostatischer Druck ein wirkungsvolles Werkzeug zur Abstimmung der funktionellen Eigenschaften organischer NLO-Kristalle ist, ohne deren Symmetrie zu verändern. Durch die Aufklärung der Mechanismen der mechanischen Versteifung, der Verstärkung des Ladungstransfers und der intermolekularen Reorganisation bestätigen die Ergebnisse die Stabilität von MBANP für Anwendungen in miniaturisierten optoelektronischen Bauteilen und bieten Orientierung für das Design der nächsten Generation molekularer Kristalle mit abstimmbaren elektromechanischen und optischen Eigenschaften.