Normalized to: Kacprzak, T.

**47** article(s) in total. **408** co-authors, from **1** to **31** common article(s). Median position in authors list is **6,0**.

[1] oai:arXiv.org:1911.08494 [pdf] - 2071705

Baryonic effects for weak lensing. Part II. Combination with X-ray data
and extended cosmologies

Comments: Accepted version (JCAP)

Submitted: **2019-11-19**, last modified: **2020-03-30**

An accurate modelling of baryonic feedback effects is required to exploit the
full potential of future weak-lensing surveys such as Euclid or LSST. In this
second paper in a series of two, we combine Euclid-like mock data of the cosmic
shear power spectrum with an eROSITA X-ray mock of the cluster gas fraction to
run a combined likelihood analysis including both cosmological and baryonic
parameters. Following the first paper of this series, the baryonic effects
(based on the baryonic correction model of Schneider et al. 2019) are included
in both the tomographic power spectrum and the covariance matrix. However, this
time we assume the more realistic case of a $\Lambda$CDM cosmology with massive
neutrinos, and we consider several extensions of the currently favoured
cosmological model. For the standard $\Lambda$CDM case, we show that including
X-ray data reduces the uncertainties on the sum of the neutrino mass by
$\sim30$ percent, while there is only a mild improvement on other parameters
such as $\Omega_m$ and $\sigma_8$. As extensions of $\Lambda$CDM, we consider
the cases of a dynamical dark energy model (wCDM), a $f(R)$ gravity model
(fRCDM), and a mixed dark matter model ($\Lambda$MDM) with both a cold and a
warm/hot dark matter component. We find that combining weak lensing with X-ray
data only leads to a mild improvement of the constraints on the additional
parameters of wCDM, while the improvement is more substantial for both fRCDM
and $\Lambda$MDM. Ignoring baryonic effects in the analysis pipeline leads to
significant false-detections of either phantom dark energy or a light
subdominant dark matter component. Overall we conclude that for all cosmologies
considered, a general parametrisation of baryonic effects is both necessary and
sufficient to obtain tight constraints on cosmological parameters.

[2] oai:arXiv.org:2001.07727 [pdf] - 2035852

Measurement of the B-band Galaxy Luminosity Function with Approximate
Bayesian Computation

Comments: Submitted to JCAP, 42 pages, 20 figures, 4 appendices

Submitted: **2020-01-21**

The galaxy Luminosity Function (LF) is a key observable for galaxy formation,
evolution studies and for cosmology. In this work, we propose a novel technique
to forward model wide-field galaxy surveys using the fast image simulator UFig
and measure the LF of galaxies in the B-band. We use Approximate Bayesian
Computation (ABC) to constrain the galaxy population model parameters of the
simulations and match data from the Canada-France-Hawaii Telescope Legacy
Survey (CFHTLS). We define a number of distance metrics between the simulated
and the survey data. By exploring the parameter space of the galaxy population
model through ABC to find the set of parameters that minimize these distance
metrics, we obtain constraints on the LFs of blue and red galaxies as a
function of redshift. We find that $\mathrm{M^*}$ fades by $\Delta
\mathrm{M}^*_{\mathrm{0.1-1.0,b}} = 0.37 \pm 0.34$ and $\Delta
\mathrm{M}^*_{\mathrm{0.1-1.0,r}} = 0.58 \pm 0.42$ magnitudes between redshift
$\mathrm{z = 1}$ and $\mathrm{z = 0.1}$ for blue and red galaxies,
respectively. We also find that $\phi^*$ for blue galaxies stays roughly
constant between redshift $\mathrm{z = 0.1}$ and $\mathrm{z=1}$, while for red
galaxies it decreases of $\sim 90\%$. We compare our results to other
measurements, finding good agreement at all redshifts, for both blue and red
galaxies. To further test our results, we compare the redshift distributions
for survey and simulated data. We use the spectroscopic redshift distribution
from the VIMOS Public Extragalactic Redshift Survey (VIPERS) and we apply the
same selection in colours and magnitudes on our simulations. We find a good
agreement between the survey and the simulated redshift distributions. We
provide best-fit values and uncertainties for the parameters of the LF. This
work offers excellent prospects for measuring other galaxy population
properties as a function of redshift using ABC.

[3] oai:arXiv.org:1908.05519 [pdf] - 2017214

Cosmological N-body simulations: a challenge for scalable generative
models

Comments:

Submitted: **2019-08-15**, last modified: **2019-12-18**

Deep generative models, such as Generative Adversarial Networks (GANs) or
Variational Autoencoders (VAs) have been demonstrated to produce images of high
visual quality. However, the existing hardware severely limits the size of the
images that can be generated. The rapid growth of high dimensional data in many
fields of science therefore poses a significant challenge for generative
models. In cosmology, the large-scale, three-dimensional matter distribution,
modeled with N-body simulations, plays a crucial role in understanding the
evolution of the universe. As these simulations are computationally very
expensive, GANs have recently generated interest as a possible method to
emulate these datasets, but they have been, so far, mostly limited to two
dimensional data. In this work, we introduce a new benchmark for the generation
of three dimensional N-body simulations, in order to stimulate new ideas in the
machine learning community and move closer to the practical use of generative
models in cosmology. As a first benchmark result, we propose a scalable GAN
approach for training a generator of N-body three-dimensional cubes. Our
technique relies on two key building blocks, (i) splitting the generation of
the high-dimensional data into smaller parts, and (ii) using a multi-scale
approach that efficiently captures global image features that might otherwise
be lost in the splitting process. We evaluate the performance of our model for
the generation of N-body samples using various statistical measures commonly
used in cosmology. Our results show that the proposed model produces samples of
high visual quality, although the statistical analysis reveals that capturing
rare features in the data poses significant problems for the generative models.
We make the data, quality evaluation routines, and the proposed GAN
architecture publicly available at https://github.com/nperraud/3DcosmoGAN

[4] oai:arXiv.org:1905.11636 [pdf] - 1975276

Effects of baryons on weak lensing peak statistics

Comments: 22 pages, 11 figures, JCAP accepted

Submitted: **2019-05-28**, last modified: **2019-10-01**

Upcoming weak-lensing surveys have the potential to become leading
cosmological probes provided all systematic effects are under control.
Recently, the ejection of gas due to feedback energy from active galactic
nuclei (AGN) has been identified as major source of uncertainty, challenging
the success of future weak-lensing probes in terms of cosmology. In this paper
we investigate the effects of baryons on the number of weak-lensing peaks in
the convergence field. Our analysis is based on full-sky convergence maps
constructed via light-cones from $N$-body simulations, and we rely on the
baryonic correction model of Schneider et al. (2019) to model the baryonic
effects on the density field. As a result we find that the baryonic effects
strongly depend on the Gaussian smoothing applied to the convergence map. For a
DES-like survey setup, a smoothing of $\theta_k\gtrsim8$ arcmin is sufficient
to keep the baryon signal below the expected statistical error. Smaller
smoothing scales lead to a significant suppression of high peaks (with
signal-to-noise above 2), while lower peaks are not affected. The situation is
more severe for a Euclid-like setup, where a smoothing of $\theta_k\gtrsim16$
arcmin is required to keep the baryonic suppression signal below the
statistical error. Smaller smoothing scales require a full modelling of
baryonic effects since both low and high peaks are strongly affected by
baryonic feedback.

[5] oai:arXiv.org:1906.03156 [pdf] - 1961599

Cosmological constraints with deep learning from KiDS-450 weak lensing
maps

Comments: 22 pages, 15 figures

Submitted: **2019-06-07**, last modified: **2019-09-16**

