We constrain the evolution of the galaxy stellar mass function from 2 2 × 105M⊙ at z = 0). The most recent studies from ZFOURGE and CANDELS also suggest flatter faint end slopes that are consistent with our results, but with a lower degree of precision. This work illustrates the strong connections between low and high-z observations when viewed through the lens of lambda cold dark matter numerical simulations.

We investigate the response of self-interacting dark matter (SIDM) halos to the growth of galaxy potentials using idealized simulations, with each run in tandem with collisionless cold dark matter (CDM). We find that if the stellar potential strongly dominates in the central parts of a galaxy, then SIDM halos can be as dense as CDM halos on observable scales. For extreme cases, core collapse can occur, leading to SIDM halos that are denser and cuspier than their CDM counterparts. If the stellar potential is not dominant, then SIDM halos retain isothermal cores with densities far below CDM predictions. When a disk is present, the inner SIDM halo becomes more flattened in the disk plane than the CDM halo. These results are in excellent quantitative agreement with the predictions of Kaplinghat et al. We also simulated a cluster halo with a central stellar distribution similar to the brightest central galaxy of the cluster A2667. An SIDM halo simulated with the cross-section over mass provides a good match to the measured dark matter (DM) density profile, while an adiabatically contracted CDM halo is denser and cuspier. The profile of the same halo simulated with is not dense enough. Our findings are in agreement with previous results that is disfavored for DM collision velocities above about 1500 km s-1. More generally, the interaction between baryonic potentials and SIDM densities offers new directions for constraining SIDM cross-sections in galaxies where baryons are dynamically important.

We analyse the cold dark matter density profiles of 54 galaxy haloes simulated with Feedback In Realistic Environments (FIRE)-2 galaxy formation physics, each resolved within 0.5 per cent of the halo virial radius. These haloes contain galaxies with masses that range from ultrafaint dwarfs (M* ~ 104.5M) to the largest spirals (M* ˜ 1011M) and have density profiles that are both cored and cuspy. We characterize our results using a new, analytic density profile that extends the standard two-parameter Einasto form to allow for a pronounced constant density core in the resolved innermost radius. With one additional core-radius parameter, rc, this three-parameter core-Einasto profile is able to characterize our feedback-impacted dark matter haloes more accurately than other three-parameter profiles proposed in the literature. To enable comparisons with observations, we provide fitting functions for rc and other profile parameters as a function of both M* and M/Mhalo. In agreement with past studies, we find that dark matter core formation is most efficient at the characteristic stellar-to-halo mass ratio M/Mhalo ˜5 × 10-3, or M* ˜ 109 M, with cores that are roughly the size of the galaxy half-light radius, rc ˜ 15kpc. Furthermore, we find no evidence for core formation at radii > 100 pc in galaxies with M*/Mhalo < 5 × 10-4 or M*-106M. For Milky Way-size galaxies, baryonic contraction often makes haloes significantly more concentrated and dense at the stellar half-light radius than DMO runs. However, even at the Milky Way scale, FIRE-2 galaxy formation still produces small dark matter cores of ˜ 0.5-2 kpc in size. Recent evidence for a ~2 kpc core in the Milky Way's dark matter halo is consistent with this expectation.

We study stellar-halo formation using six Milky-Way-mass galaxies in FIRE-2 cosmological zoom simulations. We find that 5-40 per cent of the outer (50-300 kpc) stellar halo in each system consists of in-situ stars that were born in outflows from the main galaxy. Outflow stars originate from gas accelerated by superbubble winds, which can be compressed, cool, and form co-moving stars. The majority of these stars remain bound to the halo and fall back with orbital properties similar to the rest of the stellar halo at z = 0. In the outer halo, outflow stars are more spatially homogeneous, metal-rich, and alpha-element-enhanced than the accreted stellar halo. At the solar location, up to ∼10 per cent of our kinematically identified halo stars were born in outflows; the fraction rises to as high as ∼40 per cent for the most metal-rich local halo stars ([Fe/H] >-0.5). Such stars can be retrograde and create features similar to the recently discoveredMilkyWay 'Splash' in phase space.We conclude that theMilkyWay stellar halo could contain local counterparts to stars that are observed to form in molecular outflows in distant galaxies. Searches for such a population may provide a new, near-field approach to constraining feedback and outflow physics. A stellar halo contribution from outflows is a phase-reversal of the classic halo formation scenario of Eggen, Lynden-Bell & Sandange, who suggested that halo stars formed in rapidly infalling gas clouds. Stellar outflows may be observable in direct imaging of external galaxies and could provide a source for metal-rich, extreme-velocity stars in the Milky Way.

Among the most important goals in cosmology is detecting and quantifying small (Mhalo ≃ 106-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 × 106M⊙ 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.

We explore the radial variation of star formation histories (SFHs) in dwarf galaxies simulated with Feedback In Realistic Environments (FIRE) physics. The sample contains 26 field dwarf galaxies with Mstar = 105–109 M. We find age gradients are common in our dwarfs, with older stars dominant at large radii. The strength of the gradient correlates with overall galaxy age such that earlier star formation produces a more pronounced gradient. The relation between formation time and strength of the gradient is driven by both mergers and star formation feedback. Mergers can both steepen and flatten the age gradient depending on the timing of the merger and SFHs of the merging galaxy. In galaxies without significant mergers, feedback pushes stars to the outskirts. The strength of the age gradient is determined by the subsequent evolution of the galaxy. Galaxies with weak age gradients constantly grow to z = 0, meaning that young star formation occurs at a similar radius to which older stars are heated to. In contrast, galaxies with strong age gradients tend to maintain a constant half-mass radius over time. If real galaxies have age gradients as we predict, stellar population studies that rely on sampling a limited fraction of a galaxy can give a biased view of its global SFH. Central fields can be biased young by Gyrs while outer fields are biased old. Fields positioned near the 2D half-light radius will provide the least biased measure of a dwarf galaxy’s global SFH.

We explore the origin of stellar metallicity gradients in simulated and observed dwarf galaxies. We use FIRE-2 cosmological baryonic zoom-in simulations of 26 isolated galaxies as well as existing observational data for 10 Local Group dwarf galaxies. Our simulated galaxies have stellar masses between 105.5 and 108.6 MO. Whilst gas-phase metallicty gradients are generally weak in our simulated galaxies, we find that stellar metallicity gradients are common, with central regions tending to be more metal-rich than the outer parts. The strength of the gradient is correlated with galaxy-wide median stellar age, such that galaxies with younger stellar populations have flatter gradients. Stellar metallicty gradients are set by two competing processes: (1) the steady 'puffing' of old, metal-poor stars by feedback-driven potential fluctuations and (2) the accretion of extended, metal-rich gas at late times, which fuels late-time metal-rich star formation. If recent star formation dominates, then extended, metal-rich star formation washes out pre-existing gradients from the 'puffing' process. We use published results from ten Local Group dwarf galaxies to show that a similar relationship between age and stellar metallicity-gradient strength exists among real dwarfs. This suggests that observed stellar metallicity gradients may be driven largely by the baryon/feedback cycle rather than by external environmental effects.