Deep Neural Networks (DNNs) are powerful algorithms that have been proven capable of extracting non-Gaussian information from weak lensing (WL) data sets. Understanding which features in the data determine the output of these nested, non-linear algorithms is an important but challenging task. We analyze a DNN that has been found in previous work to accurately recover cosmological parameters in simulated maps of the WL convergence ($\kappa$). We derive constraints on the cosmological parameter pair $(\Omega_m,\sigma_8)$ from a combination of three commonly used WL statistics (power spectrum, lensing peaks, and Minkowski functionals), using ray-traced simulated $\kappa$ maps. We show that the network can improve the inferred parameter constraints relative to this combination by $20\%$ even in the presence of realistic levels of shape noise. We apply a series of well established saliency methods to interpret the DNN and find that the most relevant pixels are those with extreme $\kappa$ values. For noiseless maps, regions with negative $\kappa$ account for $86-69\%$ of the attribution of the DNN output, defined as the square of the saliency in input space. In the presence of shape nose, the attribution concentrates in high convergence regions, with $36-68\%$ of the attribution in regions with $\kappa > 3 \sigma_{\kappa}$.

We performed a series of numerical experiments to quantify the sensitivity of the predictions for weak lensing statistics obtained in raytracing DM-only simulations, to two hyper-parameters that influence the accuracy as well as the computational cost of the predictions: the thickness of the lens planes used to build past light-cones and the mass resolution of the underlying DM simulation. The statistics considered are the power spectrum and a series of non-Gaussian observables, including the one-point probability density function, lensing peaks, and Minkowski functionals. Counter-intuitively, we find that using thin lens planes ($< 60~h^{-1}$Mpc on a $240~h^{-1}$Mpc simulation box) suppresses the power spectrum over a broad range of scales beyond what would be acceptable for an LSST-type survey. A mass resolution of $7.2\times 10^{11}~h^{-1}\,M_{\odot}$ per DM particle (or 256$^3$ particles in a ($240~h^{-1}$Mpc)$^3$ box) is sufficient to extract information using the power spectrum and non-Gaussian statistics from weak lensing data at angular scales down to 1 arcmin with LSST-like levels of shape noise.

The rotational kinematic Sunyaev-Zeldovich (rkSZ) signal, imprinted on the cosmic microwave background (CMB) by the gaseous halos (spinning "atmospheres") of foreground galaxies, would be a novel probe of galaxy formation. Although the signal is too weak to detect in individual galaxies, we analyze the feasibility of its statistical detection via stacking CMB data on many galaxies for which the spin orientation can be estimated spectroscopically. We use an "optimistic" model, in which fully ionized atmospheres contain the cosmic baryon fraction and spin at the halo's circular velocity $v_{\rm circ}$, and a more realistic model, based on hydrodynamical simulations, with multi-phase atmospheres spinning at a fraction of $v_{\rm circ}$. We incorporate realistic noise estimates into our analysis. Using low-redshift galaxy properties from the MaNGA spectroscopic survey (with median halo mass of $6.6\times10^{11}\,M_{\odot}$), and CMB data quality from Planck, we find that a $3\sigma$ detection would require a few$\times 10^4$ galaxies, even in the optimistic model. This is too high for current surveys, but upcoming higher-angular resolution CMB experiments will significantly reduce the requirements: stacking CMB data on galaxy spins in a $\sim10$ deg$^2$ can rule out the optimistic models, and $\approx$350 deg$^2$ will suffice for a $3\sigma$ detection with ACT. As a proof-of-concept, we stacked Planck data on the position of $\approx2,000$ MaNGA galaxies, aligned with the galaxies' projected spin, and scaled to their halos' angular size. We rule out average temperature dipoles larger than $\approx1.9\,\mu$K around field spiral galaxies.

