Breathing FIRE: How Stellar Feedback Drives Radial Migration, Rapid Size Fluctuations, and Population Gradients in Low-Mass Galaxies

We examine the effects of stellar feedback and bursty star formation on low-mass galaxies ($M_{\rm star}=2\times10^6-5\times10^{10}{\rm M_{\odot}}$) using the FIRE (Feedback in Realistic Environments) 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 $\sim1{\rm\,kpc}$ within their first $100 {\rm\,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 $\sim200 {\rm\,Myr}$, and these rapid size fluctuations can account for much of the observed scatter in radius at fixed $M_{\rm star}.$ 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 $M_{\rm star}\approx10^{7-9.6}{\rm M_{\odot}}$, 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 $\Lambda$CDM.Comment: Accepted to ApJ (820, 131) with minor revisions from v1. Figure 4 now includes dark matter. Main results in Figures 7 and 1

How To Model Supernovae in Simulations of Star and Galaxy Formation

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 >100 solar masses, 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. 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.Comment: 18 pages, 9 figures (+8 pages, 6 figures in appendices). MNRAS (updated to match published version

The origin of ultra diffuse galaxies: stellar feedback and quenching

We test if the cosmological zoom-in simulations of isolated galaxies from the FIRE project reproduce the properties of ultra diffuse galaxies. We show that stellar feedback-generated outflows that dynamically heat galactic stars, together with a passively aging stellar population after imposed quenching (from e.g. infall into a galaxy cluster), naturally reproduce the observed population of red UDGs, without the need for high spin halos or dynamical influence from their host cluster. We reproduce the range of surface brightness, radius and absolute magnitude of the observed z=0 red UDGs by quenching simulated galaxies at a range of different times. They represent a mostly uniform population of dark matter-dominated galaxies with M_star ~1e8 Msun, low metallicity and a broad range of ages. The most massive simulated UDGs require earliest quenching and are therefore the oldest. Our simulations provide a good match to the central enclosed masses and the velocity dispersions of the observed UDGs (20-50 km/s). The enclosed masses of the simulated UDGs remain largely fixed across a broad range of quenching times because the central regions of their dark matter halos complete their growth early. A typical UDG forms in a dwarf halo mass range of Mh~4e10-1e11 Msun. The most massive red UDG in our sample requires quenching at z~3 when its halo reached Mh ~ 1e11 Msun. If it, instead, continues growing in the field, by z=0 its halo mass reaches > 5e11 Msun, comparable to the halo of an L* galaxy. If our simulated dwarfs are not quenched, they evolve into bluer low-surface brightness galaxies with mass-to-light ratios 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 Ms~1e8 Msun, both in the field and in clusters.Comment: 20 pages, 13 figures; match the MNRAS accepted versio

A profile in FIRE: resolving the radial distributions of satellite galaxies in the Local Group with simulations

While many tensions between Local Group (LG) satellite galaxies and LCDM 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 Mstar > 10^5 Msun around 8 isolated Milky Way- (MW) mass host galaxies and 4 hosts in LG-like pairs. We demonstrate that these simulations resolve the survival and physical destruction of satellites with Mstar >~ 10^5 Msun. 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 disks of host galaxies. Furthermore, we quantify and provide fits to the tidal depletion of subhalos in baryonic relative to dark matter-only simulations as a function of distance. Our simulated profiles imply observational incompleteness in the LG even at Mstar >~ 10^5 Msun: 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 analogs in the SAGA survey, finding that our simulations are broadly consistent with most SAGA systems.Comment: 18 pages, 10 figures, plus appendices. Main results in figures 2, 3, and 4. Accepted versio

Towards harmonisation of chemical monitoring using avian apex predators: identification of key species for pan-European biomonitoring

Biomonitoring in raptors can be used to study long-term and large-scale changes in environmental pollution. In Europe, such monitoring is needed to assess environmental risks and outcomes of chemicals regulation, which is harmonised across the European Union. To be effective, the most appropriate sentinels need to be monitored. Our aim was to identify which European raptor species are the likely most appropriate biomonitors when pollutant quantification is based on analysing tissues. Our current study was restricted to terrestrial exposure pathways and considered four priority pollutant groups: toxic metals (lead and mercury), anticoagulant rodenticides, pesticides and medicinal products. We evaluated information on the distribution and key ecological traits (food web, foraging trait, diet, preferred habitat, and migratory behaviour) of European raptors to identify the most appropriate sentinel species. Common buzzard (Buteo buteo) and/or tawny owl (Strix aluco) proved the most suitable candidates for many of the pollutants considered. Moreover, they are abundant in Europe, enhancing the likelihood that samples can be collected. However, other species may be better sentinels for certain pollutants, such as the golden eagle (Aquila chrysaetos) for lead, the northern goshawk (Accipiter gentilis) for mercury across areas including Northern Europe, and vultures (where they occur in Europe) are likely best suited for monitoring non-steroidal anti-inflammatory drugs (NSAIDs). Overall, however, we argue the selection of candidate species for widescale monitoring of a range of pollutants can be reduced to very few raptor species. We recommend that the common buzzard and tawny owl should be the initial focus of any pan-European raptor monitoring. The lack of previous widespread monitoring using these species suggests that their utility as sentinels for environmnetal pollution has not been widely recognised. Finally, although the current study focussed on Europe, our trait-based approach for identifying raptor biomonitors can be applied to other continents and contaminants

Star formation histories of dwarf galaxies in the FIRE simulations: dependence on mass and Local Group environment

We study star formation histories (SFHs) of $\simeq500$ dwarf galaxies (stellar mass $M_\ast = 10^5 - 10^9\,M_\odot$) 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 vs. central" and "LG vs. individual MWvs. 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_\ast = 10^6 - 10^7\,M_\odot$ 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.Comment: Main text: 11 pages, 8 figures; appendices: 4 pages, 4 figures. Submitted to MNRAS; comments welcom

Probing the CGM of low-redshift dwarf galaxies using FIRE simulations

Observations of ultraviolet (UV) metal absorption lines have provided insight into the structure and composition of the circumgalactic medium (CGM) around galaxies. We compare these observations with the low-redshift (z ≤ 0.3) CGM around dwarf galaxies in high-resolution cosmological zoom-in runs in the FIRE-2 (Feedback In Realistic Environments) simulation suite. We select simulated galaxies that match the halo mass, stellar mass, and redshift of the observed samples. We produce absorption measurements using TRIDENT for UV transitions of C IV, O VI, Mg II, and Si III. The FIRE equivalent width (EW) distributions and covering fractions for the C IV ion are broadly consistent with observations inside 0.5R_(vir), but are underpredicted for O VI, Mg II, and Si III. The absorption strengths of the ions in the CGM are moderately correlated with the masses and star formation activity of the galaxies. The correlation strengths increase with the ionization potential of the ions. The structure and composition of the gas from the simulations exhibit three zones around dwarf galaxies characterized by distinct ion column densities: the discy interstellar medium, the inner CGM (the wind-dominated regime), and the outer CGM (the IGM accretion-dominated regime). We find that the outer CGM in the simulations is nearly but not quite supported by thermal pressure, so it is not in hydrostatic equilibrium, resulting in halo-scale bulk inflow and outflow motions. The net gas inflow rates are comparable to the star formation rate of the galaxy, but the bulk inflow and outflow rates are greater by an order of magnitude, with velocities comparable to the virial velocity of the halo. These roughly virial velocities (⁠∼100 km s⁻¹⁠) produce large EWs in the simulations. This supports a picture for dwarf galaxies in which the dynamics of the CGM at large scales are coupled to the small-scale star formation activity near the centre of their haloes