
Weak lensing peak counts are a powerful statistical tool for constraining
cosmological parameters. So far, this method has been applied only to surveys
with relatively small areas, up to several hundred square degrees. As future
surveys will provide weak lensing datasets with size of thousands of square
degrees, the demand on the theoretical prediction of the peak statistics will
become heightened. In particular, large simulations of increased cosmological
volume are required. In this work, we investigate the possibility of using
simulations generated with the fast ComovingLagrangian acceleration (COLA)
method, coupled to the convergence map generator Ufalcon, for predicting the
peak counts. We examine the systematics introduced by the COLA method by
comparing it with a full TreePM code. We find that for a 2000 deg$^2$ survey,
the systematic error is much smaller than the statistical error. This suggests
that the COLA method is able to generate promising theoretical predictions for
weak lensing peaks. We also examine the constraining power of various
configurations of data vectors, exploring the influence of splitting the sample
into tomographic bins and combining different smoothing scales. We find the
combination of smoothing scales to have the most constraining power, improving
the constraints on the $S_8$ amplitude parameter by at least 40% compared to a
single smoothing scale, with tomography brining only limited increase in
measurement precision.

Dark matter in the universe evolves through gravity to form a complex network
of halos, filaments, sheets and voids, that is known as the cosmic web.
Computational models of the underlying physical processes, such as classical
Nbody simulations, are extremely resource intensive, as they track the action
of gravity in an expanding universe using billions of particles as tracers of
the cosmic matter distribution. Therefore, upcoming cosmology experiments will
face a computational bottleneck that may limit the exploitation of their full
scientific potential. To address this challenge, we demonstrate the application
of a machine learning technique called Generative Adversarial Networks (GAN) to
learn models that can efficiently generate new, physically realistic
realizations of the cosmic web. Our training set is a small, representative
sample of 2D image snapshots from Nbody simulations of size 500 and 100 Mpc.
We show that the GANproduced results are qualitatively and quantitatively very
similar to the originals. Generation of a new cosmic web realization with a GAN
takes a fraction of a second, compared to the many hours needed by the Nbody
technique. We anticipate that GANs will therefore play an important role in
providing extremely fast and precise simulations of cosmic web in the era of
large cosmological surveys, such as Euclid and LSST.

Modeling the Point Spread Function (PSF) of widefield surveys is vital for
many astrophysical applications and cosmological probes including weak
gravitational lensing. The PSF smears the image of any recorded object and
therefore needs to be taken into account when inferring properties of galaxies
from astronomical images. In the case of cosmic shear, the PSF is one of the
dominant sources of systematic errors and must be treated carefully to avoid
biases in cosmological parameters. Recently, forward modeling approaches to
calibrate shear measurements within the MonteCarlo Control Loops ($MCCL$)
framework have been developed. These methods typically require simulating a
large amount of widefield images, thus, the simulations need to be very fast
yet have realistic properties in key features such as the PSF pattern. Hence,
such forward modeling approaches require a very flexible PSF model, which is
quick to evaluate and whose parameters can be estimated reliably from survey
data. We present a PSF model that meets these requirements based on a fast
deeplearning method to estimate its free parameters. We demonstrate our
approach on publicly available SDSS data. We extract the most important
features of the SDSS sample via principal component analysis. Next, we
construct our model based on perturbations of a fixed base profile, ensuring
that it captures these features. We then train a Convolutional Neural Network
to estimate the free parameters of the model from noisy images of the PSF. This
allows us to render a model image of each star, which we compare to the SDSS
stars to evaluate the performance of our method. We find that our approach is
able to accurately reproduce the SDSS PSF at the pixel level, which, due to the
speed of both the model evaluation and the parameter estimation, offers good
prospects for incorporating our method into the $MCCL$ framework.

Determining the redshift distribution $n(z)$ of galaxy samples is essential
for several cosmological probes including weak lensing. For imaging surveys,
this is usually done using photometric redshifts estimated on an
objectbyobject basis. We present a new approach for directly measuring the
global $n(z)$ of cosmological galaxy samples, including uncertainties, using
forward modeling. Our method relies on image simulations produced using UFig
(Ultra Fast Image Generator) and on ABC (Approximate Bayesian Computation)
within the $MCCL$ (MonteCarlo Control Loops) framework. The galaxy population
is modeled using parametric forms for the luminosity functions, spectral energy
distributions, sizes and radial profiles of both blue and red galaxies. We
apply exactly the same analysis to the real data and to the simulated images,
which also include instrumental and observational effects. By adjusting the
parameters of the simulations, we derive a set of acceptable models that are
statistically consistent with the data. We then apply the same cuts to the
simulations that were used to construct the target galaxy sample in the real
data. The redshifts of the galaxies in the resulting simulated samples yield a
set of $n(z)$ distributions for the acceptable models. We demonstrate the
method by determining $n(z)$ for a cosmic shear like galaxy sample from the
4band Subaru SuprimeCam data in the COSMOS field. We also complement this
imaging data with a spectroscopic calibration sample from the VVDS survey. We
compare our resulting posterior $n(z)$ distributions to the one derived from
photometric redshifts estimated using 36 photometric bands in COSMOS and find
good agreement. This offers good prospects for applying our approach to current
and future large imaging surveys.

