Showing 1–18 of 18 results for author: Heydenreich, S

The Construction of Large-scale Structure Catalogs for the Dark Energy Spectroscopic Instrument

Authors: A. J. Ross, J. Aguilar, S. Ahlen, S. Alam, A. Anand, S. Bailey, D. Bianchi, S. Brieden, D. Brooks, E. Burtin, A. Carnero Rosell, E. Chaussidon, T. Claybaugh, S. Cole, K. Dawson, A. de la Macorra, A. de Mattia, Arjun Dey, Biprateep Dey, P. Doel, K. Fanning, S. Ferraro, J. Ereza, A. Font-Ribera, J. E. Forero-Romero , et al. (59 additional authors not shown)

Abstract: We present the technical details on how large-scale structure (LSS) catalogs are constructed from redshifts measured from spectra observed by the Dark Energy Spectroscopic Instrument (DESI). The LSS catalogs provide the information needed to determine the relative number density of DESI tracers as a function of redshift and celestial coordinates and, e.g., determine clustering statistics. We produ… ▽ More We present the technical details on how large-scale structure (LSS) catalogs are constructed from redshifts measured from spectra observed by the Dark Energy Spectroscopic Instrument (DESI). The LSS catalogs provide the information needed to determine the relative number density of DESI tracers as a function of redshift and celestial coordinates and, e.g., determine clustering statistics. We produce catalogs that are weighted subsamples of the observed data, each matched to a weighted `random' catalog that forms an unclustered sampling of the probability density that DESI could have observed those data at each location. Precise knowledge of the DESI observing history and associated hardware performance allows for a determination of the DESI footprint and the number of times DESI has covered it at sub-arcsecond level precision. This enables the completeness of any DESI sample to be modeled at this same resolution. The pipeline developed to create LSS catalogs has been designed to easily allow robustness tests and enable future improvements. We describe how it allows ongoing work improving the match between galaxy and random catalogs, such as including further information when assigning redshifts to randoms, accounting for fluctuations in target density, accounting for variation in the redshift success rate, and accommodating blinding schemes. △ Less

Submitted 26 May, 2024; originally announced May 2024.

Comments: Supporting publication of DESI 2024II: Sample definitions, characteristics, and two-point clustering statistics

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'

KiDS-1000 and DES-Y1 combined: Cosmology from peak count statistics

Authors: Joachim Harnois-Deraps, Sven Heydenreich, Benjamin Giblin, Nicolas Martinet, Tilman Troester, Marika Asgari, Pierre Burger, Tiago Castro, Klaus Dolag, Catherine Heymans, Hendrik Hildebrandt, Benjamin Joachimi, Angus H. Wright

Abstract: We analyse the fourth data release of the Kilo Degree Survey (KiDS-1000) and extract cosmological parameter constraints based on the cosmic shear peak count statistics. Peaks are identified in aperture mass maps in which the filter is maximally sensitive to angular scales in the range 2-4arcmin, probing deep into the non-linear regime of structure formation. We interpret our results with a simulat… ▽ More We analyse the fourth data release of the Kilo Degree Survey (KiDS-1000) and extract cosmological parameter constraints based on the cosmic shear peak count statistics. Peaks are identified in aperture mass maps in which the filter is maximally sensitive to angular scales in the range 2-4arcmin, probing deep into the non-linear regime of structure formation. We interpret our results with a simulation-based inference pipeline, sampling over a broad $w$CDM prior volume and marginalising over uncertainties on shape calibration, photometric redshift distribution, intrinsic alignment and baryonic feedback. Our measurements constrain the structure growth parameter and the amplitude of the non-linear intrinsic alignment model to $Σ_8 \equiv σ_8\left[Ω_{\rm m}/0.3\right]^{0.60}=0.765^{+0.030}_{-0.030}$ and $A_{\rm IA}= 0.71^{+0.42}_{-0.42}$, respectively, in agreement with previous KiDS-1000 results based on two-point shear statistics. These results are robust against modelling of the non-linear physics, different scale cuts and selections of tomographic bins. The posterior is also consistent with that from the Dark Energy Survey Year-1 peak count analysis presented in Harnois-Déraps et al (2021), and hence we jointly analyse both surveys. We obtain $Σ_8^{\rm joint} \equiv σ_8\left[Ω_{\rm m}/0.3\right]^{0.57}=0.759^{+0.020}_{-0.017}$, in agreement with the Planck $w$CDM results. The shear-CMB tension on this parameter increases to $3.1σ$ when forcing $w=-1.0$, and to $4.1σ$ if comparing instead with $S_{8,Λ{\rm CDM}}^{\rm joint} = 0.736^{+0.016}_{-0.018}$, one of the tightest constraints to date on this quantity. (abridged) △ Less

Submitted 16 May, 2024; originally announced May 2024.

Comments: 25 pages, 20 figures. comments welcome

Systematic Effects in Galaxy-Galaxy Lensing with DESI

Authors: J. U. Lange, C. Blake, C. Saulder, N. Jeffrey, J. DeRose, G. Beltz-Mohrmann, N. Emas, C. Garcia-Quintero, B. Hadzhiyska, S. Heydenreich, M. Ishak, S. Joudaki, E. Jullo, A. Krolewski, A. Leauthaud, L. Medina-Varela, A. Porredon, G. Rossi, R. Ruggeri, E. Xhakaj, S. Yuan, J. Aguilar, S. Ahlen, D. Brooks, T. Claybaugh , et al. (34 additional authors not shown)

Abstract: The Dark Energy Spectroscopic Instrument (DESI) survey will measure spectroscopic redshifts for millions of galaxies across roughly $14,000 \, \mathrm{deg}^2$ of the sky. Cross-correlating targets in the DESI survey with complementary imaging surveys allows us to measure and analyze shear distortions caused by gravitational lensing in unprecedented detail. In this work, we analyze a series of mock… ▽ More The Dark Energy Spectroscopic Instrument (DESI) survey will measure spectroscopic redshifts for millions of galaxies across roughly $14,000 \, \mathrm{deg}^2$ of the sky. Cross-correlating targets in the DESI survey with complementary imaging surveys allows us to measure and analyze shear distortions caused by gravitational lensing in unprecedented detail. In this work, we analyze a series of mock catalogs with ray-traced gravitational lensing and increasing sophistication to estimate systematic effects on galaxy-galaxy lensing estimators such as the tangential shear $γ_{\mathrm{t}}$ and the excess surface density $ΔΣ$. We employ mock catalogs tailored to the specific imaging surveys overlapping with the DESI survey: the Dark Energy Survey (DES), the Hyper Suprime-Cam (HSC) survey, and the Kilo-Degree Survey (KiDS). Among others, we find that fiber incompleteness can have significant effects on galaxy-galaxy lensing estimators but can be corrected effectively by up-weighting DESI targets with fibers by the inverse of the fiber assignment probability. Similarly, we show that intrinsic alignment and lens magnification are expected to be statistically significant given the precision forecasted for the DESI year-1 data set. Our study informs several analysis choices for upcoming cross-correlation studies of DESI with DES, HSC, and KiDS. △ Less

Submitted 14 April, 2024; originally announced April 2024.

Comments: 18 pages, 13 figures

Redshift evolution and covariances for joint lensing and clustering studies with DESI Y1

Authors: Sihan Yuan, Chris Blake, Alex Krolewski, Johannes Lange, Jack Elvin-Poole, Alexie Leauthaud, Joseph DeRose, Jessica Nicole Aguilar, Steven Ahlen, Gillian Beltz-Mohrmann, David Brooks, Todd Claybaugh, Axel de la Macorra, Peter Doel, Ni Putu Audita Placida Emas, Simone Ferraro, Jaime E. Forero-Romero, Cristhian Garcia-Quintero, Enrique Gaztañaga, Satya Gontcho A Gontcho, Boryana Hadzhiyska, Sven Heydenreich, Klaus Honscheid, Mustapha Ishak, Shahab Joudaki , et al. (26 additional authors not shown)

Abstract: Galaxy-galaxy lensing (GGL) and clustering measurements from the Dark Energy Spectroscopic Instrument Year 1 (DESI Y1) dataset promise to yield unprecedented combined-probe tests of cosmology and the galaxy-halo connection. In such analyses, it is essential to identify and characterise all relevant statistical and systematic errors. In this paper, we forecast the covariances of DESI Y1 GGL+cluster… ▽ More Galaxy-galaxy lensing (GGL) and clustering measurements from the Dark Energy Spectroscopic Instrument Year 1 (DESI Y1) dataset promise to yield unprecedented combined-probe tests of cosmology and the galaxy-halo connection. In such analyses, it is essential to identify and characterise all relevant statistical and systematic errors. In this paper, we forecast the covariances of DESI Y1 GGL+clustering measurements and characterise the systematic bias due to redshift evolution in the lens samples. Focusing on the projected clustering and galaxy-galaxy lensing correlations, we compute a Gaussian analytical covariance, using a suite of N-body and log-normal simulations to characterise the effect of the survey footprint. Using the DESI One Percent Survey data, we measure the evolution of galaxy bias parameters for the DESI Luminous Red Galaxy (LRG) and Bright Galaxy Survey (BGS) samples. We find mild evolution in the LRGs in 0.4 < z < 0.8, subdominant compared to the expected statistical errors. For BGS, we find less evolution effects for brighter absolute magnitude cuts, at the cost of reduced sample size. We find that with a fiducial redshift bin width delta z = 0.1, evolution effects on GGL is negligible across all scales, all fiducial selection cuts, all fiducial redshift bins, given DESI Y1 sample size. Galaxy clustering is more sensitive to evolution due to the bias squared scaling. Nevertheless the redshift evolution effect is insignificant for clustering above the 1-halo scale of 0.1Mpc/h. For studies that wish to reliably access smaller scales, additional treatment of redshift evolution is likely needed. This study serves as a reference for GGL and clustering studies using the DESI Y1 sample △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 19 pages, 15 figures, submitted to MNRAS

KiDS-1000 cosmology: Combined second- and third-order shear statistics

Authors: Pierre A. Burger, Lucas Porth, Sven Heydenreich, Laila Linke, Niek Wielders, Peter Schneider, Marika Asgari, Tiago Castro, Klaus Dolag, Joachim Harnois-Deraps, Konrad Kuijken, Nicolas Martinet

Abstract: This paper performs the first cosmological parameter analysis of the KiDS-1000 data with second- and third-order shear statistics. This work builds on a series of papers that describe the roadmap to third-order shear statistics. We derive and test a combined model of the second-order shear statistic, namely the COSEBIs and the third-order aperture mass statistics $\langle M_\mathrm{ap}^3\rangle$ i… ▽ More This paper performs the first cosmological parameter analysis of the KiDS-1000 data with second- and third-order shear statistics. This work builds on a series of papers that describe the roadmap to third-order shear statistics. We derive and test a combined model of the second-order shear statistic, namely the COSEBIs and the third-order aperture mass statistics $\langle M_\mathrm{ap}^3\rangle$ in a tomographic set-up. We validate our pipeline with $N$-body simulations that mock the fourth Kilo Degree survey data release. To model the second- and third-order statistics, we use the latest version of \textsc{HMcode2020} for the power spectrum and \textsc{BiHalofit} for the bispectrum. Furthermore, we use an analytic description to model intrinsic alignments and hydro-dynamical simulations to model the effect of baryonic feedback processes. Lastly, we decreased the dimension of the data vector significantly by considering for the $\langle M_\mathrm{ap}^3\rangle$ part of the data vector only equal smoothing radii, making a data analysis of the fourth Kilo Degree survey data release using a combined analysis of COSEBIs third-order shear statistic possible. We first validate the accuracy of our modelling by analysing a noise-free mock data vector assuming the KiDS-1000 error budget, finding a shift in the maximum-a-posterior of the matter density parameter $ΔΩ_m< 0.02\, σ_{Ω_m}$ and of the structure growth parameter $ΔS_8 < 0.05\, σ_{S_8}$. Lastly, we performed the first KiDS-1000 cosmological analysis using a combined analysis of second- and third-order shear statistics, where we constrained $Ω_m=0.248^{+0.062}_{-0.055}$ and $S_8=σ_8\sqrt{Ω_m/0.3}=0.772\pm0.022$. The geometric average on the errors of $Ω_\mathrm{m}$ and $S_8$ of the combined statistics increased compared to the second-order statistic by 2.2. △ Less

Submitted 25 March, 2024; v1 submitted 15 September, 2023; originally announced September 2023.

Comments: 19 pages, 15 figures

Journal ref: A&A, 683, A103 (2024)

A roadmap to cosmological parameter analysis with third-order shear statistics III: Efficient estimation of third-order shear correlation functions and an application to the KiDS-1000 data

Authors: Lucas Porth, Sven Heydenreich, Pierre Burger, Laila Linke, Peter Schneider

Abstract: Third-order lensing statistics contain a wealth of cosmological information that is not captured by second-order statistics. However, the computational effort for estimating such statistics on forthcoming stage IV surveys is prohibitively expensive. We derive and validate an efficient estimation procedure for the three-point correlation function (3PCF) of polar fields such as weak lensing shear. W… ▽ More Third-order lensing statistics contain a wealth of cosmological information that is not captured by second-order statistics. However, the computational effort for estimating such statistics on forthcoming stage IV surveys is prohibitively expensive. We derive and validate an efficient estimation procedure for the three-point correlation function (3PCF) of polar fields such as weak lensing shear. We then use our approach to measure the shear 3PCF and the third-order aperture mass statistics on the KiDS-1000 survey. We construct an efficient estimator for third-order shear statistics which builds on the multipole decomposition of the 3PCF. We then validate our estimator on mock ellipticity catalogs obtained from $N$-body simulations. Finally, we apply our estimator to the KiDS-1000 data and present a measurement of the third-order aperture statistics in a tomographic setup. Our estimator provides a speedup of a factor of $\sim$ 100-1000 compared to the state-of-the-art estimation procedures. It is also able to provide accurate measurements for squeezed and folded triangle configurations without additional computational effort. We report a significant detection of the tomographic third-order aperture mass statistics in the KiDS-1000 data $(\mathrm{S/N}=6.69)$. Our estimator will make it computationally feasible to measure third-order shear statistics in forthcoming stage IV surveys. Furthermore, it can be used to construct empirical covariance matrices for such statistics. △ Less

Submitted 15 September, 2023; originally announced September 2023.

Comments: 15 pages, 8 figures + appendices

Cluster Cosmology Without Cluster Finding

Authors: Enia Xhakaj, Alexie Leauthaud, Johannes Lange, Elisabeth Krause, Andrew Hearin, Song Huang, Risa H. Wechsler, Sven Heydenreich

Abstract: We propose that observations of super-massive galaxies contain cosmological constraining power similar to conventional cluster cosmology, and we provide promising indications that the associated systematic errors are comparably easier to control. We consider a fiducial spectroscopic and stellar mass complete sample of galaxies drawn from the Dark Energy Spectroscopic Survey (DESI) and forecast how… ▽ More We propose that observations of super-massive galaxies contain cosmological constraining power similar to conventional cluster cosmology, and we provide promising indications that the associated systematic errors are comparably easier to control. We consider a fiducial spectroscopic and stellar mass complete sample of galaxies drawn from the Dark Energy Spectroscopic Survey (DESI) and forecast how constraints on Omega_m-sigma_8 from this sample will compare with those from number counts of clusters based on richness. At fixed number density, we find that massive galaxies offer similar constraints to galaxy clusters. However, a mass-complete galaxy sample from DESI has the potential to probe lower halo masses than standard optical cluster samples (which are typically limited to richness above 20 and halo mass above 10^13.5); additionally, it is straightforward to cleanly measure projected galaxy clustering for such a DESI sample, which we show can substantially improve the constraining power on Omega_m. We also compare the constraining power of stellar mass-limited samples to those from larger but mass-incomplete samples (e.g., the DESI Bright Galaxy Survey, BGS, Sample); relative to a lower number density stellar mass-limited samples, we find that a BGS-like sample improves statistical constraints by 60% for Omega_m and 40% for sigma_8, but this uses small scale information which will be harder to model for BGS. Our initial assessment of the systematics associated with supermassive galaxy cosmology yields promising results. The proposed samples have a 10% satellite fraction, but we show that cosmological constraints may be robust to the impact of satellites. These findings motivate future work to realize the potential of super-massive galaxies to probe lower halo masses than richness-based clusters and to avoid persistent systematics associated with optical cluster finding. △ Less

Submitted 26 May, 2023; originally announced June 2023.

What is the super-sample covariance? A fresh perspective for second-order shear statistics

Authors: Laila Linke, Pierre A. Burger, Sven Heydenreich, Lucas Porth, Peter Schneider

Abstract: Cosmological analyses of second-order weak lensing statistics require precise and accurate covariance estimates. These covariances are impacted by two sometimes neglected terms: A negative contribution to the Gaussian covariance due to finite survey area and the super-sample covariance (SSC) which for the power spectrum contains the impact by Fourier modes larger than the survey window. We show he… ▽ More Cosmological analyses of second-order weak lensing statistics require precise and accurate covariance estimates. These covariances are impacted by two sometimes neglected terms: A negative contribution to the Gaussian covariance due to finite survey area and the super-sample covariance (SSC) which for the power spectrum contains the impact by Fourier modes larger than the survey window. We show here that these two effects are connected and can be seen as correction terms to the "large-field-approximation", the asymptotic case of an infinitely large survey area. We describe the two terms collectively as "Finite-Field-Terms". We derive the covariance of second-order shear statistics from first principles. For this, we use an estimator in real space without relying on an estimator for the power spectrum. The resulting covariance does not scale inversely with the survey area, as naively assumed. This scaling is only correct under the large-field approximation when the contribution of the finite-field terms tends to zero. Furthermore, all parts of the covariance, not only the SSC, depend on the power- and trispectrum at all modes, including those larger than the survey. We also show that it is generally impossible to transform an estimate for the power spectrum covariance into the covariance of a real-space statistic. Such a transformation is only possible in the asymptotic case of the "large-field approximation". Additionally, we find that the total covariance of a real-space statistic can be calculated using correlation functions estimates on spatial scales smaller than the survey window. Consequently, estimating covariances of real-space statistics, in principle, does not require information on spatial scales larger than the survey area. We demonstrate that this covariance estimation method is equivalent to the standard sample covariance method. △ Less

Submitted 1 December, 2023; v1 submitted 23 February, 2023; originally announced February 2023.

Comments: 8 pages + appendix, 3 figures, forthcoming in A&A, replaced by version accepted by journal

Euclid Preparation. XXVIII. Forecasts for ten different higher-order weak lensing statistics

Authors: Euclid Collaboration, V. Ajani, M. Baldi, A. Barthelemy, A. Boyle, P. Burger, V. F. Cardone, S. Cheng, S. Codis, C. Giocoli, J. Harnois-Déraps, S. Heydenreich, V. Kansal, M. Kilbinger, L. Linke, C. Llinares, N. Martinet, C. Parroni, A. Peel, S. Pires, L. Porth, I. Tereno, C. Uhlemann, M. Vicinanza, S. Vinciguerra , et al. (189 additional authors not shown)

Abstract: Recent cosmic shear studies have shown that higher-order statistics (HOS) developed by independent teams now outperform standard two-point estimators in terms of statistical precision thanks to their sensitivity to the non-Gaussian features of large-scale structure. The aim of the Higher-Order Weak Lensing Statistics (HOWLS) project is to assess, compare, and combine the constraining power of ten… ▽ More Recent cosmic shear studies have shown that higher-order statistics (HOS) developed by independent teams now outperform standard two-point estimators in terms of statistical precision thanks to their sensitivity to the non-Gaussian features of large-scale structure. The aim of the Higher-Order Weak Lensing Statistics (HOWLS) project is to assess, compare, and combine the constraining power of ten different HOS on a common set of $Euclid$-like mocks, derived from N-body simulations. In this first paper of the HOWLS series, we computed the nontomographic ($Ω_{\rm m}$, $σ_8$) Fisher information for the one-point probability distribution function, peak counts, Minkowski functionals, Betti numbers, persistent homology Betti numbers and heatmap, and scattering transform coefficients, and we compare them to the shear and convergence two-point correlation functions in the absence of any systematic bias. We also include forecasts for three implementations of higher-order moments, but these cannot be robustly interpreted as the Gaussian likelihood assumption breaks down for these statistics. Taken individually, we find that each HOS outperforms the two-point statistics by a factor of around two in the precision of the forecasts with some variations across statistics and cosmological parameters. When combining all the HOS, this increases to a $4.5$ times improvement, highlighting the immense potential of HOS for cosmic shear cosmological analyses with $Euclid$. The data used in this analysis are publicly released with the paper. △ Less

Submitted 10 July, 2023; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: 33 pages, 24 figures, main results in Fig. 19 & Table 5, version published in A&A

Journal ref: A&A 675, A120 (2023)

A roadmap to cosmological parameter analysis with third-order shear statistics II: Analytic covariance estimate

Authors: Laila Linke, Sven Heydenreich, Pierre A. Burger, Peter Schneider

Abstract: Third-order weak lensing statistics are a promising tool for cosmological analyses since they extract cosmological information in the non-Gaussianity of the cosmic large-scale structure. However, such analyses require precise and accurate models for the covariance. In this second paper of a series on third-order weak lensing statistics, we derive and validate an analytic model for the covariance o… ▽ More Third-order weak lensing statistics are a promising tool for cosmological analyses since they extract cosmological information in the non-Gaussianity of the cosmic large-scale structure. However, such analyses require precise and accurate models for the covariance. In this second paper of a series on third-order weak lensing statistics, we derive and validate an analytic model for the covariance of the third-order aperture statistics $\langle M_\mathrm{ap}^3\rangle$. We derive the covariance model from a real-space estimator for $\langle M_\mathrm{ap}^3\rangle$. We validate the model by comparing it to estimates from simulated Gaussian random fields (GRF) and two sets of N-body simulations. Finally, we perform mock cosmological analyses with the model covariance and the simulation estimate to compare the resulting parameter constraints. We find good agreement between the model and the simulations, both for the GRF and the $N$-body simulations. The figure-of-merit in the $S_8$-$Ω_\mathrm{m}$ plane from our covariance model is within 3\% of the one obtained from the simulated covariances. We also show that our model, which is based on an estimator using convergence maps, can be used to obtain upper and lower bounds for the covariance of an estimator based on three-point shear correlation functions. This second estimator is required for realistic survey data. In our derivation, we find that the $\langle M_\mathrm{ap}^3\rangle$ covariance cannot be obtained from the bispectrum covariance and that it includes several `finite-field terms' that do not scale with the inverse survey area. Our covariance model is sufficiently accurate for analysing stage III surveys. Covariances for statistics in Fourier space cannot always be straightforwardly converted into covariance for real-space statistics. The modelling code is available at https://github.com/sheydenreich/threepoint/releases/ . △ Less

