This paper is the third in a series on tests of gravity using observations of stars and nearby dwarf galaxies. We carry out four distinct tests using published data on the kinematics and morphology of dwarf galaxies, motivated by the theoretical work of Hui et al. (2009) and Jain and Vanderplas (2011). In a wide class of gravity theories a scalar field couples to matter and provides an attractive fifth force. Due to their different self-gravity, stars and gas may respond differently to the scalar force leading to several observable deviations from standard gravity. HI gas, red giant stars and main sequence stars can be displaced relative to each other, and the stellar disk can display warps or asymmetric rotation curves aligned with external potential gradients. To distinguish the effects of modified gravity from standard astrophysical phenomena, we use a control sample of galaxies that are expected to be screened from the fifth force. In all cases we find no significant deviation from the null hypothesis of general relativity. The limits obtained from dwarf galaxies are not yet competitive with the limits from cepheids obtained in our first paper, but can be improved to probe regions of parameter space that are inaccessible using other tests. We discuss how our methodology can be applied to new radio and optical observations of nearby galaxies.

Cosmological galaxy surveys aim at mapping the largest volumes to test models with techniques such as cluster abundance, cosmic shear correlations or baryon acoustic oscillations (BAO), which are designed to be independent of galaxy bias. Here we explore an alternative route to constrain cosmology: sampling more moderate volumes with the cross-correlation of photometric and spectroscopic surveys. We consider the angular galaxy-galaxy autocorrelation in narrow redshift bins and its combination with different probes of weak gravitational lensing (WL) and redshift space distortions (RSD). Including the cross-correlation of these surveys improves by factors of a few the constraints on both the dark energy equation of state w(z) and the cosmic growth history, parametrized by \gamma. The additional information comes from using many narrow redshift bins and from galaxy bias, which is measured both with WL probes and RSD, breaking degeneracies that are present when using each method separately. We show forecasts for a joint w(z) and \gamma figure of merit using linear scales over a deep (i<24) photometric survey and a brighter (i<22.5) spectroscopic or very accurate (0.3%) photometric redshift survey. Magnification or shear in the photometric sample produce FoM that are of the same order of magnitude of those of RSD or BAO over the spectroscopic sample. However, the cross-correlation of these probes over the same area yields a FoM that is up to a factor 100 times larger. Magnification alone, without shape measurements, can also be used for these cross-correlations and can produce better results than using shear alone. For a spectroscopic follow-up survey strategy, measuring the spectra of the foreground lenses to perform this cross-correlation provides 5 times better FoM than targeting the higher redshift tail of the galaxy distribution to study BAO over a 2.5 times larger volume.

A new determination of the sound horizon scale in angular coordinates is presented. It makes use of ~ 0.6 x 10^6 Luminous Red Galaxies, selected from the Sloan Digital Sky Survey imaging data, with photometric redshifts. The analysis covers a redshift interval that goes from z=0.5 to z=0.6. We find evidence of the Baryon Acoustic Oscillations (BAO) signal at the ~ 2.3 sigma confidence level, with a value of theta_{BAO} (z=0.55) = (3.90 \pm 0.38) degrees, including systematic errors. To our understanding, this is the first direct measurement of the angular BAO scale in the galaxy distribution, and it is in agreement with previous BAO measurements. We also show how radial determinations of the BAO scale can break the degeneracy in the measurement of cosmological parameters when they are combined with BAO angular measurements. The result is also in good agreement with the WMAP7 best-fit cosmology. We obtain a value of w_0 = -1.03 \pm 0.16 for the equation of state parameter of the dark energy, Omega_M = 0.26 \pm 0.04 for the matter density, when the other parameters are fixed. We have also tested the sensitivity of current BAO measurements to a time varying dark energy equation of state, finding w_a = 0.06 \pm 0.22 if we fix all the other parameters to the WMAP7 best-fit cosmology.