Approximate Bayesian Computation (ABC) is a method to obtain a posterior
distribution without a likelihood function, using simulations and a set of
distance metrics. For that reason, it has recently been gaining popularity as
an analysis tool in cosmology and astrophysics. Its drawback, however, is a
slow convergence rate. We propose a novel method, which we call qABC, to
accelerate ABC with Quantile Regression. In this method, we create a model of
quantiles of distance measure as a function of input parameters. This model is
trained on a small number of simulations and estimates which regions of the
prior space are likely to be accepted into the posterior. Other regions are
then immediately rejected. This procedure is then repeated as more simulations
are available. We apply it to the practical problem of estimation of redshift
distribution of cosmological samples, using forward modelling developed in
previous work. The qABC method converges to nearly same posterior as the basic
ABC. It uses, however, only 20\% of the number of simulations compared to basic
ABC, achieving a fivefold gain in execution time for our problem. For other
problems the acceleration rate may vary; it depends on how close the prior is
to the final posterior. We discuss possible improvements and extensions to this
method.

Weak Gravitational Lensing is a powerful probe of the dark sector of the
Universe. One of the main challenges for this technique is the treatment of
systematics in the measurement of cosmic shear from galaxy shapes. In an
earlier work, Refregier & Amara (2014) have proposed the Monte Carlo Control
Loops (MCCL) to overcome these effects using a forward modeling approach. We
focus here on one of the control loops in this method, the task of which is the
calibration of the shear measurement. For this purpose, we first consider the
requirements on the shear systematics for a given survey and propagate them to
different systematics terms. We use two onepoint statistics to calibrate the
shear measurement and six further onepoint statistics as diagnostics. We also
propagate the systematics levels that we estimate from the onepoint functions
to the twopoint functions for the different systematic error sources. This
allows us to assess the consistency between the systematics levels measured in
different ways. To test the method, we construct synthetic sky surveys with an
area of 1,700 deg$^2$. With some simplifying assumptions, we are able to meet
the requirements on the shear calibration for this survey configuration.
Furthermore, we account for the total residual shear systematics in terms of
the contributing sources. We discuss how this MCCL framework can be applied to
current and future weak lensing surveys.

We demonstrate the potential of Deep Learning methods for measurements of
cosmological parameters from density fields, focusing on the extraction of
nonGaussian information. We consider weak lensing mass maps as our dataset. We
aim for our method to be able to distinguish between five models, which were
chosen to lie along the $\sigma_8$  $\Omega_m$ degeneracy, and have nearly the
same twopoint statistics. We design and implement a Deep Convolutional Neural
Network (DCNN) which learns the relation between five cosmological models and
the mass maps they generate. We develop a new training strategy which ensures
the good performance of the network for high levels of noise. We compare the
performance of this approach to commonly used nonGaussian statistics, namely
the skewness and kurtosis of the convergence maps. We find that our
implementation of DCNN outperforms the skewness and kurtosis statistics,
especially for high noise levels. The network maintains the mean discrimination
efficiency greater than $85\%$ even for noise levels corresponding to ground
based lensing observations, while the other statistics perform worse in this
setting, achieving efficiency less than $70\%$. This demonstrates the ability
of CNNbased methods to efficiently break the $\sigma_8$  $\Omega_m$
degeneracy with weak lensing mass maps alone. We discuss the potential of this
method to be applied to the analysis of real weak lensing data and other
datasets.

We use wavelet and curvelet transforms to extract signals of cosmic strings
from cosmic microwave background (CMB) temperature anisotropy maps, and to
study the limits on the cosmic string tension which various ongoing CMB
temperature anisotropy experiments will be able to achieve. We construct sky
maps with size and angular resolution corresponding to various experiments.
These maps contain the signals of a scaling solution of long string segments
with a given string tension $G \mu$, the contribution of the dominant Gaussian
primordial cosmological fluctuations, and pixel by pixel white noise with an
amplitude corresponding to the instrumental noise of the various experiments.
In the case that we include white noise, we find that the curvelets are more
powerful than wavelets. For maps with Planck specification, we obtain bounds on
the string tension comparable to what was obtained by the Planck collaboration.
Experiments with better angular resolution such as the South Pole Telescope
third generation (SPT3G) survey will be able to yield stronger limits. For
maps with a specification of SPT3G we find that string signals will be visible
down to a string tension of $G \mu = 1.4 \times 10^{7}$.

