Astronomy

Radio loud AGN restarting as Radio Quiet AGN?

Radio loud AGN restarting as Radio Quiet AGN?


We are searching data for your request:

Forums and discussions:
Manuals and reference books:
Data from registers:
Wait the end of the search in all databases.
Upon completion, a link will appear to access the found materials.

Is it possible for a Radio-loud object (eg., an FRI galaxy) to evolve into a Radio-quiet (eg., a Seyfert) source?

Also, can this happen via AGN restart? That is, have there been instances where the active nuclei in a Radio-loud object switches off, and then restarts as a weaker Radio-quiet source?


Quasars That Have Transitioned from Radio-quiet to Radio-loud on Decadal Timescales Revealed by VLASS and FIRST

Kristina Nyland, Dillon Z. Dong, Pallavi Patil, Mark Lacy, Sjoert Van Velzen, Amy E. Kimball, Sumit K. Sarbadhicary, Gregg Hallinan, Vivienne Baldassare, Tracy E. Clarke, Andy D. Goulding, Jenny Greene , Andrew Hughes, Namir Kassim, Magdalena Kunert-Bajraszewska, Thomas J. MacCarone, Kunal Mooley, Dipanjan Mukherjee, Wendy Peters, Leonid Petrov Show 4 others Show less Emil Polisensky, Wiphu Rujopakarn, Mark Whittle, Mattia Vaccari

Research output : Contribution to journal › Article › peer-review


Title: The Origin of Radio Emission from Radio-Quiet AGN

10^5 across the different AGN classes. The weakest radio sources, radio-quiet (RQ) AGN, are typically 1,000 times fainter than the radio-loud (RL) AGN, and represent the majority of the AGN population. In RL AGN, radio emission is essentially all produced by synchrotron emission from a relativistic jet. In contrast, in RQ AGN the absence of luminous jets allows us to probe radio emission from a wide range of possible mechanisms, from the host galaxy kpc scale down to the innermost region near the SMBHs: star formation, AGN driven wind, free-free emission from photo-ionized gas, low power jet, and the innermost accretion disc coronal activity. All these mechanisms can now be probed with unprecedented precision and spatial resolution, thanks to the current and forthcoming generation of highly sensitive radio arrays.


1. Introduction

Active galactic nuclei (AGNs) are powered by accretion on a supermassive black hole (BH). Although the first discovered AGNs were radio-loud (RL), most of AGNs are radio-quiet (RQ). The radio-loudness of RL AGNs, defined as the ratio of the radio flux to the optical flux (Kellermann et al., 1989), covers several orders of magnitude (e.g., Sikora et al., 2007 Lal and Ho, 2010) which implies very diverse jet production efficiencies. For jets powered by rotating BHs (Blandford and Znajek, 1977), the efficiency of the jet production can be related to the spread of the BH spins and the amount of magnetic fluxes accumulated in the central AGNs.

Magnetic fluxes can be developed stochastically in the innermost zones of accretion discs (Begelman and Armitage, 2014), or can be advected to the central regions of a galaxy prior to the AGN phase (Sikora et al., 2013 Sikora and Begelman, 2013). In the latter case there could be systematic differences between the properties of galaxies hosting radio-loud and radio-quiet AGNs. Our aim is to compare the properties of the host galaxies of radio-loud and radio-quiet AGNs to distinguish which of these two scenarios is more probable. If we find that RQ and RL are hosted by different galaxies, we could discard the scenario where radio jets are stochastic.

Radio-loud (jetted) and radio-quiet (non-jetted) AGNs have already been studied extensively and some systematic differences were found between these two groups of objects:

• The fraction of radio-loud objects and radio loudness increases with decreasing Eddington ratio, λ = Lbol/LEdd (e.g., Terashima and Wilson, 2003 Kratzer and Richards, 2015), but there is a large scatter in radio loudness for AGNs with similar Eddington ratio.

• The fraction of galaxies with disturbed morphology is larger in RL than in RQ AGNs (Bessiere et al., 2012 Chiaberge et al., 2015).

However, the differences listed above concern entire families of objects, but not objects that have the same accretion properties. Moreover, in catalogs of AGNs there are many objects with properties characteristic of RL objects like very massive black holes, low Eddington ratios, and disturbed morphologies, but they are radio-quiet. Therefore, we ask the question why is the efficiency of the jet production very different among otherwise similar objects?

In our project we concentrate on Type 2 (i.e., obscured) objects, to be able to study the properties of the host galaxies, with Eddington ratios λ ≥ 0.003, and we seek the differences between radio-loud and radio-quiet AGNs to check if there are any generic differences between the host galaxies of these two groups of objects that can indicate different evolution histories of jetted and non-jetted AGNs.


Radio AGN in the local universe: unification, triggering and evolution

Associated with one of the most important forms of active galactic nucleus (AGN) feedback, and showing a strong preference for giant elliptical host galaxies, radio AGN ( (L_<1.4,mathrm> > 10^<24>) W (hbox ^<-1>) ) are a key sub-class of the overall AGN population. Recently their study has benefitted dramatically from the availability of high-quality data covering the X-ray to far-IR wavelength range obtained with the current generation of ground- and space-based telescope facilities. Reflecting this progress, here I review our current state of understanding of the population of radio AGN at low and intermediate redshifts ( (z < 0.7) ), concentrating on their nuclear AGN and host galaxy properties, and covering three interlocking themes: the classification of radio AGN and its interpretation the triggering and fuelling of the jet and AGN activity and the evolution of the host galaxies. I show that much of the observed diversity in the AGN properties of radio AGN can be explained in terms of a combination of orientation/anisotropy, mass accretion rate, and variability effects. The detailed morphologies of the host galaxies are consistent with the triggering of strong-line radio galaxies (SLRG) in galaxy mergers. However, the star formation properties and cool ISM contents suggest that the triggering mergers are relatively minor in terms of their gas masses in most cases, and would not lead to major growth of the supermassive black holes and stellar bulges therefore, apart from a minority (<20 %) that show evidence for higher star formation rates and more massive cool ISM reservoirs, the SLRG represent late-time re-triggering of activity in mature giant elliptical galaxies. In contrast, the host and environmental properties of weak-line radio galaxies (WLRG) with Fanaroff–Riley class I radio morphologies are consistent with more gradual fuelling of the activity via gas accretion at low rates onto the supermassive black holes.

This is a preview of subscription content, access via your institution.


Radio loud AGN restarting as Radio Quiet AGN? - Astronomy

Before considering their detailed properties, it is important to start by defining radio AGN as a class. One approach is to use the shapes of the AGN spectral energy distributions (SEDs), for example defining radio AGN to have a minimum ratio of radio to optical luminosity (e.g. R = L5 GHz / LB > 10: Kellermann et al., 1989). However, the disadvantage of this method is that the optical AGN luminosity estimates may be affected by dust obscuration — a particular problem for narrow-line AGN — and contamination by the direct starlight of the host galaxies. There is also evidence that the radio-to-optical luminosity ratio that defines the boundary between radio-loud and radio-quiet AGN increases substantially towards low AGN luminosities and Eddington ratios (Sikora, Stawarz & Lasota, 2007, Chiaberge & Marconi, 2011). Use of a mid-IR-to-radio flux ratio (e.g. q24 = L24 µm / L1.4 GHz: Appleton et al., 2004, Ibar et al., 2008) would avoid the extinction problem however, in this case there is the potential for strong contamination of the mid-IR luminosities by star formation in the host galaxies.

An alternative to the SED-based approach is to use a single radio power cut above which an AGN can be considered as radio-loud. This is the approach adopted for this review: I define a radio AGN to have a monochromatic radio luminosity measured at 1.4 GHz of L1.4 GHz > 10 24 W Hz 𢄡 . Although this cut may seem arbitrary given that the general population of AGN shows a continuous range of radio power, it does have some physical basis in the sense that the populations of objects in the local universe that are generally considered radio-quiet in terms of AGN properties, such as starburst galaxies and classical, UV-selected Seyfert galaxies, have steeply declining luminosity functions above L1.4 GHz = 10 23 W Hz 𢄡 such objects are rare above L1.4 GHz = 10 24 W Hz 𢄡 (Meurs & Wilson, 1984, Condon, 1989, Sadler et al., 2002).

In what follows I will also sometimes draw a distinction between luminous, quasar-like AGN and their less luminous counterparts. Originally quasars were defined by their star-like appearance in optical images. However, with the detection and characterisation of the underlying host galaxies, it has become clear that quasars represent the higher luminosity end of a general AGN population that shows a continuous range of luminosities. Here I define quasars to have bolometric luminosity Lbol > 10 38 W. Roughly, this corresponds to 2 – 10 keV X-ray luminosities LX−ray > 10 37 W, optical absolute magnitudes MB < � and optical [OIII] � emission line luminosities L[OIII] > 10 35 W, depending on the precise SEDs and bolometric correction factors assumed. However, I will not distinguish between broad-line radio galaxies (BLRG: see below for description) and radio-loud quasars, because there is a significant overlap in the properties of these two groups, and some BLRG show evidence for relatively high levels of dust extinction (e.g. Osterbrock, Koski & Philips, 1976) consistent with them being partially obscured quasars I will refer to these objects collectively as BLRG/Q sources, but they also sometimes labelled broad-line objects (BLO).

In many sciences, including astronomy, the classification of a phenomenon often precedes its interpretation in terms of the underlying mechanisms. Faced with a diverse phenomenon, the process of classification is an attempt to sort things out, a first step in elucidating underlying patterns of behaviour. While the classifications themselves remain relatively fixed, their interpretation in terms of the underlying physical processes, structures, geometries etc. may not be unique and may change with time. For these reasons it is important to keep the classifications separate from the interpretations. In this sub-section I summarise the radio and optical classifications of radio AGN possible interpretations are then discussed in the following sub-sections. For ease of reference, the labels used to classify radio AGN are summarised in Table 2.

An interesting feature of classification in astronomy is that a particular object may me classified in different ways at different wavelengths. For radio AGN some of the greatest insights have been obtained when attempting to understand the relationships between their classifications at the different wavelengths. The main radio and optical classifications of radio AGN are illustrated in Figure 1.

The radio classification of radio AGN followed the development of large radio interferometers in the 1970s that were able to map the extended radio-emitting structures in some detail. It was found that extragalactic radio sources are generally double in character, with synchrotron-emitting lobes that sometimes contain compact higher surface brightness concentrations called hot spots, connected to radio cores sources in the nuclei of the galaxies by jets (e.g. Miley, 1980). The radio structures often extend on scales that are larger than the host galaxies (� kpc to 1 Mpc).

One of the main radio morphological classifications is that of Fanaroff & Riley (1974), who divided objects according to whether the distance between the highest surface brightness extended radio features on either side of the nucleus was more or less than 50% of the total diameter of the radio source, corresponding to the FRII and FRI classifications respectively (see Figure 1 for examples). Fanaroff & Riley showed that the FRII sources have higher radio powers on average than FRI sources, with the division between the two types occurring at a radio power of L178 MHz ∼ 10 26 W Hz 𢄡 (or L1.4 GHz ∼ 3 × 10 25 W Hz 𢄡 ). However, it was already recognised by Fanaroff & Riley that a minority of sources (�% in their sample) could not be unambiguously assigned to either the FRI or the FRII class. Certainly, with the higher resolution radio maps available since the early 1980s it has become clear that some sources have an apparently hybrid FRI/FRII character (see the example of Hercules A in Figure 10). Also, the radio power division is not sharp: some sources with an FRII morphology have radio powers well below L178 MHz = 10 26 W Hz 𢄡 , and vice versa. Some ambiguity in this classification may be caused by the fact that, with the availability of higher resolution radio maps, the original quantitative criterion of Fanaroff & Riley has often been replaced by a subjective assessment of whether the morphology of a source appears �ge brightened” or �ge darkened”.

