Recovering full crystallographic orientation in specimen coordinates from polarised μ-FTIR spectra

Sintesi Tecnica: Recupero dell'orientamento cristallografico completo da spettri $\mu$-FTIR polarizzati

Il Problema
La spettroscopia Fourier-transform infrared (FTIR) è il metodo standard per quantificare i tracce di ossidrile (OH) nei minerali nominalmente anidri (NAM) grazie alla sua accessibilità e velocità. Tuttavia, i cristalli birifrangenti assorbono la luce infrarossa in modo anisotropo, il che significa che l'assorbanza misurata dipende dall'orientamento del cristallo rispetto al vettore di polarizzazione. Sebbene la FTIR polarizzata possa teoricamente risolvere sia la concentrazione che l'orientamento dei difetti contenenti OH, ciò richiede la conoscenza preventiva dell'orientamento del cristallo rispetto allo strumento.

Le pratiche attuali affrontano limitazioni significative:

1. FTIR non polarizzata: La maggior parte degli studi utilizza spettri non polarizzati su cristalli orientati casualmente. Questo approccio scarta le informazioni sull'orientamento e sottostima sistematicamente le concentrazioni di OH, a meno che non vengano misurati centinaia di grani, il che è spesso impraticabile.
2. Pre-orientamento: Misurare spettri polarizzati su cristalli pre-orientati richiede l'allineamento del campione prima della misurazione, un ostacolo per i materiali a grana fine.
3. Accoppiamento EBSD: L'associazione della FTIR con la diffrazione elettronica di backscatter (EBSD) risolve il problema dell'orientamento, ma richiede due strumenti, introduce errori di allineamento (tipicamente $\pm 8^\circ$) durante il trasferimento dei campioni e media i dati di orientamento invece di mantenerli internamente coerenti con lo spettro.
4. Metodi spettrali esistenti: I precedenti tentativi di recuperare l'orientamento direttamente dagli spettri (ad es. Asimow et al., 2006) presentano limitazioni: non collocano il cristallo nel sistema di riferimento del campione (rendendo difficile il confronto con altre tecniche), si affidano a intervalli spettrali fissi che possono soffrire di non linearità del rivelatore e mancano di generalità per altre tecniche polarizzate.

Metodologia
Gli autori presentano un metodo per recuperare l'orientamento cristallografico completo di un cristallo anisotropo in coordinate di campione, espresso tramite angoli di Eulero, utilizzando spettri di $\mu$-FTIR polarizzati. L'innovazione centrale è il "metodo della sezione di lunghezza d'onda", che ricostruisce una sezione bidimensionale della risposta anisotropa del cristallo.

