The history of mass assembly of faint red galaxies in 28 galaxy clusters since z= 1.3

We measure the relative evolution of the number of bright and faint (as faint as  0.05 L *)  red galaxies in a sample of 28 clusters, out of which 16 are at  0.50 ≤ z ≤ 1.27  , all observed through a pair of filters bracketing the 4000-Å break rest frame. The abundance of red galaxies, relative to bright ones, is constant over all the studied redshift range,  0 < z < 1.3  , and rules out a differential evolution between bright and faint red galaxies as large as claimed in some past works. Faint red galaxies are largely assembled and in place at   z = 1.3  and their abundance does not depend on cluster mass, parametrized by velocity dispersion or X-ray luminosity. Our analysis, with respect to the previous one, samples a wider redshift range, minimizes systematics and put a more attention to statistical issues, keeping at the same time a large number of clusters.     

Survival of substructure within dark matter haloes

Using high-resolution cosmological N -body simulations, we investigate the survival of dark matter satellites falling into larger haloes. Satellites preserve their identity for some time after merging. We compute their loss of mass, energy and angular momentum as they are dissolved by dynamical friction, tidal forces and collisions with other satellites. We also analyse the evolution of their internal structure. Satellites with less than a few per cent of the mass of the main halo may survive for several billion years, whereas larger satellites rapidly sink into the centre of the main halo potential well and lose their identity. Penetrating encounters between satellites are frequent and may lead to significant mass loss and disruption. Only a minor fraction of cluster mass (10–15 per cent on average) is bound to substructure at most redshifts of interest. We discuss the application of these results to the survival and extent of dark matter haloes associated with galaxies in clusters, and to their interactions. We find that a minor fraction of galaxy-size dark matter haloes are disrupted by redshift z  = 0. The fraction of satellites undergoing close encounters is similar to the observed fraction of interacting or merging galaxies in clusters at moderate redshift.     

Simulating subhaloes at high redshift: merger rates, counts and types

Galaxies are believed to be in one-to-one correspondence with simulated dark matter subhaloes. We use high-resolution N -body simulations of cosmological volumes to calculate the statistical properties of subhalo (galaxy) major mergers at high redshift ( z = 0.6–5). We measure the evolution of the galaxy merger rate, finding that it is much shallower than the merger rate of dark matter host haloes at   z > 2.5  , but roughly parallels that of haloes at   z < 1.6  . We also track the detailed merger histories of individual galaxies and measure the likelihood of multiple mergers per halo or subhalo. We examine satellite merger statistics in detail: 15–35 per cent of all recently merged galaxies are satellites, and satellites are twice as likely as centrals to have had a recent major merger. Finally, we show how the differing evolution of the merger rates of haloes and galaxies leads to the evolution of the average satellite occupation per halo, noting that for a fixed halo mass, the satellite halo occupation peaks at   z ∼ 2.5  .     

The formation and survival of discs in a ΛCDM universe

We study the formation of galaxies in a Λ cold dark matter (ΛCDM) universe using high-resolution hydrodynamical simulations with a multiphase treatment of gas, cooling and feedback, focusing on the formation of discs. Our simulations follow eight isolated haloes similar in mass to the Milky Way and extracted from a large cosmological simulation without restriction on spin parameter or merger history. This allows us to investigate how the final properties of the simulated galaxies correlate with the formation histories of their haloes. We find that, at   z = 0  , none of our galaxies contains a disc with more than 20 per cent of its total stellar mass. Four of the eight galaxies nevertheless have well-formed disc components, three have dominant spheroids and very small discs, and one is a spheroidal galaxy with no disc at all. The   z = 0  spheroids are made of old stars, while discs are younger and formed from the inside-out. Neither the existence of a disc at   z = 0  nor the final disc-to-total mass ratio seems to depend on the spin parameter of the halo. Discs are formed in haloes with spin parameters as low as 0.01 and as high as 0.05; galaxies with little or no disc component span the same range in spin parameter. Except for one of the simulated galaxies, all have significant discs at   z ≳ 2  , regardless of their   z = 0  morphologies. Major mergers and instabilities which arise when accreting cold gas is misaligned with the stellar disc trigger a transfer of mass from the discs to the spheroids. In some cases, discs are destroyed, while in others, they survive or reform. This suggests that the survival probability of discs depends on the particular formation history of each galaxy. A realistic ΛCDM model will clearly require weaker star formation at high redshift and later disc assembly than occurs in our models.     

