Search

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

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, 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 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., 12 pages, 11 figures, 1 table. Submitted to MNRAS. Comments welcome

We report a tentative detection of the circumgalactic medium (CGM) of WLM, an isolated, low-mass (log$M_*/M_\odot\approx7.6$), dwarf irregular galaxy in the Local Group (LG). We analyze an HST/COS archival spectrum of a quasar sightline (PHL2525), which is 45 kpc (0.5 virial radius) from WLM and close to the Magellanic Stream (MS). Along this sightline, two ion absorbers are detected in Si II, Si III, Si IV, C II, and C IV at velocities of $\sim$-220 km s$^{-1}$ (Component v-220) and $\sim$-150 km s$^{-1}$ (Component v-150). To identify their origins, we study the position-velocity alignment of the components with WLM and the nearby MS. Near the Magellanic longitude of PHL2525, the MS-related neutral and ionized gas moves at $\lesssim-190$ km s$^{-1}$, suggesting an MS origin for Component v-220, but not for Component v-150. Because PHL2525 passes near WLM and Component v-150 is close to WLM's systemic velocity ($\sim$-132 km s$^{-1}$), it is likely that Component v-150 arises from the galaxy's CGM. This results in a total Si mass in WLM's CGM of $M_{\rm Si}^{\rm CGM}\sim(0.2-1.0)\times10^5~M_\odot$ using assumption from other COS dwarf studies. Comparing $M_{\rm Si}^{\rm CGM}$ to the total Si mass synthesized in WLM over its lifetime ($\sim$1.3$\times10^5~M_\odot$), we find $\sim$3% is locked in stars, $\sim$6% in the ISM, $\sim$15%-77% in the CGM, and the rest ($\sim$14%-76%) is likely lost beyond the virial radius. Our finding resonates with other COS dwarf galaxy studies and theoretical predictions that low-mass galaxies can easily lose metals into their CGM due to stellar feedback and shallow gravitational potential., Accepted for publication in MNRAS

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

