A non-parametric higher-order Jeans analysis method, GravSphere, was recently used to constrain the density profiles of Local Group dwarf galaxies, classifying them into those that are more likely to have an inner dark matter cusp and those that are likely to have a core (Read et al.). In this work we test this method using 31 simulated galaxies, comparable to Fornax, selected from the APOSTLE suite of cosmological hydrodynamics simulations, which include CDM and Self-Interacting Dark Matter (SIDM) cosmologies. We find that the mass profiles of the simulated dwarfs are often, but not always, well recovered by GravSphere. The less successful cases may be identified using a chi-squared diagnostic. Although the uncertainties are large in the inner regions, the inferred mass profiles are unbiased and exhibit smaller scatter than comparable Jeans methods. We find that GravSphere recovers the density profiles of simulated dwarfs below the half-light radius and down to the resolution limit of our simulations with better than 10+-30 per cent accuracy, making it a promising Jeans-based approach for modelling dark matter distributions in dwarf galaxies.

A number of Local Group dwarf galaxies are known to have two spatially segregated stellar metallicity populations, a centrally concentrated metal-rich population and a more extended metal-poor population. In this work we discuss mechanisms that lead to the formation of two spatially segregated metallicity populations. Using a set of high-resolution hydrodynamical simulations of Local Group-like environments, we select a sample of satellite and field galaxies, spanning the stellar mass range $10^6-10^9$M$_{\odot}$, that exhibit bimodality in their metallicity distributions. Among those, we identify a subsample with a strong spatial segregation in the two populations. We find three distinct mechanisms for their formation. In field dwarfs and in a small fraction of satellites, a merger causes the metal-poor stars to migrate to larger radii and encourages the available gas to sink to the centre of the dwarf. Most of the gas is subsequently blown out of the halo through star formation feedback, but the remaining gas is consumed in the formation of a metal-rich population. In the exclusive case of satellites that have retained some of their gas at infall, it is the compression of this gas by ram pressure near pericentre that triggers the formation of metal-rich stars, whilst simultaneously preventing star formation at larger radii through stripping. Additionally, in a small number of field and satellite dwarfs, interactions with gaseous filaments and other galaxies can result in the formation of a metal-rich population. Regardless of the formation mechanism, a history of mergers typically enhances the spatial segregation.

We present recalibrations of the GALFORM semi-analytical model of galaxy formation in a new N-body simulation with the Planck cosmology. The Planck Millennium simulation uses more than 128 billion particles to resolve the matter distribution in a cube of $800$ Mpc on a side, which contains more than 77 million dark matter haloes with mass greater than $2.12 \times 10^{9} h^{-1} {\rm M_{\odot}}$ at the present day. Only minor changes to a very small number of model parameters are required in the recalibration. We present predictions for the atomic hydrogen content (HI) of dark matter halos, which is a key input into the calculation of the HI intensity mapping signal expected from the large-scale structure of the Universe. We find that the HI mass $-$ halo mass relation displays a clear break at the halo mass above which AGN heating suppresses gas cooling, $\approx 3 \times 10^{11} h^{-1} M_{\rm \odot}$. Below this halo mass, the HI content of haloes is dominated by the central galaxy; above this mass it is the combined HI content of satellites that prevails. We find that the HI mass - halo mass relation changes little with redshift up to $z=3$. The bias of HI sources shows a scale dependence that gets more pronounced with increasing redshift.

We present predictions for the UV-to-mm extragalactic background light (EBL) from a recent version of the GALFORM semi-analytical model of galaxy formation which invokes a top-heavy stellar initial mass function (IMF) for galaxies undergoing dynamically-triggered bursts of star formation. We combine GALFORM with the GRASIL radiative transfer code for computing fully self-consistent UV-to-mm spectral energy distributions for each simulated galaxy, accounting for the absorption and re-emission of stellar radiation by interstellar dust. The predicted EBL is in near-perfect agreement with recent observations over the whole UV-to-mm spectrum, as is the evolution of the cosmic spectral energy distribution over the redshift range for which observations are available ($z\lesssim1$). We show that approximately 90~per~cent of the EBL is produced at $z<2$ although this shifts to higher redshifts for sub-mm wavelengths. We assess whether the top-heavy IMF in starbursts is necessary in order to reproduce the EBL at the same time as other key observables, and find that variant models with a universal solar-neighborhood IMF display poorer agreement with EBL observations over the whole UV-to-mm spectrum and fail to match the counts of galaxies in the sub-mm.

Pan-STARRS1 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 7-12 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 Pan-STARRS 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 Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI.

The Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Survey (BGS) will be a survey of bright, low redshift galaxies, which is planned to cover an area of ~14,000 sq deg in 3 passes. Each pass will cover the survey area with ~2000 pointings, each of area ~8 sq deg. The BGS is currently proposed to consist of a bright high priority sample to an r-band magnitude limit r ~ 19.5, with a fainter low priority sample to r ~ 20. The geometry of the DESI fibre positioners in the focal plane of the telescope affects the completeness of the survey, and has a non-trivial impact on clustering measurements. Using a BGS mock catalogue, we show that completeness due to fibre assignment primarily depends on the surface density of galaxies. Completeness is high (>95%) in low density regions, but very low (<10%) in the centre of massive clusters. We apply the pair inverse probability (PIP) weighting correction to clustering measurements from a BGS mock which has been through the fibre assignment algorithm. This method is only unbiased if it is possible to observe every galaxy pair. To facilitate this, we randomly promote a small fraction of the fainter sample to be high priority, and dither the set of tile positions by a small angle. We show that inverse pair weighting combined with angular upweighting provides an unbiased correction to galaxy clustering measurements for the complete 3 pass survey, and also after 1 pass, which is highly incomplete.

We apply the model relating halo concentration to formation history proposed by Ludlow et al. to merger trees generated using an algorithm based on excursion set theory. We find that while the model correctly predicts the median relation between halo concentration and mass, it underpredicts the scatter in concentration at fixed mass. Since the same model applied to N-body merger trees predicts the correct scatter, we postulate that the missing scatter is due to the lack of any environmental dependence in merger trees derived from excursion set theory. We show that a simple modification to the merger tree construction algorithm, which makes merger rates dependent on environment, can increase the scatter by the required amount, and simultaneously provide a qualitatively correct correlation between environment and formation epoch in the excursion set merger trees.

We present and analyse mock stellar catalogues that match the selection criteria and observables (including uncertainties) of the Gaia satellite data release 2 (DR2). The source are six cosmological high-resolution magneto-hydrodynamic $\Lambda$CDM zoom simulations of the formation of Milky Way analogues from the AURIGA project. Mock data are provided for stars with $V < 16$ mag, and $V < 20$ mag at $|b|>20$ degrees. The mock catalogues are made using two different methods: the public SNAPDRAGONS code, and a method based on that of Lowing et al. that preserves the phase-space distribution of the model stars. These publicly available catalogues contain 5-parameter astrometry, radial velocities, multi-band photometry, stellar parameters, dust extinction values, and uncertainties in all these quantities. In addition, we provide the gravitational potential and information on the origin of each star. By way of demonstration, we apply the mock catalogues to analyses of the young stellar disc and the stellar halo. We show that: i) the young outer stellar disc exhibits a flared distribution that is detectable in the height and vertical velocity distribution of A- and B-dwarf stars up to radii of ~15 kpc; and ii) the spin of the stellar halo out to 100 kpc can be accurately measured with Gaia DR2 RR Lyrae stars. These catalogues are well suited for comparisons with observations and should help to: i) develop and test analysis methods for the Gaia DR2 data; ii) gauge the limitations and biases of the data and iii) interpret the data in the light of theoretical predictions from realistic $ab$ $initio$ simulations of galaxy formation in the $\Lambda$CDM cosmological model.

The spherical Jeans equation (SJE) is widely used in dynamical modelling of the Milky Way (MW) halo potential. We use haloes and galaxies from the cosmological Millennium-II simulation and hydrodynamical APOSTLE simulations to investigate the performance of the SJE in recovering the underlying mass profiles of MW mass haloes. The best-fitting halo mass and concentration parameters scatter by 25% and 40% around their input values, respectively, when dark matter particles are used as tracers. This scatter becomes as large as a factor of 3 when using star particles instead. This is significantly larger than the estimated statistical uncertainty associated with the use of the SJE. The existence of correlated phase-space structures that violate the steady state assumption of the SJE as well as non-spherical geometries are the principal sources of the scatter. Binary haloes show larger scatter because they are more aspherical in shape and have a more perturbed dynamical state. Our results confirm the previous study of Wang et al. (2017) that the number of independent phase-space structures sets an intrinsic limiting precision on dynamical inferences based on the steady state assumption. Modelling with a radius-independent velocity anisotropy, or using tracers within a limited outer radius, result in significantly larger scatter, but the ensemble-averaged measurement over the whole halo sample is approximately unbiased.

The existence of two kinematically and chemically distinct stellar subpopulations in the Sculptor and Fornax dwarf galaxies offers the opportunity to constrain the density profile of their matter haloes by measuring the mass contained within the well-separated half-light radii of the two metallicity subpopulations. Walker and Penarrubia have used this approach to argue that data for these galaxies are consistent with constant-density `cores' in their inner regions and rule out `cuspy' Navarro-Frenk-White (NFW) profiles with high statistical significance, particularly in the case of Sculptor. We test the validity of these claims using dwarf galaxies in the APOSTLE (A Project Of Simulating The Local Environment) Lambda cold dark matter cosmological hydrodynamics simulations of analogues of the Local Group. These galaxies all have NFW dark matter density profiles and a subset of them develop two distinct metallicity subpopulations reminiscent of Sculptor and Fornax. We apply a method analogous to that of Walker and Penarrubia to a sample of 53 simulated dwarfs and find that this procedure often leads to a statistically significant detection of a core in the profile when in reality there is a cusp. Although multiple factors contribute to these failures, the main cause is a violation of the assumption of spherical symmetry upon which the mass estimators are based. The stellar populations of the simulated dwarfs tend to be significantly elongated and, in several cases, the two metallicity populations have different asphericity and are misaligned. As a result, a wide range of slopes of the density profile are inferred depending on the angle from which the galaxy is viewed.

We present predictions for the outcome of deep galaxy surveys with the $James$ $Webb$ $Space$ $Telescope$ ($JWST$) obtained from a physical model of galaxy formation in $\Lambda$CDM. We use the latest version of the GALFORM model, embedded within a new ($800$ Mpc)$^{3}$ dark matter only simulation with a halo mass resolution of $M_{\rm halo}>2\times10^{9}$ $h^{-1}$ M$_{\odot}$. For computing full UV-to-mm galaxy spectral energy distributions, including the absorption and emission of radiation by dust, we use the spectrophotometric radiative transfer code GRASIL. The model is calibrated to reproduce a broad range of observational data at $z\lesssim6$, and we show here that it can also predict evolution of the rest-frame far-UV luminosity function for $7\lesssim z\lesssim10$ which is in good agreement with observations. We make predictions for the evolution of the luminosity function from $z=16$ to $z=0$ in all broadband filters on the Near InfraRed Camera (NIRCam) and Mid InfraRed Instrument (MIRI) on $JWST$ and present the resulting galaxy number counts and redshift distributions. Our fiducial model predicts that $\sim1$ galaxy per field of view will be observable at $z\sim11$ for a $10^4$ s exposure with NIRCam. A variant model, which produces a higher redshift of reionization in better agreement with $Planck$ data, predicts number densities of observable galaxies $\sim5\times$ greater at this redshift. Similar observations with MIRI are predicted not to detect any galaxies at $z\gtrsim6$. We also make predictions for the effect of different exposure times on the redshift distributions of galaxies observable with $JWST$, and for the angular sizes of galaxies in $JWST$ bands.

It is now possible for hydrodynamical simulations to reproduce a representative galaxy population. Accordingly, it is timely to assess critically some of the assumptions of traditional semi-analytic galaxy formation models. We use the Eagle simulations to assess assumptions built into the Galform semi-analytic model, focussing on those relating to baryon cycling, angular momentum and feedback. We show that the assumption in Galform that newly formed stars have the same specific angular momentum as the total disc leads to a significant overestimate of the total stellar specific angular momentum of disc galaxies. In Eagle, stars form preferentially out of low specific angular momentum gas in the interstellar medium (ISM) due to the assumed gas density threshold for stars to form, leading to more realistic galaxy sizes. We find that stellar mass assembly is similar between Galform and Eagle but that the evolution of gas properties is different, with various indications that the rate of baryon cycling in Eagle is slower than is assumed in Galform. Finally, by matching individual galaxies between Eagle and Galform, we find that an artificial dependence of AGN feedback and gas infall rates on halo mass doubling events in Galform drives most of the scatter in stellar mass between individual objects. Put together our results suggest that the Galform semi-analytic model can be significantly improved in light of recent advances.

N-body dark matter simulations of structure formation in the $\Lambda$CDM model predict a population of subhalos within Galactic halos that have higher central densities than inferred for satellites of the Milky Way, a tension known as the `too big to fail' problem. Proposed solutions include baryonic effects, a smaller mass for the Milky Way halo, and warm dark matter. We test these three possibilities using a semi-analytic model of galaxy formation to generate luminosity functions for Milky Way halo-analogue satellite populations, the results of which are then coupled to the Jiang & van den Bosch model of subhalo stripping to predict the subhalo $V_\mathrm{max}$ functions for the 10 brightest satellites. We find that selecting the brightest satellites (as opposed to the most massive) and modelling the expulsion of gas by supernovae at early times increases the likelihood of generating the observed Milky Way satellite $V_\mathrm{max}$ function. The preferred halo mass is $6\times10^{11}M_{\odot}$, which has a 14 percent probability to host a $V_\mathrm{max}$ function like that of the Milky Way satellites. This probability is reduced to 8 percent for a $1.0\times10^{12}M_{\odot}$ halo and to 3 percent for a $1.4\times10^{12}M_{\odot}$ halo. We conclude that the Milky Way satellite $V_\mathrm{max}$ function is compatible with a CDM cosmology, as previously found by Sawala et al. using hydrodynamic simulations. Sterile neutrino-warm dark matter models achieve a higher degree of agreement with the observations, with a maximum 35 percent chance of generating the observed Milky Way satellite $V_\mathrm{max}$ function. However, more work is required to check that the semi-analytic stripping model is calibrated correctly in the sterile neutrino cosmology, and to check if our sterile neutrino models produce sufficient numbers of faint satellites.

Particle tagging is an efficient, but approximate, technique for using cosmological N-body simulations to model the phase-space evolution of the stellar populations predicted, for example, by a semi-analytic model of galaxy formation. We test the technique developed by Cooper et al. (which we call STINGS here) by comparing particle tags with stars in a smooth particle hydrodynamic (SPH) simulation. We focus on the spherically averaged density profile of stars accreted from satellite galaxies in a Milky Way (MW)-like system. The stellar profile in the SPH simulation can be recovered accurately by tagging dark matter (DM) particles in the same simulation according to a prescription based on the rank order of particle binding energy. Applying the same prescription to an N-body version of this simulation produces a density profile differing from that of the SPH simulation by <10 per cent on average between 1 and 200 kpc. This confirms that particle tagging can provide a faithful and robust approximation to a self-consistent hydrodynamical simulation in this regime (in contradiction to previous claims in the literature). We find only one systematic effect, likely due to the collisionless approximation, namely that massive satellites in the SPH simulation are disrupted somewhat earlier than their collisionless counterparts. In most cases this makes remarkably little difference to the spherically averaged distribution of their stellar debris. We conclude that, for galaxy formation models that do not predict strong baryonic effects on the present-day DM distribution of MW-like galaxies or their satellites, differences in stellar halo predictions associated with the treatment of star formation and feedback are much more important than those associated with the dynamical limitations of collisionless particle tagging.

Dark matter subhalos are the remnants of (incomplete) halo mergers. Identifying them and establishing their evolutionary links in the form of merger trees is one of the most important applications of cosmological simulations. The Hierachical Bound-Tracing (HBT) code identifies halos as they form and tracks their evolution as they merge, simultaneously detecting subhalos and building their merger trees. Here we present a new implementation of this approach, HBT+, that is much faster, more user friendly, and more physically complete than the original code. Applying HBT+ to cosmological simulations we show that both the subhalo mass function and the peak-mass function are well fit by similar double-Schechter functions.The ratio between the two is highest at the high mass end, reflecting the resilience of massive subhalos that experience substantial dynamical friction but limited tidal stripping. The radial distribution of the most massive subhalos is more concentrated than the universal radial distribution of lower mass subhalos. Subhalo finders that work in configuration space tend to underestimate the masses of massive subhalos, an effect that is stronger in the host centre. This may explain, at least in part, the excess of massive subhalos in galaxy cluster centres inferred from recent lensing observations. We demonstrate that the peak-mass function is a powerful diagnostic of merger tree defects, and the merger trees constructed using HBT+ do not suffer from the missing or switched links that tend to afflict merger trees constructed from more conventional halo finders. We make the HBT+ code publicly available.

We make use of a catalog of 1600 Pan-STARRS1 groups produced by the probability friends-of-friends algorithm to explore how the galaxy properties, i.e. the specific star formation rate (SSFR) and quiescent fraction, depend on stellar mass and group-centric radius. The work is the extension of Lin et al. (2014). In this work, powered by a stacking technique plus a background subtraction for contamination removal, a finer correction and more precise results are obtained than in our previous work. We find that while the quiescent fraction increases with decreasing group-centric radius the median SSFRs of star-forming galaxies in groups at fixed stellar mass drop slightly from the field toward the group center. This suggests that the major quenching process in groups is likely a fast mechanism. On the other hand, a reduction in SSFRs by ~0.2 dex is seen inside clusters as opposed to the field galaxies. If the reduction is attributed to the slow quenching effect, the slow quenching process acts dominantly in clusters. In addition, we also examine the density-color relation, where the density is defined by using a sixth-nearest neighbor approach. Comparing the quiescent fractions contributed from the density and radial effect, we find that the density effect dominates over the massive group or cluster galaxies, and the radial effect becomes more effective in less massive galaxies. The results support mergers and/or starvation as the main quenching mechanisms in the group environment, while harassment and/or starvation dominate in clusters.

Many dynamical models of the Milky Way halo require assumptions that the distribution function of a tracer population should be independent of time (i.e., a steady state distribution function) and that the underlying potential is spherical. We study the limitations of such modelling by applying a general dynamical model with minimal assumptions to a large sample of galactic haloes from cosmological $N$-body and hydrodynamical simulations. Using dark matter particles as dynamical tracers, we find that the systematic uncertainties in the measured mass and concentration parameters typically have an amplitude of 25% to 40%. When stars are used as tracers, however, the systematic uncertainties can be as large as a factor of $2-3$. The systematic uncertainties are not reduced by increasing the tracer sample size and vary stochastically from halo to halo. These systematic uncertainties are mostly driven by underestimated statistical noise caused by correlated phase-space structures that violate the steady state assumption. The number of independent phase-space structures inferred from the uncertainty level sets a limiting sample size beyond which a further increase no longer significantly improves the accuracy of dynamical inferences. The systematic uncertainty level is determined by the halo merger history, the shape and environment of the halo. Our conclusions apply generally to any spherical steady-state model.

Small distortions in the images of Einstein rings or giant arcs offer the exciting prospect of detecting dark matter haloes or subhaloes of mass below $10^9$M$_{\odot}$, most of which are too small to have made a visible galaxy. A very large number of such haloes are predicted to exist in the cold dark matter model of cosmogony; in contrast other models, such as warm dark matter, predict no haloes below a mass of this order which depends on the properties of the warm dark matter particle. Attempting to detect these small perturbers could therefore discriminate between different kinds of dark matter particles, and even rule out the cold dark matter model altogether. Globular clusters in the lens galaxy also induce distortions in the image which could, in principle, contaminate the test. Here, we investigate the population of globular clusters in six early type galaxies in the Virgo cluster. We find that the number density of globular clusters of mass $\sim10^6$M$_{\odot}$ is comparable to that of the dark matter perturbers (including subhaloes in the lens and haloes along the line-of-sight). We show that the very different degrees of mass concentration in globular clusters and dark matter haloes result in different lensing distortions. These are detectable with milli-arcsecond resolution imaging which can distinguish between globular cluster and dark matter halo signals.

Future galaxy surveys require realistic mock catalogues to understand and quantify systematics in order to make precise cosmological measurements. We present a halo lightcone catalogue and halo occupation distribution (HOD) galaxy catalogue built using the Millennium-XXL (MXXL) simulation. The halo catalogue covers the full sky, extending to z = 2 with a mass resolution of ~1e11 Msun/h . We use this to build a galaxy catalogue, which has an r-band magnitude limit of r < 20.0, with a median redshift of z~0.2. A Monte Carlo HOD method is used to assign galaxies to the halo lightcone catalogue, and we evolve the HODs to reproduce a target luminosity function; by construction, the luminosity function of galaxies in the mock is in agreement with the Sloan Digital Sky Survey (SDSS) at low redshifts and the Galaxy and Mass Assembly (GAMA) survey at high redshifts. A Monte Carlo method is used to assign a 0.1(g-r) colour to each galaxy, and the colour distribution of galaxies at different redshifts agrees with measurements from GAMA. The clustering of galaxies in the mock for galaxies in different magnitude and redshift bins is in good agreement with measurements from SDSS and GAMA, and the colour-dependent clustering is in reasonable agreement. We show that the baryon acoustic oscillation (BAO) can be measured in the mock catalogue, and the redshift space distortions (RSDs) are in agreement with measurements from SDSS, illustrating that this catalogue will be useful for upcoming surveys.

The Gaia astrometric mission may offer an unprecedented opportunity to discover new tidal streams in the Galactic halo. To test this, we apply nGC3, a great-circle-cell count method that combines position and proper motion data to identify streams, to ten mock Gaia catalogues of K giants and RR Lyrae stars constructed from cosmological simulations of Milky Way analogues. We analyse two sets of simulations, one using a combination of $N$-body and semi-analytical methods which has extremely high resolution, the other using hydro-dynamical methods, which captures the dynamics of baryons, including the formation of an in situ halo. These ten realisations of plausible Galactic merger histories allow us to assess the potential for the recovery of tidal streams in different Milky Way formation scenarios. We include the Gaia~selection function and observational errors in these mock catalogues. We find that the nGC3 method has a well-defined detection boundary in the space of stream width and projected overdensity, that can be predicted based on direct observables alone. We predict that about 4-13 dwarf galaxy streams can be detected in a typical Milky Way-mass halo with Gaia+nGC3, with an estimated efficiency of $>$80\% inside the detection boundary. The progenitors of these streams are in the mass range of the classical dwarf galaxies and may have been accreted as early as redshift $\sim3$. Finally, we analyse how different possible extensions of the Gaia mission will improve the detection of tidal streams.

The standard Halo Occupation Distribution (HOD) models were originally developed based on results from semi-analytic and hydrodynamical galaxy formation models. Those models have since progressed, in particular to include AGN feedback to match the galaxy luminosity function in a universe with the observed baryon fraction. AGN feedback affects the relationship between galaxy stellar mass and luminosity, in particular making the relationship non-monotonic. For matched number density samples, galaxies in luminosity-threshold samples occupy a different range of halo masses from those in stellar-mass-threshold samples. We find that the shapes of the HODs of luminosity-threshold samples are slightly more complicated in semi-analytic galaxy formation models that include AGN feedback than are assumed by standard HOD models. We also find that subhalo abundance matching (SHAM) does not preserve these non-standard shapes. We show that catalogues created using SHAM and the semi-analytic model Galform that have the same large-scale 2-point clustering by construction have different void probability functions (VPFs) in both real and redshift space. We find that these differences arise from the different HOD shapes, as opposed to assembly bias, which indicates that the VPF could be used to test the suitability of an HOD model with real data.

Placing bright sub-millimetre galaxies (SMGs) within the broader context of galaxy formation and evolution requires accurate measurements of their clustering, which can constrain the masses of their host dark matter halos. Recent work has shown that the clustering measurements of these galaxies may be affected by a `blending bias,' which results in the angular correlation function of the sources extracted from single-dish imaging surveys being boosted relative to that of the underlying galaxies. This is due to confusion introduced by the coarse angular resolution of the single-dish telescope and could lead to the inferred halo masses being significantly overestimated. We investigate the extent to which this bias affects the measurement of the correlation function of SMGs when it is derived via a cross-correlation with a more abundant galaxy population. We find that the blending bias is essentially the same as in the auto-correlation case and conclude that the best way to reduce its effects is to calculate the angular correlation function using SMGs in narrow redshift bins. Blending bias causes the inferred host halo masses of the SMGs to be overestimated by a factor of $\sim6$ when a redshift interval of $\delta z=3$ is used. However, this reduces to a factor of $\sim2$ for $\delta z=0.5$. The broadening of photometric redshift probability distributions with increasing redshift can therefore impart a mild halo `downsizing' effect onto the inferred host halo masses, though this trend is not as strong as seen in recent observational studies.

We compare observed far infra-red/sub-millimetre (FIR/sub-mm) galaxy spectral energy distributions (SEDs) of massive galaxies ($M_{\star}\gtrsim10^{10}$ $h^{-1}$M$_{\odot}$) derived through a stacking analysis with predictions from a new model of galaxy formation. The FIR SEDs of the model galaxies are calculated using a self-consistent model for the absorption and re-emission of radiation by interstellar dust based on radiative transfer calculations and global energy balance arguments. Galaxies are selected based on their position on the specific star formation rate (sSFR) - stellar mass ($M_{\star}$) plane. We identify a main sequence of star-forming galaxies in the model, i.e. a well defined relationship between sSFR and $M_\star$, up to redshift $z\sim6$. The scatter of this relationship evolves such that it is generally larger at higher stellar masses and higher redshifts. There is remarkable agreement between the predicted and observed average SEDs across a broad range of redshifts ($0.5\lesssim z\lesssim4$) for galaxies on the main sequence. However, the agreement is less good for starburst galaxies at $z\gtrsim2$, selected here to have elevated sSFRs$>10\times$ the main sequence value. We find that the predicted average SEDs are robust to changing the parameters of our dust model within physically plausible values. We also show that the dust temperature evolution of main sequence galaxies in the model is driven by star formation on the main sequence being more burst-dominated at higher redshifts.

