Relation between winds and jets in radio-loud AGN
Missagh Mehdipour Elisa Costantini






SRON Netherlands Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, the Netherlands
 M.Mehdipour@sron.nl








We investigate the relation between the two modes of outflow (wind and jet) in radio-loud active galactic nuclei (AGN).

For this study we have carried out a systematic and homogeneous analysis of XMM-Newton spectra of a sample of 16 suitable radio-loud Seyfert-1 AGN. 

The ionised winds in these AGN are parameterised through high-resolution X-ray spectroscopy and photoionisation modelling.

We discover a significant inverse correlation between the column density N_ H of the ionised wind and the radio-loudness parameter R of the jet.

We explore different possible explanations for this N_ H-R relation and find that ionisation, inclination, and luminosity effects are unlikely to be responsible for the observed relation.

We argue that the N_ H-R relation is rather a manifestation of the magnetic driving mechanism of the wind from the accretion disk. 

Change in the magnetic field configuration from toroidal to poloidal, powering either the wind or the jet mode of the outflow, is the most feasible explanation for the observed decline in the wind as the radio jet becomes stronger.

Our findings provide evidence for a wind-jet bimodality in radio-loud AGN and shine new light on the link between these two modes of outflow. This has far-reaching consequences for the accretion disk structure and the wind ejection mechanism.




M. Mehdipour & E. Costantini

Relation between winds and jets in radio-loud AGN

Relation between winds and jets in radio-loud AGN
    
Missagh Mehdipour
Elisa Costantini

    Received 5 February 2019 / Accepted 26 March 2019
=====================================================






§ INTRODUCTION



The growth of supermassive black holes (SMBHs) via accretion is linked to the growth of their host galaxy via star formation. This is observed as a correlation between the SMBH mass and the stellar mass of the galactic bulge (the so-called M-σ relation, e.g. ). This suggests that SMBHs and their host galaxies are co-evolved. However, the mechanism required for this co-evolution is rather uncertain. In AGN, accretion onto the SMBH liberates enough power to affect the host galaxy, but how this is transported from the very small scales close to the SMBH to the much larger scales of the host galaxy is poorly understood.

Apart from the radiation, accretion is ubiquitously accompanied by outflows (jets and/or winds), which transport matter and energy away from the nucleus into the galaxy (see e.g. the recent reviews by ). These outflows can provide an efficient coupling between SMBHs and their host galaxies. They can impact star formation, chemical enrichment of the surrounding intergalactic medium, and the cooling flows at the core of galaxy clusters (e.g. ). AGN outflows can therefore play an important role in the co-evolution of SMBHs and their host galaxies through a feedback mechanism. However, there are significant gaps in our understanding of the outflow phenomenon in AGN, which cause uncertainties in determining their role and impact in galaxy evolution.

Observations show that about 15–20% of AGN are radio-loud <cit.>, where the ratio of radio to optical emission is orders of magnitude higher than in typical radio-quiet AGN. In radio-loud AGN accretion onto the SMBH involves energy release in the form of radio synchrotron radiation from collimated jets of relativistic particles from the central engine. The mechanism for the production of the jet has been thought to be electromagnetic extraction of rotational energy from a spinning black hole via magnetic field lines that are supported by the accretion disk <cit.>. However, the connection between the accretion disk and the jet, and the formation of the observed dichotomy of radio-quiet and radio-loud AGN, are still open questions. Apart from the collimated jets, another important mode/component of AGN outflows is the winds of ionised gas. Importantly, the relationship of the ionised winds to the jet and the disk is not known, which we investigate in this paper.

While jets are generally understood to be magnetically powered, the origin and launching mechanism of the ionised winds in AGN are more uncertain. The winds may originate from either the accretion disk or the torus of the AGN, with different possible mechanisms postulated for their ejection. They could be either thermally-driven (e.g. ), or radiatively-driven (e.g. ), or magnetically-driven (e.g. ) winds. However, it is still debatable what physical conditions and factors govern the launch of outflows, and what is the dependence on the accretion flows, radiative power, and the mass and the spin of SMBHs. 

The X-ray energy band is where the ionised winds mainly reveal their diagnostic-rich spectral features. The nearby X-ray bright radio-quiet AGN have served as excellent laboratories to study these winds (e.g. ). However, there are also detections of ionised winds in radio-loud AGN (e.g. ).

In an X-ray study of a sample of Type-1 AGN using Advanced Satellite for Cosmology and Astrophysics (ASCA, ), <cit.> suggested that the four radio-loud AGN in their sample had relatively weaker absorption compared to the radio-quiet AGN. However, the number of radio-loud AGN, and the spectral resolution of ASCA, were too limited for a conclusive interpretation. More recently, in a search for ultra-fast outflows at hard X-rays using <cit.> and <cit.>, <cit.> find evidence of Fe K absorption with relativistic speeds in some radio-loud AGN. They report that the ultra-fast outflows are detected at a wide range of jet inclination angles (∼ 10 to 70), suggesting that the winds have a large opening angle and are not preferentially equatorial.

In this paper we investigate the link between the parameters of the wind and the jet in a sample of Type-1 AGN that are X-ray bright and radio loud. We systematically model all the X-ray absorption by the ionised winds in these objects using (RGS and EPIC spectroscopy), and derive the wind parameters from our homogeneous photoionisation modelling. This is the first time changes in the ionised AGN wind as a function of the jet radio power are investigated. Our characterisation of the wind in radio-loud AGN provides new information for testing theoretical models of the wind-jet-disk connection in AGN.

The structure of the paper is as follows. In Sect. <ref> we describe the selection criteria of our AGN sample, how the properties of our AGN sample were collated, as well as the reduction and processing of the data. The fitting of the X-ray spectra and photoionisation modelling of the AGN wind are reported in Sect. <ref>. The wind-jet correlation analysis from the results of our modelling are reported in Sect. <ref>. We discuss and interpret our findings in Sect. <ref>, and give concluding remarks in Sect. <ref>.

The spectral analysis and modelling presented in this paper were done using the SPEX package <cit.> v3.04.00 and its pion photoionisation model <cit.>. We use C-statistics for spectral fitting with the X-ray spectra optimally binned according to <cit.>. Errors on the derived parameters are reported at the 1σ confidence level. We assume proto-solar abundances of <cit.> throughout our modelling in this paper. In our computations we adopt the cosmological parameters H_0=70 km s^-1 Mpc^-1, Ω_Λ=0.70, and Ω_ m=0.30.






§ OBSERVATIONAL DATA







 §.§ Sample selection




