Source author record

Gong-Bo Zhao

Gong-Bo Zhao appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

102works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

102 published item(s)

preprint2022arXiv

A new consistency test for $Λ$CDM cosmology using galaxy surveys

We propose a new consistency test for the $Λ$CDM cosmology using baryonic acoustic oscillations (BAO) and redshift space distortion (RSD) measurements from galaxy redshift surveys. Specifically, we determine the peak position of $fσ_8(z)$ in redshift $z$ offered by a RSD measurement, and compare it to the one predicted by the BAO observables assuming a flat $Λ$CDM cosmology. We demonstrate this new test using the simulated data for the DESI galaxy survey, and argue that this test complements those using the background observables alone, and is less subject to systematics in the RSD analysis, compared to traditional methods using values of $fσ_8(z)$ directly.

preprint2022arXiv

Combined full shape analysis of BOSS galaxies and eBOSS quasars using an iterative emulator

Standard full-shape clustering analyses in Fourier space rely on a fixed power spectrum template, defined at the fiducial cosmology used to convert redshifts into distances, and compress the cosmological information into the Alcock-Paczynski parameters and the linear growth rate of structure. In this paper, we propose an analysis method that operates directly in the cosmology parameter space and varies the power spectrum template accordingly at each tested point. Predictions for the power spectrum multipoles from the TNS model are computed at different cosmologies in the framework of $Λ\rm{CDM}$. Applied to the final eBOSS QSO and LRG samples together with the low-z DR12 BOSS galaxy sample, our analysis results in a set of constraints on the cosmological parameters $Ω_{\rm cdm}$, $H_0$, $σ_8$, $Ω_{\rm b}$ and $n_s$. To reduce the number of computed models, we construct an iterative process to sample the likelihood surface, where each iteration consists of a Gaussian process regression. This method is validated with mocks from N-body simulations. From the combined analysis of the (e)BOSS data, we obtain the following constraints: $σ_8=0.877\pm 0.049$ and $Ω_{\rm m}=0.304^{+0.016}_{-0.010}$ without any external prior. The eBOSS quasar sample alone shows a $3.1σ$ discrepancy compared to the Planck prediction.

preprint2022arXiv

Cosmology Intertwined: A Review of the Particle Physics, Astrophysics, and Cosmology Associated with the Cosmological Tensions and Anomalies

In this paper we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant $H_0$, the $σ_8$--$S_8$ tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the $5.0\,σ$ tension between the {\it Planck} CMB estimate of the Hubble constant $H_0$ and the SH0ES collaboration measurements. After showing the $H_0$ evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the {\it Planck} CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density $Ω_m$, and the amplitude or rate of the growth of structure ($σ_8,fσ_8$). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the $H_0$--$S_8$ tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals.[Abridged]

preprint2022arXiv

Galaxy power spectrum and biasing results from the LOFAR Two-metre Sky Survey (first data release)

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing survey aiming to observe the entire Northern sky, providing an excellent opportunity to study the distribution and evolution of the large-scale structure of the Universe. The source catalogue from the public LoTSS first data release (DR1) covers 1% of the sky, and shows correlated noise or fluctuations of the flux density calibration on few degree scales. We explore the LoTSS DR1 to understand the survey systematics and data quality of this first data release. We produce catalog mocks to estimate uncertainties, and measure the angular clustering statistics of LoTSS galaxies, which fit the $Λ$CDM cosmology reasonably well. We employ a Markov chain Monte Carlo (MCMC) based Bayesian analysis to recover the best galaxy biasing scheme and multi-component source fraction for LoTSS DR1 above $1$ mJy assuming different possible redshift templates. After masking some noisy and uneven patches and with suitable flux density cuts, the LOFAR survey appears qualified for large-scale cosmological studies. The upcoming data releases from LOFAR are expected to be deeper and wider, and will therefore provide improved cosmological measurements.

preprint2022arXiv

The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Cosmological implications from multi-tracer BAO analysis with galaxies and voids

We construct cosmic void catalogues with the DIVE void finder upon SDSS BOSS DR12 and eBOSS DR16 galaxy samples with BAO reconstruction applied, and perform a joint BAO analysis using different types of galaxies and the corresponding voids. The BAO peak is evident for the galaxy-galaxy, galaxy-void, and void-void correlation functions of all datasets, including the ones cross correlating luminous red galaxy and emission line galaxy samples. Two multi-tracer BAO fitting schemes are then tested, one combining the galaxy and void correlation functions with a weight applied to voids, and the other using a single BAO dilation parameter for all clustering measurements of different tracers. Both methods produce consistent results with mock catalogues, and on average ~10 per cent improvements of the BAO statistical uncertainties are observed for all samples, compared to the results from galaxies alone. By combining the clustering of galaxies and voids, the uncertainties of BAO measurements from the SDSS data are reduced by 5 to 15 per cent, yielding 0.9, 0.8, 1.1, 2.3, and 2.9 per cent constraints on the distance $D_{_{\rm V}}(z)$, at effective redshifts of 0.38, 0.51, 0.70, 0.77, and 0.85, respectively. When combined with BAO measurements from SDSS MGS, QSO, and Ly$α$ samples, as well as the BBN results, we obtain $H_0 = 67.58 \pm 0.91\,{\rm km}\,{\rm s}^{-1}\,{\rm Mpc}^{-1}$, $Ω_{\rm m} = 0.290 \pm 0.015$, and $Ω_Λh^2 = 0.3241 \pm 0.0079$ in the flat-$Λ$CDM framework, where the 1$\,σ$ uncertainties are around 6, 6, and 17 per cent smaller respectively, compared to constraints from the corresponding anisotropic BAO measurements without voids and LRG-ELG cross correlations.

preprint2021arXiv

The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: measurement of the BAO and growth rate of structure of the emission line galaxy sample from the anisotropic power spectrum between redshift 0.6 and 1.1

We analyse the large-scale clustering in Fourier space of emission line galaxies (ELG) from the Data Release 16 of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey. The ELG sample contains 173,736 galaxies covering 1,170 square degrees in the redshift range $0.6 < z < 1.1$. We perform a BAO measurement from the post-reconstruction power spectrum monopole, and study redshift space distortions (RSD) in the first three even multipoles. Photometric variations yield fluctuations of both the angular and radial survey selection functions. Those are directly inferred from data, imposing integral constraints which we model consistently. The full data set has only a weak preference for a BAO feature ($1.4σ$). At the effective redshift $z_{\rm eff} = 0.845$ we measure $D_{\rm V}(z_{\rm eff})/r_{\rm drag} = 18.33_{-0.62}^{+0.57}$, with $D_{\rm V}$ the volume-averaged distance and $r_{\rm drag}$ the comoving sound horizon at the drag epoch. In combination with the RSD measurement, at $z_{\rm eff} = 0.85$ we find $fσ_8(z_{\rm eff}) = 0.289_{-0.096}^{+0.085}$, with $f$ the growth rate of structure and $σ_8$ the normalisation of the linear power spectrum, $D_{\rm H}(z_{\rm eff})/r_{\rm drag} = 20.0_{-2.2}^{+2.4}$ and $D_{\rm M}(z_{\rm eff})/r_{\rm drag} = 19.17 \pm 0.99$ with $D_{\rm H}$ and $D_{\rm M}$ the Hubble and comoving angular distances, respectively. These results are in agreement with those obtained in configuration space, thus allowing a consensus measurement of $fσ_8(z_{\rm eff}) = 0.315 \pm 0.095$, $D_{\rm H}(z_{\rm eff})/r_{\rm drag} = 19.6_{-2.1}^{+2.2}$ and $D_{\rm M}(z_{\rm eff})/r_{\rm drag} = 19.5 \pm 1.0$. This measurement is consistent with a flat $Λ$CDM model with Planck parameters.

preprint2021arXiv

The design of the Ali CMB Polarization Telescope receiver

Ali CMB Polarization Telescope (AliCPT-1) is the first CMB degree-scale polarimeter to be deployed on the Tibetan plateau at 5,250m above sea level. AliCPT-1 is a 90/150 GHz 72 cm aperture, two-lens refracting telescope cooled down to 4 K. Alumina lenses, 800mm in diameter, image the CMB in a 33.4° field of view on a 636mm wide focal plane. The modularized focal plane consists of dichroic polarization-sensitive Transition-Edge Sensors (TESes). Each module includes 1,704 optically active TESes fabricated on a 150mm diameter silicon wafer. Each TES array is read out with a microwave multiplexing readout system capable of a multiplexing factor up to 2,048. Such a large multiplexing factor has allowed the practical deployment of tens of thousands of detectors, enabling the design of a receiver that can operate up to 19 TES arrays for a total of 32,376 TESes. AliCPT-1 leverages the technological advancements in the detector design from multiple generations of previously successful feedhorn-coupled polarimeters, and in the instrument design from BICEP-3, but applied on a larger scale. The cryostat receiver is currently under integration and testing. During the first deployment year, the focal plane will be populated with up to 4 TES arrays. Further TES arrays will be deployed in the following years, fully populating the focal plane with 19 arrays on the fourth deployment year. Here we present the AliCPT-1 receiver design, and how the design has been optimized to meet the experimental requirements.

preprint2020arXiv

An accurate analytic mass model for lensing galaxies

We develop an analytic mass model for lensing galaxies, based on a broken power-law (BPL) density profile, which is a power-law profile with a mass deficit or surplus in the central region. Under the assumption of an elliptically symmetric surface mass distribution, the deflection angle and magnification can be evaluated analytically for this new model. We compute the theoretical prediction for various quantities, including the volume and surface mass density profiles of the galaxies, and the aperture and luminosity-weighted line-of-sight velocity dispersions, and compare them to those measured from the Illustris simulation. We find an excellent agreement between our model prediction and the simulation, which validates our modeling. The high efficiency and accuracy of our model manifests itself as a promising tool for studying properties of galaxies with strong lensing.

preprint2020arXiv

Investigating the degeneracy between modified gravity and massive neutrinos with redshift-space distortions

There is a well known degeneracy between the enhancement of the growth of large-scale structure produced by modified gravity models and the suppression due to the free-streaming of massive neutrinos at late times. This makes the matter power-spectrum alone a poor probe to distinguish between modified gravity and the concordance $Λ$CDM model when neutrino masses are not strongly constrained. In this work, we investigate the potential of using redshift-space distortions (RSD) to break this degeneracy when the modification to gravity is scale-dependent in the form of Hu-Sawicki $f(R)$. We find that if the linear growth rate can be recovered from the RSD signal, the degeneracy can be broken at the level of the dark matter field. However, this requires accurate modelling of the non-linearities in the RSD signal, and we here present an extension of the standard perturbation theory-based model for non-linear RSD that includes both Hu-Sawicki $f(R)$ modified gravity and massive neutrinos.

preprint2020arXiv

The clustering of the SDSS-IV extended Baryon Oscillation Spectroscopic Survey DR14 LRG sample: structure growth rate measurement from the anisotropic LRG correlation function in the redshift range 0.6 < z < 1.0

We analyze the anisotropic clustering of the Sloan Digital Sky Survey-IV Extended Baryon Oscillation Spectroscopic Survey (eBOSS) Luminous Red Galaxy Data Release 14 (DR14) sample combined with Baryon Oscillation Spectroscopic Survey (BOSS) CMASS sample of galaxies in the redshift range 0.6$<z<$1.0, which consists of 80,118 galaxies from eBOSS and 46,439 galaxies from the BOSS-CMASS sample. The eBOSS-CMASS Luminous Red Galaxy sample has a sky coverage of 1,844 deg$^2$, with an effective volume of 0.9 Gpc$^3$. The analysis was made in configuration space using a Legendre multipole expansion. The Redshift Space Distortion signal is modeled as a combination of the Convolution Lagrangian Perturbation Model and the Gaussian Streaming Model. We constrain the logarithmic growth of structure times the amplitude of dark matter density fluctuations, $f (z_{\rm eff})σ_8(z_{\rm eff})=0.454 \pm0.139 $, and the Alcock-Paczynski dilation scales which constraints the angular diameter distance $D_A(z_{eff})=1466.5 \pm 136.6 (r_s/r_s^{\rm fid})$ and $H(z_{\rm eff})=105.8 \pm 16 (r_s^{\rm fid}/r_s) \mathrm{km\,s^{-1}\,Mpc^{-1}}$, where $r_s$ is the sound horizon at the end of the baryon drag epoch and $r_s^{\rm fid}$ is its value in the fiducial cosmology at an effective redshift $z_{\rm eff}=0.72$. These results are in full agreement with the current $Λ$-Cold Dark Matter ($Λ$-CDM) cosmological model inferred from Planck measurements. This study is the first eBOSS LRG full-shape analysis i.e. including Redshift-Space Distortions (RSD) simultaneously with the Alcock-Paczynski (AP) effect and the Baryon Acoustic Oscillation (BAO) scale.

preprint2020arXiv

The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: BAO and RSD measurements from anisotropic clustering analysis of the Quasar Sample in configuration space between redshift 0.8 and 2.2

We measure the anisotropic clustering of the quasar sample from Data Release 16 (DR16) of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey (eBOSS). A sample of $343,708$ spectroscopically confirmed quasars between redshift $0.8<z<2.2$ are used as tracers of the underlying dark matter field. In comparison with DR14 sample, the final sample doubles the number of objects as well as the survey area. In this paper, we present the analysis in configuration space by measuring the two-point correlation function and decompose using the Legendre polynomials. For the full-shape analysis of the Legendre multipole moments, we measure the BAO distance and the growth rate of the cosmic structure. At an effective redshift of $z_{\rm eff}=1.48$, we measure the comoving angular diameter distance $D_{\rm M}(z_{\rm eff})/r_{\rm drag} = 30.66\pm0.88$, the Hubble distance $D_{\rm H}(z_{\rm eff})/r_{\rm drag} = 13.11\pm0.52$, and the growth rate $fσ_8(z_{\rm eff}) = 0.439\pm0.048$. The accuracy of these measurements is confirmed using an extensive set of mock simulations developed for the quasar sample. The uncertainties on the distance and growth rate measurements have been reduced substantially ($\sim 45\%$ and $\sim30\%$) with respect to the DR14 results. We also perform a BAO-only analysis to cross check the robustness of the methodology of the full-shape analysis. Combining our analysis with the Fourier space analysis, we arrive at $D^{\bf{c}}_{\rm M}(z_{\rm eff})/r_{\rm drag} = 30.22 \pm 0.79$, $D^{\bf{c}}_{\rm H}(z_{\rm eff})/r_{\rm drag} = 13.26 \pm 0.47$, and $fσ_8^{\bf{c}}(z_{\rm eff}) = 0.464 \pm 0.045$.