The defining characteristic of the cold dark matter (CDM) hypothesis is the presence of a very large number of low-mass haloes, too small to have made a visible galaxy. Other hypotheses for the nature of the dark matter, such as warm dark matter (WDM), predict a much smaller number of such low-mass haloes. Strong lensing systems offer the possibility of detecting small-mass haloes through the distortions they induce in the lensed image. Here we show that the main contribution to the image distortions comes from haloes along the line of sight rather than subhaloes in the lens as has normally been assumed so far. These interlopers enhance the differences between the predictions of CDM and WDM models. We derive the total perturber mass function, including both subhaloes and interlopers, and show that measurements of approximately 20 strong lens systems with a detection limit of $M_{\rm low}=10^7 h^{-1} M_{\odot}$ would distinguish (at 3 sigma) between CDM and a WDM model consisting of 7 keV sterile neutrinos such as those required to explain the recently detected 3.5 keV X-ray emission line from the centres of galaxies and clusters.

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 sterile neutrino is a viable dark matter candidate that can be produced in the early Universe via non-equilibrium processes, and would therefore possess a highly non-thermal spectrum of primordial velocities. In this paper we analyse the process of structure formation with this class of dark matter particles. To this end we construct primordial dark matter power spectra as a function of the lepton asymmetry, $L_6$, that is present in the primordial plasma and leads to resonant sterile neutrino production. We compare these power spectra with those of thermally produced dark matter particles and show that resonantly produced sterile neutrinos are much colder than their thermal relic counterparts. We also demonstrate that the shape of these power spectra is not determined by the free-streaming scale alone. We then use the power spectra as an input for semi-analytic models of galaxy formation in order to predict the number of luminous satellite galaxies in a Milky Way-like halo. By assuming that the mass of the Milky Way halo must be no more than $2\times10^{12}M_{\odot}$ (the adopted upper bound based on current astronomical observations) we are able to constrain the value of $L_6$ for $M_s\le 8$~keV. We also show that the range of $L_6$ that is in best agreement with the 3.5~keV line (if produced by decays of 7~keV sterile neutrino) requires that the Milky Way halo has a mass no smaller than $1.5\times10^{12}M_{\odot}$. Finally, we compare the power spectra obtained by direct integration of the Boltzmann equations for a non-resonantly produced sterile neutrino with the fitting formula of Viel~et~al. and find that the latter significantly underestimates the power amplitude on scales relevant to satellite galaxies.

We present a model for the dynamical evolution of subhaloes based on an approach combining numerical and analytical methods. Our method is based on tracking subhaloes in an N-body simulation up to the last point that it can be resolved, and applying an analytic prescription for its merger timescale that takes dynamical friction and tidal disruption into account. When applied to cosmological N-body simulations with mass resolutions that differ by two orders of magnitude, the technique produces halo occupation distributions that agree to within 3%.

We derive the expression for the observed redshift in the weak field limit in the observer's past light cone, including all relativistic terms up to second order in velocity. We then apply it to compute the cluster-galaxy cross-correlation functions (CGCF) using N-body simulations. The CGCF is asymmetric along the line of sight (LOS) owing to the presence of the small second order terms such as the gravitational redshift (GRedshift). We identify two systematics in the modelling of the GRedshift signal in stacked clusters. First, it is affected by the morphology of dark matter haloes and the large-scale cosmic-web. The non-spherical distribution of galaxies around the central halo and the presence of neighbouring clusters systematically reduce the GRedshift signal. This bias is approximately 20% for $M_{\rm min}\simeq 10^{14} {\rm M_{\odot}}/h$, and is more than 50% for haloes with $M_{\rm min}\simeq 2\times 10^{13} {\rm M_{\odot}}/h$ at $r>$4 Mpc/$h$. Second, the best-fit gravitational redshift profiles as well as the profiles of all other relativistic terms are found to be significantly different in velocity space compared to their real space versions. We discuss some subtleties relating to these effects in velocity space. We also find that the S/N of the GRedshift signal increases with decreasing halo mass.

Accurate modeling of galaxy formation in a hierarchical, cold dark matter universe requires the use of sufficiently high-resolution merger trees to obtain convergence in the predicted properties of galaxies. When semi-analytic galaxy formation models are applied to cosmological N-body simulation merger trees, it is often the case that those trees have insufficient resolution to give converged galaxy properties. We demonstrate a method to augment the resolution of N-body merger trees by grafting in branches of Monte Carlo merger trees with higher resolution, but which are consistent with the pre-existing branches in the N-body tree. We show that this approach leads to converged galaxy properties.

We present a new version of the GALFORM semi-analytical model of galaxy formation. This brings together several previous developments of GALFORM into a single unified model, including a different initial mass function (IMF) in quiescent star formation and in starbursts, feedback from active galactic nuclei supressing gas cooling in massive halos, and a new empirical star formation law in galaxy disks based on their molecular gas content. In addition, we have updated the cosmology, introduced a more accurate treatment of dynamical friction acting on satellite galaxies, and updated the stellar population model. The new model is able to simultaneously explain both the observed evolution of the K-band luminosity function and stellar mass function, and the number counts and redshift distribution of sub-mm galaxies selected at 850 mu. This was not previously achieved by a single physical model within the LambdaCDM framework, but requires having an IMF in starbursts that is somewhat top-heavy. The new model is tested against a wide variety of observational data covering wavelengths from the far-UV to sub-mm, and redshifts from z=0 to z=6, and is found to be generally successful. These observations include the optical and near-IR luminosity functions, HI mass function, fraction of early type galaxies, Tully-Fisher, metallicity-luminosity and size-luminosity relations at z=0, as well as far-IR number counts, and far-UV luminosity functions at z ~ 3-6. [abridged]

We compare global predictions from the EAGLE hydrodynamical simulation, and two semi-analytic (SA) models of galaxy formation, L-GALAXIES and GALFORM. All three models include the key physical processes for the formation and evolution of galaxies and their parameters are calibrated against a small number of observables at $z\approx 0$. The two SA models have been applied to merger trees constructed from the EAGLE dark matter only simulation. We find that at $z\leq 2$, both the galaxy stellar mass functions for stellar masses $M_*<10^{10.5} M_{\odot}$ and the median specific star formation rates (sSFRs) in the three models agree to better than $0.4$~dex. The evolution of the sSFR predicted by the three models closely follows the mass assembly history of dark matter haloes. In both EAGLE and L-GALAXIES there are more central passive galaxies with $M_*<10^{9.5} M_{\odot}$ than in L-GALAXIES. This difference is related to galaxies that have entered and then left a larger halo and which are treated as satellites in GALFORM. In the range $0<z<1$, the slope of the evolution of the star formation rate density in EAGLE is a factor of $\approx 1.5$ steeper than for the two SA models. The median sizes for galaxies with $M_*>10^{9.5} M_{\odot}$ differ in some instances by an order of magnitude, while the stellar mass-size relation in EAGLE is a factor of $\approx 2$ tighter than for the two SA models. Our results suggest the need for a revision of how SA models treat the effect of baryonic self-gravity on the underlying dark matter. The treatment of gas flows in the models needs to be revised based on detailed comparison with observations to understand in particular the evolution of the stellar mass-metallicity relation.

We use a suite of cosmological simulations to study the mass-concentration-redshift relation, $c({\rm M},z)$, of dark matter halos. Our simulations include standard $\Lambda$-cold dark matter (CDM) models, and additional runs with truncated power spectra, consistent with a thermal warm dark matter (WDM) scenario. We find that the mass profiles of CDM and WDM halos are self-similar and well approximated by the Einasto profile. The $c({\rm M},z)$ relation of CDM halos is monotonic: concentrations decrease with increasing virial mass at fixed redshift, and decrease with increasing redshift at fixed mass. The mass accretion histories (MAHs) of CDM halos are also scale-free, and can be used to infer concentrations directly. These results do not apply to WDM halos: their MAHs are not scale-free because of the characteristic scale imposed by the power-spectrum suppression. Further, the WDM $c({\rm M},z)$ relation is non-monotonic: concentrations peak at a mass scale dictated by the truncation scale, and decrease at higher and lower masses. We show that the assembly history of a halo can still be used to infer its concentration, provided that the total mass of its progenitors is considered (the "collapsed mass history"; CMH), rather than just that of its main ancestor. This exploits the scale-free nature of CMHs to derive a simple scaling that reproduces the mass-concentration-redshift relation of both CDM and WDM halos over a vast range of halo masses and redshifts. Our model therefore provides a robust account of the mass, redshift, cosmology and power spectrum dependence of dark matter halo concentrations.

We present predictions for the clustering of galaxies selected by their emission at far infra-red (FIR) and sub-millimetre wavelengths. This includes the first predictions for the effect of clustering biases induced by the coarse angular resolution of single-dish telescopes at these wavelengths. We combine a new version of the GALFORM model of galaxy formation with a self-consistent model for calculating the absorption and re-emission of radiation by interstellar dust. Model galaxies selected at $850$ $\mu$m reside in dark matter halos of mass $M_{\rm halo}\sim10^{11.5}-10^{12}$ $h^{-1}$ M$_{\odot}$, independent of redshift (for $0.2\lesssim z\lesssim4$) or flux (for $0.25\lesssim S_{850\mu\rm m}\lesssim4$ mJy). At $z\sim2.5$, the brightest galaxies ($S_{850\mu\rm m}>4$ mJy) exhibit a correlation length of $r_{0}=5.5_{-0.5}^{+0.3}$ $h^{-1}$ Mpc, consistent with observations. We show that these galaxies have descendants with stellar masses $M_{\star}\sim10^{11}$ $h^{-1}$ M$_{\odot}$ occupying halos spanning a broad range in mass $M_{\rm halo}\sim10^{12}-10^{14}$ $h^{-1}$ M$_{\odot}$. The FIR emissivity at shorter wavelengths ($250$, $350$ and $500$ $\mu$m) is also dominated by galaxies in the halo mass range $M_{\rm halo}\sim10^{11.5}-10^{12}$ $h^{-1}$ M$_{\odot}$, again independent of redshift (for $0.5\lesssim z\lesssim5$). We compare our predictions for the angular power spectrum of cosmic infra-red background anisotropies at these wavelengths with observations, finding agreement to within a factor of $\sim2$ over all scales and wavelengths, an improvement over earlier versions of the model. Simulating images at $850$ $\mu$m, we show that confusion effects boost the measured angular correlation function on all scales by a factor of $\sim4$. This has important consequences, potentially leading to inferred halo masses being overestimated by an order of magnitude.

The cold dark matter (CDM) cosmological model unambigously predicts that a large number of haloes should survive as subhaloes when they are accreted into a larger halo. The CDM model would be ruled out if such substructures were shown not to exist. By contrast, if the dark matter consists of warm particles (WDM), then below a threshold mass that depends on the particle mass far fewer substructures would be present. Finding subhaloes below a certain mass would then rule out warm particle masses below some value. Strong gravitational lensing provides a clean method to measure the subhalo mass function through distortions in the structure of Einstein rings and giant arcs.Using mock lensing observations constructed from high-resolution N-body simulations, we show that measurements of approximately 100 strong lens systems with a detection limit of $10^7 h^{-1} M_{\odot}$ would clearly distinguish CDM from WDM in the case where this consists of 7 keV sterile neutrinos such as those that might be responsible for the 3.5 keV X-ray emission line recently detected in galaxies and clusters.

N-body simulations suggest that the substructures that survive inside dark matter haloes follow universal distributions in mass and radial number density. We demonstrate that a simple analytical model can explain these subhalo distributions as resulting from tidal stripping which increasingly reduces the mass of subhaloes with decreasing halo-centric distance. As a starting point, the spatial distribution of subhaloes of any given infall mass is shown to be largely indistinguishable from the overall mass distribution of the host halo. Using a physically motivated statistical description of the amount of mass stripped from individual subhaloes, the model fully describes the joint distribution of subhaloes in final mass, infall mass and radius. As a result, it can be used to predict several derived distributions involving combinations of these quantities including, but not limited to, the universal subhalo mass function, the subhalo spatial distribution, the gravitational lensing profile, the dark matter annihilation radiation profile and boost factor. This model clarifies a common confusion when comparing the spatial distributions of galaxies and subhaloes, the so called "anti-bias", as a simple selection effect. We provide a Python code SubGen for populating haloes with subhaloes at http://icc.dur.ac.uk/data/

