Role of Local Structural Variation in X-ray Photoelectron Spectrum of Silicon Oxide Interfaces

Sintesi Tecnica: Il Ruolo della Variazione Strutturale Locale negli Spettri di Fotoemissione a Raggi X dell'Interfaccia Silicio-Ossido

Problematica
L'interpretazione dei dati di spettroscopia fotoelettronica a raggi X (XPS) per le interfacce silicio-ossido (SiO$_x$) su substrati di silicio si è tradizionalmente basata su un modello "basato su motivi" (motif-based). In questa visione convenzionale, i larghi spettri Si 2p (e, per estensione, Si 2s) osservati negli ossidi sottostequiometrici vengono deconvoluti in componenti discrete corrispondenti a stati di ossidazione interi specifici (Si$^{1+}$, Si$^{2+}$, Si$^{3+}$, Si$^{4+}$). Queste componenti sono spesso interpretate come rappresentazioni di distinti strati di ossido intermedi o di specifici motivi di legame locale (ad esempio, atomi di Si legati a 1, 2 o 3 atomi di ossigeno). Tuttavia, le proprietà fisiche dei film sottili di SiO$_x$ sono meno caratterizzate rispetto al SiO$_2$ bulk, e l'affidabilità dell'assegnazione degli stati chimici basata su singole componenti di livello core è sempre più messa in discussione, particolarmente quando i film di ossido si assottigliano e i contributi sottostequiometrici diventano dominanti.

Metodologia
Gli autori hanno impiegato un approccio duale combinando la ricostruzione sperimentale avanzata con simulazioni statistiche ad alte prestazioni:

1. Ricostruzione Sperimentale:

   • Campione: Un wafer di Si(100) commerciale con uno strato di SiO$_2$ cresciuto termicamente di 40 nm.
   • Processo: I dati XPS sono stati raccolti ai bordi O 1s e Si 2s mentre la superficie veniva progressivamente bombardata da un fascio di Ar$^+$ da 300 eV.
   • Algoritmo di Ricostruzione: Gli spettri grezzi sono stati trattati come segnali integrati in profondità e attenuati dallo scattering anelastico degli elettroni. Utilizzando un metodo di inversione di matrice ($S^{(obs)} = AS$), gli autori hanno ricostruito gli spettri risolti per strato ($s_j$) per singole profondità, evitando i problemi di splitting spin-orbita concentrandosi sui bordi O 1s e Si 2s.
   • Parametri: È stato utilizzato un tasso di erosione di 0,012 nm/s e un cammino libero medio degli elettroni inelastici ($\lambda_{el}$) di 3,5 nm.

2. Simulazioni Statistiche:

   • Modellazione: Sono state eseguite simulazioni di Dinamica Molecolare (MD) su celle periodiche (~15 Å) di SiO$_x$ con composizioni $x = 0, 1, 0,5, 1,0, 1,5, 2,0$.
   • Potenziale: È stato utilizzato il potenziale di machine learning MACE per le corse MD, che hanno incluso la fusione a 3000 K e l'annealing a 300 K per generare diversi ensemble strutturali.
   • Calcolo dell'Energia di Legame: Le energie di legame (BE) del livello core sono state calcolate tramite Teoria del Funzionale della Densità (DFT) con il software GPAW e il funzionale PBE.
   • Generazione Spettrale: 5.000 snapshot per composizione sono stati analizzati per generare spettri "stick" mediati per l'ensemble, che sono stati poi allargati tramite convoluzione Gaussiana (FWHM di 0,5 eV) per tenere conto degli effetti vibrazionali e di vita media.
   • Analisi del Descrittore: Per quantificare la relazione tra struttura e BE, è stato definito un descrittore strutturale locale $D(R)$ basato sul decadimento esponenziale dei contributi atomici rispetto al sito di ionizzazione.

Risultati Chiave

• Distribuzione Continua vs. Stati Discreti: Gli spettri sperimentali ricostruiti per strato e gli spettri mediati per l'ensemble simulati rivelano entrambi una distribuzione statistica continua di energie di legame piuttosto che picchi discreti. Le stechiometrie intermedie (ad esempio, SiO$_{1,0}$) mostrano un allargamento estremo, con lo spettro Si 2s per SiO$_{1,0}$ che copre oltre 5 eV.
• Origine Strutturale dell'Allargamento: L'allargamento è attribuito alla variazione continua degli ambienti atomici locali. Anche per una composizione macroscopica fissa $x$, la densità numerica locale dell'ossigeno e i motivi strutturali variano significativamente, portando a una vasta dispersione delle energie di legame del livello core.
• Raggio di Sensibilità Spettrale: Correlando le energie di legame con il descrittore strutturale $D(R)$, gli autori hanno determinato che l'ambiente spettroscopicamente rilevante si estende a un raggio di circa 5 Å (specificamente 5,2 Å per O 1s e 4,8 Å per Si 2s). Gli atomi oltre questa distanza contribuiscono in modo trascurabile (<0,1 eV) allo spostamento dell'energia di legame.
• Limiti di Rilevamento: Le simulazioni di un singolo atomo di ossigeno in una matrice di silicio hanno mostrato che anche l'ossigeno diluito può indurre un picco laterale distinto nello spettro Si 2s (spostato di ~1 eV), suggerendo che l'XPS del Si possa essere più sensibile all'ossigeno traccia rispetto all'XPS dell'O 1s in determinati contesti, a condizione che il rumore sperimentale sia basso.
• Validazione: Le simulazioni hanno riprodotto con successo i trend sperimentali, incluso lo spostamento dei picchi verso energie di legame inferiori con la diminuzione del contenuto di ossigeno e lo specifico allargamento osservato negli strati intermedi.

Significatività e Rivendicazioni
L'articolo sostiene che l'assegnazione convenzionale degli spettri XPS a stati chimici discreti (Si$^{1+}$, Si$^{2+}$, ecc.) è fondamentalmente errata per le interfacce silicio-ossido. Gli autori dimostrano che l'allargamento spettrale osservato deriva da una distribuzione statistica continua di energie di legame del livello core guidata dalla diversità strutturale locale, piuttosto che da una sovrapposizione di specie chimiche distinte.

Di conseguenza, gli autori affermano che:

1. Il Fitting Convenzionale è Qualitativo: Le procedure di fitting che tentano di identificare specifici motivi strutturali locali o stati di ossidazione discreti negli ossidi intermedi sono probabilmente inadeguate e possono portare a conclusioni errate sulle proprietà del film.
2. L'Ambiente Locale Domina: Lo spostamento del livello core è determinato dalla configurazione atomistica locale entro un raggio di ~5 Å, non solo dal numero di coordinazione del vicino più prossimo.
3. È Necessaria una Reinterpretazione: La "impronta digitale" di un particolare stato di ossidazione è sfumata dalle variazioni termiche e strutturali, rendendo necessaria una transizione da un'analisi basata su motivi discreti a una comprensione statistica ed basata sull'ensemble delle interfacce di ossido.

Questo lavoro mette in discussione l'interpretazione standard dei dati XPS degli ossidi nell'industria dei semiconduttori, suggerendo che i picchi larghi sono intrinseci alla natura statistica della struttura amorfa o disordinata dell'ossido piuttosto che indicatori di specifici strati intermedi.