Testing model predictions of the cold dark matter cosmology for the sizes, colours, morphologies and luminosities of galaxies with the SDSS

The huge size and uniformity of the Sloan Digital Sky Survey (SDSS) make possible an exacting test of current models of galaxy formation. We compare the predictions of the galform semi-analytical galaxy formation model for the luminosities, morphologies, colours and scalelengths of local galaxies. galform models the luminosity and size of the disc and bulge components of a galaxy, and so we can compute quantities which can be compared directly with SDSS observations, such as the Petrosian magnitude and the Sérsic index. We test the predictions of two published models set in the cold dark matter cosmology: the Baugh et al. model, which assumes a top-heavy initial mass function (IMF) in starbursts and superwind feedback, and the Bower et al. model, which uses active galactic nucleus feedback and a standard IMF. The Bower et al. model better reproduces the overall shape of the luminosity function, the morphology–luminosity relation and the colour bimodality observed in the SDSS data, but gives a poor match to the size–luminosity relation. The Baugh et al. model successfully predicts the size–luminosity relation for late-type galaxies. Both models fail to reproduce the sizes of bright early-type galaxies. These problems highlight the need to understand better both the role of feedback processes in determining galaxy sizes, in particular the treatment of the angular momentum of gas reheated by supernovae, and the sizes of the stellar spheroids formed by galaxy mergers and disc instabilities.     

Where do 'red and dead' early-type void galaxies come from?

Void regions of the Universe offer a special environment for studying cosmology and galaxy formation, which may expose weaknesses in our understanding of these phenomena. Although galaxies in voids are observed to be predominately gas rich, star forming and blue, a subpopulation of bright-red void galaxies can also be found, whose star formation was shutdown long ago. Are the same processes that quench star formation in denser regions of the Universe also at work in voids?
We compare the luminosity function of void galaxies in the 2dF Galaxy Redshift Survey, to those from a galaxy formation model built on the Millennium simulation. We show that a global star formation suppression mechanism in the form of low-luminosity 'radio-mode' active galactic nuclei (AGN) heating is sufficient to reproduce the observed population of void early types. Radio-mode heating is environment independent other than its dependence on dark matter halo mass, where, above a critical mass threshold of approximately   M _vir∼ 10^12.5 M_⊙  , gas cooling on to the galaxy is suppressed and star formation subsequently fades. In the Millennium simulation, the void halo mass function is shifted with respect to denser environments, but still maintains a high-mass tail above this critical threshold. In such void haloes, radio-mode heating remains efficient and red galaxies are found; collectively these galaxies match the observed space density without any modification to the model. Consequently, galaxies living in vastly different large-scale environments but hosted by haloes of similar mass are predicted to have similar properties, consistent with observations.     

The dependence of galaxy formation on cosmological parameters: can we distinguish between the WMAP1 and WMAP3 parameter sets?

We combine N -body simulations of structure growth with physical modelling of galaxy evolution to investigate whether the shift in cosmological parameters between the first- and third-year results from the Wilkinson Microwave Anisotropy Probe ( WMAP ) affects predictions for the galaxy population. Structure formation is significantly delayed in the WMAP3 cosmology, because the initial matter fluctuation amplitude is lower on the relevant scales. The decrease in dark matter clustering strength is, however, almost entirely offset by an increase in halo bias, so predictions for galaxy clustering are barely altered. In both cosmologies, several combinations of physical parameters can reproduce observed, low-redshift galaxy properties; the star formation, supernova feedback and active galactic nucleus feedback efficiencies can be played off against each other to give similar results. Models which fit observed luminosity functions predict projected two-point correlation functions which scatter by about 10–20 per cent on large scale and by larger factors on small scale, depending both on cosmology and on details of galaxy formation. Measurements of the pairwise velocity distribution prefer the WMAP1 cosmology, but careful treatment of the systematics is needed. Given present modelling uncertainties, it is not easy to distinguish between the WMAP1 and WMAP3 cosmologies on the basis of low-redshift galaxy properties. Model predictions diverge more dramatically at high redshift. Better observational data at   z > 2  will better constrain galaxy formation and perhaps also cosmological parameters.     

The fraction of galaxies that contain active nuclei and their accretion rates

We investigate the relationship between the present-day optical luminosity function of galaxies and the X-ray luminosity function of Seyfert 1s to determine the fraction of galaxies that host Seyfert 1 nuclei and their Eddington ratios. The local type 1 active galactic nuclei (AGN) X-ray luminosity function is well reproduced if ∼1 per cent of all galaxies are type 1 Seyferts which have Eddington ratios of ∼10⁻³. However, in such a model the X-ray luminosity function is completely dominated by AGN in E and S0 galaxies, contrary to the observed mix of Seyfert host galaxies. To obtain a plausible mix of AGN host galaxy morphologies requires that the most massive black holes in E and S0 galaxies accrete with lower Eddington ratios, or have a lower incidence of Seyfert activity, than the central black holes of later-type galaxies.     

