Showing 1–50 of 194 results for author: Heavens, A

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

Extreme data compression for Bayesian model comparison

Authors: Alan F. Heavens, Arrykrishna Mootoovaloo, Roberto Trotta, Elena Sellentin

Abstract: We develop extreme data compression for use in Bayesian model comparison via the MOPED algorithm, as well as more general score compression. We find that Bayes factors from data compressed with the MOPED algorithm are identical to those from their uncompressed datasets when the models are linear and the errors Gaussian. In other nonlinear cases, whether nested or not, we find negligible difference… ▽ More We develop extreme data compression for use in Bayesian model comparison via the MOPED algorithm, as well as more general score compression. We find that Bayes factors from data compressed with the MOPED algorithm are identical to those from their uncompressed datasets when the models are linear and the errors Gaussian. In other nonlinear cases, whether nested or not, we find negligible differences in the Bayes factors, and show this explicitly for the Pantheon-SH0ES supernova dataset. We also investigate the sampling properties of the Bayesian Evidence as a frequentist statistic, and find that extreme data compression reduces the sampling variance of the Evidence, but has no impact on the sampling distribution of Bayes factors. Since model comparison can be a very computationally-intensive task, MOPED extreme data compression may present significant advantages in computational time. △ Less

Submitted 13 July, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 15 pages, 5 figures. Invited paper for JCAP 20th anniversary edition

Almanac: MCMC-based signal extraction of power spectra and maps on the sphere

Authors: E. Sellentin, A. Loureiro, L. Whiteway, J. S. Lafaurie, S. T. Balan, M. Olamaie, A. H. Jaffe, A. F. Heavens

Abstract: Inference in cosmology often starts with noisy observations of random fields on the celestial sphere, such as maps of the microwave background radiation, continuous maps of cosmic structure in different wavelengths, or maps of point tracers of the cosmological fields. Almanac uses Hamiltonian Monte Carlo sampling to infer the underlying all-sky noiseless maps of cosmic structures, in multiple reds… ▽ More Inference in cosmology often starts with noisy observations of random fields on the celestial sphere, such as maps of the microwave background radiation, continuous maps of cosmic structure in different wavelengths, or maps of point tracers of the cosmological fields. Almanac uses Hamiltonian Monte Carlo sampling to infer the underlying all-sky noiseless maps of cosmic structures, in multiple redshift bins, together with their auto- and cross-power spectra. It can sample many millions of parameters, handling the highly variable signal-to-noise of typical cosmological signals, and it provides science-ready posterior data products. In the case of spin-weight 2 fields, Almanac infers $E$- and $B$-mode power spectra and parity-violating $EB$ power, and, by sampling the full posteriors rather than point estimates, it avoids the problem of $EB$-leakage. For theories with no $B$-mode signal, inferred non-zero $B$-mode power may be a useful diagnostic of systematic errors or an indication of new physics. Almanac's aim is to characterise the statistical properties of the maps, with outputs that are completely independent of the cosmological model, beyond an assumption of statistical isotropy. Inference of parameters of any particular cosmological model follows in a separate analysis stage. We demonstrate our signal extraction on a CMB-like experiment. △ Less

Submitted 29 August, 2023; v1 submitted 25 May, 2023; originally announced May 2023.

Comments: 27 pages, 18 figures. v2 accepted for publication by The Open Journal of Astrophysics with minor changes. v3 no changes, missing acknowledgement added

Field-level inference of cosmic shear with intrinsic alignments and baryons

Authors: Natalia Porqueres, Alan Heavens, Daniel Mortlock, Guilhem Lavaux, T. Lucas Makinen

Abstract: We construct a field-based Bayesian Hierarchical Model for cosmic shear that includes, for the first time, the important astrophysical systematics of intrinsic alignments and baryon feedback, in addition to a gravity model. We add to the BORG-WL framework the tidal alignment and tidal torquing model (TATT) for intrinsic alignments and compare them with the non-linear alignment (NLA) model. With sy… ▽ More We construct a field-based Bayesian Hierarchical Model for cosmic shear that includes, for the first time, the important astrophysical systematics of intrinsic alignments and baryon feedback, in addition to a gravity model. We add to the BORG-WL framework the tidal alignment and tidal torquing model (TATT) for intrinsic alignments and compare them with the non-linear alignment (NLA) model. With synthetic data, we have shown that adding intrinsic alignments and sampling the TATT parameters does not reduce the constraining power of the method and the field-based approach lifts the weak lensing degeneracy. We add baryon effects at the field level using the enthalpy gradient descent (EGD) model. This model displaces the dark matter particles without knowing whether they belong to a halo and allows for self-calibration of the model parameters, which are inferred from the data. We have also illustrated the effects of model misspecification for the baryons. The resulting model now contains the most important physical effects and is suitable for application to data. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Almanac: Weak Lensing power spectra and map inference on the masked sphere

Authors: A. Loureiro, L. Whiteway, E. Sellentin, J. S. Lafaurie, A. H. Jaffe, A. F. Heavens

Abstract: We present a field-based signal extraction of weak lensing from noisy observations on the curved and masked sky. We test the analysis on a simulated Euclid-like survey, using a Euclid-like mask and noise level. To make optimal use of the information available in such a galaxy survey, we present a Bayesian method for inferring the angular power spectra of the weak lensing fields, together with an i… ▽ More We present a field-based signal extraction of weak lensing from noisy observations on the curved and masked sky. We test the analysis on a simulated Euclid-like survey, using a Euclid-like mask and noise level. To make optimal use of the information available in such a galaxy survey, we present a Bayesian method for inferring the angular power spectra of the weak lensing fields, together with an inference of the noise-cleaned tomographic weak lensing shear and convergence (projected mass) maps. The latter can be used for field-level inference with the aim of extracting cosmological parameter information including non-gaussianity of cosmic fields. We jointly infer all-sky $E$-mode and $B$-mode tomographic auto- and cross-power spectra from the masked sky, and potentially parity-violating $EB$-mode power spectra, up to a maximum multipole of $\ell_{\rm max}=2048$. We use Hamiltonian Monte Carlo sampling, inferring simultaneously the power spectra and denoised maps with a total of $\sim 16.8$ million free parameters. The main output and natural outcome is the set of samples of the posterior, which does not suffer from leakage of power from $E$ to $B$ unless reduced to point estimates. However, such point estimates of the power spectra, the mean and most likely maps, and their variances and covariances, can be computed if desired. △ Less

Submitted 3 February, 2023; v1 submitted 24 October, 2022; originally announced October 2022.

Comments: 18 pages, 16 figures - small changes in response to helpful comments from a referee

Journal ref: The Open Journal of Astrophysics, Vol 6, 2023

Report of the Topical Group on Cosmic Probes of Fundamental Physics for for Snowmass 2021

Authors: Rana X. Adhikari, Luis A. Anchordoqui, Ke Fang, B. S. Sathyaprakash, Kirsten Tollefson, Tiffany R. Lewis, Kristi Engel, Amin Aboubrahim, Ozgur Akarsu, Yashar Akrami, Roberto Aloisio, Rafael Alves Batista, Mario Ballardini, Stefan W. Ballmer, Ellen Bechtol, David Benisty, Emanuele Berti, Simon Birrer, Alexander Bonilla, Richard Brito, Mauricio Bustamante, Robert Caldwell, Vitor Cardoso, Sukanya Chakrabarti, Thomas Y. Chen , et al. (96 additional authors not shown)

Abstract: Cosmic Probes of Fundamental Physics take two primary forms: Very high energy particles (cosmic rays, neutrinos, and gamma rays) and gravitational waves. Already today, these probes give access to fundamental physics not available by any other means, helping elucidate the underlying theory that completes the Standard Model. The last decade has witnessed a revolution of exciting discoveries such as… ▽ More Cosmic Probes of Fundamental Physics take two primary forms: Very high energy particles (cosmic rays, neutrinos, and gamma rays) and gravitational waves. Already today, these probes give access to fundamental physics not available by any other means, helping elucidate the underlying theory that completes the Standard Model. The last decade has witnessed a revolution of exciting discoveries such as the detection of high-energy neutrinos and gravitational waves. The scope for major developments in the next decades is dramatic, as we detail in this report. △ Less

Submitted 23 September, 2022; originally announced September 2022.

Comments: Report of theTopical Group on Cosmic Probes of Fundamental Physics, for the U.S. decadal Particle Physics Planning Exercise (Snowmass 2021)

The Cosmic Graph: Optimal Information Extraction from Large-Scale Structure using Catalogues

Authors: T. Lucas Makinen, Tom Charnock, Pablo Lemos, Natalia Porqueres, Alan Heavens, Benjamin D. Wandelt

Abstract: We present an implicit likelihood approach to quantifying cosmological information over discrete catalogue data, assembled as graphs. To do so, we explore cosmological parameter constraints using mock dark matter halo catalogues. We employ Information Maximising Neural Networks (IMNNs) to quantify Fisher information extraction as a function of graph representation. We a) demonstrate the high sensi… ▽ More We present an implicit likelihood approach to quantifying cosmological information over discrete catalogue data, assembled as graphs. To do so, we explore cosmological parameter constraints using mock dark matter halo catalogues. We employ Information Maximising Neural Networks (IMNNs) to quantify Fisher information extraction as a function of graph representation. We a) demonstrate the high sensitivity of modular graph structure to the underlying cosmology in the noise-free limit, b) show that graph neural network summaries automatically combine mass and clustering information through comparisons to traditional statistics, c) demonstrate that networks can still extract information when catalogues are subject to noisy survey cuts, and d) illustrate how nonlinear IMNN summaries can be used as asymptotically optimal compressed statistics for Bayesian simulation-based inference. We reduce the area of joint $Ω_m, σ_8$ parameter constraints with small ($\sim$100 object) halo catalogues by a factor of 42 over the two-point correlation function, and demonstrate that the networks automatically combine mass and clustering information. This work utilises a new IMNN implementation over graph data in Jax, which can take advantage of either numerical or auto-differentiability. We also show that graph IMNNs successfully compress simulations away from the fiducial model at which the network is fitted, indicating a promising alternative to n-point statistics in catalogue simulation-based analyses. △ Less

Submitted 22 December, 2022; v1 submitted 11 July, 2022; originally announced July 2022.

Comments: 16 pages, 10 figures. Accepted to the Open Journal of Astrophysics. We provide code and a tutorial for the analysis and relevant software at https://github.com/tlmakinen/cosmicGraphs

Field-Based Physical Inference From Peculiar Velocity Tracers

Authors: James Prideaux-Ghee, Florent Leclercq, Guilhem Lavaux, Alan Heavens, Jens Jasche

