Search

We characterize mass, momentum, energy, and metal outflow rates of multiphase galactic winds in a suite of FIRE-2 cosmological 'zoom-in' simulations from the Feedback in Realistic Environments (FIRE) project. We analyse simulations of low-mass dwarfs, intermediate-mass dwarfs, Milky Way-mass haloes, and high-redshift massive haloes. Consistent with previous work, we find that dwarfs eject about 100 times more gas from their interstellar medium (ISM) than they form in stars, while this mass 'loading factor' drops below one in massive galaxies. Most of the mass is carried by the hot phase (>105 K) in massive haloes and the warm phase (103-105 K) in dwarfs; cold outflows (

We explore the effect of including progenitor mass- and metallicity-dependent yields, supernova rates and energetics on variations in elemental abundance ratios (particularly [α/Fe]) in dwarf galaxies. To understand how the scatter and overall trends in [α/Fe] are affected by including variable metal yields from a discretely sampled initial mass function, we run FIRE simulations of a dwarf galaxy (M★(z = 0) ∼ 106 M☉) using nucleosynthetic yields from the NuGrid data base that depend on the stellar progenitor mass and metallicity. While NuGrid exhibits lower aggregate α-element production than default FIRE yields, we find that its explicit mass dependence, even when including turbulent metal diffusion, substantially widens the intrinsic scatter in the simulated [Fe/H]-[α/Fe] – a phenomenon visible in some observations of dwarf galaxies.

Recent searches for the hosts of z ∼4 damped Ly α absorbers (DLAs) have detected bright galaxies at distances of tens of kpc from the DLA. Using the FIRE-2 cosmological zoom simulations, we argue that these relatively large distances are due to a predominantly cool and neutral inner circumgalactic medium (CGM) surrounding high-redshift galaxies. The inner CGM is cool because of the short cooling time of hot gas in ≤ 1012, M\odot haloes, which implies that accretion and feedback energy are radiated quickly, while it is neutral due to high volume densities and column densities at high redshift that shield cool gas from photoionization. Our analysis predicts large DLA covering factors (≥ 50 per cent) out to impact parameters ∼0.3[(1 + z)/5]3/2Rvir from the central galaxies at z ≥ 1, equivalent to a proper distance of ∼21, M121/3 (1+z)/5)1/2\, kpc (Rvir and M12 are the halo virial radius and mass in units of 1012, M, respectively). This implies that DLA covering factors at z ∼4 may be comparable to unity out to a distance ∼10 times larger than stellar half-mass radii. A predominantly neutral inner CGM in the early universe suggests that its mass and metallicity can be directly constrained by absorption surveys, without resorting to the large ionization corrections as required for ionized CGM.

We present a suite of baryonic cosmological zoom-in simulations of self-interacting dark matter (SIDM) haloes within the 'Feedback In Realistic Environment' (FIRE) project. The three simulated haloes have virial masses of ∼ 1012 M⊙ at z = 0, and we study velocity-independent self-interaction cross sections of 1 and 10 cm2, g-1. We study star formation rates and the shape of dark matter density profiles of the parent haloes in both cold dark matter (CDM) and SIDM models. Galaxies formed in the SIDM haloes have higher star formation rates at z ≤ 1, resulting in more massive galaxies compared to the CDM simulations. While both CDM and SIDM simulations show diverse shape of the dark matter density profiles, the SIDM haloes can reach higher and more steep central densities within few kpcs compared to the CDM haloes. We identify a correlation between the build-up of the stars within the half-mass radii of the galaxies and the growth in the central dark matter densities. The thermalization process in the SIDM haloes is enhanced in the presence of a dense stellar component. Hence, SIDM haloes with highly concentrated baryonic profiles are predicted to have higher central dark matter densities than the CDM haloes. Overall, the SIDM haloes are more responsive to the presence of a massive baryonic distribution than their CDM counterparts.

We present the first set of cosmological baryonic zoom-in simulations of galaxies including dissipative self-interacting dark matter (dSIDM). These simulations utilize the Feedback In Realistic Environments galaxy formation physics, but allow the dark matter to have dissipative self-interactions analogous to standard model forces, parametrized by the self-interaction crosssection per unit mass, (σ/m), and the dimensionless degree of dissipation, 0 < fdiss < 1. We survey this parameter space, including constant and velocity-dependent cross-sections, and focus on structural and kinematic properties of dwarf galaxies with Mhalo ∼ 1010.11M⊙ and M∗ ∼ 105.8M⊙. Central density profiles (parametrized as ρ ∝ rα) of simulated dwarfs become cuspy when (σ/m)eff ≳ 0.1 cm2 g-1 (and fdiss = 0.5 as fiducial). The power-law slopes asymptote to α ≈ -1.5 in low-mass dwarfs independent of cross-section, which arises from a dark matter 'cooling flow'. Through comparisons with dark matter only simulations, we find the profile in this regime is insensitive to the inclusion of baryons. However, when (σ/m)eff ≤ 0.1 cm2 g-1, baryonic effects can produce cored density profiles comparable to non-dissipative cold dark matter (CDM) runs but at smaller radii. Simulated galaxies with (σ/m) ≳ 10 cm2 g-1 and the fiducial fdiss develop significant coherent rotation of dark matter, accompanied by halo deformation, but this is unlike the well-defined thin 'dark discs' often attributed to baryon-like dSIDM. The density profiles in this high cross-section model exhibit lower normalizations given the onset of halo deformation. For our surveyed dSIDM parameters, halo masses and galaxy stellar masses do not show appreciable difference from CDM, but dark matter kinematics and halo concentrations/shapes can differ.

