A&A 470, 531-538 (2007)
DOI: 10.1051/0004-6361:20065266

The influence of AGN nuclear parameters on the FRI/FRII dichotomy

M. Wold1,2,3 - M. Lacy2 - L. Armus2


1 - European Southern Observatory, Karl-Schwarzschild str. 2, 85748 Garching bei München, Germany
2 - Spitzer Science Center, California Institute of Technology, MS 220-6, Pasadena, CA 91125, USA
3 - Present address: Institute of Theoretical Astrophysics, PO Box 1029 Blindern, 0315 Oslo, Norway

Received 24 March 2006 / Accepted 4 May 2007

Abstract
We have investigated the influence of nuclear parameters such as black hole mass and photoionizing luminosity on the FRI/FRII transition in a sample of nearby (z<0.2) radio galaxies from the 3CR catalogue. The sample was observed with medium-resolution, optical spectroscopy and contains some galaxies with unpublished velocity dispersion measurements and emission-line fluxes. The measured velocity dispersions for the sample lie in the range 130-340 km s-1 with a mean of 216 km s-1. Using the M-$\sigma_*$ relation, we convert to black hole mass and find that the black hole mass distribution is identical for FRI and FRII galaxies, with a mean of $\approx$2.5 $\times$ $10^{8}~M_{\odot}$. We determine narrow emission-line luminosities from [O II] and [O III] in our spectra, as well as from the literature, and convert them to photoionizing luminosities under the assumption that the gas is ionized by the nuclear UV continuum. Most of the galaxies with FRI morphology and/or low-excitation emission-line spectra have progressively lower black hole masses at lower photoionizing (and jet) luminosities. This agrees with the well-known Ledlow-Owen relation which states that the radio luminosity at the FRI/FRII transition depends on the optical luminosity of the host, $L_{\rm radio}\propto L_{\rm optical}^{1.8}$, because these two luminosities relate to AGN nuclear parameters. When recasting the Ledlow-Owen relation into black hole mass versus photoionizing luminosity and jet luminosity, we find that the recasted relation describes the sample quite well. Furthermore, the FRI/FRII transition is seen to occur at approximately an order of magnitude lower luminosity relative to the Eddington luminosity than the soft-to-hard transition in X-ray binaries. This difference is also consistent with the Ledlow-Owen relation, which predicts a weak black hole mass dependence in the transition luminosity in Eddington units. We conclude that the FRI/FRII dichotomy is caused by a combination of external and nuclear factors, with the latter dominating.

Key words: galaxies: active - galaxies: nuclei - galaxies: jets - X-rays: binaries

1 Introduction

Radio galaxies are usually classifed as FRI or FRII sources depending on their radio morphology. FRIs have smooth jets emanating from the nucleus and lobes where the surface brightness decreases toward the edges. The more powerful FRIIs have faint jets and edge-brightened lobes. The classification scheme was invented by Fanaroff & Riley (1974), who discovered that the FRI/FRII transition occurs at a radio luminosity of $P_{\rm 178~MHz} \approx 10^{25}$ W Hz-1 sr-1, with almost all sources below the transition luminosity being FRIs.

Over the years, much work has been put into understanding the remarkable FRI/FRII transition. Models fall into two different groups, those explaining the morphological differences as arising because of different physical conditions in the environment in which the radio source propagates, and those seeking to explain the dichotomy as caused by fundamental AGN parameters or the jets (see Gopal-Krishna & Wiita 2000, for a summary). The deceleration models, belonging to the former group, have been successful in reproducing the appearance of both FRI and FRII jets. In this scenario, the jets are thought to start out supersonically and slow down to a tran/subsonic flow because of entrainment of plasma in the host galaxy (Laing 1994; Kaiser & Alexander 1997; Bicknell 1995; De Young 1993; Laing 1996). The other group of models explains the differences as arising in more fundamental parameters like black hole (BH) spin, accretion mode, or jet composition (Meier 2001,1999; Ghisellini & Celotti 2001; Baum et al. 1995; Reynolds et al. 1996; Rees et al. 1982; Marchesini et al. 2004).

