A multi-ion non-equilibrium solver for ionised astrophysical plasmas with arbitrary elemental abundances

TL;DR

This paper introduces a non-equilibrium multi-ion solver for astrophysical plasmas with arbitrary elemental abundances, enabling more accurate modeling of complex, dynamic ionized environments beyond equilibrium assumptions.

Contribution

It presents a novel multi-ion module for the extsc{pion} code that handles arbitrary elemental abundances and non-equilibrium processes, validated against existing models.

Findings

The module accurately reproduces equilibrium and non-equilibrium ionization states.

Charge exchange significantly affects ion abundances in cooling plasmas.

Application to nebulae demonstrates realistic spectral line predictions.

Abstract

While many astrophysical plasmas can be modelled successfully assuming ionisation and thermal equilibrium, in some cases this is not appropriate and a non-equilibrium approach is required. In nebulae around evolved stars the local elemental abundances may also strongly vary in space and time. Here we present a non-equilibrium multi-ion module developed for the fluid-dynamics code \textsc{pion}, describing the physical processes included and demonstrating its capabilities with some test calculations. A non-equilbrium ionisation solver is developed that allows arbitrary elemental abundances for neutral and ionised (but not molecular) gas, for the elements H, He, C, N, O, Ne, Si, S and Fe. Collisional ionisation and recombination, photoionisation and charge-exchange reactions are included, and ion-by-ion non-equilibrium radiative cooling is calculated based on the instantaneous ion fractions of each element. Element and ion mass-fractions are advected using passive scalars, operator-split from the microphysical processes. The module is validated by comparing with equilibrium and non-equilibrium calculations in the literature. Effects of charge exchange on ion abundances in cooling plasmas are discussed. Application to modelling shocks and photoionised H~\textsc{ii} regions is demonstrated. The time-dependent expansion of a Wolf-Rayet nebula is studied, including photoionisation and collisional processes, and spectral-line luminosities calculated for non-equilibrium and equilibrium plasma states. The multi-ion module enables simulation of ionised plasmas with spatially varying elemental abundances using self-consistent ion abundances and thermal evolution. This allows prediction of spectral lines in UV, optical, IR and X-ray even in cases where the plasma is out of ionisation equilibrium.