Jessica Cowell

Abstract:

Abstract We present the first application of marked power spectra to weak lensing data, using maps from the Subaru Hyper Suprime-Cam Year 1 (HSC-Y1) survey. Marked convergence fields, constructed by weighting the convergence field with non-linear functions of its smoothed version, are designed to encode higher-order information while remaining computationally tractable. Using simulations tailored to the HSC-Y1 data, we test three mark functions that up- or down-weight different density environments. Our results show that combining multiple types of marked auto- and cross-spectra improves constraints on the clustering amplitude parameter $S_8\equiv \sigma _8\sqrt{\Omega _{\rm m}/0.3}$ by ≈43 percnt compared to standard two-point power spectra. When applied to the HSC-Y1 data, this translates into a constraint on S8 = 0.807 ± 0.024. We assess the sensitivity of the marked power spectra to systematics, including baryonic effects, intrinsic alignment, photometric redshifts, and multiplicative shear bias. We note that some of the additional information introduced by the marked field originates from scales smaller than the scale cut, and is partly Gaussian in nature. This does not invalidate our systematic tests. These results demonstrate the promise of marked statistics as a practical and powerful tool for extracting non-Gaussian information from weak lensing surveys.

Abstract:

<jats:p>The study of third-order statistics in large-scale structure analyses has been hampered by the increased complexity of bispectrum estimators (compared to power spectra), the large dimensionality of the data vector, and the difficulty in estimating its covariance matrix. In this paper we present the filtered-squared bispectrum (FSB), an estimator of the projected bispectrum effectively consisting of the cross-correlation between the square of a field filtered on a range of scales and the original field. Within this formalism, we are able to recycle much of the infrastructure built around power spectrum measurement to construct an estimator that is both fast and robust against mode-coupling effects caused by incomplete sky observations. Furthermore, we demonstrate that the existing techniques for the estimation of analytical power spectrum covariances can be used within this formalism to calculate the bispectrum covariance at very high accuracy, naturally accounting for the most relevant Gaussian and non-Gaussian contributions in a model-independent manner.</jats:p>