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!

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

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

Using the M31 PAndromeda Cepheid sample and the HST PHAT data we obtain the largest Cepheid sample in M31 with HST data in four bands. For our analysis we consider three samples: A very homogeneous sample of Cepheids based on the PAndromeda data, the mean magnitude corrected PAndromeda sample and a sample complementing the PAndromeda sample with Cepheids from literature. The latter results in the largest catalog with 522 fundamental mode (FM) Cepheids and 102 first overtone (FO) Cepheids with F160W and F110W data and 559 FM Cepheids and 111 FO Cepheids with F814W and F475W data. The obtained dispersion of the Period-Luminosity relations (PLRs) is very small (e.g. 0.138 mag in the F160W sample I PLR). We find no broken slope in the PLRs when analyzing our entire sample, but we do identify a subsample of Cepheids that causes the broken slope. However, this effect only shows when the number of this Cepheid type makes up a significant fraction of the total sample. We also analyze the sample selection effect on the Hubble constant., Comment: 32 pages, 19 figures, 9 tables, accepted for publication in ApJ, electronic data will be available on CDS

We present results from the Wendelstein Weak Lensing (WWL) pathfinder project, in which we have observed three intermediate redshift Planck clusters of galaxies with the new 30'$\times 30$' wide field imager at the 2m Fraunhofer Telescope at Wendelstein Observatory. We investigate the presence of biases in our shear catalogues and estimate their impact on our weak lensing mass estimates. The overall calibration uncertainty depends on the cluster redshift and is below 8.1-15 per cent for $z \approx 0.27-0.77$. It will decrease with improvements on the background sample selection and the multiplicative shear bias calibration. We present the first weak lensing mass estimates for PSZ1 G109.88+27.94 and PSZ1 G139.61+24.20, two SZ-selected cluster candidates. Based on Wendelstein colors and SDSS photometry, we find that the redshift of PSZ1 G109.88+27.94 has to be corrected to $z \approx 0.77$. We investigate the influence of line-of-sight structures on the weak lensing mass estimates and find upper limits for two groups in each of the fields of PSZ1 G109.88+27.94 and PSZ1 G186.98+38.66. We compare our results to SZ and dynamical mass estimates from the literature, and in the case of PSZ1 G186.98+38.66 to previous weak lensing mass estimates. We conclude that our pathfinder project demonstrates that weak lensing cluster masses can be accurately measured with the 2m Fraunhofer Telescope., Comment: 24 pages, 26 figures. Submitted to MNRAS

We present the largest Cepheid sample in M31 based on the complete Pan-STARRS1 survey of Andromeda (PAndromeda) in the $r_{\mathrm{P1}}$ , $i_{\mathrm{P1}}$ and $g_{\mathrm{P1}}$ bands. We find 2686 Cepheids with 1662 fundamental mode Cepheids, 307 first-overtone Cepheids, 278 type II Cepheids and 439 Cepheids with undetermined Cepheid type. Using the method developed by Kodric et al. (2013) we identify Cepheids by using a three dimensional parameter space of Fourier parameters of the Cepheid light curves combined with a color cut and other selection criteria. This is an unbiased approach to identify Cepheids and results in a homogeneous Cepheid sample. The Period-Luminosity relations obtained for our sample have smaller dispersions than in our previous work. We find a broken slope that we previously observed with HST data in Kodric et al. (2015), albeit with a lower significance., Comment: 79 pages, 39 figures, 8 tables, accepted for publication in AJ, K18b is submittted to ApJ, electronic data will be available on CDS

We reconstruct the two-dimensional (2D) matter distributions in 20 high-mass galaxy clusters selected from the CLASH survey by using the new approach of performing a joint weak lensing analysis of 2D shear and azimuthally averaged magnification measurements. This combination allows for a complete analysis of the field, effectively breaking the mass-sheet degeneracy. In a Bayesian framework, we simultaneously constrain the mass profile and morphology of each individual cluster assuming an elliptical Navarro-Frenk-White halo characterized by the mass, concentration, projected axis ratio, and position angle of the projected major axis.. We find that spherical mass estimates of the clusters from azimuthally averaged weak-lensing measurements in previous work are in excellent agreement with our results from a full 2D analysis. Combining all 20 clusters in our sample, we detect the elliptical shape of weak-lensing halos at the $5\sigma$ significance level within a scale of 2Mpc$/h$. The median projected axis ratio is $0.67\pm 0.07$ at a virial mass of $M_\mathrm{vir}=(15.2\pm 2.8) \times 10^{14} M_\odot$, which is in agreement with theoretical predictions of the standard collisionless cold dark matter model. We also study misalignment statistics of the brightest cluster galaxy, X-ray, thermal Sunyaev-Zel'dovich effect, and strong-lensing morphologies with respect to the weak-lensing signal. Among the three baryonic tracers studied here, we find that the X-ray morphology is best aligned with the weak-lensing mass distribution, with a median misalignment angle of $21\pm 7$ degrees. We also conduct a stacked quadrupole shear analysis assuming that the X-ray major axis is aligned with that of the projected mass distribution. This yields a consistent axis ratio of $0.67\pm 0.10$, suggesting again a tight alignment between the intracluster gas and dark matter., Comment: Minor changes to match the version published in ApJ. One of three new companion papers of the CLUMP-3D project (I-Non Chiu et al., arXiv:1804.00676; Mauro Sereno et al., arXiv:1804.00667)