
PanSTARRS1 has carried out a set of distinct synoptic imaging sky surveys
including the $3\pi$ Steradian Survey and the Medium Deep Survey in 5 bands
($grizy_{P1}$). The mean 5$\sigma$ point source limiting sensitivities in the
stacked 3$\pi$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3,
21.4) respectively. The upper bound on the systematic uncertainty in the
photometric calibration across the sky is 712 millimag depending on the
bandpass. The systematic uncertainty of the astrometric calibration using the
Gaia frame comes from a comparison of the results with Gaia: the standard
deviation of the mean and median residuals ($ \Delta ra, \Delta dec $) are
(2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The PanSTARRS
system and the design of the PS1 surveys are described and an overview of the
resulting image and catalog data products and their basic characteristics are
described together with a summary of important results. The images, reduced
data products, and derived data products from the PanSTARRS1 surveys are
available to the community from the Mikulski Archive for Space Telescopes
(MAST) at STScI.

This work is on the Physics of the B Factories. Part A of this book contains
a brief description of the SLAC and KEK B Factories as well as their detectors,
BaBar and Belle, and data taking related issues. Part B discusses tools and
methods used by the experiments in order to obtain results. The results
themselves can be found in Part C.
Please note that version 3 on the archive is the auxiliary version of the
Physics of the B Factories book. This uses the notation alpha, beta, gamma for
the angles of the Unitarity Triangle. The nominal version uses the notation
phi_1, phi_2 and phi_3. Please cite this work as Eur. Phys. J. C74 (2014) 3026.

We report results of a study of Planck SunyaevZel'dovich effect (SZE)
selected galaxy cluster candidates using the Panoramic Survey Telescope & Rapid
Response System (PanSTARRS) imaging data. We first examine 150 Planck
confirmed galaxy clusters with spectroscopic redshifts to test our algorithm
for identifying optical counterparts and measuring their redshifts; our
redshifts have a typical accuracy of $\sigma_{z/(1+z)} \sim 0.022$ for this
sample. Using 60 random sky locations, we estimate that our chance of
contamination through a random superposition is ~ 3 per cent. We then examine
an additional 237 Planck galaxy cluster candidates that have no redshift in the
source catalogue. Of these 237 unconfirmed cluster candidates we are able to
confirm 60 galaxy clusters and measure their redshifts. A further 83 candidates
are so heavily contaminated by stars due to their location near the Galactic
plane that we do not attempt to identify counterparts. For the remaining 94
candidates we find no optical counterpart but use the depth of the PanSTARRS1
data to estimate a redshift lower limit $z_{\text{lim}(10^{15})}$ beyond which
we would not have expected to detect enough galaxies for confirmation. Scaling
from the already published Planck sample, we expect that $\sim$12 of these
unconfirmed candidates may be real clusters.

The PanSTARRS1 survey is acquiring multiepoch imaging in 5 bands (grizy)
over the entire sky north of declination 30deg (the $3\pi$ survey). In July
2011 a test area of about 70 sq.deg. was observed to the expected final depth
of the main survey. In this, the first of a series of papers targetting the
galaxy count and clustering properties of the combined multiepoch test area
data, we present a detailed investigation into the depth of the survey and the
reliability of the PanSTARRS1 analysis software. We show that the PanSTARRS1
reduction software can recover the properties of fake sources, and show good
agreement between the magnitudes measured by PanSTARRS1 and those from Sloan
Digital Sky Survey. We also examine the number of false detections apparent in
the PanSTARRS1 data. Our comparisons show that the test area survey is
somewhat deeper than the Sloan Digital Sky Survey in all bands, and, in
particular, the z band approaches the depth of the stacked Sloan Stripe 82
data.

We apply a recently developed scaling technique to the MillenniumXXL, one of
the largest cosmological Nbody simulations carried out to date 3x10^11
particles within a cube of volume ~70Gpc^3). This allows us to investigate the
cosmological parameter dependence of the mass and evolution of haloes in the
extreme highmass tail of the z=6 distribution. We assume these objects to be
likely hosts for the population of rare but ultraluminous highredshift quasars
discovered by the Sloan Digital Sky Survey. Haloes with a similar abundance to
these quasars have a median mass of 9x10^12 Msun in the currently preferred
cosmology, but do not evolve into equally extreme objects at z=0. Rather, their
descendants span the full range conventionally assigned to presentday
clusters, 6x10^13 to 2.5x10^15 Msun for this same cosmology. The masses both at
z=6 and at z=0 shift up or down by factors exceeding two if cosmological
parameters are pushed to the boundaries of the range discussed in published
interpretations of data from the WMAP satellite. The main factor determining
the future growth of a highmass z=6 halo is the mean overdensity of its
environment on scales of 7 to 14 Mpc, and descendant masses can be predicted 6
to 8 times more accurately if this density is known than if it is not. All
these features are not unique to extreme highz haloes, but are generic to
hierarchical growth. Finally, we find that extreme haloes at z=6 typically
acquired about half of their total mass in the preceding 100 Myr, implying very
large recent accretion rates which may be related to the large black hole
masses and high luminosities of the SDSS quasars.

