Subhaloes gone Notts: Subhaloes as tracers of the dark matter halo shape

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.

Subhaloes gone Notts: the clustering properties of subhaloes

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.

Are the halo occupation predictions consistent with large scale galaxy clustering?

We study how well we can reconstruct the 2-point clustering of galaxies on linear scales, as a function of mass and luminosity, using the halo occupation distribution (HOD) in several semi-analytical models (SAMs) of galaxy formation from the Millennium Simulation. We find that HOD with Friends of Friends groups can reproduce galaxy clustering better than gravitationally bound haloes. This indicates that Friends of Friends groups are more directly related to the clustering of these regions than the bound particles of the overdensities. In general we find that the reconstruction works at best to 5% accuracy: it underestimates the bias for bright galaxies. This translates to an overestimation of 50% in the halo mass when we use clustering to calibrate mass. We also found a degeneracy on the mass prediction from the clustering amplitude that affects all the masses. This effect is due to the clustering dependence on the host halo substructure, an indication of assembly bias. We show that the clustering of haloes of a given mass increases with the number of subhaloes, a result that only depends on the underlying matter distribution. As the number of galaxies increases with the number of subhaloes in SAMs, this results in a low bias for the HOD reconstruction. We expect this effect to apply to other models of galaxy formation, including the real universe, as long as the number of galaxies incresases with the number of subhaloes. We have also found that the reconstructions of galaxy bias from the HOD model fails for low mass haloes with M = 3-5x10^11 Msun/h. We find that this is because galaxy clustering is more strongly affected by assembly bias for these low masses.

Statistical Analysis of Galaxy Surveys-IV: An objective way to quantify the impact of superstructures on galaxy clustering statistics

For galaxy clustering to provide robust constraints on cosmological parameters and galaxy formation models, it is essential to make reliable estimates of the errors on clustering measurements. We present a new technique, based on a spatial Jackknife (JK) resampling, which provides an objective way to estimate errors on clustering statistics. Our approach allows us to set the appropriate size for the Jackknife subsamples. The method also provides a means to assess the impact of individual regions on the measured clustering, and thereby to establish whether or not a given galaxy catalogue is dominated by one or several large structures, preventing it to be considered as a fair sample. We apply this methodology to the two- and three-point correlation functions measured from a volume limited sample of M* galaxies drawn from data release seven of the Sloan Digital Sky Survey (SDSS). The frequency of jackknife subsample outliers in the data is shown to be consistent with that seen in large N-body simulations of clustering in the cosmological constant plus cold dark matter cosmology. We also present a comparison of the three-point correlation function in SDSS and 2dFGRS using this approach and find consistent measurements between the two samples.

Clustering of luminous red galaxies II: small scale redshift space distortions

This is the second paper of a series where we study the clustering of LRG galaxies in the latest spectroscopic SDSS data release, DR6, which has 75000 LRG galaxies covering over 1 $Gpc^3/h^3$ for $0.15<z<0.47$. Here we focus on modeling redshift space distortions in $xisp$, the 2-point correlation in separate line-of-sight and perpendicular directions, at small scales and in the line-of-sight. We show that a simple Kaiser model for the anisotropic 2-point correlation function in redshift space, convolved with a distribution of random peculiar velocities with an exponential form, can describe well the correlation of LRG at all scales. We show that to describe with accuracy the so called fingers-of-God (FOG) elongations in the radial direction, it is necessary to model the scale dependence of both bias $b$ and the pairwise rms peculiar velocity $sigma_{12}$ with the distance. We show how both quantities can be inferred from the $xisp$ data. From $r simeq 10$ Mpc/h to $r simeq 1$ Mpc/h, both the bias and $sigma_{12}$ are shown to increase by a factor of two: from $b=2$ to $b=4$ and from $sigma_{12}=400$ to 800 Km/s. The later is in good agreement, within a 5 percent accuracy in the recovered velocities, with direct velocity measurements in dark matter simulations with $Omega_m=0.25$ and $sigma_8$=0.85.

The anisotropic redshift space galaxy correlation function: detection on the BAO Ring

In a series of papers we have recently studied the clustering of LRG galaxies in the latest spectroscopic SDSS data release, which has 75000 LRG galaxies sampling 1.1 Gpc^3/h^3 to z=0.47. Here we focus on detecting a local maxima shaped as a circular ring in the bidimensional galaxy correlation function xi(pi,sigma), separated in perpendicular sigma and line-of-sight pi distances. We find a significant detection of such a peak at r ~110 Mpc/h. The overall shape and location of the ring is consistent with it originating from the recombination-epoch baryon acoustic oscillations (BAO). This agreement provides support for the current understanding of how large scale structure forms in the universe. We study the significance of such feature using large mock galaxy simulations to provide accurate errorbars.