preprint2020arXiv

The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: measurement of the BAO and growth rate of structure of the luminous red galaxy sample from the anisotropic correlation function between redshifts 0.6 and 1

We present the cosmological analysis of the configuration-space anisotropic clustering in the completed Sloan Digital Sky Survey IV (SDSS-IV) extended Baryon Oscillation Spectroscopic Survey (eBOSS) DR16 galaxy sample. This sample consists of luminous red galaxies (LRGs) spanning the redshift range $0.6 < z < 1$, at an effective redshift of $z_{\rm eff}=0.698$. It combines 174 816 eBOSS LRGs and 202 642 BOSS CMASS galaxies. We extract and model the baryon acoustic oscillations (BAO) and redshift-space distortions (RSD) features from the galaxy two-point correlation function to infer geometrical and dynamical cosmological constraints. The adopted methodology is extensively tested on a set of realistic simulations. The correlations between the inferred parameters from the BAO and full-shape correlation function analyses are estimated. This allows us to derive joint constraints on the three cosmological parameter combinations: $D_M(z)/r_d$, $D_H(z)/r_d$ and $fσ_8(z)$, where $D_M$ is the comoving angular diameter distance, $D_H$ is Hubble distance, $r_d$ is the comoving BAO scale, $f$ is the linear growth rate of structure, and $σ_8$ is the amplitude of linear matter perturbations. After combining the results with those from the parallel power spectrum analysis of Gil-Marin et al. 2020, we obtain the constraints: $D_M/r_d = 17.65 \pm 0.30$, $D_H/r_d = 19.77 \pm 0.47$, $fσ_8 = 0.473 \pm 0.044$. These measurements are consistent with a flat $Λ$CDM model with standard gravity.

preprint2020arXiv

The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Pairwise-Inverse-Probability and Angular Correction for Fibre Collisions in Clustering Measurements

The completed eBOSS catalogues contain redshifts of 344080 QSOs over 0.8<z<2.2 covering 4808 deg$^2$, 174816 LRGs over 0.6<z<1.0 covering 4242 deg$^2$ and 173736 ELGs over 0.6<z<1.1 covering 1170 deg$^2$ in order to constrain the expansion history of the Universe and the growth rate of structure through clustering measurements. Mechanical limitations of the fibre-fed spectrograph on the Sloan telescope prevent two fibres being placed closer than 62", the fibre-collision scale, in a single pass of the instrument on the sky. These `fibre collisions' strongly correlate with the intrinsic clustering of targets and can bias measurements of the two-point correlation function resulting in a systematic error on the inferred values of the cosmological parameters. We combine the new techniques of pairwise-inverse-probability weighting and the angular up-weighting to correct the clustering measurements for the effect of fibre collisions. Using mock catalogues we show that our corrections provide unbiased measurements, within data precision, of both the projected correlation function $w_p$ and the multipoles $ξ^l$ of the redshift-space correlation functions down to 0.1Mpc/h, regardless of the tracer type. We apply the corrections to the eBOSS DR16 catalogues. We find that, on scales greater than s~20Mpc/h for $ξ^l$, as used to make BAO and large-scale RSD measurements, approximate methods such as Nearest-Neighbour up-weighting are sufficiently accurate given the statistical errors of the data. Using the PIP method, for the first time for a spectroscopic program of the Sloan Digital Sky Survey we are able to successfully access the 1-halo term in the 3D clustering measurements down to ~0.1Mpc/h scales. Our results will therefore allow studies that use the small-scale clustering measurements to strengthen the constraints on both cosmological parameters and the halo-occupation distribution models.

preprint2019arXiv

Emulators for the non-linear matter power spectrum beyond $Λ$CDM

Accurate predictions for the non-linear matter power spectrum are needed to confront theory with observations in current and near future weak lensing and galaxy clustering surveys. We propose a computationally cheap method to create an emulator for modified gravity models by utilizing existing emulators for $Λ$CDM. Using a suite of $N$-body simulations we construct a fitting function for the enhancement of both the linear and non-linear matter power spectrum in the commonly studied Hu-Sawicki $f(R)$ gravity model valid for wave-numbers $k \lesssim 5-10\, h\text{Mpc}^{-1}$ and redshifts $z \lesssim 3$. We show that the cosmology dependence of this enhancement is relatively weak so that our fit, using simulations coming from only one cosmology, can be used to get accurate predictions for other cosmological parameters. We also show that the cosmology dependence can, if needed, be included by using linear theory, approximate $N$-body simulations (such as COLA) and semi-analytical tools like the halo model. Our final fit can easily be combined with any emulator or semi-analytical models for the non-linear $Λ$CDM power spectrum to accurately, and quickly, produce a non-linear power spectrum for this particular modified gravity model. The method we use can be applied to fairly cheaply construct an emulator for other modified gravity models. As an application of our fitting formula we use it to compute Fisher-forecasts for how well galaxy clustering and weak lensing in a Euclid-like survey will be at constraining modifications of gravity.

preprint2019arXiv

Generalized Brans-Dicke theories in light of evolving dark energy

The expansion history of the Universe reconstructed from a combination of recent data indicates a preference for a changing Dark Energy (DE) density. Moreover, the DE density appears to be increasing with cosmic time, with its equation of state being below -1 on average, and possibly crossing the so-called phantom divide. Scalar-tensor theories, in which the scalar field mediates a force between matter particles, offer a natural framework in which the effective DE equation of state can be less than -1 and cross the phantom barrier. We consider the generalized Brans-Dicke (GBD) class of scalar-tensor theories and reconstruct their Lagrangian given the effective DE density extracted from recent data. Then, given the reconstructed Lagrangian, we solve for the linear perturbations and investigate the characteristic signatures of these reconstructed GBD in the cosmological observables, such as the cosmic microwave background (CMB) anisotropy, the galaxy number counts, and their cross-correlations. In particular, we demonstrate that the Integrated Sachs-Wolfe (ISW) effect probed by the cross-correlation of CMB with the matter distribution can rule out scalar-tensor theories as the explanation of the observed DE dynamics independently from the laboratory and solar system fifth force constraints.

preprint2016arXiv

Acausality in Nonlocal Gravity Theory

We investigate the nonlocal gravity theory by deriving nonlocal equations of motion using the traditional variation principle in a homogeneous background. We focus on a class of models with a linear nonlocal modification term in the action. It is found that the resulting equations of motion contain the advanced Green's function, implying that there is an acausality problem. As a consequence, a divergence arises in the solutions due to contributions from the future infinity unless the Universe will go back to the radiation dominated era or become the Minkowski spacetime in the future. We also discuss the relation between the original nonlocal equations and its biscalar-tensor representation and identify the auxiliary fields with the corresponding original nonlocal terms. Finally, we show that the acusality problem cannot be avoided by any function of nonlocal terms in the action.

preprint2016arXiv

Constraining $f(R)$ Gravity Theory Using Weak Lensing Peak Statistics from the Canada-France-Hawaii-Telescope Lensing Survey

In this Letter, we report the observational constraints on the Hu-Sawicki $f(R)$ theory derived from weak lensing peak abundances, which are closely related to the mass function of massive halos. In comparison with studies using optical or x-ray clusters of galaxies, weak lensing peak analyses have the advantages of not relying on mass-baryonic observable calibrations. With observations from the Canada-France-Hawaii-Telescope Lensing Survey, our peak analyses give rise to a tight constraint on the model parameter $|f_{R0}|$ for $n=1$. The $95\%$ CL limit is $\log_{10}|f_{R0}| < -4.82$ given WMAP9 priors on $(Ω_{\rm m}, A_{\rm s})$. With Planck15 priors, the corresponding result is $\log_{10}|f_{R0}| < -5.16$.

preprint2016arXiv

Probing Theories of Gravity with Phase Space-Inferred Potentials of Galaxy Clusters

Modified theories of gravity provide us with a unique opportunity to generate innovative tests of gravity. In Chameleon f(R) gravity, the gravitational potential differs from the weak-field limit of general relativity (GR) in a mass dependent way. We develop a probe of gravity which compares high mass clusters, where Chameleon effects are weak, to low mass clusters, where the effects can be strong. We utilize the escape velocity edges in the radius/velocity phase space to infer the gravitational potential profiles on scales of 0.3-1 virial radii. We show that the escape edges of low mass clusters are enhanced compared to GR, where the magnitude of the difference depends on the background field value |fR0|. We validate our probe using N-body simulations and simulated light cone galaxy data. For a DESI (Dark Energy Spectroscopic Instrument) Bright Galaxy Sample, including observational systematics, projection effects, and cosmic variance, our test can differentiate between GR and Chameleon f(R) gravity models, |fR0| = 4e-6 (2e-6) at > 5σ (> 2σ), more than an order of magnitude better than current cluster-scale constraints.

preprint2016arXiv

Screening fifth forces in generalized Proca theories

For a massive vector field with derivative self-interactions, the breaking of the gauge invariance allows the propagation of a longitudinal mode in addition to the two transverse modes. We consider generalized Proca theories with second-order equations of motion in a curved space-time and study how the longitudinal scalar mode of the vector field gravitates on a spherically symmetric background. We show explicitly that cubic-order self-interactions lead to the suppression of the longitudinal mode through the Vainshtein mechanism. Provided that the dimensionless coupling of the interaction is not negligible, this screening mechanism is sufficiently efficient to give rise to tiny corrections to gravitational potentials consistent with solar-system tests of gravity. We also study the quartic interactions with the presence of non-minimal derivative coupling with the Ricci scalar and find the existence of solutions where the longitudinal mode completely vanishes. Finally, we discuss the case in which the effect of the quartic interactions dominates over the cubic one and show that local gravity constraints can be satisfied under a mild bound on the parameters of the theory.

preprint2016arXiv

Signatures of the Primordial Universe from Its Emptiness: Measurement of Baryon Acoustic Oscillations from Minima of the Density Field

Sound waves from the primordial fluctuations of the Universe imprinted in the large-scale structure, called baryon acoustic oscillations (BAOs), can be used as standard rulers to measure the scale of the Universe. These oscillations have already been detected in the distribution of galaxies. Here we propose to measure BAOs from the troughs (minima) of the density field. Based on two sets of accurate mock halo catalogues with and without BAOs in the seed initial conditions, we demonstrate that the BAO signal cannot be obtained from the clustering of classical disjoint voids, but is clearly detected from overlapping voids. The latter represent an estimate of all troughs of the density field. We compute them from the empty circumsphere centers constrained by tetrahedra of galaxies using Delaunay triangulation. Our theoretical models based on an unprecedented large set of detailed simulated void catalogues are remarkably well confirmed by observational data. We use the largest recently publicly available sample of luminous red galaxies from SDSS-III BOSS DR11 to unveil for the first time a >3$σ$ BAO detection from voids in observations. Since voids are nearly isotropically expanding regions, their centers represent the most quiet places in the Universe, keeping in mind the cosmos origin and providing a new promising window in the analysis of the cosmological large-scale structure from galaxy surveys.

preprint2016arXiv

Simulation tests of galaxy cluster constraints on chameleon gravity

We use two new hydrodynamical simulations of $Λ$CDM and $f(R)$ gravity to test the methodology used by Wilcox et al. 2015 (W15) in constraining the effects of a fifth force on the profiles of clusters of galaxies. We construct realistic simulated stacked weak lensing and X-ray surface brightness cluster profiles from these cosmological simulations, and then use these data projected along various lines-of-sight to test the spherical symmetry of our stacking procedure. We also test the applicability of the NFW profile to model weak lensing profiles of clusters in $f(R)$ gravity. Finally, we test the validity of the analytical model developed in W15 against the simulated profiles. Overall, we find our methodology is robust and broadly agrees with these simulated data. We also apply our full Markov Chain Monte Carlo (MCMC) analysis from W15 to our simulated X-ray and lensing profiles, providing consistent constraints on the modified gravity parameters as obtained from the real cluster data, e.g. for our $Λ$CDM simulation we obtain $|f_{\rm{R}0}| < 8.3 \times 10^{-5}$ (95% CL), which is in good agreement with the W15 measurement of $|f_{\rm{R}0}| < 6 \times 10^{-5}$. Overall, these tests confirm the power of our methodology which can now be applied to larger cluster samples available with the next generation surveys.

preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: combining correlated Gaussian posterior distributions

The cosmological information contained in anisotropic galaxy clustering measurements can often be compressed into a small number of parameters whose posterior distribution is well described by a Gaussian. We present a general methodology to combine these estimates into a single set of consensus constraints that encode the total information of the individual measurements, taking into account the full covariance between the different methods. We illustrate this technique by applying it to combine the results obtained from different clustering analyses, including measurements of the signature of baryon acoustic oscillations (BAO) and redshift-space distortions (RSD), based on a set of mock catalogues of the final SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). Our results show that the region of the parameter space allowed by the consensus constraints is smaller than that of the individual methods, highlighting the importance of performing multiple analyses on galaxy surveys even when the measurements are highly correlated. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The methodology presented here is used in Alam et al. (2016) to produce the final cosmological constraints from BOSS.

preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological implications of the configuration-space clustering wedges

