Chemical evolution of giant molecular clouds in simulations of galaxies

Alexander J. Richings, Joop Schaye

Published 2016-04-04Version 1

We present an analysis of Giant Molecular Clouds (GMCs) identified in hydrodynamic simulations of isolated, low-mass (M* ~ 10^9 M_sol) disc galaxies, with a particular focus on the evolution of molecular abundances and the implications for CO emission and the X_CO conversion factor in individual clouds. We define clouds either as regions above a density threshold n_H,min = 10 cm^-3, or using an observationally motivated velocity-integrated CO line intensity threshold of 0.25 K km s^-1. Our simulations include a non-equilibrium treatment for the chemistry of 157 species, including 20 molecules. We use a suite of runs to carefully investigate the effects of numerical resolution and pressure floors (i.e. Jeans mass limiters). We find cloud lifetimes up to ~40 Myr, with a median of 13 Myr, in agreement with observations. At ten per cent solar metallicity, young clouds (<10-15 Myr) tend to be underabundant in H2 and CO compared to chemical equilibrium, by factors of ~3 and 1-2 orders of magnitude, respectively. At solar metallicity, GMCs reach chemical equilibrium faster (within ~1 Myr), due to a higher formation rate of H2 on dust grains. We also compute CO J = 1-0 line emission from our simulated GMCs in post-processing. We find that the mean CO intensity, I_CO, is strongly suppressed at low dust extinction, A_v, and possibly saturates towards high A_v, in agreement with observations. Our simulated I_CO - A_v relation shifts towards higher A_v for higher metallicities and, to a lesser extent, for stronger UV radiation fields. At ten per cent solar metallicity, we find weaker CO emission in young clouds (<10-15 Myr), consistent with the underabundance of CO in such clouds. This is reflected in the median X_CO factor, which decreases by an order of magnitude from 0 to 15 Myr, albeit with a large scatter.