• Quadro Teorico: Il metodo si basa sull'equazione di trasmittanza (Asimow et al., 2006), che mette in relazione la trasmittanza con gli assi principali del cristallo, lo spessore del campione e l'orientamento del vettore di polarizzazione. A differenza degli approcci precedenti che risolvono l'orientamento del vettore di polarizzazione rispetto al cristallo, questo metodo risolve l'orientamento del cristallo rispetto a un sistema di riferimento esterno fisso (microscopio).
• Sistemi di Riferimento: Gli autori definiscono tre sistemi di riferimento accoppiati: il sistema del campione, il sistema del microscopio (dove il polarizzatore è fisso lungo l'asse x) e il sistema del cristallo (definito dalla cella unitaria). L'orientamento viene determinato trovando la rotazione (utilizzando gli angoli di Eulero Bunge $zxz$) che allinea i dati di trasmittanza misurati con un "involucro di trasmittanza" derivato da spettri di riferimento lungo i tre assi principali.
• Acquisizione dei Dati: Il metodo richiede spettri polarizzati raccolti su un intervallo di angoli del polarizzatore ($\alpha$) a una singola lunghezza d'onda diagnostica. Ogni coppia $(T, \alpha)$ viene convertita in un punto cartesiano sull'involucro di trasmittanza.
• Ottimizzazione: L'algoritmo minimizza il disaccordo tra i punti misurati e l'involucro teorico applicando iterativamente rotazioni di prova. Per gestire la natura non convessa della superficie di costo dell'orientamento (molti minimi locali), gli autori hanno testato e raccomandano l'uso di un ottimizzatore basato sul gradiente con molteplici tentativi iniziali casuali o un ottimizzatore di evoluzione differenziale globale.
• Selezione della Lunghezza d'Onda: Il metodo richiede una singola lunghezza d'onda in cui l'involucro di trasmittanza sia sufficientemente anisotropo (le trasmittanze principali differiscono) e in cui lo spessore del campione mantenga la trasmittanza nell'intervallo di risposta lineare del rivelatore (tipicamente 10–90% T). Gli autori dimostrano questo utilizzando la regione dell'overtone della silice (1500–2100 cm⁻¹) nell'olivina.

Contributi Chiave

1. Recupero dell'Orientamento Completo: Il metodo recupera l'orientamento cristallografico completo (angoli di Eulero) in coordinate di campione, non solo l'orientamento del vettore di polarizzazione rispetto agli assi del cristallo. Ciò consente il confronto diretto con l'EBSD e altri dati strutturali.
2. Flessibilità della Singola Lunghezza d'Onda: Operando a una singola lunghezza d'onda, il metodo evita la necessità di ampi intervalli spettrali, consentendo agli utenti di scegliere numeri d'onda che evitino la non linearità del rivelatore o problemi di linea di base.
3. Implementazione Open-Source: Gli autori forniscono FTIRkit, uno strumento Python open-source che implementa il metodo, basandosi su librerie scientifiche standard (NumPy, Pandas, SciPy, Matplotlib).
4. Validazione: Il metodo è validato attraverso simulazioni Monte Carlo e dati sperimentali da cristalli di olivina a orientamento noto (misurato tramite EBSD).

Risultati

• Prestazioni della Simulazione: Nelle simulazioni prive di rumore, il metodo raggiunge un'elevata accuratezza. Un ottimizzatore basato sul gradiente con 20 tentativi iniziali recupera gli orientamenti entro 1° dal valore reale in >99,75% dei casi. L'accuratezza angolare mediana è di circa 0,004°.
• Sensibilità al Rumore: Quando viene introdotto il rumore gaussiano, l'accuratezza migliora monotonicamente con il numero di misurazioni (dimensione del campione). Il metodo è più sensibile al rumore di trasmittanza che al rumore angolare ($\alpha$). Ridurre il rumore di trasmittanza migliora significativamente il recupero dell'orientamento.
• Validazione Sperimentale (Olivina): Il metodo è stato testato su grani di olivina da uno xenolite di Kilbourne Hole a due spessori (447,5 µm e 244 µm).
  • Per il campione più sottile, il 70% dei grani (16/23) è stato recuperato entro una soglia di disorientamento di 25° rispetto all'EBSD.
  • Per il campione più spesso, il 50% dei grani (15/30) ha soddisfatto la soglia.
  • Il successo dipende dalla specifica combinazione di numero d'onda e spessore del campione per garantire che la trasmittanza rimanga nell'intervallo di risposta lineare del rivelatore e che la forma dell'involucro sia distinta.
  • I fallimenti sono stati spesso attribuiti a grani in cui lo spettro misurato non poteva essere riprodotto dal modello teorico (indicando potenziali problemi con lo standard di riferimento, inclusioni nel campione o variazioni di spessore) piuttosto che un fallimento dell'algoritmo stesso.

Significato e Rivendicazioni
Gli autori pongono questo lavoro come una soluzione al "problema dell'orientamento" nella spettroscopia polarizzata. Affermano che, sfruttando l'anisotropia del segnale FTIR stesso, il metodo elimina la necessità di un secondo strumento (EBSD) e mantiene la coerenza interna tra lo spettro e l'orientamento.

Il documento afferma modestamente che l'approccio:

• Si applica a qualsiasi fase cristallina birifrangente, purché siano disponibili spettri di riferimento delle direzioni principali e un modello che relazioni lo spettro all'orientamento.
• Può essere esteso ad altre tecniche polarizzate (Raman, UV–Vis, IR diretta laser).
• Offre un'alternativa più flessibile agli approcci basati sull'intervallo spettrale, permettendo all'utente di evitare numeri d'onda degradati.

Gli autori non pretendono che il metodo sia perfetto per tutte le condizioni sperimentali; notano esplicitamente che il successo dipende dalla selezione di lunghezze d'onda e spessori appropriati per mantenere la risposta lineare del rivelatore e che le sfide sperimenti (movimento del campione, deformazione dell'epossidica) possono influenzare l'accuratezza dei confronti reali con l'EBSD. Il contributo primario è la dimostrazione che l'orientamento completo può essere recuperato dagli spettri polarizzati da soli, rimuovendo una barriera importante per l'analisi di routine della $\mu$-FTIR polarizzata.