Extremely metal-poor (EMP) stars in the Milky Way (MW) allow us to infer the properties of their progenitors by comparing their chemical composition to the metal yields of the first supernovae. This method is most powerful when applied to mono-enriched stars, i.e. stars that formed from gas that was enriched by only one previous supernova. We present a novel diagnostic to identify this subclass of EMP stars. We model the first generations of star formation semi-analytically, based on dark matter halo merger trees that yield MW-like halos at the present day. Radiative and chemical feedback are included self-consistently and we trace all elements up to zinc. Mono-enriched stars account for only $\sim 1\%$ of second generation stars in our fiducial model and we provide an analytical formula for this probability. We also present a novel analytical diagnostic to identify mono-enriched stars, based on the metal yields of the first supernovae. This new diagnostic allows us to derive our main results independently from the specific assumptions made regarding Pop III star formation, and we apply it to a set of observed EMP stars to demonstrate its strengths and limitations. Our results may provide selection criteria for current and future surveys and therefore contribute to a deeper understanding of EMP stars and their progenitors.

The search for metal-free stars has so far been unsuccessful, proving that if there are surviving stars from the first generation, they are rare, they have been polluted, or we have been looking in the wrong place. To predict the likely location of Population~III (Pop~III) survivors, we semi-analytically model early star formation in progenitors of Milky Way-like galaxies and their environments. We base our model on merger trees from the high-resolution dark matter only simulation suite \textit{Caterpillar}. Radiative and chemical feedback are taken into account self-consistently, based on the spatial distribution of the haloes. Our results are consistent with the non-detection of Pop III survivors in the Milky Way today. We find that possible surviving Population III stars are more common in Milky Way satellites than in the main Galaxy. In particular, low mass Milky Way satellites contain a much larger fraction of Pop~III stars than the Milky Way. Such nearby, low mass Milky Way satellites are promising targets for future attempts to find Pop~III survivors, especially for high-resolution, high signal-to-noise spectroscopic observations. We provide the probabilities for finding a Pop~III survivor in the red giant branch phase for all known Milky Way satellites to guide future observations.

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.

