Old high-redshift galaxies and primordial density fluctuation spectra

We have discovered a population of extremely red galaxies at z  ≃ 1.5 which have apparent stellar ages of ≳ 3 Gyr, based on detailed spectroscopy in the rest-frame ultraviolet. In order for galaxies to have existed at the high collapse redshifts indicated by these ages, there must be a minimum level of power in the density fluctuation spectrum on galaxy scales. This paper compares the required power with that inferred from other high-redshift populations: damped Lyα absorbers and Lyman-limit galaxies at z  ≃ 3.2. If the collapse redshifts for the old red galaxies are in the range z _c ≃ 6–8, there is general agreement between the various tracers on the required inhomogeneity on 1-Mpc scales. This level of small-scale power requires the Lyman-limit galaxies to be approximately ν ≃ 3.0 fluctuations, implying a very large bias parameter b  ≃ 6. If the collapse redshifts of the red galaxies are indeed in the range z _c = 6–8 required for power spectrum consistency, their implied ages at z  ≃ 1.5 are between 3 and 3.8 Gyr for essentially any model universe of current age 14 Gyr. The age of these objects as deduced from gravitational collapse thus provides independent support for the ages estimated from their stellar populations. Such early-forming galaxies are rare, and their contribution to the cosmological stellar density is consistent with an extrapolation to higher redshifts of the star formation rate measured at z  < 5; there is no evidence for a general era of spheroid formation at extreme redshifts.     

Spin temperatures and covering factors for H i 21-cm absorption in damped Lyman α systems

We investigate the practice of assigning high spin temperatures to damped Lyman α absorption systems (DLAs) not detected in H  i 21-cm absorption. In particular, Kanekar & Chengalur have attributed the mix of 21-cm detections and non-detections in low-redshift  ( z _abs≤ 2.04) DLAs  to a mix of spin temperatures, while the non-detections at high redshift were attributed to high spin temperatures. Below   z _abs= 0.9  , where some of the DLA host galaxy morphologies are known, we find that 21-cm absorption is normally detected towards large radio sources when the absorber is known to be associated with a large intermediate (spiral) galaxy. Furthermore, at these redshifts, only one of the six 21-cm non-detections has an optical identification and these DLAs tend to lie along the sight-lines to the largest background radio continuum sources. For these and many of the high-redshift DLAs occulting large radio continua, we therefore expect covering factors of less than the assumed/estimated value of unity. This would have the effect of introducing a range of spin temperatures considerably narrower than the current range of  Δ T _s≳ 9000 K  , while still supporting the hypothesis that the high-redshift DLA sample comprises a larger proportion of compact galaxies than the low-redshift sample.     

The changing interstellar medium of massive elliptical galaxies and cosmic evolution of radio galaxies and quasars

The recently discovered apparent dramatic expansion in the effective radii of massive elliptical galaxies from   z ≃ 2  to ≃0.1 has been interpreted in terms of either galaxy mergers or the rapid loss of cold gas due to active galactic nuclei (AGN) feedback. In examining the latter case, we have quantified the extent of the expansion, which is uncertain observationally, in terms of the star formation parameters and time of the expulsion of the cold gas. In either case, the large global decrease in stellar density should translate into a major drop in the interstellar medium density and pressure with cosmic epoch. These cosmological changes are expected to have a major influence on the gas accretion mode, which will shift from 'cold' thin disc accretion at high redshifts towards 'hot' Bondi fed Advection Dominated Accretion Flow (ADAF) accretion at low redshifts. The decline of angular momentum inflow would then lead to a spin down of the black hole, for which we have calculated more precise time-scales; a value of about 0.2 Gyr is typical for a  10⁹ M_⊙  central black hole. These results have implications for the different cosmological evolutionary patterns found for the luminosity functions of powerful and weak radio galaxies.     

The SCUBA Half Degree Extragalactic Survey (SHADES) – IX. The environment, mass and redshift dependence of star formation

