Draft version August 2, 2024 Typeset using L T E X default style in AASTeX631 A

## Do all the quasars and high-excitation radio galaxies (HERGs) in the 3CRR catalog contain a magnetically arrested disk (MAD)?

Shuang-Liang Li, 1 Wenwen Zuo, 1 and Xinwu Cao 2

1 Key Laboratory for Research in Galaxies and Cosmology, Shanghai Astronomical Observatory, Chinese Academy of Sciences, 80 Nandan Road, Shanghai 200030, People's Republic of China

2 Institute for Astronomy, School of Physics, Zhejiang University, 866 Yuhangtang Road, Hangzhou 310058, People's Republic of China

## ABSTRACT

Based on the magnetization, an accretion disk with large-scale magnetic field can be separated into either standard and normal evolution (SANE) or magnetically arrested disk (MAD), which are difficult to identify from observations. It is still unclear whether all the radio-loud active galactic nuclei (RLAGNs) with a thin disk and strong radio emissions contain a MAD. We investigate this issue by utilizing the 3CRR catalog. We compile a sample of 35 quasars and 14 high-excitation radio galaxies powered by a thin accretion disk. In order to consistently compare with the MAD sample given by Li et al. (2022), the optical-UV emissions of our sample are all detected by the Hubble Space Telescope (HST). It is found that the average X-ray luminosity ( L X ) of our sample is about 5.0 times higher than that of radio-quiet AGNs (RQAGNs) with matching optical-UV luminosity ( L UV ), in general accord with the factor of 4.5 times in MAD sample within the uncertainty. The relationship between radio (5 GHz) and X-ray (2 keV) luminosities in the 3CRR sources is also found to be consistent with the MAD sample. Furthermore, the jet efficiencies of 3CRR sources are consistent with those from the GRMHD simulations of MAD. Therefore, we suggest that probably all the quasars and at least a fraction of high-excitation radio galaxies in the 3CRR catalog, and perhaps all the RLAGNs with strong radio emissions contain a MAD.

Keywords: accretion, accretion disks - black hole physics - magnetohydrodynamics (MHD) - galaxies: active - X-rays: galaxies

## 1. INTRODUCTION

Narayan et al. (2003) suggested that with the accumulation of magnetic flux, the magnetic pressure in the inner accretion disk will be comparable or even higher than the gas pressure, leading to the destruction of symmetry in the inner parts of the accretion disk. This kind of disk is called MAD (magnetically arrested disk), where the magneto-rotational instability (MRI) is suppressed and the magnetic Rayleigh-Taylor (RT) instability will be responsible for the accretion of gas (e.g., Marshall et al. 2018; White et al. 2019). In the general relativistic magnetohydrodynamic (GRMHD) simulations, an advection-domintated accretion flow (ADAF) can be divided into two classes, i.e., the so-called standard and normal evolution (SANE) and MAD according to the magnetization of ADAF (e.g., Tchekhovskoy et al. 2011; McKinney et al. 2012; Chatterjee &amp; Narayan 2022). However, it is very difficult to discriminate between MAD and SANE scenarios from observations. Although the jet power produced in MAD is generally larger than that produced in SANE, we cannot distinguish MAD and SANE based on the jet power alone, because it depends not only on the strength of magnetic field but also on the black hole spin and mass accretion rate (Blandford &amp; Znajek 1977). Additionally, Xie &amp; Zdziarski (2019) suggested that the spectral energy distributions (SED) of SANE and MAD in an ADAF are quite similar, though with a higher radiative efficiency for MAD. Fortunately, the Faraday effect caused by the magnetic field along the direction of light propagation can help to identify MAD from SANE when the polarization information is available. By comparing the rotation measure of jets from observations and simulations of MAD,

Corresponding author: Shuang-Liang Li,Wenwen Zuo, Xinwu Cao

Yuan et al. (2022) suggested that the accretion flow of M87 should be a MAD. Unfortunately, this method cannot be extensively applied to other low-luminosity active galactic nuclei (AGNs) because their polarization is usually weak (see Yuan et al. 2022 for details). The more direct evidence for the presence of MAD was discovered recently in black hole X-ray binary MAXI J1820+070, where the radio fluxes are delayed about 8 days compared with the X-ray flux (You et al. 2023). In this process, the outburst of X-ray flux is produced by the expanding corona, resulting in the amplification of magnetic field and the formation of MAD around the time of the radio peak (8 days later).

Except for the MAD in an ADAF, it is also possible to generate a MAD in an optically thick and geometrically thin accretion disk (e.g., Morales Teixeira et al. 2018; Scepi et al. 2023). By compiling a sample of radio-loud active galaxies, Zamaninasab et al. (2014) discovered a tight correlation between the jet magnetic flux and disk luminosity, which is consistent with the theoretical prediction of a MAD in MHD simulations. A similar relationship between the jet power and the disk luminosity supporting the MAD scenario is also reported in blazars (Ghisellini et al. 2014). However, the samples of these two works include both BL Lacertae objects and flat-spectrum radio quasars (FSRQs), which may correspond to either the MAD in an ADAF (BL Lac) or the MAD in a thin disk (FSRQs), respectively. Furthermore, their physical parameters adopted (e.g., jet power, magnetic flux and disk luminosity) are based on several approximations and assumptions. It is therefore difficult to separate MAD and SANE in a thin disk from these relationships.

Another direct evidence for the MAD in a thin accretion disk comes from the obvious difference between the extreme ultraviolet (EUV) spectrum of radio-loud AGNs (RLAGNs) and the radio-quiet AGNs (RQAGNs, e.g., Telfer et al. 2002), where the RLAGNs show a deficit in EUV band compared with the RQAGNs when the rest-frame wavelength is smaller than 1100 Å. This phenomenon may indicate a missing or suppressed energy dissipation in the innermost region of accretion flow and can be qualitatively explained by the MAD model in a thin disk (Punsly 2014, 2015). For a thin disk with MAD in the innermost region, the presence of large-scale magnetic flux therein will extract a fraction of gas energy to the jet as Poynting flux, resulting in the observed EUV deficit. Furthermore, this picture can also qualitatively explain the positive relationship between the jet power normalized with the bolometric luminosity and the EUV deficit normalized to the RQAGN value shown by HST data (see Punsly 2015 for details).

