🔬 materials science

Optimization-Based Thickness Estimation and Depth-Resolved FIB-Induced Damage Characterization via Multislice Electron Ptychography

Este artículo demuestra que la optimización bayesiana conjunta del espesor de la muestra y otros parámetros dentro de la ptychografía electrónica de multicapa permite una estimación precisa y autoconsistente del espesor de la lámina local y la caracterización de la resolución de profundidad del daño amorfo inducido por FIB a partir de un único conjunto de datos 4D-STEM.

Autores originales: Keun-Yeol Park, Chuqiao Shi, Murat Tuna Pamuk, Sooyoung Cheong, Seonyu Lee, Juhui Oh, Yu-Tsun Shao, Celesta S. Chang

Publicado 2026-10-08

📖 5 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Keun-Yeol Park, Chuqiao Shi, Murat Tuna Pamuk, Sooyoung Cheong, Seonyu Lee, Juhui Oh, Yu-Tsun Shao, Celesta S. Chang

Artículo original bajo licencia CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Esta es una explicación generada por IA del artículo a continuación. No ha sido escrita ni avalada por los autores. Para mayor precisión técnica, consulte el artículo original. Leer descargo de responsabilidad completo

Para ver los detalles más pequeños del mundo, los científicos suelen utilizar un haz de electrones en lugar de luz. Debido a que los electrones son tan pequeños, pueden revelar la disposición de los átomos individuales dentro de un material. Sin embargo, cuando estos electrones pasan a través de un objeto sólido, no simplemente viajan en línea recta; rebotan en los átomos e interfieren entre sí de formas complejas. Este fenómeno, conocido como dispersión dinámica, significa que la imagen final depende en gran medida de qué tan grueso sea el objeto. Si un científico no conoce el espesor exacto de la muestra, la imagen resultante puede ser engañosa, dificultando la distincción entre una característica estructural real y un artefacto causado por la profundidad del material. Esta incertidumbre se convierte en un obstáculo importante al estudiar materiales preparados mediante un haz de iones focalizado, una técnica que utiliza un flujo de iones para tallar diminutas y finas rebanadas de materia para su examen. Aunque este método es esencial para crear muestras, el proceso de corte en sí mismo daña la superficie, creando una capa de material desordenado que se asienta sobre el cristal prístino del interior. Hasta ahora, medir el espesor de estas rebanadas y la profundidad del daño requería la preparación de una muestra especial y separada solo para cortarla y observar el lateral, un proceso que consumía mucho tiempo y que a menudo era imposible de realizar en la pieza exacta de material que se estaba estudiando.

Un equipo de investigadores ha desarrollado una forma de resolver este problema utilizando un único conjunto de datos, lo que permite determinar tanto el espesor local de una muestra como la profundidad del daño sin necesidad de un espécimen separado. Trabajaron con una fina rebanada de silicio que tenía forma de cuña, lo que significa que era más gruesa en un extremo y más delgada en el otro. Esta rebanada había sido preparada utilizando cuatro voltajes de corte finales diferentes, creando cuatro regiones distintas donde se esperaba que la superficie dañada variara. Los investigadores utilizaron una poderosa técnica de imagen llamada pticografía de multisuperficies (multislice electron ptychography), que trata el haz de electrones como una onda que pasa a través del material en capas. Para que esta técnica funcione, el modelo computacional necesita conocer el espesor total de la muestra para empezar, pero ese es precisamente el dato que los científicos intentaban encontrar. Para romper este problema circular, emplearon una estrategia matemática llamada optimización bayesiana. En lugar de adivinar el espesor, la computadora probó miles de combinaciones diferentes de espesor, enfoque y número de capas en su modelo. Buscó la combinación específica que hiciera que los patrones de electrones simulados coincidieran más estrechamente con los patrones reales capturados por el microscopio.

