To reach the full potential for the next generation of weak lensing surveys, it is necessary to mitigate the contamination of intrinsic alignments (IA) of galaxies in the observed cosmic shear signal. The self calibration (SC) of intrinsic alignments provides an independent method to measure the IA signal from the survey data and the photometric redshift information. It operates differently from the marginalization method based on IA modeling. In this work, we present the first application of SC to the KiDS450 data and the KV450 data, to split directly the intrinsic shape - galaxy density (Ig) correlation signal and the gravitational shear - galaxy density (Gg) correlation signal, using the information from photometric redshift (photo-z). We achieved a clear separation of the two signals and performed several validation tests. Our measured signals are found to be in general agreement with the KiDS450 cosmic shear best-fit cosmology, for both lensing and IA measurements. For KV450, we use partial (high-z) data, and our lensing measurements are also in good agreement with KV450 cosmic shear best-fit, however, our IA signal suggests a larger IA amplitude. We discussed the impact of photo-z quality on IA detection and several other potential systematic biases. Finally, we discuss the potential application of the information extracted for both the lensing signal and the IA signal in future surveys.

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}.

Weak gravitational lensing is one of the most promising cosmological probes to constrain dark matter, dark energy and the nature of gravity at cosmic scales. Intrinsic alignments (IA) of galaxies have been recognized as one of the most serious systematic effects facing gravitational lensing. Such alignments must be isolated and removed to obtain a pure lensing signal. Furthermore, the alignments are related to the processes of galaxy formation, so their extracted signal can help in understanding such formation processes and improving their theoretical modeling. We report in this letter the first detection of the gravitational shear--intrinsic shape (GI) correlation and the intrinsic shape--galaxy density (Ig) correlation in a photometric redshift survey using the self-calibration method. These direct measurements are made from the KiDS-450 photometric galaxy survey with a significance of 2.74$\sigma$ in the third bin for the Ig correlation, and 2.73$\sigma$ for the GI cross-correlation between the third and fourth bins. The self-calibration method uses the information available from photometric surveys without needing to specify an IA model and will play an important role in validating IA models and IA mitigation in future surveys such as LSST, Euclid and WFIRST.

Consistency between cosmological data sets is essential for ongoing and future cosmological analyses. We first investigate the questions of stability and applicability of some moment-based inconsistency measures to multiple data sets. We show that the recently introduced index of inconsistency (IOI) is numerically stable while it can be applied to multiple data sets. We use an illustrative construction of constraints as well as an example with real data sets (i.e. WMAP versus Planck) to show some limitations of the application of the Karhunen-Loeve decomposition to discordance measures. Second, we perform various consistency analyzes using IOI between multiple current data sets while \textit{working with the entire common parameter spaces}. We find current Large-Scale-Structure (LSS) data sets (Planck CMB lensing, DES lensing-clustering and SDSS RSD) all to be consistent with one another. This is found to be not the case for Planck temperature (TT) versus polarization (TE,EE) data, where moderate inconsistencies are present. Noteworthy, we find a strong inconsistency between joint LSS probes and Planck with IOI=5.27, and a moderate tension between DES and Planck with IOI=3.14. Next, using the IOI metric, we compare the Hubble constant from five independent probes. We confirm previous strong tensions between local measurement (SH0ES) and Planck as well as between H0LiCOW and Planck, but also find new strong tensions between SH0ES measurement and the joint LSS probes with IOI=6.73 (i.e. 3.7-$\sigma$ in 1D) as well as between joint LSS and combined probes SH0ES+H0LiCOW with IOI=8.59 (i.e. 4.1-$\sigma$ in 1D). Whether due to systematic effects in the data sets or problems with the underlying model, sources of these old and new tensions need to be identified and dealt with.

Some inconsistencies between cosmological observations continue to persist. The currently reported inconsistencies can be due to problems with the underlying model or systematic effects associated with some particular datasets. In an era of precision cosmology, it is important to develop proper tools to quantify the degree of these inconsistencies and to look for methods to identify their causes. This turns out to be a task that is not free of ambiguities, especially when the model is multi-dimensional, which is usually the case in cosmology. Measures currently-proposed in the literature disagree even in Gaussian cases. We discuss, with illustrative examples, some requirements that should be fulfilled in Gaussian cases by inconsistency measures and suggest a guiding definition of (in)consistency. As an example, we show that the recently-proposed index of inconsistency (IOI) meets those requirements and is in line with the (in)consistency definition proposed. Next, we examine the common practice to convert some measures, including IOI and other similar quantities, to a probability to exceed or a significance level that depends on the number of parameters. In the context of quantifying inconsistencies, we show that such a procedure can underestimate inconsistencies when there is more than one model parameter. We also discuss multiple-dataset comparisons and introduce a new tool based on the multi-dataset IOI that can identify outlying constraints when present. Comparison of constraints from various datasets can help identify the source of inconsistencies and will become important as more independent constraints will become available from ongoing and future surveys.

One of the main science goals of the Large Synoptic Survey Telescope (LSST) is to uncover the nature of cosmic acceleration. In the base analysis, possible deviations from the Lambda-Cold-Dark-Matter ($\Lambda$CDM) background evolution will be probed by fitting a $w(z)$CDM model, which allows for a redshift-dependent dark energy equation of state with $w(z)$, within general relativity (GR). A rich array of other phenomena can arise due to deviations from the standard $\Lambda$CDM+GR model though, including modifications to the growth rate of structure and lensing, and novel screening effects on non-linear scales. Concrete physical models are needed to provide consistent predictions for these (potentially small) effects, to give us the best chance of detecting them and separating them from astrophysical systematics. A complex plethora of possible models has been constructed over the past few decades, with none emerging as a particular favorite. This document prioritizes a subset of these models along with rationales for further study and inclusion into the LSST Dark Energy Science Collaboration (DESC) data analysis pipelines, based on their observational viability, theoretical plausibility, and level of theoretical development. We provide references and theoretical expressions to aid the integration of these models into DESC software and simulations, and give justifications for why other models were not prioritized. While DESC efforts are free to pursue other models, we provide here guidelines on which theories appear to have higher priority for collaboration efforts due to their perceived promise and greater instructional value.

The Core Cosmology Library (CCL) provides routines to compute basic cosmological observables to a high degree of accuracy, which have been verified with an extensive suite of validation tests. Predictions are provided for many cosmological quantities, including distances, angular power spectra, correlation functions, halo bias and the halo mass function through state-of-the-art modeling prescriptions available in the literature. Fiducial specifications for the expected galaxy distributions for the Large Synoptic Survey Telescope (LSST) are also included, together with the capability of computing redshift distributions for a user-defined photometric redshift model. A rigorous validation procedure, based on comparisons between CCL and independent software packages, allows us to establish a well-defined numerical accuracy for each predicted quantity. As a result, predictions for correlation functions of galaxy clustering, galaxy-galaxy lensing and cosmic shear are demonstrated to be within a fraction of the expected statistical uncertainty of the observables for the models and in the range of scales of interest to LSST. CCL is an open source software package written in C, with a python interface and publicly available at https://github.com/LSSTDESC/CCL.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