An interesting discovery was made by Ledlow & Owen (1996) (see also Owen & Laing 1989; Owen & White 1991; Owen & Ledlow 1994) who found that the FRI/FRII transition radio luminosity is an increasing function of the optical luminosity of the host galaxy, $L_{\rm rad} \propto L_{\rm opt}^{2}$. This means that the more optically luminous a galaxy is, the more powerful its radio source must be in order to produce FRII morphology. The Ledlow-Owen relation has been an important observation in the effort to explain radio galaxy dichotomy. However, both groups of models are able to explain, or reproduce, the Ledlow-Owen relation (Meier 1999; Gopal-Krishna & Wiita 2001; Bicknell 1995).

From the correlation between bulge luminosity and BH mass (Magorrian et al. 1998; Kormendy & Richstone 1995), we expect the host optical luminosity to scale with BH mass. Furthermore, because the radio-optical correlation (optical here means nuclear optical luminosity) for radio galaxies (Grimes et al. 2004; Willott et al. 1999; Saunders et al. 1989; Rawlings & Saunders 1991) relates radio luminosity to narrow-line luminosity, the radio luminosity can be expressed as a function of nuclear photoionizing luminosity. The Ledlow-Owen relation can thus be cast into the variables BH mass and nuclear photoionizing luminosity. Ghisellini & Celotti (2001) use this line of reasoning to argue that the Ledlow-Owen relation reflects a change in the accretion mode between the two classes because the transition occurs at a fixed ratio between BH mass and photoionizing/accretion luminosity. Recent work by Marchesini et al. (2004) also argues for a change in the accretion mode. They study a sample of radio-loud quasars and radio galaxies, finding that whereas quasars and radio galaxies have a similar distribution in BH mass, the distribution of accretion rates is bimodal with FRIs and weak-lined FRIIs on one side, and FRIIs and radio-loud quasars on the other.

In the above-mentioned work, optical luminosity is used as a measure of photoionizing luminosity. It is however still a matter of controversy whether this is applicable to FRI galaxies because it is uncertain whether they have dusty torii that can obscure the nucleus, leading to an underestimate of their true optical luminosity. One way to avoid the ambiguity in the interpretation of optical (and X-ray) nuclear luminosity is to use narrow emission lines as these are believed to arise from the narrow-line region outside the torus and be photoionized by the UV continuum from the central AGN.

In this paper we revisit the Ledlow-Owen relation for a sample of nearby FRI and FRII galaxies by utilizing stellar velocity dispersions to estimate BH masses and narrow emission-line luminosities to estimate photoionizing luminosities. We assume a cosmology with H0=70 km s-1 Mpc-1, $\Omega_{\rm m}=0.3$ and $\Omega_{\Lambda}=0.7$ unless otherwise noted. The abbreviations [O II] and [O III] are used for [O II]3727 Å and [O III]5007 throughout the paper.

2 The sample

The sample consists of sources drawn from the complete 3CRR catalogue of Laing et al. (1983). The selection criterion was $P_{178} \leq 10^{26.5}$ W Hz-1 sr-1, hence covering the FRI/FRII transition regime. Further constraints on RA and Dec gave a sample of 21 targets, ten FRIs and eleven FRIIs. Because of the redshift-luminosity correlation in the 3CRR sample, all targets lie at z < 0.2. The properties of the selected sources are listed in Table 1. Radio power at rest-frame 178 MHz was calculated from the 178 MHz flux and the spectral index as listed in the 3CRR Atlas web pages by Leahy, Bridle and Strom[*]. Radio morphology class was taken from the same web pages. In addition to morphology, we also list emission-line class as either high-excitation ("HEG'') or low-excitation ("LEG'') emission line galaxy. FRIIs are classified as HEGs if [O III] is strong compared to [O II], and as LEGs if [O III] is weak or absent (Laing et al. 1994). FRIs have low-excitation spectra, hence all the FRIs are listed as LEGs in Table 1.