Los resultados mostraron que este método podía estimar de manera fiable el espesor físico de la rebanada de silicio. Cuando los investigadores compararon las estimaciones de la computadora con mediciones independientes tomadas con una técnica diferente, encontraron que los valores eran muy cercanos, con una diferencia de solo un pequeño margen constante. Más importante aún, el método rastreó con éxito cómo cambiaba el espesor a través de la muestra en forma de cuña, confirmando que la computadora estaba midiendo una propiedad física real en lugar de simplemente encontrar un truco matemático. Una vez establecido el espesor correcto, los investigadores utilizaron los mismos datos para mirar en el interior de la rebanada. Debido a que la técnica reconstruye la muestra capa por capa, pudieron ver la transición desde la superficie desordenada y dañada hasta el cristal perfecto y ordenado en el medio. Midieron el espesor de esta capa dañada para cada una de las cuatro regiones. Los datos revelaron que el daño era más delgado donde el voltaje de corte era más bajo, midiendo 3.2 nanómetros a 2 kilovoltios, y más grueso donde el voltaje era más alto, midiendo 26.4 nanómetros a 30 kilovoltios. Estas mediciones coincidieron con hallazgos previos de estudios transversales separados, pero con la ventaja de que fueron derivadas directamente de la muestra destinada al análisis atómico.

Este enfoque ofrece un avance significativo para la microscopía cuantitativa, donde las mediciones precisas son críticas. Al permitir que los científicos determinen el espesor local y la extensión del daño superficial a partir del mismo conjunto de datos utilizado para obtener imágenes de los átomos, el método elimina la necesidad de pasos de preparación adicionales que podrían alterar la muestra. Proporciona una forma autoconsistente de comprender cómo está estructurado el material y cómo fue afectado por el proceso de preparación. Los investigadores encontraron que la capa de daño crecía de manera constante a medida que aumentaba el voltaje de corte, confirmando que los voltajes más bajos son, de hecho, mejores para preservar la estructura cristalina cerca de la superficie. Si bien persiste una pequeña diferencia sistemática entre sus estimaciones y otros métodos de medición, probablemente debido a la compleja física de cómo se dispersan los electrones, la capacidad de mapear estas variaciones dentro de una misma muestra abre la puerta a interpretaciones más precisas de las imágenes a escala atómica. Este trabajo demuestra que, con las herramientas computacionales adecuadas, el microscopio puede decirnos no solo cómo se ven los átomos, sino también qué tan grueso es el material y qué tan profundo llega el daño, todo al mismo tiempo.

Resumen Técnico: Estimación de Espesor Basada en la Optimización y Caracterización del Daño Inducido por FIB con Resolución de Profundidad mediante la Ptygrafía Electrónica Multislice

Planteamiento del Problema
La microscopía electrónica cuantitativa, particularmente en la microscopía electrónica de transmisión (TEM) y la microscopía electrónica de barrido de transmisión (STEM), depende en gran medida de un conocimiento preciso del espesor de la muestra. La dispersión dinámica provoca que el contraste de la imagen y las intensidades de difracción dependan de forma no lineal de la distancia de propagación, lo que convierte al espesor en una variable experimental determinante. En la ptygrafía electrónica multislice (MEP), el espesor de la muestra es un parámetro de entrada requerido para el modelo de propagación hacia adelante (forward model). Sin embargo, para muestras preparadas mediante fresado con haz de iones focalizados (FIB), el espesor local suele ser desconocido y varía espacialmente (por ejemplo, en laminillas en forma de cuña). Además, el espesor está acoplado con otros parámetros de reconstrucción, específicamente el desenfoque de la sonda y el número de capas (slices), lo que hace que su determinación independiente a partir de la propia reconstrucción sea no trivial. Adicionalmente, las laminillas preparadas por FIB suelen consistir en un núcleo cristalino limitado por capas superficiales amorfas causadas por la implantación de iones y el daño, las cuales, en las imágenes de proyección convencionales, se superponen con el núcleo cristalino a lo largo de la dirección del haz, dificultando la distinción de su contribución a la dispersión o la cuantificación de la profundidad del daño sin la preparación de muestras transversales separadas.

Metodología
Los autores proponen un marco que acopla la ptygrafía electrónica multislice con la optimización bayesiana (BO) para estimar conjuntamente el espesor de la muestra, el desenfoque de la sonda y el número de capas a partir de un único conjunto de datos 4D-STEM.

