Search

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.

In hierarchical structure formation, metal-poor stars in and around the Milky Way (MW) originate primarily from mergers of lower mass galaxies. A common expectation is therefore that metal-poor stars should have isotropic, dispersion-dominated orbits that do not correlate strongly with the MW disc. However, recent observations of stars in the MW show that metal-poor ([Fe/H] ≲ -2) stars are preferentially on prograde orbits with respect to the disc. Using the Feedback In Realistic Environments 2 (FIRE-2) suite of cosmological zoom-in simulations of MW/M31-mass galaxies, we investigate the prevalence and origin of prograde metal-poor stars. Almost all (11 of 12) of our simulations have metal-poor stars on preferentially prograde orbits today and throughout most of their history: we thus predict that this is a generic feature of MW/M31-mass galaxies. The typical prograde-to-retrograde ratio is ~2:1, which depends weakly on stellar metallicity at [Fe/H] ~ -1. These trends predicted by our simulations agree well with MW observations. Prograde metal-poor stars originate largely from a single Large/Small Magellanic Cloud (LMC/SMC)-mass gas-rich merger 7-12.5 Gyr ago, which deposited existing metal-poor stars and significant gas on an orbital vector that sparked the formation of and/or shaped the orientation of a long-lived stellar disc, giving rise to a prograde bias for all low-metallicity stars. We find subdominant contributions from in situ stars formed in the host galaxy before this merger, and in some cases, additional massive mergers. We find few clear correlations between any properties of our MW/M31-mass galaxies at z = 0 and the degree of this prograde bias as a result of diverse merger scenarios.

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.

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

We 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 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.

Surveys of the Milky Way (MW) and M31 enable detailed studies of stellar populations across ages and metallicities, with the goal of reconstructing formation histories across cosmic time. These surveys motivate key questions for galactic archaeology in a cosmological context: When did the main progenitor of an MW/M31-mass galaxy form, and what were the galactic building blocks that formed it? We investigate the formation times and progenitor galaxies of MW/M31-mass galaxies using the Feedback In Realistic Environments-2 cosmological simulations, including six isolated MW/M31-mass galaxies and six galaxies in Local Group (LG)-like pairs at z = 0. We examine main progenitor 'formation' based on two metrics: (1) transition from primarily ex-situ to in-situ stellar mass growth and (2) mass dominance compared to other progenitors. We find that the main progenitor of an MW/M31-mass galaxy emerged typically at z ∼3-4 ($11.6\!\!-\!\!12.2\, \rm {Gyr}$ ago), while stars in the bulge region (inner 2 kpc) at z = 0 formed primarily in a single main progenitor at z ≲ 5 (${\lesssim} \!12.6\, \rm {Gyr}$ ago). Compared with isolated hosts, the main progenitors of LG-like paired hosts emerged significantly earlier (Δz ∼2, $\Delta t\!\sim \!1.6\, \rm {Gyr}$), with ∼4× higher stellar mass at all z ≳ 4 (${\gtrsim} \!12.2\, \rm {Gyr}$ ago). This highlights the importance of environment in MW/M31-mass galaxy formation, especially at early times. On average, about 100 galaxies with $\rm {\it{ M}}_\rm {star}\!\gtrsim \!10^5\, \rm {M}_\odot$ went into building a typical MW/M31-mass system. Thus, surviving satellites represent a highly incomplete census (by ∼5×) of the progenitor population.

