
We present a method to reconstruct the initial conditions of the universe
using observed galaxy positions and luminosities under the assumption that the
luminosities can be calibrated with weak lensing to give the mean halo mass.
Our method relies on following the gradients of forward model and since the
standard way to identify halos is nondifferentiable and results in a discrete
sample of objects, we propose a framework to model the halo position and mass
field starting from the nonlinear matter field using Neural Networks. We
evaluate the performance of our model with multiple metrics. Our model is more
than $95\%$ correlated with the halomass fields up to $k\sim 0.7 {\rm h/Mpc}$
and significantly reduces the stochasticity over the Poisson shot noise. We
develop a data likelihood model that takes our modeling error and intrinsic
scatter in the halo masslight relation into account and show that a displaced
lognormal model is a good approximation to it. We optimize over the
corresponding loss function to reconstruct the initial density field and
develop an annealing procedure to speed up and improve the convergence. We
apply the method to halo number densities of $\bar{n} = 2.5\times 10^{4}
10^{3}({\rm h/Mpc})^3$, typical of current and future redshift surveys, and
recover a Gaussian initial density field, mapping all the higher order
information in the data into the power spectrum. We show that our
reconstruction improves over the standard reconstruction. For baryonic acoustic
oscillations (BAO) the gains are relatively modest because BAO is dominated by
large scales where standard reconstruction suffices. We improve upon it by
$\sim 1520\%$ in terms of error on BAO peak as estimated by Fisher analysis at
$z=0$. We expect larger gains will be achieved when applying this method to the
broadband linear power spectrum reconstruction on smaller scales.

The large scale structure (LSS) of the universe is generated by the linear
density gaussian modes, which are evolved into the observed nonlinear LSS. The
posterior surface of the modes is convex in the linear regime, leading to a
unique global maximum (MAP), but this is no longer guaranteed in the nonlinear
regime. In this paper we investigate the nature of posterior surface using the
recently developed MAP reconstruction method, with a simplified but realistic
Nbody simulation as the forward model. The reconstruction method uses
optimization with analytic gradients from backpropagation through the
simulation. For low noise cases we recover the initial conditions well into the
nonlinear regime ($k\sim 1$ h/Mpc) nearly perfectly. We show that the large
scale modes can be recovered more precisely than the linear expectation, which
we argue is a consequence of nonlinear mode coupling. For noise levels
achievable with current and planned LSS surveys the reconstruction cannot
recover very small scales due to noise. We see some evidence of nonconvexity,
specially for smaller scales where the noninjective nature of the mappings:
several very different initial conditions leading to the same near perfect
final data reconstruction. We investigate the nature of these phenomena further
using a 1d toy gravity model, where many well separated local maximas are
found to have identical data likelihood but differ in the prior. We also show
that in 1d the prior favors some solutions over the true solution, though no
clear evidence of these in 3d. Our main conclusion is that on very small
scales and for a very low noise the posterior surface is multimodal and the
global maximum may be unreachable with standard methods, while for realistic
noise levels in the context of the current and next generation LSS surveys MAP
optimization method is likely to be nearly optimal.

Fast Nbody PM simulations with a small number of time steps such as FastPM
or COLA have been remarkably successful in modeling the galaxy statistics, but
their lack of small scale force resolution and long time steps cannot give
accurate halo matter profiles or matter power spectrum. High resolution Nbody
simulations can improve on this, but lack baryonic effects, which can only be
properly included in hydro simulations. Here we present a scheme to calibrate
the fast simulations to mimic the precision of the hydrodynamic simulations or
high resolution Nbody simulations. The scheme is based on a gradient descent
of either effective gravitational potential, which mimics the short range
force, or of effective enthalpy, which mimics gas hydrodynamics and feedback.
The scheme is fast and differentiable, and can be incorporated as a
postprocessing step into any simulation. It gives very good results for the
matter power spectrum for several of the baryonic feedback and dark matter
simulations, and also gives improved dark matter halo profiles. The scheme is
even able to find the large subhalos, and increase the correlation coefficient
between the fast simulations and the high resolution Nbody or hydro
simulations. It can also be used to add baryonic effects to the high resolution
Nbody simulations. While the method has free parameters that can be calibrated
on various simulations, they can also be viewed as astrophysical nuisance
parameters describing baryonic effects that can be marginalized over during the
data analysis. In this view these parameters can be viewed as an efficient
parametrization of baryonic effects.

