Search

Massive dwarf galaxies that merge with the Milky Way on prograde orbits can be dragged into the disk plane before being completely disrupted. Such mergers can contribute to an accreted stellar disk and a dark matter disk. Here we present Nyx, a vast stellar stream in the vicinity of the Sun, which provides the first indication that such an event occurred in the Milky Way. We identify about 200 stars that have coherent radial and prograde motion in this stream using a catalogue of accreted stars built by applying deep learning methods to the Gaia data. Taken together with chemical abundance and orbital information, these results strongly favour the interpretation that Nyx is the remnant of a disrupted dwarf galaxy. Further justified by FIRE hydrodynamic simulations, we demonstrate that prograde streams like Nyx can be found in the disk plane of galaxies and identified using our methods.

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.

With Gaia Data Release 2, the astronomical community is entering a new era of multidimensional surveys of the Milky Way. This new phase-space view of our Galaxy demands new tools for comparing observations to simulations of Milky Way-mass galaxies in a cosmological context, to test the physics of both dark matter and galaxy formation. We present ananke, a framework for generating synthetic phase-space surveys from high-resolution baryonic simulations, and use it to generate a suite of synthetic surveys resembling Gaia DR2 in data structure, magnitude limits, and observational errors. We use three cosmological simulations of Milky Way-mass galaxies from the Latte suite of the Feedback In Realistic Environments project, which feature self-consistent clustering of star formation in dense molecular clouds and thin stellar/gaseous disks in live cosmological halos with satellite dwarf galaxies and stellar halos. We select three solar viewpoints from each simulation to generate nine synthetic Gaia-like surveys. We sample synthetic stars by assuming each star particle (of mass 7070 M o˙) represents a single stellar population. At each viewpoint, we compute dust extinction from the simulated gas metallicity distribution and apply a simple error model to produce a synthetic Gaia-like survey that includes both observational properties and a pointer to the generating star particle. We provide the complete simulation snapshot at z = 0 for each simulated galaxy. We describe data access points, the data model, and plans for future upgrades. These synthetic surveys provide a tool for the scientific community to test analysis methods and interpret Gaia data.

While many tensions between Local Group (LG) satellite galaxies and ⋀ cold dark matter cosmology have been alleviated through recent cosmological simulations, the spatial distribution of satellites remains an important test of physical models and physical versus numerical disruption in simulations. Using the FIRE-2 cosmological zoom-in baryonic simulations, we examine the radial distributions of satellites with M∗ > 105 M☉ around eight isolated Milky Way (MW) mass host galaxies and four hosts in LG-like pairs. We demonstrate that these simulations resolve the survival and physical destruction of satellites with M∗ ≲ 105 M☉. The simulations broadly agree with LG observations, spanning the radial profiles around the MW and M31. This agreement does not depend strongly on satellite mass, even at distances ≾100 kpc. Host-to-host variation dominates the scatter in satellite counts within 300 kpc of the hosts, while time variation dominates scatter within 50 kpc. More massive host galaxies within our sample have fewer satellites at small distances, likely because of enhanced tidal destruction of satellites via the baryonic discs of host galaxies. Furthermore, we quantify and provide fits to the tidal depletion of subhaloes in baryonic relative to dark matter-only simulations as a function of distance. Our simulated profiles imply observational incompleteness in the LG even at M∗ ≳ 105 M☉: we predict 2-10 such satellites to be discovered around the MW and possibly 6-9 around M31. To provide cosmological context, we compare our results with the radial profiles of satellites around MW analogues in the SAGA survey, finding that our simulations are broadly consistent with most SAGA systems.

