Tracking electron capture processes in classical molecular dynamics simulations for spectral line broadening in plasmas

Technical Summary: Tracking Electron Capture in Classical Molecular Dynamics for Spectral Line Broadening

Problem Statement
Accurate modeling of spectral line shapes is essential for diagnosing high-energy-density plasmas, particularly regarding Stark broadening. While computer simulations, specifically Full Molecular Dynamics (FMD), excel at replicating particle interactions and their impact on radiating atoms, they face a significant challenge in tracking electron capture (recombination) processes. In FMD simulations, distinguishing between a free electron and one that has become bound to an ion is non-trivial. Standard criteria often rely on instantaneous conditions, such as the electron-ion distance being below a threshold $a$ and the total energy being negative. However, for ions with charge $Z > 1$ and in strongly coupled plasmas, these instantaneous criteria fail. They can misidentify transient close encounters as bound states or miss bound states where electrons orbit at distances larger than the threshold, leading to unphysical oscillating electric fields in the simulation data. These errors corrupt the time-histories of the electric microfield required for calculating spectral line shapes and result in an overestimation of recombination broadening.

Methodology
The authors propose a new algorithm designed to operate within a classical FMD framework to precisely identify when an electron is captured by an ion and follows a stable orbit. The methodology introduces a global criterion that relies on the time-history of the system rather than instantaneous snapshots.

Key components of the algorithm include:

1. Potential Regularization: To handle divergences at short distances, the ion-electron interaction potential is regularized using a parabolic form for distances $r \leq a$ and a hyperbolic (Coulombic) form for $r > a$. The parameter $a$ is defined based on the ionization energy $V_i$ of the $Z-1$ charge state.
2. Time-Dependent Tracking: For each ion, the simulation tracks the $Z+1$ closest electrons at every timestep. The algorithm monitors the potential energy of the emitter ($E_p$). When $E_p$ drops below a threshold $E_{th} = -V_i$, the identity of the closest electron is recorded.
3. Stability Criterion: The algorithm checks if the identified electron remains among the $Z+1$ closest neighbors for a duration exceeding a threshold time $\tau_{bound}$. This threshold is set to $3\tau_T$, where $\tau_T$ is the characteristic orbital period of an electron at thermal velocity.
4. Energy Verification: If the time condition is met, the average total energy of the electron-ion pair over that interval is calculated. If the energy is negative, the electron is classified as bound.
5. Field Sequence Handling: Once an electron is deemed captured, the electric field sequence for that radiator is terminated (as the ion's charge state changes from $+Ze$ to $+(Z-1)e$), and the correlation function is set to zero. This prevents the inclusion of non-physical, high-frequency oscillations caused by the orbiting electron in the subsequent line shape calculations.

Key Contributions

• Algorithm for $Z \geq 1$: The work extends previous criteria (valid for $Z=1$) to ions with charges $Z \geq 1$, addressing the complexity of multiple bound electrons.
• Removal of Unphysical Artifacts: By incorporating a time-history check, the algorithm filters out transient collisions that would otherwise be mistaken for bound states, thereby eliminating unphysical oscillating electric field sequences.
• Ionization Balance Estimation: The method provides a mechanism to estimate the fractional population of different ion charge states within the plasma, offering a way to calculate the ionization balance directly from the simulation data.

Results and Validation
The authors validated the algorithm through two primary approaches:

1. Internal Consistency: The ionization balance derived from the new electron-capture detection algorithm was compared against a method based on the statistical analysis of the ion potential energy distribution. The results showed excellent agreement across a range of plasma temperatures and coupling parameters.
2. Benchmarking: The simulation results were compared with atomic kinetic calculations using the collisional-radiative code FLYCHK (which utilizes the quantum Saha equation). The classical simulation results agreed with FLYCHK within a factor of $\sim 2$, a level of agreement considered standard among different atomic kinetic codes. This validates the choice of the potential well depth parameter $V_b$ and the overall approach.

Significance and Claims
The paper claims that this detection algorithm is a pivotal tool for advancing the accuracy of spectral line broadening simulations in complex plasma environments. By correctly identifying electron capture events, the method allows for the inclusion of recombination broadening—a mechanism often overlooked in standard FMD simulations—without introducing non-physical artifacts. The authors emphasize that without such a methodology, particularly for strongly coupled plasmas and ions with $Z > 1$, line profile calculations would be inaccurate. The work establishes a robust, self-consistent framework for tracking ionization-recombination mechanisms in classical molecular dynamics, enabling more precise modeling of Stark-broadened lines. The authors note limitations regarding the applicability to heavily coupled systems (e.g., solid-density plasmas) where the regularization radius might exceed the inter-particle distance, and they do not claim to fully account for radiative recombination, noting their model represents a combined effect of electron capture and likely three-body recombination.