On the galaxy stellar mass function, the mass–metallicity relation and the implied baryonic mass function

A comparison between published field galaxy stellar mass functions (GSMFs) shows that the cosmic stellar mass density is in the range 4–8 per cent of the baryon density (assuming  Ω_b= 0.045  ). There remain significant sources of uncertainty for the dust correction and underlying stellar mass-to-light ratio even assuming a reasonable universal stellar initial mass function. We determine the   z < 0.05  GSMF using the New York University Value-Added Galaxy Catalog sample of 49 968 galaxies derived from the Sloan Digital Sky Survey and various estimates of stellar mass. The GSMF shows clear evidence for a low-mass upturn and is fitted with a double Schechter function that has  α₂≃−1.6  . At masses below  ∼10^8.5 M_⊙  , the GSMF may be significantly incomplete because of missing low-surface-brightness galaxies. One interpretation of the stellar mass–metallicity relation is that it is primarily caused by a lower fraction of available baryons converted to stars in low-mass galaxies. Using this principle, we determine a simple relationship between baryonic mass and stellar mass and present an 'implied baryonic mass function'. This function has a faint-end slope,  α₂≃−1.9  . Thus, we find evidence that the slope of the low-mass end of the galaxy mass function could plausibly be as steep as the halo mass function. We illustrate the relationship between halo baryonic mass function → galaxy baryonic mass function → GSMF. This demonstrates the requirement for peak galaxy formation efficiency at baryonic masses  ∼10¹¹ M_⊙  corresponding to a minimum in feedback effects. The baryonic-infall efficiency may have levelled off at lower masses.     

The anatomy of the NGC 5044 group – II. Stellar populations and star formation histories

The distribution of galaxy properties in groups and clusters holds important information on galaxy evolution and growth of structure in the Universe. While clusters have received appreciable attention in this regard, the role of groups as fundamental to formation of the present-day galaxy population has remained relatively unaddressed. Here, we present stellar ages, metallicities and α-element abundances derived using Lick indices for 67 spectroscopically confirmed members of the NGC 5044 galaxy group with the aim of shedding light on galaxy evolution in the context of the group environment.
We find that galaxies in the NGC 5044 group show evidence for a strong relationship between stellar mass and metallicity, consistent with their counterparts in both higher and lower mass groups and clusters. Galaxies show no clear trend of age or α-element abundance with mass, but these data form a tight sequence when fitted simultaneously in age, metallicity and stellar mass. In the context of the group environment, our data support the tidal disruption of low-mass galaxies at small group-centric radii, as evident from an apparent lack of galaxies below  ∼10⁹ M_⊙  within ∼100 kpc of the brightest group galaxy. Using a joint analysis of absorption- and emission-line metallicities, we are able to show that the star-forming galaxy population in the NGC 5044 group appears to require gas removal to explain the ∼1.5 dex offset between absorption- and emission-line metallicities observed in some cases. A comparison with other stellar population properties suggests that this gas removal is dominated by galaxy interactions with the hot intragroup medium.     

Further constraints on the evolution of Ks-selected galaxies in the GOODS/CDFS field

We have selected and analysed the properties of a sample of  2905 K_s < 21.5  galaxies in  ∼131 arcmin²  of the Great Observatories Origins Deep Survey (GOODS) Chandra Deep Field South (CDFS), to obtain further constraints on the evolution of K_s -selected galaxies with respect to the results already obtained in previous studies. We made use of the public deep multiwavelength imaging from the optical B through the infrared (IR) 4.5-μm bands, in conjunction with available spectroscopic and COMBO17 data in the CDFS, to construct an optimized redshift catalogue for our galaxy sample. We computed the K_s -band luminosity function and determined that its characteristic magnitude has a substantial brightening and a decreasing total density from   z = 0  to  〈 z 〉= 2.5  . We also analysed the colours and number density evolution of galaxies with different stellar masses. Within our sample, and in contrast to what is observed for less massive systems, the vast majority (∼85–90 per cent) of the most massive  ( M > 2.5 × 10¹¹ M_⊙)  local galaxies appear to be in place before redshift   z ∼ 1  . Around 65–70 per cent of the total assemble between redshifts   z = 1  and 3 and most of them display extremely red colours, suggesting that plausible star formation in these very massive systems should mainly proceed in obscured, short-time-scale bursts. The remaining fraction (up to ∼20 per cent) could be in place at even higher redshifts   z = 3–4  , pushing the first epoch of formation of massive galaxies beyond the limits of current near-IR surveys.     