We present a comparison between the SCUBA (Submillimetre Common User Bolometer Array) Half Degree Extragalactic Survey (SHADES) at 450 and  850 μm  in the Lockman Hole East with a deep Spitzer Space Telescope survey at  3.6–24 μm  conducted in guaranteed time. Using stacking analyses we demonstrate a striking correspondence between the galaxies contributing the submm extragalactic background light, with those likely to dominate the backgrounds at Spitzer wavelengths. Using a combination BRIzK plus Spitzer photometric redshifts, we show that at least a third of the Spitzer -identified submm galaxies at  1 < z < 1.5  appear to reside in overdensities when the density field is smoothed at 0.5–2 Mpc comoving diameters, supporting the high-redshift reversal of the local star formation–galaxy density relation. We derive the dust-shrouded cosmic star formation history of galaxies as a function of assembled stellar masses. For model stellar masses  <10¹¹ M_⊙  , this peaks at lower redshifts than the ostensible   z ∼ 2.2  maximum for submm point sources, adding to the growing consensus for 'downsizing' in star formation. Our surveys are also consistent with 'downsizing' in mass assembly. Both the mean star formation rates  〈d M _*/d t 〉  and specific star formation rates  〈(1/ M _*) d M _*/d t 〉  are in striking disagreement with some semi-analytic predictions from the Millenium Simulation. The discrepancy could either be resolved with a top-heavy initial mass function, or a significant component of the submm flux heated by the interstellar radiation field.     

On the reliability of C iv λ1549 as an abundance indicator for high-redshift star-forming galaxies

We reconsider the use of the equivalent width of C  iv  λ1549, EW(C  iv ), as an indicator of the oxygen abundance in star-forming galaxies, as proposed by Heckman et al. for nearby starbursts. We refine the local calibration of EW(C  iv ) versus log (O/H) by using a restricted wavelength window which minimizes blending with interstellar absorption lines. When applied to the stellar component only of the complex C  iv  λ1549 features in two high-redshift galaxies with good quality spectra,  MS 1512−cB58 ( z = 2.7268)  and  Q1307−BM1163 ( z = 1.4105)  , the local calibration gives values of the oxygen abundance which are in good agreement with other metallicity determinations based on nebular emission and interstellar absorption lines. Our main conclusion is that for this method to give reliable results at high redshifts, it should only be used on data of sufficiently high spectral resolution  ( R ≳ 1000)  for stellar and interstellar C  iv components to be clearly separated. Oxygen abundances will be systematically overestimated if the local calibration is applied to spectra of high-redshift galaxies obtained with the low resolving powers  ( R ≃ 200–300)  of many current wide-field surveys. It will also be necessary to understand better the causes of the scatter in the local relation, before we can be confident of inferences from it at high z .     

The evolution of 3CR radio galaxies from z1

We present the results of a comprehensive re-analysis of the images of a virtually complete sample of 28 powerful 3CR radio galaxies with redshifts 0.6< z <1.8 from the Hubble Space Telescope ( HST ) archive. Using a two-dimensional modelling technique we have derived scalelengths and absolute magnitudes for a total of 16 3CR galaxies with a median redshift of z =0.8. Our results confirm the basic conclusions of Best, Longair & Röttgering in that we also find z =1 3CR galaxies to be massive, well-evolved ellipticals, the infrared emission of which is dominated by starlight. However, we in fact find that the scalelength distribution of 3CR galaxies at z ≃1 is completely indistinguishable from that derived for their low-redshift counterparts from our own recently completed HST study of active galactic nuclei hosts at z ≃0.2. There is thus no evidence that 3CR radio galaxies at z ≃1 are dynamically different from 3CR galaxies at low redshift. Moreover, for a 10-object subsample we have determined the galaxy parameters with sufficient accuracy to demonstrate, for the first time, that the z ≃1 3CR galaxies follow a Kormendy relation that is indistinguishable from that displayed by low-redshift ellipticals if one allows for purely passive evolution. The implied rather modest level of passive evolution since z ≃1 is consistent with that predicted from spectrophotometric models provided one assumes a high formation redshift ( z ≥4) within a low-density universe. We conclude that there is no convincing evidence for significant dynamical evolution among 3CR galaxies in the redshift interval 0< z <1, and that simple passive evolution remains an acceptable interpretation of the K – z relation for powerful radio galaxies.     

The impact of H i in galaxies on 21-cm intensity fluctuations during the reionization epoch

