Theoretical ab initio Evolution of Satellite Intensity near Threshold for Cu K-shell transitions

Resumen Técnico: Evolución Teórica Ab Initio de la Intensidad de Satélite cerca del Umbral para Transiciones de la Cáscara K de Cu

Planteamiento del Problema
Los cálculos de estructura atómica teórica para sistemas complejos suelen basarse en la aproximación súbita (SA, por sus siglas en inglés), la cual asume que la ionización o la desintegración nuclear ocurren en una escala de tiempo mucho más corta que la relajación electrónica. Si bien la SA proporciona una excelente descripción de los procesos de alta energía, su validez es cuestionable cerca de los umbrales de ionización. En este régimen de baja energía, el electrón saliente posee una energía cinética baja y permanece en las proximidades del sistema atómico durante un tiempo comparativamente largo, lo que viola los supuestos fundamentales de la SA. En consecuencia, el sistema se aproxima a un límite adiabático donde la estructura electrónica evoluciona continuamente. Las observaciones experimentales previas de este efecto, particularmente en la espectroscopia de absorción de rayos X (XAS), se han visto limitadas por las restricciones de resolución. Además, mientras que existen modelos teóricos como los de Thomas y Roy para predecir la evolución de la intensidad de sacudida (líneas satélite) desde el régimen adiabático al súbito, existe la necesidad de validar estos modelos frente a datos de emisión de alta resolución y determinar si los métodos ab initio estándar pueden simular con precisión esta evolución sin parámetros de ajuste empíricos para la intensidad.

Metodología
Los autores emplearon un enfoque integral ab initio utilizando métodos de Multiconfiguración Dirac-Fock (MCDF) de vanguardia, específicamente utilizando el código MCDFGME. La metodología consistió en varias etapas clave:

1. Cálculos de Estructura Atómica: Se calcularon todas las transiciones de la cáscara K de Cu, incluyendo las tasas radiativas y no radiativas. Las funciones de onda atómicas se optimizaron mediante el método de nivel optimizado con relajación completa para las tasas radiativas, mientras que las tasas no radiativas utilizaron funciones de onda de electrones en el continuo sin relajación.
2. Determinación de la Probabilidad de Sacudida (Shake): Se calculó teóricamente la probabilidad de los procesos de sacudida de excitación (shake-up) y de ionización (shake-off). La probabilidad total de sacudida se derivó del producto interno de las funciones de onda radiales de los estados inicial y final. Los autores distinguieron entre las probabilidades de sacudida de un solo electrón (utilizando un enfoque de distribución binomial) y las probabilidades totales de sacudida, contabilizando las excitaciones monopolares.
3. Generación de Espectros Sintéticos: Se generó un espectro sintético para las transiciones Cu K$\alpha_{1,2}$ ponderando todas las rutas de decaimiento posibles mediante sus tasas y pesos estadísticos calculados.
4. Modelado de la Evolución de la Energía: Para simular la evolución de la intensidad de la línea cerca del umbral, los autores desarrollaron una función de modulación. Esto implicó la convolución del ancho Lorentziano natural de los niveles ionizados con un perfil Gaussiano que representa la distribución de energía del haz de excitación. Se utilizó un "Perfil Lorentziano Efectivo" para dar cuenta de la naturaleza asimétrica de la sección eficaz de ionización cerca del umbral.
5. Incorporación de Excitaciones Resonantes: Reconociendo que los datos de ionización por sí solos no podían explicar los resultados experimentales en la región cercana al umbral, los autores incorporaron excitaciones resonantes de las fases de óxido de Cu(I) y Cu(II) (específicamente transiciones $1s \to 3d$ y $1s \to 4p$) identificadas mediante datos de Estructura Fina de Absorción de Rayos X (XANES).
6. Ajuste y Análisis: Los resultados teóricos se compararon con los datos experimentales de Galambosi et al. [8]. Se empleó un análisis bayesiano mediante muestreo anidado (ajuste anidado) para determinar los desplazamientos de energía y los multiplicadores de ancho más probables para diferentes procesos de decaimiento (diagrama, shake-off, shake-up y excitación) para alinear la simulación con las observaciones experimentales.

