
The number of observed dwarf galaxies, with dark matter mass $\lesssim
10^{11}$ M$_{\odot}$ in the Milky Way or the Andromeda galaxy does not agree
with predictions from the successful $\Lambda$CDM paradigm. To alleviate this
problem a suppression of dark matter clustering power on very small scales has
been conjectured. However, the abundance of dark matter halos outside our
immediate neighbourhood (the Local Group) seem to agree with the
$\Lambda$CDMexpected abundance. Here we connect these problems to
observations of weak lensing cosmic shear, pointing out that cosmic shear can
make significant statements about the missing satellites problem in a
statistical way. As an example and pedagogical application we use recent
constraints on smallscales power suppression from measurements of the CFHTLenS
data. We find that, on average, in a region of $\sim $Gpc$^3$ there is no
significant smallscale power suppression. This implies that suppression of
smallscale power is not a viable solution to the `missing satellites problem'
or, alternatively, that on average in this volume there is no `missing
satellites problem' for dark matter masses $\gtrsim 5 \times 10^9$ M$_{\odot}$.
Further analysis of current and future weak lensing surveys will probe much
smaller scales, $k > 10h$ Mpc$^{1}$ corresponding roughly to masses $M < 10^9
M_{\odot}$.

We compute in detail how deviations from Einstein gravity at the inflation
energy scale could appear as nonGaussian features in the sky. To illustrate
this we use multifield $\alpha$attractor models in the framework of
supergravity to realise inflation. We find no obvious obstacle for having
choices of model's parameters that generate nonGaussian features of the
equilateral and local type at the $\cal O$(1) level in the $f_{\rm NL}$
nonGaussianity parameter, thus being potentially detectable in future
cosmological surveys. This nonGaussianity has its origin in either the
noncanonical kinetic term, the interactions of the fields in the potential or
the $\alpha$parameter, or a combination of these three. This opens up the
exciting possibility of constraining the law of gravity at energy scales close
to the Planck one.

It is broadly accepted that Supermassive Black Holes (SMBHs) are located in
the centers of most massive galaxies, although there is still no convincing
scenario for the origin of their massive seeds. It has been suggested that
primordial black holes (PBHs) of masses $\gtrsim 10^{2} M_\odot$ may provide
such seeds, which would grow to become SMBHs. We suggest an observational test
to constrain this hypothesis: gas accretion around PBHs during the cosmic dark
ages powers the emission of high energy photons which would modify the spin
temperature as measured by 21cm Intensity Mapping (IM) observations. We model
and compute their contribution to the standard skyaveraged signal and power
spectrum of 21cm IM, accounting for its substructure and angular dependence for
the first time. If PBHs exist, the skyaveraged 21cm IM signal in absorption
would be higher, while we expect an increase in the power spectrum for
$\ell~\gtrsim 10^210^3$. We also forecast PBH detectability and measurement
errors in the abundance and Eddington ratios for different fiducial parameter
configurations for various future experiments, ranging from SKA to a futuristic
radio array on the dark side of the Moon. While the SKA could provide a
detection, only a more ambitious experiment would provide accurate
measurements.

We forecast constraints on primordial nonGaussianity (PNG) and bias
parameters from measurements of galaxy power spectrum and bispectrum in future
radio continuum and optical surveys. In the galaxy bispectrum, we consider a
comprehensive list of effects, including the bias expansion for nonGaussian
initial conditions up to second order, redshift space distortions, redshift
uncertainties and theoretical errors. These effects are all combined in a
single PNG forecast for the first time. Moreover, we improve the bispectrum
modelling over previous forecasts, by accounting for trispectrum contributions.
All effects have an impact on final predicted bounds, which varies with the
type of survey. We find that the bispectrum can lead to improvements up to a
factor $\sim 5$ over bounds based on the power spectrum alone, leading to
significantly better constraints for localtype PNG, with respect to current
limits from \textit{Planck}. Future radio and photometric surveys could obtain
a measurement error of $\sigma(f_{\mathrm{NL}}^{\mathrm{loc}}) \approx 0.2$. In
the case of equilateral PNG, galaxy bispectrum can improve upon present bounds
only if significant improvements in the redshift determinations of future,
large volume, photometric or radio surveys could be achieved. For orthogonal
nonGaussianity, expected constraints are generally comparable to current ones.