Mergers and the spin of the dark matter halo are factors traditionally believed to determine the morphology of galaxies within a $\Lambda$CDM cosmology. We study this hypothesis by considering approximately 18,000 central galaxies at $z=0$ with stellar masses $M_{\ast} = 10^{9}-10^{12} \, {\rm M}_{\odot}$ selected from the Illustris cosmological hydrodynamic simulation. The fraction of accreted stars -- which measures the importance of massive, recent and dry mergers -- increases steeply with galaxy stellar mass, from less than 5 per cent in dwarfs to 80 per cent in the most massive objects, and the impact of mergers on galaxy morphology increases accordingly. For galaxies with $M_{\ast} \gtrsim 10^{11} \, {\rm M}_{\odot}$, mergers have the expected effect: if gas-poor they promote the formation of spheroidal galaxies, whereas gas-rich mergers favour the formation and survivability of massive discs. This trend, however, breaks at lower masses. For objects with $M_{\ast} \lesssim 10^{11} \, {\rm M}_{\odot}$, mergers do not seem to play any significant role in determining the morphology, with accreted stellar fractions and mean merger gas fractions that are indistinguishable between spheroidal and disc-dominated galaxies. On the other hand, halo spin correlates with morphology primarily in the least massive objects in the sample ($M_{\ast} \lesssim 10^{10} \, {\rm M}_{\odot}$), but only weakly for galaxies above that mass. Our results support a scenario where (1) mergers play a dominant role in shaping the morphology of massive galaxies, (2) halo spin is important for the morphology of dwarfs, and (3) the morphology of medium-sized galaxies -- including the Milky Way -- shows little dependence on galaxy assembly history or halo spin, at least when these two factors are considered individually.

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 compare the predictions of a semi-analytic model for ultra-compact dwarf galaxy (UCD) formation by tidal stripping to the observed properties of globular clusters (GCs) and UCDs in the Fornax and Virgo clusters. For Fornax we find the predicted number of stripped nuclei agrees very well with the excess number of GCs$+$UCDs above the GC luminosity function. GCs$+$UCDs with masses $>10^{7.3}$ M$_\odot$ are consistent with being entirely formed by tidal stripping. Stripped nuclei can also account for Virgo UCDs with masses $>10^{7.3}$ M$_\odot$ where numbers are complete by mass. For both Fornax and Virgo, the predicted velocity dispersions and radial distributions of stripped nuclei are consistent with that of UCDs within $\sim$50-100 kpc but disagree at larger distances where dispersions are too high and radial distributions too extended. Stripped nuclei are predicted to have radially biased anisotropies at all radii, agreeing with Virgo UCDs at clustercentric distances larger than 50 kpc. However, ongoing disruption is not included in our model which would cause orbits to become tangentially biased at small radii. We find the predicted metallicities and central black hole masses of stripped nuclei agree well with the metallicities and implied black hole masses of UCDs for masses $>10^{6.5}$ M$_\odot$. The predicted black hole masses also agree well with that of M60-UCD1, the first UCD with a confirmed central black hole. These results suggest that observed GC$+$UCD populations are a combination of genuine GCs and stripped nuclei, with the contribution of stripped nuclei increasing toward the high-mass end.

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 present the full public release of all data from the Illustris simulation project. Illustris is a suite of large volume, cosmological hydrodynamical simulations run with the moving-mesh code Arepo and including a comprehensive set of physical models critical for following the formation and evolution of galaxies across cosmic time. Each simulates a volume of (106.5 Mpc)^3 and self-consistently evolves five different types of resolution elements from a starting redshift of z=127 to the present day, z=0. These components are: dark matter particles, gas cells, passive gas tracers, stars and stellar wind particles, and supermassive black holes. This data release includes the snapshots at all 136 available redshifts, halo and subhalo catalogs at each snapshot, and two distinct merger trees. Six primary realizations of the Illustris volume are released, including the flagship Illustris-1 run. These include three resolution levels with the fiducial "full" baryonic physics model, and a dark matter only analog for each. In addition, we provide four distinct, high time resolution, smaller volume "subboxes". The total data volume is ~265 TB, including ~800 full volume snapshots and ~30,000 subbox snapshots. We describe the released data products as well as tools we have developed for their analysis. All data may be directly downloaded in its native HDF5 format. Additionally, we release a comprehensive, web-based API which allows programmatic access to search and data processing tasks. In both cases we provide example scripts and a getting-started guide in several languages: currently, IDL, Python, and Matlab. This paper addresses scientific issues relevant for the interpretation of the simulations, serves as a pointer to published and on-line documentation of the project, describes planned future additional data releases, and discusses technical aspects of the release.

The cumulative comoving number-density of galaxies as a function of stellar mass or central velocity dispersion is commonly used to link galaxy populations across different epochs. By assuming that galaxies preserve their number-density in time, one can infer the evolution of their properties, such as masses, sizes, and morphologies. However, this assumption does not hold in the presence of galaxy mergers or when rank ordering is broken owing to variable stellar growth rates. We present an analysis of the evolving comoving number density of galaxy populations found in the Illustris cosmological hydrodynamical simulation focused on the redshift range $0\leq z \leq 3$. Our primary results are as follows: 1) The inferred average stellar mass evolution obtained via a constant comoving number density assumption is systematically biased compared to the merger tree results at the factor of $\sim$2(4) level when tracking galaxies from redshift $z=0$ out to redshift $z=2(3)$; 2) The median number density evolution for galaxy populations tracked forward in time is shallower than for galaxy populations tracked backward in time; 3) A similar evolution in the median number density of tracked galaxy populations is found regardless of whether number density is assigned via stellar mass, stellar velocity dispersion, or dark matter halo mass; 4) Explicit tracking reveals a large diversity in galaxies' assembly histories that cannot be captured by constant number-density analyses; 5) The significant scatter in galaxy linking methods is only marginally reduced by considering a number of additional physical and observable galaxy properties as realized in our simulation. We provide fits for the forward and backward median evolution in stellar mass and number density and discuss implications of our analysis for interpreting multi-epoch galaxy property observations.

