We have observed the environments of a population of 33 heavily dust obscured, ultra-luminous, high-redshift galaxies, selected using WISE and NVSS at $z>$1.3 with the Infra-Red Array Camera on the $Spitzer$ Space Telescope over $\rm5.12\,'\times5.12\,'$ fields. Colour selections are used to quantify any potential overdensities of companion galaxies in these fields. We find no significant excess of galaxies with the standard colour selection for IRAC colours of $\rm[3.6]-[4.5]>-0.1$ consistent with galaxies at $z>$1.3 across the whole fields with respect to wide-area $Spitzer$ comparison fields, but there is a $\rm>2\sigma$ statistical excess within $\rm0.25\,'$ of the central radio-WISE galaxy. Using a colour selection of $\rm[3.6]-[4.5]>0.4$, 0.5 magnitudes redder than the standard method of selecting galaxies at $z>$1.3, we find a significant overdensity, in which $\rm76\%$ ($\rm33\%$) of the 33 fields have a surface density greater than the $\rm3\sigma$ ($\rm5\sigma$) level. There is a statistical excess of these redder galaxies within $\rm0.5\,'$, rising to a central peak $\rm\sim2$--4 times the average density. This implies that these galaxies are statistically linked to the radio-WISE selected galaxy, indicating similar structures to those traced by red galaxies around radio-loud AGN.

We investigate galaxy evolution in protoclusters using a semi-analytic model applied to the Millennium Simulation, scaled to a Planck cosmology. We show that the model reproduces the observed behaviour of the star formation history (SFH) both in protoclusters and the field. The rate of star formation peaks $\sim0.7\,{\rm Gyr}$ earlier in protoclusters than in the field and declines more rapidly afterwards. This results in protocluster galaxies forming significantly earlier: 80% of their stellar mass is already formed by $z=1.4$, but only 45% of the field stellar mass has formed by this time. The model predicts that field and protocluster galaxies have similar average specific star-formation rates (sSFR) at $z>3$, and we find evidence of an enhancement of star formation in the dense protoclusters at early times. At $z<3$, protoclusters have lower sSFRs, resulting in the disparity between the SFHs. We show that the stellar mass functions of protoclusters are top-heavy compared with the field due to the early formation of massive galaxies, and the disruption and merging of low-mass satellite galaxies in the main haloes. The fundamental cause of the different SFHs and mass functions is that dark matter haloes are biased tracers of the dark matter density field: the high density of haloes and the top-heavy halo mass function in protoclusters result in the early formation then rapid merging and quenching of galaxies. We compare our results with observations from the literature, and highlight which observables provide the most informative tests of galaxy formation.

We investigate the prevalence of AGN in the high-redshift protocluster $\rm{Cl}\,0218.3$-$0510$ at $z=1.62$. Using imaging from the Chandra X-ray Telescope, we find a large overdensity of AGN in the protocluster; a factor of $23\pm9$ times the field density of AGN. Only half of this AGN overdensity is due to the overdensity of massive galaxies in the protocluster (a factor of $11\pm2$), as we find that $17^{+6}_{-5}\%$ of massive galaxies ($M_* > 10^{10}\,\rm{M}_{\odot}$) in the protocluster host an X-ray luminous AGN, compared to $8\pm1\%$ in the field. This corresponds to an enhancement of AGN activity in massive protocluster galaxies by a factor of $2.1\pm0.7$ at $1.6\sigma$ significance. We also find that the AGN overdensity is centrally concentrated, located within 3 arcmin and most pronounced within 1 arcmin of the centre of the protocluster. Our results confirm that there is a reversal in the local anti-correlation between galaxy density and AGN activity, so there is an enhancement of AGN in high-redshift protoclusters. We compare the properties of AGN in the protocluster to the field and find no significant differences in the distributions of their stellar mass, X-ray luminosity, or hardness ratio. We therefore suggest that triggering mechanisms are similar in both environments, and that the mechanisms simply occur more frequently in denser environments.

We investigate the effects of dense environments on galaxy evolution by examining how the properties of galaxies in the z = 1.6 protocluster Cl 0218.3-0510 depend on their location. We determine galaxy properties using spectral energy distribution fitting to 14-band photometry, including data at three wavelengths that tightly bracket the Balmer and 4000A breaks of the protocluster galaxies. We find that two-thirds of the protocluster galaxies, which lie between several compact groups, are indistinguishable from field galaxies. The other third, which reside within the groups, differ significantly from the intergroup galaxies in both colour and specific star formation rate. We find that the fraction of red galaxies within the massive protocluster groups is twice that of the intergroup region. These excess red galaxies are due to enhanced fractions of both passive galaxies (1.7 times that of the intergroup region) and dusty star-forming galaxies (3 times that of the intergroup region). We infer that some protocluster galaxies are processed in the groups before the cluster collapses. These processes act to suppress star formation and change the mode of star formation from unobscured to obscured.

