In this paper, we show that Minkowski Functionals (MFs) of weak gravitational lensing (WL) convergence maps contain significant non-Gaussian, cosmology-dependent information. To do this, we use a large suite of cosmological ray-tracing N-body simulat
ions to create mock WL convergence maps, and study the cosmological information content of MFs derived from these maps. Our suite consists of 80 independent 512^3 N-body runs, covering seven different cosmologies, varying three cosmological parameters Omega_m, w, and sigma_8 one at a time, around a fiducial LambdaCDM model. In each cosmology, we use ray-tracing to create a thousand pseudo-independent 12 deg^2 convergence maps, and use these in a Monte Carlo procedure to estimate the joint confidence contours on the above three parameters. We include redshift tomography at three different source redshifts z_s=1, 1.5, 2, explore five different smoothing scales theta_G=1, 2, 3, 5, 10 arcmin, and explicitly compare and combine the MFs with the WL power spectrum. We find that the MFs capture a substantial amount of information from non-Gaussian features of convergence maps, i.e. beyond the power spectrum. The MFs are particularly well suited to break degeneracies and to constrain the dark energy equation of state parameter w (by a factor of ~ three better than from the power spectrum alone). The non-Gaussian information derives partly from the one-point function of the convergence (through V_0, the area MF), and partly through non-linear spatial information (through combining different smoothing scales for V_0, and through V_1 and V_2, the boundary length and genus MFs, respectively). In contrast to the power spectrum, the best constraints from the MFs are obtained only when multiple smoothing scales are combined.
Recent studies have shown that the number counts of convergence peaks N(kappa) in weak lensing (WL) maps, expected from large forthcoming surveys, can be a useful probe of cosmology. We follow up on this finding, and use a suite of WL convergence map
s, obtained from ray-tracing N-body simulations, to study (i) the physical origin of WL peaks with different heights, and (ii) whether the peaks contain information beyond the convergence power spectrum P_ell. In agreement with earlier work, we find that high peaks (with amplitudes >~ 3.5 sigma, where sigma is the r.m.s. of the convergence kappa) are typically dominated by a single massive halo. In contrast, medium-height peaks (~0.5-1.5 sigma) cannot be attributed to a single collapsed dark matter halo, and are instead created by the projection of multiple (typically, 4-8) halos along the line of sight, and by random galaxy shape noise. Nevertheless, these peaks dominate the sensitivity to the cosmological parameters w, sigma_8, and Omega_m. We find that the peak height distribution and its dependence on cosmology differ significantly from predictions in a Gaussian random field. We directly compute the marginalized errors on w, sigma_8, and Omega_m from the N(kappa) + P_ell combination, including redshift tomography with source galaxies at z_s=1 and z_s=2. We find that the N(kappa) + P_ell combination has approximately twice the cosmological sensitivity compared to P_ell alone. These results demonstrate that N(kappa) contains non-Gaussian information complementary to the power spectrum.