Convolutional Neural Networks (CNN) have recently been demonstrated on
synthetic data to improve upon the precision of cosmological inference. In
particular they have the potential to yield more precise cosmological
constraints from weak lensing mass maps than the two-point functions. We
present the cosmological results with a CNN from the KiDS-450 tomographic weak
lensing dataset, constraining the total matter density $\Omega_m$, the
fluctuation amplitude $\sigma_8$, and the intrinsic alignment amplitude
$A_{\rm{IA}}$. We use a grid of N-body simulations to generate a training set
of tomographic weak lensing maps. We test the robustness of the expected
constraints to various effects, such as baryonic feedback, simulation accuracy,
different value of $H_0$, or the lightcone projection technique. We train a set
of ResNet-based CNNs with varying depths to analyze sets of tomographic KiDS
mass maps divided into 20 flat regions, with applied Gaussian smoothing of
$\sigma=2.34$ arcmin. The uncertainties on shear calibration and $n(z)$ error
are marginalized in the likelihood pipeline. Following a blinding scheme, we
derive constraints of $S_8 = \sigma_8 (\Omega_m/0.3)^{0.5} =
0.777^{+0.038}_{-0.036}$ with our CNN analysis, with
$A_{\rm{IA}}=1.398^{+0.779}_{-0.724}$. We compare this result to the power
spectrum analysis on the same maps and likelihood pipeline and find an
improvement of about $30\%$ for the CNN. We discuss how our results offer
excellent prospects for the use of deep learning in future cosmological data
analysis.

[6] oai:arXiv.org:1906.01018 [pdf] - 1894429

Monte Carlo Control Loops for cosmic shear cosmology with DES Year 1

Comments: 31 pages, 15 figures, 2 tables

Submitted: **2019-06-03**

Weak lensing by large-scale structure is a powerful probe of cosmology and of
the dark universe. This cosmic shear technique relies on the accurate
measurement of the shapes and redshifts of background galaxies and requires
precise control of systematic errors. The Monte Carlo Control Loops (MCCL) is a
forward modelling method designed to tackle this problem. It relies on the
Ultra Fast Image Generator (UFig) to produce simulated images tuned to match
the target data statistically, followed by calibrations and tolerance loops. We
present the first end-to-end application of this method, on the Dark Energy
Survey (DES) Year 1 wide field imaging data. We simultaneously measure the
shear power spectrum $C_{\ell}$ and the redshift distribution $n(z)$ of the
background galaxy sample. The method includes maps of the systematic sources,
Point Spread Function (PSF), an Approximate Bayesian Computation (ABC)
inference of the simulation model parameters, a shear calibration scheme, and
the fast estimation of the covariance matrix. We find a close statistical
agreement between the simulations and the DES Y1 data using an array of
diagnostics. In a non-tomographic setting, we derive a set of $C_\ell$ and
$n(z)$ curves that encode the cosmic shear measurement, as well as the
systematic uncertainty. Following a blinding scheme, we measure the combination
of $\Omega_m$, $\sigma_8$, and intrinsic alignment amplitude $A_{\rm{IA}}$,
defined as $S_8D_{\rm{IA}} = \sigma_8(\Omega_m/0.3)^{0.5}D_{\rm{IA}}$, where
$D_{\rm{IA}}=1-0.11(A_{\rm{IA}}-1)$. We find
$S_8D_{\rm{IA}}=0.895^{+0.054}_{-0.039}$, where systematics are at the level of
roughly 60\% of the statistical errors. We discuss these results in the context
of earlier cosmic shear analyses of the DES Y1 data. Our findings indicate that
this method and its fast runtime offer good prospects for cosmic shear
measurements with future wide-field surveys.

[7] oai:arXiv.org:1801.03181 [pdf] - 1871376

The Dark Energy Survey Data Release 1

Comments: 30 pages, 20 Figures. Release page found at this url
https://des.ncsa.illinois.edu/releases/dr1

Submitted: **2018-01-09**, last modified: **2019-04-23**

We describe the first public data release of the Dark Energy Survey, DES DR1,
consisting of reduced single epoch images, coadded images, coadded source
catalogs, and associated products and services assembled over the first three
years of DES science operations. DES DR1 is based on optical/near-infrared
imaging from 345 distinct nights (August 2013 to February 2016) by the Dark
Energy Camera mounted on the 4-m Blanco telescope at Cerro Tololo
Inter-American Observatory in Chile. We release data from the DES wide-area
survey covering ~5,000 sq. deg. of the southern Galactic cap in five broad
photometric bands, grizY. DES DR1 has a median delivered point-spread function
of g = 1.12, r = 0.96, i = 0.88, z = 0.84, and Y = 0.90 arcsec FWHM, a
photometric precision of < 1% in all bands, and an astrometric precision of 151
mas. The median coadded catalog depth for a 1.95" diameter aperture at S/N = 10
is g = 24.33, r = 24.08, i = 23.44, z = 22.69, and Y = 21.44 mag. DES DR1
includes nearly 400M distinct astronomical objects detected in ~10,000 coadd
tiles of size 0.534 sq. deg. produced from ~39,000 individual exposures.
Benchmark galaxy and stellar samples contain ~310M and ~ 80M objects,
respectively, following a basic object quality selection. These data are
accessible through a range of interfaces, including query web clients, image
cutout servers, jupyter notebooks, and an interactive coadd image visualization
tool. DES DR1 constitutes the largest photometric data set to date at the
achieved depth and photometric precision.

[8] oai:arXiv.org:1810.12186 [pdf] - 1855643

DeepSphere: Efficient spherical Convolutional Neural Network with
HEALPix sampling for cosmological applications

Comments: arXiv admin note: text overlap with arXiv:astro-ph/0409513 by other
authors

Submitted: **2018-10-29**, last modified: **2019-03-26**

Convolutional Neural Networks (CNNs) are a cornerstone of the Deep Learning
toolbox and have led to many breakthroughs in Artificial Intelligence. These
networks have mostly been developed for regular Euclidean domains such as those
supporting images, audio, or video. Because of their success, CNN-based methods
are becoming increasingly popular in Cosmology. Cosmological data often comes
as spherical maps, which make the use of the traditional CNNs more complicated.
The commonly used pixelization scheme for spherical maps is the Hierarchical
Equal Area isoLatitude Pixelisation (HEALPix). We present a spherical CNN for
analysis of full and partial HEALPix maps, which we call DeepSphere. The
spherical CNN is constructed by representing the sphere as a graph. Graphs are
versatile data structures that can act as a discrete representation of a
continuous manifold. Using the graph-based representation, we define many of
the standard CNN operations, such as convolution and pooling. With filters
restricted to being radial, our convolutions are equivariant to rotation on the
sphere, and DeepSphere can be made invariant or equivariant to rotation. This
way, DeepSphere is a special case of a graph CNN, tailored to the HEALPix
sampling of the sphere. This approach is computationally more efficient than
using spherical harmonics to perform convolutions. We demonstrate the method on
a classification problem of weak lensing mass maps from two cosmological models
and compare the performance of the CNN with that of two baseline classifiers.
The results show that the performance of DeepSphere is always superior or equal
to both of these baselines. For high noise levels and for data covering only a
smaller fraction of the sphere, DeepSphere achieves typically 10% better
classification accuracy than those baselines. Finally, we show how learned
filters can be visualized to introspect the neural network.

[9] oai:arXiv.org:1708.01530 [pdf] - 1840662

Dark Energy Survey Year 1 Results: Cosmological Constraints from Galaxy
Clustering and Weak Lensing

Comments: Matches published version. Results essentially unchanged, except
updated covariance matrix leads to improved chi^2 (colored text removed)

