We propose a common terminology for use in describing both temporal merger trees and spatial structure trees for dark-matter halos. We specify a unified data format in HDF5 and provide example I/O routines in C, FORTRAN and PYTHON.

Merging haloes with similar masses (i.e., major mergers) pose significant challenges for halo finders. We compare five halo finding algorithms' (AHF, HBT, Rockstar, SubFind, and VELOCIraptor) recovery of halo properties for both isolated and cosmological major mergers. We find that halo positions and velocities are often robust, but mass biases exist for every technique. The algorithms also show strong disagreement in the prevalence and duration of major mergers, especially at high redshifts (z>1). This raises significant uncertainties for theoretical models that require major mergers for, e.g., galaxy morphology changes, size changes, or black hole growth, as well as for finding Bullet Cluster analogues. All finders not using temporal information also show host halo and subhalo relationship swaps over successive timesteps, requiring careful merger tree construction to avoid problematic mass accretion histories. We suggest that future algorithms should combine phase-space and temporal information to avoid the issues presented.

Investigating the spin parameter distribution of subhaloes in two high resolution isolated halo simulations, re- cent work by Onions et al. suggested that typical subhalo spins are consistently lower than the spin distribution found for field haloes. To further examine this puzzle, we have analyzed simulations of a cosmological volume with sufficient resolution to resolve a significant subhalo population. We confirm the result of Onions et al. and show that the typical spin of a subhalo decreases with decreasing mass and increasing proximity to the host halo center. We interpret this as the growing influence of tidal stripping in removing the outer layers, and hence the higher angular momentum particles, of the subhaloes as they move within the host potential. Investigating the redshift dependence of this effect, we find that the typical subhalo spin is smaller with decreasing redshift. This indicates a temporal evolution as expected in the tidal stripping scenario.

We study the shapes of subhalo distributions from four dark-matter-only simulations of Milky Way type haloes. Comparing the shapes derived from the subhalo distributions at high resolution to those of the underlying dark matter fields we find the former to be more triaxial if theanalysis is restricted to massive subhaloes. For three of the four analysed haloes the increased triaxiality of the distributions of massive subhaloes can be explained by a systematic effect caused by the low number of objects. Subhaloes of the fourth halo show indications for anisotropic accretion via their strong triaxial distribution and orbit alignment with respect to the dark matter field. These results are independent of the employed subhalo finder. Comparing the shape of the observed Milky Way satellite distribution to those of high-resolution subhalo samples from simulations, we find an agreement for samples of bright satellites, but significant deviations if faint satellites are included in the analysis. These deviations might result from observational incompleteness.

We study the dynamics of faint stellar substructures around the Umbrella Galaxy, NGC 4651, which hosts a dramatic system of streams and shells formed through the tidal disruption of a nucleated dwarf elliptical galaxy. We elucidate the basic characteristics of the system (colours, luminosities, stellar masses) using multi-band Subaru/Suprime-Cam images. The implied stellar mass-ratio of the ongoing merger event is about 1:50. We identify candidate kinematic tracers (globular clusters, planetary nebulae, H ii regions), and follow up a subset with Keck/DEIMOS spectroscopy to obtain velocities. We find that 15 of the tracers are likely associated with halo substructures, including the probable stream progenitor nucleus. These objects delineate a kinematically cold feature in position-velocity phase space. We model the stream using single test-particle orbits, plus a rescaled pre-existing N-body simulation. We infer a very eccentric orbit with a period of roughly 0.35 Gyr and turning points at approximately 2-4 and 40 kpc, implying a recent passage of the satellite through the disc, which may have provoked the visible disturbances in the host galaxy. This work confirms that the kinematics of low surface brightness substructures can be recovered and modeled using discrete tracers - a breakthrough that opens up a fresh avenue for unraveling the detailed physics of minor merging.

