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.

It has been established for decades that rotation curves deviate from the Newtonian gravity expectation given baryons alone below a characteristic acceleration scale $g_{\dagger}\sim 10^{-8}\,\rm{cm\,s^{-2}}$, a scale promoted to a new fundamental constant in MOND. In recent years, theoretical and observational studies have shown that the star formation efficiency (SFE) of dense gas scales with surface density, SFE $\sim \Sigma/\Sigma_{\rm crit}$ with $\Sigma_{\rm crit} \sim \langle\dot{p}/m_{\ast}\rangle/(\pi\,G)\sim 1000\,\rm{M_{\odot}\,pc^{-2}}$ (where $\langle \dot{p}/m_{\ast}\rangle$ is the momentum flux output by stellar feedback per unit stellar mass in a young stellar population). We argue that the SFE, more generally, should scale with the local gravitational acceleration, i.e. that SFE $\sim g_{\rm tot}/g_{\rm crit} \equiv (G\,M_{\rm tot}/R^{2}) / \langle\dot{p}/m_{\ast}\rangle$, where $M_{\rm tot}$ is the total gravitating mass and $g_{\rm crit}=\langle\dot{p}/m_{\ast}\rangle = \pi\,G\,\Sigma_{\rm crit} \approx 10^{-8}\,\rm{cm\,s^{-2}} \approx g_{\dagger}$. Hence the observed $g_\dagger$ may correspond to the characteristic acceleration scale above which stellar feedback cannot prevent efficient star formation, and baryons will eventually come to dominate. We further show how this may give rise to the observed acceleration scaling $g_{\rm obs}\sim(g_{\rm baryon}\,g_{\dagger})^{1/2}$ (where $g_{\rm baryon}$ is the acceleration due to baryons alone) and flat rotation curves. The derived characteristic acceleration $g_{\dagger}$ can be expressed in terms of fundamental constants (gravitational constant, proton mass, and Thomson cross section): $g_{\dagger}\sim 0.1\,G\,m_{p}/\sigma_{\rm T}$.

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.

We perform a series of controlled N-body simulations to study realizations of the recently discovered Antlia 2 galaxy in cold dark matter (CDM) and self-interacting dark matter (SIDM) scenarios. Our simulations contain six benchmark models, where we vary the initial halo concentration and the self-scattering cross section. We adopt well-motivated initial stellar and halo masses, and our fiducial orbit has a small pericenter. After evolving in the Milky Way's tidal field, the simulated galaxies experience significant mass loss and their stellar distributions expand accordingly. These tidal effects are more prominent if the initial halo concentration is lower and if the self-scattering cross section is larger. Our results show that Antlia 2-like galaxies could be realized in CDM if the halo concentration is low and the stellar distribution is diffuse at the infall time, while these conditions could be relaxed in SIDM. We also find all the simulated galaxies predict approximately the same stellar velocity dispersion after imposing selection criteria for stellar particles. This has important implications for testing dark matter models using tidally disturbed systems.

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.

With the release of Gaia DR2, it is now possible to measure the proper motions (PMs) of the lowest mass, ultra-faint satellite galaxies in the Milky Way's (MW) halo for the first time. Many of these faint satellites are posited to have been accreted as satellites of the Magellanic Clouds (MCs). Using their 6-dimensional phase space information, we calculate the orbital histories of 13 ultra-faint satellites and five classical dwarf spheroidals in a combined MW+LMC+SMC potential to determine which galaxies are dynamically associated with the MCs. These 18 galaxies are separated into four classes: i.) long-term Magellanic satellites that have been bound to the MCs for at least the last two consecutive orbits around the MCs (Carina 2, Carina 3, Horologium 1, Hydrus 1); ii.) Magellanic satellites that were recently captured by the MCs $<$ 1 Gyr ago (Reticulum 2, Phoenix 2); iii.) MW satellites that have interacted with the MCs (Sculptor 1, Tucana 3, Segue 1); and iv.) MW satellites (Aquarius 2, Canes Venatici 2, Crater 2, Draco 1, Draco 2, Hydra 2, Carina, Fornax, Ursa Minor). Results are reported for a range of MW and LMC masses. Contrary to previous work, we find no dynamical association between Carina, Fornax, and the MCs. Finally, we determine that the addition of the SMC's gravitational potential affects the longevity of satellites as members of the Magellanic system (long-term versus recently captured), but it does not change the total number of Magellanic satellites.

CMB-HD is a proposed ultra-deep (0.5 uk-arcmin), high-resolution (15 arcseconds) millimeter-wave survey over half the sky that would answer many outstanding questions in both fundamental physics of the Universe and astrophysics. This survey would be delivered in 7.5 years of observing 20,000 square degrees, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.

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.

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.

Bose-Einstein Condensate Dark Matter (BECDM; also known as Fuzzy Dark Matter) is motivated by fundamental physics and has recently received significant attention as a serious alternative to the established Cold Dark Matter (CDM) model. We perform cosmological simulations of BECDM gravitationally coupled to baryons and investigate structure formation at high redshifts ($z \gtrsim 5$) for a boson mass $m=2.5\cdot 10^{-22}~{\rm eV}$, exploring the dynamical effects of its wavelike nature on the cosmic web and the formation of first galaxies. Our BECDM simulations are directly compared to CDM as well as to simulations where the dynamical quantum potential is ignored and only the initial suppression of the power spectrum is considered -- a Warm Dark Matter-like ("WDM") model often used as a proxy for BECDM. Our simulations confirm that "WDM" is a good approximation to BECDM on large cosmological scales even in the presence of the baryonic feedback. Similarities also exist on small scales, with primordial star formation happening both in isolated haloes and continuously along cosmic filaments; the latter effect is not present in CDM. Global star formation and metal enrichment in these first galaxies are delayed in BECDM/"WDM" compared to the CDM case: in BECDM/"WDM" first stars form at $z\sim 13$/$13.5$ while in CDM star formation starts at $z\sim 35$. The signature of BECDM interference, not present in "WDM", is seen in the evolved dark matter power spectrum: although the small scale structure is initially suppressed, power on kpc scales is added at lower redshifts. Our simulations lay the groundwork for realistic simulations of galaxy formation in BECDM.

In hierarchical models of structure formation, the first galaxies form in low-mass dark matter potential wells, probing the behavior of dark matter on kiloparsec (kpc) scales. Even though these objects are below the detection threshold of current telescopes, future missions will open an observational window into this emergent world. In this Letter we investigate how the first galaxies are assembled in a `fuzzy' dark matter (FDM) cosmology where dark matter is an ultralight $\sim 10^{-22}$~eV boson and the primordial stars are expected to form along dense dark matter filaments. Using a first-of-its-kind cosmological hydrodynamical simulation, we explore the interplay between baryonic physics and unique wavelike features inherent to FDM. In our simulation, the dark matter filaments show coherent interference patterns on the boson de Broglie scale and develop cylindrical soliton-like cores which are unstable under gravity and collapse into kpc-scale spherical solitons. Features of the dark matter distribution are largely unaffected by the baryonic feedback. On the contrary, the distributions of gas and stars, which do form along the entire filament, exhibit central cores imprinted by dark matter -- a smoking gun signature of FDM.

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.

The dark energy plus cold dark matter ($\Lambda$CDM) cosmological model has been a demonstrably successful framework for predicting and explaining the large-scale structure of Universe and its evolution with time. Yet on length scales smaller than $\sim 1$ Mpc and mass scales smaller than $\sim 10^{11} M_{\odot}$, the theory faces a number of challenges. For example, the observed cores of many dark-matter dominated galaxies are both less dense and less cuspy than naively predicted in $\Lambda$CDM. The number of small galaxies and dwarf satellites in the Local Group is also far below the predicted count of low-mass dark matter halos and subhalos within similar volumes. These issues underlie the most well-documented problems with $\Lambda$CDM: Cusp/Core, Missing Satellites, and Too-Big-to-Fail. The key question is whether a better understanding of baryon physics, dark matter physics, or both will be required to meet these challenges. Other anomalies, including the observed planar and orbital configurations of Local Group satellites and the tight baryonic/dark matter scaling relations obeyed by the galaxy population, have been less thoroughly explored in the context of $\Lambda$CDM theory. Future surveys to discover faint, distant dwarf galaxies and to precisely measure their masses and density structure hold promising avenues for testing possible solutions to the small-scale challenges going forward. Observational programs to constrain or discover and characterize the number of truly dark low-mass halos are among the most important, and achievable, goals in this field over then next decade. These efforts will either further verify the $\Lambda$CDM paradigm or demand a substantial revision in our understanding of the nature of dark matter.

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.

A millimeter-wave survey over half the sky, that spans frequencies in the range of 30 to 350 GHz, and that is both an order of magnitude deeper and of higher-resolution than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. By providing such a deep, high-resolution millimeter-wave survey (about 0.5 uK-arcmin noise and 15 arcsecond resolution at 150 GHz), CMB-HD will enable major advances. It will allow 1) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10/hMpc), which probes dark matter particle properties. It will also allow 2) measurements of the thermal and kinetic Sunyaev-Zel'dovich effects on small scales to map the gas density and gas pressure profiles of halos over a wide field, which probes galaxy evolution and cluster astrophysics. In addition, CMB-HD would allow us to cross critical thresholds in fundamental physics: 3) ruling out or detecting any new, light (< 0.1eV), thermal particles, which could potentially be the dark matter, and 4) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe. Such a survey would also 5) monitor the transient sky by mapping the full observing region every few days, which opens a new window on gamma-ray bursts, novae, fast radio bursts, and variable active galactic nuclei. Moreover, CMB-HD would 6) provide a census of planets, dwarf planets, and asteroids in the outer Solar System, and 7) enable the detection of exo-Oort clouds around other solar systems, shedding light on planet formation. CMB-HD will deliver this survey in 5 years of observing half the sky, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. The telescopes will field about 2.4 million detectors (600,000 pixels) in total.

Observations of low-mass satellite galaxies in the nearby Universe point towards a strong dichotomy in their star-forming properties relative to systems with similar mass in the field. Specifically, satellite galaxies are preferentially gas poor and no longer forming stars, while their field counterparts are largely gas rich and actively forming stars. Much of the recent work to understand this dichotomy has been statistical in nature, determining not just that environmental processes are most likely responsible for quenching these low-mass systems but also that they must operate very quickly after infall onto the host system, with quenching timescales $\lesssim 2~ {\rm Gyr}$ at ${M}_{\star} \lesssim 10^{8}~{\rm M}_{\odot}$. This work utilizes the newly-available $Gaia$ DR2 proper motion measurements along with the Phat ELVIS suite of high-resolution, cosmological, zoom-in simulations to study low-mass satellite quenching around the Milky Way on an object-by-object basis. We derive constraints on the infall times for $37$ of the known low-mass satellite galaxies of the Milky Way, finding that $\gtrsim~70\%$ of the `classical' satellites of the Milky Way are consistent with the very short quenching timescales inferred from the total population in previous works. The remaining classical Milky Way satellites have quenching timescales noticeably longer, with $\tau_{\rm quench} \sim 6 - 8~{\rm Gyr}$, highlighting how detailed orbital modeling is likely necessary to understand the specifics of environmental quenching for individual satellite galaxies. Additionally, we find that the $6$ ultra-faint dwarf galaxies with publicly available $HST$-based star-formation histories are all consistent with having their star formation shut down prior to infall onto the Milky Way -- which, combined with their very early quenching times, strongly favors quenching driven by reionization.

