Richard G. Bower

Breaking the hierarchy of galaxy formation

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 the hierarchical clustering process intrinsic to the cold dark matter (CDM) cosmology. In this paper, we discuss a new implementation of the Durham semi-analytic model of galaxy formation 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 realisations of the galaxy formation process enabled by this simulation result in highly accurate statistics. After adjusting the values of the physical parameters in the model by reference to the properties of the local galaxy population, we investigate the evolution of the K-band luminosity and galaxy stellar mass functions. We 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 � 5 and a star formation rate density which rises at least out to z � 2 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 hierarchical CDM model provides a very good match to these observations.

The EAGLE project: simulating the evolution and assembly of galaxies and their environments

We introduce the Virgo Consortiums EAGLE project, a suite of hydrodynamical simulations that follow the formation of galaxies and black holes in representative volumes. We discuss the limitations of such simulations in light of their finite resolution and poorly constrained subgrid physics, and how these affect their predictive power. One major improvement is our treatment of feedback from massive stars and AGN in which thermal energy is injected into the gas without the need to turn off cooling or hydrodynamical forces, allowing winds to develop without predetermined speed or mass loading factors. Because the feedback efficiencies cannot be predicted from first principles, we calibrate them to the z~0 galaxy stellar mass function and the amplitude of the galaxy-central black hole mass relation, also taking galaxy sizes into account. The observed galaxy mass function is reproduced to ≲0.2 dex over the full mass range, 108<M∗/M⊙≲1011, a level of agreement close to that attained by semi-analytic models, and unprecedented for hydrodynamical simulations. We compare our results to a representative set of low-redshift observables not considered in the calibration, and find good agreement with the observed galaxy specific star formation rates, passive fractions, Tully-Fisher relation, total stellar luminosities of galaxy clusters, and column density distributions of intergalactic CIV and OVI. While the mass-metallicity relations for gas and stars are consistent with observations for M∗≳109M⊙, they are insufficiently steep at lower masses. The gas fractions and temperatures are too high for clusters of galaxies, but for groups these discrepancies can be resolved by adopting a higher heating temperature in the subgrid prescription for AGN feedback. EAGLE constitutes a valuable new resource for studies of galaxy formation.

What Shapes the Luminosity Function of Galaxies

We investigate the physical mechanisms that shape the luminosity function of galaxies in hierarchical clustering models. Beginning with the mass function of dark matter halos in the ΛCDM (Λ cold dark matter) cosmology, we show, in incremental steps, how gas cooling, photoionization at high redshift, feedback processes, galaxy merging, and thermal conduction affect the shape of the luminosity function. We consider three processes whereby supernovae and stellar wind energy can affect the forming galaxy: (1) the reheating of cold disk gas to the halo temperature; (2) expansion of the hot, diffuse halo gas; and (3) complete expulsion of cold disk gas from the halo. We demonstrate that while feedback of form 1 is able to flatten the faint end of the galaxy luminosity function, this process alone does not produce the sharp cutoff observed at large luminosities. Feedback of form 2 is also unable to solve the problem at the bright end of the luminosity function. The relative paucity of very bright galaxies can only be explained if cooling in massive halos is strongly suppressed. This might happen if thermal conduction near the centers of halos is very efficient, or if a substantial amount of gas is expelled from halos by process 3 above. 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 and the associated energy generation play a crucial role in limiting the amount of stars that form in the host galaxy. The models that best reproduce the galaxy luminosity function also give reasonable approximations to the Tully-Fisher relation and the galaxy autocorrelation function.

The 2dF Galaxy Redshift Survey: the environmental dependence of galaxy star formation rates near clusters

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

Galaxy bimodality versus stellar mass and environment