Submitted: **2017-08-04**, last modified: **2019-03-01**

We present cosmological results from a combined analysis of galaxy clustering
and weak gravitational lensing, using 1321 deg$^2$ of $griz$ imaging data from
the first year of the Dark Energy Survey (DES Y1). We combine three two-point
functions: (i) the cosmic shear correlation function of 26 million source
galaxies in four redshift bins, (ii) the galaxy angular autocorrelation
function of 650,000 luminous red galaxies in five redshift bins, and (iii) the
galaxy-shear cross-correlation of luminous red galaxy positions and source
galaxy shears. To demonstrate the robustness of these results, we use
independent pairs of galaxy shape, photometric redshift estimation and
validation, and likelihood analysis pipelines. To prevent confirmation bias,
the bulk of the analysis was carried out while blind to the true results; we
describe an extensive suite of systematics checks performed and passed during
this blinded phase. The data are modeled in flat $\Lambda$CDM and $w$CDM
cosmologies, marginalizing over 20 nuisance parameters, varying 6 (for
$\Lambda$CDM) or 7 (for $w$CDM) cosmological parameters including the neutrino
mass density and including the 457 $\times$ 457 element analytic covariance
matrix. We find consistent cosmological results from these three two-point
functions, and from their combination obtain $S_8 \equiv \sigma_8
(\Omega_m/0.3)^{0.5} = 0.783^{+0.021}_{-0.025}$ and $\Omega_m =
0.264^{+0.032}_{-0.019}$ for $\Lambda$CDM for $w$CDM, we find $S_8 =
0.794^{+0.029}_{-0.027}$, $\Omega_m = 0.279^{+0.043}_{-0.022}$, and
$w=-0.80^{+0.20}_{-0.22}$ at 68% CL. The precision of these DES Y1 results
rivals that from the Planck cosmic microwave background measurements, allowing
a comparison of structure in the very early and late Universe on equal terms.
Although the DES Y1 best-fit values for $S_8$ and $\Omega_m$ are lower than the
central values from Planck ...

[10] oai:arXiv.org:1802.04533 [pdf] - 1822746

Weak Lensing Analysis of SPT selected Galaxy Clusters using Dark Energy
Survey Science Verification Data

Comments: Revised version in response to referee report

Submitted: **2018-02-13**, last modified: **2018-12-12**

We present weak lensing (WL) mass constraints for a sample of massive galaxy
clusters detected by the South Pole Telescope (SPT) via the Sunyaev-Zeldovich
effect (SZE). We use $griz$ imaging data obtained from the Science Verification
(SV) phase of the Dark Energy Survey (DES) to fit the WL shear signal of 33
clusters in the redshift range $0.25 \le z \le 0.8$ with NFW profiles and to
constrain a four-parameter SPT mass-observable relation. To account for biases
in WL masses, we introduce a WL mass to true mass scaling relation described by
a mean bias and an intrinsic, log-normal scatter. We allow for correlated
scatter within the WL and SZE mass-observable relations and use simulations to
constrain priors on nuisance parameters related to bias and scatter from WL. We
constrain the normalization of the $\zeta-M_{500}$ relation,
$A_\mathrm{SZ}=12.0_{-6.7}^{+2.6}$ when using a prior on the mass slope
$B_\mathrm{SZ}$ from the latest SPT cluster cosmology analysis. Without this
prior, we recover $A_\mathrm{SZ}=10.8_{-5.2}^{+2.3}$ and
$B_\mathrm{SZ}=1.30_{-0.44}^{+0.22}$. Results in both cases imply lower cluster
masses than measured in previous work with and without WL, although the
uncertainties are large. The WL derived value of $B_\mathrm{SZ}$ is $\approx
20\%$ lower than the value preferred by the most recent SPT cluster cosmology
analysis. The method demonstrated in this work is designed to constrain cluster
masses and cosmological parameters simultaneously and will form the basis for
subsequent studies that employ the full SPT cluster sample together with the
DES data.

[11] oai:arXiv.org:1807.08732 [pdf] - 1806045

Cosmological constraints from noisy convergence maps through deep
learning

Comments: 17 pages, 12 figures

Submitted: **2018-07-23**, last modified: **2018-11-30**

Deep learning is a powerful analysis technique that has recently been
proposed as a method to constrain cosmological parameters from weak lensing
mass maps. Due to its ability to learn relevant features from the data, it is
able to extract more information from the mass maps than the commonly used
power spectrum, and thus achieve better precision for cosmological parameter
measurement. We explore the advantage of Convolutional Neural Networks (CNN)
over the power spectrum for varying levels of shape noise and different
smoothing scales applied to the maps. We compare the cosmological constraints
from the two methods in the $\Omega_M-\sigma_8$ plane for sets of 400 deg$^2$
convergence maps. We find that, for a shape noise level corresponding to 8.53
galaxies/arcmin$^2$ and the smoothing scale of $\sigma_s = 2.34$ arcmin, the
network is able to generate 45% tighter constraints. For smaller smoothing
scale of $\sigma_s = 1.17$ the improvement can reach $\sim 50 \%$, while for
larger smoothing scale of $\sigma_s = 5.85$, the improvement decreases to 19%.
The advantage generally decreases when the noise level and smoothing scales
increase. We present a new training strategy to train the neural network with
noisy data, as well as considerations for practical applications of the deep
learning approach.

[12] oai:arXiv.org:1801.09070 [pdf] - 1790606

Fast cosmic web simulations with generative adversarial networks

Comments:

Submitted: **2018-01-27**, last modified: **2018-11-29**

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
N-body 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 N-body simulations of size 500 and 100 Mpc.
We show that the GAN-generated samples are qualitatively and quantitatively
very similar to the originals. For the larger boxes of size 500 Mpc, it is very
difficult to distinguish them visually. The agreement of the power spectrum
$P_k$ is 1-2\% for most of the range, between $k=0.06$ and $k=0.4$. An
important advantage of generating cosmic web realizations with a GAN is the
considerable gains in terms of computation time. Each new sample generated by a
GAN takes a fraction of a second, compared to the many hours needed by
traditional N-body techniques. We anticipate that the use of generative models
such as 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 Large Synoptic Survey Telescope (LSST).

[13] oai:arXiv.org:1803.08461 [pdf] - 1776930

Weak lensing peak statistics in the era of large scale cosmological
surveys

Comments: 17 pages, 9 figures

Submitted: **2018-03-22**, last modified: **2018-10-31**

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 Comoving-Lagrangian 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.

[14] oai:arXiv.org:1801.08945 [pdf] - 1761469

Improving Weak Lensing Mass Map Reconstructions using Gaussian and
Sparsity Priors: Application to DES SV

Comments: 19 pages, 10 figures, MNRAS published: 15 May 2018

Submitted: **2018-01-26**, last modified: **2018-09-26**

Mapping the underlying density field, including non-visible dark matter,
using weak gravitational lensing measurements is now a standard tool in
cosmology. Due to its importance to the science results of current and upcoming
surveys, the quality of the convergence reconstruction methods should be well
understood. We compare three methods: Kaiser-Squires (KS), Wiener filter, and
GLIMPSE. KS is a direct inversion, not accounting for survey masks or noise.
The Wiener filter is well-motivated for Gaussian density fields in a Bayesian
framework. GLIMPSE uses sparsity, aiming to reconstruct non-linearities in the
density field. We compare these methods with several tests using public Dark
Energy Survey (DES) Science Verification (SV) data and realistic DES
simulations. The Wiener filter and GLIMPSE offer substantial improvements over
smoothed KS with a range of metrics. Both the Wiener filter and GLIMPSE
convergence reconstructions show a 12 per cent improvement in Pearson
correlation with the underlying truth from simulations. To compare the mapping
methods' abilities to find mass peaks, we measure the difference between peak
counts from simulated {\Lambda}CDM shear catalogues and catalogues with no mass
fluctuations (a standard data vector when inferring cosmology from peak
statistics); the maximum signal-to-noise of these peak statistics is increased
by a factor of 3.5 for the Wiener filter and 9 for GLIMPSE. With simulations we
measure the reconstruction of the harmonic phases; the phase residuals'
concentration is improved 17 per cent by GLIMPSE and 18 per cent by the Wiener
filter. The correlation between reconstructions from data and foreground
redMaPPer clusters is increased 18 per cent by the Wiener filter and 32 per
cent by GLIMPSE.

