
We compare reduced threepoint correlations $Q$ of matter, haloes (as proxies
for galaxies) and their cross correlations, measured in a total simulated
volume of $\sim100 \ (h^{1} \text{Gpc})^{3}$, to predictions from leading
order perturbation theory on a large range of scales in configuration space.
Predictions for haloes are based on the nonlocal bias model, employing linear
($b_1$) and nonlinear ($c_2$, $g_2$) bias parameters, which have been
constrained previously from the bispectrum in Fourier space. We also study
predictions from two other bias models, one local ($g_2=0$) and one in which
$c_2$ and $g_2$ are determined by $b_1$ via an approximately universal
relation. Overall, measurements and predictions agree when $Q$ is derived for
triangles with $(r_1r_2r_3)^{1/3} \gtrsim 60 h^{1}\text{Mpc}$, where $r_{13}$
are the sizes of the triangle legs. Predictions for $Q_{matter}$, based on the
linear power spectrum, show significant deviations from the measurements at the
BAO scale (given our small measurement errors), which strongly decrease when
adding a damping term or using the nonlinear power spectrum, as expected.
Predictions for $Q_{halo}$ agree best with measurements at large scales when
considering nonlocal contributions. The universal bias model works well for
haloes and might therefore be also useful for tightening constraints on $b_1$
from $Q$ in galaxy surveys. Such constraints are independent of the amplitude
of matter density fluctuation ($\sigma_8$) and hence break the degeneracy
between $b_1$ and $\sigma_8$, present in galaxy twopoint correlations.

Efficient estimators of Fourierspace statistics for large number of objects
rely on Fast Fourier Transforms (FFTs), which are affected by aliasing from
unresolved small scale modes due to the finite FFT grid. Aliasing takes the
form of a sum over images, each of them corresponding to the Fourier content
displaced by increasing multiples of the sampling frequency of the grid. These
spurious contributions limit the accuracy in the estimation of Fourierspace
statistics, and are typically ameliorated by simultaneously increasing grid
size and discarding highfrequency modes. This results in inefficient estimates
for e.g. the power spectrum when desired systematic biases are well under
percent level. We show that using interlaced grids removes odd images, which
include the dominant contribution to aliasing. In addition, we discuss the
choice of interpolation kernel used to define density perturbations on the FFT
grid and demonstrate that using higherorder interpolation kernels than the
standard Cloud in Cell algorithm results in significant reduction of the
remaining images. We show that combining fourthorder interpolation with
interlacing gives very accurate Fourier amplitudes and phases of density
perturbations. This results in power spectrum and bispectrum estimates that
have systematic biases below 0.01% all the way to the Nyquist frequency of the
grid, thus maximizing the use of unbiased Fourier coefficients for a given grid
size and greatly reducing systematics for applications to large cosmological
datasets.

