
A persistent theme in the study of dark energy is the question of whether it
really exists or not. It is often claimed hat we are miscalculating the
cosmological model by neglecting the effects associated with averaging over
largescale structures. In the Newtonian approximation this is clear: there is
no effect. Within the full relativistic picture this remains an important open
question however, owing to the complex mathematics involved. We study this
issue using particle numerical simulations which account for all relevant
relativistic effects without any problems from shell crossing. In this context
we show for the first time that the backreaction from structure can differ by
many orders of magnitude depending upon the slicing of spacetime one chooses to
average over. In the worst case, where smoothing is carried out in synchronous
spatial surfaces, the corrections can reach ten percent and more. However, when
smoothing on the constant time hypersurface of the Newtonian gauge,
backreaction contributions remain 35 orders of magnitude smaller.

We give an analytical interpretation of how subsamplebased internal
covariance estimators lead to biased estimates of the covariance, due to
underestimating the supersample covariance (SSC). This includes the jackknife
and bootstrap methods as estimators for the full survey area, and subsampling
as an estimator of the covariance of subsamples. The limitations of the
jackknife covariance have been previously presented in the literature because
it is effectively a rescaling of the covariance of the subsample area. However
we point out that subsampling is also biased, but for a different reason: the
subsamples are not independent, and the corresponding lack of power results in
SSC underprediction. We develop the formalism in the case of cluster counts
that allows the bias of each covariance estimator to be exactly predicted. We
find significant effects for a smallscale area or when a low number of
subsamples is used, with autoredshift biases ranging from 0.4% to 15% for
subsampling and from 5% to 75% for jackknife covariance estimates. The
crossredshift covariance is even more affected; biases range from 8% to 25%
for subsampling and from 50% to 90% for jackknife. Owing to the redshift
evolution of the probe, the covariances cannot be debiased by a simple
rescaling factor, and an exact debiasing has the same requirements as the full
SSC prediction. These results thus disfavour the use of internal covariance
estimators on data itself or a single simulation, leaving analytical prediction
and simulations suites as possible SSC predictors.

We have performed a detailed xray diffraction structural study of Zr under
pressure and unambiguously identify the existence of a firstorder
isostructural bcctobcc phase transition near 58 GPa. Firstprinciples quantum
molecular dynamics lattice dynamics calculations support the existence of this
phase transition, in excellent agreement with experimental results, triggered
by anharmonic effects. Our results highlight the potential ubiquity of
anharmonically driven isostructural transitions within the periodic table under
pressure and calls for followup experimental and theoretical studies.

The recent discovery of a $\gamma$ray counterpart to a gravitational wave
(GW) event has put extremely stringent constraints on the speed of
gravitational waves at the present epoch. In turn, these constraints place
strong theoretical pressure on potential modifications of gravity, essentially
allowing only the conformal sector to be active in the present Universe. In
this paper, we show that direct detection of gravitational waves from optically
identified sources can also measure or constrain the conformal sector of
modified gravity models through the time variation of the Planck mass. As a
first rough estimate, we find that the LISA satellite can measure the
dimensionless time variation of the Planck mass (the socalled parameter
$\alpha_M$) at redshift around 1.5 with an error of about 0.03 to 0.13,
depending on the assumptions concerning future observations. Stronger
constraints can be achieved once reliable distance indicators at $z>2$ are
developed, or with GW detectors that extend the capabilities of LISA, like the
proposed Big Bang Observer. We emphasize that, just like the constraints on the
gravitational speed, the bound on $\alpha_M$ is independent of the cosmological
model.