Submitted 16 February, 2023; v1 submitted 8 December, 2022; originally announced December 2022.

Comments: 21 pages + appendix, 16 figures, 1 table; replaced by version accepted to Astronomy & Astrophysics after considering Referee comments

Journal ref: A&A 672, A185 (2023)

A roadmap to cosmological parameter analysis with third-order shear statistics I: Modelling and validation

Authors: Sven Heydenreich, Laila Linke, Pierre Burger, Peter Schneider

Abstract: In this work, which is the first of a series to prepare a cosmological parameter analysis with third-order cosmic shear statistics, we model both the shear three-point correlation functions $Γ^{(i)}$ and the third-order aperture statistics $\langle\mathcal{M}_\mathrm{ap}^3\rangle$ from the BiHalofit bispectrum model and validate these statistics with a series of N-body simulations. We then inves… ▽ More In this work, which is the first of a series to prepare a cosmological parameter analysis with third-order cosmic shear statistics, we model both the shear three-point correlation functions $Γ^{(i)}$ and the third-order aperture statistics $\langle\mathcal{M}_\mathrm{ap}^3\rangle$ from the BiHalofit bispectrum model and validate these statistics with a series of N-body simulations. We then investigate how to bin the shear three-point correlation functions to achieve an unbiased estimate for third-order aperture statistics in real data. Finally, we perform a cosmological parameter analysis on KiDS1000-like mock data with second- and third-order statistics. We recover all cosmological parameters with very little bias. Furthermore, we find that a joint analysis almost doubles the constraining power on $S_8$ and increases the figure-of-merit in the $Ω_\mathrm{m}$-$σ_8$ plane by a factor of 5.9 with respect to an analysis with only second-order shear statistics. Our modelling pipeline is publicly available at https://github.com/sheydenreich/threepoint/releases/. △ Less

