Search

Uncertain feedback processes in galaxies affect the distribution of matter, currently limiting the power of weak lensing surveys. If we can identify cosmological statistics that are robust against these uncertainties, or constrain these effects by other means, then we can enhance the power of current and upcoming observations from weak lensing surveys such as DES, Euclid, the Rubin Observatory, and the Roman Space Telescope. In this work, we investigate the potential of the electron density auto-power spectrum as a robust probe of cosmology and baryonic feedback. We use a suite of (magneto-)hydrodynamic simulations from the CAMELS project and perform an idealized analysis to forecast statistical uncertainties on a limited set of cosmological and physically-motivated astrophysical parameters. We find that the electron number density auto-correlation, measurable through either kinematic Sunyaev-Zel'dovich observations or through Fast Radio Burst dispersion measures, provides tight constraints on $\Omega_{m}$ and the mean baryon fraction in intermediate-mass halos, $\bar{f}_{\mathrm{bar}}$. By obtaining an empirical measure for the associated systematic uncertainties, we find these constraints to be largely robust to differences in baryonic feedback models implemented in hydrodynamic simulations. We further discuss the main caveats associated with our analysis, and point out possible directions for future work., Comment: 31 pages, 10 figures, to be submitted to JCAP

The combination of different cosmological probes offers stringent tests of the $\Lambda$CDM model and enhanced control of systematics. For this purpose, we present an extension of the lightcone generator UFalcon first introduced in Sgier et al. 2019 (arXiv:1801.05745), enabling the simulation of a self-consistent set of maps for different cosmological probes. Each realization is generated from the same underlying simulated density field, and contains full-sky maps of different probes, namely weak lensing shear, galaxy overdensity including RSD, CMB lensing, and CMB temperature anisotropies from the ISW effect. The lightcone generation performed by UFalcon is parallelized and based on the replication of a large periodic volume simulated with the GPU-accelerated $N$-Body code PkdGrav3. The post-processing to construct the lightcones requires only a runtime of about 1 walltime-hour corresponding to about 100 CPU-hours. We use a randomization procedure to increase the number of quasi-independent full-sky UFalcon map-realizations, which enables us to compute an accurate multi-probe covariance matrix. Using this framework, we forecast cosmological parameter constraints by performing a multi-probe likelihood analysis for a combination of simulated future stage-IV-like surveys. We find that the inclusion of the cross-correlations between the probes significantly increases the information gain in the parameter constraints. We also find that the use of a non-Gaussian covariance matrix is increasingly important, as more probes and cross-correlation power spectra are included. A version of the UFalcon package currently including weak gravitational lensing is publicly available., Comment: 49 pages, 24 pictures, The UFalcon weak lensing package is available here: $\href{https://cosmology.ethz.ch/research/software-lab/UFalcon.html}{https://cosmology.ethz.ch/research/software-lab/UFalcon.html}$

Current and upcoming cosmological experiments open a new era of precision cosmology, thus demanding accurate theoretical predictions for cosmological observables. Because of the complexity of the codes delivering such predictions, reaching a high level of numerical accuracy is challenging. Among the codes already fulfilling this task, PyCosmo is a Python-based framework providing solutions to the Einstein–Boltzmann equations and accurate predictions for cosmological observables. We present the first public release of the code, which is valid in ΛCDM cosmology. The novel aspect of this version is that the user can work within a Python framework, either locally or through an online platform, called PyCosmo Hub. In this work we first describe how the observables are implemented. Then, we check the accuracy of the theoretical predictions for background quantities, power spectra and Limber and beyond-Limber angular power spectra by comparison with other codes: the Core Cosmology Library (CCL), CLASS, HMCode and iCosmo. In our analysis we quantify the agreement of PyCosmo with the other codes, for a range of cosmological models, monitored through a series of unit tests. PyCosmo, conceived as a multi-purpose cosmology calculation tool in Python, is designed to be interactive and user-friendly. The PyCosmo Hub is accessible from this link: https://cosmology.ethz.ch/research/software-lab/PyCosmo.html. On this platform the users can perform their own computations using Jupyter Notebooks without the need of installing any software, access to the results presented in this work and benefit from tutorial notebooks illustrating the usage of the code. The link above also redirects to the code release and documentation., Astronomy and Computing, 36, ISSN:2213-1337

46 pags., 12 figs., 7 tabs., n order to investigate the origin of the ongoing tension between the amplitude of matter fluctuations measured by weak lensing experiments at low redshifts and the value inferred from the cosmic microwave background anisotropies, we reconstruct the evolution of this amplitude from z ~ 2 using existing large-scale structure data. To do so, we decouple the linear growth of density inhomogeneities from the background expansion, and constrain its redshift dependence making use of a combination of 6 different data sets, including cosmic shear, galaxy clustering and CMB lensing. We analyze these data under a consistent harmonic-space angular power spectrum-based pipeline. We show that current data constrain the amplitude of fluctuations mostly in the range 0.2 < z < 0.7, where it is lower than predicted by . This difference is mostly driven by current cosmic shear data, although the growth histories reconstructed from different data combinations are consistent with each other, and we find no evidence of systematic deviations in any particular experiment. In spite of the tension with , the data are well-described by the model, albeit with a lower value of S8 ¿ ¿8(¿m/0.3)0.5. As part of our analysis, we find constraints on this parameter of S8 = 0.7781 ± 0.0094 (68% confidence level), reaching almost percent-level errors comparable with CMB measurements, and 3.4¿ away from the value found by ., AN is supported by NSF grant AST-181497. CGG, EB, EMM and PGF are supported by European Research Council Grant No: 693024 and the Beecroft Trust. CGG was also supported by the grant PGC2018-095157-B-I00 from Ministry of Science, Innovation and Universities of Spain and by the Spanish grant, partially funded by the ESF, BES-2016- 077038. PRL is supported by the grant PGC2018-095157-B-I00 from Ministry of Science, Innovation and Universities of Spain. DA is supported by the Science and Technology Facilities Council through an Ernest Rutherford Fellowship, grant reference ST/P004474. We also made extensive use of computational resources at the University of Oxford Department of Physics, funded by the John Fell Oxford University Press Research Fund.

We present a new method for joint cosmological parameter inference and cluster mass calibration from a combination of weak lensing measurements and the abundance of thermal Sunyaev-Zel'dovich (tSZ) selected galaxy clusters. We combine cluster counts with the spherical harmonic cosmic shear power spectrum and the cross-correlation between cluster overdensity and cosmic shear. These correlations constrain the cluster mass-observable relation. We model the observables using a halo model framework, including their full non-Gaussian covariance. Forecasting constraints on cosmological and mass calibration parameters for a combination of LSST cosmic shear and Simons Observatory tSZ cluster counts, we find competitive constraints for cluster cosmology, with a factor of 2 improvement in the dark energy figure of merit compared to LSST cosmic shear alone. We find most of the mass calibration information will be in the large and intermediate scales of the cross-correlation between cluster overdensity and cosmic shear. Finally, we find broadly comparable constraints to traditional analyses based on calibrating masses using stacked cluster lensing measurements, with the benefit of consistently accounting for the correlations with cosmic shear.

DES Collaboration: et al., Weak lensing by large-scale structure is a powerful probe of cosmology and of the dark universe. This cosmic shear technique relies on the accurate measurement of the shapes and redshifts of background galaxies and requires precise control of systematic errors. Monte Carlo control loops (MCCL) is a forward modeling method designed to tackle this problem. It relies on the ultra fast image generator (UFig) to produce simulated images tuned to match the target data statistically, followed by calibrations and tolerance loops. We present the first end-to-end application of this method, on the Dark Energy Survey (DES) Year 1 wide field imaging data. We simultaneously measure the shear power spectrum Cℓ and the redshift distribution n(z) of the background galaxy sample. The method includes maps of the systematic sources, point spread function (PSF), an approximate Bayesian computation (ABC) inference of the simulation model parameters, a shear calibration scheme, and a fast method to estimate the covariance matrix. We find a close statistical agreement between the simulations and the DES Y1 data using an array of diagnostics. In a nontomographic setting, we derive a set of Cℓ and n(z) curves that encode the cosmic shear measurement, as well as the systematic uncertainty. Following a blinding scheme, we measure the combination of Ωm, σ8, and intrinsic alignment amplitude AIA, defined as S8DIA=σ8(Ωm/0.3)0.5DIA, where DIA=1−0.11(AIA−1). We find S8DIA=0.895+0.054−0.039, where systematics are at the level of roughly 60% of the statistical errors. We discuss these results in the context of earlier cosmic shear analyses of the DES Y1 data. Our findings indicate that this method and its fast runtime offer good prospects for cosmic shear measurements with future wide-field surveys., Funding for the DES Projects has been provided by the U.S. Department of Energy, the U.S. National Science Foundation, the Ministry of Science and Education of Spain, the Science and Technology Facilities Council of the United Kingdom, the Higher Education Funding Council for England, the National Center for Supercomputing Applications at the University of Illinois at Urbana-Champaign, the Kavli Institute of Cosmological Physics at the University of Chicago, the Center for Cosmology and Astro-Particle Physics at the Ohio State University, the Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University, Financiadora de Estudos e Projetos, Fundação Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro, Conselho Nacional de Desenvolvimento Científico e Tecnológico and the Ministerio da Ciência, Tecnologia e Inovação, the Deutsche Forschungsgemeinschaft and the Collaborating Institutions in the Dark Energy Survey. The Collaborating Institutions are Argonne National Laboratory, the University of California at Santa Cruz, the University of Cambridge, Centro de Investigaciones Energeticas, Medioambientales y Tecnológicas-Madrid, the University of Chicago, University College London, the DESBrazil Consortium, the University of Edinburgh, the Eidgenössische Technische Hochschule (ETH) Zürich, Fermi National Accelerator Laboratory, the University of Illinois at Urbana-Champaign, the Institut de Ciencies de l’Espai (IEEC/CSIC), the Institut de Física d’Altes Energies, Lawrence Berkeley National Laboratory, the LudwigMaximilians Universität München and the associated Excellence Cluster Universe, the University of Michigan, the National Optical Astronomy Observatory, the University of Nottingham, The Ohio State University, the University of Pennsylvania, the University of Portsmouth, SLAC National Accelerator Laboratory, Stanford University, the University of Sussex, Texas A&M University, and the OzDES Membership Consortium. Based in part on observations at Cerro Tololo Inter-American Observatory, National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. The DES data management system is supported by the National Science Foundation under Grants No, AST-1138766 and No. AST-1536171. The DES participants from Spanish institutions are partially supported by MINECO under Grants No. AYA2015-71825, No. ESP2015-66861, No. FPA2015-68048, No. SEV-2016-0588, No. SEV2016-0597, and No. MDM-2015-0509, some of which include ERDF funds from the European Union. I. F. A. E. is partially funded by the CERCA program of the Generalitat de Catalunya. Research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Program (FP7/2007-2013) including ERC Grant agreements No. 240672, No. 291329, and No. 306478. We acknowledge support from the Brazilian Instituto Nacional de Ciência e Tecnologia (INCT) e-Universe (CNPq grant No. 465376/2014-2).

We present a new method to marginalize over uncertainties in redshift distributions, $N(z)$, within tomographic cosmological analyses applicable to current and upcoming photometric galaxy surveys. We allow for arbitrary deviations from the best-guess $N(z)$ governed by a general covariance matrix describing the uncertainty in our knowledge of redshift distributions. In principle, this is marginalization over hundreds or thousands of new parameters describing potential deviations as a function of redshift and tomographic bin. However, by linearly expanding the theory predictions around a fiducial model, this marginalization can be performed analytically, resulting in a modified data covariance matrix that effectively downweights the modes of the data vector that are more sensitive to redshift distribution variations. We showcase this method by applying it to the galaxy clustering measurements from the Hyper Suprime-Cam first data release. We illustrate how to marginalize over sample-variance of the calibration sample and a large general systematic uncertainty in photometric estimation methods, and explore the impact of priors imposing smoothness in the redshift distributions., Comment: 26 pages, 7 figures

We present the temperature and polarization angular power spectra of the CMB measured by the Atacama Cosmology Telescope (ACT) from 5400 deg$^2$ of the 2013-2016 survey, which covers $>$15000 deg$^2$ at 98 and 150 GHz. For this analysis we adopt a blinding strategy to help avoid confirmation bias and, related to this, show numerous checks for systematic error done before unblinding. Using the likelihood for the cosmological analysis we constrain secondary sources of anisotropy and foreground emission, and derive a "CMB-only" spectrum that extends to $\ell=4000$. At large angular scales, foreground emission at 150 GHz is $\sim$1% of TT and EE within our selected regions and consistent with that found by Planck. Using the same likelihood, we obtain the cosmological parameters for $\Lambda$CDM for the ACT data alone with a prior on the optical depth of $\tau=0.065\pm0.015$. $\Lambda$CDM is a good fit. The best-fit model has a reduced $\chi^2$ of 1.07 (PTE=0.07) with $H_0=67.9\pm1.5$ km/s/Mpc. We show that the lensing BB signal is consistent with $\Lambda$CDM and limit the celestial EB polarization angle to $\psi_P =-0.07^{\circ}\pm0.09^{\circ}$. We directly cross correlate ACT with Planck and observe generally good agreement but with some discrepancies in TE. All data on which this analysis is based will be publicly released., Comment: 44 pages, 27 figures, products available on the NASA LAMBDA website, version accepted for publication in JCAP

We present a re-analysis of cosmic shear and galaxy clustering from first-year Dark Energy Survey data (DES Y1), making use of a Hybrid Effective Field Theory (HEFT) approach to model the galaxy-matter relation on weakly non-linear scales, initially proposed in Modi et al. (2020) (arXiv:1910.07097). This allows us to explore the enhancement in cosmological constraining power enabled by extending the galaxy clustering scale range typically used in projected large-scale structure analyses. Our analysis is based on a recomputed harmonic-space data vector and covariance matrix, carefully accounting for all sources of mode-coupling, non-Gaussianity and shot noise, which allows us to provide robust goodness-of-fit measures. We use the \textsc{AbacusSummit} suite of simulations to build an emulator for the HEFT model predictions. We find that this model can explain the galaxy clustering and shear data up to wavenumbers $k_{\rm max}\sim 0.6\, {\rm Mpc}^{-1}$. We constrain $(S_8,\Omega_m) = (0.786\pm 0.020,0.273^{+0.030}_{-0.036})$ at the fiducial $k_{\rm max}\sim 0.3\, {\rm Mpc}^{-1}$, improving to $(S_8,\Omega_m) = (0.786^{+0.015}_{-0.018},0.266^{+0.024}_{-0.027})$ at $k_{\rm max}\sim 0.5\, {\rm Mpc}^{-1}$. This represents a $\sim10\%$ and $\sim35\%$ improvement on the constraints derived respectively on both parameters using a linear bias relation on a reduced scale range ($k_{\rm max}\lesssim0.15\,{\rm Mpc}^{-1}$), in spite of the 15 additional parameters involved in the HEFT model. We investigate whether HEFT can be used to constrain the Hubble parameter and find $H_0= 70.7_{-3.5}^{+3.0}\,{\rm km}\,s^{-1}\,{\rm Mpc}^{-1}$. Our constraints are investigative and subject to certain caveats discussed in the text., Comment: Prepared for JCAP; 29 pages, 7 figures, 3 tables

We present new cosmological parameter constraints for general Horndeski scalar-tensor theories, using CMB, redshift space distortion, matter power spectrum and BAO measurements from the Planck, SDSS/BOSS and 6dF surveys. We focus on theories with cosmological gravitational waves propagating at the speed of light, $c_{\rm GW} = c$, implementing and discussing several previously unaccounted for aspects in the constraint derivation for such theories, that qualitatively affect the resulting constraints. In order to ensure our conclusions are robust, we compare results for three different parametrisations of the free functions in Horndeski scalar-tensor theories, identifying several parametrisation-independent features of the constraints. We also consider models, where $c_{\rm GW} \neq c$ in cosmological settings (still allowed after GW170817 for frequency-dependent $c_{\rm GW}$) and show how this affects cosmological parameter constraints., Comment: 30 pages, 9 figures, 3 tables

We analyze the clustering of galaxies in the first public data release of the HSC Subaru Strategic Program. Despite the relatively small footprints of the observed fields, the data are an excellent proxy for the deep photometric datasets that will be acquired by LSST, and are therefore an ideal test bed for the analysis methods being implemented by the LSST DESC. We select a magnitude limited sample with $i, Comment: 65 pages, 30 figures, 5 tables, to be submitted to JCAP

Correlations between tracers of the matter density field and gravitational lensing are sensitive to the evolution of the matter power spectrum and the expansion rate across cosmic time. Appropriately defined ratios of such correlation functions, on the other hand, depend only on the angular diameter distances to the tracer objects and to the gravitational lensing source planes. Because of their simple cosmological dependence, such ratios can exploit available signal-to-noise down to small angular scales, even where directly modeling the correlation functions is difficult. We present a measurement of lensing ratios using galaxy position and lensing data from the Dark Energy Survey, and CMB lensing data from the South Pole Telescope and Planck, obtaining the highest precision lensing ratio measurements to date. Relative to the concordance $\Lambda$CDM model, we find a best fit lensing ratio amplitude of $A = 1.1 \pm 0.1$. We use the ratio measurements to generate cosmological constraints, focusing on the curvature parameter. We demonstrate that photometrically selected galaxies can be used to measure lensing ratios, and argue that future lensing ratio measurements with data from a combination of LSST and Stage-4 CMB experiments can be used to place interesting cosmological constraints, even after considering the systematic uncertainties associated with photometric redshift and galaxy shear estimation., Comment: 17 pages, 11 figures

Cosmic shear is one of the most powerful probes of Dark Energy, targeted by several current and future galaxy surveys. Lensing shear, however, is only sampled at the positions of galaxies with measured shapes in the catalog, making its associated sky window function one of the most complicated amongst all projected cosmological probes of inhomogeneities, as well as giving rise to inhomogeneous noise. Partly for this reason, cosmic shear analyses have been mostly carried out in real-space, making use of correlation functions, as opposed to Fourier-space power spectra. Since the use of power spectra can yield complementary information and has numerical advantages over real-space pipelines, it is important to develop a complete formalism describing the standard unbiased power spectrum estimators as well as their associated uncertainties. Building on previous work, this paper contains a study of the main complications associated with estimating and interpreting shear power spectra, and presents fast and accurate methods to estimate two key quantities needed for their practical usage: the noise bias and the Gaussian covariance matrix, fully accounting for survey geometry, with some of these results also applicable to other cosmological probes. We demonstrate the performance of these methods by applying them to the latest public data releases of the Hyper Suprime-Cam and the Dark Energy Survey collaborations, quantifying the presence of systematics in our measurements and the validity of the covariance matrix estimate. We make the resulting power spectra, covariance matrices, null tests and all associated data necessary for a full cosmological analysis publicly available., 42 pages, 16 figures, 4 tables, to be submitted to JCAP

The nature of dark matter, dark energy and large-scale gravity pose some of the most pressing questions in cosmology today. These fundamental questions require highly precise measurements, and a number of wide-field spectroscopic survey instruments are being designed to meet this requirement. A key component in these experiments is the development of a simulation tool to forecast science performance, define requirement flow-downs, optimize implementation, demonstrate feasibility, and prepare for exploitation. We present SPOKES (SPectrOscopic KEn Simulation), an end-to-end simulation facility for spectroscopic cosmological surveys designed to address this challenge. SPOKES is based on an integrated infrastructure, modular function organization, coherent data handling and fast data access. These key features allow reproducibility of pipeline runs, enable ease of use and provide flexibility to update functions within the pipeline. The cyclic nature of the pipeline offers the possibility to make the science output an efficient measure for design optimization and feasibility testing. We present the architecture, first science, and computational performance results of the simulation pipeline. The framework is general, but for the benchmark tests, we use the Dark Energy Spectrometer (DESpec), one of the early concepts for the upcoming project, the Dark Energy Spectroscopic Instrument (DESI). We discuss how the SPOKES framework enables a rigorous process to optimize and exploit spectroscopic survey experiments in order to derive high-precision cosmological measurements optimally.

The Quijote simulations are a set of 44,100 full N-body simulations spanning more than 7,000 cosmological models in the $\{\Omega_{\rm m}, \Omega_{\rm b}, h, n_s, \sigma_8, M_\nu, w \}$ hyperplane. At a single redshift the simulations contain more than 8.5 trillions of particles over a combined volume of 44,100 $(h^{-1}{\rm Gpc})^3$; each simulation follow the evolution of $256^3$, $512^3$ or $1024^3$ particles in a box of $1~h^{-1}{\rm Gpc}$ length. Billions of dark matter halos and cosmic voids have been identified in the simulations, whose runs required more than 35 million core hours. The Quijote simulations have been designed for two main purposes: 1) to quantify the information content on cosmological observables, and 2) to provide enough data to train machine learning algorithms. In this paper we describe the simulations and show a few of their applications. We also release the Petabyte of data generated, comprising hundreds of thousands of simulation snapshots at multiple redshifts, halo and void catalogs, together with millions of summary statistics such as power spectra, bispectra, correlation functions, marked power spectra, and estimated probability density functions., Comment: 20 pages, 15 figures. Matches published version. Simulations publicly available at https://github.com/franciscovillaescusa/Quijote-simulations