The prevalence of cooling cores in clusters of galaxies at z≈ 0.15–0.4

We present a Chandra study of 38 X-ray-luminous clusters of galaxies in the ROSAT Brightest Cluster Sample (BCS) that lie at moderate redshifts  ( z ≈ 0.15–0.4)  . Based primarily on power ratios and temperature maps, we find that the majority of clusters at moderate redshift generally have smooth, relaxed morphologies with some evidence for mild substructure perhaps indicative of recent minor merger activity. Using spatially resolved spectral analyses, we find that cool cores appear still to be common at moderate redshift. At a radius of 50 kpc, we find that at least 55 per cent of the clusters in our sample exhibit signs of mild cooling  ( t _cool < 10 Gyr)  , while in the central bin at least 34 per cent demonstrate signs of strong cooling  ( t _cool < 2 Gyr)  . These percentages are nearly identical to those found for luminous, low-redshift clusters of galaxies, indicating that there appears to be little evolution in cluster cores since   z ≈ 0.4  and suggesting that heating and cooling mechanisms may already have stabilized by this epoch. Comparing the central cooling times to catalogues of central Hα emission in BCS clusters, we find a strong correspondence between the detection of Hα and central cooling time. We also confirm a strong correlation between the central cooling time and cluster power ratios, indicating that crude morphological measures can be used as a proxy for more rigorous analysis in the face of limited signal-to-noise ratio data. Finally, we find that the central temperatures for our sample typically drop by no more than a factor of ∼3–4 from the peak cluster temperatures, similar to those of many nearby clusters.     

Substructures in hydrodynamical cluster simulations

We present predictions for the abundance and nature of extremely red objects (EROs) in the Λ cold dark matter model. EROs are red, massive galaxies observed at   z ≥ 1  and their numbers and properties pose a challenge to hierarchical galaxy formation models. We compare the predictions from two published models, one of which invokes a 'superwind' to regulate star formation in massive haloes and the other which suppresses gas cooling in haloes through 'radio-mode' active galactic nucleus (AGN) feedback. The superwind model underestimates the number counts of EROs by an order of magnitude, whereas the radio-mode AGN feedback model gives excellent agreement with the number counts and redshift distribution of EROs. In the AGN feedback model the ERO population is dominated by old, passively evolving galaxies, whereas observations favour an equal split between old galaxies and dusty starbursts. Also, the model predicts a more extended redshift distribution of passive galaxies than is observed. These comparisons suggest that star formation may be quenched too efficiently in this model.     

Richness dependence of the recent evolution of clusters of galaxies

We revisit the issue of the recent dynamical evolution of clusters of galaxies using a sample of Abell, Corwin & Olowin (ACO) clusters with   z < 0.14  , which has been selected such that it does not contain clusters with multiple velocity components nor strongly merging or interacting clusters, as revealed in X-rays. We use as proxies of the cluster dynamical state the projected cluster ellipticity, velocity dispersion and X-ray luminosity. We find indications for a recent dynamical evolution of this cluster population, which however strongly depends on the cluster richness. Poor clusters appear to be undergoing their primary phase of virialization, with their ellipticity increasing with redshift with a rate  dε/d z ≃ 2.5 ± 0.4  , while the richest clusters show an ellipticity evolution in the opposite direction (with  dε/d z ≃−1.2 ± 0.1  ), which could be due to secondary infall. When taking into account sampling effects due to the magnitude-limited nature of the ACO cluster catalogue we find no significant evolution of the cluster X-ray luminosity, while the velocity dispersion increases with decreasing redshift, independent of the cluster richness, at a rate  dσ _v /d z ≃−1700 ± 400 km s⁻¹  .     