The expansion history of the Universe can be constrained in a
cosmologyindependent way by measuring the differential age evolution of cosmic
chronometers. This yields a measurement of the Hubble parameter as a function
of redshift. The most reliable cosmic chronometers known so far are extremely
massive and passively evolving galaxies. Agedating these galaxies is, however,
a difficult task, and even a small contribution of an underlying young stellar
population ("frosting") could, in principle, affect the age estimate and its
cosmological interpretation. We present several spectral indicators to detect,
quantify and constrain such a young component in old galaxies, and study how
their combination can be used to maximize the purity of cosmic chronometers
selection. In particular, we analyze the CaII H/K ratio, and the presence (or
absence) of H$\alpha$ and [OII] emission lines, higher order Balmer absorption
lines, and UV flux; each indicator is especially sensitive to a particular age
range, allowing us to detect young components ranging between 10 Myr and 1 Gyr.
The combination of these indicators minimizes the contamination below the 1%
level, and offers a way to control the systematic error it could introduce in
the Hubble parameter determination. We show that for our previous measurements
this effect is well below the current errors. We envision that these indicators
will be instrumental in strengthening the selection criterion of cosmic
chronometers, paving the road for a robust and reliable dating of the old
population and its cosmological interpretation.

Given the important role that the galaxy bispectrum has recently acquired in
cosmology and the scale and precision of forthcoming galaxy clustering
observations, it is timely to derive the full expression of the largescale
bispectrum going beyond approximated treatments which neglect integrated terms
or higherorder bias terms or use the Limber approximation. On cosmological
scales, relativistic effects that arise from observing on the past lightcone
alter the observed galaxy number counts, therefore leaving their imprints on
Npoint correlators at all orders. In this paper we compute for the first time
the bispectrum including all general relativistic, local and integrated,
effects at second order, the tracers' bias at second order, geometric effects
as well as the primordial nonGaussianity contribution. This is timely
considering that future surveys will probe scales comparable to the horizon
where approximations widely used currently may not hold; neglecting these
effects may introduce biases in estimation of cosmological parameters as well
as primordial nonGaussianity.

The model in which Primordial Black Holes (PBHs) constitute a nonnegligible
fraction of the dark matter has (re)gained popularity after the first
detections of binary black hole mergers. Most of the observational constraints
to date have been derived assuming a single mass for all the PBHs, although
some more recent works tried to generalize constraints to the case of extended
mass functions. Here we derive a general methodology to obtain constraints for
any PBH Extended Mass Distribution (EMD) and any observables in the desired
mass range. Starting from those obtained for a monochromatic distribution, we
convert them into constraints for EMDs by using an equivalent, effective mass
$M_{\rm eq}$ that depends on the specific observable. We highlight how limits
of validity of the PBH modelling affect the EMD parameter space. Finally, we
present converted constraints on the total abundance of PBH from microlensing,
stellar distribution in ultrafaint dwarf galaxies and CMB accretion for
Lognormal and Power Law mass distributions, finding that EMD constraints are
generally stronger than monochromatic ones.

The possibility that a relevant fraction of the dark matter might be
comprised of Primordial Black Holes (PBHs) has been seriously reconsidered
after LIGO's detection of a $\sim 30 M_{\odot}$ binary black holes merger.
Despite the strong interest in the model, there is a lack of studies on
possible cosmological implications and effects on cosmological parameters
inference. We investigate correlations with the other standard cosmological
parameters using cosmic microwave background observations, finding significant
degeneracies, especially with the tilt of the primordial power spectrum and the
sound horizon at radiation drag. However, these degeneracies can be greatly
reduced with the inclusion of small scale polarization data. We also explore if
PBHs as dark matter in simple extensions of the standard $\Lambda$CDM
cosmological model induces extra degeneracies, especially between the
additional parameters and the PBH's ones. Finally, we present cosmic microwave
background constraints on the fraction of dark matter in PBHs, not only for
monochromatic PBH mass distributions but also for popular extended mass
distributions. Our results show that extended mass distribution's constraints
are tighter, but also that a considerable amount of constraining power comes
from the high$\ell$ polarization data. Moreover, we constrain the shape of
such mass distributions in terms of the correspondent constraints on the PBH
mass fraction.

