
We determine an optimized clustering statistic to be used for galaxy samples
with significant redshift uncertainty, such as those that rely on photometric
redshifts. To do so, we study the BAO information content as a function of the
orientation of galaxy clustering modes with respect to their angle to the
lineofsight (LOS). The clustering along the LOS, as observed in a
redshiftspace with significant redshift uncertainty, has contributions from
clustering modes with a range of orientations with respect to the true LOS. For
redshift uncertainty $\sigma_z \geq 0.02(1+z)$ we find that while the BAO
information is confined to transverse clustering modes in the true space, it is
spread nearly evenly in the observed space. Thus, measuring clustering in terms
of the projected separation (regardless of the LOS) is an efficient and nearly
lossless compression of the signal for $\sigma_z \geq 0.02(1+z)$. For reduced
redshift uncertainty, a more careful consideration is required. We then use
more than 1700 realizations of galaxy simulations mimicking the Dark Energy
Survey Year 1 sample to validate our analytic results and optimized analysis
procedure. We find that using the correlation function binned in projected
separation, we can achieve uncertainties that are within 10 per cent of of
those predicted by Fisher matrix forecasts. We predict that DES Y1 should
achieve a 5 per cent distance measurement using our optimized methods. We
expect the results presented here to be important for any future BAO
measurements made using photometric redshift data.

We derive constraints on feedback by active galactic nuclei (AGN) by setting
limits on their thermal SunyaevZel'dovich (SZ) imprint on the cosmic microwave
background (CMB). The amplitude of any SZ signature is small and degenerate
with the poorly known submm spectral energy distribution of the AGN host
galaxy and other unresolved dusty sources along the line of sight. Here we
break this degeneracy by combining microwave and submm data from Planck with
allsky farinfrared maps from the AKARI satellite. We first test our
measurement pipeline using the Sloan Digital Sky Survey (SDSS) redMaPPer
catalogue of galaxy clusters, finding a highly significant detection
($>$$20\sigma$) of the SZ effect together with correlated dust emission. We
then constrain the SZ signal associated with spectroscopically confirmed
quasistellar objects (QSOs) from SDSS data release 7 (DR7) and the Baryon
Oscillation Spectroscopic Survey (BOSS) DR12. We obtain a lowsignificance
($1.6\sigma$) hint of an SZ signal, pointing towards a mean thermal energy of
$\simeq 5 \times 10^{60}$ erg, lower than reported in some previous studies. A
comparison of our results with highresolution hydrodynamical simulations
including AGN feedback suggests QSO host masses of $M_{200c} \sim 4 \times
10^{12}~h^{1}M_\odot$, but with a large uncertainty. Our analysis provides no
conclusive evidence for an SZ signal specifically associated with AGN feedback.

We study the matter bispectrum of largescale structure by comparing the
predictions of different perturbative and phenomenological models with the full
threedimensional bispectrum from $N$body simulations estimated using modal
methods. We show that among the perturbative approaches, effective field theory
succeeds in extending the range of validity furthest on intermediate scales, at
the cost of free additional parameters. By studying the halo model, we show
that although it is satisfactory in the deeply nonlinear regime, it predicts a
deficit of power on intermediate scales, worsening at redshifts $z>0$. By
comparison with the $N$body bispectrum on those scales, we show that there is
a significant squeezed component underestimated in the halo model. On the basis
of these results, we propose a new threeshape model, based on the treelevel,
squeezed and constant bispectrum shapes we identified in the halo model; after
calibration this fits the simulations on all scales and redshifts of interest.
We extend this model further to primordial nonGaussianity of the local and
equilateral types by showing that the same shapes can be used to describe the
additional nonGaussian component in the matter bispectrum. This method
provides a HALOFITlike prototype of the bispectrum that could be used to
describe and test parameter dependencies and should be relevant for the
bispectrum of weak gravitational lensing and wider applications.

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.

