-
Cosmological parameters from the joint analysis of Density Split and Second Order Statistics: an Emulator based on the Halo Occupation Distribution
Authors:
Pierre A. Burger,
Enrique Paillas,
Michael J. Hudson
Abstract:
In this work, we develop a simulation-based model to predict the density split (DSS) and second-order shear and clustering statistics. A simulation-based model has the potential to model highly non-linear scales where current DSS models fail. To build this model, we use the \texttt{AbacusSummit} N-body simulation suite from which we measure all necessary statistics and train an emulator based on \…
▽ More
In this work, we develop a simulation-based model to predict the density split (DSS) and second-order shear and clustering statistics. A simulation-based model has the potential to model highly non-linear scales where current DSS models fail. To build this model, we use the \texttt{AbacusSummit} N-body simulation suite from which we measure all necessary statistics and train an emulator based on \texttt{CosmoPower}. In that context, we discuss possible improvements for future emulators to make the measurement less noisy and biased, resulting in more accurate and precise model predictions.
Regarding the emulator's accuracy, we find that the most important aspect is the average of the summary statistics over multiple-shot noise realizations of the foreground galaxies. However, these results probably depend on the chosen number density of the foreground galaxies. Regarding the parameter forecast based on preliminary LOWZxUNIONS data, we find that DSS has more constraining power to derive cosmological parameters and that the joint analysis with second-order statistics is particularly useful for extracting parameters of the galaxy-halo connection.
△ Less
Submitted 7 June, 2024;
originally announced June 2024.
-
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.
-
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.
-
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.
-
KiDS-1000 Cosmology: Constraints from density split statistics
Authors:
Pierre A. Burger,
Oliver Friedrich,
Joachim Harnois-Déraps,
Peter Schneider,
Marika Asgari,
Maciej Bilicki,
Hendrik Hildebrandt,
Angus H. Wright,
Tiago Castro,
Klaus Dolag,
Catherine Heymans,
Benjamin Joachimi,
Konrad Kuijken,
Nicolas Martinet,
HuanYuan Shan,
Tilman Tröster
Abstract:
Context. Weak lensing and clustering statistics beyond two-point functions can capture non-Gaussian information about the matter density field, thereby improving the constraints on cosmological parameters relative to the mainstream methods based on correlation functions and power spectra. Aims. This paper presents a cosmological analysis of the fourth data release of the Kilo Degree Survey based o…
▽ More
Context. Weak lensing and clustering statistics beyond two-point functions can capture non-Gaussian information about the matter density field, thereby improving the constraints on cosmological parameters relative to the mainstream methods based on correlation functions and power spectra. Aims. This paper presents a cosmological analysis of the fourth data release of the Kilo Degree Survey based on the density split statistics, which measures the mean shear profiles around regions classified according to foreground densities. The latter is constructed from a bright galaxy sample, which we further split into red and blue samples, allowing us to probe their respective connection to the underlying dark matter density. Methods. We use the state-of-the-art model of the density splitting statistics and validate its robustness against mock data infused with known systematic effects such as intrinsic galaxy alignment and baryonic feedback. Results. After marginalising over the photometric redshift uncertainty and the residual shear calibration bias, we measure for the full KiDS-bright sample a structure growth parameter of $S_8 = σ_8 \sqrt{Ω_\mathrm{m}/0.3} = 0.74^{+0.03}_{-0.02}$ that is competitive to and consistent with two-point cosmic shear results, a matter density of $Ω_\mathrm{m} = 0.28 \pm 0.02$, and a constant galaxy bias of $b = 1.32^{+0.12}_{-0.10}$.
△ Less
Submitted 10 January, 2023; v1 submitted 3 August, 2022;
originally announced August 2022.