Submitted 27 August, 2022; v1 submitted 24 August, 2022; originally announced August 2022.

Comments: updated acknowledgements

Persistent homology in cosmic shear II: A tomographic analysis of DES-Y1

Authors: Sven Heydenreich, Benjamin Brück, Pierre Burger, Joachim Harnois-Déraps, Sandra Unruh, Tiago Castro, Klaus Dolag, Nicolas Martinet

Abstract: We demonstrate how to use persistent homology for cosmological parameter inference in a tomographic cosmic shear survey. We obtain the first cosmological parameter constraints from persistent homology by applying our method to the first-year data of the Dark Energy Survey. To obtain these constraints, we analyse the topological structure of the matter distribution by extracting persistence diagr… ▽ More We demonstrate how to use persistent homology for cosmological parameter inference in a tomographic cosmic shear survey. We obtain the first cosmological parameter constraints from persistent homology by applying our method to the first-year data of the Dark Energy Survey. To obtain these constraints, we analyse the topological structure of the matter distribution by extracting persistence diagrams from signal-to-noise maps of aperture masses. This presents a natural extension to the widely used peak count statistics. Extracting the persistence diagrams from the cosmo-SLICS, a suite of $N$-body simulations with variable cosmological parameters, we interpolate the signal using Gaussian Processes and marginalise over the most relevant systematic effects, including intrinsic alignments and baryonic effects. We find for the structure growth parameter $S_8=0.747^{+0.025}_{-0.031}$, which is in full agreement with other late-time probes. We also constrain the intrinsic alignment parameter to $A=1.54\pm 0.52$, ruling out the case of no intrinsic alignments at a $3σ$-level. △ Less

