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