We investigate cosmic string networks in the Abelian Higgs model using data
from a campaign of largescale numerical simulations on lattices of up to
$4096^3$ grid points. We observe scaling or selfsimilarity of the networks
over a wide range of scales, and estimate the asymptotic values of the mean
string separation in horizon length units $\dot{\xi}$ and of the mean square
string velocity $\bar v^2$ in the continuum and large time limits. The scaling
occurs because the strings lose energy into classical radiation of the scalar
and gauge fields of the Abelian Higgs model. We quantify the energy loss with a
dimensionless radiative efficiency parameter, and show that it does not vary
significantly with lattice spacing or string separation. This implies that the
radiative energy loss underlying the scaling behaviour is not a lattice
artefact, and justifies the extrapolation of measured network properties to
large times for computations of cosmological perturbations. We also show that
the core growth method, which increases the defect core width with time to
extend the dynamic range of simulations, does not introduce significant
systematic error. We compare $\dot{\xi}$ and $\bar v^2$ to values measured in
simulations using the NambuGoto approximation, finding that the latter
underestimate the mean string separation by about 25%, and overestimate $\bar
v^2$ by about 10%. The scaling of the string separation implies that string
loops decay by the emission of massive radiation within a Hubble time in field
theory simulations, in contrast to the NambuGoto scenario which neglects this
energy loss mechanism. String loops surviving for only one Hubble time emit
much less gravitational radiation than in the NambuGoto scenario, and are
consequently subject to much weaker gravitational wave constraints on their
tension.

Some of the dark matter in the Universe is made up of massive neutrinos.
Their impact on the formation of large scale structure can be used to determine
their absolute mass scale from cosmology, but to this end accurate numerical
simulations have to be developed. Due to their relativistic nature, neutrinos
pose additional challenges when one tries to include them in Nbody simulations
that are traditionally based on Newtonian physics. Here we present the first
numerical study of massive neutrinos that uses a fully relativistic approach.
Our Nbody code, gevolution, is based on a weakfield formulation of general
relativity that naturally provides a selfconsistent framework for relativistic
particle species. This allows us to model neutrinos from first principles,
without invoking any adhoc recipes. Our simulation suite comprises some of the
largest neutrino simulations performed to date. We study the effect of massive
neutrinos on the nonlinear power spectra and the halo mass function, focusing
on the interesting mass range between 0.06 eV and 0.3 eV and including a case
for an inverted mass hierarchy.

Lensing of the CMB is now a welldeveloped probe of largescale clustering
over a broad range of redshifts. By exploiting the nonGaussian imprints of
lensing in the polarization of the CMB, the CORE mission can produce a clean
map of the lensing deflections over nearly the fullsky. The number of highS/N
modes in this map will exceed current CMB lensing maps by a factor of 40, and
the measurement will be samplevariance limited on all scales where linear
theory is valid. Here, we summarise this mission product and discuss the
science that it will enable. For example, the summed mass of neutrinos will be
determined to an accuracy of 17 meV combining CORE lensing and CMB twopoint
information with contemporaneous BAO measurements, three times smaller than the
minimum total mass allowed by neutrino oscillations. In the search for Bmode
polarization from primordial gravitational waves with CORE, lensinduced
Bmodes will dominate over instrument noise, limiting constraints on the
gravitational wave power spectrum amplitude. With lensing reconstructed by
CORE, one can "delens" the observed polarization internally, reducing the
lensing Bmode power by 60%. This improves to 70% by combining lensing and CIB
measurements from CORE, reducing the error on the gravitational wave amplitude
by 2.5 compared to no delensing (in the null hypothesis). Lensing measurements
from CORE will allow calibration of the halo masses of the 40000 galaxy
clusters that it will find, with constraints dominated by the clean
polarizationbased estimators. CORE can accurately remove Galactic emission
from CMB maps with its 19 frequency channels. We present initial findings that
show that residual Galactic foreground contamination will not be a significant
source of bias for lensing power spectrum measurements with CORE. [abridged]

