We carry out ``full-physics'' hydrodynamical simulations of galaxy formation in the normal-branch Dvali-Gabadadze-Porrati (nDGP) braneworld model using a new modified version of the {\sc Arepo} code and the IllustrisTNG galaxy formation model. We simulate two nDGP models (N5 and N1) which represent, respectively, weak and moderate departures from GR, in boxes of sizes $62\,h^{-1}{\rm Mpc}$ and $25\,h^{-1}{\rm Mpc}$ using $2\times512^3$ dark matter particles and initial gas cells. This allows us to explore, for the first time, the impact of baryonic physics on galactic scales in braneworld models of modified gravity and to make predictions on the stellar content of dark matter haloes and galaxy evolution through cosmic time in these models. We find significant differences between the GR and nDGP models in the power spectra and correlation functions of gas, stars and dark matter of up to $\sim 25$ per cent on large scales. Similar to their impact in the standard cosmological model ($\Lambda$CDM), baryonic effects can have a significant influence over the clustering of the overall matter distribution, with a sign that depends on scale. Studying the degeneracy between modified gravity and galactic feedback in these models, we find that these two physical effects on matter clustering can be cleanly disentangled, allowing for a method to accurately predict the matter power spectrum with baryonic effects included, without having to run hydrodynamical simulations. Depending on the braneworld model, we find differences compared with GR of up to $\sim15$ per cent in galaxy properties such as the stellar-to-halo-mass ratio, galaxy stellar mass function, gas fraction and star formation rate density. The amplitude of the fifth force is reduced by the presence of baryons in the very inner part of haloes, but this reduction quickly becomes negligible above $\sim0.1$ times the halo radius.

Luminous red galaxies (LRGs) are one of the key tracers of the large-scale structure of the Universe used by galaxy surveys. Hence, it is important to make accurate predictions for their properties and clustering, including the errors on these statistics. Here, we describe a novel technique which uses the semi-analytical model of galaxy formation {\sc Galform}, embedded in the high-resolution $N$-body Planck-Millennium simulation, to populate a thousand halo catalogues generated using the Parallel-PM $N$-body {\sc glam} code. Our hybrid scheme allows us to make clustering predictions on scales that cannot be modelled in the original $N$-body simulation. LRGs are selected in the redshift range $z=0.6-1$ from the {\sc Galform} output using similar colour-magnitude cuts in the $r$, $z$ and $W1$ bands to those that will be applied in the Dark Energy Spectroscopic Instrument (DESI) survey. We find that the LRG-halo connection is non-trivial, leading to the prediction of a non-standard halo occupation distribution; in particular, the occupation of central galaxies does not reach unity for the most massive haloes, and drops with increasing mass. The {\sc glam} catalogues reproduce the abundance and clustering of the LRGs predicted by {\sc Galform}, and show good agreement with recent measurements of the clustering of DESI-like LRGs using photometric redshifts. We use the \glam{} mocks to compute the covariance matrices for the two-point correlation function and power spectrum of the LRGs and their background dark matter density field, revealing important differences. We also make predictions for the linear-growth rate and the baryon acoustic oscillations distances at $z=0.6$, $0.74$ and $0.93$. All DESI-like LRG catalogues are made publicly available.

We explore how to mitigate the clustering distortions in Lyman-$\alpha$ emitters (LAEs) samples caused by the miss-identification of the Lyman-$\alpha$ (Ly$\alpha$) wavelength in their Ly$\alpha$ line profiles. We use the Ly$\alpha$ line profiles from our previous LAE theoretical model that includes radiative transfer in the interstellar and intergalactic mediums. We introduce a novel approach to measure the systemic redshift of LAEs from their Ly$\alpha$ line using neural networks. In detail, we assume that, for a fraction of the whole LAE population their systemic redshift is determined precisely through other spectral features. We then use this subset to train a neural network that predicts the Ly$\alpha$ wavelength given a Ly$\alpha$ line profile. We test two different training sets: i) the LAEs are selected homogeneously and ii) only the brightest LAEs are selected. In comparison with previous approaches in the literature, our methodology improves significantly both accuracy and precision in determining the Ly$\alpha$ wavelength. In fact, after applying our algorithm in ideal Ly$\alpha$ line profiles, we recover the clustering unperturbed down to 1cMpc/h. Then, we test the performance of our methodology in realistic Ly$\alpha$ line profiles by downgrading their quality. The machine learning techniques work well even if the Ly$\alpha$ line profile quality is decreased considerably. We conclude that LAE surveys such as HETDEX would benefit from determining with high accuracy the systemic redshift of a subpopulation and applying our methodology to estimate the systemic redshift of the rest of the galaxy sample.

The experiment to detect the global epoch of reionization signature (EDGES) collaboration reported the detection of a line at 78MHz in the sky-averaged spectrum due to neutral hydrogen (HI) 21-cm hyperfine absorption of cosmic microwave background photons at $z\sim 17$. This requires that the spin temperature of HI be coupled to the kinetic temperature of the gas at this $z$ through the scattering of Ly-$\alpha$ photons emitted by massive stars. To explain the experimental result, star formation needs to be sufficiently efficient at $z\sim 17$ and this can be used to constrain models in which small-scale structure formation is suppressed (DMF models), either due to dark matter free-streaming or non-standard inflationary dynamics. We combine simulations of structure formation with a simple recipe for star formation to investigate whether these models emit enough Ly-$\alpha$ photons to reproduce the experimental signal for reasonable values of the star formation efficiency, $f_\star$. We find that a thermal warm dark matter (WDM) model with mass $m_\mathrm{WDM}\sim 4.3\,\mathrm{keV}$ is consistent with the timing of the signal for $f_\star\lesssim 2\%$. The exponential growth of structure around $z\sim 17$ in such a model naturally generates a sharp onset of the absorption. A warmer model with $m_\mathrm{WDM}\sim 3\,\mathrm{keV}$ requires a higher star formation efficiency, $f_\star\sim 6\%$, which is a factor of few above predictions of current star formation models and observations of satellites in the Milky Way. However, uncertainties in the process of star formation at these $z$ do not allow to derive strong constrains on such models using 21-cm absorption line. The onset of the 21-cm absorption is generally faster in DMF compared to cold dark matter (CDM), unless some process significantly suppresses star formation in halos with masses below $\sim 10^8\,h^{-1}\,\mathrm{M}_\odot$.

To understand the nature of the accelerated expansion of the Universe, we need to combine constraints on the expansion rate and growth of structure. The growth rate is usually extracted from three dimensional galaxy maps by exploiting the effects of peculiar motions on galaxy clustering. However, theoretical models of the probability distribution function (PDF) of galaxy pairwise peculiar velocities are not accurate enough on small scales to reduce the error on theoretical predictions to the level required to match the precision expected for measurements from future surveys. Here, we improve the modelling of the pairwise velocity distribution by using the Skew-T PDF, which has nonzero skewness and kurtosis. Our model accurately reproduces the redshift-space multipoles (monopole, quadrupole and hexadecapole) predicted by N-body simulations, above scales of about $10\,h^{-1}{\rm Mpc}$. We illustrate how a Taylor expansion of the streaming model can reveal the contributions of the different moments to the clustering multipoles, which are independent of the shape of the velocity PDF. The Taylor expansion explains why the Gaussian streaming model works well in predicting the first two redshift-space multipoles, although the velocity PDF is non-Gaussian even on large scales. Indeed, any PDF with the correct first two moments would produce precise results for the monopole down to scales of about $10\,h^{-1}{\rm Mpc}$, and for the quadrupole down to about $30\,h^{-1}{\rm Mpc}$. An accurate model for the hexadecapole needs to include higher-order moments.

In the near future galaxy surveys will target Lyman alpha emitting galaxies (LAEs) to unveil the nature of the dark energy. It has been suggested that the observability of LAEs is coupled to the large scale properties of the intergalactic medium. Such coupling could introduce distortions into the observed clustering of LAEs, adding a new potential difficulty to the interpretation of upcoming surveys. We present a model of LAEs that incorporates Lyman-alpha radiative transfer processes in the interstellar and intergalactic medium. The model is implemented in the GALFORM semi-analytic model of galaxy of formation and evolution. We find that the radiative transfer inside galaxies produces selection effects over galaxy properties. In particular, observed LAEs tend to have low metallicities and intermediate star formation rates. At low redshift we find no evidence of a correlation between the spatial distribution of LAEs and the intergalactic medium properties. However, at high redshift the LAEs are linked to the line of sight velocity and density gradient of the intergalactic medium. The strength of the coupling depends on the outflow properties of the galaxies and redshift. This effect modifies the clustering of LAEs on large scales, adding non linear features. In particular, our model predicts modifications to the shape and position of the baryon acoustic oscillation peak. This work highlights the importance of including radiative transfer physics in the cosmological analysis of LAEs.

We model the large-scale linear galaxy bias $b_g(x,z)$ as a function of redshift $z$ and observed absolute magnitude threshold $x$ for broadband continuum emission from the far infrared to ultra-violet, as well as for prominent emission lines, such as the H$\alpha$, H$\beta$, Lya and [OII] lines. The modelling relies on the semi-analytic galaxy formation model GALFORM, run on the state-of-the-art $N$-body simulation SURFS with the Planck 2015 cosmology. We find that both the differential bias at observed absolute magnitude $x$ and the cumulative bias for magnitudes brighter than $x$ can be fitted with a five-parameter model: $b_g(x,z)=a + b(1+z)^e(1 + \exp{[(x-c)d]})$. We also find that the bias for the continuum bands follows a very similar form regardless of wavelength due to the mixing of star-forming and quiescent galaxies in a magnitude limited survey. Differences in bias only become apparent when an additional colour separation is included, which suggest extensions to this work could look at different colours at fixed magnitude limits. We test our fitting formula against observations, finding reasonable agreement with some measurements within $1\sigma$ statistical uncertainties, and highlighting areas of improvement. We provide the fitting parameters for various continuum bands, emission lines and intrinsic galaxy properties, enabling a quick estimation of the linear bias in any typical survey of large-scale structure.

Current and future cosmological surveys are targeting star-forming galaxies at z~1 with nebular emission lines. In this work we use a state-of-the-art model of galaxy formation and evolution to explore the large scale environment of star-forming emission line galaxies (ELGs). Model galaxies are selected with cuts such that the samples can be directly compared with the DEEP2, VVDS, eBOSS-SGC and DESI surveys. Their large scale environment have been classified using a velocity-shear-tensor and a tidal-tensor algorithms. Half of model ELGs live in filaments and about a third in sheets. Model ELGs in knots have the largest satellite fractions. We find that the shape of the mean halo occupation distribution of model ELGs varies widely for different large scale environments. To put these results in context, we have compared fixed number density samples of galaxies, generated by either imposing an extra cut in stellar mass, star formation rate (SFR) or L[OII] to the ELGs, or by imposing a single cut in these quantities to the total model galaxy population. For the fixed number density samples, ELGs are close to L[OII] and SFR selected samples for densities above $10^{-4.2}h^{3}{\rm Mpc}^{-3}$. ELGs with an extra cut in stellar mass to fix their number density, present differences in sheets and knots with respect to the other samples. ELGs, SFR and L[OII] selected samples with equal number density, have similar large scale bias but their clustering below separations of $1h^{-1}$Mpc is different.

Telescopes to be launched over the next decade-and-a-half, such as JWST, EUCLID, ATHENA and Lynx, promise to revolutionise the study of the high redshift Universe and greatly advance our understanding of the early stages of galaxy formation. We use a model that follows the evolution of the masses and spins of supermassive black holes (SMBHs) within a semi-analytic model of galaxy formation to make predictions for the Active Galactic Nucleus (AGN) luminosity function at $z\geq7$ in the broadband filters of JWST and EUCLID at near-infrared wavelengths, and ATHENA and Lynx at X-ray energies. The predictions of our model are relatively insensitive to the choice of seed black hole mass, except at the lowest luminosities ($L_{\mathrm{bol}}<10^{43}\mathrm{ergs^{-1}}$) and the highest redshifts ($z>10$). We predict that surveys with these different telescopes will select somewhat different samples of SMBHs, with EUCLID unveiling the most massive, highest accretion rate SMBHs, Lynx the least massive, lowest accretion rate SMBHs, and JWST and ATHENA covering objects inbetween. At $z=7$, we predict that typical detectable SMBHs will have masses, $M_{\mathrm{BH}}\sim10^{5-8}M_{\odot}$, and Eddington normalised mass accretion rates, $\dot{M}/\dot{M}_{\mathrm{Edd}}\sim0.6-2$. The SMBHs will be hosted by galaxies of stellar mass $M_{\star}\sim10^{8-10}M_{\odot}$, and dark matter haloes of mass $M_{\mathrm{halo}}\sim10^{11-12}M_{\odot}$. We predict that the detectable SMBHs at $z=10$ will have slightly smaller black holes, accreting at slightly higher Eddington normalised mass accretion rates, in slightly lower mass host galaxies compared to those at $z=7$, and reside in haloes of mass $M_{\mathrm{halo}}\sim10^{10-11}M_{\odot}$.

We investigate if, for a fixed number density of targets and redshift, there is an optimal way to select a galaxy sample in order to measure the baryon acoustic oscillation (BAO) scale, which is used as a standard ruler to constrain the cosmic expansion. Using the mock galaxy catalogue built by Smith et al. in the Millennium-XXL N-body simulation with a technique to assign galaxies to dark matter haloes based on halo occupation distribution modelling, we consider the clustering of galaxies selected by luminosity, colour and local density. We assess how well the BAO scale can be extracted by fitting a template to the power spectrum measured for each sample. We find that the BAO peak position is recovered equally well for samples defined by luminosity or colour, while there is a bias in the BAO scale recovered for samples defined by density. The BAO position is contracted to smaller scales for the densest galaxy quartile and expanded to large scales for the two least dense galaxy quartiles. For fixed galaxy number density, density-selected samples have higher uncertainties in the recovered BAO scale than luminosity- or colour-selected samples.

We study the evolution of the colour-magnitude relation for galaxies in the VIMOS Public Extragalactic Redshift Survey (VIPERS) by introducing the concept of the bright edge, and use this to derive constraints on the quenching of star formation activity in galaxies over the redshift range $0.5 < z < 1.1$. The bright-edge of the colour-magnitude diagram evolves with little dependence on galaxy colour, and therefore on the amount of star formation taking place in bright galaxies. We modelled this evolution with delayed exponential star formation histories (SFHs), to better understand the time-scale of the turn-off in star formation activity. We show that using SFHs without quenching, the transition from the blue cloud to the red sequence is too slow. This indicates that a scenario purely driven by the consumption of the gas inside each galaxy does not reproduce the observed evolution of the colour-magnitude bright edge. Among the quenching scenarios explored, the one that best matches the observations assumes that galaxies stop their star formation at a randomly selected time with a uniform distribution up to $2.5$ Gyr. We argue that quenching is required over a wide range of stellar masses. Qualitatively similar evolution of the bright edge is found in the predictions of a semi-analytical galaxy formation model, but quantitatively there are marked differences with the observations. This illustrates the utility of the bright edge as a test of galaxy formation models. The evolution changes and no longer matches the observed trend if feedback from heating by active galactic nuclei is turned off.

We present the first application of a variance-based sensitivity analysis (SA) to a model that aims to predict the evolution and properties of the whole galaxy population. SA is a well-established technique in other quantitative sciences, but is a relatively novel tool for the evaluation of astrophysical models. We perform a multi-parameter exploration of the GALFORM semi-analytic galaxy formation model, to compute how sensitive the present-day K-band luminosity function is to varying different model parameters. The parameter space is scanned using a low-discrepancy sampling technique proposed by Saltelli. We first demonstrate the usefulness of the SA approach by varying just two model parameters, one which controls supernova feedback and the other the heating of gas by AGN. The SA analysis matches our physical intuition regarding how these parameters affect the predictions for different parts of the galaxy luminosity function. We then use SA to compute Sobol' sensitivity indices varying seven model parameters, connecting the variance in the model output to the variance in the input parameters. The sensitivity is computed in luminosity bins, allowing us to probe the origin of the model predictions in detail. We discover that the SA correctly identifies the least- and most important parameters. Moreover, the SA also captures the combined responses of varying multiple parameters at the same time. Our study marks a much needed step away from the traditional "one-at-a-time" parameter variation often used in this area and improves the transparency of multi-parameter models of galaxy formation.

We investigate a conceptual modification of the halo occupation distribution approach, using the halos' present-day maximal circular velocity, $\vmax$, as an alternative to halo mass. In particular, using a semi-analytic galaxy formation model applied to the Millennium WMAP7 simulation, we explore the extent that switching to $\vmax$ as the primary halo property incorporates the effects of assembly bias into the formalism. We consider fixed number density galaxy samples ranked by stellar mass and examine the variations in the halo occupation functions with either halo concentration or formation time. We find that using $\vmax$ results in a significant reduction in the occupancy variation of the central galaxies, particularly for concentration. The satellites occupancy variation on the other hand increases in all cases. We find effectively no change in the halo clustering dependence on concentration, for fixed bins of $\vmax$ compared to fixed halo mass. Most crucially, we calculate the impact of assembly bias on galaxy clustering by comparing the amplitude of clustering to that of a shuffled galaxy sample, finding that the level of galaxy assembly bias remains largely unchanged. Our results suggest that while using $\vmax$ as a proxy for halo mass diminishes some of occupancy variations exhibited in the galaxy-halo relation, it is not able to encapsulate the effects of assembly bias potentially present in galaxy clustering. The use of other more complex halo properties, such as $\vpeak$, the peak value of $\vmax$ over the assembly history, provides some improvement and warrants further investigation.

Galaxy cluster counts in bins of mass and redshift have been shown to be a competitive probe to test cosmological models. This method requires an efficient blind detection of clusters from surveys with a well-known selection function and robust mass estimates. The Euclid wide survey will cover 15000 deg$^2$ of the sky in the optical and near-infrared bands, down to magnitude 24 in the $H$-band. The resulting data will make it possible to detect a large number of galaxy clusters spanning a wide-range of masses up to redshift $\sim 2$. This paper presents the final results of the Euclid Cluster Finder Challenge (CFC). The objective of these challenges was to select the cluster detection algorithms that best meet the requirements of the Euclid mission. The final CFC included six independent detection algorithms, based on different techniques, such as photometric redshift tomography, optimal filtering, hierarchical approach, wavelet and friend-of-friends algorithms. These algorithms were blindly applied to a mock galaxy catalog with representative Euclid-like properties. The relative performance of the algorithms was assessed by matching the resulting detections to known clusters in the simulations. Several matching procedures were tested, thus making it possible to estimate the associated systematic effects on completeness to $<3$%. All the tested algorithms are very competitive in terms of performance, with three of them reaching $>80$% completeness for a mean purity of 80% down to masses of $10^{14}$ M$_{\odot}$ and up to redshift $z=2$. Based on these results, two algorithms were selected to be implemented in the Euclid pipeline, the AMICO code, based on matched filtering, and the PZWav code, based on an adaptive wavelet approach. [abridged]

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 test different implementations of the halo occupation distribution (HOD) model to reconstruct the spatial distribution of galaxies as predicted by a publicly available semi-analytical model (SAM). We compare the measured two-point correlation functions of the HOD mock catalogues and the SAM samples to quantify the fidelity of the reconstruction. We use fixed number density galaxy samples selected according to stellar mass or star formation rate (SFR). We develop three different schemes to populate haloes with galaxies with increasing complexity, considering the scatter of the satellite HOD as an additional parameter in the modelling. We first modify the SAM output, removing assembly bias and using a standard Navarro-Frenk-White density profile for the satellite galaxies as the target to reproduce with our HOD mocks. We find that all models give similar reproductions of the two-halo contribution to the clustering signal, but there are differences in the one-halo term. In particular, the HOD mock reproductions work equally well using either the HOD of central and satellites separately or using a model that also accounts for whether or not the haloes contain a central galaxy. We find that the HOD scatter does not have an important impact on the clustering predictions for stellar mass selected samples. For SFR selections, we obtain the most accurate results assuming a negative binomial distribution for the number of satellites in a halo. The scatter in the satellites HOD is a key consideration for HOD mock catalogues that mimic ELG or SFR selected samples in future galaxy surveys.

We forecast the linear bias for H${\rm \alpha}$-emitting galaxies at high redshift. To simulate a Euclid-like and a WFIRST-like survey, we place galaxies into a large-volume dark matter halo lightcone by sampling a library of luminosity-dependent halo occupation distributions (HODs), which is constructed using a physically motivated galaxy formation model. We calibrate the dust attenuation in the lightcones such that they are able to reproduce the H{\alpha} luminosity function or the H{\alpha} cumulative number counts. The angle-averaged galaxy correlation function is computed for each survey in redshift slices of width $\Delta z = 0.2$. In each redshift bin the linear bias can be fitted with a single, scale-independent value that increases with increasing redshift. Fitting for the evolution of linear bias with redshift, we find that our Euclid-like and WFIRST-like surveys are both consistent within error with the relation $b(z) = 0.7z + 0.7$. Our bias forecasts are consistent with bias measurements from the HiZELS survey. We find that the Euclid-like and WFIRST-like surveys yield linear biases that are broadly consistent within error, most likely due to the HOD for the WFIRST-like survey having a steeper power-law slope towards larger halo masses.

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.

Understanding how Active Galactic Nuclei (AGN) evolve through cosmic time allows us to probe the physical processes that control their evolution. We use an updated model for the evolution of masses and spins of supermassive black holes (SMBHs), coupled to the latest version of the semi-analytical model of galaxy formation GALFORM using the Planck cosmology and a high resolution Millennium style dark matter simulation to make predictions for AGN and SMBH properties for $0 < z < 6$. We compare the model to the observed black hole mass function and the SMBH versus galaxy bulge mass relation at $z=0$, and compare the predicted bolometric, hard X-ray, soft X-ray and optical AGN luminosity functions to observations at $z < 6$, and find that the model is in good agreement with the observations. The model predicts that at $z<2$ and $L_{\mathrm{bol}} < 10^{43} \mathrm{ergs^{-1}}$, the AGN luminosity function is dominated by objects accreting in an Advection Dominated Accretion Flow (ADAF) disc state, while at higher redshifts and higher luminosities the dominant contribution is from objects accreting via a thin disc or at super-Eddington rates. The model also predicts that the AGN luminosity function at $z<3$ and $L_{\mathrm{bol}} < 10^{44} \mathrm{ergs^{-1}}$ is dominated by the contribution from AGN fuelled by quiescent hot halo accretion, while at higher luminosities and higher redshifts, the AGN luminosity function is dominated by the contribution from AGN fuelled by starbursts triggered by disc instabilities. We employ this model to predict the evolution of SMBH masses, Eddington ratios, and spins, finding that the median SMBH spin evolves very little for $0<z<6$.

We present a proof of concept of a new galaxy group finder method, Markov graph Clustering (MCL; Van Dongen 2000) that naturally handles probabilistic linking criteria. We introduce a new figure of merit, the variation of information statistic (VI; Meila 2003), used to optimise the free parameter(s) of the MCL algorithm. We explain that the common Friends-of-Friends (FoF) method is a subset of MCL. We test MCL in real space on a realistic mock galaxy catalogue constructed from a N-body simulation using the GALFORM model. With a fixed linking length FoF produces the best group catalogues as quantified by the VI statistic. By making the linking length sensitive to the local galaxy density, the quality of the FoF and MCL group catalogues improve significantly, with MCL being preferred over FoF due to a smaller VI value. The MCL group catalogue recovers accurately the underlying halo multiplicity function at all multiplicities. MCL provides better and more consistent group purity and halo completeness values at all multiplicities than FoF. As MCL allows for probabilistic pairwise connections, it is a promising algorithm to find galaxy groups in photometric surveys.

We use two high resolution N-body simulations, one assuming general relativity and the other the Hu-Sawicki form of $f\left(R\right)$ gravity with $\vert \bar{f}_{R} \vert = 10^{-6}$, to investigate the concentration--formation time relation of dark matter haloes. We stack haloes in logarithmically spaced mass bins to fit median density profiles and extract median formation times. At fixed mass, haloes in modified gravity are more concentrated than those in GR, especially at low masses and at low redshift, and do not follow the concentration--formation time relation seen in GR. We assess the sensitivity of the relation to how concentration and formation time are defined, as well as to the segregation of the halo population by the amount of gravitational screening. We find a clear difference between halo concentrations and assembly histories displayed in modified gravity and those in GR. Existing models for the mass--concentration--redshift relation that have gained success in cold and warm dark matter models require revision in $f\left(R\right)$ gravity.

We examine the evolution of assembly bias using a semi-analytical model of galaxy formation implemented in the Millennium-WMAP7 N-body simulation. We consider fixed number density galaxy samples ranked by stellar mass or star formation rate. We investigate how the clustering of haloes and their galaxy content depend on halo formation time and concentration, and how these relationships evolve with redshift. At $z=0$ the dependences of halo clustering on halo concentration and formation time are similar. However, at higher redshift, halo assembly bias weakens for haloes selected by age, and reverses and increases for haloes selected by concentration. The variation of the halo occupation with concentration and formation time is also similar at $z=0$ and changes at higher redshifts. In this case, the occupancy variation with halo age stays mostly constant with redshift but decreases for concentration. Finally, we look at the evolution of assembly bias reflected in the galaxy distribution by examining the galaxy correlation functions relative to those of shuffled galaxy samples which remove the occupancy variation. This correlation functions ratio monotonically decreases with larger redshift and for lower number density samples, going below unity in some cases, leading to reduced galaxy clustering. While the halo occupation functions themselves vary, the assembly bias trends are similar whether selecting galaxies by stellar mass or star formation rate. Our results provide further insight into the origin and evolution of assembly bias. Our extensive occupation function measurements and fits are publicly available and can be used to create realistic mock catalogues.

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.

