
Observations of galaxies and galaxy clusters in the local universe can
account for only $\sim\,10\%$ of the total baryon content. Cosmological
simulations predict that the `missing baryons' are spread throughout
filamentary structures in the cosmic web, forming a lowdensity gas with
temperatures of $10^510^7\,\!$K. We search for this warmhot intergalactic
medium (WHIM) by stacking the Planck Compton $y$parameter map of the thermal
SunyaevZel'dovich (tSZ) effect for 1,002,334 pairs of CMASS galaxies from the
Sloan Digital Sky Survey. We model the contribution from the galaxy halo pairs
assuming spherical symmetry, finding a residual tSZ signal at the
$2.9\mbox{$\sigma$}$ level from a stacked filament of length $10.5\,h^{1}\,\rm
Mpc$ with a Compton parameter magnitude $y=(0.6\pm0.2)\times10^{8}$. We
consider possible sources of contamination and conclude that bound gas in
haloes may contribute only up to $20\%$ of the measured filamentary signal. To
estimate the filament gas properties we measure the gravitational lensing
signal for the same sample of galaxy pairs; in combination with the tSZ signal,
this yields an inferred gas density of $\rho_{\rm b}=(5.5\pm
2.9)\times\bar{\rho_{\rm b}}$ with a temperature $T=(2.7\pm 1.7) \times
10^6\,$K. This result is consistent with the predicted WHIM properties, and
overall the filamentary gas can account for $ 11\pm 7\%$ of the total baryon
content of the Universe. We also see evidence that the gas filament extends
beyond the galaxy pair. Averaging over this longer baseline boosts the
significance of the tSZ signal and increases the associated baryon content to
$28\pm 12\%$ of the global value.

We present a public suite of weak lensing mock data, extending the Scinet
Light Cone Simulations (SLICS) to simulate crosscorrelation analyses with
different cosmological probes. These mocks include KiDS450 and LSSTlike
lensing data, cosmic microwave background lensing maps and simulated
spectroscopic surveys that emulate the GAMA, BOSS and 2dFLenS galaxy surveys.
With 817 independent realisations, our mocks are optimised for combinedprobe
covariance estimation, which we illustrate for the case of a joint measurement
involving cosmic shear, galaxygalaxy lensing and galaxy clustering from
KiDS450 and BOSS data. With their high spatial resolution, the SLICS are also
optimal for predicting the signal for novel lensing estimators, for the
validation of analysis pipelines, and for testing a range of systematic effects
such as the impact of neighbourexclusion bias on the measured tomographic
cosmic shear signal. For surveys like KiDS and DES, where the rejection of
neighbouring galaxies occurs within ~2 arcseconds, we show that the measured
cosmic shear signal will be biased low, but by less than a percent on the
angular scales that are typically used in cosmic shear analyses. The amplitude
of the neighbourexclusion bias doubles in deeper, LSSTlike data. The
simulation products described in this paper are made available at
http://slics.roe.ac.uk/.

We study projected underdensities in the cosmic galaxy density field known as
'troughs', and their overdense counterparts, which we call 'ridges'. We
identify these regions using a bright sample of foreground galaxies from the
photometric KiloDegree Survey (KiDS), specifically selected to mimic the
spectroscopic Galaxy And Mass Assembly survey (GAMA). From an independent
sample of KiDS background galaxies, we measure the weak gravitational lensing
profiles of the troughs/ridges. We quantify their lensing strength $A$ as a
function of galaxy density percentile rank $P$ and overdensity $\delta$, and
find that the skewness in the galaxy density distribution is reflected in the
total mass distribution measured by weak lensing. We interpret our results
using the mock galaxy catalogue from the Marenostrum Institut de Ci\`encies de
l'Espai (MICE) Grand Challenge lightcone simulation, and find a good agreement
with our observations. Using signaltonoise weights derived from the Scinet
LIghtCone Simulations (SLICS) mock catalogue we optimally stack the lensing
signal of KiDS troughs with an angular radius $\theta_A$ = {5,10,15,20} arcmin,
resulting in {16.8,14.9,10.13,7.55} $\sigma$ detections. Finally, we select
troughs using a volumelimited sample of galaxies, split into two redshift bins
between 0.1 < z < 0.3. For troughs/ridges with transverse comoving radius $R_A$
= 1.9 Mpc/h, we find no significant difference between the comoving $A'(P)$ and
$A'(\delta)$ relation of the low and highredshift sample. Using the MICE and
SLICS mocks we predict that trough and ridge evolution could be detected with
gravitational lensing using deeper and wider lensing surveys, such as those
from the Large Synoptic Survey Telescope and Euclid.

We present a systematic comparison of several existing and new void finding
algorithms, focusing on their potential power to test a particular class of
modified gravity models  chameleon $f(R)$ gravity. These models deviate from
standard General Relativity (GR) more strongly in lowdensity regions and thus
voids are a promising venue to test them. We use Halo Occupation Distribution
(HOD) prescriptions to populate haloes with galaxies, and tune the HOD
parameters such that the galaxy twopoint correlation functions are the same in
both f(R) and GR models. We identify both 3D voids as well as 2D underdensities
in the planeofthesky to find the same void abundance and void galaxy number
density profiles across all models, which suggests that they do not contain
much information beyond galaxy clustering. However, the underlying void dark
matter density profiles are significantly different, with f(R) voids being more
underdense than GR ones, which leads to f(R) voids having a larger tangential
shear signal than their GR analogues. We investigate the potential of each void
finder to test f(R) models with nearfuture lensing surveys such as EUCLID and
LSST. The 2D voids have the largest power to probe f(R) gravity, with a LSST
analysis of tunnel (which is a new type of 2D underdensity introduced here)
lensing distinguishing at 80 and 11$\sigma$ (statistical error) f(R) models
with $f_{R0}=10^{5}$ and $10^{6}$ from GR.

In theories of modified gravity with the chameleon screening mechanism, the
strength of the fifth force depends on environment. This induces an environment
dependence of structure formation, which differs from $\Lambda$CDM. We show
that these differences can be captured by the marked correlation function. With
the galaxy correlation functions and number densities calibrated to match
between $f(R)$ and $\Lambda$CDM models in simulations, we show that the marked
correlation functions from using either the local density or halo mass as the
marks encode extra information, which can be used to test these theories. We
discuss possible applications of these statistics in observations.

We have searched for the signature of cosmic voids in the CMB, in both the
Planck temperature and lensingconvergence maps; voids should give decrements
in both. We use zobov voids from the DR12 SDSS CMASS galaxy sample. We base our
analysis on Nbody simulations, to avoid a posteriori bias. For the first time,
we detect the signature of voids in CMB lensing: the significance is
$3.2\sigma$, close to $\Lambda$CDM in both amplitude and projected
densityprofile shape. A temperature dip is also seen, at modest significance
($2.3\sigma$), with amplitude about 6 times the prediction. This temperature
signal is induced mostly by voids with radius between 100 and 150 Mpc/h, while
the lensing signal is mostly contributed by smaller voids  as expected;
lensing relates directly to density, while ISW depends on gravitational
potential. The void abundance in observations and simulations agree, as well.
We also repeated the analysis excluding lowersignificance voids: no lensing
signal is detected, with an upper limit of about twice the $\Lambda$CDM
prediction. But the mean temperature decrement now becomes nonzero at the
$3.7\sigma$ level (similar to that found by Granett et al.), with amplitude
about 20 times the prediction. However, the observed dependence of temperature
on void size is in poor agreement with simulations, whereas the lensing results
are consistent with $\Lambda$CDM theory. Thus, the overall tension between
theory and observations does not favour nonstandard theories of gravity,
despite the hints of an enhanced amplitude for the ISW effect from voids.

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 study the spherical evolution model for voids in $\Lambda$CDM, where the
evolution of voids is governed by dark energy at an earlier time than that for
the whole universe or in overdensities. We show that the presence of dark
energy suppresses the growth of peculiar velocities, causing void
shellcrossing to occur at progressively later epochs as $\Omega_{\Lambda}$
increases. We apply the spherical model to evolve the initial conditions of
Nbody simulated voids and compare the resulting final void profiles. We find
that the model is successful in tracking the evolution of voids with radii
greater than $30 h^{1} \rm Mpc$, implying that void profiles could be used to
constrain dark energy. We find that the initial peculiar velocities of voids
play a significant role in shaping their evolution. Excluding the peculiar
velocity in the evolution model delays the time of shell crossing.

We have derived estimators for the linear growth rate of density fluctuations
using the crosscorrelation function of voids and haloes in redshift space,
both directly and in Fourier form. In linear theory, this crosscorrelation
contains only monopole and quadrupole terms. At scales greater than the void
radius, linear theory is a good match to voids traced out by haloes in Nbody
simulations; smallscale random velocities are unimportant at these radii, only
tending to cause small and often negligible elongation of the redshiftspace
crosscorrelation function near its origin. By extracting the monopole and
quadrupole from the crosscorrelation function, we measure the linear growth
rate without prior knowledge of the void profile or velocity dispersion. We
recover the linear growth parameter $\beta$ to 9% precision from an effective
volume of 3(Gpc/h)^3 using voids with radius greater than 25Mpc/h. Smaller
voids are predominantly subvoids, which may be more sensitive to the random
velocity dispersion; they introduce noise and do not help to improve the
measurement. Adding velocity dispersion as a free parameter allows us to use
information at radii as small as half of the void radius. The precision on
$\beta$ is reduced to approximately 5%. Contrary to the simple redshiftspace
distortion pattern in overdensities, voids show diverse shapes in redshift
space, and can appear either elongated or flattened along the line of sight.
This can be explained by the competing amplitudes of the local density
contrast, plus the radial velocity profile and its gradient, with the latter
two factors being determined by the cumulative density profile of voids. The
distortion pattern is therefore determined solely by the void profile and is
different for voidincloud and voidinvoid. This diversity of redshiftspace
void morphology complicates measurements of the AlcockPaczynski effect using
voids.

We investigate void properties in $f(R)$ models using Nbody simulations,
focusing on their differences from General Relativity (GR) and their
detectability. In the HuSawicki $f(R)$ modified gravity (MG) models, the halo
number density profiles of voids are not distinguishable from GR. In contrast,
the same $f(R)$ voids are more empty of dark matter, and their profiles are
steeper. This can in principle be observed by weak gravitational lensing of
voids, for which the combination of a spectroscopic redshift and a lensing
photometric redshift survey over the same sky is required. Neglecting the
lensing shape noise, the $f(R)$ model parameter amplitudes $f_{R0}=10^{5}$
and $10^{4}$ may be distinguished from GR using the lensing tangential shear
signal around voids by 4 and 8$\sigma$ for a volume of 1~(Gpc/$h$)$^3$. The
lineofsight projection of largescale structure is the main systematics that
limits the significance of this signal for the near future wide angle and deep
lensing surveys. For this reason, it is challenging to distinguish
$f_{R0}=10^{6}$ from GR. We expect that this can be overcome with larger
volume. The halo void abundance being smaller and the steepening of dark matter
void profiles in $f(R)$ models are unique features that can be combined to
break the degeneracy between $f_{R0}$ and $\sigma_8$.

We explore voids in dark matter and halo fields from simulations of
$\Lambda$CDM and HuSawicki $f(R)$ models. In $f(R)$ gravity, dark matter void
abundances are greater than that of general relativity (GR). However, when
using haloes to identify voids, the differences of void abundances become much
smaller, but can still be told apart, in principle, at the 2, 6 and 14 $\sigma$
level for the $f(R)$ model parameter amplitudes of $f_{R0}=10^{6}$,
$10^{5}$ and $10^{4}$. In contrast, the abundance of large voids found using
haloes in $f(R)$ gravity is lower than in GR. The more efficient halo formation
in underdense regions makes $f(R)$ voids less empty of haloes. This counter
intuitive result suggests that voids are not necessarily emptier in $f(R)$ if
one looks at galaxies in voids. Indeed, the halo number density profiles of
voids are not distinguishable from GR. However, the same $f(R)$ voids are more
empty of dark matter. This can in principle be observed by weak gravitational
lensing of voids, for which the combination of a spec$z$ and a photo$z$
survey over the same sky is necessary. For a volume of 1~(Gpc/$h$)$^3$,
neglecting the lensing shape noise, $f_{R0}=10^{5}$ and $10^{4}$ may be
distinguished from GR using the lensing tangential shear signal around voids by
4 and 8$\sigma$. The lineofsight projection of largescale structure is the
main systematics that limits the significance of this signal, limiting the
constraining power for $f_{R0}=10^{6}$. The halo void abundance being
smaller and the steepening of dark matter void profiles in $f(R)$ models are
unique features that can be combined to break the degeneracy between $f_{R0}$
and $\sigma_8$. The outflow of mass from void centers and velocity dispersions
are greater in $f(R)$. Model differences in velocity profiles imply potential
powerful constraints of the model in phase space and in redshift space.

We study the latetime Integrated SachsWolfe (ISW) effect in $f(R)$ gravity
using Nbody simulations. In the $f(R)$ model under study, the linear growth
rate is larger than that in general relativity (GR). This slows down the decay
of the cosmic potential and induces a smaller ISW effect on large scales.
Therefore, the $\dot\Phi$ (time derivative of the potential) power spectrum at
$k<0.1h$/Mpc is suppressed relative to that in GR. In the nonlinear regime,
relatively rapid structure formation in $f(R)$ gravity boosts the nonlinear
ISW effect relative to GR, and the $\dot\Phi$ power spectrum at $k>0.1h$/Mpc is
increased (100$\%$ greater on small scales at $z=0$). We explore the
detectability of the ISW signal via stacking supercluster and supervoids. The
differences in the corresponding ISW cold or hot spots are $\sim 20\%$ for
structures of $\sim 100$Mpc/$h$. Such differences are greater for smaller
structures, but the amplitude of the signal is lower. The high amplitude of ISW
signal detected by Granett et al. can not explained in the $f(R)$ model. We
find relatively big differences between $f(R)$ and GR in the transverse bulk
motion of matter, and discuss its detectability via the relative frequency
shifts of photons from multiple lensed images.

We measure the average temperature decrement on the cosmic microwave
background (CMB) produced by voids selected in the SDSS DR7 spectroscopic
redshift galaxy catalog, spanning redshifts $0<z<0.44$. We find an imprint of
amplitude between 2.6 and 2.9$\mu K$ as viewed through a compensated tophat
filter scaled to the radius of each void; we assess the statistical
significance of the imprint at ~2$\sigma$. We make crucial use of $N$body
simulations to calibrate our analysis. As expected, we find that large voids
produce cold spots on the CMB through the Integrated SachsWolfe (ISW) effect.
However, we also find that small voids in the halo density field produce hot
spots, because they reside in contracting, largerscale overdense regions. This
is an important effect to consider when stacking CMB imprints from voids of
different radius. We have found that the same filter radius that gives the
largest ISW signal in simulations also yields close to the largest detected
signal in the observations. However, although it is low in significance, our
measured signal is much higheramplitude than expected from ISW in the
concordance $\Lambda$CDM universe. The discrepancy is also at the ~2$\sigma$
level. We have demonstrated that our result is robust against the varying of
thresholds over a wide range.

Redshift space distortion (RSD) is a powerful way of measuring the growth of
structure and testing General Relativity, but it is limited by cosmic variance
and the degeneracy between galaxy bias b and the growth rate factor f. The
crosscorrelation of lensing shear with the galaxy density field can in
principle measure b in a manner free from cosmic variance limits, breaking the
fb degeneracy and allowing inference of the matter power spectrum from the
galaxy survey. We analyze the growth constraints from a realistic tomographic
weak lensing photoz survey combined with a spectroscopic galaxy redshift
survey over the same sky area. For sky coverage f_sky=0.5, analysis of the
transverse modes measures b to 23% accuracy per \Delta z=0.1 bin at z<1 when
\sim10 galaxies per square arcmin are measured in the lensing survey and all
halos with M>M_min=10^13 h^{1}M_solar have spectra. For the gravitational
growth parameter parameter \gamma (f=\Omega_m^{\gamma}), combining the lensing
information with RSD analysis of nontransverse modes yields accuracy
\sigma(\gamma)\sim 0.01. Adding lensing information to the RSD survey improves
\sigma(\gamma) by an amount equivalent to a 3 times (10 times) increase in RSD
survey area when the spectroscopic survey extends down to halo mass 10^13.5
(10^14) h^{1} M_solar. We also find that the \sigma(\gamma) of overlapping
surveys is equivalent to that of surveys 1.52 times larger if they are
separated on the sky. This gain is greatest when the spectroscopic mass
threshold is 10^13  10^14 h^{1} M_solar, similar to LRG surveys. The gain of
overlapping surveys is reduced for very deep or very shallow spectroscopic
surveys, but any practical surveys are more powerful when overlapped than when
separated. The gain of overlapped surveys is larger in the case when the
primordial power spectrum normalization is uncertain by >0.5%.

We examine the improvements in constraints on the linear growth factor G and
its derivative f=d ln G / dln a that are available from the combination of a
largescale galaxy redshift survey with a weak gravitational lensing survey of
background sources. In the linear perturbation theory limit, the
biasmodulation method of McDonald & Seljak allows one to distinguish the
realspace galaxy clustering from the peculiar velocity signal in each Fourier
mode. The ratio of lensing signal to galaxy clustering in transverse modes
yields the bias factor b of each galaxy subset (as per Pen 2004), hence
calibrating the conversion from galaxy realspace density to matter density in
every mode. In combination these techniques permit measure of the growth rate f
in each Fourier mode. This yields a measure of the growth rate free of sample
variance, i.e. the uncertainty in f can be reduced without bound by increasing
the number of redshifts within a finite volume. In practice, the gain from the
absence of sample variance is bounded by the limited range of bias modulation
among darkmatter halos. Nonetheless, the addition of background weak lensing
data to a redshift survey increases information on G and f by an amount
equivalent to a 10fold increase in the volume of a standard redshiftspace
distortion measurementif the lensing signal can be measured to subpercent
accuracy. This argues that a combined lensing and redshift survey over a common
lowredshift volume is a more powerful test of general relativity than an
isolated redshift survey over larger volume at high redshift. An example case
is that a survey of ~10^6 redshifts over half the sky in the redshift range
$z=0.5\pm 0.05$ can determine the growth exponent \gamma for the model
$f=\Omega_m^\gamma$ to an accuracy of $\pm 0.015$, using only modes with
k<0.1h/Mpc, but only if a weak lensing survey is conducted in concert.
[Abridged]

We derive the weight function w(M) to apply to darkmatter halos that
minimizes the stochasticity between the weighted halo distribution and its
underlying mass density field. The optimal w(M) depends on the range of masses
being used in the estimator. In Nbody simulations, the Poisson estimator is up
to 15 times noisier than the optimal. Implementation of the optimal weight
yields significantly lower stochasticity than weighting halos by their mass,
bias or equal. Optimal weighting could make cosmological tests based on the
matter power spectrum or crosscorrelations much more powerful and/or
costeffective. A volumelimited measurement of the mass power spectrum at
k=0.2h/Mpc over the entire z<1 universe could ideally be done using only 6
million redshifts of halos with mass M>6\times10^{13}h^{1}M_\odot
(1\times10^{13}) at z=0 (z=1); this is 5 times fewer than the Poisson model
predicts. Using halo occupancy distributions (HOD) we find that
uniformlyweighted catalogs of luminous red galaxies require >3 times more
redshifts than an optimallyweighted halo catalog to reconstruct the mass to
the same accuracy. While the mean HODs of galaxies above a threshold luminosity
are similar to the optimal w(M), the stochasticity of the halo occupation
degrades the mass estimator. Blue or emissionline galaxies are about 100 times
less efficient at reconstructing mass than an optimal weighting scheme. This
suggests an efficient observational approach of identifying and weighting halos
with a deep photoz survey before conducting a spectroscopic survey. The
optimal w(M) and massestimator stochasticity predicted by the standard halo
model for M>10^{12}h^{1}M_\odot are in reasonable agreement with our
measurements, with the important exceptions that the halos must be assumed to
be linearly biased samples of a "halo field" that is distinct from the mass
field. (Abridged)

We present a new method for constructing maps of the secondary temperature
fluctuations imprinted on the cosmic microwave background (CMB) radiation by
photons propagating through the evolving cosmic gravitational potential. Large
cosmological Nbody simulations are used to calculate the complete nonlinear
evolution of the peculiar gravitational potential. Tracing light rays back
through the past lightcone of a chosen observer accurately captures the
temperature perturbations generated by linear (the integrated SachsWolfe or
ISW effect) and nonlinear (the ReesSciama or RS effect) evolution. These
effects give rise to three kinds of nonlinear features in the temperature
maps. (a) In overdense regions, converging flows of matter induce cold spots of
order 100 Mpc in extent which can dominate over the ISW effect at high
redshift, and are surrounded by hot rings. (b) In underdense regions, the RS
effect enhances ISW cold spots which can be surrounded by weak hot rings. (c)
Transverse motions of large lumps of matter produce characteristic dipole
features, consisting of adjacent hot and cold spots separated by a few tens of
Megaparsecs. These nonlinear features are not easily detectable; they modulate
the ISW sky maps at about the 10 percent level. The RS effect causes the
angular power spectrum to deviate from linear theory at l~50 and generates
nonGaussianity, skewing the onepoint distribution function to negative
values. Cold spots of similar angular size, but much smaller amplitude than the
CMB cold spot reported by Cruz et al. are produced. Joint analysis of our maps
and the corresponding galaxy distribution may enable techniques to be developed
to detect these nonlinear, nonGaussian features. Our maps are available at
http://astro.dur.ac.uk/~cai/ISW

In a universe with a cosmological constant, the largescale gravitational
potential varies in time and this is, in principle, observable. Using an Nbody
simulation of a $\Lambda$CDM universe, we show that linear theory is not
sufficiently accurate to predict the power spectrum of the time derivative,
$\dot{\Phi}$, needed to compute the imprint of largescale structure on the
cosmic microwave background (CMB). The linear part of the $\dot{\Phi}$ power
spectrum (the integrated SachsWolfe effect or ISW) drops quickly as the
relative importance of $\Omega_{\Lambda}$ diminishes at high redshift, while
the nonlinear part (the ReesSciama effect or RS) evolves more slowly with
redshift. Therefore, the deviation of the total power spectrum from linear
theory occurs at larger scales at higher redshifts. The deviation occurs at
$k\sim 0.1 $ $h$ Mpc$^{1}$ at $z=0$. The crosscorrelation power spectrum of
the density $\delta$ with $\dot{\Phi}$ behaves differently to the power
spectrum of $\dot{\Phi}$. Firstly, the deviation from linear theory occurs at
smaller scales ($k\sim 1 $ $h$ Mpc$^{1}$ at $z=0$). Secondly, the correlation
becomes negative when the nonlinear effect dominates. For the
crosscorrelation power spectrum of galaxy samples with the CMB, the nonlinear
effect becomes significant at $l\sim 500$ and rapidly makes the cross power
spectrum negative. For high redshift samples, the crosscorrelation is expected
to be suppressed by $510%$ on arcminute scales. The RS effect makes a
negligible contribution to the largescale ISW crosscorrelation measurement.
However, on arcminute scales it will contaminate the expected
crosscorrelation signal induced by the SunyaevZel'dovich effect.

We describe a method for constructing mock galaxy catalogues which are well
suited for use in conjunction with large photometric surveys. We use the
semianalytic galaxy formation model of Bower et al. implemented in the
Millennium simulation. We apply our method to the specific case of the surveys
soon to commence with PS1, the first of 4 telescopes planned for the PanSTARRS
system. PS1 has 5 photometric bands (grizy), and will carry out an allsky 3pi
survey and a medium deep survey (MDS) over 84 sq.deg. We calculate the expected
magnitude limits for extended sources in the two surveys. We find that, after 3
years, the 3pi survey will have detected over 10^8 galaxies in all 5 bands, 10
million of which will lie at redshift z>0.9, while the MDS will have detected
over 10^7 galaxies with 0.5 million lying at z>2. These numbers at least double
if detection in the shallowest band, y is not required. We then evaluate the
accuracy of photometric redshifts estimated using an offtheshelf photoz
code. With the grizy bands alone it is possible to achieve an accuracy in the
3pi survey of Delta z/(1+z)~0.06 for 0.25<z<0.8, which could be reduced by
about 15% using near infrared photometry from the UKIDDS survey, but would
increase by about 25% for the deeper sample without the y band photometry. For
the MDS an accuracy of Delta z/(1+z)~0.05 is achievable for 0.02<z<1.5 using
grizy. A dramatic improvement in accuracy is possible by selecting only red
galaxies. In this case, Delta z/(1+z)~0.020.04 is achievable for ~100 million
galaxies at 0.4<z<1.1 in the 3pi survey and for 30 million galaxies in the MDS
at 0.4<z<2. We investigate the effect of using photoz in the estimate of the
baryonic acoustic oscillation scale. We find that PS1 will achieve a similar
accuracy in this estimate as a spectroscopic survey of 20 million galaxies.

The power spectrum of the twodegree Field Galaxy Redshift Survey (2dFGRS)
sample is estimated with the discrete wavelet transform (DWT) method. The DWT
power spectra within $0.04 <k< 2.3 h$Mpc$^{1}$ are measured for three
volumelimited samples defined in connective absolute magnitude bins $19 \sim
18$, $20 \sim 19$ and $21 \sim 20$. We show that the DWT power spectrum
can effectively distinguish $\Lambda$CDM models of $\sigma_8=0.84$ and
$\sigma_8=0.74$. We adopt maximum likelihood method to perform threeparameter
fitting with bias parameter $b$, pairwise velocity dispersion $\sigma_{pv}$ and
redshift distortion parameter $\beta=\Omega_m^{0.6}/b$ to the measured DWT
power spectrum. Fitting results denotes that in a $\sigma_8=0.84$ universe the
best fitted $\Omega_m$ given by the three samples are consistent in the range
$0.28 \sim 0.36$, and the best fitted $\sigma_{pv}$ are $398^{+35}_{27}$,
$475^{+37}_{29}$ and $550 \pm 20$km/s for the three samples, respectively.
However in the model of $\sigma_8=0.74$, our three samples give very different
values of $\Omega_m$. We repeat the fitting by using empirical formula of
redshift distortion. The result of the model of low $\sigma_8$ is still poor,
especially, one of the best value $\sigma_{pv}$ is as large as $10^3$km/s. The
power spectrum of 2dFGRS seems in disfavor of models with low amplitude of
density fluctuations.