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.     

The clustering of hot and cold IRAS galaxies: the redshift-space correlation function

We measure the autocorrelation function, ξ , of galaxies in the IRAS Point Source Catalogue galaxy redshift (PSC z ) survey and investigate its dependence on the far-infrared colour and absolute luminosity of the galaxies. We find that the PSC z survey correlation function can be modelled out to a scale of 10  h ⁻¹ Mpc as a power law of slope 1.30±0.04 and correlation length 4.77±0.20 . At a scale of 75  h ⁻¹ Mpc we find the value of J ₃ to be 1500±400 .
We also find that galaxies with higher 100 μm/60 μm flux ratio, corresponding to cooler dust temperatures, are more strongly clustered than warmer galaxies. Splitting the survey into three colour subsamples, we find that, between 1 and 10  h ⁻¹ Mpc, the ratio of ξ is a factor of 1.5 higher for the cooler galaxies compared with the hotter galaxies. This is consistent with the suggestion that hotter galaxies have higher star formation rates, and correspond to later-type galaxies which are less clustered than earlier types.
Using volume-limited subsamples, we find a weak variation of ξ as a function of absolute luminosity, in the sense that more luminous galaxies are less clustered than fainter galaxies. The trend is consistent with the colour dependence of ξ and the observed colour–luminosity correlation, but the large uncertainties mean that it has a low statistical significance.     

Evolution of the luminosity function and colours of galaxies in a Λ cold dark matter universe

The luminosity function of galaxies is derived from a cosmological hydrodynamic simulation of a Λ cold dark matter universe with the aid of a stellar population synthesis model. At     , the resulting B -band luminosity function has a flat faint-end slope of     with the characteristic luminosity and the normalization in fair agreement with observations, while the dark matter halo mass function is steep with a slope of     . The colour distribution of galaxies also agrees well with local observations. We also discuss the evolution of the luminosity function, and the colour distribution of galaxies from     to 5. A large evolution of the characteristic mass in the stellar mass function as a result of number evolution is compensated by luminosity evolution; the characteristic luminosity increases only by 0.8 mag from     to 2, and then declines towards higher redshift, while the B -band luminosity density continues to increase from     to 5 (but only slowly at     .     

Evolution of the u-band luminosity function from redshift 1.2 to 0

We produce and analyse u -band (  λ≈ 355  nm) luminosity functions (LFs) for the red and blue populations of galaxies using data from the Sloan Digital Sky Survey (SDSS) u -band Galaxy Survey ( u GS) and Deep Evolutionary Exploratory Probe 2 (DEEP2) survey. From a spectroscopic sample of 41 575 SDSS u GS galaxies and 24 561 DEEP2 galaxies, we produce colour magnitude diagrams and make use of the colour bimodality of galaxies to separate red and blue populations. LFs for eight redshift slices in the range  0.01 < z < 1.2  are determined using the  1/ V _max  method and fitted with Schechter functions showing that there is significant evolution in   M *  , with a brightening of 1.4 mag for the combined population. The integration of the Schechter functions yields the evolution in the u -band luminosity density (LD) out to   z ∼ 1  . By parametrizing the evolution as  ρ∝ (1 + z )^β  , we find that  β= 1.36 ± 0.2  for the combined populations and  β= 2.09 ± 0.2  for the blue population. By removing the contribution of the old stellar population to the u -band LD and correcting for dust attenuation, we estimate the evolution in the star formation rate (SFR) of the Universe to be  β_SFR= 2.5 ± 0.3  . Discrepancies between our result and higher evolution rates measured using the infrared and far-UV can be reconciled by considering possibilities such as an underestimated dust correction at high redshifts or evolution in the stellar initial mass function.     

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 dwarf galaxy population of the Coma cluster to MR=−11: a detailed description

We present the luminosity function and measurements of the scalelengths, colours and radial distribution of dwarf galaxies in the Coma cluster down to R =24. Our survey area is 674 arcmin²; this is the deepest and most detailed survey covering such a large area.
Our measurements agree with those of most previous authors at bright and intermediate magnitudes. The new results are as follows.
(1) Galaxies in the Coma cluster have a luminosity function φ( L )∝ L ^α that is steep (α∼−1.7) for −15< M_R <−11, and is shallower brighter than this. The curvature in the luminosity function at M_R ∼−15 is statistically significant.
(2) The galaxies that contribute most strongly to the luminosity function at −14< M_R <−12 have colours and scalelengths that are consistent with those of local dwarf spheroidal galaxies placed at the distance of Coma.
(3) These galaxies with −14< M_R <−12 have a colour distribution that is very strongly peaked at B − R =1.3. This is suggestive of a substantial degree of homogeneity in their star formation histories and metallicities.
(4) These galaxies with −14< M_R <−12 also appear to be more confined to the cluster core ( r ∼200 kpc) than the brighter galaxies. Alternatively, this observation may be explained in part or whole by the presence of an anomalously high number of background galaxies behind the cluster core. Velocity measurements of these galaxies would distinguish between these two possibilities.     

The luminosity function and surface brightness distribution of H i selected galaxies

We measure the     B -band optical luminosity function (LF) for galaxies selected in a blind H  i survey. The total LF of the H  i selected sample is flat, with Schechter parameters     and     , in good agreement with LFs of optically selected late-type galaxies. Bivariate distribution functions of several galaxy parameters show that the H  i density in the local Universe is more widely spread over galaxies of different size, central surface brightness and luminosity than the optical luminosity density is. The number density of very low surface brightness (LSB ) (>24.0 mag arcsec⁻²) gas-rich galaxies is considerably lower than that found in optical surveys designed to detect dim galaxies. This suggests that only a part of the population of LSB galaxies is gas-rich and that the rest must be gas-poor. However, we show that this gas-poor population must be cosmologically insignificant in baryon content. The contribution of gas-rich LSB galaxies (>23.0 mag arcsec⁻²) to the local cosmological gas and luminosity density is modest     and     per cent respectively); their contribution to Ω_matter is not well-determined, but probably <11 per cent. These values are in excellent agreement with the low redshift results from the Hubble Deep Field.     