We explore the cosmological implications of anisotropic clustering measurements in configuration space of the final galaxy samples from Data Release 12 of the SDSS-III Baryon Oscillation Spectroscopic Survey. We implement a new detailed modelling of the effects of non-linearities, galaxy bias and redshift-space distortions that can be used to extract unbiased cosmological information from our measurements for scales $s \gtrsim 20\,h^{-1}{\rm Mpc}$. We combined the galaxy clustering information from BOSS with the latest cosmic microwave background (CMB) observations and Type Ia supernovae samples and found no significant evidence for a deviation from the $Λ$CDM cosmological model. In particular, these data sets can constrain the dark energy equation of state parameter to $w_{\rm DE}=-0.996\pm0.042$ when assumed time-independent, the curvature of the Universe to $Ω_{k}=-0.0007\pm 0.0030$ and the sum of the neutrino masses to $\sum m_ν < 0.25\,{\rm eV}$ at 95 per cent CL. We explore the constraints on the growth rate of cosmic structures assuming $f(z)=Ω_{\rm m}(z)^γ$ and obtain $γ= 0.609\pm 0.079$, in good agreement with the predictions of general relativity of $γ=0.55$. We compress the information of our clustering measurements into constraints on the parameter combinations $D_{\rm V}(z)/r_{\rm d}$, $F_{\rm AP}(z)$ and $fσ_8(z)$ at the effective redshifts of $z=0.38$, $0.51$ and $0.61$ with their respective covariance matrices and find good agreement with the predictions for these parameters obtained from the best-fitting $Λ$CDM model to the CMB data from the Planck satellite. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. (2016) to produce the final cosmological constraints from BOSS.

preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: double-probe measurements from BOSS galaxy clustering \& Planck data -- towards an analysis without informative priors

We develop a new methodology called double-probe analysis with the aim of minimizing informative priors in the estimation of cosmological parameters. We extract the dark-energy-model-independent cosmological constraints from the joint data sets of Baryon Oscillation Spectroscopic Survey (BOSS) galaxy sample and Planck cosmic microwave background (CMB) measurement. We measure the mean values and covariance matrix of $\{R$, $l_a$, $Ω_b h^2$, $n_s$, $log(A_s)$, $Ω_k$, $H(z)$, $D_A(z)$, $f(z)σ_8(z)\}$, which give an efficient summary of Planck data and 2-point statistics from BOSS galaxy sample, where $R=\sqrt{Ω_m H_0^2}\,r(z_*)$, and $l_a=πr(z_*)/r_s(z_*)$, $z_*$ is the redshift at the last scattering surface, and $r(z_*)$ and $r_s(z_*)$ denote our comoving distance to $z_*$ and sound horizon at $z_*$ respectively. The advantage of this method is that we do not need to put informative priors on the cosmological parameters that galaxy clustering is not able to constrain well, i.e. $Ω_b h^2$ and $n_s$. Using our double-probe results, we obtain $Ω_m=0.304\pm0.009$, $H_0=68.2\pm0.7$, and $σ_8=0.806\pm0.014$ assuming $Λ$CDM; and $Ω_k=0.002\pm0.003$ and $w=-1.00\pm0.07$ assuming o$w$CDM. The results show no tension with the flat $Λ$CDM cosmological paradigm. By comparing with the full-likelihood analyses with fixed dark energy models, we demonstrate that the double-probe method provides robust cosmological parameter constraints which can be conveniently used to study dark energy models. We extend our study to measure the sum of neutrino mass and obtain $Σm_ν<0.10/0.22$ (68\%/95\%) assuming $Λ$CDM and $Σm_ν<0.26/0.52$ (68\%/95\%) assuming $w$CDM. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS.

preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: Observational systematics and baryon acoustic oscillations in the correlation function

We present baryon acoustic oscillation (BAO) scale measurements determined from the clustering of 1.2 million massive galaxies with redshifts 0.2 < z < 0.75 distributed over 9300 square degrees, as quantified by their redshift-space correlation function. In order to facilitate these measurements, we define, describe, and motivate the selection function for galaxies in the final data release (DR12) of the SDSS III Baryon Oscillation Spectroscopic Survey (BOSS). This includes the observational footprint, masks for image quality and Galactic extinction, and weights to account for density relationships intrinsic to the imaging and spectroscopic portions of the survey. We simulate the observed systematic trends in mock galaxy samples and demonstrate that they impart no bias on baryon acoustic oscillation (BAO) scale measurements and have a minor impact on the recovered statistical uncertainty. We measure transverse and radial BAO distance measurements in 0.2 < z < 0.5, 0.5 < z < 0.75, and (overlapping) 0.4 < z < 0.6 redshift bins. In each redshift bin, we obtain a precision that is 2.7 per cent or better on the radial distance and 1.6 per cent or better on the transverse distance. The combination of the redshift bins represents 1.8 per cent precision on the radial distance and 1.1 per cent precision on the transverse distance. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. (2016) to produce the final cosmological constraints from BOSS.

preprint2016arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: BAO measurement from the LOS-dependent power spectrum of DR12 BOSS galaxies

[abridged] We present an anisotropic analysis of the baryonic acoustic oscillation (BAO) scale in the twelfth and final data release of the Baryonic Oscillation Spectroscopic Survey (BOSS). We independently analyse the LOWZ and CMASS galaxy samples: the LOWZ sample contains contains 361 762 galaxies with an effective redshift of $z_{\rm LOWZ}=0.32$; the CMASS sample consists of 777 202 galaxies with an effective redshift of $z_{\rm CMASS}=0.57$. We extract the BAO peak position from the monopole power spectrum moment, $α_0$, and from the $μ^2$ moment, $α_2$, where $μ$ is the cosine of the angle to the line-of-sight. The $μ^2$-moment provides equivalent information to that available in the quadrupole but is simpler to analyse. After applying a reconstruction algorithm to reduce the BAO suppression by bulk motions, we measure the BAO peak position in the monopole and $μ^2$-moment, which are related to radial and angular shifts in scale. We report $H(z_{\rm LOWZ})r_s(z_d)=(11.60\pm0.60)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm LOWZ})/r_s(z_d)=6.66\pm0.16$ with a cross-correlation coefficient of $r_{HD_A}=0.41$, for the LOWZ sample; and $H(z_{\rm CMASS})r_s(z_d)=(14.56\pm0.37)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm CMASS})/r_s(z_d)=9.42\pm0.13$ with a cross-correlation coefficient of $r_{HD_A}=0.47$, for the CMASS sample. We combine these results with the measurements of the BAO peak position in the monopole and quadrupole correlation function of the same dataset \citep[][companion paper]{Cuestaetal2015} and report the consensus values: $H(z_{\rm LOWZ})r_s(z_d)=(11.63\pm0.69)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm LOWZ})/r_s(z_d)=6.67\pm0.15$ with $r_{HD_A}=0.35$ for the LOWZ sample; $H(z_{\rm CMASS})r_s(z_d)=(14.67\pm0.42)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm CMASS})/r_s(z_d)=9.47\pm0.12$ with $r_{HD_A}=0.52$ for the CMASS sample.

preprint2016arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: RSD measurement from the LOS-dependent power spectrum of DR12 BOSS galaxies

We measure and analyse the clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) relative to the line-of-sight (LOS), for LOWZ and CMASS galaxy samples drawn from the final Data Release 12 (DR12). The LOWZ sample contains 361\,762 galaxies with an effective redshift of $z_{\rm lowz}=0.32$, and the CMASS sample 777\,202 galaxies with an effective redshift of $z_{\rm cmass}=0.57$. From the power spectrum monopole and quadrupole moments around the LOS, we measure the growth of structure parameter $f$ times the amplitude of dark matter density fluctuations $σ_8$ by modeling the Redshift-Space Distortion signal. When the geometrical Alcock-Paczynski effect is also constrained from the same data, we find joint constraints on $fσ_8$, the product of the Hubble constant and the comoving sound horizon at the baryon drag epoch $H(z)r_s(z_d)$, and the angular distance parameter divided by the sound horizon $D_A(z)/r_s(z_d)$. We find $f(z_{\rm lowz})σ_8(z_{\rm lowz})=0.394\pm0.062$, $D_A(z_{\rm lowz})/r_s(z_d)=6.35\pm0.19$, $H(z_{\rm lowz})r_s(z_d)=(11.41\pm 0.56)\,{10^3\rm km}s^{-1}$ for the LOWZ sample, and $f(z_{\rm cmass})σ_8(z_{\rm cmass})=0.444\pm0.038$, $D_A(z_{\rm cmass})/r_s(z_d)=9.42\pm0.15$, $H(z_{\rm cmass})r_s(z_d)=(13.92 \pm 0.44)\, {10^3\rm km}s^{-1}$ for the CMASS sample. We find general agreement with previous BOSS DR11 measurements. Assuming the Hubble parameter and angular distance parameter are fixed at fiducial $Λ$CDM values, we find $f(z_{\rm lowz})σ_8(z_{\rm lowz})=0.485\pm0.044$ and $f(z_{\rm cmass})σ_8(z_{\rm cmass})=0.436\pm0.022$ for the LOWZ and CMASS samples, respectively.

preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$α$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

preprint2016arXiv

The DESI Experiment Part II: Instrument Design

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= λ/Δλ$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

preprint2016arXiv

The extended Baryon Oscillation Spectroscopic Survey (eBOSS): a cosmological forecast

We present a science forecast for the eBOSS survey, part of the SDSS-IV project, which is a spectroscopic survey using multiple tracers of large-scale structure, including luminous red galaxies (LRGs), emission line galaxies (ELGs) and quasars (both as a direct probe of structure and through the Ly-$α$ forest). Focusing on discrete tracers, we forecast the expected accuracy of the baryonic acoustic oscillation (BAO), the redshift-space distortion (RSD) measurements, the $f_{\rm NL}$ parameter quantifying the primordial non-Gaussianity, the dark energy and modified gravity parameters. We also use the line-of-sight clustering in the Ly-$α$ forest to constrain the total neutrino mass. We find that eBOSS LRGs ($0.6<z<1.0$) (combined with the BOSS LRGs at $z>0.6$), ELGs ($0.6<z<1.2$) and Clustering Quasars (CQs) ($0.6<z<2.2$) can achieve a precision of 1%, 2.2% and 1.6% precisions, respectively, for spherically averaged BAO distance measurements. Using the same samples, the constraint on $fσ_8$ is expected to be 2.5%, 3.3% and 2.8% respectively. For primordial non-Gaussianity, eBOSS alone can reach an accuracy of $σ(f_{\rm NL})\sim10-15$, depending on the external measurement of the galaxy bias and our ability to model large-scale systematic errors. eBOSS can at most improve the dark energy Figure of Merit (FoM) by a factor of $3$ for the Chevallier-Polarski-Linder (CPL) parametrisation, and can well constrain three eigenmodes for the general equation-of-state parameter (Abridged).

preprint2016arXiv

The SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Overview and Early Data

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

preprint2015arXiv

An efficient probe of the cosmological CPT violation

We develop an efficient method based on the linear regression algorithm to probe the cosmological CPT violation using the CMB polarisation data. We validate this method using simulated CMB data and apply it to recent CMB observations. We find that a combined data sample of BICEP1 and BOOMERanG 2003 favours a nonzero isotropic rotation angle at $2.3σ$ confidence level, ie, $Δα=-3.3 \pm1.4$ deg (68% CL) with systematics included.

preprint2015arXiv

Consistent Modified Gravity Analysis of Anisotropic Galaxy Clustering Using BOSS DR11

We analyse the clustering of cosmic large scale structure using a consistent modified gravity perturbation theory, accounting for anisotropic effects along and transverse to the line of sight. The growth factor has a particular scale dependence in f(R) gravity and we fit for the shape parameter f_{R0} simultaneously with the distance and the large scale (general relativity) limit of the growth function. Using more than 690,000 galaxies in the Baryon Oscillation Spectroscopy Survey Data Release 11, we find no evidence for extra scale dependence, with the 95\% confidence upper limit |f_{R0}| <8 \times 10^{-4}. Future clustering data, such as from the Dark Energy Spectroscopic Instrument, can use this consistent methodology to impose tighter constraints.

preprint2015arXiv

Cosmic Web and Environmental Dependence of Screening: Vainshtein vs. Chameleon

Theories which modify general relativity to explain the accelerated expansion of the Universe often use screening mechanisms to satisfy constraints on Solar System scales. We investigate the effects of the cosmic web and the local environmental density of dark matter halos on the screening properties of the Vainshtein and chameleon screening mechanisms. We compare the cosmic web morphology of dark matter particles, mass functions of dark matter halos, mass and radial dependence of screening, velocity dispersions and peculiar velocities, and environmental dependence of screening mechanisms in $f(R)$ and nDGP models. Using the ORIGAMI cosmic web identification routine we find that the Vainshtein mechanism depends on the cosmic web morphology of dark matter particles, since these are defined according to the dimensionality of their collapse, while the chameleon mechanism shows no morphology dependence. The chameleon screening of halos and their velocity dispersions depend on halo mass, and small halos and subhalos can be environmentally screened in the chameleon mechanism. On the other hand, the screening of halos in the Vainshtein mechanism does not depend on mass nor environment, and their velocity dispersions are suppressed. The peculiar velocities of halos in the Vainshtein mechanism are enhanced because screened objects can still feel the fifth force generated by external fields, while peculiar velocities of chameleon halos are suppressed when the halo centers are screened.

preprint2015arXiv

Cosmological implications of baryon acoustic oscillation (BAO) measurements

We derive constraints on cosmological parameters and tests of dark energy models from the combination of baryon acoustic oscillation (BAO) measurements with cosmic microwave background (CMB) and Type Ia supernova (SN) data. We take advantage of high-precision BAO measurements from galaxy clustering and the Ly-alpha forest (LyaF) in the BOSS survey of SDSS-III. BAO data alone yield a high confidence detection of dark energy, and in combination with the CMB angular acoustic scale they further imply a nearly flat universe. Combining BAO and SN data into an "inverse distance ladder" yields a 1.7% measurement of $H_0=67.3 \pm1.1$ km/s/Mpc. This measurement assumes standard pre-recombination physics but is insensitive to assumptions about dark energy or space curvature, so agreement with CMB-based estimates that assume a flat LCDM cosmology is an important corroboration of this minimal cosmological model. For open LCDM, our BAO+SN+CMB combination yields $Ω_m=0.301 \pm 0.008$ and curvature $Ω_k=-0.003 \pm 0.003$. When we allow more general forms of evolving dark energy, the BAO+SN+CMB parameter constraints remain consistent with flat LCDM. While the overall $χ^2$ of model fits is satisfactory, the LyaF BAO measurements are in moderate (2-2.5 sigma) tension with model predictions. Models with early dark energy that tracks the dominant energy component at high redshifts remain consistent with our constraints. Expansion history alone yields an upper limit of 0.56 eV on the summed mass of neutrino species, improving to 0.26 eV if we include Planck CMB lensing. Standard dark energy models constrained by our data predict a level of matter clustering that is high compared to most, but not all, observational estimates. (Abridged)

