
The observed power spectrum in redshift space appears distorted due to the
peculiar motion of galaxies, known as redshiftspace distortions (RSD). While
all the effects in RSD are accounted for by the simple mapping formula from
real to redshift spaces, accurately modeling redshiftspace power spectrum is
rather difficult due to the nonperturbative properties of the mapping. Still,
however, a perturbative treatment may be applied to the power spectrum at
largescales, and on top of a careful modeling of the FingerofGod effect
caused by the smallscale random motion, the redshiftspace power spectrum can
be expressed as a series of expansion which contains the higherorder
correlations of density and velocity fields. In our previous work [JCAP 8
(Aug., 2016) 050], we provide a perturbationtheory inspired model for power
spectrum in which the higherorder correlations are evaluated directly from the
cosmological $N$body simulations. Adopting a simple Gaussian ansatz for
FingerofGod effect, the model is shown to quantitatively describe the
simulation results. Here, we further push this approach, and present an
accurate power spectrum template which can be used to estimate the growth of
structure as a key to probe gravity on cosmological scales. Based on the
simulations, we first calibrate the uncertainties and systematics in the
pertrubation theory calculation in a fiducial cosmological model. Then, using
the scaling relations, the calibrated power spectrum template is applied to a
different cosmological model. We demonstrate that with our new template, the
bestfitted growth functions are shown to reproduce the fiducial values in a
good accuracy of 1 \% at $k<0.18 \hompc$ for cosmologies with different Hubble
parameters.

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.

The mapping of dark matter clustering from real space to redshift space
introduces the anisotropic property to the measured density power spectrum in
redshift space, known as the redshift space distortion effect. The mapping
formula is intrinsically nonlinear, which is complicated by the higher order
polynomials due to indefinite cross correlations between the density and
velocity fields, and the FingerofGod effect due to the randomness of the
peculiar velocity field. Whilst the full higher order polynomials remain
unknown, the other systematics can be controlled consistently within the same
order truncation in the expansion of the mapping formula, as shown in this
paper. The systematic due to the unknown nonlinear density and velocity fields
is removed by separately measuring all terms in the expansion directly using
simulations. The uncertainty caused by the velocity randomness is controlled by
splitting the FoG term into two pieces, 1) the "onepoint" FoG term being
independent of the separation vector between two different points, and 2) the
"correlated" FoG term appearing as an indefinite polynomials which is expanded
in the same order as all other perturbative polynomials. Using 100 realizations
of simulations, we find that the Gaussian FoG function with only one
scaleindependent free parameter works quite well, and that our new mapping
formulation accurately reproduces the observed 2dimensional density power
spectrum in redshift space at the smallest scales by far, up to $k\sim
0.2h$Mpc, considering the resolution of future experiments.

The anisotropic galaxy clustering of large scale structure observed by the
Baryon Oscillation Spectroscopic Survey Data Release 11 is analyzed to probe
the sum of neutrino mass in the small $m_\nu < 1eV$ limit in which the early
broadband shape determined before the last scattering surface is immune from
the variation of $m_\nu$. The signature of $m_\nu$ is imprinted on the altered
shape of the power spectrum at later epoch, which provides an opportunity to
access the nontrivial $m_\nu$ through the measured anisotropic correlation
function in redshift space (hereafter RSD instead of Redshift Space
Distortion). The nonlinear RSD corrections with massive neutrinos in the
quasi linear regime are approximately estimated using oneloop order terms
computed by tomographic linear solutions. We suggest a new approach to probe
$m_\nu$ simultaneously with all other distance measures and coherent growth
functions, exploiting this deformation of the early broadband shape of the
spectrum at later epoch. If the origin of cosmic acceleration is unknown,
$m_\nu$ is poorly determined after marginalising over all other observables.
However, we find that the measured distances and coherent growth functions are
minimally affected by the presence of mild neutrino mass. Although the standard
model of cosmic acceleration is assumed to be the cosmological constant, the
constraint on $m_\nu$ is little improved. Interestingly, the measured CMB
distance to the last scattering surface sharply slices the degeneracy between
the matter content and $m_\nu$, and the hidden $m_\nu$ is excavated to be
$m_\nu=0.19^{+0.28}_{0.17} eV$ which is different from massless neutrino more
than 68% confidence.