Clustering of luminous red galaxies I: large scale redshift space distortions

This is the first paper of a series where we study the clustering of LRG galaxies in the latest spectroscopic SDSS data release, DR6, which has 75000 LRG galaxies covering over 1 $Gpc^3/h^3$ at $0.15<z<0.47$. Here we focus on modeling redshift space distortions in $xips$, the 2-point correlation in separate line-of-sight and perpendicular directions, on large scales. % and away from the line-of-sight. We use large mock simulations to study the validity of models and errors. We show that errors in the data are dominated by a shot-noise term that is 40% larger than the Poisson error commonly used. We first use the normalized quadrupole for the whole sample (mean z=0.34) to estimate $beta=f(Omega_m)/b=0.34 pm 0.03$, where $f(Omega_m)$ is the linear velocity growth factor and $b$ is the linear bias parameter that relates galaxy to matter fluctuations on large scales. We next use the full $xips$ plane to find $Omega_{0m}= 0.245 pm 0.020$ (h=0.72) and the biased amplitude $b sigma_8 = 1.56 pm 0.09$. For standard gravity, we can combine these measurements to break degeneracies and find $sigma_8=0.85 pm 0.06$, $b=1.85 pm 0.25$ and $f(Omega_m)=0.64 pm 0.09$. We present constraints for modified theories of gravity and find that standard gravity is consistent with data as long as $0.80<sigma_8<0.92$. We also calculate the cross-correlation with WMAP5 and show how both methods to measure the growth history are complementary to constrain non-standard models of gravity. Finally, we show results for different redshift slices, including a prominent BAO peak in the monopole at different redshifts. (Abridged)

Statistical Analysis of Galaxy Surveys - I. Robust error estimation for 2-point clustering statistics

We present a test of different error estimators for 2-point clustering statistics, appropriate for present and future large galaxy redshift surveys. Using an ensemble of very large dark matter LambdaCDM N-body simulations, we compare internal error estimators (jackknife and bootstrap) to external ones (Monte-Carlo realizations). For 3-dimensional clustering statistics, we find that none of the internal error methods investigated are able to reproduce neither accurately nor robustly the errors of external estimators on 1 to 25 Mpc/h scales. The standard bootstrap overestimates the variance of xi(s) by ~40% on all scales probed, but recovers, in a robust fashion, the principal eigenvectors of the underlying covariance matrix. The jackknife returns the correct variance on large scales, but significantly overestimates it on smaller scales. This scale dependence in the jackknife affects the recovered eigenvectors, which tend to disagree on small scales with the external estimates. Our results have important implications for the use of galaxy clustering in placing constraints on cosmological parameters. For example, in a 2-parameter fit to the projected correlation function, we find that the standard bootstrap systematically overestimates the 95% confidence interval, while the jackknife method remains biased, but to a lesser extent. The scatter we find between realizations, for Gaussian statistics, implies that a 2-sigma confidence interval, as inferred from an internal estimator, could correspond in practice to anything from 1-sigma to 3-sigma. Finally, by an oversampling of sub-volumes, it is possible to obtain bootstrap variances and confidence intervals that agree with external error estimates, but it is not clear if this prescription will work for a general case.

Anisotropic Magnification Distortion of the 3D Galaxy Correlation: I. Real Space

It is well known gravitational lensing, mainly via magnification bias, modifies the observed galaxy/quasar clustering. Such discussions have largely focused on the 2D angular correlation. Here and in a companion paper (Paper II) we explore how magnification bias distorts the 3D correlation function and power spectrum, as first considered by Matsubara. The interesting point is: the distortion is anisotropic. Magnification bias preferentially enhances the observed correlation in the line-of-sight (LOS) orientation, especially on large scales. For example at LOS separation of ~100 Mpc/h, where the intrinsic galaxy-galaxy correlation is rather weak, the observed correlation can be enhanced by lensing by a factor of a few, even at a modest redshift of z ~ 0.35. The opportunity: this lensing anisotropy is distinctive, making it possible to separately measure the galaxy-galaxy, galaxy-magnification and magnification-magnification correlations, without measuring galaxy shapes. The anisotropy is distinguishable from the well known distortion due to peculiar motions, as will be discussed in Paper II. The challenge: the magnification distortion of the galaxy correlation must be accounted for in interpreting data as precision improves. For instance, the ~100 Mpc/h baryon acoustic oscillation scale in the correlation function is shifted by up to ~3% in the LOS orientation, and up to ~0.6% in the monopole, depending on the galaxy bias, redshift and number count slope. The corresponding shifts in the inferred Hubble parameter and angular diameter distance, if ignored, could significantly bias measurements of the dark energy equation of state. Lastly, magnification distortion offers a plausible explanation for the well known excess correlations seen in pencil beam surveys.