Inferring model parameters from experimental data is a grand challenge in many sciences, including cosmology. This often relies critically on high fidelity numerical simulations, which are prohibitively computationally expensive. The application of deep learning techniques to generative modeling is renewing interest in using high dimensional density estimators as computationally inexpensive emulators of fully-fledged simulations. These generative models have the potential to make a dramatic shift in the field of scientific simulations, but for that shift to happen we need to study the performance of such generators in the precision regime needed for science applications. To this end, in this work we apply Generative Adversarial Networks to the problem of generating weak lensing convergence maps. We show that our generator network produces maps that are described by, with high statistical confidence, the same summary statistics as the fully simulated maps.

Lensing peaks have been proposed as a useful statistic, containing cosmological information from non-Gaussianities that is inaccessible from traditional two-point statistics such as the power spectrum or two-point correlation functions. Here we examine constraints on cosmological parameters from weak lensing peak counts, using the publicly available data from the 154 deg$^2$ CFHTLenS survey. We utilize a new suite of ray-tracing N-body simulations on a grid of 91 cosmological models, covering broad ranges of the three parameters $\Omega_m$, $\sigma_8$, and $w$, and replicating the Galaxy sky positions, redshifts, and shape noise in the CFHTLenS observations. We then build an emulator that interpolates the power spectrum and the peak counts to an accuracy of $\leq 5\%$, and compute the likelihood in the three-dimensional parameter space ($\Omega_m$, $\sigma_8$, $w$) from both observables. We find that constraints from peak counts are comparable to those from the power spectrum, and somewhat tighter when different smoothing scales are combined. Neither observable can constrain $w$ without external data. When the power spectrum and peak counts are combined, the area of the error "banana'' in the ($\Omega_m$, $\sigma_8$) plane reduces by a factor of $\approx2$, compared to using the power spectrum alone. For a flat $\Lambda$ cold dark matter model, combining both statistics, we obtain the constraint $\sigma_8(\Omega_m/0.27)^{0.63}=0.85\substack{+0.03 \\ -0.03}$.

Weak gravitational lensing is a powerful cosmological probe, with non--Gaussian features potentially containing the majority of the information. We examine constraints on the parameter triplet $(\Omega_m,w,\sigma_8)$ from non-Gaussian features of the weak lensing convergence field, including a set of moments (up to $4^{\rm th}$ order) and Minkowski functionals, using publicly available data from the 154deg$^2$ CFHTLenS survey. We utilize a suite of ray--tracing N-body simulations spanning 91 points in $(\Omega_m,w,\sigma_8)$ parameter space, replicating the galaxy sky positions, redshifts and shape noise in the CFHTLenS catalogs. We then build an emulator that interpolates the simulated descriptors as a function of $(\Omega_m,w,\sigma_8)$, and use it to compute the likelihood function and parameter constraints. We employ a principal component analysis to reduce dimensionality and to help stabilize the constraints with respect to the number of bins used to construct each statistic. Using the full set of statistics, we find $\Sigma_8\equiv\sigma_8(\Omega_m/0.27)^{0.55}=0.75\pm0.04$ (68% C.L.), in agreement with previous values. We find that constraints on the $(\Omega_m,\sigma_8)$ doublet from the Minkowski functionals suffer a strong bias. However, high-order moments break the $(\Omega_m,\sigma_8)$ degeneracy and provide a tight constraint on these parameters with no apparent bias. The main contribution comes from quartic moments of derivatives.

Residual errors in shear measurements, after corrections for instrument systematics and atmospheric effects, can impact cosmological parameters derived from weak lensing observations. Here we combine convergence maps from our suite of ray-tracing simulations with random realizations of spurious shear. This allows us to quantify the errors and biases of the triplet $(\Omega_m,w,\sigma_8)$ derived from the power spectrum (PS), as well as from three different sets of non-Gaussian statistics of the lensing convergence field: Minkowski functionals (MF), low--order moments (LM), and peak counts (PK). Our main results are: (i) We find an order of magnitude smaller biases from the PS than in previous work. (ii) The PS and LM yield biases much smaller than the morphological statistics (MF, PK). (iii) For strictly Gaussian spurious shear with integrated amplitude as low as its current estimate of $\sigma^2_{sys}\approx 10^{-7}$, biases from the PS and LM would be unimportant even for a survey with the statistical power of LSST. However, we find that for surveys larger than $\approx 100$ deg$^2$, non-Gaussianity in the noise (not included in our analysis) will likely be important and must be quantified to assess the biases. (iv) The morphological statistics (MF,PK) introduce important biases even for Gaussian noise, which must be corrected in large surveys. The biases are in different directions in $(\Omega_m,w,\sigma_8)$ parameter space, allowing self-calibration by combining multiple statistics. Our results warrant follow-up studies with more extensive lensing simulations and more accurate spurious shear estimates.