We investigate thin and thick stellar disc formation in Milky Way-mass galaxies using 12 FIRE-2 cosmological zoom-in simulations. All simulated galaxies experience an early period of bursty star formation that transitions to a late-time steady phase of near-constant star formation. Stars formed during the late-time steady phase have more circular orbits and thin-disc-like morphology at z = 0, while stars born during the bursty phase have more radial orbits and thick-disc structure. The median age of thick-disc stars at z = 0 correlates strongly with this transition time. We also find that galaxies with an earlier transition from bursty to steady star formation have a higher thin-disc fractions at z = 0. Three of our systems have minor mergers with Large Magellanic Cloud-size satellites during the thin-disc phase. These mergers trigger short starbursts but do not destroy the thin disc nor alter broad trends between the star formation transition time and thin/thick-disc properties. If our simulations are representative of the Universe, then stellar archaeological studies of the Milky Way (or M31) provide a window into past star formation modes in the Galaxy. Current age estimates of the Galactic thick disc would suggest that the Milky Way transitioned from bursty to steady phase ~6.5 Gyr ago; prior to that time the Milky Way likely lacked a recognizable thin disc.

We examine the prevalence, longevity, and causes of planes of satellite dwarf galaxies, as observed in the Local Group. We use 14 Milky Way/Andromeda-(MW/M31) mass host galaxies from the Feedback In Realistic Environments-2 simulations. We select the 14 most massive satellites by stellar mass within ≤ 300, kpc of each host and correct for incompleteness from the foreground galactic disc when comparing to the MW. We find that MW-like planes as spatially thin and/or kinematically coherent as observed are uncommon, but they do exist in our simulations. Spatially thin planes occur in 1-2 per cent of snapshots during z = 0-0.2, and kinematically coherent planes occur in 5 per cent of snapshots. These planes are generally transient, surviving for

We use the FIRE-2 cosmological simulations to study the formation of a quasi-static, virial-temperature gas phase in the circumgalactic medium (CGM) at redshifts 0 < z < 5 and how the formation of this virialized phase affects the evolution of galactic disks. We demonstrate that when the halo mass crosses ∼1012 Me, the cooling time of shocked gas in the inner CGM (∼0.1Rvir, where Rvir is the virial radius) exceeds the local free-fall time. The inner CGM then experiences a transition from on average subvirial temperatures (T = Tvir), large pressure fluctuations, and supersonic inflow/outflow velocities to virial temperatures (T ∼ Tvir), uniform pressures, and subsonic velocities. This transition occurs when the outer CGM (∼0.5Rvir) is already subsonic and has a temperature ∼Tvir, indicating that the longer cooling times at large radii allow the outer CGM to virialize at lower halo masses than the inner CGM. This outside-in CGM virialization scenario is in contrast with inside-out scenarios commonly envisioned based on more idealized simulations. We demonstrate that inner CGM virialization coincides with abrupt changes in the central galaxy and its stellar feedback: the galaxy settles into a stable rotating disk, star formation transitions from “bursty” to “steady,” and stellar-driven galaxy-scale outflows are suppressed. Our results thus suggest that CGM virialization is initially associated with the formation of rotation-dominated thin galactic disks, rather than with the quenching of star formation as often assumed.

We explore the origin of stellar metallicity gradients in simulated and observed dwarf galaxies. We use FIRE-2 cosmological baryonic zoom-in simulations of 26 isolated galaxies as well as existing observational data for 10 Local Group dwarf galaxies. Our simulated galaxies have stellar masses between 105.5 and 108.6 MO. Whilst gas-phase metallicty gradients are generally weak in our simulated galaxies, we find that stellar metallicity gradients are common, with central regions tending to be more metal-rich than the outer parts. The strength of the gradient is correlated with galaxy-wide median stellar age, such that galaxies with younger stellar populations have flatter gradients. Stellar metallicty gradients are set by two competing processes: (1) the steady 'puffing' of old, metal-poor stars by feedback-driven potential fluctuations and (2) the accretion of extended, metal-rich gas at late times, which fuels late-time metal-rich star formation. If recent star formation dominates, then extended, metal-rich star formation washes out pre-existing gradients from the 'puffing' process. We use published results from ten Local Group dwarf galaxies to show that a similar relationship between age and stellar metallicity-gradient strength exists among real dwarfs. This suggests that observed stellar metallicity gradients may be driven largely by the baryon/feedback cycle rather than by external environmental effects.