Thermal inflation models (which feature two inflationary stages) can display damped primordial curvature power spectra on small scales which generate damped matter fluctuations. For a reasonable choice of parameters, thermal inflation models naturally predict a suppression of the matter power spectrum on galactic and sub-galactic scales, mimicking the effect of warm or interacting dark matter. Matter power spectra in these models are also characterised by an excess of power (w.r.t. the standard $\Lambda$CDM power spectrum) just below the suppression scale. By running a suite of N-body simulations we investigate the non-linear growth of structure in models of thermal inflation. We measure the non-linear matter power spectrum and extract halo statistics, such as the halo mass function, and compare these quantities with those predicted in the standard $\Lambda$CDM model and in other models with damped matter fluctuations. We find that the thermal inflation models considered here produce measurable differences in the matter power spectrum from $\Lambda$CDM at redshifts $z>5$, while the halo mass functions are appreciably different at all redshifts. The halo mass function at $z=0$ for thermal inflation displays an enhancement of around $\sim 20\%$ w.r.t. $\Lambda$CDM and a damping at lower halo masses, with the position of the enhancement depending on the value of the free parameter in the model. The enhancement in the halo mass function (w.r.t. $\Lambda$CDM ) increases with redshift, reaching $\sim 40\%$ at $z=5$. We also study the accuracy of the analytical Press-Schechter approach, using different filters to smooth the density field, to predict halo statistics for thermal inflation. We find that the predictions with the smooth-$k$ filter agree with the simulation results over a wider range of halo masses than is the case with other filters commonly used in the literature.

We study the cosmic evolution of the magnetic fields of a large sample of spiral galaxies in a cosmologically representative volume by employing a semi-analytic galaxy formation model and numerical dynamo solver in tandem. We start by deriving time- and radius-dependent galaxy properties using the galform galaxy formation model, which are then fed into the nonlinear mean-field dynamo equations. These are solved to give the large-scale (mean) field as a function of time and galactocentric radius for a thin disc, assuming axial symmetry. A simple prescription for the evolution of the small-scale (random) magnetic field component is also adopted. We find that, while most massive galaxies are predicted to have large-scale magnetic fields at redshift z=0, a significant fraction of them are expected to contain negligible large-scale field. Our model indicates that, for most of the galaxies containing large-scale magnetic fields today, the mean-field dynamo becomes active at z<3. We compute the radial profiles of pitch angle, and find broad agreement with observational data for nearby galaxies.

We analyse the two-point and marked correlation functions of haloes and galaxies in three variants of the chameleon $f(R)$ gravity model using N-body simulations, and compare to a fiducial $\Lambda$CDM model based on general relativity (GR). Using a halo occupation distribution prescription (HOD) we populate dark matter haloes with galaxies, where the HOD parameters have been tuned such that the galaxy number densities and the real-space galaxy two-point correlation functions in the modified gravity models match those in GR to within $1\sim3\%$. We test the idea that since the behaviour of gravity is dependent on environment, marked correlation functions may display a measurable difference between the models. For this we test marks based on the density field and the Newtonian gravitational potential. We find that the galaxy marked correlation function shows significant differences measured in different models on scales smaller than $r\lesssim 20~h^{-1}$ Mpc. Guided by simulations to identify a suitable mark, this approach could be used as a new probe of the accelerated expansion of the Universe.

Measurements of redshift space distortions (RSD) provide a means to test models of gravity on large-scales. We use mock galaxy catalogues constructed from large N-body simulations of standard and modified gravity models to measure galaxy clustering in redshift space. We focus our attention on two of the most representative and popular families of modified gravity models: the Hu \& Sawicki $f(R)$ gravity and the normal branch of the DGP model. The galaxy catalogues are built using a halo occupation distribution (HOD) prescription with the HOD parameters in the modified gravity models tuned to match with the number density and the real-space clustering of {\sc boss-cmass} galaxies. We employ two approaches to model RSD: the first is based on linear perturbation theory and the second models non-linear effects on small-scales by assuming standard gravity and including biasing and RSD effects. We measure the monopole to real-space correlation function ratio, the quadrupole to monopole ratio, clustering wedges and multipoles of the correlation function and use these statistics to find the constraints on the distortion parameter, $\beta$. We find that the linear model fails to reproduce the N-body simulation results and the true value of $\beta$ on scales $s < 40\Mpch$, while the non-linear modelling of RSD recovers the value of $\beta$ on the scales of interest for all models. RSD on large scales ($s\gtrsim20$-$40\Mpch$) have been found to show significant deviations from the prediction of standard gravity in the DGP models. However, the potential to use RSD to constrain $f(R)$ models is less promising, due to the different screening mechanism in this model.

The recent discovery of gravitational waves marks the culmination of a sequence of successful tests of the general theory of relativity (GR) since its formulation in 1915. Yet these tests remain confined to the scale of stellar systems or the strong gravity regime. A departure from GR on larger, cosmological scales has been advocated by the proponents of modified gravity theories as an alternative to the Cosmological Constant to account for the observed cosmic expansion history. While indistinguishable in these terms by construction, such models on the other hand yield distinct values for the linear growth rate of density perturbations and, as a consequence, for the associated galaxy peculiar velocity field. Measurements of the resulting anisotropy of galaxy clustering, when spectroscopic redshifts are used to derive distances, have thus been proposed as a powerful probe of the validity of GR on cosmological scales. However, despite significant effort in modelling such redshift space distortions, systematic errors remain comparable to current statistical uncertainties. Here, we present the results of a different forward-modelling approach, which fully exploits the sensitivity of the galaxy velocity field to modifications of GR. We use state-of-the-art, high-resolution N-body simulations of a standard GR and a compelling f(R) model, one of GR's simplest variants, to build simulated catalogues of stellar-mass-selected galaxies through a robust match to the Sloan Digital Sky Survey observations. We find that, well within the uncertainty of this technique, f(R) fails to reproduce the observed redshift-space clustering on scales 1-10 Mpc/h. Instead, the standard LCDM GR model agrees impressively well with the data. This result provides a strong confirmation, on cosmological scales, of the robustness of Einstein's general theory of relativity.

A recent comparison of the massive galaxy cluster Abell 2744 with the Millennium XXL (MXXL) N-body simulation has hinted at a tension between the observed substructure distribution and the predictions of LambdaCDM. Follow-up investigations indicated that this could be due to the contribution from the host halo and the subhalo finding algorithm used. To be independent of any subhalo finding algorithm, we therefore investigate the particle data of the MXXL simulation directly. We propose a new method to find substructures in 2D mass maps using a wavelet transform, which treats the simulation and observations equally. Using the same criteria to define a subhalo in observations and simulated data, we find three Abell 2744 analogues in the MXXL simulation. Thus the observations in Abell 2744 are in agreement with the predictions of LambdaCDM. We investigate the reasons for the discrepancy between the results obtained from the SUBFIND and full particle data analyses. We find that this is due to incompatible substructure definitions in observations and SUBFIND.

We study potential systematic effects of assembly bias on cosmological parameter constraints from redshift space distortion measurements. We use a semi-analytic galaxy formation model applied to the Millennium N-body WMAP-7 simulation to study the effects of halo assembly bias on the redshift space distortions of the galaxy correlation function. We look at the pairwise velocities of galaxies living in haloes with concentrations and ages in the upper and lower quintiles, and find that the velocity differences between these are consistent with those reported for real-space clustering analyses, i.e. samples with higher clustering also exhibit stronger infall pairwise motions. This can also be seen in the monopole and quadrupole of the redshift-space correlation function. We find that regardless of the method of measurement, the changes in the $\beta$ parameter due to different secondary halo parameters fully tracks the change in the bias Parameter. Hence, assembly bias does not introduce detectable systematics in the inferred logarithmic growth factor.

We present a gravitational lensing and X-ray analysis of a massive galaxy cluster and its surroundings. The core of MACS\,J0717.5+3745 ($M(R<1\,{\rm Mpc})\sim$\,$2$$\times$$10^{15}\,\msun$, $z$=$0.54$) is already known to contain four merging components. We show that this is surrounded by at least seven additional substructures with masses ranging from $3.8-6.5\times10^{13}\,\msun$, at projected radii $1.6$ to $4.9$\,Mpc. We compare MACS\,J0717 to mock lensing and X-ray observations of similarly rich clusters in cosmological simulations. The low gas fraction of substructures predicted by simulations turns out to match our observed values of $1$--$4\%$. Comparing our data to three similar simulated halos, we infer a typical growth rate and substructure infall velocity. That suggests MACS\,J0717 could evolve into a system similar to, but more massive than, Abell\,2744 by $z=0.31$, and into a $\sim$\,$10^{16}\,\msun$ supercluster by $z=0$. The radial distribution of infalling substructure suggests that merger events are strongly episodic; however we find that the smooth accretion of surrounding material remains the main source of mass growth even for such massive clusters.

We investigate the nonlinear evolution of structure in variants of the standard cosmological model which display damped density fluctuations relative to cold dark matter (e.g. in which cold dark matter is replaced by warm or interacting DM). Using N-body simulations, we address the question of how much information is retained from different scales in the initial linear power spectrum following the nonlinear growth of structure. We run a suite of N-body simulations with different initial linear matter power spectra to show that, once the system undergoes nonlinear evolution, the shape of the linear power spectrum at high wavenumbers does not affect the non-linear power spectrum, while it still matters for the halo mass function. Indeed, we find that linear power spectra which differ from one another only at wavenumbers larger than their half-mode wavenumber give rise to (almost) identical nonlinear power spectra at late times, regardless of the fact that they originate from different models with damped fluctuations. On the other hand, the halo mass function is more sensitive to the form of the linear power spectrum. Exploiting this result, we propose a two parameter model of the transfer function in generic damped scenarios, and show that this parametrisation works as well as the standard three parameter models for the scales on which the linear spectrum is relevant.

We present a mock catalogue for the Physics of the Accelerating Universe Survey (PAUS) and use it to quantify the competitiveness of the narrow band imaging for measuring spectral features and galaxy clustering. The mock agrees with observed number count and redshift distribution data. We demonstrate the importance of including emission lines in the narrow band fluxes. We show that PAUCam has sufficient resolution to measure the strength of the 4000\AA{} break to the nominal PAUS depth. We predict the evolution of a narrow band luminosity function and show how this can be affected by the OII emission line. We introduce new rest frame broad bands (UV and blue) that can be derived directly from the narrow band fluxes. We use these bands along with D4000 and redshift to define galaxy samples and provide predictions for galaxy clustering measurements. We show that systematic errors in the recovery of the projected clustering due to photometric redshift errors in PAUS are significantly smaller than the expected statistical errors. The galaxy clustering on two halo scales can be recovered quantatively without correction, and all qualitative trends seen in the one halo term are recovered. In this analysis mixing between samples reduces the expected contrast between the one halo clustering of red and blue galaxies and demonstrates the importance of a mock catalogue for interpreting galaxy clustering results. The mock catalogue is available on request at https://cosmohub.pic.es/home.

Lyman-{\alpha} emitters (LAEs) are a promising target to probe the large scale structure of the Universe at high redshifts, $z\gtrsim 2$. However, their detection is sensitive to radiative transfer effects that depend on local astrophysical conditions. Thus, modeling the bulk properties of this galaxy population remains challenging for theoretical models. Here we develop a physically-motivated scheme to predict LAEs in cosmological simulations. The escape of Ly{\alpha} photons is computed using a Monte Carlo radiative transfer code which outputs a Ly{\alpha} escape fraction. To speed-up the process of assigning escape fractions to individual galaxies, we employ fitting formulae that approximate the full Monte Carlo results within an accuracy of 10% for a broad range of column densities, gas metallicities and gas bulk velocities. We apply our methodology to the semi-analytical model GALFORM on a large N-body simulation. The Ly{\alpha} photons escape through an outflowing neutral gas medium, implemented assuming different geometries. This results in different predictions for the typical column density and outflow velocities of the LAE population. To understand the impact of radiative transfer on our predictions, we contrast our models against a simple abundance matching assignment. Our full models populate LAEs in less massive haloes than what is obtained with abundance matching. Overall, radiative transfer effects result in better agreement when confronting the properties of LAEs against observational measurements. This suggest that incorporating the effects of Ly{\alpha} radiative transfer in the analysis of this galaxy population, including their clustering, can be important for obtaining an unbiased interpretation of future datasets.

We propose a new filter, a smooth-$k$ space filter, to use in the Press-Schechter approach to model the dark matter halo mass function which overcomes shortcomings of other filters. We test this against the mass function measured in N-body simulations. We find that the commonly used sharp-$k$ filter fails to reproduce the behaviour of the halo mass function at low masses measured from simulations of models with a sharp truncation in the linear power spectrum. We show that the predictions with our new filter agree with the simulation results over a wider range of halo masses for both damped and undamped power spectra than is the case with the sharp-$k$ and real-space top-hat filters.

We use the semi-analytical model of galaxy formation GALFORM to characterise an indirect signature of AGN feedback in the environment of radio galaxies at high redshifts. The predicted environment of radio galaxies is denser than that of radio-quiet galaxies with the same stellar mass. This is consistent with observational results from the CARLA survey. Our model shows that the differences in environment are due to radio galaxies being hosted by dark matter haloes that are ~1.5 dex more massive than those hosting radio-quiet galaxies with the same stellar mass. By running a control-simulation in which AGN feedback is switched-off, we identify AGN feedback as the primary mechanism affecting the build-up of the stellar component of radio galaxies, thus explaining the different environment in radio galaxies and their radio-quiet counterparts. The difference in host halo mass between radio loud and radio quiet galaxies translates into different galaxies populating each environment. We predict a higher fraction of passive galaxies around radio loud galaxies compared to their radio-quiet counterparts. Furthermore, such a high fraction of passive galaxies shapes the predicted infrared luminosity function in the environment of radio galaxies in a way that is consistent with observational findings. Our results suggest that the impact of AGN feedback at high redshifts and environmental mechanisms affecting galaxies in high halo masses can be revealed by studying the environment of radio galaxies, thus providing new constraints on galaxy formation physics at high redshifts.

We investigate the impact of thermal velocities in N-body simulations of structure formation in warm dark matter models. Adopting the commonly used approach of adding thermal velocities, randomly selected from a Fermi-Dirac distribution, to the gravitationally-induced velocities of the simulation particles, we compare the matter and velocity power spectra measured from CDM and WDM simulations, in the latter case with and without thermal velocities. This prescription for adding thermal velocities introduces numerical noise into the initial conditions, which influences structure formation. At early times, the noise affects dramatically the power spectra measured from simulations with thermal velocities, with deviations of the order of $\sim \mathcal{O}(10)$ (in the matter power spectra) and of the order of $\sim \mathcal{O}(10^2)$ (in the velocity power spectra) compared to those extracted from simulations without thermal velocities. At late times, these effects are less pronounced with deviations of less than a few percent. Increasing the resolution of the N-body simulation shifts these discrepancies to higher wavenumbers. We also find that spurious haloes start to appear in simulations which include thermal velocities at a mass that is $\sim$3 times larger than in simulations without thermal velocities.

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.

Emission line galaxies (ELGs) are used in several ongoing and upcoming surveys (SDSS-IV/eBOSS, DESI) as tracers of the dark matter distribution. Using a new galaxy formation model, we explore the characteristics of [OII] emitters, which dominate optical ELG selections at $z\simeq 1$. Model [OII] emitters at $0.5<z<1.5$ are selected to mimic the DEEP2, VVDS, eBOSS and DESI surveys. The luminosity functions of model [OII] emitters are in reasonable agreement with observations. The selected [OII] emitters are hosted by haloes with $M_{\rm halo}\geq 10^{10.3}h^{-1}{\rm M}_{\odot}$, with ~90% of them being central star-forming galaxies. The predicted mean halo occupation distributions of [OII] emitters has a shape typical of that inferred for star-forming galaxies, with the contribution from central galaxies, $\langle N \rangle_{\left[OII\right]\, cen}$, being far from the canonical step function. The $\langle N \rangle_{\left[OII\right]\, cen}$ can be described as the sum of an asymmetric Gaussian for disks and a step function for spheroids, which plateaus below unity. The model [OII] emitters have a clustering bias close to unity, which is below the expectations for eBOSS and DESI ELGs. At $z\sim 1$, a comparison with observed g-band selected galaxy, which are expected to be dominated by [OII] emitters, indicates that our model produces too few [OII] emitters that are satellite galaxies. This suggests the need to revise our modelling of hot gas stripping in satellite galaxies.

We study the dependence of the galaxy content of dark matter halos on large-scale environment and halo formation time using semi-analytic galaxy models applied to the Millennium simulation. We analyze subsamples of halos at the extremes of these distributions and measure the occupation functions for the galaxies they host. We find distinct differences in these occupation functions. The main effect with environment is that central galaxies (and in one model also the satellites) in denser regions start populating lower-mass halos. A similar, but significantly stronger, trend exists with halo age, where early-forming halos are more likely to host central galaxies at lower halo mass. We discuss the origin of these trends and the connection to the stellar mass -- halo mass relation. We find that, at fixed halo mass, older halos and to some extent also halos in dense environments tend to host more massive galaxies. Additionally, we see a reverse trend for the satellite galaxies occupation where early-forming halos have fewer satellites, likely due to having more time for them to merge with the central galaxy. We describe these occupancy variations also in terms of the changes in the occupation function parameters, which can aid in constructing realistic mock galaxy catalogs. Finally, we study the corresponding galaxy auto- and cross-correlation functions of the different samples and elucidate the impact of assembly bias on galaxy clustering. Our results can inform theoretical models of assembly bias and attempts to detect it in the real universe.

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.

We incorporate the non-linear clustering of dark matter halos, as modelled by Jose et al. (2016) into the halo model to better understand the clustering of Lyman break galaxies (LBGs) in the redshift range $z=3-5$. We find that, with this change, the predicted LBG clustering increases significantly on quasi-linear scales ($0.1 \leq r\,/\,h^{-1} \,{\rm Mpc} \leq 10$) compared to that in the linear halo bias model. This in turn results in an increase in the clustering of LBGs by an order of magnitude on angular scales $5" \leq \theta \leq 100"$. Remarkably, the predictions of our new model on the whole remove the systematic discrepancy between the linear halo bias predictions and the observations. The correlation length and large scale galaxy bias of LBGs are found to be significantly higher in the non-linear halo bias model than in the linear halo bias model. The resulting two-point correlation function retains an approximate power-law form in contrast with that computed using the linear halo bias theory. We also find that the non-linear clustering of LBGs increases with increasing luminosity and redshift. Our work emphasizes the importance of using non-linear halo bias in order to model the clustering of high-z galaxies to probe the physics of galaxy formation and extract cosmological parameters reliably.

Early dark energy (EDE) models are a class of quintessence dark energy with a dynamically evolving scalar field which display a small but non-negligible amount of dark energy at the epoch of matter-radiation equality. Compared with a cosmological constant, the presence of dark energy at early times changes the cosmic expansion history and consequently the shape of the linear theory power spectrum and potentially other observables. We constrain the cosmological parameters in the EDE cosmology using recent measurements of the cosmic microwave background and baryon acoustic oscillations. The best-fitting models favour no EDE; here we consider extreme examples which are in mild tension with current observations in order to explore the observational consequences of a maximally allowed amount of EDE. We study the non-linear evolution of cosmic structure in EDE cosmologies using large volume N-body simulations. Many large-scale structure statistics are found to be very similar between the $\Lambda$ cold dark matter ($\Lambda$CDM) and EDE models. We find that the most promising way to distinguish EDE from $\Lambda$CDM is to measure the power spectrum on large scales, where differences of up to 15% are expected.

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.

The massive substructures found in Abell 2744 by Jauzac et al. (2016) present a challenge to the cold dark matter paradigm due to their number and proximity to the cluster centre. We use one of the biggest N-body simulations, the Millennium XXL, to investigate the substructure in a large sample of massive dark matter haloes. A range of effects which influence the comparison with the observations is considered, extending the preliminary evaluation carried out by Jauzac et al. (2016). There are many tens of haloes in the simulation with a total mass comparable with or larger than that of Abell 2744. However, we find no haloes with a number and distribution of massive substructures (> 5 10^13 Msun) that is close to that inferred from the observations of Abell 2744. The application of extreme value statistics suggests that we would need a simulation of at least ten times the volume of the Millennium XXL to find a single dark matter halo with a similar internal structure to Abell 2744. Explaining the distribution of massive substructures in clusters is a new hurdle for hierarchical models to negotiate, which is not weakened by appeals to baryonic physics or uncertainty over the nature of the dark matter particle.

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.

We use the halo occupation distribution (HOD) framework to characterise the predictions from two independent galaxy formation models for the galactic content of dark matter haloes and its evolution with redshift. Our galaxy samples correspond to a range of fixed number densities defined by stellar mass and span $0 \le z \le 3$. We find remarkable similarities between the model predictions. Differences arise at low galaxy number densities which are sensitive to the treatment of heating of the hot halo by active galactic nuclei. The evolution of the form of the HOD can be described in a relatively simple way, and we model each HOD parameter using its value at $z=0$ and an additional evolutionary parameter. In particular, we find that the ratio between the characteristic halo masses for hosting central and satellite galaxies can serve as a sensitive diagnostic for galaxy evolution models. Our results can be used to test and develop empirical studies of galaxy evolution and can facilitate the construction of mock galaxy catalogues for future surveys.

The formation and evolution of galaxies with low neutral atomic hydrogen (HI) masses, M$_{\rm HI}$$<$10$^{8}h^{-2}$M$_{\odot}$, are affected by host dark matter halo mass and photoionisation feedback from the UV background after the end of reionization. We study how the physical processes governing the formation of galaxies with low HI mass are imprinted on the distribution of neutral hydrogen in the Universe using the hierarchical galaxy formation model, GALFORM. We calculate the effect on the correlation function of changing the HI mass detection threshold at redshifts $0 \le z \le 0.5$. We parameterize the clustering as $\xi(r)=(r/r_{0})^{-\gamma}$ and we find that including galaxies with M$_{\rm HI}$$<$10$^{8}h^{-2}$M$_{\odot}$ increases the clustering amplitude $r_{0}$ and slope $\gamma$ compared to samples of higher HI masses. This is due to these galaxies with low HI masses typically being hosted by haloes with masses greater than 10$^{12}{h}^{-1}$M$_{\odot}$, and is in contrast to optically selected surveys for which the inclusion of faint, blue galaxies lowers the clustering amplitude. We show the HI mass function for different host dark matter halo masses and galaxy types (central or satellite) to interpret the values of $r_{0}$ and $\gamma$ of the clustering of HI-selected galaxies. We also predict the contribution of low HI mass galaxies to the 21cm intensity mapping signal. We calculate that a dark matter halo mass resolution better than $\sim$10$^{10}{h}^{-1}$M$_{\odot}$ at redshifts higher than 0.5 is required in order to predict converged 21cm brightness temperature fluctuations.

Using the liminality N-body simulations of Shi et. al., we present the first predictions for galaxy clustering in $f(R)$ gravity using subhalo abundance matching. We find that, for a given galaxy density, even for an $f(R)$ model with $f_{R0}=-10^{-6}$, for which the cold dark matter clustering is very similar to $\Lambda$CDM, the predicted clustering of galaxies in the $f(R)$ model is very different from $\Lambda$CDM. The deviation can be as large as $40\%$ for samples with mean densities close to that of $L_*$ galaxies. This large deviation is testable given the accuracy that future large-scale galaxy surveys aim to achieve. Our result demonstrates that galaxy surveys can provide a stringent test of General Relativity on cosmological scales, which is comparable to the tests from local astrophysical observations.

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 exploit the synergy between low-resolution spectroscopy and photometric redshifts to study environmental effects on galaxy evolution in slitless spectroscopic surveys from space. As a test case, we consider the future Euclid Deep survey (~40deg$^2$), which combines a slitless spectroscopic survey limited at H$\alpha$ flux $\geq5\times 10^{-17}$ erg cm$^{-2}$ s$^{-1}$ and a photometric survey limited in H-band ($H\leq26$). We use Euclid-like galaxy mock catalogues, in which we anchor the photometric redshifts to the 3D galaxy distribution of the available spectroscopic redshifts. We then estimate the local density contrast by counting objects in cylindrical cells with radius from 1 to 10 h$^{-1}$Mpc over the redshift range 0.9<z<1.8. We compare this density field with the one computed in a mock catalogue with the same depth as the Euclid Deep survey (H=26) but without redshift measurement errors. We find that our method successfully separates high from low density environments (the last from the first quintile of the density distribution), with higher efficiency at low redshift and large cell: the fraction of low density regions mistaken by high density peaks is <1% for all scales and redshifts explored, but for scales of 1 h$^{-1}$Mpc for which is a few percent. These results show that we can efficiently study environment in photometric samples if spectroscopic information is available for a smaller sample of objects that sparsely samples the same volume. We demonstrate that these studies are possible in the Euclid Deep survey, i.e. in a redshift range in which environmental effects are different from those observed in the local universe, hence providing new constraints for galaxy evolution models.