Except for the EUV spectrum, a significant difference in X-ray emission is also found between the MAD and SANE (Li et al. 2022). The authors reported that the average X-ray flux of a MAD is 4.5 times higher than that of RQAGNs, while the average X-ray flux of RLAGNs without an EUV deficit is about 2-3 times larger than that of RQAGNs (Wu et al. 2013; Gupta et al. 2018). The X-ray excess in MAD sample may originate from the increase of radiative efficiency in a corona, similar with the process that a MAD can improve the radiative efficiency of an ADAF (see Li et al. 2022 and the last section for details). This result implies that not all the RLAGNs contain a MAD. Unfortunately, such a statistically several times higher X-ray flux is not enough to discriminate between SANE and MAD given their error bars. Additionally, the method of utilizing the EUV deficit cannot be employed for most RLAGNs, as the data for EUV emission are often absent. Therefore, we turn our eyes to the radio band and find that all the objects in the MAD sample have strong radio emission (the radio-loudness are all higher than 100, see table 1 in Li et al. 2022). However, it is unclear whether all the RLAGNs with strong radio emission contain a MAD. In this work, we intend to investigate this issue by comparing the multiwavelength properties of the 3CRR sources, the most powerful radio galaxies in the universe, and the MAD sample of Li et al. (2022).

2. SAMPLE Table 1 . Sample.

| Name   | z     | type   | Γ                     | log L X [erg s - 1 Hz - 1 ]   | Ref.   | log L UV [erg s - 1 Hz - 1 ]   | Ref.   | ∆log L X [erg s - 1 Hz - 1 ]   | log L R [erg s - 1 Hz - 1 ]   | log R UV   | log M bh M /circledot   | λ     |
|--------|-------|--------|-----------------------|-------------------------------|--------|--------------------------------|--------|--------------------------------|-------------------------------|------------|-------------------------|-------|
| (1)    | (2)   | (3)    | (4)                   | (5)                           | (6)    | (7)                            | (8)    | (9)                            | (10)                          | (11)       | (12)                    | (13)  |
| 3C9    | 2.012 | Q      | 1.61 +0 . 12 - 0 . 08 | 27.56 +0 . 16 - 0 . 12        | 1      | 31.46                          | 8      | 0.14                           | 32.92                         | 1.46       | 9.8                     | -0.70 |
| 3C47   | 0.425 | Q      | 1.87 +0 . 21 - 0 . 22 | 27.31 +0 . 25 - 0 . 17        | 1      | 30.00                          | 8      | 1.01                           | 32.60                         | 2.61       | 9.2                     | -1.56 |
| 3C48   | 0.367 | Q      | 2.32 +0 . 16 - 0 . 01 | 27.16 +0 . 02 - 0 . 02        | 1      | 30.81                          | 8      | 0.24                           | 33.55                         | 2.74       | 9.2                     | -0.75 |
| 3C68.1 | 1.238 | Q      | 1.9                   | 27.21 +0 . 10 - 0 . 10        | 1      | 30.43                          | 8      | 0.57                           | 31.81                         | 1.38       | 9.9                     | -1.83 |
| 3C138  | 0.759 | Q      | 1.46 +0 . 24 - .      | 27.57 +0 . 10 - .             | 1      | 30.29                          | 8      | 1.04                           | 33.27                         | 2.98       | 8.9                     | -0.97 |

0 11

0 10

