We analyze the cold dark matter density profiles of 54 galaxy halos simulated with FIRE-2 galaxy formation physics, each resolved within $0.5\%$ of the halo virial radius. These halos contain galaxies with masses that range from ultra-faint dwarfs ($M_\star \simeq 10^{4.5} M_{\odot}$) to the largest spirals ($M_\star \simeq 10^{11} M_{\odot}$) and have density profiles that are both cored and cuspy. We characterize our results using a new analytic density profile that extends the standard Einasto form to allow for a pronounced constant-density core in the resolved innermost radius. With one additional core-radius parameter, $r_{c}$, this "core-Einasto" profile is able to characterize the shape and normalization of our feedback-impacted dark matter halos. In order to enable comparisons with observations, we provide fitting functions for $r_{c}$ and other profile parameters as a function of both $M_\star$ and $M_{\star}/M_{\rm halo}$. In agreement with similar studies done in the literature, we find that dark matter core formation is most efficient at the characteristic stellar-mass to halo-mass ratio $M_\star/M_{\rm halo} \simeq 5 \times 10^{-3}$, or $M_{\star} \sim 10^9 \, M_{\odot}$, with cores that are roughly the size of the galaxy half-light radius, $r_{c} \simeq 1-5$ kpc. Furthermore, we find no evidence for core formation at radii $\gtrsim 100\ \rm pc$ in galaxies with $M_{\star}/M_{\rm halo} < 5\times 10^{-4}$ or $M_\star \lesssim 10^6 \, M_{\odot}$. For Milky Way-size galaxies, baryonic contraction often makes halos significantly more concentrated and dense at the stellar half-light radius than dark matter only runs. However, even at the Milky Way scale, FIRE-2 galaxy formation still produces small dark matter cores of $\simeq 0.5-2$ kpc in size. Recent evidence for a ${\sim} 2$ kpc core in the Milky Way's dark matter halo is consistent with this expectation.

We examine the thermodynamic state and cooling of the low-$z$ Circum-Galactic Medium (CGM) in five FIRE-2 galaxy formation simulations of Milky Way-mass galaxies. We find that the CGM in these simulations is generally multiphase and dynamic, with a wide spectrum of largely nonlinear density perturbations sourced by the accretion of gas from the IGM and winds from both the central and satellite galaxies. We investigate the origin of the multiphase structure of the CGM with a particle tracking analysis and find that most of the low entropy gas has cooled from the hot halo as a result of thermal instability triggered by these perturbations. The ratio of cooling to free-fall timescales $t_{\rm cool}/t_{\rm ff}$ in the hot component of the CGM spans a wide range $\sim 1-100$ at a given radius, but exhibits approximately constant median values $\sim 5-20$ at all radii $0.1 R_{\rm vir} < r < R_{\rm vir}$. These are similar to the $\approx 10-20$ value typically adopted as the thermal instability threshold in "precipitation" models of the ICM. Consequently, a one-dimensional model based on the assumption of a constant $t_{\rm cool}/t_{\rm ff}$ and hydrostatic equilibrium approximately reproduces the simulation number density and entropy profiles, but only if it assumes the metallicity profiles taken directly from the simulations. We explicitly show that the $t_{\rm cool}/t_{\rm ff}$ value of a gas parcel in the hot component of the CGM does not predict its probability of cooling and subsequently accreting onto the central galaxy. This suggests that the value of $t_{\rm cool}/t_{\rm ff}$ is a poor predictor of thermal stability in gaseous halos in which large-amplitude density perturbations are prevalent.

We use the FIRE-2 cosmological simulations to study the formation of a virial temperature, quasi-static 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 discs. We demonstrate that when the halo mass crosses ~10^12 M_sun, the cooling time of shocked gas in the inner CGM (~0.1 R_vir, where R_vir is the virial radius) exceeds the local free-fall time. The inner CGM then experiences a transition from on average sub-virial temperatures (T<<T_vir), large pressure fluctuations and supersonic inflow/outflow velocities, to virial temperatures (T~T_vir), uniform pressures and subsonic velocities. This transition occurs when the outer CGM (~0.5 R_vir) is already subsonic and has a temperature ~T_vir, 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 virialization of the inner CGM coincides with abrupt changes in the properties of the central galaxy and its stellar feedback: the galaxy settles into a stable rotating disc, star formation transitions from `bursty' to `steady,' and stellar-driven galaxy-scale outflows are suppressed. Our results hence suggest that CGM virialization is initially associated with the formation of rotation-dominated thin galactic discs, rather than with the quenching of star formation as often assumed.