We present an analysis of the R ≲ 1.5 kpc core regions of seven simulated Milky Way-mass galaxies, from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulation suite, for a finely sampled period (Δt = 2.2 Myr) of 22 Myr at z ≈ 0, and compare them with star formation rate (SFR) and gas surface density observations of the Milky Way's Central Molecular Zone (CMZ). Despite not being tuned to reproduce the detailed structure of the CMZ, we find that four of these galaxies are consistent with CMZ observations at some point during this 22 Myr period. The galaxies presented here are not homogeneous in their central structures, roughly dividing into two morphological classes; (a) several of the galaxies have very asymmetric gas and SFR distributions, with intense (compact) starbursts occurring over a period of roughly 10 Myr, and structures on highly eccentric orbits through the CMZ, whereas (b) others have smoother gas and SFR distributions, with only slowly varying SFRs over the period analyzed. In class (a) centers, the orbital motion of gas and star-forming complexes across small apertures (R ≲ 150 pc, analogously |l| < 1° in the CMZ observations) contributes as much to tracers of star formation/dense gas appearing in those apertures, as the internal evolution of those structures does. These asymmetric/bursty galactic centers can simultaneously match CMZ gas and SFR observations, demonstrating that time-varying star formation can explain the CMZ's low star formation efficiency.

We present models of CO(1-0) emission from Milky-Way-mass galaxies at redshift zero in the FIRE-2 cosmological zoom-in simulations. We calculate the molecular abundances by post-processing the simulations with an equilibrium chemistry solver while accounting for the effects of local sources, and determine the emergent CO(1-0) emission using a line radiative transfer code. We find that the results depend strongly on the shielding length assumed, which, in our models, sets the attenuation of the incident UV radiation field. At the resolution of these simulations, commonly used choices for the shielding length, such as the Jeans length, result in CO abundances that are too high at a given H2 abundance. We find that a model with a distribution of shielding lengths, which has a median shielding length of ~3 pc in cold gas (T < 300 K) for both CO and H2, is able to reproduce both the observed CO(1-0) luminosity and inferred CO-to-H2 conversion factor at a given star formation rate compared with observations. We suggest that this short shielding length can be thought of as a subgrid model, which controls the amount of radiation that penetrates giant molecular clouds.

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.

Pressure balance plays a central role in models of the interstellar medium (ISM), but whether and how pressure balance is realized in a realistic multiphase ISM is not yet well understood. We address this question by using a set of FIRE-2 cosmological zoom-in simulations of Milky Way-mass disc galaxies, in which a multiphase ISM is self-consistently shaped by gravity, cooling, and stellar feedback. We analyse how gravity determines the vertical pressure profile as well as how the total ISM pressure is partitioned between different phases and components (thermal, dispersion/turbulence, and bulk flows). We show that, on average and consistent with previous more idealized simulations, the total ISM pressure balances the weight of the overlying gas. Deviations from vertical pressure balance increase with increasing galactocentric radius and with decreasing averaging scale. The different phases are in rough total pressure equilibrium with one another, but with large deviations from thermal pressure equilibrium owing to kinetic support in the cold and warm phases, which dominate the total pressure near the mid-plane. Bulk flows (e.g. inflows and fountains) are important at a few disc scale heights, while thermal pressure from hot gas dominates at larger heights. Overall, the total mid-plane pressure is well-predicted by the weight of the disc gas and we show that it also scales linearly with the star formation rate surface density (ςSFR). These results support the notion that the Kennicutt-Schmidt relation arises because ςSFR and the gas surface density (ςg) are connected via the ISM mid-plane pressure.