The estimated stellar masses of galaxies are widely used to characterize how
the galaxy population evolves over cosmic time. If stellar masses can be
estimated in a robust manner, free from any bias, global diagnostics such as
the stellar mass function can be used to constrain the physics of galaxy
formation. We explore how galaxy stellar masses, estimated by fitting
broadband spectral energy distributions (SEDs) with stellar population models,
can be biased as a result of commonly adopted assumptions for the
starformation and chemical enrichment histories, recycled fractions and dust
attenuation curves of galaxies. We apply the observational technique of
broadband SED fitting to model galaxy SEDs calculated by the theoretical
galaxy formation model GALFORM, isolating the effect of each of these
assumptions. We find that, averaged over the entire galaxy population, the
common assumption of exponentially declining starformation histories does not
adversely affect stellar mass estimation. We show that fixing the metallicity
in SED fitting or using sparsely sampled metallicity grids can introduce mass
dependent systematics into stellar mass estimates. We find that the common
assumption of a stardust geometry corresponding to a uniform foreground dust
screen can cause the stellar masses of dusty model galaxies to be significantly
underestimated. Finally, we show that stellar mass functions recovered by
applying SED fitting to model galaxies at high redshift can differ
significantly in both shape and normalization from the intrinsic mass functions
predicted by a given model. Given these differences, our methodology of using
stellar masses estimated from model galaxy SEDs offers a new, selfconsistent
way to compare model predictions with observations.

The formation of galaxies is regulated by a balance between the supply of gas
and the rate at which it is ejected. Traditional explanations of gas ejection
equate the energy required to escape the galaxy or host halo to an estimate for
the energy yield from supernovae. This yield is usually assumed to be a
constant fraction of the total available from the supernova, or is derived from
the assumption of a consistent momentum yield. By applying these ideas in the
context of a cold dark matter cosmogony, we derive a 1storder analytic
connection between these working assumptions and the expected relationship
between baryon content and galaxy circular velocity, and find that these quick
predictions straddle recent observational estimates. To examine the premises
behind these theories in more detail, we then explore their applicability to a
set of gasdynamical simulations of idealised galaxies. We show that different
premises dominate to differing degrees in the simulated outflow, depending on
the mass of the system and the resolution with which it is simulated. Using
this study to anticipate the emergent behaviour at arbitrarily high resolution,
we motivate more comprehensive analytic model which allows for the range of
velocities with which the gas may exit the system, and incorporates both
momentum and energybased constraints on the outflow. Using a trial exit
velocity distribution, this is shown to be compatible with the observed baryon
fractions in intermediatemass systems, but implies that current estimates for
lowmass systems can not be solely accounted for by supernova winds under
commonlyheld assumptions.

We characterize substructure in the simulated stellar halos of Cooper et al.
(2010) which were formed by the disruption of satellite galaxies within the
cosmological Nbody simulations of galactic halos of the Aquarius Project.
These stellar halos exhibit a wealth of tidal features: broad overdensities and
very narrow faint streams akin to those observed around the Milky Way. The
substructures are distributed anisotropically on the sky, a characteristic that
should become apparent in the next generation of photometric surveys. The
normalized RMS of the density of stars on the sky appears to be systematically
larger for our halos compared to the value estimated for the Milky Way from
main sequence turnoff stars in the Sloan Digital Sky Survey. We show that this
is likely to be due in part to contamination by faint QSOs and redder main
sequence stars, and might suggest that ~10% of the Milky Way halo stars have
formed insitu.