The configuration of the three neutrino masses can take two forms, known as
the normal and inverted hierarchies. We compute the Bayesian evidence
associated with these two hierarchies. Previous studies found a mild preference
for the normal hierarchy, and this was driven by the asymmetric manner in which
cosmological data has confined the available parameter space. Here we identify
the presence of a second asymmetry, which is imposed by data from neutrino
oscillations. By combining constraints on the squaredmass splittings with the
limit on the sum of neutrino masses of $\Sigma m_\nu < 0.13$ eV, and using a
minimally informative prior on the masses, we infer odds of 42:1 in favour of
the normal hierarchy, which is classified as "strong" in the Jeffreys' scale.
We explore how these odds may evolve in light of higher precision cosmological
data, and discuss the implications of this finding with regards to the nature
of neutrinos. Finally the individual masses are inferred to be $m_1 =
3.80^{+26.2}_{3.73} \, \text{meV}, m_2 = 8.8^{+18}_{1.2} \, \text{meV}, m_3 =
50.4^{+5.8}_{1.2} \, \text{meV}$ ($95\%$ credible intervals).

In this paper we discuss the commonlyused limiting cases, or approximations,
for twopoint cosmic shear statistics. We discuss the most prominent
assumptions in this statistic: the flatsky (small angle limit), the Limber
(Besseltodelta function limit) and the Hankel transform (large lmode limit)
approximations; that the vast majority of cosmic shear results to date have
used simultaneously. We find that the combined effect of these approximations
can suppress power by >1% on scales of l<40. A fully nonapproximated cosmic
shear study should use a sphericalsky, nonLimberapproximated power spectrum
analysis; and a transform involving Wigner smalld matrices in place of the
Hankel transform. These effects, unaccounted for, would constitute at least 11%
of the total budget for systematic effects for a power spectrum analysis of a
Euclidlike experiment; but they are unnecessary.

Obtaining highsensitivity measurements of degreescale cosmic microwave
background (CMB) polarization is the most direct path to detecting primordial
gravitational waves. Robustly recovering any primordial signal from the
dominant foreground emission will require highfidelity observations at
multiple frequencies, with excellent control of systematics. We explore the
potential for a new platform for CMB observations, the Airlander 10 hybrid air
vehicle, to perform this task. We show that the Airlander 10 platform,
operating at commercial airline altitudes, is wellsuited to mapping
frequencies above 220 GHz, which are critical for cleaning CMB maps of dust
emission. Optimizing the distribution of detectors across frequencies, we
forecast the ability of Airlander 10 to clean foregrounds of varying complexity
as a function of altitude, demonstrating its complementarity with both existing
(Planck) and ongoing (CBASS) foreground observations. This novel platform
could play a key role in defining our ultimate view of the polarized microwave
sky.