preprint2015arXiv

f(R) gravity on non-linear scales: The post-Friedmann expansion and the vector potential

Many modified gravity theories are under consideration in cosmology as the source of the accelerated expansion of the universe and linear perturbation theory, valid on the largest scales, has been examined in many of these models. However, smaller non-linear scales offer a richer phenomenology with which to constrain modified gravity theories. Here, we consider the Hu-Sawicki form of $f(R)$ gravity and apply the post-Friedmann approach to derive the leading order equations for non-linear scales, i.e. the equations valid in the Newtonian-like regime. We reproduce the standard equations for the scalar field, gravitational slip and the modified Poisson equation in a coherent framework. In addition, we derive the equation for the leading order correction to the Newtonian regime, the vector potential. We measure this vector potential from $f(R)$ N-body simulations at redshift zero and one, for two values of the $f_{R_0}$ parameter. We find that the vector potential at redshift zero in $f(R)$ gravity can be close to 50\% larger than in GR on small scales for $|f_{R_0}|=1.289\times10^{-5}$, although this is less for larger scales, earlier times and smaller values of the $f_{R_0}$ parameter. Similarly to in GR, the small amplitude of this vector potential suggests that the Newtonian approximation is highly accurate for $f(R)$ gravity, and also that the non-linear cosmological behaviour of $f(R)$ gravity can be completely described by just the scalar potentials and the $f(R)$ field.

preprint2015arXiv

Mass-concentration relation of clusters of galaxies from CFHTLenS

Based on weak lensing data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), in this paper we study the mass-concentration ($M$-$c$) relation for $\sim 200$ redMaPPer clusters in the fields. We extract the $M$-$c$ relation by measuring the density profiles of individual clusters instead of using stacked weak lensing signals. By performing Monte Carlo simulations, we demonstrate that although the signal-to-noise ratio for each individual cluster is low, the unbiased $M$-$c$ relation can still be reliably derived from a large sample of clusters by carefully taking into account the impacts of shape noise, cluster center offset, dilution effect from member or foreground galaxies, and the projection effect. Our results show that within error bars the derived $M$-$c$ relation for redMaPPer clusters is in agreement with simulation predictions. There is a weak deviation in that the halo concentrations calibrated by Monte Carlo simulations are somewhat higher than that predicted from ${\it Planck}$ cosmology.

preprint2015arXiv

Measuring redshift-space distortions with future SKA surveys

The peculiar motion of galaxies can be a particularly sensitive probe of gravitational collapse. As such, it can be used to measure the dynamics of dark matter and dark energy as well the nature of the gravitational laws at play on cosmological scales. Peculiar motions manifest themselves as an overall anisotropy in the measured clustering signal as a function of the angle to the line-of-sight, known as redshift-space distortion (RSD). Limiting factors in this measurement include our ability to model non-linear galaxy motions on small scales and the complexities of galaxy bias. The anisotropy in the measured clustering pattern in redshift-space is also driven by the unknown distance factors at the redshift in question, the Alcock-Paczynski distortion. This weakens growth rate measurements, but permits an extra geometric probe of the Hubble expansion rate. In this chapter we will briefly describe the scientific background to the RSD technique, and forecast the potential of the SKA phase 1 and the SKA2 to measure the growth rate using both galaxy catalogues and intensity mapping, assessing their competitiveness with current and future optical galaxy surveys.

preprint2015arXiv

Model-independent constraints on dark energy and modified gravity with the SKA

Employing a nonparametric approach of the principal component analysis (PCA), we forecast the future constraint on the equation of state $w(z)$ of dark energy, and on the effective Newton constant $μ(k,z)$, which parameterise the effect of modified gravity, using the planned SKA HI galaxy survey. Combining with the simulated data of Planck and Dark Energy Survey (DES), we find that SKA Phase 1 (SKA1) and SKA Phase 2 (SKA2) can well constrain $3$ and $5$ eigenmodes of $w(z)$ respectively. The errors of the best measured modes can be reduced to 0.04 and 0.023 for SKA1 and SKA2 respectively, making it possible to probe dark energy dynamics. On the other hand, SKA1 and SKA2 can constrain $7$ and $20$ eigenmodes of $μ(k,z)$ respectively within 10\% sensitivity level. Furthermore, 2 and 7 modes can be constrained within sub percent level using SKA1 and SKA2 respectively. This is a significant improvement compared to the combined datasets without SKA.

preprint2015arXiv

Modified Gravity N-body Code Comparison Project

Self-consistent ${\it N}$-body simulations of modified gravity models are a key ingredient to obtain rigorous constraints on deviations from General Relativity using large-scale structure observations. This paper provides the first detailed comparison of the results of different ${\it N}$-body codes for the $f(R)$, DGP, and Symmetron models, starting from the same initial conditions. We find that the fractional deviation of the matter power spectrum from $Λ$CDM agrees to better than $1\%$ up to $k \sim 5-10~h/{\rm Mpc}$ between the different codes. These codes are thus able to meet the stringent accuracy requirements of upcoming observational surveys. All codes are also in good agreement in their results for the velocity divergence power spectrum, halo abundances and halo profiles. We also test the quasi-static limit, which is employed in most modified gravity ${\it N}$-body codes, for the Symmetron model for which the most significant non-static effects among the models considered are expected. We conclude that this limit is a very good approximation for all of the observables considered here.

preprint2015arXiv

Reconstruction of the dark matter-vacuum energy interaction

An interaction between the vacuum energy and dark matter is an intriguing possibility which may offer a way of solving the cosmological constant problem. Adopting a general prescription for momentum exchange between the two dark components, we reconstruct $α(a)$, the temporal evolution of the coupling strength between dark matter and vacuum energy, in a nonparametric Bayesian approach using combined observational data sets from the cosmic microwave background, supernovae and large scale structure. An evolving interaction between the vacuum energy and dark matter removes some of the tensions between different data sets. However, it is not preferred over $Λ$CDM in the Bayesian sense, as improvement in the fit is not sufficient to compensate for the increase in the volume of the parameter space.

preprint2015arXiv

Searching for scalar gravitational interactions in current and future cosmological data

Modified gravity theories often contain a scalar field of gravitational strength which interacts with matter. We examine constraints on the range and the coupling strength of a scalar gravitational degree of freedom using a subset of current data that can be safely analyzed within the linear perturbation theory. Using a model-independent implementation of scalar-tensor theories in MGCAMB in terms of two functions of the scale factor describing the mass and the coupling of the scalar degree of freedom, we derive constraints on the $f(R)$, generalized chameleon, Symmetron and Dilaton models. Since most of the large scale structure data available today is from relatively low redshifts, only a limited range of observed scales is in the linear regime, leading to relatively weak constraints. We then perform a forecast for a future large scale structure survey, such as Large Synoptic Survey Telescope (LSST), which will map a significant volume at higher redshifts, and show that it will produce much stronger constraints on scalar interactions in specific models. We also perform a principal component analysis and find that future surveys should be able to provide tight constraints on several eigenmodes of the scalar mass evolution.

preprint2015arXiv

Synergy between the Large Synoptic Survey Telescope and the Square Kilometre Array

We provide an overview of the science benefits of combining information from the Square Kilometre Array (SKA) and the Large Synoptic Survey Telescope (LSST). We first summarise the capabilities and timeline of the LSST and overview its science goals. We then discuss the science questions in common between the two projects, and how they can be best addressed by combining the data from both telescopes. We describe how weak gravitational lensing and galaxy clustering studies with LSST and SKA can provide improved constraints on the causes of the cosmological acceleration. We summarise the benefits to galaxy evolution studies of combining deep optical multi-band imaging with radio observations. Finally, we discuss the excellent match between one of the most unique features of the LSST, its temporal cadence in the optical waveband, and the time resolution of the SKA.

preprint2015arXiv

The Eleventh and Twelfth Data Releases of the Sloan Digital Sky Survey: Final Data from SDSS-III

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

preprint2015arXiv

The XMM Cluster Survey: Testing chameleon gravity using the profiles of clusters

The chameleon gravity model postulates the existence of a scalar field that couples with matter to mediate a fifth force. If it exists, this fifth force would influence the hot X-ray emitting gas filling the potential wells of galaxy clusters. However, it would not influence the clusters' weak lensing signal. Therefore, by comparing X-ray and weak lensing profiles, one can place upper limits on the strength of a fifth force. This technique has been attempted before using a single, nearby cluster (Coma, $z=0.02$). Here we apply the technique to the stacked profiles of 58 clusters at higher redshifts ($0.1<z<1.2$), including 12 new to the literature, using X-ray data from the XMM Cluster Survey (XCS) and weak lensing data from the Canada France Hawaii Telescope Lensing Survey (CFHTLenS). Using a multi-parameter MCMC analysis, we constrain the two chameleon gravity parameters ($β$ and $ϕ_{\infty}$). Our fits are consistent with general relativity, not requiring a fifth force. In the special case of $f(R)$ gravity (where $β= \sqrt{1/6}$), we set an upper limit on the background field amplitude today of $|f_{\rm{R0}}| < 6 \times 10^{-5}$ (95% CL). This is one of the strongest constraints to date on $|f_{\rm{R0}}|$ on cosmological scales. We hope to improve this constraint in future by extending the study to hundreds of clusters using data from the Dark Energy Survey.

preprint2014arXiv

Dark energy imprints on the kinematic Sunyaev-Zel'dovich signal

We investigate the imprint of dark energy on the kinetic Sunyaev-Zel'dovich (kSZ) angular power spectrum on scales of $\ell=1000$ to $10000$, and find that the kSZ signal is sensitive to the dark energy parameter. For example, varying the constant $w$ by 20\% around $w=-1$ results in a $\gtrsim10\%$ change on the kSZ spectrum; changing the dark energy dynamics parametrized by $w_a$ by $\pm0.5$, a 30\% change on the kSZ spectrum is expected. We discuss the observational aspects and develop a fitting formula for the kSZ power spectrum. Finally, we discuss how the precise modeling of the post-reionization signal would help the constraints on patchy reionization signal, which is crucial for measuring the duration of reionization.

preprint2014arXiv

Modeling the nonlinear clustering in modified gravity models I: A fitting formula for matter power spectrum of f(R) gravity

Based on a suite of N-body simulations of the Hu-Sawicki model of f(R) gravity with different sets of model and cosmological parameters, we develop a new fitting formula with a numeric code, MGHalofit, to calculate the nonlinear matter power spectrum P(k) for the Hu-Sawicki model. We compare the MGHalofit predictions at various redshifts (z<=1) to the f(R) simulations and find that the accuracy on P(k) is 6% at k<1 h/Mpc and 12% at 1<k<10 h/Mpc respectively. Based on a sensitivity study of an ongoing and a future spectroscopic survey, we estimate the detectability of a signal of modified gravity described by the Hu-Sawicki model using the power spectrum up to quasi-nonlinear scales. MGHalofit is publicly available at http://icosmology.info/website/MGHalofit.html.

preprint2014arXiv

Observable physical modes of modified gravity

At linear order in cosmological perturbations, departures from the growth in the cosmological standard model can be quantified in terms of two functions of redshift $z$ and Fourier number $k$. Previous studies have performed principal component forecasts for several choices of these two functions, based on expected capabilities of upcoming large structure surveys. It is typically found that there will be many well-constrained degrees of freedom. However, not all and, probably most, of these degrees of freedom were physical if the parametrization had allowed for an arbitrary $k$-dependence. In this paper, we restrict the $k$-dependence to that allowed in local theories of gravity under the quasi-static approximation, i.e. ratios of polynomials in $k$, and identify the best constrained features in the ($z$,$k$)-dependence of the commonly considered functions $μ$ and $γ$ as measured by an LSST-like weak lensing survey. We estimate the uncertainty in the measurements of the eigenmodes of modified growth. We find that imposing the theoretical prior on $k$-dependence reduces the number of degrees of freedom and the covariance between parameters. On the other hand, imaging surveys like LSST are not as sensitive to the $z$-dependence as they are to the $k$-dependence of the modified growth functions. This trade off provides us with, more or less, the same number of well-constrained eigenmodes (with respect to our prior) as found before, but now these modes are physical.

preprint2014arXiv

Post-$Planck$ constraints on interacting vacuum energy

We present improved constraints on an interacting vacuum model using updated astronomical observations including the first data release from Planck. We consider a model with one dimensionless parameter, $α$, describing the interaction between dark matter and vacuum energy (with fixed equation of state $w=-1$). The background dynamics correspond to a generalised Chaplygin gas cosmology, but the perturbations have a zero sound speed. The tension between the value of the Hubble constant, $H_0$, determined by Planck data plus WMAP polarisation (Planck+WP) and that determined by the Hubble Space Telescope (HST) can be alleviated by energy transfer from dark matter to vacuum ($α>0$). A positive $α$ increases the allowed values of $H_0$ due to parameter degeneracy within the model using only CMB data. Combining with additional datasets of including supernova type Ia (SN Ia) and baryon acoustic oscillation (BAO), we can significantly tighten the bounds on $α$. Redshift-space distortions (RSD), which constrain the linear growth of structure, provide the tightest constraints on vacuum interaction when combined with Planck+WP, and prefer energy transfer from vacuum to dark matter ($α<0$) which suppresses the growth of structure. Using the combined datasets of Planck+WP+Union2.1+BAO+RSD, we obtain the constraint on $α$ to be $-0.083<α<-0.006$ (95% C.L.), allowing low $H_0$ consistent with the measurement from 6dF Galaxy survey. This interacting vacuum model can alleviate the tension between RSD and Planck+WP in the $Λ$CDM model for $α<0$, or between HST measurements of $H_0$ and Planck+WP for $α>0$, but not both at the same time.

preprint2014arXiv

Reconstructing primordial power spectrum using Planck and SDSS-III measurements