We forecast the scientific capabilities to improve our understanding of
cosmic inflation of CORE, a proposed CMB space satellite submitted in response
to the ESA fifth call for a mediumsize mission opportunity. The CORE satellite
will map the CMB anisotropies in temperature and polarization in 19 frequency
channels spanning the range 60600 GHz. CORE will have an aggregate noise
sensitivity of $1.7 \mu$K$\cdot \,$arcmin and an angular resolution of 5' at
200 GHz. We explore the impact of telescope size and noise sensitivity on the
inflation science return by making forecasts for several instrumental
configurations. This study assumes that the lower and higher frequency channels
suffice to remove foreground contaminations and complements other related
studies of component separation and systematic effects, which will be reported
in other papers of the series "Exploring Cosmic Origins with CORE." We forecast
the capability to determine key inflationary parameters, to lower the detection
limit for the tensortoscalar ratio down to the $10^{3}$ level, to chart the
landscape of single field slowroll inflationary models, to constrain the epoch
of reheating, thus connecting inflation to the standard radiationmatter
dominated Big Bang era, to reconstruct the primordial power spectrum, to
constrain the contribution from isocurvature perturbations to the $10^{3}$
level, to improve constraints on the cosmic string tension to a level below the
presumptive GUT scale, and to improve the current measurements of primordial
nonGaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the
models explored, CORE alone will improve significantly on the present
constraints on the physics of inflation. Its capabilities will be further
enhanced by combining with complementary future cosmological observations.

We forecast the main cosmological parameter constraints achievable with the
CORE space mission which is dedicated to mapping the polarisation of the Cosmic
Microwave Background (CMB). CORE was recently submitted in response to ESA's
fifth call for mediumsized mission proposals (M5). Here we report the results
from our presubmission study of the impact of various instrumental options, in
particular the telescope size and sensitivity level, and review the great,
transformative potential of the mission as proposed. Specifically, we assess
the impact on a broad range of fundamental parameters of our Universe as a
function of the expected CMB characteristics, with other papers in the series
focusing on controlling astrophysical and instrumental residual systematics. In
this paper, we assume that only a few central CORE frequency channels are
usable for our purpose, all others being devoted to the cleaning of
astrophysical contaminants. On the theoretical side, we assume LCDM as our
general framework and quantify the improvement provided by CORE over the
current constraints from the Planck 2015 release. We also study the joint
sensitivity of CORE and of future Baryon Acoustic Oscillation and Large Scale
Structure experiments like DESI and Euclid. Specific constraints on the physics
of inflation are presented in another paper of the series. In addition to the
six parameters of the base LCDM, which describe the matter content of a
spatially flat universe with adiabatic and scalar primordial fluctuations from
inflation, we derive the precision achievable on parameters like those
describing curvature, neutrino physics, extra light relics, primordial helium
abundance, dark matter annihilation, recombination physics, variation of
fundamental constants, dark energy, modified gravity, reionization and cosmic
birefringence. (ABRIDGED)

We reanalyse recent Cepheid data to estimate the Hubble parameter $H_0$ by
using Bayesian hyperparameters (HPs). We consider the two data sets from Riess
et al 2011 and 2016 (labelled R11 and R16, with R11 containing less than half
the data of R16) and include the available anchor distances (megamaser system
NGC4258, detached eclipsing binary distances to LMC and M31, and MW Cepheids
with parallaxes), use a weak metallicity prior and no period cut for Cepheids.
We find that part of the R11 data is downweighted by the HPs but that R16 is
mostly consistent with expectations for a Gaussian distribution, meaning that
there is no need to downweight the R16 data set. For R16, we find a value of
$H_0 = 73.75 \pm 2.11 \, \mathrm{km} \, \mathrm{s}^{1} \, \mathrm{Mpc}^{1}$
if we use HPs for all data points (including Cepheid stars, supernovae type Ia,
and the available anchor distances), which is about 2.6 $\sigma$ larger than
the Planck 2015 value of $H_0 = 67.81 \pm 0.92 \,\mathrm{km}\, \mathrm{s}^{1}
\, \mathrm{Mpc}^{1}$ and about 3.1 $\sigma$ larger than the updated Planck
2016 value $66.93 \pm 0.62 \,\mathrm{km}\, \mathrm{s}^{1} \,
\mathrm{Mpc}^{1}$. We test the effect of different assumptions, and find that
the choice of anchor distances affects the final value significantly. If we
exclude the Milky Way from the anchors, then the value of $H_0$ decreases. We
find however no evident reason to exclude the MW data. The HP method used here
avoids subjective rejection criteria for outliers and offers a way to test
datasets for unknown systematics.