Approximate methods to full Nbody simulations provide a fast and accurate
solution to the development of mock catalogues for the modeling of galaxy
clustering observables. In this paper we extend ICECOLA (Izard et al. 2016),
based on an optimized implementation of the approximate COLA method, to produce
weak lensing maps and halo catalogues in the light cone using an integrated and
self consistent approach. We show that despite the approximate dynamics, the
catalogues thus produced enable an accurate modeling of weak lensing
observables one decade beyond the characteristic scale where the growth becomes
nonlinear. In particular, we compare ICECOLA to the MICEGC Nbody simulation
for some fiducial cases representative of upcoming surveys and find that, for
sources at redshift $z=1$, their convergence power spectra agree to within one
percent up to high multipoles (i.e., of order $1000$). The corresponding shear
two point functions, $\xi_{+}$ and $\xi_{}$, yield similar accuracy down to
$2$ and $20$ arcmin respectively, while tangential shear around a $z=0.5$ lens
sample is accurate down to $4$ arcmin. We show that such accuracy is stable
against an increased angular resolution of the weak lensing maps. Hence, this
opens the possibility of using approximate methods for the joint modeling of
galaxy clustering and weak lensing observables and their covariance in ongoing
and future galaxy surveys.

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 investigate the cosmological implications of studying galaxy clustering
using a tomographic approach applied to the final BOSS DR12 galaxy sample,
including both auto and crosscorrelation functions between redshift shells.
We model the signal of the full shape of the angular correlation function,
$\omega(\theta)$, in redshift bins using stateoftheart modelling of
nonlinearities, bias and redshiftspace distortions. We present results on the
redshift evolution of the linear bias of BOSS galaxies, which cannot be
obtained with traditional methods for galaxyclustering analysis. We also
obtain constraints on cosmological parameters, combining this tomographic
analysis with measurements of the cosmic microwave background (CMB) and type Ia
supernova (SNIa). We explore a number of cosmological models, including the
standard $\Lambda$CDM model and its most interesting extensions, such as
deviations from $w_\rm{DE} = 1$, nonminimal neutrino masses, spatial
curvature and deviations from general relativity using the growthindex
$\gamma$ parametrisation. These results are, in general, comparable to the most
precise presentday constraints on cosmological parameters, and show very good
agreement with the standard model. In particular, combining CMB,
$\omega(\theta)$ and SNIa, we find a value of $w_\rm{DE}$ consistent with $1$
to a precision better than 5\% when it is assumed to be constant in time, and
better than 6\% when we also allow for a spatiallycurved Universe.

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 extract cosmological information from the anisotropic power spectrum
measurements from the recently completed Baryon Oscillation Spectroscopic
Survey (BOSS), extending the concept of clustering wedges to Fourier space.
Making use of new FFTbased estimators, we measure the power spectrum
clustering wedges of the BOSS sample by filtering out the information of
Legendre multipoles l > 4. Our modelling of these measurements is based on
novel approaches to describe nonlinear evolution, bias, and redshiftspace
distortions, which we test using synthetic catalogues based on largevolume
Nbody simulations. We are able to include smaller scales than in previous
analyses, resulting in tighter cosmological constraints. Using three
overlapping redshift bins, we measure the angular diameter distance, the Hubble
parameter, and the cosmic growth rate, and explore the cosmological
implications of our full shape clustering measurements in combination with CMB
and SN Ia data. Assuming a {\Lambda}CDM cosmology, we constrain the matter
density to {\Omega}_m = 0.311 0.010 +0.009 and the Hubble parameter to H_0 =
67.6 0.6 +0.7 km s^1 Mpc^1, at a confidence level (CL) of 68 per cent. We
also allow for nonstandard dark energy models and modifications of the growth
rate, finding good agreement with the {\Lambda}CDM paradigm. For example, we
constrain the equationofstate parameter to w = 1.019 0.039 +0.048. This
paper is part of a set that analyses the final galaxy clustering dataset from
BOSS. The measurements and likelihoods presented here are combined with others
in Alam et al. 2016 to produce the final cosmological constraints from BOSS.

We explore the cosmological implications of anisotropic clustering
measurements in configuration space of the final galaxy samples from Data
Release 12 of the SDSSIII Baryon Oscillation Spectroscopic Survey. We
implement a new detailed modelling of the effects of nonlinearities, galaxy
bias and redshiftspace distortions that can be used to extract unbiased
cosmological information from our measurements for scales $s \gtrsim
20\,h^{1}{\rm Mpc}$. We combined the galaxy clustering information from BOSS
with the latest cosmic microwave background (CMB) observations and Type Ia
supernovae samples and found no significant evidence for a deviation from the
$\Lambda$CDM cosmological model. In particular, these data sets can constrain
the dark energy equation of state parameter to $w_{\rm DE}=0.996\pm0.042$ when
assumed timeindependent, the curvature of the Universe to
$\Omega_{k}=0.0007\pm 0.0030$ and the sum of the neutrino masses to $\sum
m_{\nu} < 0.25\,{\rm eV}$ at 95 per cent CL. We explore the constraints on the
growth rate of cosmic structures assuming $f(z)=\Omega_{\rm m}(z)^\gamma$ and
obtain $\gamma = 0.609\pm 0.079$, in good agreement with the predictions of
general relativity of $\gamma=0.55$. We compress the information of our
clustering measurements into constraints on the parameter combinations $D_{\rm
V}(z)/r_{\rm d}$, $F_{\rm AP}(z)$ and $f\sigma_8(z)$ at the effective redshifts
of $z=0.38$, $0.51$ and $0.61$ with their respective covariance matrices and
find good agreement with the predictions for these parameters obtained from the
bestfitting $\Lambda$CDM model to the CMB data from the Planck satellite. This
paper is part of a set that analyses the final galaxy clustering dataset from
BOSS. The measurements and likelihoods presented here are combined with others
in Alam et al. (2016) to produce the final cosmological constraints from BOSS.