With the high-precision data from current and upcoming experiments, it becomes increasingly important to perform consistency tests of the standard cosmological model. In this work, we focus on consistency measures between different data sets and methods that allow us to assess the goodness of fit of different models. We address both of these questions using the relative entropy or Kullback-Leibler (KL) divergence [Kullback et al., 1951]. First, we revisit the relative entropy as a consistency measure between data sets and further investigate some of its key properties, such as asymmetry and path dependence. We then introduce a novel model rejection framework, which is based on the relative entropy and the posterior predictive distribution. We validate the method on several toy models and apply it to Type Ia supernovae data from the JLA and CMB constraints from Planck 2015, testing the consistency of the data with six different cosmological models., 31 pages, 10 figures, 4 tables, updated following referee's comments, matches version accepted by JCAP

Weak Gravitational Lensing is a powerful probe of the dark sector of the Universe. One of the main challenges for this technique is the treatment of systematics in the measurement of cosmic shear from galaxy shapes. In an earlier work, Refregier & Amara (2014) have proposed the Monte Carlo Control Loops (MCCL) to overcome these effects using a forward modeling approach. We focus here on one of the control loops in this method, the task of which is the calibration of the shear measurement. For this purpose, we first consider the requirements on the shear systematics for a given survey and propagate them to different systematics terms. We use two one-point statistics to calibrate the shear measurement and six further one-point statistics as diagnostics. We also propagate the systematics levels that we estimate from the one-point functions to the two-point functions for the different systematic error sources. This allows us to assess the consistency between the systematics levels measured in different ways. To test the method, we construct synthetic sky surveys with an area of 1,700 deg$^2$. With some simplifying assumptions, we are able to meet the requirements on the shear calibration for this survey configuration. Furthermore, we account for the total residual shear systematics in terms of the contributing sources. We discuss how this MCCL framework can be applied to current and future weak lensing surveys., 23 pages, 4 figures, prepared for submission to JCAP

