Measuring the accretion rate and kinetic luminosity functions of supermassive black holes

We derive accretion rate functions (ARFs) and kinetic luminosity functions (KLFs) for jet-launching supermassive black holes. The accretion rate as well as the kinetic power of an active galaxy is estimated from the radio emission of the jet. For compact low-power jets, we use the core radio emission while the jet power of high-power radio-loud quasars is estimated using the extended low-frequency emission to avoid beaming effects. We find that at low luminosities the ARF derived from the radio emission is in agreement with the measured bolometric luminosity function (BLF) of active galactic nucleus (AGN), i.e. all low-luminosity AGN launch strong jets. We present a simple model, inspired by the analogy between X-ray binaries (XRBs) and AGN, that can reproduce both the measured ARF of jet-emitting sources as well as the BLF. The model suggests that the break in power-law slope of the BLF is due to the inefficient accretion of strongly sub-Eddington sources. As our accretion measure is based on the jet power it also allows us to calculate the KLF and therefore the total kinetic power injected by jets into the ambient medium. We compare this with the kinetic power output from supernova remnants (SNRs) and XRBs, and determine its cosmological evolution.     

A sample of radio-loud active galactic nuclei in the Sloan Digital Sky Survey

A sample of 2712 radio-luminous galaxies is defined from the second data release of the Sloan Digital Sky Survey (SDSS) by cross-comparing the main spectroscopic galaxy sample with two radio surveys: the National Radio Astronomy Observatories (NRAO) Very Large Array (VLA) Sky Survey (NVSS) and the Faint Images of the Radio Sky at Twenty centimeters (FIRST) survey. The comparison is carried out in a multistage process and makes optimal use of both radio surveys by exploiting the sensitivity of the NVSS to extended and multicomponent radio sources in addition to the high angular resolution of the FIRST images. A radio source sample with 95 per cent completeness and 98.9 per cent reliability is achieved, far better than would be possible for this sample if only one of the surveys was used. The radio source sample is then divided into two classes: radio-loud active galactic nuclei (AGN) and galaxies in which the radio emission is dominated by star formation. The division is based on the location of a galaxy in the plane of 4000-Å break strength versus radio luminosity per unit stellar mass and provides a sample of 2215 radio-loud AGN and 497 star-forming galaxies brighter than 5 mJy at 1.4 GHz. A full catalogue of positions and radio properties is provided for these sources. The local radio luminosity function is then derived both for radio-loud AGN and for star-forming galaxies and is found to be in agreement with previous studies. By using the radio to far-infrared (FIR) correlation, the radio luminosity function of star-forming galaxies is also compared to the luminosity function derived in the FIR. It is found to agree well at high luminosities but less so at lower luminosities, confirming that the linearity of the radio to FIR correlation breaks down below about 10²² W Hz⁻¹ at 1.4 GHz.     

MAMBO observations at 240 GHz of optically obscured Spitzer sources: source clumps and radio activity at high redshift

Optically very faint  ( R > 25.5)  sources detected by the Spitzer Space Telescope at 24 μm represent a very interesting population at redshift   z ∼ (1.5–3)  . They exhibit strong clustering properties, implying that they are hosted by very massive haloes, and their mid-infrared emission could be powered by either dust-enshrouded star formation and/or by an obscured active galactic nucleus (AGN). We report observations carried out with the Max Planck Millimetre Bolometer (MAMBO) array at the IRAM 30-m antenna on Pico Veleta of a candidate protocluster with five optically obscured sources selected from the 24-μm Spitzer sample of the First-Look Survey. Interestingly, these sources appear to lie on a high-density filament aligned with the two radio jets of an AGN. Four out of five of the observed sources were detected. We combine these measurements with optical, infrared and radio observations to probe the nature of the candidate protocluster members. Our preliminary conclusions can be summarized as follows: the spectral energy distributions (SEDs) of all sources include both AGN and starburst contributions; the AGN contribution to the bolometric luminosities ranges between 14 and 26 per cent of the total. Such a contribution is enough for the AGN to dominate the emission at 5.8, 8 and 24 μm, while the stellar component, inferred from SED fitting, prevails at 1.25 mm and at  λ < 4.5 μ  m. The present analysis suggests a coherent interplay at high z between extended radio activity and the development of filamentary large-scale structures.     