We describe a correlation function statistic that quantifies the amount of
spatial and kinematic substructure in the stellar halo. We test this statistic
using model stellar halo realizations constructed from the Aquarius suite of
six highresolution Nbody simulations in combination with the Galform
semianalytic galaxy formation model. These simulations show considerable
scatter in the properties of stellar haloes. We find that our statistic can
distinguish between these plausible alternatives for the global structure of
the Milky Way stellar halo. We compare with observational data and show that
pencil beam surveys of ~100 tracer stars (such as the Spaghetti Survey) are not
sufficient to constrain the degree of structure in the Milky Way halo with this
statistic. Larger area surveys with >1000 tracer stars (such as BHB stars in
the Sloan Digital Sky Survey) provide much tighter constraints on comparisons
between models and data. In our simulations, we find examples of haloes with
spatial and kinematic substructure consistent with the available Milky Way
data.

We use a coupled model of the formation and evolution of galaxies and black
holes (BH) to study the evolution of active galactic nuclei (AGN) in a cold
dark matter universe. The model predicts the BH mass, spin and mass accretion
history. BH mass grows via accretion triggered by discs becoming dynamically
unstable or galaxy mergers (called the starburst mode) and accretion from
quasihydrostatic hot gas haloes (called the hothalo mode). By taking into
account AGN obscuration, we obtain a very good fit to the observed luminosity
functions (LF) of AGN (optical, soft and hard Xray, and bolometric) for a wide
range of redshifts (0<z<6). The model predicts a hierarchical build up of BH
mass, with the typical mass of actively growing BHs increasing with decreasing
redshift. Remarkably, despite this, we find downsizing in the AGN population,
in terms of the differential growth with redshift of the space density of faint
and bright AGN. This arises naturally from the interplay between the starburst
and hothalo accretion modes. The faint end of the LF is dominated by massive
BHs experiencing quiescent accretion via a thick disc, primarily during the
hothalo mode. The bright end of the LF, on the other hand, is dominated by AGN
which host BHs accreting close to or in excess of the Eddington limit during
the starburst mode. The model predicts that the comoving space density of AGN
peaks at z~3, similar to the star formation history. However, when taking into
account obscuration, the space density of faint AGN peaks at lower redshift
(z<2) than that of bright AGN (z~23). This implies that the cosmic evolution
of AGN is shaped in part by obscuration.

We track the coevolution of supermassive black holes (SMBHs) and their host
galaxies. The calculation is embedded in the GALFORM semianalytical model
which simulates the formation and evolution of galaxies in a cold dark matter
(CDM) universe. During the evolution of the host galaxy, hot and cold gas are
added to the SMBH by flows triggered by halo gas cooling, disc instabilities
and galaxy mergers. This builds up the mass and spin of the BH, and the
resulting accretion power regulates the gas cooling and subsequent star
formation. The accretion flow is assumed to form a geometrically thin cool disc
when the accretion rate exceeds 0.01\dot{M}_Edd, and a geometrically thick,
radiatively inefficient hot flow when the accretion rate falls below this
value. The resulting quasar optical luminosity function matches observations
very well, and the mass of the SMBH correlates with the mass of the galaxy
bulge as observed. The BH spin distribution depends strongly on whether the gas
in any given accretion episode remains in the same plane (prolonged accretion)
or whether, due to selfgravity, it fragments into multiple, randomly aligned
accretion episodes (chaotic accretion). In the chaotic accretion model there is
a clear correlation of spin with SMBH mass. Massive BHs (M>5\times10^8\Msun)
are hosted by giant elliptical galaxies and are rapidly spinning, while lower
mass BHs are hosted in spiral galaxies and have much lower spin. Using the
BlandfordZnajek mechanism for jet production to calculate the jet power, our
model is able to reproduce the radio loudness of radio galaxies, LINERS and
Seyferts. This is the first confirmation that a CDM galaxy formation model can
reproduce the observed phenomenology of AGN.