| 3C147       | 0.545       | Q   | 1.85 +0 . 24 - 0 . 22   | 26.91 +0 . 24 - 0 . 15        | 1   | 30.33       | 8   | 0.35      | 34.38       | 4.05      | 9.1     | -1.13       |
|-------------|-------------|-----|-------------------------|-------------------------------|-----|-------------|-----|-----------|-------------|-----------|---------|-------------|
| 3C175       | 0.768       | Q   | 1.52 +0 . 40 - 0 . 07   | 27.52 +0 . 08 - 0 . 07        | 1   | 31.03       | 8   | 0.42      | 32.68       | 1.65      | 9.9     | -1.22       |
| 3C181       | 1.382       | Q   | 1.72 +0 . 17 - 0 . 10   | 27.56 +0 . 10 - 0 . 08        | 1   | 30.84       | 8   | 0.61      | 32.66       | 1.82      | 9.6     | -1.12       |
| 3C186       | 1.063       | Q   | 1.88 +0 . 05 - 0 . 03   | 27.29 +0 . 03 - 0 . 03        | 1   | 30.89       | 8   | 0.30      | 32.80       | 1.91      | 9.5     | -0.96       |
| 3C190       | 1.197       | Q   | 1.61 +0 . 06 - 0 . 05   | 27.22 +0 . 04 - 0 . 04        | 1   | 30.51       | 8   | 0.53      | 33.60       | 3.09      | 8.7     | -0.55       |
| 3C191       | 1.952       | Q   | 1.68 +0 . 11 - 0 . 10   | 27.64 +0 . 09 - 0 . 07        | 1   | 30.75       | 8   | 0.76      | 33.82       | 3.07      | 9.7     | -1.31       |
| 3C196       | 0.871       | Q   | 1.67 +0 . 40 - 0 . 38   | 27.23 +0 . 49 - 0 . 25        | 1   | 30.48       | 8   | 0.56      | 32.28       | 1.80      | 9.6     | -1.48       |
| 3C204       | 1.112       | Q   | 2.15 +0 . 14 - 0 . 22   | 27.75 +0 . 16 - 0 . 14        | 1   | 31.07       | 8   | 0.63      | 33.10       | 2.03      | 9.5     | -0.79       |
| 3C205       | 1.534       | Q   | 1.83 +0 . 06 - 0 . 17   | 27.93 +0 . 05 - 0 . 07        | 1   | 31.38       | 8   | 0.57      | 33.28       | 1.90      | 9.6     | -0.58       |
| 3C207       | 0.684       | Q   | 1.59 +0 . 08 - 0 . 23   | 27.39 +0 . 19 - 0 . 15        | 1   | 30.22       | 8   | 0.91      | 33.91       | 3.69      | 8.5     | -0.64       |
| 3C208       | 1.109       | Q   | 1.63 +0 . 18 - 0 . 14   | 27.44 +0 . 10 - 0 . 08        | 1   | 30.99       | 8   | 0.38      | 33.37       | 2.39      | 9.4     | -0.77       |
| 3C212       | 1.049       | Q   | 1.85 +0 . 03 - 0 . 08   | 27.85 +0 . 01 - 0 . 01        | 1   | 30.02       | 8   | 1.52      | 33.79       | 3.77      | 9.2     | -1.54       |
| 3C215       | 0.411       | Q   | 1.88 +0 . 09 - 0 . 08   | 27.16 +0 . 07 - 0 . 05        | 1   | 29.65       | 8   | 1.12      | 31.92       | 2.27      | 8.3     | -1.01       |
| 3C216       | 0.668       | Q   | 1.72 +0 . 26 - 0 . 20   | 27.15 +0 . 08 - 0 . 09        | 1   | 29.85       | 8   | 0.96      | 34.20       | 4.35      | 8.6     | -1.11       |
| 3C245       | 1.029       | Q   | 1.57 +0 . 10 - 0 . 05   | 27.63 +0 . 04 - 0 . 03        | 1   | 30.19       | 8   | 1.17      | 34.55       | 4.37      | 9.4     | -1.57       |
| 3C249.1     | 0.311       | Q   | 1.92 +0 . 13 - 0 . 04   | 26.85 +0 . 08 - 0 . 08        | 1   | 30.25       | 8   | 0.35      | 32.29       | 2.04      | 9.3     | -1.40       |
| 3C254       | 0.734       | Q   | 1.93 +0 . 10 - 0 . 07   | 27.36 +0 . 02 - 0 . 02        | 1   | 30.70       | 8   | 0.52      | 32.55       | 1.84      | 9.3     | -0.95       |
| 3C263       | 0.652       | Q   | 1.83 +0 . 09 - 0 . 05   | 27.27 +0 . 01 - 0 . 01        | 1   | 30.88       | 8   | 0.30      | 33.35       | 2.47      | 9.1     | -0.58       |
| 3C268.4     | 1.400       | Q   | 1.45 +0 . 17 - 0 . 12   | 27.62 +0 . 11 - 0 . 09        | 1   | 31.20       | 8   | 0.40      | 33.59       | 2.39      | 9.8     | -0.96       |
| 3C270.1     | 1.519       | Q   | 1.57 +0 . 07 - 0 . 07   | 27.29 +0 . 03 - 0 . 03        | 1   | 30.87       | 8   | 0.83      | 34.25       | 3.38      | 9.0     | -0.49       |
| 3C275.1     | 0.557       | Q   | 1.79 +0 . 11 - 0 . 10   | 27.02 +0 . 06 - 0 . 06        | 1   | 29.77       | 8   | 0.89      | 33.11       | 3.35      | 8.3     | -0.89       |
| 3C286       | 0.849       | Q   | 2.12 +0 . 51 - 0 . 26   | 27.09 +0 . 12 - 0 . 10        | 1   | 30.75       | 8   | 0.22      | 35.15       | 4.41      | 8.5     | -0.11       |
| 3C287       | 1.055       | Q   | 1.86 +0 . 04 - 0 . 04   | 27.39 +0 . 03 - 0 . 03        | 1   | 30.59       | 8   | 0.64      | 35.10       | 4.51      | 9.6     | -1.37       |
| 3C309.1     | 0.904       | Q   | 1.57 +0 . 05 - 0 . 03   | 27.74 +0 . 05 - 0 . 03        | 1   | 30.81       | 8   | 0.82      | 34.84       | 4.03      | 9.1     | -0.65       |
| 3C325       | 0.860       | Q   | 1.45 +0 . 07 - 0 . 02   | 26.54 +0 . 02 - 0 . 01        | 1   | 29.84       | 8   | 0.36      | 31.80       | 1.96      | 9.6     | -2.12       |
| 3C334       | 0.555       | Q   | 1.85 +0 . 18 - 0 . 11   | 27.12 +0 . 04 - 0 . 04        | 1   | 30.56       | 8   | 0.39      | 33.04       | 2.48      | 9.7     | -1.50       |
| 3C336       | 0.927       | Q   | 1.93 +0 . 24 - 0 . 28   | 27.49 +0 . 14 - 0 . 11        | 1   | 30.58       | 8   | 0.74      | 32.80       | 2.22      | 9.2     | -0.98       |
| 3C345       | 0.594       | Q   | 1.69 +0 . 09 - 0 . 02   | 27.66 +0 . 02 - 0 . 02        | 1   | 30.22       | 8   | 1.18      | 35.00       | 4.78      | 9.3     | -1.44       |
| 3C351       | 0.371       | Q   | 1.90 +0 . 20 - 0 . 13   | 27.02 +0 . 17 - 0 . 09        | 1   | 30.52       | 8   | 0.32      | 31.42       | 0.90      | 9.5     | -1.33       |
| 4C16.49     | 1.296       | Q   | 1.83 +0 . 20 - 0 . 12   | 27.49 +0 . 10 - 0 . 09        | 1   | 30.59       | 8   | 0.74      | 33.02       | 2.43      | 9.8     | -1.57       |
| 3C20        | 0.174       | R   | 1.65 +0 . 21 - 0 . 19   | 25.99 ± 0 . 10                | 2   | 27.24       | 5   | 1.78      | 30.27       | 3.03      | 8.7     | -2.44       |
| 3C33.1      | 0.181       | R   | 0.89 +0 . 33 - 0 . 37   | 26.01 ± 0 . 13                | 2   | 28.37       | 6   | 0.94      | 30.98       | 2.61      | 8.6     | -1.58       |
| 3C61.1      | 0.188       | R   | 1.70                    | 25.90 ± 0 . 13                | 2   | 26.81       | 5   | 2.02      | 30.31       | 3.50      | 8.1     | -0.91       |
| 3C79        | 0.256       | R   | 1.77 +0 69 - 0 . 61     | 25.78 ± 0 . 58                | 3   | 28.87       | 7   | 0.33      | 30.82       | 1.96      | 9.0 8.7 | -1.42       |
| 3C171       | 0.24        | R   | 1.67 +0 . 20 - 0 . 23   | 25.67 ± 0 . 11                | 2   | 26.79       | 5   | 1.80      | 30.13       | 3.34      |         | -1.09       |
| 3C184.1     | 0.118       | R   | 0.57 +0 . 30 - 0 . 34   | 25.39 ± 0 . 16                | 2   | 28.50       | 6   | 0.23      | 30.28       | 1.78      | 8.5     | -1.55       |
| 3C192       | 0.0597      | R   | 1.70                    | 24.89 ± 0 . 19                | 3   | 27.33       | 5   | 0.61      | 29.80       | 2.47      | 8.7     | -2.64       |
| 3C219       | 0.175       | R   | 1.90 +0 . 19 - 0 . 18   | 26.03 ± 0 . 07                | 3   | 28.51       | 6   | 0.86      | 31.57       | 3.06      | 9.2     | -2.71       |
| 3C223       | 0.137       | R   | 0.71 +0 . 77 - 0 . 61   | 25.21 ± 0 . 60                | 3   | 27.67       | 6   | 0.67      | 30.59       | 2.92      | 8.7     | -1.81       |
| 3C234       | 0.185       | R   | 1.39 +0 . 10 - 0 . 16   | 26.21 ± 0 . 09                | 2   | 29.48       | 6   | 0.30      | 31.87 29.93 | 2.40 3.11 | 9.3 8.6 | -1.47       |
| 3C285 3C300 | 0.0794      | R R | 1.70 1.78 ± 0 . 06      | 25.30 ± 0 . 10 25.39 ± 0 . 02 | 3   | 26.82       | 6   | 1.41      | 31.23       | 3.19      | 8.8     | -3.33 -2.07 |
| 3C303       | 0.270 0.141 | R   | 1.60 ± 0 . 21           | 25.84 ± 0 . 08                | 2 3 | 28.04 28.77 | 6 6 | 0.57 0.46 | 31.84       | 3.07      | 9.0     | -2.54       |
| 3C321       | 0.0965      | R   | 1.70                    | 24.06 ± 1 . 32                | 4   | 27.28       | 6   | -0.18     | 30.80       | 3.52      | 9.1     | -3.47       |