We introduce an extension of the ELVIS project to account for the effects of the Milky Way galaxy on its subhalo population. Our simulation suite, Phat ELVIS, consists of twelve high-resolution cosmological dark matter-only (DMO) zoom simulations of Milky Way-size $\Lambda$CDM~ haloes ($M_{\rm v} = 0.7-2 \times 10^{12} \,\mathrm{M}_\odot$) along with twelve re-runs with embedded galaxy potentials grown to match the observed Milky Way disk and bulge today. The central galaxy potential destroys subhalos on orbits with small pericenters in every halo, regardless of the ratio of galaxy mass to halo mass. This has several important implications. 1) Most of the $\mathtt{Disk}$ runs have no subhaloes larger than $V_{\rm max} = 4.5$ km s$^{-1}$ within $20$ kpc and a significant lack of substructure going back $\sim 8$ Gyr, suggesting that local stream-heating signals from dark substructure will be rare. 2) The pericenter distributions of Milky Way satellites derived from $\mathit{Gaia}$ data are remarkably similar to the pericenter distributions of subhaloes in the $\mathtt{Disk}$ runs, while the DMO runs drastically over-predict galaxies with pericenters smaller than 20 kpc. 3) The enhanced destruction produces a tension opposite to that of the classic `missing satellites' problem: in order to account for ultra-faint galaxies known within $30$ kpc of the Galaxy, we must populate haloes with $V_\mathrm{peak} \simeq 7$ km s$^{-1}$ ($M \simeq 3 \times 10^{7} \,\mathrm{M}_\odot$ at infall), well below the atomic cooling limit of $V_\mathrm{peak} \simeq 16$ km s$^{-1}$ ($M \simeq 5 \times 10^{8} \,\mathrm{M}_\odot$ at infall). 4) If such tiny haloes do host ultra-faint dwarfs, this implies the existence of $\sim 1000$ satellite galaxies within 300 kpc of the Milky Way.

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.

The premise of "near-field cosmology" is that the study of nearby low-mass galaxies on a star-by-star basis has implications that extend far beyond the local Universe and include the nature of dark matter, lives and deaths of the first stars, cosmic reionization, and galaxy formation across cosmic time. In response to the Astro2020 call for science white papers, we provide a high-level synopsis of several major questions in near-field cosmology, the potential for transformative progress in the next decade and beyond, and mechanisms for achieving these goals. Our recommendations include: (i) obtaining kinematic information from large samples of individual line-of-sight and transverse velocities of stars and from spatially resolved gas observations in nearby dwarf galaxies to directly inform our understanding of dark matter interactions in the deeply nonlinear regime; and (ii) an intensive program of optical imaging and UV and optical spectroscopy of resolved stars in low-mass galaxies, including low-metallicity massive stars and HII regions (which are central to understanding faint galaxies at all earlier epochs). Realizing these aims requires a large aperture ($\gtrsim 9$m) UV/optical space telescope, a two-fold community investment in ground-based optical spectroscopy through both ELTs and a dedicated wide-field spectroscopic survey on a 10m class telescope, and radio facilities that can measure high-precision, spatially resolved rotation curves of low-mass galaxies. Finally, we suggest that a servicing mission to extend the life of HST beyond 2025 should be seriously explored.

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 contrast predictions for the high-redshift galaxy population and reionization history between cold dark matter (CDM) and an alternative self-interacting dark matter model based on the recently developed ETHOS framework that alleviates the small-scale CDM challenges within the Local Group. We perform the highest resolution hydrodynamical cosmological simulations (a 36~Mpc$^3$ volume with gas cell mass of $\sim10^5\mathrm{M}_{\odot}$ and minimum gas softening of $\sim180$~pc) within ETHOS to date -- plus a CDM counterpart -- to quantify the abundance of galaxies at high redshift and their impact on reionization. We find that ETHOS predicts galaxies with higher ultraviolet (UV) luminosities than their CDM counterparts and a faster build-up of the faint end of the UV luminosity function. These effects, however, make the optical depth to reionization less sensitive to the power spectrum cut-off: the ETHOS model differs from the CDM $\tau$ value by only 10 per cent and is consistent with Planck limits if the effective escape fraction of UV photons is 0.1-0.5. We conclude that current observations of high-redshift luminosity functions cannot differentiate between ETHOS and CDM models, but deep JWST surveys of strongly-lensed, inherently faint galaxies have the potential to test non-CDM models that offer attractive solutions to CDM's Local Group problems.

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 predominantly ancient stellar populations observed in the lowest-mass galaxies (i.e. ultra-faint dwarfs) suggest that their star formation was suppressed by reionization. Most of the well-studied ultra-faint dwarfs, however, are within the central half of the Milky Way dark matter halo, such that they are consistent with a population that was accreted at early times and thus potentially quenched via environmental processes. To study the potential role of environment in suppressing star formation on the smallest scales, we utilize the Exploring the Local Volume in Simulations (ELVIS) suite of $N$-body simulations to constrain the distribution of infall times for low-mass subhalos likely to host the ultra-faint population. For the ultra-faint satellites of the Milky Way with star-formation histories inferred from $Hubble~Space~Telescope$ imaging, we find that environment is highly unlikely to play a dominant role in quenching their star formation. Even when including the potential effects of pre-processing, there is a $\lesssim 0.1\%$ probability that environmental processes quenched all of the known ultra-faint dwarfs early enough to explain their observed star-formation histories. Instead, we argue for a mass floor in the effectiveness of satellite quenching at roughly $M_{\star} \sim 10^{5}~M_{\odot}$, below which star formation in surviving galaxies is globally suppressed by reionization. We predict a large population of quenched ultra-faint dwarfs in the Local Field ($1 < R/R_{\rm{vir}} < 2$), with as many as $\sim250$ to be discovered by future wide-field imaging surveys.

We investigate the spatial distribution of Lyman alpha emitting galaxies (LAEs) at $z \approx 2.67$, selected from the NOAO Deep Wide-Field Survey (NDWFS), using two-point statistics and topological diagnostics adopted from network science. We measure the clustering length, $r_0 \approx 4 h^{-1}$ Mpc, and the bias, $b_{LAE} = 2.2^{+0.2}_{-0.1}$. Fitting the clustering with halo occupation distribution (HOD) models results in two disparate possibilities: (1) where the fraction of central galaxies is $<$1% in halos of mass $>10^{12}$$M_\odot$; and (2) where the fraction is $\approx$20%. We refer to these two scenarios as the `Dusty Core Scenario' for Model#1 since most of central galaxies in massive halos are dead in Ly$\alpha$ emission, and the `Pristine Core Scenario' for Model#2 since the central galaxies are bright in Ly$\alpha$ emission. Traditional two-point statistics cannot distinguish between these disparate models given the current data sets. To overcome this degeneracy, we generate mock catalogs for each HOD model using a high resolution $N$-body simulation and adopt a network statistics approach, which provides excellent topological diagnostics for galaxy point distributions. We find three topological anomalies from the spatial distribution of observed LAEs, which are not reproduced by the HOD mocks. We find that Model#2 matches better all network statistics than Model#1, suggesting that the central galaxies in $> 10^{12} h^{-1} M_\odot$ halos at $z \approx 2.67$ need to be less dusty to be bright as LAEs, potentially implying some replenishing channels of pristine gas such as the cold mode accretion.

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 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 use a semi-analytic model for globular cluster (GC) formation built on dark matter merger trees to explore the relative role of formation physics and hierarchical assembly in determining the properties of GC populations. Many previous works have argued that the observed linear relation between total GC mass and halo mass points to a fundamental GC -- dark matter connection or indicates that GCs formed at very high redshift before feedback processes introduced nonlinearity in the baryon-to-dark matter mass relation. We demonstrate that at $M_{\rm vir}(z=0) \gtrsim 10^{11.5} M_{\odot}$, a constant ratio between halo mass and total GC mass is in fact an almost inevitable consequence of hierarchical assembly: by the central limit theorem, it is expected at $z=0$ independent of the GC-to-halo mass relation at the time of GC formation. The GC-to-halo mass relation at $M_{\rm vir}(z=0) < 10^{11.5} M_{\odot}$ is more sensitive to the details of the GC formation process. In our fiducial model, GC formation occurs in galaxies when the gas surface density exceeds a critical value. This model naturally predicts bimodal GC color distributions similar to those observed in nearby galaxies and reproduces the observed relation between GC system metallicity and halo mass. It predicts that the cosmic GC formation rate peaked at $z$ $\sim$ 4, too late for GCs to contribute significantly to the UV luminosity density during reionization.

Young massive star clusters spanning $\sim 10^4 - 10^8 M_\odot$ in mass have been observed to have similar surface brightness profiles. Recent hydrodynamical simulations of star cluster formation have also produced star clusters with this structure. We argue analytically that this type of mass distribution arises naturally in the relaxation from a hierarchically-clustered distribution of stars into a monolithic star cluster through hierarchical merging. We show that arbitrary initial profiles will tend to converge to a universal profile under hierarchical merging, owing to phase-space mixing obeying certain conservation constraints. We perform $N$-body simulations of a pairwise merger of model star clusters and find that mergers readily produce the shallow surface brightness profiles observed in young massive clusters. Finally, we simulate the relaxation of a hierarchically-clustered mass distribution constructed from an idealized fragmentation model. Assuming only power-law spatial and kinematic scaling relations, these numerical experiments are able to reproduce the surface density profiles of observed young massive star clusters. Thus we provide physical motivation for the structure of young massive clusters within the paradigm of hierarchical star formation. This has important implications for the structure of nascent globular clusters.

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.

A standard prediction of galaxy formation theory is that the ionizing background suppresses galaxy formation in haloes with peak circular velocities smaller than Vpeak ~ 20 km/s, rendering the majority of haloes below this scale completely dark. We use a suite of cosmological zoom simulations of Milky Way-like haloes that include central Milky Way disk galaxy potentials to investigate the relationship between subhaloes and ultrafaint galaxies. We find that there are far too few subhaloes within 50 kpc of the Milky Way that had Vpeak > 20 km/s to account for the number of ultrafaint galaxies already known within that volume today. In order to match the observed count, we must populate subhaloes down to Vpeak ~ 6 km/s with ultrafaint dwarfs. The required haloes have peak virial temperatures as low as 1,500 K, well below the atomic hydrogen cooling limit of 10^4 K. Allowing for the possibility that the Large Magellanic Cloud contributes several of the satellites within 50 kpc could potentially raise this threshold to 10 km/s (4,000 K), still below the atomic cooling limit and far below the nominal reionization threshold.

Building upon results of cosmological simulations of ultra-light scalar field dark matter (SFDM), we present a comprehensive model for the density profiles of SFDM haloes as a function of halo virial mass $M_{\rm h}$ and scalar field mass $m$. The central regions of SFDM haloes are dominated by solitons with characteristic densities that increase with increasing halo mass and asymptote to CDM-like profiles at large radii. For scalar field masses $m \sim 10^{-22}$ eV, consistent with large-scale structure observations, $M_{\rm h} \sim 10^{10} \,M_\odot$ haloes have lower core densities than their Cold Dark Matter (CDM) counterparts and this alleviates the Too Big to Fail problem (TBTF) in a regime where feedback is less effective. However, higher-mass SFDM haloes with $M_{\rm h} \sim 10^{11} \,M_\odot$ are denser than their CDM counterparts at small, observationally relevant radii. We use rotation curves of $V \sim 100$ km s$^{-1}$ galaxies from the SPARC database to show that SFDM exacerbates the cusp/core and central density problems seen in CDM at this scale. We conclude that if the conventional cosmological SFDM scaling relations are correct, then baryonic feedback is required to lower densities in SFDM haloes even more so than in CDM. This motivates cosmological and self-consistent hydrodynamic simulations of SFDM to determine whether central soliton structure can be altered by realistic feedback implementations.

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

Metal-poor globular clusters (GCs) are both numerous and ancient, which indicates that they may be important contributors to ionizing radiation in the reionization era. Starting from the observed number density and stellar mass function of old GCs at $z=0$, I compute the contribution of GCs to ultraviolet luminosity functions (UVLFs) in the high-redshift Universe ($10 \gtrsim z \gtrsim 4$). Even under absolutely minimal assumptions - no disruption of GCs and no reduction in GC stellar mass from early times to the present - GC star formation contributes non-negligibly to the UVLF at luminosities that are accessible to the Hubble Space Telescope (HST; $M_{1500} \approx -17$). If the stellar masses of GCs were significantly higher in the past, as is predicted by most models explaining GC chemical anomalies, then GCs dominate the UV emission from many galaxies in existing deep-field observations. On the other hand, it is difficult to reconcile observed UVLFS with models requiring stellar masses at birth that exceed present-day stellar masses by more than a factor of 10. The James Webb Space Telescope will be able to directly detect individual GCs at $z \sim 6$ in essentially all bright galaxies, and many galaxies below the knee of the UVLF, for most of the scenarios considered here. The properties of a subset of high-$z$ galaxies with $-19 \lesssim M_{1500} \lesssim -14$ in HST lensing fields indicate that they may actually be GCs in formation.

In the local Universe, there is a strong division in the star-forming properties of low-mass galaxies, with star formation largely ubiquitous amongst the field population while satellite systems are predominantly quenched. This dichotomy implies that environmental processes play the dominant role in suppressing star formation within this low-mass regime (${M}_{\star} \sim 10^{5.5-8}~{\rm M}_{\odot}$). As shown by observations of the Local Volume, however, there is a non-negligible population of passive systems in the field, which challenges our understanding of quenching at low masses. By applying the satellite quenching models of Fillingham et al. (2015) to subhalo populations in the Exploring the Local Volume In Simulations (ELVIS) suite, we investigate the role of environmental processes in quenching star formation within the nearby field. Using model parameters that reproduce the satellite quenched fraction in the Local Group, we predict a quenched fraction -- due solely to environmental effects -- of $\sim 0.52 \pm 0.26$ within $1< R/R_{\rm vir} < 2$ of the Milky Way and M31. This is in good agreement with current observations of the Local Volume and suggests that the majority of the passive field systems observed at these distances are quenched via environmental mechanisms. Beyond $2~R_{\rm vir}$, however, dwarf galaxy quenching becomes difficult to explain through an interaction with either the Milky Way or M31, such that more isolated, field dwarfs may be self-quenched as a result of star-formation feedback.

We present a suite of cosmological zoom-in simulations at z>5 from the Feedback In Realistic Environments project, spanning a halo mass range M_halo~10^8-10^12 M_sun at z=5. We predict the stellar mass-halo mass relation, stellar mass function, and luminosity function in several bands from z=5-12. The median stellar mass-halo mass relation does not evolve strongly at z=5-12. The faint-end slope of the luminosity function steepens with increasing redshift, as inherited from the halo mass function at these redshifts. Below z~6, the stellar mass function and ultraviolet (UV) luminosity function slightly flatten below M_star~10^4.5 M_sun (fainter than M_1500~-12), owing to the fact that star formation in low-mass halos is suppressed by the ionizing background by the end of reionization. Such flattening does not appear at higher redshifts. We provide redshift-dependent fitting functions for the SFR-M_halo, SFR-M_star, and broad-band magnitude-stellar mass relations. We derive the star formation rate density and stellar mass density at z=5-12 and show that the contribution from very faint galaxies becomes more important at z>8. Furthermore, we find that the decline in the z~6 UV luminosity function brighter than M_1500~-20 is largely due to dust attenuation. Approximately 37% (54%) of the UV luminosity from galaxies brighter than M_1500=-13 (-17) is obscured by dust at z~6. Our results broadly agree with current data and can be tested by future observations.

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.

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.

We present the reconstructed evolution of rest-frame ultra-violet (UV) luminosities of the most massive Milky Way dwarf spheroidal satellite galaxy, Fornax, and its five globular clusters (GCs) across redshift, based on analysis of the stellar fossil record and stellar population synthesis modeling. We find that (1) Fornax's (proto-)GCs can generate $10-100$ times more UV flux than the field population, despite comprising $<\sim 5\%$ of the stellar mass at the relevant redshifts; (2) due to their respective surface brightnesses, it is more likely that faint, compact sources in the Hubble Frontier Fields (HFFs) are GCs hosted by faint galaxies, than faint galaxies themselves. This may significantly complicate the construction of a galaxy UV luminosity function at $z>3$. (3) GC formation can introduce order-of-magnitude errors in abundance matching. We also find that some compact HFF objects are consistent with the reconstructed properties of Fornax's GCs at the same redshifts (e.g., surface brightness, star formation rate), suggesting we may already have detected proto-GCs in the early Universe. Finally, we discuss the prospects for improving the connections between local GCs and proto-GCs detected in the early Universe.

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.

Using high-resolution simulations from the FIRE-2 (Feedback In Realistic Environments) project, we study the effects of discreteness in stellar feedback processes on the evolution of galaxies and the properties of the interstellar medium (ISM). We specifically consider the discretization of supernovae (SNe), including hypernovae (HNe), and sampling the initial mass function (IMF). We study these processes in cosmological simulations of dwarf galaxies with $z=0$ stellar masses $M_{\ast}\sim 10^{4}-3\times10^{6}\,M_\odot$ (halo masses $\sim 10^{9}-10^{10}\,M_\odot$). We show that the discrete nature of individual SNe (as opposed to a model in which their energy/momentum deposition is continuous over time, similar to stellar winds) is crucial in generating a reasonable ISM structure and galactic winds and in regulating dwarf stellar masses. However, once SNe are discretized, accounting for the effects of IMF sampling on continuous mechanisms such as radiative feedback and stellar mass-loss (as opposed to adopting IMF-averaged rates) has weak effects on galaxy-scale properties. We also consider the effects of rare HNe events with energies $\sim 10^{53}\,{\rm erg}$. The effects of HNe are similar to the effects of clustered explosions of SNe -- which are already captured in our default simulation setup -- and do not quench star formation (provided that the HNe do not dominate the total SNe energy budget), which suggests that HNe yield products should be observable in ultra-faint dwarfs today.

A fundamental prediction of the cold dark matter cosmology is the existence of a large number of dark subhalos around galaxies, most of which should be entirely devoid of stars. Confirming the existence of dark substructures stands among the most important empirical challenges in modern cosmology: if they are found and quantified with the mass spectrum expected, then this would close the door on a vast array of competing theories. But in order for observational programs of this kind to reach fruition, we need robust predictions. Here we explore substructure predictions for lensing using galaxy lens-like hosts at z=0.2 from the Illustris simulations both in full hydrodynamics and dark matter only. We quantify substructures more massive than ~ 10^9 M_sun, comparable to current lensing detections derived from HST, Keck, and ALMA. The addition of full hydrodynamics reduces the overall subhalo mass function by about a factor of two. Even for the dark matter only runs, most (~ 85%) lines of sight through projected cylinders of size close to an Einstein radius contain no substructures larger than 10^9 M_sun. The fraction of empty sight lines rises to ~ 95% in full physics simulations. This suggests we will likely need hundreds of strong lensing systems suitable for substructure studies, as well as predictions that include the effects of baryon physics on substructure, to properly constrain cosmological models. Fortunately, the field is poised to fulfill these requirements.