Exploring the active galactic nucleus and starburst content of local ultraluminous infrared galaxies through 5–8 μm spectroscopy

We present a 5–8 μm analysis of the Spitzer spectra of 71 ultraluminous infrared galaxies (ULIRGs) with redshift   z < 0.15  , devoted to the study of the role of active galactic nuclei (AGN) and starbursts (SBs) as the power source of the extreme infrared emission. Around ∼5 μm, an AGN is much brighter (by a factor of ≈30) than an SB of equal bolometric luminosity. This allows us to detect the presence of even faint accretion-driven cores inside ULIRGs: signatures of AGN activity are found in ∼70 per cent of our sample (50/71 sources). Through a simple analytical model, we are also able to obtain a quantitative estimate of the AGN/SB contribution to the overall energy output of each source. Although the main fraction of ULIRG luminosity is confirmed to arise from star formation events, the AGN contribution is non-negligible (∼23 per cent) and is shown to increase with luminosity. The existence of a rather heterogeneous pattern in the composition and geometrical structure of the dust among ULIRGs is newly supported by the comparison between individual absorption features and continuum extinction.     

Radio sources in the 2dF Galaxy Redshift Survey – II. Local radio luminosity functions for AGN and star-forming galaxies at 1.4 GHz

We have cross-matched the 1.4-GHz NRAO VLA Sky Survey (NVSS) with the first 210 fields observed in the 2dF Galaxy Redshift Survey (2dFGRS), covering an effective area of 325 deg² (about 20 per cent of the final 2dFGRS area). This yields a set of optical spectra of 912 candidate NVSS counterparts, of which we identify 757 as genuine radio identifications – the largest and most homogeneous set of radio source spectra ever obtained. The 2dFGRS radio sources span the redshift range     to 0.438, and are a mixture of active galaxies (60 per cent) and star-forming galaxies (40 per cent). About 25 per cent of the 2dFGRS radio sources are spatially resolved by NVSS, and the sample includes three giant radio galaxies with projected linear size greater than 1 Mpc. The high quality of the 2dF spectra means we can usually distinguish unambiguously between AGN and star-forming galaxies. We make a new determination of the local radio luminosity function at 1.4 GHz for both active and star-forming galaxies, and derive a local star formation density of         .     

Giant radio sources

We present multifrequency Very Large Array (VLA) observations of two giant quasars, 0437−244 and 1025−229, from the Molonglo Complete Sample. These sources have well-defined FR II radio structure, possible one-sided jets, no significant depolarization between 1365 and 4935 MHz and low rotation measure (|RM|<20 rad m⁻²). The giant sources are defined to be those with overall projected size 1 Mpc. We have compiled a sample of about 50 known giant radio sources from the literature, and have compared some of their properties with a complete sample of 3CR radio sources of smaller sizes to investigate the evolution of giant sources, and test their consistency with the unified scheme for radio galaxies and quasars. We find an inverse correlation between the degree of core prominence and total radio luminosity, and show that the giant radio sources have similar core strengths to smaller sources of similar total luminosity. Hence their large sizes are unlikely to be caused by stronger nuclear activity. The degree of collinearity of the giant sources is also similar to that of the sample of smaller sources. The luminosity–size diagram shows that the giant sources are less luminous than our sample of smaller sized 3CR sources, consistent with evolutionary scenarios in which the giants have evolved from the smaller sources, losing energy as they expand to these large dimensions. For the smaller sources, radiative losses resulting from synchrotron radiation are more significant while for the giant sources the equipartition magnetic fields are smaller and inverse Compton loss owing to microwave background radiation is the dominant process. The radio properties of the giant radio galaxies and quasars are consistent with the unified scheme.     

Host galaxy morphologies of X-ray selected AGN: assessing the significance of different black hole fuelling mechanisms to the accretion density of the Universe at z∼ 1.

