The multiphase circumgalactic medium traced by low metal ions in EAGLE zoom simulations

We explore the circumgalactic metal content traced by commonly observed low ion absorbers, including C II, SiII, SiIII, SiIV, and MgII. 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* (M200= 1011.7- 1012.3M⊙) haloes hosting star-forming galaxies and group-sized (M200= 1012.7- 1013.3M⊙) 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 104K clumps with a small volume filling factor, and 3) a declining covering fraction as impact parameter increases from 20-160kpc. Simulated Si II, Si III, Si IV, CII, and C III column densities show good (mostly within 0.3 dex) agreement with observations, while MgII is under-predicted. Low ions trace a significant metal reservoir, ≈108M⊙, residing primarily at 10-100kpc from star-forming and passive central galaxies. These clumps preferentially 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 multiphase structure describes the inner CGM ( 0.5R200) tracing virial temperature gas around L* galaxies. Our simulations support previous ionization models indicating that cloud covering factors decline while densities and pressures show little decline with increasing impact parameter (typically < 0.3 dex from 40 to 160 kpc). We find the cool clumps have lower pressures than the ambient medium they are embedded in, and discuss that numerical effects within the hydrodynamic solver likely play a role. © 2018 The Author(s)

The Search for the Missing Baryons at Low Redshift

At low redshift, only about one-tenth of the known baryons lie in galaxies or the hot gas seen in galaxy clusters and groups. Models posit that these "missing baryons" are in gaseous form in overdense filaments that connect the much denser virialized groups and clusters. About 30% are cool (<1E5 K) and are detected in Ly alpha absorption studies, but about half is predicted to lie in the 1E5-1E7 K regime. Gas is detected in the 2-5E5 K range through OVI absorption studies (7% of the baryons) and possibly near 1E5 K from broad Ly absorption (20% of the baryons). Hotter gas (0.5-2E6 K) is detected at zero redshift by OVII and OVIII K X-ray absorption, and the OVII line strengths seem to correlate with the Galactic soft X-ray background, so it is probably produced by Galactic Halo gas, rather than a Local Group medium. There are no compelling detections of the intergalactic hot gas (0.5-10E6 K) either in absorption or emission and these upper limits are consistent with theoretical models. Claimed X-ray absorption lines are not confirmed, while most of the claims of soft emission are attributable to artifacts of background subtraction and field-flattening. The missing baryons should become detectable with moderate improvements in instrumental sensitivity.Comment: To appear in Annual Review of Astronomy and Astrophysics, Vol 45 (Sept 2007) 44 pages, including 11 figure

The nature of submillimetre galaxies in cosmological hydrodynamic simulations

We study the nature of rapidly star-forming galaxies at z= 2 in cosmological hydrodynamic simulations, and compare their properties to observations of submillimetre galaxies (SMGs). We identify simulated SMGs as the most rapidly star-forming systems that match the observed number density of SMGs. In our models, SMGs are massive galaxies sitting at the centres of large potential wells, being fed by smooth infall and gas-rich satellites at rates comparable to their star formation rates (SFRs). They are not typically undergoing major mergers that significantly boost their quiescent SFR, but they still often show complex gas morphologies and kinematics. Our simulated SMGs have stellar masses of M*∼ 1011−11.7 M⊙, SFRs of ∼180–500 M⊙ yr−1, a clustering length of ∼10 h−1 Mpc and solar metallicities. The SFRs are lower than those inferred from far-infrared (far-IR) data by ∼×3, which we suggest may owe to one or more systematic effects in the SFR calibrations. SMGs at z= 2 live in ∼1013 M⊙ haloes, and by z= 0 they mostly end up as brightest group galaxies in ∼1014 M⊙ haloes. We predict that higher M* SMGs should have on average lower specific SFRs, less disturbed morphologies and higher clustering. We also predict that deeper far-IR surveys will smoothly join SMGs on to the massive end of the SFR–M* relationship defined by lower mass z∼ 2 galaxies. Overall, our simulated rapid star-formers provide as good a match to available SMG data as merger-based scenarios, offering an alternative scenario that emerges naturally from cosmological simulations

