Cálculo de EBSD¶
EBSD (difracción de electrones retrodispersados) utiliza el mismo núcleo de Bethe/ondas de Bloch que CBED y STEM, pero el problema se plantea de otra manera. CBED y STEM son problemas de haz incidente: una onda electrónica entra desde el exterior de la muestra y se calcula la onda de salida. EBSD es un problema de dirección de salida: los electrones que han sufrido dispersión inelástica en el interior de la muestra emergen como electrones retrodispersados, y el cálculo pregunta cuánta intensidad sale en cada dirección externa.
ReciPro convierte este problema de dirección de salida en un problema ordinario de haz incidente mediante el teorema de reciprocidad. Primero calcula un master pattern en el espacio de direcciones y luego combina ese master pattern con los pesos de Monte Carlo para profundidad / energía / dirección y la geometría del detector para formar el patrón del detector.
Reformulación con el teorema de reciprocidad¶
Si la amplitud desde un punto fuente interno \(\mathbf r_n\) hacia una dirección externa \(\widehat{\mathbf s}\) se calculara directamente, sería necesario un problema de dispersión separado para cada punto fuente. Eso no es práctico.
El teorema de reciprocidad reescribe el problema de la siguiente manera: la amplitud de que un electrón que parte de \(\mathbf r_n\) aparezca en la dirección de campo lejano \(\widehat{\mathbf s}\) es igual a la amplitud, en \(\mathbf r_n\), de una onda recíproca que incide desde la dirección externa \(-\widehat{\mathbf s}\). Esta onda recíproca es una solución ordinaria de Bethe/ondas de Bloch. Escribiéndola como \(\psi_{\widehat{\mathbf s}}^{\mathrm{rec}}(\mathbf r)\), la intensidad EBSD en la dirección \(\widehat{\mathbf s}\) puede escribirse como
donde \(\sigma_n(E,z)\) es el peso para la dispersión inelástica cerca de la posición atómica \(\mathbf r_n\) en el canal de retrodispersión a energía \(E\) y profundidad \(z\). Los términos fuente se suman como intensidades, no como una suma coherente de amplitudes, porque se supone que la dispersión inelástica destruye la relación de fase entre las distintas posiciones fuente.
Master pattern¶
El master pattern de EBSD almacena la parte de difracción dinámica específica del cristal de la expresión anterior sobre una malla de direcciones. Conceptualmente,
donde \(w_n\) es el peso de fuente inelástica del lado del cristal en la posición atómica \(\mathbf r_n\). ReciPro utiliza el peso empírico
con el número atómico \(Z_n\) y la ocupación \(\mathrm{occ}_n\). Esto es independiente de la distribución de profundidad de transporte / energía producida por Monte Carlo.
En la implementación, la onda de Bloch recíproca se evalúa en cada posición atómica:
El código forma entonces la matriz de pares de ondas de Bloch
y la integral analítica sobre el espesor
de modo que el master pattern se evalúa como
En el límite degenerado en el que el denominador se acerca a cero, \(\mathcal F_{jj'}(t)\to t\).
Muestreo del espacio de direcciones¶
El master pattern no es la imagen del detector en sí; es una distribución de intensidad en el espacio de direcciones fijo al cristal. ReciPro muestrea ese espacio de direcciones con una proyección equiárea de Roşca–Lambert y almacena los hemisferios \(+Z\) y \(-Z\) como matrices planas separadas. El muestreo equiárea reduce el sesgo de densidad entre los polos y el ecuador.
En esta etapa, el master pattern depende de la estructura cristalina, la tensión de aceleración, la profundidad, la energía y el modelo de absorción. La geometría del detector, como el centro del patrón y la posición de la pantalla, aún no se ha aplicado.
Pesos de Monte Carlo y patrón del detector¶
Para obtener un patrón de detector EBSD cercano al observable experimental, el master pattern debe ponderarse según cuántos electrones retrodispersados emergen de cada profundidad, energía y dirección. Escribiendo este peso de transporte como
y usando \(\widehat{\mathbf s}(\mathbf p)\) para la dirección fija al cristal que corresponde al píxel del detector \(\mathbf p\), el patrón del detector final es
como una suma discreta sobre energía y profundidad.
La parte de Monte Carlo sigue la dispersión elástica, la dispersión inelástica, la pérdida de energía y el escape a través de la superficie de la muestra. Para los electrones retrodispersados construye distribuciones de profundidad, energía y dirección de salida. ReciPro distingue entre modelos que utilizan la última posición de dispersión inelástica y la energía inmediatamente posterior como fuente efectiva, y modelos que utilizan la profundidad de escape y la energía de escape.
Profundidad de la fuente y longitud de trayectoria¶
Por reciprocidad, el master pattern para una dirección \(\widehat{\mathbf s}\) se resuelve como una onda plana que entra a lo largo de \(-\widehat{\mathbf s}\), tomando la superficie del cristal perpendicular a \(\widehat{\mathbf s}\) para cada dirección, de modo que su coordenada de espesor \(t\) se mide a lo largo de la dirección de salida. En cambio, la profundidad de la fuente \(z\) de Monte Carlo se mide a lo largo de la normal de la superficie inclinada de la muestra. Con \(\mu=\cos\chi\), donde \(\chi\) es el ángulo entre la dirección de salida y la normal a la superficie, una fuente a profundidad \(z\) se encuentra, a lo largo de la dirección de salida, a la longitud de trayectoria
de la superficie. ReciPro ajusta la distribución de profundidad de cada bin de dirección de salida y de cada energía con una exponencial \(\propto e^{-z/\lambda_z(E)}\), la convierte en cada píxel del detector con el \(\mu\) de ese píxel a \(\lambda_t=\lambda_z/\mu\) e integra la distribución resultante, normalizada dentro de la matriz de profundidades, sobre cada intervalo de espesor del master pattern. Por tanto, en la suma anterior el índice de profundidad recorre el espesor \(t_j\), y el peso de un píxel depende de su \(\mu\) además de \(\widehat{\mathbf s}\). En la salida rasante (\(\mu\) pequeño), una fuente a la misma profundidad recorre un camino más largo, por lo que el peso se desplaza hacia láminas más gruesas.
La matriz de profundidades automática termina en el percentil 99,9 de la longitud de trayectoria \(z/\mu\) de los electrones que alcanzan el detector, redondeado hacia arriba a dos cifras significativas, y tiene 40 puntos más próximos entre sí cerca de la superficie.
Fondo de TDS y modelo de absorción¶
Los patrones EBSD contienen no solo la estructura geométrica de bandas de Kikuchi, sino también un fondo suave procedente de la dispersión térmica difusa (TDS). Cuando IncludeTDSBackground está activado, ReciPro evalúa la componente de TDS dispersada hacia el hemisferio posterior,
como una matriz de absorción \(\mu_{\mathrm{back}}\) y añade la intensidad de fondo usando la misma sumación de pares de ondas de Bloch que el master pattern. Como se reutiliza la misma solución de autovalores, el fondo de TDS añade relativamente poco coste adicional.
Cuando UseNonLocalAbsorption está activado, el potencial absorptivo no se trata solo como \(U'_{\mathbf g-\mathbf h}\), sino como una forma no local que depende de la dirección y de los pares de haces. Esto puede mejorar la precisión, pero también requiere reconstruir la matriz de absorción para las direcciones de la malla del master pattern, por lo que puede aumentar considerablemente el tiempo de cálculo.
Parámetros prácticos¶
- Número de haces: Un número demasiado bajo de haces pierde el detalle de las bandas de Kikuchi y la estructura de bandas HOLZ. Los ejes de zona de índices bajos pueden requerir varios cientos de haces.
- Matrices de profundidad y energía: Si son más gruesas que la escala de variación del peso de Monte Carlo \(W(E,z;\widehat{\mathbf s})\), los efectos de ancho de banda dependiente de la energía y de profundidad de canalización se promedian y se pierden.
- Geometría del detector: El centro del patrón, la distancia a la pantalla y la inclinación de la muestra determinan la aplicación \(\widehat{\mathbf s}(\mathbf p)\), de modo que el patrón del detector puede cambiar incluso cuando el master pattern permanece inalterado.
- Interpretación de la reciprocidad: El master pattern no es la imagen del detector. Solo se convierte en un patrón de detector tras la ponderación de Monte Carlo y la proyección al detector.
- Fondo de TDS: Actívelo para comparaciones cuantitativas del contraste de bandas. Desactívelo cuando la estructura geométrica de Kikuchi sea más fácil de inspeccionar sin el fondo suave.