Hyperluminous infrared galaxies

39 galaxies are now known, from follow-up of faint IRAS sources and from submillimetre observations of high-redshift AGN, with far-infrared luminosities >10¹³ L_⊙. 13 of these, which have been found in 60- or 850-μm surveys, form an important unbiased subsample. 12 have been found by comparison of 60-μm surveys with quasar or radio galaxy catalogues, or from infrared surveys with colour selection biased towards AGN, while a further 14 have been found through submillimetre observations of known high-redshift AGN. In this paper I argue, on the basis of detailed modelling of the spectral energy distributions of hyperluminous galaxies with accurate radiative transfer models, and from evidence of high gas mass in several cases, that the bulk of the emission from these galaxies at rest frame wavelengths ≥50 μm is caused by star formation. Even after correction for the effects of lensing, hyperluminous galaxies with emission peaking at rest frame wavelengths ≥50 μm are therefore undergoing star formation at rates >10³ M_⊙ yr⁻¹ and are strong candidates for being primeval galaxies, in the process of a major episode of star formation.     

IRAC photometric analysis and the mid-IR photometric properties of Lyman-break galaxies

We present photometric analysis of deep mid-infrared (mid-IR) observations obtained by Spitzer /IRAC covering the fields Q1422+2309, Q2233+1341, DSF2237a,b, HDFN, SSA22a,b and B20902+34, giving the number counts and the depths for each field. In a sample of 751 Lyman-break galaxies (LBGs) lying in those fields, 443, 448, 137 and 152 are identified at 3.6-, 4.5-, 5.8-, 8.0-μm IRAC bands, respectively, expanding their spectral energy distribution to rest-near-IR and revealing that LBGs display a variety of colours. Their rest-near-IR properties are rather inhomogeneous, ranging from those that are bright in IRAC bands and exhibit  [ R ]−[3.6] > 1.5  colours to those that are faint or not detected at all in IRAC bands with  [ R ]−[3.6] < 1.5  colours and these two groups of LBGs are investigated. We compare the mid-IR colours of the LBGs with the colours of star-forming galaxies and we find that LBGs have colours consistent with star-forming galaxies at   z ∼ 3  . The properties of the LBGs detected in the 8-μm IRAC band (rest-frame K band) are examined separately, showing that they exhibit redder  [ R ]−[3.6]  colours than the rest of the population and that although in general, a multiwavelength study is needed to reach more secure results, IRAC 8-μm band can be used as a diagnostic tool, to separate high z , luminous AGN-dominated objects from normal star-forming galaxies at   z ∼ 3  .     

Spitzer and ISO galaxy counts in the mid-infrared

Galaxy source counts that simultaneously fit the deep mid-infrared surveys at 24 microns and 15 microns made by the Spitzer Space Telescope and the Infrared Space Observatory ( ISO ), respectively, are presented for two phenomenological models. The models are based on starburst and luminous infrared galaxy dominated populations. Both models produce excellent fits to the counts in both wavebands and provide an explanation for the high-redshift population seen in the longer Spitzer 24-micron band supporting the hypothesis that they are luminous–ultraluminous infrared galaxies at   z = 2–3  , being the mid-infrared counterparts to the submillimetre galaxy population. The source counts are characterized by strong evolution to redshift unity, followed by less drastic evolution to higher redshift. The number–redshift distributions in both wavebands are well explained by the effect of the many mid-infrared features passing through the observation windows. The sharp upturn at around a millijansky in the 15-μm counts in particular depends critically on the distribution of mid-infrared features around 12 μm, in the assumed spectral energy distribution.     

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.     

The K-band Hubble diagram of submillimetre galaxies and hyperluminous galaxies

We present the K -band Hubble diagrams ( K – z relations) of submillimetre-selected galaxies and hyperluminous galaxies (HLIRGs). We report the discovery of a remarkably tight K – z relation of HLIRGs, indistinguishable from that of the most luminous radio galaxies. Like radio galaxies, the HLIRG K – z relation at   z ≲ 3  is consistent with a passively evolving ∼3 L _* instantaneous starburst starting from a redshift of   z ∼ 10  . In contrast, many submillimetre-selected galaxies are ≳2 mag fainter, and the population has a much larger dispersion. We argue that dust obscuration and/or a larger mass range may be responsible for this scatter. The galaxies so far proved to be hyperluminous may have been biased towards higher AGN bolometric contributions than submillimetre-selected galaxies due to the 60-μm selection of some, so the location on the K – z relation may be related to the presence of the most massive active galactic nucleus. Alternatively, a particular host galaxy mass range may be responsible for both extreme star formation and the most massive active nuclei.     

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.     

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.     

Starburst galaxies and structure in the submillimetre background towards the Hubble Deep Field

We use an 850-μm SCUBA map of the Hubble Deep Field (HDF) to study the dust properties of optically-selected starburst galaxies at high redshift. The optical/infrared (IR) data in the HDF allow a photometric redshift to be estimated for each galaxy, together with an estimate of the visible star-formation rate. The 850-μm flux density of each source provides the complementary information: the amount of hidden, dust-enshrouded star formation activity. Although the 850-μm map does not allow detection of the majority of individual sources, we show that the galaxies with the highest UV star-formation rates are detected statistically, with a flux density of about S ₈₅₀=0.2 mJy for an apparent UV star-formation rate of 1  h ⁻² M_⊙ yr⁻¹. This level of submillimetre output indicates that the total star-forming activity is on average a factor of approximately 6 times larger than the rate inferred from the UV output of these galaxies. The general population of optical starbursts is then predicted to contribute at least 25 per cent of the 850-μm background. We carry out a power-spectrum analysis of the map, which yields some evidence for angular clustering of the background source population, but at a level lower than that seen in Lyman-break galaxies. Together with other lines of argument, particularly from the NICMOS HDF data, this suggests that the 850-μm background originates over an extremely wide range of redshifts – perhaps 1≲ z ≲6.     

A new model for infrared and submillimetre counts

A new model for source counts from 8 to 1100 μm is presented, which agrees well with source-count data and the observed background spectrum. The model assumes different evolution for each of the four infrared template types used. The evolution is modified in two ways compared to my 2001 model: (i) the exponential factor is modified so that it tends to a constant value at late times and (ii) the power-law factor is modified so that it tends to zero at redshift   z _f   , rather than zero as assumed previously. I find strong evidence from the 850 and 1100 μm counts, and from the infrared background, that   z _f = 4–5  , with some preference for a value at the low end of the range, implying that star-forming galaxies at   z > 5  are not significant infrared emitters, presumably due to a low opacity in dust at these early epochs.
The model involves zero or even negative evolution for starbursts and active galactic nuclei at low redshifts (<0.2), suggesting that the era of major mergers and strong galaxy–galaxy interactions is over.