We investigate the impact of neutral hydrogen (H  i ) in galaxies on the statistics of 21-cm fluctuations using seminumerical modelling. Following the reionization of hydrogen, the H  i content of the Universe is dominated by damped absorption systems (DLAs), with a cosmic density in H  i that is observed to be constant at a level equal to ∼2 per cent of the cosmic baryon density from   z ∼ 1  to   z ∼ 5  . We show that extrapolation of this constant fraction into the reionization epoch results in a reduction in the amplitude of 21-cm fluctuations over a range of spatial scales. We further find that consideration of H  i in galaxies/DLAs reduces the prominence of the H  ii region induced shoulder in the 21-cm power spectrum (PS), and hence modifies the scale dependence of 21-cm fluctuations. We also estimate the 21-cm–galaxy cross PS and show that the cross PS changes sign on scales corresponding to the H  ii regions. From consideration of the sensitivity for forthcoming low-frequency arrays, we find that the effects of H  i in galaxies/DLAs on the statistics of 21-cm fluctuations will be significant with respect to the precision of a PS or cross PS measurement. In addition, since overdense regions are reionized first we demonstrate that the cross-correlation between galaxies and 21-cm emission changes sign at the end of the reionization era, providing an alternative avenue to pinpoint the end of reionization. The sum of our analysis indicates that the H  i content of the galaxies that reionize the universe will need to be considered in detailed modelling of the 21-cm intensity PS in order to correctly interpret measurements from forthcoming low-frequency arrays.     

The chemical evolution of the Universe – I. High column density absorbers

We construct a simple, robust model of the chemical evolution of galaxies from high to low redshift, and apply it to published observations of damped Lyman α quasar absorption line systems (DLAs). The elementary model assumes quiescent star formation and isolated galaxies (no interactions, mergers or gas flows). We consider the influence of dust and chemical gradients in the galaxies, and hence explore the selection effects in quasar surveys. We fit individual DLA systems to predict some observable properties of the absorbing galaxies, and also indicate the expected redshift behaviour of chemical element ratios involving nucleosynthetic time delays.
Despite its simplicity, our 'monolithic collapse' model gives a good account of the distribution and evolution of the metallicity and column density of DLAs, and of the evolution of the global star formation rate and gas density below redshifts z ∼3. However, from the comparison of DLA observations with our model, it is clear that star formation rates at higher redshifts ( z >3) are enhanced. Galaxy interactions and mergers, and gas flows very probably play a major role.     

The escape of ionizing photons from supernova-dominated primordial galaxies

In order to assess the contribution of Lyman break galaxies (LBGs) and Lyman α emitters (LAEs) at redshifts  3 < z < 7  to the ionization of intergalactic medium (IGM), we investigate the escape fractions of ionizing photons from supernova-dominated primordial galaxies by solving the three-dimensional (3D) radiative transfer. The model galaxy is employed from an ultra-high-resolution chemodynamic simulation of a primordial galaxy by Mori & Umemura, which well reproduces the observed properties of LAEs and LBGs. The total mass of model galaxy is  10¹¹ M_⊙  . We solve not only photoionization but also collisional ionization by shocks. In addition, according to the chemical enrichment, we incorporate the effect of dust extinction, taking the size distributions of dust into account. As a result, we find that dust extinction reduces the escape fractions by a factor of 1.5–8.5 in the LAE phase and by a factor of 2.5–11 in the LBG phase, while the collisional ionization by shocks increases the escape fractions by a factor of  ≈2  . The resultant escape fractions are 0.07–0.47 in the LAE phase and 0.06–0.17 in the LBG phase. These results are well concordant with the recent estimations derived from the flux density ratio at 1500 to 900 Å of LAEs and LBGs. Combining the resultant escape fractions with the luminosity functions of LAEs and LBGs, we find that high- z LAEs and LBGs can ionize the IGM at   z = 3–5  . However, ionizing radiation from LAEs as well as LBGs falls short of ionizing the IGM at   z > 6  . That implies that additional ionization sources may be required at   z > 6  .     

Understanding the redshift evolution of the luminosity functions of Lyman α emitters

We present a semi-analytical model of star formation which explains simultaneously the observed ultraviolet (UV) luminosity function (LF) of high-redshift Lyman break galaxies (LBGs) and LFs of Lyman α emitters. We consider both models that use the Press–Schechter (PS) and Sheth–Tormen (ST) halo mass functions to calculate the abundances of dark matter haloes. The Lyman α LFs at   z ≲ 4  are well reproduced with only ≲10 per cent of the LBGs emitting Lyman α lines with rest equivalent width greater than the limiting equivalent width of the narrow band surveys. However, the observed LF at   z > 5  can be reproduced only when we assume that nearly all LBGs are Lyman α emitters. Thus, it appears that  4 < z < 5  marks the epoch when a clear change occurs in the physical properties of the high-redshift galaxies. As Lyman α escape depends on dust and gas kinematics of the interstellar medium (ISM), this could mean that on an average the ISM at   z > 5  could be less dusty, more clumpy and having more complex velocity field. All of these will enable easier escape of the Lyman α photons. At   z > 5  , the observed Lyman α LF are well reproduced with the evolution in the halo mass function along with very minor evolution in the physical properties of high-redshift galaxies. In particular, up to   z = 6.5  , we do not see the effect of evolving intergalactic medium opacity on the Lyman α escape from these galaxies.     

