We present models of CO(1-0) emission from Milky-Way-mass galaxies at redshift zero in the FIRE-2 cosmological zoom-in simulations. We calculate the molecular abundances by post-processing the simulations with an equilibrium chemistry solver while accounting for the effects of local sources, and determine the emergent CO(1-0) emission using a line radiative transfer code. We find that the results depend strongly on the shielding length assumed, which, in our models, sets the attenuation of the incident UV radiation field. At the resolution of these simulations, commonly used choices for the shielding length, such as the Jeans length, result in CO abundances that are too high at a given H2 abundance. We find that a model with a distribution of shielding lengths, which has a median shielding length of ~3 pc in cold gas (T < 300 K) for both CO and H2, is able to reproduce both the observed CO(1-0) luminosity and inferred CO-to-H2 conversion factor at a given star formation rate compared with observations. We suggest that this short shielding length can be thought of as a subgrid model, which controls the amount of radiation that penetrates giant molecular clouds.

Pressure balance plays a central role in models of the interstellar medium (ISM), but whether and how pressure balance is realized in a realistic multiphase ISM is not yet well understood. We address this question by using a set of FIRE-2 cosmological zoom-in simulations of Milky Way-mass disc galaxies, in which a multiphase ISM is self-consistently shaped by gravity, cooling, and stellar feedback. We analyse how gravity determines the vertical pressure profile as well as how the total ISM pressure is partitioned between different phases and components (thermal, dispersion/turbulence, and bulk flows). We show that, on average and consistent with previous more idealized simulations, the total ISM pressure balances the weight of the overlying gas. Deviations from vertical pressure balance increase with increasing galactocentric radius and with decreasing averaging scale. The different phases are in rough total pressure equilibrium with one another, but with large deviations from thermal pressure equilibrium owing to kinetic support in the cold and warm phases, which dominate the total pressure near the mid-plane. Bulk flows (e.g. inflows and fountains) are important at a few disc scale heights, while thermal pressure from hot gas dominates at larger heights. Overall, the total mid-plane pressure is well-predicted by the weight of the disc gas and we show that it also scales linearly with the star formation rate surface density (ςSFR). These results support the notion that the Kennicutt-Schmidt relation arises because ςSFR and the gas surface density (ςg) are connected via the ISM mid-plane pressure.