Andrew Wetzel

MERGERS AND BULGE FORMATION IN ΛCDM: WHICH MERGERS MATTER?

We use a suite of semi-empirical models to predict the galaxy-galaxy merger rate and relative contributions to bulge growth as a function of mass (both halo and stellar), redshift, and mass ratio. The models use empirical constraints on the halo occupation distribution, evolved forward in time, to robustly identify where and when galaxy mergers occur. Together with the results of high-resolution merger simulations, this allows us to quantify the relative contributions of mergers with different properties (e.g., mass ratios, gas fractions, redshifts) to the bulge population. We compare with observational constraints, and find good agreement. We also provide useful fitting functions and make public a code to reproduce the predicted merger rates and contributions to bulge mass growth. We identify several robust conclusions. (1) Major mergers dominate the formation and assembly of ~L * bulges and the total spheroid mass density, but minor mergers contribute a non-negligible ~30%. (2) This is mass dependent: bulge formation and assembly is dominated by more minor mergers in lower-mass systems. In higher-mass systems, most bulges originally form in major mergers near ~L *, but assemble in increasingly minor mergers. (3) The minor/major contribution is also morphology dependent: higher B/T systems preferentially form in more major mergers, with B/T roughly tracing the mass ratio of the largest recent merger; lower B/T systems preferentially form in situ from minor mergers. (4) Low-mass galaxies, being gas-rich, require more mergers to reach the same B/T as high-mass systems. Gas-richness dramatically suppresses the absolute efficiency of bulge formation, but does not strongly influence the relative contribution of major versus minor mergers. (5) Absolute merger rates at fixed mass ratio increase with galaxy mass. (6) Predicted merger rates agree well with those observed in pair and morphology-selected samples, but there is evidence that some morphology-selected samples include contamination from minor mergers. (7) Predicted rates also agree with the integrated growth in bulge mass density with cosmic time, but with a factor ~2 uncertainty in both—up to half the bulge mass density could come from non-merger processes. We systematically vary the model assumptions, totaling ~103 model permutations, and quantify the resulting uncertainties. Our conclusions regarding the importance of different mergers for bulge formation are very robust to these changes. The absolute predicted merger rates are systematically uncertain at the factor ~2 level; uncertainties grow at the lowest masses and high redshifts.

What determines satellite galaxy disruption

In hierarchical structure formation, dark matter halos that merge with larger halos can persist as subhalos. These subhalos are likely hosts of visible galaxies. While the dense halo environment rapidly strips subhalos of their dark mass, the compact luminous material can remain intact for some time, making the correspondence of galaxies with severely stripped subhalos unclear. Many galaxy evolution models assume that satellite galaxies eventually merge with their central galaxy, but this ignores the possibility of satellite tidal disruption. We use a high-resolution N-body simulation of cosmological volume to explore satellite galaxy merging/disruption criteria based on dark matter subhalo dynamics. We explore the impact that satellite merging/disruption has on the Halo Occupation Distribution and radial profile of the remnants. Using abundance matching to assign stellar mass/luminosity to subhalos, we compare with observed galaxy clustering, satellite fractions, cluster satellite luminosity functions, finding that subhalos reproduce well these observables. Our results imply that satellite subhalos

On the orbits of infalling satellite haloes

The orbital properties of infalling satellite haloes set the initial conditions which control the subsequent evolution of subhaloes and the galaxies that they host, with implications for mass stripping, star formation quenching and merging. Using a high-resolution cosmological N-body simulation, we examine the orbital parameters of satellite haloes as they merge with larger host haloes, focusing primarily on orbital circularity and pericentre. We explore in detail how these orbital parameters depend on mass and redshift. Satellite orbits become more radial and plunge deeper into their host halo at higher host halo mass, but they do not significantly depend on satellite halo mass. Additionally, satellite orbits become more radial and plunge deeper into their host haloes at higher redshift. We also examine satellite velocities, finding that most satellites infall with less specific angular momentum than the host halo virial value, but that satellites are ‘hotter’ than the host virial velocity. We discuss the implications of these results to the processes of galaxy formation and evolution, and we provide fitting formulae to the mass and redshift dependence of satellite orbital circularity and pericentre.

What does Clustering Tell us About the Buildup of the Red Sequence