We started with the Quasars and Active Galactic Nuclei Catalogue (13th Ed.) of <cit.> for selecting our AGN and cross matching with observations. This catalogue provides a collection of useful information on a large number of known AGN (about 169 thousands), such as the spectral classification (e.g. Seyfert 1, BL Lac, etc.), the 6 cm (5 GHz) radio flux density F_ 6  cm, as well as the optical magnitude. From these we calculated the radio-loudness parameter R <cit.>, which is defined as the ratio of radio to optical flux density: R = F_ 6  cm / F_ opt. For the calculation of F_ opt in units of erg cm^-2 s^-1 Hz^-1 from the B magnitude (at 4400 Å) we used the relation B = -2.5 log F_ opt - 48.36 <cit.> in <cit.>. The radio-loud AGN are formally classed as those with R > 10 <cit.>.

We cross matched the AGN in <cit.> catalogue with the Serendipitous Source Catalogue 3XMM-DR6 (XMM-SSC, 2016), in order to identify the associated X-ray emission of the AGN in the <cit.> catalogue, and acquire their X-ray fluxes. For the cross matching of the two catalogues we used the CDS X-Match Service <cit.> provided by the Centre de Données astronomiques de Strasbourg (CDS). The cross matching of the two catalogues resulted in about 8800 matched X-ray sources. Following this, we proceeded to filter out those objects classed by <cit.> as BL-Lacs, highly polarised objects, LINERs, and Type-2 AGN. This left us with about 2400 Type-1 AGN remaining, which consist of Seyfert-1 galaxies and quasars. We then proceeded to select only the radio-loud AGN. Thus, for the objects with available F_ 6  cm and B data, those with R > 10 were selected. This gave us 99 radio-loud Type-1 AGN. The distribution of the X-ray flux of these AGN is shown in the histogram of Fig. <ref>. In order to model the AGN wind and derive its parameters, high-resolution X-ray spectroscopy of the absorption features in the soft X-ray band is needed. Therefore, RGS spectrum with sufficient signal-to-noise ratio (S/N) is required. We thus selected only those AGN with S/N of at least 3 in the RGS resolution element (about 70 mÅ) at the continuum level. This means those AGN with extreme neutral absorption (either Galactic or in the host galaxy), or intrinsically faint X-ray continuum, or those with insufficient exposure are effectively filtered out. The 16 remaining AGN that are suitable for the wind study are highlighted in green in the histogram of Fig. <ref>. The properties of these selected radio-loud Type-1 AGN are listed in Table <ref>.
















 §.§ Properties of the selected AGN




Following the selection of our AGN sample in Sect. <ref>, we compiled information on properties of these AGN from the literature. These properties which are listed in Table <ref>, alongside the results of our modelling of the AGN wind (Sect. <ref>), are required for our investigation of the wind-jet-disk connection. We describe below how the properties in Table <ref> were obtained.

The classification, the cosmological redshift z, and the B magnitude of all the AGN in our sample are taken from the <cit.> catalogue and references therein. The radio and X-ray fluxes were initially taken from <cit.> and the Serendipitous Source Catalogue, respectively, for the purpose of selecting the AGN (Sect. <ref>). However, these flux measurements were further refined for our study in this paper, which are provided in Table <ref> with their references. The soft (0.3–2 keV) and hard (2–10 keV) observed X-ray fluxes, and the intrinsic X-ray luminosity in Table <ref>, are derived from our spectral modelling presented in Sect. <ref>. 

In order to assess the strength of emission from the AGN dusty torus we obtained the mid-IR flux of the objects in our sample. We used mid-IR flux measurements from Multi-band Imaging Photometer (MIPS) taken in the 24 μm photometric channel. The MIPS fluxes were retrieved from the Spitzer Enhanced Imaging Products (SEIP) source list through NASA/IPAC Infrared Science Archive (IRSA). For two of the AGN in our sample (ESO 075-G041 and Mrk 896), which /MIPS observations were not available, we used mid-IR flux measurements from Infrared Astronomical Satellite (IRAS) taken in the 25 μm photometric channel. The IRAS fluxes were obtained from the IRAS Faint Source catalogue v2.0 <cit.>. The mid-IR flux (F_ mid  IR) of the AGN in our sample are provided in Table <ref>. 

We also collected from the literature the inclination angle i of the radio jet axis relative to our line of sight, and the mass of the SMBH (M_ BH) in the AGN. For one of the objects in our sample (Mrk 896), where the jet inclination angle is not available, we instead take the reported inclination angle of the accretion disk. The i and M_ BH values of our AGN sample, and their reference papers, are provided in Table <ref>. The intrinsic bolometric luminosity L_ bol of the AGN was estimated using a spectral energy distribution (SED) model. As the objects in our sample are radio-loud Seyfert-1 AGN, we use a radio-loud version of the SED model derived from a multi-wavelength campaign study of the archetypal Seyfert-1 AGN NGC 5548 <cit.>. We modified this radio-quiet SED, which spans from near-IR to hard X-rays, by extending it to lower energies using a power-law model for the synchrotron emission from the jet. In our model this power-law continuum spans from near-IR (1 μm) to radio (100 MHz) energies with a fiducial photon index of 1.8. For radio-loudness parameter R ranging from 10 to 10^3, the corresponding luminosity contribution by the power-law extension ranges from 10^42 to 10^44 . This radio-loud broadband SED model was used to estimate the bolometric luminosity L_ bol by matching the intrinsic X-ray luminosity L_ X and the radio flux F_ 6  cm of each AGN in our sample. Finally, from M_ BH and L_ bol, the Eddington ratio λ was calculated. The L_ X, L_ bol, and the λ values of the AGN in our sample are shown in Table <ref>.




 §.§ X-ray data reduction and processing