Submitted 25 April, 2022; originally announced April 2022.

Journal ref: A&A 667, A125 (2022)

Abundance of LIGO/Virgo Black Holes from Microlensing Observations of Quasars with Reverberation Mapping Size Estimates

Authors: A. Esteban-Gutiérrez, E. Mediavilla, J. Jiménez-Vicente, N. Agües-Paszkowsky, J. A. Muñoz, S. Heydenreich

Abstract: Assuming a population of Black Holes (BHs) with masses in the range inferred by LIGO/Virgo from BH mergers, we use quasar microlensing observations to estimate their abundances. We consider a mixed population of stars and BHs and the presence of a smooth dark matter component. We adopt reverberation mapping estimates of the quasar size. According to a Bayesian analysis of the measured microlensing… ▽ More Assuming a population of Black Holes (BHs) with masses in the range inferred by LIGO/Virgo from BH mergers, we use quasar microlensing observations to estimate their abundances. We consider a mixed population of stars and BHs and the presence of a smooth dark matter component. We adopt reverberation mapping estimates of the quasar size. According to a Bayesian analysis of the measured microlensing magnifications, a population of BHs with masses $\sim$ 30$M_{\odot}$ constitutes less than 0.4 % of the total matter at 68 % confidence level (less than 0.9 % at 90 % confidence). We have explored the whole mass range of LIGO/Virgo BHs finding that this upper limit ranges from 0.5 % to 0.4 % at 68 % C.L. (from 1.1 % to 0.9 % at 90 % C.L.) when the BHs mass change from 10 to 60$M_{\odot}$. We estimate a 16 % contribution from the stars, in agreement with previous studies based on a single mass population that do not consider explicitly the presence of BHs. These results are consistent with the estimates of BH abundances from the statistics of LIGO/Virgo mergers and rule out that PBHs (or any type of compact objects), in this mass range constitute a significant fraction of the dark matter. △ Less