Modified Gravity theories generally affect the Poisson equation and the
gravitational slip (effective anisotropic stress) in an observable way, that
can be parameterized by two generic functions ($\eta$ and $\mu$) of time and
space. We bin the time dependence of these functions in redshift and present
forecasts on each bin for future surveys like Euclid. We consider both Galaxy
Clustering and Weak Lensing surveys, showing the impact of the nonlinear
regime, treated with two different semianalytical approximations. In addition
to these future observables, we use a prior covariance matrix derived from the
Planck observations of the Cosmic Microwave Background. Our results show that
$\eta$ and $\mu$ in different redshift bins are significantly correlated, but
including nonlinear scales reduces or even eliminates the correlation,
breaking the degeneracy between Modified Gravity parameters and the overall
amplitude of the matter power spectrum. We further decorrelate parameters with
a Zerophase Component Analysis and identify which combinations of the Modified
Gravity parameter amplitudes, in different redshift bins, are best constrained
by future surveys. We also extend the analysis to two particular
parameterizations of the time evolution of $\mu$ and $\eta$ and consider, in
addition to Euclid, also SKA1, SKA2, DESI: we find in this case that future
surveys will be able to constrain the current values of $\eta$ and $\mu$ at the
$25\%$ level when using only linear scales (wavevector k < 0.15 h/Mpc),
depending on the specific time parameterization; sensitivity improves to about
$1\%$ when nonlinearities are included.

We present cosmic microwave background (CMB) power spectra from recent
numerical simulations of cosmic strings in the Abelian Higgs model and compare
them to CMB power spectra measured by Planck. We obtain revised constraints on
the cosmic string tension parameter $G\mu$. For example, in the $\Lambda$CDM
model with the addition of strings and no primordial tensor perturbations, we
find $G\mu < 2.0 \times 10^{7}$ at 95% confidence, about 20% lower than the
value obtained from previous simulations, which had 1/64 of the spatial volume.
We investigate the source of the difference, showing that the main cause is an
improved treatment of the string evolution across the radiationmatter
transition. The increased computational volume also makes possible to simulate
fully the physical equations of motion, in which the string cores shrink in
comoving coordinates. This, and the larger dynamic range, changes the amplitude
of the power spectra by only about 10%, demonstrating that field theory
simulations of cosmic strings have now reached the required dynamic range for
CMB calculations.

We report on the energymomentum correlators obtained with recent numerical
simulations of the Abelian Higgs model, essential for the computation of cosmic
microwave background and matter perturbations of cosmic strings. Due to
significant improvements both in raw computing power and in our parallel
simulation framework, the dynamical range of the simulations has increased
fourfold both in space and time, and for the first time we are able to
simulate strings with a constant physical width in both the radiation and
matter eras. The new simulations improve the accuracy of the measurements of
the correlation functions at the horizon scale and confirm the shape around the
peak. The normalization is slightly higher in the high wavenumber tails, due
to a small increase in the string density. We study for the first time the
behaviour of the correlators across cosmological transitions, and discover that
the correlation functions evolve adiabatically, ie the network adapts quickly
to changes in the expansion rate. We propose a new method for constructing
source functions for EinsteinBoltzmann integrators, comparing it with two
other methods previously used. The new method is more consistent, easier to
implement, and significantly more accurate.

We demonstrate the importance of including the lensing contribution in galaxy
clustering analyses with large galaxy redshift surveys. It is well known that
radial crosscorrelations between different redshift bins of galaxy surveys are
dominated by lensing. But we show here that also neglecting lensing in the
autocorrelations within one bin severely biases cosmological parameter
estimation with redshift surveys. It leads to significant shifts for several
cosmological parameters, most notably the scalar spectral index and the
neutrino mass scale. Especially the latter parameter is one of the main targets
of future galaxy surveys.