The Observation Data Files (ODFs) of the AGN in our sample were retrieved from the Science Archive (XSA) and processed using the Science Analysis System (SAS v17.0.0). The ID and the duration of the observations that we used in our study are provided in Table <ref>. For those AGN with multiple observations we have used one with sufficient S/N for line spectroscopy of the wind with RGS. The data from the RGS instruments <cit.> were processed through the rgsproc pipeline task; the source and background spectra were extracted and the response matrices were generated. We filtered out time intervals with background count rates > 0.1 count s^-1 in CCD number 9. The rgscombine task was used to stack the first-order RGS1 and RGS2 spectra of each observation. The fitted spectral range is 7–36 Å for RGS. The EPIC-pn <cit.> and EPIC-MOS <cit.> data were processed with the epproc and emproc pipeline tasks, respectively. Periods of high flaring background for EPIC-pn (> 0.40 count s^-1) and EPIC-MOS (> 0.35 count s^-1) were filtered out by applying the #XMMEA_EP and #XMMEA_EM filters, respectively. The EPIC spectra were extracted from a circular region centred on the source with a radius of 40”. The background was extracted from a nearby source-free region of the same size. The pileup was evaluated to be negligible for majority of the AGN in our sample. However, the largest pileup for the brightest object (3C 273) was found to be about 8% in MOS and 4% in pn. Nonetheless, for our purpose of the line spectroscopy of the wind in the soft band pileup is not an issue. The single and double events were selected for the pn (PATTERN <= 4) and MOS (PATTERN <= 12). Response matrices were generated for the EPIC spectra of each observation using the rmfgen and arfgen tasks. The epicspeccombine task was used for stacking the MOS1 and MOS2 spectra of each observation. Our fitted range for the EPIC spectra is from 0.3 to 10 keV. In our spectral modelling (Sect. <ref>) in , we simultaneously fit the RGS, EPIC-pn, and EPIC-MOS spectra of each observation. Our fitting takes into account differences in their instrumental flux normalisation, which ranges from about few % to 10%.






§ SPECTRAL MODELLING OF THE X-RAY DATA







 §.§ X-ray continuum




The continuum model that we applied to fit the X-ray spectra in consists of three components: a power-law component (pow), a modified black body component (mbb), and a neutral reflection component (refl). The pow component models the primary broadband continuum, while the mbb component models the commonly seen `soft X-ray excess', which is an excess emission at ≤ 2 keV in the X-ray spectra of AGN (e.g. ). The mbb model in is a blackbody component modified by coherent Compton scattering <cit.>. The refl model fits the line in the EPIC spectra, which is part of the reflection spectrum produced by the reprocessing of the primary X-ray power-law by neutral gas in the AGN. The refl model in computes the line according to <cit.>, and the Compton-reflected continuum according to <cit.>, as described in <cit.>. The normalisation and the photon index Γ of the incident power-law for refl were set to those of the observed power-law component (pow). The exponential high-energy cut-off of pow and the incident power-law for refl were fixed to a fiducial value of 200 keV. In our modelling we fitted the reflection scale s parameter of the refl model. The ionisation parameter of refl was set to zero to produce a neutral reflection component, and its abundances were kept at their default solar values. The best-fit parameters of the continuum components (pow, mbb, and refl) for each AGN in our sample are shown in Table <ref>.




 §.§ X-ray absorption by the diffuse ISM




In our spectral modelling we take into account the continuum and line absorption by the diffuse interstellar medium (ISM) in the Milky Way and the host galaxy of the AGN. This is done using the hot model in <cit.>. This model calculates the transmission of a plasma in collisional ionisation equilibrium at a given temperature, which for neutral ISM is set to the minimum temperature of the model at 0.5 eV. We applied one hot component for the Milky Way and another one for the host galaxy of the AGN. The total Galactic in our line of sight towards each object (provided in Table <ref>) was obtained from <cit.> and kept fixed in our modelling. The host galaxy was allowed to be fitted, and its best-fit value (or upper limit) for each object is shown in Table <ref>.




 §.§ Photoionisation and X-ray spectral modelling of the ionised wind in AGN




The photoionisation modelling of the ionised AGN wind was carried out in using the pion model <cit.>. We used a template SED for the calculation of the photoionisation equilibrium. The ionising SED of the Seyfert-1 AGN in our sample, extending from optical/UV to hard X-rays, can be reasonably approximated by the SED of the archetypal Seyfert-1 AGN NGC 5548, which was derived in <cit.> from continuum modelling of simultaneous multi-wavelength data. The results of the pion calculations were then inputted into the xabs model in , which calculates the absorption spectrum of a slab of photoionised gas, where all ionic column densities are linked in a physically consistent fashion according to the pion run. The parameters of the xabs model are the ionisation parameter ξ, the hydrogen column density , the covering fraction C_f, and the outflow v_ out and the turbulent σ_v velocities. The ionisation parameter ξ <cit.> is defined as ξ≡L  / n_H r^2, where L is the luminosity of the ionising source over 1–1000 Ryd (13.6 eV to 13.6 keV) in erg s^-1, n_H the hydrogen density in cm^-3, and r the distance between the ionised gas and the ionising source (in cm). The covering fraction C_f was fixed to unity in our modelling.

We started with one xabs component and fitted its , ξ, and v_ out parameters. This resulted in significant improvement to the spectral fit of all AGN in our sample. For some of the AGN we found that after fitting the first xabs component, there were still remaining absorption features at either different ionisation states or outflow velocities that were not fitted. Therefore, for those cases we fitted a second xabs component. Two xabs components were found to be sufficient to fit all the absorption features, and a third component did not improve the fit and hence was not needed. We kept σ_v of the xabs components at the default 100 in our modelling in order to limit the number of unnecessary free parameters and also because fitting σ_v did not significantly improve the fits. The best-fit parameters of the xabs model for the ionised wind in our AGN sample are reported in Table <ref>.























§ CORRELATION ANALYSIS OF THE MODELLING RESULTS




Here we examine the results of our modelling of the ionised winds in radio-loud AGN, carried out in Sect. <ref>. We look for correlations between the parameters of the wind, jet, disk, and the torus (see Tables <ref>, <ref>, and <ref>) in order to find out any physical link between them in radio-loud AGN. For our correlation analysis we calculated the Spearman's rank correlation coefficient (r_s) and the associated null hypothesis p-value probability (p_ null) from the two-sided t-test.

Figure <ref> shows the total column density of the ionised wind plotted versus the radio-loudness parameter R of the jet in our AGN sample. Interestingly, it is apparent that decreases as R increases. This inverse correlation is statistically significant with r_s = -0.83 and probability p_ null = 0.00007. The data in Fig. <ref> are fitted with a power-law function, given by


    = 8.2× 10^22 R^-0.8


where is the total column density of the ionised wind in , and the radio-loudness parameter of the jet R = F_ 6  cm / F_ opt.


There can be different possible explanations for the discovered -R relation, which we will assess and discuss in Sect. <ref>. In order to aid us identify the physical reason, we investigate the relation between the wind and other key parameters in Fig. <ref>. These parameters are the ionisation parameter ξ of the wind (top left), intrinsic X-ray luminosity of the AGN (top right), the radio jet inclination angle relative to our line of sight (middle left), the ratio of mid-IR to optical flux (middle right), SMBH mass (bottom left), and the Eddington ratio λ (bottom right). The r_s correlation coefficient and the null hypothesis probability p_ null are shown in the corner of each panel. The results show no significant correlation between of the wind and the aforementioned parameters. In Sect. <ref> we discuss the -R relation and the implications of the lack of correlation between and the other parameters.












§ DISCUSSION







 §.§ Origin of the observed NH-R relation in radio-loud AGN




In Sect. <ref> we found that there is a significant inverse correlation between the wind and the radio-loudness parameter R (Fig. <ref>). This shows that as the radio jet becomes stronger, the ionised wind becomes weaker. Importantly, we found no significant correlation between the wind and the other key parameters shown in Fig. <ref>.

The fact that is not related to the ionisation parameter ξ and the X-ray luminosity L_ X, implies that the observed -R relation is not a manifestation of changes in the ionisation state of the wind. This means the lowering of the detected in the soft X-ray band is not caused by a shift in ξ of the gas to higher ionisation states. If the -R relation was a result of changes in ξ, then an inverse relation between and ξ would be expected, which is not seen (Fig. <ref>, top left panel). This is also supported by the fact that there is no correlation between and the ionising X-ray luminosity L_ X of the AGN (Fig. <ref>, top right panel). We also find no correlation between ξ and the outflow velocity v_ out of the wind.

Another possibility for the -R relation could be the effect of the viewing angle of the AGN. According to some theoretical wind models, from high to low inclination angles (i.e. from edge-on to face-on view), the wind is expected to decrease (e.g. ). However, the inclination angle i of the objects in our sample, which are all Seyfert-1 AGN, are rather limited to a small range (10 to 40). Moreover, we do not find a correlation between and i (Fig. <ref>, middle left panel). Therefore, the inclination angle is unlikely to be the cause of the observed lowering in the wind as R increases.

One of the mechanisms for the launch of ionised winds in AGN is thermal evaporation from the AGN dusty torus (e.g. ). Moreover, such winds from the torus may be boosted by the IR radiation pressure on dust grains in the wind (e.g. ). In order to find out if there is a relation between the ionised wind and the dusty torus in our AGN sample, we checked for correlations between the wind parameters and the emission from the dusty torus, which peaks in the mid-IR energy band. Similar to the radio-loudness parameter R of the jet, the `mid-IR loudness' parameter of the torus is a measure of the relative strength of the emission from the torus compared to the emission from the accretion disk, which we define as the ratio of mid-IR to optical (B) flux density (F_ mid  IR / F_ opt). We find no significant correlation between F_ mid  IR / F_ opt and the wind parameters, including , which is shown in Fig. <ref>, middle right panel. Recent mid-IR interferometry observations have shown that significant fraction of the mid-IR emission in some AGN originates along the polar directions (e.g. ). This is thought to be associated to dusty winds from the AGN torus. However, the lack of a correlation between and the mid-IR loudness parameter, as well as the no correlation between and the inclination angle i, indicate that the AGN torus is unlikely to play a key role in the observed variation in our sample and hence the observed -R relation.

