We introduce a new version of the Integrated Software in Testing General Relativity (ISiTGR) which is a patch to the software CAMB and CosmoMC. ISiTGR is intended to test deviations from GR at cosmological scales using cosmological data sets. While doing so, it allows for various extensions to the standard flat $\Lambda$CDM model. In this new release, we have support for the following: 1) dynamical dark energy parametrizations with a constant or time-dependent equation of state; 2) a consistent implementation of anisotropic shear to model massive neutrinos throughout the full formalism; 3) multiple commonly-used parametrizations of modified growth (MG) parameters; 4) functional, binned and hybrid time- and scale-dependencies for all MG parameters; 5) spatially flat or curved backgrounds. ISiTGR is designed to allow cosmological analyses to take full advantage of ongoing and future surveys to test simultaneously or separately various extensions to the standard model. We describe here the formalism and its implementation in the CMB code, the Integrated Sachs-Wolfe (ISW) effect, and the 3x2 point statistics. Next, we apply ISiTGR to current data sets from Planck-2018, Planck-2015, Dark Energy Survey YR1 release, Baryonic Acoustic Oscillations (BAO), Redshift Space Distortions (BAO/RSD) from the BOSS Data Release 12, the 6DF Galaxy Survey and the SDSS Data Release 7 Main Galaxy Sample, and Supernova from the Pantheon compilation, joint SNLS/SDSS data analysis and the Hubble Space Telescope. We derive constraints on MG parameters for various combinations of the five features above and find that GR is consistent with current data sets in all cases. The code is made publicly available at \url{https://github.com/mishakb/ISiTGR}.

The clustering ratio is defined as the ratio between the correlation function and the variance of the smoothed overdensity field. In LCDM cosmologies not accounting for massive neutrinos, it has already been proved to be independent from bias and redshift space distortions on a range of linear scales. It therefore allows for a direct comparison of measurements (from galaxies in redshift space) to predictions (for matter in real space). In this paper we first extend the applicability of such properties of the clustering ratio to cosmologies that include massive neutrinos, by performing tests against simulated data. We then investigate the constraining power of the clustering ratio when cosmological parameters such as the total neutrino mass and the equation of state of dark energy are left free. We analyse the joint posterior distribution of the parameters that must satisfy, at the same time, the measurements of the galaxy clustering ratio in the SDSS DR12, and the angular power spectrum of temperature and polarization anisotropies of the CMB measured by the Planck satellite. We find the clustering ratio to be very sensitive to the CDM density parameter, but not very much so to the total neutrino mass. Lastly, we forecast the constraining power the clustering ratio will achieve with forthcoming surveys, predicting the amplitude of its errors in a Euclid-like galaxy survey. In this case, we find it is expected to improve the constraint at 95% level on the CDM density by 40% and on the total neutrino mass by 14%.

We present constraints on testing general relativity (GR) at cosmological scales using recent data sets and assess the impact of galaxy intrinsic alignment (IA) in the CFHTLenS lensing data on those constraints. We consider CMB temperature data from Planck, the galaxy power spectrum from WiggleZ, weak lensing tomography from the CFHTLenS, ISW-galaxy cross correlations, and BAO data from 6dF, SDSS DR7, and BOSS DR9. We use a parameterization of the modified gravity (MG) that is binned in redshift and scale, a parameterization that evolves monotonically in scale but is binned in redshift, and a functional parameterization that evolves only in redshift. We present the results in terms of the MG parameters $Q$ and $\Sigma$. We employ an IA model with an amplitude $A_{CFHTLenS}$ that is included in the parameter analysis. We find an improvement in the constraints on the MG parameters corresponding to $40-53\%$ increase on the figure of merit compared to previous studies, and GR is found consistent with the data at the $95\%$ CL. The bounds found on $A_{CFHTLenS}$ are sensitive to the MG parameterization used, and the correlations between $A_{CFHTLenS}$ and MG parameters are found to be weak to moderate. For all 3 MG parameterizations $A_{\rm CFHTLenS}$ is found consistent with zero when the whole lensing sample is used, however, when using the optimized early-type galaxy sample a significantly non-zero $A_{\rm CFHTLenS}$ is found for GR and the third MG parameterization. We find that the tensions observed in previous studies persist, and there is an indication that CMB data and lensing data prefer different values for MG parameters, particularly for the parameter $\Sigma$. The analysis of the confidence contours and probability distributions suggest that the bimodality found follows that of the known tension in the $\sigma_8$ parameter. (Abridged)

We present measurements of both scale- and time-dependent deviations from the standard gravitational field equations. These late-time modifications are introduced separately for relativistic and non-relativistic particles, by way of the parameters $G_{\rm matter}(k,z)$ and $G_{\rm light}(k,z)$ using two bins in both scale and time, with transition wavenumber $0.01$ Mpc$^{-1}$ and redshift 1. We emphasize the use of two dynamical probes to constrain this set of parameters, galaxy power spectrum multipoles and the direct peculiar velocity power spectrum, which probe fluctuations on different scales. The multipole measurements are derived from the WiggleZ and BOSS Data Release 11 CMASS galaxy redshift surveys and the velocity power spectrum is measured from the velocity sub-sample of the 6-degree Field Galaxy Survey. We combine with additional cosmological probes including baryon acoustic oscillations, Type Ia SNe, the cosmic microwave background (CMB), lensing of the CMB, and the temperature--galaxy cross-correlation. Using a Markov Chain Monte Carlo likelihood analysis, we find the inferred best-fit parameter values of $G_{\rm matter}(k,z)$ and $G_{\rm light}(k,z)$ to be consistent with the standard model at the $95\%$ confidence level. Furthermore, accounting for the Alcock-Paczynski effect, we perform joint fits for the expansion history and growth index gamma; we measure $\gamma = 0.665 \pm 0.0669$ ($68\%$ C.L) for a fixed expansion history, and $\gamma = 0.73^{+0.08}_{-0.10}$ ($68\%$ C.L) when the expansion history is allowed to deviate from $\Lambda$CDM. With a fixed expansion history the inferred value is consistent with GR at the $95\%$ C.L; alternatively, a $2\sigma$ tension is observed when the expansion history is not fixed, this tension is worsened by the combination of growth and SNe data.

In order to explain cosmic acceleration without invoking "dark" physics, we consider $f(R)$ modified gravity models, which replace the standard Einstein-Hilbert action in General Relativity with a higher derivative theory. We use data from the WiggleZ Dark Energy survey to probe the formation of structure on large scales which can place tight constraints on these models. We combine the large-scale structure data with measurements of the cosmic microwave background from the Planck surveyor. After parameterising the modification of the action using the Compton wavelength parameter $B_0$, we constrain this parameter using ISiTGR, assuming an initial non-informative log prior probability distribution of this cross-over scale. We find that the addition of the WiggleZ power spectrum provides the tightest constraints to date on $B_0$ by an order of magnitude, giving ${\rm log}_{10}(B_0) < -4.07$ at 95% confidence limit. Finally, we test whether the effect of adding the lensing amplitude $A_{\rm Lens}$ and the sum of the neutrino mass $\sum m_\nu$ is able to reconcile current tensions present in these parameters, but find $f(R)$ gravity an inadequate explanation.

Cosmological tests to distinguish between dark energy (DE) and modifications to gravity are a promising route to obtain clues on the origin of cosmic acceleration. We study here the robustness of these tests to the presence of DE density, velocity, and anisotropic stress perturbations. We find that the dispersion in the growth index parameter remains small enough to distinguish between extreme cases of DE models and some commonly used modified gravity models. The sign of the slope parameter for a redshift dependent growth index was found to be inconsistent as an additional test in extreme cases of DE models with perturbations. Next, we studied the effect of DE perturbations on the modified growth (MG) parameters that enter the perturbed Einstein equations. We find that while the dark energy perturbations affect the MG parameters, the deviations remain smaller than those due to modified gravity models. Additionally, the deviations due to DE perturbations with a non-zero effective sound speed occur at scale ranges that are completely different than those due to some modified gravity models such as the $f(R)$ models. In the case of modified gravity models with zero anisotropic stress at late times, the simultaneous determination of the effective dark energy equation of state and the MG parameters can provide the distinction between these models and DE. The growth index test was found to be the most robust to these perturbations. The scale dependence of the MG parameters in some cases of modified gravity constitute a clear-cut discriminant regardless of any DE perturbations. In summary, we find that the currently proposed cosmological tests to distinguish between DE and modified gravity are robust to DE perturbations even for extreme cases. This is certainly the case even for DE models with equations of state of DE that fall well outside of current cosmological constraints.

It is well-known that allowing for spatial curvature affects constraints on cosmological parameters such as the dark energy equation of state parameters. Here we study the effect of curvature on constraints on parameters used to test general relativity (GR) at cosmological scales, commonly known as modified growth (MG) parameters, as while current data taken in the context of the $\Lambda$CDM model points to a universe that is flat or very close to it, this constraint may not hold in modified theories of gravity. Using the latest cosmological data sets we find that MG parameters are correlated with the curvature parameter $\Omega_k$ and the constraints on the MG parameters are weakened compared to when $\Omega_k$ is not included in the parameter analysis. We next use various future simulated data sets, including cosmic microwave background, weak lensing, and Integrated Sachs-Wolfe-galaxy cross-correlations, where the fiducial model is spatially curved but we assume a flat model when fitting the MG parameters. We find the assumption of a spatially flat model on a spatially curved universe does indeed cause an artificial shift in the constraints on the MG parameters, in some cases even producing an apparent deviation from GR in the MG parameter space. For our simulated data, tension with GR begins to manifest itself for fiducial models with $\abs{\Omega_k} \geq 0.02$ and apparent deviations appear for $\abs{\Omega_k} \geq 0.05$. We find that for negatively curved models the apparent deviation is more significant. The manifestation of this apparent deviation from GR due to the assumption of spatial flatness above leads one to conclude that, when using future high-precision data to perform these tests, spatial curvature must be included in the parameter analysis along with the other core cosmological parameters and the MG parameters.

We use cosmological constraints from current data sets and a figure of merit (FoM) approach to probe any deviations from general relativity (GR) at cosmological scales. The FoM approach is used to study the constraining power of various combinations of data sets on modified gravity (MG) parameters. We use recently refined HST-COSMOS weak-lensing tomography data, ISW-galaxy cross correlations from 2MASS and SDSS LRG surveys, matter power spectrum from SDSS-DR7 (MPK), WMAP7 temperature and polarization spectra, BAO from 2DF and SDSS-DR7, and Union2 compilation of supernovae, in addition to other bounds from H_0 measurements and BBN. We use 3 parametrizations of MG parameters that enter the perturbed field equations. In order to allow for variations with redshift and scale, the first 2 parametrizations use recently suggested functional forms while the third is based on binning methods. Using the first parametrization, we find that CMB + ISW + WL provides the strongest constraints on MG parameters followed by CMB+WL or CMB+MPK+ISW. Using the second parametrization or binning methods, CMB+MPK+ISW consistently provides some of the strongest constraints. This shows that the constraints are parametrization dependent. We find that adding up current data sets does not improve consistently uncertainties on MG parameters due to tensions between best-fit MG parameters preferred by different data sets. Furthermore, some functional forms imposed by the parametrizations can lead to an exacerbation of these tensions. Next, unlike some studies that used the CFHTLS lensing data, we do not find any deviation from GR using the refined HST-COSMOS data, confirming previous claims in those studies that their result may have been due to some systematic effect. Finally, we find in all cases that the values corresponding to GR are within the 95% confidence level contours for all data set combinations. (abridged)

The testing of general relativity at cosmological scales has become a possible and timely endeavor that is not only motivated by the pressing question of cosmic acceleration but also by the proposals of some extensions to general relativity that would manifest themselves at large scales of distance. We analyze here correlations between modified gravity growth parameters and some core cosmological parameters using the latest cosmological data sets including the refined Cosmic Evolution Survey 3D weak lensing. We provide parametrized modified growth equations and their evolution. We implement known functional and binning approaches, and propose a new hybrid approach to evolve modified gravity parameters in redshift (time) and scale. The hybrid parametrization combines a binned redshift dependence and a smooth evolution in scale avoiding a jump in the matter power spectrum. The formalism developed to test the consistency of current and future data with general relativity is implemented in a package that we make publicly available and call ISiTGR (Integrated Software in Testing General Relativity), an integrated set of modified modules for the publicly available packages CosmoMC and CAMB, including a modified version of the integrated Sachs-Wolfe-galaxy cross correlation module of Ho et al and a new weak-lensing likelihood module for the refined HST-COSMOS weak lensing tomography data. We obtain parameter constraints and correlation coefficients finding that modified gravity parameters are significantly correlated with \sigma_8 and mildly correlated with \Omega_m, for all evolution methods. The degeneracies between \sigma_8 and modified gravity parameters are found to be substantial for the functional form and also for some specific bins in the hybrid and binned methods indicating that these degeneracies will need to be taken into consideration when using future high precision data.

We use current and future simulated data of the growth rate of large scale structure in combination with data from supernova, BAO, and CMB surface measurements, in order to put constraints on the growth index parameters. We use a recently proposed parameterization of the growth index that interpolates between a constant value at high redshifts and a form that accounts for redshift dependencies at small redshifts. We also suggest here another exponential parameterization with a similar behaviour. The redshift dependent parametrizations provide a sub-percent precision level to the numerical growth function, for the full redshift range. Using these redshift parameterizations or a constant growth index, we find that current available data from galaxy redshift distortions and Lyman-alpha forests is unable to put significant constraints on any of the growth parameters. For example both $\Lambda$CDM and flat DGP are allowed by current growth data. We use an MCMC analysis to study constraints from future growth data, and simulate pessimistic and moderate scenarios for the uncertainties. In both scenarios, the redshift parameterizations discussed are able to provide significant constraints and rule out models when incorrectly assumed in the analysis. The values taken by the constant part of the parameterizations as well as the redshift slopes are all found to significantly rule out an incorrect background. We also find that, for our pessimistic scenario, an assumed constant growth index over the full redshift range is unable to rule out incorrect models in all cases. This is due to the fact that the slope acts as a second discriminator at smaller redshifts and therefore provide a significant test to identify the underlying gravity theory.

The question of whether or not the cosmological constant affects the bending of light around a concentrated mass has been the subject of some recent papers. We present here a simple, specific and transparent example where $\Lambda$ bending clearly takes place, and where it is clearly neither a coordinate effect nor an aberration effect. We then show that in some recent works using perturbation theory the $\Lambda$ contribution was missed because of initial too-stringent smallness assumptions. Namely: Our method has been to insert a Kottler (Schwarzschild with $\Lambda$) vacuole into a Friedmann universe, and to calculate the total bending within the vacuole. We assume that no more bending occurs outside. It is important to observe that while the mass contribution to the bending takes place mainly quite near the lens, the $\Lambda$ bending continues throughout the vacuole. Thus if one deliberately restricts one's search for $\Lambda$ bending to the immediate neighborhood of the lens, one will not find it. Lastly, we show that the $\Lambda$ bending also follows from standard Weyl focusing, and so again, it cannot be a coordinate effect.

The growth rate of matter perturbations can be used to distinguish between different gravity theories and to distinguish between dark energy and modified gravity at cosmological scales as an explanation to the observed cosmic acceleration. We suggest here parameterizations of the growth index as functions of the redshift. The first one is given by $\gamma(a)=\tilde\gamma(a) \frac{1}{1+(a_{_{ttc}}/a)}+\gamma_{_{early}} \frac{1}{1+(a/a_{_{ttc}})}$ that interpolates between a low/intermediate redshift parameterization $\tilde\gamma(a)=\gamma_{_{late}}(a)= \gamma_0 + (1-a) \gamma_a$ and a high redshift $\gamma_{_{early}}$ constant value. For example, our interpolated form $\gamma(a)$ can be used when including the CMB to the rest of the data while the form $\gamma_{_{late}}(a)$ can be used otherwise. It is found that the parameterizations proposed achieve a fit that is better than 0.004% for the growth rate in a $\Lambda$CDM model, better than 0.014% for Quintessence-Cold-Dark-Matter (QCDM) models, and better than 0.04% for the flat Dvali-Gabadadze-Porrati (DGP) model (with $\Omega_m^0=0.27$) for the entire redshift range up to $z_{_{CMB}}$. We find that the growth index parameters $(\gamma_0,\gamma_a)$ take distinctive values for dark energy models and modified gravity models, e.g. $(0.5655,-0.02718)$ for the $\Lambda$CDM model and $(0.6418,0.06261)$ for the flat DGP model. This provides a means for future observational data to distinguish between the models.

We derive new limits on the value of the cosmological constant, $\Lambda$, based on the Einstein bending of light by systems where the lens is a distant galaxy or a cluster of galaxies. We use an amended lens equation in which the contribution of $\Lambda$ to the Einstein deflection angle is taken into account and use observations of Einstein radii around several lens systems. We use in our calculations a Schwarzschild-de Sitter vacuole exactly matched into a Friedmann-Robertson-Walker background and show that a $\Lambda$-contribution term appears in the deflection angle within the lens equation. We find that the contribution of the $\Lambda$-term to the bending angle is larger than the second-order term for many lens systems. Using these observations of bending angles, we derive new limits on the value of $\Lambda$. These limits constitute the best observational upper bound on $\Lambda$ after cosmological constraints and are only two orders of magnitude away from the value determined by those cosmological constraints.