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

Résumé technique : Transitions de phase structurelle induites par la pression dans le MBANP

Énoncé du problème
Les cristaux organiques non linéaires (NLO), spécifiquement le 2-(α-méthylbenzylamino)-5-nitropyridine (MBANP), sont essentiels pour les technologies photoniques et optoélectroniques de nouvelle génération en raison de leur haute efficacité de génération de seconde harmonique (SHG) et de leur flexibilité structurelle. Bien que le MBANP ait été largement caractérisé concernant son comportement dépendant de la température et ses interactions intermoléculaires, sa réponse à la pression hydrostatique demeure entièrement inexplorée. Comprendre cette réponse est fondamental pour élucider comment les interactions non covalentes régissent la stabilité cristalline et l'activité NLO, et est pratiquement essentiel pour évaluer la stabilité mécanique de ces matériaux dans les dispositifs miniaturisés soumis à des contraintes.

Méthodologie
Cette étude emploie une approche de premier principe basée sur la théorie de la fonctionnelle de la densité (DFT) utilisant le code CASTEP. Les calculs utilisent l'approximation du gradient généralisé (GGA) avec la fonctionnelle de Perdew–Burke–Ernzerhof (PBE), augmentée de la correction de van der Waals de Tkatchenko–Scheffler (TS) pour modéliser avec précision les interactions intermoléculaires faibles. Des pseudopotentiels conservant la norme et un ensemble de fonctions de base d'ondes planes avec une énergie de coupure de 830 eV ont été utilisés. L'étude étudie systématiquement les monocristaux de MBANP sous une pression hydrostatique allant de 0 à 20 GPa. Les analyses clés incluent :

• Évolution structurelle : Suivi des paramètres de maille et de la conformation moléculaire.
• Propriétés élastiques : Calcul des constantes élastiques ($C_{ij}$) et dérivation des modules mécaniques (module de compressibilité, de cisaillement, de Young, coefficient de Poisson) via la moyenne de Voigt-Reuss-Hill, ainsi que la vérification des critères de stabilité de Born.
• Structure électronique : Analyse des structures de bandes, de la densité d'états (PDOS), de la densité de charge, de la fonction de localisation électronique (ELF) et des mécanismes de transfert de charge.
• Interactions intermoléculaires : Utilisation de l'analyse de surface de Hirshfeld et des diagrammes d'empreintes digitales 2D pour quantifier les types de contacts.
• Propriétés vibrationnelles et optiques : Simulation des spectres Raman et calcul des métriques optiques dépendantes de la fréquence (fonction diélectrique, indice de réfraction, absorption, réflectivité).

Contributions clés et résultats

1. Transition de phase structurelle à ~14 GPa :
   L'étude identifie une transition structurelle isosymétrique induite par la pression près de 14 GPa. En dessous de ce seuil, le cristal présente une compressibilité anisotrope prononcée ($a$-axe > $c$-axe > $b$-axe). À 14 GPa, une contraction anormale de l'axe $a$ est accompagnée d'une expansion discontinue des axes $b$ et $c$ et d'un saut de l'angle $\beta$. Cela signifie un réarrangement moléculaire concerté ou une distorsion de type cisaillement, résultant en une phase à haute pression dotée de propriétés mécaniques quasi isotropes. Malgré ces changements, la symétrie du groupe d'espace ($P2_1$) est préservée.

2. Stabilité mécanique et durcissement :
   Le MBANP reste mécaniquement et dynamiquement stable sur toute la plage 0–20 GPa, confirmé par la satisfaction des critères de stabilité de Born modifiés et l'absence de fréquences imaginaires dans les spectres de dispersion de phonons. Le matériau subit un durcissement progressif, le module de compressibilité passant de 10,7 GPa à 105,2 GPa et le module de cisaillement de 3,76 GPa à 27,2 GPa. Le coefficient de Poisson ($\nu$) augmente de manière monotone de 0,343 à 0,4 dans la plage 0–12 GPa avant de se stabiliser autour de 0,38, indiquant une nature ductile. Le rapport $B/G$ présente une variation significative, augmentant de 2,84 à 4,7 entre 0 et 12 GPa, puis diminuant à 3,86 à 20 GPa, avec un point d'inflexion à 14 GPa corrélé à la transition structurelle.

3. Améliorations de la structure électronique et du transfert de charge :
   La pression améliore significativement le transfert de charge (CT) intramoléculaire le long de l'axe donneur-accepteur (méthylbenzylamino vers le groupe nitro). La bande interdite (bandgap) se rétrécit de 2,226 eV à 0 GPa à 1,056 eV à 20 GPa, pilotée par l'élargissement de la bande de valence (orbitales $\pi$ de la pyridine) et de la bande de conduction (états $\pi^*$ du nitro). L'analyse de la densité de charge révèle une augmentation de l'appauvrissement électronique au niveau de l'azote donneur et une accumulation au niveau de l'oxygène accepteur, ainsi qu'un raccourcissement de la liaison hydrogène N–H⋯O, confirmant un effet « push-pull » renforcé.

4. Réorganisation des interactions intermoléculaires :
   L'analyse de la surface de Hirshfeld quantifie un changement de dominance des interactions piloté par la pression. Les contacts de van der Waals faibles et isotropes (H⋯H et C⋯H) voient leur contribution diminuer, tandis que les interactions directionnelles fortes augmentent. Notamment, les contacts C⋯C (indicateurs d'un empilement $\pi$–$\pi$) augmentent brusquement, et des contacts C–H⋯O spécifiques se transforment en véritables liaisons hydrogène (distances < 2,6 Å) au-dessus de 14 GPa. Cette réorganisation est identifiée comme un mécanisme clé pour maintenir la stabilité sous compression extrême.

5. Réponses vibrationnelles et optiques :
   Les spectres Raman présentent un déplacement vers le bleu (blue-shift) uniforme pour tous les modes, attribué à la compression des liaisons et à l'augmentation des constantes de force. Une discontinuité distincte dans les intensités de pics et l'élargissement se produit à 14 GPa, corroborant la transition de phase. Optiquement, la constante diélectrique statique et l'indice de réfraction augmentent avec la pression (passant d'environ 3,56 à ~4,4 et de ~1,86 à ~2,1, respectivement). Cependant, le bord d'absorption et les pics principaux de la fonction diélectrique se déplacent vers des énergies plus élevées (déplacement vers le bleu), ce qui est cohérent avec la conclusion que le bord d'absorption se déplace vers des énergies plus hautes. Par conséquent, l'absorption du matériau dans la région visible s'affaiblit, entraînant une diminution de la transmittance (de ~87 % à ~73 % pour un échantillon de 100 nm) et une augmentation de la réflectivité.

Signification
Ce travail constitue la première investigation théorique complète du MBANP sous haute pression, comblant une lacune critique dans la compréhension de ses relations structure-propriété. L'étude démontre que la pression hydrostatique est un outil puissant pour ajuster les propriétés fonctionnelles des cristaux NLO organiques sans altérer leur symétrie. En élucidant les mécanismes de durcissement mécanique, d'amélioration du transfert de charge et de réorganisation intermoléculaire, les résultats confirment la stabilité du MBANP pour des applications dans les disposités optoélectroniques miniaturisés et offrent des orientations pour la conception de cristaux moléculaires de prochaine génération dotés de caractéristiques électromécaniques et optiques ajustables.