Cosmological observables show a dependence with the neutrino mass, which is
partially degenerate with parameters of extended models of gravity. We study
and explore this degeneracy in Horndeski generalized scalartensor theories of
gravity. Using forecasted cosmic microwave background and galaxy power spectrum
datasets, we find that a single parameter in the linear regime of the effective
theory dominates the correlation with the total neutrino mass. For any given
mass, a particular value of this parameter approximately cancels the power
suppression due to the neutrino mass at a given redshift. The extent of the
cancellation of this degeneracy depends on the cosmological largescale
structure data used at different redshifts. We constrain the parameters and
functions of the effective gravity theory and determine the influence of
gravity on the determination of the neutrino mass from present and future
surveys.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV groundbased dark
energy experiment that will study baryon acoustic oscillations (BAO) and the
growth of structure through redshiftspace distortions with a widearea galaxy
and quasar redshift survey. To trace the underlying dark matter distribution,
spectroscopic targets will be selected in four classes from imaging data. We
will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to
even higher redshift, DESI will target bright [O II] emission line galaxies up
to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying
dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the
Ly$\alpha$ forest absorption features in their spectra, which will be used to
trace the distribution of neutral hydrogen. When moonlight prevents efficient
observations of the faint targets of the baseline survey, DESI will conduct a
magnitudelimited Bright Galaxy Survey comprising approximately 10 million
galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy
and quasar redshifts will be obtained to measure the BAO feature and determine
the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV groundbased dark
energy experiment that will study baryon acoustic oscillations and the growth
of structure through redshiftspace distortions with a widearea galaxy and
quasar redshift survey. The DESI instrument is a roboticallyactuated,
fiberfed spectrograph capable of taking up to 5,000 simultaneous spectra over
a wavelength range from 360 nm to 980 nm. The fibers feed ten threearm
spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500,
depending on wavelength. The DESI instrument will be used to conduct a
fiveyear survey designed to cover 14,000 deg$^2$. This powerful instrument
will be installed at prime focus on the 4m Mayall telescope in Kitt Peak,
Arizona, along with a new optical corrector, which will provide a threedegree
diameter field of view. The DESI collaboration will also deliver a
spectroscopic pipeline and data management system to reduce and archive all
data for eventual public use.

Primordial features are one of the most important extensions of the Standard
Model of cosmology, providing a wealth of information on the primordial
universe, ranging from discrimination between inflation and alternative
scenarios, new particle detection, to fine structures in the inflationary
potential. We study the prospects of future largescale structure (LSS) surveys
on the detection and constraints of these features. We classify primordial
feature models into several classes, and for each class we present a simple
template of power spectrum that encodes the essential physics. We study how
well the most ambitious LSS surveys proposed to date, including both
spectroscopic and photometric surveys, will be able to improve the constraints
with respect to the current Planck data. We find that these LSS surveys will
significantly improve the experimental sensitivity on features signals that are
oscillatory in scales, due to the 3D information. For a broad range of models,
these surveys will be able to reduce the errors of the amplitudes of the
features by a factor of 5 or more, including several interesting candidates
identified in the recent Planck data. Therefore, LSS surveys offer an
impressive opportunity for primordial feature discovery in the next decade or
two. We also compare the advantages of both types of surveys.

We investigate our knowledge of early universe cosmology by exploring how
much additional energy density can be placed in different components beyond
those in the $\Lambda$CDM model. To do this we use a method to separate early
and lateuniverse information enclosed in observational data, thus markedly
reducing the modeldependency of the conclusions. We find that the 95\%
credibility regions for extra energy components of the early universe at
recombination are: nonaccelerating additional fluid density parameter
$\Omega_{\rm MR} < 0.006$ and extra radiation parameterised as extra effective
neutrino species $2.3 < N_{\rm eff} < 3.2$ when imposing flatness. Our
constraints thus show that even when analyzing the data in this largely
modelindependent way, the possibility of hiding extra energy components beyond
$\Lambda$CDM in the early universe is seriously constrained by current
observations. We also find that the standard ruler, the sound horizon at
radiation drag, can be well determined in a way that does not depend on
latetime Universe assumptions, but depends strongly on earlytime physics and
in particular on additional components that behave like radiation. We find that
the standard ruler length determined in this way is $r_{\rm s} = 147.4 \pm 0.7$
Mpc if the radiation and neutrino components are standard, but the uncertainty
increases by an order of magnitude when nonstandard dark radiation components
are allowed, to $r_{\rm s} = 150 \pm 5$ Mpc.