We develop an accurate and efficient Bayesian method to reconstruct the primordial power spectrum in a model-independent way, and apply it to the latest cosmic microwave background measurement from Planck mission, and the large scale structure observation of SDSS-III BOSS (CMASS) sample, combined with the type Ia supernovae sample (SNLS 3-year) and the measurements of baryon acoustic oscillations from SDSS-II, 6dF, and WiggleZ survey. We confirm that the scale-invariant primordial power spectrum is strongly disfavored, and a model with suppressed power on horizon scales is supported by current data. We also find that a modulation on scales $5\times10^{-4}~\textrm{Mpc}^{-1} \lesssim k \lesssim 0.01~\textrm{Mpc}^{-1}$ is mildly preferred at $2σ$ confidence level, whose origin needs further investigation.

preprint2014arXiv

The Clustering of Galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS): measuring growth rate and geometry with anisotropic clustering

We use the observed anisotropic clustering of galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 11 CMASS sample to measure the linear growth rate of structure, the Hubble expansion rate and the comoving distance scale. Our sample covers 8498 ${\rm deg}^2$ and encloses an effective volume of 6.0 ${\rm Gpc}^3$ at an effective redshift of $\bar{z} = 0.57$. We find $fσ_8 = 0.441 \pm 0.044$, $H = 93.1 \pm 3.0\ {\mathrm{km}\ \mathrm{s}^{-1} \mathrm{Mpc}^{-1}}$ and $D_{\rm A} = 1380 \pm 23\ {\rm Mpc}$ when fitting the growth and expansion rate simultaneously. When we fix the background expansion to the one predicted by spatially-flat $Λ$CDM model in agreement with recent Planck results, we find $fσ_8 = 0.447 \pm 0.028$ (6 per cent accuracy). While our measurements are generally consistent with the predictions of $Λ$CDM and General Relativity, they mildly favor models in which the strength of gravitational interactions is weaker than what is predicted by General Relativity. Combining our measurements with recent cosmic microwave background data results in tight constraints on basic cosmological parameters and deviations from the standard cosmological model. Separately varying these parameters, we find $w = -0.983 \pm 0.075$ (8 per cent accuracy) and $γ= 0.69 \pm 0.11$ (16 per cent accuracy) for the effective equation of state of dark energy and the growth rate index, respectively. Both constraints are in good agreement with the standard model values of $w=-1$ and $γ= 0.554$.

preprint2014arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the Data Release 10 and 11 galaxy samples