We present six simulations of galactic stellar haloes formed by the tidal
disruption of accreted dwarf galaxies in a fully cosmological setting. Our
model is based on the Aquarius project, a suite of high resolution Nbody
simulations of individual dark matter haloes. We tag subsets of particles in
these simulations with stellar populations predicted by the Galform
semianalytic model. Our method selfconsistently tracks the dynamical
evolution and disruption of satellites from high redshift. The luminosity
function and structural properties of surviving satellites, which agree well
with observations, suggest that this technique is appropriate. We find that
accreted stellar haloes are assembled between 1<z<7 from less than 5
significant progenitors. These progenitors are old, metalrich satellites with
stellar masses similar to the brightest Milky Way dwarf spheroidals. In
contrast to previous stellar halo simulations, we find that several of these
major contributors survive as selfbound systems to the present day. Both the
number of these significant progenitors and their infall times are inherently
stochastic. This results in great diversity among our stellar haloes, which
amplifies small differences between the formation histories of their dark halo
hosts. The masses and density/surfacebrightness profiles of the stellar haloes
are consistent with expectations from the Milky Way and M31. Each halo has a
complex structure, consisting of wellmixed components, tidal streams, shells
and other subcomponents. This structure is not adequately described by smooth
models. We find one example of an accreted thick disk. Contrasts in age and
metallicity between halo stars and those in surviving satellites are in broad
agreement with recent observations. [Abridged]

We derive constraints on cosmological parameters using the power spectrum of
galaxy clustering measured from the final twodegree field galaxy redshift
survey (2dFGRS) and a compilation of measurements of the temperature power
spectrum and temperaturepolarization crosscorrelation of the cosmic microwave
background radiation. We analyse a range of parameter sets and priors, allowing
for massive neutrinos, curvature, tensors and general dark energy models. In
all cases, the combination of datasets tightens the constraints, with the most
dramatic improvements found for the density of dark matter and the
energydensity of dark energy. If we assume a flat universe, we find a matter
density parameter of $\Omega_{\rm m}=0.237 \pm 0.020$, a baryon density
parameter of $\Omega_{\rm b} = 0.041 \pm 0.002$, a Hubble constant of
$H_{0}=74\pm2 {\rm kms}^{1}{\rm Mpc}^{1}$, a linear theory matter fluctuation
amplitude of $\sigma_{8}=0.77\pm0.05$ and a scalar spectral index of $n_{\rm
s}=0.954 \pm 0.023$ (all errors show the 68% interval). Our estimate of $n_{\rm
s}$ is only marginally consistent with the scale invariant value $n_{\rm s}=1$;
this spectrum is formally excluded at the 95% confidence level. However, the
detection of a tilt in the spectrum is sensitive to the choice of parameter
space. If we allow the equation of state of the dark energy to float, we find
$w_{\rm DE}= 0.85_{0.17}^{+0.18}$, consistent with a cosmological constant.
We also place new limits on the mass fraction of massive neutrinos: $f_{\nu} <
0.105$ at the 95% level, corresponding to $\sum m_{\nu} < 1.2$ eV.

Galactic bulges are known to harbour central black holes whose mass is
tightly correlated with the stellar mass and velocity dispersion of the bulge.
In a hierarchical universe, mergers of subgalactic units are accompanied by the
amalgamation of bulges and the likely coalescence of galactocentric black
holes. In these mergers, the beaming of gravitational radiation during the
plunge phase of the black hole collision can impart a linear momentum kick or
``gravitational recoil'' to the remnant. If large enough, this kick will eject
the remnant from the galaxy and populate intergalactic space with wandering
black holes. Using a semianalytic model of galaxy formation, we investigate
the effect of black hole ejections on the scatter in the relation between black
hole and bulge mass. We find that although not the dominant source of the
measured scatter, they do make a significant contribution and may be used to
set a constraint, v_kick<500 km/s, on the typical kick velocity, in agreement
with values found from general relativistic calculations. Even for the more
modest kick velocities implied by these calculations, we find that a
substantial number of central black holes are ejected from the progenitors of
present day galaxies, giving rise to a population of wandering intrahalo and
intergalactic black holes whose distribution we investigate in highresolution
Nbody simulations of MilkWay mass halos. We find that intergalactic black
holes make up only ~23% of the total galactic black hole mass but, within a
halo, wandering black holes can contribute up to about half of the total black
hole mass orbiting the central galaxy. Intrahalo black holes offer a natural
explanation for the compact Xray sources often seen near the centres of
galaxies and for the hyperluminous noncentral Xray source in M82.