The statistics of shear peaks have been shown to provide valuable cosmological information beyond the power spectrum, and will be an important constraint of models of cosmology with the large survey areas provided by forthcoming astronomical surveys. Surveys include masked areas due to bright stars, bad pixels etc, which must be accounted for in producing constraints on cosmology from shear maps. We advocate a forward-modeling approach, where the impact of masking (and other survey artifacts) are accounted for in the theoretical prediction of cosmological parameters, rather than removed from survey data. We use masks based on the Deep Lens Survey, and explore the impact of up to 37% of the survey area being masked on LSST and DES-scale surveys. By reconstructing maps of aperture mass, the masking effect is smoothed out, resulting in up to 14% smaller statistical uncertainties compared to simply reducing the survey area by the masked area. We show that, even in the presence of large survey masks, the bias in cosmological parameter estimation produced in the forward-modeling process is ~1%, dominated by bias caused by limited simulation volume. We also explore how this potential bias scales with survey area and find that small survey areas are more significantly impacted by the differences in cosmological structure in the data and simulated volumes, due to cosmic variance.

We compare the efficiency of moments and Minkowski functionals (MFs) in constraining the subset of cosmological parameters (Omega_m,w,sigma_8) using simulated weak lensing convergence maps. We study an analytic perturbative expansion of the MFs in terms of the moments of the convergence field and of its spatial derivatives. We show that this perturbation series breaks down on smoothing scales below 5', while it shows a good degree of convergence on larger scales (15'). Most of the cosmological distinguishing power is lost when the maps are smoothed on these larger scales. We also show that, on scales comparable to 1', where the perturbation series does not converge, cosmological constraints obtained from the MFs are approximately 1.5-2 times better than the ones obtained from the first few moments of the convergence distribution --- provided that the latter include spatial information, either from moments of gradients, or by combining multiple smoothing scales. Including either a set of these moments or the MFs can significantly tighten constraints on cosmological parameters, compared to the conventional method of using the power spectrum alone.

The weak lensing power spectrum is a powerful tool to probe cosmological parameters. Additionally, lensing peak counts contain cosmological information beyond the power spectrum. Both of these statistics can be affected by the preferential selection of source galaxies in patches of the sky with high magnification, as well as by the dilution in the source galaxy surface density in such regions. If not accounted for, these biases introduce systematic errors for cosmological measurements. Here we quantify these systematic errors, using convergence maps from a suite of ray-tracing N-body simulations. At the cut-off magnitude m of on-going and planned major weak lensing surveys, the logarithmic slope of the cumulative number counts s = dlog[n(>m)]/dlog(m) is in the range 0.1 < s < 0.5. At s = 0.2, expected in the I band for LSST, the inferred values of Omega_m, w and sigma_8 are biased by many sigma (where sigma denotes the marginalized error) and therefore the biases will need to be carefully modeled. We also find that the parameters are biased differently in the (Omega_m, w, sigma_8) parameter space when the power spectrum and when the peak counts are used. In particular, w derived from the power spectrum is less affected than w derived from peak counts, while the opposite is true for the best-constrained combination of [sigma_8 Omega_m^gamma] (with gamma=0.62 from the power spectrum and gamma = 0.48 from peak counts). This suggests that the combination of the power spectrum and peak counts can help mitigate the impact of magnification and size biases.

The statistics of peak counts in reconstructed shear maps contain information beyond the power spectrum, and can improve cosmological constraints from measurements of the power spectrum alone if systematic errors can be controlled. We study the effect of galaxy shape measurement errors on predicted cosmological constraints from the statistics of shear peak counts with the Large Synoptic Survey Telescope (LSST). We use the LSST image simulator in combination with cosmological N-body simulations to model realistic shear maps for different cosmological models. We include both galaxy shape noise and, for the first time, measurement errors on galaxy shapes. We find that the measurement errors considered have relatively little impact on the constraining power of shear peak counts for LSST.

Cosmological measurements require the calculation of nontrivial quantities over large datasets. The next generation of survey telescopes (such as DES, PanSTARRS, and LSST) will yield measurements of billions of galaxies. The scale of these datasets, and the nature of the calculations involved, make cosmological calculations ideal models for implementation on graphics processing units (GPUs). We consider two cosmological calculations, the two-point angular correlation function and the aperture mass statistic, and aim to improve the calculation time by constructing code for calculating them on the GPU. Using CUDA, we implement the two algorithms on the GPU and compare the calculation speeds to comparable code run on the CPU. We obtain a code speed-up of between 10 - 180x faster, compared to performing the same calculation on the CPU. The code has been made publicly available.

