Search

We present an analysis of \textit{HST}/COS/G160M observations of CIV in the inner circumgalactic medium (CGM) of a novel sample of eight z$\sim$0, L$\approx$L$^{\star}$ galaxies, paired with UV-bright QSOs at impact parameters ($R_\mathrm{proj}$) between 25-130 kpc. The galaxies in this stellar-mass-controlled sample (log$_{10}$M$_{\star}$/M$_{\odot}$ $\sim$ 10.2-10.9 M$_{\odot}$) host super-massive black holes (SMBHs) with dynamically-measured masses spanning log$_{10}$M$_\mathrm{BH}$/M$_{\odot}$ $\sim$ 6.8-8.4; this allows us to compare our results with models of galaxy formation where the integrated feedback history from the SMBH alters the CGM over long timescales. We find that the \ion{C}{IV} column density measurements (N$_{\rm C IV}$) (average log$_{10}$N$_{\rm C IV, CH}$ = 13.94$\pm$0.09 cm$^{-2}$) are largely consistent with existing measurements from other surveys of N$_{\rm C IV}$ in the CGM (average log$_{10}$N$_{\rm C IV, Lit}$ = 13.90$\pm$0.08 cm$^{-2}$), but do not show obvious variation as a function of the SMBH mass. In contrast, specific star-formation rate (sSFR) is highly correlated with the ionized content of the CGM. We find a large spread in sSFR for galaxies with log$_{10}$M$_\mathrm{BH}$/M$_{\odot}$ $>$ 7.0, where the CGM \ion{C}{IV} content shows clear dependence on galaxy sSFR but not M$_\mathrm{BH}$. Our results do not indicate an obvious causal link between CGM CIV and the mass of the galaxy's SMBH; however through comparisons to the EAGLE, Romulus25, $\&$ IllustrisTNG simulations, we find that our sample is likely too small to constrain such causality., Comment: Accepted to ApJ, 28 pages, 7 figures

We analyze the \ovi\ content and kinematics for 126 {\hi}-selected absorbers at $0.14 \lesssim z \lesssim 0.73$ for which the metallicities of their cool photoionized phase have been determined. We separate the absorbers into 100 strong {\lya} forest systems (SLFSs with $15 \la$\,{\colden}\,$< 16.2$) and 26 partial Lyman Limit systems (pLLSs with $16.2\le$\,{\colden}\,$\le 17.2$). The sample is drawn from the COS CGM Compendium (CCC) and has \ovi\ coverage in $\sn \geq 8$ {\it HST}/COS G130M/G160M QSO spectra, yielding a $2\sigma$ completeness level of {\coldenovi}$\,\geq 13.6$. The \ovi\ detection rates differ substantially between low-metallicity (LM; {\xh} $\leq -1.4$) and high-metallicity (HM; {\xh} $> -1.4$) SLFSs, with 20\% and 60\% detection rates, respectively. The \ovi\ detection frequency for the HM and LM pLLSs is, however, similar at $\sim$60\%. The SLFSs and pLLSs without detected \ovi\ are consistent with the absorbing gas being in a single phase, while those with \ovi\ trace multiphase gas. We show that the \ovi\ velocity widths and column densities have different distributions in LM and HM gas. We find a strong correlation between \ovi\ column density and metallicity. The strongest (\coldenovi$\,\ga 14$) and broadest {\ovi} absorbers are nearly always associated with HM absorbers, while weaker \ovi\ absorbers are found in both LM and HM absorbers. From comparisons with galaxy-selected and blind \ovi\ surveys, we conclude absorbers with \coldenovi$\,\ga 14$ most likely arise in the circumgalactic medium (CGM) of star-forming galaxies. Absorbers with weak \ovi\ likely trace the extended CGM or intergalactic medium (IGM), while those without \ovi\ likely originate in the IGM., Comment: Accepted for publication in ApJ

The circum-galactic medium (CGM) can feasibly be mapped by multiwavelength surveys covering broad swaths of the sky. With multiple large datasets becoming available in the near future, we develop a likelihood-free Deep Learning technique using convolutional neural networks (CNNs) to infer broad-scale physical properties of a galaxy's CGM and its halo mass for the first time. Using CAMELS (Cosmology and Astrophysics with MachinE Learning Simulations) data, including IllustrisTNG, SIMBA, and Astrid models, we train CNNs on Soft X-ray and 21-cm (HI) radio 2D maps to trace hot and cool gas, respectively, around galaxies, groups, and clusters. Our CNNs offer the unique ability to train and test on ''multifield'' datasets comprised of both HI and X-ray maps, providing complementary information about physical CGM properties and improved inferences. Applying eRASS:4 survey limits shows that X-ray is not powerful enough to infer individual halos with masses $\log(M_{\rm{halo}}/M_{\odot}) < 12.5$. The multifield improves the inference for all halo masses. Generally, the CNN trained and tested on Astrid (SIMBA) can most (least) accurately infer CGM properties. Cross-simulation analysis -- training on one galaxy formation model and testing on another -- highlights the challenges of developing CNNs trained on a single model to marginalize over astrophysical uncertainties and perform robust inferences on real data. The next crucial step in improving the resulting inferences on physical CGM properties hinges on our ability to interpret these deep-learning models.

We derive predictions from state-of-the-art cosmological galaxy simulations for the spatial distribution of the hot circumgalactic medium (CGM, ${\rm [0.1-1]R_{200c}}$) through its emission lines in the X-ray soft band ($[0.3-1.3]$ keV). In particular, we compare IllustrisTNG, EAGLE, and SIMBA and focus on galaxies with stellar mass $10^{10-11.6}\, \MSUN$ at $z=0$. The three simulation models return significantly different surface brightness radial profiles of prominent emission lines from ionized metals such as OVII(f), OVIII, and FeXVII as a function of galaxy mass. Likewise, the three simulations predict varying azimuthal distributions of line emission with respect to the galactic stellar planes, with IllustrisTNG predicting the strongest angular modulation of CGM physical properties at radial range ${\gtrsim0.3-0.5\,R_{200c}}$. This anisotropic signal is more prominent for higher-energy lines, where it can manifest as X-ray eROSITA-like bubbles. Despite different models of stellar and supermassive black hole (SMBH) feedback, the three simulations consistently predict a dichotomy between star-forming and quiescent galaxies at the Milky-Way and Andromeda mass range, where the former are X-ray brighter than the latter. This is a signature of SMBH-driven outflows, which are responsible for quenching star formation. Finally, we explore the prospect of testing these predictions with a microcalorimeter-based X-ray mission concept with a large field-of-view. Such a mission would probe the extended hot CGM via soft X-ray line emission, determine the physical properties of the CGM, including temperature, from the measurement of line ratios, and provide critical constraints on the efficiency and impact of SMBH feedback on the CGM., Comment: 22 pages, 16 figures. Accepted for publication in MNRAS