We present a comprehensive picture of the Cl0218.3-0510 protocluster at $z=1.623$ across 10 co-moving Mpc. Using filters that tightly bracket the Balmer and 4000 Angstrom breaks of the protocluster galaxies we obtain precise photometric redshifts resulting in a protocluster galaxy sample that is 89+/-5% complete and has a contamination of only 12+/-5%. Both star forming and quiescent protocluster galaxies are located allowing us to map the structure of the forming cluster for the first time. The protocluster contains 6 galaxy groups, the largest of which is the nascent cluster. Only a small minority of the protocluster galaxies are in the nascent cluster (11%) or in the other galaxy groups (22%), as most protocluster galaxies reside between the groups. Unobscured star forming galaxies predominantly reside between the protocluster's groups, whereas red galaxies make up a large fraction of the groups' galactic content, so observing the protocluster through only one of these types of galaxies results in a biased view of the protocluster's structure. The structure of the protocluster reveals how much mass is available for the future growth of the cluster and we use the Millennium Simulation, scaled to a Planck cosmology, to predict that Cl0218.3-0510 will evolve into a 2.7x 10^14 Msun cluster by the present day.

We show that black holes supplied with mass at hyper--Eddington rates drive outflows with mildly sub--relativistic velocities. These are $\sim 0.1 - 0.2c$ for Eddington accretion factors $\dot m_{\rm acc} \sim 10 - 100$, and $\sim 1500\,{\rm km\, s^{-1}}$ for $\dot m_{\rm acc} \sim 10^4$. Winds like this are seen in the X--ray spectra of ultraluminous sources (ULXs), strongly supporting the view that ULXs are stellar--mass compact binaries in hyper--Eddington accretion states. SS433 appears to be an extreme ULX system ($\dot m_{\rm acc} \sim 10^4$) viewed from outside the main X--ray emission cone. For less extreme Eddington factors $\dot m_{\rm acc} \sim 10 - 100$ the photospheric temperatures of the winds are $\sim 100$\, eV, consistent with the picture that the ultraluminous supersoft sources (ULSs) are ULXs seen outside the medium--energy X--ray beam, unifying the ULX/ULS populations and SS433 (actually a ULS but with photospheric emission too soft to detect). For supermassive black holes (SMBHs), feedback from hyper--Eddington accretion is significantly more powerful than the usual near--Eddington (`UFO') case, and if realised in nature would imply $M - \sigma$ masses noticeably smaller than observed. We suggest that the likely warping of the accretion disc in such cases may lead to much of the disc mass being expelled, severely reducing the incidence of such strong feedback. We show that hyper--Eddington feedback from bright ULXs can have major effects on their host galaxies. This is likely to have important consequences for the formation and survival of small galaxies.

We use a sample of 37 of the densest clusters and protoclusters across $1.3 \le z \le 3.2$ from the Clusters Around Radio-Loud AGN (CARLA) survey to study the formation of massive cluster galaxies. We use optical $i'$-band and infrared 3.6$\mu$m and 4.5$\mu$m images to statistically select sources within these protoclusters and measure their median observed colours; $\langle i'-[3.6] \rangle$. We find the abundance of massive galaxies within the protoclusters increases with decreasing redshift, suggesting these objects may form an evolutionary sequence, with the lower redshift clusters in the sample having similar properties to the descendants of the high redshift protoclusters. We find that the protocluster galaxies have an approximately unevolving observed-frame $i'-[3.6]$ colour across the examined redshift range. We compare the evolution of the $\langle i'-[3.6] \rangle$ colour of massive cluster galaxies with simplistic galaxy formation models. Taking the full cluster population into account, we show that the formation of stars within the majority of massive cluster galaxies occurs over at least 2Gyr, and peaks at $z \sim 2$-3. From the median $i'-[3.6]$ colours we cannot determine the star formation histories of individual galaxies, but their star formation must have been rapidly terminated to produce the observed red colours. Finally, we show that massive galaxies at $z>2$ must have assembled within 0.5Gyr of them forming a significant fraction of their stars. This means that few massive galaxies in $z>2$ protoclusters could have formed via dry mergers.

We explore the structures of protoclusters and their relationship with high redshift clusters using the Millennium Simulation combined with a semi-analytic model. We find that protoclusters are very extended, with 90 per cent of their mass spread across $\sim35\,h^{-1}{\rm Mpc}$ comoving at $z=2$ ($\sim30\, \rm{arcmin}$). The `main halo', which can manifest as a high redshift cluster or group, is only a minor feature of the protocluster, containing less than 20 per cent of all protocluster galaxies at $z=2$. Furthermore, many protoclusters do not contain a main halo that is massive enough to be identified as a high redshift cluster. Protoclusters exist in a range of evolutionary states at high redshift, independent of the mass they will evolve to at $z=0$. We show that the evolutionary state of a protocluster can be approximated by the mass ratio of the first and second most massive haloes within the protocluster, and the $z=0$ mass of a protocluster can be estimated to within 0.2 dex accuracy if both the mass of the main halo and the evolutionary state is known. We also investigate the biases introduced by only observing star-forming protocluster members within small fields. The star formation rate required for line-emitting galaxies to be detected is typically high, which leads to the artificial loss of low mass galaxies from the protocluster sample. This effect is stronger for observations of the centre of the protocluster, where the quenched galaxy fraction is higher. This loss of low mass galaxies, relative to the field, distorts the size of the galaxy overdensity, which in turn can contribute to errors in predicting the $z=0$ evolved mass.

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.

This article is the second in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our aim is to quantify the scatter, systematic bias and completeness of cluster masses derived from a diverse set of 25 galaxy-based methods using two contrasting mock galaxy catalogues based on a sophisticated halo occupation model and a semi-analytic model. Analysing 968 clusters, we find a wide range in the RMS errors in log M200c delivered by the different methods (0.18 to 1.08 dex, i.e., a factor of ~1.5 to 12), with abundance matching and richness methods providing the best results, irrespective of the input model assumptions. In addition, certain methods produce a significant number of catastrophic cases where the mass is under- or over-estimated by a factor greater than 10. Given the steeply falling high-mass end of the cluster mass function, we recommend that richness or abundance matching-based methods are used in conjunction with these methods as a sanity check for studies selecting high mass clusters. We see a stronger correlation of the recovered to input number of galaxies for both catalogues in comparison with the group/cluster mass, however, this does not guarantee that the correct member galaxies are being selected. We do not observe significantly higher scatter for either mock galaxy catalogues. Our results have implications for cosmological analyses that utilise the masses, richnesses, or abundances of clusters, which have different uncertainties when different methods are used.

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.

This paper is the first in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our primary aim is to test the performance of these cluster mass estimation techniques on a diverse set of models that will increase in complexity. We begin by providing participating methods with data from a simple model that delivers idealised clusters, enabling us to quantify the underlying scatter intrinsic to these mass estimation techniques. The mock catalogue is based on a Halo Occupation Distribution (HOD) model that assumes spherical Navarro, Frenk and White (NFW) haloes truncated at R_200, with no substructure nor colour segregation, and with isotropic, isothermal Maxwellian velocities. We find that, above 10^14 M_solar, recovered cluster masses are correlated with the true underlying cluster mass with an intrinsic scatter of typically a factor of two. Below 10^14 M_solar, the scatter rises as the number of member galaxies drops and rapidly approaches an order of magnitude. We find that richness-based methods deliver the lowest scatter, but it is not clear whether such accuracy may simply be the result of using an over-simplistic model to populate the galaxies in their haloes. Even when given the true cluster membership, large scatter is observed for the majority non-richness-based approaches, suggesting that mass reconstruction with a low number of dynamical tracers is inherently problematic.

We present results from a narrow-band survey of the field around the high redshift radio galaxy MRC 2104-242. We have selected Halpha emitters in a 7sq.arcmin field and compared the measured number density with that of a field sample at similar redshift. We find that MRC 2104-242 lies in an overdensity of galaxies that is 8.0 +/- 0.8 times the average density of a blank field, suggesting it resides in a large-scale structure that may eventually collapse to form a massive cluster. We find that there is more dust obscured star formation in the protocluster galaxies than in similarly selected control field galaxies and there is tentative evidence of a higher fraction of starbursting galaxies in the denser environment. However, on average we do not find a difference between the star formation rate (SFR)-mass relations of the protocluster and field galaxies and so conclude that the SFR of these galaxies at z~2.5 is governed predominantly by galaxy mass and not the host environment. We also find that the stellar mass distribution of the protocluster galaxies is skewed towards higher masses and there is a significant lack of galaxies at M < 10^10Msun within our small field of view. Based on the level of overdensity we expect to find ~22 star forming galaxies below 10^10Msun in the protocluster and do not detect any. This lack of low mass galaxies affects the level of overdensity which we detect. If we only consider high mass (M > 10^10.5Msun) galaxies, the density of the protocluster field increases to ~55 times the control field density.

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.

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.

There is strong evidence that supermassive black holes reside in all galaxies that contain a stellar spheroid and their mass is tightly correlated with properties such as stellar bulge mass and velocity dispersion. There are also strong theoretical arguments that feedback from supermassive black holes plays an important role in shaping the high mass end of the galaxy mass function, hence to accurately model galaxies we also need to model the black holes. We present a comparison of two black hole growth models implemented within a large-scale, cosmological SPH simulation including star formation and feedback. One model is a modified Bondi-Hoyle prescription that grows black holes based on the smooth density of local gas, while the other is the recently proposed Accretion Disc Particle (ADP) method. This model swallows baryonic particles that pass within an accretion radius of the black hole and adds them to a subgrid accretion disc. Black holes are then grown by material from this disc. We find that both models can reproduce local scaling relations, although the ADP model is offset from the observed relations at high black hole masses. The total black hole mass density agrees between models to within a factor of three, but both struggle to reproduce the black hole mass function. The simulated mass functions are too steep and underestimate the number of intermediate and high mass black holes. In addition, the ADP model swallows excessive amounts of material at the resolution of large-scale, cosmological simulations producing unrealistically large accretion discs. Future work needs to be performed to improve the black hole mass function within simulations. This should be done through the mass growth and feedback as they are strongly coupled and should not be treated as separate entities.

Galaxy environment is frequently discussed, but inconsistently defined. It is especially difficult to measure at high redshift where only photometric redshifts are available. With a focus on early forming proto-clusters, we use a semi-analytical model of galaxy formation to show how the environment measurement around high redshift galaxies is sensitive to both scale and metric, as well as to cluster viewing angle, evolutionary state, and the availability of either spectroscopic or photometric data. We use two types of environment metrics (nearest neighbour and fixed aperture) at a range of scales on simulated high-z clusters to see how "observed" overdensities compare to "real" overdensities. We also "observationally" identify z = 2 proto-cluster candidates in our model and track the growth histories of their parent halos through time, considering in particular their final state at z = 0. Although the measured environment of early forming clusters is critically dependent on all of the above effects (and in particular the viewing angle), we show that such clusters are very likely (< 90%) to remain overdense at z = 0, although many will no longer be among the most massive. Object to object comparisons using different methodologies and different data, however, require much more caution.

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 analyze environmental correlations using mark clustering statistics with the mock galaxy catalogue constructed by Muldrew et al. (Paper I). We find that mark correlation functions are able to detect even a small dependence of galaxy properties on the environment, quantified by the overdensity $1+\delta$, while such a small dependence would be difficult to detect by traditional methods. We then show that rank ordering the marks and using the rank as a weight is a simple way of comparing the correlation signals for different marks. With this we quantify to what extent fixed-aperture overdensities are sensitive to large-scale halo environments, nearest-neighbor overdensities are sensitive to small-scale environments within haloes, and colour is a better tracer of overdensity than is luminosity.

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.

The influence of a galaxy's environment on its evolution has been studied and compared extensively in the literature, although differing techniques are often used to define environment. Most methods fall into two broad groups: those that use nearest neighbours to probe the underlying density field and those that use fixed apertures. The differences between the two inhibit a clean comparison between analyses and leave open the possibility that, even with the same data, different properties are actually being measured. In this work we apply twenty published environment definitions to a common mock galaxy catalogue constrained to look like the local Universe. We find that nearest neighbour-based measures best probe the internal densities of high-mass haloes, while at low masses the inter-halo separation dominates and acts to smooth out local density variations. The resulting correlation also shows that nearest neighbour galaxy environment is largely independent of dark matter halo mass. Conversely, aperture-based methods that probe super-halo scales accurately identify high-density regions corresponding to high mass haloes. Both methods show how galaxies in dense environments tend to be redder, with the exception of the largest apertures, but these are the strongest at recovering the background dark matter environment. We also warn against using photometric redshifts to define environment in all but the densest regions. When considering environment there are two regimes: the 'local environment' internal to a halo best measured with nearest neighbour and 'large-scale environment' external to a halo best measured with apertures. This leads to the conclusion that there is no universal environment measure and the most suitable method depends on the scale being probed.

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.

[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).

With the ever increasing resolution of N-body simulations, accurate subhalo detection is becoming essential in the study of the formation of structure, the production of merger trees and the seeding of semi-analytic models. To investigate the state of halo finders, we compare two different approaches to detecting subhaloes; the first based on overdensities in a halo and the second being adaptive mesh refinement. A set of stable mock NFW dark matter haloes were produced and a subhalo was placed at different radii within a larger halo. SUBFIND (a Friends-of-Friends based finder) and AHF (an adaptive mesh based finder) were employed to recover the subhalo. As expected, we found that the mass of the subhalo recovered by SUBFIND has a strong dependence on the radial position and that neither halo finder can accurately recover the subhalo when it is very near the centre of the halo. This radial dependence is shown to be related to the subhalo being truncated by the background density of the halo and originates due to the subhalo being defined as an overdensity. If the subhalo size is instead determined using the peak of the circular velocity profile, a much more stable value is recovered. The downside to this is that the maximum circular velocity is a poor measure of stripping and is affected by resolution. For future halo finders to recover all the particles in a subhalo, a search of phase space will need to be introduced.