The possibility of measuring redshift space (RSD) distortions using photometric data have been recently highlighted. This effect complements and significantly alters the detectability of baryon acoustic oscillations (BAO) in photometric surveys. In this paper we present measurements of the angular correlation function of luminous red galaxies (LRGs) in the photometric catalog of the final data release (DR7) of the Sloan Digital Sky Survey II (SDSS). The sample compromise ~ 1.5 x 10^6 LRGs distributed in 0.45 < z < 0.65, with a characteristic photometric error of ~ 0.05. Our measured correlation centered at z=0.55 is in very good agreement with predictions from standard LCDM in a broad range of angular scales, $0.5^\circ < \theta < 6^\circ$. We find that the growth of structure can indeed be robustly measured, with errors matching expectations. The velocity growth rate is recovered as $f \sigma_8 = 0.53 \pm 0.42$ when no prior is imposed on the growth factor and the background geometry follows a LCDM model with WMAP7+SNIa priors. This is compatible with the corresponding General Relativity (GR) prediction $f \sigma_8 = 0.45$ for our fiducial cosmology. If we adopt a parametrization such that $f=\Omega ^\gamma_m(z)$, with $\gamma \approx 0.55$ in GR, and combine our $f\sigma_8$ measurement with the corresponding ones from spectroscopic LRGs at lower redshifts we obtain $\gamma=0.54 \pm 0.17$. In addition we find evidence for the presence of the baryon acoustic feature matching the amplitude, location and shape of LCDM predictions. The photometric BAO feature is detected with 98 % confidence level at z=0.55.

We outline how redshift-space distortions (RSD) can be measured from the angular correlation function w({\theta}), of galaxies selected from photometric surveys. The natural degeneracy between RSD and galaxy bias can be minimized by comparing results from bins with top-hat galaxy selection in redshift, and bins based on the radial position of galaxy pair centres. This comparison can also be used to test the accuracy of the photometric redshifts. The presence of RSD will be clearly detectable with the next generation of photometric redshift surveys. We show that the Dark Energy Survey (DES) will be able to measure f(z){\sigma}_8(z) to a 1{\sigma} accuracy of (17 {\times} b)%, using galaxies drawn from a single narrow redshift slice centered at z = 1. Here b is the linear bias, and f is the logarithmic rate of change of the linear growth rate with respect to the scale factor. Extending to measurements of w({\theta}) for a series of bins of width 0.02(1 + z) over 0.5 < z < 1.4 will measure {\gamma} to a 1{\sigma} accuracy of 25%, given the model f = {\Omega}_m(z)^{\gamma}, and assuming a linear bias model that evolves such that b = 0.5 + z (and fixing other cosmological parameters). The accuracy of our analytic predictions is confirmed using mock catalogs drawn from simulations conducted by the MICE collaboration.

Near future cosmology will see the advent of wide area photometric galaxy surveys, like the Dark Energy Survey (DES), that extent to high redshifts (z ~ 1 - 2) but with poor radial distance resolution. In such cases splitting the data into redshift bins and using the angular correlation function $w(\theta)$, or the $C_{\ell}$ power spectrum, will become the standard approach to extract cosmological information or to study the nature of dark energy through the Baryon Acoustic Oscillations (BAO) probe. In this work we present a detailed model for $w(\theta)$ at large scales as a function of redshift and bin width, including all relevant effects, namely nonlinear gravitational clustering, bias, redshift space distortions and photo-z uncertainties. We also present a model for the full covariance matrix characterizing the angular correlation measurements, that takes into account the same effects as for $w(\theta)$ and also the possibility of a shot-noise component and partial sky coverage. Provided with a large volume N-body simulation from the MICE collaboration we built several ensembles of mock redshift bins with a sky coverage and depth typical of forthcoming photometric surveys. The model for the angular correlation and the one for the covariance matrix agree remarkably well with the mock measurements in all configurations. The prospects for a full shape analysis of $w(\theta)$ at BAO scales in forthcoming photometric surveys such as DES are thus very encouraging.