Notes: Col. (1): Source name. Col. (2): Redshift. Col. (3): AGN type. Q: quasar; R: radio galaxy. Col. (4): X-ray photon index. Col.(5): X-ray intrinsic luminosity at 2 keV. Col. (6): References for L X . Col. (7): optical-UV luminosity at 2500 Å. Col. (8): References for L UV . Col. (9): Excess of X-ray luminosity in 3CRR sources to that in RQAGNs. Col. (10): Radio spectral luminosity at 5 GHz. Col. (11): Radio loudness defined as R UV = L /L R UV . Col. (12): black hole mass. Col. (13): Eddington ratio.

References. 1, Zhou &amp; Gu (2020); 2, Hardcastle et al. (2009); 3, Hardcastle et al. (2006); 4, Evans et al. (2006); 5, Madrid et al. (2006); 6, Baldi et al. (2010); 7, Fu &amp; Stockton (2008); 8, this work.

The RLAGNs with EUV deficit are suggested to contain a thin accretion disk with MAD in the inner region (Punsly 2014, 2015). In order to compare consistently with the MAD sample, the accretion mode of 3CRR sources is also required to be a thin disk. Therefore, we adopt the quasars and the high-excitation radio galaxies (HERGs) in 3CRR catalog in this work. For one thing, the accretion flow of quasar is widely accepted to be a radiatively efficient thin disk due to the big-blue-bump in optical-UV band, which is thought to originate from an optically thick, geometrically thin disk. For another, Buttiglione et al. (2010) suggested that a radio galaxy with excitation index ( EI = log[O III ] / H β -1 / 3 (log[N II ] / H α +log[S II ] / H α +log[O ] I / H ) α ) higher than 0.95 is a HERG, whose Eddington ratio is found to be higher than 0.01 and is suggested to comprise a thin accretion disk around the black hole (e.g., Best &amp; Heckman 2012).

3CRR catalog contains 43 quasars and 130 radio galaxies. Firstly, we pick up all the quasars and HERGs with detected EI larger than 0.95 in 3CRR catalog, where the values of EI are given by Buttiglione et al. (2010) and Hu et al. (2016). A sample of 43 quasars and 17 radio galaxies are obtained. To avoid uncertainty from different telescopes, only sources observed by Hubble Space Telescope (HST) are chosen, since the optical-UV emission of MAD sample is detected by HST (Punsly 2015; Li et al. 2022). The optical-UV data of 3CRR catalog observed with HST has been extensively explored in many previous works (e.g.,Buttiglione et al. 2010; Westhues et al. 2016; Azadi et al. 2023; Kruk et al. 2023).

In order to obtain their properties in a consistent way, we analyze the HST images of all the 43 quasars ourselves using the Python package photulis (Bradley et al. 2020). First, we estimate the background level and its standard deviation σ , and sigma clip the background data with a threshold of 2 σ . Following this, we use the task detect\_sources to detect sources and generate a segmentation map, providing information of sources including their centroid, semimajor and semi-minor axis, as well as position angle (PA), etc. Using this information, we can derive the surface brightness profile of the target with the task ellipse. By maintaining the ellipticity and PA obtained initially, we refit the isophotes in linear steps to sample the target's surface brightness profile. Finally, we determine the aperture size by the isophote at the surface brightness level of 25 mag arcsec -2 , and perform the force aperture photometry to estimate the flux and uncertainty of the target. By visually checking the fitting results and analyzing the ratio of the obtained flux to the corresponding uncertainty (SNR), we have adopted the results of 35 quasars with SNR larger than 10. The flux at different optical-UV wavelengths (eg. 6060 Å of F606W filter, 4750 Å of F475W filter) are obtained from photometry analysis of the HST images, and then de-reddened for Galactic extinction using the Cardelli et al. (1989) Milky Way reddening law and E(B-V) derived from the Schlegel et al. (1998) dust map. The fluxes at 2500 Å are then calculated by extrapolating the HST optical-UV flux and assuming a spectral index α o = -0 5 . ( f ν ∼ ν α o , e.g., Shang et al. 2011). We fit only the quasars because it is difficult to detect the weak emission lines needed to calculate EI for middle and high redshift HERGs. Additionally, it is challenging to achieve a high SNR in middle and high redshift (z &gt; 0.5) HERG observations due to their faint flux compared to background noise.