Upcoming weak lensing surveys will probe large fractions of the sky with unprecedented accuracy. To infer cosmological constraints, a large ensemble of survey simulations are required to accurately model cosmological observables and their covariances. We develop a parallelized multi-lens-plane pipeline called UFalcon, designed to generate full-sky weak lensing maps from lightcones within a minimal runtime. It makes use of L-PICOLA, an approximate numerical code, which provides a fast and accurate alternative to cosmological $N$-Body simulations. The UFalcon maps are constructed by nesting 2 simulations covering a redshift-range from $z=0.1$ to $1.5$ without replicating the simulation volume. We compute the convergence and projected overdensity maps for L-PICOLA in the lightcone or snapshot mode. The generation of such a map, including the L-PICOLA simulation, takes about 3 hours walltime on 220 cores. We use the maps to calculate the spherical harmonic power spectra, which we compare to theoretical predictions and to UFalcon results generated using the full $N$-Body code GADGET-2. We then compute the covariance matrix of the full-sky spherical harmonic power spectra using 150 UFalcon maps based on L-PICOLA in lightcone mode. We consider the PDF, the higher-order moments and the variance of the smoothed field variance to quantify the accuracy of the covariance matrix, which we find to be a few percent for scales $\ell \sim 10^2$ to $10^3$. We test the impact of this level of accuracy on cosmological constraints using an optimistic survey configuration, and find that the final results are robust to this level of uncertainty. The speed and accuracy of our developed pipeline provides a basis to also include further important features such as masking, varying noise and will allow us to compute covariance matrices for models beyond $\Lambda$CDM. [abridged], Comment: 20 pages, 7 figures