We present a census of variable stars in six M31 dwarf spheroidal satellites observed with the Hubble Space Telescope. We detect 870 RR Lyrae (RRL) stars in the fields of And I (296), II (251), III (111), XV (117), XVI (8), XXVIII (87). We also detect a total of 15 Anomalous Cepheids, three Eclipsing Binaries, and seven field RRL stars compatible with being members of the M31 halo or the Giant Stellar Stream. We derive robust and homogeneous distances to the six galaxies using different methods based on the properties of the RRL stars. Working with the up-to-date set of Period-Wesenheit ($I$, $B$ - $I$) relations published by Marconi et al., we obtain distance moduli of $\mu_0$ = [24.49, 24.16, 24.36, 24.42, 23.70, 24.43] mag (respectively), with systematic uncertainties of 0.08 mag and statistical uncertainties $<$ 0.11 mag. We have considered an enlarged sample of sixteen M31 satellites with published variability studies, and compared their pulsational observables (e.g., periods, amplitudes), with those of fifteen Milky Way satellites for which similar data are available. The properties of the (strictly old) RRL in both satellite systems do not show any significant difference. In particular, we found a strikingly similar correlation between the mean period distribution of the fundamental RRL pulsators (RRab) and the mean metallicities of the galaxies. This indicates that the old RRL progenitors were similar at the early stage in the two environments, suggesting very similar characteristics for the earliest stages of evolution of both satellite systems.

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.

I present a simple phenomenological model for the observed linear scaling of the stellar mass in old globular clusters (GCs) with $z=0$ halo mass in which the stellar mass in GCs scales linearly with progenitor halo mass at $z=6$ above a minimum halo mass for GC formation. This model reproduces the observed $M_{\rm GCs}-M_{\rm halo}$ relation at $z=0$ and results in a prediction for the minimum halo mass at $z=6$ required for hosting one GC: $M_{\rm min}(z=6)=1.07 \times 10^9\,M_{\odot}$. Translated to $z=0$, the mean threshold mass is $M_{\rm halo}(z=0) \approx 2\times 10^{10}\,M_{\odot}$. I explore the observability of GCs in the reionization era and their contribution to cosmic reionization, both of which depend sensitively on the (unknown) ratio of GC birth mass to present-day stellar mass, $\xi$. Based on current detections of $z \gtrsim 6$ objects with $M_{1500} < -17$, values of $\xi > 10$ are strongly disfavored; this, in turn, has potentially important implications for GC formation scenarios. Even for low values of $\xi$, some observed high-$z$ galaxies may actually be GCs, complicating estimates of reionization-era galaxy ultraviolet luminosity functions and constraints on dark matter models. GCs are likely important reionization sources if $5 \lesssim \xi \lesssim 10$. I also explore predictions for the fraction of accreted versus in situ GCs in the local Universe and for descendants of systems at the halo mass threshold of GC formation (dwarf galaxies). An appealing feature of the model presented here is the ability to make predictions for GC properties based solely on dark matter halo merger trees.

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 present a theoretical analysis of some unexplored aspects of relaxed Bose-Einstein condensate dark matter (BECDM) haloes. This type of ultralight bosonic scalar field dark matter is a viable alternative to the standard cold dark matter (CDM) paradigm, as it makes the same large-scale predictions as CDM and potentially overcomes CDM's small-scale problems via a galaxy-scale de Broglie wavelength. We simulate BECDM halo formation through mergers, evolved under the Schr\"odinger-Poisson equations. The formed haloes consist of a soliton core supported against gravitational collapse by the quantum pressure tensor and an asymptotic $r^{-3}$ NFW-like profile. We find a fundamental relation of the core=to-halo mass with the dimensionless invariant $\Xi \equiv \lvert E \rvert/M^3/(Gm/\hbar)^2$ or $M_{\rm c}/M \simeq 2.6 \Xi^{1/3}$, linking the soliton to global halo properties. For $r \geq 3.5 \,r_{\rm c}$ core radii, we find equipartition between potential, classical kinetic, and quantum gradient energies. The haloes also exhibit a conspicuous turbulent behavior driven by the continuous reconnection of vortex lines due to wave interference. We analyse the turbulence 1D velocity power spectrum and find a $k^{-1.1}$ power-law. This suggests the vorticity in BECDM haloes is homogeneous, similar to thermally-driven counterflow BEC systems from condensed matter physics, in contrast to a $k^{-5/3}$ Kolmogorov power-law seen in mechanically-driven quantum systems. The mode where the power spectrum peaks is approximately the soliton width, implying the soliton-sized granules carry most of the turbulent energy in BECDM haloes.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

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

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

We present a proper motion measurement for the halo globular cluster Pyxis, using HST/ACS data as the first epoch, and GeMS/GSAOI Adaptive Optics data as the second, separated by a baseline of about 5 years. This is both the first measurement of the proper motion of Pyxis and the first calibration and use of Multi-Conjugate Adaptive Optics data to measure an absolute proper motion for a faint, distant halo object. Consequently, we present our analysis of the Adaptive Optics data in detail. We obtain a proper motion of mu_alpha cos(delta)=1.09+/-0.31 mas/yr and mu_delta=0.68+/-0.29 mas/yr. From the proper motion and the line-of-sight velocity we find the orbit of Pyxis is rather eccentric with its apocenter at more than 100 kpc and its pericenter at about 30 kpc. We also investigate two literature-proposed associations for Pyxis with the recently discovered ATLAS stream and the Magellanic system. Combining our measurements with dynamical modeling and cosmological numerical simulations we find it unlikely Pyxis is associated with either system. We examine other Milky Way satellites for possible association using the orbit, eccentricity, metallicity, and age as constraints and find no likely matches in satellites down to the mass of Leo II. We propose that Pyxis probably originated in an unknown galaxy, which today is fully disrupted. Assuming that Pyxis is bound and not on a first approach, we derive a 68% lower limit on the mass of the Milky Way of 0.95*10^12 M_sun.

We present the rest-1500\AA\ UV luminosity functions (LF) for star-forming galaxies during the cosmic \textit{high noon} -- the peak of cosmic star formation rate at $1.5<z<3$. We use deep NUV imaging data obtained as part of the \textit{Hubble} Ultra-Violet Ultra Deep Field (UVUDF) program, along with existing deep optical and NIR coverage on the HUDF. We select F225W, F275W and F336W dropout samples using the Lyman break technique, along with samples in the corresponding redshift ranges selected using photometric redshifts and measure the rest-frame UV LF at $z\sim1.7,2.2,3.0$ respectively, using the modified maximum likelihood estimator. We perform simulations to quantify the survey and sample incompleteness for the UVUDF samples to correct the effective volume calculations for the LF. We select galaxies down to $M_{UV}=-15.9,-16.3,-16.8$ and fit a faint-end slope of $\alpha=-1.20^{+0.10}_{-0.13}, -1.32^{+0.10}_{-0.14}, -1.39^{+0.08}_{-0.12}$ at $1.4<z<1.9$, $1.8<z<2.6$, and $2.4<z<3.6$, respectively. We compare the star formation properties of $z\sim2$ galaxies from these UV observations with results from H\alpha\ and UV$+$IR observations. We find a lack of high SFR sources in the UV LF compared to the H\alpha\ and UV$+$IR, likely due to dusty SFGs not being properly accounted for by the generic $IRX-\beta$ relation used to correct for dust. We compute a volume-averaged UV-to-H\alpha\ ratio by \textit{abundance matching} the rest-frame UV LF and H\alpha\ LF. We find an increasing UV-to-H\alpha\ ratio towards low mass galaxies ($M_\star \lesssim 5\times10^9$ M$_\odot$). We conclude that this could be due to a larger contribution from starbursting galaxies compared to the high-mass end.

Motivated by the stellar fossil record of Local Group (LG) dwarf galaxies, we show that the star-forming ancestors of the faintest ultra-faint dwarf galaxies (UFDs; ${\rm M}_{\rm V}$ $\sim -2$ or ${\rm M}_{\star}$ $\sim 10^{2}$ at $z=0$) had ultra-violet (UV) luminosities of ${\rm M}_{\rm UV}$ $\sim -3$ to $-6$ during reionization ($z\sim6-10$). The existence of such faint galaxies has substantial implications for early epochs of galaxy formation and reionization. If the faint-end slopes of the UV luminosity functions (UVLFs) during reionization are steep ($\alpha\lesssim-2$) to ${\rm M}_{\rm UV}$ $\sim -3$, then: (i) the ancestors of UFDs produced $>50$% of UV flux from galaxies; (ii) galaxies can maintain reionization with escape fractions that are $>$2 times lower than currently-adopted values; (iii) direct HST and JWST observations may detect only $\sim10-50$% of the UV light from galaxies; (iv) the cosmic star formation history increases by $\gtrsim4-6$ at $z\gtrsim6$. Significant flux from UFDs, and resultant tensions with LG dwarf galaxy counts, are reduced if the high-redshift UVLF turns over. Independent of the UVLF shape, the existence of a large population of UFDs requires a non-zero luminosity function to ${\rm M}_{\rm UV}$ $\sim -3$ during reionization.

We confirm that the object DDO216-A1 is a substantial globular cluster at the center of Local Group galaxy DDO216 (the Pegasus dwarf irregular), using Hubble Space Telescope ACS imaging. By fitting isochrones, we find the cluster metallicity to be [M/H] = -1.6 +/-0.2, for reddening E(B-V) = 0.16 +/-0.02; the best-fit age is 12.3 +/-0.8 Gyr. There are ~30 RR Lyrae variables in the cluster; the magnitude of the fundamental mode pulsators gives a distance modulus of 24.77 +/-0.08 - identical to the host galaxy. The ratio of overtone to fundamental mode variables and their mean periods make DDO216-A1 an Oosterhoff Type I cluster. We find an I-band central surface brightness 20.85 +/-0.17 F814W mag per square arcsecond, a half-light radius of 3.1 arcsec (13.4 pc), and an absolute magnitude M814 = -7.90 +/-0.16 (approximately 10^5 solar masses). King models fit to the cluster give the core radius and concentration index, r_c = 2.1" +/-0.9" and c = 1.24 +/-0.39. The cluster is an "extended" cluster somewhat typical of some dwarf galaxies and the outer halo of the Milky Way. The cluster is projected <30 pc south of the center of DDO216, unusually central compared to most dwarf galaxy globular clusters. Analytical models of dynamical friction and tidal destruction suggest that it probably formed at a larger distance, up to ~1 kpc, and migrated inward. DDO216 has an unexceptional cluster specific frequency, S_N = 10. DDO216 is the lowest-luminosity Local Group galaxy to host a 10^5 solar mass globular cluster, and the only transition-type (dSph/dIrr) in the Local Group with a globular.

We use Local Group galaxy counts together with the ELVIS N-body simulations to explore the relationship between the scatter and slope in the stellar mass vs. halo mass relation at low masses, $M_\star \simeq 10^5 - 10^8 M_\odot$. Assuming models with log-normal scatter about a median relation of the form $M_\star \propto M_\mathrm{halo}^\alpha$, the preferred log-slope steepens from $\alpha \simeq 1.8$ in the limit of zero scatter to $\alpha \simeq 2.6$ in the case of $2$ dex of scatter in $M_\star$ at fixed halo mass. We provide fitting functions for the best-fit relations as a function of scatter, including cases where the relation becomes increasingly stochastic with decreasing mass. We show that if the scatter at fixed halo mass is large enough ($\gtrsim 1$ dex) and if the median relation is steep enough ($\alpha \gtrsim 2$), then the "too-big-to-fail" problem seen in the Local Group can be self-consistently eliminated in about $\sim 5-10\%$ of realizations. This scenario requires that the most massive subhalos host unobservable ultra-faint dwarfs fairly often; we discuss potentially observable signatures of these systems. Finally, we compare our derived constraints to recent high-resolution simulations of dwarf galaxy formation in the literature. Though simulation-to-simulation scatter in $M_\star$ at fixed $M_\mathrm{halo}$ is large among separate authors ($\sim 2$ dex), individual codes produce relations with much less scatter and usually give relations that would over-produce local galaxy counts.

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

We test the DARKexp model for relaxed, self-gravitating, collisionless systems against equilibrium dark matter halos from the Millennium-II simulation. While limited tests of DARKexp against simulations and observations have been carried out elsewhere, this is the first time the testing is done with a large sample of simulated halos spanning a factor of ~ 50 in mass, and using independent fits to density and energy distributions. We show that DARKexp, a one shape parameter family, provides very good fits to the shapes of density profiles, \rho(r), and differential energy distributions, N(E), of individual simulated halos. The best fit shape parameter $\phi_{0}$ obtained from the two types of fits are correlated, though with scatter. Our most important conclusions come from \rho(r) and N(E) that have been averaged over many halos. These show that the bulk of the deviations between DARKexp and individual Millennium-II halos come from halo-to-halo fluctuations, likely driven by substructure, and other density perturbations. The average \rho(r) and N(E) are quite smooth and follow DARKexp very closely. The only deviation that remains after averaging is small, and located at most bound energies for N(E) and smallest radii for \rho(r). Since the deviation is confined to 3-4 smoothing lengths, and is larger for low mass halos, it is likely due to numerical resolution effects.

Recent studies of galaxies in the local Universe, including those in the Local Group, find that the efficiency of environmental (or satellite) quenching increases dramatically at satellite stellar masses below ~ $10^8\ {\rm M}_{\odot}$. This suggests a physical scale where quenching transitions from a slow "starvation" mode to a rapid "stripping" mode at low masses. We investigate the plausibility of this scenario using observed HI surface density profiles for a sample of 66 nearby galaxies as inputs to analytic calculations of ram-pressure and viscous stripping. Across a broad range of host properties, we find that stripping becomes increasingly effective at $M_{*} < 10^{8-9}\ {\rm M}_{\odot}$, reproducing the critical mass scale observed. However, for canonical values of the circumgalactic medium density ($n_{\rm halo} < 10^{-3.5}$ ${\rm cm}^{-3}$), we find that stripping is not fully effective; infalling satellites are, on average, stripped of < 40 - 70% of their cold gas reservoir, which is insufficient to match observations. By including a host halo gas distribution that is clumpy and therefore contains regions of higher density, we are able to reproduce the observed HI gas fractions (and thus the high quenched fraction and short quenching timescale) of Local Group satellites, suggesting that a host halo with clumpy gas may be crucial for quenching low-mass systems in Local Group-like (and more massive) host halos.

We use particle data from the Illustris simulation, combined with individual kinematic constraints on the mass of the Milky Way (MW) at specific distances from the Galactic center, to infer the radial distribution of the MW's dark matter halo mass. Our method allows us to convert any constraint on the mass of the MW within a fixed distance to a full circular velocity profile to the MW's virial radius. As primary examples, we take two recent (and discrepant) measurements of the total mass within 50 kpc of the Galaxy and find they imply very different mass profiles and stellar masses for the Galaxy. The dark-matter-only version of the Illustris simulation enables us to compute the effects of galaxy formation on such constraints on a halo-by-halo basis; on small scales, galaxy formation enhances the density relative to dark-matter-only runs, while the total mass density is approximately 20% lower at large Galactocentric distances. We are also able to quantify how current and future constraints on the mass of the MW within specific radii will be reflected in uncertainties on its virial mass: even a measurement of M(<50 kpc) with essentially perfect precision still results in a 20% uncertainty on the virial mass of the Galaxy, while a future measurement of M(<100 kpc) with 10% errors would result in the same level of uncertainty. We expect that our technique will become even more useful as (1) better kinematic constraints become available at larger distances and (2) cosmological simulations provide even more faithful representations of the observable Universe.

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.