Lastly, we get a 3CRR sample composed of 35 quasars and 14 HERGs (see table 1), where all objects have been observed with Chandra and/or XMM-Newton telescope, so the 2 keV X-ray flux is avaible. The photon index Γ in Col. (4) and the intrinsic X-ray luminosity at 2 keV of quasars ( L X , obtained from X-ray spectral modeling) in Col. (5) are all taken from Zhou &amp; Gu (2020), while the L X of HERGs are from NED 1 (see Col. 6 for references). All the L X have been corrected for both intrinsic and Milky Way N H and K-corrected. We can find that, except for three radio galaxies (3C79, 3C223 and 3C321), the uncertainty of X-ray luminosity in other sources are all quite small. The spectral luminosity at 2500 Å( L UV ) of HERGs in Col. (7) are also from NED (see Col. (8) for detailed references), which is derived from the data at nearby frequency assuming a spectral index α o = -0 5 . (Shang et al. 2011). The ∆log L X in Col. (9) indicates the excess of X-ray luminosity above the RQAGN value, calculated from the L UV -L X relationship of Lusso et al. (2010) (for further details see Section 3). The core radio luminosity at 5 GHz in Col. (10) is taken from the online 3CRR catalog . 2 Except for 3C216, 3C325, and 3C345, the black hole mass of quasars and radio galaxies in Col. (12) are all taken from McLure et al. (2006) and Buttiglione et al. (2010), respectively. The black hole mass in 3C216, 3C325, and 3C345 are obtained from Westhues et al. (2016), Azadi et al. (2023), and Shen et al. (2011), respectively. In quasars, the broad line width of H β , C IV , and Mg II is separately adopted to calculate the black hole mass (see McLure et al. 2006 for details), while the host galaxy bulge H band luminosity is used in radio galaxies

1 http://ned.ipac.caltech.edu/

2 https://3crr.extragalactic.info/

/s32

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00321v1-figure-7a4d2b64-4c52-49a6-a810-e25722483843.jpg)

/s32

/s85/s86

Figure 1. Relationship between the X-ray luminosity L X at 2 keV and optical-UV luminosity L UV at 2500 Å in 3CRR sample, where the red circles, blue squares and blue triangles correspond to the MAD, 3CRR quasar and 3CRR HERGs, respectively. The blue line represents our best-fitting result for 3CRR sources. The green line and the red line are the best-fitting result for RQAGNs given by Lusso et al. (2010) and the relationship for MAD sources given by Li et al. (2022), respectively.

(see Buttiglione et al. 2010 for details). These two methods give consistent BH mass estimates (see Fig. 1 in Li &amp; Gu 2021). The Eddington ratio of quasars and radio galaxies in Col. (13) are also estimated by different methods. We use the bolometric correction factor provided by Runnoe et al. (2012) to calculate the bolometric luminosity of quasars, while the calculation of radio galaxy bolometric luminosity is based on the [OIII] λ 5007 luminosity (Buttiglione et al. 2010), as the correction factor for radio galaxies is not available. Indeed, [OIII] luminosity may not be as a reliable proxy for bolometric luminosity for the HERGs as it is for the quasars. This is because the [OIII] emitting region extends to sufficiently small radius (i.e., sufficiently close to the supermassive black hole) to be partly obscured by the dusty torus when the active nucleus is viewed edge-on (see, e.g., Jackson &amp; Browne 1990; Baker 1997; Netzer et al. 2006; Dicken et al. 2014).

## 3. RESULTS

L X -L UV relationship is extensively adopted to investigate the accretion physics in quasars. Lusso et al. (2010) presented a tight relationship between L X and L UV , i.e., log L X = (0 76 . ± 0 02) log . L UV + 3 51 . ± 0 64 . by compiling a type-1 RQAGN sample. This relationship can be qualitatively fitted with the disk-corona model (Lusso &amp; Risaliti 2017; Qiao &amp; Liu 2018), where the X-ray emission comes from the inverse Compton scattering of optical-UV photons from a thin disk by the hot corona. However, though the average X-ray luminosity of MAD sample is found to be roughly 4.5 times larger than that in RQAGNs, matched in L UV , their slopes are very similar (the slope is 0.75, see Li et al. 2022 for details), which may also indicate a disk-corona origin of the L X -L UV relationship in MAD objects.

Not all the RLAGNs contain a MAD since the average X-ray luminosity is only 2-3 times higher than that in RQAGNs (Gupta et al. 2018; Li &amp; Gu 2021). Therefore, in order to check whether all the 3CRR sources with a thin disk and strong radio emissions contain a MAD, we firstly compare their L X -L UV relationship with the MAD sample

/s32

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00321v1-figure-a2967b25-5529-4ad2-9659-07e0e4836521.jpg)

/s32

Figure 2. Relationship between L /L X UV and Eddington ratio λ in 3CRR sample, where the blue and green lines are the best-fitting results for the 3CRR sample and the RQAGN sample given by Lusso et al. (2010), respectively.

given by Li et al. (2022) in Fig. 1, where the blue line represents our best-fitting result. The red and green lines are the best-fitting results for MAD sample (Li et al. 2022) and RQAGNs (Lusso et al. 2010), respectively. It is found that the average X-ray flux of 3CRR sources is about 5.0 times higher than that of RQAGNs, which is consistent with the 4.5 times given in MAD sample within uncertainty. We adopt the ordinary least squares (OLS) method to fit all the relationships for both the 3CRR and the MAD samples (Isobe et al. 1990). The correlation between L X and L UV reads

$$\log L _ { X } = ( 0. 6 8 \pm 0. 0 6 ) \log L _ { U V } + 6. 5 1 \pm 2. 0 2,$$