Assessing the consistency of parameter constraints derived from different cosmological probes is an important way to test the validity of the underlying cosmological model. In an earlier work [Nicola et al., 2017], we computed constraints on cosmological parameters for $\Lambda$CDM from an integrated analysis of CMB temperature anisotropies and CMB lensing from Planck, galaxy clustering and weak lensing from SDSS, weak lensing from DES SV as well as Type Ia supernovae and Hubble parameter measurements. In this work, we extend this analysis and quantify the concordance between the derived constraints and those derived by the Planck Collaboration as well as WMAP9, SPT and ACT. As a measure for consistency, we use the Surprise statistic [Seehars et al., 2014], which is based on the relative entropy. In the framework of a flat $\Lambda$CDM cosmological model, we find all data sets to be consistent with one another at a level of less than 1$\sigma$. We highlight that the relative entropy is sensitive to inconsistencies in the models that are used in different parts of the analysis. In particular, inconsistent assumptions for the neutrino mass break its invariance on the parameter choice. When consistent model assumptions are used, the data sets considered in this work all agree with each other and $\Lambda$CDM, without evidence for tensions., Comment: 17 pages, 4 figures, 2 tables, updated following referee's comments, now includes discussion of the Riess et al., 2016 Hubble parameter measurement, matches version accepted by JCAP

