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

Technische Samenvatting: De Rol van Lokale Structurele Variatie in het Röntgenfoto-elektronenspectrum van Siliciumoxide-interfaces

Probleemstelling
De interpretatie van röntgenfoto-elektronenspectroscopie (XPS)-gegevens voor siliciumoxide (SiO$_x$)-interfaces op siliciumsubstraten heeft traditioneel geleund op een "motief-gebaseerd" model. In dit conventionele beeld worden de brede Si 2p (en bij uitbreiding Si 2s) spectra waargenomen in substoichiometrische oxiden gedeconvolueerd in discrete componenten die overeenkomen met specifieke gehele oxidatietoestanden (Si$^{1+}$, Si$^{2+}$, Si$^{3+}$, Si$^{4+}$). Deze componenten worden vaak geïnterpreteerd als representaties van specifieke intermediaire oxidelagen of specifieke lokale bindingsmotieven (bijv. Si-atomen die gebonden zijn aan 1, 2 of 3 zuurstofatomen). Echter, de fysische eigenschappen van ultradunne SiO$_x$-films zijn minder gekarakteriseerd dan die van bulk SiO$_2$, en de betrouwbaarheid van het toewijzen van chemische toestanden op basis van enkelvoudige kernniveau-componenten wordt steeds vaker in twijfel getrokken, vooral naarmate oxidefilms dunner worden en substoichiometrische bijdragen dominant worden.

Methodologie
De auteurs hanteerden een tweeledige aanpak waarbij geavanceerde experimentele reconstructie werd gecombineerd met hoogwaardige statistische simulaties:

1. Experimentele Reconstructie:

   • Monster: Een commerciële Si(100)-wafer met een 40 nm dikke thermisch gegroeide SiO$_2$-laag.
   • Proces: XPS-gegevens werden verzameld bij de O 1s en Si 2s randen terwijl het oppervlak progressief werd gesputterd met een 300 eV Ar$^+$-bundel.
   • Reconstructie-algoritme: De ruwe spectra werden behandeld als diepte-geïntegreerde signalen die worden gedempt door inelastische elektronenverstrooiing. Gebruikmakend van een matrixinversiemethode ($S^{(obs)} = AS$) hebben de auteurs de laag-opgeloste spectra ($s_j$) voor individuele dieptes gereconstrueerd, waarbij problemen met spin-orbiet-splitsing werden vermeden door zich te concentreren op de O 1s en Si 2s randen.
   • Parameters: Een etsnelheid van 0,012 nm/s en een elektronische inelastische gemiddelde vrije weglengte ($\lambda_{el}$) van 3,5 nm werden gehanteerd.

2. Statistische Simulaties:

   • Modellering: Moleculaire Dynamica (MD) simulaties werden uitgevoerd op periodieke cellen (~15 Å) van SiO$_x$ met compositiesities $x = 0, 0,1, 0,5, 1,0, 1,5, 2,0$.
   • Potentiaal: Het MACE machine learning potentiaal werd gebruikt voor de MD-runs, die het smelten bij 3000 K en het gloeien (annealing) naar 300 K om diverse structurele ensembles te genereren.
   • Bindingenergieberekening: Kernniveau-bindingsenergieën (BE) werden berekend met behulp van Density Functional Theory (DFT) met de GPAW-software en de PBE-functionaal.
   • Spectrale Generatie: 5.000 snapshots per compositie werden geanalyseerd om ensemble-gemiddelde "stick spectra" te genereren, die vervolgens werden verbreed door Gaussische convolutie (0,5 eV FWHM) om rekening te houden met vibratie- en levensduureffecten.
   • Descriptor-analyse: Om de relatie tussen structuur en BE te kwantificeren, werd een lokale structurele descriptor $D(R)$ gedefinieerd op basis van de exponentiële afname van atomaire bijdragen relatief aan de ionisatieplaats.

Belangrijkste Resultaten

• Continue Verdeling versus Discrete Toestanden: De gereconstrueerde laag-opgeloste experimentele spectra en de gesimuleerde ensemble-gemiddelde spectra onthullen beide een continue statistische verdeling van bindingsenergieën in plaats van discrete pieken. De intermediaire stoichiometrieën (bijv. SiO$_{1,0}$) vertonen extreme verbreding, waarbij het Si 2s spectrum voor SiO$_{1,0}$ meer dan 5 eV beslaat.
• Structurele Oorsprong van de Verbreding: De verbreding wordt toegeschreven aan de continue variatie van lokale atomaire omgevingen. Zelfs voor een vaste macroscopische compositie $x$ variëren de lokale zuurstofgetal-dichtheid en structurele motieven aanzienlijk, wat leidt tot een brede spreiding van kernniveau-bindingsenergieën.
• Spectrale Sensitiviteitsstraal: Door de bindingsenergieën te correleren met de structurele descriptor $D(R)$, bepaalden de auteurs dat de spectroscopisch relevante omgeving zich uitstrekt tot een straal van ongeveer 5 Å (specifiek 5,2 Å voor O 1s en 4,8 Å voor Si 2s). Atomen buiten deze afstand dragen verwaarloosbaar (<0,1 eV) bij aan de verschuiving van de bindingsenergie.
• Detectielimieten: Simulaties van een enkel zuurstofatoom in een siliciummatrix toonden aan dat zelfs sporen van zuurstof een duidelijke zijpiek in het Si 2s spectrum kan induceren (verschoven met ~1 eV), wat suggereert dat Si XPS mogelijk gevoeliger is voor sporen van zuurstof dan O 1s XPS in bepaalde contexten, mits de experimentele ruis laag is.
• Validatie: De simulaties reproduceerden succesvol de experimentele trends, inclusclusief de verschuiving van pieken naar lagere bindingsenergieën bij afnemende zuurstofgehaltes en de specifieke verbreding waargenomen in intermediaire lagen.

Betekenis en Claims
Het artikel stelt dat de conventionele toewijzing van XPS-spectra aan discrete chemische toestanden (Si$^{1+}$, Si$^{2+}$, etc.) fundamenteel gebrekkig is voor siliciumoxide-interfaces. De auteurs demonstreren dat de waargenomen spectrale verbreding voortkomt uit een continue statistische verdeling van kernniveau-bindingsenergieën, gedreven door lokale structurele diversiteit, in plaats van een superpositie van afzonderlijke chemische soorten.

Concluderend beweren de auteurs dat:

1. Conventionele Fitting is Kwalitatief: Fittingprocedures die proberen specifieke lokale structurele motieven of discrete oxidatietoestanden te identificeren in intermediaire SiO$_x$, zijn waarschijnlijk ongeschikt en kunnen leiden tot onjuiste conclusies over de eigenschappen van de film.
2. Lokale Omgeving is Dominant: De kernniveau-verschuiving wordt bepaald door de lokale atomaire configuratie binnen een straal van ~5 Å, en niet alleen door het aantal directe buren (coördinatiegetal).
3. Herinterpretatie Vereist: De "vingerafdruk" van een specifieke oxidatietoestand wordt vertroebeld door thermische en structurele variaties, wat een verschuiving noodzakelijk maakt van discrete motief-gebaseerde analyse naar een statistisch, ensemble-gebaseerd begrip van oxide-interfaces.

Dit werk daagt de standaardinterpretatie van oxide XPS-data in de halfgeleiderindustrie uit, door te suggeren dat de brede pieken inherent zijn aan het statistische karakter van de amorfe of gedisorderde oxidestructuur, in plaats van indicatoren te zijn van specifieke intermediaire lagen.