Correlation function of quasars in real and redshift space from the Sloan Digital Sky Survey Data Release 7
Abstract
We analyze the quasar twopoint correlation function (2pCF) within the redshift interval using a sample of 52303 quasars selected from the recent 7th Data Release of the Sloan Digital Sky Survey. Our approach to 2pCF uses a concept of locally Lorentz (Fermi) frame for determination of the distance between objects and permutation method of the random catalogue generation. Assuming the spatially flat cosmological model with given , we found that the realspace 2pCF is fitted well with the powerlow model within the distance range Mpc with the correlation length Mpc and the slope . The redshiftspace 2pCF is approximated with Mpc and for Mpc, and Mpc and for Mpc. For distances Mpc the parameter describing the largescale infall to density inhomogeneities is with the linear bias that marginally (within 2) agrees with the linear theory of cosmological perturbations. We discuss possibilities to obtain a statistical estimate of the random component of quasars velocities (different from the largescale infall). We note rather slight dependence of quasars velocity dispersion upon the 2pCF parameters in the region Mpc.
keywords:
cosmology: observations – largescale structure, surveys1 Introduction
Spatial distribution of quasars is one of the main sources of observational information about the largescale structure of the Universe; this is at present almost the only source of statistical information about matter inhomogeneity and clustering at cosmological redshifts. The most important results on quasars clustering were obtained using two largest quasar surveys up to date: the 2degree Field Quasar Survey with 2QZ catalogue as a result (http://www.2dfquasar.org, Croom et al. (1998)) and the Sloan Digital Sky Survey (SDSS; http://www.sdss.org), which has been completed with the 7th release (Abazadjian et al., 2009). The twopoint correlation functions (2pCF) (Peebles, 1980) of quasars are important characteristics of matter spatial inhomogeneity that may be compared to cosmological theories of structure formation. Development of 2dF and SDSS surveys gave a powerful incentive to investigation of various aspects of correlation functions of quasars (see, e.g., Croom et al. (2005); Myers et al. (2006) and references therein). In the present paper we use the 7th data release of SDSS to study 2pCF of quasars.
The main problems concerning reconstruction of the cosmological mass distribution from redshift surveys are as follows. The first problem is that the surveys of extragalactic objects give us an information only about distribution of the luminous matter which is biased relative to the dark matter (Dekel & Rees, 1987). Biasing may depend on the physical peculiarities of extragalactic objects, e.g., on morphological type (Einasto et al., 2007; Ross, Brunner & Myers, 2007), luminosity (Beisbart & Kerschner, 2000; Sorrentino, AntonuccioDelogu & Rifatto, 2006), colorindex (Coil et al., 2007), star formation rate (Owers et al., 2007) etc. and evolves with redshift (Croom et al., 2005; Myers et al., 2006, 2007; Weinstein et al., 2004; Porciani, Magliocchetti & Norberg, 2004; da Angela et al., 2008; Mountrichas et al., 2009). Commonly used supposition of linear biasing means that the density variations of certain kind of objects (in our case quasars) is proportional to that of the whole matter; therefore
(1) 
where is 2pCF of the whole matter, is 2pCF of quasars and is the bias parameter.
The second problem of the 3dimensional analysis of the matter distribution is due to the fact that the observed redshifts of extragalactic objects are ‘contaminated’ by measurement errors and nonHubble motions. The distances calculated from these redshifts without taking into account the unknown peculiar velocities are called distances in redshiftspace (in contrast to the real space). As our Universe is isotropic, the correlation function must be spherically symmetric in the realspace. But in the redshiftspace it appears to be distorted. On smaller scales the profile of galaxies 2pCF is stretched along the line of sight (‘Finger of God’ effect) due to virial velocities of objects inside the galaxy clusters (this can be neglected for quasars if we consider ), their random velocities and redshift errors. This effect is especially noticeable for quasar pairs with the projected linear separations Mpc, where the SDSS data are expected to be incomplete (Hennawi et al., 2006). On larger scales the 2pCF profile is flattened along the line of sight (‘The Bull’s eye’ or Kaiser (1987) effect) due to the gravitational infall to density inhomogeneities; this kind of redshiftspace distortion dominates on the linear scales. These effects are parameterized by lineofsight pairwise velocity dispersion and infall parameter . It is worth to note that the nonHubble motions and redshift errors are not the only sources of redshiftspace distortions. One more effect of geometric flattening that can lead to distortion of 2pCF can be due to a wrong choice of cosmological parameters , ; this provides an additional tool for estimation of these parameters by means of a geometrical test of Alcock & Paczynski (1979). The redshiftspace 2pCF of quasars can be in principle used for estimation of all parameters (, and cosmological ones) simultaneously. But due to a degeneracy between the geometric distortions and the redshiftspace distortions (see, e.g., Hoyle et al. (2002)) this problem is complicated. Hoyle et al. (2002) and da Angela et al. (2005) proposed a method that allows to break this degeneracy by means of combination of the Alcock & Paczynski (1979) test with that based on evolution of the quasar clustering amplitude, which has a different dependence on and .
As we cannot estimate the nonHubble motion of each quasar independently, the effects of the redshiftspace distortion are at present the only source of statistical information about proper velocities of quasars. On the other hand, these effects prevent direct determination of 3D 2pCF that involves lineofsight distances between objects determined from measurements. This urges us to use only the projected distances (i.e. orthogonal to the line of sight, which may be considered as independent of proper motions) to determine the projected 2pCF and then to restore the realspace 2pCF. It is well known that such reconstruction of the realspace 2pCF from projected one is mathematically illposed problem. However, most authors avoid this difficulty using a concrete functional form of 2pCF. Typically 2pCF is represented in a powerlaw form , though it is clear that 2pCF slope and correlation length may be different on different interval. For example, da Angela et al. (2005) showed that double powerlow model gives a good fit to 2pCF of quasars distribution with over scales of Mpc and over scales of Mpc for a quasar sample from 2QZ catalogue with redshifts within the range . Using quasar catalogue of spectroscopically confirmed quasars (Schneider, Hall & Richards, 2007) based on the 5th data release of SDSS Ross et al. (2009) obtained Mpc and over scales of Mpc for lowredshift () quasars, and for highredshift () quasars from the same catalogue Shen et al. (2007) obtained Mpc and over scales Mpc. Note that similar large correlation length was obtained for radiosources by Blake, Mauch & Sadler (2003) and for Xray sources by Plionis et al. (2004) and Basilakos et al. (2004); at least a part of both of these object types are considered to be AGNs. In addition to redshift dependence of the quasar clustering confirmed by Croom et al. (2005); Myers et al. (2006, 2007); Porciani, Magliocchetti & Norberg (2004); Porciani & Norberg (2006); da Angela et al. (2008), some authors (Porciani & Norberg, 2006; da Angela et al., 2008; Mountrichas et al., 2009) argue that quasar clustering depends weakly on luminosity. Hence speaking about representation of 2pCF in powerlaw form in the present paper we always mean that 2pCF is averaged over quasars with different luminosities and over certain redshift range and distance intervals described below. In particular, when dealing with the bias estimates we confine ourselves to the intervals where our results may be compared with the linear theory of matter inhomogeneity.
The most complete full investigation on nonHubble motions of quasars together with cosmological parameters was done by 2dF team for 2dF QSO Survey (see Croom et al. (2005); da Angela et al. (2005); Hoyle et al. (2002); Outram et al. (2004)). E.g. da Angela et al. (2005) obtained and (for the sample mean redshift ). In these studies the value of pairwise velocity dispersion of quasars was usually treated as fixed parameter and its value was chosen either from Hubble volume simulations or following previous results on galaxy pairwise velocity dispersion assuming that this value does not evolve strongly with time. For galaxies, we have more definite results for the velocity dispersion obtained from the redshiftspace distortions. The first estimate of Davis & Peebles (1983), km/s, remained a canonical one almost a decade. However, later Mo, Jing & Börner (1993) showed that the value of pairwise peculiar velocity dispersion depends strongly on the sample. Later this result was confirmed by values km/s, given by Marzke et al. (1995) for the Second Center for Astrophysics Galaxy Redshift Survey (CfA2) and the Southern Sky Redshift Survey (SSRS2), and km/s, given by Ratcliffe et al. (1998) for Durham/United Kingdom Schmidt Telescope (UKST) Galaxy Redshift Survey and others. Ratcliffe et al. (1998) also presented the value , which agrees well with later results, such as the results of Hawkins et al. (2003) for 2dF Galaxy Redshift Survey. A similar value, 330 km/s, was obtained for luminous red galaxies (LRG) from the 2dFSDSS LRG and QSO (2SLAQ) survey by Ross et al. (2007).
In this paper, we present new results on the 2pCF slope and correlation length, infall parameter using the last 7th Data Release of SDSS. In Section 2 we describe the quasar samples and technique for random catalogue generation we used in our study. In Sections 3 and 4 we present the results on parameters of the realspace and redshiftspace 2pCFs and techniques for their estimation. Section 4 also contains estimations of the infall parameter and quasar bias from correlation function analysis on linear scales larger than Mpc. In Section 5 we discuss possibilities to use ‘Finger of God’ effect to estimate the random quasars velocities, which does not correlate with density inhomogeneities. Because determination of random pairwise velocity dispersion of quasars needs knowledge of 2pCF in the region Mpc, which is not well studied, we considered different simplified models of 2pCF to study their outcome for . Finally, in Section 6 we sum up the results.
2 Data
2.1 The sample
Our sample is taken from the 7th Data Release of the Sloan Digital Sky Survey (Abazadjian et al., 2009), containing about 100,000 quasars. The redshift range used in our analysis is and as the sky coverage of SDSS contains a one big piece and three narrow nearequatorial ‘stripes’, we excluded these stripes to reduce boundary effects. As the data was taken from the list of all objects classified as quasars, which is not a vetted catalogue, we firstly excluded ‘bad’ objects: objects with failed astrometry and photometry data. The second step was to exclude objects with unreliable redshifts. The redshifts in SDSS are measured mainly using two techniques, emissionline measurements and crosscorrelation, either both or one of these. We excluded objects with the redshift status ‘INCONSISTENT’ and ‘FAILED’. The first status means that the redshift was measured by two techniques but the results were inconsistent. The second one means that the redshift measurement failed. The objects with redshift error were also excluded.
The third step was examination of all objects left after the previous steps by eye (photometry and spectroscopy data) using the SDSS data base. During these examination 466 objects were excluded as the ones having too faint and noisy spectra, or being fault ‘double’ objects (i.e. single objects processed twice and thus appeared twice in the list with similar coordinates), or being artifacts of observations, or being stars (we have found 5 stars in this sample).
Our final sample, which we call full, contains 52303 objects and has the mean redshift . The sky coverage of the sample and its redshift distribution are shown in Figs. 1 and 2 (solid line). The redshiftluminosity distribution of this sample is shown in Fig. 3, where the effect of magnitude limit can be seen. Here the values of absolute magnitudes in band, , are calculated using correction from Richards et al. (2006).
2.2 Sample inhomogeneity
As one can see from Fig. 1 the sky coverage by our sample in not heterogeneous: some parts seem to be significantly denser then others. This is due the fact that these parts were twice spectroscopically covered and can influence our results. First of all we note that the algorithms we use to estimate the 2pCF parameters (sections 3.1 and 4.1) are just worked out in order to reduce the selection effects due to unphysical inhomogeneity (Ivashchenko & Zhdanov, 2010; Zhdanov & Ivashchenko, 2008).
Nevertheless, to test the influence of inhomogeneity we constructed 50 homogeneous samples in the following way. The whole sky area covered by our sample was divided into patches and the density of each patch, mean density ( objects per square degree) and its rms ( objects per square degree) were calculated. All the objects from the patches with density exceeding the mean one more than 0.5rms were excluded (see Fig. 4). Then we generated 50 new homogeneous samples by putting the excluded objects in their places, but choosing (in a random way) only those making the density of the given patch equal to the mean density. The sky coveredge of one of these samples is presented in Fig. 5. Each new test sample has 35643 objects and a density objects per deg. Apparent decrease of the density with increasing declination is the result of projection of the celestial sphere onto the plane.
2.3 Tilable targets
Most quasar candidates in the SDSS are selected based on their locations in multidimentional SDSS colour space and their crossidentification with radio sources from the Faint Images of the Radio Sky at Twentycm (FIRST) survey (Richards et al., 2002). Supplementing this primary quasar sample are quasars targeted by the GALAXY, XRAY, STAR and SERENDIPITY SDSS software packages; no attempt at completeness was made for the last three categories. Thus, the sample we use is not complete and the possible effect of this incompleteness on our results has to be investigated. For this purpose we constructed the sample of the socalled ‘tilable’ objects (uniform) from our sample. According to SDSS glossary webpage, tilable targets are supposed to have as closed to uniform completeness as possible. These are targets with primTarget flag ‘QSO_HIZ’, ‘QSO_CAP’, ‘QSO_SKIRT’, ‘QSO_FIRST_CAP’, ‘QSO_FIRST_SKIRT’, ‘GALAXY_RED’, ‘GALAXY’, ‘GALAXY_BIG’, ‘GALAXY_BRIGHT_CORE’, ‘STAR_BROWN_DWARF’, and second target flag ‘HOT_STD’. The number of objects in the subsample constructed following this criterion is 37290. One can see that the sky coverage by the uniform sample presented in Fig. 6 is close to uniform and close to the sky coverage by our homogeneous samples (Fig. 5). The redshift distribution of this uniform sample is presented in Fig. 2 (dashed line).
2.4 Reddening
As the main criterion of quasar candidates selection in SDSS is based on their magnitudes and colors, that have been corrected for Galactic extinction using the maps of Schlegel, Finkbeiner & Davis (1998), any systematic errors in the reddening model can induce additional effects on clustering results (see Ross et al. (2009) for discussion of possible effects). That is why following Ross et al. (2009) we devided our sample into low reddening () sample and high reddening () one. Values of were obtained from the reddening in each SDSS band given in the catalogue and using the relative extinction for different bandpasses from Table 6 in Schlegel, Finkbeiner & Davis (1998). The numbers of objects in subsamples are 30188 and 22115 for low and high reddening correspondingly. The redshift distributions of these subsamples are presented in Fig. 2 (dashdot for lowreddening sample and dashdotdotdot for highreddening one).
2.5 Bad fields
Since the quasar selection algorithm involves colours and magnitude all photometric errors caused by e.g. photometric callibration can induce errors to this selection and result in additional inhomogeneity of the sample. For verification of these effect we constructed a good sample using information about the data quality for each field (an image in all five bands with approximate dimentions of ) from a field quality table, called RunQA. To construct the good sample we used the FieldQAll number which describes the overall quality of a field and rejected the objects with its value equal to ‘bad’. This overall quality determination, FieldQAll, is based on the (dereddened) principal colors of all the stars in the field vs. the Galactic value of this quantity, the PSF quality, the difference between Aperture and PSF magnitudes for the same stars in the field, and the seeing. We also excluded the objects with rerun number equal to zero, as there are no fields with such rerun number in FieldQAll table. The number of objects in the good sample is 50420. Its redshift distribution is presented in Fig. 2 (long dashes).
2.6 Fibre collisions
One another peculiarity of our sample is related to the so called ‘fibre collisions’ effect, which is inherent for multifibre spectrographs. It means that two spectroscopic targets on the same plate cannot be measured simultaneously if they are closer than some distance ( for SDSS); it is only possible to obtain spectra for both such objects if they are covered by more than one plate (Schneider, Hall & Richards, 2007). Angular separation corresponds to Mpc projected distance for our redshift interval . That is why we avoided this distance range when calculating the parameters of 2pCF. Of course this effect affects mostly the precision of our random velocities dispersion estimations due to the small sample of pairs with projected separation Mpc. But in this case we cannot follow e. g. Ross et al. (2009) supplementing their sample with quasar candidates from SDSS Quasar photometric catalogue with photometric redshifts (Richards et al., 2009) due to the fact that the errors of photometric redshifts are about which is much larger than the restrictions we put on our sample.
2.7 Random samples
An important point in reconstruction of the correlation function of extragalactic objects is generation of the random sample, which should reproduce a random distribution of objects with the same redshifts and angular distributions inherent to the initial sample as much as possible. For this purpose an imaging mask may be used (see e.g. Myers et al. (2006)). In the present work we used the permutation method of random catalogue generation (see e.g. Osmer (1981); Zhdanov & Surdej (2001); Croom et al. (2005)) that does not require any technical information about the survey. This method uses the random permutation of objects’ redshifts from the initial list of quasars while keeping their coordinates unchanged. As a result one obtains a randomized list of objects with the same redshift distribution and the same distribution over the sky as in the initial sample. At the same time this randomization destroys correlations between objects of the initial sample: those objects which were clustered in the initial sample appear to be not clustered in the randomized sample as they have most probably different redshifts and thus are located far from each other.
3 Realspace correlation function
3.1 The technique
As peculiar velocities of objects distort the distance between objects along the line of sight, the parameters of the realspace 2pCF of quasars can be reconstructed from the projected 2pCF. We use this standard approach, but taking into account some specific features of the permutation method of random catalogue generation.
2pCF determines the probability to find a neighbour of the object from the sample at comoving distance in a volume element as (Peebles, 1980)
(2) 
where is a number density of th object in the neighbourhood of . Here we must specify the definition of the distance between neighbours, which can be defined differently in curved spacetime. For the distance scales under consideration this may be done with a sufficient accuracy using local Lorentz frames (Misner, Thorne & Wheeler, 1973). For definiteness we use the Fermi reference frames and normal Fermi coordinates (Synge 1960; Misner et al. 1973); this is especially convenient because in case of spatially flat cosmological background this leads to distance relations analogous to that of Euclidean space (see Appendix).
We suppose that the realspace 2pCF of quasars has powerlow form for distances and , where is some characteristic comoving scale where 2pCF approaches to zero. We also make the following assumptions. First, we neglect variations of on the scales . In fact, we also considered models that take into account variations of with redshift on this scale and found the resulting difference to be not statistically significant, as it can be expected from very beginning. Second, we assume, as usual, that quasars distribution is locally isotropic at any point of the realspace, i.e. it depends only on the distance between the objects. Third, we neglect variations of redshift distribution of the sample in different sky regions.
Further we pass to projected (comoving) separations. The only thing that makes our consideration different from the standard one (Peebles, 1980) is related to our method of the random catalogue generation. When counting the number of pairs with projected separations within the range , we include only pairs with lineofsight separations that satisfy the condition (Zhdanov & Surdej, 2001; Zhdanov & Ivashchenko, 2008). (Here we denote the distance along the line of sight as to distinguish it from the number .) This number of pairs should be compared to analogous number in the randomized sample. If, e.g., we do not restrict and include all pairs with the above projected separation, there obviously would be no difference in these numbers between the initial sample and the randomized one obtained by permutation of redshifts. On the other hand, we expect that almost all the pairs that have a noticeable correlation are concentrated within the region . Note that though at this point we use the observational values of that involve contribution of the proper velocities, corresponding correction due to velocities here is negligible, because it affects only larger distances .
The total number of neighbouring quasars with the projected distances from a given quasar within the range and lineofsight separations is obtained by summing over all quasars from the sample
(3) 
where the sum is performed over all objects in the sample, . This number is twice larger than the number of pairs with separation from the above range, if the quasars and the neighbours are chosen from the same sample. In order to study boundary effects we also considered the case when the quasars are chosen from a subsample obtained from the initial sample with the exception of objects near the boundary, and the neighbours being looked for in the whole initial sample.
In case of the randomized redshifts the corresponding number of neighbours is
where is the number density in the neighbourhood of th object in the randomized sample. As the redshift distribution of objects in initial and random samples are the same, than one estimates with a relative error , where is the number of objects in the sample. Therefore
(4) 
where we extended the integration limits to infinity.
Note that possible difference of the ratio from unity may be due to variations of distribution of quasars in the different sky regions of the real survey. Corresponding correction was analyzed and found to be statistically insignificant (Zhdanov & Ivashchenko, 2008).
In fact we used not exactly Eq. (4) but corresponding relations for separations in finite intervals ; typically we used Mpc. Then using the powerlow form of 2pCF we have the following relations for fitting parameters and (cf., e.g., Shanks et al. (1980); Collins, HeydonDumbleton & MacGillivray (1989)):
(5) 
where
The errors of all the parameters in the present study have been estimated using the ‘jackknife’ method (see, e.g., Miller (1974) for review), which has been successfully applied to clustering problems in extragalactic astronomy (see, e.g., Myers et al. (2006)). We divided our initial sample into equal stripes along declination; owing to large number of objects in each stripe, the stripes may be viewed as practically independent ones (here we neglect the boundary effects). Then new subsamples have been formed: th subsample is unification of all the stripes except th stripe. Jackknife dispersion estimate of any parameter calculated in this treatment is then
(6) 
where is the corresponding value of this parameter for th subsample and is either the value of for the initial sample or it is average over subsamples.
3.2 Results
We determine the distances between quasars using the locally Lorentz frame of the first quasar (see Appendix). This method differs, e.g., from that by Hoyle et al. (2002); da Angela et al. (2005).
dist. range, Mpc  range  mean  , h Mpc  sample  authors  
1.054  full  
2.084  test  
df  2.084  test  
–  homogeneous  
1.47  0.789  uniform  present paper  
1.231  good  
0.846  low reddening  
1.025  high reddening  
1.47  –  2QZ  Porciani et al. (2004)  
1.4  –  2QZ  da Angela et al. (2005)  
1.4  –  2QZ  da Angela et al. (2005)  
–  –  SDSS DR5Q  Shen et al. (2007)  
1.27  –  SDSS DR5Q  Ross et al. (2009) 
Our sample covers only a part of the sky; thus the boundary effects can influence the result because the number of neighbours for objects lying close to the boundary is less then for those lying far from it. To estimate this influence we made calculations in two different ways. In the first one we calculated the number of neighbours for all the quasars in the sample. In the second one we took into account the neighbours of the objects lying at least 5 degrees far from the boundary (we called this test sample). Fig. 7 shows the surface number density of neighbours in the initial sample (circles) and in the random one (triangles) correspondingly, where we choose the bin size . Values are the mean numbers over 50 realizations of random sample, error bars represent 1 jackknife errors.
As can be seen from Fig. 7, in practice the surface density of neighbours for the random catalogue does not depend on the distance (except for the distance range Mpc), which is expected for a random distribution of objects. The gap at Mpc can be explained by a deficit of close pairs in the initial catalogue. Note that angular distribution of objects in our random samples is the same as in the initial one (see Sec. 2.7). This deficit is a result of the fibre collision effect (Sec. 2.6). One can see that the number densities for initial and random samples coinside up to Mpc, thus we use Mpc (or Mpc). And we also put Mpc because of the fiber collision effect. The quasar projected 2pCF divided on in logarithmic scale for the full sample is presented in Fig. 8 with the the best fit single powerlow over the range Mpc.
The same calculations were provided for 50 test samples (see Sec. 2.2). The scattering of the resulting parameters is presented in Fig. 9 with the mean values presented in the second row of Table 1, where errors are rms. Results for homogeneous, uniform, good, low reddening and high reddening samples are presented in rows of this table. For these samples the errors are jackknifes. Comparison of the projected 2pCF for the initial sample with uniform, good, low reddening and high reddening samples are presented in Figs. 1012.
3.3 Discussion
As one can see from Figs. 1012 and Table 1, the projected 2pCFs and parameters of the realspace 2pCF for all samples agree well with each other within the errors that are slightly larger for the samples with smaller number of objects. Hence, firstly, we can neglect boundary effects for our sample, that is expected due to the sample size being several times larger than the distance range we consider. Secondly, we can also neglect all the inhomogeneities of the sample caused by the nonuniformity of the sky coverage as a result of the spectroscopic survey, reddening and the conditions of observations. The last five rows of Table 1 show comparisons of our results with the results of other authors. We can see that our results agree well with early results for the 2QZ catalogue by Porciani, Magliocchetti & Norberg (2004) within 1 for and 2 for . The best agreement is with the results by Ross et al. (2009) for the 4th edition of SDSS quasar catalogue (compiled from the 5th data release of SDSS by Schneider, Hall & Richards (2007)), but note that they used the wider distance range. The 2pCF parameters from Shen et al. (2007) are included to the table only for completeness of known results on the quasars 2pCF in the realspace as they correspond to larger redshifts. Our results cannot be directly compared with that by da Angela et al. (2005) for 2QZ catalogue, though there is rather a good agreement if we average their values over both distance ranges and Mpc. Note that our redshift range is smaller, but the fraction of quasars contained within the range is small compared with those in the range. Finally, we note that comparison of 2pCF parameters on different distance intervals shows that the powerlaw form of 2pCF must be replaced by more complicated one for larger distances.
4 RedshiftSpace correlation function and infall parameter
4.1 The technique
2pCF in the redshiftspace is not sphericallysymmetric because of the effects of velocities mentioned in Sec. 1. Here we pay attention mainly to the scales larger than 10 Mpc, where we expect the density contrast to be sufficiently small. In this region the effect of bulk gravitational infall leads to flattening of redshiftspace 2pCF contours and to certain enhancement of the monopole part of this 2pCF, which is obtained by averaging over all angles. On the other hand, in this region we neglect ‘Finger of God’ effect due to random velocities.
The relation of the realspace 2pCF to the redshiftspace 2pCF is well known within the linear theory of cosmological perturbations (Kaiser, 1987; Hamilton, 1992; Matsubara & Suto, 1996). The monopole (i. e. angle averaged) part of 2pCF is proportional to the realspace 2pCF (Kaiser, 1987) and therefore it also has a powerlaw form with the same slope , but with different correlation length:
(7) 
This formula is wellgrounded in the region Mpc, where the linear theory is applicable. Here the parameter appears as some consequence of the largescale infall onto matter overdensities, is the bias parameter in Eq. (1). The function is expressed by means of the growing mode of the density contrast (Peebles, 1980). There is a simple analytical approximation (Carroll & Press, 1992) for this function in the case of the spatiallyflat CDM cosmological model:
(8) 
Therefore we may proceed in two ways: either (i) by taking the slope from the results of Sec. 3 dealing with fitting of only, or (ii) by independent determination of .
Like in Sec. 3 one can determine the total number of neighbours of the th quasar from the whole sample in a spherical layer as
where . The similar estimation for random catalogue, which is considered to represent random spatial distribution of objects with no clustering, is
Assuming and taking into account the powerlow form of the monopole part of 2pCF , we have the relation for fitting:
(9) 
In order to determine the infall parameter on the way (i) we use Eq. (9) to find , with the slope being known from the results of the Sec. 3. Then we derive by solving the Eq. (7) with known , and then we estimate bias for our redshift interval using Eq. (8). These operations deal with the region Mpc. In addition, mainly for comparison with the results of the other authors, we also considered way (ii) with independent determination of and for different separation intervals.
Note that the method we use for estimation of the infall parameter is similar to the socalled method (see e.g. Ratcliffe et al. (1998)), which uses the volume integrals . Using of these integrals allows to smooth the influence of small sizes of the samples. Our approach has the same advantage as it includes the 2pCF parameters which are the averaged characteristics of 2pCFs within the whole distance range.
4.2 Results and discussion
We estimated the 2pCF parameters within the whole range Mpc and separately for ranges Mpc and Mpc. The results are presented in Table 2 and Fig. 13. The errors are calculated using the jackknife technique. We also found the same parameters for distance range Mpc for comparison with the results of Croom et al. (2005) and Ross et al. (2009) presented in the last 3 rows of this table. As one can see our parameters and within the distance range Mpc agree well with the results of Croom et al. (2005) and Ross et al. (2009) for 2QZ catalogue and the 4th edition of the SDSS quasar catalogue (Schneider, Hall & Richards, 2007) correspondingly. And for distance range Mpc our results agree with Croom et al. (2005) within 2 for the slope and 1 for the correlation length.
dist. range, Mpc  range  mean  , h Mpc  sample  authors  

1.440  initial  
–  homogeneous  
1.473  uniform  
1.47  1.363  good  present paper  
0.779  low reddening  
1.892  high reddening  
1.060  
1.47  0.456  initial  present paper  
1.636  
1.35  –  2QZ  Croom et al. (2005)  
1.35  –  2QZ  Croom et al. (2005)  
1.27  –  SDSS DR5Q  Ross et al. (2009) 
mean  method  sample  authors  
0.55  mean of [1],[2]  2SLAQ+2QZ+SDSS DR5  Mountrichas et al. (2009)  
1.27  [3]  SDSS DR5Q  Ross et al. (2009)  
1.35  [3]  2QZ  Croom et al. (2005)  
1.40  [2]  2QZ+2SLAQ  da Angela et al. (2008)  
1.40  2QZ  da Angela et al. (2005)  
1.40  –  [2]  2QZ  da Angela et al. (2005)  
1.40  –  [3]  SDSS DR4phot  Myers et al. (2007)  
1.47  –  [4]  2QZ  Porciani, Magliocchetti & Norberg (2004)  
initial  
homogeneous  
1.47  decribed in text  uniform  present paper  
good  
low reddening  
high reddening 
Comparing the parameters for different distance ranges we can conclude that on the large scales the redshiftspace 2pCF is steeper than on the small scales while the correlation lengthes do not differ significantly within 1. This difference in the slope of the redshiftspace 2pCF of quasars on different scales was pointed out earlier (Croom et al., 2005). It is important to note that the results for the slope on Mpc interval fairly well agree with the results for in the real space, as it must be.
As discussed in Section 4.1, we have used the correlation length and the slope of the realspace 2pCF obtained in Section 3. We have then obtained the redshiftspace correlation length Mpc for the fixed on the scales Mpc and we have found and using Eqs. (7) and (8). The results are given in Table 3 that contains the obtained values of infall and bias parameters and results of other authors with indication of the method, used for their estimation. One of them (denoted as [1]), containing simultaneous determination of and cosmological model, was proposed by Hoyle et al. (2002) and includes fitting of the observed by theoretical model from Peebles (1980) (see also Ratcliffe et al. (1998)) with the relation for infall velocity from HaleSutton (1990). The improvement of this method proposed in da Angela et al. (2005) (denoted as [2]) allows to break the degeneracy between and and is based in the consideration of the quasar clustering evolution. Another approach proposed in Croom et al. (2005) (denoted as [3]) suggests the solution of the quadratic equation in with the shape of governed by the underlying dark matter distribution described by the analytic equation from Hamilton et al. (1991, 1995). In da Angela et al. (2008) and Ross et al. (2009) the last approach was applied with theoretical dark matter power spectrum from Smith et al. (2003). One more method which was used in Porciani, Magliocchetti & Norberg (2004) (denoted as [4]) includes using the relation between projected 2pCFs of the quasars and that of the matter. As distinct from a number of authors (aimed, in particular, at estimations of the dark matter halo masses) we concentrate to the comparison with the linear theory; i.e. our results on are completely based on the interval Mpc, where there is no problems with the linear theory.
For test we also provided the same calculations for redshiftspace 2pCF, and for homogeneous, uniform, good, low reddening and high reddening samples. The scattering of the resulting redshiftspace 2pCF parameters and infall parameters with the corresponding values of for 50 homogeneous samples are presented in Figs. 14 and 15. Their mean values are presented in Tables 2 and 3 together with results for other samples. For all our samples errors of the parameters are jackknifes, except the homogeneous where errors are rms. All these results agrees well (within 1.3) with the results for initial sample.
5 Random velocities dispersion
5.1 Technique
In this Section we discuss the ‘Finger of God’ effect as a tool to estimate the random component of the peculiar velocity of quasars with respect to the CMB frame, i.e. the dispersion . We neglect virialization and assume that the random velocities do not correlate with mass density variations (contrary to Kaiser (1987) effect); such assumption is quite reasonable in case of quasars. If the quasars in a pair can be considered as independent, the lineofsight pairwise velocity dispersion is , where onedimensional pairwise velocity is the difference of the peculiar velocities (above the Hubble flow) , of quasars in a pair: . The effect of random velocities typically dominates at smaller separations, where the gravitational infall is less pronounced. Here, we neglect the effect of the gravitational infall.
There are two major obstacles. Firstly, as it was mentioned above, the SDSS data are incomplete in the region of linear comoving projected separations Mpc. Nevertheless one can try to use an information from this region by testing different distributions. Here we assume that the selection effects due to the fibre collisions do not distort considerably the distribution of quasar pairs from our sample with projected comoving separations Mpc and the redshift . Another obstacle is the poor knowledge of the 2pCF form in this region. This is the main point of this Section: to study the model dependence of the results for . In order to obtain , we need 2pCF for all distances. Therefore, we used the powerlaw form of the realspace 2pCF within the comoving distance range Mpc with the parameters defined in Sec. 3 and assumed a rough estimate for :
(10) 
According to Sec. 3 we used , Mpc and , will be varied ( to avoid divergence for ). In case of 2pCF is represented by a continuous function. On account of Hennawi et al. (2006) results, the rough choice , may be considered as plausible. However, even rather wild () choices do not lead to a drastic differences of results. Namely, we saw that the results are rather stable (taking into account errors) under different choices of , .
Formula (10) allows to estimate the number of pairs with separation . Then we restricted ourselves with the region and use relation
to derive distribution of pair numbers (minus the random background) in lineofsight separations .
We considered two kinds of the pairwise velocity distribution: Gaussian and exponential (Ratcliffe et al., 1998; Hoyle et al., 2002; Hawkins et al., 2003). The pairwise velocity contributes to the change of the comoving lineofsight distance between the quasars in a pair as follows: , where with . This coefficient grows from 0.85 to 0.97 within the interval , hence we approximate it as .
First we assumed the exponential pairwise velocity distribution (Ratcliffe et al., 1998). Then we have probability density of random variable as follows
(11) 
In case of Gaussian distribution of velocities we have
(12) 
The convolution of with yields the distribution of pair numbers as a function of . We used this to fit the observed pair numbers in different bins with the projected separation Mpc after randomized background subtraction. The randomized background is calculated as described in Sec. 2.7.
5.2 Results and discussion
Because the number of pairs is not large, there is a considerable variation of results for different binnings, so we do not consider the numerical estimates as very reliable. However, we may argue that the model dependence on and parameters is rather insignificant. For illustration, we present the results (Table 4) of fitting the pair numbers (after randomized background subtraction) on account of (11) and (Table 5) on account of (12) for two different binnings having minimal jackknife dispersion estimates for . Binning 1 involves numbers of pairs with lineofsight separations from five intervals , , , , (Mpc). Binning 2 deals with separations from three intervals , , (Mpc). The bar charts in Fig. 16 and 17 we show corresponding numbers after the average background subtraction. Here we prefer to avoid separations less than 2 Mpc. Also for illustration in Fig. 16 and 17 we present the results of fitting for and different values of . As one can see from the results, the values of do not vary considerably under rather different (even unrealistic) choices of unknown parameters , : the variations are of the order of the jackknife dispersion estimate.
620140  660150  740160  
180  580190  
720160  740170  750170  
550190  580190  590190  
750165  750165  750165  
590190  590190  590190 
58090  60090  67090  
125  600110  
650100  670100  69080  
560125  590130  600120  
67090  67090  69095  
580130  600100  600120 
In fact the values in Tables 5, 6 do not represent real quasar velocity dispersion, but a superposition of two components: , where is the true velocity dispersion and is due to the redshift measurement errors. The last value corresponds to intrinsic emission lines shifts in quasars pointed out by Croom et al. (2005), as it is an estimation from emissionline or crosscorrelation technique of redshift measurement caused by the impossibility of precise determination of the line centre due to its large width. For detailed investigation of this effect see Richards et al. (2002); Shen et al. (2007). The redshift measurement errors from the SDSS data base for our sample is , which corresponds to km/s. In this case we would obtain too large value of .
Once work on the first version of the text had been finished, the paper by Hewett & Wild (2010) was appeared with the error estimate km/s for SDSS data. In this case errors dominate in . If we take the results from Table 5, which yield smaller jackknife dispersion estimate for , than the residual upper estimate of km/s seems to be more adequate.
The values of the pairwise velocity dispersion corresponding to Tables 4, 5 are comparable with that mentioned in previous studies of quasar clustering. Outram et al. (2001) used the value 400 km/s obtained from the Hubble Volume Simulations; da Angela et al. (2005) fixed the value 800 km/s and found that this value gives an adequate result for Mpc region and noted that it is dominated by the rms pairwise redshift error km/s; the same value 800 km/s was used by da Angela et al. (2008) and the close value 690 km/s was used by Croom et al. (2005). The last value was chosen as a mean one of the range 630750 km/s; this is a quadrature superposition of (i) the real pairwise velocity dispersion recalculated with redshift estimated from the galaxies pairwise velocity dispersion 500 km/s (Hawkins et al., 2003), (ii) redshift measurement error obtained from repeat observations, and (iii) velocity error due to intrinsic emission lines shifts in QSOs (Richards et al., 2002). Only Mountrichas et al. (2009) treating the velocity dispersion as a free parameter estimated its value for a sample of quasars and luminous red galaxies (LRG) as 620 km/s and 727 km/s (two values are the result of two different approaches) and noted that they are consistent with the value 728 km/s which is a quadrature superposition of the values 800 km/s for quasars and 300 km/s for LRG from previous studies.
6 Conclusions
We have performed a detailed analysis of 2pCF using the quasar sample from the SDSS DR7 data for objects with . Our main conclusions are as follows.

Using our approach for estimation of the realspace 2pCF of quasars we confirm that it is fitted well with the powerlow model in the distance range Mpc with the correlation length of Mpc and the powerlow slope . The results are in a good agreement with the results of the other authors on the same distance interval. However, this comparison suggests that the powerlow model seems to be not valid on larger intervals.

The redshiftspace 2pCF of quasars has a break and it may be approximated well with the double powerlow model with Mpc and for Mpc and Mpc and on scales Mpc with the mean values Mpc and for both distance ranges together.

For distances Mpc the parameter describing the largescale infall to density inhomogeneities is with the linear bias to be that marginally (within 2) agrees with the linear theory of cosmological perturbations.

We point out that the value of the quasars velocity dispersion varies almost within errors for a wide range of model parameters that determine 2pCF for Mpc.
We also performed investigations of our sample and showed that we can neglect boundary effects and all inhomogeneities of the sample caused by nonuniformity of the sky coverage by the spectroscopic survey, reddening and conditions of observations.
It is worth to note, that the values of from Tables 4, 5 have been calculated on account of SDSS DR7 data on pairs with the projected separations Mpc, which are not complete and the number of pairs involved in treatment of the velocity dispersion is not too large. In view of these circumstances we consider our results on rather as motivation to further studies to observe larger sample of Mpc quasar pairs with accurate redshift measurements.
Acknowledgments
This work has been partially supported by Swiss National Science Foundation (SCOPES grant No 128040).
The authors are thankful to the Sloan Digital Sky Survey team. Funding for the SDSS has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Aeronautics and Space Administration, the National Science Foundation, the US Department of Energy, the Japanese Monbukagakusho, and the Max Planck Society. The SDSS Web site is http://www.sdss.org.
The authors are also thankful to the anounimous referee for helpful comments and suggestions.
Appendix A Distances determination
We consider the spatially flat Universe with FriedmannRobertsonWalker (FRW) metric
(13) 
where and are right ascention and declination respectively, at the present epoch . In case of the flat CDM cosmological model () for an object with redshift we have
The distances in the redshiftspace are defined as if the objects were at rest with respect to the CMB. In order to define the distance from a quasar with coordinates to a neighbouring object with coordinates we proceed as follows (Ivashchenko & Zhdanov, 2010).
(a) Consider the Fermi frame (Synge 1960; Misner et al. 1973) of the quasar with the fiducial trajectory that may be interpreted as the trajectory of an observer, which is at rest with respect to the CMB.
(b) Introduce the hypersurface of constant time in corresponding normal Fermi coordinates with the origin . In our case this is the spacelike hypersurface , which is orthogonal to the fiducial trajectory at the origin.
(c) Find the invariant (proper) distance from the origin to the intersection of the neighbour world line with the hypersurface; this is uniquely determined along the spacelike geodesic connecting two points. In a general case we must suppose that the objects are not too far from each other to avoid the conjugate points.
In our case is easily determined along the geodesic on the hypersurface between the points and ; then the comoving distance is:
(14) 
where is the angle between objects on the sky. The comoving linear projected distance from to is
(15) 
and the comoving lineofsight distance from to is
(16) 
The above formulae determine distances in the redshiftspace. Transition to the realspace involves the velocities of quasars with respect to cosmological background. However, in case of the projected distances the effect of velocities is negligible and these distances obtained from observations may be used to derive parameters of the realspace 2pCF.
Footnotes
 pagerange: Correlation function of quasars in real and redshift space from Sloan Digital Sky Survey Data Release 7–A
 pubyear: 2010
References
 Abazajian K. N., AdelmanMcCarthy J. K., Agüeros M. A. et al., 2009, ApJS, 182, 543
 Alcock C., Paczynski B, 1979, Nature, 281, 358
 Basilakos S., Plionis M., Georgakakis A., Georgantopoulos I., 2005, MNRAS, 356, 183
 Beisbart C., Kerschner M., 2000, ApJ, 545, 6
 Blake C., Mauch T., Sadler E. M., 2004, MNRAS, 347, 787
 Carroll S. M., Press W. H., 1992, Annual Rev. of A&A, 30, 499
 Coil A. L., Newman J. A., Croton D. et al., 2008, ApJ, 672, 153
 Collins C. A., HeydonDumbleton N. H., MacGillivray H. T., MNRAS, 1989, 236, 7
 Croom S. M., Miller L., Shanks T. et al., AAO Newsletters, 1998, No. 86, 7
 Croom S. M., Boyle B. J., Shanks T. et al., 2005, MNRAS, 356, 415
 da Ângela J., Outram P. J., Shanks T., Boyle B. J., Croom S. M., Loaring N. S., Miller L., Smith R. J., 2005, MNRAS, 360, 1040
 da Ângela J., Shanks T., Croom S. M. et al., 2008, MNRAS, 383, 565
 Davis M. & Peebles P. J. E., 1983, ApJ, 267, 465
 Dekel A., Rees M. J., 1987, Nat, 326, 455
 Einasto M., Saar E., Einasto J. et al., 2007, preprint (0706.1126)
 HaleSutton D., 1992, PhD Thesis, University of Durham
 Hamilton A. J. S., Kumar P., Lu E., Mattheus A., 1991, ApJ, 374, L1
 Hamilton A. J. S., 1992, ApJ, 385, L5
 Hamilton A. J. S., Kumar P., Lu E., Mattheus A., 1991, ApJ, 442, L73
 Hawkins E., Madox S., Cole S. et al., 2003, MNRAS, 346, 78
 Hennawi J. F., Strauss M. A., Oguri M. et al., 2006, AJ, 131, 1
 Hewett P. S., Wild V., 2002, preprint (1003.3017)
 Hoyle F., Outram P. J., Shanks T., Boyle B. J., Croom S. M., Smith R. J., 2002, MNRAS, 332, 311
 Ivashchenko G. Yu., Zhdanov V. I., 2010, Kinematika i Fizika Nebesnykh Tel, 26, 43
 Kaiser N., MNRAS, 1987, 227, 1
 Komatsu E., Dunkley J., Nolta M. R. et al., 2009, ApJS, 180, 330
 Larson D., Dunkley J., Hinshaw G. et al., 2010, preprint (1001.4635)
 Marzke R. O., Geller M. J., da Costa L. N., Huchra J. P., AJ, 1995, 100, 477
 Matsubara T., Suto Ya., ApJ, 1996, 470, L1
 Miller R. G., 1974, Biometrica, 61, 1
 Misner J. A., Thorne K. S., Wheeler J. A. 1973, Gravitation. W.H.Freeman & Co, San Francsko
 Mo H. J., Jing Y. P., Börner J., 1993, MNRAS, 264, 825
 Mountrichas G., Sawangwit U., Shanks T., Croom S. M., Schneider D. P., Myers A. D., Pimbblet K., MNRAS, 2009, 394, 2050
 Myers A. D., Brunner R. J., Richards G. T. et al., 2006, ApJ, 638, 622
 Myers A. D., Brunner R. J., Nichol R. C., Richards G. T., Schneider D. P., Bahcall N. A., 2007, ApJ, 658, 85
 Osmer P. S., 1981, ApJ, 247, 762
 Outram J., Hoyle F., Shanks T., Boyle B. J., Croom S. M., Loaring N. S., Miller L., Smith R. J., 2001, MNRAS, 328, 174
 Outram J., Shanks T., Boyle B. J. et al., 2004, MNRAS, 348, 745
 Owers M. S., Blake C., Couch W. J., Pracy M. B., Bekki K., 2007, MNRAS, 481, 494
 Peebles P. J. E., 1980, The LargeScale Structure of the Universe, Princeton Univ. Press, Princeton, NJ
 Plionis M., Basilakos S., Georgakakis A., Georgantopoulos I., 2004, in Majica R., Maiolino R., eds, Proc. of the Guillermo Haro Conf., World Scientific, Singapore, p. 425
 Porciani C., Magliocchetti M., Norberg P., 2004, MNRAS, 355, 1010
 Porciani C., Norberg P., MNRAS, 2006, 317, 1824
 Ratcliffe A., Shanks T., Parker Q. A., Fong R., 1998, MNRAS, 296, 191
 Richards G. T., Vanden Berk D. E., Reichard T. A., Hall P. B., Schneider D. P., SubbaRao M., Thakar A. R., York D. G., 2002, AJ, 124, 1
 Richards G. T., Strauss M. A.., Fan. X. et al. 2006, AJ, 131, 2766
 Richards G. T., Myers A. D., Gray A.G. et al. 2009, ApJS, 180, 67
 Ross A. J., Brunner R. J., Myers A. D., 2007, ApJ, 665, 67
 Ross N. P., da Angela J., Shanks T. et al, 2007, MNRAS, 381, 573
 Ross N. P., Shen Y., Strauss M. A. et al., 2009, ApJ, 697, 1634
 Shen Y., Strauss M. A., Oguri M. et al., 2007, AJ, 133, 2222
 Schlegel D. J., Finkbeiner D.P., Davis M., 1998, ApJ, 500, 525
 Schneider D. P., Hall P. B., Richards G. T., 2007, AJ, 134, 102
 Shanks T., Fong R., Ellis R. S., MacGillivray H. T., 1980, MNRAS, 192, 209
 Smith R. E., Peacock J. A., Jenkins A., White S. D. M., Frenk C. S., Pearce F. R., Thomas P. A., Efstathiou G., Couchman H. M. P., 2003, MNRAS, 341, 1311
 Sorrentino G., AntonuccioDelogu V., Rifatto A., 2006, in Combes F., Palous J., eds, Proc. of the IAU Symposium N. 235, Galaxy Evolution Across the Hubble Time, p. 243
 Synge J. L., 1960, Relativity: the General Theory. NorthHolland, Amsterdam
 Weinstein M. A., Richards G. T., Schneider D. P. et al., 2004, ApJS, 155, 243
 Zehavi I., Blanton M. R., Frieman J. A. et al., 2002, ApJ, 571, 172
 Zhdanov V. I., Ivashchenko G. Yu., 2008, Kinematika i Fizika Nebesnykh Tel, 24, 1
 Zhdanov V. I., Surdej J., 2001, A&A, 372, 1