Emission-line classes for all objects were found in the literature (see table note) except for 3C 442A, which we classify as a LEG because [O III] is weaker than [O II] in our spectrum. This galaxy is not a classical FRII with well-defined hot spots (one of the lobes has FRI morphology), and the host galaxy is part of an ongoing merger (Comins & Owen 1991). The other source with non-classical radio morphology in the sample is 3C 433, highly asymmetric and probably interacting with close neighbours. The radio morphology is classified as FRI by Leahy, Bridle and Strom, but we have chosen to define it as an FRII because its optical spectrum is characteristic of an HEG and there is evidence that it harbours an obscured quasar nucleus (Fabbiano et al. 1984; Yates & Longair 1989; Lilly et al. 1985).


  \begin{figure}
\par\includegraphics[width=5.5cm,clip]{5266fig1.eps}\end{figure} Figure 1: Open circles correspond to FRIs and stars to FRIIs, with filled symbols for high-excitation emission line galaxies. The solid line marks the Ledlow-Owen relation.
Open with DEXTER

In order to compare the sources with the Ledlow-Owen relation, we plot the sample in the radio-optical plane in Fig. 1. The apparent rest-frame $R_{\rm c}$ magnitudes were converted to absolute magnitudes, and radio luminosity at 1.4 GHz calculated using the observed flux at 1.4 GHz and spectral indices as listed in Table 1. For calculating luminosities, we used the same cosmology as Ledlow & Owen (1996). The figure shows that the galaxies fall in the expected regions of the diagram, with FRI sources below and FRII sources above the Ledlow-Owen relation, but there are also some borderline cases. All galaxies classified as HEGs are well above the relation, consistent with the unified schemes (Chiaberge et al. 2002; Jackson & Wall 1999).

Table 1: The sample. $P_{\rm 178}$ in Col. (3) is total radio power (W Hz-1sr -1) at rest-frame 178 MHz. Column (4) lists spectral index between 178 and 750 MHz, and Cols. (5) and (6) the radio morphology and emission-line type. Apparent magnitude in rest-frame $R_{\rm C}$-filter and corresponding references are listed in Cols. (7) and (8). Columns (9) and (10) show velocity dispersion measurements from the literature (km s-1) and corresponding references. The last three columns list emission line fluxes in units of 10-15 erg s-1 cm-2, both from the literature and as determined from our wide slit spectrophotometry, and references.

3 Observations and data reduction

Medium resolution ($R\sim3000$-4000) optical spectra were obtained with the Double Spectrograph (Oke & Gunn 1982) on the 5 m Hale Telescope at Palomar Observatory. The spectrograph was equipped with two 1k $\times$ 1k CCDs and the D48 dichroic, dividing the red and blue side at approximately 4800 Å. We utilized the 1200 l/mm grating blazed to 7100 Å on the red CCD and to 4700 Å on the blue CCD. Two different grating angles were used depending on the galaxy redshift such that the Mg IB $\lambda$5175 Å absorption line complex fell on the center of the red CCD, and the Ca H & K  $\lambda\lambda$3934, 3969 Å doublet at the center of the blue CCD. The red and the blue CCDs had pixel scales of 0 $.\!\!^{\prime\prime}$468 and 0 $.\!\!^{\prime\prime}$624, resulting in dispersions of 0.624 and 0.864 Å/pix, respectively.

Absorption line spectra were obtained with a 1 $^{\prime\prime}$ wide slit, giving spectral resolutions of 1.33 and 1.39 Å, corresponding to velocity resolutions of $\approx$70 and $\approx$100 km s-1 in the red and the blue, respectively. In order to measure total narrow emission line fluxes, we also took spectra with a wide slit (4 $^{\prime\prime}$) of part of the sample. Slit lengths were 128 $^{\prime\prime}$. Spectra were taken by centering the slit on the optical nucleus of the galaxy and aligning it with the major axis in order to maximize the amount of light entering the slit. Each target was also observed as close to the parallactic angle as possible in order to minimize losses due to differential atmospheric refraction. Details about the observations are given in Table 2.