We develop two general methods to infer the gravitational potential of a system using steady-state tracers, i.e., tracers with a time-independent phase-space distribution. Combined with the phase-space continuity equation, the time independence implies a universal Orbital Probability Density Function (oPDF) $\mathrm{d} P(\lambda|{\rm orbit})\propto \mathrm{d} t$, where $\lambda$ is the coordinate of the particle along the orbit. The oPDF is equivalent to Jeans theorem, and is the key physical ingredient behind most dynamical modelling of steady-state tracers. In the case of a spherical potential, we develop a likelihood estimator that fits analytical potentials to the system, and a non-parametric method ("phase-mark") that reconstructs the potential profile, both assuming only the oPDF. The methods involve no extra assumptions about the tracer distribution function and can be applied to tracers with any arbitrary distribution of orbits, with possible extension to non-spherical potentials. The methods are tested on Monte Carlo samples of steady-state tracers in dark matter haloes to show that they are unbiased as well as efficient. A fully documented \textsc{C/Python} code implementing our method is freely available at a GitHub repository linked from \url{http://icc.dur.ac.uk/data/#oPDF}.

Using realistic cosmological simulations of Milky Way sized haloes, we study their dynamical state and the accuracy of inferring their mass profiles with steady-state models of dynamical tracers. We use a new method that describes the phase-space distribution of a steady-state tracer population in a spherical potential without any assumption regarding the distribution of their orbits. Applying the method to five haloes from the Aquarius $\Lambda$CDM N-body simulation, we find that dark matter particles are an accurate tracer that enables the halo mass and concentration parameters to be recovered with an accuracy of $5\%$. Assuming a potential profile of the NFW form does not significantly affect the fits in most cases, except for halo A whose density profile differs significantly from the NFW form, leading to a $30\%$ bias in the dynamically fitted parameters. The existence of substructures in the dark matter tracers only affects the fits by $\sim 1\%$. Applying the method to mock stellar haloes generated by a particle-tagging technique, we find the stars are farther from equilibrium than dark matter particles, yielding a systematic bias of $\sim 20\%$ in the inferred mass and concentration parameter. The level of systematic biases obtained from a conventional distribution function fit to stars is comparable to ours, while similar fits to DM tracers are significantly biased in contrast to our fits. In line with previous studies, the mass bias is much reduced near the tracer half-mass radius.

The Gaussianization transform has been proposed as a method to remove the issues of scale-dependent galaxy bias and nonlinearity from galaxy clustering statistics, but these benefits have yet to be thoroughly tested for realistic galaxy samples. In this paper, we test the effectiveness of the Gaussianization transform for different galaxy types by applying it to realistic simulated blue and red galaxy samples. We show that in real space, the shapes of the Gaussianized power spectra of both red and blue galaxies agree with that of the underlying dark matter, with the initial power spectrum, and with each other to smaller scales than do the statistics of the usual (untransformed) density field. However, we find that the agreement in the Gaussianized statistics breaks down in redshift space. We attribute this to the fact that red and blue galaxies exhibit very different fingers of god in redshift space. After applying a finger-of-god compression, the agreement on small scales between the Gaussianized power spectra is restored. We also compare the Gaussianization transform to the clipped galaxy density field and find that while both methods are effective in real space, they have more complicated behaviour in redshift space. Overall, we find that Gaussianization can be useful in recovering the shape of the underlying dark matter power spectrum to k ~ 0.5 h/Mpc and of the initial power spectrum to k ~ 0.4 h/Mpc in certain cases at z = 0.

We present an analysis of the predictions made by the Galform semi-analytic galaxy formation model for the evolution of the relationship between stellar mass and halo mass. We show that for the standard implementations of supernova feedback and gas reincorporation used in semi-analytic models, this relationship is predicted to evolve weakly over the redshift range 0<z<4. Modest evolution in the median stellar mass versus halo mass (SHM) relationship implicitly requires that, at fixed halo mass, the efficiency of stellar mass assembly must be almost constant with cosmic time. We show that in our model, this behaviour can be understood in simple terms as a result of a constant efficiency of gas reincorporation, and an efficiency of SNe feedback that is, on average, constant at fixed halo mass. We present a simple explanation of how feedback from active galactic nuclei (AGN) acts in our model to introduce a break in the SHM relation whose location is predicted to evolve only modestly. Finally, we show that if modifications are introduced into the model such that, for example, the gas reincorporation efficiency is no longer constant, the median SHM relation is predicted to evolve significantly over 0<z<4. Specifically, we consider modifications that allow the model to better reproduce either the evolution of the stellar mass function or the evolution of average star formation rates inferred from observations.

The mass of the dark matter halo of the Milky Way can be estimated by fitting analytical models to the phase-space distribution of dynamical tracers. We test this approach using realistic mock stellar halos constructed from the Aquarius N-body simulations of dark matter halos in the $\Lambda$CDM cosmology. We extend the standard treatment to include a Navarro-Frenk-White (NFW) potential and use a maximum likelihood method to recover the parameters describing the simulated halos from the positions and velocities of their mock halo stars. We find that the estimate of halo mass is highly correlated with the estimate of halo concentration. The best-fit halo masses within the virial radius, $R_{200}$, are biased, ranging from a 40\% underestimate to a 5\% overestimate in the best case (when the tangential velocities of the tracers are included). There are several sources of bias. Deviations from dynamical equilibrium can potentially cause significant bias; deviations from spherical symmetry are relatively less important. Fits to stars at different galactocentric radii can give different mass estimates. By contrast, the model gives good constraints on the mass within the half-mass radius of tracers even when restricted to tracers within 60kpc. The recovered velocity anisotropies of tracers, $\beta$, are biased systematically, but this does not affect other parameters if tangential velocity data are used as constraints.

We study the formation of stellar haloes in three Milky Way-mass galaxies using cosmological smoothed particle hydrodynamics simulations, focusing on the subset of halo stars that form in situ, as opposed to those accreted from satellites. In situ stars in our simulations dominate the stellar halo out to 20 kpc and account for 30-40 per cent of its total mass. We separate in situ halo stars into three straightforward, physically distinct categories according to their origin: stars scattered from the disc of the main galaxy ('heated disc'), stars formed from gas smoothly accreted on to the halo ('smooth' gas) and stars formed in streams of gas stripped from infalling satellites ('stripped' gas). We find that most belong to the stripped gas category. Those originating in smooth gas outside the disc tend to form at the same time and place as the stripped-gas population, suggesting that their formation is associated with the same gas-rich accretion events. The scattered disc star contribution is negligible overall but significant in the Solar neighbourhood, where ~90 per cent of stars on eccentric orbits once belonged to the disc. However, the distinction between halo and thick disc in this region is highly ambiguous. The chemical and kinematic properties of the different components are very similar at the present day, but the global properties of the in situ halo differ substantially between the three galaxies in our study. In our simulations, the hierarchical buildup of structure is the driving force behind not only the accreted stellar halo, but also those halo stars formed in situ.

We measure the projected 2-point correlation function of galaxies in the 180 deg$^2$ equatorial regions of the GAMA II survey, for four different redshift slices between z = 0.0 and z=0.5. To do this we further develop the Cole (2011) method of producing suitable random catalogues for the calculation of correlation functions. We find that more r-band luminous, more massive and redder galaxies are more clustered. We also find that red galaxies have stronger clustering on scales less than ~3 $h^{-1}$ Mpc. We compare to two different versions of the GALFORM galaxy formation model, Lacey et al (in prep.) and Gonzalez-Perez et al. (2014), and find that the models reproduce the trend of stronger clustering for more massive galaxies. However, the models under predict the clustering of blue galaxies, can incorrectly predict the correlation function on small scales and under predict the clustering in our sample of galaxies with ~3$L_r$ . We suggest possible avenues to explore to improve these cluster- ing predictions. The measurements presented in this paper can be used to test other galaxy formation models, and we make the measurements available online to facilitate this.

We report results of a study of Planck Sunyaev-Zel'dovich effect (SZE) selected galaxy cluster candidates using the Panoramic Survey Telescope & Rapid Response System (Pan-STARRS) 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 Pan-STARRS1 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.

We use the WISE-2MASS infrared galaxy catalog matched with Pan-STARRS1 (PS1) galaxies to search for a supervoid in the direction of the Cosmic Microwave Background Cold Spot. Our imaging catalog has median redshift $z\simeq 0.14$, and we obtain photometric redshifts from PS1 optical colours to create a tomographic map of the galaxy distribution. The radial profile centred on the Cold Spot shows a large low density region, extending over 10's of degrees. Motivated by previous Cosmic Microwave Background results, we test for underdensities within two angular radii, $5^\circ$, and $15^\circ$. The counts in photometric redshift bins show significantly low densities at high detection significance, $\gtrsim 5 \sigma$ and $\gtrsim 6 \sigma$, respectively, for the two fiducial radii. The line-of-sight position of the deepest region of the void is $z\simeq 0.15-0.25$. Our data, combined with an earlier measurement by Granett et al. 2010, are consistent with a large $R_{\rm void}=(220 \pm 50) h^{-1}Mpc $ supervoid with $\delta_{m} \simeq -0.14 \pm 0.04$ centered at $z=0.22\pm0.03$. Such a supervoid, constituting at least a $\simeq 3.3\sigma$ fluctuation in a Gaussian distribution of the $\Lambda CDM$ model, is a plausible cause for the Cold Spot.

We present distributions of orbital parameters of infalling satellites of $\Lambda$CDM haloes in the mass range $10^{12}-10^{14}$M$_\odot$, which represent the initial conditions for the subsequent evolution of substructures within the host halo. We use merger trees constructed in a high resolution cosmological N-body simulation to trace satellite haloes, and identify the time of infall. We find signficant trends in the distribution of orbital parameters with both the host halo mass and the ratio of satellite-to-host halo masses. For all host halo masses, satellites whose infall mass is a larger fraction of the host halo mass have more eccentric, radially biased orbits. At fixed satellite-to-host halo mass ratio, high mass haloes are biased towards accreting satellites on slightly more radial orbits. To charactise the orbital distributions fully requires fitting the correlated bivariate distribution of two chosen orbital parameters (e.g. radial and tangential velocity or energy and angular momentum). We provide simple fits to one choice of the bivariate distributions, which when transformed faithfully, captures the behaviour of any of the projected one-dimensional distributions.

We present a new technique for creating mock catalogues of the individual stars that make up the accreted component of stellar haloes in cosmological simulations and show how the catalogues can be used to test and interpret observational data. The catalogues are constructed from a combination of methods. A semi-analytic galaxy formation model is used to calculate the star formation history in haloes in an N-body simulation and dark matter particles are tagged with this stellar mass. The tags are converted into individual stars using a stellar population synthesis model to obtain the number density and evolutionary stage of the stars, together with a phase-space sampling method that distributes the stars while ensuring that the phase-space structure of the original N-body simulation is maintained. A set of catalogues based on the $\Lambda$CDM Aquarius simulations of Milky Way mass haloes have been created and made publicly available on a website. Two example applications are discussed that demonstrate the power and flexibility of the mock catalogues. We show how the rich stellar substructure that survives in the stellar halo precludes a simple measurement of its density profile and demonstrate explicitly how pencil-beam surveys can return almost any value for the slope of the profile. We also show that localized variations in the abundance of particular types of stars, a signature of differences in the composition of stellar populations, allow streams to be easily identified.

Recent observational evidence suggests that the coarse angular resolution ($\sim20"$ FWHM) of single-dish telescopes at sub-mm wavelengths has biased the observed galaxy number counts by blending together the sub-mm emission from multiple sub-mm galaxies (SMGs). We use lightcones computed from an updated implementation of the \galform semi-analytic model to generate $50$ mock sub-mm surveys of $0.5$ deg$^2$ at $850$ $\mu$m, taking into account the effects of the finite single-dish beam in a more accurate way than has been done previously. We find that blending of SMGs does lead to an enhancement of source extracted number counts at bright fluxes ($S_{\mathrm{850}\mu\mathrm{m}}\gtrsim1$ mJy). Typically, $\sim3{-}6$ galaxies contribute $90\%$ of the flux of an $S_{850\mu\mathrm{m}}=5$ mJy source and these blended galaxies are physically unassociated. We find that field-to-field variations are comparable to Poisson fluctuations for our $S_{850\mu\mathrm{m}}>5$ mJy SMG population, which has a median redshift $z_{50}=2.0$, but are greater than Poisson for the $S_{850\mu\mathrm{m}}>1$ mJy population ($z_{50}=2.8$). In a detailed comparison to a recent interferometric survey targeted at single-dish detected sources, we reproduce the difference between single-dish and interferometer number counts and find a median redshift ($z_{50}=2.5$) in excellent agreement with the observed value ($z_{50}=2.5\pm 0.2$). We also present predictions for single-dish survey number counts at $450$ and $1100$ $\mu$m, which show good agreement with observational data.

It has been argued that the specific star formation rates of star forming galaxies inferred from observational data decline more rapidly below z = 2 than is predicted by hierarchical galaxy formation models. We present a detailed analysis of this problem by comparing predictions from the GALFORM semi-analytic model with an extensive compilation of data on the average star formation rates of star-forming galaxies. We also use this data to infer the form of the stellar mass assembly histories of star forming galaxies. Our analysis reveals that the currently available data favour a scenario where the stellar mass assembly histories of star forming galaxies rise at early times and then fall towards the present day. In contrast, our model predicts stellar mass assembly histories that are almost flat below z = 2 for star forming galaxies, such that the predicted star formation rates can be offset with respect to the observational data by factors of up to 2-3. This disagreement can be explained by the level of coevolution between stellar and halo mass assembly that exists in contemporary galaxy formation models. In turn, this arises because the standard implementations of star formation and supernova feedback used in the models result in the efficiencies of these process remaining approximately constant over the lifetime of a given star forming galaxy. We demonstrate how a modification to the timescale for gas ejected by feedback to be reincorporated into galaxy haloes can help to reconcile the model predictions with the data.

We use a WISE-2MASS-Pan-STARRS1 galaxy catalog to search for a supervoid in the direction of the Cosmic Microwave Background Cold Spot. We obtain photometric redshifts using our multicolor data set to create a tomographic map of the galaxy distribution. The radial density profile centred on the Cold Spot shows a large low density region, extending over 10's of degrees. Motivated by previous Cosmic Microwave Background results, we test for underdensities within two angular radii, $5^\circ$, and $15^\circ$. Our data, combined with an earlier measurement by Granett et al 2010, are consistent with a large $R_{\rm void}=(192 \pm 15)h^{-1} Mpc $ $(2\sigma)$ supervoid with $\delta \simeq -0.13 \pm 0.03$ centered at $z=0.22\pm0.01$. Such a supervoid, constituting a $\sim3.5 \sigma$ fluctuation in the $\Lambda CDM$ model, is a plausible cause for the Cold Spot.

We use the Herschel ATLAS (H-ATLAS) Phase I data to study the conditional luminosity function of far-IR (250 um) selected galaxies in optically-selected galaxy groups from the GAMA spectroscopic survey, as well as environmental effects on the far-IR-to-optical colour. We applied two methods, which gave consistent results for the far-IR conditional luminosity functions. The direct matching method matches H-ATLAS sources to GAMA/SDSS galaxies, then links the optical counterparts to GAMA groups. The stacking method counts the number of far-IR sources within the projected radii of GAMA groups, subtracting the local background. We investigated the dependence of the far-IR (250 um) luminosity function on group mass in the range 10d12 < M_h < 10d14Msun/h and on redshift in the range 0 < z < 0.4, using a sample of 3000 groups containing H-ATLAS sources with GAMA redshifts over an area of 126 sq. deg. We find that the characteristic 250 um luminosity, L^*(250), increases with group mass up to M_h ~ 10d13Msun/h, but is roughly constant above this, while it increases with redshift at high group masses, but less so at low masses. We also find that the group far-IR luminosity-to-mass ratio L(250)/M_h increases with redshift and is higher in low-mass groups. We estimate that around 70% of the 250 um luminosity density in the local universe is contributed by groups with M_h > 10d12Msun/h. We also find that the far-IR-to-optical colours of H-ATLAS galaxies are independent of group mass over the range 10d12 < M_h <10d14Msun/h in the local universe. We also compare our observational results with recent semi-analytical models, and find that none of these galaxy formation model can reproduce the conditional far-IR luminosity functions of galaxy groups.

Standard inflationary hot big bang cosmology predicts small fluctuations in the Cosmic Microwave Background (CMB) with isotropic Gaussian statistics. All measurements support the standard theory, except for a few anomalies discovered in the Wilkinson Microwave Anisotropy Probe maps and confirmed recently by the Planck satellite. The Cold Spot is one of the most significant of such anomalies, and the leading explanation of it posits a large void that imprints this extremely cold area via the linear Integrated Sachs-Wolfe (ISW) effect due to the decay of gravitational potentials over cosmic time, or via the Rees-Sciama (RS) effect due to late-time non-linear evolution. Despite several observational campaigns targeting the Cold Spot region, to date no suitably large void was found at higher redshifts $z > 0.3$. Here we report the detection of an $R =(192 \pm 15) h^{-1}Mpc$ size supervoid of depth $\delta = -0.13 \pm 0.03$, and centred at redshift $z = 0.22$. This supervoid, possibly the largest ever found, is large enough to significantly affect the CMB via the non-linear RS effect, as shown in our Lemaitre-Tolman-Bondi framework. This discovery presents the first plausible explanation for any of the physical CMB anomalies, and raises the possibility that local large-scale structure could be responsible for other anomalies as well.

The velocity field of dark matter and galaxies reflects the continued action of gravity throughout cosmic history. We show that the low-order moments of the pairwise velocity distribution, $v_{12}$, are a powerful diagnostic of the laws of gravity on cosmological scales. In particular, the projected line-of-sight galaxy pairwise velocity dispersion, $\sigma_{12}(r)$, is very sensitive to the presence of modified gravity. Using a set of high-resolution N-body simulations we compute the pairwise velocity distribution and its projected line-of-sight dispersion for a class of modified gravity theories: the chameleon \fR gravity and Galileon gravity (cubic and quartic). The velocities of dark matter halos with a wide range of masses would exhibit deviations from General Relativity at the $(5-10)\sigma$ level. We examine strategies for detecting these deviations in galaxy redshift and peculiar velocity surveys. If detected, this signature would be a "smoking gun" for modified gravity.

We present a new algorithm which groups the subhaloes found in cosmological N- body simulations by structure finders such as SUBFIND into dark matter haloes whose formation histories are strictly hierarchical. One advantage of these `Dhaloes' over the commonly used friends-of-friends (FoF) haloes is that they retain their individual identity in cases when FoF haloes are artificially merged by tenuous bridges of particles or by an overlap of their outer diffuse haloes. Dhaloes are thus well-suited for modelling galaxy formation and their merger trees form the basis of the Durham semi-analytic galaxy formation model, GALFORM. Applying the Dhalo construction to the LambdaCDM Millennium-2 simulation we find that approximately 90% of Dhaloes have a one-to-one, bijective match with a corresponding FoF halo. The remaining 10% are typically secondary components of large FoF haloes. Although the mass functions of both types of haloes are similar, the mass of Dhaloes correlates much more tightly with the virial mass, M200, than FoF haloes. Approximately 80% of FoF and bijective and non-bijective Dhaloes are relaxed according to standard criteria. For these relaxed haloes all three types have similar concentration-M200 relations and, at fixed mass, the concentration distributions are described accurately by log-normal distributions.

We study the late-time Integrated Sachs-Wolfe (ISW) effect in $f(R)$ gravity using N-body 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 non-linear regime, relatively rapid structure formation in $f(R)$ gravity boosts the non-linear 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 top-hat 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 Sachs-Wolfe (ISW) effect. However, we also find that small voids in the halo density field produce hot spots, because they reside in contracting, larger-scale 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 higher-amplitude 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.

Using a large sample of field and group galaxies drawn from the Pan-STARRS1 Medium-Deep Survey, we present the specific star formation rate (SSFR) - stellar mass (M*) relation, as well as the quiescent fraction versus M* relation in different environments. We confirm that the fraction of quiescent galaxies is strongly dependent on environment at a fixed stellar mass, but that the amplitude and the slope of the star-forming sequence is similar between the field and groups: in other words, the SSFR-density relation at a fixed stellar mass is primarily driven by the change in the star-forming and quiescent fractions between different environments rather than a global suppression in the star formation rate for the star-forming population. However, when we restrict our sample to the cluster-scale environments ($M>10^{14}M_{solar}$), we find a global reduction in the SSFR of the star forming sequence of $17\%$ at 4$\sigma$ confidence as opposed to its field counterpart. After removing the stellar mass dependence of the quiescent fraction seen in field galaxies, the excess in the quiescent fraction due to the environment quenching in groups and clusters is found to increase with stellar mass. We argue that these results are in favor of galaxy mergers to be the primary environment quenching mechanism operating in galaxy groups whereas strangulation is able to reproduce the observed trend in the environment quenching efficiency and stellar mass relation seen in clusters. Our results also suggest that the relative importance between mass quenching and environment quenching depends on stellar mass -- the mass quenching plays a dominant role in producing quiescent galaxies for more massive galaxies, while less massive galaxies are quenched mostly through the environmental effect, with the transition mass around $1-2\times10^{10}M_{solar}$ in the group/cluster environment. (abridged)

Well-motivated particle physics theories predict the existence of particles (such as sterile neutrinos) which acquire non-negligible thermal velocities in the early universe. These particles could behave as warm dark matter (WDM) and generate a small-scale cutoff in the linear density power spectrum which scales approximately inversely with the particle mass. If this mass is of order a keV, the cutoff occurs on the scale of dwarf galaxies. Thus, in WDM models the abundance of small galaxies, such as the satellites that orbit in the halo of the Milky Way, depends on the mass of the warm particle. The abundance also scales with the mass of the host galactic halo. We use the \galform semi-analytic model of galaxy formation to calculate the properties of galaxies in universes in which the dark matter is warm. Using this method, we can compare the predicted satellite luminosity functions to the observed data for the Milky Way dwarf spheroidals, and determine a lower bound on the thermally produced WDM particle mass. This depends strongly on the value of the Milky Way halo mass and, to some extent, on the baryonic physics assumed; we examine both of these dependencies. For our fiducial model we find that for a particle mass of 3.3 keV (the 2$\sigma$ lower limit found by Viel et al. from a recent analysis of the Lyman-$\alpha$ forest) the Milky Way halo mass is required to be $> 1.4 \times 10^{12}$ \msun. For this same fiducial model, we also find that all WDM particle masses are ruled out (at 95% confidence) if the halo of the Milky Way has a mass smaller than $1.1 \times 10^{12}$ \msun, while if the mass of the Galactic halo is greater than 1.8 $\times 10^{12}$ \msun, only WDM particle masses larger than 2 keV are allowed.

The Pan-STARRS1 survey is acquiring multi-epoch 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 multi-epoch test area data, we present a detailed investigation into the depth of the survey and the reliability of the Pan-STARRS1 analysis software. We show that the Pan-STARRS1 reduction software can recover the properties of fake sources, and show good agreement between the magnitudes measured by Pan-STARRS1 and those from Sloan Digital Sky Survey. We also examine the number of false detections apparent in the Pan-STARRS1 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.

The Pan-STARRS1 survey is currently obtaining imaging in 5 bands (grizy) for the $3\pi$ steradian survey, one of the largest optical surveys ever conducted. The finished survey will have spatially varying depth, due to the survey strategy. This paper presents a method to correct galaxy number counts and galaxy clustering for this potential systematic based on a simplified signal to noise measurement. A star and galaxy separation method calibrated using realistic synthetic images is also presented, along with an approach to mask bright stars. By using our techniques on a ~69 sq. degree region of science verification data this paper shows PS1 measurements of the two point angular correlation function as a function of apparent magnitude agree with measurements from deeper, smaller surveys. Clustering measurements appear reliable down to a magnitude limit of rps<22.5. Additionally, stellar contamination and false detection issues are discussed and quantified. This work is the second of two papers which pave the way for the exploitation of the full $3\pi$ survey for studies of large scale structure.

We apply a recently developed scaling technique to the Millennium-XXL, one of the largest cosmological N-body 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 high-mass tail of the z=6 distribution. We assume these objects to be likely hosts for the population of rare but ultraluminous high-redshift 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 present-day 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 high-mass 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 high-z 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 broad-band spectral energy distributions (SEDs) with stellar population models, can be biased as a result of commonly adopted assumptions for the star-formation and chemical enrichment histories, recycled fractions and dust attenuation curves of galaxies. We apply the observational technique of broad-band 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 star-formation 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 star-dust 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, self-consistent way to compare model predictions with observations.

We constrain the distributions of projected radial alignment angles of satellite galaxy shapes within the Galaxy And Mass Assembly survey group catalogue. We identify the galaxy groups using spectroscopic redshifts and measure galaxy projected ellipticities from Sloan Digital Sky Survey imaging. With a sample of 3,850 groups with 13,655 satellite galaxies with high quality shape measurements, we find a less than 2-sigma signal of radial alignments in the mean projected ellipticity components and the projected position angle when using galaxy shape estimates optimized for weak lensing measurements. Our radial alignment measurement increases to greater than 3-sigma significance relative to the expectation for no alignments if we use 2-D S\'ersic model fits to define galaxy orientations. Our weak measurement of radial alignments is in conflict with predictions from dark matter N-body simulations, which we interpret as evidence for large mis-alignments of baryons and dark matter in group and cluster satellites. Within our uncertainties, that are dominated by our small sample size, we find only weak and marginally significant trends of the radial alignment angle distributions on projected distance from the group centre, host halo mass, and redshift that could be consistent with a tidal torquing mechanism for radial alignments. Using our lensing optimized shape estimators, we estimate that intrinsic alignments of galaxy group members may contribute a systematic error to the mean differential projected surface mass density of groups inferred from weak lensing observations by -1 +/- 20% at scales around 300 kpc/h from the group centre assuming a photometric redshift r.m.s. error of 10%, and given our group sample with median redshift of 0.17 and median virial masses ~10^{13} h^{-1}M_{sun}.

We use a suite of high resolution state-of-the-art N-body Dark Matter simulations of chameleon $f(R)$ gravity to study the higher order volume averaged correlation functions $\bar{\xi_n}$ together with the hierarchical $n$-th order correlation amplitudes $S_n=\bar{\xi}_n/\bar{\xi}_2^{n-1}$ and density distribution functions (PDF). We show that under the non-linear modifications of gravity the hierarchical scaling of the reduced cumulants is preserved. This is however characterised by significant changes of both the $\bar{\xi_n}$ and $S_n$'s values and their scale dependence with respect to General Relativity gravity (GR). In addition, we measure a significant increase of the non linear $\sigma_8$ parameter reaching $14, 5$ and 0.5% in excess of the GR value for the three flavours of our $f(R)$ models. We further note that the values of the reduced cumulants up to order $n=9$ are significantly increased in $f(R)$ gravity for all our models at small scales $R\simlt 30\hmpc$. In contrast the values of the hierarchical amplitudes, $S_n$'s, are smaller in $f(R)$ indicating that the modified gravity density distribution functions are deviating from the GR case. Furthermore we find that the redshift evolution of relative deviations of the $f(R)$ hierarchical correlation amplitudes is fastest at high and moderate redshifts $1\leq z \leq4$. The growth of these deviations significantly slows down in the low redshift universe. We also compute the PDFs and show that for scales below $\sim 20\hmpc$ they are significantly shifted in $f(R)$ gravity towards the low densities. Finally we discuss the implications of our theoretical predictions for measurements of the hierarchical clustering in galaxy redshifts surveys, including the important problems of the galaxy biasing and redshifts space distortions.

(Abridged) In tandem with observational datasets, we utilize realistic mock catalogs, based on a semi-analytic galaxy formation model, constructed specifically for Pan-STARRS1 Medium Deep Surveys in order to assess the performance of the Probability Friends-of-Friends (PFOF, Liu et al.) group finder, and aim to develop a grouping optimization method applicable to surveys like Pan-STARRS1. Producing mock PFOF group catalogs under a variety of photometric redshift accuracies ({\sigma}{\Delta}z/(1+zs)), we find that catalog purities and completenesses from ``good' {\sigma}{\Delta}z/(1+zs)) ~ 0.01) to ``poor' {\sigma}{\Delta}z/(1+zs)) ~ 0.07) photo-zs gradually degrade respectively from 77% and 70% to 52% and 47%. To avoid model dependency of the mock for use on observational data we apply a ``subset optimization' approach, using spectroscopic-redshift group data from the target field to train the group finder for application to that field, as an alternative method for the grouping optimization. We demonstrate this approach using these spectroscopically identified groups as the training set, i.e. zCOSMOS groups for PFOF searches within PS1 Medium Deep Field04 (PS1MD04) and DEEP2 EGS groups for searches in PS1MD07. We ultimately apply PFOF to four datasets spanning the photo-z uncertainty range from 0.01 to 0.06 in order to quantify the dependence of group recovery performance on photo-z accuracy. We find purities and completenesses calculated from observational datasets broadly agree with their mock analogues. Further tests of the PFOF algorithm are performed via matches to X-ray clusters identified within the PS1MD04 and COSMOS footprints. Across over a decade in group mass, we find PFOF groups match ~85% of X-ray clusters in COSMOS and PS1MD04, but at a lower statistical significance in the latter.

We place constraints on the matter density of the Universe and the amplitude of clustering using measurements of the galaxy two-point correlation function from the Sloan Digital Sky Survey (SDSS). We generate model predictions for different cosmologies by populating rescaled N-body simulations with galaxies using the subhalo abundance matching (SHAM) technique. We find Omega-M = 0.29 +/- 0.03 and sigma-8 = 0.86 +/- 0.04 at 68% confidence from fitting the observed two-point galaxy correlation function of galaxies brighter than M_r = -18 in a volume limited sample of galaxies obtained by the SDSS. We discuss and quantify potential sources of systematic error, and conclude that while there is scope for improving its robustness, the technique presented in this paper provides a powerful low redshift constraint on the cosmological parameters that is complementary to other commonly used methods.

The Galaxy And Mass Assembly (GAMA) survey is a multiwavelength photometric and spectroscopic survey, using the AAOmega spectrograph on the Anglo-Australian Telescope to obtain spectra for up to ~300000 galaxies over 280 square degrees, to a limiting magnitude of r_pet < 19.8 mag. The target galaxies are distributed over 0<z<0.5 with a median redshift of z~0.2, although the redshift distribution includes a small number of systems, primarily quasars, at higher redshifts, up to and beyond z=1. The redshift accuracy ranges from sigma_v~50km/s to sigma_v~100km/s depending on the signal-to-noise of the spectrum. Here we describe the GAMA spectroscopic reduction and analysis pipeline. We present the steps involved in taking the raw two-dimensional spectroscopic images through to flux-calibrated one-dimensional spectra. The resulting GAMA spectra cover an observed wavelength range of 3750<lambda<8850 A at a resolution of R~1300. The final flux calibration is typically accurate to 10-20%, although the reliability is worse at the extreme wavelength ends, and poorer in the blue than the red. We present details of the measurement of emission and absorption features in the GAMA spectra. These measurements are characterised through a variety of quality control analyses detailing the robustness and reliability of the measurements. We illustrate the quality of the measurements with a brief exploration of elementary emission line properties of the galaxies in the GAMA sample. We demonstrate the luminosity dependence of the Balmer decrement, consistent with previously published results, and explore further how Balmer decrement varies with galaxy mass and redshift. We also investigate the mass and redshift dependencies of the [NII]/Halpha vs [OIII]/Hbeta spectral diagnostic diagram, commonly used to discriminate between star forming and nuclear activity in galaxies.

We introduce a method for constructing end-to-end mock galaxy catalogues using a semi-analytical model of galaxy formation, applied to the halo merger trees extracted from a cosmological N-body simulation. The mocks that we construct are lightcone catalogues, in which a galaxy is placed according to the epoch at which it first enters the past lightcone of the observer, and incorporate the evolution of galaxy properties with cosmic time. We determine the position between the snapshot outputs at which a galaxy enters the observer's lightcone by interpolation. As an application, we consider the effectiveness of the BzK colour selection technique, which was designed to isolate galaxies in the redshift interval 1.4<z<2.5. The mock catalogue is in reasonable agreement with the observed number counts of all BzK galaxies, as well as with the observed counts of the subsample of BzKs that are star-forming galaxies. We predict that over 75 per cent of the model galaxies with K_{AB}<=23, and 1.4<z<2.5, are selected by the BzK technique. Interloper galaxies, outside the intended redshift range, are predicted to dominate bright samples of BzK galaxies (i.e. with K_{AB}<=21). Fainter K-band cuts are necessary to reduce the predicted interloper fraction. We also show that shallow B-band photometry can lead to confusion in classifying BzK galaxies as being star-forming or passively evolving. Overall, we conclude that the BzK colour selection technique is capable of providing a sample of galaxies that is representative of the 1.4<z<2.5 galaxy population.

We investigate the luminosity functions (LFs) and projected number density profiles of galactic satellites around isolated primaries of different luminosities. We measure these quantities for model satellites placed into the Millennium and Millennium II dark matter simulations by the GALFORM semi-analytic galaxy formation model for different bins of primary galaxy magnitude and we investigate their dependence on satellite luminosity. We compare our model predictions to the data of Guo et al. from the Sloan Digital Sky Survey Data Release 8 (SDSS DR8). First, we use a mock light-cone catalogue to verify that the method we used to count satellites in the SDSS DR8 is unbiased. We find that the radial distributions of model satellites are similar to those around comparable primary galaxies in the SDSS DR8, with only slight differences at low luminosities and small projected radii. However, when splitting the satellites by colour, the model and SDSS satellite systems no longer resemble one another, with many red model satellites, in contrast to the dominant blue fraction at similar luminosity in SDSS. The few model blue satellites are also significantly less centrally concentrated in the halo of their stacked primary than their SDSS counterparts. The implications of this result for the GALFORM model are discussed.

We describe a methodology to accurately compute halo mass functions, progenitor mass functions, merger rates and merger trees in non-cold dark matter universes using a self-consistent treatment of the generalized extended Press-Schechter formalism. Our approach permits rapid exploration of the subhalo population of galactic halos in dark matter models with a variety of different particle properties or universes with rolling, truncated, or more complicated power spectra. We make detailed comparisons of analytically derived mass functions and merger histories with recent warm dark matter cosmological N-body simulations, and find excellent agreement. We show that, once the accretion of smoothly distributed matter is accounted for, coarse-grained statistics such as the mass accretion history of halos can be almost indistinguishable between cold and warm dark matter cases. However, the halo mass function and progenitor mass functions differ significantly, with the warm dark matter cases being strongly suppressed below the free-streaming scale of the dark matter. We demonstrate the importance of using the correct solution for the excursion set barrier first-crossing distribution in warm dark matter - if the solution for a flat barrier is used instead the truncation of the halo mass function is much slower, leading to an overestimate of the number of low mass halos.

Star-forming galaxies which are too faint to be detected individually produce intensity fluctuations in the cosmic background light. This contribution needs to be taken into account as a foreground when using the primordial signal to constrain cosmological parameters. The extragalactic fluctuations are also interesting in their own right as they depend on the star formation history of the Universe and the way in which this connects with the formation of cosmic structure. We present a new framework which allows us to predict the occupation of dark matter haloes by star-forming galaxies and uses this information, in conjunction with an N-body simulation of structure formation, to predict the power spectrum of intensity fluctuations in the infrared background. We compute the emission from galaxies at far-infrared, millimetre and radio wavelengths. Our method gives accurate predictions for the clustering of galaxies both for the one halo and two halo terms. We illustrate our new framework using a previously published model which reproduces the number counts and redshift distribution of galaxies selected by their emission at $850\,\mu$m. Without adjusting any of the model parameters, the predictions show encouraging agreement at high frequencies and on small angular scales with recent estimates of the extragalactic fluctuations in the background made from early data analysed by the Planck Collaboration. There are, however, substantial discrepancies between the model predictions and observations on large angular scales and at low frequencies, which illustrates the usefulness of the intensity fluctuations as a constraint on galaxy formation models.

We study the spatial distribution of satellite galaxies around isolated primaries using the Sloan Digital Sky Survey (SDSS) spectroscopic and photometric galaxy catalogues. We select isolated primaries from the spectroscopic sample and search for potential satellites in the much deeper photometric sample. For specific luminosity primaries we obtain robust statistical results by stacking as many as ~50, 000 galaxy systems. We find no evidence for any anistropy in the satellite galaxy distribution relative to the major axes of the primaries. We derive accurate projected number density profiles of satellites down to 4 magnitudes fainter than their primaries. We find the normalized satellite profiles generally have a universal form and can be well fitted by projected NFW profiles. The NFW concentration parameter increases with decreasing satellite luminosity while being independent of the luminosity of the primary except for very bright primaries. The profiles of the faintest satellites show deviations from the NFW form with an excess at small galactocentric projected distances. In addition, we quantify how the radial distribution of satellites depends on the colour of the satellites and on the colour and concentration of their primaries.

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 1st-order 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 energy-based constraints on the outflow. Using a trial exit velocity distribution, this is shown to be compatible with the observed baryon fractions in intermediate-mass systems, but implies that current estimates for low-mass systems can not be solely accounted for by supernova winds under commonly-held assumptions.

We present measurements of the triaxial dark matter halo shapes and alignment correlation functions in the Millennium and Millennium-2 dark matter N-body simulations. These two simulations allow us to measure the distributions of halo shapes down to 10% of the virial radius over a halo mass range of 6E9 - 2E14 M_sun/h. We largely confirm previous results on the distributions of halo axis ratios as a function of halo mass, but we find that the median angle between halo major axes at different halo radii can vary by a factor of 2 between the Millennium-1 and 2 simulations because of the different mass resolution. Thus, error in the shape determinations from limited resolution is potentially degenerate with the misalignment of halo inner and outer shapes used to constrain Brightest Cluster Galaxy alignments in previous works. We also present simplifying parameterizations for the 3-D halo-mass alignment correlation functions that are necessary ingredients for triaxial halo models of large-scale structure and models of galaxy intrinsic alignments as contaminants for cosmic shear surveys. We measure strong alignments between halos of all masses and the surrounding dark matter overdensities out to several tens of Mpc/h, in agreement with observed shear-galaxy and cluster shape correlations. We use these measurements to forecast the contribution to the weak lensing signal around galaxy clusters from correlated mass along the line-of-sight. For prolate clusters with major axes aligned with the line-of-sight the fraction of the weak lensing signal from mass external to the cluster can be twice that predicted if the excess halo alignment correlation is assumed to be zero.

We present new deep observations of 'shell' structures in the halo of the nearby elliptical galaxy NGC 7600, alongside a movie of galaxy formation in a cold dark matter universe (available at http://www.virgo.dur.ac.uk/shell-galaxies). The movie, based on an ab initio cosmological simulation, shows how continuous accretion of clumps of dark matter and stars creates a swath of diffuse circumgalactic structures. The disruption of a massive clump on a near-radial orbit creates a complex system of transient concentric shells which bare a striking resemblance to those of NGC 7600. With the aid of the simulation we interpret NGC 7600 in the context of the CDM model.

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 N-body 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 turn-off 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 in-situ.

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 high-resolution N-body simulations in combination with the Galform semi-analytic 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 present a new algorithm to generate a random (unclustered) version of an magnitude limited observational galaxy redshift catalogue. It takes into account both galaxy evolution and the perturbing effects of large scale structure. The key to the algorithm is a maximum likelihood (ML) method for jointly estimating both the luminosity function (LF) and the overdensity as a function of redshift. The random catalogue algorithm then works by cloning each galaxy in the original catalogue, with the number of clones determined by the ML solution. Each of these cloned galaxies is then assigned a random redshift uniformly distributed over the accessible survey volume, taking account of the survey magnitude limit(s) and, optionally, both luminosity and number density evolution. The resulting random catalogues, which can be employed in traditional estimates of galaxy clustering, make fuller use of the information available in the original catalogue and hence are superior to simply fitting a functional form to the observed redshift distribution. They are particularly well suited to studies of the dependence of galaxy clustering on galaxy properties as each galaxy in the random catalogue has the same list of attributes as measured for the galaxies in the genuine catalogue. The derivation of the joint overdensity and LF estimator reveals the limit in which the ML estimate reduces to the standard 1/Vmax LF estimate, namely when one makes the prior assumption that the are no fluctuations in the radial overdensity. The new ML estimator can be viewed as a generalization of the 1/Vmax estimate in which Vmax is replaced by a density corrected Vdc,max.

We study the luminosity function of satellite galaxies around isolated primaries using the Sloan Digital Sky Survey (SDSS) spectroscopic and photometric galaxy samples. We select isolated primaries from the spectroscopic sample and search for potential satellites in the much deeper photometric sample. For primaries of similar luminosity to the Milky Way and M31, we are able to stack as many as $\sim\negthinspace 20,000$ galaxy systems to obtain robust statistical results. We derive the satellite luminosity function extending almost 8 magnitudes fainter than the primary galaxy. We also determine how the satellite luminosity function varies with the luminosity, colour and concentration of the primary. We find that, in the mean, isolated primaries of comparable luminosity to the Milky Way and M31 contain about a factor of two fewer satellites brighter than $M_V=-14$ than the average of the Milky Way and M31.

We present an algorithm for quickly generating multiple realizations of N-body simulations to be used, for example, for cosmological parameter estimation from surveys of large-scale structure. Our algorithm uses a new method to resample the large-scale (Gaussian-distributed) Fourier modes in a periodic N-body simulation box in a manner that properly accounts for the nonlinear mode-coupling between large and small scales. We find that our method for adding new large-scale mode realizations recovers the nonlinear power spectrum to sub-percent accuracy on scales larger than about half the Nyquist frequency of the simulation box. Using 20 N-body simulations, we obtain a power spectrum covariance matrix estimate that matches the estimator in Takahashi et al. (2009) (from 5000 simulations) with < 20% errors in all matrix elements. Comparing the rates of convergence, we determine that our algorithm requires ~ 8 times fewer simulations to achieve a given error tolerance in estimates of the power spectrum covariance matrix. The degree of success of our algorithm indicates that we understand the main physical processes that give rise to the correlations in the matter power spectrum. Namely, the large-scale Fourier modes modulate both the degree of structure growth through the variation in the effective local matter density and also the spatial frequency of small-scale perturbations through large-scale displacements. We expect our algorithm to be useful for noise modeling when constraining cosmological parameters from weak lensing (cosmic shear) and galaxy surveys, rescaling summary statistics of N-body simulations for new cosmological parameter values, and any applications where the influence of Fourier modes larger than the simulation size must be accounted for.

Star formation rate and accummulated stellar mass are two fundamental physical quantities that describe the evolutionary state of a forming galaxy. Two recent attempts to determine the relationship between these quantities, by interpreting a sample of star-forming galaxies at redshift of z~4, have led to opposite conclusions. We use a model galaxy population to investigate possible causes for this discrepancy and conclude that minor errors in the conversion from observables to physical quantities can lead to major misrepresentation when applied without awareness of sample selection. We also investigate, in a general way, the physical origin of the correlation between star formation rate and stellar mass within hierarchical galaxy formation theory.

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 quasi-hydrostatic hot gas haloes (called the hot-halo 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 X-ray, 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 hot-halo accretion modes. The faint end of the LF is dominated by massive BHs experiencing quiescent accretion via a thick disc, primarily during the hot-halo 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~2-3). This implies that the cosmic evolution of AGN is shaped in part by obscuration.

We measure the projected cross-correlation between low redshift (z < 0.5) far-IR selected galaxies in the SDP field of the Herschel-ATLAS (H-ATLAS) survey and optically selected galaxies from the Galaxy and Mass Assembly (GAMA) redshift survey. In order to obtain robust correlation functions, we restrict the analysis to a subset of 969 out of 6900 H-ATLAS galaxies, which have reliable optical counterparts with r<19.4 mag and well-determined spectroscopic redshifts. The overlap region between the two surveys is 12.6 sq. deg; the matched sample has a median redshift of z ~ 0.2. The cross-correlation of GAMA and H-ATLAS galaxies within this region can be fitted by a power law, with correlation length r_0 ~ 4.63 +/- 0.51 Mpc. Comparing with the corresponding auto-correlation function of GAMA galaxies within the SDP field yields a relative bias (averaged over 2-8 Mpc) of H-ATLAS and GAMA galaxies of b_H/b_G ~ 0.6. Combined with clustering measurements from previous optical studies, this indicates that most of the low redshift H-ATLAS sources are hosted by halos with masses comparable to that of the Milky Way. The correlation function appears to depend on the 250 um luminosity, L_250, with bright (median luminosity \nu L_250 ~ 1.6 x 10^10 L_sun) objects being somewhat more strongly clustered than faint (\nu L_250 ~ 4.0 x 10^9 L_sun) objects. This implies that galaxies with higher dust-obscured star formation rates are hosted by more massive halos.

We track the co-evolution of supermassive black holes (SMBHs) and their host galaxies. The calculation is embedded in the GALFORM semi-analytical 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 self-gravity, 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 Blandford-Znajek 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.

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.

Semi-analytic models are a powerful tool for studying the formation of galaxies. However, these models inevitably involve a significant number of poorly constrained parameters that must be adjusted to provide an acceptable match to the observed universe. In this paper, we set out to quantify the degree to which observational data-sets can constrain the model parameters. By revealing degeneracies in the parameter space we can hope to better understand the key physical processes probed by the data. We use novel mathematical techniques to explore the parameter space of the GALFORM semi-analytic model. We base our investigation on the Bower et al. 2006 version of GALFORM, adopting the same methodology of selecting model parameters based on an acceptable match to the local bJ and K luminosity functions. The model contains 16 parameters that are poorly constrained, and we investigate this parameter space using the Model Emulator technique, constructing a Bayesian approximation to the GALFORM model that can be rapidly evaluated at any point in parameter space. By combining successive waves of emulation, we show that only 0.26% of the initial volume is of interest for further exploration. However, within this region we show that the Bower et al. 2006 model is only one choice from an extended sub-space of model parameters that can provide equally acceptable fits. We explore the geometry of this region and begin to explore the physical connections between parameters that are exposed by this analysis. We also consider the impact of adding additional observational data to further constrain the parameter space.

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 N-body simulations of individual dark matter haloes. We tag subsets of particles in these simulations with stellar populations predicted by the Galform semi-analytic model. Our method self-consistently 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, metal-rich 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 self-bound 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/surface-brightness profiles of the stellar haloes are consistent with expectations from the Milky Way and M31. Each halo has a complex structure, consisting of well-mixed 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 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 N-body simulations are used to calculate the complete non-linear 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 Sachs-Wolfe or ISW effect) and non-linear (the Rees-Sciama or RS effect) evolution. These effects give rise to three kinds of non-linear 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 non-linear 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 non-Gaussianity, skewing the one-point 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 non-linear, non-Gaussian features. Our maps are available at http://astro.dur.ac.uk/~cai/ISW

The distribution of cold gas in dark matter haloes is driven by key processes in galaxy formation: gas cooling, galaxy mergers, star formation and reheating of gas by supernovae. We compare the predictions of four different galaxy formation models for the spatial distribution of cold gas. We find that satellite galaxies make little contribution to the abundance or clustering strength of cold gas selected samples, and are far less important than they are in optically selected samples. The halo occupation distribution function of present-day central galaxies with cold gas mass > 10^9 h^-1 Msun is peaked around a halo mass of ~ 10^11 h^-1 Msun, a scale that is set by the AGN suppression of gas cooling. The model predictions for the projected correlation function are in good agreement with measurements from the HI Parkes All-Sky Survey. We compare the effective volume of possible surveys with the Square Kilometre Array with those expected for a redshift survey in the near-infrared. Future redshift surveys using neutral hydrogen emission will be competitive with the most ambitious spectroscopic surveys planned in the near-infrared.

Recent analysis of the Milky Way's satellite galaxies reveals that these objects share a common central mass density, even though their luminosities range over five orders of magnitude. This observation can be understood in the context of galaxy formation theory by quantifying the factors which restrict the central mass density to a small range. One limit is set by the maximum mass that can collapse into a given region by the hierarchical growth of structure in the standard cold dark matter cosmology. Another limit comes from the natural thresholds which exist for gas to be able to cool and form a galaxy. The wide range of luminosities in these satellites reflect the effects of supernova feedback on the fraction of cooled baryons which are retained.

The Galaxy And Mass Assembly (GAMA) project is the latest in a tradition of large galaxy redshift surveys, and is now underway on the 3.9m Anglo-Australian Telescope at Siding Spring Observatory. GAMA is designed to map extragalactic structures on scales of 1kpc - 1Mpc in complete detail to a redshift of z~0.2, and to trace the distribution of luminous galaxies out to z~0.5. The principal science aim is to test the standard hierarchical structure formation paradigm of Cold Dark Matter (CDM) on scales of galaxy groups, pairs, discs, bulges and bars. We will measure (1) the Dark Matter Halo Mass Function (as inferred from galaxy group velocity dispersions); (2) baryonic processes, such as star formation and galaxy formation efficiency (as derived from Galaxy Stellar Mass Functions); and (3) the evolution of galaxy merger rates (via galaxy close pairs and galaxy asymmetries). Additionally, GAMA will form the central part of a new galaxy database, which aims to contain 275,000 galaxies with multi-wavelength coverage from coordinated observations with the latest international ground- and space-based facilities: GALEX, VST, VISTA, WISE, HERSCHEL, GMRT and ASKAP. Together, these data will provide increased depth (over 2 magnitudes), doubled spatial resolution (0.7"), and significantly extended wavelength coverage (UV through Far-IR to radio) over the main SDSS spectroscopic survey for five regions, each of around 50 deg^2. This database will permit detailed investigations of the structural, chemical, and dynamical properties of all galaxy types, across all environments, and over a 5 billion year timeline.

The clustering amplitude of galaxies depends on their intrinsic luminosity. We compare the properties of publicly available galaxy formation models with clustering measurements from the two-degree field galaxy redshift survey. The model predictions show the same qualitative behaviour as the data but fail to match the observations at the level of accuracy at which current measurements can be made. We demonstrate that this is due to the model producing too many satellite galaxies in massive haloes. We implement simple models to describe two new processes, satellite-satellite mergers and the tidal dissolution of satellites to investigate their impact on the predicted clustering. We find that both processes need to be included in order to produce a model which matches the observations.

The highly flattened distribution of satellite galaxies in the Milky Way presents a number of puzzles. Firstly, its polar alignment stands out from the planar alignments commonly found in other galaxies. Secondly, recent proper motion measurements reveal that the orbital angular momentum of at least 3, and possibly as many as 8, of the Milky Ways satellites point (within 30 degrees) along the axis of their flattened configuration, suggesting some form of coherent motion. In this paper we use a high resolution cosmological simulation to investigate whether this pattern conflicts with the expectations of the cold dark matter model of structure formation. We find that this seemingly unlikely set up occurs often: approximately 35% of the time we find systems in which the angular momentum of 3 individual satellites point along, or close to, the short axis of the satellite distribution. In addition, in 30% of the systems we find that the net angular momentum of the 6 best aligned satellites lies within 35 degrees of the short axis of the satellite distribution, as observed for the Milky Way.

In a universe with a cosmological constant, the large-scale gravitational potential varies in time and this is, in principle, observable. Using an N-body 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 large-scale structure on the cosmic microwave background (CMB). The linear part of the $\dot{\Phi}$ power spectrum (the integrated Sachs-Wolfe effect or ISW) drops quickly as the relative importance of $\Omega_{\Lambda}$ diminishes at high redshift, while the non-linear part (the Rees-Sciama 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 cross-correlation 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 non-linear effect dominates. For the cross-correlation power spectrum of galaxy samples with the CMB, the non-linear effect becomes significant at $l\sim 500$ and rapidly makes the cross power spectrum negative. For high redshift samples, the cross-correlation is expected to be suppressed by $5-10%$ on arcminute scales. The RS effect makes a negligible contribution to the large-scale ISW cross-correlation measurement. However, on arc-minute scales it will contaminate the expected cross-correlation signal induced by the Sunyaev-Zel'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 semi-analytic 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 Pan-STARRS system. PS1 has 5 photometric bands (grizy), and will carry out an all-sky 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 off-the-shelf photo-z 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.02-0.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 photo-z 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.

Current models of galaxy formation predict satellite galaxies in groups and clusters that are redder than observed. We investigate the effect on the colours of satellite galaxies produced by the ram pressure stripping of their hot gaseous atmospheres as the satellites orbit within their parent halo. We incorporate a model of the stripping process based on detailed hydrodynamic simulations within the Durham semi-analytic model of galaxy formation. The simulations show that the environment in groups and clusters is less aggressive than previously assumed. The main uncertainty in the model is the treatment of gas expelled by supernovae. With reasonable assumptions for the stripping of this material, we find that satellite galaxies are able to retain a significant fraction of their hot gas for several Gigayears, thereby replenishing their reservoirs of cold, star forming gas and remaining blue for a relatively long period of time. A bimodal distribution of galaxy colours, similar to that observed in SDSS data, is established and the colours of the satellite galaxies are in good agreement with the data. In addition, our model naturally accounts for the observed dependence of satellite colours on environment, from small groups to high mass clusters.

(Abridged) We use two very large cosmological simulations to study how the density profiles of relaxed LCDM dark halos depend on redshift and on halo mass. We confirm that these profiles deviate slightly but systematically from the NFW form and are better approximated by the empirical formula, $d\log \rho/d\log r \propto r^{\alpha}$, first used by Einasto to fit star counts in the Milky Way. The best-fit value of the additional shape parameter, $\alpha$, increases gradually with mass, from alpha ~ 0.16 for present-day galaxy halos to alpha~ 0.3 for the rarest and most massive clusters. Halo concentrations depend only weakly on mass at z=0, and this dependence weakens further at earlier times. At z ~3 the average concentration of relaxed halos does not vary appreciably over the mass range accessible to our simulations (M >=3e11Msun). Furthermore, in our biggest simulation, the average concentration of the most massive}, relaxed halos is constant at c200 ~3.5 to 4 for 0<=z<= 3. These results have important implications for currently planned surveys of distant clusters.

Published galaxy power spectra from the 2dFGRS and SDSS are not in good agreement. We revisit this issue by analyzing both the 2dFGRS and SDSS DR5 catalogues using essentially identical techniques. We confirm that the 2dFGRS exhibits relatively more large scale power than the SDSS, or, equivalently, SDSS has more small scale power. We demonstrate that this difference is due to the r-band selected SDSS catalogue being dominated by more strongly clustered red galaxies, which have a stronger scale dependent bias. The power spectra of galaxies of the same rest frame colours from the two surveys match well. If not accounted for, the difference between the SDSS and 2dFGRS power spectra causes a bias in the obtained constraints on cosmological parameters which is larger than the uncertainty with which they are determined. We also found that the correction developed by Cole et al.(2005) to model the distortion in the shape of the power spectrum due to non-linear evolution and scale dependent bias is not able to reconcile the constraints obtained from the 2dFGRS and SDSS power spectra. Intriguingly, the model is able to describe the differences between the 2dFGRS and the much more strongly clustered LRG sample, which exhibits greater nonlinearities. This shows that more work is needed to understand the relation between the galaxy power spectrum and the linear perturbation theory prediction for the power spectrum of matter fluctuations. It is therefore important to accurately model these effects to get precise estimates of cosmological parameters from these power spectra and from future galaxy surveys like Pan-STARRS, or the Dark Energy Survey, which will use selection criteria similar to the one of SDSS.

We present a new Monte-Carlo algorithm to generate merger trees describing the formation history of dark matter halos. The algorithm is a modification of the algorithm of Cole et al (2000) used in the GALFORM semi-analytic galaxy formation model. As such, it is based on the Extended Press-Schechter theory and so should be applicable to hierarchical models with a wide range of power spectra and cosmological models. It is tuned to be in accurate agreement with the conditional mass functions found in the analysis of merger trees extracted from the LCDM Millennium N-body simulation. We present a comparison of its predictions not only with these conditional mass functions, but also with additional statistics of the Millennium Simulation halo merger histories. In all cases we find it to be in good agreement with the Millennium Simulation and thus it should prove to be a very useful tool for semi-analytic models of galaxy formation and for modelling hierarchical structure formation in general. We have made our merger tree generation code and code to navigate the trees available at http://star-www.dur.ac.uk/~cole/merger_trees .

We identify a large sample of isolated bright galaxies and their fainter satellites in the 2dF Galaxy Redshift Survey (2dFGRS). We analyse the dynamics of ensembles of these galaxies selected according to luminosity and morphological type by stacking the positions of their satellites and estimating the velocity dispersion of the combined set. We test our methodology using realistic mock catalogues constructed from cosmological simulations. The method returns an unbiased estimate of the velocity dispersion provided that the isolation criterion is strict enough to avoid contamination and that the scatter in halo mass at fixed primary luminosity is small. Using a maximum likelihood estimator that accounts for interlopers, we determine the satellite velocity dispersion within a projected radius of 175 kpc/h. The dispersion increases with the luminosity of the primary and is larger for elliptical galaxies than for spiral galaxies of similar bJ luminosity. Calibrating the mass-velocity dispersion relation using our mock catalogues, we find a dynamical mass within 175 kpc/h of M_175 ~ 4.0^{+2.3}_{-1.5} 10^12 (L_bJ/L_*) M_sol/h for elliptical galaxies and M_175 ~ 6.3^{+6.3}_{-3.1} 10^11 (L_bJ/L_*)^1.6 Msol/h for spiral galaxies. Finally, we compare our results with recent studies and investigate their limitations using our mock catalogues.

We introduce a method to constrain general cosmological models using Baryon Acoustic Oscillation (BAO) distance measurements from galaxy samples covering different redshift ranges, and apply this method to analyse samples drawn from the SDSS and 2dFGRS. BAO are detected in the clustering of the combined 2dFGRS and SDSS main galaxy samples, and measure the distance--redshift relation at z=0.2. BAO in the clustering of the SDSS luminous red galaxies measure the distance--redshift relation at z=0.35. The observed scale of the BAO calculated from these samples and from the combined sample are jointly analysed using estimates of the correlated errors, to constrain the form of the distance measure D_V(z)=[(1+z)^2D_A^2cz/H(z)]^(1/3). Here D_A is the angular diameter distance, and H(z) is the Hubble parameter. This gives r_s/D_V(0.2)=0.1980+/-0.0058 and r_s/D_V(0.35)=0.1094+/-0.0033 (1sigma errors), with correlation coefficient of 0.39, where r_s is the comoving sound horizon scale at recombination. Matching the BAO to have the same measured scale at all redshifts then gives D_V(0.35)/D_V(0.2)=1.812+/-0.060. The recovered ratio is roughly consistent with that predicted by the higher redshift SNLS supernovae data for Lambda cosmologies, but does require slightly stronger cosmological acceleration at low redshift. If we force the cosmological model to be flat with constant w, then we find Om_m=0.249+/-0.018 and w=-1.004+/-0.089 after combining with the SNLS data, and including the WMAP measurement of the apparent acoustic horizon angle in the CMB.

We present a comparison of the statistical properties of dark matter halo merger trees extracted from the Millennium Simulation with Extended Press-Schechter (EPS) formalism and the related GALFORM Monte-Carlo method for generating ensembles of merger trees. The volume, mass resolution and output frequency make the Millennium Simulation a unique resource for the study of the hierarchical growth of structure. We construct the merger trees of present day friends-of-friends groups and calculate a variety of statistics that quantify the masses of their progenitors as a function of redshift; accretion rates; and the redshift distribution of their most recent major merger. We also look in the forward direction and quantify the present day mass distribution of halos into which high redshift progenitors of a specific mass become incorporated. We find that EPS formalism and its Monte-Carlo extension capture the qualitative behaviour of all these statistics but, as redshift increases they systematically underestimate the masses of the most massive progenitors. This shortcoming is worst for the Monte-Carlo algorithm. We present a fitting function to a scaled version of the progenitor mass distribution and show how it can be used to make more accurate predictions of both progenitor and final halo mass distributions.

We generate mock galaxy catalogues for a grid of different cosmologies, using rescaled N-body simulations in tandem with a semi-analytic model run using consistent parameters. Because we predict the galaxy bias, rather than fitting it as a nuisance parameter, we obtain an almost pure constraint on sigma_8 by comparing the projected two-point correlation function we obtain to that from the SDSS. A systematic error arises because different semi-analytic modelling assumptions allow us to fit the r-band luminosity function equally well. Combining our estimate of the error from this source with the statistical error, we find sigma_8=0.97 +/- 0.06. We obtain consistent results if we use galaxy samples with a different magnitude threshold, or if we select galaxies by b_J-band rather than r-band luminosity and compare to data from the 2dFGRS. Our estimate for sigma_8 is higher than that obtained for other analyses of galaxy data alone, and we attempt to find the source of this difference. We note that in any case, galaxy clustering data provide a very stringent constraint on galaxy formation models.

We use the Millennium Simulation (MS) to study the statistics of LCDM halo concentrations at z = 0. Our results confirm that the average halo concentration declines monotonically with mass; a power-law fits well the concentration-mass relation for over 3 decades in mass, up to the most massive objects to form in a LCDM universe (~ 10^15 h^-1 Msol). This is in clear disagreement with the predictions of the model proposed by Bullock et al. for these rare objects, and agrees better with the original predictions of Navarro, Frenk, & White. The large volume surveyed, together with the unprecedented numerical resolution of the MS, allow us to estimate with confidence the distribution of concentrations and, consequently, the abundance of systems with unusual properties. About one in a hundred cluster haloes (M200 >~ 3x10^14 h^-1 Msol) have concentrations exceeding c200 = 7.5, a result that may be used to interpret the likelihood of unusually strong massive gravitational lenses, such as Abell 1689, in the LCDM cosmogony. A similar fraction (1 in 100) of galaxy-sized haloes (M200 ~ 10^12 h^-1 Msol) have c200 < 4.5, an important constraint on models that attempt to reconcile the rotation curves of low surface-brightness galaxies by appealing to haloes of unexpectedly low concentration. We find that halo concentrations are independent of spin once haloes manifestly out of equilibrium are removed from the sample. Compared to their relaxed brethren, the concentrations of out-of-equilibrium haloes tend to be lower and to have more scatter, while their spins tend to be higher... (continue)

Published galaxy power spectra from the 2dFGRS and SDSS are not in good agreement. We revisit this issue by analyzing both the 2dFGRS and SDSS DR5 catalogues using essentially identical technqiues. We confirm that the 2dFGRS exhibits relatively more large scale power than the SDSS, or, equivalently, SDSS has more small scale power. We demonstrate that this difference is due the r-band selected SDSS catalogue being dominated by more strongly clustered red galaxies, due to these galaxies having a stronger scale dependent bias. The power spectra of galaxies of the same rest frame colours from the two surveys match well. It is therefore important to accurately model scale dependent bias to get accurate estimates of cosmological parameters from these power spectra.

We investigate the properties of satellite galaxies formed in N-body/SPH simulations of galaxy formation in the $\Lambda$CDM cosmology. The simulations include the main physical effects thought to be important in galaxy formation and, in several cases, produce realistic spiral discs. In total, a sample of 9 galaxies of luminosity comparable to the Milky Way was obtained. At magnitudes brighter than the resolution limit, $M_V=-12$, the luminosity function of the satellite galaxies in the simulations is in excellent agreement with data for the Local Group. The radial number density profile of the model satellites, as well as their gas fractions also match observations very well. In agreement with previous N-body studies, we find that the satellites tend to be distributed in highly flattened configurations whose major axis is aligned with the major axis of the (generally triaxial) dark halo. In 2 out of 3 systems with sufficiently large satellite populations, the satellite system is nearly perpendicular to the plane of the galactic disc, a configuration analogous to that observed in the Milk Way. The discs themselves are perpendicular to the minor axis of their host halos in the inner parts, and the correlation between the orientation of the galaxy and the shape of the halo persists even out to the virial radius. However, in one case the disc's minor axis ends up, at the virial radius, perpendicular to the minor axis of the halo. The angular momenta of the galaxies and their host halo tend to be well aligned.

Recent observations of the distant Universe suggest that much of the stellar mass of bright galaxies was already in place at $z>1$. This presents a challenge for models of galaxy formation because massive halos are assembled late in hierarchical cosmologies such as cold dark matter (CDM). In this paper, we discuss a new implementation of the Durham semi-analytic model in which feedback due to active galactic nuclei (AGN) is assumed to quench cooling flows in massive halos. This mechanism naturally creates a break in the local galaxy luminosity function at bright magnitudes. The model is implemented within the Millennium N-body simulation; the accurate dark matter merger trees and large number of realizations of the galaxy formation process that the simulation provides results in highly accurate statistics. After adjusting the values of the physical parameters in the model by reference to the properties of local galaxies, we use it to investigate the evolution of the K-band luminosity and galaxy stellar mass functions. We also calculate the volume averaged star formation rate density of the Universe as a function of redshift and the way in which this is apportioned amongst galaxies of different mass. The model robustly predicts a substantial population of massive galaxies out to redshift $z\sim 5$ and a star formation rate density which rises with increasing redshift in objects of all masses. Although observational data on these properties have been cited as evidence for ``anti-hierarchical'' galaxy formation, we find that when AGN feedback is taken into account, the fundamentally hierachical CDM model provides a very good match to these observations.

We study the environmental dependence of the formation epoch of dark matter haloes in the Millennium Simulation: a ten billion particle N-body simulation of standard LambdaCDM cosmology. A sensitive test of this dependence - the marked correlation function - reveals highly significant evidence that haloes of a given mass form earlier in denser regions. We define a marked cross-correlation function, which helps quantify how this effect depends upon the choice of the halo population used to define the environment. The mean halo formation redshift as a function of the local overdensity in dark matter is also well determined, and we see an especially clear dependence for galaxy-sized haloes. This contradicts one of the basic predictions of the excursion set model of structure formation, even though we see that this theory predicts other features of the distribution of halo formation epochs rather well. It also invalidates an assumption usually employed in the popular halo, or HOD, models of galaxy clustering, namely that the distribution of halo properties is a function of halo mass but not of halo environment.

We derive constraints on cosmological parameters using the power spectrum of galaxy clustering measured from the final two-degree field galaxy redshift survey (2dFGRS) and a compilation of measurements of the temperature power spectrum and temperature-polarization cross-correlation 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 energy-density 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.

We present a hierarchical galaxy formation model which can account for the number counts of sources detected through their emission at sub-millimetre wavelengths. The first stage in our approach is an ab initio calculation of the star formation histories for a representative sample of galaxies, which is carried out using the semi-analytical galaxy formation model GALFORM. These star formation histories are then input into the spectro-photometric code GRASIL, to produce a spectral energy distribution for each galaxy. Dust extinction and emission are treated self consistently in our model, without having to resort to ad-hoc assumptions about the amount of attenuation by dust or the temperature at which the dust radiates. We argue that it is necessary to modify the form of the stellar initial mass function in starbursts in order to match the observed number of sub-mm sources, if we are to retain the previous good matches enjoyed between observations and model predictions in the local universe. We also list some other observational tests that have been passed by our model.

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 semi-analytic 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 high-resolution N-body simulations of Milk-Way mass halos. We find that intergalactic black holes make up only ~2-3% 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 X-ray sources often seen near the centres of galaxies and for the hyperluminous non-central X-ray source in M82.

The abundance of galaxy clusters as a function of mass is determined using the 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue. This is used to estimate the amplitude of the matter fluctuation spectrum, parametrised by the linear theory rms density fluctuations in spheres of 8Mpc/h, sigma_8. The best-fitting value for this parameter is highly correlated with the mean matter density in the Universe, Omega_m, and is found to satisfy sigma_8=0.25Omega_m^{-0.92-4.5(Omega_m-0.22)^2}+/-10%(statistical)+/-20% (systematic) for 0.18<Omega_m<0.50, assuming that Omega_m+Omega_Lambda=1. This gives sigma_8=0.89 when Omega_m=0.25. A ~20% correction has been applied to undo the systematic bias inherent in the measurement procedure. Mock catalogues, constructed from large cosmological N-body simulations, are used to help understand and model these systematic errors. The abundance of galaxy groups as a function of group bj band luminosity is also determined. This is used in conjunction with the halo mass function, determined from simulations, to infer the variation of halo mass-to-light ratio over four orders of magnitude in halo mass. The mass-to-light ratio shows a minimum value of 100hM_sol/L_sol in the bj band at a total group luminosity of L_bJ ~ 5.10^9Lsol/h^2. Together with the observed Tully-Fisher relation, this implies that the observed rotation speed of Tully-Fisher galaxies is within ~10% of the typical circular speed of haloes hosting brightest galaxies of the same luminosity.

We present new results from a multi-wavelength model of galaxy formation, which combines a semi-analytical treatment of the formation of galaxies within the CDM framework with a sophisticated treatment of absorption and emission of radiation by dust. We find that the model, which incorporates a top-heavy IMF in bursts, agrees well with the evolution of the rest-frame far-UV luminosity function over the range z=0-6, with the IR number counts in all bands measured by SPITZER, and with the observed evolution of the mid-IR luminosity function for z=0-2.

The ROSAT-ESO Flux Limited X-ray (REFLEX) galaxy cluster survey and the 2dF Galaxy Redshift Survey (2dFGRS) respectively comprise the largest, homogeneous X-ray selected cluster catalogue and completed galaxy redshift survey. In this work we combine these two outstanding datasets in order to study the effect of the large-scale cluster environment, as traced by X-ray luminosity, on the properties of the cluster member galaxies. We measure the LX-sigma relation from the correlated dataset and find it to be consistent with recent results found in the literature. Using a sample of 19 clusters with LX>=0.36*10^44 erg s^-1 in the (0.1-2.4 keV) band, and 49 clusters with lower X-ray luminosity, we find that the fraction of early spectral type (eta<=-1.4), passively-evolving galaxies is significantly higher in the high-LX sample within R200. We extend the investigation to include composite bJ cluster luminosity functions, and find that the characteristic magnitude of the Schechter-function fit to the early-type luminosity function is fainter for the high-LX sample compared to the low-LX sample (Delta M*=0.58+/-0.14). This seems to be driven by a deficit of such galaxies with M_bJ ~ -21. In contrast, we find no significant differences between the luminosity functions of star-forming, late-type galaxies. We believe these results are consistent with a scenario in which the high-LX clusters are more dynamically evolved systems than the low-LX clusters.

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 rest-frame colour, and deducing a self-consistent treatment of k-corrections and evolution for each population. The covariance matrix for the power-spectrum 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 newly-sampled 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.)

The halo occupation distribution (HOD) describes the relation between galaxies and dark matter at the level of individual dark matter halos. The properties of galaxies residing at the centers of halos differ from those of satellite galaxies because of differences in their formation histories. Using a smoothed particle hydrodynamics (SPH) simulation and a semi-analytic (SA) galaxy formation model, we examine the separate contributions of central and satellite galaxies to the HOD, more specifically to the probability P(N|M) that a halo of virial mass M contains N galaxies of a particular class. In agreement with earlier results for dark matter subhalos, we find that the mean occupation function <N> for galaxies above a baryonic mass threshold can be approximated by a step function for central galaxies plus a power law for satellites, and that the distribution of satellite numbers is close to Poisson at fixed halo mass. For galaxy samples defined by different baryonic mass thresholds, there is a nearly linear relation between the minimum halo mass Mmin required to host a central galaxy and the mass M1 at which an average halo hosts one satellite, with M1 ~ 14 Mmin (SPH) or M1 ~ 18 Mmin (SA). The mean occupation number of young galaxies exhibits a local minimum at M ~ 10 Mmin where halos are too massive to host a young central galaxy but not massive enough to host satellites. We show that the conditional galaxy mass function at fixed halo mass cannot be described by a Schechter function because central galaxies produce a "bump" at high masses. We suggest parameterizations for the HOD and the conditional luminosity function that can be used to model observed galaxy clustering. Many of our predictions are in good agreement with recent results inferred from clustering in the Sloan Digital Sky Survey.

We investigate the metal enrichment of elliptical galaxies in the framework of hierarchical models of galaxy formation. The semi-analytical model we use in this Letter, which has been used to study the metal enrichment of the intracluster medium (ICM) by Nagashima et al., includes the effects of flows of gas and metals both into and out of galaxies and the processes of metal enrichment due to both type Ia and type II supernovae. We adopt a solar neighbourhood IMF for star formation in discs, but consider models in which starbursts have either a solar neighbourhood IMF or a top-heavy IMF. We find that the alpha-element abundance in ellipticals is consistent with observed values only if the top-heavy IMF is used. This result is consistent with our previous study on the metal enrichment of the ICM. We also discuss the abundance ratio of alpha elements to iron as a function of velocity dispersion and metallicity. We find that models with a top-heavy IMF match the alpha/Fe ratios observed in typical L_* ellipticals, but none of the models reproduce the observed increase of alpha/Fe with velocity dispersion.

The cold dark matter model has become the leading theoretical paradigm for the formation of structure in the Universe. Together with the theory of cosmic inflation, this model makes a clear prediction for the initial conditions for structure formation and predicts that structures grow hierarchically through gravitational instability. Testing this model requires that the precise measurements delivered by galaxy surveys can be compared to robust and equally precise theoretical calculations. Here we present a novel framework for the quantitative physical interpretation of such surveys. This combines the largest simulation of the growth of dark matter structure ever carried out with new techniques for following the formation and evolution of the visible components. We show that baryon-induced features in the initial conditions of the Universe are reflected in distorted form in the low-redshift galaxy distribution, an effect that can be used to constrain the nature of dark energy with next generation surveys.

The 11 known satellite galaxies within 250 kpc of the Milky Way lie close to a great circle on the sky. We use high resolution N-body simulations of galactic dark matter halos to test if this remarkable property can be understood within the context of the cold dark matter cosmology. We construct halo merger trees from the simulations and use a semianalytic model to follow the formation of satellite galaxies. We find that in all 6 of our simulations, the 11 brightest satellites are indeed distributed along thin, disk-like structures analogous to that traced by the Milky Way's satellites. This is in sharp contrast to the overall distributions of dark matter in the halo and of subhalos within it which, although triaxial, are not highly aspherical. We find that the spatial distribution of satellites is significantly different from that of the most massive subhalos but is similar to that of the subset of subhalos that had the most massive progenitors at earlier times. The elongated disk-like structure delineated by the satellites has its long axis aligned with the major axis of the dark matter halo. We interpret our results as reflecting the preferential infall of satellites along the spines of a few filaments of the cosmic web.

We use the 2dF Galaxy Redshift Survey to measure the dependence of the bJ-band galaxy luminosity function on large-scale environment, defined by density contrast in spheres of radius 8h-1Mpc, 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_bJ-5log10h <-18.5. In contrast, cluster regions have a relative excess of very bright early-type galaxies with M_bJ-5log10h < -21. Differences in the mid to faint-end population between environments are significant: at M_bJ-5log10h=-18 early and late-type 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 late-type galaxies. These observed variations, combined with our finding that the faint-end 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 sub-mm galaxies (SMGs) and Lyman-break galaxies (LBGs) in the $\Lambda$CDM cosmology. A key feature of our model is the self-consistent 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 top-heavy 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.

We investigate the metal enrichment of the intracluster medium (ICM) in the framework of hierarchical models of galaxy formation. We calculate the formation and evolution of galaxies and clusters using a semi-analytical model which includes the effects of flows of gas and metals both into and out of galaxies. For the first time in a semi-analytical model, we calculate the production of both alpha and iron-peak elements based on theoretical models for the lifetimes and ejecta of type Ia and type II supernovae (SNe Ia and SNe II). It is essential to include the long lifetimes of the SNIa progenitors in order to correctly model the evolution of the iron-peak elements. We find that if all stars form with an IMF similar to that found in the solar neighbourhood, then the metallicities of O, Mg, Si and Fe in the ICM are predicted to be 2-3 times lower than observed values. In contrast, a model (also favoured on other grounds) in which stars formed in bursts triggered by galaxy mergers have a top-heavy IMF reproduces the observed ICM abundances of O, Mg, Si and Fe. The same model predicts ratios of ICM mass to total stellar luminosity in clusters which agree well with observations. According to our model, the bulk of the metals in clusters are produced by L* and brighter galaxies. [abridged]

The 2dFGRS is used in conjunction with the 2MASSXSC to study the near-IR light and stellar mass content of the local Universe. Mock galaxy catalogues, constructed from cosmological N-body simulations and semi-analytical 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 mass-to-light ratio in the K_S band is found to increase by a factor of ~3 when going from groups with total bJ-band 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 Group-sized 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 counts-in-cells 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 non-zero `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 log-density 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 select a sample of low-redshift (z ~ 0.1) E+A galaxies from the 2dF Galaxy Redshift Survey (2dFGRS). The spectra of these objects are defined by strong hydrogen Balmer absorption lines (H-delta, H-gamma, H-beta) combined with a lack of [OII] 3727A emission, together implying a recently-truncated burst of star formation. The E+A spectrum is thus a signpost to galaxies in the process of evolution. We quantify the local environments, clustering properties and luminosity function of the E+A galaxies. We find that the environments are consistent with the ensemble of 2dFGRS galaxies: low-redshift E+A systems are located predominantly in the field, existing as isolated objects or in poor groups. However, the luminosity distribution of galaxies selected using three Balmer absorption lines H-delta-gamma-beta appears more typical of ellipticals. Indeed, morphologically these galaxies are preferentially spheroidal (E/S0) systems. In a small but significant number we find evidence for recent major mergers, such as tidal tails. We infer that major mergers are one important formation mechanism for E+A galaxies, as suggested by previous studies. At low redshift the merger probability is high in the field and low in clusters, thus these recently-formed spheroidal systems do not follow the usual morphology-density relation for ellipticals. Regarding the selection of E+A galaxies: we find that basing the Balmer-line criterion solely on H-delta absorption leads to a significant sub-population of disk systems with detectable H-alpha emission. In these objects the [OII] emission is presumably either obscured by dust or present with a low signal-to-noise ratio, whilst the (H-gamma, H-beta) absorption features are subject to emission-filling.

We measure moments of the galaxy count probability distribution function in the two-degree 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 p-point galaxy correlation functions, xibar_p, up to p=6. The hierarchical amplitudes, S_p = xibar_p/xibar_2^{p-1}, 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 non-linear 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_bJ-5logh=-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 non-linear 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 20h-1Mpc. Conversely, the reduced VPF for dark matter in a LambdaCDM universe does appear to be lognormal on small scales but deviates significantly beyond \approx 4h-1Mpc. 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 Percolation-Inferred 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 present the result of a decomposition of the 2dFGRS galaxy overdensity field into an orthonormal basis of spherical harmonics and spherical Bessel functions. Galaxies are expected to directly follow the bulk motion of the density field on large scales, so the absolute amplitude of the observed large-scale redshift-space distortions caused by this motion is expected to be independent of galaxy properties. By splitting the overdensity field into radial and angular components, we linearly model the observed distortion and obtain the cosmological constraint Omega_m^{0.6} sigma_8=0.46+/-0.06. The amplitude of the linear redshift-space distortions relative to the galaxy overdensity field is dependent on galaxy properties and, for L_* galaxies at redshift z=0, we measure beta(L_*,0)=0.58+/-0.08, and the amplitude of the overdensity fluctuations b(L_*,0) sigma_8=0.79+/-0.03, marginalising over the power spectrum shape parameters. Assuming a fixed power spectrum shape consistent with the full Fourier analysis produces very similar parameter constraints.

We use the two-degree 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 N-body 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 non-linear 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 higher-order clustering moments measured on large scales.

Complementary one-, two-, and three-dimensional tests for detecting the presence of substructure in clusters of galaxies are applied to recently obtained data from the 2dF Galaxy Redshift Survey. The sample of 25 clusters used in this study includes 16 clusters not previously investigated for substructure. Substructure is detected at or greater than the 99% CL level in at least one test for 21 of the 25 clusters studied here. From the results, it appears that low richness clusters commonly contain subclusters participating in mergers. About half of the clusters have two or more components within 0.5 h^{-1} Mpc of the cluster centroid, and at least three clusters (Abell 1139, Abell 1663, and Abell S333) exhibit velocity-position characteristics consistent with the presence of possible cluster rotation, shear, or infall dynamics. The geometry of certain features is consistent with influence by the host supercluster environments. In general, our results support the hypothesis that low richness clusters relax to structureless equilibrium states on very long dynamical time scales (if at all).

We present an analysis of the relative bias between early- and late-type galaxies in the Two-degree Field Galaxy Redshift Survey (2dFGRS). Our analysis examines the joint counts in cells between early- and late-type 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 one-point 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 power-law 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)

We develop an analytic model to calculate the rate at which galaxy disks are heated by dark matter substructures orbiting in their halos. The model takes into account the internal structure, mass function and accretion rate of satellites expected in the LambdaCDM cosmology, as well as the growth of the disk by accretion and mergers, but it ignores resonant heating of the disk and the dynamical effects of spiral arms and bars. We calibrate this model against N-body simulations and demonstrate that it is able to reproduce the N-body heating rates to within a factor of 3 in the majority of cases. Our model gives the distribution of disk scale-heights for galaxies of different luminosities. For L* spiral galaxies, it predicts a median disk thickness of only 5% of the radial scale-length if substructure is the only source of heating. The median disk thickness increases to nearly 20% of the radial scale-length when heating due to gravitational scattering of stars by molecular clouds is also included. The latter value is close to the thickness estimated observationally for the disk of the Milky Way galaxy. The distribution of disk thickness predicted by the model is also consistent with a recent observational determination for sub-L* galaxies by Bizyaev & Mitronova. Thus, the observed thickness of the stellar disks of spiral galaxies seems to be entirely compatible with the abundance of substructure in dark matter halos predicted by the standard Lambda-dominated cold dark matter model of structure formation. In an Omega_0=1 universe, our best model of galaxy formation produces similar scale-heights, a consequence of the fact that similar amounts of substructure are accreted by halos during the lifetime of the disk in Omega_0=1 and Omega_0=0.3, Lambda_0=0.7 cold dark matter cosmologies.

The construction of a catalogue of galaxy groups from the 2-degree Field Galaxy Redshift Survey (2dFGRS) is described. Groups are identified by means of a friends-of-friends percolation algorithm which has been thoroughly tested on mock versions of the 2dFGRS generated from cosmological N-body simulations. The tests suggest that the algorithm groups all galaxies that it should be grouping, with an additional 40% of interlopers. About 55% of the ~190000 galaxies considered are placed into groups containing at least two members of which ~29000 are found. Of these, ~7000 contain at least four galaxies, and these groups have a median redshift of 0.11 and a median velocity dispersion of 260km/s. This 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue represents the largest available homogeneous sample of galaxy groups. It is publicly available on the WWW.

The 2dFGRS Percolation-Inferred 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 mass-to-light 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 bj-band M/L increases by a factor of 5, whereas the rf-band 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).