Abstract: We present a Bayesian hierarchical modelling approach to reconstruct the initial cosmic matter density field constrained by peculiar velocity observations. As our approach features a model for the gravitational evolution of dark matter to connect the initial conditions to late-time observations, it reconstructs the final density and velocity fields as natural byproducts. We implement this field-ba… ▽ More We present a Bayesian hierarchical modelling approach to reconstruct the initial cosmic matter density field constrained by peculiar velocity observations. As our approach features a model for the gravitational evolution of dark matter to connect the initial conditions to late-time observations, it reconstructs the final density and velocity fields as natural byproducts. We implement this field-based physical inference approach by adapting the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm, which explores the high-dimensional posterior through the use of Hamiltonian Monte Carlo sampling. We test the self-consistency of the method using random sets of mock tracers, and assess its accuracy in a more complex scenario where peculiar velocity tracers are non-linearly evolved mock haloes. We find that our framework self-consistently infers the initial conditions, density and velocity fields, and shows some robustness to model mis-specification. As compared to the state-of-the-art approach of constrained Gaussian random fields/Wiener filtering, our method produces more accurate final density and velocity field reconstructions. It also allows us to constrain the initial conditions by peculiar velocity observations, complementing in this aspect previous field-based approaches based on other cosmological observables. △ Less

Submitted 7 December, 2022; v1 submitted 31 March, 2022; originally announced April 2022.

Comments: 23 pages, 15 figures. Accepted for publication in MNRAS

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

Authors: Elcio Abdalla, Guillermo Franco Abellán, Amin Aboubrahim, Adriano Agnello, Ozgur Akarsu, Yashar Akrami, George Alestas, Daniel Aloni, Luca Amendola, Luis A. Anchordoqui, Richard I. Anderson, Nikki Arendse, Marika Asgari, Mario Ballardini, Vernon Barger, Spyros Basilakos, Ronaldo C. Batista, Elia S. Battistelli, Richard Battye, Micol Benetti, David Benisty, Asher Berlin, Paolo de Bernardis, Emanuele Berti, Bohdan Bidenko , et al. (178 additional authors not shown)

Abstract: 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 system… ▽ More 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] △ Less

Submitted 24 April, 2022; v1 submitted 11 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021. 224 pages, 27 figures. Accepted for publication in JHEAp

Journal ref: J. High En. Astrophys. 2204, 002 (2022)

Light bending by the cosmological constant

Authors: Lingyi Hu, Alan Heavens, David Bacon

Abstract: We revisit the question of whether the cosmological constant $Λ$ affects the cosmological gravitational bending of light, by numerical integration of the geodesic equations for a Swiss cheese model consisting of a point mass and a compensated vacuole, in a Friedmann-Robertson-Walker background. We find that there is virtually no dependence of the light bending on the cosmological constant that is… ▽ More We revisit the question of whether the cosmological constant $Λ$ affects the cosmological gravitational bending of light, by numerical integration of the geodesic equations for a Swiss cheese model consisting of a point mass and a compensated vacuole, in a Friedmann-Robertson-Walker background. We find that there is virtually no dependence of the light bending on the cosmological constant that is not already accounted for in the angular diameter distances of the standard lensing equations, plus small modifications that arise because the bending is restricted to a finite region covered by the hole. The residual $Λ$ dependence for a $10^{13}\,M_{\odot}$ lens is at the level of 1 part in $10^7$, and even this might be accounted for by small changes in the hole size evolution as the photon crosses. We therefore conclude that there is no need for modification of the standard cosmological lensing equations in the presence of a cosmological constant. △ Less

Submitted 12 January, 2022; v1 submitted 20 September, 2021; originally announced September 2021.

Comments: 15 pages, 6 figures. Accepted for publication in JCAP

Matching Bayesian and frequentist coverage probabilities when using an approximate data covariance matrix

Authors: Will J. Percival, Oliver Friedrich, Elena Sellentin, Alan Heavens

Abstract: Observational astrophysics consists of making inferences about the Universe by comparing data and models. The credible intervals placed on model parameters are often as important as the maximum a posteriori probability values, as the intervals indicate concordance or discordance between models and with measurements from other data. Intermediate statistics (e.g. the power spectrum) are usually meas… ▽ More Observational astrophysics consists of making inferences about the Universe by comparing data and models. The credible intervals placed on model parameters are often as important as the maximum a posteriori probability values, as the intervals indicate concordance or discordance between models and with measurements from other data. Intermediate statistics (e.g. the power spectrum) are usually measured and inferences made by fitting models to these rather than the raw data, assuming that the likelihood for these statistics has multivariate Gaussian form. The covariance matrix used to calculate the likelihood is often estimated from simulations, such that it is itself a random variable. This is a standard problem in Bayesian statistics, which requires a prior to be placed on the true model parameters and covariance matrix, influencing the joint posterior distribution. As an alternative to the commonly-used Independence-Jeffreys prior, we introduce a prior that leads to a posterior that has approximately frequentist matching coverage. This is achieved by matching the covariance of the posterior to that of the distribution of true values of the parameters around the maximum likelihood values in repeated trials, under certain assumptions. Using this prior, credible intervals derived from a Bayesian analysis can be interpreted approximately as confidence intervals, containing the truth a certain proportion of the time for repeated trials. Linking frequentist and Bayesian approaches that have previously appeared in the astronomical literature, this offers a consistent and conservative approach for credible intervals quoted on model parameters for problems where the covariance matrix is itself an estimate. △ Less

Submitted 1 December, 2021; v1 submitted 23 August, 2021; originally announced August 2021.

Comments: 15 pages, 6 figures, version accepted by MNRAS

Lifting weak lensing degeneracies with a field-based likelihood

Authors: Natalia Porqueres, Alan Heavens, Daniel Mortlock, Guilhem Lavaux

Abstract: We present a field-based approach to the analysis of cosmic shear data to infer jointly cosmological parameters and the dark matter distribution. This forward modelling approach samples the cosmological parameters and the initial matter fluctuations, using a physical gravity model to link the primordial fluctuations to the non-linear matter distribution. Cosmological parameters are sampled and upd… ▽ More We present a field-based approach to the analysis of cosmic shear data to infer jointly cosmological parameters and the dark matter distribution. This forward modelling approach samples the cosmological parameters and the initial matter fluctuations, using a physical gravity model to link the primordial fluctuations to the non-linear matter distribution. Cosmological parameters are sampled and updated consistently through the forward model, varying (1) the initial matter power spectrum, (2) the geometry through the distance-redshift relationship, and (3) the growth of structure and light-cone effects. Our approach extracts more information from the data than methods based on two-point statistics. We find that this field-based approach lifts the strong degeneracy between the cosmological matter density, $Ω_\mathrm{m}$, and the fluctuation amplitude, $σ_8$, providing tight constraints on these parameters from weak lensing data alone. In the simulated four-bin tomographic experiment we consider, the field-based likelihood yields marginal uncertainties on $σ_8$ and $Ω_\mathrm{m}$ that are, respectively, a factor of 3 and 5 smaller than those from a two-point power spectrum analysis applied to the same underlying data. △ Less

Submitted 3 November, 2021; v1 submitted 10 August, 2021; originally announced August 2021.

Comments: Accepted in MNRAS

Quantum coherence of photons at cosmological distances

Authors: Arjun Berera, Suddhasattwa Brahma, Robert Brandenberger, Jaime Calderón-Figueroa, Alan Heavens

Abstract: We identify potential sources of decoherence for $U(1)$ gauge bosons from a cosmological standpoint. Besides interactions with different species in the cosmological medium, we also consider effects due to the expansion of the Universe, which can produce particles (especially scalars) than can potentially interact with the photon in a quantum state. We look in particular at the case of axion-like p… ▽ More We identify potential sources of decoherence for $U(1)$ gauge bosons from a cosmological standpoint. Besides interactions with different species in the cosmological medium, we also consider effects due to the expansion of the Universe, which can produce particles (especially scalars) than can potentially interact with the photon in a quantum state. We look in particular at the case of axion-like particles and their predicted decay channels in our analysis. These interactions are shown to have a negligible effect as far as decoherence goes. Interaction rates with CMB radiation or through Thomson scattering are small, so that the interstellar medium remains the biggest decoherence factor. Thus, quantum teleportation experiments with photon energies in the range $1$-$10$ keV should be feasible at cosmological distances up to the galaxy formation epoch or beyond ($z \sim 100$). △ Less

Submitted 14 July, 2021; originally announced July 2021.

Comments: 11 pages

Journal ref: Phys. Rev. D 104, 063519 (2021)

The Kaiser-Rocket effect: three decades and counting

Authors: Benedict Bahr-Kalus, Daniele Bertacca, Licia Verde, Alan Heavens

Abstract: The peculiar motion of the observer, if not accurately accounted for, is bound to induce a well-defined clustering signal in the distribution of galaxies. This signal is related to the Kaiser rocket effect. Here we examine the amplitude and form of this effect, both analytically and numerically, and discuss possible implications for the analysis and interpretation of forthcoming cosmological surve… ▽ More The peculiar motion of the observer, if not accurately accounted for, is bound to induce a well-defined clustering signal in the distribution of galaxies. This signal is related to the Kaiser rocket effect. Here we examine the amplitude and form of this effect, both analytically and numerically, and discuss possible implications for the analysis and interpretation of forthcoming cosmological surveys. For an idealistic cosmic variance dominated full-sky survey with a Gaussian selection function peaked at $z\sim 1.5$ it is a $> 5σ$ effect and it can in principle bias very significantly the inference of cosmological parameters, especially for primordial non-Gaussianity. For forthcoming surveys, with realistic masks and selection functions, the Kaiser rocket is not a significant concern for cosmological parameter inference. However, it is a systematic effect whose origin, nature and imprint on galaxy maps are well known and thus should be subtracted or mitigated. We present several approaches to do so. △ Less

Submitted 7 November, 2021; v1 submitted 1 July, 2021; originally announced July 2021.

Comments: Version matches the one accepted by JCAP

Kernel-Based Emulator for the 3D Matter Power Spectrum from CLASS

Authors: Arrykrishna Mootoovaloo, Andrew H. Jaffe, Alan F. Heavens, Florent Leclercq