We investigate the spatial clustering of dark matter halos, collapsing from $1-4 \sigma$ fluctuations, in the redshift range $0 - 5$ using N-body simulations. The halo bias of high redshift halos ($z \geq 2$) is found to be strongly non-linear and scale-dependent on quasi-linear scales that are larger than their virial radii ($0.5-10$ Mpc/h). However, at lower redshifts, the scale-dependence of non-linear bias is weaker and and is of the order of a few percent on quasi-linear scales at $z \sim 0$. We find that the redshift evolution of the scale dependent bias of dark matter halos can be expressed as a function of four physical parameters: the peak height of halos, the non-linear matter correlation function at the scale of interest, an effective power law index of the {\it rms} linear density fluctuations and the matter density of the universe at the given redshift. This suggests that the scale-dependence of halo bias is not a universal function of the dark matter power spectrum, which is commonly assumed. We provide a fitting function for the scale dependent halo bias as a function of these four parameters. Our fit reproduces the simulation results to an accuracy of better than 4 % over the redshift range $0\leq z \leq 5$. We also extend our model by expressing the non-linear bias as a function of the linear matter correlation function. It is important to incorporate our results into the clustering models of dark matter halos at any redshift, including those hosting early generations of stars and galaxies before reionization.

We present a joint optical/X-ray analysis of the massive galaxy cluster Abell 2744 (z=0.308). Our strong- and weak-lensing analysis within the central region of the cluster, i.e., at R<1Mpc from the brightest cluster galaxy, reveals eight substructures, including the main core. All of these dark-matter halos are detected with a significance of at least 5sigma and feature masses ranging from 0.5 to 1.4x10^{14}Msun within R<150kpc. Merten et al. (2011) and Medezinski et al. (2016) substructures are also detected by us. We measure a slightly higher mass for the main core component than reported previously and attribute the discrepancy to the inclusion of our tightly constrained strong-lensing mass model built on Hubble Frontier Fields data. X-ray data obtained by XMM-Newton reveal four remnant cores, one of them a new detection, and three shocks. Unlike Merten et al. (2011), we find all cores to have both dark and luminous counterparts. A comparison with clusters of similar mass in the MXXL simulations yields no objects with as many massive substructures as observed in Abell 2744, confirming that Abell 2744 is an extreme system. We stress that these properties still do not constitute a challenge to $\Lambda$CDM, as caveats apply to both the simulation and the observations: for instance, the projected mass measurements from gravitational lensing and the limited resolution of the sub-haloes finders. We discuss implications of Abell 2744 for the plausibility of different dark-matter candidates and, finally, measure a new upper limit on the self-interaction cross-section of dark matter of sigma_{DM}<1.28cm2/g(68\% CL), in good agreement with previous results from Harvey et al. (2015).

We investigate the clustering of Lyman-break galaxies (LBGs) at $z\sim4$. Using the hierarchical galaxy formation model GALFORM, we predict, for the first time using a semi-analytical model with feedback from active galactic nuclei (AGN), the angular correlation function (ACF) of LBGs and find agreement within $3\,\sigma$ with new measurements of the ACF from surveys including the Hubble eXtreme Deep Field (XDF) and CANDELS field. Our simulations confirm the conclusion reached using independent models that although the predicted ACFs reproduce the trend of increased clustering with luminosity, the dependence is less strong than observed. We find that for the detection limits of the XDF field central LBGs at $z\sim 4$ predominantly reside in haloes of mass $\sim 10^{11}-10^{12}h^{-1}M_{\rm \odot}$ and that satellites reside in larger haloes of mass $\sim 10^{12}-10^{13}h^{-1}M_{\rm \odot}$. The model predicts fewer bright satellite LBGs at $z\sim4$ than is inferred from measurements of the ACF at small scales. By analysing the halo occupation distribution (HOD) predicted by the model, we find evidence that AGN feedback affects the HOD of central LBGs in massive haloes. This is a new high-redshift test of this important feedback mechanism. We investigate the effect of photometric errors in the observations on the ACF predictions. We find that the observational uncertainty in the galaxy luminosity reduces the clustering amplitude and that this effect increases towards faint galaxies, particularly on small scales. To compare properties of model with observed LBGs this uncertainty must be considered.

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.

A variety of methods have been proposed to define and to quantify galaxy environments. While these techniques work well in general with spectroscopic redshift samples, their application to photometric redshift surveys remains uncertain. To investigate whether galaxy environments can be robustly measured with photo-z samples, we quantify how the density measured with the nearest neighbor approach is affected by photo-z uncertainties by using the Durham mock galaxy catalogs in which the 3D real-space environments and the properties of galaxies are exactly known. Furthermore, we present an optimization scheme in the choice of parameters used in the 2D projected measurements which yield the tightest correlation with respect to the 3D real-space environments. By adopting the optimized parameters in the density measurements, we show that the correlation between the 2D projected optimized density and real-space density can still be revealed, and the color-density relation is also visible out to $z \sim 0.8$ even for a photo-z uncertainty ($\sigma_{\Delta_{z}/(1+z)}$) up to 0.06. We find that at the redshift $0.3 < z < 0.5$ a deep ($i \sim 25$) photometric redshift survey with $\sigma_{\Delta_{z}/(1+z)} = 0.02$ yields a comparable performance of small-scale density measurement to a shallower $i \sim$ 22.5 spectroscopic sample with $\sim$ 10% sampling rate. Finally, we discuss the application of the local density measurements to the Pan-STARRS1 Medium Deep survey, one of the largest deep optical imaging surveys. Using data from $\sim5$ square degrees of survey area, our results show that it is possible to measure local density and to probe the color-density relation with 3$\sigma$ confidence level out to $z \sim 0.8$ in the PS-MDS. The color-density relation, however, quickly degrades for data covering smaller areas.

In the thermal dark matter (DM) paradigm, primordial interactions between DM and Standard Model particles are responsible for the observed DM relic density. In Boehm et al. (2014), we showed that weak-strength interactions between DM and radiation (photons or neutrinos) can erase small-scale density fluctuations, leading to a suppression of the matter power spectrum compared to the collisionless cold DM (CDM) model. This results in fewer DM subhaloes within Milky Way-like DM haloes, implying a reduction in the abundance of satellite galaxies. Here we use very high resolution N-body simulations to measure the dynamics of these subhaloes. We find that when interactions are included, the largest subhaloes are less concentrated than their counterparts in the collisionless CDM model and have rotation curves that match observational data, providing a new solution to the "too big to fail" problem.

The future space missions Euclid and WFIRST-AFTA will use the Halpha emission line to measure the redshifts of tens of millions of galaxies. The Halpha luminosity function at z>0.7 is one of the major sources of uncertainty in forecasting cosmological constraints from these missions. We construct unified empirical models of the Halpha luminosity function spanning the range of redshifts and line luminosities relevant to the redshift surveys proposed with Euclid and WFIRST-AFTA. By fitting to observed luminosity functions from Halpha surveys, we build three models for its evolution. Different fitting methodologies, functional forms for the luminosity function, subsets of the empirical input data, and treatment of systematic errors are considered to explore the robustness of the results. Functional forms and model parameters are provided for all three models, along with the counts and redshift distributions up to z~2.5 for a range of limiting fluxes (F_Halpha>0.5 - 3 x 10^-16 erg cm^-2 s^-1) that are relevant for future space missions. For instance, in the redshift range 0.90<z<1.8, our models predict an available galaxy density in the range 7700--13300 and 2000--4800 deg^-2 respectively at fluxes above F_Halpha>1 and 2 x 10^-16 erg cm^-2 s^-1, and 32000--48000 for F_Halpha>0.5 x 10^-16 erg cm^-2 s^-1 in the extended redshift range 0.40<z<1.8. We also consider the implications of our empirical models for the total Halpha luminosity density of the Universe, and the closely related cosmic star formation history.

The extragalactic gamma-ray and neutrino emission may have a contribution from dark matter (DM) annihilations. In the case of discrepancies between observations and standard predictions, one could infer the DM pair annihilation cross section into cosmic rays by studying the shape of the energy spectrum. So far all analyses of the extragalactic DM signal have assumed the standard cosmological model (LambdaCDM) as the underlying theory. However, there are alternative DM scenarios where the number of low-mass objects is significantly suppressed. Therefore the characteristics of the gamma-ray and neutrino emission in these models may differ from LambdaCDM as a result. Here we show that the extragalactic isotropic signal in these alternative models has a similar energy dependence to that in LambdaCDM, but the overall normalisation is reduced. The similarities between the energy spectra combined with the flux suppression could lead one to misinterpret possible evidence for models beyond LambdaCDM as being due to CDM particles annihilating with a much weaker cross section than expected.

High-redshift galaxy clusters allow us to examine galaxy formation in extreme environments. Here we compile data for 15 $z>1$ galaxy clusters to test the predictions from a state-of-the-art semi-analytical model of galaxy formation. The model gives a good match to the slope and zero-point of the cluster red sequence. The model is able to match the cluster galaxy luminosity function at faint and bright magnitudes, but underestimates the number of galaxies around the break in the cluster luminosity function. We find that simply assuming a weaker dust attenuation improves the model predictions for the cluster galaxy luminosity function, but worsens the predictions for the red sequence at bright magnitudes. Examination of the properties of the bright cluster galaxies suggests that the default dust attenuation is large due to these galaxies having large reservoirs of cold gas as well as small radii. We find that matching the luminosity function and colours of high-redshift cluster galaxies, whilst remaining consistent with local observations, poses a challenge for galaxy formation models.

We use the GALFORM semi-analytical model to study high density regions traced by radio galaxies and quasars at high redshifts. We explore the impact that baryonic physics has upon the properties of galaxies in these environments. Star-forming emission-line galaxies (Ly{\alpha} and H{\alpha} emitters) are used to probe the environments at high redshifts. Radio galaxies are predicted to be hosted by more massive haloes than quasars, and this is imprinted on the amplitude of galaxy overdensities and cross-correlation functions. We find that Ly{\alpha} radiative transfer and AGN feedback indirectly affect the clustering on small scales and also the stellar masses, star- formation rates and gas metallicities of galaxies in dense environments. We also investigate the relation between protoclusters associated with radio galaxies and quasars, and their present- day cluster descendants. The progenitors of massive clusters associated with radio galaxies and quasars allow us to determine an average protocluster size in a simple way. Overdensities within the protoclusters are found to correlate with the halo descendant masses. We present scaling relations that can be applied to observational data. By computing projection effects due to the wavelength resolution of modern spectrographs and narrow-band filters we show that the former have enough spectral resolution to map the structure of protoclusters, whereas the latter can be used to measure the clustering around radio galaxies and quasars over larger scales to determine the mass of dark matter haloes hosting them.

Redshifted 21cm measurements of the structure of ionised regions that grow during reionization promise to provide a new probe of early galaxy and structure formation. One of the challenges of modelling reionization is to account both for the sub-halo scale physics of galaxy formation and the regions of ionization on scales that are many orders of magnitude larger. To bridge this gap we first calculate the statistical relationship between ionizing luminosity and Mpc-scale overdensity using detailed models of galaxy formation computed using relatively small volume - ($\sim$100Mpc/$h$)$^{3}$, high resolution dark matter simulations. We then use a Monte-Carlo technique to apply this relationship to reionization of the intergalactic medium within large volume dark matter simulations - ($>$1Gpc/$h$)$^{3}$. The resulting simulations can be used to address the contribution of very large scale clustering of galaxies to the structure of reionization, and show that volumes larger than 500Mpc/$h$ are required to probe the largest reionization features mid-way through reionization. As an example application of our technique, we demonstrate that the predicted 21cm power spectrum amplitude and gradient could be used to determine the importance of supernovae feedback for early galaxy formation.

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.

We study lensing by voids in Cubic Galileon and Nonlocal gravity cosmologies, which are examples of theories of gravity that modify the lensing potential. We find voids in the dark matter and halo density fields of N-body simulations and compute their lensing signal analytically from the void density profiles, which we show are well fit by a simple analytical formula. In the Cubic Galileon model, the modifications to gravity inside voids are not screened and they approximately double the size of the lensing effects compared to GR. The difference is largely determined by the direct effects of the fifth force on lensing and less so by the modified density profiles. For this model, we also discuss the subtle impact on the force and lensing calculations caused by the screening effects of haloes that exist in and around voids. In the Nonlocal model, the impact of the modified density profiles and the direct modifications to lensing are comparable, but they boost the lensing signal by only $\approx 10\%$, compared with that of GR. Overall, our results suggest that lensing by voids is a promising tool to test models of gravity that modify lensing.

We demonstrate how the properties of a galaxy depend on the mass of its host dark matter subhalo, using two independent models of galaxy formation. For the cases of stellar mass and black hole mass, the median property value displays a monotonic dependence on subhalo mass. The slope of the relation changes for subhalo masses for which heating by active galactic nuclei becomes important. The median property values are predicted to be remarkably similar for central and satellite galaxies. The two models predict considerable scatter around the median property value, though the size of the scatter is model dependent. There is only modest evolution with redshift in the median galaxy property at a fixed subhalo mass. Properties such as cold gas mass and star formation rate, however, are predicted to have a complex dependence on subhalo mass. In these cases subhalo mass is not a good indicator of the value of the galaxy property. We illustrate how the predictions in the galaxy property - subhalo mass plane differ from the assumptions made in empirical models of galaxy clustering by reconstructing the model output using a subhalo abundance matching scheme. In its simplest form, abundance matching generally does not reproduce the clustering predicted by the models, typically resulting in an overprediction of the clustering signal. We show how the basic abundance matching scheme can be extended to reproduce the model predictions more faithfully, which has implications for the analysis of galaxy clustering, particularly for low abundance samples.

We explore the galaxy formation physics governing the low mass end of the HI mass function in the local Universe. Specifically, we predict the effects on the HI mass function of varying i) the strength of photoionisation feedback and the redshift of the end of the epoch of reionization, ii) the cosmology, iii) the supernovae feedback prescription, and iv) the efficiency of star formation. We find that the shape of the low-mass end of the HI mass function is most affected by the critical halo mass below which galaxy formation is suppressed by photoionisation heating of the intergalactic medium. We model the redshift dependence of this critical dark matter halo mass by requiring a match to the low-mass end of the HI mass function. The best fitting critical dark matter halo mass decreases as redshift increases in this model, corresponding to a circular velocity of $\sim 50 \, {\rm km \,s}^{-1}$ at $z=0$, $\sim 30 \, {\rm km\, s}^{-1}$ at $z \sim 1$ and $\sim 12 \, {\rm km \, s}^{-1}$ at $z=6$. We find that an evolving critical halo mass is required to explain both the shape and abundance of galaxies in the HI mass function below $M_{\rm HI} \sim 10^{8} h^{-2} {\rm M_{\odot}}$. The model makes specific predictions for the clustering strength of HI-selected galaxies with HI masses > $10^{6} h^{-2} {\rm M_{\odot}}$ and $> 10^{7} h^{-2} {\rm M_{\odot}}$ and for the relation between the HI and stellar mass contents of galaxies which will be testable with upcoming surveys with the Square Kilometre Array and its pathfinders. We conclude that measurements of the HI mass function at $z \ge 0$ will lead to an improvement in our understanding of the net effect of photoionisation feedback on galaxy formation and evolution.

We present predictions for the two-point correlation function of galaxy clustering as a function of stellar mass, computed using two new versions of the GALFORM semi-analytic galaxy formation model. These models make use of a high resolution, large volume N-body simulation, set in the WMAP7 cosmology. One model uses a universal stellar initial mass function (IMF), while the other assumes different IMFs for quiescent star formation and bursts. Particular consideration is given to how the assumptions required to estimate the stellar masses of observed galaxies (such as the choice of IMF, stellar population synthesis model and dust extinction) influence the perceived dependence of galaxy clustering on stellar mass. Broad-band spectral energy distribution fitting is carried out to estimate stellar masses for the model galaxies in the same manner as in observational studies. We show clear differences between the clustering signals computed using the true and estimated model stellar masses. As such, we highlight the importance of applying our methodology to compare theoretical models to observations. We introduce an alternative scheme for the calculation of the merger timescales for satellite galaxies in GALFORM, which takes into account the dark matter subhalo information from the simulation. This reduces the amplitude of small-scale clustering. The new merger scheme offers improved or similar agreement with observational clustering measurements, over the redshift range 0 < z < 0.7. We find reasonable agreement with clustering measurements from GAMA, but find larger discrepancies for some stellar mass ranges and separation scales with respect to measurements from SDSS and VIPERS, depending on the GALFORM model used.

We show that, in the presence of massive neutrinos, the Galileon gravity model provides a very good fit to the current CMB temperature, CMB lensing and BAO data. This model, which we dub ${\nu} \rm{Galileon}$, when assuming its stable attractor background solution, contains the same set of free parameters as $\Lambda\rm{CDM}$, although it leads to different expansion dynamics and nontrivial gravitational interactions. The data provide compelling evidence ($\gtrsim 6\sigma$) for nonzero neutrino masses, with $\Sigma m_\nu \gtrsim 0.4\ {\rm eV}$ at the $2\sigma$ level. Upcoming precision terrestrial measurements of the absolute neutrino mass scale therefore have the potential to test this model. We show that CMB lensing measurements at multipoles $l \lesssim 40$ will be able to discriminate between the ${\nu} \rm{Galileon}$ and $\Lambda\rm{CDM}$ models. Unlike $\Lambda\rm{CDM}$, the ${\nu} \rm{Galileon}$ model is consistent with local determinations of the Hubble parameter. The presence of massive neutrinos lowers the value of $\sigma_8$ substantially, despite of the enhanced gravitational strength on large scales. Unlike $\Lambda\rm{CDM}$, the ${\nu} \rm{Galileon}$ model predicts a negative ISW effect, which is difficult to reconcile with current observational limits.

We determine the concentration-mass relation of 19 X-ray selected galaxy clusters from the CLASH survey in theories of gravity that directly modify the lensing potential. We model the clusters as NFW haloes and fit their lensing signal, in the Cubic Galileon and Nonlocal gravity models, to the lensing convergence profiles of the clusters. We discuss a number of important issues that need to be taken into account, associated with the use of nonparametric and parametric lensing methods, as well as assumptions about the background cosmology. Our results show that the concentration and mass estimates in the modified gravity models are, within the errorbars, the same as in $\Lambda$CDM. This result demonstrates that, for the Nonlocal model, the modifications to gravity are too weak at the cluster redshifts, and for the Galileon model, the screening mechanism is very efficient inside the cluster radius. However, at distances $\sim \left[2-20\right] {\rm Mpc}/h$ from the cluster center, we find that the surrounding force profiles are enhanced by $\sim20-40\%$ in the Cubic Galileon model. This has an impact on dynamical mass estimates, which means that tests of gravity based on comparisons between lensing and dynamical masses can also be applied to the Cubic Galileon model.

A framework is introduced for coupling the evolution of galactic magnetic fields sustained by the mean-field dynamo with the formation and evolution of galaxies in cold dark matter cosmology. Estimates of the steady-state strength of the large-scale and turbulent magnetic fields from mean-field and fluctuation dynamo models are used together with galaxy properties predicted by semi-analytic models of galaxy formation for a population of spiral galaxies. We find that the field strength is mostly controlled by the evolving gas content of the galaxies. Thus, because of the differences in the implementation of the star formation law, feedback from supernovae and ram-pressure stripping, each of the galaxy formation models considered predicts a distribution of field strengths with unique features. The most prominent of them is the difference in typical magnetic fields strengths obtained for the satellite and central galaxy populations as well as the typical strength of the large-scale magnetic field in galaxies of different mass.

Interactions between dark matter (DM) and radiation (photons or neutrinos) in the early Universe suppress density fluctuations on small mass scales. Here we perform a thorough analysis of structure formation in the fully non-linear regime using N-body simulations for models with DM-radiation interactions and compare the results to a traditional calculation in which DM only interacts gravitationally. Significant differences arise due to the presence of interactions, in terms of the number of low-mass DM haloes and their properties, such as their spin and density profile. These differences are clearly seen even for haloes more massive than the scale on which density fluctuations are suppressed. We also show that semi-analytical descriptions of the matter distribution in the non-linear regime fail to reproduce our numerical results, emphasizing the challenge of predicting structure formation in models with physics beyond collisionless DM.

Theoretical galaxy formation models are an established and powerful tool for interpreting the astrophysical significance of observational data, particularly galaxy surveys. Such models have been utilised with great success by optical surveys such as 2dFGRS and SDSS, but their application to radio surveys of cold gas in galaxies has been limited. In this chapter we describe recent developments in the modelling of the cold gas properties in the models, and how these developments are essential if they are to be applied to cold gas surveys of the kind that will be carried out with the SKA. By linking explicitly a galaxy's star formation rate to the abundance of molecular hydrogen in the galaxy rather than cold gas abundance, as was assumed previously, the latest models reproduce naturally many of the global atomic and molecular hydrogen properties of observed galaxies. We review some of the key results of the latest models and highlight areas where further developments are necessary. We discuss also how model predictions can be most accurately compared with observational data, what challenges we expect when creating synthetic galaxy surveys in the SKA era, and how the SKA can be used to test models of dark matter.

We study the origin of the wide distribution of angles between the angular momenta of the stellar and gas components, $\alpha_{\rm G,S}$, in early-type galaxies (ETGs). We use the GALFORM model of galaxy formation, set in the $\Lambda$ cold dark matter framework, and coupled it with a Monte-Carlo simulation to follow the angular momenta flips driven by matter accretion onto haloes and galaxies. We consider a gas disk to be misaligned with respect to the stellar body if $\alpha_{\rm G,S}>30$~degrees. By assuming that the only sources of misaligments in galaxies are galaxy mergers, we place a lower limit of $2-5$ per cent on the fraction of ETGs with misaligned gas/stellar components. These low fractions are inconsistent with the observed value of $\approx 42\pm 6$ per cent in ATLAS$^{\rm 3D}$. In the more general case, in which smooth gas accretion in addition to galaxy mergers can drive misalignments, our calculation predicts that $\approx 46$ per cent of ETGs have $\alpha_{\rm G,S}>30$~degrees. In this calculation, we find correlations between $\alpha_{\rm G,S}$ and stellar mass, cold gas fraction and star formation rate, such that ETGs with high masses, low cold gas fractions and low star formation rates are more likely to display aligned cold gas and stellar components. We confirm these trends observationally for the first time using ATLAS$^{\rm 3D}$ data. We argue that the high fraction of misaligned gas discs observed in ETGs is mostly due to smooth gas accretion (e.g. cooling from the hot halo of galaxies) which takes place after most of the stellar mass of the galaxy is in place and comes misaligned with respect to the stellar component. Galaxies that have accreted most of their cold gas content prior to the time where most of the stellar mass was in place show aligned components.

We present the [OII] luminosity function measured in the redshift range 0.1<z<1.65 with unprecedented depth and accuracy. Our measurements are based on medium resolution flux-calibrated spectra of emission line galaxies with the FORS2 instrument at VLT and with the SDSS-III/BOSS spectrograph. The FORS2 spectra and the corresponding catalog containing redshifts and line fluxes are released along with this paper. In this work we use a novel method to combine the aforementioned surveys with GAMA, zCOSMOS and VVDS, which have different target selection, producing a consistent weighting scheme to derive the [OII] luminosity function. The measured luminosity function is in good agreement with previous independent estimates. The comparison with two state-of-the-art semi-analytical models is good, which is encouraging for the production of mock catalogs of [OII] flux limited surveys. We observe the bright end evolution over 8.5 Gyr: we measure the decrease of log L* from 42.4 erg/s at redshift 1.44 to 41.2 at redshift 0.165 and we find that the faint end slope flattens when redshift decreases. This measurement confirms the feasibility of the target selection of future baryonic acoustic oscillation surveys aiming at observing [OII] flux limited samples.

We exploit wide-field Ly$alpha$ imaging with Subaru to probe the environment around TN J1338-1942, a powerful radio galaxy with a >100 kpc Ly$\alpha$ halo at z=4.11. We used a sample of Ly$\alpha$ emitters (LAEs) down to $\log(L_{\rm Ly\alpha} [erg\, s^{-1}])\sim 42.8$ to measure the galaxy density around TNJ1338, compared to a control sample from a blank field taken with the same instrument. We found that TNJ1338 resides in a region with a peak overdensity of $\delta_{\rm LAE}=2.8\pm 0.5$ on scales of $8\, h^{-1}\rm Mpc$ (on the sky) and $112\, h^{-1}\rm Mpc$ (line of sight) in comoving coordinates. Adjacent to this overdensity, we found a strong underdensity where virtually no LAEs are detected. We used a semi-analytical model of LAEs derived from the Millennium Simulation to compare our results with theoretical predictions. While the theoretical density distribution is consistent with the blank field, overdense regions such as that around TNJ1338 are very rare, with a number density of $6.4\times 10^{-8}\rm Mpc^{-3}$ (comoving), corresponding to the densest < 0.4 percentile at $z\simeq 4.1$. We also found that the Ly$\alpha$ luminosity function in the TNJ1338 field differs from that in the blank field: the number of bright LAEs ($\log(L_{\rm Ly\alpha}[erg\,s^{-1}]) \gtrsim 43.3$) is enhanced, while the number of fainter LAEs is relatively suppressed. These results suggest that some powerful radio galaxies associated with Ly$\alpha$ nebulae reside in extreme overdensities on $\sim 3$--$6\, \rm Mpc$ scales, where star-formation and AGN activity may be enhanced via frequent galaxy mergers or high rates of gas accretion from the surroundings.