We present a one per cent measurement of the cosmic distance scale from the detections of the baryon acoustic oscillations in the clustering of galaxies from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). Our results come from the Data Release 11 (DR11) sample, containing nearly one million galaxies and covering approximately $8\,500$ square degrees and the redshift range $0.2<z<0.7$. We also compare these results with those from the publicly released DR9 and DR10 samples. Assuming a concordance $Λ$CDM cosmological model, the DR11 sample covers a volume of 13\,Gpc${}^3$ and is the largest region of the Universe ever surveyed at this density. We measure the correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of over $7\,σ$ in both the correlation function and power spectrum. Fitting for the position of the acoustic features measures the distance relative to the sound horizon at the drag epoch, $r_d$, which has a value of $r_{d,{\rm fid}}=149.28\,$Mpc in our fiducial cosmology. We find $D_V=(1264\pm25\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.32$ and $D_V=(2056\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.57$. At 1.0 per cent, this latter measure is the most precise distance constraint ever obtained from a galaxy survey. Separating the clustering along and transverse to the line-of-sight yields measurements at $z=0.57$ of $D_A=(1421\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ and $H=(96.8\pm3.4\,{\rm km/s/Mpc})(r_{d,{\rm fid}}/r_d)$. Our measurements of the distance scale are in good agreement with previous BAO measurements and with the predictions from cosmic microwave background data for a spatially flat cold dark matter model with a cosmological constant.

preprint2014arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: weighing the neutrino mass using the galaxy power spectrum of the CMASS sample

We measure the sum of the neutrino particle masses using the three-dimensional galaxy power spectrum of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) CMASS galaxy sample. Combined with the cosmic microwave background (CMB), supernova (SN) and additional baryonic acoustic oscillation (BAO) data, we find upper 95 percent confidence limits of the neutrino mass $Σm_ν<0.340$ eV within a flat $Λ$CDM background, and $Σm_ν<0.821$ eV, assuming a more general background cosmological model. The number of neutrino species is measured to be $N_{\rm eff}=4.308\pm0.794$ and $N_{\rm eff}=4.032^{+0.870}_{-0.894}$ for these two cases respectively. We study and quantify the effect of several factors on the neutrino measurements, including the galaxy power spectrum bias model, the effect of redshift-space distortion, the cutoff scale of the power spectrum, and the choice of additional data. The impact of neutrinos with unknown masses on other cosmological parameter measurements is investigated. The fractional matter density and the Hubble parameter are measured to be $Ω_M=0.2796\pm0.0097$, $H_0=69.72^{+0.90}_{-0.91}$ km/s/Mpc (flat $Λ$CDM) and $Ω_M=0.2798^{+0.0132}_{-0.0136}$, $H_0=73.78^{+3.16}_{-3.17}$ km/s/Mpc (more general background model). Based on a Chevallier-Polarski-Linder (CPL) parametrisation of the equation-of-state $w$ of dark energy, we find that $w=-1$ is consistent with observations, even allowing for neutrinos. Similarly, the curvature Ω_K and the running of the spectral index $α_s$ are both consistent with zero. The tensor-to-scaler ratio is constrained down to $r<0.198$ (95 percent CL, flat $Λ$ CDM) and $r<0.440$ (95 percent CL, more general background model).

preprint2014arXiv

The Vainshtein Mechanism in the Cosmic Web

We investigate the dependence of the Vainshtein screening mechanism on the cosmic web morphology of both dark matter particles and halos as determined by ORIGAMI. Unlike chameleon and symmetron screening, which come into effect in regions of high density, Vainshtein screening instead depends on the dimensionality of the system, and screened bodies can still feel external fields. ORIGAMI is well-suited to this problem because it defines morphologies according to the dimensionality of the collapsing structure and does not depend on a smoothing scale or density threshold parameter. We find that halo particles are screened while filament, wall, and void particles are unscreened, and this is independent of the particle density. However, after separating halos according to their large scale morphological environment, we find no difference in the screening properties of halos in filaments versus halos in clusters. We find that the fifth force enhancement of dark matter particles in halos is greatest well outside the virial radius. We confirm the theoretical expectation that even if the internal field is suppressed by the Vainshtein mechanism, the object still feels the fifth force generated by the external fields, by measuring peculiar velocities and velocity dispersions of halos. Finally, we investigate the morphology and gravity model dependence of halo spins, concentrations, and shapes.

preprint2013arXiv

Chameleon f(R) gravity in the virialized cluster

Current constraints on f(R) gravity from the large-scale structure are at the verge of penetrating into a region where the modified forces become nonlinearly suppressed. For a consistent treatment of observables at these scales, we study cluster quantities produced in chameleon and linearized Hu-Sawicki f(R) gravity dark matter N-body simulations. We find that the standard Navarro-Frenk-White halo density profile and the radial power law for the pseudo phase-space density provide equally good fits for f(R) clusters as they do in the Newtonian scenario. We give qualitative arguments for why this should be the case. For practical applications, we derive analytic relations, e.g., for the f(R) scalar field, the gravitational potential, and the velocity dispersion as seen within the virialized clusters. These functions are based on three degrees of freedom fitted to simulations, i.e., the characteristic density, scale, and velocity dispersion. We further analyze predictions for these fitting parameters from the gravitational collapse and the Jeans equation, which are found to agree well with the simulations. Our analytic results can be used to consistently constrain chameleon f(R) gravity with future observations on virialized cluster scales without the necessity of running a large number of simulations.

preprint2013arXiv

Cosmology with Photometrically-Classified Type Ia Supernovae from the SDSS-II Supernova Survey

We present the cosmological analysis of 752 photometrically-classified Type Ia Supernovae (SNe Ia) obtained from the full Sloan Digital Sky Survey II (SDSS-II) Supernova (SN) Survey, supplemented with host-galaxy spectroscopy from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). Our photometric-classification method is based on the SN typing technique of Sako et al. (2011), aided by host galaxy redshifts (0.05<z<0.55). SNANA simulations of our methodology estimate that we have a SN Ia typing efficiency of 70.8%, with only 3.9% contamination from core-collapse (non-Ia) SNe. We demonstrate that this level of contamination has no effect on our cosmological constraints. We quantify and correct for our selection effects (e.g., Malmquist bias) using simulations. When fitting to a flat LambdaCDM cosmological model, we find that our photometric sample alone gives omega_m=0.24+0.07-0.05 (statistical errors only). If we relax the constraint on flatness, then our sample provides competitive joint statistical constraints on omega_m and omega_lambda, comparable to those derived from the spectroscopically-confirmed three-year Supernova Legacy Survey (SNLS3). Using only our data, the statistics-only result favors an accelerating universe at 99.96% confidence. Assuming a constant wCDM cosmological model, and combining with H0, CMB and LRG data, we obtain w=-0.96+0.10-0.10, omega_m=0.29+0.02-0.02 and omega_k=0.00+0.03-0.02 (statistical errors only), which is competitive with similar spectroscopically confirmed SNe Ia analyses. Overall this comparison is re-assuring, considering the lower redshift leverage of the SDSS-II SN sample (z<0.55) and the lack of spectroscopic confirmation used herein. These results demonstrate the potential of photometrically-classified SNe Ia samples in improving cosmological constraints.

preprint2013arXiv

Exploring Vainshtein mechanism on adaptively refined meshes

There has been a lot of research interest in modified gravity theories which utilise the Vainshtein mechanism to recover standard general relativity in regions with high matter density, such as the Dvali-Gabadadze-Porrati and Galileon models. The strong nonlinearity in the field equations of these theories implies that accurate theoretical predictions could only be made using high-resolution cosmological simulations. Previously, such simulations were usually done on regular meshes, which limits both their performance and the accuracy. In this paper, we report the development of a new algorithm and code, based on ECOSMOG, that uses adaptive mesh refinements to improve the efficiency and precision in simulating the models with Vainshtein mechanism. We have made various code tests against the numerical reliability, and found consistency with previous simulations. We also studied the velocity field in the self-accelerating branch of the DGP model. The code, parallelised using MPI, is suitable for large cosmological simulations of Galileon-type modified gravity theories.

preprint2013arXiv

Modeling halo mass functions in chameleon f(R) gravity

On cosmological scales, observations of the cluster abundance currently place the strongest constraints on f(R) gravity. These constraints lie in the large-field limit, where the modifications of general relativity can correctly be modeled by setting the Compton wavelength of the scalar field to its background value. These bounds are, however, at the verge of penetrating into a regime, where the modifications become nonlinearly suppressed due to the chameleon mechanism and cannot be described by this linearized approximation. For future constraints based on observations subjected to cluster abundance, it is therefore essential to consistently model the chameleon effect. We analyze descriptions of the halo mass function in chameleon f(R) gravity using a mass- and environment-dependent spherical collapse model in combination with excursion set theory and phenomenological fits to N-body simulations in the LCDM and f(R) gravity scenarios. Our halo mass functions consistently incorporate the chameleon suppression and cosmological parameter dependencies, improving upon previous formalisms and providing an important extension to N-body simulations for the application in consistent tests of gravity with observables sensitive to the abundance of clusters.

preprint2013arXiv

Principal Component Analysis of Modified Gravity using Weak Lensing and Peculiar Velocity Measurements

We perform a principal component analysis to assess ability of future observations to measure departures from General Relativity in predictions of the Poisson and anisotropy equations on linear scales. In particular, we focus on how the measurements of redshift-space distortions (RSD) observed from spectroscopic galaxy redshift surveys will improve the constraints when combined with lensing tomographic surveys. Assuming a Euclid-like galaxy imaging and redshift survey, we find that adding the 3D information decreases the statistical uncertainty by a factor between 3 and 10 compared to the case when only observables from lensing tomographic surveys are used. We also find that the number of well-constrained modes increases by a factor between 3 and 7. Our study indicates the importance of joint galaxy imaging and redshift surveys such as SuMIRe and Euclid to give more stringent tests of the ΛCDM model and to distinguish between various modified gravity and dark energy models.

preprint2013arXiv

Relativistic corrections and non-Gaussianity in radio continuum surveys

Forthcoming radio continuum surveys will cover large volumes of the observable Universe and will reach to high redshifts, making them potentially powerful probes of dark energy, modified gravity and non-Gaussianity. Here we extend recent works by analyzing the general relativistic (GR) corrections to the angular power spectrum. These GR corrections to the standard Newtonian analysis of the power spectrum become significant on scales near and beyond the Hubble scale at each redshift. We consider the continuum surveys with LOFAR, WSRT and ASKAP, and examples of continuum surveys with the SKA. We find that the GR corrections are at most percent-level in LOFAR, WODAN and EMU surveys, but they can produce $O(10%)$ changes for high enough sensitivity SKA continuum surveys. The signal is however dominated by cosmic variance, and multiple-tracer observations will be needed to overcome this problem. The GR corrections are suppressed in continuum surveys because of the integration over redshift -- we expect that GR corrections will be enhanced for future SKA HI surveys in which the source redshifts will be known. We also provide predictions for the angular power spectra in the case where the primordial perturbations have local non-Gaussianity. We find that non-Gaussianity dominates over GR corrections, and rises above cosmic variance when $f_{\rm NL}\gtrsim5$ for SKA continuum surveys.

preprint2013arXiv

Simulating the quartic Galileon gravity model on adaptively refined meshes

We develop a numerical algorithm to solve the high-order nonlinear derivative-coupling equation associated with the quartic Galileon model, and implement it in a modified version of the RAMSES N-body code to study the effect of the Galileon field on the large-scale matter clustering. The algorithm is tested for several matter field configurations with different symmetries, and works very well. This enables us to perform the first simulations for a quartic Galileon model which provides a good fit to the cosmic microwave background (CMB) anisotropy, supernovae and baryonic acoustic oscillations (BAO) data. Our result shows that the Vainshtein mechanism in this model is very efficient in suppressing the spatial variations of the scalar field. However, the time variation of the effective Newtonian constant caused by the curvature coupling of the Galileon field cannot be suppressed by the Vainshtein mechanism. This leads to a significant weakening of the strength of gravity in high-density regions at late times, and therefore a weaker matter clustering on small scales. We also find that without the Vainshtein mechanism the model would have behaved in a completely different way, which shows the crucial role played by nonlinearities in modified gravity theories and the importance of performing self-consistent N-body simulations for these theories.

preprint2013arXiv

Systematic simulations of modified gravity: chameleon models

In this work we systematically study the linear and nonlinear structure formation in chameleon theories of modified gravity, using a generic parameterisation which describes a large class of models using only 4 parameters. For this we have modified the N-body simulation code ECOSMOG to perform a total of 65 simulations for different models and parameter values, including the default LCDM. These simulations enable us to explore a significant portion of the parameter space. We have studied the effects of modified gravity on the matter power spectrum and mass function, and found a rich and interesting phenomenology where the difference with the LCDM paradigm cannot be reproduced by a linear analysis even on scales as large as k~0.05h/Mpc, since the latter incorrectly assumes that the modification of gravity depends only on the background matter density. Our results show that the chameleon screening mechanism is significantly more efficient than other mechanisms such as the dilaton and symmetron, especially in high-density regions and at early times, and can serve as a guidance to determine the parts of the chameleon parameter space which are cosmologically interesting and thus merit further studies in the future.

preprint2013arXiv

The Clustering of Galaxies in the SDSS-III DR10 Baryon Oscillation Spectroscopic Survey: No Detectable Colour Dependence of Distance Scale or Growth Rate Measurements

We study the clustering of galaxies, as a function of their colour, from Data Release Ten (DR10) of the SDSS-III Baryon Oscillation Spectroscopic Survey. We select 122,967 galaxies with 0.43 < z < 0.7 into a "Blue" sample and 131,969 into a "Red" sample based on k+e corrected (to z=0.55) r-i colours and i band magnitudes. The samples are chosen to each contain more than 100,000 galaxies, have similar redshift distributions, and maximize the difference in clustering amplitude. The Red sample has a 40% larger bias than the Blue (b_Red/b_Blue = 1.39+-0.04), implying the Red galaxies occupy dark matter halos with an average mass that is 0.5 log Mo greater. Spherically averaged measurements of the correlation function, ξ0, and the power spectrum are used to locate the position of the baryon acoustic oscillation (BAO) feature of both samples. Using ξ0, we obtain distance scales, relative to our reference LCDM cosmology, of 1.010+-0.027 for the Red sample and 1.005+-0.031 for the Blue. After applying reconstruction, these measurements improve to 1.013+/-0.020 for the Red sample and 1.008+-0.026 for the Blue. For each sample, measurements of ξ0 and the second multipole moment, ξ2, of the anisotropic correlation function are used to determine the rate of structure growth, parameterized by fσ8. We find fσ8,Red = 0.511+-0.083, fσ8,Blue = 0.509+/-0.085, and fσ8,Cross = 0.423+-0.061 (from the cross-correlation between the Red and Blue samples). We use the covariance between the bias and growth measurements obtained from each sample and their cross-correlation to produce an optimally-combined measurement of fσ8,comb = 0.443+-0.055. In no instance do we detect significant differences in distance scale or structure growth measurements obtained from the Blue and Red samples.

preprint2012arXiv

Astrophysical Tests of Modified Gravity: A Screening Map of the Nearby Universe

Astrophysical tests of modified modified gravity theories in the nearby universe have been emphasized recently by Hui, Nicolis and Stubbs (2009) and Jain and VanderPlas (2011). A key element of such tests is the screening mechanism whereby general relativity is restored in massive halos or high density environments like the Milky Way. In chameleon theories of gravity, including all f(R) models, field dwarf galaxies may be unscreened and therefore feel an extra force, as opposed to screened galaxies. The first step to study differences between screened and unscreened galaxies is to create a 3D screening map. We use N-body simulations to test and calibrate simple approximations to determine the level of screening in galaxy catalogs. Sources of systematic errors in the screening map due to observational inaccuracies are modeled and their contamination is estimated. We then apply our methods to create a map out to 200 Mpc in the Sloan Digital Sky Survey footprint using data from the Sloan survey and other sources. In two companion papers this map will be used to carry out new tests of gravity using distance indicators and the disks of dwarf galaxies. We also make our screening map publicly available.

preprint2012arXiv

Cosmic Mach Number: A Sensitive Probe for the Growth of Structure

In this Letter, we investigate the potential power of the Cosmic Mach Number (CMN), which is the ratio between the mean velocity and the velocity dispersion of galaxies as a function of cosmic scales, to constrain cosmologies. We first measure the CMN from 5 catalogues of galaxy peculiar velocity surveys at low redshift (0.002<z<0.03), and use them to contrast cosmological models. Overall, current data is consistent with the WMAP7 LCDM model. We find that the CMN is highly sensitive to the growth of structure on scales 0.01<k<0.1 h/Mpc in Fourier space. Therefore, modified gravity models, and models with massive neutrinos, in which the structure growth generally deviates from that in the LCDM model in a scale-dependent way, can be well differentiated from the LCDM model using future CMN data.

preprint2012arXiv

Cosmological tests of General Relativity: a principal component analysis

The next generation of weak lensing surveys will trace the evolution of matter perturbations and gravitational potentials from the matter dominated epoch until today. Along with constraining the dynamics of dark energy, they will probe the relations between matter overdensities, local curvature, and the Newtonian potential. We work with two functions of time and scale to account for any modifications of these relations in the linear regime from those in the LCDM model. We perform a Principal Component Analysis (PCA) to find the eigenmodes and eigenvalues of these functions for surveys like DES and LSST. This paper builds on and significantly extends the PCA analysis of Zhao et al. (2009) in several ways. In particular, we consider the impact of some of the systematic effects expected in weak lensing surveys. We also present the PCA in terms of other choices of the two functions needed to parameterize modified growth on linear scales, and discuss their merits. We analyze the degeneracy between the modified growth functions and other cosmological parameters, paying special attention to the effective equation of state w(z). Finally, we demonstrate the utility of the PCA as an efficient data compression stage which enables one to easily derive constraints on parameters of specific models without recalculating Fisher matrices from scratch.

preprint2012arXiv

Examining the evidence for dynamical dark energy

We apply a new non-parametric Bayesian method for reconstructing the evolution history of the equation-of-state $w$ of dark energy, based on applying a correlated prior for $w(z)$, to a collection of cosmological data. We combine the latest supernova (SNLS 3-year or Union2.1), cosmic microwave background, redshift space distortion and the baryonic acoustic oscillation measurements (including BOSS, WiggleZ and 6dF) and find that the cosmological constant appears consistent with current data, but that a dynamical dark energy model which evolves from $w<-1$ at $z\sim0.25$ to $w > -1$ at higher redshift is mildly favored. Estimates of the Bayesian evidences show little preference between the cosmological constant model and the dynamical model for a range of correlated prior choices. Looking towards future data, we find that the best fit models for current data could be well distinguished from the $Λ$CDM model by observations such as Planck and Euclid-like surveys.

preprint2012arXiv

Interacting Dark Energy -- constraints and degeneracies

In standard cosmologies, dark energy interacts only gravitationally with dark matter. There could be a non-gravitational interaction in the dark sector, leading to changes in the effective DE equation of state, in the redshift dependence of the DM density and in structure formation. We use CMB, BAO and SNIa data to constrain a model where the energy transfer in the dark sector is proportional to the DE density. There are two subclasses, defined by the vanishing of momentum transfer either in the DM or the DE frame. We conduct a Markov-Chain Monte-Carlo analysis to obtain best-fit parameters. The background evolution allows large interaction strengths, and the constraints from CMB anisotropies are weak. The growth of DM density perturbations is much more sensitive to the interaction, and can deviate strongly from the standard case. However, the deviations are degenerate with galaxy bias and thus more difficult to constrain. Interestingly, the ISW signature is suppressed since the non-standard background evolution can compensate for high growth rates. We also discuss the partial degeneracy between interacting DE and modified gravity, and how this can be broken.

preprint2012arXiv

Modified Gravity Spins Up Galactic Halos

We investigate the effect of modified gravity on the specific angular momentum of galactic halos by analyzing the halo catalogs at z=0 from high-resolution N-body simulations for a f(R) gravity model that meets the solar-system constraint. It is shown that the galactic halos in the f(R) gravity model tend to acquire significantly higher specific angular momentum than those in the standard LCDM model. The largest difference in the specific angular momentum distribution between these two models occurs for the case of the isolated galactic halos with mass less than 10^{11}M_sun/h, which are likely least shielded by the chameleon screening mechanism. As the specific angular momentum of galactic halos is rather insensitive to the other cosmological parameters, it can in principle be an independent discriminator of modified gravity. We speculate a possibility of using the relative abundance of the low surface brightness galaxies (LSBGs) as a test of GR given that the formation of the LSBGs befalls in the fast spinning dark halos.

preprint2012arXiv

Redshift space distortions in f(R) gravity

We use large volume N-body simulations to predict the clustering of dark matter in redshift space in f(R) modified gravity cosmologies. This is the first time that the nonlinear matter and velocity fields have been resolved to such a high level of accuracy over a broad range of scales in this class of models. We find significant deviations from the clustering signal in standard gravity, with an enhanced boost in power on large scales and stronger damping on small scales in the f(R) models compared to GR at redshifts z<1. We measure the velocity divergence (P_θθ) and matter (P_δδ) power spectra and find a large deviation in the ratios \sqrt{P_θθ/P_δδ} and P_δθ/P_δδ, between the f(R) models and GR for 0.03<k/(h/Mpc)<0.5. In linear theory these ratios equal the growth rate of structure on large scales. Our results show that the simulated ratios agree with the growth rate for each cosmology (which is scale dependent in the case of modified gravity) only for extremely large scales, k<0.06h/Mpc at z=0. The velocity power spectrum is substantially different in the f(R) models compared to GR, suggesting that this observable is a sensitive probe of modified gravity. We demonstrate how to extract the matter and velocity power spectra from the 2D redshift space power spectrum, P(k,μ), and can recover the nonlinear matter power spectrum to within a few percent for k<0.1h/Mpc. However, the model fails to describe the shape of the 2D power spectrum demonstrating that an improved model is necessary in order to reconstruct the velocity power spectrum accurately. The same model can match the monopole moment to within 3% for GR and 10% for the f(R) cosmology at k<0.2 h/Mpc at z=1. Our results suggest that the extraction of the velocity power spectrum from future galaxy surveys is a promising method to constrain deviations from GR.

preprint2012arXiv

Systematic simulations of modified gravity: symmetron and dilaton models

We study the linear and nonlinear structure formation in the dilaton and symmetron models of modified gravity using a generic parameterisation which describes a large class of scenarios using only a few parameters, such as the coupling between the scalar field and the matter, and the range of the scalar force on very large scales. For this we have modified the N-body simulation code ECOSMOG, which is a variant of RAMSES working in modified gravity scenarios, to perform a set of 110 simulations for different models and parameter values, including the default LCDM. These simulations enable us to explore a large portion of the parameter space. We have studied the effects of modified gravity on the matter power spectrum and mass function, and found a rich and interesting phenomenology where the difference with the LCDM template cannot be reproduced by a linear analysis even on scales as large as k~0.05 h/Mpc. Our results show the full effect of screening on nonlinear structure formation and the associated deviation from LCDM. We also investigate how differences in the force mediated by the scalar field in modified gravity models lead to qualitatively different features for the nonlinear power spectrum and the halo mass function, and how varying the individual model parameters changes these observables. The differences are particularly large in the nonlinear power spectra whose shapes for f(R), dilaton and symmetron models vary greatly, and where the characteristic bump around 1 h/Mpc of f(R) models is preserved for symmetrons, whereas an increase on much smaller scales is particular to symmetrons. No bump is present for dilatons where a flattening of the power spectrum takes place on small scales. These deviations from LCDM and the differences between modified gravity models, such as dilatons and symmetrons, could be tested with future surveys.

preprint2012arXiv

Testing Einstein Gravity with Cosmic Growth and Expansion

We test Einstein gravity using cosmological observations of both expansion and structure growth, including the latest data from supernovae (Union2.1), CMB (WMAP7), weak lensing (CFHTLS) and peculiar velocity of galaxies (WiggleZ). We fit modified gravity parameters of the generalized Poisson equations simultaneously with the effective equation of state for the background evolution, exploring the covariances and model dependence. The results show that general relativity is a good fit to the combined data. Using a Pad{é} approximant form for the gravity deviations accurately captures the time and scale dependence for theories like $f(R)$ and DGP gravity, and weights high and low redshift probes fairly. For current observations, cosmic growth and expansion can be fit simultaneously with little degradation in accuracy, while removing the possibility of bias from holding one aspect fixed.

preprint2012arXiv

The Baryon Oscillation Spectroscopic Survey of SDSS-III

The Baryon Oscillation Spectroscopic Survey (BOSS) is designed to measure the scale of baryon acoustic oscillations (BAO) in the clustering of matter over a larger volume than the combined efforts of all previous spectroscopic surveys of large scale structure. BOSS uses 1.5 million luminous galaxies as faint as i=19.9 over 10,000 square degrees to measure BAO to redshifts z<0.7. Observations of neutral hydrogen in the Lyman alpha forest in more than 150,000 quasar spectra (g<22) will constrain BAO over the redshift range 2.15<z<3.5. Early results from BOSS include the first detection of the large-scale three-dimensional clustering of the Lyman alpha forest and a strong detection from the Data Release 9 data set of the BAO in the clustering of massive galaxies at an effective redshift z = 0.57. We project that BOSS will yield measurements of the angular diameter distance D_A to an accuracy of 1.0% at redshifts z=0.3 and z=0.57 and measurements of H(z) to 1.8% and 1.7% at the same redshifts. Forecasts for Lyman alpha forest constraints predict a measurement of an overall dilation factor that scales the highly degenerate D_A(z) and H^{-1}(z) parameters to an accuracy of 1.9% at z~2.5 when the survey is complete. Here, we provide an overview of the selection of spectroscopic targets, planning of observations, and analysis of data and data quality of BOSS.

preprint2012arXiv

The Clustering of Galaxies in SDSS-III DR9 Baryon Oscillation Spectroscopic Survey: Constraints on Primordial Non-Gaussianity

We analyze the density field of 264,283 galaxies observed by the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) and included in the SDSS data release nine (DR9). In total, the SDSS DR9 BOSS data includes spectroscopic redshifts for over 400,000 galaxies spread over a footprint of more than 3,000 deg^2. We measure the power spectrum of these galaxies with redshifts 0.43 < z < 0.7 in order to constrain the amount of local non-Gaussianity, f_NL,local, in the primordial density field, paying particular attention to the impact of systematic uncertainties. The BOSS galaxy density field is systematically affected by the local stellar density and this influences the ability to accurately measure f_NL,local. In the absence of any correction, we find (erroneously) that the probability that f_NL,local is greater than zero, P(f_NL,local >0), is 99.5%. After quantifying and correcting for the systematic bias and including the added uncertainty, we find -45 < f_NL,local < 195 at 95% confidence, and P(f_NL,local >0) = 91.0%. A more conservative approach assumes that we have only learned the k-dependence of the systematic bias and allows any amplitude for the systematic correction; we find that the systematic effect is not fully degenerate with that of f_NL,local, and we determine that -82 < f_NL,local < 178 (at 95% confidence) and P(f_NL,local >0) = 68%. This analysis demonstrates the importance of accounting for the impact of Galactic foregrounds on f_NL,local measurements. We outline the methods that account for these systematic biases and uncertainties. We expect our methods to yield robust constraints on f_NL,local for both our own and future large-scale-structure investigations.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Analysis of potential systematics

