A search for the submillimetre counterparts to Lyman break galaxies

We have carried out targeted submillimetre observations as part of a programme to explore the connection between the rest-frame ultraviolet and far-infrared properties of star-forming galaxies at high redshift, which is currently poorly understood. On the one hand, the Lyman break technique is very effective at selecting     galaxies. On the other, 'blank-field' imaging in the submillimetre seems to turn up sources routinely, amongst which some are star-forming galaxies at similar redshifts. Already much work has been done searching for optical identifications of objects detected using the SCUBA instrument. Here we have taken the opposite approach, performing submillimetre photometry for a sample of Lyman break galaxies, the ultraviolet properties of which imply high star formation rates. The total signal from our Lyman break sample is undetected in the submillimetre, at an rms level of ∼0.5 mJy, which implies that the population of Lyman break galaxies does not constitute a large part of the recently detected blank-field submillimetre sources. However, our one detection suggests that with reasonable SCUBA integrations we might expect to detect those few Lyman break galaxies that are far-infrared brightest.     

The local far-infrared galaxy colour–luminosity distribution: a reference for BLAST and Herschel/SPIRE submillimetre surveys

We measure the local galaxy far-infrared (FIR) 60 to 100 μm colour–luminosity distribution using an all-sky IRAS survey. This distribution is an important reference for the next generation of FIR–submillimetre surveys that have and will conduct deep extragalactic surveys at 250–500 μm. With the peak in dust-obscured star-forming activity leading to present-day giant ellipticals now believed to occur in submillimetre galaxies near   z ∼ 2.5  , these new FIR–submillimetre surveys will directly sample the spectral energy distributions of these distant objects at rest-frame FIR wavelengths similar to those at which local galaxies were observed by IRAS . We have taken care to correct for the temperature bias and the evolution effects in our IRAS 60-μm-selected sample. We verify that our colour–luminosity distribution is consistent with the measurements of the local FIR luminosity function, before applying it to the higher redshift Universe. We compare our colour–luminosity correlation with recent dust–temperature measurements of submillimetre galaxies and find evidence for pure luminosity evolution of the form  (1 + z )³  . This distribution will be useful for the development of evolutionary models for Balloon-borne Large Aperture Submillimeter Telescope (BLAST) and Spectral and Photometric Imaging Receiver (SPIRE) surveys as it provides a statistical distribution of the rest-frame dust temperatures for galaxies as a function of luminosity.     

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.     

The correlation between black hole mass and bulge velocity dispersion in hierarchical galaxy formation models

Popular models for the origin of gamma-ray bursts (GRBs) include short-lived massive stars as the progenitors of the fireballs. Hence the redshift distribution of GRBs should track the cosmic star formation rate of massive stars accurately. A significant proportion of high-mass star formation activity appears to occur in regions that are obscured from view in the optical waveband by interstellar dust. The amount of dust-enshrouded star formation activity taking place has been estimated by observing the thermal radiation from the dust that has been heated by young stars in the far-infrared and submillimetre wavebands. Here we discuss an alternative probe – the redshift distribution of GRBs. GRBs are detectable at the highest redshifts, and because gamma-rays are not absorbed by dust, the redshift distribution of GRBs should therefore be unaffected by dust extinction. At present the redshifts of GRBs can only be determined from the associated optical transient emission; however, useful information about the prevalence of dust-obscured star formation can also be obtained from the ratio of GRBs with and without an associated optical transient. Eight GRBs currently have spectroscopic redshifts. Once about a hundred redshifts are known, the population of GRBs will provide an important test of different models of the star formation history of the Universe.     

The nature of high-redshift galaxies

Using semi-analytic models of galaxy formation set within the cold dark matter (CDM) merging hierarchy, we investigate several scenarios for the nature of the high-redshift     ) Lyman-break galaxies (LBGs). We consider a 'collisional starburst' model in which bursts of star formation are triggered by galaxy–galaxy mergers, and find that a significant fraction of LBGs are predicted to be starbursts. This model reproduces the observed comoving number density of bright LBGs as a function of redshift and the observed luminosity function at     and     with a reasonable amount of dust extinction. Model galaxies at     have star formation rates, half-light radii,     colours and internal velocity dispersions that are in good agreement with the data. Global quantities such as the star formation rate density and cold gas and metal content of the Universe as a function of redshift also agree well. Two 'quiescent' models without starbursts are also investigated. In one, the star formation efficiency in galaxies remains constant with redshift, while in the other, it scales inversely with disc dynamical time, and thus increases rapidly with redshift. The first quiescent model is strongly ruled out, as it does not produce enough high-redshift galaxies once realistic dust extinction is accounted for. The second quiescent model fits marginally, but underproduces cold gas and very bright galaxies at high redshift. A general conclusion is that star formation at high redshift must be more efficient than locally. The collisional starburst model appears to accomplish this naturally without violating other observational constraints.     

The relationship between star formation rate and radio synchrotron luminosity at 0 < z < 2