Abstract: The 3D matter power spectrum, $P_δ(k,z)$ is a fundamental quantity in the analysis of cosmological data such as large-scale structure, 21cm observations, and weak lensing. Existing computer models (Boltzmann codes) such as CLASS can provide it at the expense of immoderate computational cost. In this paper, we propose a fast Bayesian method to generate the 3D matter power spectrum, for a given set… ▽ More The 3D matter power spectrum, $P_δ(k,z)$ is a fundamental quantity in the analysis of cosmological data such as large-scale structure, 21cm observations, and weak lensing. Existing computer models (Boltzmann codes) such as CLASS can provide it at the expense of immoderate computational cost. In this paper, we propose a fast Bayesian method to generate the 3D matter power spectrum, for a given set of wavenumbers, $k$ and redshifts, $z$. Our code allows one to calculate the following quantities: the linear matter power spectrum at a given redshift (the default is set to 0); the non-linear 3D matter power spectrum with/without baryon feedback; the weak lensing power spectrum. The gradient of the 3D matter power spectrum with respect to the input cosmological parameters is also returned and this is useful for Hamiltonian Monte Carlo samplers. The derivatives are also useful for Fisher matrix calculations. In our application, the emulator is accurate when evaluated at a set of cosmological parameters, drawn from the prior, with the fractional uncertainty, $ΔP_δ/P_δ$ centred on 0. It is also $\sim 300$ times faster compared to CLASS, hence making the emulator amenable to sampling cosmological and nuisance parameters in a Monte Carlo routine. In addition, once the 3D matter power spectrum is calculated, it can be used with a specific redshift distribution, $n(z)$ to calculate the weak lensing and intrinsic alignment power spectra, which can then be used to derive constraints on cosmological parameters in a weak lensing data analysis problem. The software ($\texttt{emuPK}$) can be trained with any set of points and is distributed on Github, and comes with a pre-trained set of Gaussian Process (GP) models, based on 1000 Latin Hypercube (LH) samples, which follow roughly the current priors for current weak lensing analyses. △ Less

Submitted 8 November, 2021; v1 submitted 5 May, 2021; originally announced May 2021.

Comments: 14 pages, 10 figures, Accepted for publication in Astronomy and Computing Journal

Position-Dependent Correlation Function of Weak Lensing Convergence

Authors: D. Munshi, G. Jung, T. D. Kitching, J. McEwen, M. Liguori, T. Namikawa, A. Heavens

Abstract: We provide a systematic study of the position-dependent correlation function in weak lensing convergence maps and its relation to the squeezed limit of the three-point correlation function (3PCF) using state-of-the-art numerical simulations. We relate the position-dependent correlation function to its harmonic counterpart, i.e., the position-dependent power spectrum or equivalently the integrated… ▽ More We provide a systematic study of the position-dependent correlation function in weak lensing convergence maps and its relation to the squeezed limit of the three-point correlation function (3PCF) using state-of-the-art numerical simulations. We relate the position-dependent correlation function to its harmonic counterpart, i.e., the position-dependent power spectrum or equivalently the integrated bispectrum. We use a recently proposed improved fitting function, BiHalofit, for the bispectrum to compute the theoretical predictions as a function of source redshifts. In addition to low redshift results ($z_s=1.0-2.0$), we also provide results for maps inferred from lensing of the cosmic microwave background, i.e., $z_s=1100$. We include a {\em Euclid}-type realistic survey mask and noise. In agreement with the recent studies on the position-dependent power spectrum, we find that the results from simulations are consistent with the theoretical expectations when appropriate corrections are included. Performing a rough estimate, we find that the (S/N) for the detection of the position-dependent correlation function from {\em Euclid}-type mask with $f_{sky}=0.35$, can range between $6-12$ depending on the value of the intrinsic ellipticity distribution parameter $σ_ε = 0.3-1.0$. For reconstructed $κ$ maps using an ideal CMB survey the (S/N) $\approx 1.8$. We also found that a $10\%$ deviation in $σ_8$ can be detected using IB for the optimistic case of $σ_ε=0.3$ with a (S/N) $\approx 5$. The (S/N) for such detection in case of $Ω_M$ is lower. △ Less

Submitted 17 January, 2023; v1 submitted 2 April, 2021; originally announced April 2021.

Comments: 7 pages, 7 figures (PRD in press)

On the accuracy and precision of correlation functions and field-level inference in cosmology

Authors: Florent Leclercq, Alan Heavens

Abstract: We present a comparative study of the accuracy and precision of correlation function methods and full-field inference in cosmological data analysis. To do so, we examine a Bayesian hierarchical model that predicts log-normal fields and their two-point correlation function. Although a simplified analytic model, the log-normal model produces fields that share many of the essential characteristics of… ▽ More We present a comparative study of the accuracy and precision of correlation function methods and full-field inference in cosmological data analysis. To do so, we examine a Bayesian hierarchical model that predicts log-normal fields and their two-point correlation function. Although a simplified analytic model, the log-normal model produces fields that share many of the essential characteristics of the present-day non-Gaussian cosmological density fields. We use three different statistical techniques: (i) a standard likelihood-based analysis of the two-point correlation function; (ii) a likelihood-free (simulation-based) analysis of the two-point correlation function; (iii) a field-level analysis, made possible by the more sophisticated data assimilation technique. We find that (a) standard assumptions made to write down a likelihood for correlation functions can cause significant biases, a problem that is alleviated with simulation-based inference; and (b) analysing the entire field offers considerable advantages over correlation functions, through higher accuracy, higher precision, or both. The gains depend on the degree of non-Gaussianity, but in all cases, including for weak non-Gaussianity, the advantage of analysing the full field is substantial. △ Less

Submitted 30 July, 2021; v1 submitted 6 March, 2021; originally announced March 2021.

Comments: 6+8 pages, 4+5 figures. Matches MNRAS Letters published version. Appendices provide supplementary information, including calculations of Fisher matrices. Our code and data are publicly available at https://github.com/florent-leclercq/correlations_vs_field

Journal ref: MNRAS Letters 506, L85 (2021)

The integrated angular bispectrum of weak lensing

Authors: Gabriel Jung, Toshiya Namikawa, Michele Liguori, Dipak Munshi, Alan Heavens

Abstract: We investigate three-point statistics in weak lensing convergence, through the integrated bispectrum. This statistic involves measuring power spectra in patches, and is thus easy to measure, and avoids the complexity of estimating the very large number of possible bispectrum configurations. The integrated bispectrum principally probes the squeezed limit of the bispectrum. To be useful as a set of… ▽ More We investigate three-point statistics in weak lensing convergence, through the integrated bispectrum. This statistic involves measuring power spectra in patches, and is thus easy to measure, and avoids the complexity of estimating the very large number of possible bispectrum configurations. The integrated bispectrum principally probes the squeezed limit of the bispectrum. To be useful as a set of summary statistics, accurate theoretical predictions of the signal are required, and, assuming Gaussian sampling distributions, the covariance matrix. In this paper, we investigate through simulations how accurate are theoretical formulae for both the integrated bispectrum and its covariance, finding that there a small inaccuracies in the theoretical signal, and more serious deviations in the covariance matrix, which may need to be estimated using simulations. △ Less

Submitted 12 July, 2021; v1 submitted 10 February, 2021; originally announced February 2021.

Comments: 23 pages, 12 figures. v2: minor textual corrections and additions, matches published version

Journal ref: JCAP06(2021)055

Bayesian forward modelling of cosmic shear data

Authors: Natalia Porqueres, Alan Heavens, Daniel Mortlock, Guilhem Lavaux

Abstract: We present a Bayesian hierarchical modelling approach to infer the cosmic matter density field, and the lensing and the matter power spectra, from cosmic shear data. This method uses a physical model of cosmic structure formation to infer physically plausible cosmic structures, which accounts for the non-Gaussian features of the gravitationally evolved matter distribution and light-cone effects. W… ▽ More We present a Bayesian hierarchical modelling approach to infer the cosmic matter density field, and the lensing and the matter power spectra, from cosmic shear data. This method uses a physical model of cosmic structure formation to infer physically plausible cosmic structures, which accounts for the non-Gaussian features of the gravitationally evolved matter distribution and light-cone effects. We test and validate our framework with realistic simulated shear data, demonstrating that the method recovers the unbiased matter distribution and the correct lensing and matter power spectrum. While the cosmology is fixed in this test, and the method employs a prior power spectrum, we demonstrate that the lensing results are sensitive to the true power spectrum when this differs from the prior. In this case, the density field samples are generated with a power spectrum that deviates from the prior, and the method recovers the true lensing power spectrum. The method also recovers the matter power spectrum across the sky, but as currently implemented, it cannot determine the radial power since isotropy is not imposed. In summary, our method provides physically plausible inference of the dark matter distribution from cosmic shear data, allowing us to extract information beyond the two-point statistics and exploiting the full information content of the cosmological fields. △ Less

Submitted 21 January, 2021; v1 submitted 13 November, 2020; originally announced November 2020.

Comments: 10 pages, 13 figures. Accepted for publication in MNRAS

Cosmology Intertwined IV: The Age of the Universe and its Curvature

Authors: Eleonora Di Valentino, Luis A. Anchordoqui, Ozgur Akarsu, Yacine Ali-Haimoud, Luca Amendola, Nikki Arendse, Marika Asgari, Mario Ballardini, Spyros Basilakos, Elia Battistelli, Micol Benetti, Simon Birrer, François R. Bouchet, Marco Bruni, Erminia Calabrese, David Camarena, Salvatore Capozziello, Angela Chen, Jens Chluba, Anton Chudaykin, Eoin Ó Colgáin, Francis-Yan Cyr-Racine, Paolo de Bernardis, Javier de Cruz Pérez, Jacques Delabrouille , et al. (66 additional authors not shown)