Submitted 25 April, 2022; v1 submitted 8 March, 2022; originally announced March 2022.

Comments: 14 pages, 5 figures, 1 table

Limiting the Abundance of LIGO/Virgo Black Holes with Microlensing Observations of Quasars of Finite Size

Authors: A. Esteban-Gutiérrez, E. Mediavilla, J. Jiménez-Vicente, N. Agües-Paszkowsky, J. A. Muñoz, S. Heydenreich

Abstract: We present a simple but general argument that strongly limits the abundance of Primordial Black Holes (PBHs) (or other unknown population of compact objects) with masses similar to those determined by LIGO/Virgo from BH binary mergers. We show that quasar microlensing can be very sensitive to the mass of the lenses, and that it is able to distinguish between stars and BHs of high mass, when the fi… ▽ More We present a simple but general argument that strongly limits the abundance of Primordial Black Holes (PBHs) (or other unknown population of compact objects) with masses similar to those determined by LIGO/Virgo from BH binary mergers. We show that quasar microlensing can be very sensitive to the mass of the lenses, and that it is able to distinguish between stars and BHs of high mass, when the finite size of the source is taken into account. A significant presence of massive BHs would produce frequent high flux magnifications (except for unrealistically large sources) which have been very rarely observed. On the contrary, a typical stellar population would induce flux magnifications consistent with the observations. This result excludes PBHs (or any type of compact object) in the mass range determined by LIGO/Virgo as the main dark matter constituents in the lens galaxies. △ Less