We study a suite of extremely high-resolution cosmological Feedback in Realistic Environments simulations of dwarf galaxies (Mhalo ≲ 1010 M☉), run to z = 0 with 30 M☉ resolution, sufficient (for the first time) to resolve the internal structure of individual supernovae remnants within the cooling radius. Every halo with Mhalo ≳ 108.6 M☉ is populated by a resolved stellar galaxy, suggesting very low-mass dwarfs may be ubiquitous in the field. Our ultra-faint dwarfs (UFDs; M∗ < 105 M☉) have their star formation (SF) truncated early (z ≳ 2), likely by reionization, while classical dwarfs (M∗ > 105 M☉) continue forming stars to z < 0.5. The systems have bursty star formation histories, forming most of their stars in periods of elevated SF strongly clustered in both space and time. This allows our dwarf with M∗/Mhalo > 10−4 to form a dark matter core >200 pc, while lower mass UFDs exhibit cusps down to ≲100 pc, as expected from energetic arguments. Our dwarfs with M∗ > 104 M☉ have half-mass radii (R1/2) in agreement with Local Group (LG) dwarfs (dynamical mass versus R1/2 and stellar rotation also resemble observations). The lowest mass UFDs are below surface brightness limits of current surveys but are potentially visible in next-generation surveys (e.g. LSST). The stellar metallicities are lower than in LG dwarfs; this may reflect pre-enrichment of the LG by the massive hosts or Pop-III stars. Consistency with lower resolution studies implies that our simulations are numerically robust (for a given physical model).

