Abstract
Estimates suggest that while FRII jets appear to have lifetimes constrained to hundreds of millions of years, radio galaxies with FRI jets appear to be longer lived. We illustrate the nature of this time constraint from model perspectives, showing how compatibility between theory and data match in a way suggesting a key difference between active galaxies whose engines are characterized by accretion onto co-rotating versus counter-rotating black holes. We calculate a range of timescales for counter-rotating black holes for a range of accretion rates compatible with theory which we then compare to data. The validity of these timescales constitutes the most powerful recent piece of evidence for considering counter-rotation between black holes and accretion disks in high energy astrophysics.
1 Introduction
According to the current paradigm, powerful FRII jets in radio galaxies are the product of accretion onto rapidly spinning, prograde accreting black holes (; ; ; ). The FRII classification refers to jets that are more collimated and generally more powerful than FRI jets, the latter often being subjected to entrainment from the interstellar medium (). Because prograde accretion onto a black hole can only spin black holes up, the only constraint for the longevity of the jet is the amount of accreting fuel. Radio galaxies can either be high excitation or low excitation, depending on the degree of emission line signatures or thermal versus non-thermal nature (; ; ; ; ; ). Because FRII radio galaxies are often low excitation systems, they should experience constraints on their lifetimes that are similar to those for FRI radio galaxies, such as M87, with timescales orders of magnitude longer than those associated with feeding at near Eddington rates. But this is not supported by the data. In fact, radio galaxies with FRII jets appear to have quantifiably limited timescales unlike their FRI counterparts (e.g., ; ). High excitation FRII systems, for example, are found to be limited to 10 million years (). Most recently, have compared radio galaxies with the giant radio galaxy population, including the lifetimes of jets in FRII sources, among others. This constitutes the most exhaustive quantitative analysis of FRII lifetimes and, if these results hold up to future scrutiny, a powerful constraint on the nature of jet formation and evolution in jetted active galactic nuclei (AGN). We suggest that the difference in measured timescales for powerful FRII jets compared to powerful FRI jets points to a basic difference in the nature of the two morphologies that was captured in the gap paradigm for black hole accretion and jet formation (). Whereas powerful FRII jets, in this paradigm, are produced in accreting black holes spinning in the opposite direction as the accretion disk (i.e., counter-rotation), the opposite is true for FRI jets. And since counter-rotation spins black holes down while co-rotation spins them up indefinitely, powerful FRII jets are limited in time in a way that powerful FRI jets are not. The possibility that radio galaxies with FRII jet morphology are constrained in time unlike FRI radio galaxies is, therefore, interesting in a fundamental way in high energy astrophysics. Evidence that FRII radio galaxies are constrained in time in a way that matches model predictions for the spin down timescales in both high and low excitation systems is thus exciting for understanding the nature of the longstanding puzzle behind the FRI/FRII jet dichotomy. In Section 2 we discuss the data analyzed in , describe the theory and emphasize the match between theory and data in Section 3. In Section 4 we conclude.
2 Data
explored the time evolution of 361 FRII radio galaxies from Cambridge, 3CRR, 6CE, 5C6, and 5C7 Sky Surveys and from the Bologna B2, Green Bank GB, and GB2 Surveys in order to produce a statistically relevant sample. They obtained a range of lifetimes for FRII radio galaxies which plot in their Figure 6 and which we show on the left side of Figure 1. On the right-hand side of Figure 1 we show the maximum lifetime data from the left-hand side of Figure 1 on the dynamical age of the FRII jet with jet length obtained from . In other words, we hone in on the 4 objects with oldest dynamical ages (circles) for each class of object, namely, radio galaxies, radio quasars, and giant radio galaxies and giant radio quasars. The red objects represent high excitation FRII jetted AGN, i.e., with quasar or thermal-like signatures, indicative of radiatively efficient accretion. The red circle on the left side of the vertical line has the maximum dynamical age for an FRII high excitation radio galaxy or FRII HERG. The red object on the right side of the vertical line represents a giant radio quasar. I.e., it has all the same characteristics as the red counterpart on the left except for its jet length. The object on the right is considered a giant FRII HERG. The blue objects, similarly, distinguish themselves in the same way as the red objects do, except they belong to radiatively inefficient accretion, showing an absence of thermal or emission line signatures. They are thus FRII LERG for low excitation radio galaxies. As described in Section 3, we add theoretical values to Figure 1 with diamonds associating their theoretical age (i.e., the model prescribed duration of time for the FRII jet) with the same jet length values as the observational data for ease of comparison.
FIGURE 1
The maximum lifetime observed so far for FRII HERG and FRII LERG circled in red and blue on the left panel of Figure 1, with values equal to 9 × 107 and under 3 × 108 years, respectively, are compatible with estimates from theoretical modeling, which we will show in Section 3. While the lifetimes for giant radio quasars and giant radio galaxies increase from a theoretical perspective (), their value is more uncertain because such objects are not necessarily triggered at the Eddington accretion limit. Because of this uncertainty, we do not include theoretical maximum lifetimes for giant radio quasars and giant radio galaxies. Note that if a radio quasar had a lifetime equal to 230 million years (as the giant radio quasar indicated with the red circle), it would violate theory. The theoretical limit, in fact, is shown as the red diamond.
For FRI radio galaxies, as mentioned above, no such time constraint is found. Two decades ago, evidence began to emerge suggesting that FRI jets live about an order of magnitude longer than FRII jets (). The evidence in this respect has grown (e.g., ) and it was found that LERG systems live up to order 109 years (). The constraints on FRII jets from suggest, therefore, that long-lived LERG systems are FRI. In other words, an FRI radio galaxy can live substantially longer than FRII sources. Why is that?
3 Theory
In this section we describe how the timescales above emerge from, or are compatible with, theory. The theory is anchored to the idea that counter-rotation between black holes and accretion disks give rise to FRII jets (). As counter-rotation is an unstable and less likely configuration, the majority of mergers in this paradigm funnel cold gas into the nucleus that settles into co-rotation around a spinning black hole. It is therefore the minority that end up in counter-rotation and this is environment-dependent. It is also worth noting that an engine-based difference between FRII and FRI matters for the most powerful jets and that environment makes a difference at lower jet power (see on the nature of the Owen-Ledlow diagram). We will therefore focus on the most powerful jets in the paradigm.
For the minority of configurations that end up in counter-rotating accretion disks around spinning black holes, we explore two basic evolutionary scenarios relevant for understanding the maximum possible lifetimes of FRII jets. The question we need to answer is this: How long does it take to spin a black hole down? This is because FRII jets are associated with counter-rotation. We want to find the maximum possible time for this process to then compare with the data for FRII lifetimes. Mergers yield initial conditions involving cold gas funneled into the galactic nucleus and the formation of a radiatively efficient disk accreting in counter-rotation at the Eddington limit. A subset of these counter-rotating black holes spin down at the Eddington limit while others spin down at accretion rates that begin at the Eddington limit but drop to rates as low as 10–2 the Eddington accretion rate. This range of accretion rates determines the range of jet lifetimes. We should point out that the initial spin value is crucial in determining jet lifetime. But since we are striving to determine maximum jet lifetimes, the initial spin value is assumed to be its theoretical maximum at 0.998. The drop in accretion rate depends on feedback processes, the details of which are not of present concern (see for details). What matters here is the range of accretion for counter-rotation because it determines the lifetime for a counter-rotating black hole. This timescale depends on the amount of angular momentum added to the black hole by the accreted plasma and the rate at which the plasma accretes. The angular momentum accreted is that of the gas at the inner edge of the disk, which is referred to as the marginally stable circular orbit, rms. This radial location depends both on the spin of the black hole and the orientation of the accretion disk. From the stability of circular orbits in Kerr space-time one finds that rms drops from 9 gravitational radii, rg, for a maximally spinning black hole surrounded by a counter-rotating accretion disk, to just over 1 rg for a maximally spinning black hole surrounded by a co-rotating accretion disk (e.g., ). To determine the amount of angular momentum delivered to the black hole, one sets r to rms in the following expressions for the angular momentum per unit mass as a function of Boyer-Lindquist radial coordinate r (),andwhere the ‘+’ and ‘-’ subscripts refer to the value of the angular momentum per unit mass as a function of radial coordinate for co-rotating disks and counter-rotating disks, respectively, with a the spin of the black hole and M its mass. By multiplying by accreted mass, one obtains the amount of angular momentum accreted, which changes the spin of the black hole. Because accretion rates are model-prescribed in terms of the Eddington value, our results are scale invariant, making the actual values of angular momentum irrelevant for our calculations. To see this scale invariance explicitly, we begin with the definition of the dimensionless black hole spinwhere ℒ is the angular momentum of the black hole and c is the speed of light. If accretion proceeds at the Eddington limit we set the Eddington luminosity to the luminosity in terms of the accretion rate and the disk efficiency as in Eq. (4). For our purposes we repackage all the constants into one term and relate the accretion rate to the black hole mass in Eq. (5), which can be written as in Eq. (6). The infinitesimal mass accreted onto the black hole is therefore obtained from Eq. (6) as shown in Eq. (7). We can carry out the analysis in the Newtonian limit with the magnitude of the angular momentum of the infinitesimal parcel of gas given to the black hole given in Eq. (8), with v the velocity of the parcel of gas dm at the inner edge of the disk that is supplied to the black hole and r is its radial location. From circular motion and Newton’s 2nd law we obtain Eq. (9) from which we obtain the velocity in Eq. (10).
The inner edge of the disk depends on the black hole spin parameter and has the range given in Eq. (11). For our purposes we note that r ∝ GMBH/c2 from which we get Eq. (12). For accretion at the Eddington limit, therefore, the rate at which the angular momentum of the black hole changes is shown in Eq. (13), from which we can determine how the dimensionless spin parameter of the black hole changes by using the differential form of Eq. (3) to obtain Eq. (14), from which we get Eq. (15), which is the black hole mass independent result we anticipated.
In Figure 2 we show how the angular momentum of gas that accretes onto the black hole from the marginally stable orbit depends on the value of black hole spin. We scale or normalize the angular momentum to the angular momentum at the marginally stable orbit for a black hole spinning at a = 0.998 surrounded by an accretion disk in counter-rotation. Figure 2 allows one to appreciate why spinning a high spinning black hole down takes about an order of magnitude less time than it does to spin a zero spinning black hole up to high spin, at a given accretion rate.
FIGURE 2
As gas accretes onto the black hole from rms, both black hole spin and black hole mass change. The change in the black hole mass also depends on the location of rms and can be obtained by evaluating the distribution of energy as a function of radius, i.e., energy counterparts to Eqs (1), (2) above. One finds the black hole to gain an amount of mass given by ()where m is the mass of the black hole and the added mass ∆m is given as an expression with both Newton’s constant and the speed of light equal to unity. As accretion proceeds, the marginally stable orbit decreases, acquires the value rms = 6rg when the black hole stops rotating, and then decreases further as it spins up via a co-rotating accretion disk. Our interest is only in the time to spin the black hole down to near zero spin. The crucial element to determine the timescale for spin down is the accretion rate. For a given accretion rate f, the time to build the mass by ∆m is given in Eq. (18). If the accretion rate is constant, the time is given by Eq. (19).
For an accretion rate that is the Eddington value, a rapidly spinning counter-rotating black hole spins down to zero spin in just under 8 × 106 years. Therefore, an FRII jet lives no longer than 8 × 106 years if fed at the Eddington limit. But an FRII HERG does not need to be accreting at the Eddington limit. It could accrete at 10% the Eddington limit and still be a HERG. In this case, it would spin down to zero spin in 8 × 107 years. For lower accretion rates, the timescale for spin down increases. If the FRII jet is powerful enough, it produces a strong feedback effect on the accretion flow, lowering the accretion rate, and allowing the FRII jet phase to last longer. In such cases of powerful jet feedback, the FRII jet affects the structure of the accretion disk, forcing it to evolve into an advection dominated accretion flow (ADAF). The boundary between a thin disk and an ADAF is prescribed from theory to be at 10–2 the Eddington accretion rate. Hence, at 5 × 10−2 the Eddington accretion rate, the object may still be characterized as an FRII HERG and the jet lifetime would increase to 1.6 × 108 years. We should also note that jet lifetimes are effectively limited by some threshold low spin value below which the jet may no longer be classified as an FRII if even visible. From theory, we can estimate this to be below a spin value of about 0.1 but with some uncertainty that would also depend on the environment. Overall, we can estimate that an FRII HERG lives at most about 108 years. This value appears as a red diamond on the right hand side of Figure 1.
Since the transition in cooling is not abrupt (e.g., ), the transition from cold mode accretion into an ADAF is also gradual, and the model prescribes accretion rates to barely cross the boundary into an ADAF during counter-rotation. In other words, counter-rotation may have ADAF accretion, but it has only recently entered the ADAF regime and accretion rates will therefore have values near 10–2 the Eddington accretion rate (see for details). If one assumes such a value for the accretion rate in Eq. (17), one obtains a timescale for spin down of 8 × 108 years. But, as mentioned above, the initial state is a near-Eddington accreting black hole so the timescale prescribed in the model for systems that evolve into ADAF accretion, must be lower. In short, the black hole spins halfway down from maximal spin at near the Eddington accretion rate (as an FRII HERG) which takes 4 × 106 years, followed by accretion at 10–2 the Eddington accretion rate (as an FRII LERG), requiring 4 × 108 more years to spin the black hole down to zero spin. Therefore, the model prescription for FRII lifetimes spans the range 8 × 106 to 4 × 106 + 4 × 108 years which is roughly 4 × 108 years and is our second theoretical data point in Figure 1 (the blue diamond). Because these lifetimes occur in systems with strongest FRII jet feedback, and the model prescribes the strongest jet feedback to occur more in denser environments, there is a model-prescribed environment-dependence to the maximum FRII jet lifetimes that is worth mentioning although it is not the focus of this work.
Although our focus has not been on giant radio quasars and giant radio galaxies, such objects serve an important role as guideposts, allowing us to understand better the constraints on radio quasars and radio galaxies. As re-triggered counter-rotating black holes, giant radio quasars have the opportunity to generate jets that extend beyond the kiloparsec lengths reached by their radio galaxy ancestors, and to experience longer lived FRII phases (). Accordingly, it is interesting to note that while giant radio quasars and giant radio galaxies exceed the prescribed theoretical maximum lifetimes for radio quasars (the red diamond in Figure 1), this is not true for radio quasars. If counter-rotation is not relevant to FRII jets, there is no reason for the latter to be limited in this way. In addition, if the evolution in time from quasar mode to ADAF mode (as the model prescribes) is not the way such systems change over time, there is no reason for both regular radio quasars and radio galaxies as well as their giant counterparts to have blue objects experiencing longer jet lifetimes than red ones. But it is instead required from theory.
4 Conclusion
The timescales for FRII systems obtained from theory are tantalizingly compatible with those inferred from the data as seen from the diamonds added on the right-hand side of Figure 1 to represent theoretical timescales. While spin-down timescales are constrained by the rate of accretion, this same constraint on time for FRI jets in the theory is rather weak because a black hole that feeds forever simply remains a high spinning co-rotating black hole. The real constraint, instead, is the amount of fuel. We have not gone into model details but FRI systems are late stages in the evolution of radio galaxies that were once FRII. Their accretion rates continue to drop over time and can be orders of magnitude lower than the Eddington accretion rate. As a result, FRI systems accreting in ADAF can last characteristic timescales that are on the order of the age of the Universe, making them effectively unconstrained in time. In closing, we highlight that jet dynamical lifetimes are characterized by large uncertainties (e.g., ) and that until these have been sufficiently reduced, a clear picture cannot emerge and caution should be exercised in comparing jet lifetime from the model with dynamical timescales from data. Nonetheless, if FRII lifetimes can robustly be shown to be limited to within a half billion years, it would constitute strong evidence for counter-rotating black holes in active galaxies and the evidence appears to be pointing in that direction.
While constraints on the lifetimes of FRII jetted AGN have existed for a decade or so, data has recently emerged to solidify the case that FRII and FRI jets are different in some fundamental way. We have argued over the last decade that opening the counter-rotating window for black hole accretion allows many disparate observations to come together under a simple evolutionary picture that at its heart explains the radio loud/radio quiet dichotomy. In this work we highlight the otherwise coincidental match between the lifetimes of FRII jets in quasars and radio galaxies, showing how to understand the difference in the evolution of FRII jets compared to FRI jets. Unlike FRI LERG whose jet lifetimes are effectively unconstrained, FRII jets in either LERG or HERG form, are limited to lifetimes within hundreds of millions of years due to accretion spinning black holes down, a process that is time-limited in a way that spinning black holes up is not.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
The author confirms being the sole contributor of this work and has approved it for publication.
Acknowledgments
I thank Dr. Marek Jamrozy and Dr. Dorota Koziel-Wierzbowska for sharing their expertise. In addition, I acknowledge the role of 4 referees at FrASS but thank referees 3 and 4 for pointing to the need for clarification on key points. The reason FRIIs appear to prefer less dense environments compared to FRIs is that they live longer in such environments. This resolves an interesting point raised by referee 4 that was not included in the paper because it is outside its scope. The issue is discussed in our work on X-shaped radio galaxies in 2020.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AntonucciR. (2012). Gas-phase metallicity of 27 galaxies at intermediate redshift. Astronomy Astrophysics27, 557. 10.48550/arXiv.1210.2716
2
BardeenJ. M.PressW. H.TeukolskyS. A. (1972). Rotating black holes: Locally nonrotating frames, energy extraction, and scalar synchrotron radiation. astrophysical J.178, 347. 10.1086/151796
3
BestP. N.HeckmanT. M. (2012). On the fundamental dichotomy in the local radio-AGN population: Accretion, evolution and host galaxy properties: The fundamental dichotomy in local radio-AGN. Mon Notices R Astronomical Soc. 421, 1569–1582. 10.1111/j.1365-2966.2012.20414.x
4
DabhadeP.SaikiaD. J.MahatoM. (2023). Decoding the giant extragalactic radio sources. J. Astrophys. Astr. 44, 13. 10.1007/s12036-022-09898-5
5
FanaroffB. L.RileyJ. M. (1974). The morphology of extragalactic radio sources of high and low luminosity. Mon. Notices R. Astronomical Soc. 167, 31–36P. 10.1093/mnras/167.1.31p
6
GarofaloD. (2022). Front. Astron. Space Sci. 8, 258. 10.3389/fspas.2021.785092
7
GarofaloD.EvansD. A.SambrunaR. M. (2010). Mon. Notices R. Astronomical Soc.406, 975. 10.1111/j.1365-2966.2010.16797.x
8
GarofaloD.SinghC. B.ZackA. (2018). The distribution and lifetime of powerful radio galaxies as a function of environment and redshift. SciRep8, 15097. 10.1038/s41598-018-33532-6
9
GiustiniM.ProgaD. (2019). A global view of the inner accretion and ejection flow around super massive black holes: Radiation-driven accretion disk winds in a physical context. Astronomy Astrophysics630, A94. 10.1051/0004-6361/201833810
10
HineR. G.LongairM. S. (1979). Optical spectra of 3CR radio galaxies. Mon. Notices R. Astronomical Soc. 188, 111–130. 10.1093/mnras/188.1.111
11
MacconiD.TorresiE.GrandiP.BoccardiB.VignaliC. (2020). Radio morphology–accretion mode link in Fanaroff–Riley type II low-excitation radio galaxies. Mon. Notices R. Astronomical Soc. 493, 4355–4366. 10.1093/mnras/staa560
12
MachalskiJ.Koziel-WierzbowskaD.GoyalA. (2021). An atlas of dynamical evolution models of 361 Fanaroff–Riley type II radio sources. Astrophysical J. Suppl. Ser. 255, 22. 10.3847/1538-4365/ac08a0
13
McClintockJ. E.NarayanR.DavisS. W.GouL.KulkarniA.OroszJ. A.et al (2011). Fundamental physics and cosmology with LISA. CQG. 28, 4009. 10.1088/0264-9381/28/11/114009
14
MingoB.CrostonJ. H.BestP. N.DuncanK. J.HardcastleM. J.KondapallyR.et al (2022). Mon. Notices R. Astronomical Soc. 511, 3250. 10.1093/mnras/stac140
15
MingoB.HardcastleM. J.CrostonJ. H.DickenD.EvansD. A.MorgantiR.et al (2014). An X-ray survey of the 2 Jy sample – I Is there an accretion mode dichotomy in radio-loud AGN?Mon Notices R Astronomical Soc. 440, 269–297. 10.1093/mnras/stu263
16
O’DeaC. P.DalyR. A.KharbP.FreemanK. A.BaumS. A. (2009). Physical properties of very powerful FRII radio galaxies. Astronomy Astrophysics494, 471–488. 10.1051/0004-6361:200809416
17
ParmaP.MurgiaM.de RuiterH. R.FantiR. (2002). The lives of FR I radio galaxies. New Astr. Rev. 46, 313–325. 10.1016/s1387-6473(01)00201-9
18
RaineD.ThomasE. (2009). Deforestation causes 'boom-and-bust' development in the Amazon. 2. London, United Kingdom: Imperial College Press.
19
SaripalliL.SubrahmanyanR.ThoratK.EkersR. D.HunsteadR. W.JohnstonH. M.et al (2012). Atlbs extended source sample: The evolution in radio source morphology with flux density. astrophysical J. 199, 27. 10.1088/0067-0049/199/2/27
20
SikoraM.StawarzL.LasotaJ.-P. (2007). Radio loudness of active galactic nuclei: Observational facts and theoretical implications. astrophysical J. 658, 815–828. 10.1086/511972
21
TchekhovskoyA.NarayanR.McKinneyJ. C. (2010). Black hole spin and the radio loud/quiet dichotomy of active galactic nuclei. Astrophysical J. 711, 50–63. 10.1088/0004-637x/711/1/50
22
TchekhovskoyA.McKinneyJ. C. (2012). Prograde and Retrograde Black Holes: Whose Jet is More Powerful. Mon. Notices R. Astronomical Soc. 423, L55–L59. 10.1111/j.1745-3933.2012.01256.x
23
TurnerR. J.ShabalaS. S. (2015). Clash: The concentration-mass relation of galaxy clusters. astrophysical J. 806, 59. 10.1088/0004-637x/806/1/59
24
WilsonA. S.ColbertE. J. S. (1995). Toward a theory of interstellar turbulence. II. Strong alfvenic turbulence. astrophysical J. 438, 62. 10.1086/175054
25
WojtowiczA.StawarzJ.MachalskiJ.OstoreroL. (2021). A novel method for estimating the ambient medium density around distant radio sources from their observed radio spectra. astrophysical J. 922, 197. 10.3847/1538-4357/ac116c
Summary
Keywords
counter-rotating black holes, jets from active galaxies, FRI/FRII dichotomy, black holes, radio galaxies and quasars
Citation
Garofalo D (2023) Counter-rotating black holes from FRII lifetimes. Front. Astron. Space Sci. 10:1123209. doi: 10.3389/fspas.2023.1123209
Received
13 December 2022
Accepted
28 April 2023
Published
10 May 2023
Volume
10 - 2023
Edited by
Paola Marziani, Osservatorio Astronomico di Padova (INAF), Italy
Reviewed by
Enrico Congiu, European Southern Observatory, Chile
Dragana Ilic, University of Belgrade, Serbia
Updates
Copyright
© 2023 Garofalo.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: David Garofalo, dgarofal@kennesaw.edu
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.