Wavelength calibration was accomplished by taking arc lamp spectra at every telescope pointing and standard stars were observed at the beginning and end of each night. There were some clouds at the beginning of the first night, but for the rest of the run the weather was clear and photometric, and the seeing stable at $\approx$1 $^{\prime\prime}$. The photometric zero point varied 0.1-0.2 mag between the beginning and end of the second and third night. For the first night, which started out as non-photometric, the difference was $\approx$0.7 mag. For spectrophotometry with the 4 $^{\prime\prime}$ slit, we therefore waited until the middle of the first night when the conditions were more favourable. As templates for the velocity dispersion fitting we obtained spectra of one sub giant G0 star (HD 188121) and one giant K1 star (HD 28191). In order to avoid seeing-limited spectral resolution with the 1 $^{\prime\prime}$ slit, the template star was moved across the slit during integration.

Table 2: Table of observations. For each galaxy, the slit position angle and exposure time is shown.

The spectra were reduced and extracted in a standard manner using IRAF[*] tasks in the twodspec and onedspec packages. Night sky lines from O I and from mercury in street-lamps were used to check the wavelength calibration, and an agreement of typically $\la$0.5 pixels was found, corresponding to $\la$18 and $\la$32 km s-1 in the red and blue, respectively.

Two of the spectra, 3C 83.1B and 3C 386, are dominated by light from foreground stars. For this reason we do not consider these two galaxies in our analysis. In the case of 3C 83.1B, the star lies 3 $^{\prime\prime}$ to the east of the galaxy nucleus (de Koff et al. 2000; Poulain et al. 1992), and for 3C 386 the star is superposed on the optical galaxy core (Chiaberge et al. 2002).

4 Template fitting procedure

We fit stellar templates to the galaxy spectra in order to measure velocity dispersions using the direct fitting method described by Barth et al. (2002). Galaxy and stellar template spectra were extracted using aperture diameters of  $1\hbox{$.\!\!^{\prime\prime}$ }872$ in order to cover the same physical scale in both spectral arms. A model spectrum, M(x), where $x=\ln \lambda$ is measured in the galaxy rest frame was evaluated as

\begin{displaymath}%
M(x) = \left\{ \left[T(x) \otimes G(x)\right] + C(x) \right\} \times P(x),
\end{displaymath} (1)

where $[T(x) \otimes G(x)]$ denotes the template spectrum convolved with a Gaussian of width  $\sigma_*$, C(x) is the continuum and P(x) a polynomial. The best-fit model was found by minimizing $\chi^{2}$ using a downhill simplex method implemented in the amoeba algorithm (Press & Spergel 1988).

We experimented with different functional forms of C(x) and P(x), starting with a straight line for the continuum, and a Legendre polynomial for P(x). However, since the change in x is small over the fitting range, coefficient degeneracies were easily introduced, hence we rescaled the fitting range to a new variable  $\theta
\in[0,\pi]$ and used Fourier-type expansions instead. The continuum is then expressed as

\begin{displaymath}%
C(\theta)=c_{0}+c_{1}\sin\theta,
\end{displaymath} (2)

and the polynomial as

\begin{displaymath}%
P(\theta) = p_{0} + p_{1}\cos\theta + p_{2}\cos2\theta +
p_{3}\cos3\theta + p_{4}\cos4\theta.
\end{displaymath} (3)

The sine and cosine terms are orthogonal, thereby breaking part of the degeneracy.

Fits to three different galaxies along with the stellar template used to construct the model are shown in Fig. 2. Generally, we found that the regions redward of Ca H&K and the Mg Ib absorption line complex are best suited for fitting velocity dispersions. This has also been experienced in other works, e.g. Barth et al. (2002) and Greene & Ho (2006). We therefore tried to avoid the Mg Ib region, except for in some cases where the Fe blends redward of the region were not sufficiently covered. As Barth et al., we also find excess emission at $\approx$5200 Å in several cases (caused by [N I]5199 Å emission in the host), complicating the fitting. For the red arm spectra, we therefore started several of the fits at $\approx$5220 Å. The results from the fitting and the fitting regions used are listed in Table 3.