We present a new method to measure the redshiftdependent galaxy bias by
combining information from the galaxy density field and the weak lensing field.
This method is based on Amara et al. (2012), where they use the galaxy density
field to construct a biasweighted convergence field kg. The main difference
between Amara et al. (2012) and our new implementation is that here we present
another way to measure galaxy bias using tomography instead of bias
parameterizations. The correlation between kg and the true lensing field k
allows us to measure galaxy bias using different zerolag correlations, such as
<kgk>/<kk> or <kgkg>/<kgk>. Our method measures the linear bias factor on
linear scales under the assumption of no stochasticity between galaxies and
matter. We use the MICE simulation to measure the linear galaxy bias for a
fluxlimited sample (i < 22.5) in tomographic redshift bins using this method.
This paper is the first that studies the accuracy and systematic uncertainties
associated with the implementation of the method, and the regime where it is
consistent with the linear galaxy bias defined by projected 2point correlation
functions (2PCF). We find that our method is consistent with linear bias at the
percent level for scales larger than 30 arcmin, while nonlinearities appear at
smaller scales. This measurement is a good complement to other measurements of
bias, since it does not depend strongly on sigma8 as the 2PCF measurements. We
apply this method to the Dark Energy Survey Science Verification data in a
followup paper.

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.

Next generation galaxy surveys demand the development of massive ensembles of
galaxy mocks to model the observables and their covariances, what is
computationally prohibitive using $N$body simulations. COLA is a novel method
designed to make this feasible by following an approximate dynamics but with up
to 3 orders of magnitude speedups when compared to an exact $N$body. In this
paper we investigate the optimization of the code parameters in the compromise
between computational cost and recovered accuracy in observables such as
twopoint clustering and halo abundance. We benchmark those observables with a
stateoftheart $N$body run, the MICE Grand Challenge simulation (MICEGC).
We find that using 40 time steps linearly spaced since $z_i \sim 20$, and a
force mesh resolution three times finer than that of the number of particles,
yields a matter power spectrum within $1\%$ for $k \lesssim 1\,h {\rm
Mpc}^{1}$ and a halo mass function within $5\%$ of those in the $N$body. In
turn the halo bias is accurate within $2\%$ for $k \lesssim 0.7\,h {\rm
Mpc}^{1}$ whereas, in redshift space, the halo monopole and quadrupole are
within $4\%$ for $k \lesssim 0.4\,h {\rm Mpc}^{1}$. These results hold for a
broad range in redshift ($0 < z < 1$) and for all halo mass bins investigated
($M > 10^{12.5} \, h^{1} \, {\rm M_{\odot}}$). To bring accuracy in clustering
to one percent level we study various methods that recalibrate halo masses
and/or velocities. We thus propose an optimized choice of COLA code parameters
as a powerful tool to optimally exploit future galaxy surveys.

We present a comparison of major methodologies of fast generating mock halo
or galaxy catalogues. The comparison is done for twopoint and the threepoint
clustering statistics. The reference catalogues are drawn from the BigMultiDark
Nbody simulation. Both friendoffriends (including distinct halos only) and
spherical overdensity (including distinct halos and subhalos) catalogs have
been used with the typical number density of a largevolume galaxy surveys. We
demonstrate that a proper biasing model is essential for reproducing the power
spectrum at quasilinear and even smaller scales. With respect to various
clustering statistics a methodology based on perturbation theory and a
realistic biasing model leads to very good agreement with Nbody simulations.
However, for the quadrupole of the correlation function or the power spectrum,
only the method based on semiNbody simulation could reach high accuracy (1%
level) at small scales, i.e., r<25 Mpc/h or k>0.15 h/Mpc. Full Nbody solutions
will remain indispensable to produce reference catalogues. Nevertheless, we
have demonstrated that the far more efficient approximate solvers can reach a
few percent accuracy in terms of clustering statistics at the scales
interesting for the largescale structure analysis after calibration with a few
reference Nbody calculations. This makes them useful for massive production
aimed at covariance studies, to scan large parameter spaces, and to estimate
uncertainties in data analysis techniques, such as baryon acoustic oscillation
reconstruction, redshift distortion measurements, etc.

