Explore the vast range of titles available.

Galaxies and their dark-matter haloes have posed several challenges to the dark energy plus cold dark matter ( Λ CDM) cosmological model. These discrepancies between observations and theory intensify for the lowest-mass (‘dwarf’) galaxies. Λ CDM predictions for the number, spatial distribution and internal structure of low-mass dark-matter haloes have historically been at odds with observed dwarf galaxies, but this is partially expected, because many predictions modelled only the dark-matter component. Any robust Λ CDM prediction must include, hand in hand, a model for galaxy formation to understand how baryonic matter populates and affects dark-matter haloes. In this Review, we consider the most notable challenges to Λ CDM regarding dwarf galaxies, and we discuss how recent cosmological numerical simulations have pinpointed baryonic solutions to these challenges. We identify remaining tensions, including the diversity of the inner dark-matter content, planes of satellites, stellar morphologies and star-formation quenching. Their resolution, or validation as actual problems with Λ CDM, will probably require both refining of galaxy-formation models and improving numerical accuracy in simulations. Galaxies and their dark-matter haloes have posed several challenges to the dark energy plus cold dark matter ( Λ CDM) cosmological model. This Review discusses the most notable challenges to Λ CDM regarding dwarf galaxies and the insights provided by recent cosmological numerical simulations.

Recent discoveries of a significant population of bright galaxies at cosmic dawn z≳10 have enabled critical tests of cosmological galaxy formation models. In particular, the bright end of the galaxys’ UV luminosity functions (UVLFs) appear higher than predicted by many models. Using approximately 25,000 galaxy snapshots at 8 ≤ z ≤ 12 in a suite of FIRE-2 cosmological “zoom-in” simulations from the Feedback in Realistic Environments (FIRE) project, we show that the observed abundance of UV-bright galaxies at cosmic dawn is reproduced in these simulations with a multichannel implementation of standard stellar feedback processes, without any fine-tuning. Notably, we find no need to invoke previously suggested modifications, such as a nonstandard cosmology, a top-heavy stellar initial mass function, or a strongly enhanced star formation efficiency. We contrast the UVLFs predicted by bursty star formation in these original simulations to those derived from star formation histories (SFHs) smoothed over prescribed timescales (e.g., 100 Myr). The comparison demonstrates that the strongly time-variable SFHs predicted by the FIRE simulations play a key role in correctly reproducing the observed, bright-end UVLFs at cosmic dawn: the bursty SFHs induce order-or-magnitude changes in the abundance of UV-bright (M UV ≲ −20) galaxies at z ≳ 10. The predicted bright-end UVLFs are consistent with both the spectroscopically confirmed population and the photometrically selected candidates. We also find good agreement between the predicted and observationally inferred integrated UV luminosity densities, which evolve more weakly with redshift in FIRE than suggested by some other models.

Spatial patterns of stellar elemental abundances encode rich information about a galaxy’s formation history. We analyze the radial, vertical, and azimuthal variations of metals in stars, both today and at formation, in the FIRE-2 cosmological simulations of Milky Way (MW)-mass galaxies, and we compare them with the MW. The radial gradient today is steeper (more negative) for younger stars, which agrees with the MW, although radial gradients are shallower in FIRE-2. Importantly, this age dependence was present already at birth: radial gradients today are only modestly (≲0.01 dex kpc−1) shallower than at birth. Disk vertical settling gives rise to negative vertical gradients across all stars, but vertical gradients of mono-age stellar populations are weak. Similar to the MW, vertical gradients in FIRE-2 are shallower at larger radii, but they are overall shallower in FIRE-2. This vertical dependence was present already at birth: vertical gradients today are only modestly (≲0.1 dex kpc−1) shallower than at birth. Azimuthal scatter is nearly constant with radius, and it is nearly constant with age ≲8 Gyr ago but increases for older stars. Azimuthal scatter is slightly larger (≲0.04 dex) today than at formation. Galaxies with larger azimuthal scatter have a stronger radial gradient, implying that azimuthal scatter today arises primarily from the radial redistribution of gas and stars. Overall, spatial variations of stellar metallicities show only modest differences between formation and today; spatial variations today primarily reflect the conditions of stars at birth, with spatial redistribution of stars after birth contributing secondarily.

