Partially ionized two-fluid shocks with collisional and radiative ionization and recombination - multilevel hydrogen model

Explosive phenomena are known to trigger a wealth of shocks in warm plasma environments, including the solar chromosphere and molecular clouds where the medium consists of both ionized and neutral species. Partial ionization is critical in determining the behaviour of shocks, since the ions and neutrals locally decouple, allowing for substructure to exist within the shock. Accurately modelling partially ionized shocks requires careful treatment of the ionized and neutral species, and their interactions. Here we study a partially ionized switch-off slow-mode shock using a multilevel hydrogen model with both collisional and radiative ionization and recombination rates that are implemented into the two-fluid (PIP) code, and study physical parameters that are typical of the solar chromosphere. The multilevel hydrogen model differs significantly from magnetohydrodynamic (MHD) solutions due to the macroscopic thermal energy loss during collisional ionization. In particular, the plasma temperature both post-shock and within the finite-width is significantly cooler that the post-shock MHD temperature. Furthermore, in the mid to lower chromosphere, shocks feature far greater compression than their single-fluid MHD analogues. The decreased temperature and increased compression reveal the importance of non-equilibrium ionized in the thermal evolution of shocks in partially ionized media. Since partially ionized shocks are not accurately described by the Rankine-Hugoniot shock jump conditions, it may be incorrect to use these to infer properties of lower atmospheric shocks.