[15] oai:arXiv.org:1708.01533 [pdf] - 1746405

Dark Energy Survey Year 1 Results: Weak Lensing Shape Catalogues

Comments: 36 Pages, 29 figures. This version accepted for publication in MNRAS

Submitted: **2017-08-04**, last modified: **2018-09-07**

We present two galaxy shape catalogues from the Dark Energy Survey Year 1
data set, covering 1500 square degrees with a median redshift of $0.59$. The
catalogues cover two main fields: Stripe 82, and an area overlapping the South
Pole Telescope survey region. We describe our data analysis process and in
particular our shape measurement using two independent shear measurement
pipelines, METACALIBRATION and IM3SHAPE. The METACALIBRATION catalogue uses a
Gaussian model with an innovative internal calibration scheme, and was applied
to $riz$-bands, yielding 34.8M objects. The IM3SHAPE catalogue uses a
maximum-likelihood bulge/disc model calibrated using simulations, and was
applied to $r$-band data, yielding 21.9M objects. Both catalogues pass a suite
of null tests that demonstrate their fitness for use in weak lensing science.
We estimate the 1$\sigma$ uncertainties in multiplicative shear calibration to
be $0.013$ and $0.025$ for the METACALIBRATION and IM3SHAPE catalogues,
respectively.

[16] oai:arXiv.org:1801.07615 [pdf] - 1722450

Fast Point Spread Function Modeling with Deep Learning

Comments: 25 pages, 8 figures, 1 table

Submitted: **2018-01-23**, last modified: **2018-07-25**

Modeling the Point Spread Function (PSF) of wide-field 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 Monte-Carlo Control Loops ($MCCL$)
framework have been developed. These methods typically require simulating a
large amount of wide-field 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
deep-learning 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.

[17] oai:arXiv.org:1708.01538 [pdf] - 1747812

Dark Energy Survey Year 1 Results: Cosmological Constraints from Cosmic
Shear

Comments: 32 pages, 19 figures; matches PRD referee response version

Submitted: **2017-08-04**, last modified: **2018-04-30**

We use 26 million galaxies from the Dark Energy Survey (DES) Year 1 shape
catalogs over 1321 deg$^2$ of the sky to produce the most significant
measurement of cosmic shear in a galaxy survey to date. We constrain
cosmological parameters in both the flat $\Lambda$CDM and $w$CDM models, while
also varying the neutrino mass density. These results are shown to be robust
using two independent shape catalogs, two independent \photoz\ calibration
methods, and two independent analysis pipelines in a blind analysis. We find a
3.5\% fractional uncertainty on $\sigma_8(\Omega_m/0.3)^{0.5} =
0.782^{+0.027}_{-0.027}$ at 68\% CL, which is a factor of 2.5 improvement over
the fractional constraining power of our DES Science Verification results. In
$w$CDM, we find a 4.8\% fractional uncertainty on $\sigma_8(\Omega_m/0.3)^{0.5}
= 0.777^{+0.036}_{-0.038}$ and a dark energy equation-of-state
$w=-0.95^{+0.33}_{-0.39}$. We find results that are consistent with previous
cosmic shear constraints in $\sigma_8$ -- $\Omega_m$, and see no evidence for
disagreement of our weak lensing data with data from the CMB. Finally, we find
no evidence preferring a $w$CDM model allowing $w\ne -1$. We expect further
significant improvements with subsequent years of DES data, which will more
than triple the sky coverage of our shape catalogs and double the effective
integrated exposure time per galaxy.

[18] oai:arXiv.org:1708.01535 [pdf] - 1637501

Dark Energy Survey Year 1 Results: Curved-Sky Weak Lensing Mass Map

Comments: 25 pages, 19 figures, 1 table; revision with changes implemented
according to journal referee

Submitted: **2017-08-04**, last modified: **2017-12-19**

We construct the largest curved-sky galaxy weak lensing mass map to date from
the DES first-year (DES Y1) data. The map, about 10 times larger than previous
work, is constructed over a contiguous $\approx1,500 $deg$^2$, covering a
comoving volume of $\approx10 $Gpc$^3$. The effects of masking, sampling, and
noise are tested using simulations. We generate weak lensing maps from two DES
Y1 shear catalogs, Metacalibration and Im3shape, with sources at redshift
$0.2<z<1.3,$ and in each of four bins in this range. In the highest
signal-to-noise map, the ratio between the mean signal-to-noise in the E-mode
and the B-mode map is $\sim$1.5 ($\sim$2) when smoothed with a Gaussian filter
of $\sigma_{G}=30$ (80) arcminutes. The second and third moments of the
convergence $\kappa$ in the maps are in agreement with simulations. We also
find no significant correlation of $\kappa$ with maps of potential systematic
contaminants. Finally, we demonstrate two applications of the mass maps: (1)
cross-correlation with different foreground tracers of mass and (2) exploration
of the largest peaks and voids in the maps.

[19] oai:arXiv.org:1705.05386 [pdf] - 1583439

The redshift distribution of cosmological samples: a forward modeling
approach

Comments: 26 pages, 10 figures