We probe the relationship between star formation rate (SFR) and radio synchrotron luminosity in galaxies at  0 < z < 2  within the northern Spitzer Wide-area Infrared Extragalactic survey (SWIRE) fields, in order to investigate some of the assumptions that go into calculating the star formation history of the Universe from deep radio observations. We present new 610-MHz Giant Metrewave Radio Telescope (GMRT) observations of the European Large-Area ISO Survey-North 2 (ELAIS-N2) field, and using this data, along with previous GMRT surveys carried out in the ELAIS-N1 (North 1) and Lockman Hole regions, we construct a sample of galaxies which have redshift and SFR information available from the SWIRE survey. We test whether the local relationship between SFR and radio luminosity is applicable to   z = 2  galaxies, and look for evolution in this relationship with both redshift and SFR in order to examine whether the physical processes which lead to synchrotron radiation have remained the same since the peak of star formation in the Universe. We find that the local calibration between radio luminosity and star formation can be successfully applied to radio-selected high-redshift, high-SFR galaxies, although we identify a small number of sources where this may not be the case; these sources show evidence for inaccurate estimations of their SFR, but there may also be some contribution from physical effects such as the recent onset of starburst activity, or suppression of the radio luminosity within these galaxies.     

Millimetric properties of gamma-ray burst host galaxies

We present millimetre (mm) and submillimetre (submm) photometry of a sample of five host galaxies of gamma-ray bursts (GRBs), obtained using the Max Planck Millimetre Bolometer (MAMBO2) array and Submillimetre Common-User Bolometer Array (SCUBA). These observations were obtained as part of an ongoing project to investigate the status of GRBs as indicators of star formation. Our targets include two of the most unusual GRB host galaxies, selected as likely candidate submm galaxies: the extremely red  ( R − K ≈ 5)  host of GRB 030115, and the extremely faint  ( R > 29.5)  host of GRB 020124. Neither of these galaxies is detected, but the deep upper limits for GRB 030115 impose constraints on its spectral energy distribution, requiring a warmer dust temperature than is commonly adopted for submillimetre galaxies (SMGs).
As a framework for interpreting these data, and for predicting the results of forthcoming submm surveys of Swift -derived host samples, we model the expected flux and redshift distributions based on luminosity functions of both submm galaxies and GRBs, assuming a direct proportionality between the GRB rate density and the global star formation rate density. We derive the effects of possible sources of uncertainty in these assumptions, including (1) introducing an anticorrelation between GRB rate and the global average metallicity, and (2) varying the dust temperature.     

Galaxy formation and evolution: low-surface-brightness galaxies

We investigate in detail the hypothesis that low-surface-brightness galaxies (LSBs) differ from ordinary galaxies simply because they form in haloes with large spin parameters. We compute star formation rates using the Schmidt law, assuming the same gas infall dependence on surface density as used in models of the Milky Way. We build stellar population models, predicting colours, spectra and chemical abundances. We compare our predictions with observed values of metallicity and colours for LSBs, and find excellent agreement with all observables. In particular, integrated colours, colour gradients, surface brightness and metallicity match very well to the observed values of LSBs for models with ages larger than 7 Gyr and high values (λ > 0.05) for the spin parameter of the haloes. We also compute the global star formation rate (SFR) in the Universe due to LSBs, and show that it has a flatter evolution with redshift than the corresponding SFR for normal discs. We furthermore compare the evolution in redshift of [ Zn / H ] for our models to those observed in damped Lyman α systems by Pettini et al. and show that damped Lyman α system abundances are consistent with the predicted abundances at different radii for LSBs. Finally, we show how the required late redshift of collapse of the halo may constrain the power spectrum of fluctuations.     

The far-infrared–radio relationship at high and low redshift

We use the results of the SCUBA Local Universe Galaxy Survey, a submillimetre (submm) survey of galaxies in the nearby Universe, to investigate the relationship between the far-infrared (FIR)–submm and radio emission of galaxies at both low and high redshift. At low redshift we show that the correlation between radio and FIR emission is much stronger than the correlation between radio and submm emission, which is evidence that massive stars are the source of both the FIR and radio emission. At high redshift we show that the submm sources detected by SCUBA are brighter sources of radio emission than are predicted from the properties of galaxies in the local Universe. We discuss possible reasons for the cosmic evolution of the relationship between radio and FIR emission.     

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.     

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.     

The link between SCUBA and Spitzer: cold galaxies at z≲ 1