Submitted 8 March, 2022; originally announced March 2022.

Comments: 5 pages, 1 figure

The Impact of the Mass Spectrum of Lenses in Quasar Microlensing Studies. Constraints on a Mixed Population of Primordial Black Holes and Stars

Authors: A. Esteban-Gutiérrez, N. Agües-Paszkowsky, E. Mediavilla, J. Jiménez-Vicente, J. A. Muñoz, S. Heydenreich

Abstract: We show that quasar microlensing magnification statistics induced by a population of point microlenses distributed according to a mass-spectrum can be very well approximated by that of a single-mass, "monochromatic", population. When the spatial resolution (physically defined by the source size) is small as compared with the Einstein radius, the mass of the monochromatic population matches the geo… ▽ More We show that quasar microlensing magnification statistics induced by a population of point microlenses distributed according to a mass-spectrum can be very well approximated by that of a single-mass, "monochromatic", population. When the spatial resolution (physically defined by the source size) is small as compared with the Einstein radius, the mass of the monochromatic population matches the geometric mean of the mass-spectrum. Otherwise, the best-fit mass can be larger. Taking into account the degeneracy with the geometric mean, the interpretation of quasar microlensing observations under the hypothesis of a mixed population of primordial black holes and stars, makes the existence of a significant population of intermediate mass black holes ($\sim$ 100$M_\odot$) unlikely but allows, within a two-$σ$ confidence interval, the presence of a large population ($\gtrsim 40\%$ of the total mass) of substellar black holes ($\sim$ 0.01$M_\odot$). △ Less