The redshift-space two-point correlation function of ELAIS-S1 galaxies

We investigate the clustering properties of galaxies in the recently completed ELAIS-S1 redshift survey through their spatial two-point autocorrelation function. We used a subsample of the ELAIS-S1 catalogue covering approximately 4 deg² and consisting of 148 objects selected at 15 μm with a flux >0.5 mJy and a redshift   z < 0.5  . We detected a positive signal in the correlation function that in the range of separations  1–10  h ⁻¹ Mpc  is well approximated by a power law with a slope  γ= 1.4 ± 0.25  and a correlation length   s ₀= 5.4 ± 1.2  h ⁻¹ Mpc  , at the 90 per cent significance level. This result is in good agreement with the redshift-space correlation function measured in more local samples of mid-infrared-selected galaxies such as the IRAS Point Source Catalog (PSC z ) redshift survey. This suggests a lack of significant clustering evolution of infrared-selected objects out to   z = 0.5  that is further confirmed by the consistency found between the correlation functions measured in a local  ( z < 0.2)  and a distant  (0.2 < z < 0.5)  subsample of ELAIS-S1 galaxies. We also confirm that optically selected galaxies in the local redshift surveys, especially those of the SDSS sample, are significantly more clustered than infrared objects.     

The X-ray luminosity function of AGN at z∼ 3

We combine Lyman-break colour selection with ultradeep (≳200 ks) Chandra X-ray imaging over a survey area of ∼0.35 deg² to select high-redshift active galactic nuclei (AGN). Applying careful corrections for both the optical and X-ray selection functions, the data allow us to make the most accurate determination to date of the faint end of the X-ray luminosity function (XLF) at   z ∼ 3  . Our methodology recovers a number density of X-ray sources at this redshift which is at least as high as previous surveys, demonstrating that it is an effective way of selecting high z AGN. Comparing to results at   z = 1  , we find no evidence that the faint slope of the XLF flattens at high z , but we do find significant (factor ∼3.6) negative evolution of the space density of low luminosity AGN. Combining with bright end data from very wide surveys we also see marginal evidence for continued positive evolution of the characteristic break luminosity   L _*  . Our data therefore support models of luminosity-dependent density evolution between   z = 1  and   z = 3  . A sharp upturn in the the XLF is seen at the very lowest luminosities  ( L _X≲ 10^42.5 erg s⁻¹)  , most likely due to the contribution of pure X-ray starburst galaxies at very faint fluxes.     

The infrared luminosity of the torus and the visibility of scattered broad line emission in Seyfert 2 galaxies

A number of studies have shown that the visibility of scattered broad emission lines in Seyfert 2 galaxies is strongly dependent on the IRAS     flux ratio, where those Seyfert 2 galaxies with 'warm' IRAS colours show polarized broad line emission. It is now clear that this effect is owing to the increasing dominance of the galactic rather than the active galactic nucleus (AGN) emission at 60 μm in less-luminous 'cool' Seyfert 2 galaxies. However, we present evidence that the 25-μm emission is a good measure of the AGN luminosity for most Seyfert 2 galaxies. Using this result, we show that the visibility of scattered broad line emission has a dependence on the AGN luminosity. The observations can be interpreted self-consistently if the scaleheight of the scattering zone varies with central source luminosity whilst the scaleheight of the obscuring torus is approximately constant.     

