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

Sintesi Tecnica: Transizioni di Fase Strutturale Indotte dalla Pressione in MBANP

Definizione del Problema
I cristalli organici non lineari ottici (NLO), specificamente il 2-(α-metilbenzilamino)-5-nitropiridina (MBANP), sono critici per le tecnologie fotoniche e optoelettroniche di prossima generazione grazie alla loro elevata efficienza di generazione di seconda armonica (SHG) e alla loro flessibilità strutturale. Sebbene l'MBANP sia stato ampiamente caratterizzato per quanto riguarda il suo comportamento dipendente dalla temperatura e le interazioni intermolecolari, la sua risposta alla pressione idrostatica rimane del tutto inesplorata. Comprendere questa risposta è fondamentale per chiarire come le interazioni non covalenti governino la stabilità dei cristalli e l'attività NLO, e risulta praticamente essenziale per valutare la stabilità meccanica di questi materiali in dispositivi miniaturizzati soggetti a stress.

Metodologia
Questo studio impiega un approccio basato sul primo principio tramite la Teoria del Funzionale della Densità (DFT) utilizzando il codice CA-STEP. I calcoli utilizzano l'Approssimazione del Gradiente Generalizzato (GGA) con il funzionale Perdew–Burke–Ernzerhof (PBE), aumentata dalla correzione di van der Waals Tkatchenko–Scheffler (TS) per modellare accuratamente le deboli interazioni intermolecolari. Sono stati utilizzati pseudopotenziali normo-conservativi e un set di basi a onde piane con un'energia di cutoff di 830 eV. Lo studio investiga sistematicamente i monocristalli di MBANP sotto pressione idrostatica compresa tra 0 e 20 GPa. Le analisi chiave includono:

• Evoluzione Strutturale: Monitoraggio dei parametri reticolari e della conformazione molecolare.
• Proprietà Elastiche: Calcolo delle costanti elastiche ($C_{ij}$) e derivazione dei moduli meccanici (Modulo di Bulk, di Taglio, di Young, Rapporto di Poisson) tramite la media di Voigt-Reuss-Hill, insieme ai controlli dei criteri di stabilità di Born.
• Struttura Elettronica: Analisi delle strutture a bande, della densità degli stati (PDOS), della densità di carica, della funzione di localizzazione elettronica (ELF) e dei meccanismi di trasferimento di carica.
• Interazioni Intermolecolari: Utilizzo dell'analisi della superficie di Hirshfeld e dei grafici fingerprint 2D per quantificare i tipi di contatto.
• Proprietà Vibrazionali e Ottiche: Simulazione degli spettri Raman e calcolo delle metriche ottiche dipendenti dalla frequenza (funzione dielettrica, indice di rifrazione, assorbimento, riflettività).

Contributi Chiave e Risultati

1. Transizione di Fase Strutturale a ~14 GPa:
   Lo studio identifica una transizione strutturale isosimmetrica indotta dalla pressione vicino a 14 GPa. Al di sotto di questa soglia, il cristallo esibisce una pronunciata compressibilità anisotropa ($a$-asse > $c$-asse > $b$-asse). A 14 GPa, una contrazione anomala dell' $a$-asse è accompagnata da un'espansione discontinua dei $b$- e $c$-assi e da un salto nell'angolo $\beta$. Ciò segnala un riarrangiamento molecolare concertato o una distorsione di tipo shear, che risulta in una fase ad alta pressione con proprietà meccaniche quasi isotrope. Nonostante questi cambiamenti, la simmetria del gruppo spaziale ($P2_1$) viene preservata.

2. Stabilità Meccanica e Irrigidimento:
   L'MBANP rimane meccanicamente e dinamicamente stabile nell'intervallo 0–20 GPa, come confermato dal soddisfacimento dei criteri di stabilità modificati di Born e dall'assenza di frequenze immaginarie negli spettri di dispersione fononica. Il materiale subisce un progressivo irrigidimento, con il Modulo di Bulk che aumenta da 10,7 GPa a 105,2 GPa e il Modulo di Taglio da 3,76 GPa a 27,2 GPa. Il rapporto di Poisson ($\nu$) aumenta monotonicamente da 0,343 a 0,4 nell'intervallo 0–12 GPa prima di stabilizzarsi intorno a 0,38, indicando una natura duttile. Il rapporto $B/G$ mostra una variazione significativa, aumentando da 2,84 a 4,7 tra 0 e 12 GPa, per poi diminuire a 3,86 a 20 GPa, con un punto di inflessione a 14 GPa correlato alla transizione strutturale.

3. Potenziamento del Trasferimento di Carica e della Struttura Elettronica:
   La pressione potenzia significativamente il trasferimento di carica intramolecolare (CT) lungo l'asse donatore-accettore (dal gruppo metilbenzilamminico al gruppo nitro). Il bandgap si restringe da 2,226 eV a 0 GPa a 1,056 eV a 20 GPa, guidato dall'allargamento della banda di valenza (orbitali $\pi$ della piridina) e della banda di conduzione (stati $\pi^*$ del nitro). L'analisi della densità di carica rivela un aumento della deplezione elettronica al azoto donatore e l'accumulo all'ossigeno accettore, insieme a un accorciamento della lunghezza del legame idrogeno N–H⋯O, confermando un rafforzato effetto "push-pull".

4. Riorganizzazione delle Interazioni Intermolecolari:
   L'analisi della superficie di Hirshfeld quantifica uno spostamento della dominanza delle interazioni guidato dalla pressione. I deboli contatti di van der Waals isotropi (H⋯H e C⋯H) diminuiscono nel contributo, mentre aumentano le interazioni direzionali forti. In particolare, i contatti C⋯C (indicativi dello stacking $\pi$–$\pi$) aumentano drasticamente, e specifici contatti C–H⋯O si trasformano in veri legami a idrogeno (distanze < 2,6 Å) sopra i 14 GPa. Questa riorganizzazione è identificata come un meccanismo chiave per mantenere la stabilità sotto compressione estrema.

5. Risposte Vibrazionali e Ottiche:
   Gli spettri Raman mostrano un blue-shift uniforme in tutti i modi, attribuito alla compressione dei legami e all'aumento delle costanti di forza. Si verifica una distinta discontinuità nelle intensità dei picchi e un allargamento a 14 GPa, corroborando la transizione di fase. Otticamente, la costante dielettrica statica e l'indice di rifrazione aumentano con la pressione (salendo da ~3,56 a ~4,4 e da ~1,86 a ~2,1, rispettivamente). Tuttavia, il bordo di assorbimento e i picchi principali della funzione dielettrica si spostano verso energie più elevate (blue-shift), coerentemente con la conclusione che il bordo di assorbimento si sposta verso energie più alte. Di conseguenza, l'assorbimento del materiale nella regione del visibile si indebolisce, portando a una diminuzione della trasmittanza (da ~87% a ~73% per un campione da 100 nm) e un aumento della riflettività.

Significato
Questo lavoro fornisce la prima investigazione teorica completa dell'MBANP ad alta pressione, colmando una lacuna critica nella comprensione delle relazioni struttura-proprietà. Lo studio dimostra che la pressione idrostatica agisce come uno strumento potente per modulare le proprietà funzionali dei cristalli NLO organici senza alterarne la simmetria. Elucidando i meccanismi di irrigidimento meccanico, potenziamento del trasferimento di carica e riorganizzazione intermolecolare, i risultati confermano la stabilità dell'MBANP per applicazioni in dispositivi optoelettronici miniaturizzati e offrono guida per la progettazione di cristalli molecolari di prossima generazione con caratteristiche elettromeccaniche e ottiche modulabili.