Preliminary results of the multisite time-series campaign on the sdB pulsator PG 1325+101

I present a model for the formation and evolution of a massive disk galaxy, within a growing dark halo whose mass evolves according to cosmological simulations of structure formation. The galactic evolution is simulated with a new 3D chemo-dynamical code, including dark matter, stars and a multi-phase ISM. We follow the evolution from redshift z = 4.85 until the present epoch. The energy release by massive stars and supernovae prevents a rapid collapse of the baryonic matter and delays the maximum star formation until redshift z ≈ 1. The galaxy forms radially from inside-out and vertically from top-to-bottom. The feedback of stars leads to turbulent motions and large-scale flows in the ISM. As one result the galactic disk is significantly enriched by chemical elements synthesized in bulge stars.     

The chemo-dynamical evolution of a disk galaxy

I present a model for the formation and evolution of a massive disk galaxy, within a growing dark halo whose mass evolves according to cosmological simulations of structure formation. The galactic evolution is simulated with a new three-dimensional chemo-dynamical code, including dark matter, stars and a multi-phase ISM. We follow the evolution from redshift z= 4.85 until the present epoch. The energy release by massive stars and supernovae prevents a rapid collapse of the baryonic matter and delays the maximum star formation until redshift z ≈ 1. The galaxy forms radially from inside-out and vertically from top-to-bottom. Correspondingly, the inner halo is the oldest component, followed by the outer halo, the bar/bulge, the thick and the thin disk. The bulge in the model consists of at least two stellar subpopulations, an early collapse population and a population that formed later in the bar. This revised version was published online in August 2006 with corrections to the Cover Date.     

Populating a cluster of galaxies – I. Results at

We simulate the assembly of a massive rich cluster and the formation of its constituent galaxies in a flat, low-density universe. Our most accurate model follows the collapse, the star formation history and the orbital motion of all galaxies more luminous than the Fornax dwarf spheroidal, while dark halo structure is tracked consistently throughout the cluster for all galaxies more luminous than the SMC. Within its virial radius this model contains about     dark matter particles and almost 5000 distinct dynamically resolved galaxies. Simulations of this same cluster at a variety of resolutions allow us to check explicitly for numerical convergence both of the dark matter structures produced by our new parallel N -body and substructure identification codes, and of the galaxy populations produced by the phenomenological models we use to follow cooling, star formation, feedback and stellar aging. This baryonic modelling is tuned so that our simulations reproduce the observed properties of isolated spirals outside clusters. Without further parameter adjustment our simulations then produce a luminosity function, a mass-to-light ratio, luminosity, number and velocity dispersion profiles, and a morphology–radius relation which are similar to those observed in real clusters. In particular, since our simulations follow galaxy merging explicitly, we can demonstrate that it accounts quantitatively for the observed cluster population of bulges and elliptical galaxies.     

Radial distributions of gamma-ray bursts and type Ib/c supernovae in galaxies

Data on the positions of gamma-ray bursts (GRBs) in galaxies are used to construct the radial distributions of their surface density. The gradient in GRB surface density is shown to decrease sharply at a galactocentric distance equal to the effective galactic radius. In central galactic regions, the GRB density distribution agrees with the galactic surface-brightness distribution; in outer regions, the GRB density decreases more slowly than does the surface brightness. Based on improved statistics, we analyze the radial distribution of type Ib/c supernovae. We show that it differs insignificantly from the distributions of other types of supernova and exhibits a much closer similarity to the distribution of star-forming regions than do GRBs. Although the statistics for GRBs is poor, the deviation of their distribution from the distribution of active star-forming regions in nearby galaxies seems to have been firmly established. A correlation of GRBs with the distribution of dark matter in outer galactic regions is not ruled out.     

Galaxy–dark matter correlations applied to galaxy–galaxy lensing: predictions from the semi-analytic galaxy formation models

We use semi-analytic models of galaxy formation combined with high-resolution N -body simulations to make predictions for galaxy–dark matter correlations and apply them to galaxy–galaxy lensing. We analyse cross-power spectra between the dark matter and different galaxy samples selected by luminosity, colour or star formation rate. We compare the predictions with the recent detection by the Sloan Digital Sky Survey (SDSS). We show that the correlation amplitude and the mean tangential shear depend strongly on the luminosity of the sample on scales below 1  h ⁻¹ Mpc, reflecting the correlation between the galaxy luminosity and the halo mass. The cross-correlation cannot, however, be used to infer the halo profile directly because different halo masses dominate on different scales and because not all galaxies are at the centres of the corresponding haloes. We compute the redshift evolution of the cross-correlation amplitude and compare it with those of galaxies and dark matter. We also compute the galaxy–dark matter correlation coefficient and show that it is close to unity on scales above 1  h ⁻¹ Mpc for all considered galaxy types. This would allow one to extract the bias and the dark matter power spectrum on large scales from the galaxy and galaxy–dark matter correlations.     