CENSORS: A Combined EIS–NVSS Survey Of Radio Sources – I. Sample definition, radio data and optical identifications

A new sample of radio sources, with the designated name CENSORS (A Combined EIS–NVSS Survey Of Radio Sources), has been defined by combining the National Radio Astronomy Observatory Very Large Array Sky Survey (NVSS) at 1.4 GHz with the ESO Imaging Survey (EIS) Patch D, a 3° by 2° region of sky centred at RA     , Dec. −21°00'00' (J2000). New radio observations of 199 NVSS radio sources with NVSS flux densities   S _{1.4 GHz} > 7.8 mJy  are presented, and are compared with the EIS I -band imaging observations which reach a depth of   I ∼ 23  ; optical identifications are obtained for over two-thirds of the ∼150 confirmed radio sources within the EIS field. The radio sources have a median linear size of 6 arcsec, consistent with the trend for lower flux density radio sources to be less extended. Other radio source properties, such as the lobe flux density ratios, are consistent with those of brighter radio source samples. From the optical information, 30–40 per cent of the sources are expected to lie at redshifts   z ≳ 1.5  .
One of the key goals of this survey is to accurately determine the high-redshift evolution of the radio luminosity function. These radio sources are at the ideal flux density level to achieve this goal; at redshifts   z ∼ 2  they have luminosities which are around the break of the luminosity function and so provide a much more accurate census of the radio source population at those redshifts than the existing studies of extreme, high radio power sources. Other survey goals include investigating the dual-population unification schemes for radio sources, studying the radio luminosity dependence of the evolution of radio source environments, and understanding the radio power dependence of the K – z relation for radio galaxies.     

A ROSAT survey of Wolf–Rayet galaxies — II. The extended sample

We present results from an ongoing X-ray survey of Wolf–Rayet (WR) galaxies, a class of objects believed to be very young starbursts. This paper extends the first X-ray survey of WR galaxies by Stevens &38; Strickland by studying WR galaxies identified subsequent to the original WR galaxy catalogue of Conti.   Out of a sample of 40 new WR galaxies a total of 10 have been observed with the ROSAT PSPC, and of these seven have been detected (NGC 1365, NGC 1569, I Zw 18, NGC 3353, NGC 4449, NGC 5408 and a marginal detection of NGC 2366). Of these, all are dwarf starbursts except for NGC 1365, which is a barred spiral galaxy possibly with an active nucleus. We also report on observations of the related emission-line galaxy IRAS 0833+6517.   The X-ray properties of these galaxies are broadly in line with those found for the original sample; they are X-ray overluminous compared with their blue luminosity and have thermal spectra with typically kT  ∼ 0.4 − 1.0 keV. There are some oddities: NGC 5408 is very overluminous in X-rays, even compared with other WR galaxies; I Zw 18 has a harder X-ray spectrum; NGC 1365, although thought to contain an active nucleus, has X-ray properties that are broadly similar to other WR galaxies, and we suggest that the X-ray emission from NGC 1365 is due to starburst activity.   A good correlation between X-ray and blue luminosity is found for the WR galaxy sample as a whole. However, when just dwarf galaxies are considered there is little evidence of correlation. We discuss the implications of these results on our understanding of the X-ray emission from WR galaxies and suggest that the best explanation for the X-ray activity is starburst activity from a young starburst region.     

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/     

Statistical analysis of galaxy surveys – III. The non-linear clustering of red and blue galaxies in the 2dFGRS

We present measurements of the higher order clustering of red and blue galaxies as a function of scale and luminosity made from the two-degree field galaxy redshift survey (2dFGRS). We use a counts-in-cells analysis to estimate the volume-averaged correlation functions,     , as a function of scale up to the order of   p = 5  , and also the reduced void probability function. Hierarchical amplitudes are constructed using the estimates of the correlation functions:     . We find that (i) red galaxies display stronger clustering than blue galaxies at all orders measured; (ii) red galaxies show values of   S _p   that are strongly dependent on luminosity whereas blue galaxies show no segregation in   S _p   within the errors; this is remarkable given the segregation in the variance; (iii) the linear relative bias shows the opposite trend to the hierarchical amplitudes, with little segregation for the red sequence and some segregation for the blue; (iv) faint red galaxies deviate significantly from the 'universal' negative binomial reduced void probabilities followed by all other galaxy populations. Our results show that the characteristic colour of a galaxy population reveals a unique signature in its spatial distribution. Such signatures will hopefully further elucidate the physics responsible for shaping the cosmological evolution of galaxies.     