Baryonic conversion tree: the global assembly of stars and dark matter in galaxies from the Sloan Digital Sky Survey

Using the spectroscopic sample of the Sloan Digital Sky Survey Data Release 1 (SDSS DR1), we measure how gas was transformed into stars as a function of time and stellar mass: the baryonic conversion tree (BCT). There is a clear correlation between early star formation activity and present-day stellar mass: the more massive galaxies have formed approximately 80 per cent of their stars at   z > 1  , while for the less massive ones the value is only approximately 20 per cent. By comparing the BCT with the dark matter merger tree, we find indications that star formation efficiency at   z > 1  had to be approximately a factor of two higher than today (∼10 per cent) in galaxies with present-day stellar mass larger than  2 × 10¹¹ M_⊙  , if this early star formation occurred in the main progenitor. Therefore, the λ cold dark matter (LCDM) paradigm can accommodate a large number of red objects. On the other hand, in galaxies with present-day stellar mass less than  10¹¹ M_⊙  , efficient star formation seems to have been triggered at   z ∼ 0.2  . We show that there is a characteristic mass  ( M _*∼ 10¹⁰ M_⊙)  for feedback efficiency (or lack of star formation). For galaxies with masses lower than this, feedback (or star formation suppression) is very efficient while for higher masses it is not. The BCT, determined here for the first time, should be an important observable with which to confront theoretical models of galaxy formation.     

Strong size evolution of the most massive galaxies since z∼ 2

Using the combined capabilities of the large near-infrared Palomar/DEEP-2 survey, and the superb resolution of the Advanced Camera for Surveys HST camera, we explore the size evolution of 831 very massive galaxies  ( M _⋆≥ 10¹¹ h ⁻²₇₀ M_⊙)  since   z ∼ 2  . We split our sample according to their light concentration using the Sérsic index n . At a given stellar mass, both low  ( n < 2.5)  and high  ( n > 2.5)  concentrated objects were much smaller in the past than their local massive counterparts. This evolution is particularly strong for the highly concentrated (spheroid like) objects. At   z ∼ 1.5  , massive spheroid-like objects were a factor of 4 (±0.4) smaller (i.e. almost two orders of magnitudes denser) than those we see today. These small sized, high-mass galaxies do not exist in the nearby Universe, suggesting that this population merged with other galaxies over several billion years to form the largest galaxies we see today.     

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.     

Galaxy groups at 0.3 ≤z≤ 0.55 – II. Evolution to z∼ 0

We compare deep Magellan spectroscopy of 26 groups at  0.3 ≤ z ≤ 0.55  , selected from the Canadian Network for Observational Cosmology 2 field survey, with a large sample of nearby groups from the 2PIGG catalogue. We find that the fraction of group galaxies with significant [O  ii ]λ3727 emission (≥5 Å) increases strongly with redshift, from ∼29 per cent in 2dFGRS to ∼58 per cent in CNOC2, for all galaxies brighter than  ∼ M _*+ 1.75  . This trend is parallel to the evolution of field galaxies, where the equivalent fraction of emission-line galaxies increases from ∼53 to ∼75 per cent. The fraction of emission-line galaxies in groups is lower than in the field, across the full redshift range, indicating that the history of star formation in groups is influenced by their environment. We show that the evolution required to explain the data is inconsistent with a quiescent model of galaxy evolution; instead, discrete events in which galaxies cease forming stars (truncation events) are required. We constrain the probability of truncation ( P _trunc) and find that a high value is required in a simple evolutionary scenario neglecting galaxy mergers  ( P _trunc≳ 0.3 Gyr⁻¹)  . However, without assuming significant density evolution, P _trunc is not required to be larger in groups than in the field, suggesting that the environmental dependence of star formation was embedded at redshifts   z ≳ 0.45  .     

On the evolution of the black hole: spheroid mass ratio