We propose a new method for cosmological parameters extraction using the baryon acoustic oscillation scale as a standard ruler in deep galaxy surveys with photometric determination of redshifts. The method consists in a simple empirical parametric fit to the angular 2-point correlation function w(theta). It is parametrized as a power law to describe the continuum plus a Gaussian to describe the BAO bump. The location of the Gaussian is used as the basis for the measurement of the sound horizon scale. This method, although simple, actually provides a robust estimation, since the inclusion of the power law and the use of the Gaussian removes the shifts which affect the local maximum. We discuss the effects of projection bias, non-linearities, redshift space distortions and photo-z precision, and apply our method to a mock catalog of the Dark Energy Survey, built upon a large N-body simulation provided by the MICE collaboration. We discuss the main systematic errors associated to our method and show that they are dominated by the photo-z uncertainty.

Recent publications claim that there is no convincing evidence for measurements of the baryonic acoustic (BAO) feature in galaxy samples using either monopole or radial information. Different claims seem contradictory: data is either not consistent with the BAO model or data is consistent with both the BAO model and featureless models without BAO. We investigate this point with a set of 216 realistic mock galaxy catalogs extracted from MICE7680, one of the largest volume dark matter simulation run to date, with a volume of 1300 cubical gigaparsecs. Our mocks cover similar volume, densities and bias as the real galaxies and provide 216 realizations of the Lambda or w=-1 Cold Dark Matter (wCDM) BAO model. We find that only 20% of the mocks show a statistically significant (3 sigma) preference for the true (input) wCDM BAO model as compared to a featureless (non-physical) model without BAO. Thus the volume of current galaxy samples is not yet large enough to claim that the BAO feature has been detected. Does this mean that we can not locate the BAO position? Using a simple (non optimal) algorithm we show that in 50% (100%) of the mocks we can find the BAO position within 5% (20%) of the true value. These two findings are not in contradiction: the former is about model selection, the later is about parameter fitting within a model. We conclude that current monopole and radial BAO measurements can be used as standard rulers if we assume wCDM type of models.

We apply a new model for the spherically averaged correlation function at large pair separations to the measurement of the clustering of luminous red galaxies (LRGs) made from the SDSS by Cabre and Gaztanaga(2009). Our model takes into account the form of the BAO peak and the large scale shape of the correlation function. We perform a Monte Carlo Markov chain analysis for different combinations of datasets and for different parameter sets. When used in combination with a compilation of the latest CMB measurements, the LRG clustering and the latest supernovae results give constraints on cosmological parameters which are comparable and in remarkably good agreement, resolving the tension reported in some studies. The best fitting model in the context of a flat, Lambda-CDM cosmology is specified by Omega_m=0.261+-0.013, Omega_b=0.044+-0.001, n_s=0.96+-0.01, H_0=71.6+-1.2 km/s/Mpc and sigma_8=0.80+-0.02. If we allow the time-independent dark energy equation of state parameter to vary, we find results consistent with a cosmological constant at the 5% level using all data sets: w_DE=-0.97+-0.05. The large scale structure measurements by themselves can constrain the dark energy equation of state parameter to w_DE=-1.05+-0.15, independently of CMB or supernovae data. We do not find convincing evidence for an evolving equation of state. We provide a set of "extended distance priors" that contain the most relevant information from the CMB power spectrum and the shape of the LRG correlation function which can be used to constrain dark energy models and spatial curvature. Our model should provide an accurate description of the clustering even in much larger, forthcoming surveys, such as those planned with NASA's JDEM or ESA's Euclid mission.

We study the clustering of LRG galaxies in the latest spectroscopic SDSS data releases, DR6 and DR7, which sample over 1 Gpc^3/h^3 to z=0.47. The 2-point correlation function $\xisp$ is estimated as a function of perpendicular $\sigma$ and line-of-sight $\pi$ (radial) directions. We find a significant detection of a peak at $r\simeq 110$Mpc/h, which shows as a circular ring in the $\sigma-\pi$ plane. There is also significant evidence for a peak along the radial direction whose shape is consistent with its originating from the recombination-epoch baryon acoustic oscillations (BAO). A $\xisp$ model with no radial BAO peak is disfavored at $3.2\sigma$, whereas a model with no magnification bias is disfavored at $2\sigma$. The radial data enable, for the first time, a direct measurement of the Hubble parameter $H(z)$ as a function of redshift. This is independent from earlier BAO measurements which used the spherically averaged (monopole) correlation to constrain an integral of $H(z)$. Using the BAO peak position as a standard ruler in the radial direction, we find: $H(z=0.24)= 79.69 \pm 2.32 (\pm 1.29)$ km/s/Mpc for z=0.15-0.30 and $H(z=0.43)= 86.45 \pm 3.27 (\pm 1.69)$ km/s/Mpc for $z=0.40-0.47$. The first error is a model independent statistical estimation and the second accounts for systematics both in the measurements and in the model. For the full sample, $z=0.15-0.47$, we find $H(z=0.34)= 83.80 \pm 2.96 (\pm 1.59)$ km/s/Mpc.