We develop a new method to constraint primordial nonGaussianities of the
local kind using unclustered tracers of the Large Scale Structure. We show that
in the limit of low noise, zero bias tracers yield large improvement over
standard methods, mostly due to vanishing sampling variance. We propose a
simple technique to construct such a tracer, using environmental information
obtained from the original sample, and validate our method with Nbody
simulations. Our results indicate that
$\sigma_{f_{\mathrm{NL}}^{\mathrm{loc}}}\simeq1$ can be reached using only
information on a single tracer of sufficiently high number density.

We present measurements of $E_G$, a probe of gravity from largescale
structure, using BOSS LOWZ and CMASS spectroscopic samples, with lensing
measurements from SDSS (galaxy lensing) and Planck (CMB lensing). Using SDSS
lensing and the BOSS LOWZ sample, we measure $\langle
E_G\rangle=0.37^{+0.036}_{0.032}$ (stat) $\pm 0.026$ (systematic), consistent
with the predicted value from the Planck $\Lambda$CDM model, $E_G=0.46$,
considering the impact of systematic uncertainties. Using CMB lensing, we
measure $\langle E_G\rangle=0.43^{+0.068}_{0.073}$ (stat) for LOWZ
(statistically consistent with galaxy lensing and Planck predictions) and
$\langle E_G\rangle=0.39^{+0.05}_{0.05}$ (stat) for the CMASS sample,
consistent with the Planck prediction of $E_G=0.40$ given the higher redshift
of the sample. We also study the redshift evolution of $E_G$ by splitting the
LOWZ sample into two samples based on redshift, with results being consistent
with model predictions at $2.5\sigma$ (stat) or better. We correct for the
effects of nonlinear physics and different window functions for clustering and
lensing using analytical model and simulations, and demonstrate that these
corrections are effective to $\sim 12$\%, well within our statistical
uncertainties. For the case of SDSS lensing there is an additional $\sim5\%$
combined systematic uncertainty from shear calibration and photometric redshift
uncertainties.

One of the main unsolved problems of cosmology is how to maximize the
extraction of information from nonlinear data. If the data are nonlinear the
usual approach is to employ a sequence of statistics (Npoint statistics,
counting statistics of clusters, density peaks or voids etc.), along with the
corresponding covariance matrices. However, this approach is computationally
prohibitive and has not been shown to be exhaustive in terms of information
content. Here we instead develop a Bayesian approach, expanding the likelihood
around the maximum posterior of linear modes, which we solve for using
optimization methods. By integrating out the modes using perturbative expansion
of the likelihood we construct an initial power spectrum estimator, which for a
fixed forward model contains all the cosmological information if the initial
modes are gaussian distributed. We develop a method to construct the window and
covariance matrix such that the estimator is explicitly unbiased and nearly
optimal. We then generalize the method to include the forward model parameters,
including cosmological and nuisance parameters, and primordial nongaussianity.
We apply the method in the simplified context of nonlinear structure formation,
using either simplified 2LPT dynamics or Nbody simulations as the nonlinear
mapping between linear and nonlinear density, and 2LPT dynamics in the
optimization steps used to reconstruct the initial density modes. We
demonstrate that the method gives an unbiased estimator of the initial power
spectrum, providing among other a near optimal reconstruction of linear
baryonic acoustic oscillations.

As a major source of cosmological information, galaxy clustering is
susceptible to longwavelength density and tidal fluctuations. These long modes
modulate the growth and expansion rate of local structures, shifting them in
both amplitude and scale. These effects are often named the growth and dilation
effects, respectively. In particular the dilation shifts the baryon acoustic
oscillation (BAO) peak and breaks the assumption of the AlcockPaczynski (AP)
test. This cannot be removed with reconstruction techniques because the effect
originates from long modes outside the survey. In redshift space, the long
modes generate a largescale radial peculiar velocity that affects the
redshiftspace distortion (RSD) signal. We compute the redshiftspace response
functions of the galaxy power spectrum to long density and tidal modes at
leading order in perturbation theory, including both the growth and dilation
terms. We validate these response functions against measurements from simulated
galaxy mock catalogs. As one application, long density and tidal modes beyond
the scale of a survey correlate various observables leading to an excess error
known as the supersample covariance, and thus weaken their constraining power.
We quantify the supersample effect on BAO, AP, and RSD measurements, and study
its impact on current and future surveys.