Recent progress in cosmology has relied on combining different cosmological probes. In earlier work, we implemented an integrated approach to cosmology where the probes are combined into a common framework at the map level. This has the advantage of taking full account of the correlations between the different probes, to provide a stringent test of systematics and of the validity of the cosmological model. We extend this analysis to include not only CMB temperature, galaxy clustering, weak lensing from SDSS but also CMB lensing, weak lensing from the DES SV survey, Type Ia SNe and $H_{0}$ measurements. This yields 12 auto- and cross-power spectra as well as background probes. Furthermore, we extend the treatment of systematic uncertainties. For $��$CDM, we find results that are consistent with our earlier work. Given our enlarged data set and systematics treatment, this confirms the robustness of our analysis and results. Furthermore, we find that our best-fit cosmological model gives a good fit to the data we consider with no signs of tensions within our analysis. We also find our constraints to be consistent with those found by WMAP9, SPT and ACT and the KiDS weak lensing survey. Comparing with the Planck Collaboration results, we see a broad agreement, but there are indications of a tension from the marginalized constraints in most pairs of cosmological parameters. Since our analysis includes CMB temperature Planck data at $10 < \ell < 610$, the tension appears to arise between the Planck high$-\ell$ and the other measurements. Furthermore, we find the constraints on the probe calibration parameters to be in agreement with expectations, showing that the data sets are mutually consistent. In particular, this yields a confirmation of the amplitude calibration of the weak lensing measurements from SDSS, DES SV and Planck CMB lensing from our integrated analysis. [abridged], 30 pages, 19 figures, 3 tables, updated following referee's comments, to appear in PRD