We model the linear and nonlinear growth of large scale structure in the Cubic Galileon gravity model, by running a suite of N-body cosmological simulations using the {\tt ECOSMOG} code. Our simulations include the Vainshtein screening effect, which reconciles the Cubic Galileon model with local tests of gravity. In the linear regime, the amplitude of the matter power spectrum increases by $\sim 25%$ with respect to the standard $\Lambda$CDM model today. The modified expansion rate accounts for $\sim 20%$ of this enhancement, while the fifth force is responsible for only $\sim 5%$. This is because the effective unscreened gravitational strength deviates from standard gravity only at late times, even though it can be twice as large today. In the nonlinear regime ($k \gtrsim 0.1 h\rm{Mpc}^{-1}$), the fifth force leads to only a modest increase ($\lesssim 8%$) in the clustering power on all scales due to the very efficient operation of the Vainshtein mechanism. Such a strong effect is typically not seen in other models with the same screening mechanism. The screening also results in the fifth force increasing the number density of halos by less than 10%, on all mass scales. Our results show that the screening does not ruin the validity of linear theory on large scales which anticipates very strong constraints from galaxy clustering data. We also show that, whilst the model gives an excellent match to CMB data on small angular scales ($l \gtrsim 50$), the predicted integrated Sachs-Wolf effect is in tension with Planck/WMAP results.

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.

We use 80922 galaxies in the Galaxy And Mass Assembly (GAMA) survey to measure the galaxy luminosity function (LF) in different environments over the redshift range 0.04<z<0.26. The depth and size of GAMA allows us to define samples split by colour and redshift to measure the dependence of the LF on environment, redshift and colour. We find that the LF varies smoothly with overdensity, consistent with previous results, with little environmental dependent evolution over the last 3 Gyrs. The modified GALFORM model predictions agree remarkably well with our LFs split by environment, particularly in the most overdense environments. The LFs predicted by the model for both blue and red galaxies are consistent with GAMA for the environments and luminosities at which such galaxies dominate. Discrepancies between the model and the data seen in the faint end of the LF suggest too many faint red galaxies are predicted, which is likely to be due to the over-quenching of satellite galaxies. The excess of bright blue galaxies predicted in underdense regions could be due to the implementation of AGN feedback not being sufficiently effective in the lower mass halos.

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.

The cold dark matter (CDM) model faces persistent challenges on small scales. In particular, taken at face value, the model significantly overestimates the number of satellite galaxies around the Milky Way. Attempts to solve this problem remain open to debate and have even led some to abandon CDM altogether. However, current simulations are limited by the assumption that dark matter feels only gravity. Here, we show that including interactions between CDM and radiation (photons or neutrinos) leads to a dramatic reduction in the number of satellite galaxies, alleviating the Milky Way satellite problem and indicating that physics beyond gravity may be essential to make accurate predictions of structure formation on small scales. The methodology introduced here gives constraints on dark matter interactions that are significantly improved over those from the cosmic microwave background.

We study the nonlinear growth of structure in nonlocal gravity models with the aid of N-body simulation and the spherical collapse and halo models. We focus on a model in which the inverse-squared of the d'Alembertian operator acts on the Ricci scalar in the action. For fixed cosmological parameters, this model differs from $\Lambda{\rm CDM}$ by having a lower late-time expansion rate and an enhanced and time-dependent gravitational strength ($\sim 6\%$ larger today). Compared to $\Lambda{\rm CDM}$ today, in the nonlocal model, massive haloes are slightly more abundant (by $\sim 10\%$ at $M \sim 10^{14} M_{\odot}/h$) and concentrated ($\approx 8\%$ enhancement over a range of mass scales), but their linear bias remains almost unchanged. We find that the Sheth-Tormen formalism describes the mass function and halo bias very well, with little need for recalibration of free parameters. The fitting of the halo concentrations is however essential to ensure the good performance of the halo model on small scales. For $k \gtrsim 1 h/{\rm Mpc}$, the amplitude of the nonlinear matter and velocity divergence power spectra exhibits a modest enhancement of $\sim 12\%$ to $15\%$, compared to $\Lambda{\rm CDM}$ today. This suggests that this model might only be distinguishable from $\Lambda{\rm CDM}$ by future observational missions. We point out that the absence of a screening mechanism may lead to tensions with Solar System tests due to local time variations of the gravitational strength, although this is subject to assumptions about the local time evolution of background averaged quantities.

The non-linear, scale-dependent bias in the mass distribution of galaxies and the underlying dark matter is a key systematic affecting the extraction of cosmological parameters from galaxy clustering. Using 95 million halos from the Millennium-XXL N-body simulation, we find that the mass bias is scale independent only for $k<0.1 h{\rm Mpc}^{-1}$ today ($z=0$) and for $k<0.2 h{\rm Mpc}^{-1}$ at $z=0.7$. We test analytic halo bias models against our simulation measurements and find that the model of Tinker et al. 2005 is accurate to better then 5% at $z=0$. However, the simulation results are better fit by an ellipsoidal collapse model at $z=0.7$. We highlight, for the first time, another potentially serious systematic due to a sampling bias in the halo velocity divergence power spectra which will affect the comparison between observations and any redshift space distortion model which assumes dark matter velocity statistics with no velocity bias. By measuring the velocity divergence power spectra for different sized halo samples, we find that there is a significant bias which increases with decreasing number density. This bias is approximately 20% at $k=0.1h$Mpc$^{-1}$ for a halo sample of number density $\bar{n} = 10^{-3} (h/$Mpc$)^3$ at both $z=0$ and $z=0.7$ for the velocity divergence auto power spectrum. Given the importance of redshift space distortions as a probe of dark energy and the on-going major effort to advance models for the clustering signal in redshift space, our results show this velocity bias introduces another systematic, alongside scale-dependent halo mass bias, which cannot be neglected.

We test the cosmological implications of studying galaxy clustering using a tomographic approach, by computing the galaxy two-point angular correlation function $\omega(\theta)$ in thin redshift shells using a spectroscopic-redshift galaxy survey. The advantages of this procedure are that it is not necessary to assume a fiducial cosmology in order to convert measured angular positions and redshifts into distances, and that it gives several (less accurate) measurements of the angular diameter distance $D_\rm{A}(z)$ instead of only one (more precise) measurement of the effective average distance $D_\rm{V}(z)$, which results in better constraints on the expansion history of the Universe. We test our model for $\omega(\theta)$ and its covariance matrix against a set of mock galaxy catalogues and show that this technique is able to extract unbiased cosmological constraints. Also, assuming the best-fit $\Lambda$CDM cosmology from the cosmic microwave background measurements from the Planck satellite, we forecast the result of applying this tomographic approach to the final Baryon Oscillation Spectroscopic Survey catalogue in combination with Planck for three flat cosmological models, and compare them with the expected results of the isotropic baryon acoustic oscillation (BAO) measurements post-reconstruction on the same galaxy catalogue combined with Planck. While BAOs are more accurate for constraining cosmological parameters for the standard $\Lambda$CDM model, the tomographic technique gives better results when we allow the dark energy equation of state $w_{DE}$ to deviate from $-1$, resulting in a performance similar to BAOs in the case of a constant value of $w_{DE}$, and a moderate improvement in the case of a time-dependent value of $w_{DE}$, increasing the value of the Figure-of-Merit in the $w_0-w_a$ plane up to $15\%$.

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.

We use the latest CMB data from Planck, together with BAO measurements, to constrain the full parameter space of Galileon gravity. We constrain separately the three main branches of the theory known as the Cubic, Quartic and Quintic models, and find that all yield a very good fit to these data. Unlike in $\Lambda{\rm CDM}$, the Galileon model constraints are compatible with local determinations of the Hubble parameter and predict nonzero neutrino masses at over $5\sigma$ significance. We also identify that the low-$l$ part of the CMB lensing spectrum may be able to distinguish between $\Lambda{\rm CDM}$ and Galileon models. In the Cubic model, the lensing potential deepens at late times on sub-horizon scales, which is at odds with the current observational suggestion of a positive ISW effect. Compared to $\Lambda$CDM, the Quartic and Quintic models predict less ISW power in the low-$l$ region of the CMB temperature spectrum, and as such are slightly preferred by the Planck data. We illustrate that residual local modifications to gravity in the Quartic and Quintic models may render the Cubic model as the only branch of Galileon gravity that passes Solar System tests.