We present the results of a study which uses the 3C RR sample of radio-loud active galactic nuclei to investigate the evolution of the black hole:spheroid mass ratio in the most massive early-type galaxies from  0 < z < 2  . Radio-loud unification is exploited to obtain virial (linewidth) black hole mass estimates from the 3C RR quasars, and stellar mass estimates from the 3C RR radio galaxies, thereby providing black hole and stellar mass estimates for a single population of early-type galaxies. At low redshift  ( z ≲ 1)  , the 3C RR sample is consistent with a black hole:spheroid mass ratio of   M _bh/ M _sph≃ 0.002  , in good agreement with that observed locally for quiescent galaxies of similar stellar mass  ( M _sph≃ 5 × 10¹¹ M_⊙)  . However, over the redshift interval  0 < z < 2  the 3C RR black hole:spheroid mass ratio is found to evolve as   M _bh/ M _sph∝ (1 + z )^2.07±0.76  , reaching   M _bh/ M _sph≃ 0.008  by redshift   z ≃ 2  . This evolution is found to be inconsistent with the local black hole:spheroid mass ratio remaining constant at a moderately significant level (98 per cent). If confirmed, the detection of evolution in the 3C RR black hole:spheroid mass ratio further strengthens the evidence that, at least for massive early-type galaxies, the growth of the central supermassive black hole may be completed before that of the host spheroid.     

The K-band Hubble diagram for the brightest cluster galaxies: a test of hierarchical galaxy formation models

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

On the formation of massive galaxies: a simultaneous study of number density, size and intrinsic colour evolution in GOODS

The evolution of number density, size and intrinsic colour is determined for a volume-limited sample of visually classified early-type galaxies selected from the Hubble Space Telescope /Advanced Camera for Surveys images of the Great Observatories Origins Deep Survey (GOODS) North and South fields (version 2). The sample comprises 457 galaxies over 320 arcmin² with stellar masses above  3 × 10¹⁰ M_⊙  in the redshift range  0.4 < z < 1.2  . Our data allow a simultaneous study of number density, intrinsic colour distribution and size. We find that the most massive systems  (≳3 × 10¹¹ M_⊙)  do not show any appreciable change in comoving number density or size in our data. Furthermore, when including the results from 2dF galaxy redshift survey, we find that the number density of massive early-type galaxies is consistent with no evolution between   z = 1.2  and 0, i.e. over an epoch spanning more than half of the current age of the Universe. We find large discrepancies between the predictions of semi-analytic models. Massive galaxies show very homogeneous intrinsic colour distributions, with nearly flat radial colour gradients, but with a significant negative correlation between stellar mass and colour gradient, such that red cores appear predominantly in massive galaxies. The distribution of half-light radii – when compared to   z ∼ 0  and   z > 1  samples – is compatible with the predictions of semi-analytic models relating size evolution to the amount of dissipation during major mergers.     

The early-type galaxies NGC 1407 and NGC 1400 – II. Star formation and chemical evolutionary history

We present a possible star formation and chemical evolutionary history for two early-type galaxies NGC 1407 and NGC 1400. They are the two brightest galaxies of the NGC 1407 (or Eridanus-A) group, one of the 60 groups studied as part of the Group Evolution Multi-wavelength Study.
Our analysis is based on new high signal-to-noise ratio spatially resolved integrated spectra obtained at the ESO 3.6-m telescope, out to ∼0.6 (NGC 1407) and ∼1.3 (NGC 1400) effective radii. Using Lick/IDS indices, we estimate luminosity-weighted ages, metallicities and α-element abundance ratios. Colour radial distributions from HST /ACS and Subaru Suprime-Cam multiband wide-field imaging are compared to colours predicted from spectroscopically determined ages and metallicities using single stellar population (SSP) models. The galaxies formed over half of their mass in a single short-lived burst of star formation  (≥100 M_⊙ yr⁻¹)  at redshift z ≥ 5. This likely involved an outside–in mechanism with supernova-driven galactic winds, as suggested by the flatness of the α-element radial profiles and the strong negative metallicity gradients. Our results support the predictions of the revised version of the monolithic collapse model for galaxy formation and evolution. We speculate that, since formation, the galaxies have evolved quiescently and that we are witnessing the first infall of NGC 1400 in the group.     

Black hole demographics from the relation

We analyse a sample of 32 galaxies for which a dynamical estimate of the mass of the hot stellar component, M _bulge, is available. For each of these galaxies, we calculate the mass of the central black hole, M _•, using the tight empirical correlation between M _• and bulge stellar velocity dispersion. The frequency function     is reasonably well described as a Gaussian with     and standard deviation ∼0.45; the implied mean ratio of black hole mass to bulge mass is a factor of ∼5 smaller than generally quoted in the literature. We present marginal evidence for a lower, average black hole mass fraction in more massive galaxies. The total mass density in black holes in the local Universe is estimated to be ∼     consistent with that inferred from high-redshift     active galactic nuclei.     

The galaxy luminosity–size relation and selection biases in the Hubble Ultra Deep Field