Submitted 11 November, 2020; originally announced November 2020.

Comments: 24 pages, 4 figures

Persistent homology in cosmic shear: constraining parameters with topological data analysis

Authors: Sven Heydenreich, Benjamin Brück, Joachim Harnois-Déraps

Abstract: In recent years, cosmic shear has emerged as a powerful tool to study the statistical distribution of matter in our Universe. Apart from the standard two-point correlation functions, several alternative methods like peak count statistics offer competitive results. Here we show that persistent homology, a tool from topological data analysis, can extract more cosmological information than previous m… ▽ More In recent years, cosmic shear has emerged as a powerful tool to study the statistical distribution of matter in our Universe. Apart from the standard two-point correlation functions, several alternative methods like peak count statistics offer competitive results. Here we show that persistent homology, a tool from topological data analysis, can extract more cosmological information than previous methods from the same dataset. For this, we use persistent Betti numbers to efficiently summarise the full topological structure of weak lensing aperture mass maps. This method can be seen as an extension of the peak count statistics, in which we additionally capture information about the environment surrounding the maxima. We first demonstrate the performance in a mock analysis of the KiDS+VIKING-450 data: we extract the Betti functions from a suite of $w$CDM $N$-body simulations and use these to train a Gaussian process emulator that provides rapid model predictions; we next run a Markov-Chain Monte Carlo analysis on independent mock data to infer the cosmological parameters and their uncertainty. When comparing our results, we recover the input cosmology and achieve a constraining power on $S_8 \equiv σ_8\sqrt{Ω_\mathrm{m}/0.3}$ that is 5% tighter than that of peak count statistics. Performing the same analysis on 100 deg$^2$ of Euclid-like simulations, we are able to improve the constraints on $S_8$ and $Ω_\mathrm{m}$ by 18% and 10%, respectively, while breaking some of the degeneracy between $S_8$ and the dark energy equation of state. To our knowledge, the methods presented here are the most powerful topological tools to constrain cosmological parameters with lensing data. △ Less