where the confidence level based on a Pearson test is larger than 99.9%. The slope 0 68 . ± 0 06 . of L X -L UV relationship is shallower but consistent within errors with that in MAD sample ( 0 75 . ± 0 25 . ). However, it is known that the ratio of X-ray luminosity to the UV luminosity will increase with the decreasing Eddington ratio (see e.g., Lusso et al. 2010). Therefore, we check the distribution of Eddington ratio in our sample and find it is equally distributed (Fig. 2). Our best-fitting result (the blue line) is several times higher than that from RQAGNs in Lusso et al. (2010) (the green line), and the confidence level of the relationship is 99.3%. We note that many HERGs have λ &lt; 0 01 . in our sample, which marks the division between a thin disk and a radiatively inefficient accretion flow (Yuan &amp; Narayan 2014). However, similar values for HERGs had been reported by other works with larger samples (see e.g., Fig. 11 in Buttiglione et al. 2010 and Fig. 6 in Best &amp; Heckman 2012). The reasons for this may stem from the significant uncertainties associated with calculating both the black hole mass using bulge H-band luminosity and the bolometric luminosity using [OIII] luminosity.

The traditional relationship between radio (5GHz, L R ) and X-ray (2 keV, L X ) luminosities in AGNs is also investigated in Fig. 3, where the blue line is the best-fitting result, which can be given as

$$\log L _ { R } = ( 1. 3 5 \pm 0. 1 1 ) \log L _ { X } - 3. 4 7 \pm 3. 1 1,$$

/s32

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00321v1-figure-dbc4449a-46e8-48d2-be9c-78f2d0c49820.jpg)

/s32

Figure 3. Relationship between the radio luminosity L R at 5 GHz and the X-ray luminosity L X at 2 keV, where the blue line represents the best-fitting result to the 3CRR sources.

where the confidence level based on a Pearson test is larger than 99.9%. Though the MAD sample is located on the upper right corner of 3CRR objects due to the higher luminosity, they follow the L R -L X relationship (the blue line) in 3CRR objects quite well. The results of Fig. 1 and 3 indicate that all the quasars and many HERGs in the 3CRR catalog may contain a MAD.

In addition, the GRMHD simulations suggested that the jet efficiency ( η = L BZ /Mc ˙ 2 ) in a MAD could be very large, where L BZ and ˙ M are the jet power produce by the BZ mechanism (Blandford &amp; Znajek 1977) and the mass accretion rate, respectively. The jet efficiency η can be larger than 1 in MAD around a black hole with high spin (e.g., Tchekhovskoy et al. 2011; McKinney et al. 2012), where the extra energy is suggested to be extracted from the rotating black hole. In Fig. 4, we investigate the relationship between the jet power L jet and the bolometric luminosity L bol for both the 3CRR sample and MAD sample. Here L jet ( = 0 81 log . L R +11 9 . ) is the kinetic jet power calculated with the radio core luminosity at 5 GHz (Merloni &amp; Heinz 2007). The black solid line represents the maximum jet efficiency η = 140% in a MAD with black hole spin a = 0 99 . and the dotted line is for η = 1 4% . , which roughly corresponds to a = 0 1 . since L BZ ∝ a 2 (Tchekhovskoy et al. 2011). It is found that most of the 3CRR sources and MAD sample lie between the solid and dotted lines. Therefore, the jet efficiencies of 3CRR sources are consistent with the GRMHD simulations of MAD.

## 4. DISCUSSION

In theory, the formation of MAD depends on the accumulation of large-scale magnetic fields in the inner disk region. For a standard thin disk, however, it is difficult to amplify the large-scale magnetic field due to its low radial velocity, resulting in the relatively faster diffusion of magnetic field compared with advection (e.g., Lubow et al. 1994). However, this problem may be resolved by taking the magnetically driven outflow into account, which can help to increase the radial velocity of gas by transferring the angular momentum of disk (Cao &amp; Spruit 2013; Li &amp; Begelman 2014). Furthermore, the outflow takes away not only the angular momentum, but also the energy of gas, resulting in a

/s32

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00321v1-figure-48a0efb9-0c17-493a-b05c-1c7f85bcfd7e.jpg)

/s32

Figure 4. Relationship between the jet power L jet and the bolometric luminosity L bol , where the solid line and the dotted line represent the maximum jet efficiency ( η = 140% with a = 0 99 . ) and the jet efficiency η = 1 4% . with a = 0 1 . , respectively (Tchekhovskoy et al. 2011).

far smaller radiative efficiency. The radiative efficiency given by theoretical model can be several order of magnitudes smaller than 10% in extreme cases (see e.g., Li 2014; Li &amp; Cao 2019; Zhao et al. 2023). This result can explain the phenomenon that the jet efficiency of some objects in Fig. 4 are even several times higher than 140%. The detailed physics of a disk-corona system with MAD is still unclear, though the MHD simulations do suggest a higher jet power and smaller bolometric luminosity (Tchekhovskoy et al. 2011; McKinney et al. 2012; Morales Teixeira et al. 2018). However, for an ADAF, the presence of MAD can increase its radiative efficiency (Xie &amp; Zdziarski 2019). If the properties of corona is similar with ADAF as suggested by some works (Meyer et al. 2000; Liu et al. 2007; Qiao &amp; Liu 2018), we can anticipate that the X-ray flux will also increase several times for a MAD in thin disk.

Recent works suggest that the magnetically elevated accretion disks (e.g., Begelman et al. 2022; Begelman &amp; Armitage 2023) can also lead to a MAD. In this scenario, the authors suggested that it is the strong toroidal magnetic field making an accretion disk MAD, though weak poloidal field is still necessary to trigger the dynamo process. However, the jet power in a magnetically elevated accretion disk should be much weaker compared with the standard MAD scenario, since the poloidal field responsible for the jet power is weak ( L BZ ∝ B 2 p , see e.g., Ghosh &amp; Abramowicz 1997).

