The jet production efficiency of radio galaxies can be quantified by comparison of their kinetic jet powers
P
jet
and Bondi accretion powers
P
B
. These two parameters are known to be related linearly, with the jet power resulting from the Bondi power by multiplication with an efficiency factor of order 1%. Using a recently published (Nemmen & Tchekhovskoy 2014) high-quality sample of 27 radio galaxies, I construct a
P
B
−
P
jet
diagram that includes information on optical AGN types as far as available. This diagram indicates that the jet production efficiency is a function of AGN type: Seyfert 2 galaxies seem to be systematically (with a false alarm probability of 4.3 × 10
−4
) less efficient, by about one order of magnitude, in powering jets than Seyfert 1 galaxies, LINERs, or the remaining radio galaxies. This suggests an evolutionary sequence from Sy 2s to Sy 1s and LINERs, controlled by an interplay of jets on the one hand and dust and gas in galactic nuclei on the other hand. When taking this effect into account, the
P
B
−
P
jet
relation is probably much tighter intrinsically than currently assumed.
1. INTRODUCTION
Active galactic nuclei (AGN; see, e.g., Netzer 2013 for a recent review) are supposedly powered by accretion of gas onto supermassive black holes located in the centers of probably all galaxies. A frequent, albeit as yet poorly understood, consequence of this process is the formation of collimated outflows, jets, which are emitters of synchrotron radiation especially prominent at radio wavelengths (see, e.g., Böttcher, Harris, & Krawczynski 2012 for a recent review). One road to a better understanding of jet formation and the interplay of jets and accretion flows is provided by analysis of the energetics of the jets of nearby (
100Mpc) radio galaxies which have accretion rates much smaller (≪1%) than their Eddington limits. X-ray observations of those systems provide independent estimates of (i) the kinetic jet powers
P
jet
and (ii) the Bondi accretion powers
P
B
=
Ṁ
Bc2
, with
Ṁ
B
being the Bondi accretion rate and
c
denoting the speed of light. Multiple studies (Allen et al. 2006; Balmaverde, Baldi, & Capetti 2008; Russell et al. 2013; King et al. 2013; Nemmen & Tchekhovskoy 2014; Fujita, Kawakatu, & Shlosman 2014) concluded that, within errors,
P
jet
is proportional to
P
B
. Accordingly, one finds a jet production efficiency 𝜂 from
P
jet
= 𝜂
P
B
; observationally, 𝜂 ∼ 1% (ensemble average). However, the observed
P
B
−
P
jet
relations show substantial (
0.5 dex) intrinsic scatter, indicating an as yet overlooked systematic effect. This article addresses a possible cause.
2. DATA AND ANALYSIS
I make use of the data set published recently by Nemmen & Tchekhovskoy (2014) [their
Table 1
] which has in turn been drawn from data by Balmaverde, Baldi, & Capetti (2008) and Russell et al. (2013) and which is marked by very careful selection. Jet and Bondi powers were derived from X-ray observations obtained by the Chandra space telescope. Kinetic jet powers were estimated from the ages and creation energies of jetinflated X-ray cavities. Bondi accretion rates were estimated from temperature and density profiles of hot diffuse gas. (Please refer to Nemmen & Tchekhovskoy 2014 for technical details.)
Notes: “Type” refers to optical AGN type; Sy 1/2: Seyfert 1/2 galaxies; L: LINERs; L 1/2: type 1/2 LINERs. References: Source types: NED; SIMBAD; Ho , Filippenko, & Sargent (1997); Véron-Cetty & Véron (2010). Powers: Balmaverde, Baldi, & Capetti (2008); Russell et al. (2013); Nemmen & Tchekhovskoy (2014).
I collected optical AGN types from the NASA/IPAC Extragalactic Database (NED); the SIMBAD Astronomical Database; Ho , Filippenko, & Sargent (1997); and Véron-Cetty & Véron (2010). Eventually, I found optical types for 17 of my 27 targets: two are of type Seyfert 1, six of type Seyfert 2, and nine of type Low-Ionization Nuclear Emission-Line Region galaxy (LINER – cf., Ho 2008). The final dataset is summarized in
Table 1
.
From the data of
Table 1
, I constructed a logarithmic
P
B
−
P
jet
diagram (
Figure 1
) that includes the available information on optical AGN types. The sample spans about three orders of magnitude in both jet power and accretion power. I determined the best-fit linear relation between jet and accretion power via a least-squares fit, leading to a (ensemble averaged) jet production efficiency of 𝜂 ≈ 0.6% for my sample. The intrinsic scatter (i.e., the difference in squares of r.m.s. residual and bivariate r.m.s. measurement error) is about 0.5 dex.
Kinetic jet power Pjet as function of Bondi accretion power PB using the data in Table 1. Optical AGN types are Seyfert 1 (Sy 1), Seyfert 2 (Sy 2), and LINER; RG denotes radio galaxies of unknown type. Error bars indicate 1𝜎 uncertainties. The diagonal continuous line corresponds to a linear relation with a jet production efficiency of 𝜂 ≈ 0.6%.
3. RESULTS AND DISCUSSION
Inspection of
Figure 1
reveals that the sample galaxies follow a linear relation between
P
B
and
P
jet
globally, albeit with substantial systematic scatter. As already noted by Allen et al. (2006), the presence of this relation, as well as the value found for 𝜂, indicate (i) that the accretion flows in the nuclei of radio galaxies are indeed well described by Bondi accretion in general, and (ii) that a significant fraction (about 1%) of the restmass energy of the matter streaming into the Bondi radii of these systems ends up in jets.
Taking into account the optical AGN types, a remarkable pattern emerges: five out of six Seyfert 2 galaxies (the exception being UGC 6297) are located systematically below the best-fit linear relation between
P
B
and
P
jet
, by about 0.8 dex. Noting that most of the remaining galaxies are actually located above the model line, one finds that the five Sy 2s are separate from the remainder of the sample by about one order of magnitude in terms of jet production efficiency. Even though selection effects cannot be ruled out entirely given the limited sample size, this distribution indicates that Seyfert 2 galaxies are physically different from the other sources. This discrepancy is statistically significant: there are six out of 27 galaxies located below the model line by more than their measurement errors in both
P
B
and
P
jet
, with five of them being Sy 2s. When drawing randomly six out of 27 sources, the probability of obtaining a set with at least five Sy 2s is 4.3 × 10
−4
, corresponding to a Gaussian significance of 3.5𝜎.
Within the frame of the standard viewing-angle unification scheme of AGN (Antonucci 1993; Urry & Padovani 1995), the difference between type 1 (Sy 1, LINER 1) and type 2 (Sy 2, LINER 2) galaxies arises from geometry: in case of type 2 AGN, the view into the central engine is obscured by a surrounding extended dust torus, whereas type 1 sources are seen at inclinations that permit a direct view into the central accretion zone. In contrast to this simple geometric explanation,
Figure 1
suggests that the optical type of a given AGN is partially controlled by its jet production efficiency. A possible
evolutionary
description is the following: initially, Seyfert 2 galaxies are AGN which are relatively inefficient in powering jets and whose central regions are heavily obscured by dust. At some point, a sufficient − roughly one order of magnitude − increase in jet power, and thus in the release of energy in the central accretion zone, leads to two effects. Firstly, a partial dissolution of the surrounding dust torus. Depending on the viewing geometry, this may result in the transition of a Sy 2 to a type 1 galaxy; indeed, indications for such an evolutionary sequence were already reported by Tran (2003), Koulouridis et al. (2006) and, more recently, by Villarroel & Korn (2014). Secondly, shockionization of the gas within the galactic nucleus by the jet. As already noted by Heckman (1980), this mechanism provides a natural explanation for the presence of the optical low-ionization emission lines characteristic for LINERs (see also Reynaldi & Feinstein 2013; Couto et al. 2013; Riffel et al. 2014 for recent observational evidence for shock ionization by jets). Eventually, the scenario I propose implies an evolution from Seyfert 2s to LINERs and Seyfert 1 galaxies caused by an increase in jet production efficiency. I note that the rather clear separation of the Seyfert 2s from the rest of the sample suggests that the transition in efficiency happens relatively fast compared to the characteristic life times of galactic radio jets, meaning
10
6
years (cf., e.g., Tadhunter et al. 2012).
When taking into account the apparent dependence of jet production efficiency on optical AGN type, the
P
B
−
P
jet
correlation could actually be much tighter than usually assumed. Indeed, when removing the six Sy 2 as well as the two outliers M 87 and NGC 4696 from the sample, I find an efficiency 𝜂 ≈ 1% and an intrinsic scatter about the best-fit model line
consistent with zero
. A priori, some intrinsic scatter, as well as outliers, have to be expected because (i) the efficiency of accretion flows partially depends on the angular momentum of the infalling gas and the geometry and strength of the magnetic fields in the central accretion zone (e.g., Narayan & Fabian 2011), (ii) AGN activity shows strong red-noise type temporal variability (e.g., Park & Trippe 2012, 2014; Kim & Trippe 2013), and (iii) a dependence of the AGN power on the physical properties of host galaxies, nearest-neighbor galaxies, and the large-scale background densities (Choi, Woo, & Park 2009; Hwang et al. 2012). Those effects seem to be less influential in general than currently assumed once AGN evolution is taken into account properly.
Acknowledgements
This work made use of the NASA/IPAC Extragalactic Database, the SIMBAD Astronomical Database, and the software package DPUSER developed and maintained by Thomas Ott at MPE Garching. I acknowledge financial support from the Korean National Research Foundation (NRF) via Basic Research Grant 2012-R1A1A2041387.
Bottcher M.
,
Harris D. E.
,
Krawczynski H.
2012
Relativistic Jets from Active Galactic Nuclei
Wiley-VCH
Weinheim
Choi Y.-Y.
,
Woo J.-H.
,
Park C.
2009
Environmental Dependence of Active Galactic Nucleus Activity. I. The Effects of Host Galaxy
ApJ
699
1679 -
DOI : 10.1088/0004-637X/699/2/1679
Couto G. S.
2013
Kinematics and Excitation of the Nuclear Spiral in the Active Galaxy Arp 102B
MNRAS
435
2982 -
DOI : 10.1093/mnras/stt1491
Fujita Y.
,
Kawakatu N.
,
Shlosman I.
2014
AGN Jet Power and Feedback Controlled by Bondi Accretion in Brightest Cluster Galaxies
MNRAS
1
Heckman T. M.
1980
An Optical and Radio Survey of the Nuclei of Bright Galaxies
A&A
87
152 -
Ho L. C.
,
Filippenko A. V.
,
Sargent W. L.
1997
A Search for “Dwarf” Seyfert Nuclei. III. Spectroscopic Parameters and Properties of the Host Galaxies
ApJS
112
315 -
DOI : 10.1086/313041
Kim J.-Y.
,
Trippe S.
2013
How To Monitor AGN Intra-Day Variability at 230 GHz
JKAS
46
65 -
Koulouridis E.
2006
Local and Large-Scale Environment of Seyfert Galaxies
ApJ
639
37 -
DOI : 10.1086/498421
Nemmen R. S.
,
Tchekhovskoy A.
2014
On The Efficiency of Jet Production in Radio Galaxies
MNRAS
1
13 -
Netzer H.
2013
The Physics and Evolution of Active Galactic Nuclei
Cambridge University Press
Cambridge
Reynaldi V.
,
Feinstein C.
2013
Shock Ionization in the Extended Emission-Line Region of 3C 305: the Last Piece of the (Optical) Puzzle
MNRAS
435
1350 -
DOI : 10.1093/mnras/stt1377
Riffel R. A.
2014
Feeding versus Feedback in NGC 1068 Probed with Gemini NIFS - I. Excitation
MNRAS
442
656 -
DOI : 10.1093/mnras/stu843
Russell H. R.
2013
Radiative Efficiency, Variability, and Bondi Accretion on to Massive Black Holes: the Transition from Radio AGN to Quasars in Brightest Cluster Galaxies
MNRAS
432
530 -
DOI : 10.1093/mnras/stt490
Tran H. D.
2003
The Unified Model and Evolution of Active Galaxies: Implications from a Spectropolarimetric Study
ApJ
583
632 -
DOI : 10.1086/345473
Urry M. C.
,
Padovani P.
1995
Unified Schemes for Radio-Loud Active Galactic Nuclei
PASP
107
803 -
DOI : 10.1086/133630
Villarroel B.
,
Korn A.
2014
The Different Neighbours around Type-1 and Type-2 Active Galactic Nuclei
Nat. Phys.
10
417 -
DOI : 10.1038/nphys2951