Front. Astron. Space Sci. Frontiers in Astronomy and Space Sciences Front. Astron. Space Sci. 2296-987X Frontiers Media S.A. 10.3389/fspas.2017.00022 Astronomy and Space Sciences Original Research Confrontation of the Magnetically Arrested Disc Scenario with Observations of FR II Sources Rusinek Katarzyna * Sikora Marek Nicolaus Copernicus Astronomical Center, Polish Academy of Sciences, Warsaw, Poland

Edited by: Paola Marziani, Osservatorio Astronomico di Padova (INAF), Italy

Reviewed by: Milan S. Dimitrijevic, Belgrade Astronomical Observatory, Serbia; Brian Punsly, International Center for Relativistic Astrophysics, Italy

*Correspondence: Katarzyna Rusinek krusinek@camk.edu.pl

This article was submitted to Cosmology, a section of the journal Frontiers in Astronomy and Space Sciences

 12 10 2017 2017 4 22 29 08 2017 27 09 2017 Copyright © 2017 Rusinek and Sikora. 2017 Rusinek and Sikora

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

The main aim of our work was to check whether powers of jets in FR II radio galaxies (RGs) and quasars (QSOs) can be reproduced by the Magnetically Arrested Disc (MAD) scenario. Assuming that established in the recent numerical simulations of the MAD scenario the (H/R)2 dependence of the jet production efficiency is correct, we demonstrate that in order to reproduce the observed jet powers in FR II sources: (i) accretion discs must be geometrically much thicker than the standard ones; (ii) and/or that the jet production is strongly modulated.

 active galactic nuclei (AGN) quasars relativistic jets accretion discs magnetically arrested disc (MAD) 2013/09/B/ST9/00026 2016/21/B/ST9/01620 RADIOLIFE-320745 Narodowe Centrum Nauki10.13039/501100004281 FP7 Ideas: European Research Council10.13039/100011199 1. Introduction

The primary division of Active Galactic Nuclei (AGNs) is based on their radio loudness which is defined as the ratio of radio luminosity (typically 5 GHz) to the optical luminosity (typically in B band). Minority of AGNs are known to be radio-loud (10% in case of QSOs, Kellermann et al., 1989) and their radio loudness is sometimes even up to 3–4 orders of magnitude higher than that of radio-quiet AGNs. Moreover, which was found in last years (Rawlings and Saunders, 1991; Punsly, 2007; Fernandes et al., 2011; Sikora et al., 2013), jet powers Pj of many radio galaxies reach values comparable to the accretion powers Ṁc2, where Ṁ is the mass accretion rate. In order to produce jets with such high efficiencies, ηj≡Pj/(Ṁc2)≃1, by the Blandford and Znajek (1977) mechanism, the amount of the magnetic flux required to be accumulated on the black hole (BH) is so large that it can only be maintained if it is confined by the ram pressure of the accreting plasma (Sikora and Begelman, 2013). This affects the accretion flow in such a way that the innermost portion of the accretion flow is dynamically dominated by the poloidal magnetic field causing that accretion proceeds via interchange instabilities. This scenario is called Magnetically Arrested Disc (MAD, Narayan et al., 2003; Igumenshchev, 2008; Punsly et al., 2009; Tchekhovskoy et al., 2011; McKinney et al., 2012).

Even though the MAD scenario appears to be an attractive and plausible way to explain the existence of the most powerful jets in radio-loud FR II AGNs (which as described by Fanaroff and Riley, 1974, are characterized by edge-brightened radio structures in contrast to centre-brightened RGs of class I), there are still some open problems. One of them is to establish which parameters truly decide on the MAD occurrence. Assuming that the jet production efficiency depends mostly on the BH spin, van Velzen and Falcke (2013) showed that in fact it is not a dominant parameter for powering jets. Furthermore, since the observed mean efficiency of their FR II quasars sample is much lower than maximal predicted by the MAD model, they concluded that this mechanism does not occur in these sources. On the other hand, in their recent studies, Avara et al. (2016) performed numerical simulations of thin MADs in order to investigate its efficiencies. They found that not only BH spin a but also a geometrical thickness H/R (where H is the disc height and R is the distance from the black hole) contributes to the jet power efficiency as ηj~a2(H/R)2. Their results confirm that the MAD scenario can be responsible for explaining powerful jet systems but also point out that at moderate accretion rates discs must be geometrically thicker than the standard theory predicts.