Submitted 27 July, 2020; originally announced July 2020.

Report number: CPH-GEOTOP-DNRF151

The effects of varying depth in cosmic shear surveys

Authors: Sven Heydenreich, Peter Schneider, Hendrik Hildebrandt, Marika Asgari, Catherine Heymans, Benjamin Joachimi, Konrad Kuijken, Chieh-An Lin, Tilman Tröster, Jan Luca van den Busch

Abstract: We present a semi-analytic model for the shear two-point correlation function of a cosmic shear survey with non-uniform depth. Ground-based surveys are subject to depth variations that primarily arise through varying atmospheric conditions. For a survey like the Kilo-Degree Survey (KiDS), we find that the measured depth variation increases the amplitude of the observed shear correlation function a… ▽ More We present a semi-analytic model for the shear two-point correlation function of a cosmic shear survey with non-uniform depth. Ground-based surveys are subject to depth variations that primarily arise through varying atmospheric conditions. For a survey like the Kilo-Degree Survey (KiDS), we find that the measured depth variation increases the amplitude of the observed shear correlation function at the level of a few percent out to degree-scales, relative to the assumed uniform-depth case. The impact on the inferred cosmological parameters is shown to be insignificant for a KiDS-like survey. For next-generation cosmic shear experiments, however, we conclude that variable depth should be accounted for. △ Less

Submitted 16 December, 2019; v1 submitted 24 October, 2019; originally announced October 2019.

Journal ref: A&A 634, A104 (2020)