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 cosmologi
cal 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.
We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The par
ticipating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data sets, stellar models, and merger trees. In this paper we apply them without recalibration and this leads to a wide variety of predictions for the stellar mass function, specific star formation rates, stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter is much larger than seen in previous comparison studies primarily because the codes have been used outside of their native environment within which they are well tested and calibrated. The purpose of the `nIFTy comparison of galaxy formation models is to bring together as many different galaxy formation modellers as possible and to investigate a common approach to model calibration. This paper provides a unified description for all participating models and presents the initial, uncalibrated comparison as a baseline for our future studies where we will develop a common calibration framework and address the extent to which that reduces the scatter in the model predictions seen here.
Merger tree codes are routinely used to follow the growth and merger of dark matter haloes in simulations of cosmic structure formation. Whereas in Srisawat et. al. we compared the trees built using a wide variety of such codes here we study the infl
uence of the underlying halo catalogue upon the resulting trees. We observe that the specifics of halo finding itself greatly influences the constructed merger trees. We find that the choices made to define the halo mass are of prime importance. For instance, amongst many potential options different finders select self-bound objects or spherical regions of defined overdensity, decide whether or not to include substructures within the mass returned and vary in their initial particle selection. The impact of these decisions is seen in tree length (the period of time a particularly halo can be traced back through the simulation), branching ratio (essentially the merger rate of subhalos) and mass evolution. We therefore conclude that the choice of the underlying halo finder is more relevant to the process of building merger trees than the tree builder itself. We also report on some built-in features of specific merger tree codes that (sometimes) help to improve the quality of the merger trees produced.
We present the results of a series of adiabatic hydrodynamical simulations of several quintessence models (both with a free and an interacting scalar field) in comparison to a standard LCDM cosmology. For each we use $2times1024^3$ particles in a $25
0$hMpc periodic box assuming WMAP7 cosmology. In this work we focus on the properties of haloes in the cosmic web at $z=0$. The web is classified into emph{voids}, emph{sheets}, emph{filaments} and emph{knots} depending on the eigenvalues of the velocity shear tensor, which are an excellent proxy for the underlying overdensity distribution. We find that the properties of objects classified according to their surrounding environment shows a substantial dependence on the underlying cosmology; for example, while $V_{rm max}$ shows average deviations of $approx5$ per cent across the different models when considering the full halo sample, comparing objects classified according to their environment, the size of the deviation can be as large as $20$ per cent. We also find that halo spin parameters are positively correlated to the coupling, whereas halo concentrations show the opposite behaviour. Furthermore, when studying the concentration-mass relation in different environments, we find that in all cosmologies underdense regions have a larger normalization and a shallower slope. While this behaviour is found to characterize all the models, differences in the best-fit relations are enhanced in (coupled) dark energy, thus providing a clearer prediction for this class of models.
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.
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.
We use a 64$h^{-1}$Mpc dark matter (DM) only cosmological simulation to examine the large scale orientation of haloes and substructures with respect the cosmic web. A web classification scheme based on the velocity shear tensor is used to assign to e
ach halo in the simulation a web type: knot, filament, sheet or void. Using $sim10^6$ haloes that span ~3 orders of magnitude in mass the orientation of the halos spin and the orbital angular momentum of subhaloes with respect to the eigenvectors of the shear tensor is examined. We find that the orbital angular momentum of subhaloes tends to align with the intermediate eigenvector of the velocity shear tensor for all haloes in knots, filaments and sheets. This result indicates that the kinematics of substructures located deep within the virialized regions of a halo is determined by its infall which in turn is determined by the large scale velocity shear, a surprising result given the virilaized nature of haloes. The non-random nature of subhalo accretion is thus imprinted on the angular momentum measured at z = 0. We also find that haloes spin axis is aligned with the third eigenvector of the velocity shear tensor in filaments and sheets: the halo spin axis points along filaments and lies in the plane of cosmic sheets.
The detection of extremely massive clusters at $z>1$ such as SPT-CL J0546-5345, SPT-CL J2106-5844, and XMMU J2235.3-2557 has been considered by some authors as a challenge to the standard LCDM$;$cosmology. In fact,assuming Gaussian initial conditions
, the theoretical expectation of detecting such objects is as low as $leq 1%$. In this textit{Letter} we discuss the probability of the existence of such objects in the light of the Vector Dark Energy (VDE) paradigm, showing by means of a series of $N$-body simulations that chances of detection are substantially enhanced in this non-standard framework.
Using a dark matter only Constrained Local UniversE Simulation (CLUES) we examine the existence of subhaloes that change their affiliation from one of the two prominent hosts in the Local Group (i.e. the Milky Way and the Andromeda galaxy) to the oth
er, and call these objects renegade subhaloes. In light of recent claims that the two Magellanic Clouds (MCs) may have originated from another region (or even the outskirts) of the Local Group or that they have been spawned by a major merger in the past of the Andromeda galaxy, we investigate the nature of such events. However, we cannot confirm that renegade subhaloes enter as deep into the potential well of their present host nor that they share the most simplest properties with the MCs, namely mass and relative velocity. Our simulation rather suggests that these renegade subhaloes appear to be flying past one host before being pulled into the other. A merger is not required to trigger such an event, it is rather the distinct environment of our simulated Local Group facilitating such behavior. Since just a small fraction of the full z=0 subhalo population are renegades, our study indicates that it will be intrinsically difficult to distinguish them despite clear differences in their velocity, radial distribution, shape and spin parameter distributions.
[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity
(SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).