We derive the fraction of blue galaxies in a sample of clusters at z < 0.11 and the general field at the same redshift. The value of the blue fraction is observed to depend on the luminosity limit adopted, cluster-centric radius and, more generally, local galaxy density, but it does not depend on cluster properties. Changes in the blue fraction are due to variations in the relative proportions of red and blue galaxies but the star formation rate for these two galaxy groups remains unchanged. Our results are most consistent with a model where the star formation rate declines rapidly and the blue galaxies tend to be dwarfs and do not favour mechanisms where the Butcher-Oemler effect is caused by processes specific to the cluster environment.

The clustering properties of local, S_{1.4 GHz} > 1 mJy, radio sources are investigated for a sample of 820 objects drawn from the joint use of the FIRST and 2dF Galaxy Redshift surveys. To this aim, we present 271 new bj < 19.45 spectroscopic counterparts of FIRST radio sources to be added to those already introduced in Magliocchetti et al. (2002). The two-point correlation function for the local radio population is found to be entirely consistent with estimates obtained for the whole sample of 2dFGRS galaxies. We estimate the parameters of the real-space correlation function xi(r)=(r/r_0)^{-\gamma}, r_0=6.7^{+0.9}_{-1.1} Mpc and \gamma=1.6\pm 0.1, where h=0.7 is assumed. Different results are instead obtained if we only consider sources that present signatures of AGN activity in their spectra. These objects are shown to be very strongly correlated, with r_0=10.9^{+1.0}_{-1.2} Mpc and \gamma=2\pm 0.1, a steeper slope than has been claimed in other recent works. No difference is found in the clustering properties of radio-AGNs of different radio luminosity. These results show that AGN-fuelled sources reside in dark matter halos more massive than \sim 10^{13.4} M_{\sun}},higher the corresponding figure for radio-quiet QSOs. This value can be converted into a minimum black hole mass associated with radio-loud, AGN-fuelled objects of M_{BH}^{min}\sim 10^9 M_{\sun}. The above results then suggest -at least for relatively faint radio objects -the existence of a threshold black hole mass associated with the onset of significant radio activity such as that of radio-loud AGNs; however, once the activity is triggered, there appears to be no evidence for a connection between black hole mass and level of radio output. (abridged)