We present a power spectrum analysis of the final 2dF Galaxy Redshift Survey,
employing a direct Fourier method. The sample used comprises 221,414 galaxies
with measured redshifts. We investigate in detail the modelling of the sample
selection. A new angular mask is derived, based on revisions to the photometric
calibration. The redshift selection function is determined by dividing the
survey according to restframe colour, and deducing a selfconsistent treatment
of kcorrections and evolution for each population. The covariance matrix for
the powerspectrum estimates is determined using two different approaches to
the construction of mock surveys which are used to demonstrate that the input
cosmological model can be correctly recovered. We are confident that the 2dFGRS
power spectrum can be used to infer the matter content of the universe. On
large scales, our estimated power spectrum shows evidence for the `baryon
oscillations' that are predicted in CDM models. Fitting to a CDM model,
assuming a primordial $n_{s}=1$ spectrum, $h=0.72$ and negligible neutrino
mass, the preferred parameters are $\Omega_{M} h = 0.168 \pm 0.016$ and a
baryon fraction $\Omega_{b} /\Omega_{M} = 0.185\pm0.046$ (1$\sigma$ errors).
The value of $\Omega_{M} h$ is $1\sigma$ lower than the $0.20 \pm 0.03$ in our
2001 analysis of the partially complete 2dFGRS. This shift is largely due to
the signal from the newlysampled regions of space, rather than the refinements
in the treatment of observational selection. This analysis therefore implies a
density significantly below the standard $\Omega_{M} =0.3$: in combination with
CMB data from WMAP, we infer $\Omega_{M} =0.231\pm 0.021$. (Abridged.)

We use the 2dF Galaxy Redshift Survey to measure the dependence of the
bJband galaxy luminosity function on largescale environment, defined by
density contrast in spheres of radius 8h1Mpc, and on spectral type, determined
from principal component analysis. We find that the galaxy populations at both
extremes of density differ significantly from that at the mean density. The
population in voids is dominated by late types and shows, relative to the mean,
a deficit of galaxies that becomes increasingly pronounced at magnitudes
brighter than M_bJ5log10h <18.5. In contrast, cluster regions have a relative
excess of very bright earlytype galaxies with M_bJ5log10h < 21. Differences
in the mid to faintend population between environments are significant: at
M_bJ5log10h=18 early and latetype cluster galaxies show comparable
abundances, whereas in voids the late types dominate by almost an order of
magnitude. We find that the luminosity functions measured in all density
environments, from voids to clusters, can be approximated by Schechter
functions with parameters that vary smoothly with local density, but in a
fashion which differs strikingly for early and latetype galaxies. These
observed variations, combined with our finding that the faintend slope of the
overall luminosity function depends at most weakly on density environment, may
prove to be a significant challenge for models of galaxy formation.

We present predictions for the abundance of submm galaxies (SMGs) and
Lymanbreak galaxies (LBGs) in the $\Lambda$CDM cosmology. A key feature of our
model is the selfconsistent calculation of the absorption and emission of
radiation by dust. The new model successfully matches the LBG luminosity
function, as well reproducing the properties of the local galaxy population in
the optical and IR. The model can also explain the observed galaxy number
counts at $850\mum$, but only if we assume a topheavy IMF for the stars formed
in bursts. The predicted redshift distribution of SMGs depends relatively
little on their flux over the range 1$10\mjy$, with a median value of
$z\approx 2.0$ at a flux of $5\mjy$, in very good agreement with the recent
measurement by Chapman et al The counts of SMGs are predicted to be dominated
by ongoing starbursts. However, in the model these bursts are responsible for
making only a few per cent of the stellar mass locked up in massive ellipticals
at the present day.

The 2dFGRS is used in conjunction with the 2MASSXSC to study the nearIR
light and stellar mass content of the local Universe. Mock galaxy catalogues,
constructed from cosmological Nbody simulations and semianalytical galaxy
formation models, are used to gauge the accuracy with which quantities can be
recovered. The mean luminosity densities of the Universe are found to be
rho_J=(3.57+/0.11)*10^8 h Lsol/Mpc^3 and rho_KS=(7.04+/0.23)*10^8 h
Lsol/Mpc^3 (statistical uncertainty). Using the 2PIGG catalogue, the group
dynamical masstolight ratio in the K_S band is found to increase by a factor
of ~3 when going from groups with total bJband luminosities of 3*10^10 h^2
Lsol to rich clusters, which have typical values of Upsilon_K~80 h Upsilon_sol.
Taking into account the bias introduced by uncertainties in estimating galaxy
stellar masses from luminosities, a value of Omega_* h=(0.99+/0.03)*10^3 is
measured, assuming that a Kennicutt stellar IMF is applicable to all galaxies.
The 2PIGGs are then used to study the distribution of the stellar content of
the local Universe. The three main conclusions are: (1) a slowly rising stellar
M/L_KS is found with the clusters having the largest value of ~0.6 Upsilon_sol,
(2) in contrast, the fraction of mass in stars decreases with increasing group
size, reaching ~5*10^3 h for the rich clusters, and (3) in answer to the
question posed in the title, most stellar mass is contained in Local
Groupsized objects (M~2*10^12 h^1 Msol) with only ~2% in clusters with
M>5*10^14 h^1 Msol.

It is well known that the clustering of galaxies depends on galaxy type.Such
relative bias complicates the inference of cosmological parameters from galaxy
redshift surveys, and is a challenge to theories of galaxy formation and
evolution. In this paper we perform a joint countsincells analysis on
galaxies in the 2dF Galaxy Redshift Survey, classified by both colour and
spectral type, eta, as early or late type galaxies. We fit three different
models of relative bias to the joint probability distribution of the cell
counts, assuming Poisson sampling of the galaxy density field. We investigate
the nonlinearity and stochasticity of the relative bias, with cubical cells of
side 10Mpc \leq L \leq 45Mpc (h=0.7). Exact linear bias is ruled out with high
significance on all scales. Power law bias gives a better fit, but likelihood
ratios prefer a bivariate lognormal distribution, with a nonzero
`stochasticity'  i.e. scatter that may result from physical effects on galaxy
formation other than those from the local density field. Using this model, we
measure a correlation coefficient in logdensity space (r_LN) of 0.958 for
cells of length L=10Mpc, increasing to 0.970 by L=45Mpc. This corresponds to a
stochasticity sigma_b/bhat of 0.44\pm0.02 and 0.27\pm0.05 respectively. For
smaller cells, the Poisson sampled lognormal distribution presents an
increasingly poor fit to the data, especially with regard to the fraction of
completely empty cells. We compare these trends with the predictions of
semianalytic galaxy formation models: these match the data well in terms of
overall level of stochasticity, variation with scale, and fraction of empty
cells.