We present measurements of angular cross power spectra between galaxies and
opticallyselected galaxy clusters in the final photometric sample of the Sloan
Digital Sky Survey (SDSS). We measure the auto and crosscorrelations between
galaxy and cluster samples, from which we extract the effective biases and
study the shot noise properties. We model the nonPoissonian shot noise by
introducing an effective number density of tracers and fit for this quantity.
We find that we can only describe the crosscorrelation of galaxies and galaxy
clusters, as well as the autocorrelation of galaxy clusters, on the relevant
scales using a nonPoissonian shot noise contribution.
The values of effective bias we finally measure for a volumelimited sample
are $b_{cc}=4.09 \pm 0.47$ for the cluster autocorrelation and $b_{gc}=2.15
\pm 0.09$ for the galaxycluster crosscorrelation. We find that these results
are consistent with expectations from the autocorrelations of galaxies and
clusters and are in good agreement with previous studies. The main result is
twofold: firstly we provide a measurement of the crosscorrelation of galaxies
and clusters, which can be used for further cosmological analysis, and secondly
we describe an effective treatment of the shot noise.

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

We study the matter bispectrum of the largescale structure by comparing
different perturbative and phenomenological models with measurements from
$N$body simulations obtained with a modal bispectrum estimator. Using shape
and amplitude correlators, we directly compare simulated data with theoretical
models over the full threedimensional domain of the bispectrum, for different
redshifts and scales. We review and investigate the main perturbative methods
in the literature that predict the oneloop bispectrum: standard perturbation
theory, effective field theory, resummed Lagrangian and renormalised
perturbation theory, calculating the latter also at two loops for some triangle
configurations. We find that effective field theory (EFT) succeeds in extending
the range of validity furthest into the mildly nonlinear regime, albeit at the
price of free extra parameters requiring calibration on simulations. For the
more phenomenological halo model, we confirm that despite its validity in the
deeply nonlinear regime it has a deficit of power on intermediate scales, which
worsens at higher redshifts; this issue is ameliorated, but not solved, by
combined haloperturbative models. We show from simulations that in this
transition region there is a strong squeezed bispectrum component that is
significantly underestimated in the halo model at earlier redshifts. We thus
propose a phenomenological method for alleviating this deficit, which we
develop into a simple phenomenological "threeshape" benchmark model based on
the three fundamental shapes we have obtained from studying the halo model.
When calibrated on the simulations, this threeshape benchmark model accurately
describes the bispectrum on all scales and redshifts considered, providing a
prototype bispectrum HALOFITlike methodology that could be used to describe
and test parameter dependencies.

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.

Bayesian model comparison penalizes models with more free parameters that are
allowed to vary over a wide range, and thus offers the most robust method to
decide whether some given data require new parameters. In this paper, we ask a
simple question: do current cosmological data require extensions of the
simplest singlefield inflation models? Specifically, we calculate the Bayesian
evidence of a totally anticorrelated isocurvature perturbation and a running
spectral index of the scalar curvature perturbation. These parameters are
motivated by recent claims that the observed temperature anisotropy of the
cosmic microwave background on large angular scales is too low to be compatible
with the simplest inflation models. Both a subdominant, anticorrelated cold
dark matter isocurvature component and a negative running index succeed in
lowering the largescale temperature power spectrum. We show that the
introduction of isocurvature perturbations is disfavored, whereas that of the
running spectral index is only moderately favored, even when the BICEP2 data
are included in the analysis without any foreground subtraction.

Any Dark Energy (DE) or Modified Gravity (MG) model that deviates from a
cosmological constant requires a consistent treatment of its perturbations,
which can be described in terms of an effective entropy perturbation and an
anisotropic stress. We have considered a recently proposed generic
parameterisation of DE/MG perturbations and compared it to data from the Planck
satellite and six galaxy catalogues, including temperaturegalaxy (Tg), CMB
lensinggalaxy and galaxygalaxy (gg) correlations. Combining these observables
of structure formation with tests of the background expansion allows us to
investigate the properties of DE/MG both at the background and the perturbative
level. Our constraints on DE/MG are mostly in agreement with the cosmological
constant paradigm, while we also find that the constraint on the equation of
state w (assumed to be constant) depends on the model assumed for the
perturbation evolution. We obtain $w=0.92^{+0.20}_{0.16}$ (95% CL; CMB+gg+Tg)
in the entropy perturbation scenario; in the anisotropic stress case the result
is $w=0.86^{+0.17}_{0.16}$. Including the lensing correlations shifts the
results towards higher values of w. If we include a prior on the expansion
history from recent Baryon Acoustic Oscillations (BAO) measurements, we find
that the constraints tighten closely around $w=1$, making it impossible to
measure any DE/MG perturbation evolution parameters. If, however, upcoming
observations from surveys like DES, Euclid or LSST show indications for a
deviation from a cosmological constant, our formalism will be a useful tool
towards model selection in the dark sector.

