Abstract:An accurate treatment of the multiphase interstellar medium (ISM) in hydrodynamic galaxy simulations requires that we follow not only the thermal evolution of the gas, but also the evolution of its chemical state, including its molecular chemistry, without assuming chemical (including ionisation) equilibrium. We present a reaction network that can be used to solve for this thermo-chemical evolution. Our model follows the evolution of all ionisation states of the 11 elements that dominate the cooling rate, along with important molecules such as H2 and CO, and the intermediate molecular species that are involved in their formation (20 molecules in total). We include chemical reactions on dust grains, thermal processes involving dust, cosmic ray ionisation and heating and photochemical reactions. We focus on conditions typical for the diffuse ISM, with densities of 10^-2 cm^-3 < nH < 10^4 cm^-3 and temperatures of 10^2 K < T < 10^4 K, and we consider a range of radiation fields, including no UV radiation. In this paper we consider only gas that is optically thin, while paper II considers gas that becomes shielded from the radiation field. We verify the accuracy of our model by comparing chemical abundances and cooling functions in chemical equilibrium with the photoionisation code Cloudy. We identify the major coolants in diffuse interstellar gas to be CII, SiII and FeII, along with OI and H2 at densities nH > 10^2 cm^-3. Finally, we investigate the impact of non-equilibrium chemistry on the cooling functions of isochorically or isobarically cooling gas. We find that, at T < 10^4 K, recombination lags increase the electron abundance above its equilibrium value at a given temperature, which can enhance the cooling rate by up to two orders of magnitude. The cooling gas also shows lower H2 abundances than in equilibrium, by up to an order of magnitude.

A 3D Monte Carlo Photoionization Code for Modeling Diffuse Ionized Gas

Ionization Modeling Astrophysical Gaseous Structures. I. The Optically Thin Regime

Fervent: Chemistry-coupled, ionising and non-ionising radiative feedback in magnetohydrodynamical simulations

Three-dimensional photoionization modelling of the planetary nebula NGC 3918

PHOTOPiC: Calculate photo-ionization functions and model coefficients for gas discharge simulations

Spatially-Resolved Temperature and Density Structures of Nearby HII regions

Magnetohydrodynamic Evolution of H ii Regions in Molecular Clouds: Simulation Methodology, Tests, and Uniform Media

Photochemistry and Heating/Cooling of the Multiphase Interstellar Medium with UV Radiative Transfer for Magnetohydrodynamic Simulations

Analyzing the Abundance Discrepancy Problem in HII Regions with Photoionization Modeling

A 3D Treatment of Radiative Transfer: Including Multi-wavelength Scattering; Absorption and Far Infrared Emission; and Arbitrary Geometry

Efficient Photo-heating Algorithms in Time-dependent Photo-ionization Simulations

Three-Dimensional Pic-Mc Modeling for Relativistic Electron Beam Transport Through Dense Plasma

Observing gas and dust in simulations of star formation with Monte Carlo radiation transport on Voronoi meshes

Monte Carlo Simulation of Lyman Alpha Scattering and Application to Damped Lyman Alpha Systems

Physical Processes in Star-Gas Systems

Simulating ionization feedback from young massive stars: impact of numerical resolution

$κ$monty: a Monte Carlo Compton Scattering code including non-thermal electrons

A Monte Carlo simulation of ion transport at finite temperatures

Non-equilibrium chemistry and cooling in the diffuse interstellar medium - I. Optically thin regime

Edge Theorem for Multivariable Systems

Photoionisation Feedback in a Self-Gravitating, Magnetised, Turbulent Cloud