Near-field cosmology -- using detailed observations of the Local Group and its environs to study wide-ranging questions in galaxy formation and dark matter physics -- has become a mature and rich field over the past decade. There are lingering concerns, however, that the relatively small size of the present-day Local Group ($\sim 2$ Mpc diameter) imposes insurmountable sample-variance uncertainties, limiting its broader utility. We consider the region spanned by the Local Group's progenitors at earlier times and show that it reaches $3' \approx 7$ co-moving Mpc in linear size (a volume of $\approx 350\,{\rm Mpc}^3$) at $z=7$. This size at early cosmic epochs is large enough to be representative in terms of the matter density and counts of dark matter halos with $M_{\rm vir}(z=7) \lesssim 2\times 10^{9}\,M_{\odot}$. The Local Group's stellar fossil record traces the cosmic evolution of galaxies with $10^{3} \lesssim M_{\star}(z=0) / M_{\odot} \lesssim 10^{9}$ (reaching $M_{1500} > -9$ at $z\sim7$) over a region that is comparable to or larger than the Hubble Ultra-Deep Field (HUDF) for the entire history of the Universe. It is highly complementary to the HUDF, as it probes much fainter galaxies but does not contain the intrinsically rarer, brighter sources that are detectable in the HUDF. Archaeological studies in the Local Group also provide the ability to trace the evolution of individual galaxies across time as opposed to evaluating statistical connections between temporally distinct populations. In the JWST era, resolved stellar populations will probe regions larger than the HUDF and any deep JWST fields, further enhancing the value of near-field cosmology.

The Initial Star formation and Lifetimes of Andromeda Satellites (ISLAndS) project uses Hubble Space Telescope imaging to study a representative sample of six Andromeda dSph satellite companion galaxies. The main goal of the program is to determine whether the star formation histories (SFHs) of the Andromeda dSph satellites demonstrate significant statistical differences from those of the Milky Way, which may be attributable to the different properties of their local environments. Our observations reach the oldest main sequence turn-offs, allowing a time resolution at the oldest ages of ~ 1 Gyr, which is comparable to the best achievable resolution in the MW satellites. We find that the six dSphs present a variety of SFHs that are not strictly correlated with luminosity or present distance from M31. Specifically, we find a significant range in quenching times (lookback times from 9 to 6 Gyr), but with all quenching times more than ~ 6 Gyr ago. In agreement with observations of Milky Way companions of similar mass, there is no evidence of complete quenching of star formation by the cosmic UV background responsible for reionization, but the possibility of a degree of quenching at reionization cannot be ruled out. We do not find significant differences between the SFHs of the three members of the vast, thin plane of satellites and the three off-plane dSphs. The primary difference between the SFHs of the ISLAndS dSphs and Milky Way dSph companions of similar luminosities and host distances is the absence of very late quenching (< 5 Gyr ago) dSphs in the ISLAndS sample. Thus, models that can reproduce satellite populations with and without late quenching satellites will be of extreme interest.

Based on data aquired in 13 orbits of HST time, we present a detailed evolutionary history of the M31 dSph satellite Andromeda XVI, including its life-time star formation history, the spatial distribution of its stellar populations, and the properties of its variable stars. And XVI is characterized by prolonged star formation activity from the oldest epochs until star formation was quenched ~6 Gyr ago, and, notably, only half of the mass in stars of And XVI was in place 10 Gyr ago. And XVI appears to be a low mass galaxy for which the early quenching by either reionization or starburst feedback seems highly unlikely, and thus, is most likely due to an environmental effect (e.g., an interaction), possibly connected to a late infall in the densest regions of the Local Group. Studying the star formation history as a function of galactocentric radius, we detect a mild gradient in the star formation history: the star formation activity between 6 and 8 Gyr ago is significantly stronger in the central regions than in the external regions, although the quenching age appears to be the same, within 1 Gyr. We also report the discovery of 9 RR Lyrae stars, 8 of which belong to And XVI. The RR Lyrae stars allow a new estimate of the distance, (m-M)0= 23.72+/-0.09 mag, which is marginally larger than previous estimates based on the tip of the red giant branch.

We present in this study a first analysis of the astrometric error budget of absolute astrometry relative to background galaxies using adaptive optics. We use for this analysis multi-conjugated adaptive optics (MCAO) images obtained with GeMS/GSAOI at Gemini South. We find that it is possible to obtain 0.3 mas reference precision in a random field with 1 hour on source using faint background galaxies. Systematic errors are correctable below that level, such that the overall error is approximately 0.4 mas. Because the reference sources are extended, we find it necessary to correct for the dependency of the PSF centroid on the used aperture size, which would otherwise cause an important bias. This effect needs also to be considered for Extremely Large Telescopes (ELTs). When this effect is corrected, ELTs have the potential to measure proper motions of dwarfs galaxies around M31 with 10 km/s accuracy over a baseline of 5 years.

Sterile neutrinos comprise an entire class of dark matter models that, depending on their production mechanism, can be hot, warm, or cold dark matter. We simulate the Local Group and representative volumes of the Universe in a variety of sterile neutrino models, all of which are consistent with the possible existence of a radiative decay line at ~3.5 keV. We compare models of production via resonances in the presence of a lepton asymmetry (suggested by Shi & Fuller 1999) to "thermal" models. We find that properties in the highly nonlinear regime - e.g., counts of satellites and internal properties of halos and subhalos - are insensitive to the precise fall-off in power with wavenumber, indicating that nonlinear evolution essentially washes away differences in the initial (linear) matter power spectrum. In the quasi-linear regime at higher redshifts, however, quantitative differences in the 3D matter power spectra remain, raising the possibility that such models can be tested with future observations of the Lyman-alpha forest. While many of the sterile neutrino models largely eliminate multiple small-scale issues within the Cold Dark Matter (CDM) paradigm, we show that these models may be ruled out in the near future via discoveries of additional dwarf satellites in the Local Group.

The anomalous 3.55 keV X-ray line recently detected towards a number of massive dark matter objects may be interpreted as the radiative decays of 7.1 keV mass sterile neutrino dark matter. Depending on its parameters, the sterile neutrino can range from cold to warm dark matter with small-scale suppression that differs in form from commonly-adopted thermal warm dark matter. Here, we numerically investigate the subhalo properties for 7.1 keV sterile neutrino dark matter produced via the resonant Shi-Fuller mechanism. Using accurate matter power spectra, we run cosmological zoom-in simulations of a Milky Way-sized halo and explore the abundance of massive subhalos, their radial distributions, and their internal structure. We also simulate the halo with thermal 2.0 keV warm dark matter for comparison and discuss quantitative differences. We find that the resonantly produced sterile neutrino model for the 3.55 keV line provides a good description of structures in the Local Group, including the number of satellite dwarf galaxies and their radial distribution, and largely mitigates the too-big-to-fail problem. Future searches for satellite galaxies by deep surveys, such as the Dark Energy Survey, Large Synoptic Survey Telescope, and Wide Field Infrared Survey Telescope, will be a strong direct test of warm dark matter scenarios.

We present ultra-high resolution cosmological hydrodynamic simulations of $M_*\simeq10^{4-6.3}M_{\odot}$ dwarf galaxies that form within $M_{v}=10^{9.5-10}M_{\odot}$ dark matter halos. Our simulations rely on the FIRE implementation of star formation feedback and were run with high enough force and mass resolution to directly resolve stellar and dark matter structure on the ~200 pc scales of interest for classical and ultra-faint dwarfs in the Local Group. The resultant galaxies sit on the $M_*$ vs. $M_{v}$ relation required to match the Local Group stellar mass function. They have bursty star formation histories and also form with half-light radii and metallicities that broadly match those observed for local dwarfs at the same stellar mass. We demonstrate that it is possible to create a large (~1 kpc) dark matter core in a cosmological simulation of an $M_*\simeq10^{6.5}M_{\odot}$ dwarf galaxy that resides within an $M_{v}=10^{10}M_{\odot}$ halo -- precisely the scale of interest for resolving the Too Big to Fail problem. However, these large cores are not ubiquitous and appear to correlate closely with the star formation histories of the dwarfs: dark matter cores are largest in systems that form their stars late ($z\lesssim2$), after the early epoch of cusp building mergers has ended. Our $M_*\simeq10^4M_{\odot}$ dwarf retains a cuspy dark matter halo density profile that matches almost identically that of a dark-matter only run of the same system. Despite forming in a field environment, this very low mass dwarf has observable properties that match closely to those of ultra-faint satellite galaxies of the Milky Way, including a uniformly old stellar population (>10 Gyr). Though ancient, most of the stars in our ultra-faint form after reionization; the UV field acts mainly to suppress fresh gas accretion, not to boil away gas that is already present in the proto-dwarf.

We constrain the evolution of the galaxy stellar mass function from 2 < z < 5 for galaxies with stellar masses as low as 10^5 Msun by combining star formation histories of Milky Way satellite galaxies derived from deep Hubble Space Telescope observations with merger trees from the ELVIS suite of N-body simulations. This approach extends our understanding more than two orders of magnitude lower in stellar mass than is currently possible by direct imaging. We find the faint end slopes of the mass functions to be alpha= -1.42(+0.07/-0.05) at z = 2 and alpha = -1.57^(+0.06/-0.06) at z = 5, and show the slope only weakly evolves from z = 5 to z = 0. Our findings are in stark contrast to a number of direct detection studies that suggest slopes as steep as alpha = -1.9 at these epochs. Such a steep slope would result in an order of magnitude too many luminous Milky Way satellites in a mass regime that is observationally complete (Mstar > 2*10^5 Msun at z = 0). The most recent studies from ZFOURGE and CANDELS also suggest flatter faint end slopes that are consistent with our results, but with a lower degree of precision. This work illustrates the strong connections between low and high-z observations when viewed through the lens of LCDM numerical simulations.

We infer the UV luminosities of Local Group galaxies at early cosmic times ($z \sim 2$ and $z \sim 7$) by combining stellar population synthesis modeling with star formation histories derived from deep color-magnitude diagrams constructed from Hubble Space Telescope (HST) observations. Our analysis provides a basis for understanding high-$z$ galaxies - including those that may be unobservable even with the James Webb Space Telescope (JWST) - in the context of familiar, well-studied objects in the very low-$z$ Universe. We find that, at the epoch of reionization, all Local Group dwarfs were less luminous than the faintest galaxies detectable in deep HST observations of blank fields. We predict that JWST will observe $z \sim 7$ progenitors of galaxies similar to the Large Magellanic Cloud today; however, the HST Frontier Fields initiative may already be observing such galaxies, highlighting the power of gravitational lensing. Consensus reionization models require an extrapolation of the observed blank-field luminosity function at $z \approx 7$ by at least two orders of magnitude in order to maintain reionization. This scenario requires the progenitors of the Fornax and Sagittarius dwarf spheroidal galaxies to be contributors to the ionizing background at $z \sim 7$. Combined with numerical simulations, our results argue for a break in the UV luminosity function from a faint-end slope of $\alpha \sim -2$ at $M_{\rm UV} < -13$ to $\alpha \sim -1.2$ at lower luminosities. Applied to photometric samples at lower redshifts, our analysis suggests that HST observations in lensing fields at $z \sim 2$ are capable of probing galaxies with luminosities comparable to the expected progenitor of Fornax.

The vast majority of dwarf satellites orbiting the Milky Way and M31 are quenched, while comparable galaxies in the field are gas-rich and star-forming. Assuming that this dichotomy is driven by environmental quenching, we use the ELVIS suite of N-body simulations to constrain the characteristic timescale upon which satellites must quench following infall into the virial volumes of their hosts. The high satellite quenched fraction observed in the Local Group demands an extremely short quenching timescale (~ 2 Gyr) for dwarf satellites in the mass range Mstar ~ 10^6-10^8 Msun. This quenching timescale is significantly shorter than that required to explain the quenched fraction of more massive satellites (~ 8 Gyr), both in the Local Group and in more massive host halos, suggesting a dramatic change in the dominant satellite quenching mechanism at Mstar < 10^8 Msun. Combining our work with the results of complementary analyses in the literature, we conclude that the suppression of star formation in massive satellites (Mstar ~ 10^8 - 10^11 Msun) is broadly consistent with being driven by starvation, such that the satellite quenching timescale corresponds to the cold gas depletion time. Below a critical stellar mass scale of ~ 10^8 Msun, however, the required quenching times are much shorter than the expected cold gas depletion times. Instead, quenching must act on a timescale comparable to the dynamical time of the host halo. We posit that ram-pressure stripping can naturally explain this behavior, with the critical mass (of Mstar ~ 10^8 Msun) corresponding to halos with gravitational restoring forces that are too weak to overcome the drag force encountered when moving through an extended, hot circumgalactic medium.

We compare the dynamics of satellite galaxies in the Sloan Digital Sky Survey to simple models in order to test the hypothesis that a large fraction of satellites co-rotate in coherent planes. We confirm the previously-reported excess of co-rotating satellite pairs located near diametric opposition with respect to the host, but show that this signal is unlikely to be due to rotating discs (or planes) of satellites. In particular, no overabundance of co-rotating satellites pairs is observed within $\sim 20^{\circ}-50^{\circ}$ of direct opposition, as would be expected for planar distributions inclined relative to the line-of-sight. Instead, the excess co-rotation for satellite pairs within $\sim 10^{\circ}$ of opposition is consistent with random noise associated with undersampling of an underlying isotropic velocity distribution. We conclude that at most $10\%$ of the hosts in our sample harbor co-rotating satellite planes (as traced by the luminous satellite population).

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

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.

Using the Sloan Digital Sky Survey, we examine the quenching of satellite galaxies around isolated Milky Way-like hosts in the local Universe. We find that the efficiency of satellite quenching around isolated galaxies is low and roughly constant over two orders of magnitude in satellite stellar mass ($M_{*}$ = $10^{8.5}-10^{10.5} \, M_{\odot}$), with only $\sim~20\%$ of systems quenched as a result of environmental processes. While largely independent of satellite stellar mass, satellite quenching does exhibit clear dependence on the properties of the host. We show that satellites of passive hosts are substantially more likely to be quenched than those of star-forming hosts, and we present evidence that more massive halos quench their satellites more efficiently. These results extend trends seen previously in more massive host halos and for higher satellite masses. Taken together, it appears that galaxies with stellar masses larger than about $10^{8}~M_{\odot}$ are uniformly resistant to environmental quenching, with the relative harshness of the host environment likely serving as the primary driver of satellite quenching. At lower stellar masses ($< 10^{8}~M_{\odot}$), however, observations of the Local Group suggest that the vast majority of satellite galaxies are quenched, potentially pointing towards a characteristic satellite mass scale below which quenching efficiency increases dramatically.

Reionizing the Universe with galaxies appears to require significant star formation in low-mass halos at early times, while local dwarf galaxy counts tell us that star formation has been minimal in small halos around us today. Using simple models and the ELVIS simulation suite, we show that reionization scenarios requiring appreciable star formation in halos with $M_{\rm vir} \approx 10^{8}\,M_{\odot}$ at $z=8$ are in serious tension with galaxy counts in the Local Group. This tension originates from the seemingly inescapable conclusion that 30 - 60 halos with $M_{\rm vir} > 10^{8}\,M_{\odot}$ at $z=8$ will survive to be distinct bound satellites of the Milky Way at $z = 0$. Reionization models requiring star formation in such halos will produce dozens of bound galaxies in the Milky Way's virial volume today (and 100 - 200 throughout the Local Group), each with $\gtrsim 10^{5}\,M_{\odot}$ of old stars ($\gtrsim 13$ Gyr). This exceeds the stellar mass function of classical Milky Way satellites today, even without allowing for the (significant) post-reionization star formation observed in these galaxies. One possible implication of these findings is that star formation became sharply inefficient in halos smaller than $\sim 10^9 \,M_{\odot}$ at early times, implying that the high-$z$ luminosity function must break at magnitudes brighter than is often assumed (at ${\rm M_{UV}} \approx -14$). Our results suggest that JWST (and possibly even HST with the Frontier Fields) may realistically detect the faintest galaxies that drive reionization. It remains to be seen how these results can be reconciled with the most sophisticated simulations of early galaxy formation at present, which predict substantial star formation in $M_{\rm vir} \sim 10^8 \, M_{\odot}$ halos during the epoch of reionization.