We present the clustering of galaxy clusters as a useful addition to the
common set of cosmological observables. The clustering of clusters probes the
largescale structure of the Universe, extending galaxy clustering analysis to
the highpeak, highbias regime. Clustering of galaxy clusters complements the
traditional cluster number counts and observablemass relation analyses,
significantly improving their constraining power by breaking existing
calibration degeneracies. We use the maxBCG galaxy clusters catalogue to
constrain cosmological parameters and crosscalibrate the massobservable
relation, using cluster abundances in richness bins and weaklensing mass
estimates. We then add the redshiftspace power spectrum of the sample,
including an effective modelling of the weakly nonlinear contribution and
allowing for an arbitrary photometric redshift smoothing. The inclusion of the
power spectrum data allows for an improved selfcalibration of the scaling
relation. We find that the inclusion of the power spectrum typically brings a
$\sim 50$ per cent improvement in the errors on the fluctuation amplitude
$\sigma_8$ and the matter density $\Omega_{\mathrm{m}}$. Finally, we apply this
method to constrain models of the early universe through the amount of
primordial nonGaussianity of the local type, using both the variation in the
halo mass function and the variation in the cluster bias. We find a constraint
on the amount of skewness $f_{\mathrm{NL}} = 12 \pm 157 $ ($1\sigma$) from the
cluster data alone.

We apply a new method to measure primordial nonGaussianity, using the
crosscorrelation between galaxy surveys and the CMB lensing signal to measure
galaxy bias on very large scales, where localtype primordial nonGaussianity
predicts a $k^2$ divergence. We use the CMB lensing map recently published by
the Planck collaboration, and measure its external correlations with a suite of
six galaxy catalogues spanning a broad redshift range. We then consistently
combine correlation functions to extend the recent analysis by Giannantonio et
al. (2013), where the densitydensity and the densityCMB temperature
correlations were used. Due to the intrinsic noise of the Planck lensing map,
which affects the largest scales most severely, we find that the constraints on
the galaxy bias are similar to the constraints from densityCMB temperature
correlations. Including lensing constraints only improves the previous
statistical measurement errors marginally, and we obtain $ f_{\mathrm{NL}} = 12
\pm 21 $ (1$\sigma$) from the combined data set. However, the lensing
measurements serve as an excellent test of systematic errors: we now have three
methods to measure the largescale, scaledependent bias from a galaxy survey:
autocorrelation, and crosscorrelation with both CMB temperature and lensing.
As the publicly available Planck lensing maps have had their largestscale
modes at multipoles $l<10$ removed, which are the most sensitive to the
scaledependent bias, we consider mock CMB lensing data covering all
multipoles. We find that, while the effect of $f_{\mathrm{NL}}$ indeed
increases significantly on the largest scales, so do the contributions of both
cosmic variance and the intrinsic lensing noise, so that the improvement is
small.

