Issue 
A&A
Volume 608, December 2017



Article Number  A44  
Number of page(s)  21  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201630276  
Published online  05 December 2017 
The VIMOS Public Extragalactic Redshift Survey (VIPERS)
Gravity test from the combination of redshiftspace distortions and galaxygalaxy lensing at 0.5 < z < 1.2^{⋆}
^{1} AixMarseille Univ, CNRS, LAM, Laboratoire d’Astrophysique de Marseille, UMR , 7326 Aix Marseille, France
email: sylvain.delatorre@lam.fr
^{2} INAF–Osservatorio Astronomico di Brera, Via Brera 28, 20122 Milano – via E. Bianchi 46, 23807 Merate, Italy
^{3} Dipartimento di Fisica, Università di MilanoBicocca, P.zza della Scienza 3, 20126 Milano, Italy
^{4} AixMarseille Univ, Univ. Toulon, CNRS, CPT, 13284 Marseille, France
^{5} Università degli Studi di Milano, via G. Celoria 16, 20133 Milano, Italy
^{6} INAF–Istituto di Astrofisica Spaziale e Fisica Cosmica Milano, via Bassini 15, 20133 Milano, Italy
^{7} INAF–Osservatorio Astronomico di Bologna, via Gobetti 93/3, 40129 Bologna, Italy
^{8} INAF–Osservatorio Astrofisico di Torino, 10025 Pino Torinese, Italy
^{9} Laboratoire Lagrange, UMR 7293, Université de Nice Sophia Antipolis, CNRS, Observatoire de la Côte d’Azur, 06300 Nice, France
^{10} Dipartimento di Fisica e Astronomia – Alma Mater Studiorum Università di Bologna, viale Berti Pichat 6/2, 40127 Bologna, Italy
^{11} Institute of Physics, Jan Kochanowski University, ul. Swietokrzyska 15, 25406 Kielce, Poland
^{12} National Centre for Nuclear Research, ul. Hoza 69, 00681 Warszawa, Poland
^{13} INFN, Sezione di Bologna, viale Berti Pichat 6/2, 40127 Bologna, Italy
^{14} AixMarseille Université, Jardin du Pharo, 58 bd Charles Livon, 13284 Marseille Cedex 7, France
^{15} IRAP, 9 Av. du colonel Roche, BP 44346, 31028 Toulouse Cedex 4, France
^{16} Astronomical Observatory of the Jagiellonian University, Orla 171, 30001 Cracow, Poland
^{17} School of Physics and Astronomy, University of St. Andrews, St. Andrews, KY16 9SS, UK
^{18} INAF–Istituto di Astrofisica Spaziale e Fisica Cosmica Bologna, via Gobetti 101, 40129 Bologna, Italy
^{19} INAF–Istituto di Radioastronomia, via Gobetti 101, 40129 Bologna, Italy
^{20} Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, via della Vasca Navale 84, 00146 Roma, Italy
^{21} INFN, Sezione di Roma Tre, via della Vasca Navale 84, 00146 Roma, Italy
^{22} INAF–Osservatorio Astronomicodi Roma, via Frascati 33, 00040 Monte Porzio Catone (RM), Italy
^{23} Department of Astronomy, University of Geneva, Ch. d’Ecogia 16, 1290 Versoix, Switzerland
^{24} INAF–Osservatorio Astronomico di Trieste, via G. B. Tiepolo 11, 34143 Trieste, Italy
^{25} Department of Astronomy & Physics, Saint Mary’s University, 923 Robie Street, Halifax, Nova Scotia B3H 3C3, Canada
^{26} Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK
^{27} Instituto de Física Teórica, (UAM/CSIC), Universidad Autónoma de Madrid, Cantoblanco, 28049 Madrid, Spain
^{28} Campus of International Excellence UAM+CSIC, Cantoblanco, 28049 Madrid, Spain
^{29} Instituto de Astrofísica de Andalucía (CSIC), Glorieta de la Astronomía, 18080 Granada, Spain
^{30} Departamento de Física Teórica, M8, Universidad Autónoma de Madrid, Cantoblanco, 28049 Madrid, Spain
Received: 17 December 2016
Accepted: 2 August 2017
We carry out a joint analysis of redshiftspace distortions and galaxygalaxy lensing, with the aim of measuring the growth rate of structure; this is a key quantity for understanding the nature of gravity on cosmological scales and latetime cosmic acceleration. We make use of the final VIPERS redshift survey dataset, which maps a portion of the Universe at a redshift of z ≃ 0.8, and the lensing data from the CFHTLenS survey over the same area of the sky. We build a consistent theoretical model that combines nonlinear galaxy biasing and redshiftspace distortion models, and confront it with observations. The two probes are combined in a Bayesian maximum likelihood analysis to determine the growth rate of structure at two redshifts z = 0.6 and z = 0.86. We obtain measurements of fσ_{8}(0.6) = 0.48 ± 0.12 and fσ_{8}(0.86) = 0.48 ± 0.10. The additional galaxygalaxy lensing constraint alleviates galaxy bias and σ_{8} degeneracies, providing direct measurements of f and σ_{8}: [f(0.6),σ_{8}(0.6)] = [0.93 ± 0.22,0.52 ± 0.06] and [f(0.86),σ_{8}(0.86)] = [0.99 ± 0.19,0.48 ± 0.04]. These measurements are statistically consistent with a Universe where the gravitational interactions can be described by General Relativity, although they are not yet accurate enough to rule out some commonly considered alternatives. Finally, as a complementary test we measure the gravitational slip parameter, E_{G}, for the first time at z > 0.6. We find values of E̅_{G}(0.6) = 0.16±0.09 and E̅_{G}(0.86) = 0.09±0.07, when E_{G} is averaged over scales above 3 h^{1} Mpc. We find that our E_{G} measurements exhibit slightly lower values than expected for standard relativistic gravity in a ΛCDM background, although the results are consistent within 1−2σ.
Key words: largescale structure of Universe / cosmology: observations / cosmological parameters / dark energy / galaxies: highredshift
Based on observations collected at the European Southern Observatory, Cerro Paranal, Chile, using the Very Large Telescope under programmes 182.A0886 and partly 070.A9007. Also based on observations obtained with MegaPrime/MegaCam, a joint project of CFHT and CEA/DAPNIA, at the CanadaFranceHawaii Telescope (CFHT), which is operated by the National Research Council (NRC) of Canada, the Institut National des Sciences de l’Univers of the Centre National de la Recherche Scientifique (CNRS) of France, and the University of Hawaii. This work is based in part on data products produced at TERAPIX and the Canadian Astronomy Data Centre as part of the CanadaFranceHawaii Telescope Legacy Survey, a collaborative project of NRC and CNRS. The VIPERS web site is http://www.vipers.inaf.it/
© ESO, 2017
1. Introduction
The origin of the latetime acceleration of the universal expansion is a major question in cosmology. The source of this acceleration and its associated energy density are crucial in understanding the properties of the Universe and its evolution and fate. In the standard cosmological model, this cosmic acceleration can be associated with the presence of a dark energy component, a cosmological fluid with negative pressure, which opposes the gravitational force on large scales. However, this apparent acceleration can conversely be interpreted as a failure of the standard relativistic theory of gravity. A key goal for cosmology is therefore to investigate the nature of gravity empirically. To be clear, what can potentially be falsified is the validity of Einstein’s field equations, rather than General Relativity itself; this sets a broader framework within which Einstein gravity or modified alternatives can operate.
The largescale structure of the Universe has proved to be very powerful for testing the cosmological model through the use of various observables such as the twopoint statistics of the galaxy distribution and its features (e.g. Peacock et al. 2001; Cole et al. 2005; Tegmark et al. 2004; Eisenstein et al. 2005; Guzzo et al. 2008; Percival et al. 2010; Beutler et al. 2011; Blake et al. 2012; Anderson et al. 2014; Alam et al. 2017, and references therein). In this context, a unique probe of gravitational physics is the largescale component of galaxy peculiar velocities affecting the observed galaxy distribution in redshift surveys (Guzzo et al. 2008), sensitive to the growth rate of structure f defined as dlnD/ dlna, where D and a are respectively the linear growth factor and scale factor. In turn, the growth rate of structure tells us about the strength of gravity acting on cosmological scales and is a direct prediction of gravity theories. The distortions induced by peculiar velocities in the apparent galaxy clustering, the socalled redshiftspace distortions (RSD), are a very important cosmological probe of the nature of gravity. In the last decade, they have been studied in large galaxy redshift surveys, showing a broad consistency with Λcold dark matter (ΛCDM) and General Relativity predictions (e.g. Blake et al. 2012; Beutler et al. 2012; de la Torre et al. 2013; Samushia et al. 2014; GilMarín et al. 2016; Chuang et al. 2016).
Although galaxy redshift surveys are powerful cosmological tools for understanding the geometry and the dynamics of the Universe, they are fundamentally limited by the inherent uncertainty related to the bias of galaxies, the fact that these are not faithful tracers of the underlying matter distribution. Gravitational lensing represents a powerful probe that is complementary to galaxy redshiftspace clustering. In the weak regime in particular, the statistical shape deformations of background galaxies probe the relativistic gravitational deflection of light by the projected dark matter fluctuations due to foreground largescale structure. There are several techniques associated with weak gravitational lensing; one that is particularly useful for combining with galaxy clustering is galaxygalaxy lensing. This technique consists of studying the weak deformations of background galaxies around foreground galaxies, whose associated dark matter component acts as a gravitational lens. This is particularly useful for probing the galaxymatter crosscorrelation, which in turn provides insights on the bias of foreground galaxies and the matter energy density Ω_{m}, although the projected nature of the statistic makes it insensitive to RSD. The combination of galaxygalaxy lensing with redshiftspace galaxy correlations is therefore a very promising way to study gravitational physics, given that both lensing information on background sources and spectroscopic information on foreground galaxies are available on the same field.
Beyond the determination of the growth rate of structure, one can define consistency tests of gravity that are sensitive to both the Newtonian and curvature gravitational potentials, Ψ and Φ respectively (e.g. Simpson et al. 2013). One is the gravitational slip, E_{G}, which was originally proposed by Zhang et al. (2007) and implemented by Reyes et al. (2010) in terms of the ratio between the galaxygalaxy lensing signal and the RSD parameter β = f/b times the galaxy clustering signal of the lenses. Here b is the galaxy linear bias. E_{G} effectively tests whether the Laplacian of Ψ + Φ, to which gravitational lensing is sensitive, and that of Ψ, to which galaxy peculiar velocities are sensitive, are consistent with standard gravity predictions. In the standard cosmological model, E_{G} asymptotes to Ω_{m}/f on large linear scales. A failure of this test would either imply an incorrect matter energy density or a departure from standard gravity. This test has been performed at low redshift in the SDSS survey by Reyes et al. (2010) and more recently at redshifts up to z = 0.57 by Blake et al. (2016) and Pullen et al. (2016).
The E_{G} statistic is formally defined as E_{G} = Υ_{gm}/ (βΥ_{gg}), where Υ_{gm} and Υ_{gg} are filtered versions of the realspace projected galaxymatter and galaxygalaxy correlation functions respectively, and β is the RSD parameter. In practice, its implementation involves measuring β and the ratio Υ_{gm}/ Υ_{gg} separately, to finally combine them. But since β and Υ_{gg} are extracted from the same observable, namely the anisotropic twopoint correlation function of lens galaxies, this is suboptimal and does not account for the covariance between them. In this analysis, we follow a different approach. We combine the galaxygalaxy lensing quantity Υ_{gm} and the redshiftspace anisotropic correlation function monopole and quadrupole moments ξ_{0} and ξ_{2} (from which β can be estimated) in a joint likelihood analysis, to provide constraints on f and gravity at redshifts above z = 0.6. We note that we do not include Υ_{gg} because of the redundant cosmological information shared with ξ_{0} and ξ_{2}.
The VIMOS Public Extragalactic Redshift Survey (VIPERS) is a large galaxy redshift survey probing the z ≃ 0.8 Universe with an unprecedented density of spectroscopic galaxies of 5 × 10^{3}h^{3} Mpc^{3} and covering an overall area of 23.5 deg^{2} on the sky. The prime goal of VIPERS is an accurate measurement of the growth rate of structure at redshift around unity. A first measurement has been performed using the Public Data Release 1 (PDR1), setting a reference measurement of fσ_{8} at z = 0.8 (de la Torre et al. 2013). The survey is now complete and several analyses including this one are using the final dataset to produce the VIPERS definitive growth rate of structure measurements, but following a variety of approaches. The present analysis aims at maximizing the cosmological information available and takes advantage of the overlapping lensing information provided by CFHTLenS lensing survey, to provide a precise gravity test at redshifts 0.5 <z < 1.2 by combining RSD and galaxygalaxy lensing.
The paper is organized as follows. The data are described in Sect. 2; Sect. 3 describes our methods for estimating galaxy clustering and galaxygalaxy lensing; Sect. 4 describes the theoretical modelling that is tested in Sect. 5; Sect. 6 presents how the likelihood analysis is constructed; Sect. 7 describes the cosmological results, and Sect. 8 summarizes our findings and concludes.
Throughout this analysis and if not stated otherwise, we assume a flat ΛCDM cosmological model with (Ω_{m},Ω_{b},n_{s}) = (0.3,0.045,0.96) and a Hubble constant of H_{0} = 100 h km s^{1} Mpc^{1}.
2. Data
2.1. Combined VIPERSCFHTLenS dataset
The VIPERS galaxy target sample was selected from the optical photometric catalogues of the CanadaFranceHawaii Telescope Legacy Survey Wide (CFHTLSWide, Goranova et al. 2009). VIPERS covers 23.5 deg^{2} on the sky, divided over two areas within the W1 and W4 CFHTLS fields. Galaxies are selected to a limit of i_{AB} < 22.5, applying a simple and robust gri colour preselection to efficiently remove galaxies at z < 0.5. Coupled with a highly optimized observing strategy (Scodeggio et al. 2009), this allows us to double the galaxy sampling rate in the redshift range of interest, with respect to a pure magnitudelimited sample. At the same time, the area and depth of the survey result in a relatively large volume, 5 × 10^{7}h^{3} Mpc^{3}, analogous to that of the Two Degree Field Galaxy Redshift Survey (2dFGRS) at z ≃ 0.1 (Colless et al. 2001, 2003). Such a combination of sampling rate and depth is unique amongst current redshift surveys at z > 0.5. VIPERS spectra are collected with the VIMOS multiobject spectrograph (Le Fèvre et al. 2003) at moderate resolution (R = 220) using the LR Red grism, providing a wavelength coverage of 5500−9500 Å and a redshift error corresponding to a galaxy peculiar velocity error at any redshift of σ_{vel} = 163 km s^{1}. The full VIPERS area of 23.5 deg^{2} is covered through a mosaic of 288 VIMOS pointings (192 in the W1 area, and 96 in the W4 area). A discussion of the survey data reduction and management infrastructure is presented in Garilli et al. (2014). A complete description of the survey construction, from the definition of the target sample to the actual spectra and redshift measurements, is given in the survey description paper (Guzzo et al. 2014).
The data used here correspond to the publicly released PDR2 catalogue (Scodeggio et al. 2017) that includes 86 775 galaxy spectra, with the exception of a small subset of redshifts (340 galaxies missing in the range 0.6 <z < 1.1), for which the redshift and quality flags were revised close to the release date. Concerning the analysis presented here, this has no effect. A quality flag has been assigned to each object in the process of determining their redshift from the spectrum, which quantifies the reliability of the measured redshifts. In this analysis (as with all statistical analyses presented in the parallel papers of the final science release), we use only galaxies with flags 2 to 9 inclusive, corresponding to objects with a redshift confidence level of 96.1% or larger. This has been estimated from repeated spectroscopic observations in the VIPERS fields (see Scodeggio et al. 2017). The catalogue used here, which we will refer to just as the VIPERS sample in the following, includes 76 584 galaxies with reliable redshift measurements.
In addition to the VIPERS spectroscopic sample, we make use of the public lensing data from the CanadaFranceHawaii Lensing Survey (Heymans et al. 2012), hereafter referred to as CFHTLenS. The CFHTLenS survey analysis combined weak lensing data processing with theli (Erben et al. 2013), shear measurement with lensfit (Miller et al. 2013), and photometric redshift measurement with PSFmatched photometry (Hildebrandt et al. 2012). A full systematic error analysis of the shear measurements in combination with the photometric redshifts is presented in Heymans et al. (2012), with additional error analyses of the photometric redshift measurements presented in Benjamin et al. (2013).
2.2. Sample selection
For this analysis, we define two redshift intervals covering the full volume of the VIPERS survey: 0.5 <z < 0.7 and 0.7 <z < 1.2. The number density of galaxies in the combined W1 and W4 fields is presented in Fig. 1, after correction with survey incompleteness weights w^{C} (see Sect. 3.1). It is worth emphasizing that after application of survey incompleteness corrections, the VIPERS spectroscopic sample represents a statistically unbiased subset of the parent i_{AB} < 22.5 photometric catalogue (Guzzo et al. 2014; Garilli et al. 2014; Scodeggio et al. 2017). The redshift distribution is modelled using the V_{max} method (Cole 2011; de la Torre et al. 2013) and shown with the solid curve in the figure. In this method, we randomly sample 500 times the V_{max} of each galaxy, defined as the comoving volume between the minimum and maximum redshifts where the galaxy is observable given its apparent magnitude and the magnitude limit of VIPERS, i_{AB} = 22.5. The redshift distribution thus obtained is regular and can be straightforwardly interpolated with a smooth function, showed with the solid curve in Fig. 1.
In addition to VIPERS spectroscopic galaxies, photometric galaxies from the CFHTLenS survey on the overlapping areas with VIPERS survey, have been used for the galaxygalaxy lensing. The lens sample satisfies the VIPERS selection i_{AB} < 22.5 and uses VIPERS spectroscopic redshifts when available (i.e. for about 30% of objects) or CFHTLenS maximum likelihood photometric redshifts otherwise. The sources have been selected to have i_{AB} < 24.1 and thus have a higher surface density. Sources inside the mask delimiting bad photometric areas in the CFHTLenS catalogue have been discarded. We also make use of the individual source redshift probability distribution function estimates obtained from bpz (Hildebrandt et al. 2012) as described in Sect. 3.2. Source galaxies extend above z_{phot} = 1.4 and their number density is represented with the unfilled histogram in Fig. 1.
Fig. 1 Number densities of VIPERS galaxies in the individual W1 and W4 fields and of CFHTLenS/VIPERS photometric redshift galaxies, as a function of redshift. The number densities of VIPERS galaxies are corrected for the survey incompleteness by weighting each galaxy in the counts by its associated inverse completeness weight w^{C}. The solid curve corresponds to the model n(z) used in the analysis. It was obtained by randomly sampling galaxy redshifts within their V_{max} (see text for details). 
3. Galaxy clustering and galaxygalaxy lensing estimation
3.1. Anisotropic galaxy clustering estimation
We estimate the redshiftspace galaxy clustering by measuring the twopoint statistics of the spatial distribution of galaxies in configuration space. For this we infer the anisotropic twopoint correlation function ξ(s,μ) using the Landy & Szalay (1993) estimator: (1)where GG(s,μ), GR(s,μ), and RR(s,μ) are respectively the normalized galaxygalaxy, galaxyrandom, and randomrandom number of pairs with separation (s,μ). Since we are interested in quantifying RSD effects, we have decomposed the threedimensional galaxy separation vector s into polar coordinates (s,μ), where s is the norm of the separation vector and μ is the cosine of the angle between the lineofsight and separation vector directions. This estimator minimizes the estimation variance and circumvents discreteness and finite volume effects (Landy & Szalay 1993; Hamilton 1993). A random catalogue needs to be constructed, whose aim is to accurately estimate the number density of objects in the sample. It must be an unclustered population of objects with the same radial and angular selection functions as the data. In this analysis, we use random samples with 20 times more objects than in the data to minimize the shot noise contribution in the estimated correlation functions, and the redshifts of random points are drawn randomly from the model n(z) presented in Fig. 1.
In order to study RSD, we further extract the multipole moments of the anisotropic correlation function ξ(s,μ). This approach has the main advantage of reducing the number of observables, compressing the cosmological information contained in the correlation function. This eases the estimation of the covariance matrices associated with the data. We adopt this methodology in this analysis and use the two first nonnull moments ξ_{0}(s) and ξ_{2}(s), where most of the relevant information is contained, and ignore the contributions of the more noisy subsequent orders. The multipole moments are related to ξ(s,μ) as (2)where L_{ℓ} is the Legendre polynomial of order ℓ. In practice the integration of Eq. (2) is approximated by a Riemann sum over the binned ξ(s,μ). We use a logarithmic binning in s with Δlog (s) = 0.1 and a linear binning in μ with Δμ = 0.02.
VIPERS has a complex angular selection function which has to be taken into account carefully when estimating the correlation function. This has been studied in detail for the VIPERS Public Data Release 1 (PDR1; Guzzo et al. 2014; Garilli et al. 2014) and particularly for the galaxy clustering estimation in de la Torre et al. (2013) and Marulli et al. (2013). We follow the same methodology to account for it in this analysis with only small improvements. We summarize it in the following and refer the reader to the companion paper, Pezzotta et al. (2017), for further details and tests of the method when applied to the VIPERS final dataset.
The main source of incompleteness in the survey is introduced by the VIMOS slit positioner, SSPOC, and the VIPERS onepass observational strategy. This results in an incomplete and uneven spectroscopic sampling, described in detail in Guzzo et al. (2014) and Garilli et al. (2014). In terms of galaxy clustering, the effect is to introduce an underestimation in the amplitude of the measured galaxy correlation function, which becomes scaledependent on the smallest scales. We demonstrate in de la Torre et al. (2013) that this can be corrected by weighting each galaxy in the estimation of the correlation function. For this we define a survey completeness weight, w^{C}, which is defined for each spectroscopic galaxy as well as an angular pair weight, w^{A}, which is applied only to GG pair counts. The latter is obtained from the ratio of one plus the angular correlation functions of targeted and spectroscopic galaxies, as described in de la Torre et al. (2013).
The improvements compared to the PDR1 analysis only concern the estimation of survey completeness weights w^{C}. These in fact correspond to the inverse effective sampling rate, ESR, and are defined for each galaxy as (3)where SSR, TSR are respectively the spectroscopic and target sampling rates (for details, see Guzzo et al. 2014). A significant effort has been invested in improving the estimation of the SSR and TSR. In particular the SSR, which characterizes our ability of measuring the redshifts from observed galaxy spectra, has been refined and now accounts for new galaxy property dependencies, as described in Scodeggio et al. (2017). The TSR, defined as the fraction of spectroscopically observed galaxies in the parent target catalogue, has been recomputed with better angular resolution, on rectangular apertures of 60 by 100 arcsec^{2} around spectroscopic galaxies. In order to mitigate the shot noise contribution in the galaxy counts in such small apertures, we use the Delaunay tesselation that naturally adapts to local density of points (for details, see Pezzotta et al. 2017). The accuracy of this new set of weights is tested in the next section and in Pezzotta et al. (2017).
By applying these weights we effectively upweight galaxies in the pair counts. It is important to note that the spatial distribution of the random objects is kept consistently uniform across the survey volume. The final weights assigned to GG, GR, and RR pairs combine the survey completeness and angular pair weights as where Θ_{ij}(s,μ) is equal to unity for log (s_{ij}) in [log (s)−Δlog (s)/2,log (s) + Δlog (s)/2] and μ_{ij} in [μ−Δμ/ 2,μ + Δμ/ 2], and null otherwise. We define the separation associated with each logarithmic bin as the median pair separation inside the bin. This definition is more accurate than using the bin centre, particularly at large s when the bin size is large.
One can also extract realspace clustering information from the anisotropic redshiftspace correlation function. This can be done by measuring the latter with the estimator of Eq. (1), but where the redshiftspace galaxy separation vector is decomposed in two components, r_{p} and π, respectively perpendicular and parallel to the lineofsight (Fisher et al. 1994). This decomposition allows the isolation of the effect of peculiar velocities as these modify only the component parallel to the lineofsight. This way, RSD can then be mitigated by integrating ξ(r_{p},π) over π, thus defining the projected correlation function (7)We measure w_{p}(r_{p}) using an optimal value of π_{max} = 50 h^{1} Mpc, allowing us to reduce the underestimation of the amplitude of w_{p}(r_{p}) on large scales and at the same time to avoid including noise from uncorrelated pairs with separations of π > 50 h^{1} Mpc. From the projected correlation function, one can derive the following quantity (8)where r_{0} is a cutoff radius, ρ_{c} = 3H^{2}/ (8πG) is the critical density, H(a) = ȧ/a is the Hubble parameter, and G is the gravitational constant. This quantity is equivalent to Υ_{gm}, which is measurable from galaxygalaxy lensing (see next section), but for galaxygalaxy correlations instead of galaxymatter ones. It enters the definition of the gravitational slip parameter E_{G}. In order to measure it in practice, since the logarithmic binning in r_{p} is rather large in our analysis, we interpolate w_{p}(r_{p}) using cubic spline interpolation before evaluating the integral in Eq. (8) numerically. We find that Υ_{gg} is more accurately measured with this technique than by modelling w_{p}(r_{p}) as a power law to perform the integral, as is often done (e.g. Mandelbaum et al. 2013).
3.2. Galaxygalaxy lensing estimation
We use in this analysis the weak lensing technique usually referred to as galaxygalaxy lensing, in which one infers the tangential shear of background sources γ_{t} around foreground objects (lenses) induced by the projected matter distribution in between. This quantity is sensitive to the projected crosscorrelation between lens galaxies and the underlying matter distribution. Since the shear signal is weak and the intrinsic ellipticity of galaxies is unknown, one has to average the former over a large number of foreground sources. The quantity that is effectively measured is the differential excess surface density (9)where (10)In the above equations, r_{p} is the comoving transverse distance between lens and source galaxies, D_{S}, D_{LS}, D_{L} are the angular diameter observersource, lenssource, and observerlens distances, and c is the speed of light in the vacuum.
We use the inverse varianceweighted estimator for the differential excess surface density (e.g. Mandelbaum et al. 2013): (11)where the i and j indices run over source and lens galaxies respectively, N_{S} and N_{L} are respectively the number of source and lens galaxies, e_{t,i} is the tangential ellipticity for each lenssource pair, w^{S} are statistical weights accounting for biases in the determination of background source ellipticities, and Θ_{ij}(r_{p}) is equal to unity for r_{p,ij} in [r_{p}−Δr_{p}/ 2,r + Δr_{p}/ 2] and null otherwise. The projected separation r_{p} is calculated as r_{p} = θχ_{L}, where θ and χ_{L} are respectively the angular distance between the lens and the source, and the radial comoving distance of the lens. This estimator includes an inversevariance weight for each lenssource pair , which downweights the pairs at close redshifts that contribute little to the weak lensing signal (Mandelbaum et al. 2013).
This estimator is unbiased if the redshifts of the sources are perfectly known, but here we have only photometric redshift estimates: the maximum likelihood photometric redshift and the normalized redshift probability distribution function for each source p_{s}(z). Using the maximum likelihood photometric redshift of sources in Eq. (11) and restricting the sum to pairs with z_{S} >z_{L} can possibly lead to a dilution of the signal induced by the nonnegligible probability that z_{S} <z_{L}. This effect can be mitigated by replacing in Eq. (11) by its average over the source redshift probability distribution function p_{s}(12)which leads to the following estimator (e.g. Miyatake et al. 2015; Blake et al. 2016): (13)In principle, those estimators hold in the limit where the lens redshift distribution is narrow and lens redshifts accurate (Nakajima et al. 2012). To better understand the importance of the effects introduced by an imperfect knowledge of the source and lens redshifts in the data, we perform a comparison of different estimates using Eq. (11) or Eq. (13), and various assumptions on the source and lens redshifts. This is presented in terms of the relative difference with respect to a fiducial estimate in Fig. 2. The fiducial estimate is that obtained by using Eq. (13), which includes the individual redshift probability distribution function p_{s}(z) of the sources, and for the lenses, the VIPERS spectroscopic redshift when available or the CFHTLenS maximum likelihood photometric redshift otherwise.
Fig. 2 Relative difference between various estimates of ΔΣ_{gm}, based on different assumptions for source and lens redshifts, and the fiducial estimate in the data at 0.5 <z < 0.7. The quantity shown in the figure is as a function of the projected separation r_{p}. The fiducial estimate is that obtained by using Eq. (13), which includes the individual redshift probability distribution function p_{s}(z) of the sources, and for the lenses, the VIPERS spectroscopic redshift (z_{spec}) when available or the CFHTLenS maximum likelihood photometric redshift (z_{phot}) otherwise (see text). It corresponds to the adopted estimate for the analysis. The grey shaded area represents the relative statistical error expected in the survey. 
We find that the estimate based on Eq. (11), which only uses maximum likelihood photometric redshifts for both lenses and sources, underestimate the signal on all probed scales by about 15% with respect to the fiducial case. Here, we impose z_{S} > 0.1 + z_{L}, including the additive term of 0.1 to account for typical photometric redshift errors (e.g. Coupon et al. 2015). Further including the source redshift probability distribution function through the estimator of Eq. (13) allows a slight improvement, reaching an underestimation of about 10% with respect to the fiducial case. The two previous estimates are still affected by the uncertainty on the lens redshifts, which effectively tends to dilute the overall signal. If we now use as lenses only VIPERS spectroscopic galaxies, which represents about 30% of all galaxies with i_{AB} < 22.5, we find a remarkably good agreement with the fiducial estimate. In principle, this estimate may be considered as the reference unbiased estimate, however on the largest scales probed by the data, i.e. at r_{p} = 10−20 h^{1} Mpc, the signal drops significantly. This can be imputed to the lack of sourcelens pairs induced by the reduced number of lenses, directly affecting our ability to probe the largest scales signal. However, we find that this effect can be mitigated by adding photometric lenses from the CFHTLenS catalogue, taking the maximum likelihood photometric redshifts: this corresponds to the fiducial estimate. We note that the expected statistical uncertainty, which is shown in Fig. 2 with the grey shaded area, is not negligible particularly above r_{p} = 10 h^{1} Mpc, and higher than any residual systematic effect. This test makes us confident that our fiducial estimate of ΔΣ_{gm}(r_{p}) is robust, given the expected level of statistical error in the data. Similar results are found at 0.7 <z < 1.2, leading to the same conclusions.
A nonnegligible source of systematics in weak lensing measurements is related to the measurement of background galaxy shapes. This can lead to systematic biases in the lensing measurements. The CFHTLenS collaboration has studied these extensively in Miller et al. (2013) and Heymans et al. (2012), and we follow their method to correct our measurements. We used the additive and multiplicative shear calibration corrections c and m, as well as the optimal weights w^{S} provided by lensfit, which are available in the CFHTLenS catalogue. In particular, to correct for the multiplicative bias we applied the correction factor to ΔΣ_{gm}(r_{p}) as described in Miller et al. (2013) and Velander et al. (2014). We found this correction to boost the galaxygalaxy lensing signal by about 5% independently of the scale.
For the purpose of constraining the cosmological model, it can be difficult to use ΔΣ_{gm} as its modelling is nonlinear. One of the difficulties is to model the nonlinear scales and the intrinsic mixing of smallscale nonlinear and largescale linear information (Baldauf et al. 2010). This is achievable but at the expense of introducing additional nuisance parameters in the model (e.g. Cacciato et al. 2013; More et al. 2015). An alternative approach, which we use in this analysis, consists of using a derived statistic that allows the mitigation of nonlinearities: the annular differential surface density Υ_{gm}, which is defined as (Baldauf et al. 2010) (14)This statistic removes the smallscale nonlinear contribution of ΔΣ_{gm} below a cutoff radius r_{0}. We use this quantity in our analysis and study the impact of the choice of r_{0} in Sect. 5.1.
4. Theoretical modelling
4.1. Galaxy biasing
Galaxies are not faithful tracers of the underlying matter distribution and this has to be taken into account in cosmological analyses, since cosmological models primarily predict matter observables. The modelling of galaxy biasing is simplified when focusing on large scales, where bias can be considered as linear and simply be represented as a constant multiplicative factor in front of the matter power spectrum. This is a common assumption in RSD analyses. In our case, however, the relatively small survey volume means that much of our information lies below fully linear scales; for this reason, and because of the intrinsic nonlinearities in the excess surface density ΔΣ_{gm}, additional care must be taken to model galaxy biasing. We use a nonlinear prescription for galaxy bias based on the cosmological perturbation theory that allows describing it more accurately down to translinear scales. We adopt the nonlinear nonlocal bias model of McDonald & Roy (2009) that relates the galaxy overdensity δ_{g} and matter overdensity δ as: (15)where b_{1} and b_{2} are the linear and secondorder nonlinear bias terms, b_{s2} the nonlocal bias term, s is the tidal tensor term from which nonlocality originates. The σ^{2} and ⟨ s^{2} ⟩ terms ensure the condition ⟨ δ_{g} ⟩ = 0.
4.2. Annular differential excess surface density
The galaxygalaxy lensing quantity that we observe is the differential excess surface density. It is defined as (16)where (17)and Σ_{gm}(r_{p}) is the projected surface density defined as (18)In the above equation, Ω_{m} is matter energy density and χ is the radial comoving coordinate. Υ_{gm} can be predicted from ΔΣ_{gm} by using Eq. (14) or directly from the galaxymatter crosscorrelation function as (Baldauf et al. 2010) (19)where W_{Υ}(x,r_{p},r_{0}) is the window function (Baldauf et al. 2010): (20)where Θ(x) is the Heaviside step function.
From these equations one can see explicitly that Υ_{gm} is related to the galaxymatter crosscorrelation function ξ_{gm} or crosspower spectrum P_{gm}. If we assume the biasing model of Eq. (15), P_{gm} can be written as (McDonald & Roy 2009) (21)where P_{δδ} is the nonlinear matter densitydensity power spectrum, b_{3nl} is a thirdorder nonlocal bias term, P_{lin} is the linear matter power spectrum, and P_{b2,δ}, P_{bs2,δ} are 1loop integrals given in Appendix A. In the local Lagrangian picture where one assumes no initial nonlocal bias, one can predict that the nonlocal bias terms at later time are related to b_{1} such that (Chan et al. 2012; Saito et al. 2014) We adopt these relations and our model has finally two galaxy biasing parameters: b_{1} and b_{2}, b_{1} being the standard linear bias parameter.
4.3. Redshiftspace distortions
The most general formalism describing the redshiftspace anisotropies in the power spectrum derives from writing the matter density conservation in real and redshift space (Kaiser 1987). In particular, in the planeparallel approximation that is assumed in this analysis, the anisotropic power spectrum of matter has the general compact form (Scoccimarro et al. 1999) (24)where ν = k_{∥}/k, u_{∥}(r) = −v_{∥}(r)/(faH(a)), v_{∥}(r) is the lineofsight component of the peculiar velocity, δ is the matter density field, Δu_{∥} = u_{∥}(x)−u_{∥}(x′) and r = x−x′. It is worth noting that in Fourier space, for an irrotational velocity field, is related to the divergence of the velocity field θ via . Although exact, Eq. (24) is impractical and we use the approximation proposed by Taruya et al. (2010). In the case of perfect matter tracers, the latter model takes the form (25)where D(kνσ_{v}) is a damping function, P_{δδ}, P_{δθ}, P_{θθ} are respectively the nonlinear matter densitydensity, densityvelocity divergence, and velocity divergencevelocity divergence power spectra, and σ_{v} is an effective pairwise velocity dispersion that we can fit for and then treat as a nuisance parameter. The expressions for C_{A}(k,ν,f) and C_{B}(k,ν,f) are given in Taruya et al. (2010) and de la Torre & Guzzo (2012). This phenomenological model can be seen in configuration space as a convolution of a pairwise velocity distribution, the damping function D(kμσ_{v}) that we assume to be Lorentzian in Fourier space, i.e. (26)and a term involving the density and velocity divergence correlation functions and their spherical Bessel transforms.
This model can be generalized to the case of biased tracers by including a biasing model. By introducing that of Eq. (15), one obtains for the redshiftspace galaxy power spectrum (Beutler et al. 2014; GilMarín et al. 2014) (27)where In the above equations P_{δθ} is the nonlinear matter densityvelocity divergence power spectrum, P_{lin} is the matter linear power spectrum, and P_{b2,δ}, P_{bs2,δ}, P_{b2,θ}, P_{bs2,θ}, P_{b22}, P_{b2s2}, P_{bs22}, are 1loop integrals given in Appendix A.
The final model for is obtained from its Fourier counterpart defined as (30)where (31)In the above equation, j_{ℓ} denotes the spherical Bessel functions.
The ingredients of the model are the nonlinear power spectra of density and velocity divergence at the effective redshift of the sample. These power spectra can be predicted from perturbation theory or simulations for different cosmological models. The nonlinear matter power spectrum can also be obtained to a great accuracy from semianalytical prescriptions such as HALOFIT (Smith et al. 2003), for various cosmologies. In particular, HALOFIT allows the prediction of P_{δδ} from the linear matter power spectrum and the knowledge of the scale of nonlinearity at the redshift of interest, k_{nl}(z). We note that at fixed linear matter power spectrum shape, variations of σ_{8}(z) can be straightforwardly mapped into variations of k_{nl}(z) (see Smith et al. 2003). In this analysis, the linear matter power spectrum is predicted using the CLASS Boltzmann code (Lesgourgues 2011), and we use the latest calibration of HALOFIT by Takahashi et al. (2012) to obtain P_{δδ}. To predict P_{θθ} and P_{δθ}, we use the nearly universal fitting functions of Bel et al. (in prep.) that depend on the linear power spectrum and σ_{8}(z) as where P_{lin} is the linear power spectrum and (m_{1},m_{2},n_{1},n_{2}) are free parameters calibrated on simulations. We adopt here the values (m_{1},m_{2},n_{1},n_{2}) = (1.906,2.163,2.972,2.034). These predictions for P_{θθ} and P_{δθ} are accurate at the few percent level up to k ≃ 0.7 (Bel et al., in prep.). Therefore, the overall degree of nonlinearity in P_{δδ}, P_{δθ} and P_{θθ} is solely controlled by σ_{8}(z), which is left free when fitting the model to observations.
In the model, the linear bias and growth rate parameters, b_{1} and f, are degenerate with the normalization of the matter power spectrum parameter σ_{8}. Generally with RSD, only the combination of b_{1}σ_{8} and fσ_{8} can be constrained if no assumption is made on the actual value of σ_{8}. However in the Taruya et al. (2010) model, , , and terms appear in the correction term C_{A} (see Taruya et al. 2010; de la Torre & Guzzo 2012). Accordingly, in the general case, (f,b_{1},b_{2},σ_{8},σ_{v}) are treated as separate parameters in the fit and we provide marginalized constraints on the derived fσ_{8}.
4.4. Redshift errors
Redshift errors can potentially affect the anisotropic RSD signal. In the anisotropic correlation function they have a similar effect as galaxy random motions in virialized objects: they introduce a smearing of the correlation function along the line of sight at small transverse separations. If the probability distribution function of redshift errors is known, their effect can be forward modelled by adding another multiplicative damping function in the redshiftspace galaxy power spectrum of Eq. (19). In that case, the damping function should be the Fourier transform of the error probability distribution function. We follow this approach and the final model is obtained by multiplying Eq. (19) by a Gaussian with standard deviation set to the estimated pairwise redshift dispersion of VIPERS galaxies such that the final RSD model reads (34)where is taken from Eq. (27), G is the Fourier transform of the Gaussian kernel (35)and σ_{z} is the pairwise standard deviation associated with the redshift error probability distribution function.
Fig. 3 Probability distribution function of redshift errors at 0.5 <z < 0.7 and 0.7 <z < 1.2 in the VIPERS data. This is obtained from the redshift differences of reobserved galaxies, for which there are two independent redshift measurements. The dotted and dashed curves are bestfitting Gaussians for the redshift intervals 0.5 <z < 0.7 and 0.7 <z < 1.2 respectively. 
The Gaussian form is motivated by the data themselves as shown in Fig. 3. In this figure are shown the distributions of redshift differences at 0.5 <z < 0.7 and 0.7 <z < 1.2 in VIPERS reobservations (1061 at 0.5 <z < 0.7 and 1086 at 0.7 <z < 1.2), for which we have two independent redshift measurements for the same galaxies (see Scodeggio et al. 2017). These distributions can be rather well modelled by Gaussians, and by doing so, we obtain values of σ_{z} = 1.31 × 10^{3} and σ_{z} = 1.36 × 10^{3} for the pairwise redshift standard deviations at 0.5 <z < 0.7 and 0.7 <z < 1.2 respectively. These are further converted in comoving length assuming the fiducial cosmology to enter the model in Eq. (34).
4.5. AlcockPaczynski effect
Additional distortions can arise in galaxy clustering because of the need to assume a fiducial cosmology to convert redshift and angular positions into comoving distances, and the fact that this fiducial cosmology is not necessarily the true one. This is the Alcock & Paczynski (1979) effect (AP). More specifically, since the lineofsight separations require the knowledge of the Hubble parameter, H(z), and transverse separations that of the angular diameter distance, D_{A}(z), any difference in H(z) and D_{A}(z) between the fiducial and true cosmologies, translates into an anisotropic clustering, independently of RSD. Although AP and RSD anisotropies are degenerate to some extent in the observables (Ballinger et al. 1996; Matsubara & Suto 1996), they have a fundamentally different origin: AP is sensitive to the geometry whereas RSD are sensitive to the growth of cosmological perturbations.
We follow Xu et al. (2013) and model AP distortions using the α and ϵ parameters, which characterize respectively the isotropic and anisotropic distortion components associated with AP. These are given by where quantities calculated in the fiducial cosmology are denoted with primes. Those parameters modify the scales at which the correlation function is measured such that Therefore, for the model correlation function monopole and quadrupole in a tested cosmology, the corresponding quantities in the fiducial cosmology are obtained as (Xu et al. 2013) In the case of the galaxygalaxy lensing statistic that we are considering, since it is a function of the transverse separation r_{p}, the corresponding Υ_{gm} in the fiducial cosmology is simply given by (42)
4.6. Cosmological insights from galaxy clustering and galaxygalaxy lensing
Gravitational physics on cosmological scales can be tested from measurements of the growth rate of structure, which is well measured from RSD in the galaxy clustering pattern. We have seen that in practice, the correlation function multipole moments depend not only on the growth rate of structure f, but also on the shape and amplitude σ_{8} of the matter power spectrum, the galaxy bias parameters b_{1} and b_{2}, and the pairwise velocity dispersion σ_{v}. To derive the growth rate of structure, one then needs to marginalise over those nuisances. This is of course a source of uncertainty in the determination of the growth rate of structure. Moreover, since there is a degeneracy between the amplitude of the matter power spectrum σ_{8}, the growth rate of structure f, and the linear bias parameter b_{1}, RSD alone are sensitive to the fσ_{8} and b_{1}σ_{8} parameter combinations.
On the other hand, galaxygalaxy lensing probes the realspace galaxymatter correlations that are described by the shape and amplitude σ_{8} of the matter power spectrum, the galaxy bias parameters b_{1} and b_{2}, and the matter density parameter Ω_{m}. Projected galaxygalaxy correlations are also sensitive to σ_{8}, b_{1}, and b_{2}. But by looking in detail at those dependencies, we can see that in the linear regime , while , such that by combining the two we can break the degeneracy between b_{1} and σ_{8}. We note that ξ(s,μ), from which ξ_{0} and ξ_{2} are derived, has the same parameter dependences as Υ_{gg}, except for the additional f dependence. Therefore, additional galaxygalaxy lensing information brings an independent handle on the bias parameters b_{1} and b_{2}, and the power spectrum amplitude σ_{8}, reducing the uncertainties on the growth rate of structure induced by the lack of knowledge on the bias of galaxies, as well as a supplementary sensitivity to Ω_{m}.
5. Tests on simulated data
5.1. Simulated data
To test the robustness of redshiftspace galaxy clustering, galaxygalaxy lensing, and associated error estimates, we make use of a large number of mock galaxy samples, which are designed to be a realistic match to the VIPERS final dataset. We used the mock lensing lightcones presented in Giocoli et al. (2016). These have been built upon the Big MultiDark dark matter Nbody simulation (Klypin et al. 2016), which assumes a flat ΛCDM cosmology with (Ω_{m}, Ω_{Λ}, Ω_{b},h,n,σ_{8}) = (0.307,0.693,0.0482,0.678,0.960,0.823) and covers a volume of 15.625 h^{3} Gpc^{3}. These lightcones contain the shear information associated with simulated background galaxies distributed uniformly on the sky but following the redshift distribution of CFHTLenS galaxies. More specifically, the lightcones have been built to match the effective number density and redshift distribution of the CFHTLenS lensing catalogue. We added Gaussian random errors with standard deviation to the ellipticities to mimic those in the CFHTLenS data. The size of the simulation allowed us to create 54 independent lightcones for W1 and W4, spanning the redshift range 0 <z < 2.3 (for details, see Giocoli et al. 2016).
We populate these lightcones with foreground galaxies using the halo occupation distribution (HOD) technique and apply the detailed VIPERS selection function and observational strategy. The haloes were identified in the simulation using a friendsoffriends algorithm with a relative linking length of b = 0.17 times the interparticle separation. The mass limit to which the halo catalogues are complete is 10^{11.95}h^{1}M_{⊙}. Because this limiting mass is too large to host the faintest galaxies observed with VIPERS, we use the method of de la Torre & Peacock (2013) to reconstruct haloes below the resolution limit. This method is based on stochastically resampling the halo number density field using constraints from the conditional halo mass function. For this, one needs to assume the shapes of the halo bias factor and halo mass function at masses below the resolution limit and use the analytical formulae obtained by Tinker et al. (2008, 2010). With this method we are able to populate the simulation with lowmass haloes with a sufficient accuracy to have unbiased galaxy twopoint statistics in the simulated catalogues (for details, see de la Torre et al. 2013). The minimum reconstructed halo mass we consider for the purpose of creating VIPERS mocks is 10^{10}h^{1}M_{⊙}.
In this process, we populate each halo with galaxies according to its mass, the mean number of galaxies in a halo of a given mass being given by the HOD. It is common usage to differentiate between central and satellite galaxies in haloes. While the former are put at rest at halo centres, the latter are randomly distributed within each halo according to a NFW radial profile (Navarro et al. 1996, 1997). The halo occupation function and its dependence on redshift and luminosity/stellar mass must be precisely chosen in order to obtain mock catalogues with realistic galaxy clustering properties. We calibrated the halo occupation function directly on the VIPERS data, as presented in de la Torre et al. (2013). We add velocities to the galaxies and measure their redshiftspace positions. While the central galaxies are assigned the velocity of their host halo, satellite galaxies have an additional random component for which each Cartesian velocity component is drawn from a Gaussian distribution with a standard deviation that depends on the mass of the host halo. Details about the galaxy mock catalogue construction technique are given in Appendix A of de la Torre et al. (2013).
The final step in obtaining fully realistic VIPERS mocks is to add the detailed survey selection function. We start by applying the magnitude cut i_{AB} < 22.5 and the effect of the colour selection on the radial distribution of the mocks. This is achieved by depleting the mocks at z < 0.6 so as to reproduce the VIPERS colour sampling rate (see Guzzo et al. 2014, for detail). The mock catalogues that we obtain are then similar to the parent photometric sample in the data. We next apply the slitpositioning algorithm with the same setting as for the data. This allows us to reproduce the VIPERS footprint on the sky, the smallscale angular incompleteness and the variation of TSR across the fields. Finally, a random redshift error is added to the redshifts as in the data. We are thus able to produce realistic mock galaxy catalogues that contain the detailed survey completeness function and observational biases of VIPERS, which we refer to as the ‘observed’ mock catalogues in the following.
We note that another set of VIPERS mock catalogues spanning the redshift range of 0.4 <z < 1.2 have been constructed. This set, which comprises 306 and 549 lightcones of W1 and W4 fields respectively, has not been explicitly used in this analysis, but in accompanying VIPERS PDR2 analyses (e.g. Hawken et al. 2017; Pezzotta et al. 2017; Wilson et al., in prep.; Rota et al. 2017).
5.2. Systematics on the correlation function monopole and quadrupole
Fig. 4 Relative systematic errors on the correlation function monopole (top panel) and quadrupole (bottom panel) at 0.5 <z < 0.7 and effects of target sampling rate (TSR) and angular pair weighting (w^{A}) corrections. The grey shaded areas represent the relative statistical error expected in the survey, while light grey band mark ± 1% relative uncertainties for reference. 
The mock samples are crucial for testing the redshiftspace clustering estimation in VIPERS, which is not trivial given the complex selection function of the survey. We first study the impact of the survey selection function on the measurement of the monopole and quadrupole correlation functions. We measured these quantities in the observed mocks, applying the different weights defined in Sect. 3.1, and compare them to the reference measurements obtained from the parent mocks, including VIPERS typical spectroscopic redshift errors. The relative differences in ξ_{0} and ξ_{2} as a function of separation and averaged over the mocks are shown in Figs. 4 and 5, respectively for the two samples at 0.5 <z < 0.7 and 0.7 <z < 1.2. First of all, it is clear from these figures that without any correction the spectroscopic strategy introduces biases in the estimation of the galaxy clustering. But when applying the survey completeness weights w^{C}, one can recover within a few percent the correct amplitude of the correlation functions on scales above 5 h^{1} Mpc. By further applying the angular weights w^{A}, we obtain an almost unbiased estimate of the monopole and quadrupole down to a few h^{1} Mpc. The statistical relative error induced by sample variance and estimated from the dispersion among the mock samples, is shown with the shaded area in these figures. It is important to note that it is much larger than any residual systematics over the range of scales considered. Finally, it is worth mentioning that in the quadrupole, the apparent higher level of systematics at around s = 10 h^{1} Mpc is an artefact due to the zero crossing of the functions at slightly different separations.
5.3. Systematics on the growth rate of structure
Fig. 6 Relative systematic error on fσ_{8} at 0.5 <z < 0.7 as a function of s_{min}, for different values of r_{0} (r_{0} = 1 h^{1} Mpc and r_{0} = 1.5 h^{1} Mpc) and when including or not redshift error. The error bars represent the relative statistical error associated to analysing the mean mock predictions. The shaded area shows the 1σ confidence region associated with the relative statistical error expected in VIPERS. The squares and triangles are artificially shifted along s_{min} axis to improve the clarity of the figure. 
We further study our ability to determine fσ_{8} when combining RSD and galaxygalaxy lensing measurements in a maximum likelihood analysis. For this purpose we perform several analyses of the mean RSD and galaxygalaxy lensing measurements in the observed mocks, for different minimum separations s_{min} in the correlation functions and different cutoff scale r_{0} in the annular differential excess surface density. These analyses are performed on the mean quantities to reduce the impact of statistical errors and concentrate on systematics. The precision matrix is estimated from the mocks as explained in Sect. 6, except that each element is further divided by the number of mocks to characterize the error on the mean. As an illustration, we present in this section only the case of the sample at 0.5 <z < 0.7. The sample at 0.7 <z < 1.2 provides very similar systematic levels.
Figure 6 presents the systematic errors on fσ_{8}, i.e. the relative difference of recovered values with respect to the fiducial value of the mocks, as a function of s_{min} and for r_{0} = (1 h^{1} Mpc,1.5 h^{1} Mpc). We consider rather small minimum scales and cutoff radii to explore the extent to which our modelling is robust in the translinear regime. We can see in this figure that our model allows the recovery of the fiducial value of fσ_{8} down to s_{min} = 6.3 h^{1} Mpc, with systematic errors below 5%, independently of the choice of r_{0}. In principle, values of r_{0} smaller than the typical radius of haloes hosting these galaxies may not be optimal, since the nonlinear contribution to correlations may dominate on those scales, which are more difficult to describe (Baldauf et al. 2010). However, this also depends on the galaxy type and the redshift. For VIPERS galaxies and the considered biasing model, we find that r_{0} = 1 h^{1} Mpc can be well described by our model (see also Blake et al. 2016). This can be seen in Fig. 7 where is shown the comparison between the mean mock ΔΣ_{gm} and Υ_{gm} obtained with r_{0} = (1 h^{1} Mpc,1.5 h^{1} Mpc) and the predictions of our model, when b_{2} is allowed to vary and b_{1} is fixed to its fiducial value. We can see that although the model fails to reproduce ΔΣ_{gm} on scales below about 3 h^{1} Mpc, it provides a good description of Υ_{gm} for b_{2} = −0.1.
Fig. 7 Comparison of differential excess surface density Δ_{gm} and annular differential excess surface density Υ_{gm} predictions in the mocks (points and shaded regions) and by our theoretical model (curves) at 0.5 <z < 0.7. The mock predictions correspond to the mean signal among the mock realizations (points) and its associated 1σ error (shaded region). The curves show the theoretical predictions for the fiducial parameters of the mocks, varying only the b_{2} parameter as labeled. 
We finally test the impact of redshift errors on the recovery of fσ_{8} in Fig. 6. This figure shows the relative systematic error on fσ_{8} as a function s_{min} in the case where r_{0} is fixed to 1 h^{1} Mpc and typical VIPERS redshift errors are added randomly to mock galaxy redshifts. By comparing it to the case without redshift errors, we can see that for the RSD model where a Gaussian damping term is added to account for redshift errors, the recovery of fσ_{8} is achieved without additional bias, with only a small relative bias of about 3% at s_{min} = 6.3 h^{1} Mpc and −5% above. We note that this is the ideal case where the redshift error probability distribution function is perfectly known.
Overall, these tests demonstrate that our model with s_{min} = 6.3 h^{1} Mpc and r_{0} = 1 h^{1} Mpc is robust enough to provide a precise measurement of fσ_{8} with VIPERS data, with residual systematics of the order of a few per cent only, but only representing about one fifth of expected statistical error on the measurement as shown with the grey shaded region in Fig. 6. Based on these tests, we adopt (s_{min},r_{0}) = (6.3 h^{1} Mpc,1 h^{1} Mpc) values for the following analysis.
6. Likelihood analysis and precision matrix
In order to derive cosmological parameters from the combination of RSD and galaxygalaxy lensing measurements, we perform a maximum likelihood analysis in which we define the likelihood function ℒ such that (43)where N_{p} is the number of data points in the fit, Δ is the datamodel difference vector, and is the inverse data covariance matrix. Δ is defined such that each element is Δ_{i} = d_{i}−m_{i}, where d and m are respectively the data and model prediction vectors. In our case, d is the concatenation of ξ_{0}, ξ_{2}, and Υ_{gm}, for the set of considered separations. The parameter space of the model is explored using a Monte Carlo Markov chain (MCMC) method implementing the MetropolisHastings algorithm (Metropolis et al. 1953).
A robust estimation of the inverse data covariance matrix, or precision matrix, is crucial in order to achieve realistic posterior likelihood functions of the parameters. The different bins in ξ_{0}, ξ_{2}, and Υ_{gm} are correlated to some degree and this must be allowed for in the likelihood analysis. We measure these three quantities in the 54 mocks and estimate the covariance matrix C. The generic elements of the matrix can be evaluated as (44)where N_{m} is the number of mock realizations and the indices i,j run over the data vector d elements. An unbiased estimate of the inverse covariance matrix, , is obtained as (Hartlap et al. 2007) (45)for N_{m} >N_{p}−2. The resulting inverse covariance matrix obtained from mock realizations can be noisy, depending on how large N_{m} is with respect to N_{p}. In our case, N_{m} = 54 and N_{p} = 16, which suggests the presence of a nonnegligible noise in the inverse covariance matrix. In order to reduce the level of noise, we adopt the tapering technique of Kaufman et al. (2008). This technique has been introduced in the context of cosmological analysis by Paz & Sánchez (2015). This technique relies on the assumption that largescale covariances vanish, and consists of tapering the covariance matrix around the diagonal using a specific positive and compact taper function. Contrary to other estimators such as shrinkage (e.g. Pope & Szapudi 2008), the twotapers estimator has the advantage of being unbiased. The inverse tapered covariance matrix is obtained as (46)where “°” denotes the elementwise matrix product and T is the tapering matrix. We follow Paz & Sánchez (2015) and use the tapering matrix defined as (47)where x_{i} is the ith measurement position in the data vector, and K is the taper function that we take to be a Wendland function: This taper function has one free parameter, the tapering scale T_{p}, which essentially represents the typical scale difference above which covariances are nullified.
Fig. 8 Top panel: recovered errors on fσ_{8} normalized to that obtained without tapering (equivalent to applying tapering with T_{p} = ∞), as a function of the tapering scale T_{p} used in the estimation of the precision matrix. This is obtained from the mocks at 0.5 <z < 0.7. Bottom panel: recovered maximum likelihood values for fσ_{8} and associated 1σ error as a function of the tapering scale T_{p}. 
Fig. 9 Correlation matrix (left panel) and normalized precision matrix (right panel) for galaxy clustering and galaxygalaxy lensing data in the redshift interval 0.5 <z < 0.7. These are defined as C_{ij}/ and / respectively, where C_{ij} and refer to covariance and precision matrix elements respectively. In both panels, the upper triangular matrix represents the case without tapering, while the lower panel the case with tapering. The precision matrix is normalized such that diagonal elements are unity. 
In our case, the covariance matrix is associated with three different quantities as well as two different separation types, s and r_{p}. One would then potentially need to use a combination of several taper functions, since one does not expect the largescale covariance to vanish at the same scales for all quantities. Although it may be suboptimal to use a single taper function, we still expect to increase the signaltonoise, and since the estimator is unbiased, one cannot introduce additional bias or error. We therefore decided to use a single taper function for simplicity.
In the general case, it is not straightforward to define a priori the optimal tapering scale. Paz & Sánchez (2015) introduced a simple empirical method, which consists of performing several maximum likelihood analyses of the data varying only the tapering scale, and taking as the optimal T_{p} the one that minimizes the error on the parameter of interest. We perform the same exercise on the mean mock predictions. The marginalized 1σ error on fσ_{8} as a function of T_{p} is presented in the top panel of Fig. 8. We can see that the T_{p} value that minimizes the error is around 15 h^{1} Mpc, and we adopt this value in our analysis. We also verified that the maximum likelihood values for fσ_{8} remain unchanged for any value of T_{p} as shown in the bottom panel of Fig. 8.
To illustrate the method, we present in Fig. 9 the correlation matrix and normalized precision matrix, for the combined RSD and galaxygalaxy lensing data in the redshift interval 0.5 <z < 0.7, when applying or not the tapering technique (lower and upper triangles respectively). Those matrices are defined as C_{ij}/ and / respectively, where C_{ij} and refer to covariance and precision matrix elements respectively. We can see the reduction of noise, which is particularly clear in the normalized precision matrix for most offdiagonal terms.
The tapering technique allows a significant reduction of the noise level in the precision matrix, but cannot completely remove it. The remaining noise can propagate through the likelihood analysis into derived parameter uncertainties. In order to obtain realistic confidence limits on parameters one needs to account for the additional uncertainties coming from the precision matrix estimation (Taylor & Joachimi 2014). Percival et al. (2014) showed that this additional error can be described as a rescaling of the target parameter covariance, in the case when the precision matrix is estimated with the standard estimator of Eq. (45). But the appropriate degree of rescaling is unclear when the tapering estimator is used. The improvement on the error that we find with the tapering estimator (i.e. 26.5%) is similar to or larger than what we would expect with the standard estimator using 300 mocks or more as predicted by (Dodelson & Schneider 2013; Percival et al. 2014). This gives us confidence that only a small correction, if any, would be necessary.
7. Cosmological results
The comprehensive tests of the methodology described in previous sections make us confident that we can perform a robust combined analysis of RSD and galaxygalaxy lensing with VIPERS and CFHTLenS dataset, and infer cosmology from it. We present in this section the data measurements, growth rate of structure constraints, and derived gravitational slip parameters at 0.5 <z < 0.7 and 0.7 <z < 1.2.
7.1. Galaxy clustering and galaxygalaxy lensing measurements
The correlation function measurements are performed on the full VIPERS galaxy sample in the redshift intervals 0.5 <z < 0.7 and 0.7 <z < 1.2. We select all VIPERS galaxies above the limiting magnitude of the survey, and measure the monopole and quadrupole correlation functions in both W1 and W4 fields. The combined W1+W4 measurements are obtained by summing up the pairs in the two fields, contributing to the anisotropic twopoint correlation functions ξ(s,μ), before deriving ξ_{0} and ξ_{2} from Eq. (2). The full anisotropic twopoint correlation functions are presented in Fig. 10, and the monopole and quadrupole moments in Fig. 11. In the latter figure, the individual mock measurements are superimposed, giving a visual appreciation of the error associated with these measurements in VIPERS. We can see that the combined W1+W4 monopole and quadrupole correlation function measurements enable us to probe accurately the redshiftspace galaxy clustering signal on scales below about s = 50 h^{1} Mpc.
Fig. 10 Anisotropic correlation functions of VIPERS galaxies at 0.5 <z < 0.7 (top panel) and 0.7 <z < 1.2 (bottom panel) as a function parallel and transverse to the lineofsight separations. 
Fig. 11 Monopole (circles) and quadrupole (triangles) correlation functions of VIPERS galaxies at 0.5 <z < 0.7 (top panel) and 0.7 <z < 1.2 (bottom panel). Solid curves correspond to individual mock measurements. 
The differential excess surface density measurements are obtained by combining W1 and W4 individual field measurements in a similar fashion. The lens galaxies are taken from the VIPERS catalogue or the CFHTLenS catalogue if no spectroscopic redshift is available. They are selected to have i_{AB} < 22.5 and a redshift in the intervals 0.5 <z < 0.7 and 0.7 <z < 1.2. The source galaxies are taken from the CFHTLenS catalogue and are selected to have i_{AB} < 24.1. The differential excess surface density and annular differential excess surface density measurements for r_{0} = 1 h^{1} Mpc are presented in Fig. 12. As in Fig. 11, the individual mock measurements are superimposed. We can see that with the combined W1+W4 annular differential excess surface density measurements we can reach scales up to about r_{p} = 20 h^{1} Mpc.
Fig. 12 Differential excess surface density (circles) and annular differential excess surface density (triangles) at 0.5 <z < 0.7 (top panel) and 0.7 <z < 1.2 (bottom panel). Solid curves correspond to individual mock measurements. 
Our Υ_{gm} measurements are more uncertain than the ξ_{0} and ξ_{2} ones. This is essentially related to the way the former are estimated. Weak lensing is fundamentally limited by the unknown intrinsic ellipticity of the sources, which dominates the error budget. This can be mitigated by means of a larger number of sources. Given the surface density of sources in our sample and its rather modest angular coverage of 23.5 deg^{2}, we obtain relative errors on Υ_{gm} of about 25%, estimated from the mock samples. In contrast, the typical relative error that we obtain on ξ_{0} is of 5%. Therefore, in our combined analysis of the RSD and galaxygalaxy lensing we expect Υ_{gm} to have a much lower weight in the likelihood. We finally remark that the observed Υ_{gm} tend to exhibit lower amplitudes than expected in the mock samples, in particular in the highest redshift interval. We discuss the cosmological implications of this in Sect. 7.3.
7.2. Growth of structure constraints
We perform a combined maximum likelihood analysis of the monopole, quadrupole, and annular differential excess surface density to derive constraints on the growth rate of structure at 0.5 <z < 0.7 and 0.7 <z < 1.2. The effective redshifts associated with these intervals are z = 0.6 and z = 0.86. They correspond to the average redshift of pairs contributing the most to monopole and quadrupole correlation functions in these redshift intervals (Samushia et al. 2014). The theoretical model that we use is described in Sect. 4; it depends on 11 parameters, p = (f,b_{1},b_{2},σ_{v},σ_{8},ϵ,α,Ω_{m},Ω_{m}h^{2},Ω_{b}h^{2},n_{s}). The last three describe the shape of the matter power spectrum and these are determined most accurately by CMB data. Since our galaxy clustering and weak lensing measurements cannot provide such tight constraints on these parameters, we fix them to the bestfitting Planck 2015 TT+lowP+lensing parameters (Planck Collaboration XIII 2016). Consistently, Ω_{m} is kept fixed to the Planck value in Υ_{gm}. Possible departures from those parameter values are only allowed through variations of the AP distortion parameters ϵ and α. In the following, we first consider measurements of fσ_{8}, as a derived parameter, and later study the possibility of deriving independent measurements of f and σ_{8}. All those measurements are obtained by marginalizing over the nuisance parameters: p_{n} = (b_{1},b_{2},σ_{v},ϵ,α). The adopted uniform priors on the likelihood parameters are summarized in Table 1 and the full posterior likelihood contours for the cases presented in the next section are given in Appendix B.
Adopted priors on the likelihood parameters.
Fig. 13 fσ_{8} as a function of redshift, showing VIPERS results contrasted with a compilation of recent measurements. The previous results from VVDS (Guzzo et al. 2008), SDSS LRG (Cabré & Gaztañaga 2009; Samushia et al. 2012), WiggleZ (Blake et al. 2012), 6dFGS (Beutler et al. 2012), VIPERS PDR1 (de la Torre et al. 2013), MGS (Howlett et al. 2015), FastSound (Okumura et al. 2016), BOSSLOWZ (GilMarín et al. 2016), BOSSCMASS (GilMarín et al. 2016; Chuang et al. 2016), and VIPERS PDR2 voids (Hawken et al. 2017) are shown with the different symbols (see labels). The solid curve and associated shaded area correspond to the expectations and 68% uncertainty for General Relativity in a ΛCDM background model set to TT+lowP+lensing Planck 2015 predictions (Planck Collaboration XIII 2016). 
7.2.1. fσ_{8} measurements
In our standard configuration, the linear matter power spectrum shape is fixed to the bestfitting ΛCDM model from Planck 2015 TT+lowP+lensing data (Planck Collaboration XIII 2016). AP distortion parameters are set to (ϵ,α) = (0,1) and are not allowed to vary. In this configuration we obtain fσ_{8} values of after marginalizing over other parameters. Associated reduced chisquared values are and respectively. These measurements use both RSD and galaxygalaxy lensing information. It is instructive to see the impact of adding the galaxygalaxy lensing on the measurement of fσ_{8}. Thus if we use the standard RSD approach without including galaxygalaxy lensing information, we obtain with a reduced chisquared value of for both redshifts. In that case, we fixed b_{2} = b_{s2} = b_{3nl} = 0 in the RSD model, as bias nonlinearities are negligible for VIPERS galaxies bias given the minimum scale used in the fit (Pezzotta et al. 2017). Moreover, the shape of nonlinear power spectra in the model is fixed by setting σ_{8} to its fiducial value at the effective redshift of the sample, as is commonly done (e.g. de la Torre et al. 2013). The recovered values and associated errors are very similar to the previous case. We do not find an improvement on fσ_{8} accuracy when galaxygalaxy lensing is included, in fact errors are marginally larger. This can be explained by the lower number of degrees of freedom in the RSDonly case and the significant uncertainty associated with our galaxygalaxy lensing measurements compared to the galaxy clustering ones in the VIPERS fields. In fact the real gain is on constraining f and σ_{8} separately as discussed in Sect. 7.2.4.
7.2.2. Inclusion of AlcockPaczynski distortions
As a robustness test, we relax the assumption on the shape of the linear matter power spectrum. We allow the AP distortion parameters (ϵ,α) to vary, considering flat priors on ϵ,α parameters, extending by ± 0.1 around (ϵ,α) = (0,1). After marginalizing over those parameters as well, we obtain the following fσ_{8} measurements: with reduced chisquared values of and respectively. As expected from the additional degrees of freedom introduced in the likelihood, the marginalized 68% errors on fσ_{8} are increased, although the constraints remain completely compatible with previous measurements when ϵ and α were fixed. This test thus removes any potential concern that our measurements of fσ_{8} might lack robustness though being dependent on the assumption of a ΛCDM expansion history.
7.2.3. Comparison with other measurements
In Fig. 13 we compare our fσ_{8} measurements with previous measurements from the literature, as well as predictions of the standard relativistic model for gravity. Our measurements are consistent with previous measurements at lower or similar redshifts from VVDS (Guzzo et al. 2008), SDSS LRG (Cabré & Gaztañaga 2009; Samushia et al. 2012), WiggleZ (Blake et al. 2012), 6dFGS (Beutler et al. 2012), VIPERS PDR1 (de la Torre et al. 2013), MGS (Howlett et al. 2015), FastSound (Okumura et al. 2016), BOSSLOWZ (GilMarín et al. 2016), and BOSSCMASS (GilMarín et al. 2016; Chuang et al. 2016). In particular, our measurement at z = 0.6 is compatible within 1σ with the WiggleZ z = 0.6 (Blake et al. 2012) and BOSSCMASS z = 0.57 (GilMarín et al. 2016; Chuang et al. 2016) measurements. Our results are also very close to the standard cosmological model predictions: they are consistent within 1σ with General Relativity predictions in a ΛCDM model with cosmological parameters set to Planck CMB results (Planck Collaboration XIII 2016).
These results are part of a combined effort of the VIPERS collaboration to estimate the growth rate of structure from the same data but using different complementary techniques. Specifically, in Pezzotta et al. (2017) we provide a thorough investigation of the performances of different RSD models in configuration space, using a general consistent modelling of nonlinear RSD; in Wilson et al. (in prep.) we use the clipping technique in Fourier space to minimise the impact of nonlinearities; finally in Hawken et al. (2017) we use cosmic voids as RSD tracers. In particular in Hawken et al. (2017), we make use of the void catalogue built from the VIPERS PDR2 data and resulting from the earlier work by Micheletti et al. (2014), to estimate the voidgalaxy crosscorrelation function in redshift space. By modelling its anisotropy we obtain an estimate of fσ_{8} at z = 0.73 and derive a value of , which is lower than those obtained here. However, this technique is still in its infancy, with potential systematic errors not yet fully understood. This and the other VIPERS measurements are all fully compatible within statistical errors. More discussion is presented in the specific papers.
7.2.4. f,b_{1},σ_{8} degeneracy breaking
Fig. 14 Twodimensional marginalized posterior likelihood contours for f and σ_{8} at 0.5 <z < 0.7, showing the impact of the additional galaxygalaxy lensing constraint on the f−σ_{8} degeneracy. The black curve shows the region of constant fσ_{8} associated with Planck Collaboration XIII (2016)ΛCDM + GR bestfit, while the combined (f,σ_{8}) constraint is marked with the star. 
As discussed in Sect. 4, the use of RSD in the galaxy clustering pattern allows a measurement of the parameter combinations fσ_{8}, b_{1}σ_{8}, or β = f/b_{1}. But with the additional constraint of galaxygalaxy lensing, which exhibits different parameter dependencies, we expect to be able to break the f−b_{1}−σ_{8} degeneracy inherent to galaxygalaxy correlations. We investigate this by studying the posterior likelihood contours at 68%,95%,99% for the various pairs of f, b, σ_{8} parameters in our data. This is done for the likelihood analyses presented in the previous sections, i.e. when including or not galaxygalaxy lensing. The posterior likelihood contours are presented in Figs. 14 and 15 for the two considered redshift intervals.
These figures show strong degeneracies in the f−b_{1}, f−σ_{8}, and b_{1}−σ_{8} planes when considering only RSD. In particular, we can see in the f−σ_{8} plane the distribution of the likelihood contours along the regions with constant fσ_{8}, marked with solid and dashed curves in the figures. Now with the inclusion of galaxygalaxy lensing, we can see a shrinking of the contours, in particular along the σ_{8} direction, and to a lesser extent along the b_{1} one. Galaxygalaxy lensing thus effectively provides a strong handle on the σ_{8} parameter. This allows the f−σ_{8} degeneracy to be broken and therefore leads to the possibility of a direct measurement of the growth rate of structure, f. The f−b_{1} degeneracy is also partially broken, even if the effect is milder.
We find that the f−σ_{8} degeneracy breaking is more efficient in the highredshift interval, with measurements of (f,σ_{8}) = (0.93 ± 0.22,0.52 ± 0.06) and (f,σ_{8}) = (0.99 ± 0.19,0.48 ± 0.04) at z = 0.6 and z = 0.86 respectively. These direct measurements of the growth rate of structure and σ_{8} are in agreement within 1σ with Planck ΛCDM + GR predictions, which are (f,σ_{8}) = (0.79,0.60) and (f,σ_{8}) = (0.85,0.53) respectively at z = 0.6 and z = 0.86. Planck ΛCDM + GR predictions are represented with the stars in Figs. 14 and 15. In Fig. 16, we compare our (f,σ_{8}) constraints with those from GilMarín et al. (2017), obtained by combining redshiftspace galaxy power spectrum and bispectrum information in the BOSS survey at z = 0.57. In GilMarín et al. (2017), they use the galaxy bispectrum instead of galaxygalaxy lensing to bring additional constraints on galaxy bias. Although those measurements are quite uncertain, this parameter space and how it can be used as a cosmological model diagnostic, will be very interesting to explore for nextgeneration cosmological surveys, such as Euclid, which will allow a dramatical improvement on such measurement accuracy.
Fig. 16 Joint (f,σ_{8}) constraints at different redshifts. The combined RSD and galaxygalaxy lensing posterior likelihood contours at 1σ and 2σ and those from GilMarín et al. (2017), obtained by combining redshiftspace power spectrum and bispectrum information in the BOSS survey, are presented. The solid curve and associated grey shaded area correspond to the expectations and 68% uncertainty for General Relativity in a ΛCDM background model set to TT+lowP+lensing Planck 2015 predictions (Planck Collaboration XIII 2016), as a function of redshift from z = 2 to z = 0. 
Independent measurements of σ_{8} at different redshifts also carry information about the growth rate of structure. Since σ_{8} grows with time proportionally to the growth factor, the growth rate can be written as dlnσ_{8}/ dlna. In the case of two σ_{8} measurements at a_{1} and a_{2}, as in our analysis, this equation can be approximated through finite difference by (54)By applying this to our σ_{8} measurements we obtain an additional, independent measurement of f = 0.57 ± 0.96 at the mean redshift of z = 0.73. It is clear that this type of measurement is not compelling in our dataset, but can potentially be useful as an additional constraint to be combined with direct measurements in nextgeneration cosmological surveys.
Finally, we notice in Figs. 14 and 15 that the addition of galaxygalaxy lensing constraints significantly modifies the posterior probability distribution function of the linear bias parameter, b_{1}, becoming more compact and skewed towards larger values. This means that adding galaxygalaxy lensing information reduces the uncertainties on b_{1}, and pushes its maximum likelihood value towards values that are in excellent agreement with previous linear bias estimates that are not solely based on twopoint statistics (Di Porto et al. 2016; Cappi et al. 2015; Granett et al. 2015).
7.3. Gravitational slip
In addition to the growth rate of structure, we can measure the gravitational slip parameter E_{G}. This is done by taking the ratio of the measured Υ_{gm} and Υ_{gg}, and multiplying it by β^{1}. The RSD distortion parameter β is estimated from the combined maximum likelihood analysis of the monopole and quadrupole correlation functions (the same as for the RSDonly case presented in Sect. 7.2.1). After marginalizing over nuisance parameters we obtain The 68% error on the E_{G} measurements is obtained by adding in quadrature the fractional error on Υ_{gm}/ Υ_{gg} estimated from mock samples and the fractional error on β^{1}.
The E_{G}(r_{p}) measurements are presented in Fig. 17 for the two redshift intervals under consideration, and compared with the linear predictions for ΛCDM + GR (horizontal line and associated 68% contour). We find that our measurements at z = 0.6 are compatible within 1σ with the standard model, although the central values tend to be slightly lower. We also report in this figure the averaged gravitational slip parameter over the range 3 h^{1} Mpc <r_{p} < 50 h^{1} Mpc, , obtained by Blake et al. (2016) in the similar redshift range 0.43 <z < 0.7. It is represented with a stripe in the figure, with horizontal extent corresponding to the range of r_{p} used to measure and vertical extent showing the ± 1σ error on the measurement. By averaging our E_{G} over 3 h^{1} Mpc <r_{p} < 20 h^{1} Mpc we obtain and . Our results are in good agreement with this measurement and also with that by Pullen et al. (2016) at much higher scales, which also exhibits a slightly lower value compared with ΛCDM + GR prediction. The E_{G} measurements are lower than ΛCDM + GR at r_{p} > 3 h^{1} Mpc but remain within 1−2σ, depending on the scale. At z = 0.86, the agreement with ΛCDM + GR is poorer than at lower redshift.
Fig. 17 Gravitational slip parameter as a function of scale as measured at 0.5 <z < 0.7 (top panel) and 0.7 <z < 1.2 (bottom panel). In both panels, the solid curves and associated shaded areas correspond to the expectations and 68% uncertainties for General Relativity in a ΛCDM background model set to TT+lowP+lensing Planck 2015 predictions (Planck Collaboration XIII 2016). In the top panel, the horizontal stripe shows the averaged over the range 3 h^{1} Mpc <r_{p} < 50 h^{1} Mpc obtained by Blake et al. (2016) at 0.43 <z < 0.7. E_{G} asymptotes to Ω_{m}/f in the standard model, and the simplest way of erasing the modest discrepancy with the model prediction would be to lower the density parameter. 
The origin of the tendency of our E_{G} measurements to be smaller than expected at r_{p} > 3 h^{1} Mpc remains unclear. Particularly given our statistical errors, we have to be cautious in interpreting this trend. In any case, such a result could arise as a result of residual observational systematics or a misinterpretation of the observables, rather than any breakdown of standard gravitational physics. From the construction of E_{G}, these low values of E_{G} seem to be most probably caused by the rather low measured amplitude of Υ_{gm} at r_{p} > 3 h^{1} Mpc, and seen in Fig. 12. If this discrepancy is upheld by further data, one possible interpretation is that weak lensing prefers a lower value of Ω_{m} than that determined by CMB data. It is worth noticing that a similar tension has already been identified in the CFHTLenS cosmic shear analysis of Heymans et al. (2013), as well as in the more recent analysis performed in the KiDS lensing survey (Hildebrandt et al. 2017). It is clear that this point needs to be investigated in detail in the future, in particular in the preparation of nextgeneration very large surveys combining galaxy clustering and weak lensing observables.
8. Conclusion
This paper has presented a combined analysis of RSD and galaxygalaxy lensing in the final VIPERS dataset, making use of complementary data from the CFHTLenS lensing survey over the same area. We have built a consistent theoretical model of the two observables, which includes prescriptions for nonlinear, nonlocal galaxy bias, as well as quasilinear RSD. This model has been shown to enable robust measurements of the growth rate of structure. The model robustness and adopted methodology have been tested by using a series of realistic mock surveys constructed for this purpose.
The main goal of VIPERS has been to provide an accurate measurement of the growth rate of structure using RSD in a redshift regime where the growth is not well determined. With the first data release we were able to provide an initial measurement of fσ_{8} at z = 0.8 (de la Torre et al. 2013). The final dataset increases the survey volume by a factor of 1.6, and by further adding galaxygalaxy lensing information, we have been able to provide new accurate measurements of fσ_{8} at both z = 0.6 and z = 0.86. We have found values of fσ_{8}(z = 0.6) = 0.48 ± 0.12 and fσ_{8}(z = 0.86) = 0.48 ± 0.10, which are consistent with previous measurements at lower or similar redshifts.
The additional galaxygalaxy lensing constraint and the specific treatment of σ_{8} to describe the nonlinearity level of the realspace power spectra entering the model alleviate the degeneracy between the galaxy bias parameter, σ_{8}, and f, and has allowed direct measurements of these two parameters. We have obtained values of (f,σ_{8}) = (0.93 ± 0.22,0.52 ± 0.06) and (f,σ_{8}) = (0.99 ± 0.19,0.48 ± 0.04) at z = 0.6 and z = 0.86, respectively. These measurements put new constraints on gravity at the epoch when the Universe was almost half its present age. Our measurements are statistically consistent with a Universe where the gravitational interactions between structures on cosmological scales can be described by General Relativity, although they are not yet accurate enough to rule out some commonly considered alternatives to General Relativity.
In addition to measuring the growth rate of structure, we have been able to measure the gravitational slip parameter, E_{G}, for the first time at z > 0.6. This quantity, which can be directly constructed from galaxy clustering and galaxygalaxy lensing observables, is sensitive to the growth rate of structure and mean matter density in the Universe. We have obtained averaged values of the gravitational slip parameter of and . Our E_{G} measurements are consistent within 1−2σ, although they exhibit slightly lower values than expected in the standard model for gravity in a ΛCDM background.
Overall, this analysis has demonstrated the importance of the combination of galaxy clustering in redshift space and galaxygalaxy lensing in order to probe the origin of cosmic acceleration. This combination can alleviate the inherent uncertainty related to galaxy bias in RSD analyses and provide new insights into the gravitational physics at work on cosmological scales. This analysis and adopted methodology can be seen as a proofofconcept in the context of the preparation of nextgeneration cosmological surveys such as Euclid (Laureijs et al. 2011), which will allow galaxy clustering and galaxygalaxy lensing to be combined with exquisite precision.
Acknowledgments
We acknowledge the crucial contribution of the ESO staff for the management of service observations. In particular, we are deeply grateful to M. Hilker for his constant help and support of this program. Italian participation to VIPERS has been funded by INAF through PRIN 2008, 2010, 2014 and 2015 programs. L.G. and B.R.G. acknowledge support from the European Research Council through grant No. 291521. O.L.F. acknowledges support from the European Research Council through grant No. 268107. J.A.P. acknowledges support of the European Research Council through the COSFORM ERC Advanced Research Grant (# 670193). G.D.L. acknowledges financial support from the European Research Council through grant No. 202781. R.T. acknowledges financial support from the European Research Council through grant No. 202686. A.P., K.M., and J.K. have been supported by the National Science Centre (grants UMO2012/07/B/ST9/04425 and UMO2013/09/D/ST9/04030). E.B., F.M., and L.M. acknowledge the support from grants ASIINAF I/023/12/0 and PRIN MIUR 20102011. L.M. also acknowledges financial support from PRIN INAF 2012. S.D.L.T., E.J., and M.P. acknowledge the support of the OCEVU Labex (ANR11LABX0060) and the A*MIDEX project (ANR11IDEX000102) funded by the “Investissements d’Avenir” French government program managed by the ANR. Research conducted within the scope of the HECOLS International Associated Laboratory, supported in part by the Polish NCN grant Dec2013/08/M/ST9/00664. T.M. and S.A. acknowledge financial support from the ANR Spin(e) through the French grant ANR13BS050005.
The authors gratefully acknowledge the Gauss Centre for Supercomputing e.V. (www.gausscentre.eu) and the Partnership for Advanced Supercomputing in Europe (PRACE, www.praceri.eu) for funding the MultiDark simulation project by providing computing time on the GCS Supercomputer SuperMUC at Leibniz Supercomputing Centre (LRZ, www.lrz.de).
References
 Alam, S., Ata, M., Bailey, S., et al. 2017, MNRAS, 470, 2617 [NASA ADS] [CrossRef] [Google Scholar]
 Alcock, C., & Paczynski, B. 1979, Nature, 281, 358 [NASA ADS] [CrossRef] [Google Scholar]
 Anderson, L., Aubourg, É., Bailey, S., et al. 2014, MNRAS, 441, 24 [NASA ADS] [CrossRef] [Google Scholar]
 Baldauf, T., Smith, R. E., Seljak, U., & Mandelbaum, R. 2010, Phys. Rev. D, 81, 063531 [NASA ADS] [CrossRef] [Google Scholar]
 Ballinger, W. E., Peacock, J. A., & Heavens, A. F. 1996, MNRAS, 282, 877 [NASA ADS] [CrossRef] [Google Scholar]
 Benjamin, J., Van Waerbeke, L., Heymans, C., et al. 2013, MNRAS, 431, 1547 [NASA ADS] [CrossRef] [Google Scholar]
 Bernardeau, F., Colombi, S., Gaztañaga, E., & Scoccimarro, R. 2002, Phys. Rep., 367, 1 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
 Beutler, F., Blake, C., Colless, M., et al. 2011, MNRAS, 416, 3017 [NASA ADS] [CrossRef] [Google Scholar]
 Beutler, F., Blake, C., Colless, M., et al. 2012, MNRAS, 423, 3430 [NASA ADS] [CrossRef] [Google Scholar]
 Beutler, F., Saito, S., Seo, H.J., et al. 2014, MNRAS, 443, 1065 [NASA ADS] [CrossRef] [Google Scholar]
 Blake, C., Brough, S., Colless, M., et al. 2012, MNRAS, 425, 405 [NASA ADS] [CrossRef] [Google Scholar]
 Blake, C., Joudaki, S., Heymans, C., et al. 2016, MNRAS, 456, 2806 [NASA ADS] [CrossRef] [Google Scholar]
 Cabré, A., & Gaztañaga, E. 2009, MNRAS, 393, 1183 [NASA ADS] [CrossRef] [Google Scholar]
 Cacciato, M., van den Bosch, F. C., More, S., Mo, H., & Yang, X. 2013, MNRAS, 430, 767 [NASA ADS] [CrossRef] [Google Scholar]
 Cappi, A., Marulli, F., Bel, J., et al. 2015, A&A, 579, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chan, K. C., Scoccimarro, R., & Sheth, R. K. 2012, Phys. Rev. D, 85, 083509 [NASA ADS] [CrossRef] [Google Scholar]
 Chuang, C.H., Prada, F., PellejeroIbanez, M., et al. 2016, MNRAS, 461, 3781 [NASA ADS] [CrossRef] [Google Scholar]
 Cole, S. 2011, MNRAS, 416, 739 [NASA ADS] [CrossRef] [Google Scholar]
 Cole, S., Percival, W. J., Peacock, J. A., et al. 2005, MNRAS, 362, 505 [NASA ADS] [CrossRef] [Google Scholar]
 Colless, M., Dalton, G., Maddox, S., et al. 2001, MNRAS, 328, 1039 [NASA ADS] [CrossRef] [Google Scholar]
 Colless, M., Peterson, B. A., Jackson, C., et al. 2003, ArXiv eprints [arXiv:astroph/0306581] [Google Scholar]
 Coupon, J., Arnouts, S., van Waerbeke, L., et al. 2015, MNRAS, 449, 1352 [NASA ADS] [CrossRef] [Google Scholar]
 de la Torre, S., & Guzzo, L. 2012, MNRAS, 427, 327 [NASA ADS] [CrossRef] [Google Scholar]
 de la Torre, S., & Peacock, J. A. 2013, MNRAS, 435, 743 [NASA ADS] [CrossRef] [Google Scholar]
 de la Torre, S., Guzzo, L., Peacock, J. A., et al. 2013, A&A, 557, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Di Porto, C., Branchini, E., Bel, J., et al. 2016, A&A, 594, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dodelson, S., & Schneider, M. D. 2013, Phys. Rev. D, 88, 063537 [NASA ADS] [CrossRef] [Google Scholar]
 Eisenstein, D. J., Zehavi, I., Hogg, D. W., et al. 2005, ApJ, 633, 560 [NASA ADS] [CrossRef] [Google Scholar]
 Erben, T., Hildebrandt, H., Miller, L., et al. 2013, MNRAS, 433, 2545 [NASA ADS] [CrossRef] [Google Scholar]
 Fisher, K. B., Davis, M., Strauss, M. A., Yahil, A., & Huchra, J. P. 1994, MNRAS, 267, 927 [NASA ADS] [CrossRef] [Google Scholar]
 Garilli, B., Guzzo, L., Scodeggio, M., et al. 2014, A&A, 562, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 GilMarín, H., Wagner, C., Noreña, J., Verde, L., & Percival, W. 2014, J. Cosmology Astropart. Phys., 12, 029 [CrossRef] [Google Scholar]
 GilMarín, H., Percival, W. J., Brownstein, J. R., et al. 2016, MNRAS, 460, 4188 [NASA ADS] [CrossRef] [Google Scholar]
 GilMarín, H., Percival, W. J., Verde, L., et al. 2017, MNRAS, 465, 1757 [NASA ADS] [CrossRef] [Google Scholar]
 Giocoli, C., Jullo, E., Metcalf, R. B., et al. 2016, MNRAS, 461, 209 [NASA ADS] [CrossRef] [Google Scholar]
 Goranova, Y., Hudelot, P., Magnard, F., et al. 2009, The CFHTLS T0006 Release, http://terapix.iap.fr/cplt/table_syn_T0006.html [Google Scholar]
 Goroff, M. H., Grinstein, B., Rey, S.J., & Wise, M. B. 1986, ApJ, 311, 6 [NASA ADS] [CrossRef] [Google Scholar]
 Granett, B. R., Branchini, E., Guzzo, L., et al. 2015, A&A, 583, A61 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Guzzo, L., Pierleoni, M., Meneux, B., et al. 2008, Nature, 451, 541 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Guzzo, L., Scodeggio, M., Garilli, B., et al. 2014, A&A, 566, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hamilton, A. J. S. 1993, ApJ, 417, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Hartlap, J., Simon, P., & Schneider, P. 2007, A&A, 464, 399 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hawken, A. J., Granett, B. R., Iovino, A., et al. 2017, A&A, 607, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Heymans, C., Van Waerbeke, L., Miller, L., et al. 2012, MNRAS, 427, 146 [NASA ADS] [CrossRef] [Google Scholar]
 Heymans, C., Grocutt, E., Heavens, A., et al. 2013, MNRAS, 432, 2433 [NASA ADS] [CrossRef] [Google Scholar]
 Hildebrandt, H., Erben, T., Kuijken, K., et al. 2012, MNRAS, 421, 2355 [NASA ADS] [CrossRef] [Google Scholar]
 Hildebrandt, H., Viola, M., Heymans, C., et al. 2017, MNRAS, 465, 1454 [Google Scholar]
 Howlett, C., Ross, A. J., Samushia, L., Percival, W. J., & Manera, M. 2015, MNRAS, 449, 848 [NASA ADS] [CrossRef] [Google Scholar]
 Kaiser, N. 1987, MNRAS, 227, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Kaufman, C. G., Schervish, M. J., & Nychka, D. W. 2008, J. Am. Statist. Ass., 103, 1545 [Google Scholar]
 Klypin, A., Yepes, G., Gottlöber, S., Prada, F., & Heß, S. 2016, MNRAS, 457, 4340 [NASA ADS] [CrossRef] [Google Scholar]
 Landy, S. D., & Szalay, A. S. 1993, ApJ, 412, 64 [NASA ADS] [CrossRef] [Google Scholar]
 Laureijs, R., Amiaux, J., Arduini, S., et al. 2011, ArXiv eprints [arXiv:1110.3193] [Google Scholar]
 Le Fèvre, O., Saisse, M., Mancini, D., et al. 2003, in Proc. SPIE 4841, eds. M. Iye, & A. F. M. Moorwood, 1670 [Google Scholar]
 Lesgourgues, J. 2011, ArXiv eprints [arXiv:1104.2932] [Google Scholar]
 Mandelbaum, R., Slosar, A., Baldauf, T., et al. 2013, MNRAS, 432, 1544 [NASA ADS] [CrossRef] [Google Scholar]
 Marulli, F., Bolzonella, M., Branchini, E., et al. 2013, A&A, 557, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Matsubara, T., & Suto, Y. 1996, ApJ, 470, L1 [NASA ADS] [CrossRef] [Google Scholar]
 McDonald, P., & Roy, A. 2009, J. Cosmology Astropart. Phys., 8, 020 [Google Scholar]
 Metropolis, N., Rosenbluth, A. W., Rosenbluth, M. N., Teller, A. H., & Teller, E. 1953, J. Chem. Phys., 21, 1087 [NASA ADS] [CrossRef] [Google Scholar]
 Micheletti, D., Iovino, A., Hawken, A. J., et al. 2014, A&A, 570, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Miller, L., Heymans, C., Kitching, T. D., et al. 2013, MNRAS, 429, 2858 [Google Scholar]
 Miyatake, H., More, S., Mandelbaum, R., et al. 2015, ApJ, 806, 1 [NASA ADS] [CrossRef] [Google Scholar]
 More, S., Miyatake, H., Mandelbaum, R., et al. 2015, ApJ, 806, 2 [NASA ADS] [CrossRef] [Google Scholar]
 Nakajima, R., Mandelbaum, R., Seljak, U., et al. 2012, MNRAS, 420, 3240 [NASA ADS] [Google Scholar]
 Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 563 [NASA ADS] [CrossRef] [Google Scholar]
 Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 [NASA ADS] [CrossRef] [Google Scholar]
 Okumura, T., Hikage, C., Totani, T., et al. 2016, PASJ, 68, 38 [NASA ADS] [Google Scholar]
 Paz, D. J., & Sánchez, A. G. 2015, MNRAS, 454, 4326 [NASA ADS] [CrossRef] [Google Scholar]
 Peacock, J. A., Cole, S., Norberg, P., et al. 2001, Nature, 410, 169 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Percival, W. J., Reid, B. A., Eisenstein, D. J., et al. 2010, MNRAS, 401, 2148 [NASA ADS] [CrossRef] [Google Scholar]
 Percival, W. J., Ross, A. J., Sánchez, A. G., et al. 2014, MNRAS, 439, 2531 [NASA ADS] [CrossRef] [Google Scholar]
 Pezzotta, A., de la Torre, S., Bel, J., et al. 2017, A&A, 604, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XIII. 2016, A&A, 594, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pope, A. C., & Szapudi, I. 2008, MNRAS, 389, 766 [NASA ADS] [CrossRef] [Google Scholar]
 Pullen, A. R., Alam, S., He, S., & Ho, S. 2016, MNRAS, 460, 4098 [NASA ADS] [CrossRef] [Google Scholar]
 Reyes, R., Mandelbaum, R., Seljak, U., et al. 2010, Nature, 464, 256 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Rota, S., Granett, B. R., Bel, J., et al. 2017, A&A, 601, A144 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Saito, S., Baldauf, T., Vlah, Z., et al. 2014, Phys. Rev. D, 90, 123522 [NASA ADS] [CrossRef] [Google Scholar]
 Samushia, L., Percival, W. J., & Raccanelli, A. 2012, MNRAS, 420, 2102 [NASA ADS] [CrossRef] [Google Scholar]
 Samushia, L., Reid, B. A., White, M., et al. 2014, MNRAS, 439, 3504 [NASA ADS] [CrossRef] [Google Scholar]
 Scoccimarro, R., Couchman, H. M. P., & Frieman, J. A. 1999, ApJ, 517, 531 [NASA ADS] [CrossRef] [Google Scholar]
 Scodeggio, M., Franzetti, P., Garilli, B., LeFèvre, O., & Guzzo, L. 2009, The Messenger, 135, 13 [NASA ADS] [Google Scholar]
 Scodeggio, M., Guzzo, L., Garilli, B., et al. 2017, A&A, in press, DOI: 10.1051/00046361/201630114 [Google Scholar]
 Simpson, F., Heymans, C., Parkinson, D., et al. 2013, MNRAS, 429, 2249 [NASA ADS] [CrossRef] [Google Scholar]
 Smith, R. E., Peacock, J. A., Jenkins, A., et al. 2003, MNRAS, 341, 1311 [NASA ADS] [CrossRef] [Google Scholar]
 Takahashi, R., Sato, M., Nishimichi, T., Taruya, A., & Oguri, M. 2012, ApJ, 761, 152 [NASA ADS] [CrossRef] [Google Scholar]
 Taruya, A., Nishimichi, T., & Saito, S. 2010, Phys. Rev. D, 82, 063522 [NASA ADS] [CrossRef] [Google Scholar]
 Taylor, A., & Joachimi, B. 2014, MNRAS, 442, 2728 [NASA ADS] [CrossRef] [Google Scholar]
 Tegmark, M., Blanton, M. R., Strauss, M. A., et al. 2004, ApJ, 606, 702 [NASA ADS] [CrossRef] [Google Scholar]
 Tinker, J., Kravtsov, A. V., Klypin, A., et al. 2008, ApJ, 688, 709 [NASA ADS] [CrossRef] [Google Scholar]
 Tinker, J. L., Robertson, B. E., Kravtsov, A. V., et al. 2010, ApJ, 724, 878 [NASA ADS] [CrossRef] [Google Scholar]
 Velander, M., van Uitert, E., Hoekstra, H., et al. 2014, MNRAS, 437, 2111 [NASA ADS] [CrossRef] [Google Scholar]
 Xu, X., Cuesta, A. J., Padmanabhan, N., Eisenstein, D. J., & McBride, C. K. 2013, MNRAS, 431, 2834 [NASA ADS] [CrossRef] [Google Scholar]
 Zhang, P., Liguori, M., Bean, R., & Dodelson, S. 2007, Phys. Rev. Lett., 99, 141302 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
Appendix A: Theoretical power spectra for biased tracers
This appendix presents the models describing the realspace galaxygalaxy, galaxyvelocity divergence, and galaxymatter power spectra, which enter the modelling of RSD and galaxygalaxy lensing. We adopt the nonlinear nonlocal bias model of McDonald & Roy (2009) that relates the galaxy overdensity δ_{g} and matter overdensity δ as: (A.1)In this equation, b_{1} and b_{2} are the linear and secondorder nonlinear bias terms, b_{s2} the nonlocal bias term, s is the tidal tensor term from which nonlocality originates. The σ^{2} and ⟨ s^{2} ⟩ terms ensure the condition ⟨ δ_{g} ⟩ = 0. From the bias model of Eq. (A.1) one can derive the following power spectra for galaxygalaxy, galaxyvelocity divergence (θ), and galaxymatter correlations: where (e.g. Beutler et al. 2014; GilMarín et al. 2014) In the above equations, S_{2}, F_{2}, G_{2} perturbation theory kernels are defined by (e.g. Goroff et al. 1986; Bernardeau et al. 2002)
Appendix B: Posterior likelihood contours
In this appendix are provided the posterior likelihood contours of all pairs of parameters appearing in the likelihood analyses presented in Sect. 7.2.1. Figure B.1 shows the posterior likelihood contours in the case of the RSDonly analysis, while Fig. B.2 in the case where RSD and galaxygalaxy lensing are combined. In both figures, the three types of shaded regions in each subpanel correspond to the posterior likelihood contours at 68%, 95%, and 99%.
Fig. B.1 Posterior likelihood contours for f, σ_{8}, b_{1}, and σ_{v} parameters at z = 0.6 (left panel) and z = 0.86 (right panel) in the case where RSD are considered alone (see Sect. 7.2.1). 
Fig. B.2 Posterior likelihood contours for f, σ_{8}, b_{1}, b_{2}, and σ_{v} parameters at z = 0.6 (left panel) and z = 0.86 (right panel) in the case where RSD and galaxygalaxy lensing are combined (see Sect. 7.2.1). 
All Tables
All Figures
Fig. 1 Number densities of VIPERS galaxies in the individual W1 and W4 fields and of CFHTLenS/VIPERS photometric redshift galaxies, as a function of redshift. The number densities of VIPERS galaxies are corrected for the survey incompleteness by weighting each galaxy in the counts by its associated inverse completeness weight w^{C}. The solid curve corresponds to the model n(z) used in the analysis. It was obtained by randomly sampling galaxy redshifts within their V_{max} (see text for details). 

In the text 
Fig. 2 Relative difference between various estimates of ΔΣ_{gm}, based on different assumptions for source and lens redshifts, and the fiducial estimate in the data at 0.5 <z < 0.7. The quantity shown in the figure is as a function of the projected separation r_{p}. The fiducial estimate is that obtained by using Eq. (13), which includes the individual redshift probability distribution function p_{s}(z) of the sources, and for the lenses, the VIPERS spectroscopic redshift (z_{spec}) when available or the CFHTLenS maximum likelihood photometric redshift (z_{phot}) otherwise (see text). It corresponds to the adopted estimate for the analysis. The grey shaded area represents the relative statistical error expected in the survey. 

In the text 
Fig. 3 Probability distribution function of redshift errors at 0.5 <z < 0.7 and 0.7 <z < 1.2 in the VIPERS data. This is obtained from the redshift differences of reobserved galaxies, for which there are two independent redshift measurements. The dotted and dashed curves are bestfitting Gaussians for the redshift intervals 0.5 <z < 0.7 and 0.7 <z < 1.2 respectively. 

In the text 
Fig. 4 Relative systematic errors on the correlation function monopole (top panel) and quadrupole (bottom panel) at 0.5 <z < 0.7 and effects of target sampling rate (TSR) and angular pair weighting (w^{A}) corrections. The grey shaded areas represent the relative statistical error expected in the survey, while light grey band mark ± 1% relative uncertainties for reference. 

In the text 
Fig. 5 Same as in Fig. 4 but for the redshift interval 0.7 <z < 1.2. 

In the text 
Fig. 6 Relative systematic error on fσ_{8} at 0.5 <z < 0.7 as a function of s_{min}, for different values of r_{0} (r_{0} = 1 h^{1} Mpc and r_{0} = 1.5 h^{1} Mpc) and when including or not redshift error. The error bars represent the relative statistical error associated to analysing the mean mock predictions. The shaded area shows the 1σ confidence region associated with the relative statistical error expected in VIPERS. The squares and triangles are artificially shifted along s_{min} axis to improve the clarity of the figure. 

In the text 
Fig. 7 Comparison of differential excess surface density Δ_{gm} and annular differential excess surface density Υ_{gm} predictions in the mocks (points and shaded regions) and by our theoretical model (curves) at 0.5 <z < 0.7. The mock predictions correspond to the mean signal among the mock realizations (points) and its associated 1σ error (shaded region). The curves show the theoretical predictions for the fiducial parameters of the mocks, varying only the b_{2} parameter as labeled. 

In the text 
Fig. 8 Top panel: recovered errors on fσ_{8} normalized to that obtained without tapering (equivalent to applying tapering with T_{p} = ∞), as a function of the tapering scale T_{p} used in the estimation of the precision matrix. This is obtained from the mocks at 0.5 <z < 0.7. Bottom panel: recovered maximum likelihood values for fσ_{8} and associated 1σ error as a function of the tapering scale T_{p}. 

In the text 
Fig. 9 Correlation matrix (left panel) and normalized precision matrix (right panel) for galaxy clustering and galaxygalaxy lensing data in the redshift interval 0.5 <z < 0.7. These are defined as C_{ij}/ and / respectively, where C_{ij} and refer to covariance and precision matrix elements respectively. In both panels, the upper triangular matrix represents the case without tapering, while the lower panel the case with tapering. The precision matrix is normalized such that diagonal elements are unity. 

In the text 
Fig. 10 Anisotropic correlation functions of VIPERS galaxies at 0.5 <z < 0.7 (top panel) and 0.7 <z < 1.2 (bottom panel) as a function parallel and transverse to the lineofsight separations. 

In the text 
Fig. 11 Monopole (circles) and quadrupole (triangles) correlation functions of VIPERS galaxies at 0.5 <z < 0.7 (top panel) and 0.7 <z < 1.2 (bottom panel). Solid curves correspond to individual mock measurements. 

In the text 
Fig. 12 Differential excess surface density (circles) and annular differential excess surface density (triangles) at 0.5 <z < 0.7 (top panel) and 0.7 <z < 1.2 (bottom panel). Solid curves correspond to individual mock measurements. 

In the text 
Fig. 13 fσ_{8} as a function of redshift, showing VIPERS results contrasted with a compilation of recent measurements. The previous results from VVDS (Guzzo et al. 2008), SDSS LRG (Cabré & Gaztañaga 2009; Samushia et al. 2012), WiggleZ (Blake et al. 2012), 6dFGS (Beutler et al. 2012), VIPERS PDR1 (de la Torre et al. 2013), MGS (Howlett et al. 2015), FastSound (Okumura et al. 2016), BOSSLOWZ (GilMarín et al. 2016), BOSSCMASS (GilMarín et al. 2016; Chuang et al. 2016), and VIPERS PDR2 voids (Hawken et al. 2017) are shown with the different symbols (see labels). The solid curve and associated shaded area correspond to the expectations and 68% uncertainty for General Relativity in a ΛCDM background model set to TT+lowP+lensing Planck 2015 predictions (Planck Collaboration XIII 2016). 

In the text 
Fig. 14 Twodimensional marginalized posterior likelihood contours for f and σ_{8} at 0.5 <z < 0.7, showing the impact of the additional galaxygalaxy lensing constraint on the f−σ_{8} degeneracy. The black curve shows the region of constant fσ_{8} associated with Planck Collaboration XIII (2016)ΛCDM + GR bestfit, while the combined (f,σ_{8}) constraint is marked with the star. 

In the text 
Fig. 15 Same as Fig. 14 but for the redshift interval 0.7 <z < 1.2. 

In the text 
Fig. 16 Joint (f,σ_{8}) constraints at different redshifts. The combined RSD and galaxygalaxy lensing posterior likelihood contours at 1σ and 2σ and those from GilMarín et al. (2017), obtained by combining redshiftspace power spectrum and bispectrum information in the BOSS survey, are presented. The solid curve and associated grey shaded area correspond to the expectations and 68% uncertainty for General Relativity in a ΛCDM background model set to TT+lowP+lensing Planck 2015 predictions (Planck Collaboration XIII 2016), as a function of redshift from z = 2 to z = 0. 

In the text 
Fig. 17 Gravitational slip parameter as a function of scale as measured at 0.5 <z < 0.7 (top panel) and 0.7 <z < 1.2 (bottom panel). In both panels, the solid curves and associated shaded areas correspond to the expectations and 68% uncertainties for General Relativity in a ΛCDM background model set to TT+lowP+lensing Planck 2015 predictions (Planck Collaboration XIII 2016). In the top panel, the horizontal stripe shows the averaged over the range 3 h^{1} Mpc <r_{p} < 50 h^{1} Mpc obtained by Blake et al. (2016) at 0.43 <z < 0.7. E_{G} asymptotes to Ω_{m}/f in the standard model, and the simplest way of erasing the modest discrepancy with the model prediction would be to lower the density parameter. 

In the text 
Fig. B.1 Posterior likelihood contours for f, σ_{8}, b_{1}, and σ_{v} parameters at z = 0.6 (left panel) and z = 0.86 (right panel) in the case where RSD are considered alone (see Sect. 7.2.1). 

In the text 
Fig. B.2 Posterior likelihood contours for f, σ_{8}, b_{1}, b_{2}, and σ_{v} parameters at z = 0.6 (left panel) and z = 0.86 (right panel) in the case where RSD and galaxygalaxy lensing are combined (see Sect. 7.2.1). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.