Computer simulations of the Stark effect in the helium-beta complex of krypton in ICF conditions

Resumen Técnico: Simulaciones Computacionales del Efecto Stark en el Complejo de Helio-$\beta$ de Kriptón en Condiciones de ICF

Planteamiento del Problema
La espectroscopia de experimentos de fusión por confinamiento inercial (ICF) depende de dopantes, como el kriptón, para diagnosticar la temperatura y la densidad de los plasmas en implosión. Específicamente, las transiciones de la capa K de especies tipo He (He-like), particularmente la línea He-$\beta$ (transición entre $n=3$ y $n=1$), son críticas para las mediciones de densidad. Sin embargo, la línea He-$\beta$ está acompañada de transiciones satélite tipo Li (Li-like) ($1sn'\ell'3\ell \to 1s^2n'\ell'$) que se solapan parcialmente con la línea principal. Modelar el complejo He-$\beta$ completo bajo condiciones extremas de ICF (densidades electrónicas $n_e = 10^{24}$ a $10^{25}$ cm$^{-3}$ y temperaturas $T_e = 3$ keV) es computacionalmente prohibitivo debido al gran número de niveles atómicos y de los elementos de matriz de dipolo no nulos involucrados. Si bien los modelos de simulación computacional (CSM) son las herramientas más avanzadas para calcular perfiles de Stark, su aplicación a tales formas de línea complejas ha sido un desafío. Además, estudios previos utilizando aproximaciones de "teoría estándar" (que asumen iones estáticos y electrones de impacto) sugirieron que los términos de interferencia en el operador de impacto electrónico eran despreciables para estas transiciones, un hallazgo que requería verificación bajo condiciones de simulación completa.

Metodología
Los autores presentan un estudio comparativo utilizando múltiples códigos de simulación semiclásicos para modelar el ensanchamiento Stark de la línea He-$\beta$ del Kr y sus satélites tipo Li de $n=2$ y $n=3$. Los códigos empleados incluyen:

• SIMULA y SIMULAm: Modelos donde las partículas siguen trayectorias rectas (rectilíneas), tratando a los perturbadores como independientes. SIMULAm es un código híbrido que incorpora el operador de ensanchamiento electrónico de la teoría estándar en cada paso de tiempo para reducir el costo computacional.
• SimU y SimUSP: Modelos que resuelven explícitamente la interacción entre los perturbadores y el emisor. SimU utiliza trayectorias hiperbólicas (considerando el apantallamiento de Debye) mientras que SimUSP impone trayectorias rectas para aislar el efecto de la interacción emisor-perturbador.
• DinMol: Un código de dinámica molecular completo donde todos los electrones e iones interactúan mediante un complejo campo de Coulomb multibody, permitiendo la detección de eventos de recombinación (ensanchamiento por recombinación).

El estudio evalúa cuatro niveles de complejidad de interacción de partículas: partículas independientes, interacción con el radiador, interacción de $N$ cuerpos completa e interacción de $N$ cuerpos con ensanchamiento por recombinación. Además, los autores introducen un concepto novedoso: Términos de Interferencia Equivalentes a la Simulación (SEIT). Este enfoque aisla la contribución de los términos de interferencia (elementos fuera de la diagonal en el operador de impacto) dentro del marco de simulaciones computacionales completas mediante la omisión selectiva de términos específicos en la traza de la función de correlación de dipolo.

Resultados Clave

1. Consistencia de los Códigos: Los códigos con la misma física subyacente (trayectorias de partículas independientes) producen resultados idénticos, independientemente de las diferencias en la implementación numérica. SIMULAm, el código híbrido, muestra un acuerdo notable con los CSM completos (SIMULA, SimUSP), validando su uso para acelerar las simulaciones por un factor de $\sim$50.
2. Efecto de la Interacción de Partículas: Cuando se incluye la interacción con el radiador (SimU) o se utilizan dinámicas moleculares completas (DinMol), las líneas espectrales se vuelven más estrechas en comparación con los modelos de partículas independientes. Esto se atribuye a la repulsión de los iones perturbadores por parte del emisor cargado, lo que reduce el microcampo iónico. La inclusión del ensanchamiento por recombinación en DinMol contrarresta ligeramente este estrechamiento, resultando en un perfil intermedio entre SIMULA y SimU.
3. Dependencia de la Densidad: A densidades más bajas ($10^{24}$ cm$^{-3}$), los componentes principales de He-$\beta$ se resuelven. A densidades más altas ($10^{25}$ cm$^{-3}$), la mezcla de Stark se intensifica, causando que las transiciones prohibidas (por ejemplo, $1D_2 \to 1S_0$ y $3P_2 \to 1S_0$) dominen el espectro, mientras que las líneas no perturbadas contribuyen solo a las alas. Todos los modelos siguen la escala esperada de $n_e^{2/3}$ para los anchos de línea.
4. Términos de Interferencia (SEIT): El estudio encuentra que los SEIT tienen un efecto despreciable en la forma de los satélites de $n=2$. Sin embargo, para los satélites de $n=3$, particularmente en la densidad más alta ($10^{25}$ cm$^{-3}$), los SEIT modifican significativamente la forma de la línea. Eliminar los SEIT hace que la emisión de la satélite pierda características estructurales (como hombros en el ala de alta energía) y resulte en un ensanchamiento. Esto contradice los hallazgos recientes de Gallardo-Diaz et al. [21], que sugerían que los términos de interferencia eran menores. Los autores atribuyen esta discrepancia al hecho de que, en el estudio anterior, el ensanchamiento por impacto electrónico (donde actúan los términos de interferencia) era una contribución menor en comparación con el ensanchamiento iónico y la estructura atómica; bajo condiciones donde el ensanchamiento electrónico está aislado o es dominante, los términos de interferencia son críticos.
5. Espectros Sintéticos: Al combinarse con la cinética atómica (ABAko) y el ensanchamiento Doppler térmico, los resultados de la simulación producen espectros sintéticos que difieren de las predicciones de la "teoría estándar" previa (código MERL), particularmente en el ala de alta energía de la línea He-$\beta$, debido a la inclusión de la dinámica iónica.

Significancia y Reivindicaciones
El artículo afirma proporcionar los primeros resultados de simulación computacional para la emisión de satélites tipo Li de la línea compleja de Kr He-$\beta$ bajo condiciones relevantes de ICF. El trabajo demuestra la robustez de los códigos de simulación actuales y valida el uso de enfoques híbridos (como SIMULAm) para diagnósticos eficientes. Una contribución primaria es la definición y aplicación de SEIT, que clarifica el papel de los términos de interferencia en las transiciones de satélites complejas. Los autores concluyen que, si bien omitir los términos de interferencia puede ser aceptable para los satélites de $n=2$ o cuando el ensanchamiento electrónico es despreciable, esto conduce a inexactitudes significativas para los satélites de $n=3$ a densidades altas. El estudio enfatiza que los diagnósticos futuros deben tener en cuenta estos efectos para evitar errores en la inferencia de densidad, particularmente a medida que las condiciones experimentales se mueven hacia densidades más altas o involucran transiciones de números cuánticos principales más altos. El trabajo también destaca la necesidad de ir más allá de la aproximación de dipolo en estudios futuros para capturar los desplazamientos de polarización del plasma observados en experimentos.