We present the escape fraction of hydrogen ionizing photons (fesc) from a sample of 34 high-resolution cosmological zoom-in simulations of galaxies at z ≥ 5 in the Feedback in Realistic Environments project, post-processed with a Monte Carlo radiative transfer code for ionizing radiation. Our sample consists of 8500 haloes in Mvir ∼108-1012, M⊙ (M* ∼104-1010, M⊙) at z = 5-12. We find the sample average fesc increases with halo mass for Mvir ∼108-109.5 M⊙, becomes nearly constant for 109.5-1011 M⊙, and decreases at ≳ 1011 M⊙. Equivalently, fesc increases with stellar mass up to M* ∼ 108 M⊙ and decreases at higher masses. Even applying single-star stellar population synthesis models, we find a moderate fesc ∼0.2 for galaxies at M* ∼ 108 M⊙. Nearly half of the escaped ionizing photons come from stars 1-3 Myr old and the rest from stars 3-10 Myr old. Binaries only have a modest effect, boosting fesc by ∼25-35 per cent and the number of escaped photons by 60-80 per cent. Most leaked ionizing photons are from vigorously star-forming regions that usually contain a feedback-driven kpc-scale superbubble surrounded by a dense shell. The shell is forming stars while accelerated, so new stars formed earlier in the shell are already inside the shell. Young stars in the bubble and near the edge of the shell can fully ionize some low-column-density paths pre-cleared by feedback, allowing a large fraction of their ionizing photons to escape. The decrease of fesc at the high-mass end is due to dust attenuation, while at the low-mass end, fesc decreases owing to inefficient star formation and hence feedback. At fixed mass, fesc tends to increase with redshift. Although the absolute fesc does not fully converge with resolution in our simulations, the mass- A nd redshift-dependence of fesc is likely robust. Our simulations produce sufficient ionizing photons for cosmic reionization.

We present the hitherto largest sample of gas-phase metallicity radial gradients measured at sub-kpc resolution in star-forming galaxies in the redshift range of z ∈ [1.2, 2.3]. These measurements are enabled by the synergy of slitless spectroscopy from the Hubble Space Telescope near-infrared channels and the lensing magnification from foreground galaxy clusters. Our sample consists of 76 galaxies with stellar mass ranging from 107 to 1010 M⊙, an instantaneous star formation rate in the range of [1, 100] M⊙ yr-1, and global metallicity [1/12 ,2] of solar. At a 2σ confidence level, 15/76 galaxies in our sample show negative radial gradients, whereas 7/76 show inverted gradients. Combining ours and all other metallicity gradients obtained at a similar resolution currently available in the literature, we measure a negative mass dependence of Δlog(O/H)/Δ r [dex kpc-1] = (-0.020 ± 0.007) + (-0.016 ± 0.008), with the intrinsic scatter being σ = 0.060 ± 0.006 over 4 orders of magnitude in stellar mass. Our result is consistent with strong feedback, not secular processes, being the primary governor of the chemostructural evolution of star-forming galaxies during the disk mass assembly at cosmic noon. We also find that the intrinsic scatter of metallicity gradients increases with decreasing stellar mass and increasing specific star formation rate. This increase in the intrinsic scatter is likely caused by the combined effect of cold-mode gas accretion and merger-induced starbursts, with the latter more predominant in the dwarf mass regime of M∗ ≤ 109 M⊙.

We present the first measurement of the lifetimes of giant molecular clouds (GMCs) in cosmological simulations at z = 0, using the Latte suite of FIRE-2 simulations of Milky Way (MW) mass galaxies. We track GMCs with total gas mass ≳105 M⊙ at high spatial (∼1 pc), mass (7100 M⊙), and temporal (1 Myr) resolution. Our simulated GMCs are consistent with the distribution of masses for massive GMCs in the MW and nearby galaxies. We find GMC lifetimes of 5-7 Myr, or 1-2 freefall times, on average, with less than 2 per cent of clouds living longer than 20 Myr. We find decreasing GMC lifetimes with increasing virial parameter, and weakly increasing GMC lifetimes with galactocentric radius, implying that environment affects the evolutionary cycle of GMCs. However, our GMC lifetimes show no systematic dependence on GMC mass or amount of star formation. These results are broadly consistent with inferences from the literature and provide an initial investigation into ultimately understanding the physical processes that govern GMC lifetimes in a cosmological setting.

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.

Aims. The goal of this study is to present the development of a machine learning based approach that utilizes phase space alone to separate the Gaia DR2 stars into two categories: those accreted onto the Milky Way from those that are in situ. Traditional selection methods that have been used to identify accreted stars typically rely on full 3D velocity, metallicity information, or both, which significantly reduces the number of classifiable stars. The approach advocated here is applicable to a much larger portion of Gaia DR2. Methods. A method known as "transfer learning" is shown to be effective through extensive testing on a set of mock Gaia catalogs that are based on the FIRE cosmological zoom-in hydrodynamic simulations of Milky Way-mass galaxies. The machine is first trained on simulated data using only 5D kinematics as inputs and is then further trained on a cross-matched Gaia/RAVE data set, which improves sensitivity to properties of the real Milky Way. Results. The result is a catalog that identifies ∼767 000 accreted stars within Gaia DR2. This catalog can yield empirical insights into the merger history of the Milky Way and could be used to infer properties of the dark matter distribution.