We review recent developments and results in testing general relativity (GR) at cosmological scales. The subject has witnessed rapid growth during the last two decades with the aim of addressing the question of cosmic acceleration and the dark energy associated with it. However, with the advent of precision cosmology, it has also become a well-motivated endeavor by itself to test gravitational physics at cosmic scales. We overview cosmological probes of gravity, formalisms and parameterizations for testing deviations from GR at cosmological scales, selected modified gravity (MG) theories, gravitational screening mechanisms, and computer codes developed for these tests. We then provide summaries of recent cosmological constraints on MG parameters and selected MG models. We supplement these cosmological constraints with a summary of implications from the recent binary neutron star merger event. Next, we summarize some results on MG parameter forecasts with and without astrophysical systematics that will dominate the uncertainties. The review aims at providing an overall picture of the subject and an entry point to students and researchers interested in joining the field. It can also serve as a quick reference to recent results and constraints on testing gravity at cosmological scales.

We introduce a self-calibration method that can be applied to the intrinsic ellipticity--intrinsic ellipticity (II) and gravitational shear -- intrinsic ellipticity (GI) types of intrinsic alignment of galaxies. The method combines previous self-calibration techniques with modifications to one of them in order to use auto-spectra in addition to cross-spectra between redshift bins. This allows one to use the self-calibration while preserving all the constraining power of cosmic shear from surveys. We show that the new method provides more flexibility in using various redshift bin widths. We perform cosmological parameter constraint forecast when this method is applied to the Large Synoptic Survey Telescope (LSST). Compared to the original self-calibration, we find that the new method provides further significant reduction in any residual shift in the cosmological parameters (e.g. factors of $2-4$ for the dark energy equation of state) which is promising for accurate cosmology.

Intrinsic alignments (IA) of galaxies have been recognized as one of the most serious contaminants to weak lensing. These systematics need to be isolated and mitigated in order for ongoing and future lensing surveys to reach their full potential. The IA self-calibration (SC) method was shown in previous studies to be able to reduce the GI contamination by up to a factor of 10 for the 2-point and 3-point correlations. The SC method does not require to assume an IA model in its working and can extract the GI signal from the same photo-z survey offering the possibility to test and understand structure formation scenarios and their relationship to IA models. In this paper, we study the effects of the IA SC mitigation method on the precision and accuracy of cosmological parameter constraints from future cosmic shear surveys LSST, WFIRST and Euclid. We perform analytical and numerical calculations to estimate the loss of precision and the residual bias in the best fit cosmological parameters after the self-calibration is performed. We take into account uncertainties from photometric redshifts and the galaxy bias. We find that the confidence contours are slightly inflated from applying the SC method itself while a significant increase is due to the inclusion of the photo-z uncertainties. The bias of cosmological parameters is reduced from several-$\sigma$, when IA is not corrected for, to below 1-$\sigma$ after SC is applied. These numbers are comparable to those resulting from applying the method of marginalizing over IA model parameters despite the fact that the two methods operate very differently. We conclude that implementing the SC for these future cosmic-shear surveys will not only allow one to efficiently mitigate the GI contaminant but also help to understand their modeling and link to structure formation.