We present measurements of weak gravitational lensing cosmic shear two-point statistics using Dark Energy Survey Science Verification data. We demonstrate that our results are robust to the choice of shear measurement pipeline, either ngmix or im3shape, and robust to the choice of two-point statistic, including both real and Fourier-space statistics. Our results pass a suite of null tests including tests for B-mode contamination and direct tests for any dependence of the two-point functions on a set of 16 observing conditions and galaxy properties, such as seeing, airmass, galaxy color, galaxy magnitude, etc. We furthermore use a large suite of simulations to compute the covariance matrix of the cosmic shear measurements and assign statistical significance to our null tests. We find that our covariance matrix is consistent with the halo model prediction, indicating that it has the appropriate level of halo sample variance. We compare the same jackknife procedure applied to the data and the simulations in order to search for additional sources of noise not captured by the simulations. We find no statistically significant extra sources of noise in the data. The overall detection significance with tomography for our highest source density catalog is 9.7 sigma . Cosmological constraints from the measurements in this work are presented in a companion paper [DES et al., Phys. Rev. D 94, 022001 (2016).].

Recent observational progress has led to the establishment of the standard $\Lambda$CDM model for cosmology. This development is based on different cosmological probes that are usually combined through their likelihoods at the latest stage in the analysis. We implement here an integrated scheme for cosmological probes, which are combined in a common framework starting at the map level. This treatment is necessary as the probes are generally derived from overlapping maps and are thus not independent. It also allows for a thorough test of the cosmological model and of systematics through the consistency of different physical tracers. As a first application, we combine current measurements of the Cosmic Microwave Background (CMB) from the Planck satellite, and galaxy clustering and weak lensing from SDSS. We consider the spherical harmonic power spectra of these probes including all six auto- and cross-correlations along with the associated full Gaussian covariance matrix. This provides an integrated treatment of different analyses usually performed separately including CMB anisotropies, cosmic shear, galaxy clustering, galaxy-galaxy lensing and the Integrated Sachs-Wolfe (ISW) effect with galaxy and shear tracers. We derive constraints on $\Lambda$CDM parameters that are compatible with existing constraints and highlight tensions between data sets, which become apparent in this integrated treatment. We discuss how this approach provides a complete and powerful integrated framework for probe combination and how it can be extended to include other tracers in the context of current and future wide field cosmological surveys., Comment: 29 pages, 19 figures, 3 tables, to appear in PRD, updated following referee's comments including small changes in results

Shear peak statistics has gained a lot of attention recently as a practical alternative to the two point statistics for constraining cosmological parameters. We perform a shear peak statistics analysis of the Dark Energy Survey (DES) Science Verification (SV) data, using weak gravitational lensing measurements from a 139 deg$^2$ field. We measure the abundance of peaks identified in aperture mass maps, as a function of their signal-to-noise ratio, in the signal-to-noise range $04$ would require significant corrections, which is why we do not include them in our analysis. We compare our results to the cosmological constraints from the two point analysis on the SV field and find them to be in good agreement in both the central value and its uncertainty. We discuss prospects for future peak statistics analysis with upcoming DES data., Comment: 21 pages, 14 figures, submitted to MNRAS