We study the atomic (HI) and molecular hydrogen (H2) contents of early-type galaxies (ETGs) and their gas sources using the GALFORM model of galaxy formation. This model uses a self-consistent calculation of the star formation rate (SFR), which depends on the H2 content of galaxies. We first present a new analysis of HIPASS and ATLAS3D surveys, with special emphasis on ETGs. The model predicts HI and H2 contents of ETGs in agreement with the observations from these surveys only if partial ram pressure stripping of the hot gas is included, showing that observations of neutral gas in `quenched' galaxies place stringent constraints on the treatment of the hot gas in satellites. We find that ~90% of ETGs at z=0 have neutral gas contents supplied by radiative cooling from their hot halos, 8% were supplied by gas accretion from minor mergers that took place in the last 1Gyr, while 2% were supplied by mass loss from old stars. The model predicts neutral gas fractions strongly decreasing with increasing bulge fraction. This is due to the impeded disk regeneration in ETGs, resulting from both active galactic nuclei feedback and environmental quenching by partial ram pressure stripping of the hot gas.

We investigate the performance of semi-analytical modelling of large-scale structure in Galileon gravity cosmologies using results from N-body simulations. We focus on the Cubic and Quartic Galileon models that provide a reasonable fit to CMB, SNIa and BAO data. We demonstrate that the Sheth-Tormen mass function and linear halo bias can be calibrated to provide a very good fit to our simulation results. We also find that the halo concentration-mass relation is well fitted by a power law. The nonlinear matter power spectrum computed in the halo model approach is found to be inaccurate in the mildly nonlinear regime, but captures reasonably well the effects of the Vainshtein screening mechanism on small scales. In the Cubic model, the screening mechanism hides essentially all of the effects of the fifth force inside haloes. In the case of the Quartic model, the screening mechanism leaves behind residual modifications to gravity, which make the effective gravitational strength time-varying and smaller than the standard value. Compared to normal gravity, this causes a deficiency of massive haloes and leads to a weaker matter clustering on small scales. For both models, we show that there are realistic halo occupation distributions of Luminous Red Galaxies that can match both the observed large-scale clustering amplitude and the number density of these galaxies.

We study spherical collapse in the Quartic and Quintic Covariant Galileon gravity models within the framework of the excursion set formalism. We derive the nonlinear spherically symmetric equations in the quasi-static and weak-field limits, focusing on model parameters that fit current CMB, SNIa and BAO data. We demonstrate that the equations of the Quintic model do not admit physical solutions of the fifth force in high density regions, which prevents the study of structure formation in this model. For the Quartic model, we show that the effective gravitational strength deviates from unity at late times ($z \lesssim 1$), becoming larger if the density is low, but smaller if the density is high. This shows that the Vainshtein mechanism at high densities is not enough to screen all of the modifications of gravity. This makes halos that collapse at $z \lesssim 1$ feel an overall weaker gravity, which suppresses halo formation. However, the matter density in the Quartic model is higher than in standard $\Lambda$CDM, which boosts structure formation and dominates over the effect of the weaker gravity. In the Quartic model there is a significant overabundance of high-mass halos relative to $\Lambda$CDM. Dark matter halos are also less biased than in $\Lambda$CDM, with the difference increasing appreciably with halo mass. However, our results suggest that the bias may not be small enough to fully reconcile the predicted matter power spectrum with LRG clustering data.

We study the contribution of galaxies with different properties to the global densities of star formation rate (SFR), atomic (HI) and molecular hydrogen (H2) as a function of redshift. We use the GALFORM model of galaxy formation, which is set in the LCDM framework. This model includes a self-consistent calculation of the SFR, which depends on the H2 content of galaxies. The predicted SFR density and how much of this is contributed by galaxies with different stellar masses and infrared luminosities are in agreement with observations. The model predicts a modest evolution of the HI density at z<3, which is also in agreement with the observations. The HI density is predicted to be always dominated by galaxies with SFR<1Msun/yr. This contrasts with the H2 density, which is predicted to be dominated by galaxies with SFR>10Msun/yr. Current high-redshift galaxy surveys are limited to detect carbon monoxide in galaxies with SFR>30Msun/yr, which in our model make up, at most, 20% of the H2 in the universe. In terms of stellar mass, the predicted H2 density is dominated by massive galaxies, Mstellar>10^10Msun, while the HI density is dominated by low mass galaxies, Mstellar<10^9Msun. In the context of upcoming neutral gas surveys, we suggest that the faint nature of the galaxies dominating the HI content of the Universe will hamper the identification of optical counterparts, while for H2, we expect follow up observations of molecular emission lines of already existing galaxy catalogues to be able to uncover the H2 density of the Universe.

We present a new release of the GALFORM semi-analytical model of galaxy formation and evolution, which exploits a Millennium Simulation-class N-body run performed with the WMAP7 cosmology. We use this new model to study the impact of the choice of stellar population synthesis (SPS) model on the predicted evolution of the galaxy luminosity function. The semi-analytical model is run using seven different SPS models. In each case we obtain the rest-frame luminosity function in the far-ultra-violet, optical and near-infrared (NIR) wavelength ranges. We find that both the predicted rest-frame ultra-violet and optical luminosity function are insensitive to the choice of SPS model. However, we find that the predicted evolution of the rest-frame NIR luminosity function depends strongly on the treatment of the thermally pulsating asymptotic giant branch (TP-AGB) stellar phase in the SPS models, with differences larger than a factor of 2 for model galaxies brighter than $M_{\rm AB}(K)-5$log$h<-22$ ($\sim$L$_*$ for $0\leq z\leq 1.5$). We have also explored the predicted number counts of galaxies, finding remarkable agreement between the results with different choices of SPS model, except when selecting galaxies with very red optical-NIR colours. The predicted number counts of these extremely red galaxies appear to be more affected by the treatment of star formation in disks than by the treatment of TP-AGB stars in the SPS models.

We measure the angular clustering of 33 415 extremely red objects (EROs) in the Elais-N1 field covering 5.33 deg$^{2}$, which cover the redshift range $z=0.8$ to $2$. This sample was made by merging the UKIDSS Deep eXtragalactic Survey (DXS) with the optical Subaru and Pan-STARRS PS1 datasets. We confirm the existence of a clear break in the angular correlation function at $\sim 0.02^{\circ}$ corresponding to $1 h^{-1}$ Mpc at $z\sim1$. We find that redder or brighter EROs are more clustered than bluer or fainter ones. Halo Occupation Distribution (HOD) model fits imply that the average mass of dark matter haloes which host EROs is over $10^{13} h^{-1} M_{\odot}$ and that EROs have a bias ranging from 2.7 to 3.5. Compared to EROs at $z\sim1.1$, at $z\sim1.5$ EROs have a higher bias and fewer are expected to be satellite galaxies. Furthermore, EROs reside in similar dark matter haloes to those that host $10^{11.0} M_{\odot}<M_{*}<10^{11.5} M_{\odot}$ galaxies. We compare our new measurement and HOD fits with the predictions of the GALFORM semi-analytical galaxy formation model. Overall, the clustering predicted by GALFORM gives an encouraging match to our results. However, compared to our deductions from the measurements, GALFORM puts EROs into lower mass haloes and predicts that a larger fraction of EROs are satellite galaxies. This suggests that the treatment of gas cooling may need to be revised in the model. Our analysis illustrates the potential of clustering analyses to provide observational constraints on theoretical models of galaxy formation.

We develop a numerical algorithm to solve the high-order nonlinear derivative-coupling equation associated with the quartic Galileon model, and implement it in a modified version of the RAMSES N-body code to study the effect of the Galileon field on the large-scale matter clustering. The algorithm is tested for several matter field configurations with different symmetries, and works very well. This enables us to perform the first simulations for a quartic Galileon model which provides a good fit to the cosmic microwave background (CMB) anisotropy, supernovae and baryonic acoustic oscillations (BAO) data. Our result shows that the Vainshtein mechanism in this model is very efficient in suppressing the spatial variations of the scalar field. However, the time variation of the effective Newtonian constant caused by the curvature coupling of the Galileon field cannot be suppressed by the Vainshtein mechanism. This leads to a significant weakening of the strength of gravity in high-density regions at late times, and therefore a weaker matter clustering on small scales. We also find that without the Vainshtein mechanism the model would have behaved in a completely different way, which shows the crucial role played by nonlinearities in modified gravity theories and the importance of performing self-consistent N-body simulations for these theories.

We present the first constraints on the full parameter space of the Galileon modified gravity model, considering both the cosmological parameters and the coefficients which specify the additional terms in the Lagrangian due to the Galileon field, which we call the Galileon parameters. We use the latest cosmic microwave background measurements, along with distance measurements from supernovae and baryonic acoustic oscillations, performing a Monte Carlo Markov Chain exploration of the 9-dimensional parameter space. The integrated Sachs-Wolfe signal can be very different in Galileon models compared to standard gravity, making it essential to use the full CMB data rather than the CMB distance priors. We demonstrate that meaningful constraints are only possible in the Galileon parameter space after taking advantage of a scaling degeneracy. We find that the Galileon model can fit the WMAP 9-year results better than the standard \Lambda-Cold Dark Matter model, but gives a slightly worse fit overall once lower redshift distance measurements are included. The best-fitting cosmological parameters (e.g. matter density, scalar spectral index, fluctuation amplitude) can differ by more than 2\sigma\ in the Galileon model compared with \Lambda CDM. We highlight other potential constraints of the Galileon model using galaxy clustering and weak lensing measurements.

We present a scheme to extend the halo mass resolution of N-body simulations of the hierarchical clustering of dark matter. The method uses the density field of the simulation to predict the number of sub-resolution dark matter haloes expected in different regions. The technique requires as input the abundance of haloes of a given mass and their average clustering, as expressed through the linear and higher order bias factors. These quantities can be computed analytically or, more accurately, derived from a higher resolution simulation as done here. Our method can recover the abundance and clustering in real- and redshift-space of haloes with mass below $\sim 7.5 \times 10^{13}h^{-1}M_{\odot}$ at $z=0$ to better than 10%. We demonstrate the technique by applying it to an ensemble of 50 low resolution, large-volume $N$-body simulations to compute the correlation function and covariance matrix of luminous red galaxies (LRGs). The limited resolution of the original simulations results in them resolving just two thirds of the LRG population. We extend the resolution of the simulations by a factor of 30 in halo mass in order to recover all LRGs. With existing simulations it is possible to generate a halo catalogue equivalent to that which would be obtained from a $N$-body simulation using more than 20 trillion particles; a direct simulation of this size is likely to remain unachievable for many years. Using our method it is now feasible to build the large numbers of high-resolution large volume mock galaxy catalogues required to compute the covariance matrices necessary to analyse upcoming galaxy surveys designed to probe dark energy.

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.

We present a very large high-resolution cosmological N-body simulation, the Millennium-XXL or MXXL, which uses 303 billion particles to represent the formation of dark matter structures throughout a 4.1Gpc box in a LambdaCDM cosmology. We create sky maps and identify large samples of galaxy clusters using surrogates for four different observables: richness estimated from galaxy surveys, X-ray luminosity, integrated Sunyaev-Zeldovich signal, and lensing mass. The unprecedented combination of volume and resolution allows us to explore in detail how these observables scale with each other and with cluster mass. The scatter correlates between different mass-observable relations because of common sensitivities to the internal structure, orientation and environment of clusters, as well as to line-of-sight superposition of uncorrelated structure. We show that this can account for the apparent discrepancies uncovered recently between the mean thermal SZ signals measured for optically and X-ray selected clusters by stacking data from the Planck satellite. Related systematics can also affect inferences from extreme clusters detected at high redshift. Our results illustrate that cosmological conclusions from galaxy cluster surveys depend critically on proper modelling, not only of the relevant physics, but also of the full distribution of the observables and of the selection biases induced by cluster identification procedures.

We present a dynamical model of supernova feedback which follows the evolution of pressurised bubbles driven by supernovae in a multi-phase interstellar medium (ISM). The bubbles are followed until the point of break-out into the halo, starting from an initial adiabatic phase to a radiative phase. We show that a key property which sets the fate of bubbles in the ISM is the gas surface density, through the work done by the expansion of bubbles and its role in setting the gas scaleheight. The multi-phase description of the ISM is essential, and neglecting it leads to order of magnitude differences in the predicted outflow rates. We compare our predicted mass loading and outflow velocities to observations of local and high-redshift galaxies and find good agreement over a wide range of stellar masses and velocities. With the aim of analysing the dependence of the mass loading of the outflow, beta (i.e. the ratio between the outflow and star formation rates), on galaxy properties, we embed our model in the galaxy formation simulation, GALFORM, set in the LCDM framework. We find that a dependence of beta solely on the circular velocity, as is widely assumed in the literature, is actually a poor description of the outflow rate, as large variations with redshift and galaxy properties are obtained. Moreover, we find that below a circular velocity of 80km/s the mass loading saturates. A more fundamental relation is that between beta and the gas scaleheight of the disk, hg, and the gas fraction, fgas, as beta hg^(1.1) fgas^(0.4), or the gas surface density, \Sigma_g, and the gas fraction, as beta \Sigma_g^(-0.6) fgas^(0.8). We find that using the new mass loading model leads to a shallower faint-end slope in the predicted optical and near-IR galaxy luminosity functions.

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.

Quasars represent the brightest Active Galactic Nuclei (AGN) in the Universe and are thought to indicate the location of prodigiously growing Black Holes (BHs), with luminosities as high as 10^48 erg/sec. It is often expected though that such an extremely energetic process will take place in the most massive bound structures in the dark matter (DM) distribution. We show that in contrast to this expectation, in a galaxy formation model which includes AGN feedback, quasars are predicted to live in average DM halo environments with typical masses of a few times 10^12 Msun. This fundamental prediction arises from the fact that quasar activity (i.e., BH accretion with luminosity greater than 10^46 erg/sec) is inhibited in DM haloes where AGN feedback operates. The galaxy hosts of quasars in our simulations are identified with over massive (in gas and stars) spheroidal galaxies, in which BH accretion is triggered via a galaxy merger or secular processes. We further show that the z=0 descendants of high redshift (z~6) QSOs span a wide range of morphologies, galaxy and halo masses. The z~6 BHs typically grow only by a modest factor by the present day. Remarkably, high redshift QSOs never inhabit the largest DM haloes at that time and their descendants are very seldom found in the most massive haloes at z=0. We also show that observationally it is very likely to find an enhancement in the abundance of galaxies around quasars at z~5. However, these enhancements are considerably weaker compared to the overdensities expected at the extreme peaks of the DM distribution. Thus, it is very unlikely that a quasar detected in the $z\gtrsim5$ Universe pinpoints the location of the progenitors of superclusters in the local Universe.

We present a clustering analysis of X-ray selected AGN by compiling X-ray samples from the literature and re-estimating the dark matter (DM) halo masses of AGN in a uniform manner. We find that moderate luminosity AGN (Lx(2-10 keV)=10^42-10^44 erg/sec) in the z=0-1.3 Universe are typically found in DM haloes with masses of ~10^13 Msun. We then compare our findings to the theoretical predictions of the coupled galaxy and black hole formation model GALFORM. We find good agreement when our calculation includes the hot-halo mode of accretion onto the central black hole. This type of accretion, which is additional to the common cold accretion during disk instabilities and galaxy mergers, is tightly coupled to the AGN feedback in the model. The hot-halo mode becomes prominent in DM haloes with masses greater than ~10^12.5 Msun, where AGN feedback typically operates, giving rise to a distinct class of moderate luminosity AGN that inhabit rich clusters and superclusters. Cold gas fuelling of the black hole cannot produce the observationally inferred DM halo masses of X-ray AGN. Switching off AGN feedback in the model results in a large population of luminous quasars (Lx(2-10 keV) > 10^44 erg/sec) in DM haloes with masses up to ~10^14 Msun, which is inconsistent with the observed clustering of quasars. The abundance of hot-halo AGN decreases significantly in the z~3-4 universe. At such high redshifts, the cold accretion mode is solely responsible for shaping the environment of moderate luminosity AGN. Our analysis supports two accretion modes (cold and hot) for the fuelling of supermassive black holes and strongly underlines the importance of AGN feedback in cosmological models both of galaxy formation and black hole growth.

(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.

Accurate galaxy stellar masses are crucial to better understand the physical mechanisms driving the galaxy formation process. We use synthetic star formation and metal enrichment histories predicted by the {\sc galform} galaxy formation model to investigate the precision with which various colours $(m_{a}-m_{b})$ can alone be used as diagnostics of the stellar mass-to-light ratio. As an example, we find that, at $z=0$, the {\em intrinsic} (B$_{f435w}-$V$_{f606w}$) colour can be used to determine the intrinsic rest-frame $V$-band stellar mass-to-light ratio ($\log_{10}\Gamma_{V}=\log_{10}[(M/M_{\odot})/(L_{V}/L_{V\odot})]$) with a precision of $\sigma_{lg\Gamma}\simeq 0.06$ when the initial mass function and redshift are known beforehand. While the presence of dust, assuming a universal attenuation curve, can have a systematic effect on the inferred mass-to-light ratio using a single-colour relation, this is typically small as it is often possible to choose a colour for which the dust reddening vector is approximately aligned with the $(m_{a}-m_{b})-\log_{10}\Gamma_{V}$ relation. The precision with which the stellar mass-to-light ratio can be recovered using a single colour diagnostic rivals implementations of SED fitting using more information but in which simple parameterisations of the star formation and metal enrichment histories are assumed. To facilitate the wide use of these relations, we provide the optimal observer frame colour to estimate the stellar mass-to-light ratio, along with the associated parameters, as a function of redshift ($0<z<1.5$) for two sets of commonly used filters.

Large surveys of the local Universe have shown that galaxies with different intrinsic properties, such as colour, luminosity and morphological type display a range of clustering amplitudes. Galaxies are therefore not faithful tracers of the underlying matter distribution. This modulation of galaxy clustering, called bias, contains information about the physics behind galaxy formation. It is also a systematic to be overcome before the large-scale structure of the Universe can be used as a cosmological probe. Two types of approaches have been developed to model the clustering of galaxies. The first class is empirical and filters or weights the distribution of dark matter to reproduce the measured clustering. In the second approach an attempt is made to model the physics which governs fate of baryons in order to predict the number of galaxies in dark matter haloes. I will review the development of both approaches and summarize what we have learnt about galaxy bias.

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 study the cosmology of Galileon modified gravity models in the linear perturbation regime. We derive the fully covariant and gauge invariant perturbed field equations using two different methods, which give consistent results, and solve them using a modified version of the {\tt CAMB} code. We find that, in addition to modifying the background expansion history and therefore shifting the positions of the acoustic peaks in the cosmic microwave background (CMB) power spectrum, the Galileon field can cluster strongly from early times, and causes the Weyl gravitational potential to grow, rather than decay, at late times. This leaves clear signatures in the low-$l$ CMB power spectrum through the modified integrated Sachs-Wolfe effect, strongly enhances the linear growth of matter density perturbations and makes distinctive predictions for other cosmological signals such as weak lensing and the power spectrum of density fluctuations. The quasi-static approximation is shown to work quite well from small to the near-horizon scales. We demonstrate that Galileon models display a rich phenomenology due to the large parameter space and the sensitive dependence of the model predictions on the Galileon parameters. Our results show that some Galileon models are already ruled out by present data and that future higher significance galaxy clustering, ISW and lensing measurements will place strong constraints on Galileon gravity.

We use the Millennium Simulation database to compare how different versions of the Durham and Munich semi-analytical galaxy formation models populate dark matter haloes with galaxies. The models follow the same physical processes but differ in how these are implemented. All of the models we consider use the Millennium N-body Simulation; however, the Durham and Munich groups use independent algorithms to construct halo merger histories. We compare the predicted halo occupation distributions (HODs) and correlation functions for galaxy samples defined by stellar mass, cold gas mass and star formation rate. The predictions for the HOD are remarkably similar for samples ranked by stellar mass. The predicted bias averaged over 5-25 Mpc/h is consistent between models to within 10%. At small pair separations there is a difference in the predicted clustering. This arises because the Durham models allow some satellite galaxies to merge with the central galaxy in a halo when they are still associated with resolved subhaloes. The agreement between the models is less good for samples defined by cold gas mass or star formation rate, with the small scale clustering differing by an order of magnitude, reflecting the uncertainty in the modelling of star formation. The model predictions in these cases are qualitatively similar, with a markedly shallower slope for the correlation function than is found for stellar mass selected samples and with the HOD displaying an asymmetric peak for central galaxies. We provide illustrative parametric fits to the HODs predicted by the models. Our results reveal the current limitations on how well we can predict galaxy bias in a fixed cosmology, which has implications for the interpretation of constraints on the physics of galaxy formation from galaxy clustering measurements and the ability of future galaxy surveys to measure dark energy.

The clustering of galaxies observed in future redshift surveys will provide a wealth of cosmological information. Matching the signal at different redshifts constrains the dark energy driving the acceleration of the expansion of the Universe. In tandem with these geometrical constraints, redshift-space distortions (RSD) depend on the build up of large-scale structure. As pointed out by many authors measurements of these effects are intrinsically coupled. We investigate this link, and argue that it strongly depends on the cosmological assumptions adopted when analysing data. Using representative assumptions for the parameters of the Euclid survey in order to provide a baseline future experiment, we show how the derived constraints change due to different model assumptions. We argue that even the assumption of a Friedman-Robertson-Walker (FRW) space-time is sufficient to reduce the importance of the coupling to a significant degree. Taking this idea further, we consider how the data would actually be analysed and argue that we should not expect to be able to simultaneously constrain multiple deviations from the standard $\Lambda$CDM model. We therefore consider different possible ways in which the Universe could deviate from the $\Lambda$CDM model, and show how the coupling between geometrical constraints and structure growth affects the measurement of such deviations.

Using GALFORM, a semi-analytical model of galaxy formation in the Lambda cold dark matter cosmology, we study the rest-frame ultraviolet (UV) colours of Lyman-break galaxies (LBGs) in the redshift range 2.5 < z < 10. As the impact of dust on UV luminosity can be dramatic, our model includes a self-consistent computation of dust attenuation based on a radiative transfer model. We find that intrinsically brighter galaxies suffer stronger dust attenuation than fainter ones, though the relation has a large scatter. The model predicts galaxies with UV colours consistent with the colour selection regions designed to select LBGs in observational surveys. We find that the drop-out technique that selects LBGs based on two rest-frame UV colours is robust and effective, selecting more than 70 per cent of UV bright galaxies at a given redshift. We investigate the impact on the predicted UV colours of varying selected model parameters. We find that the UV colours are most sensitive to the modelling of dust attenuation and in particular, to the extinction curve used in the radiative transfer calculation. If we assume a Milky Way dust extinction curve, the predicted UV continuum slopes are, in general, bluer than observed. However, we find that the opposite is true when using the Small Magellanic Cloud dust extinction curve. This demonstrates the strong dependence of UV colours on dust properties and highlights the inadequacy of using the UV continuum slope as a tracer of dust attenuation without any further knowledge of the galaxy inclination or dust characteristics in high redshift galaxies.

We investigate the impact of feedback - from supernovae (SNe), active galactic nuclei (AGN) and a photo-ionizing background at high redshifts - on the neutral atomic hydrogen (HI) mass function, the $b_{\rm J}$ band luminosity function, and the spatial clustering of these galaxies at $z$=0. We use a version of the semi-analytical galaxy formation model GALFORM that calculates self-consistently the amount of HI in a galaxy as a function of cosmic time and links its star formation rate to its mass of molecular hydrogen (H$_2$). We find that a systematic increase or decrease in the strength of SNe feedback leads to a systematic decrease or increase in the amplitudes of the luminosity and HI mass functions, but has little influence on their overall shapes. Varying the strength of AGN feedback influences only the numbers of the brightest or most HI massive galaxies, while the impact of varying the strength of photo-ionization feedback is restricted to changing the numbers of the faintest or least HI massive galaxies.Our results suggest that the HI mass function is a more sensitive probe of the consequences of cosmological reionization for galaxy formation than the luminosity function. We find that increasing the strength of any of the modes of feedback acts to weaken the clustering strength of galaxies, regardless of their HI-richness. In contrast, weaker AGN feedback has little effect on the clustering strength whereas weaker SNe feedback increases the clustering strength of HI-poor galaxies more strongly than HI-rich galaxies. These results indicate that forthcoming HI surveys on next generation radio telescopes such as the Square Kilometre Array and its pathfinders will be exploited most fruitfully as part of multiwavelength survey campaigns.

The latest observations of molecular gas and the atomic hydrogen content of local and high-redshift galaxies, coupled with how these correlate with star formation activity, have revolutionized our ideas about how to model star formation in a galactic context. A successful theory of galaxy formation has to explain some key facts: (i) high-redshift galaxies have higher molecular gas fractions and star formation rates than local galaxies, (ii) scaling relations show that the atomic-to-stellar mass ratio decreases with stellar mass in the local Universe, and (iii) the global abundance of atomic hydrogen evolves very weakly with time. We review how modern cosmological simulations of galaxy formation attempt to put these pieces together and highlight how approaches simultaneously solving dark matter and gas physics, and approaches first solving the dark matter N-body problem and then dealing with gas physics using semi-analytic models, differ and complement each other. We review the observable predictions, what we think we have learned so far and what still needs to be done in the simulations to allow robust testing by the new observations expected from telescopes such as ALMA, PdBI, LMT, JVLA, ASKAP, MeerKAT, SKA.

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 matter and velocity divergence power spectra in a f(R) gravity theory and their time evolution measured from several large-volume N-body simulations with varying box sizes and resolution. We find that accurate prediction of the matter power spectrum in f(R) gravity places stronger requirements on the simulation than is the case with LCDM, because of the nonlinear nature of the fifth force. Linear perturbation theory is shown to be a poor approximation for the f(R) models, except when the chameleon effect is very weak. We show that the relative differences from the fiducial LCDM model are much more pronounced in the nonlinear tail of the velocity divergence power spectrum than in the matter power spectrum, which suggests that future surveys which target the collection of peculiar velocity data will open new opportunities to constrain modified gravity theories. A close investigation of the time evolution of the power spectra shows that there is a pattern in the evolution history, which can be explained by the properties of the chameleon-type fifth force in f(R) gravity. Varying the model parameter |f_R0|, which quantifies the strength of the departure from standard gravity, mainly varies the epoch marking the onset of the fifth force, as a result of which the different f(R) models are in different stages of the same evolutionary path at any given time

The rest-frame intrinsic UV luminosity is often used as an indicator of the instantaneous star formation rate (SFR) in a galaxy. While it is in general a robust indicator of the ongoing star formation activity, the precise value of the calibration relating the UV luminosity to the SFR ($B_{\nu}$), is sensitive to various physical properties, such as the recent star formation and metal enrichment histories, along with the choice of stellar initial mass function. The distribution of these properties for the star-forming galaxy population then suggests that the adoption of a single calibration is not appropriate unless properly qualified with the uncertainties on the calibration. We investigate, with the aid of the {\sc galform} semi-analytic model of galaxy formation, the distribution of UV-SFR calibrations obtained using realistic star formation and metal enrichment histories. At $z=0$, we find that when the initial mass function is fixed (to the Kennicutt IMF), the median calibration is $B_{\rm fuv}=0.9$ where ${\rm SFR}/[{\rm M_{\odot}\,yr^{-1}}]=B_{\nu}\times 10^{-28}\times L_{\nu}/[{\rm ergs\,s^{-1}\,Hz^{-1}}]$. However, the width of the distribution $B_{\rm fuv}$ suggests that for a single object there is around a 20% {\em intrinsic} uncertainty (at $z=0$, rising to $\simeq 30%$ at $z=6$) on the star formation rate inferred from the FUV luminosity without additional constraints on the star formation history or metallicity. We also find that the median value of the calibration $B_{\rm fuv}$ is correlated with the star formation rate and redshift (at $z>3$) raising implications for the correct determination of the star formation rate from the UV.

We present predictions for the clustering of Extremely Red Objects (EROs) in a \Lambda cold dark matter Universe, using a semi-analytical galaxy formation model in combination with a cosmological N-body simulation. EROs are red, massive galaxies observed at 0.7< z < 3, and their numbers and properties have posed a challenge to hierarchical galaxy formation models. We analyse the halo occupation distribution and two-point correlation function of EROs, exploring how these quantities change with apparent magnitude, colour cut and redshift. Our model predicts a halo occupation distribution that is significantly different from that typically assumed. This is due to the inclusion of AGN feedback, which changes the slope and scatter of the luminosity-host halo mass relation above the mass where AGN feedback first becomes important. We predict that, on average, dark matter haloes with masses above 10^{13}h^{-1}Msun host at least one ERO at 1.5 < z < 2.5. Taking into account sample variance in observational estimates, the predicted angular clustering for EROs with either (R-K)>5 or (i-K)>4.5 is in reasonable agreement with current observations.

We combine the galaxy formation model GALFORM with the Photon Dominated Region code UCL_PDR to study the emission from the rotational transitions of 12CO (CO) in galaxies from z=0 to z=6 in the Lambda CDM framework. GALFORM is used to predict the molecular (H2) and atomic hydrogen (HI) gas contents of galaxies using the pressure-based empirical star formation relation of Blitz & Rosolowsky. From the predicted H2 mass and the conditions in the interstellar medium, we estimate the CO emission in the rotational transitions 1-0 to 10-9 by applying the UCL_PDR model to each galaxy. We find that deviations from the Milky-Way CO-H2 conversion factor come mainly from variations in metallicity, and in the average gas and star formation rate surface densities. In the local universe, the model predicts a CO(1-0) luminosity function (LF), CO-to-total infrared (IR) luminosity ratios for multiple CO lines and a CO spectral line energy distribution (SLED) which are in good agreement with observations of luminous and ultra-luminous IR galaxies. At high redshifts, the predicted CO SLED of the brightest IR galaxies reproduces the shape and normalization of the observed CO SLED. The model predicts little evolution in the CO-to-IR luminosity ratio for different CO transitions, in good agreement with observations up to z~5. We use this new hybrid model to explore the potential of using colour selected samples of high-redshift star-forming galaxies to characterise the evolution of the cold gas mass in galaxies through observations with the Atacama Large Millimeter Array.

Using a sample of BzK-selected galaxies at z~2 identified from the CFHT/WIRCAM near-infrared survey of GOODS-North, we discuss the relation between star formation rate (SFR), specific star formation rate (SSFR), and stellar mass (M_{*}), and the clustering of galaxies as a function of these parameters. For star-forming galaxies (sBzKs), the UV-based SFR, corrected for extinction, scales with the stellar mass as SFR ~ M_{*}^{alpha} with alpha = 0.74+/-0.20 down to M_{*} ~ 10^{9} M_{solar}, indicating a weak dependence on the stellar mass of the SSFR. We also measure the angular correlation function and hence infer the correlation length for sBzK galaxies as a function of M_{*}, SFR, and SSFR, as well as K-band apparent magnitude. We show that passive galaxies (pBzKs) are more strongly clustered than sBzK galaxies at a given stellar mass, mirroring the color-density relation seen at lower redshifts. We also find that the correlation length of sBzK galaxies ranges from 4 to 20 h^{-1}Mpc, being a strong function of M_{K}, M_{*}, and SFR. On the other hand, the clustering dependence on SSFR changes abruptly at 2x10^{-9} yr^{-1}, which is the typical value for "main sequence" star-forming galaxies at z~2. We show that the correlation length reaches a minimum at this characteristic value, and is larger for galaxies with both smaller and larger SSFRs; a dichotomy that is only marginally implied from the predictions of the semi-analytical models. Our results suggest that there are two types of environmental effects at work at z~2. Stronger clustering for relatively quiescent galaxies implies that the environment has started to play a role in quenching star formation. At the same time, stronger clustering for galaxies with elevated SSFRs ("starbursts") might be attributed to an increased efficiency for galaxy interactions and mergers in dense environments.

We use large volume N-body simulations to predict the clustering of dark matter in redshift space in f(R) modified gravity cosmologies. This is the first time that the nonlinear matter and velocity fields have been resolved to such a high level of accuracy over a broad range of scales in this class of models. We find significant deviations from the clustering signal in standard gravity, with an enhanced boost in power on large scales and stronger damping on small scales in the f(R) models compared to GR at redshifts z<1. We measure the velocity divergence (P_\theta \theta) and matter (P_\delta \delta) power spectra and find a large deviation in the ratios \sqrt{P_\theta \theta/P_\delta \delta} and P_\delta \theta/P_\delta\delta, between the f(R) models and GR for 0.03<k/(h/Mpc)<0.5. In linear theory these ratios equal the growth rate of structure on large scales. Our results show that the simulated ratios agree with the growth rate for each cosmology (which is scale dependent in the case of modified gravity) only for extremely large scales, k<0.06h/Mpc at z=0. The velocity power spectrum is substantially different in the f(R) models compared to GR, suggesting that this observable is a sensitive probe of modified gravity. We demonstrate how to extract the matter and velocity power spectra from the 2D redshift space power spectrum, P(k,\mu), and can recover the nonlinear matter power spectrum to within a few percent for k<0.1h/Mpc. However, the model fails to describe the shape of the 2D power spectrum demonstrating that an improved model is necessary in order to reconstruct the velocity power spectrum accurately. The same model can match the monopole moment to within 3% for GR and 10% for the f(R) cosmology at k<0.2 h/Mpc at z=1. Our results suggest that the extraction of the velocity power spectrum from future galaxy surveys is a promising method to constrain deviations from GR.

(abridged) We present a clustering analysis of 370 high-confidence H-alpha emitters (HAEs) at z=2.23. The HAEs are detected in the Hi-Z Emission Line Survey (HiZELS), a large-area blank field 2.121um narrowband survey using the United Kingdom Infrared Telescope (UKIRT) Wide Field Camera (WFCAM). Averaging the two-point correlation function of HAEs in two ~1 degree scale fields (United Kingdom Infrared Deep Sky Survey/Ultra Deep Survey [UDS] and Cosmological Evolution Survey [COSMOS] fields) we find a clustering amplitude equivalent to a correlation length of r_0=3.7+/-0.3 Mpc/h for galaxies with star formation rates of >7 M_sun/yr. The data are well-fitted by the expected correlation function of Cold Dark Matter, scaled by a bias factor: omega_HAE=b^2 omega_DM where b=2.4^{+0.1}_{-0.2}. The corresponding 'characteristic' mass for the halos hosting HAEs is log(M_h/[M_sun/h])=11.7+/-0.1. Comparing to the latest semi-analytic GALFORM predictions for the evolution of HAEs in a LCDM cosmology, we find broad agreement with the observations, with GALFORM predicting a HAE correlation length of ~4 Mpc/h. Motivated by this agreement, we exploit the simulations to construct a parametric model of the halo occupation distribution (HOD) of HAEs, and use this to fit the observed clustering. Our best-fitting HOD can adequately reproduce the observed angular clustering of HAEs, yielding an effective halo mass and bias in agreement with that derived from the scaled omega_DM fit, but with the relatively small sample size the current data provide a poor constraint on the HOD. Our results support the broad picture that 'typical' (~L*) star-forming galaxies have been hosted by dark matter haloes with M_h<10^12 M_sun/h since z~2, but with a broad occupation distribution and clustering that is likely to be a strong function of luminosity.

The observed ultraviolet continuum (UVC) slope is potentially a powerful diagnostic of dust obscuration in star forming galaxies. However, the intrinsic slope is also sensitive to the form of the stellar initial mass function (IMF) and to the recent star formation and metal enrichment histories of a galaxy. Using the galform semi-analytical model of galaxy formation, we investigate the intrinsic distribution of UVC slopes. For star-forming galaxies, we find that the intrinsic distribution of UVC slopes at z=0, parameterised by the power law index beta, has a standard deviation of sigma_beta=0.30. This suggests an uncertainty on the inferred UV attenuation of A_fuv=0.7$ (assuming a Calzetti attenuation curve) for an individual object, even with perfect photometry. Furthermore, we find that the intrinsic UVC slope correlates with star formation rate, intrinsic UV luminosity, stellar mass and redshift. These correlations have implications for the interpretation of trends in the observed UVC slope with these quantities irrespective of the sample size or quality of the photometry. Our results suggest that in some cases the attenuation by dust has been incorrectly estimated.

The distribution of angles subtended between pairs of galaxies and the line of sight,which is uniform in real space, is distorted by their peculiar motions, and has been proposed as a probe of cosmic expansion. We test this idea using N-body simulations of structure formation in a cold dark matter universe with a cosmological constant and in two variant cosmologies with different dark energy models. We find that the distortion of the distribution of angles is sensitive to the nature of dark energy. However, for the first time, our simulations also reveal dependences of the normalization of the distribution on both redshift and cosmology that have been neglected in previous work. This introduces systematics that severely limit the usefulness of the original method. Guided by our simulations, we devise a new, improved test of the nature of dark energy. We demonstrate that this test does not require prior knowledge of the background cosmology and that it can even distinguish between models that have the same baryonic acoustic oscillations and dark matter halo mass functions. Our technique could be applied to the completed BOSS galaxy redshift survey to constrain the expansion history of the Universe to better than 2%. The method will also produce different signals for dark energy and modified gravity cosmologies even when they have identical expansion histories, through the different peculiar velocities induced in these cases.

We study the properties of Ly-alpha emitters in a cosmological framework by computing the escape of Ly-alpha photons through galactic outflows. We combine the GALFORM semi-analytical model of galaxy formation with a Monte Carlo Ly-alpha radiative transfer code. The properties of Ly-alpha emitters at 0<z<7 are predicted using two outflow geometries: a Shell of neutral gas and a Wind ejecting material, both expanding at constant velocity. We characterise the differences in the Ly-alpha line profiles predicted by the two outflow geometries in terms of their width, asymmetry and shift from the line centre for a set of outflows with different hydrogen column densities, expansion velocities and metallicities. In general, the Ly-alpha line profile of the Shell geometry is broader and more asymmetric, and the Ly-alpha escape fraction is lower than with the Wind geometry for the same set of parameters. In order to implement the outflow geometries in the semi-analytical model GALFORM, a number of free parameters in the outflow model are set by matching the luminosity function of Ly-alpha emitters over the whole observed redshift range. The models are consistent with the observationally inferred Ly-alpha escape fractions, equivalent width distributions and with the shape of the Ly-alpha line from composite spectra. Interestingly, our predicted UV luminosity function of Ly-alpha emitters and the fraction of Ly-alpha emitters in Lyman-break galaxy samples at high redshift are in partial agreement with observations. Attenuation of the Ly-alpha line by the presence of a neutral intergalactic medium at high redshift could be responsible for this disagreement. We predict that Ly-alpha emitters constitute a subset of the galaxy population with lower metallicities, lower instantaneous star formation rates and larger sizes than the overall population at the same UV luminosity.

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

We study the evolution of the cold gas content of galaxies by splitting the interstellar medium into its atomic and molecular hydrogen components, using the galaxy formation model GALFORM in the LCDM framework. We calculate the molecular-to-atomic hydrogen mass ratio, H2/HI, in each galaxy using two different approaches; the pressure-based empirical relation of Blitz & Rosolowsky and the theoretical model of Krumholz, McKeee & Tumlinson, and apply them to consistently calculate the star formation rates of galaxies. We find that the model based on the Blitz & Rosolowsky law predicts an HI mass function, CO(1-0) luminosity function, correlations between the H2/HI ratio and stellar and cold gas mass, and infrared-CO luminosity relation in good agreement with local and high redshift observations. The HI mass function evolves weakly with redshift, with the number density of high mass galaxies decreasing with increasing redshift. In the case of the H2 mass function, the number density of massive galaxies increases strongly from z=0 to z=2, followed by weak evolution up to z=4. We also find that the H2/HI ratio of galaxies is strongly dependent on stellar and cold gas mass, and also on redshift. The slopes of the correlations between H2/HI and stellar and cold gas mass hardly evolve, but the normalisation increases by up to two orders of magnitude from z=0-8. The strong evolution in the H2 mass function and the H2/HI ratio is primarily due to the evolution in the sizes of galaxies and secondarily, in the gas fractions. The predicted cosmic density evolution of HI agrees with the observed evolution inferred from DLAs, and is dominated by low/intermediate mass halos. We find that previous theoretical studies have largely overestimated the redshift evolution of the global H2/HI ratio due to limited resolution. We predict a maximum of rho_H2/rho_HI~1.2 at z~3.5.

For galaxy clustering to provide robust constraints on cosmological parameters and galaxy formation models, it is essential to make reliable estimates of the errors on clustering measurements. We present a new technique, based on a spatial Jackknife (JK) resampling, which provides an objective way to estimate errors on clustering statistics. Our approach allows us to set the appropriate size for the Jackknife subsamples. The method also provides a means to assess the impact of individual regions on the measured clustering, and thereby to establish whether or not a given galaxy catalogue is dominated by one or several large structures, preventing it to be considered as a "fair sample". We apply this methodology to the two- and three-point correlation functions measured from a volume limited sample of M* galaxies drawn from data release seven of the Sloan Digital Sky Survey (SDSS). The frequency of jackknife subsample outliers in the data is shown to be consistent with that seen in large N-body simulations of clustering in the cosmological constant plus cold dark matter cosmology. We also present a comparison of the three-point correlation function in SDSS and 2dFGRS using this approach and find consistent measurements between the two samples.

BigBOSS is a Stage IV ground-based dark energy experiment to study baryon acoustic oscillations (BAO) and the growth of structure with a wide-area galaxy and quasar redshift survey over 14,000 square degrees. It has been conditionally accepted by NOAO in response to a call for major new instrumentation and a high-impact science program for the 4-m Mayall telescope at Kitt Peak. The BigBOSS instrument is a robotically-actuated, fiber-fed spectrograph capable of taking 5000 simultaneous spectra over a wavelength range from 340 nm to 1060 nm, with a resolution R = 3000-4800. Using data from imaging surveys that are already underway, spectroscopic targets are selected that trace the underlying dark matter distribution. In particular, targets include luminous red galaxies (LRGs) up to z = 1.0, extending the BOSS LRG survey in both redshift and survey area. To probe the universe out to even higher redshift, BigBOSS will target bright [OII] emission line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are obtained to measure the BAO feature, trace the matter power spectrum at smaller scales, and detect redshift space distortions. BigBOSS will provide additional constraints on early dark energy and on the curvature of the universe by measuring the Ly-alpha forest in the spectra of over 600,000 2.2 < z < 3.5 quasars. BigBOSS galaxy BAO measurements combined with an analysis of the broadband power, including the Ly-alpha forest in BigBOSS quasar spectra, achieves a FOM of 395 with Planck plus Stage III priors. This FOM is based on conservative assumptions for the analysis of broad band power (kmax = 0.15), and could grow to over 600 if current work allows us to push the analysis to higher wave numbers (kmax = 0.3). BigBOSS will also place constraints on theories of modified gravity and inflation, and will measure the sum of neutrino masses to 0.024 eV accuracy.

We investigate the consequences of applying different star formation laws in the galaxy formation model GALFORM. Three broad star formation laws are implemented: the empirical relations of Kennicutt and Schmidt and Blitz & Rosolowsky and the theoretical model of Krumholz, McKee & Tumlinson. These laws have no free parameters once calibrated against observations of the star formation rate (SFR) and gas surface density in nearby galaxies. We start from published models, and investigate which observables are sensitive to a change in the star formation law, without altering any other model parameters. We show that changing the star formation law (i) does not significantly affect either the star formation history of the universe or the galaxy luminosity functions in the optical and near-IR, due to an effective balance between the quiescent and burst star formation modes; (ii) greatly affects the cold gas contents of galaxies; (iii) changes the location of galaxies in the SFR versus stellar mass plane, so that a second sequence of "passive" galaxies arises, in addition to the known "active" sequence. We show that this plane can be used to discriminate between the star formation laws.

We predict the formation histories, properties and descendants of Lyman-break galaxies (LBGs) in the Lambda-CDM cosmology. In our model, which incorporates a top-heavy IMF in starbursts, we find that most LBGs are starbursts triggered by minor mergers of galaxies. The duration of the LBG phase is predicted to be quite short, ~20-60 Myr. We investigate the distributions of stellar and halo masses and morphologies for bright (L_UV > L*_UV) and faint (L_UV > 0.1 L*_UV) LBGs at z=3, z=6 and z=10 (where we classify LBGs according to their rest-frame UV luminosities relative the observed characteristic luminosity L*_UV at z \approx 3). Bright LBGs at z=3 are predicted to have median stellar masses ~ 1x10^9 Msun/h and host halo masses ~ 3x10^{11} Msun/h, and to be typically mildly disk-dominated in stellar mass. On the other hand, faint LBGs at z=10 are predicted to have median stellar masses of only ~ 1x10^7 Msun/h and host halo masses 2x10^{10} Msun/h, and to be generally bulge-dominated. Bright LBGs at z=3 evolve to present-day galaxies with median stellar mass ~ 5x^{10} Msun/h (comparable to the Milky Way), consisting of roughly equal numbers of disk- and bulge-dominated systems, and hosted by halos with median mass ~2x10^{13} Msun/h (corresponding to medium-size galaxy groups). The model predicts that 40% of Milky Way mass galaxies at the present-day have a bright LBG progenitor in the redshift range 3<z<4, while 95% have a faint LBG progenitor in the same redshift range, and 7% have a faint LBG progenitor at 10<z<11. With our multiwavelength model, we also investigate the overlap between the LBG population and that of submillimetre selected galaxies (SMGs); at z=3, only ~1% of bright LBGs are also predicted to also be bright SMGs (with an 850 mum flux in excess of 5 mJy).

We introduce an objective method to assess the probability of finding extreme events in the distribution of cold dark matter such as voids, overdensities or very high mass haloes. Our approach uses an ensemble of N-body simulations of the hierar- chical clustering of dark matter to find extreme structures. The frequency of extreme events, in our case the cell or smoothing volume with the highest count of cluster-mass dark matter haloes, is well described by a Gumbel distribution. This distribution can then be used to forecast the probability of finding even more extreme events, which would otherwise require a much larger ensemble of simulations to quantify. We use our technique to assess the chance of finding concentrations of massive clusters or super- clusters, like the two found in the two-degree field galaxy redshift survey (2dFGRS), using a counts-in-cells analysis. The Gumbel distribution gives an excellent descrip- tion of the distribution of extreme cell counts across two large ensembles of simulations covering different cosmologies, and also when measuring the clustering in both real and redshift space. We find examples of structures like those found in the 2dFGRS in the simulations. The chance of finding such structures in a volume equal to that of the 2dFGRS is around 2%.

Future galaxy surveys hope to distinguish between the dark energy and modified gravity scenarios for the accelerating expansion of the Universe using the distortion of clustering in redshift space. The aim is to model the form and size of the distortion to infer the rate at which large scale structure grows. We test this hypothesis and assess the performance of current theoretical models for the redshift space distortion using large volume N-body simulations of the gravitational instability process. We simulate competing cosmological models which have identical expansion histories - one is a quintessence dark energy model with a scalar field and the other is a modified gravity model with a time varying gravitational constant - and demonstrate that they do indeed produce different redshift space distortions. This is the first time this approach has been verified using a technique that can follow the growth of structure at the required level of accuracy. Our comparisons show that theoretical models for the redshift space distortion based on linear perturbation theory give a surprisingly poor description of the simulation results. Furthermore, the application of such models can give rise to catastrophic systematic errors leading to incorrect interpretation of the observations. We show that an improved model is able to extract the correct growth rate. Further enhancements to theoretical models of redshift space distortions, calibrated against simulations, are needed to fully exploit the forthcoming high precision clustering measurements.

We study the role of submillimetre galaxies (SMGs) in the galaxy formation process in the Lambda Cold Dark Matter cosmology. We use the Baugh et al. (2005) semi-analytical model, which matches the observed SMG number counts and redshift distribution by assuming a top-heavy initial mass function (IMF) in bursts triggered by galaxy mergers. We build galaxy merger trees and follow the evolution and properties of SMGs and their descendants. Our primary sample of model SMGs consists of galaxies which had 850 mu fluxes brighter than 5 mJy at some redshift z>1. Our model predicts that the present-day descendants of such SMGs cover a wide range of stellar masses ~ 10^{10} - 10^{12} Msun/h, with a median ~ 10^{11} Msun/h, and that more than 70% of these descendants are bulge-dominated. More than 50% of present day galaxies with stellar masses larger than 7 x 10^{11} Msun/h are predicted to be descendants of such SMGs. We find that although SMGs make an important contribution to the total star formation rate at z~2, the final stellar mass produced in the submillimetre phase contributes only 0.2% of the total present-day stellar mass, and 2% of the stellar mass of SMG descendants, in stark contrast to the popular picture in which the SMG phase marks the production of the bulk of the mass of present day massive ellipticals.

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 make a detailed investigation of the properties of Lyman-break galaxies (LBGs) in the LambdaCDM model. We present predictions for two published variants of the GALFORM semi-analytical model: the Baugh et al. (2005) model, which has star formation at high redshifts dominated by merger-driven starbursts with a top-heavy IMF, and the Bower et al. (2006) model, which has AGN feedback and a standard Solar neighbourhood IMF throughout. We show predictions for the evolution of the rest-frame far-UV luminosity function in the redshift range z=3-20, and compare with the observed luminosity functions of LBGs at z=3-10. We find that the Baugh et al. model is in excellent agreement with these observations, while the Bower et al. model predicts too many high-luminosity LBGs. Dust extinction, which is predicted self-consistently based on galaxy gas contents, metallicities and sizes, is found to have a large effect on LBG luminosities. We compare predictions for the size evolution of LBGs at different luminosities with observational data for 2<z<7, and find the Baugh et al. model to be in good agreement. We present predictions for stellar, halo and gas masses, star formation rates, circular velocities, bulge-to-disk ratios, gas and stellar metallicities and clustering bias, as functions of far-UV luminosity and redshift. We find broad consistency with current observational constraints. We then present predictions for the abundance and angular sizes of LBGs out to very high redshift (z<20), finding that planned deep surveys with JWST should detect objects out to z<15. The typical UV luminosities of galaxies are predicted to be very low at high redshifts, which has implications for detecting the galaxies responsible for reionizing the IGM; for example, at z=10, 50% of the ionizing photons are expected to be produced by galaxies fainter than M_AB(1500A)-5logh ~ -15.

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 combine the GALFORM semi-analytical model of galaxy formation, which predicts the star formation and merger histories of galaxies, the GRASIL spectro-photometric code, which calculates the spectral energy distributions (SEDs) of galaxies self-consistently including reprocessing of radiation by dust, and artificial neural networks (ANN), to investigate the clustering properties of galaxies selected by their emission at submillimetre wavelengths (SMGs). We use the Millennium Simulation to predict the spatial and angular distribution of SMGs. At redshift z = 2, we find that these galaxies are strongly clustered, with a comoving correlation length of r0 = 5.6 \pm 0.9 Mpc/h for galaxies with 850{\mu}m flux densities brighter than 5 mJy, in agreement with observations. We predict that at higher redshifts these galaxies trace denser and increasingly rarer regions of the universe. We present the predicted dependence of the clustering on luminosity, submillimetre colour, halo and total stellar masses. Interestingly, we predict tight relations between correlation length and halo and stellar masses, independent of sub-mm luminosity.

The spectral energy distribution of galaxies is a complex function of the star formation history and geometrical arrangement of stars and gas in galaxies. The computation of the radiative transfer of stellar radiation through the dust distribution is time-consuming. This aspect becomes unacceptable in particular when dealing with the predictions by semi-analytical galaxy formation models populating cosmological volumes, to be then compared with multi-wavelength surveys. Mainly for this aim, we have implemented an artificial neural network algorithm into the spectro-photometric and radiative transfer code GRASIL in order to compute the spectral energy distribution of galaxies in a short computing time. This allows to avoid the adoption of empirical templates that may have nothing to do with the mock galaxies output by models. The ANN has been implemented to compute the dust emission spectrum (the bottleneck of the computation), and separately for the star-forming molecular clouds and the diffuse dust (due to their different properties and dependencies). We have defined the input neurons effectively determining their emission, which means this implementation has a general applicability and is not linked to a particular galaxy formation model. We have trained the net for the disc and spherical geometries, and tested its performance to reproduce the SED of disc and starburst galaxies, as well as for a semi-analytical model for spheroidal galaxies. We have checked that for this model both the SEDs and the galaxy counts in the Herschel bands obtained with the ANN approximation are almost superimposed to the same quantities obtained with the full GRASIL. We conclude that this method appears robust and advantageous, and will present the application to a more complex SAM in another paper.

The anisotropy of clustering in redshift space provides a direct measure of the growth rate of large scale structure in the Universe. Future galaxy redshift surveys will make high precision measurements of these distortions, and will potentially allow us to distinguish between different scenarios for the accelerating expansion of the Universe. Accurate predictions are needed in order to distinguish between competing cosmological models. We study the distortions in the redshift space power spectrum in $\Lambda$CDM and quintessence dark energy models, using large volume N-body simulations, and test predictions for the form of the redshift space distortions. We find that the linear perturbation theory prediction by Kaiser (1987) is a poor fit to the measured distortions, even on surprisingly large scales $k \ge 0.05 h$Mpc$^{-1}$. An improved model for the redshift space power spectrum, including the non-linear velocity divergence power spectrum, is presented and agrees with the power spectra measured from the simulations up to $k \sim 0.2 h$Mpc$^{-1}$. We have found a density-velocity relation which is cosmology independent and which relates the non-linear velocity divergence spectrum to the non-linear matter power spectrum. We provide a formula which generates the non-linear velocity divergence $P(k)$ at any redshift, using only the non-linear matter power spectrum and the linear growth factor at the desired redshift. This formula is accurate to better than 5% on scales $k<0.2 h $Mpc$^{-1}$ for all the cosmological models discussed in this paper. Our results will extend the statistical power of future galaxy surveys.

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.

A space-based galaxy redshift survey would have enormous power in constraining dark energy and testing general relativity, provided that its parameters are suitably optimized. We study viable space-based galaxy redshift surveys, exploring the dependence of the Dark Energy Task Force (DETF) figure-of-merit (FoM) on redshift accuracy, redshift range, survey area, target selection, and forecast method. Fitting formulae are provided for convenience. We also consider the dependence on the information used: the full galaxy power spectrum P(k, P(k) marginalized over its shape, or just the Baryon Acoustic Oscillations (BAO). We find that the inclusion of growth rate information (extracted using redshift space distortion and galaxy clustering amplitude measurements) leads to a factor of ~ 3 improvement in the FoM, assuming general relativity is not modified. This inclusion partially compensates for the loss of information when only the BAO are used to give geometrical constraints, rather than using the full P(k) as a standard ruler. We find that a space-based galaxy redshift survey covering ~20,000 deg^2 over 0.5< z < 2 with \sigma_z/(1+z)<= 0.001 exploits a redshift range that is only easily accessible from space, extends to sufficiently low redshifts to allow both a vast 3-D map of the universe using a single tracer population, and overlaps with ground-based surveys to enable robust modeling of systematic effects. We argue that these parameters are close to their optimal values given current instrumental and practical constraints.

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 predictions for the abundance and nature of Extremely Red Objects (EROs) in the Lambda cold dark matter model. EROs are red, massive galaxies observed at z>= 1 and their numbers and properties pose a challenge to hierarchical galaxy formation models. We compare the predictions from two published models, one of which invokes a "superwind" to regulate star formation in massive haloes and the other which suppresses gas cooling in haloes through "radio-mode" AGN feedback. The superwind model underestimates the number counts of EROs by an order of magnitude, whereas the radio-mode AGN feedback model gives excellent agreement with the number counts and redshift distribution of EROs. In the AGN feedback model the ERO population is dominated by old, passively evolving galaxies, whereas observations favour an equal split between old galaxies and dusty starbursts. Also, the model predicts a more extended redshift distribution of passive galaxies than is observed. These comparisons suggest that star formation may be quenched too efficiently in this model.

We analyse whether hierarchical formation models based on Lambda cold dark matter cosmology can produce enough massive red galaxies to match observations. For this purpose, we compare with observations the predictions from two published models for the abundance and redshift distribution of Extremely Red Objects (EROs), which are red, massive galaxies observed at z >= 1. One of the models invokes a "superwind" to regulate star formation in massive haloes and the other suppresses cooling through "radio-mode" AGN feedback. The first one underestimates the number counts of EROs by an order of magnitude, whereas the radio-mode AGN feedback model gives excellent agreement with the number counts of EROs and redshift distribution of K-selected galaxies. This study highlights the need to consider AGN feedback in order to understand the formation and evolution of massive galaxies at z >= 1.

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.

We model the nonlinear growth of cosmic structure in different dark energy models, using large volume N-body simulations. We consider a range of quintessence models which feature both rapidly and slowly varying dark energy equations of state, and compare the growth of structure to that in a universe with a cosmological constant. The adoption of a quintessence model changes the expansion history of the universe, the form of the linear theory power spectrum and can alter key observables, such as the horizon scale and the distance to last scattering. The difference in structure formation can be explained to first order by the difference in growth factor at a given epoch; this scaling also accounts for the nonlinear growth at the 15% level. We find that quintessence models which feature late $(z<2)$, rapid transitions towards $w=-1$ in the equation of state, can have identical baryonic acoustic oscillation (BAO) peak positions to those in $\Lambda$CDM, despite being very different from $\Lambda$CDM both today and at high redshifts $(z \sim 1000)$. We find that a second class of models which feature non-negligible amounts of dark energy at early times cannot be distinguished from $\Lambda$CDM using measurements of the mass function or the BAO. These results highlight the need to accurately model quintessence dark energy in N-body simulations when testing cosmological probes of dynamical dark energy.

(Abridged) Accurately predicting how the cosmic abundance of neutral hydrogen evolves with redshift is a challenging problem facing modellers of galaxy formation. We investigate the predictions of four currently favoured semi-analytical galaxy formation models applied to the Millennium simulation for the mass function of cold neutral gas (atomic and molecular) in galaxies as a function of redshift, and we use these predictions to construct number counts for the next generation of all-sky neutral atomic hydrogen (HI) surveys. Despite the different implementations of the physical ingredients of galaxy formation, we find that the model predictions are broadly consistent with one another; the key differences reflect how the models treat AGN feedback and how the timescale for star formation evolves with redshift. The models produce mass functions of cold gas in galaxies that are generally in good agreement with HI surveys at $z$=0. Interestingly we find that these mass functions do not evolve significantly with redshift. Adopting a simple conversion factor for cold gas mass to HI mass that we apply to all galaxies at all redshifts, we derive mass functions of HI in galaxies from the predicted mass functions of cold gas, which we use to predict the number counts of sources likely to be detected by HI surveys on next generation radio telescopes such as the Square Kilometre Array and its pathfinders. We show also how adopting a conversion factor for cold gas mass to HI mass that varies from galaxy to galaxy impacts on number counts. In addition, we examine how the typical angular sizes of galaxies vary with redshift. Taken together, these results make clear that forthcoming HI surveys will provide important and powerful tests of theoretical galaxy formation models.

We use a model for the evolution of galaxies in the far-IR based on the LambdaCDM cosmology to make detailed predictions for upcoming cosmological surveys with the Herschel Space Observatory. We use the combined GALFORM semi-analytical galaxy formation model and GRASIL spectrophotometric code to compute galaxy SEDs including the reprocessing of radiation by dust. The model, which is the same as that in Baugh et al. (2005), assumes two different IMFs: a normal solar neighbourhood IMF for quiescent star formation in disks, and a very top-heavy IMF in starbursts triggered by galaxy mergers. We have shown previously that the top-heavy IMF appears necessary to explain the number counts and redshifts of faint sub-mm galaxies. In this paper, we present predictions for galaxy luminosity functions, number counts and redshift distributions in the Herschel imaging bands. We find that source confusion will be a serious problem in the deepest planned surveys. We also show predictions for physical properties such as star formation rates and stellar, gas and halo masses, together with fluxes at other wavelengths (from the far-UV to the radio) relevant for multi-wavelength follow-up observations. We investigate what fraction of the total IR emission from dust and of the high-mass star formation over the history of the Universe should be resolved by planned surveys with Herschel, and find a fraction ~30-50%, depending on confusion. Finally, we show that galaxies in Herschel surveys should be significantly clustered.

Future galaxy surveys will map the galaxy distribution in the redshift interval $0.5<z<2$ using near-infrared cameras and spectrographs. The primary science goal of such surveys is to constrain the nature of the dark energy by measuring the large-scale structure of the Universe. This requires a tracer of the underlying dark matter which maximizes the useful volume of the survey. We investigate two potential survey selection methods: an emission line sample based on the \ha line and a sample selected in the H-band. We present predictions for the abundance and clustering of such galaxies, using two published versions of the \galform galaxy formation model. Our models predict that \ha selected galaxies tend to avoid massive dark matter haloes and instead trace the surrounding filamentary structure; H-band selected galaxies, on the other hand, are found in the highest mass haloes. This has implications for the measurement of the rate at which fluctuations grow due to gravitational instability. We use mock catalogues to compare the effective volumes sampled by a range of survey configurations. To give just two examples: a redshift survey down to $H_{\rm AB}=22$ samples an effective volume that is $\sim 5-10$ times larger than that probed by an \ha survey with $\logfha > -15.4$; a flux limit of at least $\logfha = -16$ is required for an \ha sample to become competitive in effective volume.

Future galaxy redshift surveys aim to measure cosmological quantities from the galaxy power spectrum. A prime example is the detection of baryonic acoustic oscillations (BAOs), providing a standard ruler to measure the dark energy equation of state, w(z), to high precision. The strongest practical limitation for these experiments is how quickly accurate redshifts can be measured for sufficient galaxies to map the large-scale structure. A promising strategy is to target emission-line (i.e. star-forming) galaxies at high-redshift (z~0.5-2); not only is the space density of this population increasing out to z~2, but also emission-lines provide an efficient method of redshift determination. Motivated by the prospect of future dark energy surveys targeting H-alpha emitters at near-infrared wavelengths (i.e. z>0.5), we use the latest empirical data to model the evolution of the H-alpha luminosity function out to z~2, and thus provide predictions for the abundance of H-alpha emitters for practical limiting fluxes. We caution that the estimates presented in this work must be tempered by an efficiency factor, epsilon, giving the redshift success rate from these potential targets. For a range of practical efficiencies and limiting fluxes, we provide an estimate of nP_{0.2}, where n is the 3D galaxy number density and P_{0.2} is the galaxy power spectrum evaluated at k=0.2h/Mpc. Ideal surveys must provide nP_{0.2}>1 in order to balance shot-noise and cosmic variance errors. We show that a realistic emission-line survey (epsilon=0.5) could achieve nP_{0.2}=1 out to z~1.5 with a limiting flux of 10^{-16} erg/s/cm^{-2}. If the limiting flux is a factor 5 brighter, then this goal can only be achieved out to z~0.5, highlighting the importance of survey depth and efficiency in cosmological redshift surveys.

We study the nonlinear growth of cosmic structure in different dark energy models, using large volume N-body simulations. We consider a range of quintessence models which feature both rapidly and slowly varying dark energy equations of state, and compare the growth of structure to that in a universe with a cosmological constant. The adoption of a quintessence model changes the expansion history of the universe, the form of the linear theory power spectrum and can alter key observables, such as the horizon scale and the distance to last scattering. We incorporate these effects into our simulations in stages to isolate the impact of each on the growth of structure. The difference in structure formation can be explained to first order by the difference in growth factor at a given epoch; this scaling also accounts for the nonlinear growth at the 15% level. We find that quintessence models that are different from $\Lambda$CDM both today and at high redshifts $(z \sim 1000)$ and which feature late $(z<2)$, rapid transitions in the equation of state, can have identical baryonic acoustic oscillation (BAO) peak positions to those in $\Lambda$CDM. We find that these models have higher abundances of dark matter haloes at $z>0$ compared to $\Lambda$CDM and so measurements of the mass function should allow us to distinguish these quintessence models from a cosmological constant. However, we find that a second class of quintessence models, whose equation of state makes an early $(z>2)$ rapid transition to $w=-1$, cannot be distinguished from $\Lambda$CDM using measurements of the mass function or the BAO, even if these models have non-negligible amounts of dark energy at early times.

We apply a new model for the spherically averaged correlation function at large pair separations to the measurement of the clustering of luminous red galaxies (LRGs) made from the SDSS by Cabre and Gaztanaga(2009). Our model takes into account the form of the BAO peak and the large scale shape of the correlation function. We perform a Monte Carlo Markov chain analysis for different combinations of datasets and for different parameter sets. When used in combination with a compilation of the latest CMB measurements, the LRG clustering and the latest supernovae results give constraints on cosmological parameters which are comparable and in remarkably good agreement, resolving the tension reported in some studies. The best fitting model in the context of a flat, Lambda-CDM cosmology is specified by Omega_m=0.261+-0.013, Omega_b=0.044+-0.001, n_s=0.96+-0.01, H_0=71.6+-1.2 km/s/Mpc and sigma_8=0.80+-0.02. If we allow the time-independent dark energy equation of state parameter to vary, we find results consistent with a cosmological constant at the 5% level using all data sets: w_DE=-0.97+-0.05. The large scale structure measurements by themselves can constrain the dark energy equation of state parameter to w_DE=-1.05+-0.15, independently of CMB or supernovae data. We do not find convincing evidence for an evolving equation of state. We provide a set of "extended distance priors" that contain the most relevant information from the CMB power spectrum and the shape of the LRG correlation function which can be used to constrain dark energy models and spatial curvature. Our model should provide an accurate description of the clustering even in much larger, forthcoming surveys, such as those planned with NASA's JDEM or ESA's Euclid mission.

We introduce a new technique based on artificial neural networks which allows us to make accurate predictions for the spectral energy distributions (SEDs) of large samples of galaxies, at wavelengths ranging from the far-ultra-violet to the sub-millimetre and radio. The neural net is trained to reproduce the SEDs predicted by a hybrid code comprised of the GALFORM semi-analytical model of galaxy formation, which predicts the full star formation and galaxy merger histories, and the GRASIL spectro-photometric code, which carries out a self-consistent calculation of the SED, including absorption and emission of radiation by dust. Using a small number of galaxy properties predicted by GALFORM, the method reproduces the luminosities of galaxies in the majority of cases to within 10% of those computed directly using GRASIL. The method performs best in the sub-mm and reasonably well in the mid-infrared and the far-ultra-violet. The luminosity error introduced by the method has negligible impact on predicted statistical distributions, such as luminosity functions or colour distributions of galaxies. We use the neural net to predict the overlap between galaxies selected in the rest-frame UV and in the observer-frame sub-mm at z=2. We find that around half of the galaxies with a 850um flux above 5 mJy should have optical magnitudes brighter than R_AB < 25 mag. However, only 1% of the galaxies selected in the rest-frame UV down to R_AB < 25 mag should have 850um fluxes brighter than 5 mJy. Our technique will allow the generation of wide-angle mock catalogues of galaxies selected at rest-frame UV or mid- and far-infrared wavelengths.

The huge size and uniformity of the Sloan Digital Sky Survey makes possible an exacting test of current models of galaxy formation. We compare the predictions of the GALFORM semi-analytical galaxy formation model for the luminosities, morphologies, colours and scale-lengths of local galaxies. GALFORM models the luminosity and size of the disk and bulge components of a galaxy, and so we can compute quantities which can be compared directly with SDSS observations, such as the Petrosian magnitude and the Sersic index. We test the predictions of two published models set in the cold dark matter cosmology: the Baugh et al. (2005) model, which assumes a top-heavy initial mass function (IMF) in starbursts and superwind feedback, and the Bower et al. (2006) model, which uses AGN feedback and a standard IMF. The Bower et al model better reproduces the overall shape of the luminosity function, the morphology-luminosity relation and the colour bimodality observed in the SDSS data, but gives a poor match to the size-luminosity relation. The \Baugh et al. model successfully predicts the size-luminosity relation for late-type galaxies. Both models fail to reproduce the sizes of bright early-type galaxies. These problems highlight the need to understand better both the role of feedback processes in determining galaxy sizes, in particular the treatment of the angular momentum of gas reheated by supernovae, and the sizes of the stellar spheroids formed by galaxy mergers and disk instabilities.

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.

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.

We combine a semi-analytical model of galaxy formation with a very large simulation which follows the growth of large scale structure in a LambdaCDM universe to predict the clustering of Ly-alpha emitters. We find that the clustering strength of Ly-alpha emitters has only a weak dependence on Ly-alpha luminosity but a strong dependence on redshift. With increasing redshift, Ly-alpha emitters trace progressively rarer, higher density regions of the universe. Due to the large volume of the simulation, over 100 times bigger than any previously used for this application, we can construct mock catalogues of Ly-alpha emitters and study the sample variance of current and forthcoming surveys. We find that the number and clustering of Ly-alpha emitters in our mock catalogues are in agreement with measurements from current surveys, but that there is a considerable scatter in these quantities. We argue that a proposed survey of emitters at z=8.8 should be extended significantly in solid angle to allow a robust measurement of Ly-alpha emitter clustering.

[ABRIDGED VERSION] The future of cm and m-wave astronomy lies with the Square Kilometre Array (SKA), a telescope under development by a consortium of 17 countries. The SKA will be 50 times more sensitive than any existing radio facility. A majority of the key science for the SKA will be addressed through large-area imaging of the Universe at frequencies from 300 MHz to a few GHz. The Australian SKA Pathfinder (ASKAP) is aimed squarely in this frequency range, and achieves instantaneous wide-area imaging through the development and deployment of phase-array feed systems on parabolic reflectors. This large field-of-view makes ASKAP an unprecedented synoptic telescope poised to achieve substantial advances in SKA key science. The central core of ASKAP will be located at the Murchison Radio Observatory in inland Western Australia, one of the most radio-quiet locations on the Earth and one of the sites selected by the international community as a potential location for the SKA. Following an introductory description of ASKAP, this document contains 7 chapters describing specific science programmes for ASKAP. The combination of location, technological innovation and scientific program will ensure that ASKAP will be a world-leading radio astronomy facility, closely aligned with the scientific and technical direction of the SKA. A brief summary chapter emphasizes the point, and considers discovery space.

We use the Millennium Simulation, a large, high resolution N-body simulation of the evolution of structure in a LambdaCDM cosmology, to study the properties and fate of substructures within a large sample of dark matter haloes. We find that the subhalo mass function departs significantly from a power law at the high mass end. We also find that the radial and angular distributions of substructures depend on subhalo mass. In particular, high mass subhaloes tend to be less radially concentrated and to have angular distributions closer to the direction perpendicular to the spin of the host halo than their less massive counterparts. We find that mergers between subhaloes occur. These tend to be between substructures that were already dynamically associated before accretion into the main halo. For subhaloes larger than 0.001 times the mass of the host halo, it is more likely that the subhalo will merge with the central or main subhalo than with another subhalo larger than itself. For lower masses, subhalo-subhalo mergers become equally likely to mergers with the main subhalo. Our results have implications for the variation of galaxy properties with environment and for the treatment of mergers in galaxy formation models.

We present a test of different error estimators for 2-point clustering statistics, appropriate for present and future large galaxy redshift surveys. Using an ensemble of very large dark matter LambdaCDM N-body simulations, we compare internal error estimators (jackknife and bootstrap) to external ones (Monte-Carlo realizations). For 3-dimensional clustering statistics, we find that none of the internal error methods investigated are able to reproduce neither accurately nor robustly the errors of external estimators on 1 to 25 Mpc/h scales. The standard bootstrap overestimates the variance of xi(s) by ~40% on all scales probed, but recovers, in a robust fashion, the principal eigenvectors of the underlying covariance matrix. The jackknife returns the correct variance on large scales, but significantly overestimates it on smaller scales. This scale dependence in the jackknife affects the recovered eigenvectors, which tend to disagree on small scales with the external estimates. Our results have important implications for the use of galaxy clustering in placing constraints on cosmological parameters. For example, in a 2-parameter fit to the projected correlation function, we find that the standard bootstrap systematically overestimates the 95% confidence interval, while the jackknife method remains biased, but to a lesser extent. The scatter we find between realizations, for Gaussian statistics, implies that a 2-sigma confidence interval, as inferred from an internal estimator, could correspond in practice to anything from 1-sigma to 3-sigma. Finally, by an oversampling of sub-volumes, it is possible to obtain bootstrap variances and confidence intervals that agree with external error estimates, but it is not clear if this prescription will work for a general case.

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.

We use the combined GALFORM semi-analytical model of galaxy formation and GRASIL spectrophotometric code to investigate the properties of galaxies selected via their sub-mm emission. Our fiducial model has previously been shown to fit the properties of local ULIRGs, as well as the number counts of faint sub-mm galaxies. Here, we test the model in detail by comparing the SEDs and stellar, dynamical, gas and halo masses of sub-mm galaxies against observational data. We precisely mimic the sub-mm and radio selection function of the observations and show that the predicted far-infrared properties of model galaxies with S_850>5mJy and S_1.4>30uJy are in good agreement with observations. Although the dust emission model does not assume a single dust temperature, the far-infrared SEDs are well described by single component modified black-body spectrum with characteristic temperature 32+/-5K. We also find evidence that the observations may have uncovered evolution in the far-infrared--radio relation in ULIRGs out to z~2. We show that the predicted redshift distribution of sub-mm galaxies provides a reasonable fit to the observational data with a median redshift z=2.0, with the radio-selected subset predicted to make up approximately 75% of the population. However, the predicted K-band and mid-infrared (3--8um) flux densities of the sub-mm galaxies (and LBGs) are up to a factor 10x fainter than observed. This discrepancy may indicate that the stellar masses of the sub-mm galaxies in the model are too low: M~10^10Mo, while observations suggest more massive systems, M~10^11Mo. Finally, we discuss the potential modifications to the models which may improve the fit to the observational data. [Abridged]

We use an extremely large volume ($2.4h^{-3}{\rm Gpc}^{3}$), high resolution N-body simulation to measure the higher order clustering of dark matter haloes as a function of mass and internal structure. As a result of the large simulation volume and the use of a novel ``cross-moment'' counts-in-cells technique which suppresses discreteness noise, we are able to measure the clustering of haloes corresponding to rarer peaks than was possible in previous studies; the rarest haloes for which we measure the variance are 100 times more clustered than the dark matter. We are able to extract, for the first time, halo bias parameters from linear up to fourth order. For all orders measured, we find that the bias parameters are a strong function of mass for haloes more massive than the characteristic mass $M_{*}$. Currently, no theoretical model is able to reproduce this mass dependence closely. We find that the bias parameters also depend on the internal structure of the halo up to fourth order. For haloes more massive than $M_{*}$, we find that the more concentrated haloes are more weakly clustered than the less concentrated ones. We see no dependence of clustering on concentration for haloes with masses $M<M_{*}$; this is contrary to the trend reported in the literature when segregating haloes by their formation time. Our results are insensitive to whether haloes are labelled by the total mass returned by the friends-of-friends group finder or by the mass of the most massive substructure. This implies that our conclusions are not an artefact of the particular choice of group finding algorithm. Our results will provide important input to theoretical models of galaxy clustering.

Oscillations in the baryon-photon fluid prior to recombination imprint different signatures on the power spectrum and correlation function of matter fluctuations. The measurement of these features using galaxy surveys has been proposed as means to determine the equation of state of the dark energy. The accuracy required to achieve competitive constraints demands an extremely good understanding of systematic effects which change the baryonic acoustic oscillation (BAO) imprint. We use 50 very large volume N-body simulations to investigate the BAO signature in the two-point correlation function. The location of the BAO bump does not correspond to the sound horizon scale at the level of accuracy required by future measurements, even before any dynamical or statistical effects are considered. Careful modelling of the correlation function is therefore required to extract the cosmological information encoded on large scales. We find that the correlation function is less affected by scale dependent effects than the power spectrum. We show that a model for the correlation function proposed by Crocce & Scoccimarro (2008), based on renormalised perturbation theory, gives an essentially unbiased measurement of the dark energy equation of state. This means that information from the large scale shape of the correlation function, in addition to the form of the BAO peak, can be used to provide robust constraints on cosmological parameters. The correlation function therefore provides a better constraint on the distance scale (~50% smaller errors with no systematic bias) than the more conservative approach required when using the power spectrum (i.e. which requires amplitude and long wavelength shape information to be discarded).

We present an analysis of the small-to-intermediate scale clustering of samples of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey and the 2dF-SDSS LRG and QSO (2SLAQ) survey carefully matched to have the same rest-frame colours and luminosity. We study the spatial two-point auto-correlation function in both redshift-space and real-space of a combined sample of over 10,000 LRGs, which represent the most massive galaxies in the universe with stellar masses > 10^11 h^-1 M_sun and space densities 10^-4 h^-3 Mpc^-3. We find no significant evolution in the amplitude r_0 of the correlation function with redshift, but do see a slight decrease in the slope with increasing redshift over 0.19 < z < 0.55 and scales of 0.32 < r < 32 h^-1 Mpc. We compare our measurements with the predicted evolution of dark matter clustering and use the halo model to interpret our results. We find that our clustering measurements are inconsistent (>99.9% significance) with a passive model whereby the LRGs do not merge with one another; a model with a merger rate of 7.5 +/- 2.3% from z = 0.55 to z = 0.19 (i.e. an average rate of 2.4% Gyr^-1) provides a better fit to our observations. Our clustering and number density measurements are consistent with the hypothesis that the merged LRGs were originally central galaxies in different haloes which, following the merger of these haloes, merged to create a single Brightest Cluster Galaxy. In addition, we show that the small-scale clustering signal constrains the scatter in halo merger histories. When combined with measurements of the luminosity function, our results suggest that this scatter is sub-Poisson. While this is a generic prediction of hierarchical models, it has not been tested before.

We describe the scientific motivations, the mission concept and the instrumentation of SPACE, a class-M mission proposed for concept study at the first call of the ESA Cosmic-Vision 2015-2025 planning cycle. SPACE aims to produce the largest three-dimensional evolutionary map of the Universe over the past 10 billion years by taking near-IR spectra and measuring redshifts for more than half a billion galaxies at 0<z<2 down to AB~23 over 3\pi sr of the sky. In addition, SPACE will also target a smaller sky field, performing a deep spectroscopic survey of millions of galaxies to AB~26 and at 2<z<10+. These goals are unreachable with ground-based observations due to the ~500 times higher sky background. To achieve the main science objectives, SPACE will use a 1.5m diameter Ritchey-Chretien telescope equipped with a set of arrays of Digital Micro-mirror Devices (DMDs) covering a total field of view of 0.4 deg2, and will perform large-multiplexing multi-object spectroscopy (e.g. ~6000 targets per pointing) at a spectral resolution of R~400 as well as diffraction-limited imaging with continuous coverage from 0.8mum to 1.8mum.

Luminous red galaxies (LRGs) are much rarer and more massive than L* galaxies. Coupled with their extreme colours, LRGs therefore provide a demanding testing ground for the physics of massive galaxy formation. We present the first self-consistent predictions for the abundance and properties of LRGs in hierarchical structure formation models. We test two published models which use quite different mechanisms to suppress the formation of massive galaxies: the Bower et al. (2006) model, which invokes ``AGN-feedback'' to prevent gas from cooling in massive haloes, and the Baugh et al. (2005) model which relies upon a ``superwind'' to eject gas before it is turned into stars. Without adjusting any parameters, the Bower et al. model gives an excellent match to the observed luminosity function of LRGs in the SDSS (with a median redshift of z=0.24) and to their clustering; the Baugh et al. model is less successful in these respects. Both models fail to match the observed abundance of LRGs at z=0.5 to better than a factor of ~2. In the models, LRGs are typically bulge dominated systems with M* of ~2x10^11 h^{-1} M_sun and velocity dispersions of ~250 km s^{-1}. Around half of the stellar mass in the model LRGs is already formed by z~2.2 and is assembled into one main progenitor by z~1.5; on average, only 25% of the mass of the main progenitor is added after z~1. LRGs are predicted to be found in a wide range of halo masses, a conclusion which relies on properly taking into account the scatter in the formation histories of haloes. Remarkably, we find that the correlation function of LRGs is predicted to be a power law down to small pair separations, in excellent agreement with observational estimates. Neither the Bower et al. nor the Baugh et al. model is able to reproduce the observed radii of LRGs.

We present predictions for the evolution of the galaxy luminosity function, number counts and redshift distributions in the IR based on the Lambda-CDM cosmological model. We use the combined GALFORM semi-analytical galaxy formation model and GRASIL spectrophotometric code to compute galaxy SEDs including the reprocessing of radiation by dust. The model, which is the same as that in Baugh et al (2005), assumes two different IMFs: a normal solar neighbourhood IMF for quiescent star formation in disks, and a very top-heavy IMF in starbursts triggered by galaxy mergers. We have shown previously that the top-heavy IMF seems to be necessary to explain the number counts of faint sub-mm galaxies. We compare the model with observational data from the SPITZER Space Telescope, with the model parameters fixed at values chosen before SPITZER data became available. We find that the model matches the observed evolution in the IR remarkably well over the whole range of wavelengths probed by SPITZER. In particular, the SPITZER data show that there is strong evolution in the mid-IR galaxy luminosity function over the redshift range z ~ 0-2, and this is reproduced by our model without requiring any adjustment of parameters. On the other hand, a model with a normal IMF in starbursts predicts far too little evolution in the mid-IR luminosity function, and is therefore excluded.

We incorporate a model for black hole growth during galaxy mergers into the semi-analytical galaxy formation model based on Lambda-CDM proposed by Baugh et al. (2005). Our black hole model has one free parameter, which we set by matching the observed zeropoint of the local correlation between black hole mass and bulge luminosity. We present predictions for the evolution with redshift of the relationships between black hole mass and bulge properties. Our simulations reproduce the evolution of the optical luminosity function of quasars. We study the demographics of the black hole population and address the issue of how black holes acquire their mass. We find that the direct accretion of cold gas during starbursts is an important growth mechanism for lower mass black holes and at high redshift. On the other hand, the re-assembly of pre-existing black hole mass into larger units via merging dominates the growth of more massive black holes at low redshift. This prediction could be tested by future gravitational wave experiments. As redshift decreases, progressively less massive black holes have the highest fractional growth rates, in line with recent claims of "downsizing" in quasar activity.

We assess the detectability of baryonic acoustic oscillations (BAO) in the power spectrum of galaxies using ultra large volume N-body simulations of the hierarchical clustering of dark matter and semi-analytical modelling of galaxy formation. A step-by-step illustration is given of the various effects (nonlinear fluctuation growth, peculiar motions, nonlinear and scale dependent bias) which systematically change the form of the galaxy power spectrum on large scales from the simple prediction of linear perturbation theory. Using a new method to extract the scale of the oscillations, we nevertheless find that the BAO approach gives an unbiased estimate of the sound horizon scale. Sampling variance remains the dominant source of error despite the huge volume of our simulation box ($=2.41 h^{-3}{\rm Gpc}^{3}$). We use our results to forecast the accuracy with which forthcoming surveys will be able to measure the sound horizon scale, $s$, and, hence constrain the dark energy equation of state parameter, $w$ (with simplifying assumptions and without marginalizing over the other cosmological parameters). Pan-STARRS could potentially yield a measurement with an accuracy of $\Delta s/s = 0.5-0.7 % $ (corresponding to $\Delta w \approx 2-3% $), which is competitive with the proposed WFMOS survey ($\Delta s/s = 1% $ $\Delta w \approx 4 % $). Achieving $\Delta w \le 1% $ using BAO alone is beyond any currently commissioned project and will require an all-sky spectroscopic survey, such as would be undertaken by the SPACE mission concept under proposal to ESA.

Context: Many current and future surveys aim to detect the highest redshift (z >~ 7) sources through their Lyman-alpha (Ly-alpha) emission, using the narrow-band imaging method. However, to date the surveys have only yielded non-detections and upper limits as no survey has reached the necessary combination of depth and area to detect these very young star forming galaxies. Aims: We aim to calculate model luminosity functions and mock surveys of Ly-alpha emitters at z >~ 7 based on a variety of approaches. Methods: We calculate model luminosity functions at different redshifts based on three different approaches: a semi-analytical model based on CDM, a simple phenomenological model, and an extrapolation of observed Schechter functions at lower redshifts. The results of the first two models are compared with observations made at redshifts z ~ 5.7 and z ~ 6.5, and they are then extrapolated to higher redshift. Results: We present model luminosity functions for redshifts between z = 7 - 12.5 and give specific number predictions for future planned or possible narrow-band surveys for Ly-alpha emitters. We also investigate what constraints future observations will be able to place on the Ly-alpha luminosity function at very high redshift. Conclusion: It should be possible to observe z = 7 - 10 Ly-alpha emitters with present or near-future instruments if enough observing time is allocated. In particular, large area surveys such as ELVIS (Emission Line galaxies with VISTA Survey) will be useful in collecting a large sample. However, to get a large enough sample to constrain well the z >= 10 Ly-alpha luminosity function, instruments further in the future, such as an ELT, will be necessary.

We present measurements of the higher-order clustering of red and blue galaxies as a function of scale and luminosity made from the two-degree field galaxy redshift survey (2dFGRS). We use a counts-in-cells analysis to estimate the volume averaged correlation functions, xi_p, as a function of scale up to order p=5, and also the reduced void probability function. Hierarchical amplitudes are constructed using the estimates of the correlation functions: S_p=(xi_p/xi_2)^(p-1). We find that: 1) Red galaxies display stronger clustering than blue galaxies at all orders measured. 2) Red galaxies show values of S_p that are strongly dependent on luminosity, whereas blue galaxies show no segregation in S_p within the errors; this is remarkable given the segregation in the variance. 3) The linear relative bias shows the opposite trend to the hierarchical amplitudes, with little segregation for the red sequence and some segregation for the blue. 4) Faint red galaxies deviate significantly from the "universal" negative binomial reduced void probabilities followed by all other galaxy populations. Our results show that the characteristic colour of a galaxy population reveals a unique signature in its spatial distribution. Such signatures will hopefully further elucidate the physics responsible for shaping the cosmological evolution of galaxies.

We present predictions for the structural and photometric properties of early-type galaxies in the cold dark matter cosmology from the published semi-analytical galaxy formation models of Baugh et al. and Bower et al. These calculations were made with the GALFORM code, which tracks the evolution of the disc and bulge components of a galaxy, using a self-consistent model to compute the scalelengths. The sizes of galactic discs are determined by the conservation of the angular momentum of cooling gas. The sizes of merger remnants are computed by applying the virial theorem and conserving the binding energy of the progenitors and their orbital energy. We compare the model predictions with observational results derived from the SDSS. The model enjoys a number of notable successes, such as giving reasonable reproductions of the local Faber-Jackson relation, the velocity dispersion-age relation, and the fundamental plane relating the luminosity, velocity dispersion and effective radius of spheroids. These achievements are all the more remarkable when one bears in mind that none of the parameters have been adjusted to refine the model predictions. We study how the residuals around the fundamental plane relation depend on galaxy properties. We examine in detail the physical ingredients of the calculation of galaxy sizes. We also study the evolution of the scaling relations with redshift. However, there are some important disagreements between the predictions of the model and observations: the brightest model spheroids have effective radii smaller than observed and the zero-point of the fundamental plane shows little or no evolution with redshift in the model.

Recent measurements of the cosmic microwave background radiation (CMB), particularly when combined with other datasets, have revolutionised our knowledge of the values of the basic cosmological parameters. Here we summarize the state of play at the end of 2006, focusing on the combination of CMB measurements with the power spectrum of galaxy clustering. We compare the constraints derived from the extant CMB data circa 2005 and the final 2dFGRS galaxy power spectrum, with the results obtained when the WMAP 1-year data is replaced by the 3-year measurements (hereafter WMAP1 and WMAP3). Remarkably, the picture has changed relatively little with the arrival of WMAP3, though some aspects have been brought into much sharper focus. One notable example of this is the index of primordial scalar fluctuations, n_s. Prior to WMAP3, Sanchez et al. (2006) found that the scale invariant value of n_s = 1 was excluded at the 95% level. With WMAP3, this becomes a 3sigma result, with implications for models of inflation. We find some disagreement between the constraints on certain parameters when the 2dFGRS P(k) is replaced by the SDSS measurement. This suggests that more work is needed to understand the relation between the clustering of different types of galaxies and the linear perturbation theory prediction for the power spectrum of matter fluctuations.

Recent observational and theoretical breakthroughs make this an exciting time to be working towards understanding the physics of galaxy formation. The goal of this review is to make the principles behind the hierarchical paradigm accessible to a wide audience by providing a pedagogical introduction to modern theories of galaxy formation. I outline the ingredients of the powerful approach called semi-analytical modelling and contrast this method with numerical simulations of the gas dynamics of baryons. Semi-analytical models have enjoyed many successes, but it is the observations which the models struggle to match which mark out areas where future progress is most likely to be made; these are also reviewed.

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 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.

We present detailed predictions for the properties of Ly-alpha-emitting galaxies in the framework of the Lambda-CDM cosmology, calculated using the semi-analytical galaxy formation model GALFORM. We explore a model which assumes a top-heavy IMF in starbursts, and which has previously been shown to explain the sub-mm number counts and the luminosity function of Lyman-break galaxies at high redshift. We show that this model, with the simple assumption that a fixed fraction of Ly-alpha photons escape from each galaxy, is remarkably successful at explaining the observed luminosity function of Ly-alpha emitters over the redshift range 3<z<6.6. We also examine the distribution of Ly-alpha equivalent widths and the broad-band continuum magnitudes of emitters, which are in good agreement with the available observations. We look more deeply into the nature of Ly-alpha emitters, presenting predictions for fundamental properties such as the stellar mass and radius of the emitting galaxy and the mass of the host dark matter halo. The model predicts that the clustering of Ly-alpha emitters at high redshifts should be strongly biased relative to the dark matter, in agreement with observational estimates. We also present predictions for the luminosity function of Ly-alpha emitters at z>7, a redshift range which is starting to be be probed by near-IR surveys and using new instruments such as DAzLE.

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.

We present new results for the 3-point correlation function, \zeta, measured as a function of scale, luminosity and colour from the final version of the two-degree field galaxy redshift survey (2dFGRS). The reduced three point correlation function, Q_3 is estimated for different triangle shapes and sizes, employing a full covariance analysis. The form of Q_3 is consistent with the expectations for the \Lambda-cold dark matter model, confirming that the primary influence shaping the distribution of galaxies is gravitational instability acting on Gaussian primordial fluctuations. However, we find a clear offset in amplitude between Q_3 for galaxies and the predictions for the dark matter. We are able to rule out the scenario in which galaxies are unbiased tracers of the mass at the 9-sigma level. On weakly non-linear scales, we can interpret our results in terms of galaxy bias parameters. We find a linear bias term that is consistent with unity, b_1 = 0.93^{+0.10}_{-0.08} and a quadratic bias c_2 = b_2 /b_1 = -0.34^{+0.11}_{-0.08}. This is the first significant detection of a non-zero quadratic bias, indicating a small but important non-gravitational contribution to the three point function. Our estimate of the linear bias from the three point function is independent of the normalisation of underlying density fluctuations, so we can combine this with the measurement of the power spectrum of 2dFGRS galaxies to constrain the amplitude of matter fluctuations. We find that the rms linear theory variance in spheres of radius 8Mpc/h is \sigma_8 = 0.88^{+0.12}_{-0.10}, providing an independent confirmation of values derived from other techniques. On non-linear scales, where \xi>1, we find that Q_3 has a strong dependence on scale, colour and luminosity.

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.

Acoustic oscillations in the baryon-photon fluid leave a signature in the matter power spectrum. The overall shape of the spectrum and the wavelength of the oscillations depend upon the sound horizon scale at recombination. Using the $\Lambda$ cold dark matter Hubble Volume simulation, we show that the imprint of baryons is visible in the power spectrum of cluster-mass dark matter haloes, in spite of significant differences between the halo power spectrum and the prediction of linear perturbation theory. A measurement of the sound horizon scale can constrain the dark energy equation of state. We show that a survey of clusters at intermediate redshift ($ z\sim1 $), like the Sunyaev-Zeldovich survey proposed by the South Pole Telescope or a red sequence photometric survey with VISTA, could potentially constrain the sound horizon scale to an accuracy of $\sim 2%$, in turn fixing the ratio of the pressure of the dark energy to its density ($w$) to better than $\sim 10%$. Our approach does not require knowledge of the cluster mass, unlike those that depend upon the abundance of clusters.

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 present predictions for the abundance of Ly-$\alpha$ emitters in hierarchical structure formation models. We use the {\tt GALFORM} semi-analytical model to explore the impact on the predicted counts of varying assumptions about the escape fraction of Ly-$\alpha$ photons, the redshift at which the universe reionised and the cosmological density parameter. A model with a fixed escape fraction gives a remarkably good match to the observed counts over a wide redshift interval. We present predictions for the expected counts in a typical observation with the Multi Unit Spectroscopic Explorer instrument proposed for the Very Large Telescope.

The present-day Universe is seemingly dominated by dark energy and dark matter, but mapping the normal (baryonic) content remains vital for both astrophysics - understanding how galaxies form - and astro-particle physics - inferring properties of the dark components. The Square Kilometre Array (SKA) will provide the only means of studying the cosmic evolution of neutral Hydrogen (HI) which, alongside information on star formation from the radio continuum, is needed to understand how stars formed from gas within dark-matter over-densities and the roles of gas accretion and galaxy merging. `All hemisphere' HI redshift surveys to redshift 1.5 are feasible with wide-field-of-view realizations of the SKA and, by measuring the galaxy power spectrum in exquisite detail, will allow the first precise studies of the equation-of-state of dark energy. The SKA will be capable of other uniquely powerful cosmological studies including the measurement of the dark-matter power spectrum using weak gravitational lensing, and the precise measurement of H0 using extragalactic water masers. The SKA is likely to become the premier dark-energy-measuring machine, bringing breakthroughs in cosmology beyond those likely to be made possible by combining CMB (e.g. Planck), optical (e.g. LSST, SNAP) and other early-21st-century datasets.

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 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 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 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 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 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.

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)

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.