We analyse the observed correlation between galaxy environment and H-alpha emission line strength, using volume-limited samples and group catalogues of 24968 galaxies drawn from the 2dF Galaxy Redshift Survey (Mb<-19.5) and the Sloan Digital Sky Survey (Mr<-20.6). We characterise the environment by 1) Sigma_5, the surface number density of galaxies determined by the projected distance to the 5th nearest neighbour; and 2) rho1.1 and rho5.5, three-dimensional density estimates obtained by convolving the galaxy distribution with Gaussian kernels of dispersion 1.1 Mpc and 5.5 Mpc, respectively. We find that star-forming and quiescent galaxies form two distinct populations, as characterised by their H-alpha equivalent width, EW(Ha). The relative numbers of star-forming and quiescent galaxies varies strongly and continuously with local density. However, the distribution of EW(Ha) amongst the star-forming population is independent of environment. The fraction of star-forming galaxies shows strong sensitivity to the density on large scales, rho5.5, which is likely independent of the trend with local density, rho1.1. We use two differently-selected group catalogues to demonstrate that the correlation with galaxy density is approximately independent of group velocity dispersion, for sigma=200-1000 km/s. Even in the lowest density environments, no more than ~70 per cent of galaxies show significant H-alpha emission. Based on these results, we conclude that the present-day correlation between star formation rate and environment is a result of short-timescale mechanisms that take place preferentially at high redshift, such as starbursts induced by galaxy-galaxy interactions.

We reconstruct the underlying density field of the 2 degree Field Galaxy Redshift Survey (2dFGRS) for the redshift range 0.035<z<0.200 using the Wiener Filtering method. The Wiener Filter suppresses shot noise and accounts for selection and incompleteness effects. The method relies on prior knowledge of the 2dF power spectrum of fluctuations and the combination of matter density and bias parameters however the results are only slightly affected by changes to these parameters. We present maps of the density field in two different resolutions: 5 Mpc/h and 10 Mpc/h. We identify all major superclusters and voids in the survey. In particular, we find two large superclusters and two large local voids. A version of this paper with full set of colour maps can be found at http://www.ast.cam.ac.uk/~pirin.

We investigate the physical mechanisms that shape the luminosity function. Beginning with the mass function of dark matter halos, we show how gas cooling, photoionization, feedback, merging and thermal conduction affect the shape of the luminosity function. We consider three processes whereby supernovae can affect the forming galaxy: (1) reheating of disk gas to the halo temperature; (2) expansion of the diffuse halo gas; (3) expulsion of cold disk gas from the halo. While feedback of form (1) is able to flatten the faint end of the luminosity function, this alone does not produce the sharp cut-off observed at large luminosities. Feedback of form (2) is also unable to solve this problem. The relative paucity of very bright galaxies can only be explained if cooling in massive halos is strongly suppressed. Conduction is a promising mechanism, but an uncomfortably high efficiency is required to suppress cooling to the desired level. If, instead, superwinds are responsible for the lack of bright galaxies, then the total energy budget required to obtain a good match to the galaxy luminosity function greatly exceeds the energy available from supernova explosions. The mechanism is only viable if the formation of central supermassive black holes play a crucial role in limiting the amount of stars that form in the host galaxy. (abridged)

The Millennium Galaxy Catalogue (MGC) is a deep ($\mu_{\rm B,lim}=26$ mag arcsec$^{-2}$), wide field CCD imaging survey, covering 37.5deg$^2$ and is completely contained within the 2dFGRS and SDSS-EDR. We compare the photometry and completeness of the 2dFGRS and the SDSS-EDR with the MGC. We have also undertaken a photometric comparison to SCOS and SDSS-DR1 data. We find that $B_{MGC}-B_{2dF}=0.035$ mag with an uncertainty of 0.142 mag per galaxy, $B_{MGC}-B_{SCOS}=0.032$ mag with an uncertainty of 0.108 mag, $B_{MGC}-B_{SDSS-EDR}=0.032$ mag with an uncertainty of 0.094 mag, and $B_{MGC}-B_{SDSS-DR1}=0.039$ mag with an uncertainty of 0.086 mag. We find that high surface brightness 2dFGRS galaxies are systematically too faint. In the SDSS there is a weak non-linear scale error, which is negligible for faint galaxies. LSBGs in the SDSS are systematically fainter. We find that the 2dFGRS catalogue has 5.2% stellar contamination, 7.0% of objects are resolved into 2 or more by the MGC and is 8.7% incomplete compared to the MGC. From our all object spectroscopic survey we find that the MGC is itself misclassifying 5.6% of galaxies as stars, hence the 2dFGRS misses 14.3% of the galaxies. The SDSS-EDR galaxy catalogue has 1.3% stellar contamination and 5.3% of galaxies misclassified as stars, and is 1.8% incomplete compared to the MGC. Altogether 7.1% of the total galaxy population are missing from the SDSS-EDR catalogue from incompleteness or misclassification.

We present a detailed analysis of the two-point correlation function, from the 2dF Galaxy Redshift Survey (2dFGRS). We estimate the redshift-space correlation function, xi(s), from which we measure the redshift-space clustering length, s_0=6.82+/-0.28 Mpc/h. We also estimate the projected correlation function, Xi(sigma), and the real-space correlation function, xi(r), which can be fit by a power-law, with r_0=5.05+/-0.26Mpc/h, gamma_r=1.67+/-0.03. For r>20Mpc/h, xi drops below a power-law as is expected in the popular LCDM model. The ratio of amplitudes of the real and redshift-space correlation functions on scales of 8-30Mpc/h gives an estimate of the redshift-space distortion parameter beta. The quadrupole moment of xi on scales 30-40Mpc/h provides another estimate of beta. We also estimate the distribution function of pairwise peculiar velocities, f(v), including rigorously the effect of infall velocities, and find that it is well fit by an exponential. The accuracy of our xi measurement is sufficient to constrain a model, which simultaneously fits the shape and amplitude of xi(r) and the two redshift-space distortion effects parameterized by beta and velocity dispersion, a. We find beta=0.49+/-0.09 and a=506+/-52km/s, though the best fit values are strongly correlated. We measure the variation of the peculiar velocity dispersion with projected separation, a(sigma), and find that the shape is consistent with models and simulations. Using the constraints on bias from recent estimates, and taking account of redshift evolution, we conclude that beta(L=L*,z=0)=0.47+/-0.08, and that the present day matter density of the Universe is 0.3, consistent with other 2dFGRS estimates and independent analyses.

