Gravitational lensing and the Sunyaev–Zel'dovich effect in the millimetre/submillimetre waveband

The intensity of the cosmic microwave background radiation in the fields of clusters of galaxies is modified by inverse Compton scattering in the hot intracluster gas — the Sunyaev–Zel'dovich (SZ) effect. The effect is expected to be most pronounced at a frequency of about 350 GHz (a wavelength of about 800 μm), and has been detected in the centimetre and millimetre wavebands. In the millimetre/submillimetre waveband, the gravitationally lensed images of distant dusty star-forming galaxies in the background of the cluster are predicted to dominate the appearance of clusters on scales of several arcsec, and could confuse observations of the SZ effect at frequencies greater than about 200 GHz (wavelengths shorter than about 1.5 mm). Recent observations by Smail, Ivison &38; Blain confirm that a significant population of confusing sources is present in this waveband. Previous estimates of source confusion in observations of the millimetre/submillimetre-wave SZ effect did not include the effects of lensing by the cluster, and so the accuracy of such measurements could be lower than expected. Source subtraction may be required in order to measure the SZ effect accurately, and a careful analysis of the results of an ensemble of SZ measurements could be used to impose limits to the form of evolution of distant dusty star-forming galaxies.     

The Planck Surveyor mission and gravitational lenses

An extremely sensitive all-sky survey will be carried out in the millimetre/submillimetre waveband by the forthcoming ESA mission Planck Surveyor . The main scientific goal of the mission is to make very accurate measurements of the spatial power spectrum of primordial anisotropies in the cosmic microwave background radiation; however, hundreds of thousands of distant dusty galaxies and quasars will also be detected. These sources are much more likely to be gravitationally lensed by intervening galaxies compared with sources discovered in surveys in other wavebands. Here the number of lenses expected in the survey is estimated, and techniques for discriminating between lensed and unlensed sources are discussed. A practical strategy for this discrimination is presented, based on exploiting the remarkable sensitivity and resolving power of large ground-based millimetre/submillimetre-wave interferometer arrays. More than a thousand gravitational lenses could be detected: a sample that would be an extremely valuable resource in observational cosmology.     

Submillimetre-wave gravitational lenses and cosmology

One of the most direct routes for investigating the geometry of the Universe is provided by the numbers of strongly magnified gravitationally lensed galaxies as compared with those that are either weakly magnified or de-magnified. In the submillimetre waveband the relative abundance of strongly lensed galaxies is expected to be larger as compared with the optical or radio wavebands, both in the field and in clusters of galaxies. The predicted numbers depend on the properties of the population of faint galaxies in the submillimetre waveband, which was formerly very uncertain; however, recent observations of lensing clusters have reduced this uncertainty significantly and confirm that a large sample of galaxy–galaxy lenses could be detected and investigated using forthcoming facilities, including the FIRST and Planck Surveyor space missions and a large ground-based millimetre/submillimetre-wave interferometer array (MIA). We discuss how this sample could be used to impose limits on the values of cosmological parameters and the total density and form of evolution of the mass distribution of bound structures, even in the absence of detailed lens modelling for individual members of the sample. The effects of different world models on the form of the magnification bias expected in sensitive submillimetre-wave observations of clusters are also discussed, because an MIA could resolve and investigate images in clusters in detail.     

Submillimetre detection of the z=2.83 Lyman-break galaxy, Westphal–MM8, and implications for SCUBA2

We present confusion-limited submillimetre (submm) observations with the Submillimetre Common-User Bolometer Array (SCUBA) camera on the James Clerk Maxwell Telescope of the   z = 2.83  Lyman-break galaxy (LBG), Westphal–MM8, reaching an 850 μm sensitivity even greater than that achieved in the SCUBA map of the Hubble Deep Field region. The detection of MM8  ( S _850 μm= 1.98 ± 0.48 mJy)  , along with the literature submm detections of lensed LBGs, suggests that the LBG population may contribute significantly to the source counts of submm-selected galaxies in the 1–2 mJy regime. Additionally, submm-luminous LBGs are a viable progenitor population for the recently discovered evolved galaxies at   z ∼ 2–2.5  . These observations represent an important baseline for SCUBA2 observations which will regularly map large regions of the sky to this depth.     

Millimetre/submillimetre-wave emission-line searches for high-redshift galaxies