We measure the power spectrum of galaxy clustering in real space from the APM Galaxy Survey. We present an improved technique for the numerical inversion of Limber's equation that relates the angular clustering of galaxies to an integral over the power spectrum in three dimensions. Our approach is underpinned by a large ensemble of mock galaxy catalogues constructed from the Hubble Volume N-body simulations. The mock catalogues are used to test for systematic effects in the inversion algorithm and to estimate the errors on our measurement. We find that we can recover the power spectrum to an accuracy of better than 15% over three decades in wavenumber. A key advantage of the use of mock catalogues to infer errors is that we can apply our technique on scales for which the density fluctuations are not Gaussian, thus probing the regime that offers the best constraints on models of galaxy formation. On large scales, our measurement of the power spectrum is consistent with the shape of the mass power spectrum in the popular ``concordance'' cold dark matter model. The galaxy power spectrum on small scales is strongly affected by nonlinear evolution of density fluctuations, and, to a lesser degree, by galaxy bias. The rms variance in the galaxy distribution, when smoothed in spheres of radius 8 h^(-1)Mpc, is sigma^g_8=0.96^(+0.17)_(-0.20) and the shape of the power spectrum on large scales is described by a simple fitting formula with parameter Gamma= 0.19^(+0.13)_(-0.04) (these errors are the 1-sigma ranges for a two parameter fit). We use our measurement of the power spectrum to estimate the galaxy two point correlation function; the results are well described by a power law with correlation length r_0 = 5.9 +- 0.7 h^(-1)Mpc and slope gamma = 1.61 +- 0.06 for pair separations in the range 0.1 < r/[h^(-1)Mpc] < 20.

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.