The Astro2020 Decadal Survey has identified the mapping of the circumgalactic medium (CGM, gaseous plasma around galaxies) as a key objective. We explore the prospects for characterizing the CGM in and around nearby galaxy halos with a future, large grasp X-ray microcalorimeter. We create realistic mock observations from hydrodynamical simulations (EAGLE, IllustrisTNG, and Simba) that demonstrate a wide range of potential measurements, which will address the open questions in galaxy formation and evolution. By including all background and foreground components in our mock observations, we show why it is impossible to perform these measurements with current instruments, such as X-ray CCDs, and only microcalorimeters will allow us to distinguish the faint CGM emission from the bright Milky Way (MW) foreground emission lines. We find that individual halos of MW mass can, on average and depending on star formation rate, be traced out to large radii, around R500, and for larger galaxies even out to R200, using prominent emission lines, such as OVII, or OVIII. Furthermore, we show that emission line ratios for individual halos can reveal the radial temperature structure. Substructure measurements show that it will be possible to relate azimuthal variations to the feedback mode of the galaxy. We demonstrate the ability to construct temperature, velocity, and abundance ratio maps from spectral fitting for individual galaxy halos, which reveal rotation features, AGN outbursts, and enrichment., Comment: 41 pages, 19 figures, accepted for publication in ApJ

Dwarf galaxies are found to have lost most of their metals via feedback processes; however, there still lacks consistent assessment on the retention rate of metals in their circumgalactic medium (CGM). Here we investigate the metal content in the CGM of 45 isolated dwarf galaxies with $M_*=10^{6.5-9.5}~M_\odot$ ($M_{\rm 200m}=10^{10.0-11.5}~M_\odot$) using {\it HST}/COS. While H I (Ly$\alpha$) is ubiquitously detected ($89\%$) within the CGM, we find low detection rates ($\approx5\%-22\%$) in C II, C IV, Si II, Si III, and Si IV, largely consistent with literature values. Assuming these ions form in the cool ($T\approx10^4$ K) CGM with photoionization equilibrium, the observed H I and metal column density profiles can be best explained by an empirical model with low gas density and high volume filling factor. For a typical galaxy with $M_{\rm 200m}=10^{10.9}~M_\odot$ (median of the sample), our model predicts a cool gas mass of $M_{\rm CGM,cool}\sim10^{8.4}~M_\odot$, corresponding to $\sim2\%$ of the galaxy's baryonic budget. Assuming a metallicity of $0.3Z_\odot$, we estimate that the dwarf galaxy's cool CGM likely harbors $\sim10\%$ of the metals ever produced, with the rest either in more ionized states in the CGM or transported to the intergalactic medium. We further examine the EAGLE simulation and show that H I and low ions may arise from a dense cool medium, while C IV arises from a diffuse warmer medium. Our work provides the community with a uniform dataset on dwarf galaxies' CGM that combines our recent observations, additional archival data and literature compilation, which can be used to test various theoretical models of dwarf galaxies., Comment: Finalized version. Accepted for publication in ApJ

Known as the "Missing Baryon Problem", about one-third of baryons in the local universe remain unaccounted for. The missing baryons are thought to reside in the warm-hot intergalactic medium (WHIM) of the cosmic web filaments, which are challenging to detect. Recent Chandra X-ray observations used a novel stacking analysis and detected an OVII absorption line toward the sightline of a luminous quasar, hinting that the missing baryons may reside in the WHIM. To explore how the properties of the OVII absorption line depend on feedback physics, we compare the observational results with predictions obtained from the Cosmology and Astrophysics with MachinE Learning (CAMEL) Simulation suite. CAMELS consists of cosmological simulations with state-of-the-art supernova (SN) and active galactic nuclei (AGN) feedback models from the IllustrisTNG and SIMBA simulations, with varying strengths. We find that the simulated OVII column densities are higher in the outskirts of galaxies than in the large-scale WHIM, but they are consistently lower than those obtained in the Chandra observations, for all feedback runs. We establish that the OVII distribution is primarily sensitive to changes in the SN feedback prescription, whereas changes in the AGN feedback prescription have minimal impact. We also find significant differences in the OVII column densities between the IllustrisTNG and SIMBA runs. We conclude that the tension between the observed and simulated OVII column densities cannot be explained by the wide range of feedback models implemented in CAMELS., Comment: 11 pages, 9 figures. Matching to the published version on MNRAS

The Line Emission Mapper (LEM) is an X-ray Probe for the 2030s that will answer the outstanding questions of the Universe's structure formation. It will also provide transformative new observing capabilities for every area of astrophysics, and to heliophysics and planetary physics as well. LEM's main goal is a comprehensive look at the physics of galaxy formation, including stellar and black-hole feedback and flows of baryonic matter into and out of galaxies. These processes are best studied in X-rays, and emission-line mapping is the pressing need in this area. LEM will use a large microcalorimeter array/IFU, covering a 30x30' field with 10" angular resolution, to map the soft X-ray line emission from objects that constitute galactic ecosystems. These include supernova remnants, star-forming regions, superbubbles, galactic outflows (such as the Fermi/eROSITA bubbles in the Milky Way and their analogs in other galaxies), the Circumgalactic Medium in the Milky Way and other galaxies, and the Intergalactic Medium at the outskirts and beyond the confines of galaxies and clusters. LEM's 1-2 eV spectral resolution in the 0.2-2 keV band will make it possible to disentangle the faintest emission lines in those objects from the bright Milky Way foreground, providing groundbreaking measurements of the physics of these plasmas, from temperatures, densities, chemical composition to gas dynamics. While LEM's main focus is on galaxy formation, it will provide transformative capability for all classes of astrophysical objects, from the Earth's magnetosphere, planets and comets to the interstellar medium and X-ray binaries in nearby galaxies, AGN, and cooling gas in galaxy clusters. In addition to pointed observations, LEM will perform a shallow all-sky survey that will dramatically expand the discovery space., Comment: 18 pages. White paper for a mission concept to be submitted for the 2023 NASA Astrophysics Probes opportunity. v2: All-sky survey figure expanded, references fixed. v3: Added energy resolution measurements for prototype detector array. v4: Author list and reference fixes