Measurements of the linear growth factor $D$ at different redshifts $z$ are
key to distinguish among cosmological models. One can estimate the derivative
$dD(z)/d\ln(1+z)$ from redshift space measurements of the 3D anisotropic galaxy
twopoint correlation $\xi(z)$, but the degeneracy of its transverse (or
projected) component with galaxy bias $b$, i.e. $\xi_{\perp}(z) \propto\ D^2(z)
b^2(z)$, introduces large errors in the growth measurement. Here we present a
comparison between two methods which break this degeneracy by combining second
and thirdorder statistics. One uses the shape of the reduced threepoint
correlation and the other a combination of thirdorder one and twopoint
cumulants. These methods use the fact that, for Gaussian initial conditions and
scales larger than $20$ $h^{1}$Mpc, the reduced thirdorder matter
correlations are independent of redshift (and therefore of the growth factor)
while the thirdorder galaxy correlations depend on $b$. We use matter and halo
catalogs from the MICEGC simulation to test how well we can recover $b(z)$ and
therefore $D(z)$ with these methods in 3D real space. We also present a new
approach, which enables us to measure $D$ directly from the redshift evolution
of second and thirdorder galaxy correlations without the need of modelling
matter correlations. For haloes with masses lower than $10^{14}$
$h^{1}$M$_\odot$, we find $10%$ deviations between the different estimates of
$D$, which are comparable to current observational errors. At higher masses we
find larger differences that can probably be attributed to the breakdown of the
bias model and nonPoissonian shot noise.

Several papers have recently highlighted the possibility of measuring
redshift space distortions from angular autocorrelations of galaxies in
photometric redshift bins. In this work we extend this idea to include as
observables the crosscorrelations between redshift bins, as an additional way
of measuring radial information. We show that this extra information allows to
reduce the recovered error in the growth rate index \gamma by a factor of ~2.
Although the final error in \gamma depends on the bias and the mean photometric
accuracy of the galaxy sample, the improvement from adding crosscorrelations
is robust in different settings. Another factor of 23 improvement in the
determination of \gamma can be achieved by considering two galaxy populations
over the same photometric sky area but with different biases. This additional
gain is shown to be much larger than the one from the same populations when
observed over different areas of the sky (with twice the combined area). The
total improvement of ~5 implies that a photometric survey such as the Dark
Energy Survey should be able to recover \gamma at the 510% from the angular
clustering in linear scales of two different tracers. It can also constrain the
evolution of f(z)x\sigma_8(z) in few bins beyond z~0.80.9 at the 1015% level
perbin, compatible with recent constrains from lowerz spectroscopic surveys.
We also show how further improvement can be achieved by reducing the
photometric redshift error.