The strong accretion-powered radiation in AGN can drive outflowing winds, provided the gas is not over-ionised (e.g. ). For such radiatively-driven winds, the mass outflow rate would be dependent on the mass accretion rate Ṁ and hence the Eddington ratio λ of the AGN (e.g. ). Interestingly, an inverse relation between the radio-loudness parameter R and the Eddington ratio λ has been previously observed in radio-loud AGN (e.g. ). This relation shows that the radio power of the jet decreases with increasing Ṁ and λ. However, this relation is evident over a wide range of λ (spanning several orders of magnitude), significantly more than the sub-Eddington λ range of our AGN sample (spanning 1.7 orders of magnitude) as shown in Table <ref> and Fig. <ref> (bottom right panel). Moreover, we do not find a statistically significant correlation between R and λ for our AGN sample. Similarly, there is no significant correlation between the wind and λ, nor between and the black hole mass M_ BH, as shown in Fig. <ref> (bottom panels). This means λ of our selected AGN is effectively sampling a bin of the wider R-λ relation reported in the literature. Therefore, the above suggest that our observed -R relation is not a consequence of large Ṁ and λ variations in our AGN sample.

Since the different scenarios discussed earlier cannot provide a feasible explanation for the observed -R relation, we suggest that the most plausible explanation is associated to the only mechanism that can launch both jets and winds, i.e. the magnetic fields. Magnetic pressure generated in the accretion disk can drive winds (e.g. ), while magnetocentrifugal acceleration along the field lines can drive jets aided by the spin of the black hole (e.g. ). Thus, depending on the black hole spin and the magnetic field configuration, the magnetic power maybe distributed differently over the jet and wind modes of the outflow. We discuss the magnetic explanation for the wind-jet -R relation in Sect. <ref>.




 §.§ Wind-jet connection by the magnetic fields in radio-loud AGN




In Sect. <ref> we deduced that the -R relation is most likely the manifestation of the magnetic-powered outflow mechanism in radio-loud AGN. Interestingly, in black hole binaries (BHBs) the wind properties have been observed to change with the spectral state <cit.>. Observations show that while the wind is present in the high/soft state, it disappears in the low/hard state. The magnetic outflow mechanism has been proposed as a feasible explanation for the inverse relation between the winds and jets in BHBs <cit.>. In this scenario, the wind and the jet are driven by different configurations of the same magnetic field. In other words, the same magnetically-driven outflow is either a wind in the high/soft state, or a jet in the low/hard state (e.g. ).

The wind-jet inverse relation phenomenon in BHBs is analogous to the one discovered in radio-loud AGN by our investigation. As discussed by <cit.>, the magnetic field configuration may change from toroidal to poloidal in transition from high/soft disk-dominated state to low/hard jet-dominated state. Thus, the magnetic outflow mechanism can provide the required explanation for the observed -R relation in radio-loud AGN. This is also further supported by the study of <cit.>, where they report that black hole winds and jets are regulated across the mass scale with a common launching mechanism driving both outflows, i.e. magnetohydrodynamics for winds and magnetocentrifugal forces for jets.

An important consequence of the above scenario is the likely changing geometry of the accretion disk with increasing jet activity. Disk truncation has been discussed by <cit.> as an explanation for the jet cycle seen in 3C 120. From a multi-epoch case-study of this AGN in the radio, optical/UV, and X-rays, <cit.> find a complete disk at one epoch, and a truncated/refilling disk at another epoch. The first epoch is seen as a rising trend in the optical/UV and X-ray emission from the disk, while the second epoch is associated to a dip in the disk emission and an increase in the radio jet activity. The behaviour of the ionised wind in 3C 120 has not been investigated in <cit.>.







Finally, in Fig. <ref> we illustrate the relation between the wind and the jet in radio-loud AGN based on the findings of our investigation. In radio-loud AGN with spinning black holes, stronger poloidal magnetic field and disk truncation give rise to brighter radio jets (i.e. higher R), and hence less toroidal magnetic winds are launched from the accretion disk (i.e. lower ).