We present a new Nbody code, gevolution, for the evolution of large scale
structure in the Universe. Our code is based on a weak field expansion of
General Relativity and calculates all six metric degrees of freedom in Poisson
gauge. Nbody particles are evolved by solving the geodesic equation which we
write in terms of a canonical momentum such that it remains valid also for
relativistic particles. We validate the code by considering the Schwarzschild
solution and, in the Newtonian limit, by comparing with the Newtonian Nbody
codes Gadget2 and RAMSES. We then proceed with a simulation of large scale
structure in a Universe with massive neutrinos where we study the gravitational
slip induced by the neutrino shear stress. The code can be extended to include
different kinds of dark energy or modified gravity models and going beyond the
usually adopted quasistatic approximation. Our code is publicly available.

We use largescale cosmological observations to place constraints on the
darkmatter pressure, sound speed and viscosity, and infer a limit on the mass
of warmdarkmatter particles. Measurements of the cosmic microwave background
(CMB) anisotropies constrain the equation of state and sound speed of the dark
matter at last scattering at the per mille level. Since the redshifting of
collisionless particles universally implies that these quantities scale like
$a^{2}$ absent shell crossing, we infer that today $w_{\rm (DM)}< 10^{10.0}$,
$c_{\rm s,(DM)}^2 < 10^{10.7}$ and $c_{\rm vis, (DM)}^{2} < 10^{10.3}$ at the
$99\%$ confidence level. This very general bound can be translated to
modeldependent constraints on darkmatter models: for warm dark matter these
constraints imply $m> 70$ eV, assuming it decoupled while relativistic around
the same time as the neutrinos; for a cold relic, we show that $m>100$ eV. We
separately constrain the properties of the DM fluid on linear scales at late
times, and find upper bounds $c_{\rm s, (DM)}^2<10^{5.9}$, $c_{\rm vis,
(DM)}^{2} < 10^{5.7}$, with no detection of nondust properties for the DM.

The angular power spectrum of the cosmic infrared background (CIB) is a
sensitive probe of the local primordial bispectrum. CIB measurements are
integrated over a large volume so that the scale dependent bias from the
primordial nonGaussianity leaves a strong signal in the CIB power spectrum.
Although galactic dust dominates over the nonGaussian CIB signal, it is
possible to mitigate the dust contamination with enough frequency channels,
especially if high frequencies such as the Planck 857 GHz channel are
available. We show that, in this case, measurements of the cosmic microwave
background from future space missions should be able to probe the local
bispectrum shape down to an amplitude f_nl < 1.

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 comprehensive and updated comparison with cosmological
observations of two nonlocal modifications of gravity previously introduced by
our group, the so called RR and RT models. We implement the background
evolution and the cosmological perturbations of the models in a modified
Boltzmann code, using CLASS. We then test the nonlocal models against the {\em
Planck} 2015 TT, TE, EE and Cosmic Microwave Background (CMB) lensing data,
isotropic and anisotropic Baryonic Acoustic Oscillations (BAO) data, JLA
supernovae, $H_0$ measurements and growth rate data, and we perform Bayesian
parameter estimation. We then compare the RR, RT and $\Lambda$CDM models, using
the SavageDickey method. We find that the RT model and $\Lambda$CDM perform
equally well, while the performance of the RR model with respect to
$\Lambda$CDM depends on whether or not we include a prior on $H_0$ based on
local measurements.

Euclid is a European Space Agency medium class mission selected for launch in
2019 within the Cosmic Vision 20152025 programme. 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.

Radio interferometers suffer from the problem of missing information in their
data, due to the gaps between the antennas. This results in artifacts, such as
bright rings around sources, in the images obtained. Multiple deconvolution
algorithms have been proposed to solve this problem and produce cleaner radio
images. However, these algorithms are unable to correctly estimate
uncertainties in derived scientific parameters or to always include the effects
of instrumental errors. We propose an alternative technique called Bayesian
Inference for Radio Observations (BIRO) which uses a Bayesian statistical
framework to determine the scientific parameters and instrumental errors
simultaneously directly from the raw data, without making an image. We use a
simple simulation of Westerbork Synthesis Radio Telescope data including
pointing errors and beam parameters as instrumental effects, to demonstrate the
use of BIRO.

