Can the faint submillimetre galaxies be explained in the Λ cold dark matter model?

We present predictions for the abundance of submillimetre galaxies (SMGs) and Lyman-break galaxies (LBGs) in the Λ cold dark matter cosmology. A key feature of our model is the self-consistent calculation of the absorption and emission of radiation by dust. The new model successfully matches the LBG luminosity function, as well as reproducing the properties of the local galaxy population in the optical and infrared. The model can also explain the observed galaxy number counts at 850 μm, but only if we assume a top-heavy initial mass function for the stars formed in bursts. The predicted redshift distribution of SMGs depends relatively little on their flux over the range 1–10 mJy, with a median value of   z ≈ 2.0  at a flux of 5 mJy, in good agreement with the recent measurement by Chapman et al. The counts of SMGs are predicted to be dominated by ongoing starbursts. However, in the model these bursts are responsible for making only a few per cent of the stellar mass locked up in massive ellipticals at the present day.     

Substructures in hydrodynamical cluster simulations

We present predictions for the abundance and nature of extremely red objects (EROs) in the Λ cold dark matter model. EROs are red, massive galaxies observed at   z ≥ 1  and their numbers and properties pose a challenge to hierarchical galaxy formation models. We compare the predictions from two published models, one of which invokes a 'superwind' to regulate star formation in massive haloes and the other which suppresses gas cooling in haloes through 'radio-mode' active galactic nucleus (AGN) feedback. The superwind model underestimates the number counts of EROs by an order of magnitude, whereas the radio-mode AGN feedback model gives excellent agreement with the number counts and redshift distribution of EROs. In the AGN feedback model the ERO population is dominated by old, passively evolving galaxies, whereas observations favour an equal split between old galaxies and dusty starbursts. Also, the model predicts a more extended redshift distribution of passive galaxies than is observed. These comparisons suggest that star formation may be quenched too efficiently in this model.     

Light curves for off-centre ignition models of Type Ia supernovae

We compare observations of the high-redshift galaxy population to the predictions of the galaxy formation model of Croton et al. and De Lucia & Blaizot. This model, implemented on the Millennium Simulation of the concordance Lambda cold dark matter cosmogony, introduces 'radio mode' feedback from the central galaxies of groups and clusters in order to obtain quantitative agreement with the luminosity, colour, morphology and clustering properties of the present-day galaxy population. Here we construct deep light cone surveys in order to compare model predictions to the observed counts and redshift distributions of distant galaxies, as well as to their inferred luminosity and mass functions out to redshift 5. With the exception of the mass functions, all these properties are sensitive to modelling of dust obscuration. A simple but plausible treatment agrees moderately well with most of the data. The predicted abundance of relatively massive  (∼ M _*)  galaxies appears systematically high at high redshift, suggesting that such galaxies assemble earlier in this model than in the real Universe. An independent galaxy formation model implemented on the same simulation matches the observed mass functions slightly better, so the discrepancy probably reflects incomplete or inaccurate galaxy formation physics rather than problems with the underlying cosmogony.     

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.     

On the lack of cold dust in IRAS P09104+4109 and IRAS F15307+3252: their spectral energy distributions and implications for finding dusty AGNs at high redshift

We present upper limits on the 850-μm and 450-μm fluxes of the warm hyperluminous (bolometric luminosity     galaxies IRAS P09104+4109     and IRAS F15307+3252     , derived from measurements using the SCUBA bolometer array on the James Clerk Maxwell Telescope. Hot luminous infrared sources like these are thought to differ from more normal cold ultraluminous infrared     galaxies in that they derive most of their bolometric luminosities from dusty active galactic nuclei (AGNs) as opposed to starbursts. Such hot, dusty AGNs at high redshift are thought to be responsible for much of the mass accretion of the Universe that is in turn responsible for the formation of the supermassive black holes seen in the centres of local galaxies. The galaxy IRAS P09104+4109 is also unusual in that it is a cD galaxy in the centre of a substantial cooling-flow cluster, not an isolated interacting galaxy like most ultraluminous infrared galaxies. Previously it was known to have large amounts of hot     dust from IRAS observations. We now show that the contribution of cold dust to the bolometric luminosity is less than 3 per cent. Most ultraluminous infrared galaxies possess large amounts of cold dust, and it is now known that some cooling-flow cluster cD galaxies do as well. Yet this object, which is an extreme example of both, does not have enough cold gas to contribute significantly to the bolometric luminosity. We outline physical reasons why this could have happened. We then provide a discussion of strategies for finding hot dusty AGNs, given the limitations on submillimetre surveys implied by this work.     

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  .     

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).     

