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 using a set of FIRE-2 cosmological zoom-in simulations of Milky Way-mass disk galaxies, in which a multiphase ISM is self-consistently shaped by gravity, cooling, and stellar feedback. We analyze 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 midplane. Bulk flows (e.g., inflows and fountains) are important at a few disk scale heights, while thermal pressure from hot gas dominates at larger heights. Overall, the total midplane pressure is well-predicted by the weight of the disk gas, and we show that it also scales linearly with the star formation rate surface density (Sigma_SFR). These results support the notion that the Kennicutt-Schmidt relation arises because Sigma_SFR and the gas surface density (Sigma_g) are connected via the ISM midplane pressure.

We present an analysis of the physical and dynamical states of two sets of EAGLE zoom simulations of galaxy haloes, one at high redshift ($z=2-3$) and the other at low redshift ($z=0$), with masses of $\approx 10^{12} M_{\odot}$. Our focus is how the circumgalactic medium (CGM) of these $L^*$ star-forming galaxies change over the last 10 Gyr. We find that the high-$z$ CGM is almost equally divided between the "cool" ($T<10^5$ K) and "hot" ($T\geq 10^5$ K) phases, while the low-$z$ hot CGM phase contains $5\times$ more mass. The high-$z$ hot CGM contains 60% more metals than the cool CGM, while the low-$z$ cool CGM contains 35% more metals than the hot CGM content. The metals are evenly distributed radially between the hot and cool phases throughout the high-$z$ CGM. At high $z$, the CGM volume is dominated by hot outflows, cool gas is mainly inflowing, but cool metals are flowing outward. At low $z$, the cool metals dominate the interior and the hot metals are more prevalent at larger radii. The low-$z$ cool CGM has tangential motions consistent with rotational support out to $0.2 R_{200}$, often exhibiting $r \approx 40$ kpc disc-like structures. The low-$z$ hot CGM has several times greater angular momentum than the cool CGM, and a more flattened radial density profile than the high-$z$ hot CGM. This study verifies that, just as galaxies demonstrate significant evolutionary stages over cosmic time, the gaseous haloes surrounding them also undergo considerable changes of their own both in physical characteristics of density, temperature and metallicity, and dynamic properties of velocity and angular momentum.

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. Commonly used choices for the shielding length, such as the Jeans length, result in CO abundances that are too high at a given H$_2$ abundance. We find that a model with a distribution of shielding lengths, which has a median shielding length of $\sim 3$ pc in cold gas ($T < 300$ K) for both CO and H$_{2}$, is able to reproduce both the observed CO(1-0) luminosity and inferred CO-to-H$_{2}$ 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.

We recently used hydro-chemical simulations to demonstrate that molecular outflows observed in luminous quasars can be explained by molecule formation within the AGN wind. However, these simulations cover a limited parameter space, due to their computational cost. We have therefore developed an analytic model to follow cooling in the shocked ISM layer of an AGN wind. We explore different ambient densities ($1-10^{4} \, \rm{cm}^{-3}$), density profile slopes ($0-1.5$), AGN luminosities ($10^{44}-10^{47} \, \rm{erg} \, \rm{s}^{-1}$), and metallicities ($0.1-3 \rm{Z}_{\odot}$). The swept up gas mostly cools within ~1 Myr. Based on our previous simulations, we predict that this gas would produce observable molecular outflows. The instantaneous momentum boost initially increases as the outflow decelerates. However, it reaches a maximum of $\approx$20, due to work done against the gravitational potential. The predicted time-averaged observational estimate of the molecular outflow momentum boost reaches a maximum of $\approx1-2$, partly due to our assumed molecular fraction, 0.2, but also because the instantaneous and observational, time-averaged definitions are not equivalent. Thus recent observational estimates of order unity momentum boosts do not necessarily rule out energy-driven outflows. Finally, we find that dust grains are likely to re-form by accretion of metals after the shocked ISM layer has cooled, assuming that a small fraction of dust grains swept up after this layer has cooled are able to mix into the cool phase, and assuming that grain growth remains efficient in the presence of the strong AGN radiation field. This would enable rapid molecule formation, as assumed in our models.

