MNRAS 000,1–15 (2020) Preprint 23 December 2020 Compiled using MNRAS LATEX style file v3.0 Testing General Relativity on cosmological scales at redshift I ∼ 1.5 with quasar and CMB lensing Yucheng Zhang ,1¢ Anthony R. Pullen,1,2 Shadab Alam ,3 Sukhdeep Singh,4 Etienne Burtin,5 Chia-Hsun Chuang,6 Jiamin Hou,7 Brad W. Lyke,8 Adam D. Myers,8 Richard Neveux,5 Ashley J. Ross,9 Graziano Rossi,10 Cheng Zhao11 1 Center for Cosmology and Particle Physics, Department of Physics, New York University, 726 Broadway, New York, NY 10003, USA 2 Center for Computational Astrophysics, Flatiron Institute, New York, NY 10010, USA 3 Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh, EH9 3HJ, UK 4 Berkeley Center for Cosmological Physics, University of California, Berkeley, CA 94720, USA 5 IRFU,CEA, Université Paris-Saclay, F-91191 Gif-sur-Yvette, France 6 Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, 452 Lomita Mall, Stanford, CA 94305, USA 7 Max-Planck-Institut für Extraterrestrische Physik, Postfach 1312, Giessenbachstrasse 1, 85748 Garching bei München, Germany 8 University of Wyoming, 1000 E. University Ave., Laramie, WY 82071, USA 9 Center for Cosmology and Astro-Particle Physics, Ohio State University, Columbus, Ohio, USA 10 Department of Physics and Astronomy, Sejong University, Seoul, 143-747, Korea 11 Institute of Physics, Laboratory of Astrophysics, École Polytechnique Fédérale de Lausanne (EPFL), Observatoire de Sauverny, 1290 Versoix, Switzerland Accepted XXX. Received YYY; in original form ZZZ ABSTRACT We test general relativity (GR) at the effective redshift I¯ ∼ 1.5 by estimating the statistic 퐸퐺, a probe of gravity, on cosmological scales 19 − 190 ℎ−1Mpc. This is the highest redshift and largest scale estimation of 퐸퐺 so far. We use the quasar sample with redshifts 0.8 < I < 2.2 from Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey Data Release 16 as the large-scale structure (LSS) tracer, for which the angular power spectrum @@ 퐶ℓ and the redshift-space distortion parameter V are estimated. By cross-correlating with the Planck 2018 cosmic microwave background (CMB) lensing map, we detect the angular cross- ^@ power spectrum 퐶ℓ signal at 12 f significance. Both jackknife resampling and simulations are used to estimate the covariance matrix (CM) of 퐸퐺 at five bins covering different scales, with the later preferred for its better constraints on the covariances. We find 퐸퐺 estimates agree with the GR prediction at 1 f level over all these scales. With the CM estimated with 300 simulations, we report a best-fitting scale-averaged estimate of 퐸퐺 ¹I¯º = 0.30 ± 0.05, which is GR ¹ º in line with the GR prediction 퐸퐺 I¯ = 0.33 with Planck 2018 CMB+BAO matter density fraction Ωm = 0.31. The statistical errors of 퐸퐺 with future LSS surveys at similar redshifts will be reduced by an order of magnitude, which makes it possible to constrain modified gravity models. Key words: cosmology: theory – cosmology: observations – large-scale structure of Universe – gravitation – gravitational lensing: weak – cosmic background radiation arXiv:2007.12607v2 [astro-ph.CO] 22 Dec 2020 1 INTRODUCTION explain the cosmic expansion and acceleration (see e.g. Silvestri & Trodden 2009, and references therein), among which Λ-cold dark The expansion of the universe was first discovered by measuring matter (ΛCDM) has been regarded as the standard model for its the redshifts and relative distances of galaxies (Hubble 1929). One simplicity and success in explaining a wide range of cosmological of the milestones in cosmology in the past decades has been the observations, including the CMB (cosmic microwave background) detection of a negative deceleration parameter and hence the ac- surveys (e.g. Planck Collaboration VI 2018) and galaxy redshift sur- celerated expansion of the Universe at late times from supernovae veys (e.g. Alam et al. 2017b). ΛCDM takes general relativity (GR) observations (Riess et al. 1998; Perlmutter et al. 1999). Many the- as the true theory for gravity on both galactic and cosmological oretical models of cosmology and gravity have been proposed to scales, and assumes the existence of the cosmological constant (Λ), a special form of dark energy (DE) whose spatially uniform energy density does not evolve with cosmic expansion, and CDM, along ¢ E-mail: [email protected] © 2020 The Authors 2 Y. Zhang et al. with ordinary (baryonic) matter. Although the expansion history can order to do certain estimations and generate the simulations needed. be well described by ΛCDM-GR by fine-tuning the relative density So for now it is very difficult to design one blind test for various ratios of the energy components, the nature of dark matter (DM) and gravity models. To do a rigorous estimation of 퐸퐺 based on other DE are not well understood and their properties are hard to detect MG models, the corresponding changes have to be made for either with observations. On the other hand, some modified gravity (MG) simulations or analytic calculations (see e.g. Hojjati et al. 2011). models (see e.g. Carroll et al. 2005; Sotiriou & Faraoni 2010; Dvali The paper is organized as follows. In Section2, we review the et al. 2000), which can predict the same expansion history of the 퐸퐺 theory and describe the estimator we use. The quasar and CMB universe as ΛCDM-GR with or completely without assuming the data, simulations and jackknife resampling for the estimation of existence of DE, have been developed to challenge GR as the true covariance matrices are described in Section3. Section4 includes theory for gravity on cosmological scales. There have been some analytic models, estimators, systematics and calibrations for the great reviews of the two approaches, see e.g. Peebles & Ratra(2003) angular power spectra. Section5 describes our estimation of the for the cosmological constant and DE, Clifton et al.(2012) for MG, quasar 2-point correlation function and the maximum likelihood and Joyce et al.(2016) for a comparison. fitting of the redshift-space distortion (RSD) parameter. We present Despite the degeneracy between ΛCDM-GR and MG models all the estimates and our final results in Section6 and conclude in in explaining the cosmic expansion, their predictions of the growth Section7. of the DM large scale structure (LSS) are usually distinguishable. For our self-consistency test of GR, wherever needed, we Combining the gravitational lensing r2 ¹Ψ − Φº and the divergence assume a flat ΛCDM fiducial cosmology with Planck 2018 of the peculiar velocity \, Zhang et al.(2007) proposed a statistic CMB¸BAO parameters (Planck Collaboration VI 2018): Ωm = 2 2 퐸퐺 as a function of redshift and scale, to probe gravity on cosmo- 0.3111 ± 0.0056, Ωcℎ = 0.11933 ± 0.00091, Ωbℎ = 0.02242 ± logical scales. Lensing is related to the underlying matter overden- 0.00014, =s = 0.9665 ± 0.0038, 퐻0 = 67.66 ± 0.42, and f8 = sity X through the Poisson equation which depends on the gravity 0.8102 ± 0.0060. model (see e.g. Hojjati et al. 2011). On linear scales, \ = − 5 X, where 5 is the linear growth rate. In real surveys, instead of the DM field, the direct observables are the LSS tracers, e.g. galaxies or quasars. The distribution of these tracers is connected to the un- 2 퐸 FORMALISM AND ESTIMATOR derlying matter perturbation field with the clustering bias 1, which 퐺 varies with the physical properties of the tracers that are targeted in In this section, we briefly review the 퐸퐺 theory and describe the a particular survey. Defined as the ratio between r2 ¹Ψ − Φº and \, estimator used in this work. We assume a flat Universe described 퐸퐺 has the advantage of being independent of 1 and the variance by the perturbed Friedmann-Robertson-Walker (FRW) metric in of the matter density field f8. Interested readers may refer to Ishak conformal Newtonian gauge, (2019) for a comprehensive review on various cosmological tests of 2 h 2 2i GR, including the 퐸퐺 statistic. 3B = 0¹gº ¹1 ¸ 2Ψº3g − ¹1 ¸ 2Φº3G , (1) The estimation of 퐸퐺 requires data from both gravitational lensing and redshift surveys. Accurate estimates of tracers’ red- where Ψ and Φ are the scalar perturbations to the time and spatial shifts are necessary in order to do the 3-D clustering analysis, from components of the metric. The statistic 퐸퐺 is defined in Fourier which the growth of the structure can be probed. Thus spectro- Space (Zhang et al. 2007) as scopic redshift surveys are usually preferred. For photometric sur- " 2 # veys, Giannantonio et al.(2016) proposed a statistic 퐷 , which does r ¹Ψ − Φº 퐺 퐸퐺 ¹:,Iº = not require the estimation of the growth rate. However, this quan- −3퐻2 ¹1 ¸ Iº\ 0 : tity cannot be directly used to discriminate GR and MG models. (2) :2 ¹Ψ − Φº Using galaxy-galaxy lensing and luminous red galaxies (LRGs), = , −1 3퐻2 ¹1 ¸ Iº\ 퐸퐺 has been measured over scales . 70 ℎ Mpc at redshifts in 0 0.2 < I < 0.6 (Reyes et al. 2010; Blake et al. 2015; de la Torre et al. where 퐻 is the Hubble constant and \ = r·v/퐻¹Iº is the divergence 2017; Alam et al. 2017a; Amon et al. 2018; Singh et al. 2018; Blake 0 of the comoving peculiar velocity field. In linear perturbation theory, et al. 2020). Besides tracing the lensing signal with background \ = − 5 X, where 5 is the linear growth rate and X is the matter galaxies, Pullen et al.(2015) proposed to use the cosmic microwave perturbation. For GR, assuming no anisotropic stress (Φ = −Ψ) and background (CMB) lensing map, which allows the estimation of 2 2 using Poisson equation r Ψ = 4c퐺0 d<X, we have 퐸퐺 at higher redshifts and larger scales (Pullen et al.