The red arm spectra were easier to fit than those from the blue arm, probably because the red spectra have better signal-to-noise and containing less spectral features to fit. Also, because the galaxies contain spectral features that match better with K stars than with G stars in the chosen fitting regions, the K1 template was found to produce better fits with less scatter than the G0 template. As there is also less contamination by the AGN continuum in the red part, we chose to use the red arm spectra with the K1 template. In Fig. 3 we show a comparison between measurements from the red and the blue arm in cases where spectra from both could be fitted. The errors are seen to be larger for the blue fits, and the two agree to within typically ten per cent, except for two cases, 3C 449 and 3C 98 where there are differences at the 20-30% level.


  \begin{figure}
\par\includegraphics[width=16cm,clip]{5266fig2.eps}\end{figure} Figure 2: Examples of radio galaxy spectra with best-fit models overplotted as solid lines where an acceptable fit could be obtained. The K1 star used as a template for the fitting is shown at the top of each panel.
Open with DEXTER

Table 3: Results from the fitting. The fitting regions and signal-to-noise per pixel in the continuum is listed. Best-fit velocity dispersions normalized to the same physical aperture ( $r_{\rm ap} = 0.595~h^{-1}$ kpc) and their confidence intervals are shown. "AGN cont.'' indicates that there is too much AGN contamination in the spectrum to obtain a reliable fit.

We tested the fitting routine by constructing artificial galaxy spectra from the two stellar templates (50% K1- and 50% G0-star), broadening them by a specified amount and applying the fitting routine. Noise spectra were also added, broadened to the instrumental resolution ($\approx$45 km s-1) and scaled in order to obtain a signal-to-noise ratio of approximately 15. The result from the simulations is shown Fig. 4, where it can be seen that the velocity dispersions are generally recovered to within 10% of the input value. As expected, at dispersions lower than the instrumental resolution there is larger scatter.

Nine of the radio galaxies in the sample have velocity dispersions from the literature, as listed in column nine of Table 1. A comparison with the literature measurements is shown in Fig. 5. There is agreement to within less than 20%, except for one galaxy, 3C 388, where our estimated velocity dispersion is 36% lower than the literature value. This galaxy is one of the objects in our sample with noisier spectra, hence it has larger error bars on the velocity dispersion than the others. Part of the discrepancy could be due to our estimate and the estimate by Heckman et al. (1985) being based on different spectral regions. The blue part of the spectrum, used by Heckman et al. (1985), could be more contaminated by AGN continuum than the red part, affecting velocity dispersion estimates.

5 Discussion

   
5.1 Conversion to AGN parameters

We convert the velocity dispersions derived from the spectral fitting to a common physical scale following the scheme by Jørgensen et al. (1995), and thereafter to BH mass using the $M_{\rm
BH}$- $\sigma_*$ relation from Tremaine et al. (2002). We find that the sample spans a relatively narrow range in BH mass, from 2.2 $\times$ 107 to 1.1 $\times$ 10 $^{9}~M_{\odot}$. Assuming that we are able to recover the velocity dispersions to within an accuracy of $\approx$10% as shown by the simulations, and including a scatter of 0.3 dex in $\log M_{\rm BH}$ (Tremaine et al. 2002) around the $M_{\rm BH}-\sigma_{*}$ relation, we estimate that we are able to determine BH masses to an accuracy of approximately 0.35 dex.

Narrow emission line fluxes for [O II] and [O III] were determined from our 4 $^{\prime\prime}$ slit spectra, and some were also taken from the literature, see Cols. 11 and 12 of Table 1. Under the assumption that the narrow emission lines are photoionzied by the UV continuum from the AGN, we convert to total luminosity in the narrow-line region using the relation by Rawlings & Saunders (1991), $L_{\rm NLR}
\approx 3\left(3L_{\rm [OII]}+1.5 L_{\rm [OIII]}\right)$. In cases where only one of the two line fluxes are known (for FRIIs), we apply a relation between $L_{\rm [OII]}$ and $L_{\rm [OIII]}$ from Grimes et al. (2004), whereas for FRIs we take $L_{\rm [OII]} = L_{\rm
[OIII]}$. If there is an upper limit on the flux, we take the line flux to be equal to the limit. The narrow line flux is converted to photoionizing luminosity, $Q_{\rm phot}$, by assuming a covering factor of $\kappa=0.005$ (Willott et al. 1999) for the narrow-line gas. We estimate that the uncertainty in $\log Q_{\rm phot}$ is rougly $\pm $1 dex, which includes uncertainties related to the covering factor for the narrow-line gas and the fact that narrow lines may also be powered by shocks (Inskip et al. 2002; Dopita & Sutherland 1995).