Abstract: A precise measurement of the curvature of the Universe is of primeval importance for cosmology since it could not only confirm the paradigm of primordial inflation but also help in discriminating between different early Universe scenarios. The recent observations, while broadly consistent with a spatially flat standard $Λ$ Cold Dark Matter ($Λ$CDM) model, are showing tensions that still allow (and… ▽ More A precise measurement of the curvature of the Universe is of primeval importance for cosmology since it could not only confirm the paradigm of primordial inflation but also help in discriminating between different early Universe scenarios. The recent observations, while broadly consistent with a spatially flat standard $Λ$ Cold Dark Matter ($Λ$CDM) model, are showing tensions that still allow (and, in some cases, even suggest) a few percent deviations from a flat universe. In particular, the Planck Cosmic Microwave Background power spectra, assuming the nominal likelihood, prefer a closed universe at more than 99\% confidence level. While new physics could be in action, this anomaly may be the result of an unresolved systematic error or just a statistical fluctuation. However, since a positive curvature allows a larger age of the Universe, an accurate determination of the age of the oldest objects provides a smoking gun in confirming or falsifying the current flat $Λ$CDM model. △ Less

Submitted 13 October, 2020; v1 submitted 25 August, 2020; originally announced August 2020.

Comments: Snowmass2021 - Letter of Interest

Journal ref: Astropart. Phys. 131, 102607 (2021)

Cosmology Intertwined III: $f σ_8$ and $S_8$

Authors: Eleonora Di Valentino, Luis A. Anchordoqui, Ozgur Akarsu, Yacine Ali-Haimoud, Luca Amendola, Nikki Arendse, Marika Asgari, Mario Ballardini, Spyros Basilakos, Elia Battistelli, Micol Benetti, Simon Birrer, François R. Bouchet, Marco Bruni, Erminia Calabrese, David Camarena, Salvatore Capozziello, Angela Chen, Jens Chluba, Anton Chudaykin, Eoin Ó Colgáin, Francis-Yan Cyr-Racine, Paolo de Bernardis, Javier de Cruz Pérez, Jacques Delabrouille , et al. (67 additional authors not shown)

Abstract: The standard $Λ$ Cold Dark Matter cosmological model provides a wonderful fit to current cosmological data, but a few tensions and anomalies became statistically significant with the latest data analyses. While these anomalies could be due to the presence of systematic errors in the experiments, they could also indicate the need for new physics beyond the standard model. In this Letter of Interest… ▽ More The standard $Λ$ Cold Dark Matter cosmological model provides a wonderful fit to current cosmological data, but a few tensions and anomalies became statistically significant with the latest data analyses. While these anomalies could be due to the presence of systematic errors in the experiments, they could also indicate the need for new physics beyond the standard model. In this Letter of Interest we focus on the tension of the Planck 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 for solving this tension, and we discuss the importance of trying to fit with a single model a full array of data and not just one parameter at a time. △ Less

Submitted 13 October, 2020; v1 submitted 25 August, 2020; originally announced August 2020.

Comments: Snowmass2021 - Letter of Interest

Journal ref: Astropart. Phys. 131, 102604 (2021)

Cosmology Intertwined II: The Hubble Constant Tension

Authors: Eleonora Di Valentino, Luis A. Anchordoqui, Ozgur Akarsu, Yacine Ali-Haimoud, Luca Amendola, Nikki Arendse, Marika Asgari, Mario Ballardini, Spyros Basilakos, Elia Battistelli, Micol Benetti, Simon Birrer, François R. Bouchet, Marco Bruni, Erminia Calabrese, David Camarena, Salvatore Capozziello, Angela Chen, Jens Chluba, Anton Chudaykin, Eoin Ó Colgáin, Francis-Yan Cyr-Racine, Paolo de Bernardis, Javier de Cruz Pérez, Jacques Delabrouille , et al. (68 additional authors not shown)

Abstract: The current cosmological probes have provided a fantastic confirmation of the standard $Λ$ Cold Dark Matter cosmological model, that has been constrained with unprecedented accuracy. However, with the increase of the experimental sensitivity a few statistically significant tensions between different independent cosmological datasets emerged. While these tensions can be in portion the result of sys… ▽ More The current cosmological probes have provided a fantastic confirmation of the standard $Λ$ Cold Dark Matter cosmological model, that has been constrained with unprecedented accuracy. However, with the increase of the experimental sensitivity a few statistically significant tensions between different independent cosmological datasets emerged. While these tensions can be in portion the result of systematic errors, the persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the need for new physics. In this Letter of Interest we will focus on the $4.4σ$ tension between the Planck 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 will list a few interesting new physics models that could solve this tension and discuss how the next decade experiments will be crucial. △ Less

Submitted 13 October, 2020; v1 submitted 25 August, 2020; originally announced August 2020.

Comments: Snowmass2021 - Letter of Interest

Journal ref: Astropart. Phys. 131, 102605 (2021)

Cosmology Intertwined I: Perspectives for the Next Decade

Authors: Eleonora Di Valentino, Luis A. Anchordoqui, Ozgur Akarsu, Yacine Ali-Haimoud, Luca Amendola, Nikki Arendse, Marika Asgari, Mario Ballardini, Spyros Basilakos, Elia Battistelli, Micol Benetti, Simon Birrer, François R. Bouchet, Marco Bruni, Erminia Calabrese, David Camarena, Salvatore Capozziello, Angela Chen, Jens Chluba, Anton Chudaykin, Eoin Ó Colgáin, Francis-Yan Cyr-Racine, Paolo de Bernardis, Javier de Cruz Pérez, Jacques Delabrouille , et al. (67 additional authors not shown)

Abstract: The standard $Λ$ Cold Dark Matter cosmological model provides an amazing description of a wide range of astrophysical and astronomical data. However, there are a few big open questions, that make the standard model look like a first-order approximation to a more realistic scenario that still needs to be fully understood. In this Letter of Interest we will list a few important goals that need to be… ▽ More The standard $Λ$ Cold Dark Matter cosmological model provides an amazing description of a wide range of astrophysical and astronomical data. However, there are a few big open questions, that make the standard model look like a first-order approximation to a more realistic scenario that still needs to be fully understood. In this Letter of Interest we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances present between the different cosmological probes, as the Hubble constant $H_0$ value, the $σ_8 - S_8$ tension, and the anomalies present in the Planck results. Finally, we will give an overview of upgraded experiments and next-generation space-missions and facilities on Earth, that will be of crucial importance to address all these questions. △ Less

Submitted 13 October, 2020; v1 submitted 25 August, 2020; originally announced August 2020.

Comments: Snowmass2021 - Letter of Interest

Journal ref: Astropart. Phys. 131, 102606 (2021)

Extreme data compression while searching for new physics

Authors: Alan Heavens, Elena Sellentin, Andrew Jaffe

Abstract: Bringing a high-dimensional dataset into science-ready shape is a formidable challenge that often necessitates data compression. Compression has accordingly become a key consideration for contemporary cosmology, affecting public data releases, and reanalyses searching for new physics. However, data compression optimized for a particular model can suppress signs of new physics, or even remove them… ▽ More Bringing a high-dimensional dataset into science-ready shape is a formidable challenge that often necessitates data compression. Compression has accordingly become a key consideration for contemporary cosmology, affecting public data releases, and reanalyses searching for new physics. However, data compression optimized for a particular model can suppress signs of new physics, or even remove them altogether. We therefore provide a solution for exploring new physics \emph{during} data compression. In particular, we store additional agnostic compressed data points, selected to enable precise constraints of non-standard physics at a later date. Our procedure is based on the maximal compression of the MOPED algorithm, which optimally filters the data with respect to a baseline model. We select additional filters, based on a generalised principal component analysis, which are carefully constructed to scout for new physics at high precision and speed. We refer to the augmented set of filters as MOPED-PC. They enable an analytic computation of Bayesian evidences that may indicate the presence of new physics, and fast analytic estimates of best-fitting parameters when adopting a specific non-standard theory, without further expensive MCMC analysis. As there may be large numbers of non-standard theories, the speed of the method becomes essential. Should no new physics be found, then our approach preserves the precision of the standard parameters. As a result, we achieve very rapid and maximally precise constraints of standard and non-standard physics, with a technique that scales well to large dimensional datasets. △ Less

Submitted 18 August, 2020; v1 submitted 11 June, 2020; originally announced June 2020.

Comments: 12 pages, 12 figures. Accepted for publication in MNRAS

The distribution of dark galaxies and spin bias

Authors: Raul Jimenez, Alan F. Heavens

Abstract: In the light of the discovery of numerous (almost) dark galaxies from the ALFALAFA and LITTLE THINGS surveys, we revisit the predictions of Jimenez et al. 1997, based on the Toomre stability of rapidly-spinning gas disks. We have updated the predictions for $Λ$CDM with parameters given by Planck18, computing the expected number densities of dark objects, and their spin parameter and mass distribut… ▽ More In the light of the discovery of numerous (almost) dark galaxies from the ALFALAFA and LITTLE THINGS surveys, we revisit the predictions of Jimenez et al. 1997, based on the Toomre stability of rapidly-spinning gas disks. We have updated the predictions for $Λ$CDM with parameters given by Planck18, computing the expected number densities of dark objects, and their spin parameter and mass distributions. Comparing with the data is more challenging, but where the spins are more reliably determined, the spins are close to the threshold for disks to be stable according to the Toomre criterion, where the expected number density is highest, and reinforces the concept that there is a bias in the formation of luminous galaxies based on the spin of their parent halo. △ Less

Submitted 1 August, 2020; v1 submitted 24 May, 2020; originally announced May 2020.

Comments: matches published version in MNRAS Letters; theoretical model improved, results unchanged

Journal ref: Monthly Notices of the Royal Astronomical Society Letters. 498, L93, 2020

Parameter Inference for Weak Lensing using Gaussian Processes and MOPED

Authors: Arrykrishna Mootoovaloo, Alan F. Heavens, Andrew H. Jaffe, Florent Leclercq

Abstract: In this paper, we propose a Gaussian Process (GP) emulator for the calculation of a) tomographic weak lensing band-power spectra, and b) coefficients of summary data massively compressed with the MOPED algorithm. In the former case cosmological parameter inference is accelerated by a factor of $\sim 10$-$30$ compared to explicit calls to the Boltzmann solver CLASS when applied to KiDS-450 weak len… ▽ More In this paper, we propose a Gaussian Process (GP) emulator for the calculation of a) tomographic weak lensing band-power spectra, and b) coefficients of summary data massively compressed with the MOPED algorithm. In the former case cosmological parameter inference is accelerated by a factor of $\sim 10$-$30$ compared to explicit calls to the Boltzmann solver CLASS when applied to KiDS-450 weak lensing data. Much larger gains will come with future data, where with MOPED compression, the speed up can be up to a factor of $\sim 10^3$ when the common Limber approximation is used. Furthermore, the GP opens up the possibility of dropping the Limber approximation, without which the theoretical calculations may be unfeasibly slow. A potential advantage of GPs is that an error on the emulated function can be computed and this uncertainty incorporated into the likelihood. If speed is of the essence, then the mean of the Gaussian Process can be used and the uncertainty ignored. We compute the Kullback-Leibler divergence between the emulator likelihood and the CLASS likelihood, and on the basis of this and from analysing the uncertainties on the parameters, we find that in this case, the inclusion of the GP uncertainty does not justify the extra computational expense in the test application. For future weak lensing surveys such as Euclid and Legacy Survey of Space and Telescope (LSST), the number of summary statistics will be large, up to $\sim 10^{4}$. The speed of MOPED is determined by the number of parameters, not the number of summary data, so the gains are very large. In the non-Limber case, the speed-up can be a factor of $\sim 10^5$, provided that a fast way to compute the theoretical MOPED coefficients is available. The GP presented here provides such a fast mechanism and enables MOPED to be employed. △ Less

Submitted 15 July, 2020; v1 submitted 13 May, 2020; originally announced May 2020.