We present cosmological constraints from the Dark Energy Survey (DES) using a
combined analysis of angular clustering of red galaxies and their
crosscorrelation with weak gravitational lensing of background galaxies. We
use a 139 square degree contiguous patch of DES data from the Science
Verification (SV) period of observations. Using large scale measurements, we
constrain the matter density of the Universe as Omega_m = 0.31 +/ 0.09 and the
clustering amplitude of the matter power spectrum as sigma_8 = 0.74 +/ 0.13
after marginalizing over seven nuisance parameters and three additional
cosmological parameters. This translates into S_8 = sigma_8(Omega_m/0.3)^{0.16}
= 0.74 +/ 0.12 for our fiducial lens redshift bin at 0.35 <z< 0.5, while S_8 =
0.78 +/ 0.09 using two bins over the range 0.2 <z< 0.5. We study the
robustness of the results under changes in the data vectors, modelling and
systematics treatment, including photometric redshift and shear calibration
uncertainties, and find consistency in the derived cosmological parameters. We
show that our results are consistent with previous cosmological analyses from
DES and other data sets and conclude with a joint analysis of DES angular
clustering and galaxygalaxy lensing with Planck CMB data, Baryon Accoustic
Oscillations and Supernova type Ia measurements.

We present first results from the third GRavitational lEnsing Accuracy
Testing (GREAT3) challenge, the third in a sequence of challenges for testing
methods of inferring weak gravitational lensing shear distortions from
simulated galaxy images. GREAT3 was divided into experiments to test three
specific questions, and included simulated space and groundbased data with
constant or cosmologicallyvarying shear fields. The simplest (control)
experiment included parametric galaxies with a realistic distribution of
signaltonoise, size, and ellipticity, and a complex point spread function
(PSF). The other experiments tested the additional impact of realistic galaxy
morphology, multiple exposure imaging, and the uncertainty about a
spatiallyvarying PSF; the last two questions will be explored in Paper II. The
24 participating teams competed to estimate lensing shears to within systematic
error tolerances for upcoming StageIV dark energy surveys, making 1525
submissions overall. GREAT3 saw considerable variety and innovation in the
types of methods applied. Several teams now meet or exceed the targets in many
of the tests conducted (to within the statistical errors). We conclude that the
presence of realistic galaxy morphology in simulations changes shear
calibration biases by $\sim 1$ per cent for a wide range of methods. Other
effects such as truncation biases due to finite galaxy postage stamps, and the
impact of galaxy type as measured by the S\'{e}rsic index, are quantified for
the first time. Our results generalize previous studies regarding sensitivities
to galaxy size and signaltonoise, and to PSF properties such as seeing and
defocus. Almost all methods' results support the simple model in which additive
shear biases depend linearly on PSF ellipticity.

GALSIM is a collaborative, opensource project aimed at providing an image
simulation tool of enduring benefit to the astronomical community. It provides
a software library for generating images of astronomical objects such as stars
and galaxies in a variety of ways, efficiently handling image transformations
and operations such as convolution and rendering at high precision. We describe
the GALSIM software and its capabilities, including necessary theoretical
background. We demonstrate that the performance of GALSIM meets the stringent
requirements of high precision image analysis applications such as weak
gravitational lensing, for current datasets and for the Stage IV dark energy
surveys of the Large Synoptic Survey Telescope, ESA's Euclid mission, and
NASA's WFIRSTAFTA mission. The GALSIM project repository is public and
includes the full code history, all open and closed issues, installation
instructions, documentation, and wiki pages (including a Frequently Asked
Questions section). The GALSIM repository can be found at
https://github.com/GalSimdevelopers/GalSim .

The GRavitational lEnsing Accuracy Testing 3 (GREAT3) challenge is the third
in a series of image analysis challenges, with a goal of testing and
facilitating the development of methods for analyzing astronomical images that
will be used to measure weak gravitational lensing. This measurement requires
extremely precise estimation of very small galaxy shape distortions, in the
presence of far larger intrinsic galaxy shapes and distortions due to the
blurring kernel caused by the atmosphere, telescope optics, and instrumental
effects. The GREAT3 challenge is posed to the astronomy, machine learning, and
statistics communities, and includes tests of three specific effects that are
of immediate relevance to upcoming weak lensing surveys, two of which have
never been tested in a community challenge before. These effects include
realistically complex galaxy models based on highresolution imaging from
space; spatially varying, physicallymotivated blurring kernel; and combination
of multiple different exposures. To facilitate entry by people new to the
field, and for use as a diagnostic tool, the simulation software for the
challenge is publicly available, though the exact parameters used for the
challenge are blinded. Sample scripts to analyze the challenge data using
existing methods will also be provided. See http://great3challenge.info and
http://great3.projects.phys.ucl.ac.uk/leaderboard/ for more information.