We measure and analyse the bispectrum of the final, Data Release 12, galaxy
sample provided by the Baryon Oscillation Spectroscopic Survey, splitting by
selection algorithm into LOWZ and CMASS galaxies. The LOWZ sample contains
361\,762 galaxies with an effective redshift of $z_{\rm LOWZ}=0.32$, and the
CMASS sample 777\,202 galaxies with an effective redshift of $z_{\rm
CMASS}=0.57$. Combining the power spectrum, measured relative to the
lineofsight, with the spherically averaged bispectrum, we are able to
constrain the product of the growth of structure parameter, $f$, and the
amplitude of dark matter density fluctuations, $\sigma_8$, along with the
geometric AlcockPaczynski parameters, the product of the Hubble constant and
the comoving sound horizon at the baryon drag epoch, $H(z)r_s(z_d)$, and the
angular distance parameter divided by the sound horizon, $D_A(z)/r_s(z_d)$.
After combining prereconstruction RSD analyses of the power spectrum monopole,
quadrupole and bispectrum monopole; with postreconstruction analysis of the
BAO power spectrum monopole and quadrupole, we find $f(z_{\rm
LOWZ})\sigma_8(z_{\rm LOWZ})=0.427\pm 0.056$, $D_A(z_{\rm LOWZ})/r_s(z_d)=6.60
\pm 0.13$, $H(z_{\rm LOWZ})r_s(z_d)=(11.55\pm 0.38)10^3\,{\rm kms}^{1}$ for
the LOWZ sample, and $f(z_{\rm CMASS})\sigma_8(z_{\rm CMASS})=0.426\pm 0.029$,
$D_A(z_{\rm CMASS})/r_s(z_d)=9.39 \pm 0.10$, $H(z_{\rm
CMASS})r_s(z_d)=(14.02\pm 0.22)10^3\,{\rm kms}^{1}$ for the CMASS sample. We
find general agreement with previous BOSS DR11 and DR12 measurements. Combining
our dataset with {\it Planck15} we perform a null test of General Relativity
(GR) through the $\gamma$parametrisation finding
$\gamma=0.733^{+0.068}_{0.069}$, which is $\sim2.7\sigma$ away from the GR
predictions.

We perform a comprehensive cosmological study of the $H_0$ tension between
the direct local measurement and the modeldependent value inferred from the
Cosmic Microwave Background. With the recent measurement of $H_0$ this tension
has raised to more than $3\sigma$. We consider changes in the early time
physics without modifying the late time cosmology. We also reconstruct the late
time expansion history in a model independent way with minimal assumptions
using distances measures from Baryon Acoustic Oscillations and Type Ia
Supernovae, finding that at $z<0.6$ the recovered shape of the expansion
history is less than 5 % different than that of a standard LCDM model. These
probes also provide a model insensitive constraint on the lowredshift standard
ruler, measuring directly the combination $r_s h$ where $H_0=h \times 100$
km/s/Mpc and $r_s$ is the sound horizon at radiation drag (the standard ruler),
traditionally constrained by CMB observations. Thus $r_s$ and $H_0$ provide
absolute scales for distance measurements (anchors) at opposite ends of the
observable Universe. We calibrate the cosmic distance ladder and obtain a
modelindependent determination of the standard ruler for acoustic scale,
$r_s$. The tension in $H_0$ reflects a mismatch between our determination of
$r_s$ and its standard, CMBinferred value. Without including highl Planck CMB
polarization data (i.e., only considering the "recommended baseline" lowl
polarisation and temperature and the high l temperature data), a modification
of the earlytime physics to include a component of dark radiation with an
effective number of species around 0.4 would reconcile the CMBinferred
constraints, and the local $H_0$ and standard ruler determinations. The
inclusion of the "preliminary" highl Planck CMB polarisation data disfavours
this solution.

Assuming the existence of standard rulers, standard candles and standard
clocks, requiring only the cosmological principle, a metric theory of gravity,
a smooth expansion history, and using stateoftheart observations, we
determine the length of the "lowredshift standard ruler". The data we use are
a compilation of recent Baryon acoustic oscillation data (relying on the
standard ruler), Type 1A supernov\ae\ (as standard candles), ages of early type
galaxies (as standard clocks) and local determinations of the Hubble constant
(as a local anchor of the cosmic distance scale). In a standard $\Lambda$CDM
cosmology the "lowredshift standard ruler" coincides with the sound horizon at
radiation drag, which can also be determined in a model dependent way from
CMB observations. However, in general, the two quantities need not coincide. We
obtain constraints on the length of the lowredshift standard ruler: $r^h_{\rm
s}=101.0 \pm 2.3 h^{1}$ Mpc, when using only Type 1A supernov\ae\ and Baryon
acoustic oscillations, and $r_{\rm s}=150.0\pm 4.7 $ Mpc when using clocks to
set the Hubble normalisation, while $r_{\rm s}=141.0\pm 5.5 $ Mpc when using
the local Hubble constant determination (using both yields $r_{\rm s}=143.9\pm
3.1 $ Mpc).
The lowredshift determination of the standard ruler has an error which is
competitive with the modeldependent determination from cosmic microwave
background measurements made with the {\em Planck} satellite, which assumes it
is the sound horizon at the end of baryon drag.