The 3CRR catalog contains objects with the strongest radio emissions at 178 MHz in the sky, so the black hole spins in these objects are possibly high, leading to high jet efficiencies. Although the jet in the objects around the dotted line in Fig. 4 can be launched by a SANE with high black hole spin, other objects around the solid line should contain a MAD judging by their high jet efficiency. Furthermore, the EUV deficit in the MAD sample is difficult to explain with the SANE scenario. The origin of the X-ray luminosity excess remains debatable. In principle, both the jet core and the corona, influenced by the large-scale magnetic field, can contribute to this excess. However, if this excess were attributed to the jet core, two specific conditions must be met to maintain a similar L X -L UV slope between MAD and RQAGNs: 1) the jet core must depend strongly on L UV ; and 2) the slope of this dependence must be consistent with

the L X -L UV slope in RQAGNs (see Zhu et al. 2020 for details). Therefore, it is more natural to assume a coronal origin for the X-ray emissions. However, the jet core and corona may be indistinguishable in spatial resolution. In this context, one could also argue that the X-ray luminosity excess originates from the jet. This slope is also consistent with that observed in radio-loud quasars (see e.g., Zhu et al. 2020; Li &amp; Gu 2021).

## 5. SUMMARY

We investigate whether all RLAGNs with strong radio emission harbor a MAD by compiling a sample of quasars and HERGs from 3CRR catalog in this work. By comparing the multi-band data with the MAD sample given in Li et al. (2022), it is found that the average X-ray luminosity of 3CRR objects is higher by the same factor as the MAD objects when compared with the RQAGNs matched in L UV from Lusso et al. (2010) (Fig. 1). The L X -L UV slope of the 3CRR sample is also consistent with that of the MAD sample. Furthermore, the 3CRR sample has a larger range of L X (4 dex) than the MAD sample (1 dex), so the L X -L UV relation and slope are much better determined than those from the MAD sample. We find a positive relationship between the radio (5 GHz) luminosity L R and the 2 keV X-ray luminosity L X for the 3CRR sample (Fig. 3). The MAD sample is found to follow the L R -L X relationship quite well, which further support the same physical mechanism for MAD sample and 3CRR sources. Lastly, the jet efficiencies of the 3CRR sources (Fig. 4) are also found to be qualitatively consistent with GRMHD simulations of MAD. These results indicate that all the quasars and a fraction of HERGs (at least for z&lt;0.3) in 3CRR catalog should contain a MAD.

## ACKNOWLEDGEMENTS

We thank the reviewer for very valuable and detailed comments, which have greatly improved the quality of this manuscript. This work is supported by the NSFC (grants 12273089, 12073023, 12233007, 12361131579, 12347103), the science research grants from the China Manned Space Project with No. CMSCSST-2021-A06, and the Fundamental Research Fund for Chinese Central Universities. This work has made extensive use of the NASA/IPAC Extragalactic Database (NED) and the 3CRR catalog (Laing et al. 1983). NED is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration (NASA). This research has also made extensive use of data obtained from the Chandra Data Archive and the observations made with the NASA/ESA Hubble Space Telescope, obtained from the Hubble Legacy Archive, which is a collaboration between the Space Telescope Science Institute (STScI/NASA), the Space Telescope European Coordinating Facility (ST-ECF/ESAC/ESA) and the Canadian Astronomy Data Centre (CADC/NRC/CSA).

## REFERENCES

```
<_SQL_>
```

```

/*   ...
```

| Ghosh, P., & Abramowicz, M. A. 1997, MNRAS, 292, 887, doi: 10.1093/mnras/292.4.887                                                         |
|--------------------------------------------------------------------------------------------------------------------------------------------|
| Gupta, M., Sikora, M., Rusinek, K., & Madejski, G. M. 2018, MNRAS, 480, 2861, doi: 10.1093/mnras/sty2043                                   |
| Hardcastle, M. J., Evans, D. A., & Croston, J. H. 2006, MNRAS, 370, 1893, doi: 10.1111/j.1365-2966.2006.10615.x -. 2009, MNRAS, 396, 1929, |
| Hu, J.-F., Cao, X.-W., Chen, L., & You, B. 2016, Research in Astronomy and Astrophysics, 16, 136, doi: 10.1088/1674-4527/16/9/136          |
| Isobe, T., Feigelson, E. D., Akritas, M. G., & Babu, G. J. 1990, ApJ, 364, 104, doi: 10.1086/169390                                        |
| Jackson, N., & Browne, I. W. A. 1990, Nature, 343, 43, doi: 10.1038/343043a0                                                               |
| Kruk, S., García-Martín, P., Popescu, M., et al. 2023, Nature Astronomy, 7, 262, doi: 10.1038/s41550-023-01903-3                           |
| Laing, R. A., Riley, J. M., & Longair, M. S. 1983, MNRAS, 204, 151, doi: 10.1093/mnras/204.1.151                                           |
| Li, J., & Cao, X. 2019, ApJ, 886, 92, doi: 10.3847/1538-4357/ab4c36                                                                        |
| Li, S.-L. 2014, ApJ, 788, 71, doi: 10.1088/0004-637X/788/1/71                                                                              |
| Li, S.-L., & Begelman, M. C. 2014, ApJ, 786, 6, doi: 10.1088/0004-637X/786/1/6                                                             |
| Li, S.-L., & Gu, M. 2021, A&A, 654, A141, doi: 10.1051/0004-6361/202141301                                                                 |
| doi: 10.1051/0004-6361/202244105 B. F., Taam, R. E., Meyer-Hofmeister, E., & F.                                                            |
| Liu, Meyer, 2007, ApJ, 671, 695, doi: 10.1086/522619                                                                                       |
| Lubow, S. H., Papaloizou, J. C. B., & Pringle, J. E. 1994, MNRAS, 267, 235, doi: 10.1093/mnras/267.2.235                                   |
| Lusso, E., & Risaliti, G. 2017, A&A, 602, A79, doi: 10.1051/0004-6361/201630079                                                            |
| Lusso, E., Comastri, A., Vignali, C., et al. 2010, A&A, 512, A34, doi: 10.1051/0004-6361/200913298                                         |
| Madrid, J. P., Chiaberge, M., Floyd, D., et al. 2006, ApJS, 164, 307, doi: 10.1086/504480                                                  |
| Marshall, M. D., Avara, M. J., & McKinney, J. C. 2018, MNRAS, 478, 1837, doi: 10.1093/mnras/sty1184                                        |
| McKinney, J. C., Tchekhovskoy, A., & Blandford, R. D. 2012, MNRAS, 423, 3083, doi: 10.1111/j.1365-2966.2012.21074.x                        |
| McLure, R. J., Jarvis, M. J., Targett, T. A., Dunlop, J. S., & Best, P. N. 2006, MNRAS, 368, 1395, doi: 10.1111/j.1365-2966.2006.10228.x   |