We show that the far-IR properties of distant Luminous and UltraLuminous InfraRed Galaxies (LIRGs and ULIRGs, respectively) are on average divergent from analogous sources in the local Universe. Our analysis is based on Spitzer Multiband Imaging Photometer (MIPS) and Infrared Array Camera (IRAC) data of   L _IR > 10¹⁰ L_⊙, 70 μm  selected objects in the  0.1 < z < 2  redshift range and supported by a comparison with the IRAS Bright Galaxy Sample. The majority of the objects in our sample are described by spectral energy distributions (SEDs) which peak at longer wavelengths than local sources of equivalent total infrared luminosity. This shift in SED peak wavelength implies a noticeable change in the dust and/or star-forming properties from   z ∼ 0  to the early Universe, tending towards lower dust temperatures, indicative of strong evolution in the cold dust, 'cirrus', component. We show that these objects are potentially the missing link between the well-studied local IR-luminous galaxies, Spitzer IR populations and SCUBA sources – the   z < 1  counterparts of the cold   z > 1  SubMillimetre Galaxies (SMGs) discovered in blank-field submillimetre surveys. The Herschel Space Observatory is well placed to fully characterize the nature of these objects, as its coverage extends over a major part of the far-IR/sub-mm SED for a wide redshift range.     

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.     

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 ultraviolet-selected galaxy redshift survey: new estimates of the local star formation rate

We present the first results of an ongoing spectroscopic survey of galaxies selected in the rest frame ultraviolet (UV). The source catalogue has been constructed from a flux-limited sample of stars, galaxies and QSOs imaged at 2000 Å in Selected Area 57 with the FOCA balloon-borne imaging camera. Accurate positions for the UV sources have been obtained by matching with optical counterparts using APM scans of the Palomar Sky Survey limited at B   20.5. Here we present results derived from optical spectroscopy conducted with the WIYN telescope and the WHT for 142 faint sources. The redshift distribution for this UV-selected sample extends over 0 <  z  < 0.5, and a high fraction of the sources show intense nebular emission lines and UV–optical colours bluer than normal Hubble sequence galaxies. Such UV-selected surveys are thus a very efficient way to locate and study intermediate-redshift galaxies undergoing intense star formation. Although our sample is currently small, we derive a rest frame UV luminosity function with a steep faint-end slope consistent with that found for late-type galaxies in optical samples. However, the integrated luminosity density derived implies a volume-averaged star formation rate higher than other recent estimates, assuming a normal initial mass function. If representative of other UV fields, as suggested by UV number count studies, our data imply that the local abundance of star-forming galaxies may have been underestimated, and consequently claims for strong evolution in the global star formation rate in the range 0 <  z  < 1 overstated. An intensive study of a large UV-selected sample is likely to reveal important information on the declining population of star-forming galaxies of all types.     

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.     

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.     

Star formation history up to z= 7.4: implications for gamma-ray bursts and cosmic metallicity evolution

The current Swift sample of gamma-ray bursts (GRBs) with measured redshifts allows us to test the assumption that GRBs trace star formation in the Universe. Some authors have claimed that the rate of GRBs increases with cosmic redshift faster than the star formation rate, whose cause is not yet known. In this paper, I investigate the possibility of interpreting the observed discrepancy between the GRB rate history and the star formation rate history using cosmic metallicity evolution. I am motivated by the observation that cosmic metallicity evolves with redshift and GRBs tend to occur in low-metallicity galaxies. First, I derive a star formation history up to redshift   z = 7.4  from an updated sample of star formation rate densities. This is obtained by adding the new ultraviolet measurements of Bouwens et al. and the new ultraviolet and infrared measurements of Reddy et al. to the existing sample compiled by Hopkins & Beacom. Then, adopting a simple model for the relation between GRB production and the cosmic metallicity history as proposed by Langer & Norman, I show that the observed redshift distribution of the Swift GRBs can be reproduced with good accuracy. Although the results are limited by the small size of the GRB sample and the poorly understood selection biases in detection and localization of GRBs and in redshift determination, they suggest that GRBs trace both star formation and metallicity evolution. If the star formation history can be accurately measured with other approaches, which is presumably achievable in the near future, it will be possible to determine the cosmic metallicity evolution using the study of the redshift distribution of GRBs.     

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.     

The effects of photoionization on galaxy formation – II. Satellite galaxies in the Local Group

We use a self-consistent model of galaxy formation and the evolution of the intergalactic medium to study the effects of the reionization of the Universe at high redshift on the properties of satellite galaxies like those seen around the Milky Way. Photoionization suppresses the formation of small galaxies, so that surviving satellites are preferentially those that formed before the Universe reionized. As a result, the number of satellites expected today is about an order of magnitude smaller than the number inferred by identifying satellites with subhaloes of the same circular velocity in high-resolution simulations of the dark matter. The resulting satellite population has an abundance similar to that observed in the Local Group, although the distribution of circular velocities differs somewhat from the available data. We explore many other properties of satellite galaxies, including their gas content, metallicity and star formation rate, and find generally good agreement with available data. Our model predicts the existence of many as yet undetected satellites in the Local Group. We quantify their observability in terms of their apparent magnitude and surface brightness, and also in terms of their constituent stars. A near-complete census of the Milky Way's satellites would require imaging to   V ≈20  and to a surface brightness fainter than 26 V -band magnitudes per square arcsecond. Satellites with integrated luminosity   V =15  should contain of order 100 stars brighter than   B =26  , with central stellar densities of a few tens per square arcminute. Discovery of a large population of faint satellites would provide a strong test of current models of galaxy formation.