We analyze the clustering of red and blue galaxies from four samples spanning a redshift range of 0.4 < z < 2.0 to test the various scenarios by which galaxies evolve onto the red sequence. The data are taken from the UKIDSS Ultra Deep Survey, DEEP2, and COMBO-17. The use of clustering allows us to determine what fraction of the red sequence is made up of central galaxies and satellite galaxies. At all redshifts, including z = 0, the data are consistent with ~60% of satellite galaxies being red or quenched, implying that ~1/3 of the red sequence is comprised of satellite galaxies. More than three-fourths of red satellite galaxies were moved to the red sequence after they were accreted onto a larger halo. The constant fraction of satellite galaxies that are red yields a quenching time for satellite galaxies that depends on redshift in the same way as halo dynamical times: t Q ~ (1 + z)–1.5. In three of the four samples, the data favor a model in which red central galaxies are a random sample of all central galaxies; there is no preferred halo mass scale at which galaxies make the transition from star-forming to red and dead. The large errors on the fourth sample inhibit any conclusions. Theoretical models in which star formation is quenched above a critical halo mass are excluded by these data. A scenario in which mergers create red central galaxies imparts a weaker correlation between halo mass and central galaxy color, but even the merger scenario creates tension with red galaxy clustering at redshifts above 0.5. These results suggest that the mechanism by which central galaxies become red evolves from z = 0.5 to z = 0.

The Clustering of Massive Halos

The clustering properties of dark matter halos are a firm prediction of modern theories of structure formation. We use two large-volume, high-resolution N-body simulations to study how the correlation function of massive dark matter halos depends on their mass, assembly history, and recent merger activity. We find that halos with the lowest concentrations are currently more clustered than those of higher concentration, the size of the effect increasing with halo mass; this agrees with trends found in studies of lower mass halos. The clustering dependence on other characterizations of the full mass accretion history appears weaker than the effect with concentration. Using the integrated correlation function, marked correlation functions, and a power-law fit to the correlation function, we find evidence that halos that have recently undergone a major merger or a large mass gain have slightly enhanced clustering relative to a randomly chosen population with the same mass distribution.

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: the low redshift sample

We report on the small scale (0:5 < r < 40h 1 Mpc) clustering of 78895 massive (M 10 11:3 M ) galaxies at 0:2 < z < 0:4 from the first two years of data from the Baryon Oscillation Spectroscopic Survey (BOSS), to be released as part of SDSS Data Release 9 (DR9). We describe the sample selection, basic properties of the galaxies, and caveats for working with the data. We calculate the real- and redshift-space two-point correlation functions of these galaxies, fit these measurements using Halo Occupation Distribution (HOD) modeling within dark matter cosmological simulations, and estimate the errors using mock catalogs. These galaxies lie in massive halos, with a mean halo mass of 5:2 10 13 h 1 M , a large scale bias of 2:0, and a satellite fraction of 12 2%. Thus,

Simulating subhaloes at high redshift: merger rates, counts and types

Galaxies are believed to be in one-to-one correspondence with simulated dark matter subhaloes. We use high-resolutionN-body simulations of cosmological volumes to calculate the statistical properties of subhalo (galaxy) major mergers at high redshift (z = 0.6–5). We measure the evolution of the galaxy merger rate, finding that it is much shallower than the merger rate of dark matter host haloes at z> 2.5, but roughly parallels that of haloes at z< 1.6. We also track the detailed merger histories of individual galaxies and measure the likelihood of multiple mergers per halo or subhalo. We examine satellite merger statistics in detail: 15–35 per cent of all recently merged galaxies are satellites, and satellites are twice as likely as centrals to have had a recent major merger. Finally, we show how the differing evolution of the merger rates of haloes and galaxies leads to the evolution of the average satellite occupation per halo, noting that for a fixed halo mass, the satellite halo occupation peaks at z ∼ 2.5.

FIRE-2 Simulations: Physics versus Numerics in Galaxy Formation