We introduce FastPM, a highlyscalable approximated particle mesh Nbody
solver, which implements the particle mesh (PM) scheme enforcing correct linear
displacement (1LPT) evolution via modified kick and drift factors. Employing a
2dimensional domain decomposing scheme, FastPM scales extremely well with a
very large number of CPUs. In contrast to COmovingLAgrangian (COLA) approach,
we do not require to split the force or track separately the 2LPT solution,
reducing the code complexity and memory requirements. We compare FastPM with
different number of steps ($N_s$) and force resolution factor ($B$) against 3
benchmarks: halo mass function from Friends of Friends halo finder, halo and
dark matter power spectrum, and cross correlation coefficient (or
stochasticity), relative to a high resolution TreePM simulation. We show that
the modified time stepping scheme reduces the halo stochasticity when compared
to COLA with the same number of steps and force resolution. While increasing
$N_s$ and $B$ improves the transfer function and cross correlation coefficient,
for many applications FastPM achieves sufficient accuracy at low $N_s$ and $B$.
For example, $N_s=10$ and $B=2$ simulation provides a substantial saving (a
factor of 10) of computing time relative to $N_s=40$, $B=3$ simulation, yet the
halo benchmarks are very similar at $z=0$. We find that for abundance matched
halos the stochasticity remains low even for $N_s=5$. FastPM compares well
against less expensive schemes, being only 7 (4) times more expensive than 2LPT
initial condition generator for $N_s=10$ ($N_s=5$). Some of the applications
where FastPM can be useful are generating a large number of mocks, producing
nonlinear statistics where one varies a large number of nuisance or
cosmological parameters, or serving as part of an initial conditions solver.

We analyse the power spectrum of the Baryon Oscillation Spectroscopic Survey
(BOSS), Data Release 12 (DR12) to constrain the relative velocity effect, which
represents a potential systematic for measurements of the Baryon Acoustic
Oscillation (BAO) scale. The relative velocity effect is sourced by the
different evolution of baryon and cold dark matter perturbations before
decoupling. Our power spectrum model includes all $1$loop redshiftspace terms
corresponding to $v_{\rm bc}$ parameterised by the bias parameter $b_{v^2}$. We
also include the linear terms proportional to the relative density,
$\delta_{\rm bc}$, and relative velocity dispersion, $\theta_{\rm bc}$, which
we parameterise with the bias parameters $b^{\rm bc}_{\delta}$ and $b^{\rm
bc}_{\theta}$. Our data does not support a detection of the relative velocity
effect in any of these parameters. Combining the low and high redshift bins of
BOSS, we find limits of $b_{v^2} = 0.012 \pm 0.015\;(\pm 0.031)$, $b^{\rm
bc}_{\delta} = 1.0 \pm 2.5\;(\pm 6.2)$ and $b^{\rm bc}_{\theta} = 114 \pm
55\;(\pm 175)$ with $68\%$ ($95\%$) confidence levels. These constraints
restrict the potential systematic shift in $D_A(z)$, $H(z)$ and $f\sigma_8$,
due to the relative velocity, to $1\%$, $0.8\%$ and $2\%$, respectively. Given
the current uncertainties on the BAO measurements of BOSS these shifts
correspond to $0.53\sigma$, $0.5\sigma$ and $0.22\sigma$ for $D_A(z)$, $H(z)$
and $f\sigma_8$, respectively.