The bimodal absorption system imaging campaign (BASIC) aims to characterize the galaxy environments of a sample of 36 HI-selected partial Lyman limit systems (pLLSs) and Lyman limit systems (LLSs) in 23 QSO fields at $z \lesssim 1$. These pLLSs/LLSs provide a unique sample of absorbers with unbiased and well-constrained metallicities, allowing us to explore the origins of metal-rich and low-metallicity circumgalactic medium (CGM) at $z<1$. Here we present Keck/KCWI and VLT/MUSE observations of 11 of these QSO fields (19 pLLSs) that we combine with HST/ACS imaging to identify and characterize the absorber-associated galaxies. We find 23 unique absorber-associated galaxies, with an average of one associated galaxy per absorber. For seven absorbers, all with $<10\%$ solar metallicities, we find no associated galaxies with $\log M_\star \gtrsim 9.0$ within $\rho/R_{vir}$ and $|\Delta v|/v_{esc} \le$ 1.5 with respect to the absorber. We do not find any strong correlations between the metallicities or HI column densities of the gas and most of the galaxy properties, except for the stellar mass of the galaxies: the low-metallicity ([X/H] $\le -1.4$) systems have a probability of $0.39^{+0.16}_{-0.15}$ for having a host galaxy with $\log M_\star \ge 9.0$ within $\rho/R_{vir} \le 1.5$, while the higher metallicity absorbers have a probability of $0.78^{+0.10}_{-0.13}$. This implies metal-enriched pLLSs/LLSs at $z<1$ are typically associated with the CGM of galaxies with $\log M_\star > 9.0$, whereas low-metallicity pLLSs/LLSs are found in more diverse locations, with one population arising in the CGM of galaxies and another more broadly distributed in overdense regions of the universe. Using absorbers not associated with galaxies, we estimate the unweighted geometric mean metallicity of the intergalactic medium to be [X/H] $\lesssim -2.1$ at $z<1$, which is lower than previously estimated., Comment: Key figures: 13 and 18. Modified abstract. Accepted for publication in ApJ. Updated after referee comments, expanded subsection 4.4

Metallicity estimates of circumgalactic gas based on absorption line measurements typically require photoionization modeling to account for unseen ionization states. We explore the impact of uncertainties in the extreme ultraviolet background (EUVB) radiation on such metallicity determinations for the z < 1 circumgalactic medium (CGM). In particular, we study how uncertainties in the power-law slope of the EUV radiation, $\mathrm{\alpha_{EUVB}}$, from active galactic nuclei affect metallicity estimates in a sample of 34 absorbers with HI column densities between 15.25 < log ($\mathrm{N_{HI}}$ / $\mathrm{cm^{-2}}$) < 17.25 and measured metal ion column densities. We demonstrate the sensitivity of metallicity estimates to changes in the EUV power-law slope of active galactic nuclei, $\mathrm{\alpha_{EUVB}}$, at low redshift (z < 1), showing derived absorber metallicities increase on average by approximately 0.3 dex as the EUV slope is hardened from $\mathrm{\alpha_{EUVB}}$ = -2.0 to -1.4. We use Markov Chain Monte Carlo sampling of photoionization models with $\mathrm{\alpha_{EUVB}}$ as a free parameter to derive metallicities for these absorbers. The current sample of absorbers does not provide a robust constraint on the slope, $\mathrm{\alpha_{EUVB}}$, itself; we discuss how future analyses may provide stronger constraints. Marginalizing over the uncertainty in the slope of the background, we find the average uncertainties in the metallicity determinations increase from 0.08 dex to 0.14 dex when switching from a fixed EUVB slope to one that freely varies. Thus, we demonstrate that EUVB uncertainties can be included in ionization models while still allowing for robust metallicity inferences., Comment: 42 figures (34 in appendix), 3 tables. Accepted for publication in The Astronomical Journal

CMB-HD is a proposed millimeter-wave survey over half the sky that would be ultra-deep (0.5 uK-arcmin) and have unprecedented resolution (15 arcseconds at 150 GHz). Such a survey would answer many outstanding questions about the fundamental physics of the Universe. Major advances would be 1.) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10 h Mpc^(-1)), which probes dark matter particle properties. It will also allow 2.) measurements of the thermal and kinetic Sunyaev-Zel'dovich effects on small scales to map the gas density and velocity, another probe of cosmic structure. In addition, CMB-HD would allow us to cross critical thresholds: 3.) ruling out or detecting any new, light (< 0.1 eV) particles that were in thermal equilibrium with known particles in the early Universe, 4.) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe, and 5.) ruling out or detecting inflationary magnetic fields. CMB-HD would also provide world-leading constraints on 6.) axion-like particles, 7.) cosmic birefringence, 8.) the sum of the neutrino masses, and 9.) the dark energy equation of state. The CMB-HD survey would be delivered in 7.5 years of observing 20,000 square degrees of sky, using two new 30-meter-class off-axis crossed Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors., Comment: Contribution to Snowmass 2021. Note some text overlap with CMB-HD Astro2020 APC and RFI (arXiv:1906.10134, arXiv:2002.12714). Science case further broadened and updated

The Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project was developed to combine cosmology with astrophysics through thousands of cosmological hydrodynamic simulations and machine learning. CAMELS contains 4,233 cosmological simulations, 2,049 N-body and 2,184 state-of-the-art hydrodynamic simulations that sample a vast volume in parameter space. In this paper we present the CAMELS public data release, describing the characteristics of the CAMELS simulations and a variety of data products generated from them, including halo, subhalo, galaxy, and void catalogues, power spectra, bispectra, Lyman-$\alpha$ spectra, probability distribution functions, halo radial profiles, and X-rays photon lists. We also release over one thousand catalogues that contain billions of galaxies from CAMELS-SAM: a large collection of N-body simulations that have been combined with the Santa Cruz Semi-Analytic Model. We release all the data, comprising more than 350 terabytes and containing 143,922 snapshots, millions of halos, galaxies and summary statistics. We provide further technical details on how to access, download, read, and process the data at \url{https://camels.readthedocs.io}., Comment: 18 pages, 3 figures. More than 350 Tb of data from thousands of simulations publicly available at https://www.camel-simulations.org

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 {\it eROSITA}-observed galaxies oriented to be edge-on, we make several observationally-testable predictions for galaxies in the stellar mass range $M_\star=10^{10.7-11.2}\;$M$_{\odot}$. 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>10^6$ 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_\star$ bin, this owes to spheroids having higher stellar and halo masses within that $M_\star$ bin, and obscures the fact that both simulated populations have similar total CGM luminosities at the exact same $M_\star$. Discs have brighter emission inside 20 kpc and more steeply declining profiles with radius than spheroids. We predict that the {\it 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., Comment: 12 pages, 11 figures, 1 table. Submitted to MNRAS. Comments welcome

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~10^4 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., Comment: 61 pages, 18 figures, and 2 tables. This review article is part of the special issue "The Physical Properties of the Groups of Galaxies", edited by L. Lovisari and S. Ettori. Published in MDPI - Universe: https://www.mdpi.com/journal/universe/special_issues/PPGG

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., Comment: 16 pages, 15 figures. Submitted to MNRAS. Comments welcome. (revision submitted on May 22, 2020 fixes Fig. 3. right panel that mistakenly displayed metal mass and now correctly displays total gas mass)

We use the EAGLE (Evolution and Assembly of GaLaxies and their Environments) cosmological simulation to study the distribution of baryons, and far-ultraviolet (O VI), extreme-ultraviolet (Ne VIII) and X-ray (O VII, O VIII, Ne IX, and Fe XVII) line absorbers, around galaxies and haloes of mass $\mathrm{M}_{200c}=10^{11}$-$10^{14.5}\,\mathrm{M}_{\odot}$ at redshift 0.1. EAGLE predicts that the circumgalactic medium (CGM) contains more metals than the interstellar medium across halo masses. The ions we study here trace the warm-hot, volume-filling phase of the CGM, but are biased towards temperatures corresponding to the collisional ionization peak for each ion, and towards high metallicities. Gas well within the virial radius is mostly collisionally ionized, but around and beyond this radius, and for O VI, photoionization becomes significant. When presenting observables we work with column densities, but quantify their relation with equivalent widths by analysing virtual spectra. Virial-temperature collisional ionization equilibrium ion fractions are good predictors of column density trends with halo mass, but underestimate the diversity of ions in haloes. Halo gas dominates the highest column density absorption for X-ray lines, but lower density gas contributes to strong UV absorption lines from O VI and Ne VIII. Of the O VII (O VIII) absorbers detectable in an Athena X-IFU blind survey, we find that 41 (56) per cent arise from haloes with $\mathrm{M}_{200c}=10^{12.0}$-$10^{13.5}\,\mathrm{M}_{\odot}$. We predict that the X-IFU will detect O VII (O VIII) in 77 (46) per cent of the sightlines passing $\mathrm{M}_{\star}=10^{10.5}$-$10^{11.0}\,\mathrm{M}_{\odot}$ galaxies within 100 pkpc (59 (82) per cent for $\mathrm{M}_{\star}>10^{11.0}\,\mathrm{M}_{\odot}$). Hence, the X-IFU will probe covering fractions comparable to those detected with the Cosmic Origins Spectrograph for O VI., Comment: 22 pages, 17 figures + 6 pages references and appendices. Accepted for publication in MNRAS

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., Comment: 10 pages, 4 figures, 1 table, Accepted ApJL

CMB-HD is a proposed ultra-deep (0.5 uk-arcmin), high-resolution (15 arcseconds) millimeter-wave survey over half the sky that would answer many outstanding questions in both fundamental physics of the Universe and astrophysics. This survey would be delivered in 7.5 years of observing 20,000 square degrees, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors., Comment: Response to request for information (RFI) by the Panel of Radio, Millimeter, and Submillimeter Observations from the Ground (RMS) of the Astro2020 Decadal Survey regarding the CMB-HD APC (arXiv:1906.10134). Note some text overlap with original APC. Note also detector count and cost have been reduced by 1/3, and observing time increased by 1/3 compared to original APC; science goals expanded

Many phenomenologically successful cosmological galaxy formation simulations employ kinetic winds to model galactic outflows, a crucial ingredient in obtaining predictions that agree with various observations. Yet systematic studies of how variations in kinetic wind scalings might alter observable galaxy properties are rare. Here we employ GADGET-3 simulations to study how the baryon cycle, stellar mass function, and other galaxy and CGM predictions vary as a function of the assumed outflow speed $v_w$ and the scaling of the mass loading factor $\eta$ with velocity dispersion $\sigma$. We design our fiducial model to reproduce the measured wind properties at 25% of the virial radius from the Feedback In Realistic Environments (FIRE) simulations. We find that a strong dependence of $\eta \sim \sigma^5$ in low mass haloes with $\sigma < 106\ \mathrm{km\ s^{-1}}$ is required to match the faint end of the stellar mass functions at $z > 1$. The wind speed also has a major impact, with faster winds significantly reducing wind recycling and heating more halo gas. Both effects result in less stellar mass growth in massive haloes and impact high ionization absorption in halo gas. We cannot simultaneously match the stellar content at $z=2$ and $z=0$ within a single model, suggesting that an additional feedback source such as AGN might be required in massive galaxies at lower redshifts, but the amount needed depends strongly on assumptions regarding the outflow properties. We run a 50 $\mathrm{Mpc/h}$, $2\times576^3$ simulation with our fiducial parameters and show that it matches a range of star-forming galaxy properties at $z\sim0-2$. In closing, the results from simulations of galaxy formation are much more sensitive to small changes in the feedback implementation than to the hydrodynamic technique., Comment: 31 pages, 19 figures, accepted by MNRAS

We examine the connection between the properties of the circumgalactic medium (CGM) and the quenching and morphological evolution of central galaxies in the EAGLE and IllustrisTNG simulations. The simulations yield very different median CGM mass fractions, $f_{\rm CGM}$, as a function of halo mass, $M_{200}$, with low-mass haloes being significantly more gas-rich in IllustrisTNG than in EAGLE. Nonetheless, in both cases scatter in $f_{\rm CGM}$ at fixed $M_{200}$ is strongly correlated with the specific star formation rate and the kinematic morphology of central galaxies. The correlations are strongest for $\sim L^\star$ galaxies, corresponding to the mass scale at which AGN feedback becomes efficient. This feedback elevates the CGM cooling time, preventing gas from accreting onto the galaxy to fuel star formation, and thus establishing a preference for quenched, spheroidal galaxies to be hosted by haloes with low $f_{\rm CGM}$ for their mass. In both simulations, $f_{\rm CGM}$ correlates negatively with the host halo's intrinsic concentration, and hence with its binding energy and formation redshift, primarily because early halo formation fosters the rapid early growth of the central black hole (BH). This leads to a lower $f_{\rm CGM}$ at fixed $M_{200}$ in EAGLE because the BH reaches high accretion rates sooner, whilst in IllustrisTNG it occurs because the central BH reaches the mass threshold at which AGN feedback is assumed to switch from thermal to kinetic injection earlier. Despite these differences, there is consensus from these state-of-the-art simulations that the expulsion of efficiently-cooling gas from the CGM is a crucial step in the quenching and morphological evolution of central galaxies., Comment: 19 pages, 7 figures, Accepted for publication in MNRAS

We use high-resolution hydrodynamical simulations run with the EAGLE model of galaxy formation to study the differences between the properties of - and subsequently the lensing signal from - subhaloes of massive elliptical galaxies at redshift 0.2, in Cold and Sterile Neutrino (SN) Dark matter models. We focus on the two 7 keV SN models that bracket the range of matter power spectra compatible with resonantly-produced SN as the source of the observed 3.5 keV line. We derive an accurate parametrisation for the subhalo mass function in these two SN models relative to CDM, as well as the subhalo spatial distribution, density profile, and projected number density and the dark matter fraction in subhaloes. We create mock lensing maps from the simulated haloes to study the differences in the lensing signal in the framework of subhalo detection. We find that subhalo convergence is well described by a log-normal distribution and that signal of subhaloes in the power spectrum is lower in SN models with respect to CDM, at a level of 10 to 80 per cent, depending on the scale. However, the scatter between different projections is large and might make the use of power-spectrum studies on the typical scales of current lensing images very difficult. Moreover, in the framework of individual detections through gravitational imaging a sample of ~30 lenses with an average sensitivity of M_sub=5x10^7M_sun would be required to discriminate between CDM and the considered sterile neutrino models., Comment: 16 pages, 11 figures, accepted for publication in MNRAS

A millimeter-wave survey over half the sky, that spans frequencies in the range of 30 to 350 GHz, and that is both an order of magnitude deeper and of higher-resolution than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. By providing such a deep, high-resolution millimeter-wave survey (about 0.5 uK-arcmin noise and 15 arcsecond resolution at 150 GHz), CMB-HD will enable major advances. It will allow 1) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10/hMpc), which probes dark matter particle properties. It will also allow 2) measurements of the thermal and kinetic Sunyaev-Zel'dovich effects on small scales to map the gas density and gas pressure profiles of halos over a wide field, which probes galaxy evolution and cluster astrophysics. In addition, CMB-HD would allow us to cross critical thresholds in fundamental physics: 3) ruling out or detecting any new, light (< 0.1eV), thermal particles, which could potentially be the dark matter, and 4) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe. Such a survey would also 5) monitor the transient sky by mapping the full observing region every few days, which opens a new window on gamma-ray bursts, novae, fast radio bursts, and variable active galactic nuclei. Moreover, CMB-HD would 6) provide a census of planets, dwarf planets, and asteroids in the outer Solar System, and 7) enable the detection of exo-Oort clouds around other solar systems, shedding light on planet formation. CMB-HD will deliver this survey in 5 years of observing half the sky, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. The telescopes will field about 2.4 million detectors (600,000 pixels) in total., Comment: APC White Paper for the Astro2020 Decadal, with updated proposing team

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., Comment: 15 pages, 10 figures, 1 table, submitted to MNRAS

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., Comment: Accepted for publication in MNRAS. Minor changes w.r.t version 1 including updated values for Table 2