We describe a public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation (available at http://flathub.flatironinstitute.org/fire) from the Feedback In Realistic Environments (FIRE) project. FIRE-2 simulations achieve parsec-scale resolution to explicitly model the multiphase interstellar medium while implementing direct models for stellar evolution and feedback, including stellar winds, core-collapse and Type Ia supernovae, radiation pressure, photoionization, and photoelectric heating. We release complete snapshots from three suites of simulations. The first comprises 20 simulations that zoom in on 14 Milky Way (MW)–mass galaxies, five SMC/LMC-mass galaxies, and four lower-mass galaxies including one ultrafaint; we release 39 snapshots across z = 0–10. The second comprises four massive galaxies, with 19 snapshots across z = 1–10. Finally, a high-redshift suite comprises 22 simulations, with 11 snapshots across z = 5–10. Each simulation also includes dozens of resolved lower-mass (satellite) galaxies in its zoom-in region. Snapshots include all stored properties for all dark matter, gas, and star particles, including 11 elemental abundances for stars and gas, and formation times (ages) of star particles. We also release accompanying (sub)halo catalogs, which include galaxy properties and member star particles. For the simulations to z = 0, including all MW-mass galaxies, we release the formation coordinates and an “ex situ” flag for all star particles, pointers to track particles across snapshots, catalogs of stellar streams, and multipole basis expansions for the halo mass distributions. We describe publicly available python packages for reading and analyzing these simulations.

In the Gaia era it is increasingly apparent that traditional static, parameterized models are insufficient to describe the mass distribution of our complex, dynamically evolving Milky Way (MW). In this work, we compare different time-evolving and time-independent representations of the gravitational potentials of simulated MW-mass galaxies from the FIRE-2 suite of cosmological-baryonic simulations. Using these potentials, we calculate actions for star particles in tidal streams around three galaxies with varying merger histories at each snapshot from 7 Gyr ago to the present day. We determine the action-space coherence preserved by each model using the Kullback–Leibler divergence to gauge the degree of clustering in actions and the relative stability of the clusters over time. We find that all models produce a clustered action space for simulations with no significant mergers. However, a massive (mass ratio prior to infall more similar than 1:8) interacting galaxy not present in the model will result in mischaracterized orbits for stars most affected by the interaction. The locations of the action space clusters (i.e., the orbits of the stream stars) are only preserved by the time-evolving model, while the time-independent models can lose significant amounts of information as soon as 0.5–1 Gyr ago, even if the system does not undergo a significant merger. Our results imply that reverse-integration of stream orbits in the MW using a fixed potential is likely to give incorrect results if integrated longer than 0.5 Gyr into the past.

The physical mechanisms responsible for bar formation and destruction in galaxies remain a subject of debate. While we have gained valuable insight into how bars form and evolve from isolated idealized simulations, in the cosmological domain, galactic bars evolve in complex environments, with mergers and gas accretion events occurring in the presence of the turbulent interstellar medium with multiple star formation episodes, in addition to coupling with their host galaxies’ dark matter halos. We investigate the bar formation in 13 Milky Way–mass galaxies from the Feedback in Realistic Environments (FIRE-2) cosmological zoom-in simulations. 8 of the 13 simulated galaxies form bars at some point during their history: three from tidal interactions and five from internal evolution of the disk. The bars in FIRE-2 are generally shorter than the corotation radius (mean bar radius ∼1.53 kpc), have a wide range of pattern speeds (36–97 km s−1 kpc−1), and live for a wide range of dynamical times (2–160 bar rotations). We find that the bar formation in FIRE-2 galaxies is influenced by satellite interactions and the stellar-to-dark-matter mass ratio in the inner galaxy, but neither is a sufficient condition for bar formation. Bar formation is more likely to occur, with the bars formed being stronger and longer-lived, if the disks are kinematically cold; galaxies with high central gas fractions and/or vigorous star formation, on the other hand, tend to form weaker bars. In the case of the FIRE-2 galaxies, these properties combine to produce ellipsoidal bars with strengths A 2/A 0 ∼ 0.1–0.2.

Within the Milky Way (MW), younger stellar populations exhibit steeper (more negative) metallicity radial gradients; the origin of this trend remains debated. The FIRE-2 cosmological simulations of MW-mass galaxies show the same trend as the MW, which in FIRE-2 arises because the metallicity gradient of the interstellar medium (ISM), and thus of stars at birth, became steeper over time. We seek to understand this evolution in the context of inside-out radial growth of galaxies. Most FIRE-2 galaxies grew radially from the inside out in both gas and stars; specifically, their surface density profiles, Σ(R), became shallower over time. Combined with a realized superlinear (Kennicutt–Schmidt-like) relation between star formation and total gas density, the profile of the ratio Σstar(R)/Σgas(R) became shallower (flatter) over time. Thus, if metals stayed where they were injected into the ISM from stars, the metallicity gradient would become shallower over time, as some models predict. However, metallicity gradients in FIRE-2 became steeper over time because of the additional effects of (radial) mixing of metals in the ISM. Specifically, the velocity dispersion and net radial advection of gas declined over time, as ISM turbulence decreased and the disk settled, leading to upside-down vertical growth. In FIRE-2, this evolution in metal mixing of gas associated with upside-down growth dominates over inside-out radial growth, causing the metallicity radial gradient of the ISM and of stars at birth to become steeper over time. We argue that this reflects the ISM history of typical MW-mass galaxies and a plausible history of the MW.

We present the first detailed comparison of populations of dwarf galaxy stellar streams in cosmological simulations and the Milky Way. In particular, we compare streams identified around 13 Milky Way analogs in the FIRE-2 simulations to streams observed by the Southern Stellar Stream Spectroscopic Survey (S 5). For an accurate comparison, we produce mock Dark Energy Survey (DES) observations of the FIRE streams and estimate the detectability of their tidal tails and progenitors. The number and stellar mass distributions of detectable stellar streams is consistent between observations and simulations. However, there are discrepancies in the distributions of pericenters and apocenters, with the detectable FIRE streams, on average, forming at larger pericenters (out to >110 kpc) and surviving only at larger apocenters (≳40 kpc) than those observed in the Milky Way. We find that the population of high-stellar-mass dwarf galaxy streams in the Milky Way is incomplete. Interestingly, a large fraction of the FIRE streams would only be detected as intact satellites in DES-like observations, since their tidal tails have too low surface brightness to be detectable. We thus predict a population of yet-undetected tidal tails around Milky Way satellites, as well as a population of fully undetected low-surface-brightness stellar streams, and estimate their detectability with the Rubin Observatory. Finally, we discuss the causes and implications of the discrepancies between the stream populations in FIRE and the Milky Way, and explore future avenues for tests of satellite disruption in cosmological simulations.

We analyze the first cosmological baryonic zoom-in simulations of galaxies in dissipative self-interacting dark matter (dSIDM). The simulations utilize the FIRE-2 galaxy formation physics with the inclusion of dissipative dark matter self-interactions modeled as a constant fractional energy dissipation (f diss = 0.75). In this paper, we examine the properties of dwarf galaxies with M * ∼ 105–109 M ⊙ in both isolation and within Milky Way–mass hosts. For isolated dwarfs, we find more compact galaxy sizes and promotion of disk formation in dSIDM with (σ/m) ≤ 1 cm2 g−1. On the contrary, models with (σ/m) = 10 cm2 g−1 produce puffier stellar distributions that are in tension with the observed size–mass relation. In addition, owing to the steeper central density profiles, the subkiloparsec circular velocities of isolated dwarfs when (σ/m) ≥ 0.1 cm2 g−1 are enhanced by about a factor of 2, which are still consistent with the kinematic measurements of Local Group dwarfs but in tension with the H i rotation curves of more massive field dwarfs. Meanwhile, for satellites of Milky Way–mass hosts, the median circular velocity profiles are marginally affected by dSIDM physics, but dSIDM may help promote the structural diversity of dwarf satellites. The number of satellites is slightly enhanced in dSIDM, but the differences are small compared with the large host-to-host variations. In conclusion, the dSIDM models with (σ/m) ≳ 0.1 cm2 g−1, f diss = 0.75 are in tension in massive dwarfs (M halo ∼ 1011 M ⊙) due to circular velocity constraints. However, models with lower effective cross sections (at this halo mass/velocity scale) are still viable and can produce nontrivial observable signatures.

The third data release (DR3) of Gaia has provided a fivefold increase in the number of radial velocity measurements of stars, as well as a stark improvement in parallax and proper motion measurements. To help with studies that seek to test models and interpret Gaia DR3, we present nine Gaia synthetic surveys, based on three solar positions in three Milky Way-mass galaxies of the Latte suite of the Fire-2 cosmological simulations. These synthetic surveys match the selection function, radial velocity measurements, and photometry of Gaia DR3, adapting the code base Ananke, previously used to match the Gaia DR2 release by Sanderson et al. The synthetic surveys are publicly available and can be found at http://ananke.hub.yt/. Similarly to the previous release of Ananke, these surveys are based on cosmological simulations and thus are able to model nonequilibrium dynamical effects, making them a useful tool in testing and interpreting Gaia DR3.