We analyse a z < 0.1 galaxy sample from the Sloan Digital Sky Survey focusing on the variation in the galaxy colour bimodality with stellar mass M and projected neighbour density Σ, and on measurements of the galaxy stellar mass functions. The characteristic mass increases with environmental density from about 10 10. 6 to 10 10.9 M ⊙ (Kroupa initial mass function, H 0 = 70) for Σ in the range 0.1-10 Mpc -2 . The galaxy population naturally divides into a red and blue sequence with the locus of the sequences in colour-mass and colour-concentration indices not varying strongly with environment. The fraction of galaxies on the red sequence is determined in bins of 0.2 in log Σ and log M (12 x 13 bins). The red fraction f r generally increases continuously in both Σ and M such that there is a unified relation: f, = F(Σ, M). Two simple functions are proposed which provide good fits to the data. These data are compared with analogous quantities in semi-analytical models based on the Millennium N-body simulation: the Bower et al. and Croton et al. models that incorporate active galactic nucleus feedback. Both models predict a strong dependence of the red fraction on stellar mass and environment that is qualitatively similar to the observations. However, a quantitative comparison shows that the Bower et al. model is a significantly better match; this appears to be due to the different treatment of feedback in central galaxies.

Gone with the Wind: The Origin of S0 Galaxies in Clusters

We present three-dimensional, high-resolution hydrodynamical simulations of the interaction between the hot ionized intracluster medium and the cold interstellar medium of spiral galaxies. Ram pressure and turbulent/viscous stripping remove 100% of the atomic hydrogen content of luminous galaxies like the Milky Way within 100 million years. These mechanisms naturally account for the morphology of S0 galaxies and the rapid truncation of star formation implied by spectroscopic observations, as well as a host of observational data on the neutral hydrogen (HI) morphology of galaxies in clusters.

Galaxy ecology: groups and low-density environments in the SDSS and 2dFGRS

We analyse the observed correlation between galaxy environment and Halpha emission-line strength, using volume-limited samples and group catalogues of 24 968 galaxies at 0.05 < z < 0.095, drawn from the 2dF Galaxy Redshift Survey (M-bJ < -19.5) and the Sloan Digital Sky Survey (M-r < -20.6). We characterize the environment by: (1) Sigma(5), the surface number density of galaxies determined by the projected distance to the fifth nearest neighbour; and (2) rho(1.1) and rho(5.5), three-dimensional density estimates obtained by convolving the galaxy distribution with Gaussian kernels of dispersion 1.1 and 5.5 Mpc, respectively. We find that star-forming and quiescent galaxies form two distinct populations, as characterized by their H equivalent width, W-0(Halpha). The relative numbers of star-forming and quiescent galaxies vary strongly and continuously with local density. However, the distribution of W-0(Halpha) amongst the star-forming population is independent of environment. The fraction of star-forming galaxies shows strong sensitivity to the density on large scales, rho(5.5), which is likely independent of the trend with local density, rho(1.1). We use two differently selected group catalogues to demonstrate that the correlation with galaxy density is approximately independent of group velocity dispersion, for sigma = 200-1000 km s(-1). Even in the lowest-density environments, no more than similar to70 per cent of galaxies show significant Halpha emission. Based on these results, we conclude that the present-day correlation between star formation rate and environment is a result of short-time-scale mechanisms that take place preferentially at high redshift, such as starbursts induced by galaxy-galaxy interactions.

The Bimodal Galaxy Color Distribution: Dependence on Luminosity and Environment

We analyze the u - r color distribution of 24,346 galaxies with Mr ≤ -18 and z < 0.08, drawn from the Sloan Digital Sky Survey first data release, as a function of luminosity and environment. The color distribution is well fitted with two Gaussian distributions, which we use to divide the sample into a blue and red population. At fixed luminosity, the mean color of the blue (red) distribution is nearly independent of environment, with a weakly significant (~3 σ) detection of a trend for colors to become redder by 0.1-0.14 (0.03-0.06) mag with a factor of ~100 increase in local density, as characterized by the surface density of galaxies within a ±1000 km s-1 redshift slice. In contrast, at fixed luminosity the fraction of galaxies in the red distribution is a strong function of local density, increasing from ~10%-30% of the population in the lowest density environments to ~70% at the highest densities. The strength of this trend is similar for both the brightest (-23 < Mr < -22) and faintest (-19 < Mr < -18) galaxies in our sample. The fraction of red galaxies within the virialized regions of clusters shows no significant dependence on velocity dispersion. Even at the lowest densities explored, a substantial population of red galaxies exists, which might be fossil groups. We propose that most star-forming galaxies today evolve at a rate that is determined primarily by their intrinsic properties and independent of their environment. Any environmentally triggered transformations from blue to red colors must occur either on a short timescale or preferentially at high redshift to preserve the simple Gaussian nature of the color distribution. The mechanism must be effective for both bright and faint galaxies.