The redshifted spectral line radiation emitted from both atomic fine-structure and molecular rotational transitions in the interstellar medium (ISM) of high-redshift galaxies can be detected in the centimetre, millimetre and submillimetre wavebands. Here we predict the counts of galaxies detectable in an array of molecular and atomic lines. This calculation requires a reasonable knowledge of both the surface density of these galaxies on the sky, and the physical conditions in their ISM. The surface density is constrained using the results of submillimetre-wave continuum surveys. Follow-up OVRO Millimeter Array observations of two of the galaxies detected in the dust continuum have provided direct measurements of CO rotational line emission at redshifts of 2.56 and 2.81. Based on these direct high-redshift observations and on models of the ISM that are constrained by observations of low-redshift ultraluminous infrared galaxies, we predict the surface density of line-emitting galaxies as a function of line flux density and observing frequency. We incorporate the sensitivities and mapping speeds of existing and future millimetre/submillimetre-wave telescopes and spectrographs, and so assess the prospects for blank-field surveys to detect this line emission from gas-rich high-redshift galaxies.     

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.     

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.     

Millimetre and submillimetre τ relations

Relations are presented between the optical depth of the atmosphere  (τ)  in the millimetre and submillimetre windows at 350, 450, 850 and 1250 μm. These relations are derived by an analysis of 849 skydips at 850 μm performed with the submillimetre camera SCUBA at the James Clerk Maxwell Telescope (JCMT) and polynomial fits to over 19 000 measurements of a tipping radiometer operating at 350 μm. The relations allow for an improved calibration of millimetre and submillimetre observations on Mauna Kea, Hawai'i, and for better monitoring of atmospheric conditions.     

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 effect of lensing on the identification of SCUBA galaxies

Spectroscopic surveys of submillimetre(submm)-selected sources have uncovered optically bright galaxies at   z ≲ 1  close to the positions of several submm sources. Naive statistical analyses suggest that these galaxies are associated with the submm emission. However, in some cases, it is difficult to understand this association given the relatively modest redshifts and unprepossessing spectral characteristics of the galaxies. These are in stark contrast to those expected from the massive dust-enshrouded starbursts and AGN thought to power the bulk of the bright submm population. We present new observations of optically bright counterparts to two luminous submm sources, along with a compilation of previously proposed optically bright counterparts with   z ≲ 1  . We suggest that the majority of these associations between bright galaxies and submm sources may be as a result of the action of the foreground galaxies as gravitational lenses on the much fainter and more distant submm sources. We discuss the implications of this conclusion for our understanding of the SCUBA population.     

Evolutionary constraints from infrared source counts

The backward evolution approach to modelling galaxy source counts is re-visited in the wake of the numerous results and revelations from the Infrared Space Observatory ( ISO ), the Submillimetre Common User Bolometer Array (SCUBA) and the detections and measurements of the cosmic extragalactic background light. Using the framework of the Pearson & Rowan-Robinson galaxy evolution model, the observed source counts and background measurements are used to constrain the evolution in the galaxy population. It is found that a strong evolution in both density and luminosity of the high-luminosity tail of the infrared (IR) luminosity function, interpreted as the ultraluminous galaxies discovered first by IRAS and later elevated in status by SCUBA and ISO , can account for the source counts from 15 μm (where it matches the undulations in the integral counts and the hump in the differential counts extremely well) to the submillimetre region, as well as explain the peak in the cosmic infrared background at ∼140 μm. The submillimetre counts are interpreted as the superposition of two separate populations comprising ultraluminous infrared galaxies (ULIGs) at the brighter submillimetre fluxes and starburst galaxies at fluxes fainter than ∼2 mJy. In this scenario the high-redshift ULIGs are tenuously interpreted as the progenitors of today's giant elliptical (gE) galaxies.
All the source count models can be accessed via the world wide web at the URL http://www.ir.isas.ac.jp/~cpp/counts/     

The SCUBA 8-mJy survey – II. Multiwavelength analysis of bright submillimetre sources