In order to estimate the average kinetic power transported by the jets to the lobes, we use the relation of Willott et al. (1999), $Q_{\rm jet} =
3$ $\times$ $10^{38} L_{\rm 151}^{6/7}$ Watts, where  $L_{\rm 151}$ has units of 1028 W Hz-1 sr-1. The radio luminosity at 151 MHz for the sources in our sample was taken from Grime's compilation[*]. Given that Willott et al.'s relation applies to a "typical'' radio galaxy of median age, $\pm $1 dex is a reasonable estimate for the uncertainty in  $\log Q_{\rm jet}$.


  \begin{figure}
\par\includegraphics[width=5.8cm,clip]{5266fig3.eps}\end{figure} Figure 3: Comparison of red and blue velocity dispersions. Error bars represent 95% confidence intervals.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=5.55cm,clip]{5266fig4.eps}\end{figure} Figure 4: Results from simulations with artificial galaxy spectra. The input velocity dispersions, $\sigma _{\rm in}$, are recovered ( $\sigma _{\rm out}$) to within 10%.
Open with DEXTER

We also derive the dimensionless ratio between the jet luminosity and the photoionizing luminosity, $\log (Q_{\rm jet}/Q_{\rm phot})$, and estimate an uncertainty in this parameter of approximately $\pm $1.5 dex. The $Q_{\rm jet}/Q_{\rm phot}$ ratio can be seen as the ratio of kinetic to radiative energy output by the AGN, and indicates how effectively jets are formed. A small ratio implies that the energy output is dominated by thermal radiation from an accretion disk, whereas a large ratio indicates that most of the energy goes into forming jets, i.e. to kinetic energy, and that the accretion process thus may be radiatively inefficient.

   
5.2 The Ledlow-Owen relation

The parameters $Q_{\rm jet}$, $Q_{\rm phot}$ and their ratio are plotted as a function of BH mass in Fig. 6. The line of separation between FRIs and FRIIs, $L_{\rm radio}\propto L_{\rm optical}^{1.8}$ (Ledlow & Owen 1996), can be written as

 \begin{displaymath}%
\log\left(\frac{P_{\rm 1.4}}{{\rm W~Hz^{-1}}}\right) = -0.67 M_{R} + 10.13,
\end{displaymath} (4)

where MR is the optical absolute magnitude of the host and  $P_{\rm 1.4}$ is the radio power at 1.4 GHz. By using the relation by McLure & Dunlop (2001) which relates MR to  $M_{\rm
BH}$ and the relation from Willott et al. (1999) linking $P_{\rm 1.4}$ to  $Q_{\rm jet}$, we obtain

 \begin{displaymath}%
\left(\frac{Q_{\rm jet}}{{\rm W}}\right) = 6.17\times10^{27}\left(\frac{M_{\rm BH}}{M_{\odot}}\right)^{0.92}\cdot
\end{displaymath} (5)

This translates into $Q_{\rm jet} \simeq 5$ $\times$ $10^{-4} L_{\rm Edd}$, implying that the FRI/FRII transition occurs at a ratio of $\simeq$0.0005 between jet power and Eddington luminosity if $M_{\rm
BH}\approx10^{8}~M_{\odot}$. Ghisellini & Celotti (2001) derive a ratio of 0.015, but this is probably erroneous[*]. The left-hand panel of Fig. 6 shows $Q_{\rm jet}$ as a function of  $M_{\rm
BH}$, and the Ledlow-Owen relation in the form $Q_{\rm jet} = 5$ $\times$ $10^{-4} M_{\rm BH}$ is indicated with a solid line and seems to fit the data well.


  \begin{figure}
\par\includegraphics[width=5.55cm,clip]{5266fig5.eps}\end{figure} Figure 5: Comparison with dispersions from the literature.
Open with DEXTER

