Stuart I. Muldrew

Haloes gone MAD: The Halo-Finder Comparison Project

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, spherical-overdensity and phase-space-based algorithms. We

Measures of galaxy environment – I. What is ‘environment’?

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.

Subhaloes going Notts: the subhalo-finder comparison project

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 in the presence and location of substructure and even for properties relating to the inner 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 the basic properties (mass and 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.

The Accuracy of Subhalo Detection

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 Navarro-Frenk-White (NFW) dark matter haloes was 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.

Galaxy cluster mass reconstruction project – I. Methods and first results on galaxy-based techniques

This paper is the first in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilize 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 idealized 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(aS (TM)), recovered cluster masses are correlated with the true underlying cluster mass with an intrinsic scatter of typically a factor of 2. Below 10(14)M(aS (TM)), 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.

What are protoclusters? – Defining high-redshift galaxy clusters and protoclusters

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∼35 h−1 Mpc commoving at z =2 (∼30 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 are 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.

A z = 2.5 protocluster associated with the radio galaxy MRC 2104−242: star formation and differing mass functions in dense environments

We present results from a narrow-band survey of the field around the high-redshift radio galaxy MRC 2104−242. We have selected Hα emitters in a 7 arcmin2 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 1010.5M_) galaxies, the density of the protocluster field increases to ∼55 times the control field density.

Galaxy Cluster Mass Reconstruction Project – II. Quantifying scatter and bias using contrasting mock catalogues

We would like to acknowledge funding from the Science and Technology Facilities Council (STFC). DC would like to thank the Australian Research Council for receipt of a QEII Research Fellowship. The Dark Cosmology Centre is funded by the Danish National Research Foundation. The authors would like to express special thanks to the Instituto de Fisica Teorica (IFT-UAM/CSIC in Madrid) for its hospitality and support, via the Centro de Excelencia Severo Ochoa Program under Grant No. SEV-2012-0249, during the three week workshop “nIFTy Cosmology” where this work developed. We further acknowledge the financial support of the University of Western 2014 Australia Research Collaboration Award for “Fast Approximate Synthetic Universes for the SKA”, the ARC Centre of Excellence for All Sky Astrophysics (CAASTRO) grant number CE110001020, and the two ARC Discovery Projects DP130100117 and DP140100198. We also recognise support from the Universidad Autonoma de Madrid (UAM) for the workshop infrastructure. RAS acknowledges support from the NSF grant AST-1055081. CS acknowledges support from the European Research Council under FP7 grant number 279396. SIM acknowledges the support of the STFC consolidated grant (ST/K001000/1) to the astrophysics group at the University of Leicester. ET acknowledge the support from the ESF grant IUT40-2.

Subhaloes gone Notts: spin across subhaloes and finders

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 10 different substructure finders, which span a wide range of methodologies from simple overdensity in configuration space to full 6D 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 substructures the nearer they are to the host haloes, 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.

Major mergers going Notts: challenges for modern halo finders

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.