Numerical simulations are a versatile tool providing insight into the
complicated process of structure formation in cosmology. This process is mainly
governed by gravity, which is the dominant force on large scales. To date, a
century after the formulation of general relativity, numerical codes for
structure formation still employ Newton's law of gravitation. This
approximation relies on the two assumptions that gravitational fields are weak
and that they are only sourced by nonrelativistic matter. While the former
appears well justified on cosmological scales, the latter imposes restrictions
on the nature of the "dark" components of the Universe (dark matter and dark
energy) which are, however, poorly understood. Here we present the first
simulations of cosmic structure formation using equations consistently derived
from general relativity. We study in detail the small relativistic effects for
a standard {\Lambda}CDM cosmology which cannot be obtained within a purely
Newtonian framework. Our particlemesh Nbody code computes all six degrees of
freedom of the metric and consistently solves the geodesic equation for
particles, taking into account the relativistic potentials and the
framedragging force. This conceptually clean approach is very general and can
be applied to various settings where the Newtonian approximation fails or
becomes inaccurate, ranging from simulations of models with dynamical dark
energy or warm/hot dark matter to core collapse supernova explosions.

From the galaxy power spectrum in redshift space, we derive semianalytical
results on the generic degeneracy of galaxy clustering measurements. Defining
the observables $\bar{A}= Gb\sigma_8$ and $\bar{R} = Gf\sigma_8$, (being $G$
the growth function, $b$ the bias, $f$ the growth rate, and $\sigma_8$ the
amplitude of the power spectrum), we perform a Fisher matrix formalism to
forecast the expected precision of these quantities for a Euclidlike survey.
Among the results we found that galaxy surveys have generically a slightly
negative correlation between $\bar{A}$ and $\bar{R}$, and they can always
measure $\bar{R}$ about 3.7 to 4.7 times better than $\bar{A}$.

New telescopes like the Square Kilometre Array (SKA) will push into a new
sensitivity regime and expose systematics, such as directiondependent effects,
that could previously be ignored. Current methods for handling such systematics
rely on alternating best estimates of instrumental calibration and models of
the underlying sky, which can lead to inadequate uncertainty estimates and
biased results because any correlations between parameters are ignored. These
deconvolution algorithms produce a single image that is assumed to be a true
representation of the sky, when in fact it is just one realization of an
infinite ensemble of images compatible with the noise in the data. In contrast,
here we report a Bayesian formalism that simultaneously infers both systematics
and science. Our technique, Bayesian Inference for Radio Observations (BIRO),
determines all parameters directly from the raw data, bypassing imagemaking
entirely, by sampling from the joint posterior probability distribution. This
enables it to derive both correlations and accurate uncertainties, making use
of the flexible software MEQTREES to model the sky and telescope
simultaneously. We demonstrate BIRO with two simulated sets of Westerbork
Synthesis Radio Telescope data sets. In the first, we perform joint estimates
of 103 scientific (flux densities of sources) and instrumental (pointing
errors, beamwidth and noise) parameters. In the second example, we perform
source separation with BIRO. Using the Bayesian evidence, we can accurately
select between a single point source, two point sources and an extended
Gaussian source, allowing for 'superresolution' on scales much smaller than
the synthesized beam.

Recent work has demonstrated that it is important to constrain the dynamics
of cosmological perturbations, in addition to the evolution of the background,
if we want to distinguish among different models of the dark sector. Especially
the anisotropic stress of the (possibly effective) dark energy fluid has been
shown to be an important discriminator between modified gravity and dark energy
models. In this paper we use approximate analytical solutions of the
perturbation equations in the presence of viscosity to study how the
anisotropic stress affects the weak lensing and galaxy power spectrum. We then
forecast how sensitive the photometric and spectroscopic Euclid surveys will be
to both the speed of sound and the viscosity of our effective dark energy fluid
when using weak lensing tomography and the galaxy power spectrum. We find that
Euclid alone can only constrain models with very small speed of sound and
viscosity, while it will need the help of other observables in order to give
interesting constraints on models with a sound speed close to one. This
conclusion is also supported by the expected Bayes factor between models.