The 2dF Galaxy Redshift Survey (2dFGRS) has obtained spectra for 245591 sources, mainly galaxies, brighter than a nominal extinction-corrected magnitude limit of b_J=19.45. Reliable redshifts were measured for 221414 galaxies. The galaxies are selected from the extended APM Galaxy Survey and cover an area of approximately 1500 square degrees in three regions: an NGP strip, an SGP strip and random fields scattered around the SGP strip. This paper describes the 2dFGRS final data release of 30 June 2003 and complements Colless et al. (2001), which described the survey and the initial 100k data release. The 2dFGRS database and full documentation are available on the WWW at http://www.mso.anu.edu.au/2dFGRS/

The halo occupation distribution (HOD) describes the bias between galaxies and dark matter by specifying (a) the probability P(N|M) that a halo of virial mass M contains N galaxies of a particular class and (b) the relative distributions of galaxies and dark matter within halos. We calculate predicted HODs for a Lambda-CDM cosmological model using an SPH hydrodynamic simulation and a semi-analytic (SA) galaxy formation model. Although the two methods predict different galaxy mass functions, their HOD predictions agree remarkably well. For mass-selected samples, the mean occupation <N(M)> exhibits a sharp cutoff at low halo masses, a slowly rising plateau for <N>~1-2, and a more steeply rising high occupancy regime. At low <N>, the mean pair and triple counts are well below Poisson expectations, with important consequences for small scale behavior of 2- and 3-point correlation functions. The HOD depends strongly on galaxy age, with high mass halos populated mainly by old galaxies and low mass halos by young galaxies. The SPH simulation supports several simplifying assumptions about HOD bias: the most massive galaxy in a halo usually lies close to the center and moves near the halo's mean velocity; satellite galaxies have the same radial profile and velocity dispersion as the dark matter; and the mean occupation at fixed halo mass is independent of the halo's larger scale environment. By applying the SPH and SA HODs to a large volume N-body simulation, we show that both methods predict slight, observable departures from a power-law galaxy correlation function. The predicted HODs are closely tied to the underlying galaxy formation physics, they offer useful guidance to theoretical models of galaxy clustering, and they will be tested empirically by ongoing analyses of galaxy redshift surveys. (Shortened)

We compare the properties of local spiral galaxies with the predictions of the Cole et al. semi-analytic model of hierarchical galaxy formation, in order to gain insight into the baryonic processes that were responsible for shaping these galaxies. On the whole, the models reproduce the properties of present-day spirals (such as SFRs, SFHs, scale sizes, metallicities, and gas fractions) rather well. In particular, we find that once the effects of dust and variations in stellar populations have been taken into account, published spiral galaxy scale-size distributions derived from optical data (with logarithmic widths ~0.3) can be reconciled with the width of the stellar mass scale-size distribution predicted by the semi-analytic model (~0.5). There are some illuminating discrepancies between the observations and the model predictions. The model colours of luminous spiral galaxies are somewhat too blue and those of faint galaxies somewhat too red, most likely indicating shortcomings in the way that gas is accreted by spiral galaxies. Furthermore, the model produces too few luminous spiral galaxies. These difficulties could be resolved by altering the way in which gas cooling is treated or, perhaps, by adopting a higher baryon fraction and invoking galactic `superwinds.' Secondly, stellar M/Ls are found to be as high as observations permit. Yet, typically 60% of the mass in the inner half-light radius of the model galaxies is dark. This results in an offset between the model and observed spiral galaxy luminosity-linewidth relation. This could be resolved by substantially reducing the mass of baryons which make it into a galaxy disc (with an attendant decrease in stellar M/L), or by modifying the assumed dark matter profile to include less dark matter in the inner parts. [Abridged]

We have calculated the two-point correlation functions in redshift space, xi(sigma,pi), for galaxies of different spectral types in the 2dF Galaxy Redshift Survey. Using these correlation functions we are able to estimate values of the linear redshift-space distortion parameter, beta = Omega_m^0.6/b, the pairwise velocity dispersion, a, and the real-space correlation function, xi(r), for galaxies with both relatively low star-formation rates (for which the present rate of star formation is less than 10% of its past averaged value) and galaxies with higher current star-formation activity. At small separations, the real-space clustering of passive galaxies is very much stronger than that of the more actively star-forming galaxies; the correlation-function slopes are respectively 1.93 and 1.50, and the relative bias between the two classes is a declining function of radius. On scales larger than 10 h^-1 Mpc there is evidence that the relative bias tends to a constant, b(passive)/b(active) ~ 1. This result is consistent with the similar degrees of redshift-space distortions seen in the correlation functions of the two classes -- the contours of xi(sigma,pi) require beta(active)=0.49+/-0.13, and beta(passive)=0.48+/-0.14. The pairwise velocity dispersion is highly correlated with beta. However, despite this a significant difference is seen between the two classes. Over the range 8-20 h^-1 Mpc, the pairwise velocity dispersion has mean values 416+/-76 km/s and 612+/-92 km/s for the active and passive galaxy samples respectively. This is consistent with the expectation from morphological segregation, in which passively evolving galaxies preferentially inhabit the cores of high-mass virialised regions.

We have determined the composite luminosity function (LF) for galaxies in 60 clusters from the 2dF Galaxy Redshift Survey. The LF spans the range $-22.5<M_{b_{\rm J}}<-15$, and is well-fitted by a Schechter function with ${M_{b_{\rm J}}}^{*}=-20.07\pm0.07$ and $\alpha=-1.28\pm0.03$ ($H_0$=100 km s$^{-1}$ Mpc$^{-1}$, $\Omega_M$=0.3, $\Omega_\Lambda$=0.7). It differs significantly from the field LF of \cite{mad02}, having a characteristic magnitude that is approximately 0.3 mag brighter and a faint-end slope that is approximately 0.1 steeper. There is no evidence for variations in the LF across a wide range of cluster properties. However the LF of early-type galaxies in clusters is both brighter and steeper than its field counterpart. The differences between the field and cluster LFs for the various spectral types can be qualitatively explained by the suppression of star formation in the dense cluster environment, together with mergers to produce the brightest early-type galaxies.

Using semi-analytic modeling techniques, we calculate the luminosity function of galaxy populations residing in cold dark matter halos of different mass. We pay particular attention to the influence of the reionization of the Universe on the number of faint galaxies and to the effects of dynamical friction and tidal limitation of satellites on the number of bright galaxies. We find substantial differences in the shapes of the galaxy luminosity functions in halos of different mass which reflect generic features of the cold dark matter model of galaxy formation and thus offer the opportunity to test it. We then consider how the individual halo luminosity functions combine together to produce the global luminosity function. Surprisingly, the global function ends up having a shallower faint end slope than those of the constituent halo luminosity functions. We compare our model predictions with the limited datasets compiled by Trentham & Hodgkin. We find good agreement with the luminosity functions measured in the Virgo and Coma clusters but significant disagreement with the luminosity functions measured in the Local Group and Ursa Minor cluster. We speculate on possible inadequacies in our modeling and in the existing observational samples. The luminosity functions of galaxies in groups and clusters identified in the 2dF and SDSS galaxy redshift surveys offer the prospect of testing galaxy formation models in detail.

We compare the results of two techniques used to calculate the evolution of cooling gas during galaxy formation: Smooth Particle Hydrodynamics (SPH) simulations and semi-analytic modelling. We improve upon the earlier statistical comparison of Benson et al. by taking halo merger histories from the dark matter component of the SPH simulation, which allows us to compare the evolution of galaxies on an object-by-object basis in the two treatments. We use a ``stripped-down'' version of the semi-analytic model described by Helly et al. which includes only shock heating and radiative cooling of gas and which is adjusted to mimic the resolution and other parameters of a comparison SPH simulation as closely as possible. We compare the total mass of gas that cools in halos of different mass as a function of redshift as well as the masses and spatial distribution of individual ``galaxies.'' At redshift z=0, the cooled gas mass in well-resolved halos agrees remarkably well (to better than ~20%) in the SPH simulation and stripped-down semi-analytic model. At high redshift, resolution effects in the simulation become increasingly important and, as a result, more gas tends to cool in low mass halos in the SPH simulation than in the semi-analytic model. The cold gas mass function of individual galaxies in the two treatments at z=0 also agrees very well and, when the effects of mergers are accounted for, the masses of individual galaxies and their 2-point correlation functions are also in excellent agreement in the two treatments. Thus, our comparison confirms and extends the earlier conclusion of Benson et al. that SPH simulations and semi-analytic models give consistent results for the evolution of cooling galactic gas.

We develop a hybrid galaxy formation model which uses outputs from an N-body simulation to follow the merger histories (or ``merger trees'') of dark matter halos and treats baryonic processes, such as the cooling of gas within halos and subsequent star formation, using the semi-analytic model of Cole et al. We compare this hybrid model to an otherwise identical model which utilises merger tree realisations generated by a Monte-Carlo algorithm and find that, apart from the limited mass resolution imposed by the N-body particle mass, the only significant differences between the models are due to the known discrepancy between the distribution of halo progenitor masses predicted by extended Press Schechter theory and that found in N-body simulations. We investigate the effect of limited mass resolution on the hybrid model by comparing to a purely semi-analytic model with greatly improved mass resolution. We find that the mass resolution of the simulation we use, which has a particle mass of 1.4x10^10 Msolar/h, is insufficient to produce a reasonable luminosity function for galaxies with magnitudes in the b_J band fainter than -17.

We constrain flat cosmological models with a joint likelihood analysis of a new compilation of data from the cosmic microwave background (CMB) and from the 2dF Galaxy Redshift Survey (2dFGRS). Fitting the CMB alone yields a known degeneracy between the Hubble constant h and the matter density Omega_m, which arises mainly from preserving the location of the peaks in the angular power spectrum. This `horizon-angle degeneracy' is considered in some detail and shown to follow a simple relation Omega_m h^{3.4} = constant. Adding the 2dFGRS power spectrum constrains Omega_m h and breaks the degeneracy. If tensor anisotropies are assumed to be negligible, we obtain values for the Hubble constant h=0.665 +/- 0.047, the matter density Omega_m=0.313 +/- 0.055, and the physical CDM and baryon densities Omega_c h^2 = 0.115 +/- 0.009, Omega_b h^2 = 0.022 +/- 0.002 (standard rms errors). Including a possible tensor component causes very little change to these figures; we set a upper limit to the tensor-to-scalar ratio of r<0.7 at 95% confidence. We then show how these data can be used to constrain the equation of state of the vacuum, and find w<-0.52 at 95% confidence. The preferred cosmological model is thus very well specified, and we discuss the precision with which future CMB data can be predicted, given the model assumptions. The 2dFGRS power-spectrum data and covariance matrix, and the CMB data compilation used here, are available from http://www.roe.ac.uk/~wjp/

We constrain f_nu = Omega_nu / Omega_m, the fractional contribution of neutrinos to the total mass density in the Universe, by comparing the power spectrum of fluctuations derived from the 2dF Galaxy Redshift Survey with power spectra for models with four components: baryons, cold dark matter, massive neutrinos and a cosmological constant. Adding constraints from independent cosmological probes we find f_nu < 0.13 (at 95% confidence) for a prior of 0.1< Omega_m <0.5, and assuming the scalar spectral index n=1. This translates to an upper limit on the total neutrino mass and m_nu,tot < 1.8 eV for "concordance" values of Omega_m and the Hubble constant. Very similar results are obtained with a prior on Omega_m from Type Ia supernovae surveys, and with marginalization over n.

The 2dF Galaxy Redshift Survey is the first to observe more than 100,000 redshifts. This allows precise measurements of many of the key statistics of galaxy clustering, in particular redshift-space distortions and the large-scale power spectrum. This paper presents the current 2dFGRS results in these areas. Redshift-space distortions are detected with a high degree of significance, confirming the detailed Kaiser distortion from large-scale infall velocities, and measuring the distortion parameter beta equiv Omega_m^{0.6}/b = 0.43 +- 0.07. The power spectrum is measured to < 10% accuracy for k > 0.02 h Mpc^{-1}, and is well fitted by a CDM model with Omega_m h = 0.20 +- 0.03 and a baryon fraction of 0.15 +- 0.07. A joint analysis with CMB data requires Omega_m = 0.29 +- 0.05, assuming scalar fluctuations, but no priors on other parameters. Two methods are used to determine the large-scale bias parameter: an internal bispectrum analysis yields b = 1.04 +- 0.11, in very good agreement with the b = 1.10 +- 0.08 obtained from a joint 2dFGRS+CMB analysis, again assuming scalar fluctuations. These figures refer to galaxies of approximate luminosity 2L^*; luminosity dependence of clustering is detected at high significance, and is well described by b/b^* = 0.85 + 0.15(L/L^*).

We compare the amplitudes of fluctuations probed by the 2dF Galaxy Redshift Survey and by the latest measurements of the Cosmic Microwave Background anisotropies. By combining the 2dFGRS and CMB data we find the linear-theory rms mass fluctuations in 8 Mpc/h spheres to be sigma_8 = 0.73 +-0.05 (after marginalization over the matter density parameter Omega_m and three other free parameters). This normalization is lower than the COBE normalization and previous estimates from cluster abundance, but it is in agreement with some revised cluster abundance determinations. We also estimate the scale-independent bias parameter of present-epoch L_s = 1.9L_* APM-selected galaxies to be b(L_s,z=0) = 1.10 +- 0.08 on comoving scales of 0.02 < k < 0.15 h/Mpc. If luminosity segregation operates on these scales, L_* galaxies would be almost un-biased, b(L_*,z=0) = 0.96. These results are derived by assuming a flat Lambda-CDM Universe, and by marginalizing over other free parameters and fixing the spectral index n=1 and the optical depth due to reionization tau=0. We also study the best fit pair (Omega_m,b), and the robustness of the results to varying n and tau. Various modelling corrections can each change the resulting b by 5-15 per cent. The results are compared with other independent measurements from the 2dFGRS itself, and from the SDSS, cluster abundance and cosmic shear.

We compute the bispectrum of the 2dF Galaxy Redshift Survey (2dFGRS) and use it to measure the bias parameter of the galaxies. This parameter quantifies the strength of clustering of the galaxies relative to the mass in the Universe. By analysing 80 million triangle configurations in the wavenumber range 0.1 < k < 0.5 h/Mpc (i.e. on scales roughly between 5 and 30 Mpc/h) we find that the linear bias parameter is consistent with unity: b_1=1.04 pm 0.11, and the quadratic (nonlinear) bias is consistent with zero: b_2=-0.054 pm 0.08. Thus, at least on large scales, optically-selected galaxies do indeed trace the underlying mass distribution. The bias parameter can be combined with the 2dFGRS measurement of the redshift distortion parameter beta = Omega_m^{0.6}/b_1, to yield Omega_m = 0.27 pm 0.06 for the matter density of the Universe, a result which is determined entirely from this survey, independently of other datasets. Our measurement of the matter density of the Universe should be interpreted as Omega_m at the effective redshift of the survey (z=0.17).

We use more than 110500 galaxies from the 2dF galaxy redshift survey (2dFGRS) to estimate the b_J-band galaxy luminosity function at redshift z=0, taking account of evolution, the distribution of magnitude measurement errors and small corrections for incompletenessin the galaxy catalogue. Throughout the interval -16.5>M- 5log h>-22, the luminosity function is accurately described by a Schechter function with M* -5log h =-19.66+/-0.07, alpha=-1.21+/-0.03 and phistar=(1.61+/-0.08) 10^{-2} h^3/Mpc^3, giving an integrated luminosity density of rho_L=(1.82+/-0.17) 10^8 h L_sol/Mpc^3 (assuming an Omega_0=0.3, Lambda_0=0.7 cosmology). The quoted errors have contributions from the accuracy of the photometric zeropoint, large scale structure in the galaxy distribution and, importantly, from the uncertainty in the appropriate evolutionary corrections. Our luminosity function is in excellent agreement with, but has much smaller statistical errors than an estimate from the Sloan Digital Sky Survey (SDSS) data when the SDSS data are accurately translated to the b_J-band and the luminosity functions are normalized in the same way. We use the luminosity function, along with maps describing the redshift completeness of the current 2dFGRS catalogue, and its weak dependence on apparent magnitude, to define a complete description of the 2dFGRS selection function. Details and tests of the calibration of the 2dFGRS photometric parent catalogue are also presented.

We have measured the equivalent width of the H-alpha emission line for 11006 galaxies brighter than M_b=-19 (LCDM) at 0.05<z<0.1 in the 2dF Galaxy Redshift Survey (2dF), in the fields of seventeen known galaxy clusters. The limited redshift range ensures that our results are insensitive to aperture bias, and to residuals from night sky emission lines. We use these measurements to trace mustar, the star formation rate normalized to Lstar, as a function of distance from the cluster centre, and local projected galaxy density. We find that the distribution of mustar steadily skews toward larger values with increasing distance from the cluster centre, converging to the field distribution at distances greater than ~3 times the virial radius. A correlation between star formation rate and local projected density is also found, which is independent of cluster velocity dispersion and disappears at projected densities below ~1 galaxy (brighter than M_b=-19) per Mpc^2. This characteristic scale corresponds approximately to the mean density at the cluster virial radius. The same correlation holds for galaxies more than two virial radii from the cluster centre. We conclude that environmental influences on galaxy properties are not restricted to cluster cores, but are effective in all groups where the density exceeds this critical value. The present day abundance of such systems, and the strong evolution of this abundance, makes it likely that hierarchical growth of structure plays a significant role in decreasing the global average star formation rate. Finally, the low star formation rates well beyond the virialised cluster rule out severe physical processes, such as ram pressure stripping of disk gas, as being completely responsible for the variations in galaxy properties with environment.

Advances in extragalactic astronomy have prompted the development of increasingly realistic models which aim to describe the formation and evolution of galaxies. We review the philosophy behind one such technique, called semi-analytic modelling, and explain the relation between this approach and direct simulations of gas dynamics. Finally, we present model predictions for the evolution of the stellar mass of galaxies in a universe in which structure formation is hierarchical.

We present the bivariate brightness distribution (BBD) for the 2dF Galaxy Redshift Survey (2dFGRS) based on a preliminary subsample of 45,000 galaxies. The BBD is an extension of the galaxy luminosity function incorporating surface brightness information. It allows the measurement of the local luminosity density, j_B, and the galaxy luminosity and surface brightness distributions while accounting for surface brightness selection biases. The recovered 2dFGRS BBD shows a strong surface brightness-luminosity relation (M_B~2.4\mu_e). The luminosity-density is dominated by normal galaxies and the luminosity-density peak lies away from the selection boundaries implying that the 2dFGRS is complete and that luminous low surface brightness galaxies are rare. The final value we derive for the local luminosity-density, inclusive of surface brightness corrections, is: j_B=2.49+/-0.20x10^8 h L_solar Mpc^-3. Representative Schechter function parameters are: M*=-19.75+/-0.05, phi*=2.02+/-0.02x10^-2 and alpha=-1.09+/-0.03. Extending the conventional methodology to incorporate surface brightness selection effects has resulted in an increase in the luminosity-density of 37%.

We calculate the optical bj luminosity function of the 2dF Galaxy Redshift Survey (2dFGRS) for different subsets defined by their spectral properties. These spectrally selected subsets are defined using a new parameter, eta, which is a linear combination of the first two projections derived from a Principal Component Analysis. This parameter eta identifies the average emission and absorption line strength in the galaxy rest-frame spectrum and hence is a useful indicator of the present star formation. We use a total of 75,000 galaxies in our calculations, chosen from a sample of high signal-to-noise ratio, low redshift galaxies observed before January 2001. We find that there is a systematic steepening of the faint end slope (alpha) as one moves from passive (alpha = -0.54) to active (alpha = -1.50) star-forming galaxies, and that there is also a corresponding faintening of the rest-frame characteristic magnitude M* - 5\log_10(h) (from -19.6 to -19.2). We also show that the Schechter function provides a poor fit to the quiescent (Type 1) LF for very faint galaxies (M - 5log_10(h) fainter than -16.0), perhaps suggesting the presence of a significant dwarf population. The luminosity functions presented here give a precise confirmation of the trends seen previously in a much smaller preliminary 2dFGRS sample, and in other surveys. We also present a new procedure for determining self-consistent K-corrections and investigate possible fibre-aperture biases.

We use redshift determinations and spectral analysis of galaxies in the 2dF Galaxy Redshift Survey to study the properties of local radio sources with $S\ge 1$ mJy. 557 objects drawn from the FIRST survey, corresponding to 2.3 per cent of the total radio sample, are found in the 2dFGRS cataloge within the area $9^h 48^m \simlt {\rm RA}({\rm 2000}) \simlt 14^h 32^m$ and $-2.77^\circ \simlt {\rm dec}({\rm 2000}) \simlt 2.25^\circ$, down to a magnitude limit $b_J= 19.45$. The excellent quality of 2dF spectra allows us to divide these sources into classes, according to their optical spectra. Radio galaxies (i.e. absorption systems, mainly belonging to the class of FR I sources, whose spectra may or may not show emission lines due to AGN activity) make up 63 per cent of the spectroscopic sample; they are characterized by relatively high radio-to-optical ratios, red colours and luminosities $10^{21}\simlt P_{1.4{\rm GHz}} W Hz^${-1}$ sr$^{-1}$ \simlt 10^{24}$. Actively star-forming galaxies instead contribute with 30 per cent to the total sample. These objects -- principally found at low redshifts ($z \simlt 0.1$) -- in general show low radio-to-optical ratios, blue colours and low radio luminosities. We also found 18 Seyfert 2 galaxies and 4 Seyfert 1's. Analysis of the local radio luminosity function shows that radio galaxies are well described by models that assume pure luminosity evolution, at least down to radio powers $P_{1.4{\rm GHz}}\simlt 10^{20.5} W Hz$^{-1}$ sr$^{-1}$ Late-type galaxies, whose relative contribution to the RLF is found to be lower than was predicted by previous work, present a luminosity function which is compatible with the IRAS galaxy LF, and therefore plausibly constitute the radio counterpart of the dusty spirals and starbursts that dominate the counts at 60 $\mu$m.

We have carried out a study of known clusters within the 2dF Galaxy Redshift Survey (2dFGRS) observed areas and have identified 431 Abell, 173 APM and 343 EDCC clusters. Precise redshifts, velocity dispersions and new centroids have been measured for the majority of these objects, and this information has been used to study the completeness of these catalogues, the level of contamination from foreground and background structures along the cluster's line of sight, the space density of the clusters as a function of redshift, and their velocity dispersion distributions. We find that the Abell and EDCC catalogues are contaminated at the level of about 10%, whereas the APM catalogue suffers only 5% contamination. If we use the original catalog centroids, the level of contamination rises to approximately 15% for the Abell and EDCC catalogues, showing that the presence of foreground and background groups may alter the richness of clusters in these catalogues. There is a deficiency of clusters at $z \sim 0.05$ that may correspond to a large underdensity in the Southern hemisphere. From the cumulative distribution of velocity dispersions for these clusters, we derive an upper limit to the space density of $\sigma > 1000 \kms$ clusters of $3.6 \times 10^{-6} \hdens$. This result is used to constrain models for structure formation; our data favour low-density cosmologies, subject to the usual assumptions concerning the shape and normalization of the power spectrum.

We have cross-matched the 1.4 GHz NRAO VLA Sky Survey (NVSS) with the first 210 fields observed in the 2dF Galaxy Redshift Survey (2dFGRS), covering an effective area of 325 square degrees (about 20% of the final 2dFGRS area). This yields a set of optical spectra of 912 candidate NVSS counterparts, of which we identify 757 as genuine radio IDs - the largest and most homogeneous set of radio-source spectra ever obtained. The 2dFGRS radio sources span the redshift range z=0.005 to 0.438, and are a mixture of active galaxies (60%) and star-forming galaxies (40%). About 25% of the 2dFGRS radio sources are spatially resolved by NVSS, and the sample includes three giant radio galaxies with projected linear size greater than 1 Mpc. The high quality of the 2dF spectra means we can usually distinguish unambiguously between AGN and star-forming galaxies. We have made a new determination of the local radio luminosity function at 1.4 GHz for both active and star-forming galaxies, and derive a local star-formation density of 0.022+/-0.004 solar masses per year per cubic Mpc. (Ho=50 km/s/Mpc).

