
The Planck full mission cosmic microwave background(CMB) temperature and
Emode polarization maps are analysed to obtain constraints on primordial
nonGaussianity(NG). Using three classes of optimal bispectrum estimators 
separable templatefitting (KSW), binned, and modal  we obtain consistent
values for the local, equilateral, and orthogonal bispectrum amplitudes,
quoting as our final result from temperature alone fNL^local=2.5+\5.7,
fNL^equil=16+\70 and fNL^ortho=34+\33(68%CL). Combining temperature and
polarization data we obtain fNL^local=0.8+\5.0, fNL^equil=4+\43 and
fNL^ortho=26+\21 (68%CL). The results are based on crossvalidation of these
estimators on simulations, are stable across component separation techniques,
pass an extensive suite of tests, and are consistent with Minkowski functionals
based measurements. The effect of timedomain deglitching systematics on the
bispectrum is negligible. In spite of these test outcomes we conservatively
label the results including polarization data as preliminary, owing to a known
mismatch of the noise model in simulations and the data. Beyond fNL estimates,
we present modelindependent reconstructions of the CMB bispectrum and derive
constraints on early universe scenarios that generate NG, including general
singlefield and axion inflation, initial state modifications, parityviolating
tensor bispectra, and directionally dependent vector models. We also present a
wide survey of scaledependent oscillatory bispectra, and we look for
isocurvature NG. Our constraint on the local primordial trispectrum amplitude
is gNL^local=(9.0+\7.7)x10^4 (68%CL), and we perform an analysis of
additional trispectrum shapes. The global picture is one of consistency with
the premises of the LambdaCDM cosmology, namely that the structure we observe
today was sourced by adiabatic, passive, Gaussian, and primordial seed
perturbations.[abridged]

We study the implications of Planck data for models of dark energy (DE) and
modified gravity (MG), beyond the cosmological constant scenario. We start with
cases where the DE only directly affects the background evolution, considering
Taylor expansions of the equation of state, principal component analysis and
parameterizations related to the potential of a minimally coupled DE scalar
field. When estimating the density of DE at early times, we significantly
improve present constraints. We then move to general parameterizations of the
DE or MG perturbations that encompass both effective field theories and the
phenomenology of gravitational potentials in MG models. Lastly, we test a range
of specific models, such as kessence, f(R) theories and coupled DE. In
addition to the latest Planck data, for our main analyses we use baryonic
acoustic oscillations, typeIa supernovae and local measurements of the Hubble
constant. We further show the impact of measurements of the cosmological
perturbations, such as redshiftspace distortions and weak gravitational
lensing. These additional probes are important tools for testing MG models and
for breaking degeneracies that are still present in the combination of Planck
and background data sets. All results that include only background
parameterizations are in agreement with LCDM. When testing models that also
change perturbations (even when the background is fixed to LCDM), some tensions
appear in a few scenarios: the maximum one found is \sim 2 sigma for Planck
TT+lowP when parameterizing observables related to the gravitational potentials
with a chosen time dependence; the tension increases to at most 3 sigma when
external data sets are included. It however disappears when including CMB
lensing.

Lensing of the Cosmic Microwave Background (CMB) changes the morphology of
pattern of temperature fluctuations, so topological descriptors such as
Minkowski Functionals can probe the gravity model responsible for the lensing.
We show how the recently introduced twototwo and threetoone kurtspectra
(and their associated correlation functions), which depend on the power
spectrum of the lensing potential, can be used to probe modified gravity
theories such as $f({R})$ theories of gravity and quintessence models. We also
investigate models based on effective field theory, which include the
constant$\Omega$ model, and lowenergy Ho\vrava theories. Estimates of the
cumulative signaltonoise for detection of lensinginduced morphology changes,
reaches ${\cal O}(10^3)$ for the future planned CMB polarization mission
COrE$^{+}$. Assuming foreground removal is possible to $\ell_{max}=3000$, we
show that many modified gravity theories can be rejected with a high level of
significance, making this technique comparable in power to galaxy weak lensing
or redshift surveys. These topological estimators are also useful in
distinguishing {\em lensing} from other scattering secondaries at the level of
the fourpoint function or trispectrum. Examples include the kinetic
SunyaevZel'dovich (kSZ) effect which shares, with lensing, a lack of spectral
distortion. We also discuss the complication of foreground contamination from
unsubtracted point sources.