Comments: 16 pages, 11 figures (Accepted for publication in MNRAS)

Perfectly parallel cosmological simulations using spatial comoving Lagrangian acceleration

Authors: Florent Leclercq, Baptiste Faure, Guilhem Lavaux, Benjamin D. Wandelt, Andrew H. Jaffe, Alan F. Heavens, Will J. Percival, Camille Noûs

Abstract: Existing cosmological simulation methods lack a high degree of parallelism due to the long-range nature of the gravitational force, which limits the size of simulations that can be run at high resolution. To solve this problem, we propose a new, perfectly parallel approach to simulate cosmic structure formation, which is based on the spatial COmoving Lagrangian Acceleration (sCOLA) framework. Buil… ▽ More Existing cosmological simulation methods lack a high degree of parallelism due to the long-range nature of the gravitational force, which limits the size of simulations that can be run at high resolution. To solve this problem, we propose a new, perfectly parallel approach to simulate cosmic structure formation, which is based on the spatial COmoving Lagrangian Acceleration (sCOLA) framework. Building upon a hybrid analytical and numerical description of particles' trajectories, our algorithm allows for an efficient tiling of a cosmological volume, where the dynamics within each tile is computed independently. As a consequence, the degree of parallelism is equal to the number of tiles. We optimised the accuracy of sCOLA through the use of a buffer region around tiles and of appropriate Dirichlet boundary conditions around sCOLA boxes. As a result, we show that cosmological simulations at the degree of accuracy required for the analysis of the next generation of surveys can be run in drastically reduced wall-clock times and with very low memory requirements. The perfect scalability of our algorithm unlocks profoundly new possibilities for computing larger cosmological simulations at high resolution, taking advantage of a variety of hardware architectures. △ Less

Submitted 16 September, 2022; v1 submitted 10 March, 2020; originally announced March 2020.

Comments: 24 pages, 8 figures, 2 tables. Typos corrected with respect to A&A published version. The code is publicly available at https://bitbucket.org/florent-leclercq/simbelmyne/

Journal ref: A&A 639, A91 (2020)

Gaussian Mixture Models for Blended Photometric Redshifts

Authors: Daniel M. Jones, Alan F. Heavens

Abstract: Future cosmological galaxy surveys such as the Large Synoptic Survey Telescope (LSST) will photometrically observe very large numbers of galaxies. Without spectroscopy, the redshifts required for the analysis of these data will need to be inferred using photometric redshift techniques that are scalable to large sample sizes. The high number density of sources will also mean that around half are bl… ▽ More Future cosmological galaxy surveys such as the Large Synoptic Survey Telescope (LSST) will photometrically observe very large numbers of galaxies. Without spectroscopy, the redshifts required for the analysis of these data will need to be inferred using photometric redshift techniques that are scalable to large sample sizes. The high number density of sources will also mean that around half are blended. We present a Bayesian photometric redshift method for blended sources that uses Gaussian mixture models to learn the joint flux-redshift distribution from a set of unblended training galaxies, and Bayesian model comparison to infer the number of galaxies comprising a blended source. The use of Gaussian mixture models renders both of these applications computationally efficient and therefore suitable for upcoming galaxy surveys. △ Less

Submitted 4 October, 2019; v1 submitted 24 July, 2019; originally announced July 2019.

Comments: 23 pages, updated to version accepted in MNRAS

Measuring the Homogeneity of the Universe Using Polarization Drift

Authors: Raul Jimenez, Roy Maartens, Ali Rida Khalifeh, Robert R. Caldwell, Alan F. Heavens, Licia Verde

Abstract: We propose a method to probe the homogeneity of a general universe, without assuming symmetry. We show that isotropy can be tested at remote locations on the past lightcone by comparing the line-of-sight and transverse expansion rates, using the time dependence of the polarization of Cosmic Microwave Background photons that have been inverse-Compton scattered by the hot gas in massive clusters of… ▽ More We propose a method to probe the homogeneity of a general universe, without assuming symmetry. We show that isotropy can be tested at remote locations on the past lightcone by comparing the line-of-sight and transverse expansion rates, using the time dependence of the polarization of Cosmic Microwave Background photons that have been inverse-Compton scattered by the hot gas in massive clusters of galaxies. This probes a combination of remote transverse and parallel components of the expansion rate of the metric, and we may use radial baryon acoustic oscillations or cosmic clocks to measure the parallel expansion rate. Thus we can test remote isotropy, which is a key requirement of a homogeneous universe. We provide explicit formulas that connect observables and properties of the metric. △ Less

Submitted 20 May, 2019; v1 submitted 28 February, 2019; originally announced February 2019.

Comments: Matches accepted version to JCAP; conclusions unchanged; now includes detailed section of the expected S/N of the signal

Journal ref: JCAP05(2019)048

Primordial power spectrum and cosmology from black-box galaxy surveys

Authors: Florent Leclercq, Wolfgang Enzi, Jens Jasche, Alan Heavens

Abstract: We propose a new, likelihood-free approach to inferring the primordial matter power spectrum and cosmological parameters from arbitrarily complex forward models of galaxy surveys where all relevant statistics can be determined from numerical simulations, i.e. black-boxes. Our approach, which we call simulator expansion for likelihood-free inference (SELFI), builds upon approximate Bayesian computa… ▽ More We propose a new, likelihood-free approach to inferring the primordial matter power spectrum and cosmological parameters from arbitrarily complex forward models of galaxy surveys where all relevant statistics can be determined from numerical simulations, i.e. black-boxes. Our approach, which we call simulator expansion for likelihood-free inference (SELFI), builds upon approximate Bayesian computation using a novel effective likelihood, and upon the linearisation of black-box models around an expansion point. Consequently, we obtain simple "filter equations" for an effective posterior of the primordial power spectrum, and a straightforward scheme for cosmological parameter inference. We demonstrate that the workload is computationally tractable, fixed a priori, and perfectly parallel. As a proof of concept, we apply our framework to a realistic synthetic galaxy survey, with a data model accounting for physical structure formation and incomplete and noisy galaxy observations. In doing so, we show that the use of non-linear numerical models allows the galaxy power spectrum to be safely fitted up to at least $k_\mathrm{max} = 0.5$ $h$/Mpc, outperforming state-of-the-art backward-modelling techniques by a factor of $\sim 5$ in the number of modes used. The result is an unbiased inference of the primordial matter power spectrum across the entire range of scales considered, including a high-fidelity reconstruction of baryon acoustic oscillations. It translates into an unbiased and robust inference of cosmological parameters. Our results pave the path towards easy applications of likelihood-free simulation-based inference in cosmology. We have made our code pySELFI and our data products publicly available at http://pyselfi.florent-leclercq.eu. △ Less

Submitted 8 October, 2019; v1 submitted 26 February, 2019; originally announced February 2019.

Comments: 22 pages, 12 figures. Matches MNRAS published version. The main results are equations (25), (26), and figure 10. The pySELFI code is publicly available on the first author's website (currently http://pyselfi.florent-leclercq.eu) and on GitHub (https://github.com/florent-leclercq/pyselfi/)

Fast Sampling from Wiener Posteriors for Image Data with Dataflow Engines

Authors: Niall Jeffrey, Alan F. Heavens, Philip D. Fortio

Abstract: We use Dataflow Engines (DFE) to construct an efficient Wiener filter of noisy and incomplete image data, and to quickly draw probabilistic samples of the compatible true underlying images from the Wiener posterior. Dataflow computing is a powerful approach using reconfigurable hardware, which can be deeply pipelined and is intrinsically parallel. The unique Wiener-filtered image is the minimum-va… ▽ More We use Dataflow Engines (DFE) to construct an efficient Wiener filter of noisy and incomplete image data, and to quickly draw probabilistic samples of the compatible true underlying images from the Wiener posterior. Dataflow computing is a powerful approach using reconfigurable hardware, which can be deeply pipelined and is intrinsically parallel. The unique Wiener-filtered image is the minimum-variance linear estimate of the true image (if the signal and noise covariances are known) and the most probable true image (if the signal and noise are Gaussian distributed). However, many images are compatible with the data with different probabilities, given by the analytic posterior probability distribution referred to as the Wiener posterior. The DFE code also draws large numbers of samples of true images from this posterior, which allows for further statistical analysis. Naive computation of the Wiener-filtered image is impractical for large datasets, as it scales as $n^3$, where $n$ is the number of pixels. We use a messenger field algorithm, which is well suited to a DFE implementation, to draw samples from the Wiener posterior, that is, with the correct probability we draw samples of noiseless images that are compatible with the observed noisy image. The Wiener-filtered image can be obtained by a trivial modification of the algorithm. We demonstrate a lower bound on the speed-up, from drawing 10$^5$ samples of a 128$^2$ image, of 11.3 ${\pm}$ 0.8 with 8 DFEs in a 1U MPC-X box when compared with a 1U server presenting 32 CPU threads. We also discuss a potential application in astronomy, to provide better dark matter maps and improved determination of the parameters of the Universe. △ Less

Submitted 5 October, 2018; originally announced October 2018.

Comments: 20 pages, 5 figures, accepted by Astronomy and Computing

The gravitational and lensing-ISW bispectrum of 21cm radiation

Authors: Claude J. Schmit, Alan F. Heavens, Jonathan R. Pritchard

Abstract: Cosmic Microwave Background experiments from COBE to Planck, have launched cosmology into an era of precision science, where many cosmological parameters are now determined to the percent level. Next generation telescopes, focussing on the cosmological 21cm signal from neutral hydrogen, will probe enormous volumes in the low-redshift Universe, and have the potential to determine dark energy proper… ▽ More Cosmic Microwave Background experiments from COBE to Planck, have launched cosmology into an era of precision science, where many cosmological parameters are now determined to the percent level. Next generation telescopes, focussing on the cosmological 21cm signal from neutral hydrogen, will probe enormous volumes in the low-redshift Universe, and have the potential to determine dark energy properties and test modifications of Einstein's gravity. We study the 21cm bispectrum due to gravitational collapse as well as the contribution by line of sight perturbations in the form of the lensing-ISW bispectrum at low-redshifts ($z \sim 0.35-3$), targeted by upcoming neutral hydrogen intensity mapping experiments. We compute the expected bispectrum amplitudes and use a Fisher forecast model to compare power spectrum and bispectrum observations of intensity mapping surveys by CHIME, MeerKAT and SKA-mid. We find that combined power spectrum and bispectrum observations have the potential to decrease errors on the cosmological parameters by an order of magnitude compared to Planck. Finally, we compute the contribution of the lensing-ISW bispectrum, and find that, unlike for the cosmic microwave background analyses, it can safely be ignored for 21cm bispectrum observations. △ Less