A large fraction of the baryons and most of the metals in the Universe are unaccounted for. They likely lie in extended galaxy halos, galaxy groups, and the cosmic web, and measuring their nature is essential to understanding galaxy formation. These environments have virial temperatures >10^5.5 K, so the gas should be visible in X-rays. Here we show the breakthrough capabilities of grating spectrometers to 1) detect these reservoirs of hidden metals and mass, and 2) quantify hot gas flows, turbulence, and rotation around the Milky Way and external galaxies. Grating spectrometers are essential instruments for future X-ray missions, and existing technologies provide 50-1500-fold higher throughput compared to current orbiting instruments., Comment: Submitted to the Astro2020 Decadal (6 pages, 5 figures)

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., Comment: White paper submitted to the Astro2020 Decadal Survey (5 pages, 2 figures)

The past decade has seen an explosion of discoveries and new insights into the diffuse gas within galaxies, galaxy clusters, and the filaments composing the Cosmic Web. A new decade will bring fresh opportunities to further this progress towards developing a comprehensive view of the composition, thermal state, and physical processes of diffuse gas in the Universe. Ultraviolet (UV) spectroscopy, probing diffuse 10^4-10^6 K gas at high spectral resolution, is uniquely poised to (1) witness environmental galaxy quenching processes in action, such as strangulation and tidal- and ram-pressure stripping, (2) directly account for the baryon content of galaxy clusters in the cold-warm (T<10^6 K) gas, (3) determine the phase structure and kinematics of gas participating in the equilibrium-regulating exchange of energy at the cores of galaxy clusters, and (4) map cold streams and filaments of the Cosmic Web that feed galaxies and clusters. With a substantial UV undertaking beyond the Hubble Space Telescope, all of the above would be achievable over the entire epoch of galaxy cluster formation. Such capabilities, coupled with already-planned advancements at other wavelengths, will transform extragalactic astronomy by revealing the dominant formation and growth mechanisms of gaseous halos over the mass spectrum, settling the debate between early- and late-time metal enrichment scenarios, and revealing how the ecosystems in which galaxies reside ultimately facilitate their demise., Comment: White Paper submitted to the Astro2020 Decadal Survey

The past decade began with the first light of ALMA and will end at the start of the new era of hyperdimensional astrophysics. Our community-wide movement toward highly multiwavelength and multidimensional datasets has enabled immense progress in each science frontier identified by the 2010 Decadal Survey, particularly with regard to black hole feedback and the cycle of baryons in galaxies. Facilities like ALMA and the next generation of integral field unit (IFU) spectrographs together enable mapping the physical conditions and kinematics of warm ionized and cold molecular gas in galaxies in unprecedented detail (Fig. 1). JWST's launch at the start of the coming decade will push this capability to the rest-frame UV at redshifts z > 6, mapping the birth of stars in the first galaxies at cosmic dawn. Understanding of their subsequent evolution, however, now awaits an ability to map the processes that transform galaxies directly, rather than the consequences of those processes in isolation. In this paper, we argue that doing so requires an equivalent revolution in spatially resolved spectroscopy for the hot plasma that pervades galaxies, the atmospheres in which they reside, and the winds that are the engines of their evolution., Comment: A Science White Paper submitted to the Astro2020 Decadal Survey (7 papers, 4 figures)