We study the covariance properties of real space correlation function
estimators  primarily galaxyshear correlations, or galaxygalaxy lensing 
using SDSS data for both shear catalogs and lenses (specifically the BOSS LOWZ
sample). Using mock catalogs of lenses and sources, we disentangle the various
contributions to the covariance matrix and compare them with a simple
analytical model. We show that not subtracting the lensing measurement around
random points from the measurement around the lens sample is equivalent to
performing the measurement using the lens density field instead of the lens
overdensity field. While the measurement using the lens density field is
unbiased (in the absence of systematics), its error is significantly larger due
to an additional term in the covariance. Therefore, this subtraction should be
performed regardless of its beneficial effects on systematics. Comparing the
error estimates from data and mocks for estimators that involve the
overdensity, we find that the errors are dominated by the shape noise and lens
clustering, that empirically estimated covariances (jackknife and standard
deviation across mocks) are consistent with theoretical estimates, and that
both the connected parts of the 4point function and the supersample
covariance can be neglected for the current levels of noise. While the
tradeoff between different terms in the covariance depends on the survey
configuration (area, source number density), the diagnostics that we use in
this work should be useful for future works to test their
empiricallydetermined covariances.

We present a new model for the redshiftspace power spectrum of galaxies and
demonstrate its accuracy in modeling the monopole, quadrupole, and hexadecapole
of the galaxy density field down to scales of $k = 0.4 \ h\mathrm{Mpc}^{1}$.
The model describes the clustering of galaxies in the context of a halo model
and the clustering of the underlying halos in redshift space using a
combination of Eulerian perturbation theory and $N$body simulations. The
modeling of redshiftspace distortions is done using the socalled distribution
function approach. The final model has 13 free parameters, and each parameter
is physically motivated rather than a nuisance parameter, which allows the use
of wellmotivated priors. We account for the FingerofGod effect from centrals
and both isolated and nonisolated satellites rather than using a single
velocity dispersion to describe the combined effect. We test and validate the
accuracy of the model on several sets of highfidelity $N$body simulations, as
well as realistic mock catalogs designed to simulate the BOSS DR12 CMASS data
set. The suite of simulations covers a range of cosmologies and galaxy bias
models, providing a rigorous test of the level of theoretical systematics
present in the model. The level of bias in the recovered values of $f \sigma_8$
is found to be small. When including scales to $k = 0.4 \ h\mathrm{Mpc}^{1}$,
we find 1530\% gains in the statistical precision of $f \sigma_8$ relative to
$k = 0.2 \ h\mathrm{Mpc}^{1}$ and a roughly 1015\% improvement for the
perpendicular AlcockPaczynski parameter $\alpha_\perp$. Using the BOSS DR12
CMASS mocks as a benchmark for comparison, we estimate an uncertainty on $f
\sigma_8$ that is $\sim$1020\% larger than other similar Fourierspace RSD
models in the literature that use $k \leq 0.2 \ h\mathrm{Mpc}^{1}$, suggesting
that these models likely have a toolimited parametrization.

Measurements of lineofsight dependent clustering via the galaxy power
spectrum's multipole moments constitute a powerful tool for testing theoretical
models in largescale structure. Recent work shows that this measurement,
including a moving lineofsight, can be accelerated using Fast Fourier
Transforms (FFTs) by decomposing the Legendre polynomials into products of
Cartesian vectors. Here, we present a faster, optimal means of using FFTs for
this measurement. We avoid redundancy present in the Cartesian decomposition by
using a spherical harmonic decomposition of the Legendre polynomials.
Consequently, our method is substantially faster: a given multipole of order
$\ell$ requires only $2\ell+1$ FFTs rather than the $(\ell+1)(\ell+2)/2$ FFTs
of the Cartesian approach. For the hexadecapole ($\ell = 4$), this translates
to $40\%$ fewer FFTs, with increased savings for higher $\ell$. The reduction
in wallclock time enables the calculation of finelybinned wedges in
$P(k,\mu)$, obtained by computing multipoles up to a large $\ell_{\rm max}$ and
combining them. This transformation has a number of advantages. We demonstrate
that by using nonuniform bins in $\mu$, we can isolate planeofsky (angular)
systematics to a narrow bin at $\mu \simeq 0$ while eliminating the
contamination from all other bins. We also show that the covariance matrix of
clustering wedges binned uniformly in $\mu$ becomes illconditioned when
combining multipoles up to large values of $\ell_{\rm max}$, but that the
problem can be avoided with nonuniform binning. As an example, we present
results using $\ell_{\rm max}=16$, for which our procedure requires a factor of
3.4 fewer FFTs than the Cartesian method, while removing the first $\mu$ bin
leads only to a 7% increase in statistical error on $f \sigma_8$, as compared
to a 54% increase with $\ell_{\rm max}=4$.