Photochemical evolution of elliptical galaxies – II. The impact of merging-induced starbursts

The effects of late gas accretion episodes and subsequent merger-induced starbursts on the photochemical evolution of elliptical galaxies are studied and compared to the picture of galaxy formation occurring at high redshift with a unique and intense starburst modulated by a very short infall, as suggested by Pipino and Matteucci in Paper I. By means of the comparison with the colour–magnitude relations (CMRs) and the  [〈Mg/Fe〉 _V ]–σ  relation observed in ellipticals, we conclude that either bursts involving a gas mass comparable to the mass already transformed into stars during the first episode of star formation (SF) and occurring at any redshift, or bursts occurring at low redshift (i.e. z ≤ 0.2) and with a large range of accreted mass, are ruled out. These models fail in matching the above relations even if the initial infalling hypothesis is relaxed, and the galaxies form either by means of more complicated SF histories or by means of the classical monolithic model. On the other hand, galaxies accreting a small amount of gas at high redshift (i.e. z ≥ 3) produce a spread in the model results, with respect to the best model of Paper I, which is consistent with the observational scatter of the CMRs, although there is only marginal agreement with the  [〈Mg/Fe〉 _V ]–σ  relation. Therefore, only small perturbations to the standard scenario seem to be allowed. We stress that the strongest constraints to galaxy-formation mechanisms are represented by the chemical abundances, whereas the colours can be reproduced under several different hypotheses.     

A model of supernova feedback in galaxy formation

A model of supernova feedback during disc galaxy formation is developed. The model incorporates infall of cooling gas from a halo, and outflow of hot gas from a multiphase interstellar medium (ISM). The star formation rate is determined by balancing the energy dissipated in collisions between cold gas clouds with that supplied by supernovae in a disc marginally unstable to axisymmetric instabilities. Hot gas is created by thermal evaporation of cold gas clouds in supernova remnants, and criteria are derived to estimate the characteristic temperature and density of the hot component and hence the net mass outflow rate. A number of refinements of the model are investigated, including a simple model of a galactic fountain, the response of the cold component to the pressure of the hot gas, pressure-induced star formation and chemical evolution. The main conclusion of this paper is that low rates of star formation can expel a large fraction of the gas from a dwarf galaxy. For example, a galaxy with circular speed 50 km s¹ can expel 6080 per cent of its gas over a time-scale of 1 Gyr, with a star formation rate that never exceeds 0.1 M yr¹. Effective feedback can therefore take place in a quiescent mode and does not require strong bursts of star formation. Even a large galaxy, such as the Milky Way, might have lost as much as 20 per cent of its mass in a supernova-driven wind. The models developed here suggest that dwarf galaxies at high redshifts will have low average star formation rates and may contain extended gaseous discs of largely unprocessed gas. Such extended gaseous discs might explain the numbers, metallicities and metallicity dispersions of damped Lyman systems.     

The contribution of faint blue galaxies to the submillimetre counts and background

Observations in the submillimetre (submm) waveband have recently revealed a new population of luminous sources. These are proposed to lie at high redshift and to be optically faint because of their high intrinsic dust obscuration. The presence of dust has been previously invoked in optical galaxy count models which use the Bruzual & Charlot evolution models with an exponential τ =9 Gyr star formation rate (SFR) for spirals, and these fit the count data well from U to K . We now show that by using either a 1/ λ or Calzetti absorption law for the dust and re-distributing the evolved spiral galaxy ultraviolet (UV) radiation into the far-infrared (FIR), these models can account for all of the 'faint' ( 1 mJy) 850-μm galaxy counts, but fail to fit 'bright' ( 2 mJy) sources, indicating that another explanation for the submm counts may apply at brighter fluxes, e.g., quasi-stellar objects (QSOs) or ultraluminous infrared galaxies (ULIRGs). We find that the main contribution to the faint, submm number counts is in the redshift range 0.5< z <3, peaking at z ≈1.8. The above model, using either dust law, can also explain a significant proportion of the extragalactic background at 850 μm, as well as producing a reasonable fit to the bright 60-μm IRAS counts.     