A fundamental problem in our understanding of low-mass galaxy evolution

Recent studies have found a dramatic difference between the observed number density evolution of low-mass galaxies and that predicted by semi-analytic models. Whilst models accurately reproduce the z= 0 number density, they require that the evolution occurs rapidly at early times, which is incompatible with the strong late evolution found in observational results. We report here the same discrepancy in two state-of-the-art cosmological hydrodynamical simulations, which is evidence that the problem is fundamental. We search for the underlying cause of this problem using two complementary methods. First, we consider a narrow range in stellar mass of log (Mstar/(h−2M_)) = 9–9.5 and look for evidence of a different history of today’s low-mass galaxies in models and observations. We find that the exclusion of satellite galaxies from the analysis brings the median ages and star formation rates of galaxies into reasonable agreement. However, the models yield too few young, strongly star-forming galaxies. Secondly, we construct a toy model to link the observed evolution of specific star formation rates with the evolution of the galaxy stellar mass function. We infer from this model that a key problem in both semi-analytic and hydrodynamical models is the presence of a positive instead of a negative correlation between specific star formation rate and stellar mass. A similar positive correlation is found between the specific dark matter halo accretion rate and the halo mass, indicating that model galaxies are growing in a way that follows the growth of their host haloes too closely. It therefore appears necessary to find a mechanism that decouples the growth of low-mass galaxies, which occurs primarily at late times, from the growth of their host haloes, which occurs primarily at early times. We argue that the current form of star formation-driven feedback implemented in most galaxy formation models is unlikely to achieve this goal, owing to its fundamental dependence on host halo mass and time

Feedback from supermassive black holes transforms centrals into passive galaxies by ejecting circumgalactic gas

Davies et al. (2019) established that for L^* galaxies the fraction of baryons in the circumgalactic medium (CGM) is inversely correlated with the mass of their central supermassive black holes (BHs) in the EAGLE hydrodynamic simulation. The interpretation is that, over time, a more massive BH has provided more energy to transport baryons beyond the virial radius, which additionally reduces gas accretion and star formation. We continue this research by focusing on the relationship between the 1) BH masses, 2) physical and observational properties of the CGM, and 3) galaxy colours for Milky Way-mass systems. The ratio of the cumulative BH feedback energy over the gaseous halo binding energy is a strong predictor of the CGM gas content, with BHs injecting >~10x the binding energy resulting in gas-poor haloes. Observable tracers of the CGM, including CIV, OVI, and HI absorption line measurements, are found to be effective tracers of the total z~0 CGM halo mass. We use high-cadence simulation outputs to demonstrate that BH feedback pushes baryons beyond the virial radius within 100 Myr timescales, but that CGM metal tracers take longer (0.5-2.5 Gyr) to respond. Secular evolution of galaxies results in blue, star-forming or red, passive populations depending on the cumulative feedback from BHs. The reddest quartile of galaxies with M_*=10^{10.2-10.7} M_solar (median u-r = 2.28) has a CGM mass that is 2.5x lower than the bluest quartile (u-r=1.59). We propose strategies for observing the predicted lower CGM column densities and covering fractions around galaxies hosting more massive BHs using the Cosmic Origins Spectrograph on Hubble

Simulating Groups and the IntraGroup Medium: The Surprisingly Complex and Rich Middle Ground between Clusters and Galaxies