The evolution of the Lyα forest effective optical depth following He ii reionization

Three independent observational studies have now detected a narrow  (Δ z ≃ 0.5)  dip centred at   z = 3.2  in the otherwise smooth redshift evolution of the Lyα forest effective optical depth. This feature has previously been interpreted as an indirect signature of rapid photoheating in the intergalactic medium (IGM) during the epoch of He  ii reionization. We examine this interpretation using a semi-analytic model of inhomogeneous He  ii reionization and high-resolution hydrodynamical simulations of the Lyα forest. We instead find that a rapid  (Δ z ≃ 0.2)  boost to the IGM temperature  (Δ T ≃ 10⁴ K)  beginning at   z = 3.4  produces a well understood and generic evolution in the Lyα effective optical depth, where a sudden reduction in the opacity is followed by a gradual, monotonic recovery driven largely by adiabatic cooling in the low-density IGM. This behaviour is inconsistent with the narrow feature in the observational data. If photoheating during He  ii reionization is instead extended over several redshift units, as recent theoretical studies suggest, then the Lyα opacity will evolve smoothly with redshift. We conclude that the sharp dip observed in the Lyα forest effective optical depth is instead most likely due to a narrow peak in the hydrogen photoionization rate around   z = 3.2  , and suggest that it may arise from the modulation of either reprocessed radiation during He  ii reionization, or the opacity of Lyman limit systems.     

The evolution of the cluster X-ray scaling relations in the Wide Angle ROSAT Pointed Survey sample at 0.6 < z < 1.0

The X-ray properties of a sample of 11 high-redshift  (0.6 < z < 1.0)  clusters observed with Chandra and/or XMM–Newton are used to investigate the evolution of the cluster scaling relations. The observed evolution in the normalization of the   L – T , M – T , M _g– T   and M – L relations is consistent with simple self-similar predictions, in which the properties of clusters reflect the properties of the Universe at their redshift of observation. Under the assumption that the model of self-similar evolution is correct and that the local systems formed via a single spherical collapse, the high-redshift L – T relation is consistent with the high- z clusters having virialized at a significantly higher redshift than the local systems. The data are also consistent with the more realistic scenario of clusters forming via the continuous accretion of material.
The slope of the L – T relation at high redshift  ( B = 3.32 ± 0.37)  is consistent with the local relation, and significantly steeper than the self-similar prediction of   B = 2  . This suggests that the same non-gravitational processes are responsible for steepening the local and high- z relations, possibly occurring universally at   z ≳ 1  or in the early stages of the cluster formation, prior to their observation.
The properties of the intracluster medium at high redshift are found to be similar to those in the local Universe. The mean surface-brightness profile slope for the sample is  β= 0.66 ± 0.05  , the mean gas mass fractions within   R _{2500( z )}  and   R _{200( z )}  are  0.069 ± 0.012  and  0.11 ± 0.02  , respectively, and the mean metallicity of the sample is  0.28 ± 0.11 Z_⊙  .     

The low-redshift evolution of the Lyman-β forest

The low-redshift evolution of the intergalactic medium is investigated using hydrodynamic cosmological simulations. The assumed cosmological model is a critical density cold dark matter universe. The imposed uniform background of ionizing radiation has the amplitude, shape and redshift evolution as computed from the observed quasar luminosity function by Haardt &38; Madau. We have analysed simulated Lyman-α spectra using Voigt-profile fitting, mimicking the procedure with which quasar spectra are analysed. Our simulations reproduce the observed evolution of the number of Lyman-α absorption lines over the whole observed interval of z  = 0.5 to 4. In particular, our simulations show that the decrease in the rate of evolution of Lyman-α absorption lines at z  ≤ 2, as observed by the Hubble Space Telescope , can be explained by the steep decline in the photoionizing background resulting from the rapid decline in quasar numbers at low redshift.     