Galaxy groups at 0.3 ≤z≤ 0.55 – I. Group properties

The evolution of galaxies in groups may have important implications for the evolution of the star formation history of the Universe, since many processes which operate in groups may suppress star formation and the fraction of galaxies in bound groups grows rapidly between   z = 1  and the present day. In this paper, we present an investigation of the properties of galaxies in galaxy groups at intermediate redshift  ( z ∼ 0.4)  . The groups were selected from the Canadian Network for Observational Cosmology Redshift Survey (CNOC2) redshift survey as described by Carlberg et al., with further spectroscopic follow-up undertaken at the Magellan telescope in order to improve the completeness and depth of the sample. We present the data for the individual groups, and find no clear trend in the fraction of passive galaxies with group velocity dispersion and group concentration. We stack the galaxy groups in order to compare the properties of group galaxies with those of field galaxies at the same redshift. The groups contain a larger fraction of passive galaxies than the field, this trend being particularly clear for galaxies brighter than   M _{B J} < −20  in the higher velocity dispersion groups. In addition, we see evidence for an excess of bright passive galaxies in the groups relative to the field. In contrast, the luminosity functions of the star-forming galaxies in the groups and the field are consistent. These trends are qualitatively consistent with the differences between group and field galaxies seen in the local Universe.     

An Analytic Model of Galactic Winds and Mass Outflows

Galactic winds and mass outflows are observed both in nearby starburst galaxies and in high-redshift star-forming galaxies. We develop a simple analytic model to understand the observed superwind phenomenon with a discussion of the model uncertainties. Our model is built upon the model of McKee & Ostriker for the interstellar medium. It allows one to predict how properties of a superwind, such as wind velocity and mass outflow rate, are related to properties of its starforming host galaxy, such as size, gas density and star formation rate. The model predicts a threshold of star formation rate density for the generation of observable galactic winds. Galaxies with more concentrated star formation activities produce superwinds with higher velocities. The predicted mass outflow rates are comparable to (or slightly larger than) the corresponding star formation rates. We apply our model to both local starburst galaxies and high-redshift Lyman break galaxies, and find its predictions to be in good agreement with current observations. Our model is simple and so can be easily incorporated into numerical simulations and semi-analytical models of galaxy formation.     

Effects of galaxy interactions in different environments

We analyse star formation rates (SFRs) derived from photometric and spectroscopic data of galaxies in pairs in different environments using the 2-degree field galaxy redshift survey (2dFGRS) and the Sloan digital sky survey (SDSS). The two samples comprise several thousand pairs, suitable to explore into detail the dependence of star formation activity in pairs on orbital parameters and global environment. We use the projected galaxy density derived from the fifth brightest neighbour of each galaxy, with a convenient luminosity threshold to characterize environment in both surveys in a consistent way. Star formation activity is derived through the η parameter in 2dFGRS and through the SFR normalized to the total mass in stars,  SFR/ M *  , given by Brinchmann et al. in the SDSS-second data release (SDSS-DR2). For both galaxy pair catalogs, the star formation birth rate parameter is a strong function of the global environment and orbital parameters. Our analysis on SDSS pairs confirms previous results found with the 2dFGRS where suitable thresholds for the star formation activity induced by interactions are estimated at a projected distance   r _p= 100  h ⁻¹ kpc  and a relative velocity  Δ V = 350 km s⁻¹  . We observe that galaxy interactions are more effective at triggering important star formation activity in low- and moderate-density environments with respect to the control sample of galaxies without a close companion. Although close pairs have a larger fraction of actively star-forming galaxies, they also exhibit a greater fraction of red galaxies with respect to those systems without a close companion, an effect that may indicate that dust stirred up during encounters could affect colours and, partially, obscure tidally induced star formation.     

A submillimetre survey of the star formation history of radio galaxies