Proximity zone fossils (PZFs) are ionization signatures around recently active galactic nuclei (AGN) where metal species in the circumgalactic medium remain over-ionized after the AGN has shut-off due to their long recombination timescales. We explore cosmological zoom hydrodynamic simulations using the EAGLE model paired with a non-equilibrium ionization and cooling module including time-variable AGN radiation to model PZFs around star-forming, disk galaxies in the z~0.2 Universe. Previous simulations typically under-estimated the O VI content of galactic haloes, but we show that plausible PZF models increase O VI column densities by 2-3x to achieve the levels observed around COS-Halos star-forming galaxies out to 150 kpc. Models with AGN bolometric luminosities >~10^43.6 erg s^-1, duty cycle fractions <~10%, and AGN lifetimes <~10^6 yr are the most promising, because their super-massive black holes grow at the cosmologically expected rate and they mostly appear as inactive AGN, consistent with COS-Halos. The central requirement is that the typical star-forming galaxy hosted an active AGN within a timescale comparable to the recombination time of a high metal ion, which for circumgalactic O VI is 10^7 years. H I, by contrast, returns to equilibrium much more rapidly due to its low neutral fraction and does not show a significant PZF effect. O VI absorption features originating from PZFs appear narrow, indicating photo-ionization, and are often well-aligned with lower metal ion species. PZFs are highly likely to affect the physical interpretation of circumgalactic high ionization metal lines if, as expected, normal galaxies host flickering AGN.

We explore the origin of fast molecular outflows that have been observed in Active Galactic Nuclei (AGN). Previous numerical studies have shown that it is difficult to create such an outflow by accelerating existing molecular clouds in the host galaxy, as the clouds will be destroyed before they can reach the high velocities that are observed. In this work, we consider an alternative scenario where molecules form in-situ within the AGN outflow. We present a series of hydro-chemical simulations of an isotropic AGN wind interacting with a uniform medium. We follow the time-dependent chemistry of 157 species, including 20 molecules, to determine whether molecules can form rapidly enough to produce the observed molecular outflows. We find H$_2$ outflow rates up to 140 M$_\odot$ yr$^{-1}$, which is sensitive to density, AGN luminosity, and metallicity. We compute emission and absorption lines of CO, OH and warm (a few hundred K) H$_2$ from the simulations in post-processing. The CO-derived outflow rates and OH absorption strengths at solar metallicity agree with observations, although the maximum line of sight velocities from the model CO spectra are a factor $\approx$2 lower than is observed. We derive a CO (1-0) to H$_2$ conversion factor of $\alpha_{\rm{CO} (1-0)}$ = 0.13 M$_\odot$ (K km s$^{-1}$ pc$^2$)$^{-1}$, 6 times lower than is commonly assumed in observations of such systems. We find strong emission from the mid-infrared lines of H$_2$. The mass of H$_2$ traced by this infrared emission is within a few per cent of the total H$_2$ mass. This H$_2$ emission may be observable by JWST.

We explore the circumgalactic metal content traced by commonly observed low ion absorbers, including C II, Si II, Si III, Si IV, and Mg II. We use a set of cosmological hydrodynamical zoom simulations run with the EAGLE model and including a non-equilibrium ionization and cooling module that follows 136 ions. The simulations of z~0.2 L* (M_200=10^11.7-10^12.3 Msol) haloes hosting star-forming galaxies and group-sized (M_200=10^12.7-10^13.3 Msol) haloes hosting mainly passive galaxies reproduce key trends observed by the COS-Halos survey-- low ion column densities show 1) little dependence on galaxy specific star formation rate, 2) a patchy covering fraction indicative of 10^4 K clouds with a small volume filling factor, and 3) a declining covering fraction as impact parameter increases from 20-160 kpc. Simulated Si II, Si III, Si IV, C II, and C III column densities show good agreement with observations, while Mg II is under-predicted. Low ions trace a significant metal reservoir, ~10^8 Msol, residing primarily at 10-100 kpc from star-forming and passive central galaxies. These clouds tend to flow inwards and most will accrete onto the central galaxy within the next several Gyr, while a small fraction are entrained in strong outflows. A two-phase structure describes the inner CGM (<0.5 R_200) with low-ion metal clouds surrounded by a hot, ambient medium. This cool phase is separate from the O VI observed by COS-Halos, which arises from the outer CGM (>0.5 R_200) tracing virial temperature gas around L* galaxies. Physical parameters derived from standard photo-ionization modelling of observed column densities (e.g. aligned Si II/Si III absorbers) are validated against our simulations. Our simulations therefore support previous ionization models indicating that cloud covering factors decline while densities and pressures show little variation with increasing impact parameter.

