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

Arun Mathew,Jonathan Mackey,Maggie Celeste,Thomas J. Haworth,Garrelt Mellema
2024-09-28
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.
Astrophysics of Galaxies,Computational Physics
What problem does this paper attempt to address?