The Feedback In Realistic Environments (FIRE) project explores feedback in cosmological galaxy formation simulations. Previous FIRE simulations used an identical source code (“FIRE-1”) for consistency. Motivated by the development of more accurate numerics – including hydrodynamic solvers, gravitational softening, and supernova coupling algorithms – and exploration of new physics (e.g. magnetic fields), we introduce “FIRE-2”, an updated numerical implementation of FIRE physics for the GIZMO code. We run a suite of simulations and compare against FIRE-1: overall, FIRE-2 improvements do not qualitatively change galaxy-scale properties. We pursue an extensive study of numerics versus physics. Details of the star-formation algorithm, cooling physics, and chemistry have weak effects, provided that we include metal-line cooling and star formation occurs at higher-than-mean densities. We present new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are robust to numerics we test, provided: (1) Toomre masses are resolved; (2) feedback coupling ensures conservation, and (3) individual supernovae are time-resolved. Stellar masses and profiles are most robust to resolution, followed by metal abundances and morphologies, followed by properties of winds and circum-galactic media (CGM). Central (∼kpc) mass concentrations in massive (>L*) galaxies are sensitive to numerics (via trapping/recycling of winds in hot halos). Multiple feedback mechanisms play key roles: supernovae regulate stellar masses/winds; stellar mass-loss fuels late star formation; radiative feedback suppresses accretion onto dwarfs and instantaneous star formation in disks. We provide all initial conditions and numerical algorithms used.

SATELLITE DWARF GALAXIES IN A HIERARCHICAL UNIVERSE: THE PREVALENCE OF DWARF-DWARF MAJOR MERGERS

Mergers are a common phenomenon in hierarchical structure formation, especially for massive galaxies and clusters, but their importance for dwarf galaxies in the Local Group remains poorly understood. We investigate the frequency of major mergers between dwarf galaxies in the Local Group using the ELVIS suite of cosmological zoom-in dissipationless simulations of Milky Way- and M31-like host halos. We find that ~10% of satellite dwarf galaxies with Mstar > 10^6 M_☉ that are within the host virial radius experienced a major merger of stellar mass ratio closer than 0.1 since z = 1, with a lower fraction for lower mass dwarf galaxies. Recent merger remnants are biased toward larger radial distance and more recent virial infall times, because most recent mergers occurred shortly before crossing within the virial radius of the host halo. Satellite–satellite mergers also occur within the host halo after virial infall, catalyzed by the large fraction of dwarf galaxies that fell in as part of a group. The merger fraction doubles for dwarf galaxies outside of the host virial radius, so the most distant dwarf galaxies in the Local Group are the most likely to have experienced a recent major merger. We discuss the implications of these results on observable dwarf merger remnants, their star formation histories, the gas content of mergers, and massive black holes in dwarf galaxies.

Not so lumpy after all: modelling the depletion of dark matter subhaloes by Milky Way-like galaxies

Among the most important goals in cosmology is detecting and quantifying small (M_(halo)≃10^(6−9) M⊙) dark matter (DM) subhaloes. Current probes around the Milky Way (MW) are most sensitive to such substructure within ∼20 kpc of the halo centre, where the galaxy contributes significantly to the potential. We explore the effects of baryons on subhalo populations in ΛCDM using cosmological zoom-in baryonic simulations of MW-mass haloes from the Latte simulation suite, part of the Feedback In Realistic Environments (FIRE) project. Specifically, we compare simulations of the same two haloes run using (1) DM-only (DMO), (2) full baryonic physics and (3) DM with an embedded disc potential grown to match the FIRE simulation. Relative to baryonic simulations, DMO simulations contain ∼2 × as many subhaloes within 100 kpc of the halo centre; this excess is ≳5 × within 25 kpc. At z = 0, the baryonic simulations are completely devoid of subhaloes down to 3×10^6M⊙ within 15 kpc of the MW-mass galaxy, and fewer than 20 surviving subhaloes have orbital pericentres <20 kpc. Despite the complexities of baryonic physics, the simple addition of an embedded central disc potential to DMO simulations reproduces this subhalo depletion, including trends with radius, remarkably well. Thus, the additional tidal field from the central galaxy is the primary cause of subhalo depletion. Subhaloes on radial orbits that pass close to the central galaxy are preferentially destroyed, causing the surviving population to have tangentially biased orbits compared to DMO predictions. Our method of embedding a potential in DMO simulations provides a fast and accurate alternative to full baryonic simulations, thus enabling suites of cosmological simulations that can provide accurate and statistical predictions of substructure populations.