We study the effect of a fluctuating active galactic nucleus (AGN) on the abundance of circumgalactic OVI in galaxies selected from the EAGLE simulations. We follow the time-variable OVI abundance in post-processing around four galaxies - two at $z=0.1$ with stellar masses of $M_{\ast} \sim 10^{10}$ M$_{\odot}$ and $M_{\ast} \sim 10^{11}$ M$_{\odot}$, and two at $z=3$ with similar stellar masses - out to impact parameters of twice their virial radii, implementing a fluctuating central source of ionizing radiation. Due to delayed recombination, the AGN leave significant `AGN proximity zone fossils' around all four galaxies, where OVI and other metal ions are out of ionization equilibrium for several megayears after the AGN fade. The column density of OVI is typically enhanced by $\approx 0.3-1.0$ dex at impact parameters within $0.3R_{\rm vir}$, and by $\approx 0.06-0.2$ dex at $2R_{\rm vir}$, thereby also enhancing the covering fraction of OVI above a given column density threshold. The fossil effect tends to increase with increasing AGN luminosity, and towards shorter AGN lifetimes and larger AGN duty cycle fractions. In the limit of short AGN lifetimes, the effect converges to that of a continuous AGN with a luminosity of $(f_{\rm duty}/100\%)$ times the AGN luminosity. We also find significant fossil effects for other metal ions, where low-ionization state ions are decreased (SiIV, CIV at $z=3$) and high-ionization state ions are increased (CIV at $z=0.1$, NeVIII, MgX). Using observationally motivated AGN parameters, we predict AGN proximity zone fossils to be ubiquitous around $M_{\ast} \sim 10^{10-11}$ M$_{\odot}$ galaxies, and to affect observations of metals in the circumgalactic medium at both low and high redshifts.

We introduce a series of 20 cosmological hydrodynamical simulations of Lstar (M_200 =10^11.7 - 10^12.3 Msol) and group-sized (M_200 = 10^12.7 - 10^13.3 Msol) haloes run with the model used for the EAGLE project, which additionally includes a non-equilibrium ionization and cooling module that follows 136 ions. The simulations reproduce the observed correlation, revealed by COS-Halos at z~0.2, between O VI column density at impact parameters b < 150 kpc and the specific star formation rate (sSFR=SFR/Mstar) of the central galaxy at z~0.2. We find that the column density of circumgalactic O VI is maximal in the haloes associated with Lstar galaxies, because their virial temperatures are close to the temperature at which the ionization fraction of O VI peaks (T~10^5.5 K). The higher virial temperature of group haloes (> 10^6 K) promotes oxygen to higher ionization states, suppressing the O VI column density. The observed NO VI-sSFR correlation therefore does not imply a causal link, but reflects the changing characteristic ionization state of oxygen as halo mass is increased. In spite of the mass-dependence of the oxygen ionization state, the most abundant circumgalactic oxygen ion in both Lstar and group haloes is O VII; O VI accounts for only 0.1% of the oxygen in group haloes and 0.9-1.3% with Lstar haloes. Nonetheless, the metals traced by O VI absorbers represent a fossil record of the feedback history of galaxies over a Hubble time; their characteristic epoch of ejection corresponds to z > 1 and much of the ejected metal mass resides beyond the virial radius of galaxies. For both Lstar and group galaxies, more of the oxygen produced and released by stars resides in the circumgalactic medium (within twice the virial radius) than in the stars and ISM of the galaxy.

We present an analysis of Giant Molecular Clouds (GMCs) within hydrodynamic simulations of isolated, low-mass (M* ~ 10^9 M_sol) disc galaxies. We study 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 CO intensity threshold of 0.25 K km s^-1. Our simulations include a non-equilibrium chemical model with 157 species, including 20 molecules. We also investigate the effects of resolution and pressure floors (i.e. Jeans limiters). We find cloud lifetimes up to ~40 Myr, with a median of 13 Myr, in agreement with observations. At one tenth solar metallicity, young clouds (<10-15 Myr) are 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). We also compute CO emission from individual clouds. 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. The I_CO - A_v relation shifts towards higher A_v for higher metallicities and, to a lesser extent, for stronger UV radiation. At one tenth solar metallicity, CO emission is weaker in young clouds (<10-15 Myr), consistent with the underabundance of CO. Consequently, X_CO decreases by an order of magnitude from 0 to 15 Myr, albeit with a large scatter.

We present a series of hydrodynamic simulations of isolated galaxies with stellar mass of $10^{9} \, \rm{M}_{\odot}$. The models use a resolution of $750 \, \rm{M}_{\odot}$ per particle and include a treatment for the full non-equilibrium chemical evolution of ions and molecules (157 species in total), along with gas cooling rates computed self-consistently using the non-equilibrium abundances. We compare these to simulations evolved using cooling rates calculated assuming chemical (including ionisation) equilibrium, and we consider a wide range of metallicities and UV radiation fields, including a local prescription for self-shielding by gas and dust. We find higher star formation rates and stronger outflows at higher metallicity and for weaker radiation fields, as gas can more easily cool to a cold (few hundred Kelvin) star forming phase under such conditions. Contrary to variations in the metallicity and the radiation field, non-equilibrium chemistry generally has no strong effect on the total star formation rates or outflow properties. However, it is important for modelling molecular outflows. For example, the mass of H$_{2}$ outflowing with velocities $> 50 \, \rm{km} \, \rm{s}^{-1}$ is enhanced by a factor $\sim 20$ in non-equilibrium. We also compute the observable line emission from CII and CO. Both are stronger at higher metallicity, while CII and CO emission are higher for stronger and weaker radiation fields respectively. We find that CII is generally unaffected by non-equilibrium chemistry. However, emission from CO varies by a factor of $\sim 2 - 4$. This has implications for the mean $X_{\rm{CO}}$ conversion factor between CO emission and H$_{2}$ column density, which we find is lowered by up to a factor $\sim 2.3$ in non-equilibrium, and for the fraction of CO-dark molecular gas.

We extend the non-equilibrium model for the chemical and thermal evolution of diffuse interstellar gas presented in Richings et al. (2014) to account for shielding from the UV radiation field. We attenuate the photochemical rates by dust and by gas, including absorption by HI, H2, HeI, HeII and CO where appropriate. We then use this model to investigate the dominant cooling and heating processes in interstellar gas as it becomes shielded from the UV radiation. We consider a one-dimensional plane-parallel slab of gas irradiated by the interstellar radiation field, either at constant density and temperature or in thermal and pressure equilibrium. The dominant thermal processes tend to form three distinct regions in the clouds. At low column densities cooling is dominated by ionised metals such as SiII, FeII, FeIII and CII, which are balanced by photoheating, primarily from HI. Once the hydrogen-ionising radiation becomes attenuated by neutral hydrogen, photoelectric dust heating dominates, while CII becomes dominant for cooling. Finally, dust shielding triggers the formation of CO and suppresses photoelectric heating. The dominant coolants in this fully shielded region are H2 and CO. The column density of the HI-H2 transition predicted by our model is lower at higher density (or at higher pressure for gas clouds in pressure equilibrium) and at higher metallicity, in agreement with previous PDR models. We also compare the HI-H2 transition in our model to two prescriptions for molecular hydrogen formation that have been implemented in hydrodynamic simulations.

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.

Recent results indicate a correlation between nuclear radio-loudness of active galaxies and their central stellar surface-brightness profiles, in that `core' galaxies (with inner logarithmic slope {\gamma}<0.3) are significantly more radio loud than `power-law' galaxies ({\gamma}>0.5). This connection, which indicates possible links between radio-loudness and galaxy formation history (e.g. through black hole spin) has so far only been confirmed for a radio-selected sample of galaxies. Furthermore, it has since been shown that the Nuker law, which was used to parameterise the brightness profiles in these studies, gives a poor description of the brightness profile. Here, we present an analysis of the central surface brightness profiles of the active galaxies of Hubble Type T<3, that were identified by the optically-selected Palomar spectroscopic survey of nearby galaxies. We fit the brightness profiles using Sersic, Core-Sersic and, where necessary, Double-Sersic models, which we fit to the semi-major axis brightness profiles extracted from high resolution images of the galaxies from the Hubble Space Telescope (HST). We use these fits to classify the galaxies as `Core', `Sersic' or `Double-Sersic' and compare this classification with the properties of the Active Galactic Nuclei (AGNs). We find that AGN hosted in Core galaxies are generally more radio-loud than those hosted in Sersic galaxies, although there is a large overlap between the two subsamples. The correlation between radio-loudness and brightness profile can partly be explained by a correlation between radio-loudness and black hole mass. Additionally, there is a significant (99 per cent confidence) partial correlation between radio-loudness and the Core/Sersic classification of the host galaxy, which lends support to the previous results based on the radio-selected sample.

