How low does it go? Too few Galactic satellites with standard reionization quenching
ABSTRACT A standard prediction of galaxy formation theory is that the ionizing background suppresses galaxy formation in haloes with peak circular velocities smaller than $V_{\rm peak}\simeq 20 \, \rm km \, s^{-1}$, rendering the majority of haloes below this scale completely dark. We use a suite of cosmological zoom simulations of Milky Way-like haloes that include central Milky Way disc galaxy potentials to investigate the relationship between subhaloes and ultrafaint galaxies. We find that there are far too few subhaloes within 50 kpc of the Milky Way that had $V_{\rm peak}\gtrsim 20\, \rm km \, s^{-1}$ to account for the number of ultrafaint galaxies already known within that volume today. In order to match the observed count, we must populate subhaloes down to $V_{\rm peak}\simeq 6\, \rm km \, s^{-1}$ with ultrafaint dwarfs. The required haloes have peak virial temperatures as low as 1500 K, well below the atomic hydrogen cooling limit of 104 K. Allowing for the possibility that the Large Magellanic Cloud contributes several of the satellites within 50 kpc could potentially raise this threshold to $10\, \rm km \, s^{-1}$ (4000 K), still below the atomic cooling limit and far below the nominal reionization threshold.
Authors:
; ; ; ; ;
Award ID(s):
Publication Date:
NSF-PAR ID:
10174837
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
488
Issue:
4
Page Range or eLocation-ID:
4585 to 4595
ISSN:
0035-8711
1. ABSTRACT 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{{\ \rm per\ cent}}$ of the halo virial radius. These haloes contain galaxies with masses that range from ultrafaint dwarfs ($M_\star \simeq 10^{4.5}\, \mathrm{M}_{\odot }$) to the largest spirals ($M_\star \simeq 10^{11}\, \mathrm{M}_{\odot }$) 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_{\star } \sim 10^9 \, \mathrm{M}_{\odot }$, with cores that are roughly the size of the galaxy half-light radius, rc ≃ 1−5 kpc. Furthermore,more »
2. ABSTRACT We revisit the question of ‘hot mode’ versus ‘cold mode’ accretion on to galaxies using steady-state cooling flow solutions and idealized 3D hydrodynamic simulations. We demonstrate that for the hot accretion mode to exist, the cooling time is required to be longer than the free-fall time near the radius where the gas is rotationally supported, Rcirc, i.e. the existence of the hot mode depends on physical conditions at the galaxy scale rather than on physical conditions at the halo scale. When allowing for the depletion of the halo baryon fraction relative to the cosmic mean, the longer cooling times imply that a virialized gaseous halo may form in halo masses below the threshold of $\sim 10^{12}\, {\rm M_{\odot }}$ derived for baryon-complete haloes. We show that for any halo mass there is a maximum accretion rate for which the gas is virialized throughout the halo and can accrete via the hot mode of ${\dot{M}}_{\rm crit}\approx 0.7(v_{\rm c}/100\, \rm km\ s^{-1})^{5.4}(R_{\rm circ}/10\, {\rm kpc})(Z/\, {\rm Z_{\odot }})^{-0.9}\, {\rm M_{\odot }}\, {\rm yr}^{-1}$, where Z and vc are the metallicity and circular velocity measured at Rcirc. For accretion rates $\gtrsim {\dot{M}}_{\rm crit}$ the volume-filling gas phase can in principle be ‘transonic’ –more »
3. ABSTRACT A possibility of DM being multicomponent has a strong implication on resolving decades-long known cosmological problems on small scale. In addition to elastic scattering, the model allows for inelastic interactions, which can be characterized by a ‘velocity kick’ parameter. The simplest 2cDM model with cross-section $0.01\lesssim \sigma /m\lt 1\, \textrm {cm}^{2}{ \rm g}^{-1}$ and the kick velocity $V_{\mathrm{ k}}\simeq 100\, \rm {km\, s}^{-1}$ have been shown to robustly resolve the missing satellites, core-cusp, and too-big-to-fail problems in N-body cosmological simulations tested on Milky Way (MW)-like haloes of a virial mass ${\sim}5 \times 10^{11}\, {\rm M_{\odot }}$ (Papers I & II). With the aim of further constraining the parameter space available for the 2cDM model, we extend our analysis to dwarf and galaxy cluster haloes with their virial mass of ∼107−108 and ${\sim}10^{13} - 10^{14}\, {\rm M_{\odot }}$, respectively. We find that σ0/m ≳ 0.1 cm2g−1 is preferentially disfavoured for both dwarfs and galaxy cluster haloes in comparison with observations, while σ0/m = 0.001 cm2g−1 causes little perceptible difference from that of the CDM counterpart for most of the cross-section’s velocity dependence studied in this work. Our main result is that within the reasonable set of parameters, the 2cDM model can successfully explain themore »
We perform high-resolution simulations of an MW-like galaxy in a self-interacting cold dark matter model with elastic cross-section over mass of $1~\rm cm^2\, g^{-1}$ (SIDM) and compare to a model without self-interactions (CDM). We run our simulations with and without a time-dependent embedded potential to capture effects of the baryonic disc and bulge contributions. The CDM and SIDM simulations with the embedded baryonic potential exhibit remarkably similar host halo profiles, subhalo abundances, and radial distributions within the virial radius. The SIDM host halo is denser in the centre than the CDM host and has no discernible core, in sharp contrast to the case without the baryonic potential (core size ${\sim}7 \, \rm kpc$). The most massive subhaloes (with $V_{\mathrm{peak}}\gt 20 \, \rm km\, s^{-1}$) in our SIDM simulations, expected to host the classical satellite galaxies, have density profiles that are less dense than their CDM analogues at radii less than 500 pc but the deviation diminishes for less massive subhaloes. With the baryonic potential included in the CDM and SIDM simulations, the most massive subhaloes do not display the too-big-to-fail problem. However, the least dense among the massive subhaloes in both these simulations tend to have the smallest pericenter values,more »