The Ledlow-Owen relation can also be written in terms of photoionizing luminosity through the radio-optical correlation for radio galaxies (Willott et al. 1999). Doing this, we find $Q_{\rm phot} \simeq 0.0037
L_{\rm Edd}$ for the typical BH mass in our sample. (Ghisellini & Celotti 2001 derive $Q_{\rm phot} \simeq 0.005 L_{\rm Edd}$). The solid line in the middle panel of Fig. 6 shows the Ledlow-Owen relation recast in the two variables  $Q_{\rm phot}$ and  $M_{\rm
BH}$and can be seen to fit the transition region quite well. However, the separation of FRIs and FRIIs appears to be cleaner in the diagram of $\log Q_{\rm phot}$ vs. $\log M_{\rm BH}$, probably because $Q_{\rm jet}$ is more influenced by external factors, such as conditions in the interstellar medium of the host. If the accretion efficiency is approximately 10%, and assumed to be in the form of a standard optically thick, geometrically thin accretion disk (Shakura & Sunyaev 1973), the relation $Q_{\rm phot} \approx 0.004 L_{\rm Edd}$indicates that the FRI/FRII transition occurs at $\sim$4% of the Eddington rate.

The transition can also be described as occurring at an almost fixed ratio between $Q_{\rm jet}$ and  $Q_{\rm phot}$, i.e. at a fixed "radio loudness''. Using the relations above, we find $Q_{\rm
jet}/Q_{\rm phot} \approx 0.1$, shown as a solid line in the rightmost panel of Fig. 6, but the sources do not separate as well in this diagram as in the other two.


  \begin{figure}
\par\includegraphics[width=15.9cm,clip]{5266fig6.eps}\end{figure} Figure 6: From left to right, jet power, photoionizing luminosity and their ratio as a function of BH mass. An uncertainty of $\pm $1 dex in $Q_{\rm jet}$ and $Q_{\rm phot}$ has been assumed, and the 95% confidence interval on BH mass is used. Symbols are as in Fig. 1.
Open with DEXTER

Both AGN and X-ray binaries (XRBs) are associated with accretion onto compact objects and the production of jets. The two have therefore been compared at several occations (e.g. McHardy et al. 2006; Merloni et al. 2003; Maccarone 2003; Falcke et al. 2004). XRBs are characterized by three different spectral states (see e.g. Gallo et al. 2003). The low/hard state has a hard power-law spectrum, a weak thermal component from an accretion disk/flow, and steady jet emission is seen at radio wavelengths. FRIs have been suggested as the high BH mass analogue of the low/hard state in XRBs (see e.g. Merloni et al. 2003; Falcke et al. 2004) and having accretion flows with low accretion rates and inefficient cooling, such as advection-dominated accretion flows (Narayan & Yi 1994; Rees et al. 1982). FRIIs, on the other hand, are thought of as analogous to the high/soft and very high state XRBs with accretion disks modeled as standard, optically thick, geometrically thin disks (Shakura & Sunyaev 1973). XRBs in the high/soft state show no strong or steady radio emission, and jet formation is thought to be inhibited. In the so-called very high state, thought to occur at high accretion rates, XRBs show transient jet phenomena.

If, when we derive the relations in Eqs. (4) and (5), keep the BH mass dependence instead of substituting the mean black hole mass of the sample, we find that $Q_{\rm phot}/L_{\rm
Edd} \propto (M_{\rm BH}/M_{\odot})^{-0.14}$ and $Q_{\rm jet}/L_{\rm
Edd} \propto (M_{\rm BH}/M_{\odot})^{-0.08}$. Hence both the photoionizing and jet luminosity in Eddington units at the FR/FRII transition is weakly dependent on BH mass. It is too weak to be noticeable in a sample of radio galaxies, but when extrapolated to XRBs there should be roughly a factor of 10 difference. From the Ledlow-Owen relation we therefore expect that the FRI/FRII transition luminosity in Eddington units is about a factor of 10 lower than for spectral state transition in XRBs. This does indeed seem to be observed. For the FRI/FRII transition we observe transition luminosities of $\approx$0.4% of the Eddington luminosity, see Fig. 6. The transition luminosity of XRBs from the high/soft to the low/hard state is $\approx$2% of the Eddington luminosity (may be up to four times higher for the low/hard to high/soft transition) (Maccarone et al. 2003). The concurrence of observed transition luminosities with those predicted from radio galaxy scaling relations may indicate that the transition between low/hard and high/soft states in XRBs is indeed similar in nature to the FRI/FRII transition. Maccarone (2003) discusses that neutron star XRBs may have higher transition luminosities than BH XRBs, also indicating that the transition luminosity may depend on the mass of the compact object.

