Push it to the limit: Local Group constraints on high-redshift stellar mass functions for Mstar > 10^5 Msun

We constrain the evolution of the galaxy stellar mass function from 2 < z < 5 for galaxies with stellar masses as low as 10^5 Msun by combining star formation histories of Milky Way satellite galaxies derived from deep Hubble Space Telescope observations with merger trees from the ELVIS suite of N-body simulations. This approach extends our understanding more than two orders of magnitude lower in stellar mass than is currently possible by direct imaging. We find the faint end slopes of the mass functions to be alpha= -1.42(+0.07/-0.05) at z = 2 and alpha = -1.57^(+0.06/-0.06) at z = 5, and show the slope only weakly evolves from z = 5 to z = 0. Our findings are in stark contrast to a number of direct detection studies that suggest slopes as steep as alpha = -1.9 at these epochs. Such a steep slope would result in an order of magnitude too many luminous Milky Way satellites in a mass regime that is observationally complete (Mstar > 2*10^5 Msun 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 LCDM numerical simulations

Through a Smoother Lens: An expected absence of LCDM substructure detections from hydrodynamic and dark matter only simulations

A fundamental prediction of the cold dark matter cosmology is the existence of a large number of dark subhalos around galaxies, most of which should be entirely devoid of stars. Confirming the existence of dark substructures stands among the most important empirical challenges in modern cosmology: if they are found and quantified with the mass spectrum expected, then this would close the door on a vast array of competing theories. But in order for observational programs of this kind to reach fruition, we need robust predictions. Here we explore substructure predictions for lensing using galaxy lens-like hosts at z=0.2 from the Illustris simulations both in full hydrodynamics and dark matter only. We quantify substructures more massive than ~ 10^9 M_sun, comparable to current lensing detections derived from HST, Keck, and ALMA. The addition of full hydrodynamics reduces the overall subhalo mass function by about a factor of two. Even for the dark matter only runs, most (~ 85%) lines of sight through projected cylinders of size close to an Einstein radius contain no substructures larger than 10^9 M_sun. The fraction of empty sight lines rises to ~ 95% in full physics simulations. This suggests we will likely need hundreds of strong lensing systems suitable for substructure studies, as well as predictions that include the effects of baryon physics on substructure, to properly constrain cosmological models. Fortunately, the field is poised to fulfill these requirements.Comment: 11 pages, 9 figure

A Predicted Correlation Between Age Gradient and Star Formation History in FIRE Dwarf Galaxies

We explore the radial variation of star formation histories in dwarf galaxies simulated with Feedback In Realistic Environments (FIRE) physics. The sample contains 9 low-mass field dwarfs with M_ star = 10^5 - 10^7 M_sun from previous FIRE results, and a new suite of 17 higher mass field dwarfs with M_star = 10^7 - 10^9 M_sun introduced here. We find that 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 star formation history of the merging galaxy. In galaxies without significant mergers, early feedback pushes stars to the outskirts at early times. Interestingly, among galaxies without mergers, those with large dark matter cores have flatter age gradients because these galaxies have more late-time feedback. 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 star formation history. We show that central fields can be biased young by a few 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 star formation history.Comment: 13 pages, 8 figures. Submitted to MNRAS, comments welcom

A Testable Conspiracy: Simulating Baryonic Effects on Self-interacting Dark Matter Halos

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 σ/m =0.1 cm^2 g^(-1) 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 σ/m = 0.5 cm^2 g^(-1) is not dense enough. Our findings are in agreement with previous results that σ/m\ ≳ 0.1 cm^2 g^(1-) 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

Is it possible to reconcile extragalactic IMF variations with a universal Milky Way IMF?

One of the most robust observations of the stellar initial mass function (IMF) is its near-universality in the Milky Way and neighbouring galaxies. But recent observations of early-type galaxies can be interpreted to imply a ‘bottom-heavy’ IMF, while others of ultrafaint dwarfs could imply a ‘top-heavy’ IMF. This would impose powerful constraints on star formation models. We explore what sort of ‘cloud-scale’ IMF models could possibly satisfy these constraints. We utilize simulated galaxies that reproduce (broadly) the observed galaxy properties, while they also provide the detailed star formation history and properties of each progenitor star-forming cloud. We then consider generic models where the characteristic mass of the IMF is some arbitrary power-law function of progenitor cloud properties, along with well-known literature IMF models which scale with Jeans mass, ‘turbulent Bonnor–Ebert mass’, temperature, the opacity limit, metallicity, or the ‘protostellar heating mass’. We show that no IMF models currently in the literature – nor any model where the turnover mass is an arbitrary power-law function of a combination of cloud temperature/density/size/metallicity/velocity dispersion/magnetic field – can reproduce the claimed IMF variation in ellipticals or dwarfs without severely violating observational constraints in the Milky Way. Specifically, they predict too much variation in the ‘extreme’ environments of the Galaxy compared to that observed. Either the IMF varies in a more complicated manner, or alternative interpretations of the extragalactic observations must be explored

Star formation histories of dwarf galaxies in the FIRE simulations: dependence on mass and Local Group environment

We study star formation histories (SFHs) of $\simeq500$ dwarf galaxies (stellar mass $M_\ast = 10^5 - 10^9\,M_\odot$) from FIRE-2 cosmological zoom-in simulations. We compare dwarfs around individual Milky Way (MW)-mass galaxies, dwarfs in Local Group (LG)-like environments, and true field (i.e. isolated) dwarf galaxies. We reproduce observed trends wherein higher-mass dwarfs quench later (if at all), regardless of environment. We also identify differences between the environments, both in terms of "satellite vs. central" and "LG vs. individual MWvs. isolated dwarf central." Around the individual MW-mass hosts, we recover the result expected from environmental quenching: central galaxies in the "near field" have more extended SFHs than their satellite counterparts, with the former more closely resemble isolated ("true field") dwarfs (though near-field centrals are still somewhat earlier forming). However, this difference is muted in the LG-like environments, where both near-field centrals and satellites have similar SFHs, which resemble satellites of single MW-mass hosts. This distinction is strongest for $M_\ast = 10^6 - 10^7\,M_\odot$ but exists at other masses. Our results suggest that the paired halo nature of the LG may regulate star formation in dwarf galaxies even beyond the virial radii of the MW and Andromeda. Caution is needed when comparing zoom-in simulations targeting isolated dwarf galaxies against observed dwarf galaxies in the LG.Comment: Main text: 11 pages, 8 figures; appendices: 4 pages, 4 figures. Submitted to MNRAS; comments welcom