We examine systematically the (in)consistency between cosmological constraints as obtained from various current data sets of the expansion history, Large Scale Structure (LSS), and Cosmic Microwave Background (CMB) from Planck. We run (dis)concordance tests within each set and across the sets using a recently introduced index of inconsistency (IOI) capable of dissecting inconsistencies between two or more data sets. First, we compare the constraints on $H_0$ from five different methods and find that the IOI drops from 2.85 to 0.88 (on Jeffreys' scales) when the local $H_0$ measurements is removed. This seems to indicate that the local measurement is an outlier, thus favoring a systematics-based explanation. We find a moderate inconsistency (IOI=2.61) between Planck temperature and polarization. We find that current LSS data sets including WiggleZ, SDSS RSD, CFHTLenS, CMB lensing and SZ cluster count, are consistent one with another and when all combined. However, we find a persistent moderate inconsistency between Planck and individual or combined LSS probes. For Planck TT+lowTEB versus individual LSS probes, the IOI spans the range 2.92--3.72 and increases to 3.44--4.20 when the polarization data is added in. The joint LSS versus the combined Planck temperature and polarization has an IOI of 2.83 in the most conservative case. But if Planck lowTEB is added to the joint LSS to constrain $\tau$ and break degeneracies, the inconsistency between Planck and joint LSS data increases to the high-end of the moderate range with IOI=4.81. Whether due to systematic effects in the data or to the underlying model, these inconsistencies need to be resolved. Finally, we perform forecast calculations using LSST and find that the discordance between Planck and future LSS data, if it persists as present, can rise up to a high IOI of 17, thus falling in the strong range of inconsistency. (Abridged).

The continuous progress toward more precise cosmological surveys and experiments has galvanized recent interest into consistency tests on cosmological parameters and models. At the heart of this effort is quantifying the degree of inconsistency between two or more cosmological data sets. We introduce an intuitive moment-based measure we call the \textit{index of inconsistency} (IOI) and show that it is sensitive to the separation of the means, the size of the constraint ellipsoids, and their orientations in the parameter space. We find that it tracks accurately the inconsistencies when present. Next, we show that parameter marginalization can cause a loss of information on the inconsistency between two experiments and we quantify such a loss using the drop in IOI. In order to zoom on a given parameter, we define the relative residual IOI and the relative drop in IOI. While these two quantities can provide insights on the parameters that are most responsible for inconsistencies, we find that the full IOI applied to the whole parameter spaces is what must be used to correctly reflect the degree of inconsistency between two experiments. We discuss various properties of IOI, provide its eigenmode decomposition, and compare it to other measures of discordance. Finally, we apply IOI to current geometry data sets (i.e. an improved Supernovae Type Ia compilation, baryon acoustic oscillations from 6dF, SDSS MGS and Lyman-$\alpha$ forest, and high-$\ell$ CMB temperature data from Planck-2015) versus growth data sets (i.e. Redshift Space Distortions from WiggleZ and SDSS, Weak Lensing from CFHTLenS, CMB Lensing, Sunyav-Zeldovich effect, and low-$\ell$ CMB temperature and polarization data from Planck-2015). We find that a persistent inconsistency is present between the two data sets. This could reflect the presence of systematics in the data or inconsistencies in the underlying model.

Primordial gravitational waves constitute a promising probe of the very early Universe and the laws of gravity. We study in this work changes to tensor-mode perturbations (TMPs) that can arise in various proposed modified gravity (MG) theories. These include additional friction effects, nonstandard dispersion relations involving a massive graviton, a modified speed, and a small-scale modification. We introduce a physically motivated parameterization of these effects and use current available data to obtain exclusion regions in the parameter spaces. Taking into account the foreground subtraction, we then perform a forecast analysis focusing on tensor-mode MG parameters as constrained by future experiments COrE, Stage-IV and PIXIE. For a fiducial value of the tensor-to-scalar ratio $r=0.01$, we find that an additional friction of $3.5\sim4.5\%$ compared to GR will be detected at $3$-$\sigma$ by these experiments, while a decrease in friction will be more difficult to detect. The speed of gravitational waves needs to be by $5\sim15\%$ different from the speed of light for detection. We find that the minimum detectable graviton mass is about $7.8\sim9.7\times 10^{-33}$ eV, which is of the same order of magnitude as the graviton mass that allows massive gravity theories to produce cosmic acceleration. Finally, we study TMPs in MG during inflation using our parameterization. We find that the tensor spectral index is related to $r$ and the friction parameter $\nu_0$ by $n_T=-3\nu_0-r/8$. Assuming that the friction parameter is unchanged throughout the history of the Universe and is much larger than $r$, future experiments will be able to distinguish this MG consistency relation from the standard inflation one, and thus can be used as a further test of MG. In summary, TMPs and cosmic-microwave-background B-mode polarization provide a complementary avenue to test gravity theories. (Abridged)

One of the outstanding problems in general relativistic cosmology is that of the averaging. That is, how the lumpy universe that we observe at small scales averages out to a smooth Friedmann-Lemaitre-Robertson-Walker (FLRW) model. The root of the problem is that averaging does not commute with the Einstein equations that govern the dynamics of the model. This leads to the well-know question of backreaction in cosmology. In this work, we approach the problem using the covariant framework of Macroscopic Gravity (MG). We use its cosmological solution with a flat FLRW macroscopic background where the result of averaging cosmic inhomogeneities has been encapsulated into a backreaction density parameter denoted $\Omega_\mathcal{A}$. We constrain this averaged universe using available cosmological data sets of expansion and growth including, for the first time, a full CMB analysis from Planck temperature anisotropy and polarization data, the supernovae data from Union 2.1, the galaxy power spectrum from WiggleZ, the weak lensing tomography shear-shear cross correlations from the CFHTLenS survey and the baryonic acoustic oscillation data from 6Df, SDSS DR7 and BOSS DR9. We find that $-0.0155 \le \Omega_\mathcal{A} \le 0$ (at the 68\% CL) thus providing a tight upper-bound on the backreaction term. We also find that the term is strongly correlated with cosmological parameters such $\Omega_\Lambda$, $\sigma_8$ and $H_0$. While small, a backreaction density parameter of a few percent should be kept in consideration along with other systematics for precision cosmology.

The scenario consistent with a wealth of observations for the missing mass problem is that of weakly interacting dark matter particles. However, arguments or proposals for a Newtonian or relativistic modified gravity scenario continue to be made. A distinguishing characteristic between the two scenarios is that dark matter particles can produce a gravitational effect, in principle, without the need of baryons while this is not the case for the modified gravity scenario where such an effect must be correlated with the amount of baryonic matter. We consider here ultra-faint dwarf (UFD) galaxies as a promising arena to test the two scenarios based on the above assertion. We compare the correlation of the luminosity with the velocity dispersion between samples of UFD and non-UFD galaxies, finding a significant loss of correlation for UFD galaxies. For example, we find for 28 non-UFD galaxies a strong correlation coefficient of -0.688 which drops to -0.077 for the 23 UFD galaxies. Incoming and future data will determine whether the observed stochasticity for UFD galaxies is physical or due to systematics in the data. Such a loss of correlation (if it is to persist) is possible and consistent with the dark matter scenario for UFD galaxies but would constitute a new challenge for the modified gravity scenario.

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 investigate the effect of averaging inhomogeneities on expansion and large-scale structure growth observables using the exact and covariant framework of Macroscopic Gravity (MG). It is well-known that applying the Einstein's equations and spatial averaging do not commute and lead to the averaging problem. For the MG formalism applied to the Friedmann-Lemaitre-Robertson-Walker (FLRW) metric, this gives an extra dynamical term encapsulated as an averaging density parameter denoted $\Omega_A$. An exact isotropic cosmological solution of MG for the flat FLRW metric is already known in the literature, we derive here an anisotropic exact solution. Using the isotropic solution, we compare the expansion history to current data of distances to supernovae, Baryon Acoustic Oscillations, CMB last scattering surface, and Hubble constant measurements, and find $-0.05 \le \Omega_A \le 0.07$ (at the 95% CL). For the flat metric case this reduces to $-0.03 \le \Omega_A \le 0.05$. We also find that the inclusion of this term in the fits can shift the values of the usual cosmological parameters by a few to several percents. Next, we derive an equation for the growth rate of large scale structure in MG that includes a term due to the averaging and compare it to that of the LCDM concordance model. We find that an $\Omega_A$ of an amplitude range within the bounds above leads to a relative deviation of the growth from that of the LCDM of up to 2-4% at late times. Thus, the shift in the growth could be of comparable amplitude to that caused by similar changes in cosmological parameters like the dark energy density parameter or its equation of state. The effect could also be comparable in amplitude to some systematic effects considered for future surveys. This indicates that the averaging term and its possible effect need to be tightly constrained in future precision cosmological studies. (Abridged)

We study effects of inhomogeneities on distance measures in an exact relativistic Swiss-cheese model of the universe, focusing on the distance modulus. The model has LCDM background dynamics, and the `holes' are non-symmetric structures described by the Szekeres metric. The Szekeres exact solution of Einstein's equations, which is inhomogeneous and anisotropic, allows us to capture potentially relevant effects on light propagation due to nontrivial evolution of structures in an exact framework. Light beams traversing a single Szekeres structure in different ways can experience either magnification or demagnification, depending on the particular path, and we explore a small additional effect due time evolution of the structures. We study the probability distributions of $\Delta\mu=\mu_{LCDM}-\mu_{SC}$ for sources at different redshifts in various Swiss-cheese constructions, where the light beams travel through a large number of randomly oriented holes with random impact parameters. We find for $\Delta\mu$ the dispersions $0.004\le \sigma_{\Delta\mu} \le 0.008$ mag for sources with redshifts $1.0\le z \le 1.5$, which are smaller than the intrinsic dispersion of, for example, magnitudes of type Ia supernovae. The shapes of the distributions we obtain for our Swiss-cheese constructions are peculiar in the sense that they are not consistently skewed toward the demagnification side, as they are in analyses of lensing in cosmological simulations, reflecting a limitation of these constructions. This could be the result of requiring the continuity of Einstein's equations throughout the overall spacetime patchwork, which imposes the condition that compensating overdense shells must accompany the underdense void regions in the holes. The possibility to explore other uses of these constructions that could circumvent this limitation and lead to different statistics remains open. (Abridged)

The wealth of incoming and future cosmological observations will allow us to map out the structure and evolution of the observable universe to an unprecedented level of precision. Among these observations is the weak gravitational lensing of galaxies, e.g., cosmic shear that measures the minute distortions of background galaxy images by intervening cosmic structure. Weak lensing and cosmic shear promise to be a powerful probe of astrophysics and cosmology, constraining models of dark energy, measuring the evolution of structure in the universe, and testing theories of gravity on cosmic scales. However, the intrinsic alignment of galaxies -- their shape and orientation before being lensed -- may pose a great challenge to the use of weak gravitational lensing as an accurate cosmological probe, and has been identified as one of the primary physical systematic biases in cosmic shear studies. Correlations between this intrinsic alignment and the lensing signal can persist even for large physical separations, and isolating the effect of intrinsic alignment from weak lensing is not trivial. A great deal of work in the last two decades has been devoted to understanding and characterizing this intrinsic alignment, which is also a direct and complementary probe of structure formation and evolution in its own right. In this review, we report in a systematic way the state of our understanding of the intrinsic alignment of galaxies, with a particular emphasis on its large-scale impact on weak lensing measurements and methods for its isolation or mitigation. (Abridged)

The study of relativistic, higher order and nonlinear effects has become necessary in recent years in the pursuit of precision cosmology. We develop and apply here a framework to study gravitational lensing in exact models in general relativity that are not restricted to homogeneity and isotropy, and where full nonlinearity and relativistic effects are included. We apply the framework to a specific, anisotropic galaxy cluster model which is based on a modified NFW halo density profile and described by the Szekeres metric. We examine the effects of increasing levels of anisotropy in the galaxy cluster on lensing observables like the convergence and shear for various lensing geometries, finding a strong nonlinear response in both the convergence and shear for rays passing through anisotropic regions of the cluster. Deviation from the expected values in a spherically symmetric structure are asymmetric with respect to path direction and thus will persist as a statistical effect when averaged over some ensemble of such clusters. The resulting relative difference in various geometries can be as large as approximately 2%, 8%, and 24% in the measure of convergence for levels of anisotropy of 5%, 10%, and 15%, respectively, as a fraction of total cluster mass. For the total magnitude of shear, the relative difference can grow near the center of the structure to be as large as 15%, 32%, and 44% for the same levels of anisotropy, averaged over the two extreme geometries. The convergence is impacted most strongly for rays which pass in directions along the axis of maximum dipole anisotropy in the structure, while the shear is most strongly impacted for rays which pass in directions orthogonal to this axis, as expected. These effects due to anisotropic structures will affect lensing measurements and must be fully examined in an era of precision cosmology.

We introduce here a cross-correlation term between CMB lensing and galaxy intrinsic alignment, noted hereby as $\phi I$. This effect acts as a contaminant to the cross-correlation between CMB lensing and galaxy lensing. The latter cross-correlation has recently been detected for the first time, and measurements will greatly improve as the area of overlap between galaxy and CMB surveys increases and measurements of the CMB polarization become more significant. This will constitute a powerful probe for studying the structure and evolution of the universe. The magnitude of the $\phi I$ term is found to be about 15% of the pure CMB lensing-galaxy lensing component and acts to reduce the magnitude of its measured spectrum. This offset in the spectrum will strongly impact its use for precision cosmological study if left unmitigated. We also propose here a method to calibrate this $\phi I$ contamination through use of a scaling relation that allows one to reduce the impact of $\phi I$ by a factor of 20 or more in all redshift bins, which would reduce its magnitude down to detection limits in almost all cases. This will allow the full use of this probe for precision cosmology.

Probes of cosmic expansion constitute the main basis for arguments to support or refute a possible apparent acceleration due to different expansion rates in the universe as described by inhomogeneous cosmological models. We present in this Letter a separate argument based on results from an analysis of the growth rate of large-scale structure in the universe as modeled by the inhomogeneous cosmological models of Szekeres. We use the models with no assumptions of spherical or axial symmetries. We find that while the Szekeres models can fit very well the observed expansion history without a $\Lambda$, they fail to produce the observed late-time suppression in the growth unless $\Lambda$ is added to the dynamics. A simultaneous fit to the supernova and growth factor data shows that the cold dark matter model with a cosmological constant ($\Lambda$CDM) provides consistency with the data at a confidence level of 99.65% while the Szekeres model without $\Lambda$ achieves only a 60.46% level. When the data sets are considered separately, the Szekeres with no $\Lambda$ fits the supernova data as well as the $\Lambda$CDM does, but provides a very poor fit to the growth data with only 31.31% consistency level compared to 99.99% for the $\Lambda$CDM. This absence of late-time growth suppression in inhomogeneous models without a $\Lambda$ is consolidated by a physical explanation.

We study the effects and implications of anisotropies at the scale of galaxy clusters by building an exact general relativistic model of a cluster using the inhomogeneous and anisotropic Szekeres metric. The model is built from a modified Navarro-Frenk-White (NFW) density profile. We compare this to a corresponding spherically symmetric structure in the Lemaitre-Tolman (LT) model and quantify the impact of introducing varying levels of anisotropy. We examine two physical measures of gravitational infall -- the growth rate of density and the velocity of the source dust in the model. We introduce a generalization of the LT dust velocity profile for the Szekeres metric and demonstrate its consistency with the growth rate of density. We find that the growth rate of density in one substructure increases by 0.5%, 1.5%, and 3.75% for 5%, 10%, and 15% levels of introduced anisotropy, which is measured as the fractional displaced mass relative to the spherically symmetric case. The infall velocity of the dust is found to increase by 2.5, 10, and 20 km/s (0.5%, 2%, and 4.5%), respectively, for the same three levels of anisotropy. This response to the anisotropy in a structure is found to be strongly nonlinear with respect to the strength of anisotropy. These relative velocities correspond to an equivalent increase in the total mass of the spherically symmetric structure of 1%, 3.8%, and 8.4%, indicating that not accounting for the presence of anisotropic mass distributions in cluster models can strongly bias the determination of physical properties like the total mass.

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 the Szekeres inhomogeneous cosmological models to study the growth of large-scale structure in the universe including nonzero spatial curvature and a cosmological constant. In particular, we use the Goode and Wainwright formulation, as in this form the models can be considered to represent exact nonlinear perturbations of an averaged background. We identify a density contrast in both classes I and II of the models, for which we derive growth evolution equations. By including Lambda, the time evolution of the density contrast as well as kinematic quantities can be tracked through the matter- and Lambda-dominated cosmic eras up to the present and into the future. In various models of class I and class II, the growth rate is found to be stronger than that of the LCDM cosmology, and it is suppressed at later times due to the presence of Lambda. We find that there are Szekeres models able to provide a growth history similar to that of LCDM while requiring less matter content and nonzero curvature, which speaks to the importance of including the effects of large-scale inhomogeneities in analyzing the growth of large-scale structure. Using data for the growth factor f from redshift space distortions and the Lyman-alpha forest, we obtain best fit parameters for class II models and compare their ability to match observations with LCDM. We find that there is negligible difference between best fit Szekeres models with no priors and those for LCDM, both including and excluding Lyman-alpha data. We also find that the growth index gamma parametrization cannot be applied in a simple way to the growth in Szekeres models, so a direct comparison of the function f to the data is performed. We conclude that the Szekeres models can provide an exact framework for the analysis of large-scale growth data that includes inhomogeneities and allows for different interpretations of observations. (abridged)

We extend the 3-point intrinsic alignment self-calibration technique to the gravitational shear-intrinsic ellipticity-intrinsic ellipticity (GII) bispectrum. The proposed technique will allow the measurement and removal of the GII intrinsic alignment contamination from the cross-correlation weak lensing signal. While significantly decreased from using cross-correlations instead of auto-correlation in a single photo-z bin, the GII contamination persists in adjacent photo-z bins and must be accounted for and removed from the lensing signal. We relate the GII and galaxy density-intrinsic ellipticity-intrinsic ellipticity (gII) bispectra through use of the galaxy bias, and develop the estimator necessary to isolate the gII bispectrum from observations. We find that the GII self-calibration technique performs at a level comparable to that of the gravitational shear-gravitational shear-intrinsic ellipticity correlation (GGI) self-calibration technique, with measurement error introduced through the gII estimator generally negligible when compared to minimum survey error. The accuracy of the relationship between the GII and gII bispectra typically allows the GII self-calibration to reduce the GII contamination by a factor of 10 or more for all adjacent photo-z bin combinations at $\ell>300$. For larger scales, we find that the GII contamination can be reduced by a factor of 3-5 or more. The GII self-calibration technique is complementary to the existing GGI self-calibration technique, which together will allow the total intrinsic alignment cross-correlation signal in 3-point weak lensing to be measured and removed.

The weak lensing signal (cosmic shear) has been shown to be strongly contaminated by the various types of galaxy intrinsic alignment (IA) correlations, which poses a barrier to precision weak lensing measurements. The redshift dependence of the IA signal has been used at the 2-point level to reduce this contamination by only measuring cross-correlations between large redshift bins, which significantly reduces the galaxy intrinsic ellipticity - intrinsic ellipticity (II) correlation. A self-calibration technique based on the redshift dependences of the IA correlations has also been proposed as a means to remove the 2-point IA contamination from the lensing signal. We explore here the redshift dependences of the IA and lensing bispectra in order to propose a self-calibration of the IA auto-correlations at the 3-point level (i.e. GGI, GII, and III), which can be well understood without the assumption of any particular IA model. We find that future weak lensing surveys will be able to measure the distinctive IA redshift dependence over ranges of $|\Delta z^P|\le 0.2$. Using conservative estimates of photo-z accuracy, we describe the 3-point self-calibration technique for the total IA signal, which can be accomplished through lensing tomography of photo-z bin size $\sim 0.01$. We find that the 3-point self-calibration can function at the accuracy of the 2-point technique with modest constraints in redshift separation. This allows the 3-point IA auto-correlation self-calibration technique proposed here to significantly reduce the contamination of the IA contamination to the weak lensing bispectrum.

This study belongs to a series devoted to using Szekeres inhomogeneous models to develop a theoretical framework where observations can be investigated with a wider range of possible interpretations. We look here into the growth of large-scale structure in the models. The Szekeres models are exact solutions to Einstein's equations that were originally derived with no symmetries. We use a formulation of the models that is due to Goode and Wainwright, who considered the models as exact perturbations of an FLRW background. Using the Raychaudhuri equation, we write for the two classes of the models, exact growth equations in terms of the under/overdensity and measurable cosmological parameters. The new equations in the overdensity split into two informative parts. The first part, while exact, is identical to the growth equation in the usual linearly perturbed FLRW models, while the second part constitutes exact non-linear perturbations. We integrate numerically the full exact growth rate equations for the flat and curved cases. We find that for the matter-dominated era, the Szekeres growth rate is up to a factor of three to five stronger than the usual linearly perturbed FLRW cases, reflecting the effect of exact Szekeres non-linear perturbations. The growth is also stronger than that of the non-linear spherical collapse model, and the difference between the two increases with time. This highlights the distinction when we use general inhomogeneous models where shear and a tidal gravitational field are present and contribute to the gravitational clustering. Additionally, it is worth observing that the enhancement of the growth found in the Szekeres models during the matter-dominated era could suggest a substitute to the argument that dark matter is needed when using FLRW models to explain the enhanced growth and resulting large-scale structures that we observe today (abridged)

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 intrinsic alignment (IA) of galaxies has been shown to be a significant barrier to precision cosmic shear measurements. (Zhang, 2010, ApJ, 720, 1090) proposed a self-calibration technique for the power spectrum to calculate the induced gravitational shear-galaxy intrinsic ellipticity correlation (GI) in weak lensing surveys with photo-z measurements which is expected to reduce the IA contamination by at least a factor of 10 for currently proposed surveys. We confirm this using an independent analysis and propose an expansion to the self-calibration technique for the bispectrum in order to calculate the dominant IA gravitational shear-gravitational shear-intrinsic ellipticity correlation (GGI) contamination. We first establish an estimator to extract the galaxy density-density-intrinsic ellipticity (ggI) correlation from the galaxy ellipticity-density-density measurement for a photo-z galaxy sample. We then develop a relation between the GGI and ggI bispectra, which allows for the estimation and removal of the GGI correlation from the cosmic shear signal. We explore the performance of these two methods, compare to other possible sources of error, and show that the GGI self-calibration technique can potentially reduce the IA contamination by up to a factor of 5-10 for all but a few bin choices, thus reducing the contamination to the percent level. The self-calibration is less accurate for adjacent bins, but still allows for a factor of three reduction in the IA contamination. The self-calibration thus promises to be an efficient technique to isolate both the 2-point and 3-point intrinsic alignment signals from weak lensing measurements.

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.

The Szekeres inhomogeneous models can be used to model the true lumpy universe that we observe. This family of exact solutions to Einstein's equations was originally derived with a general metric that has no symmetries. In this work, we develop and use a framework to integrate the angular diameter and luminosity distances in the general Szekeres models. We use the affine null geodesic equations in order to derive a set of first-order ordinary differential equations that can be integrated numerically to calculate the partial derivatives of the null vector components. These equations allow the integration in all generality of the distances in the Szekeres models and some examples are given. The redshift is determined from simultaneous integration of the null geodesic equations. This work does not assume spherical or axial symmetry, and the results will be useful for comparisons of the general Szekeres inhomogeneous models to current and future cosmological data.

This review surveys some recent developments concerning the effect of the cosmological constant on the bending of light by a spherical mass in Kottler (Schwarzchild-de Sitter) spacetime. Some proposals of how such an effect may be put into a setting of gravitational lensing in cosmology are also discussed. The picture that emerges from this review is that it seems fair to assert that the contribution of $\Lambda$ to the bending of light has by now been well established, while putting the $\Lambda$ light-bending terms into a cosmological context is still subject to some interpretation and requires further work and clarification.

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.

We consider recently proposed higher order gravity models where the action is built from the Einstein-Hilbert action plus a function f(G) of the Gauss-Bonnet invariant. The models were previously shown to pass physical acceptability conditions as well as solar system tests. In this paper, we compare the models to combined data sets of supernovae, baryon acoustic oscillations, and constraints from the CMB surface of last scattering. We find that the models provide fits to the data that are close to those of the LCDM concordance model. The results provide a pool of higher order gravity models that pass these tests and need to be compared to constraints from large scale structure and full CMB analysis.

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.

We compare higher order gravity models to observational constraints from magnitude-redshift supernova data, distance to the last scattering surface of the CMB, and Baryon Acoustic Oscillations. We follow a recently proposed systematic approach to higher order gravity models based on minimal sets of curvature invariants, and select models that pass some physical acceptability conditions (free of ghost instabilities, real and positive propagation speeds, and free of separatrices). Models that satisfy these physical and observational constraints are found in this analysis and do provide fits to the data that are very close to those of the LCDM concordance model. However, we find that the limitation of the models considered here comes from the presence of superluminal mode propagations for the constrained parameter space of the models.

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.

The growth rate of matter perturbation and the expansion rate of the Universe can be used to distinguish modified gravity and dark energy models in explaining cosmic acceleration. We explore here the inclusion of spatial curvature into the growth factor. We expand previous results using the approximation $\Omega_{m}^\gamma$ and then suggest a new form, $f_a=\Omega_m^\gamma+(\gamma-4/7)\Omega_k$, as an approximation for the growth factor when the curvature $\Omega_k$ is not negligible, and where the growth index $\gamma$ is usually model dependent. The expression recovers the standard results for the curved and flat $\Lambda$CDM and Dvali-Gabadadze-Porrati models. Using the best fit values of $\Omega_{m0}$ and $\Omega_{k0}$ to the expansion/distance measurements from Type Ia supernovae, baryon acoustic oscillation, WMAP5, and $H(z)$ data, we fit the growth index parameter to current growth factor data and obtain $\gamma_{\Lambda}(\Omega_{k} \not= 0) = 0.65^{+0.17}_{-0.15}$ and $\gamma_{DGP}(\Omega_{k} \not= 0) = 0.53^{+0.14}_{-0.12}$. For the $\Lambda$CDM model, the 1-$\sigma$ observational bounds are found consistent with theoretical value, unlike the case for the Dvali-Gabadadze-Porrati model. We also find that the current data we used is not enough to put significant constraints when the 3 parameters in $f_a$ are fit simultaneously. Importantly, we find that, in the presence of curvature, the analytical expression proposed for $f_a$ provides a better fit to the growth factor than other forms and should be useful for future high precision missions and studies.

We use the Szekeres inhomogeneous relativistic models in order to fit supernova combined data sets. We show that with a choice of the spatial curvature function that is guided by current observations, the models fit the supernova data almost as well as the LCDM model without requiring a dark energy component. The Szekeres models were originally derived as an exact solution to Einstein's equations with a general metric that has no symmetries and are regarded as good candidates to model the true lumpy universe that we observe. The null geodesics in these models are not radial. The best fit model found is also consistent with the requirement of spatial flatness at CMB scales. The first results presented here seem to encourage further investigations of apparent acceleration using various inhomogeneous models and other constraints from CMB and large structure need to be explored next.

Higher-order gravity models have been recently the subject of much attention in the context of cosmic acceleration. These models are derived by adding various curvature invariants to the Einstein-Hilbert action. Several studies showed that these models can have late-time self-acceleration and could, in some cases, fit various observational constraints. In view of the infinite spectrum of invariants that could be built from curvature tensors, we propose here a method based on minimal sets of independent invariants as a possible route for a systematic approach to these models. We explore a connection made between theorems on bases in invariants theory in relativity and higher-order cosmological models. A dynamical system analysis is performed and some models with accelerating attractors are discussed. The asymptotic behavior of the models is also studied using analytical calculations.

The contribution of the cosmological constant to the deflection angle and the time delays are derived from the integration of the gravitational potential as well as from Fermat's Principle. The findings are in agreement with recent results using exact solutions to Einstein's equations and reproduce precisely the new $\Lambda$-term in the bending angle and the lens equation. The consequences on time delay expressions are explored. While it is known that $\Lambda$ contributes to the gravitational time delay, it is shown here that a new $\Lambda$-term appears in the geometrical time delay as well. Although these newly derived terms are perhaps small for current observations, they do not cancel out as previously claimed. Moreover, as shown before, at galaxy cluster scale, the $\Lambda$ contribution can be larger than the second-order term in the Einstein deflection angle for several cluster lens systems.

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.

We expand previous work on an inverse approach to Einstein Field Equations where we include fluids with energy flux and consider the vanishing of the anisotropic stress tensor. We consider the approach using warped product spacetimes of class $B_1$. Although restricted, these spacetimes include many exact solutions of interest to compact object studies and to cosmological models studies. The question explored here is as follows: given a spacetime metric, what fluid flow (timelike congruence), if any, could generate the spacetime via Einstein's equations. We calculate the flow from the condition of a vanishing anisotropic stress tensor and give results in terms of the metric functions in the three canonical types of coordinates. A condition for perfect fluid sources is also provided. The framework developed is algorithmic and suited for the study and validation of exact solutions using computer algebra systems. The framework can be applied to solutions in comoving and non-comoving frames of reference, and examples in different types of coordinates are worked out.

Correlations between intrinsic shear and the density field on large scales, a potentially important contaminant for cosmic shear surveys, have been robustly detected at low redshifts with bright galaxies in SDSS data. Here we present a more detailed characterization of this effect, which can cause anti-correlations between gravitational lensing shear and intrinsic ellipticity (GI correlations). This measurement uses 36278 Luminous Red Galaxies (LRGs) from the SDSS spectroscopic sample with 0.15<z<0.35, split by redshift and luminosity; 7758 LRGs from the 2SLAQ Survey at 0.4<z<0.8; and a variety of other SDSS samples from previous, related work. We find >3sigma detections of the effect for all galaxy subsamples within the SDSS LRG sample; for the 2SLAQ sample, we find a 2sigma detection for a bright subsample, and no detection for a fainter subsample. Fitting formulae are provided for the scaling of the GI correlations with luminosity, transverse separation, and redshift. We estimate contamination in the measurement of sigma_8 for future cosmic shear surveys on the basis of the fitted dependence of GI correlations on galaxy properties. We find contamination to the power spectrum ranging from -1.5 (optimistic) to -33 per cent (pessimistic) for a toy cosmic shear survey using all galaxies to a depth of R=24 using scales l~500. This corresponds to a bias in sigma_8 of Delta sigma_8=-0.004 (optimistic), -0.02 (central), or -0.10 (pessimistic). We provide a prescription for inclusion of this error in cosmological parameter estimation codes. The principal uncertainty is in the treatment of the L<=L* blue galaxies. Characterization of the tidal alignments of these galaxies, especially at redshifts relevant for cosmic shear, should be a high priority for the cosmic shear community. (Abridged)

Associated with the cosmic acceleration are the old and new cosmological constant problems, recently put into the more general context of the dark energy problem. In broad terms, the old problem is related to an unexpected order of magnitude of this component while the new problem is related to this magnitude being of the same order of the matter energy density during the present epoch of cosmic evolution. Current plans to measure the equation of state certainly constitute an important approach; however, as we discuss, this approach is faced with serious feasibility challenges and is limited in the type of conclusive answers it could provide. Therefore, is it really too early to seek actively for new tests and approaches to these problems? In view of the difficulty of this endeavor, we argue here that a good place to start is by questioning some of the assumptions underlying the formulation of these problems and finding new ways to put this questioning to the test. First, we calculate how much fine tuning the cosmic coincidence problem represents. Next, we discuss the potential of some cosmological probes such as gravitational lensing to identify novel tests to probe dark energy questions and assumptions. Then, motivated by some theorems in General Relativity, we discuss if the full identification of the cosmological constant with vacuum energy is unquestionable. Also, we point out the relevance of experiments at the interface of astrophysics and quantum field theory, such as the Casimir effect in gravitational and cosmological contexts. We conclude that challenging some of the assumptions underlying the cosmological constant problems and putting them to the test may prove useful and necessary to make progress on these questions. (Abridged)

We study the effect of the cosmological constant $\Lambda$ on the bending of light by a concentrated spherically symmetric mass. Contrarily to previous claims, we show that when the Schwarzschild-de Sitter geometry is taken into account, $\Lambda$ does indeed contribute to the bending.

If general relativity is the correct theory of physics on large scales, then there is a differential equation that relates the Hubble expansion function, inferred from measurements of angular diameter distance and luminosity distance, to the growth rate of large scale structure. For a dark energy fluid without couplings or an unusual sound speed, deviations from this consistency relationship could be the signature of modified gravity on cosmological scales. We propose a procedure based on this consistency relation in order to distinguish between some dark energy models and modified gravity models. The procedure uses different combinations of cosmological observations and is able to find inconsistencies when present. As an example, we apply the procedure to a universe described by a recently proposed 5-dimensional modified gravity model. We show that this leads to an inconsistency within the dark energy parameter space detectable by future experiments.

The power spectrum of weak lensing shear caused by large-scale structure is an emerging tool for precision cosmology, in particular for measuring the effects of dark energy on the growth of structure at low redshift. One potential source of systematic error is intrinsic alignments of ellipticities of neighbouring galaxies (II correlation) that could mimic the correlations due to lensing. A related possibility pointed out by Hirata and Seljak (2004) is correlation between the intrinsic ellipticities of galaxies and the density field responsible for gravitational lensing shear (GI correlation). We present constraints on both the II and GI correlations using 265 908 spectroscopic galaxies from the SDSS, and using galaxies as tracers of the mass in the case of the GI analysis. The availability of redshifts in the SDSS allows us to select galaxies at small radial separations, which both reduces noise in the intrinsic alignment measurement and suppresses galaxy- galaxy lensing (which otherwise swamps the GI correlation). While we find no detection of the II correlation, our results are nonetheless statistically consistent with recent detections found using the SuperCOSMOS survey. In contrast, we have a clear detection of GI correlation in galaxies brighter than L* that persists to the largest scales probed (60 Mpc/h) and with a sign predicted by theoretical models. This correlation could cause the existing lensing surveys at z~1 to underestimate the linear amplitude of fluctuations by as much as 20% depending on the source sample used, while for surveys at z~0.5 the underestimation may reach 30%. (Abridged.)

We study future constraints on dark energy parameters determined from several combinations of CMB experiments, supernova data, and weak lensing surveys with and without tomography. In this analysis, we look in particular for combinations that will bring the uncertainties to a level of precision tight enough (a few percent) to answer decisively some of the dark energy questions. We probe the dark energy using two variants of its equation of state, and its energy density.We consider a set of 13 cosmological and systematic parameters, and assume reasonable priors on the lensing and supernova systematics. We consider various lensing surveys: a wide survey with f_{sky}=0.7, and with 2 (WLT2) and 5 (WLT5) tomographic bins; a deep survey with 10 bins (WLT10). The constraints found from Planck, 2000 supernovae with z_max=0.8, and WLT2 are: {sigma(w_0)=0.086, sigma(w_1)=0.069}, {sigma(w_0)=0.088, sigma(w_a)=0.11}, and {sigma(E_1)=0.029, sigma(E_2)=0.065}. With 5 bins, we find {sigma(w_0)=0.04, sigma(w_1)=0.034}, {sigma(w_0)=0.041, sigma(w_a)=0.056}, and {sigma(E_1)=0.012, sigma(E_2)=0.049}. Finally, we find from Planck, 2000 supernovae with z_max=1.5, and WLT10 with f_{sky}=0.1: {sigma(w_0)=0.032, sigma(w_1)=0.027}, {sigma(w_0)=0.033, sigma(w_a)=0.040}, and {sigma(E_1)=0.01, sigma(E_2)=0.04}. Although some worries remain about other systematics, our study shows that after the combination of the 3 probes, lensing tomography with many redshift bins and large coverages of the sky has the potential to add key improvements to the dark energy parameter constraints. However, the requirement for very ambitious and sophisticated surveys in order to achieve some of the constraints or to improve them suggests the need for new tests to probe the nature of dark energy in addition to constraining its equation of state. (Abriged)

We consider how well the dark energy equation of state $w$ as a function of red shift $z$ will be measured using current and anticipated experiments. We use a procedure which takes fair account of the uncertainties in the functional dependence of $w$ on $z$, as well as the parameter degeneracies, and avoids the use of strong prior constraints. We apply the procedure to current data from WMAP, SDSS, and the supernova searches, and obtain results that are consistent with other analyses using different combinations of data sets. The effects of systematic experimental errors and variations in the analysis technique are discussed. Next, we use the same procedure to forecast the dark energy constraints achieveable by the end of the decade, assuming 8 years of WMAP data and realistic projections for ground-based measurements of supernovae and weak lensing. We find the $2 \sigma$ constraints on the current value of $w$ to be $\Delta w_0 (2 \sigma) = 0.20$, and on $dw/dz$ (between $z=0$ and $z=1$) to be $\Delta w_1 (2 \sigma)=0.37$. Finally, we compare these limits to other projections in the literature. Most show only a modest improvement; others show a more substantial improvement, but there are serious concerns about systematics. The remaining uncertainty still allows a significant span of competing dark energy models. Most likely, new kinds of measurements, or experiments more sophisticated than those currently planned, are needed to reveal the true nature of dark energy.

Weak lensing is a potentially robust and model-independent cosmological probe, but its accuracy is dependent on knowledge of the redshift distribution of the source galaxies used. The most robust way to determine the redshift distribution is via spectroscopy of a subsample of the source galaxies. We forecast constraints from combining CMB anisotropies with cosmic shear using a spectroscopically determined redshift distribution, varying the number of spectra $N_{spec}$ obtained from 64 to $\infty$. The source redshift distribution is expanded in a Fourier series, and the amplitudes of each mode are considered as parameters to be constrained via both the spectroscopic and weak lensing data. We assume independent source redshifts, and consider in what circumstances this is a good approximation (the sources are clustered and for narrow spectroscopic surveys with many objects this results in the redshifts being correlated). It is found that for the surveys considered and for a prior of 0.04 on the calibration parameters, the addition of redshift information make significant improvements on the constraints on the cosmological parameters; however, beyond $N_{spec}\sim$few$\times 10^3$ the addition of further spectra will make only a very small improvement to the cosmological parameters. We find that a better calibration makes large $N_{spec}$ more useful. Using an eigenvector analysis, we find that the improvement continues with even higher $N_{spec}$, but not in directions that dominate the uncertainties on the standard cosmological parameters.

We use null spherical (observational) coordinates to describe a class of inhomogeneous cosmological models. The proposed cosmological construction is based on the observer past null cone. A known difficulty in using inhomogeneous models is that the null geodesic equation is not integrable in general. Our choice of null coordinates solves the radial ingoing null geodesic by construction. Furthermore, we use an approach where the velocity field is uniquely calculated from the metric rather than put in by hand. Conveniently, this allows us to explore models in a non-comoving frame of reference. In this frame, we find that the velocity field has shear, acceleration and expansion rate in general. We show that a comoving frame is not compatible with expanding perfect fluid models in the coordinates proposed and dust models are simply not possible. We describe the models in a non-comoving frame. We use the dust models in a non-comoving frame to outline a fitting procedure.

We use statistical inference theory to explore the constraints from future galaxy weak lensing (cosmic shear) surveys combined with the current CMB constraints on cosmological parameters, focusing particularly on the running of the spectral index of the primordial scalar power spectrum, $\alpha_s$. Recent papers have drawn attention to the possibility of measuring $\alpha_s$ by combining the CMB with galaxy clustering and/or the Lyman-$\alpha$ forest. Weak lensing combined with the CMB provides an alternative probe of the primordial power spectrum. We run a series of simulations with variable runnings and compare them to semi-analytic non-linear mappings to test their validity for our calculations. We find that a ``Reference'' cosmic shear survey with $f_{sky}=0.01$ and $6.6\times 10^8$ galaxies per steradian can reduce the uncertainty on $n_s$ and $\alpha_s$ by roughly a factor of 2 relative to the CMB alone. We investigate the effect of shear calibration biases on lensing by including the calibration factor as a parameter, and show that for our Reference Survey, the precision of cosmological parameter determination is only slightly degraded even if the amplitude calibration is uncertain by as much as 5%. We conclude that in the near future weak lensing surveys can supplement the CMB observations to constrain the primordial power spectrum.

We show that a flow (timelike congruence) in any type $B_{1}$ warped product spacetime is uniquely and algorithmically determined by the condition of zero flux. (Though restricted, these spaces include many cases of interest.) The flow is written out explicitly for canonical representations of the spacetimes. With the flow determined, we explore an inverse approach to Einstein's equations where a phenomenological fluid interpretation of a spacetime follows directly from the metric irrespective of the choice of coordinates. This approach is pursued for fluids with anisotropic pressure and shear viscosity. In certain degenerate cases this interpretation is shown to be generically not unique. The framework developed allows the study of exact solutions in any frame without transformations. We provide a number of examples, in various coordinates, including spacetimes with and without unique interpretations. The results and algorithmic procedure developed are implemented as a computer algebra program called GRSource.

The Tolman VII solution is an exact static spherically symmetric perfect fluid solution of Einstein's equations that exhibits a surprisingly good approximation to a neutron star. We show that this solution exhibits trapped null orbits in a causal region even for a tenuity (total radius to mass ratio) $> 3$. In this region the dynamical part of the potential for axial w - modes dominates over the centrifugal part.

An explicit necessary condition for the occurrence of resonance scattering of axial gravitational waves, along with the internal trapping of null geodesics, is proposed for static spherically symmetric perfect fluid solutions to Einstein's equations. Some exact inhomogeneous solutions which exhibit this trapping are given with special attention to boundary conditions and the physical acceptability of the space times. In terms of the tenuity ($\alpha = R/M$ at the boundary) all the examples given lie in the narrow range $2.8 < \alpha < 2.9$. The tenuity can be raised to more interesting values by the addition of an envelope without altering the trapping.

GRLite and GRTensorJ are first and second generation graphical user interfaces to the computer algebra system GRTensorII. Current development centers on GRTensorJ, which provides fully customizable symbolic procedures that reduce many complex calculations to "elementary functions". Although still in development, GRTensorJ, which is now available (free of charge) over the internet, is sufficiently advanced to be of interest to researchers in general relativity and related fields.