Submitted 10 December, 2018; v1 submitted 1 October, 2018; originally announced October 2018.

Comments: 19 pages, 10 figures

Bayesian Evidence against Harrison-Zel'dovich spectrum in tension cosmology

Authors: Eleonora Di Valentino, Alessandro Melchiorri, Yabebal Fantaye, Alan Heavens

Abstract: Current cosmological constraints on the scalar spectral index of primordial fluctuations $n_{\rm s}$ in the $Λ$CDM model have excluded the minimal scale-invariant Harrison-Zel'dovich model ($n_{\rm s}=1$; hereafter HZ) at high significance, providing support for inflation. In recent years, however, some tensions have emerged between different cosmological datasets that, if not due to systematics,… ▽ More Current cosmological constraints on the scalar spectral index of primordial fluctuations $n_{\rm s}$ in the $Λ$CDM model have excluded the minimal scale-invariant Harrison-Zel'dovich model ($n_{\rm s}=1$; hereafter HZ) at high significance, providing support for inflation. In recent years, however, some tensions have emerged between different cosmological datasets that, if not due to systematics, could indicate the presence of new physics beyond the $Λ$CDM model. In the light of these developments, we evaluate the Bayesian evidence against HZ in different data combinations and model extensions. Considering only the Planck temperature data, we find inconclusive evidence against HZ when including variations in the neutrino number $N_{\rm eff}$ and/or the Helium abundance $Y_{\rm He}$. Adding the Planck polarization data, on the other hand, yields strong evidence against HZ in the extensions we considered. Perhaps most interestingly, Planck temperature data combined with local measurements of the Hubble constant give as the most probable model an HZ spectrum, with additional neutrinos. However, with the inclusion of polarisation, standard $Λ$CDM is once again preferred, but the HZ model with extra neutrinos is not strongly disfavored. The possibility of fully ruling out the HZ spectrum is therefore ultimately connected with the solution to current tensions between cosmological datasets. If these tensions are confirmed by future data, then new physical mechanisms could be at work and an HZ spectrum could still offer a valid alternative. △ Less

Submitted 28 August, 2018; originally announced August 2018.

Comments: 14 pages, 1 figure

Journal ref: Phys. Rev. D 98, 063508, 2018

Bayesian photometric redshifts of blended sources

Authors: Daniel M. Jones, Alan F. Heavens

Abstract: Photometric redshifts are necessary for enabling large-scale multicolour galaxy surveys to interpret their data and constrain cosmological parameters. While the increased depth of future surveys such as the Large Synoptic Survey Telescope (LSST) will produce higher precision constraints, it will also increase the fraction of sources that are blended. In this paper, we present a Bayesian photometri… ▽ More Photometric redshifts are necessary for enabling large-scale multicolour galaxy surveys to interpret their data and constrain cosmological parameters. While the increased depth of future surveys such as the Large Synoptic Survey Telescope (LSST) will produce higher precision constraints, it will also increase the fraction of sources that are blended. In this paper, we present a Bayesian photometric redshift method for blended sources with an arbitrary number of intrinsic components. This method generalises existing template-based Bayesian photometric redshift (BPZ) methods, and produces joint posterior distributions for the component redshifts that allow uncertainties to be propagated in a principled way. Using Bayesian model comparison, we infer the probability that a source is blended and the number of components that it contains. We extend our formalism to the case where sources are blended in some bands and resolved in others. Applying this to the combination of LSST- and Euclid-like surveys, we find that the addition of resolved photometry results in a significant improvement in the reduction of outliers over the fully-blended case. We make available blendz, a Python implementation of our method. △ Less

Submitted 30 November, 2018; v1 submitted 8 August, 2018; originally announced August 2018.

Comments: Added section on LSST+Euclid. Accepted in MNRAS. Python implementation available at http://blendz.readthedocs.io/

Objective Bayesian analysis of neutrino masses and hierarchy

Authors: Alan F. Heavens, Elena Sellentin

Abstract: Given the precision of current neutrino data, priors still impact noticeably the constraints on neutrino masses and their hierarchy. To avoid our understanding of neutrinos being driven by prior assumptions, we construct a prior that is mathematically minimally informative. Using the constructed uninformative prior, we find that the normal hierarchy is favoured but with inconclusive posterior odds… ▽ More Given the precision of current neutrino data, priors still impact noticeably the constraints on neutrino masses and their hierarchy. To avoid our understanding of neutrinos being driven by prior assumptions, we construct a prior that is mathematically minimally informative. Using the constructed uninformative prior, we find that the normal hierarchy is favoured but with inconclusive posterior odds of 5.1:1. Better data is hence needed before the neutrino masses and their hierarchy can be well constrained. We find that the next decade of cosmological data should provide conclusive evidence if the normal hierarchy with negligible minimum mass is correct, and if the uncertainty in the sum of neutrino masses drops below 0.025 eV. On the other hand, if neutrinos obey the inverted hierarchy, achieving strong evidence will be difficult with the same uncertainties. Our uninformative prior was constructed from principles of the Objective Bayesian approach. The prior is called a reference prior and is minimally informative in the specific sense that the information gain after collection of data is maximised. The prior is computed for the combination of neutrino oscillation data and cosmological data and still applies if the data improve. △ Less

Submitted 6 April, 2018; v1 submitted 26 February, 2018; originally announced February 2018.

Comments: 15 pages. Minor changes to match accepted version in JCAP

On the use of the Edgeworth expansion in cosmology I: how to foresee and evade its pitfalls

Authors: Elena Sellentin, Andrew H. Jaffe, Alan F. Heavens

Abstract: Non-linear gravitational collapse introduces non-Gaussian statistics into the matter fields of the late Universe. As the large-scale structure is the target of current and future observational campaigns, one would ideally like to have the full probability density function of these non-Gaussian fields. The only viable way we see to achieve this analytically, at least approximately and in the near f… ▽ More Non-linear gravitational collapse introduces non-Gaussian statistics into the matter fields of the late Universe. As the large-scale structure is the target of current and future observational campaigns, one would ideally like to have the full probability density function of these non-Gaussian fields. The only viable way we see to achieve this analytically, at least approximately and in the near future, is via the Edgeworth expansion. We hence rederive this expansion for Fourier modes of non-Gaussian fields and then continue by putting it into a wider statistical context than previously done. We show that in its original form, the Edgeworth expansion only works if the non-Gaussian signal is averaged away. This is counterproductive, since we target the parameter-dependent non-Gaussianities as a signal of interest. We hence alter the analysis at the decisive step and now provide a roadmap towards a controlled and unadulterated analysis of non-Gaussianities in structure formation (with the Edgeworth expansion). Our central result is that, although the Edgeworth expansion has pathological properties, these can be predicted and avoided in a careful manner. We also show that, despite the non-Gaussianity coupling all modes, the Edgeworth series may be applied to any desired subset of modes, since this is equivalent (to the level of the approximation) to marginalising over the exlcuded modes. In this first paper of a series, we restrict ourselves to the sampling properties of the Edgeworth expansion, i.e.~how faithfully it reproduces the distribution of non-Gaussian data. A follow-up paper will detail its Bayesian use, when parameters are to be inferred. △ Less

Submitted 11 September, 2017; originally announced September 2017.

Comments: 25 pages, 7 figures

Massive data compression for parameter-dependent covariance matrices

Authors: Alan Heavens, Elena Sellentin, Damien de Mijolla, Alvise Vianello

Abstract: We show how the massive data compression algorithm MOPED can be used to reduce, by orders of magnitude, the number of simulated datasets that are required to estimate the covariance matrix required for the analysis of gaussian-distributed data. This is relevant when the covariance matrix cannot be calculated directly. The compression is especially valuable when the covariance matrix varies with th… ▽ More We show how the massive data compression algorithm MOPED can be used to reduce, by orders of magnitude, the number of simulated datasets that are required to estimate the covariance matrix required for the analysis of gaussian-distributed data. This is relevant when the covariance matrix cannot be calculated directly. The compression is especially valuable when the covariance matrix varies with the model parameters. In this case, it may be prohibitively expensive to run enough simulations to estimate the full covariance matrix throughout the parameter space. This compression may be particularly valuable for the next-generation of weak lensing surveys, such as proposed for Euclid and LSST, for which the number of summary data (such as band power or shear correlation estimates) is very large, $\sim 10^4$, due to the large number of tomographic redshift bins that the data will be divided into. In the pessimistic case where the covariance matrix is estimated separately for all points in an MCMC analysis, this may require an unfeasible $10^9$ simulations. We show here that MOPED can reduce this number by a factor of 1000, or a factor of $\sim 10^6$ if some regularity in the covariance matrix is assumed, reducing the number of simulations required to a manageable $10^3$, making an otherwise intractable analysis feasible. △ Less

Submitted 5 September, 2017; v1 submitted 20 July, 2017; originally announced July 2017.

Comments: 8 pages. Accepted by MNRAS

On the insufficiency of arbitrarily precise covariance matrices: non-Gaussian weak lensing likelihoods

Authors: Elena Sellentin, Alan F. Heavens

Abstract: We investigate whether a Gaussian likelihood, as routinely assumed in the analysis of cosmological data, is supported by simulated survey data. We define test statistics, based on a novel method that first destroys Gaussian correlations in a dataset, and then measures the non-Gaussian correlations that remain. This procedure flags pairs of datapoints which depend on each other in a non-Gaussian fa… ▽ More We investigate whether a Gaussian likelihood, as routinely assumed in the analysis of cosmological data, is supported by simulated survey data. We define test statistics, based on a novel method that first destroys Gaussian correlations in a dataset, and then measures the non-Gaussian correlations that remain. This procedure flags pairs of datapoints which depend on each other in a non-Gaussian fashion, and thereby identifies where the assumption of a Gaussian likelihood breaks down. Using this diagnostic, we find that non-Gaussian correlations in the CFHTLenS cosmic shear correlation functions are significant. With a simple exclusion of the most contaminated datapoints, the posterior for $s_8$ is shifted without broadening, but we find no significant reduction in the tension with $s_8$ derived from Planck Cosmic Microwave Background data. However, we also show that the one-point distributions of the correlation statistics are noticeably skewed, such that sound weak lensing data sets are intrinsically likely to lead to a systematically low lensing amplitude being inferred. The detected non-Gaussianities get larger with increasing angular scale such that for future wide-angle surveys such as Euclid or LSST, with their very small statistical errors, the large-scale modes are expected to be increasingly affected. The shifts in posteriors may then not be negligible and we recommend that these diagnostic tests be run as part of future analyses. △ Less