We use morphological information of X-ray selected active galactic nuclei (AGN) hosts to set limits on the fraction of the accretion density of the Universe at   z ≈ 1  that is not likely to be associated with major mergers. Deep X-ray observations are combined with high-resolution optical data from the Hubble Space Telescope in the All-wavelength Extended Groth strip International Survey, Great Observatories Origins Deep Survey (GOODS) North and GOODS South fields to explore the morphological breakdown of X-ray sources in the redshift interval  0.5 < z < 1.3  . The sample is split into discs, early-type bulge-dominated galaxies, peculiar systems and point sources in which the nuclear source outshines the host galaxy. The X-ray luminosity function and luminosity density of AGN at   z ≈ 1  are then calculated as a function of morphological type. We find that disc-dominated hosts contribute  30 ± 9  per cent to the total AGN space density and  23 ± 6  per cent to the luminosity density at   z ≈ 1  . We argue that AGN in disc galaxies are most likely fuelled not by major merger events but by minor interactions or internal instabilities. We find evidence that these mechanisms may be more efficient in producing luminous AGN     compared to predictions for the stochastic fuelling of massive black holes in disc galaxies.     

High-redshift obscured quasars: radio emission at sub-kiloparsec scales

The radio properties of 11 obscured 'radio-intermediate' quasars at redshifts   z ≳ 2  have been investigated using the European Very-Long-Baseline-Interferometry Network (EVN) at 1.66 GHz. A sensitivity of  ∼25 μJy per 17 × 14 mas² beam  was achieved, and in seven out of 11 sources unresolved radio emission was securely detected. The detected radio emission of each source accounts for ∼30–100 per cent of the total source flux density. The physical extent of this emission is ≲150 pc, and the derived properties indicate that this emission originates from an active galactic nucleus (AGN). The missing flux density is difficult to account for by star formation alone, so radio components associated with jets of physical size ≳150 pc and ≲40 kpc are likely to be present in most of the sources. Amongst the observed sample steep, flat, gigahertz-peaked and compact-steep spectrum sources are all present. Hence, as well as extended and compact jets, examples of beamed jets are also inferred, suggesting that in these sources, the obscuration must be due to dust in the host galaxy, rather than the torus invoked by the unified schemes. Comparing the total to core (≲150 pc) radio luminosities of this sample with different types of AGN suggests that this sample of   z ≳ 2  radio-intermediate obscured quasars shows radio properties that are more similar to those of the high-radio-luminosity end of the low-redshift radio-quiet quasar population than those of Fanaroff–Riley type I (FR I) radio galaxies. This conclusion may reflect intrinsic differences, but could be strongly influenced by the increasing effect of inverse-Compton cooling of extended radio jets at high redshift.     

Hubble Space Telescope imaging survey of sub-mJy star-forming galaxies – I. Morphologies at z ∼ 0.2

We present the first results of our Hubble Space Telescope HST WFPC2 F814W snapshot imaging survey, targeting virtually all sub-mJy decimetric radio-selected star-forming galaxies. The radio selection at ∼1 GHz is free from extinction effects and the radio luminosities are largely unaffected by AGN contamination, making these galaxies ideal tracers of the cosmic star formation history. A subsample of four targets is presented here, selected at 1.4 GHz from the spectroscopically homogenous and complete samples of Benn et al. and Hopkins et al. The redshifts are confined to a narrow range around z ∼0.2, to avoid differential evolution, with a radio luminosity close to L ∗ where the galaxies dominate the comoving volume-averaged star formation rate. We find clearly disturbed morphologies resembling those of ultraluminous infrared galaxies, indicating that galaxy interactions may be the dominant mechanism for triggering star formation at these epochs. The morphologies are also clearly different from those of coeval quasars and radio galaxies, as found in star-forming galaxies selected at other wavelengths. This may prove challenging for models that propose direct causal links between AGN evolution and the cosmic star formation history at these epochs. The asymmetries are typically much larger than seen in the Canada–France Redshift Survey at similar redshifts, optical luminosities and H α -derived star formation rates, indicating the possible existence of an obscuration-related morphological bias in such samples.     