Although the FRI/FRII classification is the most common, other radio morphological classifications are also used 2 . For example, some radio sources resemble double-lobed, edge-brightened FRII sources, with similar steep-spectrum radio synchrotron emission, but have much smaller linear diameters. The class of compact steep-spectrum (CSS) sources has been defined (see Fanti et al., 1990, O'Dea, 1998) to include steep spectrum radio sources with total diameters D ≤ 20 kpc. Further diversity in the radio classification is added by objects that, rather than showing the classical double-lobed morphology, exhibit a core-jet or core-halo structure (e.g. Antonucci & Ulvestad, 1985a). Finally, some extended double-lobed radio sources have been classified as �t doubles” (FD) based on the fact that their lobes are unusually extended in the direction perpendicular to the axis defined by the radio jets (Owen & Laing, 1989). Such sources are also sometimes labelled “relaxed doubles” (RD) on the basis that their radio lobes lack bright hotspots and jets (Leahy, 1993). However, both the fat and the relaxed classifications lack clear, quantitative definitions.

In addition to the morphological classifications, radio sources are also classified on the basis of their long wavelength SED shapes, often quantified in terms of the spectral index (α) of a power law fitted to their radio spectra when expressed in frequency units (Fν∝ ν −α ). Flat- and steep-spectrum radio sources are defined to have spectral indices smaller or greater than a particular limiting value (typically αlim ∼ 0.5, but definitions vary). It is notable that, whereas steep-spectrum radio sources are commonly associated with FRI, FRII or CSS radio morphologies, flat-spectrum sources are generally core dominated, appearing as core-jet or core-halo sources in high resolution very long baseline interferometry (VLBI) observations (e.g. Antonucci & Ulvestad, 1985a). However, some radio sources have more complex radio spectra that cannot be approximated using a simple power-law fit. In particular, the GHz-peaked sources (GPS) have radio spectra that peak at GHz frequencies and decline at higher and lower frequencies (Bolton, Gardner & Mackey, 1963) such spectra are often attributed to synchrotron self-absorption effects (Kellermann, 1966). Morphologically, GPS sources often show double-lobed radio structures that are even more compact (D ≤ 1 kpc) than those of CSS sources 3

Most optical classifications of radio AGN rest on their emission line properties. In the 1970s Osterbrock and colleagues identified a host of narrow emission lines in several radio galaxies that they used to measure key properties of the narrow line region (NLR) such as reddening, density and temperature (e.g. Osterbrock & Miller, 1975). Paralleling the classification of Seyfert galaxies into Seyfert type 1 and Seyfert type 2, based on respectively the presence or absence of broad (FWHM > 2000 km s 𢄡 ) permitted lines in the spectra (Khachikian & Weedman, 1971), they also divided the radio galaxies into two types: broad-line radio galaxies (BLRG Osterbrock, Koski & Philips, 1976) and narrow-line radio galaxies (NLRG: Costero & Osterbrock, 1977), with BLRG showing spectral similarities to quasars, but having less luminous AGN. Collectively the NLRG and BLRG/Q spectroscopic classes can be labelled as strong-line radio galaxies (SLRG).

By the late 1970s, however, it became clear that the NLRG/BLRG/Q classification does not capture the full spectral diversity of the radio AGN population: in a study of 3CR radio galaxies Hine & Longair (1979) recognised a new class of objects (labelled 𠇌lass B”) that have “only the absorption spectra typical of giant elliptical galaxies or else very weak [OII] �”. Understanding these objects, and how they relate to the SLRG, has been a major focus of research on radio AGN in the last 20 years over this period improved spectroscopic data has allowed more quantitative classification schemes to be developed. Such objects are now identified in samples of radio AGN on the basis of their low [OIII] � equivalent widths (e.g. the weak-line radio galaxy [WLRG] classification of Tadhunter et al. 1998) or their low excitation/ionization emission line spectra (e.g. the low excitation galaxy [LEG] class of Buttiglione et al. 2010), or one or other of [OIII] equivalent width and excitation/ionization criteria (Laing et al., 1994, Jackson & Rawlings, 1997, Best & Heckman, 2012).

The different methods used to classify the Class B/WLRG/LEG objects have the potential to lead to ambiguity. Fortunately, most objects classified as WLRG on the basis of their low equivalent width [OIII] emission would be also be classified as LEG on the basis of the low excitation/ionization of their emission line spectra, as quantified using emission line ratios and diagnostic diagrams. Indeed, the [OIII] equivalent width is strongly correlated with excitation/ionization state in the radio AGN population (see Figure 1 of Best et al. 2012). However, the LEG and WLRG classifications are not exactly the same, and a minority of objects classified as LEG are not classified as WLRG. For example, in the study of Buttiglione et al. (2010) a few of the objects classified as LEGs on the basis of their low excitation/ionization emission line spectra have relatively strong, high equivalent width [OIII] emission lines (e.g. 3C84 4 , 3C153, 3C196.1, 3C349).

The data in Buttiglione et al. (2009, 2011) can be used to make a more quantitative comparison between the classification schemes. Taking the 99 3CR radio AGN with redshifts z < 0.3 in the Buttiglione et al. (2010, 2011) study with secure spectroscopic classifications, 2% are classified as star forming objects, 46% as LEGs 5 , 33% as high excitation galaxies (HEGs), and 18% as broad line objects (BLOs). In comparison, according to the WLRG/NLRG/BLRG/Q classification scheme of Tadhunter et al. (1998), 32% of the same objects would be classified as WLRG, 48% as NLRG and 19% as BLRG/Q objects. Overall, 65% of the 43 objects classified as LEGs by Buttiglione et al. would be classified as WLRG based on their [OIII] equivalent withs.

An advantage of the Buttiglione et al. (2010) scheme is that it is based on a clear, quantitative criterion that is motivated by the positions of the points on emission line diagnostic diagrams (see Figure 2). In this scheme LEGs and HEGs are classified according to their excitation indices (EI):

where [OIII] / Hβ represents the ratio of the flux of the [OIII] � line to that of the Hβ line, while [NII] / Hα, [SII] / Hα and [OI] / Hα represent respectively the ratios of the fluxes of [NII] �, [SII] λ�,6731 and [OI] � lines to that of Hα LEGs are defined to have EI < 0.95 6 , while HEGs are defined to have EI > 0.95. Importantly, there is evidence that the distribution of excitation class in the sample is bimodal, with one peak representing the HEGs/BLO and the other the LEGs. However, a disadvantage of this method is that it requires accurate measurements of all the lines involved in the excitation index. For objects at the lower end of the [OIII] equivalent width (EW) distribution this information is difficult to obtain, since several (or all) of the emission lines are undetected or have flux measurements with large uncertainties. For this reason, 14 (12%) of the 113 3CR objects with redshifts z < 0.3 in the full sample of Buttiglione et al. (2010, 2011), lack spectroscopic classifications. The latter objects would be classified as WLRG in the Tadhunter et al. (1998) scheme, leading to 43% of the full sample of Buttiglione et al. (2009, 2010, 2011) being classified as WLRG.

Therefore a major advantage of the alternative technique of using an upper limiting [OIII] EW is that it captures the objects with low EW or undetected [OIII] that cannot be classified on the basis of excitation class. On the other hand, using a particular [OIII] equivalent width limit (EW[OIII] < 10 Å in the case of Tadhunter et al. 1998) is somewhat arbitrary because, while there is some evidence for bi-modality in the [OIII] emission line luminosities and equivalent withs for the high radio power end of the 2Jy sample, at low radio powers the distribution [OIII] equivalent widths appears continuous. Also, the equivalent width of [OIII] depends not only on the [OIII] emission line flux but also on the flux of the underlying stellar continuum. The latter may vary from object to object (for a given [OIII] luminosity), and be subject to aperture and reddening effects.

An important general point about the LEG and WLRG classes is that, while their emission line ratios are similar to the those of the LINER class of radio-quiet AGN (Heckman, 1980), there are differences in terms of emission line luminosity. As pointed out by Capetti et al. (2013), many of the LEGs in the 3CR sample have emission line luminosities that are comparable with those of Seyfert galaxies, but much higher than those of LINERs the same is true of many of the 2Jy radio AGN classified as WLRG by Tadhunter et al. (1998).

One further classification method is based on optical variability. A subset of objects — labelled blazars — exhibit extreme variability at optical wavelengths. This class includes BL Lac objects and optically violently variable (OVV) quasars, with the two groups separated on the basis that BL Lac objects show emission lines that are weak relative to the underlying featureless continuum (e.g. EW < 5Å : Urry & Padovani 1995), whereas OVV quasars exhibit broad emission lines as well as a strong non-stellar continuum a relatively high degree of optical polarisation is also common in both of the latter groups. Moreover, most such objects are flat-spectrum radio sources with core/halo or core/jet morphologies.

Finally, it is important to consider the degree to which the various optical and radio classifications are correlated. Table 3 shows the frequency of the different optical spectroscopic classifications for FRI, FRII, CSS/GPS, and hybrid FRI/FRII or uncertain radio morphological classifications in both the full z < 0.3 3CR sample of Buttiglione et al. (2009, 2011: 113 objects) and the full z < 0.7 2Jy sample of Tadhunter et al. (1998: 66 objects) 7 . The most striking feature of this comparison is that all the objects unambiguously classified as FRI radio sources are WLRG none are classified as BLRG/Q 8 . In addition, all SLRG have FRII or CSS/GPS radio morphologies, apart from a small minority that have hybrid FRI/FRII or ambiguous morphologies. It is also clear that, while the majority of FRII sources are SLRG, a significant minority (∼ 24%) are WLRG.

The orientation-based unified schemes were developed following the recognition in the late 1970s and early 1980s that the continuum emission from AGN is highly anisotropic. In the case of radio AGN, there are two types of anisotropy to consider. First, the beaming of the non-thermal synchrotron emission due to the bulk relativistic motions of the jets in the core regions of the sources. Direct evidence for the bulk relativistic motions and beaming is provided by the detection of apparently superluminal velocities in VLBI observations of the radio cores of blazar-like objects (Cohen et al., 1977), and the detection of a polarisation asymmetries in the radio lobes of radio-loud quasars that correlate with sidedness of the inner radio jets (Garrington et al., 1988, Laing, 1988). Second, the blocking effect of dust and gas in the circum-nuclear dust structures that are often characterised as tori. The presence of this form of anisotropy is supported by detection of significant optical polarisation that is aligned perpendicular to the radio axes in nearby radio galaxies (Antonucci, 1982, Antonucci, 1984), the detection of broad emission lines in the polarised intensity spectra of nearby Seyfert galaxies and NLRG (Antonucci, 1984, Antonucci & Miller, 1985b, Ogle et al., 1997, Cohen et al., 1999), and the detection of “ionisation cones” in narrow-band emission line images of Seyfert galaxies (e.g. Pogge, 1988, Tadhunter & Tsvetanov, 1989b) and radio galaxies (Jackson, Tadhunter & Sparks, 1998). The effect of such blocking is to heavily attenuate both the continuum and broad line region (BLR) emission of the AGN at UV and optical wavelengths in objects for which the radio jet/torus axis is pointing at a large angle to the line of sight.

The anisotropic AGN continuum emission can help to explain some of the observed diversity in the properties of the radio AGN. Perhaps the most ambitious anisotropy-based unified scheme for radio AGN involves both forms of anisotropy, and attempts to explain the relationship between steep-spectrum radio-loud quasars (SSRQ), flat-spectrum radio-loud quasars (FSRQ), and radio galaxies, (Barthel, 1989). This scheme is shown schematically in Figure 3(a). The basic idea is that, for a given radio power, the different classes of radio sources are all drawn from the same parent population with similar central AGN properties, and the differences between their radio and optical properties are explained in terms of anisotropy and orientation. In this scheme, the relative proportions of FSRQ, SSRQ/BLRG and NLRG in a particular radio AGN sample is set by the opening half-angles of both the torus (θtor) and the beaming cones of the relativistic jets (θbeam ∼ Γ 𢄡 , where Γ is the Lorentz factor of the bulk relativistic motions in the jet). Recent analysis of the radio and optical properties of the HEG and BLO with FRII morphologies in the z < 0.3 3CR sample is consistent with θtor = 50 ± 5 degrees and Γ ∼ 3 — 5 (Baldi et al., 2013).

One way to test such schemes is to compare the statistical properties of the different classes of radio AGN in flux-limited samples selected in well-defined redshift/radio luminosity ranges. For example, Barthel (1989) compared the distributions of the measured (i.e. projected) linear sizes of the radio sources of the radio galaxies and the quasars in his sample of intermediate- to high-redshift (0.5 < z < 1) radio galaxies, and found that the radio galaxies have larger diameter radio sources on average than the quasars, as expected in the case that we have a foreshortened view of the radio sources of the quasars. However, attempts to repeat the linear diameter test on other samples of radio sources, albeit samples selected in different redshift ranges and/or with different radio flux limits, have failed to find significant differences between the linear size distributions of radio galaxies and quasars (Singal et al., 1993, Kapahi et al., 1995). Note that, this lack of difference between the linear size distributions is not necessarily evidence against the orientation-based unified schemes, but could rather reflect selectional biases in the samples, incompleteness in the optical classifications of the targets, mixing up different populations of radio sources, or the effects of the radio power evolution of the sources coupled with a luminosity dependence in the opening angle of the torus (e.g. Gopal-Krishna, Kulkarni & Wiita, 1996).

Another statistical test involves comparing the radio powers of the extended radio lobes of powerful radio galaxies and quasars. In this case, since the radio lobes are not thought to be undergoing bulk relativistic motions that would cause anisotropy in their radio emission, we would expect the radio galaxies and quasars to have similar distributions of lobe power. In general, this is borne out by the observations (e.g. Urry & Padovani, 1995). Clearly, a major advantage of radio AGN for tests of the unified schemes is the isotropy of the emission from their radio lobes: samples selected on the basis of their low frequency radio emission — assumed to be lobe-dominated — can be considered orientation-independent.

The recent dramatic improvement of the optical and infrared data for samples of nearby radio sources has also allowed tests of orientation-based unification based on the optical and mid-IR emission lines and mid-IR continuum: if BLRG/Q and NLRG objects are drawn from the same parent population, on average the two groups should show similar luminosities in their narrow emission lines and mid-IR continua, assuming that the latter are emitted isotropically (i.e. they do not suffer attenuation by the circum-nuclear dust).

As one of the brightest optical emission lines that represents the high ionization conditions typical of the NLR, the [OIII]� forbidden line was the first to be used in this way. Early results suggested that the quasars are up to an order of magnitude more luminous in [OIII] than NLRG (Jackson & Browne, 1990). In contrast, the lower ionisation [OII] � line — which is likely to be emitted on larger scales than the [OIII] — showed no significant difference between the two groups (Hes, Barthel & Fosbury, 1993). Therefore, rather than providing evidence against the unified schemes, it was proposed that all, or part, of the [OIII]-emitting NLR is emitted on a relatively small scale and is subject to attenuation by the circum-nuclear dust in NLRG. This is supported by the observation of variability in the [OIII] emission lines of the BLRG 3C390.3, which provides evidence that much of the [OIII] emission in that object is emitted on a scale r < 10 pc (Clavel & Wamsteker, 1987, Zheng et al., 1995). However, some of the early studies that compared the [OIII] emission line luminosities of radio galaxies and quasars failed to distinguish between WLRG and SLRG. If WLRG represent a separate class of radio AGN with intrinsically lower luminosity AGN (see discussion in section 2.4 below), this could lead to the apparent differences between the [OIII] luminosities of narrow- and broad-line objects being exaggerated (Laing et al., 1994).

Deep mid-IR spectroscopy observations taken with the IRS spectrograph on the Spitzer satellite have deepened our understanding of the degree of attenuation suffered by the [OIII] narrow line and the mid-IR continuum emission (Haas et al., 2005, Dicken et al., 2014). Figure 4 shows example Spitzer/IRS spectra for radio AGN from the 2Jy sample and illustrates the diversity in the mid-IR spectra of radio AGN. Such observations have enabled accurate measurements to be made of the mid-IR [NeIII] 𻬕.6µm and [OIV] 𻬥.9µm fine structure lines, which are much less likely to suffer dust extinction than the optical forbidden lines. Figure 5 shows the [OIII] emission line and 24µm continuum luminosities plotted against the [OIV] 𻬥.9µm luminosity for complete sub-samples of 3CRR and 2Jy SLRG, with NLRG and BLRG/Q indicated by different symbols. It is clear that, on average, the NLRG fall below the correlations defined by the BLRG/Q. Under the reasonable assumption that the [OIV] emission is not significantly attenuated, this suggests that both the [OIII] and 24µm emission suffer a similar mild attenuation by factor of 𢏂 — 3 (Dicken et al., 2014) 9 . This is similar to the level of [OIII] extinction deduced by Baum et al. (2010) and LaMassa et al. (2010) for samples of nearby Seyfert 2 galaxies using Spitzer/IRS data.

In fact, it is remarkable that the degree of attenuation is apparently so similar for the [OIII] and 24µm emission in the radio AGN, given that the equivalent dust extinctions expressed in magnitudes in the V-band, are Av ∼ 1 magnitudes and Av ∼ 20 magnitudes for the [OIII] emission line and 24µm continuum respectively. Therefore, if the results are interpreted solely in terms of extinction, the [OIII] and 24µm emission must be extinguished by different dust structures, or by different parts of the same dust structure. For example, the 24µm continuum could be extinguished by the dust in the compact inner parts of the torus, and the [OIII] by larger-scale dust, perhaps in the outer parts of the torus, or a kpc-scale dust lane in the host galaxy. However, it is not possible to entirely rule out the possibility that, rather than being due to extinction, the difference between the 24µm luminosities of the NLRG and BLRG/Q objects at fixed [OIV] luminosity is due to contamination by beamed synchrotron emission from the radio cores, given that the degree of synchrotron contamination would be expected to be larger in the case of the BLRG/Q objects, because their jets are pointing closer to the line of sight.

An alternative technique involves using the high spatial resolution of the HST to investigate the rate of detection of compact core sources at both optical and infrared wavelengths in samples of nearby radio sources. The bright AGN nuclei are expected to be heavily attenuated at optical wavelengths by dust in the torus for objects observed as NLRG. However, due to the decrease in dust extinction towards longer wavelengths, the torus will become more transparent in the near-IR, potentially allowing the direct detection of the compact nuclei. Early attempts to detect the compact near-IR nuclei in radio galaxies using ground-based observations produced some promising results for a few individual radio galaxies (Djorgovski et al., 1991, Simpson, Ward & Wilson, 1995), but were hampered by the difficulty of separating the AGN nuclei from the starlight, given the relatively modest spatial resolution. Therefore the HST, with its order of magnitude better spatial resolution, has revolutionised this area (see Baldi et al., 2010, Ramírez et al., 2014a). Figure 6 shows the detection rate of AGN nuclei in low redshift (z < 0.11) NLRG as a function of wavelength from the optical to the mid-IR 10 . As expected, the rate of detection of AGN nuclei in the NLRG sources is low at optical wavelengths (< 30%), although interestingly not zero. The detection rate for the NLRG then rises with increasing wavelength through the near-IR to the mid-IR bands, reaching 80% at 2µm and 95% at 8µm (Ramírez et al., 2014a). This behaviour is entirely as expected on the basis of the standard orientation-based unified schemes. The infra-red fluxes and SED shapes deduced for the nuclei of NLRG are consistent with degrees of dust extinction in the range 10 < Av < 200 magnitudes (Tadhunter et al., 1999, Ramírez et al., 2014a).

At this stage it is important to add a caveat about the nature of the unresolved near-IR cores detected in NLRG. The near-IR light of type 1 AGN is dominated by the thermal emission of dust at close to the sublimation temperature in the inner parts of the obscuring torus, rather than accretion disk emission. Therefore it is possible that the unresolved nuclei detected in NLRG at near-IR wavelengths represent directly transmitted dust emission from the inner parts of the torus. However, alternative possibilities include non-thermal emission from the inner synchrotron jets and light scattered by dust in the near-nuclear regions (e.g. the far wall of the torus). Unfortunately, the recent detection of high degrees of linear polarisation in the unresolved near- and mid-IR nuclei of some NLRG (Tadhunter et al., 2000, Ramírez et al., 2009, Ramírez et al., 2014b, Lopez-Rodriguez et al., 2014) does not entirely resolve this issue, because the synchrotron, scattered AGN and transmitted AGN mechanisms could all produce significant IR polarisation (in the latter case via the dichroic extinction effect of aligned dust grains in the torus: Ramirez et al. 2009, 2014b). Note that if the near-IR emission of NLRG were dominated by synchrotron or scattered AGN emission, the level of extinction to the AGN would be higher than estimated based on the assumption that all the near-IR core emission is directly transmitted AGN light.

It is also possible to search for direct AGN continuum emission from NLRG at harder X-ray energies (> 2 keV), where the level of absorption by gas in the torus is reduced compared to that at EUV and soft-X-ray wavelengths. Again this is an area that has benefitted from improvements in technology in the last 20 years, with the launch of the Chandra and XMM satellites providing improved spatial resolution and sensitivity compared with previous X-ray satellites. One complication of the X-ray studies is that emission from the inner non-thermal jets — which are likely to suffer less attenuation than that of the AGN themselves — may contribute to the X-ray core fluxes. However, by explicitly accounting for the non-thermal component in the modelling of the X-ray spectra, it has proved possible to detect the attenuated X-ray emission from AGN in several NLRG: absorbed X-ray AGN components have been detected in 81% of the 21 NLRG 11 at z < 0.3 in the sample of 3CRR objects studied by Hardcastle et al. (2009), and in 84% of 19 NLRG in the 0.05 < z < 0.7 2Jy sample studied by Mingo et al. (2014). Thus the X-ray results are consistent with the near- to mid-IR results in the sense that they show the presence of heavily attenuated AGN is a high proportion of NLRG. However, for objects with both X-ray and near-IR HST observations, the attenuating HI columns are generally higher than predicted from the levels of near-IR dust extinction assuming the standard Galactic dust-to-gas-ratio (Ramírez et al., 2014a). This apparent inconsistency can be explained if AGN have a non-standard (order-of-magnitude higher) gas-to-dust ratio due to the AGN radiation field destroying dust grains in the inner part of the torus (Maiolino et al., 2001).

Despite the success of the X-ray and infrared imaging observations in detecting intrinsically luminous, obscured AGN in several NLRG, such observations do not allow the detailed spectra of the obscured AGN to be determined — an important step in demonstrating that the nuclei have quasar-like properties. Therefore, near-IR spectroscopy and optical spectropolarimetry observations — which have the potential to detect broad permitted lines characteristic of quasars in transmitted or scattered light — provide an important complement to the X-ray and infrared imaging observations.

Attempts have been made to detect the directly transmitted broad Paα emission using K-band spectroscopy of NLRG. However, despite claims of the detections in a few nearby NLRG (Hill et al., 1996), the observations have a low S/N, and this technique suffers from the fact that the infrared nuclei tend to be faint relative to the starlight of the cores of the host galaxies in most objects, and the broad Paα in typical type 1 AGN spectra has a low equivalent width compared with the optical Balmer lines (see discussion in Ramírez et al., 2009). Together, these factors make the detection of broad Paα in NLRG challenging, even with 8m-class telescopes.

Far more successful have been optical spectropolarimetry observations that use the alternative technique of detecting the broad Balmer lines in scattered light: there are now convincing detections of the scattered broad Hα lines in five NLRG in the local universe (Antonucci, 1984, Ogle et al., 1997, Cohen et al., 1999). Such observations provide the most direct evidence to support the orientation-based unified schemes for SLRG, because they demonstrate that individual NLRG have nuclei with the spectral characteristics of quasars.

Although scattered broad lines have so far been detected in a only handful of the dozens of NLRG in the 3CR and 2Jy samples at low redshifts (z < 0.2), the relatively low detection rate cannot be taken as strong evidence against the unified schemes. This is because several factors can confound the detection scattered broad lines, including the strong dilution of the polarised light by the starlight of the host galaxies, geometrical dilution of the polarization, a lack of scattering dust in the NLR, and illuminating AGN that are at the lower end of the intrinsic luminosity range for a given radio power. Moreover, many of the 3CR and 2Jy NLRG still lack deep spectropolarimetry observations. Unfortunately, it is unlikely that there will be rapid progress in this area in the near future, since making the requisite spectropolarimetry observations of a large sample of NLRG would be prohibitively expensive in observing time with the current generation of 8m telescopes.

Taken together, the statistical results from comparisons between the optical, mid-IR and radio properties of BLRG/Q and NLRG, the detection of highly attenuated AGN nuclei at X-ray, near-IR and mid-IR wavelengths in a high proportion of nearby NLRG, and the detection of polarised broad lines in scattered light in some NLRG, provide compelling evidence that the orientation-based unified schemes for SLRG work to first order: the data are consistent with the idea that all NLRG contain BLRG/Q nuclei that are obscured along our direct line of sight by circum-nuclear dust. The unification debate for SLRG now centres on the geometry of the central obscuring region (e.g. smooth torus, clumpy torus, warped disk), and also on whether the properties of the obscuring structures change with luminosity and redshift (e.g. Lawrence, 1991, Lawrence & Elvis, 2010, Elvis, 2012).

While the orientation-based unified schemes are successful in explaining the relationship between NLRG and BLRG/Q objects amongst the SLRG, they cannot readily explain the relationship between WLRG and SLRG, or between FRI and FRII sources. Although it has been proposed that the WLRG objects might be SLRG radio galaxies in which the NLR is unusually heavily obscured by circumnuclear dust (e.g. Cao & Rawlings, 2004), the Spitzer results presented in Figure 7 demonstrate that this cannot be the case: a heavily obscured SLRG nucleus would be expected to radiate strongly in the mid-IR continuum and emission lines, whereas the WLRG have weak mid-IR continuum and [OIV] line emission (Hardcastle et al., 2009, Dicken et al., 2014) many of the WLRG are also weak at far-IR wavelengths (Dicken et al., 2009) 12 , despite the fact that their continua are often dominated by non-thermal emission (Dicken et al., 2008, Leipski et al., 2009, van der Wolk et al., 2010). Similarly, the relationship between FRI and FRII sources cannot be explained in terms of orientation, because it is not possible to “hide” the strong radio lobe and hot spot emission typical of FRII sources in the FRI sources.

Noting that FRI sources are almost invariably associated with WLRG optical spectra, and that BL Lac objects also show low equivalent emission lines, an alternative orientation-based unified scheme (see Figure 3(b)) has been proposed to explain the link between these two classes of objects (Urry, Padovani & Stickel, 1991, Urry & Padovani, 1995). In this case, the required anisotropy in the optical continuum is produced entirely through the beaming effect caused by the bulk relativistic motions of the inner jets the presence or otherwise of a central obscuring torus in such objects is controversial. The beaming that is such a key part of this scheme explains why BL Lac objects — viewed with the line of sight close to the direction of the radio axis — have core-dominated radio morphologies and relatively flat radio spectra beamed jet emission also explains the relatively strong, point-like, and highly polarised optical continuum sources that are observed at optical wavelengths in such objects.

Interestingly, HST observations also reveal point-like nuclei at optical wavelengths in a high proportion of WLRG/FRI sources (see Figure 6 Chiaberge, Capetti & Celotti, 1999, Verdoes Kleijn et al., 1999, Capetti et al., 2002). The strong correlations between the luminosities of the optical nuclei and those of the radio cores (Chiaberge, Capetti & Celotti, 1999), the optical to X-ray SED spectral indices (Balmaverde, Capetti & Grandi, 2006), and the detection of high degrees of optical polarization (Capetti et al., 2007) are all consistent with the idea that the optical nuclei in the WLRG/FRI objects represent synchrotron emission from the inner jets. At first sight the detection of the optical nuclei in FRI sources might seem to go against the orientation-based unified schemes, since the WLRG/FRI sources are considered to be the unbeamed counterparts of the BL Lac objects. However, even in the case that the radiation from the inner jets is strongly beamed, it is expected that some radiation will be detected from the jets in objects whose jets are pointing at a large angle to the line of sight. Moreover, the inner synchrotron jets could extend beyond any central obscuring tori that would otherwise block the emission of their nuclei at X-ray and optical wavelengths. Finally, as already noted, it is not clear whether the obscuring tori are present in the nuclei of typical WLRG/FRI objects.

The unification of BL Lac objects with FRI radio galaxies in the context of the orientation-based scheme shown in Figure 3(b) can be tested in a similar way to unified scheme for SLRG: by statistical comparison of properties that are considered to be isotropic. In this case the results are mixed: while the extended radio lobe luminosities of BL Lac objects and FRI radio galaxies appear to be similar (Urry & Padovani, 1995), on average the [OIII] emission line luminosities of the BL Lac objects are significantly higher when the comparison is made for low redshift samples (z ≤ 0.2: Urry & Padovani, 1995, Wills et al., 2004) 13 . These results can be reconciled if there is sufficient circum-nuclear dust to significantly obscure the [OIII] emission in the WLRG/FRI objects. Certainly some WLRG/FRI objects do have kpc-scale dust lanes (de Koff et al., 2000). On the other hand, as we have already seen, most WLRG/FRI objects lack evidence for central obscuring tori.

A further relevant result is that broad permitted lines characteristic of type 1 AGN have been detected in some BL Lac objects — most notably Bl Lac itself (Corbett et al., 1996) and PKS0521-36 (Ulrich, 1981) — whereas, do date, there have been no convincing detections of broad permitted lines in any of the FRI objects in the 3CR and 2Jy samples 14 . Even allowing for some obscuring dust in the circum-nuclear regions, by analogy with the SLRG unification scheme, we might expect a significant fraction of the WLRG/FRI objects — those observed at intermediate angles (the equivalent of the steep spectrum radio-loud quasars) — to show broad lines if these objects are truly the unbeamed counterparts of the BL Lac objects, and most BL Lac objects have a BLR in their nucleus. However, the result that FRI sources lack broad lines should be treated with some caution, because the broad lines might be difficult to detect against the strong stellar emission of the cores of the host elliptical galaxies, especially if the AGN luminosity (including the broad-line emission) scales with the extended radio luminosity (Rawlings & Saunders, 1991), such that the relatively low radio luminosity FRI sources have correspondingly low luminosity BLR emission.

Taking the results on FRI/BL Lac unification together, it seems unlikely that all BL Lac objects would appear as FRI sources if they were observed with their radio jets pointing at a large angle to the line of sight however, it is entirely plausible that a subset of such objects (i.e. those with lower luminosity [OIII] emission) can be unified with FRI sources in this way.

This leaves the WLRG/FRII sources as true misfits in terms of the orientation-based unified schemes: they cannot be readily unified with WLRG/FRI sources because of the luminosities and morphologies of their extended radio emission, yet they cannot be unified with SLRG/FRII sources because of their low [OIII] emission line luminosities, lack of broad line emission, and low mid-IR continuum and emission line luminosities. Note that the latter mid-IR properties rule out the idea that the WLRG/FRII sources contain luminous AGN that are unusually heavily obscured, just as they do in the case of the WLRG/FRI sources.

Overall, it is clear that the orientation-based unified schemes are successful at explaining the relationship between NLRG and BLRG/Q objects in the case of SLRG, and less certainly between FRI radio galaxies and BL Lac objects in the case of WLRG. However, they cannot by themselves explain the full diversity of properties of the radio AGN population. In the following sections I discuss other factors that are important in determining the observed properties of radio AGN many of these have come to light in the last 20 years.

Aside from anisotropy/orientation effects, it is natural to consider the possibility that the varied manifestations of the radio AGN phenomenon are related to the mechanisms that produce the AGN radiation and jets in the nuclear regions — the so-called �ntral engines” of the AGN — in particular, changes in rate or mode of accretion of material by the central supermassive black holes.

In this context, it is significant that there is a strong correlation between the radio morphological properties (e.g. Fanaroff & Riley types I and II) and the optical spectroscopic properties (e.g. WLRG/SLRG) of radio AGN (see Table 3). As we have seen, FRI radio sources are almost invariably associated with WLRG spectra, and SLRG are almost invariably associated with FRII or CSS/GPS radio sources only the WLRG/FRII objects break the trend. This strongly suggests that radio morphological and optical spectroscopic properties are linked primarily through the natures of the central engines, rather than through external environmental factors. Note, however, that environmental factors may nonetheless affect the radio properties at some level. For example, entrainment of the relatively dense, hot ISM found at the centres of giant elliptical galaxies and clusters of galaxies may be required to explain the detailed properties of FRI jets (e.g. Laing & Bridle, 2014) and the confinement effect of a hot, dense ISM may act to boost the observed radio luminosity for a given intrinsic mechanical jet power (Barthel & Arnaud 1996, but see Ramos Almeida et al. 2013).

Further important clues are provided by examining the links between the levels of AGN activity, the host galaxy absolute magnitudes, and the radio morphological or optical spectroscopic classifications. Such links were first discussed in detail by Owen (1993) and Ledlow & Owen (1996) who plotted monochromatic radio power against host galaxy absolute magnitude and identified FRI and FRII sources with separate symbols (see Figure 8). Essentially, their plot demonstrates that the radio power that defines the break between FRI and FRII sources increases with increasing host galaxy luminosity: FRI sources hosted by luminous (and more massive) elliptical galaxies can have higher radio powers than those hosted by less luminous galaxies, without breaching the FRI/FRII radio power limit.

A possible explanation for the trend shown in Figure 8 in terms of the accretion processes close to the supermassive black holes was suggested by Ghisellini & Celotti (2001) who noted that the radio power could be taken as a proxy for the overall level of AGN activity, while the host galaxy absolute magnitude is related to the mass of the central supermassive black hole (e.g. Kormendy & Ho, 2013). Therefore, Figure 8 could be recast as a plot of AGN power against the black hole mass, with the dividing line between FRI and FRII sources representing a fixed ratio of the AGN power to the Eddington luminosity of the black hole. This is the first suggestion that the division between FRI and FRII sources is due to an �ington switch”: a change in the nature of the accretion flow from a radiatively inefficient accretion flow (RIAF) that is geometrically thick but optically thin (FRI objects) to a standard geometrically thin but optically thick accretion disk (FRII objects) at a critical ratio of the jet power (Qjet) to the Eddington luminosity of the central black hole (Ledd). The critical ratio derived by Ghisellini & Celotti (2001) was in the range 10 𢄣 < Qjet / Ledd < 10 𢄢 , although a more recent study by Wold, Lacy & Armus (2007) that used host galaxy velocity dispersion rather than absolute magnitude as a proxy for black hole mass found Qjet / Ledd ∼ 5 × 10 𢄤 .

Figure 8 can also be redrawn as a plot of emission line luminosity against host galaxy absolute magnitude, this time identifying objects by their optical spectroscopic classifications rather than their radio morphologies (Buttiglione et al., 2010). The results for the 3CR sample are shown in Figure 9. Just as there is a clear dividing line between FRI and FRII radio galaxies in Figure 8, there is a clear demarcation between HEG and LEG objects in Figure 9, with some evidence that the emission line luminosity that represents the boundary between the two classifications increases with host galaxy luminosity. Again, since the [OIII] luminosity is a good indicator of overall AGN power (e.g. Heckman et al., 2005), this trend can be interpreted in terms of an Eddington switch: LEGs are associated with radiatively inefficient accretion, and the HEG/BLOs with a standard thin accretion disk, with the transition between the two types occurring at Lion / Ledd ∼ 10 𢄣 , where Lion is the ionising luminosity of the AGN. This is also consistent with the recent results of both Best & Heckman (2012) for SDSS-selected radio AGN and Mingo et al. (2014) for the 2Jy and 3CR samples, who find a clear division between the Eddington ratios of low- and high-excitation objects at (Lbol + Qjet) / Ledd ∼ 10 𢄢 , where Lbol is the bolometric (radiative) luminosity of the AGN.

Interestingly, an Eddington ratio of � 𢄢 also represents the division between different accretion states in X-ray binary systems: from the low/hard to the high/soft states (Macarrone, 2003). Therefore, in terms of interpreting diagrams like Figs 8 and 9, analogies have often been drawn with X-ray binary systems (e.g. Falcke, Körding & Markoff, 2004, Körding, Jester & Fender, 2006), and it has been proposed that X-ray binaries and AGN together form part of a main sequence in black hole activity (Merloni, Heinz & di Matteo, 2003). A potential problem with this picture is that the high/soft state in X-ray binary systems is not commonly associated with the formation of non-thermal jets. However, Nipoti, Blundell & Binney (2005) have linked radio AGN with a transitory jet flaring stage of activity as sources pass between the low/hard and high/soft states (see also Körding, Jester & Fender, 2006), and radio-quiet AGN with the more usual non-flaring high/soft state. In this case, the transitory nature of jet flaring stage could be consistent with the fact that radio-loud quasars represent < 10% of the full quasar population: perhaps objects cycle through radio-loud and radio-quiet phases as part a particular quasar triggering event (but see Bessiere et al., 2012).

An important debate concerns whether the different radio AGN classifications, in particular the division between WLRG/LEG and SLRG/HEG, are related to the mode rather than the rate of accretion. On the basis that the powers of the jets of LEGs in the 3CR sample are consistent with fuelling via Bondi accretion of the hot ISM from the X-ray haloes of the objects, whereas HEGs are too powerful to be explained in this way, Hardcastle, Evans & Croston (2007) proposed that LEGs are associated with hot mode accretion, whereas HEGs are associated with cold mode accretion which has the capability to fuel the AGN at the requisite higher accretion rates in this model the division between FRI and FRII sources is explained in terms of environmental effects. Buttiglione et al. (2010) also favoured the mode of accretion as being the dominant factor, arguing that any hot gas being accreted will not be able to cool to form a thin accretion disk and will therefore be associated with WLRG/LEG nuclear activity in radio AGN. On the other hand, cool gas in the process of being accreted will naturally form a standard accretion disk, leading to HEG-style nuclear activity.

Despite these arguments, it seems unlikely that all WLRG/LEG are fuelled by hot gas. Indeed, many have plentiful cold gas in their nuclear regions, as evidenced by nuclear dust lanes (de Koff et al., 2000), at least one shows evidence at mid-IR wavelengths for a compact, warm dust structure close to the nucleus (PKS0043-42: Ramos Almeida et al., 2011b), and the nuclear NLR of several WLRG/FRI sources — albeit of low luminosity compared to SLRG — appear compact and centred on nuclei in HST long-slit and imaging observations (Verdoes Kleijn et al., 1999, Capetti, Verdoes Keijn & Chiaberge, 2005), thus providing evidence that there is at least some warm gas on a sub-100 pc scale in these objects, even if the gas masses are relatively low. Therefore, it seems more plausible that accretion rate is the dominant factor, and that some WLRG are fuelled by cold mode accretion at a relatively low rate, even if the majority are fuelled by hot mode accretion (see also Best & Heckman, 2012).

Until recently, most studies of hot mode accretion assumed that the hot gas accretes as a spherical, unperturbed flow at the Bondi rate (Bondi, 1952). However, recent simulations that include realistic prescriptions for cooling, heating and turbulence, provide evidence that much of the hot gas in the nuclear regions on a scale of 𢏀.1 — 1 kpc can condense into filaments of cool gas that fall near-radially towards the nucleus, with cloud-cloud collisions helping to dissipate any angular momentum in the gas (e.g. Gaspari, Peng & Oh, 2013, Gaspari, Brighenti & Temi, 2015). As a result, the accretion rates associated with such chaotic cold accretion (CCA) can be orders of magnitude higher than the Bondi rate — sufficient to trigger a quasar-like AGN. The CCA simulations — which clearly blur the distinction between hot and cold mode accretion — also predict the presence of large masses of cold gas and the formation of clumpy torus-like structures close to the AGN. Since these predictions are somewhat at odds with what we observe in the majority of WLRG/FRI objects, the applicability of CCA models to such objects is currently uncertain. Furthermore, if the CCA mechanism is as effective as the models suggest, we would expect luminous, quasar-like AGN to be common in rich clusters of galaxies, where the densities of hot ISM are relatively high however, this appears inconsistent with the measured environments of most SLRG (see section 3.5).

On the basis of diagrams such as Figures 8 and Figure 9 it might seem plausible that the radio morphological and optical spectroscopic classifications of radio AGN are linked through the nature of the accretion onto the black holes via the Eddington switch i.e. FRI objects are always WLRG and are associated with low rates accretion and a radiatively inefficient accretion flow, whereas FRII objects are always SLRG and are associated with thin, radiatively efficient accretion disks. However, the WLRG/FRII sources do not fit into this picture. While the reason for this apparent discrepancy might be related to some aspect of the accretion flow onto the black hole that can generate powerful jets and an FRII radio morphology but not a SLRG nucleus, or perhaps to some particular combination of environmental and nuclear accretion factors, an alternative possibility is that it is due do the intermittency or even switching off the nuclear fuel supply. This possibility will be explored in the next subsection.

All types of broad-line AGN are known to vary by up to a factor of a few at optical wavelengths on timescales of weeks to years (Matthews & Sandage, 1963, Fitch et al., 1967). Radio AGN are no exception, and there are now detailed studies of bright, nearby BLRG that track the variability of both the broad lines and continuum, and use the time lags between the two types of emission to estimate the size of the broad-line region (e.g. Dietrich et al., 2012).

Apart from this “normal” AGN variability, which is likely to be due to relatively minor changes in the accretion flow that lead to a temporary increase or decrease in luminosity of the thermal accretion disk emission and hence the flux of ionising photons in the BLR, it is also possible that there are longer-time-scale variations in the fuelling of the AGN that lead to larger amplitude (i.e. factor of 3 or more) variations. There are two situations to consider here: high-amplitude intermittency in the fuel supply within a particular cycle of SLRG activity and the switch-off phase at the end of a SLRG cycle, in which a longer-term scarcity of fuel causes the central AGN to shut down completely or enter a lower activity state associated with a low Eddington ratio.

Evidence for intermittency is provided by the radio structures of some radio sources. For example, the radio structure of the WLRG Hercules A (3C348: see Figure 10), shows hybrid FRI/FRII characteristics, including a series of bubble structures that may represent successive phases of high jet activity in the lifecycle of the source. Other objects show evidence for more than one phase of activity in the form of 𠇍ouble-double” radio structures (Schoenmakers et al., 2000), or high surface brightness, compact inner double structures combined with much larger-scale and lower-surface-brightness outer structures (e.g. PKS1345+12: Stanghellini et al., 2005).

The observational consequences of high-amplitude changes in the fuel supply depend on the timescales of the changes, and on the AGN components being considered. Observations show that the BLR, torus, NLR, and large-scale radio-emitting lobes of radio AGN have typical radial scales of rBLR ∼ 0.01 − 1 light year, rtorus ∼ 0.1 − 100 pc, rnlr ∼ 0.001 − 3 kpc 15 , rlobe ∼ 0.05 − 1 Mpc, corresponding to light crossing times of τBLR ∼ 0.01 − 1 yr, τtorus ∼ 0.3 − 300 yr, τnlr ∼ 3 − 10,000 yr and τlobe ∼ 150,000 − 3 × 10 6 yr respectively. Thus, when an AGN finally switches off due to the exhaustion of its fuel supply, this change will lead to a substantial decrease in the BLR and torus emission on a timescale of weeks to decades, in the NLR emission on a timescale of decades to thousands of years 16 , and in the radio emission on a timescale of a hundred thousand to millions of years. It is notable that, due to their larger scales, the radio components access much longer time scales that the other components.

The different timescales of variability for the emission components on different scales in radio AGN could help to explain the apparent anomaly of the WLRG/FRII sources: perhaps these are objects in which the AGN has switched off or entered a low activity phase for a timescale of hundreds to thousands of years, such that there has been a substantial decrease in the nuclear NLR emission (hence the WLRG classification), but this information has yet to reach the hotspots in the radio lobes of the sources. This possibility has been discussed by Buttiglione et al. (2010) and Tadhunter et al. (2012) Capetti et al. (2011), Capetti et al. (2013) have also considered the possibility that WLRG with unusually low excitation emission line spectra — the so-called extreme low excitation radio galaxies (ELEG) — represent objects caught in the act of the NLR switching off. The advantage of associating the WLRG/FRII sources with the switch-off phase is that it explains their ambiguous status in the unified schemes. Moreover, in the context of explanations of both radio morphological and optical spectroscopic classifications in terms of an Eddington switch, it allows the link between FRII and SLRG, and between FRI and WLRG to be preserved, without the need to explain the WLRG/FRII sources in terms of special accretion flow physics.

It is also possible that the radio AGN classified as relaxed or fat doubles (Owen & Laing, 1989), which lack radio hotspots, but nonetheless have FRII-like radio lobes (e.g. 3C310, 3C314.1, 3C386 in the 3CR sample), represent a post-switch-off stage in the evolution of the radio sources, in which the hot spots have already responded to the decline in jet activity, but the relic lobes have not yet faded below the flux limoit of the sample. Significantly, ∼ 48% of all the 25 WLRG/FRII sources, but only ∼ 11% of all the 80 SLRG/FRII sources, in the 3CR and 2Jy samples of Dicken et al. (2009) and Buttiglione et al. (2009, 2011) could be classified as relaxed or fat doubles. It is also notable that spectral ageing studies of the lobes of at least two WLRG/FRII objects with relaxed radio lobes suggest that they represent relic radio sources in which the jets switched off ∼ 6 – 20 Myr ago (Mazzotta et al., 2004, Harwood, Hardcastle & Croston, 2015).

The switch-off phase must occur at the ends of the activity cycles of all powerful radio sources. However, it is also possible that the nuclear AGN/jet activity is intermittent, and switches off for timescales τnlr < toff < τlobe, within a full activity cycle. In this case, the fraction of FRII that are WLRG sets an upper limit on the proportion of time that the radio sources are “off” for timescales longer than τnlr within the cycle: foff ≤ 0.25 for the full 2Jy and 3CR samples 17 , and foff ≤ 0.1 for the 29 sources at the higher radio power end of the 2Jy sample (P5 GHz > 10 26 W Hz 𢄡 ). These estimates are upper limits on the in-cycle intermittency because some WLRG/FRII sources may represent objects observed in the switch-off phase, and there may be other explanations for the WLRG/FRII phenomenon. We deduce from this that the fuel supply must be steady, and the AGN “on”, for the overwhelming majority of activity cycle. This is remarkable given that the rate of infall of gas to the nuclear regions may be clumpy rather than smooth, and the feedback effect of the jets and winds associated with the AGN may disrupt the fuel supply.

Note that, based on the existing spatially integrated spectra of radio galaxies, it is not possible to rule out the existence high-amplitude intermittency on shorter timescales (toff < τnlr). Indeed, Inskip et al. (2007) have argued that the 2Jy BLRG PKS1932-46 represents a case in which the AGN has recently entered a low activity phase, but this is not yet reflected in its NLR properties, because its [OIII] and mid-IR narrow emission line luminosities are much higher than expected from its relatively low X-ray, near-IR and mid-IR continuum luminosities, and Hα broad-emission-line luminosity. There are also several reports in the literature of 𠇌hanging state” AGN in which broad lines near-disappeared, or appeared, on timescale of years to decades (e.g. Penston & Pérez, 1984, LaMassa et al., 2015, MacLeod et al., 2015, and references therein).

Before concluding that all WLRG/FRII sources represent powerful radio AGN in which the nuclear activity has switched off, either permanently of temporarily, it is important to consider some objections to this interpretation. First, it has been noted that the radio cores of WLRG/FRII sources are stronger than expected in the case that the nuclear activity has switched off (Buttiglione et al., 2010). Second, some WLRG/FRII appear to be in richer environments and have more luminous host galaxies than their SLRG/FRII counterparts, suggesting that the gaseous environments of the radio sources may play a role in determining whether a source appears as a WLRG/FRII (Hardcastle, Evans & Croston, 2007, Ramos Almeida et al., 2013, Ineson et al., 2015). However, neither of these objections is insurmountable. For example, in the case of the radio cores, at least some of the WLRG/FRII sources have radio cores that are significantly weaker relative to their extended radio emission than in SLRG/FRII sources examples from the 2Jy sample include PKS0043-42, PKS0347+05, PKS1648+05 and PKS2211-17, all of which have core-to-extended radio flux ratios measured at 2.3 GHz of R2.3 GHz < 10 𢄣 compared with a median of R2.3 GHz = 3 × 10 𢄣 for the NLRG/FRII sources in the same sample (Morganti et al., 1997a). In addition, FRI radio galaxies have higher core/extended radio flux ratios on average than FRII sources (Morganti et al., 1997a), so if a SLRG/FRII source were to drop down to a lower level of nuclear activity consistent with it eventually becoming a WLRG/FRI source, this change would not necessarily have a dramatic effect on the radio core flux.

Considering the environments and host galaxies, it is notable that, even if WLRG/FRII sources are associated with richer environments and more luminous host galaxies on average than SLRG/FRII sources, there is a considerable range in, and overlap between, the environments and host galaxy properties of the two groups (Ramos Almeida et al., 2013, Ineson et al., 2015). It is also possible that the level of high amplitude intermittency depends on the environment and host galaxy properties (e.g. due to the nature of the fuel supply), such that FRII sources in rich environments and with massive host galaxies are more likely to be observed in a low activity (WLRG) phase. Alternatively, the confinement effect of the dense hot gaseous haloes in rich galaxy environments may lead to the relic radio sources remaining visible for longer than they would in lower density environments (see discussion in Murgia et al., 2011), thus making in more likely that WLRG/FRII sources will be observed in rich environments.

So far, this review has concentrated on the results for samples of nearby radio sources selected in bright, flux-limited radio surveys such as the 3CR and the 2Jy. Because such surveys select relatively high luminosity radio sources, which are rare in the local universe, the sample sizes are modest. However, cross correlation of the SDSS optical and FIRST/NVSS radio catalogues has recently allowed the selection of much larger samples of nearby radio sources (z < 0.2) down to lower flux densities and radio luminosities (Best et al., 2005, Best & Heckman, 2012), albeit lacking the detailed radio and optical morphological, X-ray and infrared spectral information available for the brighter samples. It is interesting to consider how the sources selected in the bright radio samples compare with those in the SDSS-selected samples.

The homogeneous spectral data available for these SDSS-selected samples has allowed the radio luminosity functions to be derived separately for the LEG and HEG populations. The results are shown in Figure 11, which is taken from Best & Heckman (2012). While the HEG population starts to dominate the radio source population at the highest radio luminsoities, and the LEG population dominates at lower radio luminosities, it is interesting that both types of objects are present in the SDSS sample at all radio luminosities.

Although detailed radio morphological information is lacking for the SDSS-selected samples, based on the results for the 3CR and 2Jy samples, it seems highly likely that a large fraction of the LEG population at higher radio luminosities (L1.4 GHz > 10 24 W Hz 𢄡 ) are WLRG/FRI, with perhaps an admixture of WLRG/FRII sources, whereas the majority of the high radio power HEGs are SLRG/FRII sources. Moreover, the majority of these high radio luminosity sources are likely to be hosted by giant elliptical galaxies (see section 3.1 below). However, considerable uncertainty surrounds the natures of the HEG and LEG radio sources and their host galaxies at lower radio luminosities: for a given optical spectral classification, do the lower power radio sources have the same radio and optical morphologies as their higher radio power counterparts? Do the luminosity functions for the HEGs and LEGs each represent uniform populations of objects?

Concentrating first on the LEG objects, there is already evidence that the natures of the radio sources associated with the LEGs change dramatically with radio luminosity. For example, a recent radio morphological study of the low radio luminosity end of the LEG population (L1.4 GHz < 10 24 W Hz 𢄡 ) by Baldi et al. (2015) shows that most of the lower radio luminosity sources are compact and lack the prominent jets and diffuse lobes typical of more powerful LEG objects. Baldi et al. (2015) have labelled such sources FR0 radio galaxies, in recognition of the fact that they fall outside the usual FRI/FRII classification scheme (see also Ghisellini, 2011, Sadler et al., 2014). However, the host galaxies of most of the FR0 objects appear to be massive, early-type galaxies, just like their more powerful LEG counterparts.

The situation for the HEG sources is less clear. While the high radio luminosity end of the HEG population is dominated by FRII radio sources that are hosted by giant elliptical galaxies, the lower power end of the luminosity function of the HEGs overlaps with the high power end of the radio luminosity function of Seyfert galaxies (Meurs & Wilson, 1984), which are more likely to be hosted by spiral galaxies. This raises intriguing questions about the makeup of the HEG population: is there a gradual transition between spiral and elliptical host galaxies from low to high radio powers along the HEG sequence? How do the radio morphologies change along the HEG sequence?

2 See http://www.jb.man.ac.uk/atlas for a more detailed discussion of radio classifications schemes. Back.

3 Since the GPS and related CSS sources show strong morphological similarities with their more extended counterparts (O'Dea, 1998, Tzioumis et al., 2002), and are thought to represent radio sources in a young evolutionary phase (Fanti et al., 1995, Owsianik, Conway & Polatidis, 1998, Polatidis & Conway, 2003), in most of this review I will not distinguish them from the extended sources. However, note that the detection of GSP/CSS sources in flux-limited samples may be affected by selection effects related to the strong interactions between the compact radio sources and dense ISM in the central regions of the galaxies (Tadhunter et al., 2011, Morganti et al., 2011, Dicken et al., 2012). Back.

4 N.B. This object would also be classified as a BLRG on the basis of its broad Balmer emission lines. Back.

5 This includes the 3 objects classified by Capetti et al. (2011) as extreme low excitation radio galaxies (ELEGs). Back.

6 Buttiglione et al. (2010) and Capetti et al. (2013) also distinguish a class of radio AGN showing extremely low excitation emission line spectra that they label extreme low excitation galaxies (ELEGs). However, in what follows I make no distinction between LEGs and ELEGs and label them collectively as LEGs. Back.

7 Note that this sample extends to lower redshifts than the 0.05 < z < 0.7 2Jy sample of Dicken et al. (2009) described in Table 1. Back.

8 One object deserves particular mention here: 3C84 (also known as NGC1275). This object was originally classified as an FRI by Fanaroff & Riley (1974), but shows broad Balmer emission that lead to a BLRG optical classification. Also, as expected for a BLRG, it has luminous narrow lines, albeit of low ionization (hence the LEG classification of 3C84 by Buttilgione et al. 2010). On the basis of these properties, 3C84 could be considered as the only FRI object in the combined 3CR and 2Jy sample classified as a SLRG. However, its radio structure is highly peculiar, with a strong, highly variable flat spectrum core, an inner steep spectrum double structure, and an outer halo (Pedlar et al., 1990). Therefore, its radio morphological classification must be regarded as uncertain it is certainly not a typical FRI source. Back.

9 Note that Haas et al. (2005) derive a larger attenuation factor of ∼ 7 based on a smaller, more heterogeneous sample Back.

10 Note that the mid-IR AGN detection rates are based on photometric measurements with the IRAC instrument on the Spitzer satellite. In this case, rather than using high spatial resolution to detect the compact nuclei, the AGN are detected as an excess in the mid-IR continuum over the flux in the starlight predicted on the basis of an extrapolation of the near-IR starlight flux, with the latter measured from HST observations using the same aperture as used for the Spitzer/IRAC measurements (see Ramírez et al., 2014a). Back.

11 Using updated optical spectral classifications from the data in Buttiglione et al. (2009, 2010) and other sources. Back.

12 The exceptions are WLRG that have a strong starburst heated dust component or substantial contamination of their far-IR emission by non-thermal jet emission. Back.

13 This is consistent with the idea that the relatively low equivalent widths of the optical emission lines in BL Lac objects is due to their strongly beamed optical continuum emission, rather than their low emission line luminosities. Back.

14 Discounting the 3C84, which, as noted above, has a peculiar radio morphology. Back.

15 This is the typical scale on which the narrow-line emitting gas accessed by the nuclear apertures used for the spectroscopic classification is emitted. Many radio galaxies show lower-surface-brightness extended emission line nebulosities on much larger scales (up to ∼ 100 kpc: Baum et al., 1988, Tadhunter, Fosbury & Quinn, 1989a), but these are not used for the optical spectral classifications of the sources. Back.

16 Note that, for a typical cloud in the NLR, the time taken for the [OIII] emission to fade following a substantial decline in the ionizing flux is negligible compared with the light crossing time of the NLR (Capetti et al., 2013). Back.

17 Note that WLRG with hybrid FRII/FRI radio morphologies were counted as WLRG/FRII sources when making this estimate. Back. *****


4. Baldwin Effect

4.1. Wλ vs. Lc

The relation between the emission line equivalent width (Wλ) and the continuum emission luminosity (Lc) is given by Baldwin (1977),

The slope β is found to be negative for RQ AGN (Baldwin, 1977 Shields, 2007), i.e., the equivalent width of the emission line (or the contrast between the line and continuum luminosities) decreases toward large continuum luminosities.

Equation (1) can be transformed to a relation between the total line luminosity Lline and the monochromatic continuum luminosity λLλLc measured at a certain wavelength λ,

4.2. Contribution of Non-thermal Emission to the Baldwin Effect

The optical continuum emission in RQ AGN is assumed to be isotropic and generated in the accretion disk, so that the continuum luminosity is L c = L disk RQ . On the other hand, for FSRQ the optical continuum emission has two components, the thermal emission from the disk ( L disk BL ) and the beamed non-thermal emission from the relativistic jet (Ljet), i.e., L c = L disk BL + L jet . It was assumed that the main contribution to the broad line emission is attributed to the disk thermal emission, while the beamed emission from the jet is produced beyond the BLR and, hence, has no contribution to the Broad Line (BL) emission.

A simulation was performed in order to estimate the statistical properties of the RQ samples studied by Greene and Ho (2005) and Shen et al. (2011). The purpose is to compare them to the RL sample presented in this paper. First. a distribution of continuum luminosity was taken (see Figure 1), based on the luminosity distributions observed for our sample of FSRQ. For all three lines a gaussian distribution represents well the data. Then, using the line luminosity - continuum luminosity relations described in the afore mentioned papers, and the scatter obtained for these relationships, representative line luminosities of the sample of RQ AGN were generated. Using these line luminosities and the assumed distribution of continuum luminosities, the equivalent widths for the simulated RQ sample were calculated. The number of simulated values in each case is equal to the number of data points in the original RQ samples. Figure 2 shows the simulation results for the Hβ, Mg II, and CIV lines.

Figure 1. Distributions of continuum luminosity at 3000 Å obtained for the observed FSRQ sample (left panel) and for the simulated RQ sample (right panel).

Figure 2. Relations obtained for the Hβ (left panel), Mg II (middle panel), and CIV (right panel) emission lines. Top row: Lline-Lc plot with simulated data. Bottom row: EW-Lc plot with simulated data.

After generating the data, it was separated in order to match the continuum luminosity range on our sample of FSRQ. Using the simulated data that fall in our observed range, the mean and standard deviation of the equivalent width for the three lines Hβ, C IV and Mg II, were calculated. It is worth mentioning that for C IV and Mg II, all the simulated data fall inside the desired continuum luminosity ranges.

Then, a linear least-squares algorithm in one dimension was applied, using the IDL task FITEXY 6 , to perform a linear fitting to the Lc - Wλ relation to obtain the Baldwin Effect of the RQ sample.

From these simulations, the next conclusions were drawn:

• The slope and uncertainty on the simulated Lc - Wλ relation remains unchanged, regardless of the input continuum luminosity distribution used.

• The changes on the intercept and uncertainty on the simulated Lc - Wλ relation, with changes on the input continuum luminosity distribution are negligible.

• The mean and standard deviation of the simulated equivalent widths can change drastically, depending on the continuum luminosity distribution used as input.


Radio loud AGN restarting as Radio Quiet AGN? - Astronomy

Active Galactic Nuclei (AGN) is a term encompassing a variety of energetic phenomena in galactic centers which are thought to be powered directly or indirectly by accretion of matter onto central supermassive black holes. There are many reviews of this field, including several of book length. Perhaps the most technical and broad book-length review is that of J. Krolik, which was published in 1998. The classic text on spectral analysis of gaseous nebulae and AGN has been updated in 2006 (Ferland and Osterbrock) it's a wonderful book, but it is biased towards the optical and ultraviolet regions of the spectrum. A new edition (2009) of An Introduction to Radio Astronomy by Burke and Graham-Smith is also noteworthy.

Few theoretical predictions have been borne out in this field, and understanding is still semi-quantitative at best (Alloin et al. 1985 Antonucci 1988, 2002a Courvoisier and Clavel 1991 Koratkar and Blaes 1999 Blaes 2007). This review will refer to some fairly general theoretical ideas, but it will mainly organize some observational information on radio galaxy central engines that has become clear over recent years. It will not include a general introduction to AGN (see above references), or even to orientation effects in AGN, but will address the nature of the central engines in radio galaxies via thin electromagnetic output, which is a topic tied up as a practical matter with orientation effects on observations.

The general topic of orientation effects ("Unified Models") is reviewed in detail in Antonucci 1993. The material in that review is almost entirely "still true." However, it has been updated and elucidated in several more recent (but generally narrower) reviews (Urry & Padovani 1995 Dopita 1997 Cohen et al. 1999 Wills 1999 Axon 2001 Tadhunter 2008).

In a nutshell, prior to the mid-1980s, it seemed that radio loud Quasars (with their powerful thermal optical/UV light) were quite distinct from radio galaxies the latter are analogous to Quasars in general radio properties, but apparently lacked the strong optical/UV electromagnetic luminosity (Big Blue Bump). In a landmark review of bright extragalactic radio sources, Begelman, Blandford and Rees (1984) posited that "The ratio of mass accretion rate to the mass of the hole may determine whether a [radio loud] active galactic nucleus will be primarily a thermal emitter like an optical Quasar or a nonthermal object like a radio galaxy."

Through optical spectropolarimetry and other means, it was subsequently determined that many radio galaxies, especially the most powerful ones, with strong high-ionization narrow emission lines, actually harbor hidden Quasars surrounded by opaque dusty tori, so that the observational appearance depends on the inclination of the radio axis to the line of sight. But now we know that for many radio galaxies, any hidden Quasar must be very weak.

The radio quiet and radio loud objects with high ionization are virtually all visible (Type 1) or hidden (Type 2) Seyferts or Quasars (e.g., Antonucci 2002b). At lower radio luminosities of all radio types, we find mostly LINERs 6 (Low Ionization Nuclear Emission Regions). A recent and comprehensive review of LINERs is that of Ho (2008).

In the low-redshift universe, there is a near-perfect correspondence between "radio loud" objects - with L (1.5 GHz) loosely extending from perhaps

10 28 - 10 36 erg sec -1 Hz -1 - and Elliptical hosts. (This paper uses H0 = 70 km sec -1 Mpc -1 , matter = 0.3, and = 0.7.) For a fiducial L

-1 spectrum, for which there is equal power per logarithmic frequency interval - also referred to as "per dex" - the parameter L gives the power integrated over an interval of 0.30 dex (powers of 10). Thus if one integrates such a spectrum over the "radio region" 30 MHz-300 GHz, the corresponding luminosity is 13.3 times as great a look at various radio AGN in the NASA Extragalactic Database (http://nedwww.ipac.caltech.edu) shows that the luminosity per dex of the dominant optically thin synchrotron component tends to be lower outside this range, so depending on the application, one may consider this as a crude "radio bolometric correction." Resulting radio powers can exceed 1 × 10 45 erg/s. Estimates of energy tied up in 100-kpc scale radio lobes are as high as 1 × 10 61 ergs or more, even using the particle/magnetic field minimum-energy assumption and assuming a lack of a dominant proton contribution. (The minimum-energy assumption posits that energy is apportioned between relativistic electrons and magnetic field in such a way as to minimize lobe energy content for a given synchrotron luminosity.)

Unified models assert that certain AGN classes differ only in orientation with respect to the line of sight. These models comprise two separate (though interacting) assertions, as illustrated in Fig. 1. The first to be recognized historically is the effect of relativistic beaming (aberration causing anisotropy in the observed frame) in the powerful synchrotron jets which feed particles and magnetic energy into the radio lobes. When seen at low inclinations, beaming amplifies and speeds up "core" (subparsec scale jet, usual synchrotron-self-absorbed) radio flux variability and (apparent faster-than-light) "superluminal motion." Thus a special fortuitous orientation of a nearly axisymmetric object leads to a different observational category (Blazars). The same objects, seen at higher inclination, are ordinary radio-loud galaxies and Quasars (Blandford et al 1984 Antonucci and Ulvestad 1985 Kollgaard et al 1992, etc.). We can call this the beaming unified model (or more loosely, just the beam model, though that term usually connotes some connection with actual physics).

At low redshift, radio-quiet (but not silent) AGN lie in Spiral hosts, and go by the name of Seyfert galaxies. They rarely show detectable motions in their weak radio jets - and when they do show motions, the apparent speed is usually much less than the speed of light. More luminous radio quiet objects are called Radio Quiet Quasars, or historically, Quasistellar Objects (QSOs).

But both radio quiet and many radio loud AGN widely exhibit another kind of orientation unification: many well-studied objects include energetically dominant continuum components in the optical-ultraviolet region, referred to as the Big Blue Bump, and widely attributed to thermal radiation from optical thick accretion flows. (Some confirmation of the latter can be found in Kishimoto et al 2004, 2005, 2008.) They are almost always accompanied by broad (5000-10,000 km/s) permitted emission lines. Both these components reside inside optically opaque dusty structures which to zeroth order have the shadowing properties of tori, and these structures are referred to loosely as the "AGN torus." In some cases we have a direct (low-inclination, polar) view of these compact components (in Quasars, Seyfert 1 galaxies, and Broad Line Radio Galaxies). In many other cases (high inclination, equatorial views), we can be sure these two components are present inside the tori using the technique of optical spectropolarimetry, which literally allows us to see the nuclei "from above," using ambient gas and dust as natural periscopic mirrors.

This wonderful trick of optical spectrapolarimetry (e.g. Antonucci 1982, 1983, 1984 Antonucci and Miller 1985) uses the polarization property of scattered light to separate the spectrum of hidden sources from any sources of direct light. The above references show that the Big Blue Bump and the Broad Line Region are present but hidden from direct view. The scattering polarization position angle indicates that the photons emerging from these components can only escape from the hidden nuclei if moving along the radio jet (and lobe) axis, so that the other directions must be blocked by an obscuring torus. This is illustrated in Fig. 2, from Tran et al 1995, illustrating the total and polarized fluxes, for 3CR234. The polarized photons (scattered light) contribute to both plots of course, but show up with great contrast when plotted alone.

If all objects fit in with the above descriptions of orientation unification, we would be left with as few as two independent types of "central engines" on relativistic scales, those for radio quiet and those for radio loud AGN. But we now know that many radio galaxies lack powerful visible or hidden Type 1 engines (Big Blue Bump, Broad Line Region, copious accretion flows). They occur frequently in some regions of parameter space (as delineated for example by redshift and radio flux or luminosity), and not in others. The purpose of this review is to gather the multiwavelength evidence for the last two statements, and emphasize the dependence of the distribution of central engine types on parameter space, which is the key to avoiding many errors and much confusion and even discord. We will see that a great deal of self-consistent information is available on the occurrence of the two types of radio galaxy/Quasar central engines.

Within the radio loud AGN arena, several types are also distinguished based not on central engine properties but on extended radio morphology. Unfortunately, they do not correspond very closely to the optical catagories!

The most luminous giant ( 100 kpc) double radio sources are designated FR II ("Classical Double"), for Fanaroff and Riley (1974) those authors noted that a whole suite of properties change together fairly suddenly over a critical radio luminosity (referring to the roughly isotropic diffuse emission). Sources above

2 × 10 32 erg/sec/Hz at 1.5 GHz show edge-brightening and hot spots, where the radio jets impinge on an external medium, and shocks partially convert bulk kinetic energy to particle and field energy. They also tend to have strong side-to-side asymmetry (generally attributed to relativistic beaming) of the jets over scales from the relativistic region up to tens of kpc. The lower luminosity giant (also 100 kpc) objects (FR I galaxies) also have strong side-to-side jet asymmetry, but only on 1-1000 pc scales in most cases. An important refinement to the FR classification scheme is that the dependence of the exact radio luminosity cutoff depends on the optical luminosity of the host galaxy (Owen and Ledlow 1994, but see Best 2009, Fig. 4a).

The "FR" types of radio galaxy generally have sizes of 25-1000 kpc, but large populations of smaller sources exist, and they can still be very powerful. (The small sizes mean that their lobe energy content is very much lower, however.) They are denoted in an inconsistent way, according to their means of discovery. Optically thin, steep-spectrum radio sources which were historically unresolved on arcminute scales, were (and are) called Compact Steep Spectrum (CSS) sources to distinguish them from opaque beamed synchrotron cores their spectra peak at

100 MHz, and they are generally defined to be in the size range 1 kpc - 15 or 25 kpc. Sometimes they are crudely defined to be sources smaller than typical host galaxies. Sources whose extent is less than

1kpc, with even more compact substructure, are often dominated by synchrotron components which are self-absorbed up to

GHz frequencies 7 , and are called Gigahertz-Peaked-Spectrum sources (GPS). (Even tinier sources are being sought by selecting for self-absorption peaks at even higher frequencies.) These classes are reviewed by O'Dea (1998). Since that paper was written, much evidence has accumulated from VLBI proper motions that the sources are small because they are very young (

1000-100,000 years! 8 ). In at least some cases it is known from faint extended emission that these very young ages refer only to a recent phase of activity however. Statistically, only a very small fraction of small and very short-lived sources can grow to be huge bright long-lived sources.

The relevant properties of classes will be discussed in turn, generally starting with the radio data and proceeding upwards in frequency. We will discover that some FR II radio galaxies at z 0.5-1.0 lack powerful hidden Quasars. These objects may have hidden Type 1 nuclei, but they are constrained to be much weaker than those of the "matched" 9 visible Quasars, and thus they are not "unified" (identified) with them via orientation with respect to the line of sight. At low redshifts (z <

0.5) it is probable that only a minority of FR II radio galaxies in the 3CR catalog host hidden Quasars.

Next we will tackle the less powerful (FR I) giant radio galaxies, which are by selection nearby in almost all cases. The 3CR catalog flux cutoff of 10 Jy at 178 MHz (Laing, Riley and Longair 1983) corresponds to the FR I vs. II radio luminosity separation at z

0.2. Most (but by no means all) of these objects have nuclear spectral energy distributions dominated by synchrotron radiation, with no evidence for visible or hidden "Type 1" central engines.

Finally the small, young GPS and CSS sources will be discussed. Very recent information from the ISO and Spitzer infrared satellites has greatly increased our knowledge of "shadowing unification" 10 at various radio luminosities.

A major caveat of this paper, and of this field, is that most of the information derives from the brightest radio sources, especially those in the 3CR catalog, so no implication is made for unexplored regions of parameter space! Another major caveat is that while we discuss thermal vs. nonthermal galaxies 11 , relatively little evidence is presented that any parameter is bimodal, so that there could be a continuum of properties.

Radiation absorbed by the dusty torus is largely reradiated as infrared, and many studies have concluded that in reasonably luminous AGN (so that the IR is not dominated by a normal host galaxy), at least the near- and mid-infrared reradiation ( 1-40 microns) is dominated by reprocessed nuclear optical/UV/X-ray light. In specific populations, the entire IR seems to be radiated by the torus, because the colors are warm throughout, and there is evidence for only weak star formation (e.g. PAHs) or synchrotron radiation (e.g., strong radio-mm emission). Thus the infrared reradiation of nuclear light can potentially be used as a calorimetric indicator for the luminosity of any hidden AGN.

There are at least two ambiguities in using this infrared emission as an AGN calorimeter. The first is that some of the nuclear radiation may not reach the torus, either because of the intrinsic latitude-dependence expected for many models of the Big Blue Bump (e.g. Netzer 1985), or preferentially planar Broad Line Region absorption (Maiolino et al 2001c Gaskell et al 2007). Both of these effects add noise to the dust reradiation calorimeter, and tend to make hidden Quasars look dimmer than visible Quasars in the infrared for a given opening angle. Nevertheless most studies show infrared luminosity about as expected, from detailed studies of individual objects (Carleton et al 1984 Storchi-Bergmann et al 1992) and for populations of Type 1 and Type 2 objects within isotropically selected samples (Keel et al 1994). Carleton et al (1984, see Fig. 3) shows how the infrared calorimeter works, based on the first spectropolarimetric hidden AGN, 3CR234.

Also, although the covering factor deduced from dividing infrared re-emission by Big Blue Bump luminosities are therefore lower limits, they tend to be high (

0.1-1) so that they can't be too far off! 12

Another concern with the infrared calorimeter is the expected anisotropy of the thermal dust emission due to the large dust column densities (Pier and Krolik 1992, 1993). Many Type 2 AGN have X-ray columns of 1 × 10 24 cm -2 absorption of mid-IR lines, molecular maps, and the great difference in the average X-ray columns between Type 1 and Type 2 AGN suggest that a commensurate dust extinction is present. (See Maiolino et al 2001a, b, c for arguments which affect this line of reasoning quantitatively but not qualitatively.)

There is a limit on the anisotropy of the ratio [O III] 5007 / F (60µ) in Seyferts from Fig. 3 of Keel et al 1994. The figure shows that in their well-selected sample (60µ flux with a mild 25µ - 60µ warmth criterion), Type 1 and Type 2 objects have indistinguishable distributions of L(60µ), L[5007], and of course their ratio. The 5007 line is produced outside the torus in most objects like these: it doesn't appear in polarized flux along with the broad emission lines and Big Blue Bump. Thus it's not significantly hidden inside the torus, is likely quite isotropic in this parameter space, and is nearly isotropically selected 13 , fairly powerful Seyferts.

Torus models do predict that the optical depths will be small in the far-infrared in general. For AGN-dominated infrared SEDs, that means there is an elegant and detailed method of deriving the degree and wavelength-dependence of the dust emission anisotropy. For isotropically selected samples, we can divide composite or representative Type 1 SEDs by those for Type 2, tying them together at 60µ.

In order to make the infrared calorimeter accurate, we need to account for the common anisotropy of the near- and mid-thermal dust emission. There is general agreement on near isotropy past

30µ, as was first predicted by Pier and Krolik (1992, 1993). The main basis for their prediction of anisotropy at shorter wavelengths was that many Seyfert 2's are Compton-thick, including NGC 1068 - which is completely opaque and for which N(H) 10 25 according to the X-ray spectrum (e.g., Pounds and Vaughan 2006). In Galactic dusty gas this corresponds to A(V)

1000. (See e.g., Maiolino et al 2001a, b for a quantitative correction, which however doesn't greatly affect the discussion of gas vs. dust columns below.) As an aside, I do think that dust-free atomic gas can contribute to the X-ray absorption in some cases (e.g., Risaliti et al 2011 Antonucci et al 2004), but the X-ray columns of Type 2 AGN are on average 100 × those of Type 1, so most of the column is connected with the Type 2 classification. Also in some cases, the dust column can be constrained to be similar to the very high X-ray columns. For example, Lutz et al (2000) used the lack of Pf- at 7.46µ to show that A(V) > 50 in NGC 1068.

Recall that one could pick spectral Type 1 (broad-line) radio loud Quasars and radio galaxies by some fairly isotropic AGN-related luminosity (e.g., lobe-power), and divide the two spectra to produce a spectrum of diminuation from anisotropy. We were able to do this for the z > 1 3CR (Fig. 4) because all of the radio galaxies seem to host hidden Quasars, but only out to

15µ in the rest frame. The correction for anisotropy is around a factor of 10 in the near-IR, 1.5-2 on both sides of the silicate feature, and about three in the silicate feature, but varies significantly from object to object. As expected, it flattens out at

Recall that torus shape was inferred from the high 14 typical polarization of the reflected light in hidden-Broad Line Region objects, which is generally perpendicular to the radio axis, meaning that photons can only escape the nucleus to scatter into the line of sight if they leave the nucleus in the polar directions.

The obscuring dusty tori invoked in the shadowing aspect of the unified model are active, not just passive, components. Modulo factors of order unity for geometry and dust cloud albedo, the tori will reprocess almost all of the incident Big Blue Bump/Broad Line Region luminosity into the infrared. Thus the ratio of re-emitted infrared emission to Big Blue Bump (and often considerable absorbed X-ray) emission tells us the approximate covering factor of the dusty gas which is idealized as a torus shape. 15

The covering fraction should agree with the fraction of Type 2 (hidden) nuclei in a sample which is selected by any isotropic AGN property. What are the results of this comparison? I'm crudely leaving out intermediate types, and the substantial minority of Seyfert 2s whose Type 1 nuclei are hidden by dust in the host galaxy plane (Keel 1980 Lawrence and Elvis 1982) rather than by a nuclear torus. Keel et al (1994) found 80 far-IR selected Seyfert 1s and 141 Seyfert 2s (plus some H II galaxies that they weeded out). So one expects a dust covering factor of

two-thirds in this parameter space. While the paper describes dust and gas covering factors as "broadly consistent" with the model, it would still be very valuable at this point to take advantage of this sample and actually measure covering factor for the Type 1s based on existing spectral energy distributions.

Figure 4 (kindly supplied by C. Leipski and S. Hoenig) shows the quotient SED for z > 1 individual Quasars and the radio galaxy composite, which however only covers up to 15 microns in the rest frame. 16 Again this plot purports to give us directly the anisotropy as a function of wavelength.

6 A few LINERs have strong broad lines, e.g., Filippenko and Halpern 1984. Many LINERs have very inconspicuous broad H components, but it isn't clear to me that they are strictly analogous to those in low-luminosity Seyferts. For example, we do not know whether they vary rapidly (L. Ho, 2011, private communication). Back.

7 It is proposed by Begelman (1999) that the weak fluxes at low frequencies result from free-free rather than inverse synchrotron absorption. Encouraging follow-up work can be found in Stawarz et al 2008 and Ostorero et al 2010. Back.

8 This material was reviewed recently by Giroletti 2008. Back.

9 In radio flux and redshift. Back.

10 As a reminder, unification of broad line and narrow line AGN by orientation of a toroidal nuclear obscurer is here called "shadowing unification." Ascription of superluminal motion and other relativisitic effects in subpopulations to orientation is called "beaming unification." Back.

11 Recall that thermal vs. non-thermal refers not to the radio emission itself, but to the presence of an energetically dominant optical/UV source thought to arise from an accretion disk. Back.

12 Although Maiolino et al 2007 didn't actually integrate over the SED and should therefore be viewed with much caution, it is interesting that they actually find some covering factors nominally greater than one for low-moderate luminosity objects. This was predicted qualitatively by Gaskell et al (2004), who argued that a preponderance of large grains leads to less apparent absorption (spectral curvature) in the optical/UV region than otherwise, and thus one may fall into the trap of inferring very low extinction. The inner torus is expected to lack small grains, most robustly perhaps on the grounds that the torus sublimation radius must be larger for small grains according to radiative equilibrium. The lack of small grains at the sublimation radius is supported by the data of Suganuma et al 2006: see discussion in Kishimoto et al 2007. Back.

13 Isotropic selection means selection on an isotropic AGN property, which avoids powerful biases in comparisons between classes. It is essential in order to produce intelligible results. Back.

14 A common error is to use the percent polarization of the continuum (after correcting for bulge or Elliptical starlight) as the scattering polarization: that is usually contaminated by hot stars, and only the broad lines themselves can be used to find, or more often place a lower limit on, the percent polarization of the scattered nuclear light (Antonucci 2002b). Back.

15 There is at least one obscured Quasar with a thin obscuring disk - the radio-loud, mini-MgII BAL OI287 (Goodrich and Miller 1988 Rudy and Schmidt 1988 Ulvestad and Antonucci 1988 Antonucci, Kinney, and Hurt 1993) however, even that one seems to have an upturn longward of 1µ at least according to the NED figures. Back.

16 Buchanan et al (2006) does this type of division for a 12µ (better than optical, but not ideal) selected sample of Seyfert galaxies. As expected, their anisotropy curve flattens at long wavelengths, but oddly not at a value of unity. Back. *****


Title: FRESH ACTIVITY IN OLD SYSTEMS: RADIO AGNs IN FOSSIL GROUPS OF GALAXIES

We present the first systematic 1.4 GHz Very Large Array radio continuum survey of fossil galaxy group candidates. These are virialized systems believed to have assembled over a gigayear in the past through the merging of galaxy group members into a single, isolated, massive elliptical galaxy and featuring an extended hot X-ray halo. We use new photometric and spectroscopic data from Sloan Digital Sky Survey Data Release 7 to determine that three of the candidates are clearly not fossil groups. Of the remaining 30 candidates, 67% contain a radio-loud (L > 10 W Hz) active galactic nucleus (AGN) at the center of their dominant elliptical galaxy. We find a weak correlation between the radio luminosity of the AGN and the X-ray luminosity of the halo suggesting that the AGN contributes to energy deposition into the intragroup medium. We only find a correlation between the radio and optical luminosity of the central elliptical galaxy when we include X-ray-selected, elliptically dominated non-fossil groups, indicating a weak relationship between AGN strength and the mass assembly history of the groups. The dominant elliptical galaxy of fossil groups is on average roughly an order of magnitude more luminous than normal group ellipticalmore » galaxies in optical, X-ray, and radio luminosities and our findings are consistent with previous results that the radio-loud fraction in elliptical galaxies is linked to the stellar mass of a population. The current level of activity in fossil groups suggests that AGN fueling continues long after the last major merger. We discuss several possibilities for fueling the AGN at the present epoch. « less


Gravitational Lenses

The key to understanding gravitational lensing is to recall that photons follow geodesics and mass curves spacetime. Thus mass concentrations act like optical lenses. When light passes through a lens, the change in the index of refraction between air (or vacuum) and glass causes the light to bend. The index of refraction of a material is the ratio of the speed of light in a vacuum to the speed of light in the material:

Thus n is always greater than one.

To understand how mass causes analogous behavior, recall the expression for the coordinate speed of light in the radial direction for the Schwartzchild metric:

This can be rearranged to give an approximate expression for n:

For a mass of 10e11 Solar masses, and close-approach of

10 kpc, this gives n-1

10e-6. This implies the deviations are small. As they must be (or we would see stuff like this all the time).

If one studies the problem in a way analogous to the treatment of gravitational deflection of starlight by the Sun, one can derive a relationship for an idealized point-lens that indicates there should be two images.

If one examines the quadratic relation for the two image seperations, one finds that the seperations can be used to determine an estimate for the lensing mass (assuming some value for the Hubble constant). Note that if the background source is intrinsically variable, one can in principle then use the observed lag between the variability of the two sources to directly measure the path lengths from us to the source. This means one can work out an estimate for the Hubble constant for such systems. As QSOs are intrinsically variable, this is a method of some interest. Unfortunately, QSO variability is messy. Also, to really make this method work, one must understand the mass distribution of the lens. An extended source lens behaves much differently than an idealized point-source lens.

If the background source is exactly on line with the lens, one gets an Einstein Ring:

In reality, such rings occur whenever the alignment angle is less than the Einstein ring radius (above).


Watch the video: Webinar on Active Galactic Nuclei:The energetic centers of galaxies (January 2023).