The X-ray emission in post-merger ellipticals

The evolution in X-ray properties of early-type galaxies is largely unconstrained. In particular, little is known about how, and if, remnants of mergers generate hot gas haloes. Here we examine the relationship between X-ray luminosity and galaxy age for a sample of early-type galaxies. Comparing normalized X-ray luminosity to three different age indicators, we find that L _X L _B increases with age, suggesting an increase in X-ray halo mass with time after the last major star formation episode of a galaxy. The long-term nature of this trend, which appears to continue across the full age range of our sample, poses a challenge for many models of hot halo formation. We conclude that models involving a declining rate of type Ia supernovae, and a transition from outflow to inflow of the gas originally lost by galactic stars, offer the most promising explanation for the observed evolution in X-ray luminosity.     

An indicator of galactic violence

An indicator of galactic violence is proposed, nominally equal to the total galactic mass divided by the rim radius. Four models of galaxy formation are outlined; and for each, the dependence of the rim radius on the galactic mass and the source of angular momentum is examined. If the violence indicator is small, then a normal galaxy results. If the violence indicator is large, then a galaxy with a very massive central black hole is produced (possibly a quasar). A rim radius of zero would indicate the ultimate violence and would lead directly to a non-rotating black hole of galactic mass.Los Alamos National Laboratory is operated by the University of California for the U.S. Dept. of Energy under Contract W-7405-ENG-36.     

Hydrodynamical simulations of cluster formation with central AGN heating

We analyse a hydrodynamical simulation model for the recurrent heating of the central intra-cluster medium (ICM) by active galactic nuclei (AGN). Besides the self-gravity of the dark matter and gas components, our approach includes the radiative cooling and photoheating of the gas, as well as a subresolution multiphase model for star formation and supernova feedback. Additionally, we incorporate a periodic heating mechanism in the form of hot, buoyant bubbles, injected into the intragalactic medium (IGM) during the active phases of the accreting central AGN. We use simulations of isolated cluster haloes of different masses to study the bubble dynamics and heat transport into the IGM. We also apply our model to self-consistent cosmological simulations of the formation of galaxy clusters with a range of masses. Our numerical schemes explore a variety of different assumptions for the spatial configuration of AGN-driven bubbles, for their duty cycles and for the energy injection mechanism, in order to obtain better constraints on the underlying physical picture. We argue that AGN heating can substantially affect the properties of both the stellar and gaseous components of clusters of galaxies. Most importantly, it alters the properties of the central dominant (cD) galaxy by reducing the mass deposition rate of freshly cooled gas out of the ICM, thereby offering an energetically plausible solution to the cooling-flow problem. At the same time, this leads to reduced or eliminated star formation in the central cD galaxy, giving it red stellar colours as observed.     

GABE: Galaxy Assembly with Binary Evolution

We developed a new semi-analytic galaxy formation model: Galaxy Assembly with Binary Evolution(GABE). For the first time, we introduce binary evolution into semi-analytic models of galaxy formation by using the Yunnan-II stellar population synthesis model, which includes various binary interactions. When implementing our galaxy formation model onto the merger trees extracted from the Millennium simulation, it can reproduce a large body of observational results. We find that in the local universe, the model including binary evolution reduces the luminosity at optical and infrared wavelengths slightly, but it increases the luminosity at ultraviolet wavelengths significantly, especially in FUV band. The resulting luminosity function does not change very much over SDSS optical bands and infrared band, but the predicted colors are bluer, especially when the FUV band is under consideration. The new model allows us to explore the physics of various high energy events related to the remnants of binary stars, such as type Ia supernovae, short gamma-ray bursts and gravitational wave events, and their relation with host galaxies in a cosmological context.     

The properties of Lyα emitting galaxies in hierarchical galaxy formation models

We present detailed predictions for the properties of Lyα-emitting galaxies in the framework of the Λ cold dark matter cosmology, calculated using the semi-analytical galaxy formation model galform . We explore a model that assumes a top-heavy initial mass function in starbursts and that has previously been shown to explain the sub-millimetre number counts and the luminosity function of Lyman-break galaxies at high redshift. We show that this model, with the simple assumption that a fixed fraction of Lyα photons escape from each galaxy, is remarkably successful at explaining the observed luminosity function of Lyα emitters (LAEs) over the redshift range  3 < z < 6.6  . We also examine the distribution of Lyα equivalent widths and the broad-band continuum magnitudes of emitters, which are in good agreement with the available observations. We look more deeply into the nature of LAEs, presenting predictions for fundamental properties such as the stellar mass and radius of the emitting galaxy and the mass of the host dark matter halo. The model predicts that the clustering of LAEs at high redshifts should be strongly biased relative to the dark matter, in agreement with observational estimates. We also present predictions for the luminosity function of LAEs at   z > 7  , a redshift range that is starting to be be probed by near-infrared surveys and using new instruments such as the Dark Ages Z Lyman Explorer (DAzLE).     

