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 characte
ristics 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.
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 an
d 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 halos accretion history.
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 i
nto 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.