In this study we considered only radio-loud AGN for investigating the wind-jet connection. However, in radio-quiet AGN, where winds are prominently present, there is also radio emission, albeit much fainter than in radio-loud AGN. At the low threshold limit of radio-loudness parameter (R ≈ 10), we find that the wind reaches about 10^22 (Fig. <ref>), which is a typical seen in radio-quiet Seyfert-1 AGN (see e.g. ). In radio-quiet AGN there is also a wide range of wind (e.g. ), which its association to the jet activity is yet to be established. This is however challenging because the nature and origin of radio emission in radio-quiet AGN is rather uncertain. There are different possible explanations given in the literature, such as coronal disk emission, star formation, emission from outflows including low-power jets, or perhaps a combination of them (see e.g.  and references therein). Therefore, while ionised winds are commonly seen in radio-quiet AGN, investigating their link to the jet activity is complex since multiple processes may contribute to the faint radio emission in radio-quiet AGN.

The X-ray spectroscopy of ionised winds in radio-loud AGN is currently feasible for a limited number of X-ray bright sources as presented in this study (Fig. <ref>). The upcoming X-ray observatory <cit.>, with its unprecedented X-ray sensitivity and energy resolution, enables us to extend the wind-jet study to a larger population of AGN that are not bright enough for the current X-ray observatories. It will also help us to simultaneously constrain all components of the wind using high-resolution spectroscopy in the soft and hard X-ray bands. The increase in the sample size will provide a more general characterisation of the wind-jet-disk connection in radio-loud AGN. 






§ CONCLUSIONS




In this paper we have parameterised the ionised winds in a sample of 16 radio-loud AGN by carrying out high-resolution X-ray spectroscopy and photoionisation modelling. The results of our study shine new light on the relation between the two modes of outflow (wind and jet) in AGN. From the findings of our investigation we conclude the following. 





  * We discover that in radio-loud AGN there is a significant inverse correlation between the column density of the ionised wind and the radio-loudness parameter R of the jet. The wind , ranging over ∼ 10^20–10^22 , is seen to decrease as a power-law function with R, ranging over ∼ 10–10^3. There is no significant correlation between R and the other wind parameters, namely the ionisation parameter ξ and the outflow velocity v_ out.  



  * The results of our study show that the observed -R relation is not caused by inclination or ionisation effects. There is also no association to the black hole mass and the Eddington luminosity ratio of the AGN. The -R relation is also independent of the AGN dusty torus.



  * The most feasible explanation for the origin of the observed -R relation is the magnetic outflow mechanism in radio-loud AGN. Depending on the black hole spin and the magnetic field configuration (toroidal or poloidal), the magnetic power drives the wind or the jet mode of the outflow differently. Disk truncation is also a likely consequence of this scenario as the jet becomes stronger and the wind becomes weaker.



  * The observed wind-jet bimodality in radio-loud AGN is analogous to that previously found in stellar-mass black holes. Our results point towards the magnetohydrodynamic mechanism for the ejection of ionised winds from the accretion disk in AGN.









M.M. and E.C. are supported by the Netherlands Organisation for Scientific Research (NWO) through the Innovational Research Incentives Scheme Vidi grant 639.042.525. This work is based on observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and the USA (NASA). This research made use of the cross-match service provided by CDS, Strasbourg. We thank R. Morganti for the useful discussions. We thank the anonymous referee for the constructive comments.







[Ballantyne(2005)]Ball05b
Ballantyne, D. R. 2005, , 362, 1183

[Ballantyne & Fabian(2005)]Ball05a
Ballantyne, D. R. & Fabian, A. C. 2005, , 622, L97

[Begelman et al.(1983)Begelman, McKee, & Shields]Bege83
Begelman, M. C., McKee, C. F., & Shields, G. A. 1983, , 271, 70

[Behar et al.(2017)Behar, Peretz, Kriss, Kaastra, Arav,
  Bianchi, Branduardi-Raymont, Cappi, Costantini, De Marco, Di
  Gesu, Ebrero, Kaspi, Mehdipour, Paltani, Petrucci, Ponti, &
  Ursini]Beh17
Behar, E., Peretz, U., Kriss, G. A., et al. 2017, , 601, A17

[Bicay et al.(1995)Bicay, Kojoian, Seal, Dickinson, &
  Malkan]Bica95
Bicay, M. D., Kojoian, G., Seal, J., Dickinson, D. F., & Malkan,
  M. A. 1995, The Astrophysical Journal Supplement Series, 98, 369

[Blandford & Znajek(1977)]Blan77
Blandford, R. D. & Znajek, R. L. 1977, , 179, 433

[Blandford & Payne(1982)]Blan82
Blandford, R. D. & Payne, D. G. 1982, , 199, 883

[Blustin et al.(2005)Blustin, Page, Fuerst,
  Branduardi-Raymont, & Ashton]Blu05
Blustin, A. J., Page, M. J., Fuerst, S. V., Branduardi-Raymont, G., &
  Ashton, C. E. 2005, , 431, 111

[de Plaa et al.(2004)de Plaa, Kaastra, Tamura,
  Pointecouteau, Mendez, & Peterson]dePl04
de Plaa, J., Kaastra, J. S., Tamura, T., et al. 2004, , 423, 49

[den Herder et al.(2001)den Herder, Brinkman, Kahn,
  Branduardi-Raymont, Thomsen, , et al.]denH01
den Herder, J. W., Brinkman, A. C., Kahn, S. M., et al. 2001, ,
  365, L7

[Di Gesu & Costantini(2016)]DiGe16
Di Gesu, L. & Costantini, E. 2016, , 594, A88

[Dodson et al.(2008)Dodson, Fomalont, Wiik, Horiuchi,
  Hirabayashi, Edwards, Murata, Asaki, Moellenbrock, Scott,
  Taylor, Gurvits, Paragi, Frey, Shen, Lovell, Tingay, Rioja,
  Fodor, Lister, Mosoni, Coldwell, Piner, & Yang]Dods08
Dodson, R., Fomalont, E. B., Wiik, K., et al. 2008, , 175, 314

[Dorodnitsyn et al.(2011)Dorodnitsyn, Bisnovatyi-Kogan, &
  Kallman]Doro11
Dorodnitsyn, A., Bisnovatyi-Kogan, G. S., & Kallman, T. 2011, , 741,
  29

[Doroshenko et al.(2012)Doroshenko, Sergeev, Klimanov,
  Pronik, & Efimov]Doro12
