Source author record

Benjamin D. Oppenheimer

Benjamin D. Oppenheimer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA hep-ph astro-ph.HE astro-ph astro-ph.IM hep-ex

Catalog footprint

What is connected

55works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2023arXiv

X-ray Absorption Lines in the Warm-Hot Intergalactic Medium: Probing Chandra observations with the CAMEL simulations

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.

preprint2022arXiv

Snowmass2021 CMB-HD White Paper

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.

preprint2022arXiv

The COS CGM Compendium. IV. Effects of Varying Ionization Backgrounds on Metallicity Determinations in the z < 1 Circumgalactic Medium

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{α_{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{α_{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{α_{EUVB}}$ = -2.0 to -1.4. We use Markov Chain Monte Carlo sampling of photoionization models with $\mathrm{α_{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{α_{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.

preprint2021arXiv

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

We investigate how the X-ray circumgalactic medium (CGM) of present-day galaxies depends on galaxy morphology and azimuthal angle using mock observations generated from the EAGLE cosmological hydrodynamic simulation. By creating mock stacks of {\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.

preprint2020arXiv

CMB-HD: Astro2020 RFI Response

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.

preprint2020arXiv

EAGLE and Illustris-TNG predictions for resolved eROSITA X-ray observations of the circumgalactic medium around normal galaxies

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

preprint2020arXiv

The Impact of Wind Scalings on Stellar Growth and the Baryon Cycle in Cosmological Simulations

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 $η$ with velocity dispersion $σ$. 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 $η\sim σ^5$ in low mass haloes with $σ< 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.

preprint2020arXiv

The Robustness of Cosmological Hydrodynamic Simulation Predictions to Changes in Numerics and Cooling Physics

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.

preprint2020arXiv

The warm-hot circumgalactic medium around EAGLE-simulation galaxies and its detection prospects with X-ray and UV line absorption

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.

preprint2019arXiv

CMB-HD: An Ultra-Deep, High-Resolution Millimeter-Wave Survey Over Half the Sky

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.

preprint2019arXiv

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

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

preprint2019arXiv

The COS CGM Compendium. III: Metallicity and Physical Properties of the Cool Circumgalactic Medium at z<1

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<log N(HI)<16.2). The sample is drawn from our COS circumgalactic medium (CGM) compendium (CCC), an ultraviolet survey of HI-selected circumgalactic gas around z<1 galaxies that targets 262 absorbers with 15<log N(HI)<19. We show that the metallicity probability distribution function of the SLFSs at z<1 is unimodal, skewed to low metallicities with a mean and median of [X/H]=-1.47$ and -1.18 dex. Very metal-poor gas with [X/H]<-1.4 represents about half of the population of absorbers with 15<log N(HI)<18. Thus, there are important reservoirs of primitive (though not pristine) gas around z<1 galaxies. The photoionized gas around z<1 galaxies is highly inhomogeneous based on the wide range of metallicities observed (-3<[X/H]<+0.4) and that there are large metallicity variations (factors of 2 to 25) for most of the closely-spaced absorbers (Dv<300 km/s) along the same sightlines. These absorbers show a complex evolution with redshift and HI column density, and we identify subtle cosmic evolution effects that affect the interpretation of metallicity distributions and comparison with other of absorbers samples. We discuss the physical conditions and cosmic baryon and metal budgets of the CCC absorbers. Finally, we compare the CCC results to recent cosmological zoom simulations and explore the origins of the 15<log N(HI)<19 absorbers within the EAGLE high-resolution simulations.

preprint2019arXiv

The lensing properties of subhaloes in massive elliptical galaxies in sterile neutrino cosmologies

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.

preprint2016arXiv

Bimodality of low-redshift circumgalactic O VI in non-equilibrium EAGLE zoom simulations

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.

preprint2016arXiv

Cosmic distribution of highly ionized metals and their physical conditions in the EAGLE simulations

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.

preprint2015arXiv

Baryon Cycling in the Low-Redshift Circumgalactic Medium: A Comparison of Simulations to the COS-Halos Survey

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.

preprint2015arXiv

Nearby galaxy filaments and the Lya forest: confronting simulations and the UV background with observations

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.

preprint2015arXiv

The Growth and Enrichment of the Intragroup Gas

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.

preprint2015arXiv

The impact of environment and mergers on the HI content of galaxies in hydrodynamic simulations

We quantitatively examine the effects of accretion and environment on the HI content of galaxies within a cosmological hydrodynamic simulation that reproduces basic observed trends of HI in galaxies. We show that our model broadly reproduces the observed scatter in HI at a given stellar mass as quantified by the HI mass function in bins of stellar mass, as well as the HI richness versus local galaxy density. This shows that the predicted HI fluctuations and environmental effects are roughly consistent with data with few minor discrepancies. For satellite galaxies in >= 10^12M_* halos, the HI richness distribution is bimodal and drops towards the largest halo masses. The depletion rate of HI once a galaxy enters a more massive halo is more rapid at higher halo mass, in contrast to the specific star formation rate which shows much less variation in the attenuation rate versus halo mass. This suggests that, up to halo mass scales probed here (<= 10^14M_*), star formation is mainly attenuated by starvation, but HI is additionally removed by stripping once a hot gaseous halo is present. In low mass halos, the HI richness of satellites is independent of radius, while in high mass halos they become gas-poor towards the center. We show that the gas fraction of satellite and central galaxies decreases from z=5 to z=0,tracking each other until z~1 after which the satellites' HI content drops much more quickly, particularly for the highest halo masses. Mergers somewhat increase the HI richness and its scatter about the mean relation, tracking the metallicity in a way consistent with arising from inflow fluctuations, while star formation is significantly boosted relative to HI.

preprint2015arXiv

Torque-Limited Growth of Massive Black Holes in Galaxies Across Cosmic Time

We combine cosmological hydrodynamic simulations with analytic models to evaluate the role of galaxy-scale gravitational torques on the evolution of massive black holes at the centers of star-forming galaxies. We confirm and extend our earlier results to show that torque-limited growth yields black holes and host galaxies evolving on average along the Mbh-Mbulge relation from early times down to z = 0 and that convergence onto the scaling relation occurs independent of the initial conditions and with no need for mass averaging through mergers or additional self-regulation processes. Smooth accretion dominates the long-term evolution, with black hole mergers with mass ratios >1:5 representing typically a small fraction of the total growth. Winds from the accretion disk are required to eject significant mass to suppress black hole growth, but there is no need for coupling this wind to galactic-scale gas to regulate black holes in a non-linear feedback loop. Torque-limited growth yields a close-to-linear relation for the star formation rate and the black hole accretion rate averaged over galaxy evolution time scales. However, the SFR-AGN connection has significant scatter owing to strong variability of black hole accretion at all resolved time scales. Eddington ratios can be described by a broad lognormal distribution with median value evolving roughly as (1 + z)^1.9, suggesting a main sequence for black hole growth similar to the cosmic evolution of specific SFRs. Our results offer an attractive scenario consistent with available observations in which cosmological gas infall and transport of angular momentum in the galaxy by gravitational instabilities regulate the long-term co-evolution of black holes and star-forming galaxies.

preprint2014arXiv

A Budget and Accounting of Metals at z~0: Results from the COS-Halos Survey

We present a budget and accounting of metals in and around star-forming galaxies at $z\sim 0$. We combine empirically derived star formation histories with updated supernova and AGB yields and rates to estimate the total mass of metals produced by galaxies with present-day stellar mass of $10^{9.3}$--$10^{11.6} M_{\odot}$. On the accounting side of the ledger, we show that a surprisingly constant 20--25% mass fraction of produced metals remain in galaxies' stars, interstellar gas and interstellar dust, with little dependence of this fraction on the galaxy stellar mass (omitting those metals immediately locked up in remnants). Thus, the bulk of metals are outside of galaxies, produced in the progenitors of today's $L^*$ galaxies. The COS-Halos survey is uniquely able to measure the mass of metals in the circumgalactic medium (to impact parameters of $< 150$ kpc) of low-redshift $\sim L^*$ galaxies. Using these data, we map the distribution of CGM metals as traced by both the highly ionized OVI ion and a suite of low-ionization species; combined with constraints on circumgalactic dust and hotter X-ray emitting gas out to similar impact parameters, we show that $\sim 40$% of metals produced by $M_{\star}\sim 10^{10}M_{\odot}$ galaxies can be easily accounted for out to 150 kpc. With the current data, we cannot rule out a constant mass of metals within this fixed physical radius. This census provides a crucial boundary condition for the eventual fate of metals in galaxy evolution models.

preprint2014arXiv

Cosmological Zoom Simulations of z = 2 Galaxies: The Impact of Galactic Outflows

We use high-resolution cosmological zoom simulations with ~200 pc resolution at z = 2 and various prescriptions for galactic outflows in order to explore the impact of winds on the morphological, dynamical, and structural properties of eight individual galaxies with halo masses ~ 10^11--2x10^12 Msun at z = 2. We present a detailed comparison to spatially and spectrally resolved Hα and other observations of z ~ 2 galaxies. We find that simulations without winds produce massive, compact galaxies with low gas fractions, super-solar metallicities, high bulge fractions, and much of the star formation concentrated within the inner kpc. Strong winds are required to maintain high gas fractions, redistribute star-forming gas over larger scales, and increase the velocity dispersion of simulated galaxies, more in agreement with the large, extended, turbulent disks typical of high-redshift star-forming galaxies. Winds also suppress early star formation to produce high-redshift cosmic star formation efficiencies in better agreement with observations. Sizes, rotation velocities, and velocity dispersions all scale with stellar mass in accord with observations. Our simulations produce a diversity of morphological characteristics - among our three most massive galaxies, we find a quiescent grand-design spiral, a very compact star-forming galaxy, and a clumpy disk undergoing a minor merger; the clumps are evident in Hα but not in the stars. Rotation curves are generally slowly rising, particularly when calculated using azimuthal velocities rather than enclosed mass. Our results are broadly resolution-converged. These results show that cosmological simulations including outflows can produce disk galaxies similar to those observed during the peak epoch of cosmic galaxy growth.

preprint2014arXiv

Non-equilibrium chemistry and cooling in the diffuse interstellar medium - I. Optically thin regime

An accurate treatment of the multiphase interstellar medium (ISM) in hydrodynamic galaxy simulations requires that we follow not only the thermal evolution of the gas, but also the evolution of its chemical state, including its molecular chemistry, without assuming chemical (including ionisation) equilibrium. We present a reaction network that can be used to solve for this thermo-chemical evolution. Our model follows the evolution of all ionisation states of the 11 elements that dominate the cooling rate, along with important molecules such as H2 and CO, and the intermediate molecular species that are involved in their formation (20 molecules in total). We include chemical reactions on dust grains, thermal processes involving dust, cosmic ray ionisation and heating and photochemical reactions. We focus on conditions typical for the diffuse ISM, with densities of 10^-2 cm^-3 < nH < 10^4 cm^-3 and temperatures of 10^2 K < T < 10^4 K, and we consider a range of radiation fields, including no UV radiation. In this paper we consider only gas that is optically thin, while paper II considers gas that becomes shielded from the radiation field. We verify the accuracy of our model by comparing chemical abundances and cooling functions in chemical equilibrium with the photoionisation code Cloudy. We identify the major coolants in diffuse interstellar gas to be CII, SiII and FeII, along with OI and H2 at densities nH > 10^2 cm^-3. Finally, we investigate the impact of non-equilibrium chemistry on the cooling functions of isochorically or isobarically cooling gas. We find that, at T < 10^4 K, recombination lags increase the electron abundance above its equilibrium value at a given temperature, which can enhance the cooling rate by up to two orders of magnitude. The cooling gas also shows lower H2 abundances than in equilibrium, by up to an order of magnitude.

preprint2014arXiv

Non-equilibrium chemistry and cooling in the diffuse interstellar medium - II. Shielded gas

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

preprint2014arXiv

Parametrising Star Formation Histories

We examine the star formation histories (SFHs) of galaxies in smoothed particle hydrodynamics (SPH) simulations, compare them to parametric models that are commonly used in fitting observed galaxy spectral energy distributions, and examine the efficacy of these parametric models as practical tools for recovering the physical parameters of galaxies. The commonly used tau-model, with SFR ~ exp(-t/tau), provides a poor match to the SFH of our SPH galaxies, with a mismatch between early and late star formation that leads to systematic errors in predicting colours and stellar mass-to-light ratios. A one-parameter lin-exp model, with SFR ~ t*exp(-t/tau), is much more successful on average, but it fails to match the late-time behavior of the bluest, most actively star-forming galaxies and the passive, "red and dead" galaxies. We introduce a 4-parameter model, which transitions from lin-exp to a linear ramp after a transition time, which describes our simulated galaxies very well. We test the ability of these parametrised models to recover (at z=0, 0.5, and 1) the stellar mass-to-light ratios, specific star formation rates, and stellar population ages from the galaxy colours, computed from the full SPH star formation histories using the FSPS code of Conroy et al. (2009). Fits with tau-models systematically overestimate M/L by ~ 0.2 dex, overestimate population ages by ~ 1-2 Gyr, and underestimate sSFR by ~ 0.05 dex. Fits with lin-exp are less biased on average, but the 4-parameter model yields the best results for the full range of galaxies. Marginalizing over the free parameters of the 4-parameter model leads to slightly larger statistical errors than 1-parameter fits but essentially removes all systematic biases, so this is our recommended procedure for fitting real galaxies.

preprint2014arXiv

The COS-Dwarfs Survey: The Carbon Reservoir Around sub-L* Galaxies

We report new observations of circumgalactic gas from the COS-Dwarfs survey, a systematic investigation of the gaseous halos around 43 low-mass z $\leq$ 0.1 galaxies using background QSOs observed with the Cosmic Origins Spectrograph. From the projected 1D and 2D distribution of C IV absorption, we find that C IV absorption is detected out to ~ 0.5 R$_{vir}$ of the host galaxies. The C IV absorption strength falls off radially as a power law and beyond 0.5 R$_{vir}$, no C IV absorption is detected above our sensitivity limit of ~ 50-100 m$Å$. We find a tentative correlation between detected C IV absorption strength and star formation, paralleling the strong correlation seen in highly ionized oxygen for L~L* galaxies by the COS-Halos survey. The data imply a large carbon reservoir in the CGM of these galaxies, corresponding to a minimum carbon mass of $\gtrsim$ 1.2$\times 10^6$ $M_\odot$ out to ~ 110 kpc. This mass is comparable to the carbon mass in the ISM and more than the carbon mass currently in stars of these galaxies. The C IV absorption seen around these sub-L* galaxies can account for almost two-thirds of all $W_r$> 100 m$Å$ C IV absorption detected at low z. Comparing the C IV covering fraction with hydrodynamical simulations, we find that an energy-driven wind model is consistent with the observations whereas a wind model of constant velocity fails to reproduce the CGM or the galaxy properties.

preprint2014arXiv

The COS-Halos Survey: Physical Conditions and Baryonic Mass in the Low-Redshift Circumgalactic Medium

We analyze the physical conditions of the cool, photoionized (T $\sim 10^4$ K) circumgalactic medium (CGM) using the COS-Halos suite of gas column density measurements for 44 gaseous halos within 160 kpc of $L \sim L^*$ galaxies at $z \sim 0.2$. These data are well described by simple photoionization models, with the gas highly ionized (n$_{\rm HII}$/n$_{\rm H} \gtrsim 99\%$) by the extragalactic ultraviolet background (EUVB). Scaling by estimates for the virial radius, R$_{\rm vir}$, we show that the ionization state (tracked by the dimensionless ionization parameter, U) increases with distance from the host galaxy. The ionization parameters imply a decreasing volume density profile n$_{\rm H}$ = (10$^{-4.2 \pm 0.25}$)(R/R$_{\rm vir})^{-0.8\pm0.3}$. Our derived gas volume densities are several orders of magnitude lower than predictions from standard two-phase models with a cool medium in pressure equilibrium with a hot, coronal medium expected in virialized halos at this mass scale. Applying the ionization corrections to the HI column densities, we estimate a lower limit to the cool gas mass M$_{\rm CGM}^{\rm cool} > 6.5 \times 10^{10}$ M$_{\odot}$ for the volume within R $<$ R$_{\rm vir}$. Allowing for an additional warm-hot, OVI-traced phase, the CGM accounts for at least half of the baryons purported to be missing from dark matter halos at the 10$^{12}$ M$_{\odot}$ scale.

preprint2014arXiv

The Photon Underproduction Crisis

We examine the statistics of the low-redshift Lyman-alpha forest from smoothed particle hydrodynamic simulations in light of recent improvements in the estimated evolution of the cosmic ultraviolet background (UVB) and recent observations from the Cosmic Origins Spectrograph (COS). We find that the value of the metagalactic photoionization rate required by our simulations to match the observed properties of the low-redshift Lyman-alpha forest is a factor of 5 larger than the value predicted by state-of-the art models for the evolution of this quantity. This mismatch results in the mean flux decrement of the Lyman-alpha forest being underpredicted by at least a factor of 2 (a 10-sigma discrepancy with observations) and a column density distribution of Lyman-alpha forest absorbers systematically and significantly elevated compared to observations over nearly two decades in column density. We examine potential resolutions to this mismatch and find that either conventional sources of ionizing photons (galaxies and quasars) must be significantly elevated relative to current observational estimates or our theoretical understanding of the low-redshift universe is in need of substantial revision.

preprint2013arXiv

AGN proximity zone fossils and the delayed recombination of metal lines

We model the time-dependent evolution of metal-enriched intergalactic and circumgalactic gas exposed to the fluctuating radiation field from an active galactic nucleus (AGN). We consider diffuse gas densities (n_H=10^-5-10^-3 cm^-3) exposed to the extra-galactic background (EGB) and initially in thermal equilibrium (T \sim 10^4-10^4.5 K). Once the proximate AGN field turns on, additional photo-ionisation rapidly ionises the HI and metals. The enhanced AGN radiation field turns off after a typical AGN lifetime (tau_AGN=1-20 Myr) and the field returns to the EGB intensity, but the metals remain out of ionisation equilibrium for timescales that can significantly exceed tau_AGN. We define this phase as the AGN proximity zone "fossil" phase and show that high ionisation stages (e.g. OVI, NeVIII, MgX) are in general enhanced, while the abundances of low ions are reduced. In contrast, HI re-equilibrates rapidly (<<tau_AGN) owing to its low neutral fraction at diffuse densities. We demonstrate that metal column densities of intervening gas observed in absorption in quasar sight lines are significantly affected by delayed recombination for a wide range of densities, metallicities, and AGN strengths, lifetimes, and duty cycles. We model the exceptionally strong z=0.9 NeVIII absorbers observed by Tripp et al. (2011) as arising in a possible fossil zone or near a recently turned-on AGN and we demonstrate that at low redshift even moderate strength AGN could significantly enhance the high-ion metal columns in the circumgalactic media of galaxies observed without active AGN. Fossil proximity zones may be very important during the quasar era, z \sim 2-5. AGN proximity zone fossils allow a whole new class of non-equilibrium solutions that may be applicable to a large fraction of observed metal absorbers and which could potentially change the inferred physical conditions and masses of diffuse gases.

preprint2013arXiv

Hydrogen and Metal Line Absorption Around Low-Redshift Galaxies in Cosmological Hydrodynamic Simulations

We study the physical conditions of the circum-galactic medium (CGM) around z=0.25 galaxies as traced by HI and metal line absorption, using cosmological hydrodynamic simulations that include galactic outflows. Using lines of sight targeted at impact parameters from 10 kpc to 1 Mpc around galaxies with halo masses from 10^11-10^13 M_solar, we study the physical conditions and their variation with impact parameter b and line-of-sight velocity delta v in the CGM as traced by HI, MgII, SiIV, CIV, OVI, and NeVIII absorbers. All ions show a strong excess of absorption near galaxies compared to random lines of sight. The excess continues beyond 1 Mpc, reflecting the correlation of metal absorption with large-scale structure. Absorption is particularly enhanced within about v<300 km/sec and roughly 300 kpc of galaxies (with distances somewhat larger for the highest ion), approximately delineating the CGM; this range contains the majority of global metal absorption. Low ions like MgII and SiIV predominantly arise in denser gas closer to galaxies and drop more rapidly with b, while high ions OVI and NeVIII trace more diffusely distributed gas with a comparatively flat radial profile; CIV is intermediate. All ions predominantly trace T~10^4-4.5 K photo-ionised gas at all b, but when hot CGM gas is present (mostly in larger halos), we see strong collisionally-ionised OVI and NeVIII at b <= 100 kpc. Larger halo masses generally produce more absorption, though overall the trends are not as strong as that with impact parameter. These findings arise using our favoured outflow scalings as expected for momentum-driven winds; with no winds, the CGM gas remains mostly unenriched, while our outflow model with a constant velocity and mass loading factor produce hotter, more widely dispersed metals.

preprint2013arXiv

Non-Equilibirum ionisation and cooling of metal-enriched gas in the presence of a photo-ionisation background

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)

preprint2013arXiv

The COS-Halos Survey: Rationale, Design, and A Census of Circumgalactic Neutral Hydrogen

We present the design and methods of the COS-Halos survey, a systematic investigation of the gaseous halos of 44 z = 0.15-0.35 galaxies using background QSOs observed with the Cosmic Origins Spectrograph aboard the Hubble Space Telescope. This survey has yielded 39 spectra of z_em ~ 0.5 QSOs with S/N ~ 10-15 per resolution element. The QSO sightlines pass within 150 physical kpc of the galaxies, which span early and late types over stellar mass log M* / Msun= 9.5 - 11.5. We find that the CGM exhibits strong HI, averaging 1 Ang in Lya equivalent width out to 150 kpc, with 100% covering fraction for star-forming galaxies and 75% covering for passive galaxies. We find good agreement in column densities between this survey and previous studies over similar range of impact parameter. There is weak evidence for a difference between early- and late-type galaxies in the strength and distribution of HI. Kinematics indicate that the detected material is bound to the host galaxy, such that >~90% of the detected column density is confined within +/-200 km s^-1 of the galaxies. This material generally exists well below the halo virial temperatures at T<~ 10^5 K. We evaluate a number of possible origin scenarios for the detected material, and in the end favor a simple model in which the bulk of the detected HI arises in a bound, cool, low-density photoionized diffuse medium that is generic to all L* galaxies and may harbor a total gaseous mass comparable to galactic stellar masses.

preprint2013arXiv

The effect of metal enrichment and galactic winds on galaxy formation in cosmological zoom simulations

We investigate the differential effects of metal cooling and galactic stellar winds on the cosmological formation of individual galaxies with three sets of cosmological, hydrodynamical zoom simulations of 45 halos in the mass range 10^11<M_halo<10^13M_sun. Models including both galactic winds and metal cooling (i) suppress early star formation at z>1 and predict reasonable star formation histories, (ii) produce galaxies with high cold gas fractions (30-60 per cent) at high redshift, (iii) significantly reduce the galaxy formation efficiencies for halos (M_halo<10^12M_sun) at all redshifts in agreement with observational and abundance matching constraints, (iv) result in high-redshift galaxies with reduced circular velocities matching the observed Tully-Fisher relation at z~2, and (v) significantly increase the sizes of low-mass galaxies (M_stellar<3x10^10M_sun) at high redshift resulting in a weak size evolution - a trend in agreement with observations. However, the low redshift (z<0.5) star formation rates of massive galaxies are higher than observed (up to ten times). No tested model predicts the observed size evolution for low-mass and high-mass galaxies simultaneously. Due to the delayed onset of star formation in the wind models, the metal enrichment of gas and stars is delayed and agrees well with observational constraints. Metal cooling and stellar winds are both found to increase the ratio of in situ formed to accreted stars - the relative importance of dissipative vs. dissipationless assembly. For halo masses below ~10^12M_sun, this is mainly caused by less stellar accretion and compares well to predictions from semi-analytical models but still differs from abundance matching models. For higher masses, the fraction of in situ stars is over-predicted due to the unrealistically high star formation rates at low redshifts.

preprint2013arXiv

The High-Ion Content and Kinematics of Low-Redshift Lyman Limit Systems

We study the high-ionization phase and kinematics of the circumgalactic medium around low-redshift galaxies using a sample of 23 Lyman Limit Systems (LLSs) at 0.08<z<0.93 observed with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope. In Lehner et al. (2013), we recently showed that low-z LLSs have a bimodal metallicity distribution. Here we extend that analysis to search for differences between the high-ion and kinematic properties of the metal-poor and metal-rich branches. We find that metal-rich LLSs tend to show higher O VI columns and broader O VI profiles than metal-poor LLSs. The total H I line width (dv90 statistic) in LLSs is not correlated with metallicity, indicating that the H I kinematics alone cannot be used to distinguish inflow from outflow and gas recycling. Among the 17 LLSs with O VI detections, all but two show evidence of kinematic sub-structure, in the form of O VI-H I centroid offsets, multiple components, or both. Using various scenarios for how the metallicity in the high-ion and low-ion phases of each LLS compare, we constrain the ionized hydrogen column in the O VI phase to lie in the range log N(H II)~17.6-20. The O VI phase of LLSs is a substantial baryon reservoir, with M(high-ion)~10^{8.5-10.9}(r/150 kpc)^2 solar masses, similar to the mass in the low-ion phase. Accounting for the O VI phase approximately doubles the contribution of low-z LLSs to the cosmic baryon budget.

preprint2013arXiv

The Neutral Hydrogen Content of Galaxies in Cosmological Hydrodynamic Simulations

We examine the global HI properties of galaxies in quarter-billion particle cosmological simulations using Gadget-2, focusing on how galactic outflows impact HI content. We consider four outflow models, including a new one (ezw) motivated by recent interstellar medium simulations in which the wind speed and mass loading factor scale as expected for momentum-driven outflows for larger galaxies and energy-driven outflows for dwarfs (sigma<75 km/s). To obtain predicted HI masses, we employ a simple but effective local correction for particle self-shielding, and an observationally-constrained transition from neutral to molecular hydrogen. Our ezw simulation produces an HI mass function whose faint-end slope of -1.3 agrees well with observations from the ALFALFA survey; other models agree less well. Satellite galaxies have a bimodal distribution in HI fraction versus halo mass, with smaller satellites and/or those in larger halos more often being HI-deficient. At a given stellar mass, HI content correlates with star formation rate and inversely correlates with metallicity, as expected if driven by stochasticity in the accretion rate. To higher redshifts, massive HI galaxies disappear and the mass function steepens. The global cosmic HI density conspires to remain fairly constant from z~5-0, but the relative contribution from smaller galaxies increases with redshift.

preprint2013arXiv

Tracing Inflows and Outflows with Absorption Lines in Circumgalactic Gas

We examine how HI and metal absorption lines within low-redshift galaxy halos trace the dynamical state of circumgalactic gas, using cosmological hydrodynamic simulations that include a well-vetted heuristic model for galactic outflows. We categorize inflowing, outflowing, and ambient gas based on its history and fate as tracked in our simulation. Following our earlier work showing that the ionisation level of absorbers was a primary factor in determining the physical conditions of absorbing gas, we show here that it is also a governing factor for its dynamical state. Low-ionisation metal absorbers (e.g. MgII) tend to arise in gas that will fall onto galaxies within several Gyr, while high-ionisation metal absorbers (e.g. OVI) generally trace material that was deposited by outflows many Gyr ago. Inflowing gas is dominated by enriched material that was previously ejected in an outflow, hence accretion at low redshifts is typically substantially enriched. Recycling wind material is preferentially found closer to galaxies, and is more dominant in lower-mass halos since high-mass halos have more hot gas that is able to support itself against infall. Low-mass halos also tend to re-eject more of their accreted material, owing to our outflow prescription that employs higher mass loading factors for lower-mass galaxies. Typical HI absorbers trace unenriched ambient material that is not participating in the baryon cycle, but stronger HI absorbers arise in cool, enriched inflowing gas. Instantaneous radial velocity measures of absorbers are generally poor at distinguishing between inflowing and outflowing gas, except in the case of very recent outflows. These results suggest that probing halo gas using a range of absorbers can provide detailed information about the amount and physical conditions of material that is participating in the baryon cycle.

preprint2012arXiv

A fundamental problem in our understanding of low mass galaxy evolution

Recent studies have found a dramatic difference between the observed number density evolution of low mass galaxies and that predicted by semi-analytic models. 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]

preprint2012arXiv

Not Dead Yet: Cool Circumgalactic Gas in the Halos of Early Type Galaxies

We report new observations of circumgalactic gas in the halos of early type galaxies obtained by the COS-Halos Survey with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope. We find that detections of HI surrounding early type galaxies are typically as common and strong as around star-forming galaxies, implying that the total mass of circumgalactic material is comparable in the two populations. For early type galaxies, the covering fraction for HI absorption above 10^16 cm^2 is ~40-50% within ~150 kpc. Line widths and kinematics of the detected material show it to be cold (T ~< 10^5 K) in comparison to the virial temperature of the host halos. The implied masses of cool, photoionized CGM baryons may be up to 10^9 --- 10^11 Msun. Contrary to some theoretical expectations, strong halo HI absorbers do not disappear as part of the quenching of star-formation. Even passive galaxies retain significant reservoirs of halo baryons which could replenish the interstellar gas reservoir and eventually form stars. This halo gas may feed the diffuse and molecular gas that is frequently observed inside ETGs.

preprint2012arXiv

Testing Subhalo Abundance Matching in Cosmological Smoothed Particle Hydrodynamics Simulations

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 < r < 10 Mpc/h. For the wind simulation, agreement is better than 15% at r > 2 Mpc/h, but overpopulation of massive halos increases the correlation function by a factor of ~2.5 on small scales.

preprint2011arXiv

An Analytic Model for the Evolution of the Stellar, Gas, and Metal Content of Galaxies

We present an analytic formalism that describes the evolution of the stellar, gas, and metal content of galaxies. It is based on the idea, inspired by hydrodynamic simulations, that galaxies live in a slowly-evolving equilibrium between inflow, outflow, and star formation. We argue that this formalism broadly captures the behavior of galaxy properties evolving in simulations. The resulting equilibrium equations for the star formation rate, gas fraction, and metallicity depend on three key free parameters that represent ejective feedback, preventive feedback, and re-accretion of ejected material. We schematically describe how these parameters are constrained by models and observations. Galaxies perturbed off the equilibrium relations owing to inflow stochasticity tend to be driven back towards equilibrium, such that deviations in star formation rate at a given mass are correlated with gas fraction and anti-correlated with metallicity. After an early gas accumulation epoch, quiescently star-forming galaxies are expected to be in equilibrium over most of cosmic time. The equilibrium model provides a simple intuitive framework for understanding the cosmic evolution of galaxy properties, and centrally features the cycle of baryons between galaxies and surrounding gas as the driver of galaxy growth.

preprint2011arXiv

Galaxy Evolution in Cosmological Simulations With Outflows I: Stellar Masses and Star Formation Rates

We examine the growth of the stellar content of galaxies from z=3-0 in cosmological hydrodynamic simulations incorporating parameterised galactic outflows. Without outflows, galaxies overproduce stellar masses (M*) and star formation rates (SFRs) compared to observations. Winds introduce a three-tier form for the galaxy stellar mass and star formation rate functions, where the middle tier depends on differential (i.e. mass-dependent) recycling of ejected wind material back into galaxies. A tight M*-SFR relation is a generic outcome of all these simulations, and its evolution is well-described as being powered by cold accretion, although current observations at z>2 suggest that star formation in small early galaxies must be highly suppressed. Roughly one-third of z=0 galaxies at masses below M^* are satellites, and star formation in satellites is not much burstier than in centrals. All models fail to suppress star formation and stellar mass growth in massive galaxies at z<2, indicating the need for an external quenching mechanism such as black hole feedback. All models also fail to produce dwarfs as young and rapidly star-forming as observed. An outflow model following scalings expected for momentum-driven winds broadly matches observed galaxy evolution around M^* from z=0-3, which is a significant success since these galaxies dominate cosmic star formation, but the failures at higher and lower masses highlight the challenges still faced by this class of models. We argue that central star-forming galaxies are well-described as living in a slowly-evolving equilibrium between inflows from gravity and recycled winds, star formation, and strong and ubiquitous outflows that regulate how much inflow forms into stars. Star-forming galaxy evolution is thus primarily governed by the continual cycling of baryons between galaxies and intergalactic gas.

preprint2011arXiv

Galaxy Evolution in Cosmological Simulations with Outflows II: Metallicities and Gas Fractions

We use cosmological hydrodynamic simulations to investigate how inflows, star formation, and outflows govern the the gaseous and metal content of galaxies. In our simulations, galaxy metallicities are established by a balance between inflows and outflows as governed by the mass outflow rate, implying that the mass-metallicity relation reflects how the outflow rate varies with stellar mass (M*). Gas content is set by a competition between inflow into and gas consumption within the ISM, the latter being governed by the SF law, while the former is impacted by both wind recycling and preventive feedback. Stochastic variations in the inflow rate move galaxies off the equilibrium M*-Z and Z*-fgas relations in a manner correlated with star formation rate, and the scatter is set by the timescale to re-equilibrate. The evolution of both relations from z=3-0 is slow, as individual galaxies tend to evolve mostly along the relations. Gas fractions at a given M* slowly decrease with time because the cosmic inflow rate diminishes faster than the consumption rate, while metallicities slowly increase as infalling gas becomes more enriched. Observations from z~3-0 are better matched by simulations employing momentum-driven wind scalings rather than constant wind speeds, but all models predict too low gas fractions at low masses and too high metallicities at high M*. All our models reproduce observed second-parameter trends of the mass-metallicity relation with star formation rate and environment, indicating that these are a consequence of equilibrium and not feedback. Overall, the analytical framework of our equilibrium scenario broadly captures the relevant physics establishing the galaxy gas and metal content in simulations, which suggests that the cycle of baryonic inflows and outflows centrally governs the cosmic evolution of these properties in typical star-forming galaxies.

preprint2011arXiv

Intergalactic Dust Extinction in Hydrodynamic Cosmological Simulations

Recently Menard et al. detected a subtle but systematic change in the mean color of quasars as a function of their projected separation from foreground galaxies, extending to comoving separations of ~10Mpc/h, which they interpret as a signature of reddening by intergalactic dust. We present theoretical models of this remarkable observation, using SPH cosmological simulations of a (50Mpc/h)^3 volume. Our primary model uses a simulation with galactic winds and assumes that dust traces the intergalactic metals. The predicted galaxy-dust correlation function is similar in form to the galaxy-mass correlation function, and reproducing the MSFR data requires a dust-to-metal mass ratio of 0.24, about half the value in the Galactic ISM. Roughly half of the reddening arises in dust that is more than 100Kpc/h from the nearest massive galaxy. We also examine a simulation with no galactic winds, which predicts a much smaller fraction of intergalactic metals (3% vs. 35%) and therefore requires an unphysical dust-to-metal ratio of 2.18 to reproduce the MSFR data. In both models, the signal is dominated by sightlines with E(g-i)=0.001-0.1. The no-wind simulation can be reconciled with the data if we also allow reddening to arise in galaxies up to several x 10^10 Msun. The wind model predicts a mean visual extinction of A_V ~0.0133 mag out to z=0.5, with a sightline-to-sightline dispersion similar to the mean, which could be significant for future supernova cosmology studies. Reproducing the MSFR results in these simulations requires that a large fraction of ISM dust survive its expulsion from galaxies and its residence in the intergalactic medium. Future observational studies that provide higher precision and measure the dependence on galaxy type and environment will allow detailed tests for models of enriched galactic outflows and the survival of IG dust.

preprint2011arXiv

Short-lived star-forming giant clumps in cosmological simulations of z~2 disks

Many observed massive star-forming z\approx2 galaxies are large disks that exhibit irregular morphologies, with \sim1kpc, \sim10^(8-10)Msun clumps. We present the largest sample to date of high-resolution cosmological SPH simulations that zoom-in on the formation of individual M*\sim10^(10.5)Msun galaxies in \sim10^(12)Msun halos at z\approx2. Our code includes strong stellar feedback parameterized as momentum-driven galactic winds. This model reproduces many characteristic features of this observed class of galaxies, such as their clumpy morphologies, smooth and monotonic velocity gradients, high gas fractions (f_g\sim50%) and high specific star-formation rates (\gtrsim1Gyr^(-1)). In accord with recent models, giant clumps (Mclump\sim(5x10^8-10^9)Msun) form in-situ via gravitational instabilities. However, the galactic winds are critical for their subsequent evolution. The giant clumps we obtain are short-lived and are disrupted by wind-driven mass loss. They do not virialise or migrate to the galaxy centers as suggested in recent work neglecting strong winds. By phenomenologically implementing the winds that are observed from high-redshift galaxies and in particular from individual clumps, our simulations reproduce well new observational constraints on clump kinematics and clump ages. In particular, the observation that older clumps appear closer to their galaxy centers is reproduced in our simulations, as a result of inside-out formation of the disks rather than inward clump migration.

preprint2011arXiv

The intergalactic medium over the last 10 billion years II: Metal-line absorption and physical conditions

We investigate the metallicity evolution and content of the intergalactic medium (IGM) and galactic halo gas from z=2->0 using 110-million particle cosmological hydrodynamic simulations. We focus on the detectability and physical properties of UV resonance metal-line absorbers observable with Hubble's Cosmic Origins Spectrograph (COS). We confirm that galactic superwind outflows are required to enrich the IGM to observed levels down to z=0 using three wind prescriptions contrasted to a no-wind simulation. Our favoured momentum-conserved wind prescription deposits metals closer to galaxies owing to its moderate energy input, while the more energetic constant wind model enriches the warm-hot IGM 6.4x more. Despite these significant differences, all wind models produce metal-line statistics within a factor of two of existing observations. This is because OVI, CIV, SiIV, and NeVIII absorbers primarily arise from T<10^5 K, photo-ionised gas that is enriched to similar levels in the three feedback schemes. OVI absorbers trace the diffuse phase with delta<100, which is enriched to ~1/50 Zsol at z=0, although the absorbers themselves usually exceed 0.3 Zsol and arise from inhomogeneously distributed, un-mixed winds. CIV and SiIV absorbers trace primarily T~10^4 K gas inside haloes. We predict COS will observe a population of NeVIII photo-ionised absorbers tracing T<10^5 K, delta~10 gas. MgX and SiXII are rarely detected in COS S/N=30 simulated sight lines although detections trace T=10^(6-7) K halo gas. In general, the IGM is enriched in an outside-in manner, where wind-blown metals released at higher redshift reach lower overdensities, resulting in higher ionisation species tracing lower-density, older metals. At z=0, the 90% of baryons outside of galaxies are enriched to 0.096 Zsol, but the 65% of unbound baryons in the IGM have 0.018 Zsol and contain only 4% of all metals.

preprint2011arXiv

The Large, Oxygen-Rich Halos of Star-Forming Galaxies Are A Major Reservoir of Galactic Metals

The circumgalactic medium (CGM) is fed by galaxy outflows and accretion of intergalactic gas, but its mass, heavy element enrichment, and relation to galaxy properties are poorly constrained by observations. In a survey of the outskirts of 42 galaxies with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope, we detected ubiquitous, large (150 kiloparsec) halos of ionized oxygen surrounding star-forming galaxies, but we find much less ionized oxygen around galaxies with little or no star formation. This ionized CGM contains a substantial mass of heavy elements and gas, perhaps far exceeding the reservoirs of gas in the galaxies themselves. It is a basic component of nearly all star-forming galaxies that is removed or transformed during the quenching of star formation and the transition to passive evolution.

preprint2010arXiv

A Search for Oxygen in the Low-Density Lyman-alpha Forest Using the Sloan Digital Sky Survey

We use 2167 Sloan Digital Sky Survey (SDSS) quasar spectra to search for low-density oxygen in the Intergalactic Medium (IGM). Oxygen absorption is detected on a pixel-by-pixel basis by its correlation with Lyman-alpha forest absorption. We have developed a novel Locally Calibrated Pixel (LCP) search method that uses adjacent regions of the spectrum to calibrate interlopers and spectral artifacts, which would otherwise limit the measurement of OVI absorption. Despite the challenges presented by searching for weak OVI within the Lyman-alpha forest in spectra of moderate resolution and signal-to-noise, we find a highly significant detection of absorption by oxygen at 2.7 < z < 3.2 (the null hypothesis has a chi^2=80 for 9 data points). We interpret our results using synthetic spectra generated from a lognormal density field assuming a mixed quasar-galaxy photoionizing background (Haardt & Madau 2001) and that it dominates the ionization fraction of detected OVI. The LCP search data can be fit by a constant metallicity model with [O/H] = -2.15_(-0.09)^(+0.07), but also by models in which low-density regions are unenriched and higher density regions have a higher metallicity. The density-dependent enrichment model by Aguirre et al. (2008) is also an acceptable fit. All our successful models have similar mass-weighted oxygen abundance, corresponding to [<O/H>_MW] = -2.45+-0.06. This result can be used to find the cosmic oxygen density in the Lyman-alpha forest, Omega_(Oxy, IGM) = 1.4(+-0.2)x10^(-6) = 3x10^(-4) Omega_b. This is the tightest constraint on the mass-weighted mean oxygen abundance and the cosmic oxygen density in the Lyman-alpha forest to date and indicates that it contains approximately 16% of metals produced by star formation (Bouché et al. 2008) up to z = 3.

preprint2010arXiv

Feedback and Recycled Wind Accretion: Assembling the z=0 Galaxy Mass Function

We analyse cosmological hydrodynamic simulations that include observationally-constrained prescriptions for galactic outflows. If these simulated winds accurately represent winds in the real Universe, then material previously ejected in winds provides the dominant source of gas infall for new star formation at redshifts z<1. This recycled wind accretion, or wind mode, provides a third physically distinct accretion channel in addition to the "hot" and "cold" modes emphasised in recent theoretical studies. Because of the interaction between outflows and gas in and around halos, the recycling timescale of wind material (t_rec) is shorter in higher-mass systems, which reside in denser gaseous environments. In these simulations, this differential recycling plays a central role in shaping the present-day galaxy stellar mass function (GSMF). If we remove all particles that were ever ejected in a wind, then the predicted GSMFs are much steeper than observed; galaxy masses are suppressed both by the direct removal of gas and by the hydrodynamic heating of their surroundings, which reduces subsequent infall. With wind recycling included, the simulation that incorporates our favoured momentum-driven wind scalings reproduces the observed GSMF for stellar masses 10^9 < M < 5x10^10 Msolar. At higher masses, wind recycling leads to excessive galaxy masses and excessive star formation rates relative to observations. In these massive systems, some quenching mechanism must suppress the re-accretion of gas ejected from star-forming galaxies. In short, as has long been anticipated, the form of the GSMF is governed by outflows; the unexpected twist here for our simulated winds is that it is not primarily the ejection of material but how the ejected material is re-accreted that governs the GSMF.

preprint2010arXiv

Galactic outflows and the kinematics of damped Lyman alpha absorbers

The kinematics of damped Lyman alpha absorbers (DLAs) are difficult to reproduce in hierarchical galaxy formation models, particularly the preponderance of wide systems. We investigate DLA kinematics at z=3 using high-resolution cosmological hydrodynamical simulations that include a heuristic model for galactic outflows. Without outflows, our simulations fail to yield enough wide DLAs, as in previous studies. With outflows, predicted DLA kinematics are in much better agreement with observations. Comparing two outflow models, we find that a model based on momentum-driven wind scalings provides the best match to the observed DLA kinematic statistics of Prochaska & Wolfe. In this model, DLAs typically arise a few kpc away from galaxies that would be identified in emission. Narrow DLAs can arise from any halo and galaxy mass, but wide ones only arise in halos with mass >10^11 Mo, from either large central or small satellite galaxies. This implies that the success of this outflow model originates from being most efficient at pushing gas out from small satellite galaxies living in larger halos. This increases the cross-section for large halos relative to smaller ones, thereby yielding wider kinematics. Our simulations do not include radiative transfer effects or detailed metal tracking, and outflows are modeled heuristically, but they strongly suggest that galactic outflows are central to understanding DLA kinematics. An interesting consequence is that DLA kinematics may place constraints on the nature and efficiency of gas ejection from high-z galaxies.

preprint2010arXiv

Smoothly-Rising Star Formation Histories During the Reionization Epoch

Cosmological hydrodynamic simulations robustly predict that high-redshift galaxy star formation histories (SFHs) are smoothly-rising and vary with mass only by a scale factor. We use our latest simulations to test whether this scenario can account for recent observations at z>=6 from WFC3/IR, NICMOS, and IRAC. Our simulations broadly reproduce the observed ultraviolet (UV) luminosity functions and stellar mass densities and their evolution at z=6-8, all of which are nontrivial tests of the mean SFH. In agreement with observations, simulated galaxies possess blue UV continua owing to young ages (50-150 Myr), low metallicities (0.1-0.5 Zsun), and low dust columns (E(B-V) <= 0.05). Observations imply a near-unity slope in the stellar mass--star formation rate relation at all z=6-8, confirming the prediction that SFH shapes are invariant. Current surveys detect the majority of galaxies with stellar masses exceeding 10^9 Msun and few galaxies less massive than 10^{8.5} Msun, implying that they probe no more than the brightest 30% of the complete star formation and stellar mass densities at z>=6. Finally, we demonstrate that there is no conflict between smoothly-rising SFHs and recent clustering observations. This is because momentum-driven outflows suppress star formation in low-mass halos, leading to overall occupancies of 0.2-0.4 even though the star formation duty cycle is one. This leads to many interesting predictions at z>=4, among them that (1) optically-selected and UV-selected samples largely overlap; (2) few galaxies exhibit significantly suppressed specific star formation rates; and (3) occupancy is constant or increasing with decreasing luminosity. These predictions are in tentative agreement with current observations, but further analysis of existing and upcoming data sets is required in order to test them more thoroughly. (abridged)

preprint2010arXiv

The intergalactic medium over the last 10 billion years I: Lyman alpha absorption and physical conditions

The intergalactic medium (IGM) is the dominant reservoir of baryons at all cosmic epochs. We investigate the evolution of the IGM from z=2-0 in 48 Mpc/h, 110-million particle cosmological hydrodynamic simulations using three prescriptions for galactic outflows. We focus on the evolution of IGM physical properties, and how such properties are traced by Ly-alpha absorption as detectable using HST/COS. Our results broadly confirm the canonical picture that most Ly-alpha absorbers arise from highly ionized gas tracing filamentary large-scale structure. Growth of structure causes gas to move from the diffuse photoionized IGM into other cosmic phases, namely stars, cold and hot gas within galaxy halos, and the unbound and shock-heated warm-hot intergalactic medium (WHIM). By today, baryons are roughly equally divided between bound phases (35%), the diffuse IGM (41%), and the WHIM (24%). Here we (re)define the WHIM as gas with overdensities lower than that in halos and temperatures >10^5 K, in order to more closely align it with "missing baryons". When we tune our photoionizing background to match the observed evolution of the Ly-alpha mean flux decrement, we obtain a line count evolution that broadly agrees with available data. We predict a column density distribution slope of -1.70 for our favored momentum-driven wind model, in agreement with recent observations, and it becomes shallower with redshift. With improved statistics, the frequency of strong lines can be a valuable diagnostic of outflows, and our favored wind model matches existing data best among our models. The relationship between column density and physical density is fairly tight from z=2-0, and evolves as rho N_HI^0.74 10^(-0.37z) for diffuse absorbers. Linewidths only loosely reflect the temperature of the absorbing gas, which will hamper attempts to quantify the WHIM using broad Ly-alpha absorbers. [Abridged]

preprint2010arXiv

The Nature of Sub-millimetre Galaxies in Cosmological Hydrodynamic Simulations

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

preprint2009arXiv

The Simulated HI Sky at low redshift

Observations of intergalactic neutral hydrogen can provide a wealth of information about structure and galaxy formation, potentially tracing accretion and feedback processes on Mpc scales. Below a column density of NHI ~ 10^19 cm-2, the "edge" or typical observational limit for HI emission from galaxies, simulations predict a cosmic web of extended emission and filamentary structures. We study the distribution of neutral hydrogen and its 21cm emission properties in a cosmological hydrodynamic simulation, to gain more insights into the distribution of HI below NHI ~ 10^19 cm-2. Such Lyman Limit systems are expected to trace out the cosmic web, and are relatively unexplored. Beginning with a 32 h^-1 Mpc simulation, we extract the neutral hydrogen component by determining the neutral fraction, including a post-processed correction for self-shielding based on the thermal pressure. We take into account molecular hydrogen, assuming an average density ratio Omega_H2 / Omega_HI = 0.3 at z = 0. The statistical properties of the HI emission are compared with observations, to assess the reliability of the simulation. The simulated HI distribution robustly describes the full column density range between NHI ~ 10^14 and NHI ~ 10^21 cm-2 and agrees very well with available measurements from observations. Furthermore there is good correspondence in the statistics when looking at the two-point correlation function and the HI mass function. The reconstructed maps are used to simulate observations of existing and future telescopes by adding noise and taking account of the sensitivity of the telescopes. The general agreement in statistical properties of HI suggests that neutral hydrogen as modeled in this hydrodynamic simulation is a fair representation of that in the Universe. (abridged)

preprint2009arXiv

Tracing the Reionization-Epoch Intergalactic Medium with Metal Absorption Lines

IGM metal absorption lines observed in z>6 spectra offer the opportunity to probe early feedback processes, the nature of enriching sources, and the topology of reionization. We run high-resolution cosmological simulations including galactic outflows to study the observability and physical properties of 5 ions (C II, C IV, O I, Si II, Si IV) in absorption between z=8->5. We apply three cases for ionization conditions: Fully neutral, fully reionized, and a patchy model based on the flux from the nearest galaxy. We find that our simulations broadly fit available z~5-6 IGM metal-line data, although all observations cannot be accommodated with a single ionization condition. Variations in O I absorbers among sight lines seen by Becker et al. (2006) suggest significant neutral IGM patches down to z~6. Strong C IV absorbers at z~6 may be the result of ionization by their parent galaxy. Our outflows have typical speeds of ~200 km/s and mass loading factors of ~6. Such high mass loading is critical for enriching the IGM to the observed levels while curtailing star formation to match the observed z~6 rest-frame UV luminosity function. The volume filling factor of metals increases during this epoch, but only reaches ~1% for Z>10^(-3) Zsolar by z=5. C IV is an ideal tracer of IGM metals at z~5-6, with dropping global ionization fractions to either higher or lower redshifts. This results in a strongly increasing global Omega(C IV) from z=8->5, in contrast to its relative constancy from z=5->2. Our simulations do not support widespread early IGM enrichment from e.g. Pop III stars. High-z absorbers arise from metals on their first outward journey from galaxies, at distances less than 50 kpc. The galaxies responsible for early IGM enrichment have typical M*=10^(7.0-8.5) Msolar.

preprint2008arXiv

The Flux Auto- and Cross-Correlation of the Lyman-alpha Forest. II. Modelling Anisotropies with Cosmological Hydrodynamic Simulations

The isotropy of the Lyman-alpha forest in real-space uniquely provides a measurement of cosmic geometry at z > 2. The angular diameter distance for which the correlation function along the line of sight and in the transverse direction agree corresponds to the correct cosmological model. However, the Lyman-alpha forest is observed in redshift-space where distortions due to Hubble expansion, bulk flows, and thermal broadening introduce anisotropy. Similarly, a spectrograph's line spread function affects the autocorrelation and cross-correlation differently. In this the second paper of a series on using the Lyman-alpha forest observed in pairs of QSOs for a new application of the Alcock-Paczynski (AP) test, these anisotropies and related sources of potential systematic error are investigated with cosmological hydrodynamic simulations. Three prescriptions for galactic outflow were compared and found to have only a marginal effect on the Lyman-alpha flux correlation (which changed by at most 7% with use of the currently favored variable-momentum wind model vs. no winds at all). An approximate solution for obtaining the zero-lag cross-correlation corresponding to arbitrary spectral resolution directly from the zero-lag cross-correlation computed at full-resolution (good to within 2% at the scales of interest) is presented. Uncertainty in the observationally determined mean flux decrement of the Lyman-alpha forest was found to be the dominant source of systematic error; however, this is reduced significantly when considering correlation ratios. We describe a simple scheme for implementing our results, while mitigating systematic errors, in the context of a future application of the AP test.

Benjamin D. Oppenheimer

What is connected

Connect this record

See the researcher in context

Building this map preview

55 published item(s)

X-ray Absorption Lines in the Warm-Hot Intergalactic Medium: Probing Chandra observations with the CAMEL simulations

Snowmass2021 CMB-HD White Paper

The COS CGM Compendium. IV. Effects of Varying Ionization Backgrounds on Metallicity Determinations in the z < 1 Circumgalactic Medium

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

CMB-HD: Astro2020 RFI Response

EAGLE and Illustris-TNG predictions for resolved eROSITA X-ray observations of the circumgalactic medium around normal galaxies

The Impact of Wind Scalings on Stellar Growth and the Baryon Cycle in Cosmological Simulations

The Robustness of Cosmological Hydrodynamic Simulation Predictions to Changes in Numerics and Cooling Physics

The warm-hot circumgalactic medium around EAGLE-simulation galaxies and its detection prospects with X-ray and UV line absorption

CMB-HD: An Ultra-Deep, High-Resolution Millimeter-Wave Survey Over Half the Sky

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

The COS CGM Compendium. III: Metallicity and Physical Properties of the Cool Circumgalactic Medium at z<1

The lensing properties of subhaloes in massive elliptical galaxies in sterile neutrino cosmologies

Bimodality of low-redshift circumgalactic O VI in non-equilibrium EAGLE zoom simulations

Cosmic distribution of highly ionized metals and their physical conditions in the EAGLE simulations

Baryon Cycling in the Low-Redshift Circumgalactic Medium: A Comparison of Simulations to the COS-Halos Survey

Nearby galaxy filaments and the Lya forest: confronting simulations and the UV background with observations

The Growth and Enrichment of the Intragroup Gas

The impact of environment and mergers on the HI content of galaxies in hydrodynamic simulations

Torque-Limited Growth of Massive Black Holes in Galaxies Across Cosmic Time

A Budget and Accounting of Metals at z~0: Results from the COS-Halos Survey

Cosmological Zoom Simulations of z = 2 Galaxies: The Impact of Galactic Outflows

Non-equilibrium chemistry and cooling in the diffuse interstellar medium - I. Optically thin regime

Non-equilibrium chemistry and cooling in the diffuse interstellar medium - II. Shielded gas

Parametrising Star Formation Histories

The COS-Dwarfs Survey: The Carbon Reservoir Around sub-L* Galaxies

The COS-Halos Survey: Physical Conditions and Baryonic Mass in the Low-Redshift Circumgalactic Medium

The Photon Underproduction Crisis

AGN proximity zone fossils and the delayed recombination of metal lines

Hydrogen and Metal Line Absorption Around Low-Redshift Galaxies in Cosmological Hydrodynamic Simulations

Non-Equilibirum ionisation and cooling of metal-enriched gas in the presence of a photo-ionisation background

The COS-Halos Survey: Rationale, Design, and A Census of Circumgalactic Neutral Hydrogen

The effect of metal enrichment and galactic winds on galaxy formation in cosmological zoom simulations

The High-Ion Content and Kinematics of Low-Redshift Lyman Limit Systems

The Neutral Hydrogen Content of Galaxies in Cosmological Hydrodynamic Simulations

Tracing Inflows and Outflows with Absorption Lines in Circumgalactic Gas

A fundamental problem in our understanding of low mass galaxy evolution

Not Dead Yet: Cool Circumgalactic Gas in the Halos of Early Type Galaxies

Testing Subhalo Abundance Matching in Cosmological Smoothed Particle Hydrodynamics Simulations

An Analytic Model for the Evolution of the Stellar, Gas, and Metal Content of Galaxies

Galaxy Evolution in Cosmological Simulations With Outflows I: Stellar Masses and Star Formation Rates

Galaxy Evolution in Cosmological Simulations with Outflows II: Metallicities and Gas Fractions

Intergalactic Dust Extinction in Hydrodynamic Cosmological Simulations

Short-lived star-forming giant clumps in cosmological simulations of z~2 disks

The intergalactic medium over the last 10 billion years II: Metal-line absorption and physical conditions

The Large, Oxygen-Rich Halos of Star-Forming Galaxies Are A Major Reservoir of Galactic Metals

A Search for Oxygen in the Low-Density Lyman-alpha Forest Using the Sloan Digital Sky Survey

Feedback and Recycled Wind Accretion: Assembling the z=0 Galaxy Mass Function

Galactic outflows and the kinematics of damped Lyman alpha absorbers

Smoothly-Rising Star Formation Histories During the Reionization Epoch

The intergalactic medium over the last 10 billion years I: Lyman alpha absorption and physical conditions

The Nature of Sub-millimetre Galaxies in Cosmological Hydrodynamic Simulations

The Simulated HI Sky at low redshift

Tracing the Reionization-Epoch Intergalactic Medium with Metal Absorption Lines

The Flux Auto- and Cross-Correlation of the Lyman-alpha Forest. II. Modelling Anisotropies with Cosmological Hydrodynamic Simulations