Weak gravitational lensing is one of the most promising cosmological probes of the late universe. Several large ongoing (DES, KiDS, HSC) and planned (LSST, EUCLID, WFIRST) astronomical surveys attempt to collect even deeper and larger scale data on weak lensing. Due to gravitational collapse, the distribution of dark matter is non-Gaussian on small scales. However, observations are typically evaluated through the two-point correlation function of galaxy shear, which does not capture non-Gaussian features of the lensing maps. Previous studies attempted to extract non-Gaussian information from weak lensing observations through several higher-order statistics such as the three-point correlation function, peak counts or Minkowski-functionals. Deep convolutional neural networks (CNN) emerged in the field of computer vision with tremendous success, and they offer a new and very promising framework to extract information from 2 or 3-dimensional astronomical data sets, confirmed by recent studies on weak lensing. We show that a CNN is able to yield significantly stricter constraints of ($\sigma_8, \Omega_m$) cosmological parameters than the power spectrum using convergence maps generated by full N-body simulations and ray-tracing, at angular scales and shape noise levels relevant for future observations. In a scenario mimicking LSST or Euclid, the CNN yields 2.4-2.8 times smaller credible contours than the power spectrum, and 3.5-4.2 times smaller at noise levels corresponding to a deep space survey such as WFIRST. We also show that at shape noise levels achievable in future space surveys the CNN yields 1.4-2.1 times smaller contours than peak counts, a higher-order statistic capable of extracting non-Gaussian information from weak lensing maps.

The presence of massive neutrinos affects structure formation, leaving imprints on large-scale structure observables such as the weak lensing field. The common lensing analyses with two-point statistics are insensitive to the large amount of non-Gaussian information in the density field. We investigate non-Gaussian tools, in particular the Minkowski Functionals (MFs)---morphological descriptors including area, perimeter, and genus---in an attempt to recover the higher-order information. We use convergence maps from the Cosmological Massive Neutrino Simulations (\texttt{MassiveNus}) and assume galaxy noise, density, and redshift distribution for an LSST-like survey. We show that MFs are sensitive to the neutrino mass sum, and the sensitivity is redshift dependent and is non-Gaussian. We find that redshift tomography significantly improves the constraints on neutrino mass for MFs, compared to the improvements for the power spectrum. We attribute this to the stronger redshift dependence of neutrino effects on small scales. We then build an emulator to model the power spectrum and MFs, and study the constraints on $[M_{\nu}$, $\Omega_{m}$, $A_{s}]$ from the power spectrum, MFs, and their combination. We show that MFs significantly outperform the power spectrum in constraining neutrino mass, by more than a factor of four. However, a thorough study of the impact from systematics such as baryon physics and galaxy shape and redshift biases will be important to realize the full potential of MFs.

Massive cosmic neutrinos change the structure formation history by suppressing perturbations on small scales. Weak lensing data from galaxy surveys probe the structure evolution and thereby can be used to constrain the total mass of the three active neutrinos. However, much of the information is at small scales where the dynamics are nonlinear. Traditional approaches with second order statistics thus fail to fully extract the information in the lensing field. In this paper, we study constraints on the neutrino mass sum using lensing peak counts, a statistic which captures non-Gaussian information beyond the second order. We use the ray-traced weak lensing mocks from the Cosmological Massive Neutrino Simulations (MassiveNuS), and apply LSST-like noise. We discuss the effects of redshift tomography, multipole cutoff $\ell_{\rm max}$ for the power spectrum, smoothing scale for the peak counts, and constraints from peaks of different heights. We find that combining peak counts with the traditional lensing power spectrum can improve the constraint on neutrino mass sum, $\Omega_m$, and $A_s$ by 39%, 32%, and 60%, respectively, over that from the power spectrum alone.

Weak lensing maps contain information beyond two-point statistics on small scales. Much recent work has tried to extract this information through a range of different observables or via nonlinear transformations of the lensing field. Here we train and apply a 2D convolutional neural network to simulated noiseless lensing maps covering 96 different cosmological models over a range of {$\Omega_m,\sigma_8$}. Using the area of the confidence contour in the {$\Omega_m,\sigma_8$} plane as a figure-of-merit, derived from simulated convergence maps smoothed on a scale of 1.0 arcmin, we show that the neural network yields $\approx 5 \times$ tighter constraints than the power spectrum, and $\approx 4 \times$ tighter than the lensing peaks. Such gains illustrate the extent to which weak lensing data encode cosmological information not accessible to the power spectrum or even other, non-Gaussian statistics such as lensing peaks.