The use of large, Xray selected galaxy cluster catalogues for cosmological
analyses requires a thorough understanding of the Xray mass estimates. Weak
gravitational lensing is an ideal method to shed light on such issues, due to
its insensitivity to the cluster dynamical state. We perform a weak lensing
calibration of 166 galaxy clusters from the REFLEX and BCS cluster catalogue
and compare our results to the Xray masses based on scaled luminosities from
that catalogue. To interpret the weak lensing signal in terms of cluster
masses, we compare the lensing signal to simple theoretical NavarroFrenkWhite
models and to simulated cluster lensing profiles, including complications such
as cluster substructure, projected largescale structure, and Eddington bias.
We find evidence of underestimation in the Xray masses, as expected, with
$\langle M_{\mathrm{X}}/M_{\mathrm{WL}}\rangle = 0.75 \pm 0.07$ stat. $\pm
0.05$ sys. for our bestfit model. The biases in cosmological parameters in a
typical cluster abundance measurement that ignores this mass bias will
typically exceed the statistical errors.

Thermal SunyaevZeldovich (tSZ) power spectrum is one of the most sensitive
methods to constrain cosmological parameters, scaling as the amplitude
$\sigma_8^8$. It is determined by the integral over the halo mass function
multiplied by the total pressure content of clusters, and further convolved by
the cluster gas pressure profile. It has been shown that various feedback
effects can change significantly the pressure profile, possibly even pushing
the gas out to the virial radius and beyond, strongly affecting the tSZ power
spectrum at high $l$. Energetics arguments and SZhalo mass scaling relations
suggest feedback is unlikely to significantly change the total pressure
content, making low $l$ tSZ power spectrum more robust against feedback
effects. Furthermore, the separation between the cosmic infrared background
(CIB) and tSZ is more reliable at low $l$. Low $l$ modes are however probing
very small volumes, giving rise to very large nongaussian sampling variance
errors. By computing the trispectrum contribution we identify $90<l<350$ as the
minimum variance scale where the combined error is minimized. Using the
measurement at this $l$ provides constraints that are almost as strong as using
the entire power spectrum. We extend the work by the Planck collaboration by
including a full trispectrum and feedback effects in the analysis. We perform a
MarkovChain Monte Carlo analysis over the two dimensional parameter space and
find constraints on $\sigma_8$ by marginalizing over the feedback nuisance
parameter.We obtain $\sigma_8
=0.820^{+0.021}_{0.009}\left(\Omega_m/0.31\right)^{0.4}$, when fixing other
parameters to Planck cosmology values. Our results suggest it is possible to
obtain competitive cosmological constraints from tSZ without cluster redshift
information, and that the current tSZ power spectrum shows no evidence for a
low amplitude of $\sigma_8$.

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 several methods to accurately estimate Lagrangian bias parameters,
in particular the quadratic terms, both the local and the non local ones and
show the first clear evidence for the latter. Using Fourier space correlations,
we also show for the first time, the scale dependence of the quadratic and
nonlocal bias coefficients. We fit for the scale dependence of linear bias and
demonstrate the validity of a consistency relation between linear bias
parameters. Furthermore we employ real space estimators, using both
crosscorrelations and the PeakBackground Split argument. This is the first
time the latter is used to measure anisotropic bias coefficients. We find good
agreement among the methods, and also good agreement for local bias with
ESP$\tau$ theory predictions. Possible relations among the different bias
parameters are exploited. Finally, we also show how including higher order bias
reduces the magnitude and scale dependence of stochasticity of the halo field.