We report the formation of bound star clusters in a sample of high-resolution cosmological zoom-in simulations of z ≥ 5 galaxies from the Feedback In Realistic Environments project. We find that bound clusters preferentially form in high-pressure clouds with gas surface densities over 104M⊙ pc-2, where the cloud-scale star formation efficiency is near unity and young stars born in these regions are gravitationally bound at birth. These high-pressure clouds are compressed by feedback-driven winds and/or collisions of smaller clouds/gas streams in highly gas-rich, turbulent environments. The newly formed clusters follow a power-law mass function of dN/dM ~ M-2. The cluster formation efficiency is similar across galaxies with stellarmasses of~107-1010M⊙ at z ≥ 5. The age spread of cluster stars is typically a few Myr and increases with cluster mass. The metallicity dispersion of cluster members is ~0.08 dex in [Z/H] and does not depend on clustermass significantly. Our findings support the scenario that present-day old globular clusters (GCs) were formed during relatively normal star formation in high-redshift galaxies. Simulations with a stricter/looser star formation model form a factor of a few more/fewer bound clusters per stellar mass formed, while the shape of the mass function is unchanged. Simulations with a lower local star formation efficiency form more stars in bound clusters. The simulated clusters are larger than observed GCs due to finite resolution. Our simulations are among the first cosmological simulations that form bound clusters self-consistently in a wide range of high-redshift galaxies.

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.