We present the first comparison between the lifetime star formation histories (SFHs) of M31 and Milky Way (MW) satellites. Using the Advanced Camera for Surveys aboard the Hubble Space Telescope, we obtained deep optical imaging of Andromeda II (M$_{V} = -$12.0; log(M$_{\star}$/M$_{\odot}$) $\sim$ 6.7) and Andromeda XVI (M$_{V} = -$7.5; log(M$_{\star}$/M$_{\odot}$) $\sim$ 4.9) yielding color-magnitude diagrams that extend at least 1 magnitude below the oldest main sequence turnoff, and are similar in quality to those available for the MW companions. And II and And XVI show strikingly similar SFHs: both formed 50-70% of their total stellar mass between 12.5 and 5 Gyr ago (z$\sim$5-0.5) and both were abruptly quenched $\sim$ 5 Gyr ago (z$\sim$0.5). The predominance of intermediate age populations in And XVI makes it qualitatively different from faint companions of the MW and clearly not a pre-reionization fossil. Neither And II nor And XVI appears to have a clear analog among MW companions, and the degree of similarity in the SFHs of And II and And XVI is not seen among comparably faint-luminous pairs of MW satellites. These findings provide hints that satellite galaxy evolution may vary substantially among hosts of similar stellar mass. Although comparably deep observations of more M31 satellites are needed to further explore this hypothesis, our results underline the need for caution when interpreting satellite galaxies of an individual system in a broader cosmological context.

We study dwarf satellite galaxy quenching using observations from the Geha et al. (2012) NSA/SDSS catalog together with LCDM cosmological simulations to facilitate selection and interpretation. We show that fewer than 30% of dwarfs (M* ~ 10^8.5-10^9.5 Msun) identified as satellites within massive host halos (Mhost ~ 10^12.5-10^14 Msun) are quenched, in spite of the expectation from simulations that half of them should have been accreted more than 6 Gyr ago. We conclude that whatever the action triggering environmental quenching of dwarf satellites, the process must be highly inefficient. We investigate a series of simple, one-parameter quenching models in order to understand what is required to explain the low quenched fraction and conclude that either the quenching timescale is very long (> 9.5 Gyr, a "slow starvation" scenario) or that the environmental trigger is not well matched to accretion within the virial volume. We discuss these results in light of the fact that most of the low mass dwarf satellites in the Local Group are quenched, a seeming contradiction that could point to a characteristic mass scale for satellite quenching.

We compare the dynamical masses of dwarf galaxies in the Local Group (LG) to the predicted masses of halos in the ELVIS suite of $\Lambda$CDM simulations, a sample of 48 Galaxy-size hosts, 24 of which are in paired configuration similar to the LG. We enumerate unaccounted-for dense halos ($V_\mathrm{max} \gtrsim 25$ km s$^{-1}$) in these volumes that at some point in their histories were massive enough to have formed stars in the presence of an ionizing background ($V_\mathrm{peak} > 30$ km s$^{-1}$). Within 300 kpc of the Milky Way, the number of unaccounted-for massive halos ranges from 2 - 25 over our full sample. Moreover, this "too big to fail" count grows as we extend our comparison to the outer regions of the Local Group: within 1.2 Mpc of either giant we find that there are 12-40 unaccounted-for massive halos. This count excludes volumes within 300 kpc of both the MW and M31, and thus should be largely unaffected by any baryonically-induced environmental processes. According to abundance matching -- specifically abundance matching that reproduces the Local Group stellar mass function -- all of these missing massive systems should have been quite bright, with $M_\star > 10^6M_\odot$. Finally, we use the predicted density structure of outer LG dark matter halos together with observed dwarf galaxy masses to derive an $M_\star-V_\mathrm{max}$ relation for LG galaxies that are outside the virial regions of either giant. We find that there is no obvious trend in the relation over three orders of magnitude in stellar mass (a "common mass" relation), from $M_\star \sim 10^8 - 10^5 M_\odot$, in drastic conflict with the tight relation expected for halos that are unaffected by reionization. Solutions to the too big to fail problem that rely on ram pressure stripping, tidal effects, or statistical flukes appear less likely in the face of these results.

We compare predictions of cooled masses and cooling rates from three stripped-down Semi-Analytic Models (SAMs) of galaxy formation with the results of N-body+SPH simulations with gas particle mass of 3.9x10^6 Msun, where radiative cooling of a gas of primordial composition is implemented. We also run a simulation where cooling is switched on at redshift ~2, in order to test cooling models in a regime in which their approximations are expected to be valid. We confirm that cooling models implemented in SAMs are able to predict the amount of cooled mass at z=0 to within ~20 per cent. However, some relevant discrepancies are found. (i) When the contribution from poorly resolved halos is subtracted out, SAMs tend to under-predict by ~30 per cent the mass that cools in the infall-dominated regime. (ii) At large halo masses SAMs tend to over-predict cooling rates, though the numerical result may be affected by the use of SPH. (iii) As found in our previous work, cooling rates are found to be significantly affected by model details: simulations disfavour models with large cores and with quenching of cooling at major mergers. (iv) When cooling is switched on at z~2, cold gas accumulates very quickly in the simulated halos. This accumulation is reproduced by SAMs with varying degrees of accuracy.

We compare the kinematics of M31's satellite galaxies to the mass profiles of the subhaloes they are expected to inhabit in LCDM. We consider the most massive subhaloes of an approximately M31-sized halo, following the assumption of a monotonic galaxy luminosity-to-subhalo mass mapping. While this abundance matching relation is consistent with the kinematic data for galaxies down to the luminosity of the bright satellites of the Milky Way and M31, it is not consistent with kinematic data for fainter dwarf galaxies (those with L <~ 10^8 Lsun). Comparing the kinematics of M31's dwarf Spheroidal (dSph) satellites to the subhaloes reveals that M31's dSph satellites are too low density to be consistent with the subhaloes' mass profiles. A similar discrepancy has been reported between Milky Way dSphs and their predicted subhaloes, the "too big to fail" problem. By contrast, total mass profiles of the dwarf Elliptical (and similarly bright) satellites are consistent with the subhaloes. However, they suffer from large systematic uncertainties in their dark matter content because of substantial (and potentially dominant) contributions from baryons within their half-light radii.

We present an analysis of the star formation history (SFH) of a field near the half light radius in the Local Group dwarf irregular galaxy IC 1613 based on deep Hubble Space Telescope Advanced Camera for Surveys imaging. Our observations reach the oldest main sequence turn-off, allowing a time resolution at the oldest ages of ~1 Gyr. Our analysis shows that the SFH of the observed field in IC 1613 is consistent with being constant over the entire lifetime of the galaxy. These observations rule out an early dominant episode of star formation in IC 1613. We compare the SFH of IC 1613 with expectations from cosmological models. Since most of the mass is in place at early times for low mass halos, a naive expectation is that most of the star formation should have taken place at early times. Models in which star formation follows mass accretion result in too many stars formed early and gas mass fractions which are too low today (the "over-cooling problem"). The depth of the present photometry of IC 1613 shows that, at a resolution of ~1 Gyr, the star formation rate is consistent with being constant, at even the earliest times, which is difficult to achieve in models where star formation follows mass assembly.

We measured velocities of 862 individual red giant stars in seven isolated dwarf galaxies in the Local Group: NGC 6822, IC 1613, VV 124 (UGC 4879), the Pegasus dwarf irregular galaxy (DDO 216), Leo A, Cetus, and Aquarius (DDO 210). We also computed velocity dispersions, taking into account the measurement uncertainties on individual stars. None of the isolated galaxies is denser than the densest Local Group satellite galaxy. Furthermore, the isolated dwarf galaxies have no obvious distinction in the velocity dispersion--half-light radius plane from the satellite galaxies of the Milky Way and M31. The similarity of the isolated and satellite galaxies' dynamics and structural parameters imposes limitations on environmental solutions to the too-big-to-fail problem, wherein there are fewer dense dwarf satellite galaxies than would be expected from cold dark matter simulations. This data set also has many other applications for dwarf galaxy evolution, including the transformation of dwarf irregular into dwarf spheroidal galaxies. We intend to explore these issues in future work.

We introduce a set of high-resolution dissipationless simulations that model the Local Group (LG) in a cosmological context: Exploring the Local Volume in Simulations (ELVIS). The suite contains 48 Galaxy-size halos, each within high-resolution volumes that span 2-5 Mpc in size, and each resolving thousands of systems with masses below the atomic cooling limit. Half of the ELVIS galaxy halos are in paired configurations similar to the Milky Way (MW) and M31; the other half are isolated, mass-matched analogs. We find no difference in the abundance or kinematics of substructure within the virial radii of isolated versus paired hosts. On Mpc scales, however, LG-like pairs average almost twice as many companions and the velocity field is kinematically hotter and more complex. We present a refined abundance matching relation between stellar mass and halo mass that reproduces the observed satellite stellar mass functions of the MW and M31 down to the regime where incompleteness is an issue, $M_\star \sim 5\times 10^5 \, M_\odot$. Within a larger region spanning approximately 3 Mpc, the same relation predicts that there should be $\sim$ 1000 galaxies with $M_\star > 10^{3}\,M_\odot$ awaiting discovery. We show that up to 50% of halos within 1 Mpc of the MW or M31 could be systems that have previously been within the virial radius of either giant. By associating never-accreted halos with gas-rich dwarfs, we show that there are plausibly 50 undiscovered dwarf galaxies with HI masses $> 10^5\,M_\odot$ within the Local Volume. The radial velocity distribution of these predicted gas-rich dwarfs can be used to inform follow-up searches based on ultra-compact high-velocity clouds found in the ALFALFA survey.

We use the Millennium Simulation series to study the relation between the accretion history (MAH) and mass profile of cold dark matter halos. We find that the mean density within the scale radius, r_{-2} (where the halo density profile has isothermal slope), is directly proportional to the critical density of the Universe at the time when the main progenitor's virial mass equals the mass enclosed within r_{-2}. Scaled to these characteristic values of mass and density, the mean MAH, expressed in terms of the critical density of the Universe, M(\rho_{crit}(z)), resembles that of the enclosed density profile, M(<\rho >), at z=0. Both follow closely the NFW profile, suggesting that the similarity of halo mass profiles originates from the mass-independence of halo MAHs. Support for this interpretation is provided by outlier halos whose accretion histories deviate from the NFW shape; their mass profiles show correlated deviations from NFW and are better approximated by Einasto profiles. Fitting both M(<\rho >) and M(\rho_{crit}) with either NFW or Einasto profiles yield concentration and shape parameters that are correlated, confirming and extending earlier work linking the concentration of a halo with its accretion history. These correlations also confirm that halo structure is insensitive to initial conditions: only halos whose accretion histories differ greatly from the NFW shape show noticeable deviations from NFW in their mass profiles. As a result, the NFW profile provides acceptable fits to hot dark matter halos, which do not form hierarchically, and for fluctuation power spectra other than CDM. Our findings, however, predict a subtle but systematic dependence of mass profile shape on accretion history which, if confirmed, would provide strong support for the link between accretion history and halo structure we propose here.

We use the Millennium Simulation series to investigate the mass and redshift dependence of the concentration of equilibrium cold dark matter (CDM) halos. We extend earlier work on the relation between halo mass profiles and assembly histories to show how the latter may be used to predict concentrations for halos of all masses and at any redshift. Our results clarify the link between concentration and the ``collapse redshift'' of a halo as well as why concentration depends on mass and redshift solely through the dimensionless ``peak height'' mass parameter, $\nu(M,z)=\delta_{\rm crit}(z)/\sigma(M,z)$. We combine these results with analytic mass accretion histories to extrapolate the $c(M,z)$ relations to mass regimes difficult to reach through direct simulation. Our model predicts that, at given $z$, $c(M)$ should deviate systematically from a simple power law at high masses, where concentrations approach a constant value, and at low masses, where concentrations are substantially lower than expected from extrapolating published empirical fits. This correction may reduce the expected self-annihilation boost factor from substructure by about one order of magnitude. The model also reproduces the $c(M,z)$ dependence on cosmological parameters reported in earlier work, and thus provides a simple and robust account of the relation between cosmology and the mass-concentration-redshift relation of CDM halos.

We examine the star formation properties of bright (~0.1 L*) satellites around isolated ~L* hosts in the local Universe using spectroscopically confirmed systems in the Sloan Digital Sky Survey DR7. Our selection method is carefully designed with the aid of N-body simulations to avoid groups and clusters. We find that satellites are significantly more likely to be quenched than a stellar mass-matched sample of isolated galaxies. Remarkably, this quenching occurs only for satellites of hosts that are themselves quenched: while star formation is unaffected in the satellites of star-forming hosts, satellites around quiescent hosts are more than twice as likely to be quenched than stellar-mass matched field samples. One implication of this is that whatever shuts down star formation in isolated, passive L* galaxies also plays at least an indirect role in quenching star formation in their bright satellites. The previously-reported tendency for "galactic conformity" in color/morphology may be a by-product of this host-specific quenching dichotomy. The S\'ersic indices of quenched satellites are statistically identical to those of field galaxies with the same specific star formation rates, suggesting that environmental and secular quenching give rise to the same morphological structure. By studying the distribution of pairwise velocities between the hosts and satellites, we find dynamical evidence that passive host galaxies reside in dark matter halos that are ~45% more massive than those of star-forming host galaxies of the same stellar mass. We emphasize that even around passive hosts, the mere fact that galaxies become satellites does not typically result in star formation quenching: we find that only ~30% of ~0.1 L* galaxies that fall in from the field are quenched around passive hosts, compared with ~0% around star forming hosts.

The Hubble Space Telescope (HST) has proven to be uniquely suited for the measurement of proper motions (PMs) of stars and galaxies in the nearby Universe. Here we summarize the main results and ongoing studies of the HSTPROMO collaboration, which over the past decade has executed some two dozen observational and theoretical HST projects on this topic. This is continuing to revolutionize our dynamical understanding of many objects, including: globular clusters; young star clusters; stars and stellar streams in the Milky Way halo; Local Group galaxies, including dwarf satellite galaxies, the Magellanic Clouds, and the Andromeda galaxy; and AGN Black Hole Jets.

We have combined the semi-analytic galaxy formation model of Guo et al. (2011) with the particle-tagging technique of Cooper et al. (2010) to predict galaxy surface brightness profiles in a representative sample of ~1900 massive dark matter haloes (10^12--10^14 M_sol) from the Millennium II Lambda-CDM N-body simulation. Here we present our method and basic results focusing on the outer regions of galaxies, consisting of stars accreted in mergers. These simulations cover scales from the stellar haloes of Milky Way-like galaxies to the 'cD envelopes' of groups and clusters, and resolve low surface brightness substructure such as tidal streams. We find that the surface density of accreted stellar mass around the central galaxies of dark matter haloes is well described by a Sersic profile, the radial scale and amplitude of which vary systematically with halo mass (M_200). The total stellar mass surface density profile breaks at the radius where accreted stars start to dominate over stars formed in the galaxy itself. This break disappears with increasing M_200 because accreted stars contribute more of the total mass of galaxies, and is less distinct when the same galaxies are averaged in bins of stellar mass, because of scatter in the relation between M_star and M_200. To test our model we have derived average stellar mass surface density profiles for massive galaxies at z~0.08 by stacking SDSS images. Our model agrees well with these stacked profiles and with other data from the literature, and makes predictions that can be more rigorously tested by future surveys that extend the analysis of the outer structure of galaxies to fainter isophotes. We conclude that it is likely that the outer structure of the spheroidal components of galaxies is largely determined by collisionless merging during their hierarchical assembly