The Dark Energy Spectroscopic Instrument (DESI), a multiplexed fiberfed
spectrograph, is a StageIV groundbased dark energy experiment aiming to
measure redshifts for 29 million EmissionLine Galaxies (ELG), 4 million
Luminous Red Galaxies (LRG), and 2 million QuasiStellar Objects (QSO). The
survey design includes a pattern of tiling on the sky and the locations of the
fiber positioners in the focal plane of the telescope, with the observation
strategy determined by a fiber assignment algorithm that optimizes the
allocation of fibers to targets. This strategy allows a given region to be
covered on average five times for a fiveyear survey, but with coverage varying
between zero and twelve, which imprints a spatiallydependent pattern on the
galaxy clustering. We investigate the systematic effects of the fiber
assignment coverage on the anisotropic galaxy clustering of ELGs and show that,
in the absence of any corrections, it leads to discrepancies of order ten
percent on large scales for the power spectrum multipoles. We introduce a
method where objects in a random catalog are assigned a coverage, and the mean
density is separately computed for each coverage factor. We show that this
method reduces, but does not eliminate the effect. We next investigate the
angular dependence of the contaminated signal, arguing that it is mostly
localized to purely transverse modes. We demonstrate that the cleanest way to
remove the contaminating signal is to perform an analysis of the anisotropic
power spectrum $P(k,\mu)$ and remove the lowest $\mu$ bin, leaving $\mu>0$
modes accurate at the fewpercent level. Here, $\mu$ is the cosine of the angle
between the lineofsight and the direction of $\vec{k}$. We also investigate
two alternative definitions of the random catalog and show they are comparable
but less effective than the coverage randoms method.

We evaluate the covariance matrix of the matter power spectrum using
perturbation theory up to dominant terms at 1loop order and compare it to
numerical simulations. We decompose the covariance matrix into the disconnected
(Gaussian) part, trispectrum from the modes outside the survey (beat coupling
or supersample variance), and trispectrum from the modes inside the survey,
and show how the different components contribute to the overall covariance
matrix. We find the agreement with the simulations is at a 10\% level up to $k
\sim 1 h {\rm Mpc^{1}}$. We show that all the connected components are
dominated by the largescale modes ($k<0.1 h {\rm Mpc^{1}}$), regardless of
the value of the wavevectors $k,\, k'$ of the covariance matrix, suggesting
that one must be careful in applying the jackknife or bootstrap methods to the
covariance matrix. We perform an eigenmode decomposition of the connected part
of the covariance matrix, showing that at higher $k$ it is dominated by a
single eigenmode. The full covariance matrix can be approximated as the
disconnected part only, with the connected part being treated as an external
nuisance parameter with a known scale dependence, and a known prior on its
variance for a given survey volume. Finally, we provide a prescription for how
to evaluate the covariance matrix from small box simulations without the need
to simulate large volumes.

A strong instrumentation and detector R&D program has enabled the current
generation of cosmic frontier surveys. A small investment in R&D will continue
to pay dividends and enable new probes to investigate the accelerated expansion
of the universe. Instrumentation and detector R&D provide critical training
opportunities for future generations of experimentalists, skills that are
important across the entire Department of Energy High Energy Physics program.

We explore the relationship between the nonlinear matter power spectrum and
the various Lagrangian and Standard Perturbation Theories (LPT and SPT). We
first look at it in the context of one dimensional (1d) dynamics, where 1LPT
is exact at the perturbative level and one can exactly resum the SPT series
into the 1LPT power spectrum. Shell crossings lead to nonperturbative effects,
and the PT ignorance can be quantified in terms of their ratio, which is also
the transfer function squared in the absence of stochasticity. At the order of
PT we work, this parametrization is equivalent to the results of effective
field theory (EFT), and can thus be expanded in terms of the same parameters.
We find that its radius of convergence is larger than the SPT loop expansion.
The same EFT parametrization applies to all SPT loop terms and, if
stochasticity can be ignored, to all Npoint correlators. In 3d, the LPT
structure is considerably more complicated, and we find that LPT models with
parametrization motivated by the EFT exhibit running with $k$ and that SPT is
generally a better choice. Since these transfer function expansions contain
free parameters that change with cosmological model their usefulness for
broadband power is unclear. For this reason we test the predictions of these
models on baryonic acoustic oscillations (BAO) and other primordial
oscillations, including string monodromy models, for which we ran a series of
simulations with and without oscillations. Most models are successful in
predicting oscillations beyond their corresponding PT versions, confirming the
basic validity of the model.