| Merloni, A., & Heinz, S. 2007, MNRAS, 381, 589, doi: 10.1111/j.1365-2966.2007.12253.x                                    |
|--------------------------------------------------------------------------------------------------------------------------|
| Meyer, F., Liu, B. F., & Meyer-Hofmeister, E. 2000, A&A, 361, 175, doi: 10.48550/arXiv.astro-ph/0007091                  |
| Morales Teixeira, D., Avara, M. J., & McKinney, J. C. 2018, MNRAS, 480, 3547, doi: 10.1093/mnras/sty2044                 |
| Narayan, R., Igumenshchev, I. V., & Abramowicz, M. A. 2003, PASJ, 55, L69, doi: 10.1093/pasj/55.6.L69                    |
| Netzer, H., Mainieri, V., Rosati, P., & Trakhtenbrot, B. 2006, A&A, 453, 525, doi: 10.1051/0004-6361:20054203            |
| Punsly, B. 2014, ApJL, 797, L33, doi: 10.1088/2041-8205/797/2/L33 -. 2015, ApJ, 806, 47, doi: 10.1088/0004-637X/806/1/47 |
| Qiao, E., & Liu, B. F. 2018, MNRAS, 477, 210, doi: 10.1093/mnras/sty652                                                  |
| Runnoe, J. C., Brotherton, M. S., & Shang, Z. 2012, MNRAS, 422, 478, doi: 10.1111/j.1365-2966.2012.20620.x               |
| Scepi, N., Begelman, M. C., & Dexter, J. 2023, arXiv e-prints, arXiv:2302.10226, doi: 10.48550/arXiv.2302.10226          |
| Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525, doi: 10.1086/305772                                 |
| Shang, Z., Brotherton, M. S., Wills, B. J., et al. 2011, ApJS, 196, 2, doi: 10.1088/0067-0049/196/1/2                    |
| Shen, Y., Richards, G. T., Strauss, M. A., et al. 2011, ApJS, 194, 45, doi: 10.1088/0067-0049/194/2/45                   |
| Tchekhovskoy, A., Narayan, R., & McKinney, J. C. 2011, MNRAS, 418, L79, doi: 10.1111/j.1745-3933.2011.01147.x            |
| Telfer, R. C., Zheng, W., Kriss, G. A., & Davidsen, A. F. 2002, ApJ, 565, 773, doi: 10.1086/324689                       |
| Westhues, C., Haas, M., Barthel, P., et al. 2016, AJ, 151, 120, doi: 10.3847/0004-6256/151/5/120                         |
| White, C. J., Stone, J. M., & Quataert, E. 2019, ApJ, 874, 168, doi: 10.3847/1538-4357/ab0c0c                            |
| Wu, J., Brandt, W. N., Miller, B. P., et al. 2013, ApJ, 763, 109, doi: 10.1088/0004-637X/763/2/109                       |
| Xie, F.-G., & Zdziarski, A. A. 2019, ApJ, 887, 167, doi: 10.3847/1538-4357/ab5848                                        |
| You, B., Cao, X., Yan, Z., et al. 2023, Science, 381, 961, doi: 10.1126/science.abo4504                                  |
| Yuan, F., & Narayan, R. 2014, ARA&A, 52, 529, doi: 10.1146/annurev-astro-082812-141003                                   |
| Yuan, F., Wang, H., & Yang, H. 2022, ApJ, 924, 124, doi: 10.3847/1538-4357/ac4714                                        |
| Zamaninasab, M., Clausen-Brown, E., Savolainen, T., & Tchekhovskoy, A. 2014, Nature, 510, 126, doi: 10.1038/nature13399  |
| Zhao, Y., Yang, X.-H., Xue, L., & Li, S.-L. 2023, MNRAS, 526, 862, doi: 10.1093/mnras/stad2816                           |

Zhou, M., &amp; Gu, M. 2020, ApJ, 893, 39, doi: 10.3847/1538-4357/ab7dca

Zhu, S. F., Brandt, W. N., Luo, B., et al. 2020, MNRAS, 496, 245, doi: 10.1093/mnras/staa1411
Shuang-Liang Li Wenwen Zuo Xinwu Cao
Do all the quasars and high-excitation radio galaxies (HERGs) in the 3CRR catalog contain a magnetically arrested disk (MAD)?
Based on the magnetization, an accretion disk with large-scale magnetic field
can be separated into either standard and normal evolution (SANE) or
magnetically arrested disk (MAD), which are difficult to identify from
observations. It is still unclear whether all the radio-loud active galactic
nuclei (RLAGNs) with a thin disk and strong radio emissions contain a MAD. We
investigate this issue by utilizing the 3CRR catalog. We compile a sample of 35
quasars and 14 high-excitation radio galaxies powered by a thin accretion disk.
In order to consistently compare with the MAD sample given by Li et al. (2022),
the optical-UV emissions of our sample are all detected by the Hubble Space
Telescope (HST). It is found that the average X-ray luminosity ($L_{\rm X}$) of
our sample is about 5.0 times higher than that of radio-quiet AGNs (RQAGNs)
with matching optical-UV luminosity ($L_{\rm UV}$), in general accord with the
factor of 4.5 times in MAD sample within the uncertainty. The relationship
between radio (5~GHz) and X-ray (2 keV) luminosities in the 3CRR sources is
also found to be consistent with the MAD sample. Furthermore, the jet
efficiencies of 3CRR sources are consistent with those from the GRMHD
simulations of MAD. Therefore, we suggest that probably all the quasars and at
least a fraction of high-excitation radio galaxies in the 3CRR catalog, and
perhaps all the RLAGNs with strong radio emissions contain a MAD.