We investigate the dependence of galaxy clustering on luminosity and spectral type using the 2dF Galaxy Redshift Survey (2dFGRS). Spectral types are assigned using the principal component analysis of Madgwick et al. We divide the sample into two broad spectral classes: galaxies with strong emission lines (`late-types'), and more quiescent galaxies (`early-types'). We measure the clustering in real space, free from any distortion of the clustering pattern due to peculiar velocities, for a series of volume-limited samples. The projected correlation functions of both spectral types are well described by a power law for transverse separations in the range 2 < (sigma/Mpc/h) < 15, with a marginally steeper slope for early-types than late-types. Both early and late types have approximately the same dependence of clustering strength on luminosity, with the clustering amplitude increasing by a factor of ~2.5 between L* and 4 L*. At all luminosities, however, the correlation function amplitude for the early-types is ~50% higher than that of the late-types. These results support the view that luminosity, and not type, is the dominant factor in determining how the clustering strength of the whole galaxy population varies with luminosity.

We present the first results on the history of star formation in the Universe based on the `cosmic spectrum', in particular, the volume-averaged, luminosity-weighted, stellar absorption line spectrum of present day galaxies from the 2dFGRS. This method is novel in that unlike previous studies it is not an estimator based on total luminosity density. The cosmic spectrum is fitted with models of population synthesis, tracing the history of star formation prior to the epoch of the observed galaxies, using a method we have developed which decouples continuum and spectral-line variations and is robust against spectrophotometric uncertainties. The cosmic spectrum can only be fitted with models incorporating chemical evolution and indicates there was a peak of star-formation rate in the past of at least three times the current value and that the increase back to z=1, assuming it scales as (1+z)^beta, has a strong upper limit of beta<5. We find in the general case there is some model degeneracy between star formation at low and high redshift. However, if we incorporate previous work on star formation at z<1 we can put strong upper limits on the SFR at z>1: e.g., if beta>2 then the SFR for 1<z<5 scales as (1+z)^alpha with alpha<2. This is equivalent to stating that no more than 80% of stars in the Universe formed at z>1. Our results are consistent with the best-fit results from compilations of cosmic SFR estimates based on UV luminosity density, which give 1.8<beta<2.9 and -1.0<alpha<0.7, and are also consistent with estimates of Omega_stars based on the K-band luminosity density.

We combine the 2MASS extended source catalogue and the 2dFGRS to produce an IR selected galaxy catalogue with 17,173 measured redshifts. We use this extensive dataset to estimate the J and K-band galaxy luminosity functions. The LFs are fairly well fit by Schechter functions with J: M*-5log h= -22.36+/-0.02, alpha= -0.93+/-0.04, Phi=0.0104+/-0.0016 h^3/Mpc^3 and K: M*-5log h= -23.44+/-0.03, alpha=-0.96+/-0.05, Phi=0.0108+/-0.0016 h^3/Mpc^3 (2MASS Kron magnitudes). These parameters assume a cosmological model with Omega=0.3 and Lambda=0.7. With datasets of this size, systematic rather than random errors are the dominant source of uncertainty in the determination of the LF. We carry out a careful investigation of possible systematic effects in our data. The surface brightness distribution of the sample shows no evidence that significant numbers of low surface brightness or compact galaxies are missed by the survey. We estimate the present-day distributions of B-K and J-K colours as a function of absolute magnitude and use models of the galaxy stellar populations, constrained by the observed optical and infrared colours, to infer the galaxy stellar mass function. Integrated over all galaxy masses, this yields a total mass fraction in stars (in units of the critical mass density) of Omega_*.h= (1.6+/-0.24)/10^3 for a Kennicutt IMF and Omega_*.h= (2.9+/-0.43)/10^3 for a Salpeter IMF. These values agree with those inferred from observational estimates of the star formation history of the universe provided that dust extinction corrections are modest.

The 2dF Galaxy Redshift Survey (2dFGRS) is designed to measure redshifts for approximately 250000 galaxies. This paper describes the survey design, the spectroscopic observations, the redshift measurements and the survey database. The 2dFGRS uses the 2dF multi-fibre spectrograph on the Anglo-Australian Telescope, which is capable of observing 400 objects simultaneously over a 2-degree diameter field. The source catalogue for the survey is a revised and extended version of the APM galaxy catalogue, and the targets are galaxies with extinction-corrected magnitudes brighter than b_J=19.45. The main survey regions are two declination strips, one in the southern Galactic hemisphere spanning 80deg x 15deg around the SGP, and the other in the northern Galactic hemisphere spanning 75deg x 10deg along the celestial equator; in addition, there are 99 fields spread over the southern Galactic cap. The survey covers 2000 sq.deg and has a median depth of z=0.11. Adaptive tiling is used to give a highly uniform sampling rate of 93% over the whole survey region. Redshifts are measured from spectra covering 3600A-8000A at a two-pixel resolution of 9.0A and a median S/N of 13 per pixel. All redshift identifications are visually checked and assigned a quality parameter Q in the range 1-5; Q>=3 redshifts are 98.4% reliable and have an rms uncertainty of 85 km/s. The overall redshift completeness for Q>=3 redshifts is 91.8%, but this varies with magnitude from 99% for the brightest galaxies to 90% for objects at the survey limit. The 2dFGRS database is available on the WWW at http://www.mso.anu.edu.au/2dFGRS

We investigate the dependence of the strength of galaxy clustering on intrinsic luminosity using the Anglo-Australian two degree field galaxy redshift survey (2dFGRS). The 2dFGRS is over an order of magnitude larger than previous redshift surveys used to address this issue. We measure the projected two-point correlation function of galaxies in a series of volume-limited samples. The projected correlation function is free from any distortion of the clustering pattern induced by peculiar motions and is well described by a power-law in pair separation over the range 0.1 < r /h Mpc < 10. The clustering of L* galaxies in real space is well fit by a correlation length r0 = 4.9 +/- 0.3 /h Mpc and power-law slope gamma = 1.71 +/- 0.06. The clustering amplitude increases slowly with absolute magnitude for galaxies fainter than M*, but rises more strongly at higher luminosities. At low luminosities, our results agree with measurements from the SSRS2 by Benoist et al. However, we find a weaker dependence of clustering strength on luminosity at the highest luminosities. The correlation function amplitude increases by a factor of 4.0 between $M_{b_{J}} -5\log_{10}h = -18$ and -22.5, and the most luminous galaxies are 3.0 times more strongly clustered than L* galaxies. The power-law slope of the correlation function shows remarkably little variation for samples spanning a factor of 20 in luminosity. Our measurements are in very good agreement with the predictions of the hierarchical galaxy formation models of Benson et al.

We perform a joint likelihood analysis of the power spectra of the 2dF Galaxy Redshift Survey (2dFGRS) and the cosmic microwave background (CMB) anisotropies under the assumptions that the initial fluctuations were adiabatic, Gaussian and well described by power laws with scalar and tensor indices of n_s and n_t. On its own, the 2dFGRS sets tight limits on the parameter combination Omega_m h, but relatively weak limits on the fraction of the cosmic matter density in baryons Omega_b/Omega_m. The CMB anisotropy data alone set poor constraints on the cosmological constant and Hubble constant because of a `geometrical degeneracy' among parameters. Furthermore, if tensor modes are allowed, the CMB data allow a wide range of values for the physical densities in baryons and cold dark matter. Combining the CMB and 2dFGRS data sets helps to break both the geometrical and tensor mode degeneracies. The values of the parameters derived here are consistent with the predictions of the simplest models of inflation, with the baryon density derived from primordial nucleosynthesis and with direct measurements of the Hubble parameter. In particular, we find strong evidence for a positive cosmological constant with a pm 2sigma range of 0.65 < Omega_Lambda < 0.85, completely independently of constraints on Omega_\Lambda derived from Type Ia supernovae.

We use a self-consistent model of galaxy formation and the evolution of the intergalactic medium to study the effects of the reionization of the universe at high redshift on the properties of satellite galaxies like those seen around the Milky Way. Photoionization suppresses the formation of small galaxies, so that surviving satellites are preferentially those that formed before the universe reionized. As a result, the number of satellites expected today is about an order of magnitude smaller than the number inferred by identifying satellites with subhalos in high-resolution simulations of the dark matter. The resulting satellite population has an abundance and a distribution of circular velocities similar to those observed in the Local Group. We explore many other properties of satellite galaxies, including their gas content, metallicity and star formation rate, and find generally good agreement with available data. Our model predicts the existence of many as yet undetected satellites in the Local Group. We quantify their observability in terms of their apparent magnitude and surface brightness and also in terms of their constituent stars. A near-complete census of the Milky Way's satellites would require imaging to V~20 and to a surface brightness fainter than 26 V-band magnitudes per square arcsecond. Satellites with integrated luminosity V=15 should contain of order 100 stars brighter than B=26, with central stellar densities of a few tens per square arcminute. Discovery of a large population of faint satellites would provide a strong test of current models of galaxy formation.

We develop a coupled model for the evolution of the global properties of the intergalactic medium (IGM) and the formation of galaxies, in the presence of a photoionizing background due to stars and quasars. We use this model to predict the thermodynamic history of the IGM when photoionized by galaxies forming in a cold dark matter (CDM) universe. The evolution of the galaxies is calculated using a semi-analytical model, including a detailed treatment of the effects of tidal stripping and dynamical friction on satellite galaxies orbiting inside larger dark matter halos. We include in the model the negative feedback on galaxy formation from the photoionizing background. Photoionization inhibits galaxy formation in low-mass dark matter halos in two ways: (i) heating of the IGM and inhibition of the collapse of gas into dark halos by the IGM pressure, and (ii) reduction in the radiative cooling of gas within halos. The result of our method is a self-consistent model of galaxy formation and the IGM. The IGM is reheated twice (during reionization of HI and HeII), and we find that the star formation rate per unit volume is slightly suppressed after each episode of reheating. We find that galaxies brighter than L_* are mostly unaffected by reionization, while the abundance of faint galaxies is significantly reduced, leading to present-day galaxy luminosity functions with shallow faint end slopes, in good agreement with recent observational data. Reionization also affects other properties of these faint galaxies, in a readily understandable way.

The 2dF Galaxy Redshift Survey has now measured in excess of 160000 galaxy redshifts. This paper presents the power spectrum of the galaxy distribution, calculated using a direct FFT-based technique. We argue that, within the k-space region 0.02<k<0.15 h Mpc^-1, the shape of this spectrum should be close to that of the linear density perturbations convolved with the window function of the survey. This window function and its convolving effect on the power spectrum estimate are analyzed in detail. By convolving model spectra, we are able to fit the power-spectrum data and provide a measure of the matter content of the universe. Our results show that models containing baryon oscillations are mildly preferred over featureless power spectra. Analysis of the data yields 68% confidence limits on the total matter density times the Hubble parameter \Omega_m h = 0.20 +/- 0.03, and the baryon fraction \Omega_b/\Omega_m = 0.15 +/- 0.07, assuming scale-invariant primordial fluctuations.

There is now compelling evidence in favour of the hierarchical structure formation paradigm. Semi-analytic modelling is a powerful tool which allows the formation and evolution of galaxies to be followed in a hierarchical framework. We review some of the latest developments in this area before discussing how such models can help us to interpret observations of the high redshift Universe.

The large-scale structure in the distribution of galaxies is thought to arise from the gravitational instability of small fluctuations in the initial density field of the universe. A key test of this hypothesis is that superclusters of galaxies in the process of formation should generate systematic infall of other galaxies. This would be evident in the pattern of recessional velocities, causing an anisotropy in the inferred spatial clustering of galaxies. Here we report a precise measurement of this clustering, using the redshifts of more than 141,000 galaxies from the two-degree-field galaxy redshift survey. We determine the parameter beta = Omega^{0.6}/b = 0.43 +- 0.07, where Omega is the total mass-density parameter and b is a measure of the `bias' of the luminous galaxies in the survey. Combined with the anisotropy of the cosmic microwave background, our results favour a low-density universe with Omega approximately 0.3.

We follow the evolution of the galaxy population in a Lambda-CDM cosmology by means of high-resolution N-body simulations in which the formation of galaxies and their observable properties are calculated using a semi-analytic model. We display images of the spatial distribution of galaxies in the simulations that illustrate its evolution and provide a qualitative understanding of the processes responsible for various biases that develop. We consider three specific statistical measures of clustering at z=1 and z=0: the correlation length (in real- and redshift-space) of galaxies of different luminosity, the morphology-density relation and the genus curve of the topology of galaxy isodensity surfaces. For galaxies with luminosity below L*, the z=0 correlation length depends very little on the luminosity of the sample, but for brighter galaxies it increases very rapidly, reaching values in excess of 10Mpc/h. The "accelerated" dynamical evolution experienced by galaxies in rich clusters, which is partly responsible for this effect, also results in a strong morphology-density relation. Remarkably, this relation is already well-established at z=1. The genus curves of the galaxies are significantly different from the genus curves of the dark matter but this is not due to genuine topological differences but rather to the sparse sampling of the density field provided by galaxies. The predictions of our model at z=0 will be tested by forthcoming data from the 2dF and Sloan galaxy surveys, and those at z=1 by the DEEP and VIRMOS surveys.

We present a new model for the X-ray properties of the intracluster medium that explicitly includes heating of the gas by the energy released during the evolution of cluster galaxies. We calculate the evolution of clusters by combining the semi-analytic model of galaxy formation of Cole et al. with a simple model for the radial profile of the intracluster gas. We focus on the cluster X-ray luminosity function and on the relation between X-ray temperature and luminosity. We show that cooling alone is not enough to account for the flatness of the observed T-L relation or for the lack of strong redshift evolution in the observed X-ray luminosity function. Gas heating, on the other hand, can solve these two problems: in the \Lambda-CDM cosmology, our model reproduces fairly well the T-L relation and the X-ray luminosity function. Furthermore, it predicts only weak evolution in these two properties out to z=0.5, in agreement with recent observational data. A successful model requires an energy input of 1--2 x 10^49 ergs per solar mass of stars formed. This is comparable to the total energy released by the supernovae associated with the formation of the cluster galaxies. Thus, unless the transfer of supernovae energy to the intracluster gas is very (perhaps unrealistically) efficient, additional sources of energy, such as mechanical energy from AGN winds are required. However, the amplification of an initial energy input by the response of the intracluster medium to protocluster mergers might ease the energy requirements. Our model makes definite predictions for the X-ray properties of groups and clusters at high redshift which may soon be tested with data from the Chandra and Newton satellites.

We combine data from a number of N-body simulations to predict the abundance of dark halos in Cold Dark Matter universes over more than 4 orders of magnitude in mass. A comparison of different simulations suggests that the dominant uncertainty in our results is systematic and is smaller than 10--30% at all masses, depending on the halo definition used. In particular, our ``Hubble Volume'' simulations of \tcdm and \lcdm cosmologies allow the abundance of massive clusters to be predicted with uncertainties well below those expected in all currently planned observational surveys. We show that for a range of CDM cosmologies and for a suitable halo definition, the simulated mass function is almost independent of epoch, of cosmological parameters, and of initial power spectrum when expressed in appropriate variables. This universality is of exactly the kind predicted by the familiar Press-Schechter model, although this model predicts a mass function shape which differs from our numerical results, overestimating the abundance of ``typical'' halos and underestimating that of massive systems.

We use a combination of N-body simulations of the hierarchical clustering of dark matter and semi-analytic modelling of the physics of galaxy formation to probe the relationship between the galaxy distribution and the mass distribution. The pairwise line-of-sight velocity dispersion of galaxies is sensitive to the number of galaxies present in halos of different mass. The galaxy velocity dispersion is ~40% lower than that of the dark matter and is in reasonable agreement with the values measured in the Las Campanas redshift survey by Jing et al. over two decades in pair separation. The origin of this offset is statistical rather than dynamical, and depends upon the relative efficiency of galaxy formation in dark matter halos of different mass. Although the model galaxies and the dark matter have markedly different correlation functions in real space, such biases conspire to cause the redshift space correlation functions to be remarkably similar to each other. The predicted redshift-space galaxy correlation function agrees well with observations. We find no evidence in the model for a dependence of clustering strength on intrinsic galaxy luminosity, unless extremely bright galaxies, two magnitudes brighter than L*, are considered. However, there are significant differences when model galaxies are selected either by morphology or by colour. Early type or red galaxies show a much stronger clustering amplitude than late type or blue galaxies, particularly on small scales, again in good agreement with observations.

We examine the present-day galaxy disk colors and star formation rates (SFRs) in the semi-analytic model of Cole et al. (2000). We find that the fiducial model is a good match to Kennicutt's (1998) observed global star formation law, the color-based metallicity-magnitude correlation and the colors of low-luminosity galaxies. The main limitation of the model, from the point of view of present-day spirals, is that the optical colors of the disks of very luminous spiral galaxies are too blue, even after accounting for the effects of dust.

Recently (Granato, Lacey, Silva et al. 2000, astro-ph/0001308) we have combined our spectrophotometric galaxy evolution code which includes dust reprocessing (GRASIL, Silva et al. 1998) with semi-analytical galaxy formation models (GALFORM, Cole et al. 1999). One of the most characteristic features of the former is that the dust is divided in two main phases: molecular cloud complexes, where stars are assumed to be born, and the diffuse interstellar medium. As a consequence, stellar populations of different ages have different geometrical relationships with the two phases, which is essential in understanding several observed properties of galaxies, in particular those undergoing major episodes of star formation at any redshift. Indeed, our merged GRASIL+GALFORM model reproduces fairly well the SEDs of normal spirals and starbursts from the far-UV to the sub-mm and their internal extinction properties. In particular in the model the observed starburst attenuation law (Calzetti 1999) is accounted for as an effect of geometry of stars and dust, and has nothing to do with the optical properties of dust grains.

We describe the GALFORM semi-analytic model for calculating the formation and evolution of galaxies in hierarchical models. It improves upon, and extends, the Cole et al 1994 model. The model employs a new Monte-Carlo algorithm to follow the merging evolution of dark matter halos with arbitrary mass resolution. It incorporates realistic descriptions of the density profiles of dark matter halos and their gas content; follows the chemical evolution of gas and stars, and the associated production of dust; and includes a detailed calculation of the sizes of disks and spheroids. Wherever possible, our prescriptions for modelling individual physical processes are based on results of numerical simulations. We apply our methods to the LCDM cosmology (Omega_0=0.3, Lambda_0=0.7), and find good agreement with a wide range of properties of the local galaxy population: the B-band and K-band luminosity functions, the distribution of colours for the population as a whole, the ratio of ellipticals to spirals, the distribution of disk sizes, and the current cold gas content of disks. (Abridged)

We present a new evolutionary model for the far-UV to sub-mm properties of the galaxy population. This combines a semi-analytic galaxy formation model based on hierarchical clustering (GALFORM) with a spectro-photometric code which includes dust reprocessing (GRASIL). The former provides the star formation and metal enrichment histories, together with the gas mass and various geometrical parameters, for a representative sample of galaxies formed in different density environments. These quantities allow us to model the SEDs of galaxies, taking into account stellar emission and also dust extinction and re-emission. Two phases are considered for the dust: molecular cloud complexes, where stars are assumed to be born, and the diffuse interstellar medium. The model includes both galaxies forming stars quiescently in disks, and starbursts triggered by galaxy mergers. We test our models against the observed spectro-photometric properties of galaxies in the local Universe. The models reproduce fairly well the SEDs of normal spirals and starbursts, and their internal extinction properties. The starbursts follow the observed relationship between the FIR to UV luminosity ratio and the slope of the UV continuum. They also reproduce the observed starburst attenuation law (Calzetti et al 99). This result is remarkable, because we use a dust mixture which reproduces the Milky Way extinction law. It suggests that the observed attenuation law is related to the geometry of the stars and dust. We compute galaxy luminosity functions over our wide range of wavelengths, which turn out to be in good agreement with observational data. The UV continuum turns out to be a poor star formation indicator for our models, whilst the infrared luminosity is much more reliable.

We have used a combination of high resolution cosmological N-body simulations and semi-analytic modelling of galaxy formation to investigate the processes that determine the spatial distribution of galaxies in cold dark matter (CDM) models. The galaxy distribution depends sensitively on the efficiency with which galaxies form in halos of different mass. In small mass halos, galaxy formation is inhibited by the reheating of cooled gas by feedback processes, whereas in large mass halos, it is inhibited by the long cooling time of the gas. As a result, the mass-to-light ratio of halos has a deep minimum at the halo mass associated with L* galaxies. This leads to a scale dependent bias in the distribution of galaxies relative to the distribution of mass. On large scales, the bias in the galaxy distribution is related in a simple way to the bias in the distribution of massive halos. On small scales, the correlation function is determined by the interplay between various effects including the spatial exclusion of dark matter halos, the distribution function of the number of galaxies occupying a single dark matter halo and, to a lesser extent, dynamical friction. Remarkably, these processes conspire to produce a correlation function in a flat, Omega_0=0.3, CDM model that is close to a power-law over nearly four orders of magnitude in amplitude. This model agrees well with the correlation function of galaxies measured in the APM survey. On small scales, the model galaxies are less strongly clustered than the dark matter whereas on large scales, they trace the occupied halos. Our clustering predictions are robust to changes in the parameters of the galaxy formation model, provided only those models that match the bright end of the galaxy luminosity function are considered. (abridged)

The main ingredients of recent semi-analytic models of galaxy formation are summarised. We present predictions for the galaxy clustering properties of a well specified LCDM model whose parameters are constrained by observed local galaxy properties. We present preliminary predictions for evolution of clustering that can be probed with deep pencil beam surveys.

Given the failure of existing models for redshift-space distortions to provide a highly accurate measure of the beta-parameter, and the ability of forthcoming surveys to obtain data with very low random errors, it becomes necessary to develop better models for these distortions. Here we review the failures of the commonly-used velocity dispersion models and present an empirical method for extracting beta from the quadrupole statistic that has little systematic offset over a wide range of beta and cosmologies. This empirical model is then applied to an ensemble of mock 2dF southern strip surveys, to illustrate the technique and see how accurately we can recover the true value of beta. We compare this treatment with the error we expect to find due only to the finite volume of the survey. We find that non-linear effects reduce the range of scales over which beta can be fitted, and introduce covariances between nearby modes in addition to those introduced by the convolution with the survey window function. The result is that we are only able to constrain beta to a 1-sigma accuracy of 25% (beta=0.55 +/- 0.14 for the cosmological model considered). We explore one possible means of reducing this error, that of cluster collapse, and show that accurate application of this method can greatly reduce the effect of non-linearities, improving the determination of beta. We conclude by demonstrating that, when the contaminating effects of clusters are dealt with, this simple analysis of the full 2dF survey yields beta=0.55 +/- 0.04. For this model this represents a determination of beta to an accuracy of 8% and hence an important constraint on the cosmological density parameter Omega_0.

The formation and evolution of galaxies can be followed in the context of cosmological structure formation using the technique of semi-analytic modelling. We give a brief outline of the features incorporated into the semi-analytic model of Cole etal (1999). We present two examples of model predictions that can be tested using photometric redshift techniques. The first prediction, of the star formation history of the universe, has already been shown to be in broad argeement with the observational estimates. The second prediction, of the evolution of galaxy clustering with redshift, will be addressed with some of the forthcoming deep, multi-filter imaging surveys discussed at this meeting.

Measurements of galaxy clustering are now becoming possible over a range of redshifts out to z=3. We use a semi-analytic model of galaxy formation to compute the expected evolution of the galaxy correlation function with redshift. We illustrate how the degree of clustering evolution is sensitive to the details of the sample selection. For a fixed apparent magnitude limit, galaxies selected at higher redshifts are located in progressively rarer dark matter haloes, compared with the general population of galaxies in place at each redshift. As a result these galaxies are highly biased tracers of the underlying dark matter distribution and exhibit stronger clustering than the dark matter. In general, the correlation length measured in comoving units, decreases at first with increasing redshift, before increasing again at higher redshift. We show that the epsilon-model often used to interpret the angular correlation function of faint galaxies gives an inadequate description of the evolution of clustering, and offers no physical insight into the clustering process. We compare our predictions with those of a simple, popular model in which a one-to-one correspondence between galaxies and dark halos is assumed. Qualitatively, this model reproduces the correct evolutionary behaviour at high redshift, but the quantitative results can be significantly in error. Our theoretical expectations are in good agreement with the high redshift clustering data of Carlberg etal and Postman etal but are higher than the measurements of Le Fevre etal.

We combine information on cosmological parameters from cluster abundances, CMB primordial anisotropies and the IRAS 1.2 Jy galaxy redshift survey. We take as free parameters the present values of the total matter density of the universe, Omega_m, the Hubble parameter, h, sigma_8, and the IRAS biasing factor, b_IRAS. We assume that the universe is spatially flat, with a cosmological constant, and that structure formed from adiabatic initial fluctuations with a Harrison-Zel'dovich power spectrum. The nucleosynthesis value Omega_b=0.019/h^2 is adopted. We use the full three- and four- dimensional likelihood functions for each data set and marginalise these to two- and one- dimensional distributions in a Bayesian way, integrating over the other parameters. It is shown that the three data sets are in excellent agreement, with a best fit point of Omega_m= 1-Omega_Lambda=0.36, h=0.54, sigma_8=0.74, and b_IRAS=1.08. This point is within one sigma of the minimum for each data set alone. Pairs of these data sets have their degeneracies in sufficiently different directions that using only two data sets at a time is sufficient to place good constraints on the cosmological parameters. We show that the results from each of the three possible pairings of the data are also in good agreement. Finally, we combine all three data sets to obtain marginalised 68 per cent confidence intervals of 0.30<Omega_m<0.43, 0.48<h<0.59, 0.69<sigma_8<0.79, and 1.01<b_IRAS<1.16. For the best fit parameters the CMB quadrupole is Q_rms-ps=18.0 uK, Gamma=0.15, Omega_b=0.066 and the age of the universe is 16.7 Gyr.

We describe the 2dF Galaxy Redshift Survey (2dFGRS), and the current status of the observations. In this exploratory paper, we apply a Principal Component Analysis to a preliminary sample of 5869 galaxy spectra and use the two most significant components to split the sample into five spectral classes. These classes are defined by considering visual classifications of a subset of the 2dF spectra, and also by comparing to high quality spectra of local galaxies. We calculate a luminosity function for each of the different classes and find that later-type galaxies have a fainter characteristic magnitude, and a steeper faint-end slope. For the whole sample we find M*=-19.7 (for Omega=1, H_0=100 km/sec/Mpc), alpha=-1.3, phi*=0.017. For class 1 (`early-type') we find M*=-19.6, alpha=-0.7, while for class 5 (`late-type') we find M*=-19.0, alpha=-1.7. The derived 2dF luminosity functions agree well with other recent luminosity function estimates.

