The highly r-process enhanced (r-II) metal-poor halo stars we observe today could play a key role in understanding early ultra-faint dwarf galaxies, the smallest building blocks of the Milky Way. If a significant fraction of metal-poor r-II halo stars originated in the ultra-faint dwarf galaxies that merged to help form the Milky Way, observations of r-II stars could help us study these now-destroyed systems and probe the formation history of our Galaxy. To conduct our initial investigation into this possible connection, we use high-resolution cosmological simulations of Milky-Way-mass galaxies from the Caterpillar suite in combination with a simple, empirically motivated treatment of r-process enrichment. We determine the fraction of metal-poor halo stars that could have formed from highly r-process enhanced gas in now-destroyed low-mass ultra-faint dwarf galaxies, the simulated r-II fraction, and compare it to the "as observed" r-II fraction. We find that the simulated fraction, f_{r-II,sim} ~ 1-2%, can account for around half of the "as observed" fraction, f_{r-II,obs} ~ 2-4%. The "as observed" fraction likely overrepresents the fraction of r-II stars due to incomplete sampling, though, meaning f_{r-II,sim} likely accounts for more than half of the true f_{r-II,obs}. Further considering some parameter variations and scatter between individual simulations, the simulated fraction can account for around 20-80% of the "as observed" fraction.

Triangulum (M33) is the most massive satellite galaxy of Andromeda (M31), with a stellar mass of about $3\times10^9\; M_{\odot}$. Based on abundance matching techniques, M33's total mass at infall is estimated to be of order $10^{11}\; M_{\odot}$. $\Lambda$CDM simulations predict that M33-mass halos host several of their own low mass satellite companions, yet only one candidate M33 satellite has been detected in deep photometric surveys to date. This `satellites of satellites' hierarchy has recently been explored in the context of the dwarf galaxies discovered around the Milky Way's Magellanic Clouds in the Dark Energy Survey. Here we discuss the number of satellite galaxies predicted to reside within the virial radius ($\sim$160 kpc) of M33 based on $\Lambda$CDM simulations. We also calculate the expected number of satellite detections in N fields of data using various ground--based imagers. Finally, we discuss how these satellite population predictions change as a function of M33's past orbital history. If M33 is on its first infall into M31's halo, its proposed satellites are expected to remain bound to M33 today. However, if M33 experienced a recent tidal interaction with M31, the number of surviving satellites depends strongly on the distance achieved at pericenter due to the effects of tidal stripping. We conclude that a survey extending to $\sim$100 kpc around M33 would be sufficient to constrain its orbital history and a majority of its satellite population. In the era of WFIRST, surveys of this size will be considered small observing programs.

The smallest satellites of the Milky Way ceased forming stars during the epoch of reionization and thus provide archaeological access to galaxy formation at $z>6$. Numerical studies of these ultra-faint dwarf galaxies (UFDs) require expensive cosmological simulations with high mass resolution that are carried out down to $z=0$. However, if we are able to statistically identify UFD host progenitors at high redshifts \emph{with relatively high probabilities}, we can avoid this high computational cost. To find such candidates, we analyze the merger trees of Milky Way type halos from the high-resolution ${\it Caterpillar}$ suite of dark matter only simulations. Satellite UFD hosts at $z=0 $ are identified based on four different abundance matching (AM) techniques. All the halos at high redshifts are traced forward in time in order to compute the probability of surviving as satellite UFDs today. Our results show that selecting potential UFD progenitors based solely on their mass at z=12 (8) results in a 10\% (20\%) chance of obtaining a surviving UFD at $z=0$ in three of the AM techniques we adopted. We find that the progenitors of surviving satellite UFDs have lower virial ratios ($\eta$), and are preferentially located at large distances from the main MW progenitor, while they show no correlation with concentration parameter. Halos with favorable locations and virial ratios are $\approx 3$ times more likely to survive as satellite UFD candidates at $z=0.$