Submitted 25 September, 2017; v1 submitted 14 July, 2017; originally announced July 2017.

Comments: Replacement to match accepted MNRAS version

Mapping weak lensing distortions in the Kerr metric

Authors: Arianna I. Renzini, Carlo R. Contaldi, Alan Heavens

Abstract: Einstein's theory of General Relativity implies that energy, i.e. matter, curves space-time and thus deforms lightlike geodesics, giving rise to gravitational lensing. This phenomenon is well understood in the case of the Schwarzschild metric, and has been accurately described in the past; however, lensing in the Kerr space-time has received less attention in the literature despite potential pract… ▽ More Einstein's theory of General Relativity implies that energy, i.e. matter, curves space-time and thus deforms lightlike geodesics, giving rise to gravitational lensing. This phenomenon is well understood in the case of the Schwarzschild metric, and has been accurately described in the past; however, lensing in the Kerr space-time has received less attention in the literature despite potential practical observational applications. In particular, lensing in such space is not expressible as the gradient of a scalar potential and as such is a source of curl-like signatures and an asymmetric shear pattern. In this paper, we develop a differentiable lensing map in the Kerr metric, reworking and extending previous approaches. By using standard tools of weak gravitational lensing, we isolate and quantify the distortion that is uniquely induced by the presence of angular momentum in the metric. We apply this framework to the distortion induced by a Kerr-like foreground object on a distribution of background of sources. We verify that the new unique lensing signature is orders of magnitude below current observational bounds for a range of lens configurations. △ Less

Submitted 13 June, 2017; originally announced June 2017.

Comments: 16 pages, 11 figures

Marginal Likelihoods from Monte Carlo Markov Chains

Authors: Alan Heavens, Yabebal Fantaye, Arrykrishna Mootoovaloo, Hans Eggers, Zafiirah Hosenie, Steve Kroon, Elena Sellentin

Abstract: In this paper, we present a method for computing the marginal likelihood, also known as the model likelihood or Bayesian evidence, from Markov Chain Monte Carlo (MCMC), or other sampled posterior distributions. In order to do this, one needs to be able to estimate the density of points in parameter space, and this can be challenging in high numbers of dimensions. Here we present a Bayesian analysi… ▽ More In this paper, we present a method for computing the marginal likelihood, also known as the model likelihood or Bayesian evidence, from Markov Chain Monte Carlo (MCMC), or other sampled posterior distributions. In order to do this, one needs to be able to estimate the density of points in parameter space, and this can be challenging in high numbers of dimensions. Here we present a Bayesian analysis, where we obtain the posterior for the marginal likelihood, using $k$th nearest-neighbour distances in parameter space, using the Mahalanobis distance metric, under the assumption that the points in the chain (thinned if required) are independent. We generalise the algorithm to apply to importance-sampled chains, where each point is assigned a weight. We illustrate this with an idealised posterior of known form with an analytic marginal likelihood, and show that for chains of length $\sim 10^5$ points, the technique is effective for parameter spaces with up to $\sim 20$ dimensions. We also argue that $k=1$ is the optimal choice, and discuss failure modes for the algorithm. In a companion paper (Heavens et al. 2017) we apply the technique to the main MCMC chains from the 2015 Planck analysis of cosmic background radiation data, to infer that quantitatively the simplest 6-parameter flat $Λ$CDM standard model of cosmology is preferred over all extensions considered. △ Less

Submitted 11 April, 2017; originally announced April 2017.

No evidence for extensions to the standard cosmological model

Authors: Alan Heavens, Yabebal Fantaye, Elena Sellentin, Hans Eggers, Zafiirah Hosenie, Steve Kroon, Arrykrishna Mootoovaloo

Abstract: We compute the Bayesian Evidence for models considered in the main analysis of Planck cosmic microwave background data. By utilising carefully-defined nearest-neighbour distances in parameter space, we reuse the Monte Carlo Markov Chains already produced for parameter inference to compute Bayes factors $B$ for many different model-dataset combinations. Standard 6-parameter flat $Λ$CDM model is fav… ▽ More We compute the Bayesian Evidence for models considered in the main analysis of Planck cosmic microwave background data. By utilising carefully-defined nearest-neighbour distances in parameter space, we reuse the Monte Carlo Markov Chains already produced for parameter inference to compute Bayes factors $B$ for many different model-dataset combinations. Standard 6-parameter flat $Λ$CDM model is favoured over all other models considered, with curvature being mildly favoured only when CMB lensing is not included. Many alternative models are strongly disfavoured by the data, including primordial correlated isocurvature models ($\ln B=-7.8$), non-zero scalar-to-tensor ratio ($\ln B=-4.3$), running of the spectral index ($\ln B = -4.7$), curvature ($\ln B=-3.6$), non-standard numbers of neutrinos ($\ln B=-3.1$), non-standard neutrino masses ($\ln B=-3.2$), non-standard lensing potential ($\ln B=-4.6$), evolving dark energy ($\ln B=-3.2$), sterile neutrinos ($\ln B=-6.9$), and extra sterile neutrinos with a non-zero scalar-to-tensor ratio ($\ln B=-10.8$). Other models are less strongly disfavoured with respect to flat $Λ$CDM. As with all analyses based on Bayesian Evidence, the final numbers depend on the widths of the parameter priors. We adopt the priors used in the Planck analysis, while performing a prior sensitivity analysis. Our quantitative conclusion is that extensions beyond the standard cosmological model are disfavoured by Planck data. Only when newer Hubble constant measurements are included does $Λ$CDM become disfavoured, and only mildly, compared with a dynamical dark energy model ($\ln B\sim +2$). △ Less

Submitted 9 August, 2017; v1 submitted 11 April, 2017; originally announced April 2017.

Comments: 5 pages. Accepted for publication in PRL. Effect of inclusion of recent H0 measurements is added

Journal ref: Phys. Rev. Lett. 119, 101301 (2017)

Perturbation Theory for BAO reconstructed fields: one-loop results in real-space matter density field

Authors: Chiaki Hikage, Kazuya Koyama, Alan Heavens

Abstract: We compute the power spectrum at one-loop order in standard perturbation theory for the matter density field to which a standard Lagrangian Baryonic acoustic oscillation (BAO) reconstruction technique is applied. The BAO reconstruction method corrects the bulk motion associated with the gravitational evolution using the inverse Zel'dovich approximation (ZA) for the smoothed density field. We find… ▽ More We compute the power spectrum at one-loop order in standard perturbation theory for the matter density field to which a standard Lagrangian Baryonic acoustic oscillation (BAO) reconstruction technique is applied. The BAO reconstruction method corrects the bulk motion associated with the gravitational evolution using the inverse Zel'dovich approximation (ZA) for the smoothed density field. We find that the overall amplitude of one-loop contributions in the matter power spectrum substantially decrease after reconstruction. The reconstructed power spectrum thereby approaches the initial linear spectrum when the smoothed density field is close enough to linear, i.e., the smoothing scale $R_s$ larger than around 10$h^{-1}$Mpc. On smaller $R_s$,however, the deviation from the linear spectrum becomes significant on large scales ($k\lt R_s^{-1}$) due to the nonlinearity in the smoothed density field, and the reconstruction is inaccurate. Compared with N-body simulations, we show that the reconstructed power spectrum at one loop order agrees with simulations better than the unreconstructed power spectrum. We also calculate the tree-level bispectrum in standard perturbation theory to investigate non-Gaussianity in the reconstructed matter density field. We show that the amplitude of the bispectrum significantly decreases for small $k$ after reconstruction and that the tree-level bispectrum agrees well with N-body results in the weakly nonlinear regime. △ Less

Submitted 21 July, 2017; v1 submitted 22 March, 2017; originally announced March 2017.

Comments: 18 pages, 7 figures, accepted for publications in PRD

Journal ref: Phys. Rev. D 96, 043513 (2017)

The Pan-STARRS1 Surveys

Authors: K. C. Chambers, E. A. Magnier, N. Metcalfe, H. A. Flewelling, M. E. Huber, C. Z. Waters, L. Denneau, P. W. Draper, D. Farrow, D. P. Finkbeiner, C. Holmberg, J. Koppenhoefer, P. A. Price, A. Rest, R. P. Saglia, E. F. Schlafly, S. J. Smartt, W. Sweeney, R. J. Wainscoat, W. S. Burgett, S. Chastel, T. Grav, J. N. Heasley, K. W. Hodapp, R. Jedicke , et al. (101 additional authors not shown)

Abstract: Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3π$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$σ$ point source limiting sensitivities in the stacked 3$π$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across… ▽ More Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3π$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$σ$ point source limiting sensitivities in the stacked 3$π$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ Δra, Δdec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI. △ Less

Submitted 28 January, 2019; v1 submitted 16 December, 2016; originally announced December 2016.

Comments: 38 pages, 29 figures, 12 tables

Unequal-Time Correlators for Cosmology

Authors: T. D. Kitching, A. F. Heavens

Abstract: Measurements of the power spectrum from large-scale structure surveys have to date assumed an equal-time approximation, where the full cross-correlation power spectrum of the matter density field evaluated at different times (or distances) has been approximated either by the power spectrum at a fixed time, or in an improved fashion, by a geometric mean $P(k; r_1, r_2)=[P(k; r_1) P(k; r_2)]^{1/2}$.… ▽ More Measurements of the power spectrum from large-scale structure surveys have to date assumed an equal-time approximation, where the full cross-correlation power spectrum of the matter density field evaluated at different times (or distances) has been approximated either by the power spectrum at a fixed time, or in an improved fashion, by a geometric mean $P(k; r_1, r_2)=[P(k; r_1) P(k; r_2)]^{1/2}$. In this paper we investigate the expected impact of the geometric mean ansatz, and present an application in assessing the impact on weak gravitational lensing cosmological parameter inference, using a perturbative unequal-time correlator. As one might expect, we find that the impact of this assumption is greatest at large separations in redshift $Δz > 0.3$ where the change in the amplitude of the matter power spectrum can be as much as $10$ percent for $k > 5h$Mpc$^{-1}$. However, of more concern is that the corrections for small separations, where the clustering is not close to zero, may not be negligibly small. In particular, we find that for a Euclid- or LSST-like weak lensing experiment the assumption of equal-time correlators may result in biased predictions of the cosmic shear power spectrum, and that the impact is strongly dependent on the amplitude of the intrinsic alignment signal. To compute unequal-time correlations to sufficient accuracy will require advances in either perturbation theory to high $k$-modes, or extensive use of simulations. △ Less

Submitted 2 December, 2016; originally announced December 2016.

Comments: 11 pages, 7 figures. Submitted to PRD

Journal ref: Phys. Rev. D 95, 063522 (2017)