ΛCDM predictions for galaxy protoclusters – I. The relation between galaxies, protoclusters and quasars at z∼ 6

Motivated by recent observational studies of the environment of   z ∼ 6  QSOs, we have used the Millennium Run (MR) simulations to construct a very large  (∼4°× 4°)  mock redshift survey of star-forming galaxies at   z ∼ 6  . We use this simulated survey to study the relation between density enhancements in the distribution of i ₇₇₅-dropouts and Lyα emitters, and their relation to the most massive haloes and protocluster regions at   z ∼ 6  . Our simulation predicts significant variations in surface density across the sky with some voids and filaments extending over scales of 1°, much larger than probed by current surveys. Approximately one-third of all   z ∼ 6  haloes hosting i -dropouts brighter than   z = 26.5  mag  (≈ M *_{UV, z =6})  become part of   z = 0  galaxy clusters. i -dropouts associated with protocluster regions are found in regions where the surface density is enhanced on scales ranging from a few to several tens of arcminutes on the sky. We analyse two structures of i -dropouts and Lyα emitters observed with the Subaru Telescope and show that these structures must be the seeds of massive clusters in formation. In striking contrast, six   z ∼ 6  QSO fields observed with Hubble Space Telescope show no significant enhancements in their i ₇₇₅-dropout number counts. With the present data, we cannot rule out the QSOs being hosted by the most massive haloes. However, neither can we confirm this widely used assumption. We conclude by giving detailed recommendations for the interpretation and planning of observations by current and future ground- and space-based instruments that will shed new light on questions related to the large-scale structure at   z ∼ 6  .     

The impact of dust on the scaling properties of galaxy clusters

We investigate the effect of dust on the scaling properties of galaxy clusters based on hydrodynamic N -body simulations of structure formation. We have simulated five dust models plus radiative cooling and adiabatic models using the same initial conditions for all runs. The numerical implementation of dust was based on the analytical computations of Montier & Giard. We set up dust simulations to cover different combinations of dust parameters that make evident the effects of size and abundance of dust grains. Comparing our radiative plus dust cooling runs with a purely radiative cooling simulation, we find that dust has an impact on cluster scaling relations. It mainly affects the normalization of the scalings (and their evolution), whereas it introduces no significant differences in their slopes. The strength of the effect critically depends on the dust abundance and grain size parameters as well as on the cluster scaling. Indeed, cooling due to dust is effective in the cluster regime and has a stronger effect on the 'baryon driven' statistical properties of clusters such as   L _X– M , Y – M , S – M   scaling relations. Major differences, relative to the radiative cooling model, are as high as 25 per cent for the   L _X– M   normalization, and about 10 per cent for the Y – M and S – M normalizations at redshift zero. On the other hand, we find that dust has almost no impact on the 'dark matter driven'   T _mw– M   scaling relation. The effects are found to be dependent in equal parts on both dust abundances and grain size distributions for the scalings investigated in this paper. Higher dust abundances and smaller grain sizes cause larger departures from the radiative cooling (i.e. with no dust) model.     

The evolution of [O ii] emission from cluster galaxies

