Within lambda cold dark matter (CDM), dwarf galaxies like the Large Magellanic Cloud (LMC) are expected to host numerous dark matter subhaloes, several of which should host faint dwarf companions. Recent Gaia proper motions confirm new members of the LMC system in addition to the previously known SMC, including two classical dwarf galaxies (M∗ > 105 M; Carina and Fornax) as well as several ultrafaint dwarfs (Car2, Car3, Hor1, and Hyd1). We use the Feedback In Realistic Environments (FIRE) simulations to study the dark and luminous (down to ultrafaint masses, M∗ ∼6×103 M) substructure population of isolated LMC-mass hosts (M200m = 1–3×1011 M) and place the Gaia + DES results in a cosmological context. By comparing number counts of subhaloes in simulations with and without baryons, we find that, within 0.2 r200m, LMC-mass hosts deplete ∼30 per cent of their substructure, significantly lower than the ∼70 per cent of substructure depleted by Milky Way (MW) mass hosts. For our highest resolution runs (mbary = 880 M), ∼ 5–10 subhaloes form galaxies with M∗ ≥ 104 M , in agreement with the seven observationally inferred pre-infall LMC companions. However, we find steeper simulated luminosity functions than observed, hinting at observation incompleteness at the faint end. The predicted DM content for classical satellites in FIRE agrees with observed estimates for Carina and Fornax, supporting the case for an LMC association. We predict that tidal stripping within the LMC potential lowers the inner dark matter density of ultrafaint companions of the LMC. Thus, in addition to their orbital consistency, the low densities of dwarfs Car2, Hyd1, and Hyd2 reinforce their likelihood of Magellanic association.

Understanding the rate at which stars form is central to studies of galaxy formation. Observationally, the star formation rates (SFRs) of galaxies are measured using the luminosity in different frequency bands, often under the assumption of a time-steady SFR in the recent past. We use star formation histories (SFHs) extracted from cosmological simulations of star-forming galaxies from the FIRE project to analyse the time-scales to which the H α and far-ultraviolet (FUV) continuum SFR indicators are sensitive. In these simulations, the SFRs are highly time variable for all galaxies at high redshift, and continue to be bursty to z = 0 in dwarf galaxies. When FIRE SFHs are partitioned into their bursty and time-steady phases, the best-fitting FUV time-scale fluctuates from its ∼10 Myr value when the SFR is time-steady to ≥100 Myr immediately following particularly extreme bursts of star formation during the bursty phase. On the other hand, the best-fitting averaging time-scale for H α is generally insensitive to the SFR variability in the FIRE simulations and remains ∼5 Myr at all times. These time-scales are shorter than the 100 and 10 Myr time-scales sometimes assumed in the literature for FUV and H α, respectively, because while the FUV emission persists for stellar populations older than 100 Myr, the time-dependent luminosities are strongly dominated by younger stars. Our results confirm that the ratio of SFRs inferred using H α versus FUV can be used to probe the burstiness of star formation in galaxies.

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.