Galaxy spectra are essential to probe the spatial distribution of galaxies in our Universe. To better interpret current and future spectroscopic galaxy redshift surveys, it is important to be able to simulate these data sets. We describe Uspec, a forward modeling tool to generate galaxy spectra taking into account some intrinsic galaxy properties as well as instrumental responses of a given telescope. The model for the intrinsic properties of the galaxy population, i.e., the luminosity functions, and size and spectral coefficients distribu- tions, was developed in an earlier work for broad-band imaging surveys [1], and we now aim to test the model further using spectroscopic data. We apply Uspec to the SDSS/CMASS sample of Luminous Red Galaxies (LRGs). We construct selection cuts that match those used to build this LRG sample, which we then apply to data and simulations in the same way. The resulting real and simulated average spectra show a good statistical agreement overall, with residual differences likely coming from a bluer galaxy population of the simulated sam- ple. We also do not explore the impact of non-solar element ratios in our simulations. For a quantitative comparison, we perform Principal Component Analysis (PCA) of the sets of spectra. By comparing the PCs constructed from simulations and data, we find good agree- ment for all components. The distributions of the eigencoefficients also show an appreciable overlap. We are therefore able to properly simulate the LRG sample taking into account the SDSS/BOSS instrumental responses. The differences between the two samples can be ascribed to the intrinsic properties of the simulated galaxy population, which can be reduced by further improvements of our modelling method in the future. We discuss how these results can be useful for the forward modeling of upcoming large spectroscopic surveys., Comment: 32 pages, 14 figures, accepted by JCAP

In light of the growing number of cosmological observations, it is important to develop versatile tools to quantify the constraining power and consistency of cosmological probes. Originally motivated from information theory, we use the relative entropy to compute the information gained by Bayesian updates in units of bits. This measure quantifies both the improvement in precision and the 'surprise', i.e. the tension arising from shifts in central values. Our starting point is a WMAP9 prior which we update with observations of the distance ladder, supernovae (SNe), baryon acoustic oscillations (BAO), and weak lensing as well as the 2015 Planck release. We consider the parameters of the flat $\Lambda$CDM concordance model and some of its extensions which include curvature and Dark Energy equation of state parameter $w$. We find that, relative to WMAP9 and within these model spaces, the probes that have provided the greatest gains are Planck (10 bits), followed by BAO surveys (5.1 bits) and SNe experiments (3.1 bits). The other cosmological probes, including weak lensing (1.7 bits) and {$\rm H_0$} measures (1.7 bits), have contributed information but at a lower level. Furthermore, we do not find any significant surprise when updating the constraints of WMAP9 with any of the other experiments, meaning that they are consistent with WMAP9. However, when we choose Planck15 as the prior, we find that, accounting for the full multi-dimensionality of the parameter space, the weak lensing measurements of CFHTLenS produce a large surprise of 4.4 bits which is statistically significant at the 8 $\sigma$ level. We discuss how the relative entropy provides a versatile and robust framework to compare cosmological probes in the context of current and future surveys., Comment: 26 pages, 5 figures

Determining the redshift distribution $n(z)$ of galaxy samples is essential for several cosmological probes including weak lensing. For imaging surveys, this is usually done using photometric redshifts estimated on an object-by-object basis. We present a new approach for directly measuring the global $n(z)$ of cosmological galaxy samples, including uncertainties, using forward modeling. Our method relies on image simulations produced using UFig (Ultra Fast Image Generator) and on ABC (Approximate Bayesian Computation) within the $MCCL$ (Monte-Carlo Control Loops) framework. The galaxy population is modeled using parametric forms for the luminosity functions, spectral energy distributions, sizes and radial profiles of both blue and red galaxies. We apply exactly the same analysis to the real data and to the simulated images, which also include instrumental and observational effects. By adjusting the parameters of the simulations, we derive a set of acceptable models that are statistically consistent with the data. We then apply the same cuts to the simulations that were used to construct the target galaxy sample in the real data. The redshifts of the galaxies in the resulting simulated samples yield a set of $n(z)$ distributions for the acceptable models. We demonstrate the method by determining $n(z)$ for a cosmic shear like galaxy sample from the 4-band Subaru Suprime-Cam data in the COSMOS field. We also complement this imaging data with a spectroscopic calibration sample from the VVDS survey. We compare our resulting posterior $n(z)$ distributions to the one derived from photometric redshifts estimated using 36 photometric bands in COSMOS and find good agreement. This offers good prospects for applying our approach to current and future large imaging surveys., Comment: 26 pages, 10 figures

