Search

Modern cosmological research in large scale structure has witnessed an increasing number of applications of machine learning methods. Among them, Convolutional Neural Networks (CNNs) have received substantial attention due to their outstanding performance in image classification, cosmological parameter inference and various other tasks. However, many models which make use of CNNs are criticized as "black boxes" due to the difficulties in relating their outputs intuitively and quantitatively to the cosmological fields under investigation. To overcome this challenge, we present the Cosmological Correlator Convolutional Neural Network (C3NN) -- a fusion of CNN architecture with the framework of cosmological N-point correlation functions (NPCFs). We demonstrate that the output of this model can be expressed explicitly in terms of the analytically tractable NPCFs. Together with other auxiliary algorithms, we are able to open the "black box" by quantitatively ranking different orders of the interpretable convolution outputs based on their contribution to classification tasks. As a proof of concept, we demonstrate this by applying our framework to a series of binary classification tasks using Gaussian and Log-normal random fields and relating its outputs to the analytical NPCFs describing the two fields. Furthermore, we exhibit the model's ability to distinguish different dark energy scenarios ($w_0=-0.95$ and $-1.05$) using N-body simulated weak lensing convergence maps and discuss the physical implications coming from their interpretability. With these tests, we show that C3NN combines advanced aspects of machine learning architectures with the framework of cosmological NPCFs, thereby making it an exciting tool with the potential to extract physical insights in a robust and explainable way from observational data., Comment: 19 pages, 8 figures, 5 tables; Comments are welcome!

The last few years have seen the development of a promising theoretical framework for statistics of the cosmic large-scale structure -- the theory of large deviations (LDT) for modelling weak-lensing one-point statistics in the mildly non-linear regime. The goal of this series of papers is to make the leap and lay out the steps to perform an actual data analysis with this theoretical tool. Building upon the LDT framework, in this work (Paper I) we demonstrate how to accurately model the Probability Distribution Function (PDF) of a reconstructed Kaiser-Squires convergence field under a realistic mask, that of the third data release of the Dark Energy Survey (DES). We also present how weak lensing systematics and higher-order lensing corrections due to intrinsic alignments, shear biases, photo-$z$ errors and baryonic feedback can be incorporated in the modelling of the reconstructed convergence PDF. In an upcoming work (Paper II) we will then demonstrate the robustness of our modelling through simulated likelihood analyses, the final step required before applying our method to actual data., Comment: 39 pages, 9 figures and 2 tables. Accepted for publication in JCAP, comments welcome

We present the integrated 3-point correlation functions (3PCF) involving both the cosmic shear and the galaxy density fields. These are a set of higher-order statistics that describe the modulation of local 2-point correlation functions (2PCF) by large-scale features in the fields, and which are easy to measure from galaxy imaging surveys. Based on previous works on the shear-only integrated 3PCF, we develop the theoretical framework for modelling 5 new statistics involving the galaxy field and its cross-correlations with cosmic shear. Using realistic galaxy and cosmic shear mocks from simulations, we determine the regime of validity of our models based on leading-order standard perturbation theory with an MCMC analysis that recovers unbiased constraints of the amplitude of fluctuations parameter $A_s$ and the linear and quadratic galaxy bias parameters $b_1$ and $b_2$. Using Fisher matrix forecasts for a DES-Y3-like survey, relative to baseline analyses with conventional 3$\times$2PCFs, we find that the addition of the shear-only integrated 3PCF can improve cosmological parameter constraints by $20-40\%$. The subsequent addition of the new statistics introduced in this paper can lead to further improvements of $10-20\%$, even when utilizing only conservatively large scales where the tree-level models are valid. Our results motivate future work on the galaxy and shear integrated 3PCFs, which offer a practical way to extend standard analyses based on 3$\times$2PCFs to systematically probe the non-Gaussian information content of cosmic density fields., Comment: 19 pages, 8 figures + appendix. Comments are welcome!

The integrated shear 3-point correlation function $\zeta_{\pm}$ measures the correlation between the local shear 2-point function $\xi_{\pm}$ and the 1-point shear aperture mass in patches of the sky. Unlike other higher-order statistics, $\zeta_{\pm}$ can be efficiently measured from cosmic shear data, and it admits accurate theory predictions on a wide range of scales as a function of cosmological and baryonic feedback parameters. Here, we develop and test a likelihood analysis pipeline for cosmological constraints using $\zeta_{\pm}$. We incorporate treatment of systematic effects from photometric redshift uncertainties, shear calibration bias and galaxy intrinsic alignments. We also develop an accurate neural-network emulator for fast theory predictions in MCMC parameter inference analyses. We test our pipeline using realistic cosmic shear maps based on $N$-body simulations with a DES Y3-like footprint, mask and source tomographic bins, finding unbiased parameter constraints. Relative to $\xi_{\pm}$-only, adding $\zeta_{\pm}$ can lead to $\approx 10-25\%$ improvements on the constraints of parameters like $A_s$ (or $\sigma_8$) and $w_0$. We find no evidence in $\xi_{\pm} + \zeta_{\pm}$ constraints of a significant mitigation of the impact of systematics. We also investigate the impact of the size of the apertures where $\zeta_{\pm}$ is measured, and of the strategy to estimate the covariance matrix ($N$-body vs. lognormal). Our analysis solidifies the strong potential of the $\zeta_{\pm}$ statistic and puts forward a pipeline that can be readily used to improve cosmological constraints using real cosmic shear data., Comment: 21 pages, 11 figures, 3 tables. Comments welcome