Submitted: **2017-05-15**, last modified: **2017-08-11**

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
object-by-object 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$ (Monte-Carlo 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
4-band Subaru Suprime-Cam 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.

[20] oai:arXiv.org:1708.01534 [pdf] - 1648606

Dark Energy Survey Year 1 Results: The Impact of Galaxy Neighbours on
Weak Lensing Cosmology with im3shape

Comments: 20 pages, submitted to MNRAS

Submitted: **2017-08-04**

We use a suite of simulated images based on Year 1 of the Dark Energy Survey
to explore the impact of galaxy neighbours on shape measurement and shear
cosmology. The hoopoe image simulations include realistic blending, galaxy
positions, and spatial variations in depth and PSF properties. Using the
im3shape maximum-likelihood shape measurement code, we identify four mechanisms
by which neighbours can have a non-negligible influence on shear estimation.
These effects, if ignored, would contribute a net multiplicative bias of $m
\sim 0.03 - 0.09$ in the DES Y1 im3shape catalogue, though the precise impact
will be dependent on both the measurement code and the selection cuts applied.
This can be reduced to percentage level or less by removing objects with close
neighbours, at a cost to the effective number density of galaxies
$n_\mathrm{eff}$ of 30%. We use the cosmological inference pipeline of DES Y1
to explore the cosmological implications of neighbour bias and show that
omitting blending from the calibration simulation for DES Y1 would bias the
inferred clustering amplitude $S_8\equiv \sigma_8 (\Omega _\mathrm{m}
/0.3)^{0.5}$ by $2 \sigma$ towards low values. Finally, we use the hoopoe
simulations to test the effect of neighbour-induced spatial correlations in the
multiplicative bias. We find the impact on the recovered $S_8$ of ignoring such
correlations to be subdominant to statistical error at the current level of
precision.

[21] oai:arXiv.org:1707.07498 [pdf] - 1641256

Accelerating Approximate Bayesian Computation with Quantile Regression:
Application to Cosmological Redshift Distributions

Comments: 10 pages, 5 figures, prepared for submission to JCAP

Submitted: **2017-07-24**, last modified: **2017-07-25**

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.

[22] oai:arXiv.org:1707.06233 [pdf] - 1732528

Cosmic shear calibration with forward modeling

Comments: 23 pages, 4 figures, prepared for submission to JCAP

Submitted: **2017-07-19**

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 one-point statistics to calibrate the
shear measurement and six further one-point statistics as diagnostics. We also
propagate the systematics levels that we estimate from the one-point functions
to the two-point 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.

[23] oai:arXiv.org:1707.05167 [pdf] - 1586002

Cosmological model discrimination with Deep Learning

Comments: 21 pages, 9 figures, prepared for submission to JCAP

Submitted: **2017-07-17**, last modified: **2017-07-18**

We demonstrate the potential of Deep Learning methods for measurements of
cosmological parameters from density fields, focusing on the extraction of
non-Gaussian 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 two-point 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 non-Gaussian 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 CNN-based 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.

[24] oai:arXiv.org:1507.05552 [pdf] - 1572642

Cosmology from Cosmic Shear with DES Science Verification Data

Comments: Finally updating to the published version. 20 pages, 12 figures.
Additional information at http://deswl.github.io/

Submitted: **2015-07-20**, last modified: **2017-05-03**

We present the first constraints on cosmology from the Dark Energy Survey
(DES), using weak lensing measurements from the preliminary Science
Verification (SV) data. We use 139 square degrees of SV data, which is less
than 3\% of the full DES survey area. Using cosmic shear 2-point measurements
over three redshift bins we find $\sigma_8 (\Omega_{\rm m}/0.3)^{0.5} = 0.81
\pm 0.06$ (68\% confidence), after marginalising over 7 systematics parameters
and 3 other cosmological parameters. We examine the robustness of our results
to the choice of data vector and systematics assumed, and find them to be
stable. About $20$\% of our error bar comes from marginalising over shear and
photometric redshift calibration uncertainties. The current state-of-the-art
cosmic shear measurements from CFHTLenS are mildly discrepant with the
cosmological constraints from Planck CMB data; our results are consistent with
both datasets. Our uncertainties are $\sim$30\% larger than those from CFHTLenS
when we carry out a comparable analysis of the two datasets, which we attribute
largely to the lower number density of our shear catalogue. We investigate
constraints on dark energy and find that, with this small fraction of the full
survey, the DES SV constraints make negligible impact on the Planck
constraints. The moderate disagreement between the CFHTLenS and Planck values
of $\sigma_8 (\Omega_{\rm m}/0.3)^{0.5}$ is present regardless of the value of
$w$.

[25] oai:arXiv.org:1610.06890 [pdf] - 1580445

Weak-lensing mass calibration of redMaPPer galaxy clusters in Dark
Energy Survey Science Verification data

Comments: 22 pages, 14 figures, accepted by MNRAS

Submitted: **2016-10-21**, last modified: **2017-04-28**

We use weak-lensing shear measurements to determine the mean mass of
optically selected galaxy clusters in Dark Energy Survey Science Verification
data. In a blinded analysis, we split the sample of more than 8,000 redMaPPer
clusters into 15 subsets, spanning ranges in the richness parameter $5 \leq
\lambda \leq 180$ and redshift $0.2 \leq z \leq 0.8$, and fit the averaged mass
density contrast profiles with a model that accounts for seven distinct sources
of systematic uncertainty: shear measurement and photometric redshift errors;
cluster-member contamination; miscentering; deviations from the NFW halo
profile; halo triaxiality; and line-of-sight projections. We combine the
inferred cluster masses to estimate the joint scaling relation between mass,
richness and redshift, $\mathcal{M}(\lambda,z) \varpropto M_0 \lambda^{F}
(1+z)^{G}$. We find $M_0 \equiv \langle
M_{200\mathrm{m}}\,|\,\lambda=30,z=0.5\rangle=\left[ 2.35 \pm 0.22\ \rm{(stat)}
\pm 0.12\ \rm{(sys)} \right] \cdot 10^{14}\ M_\odot$, with $F =
1.12\,\pm\,0.20\ \rm{(stat)}\, \pm\, 0.06\ \rm{(sys)}$ and $G = 0.18\,\pm\,
0.75\ \rm{(stat)}\, \pm\, 0.24\ \rm{(sys)}$. The amplitude of the mass-richness
relation is in excellent agreement with the weak-lensing calibration of
redMaPPer clusters in SDSS by Simet et al. (2016) and with the Saro et al.
(2015) calibration based on abundance matching of SPT-detected clusters. Our
results extend the redshift range over which the mass-richness relation of
redMaPPer clusters has been calibrated with weak lensing from $z\leq 0.3$ to
$z\leq0.8$. Calibration uncertainties of shear measurements and photometric
redshift estimates dominate our systematic error budget and require substantial
improvements for forthcoming studies.

[26] oai:arXiv.org:1608.00004 [pdf] - 1580104

Searching for Cosmic Strings in CMB Anisotropy Maps using Wavelets and
Curvelets

Comments: 13 pages, 7 figures, 3 tables

Submitted: **2016-07-29**, last modified: **2016-09-28**

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 (SPT-3G) survey will be able to yield stronger limits. For
maps with a specification of SPT-3G we find that string signals will be visible
down to a string tension of $G \mu = 1.4 \times 10^{-7}$.

[27] oai:arXiv.org:1608.01838 [pdf] - 1531199

Inference from the small scales of cosmic shear with current and future
Dark Energy Survey data

Comments: 14 pages, 10 figures

Submitted: **2016-08-05**

Cosmic shear is sensitive to fluctuations in the cosmological matter density
field, including on small physical scales, where matter clustering is affected
by baryonic physics in galaxies and galaxy clusters, such as star formation,
supernovae feedback and AGN feedback. While muddying any cosmological
information that is contained in small scale cosmic shear measurements, this
does mean that cosmic shear has the potential to constrain baryonic physics and
galaxy formation. We perform an analysis of the Dark Energy Survey (DES)
Science Verification (SV) cosmic shear measurements, now extended to smaller
scales, and using the Mead et al. 2015 halo model to account for baryonic
feedback. While the SV data has limited statistical power, we demonstrate using
a simulated likelihood analysis that the final DES data will have the
statistical power to differentiate among baryonic feedback scenarios. We also
explore some of the difficulties in interpreting the small scales in cosmic
shear measurements, presenting estimates of the size of several other
systematic effects that make inference from small scales difficult, including
uncertainty in the modelling of intrinsic alignment on nonlinear scales,
`lensing bias', and shape measurement selection effects. For the latter two, we
make use of novel image simulations. While future cosmic shear datasets have
the statistical power to constrain baryonic feedback scenarios, there are
several systematic effects that require improved treatments, in order to make
robust conclusions about baryonic feedback.

[28] oai:arXiv.org:1507.05598 [pdf] - 1445264

Cosmic Shear Measurements with DES Science Verification Data

Comments: measurements and covariance matrices in machine readable format are
available as ancillary data on the arXiv; high-resolution versions of figures
can be downloaded from http://deswl.github.io

Submitted: **2015-07-20**, last modified: **2016-07-27**

We present measurements of weak gravitational lensing cosmic shear two-point
statistics using Dark Energy Survey Science Verification data. We demonstrate
that our results are robust to the choice of shear measurement pipeline, either
ngmix or im3shape, and robust to the choice of two-point statistic, including
both real and Fourier-space statistics. Our results pass a suite of null tests
including tests for B-mode contamination and direct tests for any dependence of
the two-point functions on a set of 16 observing conditions and galaxy
properties, such as seeing, airmass, galaxy color, galaxy magnitude, etc. We
furthermore use a large suite of simulations to compute the covariance matrix
of the cosmic shear measurements and assign statistical significance to our
null tests. We find that our covariance matrix is consistent with the halo
model prediction, indicating that it has the appropriate level of halo sample
variance. We compare the same jackknife procedure applied to the data and the
simulations in order to search for additional sources of noise not captured by
the simulations. We find no statistically significant extra sources of noise in
the data. The overall detection significance with tomography for our highest
source density catalog is 9.7sigma. Cosmological constraints from the
measurements in this work are presented in a companion paper (DES et al. 2015).

[29] oai:arXiv.org:1601.00405 [pdf] - 1415044

Galaxy bias from the Dark Energy Survey Science Verification data:
combining galaxy density maps and weak lensing maps

Comments: 14 pages, 9 figures, 2 table; to be published in MNRAS; update
title/affiliation to match journal version

Submitted: **2016-01-04**, last modified: **2016-04-27**

We measure the redshift evolution of galaxy bias for a magnitude-limited
galaxy sample by combining the galaxy density maps and weak lensing shear maps
for a $\sim$116 deg$^{2}$ area of the Dark Energy Survey (DES) Science
Verification data. This method was first developed in Amara et al. (2012) and
later re-examined in a companion paper (Pujol et al. 2016) with rigorous
simulation tests and analytical treatment of tomographic measurements. In this
work we apply this method to the DES SV data and measure the galaxy bias for a
i$<$22.5 galaxy sample. We find the galaxy bias and 1$\sigma$ error bars in 4
photometric redshift bins to be 1.12$\pm$0.19 (z=0.2-0.4), 0.97$\pm$0.15
(z=0.4-0.6), 1.38$\pm$0.39 (z=0.6-0.8)), and 1.45$\pm$0.56 (z=0.8-1.0). These
measurements are consistent at the 2$\sigma$ level with measurements on the
same dataset using galaxy clustering and cross-correlation of galaxies with CMB
lensing, with most of the redshift bins consistent within the 1{\sigma} error
bars. In addition, our method provides the only $\sigma_8$-independent
constraint among the three. We forward-model the main observational effects
using mock galaxy catalogs by including shape noise, photo-z errors and masking
effects. We show that our bias measurement from the data is consistent with
that expected from simulations. With the forthcoming full DES data set, we
expect this method to provide additional constraints on the galaxy bias
measurement from more traditional methods. Furthermore, in the process of our
measurement, we build up a 3D mass map that allows further exploration of the
dark matter distribution and its relation to galaxy evolution.

[30] oai:arXiv.org:1604.07871 [pdf] - 1494677

Cosmology from large scale galaxy clustering and galaxy-galaxy lensing
with Dark Energy Survey Science Verification data

Comments: 18 pages, 13 figures, submitted to MNRAS

Submitted: **2016-04-26**

We present cosmological constraints from the Dark Energy Survey (DES) using a
combined analysis of angular clustering of red galaxies and their
cross-correlation 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 galaxy-galaxy lensing with Planck CMB data, Baryon Accoustic
Oscillations and Supernova type Ia measurements.

[31] oai:arXiv.org:1507.05603 [pdf] - 1411151

The DES Science Verification Weak Lensing Shear Catalogues

Comments: Accepted by MNRAS; 38 pages, 29 figures; v3: minor edits based on
referee's comments, switched to mnras style, added figure 8, updated info
about released catalogs

Submitted: **2015-07-20**, last modified: **2016-04-25**

We present weak lensing shear catalogues for 139 square degrees of data taken
during the Science Verification (SV) time for the new Dark Energy Camera
(DECam) being used for the Dark Energy Survey (DES). We describe our object
selection, point spread function estimation and shear measurement procedures
using two independent shear pipelines, IM3SHAPE and NGMIX, which produce
catalogues of 2.12 million and 3.44 million galaxies respectively. We detail a
set of null tests for the shear measurements and find that they pass the
requirements for systematic errors at the level necessary for weak lensing
science applications using the SV data. We also discuss some of the planned
algorithmic improvements that will be necessary to produce sufficiently
accurate shear catalogues for the full 5-year DES, which is expected to cover
5000 square degrees.

[32] oai:arXiv.org:1603.05040 [pdf] - 1470597

Cosmology constraints from shear peak statistics in Dark Energy Survey
Science Verification data

Comments: 21 pages, 14 figures, submitted to MNRAS

Submitted: **2016-03-16**

Shear peak statistics has gained a lot of attention recently as a practical
alternative to the two point statistics for constraining cosmological
parameters. We perform a shear peak statistics analysis of the Dark Energy
Survey (DES) Science Verification (SV) data, using weak gravitational lensing
measurements from a 139 deg$^2$ field. We measure the abundance of peaks
identified in aperture mass maps, as a function of their signal-to-noise ratio,
in the signal-to-noise range $0<\mathcal S / \mathcal N<4$. To predict the peak
counts as a function of cosmological parameters we use a suite of $N$-body
simulations spanning 158 models with varying $\Omega_{\rm m}$ and $\sigma_8$,
fixing $w = -1$, $\Omega_{\rm b} = 0.04$, $h = 0.7$ and $n_s=1$, to which we
have applied the DES SV mask and redshift distribution. In our fiducial
analysis we measure $\sigma_{8}(\Omega_{\rm m}/0.3)^{0.6}=0.77 \pm 0.07$, after
marginalising over the shear multiplicative bias and the error on the mean
redshift of the galaxy sample. We introduce models of intrinsic alignments,
blending, and source contamination by cluster members. These models indicate
that peaks with $\mathcal S / \mathcal N>4$ would require significant
corrections, which is why we do not include them in our analysis. We compare
our results to the cosmological constraints from the two point analysis on the
SV field and find them to be in good agreement in both the central value and
its uncertainty. We discuss prospects for future peak statistics analysis with
upcoming DES data.

[33] oai:arXiv.org:1512.04535 [pdf] - 1378743

Cross-correlation of gravitational lensing from DES Science Verification
data with SPT and Planck lensing

Comments: 13 pages, 6 figures

Submitted: **2015-12-14**

We measure the cross-correlation between weak lensing of galaxy images and of
the cosmic microwave background (CMB). The effects of gravitational lensing on
different sources will be correlated if the lensing is caused by the same mass
fluctuations. We use galaxy shape measurements from 139 deg$^{2}$ of the Dark
Energy Survey (DES) Science Verification data and overlapping CMB lensing from
the South Pole Telescope (SPT) and Planck. The DES source galaxies have a
median redshift of $z_{\rm med} {\sim} 0.7$, while the CMB lensing kernel is
broad and peaks at $z{\sim}2$. The resulting cross-correlation is maximally
sensitive to mass fluctuations at $z{\sim}0.44$. Assuming the Planck 2015
best-fit cosmology, the amplitude of the DES$\times$SPT cross-power is found to
be $A = 0.88 \pm 0.30$ and that from DES$\times$Planck to be $A = 0.86 \pm
0.39$, where $A=1$ corresponds to the theoretical prediction. These are
consistent with the expected signal and correspond to significances of $2.9
\sigma$ and $2.2 \sigma$ respectively. We demonstrate that our results are
robust to a number of important systematic effects including the shear
measurement method, estimator choice, photometric redshift uncertainty and CMB
lensing systematics. Significant intrinsic alignment of galaxy shapes would
increase the cross-correlation signal inferred from the data; we calculate a
value of $A = 1.08 \pm 0.36$ for DES$\times$SPT when we correct the
observations with a simple IA model. With three measurements of this
cross-correlation now existing in the literature, there is not yet reliable
evidence for any deviation from the expected LCDM level of cross-correlation,
given the size of the statistical uncertainties and the significant impact of
systematic errors, particularly IAs. We provide forecasts for the expected
signal-to-noise of the combination of the five-year DES survey and SPT-3G.

[34] oai:arXiv.org:1507.05090 [pdf] - 1324854

Weak lensing by galaxy troughs in DES Science Verification data

Comments: 15 pages, 8 figures; matches accepted version; high-resolution
versions of figures can be downloaded from http://deswl.github.io

Submitted: **2015-07-17**, last modified: **2015-12-08**

We measure the weak lensing shear around galaxy troughs, i.e. the radial
alignment of background galaxies relative to underdensities in projections of
the foreground galaxy field over a wide range of redshift in Science
Verification data from the Dark Energy Survey. Our detection of the shear
signal is highly significant (10 to 15$\sigma$ for the smallest angular scales)
for troughs with the redshift range z in [0.2,0.5] of the projected galaxy
field and angular diameters of 10 arcmin...1{\deg}. These measurements probe
the connection between the galaxy, matter density, and convergence fields. By
assuming galaxies are biased tracers of the matter density with Poissonian
noise, we find agreement of our measurements with predictions in a fiducial
Lambda cold dark matter model. The prediction for the lensing signal on large
trough scales is virtually independent of the details of the underlying model
for the connection of galaxies and matter. Our comparison of the shear around
troughs with that around cylinders with large galaxy counts is consistent with
a symmetry between galaxy and matter over- and underdensities. In addition, we
measure the two-point angular correlation of troughs with galaxies which, in
contrast to the lensing signal, is sensitive to galaxy bias on all scales. The
lensing signal of troughs and their clustering with galaxies is therefore a
promising probe of the statistical properties of matter underdensities and
their connection to the galaxy field.

[35] oai:arXiv.org:1507.05909 [pdf] - 1475231

Redshift distributions of galaxies in the DES Science Verification shear
catalogue and implications for weak lensing

Comments: high-resolution versions of figures can be downloaded from
http://deswl.github.io

Submitted: **2015-07-21**, last modified: **2015-07-23**

We present photometric redshift estimates for galaxies used in the weak
lensing analysis of the Dark Energy Survey Science Verification (DES SV) data.
Four model- or machine learning-based photometric redshift methods -- ANNZ2,
BPZ calibrated against BCC-Ufig simulations, SkyNet, and TPZ -- are analysed.
For training, calibration, and testing of these methods, we construct a
catalogue of spectroscopically confirmed galaxies matched against DES SV data.
The performance of the methods is evaluated against the matched spectroscopic
catalogue, focusing on metrics relevant for weak lensing analyses, with
additional validation against COSMOS photo-zs. From the galaxies in the DES SV
shear catalogue, which have mean redshift $0.72\pm0.01$ over the range
$0.3<z<1.3$, we construct three tomographic bins with means of $z=\{0.45, 0.67,
1.00\}$. These bins each have systematic uncertainties $\delta z \lesssim 0.05$
in the mean of the fiducial SkyNet photo-z $n(z)$. We propagate the errors in
the redshift distributions through to their impact on cosmological parameters
estimated with cosmic shear, and find that they cause shifts in the value of
$\sigma_8$ of approx. 3%. This shift is within the one sigma statistical errors
on $\sigma_8$ for the DES SV shear catalog. We further study the potential
impact of systematic differences on the critical surface density,
$\Sigma_{\mathrm{crit}}$, finding levels of bias safely less than the
statistical power of DES SV data. We recommend a final Gaussian prior for the
photo-z bias in the mean of $n(z)$ of width $0.05$ for each of the three
tomographic bins, and show that this is a sufficient bias model for the
corresponding cosmology analysis.

[36] oai:arXiv.org:1505.01871 [pdf] - 1255013

Wide-Field Lensing Mass Maps from DES Science Verification Data

Comments: 7 pages, 3 figures, 1 table; accepted to PRL

Submitted: **2015-05-07**, last modified: **2015-07-20**

We present a mass map reconstructed from weak gravitational lensing shear
measurements over 139 sq. deg from the Dark Energy Survey (DES) Science
Verification data. The mass map probes both luminous and dark matter, thus
providing a tool for studying cosmology. We find good agreement between the
mass map and the distribution of massive galaxy clusters identified using a
red-sequence cluster finder. Potential candidates for super-clusters and voids
are identified using these maps. We measure the cross-correlation between the
mass map and a magnitude-limited foreground galaxy sample and find a detection
at the 5-7 sigma level on a large range of scales. These measurements are
consistent with simulated galaxy catalogs based on LCDM N-body simulations,
suggesting low systematics uncertainties in the map. We summarize our key
findings in this letter; the detailed methodology and tests for systematics are
presented in a companion paper.

[37] oai:arXiv.org:1504.03002 [pdf] - 1248536

Wide-Field Lensing Mass Maps from DES Science Verification Data:
Methodology and Detailed Analysis

Comments: 21 pages, 14 figures, 2 tables; accepted to PRD

Submitted: **2015-04-12**, last modified: **2015-07-20**

Weak gravitational lensing allows one to reconstruct the spatial distribution
of the projected mass density across the sky. These "mass maps" provide a
powerful tool for studying cosmology as they probe both luminous and dark
matter. In this paper, we present a weak lensing mass map reconstructed from
shear measurements in a 139 sq. deg area from the Dark Energy Survey (DES)
Science Verification (SV) data. We compare the distribution of mass with that
of the foreground distribution of galaxies and clusters. The overdensities in
the reconstructed map correlate well with the distribution of optically
detected clusters. We demonstrate that candidate superclusters and voids along
the line of sight can be identified, exploiting the tight scatter of the
cluster photometric redshifts. We cross-correlate the mass map with a
foreground magnitude-limited galaxy sample from the same data. Our measurement
gives results consistent with mock catalogs from N-body simulations that
include the primary sources of statistical uncertainties in the galaxy,
lensing, and photo-z catalogs. The statistical significance of the
cross-correlation is at the 6.8-sigma level with 20 arcminute smoothing. A
major goal of this study is to investigate systematic effects arising from a
variety of sources, including PSF and photo-z uncertainties. We make maps
derived from twenty variables that may characterize systematics and find the
principal components. We find that the contribution of systematics to the
lensing mass maps is generally within measurement uncertainties. In this work,
we analyze less than 3% of the final area that will be mapped by the DES; the
tools and analysis techniques developed in this paper can be applied to
forthcoming larger datasets from the survey.

[38] oai:arXiv.org:1412.1825 [pdf] - 1005617

GREAT3 results I: systematic errors in shear estimation and the impact
of real galaxy morphology

Comments: 32 pages + 15 pages of technical appendices; 28 figures; submitted to
MNRAS; latest version has minor updates in presentation of 4 figures, no
changes in content or conclusions

Submitted: **2014-12-04**, last modified: **2015-04-02**

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 ground-based data with
constant or cosmologically-varying shear fields. The simplest (control)
experiment included parametric galaxies with a realistic distribution of
signal-to-noise, 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
spatially-varying 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 Stage-IV 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 signal-to-noise, 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.

[39] oai:arXiv.org:1405.4285 [pdf] - 996594

Mass and galaxy distributions of four massive galaxy clusters from Dark
Energy Survey Science Verification data

Comments: accepted by MNRAS; high-resolution versions of figures can be
downloaded from
https://www.physics.ohio-state.edu/~melchior.12/research/des_sv_clusters.html

Submitted: **2014-05-16**, last modified: **2015-02-28**

We measure the weak-lensing masses and galaxy distributions of four massive
galaxy clusters observed during the Science Verification phase of the Dark
Energy Survey. This pathfinder study is meant to 1) validate the DECam imager
for the task of measuring weak-lensing shapes, and 2) utilize DECam's large
field of view to map out the clusters and their environments over 90 arcmin. We
conduct a series of rigorous tests on astrometry, photometry, image quality,
PSF modeling, and shear measurement accuracy to single out flaws in the data
and also to identify the optimal data processing steps and parameters. We find
Science Verification data from DECam to be suitable for the lensing analysis
described in this paper. The PSF is generally well-behaved, but the modeling is
rendered difficult by a flux-dependent PSF width and ellipticity. We employ
photometric redshifts to distinguish between foreground and background
galaxies, and a red-sequence cluster finder to provide cluster richness
estimates and cluster-galaxy distributions. By fitting NFW profiles to the
clusters in this study, we determine weak-lensing masses that are in agreement
with previous work. For Abell 3261, we provide the first estimates of redshift,
weak-lensing mass, and richness. In addition, the cluster-galaxy distributions
indicate the presence of filamentary structures attached to 1E 0657-56 and RXC
J2248.7-4431, stretching out as far as 1 degree (approximately 20 Mpc),
showcasing the potential of DECam and DES for detailed studies of degree-scale
features on the sky.

[40] oai:arXiv.org:1407.7676 [pdf] - 935774

GalSim: The modular galaxy image simulation toolkit

Comments: 38 pages, 3 tables, 8 figures, accepted for publication in Astronomy
& Computing

Submitted: **2014-07-29**, last modified: **2015-02-15**

GALSIM is a collaborative, open-source 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 WFIRST-AFTA 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/GalSim-developers/GalSim .

[41] oai:arXiv.org:1308.4982 [pdf] - 1178725

The Third Gravitational Lensing Accuracy Testing (GREAT3) Challenge
Handbook

Comments: 30 pages, 13 figures, submitted for publication, with minor edits
(v2) to address comments from the anonymous referee. Simulated data are
available for download and participants can find more information at
http://great3.projects.phys.ucl.ac.uk/leaderboard/

Submitted: **2013-08-22**, last modified: **2014-02-28**

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 high-resolution imaging from
space; spatially varying, physically-motivated 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.

[42] oai:arXiv.org:1304.6673 [pdf] - 765729

Gaussian Process for star and planet characterisation

Comments: 31 pages, 9 figures, 2 tables, submitted to ApJ

Submitted: **2013-04-24**, last modified: **2014-01-02**

The study of exoplanetary atmospheres epitomises a continuous quest for
higher accuracy measurements. Systematic effects and noise associated with both
the stellar activity and the instrument can bias the results and thus limit the
precision of the analysis. To reach a high photometric and spectroscopic
precision, it is therefore essential to correct for these effects. We present
here a novel non-parametric approach, named Gaussian Process method for Star
Characterization (GPSC), to remove effects of stellar activity and instrumental
systematics on planetary signals, with a view to preserve the atmospheric
contribution which can be as small as 10$^{-4}$ or even 10$^{-5}$ the flux of
the star. We applied our method to data recorded with Kepler, focussing on a
sample of lightcurves with different effective temperatures and flux
modulations. We found that GPSC can very effectively correct for the short and
long term stellar activity and instrumental systematics. Additionally we run
the GPSC on both real and simulated transit data, finding transit depths
consistent with the original ones. Consequently we considered 10 hours of
continuous observations: daily, every other day and weekly, and we used the
GPSC to reconstruct the lightcurves. When data are recorded more frequently
than once every five days we found that our approach is able to extrapolate the
stellar flux at the 10$^{-4}$ level compared to the full stellar flux. These
results show a great potential of GPSC to isolate the relevant astrophysical
signal and achieve the precision needed for the correction of short and long
term stellar activity.

[43] oai:arXiv.org:1308.4663 [pdf] - 1542622

S\'{e}rsic galaxy models in weak lensing shape measurement: model bias,
noise bias and their interaction

Comments: 11 pages, 3 figures

Submitted: **2013-08-21**

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 non-linear
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.

[44] oai:arXiv.org:1202.5254 [pdf] - 630580

Image Analysis for Cosmology: Results from the GREAT10 Galaxy Challenge

Comments: Accepted to MNRAS

Submitted: **2012-02-23**, last modified: **2013-02-25**

In this paper we present results from the weak lensing shape measurement
GRavitational lEnsing Accuracy Testing 2010 (GREAT10) Galaxy Challenge. This
marks an order of magnitude step change in the level of scrutiny employed in
weak lensing shape measurement analysis. We provide descriptions of each method
tested and include 10 evaluation metrics over 24 simulation branches. GREAT10
was the first shape measurement challenge to include variable fields; both the
shear field and the Point Spread Function (PSF) vary across the images in a
realistic manner. The variable fields enable a variety of metrics that are
inaccessible to constant shear simulations including a direct measure of the
impact of shape measurement inaccuracies, and the impact of PSF size and
ellipticity, on the shear power spectrum. To assess the impact of shape
measurement bias for cosmic shear we present a general pseudo-Cl formalism,
that propagates spatially varying systematics in cosmic shear through to power
spectrum estimates. We also show how one-point estimators of bias can be
extracted from variable shear simulations. The GREAT10 Galaxy Challenge
received 95 submissions and saw a factor of 3 improvement in the accuracy
achieved by shape measurement methods. The best methods achieve sub-percent
average biases. We find a strong dependence in accuracy as a function of
signal-to-noise, and indications of a weak dependence on galaxy type and size.
Some requirements for the most ambitious cosmic shear experiments are met above
a signal-to-noise ratio of 20. These results have the caveat that the simulated
PSF was a ground-based PSF. Our results are a snapshot of the accuracy of
current shape measurement methods and are a benchmark upon which improvement
can continue. This provides a foundation for a better understanding of the
strengths and limitations of shape measurement methods.

[45] oai:arXiv.org:1302.0183 [pdf] - 1159437

IM3SHAPE: A maximum-likelihood galaxy shear measurement code for cosmic
gravitational lensing

Comments: 16 pages including appendices. Figure 8 is the money plot

Submitted: **2013-02-01**, last modified: **2013-02-04**

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 bulge-plus-disc 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 noise-free
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 sub-two percent level multiplicative bias, which is suitable for
analysis of current surveys. We fall short of the sub-percent level requirement
for upcoming surveys, which we attribute to a combination of noise bias and the
mis-match 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

[46] oai:arXiv.org:1203.5049 [pdf] - 1117484

Measurement and Calibration of Noise Bias in Weak Lensing Galaxy Shape
Estimation

Comments: 12 pages, 4 figures, submitted to MNRAS

Submitted: **2012-03-22**

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
Maximum-Likelihood forward model-fitting approach using a sum of co-elliptical
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 1-10%
on measured shears, which is significant for current and future lensing
surveys. We evaluate a simulation-based 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.

[47] oai:arXiv.org:1203.5050 [pdf] - 1117485

Noise bias in weak lensing shape measurements

Comments: 8 pages, 3 figures, submitted to MNRAS

Submitted: **2012-03-22**

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 6-parameter elliptical Gaussian model. We find that the
noise bias is generically of the order of the inverse-squared signal-to-noise
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
all-sky weak lensing surveys. We discuss possible ways to circumvent this
effect, including a calibration method using simulations discussed in an
associated paper.