Cosmic shear is a powerful probe of cosmological parameters, but its
potential can be fully utilised only if galaxy shapes are measured with great
accuracy. Two major effects have been identified which are likely to account
for most of the bias for maximum likelihood methods in recent shear measurement
challenges. Model bias occurs when the true galaxy shape is not well
represented by the fitted model. Noise bias occurs due to the nonlinear
relationship between image pixels and galaxy shape. In this paper we
investigate the potential interplay between these two effects when an imperfect
model is used in the presence of high noise. We present analytical expressions
for this bias, which depends on the residual difference between the model and
real data. They can lead to biases not accounted for in previous calibration
schemes. By measuring the model bias, noise bias and their interaction, we
provide a complete statistical framework for measuring galaxy shapes with model
fitting methods from GRavitational lEnsing Accuracy Testing (GREAT) like
images. We demonstrate the noise and model interaction bias using a simple toy
model, which indicates that this effect can potentially be significant. Using
real galaxy images from the Cosmological Evolution Survey (COSMOS) we quantify
the strength of the model bias, noise bias and their interaction. We find that
the interaction term is often a similar size to the model bias term, and is
smaller than the requirements of the current and shortly upcoming galaxy
surveys.

We present and describe im3shape, a new publicly available galaxy shape
measurement code for weak gravitational lensing shear. im3shape performs a
maximum likelihood fit of a bulgeplusdisc galaxy model to noisy images,
incorporating an applied point spread function. We detail challenges faced and
choices made in its design and implementation, and then discuss various
limitations that affect this and other maximum likelihood methods. We assess
the bias arising from fitting an incorrect galaxy model using simple noisefree
images and find that it should not be a concern for current cosmic shear
surveys.
We test im3shape on the GREAT08 Challenge image simulations, and meet the
requirements for upcoming cosmic shear surveys in the case that the simulations
are encompassed by the fitted model, using a simple correction for image noise
bias. For the fiducial branch of GREAT08 we obtain a negligible additive shear
bias and subtwo percent level multiplicative bias, which is suitable for
analysis of current surveys. We fall short of the subpercent level requirement
for upcoming surveys, which we attribute to a combination of noise bias and the
mismatch between our galaxy model and the model used in the GREAT08
simulations. We meet the requirements for current surveys across all branches
of GREAT08, except those with small or high noise galaxies, which we would cut
from our analysis. Using the GREAT08 metric we we obtain a score of Q=717 for
the usable branches, relative to the goal of Q=1000 for future experiments. The
code is freely available from https://bitbucket.org/joezuntz/im3shape

Weak gravitational lensing has the potential to constrain cosmological
parameters to high precision. However, as shown by the Shear TEsting Programmes
(STEP) and GRavitational lEnsing Accuracy Testing (GREAT) Challenges, measuring
galaxy shears is a nontrivial task: various methods introduce different
systematic biases which have to be accounted for. We investigate how pixel
noise on the image affects the bias on shear estimates from a
MaximumLikelihood forward modelfitting approach using a sum of coelliptical
S\'{e}rsic profiles, in complement to the theoretical approach of an an
associated paper. We evaluate the bias using a simple but realistic galaxy
model and find that the effects of noise alone can cause biases of order 110%
on measured shears, which is significant for current and future lensing
surveys. We evaluate a simulationbased calibration method to create a bias
model as a function of galaxy properties and observing conditions. This model
is then used to correct the simulated measurements. We demonstrate that this
method can effectively reduce noise bias so that shear measurement reaches the
level of accuracy required for estimating cosmic shear in upcoming lensing
surveys.

Weak lensing experiments are a powerful probe of cosmology through their
measurement of the mass distribution of the universe. A challenge for this
technique is to control systematic errors that occur when measuring the shapes
of distant galaxies. In this paper we investigate noise bias, a systematic
error that arises from second order noise terms in the shape measurement
process. We first derive analytical expressions for the bias of general Maximum
Likelihood Estimators (MLEs) in the presence of additive noise. We then find
analytical expressions for a simplified toy model in which galaxies are modeled
and fitted with a Gaussian with its size as a single free parameter. Even for
this very simple case we find a significant effect. We also extend our analysis
to a more realistic 6parameter elliptical Gaussian model. We find that the
noise bias is generically of the order of the inversesquared signaltonoise
ratio (SNR) of the galaxies and is thus of the order of a percent for galaxies
of SNR of 10, i.e. comparable to the weak lensing shear signal. This is nearly
two orders of magnitude greater than the systematics requirements for future
allsky weak lensing surveys. We discuss possible ways to circumvent this
effect, including a calibration method using simulations discussed in an
associated paper.