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

Resumen Técnico: El papel de la variación estructural local en el espectro de fotoelectrones de rayos X de las interfaces de óxido de silicio

Planteamiento del problema
La interpretación de los datos de espectroscopía de fotoelectrones de rayos X (XPS) para las interfaces de óxido de silicio (SiO$_x$) sobre sustratos de silicio se ha basado tradicionalmente en un modelo "basado en motivos". En esta visión convencional, los espectros anchos de Si 2p (y, por extensión, de Si 2s) observados en óxidos subestequiométricos se deconvolucionan en componentes discretos correspondientes a estados de oxidación enteros específicos (Si$^{1+}$, Si$^{2+}$, Si$^{3+}$, Si$^{4+}$). Estos componentes se interpretan a menudo como representaciones de capas de óxido intermedias distintas o motivos de enlace local específicos (por ejemplo, átomos de Si unidos a 1, 2 o 3 átomos de oxígeno). Sin embargo, las propiedades físicas de las películas de SiO$_x$ ultradelgadas están menos caracterizadas que las del SiO$_2$ voluminoso, y la fiabilidad de la asignación de estados químicos basados en componentes de nivel de núcleo únicos es cada vez más cuestionada, particularmente a medida que las películas de óxido se adelgazan y las contribuciones subestequiométricas se vuelven dominantes.

Metodología
Los autores emplearon un enfoque dual que combina la reconstrucción experimental avanzada con simulaciones estadísticas de alto rendimiento:

1. Reconstrucción Experimental:

   • Muestra: Una oblea de Si(100) comercial con una capa de SiO$_2$ crecida térmicamente de 40 nm.
   • Proceso: Se recolectaron datos de XPS en los bordes O 1s y Si 2s mientras se erosionaba progresivamente la superficie con un haz de Ar$^+$ de 300 eV.
   • Algoritmo de Reconstrucción: Los espectros brutos fueron tratados como señales integradas en la profundidad y atenuadas por la dispersión inelástica de electrones. Utilizando un método de inversión de matriz ($S^{(obs)} = AS$), los autores reconstruyeron los espectros resueltos por capa ($s_j$) para profundidades individuales, evitando problemas de desdoblamiento de espín-órbita al centrarse en los bordes O 1s y Si 2s.
   • Parámetros: Se utilizó una tasa de erosión de 0.012 nm/s y un camino libre medio inelástico de electrones ($\lambda_{el}$) de 3.5 nm.

2. Simulaciones Estadísticas:

   • Modelado: Se realizaron simulaciones de Dinámica Molecular (MD) en celdas periódicas (~15 Å) de SiO$_x$ con composiciones $x = 0, 0.1, 0.5, 1.0, 1.5, 2.0$.
   • Potencial: Se utilizó el potencial de aprendizaje automático MACE para las corridas de MD, que incluyeron fusión a 3000 K y recocido a 300 K para generar diversos conjuntos estructurales.
   • Cálculo de la Energía de Ligadura: Las energías de ligadura (BE) de los niveles de núcleo se calcularon mediante la Teoría del Funcional de la Densidad (DFT) utilizando el software GPAW y el funcional PBE.
   • Generación de Espectros: Se analizaron 5,000 instantáneas por composición para generar espectros de varillas (stick spectra) promediados por conjunto, los cuales fueron ensanchados mediante convolución Gaussiana (FWHM de 0.5 eV) para dar cuenta de los efectos vibracionales y de tiempo de vida.
   • Análisis de Descriptores: Para cuantificar la relación entre estructura y BE, se definió un descriptor estructural local $D(R)$ basado en el decaimiento exponencial de las contribuciones atómicas respecto al sitio de ionización.

Resultados Clave

• Distribución Continua vs. Estados Discretos: Los espectros experimentales reconstruidos por capa y los espectros simulados promediados por conjunto revelan ambos una distribución estadística continua de energías de ligadura en lugar de picos discretos. Las estequiometrías intermedias (por ejemplo, SiO$_{1.0}$) exhiben un ensanchamiento extremo, donde el espectro Si 2s para SiO$_{1.0}$ abarca más de 5 eV.
• Origen Estructural del Ensanchamiento: El ensanchamiento se atribuye a la variación continua de los entornos atómicos locales. Incluso para una composición macroscópica fija $x$, la densidad de número de átomos de oxígeno local y los motivos estructurales varían significativamente, lo que conduce a una amplia dispersión de las energías de ligadura de los niveles de núcleo.
• Radio de Sensibilidad Espectroscópica: Al correlacionar las energías de ligadura con el descriptor estructural $D(R)$, los autores determinaron que el entorno espectroscópicamente relevante se extiende a un radio de aproximadamente 5 Å (específicamente 5.2 Å para O 1s y 4.8 Å para Si 2s). Los átomos más allá de esta distancia contribuyen de manera insignificante (<0.1 eV) al desplazamiento de la energía de ligadura.
• Límites de Detección: Las simulaciones de un solo átomo de oxígeno en una matriz de silicio mostraron que incluso el oxígeno diluido puede inducir un pico lateral distinto en el espectro Si 2s (desplazado por ~1 eV), lo que sugiere que el XPS de Si puede ser más sensible al oxígeno traza que el XPS de O 1s en ciertos contextos, siempre que el ruido experimental sea bajo.
• Validación: Las simulaciones reprodujeron con éxito las tendencias experimentales, incluyendo el desplazamiento de los picos hacia energías de ligadura más bajas al disminuir el contenido de oxígeno, así como el ensanchamiento específico observado en las capas intermedias.

Significancia y Reivindicaciones
El artículo argumenta que la asignación convencional de los espectros XPS a estados químicos discretos (Si$^{1+}$, Si$^{2+}$, etc.) es fundamentalmente errónea para las interfaces de óxido de silicio. Los autores demuestran que el ensanchamiento espectral observado surge de una distribución estadística continua de energías de ligadura de los niveles de núcleo impulsada por la diversidad estructural local, en lugar de una superposición de especies químicas distintas.

Consecuentemente, los autores afirman que:

1. El Ajuste Convencional es Cualitativo: Los procedimientos de ajuste que intentan identificar motivos estructurales locales específicos o estados de oxidación discretos en SiO$_x$ intermedios son probablemente inadecuados y pueden conducir a conclusiones incorrectas sobre las propiedades de la película.
2. El Entorno Local Domina: El desplazamiento del nivel de núcleo está determinado por la configuración atómica local dentro de un radio de ~5 Å, no solo por el número de coordinación del vecino más cercano.
3. Se Requiere una Reinterpretación: La "huella dactilar" de un estado de oxidación específico se desdibuja por las variaciones térmicas y estructurales, lo que requiere un cambio de un análisis basado en motivos discretos a un entendimiento estadístico y basado en conjuntos de las interfaces de óxido.

Este trabajo desafía la interpretación estándar de los datos de XPS de óxidos en la industria de los semiconductores, sugiriendo que los picos anchos son intrínsecos a la naturaleza estadística de la estructura amorfa o desordenada del óxido, en lugar de ser indicadores de capas intermedias específicas.