Spectroscopic redshift surveys offer great prospects for constraining the dark sector in cosmology. Future surveys will however be both deep and wide and will thus require an analysis in 3-dimensional spherical geometry. We review and compare several methods which have been proposed in the literature for this purpose, focusing in particular on implementations of the spherical harmonic tomography (SHT) power spectrum $C^{i j}_{l}$ and the spherical Fourier Bessel (SFB) power spectrum $C_{l} (k, k')$. Using a Fisher analysis, we compare the forecasted constraints on cosmological parameters using these statistics. These constraints typically rely on approximations such as the Limber approximation and make specific choices in the numerical implementation of each statistic. Using a series of toy models, we explore the applicability of these approximations and study the sensitivity of the SHT and SFB statistics to the details of their implementation. In particular, we show that overlapping redshift bins may improve cosmological constraints using the SHT statistic when the number of bins is small, and that the SFB constraints are quite robust to changes in the assumed distance-redshift relation. We also find that the SHT can be tailored to be more sensitive to modes at redshifts close to the survey boundary, while the SFB appears better suited to capture information beyond the smooth shape of the power spectrum. In this context, we discuss the pros and cons of the different techniques and their impact on the design and analysis of future wide field spectroscopic surveys.

Monthly Notices of the Royal Astronomical Society, 452 (4), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

Optical imaging surveys measure both the galaxy density and the gravitational lensing-induced shear fields across the sky. Recently, the Dark Energy Survey (DES) collaboration used a joint fit to two-point correlations between these observables to place tight constraints on cosmology (DES Collaboration et al. 2017). In this work, we develop the methodology to extend the DES year one joint probes analysis to include cross-correlations of the optical survey observables with gravitational lensing of the cosmic microwave background (CMB) as measured by the South Pole Telescope (SPT) and Planck. Using simulated analyses, we show how the resulting set of five two-point functions increases the robustness of the cosmological constraints to systematic errors in galaxy lensing shear calibration. Additionally, we show that contamination of the SPT+Planck CMB lensing map by the thermal Sunyaev-Zel'dovich effect is a potentially large source of systematic error for two-point function analyses, but show that it can be reduced to acceptable levels in our analysis by masking clusters of galaxies and imposing angular scale cuts on the two-point functions. The methodology developed here will be applied to the analysis of data from the DES, the SPT, and Planck in a companion work., Comment: 21 pages, 11 figures; matches version resubmitted to journal

We measure the cross-correlation between weak lensing of galaxy images and of the cosmic microwave background (CMB). The effects of gravitational lensing on different sources will be correlated if the lensing is caused by the same mass fluctuations. We use galaxy shape measurements from 139 deg$^{2}$ of the Dark Energy Survey (DES) Science Verification data and overlapping CMB lensing from the South Pole Telescope (SPT) and Planck. The DES source galaxies have a median redshift of $z_{\rm med} {\sim} 0.7$, while the CMB lensing kernel is broad and peaks at $z{\sim}2$. The resulting cross-correlation is maximally sensitive to mass fluctuations at $z{\sim}0.44$. Assuming the Planck 2015 best-fit cosmology, the amplitude of the DES$\times$SPT cross-power is found to be $A = 0.88 \pm 0.30$ and that from DES$\times$Planck to be $A = 0.86 \pm 0.39$, where $A=1$ corresponds to the theoretical prediction. These are consistent with the expected signal and correspond to significances of $2.9 \sigma$ and $2.2 \sigma$ respectively. We demonstrate that our results are robust to a number of important systematic effects including the shear measurement method, estimator choice, photometric redshift uncertainty and CMB lensing systematics. Significant intrinsic alignment of galaxy shapes would increase the cross-correlation signal inferred from the data; we calculate a value of $A = 1.08 \pm 0.36$ for DES$\times$SPT when we correct the observations with a simple IA model. With three measurements of this cross-correlation now existing in the literature, there is not yet reliable evidence for any deviation from the expected LCDM level of cross-correlation, given the size of the statistical uncertainties and the significant impact of systematic errors, particularly IAs. We provide forecasts for the expected signal-to-noise of the combination of the five-year DES survey and SPT-3G., Comment: 13 pages, 6 figures

We introduce The Novel Probes Project, an initiative to advance the field of astrophysical tests of the dark sector by creating a forum that connects observers and theorists. This review focuses on tests of gravity and is intended to be of use primarily to observers, but also to theorists with interest in the development of experimental tests. It is twinned with a separate review on tests of dark matter self-interactions (Adhikari et al., in prep.). Our focus is on astrophysical probes of gravity in the weak-field regime, ranging from stars to quasilinear cosmological scales. These are complementary to both strong-field tests and background and linear probes in cosmology. In particular, the nonlinear screening mechanisms that are an integral part of viable modified gravity models lead to characteristic signals specifically on astrophysical scales. The constraining power of these signals is not limited by cosmic variance, but comes with the challenge of building robust theoretical models of the nonlinear dynamics of stars, galaxies, clusters and large scale structure. In this review we lay the groundwork for a thorough exploration of the astrophysical regime with an eye to using the current and next generation of observations for tests of gravity. We begin by setting the scene for how theories beyond General Relativity are expected to behave, focusing primarily on screened fifth forces. We describe the analytic and numerical techniques for exploring the pertinent astrophysical systems, as well as the signatures of modified gravity. With these in hand we present a range of observational tests, and discuss prospects for future measurements and theoretical developments., 108 pages; matches Reviews of Modern Physics accepted version; project website at https://www.novelprobes.org/