Given the above we checked whether jet powers in FR II radio galaxies and quasars can be reproduced by the MAD scenario. Our studies were thoroughly described in Rusinek et al. (2017). In this paper we briefly outline samples we used and discuss the origins of results we obtained, especially focusing on the importance of geometrical thickness of accretion discs and modulation of jet production.

 2. Samples and analysis

In order to adequately assess the distribution of radio galaxies and quasars in the Pj/Ld – λEdd plane (where Ld is the bolometric disc luminosity and λEdd is the Eddington ratio) we have combined four different samples of these sources. We used radio and optical data to calculate necessary values. Including both low- as well as high-redshift objects allowed us to check if cosmological evolution of jet production may have an impact on our final results.

2.1. Samples

FR II Narrow-Line Radio Galaxies (NLRGs) were extracted from the sample of z < 0.4 radio galaxies with extended radio structure selected by Sikora et al. (2013). The objects are taken from Cambridge catalogues and matched with the SDSS (Sloan Digital Sky Survey), FIRST (Faint Images of the Radio Sky at 20 cm) and NVSS (National Radio Astronomy Observatory (NRAO) Very Large Array (VLA) Sky Survey) catalogues. Due to the available optical data, this sample contains 152 sources;

FR II quasar sample was obtained by van Velzen et al. (2015) on the selection of double-lobed radio sources from the FIRST survey catalog, and cross-matching with SDSS quasars. The BH masses and Eddington ratios, when available, were taken from Shen et al. (2011) thereby reducing the sample from 458 to 414 objects;

FR II NLRGs selected in 0.9 < z < 1.1 were taken from Fernandes et al. (2011, 2015). The main reason of adding these sources was to verify how much the incompletness of very massive BHs in the local Universe affects the jet production efficiency FR II NLRGs (the first sample). This sample contains 27 objects;

Low-redshift Broad-Line Radio Galaxies (BLRGs) and radio-loud quasars (RLQs) were used by Sikora et al. (2007) to study radio loudness of these objects. We decided to add this sample (86 sources) to check if the incompletness of SDSS quasars at moderate accretion rates (the second sample) can significantly influence the average value of Pj/Ld of the FR II quasar sample.

 2.2. Jet production efficiency

For each of mentioned in the previous paragraph samples we calculated the following properties:

– The jet power Pj was obtained from the different radio luminosities (from the range of 151 MHz to 5 GHz) according to the formula based on calorimetry of radio lobes which was proposed by Willott et al. (1999) as Pj[ergs-1]=5.0×1022(f/10)3/2(L1.4[WHz-1])6/7, where f is the parameter accounting for errors in the model assumptions. Its value (typically between 10 and 20, Blundell and Rawlings, 2000) is established based on comparing jet powers of luminous FR II sources calculated from the model of hotspots Godfrey and Shabala (2013) and the one provided by Willott et al. (1999), and is adpoted as f = 10;

– The bolometric disc luminosity Ld which is related to the line or optical/IR luminosities via respective bolometric corrections. It is used as a proxy of the accretion power Ṁc2=Ld/ϵd where ϵd is the disc radiation efficiency depending on BH spin. In standard disc theory it is assumed to be ϵd = 0.1;

– The Eddington luminosity LEdd which is necessary to establish the Eddington ratio λEdd ≡ Ld/LEdd describing the accretion rate, LEdd ∝ MBH where MBH is the BH mass.

We have adopted the Λ cold dark matter (ΛCDM) cosmology with H0=70km s-1, Ωm = 0.3 and ΩΛ = 0.7. Uncertainties of Pj/Ld and λEdd, calculated as standard deviations of ratios of two independently determined quantities, are not bigger than 0.4 dex for each of the samples. All details about samples and methods used to estimate our calculations are described in Rusinek et al. (2017).

 2.3. Results