We present the basic features and preliminary results of the interface between our spectro-photometric model GRASIL (that calculates galactic SED from the UV to the sub-mm with a detailed computation of dust extinction and thermal reemission) with the semi-analytical galaxy formation model GALFORM (that computes galaxy formation and evolution in the hierarchical scenario, providing the star formation history as an input to our model). With these two models we are able to synthetize simulated samples of a few thousands galaxies suited for statistical studies of galaxy properties to investigate on galaxy formation and evolution. We find good agreement with the available data of SED and luminosity functions.

Semi-analytical models of galaxy formation based on hierarchical clustering now make a wide range of predictions for observable properties of galaxies at low and high redshift. This article concentrates on 2 aspects: (1) Self-consistent modelling of dust absorption predicts a mean UV extinction A_{UV} ~ 1 mag, depending only weakly on redshift, and similar to observational estimates. (2) The models predict that the Lyman-break galaxies found at z ~ 3 should be strongly clustered with a comoving correlation length r_0 = 4-7 Mpc/h (depending on the cosmology), in good agreement with subsequent observational determinations.

We model the two-point correlation function of galaxies in a cold dark matter Universe by combining two powerful theoretical tools - dissipationless N-body simulations of dark matter clustering (specifically the GIF simulations carried out by MPIA and the Virgo Consortium) and semi-analytic modelling of galaxy formation. We construct catalogues of galaxies containing a wide range of information for each galaxy, including magnitudes in various bands, star formation rates, disk and bulge sizes (from the semi-analytic model described by Cole et al 1998) and, importantly, spatial positions and peculiar velocities from the N-body simulation. We then use this information to study the clustering properties of galaxies.

Any realistic theory of galaxy formation must be set in the context of a model for the formation of structure in the universe. We describe a powerful approach -- semi-analytic modelling -- that combines a set of simple rules describing the gas processes involved in galaxy formation with a scheme to follow the hierarchical growth of dark matter haloes. Surprisingly few free parameters are required to specify the model, and these are fixed with reference to a subset of local observational data. The model produces the full star formation history of a galaxy, allowing a wide range of predictions to be made. We review some of the successes of the models, namely the star formation history of the universe and the evolution of galaxy clustering, before focusing our attention on early-type galaxies. We discuss the observational evidence against the classical picture in which early-type galaxies form at some arbitrarily high redshift in a single monolithic collapse and burst of star formation. The alternative scenario in which spheroidal systems are formed by the merger of disk galaxies is outlined. We review some of the predictions of this model, namely the colour-magnitude relation, the faint counts and the evolution of cluster membership.

We present a comprehensive set of mock 2dF and SDSS galaxy redshift surveys constructed from a set of large, high-resolution cosmological N-body simulations. The radial selection functions and geometrical limits of the catalogues mimic those of the genuine surveys. The catalogues span a wide range of cosmologies, including both open and flat universes. In all the models the galaxy distributions are biased so as to approximately reproduce the observed galaxy correlation function on scales of 1--10 Mpc/h In some cases models with a variety of different biasing prescriptions are included. All the mock catalogues are publically available at http://star-www.dur.ac.uk/~cole/mocks/main.html . We expect these catalogues to be a valuable aid in the development of the new algorithms and statistics that will be used to analyse the 2dF and SDSS surveys when they are completed in the next few years. Mock catalogues of the PSCZ survey of IRAS galaxies are also available at the same WWW location.

The high redshift observations of galaxies now becoming available from the Hubble Space Telescope and from large ground based telescopes are opening fresh windows on galaxy formation. Semianalytic models of galaxy formation provide us with a powerful tool to interpret and understand these exciting new data. In this review, we explain the philosophy behind this class of model and outline some of their remarkable successes, focussing our attention on the formation of elliptical galaxies and on the properties of galaxies at high redshift. Now that the recent discovery of star forming galaxies at z=3 has made possible the construction of the cosmic star formation history, which is in good agreement with our model predictions, it appears that a coherent, broadbrush picture of galaxy formation is beginning to emerge.

We use a semi-analytic model of galaxy formation in hierarchical clustering theories to interpret recent data on galaxy formation and evolution, focussing primarily on the recently discovered population of Lyman-break galaxies at $z\simeq 3$. For a variety of cold dark matter (CDM) cosmologies we construct mock galaxy catalogues subject to identical selection criteria to those applied to the real data. We find that the expected number of Lyman-break galaxies is very sensitive to the assumed stellar initial mass function and to the normalization of the primordial power spectrum. For reasonable choices of these and other model parameters, it is possible to reproduce the observed abundance of Lyman-break galaxies in CDM models with $\Omega_0=1$ and with $\Omega_0<1$. These galaxies generally form from rare peaks at high redshift and, as a result, their spatial distribution is strongly biased, with a typical bias parameter, $b\simeq 4$, and a comoving correlation length, $r_0\simeq 4 \mpc$. The typical sizes of these galaxies, $\sim 0.5\kpc$, are substantially smaller than those of present day bright galaxies. In combination with data at lower redshifts, the Lyman-break galaxies can be used to trace the cosmic star formation history. Most galaxies (in our models and in the data) never experience star formation rates in excess of a few solar masses per year. Our models predict that even at $z=5$, the integrated star formation rate is similar to that measured locally, although less than 1% of all the stars have formed prior to this redshift. These galaxies mark the onset of the epoch of galaxy formation that continues into the present day. The basic ingredients of a consistent picture of galaxy formation may well be now in place.

The discovery of a population of young galaxies at an epoch when the universe was about one tenth of its current age has shed new light on the question of when and how galaxies formed. Within the context of popular models this is the population of primeval galaxies that built themselves up to the size of present--day galaxies through the process of repeated mergers. But the recent detection of a large concentration of these primeval galaxies (Steidel et al. 1997) appears to be incompatible with hierarchical clustering models, which generally predict that clusters of this size are fully formed later in time. Here we use a combination of two powerful theoretical techniques --semi-analytic modelling and N-body simulations-- to show that such large concentrations should be quite common in a universe dominated by cold dark matter, and that they are the progenitors of the rich galaxy clusters seen today. We predict the clustering properties of primeval galaxies which should, when compared with data that will be collected in the near future, test our understanding of galaxy formation within the framework of a universe dominated by cold dark matter.

The evolution of the galaxy cluster abundance depends sensitively on the value of the cosmological density parameter, Omega_0. Recent ASCA data are used to quantify this evolution as measured by the X-ray temperature function. A chi^2 minimisation fit to the cumulative temperature function, as well as a maximum likelihood estimate (which requires additional assumptions about cluster luminosities), lead to the estimate Omega_0 \approx 0.45+/-0.2 (1-sigma statistical error). Various systematic uncertainties are considered, none of which enhance significantly the probability that Omega_0=1. These conclusions hold for models with or without a cosmological constant. The statistical uncertainties are at least as large as the individual systematic errors that have been considered here, suggesting that additional temperature measurements of distant clusters will allow an improvement in this estimate. An alternative method that uses the highest redshift clusters to place an upper limit on Omega_0 is also presented and tentatively applied, with the result that Omega_0=1 can be ruled out at the 98 per cent confidence level. Whilst this method does not require a well-defined statistical sample of distant clusters, there are still modelling uncertainties that preclude a firmer conclusion at this time.

I describe a semi-analytical model for the formation and evolution of galaxies in hierarchical clustering models, and its predictions for the properties of the galaxy population at high z. The predictions are found to agree well with the observed properties of the Lyman break galaxies found at z sim 3 by Steidel et al. The models predict that the star formation rate per comoving volume should have peaked at z sim 1-2, which also agrees well with recent observational data.

We use a set of large, high-resolution cosmological N-body simulations to examine the redshift-space distortions of galaxy clustering on scales of order 10-200h^{-1} Mpc. Galaxy redshift surveys currently in progress will, on completion, allow us to measure the quadrupole distortion in the 2-point correlation function, \xi(\sigma,\pi), or its Fourier transform, the power spectrum, P(k,\mu), to a high degree of accuracy. On these scales we typically find a positive quadrupole, as expected for coherent infall onto overdense regions and outflow from underdense regions, but the distortion is substantially weaker than that predicted by pure linear theory. We assess two models that may be regarded as refinements to linear theory, the Zel'dovich approximation and a dispersion model in which the non-linear velocities generated by the formation of virialized groups and clusters are treated as random perturbations to the velocities predicted by linear theory. We find that neither provides an adequate physical description of the clustering pattern. If used to model redshift spacedistortions on scales for 10<\lambda <200 h^{-1}Mpc the estimated value of \beta (\beta=f(\Omega_0)/b where f(\Omega_0) ~ \Omega_0^{0.6} and b is the galaxy bias parameter) is liable to systematic errors of order ten per cent or more. We discuss how such systematics can be avoided by i) development of a more complete model of redshift distortions and ii) the direct use of galaxy catalogues generated from non-linear N-body simulations.

We study the clustering of the mass distribution in COBE-normalized open and flat CDM models using large N-body simulations. With an age of the universe of 14 Gyr (12 Gyr) for the flat (open) models and a baryon density fixed by nucleosynthesis constraints, the observed abundance of rich galaxy clusters leads to tight constraints on the density parameter; 0.25<Omega_0<0.4 (0.4<Omega_0<0.5) for flat (open) models. The upper limits on Omega_0 can be relaxed if one lowers the H_0 and increases the age of the universe, but h<0.25 is required to allow Omega_0=1. The constraints also change if one allows tilted primordial power spectra. An acceptable Omega_0=1 model with h=0.5 can be constructed, with a tilt of n=0.8. We compare the evolved mass correlation functions of the most promising models with those of APM galaxies. The flat models have steep correlation functions at small scales and require galaxies to be antibiased on scales r<8 Mpc/h. The open models require little or no antibias on small scales and a positive bias on large scales. The lack of a positive bias on small scales is difficult to reconcile with the virial analysis of galaxy clusters which, if Omega_0>0.2, implies that galaxies are overabundant in clusters relative to the field. The tilted Omega_0=1 model, on the other hand, does require that galaxies be positively biased on all scales. We also compute the topology of isodensity contours in these models, obtaining theoretical predictions that are less sensitive to galaxy bias.

We review recent developments in theoretical studies of galaxy formation and evolution. In combination with new data from HST, Keck and other large telescopes, numerical and semi-analytic modelling is beginning to build up a coherent picture of galaxy formation. We summarize the current status of modelling of various galactic properties such as the structure of dark matter halos, the galaxy luminosity function, the Tully-Fisher relation, the colour-magnitude relation for ellipticals, the gross morphological properties of galaxies and the counts of faint galaxies as a function of magnitude, redshift and morphology. Many of these properties can be explained, at least at some level, within a broad class of CDM cosmologies, but a number of fundamental issues remain unresolved. We use our semi-analytic model of galaxy formation to interpret the evolutionary status of the Lyman-break galaxies at $z\simeq 3-3.5$ recently discovered by Steidel et al. The abundance and global properties of these objects are compatible with model predictions in a variety of CDM cosmologies, including the standard version. All these models predict mild evolution in the distribution of star formation rates which peaks at around $z\simeq 1$, but is never much larger than it is at present. The Steidel et al. Lyman-break galaxies are among the very first objects in which appreciable star formation is taking place; they thus signal the onset of galaxy formation. We present three example evolutionary histories of Lyman-break galaxies which illustrate that these objects are the precursors of present day, normal, bright ellipticals and spirals.

Tormen and Bertschinger have presented an algorithm which allows the dynamic range of N-body simulations to be extended by adding long-wavelength power to an evolved N-body simulation. This procedure is of considerable interest as it will enable mock galaxy catalogues to be constructed with volumes as large as those of the next generation of galaxy redshift surveys. Their algorithm, however, neglects the coupling between long-wavelength linear modes and short-wavelength non-linear modes. The growth of structure on small scales is coupled to the amplitude of long-wavelength density perturbations via their effect on the local value of the density parameter Omega_0.The effect of neglecting this coupling is quantified using a set of specially tailored N-body simulations. It is shown that the large-scale clustering of objects defined in the evolved density field such as galaxy clusters is strongly underestimated by their algorithm. An adaptation to their algorithm is proposed that, at the expense of additional complexity, remedies the shortcomings of the original one. Methods of constructing biased mock galaxy catalogues which utilise the basic algorithm of Tormen and Bertschinger, but avoid the pitfalls are discussed.

The unprecedented resolution of the refurbished Wide Field and Planetary Camera 2 (WFPC2) on the Hubble Space Telescope (HST) has led to major advances in our understanding of galaxy formation. The high image quality in the Medium Deep Survey and Hubble Deep Field has made it possible, for the first time, to classify faint distant galaxies according to morphological type. These observations have revealed a large population of galaxies classed as irregulars or which show signs of recent merger activity. Their abundance rises steeply with apparent magnitude, providing a likely explanation for the large number of blue galaxies seen at faint magnitudes. We demonstrate that such a population arises naturally in a model in which structure forms hierarchically and which is dynamically dominated by cold dark matter. The number counts of irregular, spiral and elliptical galaxies as a function of magnitude seen in the HST data are well reproduced in this model.We present detailed predictions for the outcome of spectroscopic follow-up observations of the HST surveys. By measuring the redshift distributions of faint galaxies of different morphological types, these programmes will provide a test of the hierarchical galaxy formation paradigm and might distinguish between models with different cosmological parameters.

We use a set of large cosmological N-body simulations to study the internal structure of dark matter haloes which form in scale-free models. We find that the radius r_178 corresponding to a mean interior overdensity of 178 accurately delineates the quasi-static halo interior from the surrounding infalling material, in agreement with the simple spherical collapse model. The interior velocity dispersion correlates with mass, again in good agreement with the spherical collapse model. Interior to the virial radius r_178, the spherical averaged density, circular velocity and velocity dispersion profiles are well fit by a simple 2-parameter analytic model proposed by Navarro etal (1995). This model has density going as 1/r at small radii, steepening to 1/r^3 at large radii, and fits our haloes to the resolution limit of the simulations. The two model parameters, scalelength and mass, are tightly correlated. Lower mass haloes are more centrally concentrated, and so have scalelengths which are a smaller fraction of their virial radius than those of their higher mass counterparts. This reflects the earlier formation times of low mass haloes. The haloes are moderately aspherical, with typical axial ratios 1:0.8:0.65 at their virial radii. These shapes are maintained by an anisotropic velocity dispersion tensor. The median value of the spin parameter is lambda=0.04, with a weak trend for lower lambda at higher halo mass. We also investigate how the halo properties depend on the algorithm used to identify them, using both friends-of-friends and spherical overdensity methods. We find that for groups selected at mean overdensities 100-400 by either method the properties are insensitive to how the haloes are selected, if the halo centre is taken as the position of the most bound particle

The correlation function of galaxy clusters has often been used as a test of cosmological models. A number of assumptions are implicit in the comparison of theoretical expectations to data. Here we use an ensemble of ten large N-body simulations of the standard cold dark matter cosmology to investigate how cluster selection criteria and other uncertain factors influence the cluster correlation function. We consider the effects of varying the definition of a cluster, the mean number density (or equivalently the richness) in a catalogue, and the assumed normalisation of the model; we also examine the importance of redshift space distortions. We implement five different group-finding algorithms and construct cluster catalogues defined by mass, velocity dispersion or a measure of X-ray luminosity. We find that different cluster catalogues yield correlation functions which can differ by substantially more than the statistical errors in any one determination. For a given cluster selection criteria, the correlation length typically varies by $\sim 20\%$ in catalogues spanning the range of intercluster separations covered by the APM and Abell (richness class $\gsim 1$) catalogues. Distortions produced by peculiar velocities in redshift space enhance the correlation function at large separations and lead to a larger clustering length in redshift space than in real space. The sensitivity of the cluster correlation function to various uncertain model assumptions substantially weakens previous conclusions based on the comparison of model predictions with real data. Detailed modelling of cluster selection procedures including the effects of selecting from projected galaxy catalogues is required before the cluster correlation function can be regarded as a high precision constraint on cosmological models.

We present a model for the broad morphological distinction between the disk and spheroidal components of galaxies. Elaborating on the hierarchical clustering scheme of galaxy formation proposed by Cole et al., we assume that galaxies form stars quiescently in a disk until they are disrupted into a spheroidal configuration by mergers. We calculate formation and merging histories, and the evolution in colour, luminosity and morphology of the galaxy populations in different environments. Roughly $ 50\ %$ of ellipticals in our model have undergone a major merger since $z = 0.5$, yet in spite of this we find that cluster ellipticals have colour-magnitude diagrams with remarkably small scatter. The morphological mix of galaxies that become rich cluster members at high redshift is dominated by spiral galaxies, due to the long timescale for galaxy mergers compared with the timescale for cluster assembly at high redshift. The assembly of low redshift clusters is slower, allowing more galaxy mergers to occur in the progenitor halos. As a result $z=0$ rich clusters become E/S0 dominated and we find a ``Butcher-Oemler'' effect that becomes weaker for poorer groups at high redshift. The field luminosity function of red galaxies shows little evolution out to $z\simeq 1$ and the reddest galaxies at these redshifts are as bright as their local counterparts. The blue luminosity function, on the other hand, evolves rapidly with redshift, increasing its characteristic luminosity and becoming steeper at the faint end. These trends are similar to those recently observed in the Canada-France Redshift Survey.

The morphology of galaxy clusters reflects the epoch at which they formed and hence depends on the value of the mean cosmological density, Omega. Recent studies have shown that the distribution of dark matter in clusters can be mapped from analysis of the small distortions in the shapes of background galaxies induced by weak gravitational lensing in the cluster potential. We construct new statistics to quantify the morphology of clusters which are insensitive to limitations in the mass reconstruction procedure. By simulating weak gravitational lensing in artificial clusters grown in numerical simulations of the formation of clusters in three different cosmologies, we obtain distributions of a quadrupole statistic which measures global deviations from spherical symmetry in a cluster. These distributions are very sensitive to the value of Omega_0 and, as a result, lensing observations of a small number of clusters should be sufficient to place broad constraints on Omega_{0} and certainly to distinguish between the extreme values of 0.2 and 1.

Kaiser & Squires have proposed a technique for mapping the dark matter in galaxy clusters using the coherent weak distortion of background galaxy images caused by gravitational lensing. We investigate the effectiveness of this technique under controlled conditions by creating simulated CCD frames containing galaxies lensed by a model cluster, measuring the resulting galaxy shapes, and comparing the reconstructed mass distribution with the original. Typically, the reconstructed surface density is diminished in magnitude when compared to the original. The main cause of this reduced signal is the blurring of galaxy images by atmospheric seeing, but the overall factor by which the reconstructed surface density is reduced depends also on the signal-to-noise ratio in the CCD frame and on both the sizes of galaxy images and the magnitude limit of the sample that is analysed. We propose a method for estimating a multiplicative compensation factor. We test our technique using a lensing cluster drawn from a cosmological N-body simulation with a variety of realistic background galaxy populations and observing conditions. We find that typically the compensation factor is appreciable, 1.4<f<2.2, and varies considerably depending on the observing conditions and sample selection. Our method yields a compensation factor which when used to correct the surface density estimates produces values that are in good agreement with those of the original cluster. Thus weak lensing observations when calibrated using this method yield not only accurate maps of the cluster morphology but also quantitative estimates of the cluster mass distribution.

The population of rich galaxy clusters evolves much more rapidly in a universe with critical density than one with low density, thus offering the possibility of determining the cosmological density parameter, Omega_0. We quantify this evolution using the Press-Schechter formalism which we extend to flat models with a cosmological constant. Using new large N-body simulations, we verify that this formalism accurately predicts the abundance of rich clusters as a function of redshift in various cosmologies. We normalise the models by comparing them to the local abundance of clusters as a function of their X-ray temperature which we rederive from data compiled by Henry & Arnaud. This gives values of the rms density fluctuation in spheres of radius 8 Mpc/h of sigma_8 = (0.50+/- 0.04) Omega_0^{-0.47+0.10 Omega_0} if Lambda_0=0 and sigma_8 = (0.50 +/- 0.04) Omega_0^{-0.53+0.13 Omega_0} if Lambda_0=1-Omega_0. These values depend very weakly on the shape of the power spectrum. We then examine how the distributions of mass, X-ray temperature and Sunyaev-Zel'dovich decrement evolve as a function of Omega_0. We present the expected distributions at z=0.33 and z=0.5 and the predicted number counts of the largest clusters. We find that even at z=0.33, these distributions depend very strongly on Omega_0 and only weakly on Lambda_0. For example, at this redshift, we expect 20 times as many clusters per comoving volume with M>3.5 10^{14} Msol/h and 5 times as many clusters with kT>5 keV if Omega_0=0.3 than if Omega_0=1. The splitting in the integrated counts is enhanced by the larger volume element in low Omega_0 models. There is therefore a real prospect of estimating Omega_0 from forthcoming surveys of intermediate redshift clusters that will determine their masses, X-ray temperatures or SZ decrements.

We give a brief review of current theoretical work in galaxy formation. Recent results from N-body and N-body/hydrodynamic simulations, and from semianalytic modelling are discussed. We present updated versions of some figures from Cole et al (1994). In particular, we show the effect of using the revised stellar population synthesis model of Bruzual and Charlot, which results in a much better match to the observed colour distribution of galaxies than before. We also compare the model output with recently available data on the galaxy luminosity function and the redshift distribution of galaxies in the B and K bands. The form of the Tully-Fisher relation at high redshift predicted by our semi-analytic scheme for galaxy formation is given.

We measure the anisotropy of the redshift-space power spectrum in the 1.2-Jy and QDOT redshift surveys of IRAS-selected galaxies. On large scales, this anisotropy is caused by coherent peculiar motions, and gravitational instability theory predicts a distortion of the power spectrum that depends only on the ratio $\beta \equiv f(\Omega)/b \approx \Omega^{0.6}/b$, where Omega is the cosmological density parameter and $b$ is the bias parameter. On small scales, the distortion is dominated by the random velocity dispersion in non-linear structures. We fit the observed anisotropy with an analytic model that incorporates two parameters, beta, and a small-scale velocity dispersion sigma_v. Tests on N-body simulations show that this model recovers beta quite accurately on the scales accessible to the existing IRAS redshift surveys. Applying our procedure to the 1.2-Jy and QDOT surveys, we find beta=0.52 +/- 0.13 and beta=0.54 +/- 0.3, respectively. These results imply Omega approximately 0.35 if galaxies trace mass, or a bias factor of about 2 if Omega=1.

We predict the observable properties of the galaxy population in popular hierarchical models of galaxy formation. We employ a detailed semianalytic procedure which incorporates the formation and merging of dark matter halos, the shock heating and radiative cooling of gas, self-regulated star formation, the merging of galaxies within dark matter halos, and the spectral evolution of the stellar populations. We contrast the standard CDM cosmogony with variants of the CDM model having either a low value of H_0, or a low value of Omega with or without a cosmological constant. In addition, we compare galaxy formation in these CDM universes with a CHDM model. We find that although the models have some success in remedying the shortcomings of the standard CDM cosmogony, none of these new models produce broad agreement with the whole range of observations. Although the low-Omega and Omega+Lambda=1 CDM models reduce the discrepancy between the predicted and observed Tully-Fisher relations (the main weakness of galaxy formation in standard CDM), these models predict an inverted colour-magnitude relation and do not produce an exponential cut-off at the bright end of the galaxy luminosity function. All of our models predict recent star formation and exhibit galaxy colours bluer than observed, but this problem is far more severe in the CHDM model which produces colours about two magnitudes too blue in B-K. Unlike in the variants of the CDM model in the CHDM case this result is not dependent on our model of stellar feedback, but is instead directly caused by the late epoch of structure formation in this model.

We have made a detailed comparison of the results of large N-body simulations with the analytical description of the merging histories of dark matter halos presented in Lacey & Cole 1993, which is based on an extension of the Press- Schechter method (Bond etal 1991,Bower 1991). We find the analytical predictions for the halo mass function, merger rates and formation times to be remarkably accurate. The analytical model is however expected to apply for arbitrary Omega and general power spectra. Dark matter halos were identified in the simulations using two different methods. For halos selected at mean overdensities 100-200, the analytical mass function was found to provide a good fit to the simulations with a collapse threshold close to that predicted by the spherical collapse model, with a typical error of <30% over a range of 1000 in mass, which is the full dynamic range of our N-body simulations. This was insensitive to the type of filtering used. Over a range of 100 - 1000 in mass, there was also good agreement with the analytical predictions for merger rates, including their dependence on the masses of the two halos involved and the time interval being considered, and for formation times, including the dependence on halo mass and formation epoch. The analytical Press-Schechter mass function and its extension to halo lifetimes and merger rates thus provide a very useful description of the growth of dark matter halos through hierarchical clustering, and should providea valuable tool in studies of the formation and evolution of galaxies and galaxy clusters.

The peculiar velocities of galaxies distort the pattern of galaxy clustering in redshift space, making the redshift space power spectrum anisotropic. In the linear regime, the strength of this distortion depends only on the ratio $\beta \equiv f(\Omega)/b \approx \Omega^{0.6}/b$, where $\Omega$ is the cosmological density parameter and $b$ is the bias parameter. We derive a linear theory estimator for $\beta$ based on the harmonic moments of the redshift space power spectrum. Using N-body simulations, we examine the impact of non-linear gravitational clustering on the power spectrum anisotropy and on our $\beta$-estimator. Non-linear effects can be important out to wavelengths $\lambda \sim 50$Mpc/h or larger; in most cases, they lower the quadrupole moment of the power spectrum and thereby depress the estimate of $\beta$ below the true value. With a sufficiently large redshift survey, the scaling of non-linear effects may allow separate determinations of $\Omega$ and $b$. We describe a practical technique for measuring the anisotropy of the power spectrum from galaxy redshift surveys, and we test the technique on mock catalogues drawn from the N-body simulations. Preliminary application of our methods to the 1.2 Jy IRAS galaxy survey yields $\beta_{est} \sim 0.3-0.4 $ at wavelengths $\lambda \sim 30-40$Mpc/h . Non-linear effects remain important at these scales, so this estimate of $\beta$ is probably lower than the true value.