We present the strongest robust constraints on primordial nonGaussianity
(PNG) from currently available galaxy surveys, combining largescale clustering
measurements and their crosscorrelations with the cosmic microwave background.
We update the data sets used by Giannantonio et al. (2012), and broaden that
analysis to include the full set of twopoint correlation functions between all
surveys. In order to obtain the most reliable constraints on PNG, we advocate
the use of the crosscorrelations between the catalogs as a robust estimator
and we perform an extended analysis of the possible systematics to reduce their
impact on the results. To minimize the impact of stellar contamination in our
luminous red galaxy (LRG) sample, we use the recent Baryon Oscillations
Spectroscopic Survey catalog of Ross et al. (2011). We also find evidence for a
new systematic in the NVSS radio galaxy survey similar to, but smaller than,
the known declinationdependent issue; this is difficult to remove without
affecting the inferred PNG signal, and thus we do not include the NVSS
autocorrelation function in our analyses. We find no evidence of primordial
nonGaussianity; for the localtype configuration we obtain for the skewness
parameter $ 36 < f_{\mathrm{NL}} < 45 $ at 95 % c.l. ($5 \pm 21$ at $1\sigma$)
when using the most conservative part of our data set, improving previous
results; we also find no evidence for significant kurtosis, parameterized by
$g_{\mathrm{NL}}$. In addition to PNG, we simultaneously constrain dark energy
and find that it is required with a form consistent with a cosmological
constant.

We revisit the state of the integrated SachsWolfe (ISW) effect measurements
in light of newly available data and address criticisms about the measurements
which have recently been raised. We update the data set previously assembled by
Giannantonio et al. to include new data releases for both the cosmic microwave
background (CMB) and the largescale structure (LSS) of the Universe. We find
that our updated results are consistent with previous measurements. By fitting
a single template amplitude, we now obtain a combined significance of the ISW
detection at the 4.4 sigma level, which fluctuates by 0.4 sigma when
alternative data cuts and analysis assumptions are considered. We also make new
tests for systematic contaminations of the data, focusing in particular on the
issues raised by Sawangwit et al. Amongst them, we address the rotation test,
which aims at checking for possible systematics by correlating pairs of
randomly rotated maps. We find results consistent with the expected data
covariance, no evidence for enhanced correlation on any preferred axis of
rotation, and therefore no indication of any additional systematic
contamination. We publicly release the results, the covariance matrix, and the
sky maps used to obtain them.

We present a comparison of Fisher matrix forecasts for cosmological probes
with Monte Carlo Markov Chain (MCMC) posterior likelihood estimation methods.
We analyse the performance of future Dark Energy Task Force (DETF) stageIII
and stage IV darkenergy surveys using supernovae, baryon acoustic
oscillations and weak lensing as probes. We concentrate in particular on the
darkenergy equation of state parameters w0 and wa. For purely geometrical
probes, and especially when marginalising over wa, we find considerable
disagreement between the two methods, since in this case the Fisher matrix can
not reproduce the highly nonelliptical shape of the likelihood function. More
quantitatively, the Fisher method underestimates the marginalized errors for
purely geometrical probes between 30%70%. For cases including structure
formation such as weak lensing, we find that the posterior probability contours
from the Fisher matrix estimation are in good agreement with the MCMC contours
and the forecasted errors only changing on the 5% level. We then explore
nonlinear transformations resulting in physicallymotivated parameters and
investigate whether these parameterisations exhibit a Gaussian behaviour. We
conclude that for the purely geometrical probes and, more generally, in cases
where it is not known whether the likelihood is close to Gaussian, the Fisher
matrix is not the appropriate tool to produce reliable forecasts.

We study the constraining power on primordial nonGaussianity of future
surveys of the largescale structure of the Universe for both nearterm surveys
(such as the Dark Energy Survey  DES) as well as longer term projects such as
Euclid and WFIRST. Specifically we perform a Fisher matrix analysis forecast
for such surveys, using DESlike and Euclidlike configurations as examples,
and take account of any expected photometric and spectroscopic data. We focus
on twopoint statistics and we consider three observables: the 3D galaxy power
spectrum in redshift space, the angular galaxy power spectrum, and the
projected weaklensing shear power spectrum. We study the effects of adding a
few extra parameters to the basic LCDM set. We include the two standard
parameters to model the current value for the dark energy equation of state and
its time derivative, w_0, w_a, and we account for the possibility of primordial
nonGaussianity of the local, equilateral and orthogonal types, of parameter
fNL and, optionally, of spectral index n_fNL. We present forecasted constraints
on these parameters using the different observational probes. We show that
accounting for models that include primordial nonGaussianity does not degrade
the constraint on the standard LCDM set nor on the darkenergy equation of
state. By combining the weak lensing data and the information on projected
galaxy clustering, consistently including all twopoint functions and their
covariance, we find forecasted marginalised errors sigma (fNL) ~ 3, sigma
(n_fNL) ~ 0.12 from a Euclidlike survey for the local shape of primordial
nonGaussianity, while the orthogonal and equilateral constraints are weakened
for the galaxy clustering case, due to the weaker scaledependence of the bias.
In the lensing case, the constraints remain instead similar in all
configurations.