We report a detailed spectral analysis of the population of low-mass X-ray binaries (LMXBs) detected in the elliptical galaxy NGC~4278 with Chandra. Seven luminous sources were studied individually, four in globular clusters (GCs), and three in the stellar field. The range of (0.3-8 keV) $L_X$ for these sources suggests that they may be black hole binaries (BHBs). Comparison of our results with simulations allows us to discriminate between disk and power-law dominated emission, pointing to spectral/luminosity variability, reminiscent of Galactic BHBs. The BH masses derived from a comparison of our spectral results with the $L_X \sim T^4_{in}$ relation of Galactic BHBs are in the 5-15 $M_{\odot}$ range, as observed in the Milky Way. The analysis of joint spectra of sources selected in three luminosity ranges suggests that while the high luminosity sources have prominent thermal disk emission components, power-law components are likely to be important in the mid and low-luminosity spectra. Comparing low-luminosity average spectra, we find a relatively larger $N_H$ in the GC spectrum; we speculate that this may point to either a metallicity effect, or to intrinsic physical differences between field and GC accreting binaries. Analysis of average sample properties uncover a previously unreported $L_X - R_G$ correlation (where $R_G$ is the galactocentric radius) in the GC-LMXB sample, implying richer LMXB populations in more central GCs. No such trend is seen in the field LMXB sample. We can exclude that the GC $L_X - R_G$ correlation is the by-product of a luminosity effect, and suggest that it may be related to the presence of more compact GCs at smaller galactocentric radii, fostering more efficient binary formation.

