
Cosmological parameter estimation is entering a new era. Large collaborations
need to coordinate highstakes analyses using multiple methods; furthermore
such analyses have grown in complexity due to sophisticated models of cosmology
and systematic uncertainties. In this paper we argue that modularity is the key
to addressing these challenges: calculations should be broken up into
interchangeable modular units with inputs and outputs clearly defined. We
present a new framework for cosmological parameter estimation, CosmoSIS,
designed to connect together, share, and advance development of inference tools
across the community. We describe the modules already available in CosmoSIS,
including CAMB, Planck, cosmic shear calculations, and a suite of samplers. We
illustrate it using demonstration code that you can run outofthebox with the
installer available at http://bitbucket.org/joezuntz/cosmosis

We construct forecasts for cosmological parameter constraints from weak
gravitational lensing surveys involving the Square Kilometre Array (SKA).
Considering matter content, dark energy and modified gravity parameters, we
show that the first phase of the SKA (SKA1) can be competitive with other Stage
III experiments such as the Dark Energy Survey (DES) and that the full SKA
(SKA2) can potentially form tighter constraints than Stage IV optical weak
lensing experiments, such as those that will be conducted with LSST,
WFIRSTAFTA or Euclidlike facilities. Using weak lensing alone, going from
SKA1 to SKA2 represents improvements by factors of $\sim10$ in matter, $\sim10$
in dark energy and $\sim5$ in modified gravity parameters. We also show, for
the first time, the powerful result that comparably tight constraints (within
$\sim5\%$) for both Stage III and Stage IV experiments, can be gained from
crosscorrelating shear maps between the optical and radio wavebands, a process
which can also eliminate a number of potential sources of systematic errors
which can otherwise limit the utility of weak lensing cosmology.

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.

Approximately half of the Universe's baryons are in a form that has been hard
to detect directly. However, the missing component can be traced through the
crosscorrelation of the thermal SunyaevZeldovich (tSZ) effect with weak
gravitational lensing. We build a model for this correlation and use it to
constrain the extended baryon component, employing data from the Canada France
Hawaii Lensing Survey and the {\it Planck\/} satellite. The measured
correlation function is consistent with an isothermal $\beta$model for the
halo gas pressure profile, and the 1 and 2halo terms are both detected at the
4$\sigma$ level. In addition, we measure the hydrostatic mass bias
$(1b)=0.79^{+0.07}_{0.10}$, which is consistent with numerical simulation
results and the constraints from Xray observations. The effective temperature
of the gas is found to be in the range ($7\times10^{5}$$3 \times10^{8}$)\,K,
with approximately $50\%$ of the baryons appearing to lie beyond the virial
radius of the halos, consistent with current expectations for the warmhot
intergalactic medium.

We examine the level of agreement between low redshift weak lensing data and
the CMB using measurements from the CFHTLenS and Planck+WMAP polarization. We
perform an independent analysis of the CFHTLenS six bin tomography results of
Heymans et al. (2013). We extend their systematics treatment and find the
cosmological constraints to be relatively robust to the choice of nonlinear
modeling, extension to the intrinsic alignment model and inclusion of baryons.
We find that the 90% confidence contours of CFHTLenS and Planck+WP do not
overlap even in the full 6dimensional parameter space of $\Lambda$CDM, so the
two datasets are discrepant. Allowing a massive active neutrino or tensor modes
does not significantly resolve the disagreement in the full ndimensional
parameter space. Our results differ from some in the literature because we use
the full tomographic information in the weak lensing data and marginalize over
systematics. We note that adding a sterile neutrino to $\Lambda$CDM does bring
the 8dimensional 64% contours to overlap, mainly due to the extra effective
number of neutrino species, which we find to be 0.84 $\pm$ 0.35 (68%) greater
than standard on combining the datasets. We discuss why this is not a
completely satisfactory resolution, leaving open the possibility of other new
physics or observational systematics as contributing factors. We provide
updated cosmology fitting functions for the CFHTLenS constraints and discuss
the differences from ones used in the literature.

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.