We present the 3-point function \xi_3 and Q_3=\xi_3/\xi_2^2 for a spectroscopic sample of luminous red galaxies (LRG) from SDSS DR6 & DR7. We find a strong (S/N$>$6) detection of $Q_3$ on scales of 55-125 Mpc/h, with a well defined peak around 105 Mpc/h in all \xi_2, \xi_3 and Q_3, in excellent agreement with the predicted shape and location of the imprint of the baryon acoustic oscillations (BAO). We use very large simulations to asses and test the significance of our measurement. Models without the BAO peak are ruled out by the $Q_3$ data with 99% confidence. Our measurements show the expected shape for Q_3 as a function of the triangular configuration. This provides a first direct measurement of the non-linear mode coupling coefficients of density and velocity fluctuations which, on these large scales, are independent of cosmic time, the amplitude of fluctuations or cosmological parameters. The location of the BAO peak in the data indicates \Omega_m =0.28 \pm 0.05 and \Omega_B=0.079 \pm 0.025 (h=0.70) after marginalization over spectral index (n_s=0.8-1.2) linear b_1 and quadratic c_2 bias,which are found to be in the range: b_1=1.7-2.2 and c_2=0.75-3.55. The data allows a hierarchical contribution from primordial non-Gaussianities in the range Q_3=0.55-3.35. These constraints are independent and complementary to the ones that can be obtained using the 2-point function, which are presented in a separate paper. This is the first detection of the shape of $Q_3$ on BAO scales, but our errors are shot-noise dominated and the SDSS volume is still relatively small, so there is ample room for future improvement in this type of measurements.

This is the second paper of a series where we study the clustering of LRG galaxies in the latest spectroscopic SDSS data release, DR6, which has 75000 LRG galaxies covering over 1 $Gpc^3/h^3$ for $0.15<z<0.47$. Here we focus on modeling redshift space distortions in $\xisp$, the 2-point correlation in separate line-of-sight and perpendicular directions, at small scales and in the line-of-sight. We show that a simple Kaiser model for the anisotropic 2-point correlation function in redshift space, convolved with a distribution of random peculiar velocities with an exponential form, can describe well the correlation of LRG at all scales. We show that to describe with accuracy the so called "fingers-of-God" (FOG) elongations in the radial direction, it is necessary to model the scale dependence of both bias $b$ and the pairwise rms peculiar velocity $\sigma_{12}$ with the distance. We show how both quantities can be inferred from the $\xisp$ data. From $r \simeq 10$ Mpc/h to $r \simeq 1$ Mpc/h, both the bias and $\sigma_{12}$ are shown to increase by a factor of two: from $b=2$ to $b=4$ and from $\sigma_{12}=400$ to 800 Km/s. The later is in good agreement, within a 5 percent accuracy in the recovered velocities, with direct velocity measurements in dark matter simulations with $\Omega_m=0.25$ and $\sigma_8$=0.85.

In a series of papers we have recently studied the clustering of LRG galaxies in the latest spectroscopic SDSS data release, which has 75000 LRG galaxies sampling 1.1 Gpc^3/h^3 to z=0.47. Here we focus on detecting a local maxima shaped as a circular ring in the bidimensional galaxy correlation function \xi(pi,sigma), separated in perpendicular \sigma and line-of-sight \pi distances. We find a significant detection of such a peak at r ~110 Mpc/h. The overall shape and location of the ring is consistent with it originating from the recombination-epoch baryon acoustic oscillations (BAO). This agreement provides support for the current understanding of how large scale structure forms in the universe. We study the significance of such feature using large mock galaxy simulations to provide accurate errorbars.