A recent surge in the discovery of new ultrafaint dwarf satellites of the Milky Way has inspired the idea of searching for faint satellites, $10^3\, \mathrm{M_{\odot}} < M_* < 10^6\, \mathrm{M_{\odot}}$, around less massive field galaxies in the Local Group. Such satellites would be subject to weaker environmental influences than Milky Way satellites, and could lead to new insights on low mass galaxy formation. In this paper, we predict the number of luminous satellites expected around field dwarf galaxies by applying several abundance matching models and a reionization model to the dark-matter only \textit{Caterpillar} simulation suite. For three of the four abundance matching models used, we find a $>99\%$ chance that at least one satellite with stellar mass $M_*> 10^5\, \mathrm{M_{\odot}}$ exists around the combined five Local Group field dwarf galaxies with the largest stellar mass. When considering satellites with $M_*> 10^4\, \mathrm{M_{\odot}}$, we predict a combined $5-25$ satellites for the five largest field dwarfs, and $10-50$ for the whole Local Group field dwarf population. Because of the relatively small number of predicted dwarfs, and their extended spatial distribution, a large fraction each Local Group dwarf's virial volume will need to be surveyed to guarantee discoveries. We compute the predicted number of satellites in a given field of view of specific Local Group galaxies, as a function of minimum satellite luminosity, and explicitly obtain such values for the Solitary Local dwarfs survey. Uncertainties in abundance matching and reionization models are large, implying that comprehensive searches could lead to refinements of both models.

Recent discovery of many dwarf satellite galaxies in the direction of the Small and Large Magellanic Clouds (SMC and LMC) provokes questions of their origins, and what they can reveal about galaxy evolution theory. Here, we predict the satellite stellar mass function of Magellanic Cloud-mass host galaxies using abundance matching and reionization models applied to the \textit{Caterpillar} simulations. Specifically focusing on the volume within $50$~kpc of the LMC, we predict a mean of 4-8 satellites with stellar mass $M_* > 10^4\, \mathrm{M_\odot}$, and 3-4 satellites with $80 < M_* \leq 3000\, \mathrm{M_\odot}$. Surprisingly, all $12$ currently known satellite candidates have stellar masses of $80 < M_* \leq 3000\, \mathrm{M_\odot}$. Reconciling the dearth of large satellites and profusion of small satellites is challenging and may require a combination of a major modification of the \mstarmhalo relationship (steep, but with an abrupt flattening at $10^3\, \mathrm{M_\odot}$), late reionization for the Local Group ($z_{\rm{reion}} \lesssim 9$ preferred), and/or strong tidal stripping. We can more robustly predict that $\sim 53\%$ of satellites within this volume were accreted together with the LMC and SMC, and $\sim 47\%$ were only ever Milky Way satellites. Observing satellites of isolated LMC-sized field galaxies is essential to placing the LMC in context, and to better constrain the \mstarmhalo relationship. Modeling known LMC-sized galaxies within $8$ Mpc, we predict 1-6 (2-12) satellites with $M_* > 10^5\, \mathrm{M_\odot}$ ($M_* > 10^4\, \mathrm{M_\odot}$) within the virial volume of each, and 1-3 (1-7) within a single $1.5^{\circ}$ diameter field of view, making their discovery likely.

We use 30 high-resolution dark matter halos of the $Caterpillar$ simulation suite to probe the first stars and galaxies of Milky Way-mass systems. We quantify the environment of the high-$z$ progenitors of the Milky Way and connect them to the properties of the host and satellites today. We identify the formation sites of the first generation of Population III (Pop III) stars ($z$ ~ 25) and first galaxies ($z$ ~ 22) with several different models based on a minimum halo mass including a simple model for Lyman-Werner feedback. Through this method we find approximately 23,000 $\pm$ 5,000 Pop III potentially star-forming sites per Milky Way-mass host, though this number is drastically reduced to ~550 star-forming sites when Lyman-Werner feedback is included, as it has critical effects at these length scales. The majority of these halos identified form in isolation (96% at $z$ = 15) and are not subject to external enrichment by neighboring halos (median separation ~1 pkpc at $z$ = 15), though half merge with a system larger than themselves within 1.5 Gyrs. Approximately 55% of the entire population has merged into the host halo by $z$ = 0. Using particle tagging, we additionally trace the Pop III remnant population to $z$ = 0 and find an order of magnitude scatter in their number density at small (i.e. r $<$ 5 kpc) and large (i.e. r $>$ 50 kpc) galactocentric radii at $z$ = 0. Using our large number of realizations, we provide fitting functions for determining the number of progenitor minihalo and atomic cooling halo systems that present-day dwarf galaxies and the Magellanic cloud system might have accreted since their formation. We demonstrate that observed dwarf galaxies with stellar masses below 10$^{4.6}$ M$_{\odot}$ are unlikely to have merged with any other star-forming systems.