Doroshenko, V. T., Sergeev, S. G., Klimanov, S. A., Pronik, V. I., &
  Efimov, Y. S. 2012, , 426, 416

[Fabian(2012)]Fabi12
Fabian, A. C. 2012, , 50, 455

[Ferrarese & Merritt(2000)]Ferr00
Ferrarese, L. & Merritt, D. 2000, , 539, L9

[Foschini et al.(2015)Foschini, Berton, Caccianiga, Ciroi,
  Cracco, Peterson, Angelakis, Braito, Fuhrmann, Gallo, Grupe,
  Järvelä, Kaufmann, Komossa, Kovalev, Lähteenmäki,
  Lisakov, Lister, Mathur, Richards, Romano, Sievers,
  Tagliaferri, Tammi, Tibolla, Tornikoski, Vercellone, La Mura,
  Maraschi, & Rafanelli]Fosc15
Foschini, L., Berton, M., Caccianiga, A., et al. 2015, , 575, A13

[Fukumura et al.(2010)Fukumura, Kazanas, Contopoulos, &
  Behar]Fuku10
Fukumura, K., Kazanas, D., Contopoulos, I., & Behar, E. 2010, ,
  715, 636

[Gear et al.(1994)Gear, Stevens, Hughes, Litchfield,
  Robson, Terasranta, Valtaoja, Steppe, Aller, & Aller]Gear94
Gear, W. K., Stevens, J. A., Hughes, D. H., et al. 1994, , 267,
  167

[Giovannini et al.(2001)Giovannini, Cotton, Feretti, Lara,
  & Venturi]Giov01
Giovannini, G., Cotton, W. D., Feretti, L., Lara, L., & Venturi, T.
  2001, , 552, 508

[Grier et al.(2017)Grier, Pancoast, Barth, Fausnaugh,
  Brewer, Treu, & Peterson]Grie17
Grier, C. J., Pancoast, A., Barth, A. J., et al. 2017, , 849, 146

[Higginbottom et al.(2014)Higginbottom, Proga, Knigge,
  Long, Matthews, & Sim]Higg14
Higginbottom, N., Proga, D., Knigge, C., et al. 2014, , 789, 19

[Ho(2002)]Ho02
Ho, L. C. 2002, , 564, 120

[Hogan et al.(2011)Hogan, Lister, Kharb, Marshall, &
  Cooper]Hoga11
Hogan, B. S., Lister, M. L., Kharb, P., Marshall, H. L., & Cooper,
  N. J. 2011, , 730, 92

[Hönig et al.(2013)Hönig, Kishimoto, Tristram,
  Prieto, Gandhi, Asmus, Antonucci, Burtscher, Duschl, &
  Weigelt]Honi13
Hönig, S. F., Kishimoto, M., Tristram, K. R. W., et al. 2013, ,
  771, 87

[Jansen et al.(2001)Jansen, Lumb, Altieri, Clavel, Ehle,
  , et al.]Jans01
Jansen, F., Lumb, D., Altieri, B., et al. 2001, , 365, L1

[Jorstad et al.(2005)Jorstad, Marscher, Lister, Stirling,
  Cawthorne, Gear, Gómez, Stevens, Smith, Forster, &
  Robson]Jors05
Jorstad, S. G., Marscher, A. P., Lister, M. L., et al. 2005, , 130,
  1418

[Kaastra & Barr(1989)]Kaa89
Kaastra, J. S. & Barr, P. 1989, , 226, 59

[Kaastra et al.(1996)Kaastra, Mewe, &
  Nieuwenhuijzen]Kaa96
Kaastra, J. S., Mewe, R., & Nieuwenhuijzen, H. 1996, in UV and X-ray
  Spectroscopy of Astrophysical and Laboratory Plasmas, ed. K. Yamashita &
  T. Watanabe, 411–414

[Kaastra et al.(2000)Kaastra, Mewe, Liedahl, Komossa, &
  Brinkman]Kaa00
Kaastra, J. S., Mewe, R., Liedahl, D. A., Komossa, S., & Brinkman,
  A. C. 2000, , 354, L83

[Kaastra et al.(2012)Kaastra, Detmers, Mehdipour, Arav,
  Behar, Bianchi, Branduardi-Raymont, Cappi, Costantini, Ebrero,
  Kriss, Paltani, Petrucci, Pinto, Ponti, Steenbrugge, & de
  Vries]Kaas12
Kaastra, J. S., Detmers, R. G., Mehdipour, M., et al. 2012, , 539,
  A117

[Kaastra & Bleeker(2016)]Kaas16
Kaastra, J. S. & Bleeker, J. A. M. 2016, , 587, A151

[Kaastra et al.(2017)Kaastra, Raassen, de Plaa, & Gu]Kaas17
Kaastra, J. S., Raassen, A. J. J., de Plaa, J., & Gu, L. 2017, SPEX X-ray
  spectral fitting package, https://doi.org/10.5281/zenodo.2272992

[Karamanavis(2015)]Karam15
Karamanavis, V. . 2015, PhD thesis, Max-Planck-Institut für
  Radioastronomie

[Kellermann et al.(1989)Kellermann, Sramek, Schmidt,
  Shaffer, & Green]Kell89
Kellermann, K. I., Sramek, R., Schmidt, M., Shaffer, D. B., & Green,
  R. 1989, , 98, 1195

[Kharb et al.(2006)Kharb, O'Dea, Baum, Colbert, &
  Xu]Khar06
Kharb, P., O'Dea, C. P., Baum, S. A., Colbert, E. J. M., & Xu, C.
  2006, , 652, 177

[King et al.(2013)King, Miller, Raymond, Fabian,
  Reynolds, Gültekin, Cackett, Allen, Proga, &
  Kallman]King13
King, A. L., Miller, J. M., Raymond, J., et al. 2013, , 762, 103

[King & Pounds(2015)]King15
King, A. & Pounds, K. 2015, , 53, 115

[Königl & Kartje(1994)]Koni94
Königl, A. & Kartje, J. F. 1994, , 434, 446

[Krolik et al.(1981)Krolik, McKee, & Tarter]Kro81
Krolik, J. H., McKee, C. F., & Tarter, C. B. 1981, , 249, 422

[Krolik & Kriss(2001)]Krol01
Krolik, J. H. & Kriss, G. A. 2001, , 561, 684

[Kuehr et al.(1981)Kuehr, Pauliny-Toth, Witzel, &
  Schmidt]Kueh81
Kuehr, H., Pauliny-Toth, I. I. K., Witzel, A., & Schmidt, J. 1981,
  , 86, 854

[Kurosawa & Proga(2009)]Kuro09
Kurosawa, R. & Proga, D. 2009, , 397, 1791