Massive dark matter haloes around bright isolated galaxies in the 2dFGRS

We identify a large sample of isolated bright galaxies and their fainter satellites in the 2dF Galaxy Redshift Survey (2dFGRS). We analyse the dynamics of ensembles of these galaxies selected according to luminosity and morphological type by stacking the positions of their satellites and estimating the velocity dispersion of the combined set. We test our methodology using realistic mock catalogues constructed from cosmological simulations. The method returns an unbiased estimate of the velocity dispersion provided that the isolation criterion is strict enough to avoid contamination and that the scatter in halo mass at fixed primary luminosity is small. Using a maximum likelihood estimator that accounts for interlopers, we determine the satellite velocity dispersion within a projected radius of 175  h ⁻¹ kpc. The dispersion increases with the luminosity of the primary and is larger for elliptical galaxies than for spiral galaxies of similar b _J luminosity. Calibrating the mass–velocity dispersion relation using our mock catalogues, we find a dynamical mass within 175  h ⁻¹ kpc of     for elliptical galaxies and     for spiral galaxies. Finally, we compare our results with recent studies and investigate their limitations using our mock catalogues.     

Modelling the number counts of early-type galaxies by pure luminosity evolution

In this paper, we explore the plausible luminosity evolution of early-type galaxies in different cosmological models by constructing a set of pure luminosity evolution (PLE) models via the choices of the star-formation rate (SFR) parameters and formation redshift z _f of galaxies, with the observational constraints derived from the Hubble Space Telescope ( HST  ) morphological number counts for elliptical and S0 galaxies of the Medium Deep Survey (MDS) and the Hubble Deep Field (HDF). We find that the number counts of early-type galaxies can be explained by the pure luminosity evolution models, without invoking exotic scenarios such as merging or introducing an additional population, but the evolution should be nearly passive, with a high z _f assumed. The conclusion is valid in all of the three cosmological models we adopt in this paper. We also present the redshift distributions for three bins of observed magnitudes in the F814w passband, to show the redshift at which the objects that dominate the counts at a given magnitude may be found. The predictions of the redshift distribution of 22.5 <  b _j  < 24.0 are also presented for comparison with future data.     

Radio galaxies in the 2SLAQ Luminous Red Galaxy survey – II. The stellar populations of radio-loud and radio-quiet LRGs

We present an analysis of the optical spectra of a volume-limited sample of 375 radio galaxies at redshift  0.4 < z < 0.7  from the 2dF-SDSS (Sloan Digital Sky Survey) Luminous Red Galaxy (LRG) and QSO (quasi-stellar object) (2SLAQ) redshift survey. We investigate the evolution of the stellar populations and emission-line properties of these galaxies. By constructing composite spectra and comparing with a matched sample of radio-quiet sources from the same survey, we also investigate the effect on the galaxy of the presence of an active nucleus.
The composite spectra, binned by redshift and radio luminosity, all require two components to describe them, which we interpret as an old and a younger population. We found no evolution with redshift of the age of the younger population in radio galaxies, nor were they different from the radio-quiet comparison sample. Similarly, there is no correlation with radio power, with the exception that the most powerful radio sources  ( P _1.4 > 10²⁶  W Hz⁻¹) have younger stars and stronger emission lines than the less powerful sources. This suggests that we have located the threshold in radio power where strong emission lines 'switch on', at radio powers of around 10²⁶ W Hz⁻¹. Except for the very powerful radio galaxies, the presence of a currently active radio active galactic nucleus (AGN) does not appear to be correlated with any change in the observed stellar population of a luminous red galaxy at   z ∼ 0.5  .     

The European Large Area ISO Survey – III. 90-μm extragalactic source counts

We present results and source counts at 90 μm extracted from the preliminary analysis of the European Large Area ISO Survey (ELAIS). The survey covered about 12 deg² of the sky in four main areas and was carried out with the ISOPHOT instrument onboard the Infrared Space Observatory ( ISO ). The survey is at least an order of magnitude deeper than the IRAS 100-μm survey and is expected to provide constraints on the formation and evolution of galaxies. The majority of the detected sources are associated with galaxies on optical images. In some cases the optical associations are interacting pairs or small groups of galaxies, suggesting that the sample may include a significant fraction of luminous infrared galaxies. The source counts extracted from a reliable subset of the detected sources are in agreement with strongly evolving models of the starburst galaxy population.