We present an analysis of the crosscorrelation between the CMB and the
largescale structure (LSS) of the Universe in Unified Dark Matter (UDM) scalar
field cosmologies. We work out the predicted crosscorrelation function in UDM
models, which depends on the speed of sound of the unified component, and
compare it with observations from six galaxy catalogues (NVSS, HEAO, 2MASS, and
SDSS main galaxies, luminous red galaxies, and quasars). We sample the value of
the speed of sound and perform a likelihood analysis, finding that the UDM
model is as likely as the LambdaCDM, and is compatible with observations for a
range of values of c_\infinity (the value of the sound speed at late times) on
which structure formation depends. In particular, we obtain an upper bound of
c_\infinity^2 \leq 0.009 at 95% confidence level, meaning that the LambdaCDM
model, for which c_\infinity^2 = 0, is a good fit to the data, while the
posterior probability distribution peaks at the value c_\infinity^2=10^(4) .
Finally, we study the time dependence of the deviation from LambdaCDM via a
tomographic analysis using a mock redshift distribution and we find that the
largest deviation is for lowredshift sources, suggesting that future lowz
surveys will be best suited to constrain UDM models.

A large submm survey with Herschel will enable many exciting science
opportunities, especially in an era of widefield optical and radio surveys and
high resolution cosmic microwave background experiments. The HerschelSPIRE
Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250,
350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3
million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion
noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere
(3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which
have extensive multiwavelength coverage and are easily accessible from ALMA.
Two thirds of the of the sources are expected to be at z > 1, one third at z >
2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary
anisotropy studies with Planck and groundbased CMB data. (c) Find at least
1200 strongly lensed bright submm sources leading to a 2% test of general
relativity. (d) Identify 200 protocluster regions at z of 2 and perform an
unbiased study of the environmental dependence of star formation. (e) Perform
an unbiased survey for star formation and dust at high Galactic latitude and
make a census of debris disks and dust around AGB stars and white dwarfs.

We study the effects of modified theories of gravity on the cosmic microwave
background (CMB) anisotropies power spectrum, and in particular on its large
scales, where the integrated SachsWolfe (ISW) effect is important. Starting
with a general parametrisation, we then specialise to f(R) theories and
theories with Yukawatype interactions between dark matter particles. In these
models, the evolution of the metric potentials is altered, and the contribution
to the ISW effect can differ significantly from that in the standard model of
cosmology. We proceed to compare these predictions with observational data for
the CMB and the ISW, performing a full Monte Carlo Markov chain (MCMC)
analysis. In the case of f(R) theories, the result is an upper limit on the
lengthscale associated to the extra scalar degree of freedom characterising
these theories. With the addition of data from the Hubble diagram of Type Ia
supernovae, we obtain an upper limit on the lengthscale of the theory of B_0 <
0.4, or correspondingly \lambda_1 < 1900 Mpc/h at 95% c.l. improving previous
CMB constraints. For Yukawatype models we get a bound on the coupling 0.75 <
\beta_1 < 1.25 at the 95% c.l. We also discuss the implications of the assumed
priors on the estimation of modified gravity parameters, showing that a
marginally less conservative choice improves the f(R) constraints to \lambda_1
< 1400 Mpc/h, corresponding to B_0 < 0.2 at 95% c.l.

We study structure formation in the presence of primordial nonGaussianity of
the local type with parameters f_NL and g_NL. We show that the distribution of
darkmatter halos is naturally described by a multivariate bias scheme where
the halo overdensity depends not only on the underlying matter density
fluctuation delta, but also on the Gaussian part of the primordial
gravitational potential phi. This corresponds to a nonlocal bias scheme in
terms of delta only. We derive the coefficients of the bias expansion as a
function of the halo mass by applying the peakbackground split to common
parametrizations for the halo mass function in the nonGaussian scenario. We
then compute the halo power spectrum and halomatter cross spectrum in the
framework of Eulerian perturbation theory up to third order. Comparing our
results against Nbody simulations, we find that our model accurately describes
the numerical data for wavenumbers k < 0.10.3 h/Mpc depending on redshift and
halo mass. In our multivariate approach, perturbations in the halo counts trace
phi on large scales and this explains why the halo and matter power spectra
show different asymptotic trends for k > 0. This strongly scaledependent bias
originates from terms at leading order in our expansion. This is different from
what happens using the standard univariate local bias where the scaledependent
terms come from badly behaved higherorder corrections. On the other hand, our
biasing scheme reduces to the usual local bias on smaller scales where phi is
typically much smaller than the density perturbations. We finally discuss the
halo bispectrum in the context of multivariate biasing and show that, due to
its strong scale and shape dependence, it is a powerful tool for the detection
of primordial nonGaussianity from future galaxy surveys.

We present posterior likelihoods and Bayesian model selection analysis for
generalized cosmological models where the primordial perturbations include
correlated adiabatic and cold dark matter isocurvature components. We perform
nested sampling with flat and, for the first time, curved spatial geometries of
the Universe, using data from the cosmic microwave background (CMB)
anisotropies, the Union supernovae (SN) sample and a combined measurement of
the integrated SachsWolfe (ISW) effect. The CMB alone favors a 3% (positively
correlated) isocurvature contribution in both the flat and curved cases. The
nonadiabatic contribution to the observed CMB temperature variance is 0 <
alpha_T < 7% at 98% CL in the curved case. In the flat case, combining the CMB
with SN data artificially biases the result towards the pure adiabatic LCDM
concordance model, whereas in the curved case the favored level of
nonadiabaticity stays at 3% level with all combinations of data. However, the
ratio of Bayes factors, or Delta ln(evidence), is more than 5 points in favor
of the flat adiabatic LCDM model, which suggests that the inclusion of the 5
extra parameters of the curved isocurvature model is not supported by the
current data. The results are very sensitive to the second and third acoustic
peak regions in the CMB temperature angular power: therefore a careful
calibration of these data will be required before drawing decisive conclusions
on the nature of primordial perturbations. Finally, we point out that the odds
for the flat nonadiabatic model are 1:3 compared to the curved adiabatic
model. This may suggest that it is not much less motivated to extend the
concordance model with 4 isocurvature degrees of freedom than it is to study
the spatially curved adiabatic model.

I present to this conference our latest measurements of the integrated
SachsWolfe (ISW) effect. After a brief review of the reasons for which this
effect arises and of the technique to detect it by crosscorrelating the cosmic
microwave background (CMB) with the large scale structure of the Universe
(LSS), I describe the current state of the art measurement. This is obtained
from a combined analysis of six different galaxy datasets, and has a
significance level of ~ 4.5 sigma. I then describe the cosmological
implications, which show agreement with a flat LCDM model with Omega_m = 0.20
+0.19 0.11 at 95% confidence level. I finally show how these data can be used
to constrain modified gravity theories, focusing in particular on the
DvaliGabadazePorrati (DGP) model.

We present a global measurement of the integrated SachsWolfe (ISW) effect
obtained by crosscorrelating all relevant large scale galaxy data sets with
the cosmic microwave background radiation map provided by the Wilkinson
Microwave Anisotropy Probe. With these measurements, the overall ISW signal is
detected at the ~ 4.5 sigma level. We also examine the cosmological
implications of these measurements, particularly the dark energy equation of
state w, its sound speed, and the overall curvature of the Universe. The flat
LCDM model is a good fit to the data and, assuming this model, we find that the
ISW data constrain Omega_m = 0.20 +0.19 0.11 at the 95% confidence level. When
we combine our ISW results with the latest baryon oscillation and supernovae
measurements, we find that the result is still consistent with a flat LCDM
model with w = 1 out to redshifts z > 1.