We study the spatially resolved (sub-kpc) gas velocity dispersion (σ)-star formation rate (SFR) relation in the FIRE-2 (Feedback in Realistic Environments) cosmological simulations. We specifically focus on Milky Way-mass disc galaxies at late times (z ≈ 0). In agreement with observations, we find a relatively flat relationship, with σ ≈ 15-30 km s−1 in neutral gas across 3 dex in SFRs. We show that higher dense gas fractions (ratios of dense gas to neutral gas) and SFRs are correlated at constant σ. Similarly, lower gas fractions (ratios of gas to stellar mass) are correlated with higher σ at constant SFR. The limits of the σ-ΣSFR relation correspond to the onset of strong outflows. We see evidence of 'on-off' cycles of star formation in the simulations, corresponding to feedback injection time-scales of 10-100 Myr, where SFRs oscillate about equilibrium SFR predictions. Finally, SFRs and velocity dispersions in the simulations agree well with feedback-regulated and marginally stable gas disc (Toomre's Q = 1) model predictions, and the simulation data effectively rule out models assuming that gas turns into stars at (low) constant efficiency (i.e. 1 per cent per free-fall time). And although the simulation data do not entirely exclude gas accretion/gravitationally powered turbulence as a driver of σ, it appears to be subdominant to stellar feedback in the simulated galaxy discs at z ≈ 0.

We study stellar-halo formation using six Milky-Way-mass galaxies in FIRE-2 cosmological zoom simulations. We find that 5-40 per cent of the outer (50-300 kpc) stellar halo in each system consists of in-situ stars that were born in outflows from the main galaxy. Outflow stars originate from gas accelerated by superbubble winds, which can be compressed, cool, and form co-moving stars. The majority of these stars remain bound to the halo and fall back with orbital properties similar to the rest of the stellar halo at z = 0. In the outer halo, outflow stars are more spatially homogeneous, metal-rich, and alpha-element-enhanced than the accreted stellar halo. At the solar location, up to ∼10 per cent of our kinematically identified halo stars were born in outflows; the fraction rises to as high as ∼40 per cent for the most metal-rich local halo stars ([Fe/H] >-0.5). Such stars can be retrograde and create features similar to the recently discoveredMilkyWay 'Splash' in phase space.We conclude that theMilkyWay stellar halo could contain local counterparts to stars that are observed to form in molecular outflows in distant galaxies. Searches for such a population may provide a new, near-field approach to constraining feedback and outflow physics. A stellar halo contribution from outflows is a phase-reversal of the classic halo formation scenario of Eggen, Lynden-Bell & Sandange, who suggested that halo stars formed in rapidly infalling gas clouds. Stellar outflows may be observable in direct imaging of external galaxies and could provide a source for metal-rich, extreme-velocity stars in the Milky Way.

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.

Radiative feedback (RFB) from stars plays a key role in galaxies, but remains poorly understood. We explore this using high-resolution, multifrequency radiation-hydrodynamics (RHD) simulations from the Feedback In Realistic Environments (FIRE) project. We study ultrafaint dwarf through Milky Way mass scales, including H+He photoionization; photoelectric, Lyman Werner, Compton, and dust heating; and single+multiple scattering radiation pressure (RP). We compare distinct numerical algorithms: ray-based LEBRON (exact when optically thin) and moments-based M1 (exact when optically thick). The most important RFB channels on galaxy scales are photoionization heating and single-scattering RP: in all galaxies, most ionizing/far-UV luminosity (∼1/2 of lifetime-integrated bolometric) is absorbed. In dwarfs, the most important effect is photoionization heating from the UV background suppressing accretion. In MW-mass galaxies, metagalactic backgrounds have negligible effects; but local photoionization and single-scattering RP contribute to regulating the galactic star formation efficiency and lowering central densities. Without some RFB (or other 'rapid' FB), resolved GMCs convert too-efficiently into stars, making galaxies dominated by hyperdense, bound star clusters. This makes star formation more violent and 'bursty' when SNe explode in these hyperclustered objects: thus, including RFB 'smoothes' SFHs. These conclusions are robust to RHD methods, but M1 produces somewhat stronger effects. Like in previous FIRE simulations, IR multiple-scattering is rare (negligible in dwarfs, ∼ 10 per cent of RP in massive galaxies): absorption occurs primarily in 'normal' GMCs with AV ∼ 1.

We present the star formation history (SFH) of the isolated (D ∼ 970 kpc) Local Group dwarf galaxy Wolf-Lundmark-Melotte (WLM) measured from colour-magnitude diagrams (CMDs) constructed from deep Hubble Space Telescope imaging. Our observations include a central (0.5 rh) and outer field (0.7 rh) that reach below the oldest main-sequence turn-off. WLM has no early dominant episode of star formation: 20 per cent of its stellar mass formed by ∼12.5 Gyr ago (z ∼ 5). It also has an SFR that rises to the present with 50 per cent of the stellar mass within the most recent 5 Gyr (z < 0.7). There is evidence of a strong age gradient: the mean age of the outer field is 5 Gyr older than the inner field despite being only 0.4 kpc apart. Some models suggest such steep gradients are associated with strong stellar feedback and dark-matter core creation. The SFHs of real isolated dwarf galaxies and those from the Feedback in Realistic Environment suite are in good agreement for M*(z = 0) ∼ 107-109M☉, but in worse agreement at lower masses (M*(z = 0) ∼ 105-107 M☉). These differences may be explainable by systematics in the models (e.g. reionization model) and/or observations (HST field placement). We suggest that a coordinated effort to get deep CMDs between HST/JWST (crowded central fields) and WFIRST (wide-area halo coverage) is the optimal path for measuring global SFHs of isolated dwarf galaxies.

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 explore the radial variation of star formation histories (SFHs) in dwarf galaxies simulated with Feedback In Realistic Environments (FIRE) physics. The sample contains 26 field dwarf galaxies with Mstar = 105–109 M. We find 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 SFHs of the merging galaxy. In galaxies without significant mergers, feedback pushes stars to the outskirts. The strength of the age gradient is determined by the subsequent evolution of the galaxy. Galaxies with weak age gradients constantly grow to z = 0, meaning that young star formation occurs at a similar radius to which older stars are heated to. In contrast, galaxies with strong age gradients tend to maintain a constant half-mass radius over time. 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 SFH. Central fields can be biased young by 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 SFH.

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.

One important result from recent large integral field spectrograph (IFS) surveys is that the intrinsic velocity dispersion of galaxies traced by star-forming gas increases with redshift. Massive, rotation-dominated discs are already in place at z ∼ 2, but they are dynamically hotter than spiral galaxies in the local Universe. Although several plausible mechanisms for this elevated velocity dispersion (e.g. star formation feedback, elevated gas supply, or more frequent galaxy interactions) have been proposed, the fundamental driver of the velocity dispersion enhancement at high redshift remains unclear. We investigate the origin of this kinematic evolution using a suite of cosmological simulations from the FIRE (Feedback In Realistic Environments) project. Although IFS surveys generally cover a wider range of stellar masses than in these simulations, the simulated galaxies show trends between intrinsic velocity dispersion (σ intr ), SFR, and z in agreement with observations. In both observations and simulations, galaxies on the star-forming main sequence have median σ intr values that increase from z ∼ 0 to z ∼ 1–1.5, but this increasing trend is less evident at higher redshift. In the FIRE simulations, σ intr can vary significantly on time-scales of 100 Myr. These variations closely mirror the time evolution of the SFR and gas inflow rate (M gas ). By cross-correlating pairs of σ intr, M gas, and SFR, we show that increased gas inflow leads to subsequent enhanced star formation, and enhancements in σ intr tend to temporally coincide with increases in M gas and SFR.

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 investigate the merger histories of isolated dwarf galaxies based on a suite of 15 highresolution cosmological zoom-in simulations, all with masses of Mhalo ≈ 1010M⊙ (and M* ~ 105 - 107M⊙) at z= 0, from the Feedback in Realistic Environments project. The stellar populations of these dwarf galaxies at z= 0 are formed essentially entirely 'in situ': over 90 per cent of the stellar mass is formed in the main progenitor in all but two cases, and all 15 of the galaxies have > 70 per cent of their stellar mass formed in situ. Virtually all galaxy mergers occur prior to z~ 3, meaning that accreted stellar populations are ancient. On average, our simulated dwarfs undergo five galaxy mergers in their lifetimes, with typical pre-merger galaxy mass ratios that are less than 1:10. This merger frequency is generally comparable to what has been found in dissipationless simulations when coupled with abundance matching. Two of the simulated dwarfs have a luminous satellite companion at z= 0. These ultra-faint dwarfs lie at or below current detectability thresholds but are intriguing targets for nextgeneration facilities. The small contribution of accreted stars makes it extremely difficult to discern the effects of mergers in the vast majority of dwarfs either photometrically or using resolved-star colour-magnitude diagrams (CMDs). The important implication for near-field cosmology is that star formation histories (SFHs) of comparably massive galaxies derived from resolved CMDs should trace the build-up of stellar mass in one main system across cosmic time as opposed to reflecting the contributions of many individual SFHs of merged dwarfs.

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 implementation of mechanical feedback from supernovae (SNe) and stellar mass loss in galaxy simulations, within the Feedback In Realistic Environments (FIRE) project. We present the FIRE-2 algorithm for coupling mechanical feedback, which can be applied to any hydrodynamics method (e.g. fixed-grid, moving-mesh, and mesh-less methods), and black hole as well as stellar feedback. This algorithm ensures manifest conservation of mass, energy, and momentum, and avoids imprinting 'preferred directions' on the ejecta. We show that it is critical to incorporate both momentum and thermal energy of mechanical ejecta in a self-consistent manner, accounting for SNe cooling radii when they are not resolved. Using idealized simulations of single SN explosions, we show that the FIRE-2 algorithm, independent of resolution, reproduces converged solutions in both energy and momentum. In contrast, common 'fully thermal' (energy-dump) or 'fully kinetic' (particle-kicking) schemes in the literature depend strongly on resolution: when applied at mass resolution ≳100M⊙, they diverge by orders of magnitude from the converged solution. In galaxy-formation simulations, this divergence leads to orders-of-magnitude differences in galaxy properties, unless those models are adjusted in a resolution-dependent way. We show that all models that individually time-resolve SNe converge to the FIRE-2 solution at sufficiently high resolution (< 100M⊙). However, in both idealized single-SN simulations and cosmological galaxy-formation simulations, the FIRE-2 algorithm converges much faster than other sub-grid models without re-tuning parameters.

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.

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 stellar metallicity distribution functions (MDFs), including Fe and a-element abundances, in dwarf galaxies from the Feedback in Realistic Environment (FIRE) project. We examine both isolated dwarf galaxies and those that are satellites of a MilkyWay-mass galaxy. In particular, we study the effects of including a sub-grid turbulent model for the diffusion of metals in gas. Simulations that include diffusion have narrower MDFs and abundance ratio distributions, because diffusion drives individual gas and star particles towards the average metallicity. This effect provides significantly better agreement with observed abundance distributions in dwarf galaxies in the Local Group, including small intrinsic scatter in [α/Fe] versus [Fe/H] of ≲0.1 dex. This small intrinsic scatter arises in our simulations because the interstellar medium in dwarf galaxies is well mixed at nearly all cosmic times, such that stars that form at a given time have similar abundances to ≲0.1 dex. Thus, most of the scatter in abundances at z = 0 arises from redshift evolution and not from instantaneous scatter in the ISM. We find similar MDF widths and intrinsic scatter for satellite and isolated dwarf galaxies, which suggests that environmental effects play a minor role compared with internal chemical evolution in our simulations. Overall, with the inclusion of metal diffusion, our simulations reproduce abundance distribution widths of observed low-mass galaxies, enabling detailed studies of chemical evolution in galaxy formation.

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.

We introduce massive black holes (BHs) in the Feedback In Realistic Environments (FIRE) project and perform high-resolution cosmological hydrodynamic simulations of quasar-mass haloes [Mhalo(z = 2) ≈ 1012.5M⊙] down to z = 1. These simulations model stellar feedback by supernovae, stellar winds and radiation, and BH growth using a gravitational torque-based prescription tied to the resolved properties of galactic nuclei.We do not include BH feedback. We show that early BH growth occurs through short (≲1Myr) accretion episodes that can reach or even exceed the Eddington rate. In this regime, BH growth is limited by bursty stellar feedback continuously evacuating gas from galactic nuclei, and BHs remain undermassive in low-mass galaxies relative to the local MBH-Mbulgerelation. BH growth is more efficient at later times, when the nuclear stellar potential retains a significant gas reservoir, star formation becomes less bursty and galaxies settle into a more ordered state. BHs rapidly converge on to the observed scaling relations when the host reaches Mbulge ~1010M⊙. We show that resolving the effects of stellar feedback on the gas supply in the inner ~100 pc of galaxies is necessary to accurately capture the growth of central BHs. Our simulations imply that bursty stellar feedback has important implications for BH-galaxy relations, AGN demographics and time variability, the formation of early quasars and massive BH mergers.

We present a suite of 15 cosmological zoom-in simulations of isolated dark matter haloes, all with masses of Mhalo ≈ 1010M⊙ at z = 0, in order to understand the relationship among halo assembly, galaxy formation and feedback's effects on the central density structure in dwarf galaxies. These simulations are part of the Feedback in Realistic Environments (FIRE) project and are performed at extremely high resolution (mbaryon = 500M⊙, mdm = 2500M⊙). The resultant galaxies have stellar masses that are consistent with rough abundance matching estimates, coinciding with the faintest galaxies that can be seen beyond the virial radius of the Milky Way (M*/M⊙ ≈ 105 - 107). This non-negligible spread in stellar mass at z = 0 in haloes within a narrow range of virial masses is strongly correlated with central halo density or maximum circular velocity Vmax, both of which are tightly linked to halo formation time. Much of this dependence of M* on a second parameter (beyond Mhalo) is a direct consequence of the Mhalo ~ 1010M⊙ mass scale coinciding with the threshold for strong reionization suppression: the densest, earliest-forming haloes remain above the UV-suppression scale throughout their histories while late-forming systems fall below the UV-suppression scale over longer periods and form fewer stars as a result. In fact, the latest-forming, lowest-concentration halo in our suite fails to form any stars. Haloes that form galaxies with M ≳ 2 × 106 M⊙ have reduced central densities relative to dark-matter-only simulations, and the radial extent of the density modifications is well-approximated by the galaxy half-mass radius r1/2. Lower-mass galaxies do not modify their host dark matter haloes at the mass scale studied here. This apparent stellar mass threshold of M ≈ 2 × 106-2 × 10-4 Mhalo is broadly consistent with previous work and provides a testable prediction of FIRE feedback models in Λcold dark matter.

With the upcoming release of the Gaia catalog and the many multiplexed spectroscopic surveys on the horizon, we are rapidly moving into a new data-driven era in the study of the Milky Way's stellar halo. When combined, these data sets will give us a many-dimensional view of stars in accreted structures in the halo that includes both dynamical information about their orbits and chemical information about their formation histories. Using simulated data from the state-of-the-art Latte simulations of Milky-Way-like galaxies, which include hydrodynamics, feedback, and chemical evolution in a cosmological setting, we demonstrate that while dynamical information alone can be used to constrain models of the Galactic mass distribution in the halo, including the extra dimensions provided by chemical abundances can improve these constraints as well as assist in untangling different accreted components.

We use the first Gaia data release, combined with the RAVE and APOGEE spectroscopic surveys, to investigate the origin of halo stars ≲3 kpc within kpc from the Sun. We identify halo stars kinematically as moving at a relative speed of at least 220 km s-1 with respect to the local standard of rest. These stars are generally less metal-rich than the disk, but surprisingly, half of our halo sample is comprised of stars with [Fe H] >-1. The orbital directions of these metal-rich halo stars are preferentially aligned with the disk rotation, in sharp contrast with the intrinsically isotropic orbital distribution of the metal-poor halo stars. We find similar properties in the Latte cosmological zoom-in simulation of a Milky Way-like galaxy from the FIRE project. In Latte, metal-rich halo stars formed primarily inside of the solar circle, whereas lower-metallicity halo stars preferentially formed at larger distances (extending beyond the virial radius). This suggests that metal-rich halo stars in the solar neighborhood actually formed in situ within the Galactic disk, rather than having been accreted from satellite systems. These stars, currently on halo-like orbits, therefore have likely undergone substantial radial migration/heating.

We study the structure, age and metallicity gradients, and dynamical evolution using a cosmological zoom-in simulation of a Milky Way-mass galaxy from the Feedback in Realistic Environments project. In the simulation, stars older than 6 Gyr were formed in a chaotic, bursty mode and have the largest vertical scaleheights (1.5-2.5 kpc) by z = 0, while stars younger than 6 Gyr were formed in a relatively calm, stable disc. The vertical scaleheight increases with stellar age at all radii, because (1) stars that formed earlier were thicker 'at birth', and (2) stars were kinematically heated to an even thicker distribution after formation. Stars of the same age are thicker in the outer disc than in the inner disc (flaring). These lead to positive vertical age gradients and negative radial age gradients. The radial metallicity gradient is negative at the mid-plane, flattens at larger disc height |Z|, and turns positive above |Z| ~ 1.5 kpc. The vertical metallicity gradient is negative at all radii, but is steeper at smaller radii. These trends broadly agree with observations in the Milky Way and can be naturally understood from the age gradients. The vertical stellar density profile can be well described by two components, with scaleheights 200-500 pc and 1-1.5 kpc, respectively. The thick component is a mix of stars older than 4 Gyr, which formed through a combination of several mechanisms. Our results also demonstrate that it is possible to form a thin disc in cosmological simulations even with a strong stellar feedback.