The European Space Agency's Planck satellite, dedicated to studying the early
Universe and its subsequent evolution, was launched 14~May 2009 and scanned the
microwave and submillimetre sky continuously between 12~August 2009 and
23~October 2013. In February~2015, ESA and the Planck Collaboration released
the second set of cosmology products based on data from the entire Planck
mission, including both temperature and polarization, along with a set of
scientific and technical papers and a webbased explanatory supplement. This
paper gives an overview of the main characteristics of the data and the data
products in the release, as well as the associated cosmological and
astrophysical science results and papers. The science products include maps of
the cosmic microwave background (CMB), the thermal SunyaevZeldovich effect,
and diffuse foregrounds in temperature and polarization, catalogues of compact
Galactic and extragalactic sources (including separate catalogues of
SunyaevZeldovich clusters and Galactic cold clumps), and extensive simulations
of signals and noise used in assessing the performance of the analysis methods
and assessment of uncertainties. The likelihood code used to assess
cosmological models against the Planck data are described, as well as a CMB
lensing likelihood. Scientific results include cosmological parameters deriving
from CMB power spectra, gravitational lensing, and cluster counts, as well as
constraints on inflation, nonGaussianity, primordial magnetic fields, dark
energy, and modified gravity.

The ESA's Planck satellite, dedicated to studying the early Universe and its
subsequent evolution, was launched 14 May 2009 and has been scanning the
microwave and submillimetre sky continuously since 12 August 2009. This paper
gives an overview of the mission and its performance, the processing, analysis,
and characteristics of the data, the scientific results, and the science data
products and papers in the release. The science products include maps of the
CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and
extragalactic sources, and a list of sources detected through the SZ effect.
The likelihood code used to assess cosmological models against the Planck data
and a lensing likelihood are described. Scientific results include robust
support for the standard sixparameter LCDM model of cosmology and improved
measurements of its parameters, including a highly significant deviation from
scale invariance of the primordial power spectrum. The Planck values for these
parameters and others derived from them are significantly different from those
previously determined. Several largescale anomalies in the temperature
distribution of the CMB, first detected by WMAP, are confirmed with higher
confidence. Planck sets new limits on the number and mass of neutrinos, and has
measured gravitational lensing of CMB anisotropies at greater than 25 sigma.
Planck finds no evidence for nonGaussianity in the CMB. Planck's results agree
well with results from the measurements of baryon acoustic oscillations. Planck
finds a lower Hubble constant than found in some more local measures. Some
tension is also present between the amplitude of matter fluctuations derived
from CMB data and that derived from SZ data. The Planck and WMAP power spectra
are offset from each other by an average level of about 2% around the first
acoustic peak.

We show how crosscorrelating a high redshift external tracer field, such as
the 21cm neutral hydrogen distribution and product maps involving Cosmic
Microwave Background (CMB) temperature and polarisation fields, that probe
mixed bispectrum involving these fields, can help to determine the reionization
history of the Universe, beyond what can be achieved from crossspectrum
analysis. Taking clues from recent studies for the detection of primordial
nonGaussianity (Munshi & Heavens 2010), we develop a set of estimators that
can study reionization using a power spectrum associated with the bispectrum
(or skewspectrum). We use the matched filtering inherent in this method to
investigate different reionization histories. We check to what extent they can
be used to rule out various models of reionization and study cross
contamination from different sources such as the lensing of the CMB. The
estimators can be finetuned to optimize study of a specific reionization
history. We consider three different types of tracers in our study, namely:
protogalaxies; 21cm maps of neutral hydrogen; quasars. We also consider four
alternative models of reionization. We find that the cumulative signaltonoise
(S/N) for detection at $\ell_{max}=2000$ can reach $\mathcal{O}(70)$ for cosmic
variance limited allsky experiments. Combining 100GHz, 143GHz and 217GHz
channels of the Planck experiment, we find that the (S/N) lies in the range
$\mathcal{O}(5)\mathcal{O}(35)$. The (S/N) depends on the specific choice of a
tracer field, and multiple tracers can be effectively used to map out the
entire reionization history with reasonable S/N. Contamination from weak
lensing is investigated and found to be negligible, and the effects of Thomson
scattering from patchy reionizationare also considered.