We apply the semi-analytical galaxy formation code of Cole et al. to investigate the formation of globular cluster (GC) systems in hierarchical clustering scenarios. The nature of the model allows us to investigate the properties of GC systems and their parent galaxies within a cosmological framework, over a wide dynamic range of mass and time resolution. Assuming GCs form during mergers of gaseous systems, the metal-rich peak of the classical 'bimodal' metallicity distribution of GCs naturally falls out of our model, where such merging occurs over a wide range of redshifts. The physical origin of old, metal-poor GCs (the metal-poor peak) is harder to understand, since their formation must be decoupled from the ongoing star formation in these systems at high redshift (z~5). Within the context of semi-analytic models in general, a possible solution lies in a cut-off in the GC formation efficiency at a characteristic local star formation rate.

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 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 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 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 describe a scheme for the formation of globular cluster systems in early-type galaxies using a semi-analytic model of galaxy formation. Operating within a Lambda CDM cosmology, we assume that metal-poor globular clusters are formed at high-redshift in pre-galactic fragments, and that the subsequent gas-rich merging of these fragments leads to the formation of metal-rich clusters. We compare our results with contemporary data, and look at the particular case of the globular cluster and stellar metallicity distribution function of the nearby elliptical galaxy Centaurus A.

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 have investigated the formation of globular cluster (GC) systems in the fiducial semi-analytic galaxy formation model of Cole et al. We assume that GCs are formed at high-redshift (z > 5) in proto-galactic fragments, and during the subsequent gas-rich merging of these fragments. We find that the mean ages of the metal-rich GCs are dependent upon host galaxy luminosity and halo circular velocity, whereas the mean ages of the metal-poor GCs are not. We find that gaseous merging of the proto-galactic fragments leads to a significant age-spread amongst the metal-rich GCs, which increases for lower-luminosity galaxies. We find the total GC populations scale with host galaxy luminosity as L^1.25, due to the systematic increase in the mass-to-light ratios (with mass) of the galaxy haloes in the model. A comparison between the luminosity growth of the model ellipticals and their GC formation indicates that mergers do not significantly effect S_N at z < 2. We conclude that gaseous merging, the bulk of which occurs at 1 < z < 4 in the Lambda-CDM model, leads to the formation of the metal-rich peak of the GC systems of elliptical galaxies. We suggest that the formation and subsequent truncation of the metal-poor GCs is closely related to the star formation rate in the proto-galactic fragments, which may have been significantly higher at early times (abridged).

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 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 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 test an analytic model for the two-point correlations of galaxy clusters in redshift space using the Hubble Volume N-body simulations. The correlation function of clusters shows no enhancement along the line of sight, due to the lack of any virialised structures in the cluster distribution. However, the distortion of the clustering pattern due to coherent bulk motions is clearly visible. The distribution of cluster peculiar motions is well described by a Gaussian, except in the extreme high velocity tails. The simulations produce a small but significant number of clusters with large peculiar motions. The form of the redshift space power spectrum is strongly influenced by errors in measured cluster redshifts in extant surveys. When these errors are taken into account, the model reproduces the power spectrum recovered from the simulation to an accuracy of 15% or better over a decade in wavenumber. We compare our analytic predictions with the power spectrum measured from t! he! APM cluster redshift survey. The cluster power spectrum constrains the amplitude of density fluctuations, as measured by the linear rms variance in spheres of radius 8 h^(-1) Mpc, denoted by sigma_8. When combined with the constraints on sigma_8 and the density parameter Omega derived from the local abundance of clusters, we find a best fitting cold dark matter model with sigma_8~1.25 and Omega~0.2, for a power spectrum shape that matches that measured for galaxies. However, for the best fitting value of Omega and given the value of Hubble's constant from recent measurements, the assumed shape of the power spectrum is incompatible with the most readily motivated predictions from the cold dark matter paradigm.

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.

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.