[Laha et al.(2014)Laha, Guainazzi, Dewangan, Chakravorty,
  & Kembhavi]Laha14
Laha, S., Guainazzi, M., Dewangan, G. C., Chakravorty, S., &
  Kembhavi, A. K. 2014, , 441, 2613

[Laurent-Muehleisen et al.(1997)Laurent-Muehleisen, Kollgaard,
  Ryan, Feigelson, Brinkmann, & Siebert]Laur97
Laurent-Muehleisen, S. A., Kollgaard, R. I., Ryan, P. J., et al. 1997,
  , 122, 235

[Lodders et al.(2009)Lodders, Palme, & Gail]Lod09
Lodders, K., Palme, H., & Gail, H.-P. 2009, Landolt Börnstein, 44

[Lohfink et al.(2013)Lohfink, Reynolds, Jorstad, Marscher,
  Miller, Aller, Aller, Brenneman, Fabian, Miller, Mushotzky,
  Nowak, & Tombesi]Lohf13
Lohfink, A. M., Reynolds, C. S., Jorstad, S. G., et al. 2013, ,
  772, 83

[Magdziarz & Zdziarski(1995)]Magd95
Magdziarz, P. & Zdziarski, A. A. 1995, , 273, 837

[Marchesini et al.(2004)Marchesini, Celotti, &
  Ferrarese]Marc04
Marchesini, D., Celotti, A., & Ferrarese, L. 2004, , 351, 733

[Mehdipour et al.(2011)Mehdipour, Branduardi-Raymont,
  Kaastra, Petrucci, Kriss, , et al.]Meh11
Mehdipour, M., Branduardi-Raymont, G., Kaastra, J. S., et al. 2011,
  , 534, A39

[Mehdipour et al.(2015)Mehdipour, Kaastra, Kriss, Cappi,
  Petrucci, Steenbrugge, Arav, Behar, Bianchi, Boissay,
  Branduardi-Raymont, Costantini, Ebrero, Di Gesu, Harrison, Kaspi,
  De Marco, Matt, Paltani, Peterson, Ponti, Pozo Nuñez, De
  Rosa, Ursini, de Vries, Walton, & Whewell]Meh15a
Mehdipour, M., Kaastra, J. S., Kriss, G. A., et al. 2015, , 575,
  A22

[Mehdipour et al.(2016)Mehdipour, Kaastra, &
  Kallman]Meh16b
Mehdipour, M., Kaastra, J. S., & Kallman, T. 2016, , 596, A65

[Mehdipour & Costantini(2018)]Meh18b
Mehdipour, M. & Costantini, E. 2018, , 619, A20

[Miller et al.(2012)Miller, Raymond, Fabian, Reynolds,
  King, Kallman, Cackett, van der Klis, & Steeghs]Mill12
Miller, J. M., Raymond, J., Fabian, A. C., et al. 2012, , 759, L6

[Mingo et al.(2014)Mingo, Hardcastle, Croston, Dicken,
  Evans, Morganti, & Tadhunter]Ming14
Mingo, B., Hardcastle, M. J., Croston, J. H., et al. 2014, , 440,
  269

[Mitsuda et al.(2007)Mitsuda, Bautz, Inoue, Kelley,
  Koyama, Kunieda, Makishima, Ogawara, Petre, Takahashi, Tsunemi,
  White, Anabuki, Angelini, Arnaud, Awaki, Bamba, Boyce, Brown,
  Chan, Cottam, Dotani, Doty, Ebisawa, Ezoe, Fabian, Figueroa,
  Fujimoto, Fukazawa, Furusho, Furuzawa, Gendreau, Griffiths,
  Haba, Hamaguchi, Harrus, Hasinger, Hatsukade, Hayashida, Henry,
  Hiraga, Holt, Hornschemeier, Hughes, Hwang, Ishida, Ishisaki,
  Isobe, Itoh, Iyomoto, Kahn, Kamae, Katagiri, Kataoka,
  Katayama, Kawai, Kilbourne, Kinugasa, Kissel, Kitamoto, Kohama,
  Kohmura, Kokubun, Kotani, Kotoku, Kubota, Madejski, Maeda,
  Makino, Markowitz, Matsumoto, Matsumoto, Matsuoka, Matsushita,
  McCammon, Mihara, Misaki, Miyata, Mizuno, Mori, Mori, Morii,
  Moseley, Mukai, Murakami, Murakami, Mushotzky, Nagase, Namiki,
  Negoro, Nakazawa, Nousek, Okajima, Ogasaka, Ohashi, Oshima,
  Ota, Ozaki, Ozawa, Parmar, Pence, Porter, Reeves, Ricker,
  Sakurai, Sanders, Senda, Serlemitsos, Shibata, Soong, Smith,
  Suzuki, Szymkowiak, Takahashi, Tamagawa, Tamura, Tamura,
  Tanaka, Tashiro, Tawara, Terada, Terashima, Tomida, Torii,
  Tsuboi, Tsujimoto, Tsuru, Turner, Ueda, Ueno, Ueno, Uno,
  Urata, Watanabe, Yamamoto, Yamaoka, Yamasaki, Yamashita,
  Yamauchi, Yamauchi, Yaqoob, Yonetoku, & Yoshida]Mits07
Mitsuda, K., Bautz, M., Inoue, H., et al. 2007, , 59, S1

[Morganti(2017)]Morg17
Morganti, R. 2017, Frontiers in Astronomy and Space Sciences, 4, 42

[Moshir et al.(1990)Moshir, Copan, Conrow, McCallon,
  Hacking, Gregorich, Rohrbach, Melnyk, Rice, & Fullmer]Mosh90
Moshir, M., Copan, G., Conrow, T., et al. 1990, in IRAS Faint Source
  Catalogue, version 2.0 (1990)

[Nandra et al.(2013)Nandra, Barret, Barcons, Fabian, den
  Herder, Piro, Watson, Adami, Aird, Afonso, & et al.]Nand13
Nandra, K., Barret, D., Barcons, X., et al. 2013, ArXiv e-prints
  [[arXiv]1306.2307]

[Neilsen & Lee(2009)]Neil09
Neilsen, J. & Lee, J. C. 2009, , 458, 481

[Ohsuga et al.(2009)Ohsuga, Mineshige, Mori, &
  Kato]Ohsu09
Ohsuga, K., Mineshige, S., Mori, M., & Kato, Y. 2009, , 61, L7

[Page et al.(2003)Page, O'Brien, Reeves, &
  Breeveld]Page03
Page, K. L., O'Brien, P. T., Reeves, J. N., & Breeveld, A. A. 2003,
  , 340, 1052