We use the Millennium II cosmological simulation combined with the semi-analytic galaxy formation model of Guo et al. (2011) to predict the contribution of galactic nuclei formed by the tidal stripping of nucleated dwarf galaxies to globular cluster (GC) and ultra-compact dwarf galaxy (UCD) populations of galaxies. We follow the merger trees of galaxies in clusters back in time and determine the absolute number and stellar masses of disrupted galaxies. We assume that at all times nuclei have a distribution in nucleus-to-galaxy mass and nucleation fraction of galaxies similar to that observed in the present day universe. Our results show stripped nuclei follow a mass function $N(M) \sim M^{-1.5}$ in the mass range $10^6 < M/M_\odot < 10^8$, significantly flatter than found for globular clusters. The contribution of stripped nuclei will therefore be most important among high-mass GCs and UCDs. For the Milky Way we predict between 1 and 3 star clusters more massive than $10^5 M_\odot$ come from tidally disrupted dwarf galaxies, with the most massive cluster formed having a typical mass of a few times $10^6 M_\odot$, like omega Centauri. For a galaxy cluster with a mass $7 \times 10^{13} M_\odot$, similar to Fornax, we predict $\sim$19 UCDs more massive than $2\times10^6 M_\odot$ and $\sim$9 UCDs more massive than $10^7 M_\odot$ within a projected distance of 300 kpc come from tidally stripped dwarf galaxies. The observed number of UCDs are $\sim$200 and 23, respectively. We conclude that most UCDs in galaxy clusters are probably simply the high mass end of the GC mass function.

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.

We present detailed radiative transfer simulations of the reionization history of the Milky Way by metal-poor globular clusters. We identify potential metal-poor globular cluster candidates within the Aquarius simulation using dark matter halo velocity dispersions. We calculate the local ionization fields via a photon-conserving, three dimensional non-equilibrium chemistry code and allow the model to propagate through to the present day. The key feature of the model is that globular cluster formation is suppressed if the local gas is ionized. We find that our spatial treatment of the ionization field leads to drastically different numbers and spatial distributions when compared to models where globular cluster formation is simply truncated at a given redshift. We find that it is possible for metal-poor globular clusters to have formed via the dark matter halo formation channel as our secondary model (delayed formation) combined with truncation at z = 10 produces radial distributions statistically consistent with that of the Milky Way metal-poor globular clusters. If globular clusters do indeed form within high-redshift dark matter halos, if only in-part, their contributions to the reionization of the local (i.e. 2^3 h^-3 Mpc^3 centred on the host galaxy) volume and mass by redshift 10 could be as high as 98% and 90%, respectively. In our photon poorest model, this contribution drops to 60% and 50%. The surviving clusters in all models have a narrow average age range (mean = 13.34 Gyr, \sigma = 0.04 Gyr) consistent with current ages estimates of the Milky Way metal-poor globular clusters. We also test a simple dynamical destruction model and estimate that ~60% of all metal-poor globular clusters formed at high redshift have since been destroyed via tidal interactions with the host galaxy.

The Aquarius project is the first simulation that can resolve the full mass range of potential globular cluster formation sites. With a particle mass $m_\mathrm{p}=1.4 \times 10^4$\Msun, Aquarius yields more than 100 million particles within the virial radius of the central halo which has a mass of $1.8 \times 10^{12}$\Msun, similar to that of the Milky Way. With this particle mass, dark matter concentrations (haloes) as small as 10$^6$ M$_\odot$ will contain a minimum of 100 particles.Here, we use this simulation to test a model of metal-poor globular cluster formation based on collapse physics. In our model, globular clusters form when the virial temperatures of haloes first exceed $10^4$ K as this is when electronic transitions allow the gas to cool efficiently. We calculate the ionising flux from the stars in these first clusters and stop the formation of new clusters when all the baryonic gas of the galaxy is ionised. This is achieved by adopting reasonable values for the star formation efficiencies and escape fraction of ionising photons. The model is successful in that it predicts ages (peak age $\sim$ 13.3 Gyrs) and a spatial distribution of metal-poor globular clusters which are consistent with the observed populations in the Milky Way. We then test a scenario of metal-rich cluster formation by examining mergers that trigger star formation within central gas disks. This results in younger ($\sim$ 7--13.3 Gyrs), more centrally-located clusters (40 metal-rich GCs within 18 kpc from the centre of the host) which are consistent with the Galactic metal-rich population. We test an alternate model in which metal-rich globular clusters form in dwarf galaxies that become stripped as they merge with the main halo. We find that this process is inconsistent with observed metal-rich globulars.