The clustering of galaxy clusters in cosmological models with non-Gaussian initial conditions: predictions for future surveys

We predict the biasing and clustering properties of galaxy clusters that are expected to be observed in the catalogues produced by two forthcoming X-ray and Sunyaev–Zel'dovich effect surveys. We study a set of flat cosmological models where the primordial density probability distribution shows deviations from Gaussianity in agreement with current observational bounds form the background radiation. We consider both local and equilateral shapes for the primordial bispectrum in non-Gaussian models. The two catalogues investigated are those produced by the e ROSITA wide survey and from a survey based on South Pole Telescope observations. It turns out that both the bias and observed power spectrum of galaxy clusters are severely affected in non-Gaussian models with local shape of the primordial bispectrum, especially at large scales. On the other hand, models with equilateral shape of the primordial bispectrum show only a mild effect at all scales, that is difficult to be detected with clustering observations. Between the two catalogues, the one performing better is the e ROSITA one, since it contains only the largest masses that are more sensitive to primordial non-Gaussianity.     

Constraints on σ8 from galaxy clustering in N-body simulations and semi-analytic models

We generate mock galaxy catalogues for a grid of different cosmologies, using rescaled N -body simulations in tandem with a semi-analytic model run using consistent parameters. Because we predict the galaxy bias, rather than fitting it as a nuisance parameter, we obtain an almost pure constraint on σ₈ by comparing the projected two-point correlation function we obtain to that from the Sloan Digital Sky Survey (SDSS). A systematic error arises because different semi-analytic modelling assumptions allow us to fit the r -band luminosity function equally well. Combining our estimate of the error from this source with the statistical error, we find  σ₈= 0.97 ± 0.06  . We obtain consistent results if we use galaxy samples with a different magnitude threshold, or if we select galaxies by b _J-band rather than r -band luminosity and compare to data from the 2dF Galaxy Redshift Survey (2dFGRS). Our estimate for σ₈ is higher than that obtained for other analyses of galaxy data alone, and we attempt to find the source of this difference. We note that in any case, galaxy clustering data provide a very stringent constraint on galaxy formation models.     

MoMaF: the Mock Map Facility