Galaxy groups are more than an intermediate scale between clusters and halos hosting individual galaxies, they are crucial laboratories capable of testing a range of astrophysics from how galaxies form and evolve to large scale structure (LSS) statistics for cosmology. Cosmological hydrodynamic simulations of groups on various scales offer an unparalleled testing ground for astrophysical theories. Widely used cosmological simulations with ∼(100 Mpc)3 volumes contain statistical samples of groups that provide important tests of galaxy evolution influenced by environmental processes. Larger volumes capable of reproducing LSS while following the redistribution of baryons by cooling and feedback are essential tools necessary to constrain cosmological parameters. Higher resolution simulations can currently model satellite interactions, the processing of cool (T ∼ 104 K) multi-phase gas, and non-thermal physics including turbulence, magnetic fields, and cosmic ray transport. We review simulation results regarding the gas and stellar contents of groups, cooling flows and the relation to the central galaxy, the formation and processing of multi-phase gas, satellite interactions with the intragroup medium, and the impact of groups for cosmological parameter estimation. Cosmological simulations provide evolutionarily consistent predictions of these observationally difficult-to-define objects, and have untapped potential to accurately model their gaseous, stellar, and dark matter distributions

The abundance and physical properties of O VII and O VIII X-ray absorption systems in the EAGLE simulations

We use the EAGLE cosmological, hydrodynamical simulations to predict the column density and equivalent width distributions of intergalactic O VII ($E=574$ eV) and O VIII ($E=654$ eV) absorbers at low redshift. These two ions are predicted to account for 40% of the gas-phase oxygen, which implies that they are key tracers of cosmic metals. We find that their column density distributions evolve little at observable column densities from redshift 1 to 0, and that they are sensitive to AGN feedback, which strongly reduces the number of strong (column density $N \gtrsim 10^{16} \, \mathrm{cm}^{-2})$ absorbers. The distributions have a break at $N \sim 10^{16} \, \mathrm{cm}^{-2}$, corresponding to overdensities of $\sim 10^{2}$, likely caused by the transition from sheet/filament to halo gas. Absorption systems with $N \gtrsim 10^{16} \mathrm{cm}^{-2}$ are dominated by collisionally ionized O VII and O VIII, while the ionization state of oxygen at lower column densities is also influenced by photoionization. At these high column densities, O VII and O VIII arising in the same structures probe systematically different gas temperatures, meaning their line ratio does not translate into a simple estimate of temperature. While O VII and O VIII column densities and covering fractions correlate poorly with the H I column density at $N_{\mathrm{H \, I}} \gtrsim 10^{15} \, \mathrm{cm}^{-2}$, O VII and O VIII column densities are higher in this regime than at the more common, lower H I column densities. The column densities of O VI and especially Ne VIII, which have strong absorption lines in the UV, are good predictors of the strengths of O VII and O VIII absorption and can hence aid in the detection of the X-ray lines

Predictions for the X-ray circumgalactic medium of edge-on discs and spheroids

We investigate how the X-ray circumgalactic medium (CGM) of present-day galaxies depends on galaxy morphology and azimuthal angle using mock observations generated from the EAGLE cosmological hydrodynamic simulation. By creating mock stacks of eROSITA-observed galaxies oriented to be edge-on, we make several observationally-testable predictions for galaxies in the stellar mass range M⋆ = 1010.7 − 11.2 M⊙. The soft X-ray CGM of disc galaxies is between 60 and 100% brighter along the semi-major axis compared to the semi-minor axis, between 10-30 kpc. This azimuthal dependence is a consequence of the hot (T &amp;gt; 106 K) CGM being non-spherical: specifically it is flattened along the minor axis such that denser and more luminous gas resides in the disc plane and co-rotates with the galaxy. Outflows enrich and heat the CGM preferentially perpendicular to the disc, but we do not find an observationally-detectable signature along the semi-minor axis. Spheroidal galaxies have hotter CGMs than disc galaxies related to spheroids residing at higher halos masses, which may be measurable through hardness ratios spanning the 0.2 − 1.5 keV band. While spheroids appear to have brighter CGMs than discs for the selected fixed M⋆ bin, this owes to spheroids having higher stellar and halo masses within that M⋆ bin, and obscures the fact that both simulated populations have similar total CGM luminosities at the exact same M⋆. Discs have brighter emission inside 20 kpc and more steeply declining profiles with radius than spheroids. We predict that the eROSITA 4-year all-sky survey should detect many of the signatures we predict here, although targeted follow-up observations of highly inclined nearby discs after the survey may be necessary to observe some of our azimuthally-dependent predictions

