Data analysis methods for powder x-ray diffraction intensity under laser-driven dynamic compression at Omega and NIF laser facilities

Sintesi Tecnica: Metodi di Analisi dei Dati per l'Intensità della Diffrazione di Raggi X su Polvere sotto Compressione Dinamica Indotta da Laser

Definizione del Problema
La diffrazione di raggi X su polvere (PXRD) sotto compressione dinamica indotta da laser presso strutture come Omega e la National Ignition Facility (NIF) è uno strumento critico per investigare il comportamento dei materiali sotto pressione, temperatura e velocità di deformazione estreme. Sebbene piattaforme robuste come la Powder X-Ray Diffraction Image Plate (PXRDIP) presso Omega e la TARDIS (TARget Diffraction In Situ) presso NIF abbiano documentato con successo la presenza di picchi di Bragg e inferito la densità, determinare la fedeltà dell'intensità PXRD rimane una sfida significativa.

Le difficoltà primarie includono:

1. Geometria Complessa: Determinare accuratamente le posizioni relative della sorgente di raggi X, del campione, dei materiali di riferimento e dei pannelli delle lastre di immagine (image plate), considerando anche l'orientamento non standard dei detector.
2. Rumore di Fondo: I dati grezzi sono spesso dominati da background spazialmente variabili generati dall'ablazione laser ad alta intensità e dalle sorgenti di plasma a raggi X.
3. Fedeltà del Segnale: Gli studi precedenti si sono concentrati ampiamente sulle posizioni dei picchi; tuttavia, vincolare accuratamente l'intensità XRD è necessario per studiare le trasformazioni di fase, la coesistenza di fasi e transizioni sottili (ad esempio, verso il comportamento superionico).
4. Effetti Termici e di Densità: Confrontare la PXRD di materiali compressi per shock con materiali di riferimento ambientali richiede di tenere conto dello smorzamento termico (effetti Debye-Waller) e delle variazioni di densità.

Metodologia
Gli autori presentano un workflow di analisi dati raffinato, implementato nel pacchetto ANALYZEPXRDIP (scritto in IGORPRO), per migliorare la determinazione dell'intensità. La metodologia è illustrata utilizzando dati di diamante sottoposto a shock raccolti vicino a 700 GPa presso la struttura laser Omega EP. Le fasi chiave includono:

• Calibrazione Geometrica: Viene utilizzato un approccio di calibrazione iterativo per determinare le posizioni precise della sorgente di raggi X, del pacchetto campione e dei pannelli delle lastre di immagine. Ciò utilizza il segnale di diffrazione del foro stenopeico (pinhole) di collimazione (o uno strato di riferimento non compresso) come riferimento in situ, anziché affidarsi a standard esterni come la polvere di $\text{CeO}_2$.
• Proiezione dei Dati: I dati grezzi della lastra di immagine vengono proiettati in array 2D di angolo di diffrazione ($2\theta$) rispetto all'angolo azimutale ($\phi$) e al raggio di scattering ($R_{ZA}$), tenendo conto della specifica geometria dove il volume di scattering varia attraverso il detector.
• Correzioni di Intensità (Pre-Sottrazione del Background):
  • Area Effettiva del Pixel: Correzione per la distanza non uniforme campione-detector e l'angolo di incidenza sulle lastre di immagine.
  • Area del Campione Rilevata: Scalatura dell'intensità in base al volume di scattering effettivo del campione visibile a ciascun punto del detector, definito dalla geometria del pinhole e dall'angolo di attacco.
  • Modellazione del Segnale del Pinhole: Stima dell'intensità XRD del pinhole basata sul volume di materiale del pinhole che contribuisce alla diffrazione in ogni punto del detector.
• Sottrazione del Background: Un algoritmo 2D SNIP (Statistics-sensitive Non-linear Iterative Peak-clipping) con kernel a $2\theta$ costante viene applicato per rimuovere il background intenso generato dal drive laser e dal plasma.
• Correzioni Post-Sottrazione:
  • Sensibilità del Detector: Correzione per la sensibilità dipendente dall'energia delle lastre di immagine e dall'angolo di incidenza.
  • Attenuazione: Compensazione dell'attenuazione dei raggi X attraverso i filtri, il pacchetto campione (ablatore, scudi, finestre) e lo strato di copertura della lastra di immagine.
• Scalatura dell'Intensità e Analisi:
  • Scalatura del Volume: Scalatura delle aree dei picchi tramite lo spessore effettivo delle frazioni di campione shockate ($d_S$) e non shockate ($d_U$), derivate da dati VISAR o simulazioni idrodinamiche.
  • Scalatura della Luminosità della Sorgente: Normalizzazione dell'intensità tramite la luminosità della sorgente di raggi X, stimata tramite l'energia del laser e l'efficienza di conversione o misurata tramite spettrometri.
  • Scalatura del Riferimento: Utilizzo della porzione non shockata del campione o del materiale del pinhole come riferimento interno.
  • Correzione dello Smorzamento Termico: Applicazione di un approccio Debye-Waller per tenere conto della riduzione dell'intensità XRD dovuta alle alte temperature e densità del campione shockato rispetto al riferimento ambiente. Ciò comporta il calcolo dello spostamento interatomico medio al quadrato ($u^2$) utilizzando un modello Debye.

Risultati Chiave
Il paper applica questo workflow a un esperimento di diamante sottoposto a shock vicino a 700 GPa (circa 1 TPa in studi correlati).

• Il workflow separa con successo i segnali del diamante shockato, del diamante non shockato e del pinhole di platino.
• Applicando la correzione Debye-Waller e la scalatura del volume, gli autori hanno calcolato la frazione solida ($SF$) del campione shockato.
• Per il dataset specifico analizzato, la stima sperimentale ha fornito una frazione solida di $SF = 0,6 \pm 0,05$.
• Gli autori osservano che i modelli teorici che assumono uno stato di Hugoniot idrostatico suggeriscono una frazione solida più vicina a 0,9, evidenziando la difficoltà nel raggiungere l'accuratezza assoluta dell'intensità a queste strutture.

Significato e Rivendicazioni
Il documento sostiene che questi metodi raffinati offrono "opportunità promettenti" per andare oltre la semplice identificazione dei picchi verso l'analisi quantitativa delle trasformazioni di fase e della coesistenza. Nello specifico:

• Miglioramento della Fedeltà: I metodi consentono una migliore costrizione dell'intensità XRD, essenziale per identificare cambiamenti di fase sottili, come la transizione al comportamento superionico.
• Referenziazione In-situ: L'uso del pinhole o degli strati di campione non shockati come riferimenti in situ fornisce un modo robusto per calibrare geometria e intensità senza standard esterni.
• Trasferibilità: Sebbene focalizzato sulla piattaforma Omega PXRDIP, gli autori affermano che i metodi sono direttamente applicabili alla piattaforma TARDIS presso NIF e sono trasferibili al Laser Fusion Research Center (LFRC) e altre strutture come LULI2000 e ORION con modifiche minori.
• Modestia sulle Limitazioni: Gli autori riconoscono esplicitamente le limitazioni, tra cui il potenziale bias dell'algoritmo SNIP contro i segnali ampi, l'omissione degli effetti di orientazione preferita (texture), l'incapacità di quantificare lo scattering termico diffuso a causa della rimozione del background e la semplificazione dell'uso di una singola temperatura di Debye per stati fononici complessi. Concludono che, sebbene il framework migliori la determinazione dell'intensità, le discrepanze tra le stime sperimentali e la teoria (come visto nel calcolo della frazione solida) indicano che è necessaria ulteriore evoluzione per una quantificazione assoluta.