Accurate galaxy stellar masses are crucial to better understand the physical
mechanisms driving the galaxy formation process. We use synthetic star
formation and metal enrichment histories predicted by the {\sc galform} galaxy
formation model to investigate the precision with which various colours
$(m_{a}m_{b})$ can alone be used as diagnostics of the stellar masstolight
ratio. As an example, we find that, at $z=0$, the {\em intrinsic}
(B$_{f435w}$V$_{f606w}$) colour can be used to determine the intrinsic
restframe $V$band stellar masstolight ratio
($\log_{10}\Gamma_{V}=\log_{10}[(M/M_{\odot})/(L_{V}/L_{V\odot})]$) with a
precision of $\sigma_{lg\Gamma}\simeq 0.06$ when the initial mass function and
redshift are known beforehand. While the presence of dust, assuming a universal
attenuation curve, can have a systematic effect on the inferred masstolight
ratio using a singlecolour relation, this is typically small as it is often
possible to choose a colour for which the dust reddening vector is
approximately aligned with the $(m_{a}m_{b})\log_{10}\Gamma_{V}$ relation.
The precision with which the stellar masstolight ratio can be recovered using
a single colour diagnostic rivals implementations of SED fitting using more
information but in which simple parameterisations of the star formation and
metal enrichment histories are assumed. To facilitate the wide use of these
relations, we provide the optimal observer frame colour to estimate the stellar
masstolight ratio, along with the associated parameters, as a function of
redshift ($0<z<1.5$) for two sets of commonly used filters.

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

There are a number of approaches to testing General Relativity (GR) on linear
scales using parameterized frameworks for modifying cosmological perturbation
theory. It is sometimes assumed that the details of any given parameterization
are unimportant if one uses it as a diagnostic for deviations from GR. In this
brief report we argue that this is not necessarily so. First we show that
adopting alternative combinations of modifications to the field equations
significantly changes the constraints that one obtains. In addition, we show
that using a parameterization with insufficient freedom significantly tightens
the apparent theoretical constraints. Fundamentally we argue that it is almost
never appropriate to consider modifications to the perturbed Einstein equations
as being constraints on the effective gravitational constant, for example, in
the same sense that solar system constraints are. The only consistent
modifications are either those that grant neartotal freedom, as in
decomposition methods, or ones which map directly to a particular part of
theory space.

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.