We present cosmological results from the final galaxy clustering data set of
the Baryon Oscillation Spectroscopic Survey, part of the Sloan Digital Sky
Survey III. Our combined galaxy sample comprises 1.2 million massive galaxies
over an effective area of 9329 deg^2 and volume of 18.7 Gpc^3, divided into
three partially overlapping redshift slices centred at effective redshifts
0.38, 0.51, and 0.61. We measure the angular diameter distance DM and Hubble
parameter H from the baryon acoustic oscillation (BAO) method after applying
reconstruction to reduce nonlinear effects on the BAO feature. Using the
anisotropic clustering of the prereconstruction density field, we measure the
product DM*H from the AlcockPaczynski (AP) effect and the growth of structure,
quantified by f{\sigma}8(z), from redshiftspace distortions (RSD). We combine
measurements presented in seven companion papers into a set of consensus values
and likelihoods, obtaining constraints that are tighter and more robust than
those from any one method. Combined with Planck 2015 cosmic microwave
background measurements, our distance scale measurements simultaneously imply
curvature {\Omega}_K =0.0003+/0.0026 and a dark energy equation of state
parameter w = 1.01+/0.06, in strong affirmation of the spatially flat cold
dark matter model with a cosmological constant ({\Lambda}CDM). Our RSD
measurements of f{\sigma}_8, at 6 per cent precision, are similarly consistent
with this model. When combined with supernova Ia data, we find H0 = 67.3+/1.0
km/s/Mpc even for our most general dark energy model, in tension with some
direct measurements. Adding extra relativistic species as a degree of freedom
loosens the constraint only slightly, to H0 = 67.8+/1.2 km/s/Mpc. Assuming
flat {\Lambda}CDM we find {\Omega}_m = 0.310+/0.005 and H0 = 67.6+/0.5
km/s/Mpc, and we find a 95% upper limit of 0.16 eV/c^2 on the neutrino mass
sum.

We demonstrate that a scalar field is unable to reverse its direction of
motion while continuously exchanging energy with another fluid. If the rate of
transfer is modulated by the scalar's acceleration, the field can undergo a
rapid process of freezing, despite being displaced from the local minimum of
its potential. This enables dark energy to form from any potential, regardless
of its shape. The field's equation of state mimicks that of a cosmological
constant. We present a physically motivated realisation in the form of a
derivative neutrinomajoron coupling. Coherent motions, which form only once
the neutrinos become nonrelativistic, could be responsible for instigating the
freezing process. This would provide a natural resolution to the dark energy
coincidence problem, while avoiding the dynamical instabilities associated with
massvarying neutrino models. Finally we discuss possible means by which this
model could be experimentally verified.

Recent anomalies found in cosmological datasets such as the low multipoles of
the Cosmic Microwave Background or the low redshift amplitude and growth of
clustering measured by e.g., abundance of galaxy clusters and redshift space
distortions in galaxy surveys, have motivated explorations of models beyond
standard $\Lambda$CDM. Of particular interest are models where general
relativity (GR) is modified on large cosmological scales. Here we consider
deviations from $\Lambda$CDM+GR within the context of Horndeski gravity, which
is the most general theory of gravity with second derivatives in the equations
of motion. We adopt a parametrization in which the four additional Horndeski
functions of time $\alpha_i(t)$ are proportional to the cosmological density of
dark energy $\Omega_{DE}(t)$. Constraints on this extended parameter space
using a suite of stateofthe art cosmological observations are presented for
the first time. Although the theory is able to accommodate the low multipoles
of the Cosmic Microwave Background and the low amplitude of fluctuations from
redshift space distortions, we find no significant tension with $\Lambda$CDM+GR
when performing a global fit to recent cosmological data and thus there is no
evidence against $\Lambda$CDM+GR from an analysis of the value of the Bayesian
evidence ratio of the modified gravity models with respect to $\Lambda$CDM,
despite introducing extra parameters. The posterior distribution of these extra
parameters that we derive return strong constraints on any possible deviations
from $\Lambda$CDM+GR in the context of Horndeski gravity. We illustrate how our
results can be applied to a more general frameworks of modified gravity models.