We use the optimised skewspectrum as well as the skewspectra associated
with the Minkowski Functionals (MFs) to test the possibility of using the
crosscorrelation of the Integrated SachsWolfe effect (ISW) and lensing of the
cosmic microwave background (CMB) radiation to detect deviations in the theory
of gravity away from General Relativity (GR). We find that the although both
statistics can put constraints on modified gravity, the optimised skewspectra
are especially sensitive to the parameter $\rm B_0$ that denotes the the {\em
Compton wavelength} of the scalaron at the present epoch. We investigate three
modified gravity theories, namely: the PostParametrised Friedmanian (PPF)
formalism; the HuSawicki (HS) model; and the BertschingerZukin (BZ)
formalism. Employing a likelihood analysis for an experimental setup similar to
ESA's Planck mission, we find that, assuming GR to be the correct model, we
expect the constraints from the first two skewspectra, $S_{\ell}^{(0)}$ and
$S_{\ell}^{(1)}$, to be the same: $\rm B_0<0.45$ at $95%$ confidence level
(CL), and $\rm B_0<0.67$ at $99%$ CL in the BZ model. The third skewspectrum
does not give any meaningful constraint. We find that the optimal skewspectrum
provides much more powerful constraint, giving $\rm B_0<0.071$ at $95%$ CL and
$\rm B_0<0.15$ at $99%$ CL, which is essentially identical to what can be
achieved using the full bispectrum.

The Planck nominal mission cosmic microwave background (CMB) maps yield
unprecedented constraints on primordial nonGaussianity (NG). Using three
optimal bispectrum estimators, separable templatefitting (KSW), binned, and
modal, we obtain consistent values for the primordial local, equilateral, and
orthogonal bispectrum amplitudes, quoting as our final result fNL^local=
2.7+/5.8, fNL^equil= 42+/75, and fNL^ortho= 25+\39 (68% CL statistical).
NG is detected in the data; using skewC_l statistics we find a nonzero
bispectrum from residual point sources, and the ISWlensing bispectrum at a
level expected in the LambdaCDM scenario. The results are based on
comprehensive crossvalidation of these estimators on Gaussian and nonGaussian
simulations, are stable across component separation techniques, pass an
extensive suite of tests, and are confirmed by skewC_l, wavelet bispectrum and
Minkowski functional estimators. Beyond estimates of individual shape
amplitudes, we present modelindependent, 3dimensional reconstructions of the
Planck CMB bispectrum and thus derive constraints on earlyUniverse scenarios
that generate primordial NG, including general singlefield models of
inflation, excited initial states (nonBunchDavies vacua), and
directionallydependent vector models. We provide an initial survey of
scaledependent feature and resonance models. These results bound both general
singlefield and multifield model parameter ranges, such as the speed of
sound, c_s \geq 0.02 (95% CL), in an effective field theory parametrization,
and the curvaton decay fraction r_D \geq 0.15 (95% CL). The Planck data
significantly limit the viable parameter space of the ekpyrotic/cyclic
scenarios. The amplitude of the 4point function in the local model tauNL <
2800 (95% CL). These constraints represent the highest precision tests to date
of physical mechanisms for the origin of cosmic structure.

Euclid is a spacebased survey mission from the European Space Agency
designed to understand the origin of the Universe's accelerating expansion. It
will use cosmological probes to investigate the nature of dark energy, dark
matter and gravity by tracking their observational signatures on the geometry
of the universe and on the cosmic history of structure formation. The mission
is optimised for two independent primary cosmological probes: Weak
gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid
payload consists of a 1.2 m Korsch telescope designed to provide a large field
of view. It carries two instruments with a common fieldofview of ~0.54 deg2:
the visual imager (VIS) and the near infrared instrument (NISP) which contains
a slitless spectrometer and a three bands photometer. The Euclid wide survey
will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20
deg2 deep fields. For WL, Euclid measures the shapes of 3040 resolved galaxies
per arcmin2 in one broad visible R+I+Z band (550920 nm). The photometric
redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are
derived from three additional Euclid NIR bands (Y, J, H in the range 0.922.0
micron), complemented by ground based photometry in visible bands derived from
public data or through engaged collaborations. The BAO are determined from a
spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless
spectrometer, with spectral resolution ~250, predominantly detects Ha emission
line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision
20152025 programme, with a foreseen launch date in 2019. This report (also
known as the Euclid Red Book) describes the outcome of the Phase A study.

