Benjamin D. Oppenheimer

Cosmological simulations of intergalactic medium enrichment from galactic outflows

We investigate models of self-consistent chemical enrichment of the intergalactic medium (IGM) from z = 6.0 → 1.5, based on hydrodynamic simulations of structure formation that explicitly incorporate outflows from star-forming galaxies. Our main result is that outflow parametrizations derived from observations of local starburst galaxies, in particular momentum-driven wind scenarios, provide the best agreement with observations of C iv absorption at z ∼ 2-5. Such models sufficiently enrich the high-z IGM to produce a global mass density of Civ absorbers that is relatively invariant from z = 5.5 → 1.5, in agreement with observations. This occurs despite continual IGM enrichment causing an increase in volume-averaged metallicity by ∼ x 5-10 over this redshift range, because energy input accompanying the enriching outflows causes a drop in the global ionization fraction of C iv. Comparisons to observed C IV column density and linewidth distributions and C IV-based pixel optical depth ratios provide significant constraints on wind models. Our best-fitting outflow models show mean IGM temperatures only slightly above our no-outflow case, metal filling factors of just a few per cent with volume-weighted metallicities around 10 -3 at z ∼ 3, significant amounts of collisionally ionized C iv absorption and a metallicity-density relationship that rises rapidly at low overdensities and flattens at higher ones. In general, we find that outflow speeds must be high enough to enrich the low-density IGM at early times but low enough not to overheat it, and concurrently must significantly suppress early star formation while still producing enough early metals. It is therefore non-trivial that locally calibrated momentum-driven wind scenarios naturally yield the desired strength and evolution of outflows, and suggest that such models represent a significant step towards understanding the impact of galactic outflows on galaxies and the IGM across cosmic time.

Mass, Metal, and Energy Feedback in Cosmological Simulations

Using GADGET-2 cosmological hydrodynamic simulations including an observationally constrained model for galactic outflows, we investigate how feedback from star formation distributes mass, metals, and energy on cosmic scales from z = 6 → 0. We include instantaneous enrichment from Type II supernovae (SNe), as well as delayed enrichment from Type Ia SNe and stellar [asymptotic giant branch (AGB)] mass loss, and we individually track carbon, oxygen, silicon and iron using the latest yields. Following on the success of the momentum-driven wind scalings, we improve our implementation by using an on-the-fly galaxy finder to derive wind properties based on host galaxy masses. By tracking wind particles in a suite of simulations, we find: (1) wind material re-accretes on to a galaxy (usually the same one it left) on a recycling time-scale that varies inversely with galaxy mass (e.g. <1 Gyr for L ∗ galaxies at

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 behaviour 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 reaccretion 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 anticorrelated 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.

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

Observations with the Hubble Space Telescope show that halos of ionized gas are common around star-forming galaxies. 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; we found 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. Our data indicate that 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.

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 within a hierarchical structure formation context. 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. Gas content, meanwhile, is set by a competition between inflow into and gas consumption within the interstellar medium, the latter being governed by the star formation law, while the former is impacted by both wind recycling and preventive feedback. Stochastic variations in the inflow rate move galaxies off the equilibrium mass-metallicity and mass-gas fraction 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 stellar mass 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 masses. 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.

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 parametrized 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 the 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 the star formation in small early galaxies must be highly suppressed. Roughly, one-third of z = 0 galaxies at masses below M * are satellites and the star formation in satellites is not much burstier than in centrals. All models fail to suppress the 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 the observed galaxy evolution around M * from z = 0 to 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.

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 (� < 75 kms 1 ). 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.

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

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

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

Many observed massive star-forming z 2 galaxies are large disks that exhibit irregular morphologies, with 1 kpc, 108-1010M⊙ clumps. We present the largest sample to date of high-resolution cosmological smoothed particle hydrodynamics simulations that zoom-in on the formation of individual M * 1010.5M⊙ galaxies in 1012M⊙ halos at z 2. 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 50%), and high specific star formation rates (1 Gyr–1). In accord with recent models, giant clumps (M clump (5 × 108-109)M⊙) 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 virialize 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.

The nature and origin of low-redshift O vi absorbers

The O VI ion observed in quasar absorption-line spectra is the most accessible tracer of the cosmic metal distribution in the low-redshift (z < 0.5) intergalactic medium (IGM). We explore the nature and origin of O VI absorbers using cosmological hydrodynamic simulations including galactic outflows with a range of strengths. We consider the effects of ionization background variations, non-equilibrium ionization and cooling, uniform metallicity and small-scale (sub-resolution) turbulence. Our main results are as follows. (1) IGM O VI is predominantly photoionized with T ≈ 10 4.2±0.2 K. A key reason for this is that O VI absorbers preferentially trace overenriched (by ∼ x 5) regions of the IGM at a given density, which enhances metal-line cooling such that absorbers can cool to photo-ionized temperatures within a Hubble time. As such, O VI is not a good tracer of the warm-hot intergalactic medium. (2) The predicted O VI properties fit observables if and only if sub-resolution turbulence is added, regardless of any other model variations. The required turbulence increases with O VI absorber strength. Stronger absorbers arise from more recent outflows, so qualitatively this can be understood if IGM turbulence dissipates on the order of a Hubble time. The amount of turbulence is consistent with other examples of turbulence observed in the IGM and galactic haloes. (3) Metals traced by O VI and H 1 do not trace exactly the same baryons, but reside in the same large-scale structure. Our simulations reproduce observed alignment statistics between O VI and H I, yet aligned absorbers typically have O VI arising from cooler gas, and for stronger absorbers lower densities, than H I. Owing to peculiar velocities dominating the line structure, coincident absorption often arises from spatially distinct gas. (4) Photo-ionized O VI traces gas in a variety of environments, and is not directly associated with the nearest galaxy, though is typically nearest to ∼0.1 L * galaxies. Weaker O VI components trace some of the oldest cosmic metals. (5) Very strong absorbers (EW ≤100 mA) are more likely to be collisionally ionized, tracing more recent enrichment (≥2 Gyr) within or near galactic haloes.