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

We present an empirical approach for interpreting gravitational wave signals of binary black hole mergers under the assumption that the underlying black hole population is sourced by remnants of stellar evolution. Using the observed relationship between galaxy mass and stellar metallicity, we predict the black hole count as a function of galaxy stellar mass. We show, for example, that a galaxy like the Milky Way should host millions of $\sim 30~M_\odot$ black holes and dwarf satellite galaxies like Draco should host $\sim 100$ such remnants, with weak dependence on the assumed IMF and stellar evolution model. Most low-mass black holes ($\sim10 M_\odot$) typically reside within massive galaxies ($M_\star \simeq 10^{11} M_\odot$) while massive black holes ($\sim 50~M_\odot$) typically reside within dwarf galaxies ($M_\odot \simeq 10^9 M_\odot$) today. If roughly $1\%$ of black holes are involved in a binary black hole merger, then the reported merger rate densities from Advanced LIGO can be accommodated for a range of merger timescales, and the detection of mergers with $> 50~M_\odot$ black holes should be expected within the next decade. Identifying the host galaxy population of the mergers provides a way to constrain both the binary neutron star or black hole formation efficiencies and the merger timescale distributions; these events would be primarily localized in dwarf galaxies if the merger timescale is short compared to the age of the universe and in massive galaxies otherwise. As more mergers are detected, the prospect of identifying the host galaxy population, either directly through the detection of electromagnetic counterparts of binary neutron star mergers or indirectly through the anisotropy of the events, will become a realistic possibility.

We present a suite of 15 cosmological zoom-in simulations of isolated dark matter halos, all with masses of $M_{\rm halo} \approx 10^{10}\,{\rm M}_\odot$ at $z=0$, in order to understand the relationship between 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. 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_\star/{\rm M}_\odot\approx 10^5-10^7$). This non-negligible spread in stellar mass at $z=0$ in halos within a narrow range of virial masses is strongly correlated with central halo density or maximum circular velocity $V_{\rm max}$. Much of this dependence of $M_\star$ on a second parameter (beyond $M_{\rm halo}$) is a direct consequence of the $M_{\rm halo}\sim10^{10}\,{\rm M}_\odot$ mass scale coinciding with the threshold for strong reionization suppression: the densest, earliest-forming halos 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. Halos that form galaxies with $M_\star\gtrsim2\times10^{6}\,{\rm M}_\odot$ 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 $r_{1/2}$. This apparent stellar mass threshold of $M_\star \approx 2\times 10^{6} \approx 2\times 10^{-4} \,M_{\rm halo}$ is broadly consistent with previous work and provides a testable prediction of FIRE feedback models in LCDM.

We compare a suite of four simulated dwarf galaxies formed in 10$^{10} M_{\odot}$ haloes of collisionless Cold Dark Matter (CDM) with galaxies simulated in the same haloes with an identical galaxy formation model but a non-zero cross-section for dark matter self-interactions. These cosmological zoom-in simulations are part of the Feedback In Realistic Environments (FIRE) project and utilize the FIRE-2 model for hydrodynamics and galaxy formation physics. We find the stellar masses of the galaxies formed in Self-Interacting Dark Matter (SIDM) with $\sigma/m= 1\, cm^2/g$ are very similar to those in CDM (spanning $M_{\star} \approx 10^{5.7 - 7.0} M_{\odot}$) and all runs lie on a similar stellar mass -- size relation. The logarithmic dark matter density slope ($\alpha=d\log \rho / d\log r$) in the central $250-500$ pc remains steeper than $\alpha= -0.8$ for the CDM-Hydro simulations with stellar mass $M_{\star} \sim 10^{6.6} M_{\odot}$ and core-like in the most massive galaxy. In contrast, every SIDM hydrodynamic simulation yields a flatter profile, with $\alpha >-0.4$. Moreover, the central density profiles predicted in SIDM runs without baryons are similar to the SIDM runs that include FIRE-2 baryonic physics. Thus, SIDM appears to be much more robust to the inclusion of (potentially uncertain) baryonic physics than CDM on this mass scale, suggesting SIDM will be easier to falsify than CDM using low-mass galaxies. Our FIRE simulations predict that galaxies less massive than $M_{\star} < 3 \times 10^6 M_{\odot}$ provide potentially ideal targets for discriminating models, with SIDM producing substantial cores in such tiny galaxies and CDM producing cusps.