The dependence of Pj/Ld on the λEdd for all four samples is presented in Figure 1. The distribution of all objects confirms the trend of the increase of Pj/Ld with decreasing Eddington ratio (Sikora et al., 2007). All samples are pretty consistent with each other: z < 0.4 FR II NLRGs and 0.9 < z < 1.1 NLRGs samples have similar median of Pj/Ld at moderate accretion rates (which corresponds to λEdd > 0.003); FR II quasars and BLRGs+RLQs samples overlap each other.

The dependence of Pj/Ld on the λEdd for all samples: z < 0.4 FR II NLRGs as green crosses; FR II quasars as grey dots; 0.9 < z < 1.1 NLRGs as empty circles; BLRGs+RLQs as black triangles. The horizontal dashed line corresponds to Pj=(ϵd/0.1)Ṁc2. The inclined dashed line is derived from the Equation (2) for a = 0.5 and assuming a correct normalization. The presence of objects in the shaded region of the plot calls into question our underlying assumptions. Purple arrows show how the modulation of jet production affects the Pj/Ld − λEdd distribution.

The horizontal dashed line marks the upper limits for Pj/Ld predicted by the MAD model for the maximal BH spin and adopting ϵd = 0.1. The inclined dashed line results from the dependence of the jet production efficiency on H/R given as

ηj≃4a2(1+0.3a1+2(H/R)4)2(H/R)2

(Avara et al., 2016). For 0.03 < λEdd < 1.0 discs have H/R ≥ 1 and then such a dependence disappears. The existence of objects in the shaded area, which represents a significant fraction of all studied sources, contradicts with the above predictions and challenges current jet production theories.

 3. Discussion

Depending on the Eddington ratio few different explanations for obtained by us jet powers can be proposed. The transition from Radiatively Inefficient, optically thin Accretion Flows (RIAF) to the standard, optically thick accretion discs occurs at λEdd ≃ 0.01 (Best and Heckman, 2012). At lower Eddington ratios the disc radiative efficiency ϵd is expected to be much lower than usually assumed 0.1 and this can explain the presence of objects with Pj/Ld > 10 in this area. At higher Eddington ratios there are still some sources with Pj/Ld > 10. All of them belong to the 0.9 < z < 1.1 NLRGs sample for which the bolometric disc luminosity was calculated from the mid-IR data. As it was pointed out by Ogle et al. (2006) this method is not very accurate and therefore discs luminosities of sources taken from Fernandes et al. (2011) may be overestimated.

The presence of objects in the shaded area corresponding to moderate accretion ratios and above the sloping dashed line may indicate that the accretion discs are thicker than the standard ones or/and that the jet production is modulated. However it should be noted here that the extension of the shaded area towards Pj/Ld ≪ 1 on the Figure 1 may be inappropriate by noting that the dependence of the jet production efficiency on geometrical thickness of the disk given by Equation (2) can be overestimated resulting from the very approximate treatment of the radiative transfer in Avara et al. (2016) simulations and by not including contribution to the vertical pressure in the MAD zone from toroidal magnetic fields.

3.1. Thicker accretion discs

Assuming maximally rotating black hole with a = 1 and producing radiation at a rate λEdd ~ 0.01, the standard accretion disc model (Novikov and Thorne, 1973; Laor and Netzer, 1989) predicts that the maximal value of geometrical thickness of the disc is H/R ~ 0.4 which gives the jet production efficiency equal to 0.01 (from the Equation 2). This value is much lower than the median value of FR II quasars sample as well as its upper bound in the Pj/Ld − λEdd plot. This discrepancy can be explained if instead of optically thick, geometrically thin standard disc the thicker ones will be considered. These kind of discs have been proposed to avoid gravitational and thermal instabilities in gas and radiation pressure supported discs and they can be formed in presence of strong toroidal magnetic fields (e.g., Begelman and Pringle, 2007; Sądowski, 2016).