We use the Hubble Ultra Deep Field to study the galaxy luminosity–size  ( M – R _e )  distribution. With a careful analysis of selection effects due to both detection completeness and measurement reliability, we identify bias-free regions in the   M – R _e   plane for a series of volume-limited samples. By comparison to a nearby survey also having well-defined selection limits, namely the Millennium Galaxy Catalogue, we present clear evidence for evolution in surface brightness since   z ∼ 0.7  . Specifically, we demonstrate that the mean, rest-frame B -band  〈μ〉 _e   for galaxies in a sample spanning 8 mag in luminosity between   M _B =−22  and −14 mag increases by ∼1.0 mag arcsec⁻² from   z ∼ 0.1  to 0.7. We also highlight the importance of considering surface brightness-dependent measurement biases in addition to incompleteness biases. In particular, the increasing, systematic underestimation of Kron fluxes towards low surface brightnesses may cause diffuse, yet luminous, systems to be mistaken for faint, compact objects.     

Dry mergers: a crucial test for galaxy formation

We investigate the role that dry mergers play in the build-up of massive galaxies within the cold dark matter paradigm. Implementing an empirical shut-off mass scale for star formation, we find a nearly constant dry merger rate of  ∼6 × 10⁻⁵ Mpc⁻³ Gyr⁻¹  at   z ≤ 1  and a steep decline at larger z . Less than half of these mergers are between two galaxies that are morphologically classified as early-types, and the other half is mostly between an early- and late-type galaxy. Latter are prime candidates for the origin of tidal features around red elliptical galaxies. The introduction of a transition mass scale for star formation has a strong impact on the evolution of galaxies, allowing them to grow above a characteristic mass scale of   M _{*, c} ∼ 6.3 × 10¹⁰ M_⊙  by mergers only. As a consequence of this transition, we find that around   M _{*, c}   , the fraction of 1:1 mergers is enhanced with respect to unequal mass major mergers. This suggests that it is possible to detect the existence of a transition mass scale by measuring the relative contribution of equal mass mergers to unequal mass mergers as a function of galaxy mass. The evolution of the high-mass end of the luminosity function is mainly driven by dry mergers at low z . We however find that only 10–20 per cent of galaxies more massive than   M _{*, c}   experience dry major mergers within their last Gyr at any given redshift   z ≤ 1  .     

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.     

The environmental dependence of the chemical properties of star-forming galaxies

We use a  0.040 < z < 0.085  sample of 37 866 star-forming galaxies from the Fourth Data Release of the Sloan Digital Sky Survey to investigate the dependence of gas-phase chemical properties on stellar mass and environment. The local density, determined from the projected distances to the fourth and fifth nearest neighbours, is used as an environment indicator. Considering environments ranging from voids, i.e.  log Σ≲−0.8  , to the periphery of galaxy clusters, i.e.  log Σ≈ 0.8  , we find no dependence of the relationship between galaxy stellar mass and gas-phase oxygen abundance, along with its associated scatter, on local galaxy density. However, the star-forming gas in galaxies shows a marginal increase in the chemical enrichment level at a fixed stellar mass in denser environments. Compared with galaxies of similar stellar mass in low-density environments, they are enhanced by a few per cent for massive galaxies to about 20 per cent for galaxies with stellar masses  ≲10^9.5 M_⊙  . These results imply that the evolution of star-forming galaxies is driven primarily by their intrinsic properties and is largely independent of their environment over a large range of local galaxy density.     

The evolution of galactic discs

We use recent observations of high-redshift galaxies to study the evolution of galactic discs over the redshift range 0 <  z ≲1. The data are inconsistent with models in which discs were already assembled at z  = 1 and have evolved only in luminosity since that time. Assuming that disc properties change with redshift as powers of 1 +   z and analysing the observations assuming an Einstein–de Sitter universe, we find that for given rotation speed, disc scalelength decreases with z as ∼ (1 +  z )⁻¹, total B -band mass-to-light ratio decreases with z as ∼ (1 +  z )⁻¹, and disc luminosity (again in B ) depends only weakly on z . These scalings are consistent with current data on the evolution of disc galaxy abundance as a function of size and luminosity. Both the scalings and the abundance evolution are close to the predictions of hierarchical models for galaxy formation. If different cosmogonies are compared, the observed evolution in disc size and disc abundance favours a flat low-Ω₀ universe over an Einstein–de Sitter universe.     