If FRIIs are analogues to the very high state XRBs with thermally unstable and radiation pressure dominated accretion disks (Meier 2001), the FRI/FRII transition could correspond to a switch between the low/hard and the very high state, without going through the high/soft state. Meier (2001) offers a solution which has the transition luminosity relative to the Eddington luminosity proportional to $M_{\rm BH}^{-1/8}$ (see also Merloni 2003; Chen et al. 1995). There is some inconsistency in the implied transition luminosities in this scenario because the very high state in XRBs is thought to occur at a (bolometric) luminosity of 20-30% of the Eddington luminosity (Meier 2001; Maccarone et al. 2003), whereas we find <1% for the radio galaxies. The discrepancy could be due to our different methods of estimating the bolometric/photoionizing luminosity, but another possibility is that the physical conditions (e.g. density) in the accretion flow close to the BH are different for radio galaxies and XRBs.

Sources that do not follow the Ledlow-Owen relation may be the ones in clusters or high density environments (Marchesini et al. 2004). Besides external density in this case probably being a dominant factor for the radio morphology, another possible explanation is that mergers or collisions in cluster galaxies may have altered the spin of the BH (Natarajan & Pringle 1998) and BH spin is thought to be important for the formation of jets.

The Ledlow-Owen relation nevertheless seems to be an indication that the FRI/FRII morphology is largely determined by nuclear AGN parameters. In hindsight it is perhaps not surprising that the empirical scaling relations used by Bicknell (1995) to explain the Ledlow-Owen relation relate to the central stellar velocity dispersion.

6 Summary and conclusions

We have investigated BH masses and narrow emission line luminosities in a sample of FRI and FRII radio galaxies. BH masses were estimated via the $M_{\rm
BH}$-$\sigma_*$ relation by determining central stellar velocity dispersions from medium resolution spectra, and photoionizing luminosities were estimated from spectrophotometry of the narrow [O II]3727 Å and [O III]5007 Å emission lines. The measured velocity dispersions span the range 130-340 km s-1with a mean of $\approx$215 km s-1. Including the scatter in the $M_{\rm
BH}$-$\sigma_*$ relation, we estimate that we are able to determine BH masses with an accuracy of $\approx$0.35 dex. We find FRIs and FRIIs to have the similar BH mass distributions, confirming previous results (e.g. Marchesini et al. 2004; McLure et al. 2004).

The data fit the Ledlow-Owen relation well when it is rewritten in terms of jet luminosity, photoionizing luminosity and BH mass (Ghisellini & Celotti 2001). The separation of FRIs and FRIIs appears cleaner in the $Q_{\rm phot}$- $M_{\rm
BH}$ plane than in the $Q_{\rm jet}$- $M_{\rm
BH}$ plane, indicating that FRIs and FRIIs are more easily separated in terms of nuclear parameters than in terms of parameters such as $Q_{\rm jet}$ which may be more influenced by external factors. However, it is still possible that the FR/FRII transition may be caused by a mixture of both nuclear and external factors.

We discuss whether the FRI/FRII transition can be compared to spectral state transitions in XRBs, and highlight the fact that the Ledlow-Owen relation predicts a BH mass dependence in the transition luminosity in Eddington units. This implies that radio galaxies change from FRI to FRII at luminosities (in Eddington units) that are a factor of 10 lower that those of XRBs.

Acknowledgements
The authours thank the referee for comments which helped to improve the original manuscript. This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. The authors would also like to thank the maintainers of the DRAGN web page (http://www.jb.man.ac.uk/altas) which was used in this work.

References

 

Copyright ESO 2007