The H i opacity of the intergalactic medium at redshifts 1.6 < z < 3.2

We use high-quality echelle spectra of 24 quasi-stellar objects to provide a calibrated measurement of the total amount of Lyα forest absorption (DA) over the redshift range  2.2 < z < 3.2  . Our measurement of DA excludes absorption from metal lines or the Lyα lines of Lyman-limit systems and damped Lyα systems. We use artificial spectra with realistic flux calibration errors to show that we are able to place continuum levels that are accurate to better than 1 per cent. When we combine our results with our previous results between  1.6 < z < 2.2  , we find that the redshift evolution of DA is well described over  1.6 < z < 3.2  as   A (1 + z )^γ  , where   A = 0.0062  and  γ= 2.75  . We detect no significant deviations from a smooth power-law evolution over the redshift range studied. We find less H  i absorption than expected at   z = 3  , implying that the ultraviolet background is  ∼40  per cent higher than expected. Our data appears to be consistent with an H  i ionization rate of  Γ∼ 1.4 × 10⁻¹² s⁻¹  .     

Faint Lyman-break galaxies as a crucial test for galaxy formation models

It has recently been shown that galaxy formation models within the Λ cold dark matter cosmology predict that, compared to the observed population, small galaxies (with stellar masses  <10¹¹ M_⊙  ) form too early, are too passive since   z ∼ 3  and host too old stellar populations at   z = 0  . We then expect an overproduction of small galaxies at   z ≳ 4  that should be visible as an excess of faint Lyman-break galaxies. To check whether this excess is present, we use the morgana galaxy formation model and grasil spectrophotometric  +  radiative transfer code to generate mock catalogues of deep fields observed with Hubble Space Telescope Advanced Camera for Surveys. We add observational noise and the effect of Lyman α emission, and perform colour–colour selections to identify Lyman-break galaxies. The resulting mock candidates have plausible properties that closely resemble those of observed galaxies. We are able to reproduce the evolution of the bright tail of the luminosity function of Lyman-break galaxies (with a possible underestimate of the number of the brightest i -dropouts), but uncertainties and degeneracies in dust absorption parameters do not allow to give strong constraints to the model. Besides, our model shows a clear excess with respect to observations of faint Lyman-break galaxies, especially of   z ₈₅₀∼ 27 V   -dropouts at   z ∼ 5  . We quantify the properties of these 'excess' galaxies and discuss the implications: these galaxies are hosted in dark matter haloes with circular velocities in excess of 100 km s⁻¹, and their suppression may require a deep rethinking of stellar feedback processes taking place in galaxy formation.     

Hubble Space Telescope imaging survey of sub-mJy star-forming galaxies – I. Morphologies at z ∼ 0.2

We present the first results of our Hubble Space Telescope HST WFPC2 F814W snapshot imaging survey, targeting virtually all sub-mJy decimetric radio-selected star-forming galaxies. The radio selection at ∼1 GHz is free from extinction effects and the radio luminosities are largely unaffected by AGN contamination, making these galaxies ideal tracers of the cosmic star formation history. A subsample of four targets is presented here, selected at 1.4 GHz from the spectroscopically homogenous and complete samples of Benn et al. and Hopkins et al. The redshifts are confined to a narrow range around z ∼0.2, to avoid differential evolution, with a radio luminosity close to L ∗ where the galaxies dominate the comoving volume-averaged star formation rate. We find clearly disturbed morphologies resembling those of ultraluminous infrared galaxies, indicating that galaxy interactions may be the dominant mechanism for triggering star formation at these epochs. The morphologies are also clearly different from those of coeval quasars and radio galaxies, as found in star-forming galaxies selected at other wavelengths. This may prove challenging for models that propose direct causal links between AGN evolution and the cosmic star formation history at these epochs. The asymmetries are typically much larger than seen in the Canada–France Redshift Survey at similar redshifts, optical luminosities and H α -derived star formation rates, indicating the possible existence of an obscuration-related morphological bias in such samples.     