We present scaledependent measurements of the normalised growth rate of
structure $f\sigma_{8}(k, z=0)$ using only the peculiar motions of galaxies. We
use data from the 6degree Field Galaxy Survey velocity sample (6dFGSv)
together with a newlycompiled sample of lowredshift $(z < 0.07)$ type Ia
supernovae. We constrain the growth rate in a series of $\Delta k \sim 0.03
h{\rm Mpc^{1}}$ bins to $\sim35\%$ precision, including a measurement on
scales $>300 h^{1}{\rm Mpc}$, which represents one of the largestscale growth
rate measurement to date. We find no evidence for a scale dependence in the
growth rate, or any statistically significant variation from the growth rate as
predicted by the {\it Planck} cosmology. Bringing all the scales together, we
determine the normalised growth rate at $z=0$ to $\sim15\%$ in a manner {\it
independent} of galaxy bias and in excellent agreement with the constraint from
the measurements of redshiftspace distortions from 6dFGS. We pay particular
attention to systematic errors. We point out that the intrinsic scatter present
in FundamentalPlane and TullyFisher relations is only Gaussian in logarithmic
distance units; wrongly assuming it is Gaussian in linear (velocity) units can
bias cosmological constraints. We also analytically marginalise over zeropoint
errors in distance indicators, validate the accuracy of all our constraints
using numerical simulations, and demonstrate how to combine different
(correlated) velocity surveys using a matrix `hyperparameter' analysis.
Current and forthcoming peculiar velocity surveys will allow us to understand
in detail the growth of structure in the lowredshift universe, providing
strong constraints on the nature of dark energy.

Future largescale structure surveys of the Universe will aim to constrain
the cosmological model and the true nature of dark energy with unprecedented
accuracy. In order for these surveys to achieve their designed goals, they will
require predictions for the nonlinear matter power spectrum to subpercent
accuracy. Through the use of a large ensemble of cosmological Nbody
simulations, we demonstrate that if we do not understand the uncertainties
associated with simulating structure formation, i.e. knowledge of the `true'
simulation parameters, and simply seek to marginalize over them, then the
constraining power of such future surveys can be significantly reduced.
However, for the parameters {n_s, h, Om_b, Om_m}, this effect can be largely
mitigated by adding the information from a CMB experiment, like Planck. In
contrast, for the amplitude of fluctuations sigma8 and the timeevolving
equation of state of dark energy {w_0, w_a}, the mitigation is mild. On
marginalizing over the simulation parameters, we find that the darkenergy
figure of merit can be degraded by ~2. This is likely an optimistic assessment,
since we do not take into account other important simulation parameters. A
caveat is our assumption that the Hessian of the likelihood function does not
vary significantly when moving from our adopted to the 'true' simulation
parameter set. This paper therefore provides strong motivation for rigorous
convergence testing of Nbody codes to meet the future challenges of precision
cosmology.