Ram pressure stripping the hot gaseous haloes of galaxies in groups and clusters

We use a large suite of carefully controlled full hydrodynamic simulations to study the ram pressure stripping of the hot gaseous haloes of galaxies as they fall into massive groups and clusters. The sensitivity of the results to the orbit, total galaxy mass, and galaxy structural properties is explored. For typical structural and orbital parameters, we find that ∼30 per cent of the initial hot galactic halo gas can remain in place after 10 Gyr. We propose a physically simple analytic model that describes the stripping seen in the simulations remarkably well. The model is analogous to the original formulation of Gunn & Gott, except that it is appropriate for the case of a spherical (hot) gas distribution (as opposed to a face-on cold disc) and takes into account that stripping is not instantaneous but occurs on a characteristic time-scale. The model reproduces the results of the simulations to within ≈10 per cent at almost all times for all the orbits, mass ratios, and galaxy structural properties we have explored. The one exception involves unlikely systems where the orbit of the galaxy is highly non-radial and its mass exceeds about 10 per cent of the group or cluster into which it is falling (in which case the model underpredicts the stripping following pericentric passage). The proposed model has several interesting applications, including modelling the ram pressure stripping of both observed and cosmologically simulated galaxies and as a way to improve present semi-analytic models of galaxy formation. One immediate consequence is that the colours and morphologies of satellite galaxies in groups and clusters will differ significantly from those predicted with the standard assumption of complete stripping of the hot coronae.

The case for AGN feedback in galaxy groups

The relatively recent insight that energy input from supermassive black holes (BHs) can have a substantial effect on the star formation rates (SFRs) of galaxies motivates us to examine the effects of BH feedback on the scale of galaxy groups. At present, groups contain most of the galaxies and a significant fraction of the overall baryon content of the universe and, along with massive clusters, they represent the only systems for which it is possible to measure both the stellar and gaseous baryonic components directly. To explore the effects of BH feedback on groups, we analyse two high resolution cosmological hydrodynamic simulations from the OverWhelmingly Large Simulations (OWLS) project. While both include galactic winds driven by supernovae, only one of the models includes feedback from accreting BHs. We compare the properties of the simulated galaxy groups to a wide range of observational data, including the entropy and temperature profiles of the intragroup medium, hot gas mass fractions, the luminosity temperature and mass temperature scaling relations, the K-band luminosity of the group and its central brightest galaxy (CBG), star formation rates and ages of the CBG, and gas- and stellar-phase metallicities. Both runs yield entropy distributions similar to the data, while the run without AGN feedback yields highly peaked temperature profiles, in discord with the observations. Energy input from supermassive BHs significantly reduces the gas mass fractions of galaxy groups with masses less than a few times 10 14 M� , yielding a gas mass fraction and X-ray luminosity scaling with system temperature that is in excellent agreement with the data, although the detailed scatter in the L T relation is not quite correct. The run without AGN feedback suffers from the well known overcooling problem — the resulting stellar mass fractions are several times larger than observed and present-day cooling flows operate uninhibitedly. By contrast, the run that includes BH feedback yields stellar mass fractions, SFRs, and stellar age distributions in excellent agreement with current estimates, thus resolving the long standing ‘cooling crisis’ of simulations on the scale of groups. Both runs yield very similar gas-phase metal abundance profiles that match X-ray measurements, but they predict very different stellar metallicities. Based on the above, galaxy groups provide a compelling case that feedback from supermassive BHs is a crucial ingredient in the formation of massive galaxies.