Different approach assumes the existence of moderately hot, optically thick and massive layer on top of the relatively colder, accretion disc which stays geometrically thin and optically thick as it was originally developed by Shakura and Sunyaev (1973). Różańska et al. (2015), and Begelman et al. (2015) proved that the model of heavly, viscously driven corona is real if the disc/corona system is stabilized by either strong magnetic fields or vertical outflows. The massive, dense coronas can also arise during the transition from Shakura and Sunyaev (1973) discs to Advection-Dominated Acretion Flows (ADAFs, Abramowicz et al., 1995; Narayan and Yi, 1995) which the transition coincides with the Luminous Hot Accretion Flow (LHAF, Yuan, 2003; Yuan and Narayan, 2014).

 3.2. Modulation of jet production

Another worth considering idea which may explain the visible Pj/Ld − λEdd distribution is taking into account the variability of accretion rate. The method we used to establish jet powers is based on the total energy content of the radio lobes. Their lifetime is long (even up to 108 years, Komissarov and Gubanov, 1994) so the only possibility of observing variations of jet powers is connected with the hotspot luminosity which may vary on much shorter time-scales. Hotspots however make a small contribution to the total radio luminosity (Mullin et al., 2008) as a result of which jet powers from integrated lobe luminosities are not significantly affected by the variability of the accretion rate. This property may be noticed in the disc luminosity though, which is a direct measure of instantaneous accretion rate. Taken together these findings implicate a modulation of both, the “apparent” jet production efficiency as well as the Eddington ratio in such a way that with decreasing λEdd the Pj/Ld increases. This causes the stretching with the slope −1 on the plot which is presented as the set of arrows on the Figure 1. just as the sources are distributed. Variability of accretion rate may be naturally caused by the viscous instabilities in accretion discs (Janiuk et al., 2002; Janiuk and Czerny, 2011).

 4. Summary

In this proceeding we argue that jets in FR II radio galaxies and quasars accreting at moderate accretion rates are not only powered by MAD scenario but they are also much more powerful than this mechanism predicts. We highlight two possible reasons for this discrepancy. The first one indicates that the geometrical thickness of accretion flows can have much bigger impact on jets formation and evolution than it was originally considered. Here we point out that the future work should concentrate on the existence of thicker than the standard accretion discs as well as on systems with massive, hot coronas. The second solution focuses on the variability of accretion rate which can significantly influence observed Pj/Ld − λEdd distribution. Since estimated by us jet powers are based on the calorimetry on the radio lobes and thereby they are averaged over the lifetime of the source, an important issue to resolve for future studies is to check if using different methods for calculating jet powers such as X-ray cavities (Cavagnolo et al., 2010; Nemmen and Tchekhovskoy, 2015) or model of hotspots (Godfrey and Shabala, 2013) such a trend of decreasing λEdd with increasing Pj/Ld will be still visible.

 Author contributions

Both authors made substantial, direct and intellectual contribution to the work. KR collected and analyzed presented data. MS formulated the original problem, provided direction and guidance. Both authors contributed to the discussion of the obtained results.

Conflict of interest statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

