How Do Galaxies Get Their Gas?

We examine the temperature history of gas accreted by forming galaxies in SPH simulations. About half the gas shock heats to roughly the virial temperature of the galaxy potential well before cooling, condensing, and forming stars, but the other half radiates its acquired gravitational energy at much lower temperatures, typically T<10^5 K, and the histogram of maximum gas temperatures is clearly bimodal. The "cold mode" of gas accretion dominates for low mass galaxies (M_baryon < 10^{10.3}Msun or M_halo < 10^{11.4}Msun), while the conventional "hot mode" dominates the growth of high mass systems. Cold accretion is often directed along filaments, allowing galaxies to efficiently draw gas from large distances, while hot accretion is quasi-spherical. The galaxy and halo mass dependence leads to redshift and environment dependence of cold and hot accretion rates, with cold mode dominating at high redshift and in low density regions today, and hot mode dominating in group and cluster environments at low redshift. Star formation rates closely track accretion rates, and we discuss the physics behind the observed environment and redshift dependence of galactic scale star formation. If we allowed hot accretion to be suppressed by conduction or AGN feedback, then the simulation predictions would change in interesting ways, perhaps resolving conflicts with the colors of ellipticals and the cutoff of the galaxy luminosity function. The transition between cold and hot accretion at M_h ~ 10^{11.4}Msun is similar to that found by Birnboim & Dekel (2003) using 1-d simulations and analytic arguments. The corresponding baryonic mass is tantalizingly close to the scale at which Kauffmann et al. (2003) find a marked shift in galaxy properties. We speculate on connections between these theoretical and observational transitions.Comment: 1 figure added, Appendix discussing SAMs added, some text changes. Matches the version accepted by MNRAS. 31 pages (MNRAS style), 21 figures,For high resolution version of the paper (highly recommended) follow http://www.astro.umass.edu/~keres/paper/ms2.ps.g

The formation of massive, quiescent galaxies at cosmic noon

The cosmic noon (z~1.5-3) marked a period of vigorous star formation for most galaxies. However, about a third of the more massive galaxies at those times were quiescent in the sense that their observed stellar populations are inconsistent with rapid star formation. The reduced star formation activity is often attributed to gaseous outflows driven by feedback from supermassive black holes, but the impact of black hole feedback on galaxies in the young Universe is not yet definitively established. We analyze the origin of quiescent galaxies with the help of ultra-high resolution, cosmological simulations that include feedback from stars but do not model the uncertain consequences of black hole feedback. We show that dark matter halos with specific accretion rates below ~0.25-0.4 per Gyr preferentially host galaxies with reduced star formation rates and red broad-band colors. The fraction of such halos in large dark matter only simulations matches the observed fraction of massive quiescent galaxies (~10^10-10^11 Msun). This strongly suggests that halo accretion rate is the key parameter determining which massive galaxies at z~1.5-3 become quiescent. Empirical models that connect galaxy and halo evolution, such as halo occupation distribution or abundance matching models, assume a tight link between galaxy properties and the masses of their parent halos. These models will benefit from adding the specific accretion rate of halos as a second model parameter.Comment: 5 pages, 5 figures, to appear in MNRAS Letter

Galactic r-process enrichment by neutron star mergers in cosmological simulations of a Milky Way-mass galaxy

We quantify the stellar abundances of neutron-rich r-process nuclei in cosmological zoom-in simulations of a Milky Way-mass galaxy from the Feedback In Realistic Environments project. The galaxy is enriched with r-process elements by binary neutron star (NS) mergers and with iron and other metals by supernovae. These calculations include key hydrodynamic mixing processes not present in standard semi-analytic chemical evolution models, such as galactic winds and hydrodynamic flows associated with structure formation. We explore a range of models for the rate and delay time of NS mergers, intended to roughly bracket the wide range of models consistent with current observational constraints. We show that NS mergers can produce [r-process/Fe] abundance ratios and scatter that appear reasonably consistent with observational constraints. At low metallicity, [Fe/H]<-2, we predict there is a wide range of stellar r-process abundance ratios, with both supersolar and subsolar abundances. Low-metallicity stars or stars that are outliers in their r-process abundance ratios are, on average, formed at high redshift and located at large galactocentric radius. Because NS mergers are rare, our results are not fully converged with respect to resolution, particularly at low metallicity. However, the uncertain rate and delay time distribution of NS mergers introduces an uncertainty in the r-process abundances comparable to that due to finite numerical resolution. Overall, our results are consistent with NS mergers being the source of most of the r-process nuclei in the Universe.Comment: Accepted for publication in MNRAS, 10 pages and 4 figures. Revised version: minor change

Reconciling dwarf galaxies with LCDM cosmology: Simulating a realistic population of satellites around a Milky Way-mass galaxy