The observed central densities of Milky Way dwarf spheroidal galaxies (dSphs) are significantly lower than the densities of the largest (Vmax about 35 km/s) subhalos found in dissipationless simulations of Galaxy-size dark matter hosts. One possible explanation is that gas removal from feedback can lower core densities enough to match observations. We model the dynamical effects of supernova feedback through the use of a time-varying central potential in high resolution, idealized numerical simulations and explore the resulting impact on the mass distributions of dwarf dark matter halos. We find that in order to match the observed central masses of M_star about 10^6 M_sun dSphs, the energy equivalent of more than 40,000 supernovae must be delivered with 100% efficiency directly to the dark matter. This energy requirement exceeds the number of supernovae that have ever exploded in most dSphs for typical initial mass functions. We also find that, per unit energy delivered and per cumulative mass removed from the galaxy, single blow-out events are more effective than repeated small bursts in reducing central dark matter densities. We conclude that it is unlikely that supernova feedback alone can solve the "Too Big to Fail" problem for Milky Way subhalos.

Stellar archeology shows that massive elliptical galaxies today formed rapidly about ten billion years ago with star formation rates above several hundreds solar masses per year (M_sun/yr). Their progenitors are likely the sub-millimeter-bright galaxies (SMGs) at redshifts (z) greater than 2. While SMGs' mean molecular gas mass of 5x10^10 M_sun can explain the formation of typical elliptical galaxies, it is inadequate to form ellipticals that already have stellar masses above 2x10^11 M_sun at z ~ 2. Here we report multi-wavelength high-resolution observations of a rare merger of two massive SMGs at z = 2.3. The system is currently forming stars at a tremendous rate of 2,000 M_sun/yr. With a star formation efficiency an order-of-magnitude greater than that of normal galaxies, it will quench the star formation by exhausting the gas reservoir in only ~200 million years. At a projected separation of 19 kiloparsecs, the two massive starbursts are about to merge and form a passive elliptical galaxy with a stellar mass of ~4x10^11 M_sun. Our observations show that gas-rich major galaxy mergers, concurrent with intense star formation, can form the most massive elliptical galaxies by z ~ 1.5.

Segue 2, discovered by Belokurov et al. (2009), is a galaxy with a luminosity of only 900 L_sun. We present Keck/DEIMOS spectroscopy of 25 members of Segue 2--a threefold increase in spectroscopic sample size. The velocity dispersion is too small to be measured with our data. The upper limit with 90% (95%) confidence is sigma_v < 2.2 (2.6) km/s, the most stringent limit for any galaxy. The corresponding limit on the mass within the 3-D half-light radius (46 pc) is M_1/2 < 1.5 (2.1) x 10^5 M_sun. Segue 2 is the least massive galaxy known. We identify Segue 2 as a galaxy rather than a star cluster based the wide dispersion in [Fe/H] (from -2.85 to -1.33) among the member stars. The stars' [alpha/Fe] ratios decline with increasing [Fe/H], indicating that Segue 2 retained Type Ia supernova ejecta despite its presently small mass and that star formation lasted for at least 100 Myr. The mean metallicity, <[Fe/H]> = -2.22 +/- 0.13 (about the same as the Ursa Minor galaxy, 330 times more luminous than Segue 2), is higher than expected from the luminosity-metallicity relation defined by more luminous dwarf galaxy satellites of the Milky Way. Segue 2 may be the barest remnant of a tidally stripped, Ursa Minor-sized galaxy. If so, it is the best example of an ultra-faint dwarf galaxy that came to be ultra-faint through tidal stripping. Alternatively, Segue 2 could have been born in a very low-mass dark matter subhalo (v_max < 10 km/s), below the atomic hydrogen cooling limit.

We compare spherically-averaged radial number counts of bright (> 10^5 Lsun) dwarf satellite galaxies within 400 kpc of the Milky Way (MW) and M31 and find that the MW satellites are much more centrally concentrated. Remarkably, the two satellite systems are almost identical within the central 100 kpc, while M31 satellites outnumber MW satellites by about a factor of four at deprojected distances spanning 100 - 400 kpc. We compare the observed distributions to those predicted for LCDM suhbalos using a suite of 44 high-resolution ~10^12 halo zoom simulations, 22 of which are in pairs like the MW and M31. We find that the radial distribution of satellites around M31 is fairly typical of those predicted for subhalos, while the Milky Way's distribution is more centrally concentrated that any of our simulated LCDM halos. One possible explanation is that our census is bright (> 10^5 Lsun) MW dwarf galaxies is significantly incomplete beyond ~ 100 kpc of the Sun. If there were ~8 - 20 more bright dwarfs orbiting undetected at 100 - 400 kpc, then the Milky Way's radial distribution would fall within the range expected from subhalo distributions and alos look very much like the known M31 system. We use our simulations to demonstrate that there is enough area left unexplored by the Sloan Digital Sky Survey and its extensions that the discovery of ~10 new bright dwarfs is not implausible given the expected range of angular anisotropy of subhalos in the sky.

We present the first absolute proper motion measurement of Leo I, based on two epochs of HST ACS/WFC images separated by ~5 years. The average shift of Leo I stars with respect to ~100 background galaxies implies a proper motion of (mu_W, mu_N) = (0.1140 +/- 0.0295, -0.1256 +/- 0.0293) mas/yr. The implied Galactocentric velocity vector, corrected for the reflex motion of the Sun, has radial and tangential components V_rad = 167.9 +/- 2.8 km/s and V_tan = 101.0 +/- 34.4 km/s, respectively. We study the detailed orbital history of Leo I by solving its equations of motion backward in time for a range of plausible mass models for the Milky Way and its surrounding galaxies. Leo I entered the Milky Way virial radius 2.33 +/- 0.21 Gyr ago, most likely on its first infall. It had a pericentric approach 1.05 +/- 0.09 Gyr ago at a Galactocentric distance of 91 +/- 36 kpc. We associate these time scales with characteristic time scales in Leo I's star formation history, which shows an enhanced star formation activity ~2 Gyr ago and quenching ~1 Gyr ago. There is no indication from our calculations that other galaxies have significantly influenced Leo I's orbit, although there is a small probability that it may have interacted with either Ursa Minor or Leo II within the last ~1 Gyr. For most plausible Milky Way masses, the observed velocity implies that Leo I is bound to the Milky Way. However, it may not be appropriate to include it in models of the Milky Way satellite population that assume dynamical equilibrium, given its recent infall. Solution of the complete (non-radial) timing equations for the Leo I orbit implies a Milky Way mass M_MW,vir = 3.15 (-1.36, +1.58) x 10^12 Msun, with the large uncertainty dominated by cosmic scatter. In a companion paper, we compare the new observations to the properties of Leo I subhalo analogs extracted from cosmological simulations.

We combine our Hubble Space Telescope measurement of the proper motion of the Leo I dwarf spheroidal galaxy (presented in a companion paper) with the highest resolution numerical simulations of Galaxy-size dark matter halos in existence to constrain the mass of the Milky Way's dark matter halo (M_MW). Despite Leo I's large Galacto-centric space velocity (200 km/s) and distance (261 kpc), we show that it is extremely unlikely to be unbound if Galactic satellites are associated with dark matter substructure, as 99.9% of subhalos in the simulations are bound to their host. The observed position and velocity of Leo I strongly disfavor a low mass Milky Way: if we assume that Leo I is the least bound of the Milky Way's classical satellites, then we find that M_MW > 10^{12} M_sun at 95% confidence for a variety of Bayesian priors on M_MW. In lower mass halos, it is vanishingly rare to find subhalos at 261 kpc moving as fast as Leo I. Should an additional classical satellite be found to be less bound than Leo I, this lower limit on M_MW would increase by 30%. Imposing a mass weighted LCDM prior, we find a median Milky Way virial mass of M_MW=1.6 x 10^{12} M_sun, with a 90% confidence interval of [1.0-2.4] x 10^{12} M_sun. We also confirm a strong correlation between subhalo infall time and orbital energy in the simulations and show that proper motions can aid significantly in interpreting the infall times and orbital histories of satellites.

Using the technique of Angulo & White (2010) we scale the Millennium and Millennium-II simulations of structure growth in a LCDM universe from the cosmological parameters with which they were carried out (based on first-year results from the Wilkinson Microwave Anisotropy Probe, WMAP1) to parameters consistent with the seven-year WMAP data (WMAP7). We implement semi-analytic galaxy formation modelling on both simulations in both cosmologies to investigate how the formation, evolution and clustering of galaxies are predicted to vary with cosmological parameters. The increased matter density Omega_m and decreased linear fluctuation amplitude sigma8 in WMAP7 have compensating effects, so that the abundance and clustering of dark halos are predicted to be very similar to those in WMAP1 for z <= 3. As a result, local galaxy properties can be reproduced equally well in the two cosmologies by slightly altering galaxy formation parameters. The evolution of the galaxy populations is then also similar. In WMAP7, structure forms slightly later. This shifts the peak in cosmic star formation rate to lower redshift, resulting in slightly bluer galaxies at z=0. Nevertheless, the model still predicts more passive low-mass galaxies than are observed. For rp< 1Mpc, the z=0 clustering of low-mass galaxies is weaker for WMAP7 than for WMAP1 and closer to that observed, but the two cosmologies give very similar results for more massive galaxies and on large scales. At z>1 galaxies are predicted to be more strongly clustered for WMAP7. Differences in galaxy properties, including, clustering, in these two cosmologies are rather small up to redshift 3. Given that there are still considerable residual uncertainties in galaxy formation models, it is very difficult to distinguish WMAP1 from WMAP7 through observations of galaxy properties or their evolution.

We use the Millennium Simulation series to study how the dynamical state of dark matter halos affects the relation between mass and concentration. We find that a large fraction of massive systems are identified when they are substantially out of equilibrium and in a particular phase of their dynamical evolution: the more massive the halo, the more likely it is found at a transient stage of high concentration. This state reflects the recent assembly of massive halos and corresponds to the first pericentric passage of recently-accreted material when, before virialization, the kinetic and potential energies reach maximum and minimum values, respectively. This result explains the puzzling upturn in the mass-concentration relation reported in recent work for massive halos; indeed, the upturn disappears when only dynamically-relaxed systems are considered in the analysis. Our results warn against applying simple equilibrium models to describe the structure of rare, massive galaxy clusters and urges caution when extrapolating scaling laws calibrated on lower-mass systems, where such deviations from equilibrium are less common. The evolving dynamical state of galaxy clusters ought to be carefully taken into account if cluster studies are to provide precise cosmological constraints.

We study the impact of star formation and stellar feedback prescriptions on galaxy properties predicted by means of "stripped-down" versions of independently developed semi-analytic models (SAMs). These include cooling, star formation, feedback from supernovae (SNe) and simplified prescriptions for galaxy merging, but no chemical evolution, disc instabilities or AGN feedback. We run these versions on identical samples of dark matter (DM) haloes extracted from high-resolution N-body simulations in order to perform both statistical analysis and object-by-object comparisons. We compare our results with previous work based on stripped-down versions of the same SAMs including only gas cooling, and show that all feedback models provide coherent modifications in the distribution of baryons between the various gas phases. In particular, we find that the predicted hot gas fractions are considerably increased by up to a factor of three, while the corresponding cold gas fractions are correspondingly decreased, and a significant amount of mass is ejected from the DM halo. Nonetheless, we also find relevant differences in the predicted properties of model galaxies among the three SAMs: these deviations are more relevant at mass scales comparable to that of our own Galaxy, and are reduced at larger masses, confirming the varying impact of stellar feedback at different mass scales. We also check the effect of enhanced star formation events (i.e. starbursts modes), defined in connection with galaxy mergers. We find that, in general, these episodes have a limited impact in the overall star formation histories of model galaxies, even in massive DM halos where merger-driven star formation has often been considered very important.

The Milky Way appears to be missing baryons, as the observed mass in stars and gas is well below the cosmic mean. One possibility is that a substantial fraction of the Galaxy's baryons are embedded within an extended, million-degree hot halo, an idea supported indirectly by observations of warm gas clouds in the halo and gas-free dwarf spheroidal satellites. X-ray observations have established that hot gas does exist in our Galaxy beyond the local hot bubble; however, it may be distributed in a hot disk configuration. Moreover, recent investigations into the X-ray constraints have suggested that any Galactic corona must be insignificant. Here we re-examine the observational data, particularly in the X-ray and radio bands, in order to determine whether it is possible for a substantial fraction of the Galaxy's baryons to exist in ~ 10^6 K gas. In agreement with past studies, we find that a baryonically closed halo is clearly ruled out if one assumes that the hot corona is distributed with a cuspy NFW profile. However, if the hot corona of the galaxy is in an extended, low-density distribution with a large central core, as expected for an adiabatic gas in hydrostatic equilibrium, then it may contain up to 10^11 M_sun of material, possibly accounting for all of the missing Galactic baryons. We briefly discuss some potential avenues for discriminating between a massive, extended hot halo and a local hot disk.

We use the Millennium I and II cosmological simulations to revisit the impact of mergers in the growth of bulges in central galaxies in the LCDM scenario. We seed galaxies within the growing CDM haloes using semi-empirical relations to assign stellar and gaseous masses, and an analytic treatment to estimate the transfer of stellar mass to the bulge of the remnant after a galaxy merger. We find that this model roughly reproduces the observed correlation between the bulge-to-total (B/T) mass ratio and stellar mass in present-day central galaxies as well as their observed demographics, although low-mass B/T<0.1 (bulgeless) galaxies might be scarce relative to the observed abundance. In our merger-driven scenario, bulges have a composite population made of (i) stars acquired from infalling satellites, (ii) stars transferred from the primary disc due to merger-induced perturbations, and (iii) newly formed stars in starbursts triggered by mergers. We find that (i) and (ii) are the main channels of mass assembly, with the first being dominant for massive galaxies, creating large bulges with different stellar populations than those of the inner discs, while the second is dominant for intermediate/low-mass galaxies creating small bulges with similar stellar populations to the inner discs. We associate the dominion of the first (second) channel to classical (pseudo) bulges, and compare the predicted fractions to observations. We remark that our treatment does not include other mechanisms of bulge growth such as intrinsic secular disc instabilities or misaligned gas accretion. We find that the evolution of the stellar and gaseous contents of the satellite as it moves towards the central galaxy is a key ingredient in setting the morphology of the remnant, and that a good match to the observed bulge demographics occurs when this evolution proceeds closely to that of the central galaxy.

We use the Aquarius simulations to show that the most massive subhalos in galaxy-mass dark matter halos in LCDM are grossly inconsistent with the dynamics of the brightest Milky Way dwarf spheroidal galaxies. While the best-fitting hosts of the dwarf spheroidals all have 12 < Vmax < 25 km/s, LCDM simulations predict at least ten subhalos with Vmax > 25 km/s. These subhalos are also among the most massive at earlier times, and significantly exceed the UV suppression mass back to z ~ 10. No LCDM-based model of the satellite population of the Milky Way explains this result. The problem lies in the satellites' densities: it is straightforward to match the observed Milky Way luminosity function, but doing so requires the dwarf spheroidals to have dark matter halos that are a factor of ~5 more massive than is observed. Independent of the difficulty in explaining the absence of these dense, massive subhalos, there is a basic tension between the derived properties of the bright Milky Way dwarf spheroidals and LCDM expectations. The inferred infall masses of these galaxies are all approximately equal and are much lower than standard LCDM predictions for systems with their luminosities. Consequently, their implied star formation efficiencies span over two orders of magnitude, from 0.2% to 20% of baryons converted into stars, in stark contrast with expectations gleaned from more massive galaxies. We explore possible solutions to these problems within the context of LCDM and find them to be unconvincing. In particular, we use controlled simulations to demonstrate that the small stellar masses of the bright dwarf spheroidals make supernova feedback an unlikely explanation for their low inferred densities.