A unifying view of the spectral energy distributions of blazars

We collect data at well-sampled frequencies from the radio to the γ-ray range for the following three complete samples of blazars: the Slew survey, the 1-Jy samples of BL Lacs and the 2-Jy sample of flat-spectrum radio-loud quasars (FSRQs). The fraction of objects detected in γ-rays ( E  ≳ 100 MeV) is ∼ 17, 26 and 40 per cent in the three samples respectively. Except for the Slew survey sample, γ-ray detected sources do not differ either from other sources in each sample, or from all the γ-ray detected sources, in terms of the distributions of redshift, radio and X-ray luminosities or of the broad-band spectral indices (radio to optical and radio to X-ray). We compute average spectral energy distributions (SEDs) from radio to γ-rays for each complete sample and for groups of blazars binned according to radio luminosity, irrespective of the original classification as BL Lac or FSRQ. The resulting SEDs show a remarkable continuity in that (i) the first peak occurs in different frequency ranges for different samples/luminosity classes, with most luminous sources peaking at lower frequencies; (ii) the peak frequency of the γ-ray component correlates with the peak frequency of the lower energy one; (iii) the luminosity ratio between the high and low frequency components increases with bolometric luminosity. The continuity of properties among different classes of sources and the systematic trends of the SEDs as a function of luminosity favour a unified view of the blazar phenomenon: a single parameter, related to luminosity, seems to govern the physical properties and radiation mechanisms in the relativistic jets present in BL Lac objects as well as in FSRQs. The general implications of this unified scheme are discussed while a detailed theoretical analysis, based on fitting continuum models to the individual spectra of most γ-ray blazars, is presented in a separate paper.     

Cooling flows as a calorimeter of active galactic nucleus mechanical power

The assumption that radiative cooling of gas in the centres of galaxy clusters is approximately balanced by energy input from a central supermassive black hole implies that the observed X-ray luminosity of the cooling flow region sets a lower limit on active galactic nucleus (AGN) mechanical power. The conversion efficiency of the mechanical power of the AGN into gas heating is uncertain, but we argue that it can be high even in the absence of strong shocks. These arguments inevitably lead to the conclusion that the time-averaged mechanical power of AGNs in cooling flows is much higher than the bolometric luminosity of these objects observed currently.
The energy balance between cooling losses and AGN mechanical power requires some feedback mechanism. We consider a toy model in which the accretion rate on to a black hole is set by the classic Bondi formula. Application of this model to the best studied case of M87 suggests that accretion proceeds at approximately the Bondi rate down to a few gravitational radii with most of the power (at the level of a few per cent of the rest mass) being carried away by an outflow.     

A MERLIN neutral hydrogen absorption study of the luminous infrared merger NGC 6240

We present neutral hydrogen absorption observations of the luminous infrared merger NGC 6240 using MERLIN with a resolution of 0.2 arcsec. Broad absorption (a few hundred km s⁻¹) has been found against two compact radio sources within the central kpc providing dynamical information about the neutral gas components in front of these sources. A narrow absorption component is also detected superimposed upon this broad absorption and additionally against some of the extended L -band continuum. From these results we deduce that the broad component is a result of absorption by a highly disturbed disc-like structure of neutral gas aligned along the position angle of the two compact radio sources, similar to the model previously proposed by Tacconi et al. at the end of the last century based on spectral CO emission data. The narrow component is likely to arise from absorption by less disturbed neutral gas at much larger scales within the system.
Continuum observations presented here at 1.4 and 5 GHz support the view that NGC 6240 contains a double nucleus resulting from a galactic merger event and show these as two compact radio sources separated by 1.52 arcsec. We have also applied luminosity and morphological considerations to the continuum results to determine the most feasible source of radio emission for this luminous merger galaxy. We conclude that the most likely source of the radio flux found in NGC 6240 is a combination of starburst emission from radio luminous supernova remnants, similar to those found in Arp 220, and emission from a weak AGN probably triggered by a merger event.     

On the redshift cut-off for steep-spectrum radio sources