The Limits of Cosmic Shear

Authors: Thomas D. Kitching, Justin Alsing, Alan F. Heavens, Raul Jimenez, Jason D. McEwen, Licia Verde

Abstract: In this paper we discuss the commonly-used limiting cases, or approximations, for two-point cosmic shear statistics. We discuss the most prominent assumptions in this statistic: the flat-sky (small angle limit), the Limber (Bessel-to-delta function limit) and the Hankel transform (large l-mode limit) approximations; that the vast majority of cosmic shear results to date have used simultaneously. W… ▽ More In this paper we discuss the commonly-used limiting cases, or approximations, for two-point cosmic shear statistics. We discuss the most prominent assumptions in this statistic: the flat-sky (small angle limit), the Limber (Bessel-to-delta function limit) and the Hankel transform (large l-mode limit) approximations; that the vast majority of cosmic shear results to date have used simultaneously. We find that the combined effect of these approximations can suppress power by >1% on scales of l<40. A fully non-approximated cosmic shear study should use a spherical-sky, non-Limber-approximated power spectrum analysis; and a transform involving Wigner small-d matrices in place of the Hankel transform. These effects, unaccounted for, would constitute at least 11% of the total budget for systematic effects for a power spectrum analysis of a Euclid-like experiment; but they are unnecessary. △ Less

Submitted 3 May, 2017; v1 submitted 15 November, 2016; originally announced November 2016.

Comments: 15 pages, accepted to MNRAS

Journal ref: Mon Not R Astron Soc (2017) 469 (3): 2737-2749

Early Cosmology Constrained

Authors: Licia Verde, Emilio Bellini, Cassio Pigozzo, Alan F. Heavens, Raul Jimenez

Abstract: We investigate our knowledge of early universe cosmology by exploring how much additional energy density can be placed in different components beyond those in the $Λ$CDM model. To do this we use a method to separate early- and late-universe information enclosed in observational data, thus markedly reducing the model-dependency of the conclusions. We find that the 95\% credibility regions for extra… ▽ More We investigate our knowledge of early universe cosmology by exploring how much additional energy density can be placed in different components beyond those in the $Λ$CDM model. To do this we use a method to separate early- and late-universe information enclosed in observational data, thus markedly reducing the model-dependency of the conclusions. We find that the 95\% credibility regions for extra energy components of the early universe at recombination are: non-accelerating additional fluid density parameter $Ω_{\rm MR} < 0.006$ and extra radiation parameterised as extra effective neutrino species $2.3 < N_{\rm eff} < 3.2$ when imposing flatness. Our constraints thus show that even when analyzing the data in this largely model-independent way, the possibility of hiding extra energy components beyond $Λ$CDM in the early universe is seriously constrained by current observations. We also find that the standard ruler, the sound horizon at radiation drag, can be well determined in a way that does not depend on late-time Universe assumptions, but depends strongly on early-time physics and in particular on additional components that behave like radiation. We find that the standard ruler length determined in this way is $r_{\rm s} = 147.4 \pm 0.7$ Mpc if the radiation and neutrino components are standard, but the uncertainty increases by an order of magnitude when non-standard dark radiation components are allowed, to $r_{\rm s} = 150 \pm 5$ Mpc. △ Less

Submitted 1 November, 2016; originally announced November 2016.

Comments: Submitted to JCAP

Journal ref: Journal of Cosmology and Astroparticle Physics, Issue 04, article id. 023 (2017)

Quantifying lost information due to covariance matrix estimation in parameter inference

Authors: Elena Sellentin, Alan F. Heavens

Abstract: Parameter inference with an estimated covariance matrix systematically loses information due to the remaining uncertainty of the covariance matrix. Here, we quantify this loss of precision and develop a framework to hypothetically restore it, which allows to judge how far away a given analysis is from the ideal case of a known covariance matrix. We point out that it is insufficient to estimate thi… ▽ More Parameter inference with an estimated covariance matrix systematically loses information due to the remaining uncertainty of the covariance matrix. Here, we quantify this loss of precision and develop a framework to hypothetically restore it, which allows to judge how far away a given analysis is from the ideal case of a known covariance matrix. We point out that it is insufficient to estimate this loss by debiasing a Fisher matrix as previously done, due to a fundamental inequality that describes how biases arise in non-linear functions. We therefore develop direct estimators for parameter credibility contours and the figure of merit. We apply our results to DES Science Verification weak lensing data, detecting a 10% loss of information that increases their credibility contours. No significant loss of information is found for KiDS. For a Euclid-like survey, with about 10 nuisance parameters we find that 2900 simulations are sufficient to limit the systematically lost information to 1%, with an additional uncertainty of about 2%. Without any nuisance parameters 1900 simulations are sufficient to only lose 1% of information. We also derive an estimator for the Fisher matrix of the unknown true covariance matrix, two estimators of its inverse with different physical meanings, and an estimator for the optimally achievable figure of merit. The formalism here quantifies the gains to be made by running more simulated datasets, allowing decisions to be made about numbers of simulations in an informed way. △ Less

Submitted 15 March, 2017; v1 submitted 2 September, 2016; originally announced September 2016.

Comments: Figure 6 holds for all surveys; replacement to match published version

The length of the low-redshift standard ruler

Authors: Licia Verde, Jose Luis Bernal, Alan F. Heavens, Raul Jimenez

Abstract: Assuming the existence of standard rulers, standard candles and standard clocks, requiring only the cosmological principle, a metric theory of gravity, a smooth expansion history, and using state-of-the-art observations, we determine the length of the "low-redshift standard ruler". The data we use are a compilation of recent Baryon acoustic oscillation data (relying on the standard ruler), Type 1A… ▽ More Assuming the existence of standard rulers, standard candles and standard clocks, requiring only the cosmological principle, a metric theory of gravity, a smooth expansion history, and using state-of-the-art observations, we determine the length of the "low-redshift standard ruler". The data we use are a compilation of recent Baryon acoustic oscillation data (relying on the standard ruler), Type 1A supernovæ (as standard candles), ages of early type galaxies (as standard clocks) and local determinations of the Hubble constant (as a local anchor of the cosmic distance scale). In a standard $Λ$CDM cosmology the "low-redshift standard ruler" coincides with the sound horizon at radiation drag, which can also be determined --in a model dependent way-- from CMB observations. However, in general, the two quantities need not coincide. We obtain constraints on the length of the low-redshift standard ruler: $r^h_{\rm s}=101.0 \pm 2.3 h^{-1}$ Mpc, when using only Type 1A supernovæ and Baryon acoustic oscillations, and $r_{\rm s}=150.0\pm 4.7 $ Mpc when using clocks to set the Hubble normalisation, while $r_{\rm s}=141.0\pm 5.5 $ Mpc when using the local Hubble constant determination (using both yields $r_{\rm s}=143.9\pm 3.1 $ Mpc). The low-redshift determination of the standard ruler has an error which is competitive with the model-dependent determination from cosmic microwave background measurements made with the {\em Planck} satellite, which assumes it is the sound horizon at the end of baryon drag. △ Less

Submitted 18 July, 2016; originally announced July 2016.

Comments: MCMC output available on request; submitted to MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 467, Issue 1, p.731-736 (2017)

Cosmological parameters, shear maps and power spectra from CFHTLenS using Bayesian hierarchical inference

Authors: Justin Alsing, Alan F. Heavens, Andrew H. Jaffe

Abstract: We apply two Bayesian hierarchical inference schemes to infer shear power spectra, shear maps and cosmological parameters from the CFHTLenS weak lensing survey - the first application of this method to data. In the first approach, we sample the joint posterior distribution of the shear maps and power spectra by Gibbs sampling, with minimal model assumptions. In the second approach, we sample the j… ▽ More We apply two Bayesian hierarchical inference schemes to infer shear power spectra, shear maps and cosmological parameters from the CFHTLenS weak lensing survey - the first application of this method to data. In the first approach, we sample the joint posterior distribution of the shear maps and power spectra by Gibbs sampling, with minimal model assumptions. In the second approach, we sample the joint posterior of the shear maps and cosmological parameters, providing a new, accurate and principled approach to cosmological parameter inference from cosmic shear data. As a first demonstration on data we perform a 2-bin tomographic analysis to constrain cosmological parameters and investigate the possibility of photometric redshift bias in the CFHTLenS data. Under the baseline $Λ$CDM model we constrain $S_8 = σ_8(Ω_\mathrm{m}/0.3)^{0.5} = 0.67 ^{\scriptscriptstyle+ 0.03 }_{\scriptscriptstyle- 0.03 }$ $(68\%)$, consistent with previous CFHTLenS analysis but in tension with Planck. Adding neutrino mass as a free parameter we are able to constrain $\sum m_ν< 4.6\mathrm{eV}$ (95%) using CFHTLenS data alone. Including a linear redshift dependent photo-$z$ bias $Δz = p_2(z - p_1)$, we find $p_1=-0.25 ^{\scriptscriptstyle+ 0.53 }_{\scriptscriptstyle- 0.60 }$ and $p_2 = -0.15 ^{\scriptscriptstyle+ 0.17 }_{\scriptscriptstyle- 0.15 }$, and tension with Planck is only alleviated under very conservative prior assumptions. Neither the non-minimal neutrino mass or photo-$z$ bias models are significantly preferred by the CFHTLenS (2-bin tomography) data. △ Less

Submitted 9 May, 2017; v1 submitted 30 June, 2016; originally announced July 2016.

Comments: Matches accepted version

Generalisations of Fisher Matrices

Authors: Alan Heavens

Abstract: Fisher matrices play an important role in experimental design and in data analysis. Their primary role is to make predictions for the inference of model parameters - both their errors and covariances. In this short review, I outline a number of extensions to the simple Fisher matrix formalism, covering a number of recent developments in the field. These are: (a) situations where the data (in the f… ▽ More Fisher matrices play an important role in experimental design and in data analysis. Their primary role is to make predictions for the inference of model parameters - both their errors and covariances. In this short review, I outline a number of extensions to the simple Fisher matrix formalism, covering a number of recent developments in the field. These are: (a) situations where the data (in the form of (x,y) pairs) have errors in both x and y; (b) modifications to parameter inference in the presence of systematic errors, or through fixing the values of some model parameters; (c) Derivative Approximation for LIkelihoods (DALI) - higher-order expansions of the likelihood surface, going beyond the Gaussian shape approximation; (d) extensions of the Fisher-like formalism, to treat model selection problems with Bayesian evidence. △ Less

Submitted 21 June, 2016; originally announced June 2016.

Comments: Invited review article for Entropy special issue on 'Applications of Fisher Information in Sciences'. Accepted version