We present an analysis of the hot interstellar medium (ISM) in the spiral galaxy NGC 4490, which is interacting with the irregular galaxy NGC 4485, using ~100ks of Chandra ACIS-S observations. The high angular resolution of Chandra enables us to remove discrete sources and perform spatially resolved spectroscopy for the star forming regions and associated outflows, allowing us to look at how the physical properties of the hot ISM such as temperature, hydrogen column density and metal abundances vary throughout these galaxies. We find temperatures of >0.41 keV and 0.85 +0.59/-0.12 keV, electron densities of >1.87 eta^(-1/2) x 10^(-3) cm^(-3) and 0.21 +0.03/-0.04 eta^(-1/2) x 10^(-3) cm^(-3), and hot gas masses of >1.1 eta^(1/2) x 10^7 M_{\odot} and ~3.7 eta^(1/2) x 10^7 M_{\odot} in the plane and halo of NGC 4490 respectively, where eta is the filling factor of the hot gas. The abundance ratios of Ne, Mg and Si with respect to Fe are found to be consistent with those predicted by theoretical models of type II supernovae. The thermal energy in the hot ISM is ~5% of the total mechanical energy input from supernovae, so it is likely that the hot ISM has been enriched and heated by type II supernovae. The X-ray emission is anticorrelated with the H-alpha and mid-infrared emission, suggesting that the hot gas is bounded by filaments of cooler ionized hydrogen mixed with warm dust.