We use three samples (3CRR, 6CE and 6C*) selected at low radio frequency to constrain the cosmic evolution in the radio luminosity function (RLF) for the 'most luminous' steep-spectrum radio sources. Though intrinsically rare, such sources give the largest possible baseline in redshift for the complete flux-density-limited samples currently available. Using parametric models to describe the RLF, incorporating distributions in radio spectral shape and linear size, as well as the usual luminosity and redshift, we find that the data are consistent with a constant comoving space density between     and     . We find that this model is favoured over a model with similar evolutionary behaviour to that of optically selected quasars (i.e., a roughly Gaussian distribution in redshift) with probability ratios of     and     for spatially flat cosmologies with     and     respectively. Within the uncertainties, this evolutionary behaviour may be reconciled with the shallow decline preferred for the comoving space density of flat-spectrum sources by Dunlop & Peacock and Jarvis & Rawlings, in line with the expectations of unified schemes.     

Evidence of a link between the evolution of clusters and their AGN fraction

We discuss the optical properties, X-ray detections and active galactic nucleus (AGN) populations of four clusters at   z ∼ 1  in the Subaru–XMM Deep Field (SXDF). The velocity distribution and plausible extended X-ray detections are examined, as well as the number of X-ray point sources and radio sources associated with the clusters. We find that the two clusters that appear virialized and have an extended X-ray detection contain few, if any, AGN, whereas the two pre-virialized clusters have a large AGN population. This constitutes evidence that the AGN fraction in clusters is linked to the clusters' evolutionary stage. The number of X-ray AGN in the pre-virialized clusters is consistent with an overdensity of a factor of ∼200; the radio AGN appear to be clustered with a factor of 3 to 6 higher. The median K -band luminosities of   L_K = 1.7 ± 0.7 L *  for the X-ray sources and   L_K = 2.3 ± 0.1 L *  for the radio sources support the theory that these AGN are triggered by galaxy interaction and merging events in sub-groups with low internal velocity distributions, which make up the cluster environment in a pre-virialization evolutionary stage.     

XMM–Newton surveys of the Canada–France Redshift Survey Fields – III. The environments of X-ray selected active galactic nuclei at 0.4 < z < 0.6

The environmental properties of a sample of 31 hard X-ray selected active galactic nuclei (AGN) are investigated, from scales of 500 kpc down to 30 kpc, and are compared to a control sample of inactive galaxies. All the AGN lie in the redshift range  0.4 < z < 0.6  . The accretion luminosity density of the Universe peaks close to this redshift range, and the AGN in the sample have X-ray luminosities close to the knee in the hard X-ray luminosity function, making them representative of the population that dominated this important phase of energy conversion.
Using both the spatial clustering amplitude and near-neighbour counts, it is found that the AGN have environments that are indistinguishable from normal, inactive galaxies over the same redshift range and with similar optical properties. Typically, the environments are of subcluster richness, in contrast to similar studies of high- z quasars, which are often found in clusters with comparable richness to the Abell   R ≥ 0  clusters.
It is suggested that minor mergers with low-mass companions are a likely candidate for the mechanism by which these modest luminosity AGN are fuelled.     

The bulk kinetic power of radio jets in active galactic nuclei

Based on the Königl's inhomogeneous jet model, we estimate the jet parameters, such as bulk Lorentz factor Γ, viewing angle θ and electron number density n _e from radio very long-baseline interferometry and X-ray data for a sample of active galactic nuclei (AGNs) assuming that the X-rays are from the jet rather than the intracluster gas. The bulk kinetic power of jets is then calculated using the derived jet parameters. We find a strong correlation between the total luminosity of broad emission lines and the bulk kinetic power of the jets. This result supports the scenario that the accretion process is tightly linked with the radio jets, though how the disc and jet are coupled is not revealed by present correlation analysis. Moreover, we find a significant correlation between the bulk kinetic power and radio extended luminosity. This implies that the emission from the radio lobes is closely related with the energy flux transported through jets from the central part of AGNs.     

Doppler boosting, superluminal motion, and the kinematics of AGN jets