There is a distinct possibility that current and future cosmological data can
be used to constrain Einstein's theory of gravity on the very largest scales.
To be able to do this in a modelindependent way, it makes sense to work with a
general parameterization of modified gravity. Such an approach would be
analogous to the Parameterized PostNewtonian (PPN) approach which is used on
the scale of the Solar System. A few such parameterizations have been proposed
and preliminary constraints have been obtained. We show that the majority of
such parameterizations are only exactly applicable in the quasistatic regime.
On larger scales they fail to encapsulate the full behaviour of typical models
currently under consideration. We suggest that it may be possible to capture
the additions to the `Parameterized PostFriedmann' (PPF) formalism by treating
them akin to fluid perturbations.

It is well known that observations of the cosmic microwave background (CMB)
are highly sensitive to the spatial curvature of the Universe, k. Here we find
that what is in fact being tightly constrained by small angle fluctuations is
spatial curvature near the surface of last scattering, and that if we allow k
to be a function of position, rather than taking a constant value everywhere,
then considerable spatial curvature is permissible within our own locale. This
result is of interest for the giant void models that attempt to explain the
supernovae observations without Dark Energy. We find voids models with a
homogeneous big bang can be compatible with the observed small angle CMB, but
only if they exist in a positively curved universe. To be compatible with local
measurements of H_0, however, we find that a radially varying bang time is
required.

Statistically anomalous signals in the microwave background have been
extensively studied in general in multipole space, and in real space mainly for
circular and other simple patterns. In this paper we search for a range of
nontrivial patterns in the temperature data from WMAP 7year observations. We
find a very significant detection of a number of such features and discuss
their consequences for the essential character of the cosmos.

We analyze correlations between the first letter of the name of an author and
the number of citations their papers receive. We look at simple mean counts,
numbers of highlycited papers, and normalized hindices, by letter. To our
surprise, we conclude that orthographically senior authors produce a better
body of work than their colleagues, despite some evidence of discrimination
against them.

We describe and discuss the application of Gumbel statistics, which model
extreme events, to WMAP 5year measurements of the cosmic microwave background.
We find that temperature extrema of the CMB are well modelled by the Gumbel
formalism and describe tests for Gaussianity that the approach can provide.
Comparison to simulations reveals Gumbel statistics to have only weak
discriminatory power for the conventional statistic: $f_{NL}<1000$, though it
may probe other regimes of nonGaussianity. Tests based on hemispheric cuts
reveal interesting alignment with other reported CMB anomalies. The approach
has the advantage of model independence and may find further utility with
smaller scale data.

We consider the effects of large structures in the Universe on the Hubble
diagram. This problem is treated nonlinearly by considering a Swiss Cheese
model of the Universe in which underdense voids are represented as negatively
curved regions of spacetime. Exact expressions for luminosity distances and
redshifts are used to investigate the nonlinear effects of structure on the
magnitudes of astrophysical sources. It is found that the intervening voids we
consider, between the observer and source, produce changes in apparent
magnitude of less than 0.012. Sources inside voids, however, can be affected
considerably at redshifts below z~0.5. By averaging observable quantities over
many randomly generated distributions of voids we find that the presence of
these structures has the effect of introducing a dispersion around the mean,
which itself can be displaced the background value. Observers in an
inhomogeneous universe, who take averages of observables along many different
lines of sight, may then introduce systematic biases, and underestimate
errors, if these effects are not taken into account. Estimates of the potential
size of these effects are made using data from simulated largescale structure.

Observations by the Wilkinson Microwave Anisotropy Probe (WMAP) satellite
have identified an excess of microwave emission from the centre of the Milky
Way. It has been suggested that this WMAP haze emission could potentially be
synchrotron emission from relativistic electrons and positrons produced in the
annihilations of one (or more) species of dark matter particles. In this paper
we recalculate the intensity and morphology of the WMAP haze using a
multilinear regression involving fullsky templates of the dominant forms of
galactic foreground emission, using two different CMB sky signal estimators.
The first estimator is a posterior mean CMB map, marginalized over a general
foreground model using a Gibbs sampling technique, and the other is the ILC map
produced by the WMAP team. Earlier analyses of the WMAP haze used the ILC map,
which is more contaminated by galactic foregrounds than the Gibbs map. In
either case, we reconfirm earlier results that a statistically significant
residual emission remains after foreground subtraction that is concentrated
around the galactic centre. However, we find that the significance of this
emission can be significantly reduced by allowing for a subtle spatial
variation in the frequency dependence of soft synchrotron emission in the inner
and outer parts of the galaxy. We also reinvestigate the prospect of a
neutralino dark matter interpretation of the origin of the haze, and find that
significant boosting in the dark matter annihilation rate is required, relative
to that obtained with a smooth galactic dark matter distribution, in order to
reproduce the inferred residual emission, contrary to that deduced in several
recent studies.

Residual star formation at late times in earlytype galaxies and their
progenitors must be suppressed in order to explain the population of red,
passively evolving systems we see today. Likewise, residual or newly accreted
reservoirs of molecular gas that are fuelling star formation must be destroyed.
This suppression of star formation in earlytype galaxies is now commonly
attributed to AGN feedback wherein the reservoir of gas is heated and expelled
during a phase of accretion onto the central supermassive black hole. However,
direct observational evidence for a link between the destruction of this
molecular gas and an AGN phase has been missing so far. We present new
mmwavelength observations from the IRAM 30m telescope of a sample of low
redshift SDSS earlytype galaxies currently undergoing this process of
quenching of latetime star formation. Our observations show that the
disappearance of the molecular gas coincides within less than 100 Myr with the
onset of accretion onto the black hole and is too rapid to be due to star
formation alone. Since our sample galaxies are not associated to powerful
quasar activity or radio jets, we conclude that lowluminosity AGN episodes are
sufficient to suppress residual star formation in earlytype galaxies. This
`suppression mode' of AGN feedback is very different from the `truncation mode'
linked to powerful quasar activity during early phases of galaxy formation.

We discuss Cosmic Microwave Background constraints on the causal set theory
of quantum gravity, which has made testable predictions about the nature of
dark energy. We flesh out previously discussed heuristic constraints by showing
how the power spectrum of causal set dark energy fluctuations can be found from
the overlap volumes of past light cones of points in the universe. Using a
modified Boltzmann code we put constraints on the single parameter of the
theory that are somewhat stronger than previous ones. We conclude that causal
set theory cannot explain latetime acceleration without radical alterations to
General Relativity.