Cosmic surveys provide crucial information about high energy physics
including strong evidence for dark energy, dark matter, and inflation. Ongoing
and upcoming surveys will start to identify the underlying physics of these new
phenomena, including tight constraints on the equation of state of dark energy,
the viability of modified gravity, the existence of extra light species, the
masses of the neutrinos, and the potential of the field that drove inflation.
Even after the Stage IV experiments, DESI and LSST, complete their surveys,
there will still be much information left in the sky. This additional
information will enable us to understand the physics underlying the dark
universe at an even deeper level and, in case Stage IV surveys find hints for
physics beyond the current Standard Model of Cosmology, to revolutionize our
current view of the universe. There are many ideas for how best to supplement
and aid DESI and LSST in order to access some of this remaining information and
how surveys beyond Stage IV can fully exploit this regime. These ideas flow to
potential projects that could start construction in the 2020's.

Amongst standard model parameters that are constrained by cosmic microwave
background (CMB) observations, the optical depth $\tau$ stands out as a
nuisance parameter. While $\tau$ provides some crude limits on reionization, it
also degrades constraints on other cosmological parameters. Here we explore how
21 cm cosmologyas a direct probe of reionizationcan be used to
independently predict $\tau$ in an effort to improve CMB parameter constraints.
We develop two complementary schemes for doing so. The first uses 21 cm power
spectrum observations in conjunction with semianalytic simulations to predict
$\tau$. The other uses global 21 cm measurements to directly constrain low
redshift (postreheating) contributions to $\tau$ in a relatively
modelindependent way. Forecasting the performance of the upcoming Hydrogen
Epoch of Reionization Array, we find that significant reductions in the errors
on $\tau$ can be achieved. These results are particularly effective at breaking
the CMB degeneracy between $\tau$ and the amplitude of the primordial
fluctuation spectrum $A_s$, with errors on $\ln (10^{10} A_s)$ reduced by up to
a factor of four. Stage 4 CMB constraints on the neutrino mass sum are also
improved, with errors potentially reduced to $12\,\textrm{meV}$ regardless of
whether CMB experiments can precisely measure the reionization bump in
polarization power spectra. Observations of the 21 cm line are therefore
capable of improving not only our understanding of reionization astrophysics,
but also of cosmology in general.

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct
novel cosmological observations using the BOSS spectrograph at Apache Point
Observatory. Observations will be simultaneous with the Time Domain
Spectroscopic Survey (TDSS) designed for variability studies and the
Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for
studies of Xray sources. eBOSS will use four different tracers to measure the
distanceredshift relation with baryon acoustic oscillations (BAO). Using more
than 250,000 new, spectroscopically confirmed luminous red galaxies at a median
redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to
an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the
z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy
redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and
H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000
spectroscopicallyconfirmed quasars will provide the first BAO distance
measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8%
and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars
and reobservation of 60,000 quasars known from BOSS, we will obtain new
Lymanalpha forest measurements at redshifts z>2.1; these new data will enhance
the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS.
Furthermore, eBOSS will provide improved tests of General Relativity on
cosmological scales through redshiftspace distortion measurements, improved
tests for nonGaussianity in the primordial density field, and new constraints
on the summed mass of all neutrino species. Here, we provide an overview of the
cosmological goals, spectroscopic target sample, demonstration of spectral
quality from early data, and projected cosmological constraints from eBOSS.

The standard formalism for the coevolution of halos and dark matter predicts
that any initial halo velocity bias rapidly decays to zero. We argue that, when
the purpose is to compute statistics like power spectra etc., the coupling in
the momentum conservation equation for the biased tracers must be modified. Our
new formulation predicts the constancy in time of any statistical halo velocity
bias present in the initial conditions, in agreement with peak theory. We test
this prediction by studying the evolution of a conserved halo population in
Nbody simulations. We establish that the initial simulated halo density and
velocity statistics show distinct features of the peak model and, thus, deviate
from the simple local Lagrangian bias. We demonstrate, for the first time, that
the time evolution of their velocity is in tension with the rapid decay
expected in the standard approach.