We analyze the density field of galaxies observed by the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) included in the SDSS Data Release Nine (DR9). DR9 includes spectroscopic redshifts for over 400,000 galaxies spread over a footprint of 3,275 deg^2. We identify, characterize, and mitigate the impact of sources of systematic uncertainty on large-scale clustering measurements, both for angular moments of the redshift-space correlation function and the spherically averaged power spectrum, P(k), in order to ensure that robust cosmological constraints will be obtained from these data. A correlation between the projected density of stars and the higher redshift (0.43 < z < 0.7) galaxy sample (the `CMASS' sample) due to imaging systematics imparts a systematic error that is larger than the statistical error of the clustering measurements at scales s > 120h^-1Mpc or k < 0.01hMpc^-1. We find that these errors can be ameliorated by weighting galaxies based on their surface brightness and the local stellar density. We use mock galaxy catalogs that simulate the CMASS selection function to determine that randomly selecting galaxy redshifts in order to simulate the radial selection function of a random sample imparts the least systematic error on correlation function measurements and that this systematic error is negligible for the spherically averaged correlation function. The methods we recommend for the calculation of clustering measurements using the CMASS sample are adopted in companion papers that locate the position of the baryon acoustic oscillation feature (Anderson et al. 2012), constrain cosmological models using the full shape of the correlation function (Sanchez et al. 2012), and measure the rate of structure growth (Reid et al. 2012). (abridged)

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the Data Release 9 Spectroscopic Galaxy Sample

We present measurements of galaxy clustering from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). These use the Data Release 9 (DR9) CMASS sample, which contains 264,283 massive galaxies covering 3275 square degrees with an effective redshift z=0.57 and redshift range 0.43 < z < 0.7. Assuming a concordance Lambda-CDM cosmological model, this sample covers an effective volume of 2.2 Gpc^3, and represents the largest sample of the Universe ever surveyed at this density, n = 3 x 10^-4 h^-3 Mpc^3. We measure the angle-averaged galaxy correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of 5σin both the correlation function and power spectrum. Combining with the SDSS-II Luminous Red Galaxy Sample, the detection significance increases to 6.7σ. Fitting for the position of the acoustic features measures the distance to z=0.57 relative to the sound horizon DV /rs = 13.67 +/- 0.22 at z=0.57. Assuming a fiducial sound horizon of 153.19 Mpc, which matches cosmic microwave background constraints, this corresponds to a distance DV(z=0.57) = 2094 +/- 34 Mpc. At 1.7 per cent, this is the most precise distance constraint ever obtained from a galaxy survey. We place this result alongside previous BAO measurements in a cosmological distance ladder and find excellent agreement with the current supernova measurements. We use these distance measurements to constrain various cosmological models, finding continuing support for a flat Universe with a cosmological constant.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: measurements of the growth of structure and expansion rate at z=0.57 from anisotropic clustering

We analyze the anisotropic clustering of massive galaxies from the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) sample, which consists of 264,283 galaxies in the redshift range 0.43 < z < 0.7 spanning 3,275 square degrees. Both peculiar velocities and errors in the assumed redshift-distance relation ("Alcock-Paczynski effect") generate correlations between clustering amplitude and orientation with respect to the line-of-sight. Together with the sharp baryon acoustic oscillation (BAO) standard ruler, our measurements of the broadband shape of the monopole and quadrupole correlation functions simultaneously constrain the comoving angular diameter distance (2190 +/- 61 Mpc) to z=0.57, the Hubble expansion rate at z=0.57 (92.4 +/- 4.5 km/s/Mpc), and the growth rate of structure at that same redshift (d sigma8/d ln a = 0.43 +/- 0.069). Our analysis provides the best current direct determination of both DA and H in galaxy clustering data using this technique. If we further assume a LCDM expansion history, our growth constraint tightens to d sigma8/d ln a = 0.415 +/- 0.034. In combination with the cosmic microwave background, our measurements of DA, H, and growth all separately require dark energy at z > 0.57, and when combined imply Ω_Λ = 0.74 +/- 0.016, independent of the Universe's evolution at z<0.57. In our companion paper (Samushia et al. prep), we explore further cosmological implications of these observations.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: measuring structure growth using passive galaxies

We explore the benefits of using a passively evolving population of galaxies to measure the evolution of the rate of structure growth between z=0.25 and z=0.65 by combining data from the SDSS-I/II and SDSS-III surveys. The large-scale linear bias of a population of dynamically passive galaxies, which we select from both surveys, is easily modeled. Knowing the bias evolution breaks degeneracies inherent to other methodologies, and decreases the uncertainty in measurements of the rate of structure growth and the normalization of the galaxy power-spectrum by up to a factor of two. If we translate our measurements into a constraint on sigma_8(z=0) assuming a concordance cosmological model and General Relativity (GR), we find that using a bias model improves our uncertainty by a factor of nearly 1.5. Our results are consistent with a flat Lambda Cold Dark Matter model and with GR.

preprint2012arXiv

The Clustering of Galaxies in the SDSS-III DR9 Baryon Oscillation Spectroscopic Survey: Testing Deviations from $Λ$ and General Relativity using anisotropic clustering of galaxies

We use the joint measurement of geometry and growth from anisotropic galaxy clustering in the Baryon Oscillation Spectroscopic Survey Data Release 9 CMASS sample reported by Reid et al. to constrain dark energy properties and possible deviations from the General Relativity. Assuming GR and taking a prior on the linear matter power spectrum at high redshift from the cosmic microwave background (CMB), anisotropic clustering of the CMASS DR9 galaxies alone constrains $Ω_{\rm m} = 0.308 \pm 0.022$ and $100Ω_{\rm k} = 5.9 \pm 4.8$ for $w = -1$, or $w = -0.91 \pm 0.12$ for $Ω_k = 0$. When combined with the full CMB likelihood, the addition of the anisotropic clustering measurements to the spherically-averaged BAO location increases the constraining power on dark energy by a factor of 4 in a flat CDM cosmology with constant dark energy equation of state $w$ (giving $w = -0.87 \pm 0.05$). This impressive gain depends on our measurement of both the growth of structure and Alcock-Paczynski effect, and is not realised when marginalising over the amplitude of redshift space distortions. Combining with both the CMB and Supernovae Type Ia (SNeIa), we find $Ω_{\rm m} = 0.281 \pm 0.014$ and $1000Ω_{\rm k}=-9.2\pm5.0$ for $w = -1$, or $w_0 = -1.13 \pm 0.12$ and $w_{\rm a}=0.65 \pm 0.36$ assuming $Ω_k = 0$. Finally, when a $Λ$CDM background expansion is assumed, the combination of our estimate of the growth rate with previous growth measurements provides tight constraints on the parameters describing possible deviations from GR giving $γ= 0.64 \pm 0.05$. For one parameter extensions of the flat $Λ$CDM model, we find a $\sim 2σ$ preference either for $w > -1$ or slower growth than in GR. However, the data is fully consistent with the concordance model, and the evidence for these additional parameters is weaker than $2σ$.

preprint2012arXiv

The Ninth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Baryon Oscillation Spectroscopic Survey

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

preprint2012arXiv

The nonlinear matter and velocity power spectra in f(R) gravity

We study the matter and velocity divergence power spectra in a f(R) gravity theory and their time evolution measured from several large-volume N-body simulations with varying box sizes and resolution. We find that accurate prediction of the matter power spectrum in f(R) gravity places stronger requirements on the simulation than is the case with LCDM, because of the nonlinear nature of the fifth force. Linear perturbation theory is shown to be a poor approximation for the f(R) models, except when the chameleon effect is very weak. We show that the relative differences from the fiducial LCDM model are much more pronounced in the nonlinear tail of the velocity divergence power spectrum than in the matter power spectrum, which suggests that future surveys which target the collection of peculiar velocity data will open new opportunities to constrain modified gravity theories. A close investigation of the time evolution of the power spectra shows that there is a pattern in the evolution history, which can be explained by the properties of the chameleon-type fifth force in f(R) gravity. Varying the model parameter |f_R0|, which quantifies the strength of the departure from standard gravity, mainly varies the epoch marking the onset of the fifth force, as a result of which the different f(R) models are in different stages of the same evolutionary path at any given time

preprint2011arXiv

Cosmological Measurements with Forthcoming Radio Continuum Surveys

We present forecasts for constraints on cosmological models which can be obtained by forthcoming radio continuum surveys: the wide surveys with the LOw Frequency ARray (LOFAR), Australian Square Kilometre Array Pathfinder (ASKAP) and the Westerbork Observations of the Deep APERTIF Northern sky (WODAN). We use simulated catalogues appropriate to the planned surveys to predict measurements obtained with the source auto-correlation, the cross-correlation between radio sources and CMB maps (the Integrated Sachs-Wolfe effect), the cross-correlation of radio sources with foreground objects due to cosmic magnification, and a joint analysis together with the CMB power spectrum and supernovae. We show that near future radio surveys will bring complementary measurements to other experiments, probing different cosmological volumes, and having different systematics. Our results show that the unprecedented sky coverage of these surveys combined should provide the most significant measurement yet of the Integrated Sachs-Wolfe effect. In addition, we show that using the ISW effect will significantly tighten constraints on modified gravity parameters, while the best measurements of dark energy models will come from galaxy auto-correlation function analyses. Using the combination of EMU and WODAN to provide a full sky survey, it will be possible to measure the dark energy parameters with an uncertainty of \{$σ(w_0) = 0.05$, $σ(w_a) = 0.12$\} and the modified gravity parameters \{$σ(η_0) = 0.10$, $σ(μ_0) = 0.05$\}, assuming Planck CMB+SN(current data) priors. Finally, we show that radio surveys would detect a primordial non-Gaussianity of $f_{\rm NL}$ = 8 at 1-$σ$ and we briefly discuss other promising probes.

preprint2011arXiv

ECOSMOG: An Efficient Code for Simulating Modified Gravity

We introduce a new code, ECOSMOG, to run N-body simulations for a wide class of modified gravity and dynamical dark energy theories. These theories generally have one or more new dynamical degrees of freedom, the dynamics of which are governed by their (usually rather nonlinear) equations of motion. Solving these non-linear equations has been a great challenge in cosmology. Our code is based on the RAMSES code, which solves the Poisson equation on adaptively refined meshes to gain high resolutions in the high-density regions. We have added a solver for the extra degree(s) of freedom and performed numerous tests for the f(R) gravity model as an example to show its reliability. We find that much higher efficiency could be achieved compared with other existing mesh/grid-based codes thanks to two new features of the present code: (1) the efficient parallelisation and (2) the usage of the multigrid relaxation to solve the extra equation(s) on both the regular domain grid and refinements, giving much faster convergence even under much more stringent convergence criteria. This code is designed for performing high-accuracy, high-resolution and large-volume cosmological simulations for modified gravity and general dark energy theories, which can be utilised to test gravity and the dark energy hypothesis using the upcoming and future deep and high-resolution galaxy surveys.

preprint2011arXiv

Fables of reconstruction: controlling bias in the dark energy equation of state

We develop an efficient, non-parametric Bayesian method for reconstructing the time evolution of the dark energy equation of state w(z) from observational data. Of particular importance is the choice of prior, which must be chosen carefully to minimise variance and bias in the reconstruction. Using a principal component analysis, we show how a correlated prior can be used to create a smooth reconstruction and also avoid bias in the mean behaviour of w(z). We test our method using Wiener reconstructions based on Fisher matrix projections, and also against more realistic MCMC analyses of simulated data sets for Planck and a future space-based dark energy mission. While the accuracy of our reconstruction depends on the smoothness of the assumed w(z), the relative error for typical dark energy models is <10% out to redshift z=1.5.

preprint2011arXiv

Halos and Voids in f(R) Gravity

In this paper, we study the distribution of dark matter halos and voids using high resolution simulations in f(R) gravity models with the chameleon mechanism to screen the fifth force in dense environment. For dark matter halos, we show that the semi-analytic thin shell condition, with a suitably-defined environment, provides a good approximation to describe the mass and environmental dependence of the screening of the fifth force in halos. Due to stronger gravity, there are far more massive halos and large voids in f(R) models compared with the ΛCDM model. The numbers of voids with an effective radius of 15Mpc/h are twice and four times as many as those in ΛCDM for f(R) models with |f_{R0}|=1e-5 and 1e-4 respectively. This provides a new means to test the models using the upcoming observational data. We also find that halos inside voids are all unscreened in our simulations, which are ideal objects for the gravity test.

preprint2011arXiv

How to optimally parametrize deviations from General Relativity in the evolution of cosmological perturbations

The next generation of weak lensing surveys will trace the growth of large scale perturbations through a sequence of epochs, offering an opportunity to test General Relativity (GR) on cosmological scales. We review in detail the parametrization used in MGCAMB to describe the modified growth expected in alternative theories of gravity and generalized dark energy models. We highlight its advantages and examine several theoretical aspects. In particular, we show that the same set of equations can be consistently used on super-horizon and sub-horizon linear scales. We also emphasize the sensitivity of data to scale-dependent features in the growth pattern, and propose using Principal Component Analysis to converge on a practical set of parameters which is most likely to detect departures from GR. The connection with other parametrizations is also discussed.

preprint2011arXiv

Testing General Relativity using the Environmental Dependence of Dark Matter Halos

In this Letter, we investigate the environmental dependence of dark matter halos in theories that attempt to explain the accelerated expansion of the Universe by modifying general relativity (GR). Using high-resolution N-body simulations in f(R) gravity models which recover GR in dense environments by virtue of the chameleon mechanism, we find a strong environmentally-dependent difference between the lensing mass and dynamical mass estimates of dark matter halos. This environmental dependence of the halo properties can be used as a smoking gun to test GR observationally.

preprint2011arXiv

Testing gravity with CAMB and CosmoMC

We introduce a patch to the commonly used public codes CAMB and CosmoMC that allows the user to implement a general modification of the equations describing the growth of cosmological perturbations, while preserving the covariant conservation of the energy-momentum. This patch replaces the previously publicly released code MGCAMB, while also extending it in several ways. The new version removes the limitation of late-time-only modifications to the perturbed Einstein equations, and includes several parametrization introduced in the literature. To demonstrate the use of the patch, we obtain joint constraints on the neutrino mass and parameters of a scalar-tensor gravity model from CMB, SNe and ISW data as measured from the correlation of CMB with large scale structure.

preprint2010arXiv

A Semi-analytic Ray-tracing Algorithm for Weak Lensing

We propose a new ray-tracing algorithm to measure the weak lensing shear and convergence fields directly from N-body simulations. We calculate the deflection of the light rays lensed by the 3-D mass density field or gravitational potential along the line of sight on a grid-by-grid basis, rather than using the projected 2-D lens planes. Our algorithm uses simple analytic formulae instead of numerical integrations in the computation of the projected density field along the line of sight, and so is computationally efficient, accurate and straightforward to implement. This will prove valuable in the interpretation of data from the next generation of surveys that will image many thousands of square degrees of sky.

preprint2010arXiv

Complementarity of Weak Lensing and Peculiar Velocity Measurements in Testing General Relativity

We explore the complementarity of weak lensing and galaxy peculiar velocity measurements to better constrain modifications to General Relativity. We find no evidence for deviations from GR on cosmological scales from a combination of peculiar velocity measurements (for Luminous Red Galaxies in the Sloan Digital Sky Survey) with weak lensing measurements (from the CFHT Legacy Survey). We provide a Fisher error forecast for a Euclid-like space-based survey including both lensing and peculiar velocity measurements, and show that the expected constraints on modified gravity will be at least an order of magnitude better than with present data, i.e. we will obtain 5% errors on the modified gravity parametrization described here. We also present a model--independent method for constraining modified gravity parameters using tomographic peculiar velocity information, and apply this methodology to the present dataset.

preprint2010arXiv

Detecting Features in the Dark Energy Equation of State: A Wavelet Approach

We study the utility of wavelets for detecting the redshift evolution of the dark energy equation of state w(z) from the combination of supernovae, CMB and BAO data. We show that local features in w, such as bumps, can be detected efficiently using wavelets. To demonstrate, we first generate a mock supernovae (SNe) data sample for a SNAP-like survey with a bump feature in w(z) hidden in, then successfully discover it by performing a blind wavelet analysis. We also apply our method to analyze the recently released "Constitution" SNe data, combined with WMAP and BAO from SDSS, and find weak hints of dark energy dynamics. Namely, we find that models with w(z) < -1 for 0.2 < z < 0.5, and w(z)> -1 for 0.5 < z <1, are mildly favored at 95% confidence level. This is in good agreement with several recent studies using other methods, such as redshift binning with principal component analysis (PCA) (e.g. Zhao and Zhang, arXiv:0908.1568)

preprint2010arXiv

N-body Simulations for f(R) Gravity using a Self-adaptive Particle-Mesh Code

We perform high resolution N-body simulations for f(R) gravity based on a self-adaptive particle- mesh code MLAPM. The Chameleon mechanism that recovers General Relativity on small scales is fully taken into account by self-consistently solving the non-linear equation for the scalar field. We independently confirm the previous simulation results, including the matter power spectrum, halo mass function and density profiles, obtained by Oyaizu et al. (Phys.Rev.D 78, 123524, 2008) and Schmidt et al. (Phys.Rev.D 79, 083518, 2009), and extend the resolution up to k~20 h/Mpc for the measurement of the matter power spectrum. Based on our simulation results, we discuss how the Chameleon mechanism affects the clustering of dark matter and halos on full non-linear scales.

preprint2010arXiv

Probing Dark Energy Dynamics from Current and Future Cosmological Observations

We report the constraints on the dark energy equation-of-state w(z) using the latest 'Constitution' SNe sample combined with the WMAP5 and SDSS data. Based on the localized principal component analysis and the model selection criteria, we find that the LCDM model is generally consistent with the current data, yet there exists weak hint of the possible dynamics of dark energy. In particular, a model predicting w(z)<-1 at z\in[0.25,0.5) and w(z)>-1 at z\in[0.5,0.75), which means that w(z) crosses -1 in the range of z\in[0.25,0.75), is mildly favored at 95% confidence level. Given the best fit model for current data as a fiducial model, we make future forecast from the joint data sets of JDEM, Planck and LSST, and we find that the future surveys can reduce the error bars on the w bins by roughly a factor of 10 for a 5-w-bin model.

preprint2010arXiv

Probing Dark Energy with the Kunlun Dark Universe Survey Telescope

Dark energy is an important science driver of many upcoming large-scale surveys. With small, stable seeing and low thermal infrared background, Dome A, Antarctica, offers a unique opportunity for shedding light on fundamental questions about the universe. We show that a deep, high-resolution imaging survey of 10,000 square degrees in \emph{ugrizyJH} bands can provide competitive constraints on dark energy equation of state parameters using type Ia supernovae, baryon acoustic oscillations, and weak lensing techniques. Such a survey may be partially achieved with a coordinated effort of the Kunlun Dark Universe Survey Telescope (KDUST) in \emph{yJH} bands over 5000--10,000 deg$^2$ and the Large Synoptic Survey Telescope in \emph{ugrizy} bands over the same area. Moreover, the joint survey can take advantage of the high-resolution imaging at Dome A to further tighten the constraints on dark energy and to measure dark matter properties with strong lensing as well as galaxy--galaxy weak lensing.

preprint2010arXiv

Probing modifications of General Relativity using current cosmological observations

We test General Relativity (GR) using current cosmological data: the cosmic microwave background (CMB) from WMAP5 (Komatsu et al. 2009), the integrated Sachs-Wolfe (ISW) effect from the cross-correlation of the CMB with six galaxy catalogs (Giannantonio et al. 2008), a compilation of supernovae Type Ia (SNe) including the latest SDSS SNe (Kessler et al. 2009), and part of the weak lensing (WL) data from CFHTLS (Fu et al. 2008, Kilbinger et al. 2009) that probe linear and mildly non-linear scales. We first test a model where the effective Newton's constant, mu, and the ratio of the two gravitational potentials, eta, transit from the GR value to another constant at late times; in this case, we find that standard GR is fully consistent with the combined data. The strongest constraint comes from the ISW effect which would arise from this gravitational transition; the observed ISW signal imposes a tight constraint on a combination of mu and eta that characterizes the lensing potential. Next, we consider four pixels in time and space for each function mu and eta, and perform a Principal Component Analysis (PCA) finding that seven of the resulting eight eigenmodes are consistent with GR within the errors. Only one eigenmode shows a 2-sigma deviation from the GR prediction, which is likely to be due to a systematic effect. However, the detection of such a deviation demonstrates the power of our time- and scale-dependent PCA methodology when combining observations of structure formation and expansion history to test GR.

preprint2009arXiv

Cosmological Tests of General Relativity with Future Tomographic Surveys

Future weak lensing surveys will map the evolution of matter perturbations and gravitational potentials, yielding a new test of general relativity on cosmic scales. They will probe the relations between matter overdensities, local curvature, and the Newtonian potential. These relations can be modified in alternative gravity theories or by the effects of massive neutrinos or exotic dark energy fluids. We introduce two functions of time and scale which account for any such modifications in the linear regime. We use a principal component analysis to find the eigenmodes of these functions that cosmological data will constrain. The number of constrained modes gives a model-independent forecast of how many parameters describing deviations from general relativity could be constrained, along with $w(z)$. The modes' scale and time dependence tell us which theoretical models will be better tested.

preprint2009arXiv

Investigating dark energy experiments with principal components

We use a principal component approach to contrast different kinds of probes of dark energy, and to emphasize how an array of probes can work together to constrain an arbitrary equation of state history w(z). We pay particular attention to the role of the priors in assessing the information content of experiments and propose using an explicit prior on the degree of smoothness of w(z) that is independent of the binning scheme. We also show how a figure of merit based on the mean squared error probes the number of new modes constrained by a data set, and use it to examine how informative various experiments will be in constraining the evolution of dark energy.

preprint2008arXiv

On using the WMAP distance priors in constraining the time evolving equation of state of dark energy

Recently, the WMAP group has published their five-year data and considered the constraints on the time evolving equation of state of dark energy for the first time from the WMAP distance information. In this paper, we study the effectiveness of the usage of these distance information and find that these compressed CMB information can give similar constraints on dark energy parameters compared with the full CMB power spectrum if dark energy perturbations are included, however, once incorrectly neglecting the dark energy perturbations, the difference of the results are sizable.

preprint2008arXiv

Overcoming the Circular Problem for γ-ray Bursts in Cosmological Global Fitting Analysis

Due to the lack of low redshift long Gamma-Ray Bursts (GRBs), the circular problem has been a severe obstacle for using GRBs as cosmological candles. In this paper, we present a new method to deal with such a problem in MCMC global fitting analysis. Assuming that a certain type of correlations between different observables exists in a subsample of GRBs, for the parameters involved in the correlation relation, we treat them as free parameters and determine them simultaneously with cosmological parameters through MCMC analysis on GRB data together with other observational data. Then the circular problem is naturally eliminated in this procedure. We take the Ghirlanda relation as an example while keeping in mind the debate about its physical validity. Together with SNe Ia, WMAP and SDSS data, we include 27 GRBs with the reported Ghirlanda relation in our study, and perform MCMC global fitting. We consider the $Λ$CDM model and dynamical dark energy models. In each case, in addition to the constraints on the relevant cosmological parameters, we obtain the best fit values as well as the distributions of the correlation parameters $A$ and $C$. We find that the observational data sets other than GRBs can affect $A$ and $C$ considerably through their degeneracies with the cosmological parameters. The results on $A$ and $C$ for different cosmological models are in well agreement within $1σ$ range. The best fit value of $A$ in all models being analyzed is $A\sim 1.53$ with $σ\sim 0.08$. For $C$, we have the best value in the range of $0.94-0.98$ with $σ\sim 0.1$. It is also noted that the distributions of $A$ and $C$ are generally broader than the priors used in many studies in literature. (Abriged)

preprint2008arXiv

Probing for the Cosmological Parameters with PLANCK Measurement

We investigate the constraints on cosmological parameters especially for EoS of dark energy, inflationary parameters, neutrino mass and curvature of universe using simulated Planck data. Firstly we determine cosmological parameters with current observations including ESSENCE, WMAP3, Boomerang-2K2, CBI, VSA, ACBAR, SDSS LRG and 2dFGRS, and take best-fit model as the fiducial model in simulations. In simulations we pay attention to the effects of dynamical dark energy in determination of cosmological parameters. We add simulated SNAP data to do all the simulations. Using present data, we find Quintom dark energy model is mildly favored while \LambdaCDM remains a good fit. In the framework of dynamical dark energy, the constraints on inflationary parameters, m_ν and Ω_{K} become weak compared with the constraints in \LambdaCDM. Intriguingly, we find that the inflationary models with a "blue" tilt, which are excluded about 2σin \LambdaCDM model, are well within 2σregion with the presence of the dynamics of dark energy. The upper limits of neutrino mass are weakened by a factor of 2 (95% C.L.), say, m_ν<1.59 eV and m_ν<1.53 eV for two forms of parametrization of the equation of state of dark energy. The flat universe is a good fit to the current data, namely, |Ω_{K}|<0.03 (95% C.L.). With the simulated Planck and SNAP data, dynamical dark energy and \LambdaCDM might be distinguished at 4σ. And uncertainties of inflationary parameters, m_ν and Ω_{K} can be reduced obviously. We also constrain the rotation angle Δα, denoting possible cosmological CPT violation, with simulated Planck and CMBpol data and find that our results are much more stringent than current constraint and will verify cosmological CPT symmetry with a higher precision. (Abridged)

preprint2008arXiv

Testing CPT Symmetry with CMB Measurements: Update after WMAP5

In this Letter we update our previous result on the test of CPT symmetry with Cosmic Microwave Background (CMB) measurements. A CPT violating interaction in the photon sector ${\cal L}_{cs} \sim p_μA_ν\tilde F^{μν}$ gives rise to a rotation of the polarization vectors of the propagating CMB photons. Recently the WMAP group used the newly released polarization data of WMAP5 to measure this rotation angle $Δα$ and obtained $Δα=-1.7\pm2.1$ deg ($1σ$). However, in their analysis the BOOMERanG 2003 data is not included. Here we revisit this issue by combining the full data of WMAP5 and BOOMERanG 2003 angular power spectra for the measurement of this rotation angle $Δα$ and find that $Δα=-2.6\pm1.9$ deg at a 68% confidence level.

preprint2007arXiv

Probing for variation of neutrino mass with current observations

With the latest astronomical data including Cosmic Microwave Background (WMAP three year, CBI, ACBAR, VSA), Type Ia Supernova ("gold sample"), Galaxy Clustering (SDSS 3-D matter power, Lyman-$α$ forest and Baryon Acoustic Oscillating (BAO)), we make a global fitting to constrain the mass varying neutrinos. We find that the parameter $δ$, denoting time evolving of neutrino mass, is weakly constrained and the neutrino mass limit today can be relaxed at least by a factor of two. Adding data of $0\nu2β$ decay of Heidelberg-Moscow experiment to our analysis, we find that $δ$ can be constrained tightly and mass varying neutrinos are favored at about 99.7% confidence level.