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

Resumo Técnico: Simulações Computacionais do Efeito Stark no Complexo de Hélio-$\beta$ do Criptônio em Condições de ICF

Definição do Problema
A espectroscopia de experimentos de fusão por confinamento inercial (ICF) depende de dopantes, como o criptônio, para diagnosticar a temperatura e a densidade de plasmas em implosão. Especificamente, as transições da camada K de espécies tipo He, particularmente a linha He-$\beta$ (transição entre $n=3$ e $n=1$), são críticas para medições de densidade. No entanto, a linha He-$\beta$ é acompanhada por transições satélite do tipo Li ($1sn'\ell'3\ell \to 1s^2n'\ell'$), que se sobrepõem parcialmente à linha principal. Modelar todo o complexo He-$\beta$ sob condições extremas de ICF (densidades eletrônicas $n_e = 10^{24}$ a $10^{25}$ cm$^{-3}$ e temperaturas $T_e = 3$ keV) é computacionalmente proibitivo devido ao grande número de níveis atômicos e elementos de matriz de dipolo não nulos envolvidos. Embora os modelos de simulação computacional (CSM) sejam as ferramentas mais avançadas para o cálculo de perfis Stark, sua aplicação a tais formas de linha complexas tem sido desafiadora. Além disso, estudos anteriores usando aproximações de "teoria padrão" (assumindo íons estáticos e elétrons de impacto) sugeriram que os termos de interferência no operador de impacto eletrônico eram negligenciáveis para essas transições, um achado que exigia verificação sob condições de simulação completa.

Metodologia
Os autores apresentam um estudo comparativo utilizando múltiplos códigos de simulação computacional semiclássicos para modelar o alargamento Stark da linha He-$\beta$ do Kr e suas satélites Li-like de $n=2$ e $n=3$. Os códigos empregados incluem:

• SIMULA e SIMULAm: Modelos onde as partículas seguem trajetórias retilíneas, tratando os perturbadores como independentes. O SIMULAm é um código híbrido que incorpora o operador de alargamento eletrônico da teoria padrão em cada passo de tempo para reduzir o custo computacional.
• SimU e SimUSP: Modelos que resolvem explicitamente a interação entre os perturbadores e o emissor. O SimU utiliza trajetórias hiperbólicas (levando em conta o blindagem de Debye) enquanto o SimUSP impõe caminhos retilíneos para isolar o efeito da interação emissor-perturbador.
• DinMol: Um código de dinâmica molecular completo onde todos os elétrons e íons interagem via um campo de Coulomb multibody complexo, permitindo a detecção de eventos de recombinação (alargamento por recombinação).

O estudo avalia quatro níveis de complexidade de interação de partículas: partículas independentes, interação com o radiador, interação $N$-body completa e interação $N$-body com alargamento por recombinação. Adicionalmente, os autores introduzem um conceito inovador: Termos de Interferência Equivalentes à Simulação (SEIT). Esta abordagem isola a contribuição dos termos de interferência (elementos fora da diagonal no operador de impacto) dentro do framework de simulações computacionais completas, ao omitir seletivamente termos específicos no traço da função de correlação de dipolo.

Principais Resultados

1. Consistência dos Códigos: Códigos com a mesma física subjacente (trajetórias de partículas independentes) produzem resultados idênticos, independentemente das diferenças de implementação numérica. O SIMULAm, o código híbrido, mostra um acordo notável com os CSMs completos (SIMULA, SimUSP), validando seu uso para acelerar as simulações em um fator de $\sim$50.
2. Efeito da Interação de Partículas: Quando a interação com o radiador é incluída (SimU) ou quando a dinâmica molecular completa é usada (DinMol), as linhas espectrais tornam-se mais estreitas em comparação com modelos de partículas independentes. Isso é atribuído à repulsão de íons perturbadores pelo emissor carregado, o que reduz o microcampo iônico. A inclusão do alargamento por recombinação no DinMol contrabalanceia levemente esse estreitamento, resultando em um perfil intermediário entre o SIMULA e o SimU.
3. Dependência de Densidade: Em densidades mais baixas ($10^{24}$ cm$^{-3}$), os componentes principais da He-$\beta$ são resolvidos. Em densidades mais altas ($10^{25}$ cm$^{-3}$), o mixing Stark se intensifica, fazendo com que transições proibidas (ex: $1D_2 \to 1S_0$ e $3P_2 \to 1S_0$) dominem o espectro, enquanto as linhas não perturbadas contribuem apenas para as asas. Todos os modelos seguem a escala esperada de $n_e^{2/3}$ para as larguras de linha.
4. Termos de Interferência (SEIT): O estudo descobre que os SEIT têm um efeito negligenciável na forma das satélites de $n=2$. No entanto, para as satélites de $n=3$, particularmente na densidade mais alta ($10^{25}$ cm$^{-3}$), os SEIT modificam significativamente a forma da linha. Remover os SEIT faz com que a emissão da satélite perca características estruturais (como ombros na asa de alta energia) e resulte em alargamento. Isso contradiz achados recentes de Gallardo-Diaz et al. [21], que sugeriram que os termos de interferência eram menores; os autores atribuem essa discrepância ao fato de que, no estudo anterior, o alargamento por impacto eletrônico (onde os termos de interferência atuam) era uma contribuição menor comparada ao alargamento iônico e à estrutura atômica; sob condições onde o alargamento eletrônico é isolado ou dominante, os termos de interferência são críticos.
5. Espectros Sintéticos: Quando combinados com cinética atômica (ABAko) e alargamento Doppler térmico, os resultados de simulação produzem espectros sintéticos que diferem das previsões da "teoria padrão" anteriores (código MERL), particularmente na asa de alta energia da linha He-$\beta$, devido à inclusão da dinâmica iônica.

Significância e Alegações
O artigo afirma fornecer os primeiros resultados de simulação computacional para a emissão de satélite Li-like do complexo de He-$\beta$ do Kr sob condições relevantes para ICF. O trabalho demonstra a robustez dos atuais códigos de simulação e valida o uso de abordagens híbridas (como o SIMULAm) para diagnósticos eficientes. Uma contribuição primária é a definição e aplicação de SEIT, que esclarece o papel dos termos de interferência em transições de satélite complexas. Os autores concluem que, embora omitir termos de interferência possa ser aceitável para satélites de $n=2$ ou quando o alargamento eletrônico é negligenciável, isso leva a imprecisões significativas para satélites de $n=3$ em altas densidades. O estudo enfatiza que diagnósticos futuros devem considerar esses efeitos para evitar erros na inferência de densidade, particularmente conforme as condições experimentais se movem para densidades mais altas ou envolvem transições de números quânticos principais mais elevados. O trabalho também destaca a necessidade de ir além da aproximação de dipolo em estudos futuros para capturar deslocamentos de polarização de plasma observados em experimentos.