We combine the high-resolution Aquarius simulations with three-dimensional models of reionization based on the initial density field of the Aquarius parent simulation, Millennium-II, to study the impact of patchy reionization on the faint satellite population of Milky Way halos. Because the Aquarius suite consists of zoom-in simulations of halos in the Millennium-II volume, we follow the formation of substructure and the growth of reionization bubbles due to the larger environment simultaneously, and thereby determine the reionization redshifts of satellite candidates. We do this for four different reionization models, and also compare results to instantaneous reionization. Using a simple procedure for selecting satellites and assigning luminosities in the simulations, we compare the resulting satellite populations. We find that the overall number of satellites depends sensitively on the reionization model, with a factor of 3-4 variation between the four models for a given host halo, although the difference is entirely in the population of faint satellites (M_V > -10). In addition, we find that for a given reionization model the total number of satellites differs by 10%-20% between the patchy and homogeneous scenarios, provided that the redshift is chosen appropriately for the instantaneous case. However, the halo-halo scatter from the six Aquarius halos is large, up to a factor of 2-3, and so is comparable to the difference between reionization scenarios. In order to use the population of faint dwarf galaxies around the Milky Way as a probe of the local reionization history, then, it is necessary to first better understand the general distribution of substructure around Milky Way-mass halos.

Dark matter halo merger trees are now routinely extracted from cosmological simulations of structure formation. These trees are frequently used as inputs to semi-analytic models of galaxy formation to provide the backbone within which galaxy formation takes place. By necessity, these merger trees are constructed from a finite set of discrete "snapshots" of the N-body simulation and so have a limited temporal resolution. To date, there has been little consideration of how this temporal resolution affects the properties of galaxies formed within these trees. In particular, the question of how many snapshots are needed to achieve convergence in galaxy properties has not be answered. Therefore, we study the convergence in the stellar and total baryonic masses of galaxies, distribution of merger times, stellar mass functions and star formation rates in the Galacticus model of galaxy formation as a function of the number of "snapshot" times used to represent dark matter halo merger trees. When utilizing snapshots between z=20 and z=0, we find that at least 128 snapshots are required to achieve convergence to within 5% for galaxy masses. This convergence is obtained for mean quantities averaged over large samples of galaxies - significant variance for individual galaxies remains even when using very large numbers of snapshots. We find only weak dependence of the rate of convergence on the distribution of snapshots in time - snapshots spaced uniformly in the expansion factor, uniformly in the logarithm of expansion factor or uniformly in the logarithm of critical overdensity for collapse work equally well in almost all cases. We provide input parameters to Galacticus which allow this type of convergence study to be tuned to other simulations and to be carried out for other galaxy properties.

We develop a new empirical methodology to study the relation between the stellar mass of galaxies and the mass of their host subhaloes. Our approach is similar to abundance matching, and is based on assigning a stellar mass to each subhalo within a large cosmological N-body simulation. However, we significantly extend the method to include a different treatment for central and satellite galaxies, allowing the stellar mass of satellite galaxies to depend on both the host subhalo mass, and on the halo mass. Unlike in previous studies, our models are constrained by two observations: the stellar mass function of galaxies, and their auto-correlation functions (CFs). We present results for ~10^6 different successful models, showing that the uncertainty in the host subhalo mass reaches a factor of ~10 for a given stellar mass. Our results thus indicate that the host subhalo mass of central and satellite galaxies is poorly constrained by using information coming solely from the abundance and CFs of galaxies. In addition, we demonstrate that the fraction of stellar mass relative to the universal baryon fraction can vary between different models, and can reach ~0.6 for a specific population of low mass galaxies. We conclude that additional observational constraints are needed, in order to better constrain the mass relation between haloes and galaxies. These might include weak lensing, satellite kinematics, or measures of clustering other than the CFs. Moreover, CFs based on larger surveys with a better estimate of the systematic uncertainties are needed.

We use a volume-limited spectroscopic sample of isolated galaxies in the Sloan Digital Sky Survey (SDSS) to investigate the frequency and radial distribution of luminous (M_r <~ -18.3) satellites like the Large Magellanic Cloud (LMC) around ~L* Milky Way analogs and compare our results object-by-object to LCDM predictions based on abundance matching in simulations. We show that 12% of Milky Way-like galaxies host an LMC-like satellite within 75 kpc (projected), and 42 % within 250 kpc (projected). This implies ~10% have a satellite within the distance of the LMC, and ~40% of L* galaxies host a bright satellite within the virialized extent of their dark matter halos. Remarkably, the simulation reproduces the observed frequency, radial dependence, velocity distribution, and luminosity function of observed secondaries exceptionally well, suggesting that LCDM provides an accurate reproduction of the observed Universe to galaxies as faint as L~10^9 Lsun on ~50 kpc scales. When stacked, the observed projected pairwise velocity dispersion of these satellites is sigma~160 km/s, in agreement with abundance-matching expectations for their host halo masses. Finally, bright satellites around L* primaries are significantly redder than typical galaxies in their luminosity range, indicating that environmental quenching is operating within galaxy-size dark matter halos that typically contain only a single bright satellite. This redness trend is in stark contrast to the Milky Way's LMC, which is unusually blue even for a field galaxy. We suggest that the LMC's discrepant color might be further evidence that it is undergoing a triggered star-formation event upon first infall.

We show that dissipationless LCDM simulations predict that the majority of the most massive subhaloes of the Milky Way are too dense to host any of its bright satellites (L_V > 10^5 L_sun). These dark subhaloes have circular velocities at infall of 30-70 km/s and infall masses of [0.2-4] x 10^10 M_sun. Unless the Milky Way is a statistical anomaly, this implies that galaxy formation becomes effectively stochastic at these masses. This is in marked contrast to the well-established monotonic relation between galaxy luminosity and halo circular velocity (or halo mass) for more massive haloes. We show that at least two (and typically four) of these massive dark subhaloes are expected to produce a larger dark matter annihilation flux than Draco. It may be possible to circumvent these conclusions if baryonic feedback in dwarf satellites or different dark matter physics can reduce the central densities of massive subhaloes by order unity on a scale of 0.3 - 1 kpc.

We examine Milky Way-Magellanic Cloud systems selected from the Millennium-II Simulation in order to place the orbits of the Magellanic Clouds in a cosmological context. Our analysis shows that satellites massive enough to be LMC analogs are typically accreted at late times. Moreover, those that are accreted at early times and survive to the present have orbital properties that are discrepant with those observed for the LMC. The high velocity of the LMC, coupled with the dearth of unbound orbits seen in the simulation, argues that the mass of the MW's halo is unlikely to be less than 2 x 10^12 Msun. This conclusion is further supported by statistics of halos hosting satellites with masses, velocities, and separations comparable to those of the LMC. We further show that: (1) LMC and SMC-mass objects are not particularly uncommon in MW-mass halos; (2) the apparently high angular momentum of the LMC is not cosmologically unusual; and (3) it is rare for a MW halo to host a LMC-SMC binary system at z=0, but high speed binary pairs accreted at late times are possible. Based on these results, we conclude that the LMC was accreted within the past four Gyr and is currently making its first pericentric passage about the MW.

Recent studies emphasize that an empirical relation between the stellar mass of galaxies and the mass of their host dark matter subhaloes can predict the clustering of galaxies and its evolution with cosmic time. In this paper we study the assumptions made by this methodology using a semi-analytical model (SAM). To this end, we randomly swap between the locations of model galaxies within a narrow range of subhalo mass (M_infall). We find that shuffled samples of galaxies have different auto-correlation functions in comparison with the original model galaxies. This difference is significant even if central and satellite galaxies are allowed to follow a different relation between M_infall and stellar mass, and can reach a factor of 2 for massive galaxies at redshift zero. We analyze three features within SAMs that contribute to this effect: a) The relation between stellar mass and subhalo mass evolves with redshift for central galaxies, affecting satellite galaxies at the time of infall. b) The stellar mass of galaxies falling into groups and clusters at high redshift is different from the mass of central galaxies at the same time. c) The stellar mass growth for satellite galaxies after infall can be significant and depends on the infall redshift and the group mass. We show that the above is true for differing SAMs, and that the effect is sensitive to the treatment of dynamical friction and stripping of gas in satellite galaxies. We find that by using the FoF group mass at redshift zero in addition to M_infall, an empirical model is able to accurately reproduce the clustering properties of galaxies. On the other hand, using the infall redshift as a second parameter does not yield as good results because it is less correlated with stellar mass. Our analysis indicates that environmental processes are important for modeling the clustering and abundance of galaxies. (Abridged)

Cosmological N-body simulations indicate that the spherically-averaged density profiles of cold dark matter halos are accurately described by Einasto profiles, where the logarithmic slope is a power-law of adjustable exponent, \gamma =dln\rho /dlnr ~ r^\alpha $. The pseudo-phase-space density (PPSD) profiles of CDM halos also show remarkable regularity, and are well approximated by simple power laws, Q(r)=\rho /\sigma ^3 ~ r^-\chi . We show that this is expected from dynamical equilibrium considerations, since Jeans' equations predict that the pseudo-phase-space density profiles of Einasto halos should resemble power laws over a wide range of radii. For the values of \alpha typical of CDM halos, the inner Q profiles of equilibrium halos deviate significantly from a power law only very close to the center, and simulations of extremely high-resolution would be needed to detect such deviations unambiguously. We use an ensemble of halos drawn from the Millennium-II simulation to study which of these two alternatives describe best the mass profile of CDM halos. Our analysis indicates that at the resolution of the best available simulations, both Einasto and power-law PPSD profiles (with adjustable exponents \alpha and \chi, respectively) provide equally acceptable fits to the simulations. A full account of the structure of CDM halos requires understanding how the shape parameters that characterize departures from self-similarity, like \alpha or \chi, are determined by evolutionary history, environment or initial conditions.

We apply updated semi-analytic galaxy formation models simultaneously to the stored halo/subhalo merger trees of the Millennium and Millennium-II simulations. These differ by a factor of 125 in mass resolution, allowing explicit testing of resolution effects on predicted galaxy properties. We have revised the treatments of the transition between the rapid infall and cooling flow regimes of gas accretion, of the sizes of bulges and of gaseous and stellar disks, of supernova feedback, of the transition between central and satellite status as galaxies fall into larger systems, and of gas and star stripping once they become satellites. Plausible values of efficiency and scaling parameters yield an excellent fit not only to the observed abundance of low-redshift galaxies over 5 orders of magnitude in stellar mass and 9 magnitudes in luminosity, but also to the observed abundance of Milky Way satellites. This suggests that reionisation effects may not be needed to solve the "missing satellite" problem except, perhaps, for the faintest objects. The same model matches the observed large-scale clustering of galaxies as a function of stellar mass and colour. The fit remains excellent down to ~30kpc for massive galaxies. For M* < 6 x 10^10Msun, however, the model overpredicts clustering at scales below 1 Mpc, suggesting that the sigma_8 adopted in the simulations (0.9) is too high. Galaxy distributions within rich clusters agree between the simulations and match those observed, but only if galaxies without dark matter subhalos (so-called orphans) are included. Our model predicts a larger passive fraction among low-mass galaxies than is observed, as well as an overabundance of ~10^10Msun galaxies beyond z~0.6, reflecting deficiencies in the way star-formation rates are modelled.

We study the population statistics of the surviving subhaloes of LCDM dark matter haloes using a set of very high resolution N-body simulations. These include both simulations of representative regions of the Universe and ultra-high resolution resimulations of individual dark matter haloes. We find that more massive haloes tend to have a larger mass fraction in subhaloes. For example, cluster size haloes typically have 7.5 percent of their mass within R200 in substructures of fractional mass larger than 1e-5, which is 25 percent higher than galactic haloes. There is, however, a large variance in the subhalo mass fraction from halo to halo, whereas the subhalo abundance shows much higher regularity. For dark matter haloes of fixed mass, the subhalo abundance decreases by 30 percent between redshift 2 and 0. The subhalo abundance function correlates with the host halo concentration parameter and formation redshift. However, the intrinsic scatter is not significantly reduced for narrow ranges of concentration parameter or formation redshift, showing that they are not the dominant parameters that determine the subhalo abundance in a halo.

We present an analysis of the distribution of structural properties for Milky Way-mass halos in the Millennium-II Simulation (MS-II). This simulation of structure formation within the standard LCDM cosmology contains thousands of Milky Way-mass halos and has sufficient resolution to properly resolve many subhalos per host. It thus provides a major improvement in the statistical power available to explore the distribution of internal structure for halos of this mass. In addition, the MS-II contains lower resolution versions of the Aquarius Project halos, allowing us to compare our results to simulations of six halos at a much higher resolution. We study the distributions of mass assembly histories, of subhalo mass functions and accretion times, and of merger and stripping histories for subhalos capable of impacting disks at the centers of halos. We show that subhalo abundances are not well-described by Poisson statistics at low mass, but rather are dominated by intrinsic scatter. Using the masses of subhalos at infall and the abundance-matching assumption, there is less than a 10% chance that a Milky Way halo with M_vir =10^12 M_sun will host two galaxies as bright as the Magellanic Clouds. This probability rises to ~25% for a halo with M_vir=2.5 x 10^12 M_sun. The statistics relevant for disk heating are very sensitive to the mass range that is considered relevant. Mergers with infall mass : redshift zero virial mass greater than 1:30 could well impact a central galactic disk and are a near inevitability since z=2, whereas only half of all halos have had a merger with infall mass : redshift zero virial mass greater than 1:10 over this same period.

We construct merger trees of dark matter haloes and quantify their merger rates and mass growth rates using the joint dataset from the Millennium and Millennium-II simulations. The finer resolution of the Millennium-II Simulation has allowed us to extend our earlier analysis of halo merger statistics to an unprecedentedly wide range of descendant halo mass (10^10 < M0 < 10^15 Msun), progenitor mass ratio (10^-5 < xi < 1), and redshift (0 < z < 15). We update our earlier fitting form for the mean merger rate per halo as a function of M_0, xi, and z. The overall behavior of this quantity is unchanged: the rate per unit redshift is nearly independent of z out to z~15; the dependence on halo mass is weak (M0^0.13); and it is nearly a power law in the progenitor mass ratio (xi^-2). We also present a simple and accurate fitting formula for the mean mass growth rate of haloes as a function of mass and redshift. This mean rate is 46 Msun/yr for 10^12 Msun haloes at z=0, and it increases with mass as M^{1.1} and with redshift as (1+z)^2.5 (for z > 1). When the fit for the mean mass growth rate is integrated over a halo's history, we find excellent match to the mean mass assembly histories of the simulated haloes. By combining merger rates and mass assembly histories, we present results for the number of mergers over a halo's history and the statistics of the redshift of the last major merger.