The chemo-dynamics of galaxy halos beyond the Local Group may now be mapped out through the use of globular clusters and planetary nebulae as bright tracer objects, along with deep multi-slit spectroscopy of the integrated stellar light. We present results from surveying nearby early-type galaxies, including evidence for kinematically distinct halos that may reflect two-phase galaxy assembly. We also demonstrate the utility of the tracer approach in measuring the kinematics of stellar substructures around the Umbrella Galaxy, which allow us to reconstruct the progenitor properties and stream orbit.

We present a study of the substructure finder dependence of subhalo clustering in the Aquarius Simulation. We run 11 different subhalo finders on the haloes of the Aquarius Simulation and we study their differences in the density profile, mass fraction and 2-point correlation function of subhaloes in haloes. We also study the mass and vmax dependence of subhalo clustering. As the Aquarius Simulation has been run at different resolutions, we study the convergence with higher resolutions. We find that the agreement between finders is at around the 10% level inside R200 and at intermediate resolutions when a mass threshold is applied, and better than 5% when vmax is restricted instead of mass. However, some discrepancies appear in the highest resolution, underlined by an observed resolution dependence of subhalo clustering. This dependence is stronger for the smallest subhaloes, which are more clustered in the highest resolution, due to the detection of subhaloes within subhaloes (the sub-subhalo term). This effect modifies the mass dependence of clustering in the highest resolutions. We discuss implications of our results for models of subhalo clustering and their relation with galaxy clustering.

Tidal streams are a powerful probe of the Milky Way (MW) potential shape. In this paper, we introduce a simple test particle method to fit stream data, using a Markov Chain Monte Carlo technique to marginalise over uncertainties in the progenitor's orbit and the Milky Way halo shape parameters. Applying it to mock data of thin streams in the MW halo, we show that, even for very cold streams, stream-orbit offsets - not modelled in our simple method - introduce systematic biases in the recovered shape parameters. For the streams that we consider, and our particular choice of potential parameterisation, these errors are of order ~20% on the halo flattening parameters. However, larger systematic errors can arise for more general streams and potentials; such offsets need to be correctly modelled in order to obtain an unbiased recovery of the underlying potential. Assessing which of the known Milky Way streams are most constraining, we find NGC 5466 and Pal 5 are the most promising candidates. These form an interesting pair as their orbital planes are both approximately perpendicular to each other and to the disc, giving optimal constraints on the MW halo shape. We show that - while with current data their constraints on potential parameters are poor - good radial velocity data along the Pal 5 stream will provide constraints on qz - the flattening perpendicular to the disc. Furthermore, as discussed in a companion paper, NGC 5466 can provide rather strong constraints on the MW halo shape parameters, if the tentative evidence for a departure from the smooth orbit towards its western edge is confirmed.

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

With the ever increasing size and complexity of fully self-consistent simulations of galaxy formation within the framework of the cosmic web, the demands upon object finders for these simulations has simultaneously grown. To this extent we initiated the Halo Finder Comparison Project that gathered together all the experts in the field and has so far led to two comparison papers, one for dark matter field haloes (Knebe et al. 2011), and one for dark matter subhaloes (Onions et al. 2012). However, as state-of-the-art simulation codes are perfectly capable of not only following the formation and evolution of dark matter but also account for baryonic physics (e.g. hydrodynamics, star formation, feedback) object finders should also be capable of taking these additional processes into consideration. Here we report on a comparison of codes as applied to the Constrained Local UniversE Simulation (CLUES) of the formation of the Local Group which incorporates much of the physics relevant for galaxy formation. We compare both the properties of the three main galaxies in the simulation (representing the MW, M31, and M33) as well as their satellite populations for a variety of halo finders ranging from phase-space to velocity-space to spherical overdensity based codes, including also a mere baryonic object finder. We obtain agreement amongst codes comparable to (if not better than) our previous comparisons, at least for the total, dark, and stellar components of the objects. However, the diffuse gas content of the haloes shows great disparity, especially for low-mass satellite galaxies. This is primarily due to differences in the treatment of the thermal energy during the unbinding procedure. We acknowledge that the handling of gas in halo finders is something that needs to be dealt with carefully, and the precise treatment may depend sensitively upon the scientific problem being studied.