We use a suite of high-resolution cosmological dwarf galaxy simulations to test the accuracy of commonly-used mass estimators from Walker et al.(2009) and Wolf et al.(2010), both of which depend on the observed line-of-sight velocity dispersion and the 2D half-light radius of the galaxy, $Re$. The simulations are part of the the Feedback in Realistic Environments (FIRE) project and include twelve systems with stellar masses spanning $10^{5} - 10^{7} M_{\odot}$ that have structural and kinematic properties similar to those of observed dispersion-supported dwarfs. Both estimators are found to be quite accurate: $M_{Wolf}/M_{true} = 0.98^{+0.19}_{-0.12}$ and $M_{Walker}/M_{true} =1.07^{+0.21}_{-0.15}$, with errors reflecting the 68% range over all simulations. The excellent performance of these estimators is remarkable given that they each assume spherical symmetry, a supposition that is broken in our simulated galaxies. Though our dwarfs have negligible rotation support, their 3D stellar distributions are flattened, with short-to-long axis ratios $ c/a \simeq 0.4-0.7$. The accuracy of the estimators shows no trend with asphericity. Our simulated galaxies have sphericalized stellar profiles in 3D that follow a nearly universal form, one that transitions from a core at small radius to a steep fall-off $\propto r^{-4.2}$ at large $r$, they are well fit by S\'ersic profiles in projection. We find that the most important empirical quantity affecting mass estimator accuracy is $Re$ . Determining $Re$ by an analytic fit to the surface density profile produces a better estimated mass than if the half-light radius is determined via direct summation.

We perform a systematic Bayesian analysis of rotation vs. dispersion support ($v_{\rm rot} / \sigma$) in $40$ dwarf galaxies throughout the Local Volume (LV) over a stellar mass range $10^{3.5} M_{\rm \odot} < M_{\star} < 10^8 M_{\rm \odot}$. We find that the stars in $\sim 80\%$ of the LV dwarf galaxies studied -- both satellites and isolated systems -- are dispersion-supported. In particular, we show that $6/10$ *isolated* dwarfs in our sample have $v_{\rm rot} / \sigma < 1.0$. All have $v_{\rm rot} / \sigma \lesssim 2.0$. These results challenge the traditional view that the stars in gas-rich dwarf irregulars (dIrrs) are distributed in cold, rotationally-supported stellar disks, while gas-poor dwarf spheroidals (dSphs) are kinematically distinct in having dispersion-supported stars. We see no clear trend between $v_{\rm rot} / \sigma$ and distance to the closest $\rm L_{\star}$ galaxy, nor between $v_{\rm rot} / \sigma$ and $M_{\star}$ within our mass range. We apply the same Bayesian analysis to four FIRE hydrodynamic zoom-in simulations of isolated dwarf galaxies ($10^9 M_{\odot} < M_{\rm vir} < 10^{10} M_{\rm \odot}$) and show that the simulated *isolated* dIrr galaxies have stellar ellipticities and stellar $v_{\rm rot} / \sigma$ ratios that are consistent with the observed population of dIrrs *and* dSphs without the need to subject these dwarfs to any external perturbations or tidal forces. We posit that most dwarf galaxies form as puffy, dispersion-dominated systems, rather than cold, angular momentum-supported disks. If this is the case, then transforming a dIrr into a dSph may require little more than removing its gas.