We investigate the evolution of the star formation rate in cluster galaxies. We complement data from the Canadian Network for Observational Cosmology 1 (CNOC1) cluster survey  (0.15 < z < 0.6)  with measurements from galaxy clusters in the Two-degree Field (2dF) galaxy redshift survey  (0.05 < z < 0.1)  and measurements from recently published work on higher-redshift clusters, up to almost   z = 1  . We focus our attention on galaxies in the cluster core, i.e. galaxies with   r < 0.7  h ⁻¹₇₀ Mpc  . Averaging over clusters in redshift bins, we find that the fraction of galaxies with strong [O  ii ] emission is ≲20 per cent in cluster cores, and the fraction evolves little with redshift. In contrast, field galaxies from the survey show a very strong increase over the same redshift range. It thus appears that the environment in the cores of rich clusters is hostile to star formation at all the redshifts studied. We compare this result with the evolution of the colours of galaxies in cluster cores, first reported by Butcher and Oemler. Using the same galaxies for our analysis of the [O  ii ] emission, we confirm that the fraction of blue galaxies, which are defined as galaxies 0.2 mag bluer in the rest-frame B – V than the red sequence of each cluster, increases strongly with redshift. Because the colours of galaxies retain a memory of their recent star formation history, while emission from the [O  ii ] line does not, we suggest that these two results can best be reconciled if the rate at which the clusters are being assembled is higher in the past, and the galaxies from which it is being assembled are typically bluer.     

Simulating the formation of a protocluster at z∼ 2

We present results from two high-resolution hydrodynamical simulations of protocluster regions at   z ≃ 2.1  . The simulations have been compared to observational results for the so-called Spiderweb galaxy system, the core of a putative protocluster region at   z = 2.16  , found around a radio galaxy. The simulated regions have been chosen so as to form a poor cluster with   M ₂₀₀≃ 10¹⁴  h ⁻¹ M_⊙  (C1) and a rich cluster with   M ₂₀₀≃ 2 × 10¹⁵  h ⁻¹ M_⊙  (C2) at   z = 0  . The simulated protoclusters show evidence of ongoing assembly of a dominating central galaxy. The stellar mass of the brightest cluster galaxy of the C2 system is in excess with respect to observational estimates for the Spiderweb galaxy, with a total star formation rate which is also larger than indicated by observations. We find that the projected velocities of galaxies in the C2 cluster are consistent with observations, while those measured for the poorer cluster C1 are too low compared with the observed velocities. We argue that the Spiderweb complex resembles the high-redshift progenitor of a rich galaxy cluster. Our results indicate that the included supernovae feedback is not enough to suppress star formation in these systems, supporting the need of introducing active galactic nuclei feedback. According to our simulations, a diffuse atmosphere of hot gas in hydrostatic equilibrium should already be present at this redshift, and enriched at a level comparable to that of nearby galaxy clusters. The presence of this gas should be detectable with future deep X-ray observations.     

Angular momentum and clustering properties of early dark matter haloes

In this paper, we study the angular momentum properties of simulated dark matter haloes at high redshifts that likely host the first stars in the Universe. Calculating the spin distributions of these  10⁶– 10⁷ M_⊙  haloes in redshift slices from   z = 15  to 6, we find that they are well fit by a lognormal distribution as is found for lower redshift and more massive haloes in earlier work. We find that both the mean value of the spin and dispersion are largely unchanged with redshift for all haloes. Our key result is that subsamples of low- and high-spin, 10⁶ and  10⁷ M_⊙  , haloes show difference in clustering strength. In both mass bins, higher spin haloes are more strongly clustered in concordance with a tidal torquing picture for the growth of angular momentum in dark matter haloes in the cold dark matter paradigm.     

Recent arrival of faint cluster galaxies on the red sequence: luminosity functions from 119 deg2 of CFHTLS