We present a suite of FIRE-2 cosmological zoom-in simulations of isolated field dwarf galaxies, all with masses of Mhalo ≈ 1010 M at z = 0, across a range of dark matter models. For the first time, we compare how both self-interacting dark matter (SIDM) and/or warm dark matter (WDM) models affect the assembly histories as well as the central density structure in fully hydrodynamical simulations of dwarfs. Dwarfs with smaller stellar half-mass radii (r1/2 < 500 pc) have lower σ/Vmax ratios, reinforcing the idea that smaller dwarfs may reside in haloes that are more massive than is naively expected. The majority of dwarfs simulated with self-interactions actually experience contraction of their inner density profiles with the addition of baryons relative to the cores produced in dark-matter-only runs, though the simulated dwarfs are always less centrally dense than in CDM. The V1/2–r1/2 relation across all simulations is generally consistent with observations of Local Field dwarfs, though compact objects such as Tucana provide a unique challenge. Overall, the inclusion of baryons substantially reduces any distinct signatures of dark matter physics in the observable properties of dwarf galaxies. Spatially resolved rotation curves in the central regions (

We study star formation histories (SFHs) of 500 dwarf galaxies (stellar mass M∗ = 105–109 M) 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 versus central’ and ‘LG versus individual MW versus 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∗ = 106–107 M 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.

The Gaia era opens new possibilities for discovering the remnants of disrupted satellite galaxies in the solar neighborhood. If the population of local accreted stars is correlated with the dark matter sourced by the same mergers, one can then map the dark matter distribution directly. Using two cosmological zoom-in hydrodynamic simulations of Milky-Way-mass galaxies from the Latte suite of the Fire-2 simulations, we find a strong correlation between the velocity distribution of stars and dark matter at the solar circle that were accreted from luminous satellites. This correspondence holds for dark matter that is either relaxed or in a kinematic substructure called debris flow, and is consistent between two simulated hosts with different merger histories. The correspondence is more problematic for streams because of possible spatial offsets between the dark matter and stars. We demonstrate how to reconstruct the dark matter velocity distribution from the observed properties of the accreted stellar population by properly accounting for the ratio of stars to dark matter contributed by individual mergers. This procedure does not account for the dark matter that originates from nonluminous satellites, which may constitute a nontrivial fraction of the local contribution. After validating this method using the Fire-2 simulations, we apply it to the Milky Way and use it to recover the dark matter velocity distribution associated with the recently discovered stellar debris field in the solar neighborhood. Based on results from Gaia, we estimate that of the local dark matter that is accreted from luminous mergers is in debris flow.

We use a particle tracking analysis to study the origins of the circumgalactic medium (CGM), separating it into (1) accretion from the intergalactic medium (IGM), (2) wind from the central galaxy, and (3) gas ejected from other galaxies. Our sample consists of 21 FIRE-2 simulations, spanning the halo mass range Mh ∼ 1010–1012 M, and we focus on z = 0.25 and z = 2. Owing to strong stellar feedback, only ∼L haloes retain a baryon mass 50 per cent of their cosmic budget. Metals are more efficiently retained by haloes, with a retention fraction 50 per cent. Across all masses and redshifts analysed 60 per cent of the CGM mass originates as IGM accretion (some of which is associated with infalling haloes). Overall, the second most important contribution is wind from the central galaxy, though gas ejected or stripped from satellites can contribute a comparable mass in ∼L haloes. Gas can persist in the CGM for billions of years, resulting in well mixed-halo gas. Sightlines through the CGM are therefore likely to intersect gas of multiple origins. For low-redshift ∼L haloes, cool gas (T < 104.7 K) is distributed on average preferentially along the galaxy plane, however with strong halo-to-halo variability. The metallicity of IGM accretion is systematically lower than the metallicity of winds (typically by 1 dex), although CGM and IGM metallicities depend significantly on the treatment of subgrid metal diffusion. Our results highlight the multiple physical mechanisms that contribute to the CGM and will inform observational efforts to develop a cohesive picture.

We present a new set of high-resolution hydrodynamic cosmological zoom-in simulations that apply the Feedback In Realistic Environments physics to both Local Group (LG)-like and isolated MilkyWay (MW)-like volumes (10 host systems in total with a baryonic particle mass ≃3500-7000M⊙).We study the stellar mass functions, circular velocity or mass profiles, and velocity dispersions of the dwarf galaxy populations. The simulations reproduce the stellar mass function and central densities ofMWsatellite dwarfs forM∗≥105.5M⊙ and predict the existence of ∼3 unidentified galaxies with M∗∼ 105M⊙ within 300 kpc of the MW. Overall, we find no evidence for the classical missing satellites or too-big-to-fail (TBTF) problems for satellite galaxies in our sample. Among the satellites, TBTF is resolved primarily by subhalo disruption and overall mass-loss; central density profiles of subhaloes are of secondary importance. For non-satellite galaxies, our LG-like simulations predict as many as ∼10 as-ofyet unseen galaxies at distances 0.3-1Mpc from both hosts, with M∗≃ 105-6M⊙ (in haloes with Vmax ∼ 20 km s-1), albeit with large halo-to-halo variance. None of our simulations produces a compact, baryon-dominated, high-density dwarf elliptical-type galaxy (with Vcirc ≳ 35 km s-1 at r < 1 kpc), of which six may appear in the LG (but none in the MW). It may therefore remain a challenge to reproduce the full diversity of the dwarf population, including both the highest and lowest density systems.

We present the late-time evolution of m12m, a cosmological simulation of a Milky Way-like galaxy from the FIRE project. The simulation forms a bar after redshift z = 0.2. We show that the evolution of the model exhibits behaviours typical of kinematic fractionation, with a bar weaker in older populations, an X-shape traced by the younger, metal-rich populations, and a prominent X-shape in the edge-on mean metallicity map. Because of the late formation of the bar in m12m, stars forming after 10 Gyr (z = 0.34) significantly contaminate the bulge, at a level higher than is observed at high latitudes in the Milky Way, implying that its bar cannot have formed as late as in m12m. We also study the model's vertex deviation of the velocity ellipsoid as a function of stellar metallicity and age in the equivalent of Baade's Window. The formation of the bar leads to a non-zero vertex deviation. We find that metal-rich stars have a large vertex deviation (∼40◦), which becomes negligible for metal-poor stars, a trend also found in the Milky Way, despite not matching in detail. We demonstrate that the vertex deviation also varies with stellar age and is large for stars as old as 9 Gyr, while 13 Gyr old stars have negligible vertex deviation. When we exclude stars that have been accreted, the vertex deviation is not significantly changed, demonstrating that the observed variation of vertex deviation with metallicity is not necessarily due to an accreted population.

We study the impact of a warm dark matter (WDM) cosmology on dwarf galaxy formation through a suite of cosmological hydrodynamical zoom-in simulations of Mhalo ≈ 1010 M dark matter haloes as part of the Feedback in Realistic Environments (FIRE) project. A main focus of this paper is to evaluate the combined effects of dark matter physics and stellar feedback on the well-known small-scale issues found in cold dark matter (CDM) models. We find that the z = 0 stellar mass of a galaxy is strongly correlated with the central density of its host dark matter halo at the time of formation, zf, in both CDM and WDM models. WDM haloes follow the same M(z = 0)–Vmax(zf) relation as in CDM, but they form later, are less centrally dense, and therefore contain galaxies that are less massive than their CDM counterparts. As a result, the impact of baryonic effects on the central gravitational potential is typically diminished relative to CDM. However, the combination of delayed formation in WDM and energy input from stellar feedback results in dark matter profiles with lower overall densities. The WDM galaxies studied here have a wider diversity of star formation histories (SFHs) than the same systems simulated in CDM, and the two lowest M WDM galaxies form all of their stars at late times. The discovery of young ultrafaint dwarf galaxies with no ancient star formation – which do not exist in our CDM simulations – would therefore provide evidence in support of WDM.

We use cosmological hydrodynamical simulations of Milky Way-mass galaxies from the FIRE project to evaluate various strategies for estimating the mass of a galaxy's stellar halo from deep, integrated-light images. We find good agreement with integrated-light observations if we mimic observational methods to measure the mass of the stellar halo by selecting regions of an image via projected radius relative to the disk scale length or by their surface density in stellar mass. However, these observational methods systematically underestimate the accreted stellar component, defined in our (and most) simulations as the mass of stars formed outside of the host galaxy, by up to a factor of 10, since the accreted component is centrally concentrated and therefore substantially obscured by the galactic disk. Furthermore, these observational methods introduce spurious dependencies of the estimated accreted stellar component on the stellar mass and size of galaxies that can obscure the trends in accreted stellar mass predicted by cosmological simulations, since we find that in our simulations, the size and shape of the central galaxy are not strongly correlated with the assembly history of the accreted stellar halo. This effect persists whether galaxies are viewed edge-on or face-on. We show that metallicity or color information may provide a way to more cleanly delineate in observations the regions dominated by accreted stars. Absent additional data, we caution that estimates of the mass of the accreted stellar component from single-band images alone should be taken as lower limits.

Binary black holes are the primary endpoint of massive stars. Their properties provide a unique opportunity to constrain binary evolution, which remains poorly understood. We predict the main properties of binary black holes and their merger products in/around the Milky Way. We present the first combination of a high-resolution cosmological simulation of a Milky Way-mass galaxy with a binary population synthesis model in this context. The hydrodynamic simulation, taken from the FIRE project, provides a cosmologically realistic star formation history for the galaxy, its stellar halo, and satellites. During post-processing, we apply a metallicity-dependent evolutionary model to the star particles to produce individual binary black holes. We find that 7 × 105 binary black holes have merged in the model Milky Way, and 1.2 × 106 binaries are still present, with a mean mass of 28M⊙. Because the black hole progenitors are strongly biased towards low-metallicity stars, half reside in the stellar halo and satellites and a third were formed outside the main galaxy. The numbers and mass distribution of the merged systems is broadly compatible with the LIGO/Virgo detections. Our simplified binary evolution models predict that LISA will detect more than 20 binary black holes, but that electromagnetic observations will be challenging. Our method will allow for constraints on the evolution of massive binaries based on comparisons between observations of compact objects and the predictions of varying binary evolution models. We provide online data of our star formation model and binary black hole distribution.

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. Central (~kpc) mass concentrations in massive (> L*) galaxies are sensitive to numerics (via trapping/recycling of winds in hot haloes). Multiple feedback mechanisms play key roles: supernovae regulate stellar masses/winds; stellar mass-loss fuels late star formation; radiative feedback suppresses accretion on to dwarfs and instantaneous star formation in discs. We provide all initial conditions and numerical algorithms used.

We test if the cosmological zoom-in simulations of isolated galaxies from the FIRE project reproduce the properties of ultra diffuse galaxies (UDGs). We show that outflows that dynamically heat galactic stars, together with a passively aging stellar population after imposed quenching, naturally reproduce the observed population of red UDGs, without the need for high spin haloes, or dynamical influence from their host cluster. We reproduce the range of surface brightness, radius, and absolute magnitude of the observed red UDGs by quenching simulated galaxies at a range of different times. They represent a mostly uniform population of dark matter-dominated dwarf galaxies with M * ~ 108 M☉, low metallicity, and a broad range of ages; the more massive the UDGs, the older they are. The most massive red UDG in our sample(M * ~ 3 × 108 M☉) requires quenching at z ~ 3 when its halo reached M h ~ 1011M☉. Our simulated UDGs form with normal stellar-to-halo ratios and match the central enclosed masses and the velocity dispersions of the observed UDGs. Enclosed masses remain largely fixed across a broad range of quenching times because the central regions of their dark matter haloes complete their growth early. If our simulated dwarfs are not quenched, they evolve into bluer low surface brightness galaxies with M/L similar to observed field dwarfs. While our simulation sample covers a limited range of formation histories and halo masses, we predict that UDG is a common, and perhaps even dominant, galaxy type around M * ~ 108 M☉, both in the field and in clusters.

We study the morphologies and sizes of galaxies at z ≥ 5 using high-resolution cosmological zoom-in simulations from the Feedback In Realistic Environments project. The galaxies show a variety of morphologies, from compact to clumpy to irregular. The simulated galaxies have more extended morphologies and larger sizes when measured using rest-frame optical B-band light than rest-frameUVlight; sizes measured from stellarmass surface density are even larger. The UV morphologies are usually dominated by several small, bright young stellar clumps that are not always associated with significant stellarmass. The B-band light traces stellarmass better than the UV, but it can also be biased by the bright clumps. At all redshifts, galaxy size correlates with stellar mass/luminosity with large scatter. The half-light radii range from 0.01 to 0.2 arcsec (0.05-1 kpc physical) at fixed magnitude. At z ≥ 5, the size of galaxies at fixed stellar mass/luminosity evolves as (1 + z)-m, with m ~ 1-2. For galaxies less massive than M* ~ 108M⊙, the ratio of the half-mass radius to the halo virial radius is ~10 per cent and does not evolve significantly at z = 5-10; this ratio is typically 1-5 per cent for more massive galaxies. A galaxy's 'observed' size decreases dramatically at shallower surface brightness limits. This effect may account for the extremely small sizes of z ≥ 5 galaxies measured in the Hubble Frontier Fields. We provide predictions for the cumulative light distribution as a function of surface brightness for typical galaxies at z = 6.

We identify subhalos in dark matter-only (DMO) zoom-in simulations that are likely to be disrupted due to baryonic effects by using a random forest classifier trained on two hydrodynamic simulations of Milky Way (MW)-mass host halos from the Latte suite of the Feedback in Realistic Environments (FIRE) project. We train our classifier using five properties of each disrupted and surviving subhalo: pericentric distance and scale factor at first pericentric passage after accretion and scale factor, virial mass, and maximum circular velocity at accretion. Our five-property classifier identifies disrupted subhalos in the FIRE simulations with an 85% out-of-bag classification score. We predict surviving subhalo populations in DMO simulations of the FIRE host halos, finding excellent agreement with the hydrodynamic results; in particular, our classifier outperforms DMO zoom-in simulations that include the gravitational potential of the central galactic disk in each hydrodynamic simulation, indicating that it captures both the dynamical effects of a central disk and additional baryonic physics. We also predict surviving subhalo populations for a suite of DMO zoom-in simulations of MW-mass host halos, finding that baryons impact each system consistently and that the predicted amount of subhalo disruption is larger than the host-to-host scatter among the subhalo populations. Although the small size and specific baryonic physics prescription of our training set limits the generality of our results, our work suggests that machine-learning classification algorithms trained on hydrodynamic zoom-in simulations can efficiently predict realistic subhalo populations.

Using a state-of-the-art cosmological simulation of merging proto-galaxies at high redshift from the FIRE project, with explicit treatments of star formation and stellar feedback in the interstellar medium, we investigate the formation of star clusters and examine one of the formation hypotheses of present-day metal-poor globular clusters. We find that frequent mergers in high-redshift proto-galaxies could provide a fertile environment to produce long-lasting bound star clusters. The violent merger event disturbs the gravitational potential and pushes a large gas mass of ≳ 105-6M⊙ collectively to high density, at which point it rapidly turns into stars before stellar feedback can stop star formation. The high dynamic range of the reported simulation is critical in realizing such dense star-forming clouds with a small dynamical time-scale, tff ≲ 3 Myr, shorter than most stellar feedback time-scales. Our simulation then allows us to trace how clusters could become virialized and tightly bound to survive for up to ~420 Myr till the end of the simulation. Because the cluster's tightly bound core was formed in one short burst, and the nearby older stars originally grouped with the cluster tend to be preferentially removed, at the end of the simulation the cluster has a small age spread.

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 study the z = 0 gas kinematics, morphology and angular momentum content of isolated galaxies in a suite of cosmological zoom-in simulations from the FIRE project spanning Mstar = 106-11M⊙. Gas becomes increasingly rotationally supported with increasing galaxy mass. In the lowest mass galaxies (Mstar < 108M⊙), gas fails to form a morphological disc and is primarily dispersion and pressure supported. At intermediate masses (Mstar = 108-10M⊙), galaxies display a wide range of gas kinematics and morphologies, from thin, rotating discs to irregular spheroids with negligible net rotation. All the high-mass (Mstar = 1010-11M⊙) galaxies form rotationally supported gas discs. Many of the haloes whose galaxies fail to form discs harbour high angular momentum gas in their circumgalactic medium. The ratio of the specific angular momentum of gas in the central galaxy to that of the dark matter halo increases significantly with galaxy mass, from 〈jgas〉/〈jDM〉 ~ 0.1 at Mstar = 106-7M⊙ to 〈jgas〉/〈jDM〉 ~ 2 at Mstar = 1010-11M⊙. The reduced rotational support in the lowest mass galaxies owes to (a) stellar feedback and the UV background suppressing the accretion of high angular momentum gas at late times, and (b) stellar feedback driving large non-circular gas motions. We broadly reproduce the observed scaling relations between galaxy mass, gas rotation velocity, size and angular momentum, but may somewhat underpredict the incidence of disky, high angular momentum galaxies at the lowest observed masses (Mstar = (106- 2 × 107) M⊙). Stars form preferentially from low angular momentum gas near the galactic centre and are less rotationally supported than gas. The common assumption that stars follow the same rotation curve as gas thus substantially overestimates the simulated galaxies' stellar angular momentum, particularly at low masses.

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.

Motivated by the recent discovery of several dwarfs near the Large Magellanic Cloud (LMC), we study the accretion of massive satellites onto Milky Way (MW)/M31-like haloes using the ELVIS suite of N-body simulations. We identify 25 surviving LMC-mass subhaloes, and investigate the lower-mass satellites that were associated with these subhaloes before they fell into the MW/M31 haloes. Typically, 7 per cent of the overall z = 0 satellite population of MW/M31 haloes were in a surviving LMC-group before falling into the MW/M31 halo. This fraction can vary between 1 and 25 per cent, being higher for groups with higher mass and/or more recent infall times. Groups of satellites disperse rapidly in phase space after infall, and their distances and velocities relative to the group centre become statistically similar to the overall satellite population after 4-8 Gyr. We quantify the likelihood that satellites were associated with an LMC-mass group as a function of both distance and velocity relative to the LMC at z = 0. The close proximity in distance of the nine Dark Energy Survey candidate dwarf galaxies to the LMC suggest that ~2-4 are likely associated with the LMC. Furthermore, if several of these dwarfs are genuine members, then the LMC-group probably fell into the MW very recently, ≲2 Gyr ago. If the connection with the LMC is established with follow-up velocity measurements, these 'satellites of satellites' represent prime candidates to study the effects of group pre-processing on lower mass dwarfs.

In the Local Group (LG), almost all satellite dwarf galaxies that are within the virial radius of the Milky Way (MW) and Andromeda (M31) exhibit strong environmental influence. The orbital histories of these satellites provide the key to understanding the role of the MW/M31 halo, lower-mass groups, and cosmic reionization on the evolution of dwarf galaxies. We examine the virial-infall histories of satellites with Mstar = 103-9 M⊙ using the Exploring the Local Volume in Simulations suite of cosmological zoom-in dissipationless simulations of 48 MW/M31-like halos. Satellites at z = 0 fell into the MW/M31 halos typically 5-8 Gyr ago at z = 0.5-1. However, they first fell into any host halo typically 7-10 Gyr ago at z = 0.7-1.5. This difference arises because many satellites experienced "group preprocessing" in another host halo, typically of Mvir ∼ 1010-12 M⊙, before falling into the MW/M31 halos. Satellites with lower mass and/or those closer to the MW/M31 fell in earlier and are more likely to have experienced group preprocessing; half of all satellites with Mstar < 106 M⊙were preprocessed in a group. Infalling groups also drive most satellite-satellite mergers within the MW/M31 halos. Finally, none of the surviving satellites at z = 0 were within the virial radius of their MW/M31 halo during reionization (z > 6), and only

We investigate the effect of self-interacting dark matter (SIDM) on the density profiles of Vmax ≃ 40km s-1 isolated dwarf dark matter haloes-the scale of relevance for the too big to fail problem (TBTF)-using very high resolution cosmological zoom simulations. Each halo has millions of particles within its virial radius. We find that SIDM models with cross-sections per unitmass spanning the range s/m=0.5-50 cm2 g-1 alleviate TBTF and produce constantdensity cores of size 300-1000 pc, comparable to the half-light radii of M* ~ 105-7M⊙ dwarfs. The largest, lowest density cores develop for cross-sections in the middle of this range, σ/m ~ 5-10 cm2 g-1. Our largest SIDM cross-section run (σ/m = 50 cm2 g-1) develops a slightly denser core owing to mild core-collapse behaviour, but it remains less dense than the cold dark matter case and retains a constant-density core profile. Our work suggests that SIDM cross-sections as large or larger than 50 cm2 g-1 remain viable on velocity scales of dwarf galaxies (vrms ~ 40 km s-1). The range of SIDM cross-sections that alleviate TBTF and the cusp/core problem spans at least two orders of magnitude and therefore need not be particularly fine-tuned.

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 M star > 106 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.

The BICEP2 results, when interpreted as a gravitational wave signal and combined with other cosmic microwave background data, suggest a roll-off in power towards small scales in the primordial matter power spectrum. Among the simplest possibilities is a running of the spectral index. Here we show that the preferred level of running alleviates small-scale issues within the ΛCDM model, more so even than viable WDM models. We use cosmological zoom-in simulations of a MilkyWay-sized halo along with full-box simulations to compare predictions among four separate cosmologies: a BICEP2-inspired running index model (αs = -0.024), two fixed-tilt ΛCDM models motivated by Planck, and a 2.6 keV thermal WDM model. We find that the running BICEP2 model reduces the central densities of large dwarf-sized haloes (Vmax ~ 30-80 km s-1) and alleviates the too-big-to-fail problem significantly compared to our adopted Planck andWDMcases. Further, the BICEP2 model suppresses the count of small subhaloes by ~50 per cent relative to Planck models, and yields a significantly lower 'boost' factor for dark matter annihilation signals. Our findings highlight the need to understand the shape of the primordial power spectrum in order to correctly interpret small-scale data.

We show that the canonical oscillation-based (nonresonant) production of sterile neutrino dark matter is inconsistent at >99% confidence with observations of galaxies in the Local Group. We set lower limits on the nonresonant sterile neutrino mass of 2.5 keV (equivalent to 0.7 keV thermal mass) using phase-space densities derived for dwarf satellite galaxies of the Milky Way as well as limits of 8.8 keV (equivalent to 1.8 keV thermal mass) based on subhalo counts of N-body simulations of M 31 analogs. Combined with improved upper mass limits derived from significantly deeper x-ray data of M 31 with full consideration for background variations, we show that there remains little room for nonresonant production if sterile neutrinos are to explain 100% of the dark matter abundance. Resonant and nonoscillation sterile neutrino production remain viable mechanisms for generating sufficient dark matter sterile neutrinos. © 2014 American Physical Society.

The predominantly ancient stellar populations observed in the lowest mass galaxies (i.e. ultra-faint dwarfs) suggest that their star formation was suppressed by reionization. Most of the well-studied ultra-faint dwarfs, however, are within the central half of the Milky Way dark matter halo, such that they are consistent with a population that was accreted at early times and thus potentially quenched via environmental processes. To study the potential role of environment in suppressing star formation on the smallest scales, we utilize the Exploring the Local Volume in Simulations suite of N -body simulations to constrain the distribution of infall times for low-mass subhaloes likely to host the ultra-faint population. For the ultra-faint satellites of the Milky Way with star formation histories inferred from Hubble Space Telescope imaging, we find that environment is highly unlikely to play a dominant role in quenching their star formation. Even when including the potential effects of pre-processing, there is a |${\lesssim } 0.1{{\ \rm per\ cent}}$| probability that environmental processes quenched all of the known ultra-faint dwarfs early enough to explain their observed star formation histories. Instead, we argue for a mass floor in the effectiveness of satellite quenching at roughly |${M}_{\star }\sim 10^{5}\, {\rm M}_{\odot }$|⁠, below which star formation in surviving galaxies is globally suppressed by reionization. We predict a large population of quenched ultra-faint dwarfs in the Local Field (1 < R / R vir < 2), with as many as ∼250 to be discovered by future wide-field imaging surveys. [ABSTRACT FROM AUTHOR]

We use hydrodynamic cosmological zoom-in simulations from the Feedback in Realistic Environments project to explore the morphologies and kinematics of 15 Milky Way (MW)-mass galaxies. Our sample ranges from compact, bulge-dominated systems with 90 per cent of their stellar mass within 2.5 kpc to well-ordered discs that reach ≳15 kpc. The gas in our galaxies always forms a thin, rotation-supported disc at z = 0, with sizes primarily determined by the gas mass. For stars, we quantify kinematics and morphology both via the fraction of stars on disc-like orbits and with the radial extent of the stellar disc. In this mass range, stellar morphology and kinematics are poorly correlated with the properties of the halo available from dark matter-only simulations (halo merger history, spin, or formation time). They more strongly correlate with the gaseous histories of the galaxies: those that maintain a high gas mass in the disc after z ~ 1 develop well-ordered stellar discs. The best predictor of morphology we identify is the spin of the gas in the halo at the time the galaxy formed 1/2 of its stars (i.e. the gas that builds the galaxy). High- z mergers, before a hot halo emerges, produce some of the most massive bulges in the sample (from compact discs in gas-rich mergers), while later-forming bulges typically originate from internal processes, as satellites are stripped of gas before the galaxies merge. Moreover, most stars in z = 0 MW-mass galaxies (even z = 0 bulge stars) form in a disc: ≳60-90 per cent of stars begin their lives rotationally supported.

The oldest stars in the Milky Way (MW) bear imprints of the Galaxy's early assembly history. We use FIRE cosmological zoom-in simulations of three MW-mass disc galaxies to study the spatial distribution, chemistry, and kinematics of the oldest surviving stars ( z form ≳ 5) in MW-like galaxies. We predict the oldest stars to be less centrally concentrated at z = 0 than stars formed at later times as a result of two processes. First, the majority of the oldest stars are not formed in situ but are accreted during hierarchical assembly. These ex situ stars are deposited on dispersion-supported, halo-like orbits but dominate over old stars formed in situ in the solar neighbourhood, and in some simulations, even in the galactic centre. Secondly, old stars formed in situ are driven outwards by bursty star formation and energetic feedback processes that create a time-varying gravitational potential at z ≳ 2, similar to the process that creates dark matter cores and expands stellar orbits in bursty dwarf galaxies. The total fraction of stars that are ancient is more than an order of magnitude higher for sight lines away from the bulge and inner halo than for inward-looking sight lines. Although the task of identifying specific stars as ancient remains challenging, we anticipate that million-star spectral surveys and photometric surveys targeting metal-poor stars already include hundreds of stars formed before z = 5. We predict most of these targets to have higher metallicity (-3 < [Fe/H] < -2) than the most extreme metal-poor stars.

In the local Universe, there is a strong division in the star-forming properties of low-mass galaxies, with star formation largely ubiquitous amongst the field population while satellite systems are predominantly quenched. This dichotomy implies that environmental processes play the dominant role in suppressing star formation within this low-mass regime ( M ★ ~ 10 5.5-8 M ☉ ). As shown by observations of the Local Volume, however, there is a non-negligible population of passive systems in the field, which challenges our understanding of quenching at low masses. By applying the satellite quenching models of Fillingham et al. (2015) to subhalo populations in the Exploring the Local Volume In Simulations suite, we investigate the role of environmental processes in quenching star formation within the nearby field. Using model parameters that reproduce the satellite quenched fraction in the Local Group, we predict a quenched fraction - due solely to environmental effects - of ~0.52 ± 0.26 within 1 < R/R vir < 2 of the Milky Way and M31. This is in good agreement with current observations of the Local Volume and suggests that the majority of the passive field systems observed at these distances are quenched via environmental mechanisms. Beyond 2 R vir , however, dwarf galaxy quenching becomes difficult to explain through an interaction with either the Milky Way or M31, such that more isolated, field dwarfs may be self-quenched as a result of star-formation feedback.