Feast and Famine: regulation of black hole growth in low-redshift galaxies

We analyse the observed distribution of Eddington ratios  ( L / L _Edd)  as a function of supermassive black hole mass for a large sample of nearby galaxies drawn from the Sloan Digital Sky Survey. We demonstrate that there are two distinct regimes of black hole growth in nearby galaxies. The first is associated with galaxies with significant star formation [   M _*/star formation rate (SFR) ∼  a Hubble time] in their central kiloparsec regions, and is characterized by a broad lognormal distribution of accretion rates peaked at a few per cent of the Eddington limit. In this regime, the Eddington ratio distribution is independent of the mass of the black hole and shows little dependence on the central stellar population of the galaxy. The second regime is associated with galaxies with old central stellar populations (   M _*/SFR ≫  a Hubble time), and is characterized by a power-law distribution function of Eddington ratios. In this regime, the time-averaged mass accretion rate on to black holes is proportional to the mass of stars in the galaxy bulge, with a constant of proportionality that depends on the mean stellar age of the stars. This result is once again independent of black hole mass. We show that both the slope of the power law and the decrease in the accretion rate on to black holes in old galaxies are consistent with population synthesis model predictions of the decline in stellar mass loss rates as a function of mean stellar age. Our results lead to a very simple picture of black hole growth in the local Universe. If the supply of cold gas in a galaxy bulge is plentiful, the black hole regulates its own growth at a rate that does not further depend on the properties of the interstellar medium. Once the gas runs out, black hole growth is regulated by the rate at which evolved stars lose their mass.     

Radio galaxies in the 2SLAQ Luminous Red Galaxy survey – II. The stellar populations of radio-loud and radio-quiet LRGs

We present an analysis of the optical spectra of a volume-limited sample of 375 radio galaxies at redshift  0.4 < z < 0.7  from the 2dF-SDSS (Sloan Digital Sky Survey) Luminous Red Galaxy (LRG) and QSO (quasi-stellar object) (2SLAQ) redshift survey. We investigate the evolution of the stellar populations and emission-line properties of these galaxies. By constructing composite spectra and comparing with a matched sample of radio-quiet sources from the same survey, we also investigate the effect on the galaxy of the presence of an active nucleus.
The composite spectra, binned by redshift and radio luminosity, all require two components to describe them, which we interpret as an old and a younger population. We found no evolution with redshift of the age of the younger population in radio galaxies, nor were they different from the radio-quiet comparison sample. Similarly, there is no correlation with radio power, with the exception that the most powerful radio sources  ( P _1.4 > 10²⁶  W Hz⁻¹) have younger stars and stronger emission lines than the less powerful sources. This suggests that we have located the threshold in radio power where strong emission lines 'switch on', at radio powers of around 10²⁶ W Hz⁻¹. Except for the very powerful radio galaxies, the presence of a currently active radio active galactic nucleus (AGN) does not appear to be correlated with any change in the observed stellar population of a luminous red galaxy at   z ∼ 0.5  .     

The imprint of cosmic reionization on galaxy clustering

We consider the effect of reionization on the clustering properties of galaxy samples at intermediate redshifts ( z ∼ 0.3–5.5). Current models for the reionization of intergalactic hydrogen predict that overdense regions will be reionized early, thus delaying the build-up of stellar mass in the progenitors of massive lower redshift galaxies. As a result, the stellar populations observed in intermediate-redshift galaxies are somewhat younger and hence brighter in overdense regions of the Universe. Galaxy surveys would therefore be sensitive to galaxies with a somewhat lower dark matter mass in overdense regions. The corresponding increase in the observed number density of galaxies can be parametrized as a galaxy bias due to reionization. We model this process using merger trees combined with a stellar synthesis code. Our model demonstrates that reionization has a significant effect on the clustering properties of galaxy samples that are selected based on their star formation properties. The bias correction in Lyman-break galaxies (including those in proposed baryonic oscillation surveys at z < 1) is at the level of 10–20 per cent for a halo mass of  10¹² M_⊙  , leading to corrections factors of 1.5–2 in the halo mass inferred from measurements of clustering length. The reionization of helium could also lead to a sharp increase in the amplitude of the galaxy correlation function at z ∼ 3. We find that the reionization bias is approximately independent of scale and halo mass. However, since the traditional galaxy bias is mass dependent, the reionization bias becomes relatively more important for lower mass systems. The correction to the bias due to reionization is very small in surveys of luminous red galaxies at z < 1.