We study how to recover the full 3D clustering information of P(\vec{k},z),
including redshift space distortions (RSD), from 2D tomography using the
angular auto and cross spectra of different redshift bins C_\ell(z,z'). We
focus on quasilinear scales where the minimum scale \lambda_{min} or
corresponding maximum wavenumber k_{max}= 2\pi/\lambda_{min} is targeted to be
between k_{max}={0.050.2} h/Mpc. For spectroscopic surveys, we find that we
can recover the full 3D clustering information when the redshift bin width
\Delta z used in the 2D tomography is similar to the targeted minimum scale,
i.e. \Delta z ~ {0.60.8} \lambda_{min} H(z)/c which corresponds to \Delta z ~
0.010.05 for z<1. This value of \Delta z is optimal in the sense that larger
values of \Delta z lose information, while smaller values violate our minimum
scale requirement. For a narrowband photometric survey, with photoz error
\sigma_z=0.004, we find almost identical results to the spectroscopic survey
because the photoz error is smaller than the optimal bin width \sigma_z<\Delta
z. For a typical broadband photometric survey with \sigma_z=0.1, we have that
\sigma_z>\Delta z and most radial information is intrinsically lost. The
remaining information can be recovered from the 2D tomography if we use \Delta
z ~ 2\sigma_z. While 3D and 2D analysis are shown here to be equivalent, the
advantage of using angular positions and redshifts is that we do not need a
fiducial cosmology to convert to 3D coordinates. This avoids assumptions and
marginalization over the fiducial model. In addition, it becomes straight
forward to combine RSD, clustering and weak lensing in 2D space.

We put forward and test a simple description of multipoint propagators (MP),
which serve as buildingblocks to calculate the nonlinear matter power
spectrum. On large scales these propagators reduce to the wellknown kernels in
standard perturbation theory, while at smaller scales they are suppresed due to
nonlinear couplings. Through extensive testing with numerical simulations we
find that this decay is characterized by the same damping scale for both two
and threepoint propagators. In turn this transition can be well modeled with
resummation results that exponentiate oneloop computations. For the first
time, we measure the four components of the nonlinear (twopoint) propagator
using dedicated simulations started from two independent random Gaussian fields
for positions and velocities, verifying in detail the fundamentals of
propagator resummation.
We use these results to develop an implementation of the MPexpansion for the
nonlinear power spectrum that only requires seconds to evaluate at BAO scales.
To test it we construct six suites of large numerical simulations with
different cosmologies. From these and
LasDamas runs we show that the nonlinear power spectrum can be described at
the ~ 2% level at BAO scales for redshifts in the range [02.5]. We make a
public release of the MPTbreeze code with the hope that it can be useful to the
community.

We present a new scheme for the general computation of cosmic propagators
that allow to interpolate between standard perturbative results at lowk and
their expected largek resummed behavior. This scheme is applicable to any
multipoint propagator and allows the matching of perturbative lowk
calculations to any number of loops to their largek behavior, and can
potentially be applied in case of nonstandard cosmological scenarios such as
those with nonGaussian initial conditions. The validity of our proposal is
checked against previous prescriptions and measurements in numerical
simulations showing a remarkably good agreement. Such a generic prescription
for multipoint propagators provides the necessary building blocks for the
computation of polyspectra in the context of the socalled Gammaexpansion
introduced by Bernardeau et al. (2008). As a concrete application we present a
consistent calculation of the matter bispectrum at oneloop order.

We present measurements of the bispectrum of dark matter halos in numerical
simulations with nonGaussian initial conditions of the local type. We show, in
the first place, that the overall effect of primordial nonGaussianity on the
halo bispectrum is larger than on the halo power spectrum when all measurable
configurations are taken into account. We then compare our measurements with a
treelevel perturbative prediction finding good agreement at large scale when
the constant Gaussian bias parameter, both linear and quadratic, and their
constant nonGaussian corrections are fitted for. The bestfit values of the
Gaussian bias factors and their nonGaussian, scaleindependent corrections are
in qualitative agreement with the peakbackground split expectations. In
particular, we show that the effect of nonGaussian initial conditions on
squeezed configurations is fairly large (up to 30% for f_NL=100 at redshift
z=0.5) and results from contributions of similar amplitude induced by the
initial matter bispectrum, scaledependent bias corrections as well as from
nonlinear matter bispectrum corrections. We show, in addition, that effects at
second order in f_NL are irrelevant for the range of values allowed by CMB and
galaxy power spectrum measurements, at least on the scales probed by our
simulations. Finally, we present a Fisher matrix analysis to assess the
possibility of constraining primordial nonGaussianity with future measurements
of the galaxy bispectrum. We find that a survey with a volume of about 10 cubic
Gpc at mean redshift z ~ 1 could provide an error on f_NL of the order of a
few. This shows the relevance of a joint analysis of galaxy power spectrum and
bispectrum in future redshift surveys.

The possibility of measuring redshift space (RSD) distortions using
photometric data have been recently highlighted. This effect complements and
significantly alters the detectability of baryon acoustic oscillations (BAO) in
photometric surveys. In this paper we present measurements of the angular
correlation function of luminous red galaxies (LRGs) in the photometric catalog
of the final data release (DR7) of the Sloan Digital Sky Survey II (SDSS). The
sample compromise ~ 1.5 x 10^6 LRGs distributed in 0.45 < z < 0.65, with a
characteristic photometric error of ~ 0.05. Our measured correlation centered
at z=0.55 is in very good agreement with predictions from standard LCDM in a
broad range of angular scales, $0.5^\circ < \theta < 6^\circ$. We find that the
growth of structure can indeed be robustly measured, with errors matching
expectations. The velocity growth rate is recovered as $f \sigma_8 = 0.53 \pm
0.42$ when no prior is imposed on the growth factor and the background geometry
follows a LCDM model with WMAP7+SNIa priors. This is compatible with the
corresponding General Relativity (GR) prediction $f \sigma_8 = 0.45$ for our
fiducial cosmology. If we adopt a parametrization such that $f=\Omega
^\gamma_m(z)$, with $\gamma \approx 0.55$ in GR, and combine our $f\sigma_8$
measurement with the corresponding ones from spectroscopic LRGs at lower
redshifts we obtain $\gamma=0.54 \pm 0.17$. In addition we find evidence for
the presence of the baryon acoustic feature matching the amplitude, location
and shape of LCDM predictions. The photometric BAO feature is detected with 98
% confidence level at z=0.55.

Near future cosmology will see the advent of wide area photometric galaxy
surveys, like the Dark Energy Survey (DES), that extent to high redshifts (z ~
1  2) but with poor radial distance resolution. In such cases splitting the
data into redshift bins and using the angular correlation function $w(\theta)$,
or the $C_{\ell}$ power spectrum, will become the standard approach to extract
cosmological information or to study the nature of dark energy through the
Baryon Acoustic Oscillations (BAO) probe. In this work we present a detailed
model for $w(\theta)$ at large scales as a function of redshift and bin width,
including all relevant effects, namely nonlinear gravitational clustering,
bias, redshift space distortions and photoz uncertainties. We also present a
model for the full covariance matrix characterizing the angular correlation
measurements, that takes into account the same effects as for $w(\theta)$ and
also the possibility of a shotnoise component and partial sky coverage.
Provided with a large volume Nbody simulation from the MICE collaboration we
built several ensembles of mock redshift bins with a sky coverage and depth
typical of forthcoming photometric surveys. The model for the angular
correlation and the one for the covariance matrix agree remarkably well with
the mock measurements in all configurations. The prospects for a full shape
analysis of $w(\theta)$ at BAO scales in forthcoming photometric surveys such
as DES are thus very encouraging.

Models of galaxy and halo clustering commonly assume that the tracers can be
treated as a continuous field locally biased with respect to the underlying
mass distribution. In the peak model pioneered by BBKS, one considers instead
density maxima of the initial, Gaussian mass density field as an approximation
to the formation site of virialized objects. In this paper, the peak model is
extended in two ways to improve its predictive accuracy. Firstly, we derive the
twopoint correlation function of initial density peaks up to second order and
demonstrate that a peakbackground split approach can be applied to obtain the
kindependent and kdependent peak bias factors at all orders. Secondly, we
explore the gravitational evolution of the peak correlation function within the
Zel'dovich approximation. We show that the local (Lagrangian) bias approach
emerges as a special case of the peak model, in which all bias parameters are
scaleindependent and there is no statistical velocity bias. We apply our
formulae to study how the Lagrangian peak biasing, the diffusion due to large
scale flows and the modecoupling due to nonlocal interactions affect the scale
dependence of bias from small separations up to the baryon acoustic oscillation
(BAO) scale. For 2sigma density peaks collapsing at z=0.3, our model predicts
a ~ 5% residual scaledependent bias around the acoustic scale that arises
mostly from firstorder Lagrangian peak biasing (as opposed to secondorder
gravity modecoupling). We also search for a scale dependence of bias in the
large scale autocorrelation of massive halos extracted from a very large
Nbody simulation provided by the MICE collaboration. For halos with mass
M>10^{14}Msun/h, our measurements demonstrate a scaledependent bias across the
BAO feature which is very well reproduced by a prediction based on the peak
model.

We show here how Renormalized Perturbation Theory (RPT) calculations applied
to the quasilinear growth of the largescale structure can be carried on in
presence of primordial nonGaussian (PNG) initial conditions. It is explicitly
demonstrated that the series reordering scheme proposed in Bernardeau, Crocce
and Scoccimarro (2008) is preserved for nonGaussian initial conditions. This
scheme applies to the power spectrum and higher order spectra and is based on a
reorganization of the contributing terms into sum of products of multipoint
propagators. In case of PNG new contributing terms appear, the importance of
which is discussed in the context of current PNG models. The properties of the
building blocks of such resummation schemes, the multipoint propagators, are
then investigated. It is first remarked that their expressions are left
unchanged at oneloop order irrespectively of statistical properties of the
initial field. We furthermore show that the highmomemtum limit of each of
these propagators can be explicitly computed even for arbitrary initial
conditions. They are found to be damped by an exponential cutoff whose
expression is directly related to the moment generating function of the
onedimensional displacement field. This extends what had been established for
multipoint propagators for Gaussian initial conditions. Numerical forms of the
cutoff are shown for the socalled local model of PNG.