The non-zero mass of neutrinos suppresses the growth of cosmic structure on small scales. Since the level of suppression depends on the sum of the masses of the three active neutrino species, the evolution of large-scale structure is a promising tool to constrain the total mass of neutrinos and possibly shed light on the mass hierarchy. In this work, we investigate these effects via a large suite of N-body simulations that include massive neutrinos using an analytic linear-response approximation: the Cosmological Massive Neutrino Simulations (MassiveNuS). The simulations include the effects of radiation on the background expansion, as well as the clustering of neutrinos in response to the nonlinear dark matter evolution. We allow three cosmological parameters to vary: the neutrino mass sum M_nu in the range of 0-0.6 eV, the total matter density Omega_m, and the primordial power spectrum amplitude A_s. The rms density fluctuation in spheres of 8 comoving Mpc/h (sigma_8) is a derived parameter as a result. Our data products include N-body snapshots, halo catalogues, merger trees, ray- traced galaxy lensing convergence maps for four source redshift planes between z_s=1-2.5, and ray-traced cosmic microwave background lensing convergence maps. We describe the simulation procedures and code validation in this paper. The data are publicly available at http://columbialensing.org.

We explore the sensitivity of weak lensing observables to the expansion history of the universe and to the growth of cosmic structures, as well as the relative contribution of both effects to constraining cosmological parameters. We utilize ray-tracing dark-matter-only N-body simulations and validate our technique by comparing our results for the convergence power spectrum with analytic results from past studies. We then extend our analysis to non-Gaussian observables which cannot be easily treated analytically. We study the convergence (equilateral) bispectrum and two topological observables, lensing peaks and Minkowski functionals, focusing on their sensitivity to the matter density $\Omega_m$ and the dark energy equation of state $w$. We find that a cancelation between the geometry and growth effects is a common feature for all observables, and exists at the map level. It weakens the overall sensitivity by up to a factor of 3 and 1.5 for $w$ and $\Omega_m$, respectively, with the bispectrum worst affected. However, combining geometry and growth information alleviates the degeneracy between $\Omega_m$ and $w$ from either effect alone. As a result, the magnitude of marginalized errors remain similar to those obtained from growth-only effects, but with the correlation between the two parameters switching sign. These results shed light on the origin of cosmology-sensitivity of non-Gaussian statistics, and should be useful in optimizing combinations of observables.

We have investigated a recently proposed halo-based model, Camelus, for predicting weak-lensing peak counts, and compared its results over a collection of 162 cosmologies with those from N-body simulations. While counts from both models agree for peaks with $\mathcal{S/N}>1$ (where $\mathcal{S/N}$ is the ratio of the peak height to the r.m.s. shape noise), we find $\approx 50\%$ fewer counts for peaks near $\mathcal{S/N}=0$ and significantly higher counts in the negative $\mathcal{S/N}$ tail. Adding shape noise reduces the differences to within $20\%$ for all cosmologies. We also found larger covariances that are more sensitive to cosmological parameters. As a result, credibility regions in the $\{\Omega_m, \sigma_8\}$ are $\approx 30\%$ larger. Even though the credible contours are commensurate, each model draws its predictive power from different types of peaks. Low peaks, especially those with $2<\mathcal{S/N}<3$, convey important cosmological information in N-body data, as shown in \cite{DietrichHartlap, Kratochvil2010}, but \textsc{Camelus} constrains cosmology almost exclusively from high significance peaks $(\mathcal{S/N}>3)$. Our results confirm the importance of using a cosmology-dependent covariance with at least a 14\% improvement in parameter constraints. We identified the covariance estimation as the main driver behind differences in inference, and suggest possible ways to make Camelus even more useful as a highly accurate peak count emulator.