We use stripped-down versions of three semi-analytic galaxy formation models to study the influence of different assumptions about gas cooling and galaxy mergers. By running the three models on identical sets of merger trees extracted from high-resolution cosmological N-body simulations, we are able to perform both statistical analyses and halo-by-halo comparisons. Our study demonstrates that there is a good statistical agreement between the three models used here, when operating on the same merger trees, reflecting a general agreement in the underlying framework for semi-analytic models. We also show, however, that various assumptions that are commonly adopted to treat gas cooling and galaxy mergers can lead to significantly different results, at least in some regimes. In particular, we find that the different models adopted for gas cooling lead to similar results for mass scales comparable to that of our own Galaxy. Significant differences, however, arise at larger mass scales. These are largely (but not entirely) due to different treatments of the `rapid cooling' regime, and different assumptions about the hot gas distribution. At this mass regime, the predicted cooling rates can differ up to about one order of magnitude, with important implications on the relative weight that these models give to AGN feedback in order to counter-act excessive gas condensation in relatively massive haloes at low redshift. Different assumptions in the modelling of galaxy mergers can also result in significant differences in the timings of mergers, with important consequences for the formation and evolution of massive galaxies.

We discuss how the effective radius Phi(Re) function (ERF) recently worked out by Bernardi et al. (2009) represents a new testbed to improve the current understanding of Semi-analytic Models of Galaxy formation. In particular, we here show that a detailed hierarchical model of structure formation can broadly reproduce the correct peak in the size distribution of local early-type galaxies, although it significantly overpredicts the number of very compact and very large galaxies. This in turn is reflected in the predicted size-mass relation, much flatter than the observed one, due to too large (~3 kpc) low-mass galaxies (<10^11 \msun), and to a non-negligible fraction of compact (< 0.5-1 kpc) and massive galaxies (> 10^11 \msun). We also find that the latter discrepancy is smaller than previously claimed, and limited to only ultracompact (Re < 0.5 kpc) galaxies when considering elliptical-dominated samples. We explore several causes behind these effects. We conclude that the former problem might be linked to the initial conditions, given that large and low-mass galaxies are present at all epochs in the model. The survival of compact and massive galaxies might instead be linked to their very old ages and peculiar merger histories. Overall, knowledge of the galactic stellar mass {\em and} size distributions allows a better understanding of where and how to improve models.

If dark matter is composed of neutralinos, one of the most exciting prospects for its detection lies in observations of the gamma-ray radiation created in pair annihilations between neutralinos, a process that may contribute significantly to the extragalactic gamma-ray background (EGB) radiation. We here use the high-resolution Millennium-II simulation of cosmic structure formation to produce the first full-sky maps of the expected radiation coming from extragalactic dark matter structures. Our map making procedure takes into account the total gamma-ray luminosity from all haloes and their subhaloes, and includes corrections for unresolved components of the emission as well as an extrapolation to the damping scale limit of neutralinos. Our analysis also includes a proper normalization of the signal according to a specific supersymmetric model based on minimal supergravity. The new simulated maps allow a study of the angular power spectrum of the gamma-ray background from dark matter annihilation, which has distinctive features associated with the nature of the annihilation process and may be detectable in forthcoming observations by the recently launched FERMI satellite. Our results are in broad agreement with analytic models for the gamma-ray background, but they also include higher-order correlations not readily accessible in analytic calculations and, in addition, provide detailed spectral information for each pixel. In particular, we find that difference maps at different energies can reveal cosmic large-scale structure at low and intermediate redshifts. If the intrinsic emission spectrum is characterized by an emission peak, cosmological tomography with gamma ray annihilation radiation is in principle possible.

If dark matter is composed of neutralinos, the gamma-ray radiation produced in their annihilation offers an attractive possibility for dark matter detection. This process may contribute significantly to the extragalactic gamma-ray background (EGB) radiation, which is being measured by the FERMI satellite with unprecedented sensitivity. Using the high-resolution Millennium-II simulation of cosmic structure formation we have produced the first full-sky maps of the expected contribution of dark matter annihilation to the EGB radiation. Our maps include a proper normalization of the signal according to a specific supersymmetric model based on minimal supergravity. The new simulated maps allow a study of the angular power spectrum of the gamma-ray background from dark matter annihilation, which has distinctive features associated with the nature of the annihilation process. Our results are in broad agreement with analytic models for the gamma-ray background, but they also include higher-order correlations not readily accessible in analytic calculations and, in addition, provide detailed spectral information for each pixel. In particular, we find that color maps combining different energies can reveal the cosmic large-scale structure at low and intermediate redshifts.

For any assumed stellar Initial Mass Function, the Sloan Digital Sky Survey (SDSS) gives a precise determination of the stellar mass function of galaxies for 10^8 M_sun < M_* < 10^12 M_sun. Within the concordance LCDM cosmology, the Millennium simulations give a precise halo mass function for all halos within which galaxies can form. Under the plausible hypothesis that the stellar mass of a galaxy is an increasing function of the maximum mass ever attained by its halo, these combine to give halo mass as a function of stellar mass. The result agrees quite well with observational estimates of mean halo mass as a function of stellar mass from stacking analyses of the gravitational lensing signal and the satellite dynamics of SDSS galaxies. For M_* ~ 5.5 x 10^10 M_sun, the stellar mass usually assumed for the Milky Way, the implied halo mass is ~ 2 x 10^12 M_sun, consistent with most recent direct estimates and inferences from the MW/M31 Timing Argument. The fraction of the baryons associated with each halo which are present as stars in its central galaxy reaches a maximum of 20% at masses somewhat below that of the Milky Way, and falls rapidly at both higher and lower masses. These conversion efficiencies are lower than in almost all recent high-resolution simulations of galaxy formation, showing that these are not yet viable models for the formation of typical members of the galaxy population. When inserted in the Millennium-II Simulation, our derived relation between stellar mass and halo mass predicts a stellar mass autocorrelation function in excellent agreement with that measured directly in the SDSS. The implied Tully-Fisher relation also appears consistent with observation, suggesting that galaxy luminosity functions and Tully-Fisher relations can be reproduced simultaneously in a LCDM cosmology.

We present the Millennium-II Simulation (MS-II), a very large N-body simulation of dark matter evolution in the concordance LCDM cosmology. The MS-II assumes the same cosmological parameters and uses the same particle number and output data structure as the original Millennium Simulation (MS), but was carried out in a periodic cube one-fifth the size (100 Mpc/h) with 5 times better spatial resolution (a Plummer equivalent softening of 1.0 kpc/h) and with 125 times better mass resolution (a particle mass of 6.9 \times 10^6 Msun/h). By comparing results at MS and MS-II resolution, we demonstrate excellent convergence in dark matter statistics such as the halo mass function, the subhalo abundance distribution, the mass dependence of halo formation times, the linear and nonlinear autocorrelations and power spectra, and halo assembly bias. Together, the two simulations provide precise results for such statistics over an unprecedented range of scales, from halos similar to those hosting Local Group dwarf spheroidal galaxies to halos corresponding to the richest galaxy clusters. The "Milky Way" halos of the Aquarius Project were selected from a lower resolution version of the MS-II and were then resimulated at much higher resolution. As a result, they are present in the MS-II along with thousands of other similar mass halos. A comparison of their assembly histories in the MS-II and in resimulations of 1000 times better resolution shows detailed agreement over a factor of 100 in mass growth. We publicly release halo catalogs and assembly trees for the MS-II in the same format within the same archive as those already released for the MS.

The timescale for galaxies within merging dark matter halos to merge with each other is an important ingredient in galaxy formation models. Accurate estimates of merging timescales are required for predictions of astrophysical quantities such as black hole binary merger rates, the build-up of stellar mass in central galaxies, and the statistical properties of satellite galaxies within dark matter halos. In this paper, we study the merging timescales of extended dark matter halos using N-body simulations. We compare these results to standard estimates based on the Chandrasekhar theory of dynamical friction. We find that these standard predictions for merging timescales, which are often used in semi-analytic galaxy formation models, are systematically shorter than those found in simulations. The discrepancy is approximately a factor of 1.7 for $M_sat/M_host \approx 0.1$ and becomes larger for more disparate satellite-to-host mass ratios, reaching a factor of $\sim 3.3$ for $M_sat/M_host\approx 0.01$. Based on our simulations, we propose a new, easily implementable fitting formula that accurately predicts the timescale for an extended satellite to sink from the virial radius of a host halo down to the halo's center for a wide range of $M_sat/M_host$ and orbits. Including a central bulge in each galaxy changes the merging timescale by $\la 10%$. To highlight one concrete application of our results, we show that merging timescales often used in the literature overestimate the growth of stellar mass by satellite accretion by $\approx 40 %$, with the extra mass gained in low mass ratio mergers.

Minor mergers of galaxies are expected to be common in a hierarchical cosmology such as $\Lambda$CDM and have the potential to significantly affect galactic structure. In this paper we dissect the case-by-case outcome from a set of numerical simulations of a single satellite elliptical galaxy accreting onto a massive elliptical galaxy. We take care to explore cosmologically relevant orbital parameters and to set up realistic initial galaxy models that include all three relevant dynamical components: dark matter halos, stellar bulges, and central massive black holes. The effects of several different parameters are considered, including orbital energy and angular momentum, satellite density and inner density profile, satellite-to-host mass ratio, and presence of a black hole at the center of the host. Black holes play a crucial role in protecting the shallow stellar cores of the hosts, as satellites merging onto a host with a central black hole are more strongly disrupted than those merging onto hosts without black holes. Orbital parameters play an important role in determining the degree of disruption: satellites on less bound or more eccentric orbits are more easily destroyed than those on more bound or more circular orbits as a result of an increased number of pericentric passages and greater cumulative effects of gravitational shocking and tidal stripping. In addition, satellites with densities typical of faint elliptical galaxies are disrupted relatively easily, while denser satellites can survive much better in the tidal field of the host. Over the range of parameters explored, we find that the accretion of a single satellite elliptical galaxy can result in a broad variety of changes, in both signs, in the surface brightness profile and color of the central part of an elliptical galaxy.

Several recent observations suggest that gas-poor (dissipationless) mergers of elliptical galaxies contribute significantly to the build-up of the massive end of the red sequence. We perform a series of major merger simulations to investigate the spatial and velocity structure of the remnants of such mergers. Regardless of orbital energy or angular momentum, we find that the stellar remnants lie on the fundamental plane defined by their progenitors, a result of virial equilibrium with a small tilt due to an increasing central dark matter fraction. However, the locations of merger remnants in the projections of the fundamental plane -- the Faber-Jackson and R_e-M_* relations -- depend strongly on the merger orbit, and the relations steepen significantly from the canonical scalings (L sigma^4 and R_e M_*^0.6) for mergers on radial orbits. Our results imply that the projections of the fundamental plane -- but not necessarily the plane itself -- provide a powerful way of investigating the assembly history of massive elliptical galaxies, including the brightest cluster galaxies at or near the centers of galaxy clusters. We argue that most massive ellipticals are formed by anisotropic merging and that their fundamental plane projections should thus differ noticeably from those of lower mass ellipticals even though they should lie on the same fundamental plane. Current observations are consistent with this conclusion. The steepening in the L-sigma relation for luminous ellipticals may also be reflected in a corresponding steepening in the M_BH-sigma relation for massive black holes.

We carry out numerical simulations of dissipationless major mergers of elliptical galaxies using initial galaxy models that consist of a dark matter halo and a stellar bulge with properties consistent with the observed fundamental plane. By varying the density profile of the dark matter halo (standard NFW versus adiabatically contracted NFW), the global stellar to dark matter mass ratio, and the orbit of the merging galaxies, we are able to assess the impact of each of these factors on the structure of the merger remnant. Our results indicate that the properties of the remnant bulge depend primarily on the angular momentum and energy of the orbit; for a cosmologically motivated orbit, the effective radius and velocity dispersion of the remnant bulge remain approximately on the fundamental plane. This indicates that the observed properties of elliptical galaxies are consistent with significant growth via late dissipationless mergers. We also find that the dark matter fraction within the effective radius of our remnants increases after the merger, consistent with the hypothesis that the tilt of the fundamental plane from the virial theorem is due to a varying dark matter fraction as a function of galaxy mass.

Anisotropic gravitational radiation from a coalescing black hole binary can impart a recoil velocity of up to several hundred km/s to the remnant black hole. We examine the effects of recoiling massive black holes on their host stellar bulges, both for holes that escape their host and those that return to the galactic center via dynamical friction. We show that removal of a black hole via radiation recoil generally results in a rapidly-formed central core in the stellar system, with the effect being largest when the hole stays bound to the bulge and the recoil velocity is comparable to the bulge velocity dispersion. Black hole recoil therefore provides a mechanism for producing cores in some early type galaxies, but it is expected to be most efficient in faint ellipticals that are known to have steep density profiles. We argue that these results may hint at a significant role for gas in facilitating the coalescence of binary black holes in faint (power-law) early-type galaxies.

Much recent discussion about dark matter has been centered on two seemingly independent problems: the abundance of substructure in dark matter halos, and the cuspiness of the halos' inner density profile. We explore possible connections between the two problems by studying the gravitational scattering effects due to subhalos on the phase-space distribution of dark matter particles in the main halos. Our series of controlled numerical experiments indicates that the number and mass density of subhalos can be high enough to cause the collisionless dark matter particles in the inner part of a main halo to diffuse, flattening the main halo's inner cusp within a few dynamical times. Depending on the masses and concentration of the subhalos, the inner density profile of the whole system (main plus sub halos) can either steepen or flatten. Subhalo accretion can therefore introduce significant scatter in the inner density profiles of dark matter halos, offering a possible explanation for the range of profiles seen in both observations and cosmological simulations.

We present the results from a series of collisionless N-body simulations of major mergers of galaxy dark matter haloes with density profiles having either inner cusps or cores. Our simulations range from 2x10^5 to 10^7 particles, allowing us to probe the phase-space distribution of dark matter particles in the innermost regions (less than 0.005 virial radii) of cold dark matter haloes, a subject of much recent debate. We find that a major merger of two cored haloes yields a cored halo and does not result in a cuspy profile seen in many cosmological simulations. This result is unchanged if we consider mergers with parent mass ratios of 3:1 instead of 1:1. Mergers of a cuspy halo with either a cored halo or a second cuspy halo of equal mass, on the other hand, produce cuspy haloes with a slightly reduced inner logarithmic slope. Cuspy haloes, once formed, therefore appear resilient to major mergers. We find the velocity structure of the remnants to be mildly anisotropic, with a Maxwellian velocity distribution near the centre but not in the outer portions of the final haloes. Violent relaxation is effective only during the early phase of mergers, with phase mixing likely to be the dominant relaxation process at late times.

Since their discovery 30 years ago, pulsars have been understood to be neutron stars (NSs) born rotating rapidly (~ 10-100 ms). These neutron stars are thought to be created in supernova explosions involving massive stars, which give rise to expanding supernova remnants (SNRs). With over 220 Galactic SNRs known (Green 1998) and over 1200 radio pulsars detected (Camilo et al. 2000), it is quite surprising that few associations between the two populations have been identified with any certainty. Here we report the discovery of a remarkable 0.3 sec X-ray pulsar, PSR J1846-0258, associated with the supernova remnant Kes 75. With a characteristic age of only 723 yr, consistent with the age of Kes 75, PSR J1846-0258 is the youngest pulsar yet discovered and is being rapidly spun down by torques from a large magnetic dipole of strength ~ 5E13 G, just above the so-called quantum critical field. PSR J1846-0258 resides in this transitional regime where the magnetic field is hypothesized to separate the regular pulsars from the so-called magnetars. PSR J1846-0258 is evidently a Crab-like pulsar, however, its period, spin-down rate, spin-down conversion efficiency, are each an order-of-magnitude greater, likely the result of its extreme magnetic field.