This is the first paper of a series where we study the clustering of LRG galaxies in the latest spectroscopic SDSS data release, DR6, which has 75000 LRG galaxies covering over 1 $Gpc^3/h^3$ at $0.15<z<0.47$. Here we focus on modeling redshift space distortions in $\xips$, the 2-point correlation in separate line-of-sight and perpendicular directions, on large scales. % and away from the line-of-sight. We use large mock simulations to study the validity of models and errors. We show that errors in the data are dominated by a shot-noise term that is 40% larger than the Poisson error commonly used. We first use the normalized quadrupole for the whole sample (mean z=0.34) to estimate $\beta=f(\Omega_m)/b=0.34 \pm 0.03$, where $f(\Omega_m)$ is the linear velocity growth factor and $b$ is the linear bias parameter that relates galaxy to matter fluctuations on large scales. We next use the full $\xips$ plane to find $\Omega_{0m}= 0.245 \pm 0.020$ (h=0.72) and the biased amplitude $b \sigma_8 = 1.56 \pm 0.09$. For standard gravity, we can combine these measurements to break degeneracies and find $\sigma_8=0.85 \pm 0.06$, $b=1.85 \pm 0.25$ and $f(\Omega_m)=0.64 \pm 0.09$. We present constraints for modified theories of gravity and find that standard gravity is consistent with data as long as $0.80<\sigma_8<0.92$. We also calculate the cross-correlation with WMAP5 and show how both methods to measure the growth history are complementary to constrain non-standard models of gravity. Finally, we show results for different redshift slices, including a prominent BAO peak in the monopole at different redshifts. (Abridged)

Constraining cosmological parameters from measurements of the Integrated Sachs-Wolfe effect requires developing robust and accurate methods for computing statistical errors in the cross-correlation between maps. This paper presents a detailed comparison of such error estimation applied to the case of cross-correlation of Cosmic Microwave Background (CMB) and large-scale structure data. We compare theoretical models for error estimation with montecarlo simulations where both the galaxy and the CMB maps vary around a fiducial auto-correlation and cross-correlation model which agrees well with the current concordance LCDM cosmology. Our analysis compares estimators both in harmonic and configuration (or real) space, quantifies the accuracy of the error analysis and discuss the impact of partial sky survey area and the choice of input fiducial model on dark-energy constraints. We show that purely analytic approaches yield accurate errors even in surveys that cover only 10% of the sky and that parameter constraints strongly depend on the fiducial model employed. Alternatively, we discuss the advantages and limitations of error estimators that can be directly applied to data. In particular, we show that errors and covariances from the Jack-Knife method agree well with the theoretical approaches and simulations. We also introduce a novel method in real space that is computationally efficient and can be applied to real data and realistic survey geometries. Finally, we present a number of new findings and prescriptions that can be useful for analysis of real data and forecasts, and present a critical summary of the analyses done to date.

We cross-correlate the third-year WMAP data with galaxy samples extracted from the SDSS DR4 (SDSS4) covering 13% of the sky, increasing by a factor of 3.7 the volume sampled in previous analyses. The new measurements confirm a positive cross-correlation with higher significance (total signal-to-noise of about 4.7). The correlation as a function of angular scale is well fitted by the integrated Sachs-Wolfe (ISW) effect for LCDM flat FRW models with a cosmological constant. The combined analysis of different samples gives Omega_L=0.80-0.85$ (68% Confidence Level, CL) or $0.77-0.86$ (95% CL). We find similar best fit values for Omega_L for different galaxy samples with median redshifts of z ~0.3 and z ~0.5, indicating that the data scale with redshift as predicted by the LCDM cosmology (with equation of state parameter w=-1). This agreement is not trivial, but can not yet be used to break the degeneracy constraints in the w versus Omega_L plane using only the ISW data.