We characterize the metallicities and physical properties of cool, photoionized gas in a sample of 152 z<1 strong Lya forest systems (SLFSs, absorbers with 15

A small survey of the UV-absorbing gas in 12 low-$z$ galaxy groups has been conducted using the Cosmic Origins Spectrograph (COS) on-board the Hubble Space Telescope (HST). Targets were selected from a large, homogeneously-selected sample of groups found in the Sloan Digital Sky Survey (SDSS). A critical selection criterion excluded sight lines that pass close ($<1.5$ virial radii) to a group galaxy, to ensure absorber association with the group as a whole. Deeper galaxy redshift observations are used both to search for closer galaxies and also to characterize these $10^{13.5}$ to $10^{14.5} M_{\odot}$ groups, the most massive of which are highly-virialized with numerous early-type galaxies (ETGs). This sample also includes two spiral-rich groups, not yet fully-virialized. At group-centric impact parameters of 0.3-2 Mpc, these $\mathrm{S/N}=15$-30 spectra detected HI absorption in 7 of 12 groups; high (OVI) and low (SiIII) ion metal lines are present in 2/3 of the absorption components. None of the three most highly-virialized, ETG-dominated groups are detected in absorption. Covering fractions $\gtrsim50$% are seen at all impact parameters probed, but do not require large filling factors despite an enormous extent. Unlike halo clouds in individual galaxies, group absorbers have radial velocities which are too low to escape the group potential well without doubt. This suggests that these groups are "closed boxes" for galactic evolution in the current epoch. Evidence is presented that the cool and warm group absorbers are not a pervasive intra-group medium (IGrM), requiring a hotter ($T\sim10^6$ to $10^7$ K) IGrM to be present to close the baryon accounting., Comment: Resubmitted to ApJS after first review; 82 pages (27 for main text, rest are Appendices and supplemental figures and tables), 47 figures, 21 tables

We present the results from our COS circumgalactic medium (CGM) compendium (CCC), a survey of the CGM at z<1 using HI-selected absorbers with 1519 have mostly -1<[X/H]<0 at z<1. The metal-enriched gas probed by pLLSs and LLSs confirms that galaxies that have been enriching their CGM over billions of years. Surprisingly, despite this enrichment, there is also abundant metal-poor CGM gas (41-59% of the pLLSs have [X/H]<-1.4), representing a reservoir of near-pristine gas around z<1 galaxies. We compare our empirical results to recent cosmological zoom simulations, finding some discrepancies, including an overabundance of metal-enriched CGM gas in simulations., Comment: Resubmitted to ApJ after the first referee report

We test and improve the numerical schemes in our smoothed particle hydrodynamics (SPH) code for cosmological simulations, including the pressure-entropy formulation (PESPH), a time-dependent artificial viscosity, a refined timestep criterion, and metal-line cooling that accounts for photoionisation in the presence of a recently refined Haardt \& Madau (2012) model of the ionising background. The PESPH algorithm effectively removes the artificial surface tension present in the traditional SPH formulation, and in our test simulations it produces better qualitative agreement with mesh-code results for Kelvin-Helmholtz instability and cold cloud disruption. Using a set of cosmological simulations, we examine many of the quantities we have studied in previous work. Results for galaxy stellar and HI mass functions, star formation histories, galaxy scaling relations, and statistics of the Ly$\alpha$ forest are robust to the changes in numerics and microphysics. As in our previous simulations, cold gas accretion dominates the growth of high-redshift galaxies and of low mass galaxies at low redshift, and recycling of winds dominates the growth of massive galaxies at low redshift. However, the PESPH simulation removes spurious cold clumps seen in our earlier simulations, and the accretion rate of hot gas increases by up to an order of magnitude at some redshifts. The new numerical model also influences the distribution of metals among gas phases, leading to considerable differences in the statistics of some metal absorption lines, most notably NeVIII., Comment: 29 pages, 25 figures, accepted by MNRAS

We examine the origin of scatter in the relationship between the gas fraction and mass of dark matter haloes hosting present-day $\sim L^\star$ central galaxies in the EAGLE simulations. The scatter is uncorrelated with the accretion rate of the central galaxy's black hole (BH), but correlates strongly and negatively with the BH's mass, implicating differences in the expulsion of gas by active galactic nucleus feedback, throughout the assembly of the halo, as the main cause of scatter. Haloes whose central galaxies host undermassive BHs also tend to retain a higher gas fraction, and exhibit elevated star formation rates (SFRs). Diversity in the mass of central BHs stems primarily from diversity in the dark matter halo binding energy, as these quantities are strongly and positively correlated at fixed halo mass, such that $\sim L^\star$ galaxies hosted by haloes that are more (less) tightly-bound develop central BHs that are more (less) massive than is typical for their halo mass. Variations in the halo gas fraction at fixed halo mass are reflected in both the soft X-ray luminosity and thermal Sunyaev-Zel'dovich flux, suggesting that the prediction of a strong coupling between the properties of galaxies and their halo gas fractions can be tested with measurements of these diagnostics for galaxies with diverse SFRs but similar halo masses., Comment: 12 pages, 5 figures, accepted for publication in MNRAS

We present a neutral hydrogen-selected absorption-line survey of gas with HI column densities 152 orders of magnitude at any given N(HI), indicating a wide range of metallicities (from solar to <1/100 solar). In the range 16.2

Active galactic nuclei (AGN) are thought to play a critical role in shaping galaxies, but their effect on the circumgalactic medium (CGM) is not well studied. We present results from the COS-AGN survey: 19 quasar sightlines that probe the CGM of 20 optically-selected AGN host galaxies with impact parameters $80 < \rho_{imp} < 300$ kpc. Absorption lines from a variety of species are measured and compared to a stellar mass and impact parameter matched sample of sightlines through non-AGN galaxies. Amongst the observed species in the COS-AGN sample (HI, CII, SiII, SiIII, CIV, SiIV, NV), only Ly$\alpha$ shows a high covering fraction ($94^{+6}_{-23}$% for rest-frame equivalent widths EW $> 124$ m\AA) whilst many of the metal ions are not detected in individual sightlines. A sightline-by-sightline comparison between COS-AGN and the control sample yields no significant difference in EW distribution. However, stacked spectra of the COS-AGN and control samples show significant (> 3 sigma) enhancements in the EW of both Ly$\alpha$ and SiIII at impact parameters $> 164$ kpc by a factor of $+0.45\pm0.05$ dex and $> +0.75$ dex respectively. The lack of detections of both high-ionization species near the AGN and strong kinematic offsets between the absorption systemic galaxy redshifts indicates that neither the AGN's ionization nor its outflows are the origin of these differences. Instead, we suggest the observed differences could result from either AGN hosts residing in haloes with intrinsically distinct gas properties, or that their CGM has been affected by a previous event, such as a starburst, which may also have fuelled the nuclear activity., Comment: Accepted for publication in MNRAS

We present the Technicolor Dawn simulations, a suite of cosmological radiation-hydrodynamic simulations of the first 1.2 billion years. By modeling a spatially-inhomogeneous UVB on-the-fly with 24 frequencies and resolving dark matter halos down to $10^8 M_\odot$ within 12 $h^{-1}$ Mpc volumes, our simulations unify observations of the intergalactic and circumgalactic media, galaxies, and reionization into a common framework. The only empirically-tuned parameter, the fraction $f_{\mathrm{esc,gal}}(z)$ of ionizing photons that escape the interstellar medium, is adjusted to match observations of the Lyman-$\alpha$ forest and the cosmic microwave background. With this single calibration, our simulations reproduce the history of reionization; the stellar mass-star formation rate relation of galaxies; the number density and metallicity of damped Lyman-$\alpha$ absorbers (DLAs) at $z\sim5$; the abundance of weak metal absorbers; the ultraviolet background (UVB) amplitude; and the Lyman-$\alpha$ flux power spectrum at $z=5.4$. The galaxy stellar mass and UV luminosity functions are underproduced by $\leq2\times$, suggesting an overly vigorous feedback model. The mean transmission in the Lyman-$\alpha$ forest is underproduced at $z<6$, indicating tension between measurements of the UVB amplitude and Lyman-$\alpha$ transmission. The observed SiIV column density distribution is reasonably well-reproduced ($\sim 1\sigma$ low). By contrast, CIV remains significantly underproduced despite being boosted by an intense $>4$ Ryd UVB. Solving this problem by increasing metal yields would overproduce both weak absorbers and DLA metallicities. Instead, the observed strength of high-ionization emission from high-redshift galaxies and absorption from their environments suggest that the ionizing flux from conventional stellar population models is too soft., Comment: 24 pages, 17 figures, accepted to MNRAS

We use observations of spatially-aligned C II, C IV, Si II, Si IV, and O I absorbers to probe the slope and intensity of the ultraviolet background (UVB) at $z \sim 6$. We accomplish this by comparing observations with predictions from a cosmological hydrodynamic simulation using three trial UVBs applied in post-processing: a spectrally soft, fluctuating UVB calculated using multi-frequency radiative transfer; a soft, spatially-uniform UVB; and a hard, spatially-uniform "quasars-only" model. When considering our paired high-ionization absorbers (C IV/Si IV), the observed statistics strongly prefer the hard, spatially-uniform UVB. This echoes recent findings that cosmological simulations generically underproduce strong C IV absorbers at $z>5$. A single low/high ionization pair (Si II/Si IV), by contrast, shows a preference for the HM12 UVB, while two more (C II/C IV and O I/C IV) show no preference for any of the three UVBs. Despite this, future observations of specific absorbers, particularly Si IV/C IV, with next-generation telescopes probing to lower column densities should yield tighter constraints on the UVB., Comment: 13 pages, 4 figures, accepted by MNRAS

Hydrostatic equilibrium (HSE), where the thermal pressure gradient balances the force of gravity, is tested across a range of simulated EAGLE haloes from Milky Way L* haloes (M_200~10^12 Msol) to cluster scales. Clusters (M_200>=10^14 Msol) reproduce previous results with thermal pressure responsible for ~90% of the support against gravity, but this fraction drops for group-sized haloes (M_200~10^13 Msol) and is even lower (40-70%) for L* haloes between 0.1-0.3 R_200. Energy from feedback grows relative to the binding energy of a halo toward lower mass resulting in greater deviations from HSE. Tangential motions comprise the largest deviation from HSE in L* haloes indicating that the hot circumgalactic medium (CGM) has significant sub- centrifugal rotation and angular momentum spin parameters 2-3x higher than the dark matter spin parameters. Thermal feedback can buoyantly rise to the outer CGM of M_200<10^12 Msol haloes, both moving baryons beyond R_200 and feeding uncorrelated tangential motions. The resulting hot halo density and rotation profiles show promising agreement with X-ray observations of the inner Milky Way halo, and we discuss future observational prospects to detect spinning hot haloes around other galaxies. Acceleration and radial streaming motions also comprise significant deviations from HSE, especially net outward accelerations seen in L* and group haloes indicating active feedback. Black hole feedback acts in a preventative manner during the later growth of group haloes, applying significant accelerations via shocks that do not feed tangential motions. We argue that HSE is a poor assumption for the CGM, especially in the inner regions, and rotating baryonic hot haloes are a critical consideration for analytic models of the CGM., Comment: Accepted to MNRAS (15 pages, 6 figures) Revisions include more simulation clarifications, clearer nomenclature, and some equation corrections

The metallicity of strong HI systems, spanning from damped Lyman-alpha absorbers (DLAs) to Lyman-limit systems (LLSs) is explored between z = 5->0 using the EAGLE high-resolution cosmological hydrodynamic simulation of galaxy formation. The metallicities of LLSs and DLAs steadily increase with time in agreement with observations. DLAs are more metal rich than LLSs, although the metallicities in the LLS column density range (NHI = 10^17 -10^20 cm^-2) are relatively flat, evolving from a median HI-weighted metallicity of Z<10^-2 Zsol at z = 3 to ~10^-0.5 Zsol by z = 0. The metal content of HI systems tracks the increasing stellar content of the Universe, holding ~5% of the integrated total metals released from stars at z = 0. We also consider partial LLS (pLLS, NHI = 10^16-10^17 cm^-2) metallicities, and find good agreement with Wotta et al. (2016) for the fraction of systems above (37%) and below (63%) 0.1 Zsol. We also find a large dispersion of pLLS metallicities, although we do not reproduce the observed metallicity bimodality and instead we make the prediction that a larger sample will yield more pLLSs around 0.1 Zsol. We under-predict the median metallicity of strong LLSs, and predict a population of Z < 10^-3 Zsol DLAs at z > 3 that are not observed, which may indicate more widespread early enrichment in the real Universe compared to EAGLE., Comment: 8 pages, 6 figures, accepted by MNRAS

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., Comment: 26 pages, 19 figures, submitted to MNRAS

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., Comment: 17 pages, 12 figures, accepted in MNRAS, companion website: http://noneq.strw.leidenuniv.nl/PZF/

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., Comment: 17 pages, 11 figures (not including appendix), accepted for publication in MNRAS

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., Comment: 25 pages, 14 figures, 2 tables. Accepted in MNRAS

The thermal and chemical properties of the hot diffuse intragroup medium (IGrM) provide important constraints on the feedback processes associated with massive galaxy formation and evolution. Here we explore these constraints via a detailed analysis of the global properties of simulated z<3 galaxy groups from a cosmological simulation that includes a well-constrained prescription for stellar/supernovae-powered galactic outflows but no AGN feedback. Our aims are to (a) establish a baseline against which we will compare future models; (b) identify model successes due to stellar/supernovae-powered outflows; and (c) pinpoint features that signal the need for, and constrain the nature of, AGN feedback. Our simulation successfully reproduces key observed z=0 group IGrM properties, including the various X-ray Lx - Tx - entropy scaling relations, for all but the most massive groups. The z<1 redshift evolution of these also agree with observations. Contrary to expectations, the simulated groups' IGrM does not suffer catastrophic cooling. Yet, the z=0 group stellar mass is ~ 2X too large. This is due to the build-up of cold gas in the massive galaxies before they are incorporated inside groups. This in turn indicates that other feedback mechanisms must activate in real galaxies once their stellar masses grow to a few X 10^{10} M_sun. We show that these must be powerful enough to expel a significant fraction of the gas from the galactic halos. Gentle maintenance-mode (quenching) AGN feedback, as seen in galaxy clusters, will not do. Just as importantly, we find that the stellar/supernovae-powered winds are essential for understanding the IGrM metal abundances. Our simulation is able to reproduce the observed relationship between the global IGrM iron and silicon abundance and the group X-ray temperature, and these results ought to be relatively insensitive to the addition of AGN feedback., Comment: 26 pages, 16 figures, accepted by MNRAS

We study the distribution and evolution of highly ionised intergalactic metals in the Evolution and Assembly of Galaxies and their Environment (EAGLE) cosmological, hydrodynamical simulations. EAGLE has been shown to reproduce a wide range of galaxy properties while its subgrid feedback was calibrated without considering gas properties. We compare the predictions for the column density distribution functions (CDDFs) and cosmic densities of SiIV, CIV, NV, OVI and NeVIII absorbers with observations at redshift z = 0 to ~ 6 and find reasonable agreement, although there are some differences. We show that the typical physical densities of the absorbing gas increase with column density and redshift, but decrease with the ionization energy of the absorbing ion. The typical metallicity increases with both column density and time. The fraction of collisionally ionized metal absorbers increases with time and ionization energy. While our results show little sensitivity to the presence or absence of AGN feedback, increasing/decreasing the efficiency of stellar feedback by a factor of two substantially decreases/increases the CDDFs and the cosmic densities of the metal ions. We show that the impact of the efficiency of stellar feedback on the CDDFs and cosmic densities is largely due to its effect on the metal production rate. However, the temperatures of the metal absorbers, particularly those of strong OVI, are directly sensitive to the strength of the feedback., Comment: 25 pages, 20 figures, 2 tables; Accepted for publication in MNRAS after minor changes

Simulations of the formation of large-scale structure predict that dark matter, low density highly ionized gas, and galaxies form 10 40 Mpc scale filaments. These structure are easily recognized in the distribution of galaxies, but have not been directly observed in the distribution of the gas. We use Ly-alpha absorption lines in the spectra of 24 AGN to present a new way to probe these filaments. We use a new catalogue of nearby (cz<10,000 km/s) galaxies, complete down to a luminosity of about 0.05 L* for the region of space analyzed here. Using HST spectra of 24 AGN we sample the gas associated with a 30x5 Mpc galaxy filament at cz~3500 km/s. All of our sightlines pass outside the virial radius of any known filament galaxy. Within 500 kpc of the filament axis the detection rate is ~80%, while no detections are seen more than 2.1 Mpc from the filament. The width of the Lya lines correlates with filament impact parameter and the four BLAs in our sample all occur within 400 kpc of the filament axis, indicating increased temperature and/or turbulence. Comparing to simulations, we find that the recent Haardt & Madau (2012) extragalactic ionizing background predicts a factor 3-5 too few ionizing photons. Using a more intense radiation field matches the hydrogen density profile within 2.1 Mpc of the filament axis, but the simulations still overpredict the detection rate between 2.1 and 5 Mpc from the axis., Comment: 31 pages, 28 figures, replaced original with updated version adding two figures and one subsection

We analyze the low-redshift (z~0.2) circumgalactic medium by comparing absorption-line data from the COS-Halos Survey to absorption around a matched galaxy sample from two cosmological hydrodynamic simulations. The models include different prescriptions for galactic outflows, namely hybrid energy/momentum driven wind (ezw), and constant winds (cw). We extract for comparison direct observables including equivalent widths, covering factors, ion ratios, and kinematics. Both wind models are generally in good agreement with these observations for HI and certain low ionization metal lines, but show poorer agreement with higher ionization metal lines including SiIII and OVI that are well-observed by COS-Halos. These discrepancies suggest that both wind models predict too much cool, metal-enriched gas and not enough hot gas, and/or that the metals are not sufficiently well-mixed. This may reflect our model assumption of ejecting outflows as cool and unmixing gas. Our ezw simulation includes a heuristic prescription to quench massive galaxies by super-heating its ISM gas, which we show yields sufficient low ionisation absorption to be broadly consistent with observations, but also substantial OVI absorption that is inconsistent with data, suggesting that gas around quenched galaxies in the real Universe does not cool. At impact parameters of <50 kpc, recycling winds dominate the absorption of low ions and even HI, while OVI almost always arises from metals ejected longer than 1 Gyr ago. The similarity between the wind models is surprising, since we show that they differ substantially in their predicted amount and phase distribution of halo gas. We show that this similarity owes mainly to our comparison here at fixed stellar mass rather than at fixed halo mass in our previous works, which suggests that CGM properties are more closely tied to the stellar mass of galaxies rather than halo mass., Comment: 21 pages, 11 figures, submitted to MNRAS