EAGLE and Illustris-TNG Predictions for Resolved eROSITA X-Ray Observations of the Circumgalactic Medium around Normal Galaxies

We simulate stacked observations of nearby hot X-ray coronae associated with galaxies in the EAGLE and Illustris-TNG hydrodynamic simulations. A forward modeling pipeline is developed to predict 4-year eROSITA observations and stacked image analysis, including the effects of instrumental and astrophysical backgrounds. We propose an experiment to stack z~0.01 galaxies separated by specific star-formation rate (sSFR) to examine how the hot (T>=10^6 K) circumgalactic medium (CGM) differs for high- and low-sSFR galaxies. The simulations indicate that the hot CGM of low-mass (M_*~10^{10.5} Msol), high-sSFR (defined as the top one-third ranked by sSFR) central galaxies will be detectable to a galactocentric radius r~30-50 kpc. Both simulations predict lower luminosities at fixed stellar mass for the low-sSFR galaxies (the lower third of sSFR) with Illustris-TNG predicting 3x brighter coronae around high-sSFR galaxies than EAGLE. Both simulations predict detectable emission out to r~150-200 kpc for stacks centered on high-mass (M_*~10^{11.0} Msol) galaxies, with EAGLE predicting brighter X-ray halos. The extended soft X-ray luminosity correlates strongly and positively with the mass of circumgalactic gas within the virial radius (f_{CGM}). Prior analyses of both simulations have established that f_{CGM} is reduced by expulsive feedback driven mainly by black hole growth, which quenches galaxy growth by inhibiting replenishment of the ISM. Both simulations predict that eROSITA stacks should not only conclusively detect and resolve the hot CGM around L^* galaxies for the first time, but provide a powerful probe of how the baryon cycle operates, for which there remains an absence of consensus between state-of-the-art simulations

Imprint of Drivers of Galaxy Formation in the Circumgalactic Medium

The majority of baryons reside beyond the optical extent of a galaxy in the circumgalactic and intergalactic media (CGM/IGM). Gaseous halos are inextricably linked to the appearance of their host galaxies through a complex story of accretion, feedback, and continual recycling. The energetic processes, which define the state of gas in the CGM, are the same ones that 1) regulate stellar growth so that it is not over-efficient, and 2) create the diversity of today's galaxy colors, SFRs, and morphologies spanning Hubble's Tuning Fork Diagram. They work in concert to set the speed of growth on the star-forming Main Sequence, transform a galaxy across the Green Valley, and maintain a galaxy's quenched appearance on the Red Sequence. Most baryons in halos more massive than 10^12 Msolar along with their high-energy physics and dynamics remain invisible because that gas is heated above the UV ionization states. We argue that information on many of the essential drivers of galaxy evolution is primarily contained in this "missing" hot gas phase. Completing the picture of galaxy formation requires uncovering the physical mechanisms behind stellar and SMBH feedback driving mass, metals, and energy into the CGM. By opening galactic hot halos to new wavebands, we not only obtain fossil imprints of >13 Gyrs of evolution, but observe on-going hot-mode accretion, the deposition of superwind outflows into the CGM, and the re-arrangement of baryons by SMBH feedback. A description of the flows of mass, metals, and energy will only be complete by observing the thermodynamic states, chemical compositions, structure, and dynamics of T>=10^6 K halos. These measurements are uniquely possible with a next-generation X-ray observatory if it provides the sensitivity to detect faint CGM emission, spectroscopic power to measure absorption lines and gas motions, and high spatial resolution to resolve structures