The global star formation rate has decreased significantly since   z ∼ 1  , for reasons that are not well understood. Red-sequence galaxies, dominating in galaxy clusters, represent the population that have had their star formation shut off, and may therefore be the key to this problem. In this work, we select 127 rich galaxy clusters at  0.17 ≤ z ≤ 0.36  , from 119 deg² of the Canada–France–Hawaii Telescope Legacy Survey (CFHTLS) optical imaging data, and construct the r '-band red-sequence luminosity functions (LFs). We show that the faint end of the LF is very sensitive to how red-sequence galaxies are selected, and an optimal way to minimize the contamination from the blue cloud is to mirror galaxies on the redder side of the colour–magnitude relation. The LFs of our sample have a significant inflexion centred at     , suggesting a mixture of two populations. Combining our survey with low-redshift samples constructed from the Sloan Digital Sky Survey, we show that there is no strong evolution of the faint end of the LF (or the red-sequence dwarf-to-giant ratio) over the redshift range  0.2 ≲ z ≲ 0.4  , but from   z ∼ 0.2  to ∼0 the relative number of red-sequence dwarf galaxies has increased by a factor of ∼3, implying a significant build-up of the faint end of the cluster red sequence over the last 2.5 Gyr.     

Tracing the nature of dark energy with galaxy distribution

Dynamical dark energy (DE) is a viable alternative to the cosmological constant. Constructing tests to discriminate between Λ and dynamical DE models is difficult, however, because the differences are not large. In this paper we explore tests based on the galaxy mass function, the void probability function (VPF), and the number of galaxy clusters. At high z , the number density of clusters shows large differences between DE models, but geometrical factors reduce the differences substantially. We find that detecting a model dependence in the cluster redshift distribution is a significant challenge. We show that the galaxy redshift distribution is potentially a more sensitive characteristic. We do this by populating dark matter haloes in N -body simulations with galaxies using well-tested halo occupation distributions. We also estimate the VPF and find that samples with the same angular surface density of galaxies, in different models, exhibition almost model-independent VPF which therefore cannot be used as a test for DE. Once again, geometry and cosmic evolution compensate each other. By comparing VPFs for samples with fixed galaxy mass limits, we find measurable differences.     

The WARPS survey – IV. The X-ray luminosity–temperature relation of high-redshift galaxy clusters

We present a measurement of the cluster X-ray luminosity–temperature ( L – T ) relation out to high redshift ( z ∼0.8). Combined ROSAT PSPC spectra of 91 galaxy clusters detected in the Wide Angle ROSAT Pointed Survey (WARPS) are simultaneously fitted in redshift and luminosity bins. The resulting temperature and luminosity measurements of these bins, which occupy a region of the high-redshift L – T relation not previously sampled, are compared with existing measurements at low redshift in order to constrain the evolution of the L – T relation. We find the best fit to low-redshift ( z <0.2) cluster data, at T >1 keV, to be L ∝ T ^3.15±0.06. Our data are consistent with no evolution in the normalization of the L – T relation up to z ∼0.8. Combining our results with ASCA measurements taken from the literature, we find η =0.19±0.38 (for Ω₀=1, with 1 σ errors) where L _Bol∝(1+ z ) ^η T ^3.15, or η =0.60±0.38 for Ω₀=0.3. This lack of evolution is considered in terms of the entropy-driven evolution of clusters. Further implications for cosmological constraints are also discussed.     

Galaxy groups in the 2dF Galaxy Redshift Survey: the catalogue

We construct a galaxy groups catalogue from the public 100-K data release of the 2dF Galaxy Redshift Survey. The group identification is carried out using a slightly modified version of the group-finding algorithm developed by Huchra & Geller. Several tests using mock catalogues allow us to find the optimal conditions to increase the reliability of the final group sample. A minimum number of four members, an outer number density enhancement of 80 and a linking radial cut-off of 200 km s⁻¹ are the best obtained values from the analysis. Using these parameters, approximately 90 per cent of groups identified in real space have a redshift space counterpart. On the other hand, the level of contamination in redshift space reaches 30 per cent, including ∼6 per cent of artificial groups and ∼24 per cent of groups associated with binaries or triplets in real space. The final sample comprises 2209 galaxy groups covering the sky region described by Colless et al.spanning over the redshift range of  0.003 ≤ z ≤ 0.25  with a mean redshift of 0.1.