Recent zoom-in cosmological simulations have shown that stellar feedback can flatten the inner density profile of the dark matter halo in low-mass galaxies. A correlation between the stellar/gas velocity dispersion (σ star, σ gas) and the specific star formation rate (sSFR) is predicted as an observational test of the role of stellar feedback in re-shaping the dark matter density profile. In this work we test the validity of this prediction by studying a sample of star-forming galaxies at 0.6 < z < 1.0 from the LEGA-C survey, which provides high signal-to-noise measurements of stellar and gas kinematics. We find that a weak but significant correlation between σ star (and σ gas) and sSFR indeed exists for galaxies in the lowest mass bin (M ∗ ∼ 1010 M o˙). This correlation, albeit with a ∼35% scatter, holds for different tracers of star formation, and becomes stronger with redshift. This result generally agrees with the picture that at higher redshifts star formation rate was generally higher, and galaxies at M ∗ ≲ 1010 M o˙ have not yet settled into a disk. As a consequence, they have shallower gravitational potentials more easily perturbed by stellar feedback. The observed correlation between σ star (and σ gas) and sSFR supports the scenario predicted by cosmological simulations, in which feedback-driven outflows cause fluctuations in the gravitation potential that flatten the density profiles of low-mass galaxies.

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 analyse the different fates of the circumgalactic medium (CGM) in FIRE-2 cosmological simulations, focusing on the redshifts z = 0.25 and 2 representative of recent surveys. Our analysis includes 21 zoom-in simulations covering the halo mass range Mh(z = 0) ∼ 1010-1012 M⊙. We analyse both where the gas ends up after first leaving the CGM (its 'proximate' fate) and its location at z = 0 (its 'ultimate' fate). Of the CGM at z = 2, about half is found in the ISM or stars of the central galaxy by z = 0 in Mh(z = 2) ∼ 5 × 1011 M⊙ haloes, but most of the CGM in lower mass haloes ends up in the intergalactic medium (IGM). This is so even though most of the CGM in Mh(z = 2) ∼ 5 × 1010 M⊙ haloes first accretes on to the central galaxy before being ejected into the IGM. On the other hand, most of the CGM mass at z = 0.25 remains in the CGM by z = 0 at all halo masses analysed. Of the CGM gas that subsequently accretes on to the central galaxy in the progenitors of Mh(z = 0) ∼ 1012 M⊙ haloes, most of it is cool (T ∼ 104 K) at z = 2 but hot (∼Tvir) at z ∼ 0.25, consistent with the expected transition from cold mode to hot mode accretion. Despite the transition in accretion mode, at both z =0.25 and 2 ≥ 80 per cent of the cool gas inMh ≥ 1011 M⊙ haloes will accrete on to a galaxy. We find that the metallicity of CGM gas is typically a poor predictor of both its proximate and ultimate fates. This is because there is in general little correlation between the origin of CGM gas and its fate owing to substantial mixing while in the CGM.

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 (

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

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.

We investigate the lowest-mass quiescent galaxies known to exist in isolated environments (M* = 109.0-9.5 M⊙; 1.5 Mpc from a more massive galaxy). This population may represent the lowest stellar mass galaxies in which internal feedback quenches galaxy-wide star formation. We present a Keck/Echelle Spectrograph and Imager longslit spectroscopy for 27 isolated galaxies in this regime (20 quiescent galaxies and 7 star-forming galaxies). We measure emission line strengths as a function of radius and place galaxies on the Baldwin Phillips Terlevich (BPT) diagram. Remarkably, 16 of 20 quiescent galaxies in our sample host central active galactic nucleus (AGN)- like line ratios. Only five of these quiescent galaxies were identified as AGN-like in the Sloan Digital Sky Survey due to a lower spatial resolution and signal-to-noise ratio. We find that many of the quiescent galaxies in our sample have spatially extended emission across the non-star-forming regions of BPT-space. While quenched galaxies in denser environments in this mass range often show no evidence for AGN activity, a significant fraction of quiescent galaxies in isolation host AGNs despite their overall passive appearances.

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 introduce the OSIRIS Lens-Amplified Survey (OLAS), a kinematic survey of gravitationally lensed galaxies at cosmic noon taken with Keck adaptive optics. In this paper, we present spatially resolved spectroscopy and nebular emission kinematic maps for 17 star-forming galaxies with stellar masses 8 < log(M ∗/M o) < 9.8 and redshifts 1.2 < z < 2.3. OLAS is designed to probe the stellar mass (M ∗) and specific star formation rate (sSFR) range where simulations suggest that stellar feedback is most effective at driving gaseous outflows that create galaxy-wide potential fluctuations, which can generate dark matter cores. We compare our kinematic data with the trend among sSFR, M ∗, and Hα velocity dispersion, σ, from the Feedback In Realistic Environments (FIRE) simulations. Our observations reveal a correlation between sSFR and σ at fixed M ∗ that is similar to the trend predicted by simulations: feedback from star formation drives star-forming gas and newly formed stars into more dispersion-dominated orbits. The observed magnitude of this effect is in good agreement with the FIRE simulations, in which feedback alters the central density profiles of low-mass galaxies, converting dark matter cusps into cores over time. Our data support the scenario that stellar feedback drives gaseous outflows and potential fluctuations, which in turn drive dark matter core formation in dwarf galaxies.

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 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 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.

In low-mass galaxies, stellar feedback can drive gas outflows that generate non-equilibrium fluctuations in the gravitational potential. Using cosmological zoom-in baryonic simulations from the Feedback in Realistic Environments project, we investigate how these fluctuations affect stellar kinematics and the reliability of Jeans dynamical modeling in low-mass galaxies. We find that stellar velocity dispersion and anisotropy profiles fluctuate significantly over the course of galaxies' starburst cycles. We therefore predict an observable correlation between star formation rate and stellar kinematics: dwarf galaxies with higher recent star formation rates should have systemically higher stellar velocity dispersions. This prediction provides an observational test of the role of stellar feedback in regulating both stellar and dark-matter densities in dwarf galaxies. We find that Jeans modeling, which treats galaxies as virialized systems in dynamical equilibrium, overestimates a galaxy's dynamical mass during periods of post-starburst gas outflow and underestimates it during periods of net inflow. Short-timescale potential fluctuations lead to typical errors of ∼20% in dynamical mass estimates, even if full three-dimensional stellar kinematics - including the orbital anisotropy - are known exactly. When orbital anisotropy is not known a priori, typical mass errors arising from non-equilibrium fluctuations in the potential are larger than those arising from the mass-anisotropy degeneracy. However, Jeans modeling alone cannot reliably constrain the orbital anisotropy, and problematically, it often favors anisotropy models that do not reflect the true profile. If galaxies completely lose their gas and cease forming stars, fluctuations in the potential subside, and Jeans modeling becomes much more reliable.

We examine the effects of stellar feedback and bursty star formation on low-mass galaxies (Mstar = 2 ×106 - 5 ×1010 Mo) using the Feedback in Realistic Environments (FIRE) simulations. While previous studies emphasized the impact of feedback on dark matter profiles, we investigate the impact on the stellar component: kinematics, radial migration, size evolution, and population gradients. Feedback-driven outflows/inflows drive significant radial stellar migration over both short and long timescales via two processes: (1) outflowing/infalling gas can remain star-forming, producing young stars that migrate ∼1 kpc within their first 100 Myr, and (2) gas outflows/inflows drive strong fluctuations in the global potential, transferring energy to all stars. These processes produce several dramatic effects. First, galaxies' effective radii can fluctuate by factors of >2 over ∼200 Myr, and these rapid size fluctuations can account for much of the observed scatter in the radius at fixed Mstar. Second, the cumulative effects of many outflow/infall episodes steadily heat stellar orbits, causing old stars to migrate outward most strongly. This age-dependent radial migration mixes - and even inverts - intrinsic age and metallicity gradients. Thus, the galactic-archaeology approach of calculating radial star formation histories from stellar populations at z = 0 can be severely biased. These effects are strongest at Mstar ≈ 107-9.6 Mo, the same regime where feedback most efficiently cores galaxies. Thus, detailed measurements of stellar kinematics in low-mass galaxies can strongly constrain feedback models and test baryonic solutions to small-scale problems in ΛCDM.

In preparation for the release of the astrometric orbits of Gaia , Shahaf et al. (2019) proposed a triage technique to identify astrometric binaries with compact companions based on their astrometric semimajor axis, parallax, and primary mass. The technique requires the knowledge of the appropriate mass–luminosity relation to rule out single or close-binary main-sequence companions. The recent publication of the Gaia DR3 astrometric orbits used a schematic version of this approach, identifying 735 astrometric binaries that might have compact companions. In this communication, we return to the triage of the DR3 astrometric binaries with more careful analysis, estimating the probability for its astrometric secondary to be a compact object or a main-sequence close binary. We compile a sample of 177 systems with highly probable non-luminous massive companions, which is smaller but cleaner than the sample reported in Gaia DR3. The new sample includes eight candidates to be black-hole systems with compact-object masses larger than 2.4 M⊙. The orbital–eccentricity–secondary–mass diagram of the other 169 systems suggests a tentative separation between the white-dwarf and the neutron-star binaries. Most white-dwarf binaries are characterized by small eccentricities of about 0.1 and masses of 0.6 M⊙, while the neutron star binaries display typical eccentricities of 0.4 and masses of 1.3 M⊙. [ABSTRACT FROM AUTHOR]