We measure moments of the galaxy count probability distribution function in
the twodegree field galaxy redshift survey (2dFGRS). The survey is divided
into volume limited subsamples in order to examine the dependence of the higher
order clustering on galaxy luminosity. We demonstrate the hierarchical scaling
of the averaged ppoint galaxy correlation functions, xibar_p, up to p=6. The
hierarchical amplitudes, S_p = xibar_p/xibar_2^{p1}, are approximately
independent of the cell radius used to smooth the galaxy distribution on small
to medium scales. On larger scales we find the higher order moments can be
strongly affected by the presence of rare, massive superstructures in the
galaxy distribution. The skewness S_3 has a weak dependence on luminosity,
approximated by a linear dependence on log luminosity. We discuss the
implications of our results for simple models of linear and nonlinear bias
that relate the galaxy distribution to the underlying mass.

We measure the redshift space reduced void probability function (VPF) for
2dFGRS volume limited galaxy samples covering the absolute magnitude range
M_bJ5logh=18 to 22. Theoretically, the VPF connects the distribution of
voids to the moments of galaxy clustering of all orders, and can be used to
discriminate clustering models in the weakly nonlinear regime. The reduced VPF
measured from the 2dFGRS is in excellent agreement with the paradigm of
hierarchical scaling of the galaxy clustering moments. The accuracy of our
measurement is such that we can rule out, at a very high significance, popular
models for galaxy clustering, including the lognormal distribution. We
demonstrate that the negative binomial model gives a very good approximation to
the 2dFGRS data over a wide range of scales, out to at least 20h1Mpc.
Conversely, the reduced VPF for dark matter in a LambdaCDM universe does appear
to be lognormal on small scales but deviates significantly beyond \approx
4h1Mpc. We find little dependence of the 2dFGRS reduced VPF on galaxy
luminosity. Our results hold independently in both the north and south Galactic
pole survey regions.