We study the escape fraction of ionizing photons (f_esc) in two cosmological zoom-in simulations of galaxies in the reionization era with halo mass M_halo~10^10 and 10^11 M_sun (stellar mass M*~10^7 and 10^9 M_sun) at z=5 from the Feedback in Realistic Environments project. These simulations explicitly resolve the formation of proto-globular clusters (GCs) self-consistently, where 17-39% of stars form in bound clusters during starbursts. Using post-processing Monte Carlo radiative transfer calculations of ionizing radiation, we compute f_esc from cluster stars and non-cluster stars formed during a starburst over ~100 Myr in each galaxy. We find that the averaged f_esc over the lifetime of a star particle follows a similar distribution for cluster stars and non-cluster stars. Clusters tend to have low f_esc in the first few Myrs, presumably because they form preferentially in more extreme environments with high optical depths; the f_esc increases later as feedback starts to disrupt the natal cloud. On the other hand, non-cluster stars formed between cluster complexes or in the compressed shell at the front of a superbubble can also have high f_esc. We find that cluster stars on average have comparable f_esc to non-cluster stars. This result is robust across several star formation models in our simulations. Our results suggest that the fraction of ionizing photons from proto-GCs to cosmic reionization is comparable to the cluster formation efficiency in high-redshift galaxies and hence proto-GCs likely contribute an appreciable fraction of photons but are not the dominant sources for reionization.

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 ($\sigma_{star}$, $\sigma_{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 $\sigma_{star}$ (and $\sigma_{gas}$) and sSFR indeed exists for galaxies in the lowest mass bin (M$_\ast\sim10^{10}\,$M$_\odot$). This correlation, albeit with a $\sim$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$_\ast\lesssim10^{10}\,$M$_\odot$ 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 $\sigma_{star}$ (and $\sigma_{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.

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 using a set of FIRE-2 cosmological zoom-in simulations of Milky Way-mass disk galaxies, in which a multiphase ISM is self-consistently shaped by gravity, cooling, and stellar feedback. We analyze 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 midplane. Bulk flows (e.g., inflows and fountains) are important at a few disk scale heights, while thermal pressure from hot gas dominates at larger heights. Overall, the total midplane pressure is well-predicted by the weight of the disk gas, and we show that it also scales linearly with the star formation rate surface density (Sigma_SFR). These results support the notion that the Kennicutt-Schmidt relation arises because Sigma_SFR and the gas surface density (Sigma_g) are connected via the ISM midplane pressure.

The diffuse plasma inside clusters of galaxies has X-ray emitting temperatures of a few keV. The physical mechanisms that heat this intracluster medium (ICM) to such temperatures include the accretion shock at the periphery of a galaxy cluster, the shocks driven by merger events, as well as a somewhat overlooked mechanism -- the dissipation of intracluster turbulent motions. We study the relative role of these heating mechanisms using galaxy clusters in Lagrangian tracer particle re-simulations of the Omega500 cosmological simulation. We adopt a novel analysis method of decomposing the temperature increase at each time step into the contribution from dissipative heating and that from adiabatic heating. In the high-resolution spatial-temporal map of these heating rates, merger tracks are clearly visible, demonstrating the dominant role of merger events in heating the ICM. The dissipative heating contributed by each merger event is extended in time and also occurs in the rarefaction regions, suggesting the importance of heating by the dissipation of merger-induced turbulence. Quantitative analysis shows that turbulence heating, rather than direct heating at merger shocks, dominates the temperature increase of the ICM, especially at inner radii $r < r_{\rm 500c}$. In addition, we find that many merger shocks can propagate with almost constant velocity to very large radii $r \gg r_{\rm 500c}$, some even reach and join with the accretion shock and becoming the outer boundary of the ICM. Altogether, these results suggest that the ICM is heated more in an `inside-out' fashion rather than `outside-in' as depicted in the classical smooth accretion picture.

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. 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. The result is a catalog that identifies around 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.

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 (FIRE) 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_{\odot}$) 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.

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.

We present the escape fraction of hydrogen ionizing photons (f_esc) 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 halos in M_vir~10^8--10^{12} M_sun (M_star~10^4--10^{10} M_sun) at z=5--12. We find the sample average <f_esc> increases with halo mass for M_vir~10^8--10^{9.5} M_sun, becomes nearly constant for M_vir~10^{9.5}--10^{11} M_sun, and decreases at M_vir>10^{11} M_sun. Equivalently, <f_esc> increases with stellar mass up to M_star~10^8 M_sun and decreases at higher masses. Even applying single-star stellar population synthesis models, we find a moderate <f_esc>~0.2 for galaxies at M_star~10^8 M_sun. 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 <f_esc> by ~25--35% and the number of escaped photons by 60--80%. 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 <f_esc> at the high-mass end is due to dust attenuation, while at the low-mass end, <f_esc> decreases owing to inefficient star formation (and hence feedback). At fixed mass, <f_esc> tends to increase with redshift. Our simulations produce sufficient ionizing photons for cosmic reionization.

Giant molecular clouds (GMCs) are well-studied in the local Universe, however, exactly how their properties vary during galaxy evolution is poorly understood due to challenging resolution requirements, both observational and computational. We present the first time-dependent analysis of giant molecular clouds in a Milky Way-like galaxy and an LMC-like dwarf galaxy of the FIRE-2 (Feedback In Realistic Environments) simulation suite, which have sufficient resolution to predict the bulk properties of GMCs in cosmological galaxy formation self-consistently. We show explicitly that the majority of star formation outside the galactic center occurs within self-gravitating gas structures that have properties consistent with observed bound GMCs. We find that the typical cloud bulk properties such as mass and surface density do not vary more than a factor of 2 in any systematic way after the first Gyr of cosmic evolution within a given galaxy from its progenitor. While the median properties are constant, the tails of the distributions can briefly undergo drastic changes, which can produce very massive and dense self-gravitating gas clouds. Once the galaxy forms, we identify only two systematic trends in bulk properties over cosmic time: a steady increase in metallicity produced by previous stellar populations and a weak decrease in bulk cloud temperatures. With the exception of metallicity we find no significant differences in cloud properties between the Milky Way-like and dwarf galaxies. These results have important implications for cosmological star and star cluster formation and put especially strong constraints on theories relating the stellar initial mass function to cloud properties.

We report the formation of bound star clusters in a sample of high-resolution cosmological zoom-in simulations of z>5 galaxies from the FIRE project. We find that bound clusters preferentially form in high-pressure clouds with gas surface densities over 10^4 Msun 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 stellar masses of ~10^7-10^10 Msun at z>5. The age spread of cluster stars is typically a few Myrs and increases with cluster mass. The metallicity dispersion of cluster members is ~0.08 dex in [Z/H] and does not depend on cluster mass 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.

Radiative feedback (RFB) from stars plays a key role in galaxies, but remains poorly-understood. We explore this using high-resolution, multi-frequency radiation-hydrodynamics (RHD) simulations from the Feedback In Realistic Environments (FIRE) project. We study ultra-faint dwarf through Milky Way mass scales, including H+He photo-ionization; photo-electric, 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 photo-ionization 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 photo-ionization heating from the UV background suppressing accretion. In MW-mass galaxies, meta-galactic backgrounds have negligible effects; but local photo-ionization 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 hyper-dense, bound star clusters. This makes star formation more violent and 'bursty' when SNe explode in these hyper-clustered 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% of RP in massive galaxies): absorption occurs primarily in 'normal' GMCs with A_v~1.

The correlation between galaxies' integrated stellar masses and star formation rates (the `star formation main sequence'; SFMS) is a well-established scaling relation. Recently, surveys have found a relationship between the star formation rate and stellar mass surface densities on kpc and sub-kpc scales (the `resolved SFMS'; rSFMS). In this work, we demonstrate that the rSFMS emerges naturally in FIRE-2 zoom-in simulations of Milky Way-mass galaxies. We make SFR and stellar mass maps of the simulated galaxies at a variety of spatial resolutions and star formation averaging time-scales and fit the rSFMS using multiple methods from the literature. While the absolute value of the SFMS slope depends on the fitting method, the slope is steeper for longer star formation time-scales and lower spatial resolutions regardless of the fitting method employed. We present a toy model that quantitatively captures the dependence of the simulated galaxies' rSFMS slope on spatial resolution and use it to illustrate how this dependence can be used to constrain the characteristic mass of star-forming clumps.

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. Commonly used choices for the shielding length, such as the Jeans length, result in CO abundances that are too high at a given H$_2$ abundance. We find that a model with a distribution of shielding lengths, which has a median shielding length of $\sim 3$ pc in cold gas ($T < 300$ K) for both CO and H$_{2}$, is able to reproduce both the observed CO(1-0) luminosity and inferred CO-to-H$_{2}$ 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.

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 a 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 FIRE-2 cosmological simulations, including 6 isolated MW/M31-mass galaxies and 6 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 a MW/M31-mass galaxy emerged typically at z ~ 3-4 (11.6-12.2 Gyr ago), while stars in the bulge region (inner 2 kpc) at z = 0 formed primarily in a single main progenitor at z < 5 (< 12.6 Gyr ago). Compared with isolated hosts, the main progenitors of LG-like paired hosts emerged significantly earlier (\Delta z ~ 2, \Delta t ~ 1.6 Gyr), with ~ 4x higher stellar mass at all z > 4 (> 12.2 Gyr ago). This highlights the importance of environment in MW/M31-mass galaxy formation, especially at early times. Overall, about 100 galaxies with M_star > 10^5 M_sun formed a typical MW/M31-mass system. Thus, surviving satellites represent a highly incomplete census (by ~ 5x) of the progenitor population.

The past two decades have seen a tremendous investment in observational facilities that promise to reveal new and unprecedented discoveries about the universe. In comparison, the investment in theoretical work is completely dwarfed, even though theory plays a crucial role in the interpretation of these observations, predicting new types of phenomena, and informing observing strategies. In this white paper, we argue that in order to reach the promised critical breakthroughs in astrophysics over the next decade and well beyond, the national agencies must take a serious approach to investment in theoretical astrophysics research. We discuss the role of theory in shaping our understanding of the universe, and then we provide a multi-level strategy, from the grassroots to the national, to address the current underinvestment in theory relative to observational work.

We study stellar-halo formation using six Milky Way-mass galaxies in FIRE-2 cosmological zoom simulations. We find that $5-40\%$ of the outer ($50-300$ kpc) stellar halo in each system consists of $\textit{in-situ}$ stars that were born in outflows from the main galaxy. Outflow stars originate from gas accelerated by super-bubble 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 $\sim 10 \%$ of our kinematically-identified halo stars were born in outflows; the fraction rises to as high as $\sim 40\%$ for the most metal-rich local halo stars ([Fe/H] $> -0.5$). We conclude that the Milky Way 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 $\textit{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.

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-mass galaxies. We track GMCs with total gas mass $\gtrsim 10^5$ M$_\odot$ at high spatial ($\sim1$ pc), mass ($7100$ M$_{\odot}$), and temporal (1 Myr) resolution. Our simulated GMCs are consistent with the distribution of masses for massive GMCs in the Milky Way and nearby galaxies. We find GMC lifetimes of $5-7$ Myr, or 1-2 freefall times, on average, with less than 1\% of clouds living longer than 20 Myr. We find 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 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 ten, 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 is 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.

We present constraints on the timescale of star formation variability and the correlation between star formation and host halo accretion histories of star-forming (SF) central galaxies from the scatter of the stellar-to-halo mass relation (SHMR). SF galaxies are found to have a tight relationship between their star formation rates and stellar masses on the so-called "star-forming sequence" (SFS), which characterizes their SF histories and $M_*$ growths. Meanwhile, observed SHMR constraints connect $M_*$ growth to halo accretion history. Combining these observed trends with a cosmological $N$-body simulation, we present flexible models that track the SF, $M_*$, and halo accretion histories of SF central galaxies at $z<1$ while reproducing the observed stellar mass function and SFS in SDSS DR7. We find that the scatter in SHMR at $M_h=10^{12}M_\odot$, $\sigma_{M_*|M_h=10^{12}M_\odot}$, is sensitive to the timescale of star formation variability, $t_{\rm duty}$, and the correlation coefficient, $r$, between SF and halo accretion histories: shorter $t_{\rm duty}$ and higher $r$ both result in tighter $\sigma_{M_*|M_h=10^{12}M_\odot}$. To reproduce a constant 0.2 dex scatter over $z=1-0$, our models require $t_{\rm duty}\leq1.5$ Gyr for $r=0.99$ or $r>0.6$ for $t_{\rm duty}=0.1$ Gyr. For $r\sim0.6$, as found in the literature, $t_{\rm duty}<0.2$ Gyr is necessary. Meanwhile, to reproduce the tightening of $\sigma_{M_*|M_h=10^{12}M_\odot}=0.35$ to 0.2 dex from $z=1-0$ in hydrodynamical simulations, our models require $t_{\rm duty}=0.1$ Gyr for $r>0.5$. Although, the lack of consensus on $\sigma_{M_*|M_h=10^{12}M_\odot}$ at $M_h=10^{12}M_\odot$ and at $z=1$ from observations and galaxy formation models remains the main bottleneck in precisely constraining $r$ and $t_{\rm duty}$, we demonstrate that SHMR can be used to constrain the SF and host halo accretion histories of SF central galaxies.

We analyze the different fates of the circumgalactic medium (CGM) in FIRE-2 cosmological simulations, focusing on the redshifts z=0.25 and z=2 representative of recent surveys. Our analysis includes 21 zoom-in simulations covering the halo mass range Mh(z=0) ~ 10^10 - 10^12 Msun. We analyze both where the gas ends up after first leaving the CGM (its "proximate" fate), as well as 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) ~ 5e11 Msun halos, but most of the CGM in lower-mass halos ends up in the IGM. This is so even though most of the CGM in M_h(z=2) ~ 5e10 Msun halos first accretes onto 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 analyzed. Of the CGM gas that subsequently accretes onto the central galaxy in the progenitors of Mh(z=0) ~10^12 Msun halos, most of it is cool (T~10^4 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 z=2 >~80% of the cool gas in Mh >~ 10^11 Msun halos will accrete onto 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 $M_\mathrm{halo} \approx 10^{10}\,$M$_\odot$ 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 (r$_{1/2}<500$ pc) have lower $\sigma_\star/V_\mathrm{max}$ ratios, reinforcing the idea that smaller dwarfs may reside in halos 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 $\Lambda$CDM. The V$_{1/2}-$r$_{1/2}$ relation across all simulations is generally consistent with observations of Local Field dwarfs, though compact objects such as Tucana provide a unique challenge. Spatially-resolved rotation curves in the central regions ($<400$ pc) of small dwarfs could provide a way to distinguish between CDM, WDM, and SIDM, however: at the masses probed in this simulation suite, cored density profiles in dwarfs with small r$_{1/2}$ values can only originate from dark matter self-interactions.

We present the star formation history (SFH) of the isolated (D~970 kpc) Local Group dwarf galaxy WLM measured from color-magnitude diagrams constructed from deep Hubble Space Telescope imaging. Our observations include a central (0.5 $r_h$) and outer field (0.7 $ r_h$) that reach below the oldest main sequence turnoff. WLM has no early dominant episode of star formation: 20% of its stellar mass formed by ~12.5 Gyr ago (z~5). It also has an SFR that rises to the present with 50% 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 the Feedback In Realistic Environment suite are in good agreement for $M_{\star}(z=0) \sim 10^7-10^9 M_{\odot}$, but in worse agreement at lower masses ($M_{\star}(z=0) \sim 10^5-10^7 M_{\odot}$). 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.

Within $\Lambda$CDM, dwarf galaxies like the Large Magellanic Cloud (LMC) are expected to host numerous dark matter subhalos, 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_\ast$ > $10^5$ M$_{\odot}$; Carina and Fornax) as well as several ultra-faint 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_\ast$ ~ $6\times10^3$ M$_{\odot}$) substructure population of isolated LMC-mass hosts ($M_\text{200}$ = $1-3\times10^{11}$ M$_{\odot}$) and place the Gaia + DES results in a cosmological context. By comparing number counts of subhalos in simulations with and without baryons, we find that, within 0.2 $r_\text{200}$, LMC-mass hosts deplete ~30% of their substructure, significantly lower than the ~70% of substructure depleted by Milky Way (MW) mass hosts. For our highest resolution runs ($m_\text{bary}$ = 880 M$_{\odot}$), ~5-10 subhalos form galaxies with $M_\ast$ > $10^4$ M$_{\odot}$, in agreement with the 7 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 ultra faint 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.

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. We present evidence for Nyx, a vast new stellar stream in the vicinity of the Sun, that may provide the first indication that such an event occurred in the Milky Way. We identify about 500 stars that have coherent radial and prograde motion in this stream using a catalog of accreted stars built by applying deep learning methods to the second Gaia data release. Nyx is concentrated within $\pm 2$ kpc of the Galactic midplane and spans the full radial range studied (6.5-9.5 kpc). The kinematics of Nyx stars are distinct from those of both the thin and thick disk. In particular, its rotational speed lags the disk by $\sim 80$ km/s and its stars follow more eccentric orbits. A small number of Nyx stars have chemical abundances or inferred ages; from these, we deduce that Nyx stars have a peak metallicity of [Fe/H] $\sim -0.5$ and ages $\sim $10-13 Gyr. Taken together with the kinematic observations, these results strongly favor the interpretation that Nyx is the remnant of a disrupted dwarf galaxy. To further justify this interpretation, we explicitly demonstrate that metal-rich, prograde streams like Nyx can be found in the disk plane of Milky Way-like galaxies using the FIRE hydrodynamic simulations. Future spectroscopic studies will be able to validate whether Nyx stars originate from a single progenitor.

In response to the need for the Astro2020 Decadal Survey to explicitly engage early career astronomers, the National Academies of Sciences, Engineering, and Medicine hosted the Early Career Astronomer and Astrophysicist Focus Session (ECFS) on October 8-9, 2018 under the auspices of Committee of Astronomy and Astrophysics. The meeting was attended by fifty six pre-tenure faculty, research scientists, postdoctoral scholars, and senior graduate students, as well as eight former decadal survey committee members, who acted as facilitators. The event was designed to educate early career astronomers about the decadal survey process, to solicit their feedback on the role that early career astronomers should play in Astro2020, and to provide a forum for the discussion of a wide range of topics regarding the astrophysics career path. This white paper presents highlights and themes that emerged during two days of discussion. In Section 1, we discuss concerns that emerged regarding the coming decade and the astrophysics career path, as well as specific recommendations from participants regarding how to address them. We have organized these concerns and suggestions into five broad themes. These include (sequentially): (1) adequately training astronomers in the statistical and computational techniques necessary in an era of "big data", (2) responses to the growth of collaborations and telescopes, (3) concerns about the adequacy of graduate and postdoctoral training, (4) the need for improvements in equity and inclusion in astronomy, and (5) smoothing and facilitating transitions between early career stages. Section 2 is focused on ideas regarding the decadal survey itself, including: incorporating early career voices, ensuring diverse input from a variety of stakeholders, and successfully and broadly disseminating the results of the survey.

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 log(Mh/Msun) ~ 10-12 , and we focus on z=0.25 and z=2. Owing to strong stellar feedback, only ~L* halos retain a baryon mass >~50% of their cosmic budget. Metals are more efficiently retained by halos, with a retention fraction >~50%. Across all masses and redshifts analyzed >~60% of the CGM mass originates as IGM accretion (some of which is associated with infalling halos). 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* halos. Gas can persist in the CGM for billions of years, resulting in well-mixed halo gas. Sight lines through the CGM are therefore likely to intersect gas of multiple origins. For low-redshift ~L* halos, cool gas (T<10^4.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_{\odot}$) < 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 between sSFR, $M_*$ and H$\alpha$ velocity dispersion, $\sigma$, from the Feedback In Realistic Environments (FIRE) simulations. Our observations reveal a correlation between sSFR and sigma 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 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 degrees), 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.

The Great Andromeda Galaxy (M31) is the nexus of the near-far galaxy evolution connection and a principal data point for near-field cosmology. Due to its proximity (780 kpc), M31 can be resolved into individual stars like the Milky Way (MW). Unlike the MW, we have the advantage of a global view of M31, enabling M31 to be observed with techniques that also apply to more distant galaxies. Moreover, recent evidence suggests that M31 may have survived a major merger within the last several Gyr, shaping the morphology of its stellar halo and triggering a starburst, while leaving the stellar disk largely intact. The MW and M31 thus provide complementary opportunities for in-depth studies of the disks, halos, and satellites of L* galaxies. Our understanding of the M31 system will be transformed in the 2020s if they include wide field facilities for both photometry (HST-like sensitivity and resolution) and spectroscopy (10-m class telescope, >1 sq. deg. field, highly multiplexed, R~ 3000 to 6000). We focus here on the power of these facilities to constrain the past, present, and future merger history of M31, via chemo-dynamical analyses and star formation histories of phase-mixed stars accreted at early times, as well as stars in surviving tidal debris features, M31's extended disk, and intact satellite galaxies that will eventually be tidally incorporated into the halo. This will yield an unprecedented view of the hierarchical formation of the M31 system and the subhalos that built it into the L* galaxy we observe today.

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.

We investigate the lowest-mass quiescent galaxies known to exist in isolated environments ($\mathrm{M^* = 10^{9.0-9.5} \ M_\odot}$; 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 Keck/ESI long-slit 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 AGN-like line ratios. Only 5 of these quiescent galaxies were identified as AGN-like in SDSS due to lower spatial resolution and signal-to-noise. We find that many of the quiescent galaxies in our sample have spatially-extended emission across the non-SF regions of BPT-space. When considering only the central 1$^{\prime\prime}$, we identify a tight relationship between distance from the BPT star-forming sequence and host galaxy stellar age as traced by $\mathrm{D_n4000}$, such that older stellar ages are associated with larger distances from the star-forming locus. Our results suggest that the presence of hard ionizing radiation (AGN-like line ratios) is intrinsically tied to the quenching of what may be the lowest-mass self-quenched galaxies.

We explore the radial variation of star formation histories in dwarf galaxies simulated with Feedback In Realistic Environments (FIRE) physics. The sample contains 9 low-mass field dwarfs with M_ star = 10^5 - 10^7 M_sun from previous FIRE results, and a new suite of 17 higher mass field dwarfs with M_star = 10^7 - 10^9 M_sun introduced here. We find that 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 star formation history of the merging galaxy. In galaxies without significant mergers, early feedback pushes stars to the outskirts at early times. Interestingly, among galaxies without mergers, those with large dark matter cores have flatter age gradients because these galaxies have more late-time feedback. 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 star formation history. We show that central fields can be biased young by a few 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 star formation history.

We study a suite of extremely high-resolution cosmological FIRE simulations of dwarf galaxies ($M_{\rm halo} \lesssim 10^{10}$$M_{\odot}$), run to $z=0$ with $30 M_{\odot}$ resolution, sufficient (for the first time) to resolve the internal structure of individual supernovae remnants within the cooling radius. Every halo with $M_{\rm halo} \gtrsim 10^{8.6} M_{\odot}$ is populated by a resolved {\em stellar} galaxy, suggesting very low-mass dwarfs may be ubiquitous in the field. Our ultra-faint dwarfs (UFDs; $M_{\ast}<10^{5}\,M_{\odot}$) have their star formation truncated early ($z\gtrsim2$), likely by reionization, while classical dwarfs ($M_{\ast}>10^{5} M_{\odot}$) continue forming stars to $z<0.5$. The systems have bursty star formation (SF) histories, forming most of their stars in periods of elevated SF strongly clustered in both space and time. This allows our dwarf with $M_{\ast}/M_{\rm halo} > 10^{-4}$ to form a dark matter core $>200$pc, while lower-mass UFDs exhibit cusps down to $\lesssim100$pc, as expected from energetic arguments. Our dwarfs with $M_{\ast}>10^{4}\,M_{\odot}$ have half-mass radii ($R_{\rm 1/2}$) in agreement with Local Group (LG) dwarfs; dynamical mass vs. $R_{1/2}$ and the degree of rotational support 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).

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

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 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 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. After demonstrating 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 $42 ^{+26}_{-22}\%$ of the local dark matter that is accreted from luminous mergers is in debris flow.

The Halo Assembly in Lambda-CDM: Observations in 7 Dimensions (HALO7D) dataset consists of Keck II/DEIMOS spectroscopy and Hubble Space Telescope-measured proper motions of Milky Way (MW) halo main sequence turnoff stars in the CANDELS fields. In this paper, the second in the HALO7D series, we present the proper motions for the HALO7D sample. We discuss our measurement methodology, which makes use of a Bayesian mixture modeling approach for creating the stationary reference frame of distant galaxies. Using the 3D kinematic HALO7D sample, we estimate the parameters of the halo velocity ellipsoid, $\langle v_{\phi} \rangle, \sigma_r, \sigma_{\phi}, \sigma_{\theta}$, and the velocity anisotropy $\beta$. Using the full HALO7D sample, we find $\beta=0.63 \pm 0.05$ at $\langle r \rangle =24$ kpc. We also estimate the ellipsoid parameters for our sample split into three apparent magnitude bins; the posterior medians for these estimates of $\beta$, while consistent with one another, increase as a function of mean sample distance. Finally, we estimate $\beta$ in each of the individual HALO7D fields. We find that the velocity anisotropy $\beta$ can vary from field to field, which suggests that the halo is not phase mixed at $\langle r \rangle =24$ kpc. We explore the $\beta$ variation across the skies of two stellar halos from the \textit{Latte} suite of FIRE-2 simulations, finding that both simulated galaxies show $\beta$ variation over a similar range to the variation observed across the four HALO7D fields. The accretion histories of the two simulated galaxies result in different $\beta$ variation patterns; spatially mapping $\beta$ is thus a way forward in characterizing the accretion history of the Galaxy.

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.

Binary black holes are the primary endpoint of massive stellar evolution. Their properties provide a unique opportunity to constrain binary evolution, which is still poorly understood. In this paper, we predict the inventory of binary black holes and their merger products in/around the Milky Way, and detail their main properties. We present the first combination of a high-resolution cosmological simulation of a Milky Way-mass galaxy with a binary population synthesis model. The hydrodynamic simulation, taken from the FIRE project, provides a cosmologically realistic star formation history for the galaxy and its stellar halo and satellites. We apply a metallicity-dependent evolutionary model to the star particles to produce individual binary black holes. We find that a million binary black holes have merged in the model Milky Way, and 3 million binaries are still present, with an average mass of 28 Msun per binary. Because the black hole progenitors are biased towards low metallicity stars, half reside in the stellar halo and satellites and 40 per cent of the binaries were formed outside the main galaxy. This trend increases with the masses of the black holes. The numbers and mass distribution of the merged systems is compatible with the LIGO/Virgo detections. Observations of these black holes will be challenging, both with electromagnetic methods and LISA. We find that a cosmologically realistic star formation history, with self-consistent metal enrichment and Galactic accretion history, are key ingredients for determining binary black hole rates that can be compared with observations to constrain massive binary evolution.

The oldest stars in the Milky Way (MW) bear imprints of the Galaxy's early assembly history. We use FIRE cosmological zoom-in simulations of three MW-mass disk galaxies to study the spatial distribution, chemistry, and kinematics of the oldest surviving stars ($z_{\rm form} \gtrsim 5$) in MW-like galaxies. We predict the oldest stars to be less centrally concentrated at $z=0$ than stars formed at later times as a result of two processes. First, the majority of the oldest stars are not formed $\textit{in situ}$ but are accreted during hierarchical assembly. These $\textit{ex situ}$ stars are deposited on dispersion-supported, halo-like orbits but dominate over old stars formed $\textit{in situ}$ in the solar neighborhood, and in some simulations, even in the galactic center. Secondly, old stars formed $\textit{in situ}$ are driven outwards by bursty star formation and energetic feedback processes that create a time-varying gravitational potential at $z\gtrsim 2$, similar to the process that creates dark matter cores and expands stellar orbits in bursty dwarf galaxies. The total fraction of stars that are ancient is more than an order of magnitude higher for sight lines $\textit{away}$ from the bulge and inner halo than for inward-looking sight lines. Although the task of identifying specific stars as ancient remains challenging, we anticipate that million-star spectral surveys and photometric surveys targeting metal-poor stars already include hundreds of stars formed before $z=5$. We predict most of these targets to have higher metallicity ($-3 < \rm [Fe/H] < -2$) than the most extreme metal-poor stars.

We investigate stellar metallicity distribution functions (MDFs), including Fe and ${\alpha}$-element abundances, in dwarf galaxies from the Feedback in Realistic Environments (FIRE) project. We examine both isolated dwarf galaxies and those that are satellites of a Milky Way-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 toward the average metallicity. This effect provides significantly better agreement with observed abundance distributions of dwarf galaxies in the Local Group, including the small intrinsic scatter in [${\alpha}$/Fe] vs. [Fe/H] (less than 0.1 dex). This small intrinsic scatter arises in our simulations because the interstellar medium (ISM) in dwarf galaxies is well-mixed at nearly all cosmic times, such that stars that form at a given time have similar abundances to within 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 kinematics of galaxies within massive clusters, as a probe of the physics of star-formation quenching within clusters. Using N-body simulations, we argue that satellite kinematics provide information about galaxy infall that is complementary to the (instantaneous) spatial distribution of satellites. Comparing the simulation results with measurements of real cluster galaxies, we find evidence that the kinematics of red (quiescent) satellite galaxies are consistent with earlier infall times than that of blue (star-forming) satellites.

We present a new set of high-resolution hydrodynamic cosmological zoom-in simulations that apply the Feedback In Realistic Environments (FIRE) physics to both Local Group (LG)-like and isolated Milky Way (MW)-like volumes (ten host systems in total with baryonic particle mass $\simeq 3,500-7,000\,M_\odot$). 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 of MW satellite dwarfs for $M_\ast \geq10^{5.5}\,M_\odot$ and predict the existence of $\sim3$ unidentified galaxies with $M_\ast\sim10^5\,M_\odot$ 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 subhalos are of secondary importance. For non-satellite galaxies, our LG-like simulations predict as many as $\sim10$ as-of-yet unseen galaxies at distances $0.3-1$ Mpc from both hosts, with $M_\ast\simeq10^{5-6}\,M_\odot$ (in halos with $V_\mathrm{max}\sim 20~\mathrm{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 $V_\mathrm{circ} \gtrsim 35~\mathrm{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.

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, SFR, and redshift in agreement with observations. In both the observed and simulated galaxies, intrinsic velocity dispersion is positively correlated with SFR. Intrinsic velocity dispersion increases with redshift out to z~1 and then flattens beyond that. In the FIRE simulations, intrinsic velocity dispersion can vary significantly on timescales of <100 Myr. These variations closely mirror the time evolution of the SFR and gas inflow rate. By cross-correlating pairs of intrinsic velocity dispersion, gas inflow rate, and SFR, we show that increased gas inflow leads to subsequent enhanced star formation, and enhancements in intrinsic velocity dispersion tend to temporally coincide with increases in gas inflow rate and SFR.

We identify subhalos in dark matter-only (DMO) zoom-in simulations that are likely to be disrupted due to baryonic effects by using a random forest classifier trained on two hydrodynamic simulations of Milky Way (MW)-mass host halos from the Latte suite of the Feedback in Realistic Environments (FIRE) project. We train our classifier using five properties of each disrupted and surviving subhalo: pericentric distance and scale factor at first pericentric passage after accretion, and scale factor, virial mass, and maximum circular velocity at accretion. Our five-property classifier identifies disrupted subhalos in the FIRE simulations with an $85\%$ out-of-bag classification score. We predict surviving subhalo populations in DMO simulations of the FIRE host halos, finding excellent agreement with the hydrodynamic results; in particular, our classifier outperforms DMO zoom-in simulations that include the gravitational potential of the central galactic disk in each hydrodynamic simulation, indicating that it captures both the dynamical effects of a central disk and additional baryonic physics. We also predict surviving subhalo populations for a suite of DMO zoom-in simulations of MW-mass host halos, finding that baryons impact each system consistently and that the predicted amount of subhalo disruption is larger than the host-to-host scatter among the subhalo populations. Although the small size and specific baryonic physics prescription of our training set limits the generality of our results, our work suggests that machine-learning classification algorithms trained on hydrodynamic zoom-in simulations can efficiently predict realistic subhalo populations.

We present the first proper motion measurement for an ultra-faint dwarf spheroidal galaxy, Segue 1, using SDSS and LBC data as the first and second epochs separated by a baseline of $\sim 10$ years. We obtain a motion of $\mu_{\alpha}\,\cos(\delta) = -0.37\pm0.57$ mas yr$^{-1}$ and $\mu_{\delta} =-3.39\pm0.58$ mas yr$^{-1}$. Combining this with the known line-of-sight velocity, this corresponds to a Galactocentric V$_\mathrm{rad}=84\pm9$ and V$_\mathrm{tan}=164^{+66}_{-55}$ km s$^{-1}$. Applying Milky Way halo masses between 0.8 to 1.6$\times 10^{12}$ M$_\odot$ results in an apocenter at 33.9$^{+21.7}_{-7.4}$ kpc and pericenter at 15.4$^{+10.1}_{-9.0}$ kpc from the Galactic center, indicating Segue~1 is rather tightly bound to the Milky Way. Since neither the orbital pole of Segue 1 nor its distance to the Milky Way is similar to the more massive classical dwarfs, it is very unlikely that Segue 1 was once a satellite of a massive known galaxy. Using cosmological zoom-in simulations of Milky Way-mass galaxies, we identify subhalos on similar orbits as Segue~1, which imply the following orbital properties: a median first infall 8.1$^{+3.6}_{-4.3}$ Gyrs ago, a median of 4 pericentric passages since then and a pericenter of 22.8$^{+4.7}_{-4.8}$ kpc. This is slightly larger than the pericenter derived directly from Segue 1 and Milky Way parameters, because galaxies with a small pericenter are more likely to be destroyed. Of the surviving subhalo analogs only 27\% were previously a satellite of a more massive dwarf galaxy (that is now destroyed), thus Segue 1 is more likely to have been accreted on its own.

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 $M_{\rm halo} \approx10^{10}\,M_{\odot}$ dark matter halos 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, $z_{\rm f}$, in both CDM and WDM models. WDM halos follow the same $M_{\star}(z=0)-V_{\rm max}(z_{\rm f})$ 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_{\star}$ WDM galaxies form all of their stars at late times. The discovery of young ultra-faint dwarf galaxies with no ancient star formation -- which do not exist in our CDM simulations -- would therefore provide evidence in support of WDM.

The shape of a galaxy's spatially unresolved, globally integrated 21-cm emission line depends on its internal gas kinematics: galaxies with rotation-supported gas disks produce double-horned profiles with steep wings, while galaxies with dispersion-supported gas produce Gaussian-like profiles with sloped wings. Using mock observations of simulated galaxies from the FIRE project, we show that one can therefore constrain a galaxy's gas kinematics from its unresolved 21-cm line profile. In particular, we find that the kurtosis of the 21-cm line increases with decreasing $V/\sigma$, and that this trend is robust across a wide range of masses, signal-to-noise ratios, and inclinations. We then quantify the shapes of 21-cm line profiles from a morphologically unbiased sample of $\sim$2000 low-redshift, HI-detected galaxies with $M_{\rm star} = 10^{7-11} M_{\odot}$ and compare to the simulated galaxies. At $M_{\rm star} \gtrsim 10^{10} M_{\odot}$, both the observed and simulated galaxies produce double-horned profiles with low kurtosis and steep wings, consistent with rotation-supported disks. Both the observed and simulated line profiles become more Gaussian-like (higher kurtosis and less-steep wings) at lower masses, indicating increased dispersion support. However, the simulated galaxies transition from rotation to dispersion support more strongly: at $M_{\rm star} = 10^{8-10}M_{\odot}$, most of the simulations produce more Gaussian-like profiles than typical observed galaxies with similar mass, indicating that gas in the low-mass simulated galaxies is, on average, overly dispersion-supported. Most of the lower-mass simulated galaxies also have somewhat lower gas fractions than the median of the observed population. The simulations nevertheless reproduce the observed line-width baryonic Tully-Fisher relation, which is insensitive to rotation vs. dispersion support.

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 hypothesis 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 ~> 1e5-6 Msun 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 timescale, t_ff <~ 3 Myr, shorter than most stellar feedback timescales. 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 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-frame UV light; sizes measured from stellar mass 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 stellar mass. The B-band light traces stellar mass 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_star~10^8 M_sun, the ratio of the half-mass radius to the halo virial radius is ~10% and does not evolve significantly at z=5-10; this ratio is typically 1-5% 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.

We investigate the merger histories of isolated dwarf galaxies based on a suite of 15 high-resolution cosmological zoom-in simulations, all with masses of $M_{\rm halo} \approx 10^{10}\,{\rm M}_{\odot}$ (and M$_\star\sim10^5-10^7\,{\rm M}_{\odot}$) at $z=0$, from the Feedback in Realistic Environments (FIRE) project. The stellar populations of these dwarf galaxies at $z=0$ are formed essentially entirely "in situ": over 90$\%$ of the stellar mass is formed in the main progenitor in all but two cases, and all 15 of the galaxies have >70$\%$ of their stellar mass formed in situ. Virtually all galaxy mergers occur prior to $z\sim3$, meaning that accreted stellar populations are ancient. On average, our simulated dwarfs undergo 5 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 next-generation facilities. The small contribution of accreted stars make it extremely difficult to discern the effects of mergers in the vast majority of dwarfs either photometrically or using resolved-star color-magnitude diagrams (CMDs). The important implication for near-field cosmology is that star formation histories 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 star formation histories of merged dwarfs.

We use hydrodynamic cosmological zoom-in simulations from the FIRE project to explore the morphologies and kinematics of fifteen Milky Way (MW)-mass galaxies. Our sample ranges from compact, bulge-dominated systems with 90% of their stellar mass within 2.5 kpc to well-ordered disks that reach $\gtrsim15$ kpc. The gas in our galaxies always forms a thin, rotation-supported disk at $z=0$, with sizes primarily determined by the gas mass. For stars, we quantify kinematics and morphology both via the fraction of stars on disk-like orbits and with the radial extent of the stellar disk. In this mass range, stellar morphology and kinematics are poorly correlated with the properties of the halo available from dark matter-only simulations (halo merger history, spin, or formation time). They more strongly correlate with the gaseous histories of the galaxies: those that maintain a high gas mass in the disk after $z\sim1$ develop well-ordered stellar disks. The best predictor of morphology we identify is the spin of the gas in the halo at the time the galaxy formed 1/2 of its stars (i.e. the gas that builds the galaxy). High-$z$ mergers, before a hot halo emerges, produce some of the most massive bulges in the sample (from compact disks in gas-rich mergers), while later-forming bulges typically originate from internal processes, as satellites are stripped of gas before the galaxies merge. Moreover, most stars in $z=0$ MW-mass galaxies (even $z=0$ bulge stars) form in a disk: $\gtrsim$60-90% of stars begin their lives rotationally supported.

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 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 $M_{\star}=10^{6-11}M_{\odot}$. Gas becomes increasingly rotationally supported with increasing galaxy mass. In the lowest-mass galaxies ($M_{\star}<10^{8}M_{\odot}$), gas fails to form a morphological disk and is primarily dispersion and pressure supported. At intermediate masses ($M_{\star}=10^{8-10}M_{\odot}$), galaxies display a wide range of gas kinematics and morphologies, from thin, rotating disks, to irregular spheroids with negligible net rotation. All the high-mass ($M_{\star}=10^{10-11}M_{\odot}$) galaxies form rotationally supported gas disks. Many of the halos whose galaxies fail to form disks harbor 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 $j_{\rm gas}/j_{\rm DM}\sim0.1$ at $M_{\star}=10^{6-7}M_{\odot}$ to $j_{\rm gas}/j_{\rm DM}\sim2$ at $M_{\star}=10^{10-11}M_{\odot}$. 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 ($M_{\star}=(10^{6}-2\times10^{7})M_{\odot}$). In our simulations, stars are uniformly less rotationally supported than gas. The common assumption that stars follow the same rotation curve as gas thus substantially overestimates galaxies' stellar angular momentum, particularly at low masses.

Dwarf galaxies are known to have remarkably low star formation efficiency due to strong feedback. Adopting the dwarf galaxies of the Milky Way as a laboratory, we explore a flexible semi-analytic galaxy formation model to understand how the feedback processes shape the satellite galaxies of the Milky Way. Using Markov-Chain Monte-Carlo, we exhaustively search a large parameter space of the model and rigorously show that the general wisdom of strong outflows as the primary feedback mechanism cannot simultaneously explain the stellar mass function and the mass--metallicity relation of the Milky Way satellites. An extended model that assumes that a fraction of baryons is prevented from collapsing into low-mass halos in the first place can be accurately constrained to simultaneously reproduce those observations. The inference suggests that two different physical mechanisms are needed to explain the two different data sets. In particular, moderate outflows with weak halo mass dependence are needed to explain the mass--metallicity relation, and prevention of baryons falling into shallow gravitational potentials of low-mass halos (e.g. "pre-heating") is needed to explain the low stellar mass fraction for a given subhalo mass.

Among the most important goals in cosmology is detecting and quantifying small ($M_{\rm halo}\simeq10^{6-9}~\mathrm{M}_\odot$) dark matter (DM) subhalos. Current probes around the Milky Way (MW) are most sensitive to such substructure within $\sim20$ kpc of the halo center, where the galaxy contributes significantly to the potential. We explore the effects of baryons on subhalo populations in $\Lambda$CDM using cosmological zoom-in baryonic simulations of MW-mass halos from the Latte simulation suite, part of the Feedback In Realistic Environments (FIRE) project. Specifically, we compare simulations of the same two halos run using (1) DM-only (DMO), (2) full baryonic physics, and (3) DM with an embedded disk potential grown to match the FIRE simulation. Relative to baryonic simulations, DMO simulations contain $\sim2\times$ as many subhalos within 100 kpc of the halo center; this excess is $\gtrsim5\times$ within 25 kpc. At $z=0$, the baryonic simulations are completely devoid of subhalos down to $3\times10^6~\mathrm{M}_\odot$ within $15$ kpc of the MW-mass galaxy, and fewer than 20 surviving subhalos have orbital pericenters <20 kpc. Despite the complexities of baryonic physics, the simple addition of an embedded central disk potential to DMO simulations reproduces this subhalo depletion, including trends with radius, remarkably well. Thus, the additional tidal field from the central galaxy is the primary cause of subhalo depletion. Subhalos on radial orbits that pass close to the central galaxy are preferentially destroyed, causing the surviving subhalo population to have tangentially biased orbits compared to DMO predictions. Our method of embedding a disk potential in DMO simulations provides a fast and accurate alternative to full baryonic simulations, thus enabling suites of cosmological simulations that can provide accurate and statistical predictions of substructure populations.

We introduce massive black holes (BHs) in the Feedback In Realistic Environments project and perform high-resolution cosmological hydrodynamic simulations of quasar-mass halos ($M_{\rm halo}(z=2) \approx 10^{12.5}\,\rm{M}_{\odot}$) 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 resolved properties of galactic nuclei. We do not include BH feedback. We show that early BH growth occurs through short ($\lesssim 1\,$Myr) 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 under-massive relative to the local $M_{\rm BH}$-$M_{\rm bulge}$ relation. 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, with BHs rapidly converging onto the scaling relation when the host reaches $M_{\rm bulge} \sim 10^{10}\,\rm{M}_{\odot}$. Our results are not sensitive to the details of the accretion model so long as BH growth is tied to the gas content within $\sim 100\,$pc of the BH. Our simulations imply that bursty stellar feedback has strong implications for BH and AGN demographics, especially in the early Universe and for low-mass galaxies.

We describe a new Large Program in progress on the Gemini North and South telescopes: Gemini Observations of Galaxies in Rich Early Environments (GOGREEN). This is an imaging and deep spectroscopic survey of 21 galaxy systems at $1<z<1.5$, selected to span a factor $>10$ in halo mass. The scientific objectives include measuring the role of environment in the evolution of low-mass galaxies, and measuring the dynamics and stellar contents of their host haloes. The targets are selected from the SpARCS, SPT, COSMOS and SXDS surveys, to be the evolutionary counterparts of today's clusters and groups. The new red-sensitive Hamamatsu detectors on GMOS, coupled with the nod-and-shuffle sky subtraction, allow simultaneous wavelength coverage over $\lambda\sim 0.6$--$1.05\mu$m, and this enables a homogeneous and statistically complete redshift survey of galaxies of all types. The spectroscopic sample targets galaxies with AB magnitudes $z^{\prime}<24.25$ and [3.6]$\mu$m$<22.5$, and is therefore statistically complete for stellar masses $M_\ast\gtrsim10^{10.3}M_\odot$, for all galaxy types and over the entire redshift range. Deep, multiwavelength imaging has been acquired over larger fields for most systems, spanning $u$ through $K$, in addition to deep IRAC imaging at 3.6$\mu$m. The spectroscopy is $\sim 50$ per cent complete as of semester 17A, and we anticipate a final sample of $\sim 500$ new cluster members. Combined with existing spectroscopy on the brighter galaxies from GCLASS, SPT and other sources, GOGREEN will be a large legacy cluster and field galaxy sample at this redshift that spectroscopically covers a wide range in stellar mass, halo mass, and clustercentric radius.

We measure how the properties of star-forming central galaxies correlate with large-scale environment, $\delta$, measured on $10$Mpc/h scales. We use group catalogs to isolate a robust sample of central galaxies with high purity and completeness. The properties we investigate are star formation rate (SFR), exponential disk scale length $R_{\rm exp}$, and Sersic index of the light profile, $n$. We find that, at all stellar masses, there is an inverse correlation between SFR and $\delta$, meaning that above-average star forming centrals live in underdense regions. For $n$ and $R_{\rm exp}$, there is no correlation with $\delta$ at $M_{\rm star}\lesssim 10^{10.5}$ $M_\odot$, but at higher masses there are positive correlations; a weak correlation with $R_{\rm exp}$ and a strong correlation with $n$. These data are evidence of assembly bias within the star-forming population. The results for SFR are consistent with a model in which SFR correlates with present-day halo accretion rate, $\dot{M}_h$. In this model, galaxies are assigned to halos using the abundance matching ansatz, which maps galaxy stellar mass onto halo mass. At fixed halo mass, SFR is assigned to galaxies using the same approach,but $\dot{M}_h$ is used to map onto SFR. The best-fit model requires some scatter in the $\dot{M}_h$-SFR relation. The $R_{\rm exp}$ and $n$ measurements are consistent with a model in which these quantities are correlated with the spin parameter of the halo, $\lambda$. Halo spin does not correlate with $\delta$ at low halo masses, but for higher mass halos, high-spin halos live in higher density environments at fixed $M_h$. Put together with the earlier installments of this series, these data demonstrate that quenching processes have limited correlation with halo formation history, but the growth of active galaxies, as well as other detailed properties, are influenced by the details of halo assembly.

We use the first Gaia data release, combined with RAVE and APOGEE spectroscopic surveys, to investigate the origin of halo stars within <~3 kpc from the Sun. We identify halo stars kinematically, as moving with a relative speed of at least 220 km/s with respect to the local standard of rest. These stars are in general more metal-poor 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 isotropic orbital distribution of the more 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, while 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 in fact 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.

Using group catalogs from the SDSS DR7, we attempt to measure galactic conformity in the local universe. We measure the quenched fraction of neighbor galaxies around isolated primary galaxies, dividing the isolated sample into star-forming and quiescent objects. We restrict our measurements to scales $>1$ Mpc to probe the correlations between the formation histories of distinct halos. Over the stellar mass range $10^{9.7} \le M_\ast/M_\odot \le 10^{10.9}$, we find minimal statistical evidence for conformity. We further compare these data to predictions of the halo age-matching model, in which the oldest galaxies are associated with the oldest halos at fixed $M_\ast$. For models with strong correlations between halo and stellar age, the conformity signal is too large to be consistent with the data. For weaker implementations of age-matching, galactic conformity is not a sensitive diagnostic of halo assembly bias, and would not produce a detectable signal in SDSS data. We reproduce the results of Kauffmann et al 2013, in which the star formation rates of neighbor galaxies are significantly reduced around primary galaxies when the primaries are themselves low star formers. However, we find this result is mainly driven by contamination in the isolation criterion, when using our group catalog to remove the small fraction of satellite galaxies in the sample, the conformity signal largely goes away. Lastly, we show that small conformity signals, i.e., 2-5% differences in the quenched fractions of neighbor galaxies, can be produced by mechanisms other than halo assembly bias. For example, if passive galaxies occupy more massive halos than star forming galaxies of the same stellar mass, a conformity signal that is consistent with recent measurements from PRIMUS (Berti et al 2016) can be produced.

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 scale heights (1.5-2.5 kpc) by z=0, while stars younger than 6 Gyr were formed in a relatively calm, stable disk. The vertical scale height 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 disk than in the inner disk (flaring). These lead to positive vertical age gradients and negative radial age gradients. The radial metallicity gradient is neg- ative at the mid-plane, flattens at larger disk height |Z|, and turns positive above |Z|~1.5kpc. 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 scale heights 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 disk in cosmological simulations even with strong stellar feedback.

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 (FIRE) 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 $\sim 20\%$ in dynamical mass estimates, even if full 3-dimensional stellar kinematics -- including the orbital anisotropy -- are known exactly. When orbital anisotropy is $\textit{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 $\textit{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.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

Central galaxies make up the majority of the galaxy population, including the majority of the quiescent population at $\mathcal{M}_* > 10^{10}\mathrm{M}_\odot$. Thus, the mechanism(s) responsible for quenching central galaxies plays a crucial role in galaxy evolution as whole. We combine a high resolution cosmological $N$-body simulation with observed evolutionary trends of the "star formation main sequence," quiescent fraction, and stellar mass function at $z < 1$ to construct a model that statistically tracks the star formation histories and quenching of central galaxies. Comparing this model to the distribution of central galaxy star formation rates in a group catalog of the SDSS Data Release 7, we constrain the timescales over which physical processes cease star formation in central galaxies. Over the stellar mass range $10^{9.5}$ to $10^{11} \mathrm{M}_\odot$ we infer quenching e-folding times that span $1.5$ to $0.5\; \mathrm{Gyr}$ with more massive central galaxies quenching faster. For $\mathcal{M}_* = 10^{10.5}\mathrm{M}_\odot$, this implies a total migration time of $\sim 4~\mathrm{Gyrs}$ from the star formation main sequence to quiescence. Compared to satellites, central galaxies take $\sim 2~\mathrm{Gyrs}$ longer to quench their star formation, suggesting that different mechanisms are responsible for quenching centrals versus satellites. Finally, the central galaxy quenching timescale we infer provides key constraints for proposed star formation quenching mechanisms. Our timescale is generally consistent with gas depletion timescales predicted by quenching through strangulation. However, the exact physical mechanism(s) responsible for this still remain unclear.

We test whether halo age and galaxy age are correlated at fixed halo and galaxy mass. The formation histories, and thus ages, of dark matter halos correlate with their large-scale density $\rho$, an effect known as assembly bias. We test whether this correlation extends to galaxies by measuring the dependence of galaxy stellar age on $\rho$. To clarify the comparison between theory and observation, and to remove the strong environmental effects on satellites, we use galaxy group catalogs to identify central galaxies and measure their quenched fraction, $f_Q$, as a function of large-scale environment. Models that match halo age to central galaxy age predict a strong positive correlation between $f_Q$ and $\rho$. However, we show that the amplitude of this effect depends on the definition of halo age: assembly bias is significantly reduced when removing the effects of splashback halos---those halos that are central but have passed through a larger halo or experienced strong tidal encounters. Defining age using halo mass at its peak value rather than current mass removes these effects. In SDSS data, at M$_{\rm gal}\gtrsim 10^{10.0}$ M_sol/h$^2$, there is a $\sim 5\%$ increase in $f_Q$ from low to high densities, which is in agreement with predictions of dark matter halos using peak halo mass. At lower stellar mass there is little to no correlation of $f_Q$ with $\rho$. For these galaxies, age-matching is inconsistent with the data across the wide range the halo formation metrics that we tested. This implies that halo formation history has a small but statistically significant impact on quenching of star formation at high masses, while the quenching process in low-mass central galaxies is uncorrelated with halo formation history.

Low-mass "dwarf" galaxies represent the most significant challenges to the cold dark matter (CDM) model of cosmological structure formation. Because these faint galaxies are (best) observed within the Local Group (LG) of the Milky Way (MW) and Andromeda (M31), understanding their formation in such an environment is critical. We present first results from the Latte Project: the Milky Way on FIRE (Feedback in Realistic Environments). This simulation models the formation of a MW-mass galaxy to z = 0 within LCDM cosmology, including dark matter, gas, and stars at unprecedented resolution: baryon particle mass of 7070 Msun with gas kernel/softening that adapts down to 1 pc (with a median of 25 - 60 pc at z = 0). Latte was simulated using the GIZMO code with a mesh-free method for accurate hydrodynamics and the FIRE-2 model for star formation and explicit feedback within a multi-phase interstellar medium. For the first time, Latte self-consistently resolves the spatial scales corresponding to half-light radii of dwarf galaxies that form around a MW-mass host down to Mstar > 10^5 Msun. Latte's population of dwarf galaxies agrees with the LG across a broad range of properties: (1) distributions of stellar masses and stellar velocity dispersions (dynamical masses), including their joint relation; (2) the mass-metallicity relation; and (3) a diverse range of star-formation histories, including their mass dependence. Thus, Latte produces a realistic population of dwarf galaxies at Mstar > 10^5 Msun that does not suffer from the "missing satellites" or "too big to fail" problems of small-scale structure formation. We conclude that baryonic physics can reconcile observed dwarf galaxies with standard LCDM cosmology.

Many properties of the Milky Way's dark matter halo, including its mass assembly history, concentration, and subhalo population, remain poorly constrained. We explore the connection between these properties of the Milky Way and its satellite galaxy population, especially the implication of the presence of the Magellanic Clouds for the properties of the Milky Way halo. Using a suite of high-resolution $N$-body simulations of Milky Way-mass halos with a fixed final Mvir ~ 10^{12.1}Msun, we find that the presence of Magellanic Cloud-like satellites strongly correlates with the assembly history, concentration, and subhalo population of the host halo, such that Milky Way-mass systems with Magellanic Clouds have lower concentration, more rapid recent accretion, and more massive subhalos than typical halos of the same mass. Using a flexible semi-analytic galaxy formation model that is tuned to reproduce the stellar mass function of the classical dwarf galaxies of the Milky Way with Markov-Chain Monte-Carlo, we show that adopting host halos with different mass-assembly histories and concentrations can lead to different best-fit models for galaxy-formation physics, especially for the strength of feedback. These biases arise because the presence of the Magellanic Clouds boosts the overall population of high-mass subhalos, thus requiring a different stellar-mass-to-halo-mass ratio to match the data. These biases also lead to significant differences in the mass--metallicity relation, the kinematics of low-mass satellites, the number counts of small satellites associated with the Magellanic Clouds, and the stellar mass of Milky Way itself. Observations of these galaxy properties can thus provide useful constraints on the properties of the Milky Way halo.

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.

Motivated by the recent discovery of several dwarfs near the Large Magellanic Cloud (LMC), we study the accretion of massive satellites onto Milky Way (MW)/M31-like halos using the ELVIS suite of N-body simulations. We identify 25 surviving LMC-mass subhalos, and investigate the lower-mass satellites that were associated with these subhalos before they fell into the MW/M31 halos. Typically, 7% of the overall z=0 satellite population of MW/M31 halos were in a surviving LMC-group before falling into the MW/M31 halo. This fraction can vary between 1% and 25%, being higher for groups with higher mass and/or more recent infall times. Groups of satellites disperse rapidly in phase space after infall, and their distances and velocities relative to the group center become statistically similar to the overall satellite population after 4-8 Gyr. We quantify the likelihood that satellites were associated with an LMC-mass group as a function of both distance and velocity relative to the LMC at z=0. The close proximity in distance of the nine Dark Energy Survey candidate dwarf galaxies to the LMC suggest that ~2-4 are likely associated with the LMC. Furthermore, if several of these dwarfs are genuine members, then the LMC-group probably fell into the MW very recently, <2 Gyr ago. If the connection with the LMC is established with follow-up velocity measurements, these "satellites of satellites" represent prime candidates to study the affects of group pre-processing on lower mass dwarfs.

In the Local Group, nearly all of the dwarf galaxies (M_star < 10^9 M_sun) that are satellites within 300 kpc (the virial radius) of the Milky Way (MW) and Andromeda (M31) have quiescent star formation and little-to-no cold gas. This contrasts strongly with comparatively isolated dwarf galaxies, which are almost all actively star-forming and gas-rich. This near dichotomy implies a rapid transformation of satellite dwarf galaxies after falling into the halos of the MW or M31. We combine the observed quiescent fractions for satellites of the MW and M31 with the infall times of satellites from the Exploring the Local Volume in Simulations (ELVIS) suite of cosmological zoom-in simulations to determine the typical timescales over which environmental processes within the MW/M31 halos remove gas and quench star formation in low-mass satellite galaxies. The quenching timescales for satellites with M_star < 10^8 M_sun are short, < 2 Gyr, and quenching is more rapid at lower M_star. These satellite quenching timescales can be 1 - 2 Gyr longer if one includes the time that satellites were environmentally preprocessed by low-mass groups prior to MW/M31 infall. We compare with quenching timescales for more massive satellites from previous works to synthesize the nature of satellite galaxy quenching across the observable range of M_star = 10^{3-11} M_sun. The satellite quenching timescale increases rapidly with satellite M_star, peaking at ~9.5 Gyr for M_star ~ 10^9 M_sun, and the timescale rapidly decreases at higher M_star to < 5 Gyr at M_star > 5 x 10^9 M_sun. Overall, galaxies with M_star ~ 10^9 M_sun, similar to the Magellanic Clouds, exhibit the longest quenching timescales, regardless of environmental or internal mechanisms.

The cosmic accretion of both dark matter and baryons into halos typically is measured using some evolving virial relation, but recent work suggests that most halo growth at late cosmic time (z < 2) is not physical but is rather the by-product of an evolving virial radius ("pseudo-evolution"). Using Omega25, a suite of cosmological simulations that incorporate both dark matter and gas dynamics with differing treatments of gas cooling, star formation, and thermal feedback, we systematically explore the physics that governs cosmic accretion into halos and their galaxies. Physically meaningful cosmic accretion of both dark matter and baryons occurs at z > 1 across our halo mass range: M_200m = 10^{11-14} M_sun. However, dark matter, because it is dissipationless, is deposited (in a time-average sense) at > R_200m(z) in a shell-like manner, such that dark-matter mass and density experience little-to-no physical growth at any radius within a halo at z < 1. In contrast, gas, because it is able to cool radiatively, experiences significant accretion at all radii, at a rate that roughly tracks the accretion rate at R_200m, at all redshifts. Infalling gas starts to decouple from dark matter at ~2 R_200m and continues to accrete to smaller radii until the onset of strong angular-momentum support at ~ 0.1 R_200m. Thus, while the growth of dark matter is subject to pseudo-evolution, the growth of baryons is not. The fact that the accretion rate of gas on galactic scales tracks the accretion rate near R_200m provides insight into the tight relations between the masses/sizes of galaxies and those of their host halos across cosmic time.

In the Local Group (LG), almost all satellite dwarf galaxies that are within the virial radius of the Milky Way (MW) and Andromeda (M31) exhibit strong environmental influence. The orbital histories of these satellites provide the key to understanding the role of the MW/M31 halo, lower-mass groups, and cosmic reionization on the evolution of dwarf galaxies. We examine the virial-infall histories of satellites with M_star = 10^{3-9} M_sun using the ELVIS suite of cosmological zoom-in dissipationless simulations of 48 MW/M31-like halos. Satellites at z = 0 fell into the MW/M31 halos typically 5 - 8 Gyr ago at z = 0.5 - 1. However, they first fell into any host halo typically 7 - 10 Gyr ago at z = 0.7 - 1.5. This difference arises because many satellites experienced "group preprocessing" in another host halo, typically of M_vir ~ 10^{10-12} M_sun, before falling into the MW/M31 halos. Satellites with lower mass and/or those closer to the MW/M31 fell in earlier and are more likely to have experienced group preprocessing; half of all satellites with M_star < 10^6 M_sun were preprocessed in a group. Infalling groups also drive most satellite-satellite mergers within the MW/M31 halos. Finally, none of the surviving satellites at z = 0 were within the virial radius of their MW/M31 halo during reionization (z > 6), and only < 4% were satellites of any other host halo during reionization. Thus, effects of cosmic reionization versus host-halo environment on the formation histories of surviving dwarf galaxies in the LG occurred at distinct epochs, separated typically by 2 - 4 Gyr, so they are separable theoretically and, in principle, observationally.

High-precision astrometry throughout the Local Group is a unique capability of the Hubble Space Telescope (HST), with potential for transformative science, including constraining the nature of dark matter, probing the epoch of reionization, and understanding key physics of galaxy evolution. While Gaia will provide unparalleled astrometric precision for bright stars in the inner halo of the Milky Way, HST is the only current mission capable of measuring accurate proper motions for systems at greater distances (> 80 kpc), which represents the vast majority of galaxies in the Local Group. The next generation of proper-motion measurements will require long time baselines, spanning many years to decades and possibly multiple telescopes, combining HST with the James Webb Space Telescope (JWST) or the Wide-Field Infrared Survey Telescope (WFIRST). However, the current HST allocation process is not conducive to such multi-cycle/multi-mission science, which will bear fruit primarily over many years. We propose an HST astrometry initiative to enable long-time-baseline, multi-mission science, which we suggest could be used to provide comprehensive kinematic measurements of all dwarf galaxies and high surface-density stellar streams in the Local Group with HST's Advanced Camera for Surveys (ACS) or Wide Field Camera 3 (WFC3). Such an initiative not only would produce forefront scientific results within the next 5 years of HST's life, but also would serve as a critical anchor point for future missions to obtain unprecedented astrometric accuracy, ensuring that HST leaves a unique and lasting legacy for decades to come.

Mergers are a common phenomenon in hierarchical structure formation, especially for massive galaxies and clusters, but their importance for dwarf galaxies in the Local Group remains poorly understood. We investigate the frequency of major mergers between dwarf galaxies in the Local Group using the ELVIS suite of cosmological zoom-in dissipationless simulations of Milky Way- and M31-like host halos. We find that ~10% of satellite dwarf galaxies with M_star > 10^6 M_sun that are within the host virial radius experienced a major merger of stellar mass ratio closer than 0.1 since z = 1, with a lower fraction for lower mass dwarf galaxies. Recent merger remnants are biased towards larger radial distance and more recent virial infall times, because most recent mergers occurred shortly before crossing within the virial radius of the host halo. Satellite-satellite mergers also occur within the host halo after virial infall, catalyzed by the large fraction of dwarf galaxies that fell in as part of a group. The merger fraction doubles for dwarf galaxies outside of the host virial radius, so the most distant dwarf galaxies in the Local Group are the most likely to have experienced a recent major merger. We discuss the implications of these results on observable dwarf merger remnants, their star formation histories, the gas content of mergers, and massive black holes in dwarf galaxies.

Galaxies that are several virial radii beyond groups/clusters show preferentially quiescent star formation rates. Using a galaxy group/cluster catalog from the Sloan Digital Sky Survey, together with a cosmological N-body simulation, we examine the origin of this environmental quenching beyond the virial radius. Accounting for the clustering of groups/clusters, we show that central galaxies show enhanced SFR quenching out to 2.5 virial radii beyond groups/clusters, and we demonstrate that this extended environmental enhancement can be explained simply by 'ejected' satellite galaxies that orbit beyond their host halo's virial radius. We show that ejected satellites typically orbit for several Gyr beyond the virial radius before falling back in, and thus they compose up to 40% of all central galaxies near groups/clusters. We show that a model in which ejected satellites experience the same SFR quenching as satellites within a host halo can explain essentially all environmental dependence of galaxy quenching. Furthermore, ejected satellites (continue to) lose significant halo mass, an effect that is potentially observable via gravitational lensing. The SFRs/colors and stellar-to-halo masses of ejected satellites highlight the importance of environmental history and present challenges to models of galaxy occupation that ignore such history.

We investigate the use of the halo mass-gap statistic --- defined as the logarithmic difference in mass between the host halo and its most massive satellite subhalo --- as a probe of halo age and concentration. A cosmological N-body simulation is used to study N ~25, 000 group/cluster sized halos in the mass range 10^12.5 < M_halo/M_sol < 10^14.5. In agreement with previous work, we find that halo mass-gap is related to halo formation time and concentration. On average, older and more highly concentrated halos have larger halo mass-gaps, and this trend is stronger than the mass-concentration relation over a similar dynamic range. However, there is a large amount of scatter owing to the transitory nature of the satellite subhalo population, which limits the use of the halo mass-gap statistic on an object-by-object basis. For example, we find that 20% of very large halo mass-gap systems (akin to "fossil groups") are young, and have likely experienced a recent merger between a massive satellite subhalo and the central subhalo. We relate halo mass-gap to the observable stellar mass-gap via abundance matching. Using a galaxy group catalog constructed from the Sloan Digital Sky Survey Data Release 7, we find that the star formation and structural properties of galaxies at fixed mass show no trend with stellar mass-gap. This is despite a variation in halo-age of ~2.5 Gyr over ~1.2 dex in stellar mass-gap. Thus, we find no evidence to suggest that the halo formation history significantly affects galaxy properties.

Satellite galaxies in groups and clusters are more likely to have low star formation rates (SFR) and lie on the red-sequence than central (field) galaxies. Using galaxy group/cluster catalogs from SDSS DR7, together with a cosmological N-body simulation to track satellite orbits, we examine the star formation histories and quenching timescales of satellites of M_star > 5 x 10^9 M_sun at z=0. We first explore satellite infall histories: group preprocessing and ejected orbits are critical aspects of satellite evolution, and properly accounting for these, satellite infall typically occurred at z~0.5, or ~5 Gyr ago. To obtain accurate initial conditions for the SFRs of satellites at their time of first infall, we construct an empirical parametrization for the evolution of central galaxy SFRs and quiescent fractions. With this, we constrain the importance and efficiency of satellite quenching as a function of satellite and host halo mass, finding that satellite quenching is the dominant process for building up all quiescent galaxies at M_star < 10^10 M_sun. We then constrain satellite star formation histories, finding a 'delayed-then-rapid' quenching scenario: satellite SFRs evolve unaffected for 2-4 Gyr after infall, after which star formation quenches rapidly, with an e-folding time of < 0.8 Gyr. These quenching timescales are shorter for more massive satellites but do not depend on host halo mass: the observed increase in satellite quiescent fraction with halo mass arises simply because of satellites quenching in a lower mass group prior to infall (group preprocessing), which is responsible for up to half of quenched satellites in massive clusters. Because of the long time delay before quenching starts, satellites experience significant stellar mass growth after infall, nearly identical to central galaxies. This fact provides key physical insight into the subhalo abundance matching method.

We report on the small scale (0.5<r<40h^-1 Mpc) clustering of 78895 massive (M*~10^11.3M_sun) galaxies at 0.2<z<0.4 from the first two years of data from the Baryon Oscillation Spectroscopic Survey (BOSS), to be released as part of SDSS Data Release 9 (DR9). We describe the sample selection, basic properties of the galaxies, and caveats for working with the data. We calculate the real- and redshift-space two-point correlation functions of these galaxies, fit these measurements using Halo Occupation Distribution (HOD) modeling within dark matter cosmological simulations, and estimate the errors using mock catalogs. These galaxies lie in massive halos, with a mean halo mass of 5.2x10^13 h^-1 M_sun, a large scale bias of ~2.0, and a satellite fraction of 12+/-2%. Thus, these galaxies occupy halos with average masses in between those of the higher redshift BOSS CMASS sample and the original SDSS I/II LRG sample.

Using galaxy group/cluster catalogs created from the Sloan Digital Sky Survey Data Release 7, we examine in detail the specific star formation rate (SSFR) distribution of satellite galaxies and its dependence on stellar mass, host halo mass, and halo-centric radius. All galaxies, regardless of central-satellite designation, exhibit a similar bimodal SSFR distribution, with a strong break at SSFR ~ 10^-11 yr^-1 and the same high SSFR peak; in no regime is there ever an excess of galaxies in the `green valley'. Satellite galaxies are simply more likely to lie on the quenched (`red sequence') side of the SSFR distribution. Furthermore, the satellite quenched fraction excess above the field galaxy value is nearly independent of galaxy stellar mass. An enhanced quenched fraction for satellites persists in groups with halo masses down to 3 x 10^11 Msol and increases strongly with halo mass and toward halo center. We find no detectable quenching enhancement for galaxies beyond ~2R_vir around massive clusters once these galaxies have been decomposed into centrals and satellites. These trends imply that (1) galaxies experience no significant environmental effects until they cross within ~R_vir of a more massive host halo, (2) after this, star formation in active satellites continues to evolve in the same manner as active central galaxies for several Gyrs, and (3) once begun, satellite star formation quenching occurs rapidly. These results place strong constraints on satellite-specific quenching mechanisms, as we will discuss further in companion papers.

The properties of dark matter halos, including mass growth, correlate with larger scale environment at fixed mass, an effect known as assembly bias. However, whether this environmental dependence manifests itself in galaxy properties remains unclear. We apply a group-finding algorithm to DR7 of the SDSS to estimate the halo mass of each galaxy and to decompose galaxies into those that exist at the centers of distinct halos and those that orbit as satellites within larger halos. Using the 4000-A break as a measure of star formation history, we examine the correlation between the quenched fraction of galaxies, f_q, and large-scale environment, rho. At all galaxy magnitudes, there is a positive, monotonic relationship between f_q and rho. We use the group catalog to decompose this correlation into the contribution from central and satellite galaxies as a function of halo mass. Because satellites are more likely to be quenched than central galaxies, the observed f_q-rho correlation is primarily due to variations of the halo mass function with environment, which causes a larger fraction of satellite galaxies at high rho. For low-mass central galaxies (Mgal <~ 10^10.0 Msol/h^2), there is no correlation between f_q and rho. These results are inconsistent with the strong assembly bias of dark matter halos seen in this mass regime if recent galaxy growth at all correlates with recent halo growth, as we demonstrate through a high resolution N-body simulation. We also find that the mean stellar age of quenched central galaxies is independent of rho at fixed Mgal, while the formation times of low mass halos vary significantly. We conclude that the processes that halt the star formation of low mass central galaxies are not correlated to the formation histories of their host halos, and old galaxies do not reside preferentially in old halos. (Abridged)

The orbital properties of infalling satellite halos set the initial conditions which control the subsequent evolution of subhalos and the galaxies that they host, with implications for mass stripping, star formation quenching, and merging. Using a high-resolution, cosmological N-body simulation, I examine the orbital parameters of satellite halos as they merge with larger host halos, focusing primarily on orbital circularity and pericenter. I explore in detail how these orbital parameters depend on mass and redshift. Satellite orbits become more radial and plunge deeper into their host halo at higher host halo mass, but they do not significantly depend on satellite halo mass. Additionally, satellite orbits become more radial and plunge deeper into their host halos at higher redshift. I also examine satellite velocities, finding that most satellites infall with less specific angular momentum than the host halo virial value, but that satellites are `hotter' than the host virial velocity. I discuss the implications of these results to the processes of galaxy formation and evolution, and I provide fitting formulas to the mass and redshift dependence of satellite orbital circularity and pericenter.

Different methodologies lead to order-of-magnitude variations in predicted galaxy merger rates. We examine and quantify the dominant uncertainties. Different halo merger rates and subhalo 'destruction' rates agree to within a factor ~2 given proper care in definitions. If however (sub)halo masses are not appropriately defined or are under-resolved, the major merger rate can be dramatically suppressed. The dominant differences in galaxy merger rates owe to baryonic physics. Hydrodynamic simulations without feedback and older models that do not agree with the observed galaxy mass function propagate factor ~5 bias in the resulting merger rates. However, if the model matches the galaxy mass function, properties of central galaxies are sufficiently converged to give small differences in merger rates. But variations in baryonic physics of satellites also have dramatic effects. The known problem of satellite 'over-quenching' in most semi-analytic models (SAMs), whereby SAM satellites are too efficiently stripped of gas, could lead to order-of-magnitude under-estimates of merger rates for low-mass, gas-rich galaxies. Fixing the satellite properties to observations tends to predict higher merger rates, but with factor ~2 empirical uncertainties. Choice of mass ratio definition matters: at low masses, most true major mergers (in baryonic/dynamical galaxy mass) will appear to be minor mergers in their stellar or luminosity mass ratio. Observations and models using these criteria may underestimate major merger rates by factors ~5. Orbital parameters and gas fractions also introduce factor ~3 differences in amount of bulge formed by mergers, even for fixed mass ratio encounters.

We use a suite of semi-empirical models to predict galaxy merger rates and contributions to bulge growth as functions of merger mass, redshift, and mass ratio. The models use empirical halo occupation constraints to identify mergers, together with high-resolution simulations to quantify how mergers with different properties contribute to the bulge population. We find good agreement with a variety of observational constraints, and provide fitting functions for merger rates and contributions to bulge growth. We identify several robust conclusions. (1) Major mergers dominate formation and assembly of L* bulges and the spheroid mass density, minor mergers contribute ~30%. (2) This is mass-dependent: bulge formation is dominated by more minor mergers in lower-mass systems. At higher masses, bulges form in major mergers near L*, but subsequently assemble in minor mergers. (3) The minor/major contribution is also morphology-dependent: higher B/T systems form in more major mergers, lower B/T systems form in situ from minor mergers. (4) Low-mass galaxies, being gas-rich, require more major mergers to reach the same B/T as high-mass systems. (5) Absolute merger rates increase with galaxy mass. (6) Predicted rates agree well with observations, but suggest that some morphology-selected samples include contamination from minor mergers. (7) Predicted rates agree with integrated growth in bulge mass with cosmic time, but with factor ~2 uncertainty - half the bulge mass density could come from non-mergers. We consider ~1000 model variations and quantify resulting uncertainties. Conclusions regarding the major/minor contribution to bulge growth are very robust, absolute merger rates have systematic factor ~2 uncertainties.

In hierarchical structure formation, dark matter halos that merge with larger halos can persist as subhalos. These subhalos are likely hosts of visible galaxies. While the dense halo environment rapidly strips subhalos of their dark mass, the compact luminous material can remain intact for some time, making the correspondence of galaxies with severely stripped subhalos unclear. Many galaxy evolution models assume that satellite galaxies eventually merge with their central galaxy, but this ignores the possibility of satellite tidal disruption. We use a high-resolution $N$-body simulation of cosmological volume to explore satellite galaxy merging/disruption criteria based on dark matter subhalo dynamics. We explore the impact that satellite merging/disruption has on the Halo Occupation Distribution and radial profile of the remnants. Using abundance matching to assign stellar mass/luminosity to subhalos, we compare with observed galaxy clustering, satellite fractions, cluster satellite luminosity functions, finding that subhalos reproduce well these observables on all scales. Our results imply that satellite subhalos corresponding to >0.2L* galaxies must be well-resolved down to 1-3% of their mass at infall to robustly trace the galaxy population. We also explore a simple analytic model based on dynamical friction for satellite galaxy infall, finding good agreement with our subhalo catalog and observations.

We analyze the clustering of red and blue galaxies from four samples spanning a redshift range of 0.4<z<2.0 to test the various scenarios by which galaxies evolve onto the red sequence. The data are taken from the UKIDSS Ultra Deep Survey, DEEP2, and COMBO-17. The use of clustering allows us to determine what fraction of the red sequence is made up of central galaxies and satellite galaxies. At all redshifts, including z=0, the data are consistent with ~60% of satellite galaxies being red or quenched, implying that ~1/3 of the red sequence is comprised of satellite galaxies. More than three-fourths of red satellite galaxies were moved to the red sequence after they were accreted onto a larger halo. The constant fraction of satellite galaxies that are red yields a quenching time for satellite galaxies that depends on redshift in the same way as halo dynamical times; t_Q ~ (1+z)^{-1.5}. In three of the four samples, the data favor a model in which red central galaxies are a random sample of all central galaxies; there is no preferred halo mass scale at which galaxies make the transition from star-forming to red and dead. The large errors on the fourth sample inhibit any conclusions. Theoretical models in which star formation is quenched above a critical halo mass are excluded by these data. A scenario in which mergers create red central galaxies imparts a weaker correlation between halo mass and central galaxy color, but even the merger scenario creates tension with red galaxy clustering at redshifts above 0.5. These results suggest that the mechanism by which central galaxies become red evolves from z=0.5 to z=0.

While both society and astronomy have evolved greatly over the past fifty years, the academic institutions and incentives that shape our field have remained largely stagnant. As a result, the astronomical community is faced with several major challenges, including: (1) the training that we provide does not align with the skills that future astronomers will need, (2) the postdoctoral phase is becoming increasingly demanding and demoralizing, and (3) our jobs are increasingly unfriendly to families with children. Solving these problems will require conscious engineering of our profession. Fortunately, this Decadal Review offers the opportunity to revise outmoded practices to be more effective and equitable. The highest priority of the Subcommittee on the State of the Profession should be to recommend specific, funded activities that will ensure the field meets the challenges we describe.

High-resolution simulations of cosmological structure formation indicate that dark matter substructure in dense environments, like groups and clusters, may survive for a long time. These dark matter subhalos are the likely hosts of galaxies. We examine the small-scale spatial clustering of subhalo major mergers at high redshift using high-resolution N-body simulations of cosmological volumes. Recently merged, massive subhalos exhibit enhanced clustering on scales ~100-300 kpc/h, relative to all subhalos of the same infall mass, for a short time after a major merger (< 500 Myr). The small-scale clustering enhancement is smaller for lower mass subhalos, which also show a deficit on scales just beyond the excess. Halos hosting recent subhalo mergers tend to have more subhalos; for massive subhalos the excess is stronger and it tends to increase for the most massive host halos. The subhalo merger fraction is independent of halo mass for the scales we probe. In terms of satellite and central subhalos, the merger increase in small-scale clustering for massive subhalos arises from recently merged massive central subhalos having an enhanced satellite population. Our mergers are defined via their parent infall mass ratios. Subhalos experiencing major mass gains also exhibit a small-scale clustering enhancement, but these correspond to two-body interactions leading to two final subhalos, rather than subhalo coalescence.

Galaxies are believed to be in one-to-one correspondence with simulated dark matter subhalos. We use high-resolution N-body simulations of cosmological volumes to calculate the statistical properties of subhalo (galaxy) major mergers at high redshift (z=0.6-5). We measure the evolution of the galaxy merger rate, finding that it is much shallower than the merger rate of dark matter host halos at z>2.5, but roughly parallels that of halos at z<1.6. We also track the detailed merger histories of individual galaxies and measure the likelihood of multiple mergers per halo or subhalo. We examine satellite merger statistics in detail: 15%-35% of all recently merged galaxies are satellites and satellites are twice as likely as centrals to have had a recent major merger. Finally, we show how the differing evolution of the merger rates of halos and galaxies leads to the evolution of the average satellite occupation per halo, noting that for a fixed halo mass, the satellite halo occupation peaks at z~2.5.

Galaxy cluster merger statistics are an important component in understanding the formation of large-scale structure. Unfortunately, it is difficult to study merger properties and evolution directly because the identification of cluster mergers in observations is problematic. We use large N-body simulations to study the statistical properties of massive halo mergers, specifically investigating the utility of close halo pairs as proxies for mergers. We examine the relationship between pairs and mergers for a wide range of merger timescales, halo masses, and redshifts (0<z<1). We also quantify the utility of pairs in measuring merger bias. While pairs at very small separations will reliably merge, these constitute a small fraction of the total merger population. Thus, pairs do not provide a reliable direct proxy to the total merger population. We do find an intriguing universality in the relation between close pairs and mergers, which in principle could allow for an estimate of the statistical merger rate from the pair fraction within a scaled separation, but including the effects of redshift space distortions strongly degrades this relation. We find similar behavior for galaxy-mass halos, making our results applicable to field galaxy mergers at high redshift. We investigate how the halo merger rate can be statistically described by the halo mass function via the merger kernel (coagulation), finding an interesting environmental dependence of merging: halos within the mass resolution of our simulations merge less efficiently in overdense environments. Specifically, halo pairs with separations less than a few Mpc/h are more likely to merge in underdense environments; at larger separations, pairs are more likely to merge in overdense environments.

The clustering properties of dark matter halos are a firm prediction of modern theories of structure formation. We use two large volume, high-resolution N-body simulations to study how the correlation function of massive dark matter halos depends upon their mass and formation history. We find that halos with the lowest concentrations are presently more clustered than those of higher concentration, the size of the effect increasing with halo mass; this agrees with trends found in studies of lower mass halos. The clustering dependence on other characterizations of the full mass accretion history appears weaker than the effect with concentration. Using the integrated correlation function, marked correlation functions, and a power-law fit to the correlation function, we find evidence that halos which have recently undergone a major merger or a large mass gain have slightly enhanced clustering relative to a randomly chosen population with the same mass distribution.