We present the results of the first major systematic submillimetre survey of radio galaxies spanning the redshift range 1< z <5. The primary aim of this work is to elucidate the star formation history of this sub class of elliptical galaxies by tracing the cosmological evolution of dust mass. Using SCUBA on the JCMT, we have obtained 850-μm photometry of 47 radio galaxies to a consistent rms depth of 1 mJy, and have detected dust emission in 14 cases. The radio galaxy targets have been selected from a series of low-frequency radio surveys of increasing depth (3CRR, 6CE, etc.), in order to allow us to separate the effects of increasing redshift and increasing radio power on submillimetre luminosity. Although the dynamic range of our study is inevitably small, we find clear evidence that the typical submillimetre luminosity (and hence dust mass) of a powerful radio galaxy is a strongly increasing function of redshift; the detection rate rises from ≃15 per cent at z <2.5 to ≳75 per cent at z >2.5, and the average submillimetre luminosity rises at a rate ∝(1+ z )³ out to z ≃4. Moreover, our extensive sample allows us to argue that this behaviour is not driven by underlying correlations with other radio galaxy properties such as radio power, radio spectral index, or radio source size/age. Although radio selection may introduce other more subtle biases, the redshift distribution of our detected objects is in fact consistent with the most recent estimates of the redshift distribution of comparably bright submillimetre sources discovered in blank field surveys. The evolution of submillimetre luminosity found here for radio galaxies may thus be representative of massive ellipticals in general.     

A census of metals at high and low redshift and the connection between submillimetre sources and spheroid formation

Deep surveys in many wavebands have shown that the rate at which stars were forming was at least a factor of 10 higher at redshifts >1 than today. Heavy elements ('metals') are produced by stars, and the star formation history deduced by these surveys implies that a significant fraction of all metals in the Universe today should already exist at   z ∼ 2–3  . However, only 10 per cent of the total metals expected to exist at this redshift have so far been accounted for (in damped Lyman α absorbers and the Lyman forest). In this paper, we use the results of submillimetre surveys of the local and high-redshift Universe to show that there was much more dust in galaxies in the past. We find that a large proportion of the missing metals are traced by this dust, bringing the metals implied from the star formation history and observations into agreement. We also show that the observed distribution of dust masses at high redshift can be reproduced remarkably well by a simple model for the evolution of dust in spheroids, suggesting that the descendants of the dusty galaxies found in deep submillimetre surveys are the relatively dust-free spiral bulges and ellipticals in the Universe today.     

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.     

Evolutionary synthesis models for galaxy transformation in clusters

The galaxy population in rich local galaxy clusters shows a ratio of one quarter elliptical galaxies, two quarters S0 galaxies, and one quarter spiral galaxies. Observations of clusters at redshift 0.5 show a perspicuously different ratio, the dominant galaxy type are spiral galaxies with a fraction of two quarters while the number of S0 galaxies decreases to a fraction of one quarter (Dressler et al. 1997). This shows an evolution of the galaxy population in clusters since redshift 0.5 and it has been suspected that galaxy transformation processes during the infall into a cluster are responsible for this change. These could be merging, starburst or ram-pressure stripping. We use our evolutionary synthesis models to describe various possible effects of those interactions on the star formation of spiral galaxies infalling into clusters. We study the effects of starbursts of various strengths as well as of the truncation of star formation at various epochs on the color and luminosity evolution of model galaxies of various spectral types. As a first application we present the comparison of our models with observed properties of the local S0 galaxy population to constrain possible S0 formation mechanisms in clusters. Application to other types of galaxies is planned for the future. This revised version was published online in August 2006 with corrections to the Cover Date.     

Confirming a population of hot-dust dominated, star-forming, ultraluminous galaxies at high redshift

We identify eight   z > 1  radio sources undetected at 850 μm but robustly detected at 70 μm, confirming that they represent ultraluminous infrared galaxies (ULIRGs) with hotter dust temperatures  (〈 T _d〉= 52 ± 10 K)  than submillimetre galaxies (SMGs) at similar luminosities and redshifts. These galaxies share many properties with SMGs: ultraviolet spectra consistent with starbursts, high stellar masses and radio luminosities. We can attribute their radio emission to star formation since high-resolution Multi-Element Radio Linked Interferometer Network (MERLIN) radio maps show extended emission regions (with characteristic radii of 2–3 kpc), which are unlikely to be generated by active galactic nucleus (AGN) activity. These observations provide the first direct confirmation of hot, dusty ULIRGs which are missed by current submillimetre surveys. They have significant implications for future observations from the Herschel Space Observatory and Submillimetre Common-User Bolometer Array 2 (SCUBA2), which will select high-redshift luminous galaxies with less selection biases.     

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  .