A method to measure the mass of damped Lyα absorber host galaxies using fluctuations in 21-cm emission

Observations of damped Lyα absorbers (DLAs) indicate that the fraction of hydrogen in its neutral form (H  i ) is significant by mass at all redshifts. This gas represents the reservoir of material that is available for star formation at late times. As a result, observational identification of the systems in which this neutral hydrogen resides is an important missing ingredient in models of galaxy formation. Precise identification of DLA host mass via traditional clustering studies is not practical owing to the small numbers of known systems being spread across sparsely distributed sightlines. However, following the completion of re-ionization, 21-cm surface brightness fluctuations will be dominated by neutral hydrogen in DLAs. No individual DLAs could be detected in 21-cm emission. Rather, observations of these fluctuations will measure the combined clustering signal from all DLAs within a large volume. We show that measurement of the spherically averaged power spectrum of 21-cm intensity fluctuations due to DLAs could be used to measure the galaxy bias for DLA host galaxies when combined with an independent measurement of the cosmological H  i mass density from quasar absorption studies. Utilizing this technique, the low-frequency arrays now under construction could measure the characteristic DLA host mass with a statistical precision as low as 0.3 dex at z ≳ 4. In addition, high signal-to-noise ratio observations of the peculiar-motion-induced anisotropy of the power spectrum would facilitate measurement of both the DLA host mass and the cosmic H  i density directly from 21-cm fluctuations. By exploiting this anisotropy, a second generation of low-frequency arrays with an order of magnitude increase in collecting area could measure the values of cosmic H  i density and DLA host mass, with uncertainties of a few per cent and a few tens of per cent, respectively.     

The galaxy environment of a quasar at z= 1.226: a possible cluster merger

We have conducted ultra-deep optical and deep near-infrared observations of a field around the z =1.226 radio-quiet quasar 104420.8+055739 from the Clowes–Campusano LQG of 18 quasars at z ∼1.3, in search of associated galaxy clustering. Galaxies at these redshifts are distinguished by their extremely red colours, with I − K >3.75, and we find a factor ∼11 overdensity of such galaxies in a 2.25×2.25 arcmin² field centred on the quasar. In particular, we find 15–18 galaxies that have colours consistent with being a population of passively evolving massive ellipticals at the quasar redshift. They form 'fingers' in the V − K K , I − K K colour–magnitude plots at V − K ≃6.9, I − K ≃4.3 comparable to the red sequences observed in other z ≃1.2 clusters. We find suggestive evidence for substructure among the red sequence galaxies in the K image, in the form of two compact groups, 40 arcsec to the north, and 60 arcsec to the south-east of the quasar. An examination of the wider optical images indicates that this substructure is significant, and that the clustering extends to form a large-scale structure 2–3  h ⁻¹ Mpc across. We find evidence for a high (≳50 per cent) fraction of blue galaxies in this system, in the form of 15–20 'red outlier' galaxies with I − K >3.75 and V − I <2.00, which we suggest are dusty, star-forming galaxies at the quasar redshift. Within 30 arcsec of the quasar we find a concentration of blue ( V − I <1) galaxies in a band that bisects the two groups of red sequence galaxies. This band of blue galaxies is presumed to correspond to a region of enhanced star formation. We explain this distribution of galaxies as the early stages of a cluster merger which has triggered both the star formation and the quasar.     

K-band imaging of strong Ca ii-absorber host galaxies at z∼1

We present K -band imaging of fields around 30 strong Ca  ii absorption-line systems, at  0.7 < z < 1.2  , three of which are confirmed damped Lyman α systems. A significant excess of galaxies is found within 6.0 arcsec (≃50 kpc) from the absorber line of sight. The excess galaxies are preferentially luminous compared to the population of field galaxies. A model in which field galaxies possess a luminosity-dependent cross-section for Ca  ii absorption of the form  ( L / L *)^0.7  reproduces the observations well. The luminosity-dependent cross-section for the Ca  ii absorbers appears to be significantly stronger than the established  ( L / L *)^0.4  dependence for Mg  ii absorbers. The associated galaxies lie at large physical distances from the Ca  ii -absorbing gas; we find a mean impact parameter of 24 kpc  ( H ₀= 70 km s⁻¹ Mpc⁻¹)  . Combined with the observed number density of Ca  ii absorbers the large physical separations result in an inferred filling factor of only ∼10 per cent. The physical origin of the strong Ca  ii absorption remains unclear, possible explanations vary from very extended discs of the luminous galaxies to associated dwarf galaxy neighbours, remnants of outflows from the luminous galaxies, or tidal debris from cannibalism of smaller galaxies.     