The spatial segregation between dSphs and dIrrs in the Local Group has long been regarded as evidence of an interaction with their host galaxies. In this paper, we assume that ram-pressure stripping is the dominant mechanism that removed gas from the dSphs and we use this to derive a lower bound on the density of the corona of the Milky Way at large distances (50-90 kpc) from the Galactic centre. At the same time, we derive an upper bound by demanding that the interstellar medium of the dSphs is in pressure equilibrium with the hot corona. We consider two dwarfs (Sextans and Carina) with well-determined orbits and star formation histories. Our approach introduces several novel features: we use the measured star formation histories of the dwarfs to derive the time at which they last lost their gas, and (via a modified version of the Kennicutt-Schmidt relation) their internal gas density at that time; we use a large suite of 2D hydrodynamical simulations to model the gas stripping; and we include supernova feedback tied to the gas content. Despite having very different orbits and star formation histories, we find results for the two dSphs that are in excellent agreement with one another. We derive an average particle density of the corona of the Milky Way at 50-90 kpc in the range 1.3-3.6 10^{-4} cm^{-3}. Including additional constraints from X-ray emission limits and pulsar dispersion measurements, we extrapolate Galactic coronal density profiles and we estimate the fraction of baryons that can exist within the virial radius of the Milky Way. For an isothermal corona (T=1.8 10^6 K) this is small, 10-20 % of the universal baryon fraction. Only a hot (T=3 10^6 K) and adiabatic corona can contain all of the Galaxy's missing baryons. Models for the Milky Way must explain why its corona is in a hot adiabatic thermal state or why a large fraction of its baryons lie beyond the virial radius.

While various codes exist to systematically and robustly find haloes and subhaloes in cosmological simulations (Knebe et al., 2011, Onions et al., 2012), this is the first work to introduce and rigorously test codes that find tidal debris (streams and other unbound substructure) in fully cosmological simulations of structure formation. We use one tracking and three non-tracking codes to identify substructure (bound and unbound) in a Milky Way type simulation from the Aquarius suite (Springel et al., 2008) and post-process their output with a common pipeline to determine the properties of these substructures in a uniform way. By using output from a fully cosmological simulation, we also take a step beyond previous studies of tidal debris that have used simple toy models. We find that both tracking and non-tracking codes agree well on the identification of subhaloes and more importantly, the {\em unbound tidal features} associated with them. The distributions of basic properties of the total substructure distribution (mass, velocity dispersion, position) are recovered with a scatter of $\sim20%$. Using the tracking code as our reference, we show that the non-tracking codes identify complex tidal debris with purities of $\sim40%$. Analysing the results of the substructure finders, we find that the general distribution of {\em substructures} differ significantly from the distribution of bound {\em subhaloes}. Most importantly, both bound and unbound {\em substructures} together constitute $\sim18%$ of the host halo mass, which is a factor of $\sim2$ higher than the fraction in self-bound {\em subhaloes}. However, this result is restricted by the remaining challenge to cleanly define when an unbound structure has become part of the host halo. Nevertheless, the more general substructure distribution provides a more complete picture of a halo's accretion history.

We present a study of a comparison of spin distributions of subhaloes found associated with a host halo. The subhaloes are found within two cosmological simulation families of Milky Way-like galaxies, namely the Aquarius and GHALO simulations. These two simulations use different gravity codes and cosmologies. We employ ten different substructure finders, which span a wide range of methodologies from simple overdensity in configuration space to full 6-d phase space analysis of particles.We subject the results to a common post-processing pipeline to analyse the results in a consistent manner, recovering the dimensionless spin parameter. We find that spin distribution is an excellent indicator of how well the removal of background particles (unbinding) has been carried out. We also find that the spin distribution decreases for substructure the nearer they are to the host halo's, and that the value of the spin parameter rises with enclosed mass towards the edge of the substructure. Finally subhaloes are less rotationally supported than field haloes, with the peak of the spin distribution having a lower spin parameter.