We present the Mock Map Facility, a powerful tool for converting theoretical outputs of hierarchical galaxy formation models into catalogues of virtual observations. The general principle is straightforward: mock observing cones can be generated using semi-analytically post-processed snapshots of cosmological N -body simulations. These cones can then be projected to synthesize mock sky images. To this end, the paper describes in detail an efficient technique for creating such mock cones and images from the galaxies in cosmological simulations ( galics ) semi-analytic model, providing the reader with an accurate quantification of the artefacts it introduces at every step. We show that replication effects introduce a negative bias on the clustering signal – typically peaking at less than 10 per cent around the correlation length. We also thoroughly discuss how the clustering signal is affected by finite-volume effects, and show that it vanishes at scales larger than approximately one-tenth of the simulation box size. For the purpose of analysing our method, we show that number counts and redshift distributions obtained with galics / momaf compare well with K -band observations and the two-degree field galaxy redshift survey. Given finite-volume effects, we also show that the model can reproduce the automatic plate measuring machine angular correlation function. The momaf results discussed here are made publicly available to the astronomical community through a public data base. Moreover, a user-friendly Web interface ( http://galics.iap.fr ) allows any user to recover her/his own favourite galaxy samples through simple SQL queries. The flexibility of this tool should permit a variety of uses ranging from extensive comparisons between real observations and those predicted by hierarchical models of galaxy formation, to the preparation of observing strategies for deep surveys and tests of data processing pipelines.     

The three-dimensional power spectrum of dark and luminous matter from the VIRMOS-DESCART cosmic shear survey

We present the first optimal power spectrum estimation and three-dimensional deprojections for the dark and luminous matter and their cross-correlations. The results are obtained using a new optimal fast estimator, deprojected using minimum variance and Singular Value Decomposition (SVD) techniques. We show the resulting 3D power spectra for dark matter and galaxies, and their covariance for the VIRMOS-DESCART weak lensing shear and galaxy data. The survey is most sensitive to non-linear scales   k _NL∼ 1 h Mpc⁻¹  . On these scales, our 3D power spectrum of dark matter is in good agreement with the RCS 3D power spectrum found by Tegmark & Zaldarriaga. Our galaxy power is similar to that found by the 2MASS survey, and larger than that of SDSS, APM and RCS, consistent with the expected difference in galaxy population.
We find an average bias   b = 1.24 ± 0.18  for the I -selected galaxies, and a cross-correlation coefficient   r = 0.75 ± 0.23  . Together with the power spectra, these results optimally encode the entire two point information about dark matter and galaxies, including galaxy–galaxy lensing. We address some of the implications regarding galaxy haloes and mass-to-light ratios. The best-fitting 'halo' parameter   h ≡ r / b = 0.57 ± 0.16  , suggesting that dynamical masses estimated using galaxies systematically underestimate total mass.
Ongoing surveys, such as the Canada–France–Hawaii Telescope Legacy Survey, will significantly improve on the dynamic range, and future photometric redshift catalogues will allow tomography along the same principles.     

The influence of halo assembly on galaxies and galaxy groups

In this paper, we study the variations of group galaxy properties according to the assembly history in Sloan Digital Sky Survey Data Release 6 (SDSS-DR6) selected groups. Using mock SDSS group catalogues, we find two suitable indicators of group formation time: (i) the isolation of the group, defined as the distance to the nearest neighbour in terms of its virial radius and (ii) the concentration, measured as the group inner density calculated using the fifth nearest bright galaxy to the group centre. Groups within narrow ranges of mass in the mock catalogue show increasing group age with isolation and concentration. However, in the observational data the stellar age, as indicated by the spectral type, only shows a correlation with concentration.
We study groups of similar mass and different assembly history, finding important differences in their galaxy population. Particularly, in high-mass SDSS groups, the number of members, mass-to-light ratios, red galaxy fractions and the magnitude difference between the brightest and second-brightest group galaxies, show different trends as a function of isolation and concentration, even when it is expected that the latter two quantities correlate with group age. Conversely, low-mass SDSS groups appear to be less sensitive to their assembly history.
The correlations detected in the SDSS are not consistent with the trends measured in the mock catalogues. However, discrepancies can be explained in terms of the disagreement found in the age-isolation trends, suggesting that the model might be overestimating the effects of environment. We discuss how the modelling of the cold gas in satellite galaxies could be responsible for this problem. These results can be used to improve our understanding of the evolution of galaxies in high-density environments.     

The Hubble constant from galaxy lenses: impacts of triaxiality and model degeneracies

The Hubble constant can be constrained using the time delays between multiple images of gravitationally lensed sources. In some notable cases, typical lensing analyses assuming isothermal galaxy density profiles produce low values for the Hubble constant, inconsistent with the result of the HST Key Project  (72 ± 8 km s⁻¹ Mpc⁻¹)  . Possible systematics in the values of the Hubble constant derived from galaxy lensing systems can result from a number of factors, for example, neglect of environmental effects, assumption of isothermality, or contamination by line-of-sight structures. One additional potentially important factor is the triaxial structure of the lensing galaxy halo; most lens models account for halo shape simply by perturbing the projected spherical lensing potential, an approximation that is often necessary but that is inadequate at the levels of triaxiality predicted in the cold dark matter paradigm. To quantify the potential error introduced by this assumption in estimates of the Hubble parameter, we strongly lens a distant galaxy through a sample of triaxial softened isothermal haloes and use an Markov Chain Monte Carlo method to constrain the lensing halo profile and the Hubble parameter from the resulting multiple image systems. We explore the major degeneracies between the Hubble parameter and several parameters of the lensing model, finding that without a way to accurately break these degeneracies accurate estimates of the Hubble parameter are not possible. Crucially, we find that triaxiality does not significantly bias estimates of the Hubble constant, and offer an analytic explanation for this behaviour in the case of isothermal profiles. Neglected triaxial halo shape cannot contribute to the low Hubble constant values derived in a number of galaxy lens systems.     

Generalized cumulant correlators and hierarchical clustering

The cumulant correlators, C _pq , are statistical quantities that generalize the better-known S _p parameters; the former are obtained from the two-point probability distribution function of the density fluctuations while the latter describe only the one-point distribution. If galaxy clustering develops from Gaussian initial fluctuations and a small-angle approximation is adopted, standard perturbative methods suggest a particular hierarchical relationship of the C _pq for projected clustering data, such as that obtained from the Automatic Plate Measuring (APM) survey. We establish the usefulness of the two-point cumulants for describing hierarchical clustering by comparing such calculations against available measurements from projected catalogues, finding very good agreement. We extend the idea of cumulant correlators to multipoint generalized cumulant correlators (related to the higher-order correlation functions). We extend previous studies in the highly non-linear regime to express the generalized cumulant correlators in terms of the underlying 'tree amplitudes' of hierarchical scaling models. Such considerations lead to a technique for determining these hierarchical amplitudes, to arbitrary order, from galaxy catalogues and numerical simulations. Knowledge of these amplitudes yields important clues about the phenomenology of gravitational clustering. For instance, we show that a three-point cumulant correlator can be used to separate the tree amplitudes up to sixth order. We also combine the particular hierarchical Ansatz of Bernardeau & Schaeffer with extended and hyper-extended perturbation theory to infer values of the tree amplitudes in the highly non-linear regime.     

Weak lensing from strong clustering

We investigate the effect of weak gravitational lensing in the limit of small angular scales where projected galaxy clustering is strongly non-linear. This is the regime likely to be probed by future weak lensing surveys. We use well-motivated hierarchical scaling arguments and the plane-parallel approximation to study multi-point statistical properties of the convergence field. These statistics can be used to compute the vertex amplitudes in tree models of hierarchical clustering; these can be compared with similar measurements from galaxy surveys, leading to a powerful probe of galaxy bias.     

Apparent and actual galaxy cluster temperatures

The redshift evolution of the galaxy cluster temperature function is a powerful probe of cosmology. However, its determination requires the measurement of redshifts for all clusters in a catalogue, which is likely to prove challenging for large catalogues expected from XMM-Newton , which may contain of the order of 2000 clusters with measurable temperatures, distributed around the sky. In this paper we study the apparent cluster temperature, which can be obtained without cluster redshifts. We show that the apparent temperature function itself is of limited use in constraining cosmology, and so concentrate our focus on studying how apparent temperatures can be combined with other X-ray information to constrain the cluster redshift. We also briefly study the circumstances under which the non-thermal spectral features can provide redshift information.     

Large-scale bias in the Universe — II. Redshift-space bispectrum

The determination of the density parameter Ω₀ from the large-scale distribution of galaxies is one of the major goals of modern cosmology. However, if galaxies are biased tracers of the underlying mass distribution, linear perturbation theory leads to a degeneracy between Ω₀ and the linear bias parameter b , and the density parameter cannot be estimated. In Matarrese, Verde &38; Heavens we developed a method based on second-order perturbation theory to use the bispectrum to lift this degeneracy by measuring the bias parameter in an Ω₀-independent way. The formalism was developed assuming that one has perfect information on the positions of galaxies in three dimensions. In galaxy redshift surveys, the three-dimensional information is imperfect, because of the contaminating effects of peculiar velocities, and the resulting clustering pattern in redshift space is distorted. In this paper we combine second-order perturbation theory with a model for collapsed, virialized structures, to extend the method to redshift space, and demonstrate that the method should be successful in determining with reasonable accuracy the bias parameter from state-of-the-art surveys such as the Anglo-Australian 2 degree Field Survey and the Sloan Digital Sky Survey.     

Parameter information from non-linear cosmological fields

We develop a general formalism for analysing parameter information from non-Gaussian cosmic fields. The method can be adapted to include the non-linear effects in galaxy redshift surveys, weak lensing surveys and cosmic velocity field surveys as part of parameter estimation. It can also be used as a test of non-Gaussianity of the cosmic microwave background. Generalizing maximum-likelihood analysis to second order, we calculate the non-linear Fisher information matrix and likelihood surfaces in parameter space. To this order we find that the information content is always increased by including non-linearity. Our methods are applied to a realistic model of a galaxy redshift survey, including non-linear evolution, galaxy bias, shot-noise and redshift-space distortions to second order. We find that including non-linearities allows all of the degeneracies between parameters to be lifted. Marginalized parameter uncertainties of a few per cent will then be obtainable using forthcoming galaxy redshift surveys.     

Probing the gravity-induced bias with weak lensing: test of analytical results against simulations

Future weak lensing surveys will directly probe the density fluctuation in the Universe. Recent studies have shown how the statistics of the weak lensing convergence field is related to the statistics of collapsed objects. Extending earlier analytical results on the probability distribution function of the convergence field, we show that the bias associated with the convergence field can be directly related to the bias associated with the statistics of underlying overdense objects. This will provide us with a direct method to study the gravity-induced bias in galaxy clustering. Based on our analytical results, which use the hierarchical Ansatz for non-linear clustering, we study how such a bias depends on the smoothing angle and the source redshift. We compare our analytical results with ray-tracing experiments through N -body simulations of four different realistic cosmological scenarios, and find a very good match. Our study shows that the bias in the convergence map strongly depends on the background geometry and hence can help us in distinguishing different cosmological models in addition to improving our understanding of the gravity-induced bias in galaxy clustering.     

The spatial evolution of stellar structures in the Large Magellanic Cloud

We present an analysis of the spatial distribution of various stellar populations within the Large Magellanic Cloud (LMC). We combine mid-infrared selected young stellar objects, optically selected samples with mean ages between ∼9 and ∼1000 Myr and existing stellar cluster catalogues to investigate how stellar structures form and evolve within the LMC. For the analysis we use Fractured Minimum Spanning Trees, the statistical Q parameter and the two-point correlation function. Restricting our analysis to young massive (OB) stars, we confirm our results obtained for M33, namely that the luminosity function of the groups is well described by a power law with index −2, and that there is no characteristic length-scale of star-forming regions. We find that stars in the LMC are born with a large amount of substructure, consistent with a two-dimensional fractal distribution with dimension     and evolve towards a uniform distribution on a time-scale of ∼175 Myr. This is comparable to the crossing time of the galaxy, and we suggest that stellar structure, regardless of spatial scale, will be eliminated in a crossing time. This may explain the smooth distribution of stars in massive/dense young clusters in the Galaxy, while other, less massive, clusters still display large amounts of structure at similar ages. By comparing the stellar and star cluster distributions and evolving time-scales, we show that infant mortality of clusters (or 'popping clusters') has a negligible influence on the galactic structure. Finally, we quantify the influence of the elongation, differential extinction and contamination of a population on the measured Q value.     

Scale-dependent galaxy bias in the Sloan Digital Sky Survey as a function of luminosity and colour

It has been known for a long time that the clustering of galaxies changes as a function of galaxy type. This galaxy bias acts as a hindrance to the extraction of cosmological information from the galaxy power spectrum or correlation function. Theoretical arguments show that a change in the amplitude of the clustering between galaxies and mass on large scales is unavoidable, but cosmological information can be easily extracted from the shape of the power spectrum or correlation function if this bias is independent of scale. Scale-dependent bias is generally small on large scales,   k < 0.1  h  Mpc⁻¹  , but on smaller scales can affect the recovery of  Ω_m h   from the measured shape of the clustering signal, and have a small effect on the Baryon Acoustic Oscillations. In this paper, we investigate the transition from scale-independent to scale-dependent galaxy bias as a function of galaxy population. We use the Sloan Digital Sky Survey Data Release 5 sample to fit various models, which attempt to parametrize the turn-off from scale-independent behaviour. For blue galaxies, we find that the strength of the turn-off is strongly dependent on galaxy luminosity, with stronger scale-dependent bias on larger scales for more luminous galaxies. For red galaxies, the scale dependence is a weaker function of luminosity. Such trends need to be modelled in order to optimally extract the information available in future surveys, and can help with the design of such surveys.     

-ray wakes as probes of galaxy cluster dynamics

If a galaxy resides in a cluster, then its passage through the pervasive intracluster medium will produce a detectable signature in the X-ray emission from the cluster. Such features have now been detected in a number of systems. The simplest kinematic information that can be extracted from this signature is the galaxy's direction of motion on the plane of the sky. This paper explores the constraints on cluster dynamics that could be derived from such information. In particular, we show that it is possible to define a projected anisotropy parameter, B ( r ), which is directly analogous to the usual orbital anisotropy parameter. We describe an estimator for this quantity, ( R ), which can be derived in a robust and straightforward manner. We present a simple dynamical model for a cluster consisting of a Michie distribution function of galaxies orbiting in a truncated singular isothermal sphere potential. Using this model, we demonstrate the ambiguity between the distribution of mass and the distribution of galaxy orbits when interpreting the traditional measures of cluster kinematics (the projected density of galaxies and their line-of-sight velocity dispersion). As an example, we show how two very different dynamical models can fit the kinematic properties of the Coma cluster. We demonstrate that the measurement of using a relatively small sample of wake directions ( N _wake≈50) would provide an effective mechanism for lifting this degeneracy. Thus, by combining X-ray measurements of wake directions with number counts and line-of-sight velocities derived from optical data, it will prove possible to measure both the orbit distribution and the form of the gravitational potential in clusters of galaxies. The requisite X-ray observations lie within reach of the soon-to-be-launched AXAF satellite.     

Searching for the scale of homogeneity

We introduce a statistical quantity, known as the K function, related to the integral of the two-point correlation function. It gives us straightforward information about the scale where clustering dominates and the scale at which homogeneity is reached. We evaluate the correlation dimension, D ₂, as the local slope of the log–log plot of the K function. We apply this statistic to several stochastic point fields, to three numerical simulations describing the distribution of clusters and finally to real galaxy redshift surveys. Four different galaxy catalogues have been analysed using this technique: the Center for Astrophysics I, the Perseus–Pisces redshift surveys (these two lying in our local neighbourhood), the Stromlo–APM and the 1.2-Jy IRAS redshift surveys (these two encompassing a larger volume). In all cases, this cumulant quantity shows the fingerprint of the transition to homogeneity. The reliability of the estimates is clearly demonstrated by the results from controllable point sets, such as the segment Cox processes. In the cluster distribution models, as well as in the real galaxy catalogues, we never see long plateaus when plotting D ₂ as a function of the scale, leaving no hope for unbounded fractal distributions.     

Peculiar velocities and the mean density parameter

We study the peculiar velocity field inferred from the Mark III spirals using a new method of analysis. We estimate optimal values of Tully–Fisher scatter and zero-point offset, and we derive the three-dimensional rms peculiar velocity ( σ _v ) of the galaxies in the samples analysed. We check our statistical analysis using mock catalogues derived from numerical simulations of cold dark matter (CDM) models considering measurement uncertainties and sampling variations. Our best determination for the observations is σ _v =(660±50) km s⁻¹. We use the linear theory relation between σ _v , the density parameter Ω, and the galaxy correlation function ξ ( r ) to infer the quantity     , where b is the linear bias parameter of optical galaxies and the uncertainties correspond to bootstrap resampling and an estimated cosmic variance added in quadrature. Our findings are consistent with the results of cluster abundances and redshift-space distortion of the two-point correlation function. These statistical measurements suggest a low value of the density parameter Ω∼0.4 if optical galaxies are not strongly biased tracers of mass.