We investigate the response of self-interacting dark matter (SIDM) halos to the growth of galaxy potentials using idealized simulations, each run in tandem with standard collisionless Cold Dark Matter (CDM). We find a greater diversity in the SIDM halo profiles compared to the CDM halo profiles. If the stellar gravitational potential strongly dominates in the central parts of a galaxy, then SIDM halos can be as dense as CDM halos on observable scales. For extreme cases with highly compact disks core collapse can occur, leading to SIDM halos that are denser and cuspier than their CDM counterparts. If the stellar potential is not dominant, then SIDM halos retain constant density cores with densities far below CDM predictions. When a disk potential is present, the inner SIDM halo becomes \em{more flattened} in the disk plane than the CDM halo. These results are in excellent quantitative agreement with the predictions of Kaplinghat et al. (2014). We also simulated a galaxy cluster halo with a central stellar distribution similar to the brightest central galaxy of the cluster A2667. A SIDM halo simulated with cross section over mass $\sigma/m = 0.1\ \mathrm{cm^2 g^{-1}}$ provides a good match to the measured dark matter density profile of A2667, while an adiabatically-contracted CDM halo is denser and cuspier. The cored profile of the same halo simulated with $\sigma/m = 0.5\ \mathrm{cm^2 g^{-1}}$ is not dense enough to match A2667. Our findings are in agreement with previous results that $\sigma/m \gtrsim 0.1\ \mathrm{cm^2 g^{-1}}$ is disfavored for dark matter collision velocities in excess of about 1500 km/s. More generally, the predictive cross-talk between baryonic potentials and SIDM density distributions offers new directions for constraining SIDM cross sections in massive galaxies where baryons are dynamically important.

We present FIRE/Gizmo hydrodynamic zoom-in simulations of isolated dark matter halos, two each at the mass of classical dwarf galaxies ($M_{\rm vir} \simeq 10^{10} M_{\odot}$) and ultra-faint galaxies ($M_{\rm vir} \simeq 10^9 M_{\odot}$), and with two feedback implementations. The resultant central galaxies lie on an extrapolated abundance matching relation from $M_{\star} \simeq 10^6$ to $10^4 M_{\odot}$ without a break. Every host is filled with subhalos, many of which form stars. Our dwarfs with $M_{\star} \simeq 10^6 M_{\odot}$ each have 1-2 well-resolved satellites with $M_{\star} = 3-200 \times 10^3 M_{\odot}$. Even our isolated ultra-faint galaxies have star-forming subhalos. If this is representative, dwarf galaxies throughout the universe should commonly host tiny satellite galaxies of their own. We combine our results with the ELVIS simulations to show that targeting $\sim 50~ \rm kpc$ regions around nearby isolated dwarfs could increase the chances of discovering ultra-faint galaxies by $\sim 35\%$ compared to random halo pointings, and specifically identify the region around the Phoenix dwarf galaxy as a good potential target. The well-resolved ultra-faint galaxies in our simulations ($M_{\star} \simeq 3 - 30 \times 10^3 M_{\odot}$) form within $M_{\rm peak} \simeq 0.5 - 3 \times 10^9 M_{\odot}$ halos. Each has a uniformly ancient stellar population ($ > 10~ \rm Gyr$) owing to reionization-related quenching. More massive systems, in contrast, all have late-time star formation. Our results suggest that $M_{\rm halo} \simeq 5 \times 10^9 M_{\odot}$ is a probable dividing line between halos hosting reionization "fossils" and those hosting dwarfs that can continue to form stars in isolation after reionization.

We investigate the effect of self-interacting dark matter (SIDM) on the density profiles of $V_{\rm max} \simeq 40~km~s^{-1}$ isolated dwarf dark matter halos -- the scale of relevance for the too big to fail problem (TBTF) -- using very high-resolution cosmological zoom simulations. Each halo has millions of particles within its virial radius. We find that SIDM models with cross sections per unit mass spanning the range \sigma/m = $0.5 - 50$ $cm^2~g^{-1}$ alleviate TBTF and produce constant density cores of size 300-1000 pc, comparable to the half-light radii of $M_\star$ ~ $10^{5-7}$ $M_\odot$ dwarfs. The largest, lowest density cores develop for cross sections in the middle of this range, \sigma/m ~ $5-10~cm^2~g^{-1}$. Our largest SIDM cross section run (\sigma/m = $50~cm^2~g^{-1}$) develops a slightly denser core owing to mild core-collapse behavior, but it remains less dense than the CDM case and retains a constant density core profile. Our work suggests that SIDM cross sections as large or larger than $50~cm^2~g^{-1}$ remain viable on velocity scales of dwarf galaxies ($v_{\rm rms}$ ~ $40~km~s^{-1}$). The range of SIDM cross sections that alleviate TBTF and the cusp/core problem spans at least two orders of magnitude and therefore need not be particularly fine-tuned.