Euclid is a European Space Agency medium class mission selected for launch in
2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to
understand the origin of the accelerated expansion of the universe. Euclid will
explore the expansion history of the universe and the evolution of cosmic
structures by measuring shapes and redshifts of galaxies as well as the
distribution of clusters of galaxies over a large fraction of the sky. Although
the main driver for Euclid is the nature of dark energy, Euclid science covers
a vast range of topics, from cosmology to galaxy evolution to planetary
research. In this review we focus on cosmology and fundamental physics, with a
strong emphasis on science beyond the current standard models. We discuss five
broad topics: dark energy and modified gravity, dark matter, initial
conditions, basic assumptions and questions of methodology in the data
analysis. This review has been planned and carried out within Euclid's Theory
Working Group and is meant to provide a guide to the scientific themes that
will underlie the activity of the group during the preparation of the Euclid
mission.

We present a minimally parametric, model independent reconstruction of the
shape of the primordial power spectrum. Our smoothing spline technique is
wellsuited to search for smooth features such as deviations from scale
invariance, and deviations from a power law such as running of the spectral
index or smallscale power suppression. We use a comprehensive set of the
stateof the art cosmological data: {\it Planck} observations of the
temperature and polarisation anisotropies of the cosmic microwave background,
WiggleZ and Sloan Digital Sky Survey Data Release 7 galaxy power spectra and
the CanadaFranceHawaii Lensing Survey correlation function. This
reconstruction strongly supports the evidence for a power law primordial power
spectrum with a red tilt and disfavours deviations from a power law power
spectrum including smallscale power suppression such as that induced by
significantly massive neutrinos. This offers a powerful confirmation of the
inflationary paradigm, justifying the adoption of the inflationary prior in
cosmological analyses.

Deriving the expansion history of the Universe is a major goal of modern
cosmology. To date, the most accurate measurements have been obtained with Type
Ia Supernovae and Baryon Acoustic Oscillations, providing evidence for the
existence of a transition epoch at which the expansion rate changes from
decelerated to accelerated. However, these results have been obtained within
the framework of specific cosmological models that must be implicitly or
explicitly assumed in the measurement. It is therefore crucial to obtain
measurements of the accelerated expansion of the Universe independently of
assumptions on cosmological models. Here we exploit the unprecedented
statistics provided by the Baryon Oscillation Spectroscopic Survey (BOSS) Data
Release 9 to provide new constraints on the Hubble parameter $H(z)$ using the
em cosmic chronometers approach. We extract a sample of more than 130000 of the
most massive and passively evolving galaxies, obtaining five new
cosmologyindependent $H(z)$ measurements in the redshift range $0.3<z<0.5$,
with an accuracy of $\sim$1116\% incorporating both statistical and systematic
errors. Once combined, these measurements yield a 6\% accuracy constraint of
$H(z=0.4293)=91.8\pm5.3$ km/s/Mpc. The new data are crucial to provide the
first cosmologyindependent determination of the transition redshift at high
statistical significance, measuring $z_{t}=0.4\pm0.1$, and to significantly
disfavor the null hypothesis of no transition between decelerated and
accelerated expansion at 99.9\% confidence level. This analysis highlights the
wide potential of the cosmic chronometers approach: it permits to derive
constraints on the expansion history of the Universe with results competitive
with standard probes, and most importantly, being the estimates independent of
the cosmological model, it can constrain cosmologies beyond and including the
$\Lambda$CDM model.