Search

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., Comment: 14 pages, 9 figures, accepted for publication in MNRAS

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., Comment: 9 pages, 8 figures, Accepted to MNRAS (24th May 2017)

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., Comment: 7 pages, 4 figures, accepted for publication in MNRAS

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., Comment: 14 pages, 14 figures, accepted to MNRAS

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., Comment: Figs. 2, 4, and 7 show the main issues. This project was initiated at the Subhaloes Going Notts conference (http://popia.ft.uam.es/SubhaloesGoingNotts/Home.html). MNRAS submitted

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., Comment: 10 pages; 5 figures

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., Comment: 15 pages, 11 figures, 2 tables. 1 figure removed to reduce paper extension, shorter discussion, accepted for publication in MNRAS

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., Comment: 19 pages, 15 figures

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., Comment: 10 pages, 7 figures, submitted to MNRAS, comments welcome

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., Comment: 12 pages, 9 figures, accepted to MNRAS

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., Comment: 19 pages, 12 figures, accepted for publication in MNRAS

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., Comment: 28 pages containing 13 figures & 4 tables + 9 pages appendix containing another 4 tables + 4 pages of references, accepted for publication in MNRAS

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., Comment: 10 pages, 14 figures and 3 tables. Accepted for publication in MNRAS

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., Comment: 12 pages, 15 figures; published in MNRAS. The mock galaxy catalogues and environment measures used in this paper and Paper I (arXiv:1109.6328) are available here: http://tao.it.swin.edu.au/mock-galaxy-factory/precomputed-mock-catalogues/

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., Comment: 14 interesting pages, 17 beautiful figures, and 2 informative tables accepted for publication in MNRAS (matches published version)

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., Comment: 16 pages, 7 figures, 2 tables, Accepted for MNRAS

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., Comment: 14 pages, 9 figures, 1 table, published in MNRAS

[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)., Comment: 27 interesting pages, 20 beautiful figures, and 4 informative tables accepted for publication in MNRAS. The high-resolution version of the paper as well as all the test cases and analysis can be found at the web site http://popia.ft.uam.es/HaloesGoingMAD

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., Comment: 9 pages, 7 figures, accepted for publication in MNRAS