Weak gravitational lensing is responsible for the shearing and magnification
of the images of highredshift sources due to the presence of intervening
matter. The distortions are due to fluctuations in the gravitational potential,
and are directly related to the distribution of matter and to the geometry and
dynamics of the Universe. As a consequence, weak gravitational lensing offers
unique possibilities for probing the Dark Matter and Dark Energy in the
Universe. In this review, we summarise the theoretical and observational state
of the subject, focussing on the statistical aspects of weak lensing, and
consider the prospects for weak lensing surveys in the future. Weak
gravitational lensing surveys are complementary to both galaxy surveys and
cosmic microwave background (CMB) observations as they probe the unbiased
nonlinear matter power spectrum at modest redshifts. Most of the cosmological
parameters are accurately estimated from CMB and largescale galaxy surveys, so
the focus of attention is shifting to understanding the nature of Dark Matter
and Dark Energy. On the theoretical side, recent advances in the use of 3D
information of the sources from photometric redshifts promise greater
statistical power, and these are further enhanced by the use of statistics
beyond twopoint quantities such as the power spectrum. The use of 3D
information also alleviates difficulties arising from physical effects such as
the intrinsic alignment of galaxies, which can mimic weak lensing to some
extent. (Abridged)

We use the LUQAS sample of 27 high resolution, high signaltonoise QSO
absorption spectra (Kim et al. 2004) and the results from Croft et al. (2002)
together with a suite of highresolution hydrodynamical simulations run with
the GADGETII code, to infer the linear dark matter power spectrum on scales of
\~ 0.330 Mpc at z=2.125 and at z=2.72.

We present a determination of the bispectrum of the flux in the Lymanalpha
forest of QSO absorption spectra obtained from LUQAS which consists of spectra
observed with the high resolution Echelle spectrograph UVES. Typical errors on
the observed bispectrum as obtained from a jackknife estimator are ~ 50%. For
wavenumbers in the range 0.03 s/km <k < 0.1 s/km the observed bispectrum agrees
within the errors with that of the synthetic absorption spectra obtained from
numerical hydrosimulations of a LCDM model with and without feedback from star
formation. Including galactic feedback changes the bispectrum by less than 10%.
At smaller wavenumbers the associated metal absorption lines contribute about
50% to the bispectrum and the observed bispectrum exceeds that of the
simulations. At wavenumbers k< 0.03 (s/km) secondorder perturbation theory
applied to the flux spectrum gives a reasonable (errors smaller than 30%)
approximation to the bispectra of observed and simulated absorption spectra.
The bispectrum of the observed absorption spectra also agrees, within the
errors, with that of a randomized set of absorption spectra where a random
shift in wavelength has been added to absorption lines identified with VPFIT.
This suggests that for a sample of the size presented here, the errors on the
bispectrum are too large to discriminate between models with very different 3D
distribution of Lymanalpha absorption. If it were possible to substantially
reduce these errors for larger samples of absorption spectra, the bispectrum
might become an important statistical tool for probing the growth of
gravitational structure in the Universe at redshift z>2.

We describe different methods for estimating the bispectrum of Cosmic
Microwave Background data. In particular we construct a minimum variance
estimator for the flatsky limit and compare results with previouslystudied
frequentist methods. Application to the MAXIMA dataset shows consistency with
primordial Gaussianity. Weak quadratic nonGaussianity is characterised by a
tunable parameter $f_{NL}$, corresponding to nonGaussianity at a level $\sim
10^{5}f_{NL}$ (ratio of nonGaussian to Gaussian terms), and we find limits of
$f_{NL}<950$ for the minimumvariance estimator and $f_{NL}<1650$ for the
usual frequentist estimator. These are the tightest limits on primordial
nonGaussianity which include the full effects of the radiation transfer
function.

We investigate the influence of the star formation rate on the TullyFisher
relation. We find that a simple model which combines the empiricallydetermined
starformation rate with the expected properties of galaxy halos provides a
remarkably good fit to the absolute magnituderotation speed correlation. We
find that the powerlaw nature, its slope, normalisation and scatter, are all
readily accounted for if the Universe has a low density parameter, with or
without a cosmological constant and disks are assembled at $z \sim 1  1.5$.
Moreover, this agreement is found simultaneously in 4 wavebands. An Einsteinde
Sitter Universe produces disks which are too faint unless the disks are
assembled at $z\sim 0.5$. The scatter in the relation is due to a combination
of the expected range of spin parameters of the halos and the range of
formation redshifts. The source of the scatter opens up possibilities of a
better galaxy distance indicator, if spectroscopic observations of globular
clusters can be used to determine the halo rotation.