The first galaxies: assembly, cooling and the onset of turbulence

We investigate the properties of the first galaxies at   z ≳ 10  with highly resolved numerical simulations, starting from cosmological initial conditions and taking into account all relevant primordial chemistry and cooling. A first galaxy is characterized by the onset of atomic hydrogen cooling, once the virial temperature exceeds  ≃10⁴ K  , and its ability to retain photoheated gas. We follow the complex accretion and star formation history of a  ≃5 × 10⁷ M_⊙  system by means of a detailed merger tree and derive an upper limit on the number of Population III (Pop III) stars formed prior to its assembly. We investigate the thermal and chemical evolution of infalling gas and find that partial ionization at temperatures  ≳10⁴ K  catalyses the formation of  H₂  and hydrogen deuteride, allowing the gas to cool to the temperature of the cosmic microwave background. Depending on the strength of radiative and chemical feedback, primordial star formation might be dominated by intermediate-mass Pop III stars formed during the assembly of the first galaxies. Accretion on to the nascent galaxy begins with hot accretion, where gas is accreted directly from the intergalactic medium and shock heated to the virial temperature, but is quickly accompanied by a phase of cold accretion, where the gas cools in filaments before flowing into the parent halo with high velocities. The latter drives supersonic turbulence at the centre of the galaxy and could lead to very efficient chemical mixing. The onset of turbulence in the first galaxies thus likely marks the transition to Pop II star formation.     

An excess of damped Lyman α galaxies near quasi-stellar objects

We present a sample of 33 damped Lyman α systems (DLAs) discovered in the Sloan Digital Sky Survey (SDSS) whose absorption redshifts ( z _abs) are within 6000 km s⁻¹ of the quasi-stellar object's (QSO) systemic redshift ( z _sys). Our sample is based on  731 2.5 < z _sys < 4.5  non-broad absorption line (non-BAL) QSOs from Data Release 3 (DR3) of the SDSS. We estimate that our search is ≈100 per cent complete for absorbers with N (H  i )  ≥ 2 × 10²⁰ cm⁻²  . The derived number density of DLAs per unit redshift, n ( z ), within  Δ v < 6000 km s⁻¹  is higher (3.5σ significance) by almost a factor of 2 than that of intervening absorbers observed in the SDSS DR3, i.e. there is evidence for an overdensity of galaxies near the QSOs. This provides a physical motivation for excluding DLAs at small velocity separations in surveys of intervening 'field' DLAs. In addition, we find that the overdensity of proximate DLAs is independent of the radio-loudness of the QSO, consistent with the environments of radio-loud and radio-quiet QSOs being similar.     

Clustering of galaxies at 3.6 μm in the Spitzer Wide-area Infrared Extragalactic legacy survey

We investigate the clustering of galaxies selected in the 3.6 μm band of the Spitzer Wide-area Infrared Extragalactic (SWIRE) legacy survey. The angular two-point correlation function is calculated for 11 samples with flux limits of S _3.6≥ 4–400 μJy, over an 8 deg² field. The angular clustering strength is measured at >5σ significance at all flux limits, with amplitudes of A = (0.49–29) × 10⁻³ at 1°, for a power-law model, A θ^−0.8. We estimate the redshift distributions of the samples using phenomological models, simulations and photometric redshifts, and so derive the spatial correlation lengths. We compare our results with the Galaxies In Cosmological Simulations (GalICS) models of galaxy evolution and with parametrized models of clustering evolution. The GalICS simulations are consistent with our angular correlation functions, but fail to match the spatial clustering inferred from the phenomological models or the photometric redshifts. We find that the uncertainties in the redshift distributions of our samples dominate the statistical errors in our estimates of the spatial clustering. At low redshifts (median z ≤ 0.5), the comoving correlation length is approximately constant,   r ₀= 6.1 ± 0.5  h ⁻¹  Mpc, and then decreases with increasing redshift to a value of 2.9 ± 0.3  h ⁻¹ Mpc for the faintest sample, for which the median redshift is z ∼ 1. We suggest that this trend can be attributed to a decrease in the average galaxy and halo mass in the fainter flux-limited samples, corresponding to changes in the relative numbers of early- and late-type galaxies. However, we cannot rule out strong evolution of the correlation length over  0.5 < z < 1  .