Broad reprocessed Balmer emission from warped accretion discs

Double-peaked broad emission lines in active galactic nuclei are generally considered to be formed in an accretion disc. In this paper, we compute the profiles of reprocessing emission lines from a relativistic, warped accretion disc around a black hole in order to explore the possibility that certain asymmetries in the double-peaked emission-line profile which cannot be explained by a circular Keplerian disc may be induced by disc warping. The disc warping also provides a solution for the energy budget in the emission-line region because it increases the solid angle of the outer disc portion subtended to the inner portion of the disc. We adopted a parametrized disc geometry and a central point-like source of ionizing radiation to capture the main characteristics of the emission-line profile from such discs. We find that the ratio between the blue and red peaks of the line profiles becoming less than unity can be naturally predicted by a twisted warped disc, and a third peak can be produced in some cases. We show that disc warping can reproduce the main features of multipeaked line profiles of four active galactic nuclei from the Sloan Digital Sky Survey.     

Supersonic motions of galaxies in clusters

We study motions of galaxies in galaxy clusters formed in the concordance Λ cold dark matter cosmology. We use high-resolution cosmological simulations that follow the dynamics of dark matter and gas and include various physical processes critical for galaxy formation: gas cooling, heating and star formation. Analysing the motions of galaxies and the properties of intracluster gas in a sample of eight simulated clusters at z = 0, we study the velocity dispersion profiles of the dark matter, gas and galaxies. We measure the mean velocity of galaxy motions and gas sound speed as a function of radius and calculate the average Mach number of galaxy motions. The simulations show that galaxies, on average, move supersonically with the average Mach number of ≈1.4, approximately independent of the cluster-centric radius. The supersonic motions of galaxies may potentially provide an important source of heating for the intracluster gas by driving weak shocks and via dynamical friction, although these heating processes appear to be inefficient in our simulations. We also find that galaxies move slightly faster than the dark matter particles. The magnitude of the velocity bias,   b _v ≈ 1.1  , is, however, smaller than the bias estimated for subhaloes in dissipationless simulations. Interestingly, we find velocity bias in the tangential component of the velocity dispersion, but not in the radial component. Finally, we find significant random bulk motions of gas. The typical gas velocities are of order ≈20–30 per cent of the gas sound speed. These random motions provide about 10 per cent of the total pressure support in our simulated clusters. The non-thermal pressure support, if neglected, will bias measurements of the total mass in the hydrostatic analyses of the X-ray cluster observations.     

Adiabatic scaling relations of galaxy clusters

The aim of this work is to show that, contrary to popular belief, galaxy clusters are not expected to be self-similar, even when the only energy sources available are gravity and shock-wave heating. In particular, we investigate the scaling relations between mass, luminosity and temperature of galaxy groups and clusters in the absence of radiative processes. Theoretical expectations are derived from a polytropic model of the intracluster medium and compared with the results of high-resolution adiabatic gasdynamical simulations. It is shown that, in addition to the well-known relation between the mass and concentration of the dark matter halo, the effective polytropic index of the gas also varies systematically with cluster mass, and therefore neither the dark matter nor the gas profiles are exactly self-similar. It is remarkable, though, that the effects of concentration and polytropic index tend to cancel each other, leading to scaling relations whose logarithmic slopes roughly match the predictions of the most-basic self-similar models. We provide a phenomenological fit to the relation between polytropic index and concentration, as well as a self-consistent scheme to derive the non-linear scaling relations expected for any cosmology and the best-fitting normalizations of the M – T , L – T and F – T relations appropriate for a Λ cold dark matter universe. The predicted scaling relations reproduce observational data reasonably well for massive clusters, where the effects of cooling and star formation are expected to play a minor role.     

The effects of photoionization on galaxy formation – I. Model and results at z=0

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 haloes. We include in the model the negative feedback on galaxy formation from the photoionizing background. Photoionization inhibits galaxy formation in low-mass dark matter haloes in two ways: (i) heating of the IGM and inhibition of the collapse of gas into dark haloes by the IGM pressure, and (ii) reduction in the rate of radiative cooling of gas within haloes. The result of our method is a self-consistent model of galaxy formation and the IGM. The IGM is reheated twice (during reionization of H  i and He  ii ), 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.