Low-mass "dwarf" galaxies represent the most significant challenges to the cold dark matter (CDM) model of cosmological structure formation. Because these faint galaxies are (best) observed within the Local Group (LG) of the Milky Way (MW) and Andromeda (M31), understanding their formation in such an environment is critical. We present first results from the Latte Project: the Milky Way on FIRE (Feedback in Realistic Environments). This simulation models the formation of a MW-mass galaxy to z = 0 within LCDM cosmology, including dark matter, gas, and stars at unprecedented resolution: baryon particle mass of 7070 Msun with gas kernel/softening that adapts down to 1 pc (with a median of 25 - 60 pc at z = 0). Latte was simulated using the GIZMO code with a mesh-free method for accurate hydrodynamics and the FIRE-2 model for star formation and explicit feedback within a multi-phase interstellar medium. For the first time, Latte self-consistently resolves the spatial scales corresponding to half-light radii of dwarf galaxies that form around a MW-mass host down to Mstar > 10^5 Msun. Latte's population of dwarf galaxies agrees with the LG across a broad range of properties: (1) distributions of stellar masses and stellar velocity dispersions (dynamical masses), including their joint relation; (2) the mass-metallicity relation; and (3) a diverse range of star-formation histories, including their mass dependence. Thus, Latte produces a realistic population of dwarf galaxies at Mstar > 10^5 Msun that does not suffer from the "missing satellites" or "too big to fail" problems of small-scale structure formation. We conclude that baryonic physics can reconcile observed dwarf galaxies with standard LCDM cosmology.Comment: 7 pages, 5 figures. Accepted for publication in ApJ Letters. Several updates, including: (1) fixed a bug in halo finder, now identifies 13 satellite galaxies and more subhalos in the baryonic simulation; (2) fixed a minor bug in the feedback coupling and reran the simulation, resulting in a somewhat lower-mass host galaxy; (3) Fig 2 now shows stellar velocity dispersion profiles of satellite

The Origin and Evolution of the Galaxy Mass-Metallicity Relation

We use high-resolution cosmological zoom-in simulations from the Feedback in Realistic Environment (FIRE) project to study the galaxy mass-metallicity relations (MZR) from z=0-6. These simulations include explicit models of the multi-phase ISM, star formation, and stellar feedback. The simulations cover halo masses Mhalo=10^9-10^13 Msun and stellar mass Mstar=10^4-10^11 Msun at z=0 and have been shown to produce many observed galaxy properties from z=0-6. For the first time, our simulations agree reasonably well with the observed mass-metallicity relations at z=0-3 for a broad range of galaxy masses. We predict the evolution of the MZR from z=0-6 as log(Zgas/Zsun)=12+log(O/H)-9.0=0.35[log(Mstar/Msun)-10]+0.93 exp(-0.43 z)-1.05 and log(Zstar/Zsun)=[Fe/H]-0.2=0.40[log(Mstar/Msun)-10]+0.67 exp(-0.50 z)-1.04, for gas-phase and stellar metallicity, respectively. Our simulations suggest that the evolution of MZR is associated with the evolution of stellar/gas mass fractions at different redshifts, indicating the existence of a universal metallicity relation between stellar mass, gas mass, and metallicities. In our simulations, galaxies above Mstar=10^6 Msun are able to retain a large fraction of their metals inside the halo, because metal-rich winds fail to escape completely and are recycled into the galaxy. This resolves a long-standing discrepancy between "sub-grid" wind models (and semi-analytic models) and observations, where common sub-grid models cannot simultaneously reproduce the MZR and the stellar mass functions.Comment: 17 pages, 14 figures, re-submitted to MNRAS after revisions on referee comment

Accretion, feedback and galaxy bimodality: a comparison of the GalICS semi-analytic model and cosmological SPH simulations

We compare the galaxy population of an SPH simulation to those predicted by the GalICS semi-analytic model and a stripped down version without supernova and AGN feedback. The SPH simulation and the no-feedback GalICS model make similar predictions for the baryonic mass functions of galaxies and for the dependence of these mass functions on environment and redshift. The two methods also make similar predictions for the galaxy content of dark matter haloes as a function of halo mass and for the gas accretion history of galaxies. Both the SPH and no-feedback GalICS models predict a bimodal galaxy population at z=0. The "red'' sequence of gas poor, old galaxies is populated mainly by satellite systems while, contrary to observations, the central galaxies of massive haloes lie on the "blue'' star-forming sequence as a result of continuing hot gas accretion at late times. Furthermore, both models overpredict the observed baryonic mass function, especially at the high mass end. In the full GalICS model, supernova-driven outflows reduce the masses of low and intermediate mass galaxies by about a factor of two. AGN feedback suppresses gas cooling in large haloes, producing a sharp cut-off in the baryonic mass function and moving the central galaxies of these massive haloes to the red sequence. Our results imply that the observational failings of the SPH simulation and the no-feedback GalICS model are a consequence of missing input physics rather than computational inaccuracies, that truncating gas accretion by satellite galaxies automatically produces a bimodal galaxy distribution with a red sequence, but that explaining the red colours of the most massive galaxies requires a mechanism like AGN feedback that suppresses the accretion onto central galaxies in large haloes.Comment: 17 pages, 11 figures, submitted to MNRA