The observed power spectrum in redshift space appears distorted due to the peculiar motion of galaxies, known as redshift-space distortions (RSD). While all the effects in RSD are accounted for by the simple mapping formula from real to redshift spaces, accurately modeling redshift-space power spectrum is rather difficult due to the non-perturbative properties of the mapping. Still, however, a perturbative treatment may be applied to the power spectrum at large-scales, and on top of a careful modeling of the Finger-of-God effect caused by the small-scale random motion, the redshift-space power spectrum can be expressed as a series of expansion which contains the higher-order correlations of density and velocity fields. In our previous work [JCAP 8 (Aug., 2016) 050], we provide a perturbation-theory inspired model for power spectrum in which the higher-order correlations are evaluated directly from the cosmological $N$-body simulations. Adopting a simple Gaussian ansatz for Finger-of-God 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 best-fitted 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 ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area 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 magnitude-limited 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 ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree 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 non-linear, which is complicated by the higher order polynomials due to indefinite cross correlations between the density and velocity fields, and the Finger-of-God 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 non-linear 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 "one-point" 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 scale-independent free parameter works quite well, and that our new mapping formulation accurately reproduces the observed 2-dimensional 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 non--trivial $m_\nu$ through the measured anisotropic correlation function in redshift space (hereafter RSD instead of Redshift Space Distortion). The non--linear RSD corrections with massive neutrinos in the quasi linear regime are approximately estimated using one-loop 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 non-linear smearing effects such as non-linear 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 Alcock--Paczynski effect and redshift-space distortions. While the improved theoretical models have been proposed and developed to describe the apparent anisotropic clustering at weakly non-linear scales, the applicability of these models is still limited in the presence of the non--perturbative 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 (Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer) ( http://spherex.caltech.edu ) is a proposed all-sky 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 entire-sky. In this paper, we detail the extra-galactic 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 Alcock-Paczynski 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 line-of-sight, 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 Alcock--Paczynski method using BOSS--like 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 two-dimensional 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 wide--deep 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 density-density and velocity-velocity 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/quasi-linear 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 Friedmann-Lematre-Robertson-Walker (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 uniform-conformal 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 density-density and velocity-velocity 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 cross-check 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 pre-assumption about galaxy bias.

We present new measurements of the coherent motion of galaxies based on observations of the large-scale redshift-space distortions seen in the two-dimensional two-point 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 late-time acceleration of the expansion of the Universe. We measure a one-dimensional velocity dispersion of galaxies on large-scales 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 WMAP7-normalised LCDM Universe and inconsistent at confidence of 3.8\sigma with a Dvali-Gabadadze-Porrati (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 Sunyaev-Zeldovich 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 higher-resolution ground-based data. The multi-frequency 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 Euclid-like space-based 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 model--independent method for constraining modified gravity parameters using tomographic peculiar velocity information, and apply this methodology to the present dataset.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE 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 multi-wavelength 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 ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster 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 density--density and velocity--velocity spectra out of observed spectra in redshift space. The measured spectra of the distribution of halos from redshift distorted mock map are binned into 2--dimensional 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 velocity--velocity 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 two-dimensional two-point correlation function. We further propose the one-dimensional velocity dispersion (v_p) of the bulk flow as a complementary measure of redshift-space distortions, which is model-independent 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 WMAP5-normalised 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 re-construct the perturbed potential, as well as estimating the unbiased matter density fluctuations and scale--independent bias.

Next generation surveys will observe the large-scale structure of the Universe with unprecedented accuracy. This will enable us to test the relationships between matter over-densities, the curvature perturbation and the Newtonian potential. Any large-distance 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 over-density, and second, the Newtonian potential. We investigate the predictions for these functions in Brans-Dicke 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 non-interacting 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 redshift-space 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 redshift-space 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 root-mean-square mass fluctuation in spheres with radius $8h^{-1}$Mpc). We argue that redshift-space 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.