Contribuciones Clave

• Validación de Métodos Ab Initio: El estudio demuestra que los métodos ab initio estándar de vanguardia, sin parámetros de ajuste de intensidad para la ionización, pueden simular con éxito la evolución de la intensidad de las líneas satélite cerca de los umbrales de ionización.
• Cuantificación de las Probabilidades de Sacudida: Los autores calcularon una probabilidad de sacudida total para la ionización de la cáscara K de Cu de aproximadamente 25.97%, con un desglose de $\approx$13.08% para shake-off y $\approx$12.89% para shake-up.
• Identificación de los Orígenes del Umbral: El trabajo identifica que la intensidad de la satélite "por debajo del umbral" observada en experimentos no se origina únicamente en los procesos de sacudida, sino que está significativamente influenciada por excitaciones resonantes en las fases de óxido de Cu(I) ($1s \to 4p$) y Cu(II) ($1s \to 3d$).
• Refinamiento de Modelos Teóricos: El artículo proporciona una prueba rigurosa de los modelos de Thomas y Roy. Confirma la forma funcional del modelo de Thomas para la transición adiabática-a-súbita y destaca las discrepancias en los parámetros del modelo de Roy cuando se ajustan a los datos experimentales, sugiriendo que los datos experimentales contienen procesos físicos (como efectos de estado sólido) no capturados totalmente por los modelos atómicos simples.

Resultados

• Acuerdo con el Experimento: Las relaciones de intensidad simuladas de las líneas satélite respecto a la intensidad total de K$\alpha_{1,2}$ mostraron un buen acuerdo con los datos experimentales para energías de excitación por encima de 9050 eV, donde todos los canales de decaimiento están abiertos. La relación simulada a 10 keV resultó ser $0.262 \pm 0.00004$, consistente con la probabilidad de sacudida total calculada en la aproximación súbita.
• Comportamiento del Umbral: Por debajo de 9050 eV, la simulación de solo ionización falló en reproducir el aumento experimental en la intensidad de la satélite. Sin embargo, al incluir las excitaciones resonantes de los óxidos de Cu, la simulación logró reproducir con éxito los datos experimentales en el rango de 8975–9060 eV.
• Ajuste de Modelos: El ajuste del modelo de Thomas a la simulación arrojó parámetros ($\Delta E$, $r$, $\mu_\infty$) que estuvieron en buen acuerdo tanto con el ajuste experimental como con los cálculos teóricos. El ajuste del modelo de Roy mostró una proximidad razonable a los valores teóricos, pero desviaciones significativas en los parámetros de energía de enlace, lo que indica las limitaciones del modelo para capturar la complejidad total del espectro experimental.
• Características Espectrales: La inclusión de transiciones de sacudida de orden superior (por ejemplo, shake-off de $3s$, $3p$, $2p$ y shake-up de $3d$) alteró significativamente la forma de la línea de la satélite teórica en comparación con cálculos previos de Dirac-Fock, resaltando la necesidad de cálculos detallados de sacudida para la metrología de alta resolución.

Significancia
El artículo afirma que su principal significancia radica en demostrar que los cálculos teóricos estándar pueden predecir la evolución media de la intensidad espectral cerca de los umbrales de ionización sin depender de ajustes de intensidad empíricos. Al separar con éxito la intensidad de la satélite (originada en procesos de sacudida) de la intensidad directa de diagrama e identificar la contribución de las excitaciones resonantes en las fases de óxido, los autores proporcionan una imagen física más completa de la región cercana al umbral. El trabajo subraya la importancia de los cálculos detallados de shake-up y shake-off para la espectroscopia de rayos X de vanguardia y sugiere que las discrepancias entre la teoría y el experimento en este régimen a menudo pueden atribuirse a efectos de estado sólido (como la presencia de óxidos) y a las formas de línea específicas utilizadas en el análisis de datos, más que a un fallo de los modelos de física atómica subyacentes.