Photodissociation regions (PDRs) are ubiquitously present in the interstellar medium (ISM) and consist of predominantly neutral gas and dust illuminated by far-ultraviolet (FUV) radiation (6 < hv < 13.6 eV). Studies of PDRs allow us to understand the effects of FUV photons on the chemistry and structure of the neutral ISM in galaxies, as well as diagnosing the conditions within star-forming regions. PDRs are responsible for most of the infrared radiation from galaxies. Solving numerically the astrochemical network that occurs in PDRs is a powerful way to estimate the abundances of species and simulate the line emission of the observed galactic and extragalactic clouds.
I am the main developer of the astrochemical code 3D-PDR. This is an open-source OpenMP-parallelized 3D code. It is specifically designed to treat PDRs of arbitrary density distributions. It uses the Large Velocity Gradient approximation to compute the cooling functions, using a HEALPIX-based ray-tracing scheme. The code solves for the thermal balance self-consistently with the chemistry, to determine the gas temperature. It includes heating due to photoionisation and photodissociation reactions, in addition to the standard gas-phase chemistry and it accounts for self-shielding of H2 and CO against photodissociation. Suprathermal formation of CO via CH+ is included.