The cosmic distance can be precisely determined using a `standard ruler'
imprinted by primordial baryon acoustic oscillation (hereafter BAO) in the
early Universe. The BAO at the targeted epoch is observed by analysing galaxy
clustering in redshift space (hereafter RSD) for which a theoretical
formulation is not yet fully understood, and thus makes this methodology
unsatisfactory. The BAO analysis following a full RSD modelling is contaminated
by systematic uncertainties due to a nonlinear smearing effects such as
nonlinear corrections and by random viral velocity of galaxies. However, the
BAO can be probed independently of RSD contamination using the BAO peak
positions located in the 2D anisotropic correlation function. A new methodology
is presented to measure peak positions, to test whether it is also contaminated
by the same systematics in RSD, and to provide the radial and transverse cosmic
distances determined by the 2D BAO peak positions. We find that in our model
independent anisotropic clustering analysis we can obtain about $2\%$ and $5\%$
constraints on $D_A$ and $H^{1}$ respectively with current BOSS data, which is
competitive with other analysis.

We analyse the clustering of cosmic large scale structure using a consistent
modified gravity perturbation theory, accounting for anisotropic effects along
and transverse to the line of sight. The growth factor has a particular scale
dependence in f(R) gravity and we fit for the shape parameter f_{R0}
simultaneously with the distance and the large scale (general relativity) limit
of the growth function. Using more than 690,000 galaxies in the Baryon
Oscillation Spectroscopy Survey Data Release 11, we find no evidence for extra
scale dependence, with the 95\% confidence upper limit f_{R0} <8 \times
10^{4}. Future clustering data, such as from the Dark Energy Spectroscopic
Instrument, can use this consistent methodology to impose tighter constraints.

The apparent anisotropies of the galaxy clustering in observable redshift
space provide a unique opportunity to simultaneously probe cosmic expansion and
gravity on cosmological scales via the AlcockPaczynski effect and
redshiftspace distortions. While the improved theoretical models have been
proposed and developed to describe the apparent anisotropic clustering at
weakly nonlinear scales, the applicability of these models is still limited in
the presence of the nonperturbative smearing effect caused by the randomness
of the relative velocities. Although the cosmological constraint from the
anisotropic clustering will be certainly improved with a more elaborate
theoretical model, we here consider an alternative way by using the statistical
power of both the power spectrum and bispectrum at large scales. Based on the
Fisher matrix analysis, we estimate the benefit of combining the power spectra
and bispectra, finding that the constraints on the cosmic expansion and growth
of structure will be improved by a factor of two. This compensates for the loss
of constraining power using the power spectrum alone due to the randomness of
the relative velocities.

SPHEREx (SpectroPhotometer for the History of the Universe, Epoch of
Reionization, and Ices Explorer) ( http://spherex.caltech.edu ) is a proposed
allsky spectroscopic survey satellite designed to address all three science
goals in NASA's Astrophysics Division: probe the origin and destiny of our
Universe; explore whether planets around other stars could harbor life; and
explore the origin and evolution of galaxies. SPHEREx will scan a series of
Linear Variable Filters systematically across the entire sky. The SPHEREx data
set will contain R=40 spectra fir 0.75$<\lambda<$4.1$\mu$m and R=150 spectra
for 4.1$<\lambda<$4.8$\mu$m for every 6.2 arc second pixel over the entiresky.
In this paper, we detail the extragalactic and cosmological studies SPHEREx
will enable and present detailed systematic effect evaluations. We also outline
the Ice and Galaxy Evolution Investigations.

We analyze the clustering of large scale structure in the Universe in a model
independent method, accounting for anisotropic effects along and transverse to
the line of sight. The Baryon Oscillation Spectroscopy Survey Data Release 11
provides a large sample of 690,000 galaxies, allowing determination of the
Hubble expansion H, angular distance D_A, and growth rate G_T at an effective
redshift of z=0.57. After careful bias and convergence studies of the effects
from small scale clustering, we find that cutting transverse separations below
40 Mpc/h delivers robust results while smaller scale data leads to a bias due
to unmodelled nonlinear and velocity effects. The converged results are in
agreement with concordance LCDM cosmology, general relativity, and minimal
neutrino mass, all within the 68% confidence level. We also present results
separately for the northern and southern hemisphere sky, finding a slight
tension in the growth rate  potentially a signature of anisotropic stress, or
just covariance with small scale velocities  but within 68% CL.

Baryon acoustic oscillations (BAO), known as one of the largest cosmological
objects, is now recognized as standard cosmological tool to measure geometric
distances via the AlcockPaczynski effect, by which the observed BAO exhibits
characteristic anisotropies in addition to the redshift distortions. This
implies that once we know the correct distances to the observed BAO, the tip
points of baryon acoustic peaks in the anisotropic correlation function of
galaxies, $\xi(\sigma,\pi)$, can form a great circle (hereafter 2D BAO circle)
in the $\sigma$ and $\pi$ plane, where $\sigma$ and $\pi$ are the separation of
galaxy pair parallel and perpendicular to the lineofsight, respectively. This
2D BAO circle remains unchanged under the variations of the unknown galaxy bias
and/or coherent motion, while it varies transversely and radially with respect
to the variations of $D_A$ and $H^{1}$, respectively. Hereby the ratio between
transverse distance $D_A$ and the radial distance $H^{1}$ reproduces the
intrinsic shape of 2D BAO circle, which is {\it a priori} given by the known
broadband shape of spectra. All BAO peaks of $\xi(\sigma,\pi)$ are precisely
calculated with the improved theoretical model of redshift distortion. We test
this broadband AlcockPaczynski method using BOSSlike mock catalogues. The
transverse and radial distances are probed in precision of several percentage
fractional errors, and the coherent motion is observed to match with the
fiducial values accurately.

Our observations of the Universe are fundamentally anisotropic, with data
from galaxies separated transverse to the line of sight coming from the same
epoch while that from galaxies separated parallel to the line of sight coming
from different times. Moreover, galaxy velocities along the line of sight
change their redshift, giving redshift space distortions. We perform a full
twodimensional anisotropy analysis of galaxy clustering data, fitting in a
substantially model independent manner the angular diameter distance D_A,
Hubble parameter H, and growth rate ddelta/dln a without assuming a dark energy
model. The results demonstrate consistency with LCDM expansion and growth,
hence also testing general relativity. We also point out the interpretation
dependence of the effective redshift z_eff, and its cosmological impact for
next generation surveys.

The cosmological power spectrum of the coherent matter flow is measured
exploiting an improved prescription for the apparent anisotropic clustering
pattern in redshift space. New statistical analysis is presented to provide an
optimal observational platform to link the improved redshift distortion
theoretical model to future real datasets. The statistical power as well as
robustness of our method are tested against 60 realizations of 8 Gpc/h^3 dark
matter simulation maps mocking the precision level of upcoming widedeep
surveys. We showed that we can accurately extract the velocity power spectrum
up to quasi linear scales of k~0.1 h/Mpc at z = 0.35 and up to k~0.15 h/Mpc at
higher redshifts within a couple of percentage precision level. Our
understanding of redshift space distortion is proved to be appropriate for
precision cosmology, and our statistical method will guide us to righteous path
to meet the real world.

We investigate the direct determination of expansion history using redshift
distortions without plugging into detailed cosmological parameters. The
observed spectra in redshift space include a mixture of information:
fluctuations of densitydensity and velocityvelocity spectra, and distance
measures of perpendicular and parallel components to the line of sight.
Unfortunately it is hard to measure all the components simultaneously without
any specific prior assumption. Common prior assumptions include a
linear/quasilinear model of redshift distortions or a model for the shape of
the power spectra, which eventually breaks down on small scales at later epochs
where nonlinear structure formation disturbs coherent growth. The degeneracy
breaking, between the effect of cosmic distances and redshift distortions for
example, depends on the prior we assume. An alternative approach is to utilize
the cosmological principle inscribed in the heart of the
FriedmannLematreRobertsonWalker (hereafter FLRW) universe, that is, the
specific relation between the angular diameter distance and the Hubble
parameter, in this degeneracy breaking. We show that utilizing this FLRW prior
early in the step of distinguishing the distance effect from redshift
distortions helps us improve the detectability of power spectra and distance
measures with no leaning on a combination of other experiments.

We show that in the single component situation all perturbation variables in
the comoving gauge are conformally invariant to all perturbation orders.
Generally we identify a special time slicing, the uniformconformal
transformation slicing, where all perturbations are again conformally invariant
to all perturbation orders. We apply this result to the delta N formalism, and
show its conformal invariance.

New statistical method is proposed to coherently combine Baryon Acoustic
Oscillation statistics (BAO) and peculiar velocity measurements exploiting
decomposed densitydensity and velocityvelocity spectra in real space from the
observed redshift distortions in redshift space, 1) to achieve stronger dark
energy constraints, \sigma(w)=0.06 and \sigma(w_a)=0.20, which are enhanced
from BAO or velocity measurements alone, and 2) to crosscheck consistency of
dark energy constraints from two different approaches; BAO as geometrical
measurements and peculiar velocity as large scale structure formation
observables. In addition to those advantages, as power spectra decomposition
procedure is free from uncertainty of galaxy bias, this simultaneous fitting is
an optimal method to extract cosmological parameters without any preassumption
about galaxy bias.

We present new measurements of the coherent motion of galaxies based on
observations of the largescale redshiftspace distortions seen in the
twodimensional twopoint correlation function of Luminous Red Galaxies in Data
Release Seven of the Sloan Digital Sky Survey. We have developed a new
methodology for estimating these coherent motions, which is less dependent on
the details of galaxy bias and of the cosmological model to explain the
latetime acceleration of the expansion of the Universe. We measure a
onedimensional velocity dispersion of galaxies on largescales of
\sigma_v=3.01^{+0.45}_{0.46} Mpc/h and \sigma_v=3.69^{+0.47}_{0.47} \mpcoh$
at a mean redshift of z=0.25 and 0.38 respectively. These values are fully
consistent with predictions for a WMAP7normalised LCDM Universe and
inconsistent at confidence of 3.8\sigma with a DvaliGabadadzePorrati (DGP)
model for the Universe. We can convert the units of these $\sigma_v$
measurements to 270^{+40}_{41} km/s and 320^{+41}_{41} km/s respectively
(assuming a $\Lambda$CDM universe), which are lower that expected based on
recent low redshift (z<0.2) measurements of the peculiar velocity field (or
"bulk flows"). It is difficult to directly compare these measurements as they
cover different redshift ranges and different areas of the sky. However, one
possible cosmological explanation for this discrepancy is that our Galaxy is
located in unusually over, or under, dense region of the Universe.

Using the latest physical modeling and constrained by the most recent data,
we develop a phenomenological parameterized model of the contributions to
intensity and polarization maps at millimeter wavelengths from external
galaxies and SunyaevZeldovich effects. We find such modeling to be necessary
for estimation of cosmological parameters from Planck data. For example,
ignoring the clustering of the infrared background would result in a bias in
n_s of 7 sigma. We show that the simultaneous marginalization over a full
foreground model can eliminate such biases, while increasing the statistical
uncertainty in cosmological parameters by less than 20%. The small increases in
uncertainty can be significantly reduced with the inclusion of
higherresolution groundbased data.
The multifrequency analysis we employ involves modeling 46 total power
spectra and marginalization over 17 foreground parameters. We show that we can
also reduce the data to a best estimate of the CMB power spectra, and just two
principal components (with constrained amplitudes) describing residual
foreground contamination.

We explore the complementarity of weak lensing and galaxy peculiar velocity
measurements to better constrain modifications to General Relativity. We find
no evidence for deviations from GR on cosmological scales from a combination of
peculiar velocity measurements (for Luminous Red Galaxies in the Sloan Digital
Sky Survey) with weak lensing measurements (from the CFHT Legacy Survey). We
provide a Fisher error forecast for a Euclidlike spacebased survey including
both lensing and peculiar velocity measurements, and show that the expected
constraints on modified gravity will be at least an order of magnitude better
than with present data, i.e. we will obtain 5% errors on the modified gravity
parametrization described here. We also present a modelindependent method for
constraining modified gravity parameters using tomographic peculiar velocity
information, and apply this methodology to the present dataset.

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 develop an optimized technique to extract densitydensity and
velocityvelocity spectra out of observed spectra in redshift space. The
measured spectra of the distribution of halos from redshift distorted mock map
are binned into 2dimensional coordinates in Fourier space so as to be
decomposed into both spectra using angular projection dependence. With the
threshold limit introduced to minimize nonlinear suppression, the decomposed
velocityvelocity spectra are reasonably well measured up to scale k=0.07
h/Mpc, and the measured variances using our method are consistent with errors
predicted from a Fisher matrix analysis. The detectability is extendable to
k\sim 0.1 h/Mpc with more conservative bounds at the cost of weakened
constraint.

We present here a new parameterization for the bulk motions of galaxies and
clusters (in the linear regime) that can be measured statistically from the
shape and amplitude of the twodimensional twopoint correlation function. We
further propose the onedimensional velocity dispersion (v_p) of the bulk flow
as a complementary measure of redshiftspace distortions, which is
modelindependent and not dependent on the normalisation method. As a
demonstration, we have applied our new methodology to the C4 cluster catalogue
constructed from Data Release Three (DR3) of the Sloan Digital Sky Survey. We
find v_p=270^{+433}km/s (also consistent with v_p=0) for this cluster sample
(at z=0.1), which is in agreement with that predicted for a WMAP5normalised
LCDM model (i.e., v_p(LCDM=203km/s). This measurement does not lend support to
recent claims of excessive bulk motions (\simeq1000 km/s) which appear in
conflict with LCDM, although our large statistical error cannot rule them out.
From the measured coherent evolution of v_p, we develop a technique to
reconstruct the perturbed potential, as well as estimating the unbiased matter
density fluctuations and scaleindependent bias.

Next generation surveys will observe the largescale structure of the
Universe with unprecedented accuracy. This will enable us to test the
relationships between matter overdensities, the curvature perturbation and the
Newtonian potential. Any largedistance modification of gravity or exotic
nature of dark energy modifies these relationships as compared to those
predicted in the standard smooth dark energy model based on General Relativity.
In linear theory of structure growth such modifications are often parameterised
by virtue of two functions of space and time that enter the relation of the
curvature perturbation to, first, the matter overdensity, and second, the
Newtonian potential. We investigate the predictions for these functions in
BransDicke theory, clustering dark energy models and interacting dark energy
models. We find that each theory has a distinct path in the parameter space of
modified growth. Understanding these theoretical priors on the
parameterisations of modified growth is essential to reveal the nature of
cosmic acceleration with the help of upcoming observations of structure
formations.

Dark energy in General Relativity is typically noninteracting with other
matter. However, it is possible that the dark energy interacts with the dark
matter, and in this case, the dark matter can violate the universality of free
fall (the weak equivalence principle). We show that some forms of the dark
sector interaction do not violate weak equivalence. For those interactions that
do violate weak equivalence, there are no available laboratory experiments to
probe this violation for dark matter. But cosmology provides a test for
violations of the equivalence principle between dark matter and baryons  via
a test for consistency of the observed galaxy velocities with the Euler
equation.

Measuring the statistics of galaxy peculiar velocities using redshiftspace
distortions is an excellent way of probing the history of structure formation.
Because galaxies are expected to act as test particles within the flow of
matter, this method avoids uncertainties due to an unknown galaxy density bias.
We show that the parameter combination measured by redshiftspace distortions,
$f\sigma_8^{\rm mass}$ provides a good test of dark energy models, even without
the knowledge of bias or $\sigma_8^{\rm mass}$ required to extract $f$ from
this measurement (here $f$ is the logarithmic derivative of the linear growth
rate, and $\sigma_8^{\rm mass}$ is the rootmeansquare mass fluctuation in
spheres with radius $8h^{1}$Mpc). We argue that redshiftspace distortion
measurements will help to determine the physics behind the cosmic acceleration,
testing whether it is related to dark energy or modified gravity, and will
provide an opportunity to test possible dark energy clumping or coupling
between dark energy and dark matter. If we can measure galaxy bias in addition,
simultaneous measurement of both the overdensity and velocity fields can be
used to test the validity of equivalence principle, through the continuity
equation.