Preparación de la Muestra: Se preparó una laminilla de silicio (Si) en forma de cuña a partir de una oblea de Si (001) utilizando un FIB-SEM de doble haz. La laminilla se dividió en cuatro regiones, cada una sometida a un voltaje de fresado final de 2, 5, 8 o 30 kV para variar sistemáticamente la extensión de la amorfización inducida por el FIB.
Adquisición de Datos: Se adquirieron conjuntos de datos 4D-STEM a 200 kV utilizando un detector de electrones directo. Los conjuntos de datos cubrieron un campo de visión de 9.2 × 9.2 nm² con un semianulo de convergencia de la sonda de 30 mrad.
Optimización Bayesiana (BO): Se empleó un motor de BO para optimizar conjuntamente tres parámetros: el espesor total de la muestra ( $t$ ), el desenfoque de la sonda ( $\Delta f$ ) y el número de capas ( $N$ ). La función objetivo se definió como el logaritmo negativo de la diferencia de la raíz cuadrática media entre las amplitudes de difracción medidas y simuladas. Para asegurar que la estimación del espesor fuera físicamente significativa y no un artefacto numérico, se desactivaron múltiples estados de la sonda y la relajación ortogonal de la sonda durante la optimización para evitar que otros grados de libertad compensaran un espesor incorrecto.
Caracterización con Resolución de Profundidad: Utilizando los parámetros optimizados, se realizaron reconstrucciones MEP completas. Los volúmenes 3D resultantes se analizaron para separar el interior cristalino de las capas superficiales amorfas inducidas por el FIB. El espesor de las capas amorfas se cuantificó ajustando los perfiles de fase versus profundidad de las columnas atómicas con una función de error doble simétrica para identificar las interfaces amorfo-cristalinas.
Validación: Las estimaciones de espesor derivadas del método BO-MEP se compararon con mediciones independientes obtenidas mediante espectroscopía de pérdida de energía de electrones (EELS) utilizando el método de la relación de logaritmos (log-ratio method).

Resultados Clave

Estimación de Espesor: La optimización conjunta produjo un mínimo bien definido y único para el espesor, el desenfoque y el número de capas. Al aplicarse a diferentes subregiones de la laminilla en forma de cuña, los valores de espesor optimizados reprodujeron la tendencia espacial esperada de aumento de espesor de arriba hacia abajo.
Comparación con EELS: Las estimaciones de espesor por BO-MEP mostraron un desfase sistemático del 10–17% respecto a las mediciones de EELS (subestimando el espesor). Sin embargo, los autores señalan que el espesor absoluto de EELS depende del camino libre medio inelástico asumido, el cual tiene una precisión de aproximadamente ±20%. A pesar del desfase, la capacidad del método para reproducir las variaciones locales de espesor independientemente de la configuración de la reconstrucción respalda su interpretación como una estimación física.
Robustez: Se encontró que el espesor optimizado es robusto frente a variaciones en la fuerza de regularización de la capa, mientras que el número óptimo de capas se desplazó para mantener un muestreo axial consistente.
Caracterización del Daño: Las reconstrucciones con resolución de profundidad lograron separar las capas superficiales amorfas del núcleo cristalino. Los espesores de las capas amorfas cuantificados aumentaron monótonamente con el voltaje de fresado: 3.2 nm (2 kV), 5.0 nm (5 kV), 7.0 nm (8 kV) y 26.4 nm (30 kV). Estos valores concuerdan razonablemente con mediciones transversales reportadas previamente, aunque la región de 30 kV mostró mayor variabilidad debido a la complejidad de la muestra más gruesa y la región distorsionada extendida.

Significado y Reivindicaciones
El artículo sostiene que la aplicación de MEP acoplada con optimización bayesiana sobre un único conjunto de datos 4D-STEM proporciona una vía hacia un análisis cuantitativo auticonsistente y consciente del espesor. La importancia principal reside en dos áreas:

Interpretabilidad Física: El estudio demuestra que el espesor seleccionado por la optimización no es meramente un parámetro numérico que minimiza el error de reconstrucción, sino que puede interpretarse como una estimación significativa del espesor físico local de la muestra, siempre que se utilice la optimización conjunta para resolver el acoplamiento de parámetros.
Evaluación Integrada del Daño: El método permite la determinación simultánea del espesor local de la muestra y la caracterización con resolución de profundidad del daño inducido por FIB (espesor de la capa amorfa) directamente desde la misma laminilla utilizada para el análisis a resolución atómica. Esto elimina la necesidad de preparar muestras transversales separadas para evaluar la profundidad del daño o el espesor total.

Los autores posicionan este enfoque como una herramienta de refinamiento que requiere solo una estimación previa aproximada del espesor para definir los límites de búsqueda, lo que lo hace adecuado para su integración en flujos de trabajo rutinarios de MEP para restringir simulaciones y mejorar la interpretación de mediciones fuertemente afectadas por la dispersión dinámica.

Más como este