Stellar streams provide unique probes of galactic potentials, with the longer streams normally providing the cleaner measurements. In this paper, we show an example of a short tidal stream that is particularly sensitive to the shape of the Milky Way's dark matter halo: the globular cluster tidal stream NGC 5466. This stream has an interesting deviation from a smooth orbit at its western edge. We show that such a deviation favours an underlying oblate or triaxial halo (irrespective of plausible variations in the Milky Way disc properties and the specific halo parametrisation chosen); spherical or prolate halo shapes can be excluded at a high confidence level. Therefore, more extensive data sets along the NGC 5466 tidal stream promise strong constraints on the Milky Way halo shape.

We present a detailed comparison of the substructure properties of a single Milky Way sized dark matter halo from the Aquarius suite at five different resolutions, as identified by a variety of different (sub-)halo finders for simulations of cosmic structure formation. These finders span a wide range of techniques and methodologies to extract and quantify substructures within a larger non-homogeneous background density (e.g. a host halo). This includes real-space, phase-space, velocity-space and time- space based finders, as well as finders employing a Voronoi tessellation, friends-of-friends techniques, or refined meshes as the starting point for locating substructure.A common post-processing pipeline was used to uniformly analyse the particle lists provided by each finder. We extract quantitative and comparable measures for the subhaloes, primarily focusing on mass and the peak of the rotation curve for this particular study. We find that all of the finders agree extremely well on the presence and location of substructure and even for properties relating to the inner part part of the subhalo (e.g. the maximum value of the rotation curve). For properties that rely on particles near the outer edge of the subhalo the agreement is at around the 20 per cent level. We find that basic properties (mass, maximum circular velocity) of a subhalo can be reliably recovered if the subhalo contains more than 100 particles although its presence can be reliably inferred for a lower particle number limit of 20. We finally note that the logarithmic slope of the subhalo cumulative number count is remarkably consistent and <1 for all the finders that reached high resolution. If correct, this would indicate that the larger and more massive, respectively, substructures are the most dynamically interesting and that higher levels of the (sub-)subhalo hierarchy become progressively less important.

We calculate orbits for the Milky Way dwarf galaxies with proper motions, and compare these to subhalo orbits in a high resolution cosmological simulation. We use the simulation data to assess how well orbits may be recovered in the face of measurement errors, a time varying triaxial gravitational potential, and satellite-satellite interactions. For present measurement uncertainties, we recover the apocentre r_a and pericentre r_p to ~40%. With improved data from the Gaia satellite we should be able to recover r_a and r_p to ~14%, respectively. However, recovering the 3D positions and orbital phase of satellites over several orbits is more challenging. This owes primarily to the non-sphericity of the potential and satellite interactions during group infall. Dynamical friction, satellite mass loss and the mass evolution of the main halo play a more minor role in the uncertainties. We apply our technique to nine Milky Way dwarfs with observed proper motions. We show that their mean apocentre is lower than the mean of the most massive subhalos in our cosmological simulation, but consistent with the most massive subhalos that form before z=10. This lends further support to the idea that the Milky Way's dwarfs formed before reionisation.

We calculate orbits for the Milky Way dwarf galaxies with proper motions, and compare these to subhalo orbits in a high resolution cosmological simulation. We use this same simulation to assess how well are able to recover orbits in the face of measurement errors, a time varying triaxial gravitational potential, and satellite-satellite interactions. We find that, for present measurement uncertainties, we are able to recover the apocentre r_a and pericentre r_p to ~ 40%. However, even with better data the non-sphericity of the potential and satellite interactions during group infall make the orbital recovery more challenging. Dynamical friction, satellite mass loss and the mass evolution of the main halo play a more minor role. We apply our technique to nine Milky Way dwarfs with observed proper motions. We show that their mean apocentre is consistent with the most massive subhalos that form before z=10, lending support to the idea that the Milky Way dwarfs formed before reionisation.