Search

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., Comment: 10 pages. To appear in ``From Z-Machines to ALMA: (Sub)millimeter Spectroscopy of Galaxies'', ASP Conference Series, eds. A. J. Baker, J. Glenn, A. I. Harris, J. G. Mangum, M. S. Yun

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

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., Comment: 11 pages, 8 colour figures. Accepted for publication in MNRAS. Galaxy catalogues are available for down load from http://www.icc.dur.ac.uk/

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

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., Comment: 6 pages, 3 figures, to appear in proceedings of the conference "The fabulous destiny of galaxies: bridging past and present", 20-24 June 2005, Marseille, France

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., Comment: MNRAS, in press. Minor revision after referee's report. 22 pages, 18 plots. Colour figures for talks (including additional plots which do not appear in the paper) can be downloaded from http://star-www.dur.ac.uk/~cmb/Params/

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.), Comment: Updated to match the in press MNRAS version. P(k) data and covariance matrix available at http://www.mso.anu.edu.au/2dFGRS/Public/Release/PowSpec

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., Comment: 15 pages, 14 figures (1 colour), submitted to MNRAS

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., Comment: 14 pages, 12 figures, accepted MNRAS

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., Comment: Minor revisions to match version published in MNRAS. Colour versions of selected plots are available at http://star-www.dur.ac.uk/~cmb/SCUBAPLOTS/

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., Comment: 24 pages, 13 figures, LaTex, minor changes to match version accepted to MNRAS. Version with high resolution and colour figures available from http://www.ast.cam.ac.uk/~vw/stocbias.html

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), Comment: 19 pages, 17 figures, submitted to MNRAS

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., Comment: 16 pages. Submitted to MNRAS

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., Comment: 15 pages, 8 figures (1 in colour), MNRAS, 348, 866 (as before)

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)., Comment: Same as before

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., Comment: 14 pages, 10 figures. Slight revisions to match version published in MNRAS

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., Comment: 5 pages, 4 figures. Slight revisions to match version published in MNRAS

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., Comment: 10 pages, 6 figures. Dark matter comparison added. Accepted MNRAS

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

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., Comment: 23 pages, replaced with accepted version (MNRAS), minor changes and additional discussion

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., Comment: 10 pages, 7 figures. Submitted to MNRAS

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), Comment: 10 pages, 5 figures. Submitted to ApJ

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

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., Comment: 19 pages, revised following referee's report, and accepted by MNRAS. Higher resolution Figures, an animated version of Figure 12 and a colour version of Figure 22 are available from http://www.nottingham.ac.uk/~ppxeh/

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., Comment: 15 pages, 6 figures, submitted to MNRAS

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^*)., Comment: 20 Pages. To appear in "A New Era in Cosmology" (ASP Conference Proceedings), eds T. Shanks and N. Metcalfe. Replaced with minor corrections to match revised 2dFGRS bias paper

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., Comment: 5 pages, 3 figures, version accepted for publication in PRL. Slightly extended discussion of the choice of priors

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., Comment: To appear in the proceedings `The Mass of Galaxies at Low and High Redshift', Venice 2001, eds. R. Bender, A. Renzini

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)., Comment: 9 pages, 4 figures, MNRAS accepted (minor changes to match the accepted version)

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., Comment: accepted by MNRAS after minor revision. 13 pages, 10 figures

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., Comment: Revised to include additional material, requested by the MNRAS referee, on the photometric calibration of the 2dFGRS and surface brightness dependencies. Full resolution figures and additional material at http://star-www.dur.ac.uk/~cole/2df_lf_selfun

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 11. Our results are consistent with the best-fit results from compilations of cosmic SFR estimates based on UV luminosity density, which give 1.8

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., Comment: 14 pages, 11 figures, submitted to MNRAS

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., Comment: 22 pages, 11 figures, submitted to MNRAS

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., Comment: 14 pages, 13 figures. Accepted for publication in MNRAS

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

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

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

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., Comment: Revised to matched published version: 2001, MNRAS, 328, 64. Colour versions of selected figures for use in talks are available at http://star-www.dur.ac.uk/~irpn/public_2dFGRS/2dF_xir/PaperI/

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., Comment: 10 pages. Invited review to appear in the proceedings of the workshop on "QSO hosts and their environments", IAA, Granada, 10-12 Jan 2001

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., Comment: Nature 410, 169-173 (2001)

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., Comment: 18 pages, 9 figures. Submitted to MNRAS. High-resolution copies of the colour figures can be found at http://www.astro.caltech.edu/~abenson/Mocks/mocks.html

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., Comment: 4 pages with 1 figures in Latex-Kluwer style. To be published in the proceedings of the Granada Euroconference "The Evolution of Galaxies.I-Observational Clues"

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., Comment: 11 pages, Latex. Accepted by MNRAS (revised text includes a clearer explanation of the method used)

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

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., Comment: 31 pages, 17 figures, revised version to appear on 20 October 2000 issue of ApJ

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., Comment: 14 pages, 14 figures, to appear in MNRAS, discussion of number of galaxies per halo added - conclusions unchanged

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., Comment: 6 pages, 3 figures, uses paspconf.sty. To appear in Proceedings of OCIW Workshop (April 1999) on Photometric Redshifts and High Redshift Galaxies

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

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