We discuss results from a decade long program to study the fine-scale structure and the kinematics of relativistic AGN jets with the aim of better understanding the acceleration and collimation of the relativistic plasma forming AGN jets. From the observed distribution of brightness temperature, apparent velocity, flux density, time variability, and apparent luminosity, the intrinsic properties of the jets including Lorentz factor, luminosity, orientation, and brightness temperature are discussed. Special attention is given to the jet in M87, which has been studied over a wide range of wavelengths and which, due to its proximity, is observed with excellent spatial resolution. Most radio jets appear quite linear, but we also observe curved non-linear jets and non-radial motions. Sometimes, different features in a given jet appear to follow the same curved path but there is evidence for ballistic trajectories as well. The data are best fit with a distribution of Lorentz factors extending up to γ∼30 and intrinsic luminosity up to ∼10²⁶ W Hz⁻¹. In general, gamma-ray quasars may have somewhat larger Lorentz factors than non gamma-ray quasars. Initially the observed brightness temperature near the base of the jet extend up to ∼5×10¹³ K which is well in excess of the inverse Compton limit and corresponds to a large excess of particle energy over magnetic energy. However, more typically, the observed brightness temperatures are ∼2×10¹¹ K, i.e., closer to equipartition.     

A diffuse bubble-like radio-halo source MRC 0116+111: imprint of AGN feedback in a low-mass cluster of galaxies

We present detailed observations of MRC 0116+111, revealing a luminous, miniradio halo of ∼240-kpc diameter located at the centre of a cluster of galaxies at redshift   z = 0.131  . Our optical and multiwavelength Giant Metrewave Radio Telescope and Very Large Array radio observations reveal a highly unusual radio source: showing a pair of giant (∼100-kpc diameter) bubble-like diffuse structures, that are about three times larger than the analogous extended radio emission observed in M87 – the dominant central radio galaxy in the Virgo cluster. However, in MRC 0116+111 we do not detect any ongoing active galactic nucleus (AGN) activity, such as a compact core or active radio jets feeding the plasma bubbles. The radio emitting relativistic particles and magnetic fields were probably seeded in the past by a pair of radio jets originating in the AGN of the central cD galaxy. The extremely steep high-frequency radio spectrum of the north-western bubble, located ∼100 kpc from cluster centre, indicates radiation losses, possibly because having detached, it is rising buoyantly and moving away into the putative hot intracluster medium. The other bubble, closer to the cluster centre, shows signs of ongoing particle re-acceleration. We estimate that the radio jets which inflated these two bubbles might have also fed enough energy into the intracluster medium to create an enormous system of cavities and shock fronts, and to drive a massive outflow from the AGN, which could counter-balance and even quench a cooling flow. Therefore, this source presents an excellent opportunity to understand the energetics and the dynamical evolution of radio jet inflated plasma bubbles in the hot cluster atmosphere.     

Radio galaxies with a 'double-double' morphology – II. The evolution of double-double radio galaxies and implications for the alignment effect in FRII sources

A double-double radio galaxy (DDRG) is defined as consisting of a pair of double radio sources with a common centre. In this paper we present an analytical model in which the peculiar radio structure of DDRGs is caused by an interruption of the jet flow in the central AGN. The new jets emerging from the restarted AGN give rise to an inner source structure within the region of the old, outer cocoon. Standard models of the evolution of FRII sources predict gas densities within the region of the old cocoon that are insufficient to explain the observed properties of the inner source structure. Therefore additional material must have passed from the environment of the source through the bow shock surrounding the outer source structure into the cocoon. We propose that this material is warm clouds (∼10⁴ K) of gas embedded in the hot IGM which are eventually dispersed over the cocoon volume by surface instabilities induced by the passage of cocoon material. The derived lower limits for the volume filling factors of these clouds are in good agreement with results obtained from optical observations. The long time-scales for the dispersion of the clouds (∼10⁷ yr) are consistent with the apparently exclusive occurrence of the DDRG phenomenon in large (≳700 kpc) radio sources, and with the observed correlation of the strength of the optical/UV alignment effect in z ∼1 FRII sources with their linear size.