We present the results of a multiwavelength study of the 19 most significant submillimetre (submm) sources detected in the SCUBA 8-mJy survey. As described in Scott et al. , this survey covers ≃260 arcmin² using the submillimetre camera SCUBA, to a limiting source detection limit   S _850 μm≃8 mJy  . One advantage of this relatively bright flux-density limit is that accurate astrometric positions are potentially achievable for every source using existing radio and/or millimetre-wave interferometers. However, an associated advantage is that spectral energy distribution (SED) based redshift constraints should be more powerful than in fainter submm surveys. Here we therefore exploit the parallel SCUBA 450-μm data, in combination with existing radio and Infrared Space Observatory ( ISO ) data at longer and shorter wavelengths to set constraints on the redshift of each source. We also analyse new and existing optical and near-infrared imaging of our SCUBA survey fields to select potential identifications consistent with these constraints. Our derived SED-based redshift constraints, and the lack of statistically significant associations with even moderately bright galaxies allow us to conclude that all 19 sources lie at   z >1  , and at least half of them apparently lie at   z >2  .     

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.     

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.     

The temperatures of dust-enshrouded active galactic nuclei

A high density of massive dark objects (MDOs), probably supermassive black holes, in the centres of nearby galaxies has been inferred from recent observations. There are various indications that much of the accretion responsible for producing these objects took place in dust-enshrouded active galactic nuclei (AGNs). If so, then measurements of the intensity of background radiation and the source counts in the far-infrared and submillimetre wavebands constrain the temperature of dust in these AGNs. An additional constraint comes from the hard X-ray background, if this is produced by accretion. One possibility is that the dust shrouds surrounding the accreting AGNs are cold, about 30 K. In this event, the dusty AGNs could be some subset of the population of luminous distant sources discovered at 850 μm using the SCUBA array on the James Clerk Maxwell Telescope, as proposed by Almaini, Lawrence & Boyle. An alternative is that the dust shrouds surrounding the accreting AGNs are much hotter (>60 K). These values are closer to the dust temperatures of a number of well-studied low-redshift ultraluminous galaxies that are thought to derive their power from accretion. If the local MDO density is close to the maximum permitted, then cold sources cannot produce this density without the submillimetre background being overproduced if they accrete at high radiative efficiency, and thus a hot population is required. If the dust-enshrouded accretion occurred at similar redshifts to that taking place in unobscured optical quasars, then a significant fraction of the far-infrared background radiation measured by COBE at 140 μm, but very little of the submillimetre background at 850 μm, may have been produced by hot dust-enshrouded AGNs which may have already been seen in recent Chandra X-ray surveys.     

The diversity of SCUBA-selected galaxies

We present extensive observations of a sample of distant, submillimetre (submm) galaxies detected in the field of the massive cluster lens Abell 1835, using the Submm Common-User Bolometer Array (SCUBA). Taken in conjunction with earlier observations of other submm-selected sources, we now have detailed, multiwavelength observations of seven examples of the submm population, having exploited the combination of achromatic amplification by cluster lenses and lavish archival data sets. These sources, all clearly at z ≳1, illustrate the wide range in the radio and optical properties of distant submm-selected galaxies. We include detailed observations of the first candidate 'pure' starburst submm galaxy at high redshift, a z =2.56 interacting galaxy which shows no obvious sign of hosting an AGN. The remaining sources have varying degrees of inferred AGN activity (three out of seven of the most luminous show some evidence of the presence of an AGN), although even when an AGN is obviously present it is still not apparent whether reprocessed radiation from this source dominates the submm emission. In contrast with the variation in the spectral properties, we see relatively homogeneous morphologies for the population, with a large fraction of merging or interacting systems. Our study shows that virtually identical spectral energy distributions are seen for galaxies that exhibit strikingly different optical/UV spectral-line characteristics. We conclude that standard optical/UV spectral classifications are misleading when applied to distant, highly obscured galaxies, and that we must seek other means of determining the various contributions to the overall energy budget of submm galaxies and hence to the far-infrared extragalactic background.     

Extended dust emission in NGC 7465

We present SCUBA 850-μm, JCMT  CO( J =2→1)  , B -band imaging and VLA H  i observations of the NGC 7465/4/3 group of galaxies. The 850-μm emission associated with NGC 7465 extends to at least ∼2 R ₂₅ and is well correlated with the H  i . We investigate a range of possible mechanisms by which dust beyond R ₂₅ may be heated to give the observed extended submillimetre emission. By modelling the dust heating by stars in two extreme geometries, we fail to find any reasonable star formation scenario that is consistent with both the 850-μm and optical data. Furthermore, we do not detect any  CO( J =2→1)  emission coincident with the extended dust and atomic gas as would be expected if significant star formation were occurring. We show that shock-heating of dust via cloud–cloud collisions in the stripped interstellar medium of NGC 7465 could be sufficient to explain the extended 850-μm emission and lack of optical emission in the stripped gas, and suggest that cloud–cloud collisions may be an important dust heating mechanism in gas-rich systems.     