We are particularly grateful to Dorota Kozieł-Wierzbowska and Leith Godfrey for their significant contribution in the analysis of the results and a crucial discussion on their possible explanations. We also thank Greg Madejski for helpful comments.

 References Abramowicz M. A. Chen X. Kato S. Lasota J.-P. Regev O. (1995). Thermal equilibria of accretion disks. Astrophys. J. 438, L37–L39. 10.1086/187709 Avara M. J. McKinney J. C. Reynolds C. S. (2016). Efficiency of thin magnetically arrested discs around black holes. Mon. Not. R. Astron. Soc. 462:636. 10.1093/mnras/stw1643 Begelman M. C. Pringle J. E. (2007). Accretion discs with strong toroidal magnetic fields. Mon. Not. R. Astron. Soc. 375, 1070–1076. 10.1111/j.1365-2966.2006.11372.x Begelman M. C. Armitage P. J. Reynolds C. S. (2015). Accretion disk dynamo as the trigger for X-ray binary state transitions. Astrophys. J. 809:118. 10.1088/0004-637X/809/2/118 Best P. N. Heckman T. M. (2012). On the fundamental dichotomy in the local radio-AGN population: accretion, evolution and host galaxy properties. Mon. Not. R. Astron. Soc. 421, 1569–1582. 10.1111/j.1365-2966.2012.20414.x Blandford R. D. Znajek R. L. (1977). Electromagnetic extraction of energy from Kerr black holes. Mon. Not. R. Astron. Soc. 179, 433–456. 10.1093/mnras/179.3.433 Blundell K. M. Rawlings S. (2000). The spectra and energies of classical double radio lobes. Astron. J. 119, 1111–1122. 10.1086/301254 Cavagnolo K. W. McNamara B. R. Nulsen P. E. J. Carilli C. L. Jones C. Birzan L. (2010). A relationship between AGN jet power and radio power. Astrophys. J. 720, 1066–1072. 10.1088/0004-637X/720/2/1066 Fanaroff B. L. Riley J. M. (1974). The morphology of extragalactic radio sources of high and low luminosity. Mon. Not. R. Astron. Soc. 167:31P. 10.1093/mnras/167.1.31P Fernandes C. A. C. Jarvis M. J. Rawlings S. Martinez-Sansigre A. Hatziminaoglou E. Lacy M. . (2011). Evidence for a maximum jet efficiency for the most powerful radio galaxies. Mon. Not. R. Astron. Soc. 411, 1909–1916. 10.1111/j.1365-2966.2010.17820.x Fernandes C. A. C. Jarvis M. J. Martínez-Sansigre A. Rawlings S. Afonso J. Hardcastle M. J. . (2015). Black hole masses, accretion rates and hot- and cold-mode accretion in radio galaxies at z ~ 1. Mon. Not. R. Astron. Soc. 447, 1184–1203. 10.1093/mnras/stu2517 Godfrey L. E. H. Shabala S. S. (2013). AGN jet kinetic power and the energy budget of radio galaxy lobes. Astrophys. J. 767:12. 10.1088/0004-637X/767/1/12 Igumenshchev I. V. (2008). Magnetically arrested disks and the origin of poynting jets: a numerical study. Astrophys. J. 677, 317–326. 10.1086/529025 Janiuk A. Czerny B. (2011). On different types of instabilities in black hole accretion discs: implications for X-ray binaries and active galactic nuclei. Mon. Not. R. Astron. Soc. 414, 2186–2194. 10.1111/j.1365-2966.2011.18544.x Janiuk A. Czerny B. Siemiginowska A. (2002). Radiation pressure instability driven variability in the accreting black holes. Astrophys. J. 576, 908–922. 10.1086/341804 Kellermann K. I. Sramek R. Schmidt M. Shaffer D. B. Green R. (1989). VLA observations of objects in the Palomar Bright Quasar Survey. Astron. J. 98, 1195–1207. 10.1086/115207 Komissarov S. S. Gubanov A. G. (1994). Relic radio galaxies: evolution of synchrotron spectrum. Astron. Astrophys. 285, 27–43 Laor A. Netzer H. (1989). Massive thin accretion discs. I - Calculated spectra. Mon. Not. R. Astron. Soc. 238, 897–916. 10.1093/mnras/238.3.897 McKinney J. C. Tchekhovskoy A. Blandford R. D. (2012). General relativistic magnetohydrodynamic simulations of magnetically choked accretion flows around black holes. Mon. Not. R. Astron. Soc. 423, 3083–3117. 10.1111/j.1365-2966.2012.21074.x Mullin L. M. Riley J. M. Hardcastle M. J. (2008). Observed properties of FRII quasars and radio galaxies at z < 1.0. Mon. Not. R. Astron. Soc. 390, 595–621. 10.1111/j.1365-2966.2008.13534.x Narayan R. Yi I. (1995). Advection-dominated accretion: self-similarity and bipolar outflows. Astrophys. J. 444, 231–243. 10.1086/175599 Narayan R. Igumenshchev I. V. Abramowicz M. A. (2003). Magnetically arrested disk: an energetically efficient accretion flow. Publ. Astron. Soc. Japan 55, L69–L72. 10.1093/pasj/55.6.L69 Nemmen R. S. Tchekhovskoy A. (2015). On the efficiency of jet production in radio galaxies. Mon. Not. R. Astron. Soc. 449, 316–327. 10.1093/mnras/stv260 Novikov I. D. Thorne K. S. (1973). Astrophysics of black holes, in Black Holes (Les Astres Occlus), eds Dewitt C. Dewitt B. S. (New York, NY: Gordon and Breach Science Publishers), 343–450. Ogle P. Whysong D. Antonucci R. (2006). Spitzer reveals hidden quasar nuclei in some powerful FR II radio galaxies. Astrophys. J. 647:161. 10.1086/505337 Punsly B. (2007). Kinetically dominated FRII radio sources. Mon. Not. R. Astron. Soc. 374, L10–L14. 10.1111/j.1745-3933.2006.00254.x Punsly B. Igumenshchev I. V. Hirose S. (2009). Three-dimensional simulations of vertical magnetic flux in the immediate vicinity of black holes. Astrophys. J. 704, 1065–1085. 10.1088/0004-637X/704/2/1065 Rawlings S. Saunders R. (1991). Evidence for a common central-engine mechanism in all extragalactic radio sources. Nature 349, 138–140. 10.1038/349138a0 Różańska A. Malzac J. Belmont R. Czerny B. Petrucci P.-O. (2015). Warm and optically thick dissipative coronae above accretion disks. Astron. Astrophys. 580:A77. 10.1051/0004-6361/201526288 Rusinek K. Sikora M. Kozieł-Wierzbowska D. Godfrey L. (2017). On the efficiency of jet production in FR II radio galaxies and quasars. Mon. Not. R. Astron. Soc. 466, 2294–2301. 10.1093/mnras/stw3330 Sądowski A. (2016). Thin accretion discs are stabilized by a strong magnetic field. Mon. Not. R. Astron. Soc. 459, 4397–4407. 10.1093/mnras/stw913 Shakura N. I. Sunyaev R. A. (1973). Black holes in binary systems. Observational appearance. Astron. Astrophys. 24, 337–355 Shen Y. Richards G. T. Strauss M. A. Hall P. B. Schneider D. P. Snedden S. . (2011). A catalog of quasar properties from sloan digital sky survey data release 7. Astrophys. J. 194:45. 10.1088/0067-0049/194/2/45 Sikora M. Begelman M. C. (2013). Magnetic flux paradigm for radio loudness of active galactic nuclei. Astrophys. J. 764:L24. 10.1088/2041-8205/764/2/L24 Sikora M. Stawarz Ł. Lasota J.-P. (2007). Radio loudness of active galactic nuclei: observational facts and theoretical implications. Astrophys. J. 658, 815–828. 10.1086/511972 Sikora M. Stasińska G. Kozieł-Wierzbowska D. Madejski G. M. Asari N. V. (2013). Constraining jet production scenarios by studies of narrow-line radio galaxies. Astrophys. J. 765:62. 10.1088/0004-637X/765/1/62 Tchekhovskoy A. Narayan R. McKinney J. C. (2011). Efficient generation of jets from magnetically arrested accretion on a rapidly spinning black hole. Mon. Not. R. Astron. Soc. 418, L79–L83. 10.1111/j.1745-3933.2011.01147.x Yuan F. (2003). Luminous hot accretion flows: thermal equilibrium curve and thermal stability. Astrophys. J. 594, L99–L102. 10.1086/378666 Yuan F. Narayan R. (2014). Hot accretion flows around black holes. Annu. Rev. Astron. Astrophys. 52, 529–588. 10.1146/annurev-astro-082812-141003 van Velzen S. Falcke H. (2013). The contribution of spin to jet-disk coupling in black holes. Astron. Astrophys. 557:L7. 10.1051/0004-6361/201322127 van Velzen S. Falcke H. Körding E. (2015). Nature and evolution of powerful radio galaxies at z ~ 1 and their link with the quasar luminosity function. Mon. Not. R. Astron. Soc. 446, 2985–3001. 10.1093/mnras/stu2213 Willott C. J. Rawlings S. Blundell K. M. Lacy M. (1999). The emission line-radio correlation for radio sources using the 7C Redshift Survey. Mon. Not. R. Astron. Soc. 309, 1017–1033. 10.1046/j.1365-8711.1999.02907.x

Funding. The research leading to these results has received funding from the Polish National Science Centre grant 2016/21/B/ST9/01620.