The integrated shear 3-point correlation function $\zeta_{\pm}$ is a higher-order statistic of the cosmic shear field that describes the modulation of the 2-point correlation function $\xi_{\pm}$ by long-wavelength features in the field. Here, we introduce a new theoretical model to calculate $\zeta_{\pm}$ that is accurate on small angular scales, and that allows to take baryonic feedback effects into account. Our model builds on the realization that the small-scale $\zeta_{\pm}$ is dominated by the nonlinear matter bispectrum in the squeezed limit, which can be evaluated accurately using the nonlinear matter power spectrum and its first-order response functions to density and tidal field perturbations. We demonstrate the accuracy of our model by showing that it reproduces the small-scale $\zeta_{\pm}$ measured in simulated cosmic shear maps. The impact of baryonic feedback enters effectively only through the corresponding impact on the nonlinear matter power spectrum, thereby permitting to account for these astrophysical effects on $\zeta_{\pm}$ similarly to how they are currently accounted for on $\xi_{\pm}$. Using a simple idealized Fisher matrix forecast for a DES-like survey we find that, compared to $\xi_{\pm}$, a combined $\xi_{\pm}\ \&\ \zeta_{\pm}$ analysis can lead to improvements of order $20-40\%$ on the constraints of cosmological parameters such as $\sigma_8$ or the dark energy equation of state parameter $w_0$. We find similar levels of improvement on the constraints of the baryonic feedback parameters, which strengthens the prospects for cosmic shear data to obtain tight constraints not only on cosmology but also on astrophysical feedback models. These are encouraging results that motivate future works on the integrated shear 3-point correlation function towards applications to real survey data., Comment: 15 pages, 9 figures; v2 matches the published version in MNRAS

We study the connection of matter density and its tracers from the PDF perspective. One aspect of this connection is the conditional expectation value $\langle \delta_{\mathrm{tracer}}|\delta_m\rangle$ when averaging both tracer and matter density over some scale. We present a new way to incorporate a Lagrangian bias expansion of this expectation value into standard frameworks for modelling the PDF of density fluctuations and counts-in-cells statistics. Using N-body simulations and mock galaxy catalogs we confirm the accuracy of this expansion and compare it to the more commonly used Eulerian parametrization. For halos hosting typical luminous red galaxies, the Lagrangian model provides a significantly better description of $\langle \delta_{\mathrm{tracer}}|\delta_m\rangle$ at second order in perturbations. A second aspect of the matter-tracer connection is shot-noise, \ie the scatter of tracer density around $\langle \delta_{\mathrm{tracer}}|\delta_m\rangle$. It is well known that this noise can be significantly non-Poissonian and we validate the performance of a more general, two-parameter shot-noise model for different tracers and simulations. Both parts of our analysis are meant to pave the way for forthcoming applications to survey data., Comment: 20 pages, submitted to MNRAS, numerical tools available at https://github.com/OliverFHD/CosMomentum

We present the integrated 3-point shear correlation function $i\zeta_{\pm}$ -- a higher-order statistic of the cosmic shear field -- which can be directly estimated in wide-area weak lensing surveys without measuring the full 3-point shear correlation function, making this a practical and complementary tool to 2-point statistics for weak lensing cosmology. We define it as the 1-point aperture mass statistic $M_{\mathrm{ap}}$ measured at different locations on the shear field correlated with the corresponding local 2-point shear correlation function $\xi_{\pm}$. Building upon existing work on the integrated bispectrum of the weak lensing convergence field, we present a theoretical framework for computing the integrated 3-point function in real space for any projected field within the flat-sky approximation and apply it to cosmic shear. Using analytical formulae for the non-linear matter power spectrum and bispectrum, we model $i\zeta_{\pm}$ and validate it on N-body simulations within the uncertainties expected from the sixth year cosmic shear data of the Dark Energy Survey. We also explore the Fisher information content of $i\zeta_{\pm}$ and perform a joint analysis with $\xi_{\pm}$ for two tomographic source redshift bins with realistic shape-noise to analyse its power in constraining cosmological parameters. We find that the joint analysis of $\xi_{\pm}$ and $i\zeta_{\pm}$ has the potential to considerably improve parameter constraints from $\xi_{\pm}$ alone, and can be particularly useful in improving the figure of merit of the dynamical dark energy equation of state parameters from cosmic shear data., Comment: Accepted for publication in MNRAS; v2 matches the accepted manuscript; 18 pages + appendix