We measure the clustering of galaxy groups in the 2dFGRS PercolationInferred
Galaxy Group (2PIGG) catalogue. The 2PIGG sample has 29,000 groups with at
least two members. The clustering amplitude of the full 2PIGG catalogue is
weaker than that of 2dFGRS galaxies, in agreement with theoretical predictions.
We have subdivided the 2PIGG catalogue into samples that span a factor of 25 in
median total luminosity. Our correlation function measurements span an
unprecedented range of clustering strengths, connecting the regimes probed by
groups fainter than L* galaxies and rich clusters. There is a steady increase
in clustering strength with group luminosity; the most luminous groups are ten
times more strongly clustered than the full 2PIGG catalogue. We demonstrate
that the 2PIGG results are in very good agreement with the clustering of groups
expected in the LCDM model.

We use the twodegree field Galaxy Redshift Survey (2dFGRS) to test the
hierarchical scaling hypothesis: namely, that the $p$point galaxy correlation
functions can be written in terms of the two point correlation function or
variance. This scaling is expected if an initially Gaussian distribution of
density fluctuations evolves under the action of gravitational instability. We
measure the volume averaged $p$point correlation functions using a counts in
cells technique applied to a volume limited sample of 44,931 $L_*$ galaxies. We
demonstrate that $L_{*}$ galaxies display hierarchical clustering up to order
$p=6$ in redshift space. The variance measured for $L_{*}$ galaxies is in
excellent agreement with the predictions from a $\Lambda$cold dark matter
Nbody simulation. This applies to all cell radii considered,
$0.3<(R/h^{1}{\rm Mpc})<30$. However, the higher order correlation functions
of $L_*$ galaxies have a significantly smaller amplitude than is predicted for
the dark matter for $R<10h^{1}$Mpc. This disagreement implies that a
nonlinear bias exists between the dark matter and $L_*$ galaxies on these
scales. We also show that the presence of two rare, massive superclusters in
the 2dFGRS has an impact on the higherorder clustering moments measured on
large scales.

We present an analysis of the relative bias between early and latetype
galaxies in the Twodegree Field Galaxy Redshift Survey (2dFGRS). Our analysis
examines the joint counts in cells between early and latetype galaxies, using
approximately cubical cells with sides ranging from 7h^{1}Mpc to 42h^{1}Mpc.
We measure the variance of the counts in cells using the method of Efstathiou
et al. (1990), which we find requires a correction for a finite volume effect.
We fit lognormal models to the onepoint density distribution and develop
methods of dealing with biases in the recovered variances resulting from this
technique. We directly fit deterministic models for the joint density
distribution function, f(delta_E,delta_L), to the joint counts in cells using a
maximum likelihood technique. Our results are consistent with a scale invariant
relative bias factor on all scales studied. Linear bias is ruled out on scales
less than l=28h^{1}Mpc. A powerlaw bias model is a significantly better fit
to the data on all but the largest scales studied; the relative goodness of fit
of this model as compared to that of the linear bias model suggests that any
nonlinearity is negligible for l>~40h^{1}Mpc, consistent with the expectation
from theory that the bias should become linear on large scales. (abridged)

The 2dFGRS PercolationInferred Galaxy Group (2PIGG) catalogue of ~29000
objects is used to study the luminous content of galaxy systems of various
sizes. Mock galaxy catalogues constructed from cosmological simulations are
used to gauge the accuracy with which intrinsic group properties can be
recovered. A Schechter function is found to be a reasonable fit to the galaxy
luminosity functions in groups of different mass in the real data. The
characteristic luminosity L* is larger for more massive groups. However, the
mock data show that the shape of the recovered luminosity function is expected
to differ from the true shape, and this must be allowed for when interpreting
the data. The variation of halo masstolight ratio with group size is studied
in both these wavebands. A robust trend of increasing M/L with increasing group
luminosity is found in the 2PIGG data. From groups with L_bj=10^{10}Lsol to
those 100 times more luminous, the typical bjband M/L increases by a factor of
5, whereas the rfband M/L grows by a factor of 3.5. These trends agree well
with the simulations, which also predict a minimum M/L on a scale corresponding
to the Local Group. Our data indicate that if such a minimum exists, then it
must occur at L<~10^{10}Lsol, below the range accurately probed by the 2PIGG
catalogue. According to the mock data, the bj M/Ls of the largest groups are
expected to be approximately 1.1 times the global value. Assuming that this
correction applies to the real data yields an estimate of Omega_m=0.26+/0.03
(statistical).