Peaks in two-dimensional weak lensing (WL) maps contain significant cosmological information, complementary to the WL power spectrum. This has recently been demonstrated using N-body simulations which neglect baryonic effects. Here we employ ray-tracing N-body simulations in which we manually steepen the density profile of each dark matter halo, mimicking the cooling and concentration of baryons into dark matter potential wells. We find, in agreement with previous works, that this causes a significant increase in the amplitude of the WL power spectrum on small scales (spherical harmonic index l>1,000). We then study the impact of the halo concentration increase on the peak counts, and find the following. (i) Low peaks (with convergence 0.02 < kappa_peak < 0.08), remain nearly unaffected. These peaks are created by a constellation of several halos with low masses (10^12-10^13 M_sun) and large angular offsets from the peak center (> 0.5 R_vir); as a result, they are insensitive to the central halo density profiles. These peaks contain most of the cosmological information, and thus provide an unusually sensitive and unbiased probe. (ii) The number of high peaks (with convergence kappa_peak > 0.08) is increased. However, when the baryon effects are neglected in cosmological parameter estimation, then the high peaks lead to a modest bias, comparable to that from the power spectrum on relatively large-scales (l<2000), and much smaller than the bias from the power spectrum on smaller scales (l>2,000). (iii) In the 3D parameter space (sigma_8, Omega_m, w), the biases from the high peaks and the power spectra are in different directions. This suggests the possibility of "self-calibration": the combination of peak counts and power spectrum can simultaneously constrain baryonic physics and cosmological parameters.

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 simulations 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 maps, 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.

We propose counting peaks in weak lensing (WL) maps, as a function of their height, to probe models of dark energy and to constrain cosmological parameters. Because peaks can be identified in two-dimensional WL maps directly, they can provide constraints that are free from potential selection effects and biases involved in identifying and determining the masses of galaxy clusters. We have run cosmological N-body simulations to produce WL convergence maps in three models with different constant values of the dark energy equation of state parameter, w=-0.8, -1, and -1.2, with a fixed normalization of the primordial power spectrum (corresponding to present-day normalizations of sigma8=0.742, 0.798, and 0.839, respectively). By comparing the number of WL peaks in 8 convergence bins in the range of -0.1 < kappa < 0.2, in multiple realizations of a single simulated 3x3 degree field, we show that the first (last) pair of models can be distinguished at the 95% (85%) confidence level. A survey with depth and area (20,000 sq. degrees), comparable to those expected from LSST, should have a factor of approx. 50 better parameter sensitivity. We find that relatively low-amplitude peaks (kappa = 0.03), which typically do not correspond to a single collapsed halo along the line of sight, account for most of this sensitivity. We study a range of smoothing scales and source galaxy redshifts (z_s). With a fixed source galaxy density of 15/arcmin^2, the best results are provided by the smallest scale we can reliably simulate, 1 arcminute, and z_s=2 provides substantially better sensitivity than z_s< 1.5.

During the last ten years a detailed investigation of preheating was performed for chaotic inflation and for hybrid inflation. However, nonperturbative effects during reheating in the new inflation scenario remained practically unexplored. We do a full analysis of preheating in new inflation, using a combination of analytical and numerical methods. We find that the decay of the homogeneous component of the inflaton field and the resulting process of spontaneous symmetry breaking in the simplest models of new inflation usually occurs almost instantly: for the new inflation on the GUT scale it takes only about 5 oscillations of the field distribution. The decay of the homogeneous inflaton field is so efficient because of a combined effect of tachyonic preheating and parametric resonance. At that stage, the homogeneous oscillating inflaton field decays into a collection of waves of the inflaton field, with a typical wavelength of the order of the inverse inflaton mass. This stage usually is followed by a long stage of decay of the inflaton field into other particles, which can be described by the perturbative approach to reheating after inflation. The resulting reheating temperature typically is rather low.

In a broad class of dark energy models, the universe may collapse within a finite time t_c. Here we study a representative model of dark energy with a linear potential, V(\phi)=V_0(1+\alpha\phi). This model is the simplest doomsday model, in which the universe collapses rather quickly after it stops expanding. Observational data from type Ia supernovae (SNe Ia), cosmic microwave background anisotropy (CMB), and large scale structure (LSS) are complementary in constraining dark energy models. Using the new SN Ia data (Riess sample), the CMB data from WMAP, and the LSS data from 2dF, we find that the collapse time of the universe is t_c > 42 (24) gigayears from today at 68% (95%) confidence.

It may be difficult to single out the best model of dark energy on the basis of the existing and planned cosmological observations, because many different models can lead to similar observational consequences. However, each particular model can be studied and either found consistent with observations or ruled out. In this paper, we concentrate on the possibility to test and rule out the simplest and by far the most popular of the models of dark energy, the theory described by general relativity with positive vacuum energy (the cosmological constant). We evaluate the conditions under which this model could be ruled out by the future observations made by the Supernova/Acceleration Probe SNAP (both for supernovae and weak lensing) and by the Planck Surveyor cosmic microwave background satellite.

Recently it was found, in a broad class of models, that the dark energy density may change its sign during the evolution of the universe. This may lead to a global collapse of the universe within the time t_c ~ 10^{10}-10^{11} years. Our goal is to find what bounds on the future lifetime of the universe can be placed by the next generation of cosmological observations. As an example, we investigate the simplest model of dark energy with a linear potential V(\phi) =V_0(1+\alpha\phi). This model can describe the present stage of acceleration of the universe if \alpha is small enough. However, eventually the field \phi rolls down, V(\phi) becomes negative, and the universe collapses. The existing observational data indicate that the universe described by this model will collapse not earlier than t_c > 10 billion years from the present moment. We show that the data from SNAP and Planck satellites may extend the bound on the "doomsday" time to t_c > 40 billion years at the 95% confidence level.