We describe an instrument concept that will provide simultaneous spectra for a million spatial samples on the sky. With the proposed field of view and spectral resolution, it will be able to measure redshifts and line strengths for around 2-4000 z~3-7 galaxies in a 16 night campaign. The main science driver is to obtain a complete census of the star formation properties of galaxies with 2.5<z<6.7 as a function of luminosity, environment, morphology and redshift. This survey will also allow us to study the evolution of the 3 dimensional power spectrum and test the connection between these early galaxies and the large scale structure of the Lyman Alpha forest. Additional science drivers include [OII] emitters with 0.13<z<1.52, abundance gradients in nearby galaxies, searches for planetary nebulae in nearby galaxy groups and clusters, lensed emission line objects, and spectroscopic searches for objects with continuum breaks. We expect that this instrument will make an ideal partner to the HST ACS, where ACS morphology and broadband colour information can complement the MEIFU spectroscopy and help remove any degeneracies in line identification.

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 present a counts-in-cells analysis of clustering in the optically selected Durham/UKST and Stromlo-APM Galaxy Redshift Surveys. Minimum variance estimates of the second moment, skewness (S_3) and kurtosis (S_4) of the count probability distribution are extracted from a series of volume limited samples of varying radial depth. The corresponding theoretical error calculation takes into account all sources of statistical error on the measurement of the moments, and is in good agreement with the dispersion over mock redshift catalogues. The errors that we find on S_3 and S_4 are larger than those quoted in previous studies, in spite of the fact that the surveys we consider cover larger volumes. S_3 varies little with cell size, with values in the range 1.8-2.2 and errors ~20%, for cubical cells of side 3-20 h^{-1}Mpc. Direct measurements of S_3 are possible out to ~35 h^{-1}Mpc, though with larger errors. A significant determination of S_4 is only possible for one scale, l=6 h^{-1}Mpc, with S_4~5. We compare our results with theoretical predictions from N-body simulations of cold dark matter universes. Qualitatively, the skewness of the dark matter has the same form as that of the galaxies. However, the amplitude of the galaxy S_3 is lower than that predicted for the dark matter. Our measurements of S_3 are consistent with the predictions of a simple model in which initially Gaussian fluctuations in the dark matter evolve gravitationally, if a second order bias term is specified, in addition to the traditional linear bias, in order to describe the relation between the distribution of galaxies and dark matter.

We compare the statistical properties of galaxies found in two different models of hierarchical galaxy formation: the semi-analytic model of Cole et al. and the smoothed particle hydrodynamics (SPH) simulations of Pearce et al. Using a `stripped-down' version of the semi-analytic model which mimics the resolution of the SPH simulations and excludes physical processes not included in them, we find that the two models produce an ensemble of galaxies with remarkably similar properties, although there are some differences in the gas cooling rates and in the number of galaxies that populate halos of different mass. The full semi-analytic model, which has effectively no resolution limit and includes a treatment of star formation and supernovae feedback, produces somewhat different (but readily understandable) results. Agreement is particularly good for the present-day global fractions of hot gas, cold dense (i.e. galactic) gas and uncollapsed gas, for which the SPH and stripped-down semi-analytic calculations differ by at most 25%. In the most massive halos, the stripped-down semi-analytic model predicts, on the whole, up to 50% less gas in galaxies than is seen in the SPH simulations. The two techniques apportion this cold gas somewhat differently amongst galaxies in a given halo. This difference can be tracked down to the greater cooling rate in massive halos in the SPH simulation compared to the semi-analytic model. (abridged)

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.

We present the power spectrum analysis of clustering in the Durham/UKST Galaxy Redshift Survey, which covers 1450 square degrees and consists of 2501 galaxy redshifts sampled at a rate of 1 in 3 to b_J = 17. The fluctuations that we measure can be expressed as sigma_{8} = 1.01 +/- 0.17 . We find remarkably good agreement between the power spectrum measured for the Durham/UKST Survey and those from other optical studies, in spite of the very different sampling rates, magnitude limits and survey geometries employed, on scales up to lambda = (2 pi)/k = 80/h Mpc. On scales larger than this we find good agreement with the power measured from the Stromlo-APM Survey (Tadros & Efstathiou), but find more power than estimated from the Las Campanas Redshift Survey (Lin et al). We apply a simple model for redshift space distortions to the real space APM Galaxy Survey power spectrum and find a shape and amplitude that is in very good agreement with the power spectrum of the Durham/UKST Survey. This implies beta=Omega^{0.6}/b=0.60 +/- 0.35, where Omega is the density parameter and b is the bias between fluctuations in the galaxy and mass distributions, and also suggests a one dimensional velocity dispersion of sigma = 320 +/-140 km/s. We find that for any choice of normalisation, the standard CDM model has a shape that cannot be reconciled with the Durham/UKST Survey power spectrum, unless either unacceptably high values of the one dimensional velocity dispersion are adopted or the assumption that bias is constant is invalid on scales greater than 20/h Mpc. Over the range of wavenumbers for which we have a robust measurement of the power spectrum, we find the best agreement is obtained for a critical density CDM model in which the shape of the power spectrum is modified. (Abridged)

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 study the evolution of galaxy satellites with high resolution N-body simulations. Satellites are modeled as replicas of typical low and high surface brightness galaxies (LSBs and HSBs). Encounters on high eccentricity orbits (as typical in hierarchical models of galaxy formation) strip LSBs of most of their stars and tend to decrease their surface brightness. On the contrary, bar instability in HSBs leads to substantial loss of angular momentum of the stellar component and to an increase of central surface brightness. In both cases the remnant resembles a spheroidal galaxy with an exponential surface brightness profile. A simple modeling of color evolution and interactions driven star formation gives M/L ratios for the remnants that are roughly consistent with observations. These results suggest an evolutionary scenario for the dwarf galaxies in our Local Group, faint dSphs being the descendents of LSBs and brighter dSphs/dEs being the final state of HSB satellites.

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 analyse parallel N-body simulations of three Cold Dark Matter (CDM) universes to study the abundance and clustering of galaxy clusters. The simulations cover a volume comparable to the forthcoming SDSS. We are able to make robust measurements of cluster properties to a redshift larger than unity. We extract halos using two independent, public domain group finders (FOF & HOP) and find consistent results. The correlation function of clusters is in very good agreement with a simple analytic prescription based upon a Lagrangian biasing scheme developed by Mo & White (1996) and the Press-Schechter (PS) formalism for the mass function. The R_0--D_c relation for the open CDM model is in good agreement with the results from the APM Cluster Survey. The SCDM universe shows a robust deviation in the shape and evolution of the mass function when compared with that predicted by the PS formalism. Critical models with a low sigma_8 normalization or small shape parameter Gamma show an excess of massive clusters compared with the PS prediction. When cluster normalized, the SCDM universe at z =1 contains 10 times more clusters with temperatures greater than 7keV, compared with the PS prediction. The agreement between the analytic and N-body mass functions of SCDM can be improved if the value of the delta_c (the extrapolated linear theory threshold for collapse) is revised to be $ delta_c(z) = 1.685[(0.7/sigma_8)(1+z)]^{-0.125}. Our best estimate for the amplitude of fluctuations inferred from the local cluster abundance for SCDM is sigma_{8} = 0.5 \pm 0.04. However, the discrepancy between the temperature function predicted in a critical density universe and that observed at z=0.33 (Henry et al. 1998) remains. (abridged)

The APM Galaxy Survey maps the angular positions of more than one million galaxies to $b_{J} = 20$, covering a volume comparable to the forthcoming Sloan Digital Sky Survey. A numerical algorithm has been developed to estimate the power spectrum in three dimensions, using the angular clustering and a model for the redshift distribution of APM Galaxies. The power spectrum obtained is free from distortions of the pattern of clustering caused by the peculiar motions of galaxies. We discuss the uncertainties in the estimated power spectrum and describe tests of the algorithm using large numerical simulations. The APM Galaxy power spectrum shows an inflection at a wavenumber $k \sim 0.15 h/Mpc$, with evidence for a peak or turnover in the range $ k \sim 0.03-0.06 h/Mpc$. These features can place strong constraints on Cold Dark Matter models for structure formation.

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.

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 K-band Hubble diagram for a sample of brightest cluster galaxies (BCGs) in the redshift range 0<z<1 shows a very small scatter (0.3 magnitudes r.m.s). The BCGs exhibit very little luminosity evolution in this redshift range: if q_0=0.0 we detect no luminosity evolution; for q_0=0.5 we measure a small negative evolution. If the mass in stars of these galaxies had remained constant over this period of time, substantial positive luminosity evolution would be expected: BCGs should have been brighter in the past since their stars were younger. This suggests that the stellar mass of the BCGs has been assembled over time through merging and accretion. We estimate that the stellar mass in a typical BCG has grown by a factor ~ 2 since z=1 if q_0=0.0 or by factor ~ 4 if q_0=0.5. These results are in remarkably good agreement with the predictions of semi-analytic models of galaxy formation and evolution set in the context of a hierarchical scenario for structure formation.

We analyse the K-band Hubble diagram for a sample of brightest cluster galaxies (BCGs) in the redshift range 0<z<1. In good agreement with earlier studies, we confirm that the scatter in the absolute magnitudes of the galaxies is small (0.3 magnitudes). The BCGs exhibit very little luminosity evolution in this redshift range: if q_0=0.0 we detect no luminosity evolution; for q_0=0.5 we measure a small negative evolution (i.e., BCGs were about 0.5 magnitudes fainter at z=1 than today). If the mass in stars of these galaxies had remained constant over this period of time, substantial positive luminosity evolution would be expected: BCGs should have been brighter in the past since their stars were younger. A likely explanation for the observed zero or negative evolution is that the stellar mass of the BCGs has been assembled over time through merging and accretion, as expected in hierarchical models of galaxy formation. The colour evolution of the BCGs is consistent with that of an old stellar population (z_form >2) that is evolving passively. We can thus use evolutionary population synthesis models to estimate the rate of growth in stellar mass for these systems. We find that the stellar mass in a typical BCG has grown by a factor ~ 2 since z=1 if q_0=0.0 or by factor ~ 4 if q_0=0.5. These results are in good agreement with the predictions of semi-analytic models of galaxy formation and evolution set in the context of a hierarchical scenario for structure formation. The models predict a scatter in the luminosities of the BCGs that is somewhat larger than the observed one, but that depends on the criterion used to select the model clusters.

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 produce mock angular catalogues from simulations with different initial power spectra to test methods that recover measures of clustering in three dimensions, such as the power spectrum, variance and higher order cumulants. We find that the statistical properties derived from the angular mock catalogues are in good agreement with the intrinsic clustering in the simulations. In particular, we concentrate on the detailed predictions for the shape of the power spectrum, $P(k)$. We find that there is good evidence for a break in the galaxy $P(k)$ at scales between $ 0.02 < k < 0.06 ~\impc$ using an inversion technique applied to the angular correlation function measured from the APM Galaxy Survey. For variants on the standard Cold Dark Matter model, a fit at the location of the break implies $\Omega h= 0.45 \pm 0.10$, where $\Omega$ is the ratio of the total matter density to the critical density and Hubble's constant is parameterised as $H_{0}= 100~hkm ~s^{-1}Mpc^{-1}$. On slightly smaller, though still quasi-linear scales, there is a feature in the APM power spectrum where the local slope changes appreciably, with the best match to CDM models obtained for $\Omega h \simeq 0.2$. Hence the location and narrowness of the break in the APM power spectrum combined with the rapid change in its slope on quasi-linear scales cannot be matched by any variant of CDM, including models that have a non-zero cosmological constant or a tilt to the slope of the primordial $P(k)$. These results are independent of the overall normalization of the CDM models or any simple bias that exists between the galaxy and mass distributions.

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.

We present the first measurement of the angular correlation function in a $K$-selected galaxy survey, from two fields covering almost 10 square degrees. The angular correlation function at small angles is well described by a $\theta^{-0.8}$ power law, as for optically selected samples. The clustering amplitude is reduced by a factor of $\sim 4$ between galaxies with $K<15$ and a fainter magnitude slice with $15<K<16$. This allows us to place constraints upon the redshift distribution of the galaxies and their spatial correlation function. We find no clear evidence for a change in clustering amplitude when galaxies are selected by their observed $B-K$ colours.

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.

The null hypothesis that the three dimensional power spectrum measured from the APM Survey is consistent with the one dimensional power spectrum measured from the pencil beams surveys of Broadhurst et al and of Szalay et al is tested. The external estimates of the mean power that we make are sensitive to details of the model for the survey geometry and to the assumed level of the distortion of the pattern of galaxy clustering caused by peculiar motions. We find that the measured 3D clustering of galaxies can account for the presence of peaks in the one dimensional power spectrum, but is less successful in recovering the detailed appearence of the observations. We find no strong evidence for any additional large scale structure in the deep pencil beams beyond that recovered from the APM Survey. We conclude therefore that it is unlikey that large scale structure can be responsible for the steep local number counts of bright galaxies.

We compare the accuracy of published formulae that transform the linear perturbation theory power spectrum into the nonlinear regime against the results of an ensemble of large N-body simulations. We find that the modified transformation given by Jain, Mo and White (1995) performs well for initial fluctuations that have a power law power spectrum, but is less succesful for scale dependent spectra. The inverse transformation to obtain the linear power spectrum is applied to the APM Galaxy Survey power spectrum measured by Baugh \& Efstathiou (1993). The resulting linear spectrum is used to generate the initial density fluctuations in an N-body simulation, which is evolved to match the measured APM amplitude on large scales. We find very good agreement between the final power spectrum of the simulation and the measured APM power spectrum. However, the higher moments for the particle distribution only match the ones recovered from the APM Survey on large scales, $R \simgt 10 \mpc$. On small scales, $R \simlt 10 \mpc$, the APM estimations give smaller amplitudes, indicating that non-gravitational effects, such as biasing, are important on those scales. Our approach can be used to constrain a model of how light from galaxies traces the underlying mass distribution.

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.

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 present a determination of the real space correlation function for galaxies in the APM Survey, with magnitudes in the range $17 < b_J < 20$. We recover a power law form for the correlation function $\xi(r) = (r/4.5)^{-1.7}$ on scales $r < 4 Mpc/h$, for a flat universe and clustering that is fixed in comoving coordinates.If we assume that clustering evolves according to linear perturbation theory, we obtain $\xi(r) = (r/5.25)^{-1.7}$. We recover a shoulder feature in $\xi(r)$ for scales $4 < r < 25 Mpc/h$, with the correlation function rising above the power law, before becoming consistent with zero around $r > 40 Mpc/h$. The form we obtain for $\xi(r)$ is smooth on large scales, allowing a determination of the distortions in the correlation function in redshift space, caused by the peculiar motions of galaxies to be made. Comparing the real space correlation function with the measurement of the redshift space correlation function by Loveday et al (1995), we find that a value for $omega^{0.6}/b$ (where $\omega$ is the density parameter of the universe and b is the bias factor between the density of light and mass) of unity is ruled out at the 2 sigma level.

We test the hypothesis that elliptical galaxies are formed by violent mergers in a universe with hierarchical structure formation. Within the framework of a semi-analytic scheme for galaxy formation, we predict the distribution of morphological types with redshift and compare the colours of elliptical galaxies in different environments.

We present a comparison of the predictions of perturbation theory for the hierarchical J-order correlation amplitudes $S_J=\xibar_J/\xibar_2^{J-1}$, $J=3-10$, with the results of large numerical simulations of gravitational clustering. We consider two different initial power spectra of density fluctuations, one being flatter than the other. The clustering amplitudes $S_J$ measured in these models are different at all scales. In each case, the perturbation theory predictions give an excellent agreement with the simulations on scales for which the variance is approximately linear, i.e. for scales on which $\xibar_2 \simlt 1$. For cells of radius $R \simgt 30 \Mpc$, the sampling variance arising from the finite size of our simulations dominates the results. We also find that $S_J$ are roughly independent of time, $\Omega$, or $\lambda$, as expected. On comparing with the observations, one can use these results to discriminate between models for structure formation.

We present a detailed comparison of the predictions of perturbation theory for the averaged J-point correlation functions, $\xibar_J$, with the results of numerical simulations of gravitational clustering. We have carried out a systematic analysis of this method using ensembles of simulations with different numbers of particles, different box sizes and using different particle arrangements and clustering amplitudes in the initial conditions. We estimate $\xibar_J$, for $J=2-10$, from moments of counts-in-cells. We find significant non-linear effects in the variance, $J=2$, even at scales as large as $R \sim 30 \Mpc$. Perturbation theory gives remarkable agreement at large scales, where $\xibar_2 \simlt 1$, with the measured hierarchical amplitudes $S_J=\xibar_J/\xibar_2^{J-1}$. We have followed the evolution of $\xibar_J$ in time and find that for a change in the effective redshift of at least $\Delta z \simeq 2$ the amplitudes $S_J$ remain unchanged, despite the fact that the $\xibar_J$ have evolved by large factors, $\simeq 10^{J-1}$. We illustrate how these results can be applied to interpret the clustering in galaxy surveys and conclude that the observed hierarchical pattern in the APM is compatible with gravitational evolution in unbiased, initially Gaussian, models.