[Paltani & Türler(2005)]Palt05
Paltani, S. & Türler, M. 2005, , 435, 811

[Panessa et al.(2019)Panessa, Baldi, Laor, Padovani,
  Behar, & McHardy]Pane19
Panessa, F., Baldi, R. D., Laor, A., et al. 2019, arXiv e-prints
  [[arXiv]1902.05917]

[Peterson et al.(2004)Peterson, Ferrarese, Gilbert, Kaspi,
  Malkan, , et al.]Pet04
Peterson, B. M., Ferrarese, L., Gilbert, K. M., et al. 2004, , 613,
  682

[Pineau et al.(2011)Pineau, Boch, & Derriere]Pine11
Pineau, F.-X., Boch, T., & Derriere, S. 2011, in Astronomical Society of
  the Pacific Conference Series, Vol. 442, Astronomical Data Analysis Software
  and Systems XX, ed. I. N. Evans, A. Accomazzi, D. J. Mink, & A. H.
  Rots, 85

[Piotrovich et al.(2015)Piotrovich, Gnedin, Silant'ev,
  Natsvlishvili, & Buliga]Piot15
Piotrovich, M. Y., Gnedin, Y. N., Silant'ev, N. A., Natsvlishvili,
  T. M., & Buliga, S. D. 2015, , 454, 1157

[Ponti et al.(2012)Ponti, Fender, Begelman, Dunn,
  Neilsen, & Coriat]Pont12
Ponti, G., Fender, R. P., Begelman, M. C., et al. 2012, , 422,
  L11

[Proga & Kallman(2004)]Prog04
Proga, D. & Kallman, T. R. 2004, , 616, 688

[Proga(2007)]Prog07
Proga, D. 2007, , 661, 693

[Reeves et al.(2009)Reeves, Sambruna, Braito, &
  Eracleous]Reev09b
Reeves, J. N., Sambruna, R. M., Braito, V., & Eracleous, M. 2009,
  , 702, L187

[Reynolds(1997)]Reyn97
Reynolds, C. S. 1997, , 286, 513

[Rokaki et al.(2003)Rokaki, Lawrence, Economou, &
  Mastichiadis]Roka03
Rokaki, E., Lawrence, A., Economou, F., & Mastichiadis, A. 2003,
  , 340, 1298

[Schmidt & Green(1983)]Schm83
Schmidt, M. & Green, R. F. 1983, , 269, 352

[Shimmins & Bolton(1972)]Shim72
Shimmins, A. J. & Bolton, J. G. 1972, Australian Journal of Physics
  Astrophysical Supplement, 23, 1

[Sikora et al.(2007)Sikora, Stawarz, & Lasota]Siko07
Sikora, M., Stawarz, Ł., & Lasota, J.-P. 2007, , 658, 815

[Stawarz(2004)]Staw04
Stawarz, Ł. 2004, , 613, 119

[Steenbrugge et al.(2005)Steenbrugge, Kaastra, Crenshaw,
  Kraemer, Arav, George, Liedahl, van der Meer, Paerels, Turner,
  & Yaqoob]Stee05
Steenbrugge, K. C., Kaastra, J. S., Crenshaw, D. M., et al. 2005, ,
  434, 569

[Strüder et al.(2001)Strüder, Briel, Dennerl,
  Hartmann, Kendziorra, , et al.]Stru01
Strüder, L., Briel, U., Dennerl, K., et al. 2001, , 365, L18

[Tanaka et al.(1994)Tanaka, Inoue, & Holt]Tana94
Tanaka, Y., Inoue, H., & Holt, S. S. 1994, , 46, L37

[Tarter et al.(1969)Tarter, Tucker, & Salpeter]Tar69
Tarter, C. B., Tucker, W. H., & Salpeter, E. E. 1969, , 156, 943

[Tombesi et al.(2014)Tombesi, Tazaki, Mushotzky, Ueda,
  Cappi, Gofford, Reeves, & Guainazzi]Tomb14
Tombesi, F., Tazaki, F., Mushotzky, R. F., et al. 2014, , 443,
  2154

[Torresi et al.(2012)Torresi, Grandi, Costantini, &
  Palumbo]Torr12
Torresi, E., Grandi, P., Costantini, E., & Palumbo, G. G. C. 2012,
  , 419, 321

[Tristram et al.(2014)Tristram, Burtscher, Jaffe,
  Meisenheimer, Hönig, Kishimoto, Schartmann, &
  Weigelt]Tris14
Tristram, K. R. W., Burtscher, L., Jaffe, W., et al. 2014, , 563,
  A82

[Turner et al.(2001)Turner, Abbey, Arnaud, Balasini,
  Barbera, , et al.]Turn01
Turner, M. J. L., Abbey, A., Arnaud, M., et al. 2001, , 365, L27

[Urry & Padovani(1995)]Urry95
Urry, C. M. & Padovani, P. 1995, , 107, 803

[Vasudevan & Fabian(2007)]Vasu07
Vasudevan, R. V. & Fabian, A. C. 2007, , 381, 1235

[Veilleux et al.(2009)Veilleux, Kim, Rupke, Peng,
  Tacconi, Genzel, Lutz, Sturm, Contursi, Schweitzer, Dasyra,
  Ho, Sanders, & Burkert]Veil09
Veilleux, S., Kim, D.-C., Rupke, D. S. N., et al. 2009, , 701, 587

[Véron-Cetty & Véron(2010)]Vero10
Véron-Cetty, M. P. & Véron, P. 2010, , 518, A10

[Willingale et al.(2013)Willingale, Starling, Beardmore,
  Tanvir, & O'Brien]Will13
Willingale, R., Starling, R. L. C., Beardmore, A. P., Tanvir, N. R., &
  O'Brien, P. T. 2013, , 431, 394

[Woo & Urry(2002)]Woo02
Woo, J.-H. & Urry, C. M. 2002, , 579, 530

[Worrall et al.(2012)Worrall, Birkinshaw, Young, Momtahan,
  Fosbury, Morganti, Tadhunter, & Verdoes Kleijn]Worr12
Worrall, D. M., Birkinshaw, M., Young, A. J., et al. 2012, , 424,
  1346

[Wright & Otrupcek(1990)]Wrig90
Wright, A. & Otrupcek, R. 1990, PKS Catalog (1990), 0

[Zycki & Czerny(1994)]Zyck94
Zycki, P. T. & Czerny, B. 1994, , 266, 653

[Zycki et al.(1999)Zycki, Done, & Smith]Zyck99
Zycki, P. T., Done, C., & Smith, D. A. 1999, , 305, 231