Exploiting magnification bias in ultradeep submillimetre-wave surveys using ALMA

The surface density of populations of galaxies with steep/shallow source counts is increased/decreased by gravitational lensing magnification. These effects are usually called 'magnification bias' and 'depletion', respectively. However, if sources are demagnified by lensing, then the situation is reversed, and the detectable surface density of galaxies with a shallow source count, as expected at the faintest flux densities, is increased. In general, demagnified sources are difficult to detect and study: exquisite subarcsec angular resolution and surface brightness sensitivity are required, and emission from the lensing object must not dominate the image. These unusual conditions are expected to be satisfied for observations made of the dense swarm of demagnified images that could form very close to the line of sight through the centre of a rich cluster of galaxies using the forthcoming submillimetre-wave Atacama Large Millimeter Array (ALMA) interferometer. The demagnified images of most of the background galaxies lying within about 1 arcmin of a rich cluster of galaxies could be detected in a single 18-arcsec-diameter ALMA field centred on the cluster core, providing an effective increase in the ALMA field of view. This technique could allow a representative sample of faint,  10–100 μJy  submillimetre galaxies to be detected several times more rapidly than in a blank field.     

Confusion limit resulting from galaxies: using the Infrared Array Camera on board SIRTF

Recent ISO data have allowed, for the first time, observationally based estimates for source confusion in mid-infrared surveys. We use the extragalactic source counts from ISOCAM in conjunction with K -band counts to predict the confusion resulting from galaxies in deep mid-infrared observations. We specifically concentrate on the near-future Space Infrared Telescope Facility ( SIRTF ) mission, and calculate expected confusion for the Infrared Array Camera (IRAC) on board SIRTF . A defining scientific goal of the IRAC instrument will be the study of high-redshift galaxies using a deep, confusion-limited wide-field survey at 3–10 μm . A deep survey can reach 3-μJy sources with reasonable confidence in the shorter wavelength IRAC bands. Truly confusion-limited images with the 8 μm will be difficult to obtain because of practical time constraints, unless infrared galaxies exhibit very strong evolution beyond the deepest current observations. We find L * galaxies to be detectable to z =3–3.5 at 8 μm, which is slightly more pessimistic than found in 1999 by Simpson & Eisenhardt.     

Millimetre continuum observations of southern massive star formation regions – II. SCUBA observations of cold cores and the dust grain emissivity index (β)

We report the results of a submillimetre continuum emission survey targeted towards 78 star formation regions, 72 of which are devoid of methanol maser and UC H  ii  regions, identified in the Swedish ESO Submillimetre Telescope (SEST)/SEST IMaging Bolometer Array (SIMBA) millimetre continuum survey of Hill et al. At least 45 per cent of the latter sources, dubbed 'mm-only', detected in this survey are also devoid of the mid-infrared MSX  emission. The 450- and 850-μm continuum emission was mapped using the Submillimetre Common User Bolometer Array (SCUBA) instrument on the James Clerk Maxwell Telescope (JCMT). Emission is detected towards 97 per cent of the 78 sources targeted as well as towards 28 other SIMBA sources lying in the SCUBA fields.
In total, we have identified 212 cores in this submillimetre survey, including 106 previously known from the SIMBA survey. Of the remaining 106 sources, 53 result from resolving a SIMBA source into multiple submillimetre components, whilst the other 53 sources are submillimetre cores, not seen in the SIMBA. Additionally, we have identified two further mm-only sources in the SIMBA images. Of the total 405 sources identified in the SIMBA survey, 255 are only seen at millimetre wavelengths.
We concatenate the results from four (sub)millimetre continuum surveys of massive star formation, together with the Galactic plane map of Pierce-Price et al. in order to determine the dust grain emissivity index β for each of the sources in the SIMBA source list. We examine the value of β with respect to temperature, as well as for the source classes identified in the SIMBA survey, for variation of this index. Our results indicate that β is typically 2, which is consistent with previous determinations in the literature, but for a considerably larger sample than previous work.