We investigate the effects of self-interacting dark matter (SIDM) on the tidal stripping and evaporation of satellite galaxies in a Milky Way-like host. We use a suite of five zoom-in, dark-matter-only simulations, two with velocity-independent SIDM cross sections, two with velocity-dependent SIDM cross sections, and one cold dark matter simulation for comparison. After carefully assigning stellar mass to satellites at infall, we find that stars are stripped at a higher rate in SIDM than in CDM. In contrast, the total bound dark matter mass loss rate is minimally affected, with subhalo evaporation having negligible effects on satellites for viable SIDM models. Centrally located stars in SIDM haloes disperse out to larger radii as cores grow. Consequently, the half-light radius of satellites increases, stars become more vulnerable to tidal stripping, and the stellar mass function is suppressed. We find that the ratio of core radius to tidal radius accurately predicts the relative strength of enhanced SIDM stellar stripping. Velocity-independent SIDM models show a modest increase in the stellar stripping effect with satellite mass, whereas velocity-dependent SIDM models show a large increase in this effect towards lower masses, making observations of ultra-faint dwarfs prime targets for distinguishing between and constraining SIDM models. Due to small cores in the largest satellites of velocity-dependent SIDM, no identifiable imprint is left on the all-sky properties of the stellar halo. While our results focus on SIDM, the main physical mechanism of enhanced tidal stripping of stars apply similarly to satellites with cores formed via other means.

We present the largest number of Milky Way sized dark matter halos simulated at very high mass ($\sim$$10^4$ M$_\odot$/particle) and temporal resolution ($\sim$5 Myrs/snapshot) done to date, quadrupling what is currently available in the literature. This initial suite consists of the first 24 halos of the $Caterpillar$ $Project$ (www.caterpillarproject.org) whose project goal of 60 - 70 halos will be made public when complete. We resolve $\sim$20,000 gravitationally bound subhalos within the virial radius of each host halo. Over the ranges set by our spatial resolution our convergence is excellent and improvements were made upon current state-of-the-art halo finders to better identify substructure at such high resolutions (e.g., on average we recover $\sim$4 subhalos in each host halo above 10$^8$ M$_\odot$ which would have otherwise not been found using conventional methods). For our relaxed halos, the inner profiles are reasonably fit by Einasto profiles ($\alpha$ = 0.169 $\pm$ 0.023) though this depends on the relaxed nature and assembly history of a given halo. Averaging over all halos, the substructure mass fraction is $f_{m,subs} = 0.121 \pm 0.041$, and mass function slope is d$N$/d$M\propto M^{-1.88 \pm 0.10}$ though we find scatter in the normalizations for fixed halo mass due to more concentrated hosts having less subhalos at fixed subhalo mass. There are no biases stemming from Lagrangian volume selection as all Lagrangian volume types are included in our sample. Our detailed contamination study of 264 low resolution halos has resulted in obtaining very large and unprecedented, high-resolution regions around our host halos for our target resolution (sphere of radius $\sim$$1.4 \pm 0.4$ Mpc) allowing for accurate studies of low mass dwarf galaxies at large galactocentric radii and the very first stellar systems at high redshift ($z \geq$ 10).

We investigate how different cosmological parameters, such as those delivered by the WMAP and Planck missions, affect the nature and evolution of dark matter halo substructure. We use a series of flat $\Lambda$ cold dark matter ($\Lambda$CDM) cosmological $N$-body simulations of structure formation, each with a different power spectrum but the same initial white noise field. Our fiducial simulation is based on parameters from the WMAP 7th year cosmology. We then systematically vary the spectral index, $n_s$, matter density, $\Omega_M$, and normalization of the power spectrum, $\sigma_8$, for 7 unique simulations. Across these, we study variations in the subhalo mass function, mass fraction, maximum circular velocity function, spatial distribution, concentration, formation times, accretion times, and peak mass. We eliminate dependence of subhalo properties on host halo mass and average over many hosts to reduce variance. While the "same" subhalos from identical initial overdensity peaks in higher $\sigma_8, n_s$, and $\Omega_m$ simulations accrete earlier and end up less massive and closer to the halo center at $z=0$, the process of continuous subhalo accretion and destruction leads to a steady state distribution of these properties across all subhalos in a given host. This steady state mechanism eliminates cosmological dependence on all properties listed above except subhalo concentration and $V_{max}$, which remain greater for higher $\sigma_8, n_s$ and $\Omega_m$ simulations, and subhalo formation time, which remains earlier. We also find that the numerical technique for computing scale radius and the halo finder used can significantly affect the concentration-mass relationship computed for a simulation.