We present a study of the galaxy environment of 9 strong HI+CIV absorption line systems ($16.2, Comment: 19 pages (+2 pages appendices), 19 figures (+4 in appendices). Submitted to MNRAS

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 < ��_{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$��$ shows a high covering fraction ($94^{+6}_{-23}$% for rest-frame equivalent widths EW $> 124$ m��) 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$��$ 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., 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., 16 pages, 11 figures, accepted for publication in 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

Recent analysis of strongly-lensed sources in the Hubble Frontier Fields indicates that the rest-frame UV luminosity function of galaxies at $z=$6--8 rises as a power law down to $M_\mathrm{UV}=-15$, and possibly as faint as -12.5. We use predictions from a cosmological radiation hydrodynamic simulation to map these luminosities onto physical space, constraining the minimum dark matter halo mass and stellar mass that the Frontier Fields probe. While previously-published theoretical studies have suggested or assumed that early star formation was suppressed in halos less massive than $10^9$--$10^{11} M_\odot$, we find that recent observations demand vigorous star formation in halos at least as massive as (3.1, 5.6, 10.5)$\times10^9 M_\odot$ at $z=(6,7,8)$. Likewise, we find that Frontier Fields observations probe down to stellar masses of (8.1, 18, 32)$\times10^6 M_\odot$; that is, they are observing the likely progenitors of analogues to Local Group dwarfs such as Pegasus and M32. Our simulations yield somewhat different constraints than two complementary models that have been invoked in similar analyses, emphasizing the need for further observational constraints on the galaxy-halo connection., Comment: 8 pages, 5 figures, re-submitted to MNRAS after incorporating referee's comments

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, 26 pages, 16 figures, accepted by MNRAS

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 ($, 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 z19 have mostly -1, 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$��$ 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., 29 pages, 25 figures, accepted by 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 $z4$ 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., 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, 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 3 that are not observed, which may indicate more widespread early enrichment in the real Universe compared to EAGLE., 8 pages, 6 figures, accepted by MNRAS

Observations suggest that CII was more abundant than CIV in the intergalactic medium towards the end of the hydrogen reionization epoch. This transition provides a unique opportunity to study the enrichment history of intergalactic gas and the growth of the ionizing background (UVB) at early times. We study how carbon absorption evolves from z=10-5 using a cosmological hydrodynamic simulation that includes a self-consistent multifrequency UVB as well as a well-constrained model for galactic outflows to disperse metals. Our predicted UVB is within 2-4 times that of Haardt & Madau (2012), which is fair agreement given the uncertainties. Nonetheless, we use a calibration in post-processing to account for Lyman-alpha forest measurements while preserving the predicted spectral slope and inhomogeneity. The UVB fluctuates spatially in such a way that it always exceeds the volume average in regions where metals are found. This implies both that a spatially-uniform UVB is a poor approximation and that metal absorption is not sensitive to the epoch when HII regions overlap globally even at column densites of 10^{12} cm^{-2}. We find, consistent with observations, that the CII mass fraction drops to low redshift while CIV rises owing the combined effects of a growing UVB and continued addition of carbon in low-density regions. This is mimicked in absorption statistics, which broadly agree with observations at z=6-3 while predicting that the absorber column density distributions rise steeply to the lowest observable columns. Our model reproduces the large observed scatter in the number of low-ionization absorbers per sightline, implying that the scatter does not indicate a partially-neutral Universe at z=6., 16 pages, 14 figures, accepted 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, 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., 17 pages, 11 figures (not including appendix), accepted for publication in MNRAS

In preparation for a Hubble Space Telescope (HST) observing project using the Cosmic Origins Spectrograph (COS), the positions of all AGN targets having high-S/N far-UV G130M spectra were cross-correlated with a large catalog of low-redshift galaxy groups homogeneously selected from the spectroscopic sample of the Sloan Digital Sky Survey (SDSS). Searching for targets behind only those groups at z = 0.1-0.2 (which places the OVI doublet in the wavelength region of peak COS sensitivity) we identified only one potential S/N = 15-20 target, FBQS 1010+3003. An OVI-only absorber was found in its G130M spectrum at z = 0.11326, close to the redshift of a foreground small group of luminous galaxies at z = 0.11685. Because there is no associated Lyalpha absorption, any characterization of this absorber is necessarily minimal; however, the OVI detection likely traces "warm" gas in collisional ionization equilibrium at T ~ 300,000 K. While this discovery is consistent with being interface gas between cooler, photoionized clouds and a hotter intra-group medium, it could also be warm, interface gas associated with the circum-galactic medium (CGM) of the single closest galaxy. In this case a detailed analysis of the galaxy distribution (complete to 0.2 L*) strongly favors the individual galaxy association. This analysis highlights the necessity of both high-S/N > 20 COS data and a deep galaxy redshift survey of the region in order to test more rigorously the association of OVI-absorbing gas with a galaxy group. A Cycle 23 HST/COS program currently is targeting 10 UV-bright AGN behind 12 low-redshift galaxy groups to test the warm, group gas hypothesis., 11 pages, 3 figures, submitted to ApJ

We use a radiation hydrodynamic simulation of the hydrogen reionization epoch to study OI absorbers at z~6. The intergalactic medium (IGM) is reionized before it is enriched, hence OI absorption originates within dark matter halos. The predicted abundance of OI absorbers is in reasonable agreement with observations. At z=10, roughly 70% of sightlines through atomically-cooled halos encounter a visible (N_OI > 10^14 cm^-2) column. Reionization ionizes and removes gas from halos less massive than 10^8.4 M_0, but 20% of sightlines through more massive halos encounter visible columns even at z=5. The mass scale of absorber host halos is 10-100 times smaller than the halos of Lyman break galaxies and Lyman-alpha emitters, hence absorption probes the dominant ionizing sources more directly. OI absorbers have neutral hydrogen columns of 10^19-10^21 cm^-2, suggesting a close resemblance between objects selected in OI and HI absorption. Finally, the absorption in the foreground of the z=7.085 quasar ULASJ1120+0641 cannot originate in a dark matter halo because halo gas at the observed HI column density is enriched enough to violate the upper limits on the OI column. By contrast, gas at less than one third the cosmic mean density satisfies the constraints. Hence the foreground absorption likely originates in the IGM., 21 pages, 15 figures, re-submitted to MNRAS after many improvements suggested by referee, results unchanged

Simulations of the formation of galaxies, as well as ionisation models used to interpret observations of quasar absorption lines, generally either assume ionisation equilibrium or ignore the presence of the extra-galactic background (EGB) radiation. We introduce a method to compute the non-equilibrium ionisation and cooling of diffuse gas exposed to the EGB. Our method iterates the ionisation states of the 11 elements that dominate the cooling and uses tabulated ion-by-ion cooling and photo-heating efficiencies to update the temperature of the gas. Our reaction network includes radiative and di-electric recombination, collisional ionisation, photo-ionisation, Auger ionisation, and charge transfer. We verify that our method reproduces published results for collisional equilibrium, collisional non-equilibrium, and photo-ionised equilibrium. Non-equilibrium effects can become very important in cooling gas, particularly below 10^6 K. Photo-ionisation and non-equilibrium effects both tend to boost the degree of ionisation and to reduce cooling efficiencies. The effect of the EGB is larger for lower densities. Hence, photo-ionisation affects cooling more under isochoric than under isobaric conditions. Non-equilibrium effects are smaller in the presence of the EGB and are thus overestimated when using collisional-only processes. The inclusion of the EGB alters the observational diagnostics of diffuse, metal-enriched gas (e.g. metal absorption lines) even more significantly than the cooling efficiencies. We argue that the cooling efficiency should be considered if ionisation models are used to infer physical conditions from observed line ratios, as the a priori probability of observing gas is lower if its cooling time is shorter. We provide on-line tables of ionisation fractions and cooling efficiencies for equilibrium and non-equilibrium scenarios. (abridged), Acceptied to MNRAS, 16 figures, paper website- http://noneq.strw.leidenuniv.nl/

The total contribution of diffuse halo gas to the galaxy baryon budget strongly depends on its dominant ionization state. In this paper, we address the physical conditions in the highly ionized circumgalactic medium (CGM) traced by OVI absorption lines observed in COS-Halos spectra. We analyze the observed ionic column densities, absorption-line widths and relative velocities, along with the ratios of NV/OVI for 39 fitted Voigt profile components of OVI. We compare these quantities with the predictions given by a wide range of ionization models. Photoionization models that include only extragalactic UV background radiation are ruled out; conservatively, the upper limits to NV/OVI and measurements of N$_{\rm OVI}$ imply unphysically large path lengths > 100 kpc. Furthermore, very broad OVI absorption (b > 40 km s$^{-1}$) is a defining characteristic of the CGM of star-forming L* galaxies. We highlight two possible origins for the bulk of the observed OVI: (1) highly structured gas clouds photoionized primarily by local high energy sources or (2) gas radiatively cooling on large scales behind a supersonic wind. Approximately 20% of circumgalactic OVI does not align with any low-ionization state gas within $\pm$50 km s$^{-1}$ and is found only in halos with M$_{\rm halo}$ < 10$^{12}$ M$_{\odot}$. We suggest that this type of unmatched OVI absorption traces the hot corona itself at a characteristic temperature of 10$^{5.5}$ K. We discuss the implications of these very distinct physical origins for the dynamical state, gas cooling rates, and total baryonic content of L* gaseous halos., 25 pages, 13 figures. Accepted for publication in ApJ October 6, 2016

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., 25 pages, 20 figures, 2 tables; Accepted for publication in MNRAS after minor changes

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

Recent studies have found a dramatic difference between the observed number density evolution of low mass galaxies and that predicted by semi-analytic models. While 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 observationally. 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. Firstly, we look for evidence of a different history of today's low mass galaxies in models and observations and we find that 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 (sSFR) 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 sSFR 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. [Abridged]

Subhalo abundance matching (SHAM) is a technique for populating simulated dark matter distributions with galaxies, assuming a monotonic relation between a galaxy's stellar mass or luminosity and the mass of its parent dark matter halo or subhalo. We examine the accuracy of SHAM in two cosmological SPH simulations, one of which includes momentum-driven winds. The SPH simulations indeed show a nearly monotonic relation between stellar mass and halo mass provided that, for satellite galaxies, we use the mass of the subhalo at the epoch when it became a satellite. In each simulation, the median relation for central and satellite galaxies is nearly identical, though a somewhat larger fraction of satellites are outliers. SHAM-assigned masses (at z=0-2), luminosities (R-band at z=0), or star formation rates (at z=2) have a 68% scatter of 0.09-0.15 dex relative to the true simulation values. When we apply SHAM to the subhalo population of collisionless N-body simulation with the same initial conditions as the SPH runs, we find generally good agreement for the halo occupation distributions and halo radial profiles of galaxy samples defined by thresholds in stellar mass. However, because a small fraction of SPH galaxies suffer severe stellar mass loss after becoming satellites, SHAM slightly overpopulates high mass halos; this effect is more significant for the wind simulation, which produces galaxies that are less massive and more fragile. SHAM recovers the two-point correlation function of the SPH galaxies in the no-wind simulation to better than 10% at scales 0.1 2 Mpc/h, but overpopulation of massive halos increases the correlation function by a factor of ~2.5 on small scales.