X-Ray Weak AGNs from Super-Eddington Accretion onto Infant Black Holes
Abstract
A simple model for the X-ray weakness of JWST-selected broad-line AGNs is proposed under the assumption that the majority of these sources are fed at super-Eddington accretion rates. In these conditions, the hot inner corona above the geometrically thin disk that is responsible for the emission of X-rays in “normal” AGNs will be embedded instead in a funnel-like reflection geometry. The coronal plasma will Compton upscatter optical/UV photons from the underlying thick disk as well as the surrounding funnel walls, and the high soft-photon energy density will cool down the plasma to temperatures in the range 30–40 keV. The resulting X-ray spectra are predicted to be extremely soft, with power-law photon indices 2.8–4.0, making high- super-Eddington AGNs largely undetectable by Chandra.
1 Introduction
Deep surveys with the James Webb Space Telescope (JWST) have revealed an emergent, large population of moderate-luminosity, broad-line active galactic nuclei (AGNs) at powered by accretion onto early massive black holes (MBHs) (see, e.g., Kocevski et al., 2023; Harikane et al., 2023; Maiolino et al., 2023). These MBHs appear to be overmassive (but see Ananna et al., 2024) compared to their galaxy hosts and to grow while being extremely weak in X-rays (e.g. Maiolino et al., 2024; Yue et al., 2024). It has already been suggested that accretion at super-Eddington rates may answer many of the theoretical challenges posed by these and other observations. Some of the attractive features of this scenario include the very rapid grow of light seeds occurring during short-lived supercritical episodes (Madau et al., 2014; Volonteri et al., 2015; Pezzulli et al., 2016), as well as the artificial bias in black hole mass estimates induced by the anisotropic radiation field emitted by thick accretion tori (e.g. King, 2024; Lupi et al., 2024).
Supercritical accretion flows are known to form geometrically and optically thick disks dominated by radiation pressure. The large thickness of the disk naturally collimates radiation and produces a highly super-Eddington photon flux along the rotation axis (for a review, see Abramowicz & Fragile, 2013). In this regime, the inner, hot and luminous funnel region remains hidden from view, visible only at small viewing angles from the rotation axis of the system (e.g., Sikora, 1981; Madau, 1988; Sadowski et al., 2014; Jiang et al., 2014; Ogawa et al., 2017). Pacucci & Narayan (2024) have recently argued that the diminishing X-ray contribution predicted at high inclination angles from the pole may offer an explanation for the X-ray weakness of the “little red dots”, an enigmatic subcomponent of JWST-selected broad-line AGNs (e.g. Greene et al., 2024; Kokorev et al., 2024; Matthee et al., 2024).
In this Letter we propose a simple model for the X-ray weakness of JWST-selected AGNs. At super-Eddington rates, the hot corona above the inner accretion disks that is thought to be responsible for the emission of X-rays in “normal” AGNs (e.g., Haardt & Maraschi, 1991) is embedded in a funnel-like reflection geometry. The coronal plasma upscatters soft photons from the underlying disk as well as the surrounding funnel walls, and the high UV energy density cools down the plasma to temperatures below 40 keV. We shall find that the resulting X-ray spectra are extremely soft, with hard X-ray bolometric corrections that can be two orders of magnitude larger than those of standard AGNs.
2 Supercritical Accretion Flows
The shape and luminosity of low-viscosity, rotating, radiation-pressure supported accretion flows around black holes have been computed by several authors under a number of simplifying assumptions (e.g., Paczyńsky & Wiita, 1980; Wiita, 1982; Wielgus et al., 2016). These axisymmetric, , semi-analytical “thick disk” models differ from the standard “slim disk” solution (Abramowicz et al., 1988; Sadowski, 2009; Wang et al., 2014; Lasota et al., 2016), in that they are geometrically thick and radiatively efficient. We choose and briefly summarize below the thick disk formulation because numerical simulations of super-Eddington accretion flows have shown that vertical advection of radiation caused by magnetic buoyancy transports energy faster than radiative diffusion (Jiang et al., 2014, 2019). This effect allows photons to escape from the surface of the thick disk before being trapped and radially advected into the hole and undermines the underluminous – at a given – slim disk solution.
Let us adopt a cylindrical coordinate system ( centered on a Schwarzschild black hole of mass , and use a pseudo-Newtonian potential to mimic general relativistic effects, , where is the spherical radius and is the Schwarzschild radius. A physically realistic specific angular momentum distribution of the form (Paczyńsky & Wiita, 1980; Wiita, 1982)
(1) |
completely determines the shape of a thick disk (in the case of a polytropic gas depends only on the position coordinate ). The disk half-thickness can be shown to be given by
(2) |
where . The intersections of with the Keplerian angular momentum distribution provides the inner edge of the torus (), a pressure maximum (), and an outer transition radius () where the thick solution matches on to a thin accretion disk. The inner radius lies between the marginally bound and marginally stable orbits; its location determines the efficiency of the conversion of accreted matter into radiation.
Radiation is emitted from the photosphere at the local Eddington rate
(3) |
where is the effective gravity vector perpendicular to the surface of the disk, and is the electron scattering opacity. The total luminosity radiated by the torus, , is calculated by integrating the emitted flux over the two disk surfaces,
(4) |
where is the disk area element. In the case of large tori, the total luminosity generated by viscosity is simply given by
(5) |
where is the specific binding energy at the inner edge. Global energy conservation requires that the total energy gain be compensated by radiative losses, i.e. . By equating Equations (4) and (5) we then finally calculate .
\topruleParameter | Values for Models A and B |
---|---|
2.2, 2.3 | |
0.04792, 0.063301 | |
320, 163 | |
0.035, 0.044 | |
4.3, 3.1 | |
12.5, 7.1 | |
35∘, 44∘ |
The properties of two representative super-Eddington disk models are summarized in Table 1. Values for the disk inner edge and the constant in Equation (1) were chosen following Paczyńsky & Wiita (1980). As decreases, the disk grows larger and fatter (narrower funnel) and the accretion flow becomes more super-Eddington. The locally-generated radiation flux peaks at a few Schwarzschild radii and then drops as in the inner funnel. The total isotropic luminosity radiated in Models A and B is and 3.1, corresponding to and 7.1, respectively. Here, we have defined the critical accretion rate, , for powering the Eddington luminosity assuming a 10% radiative efficiency. The shape of the thick portion of the disks (one quadrant only) is outlined in Figure 1, while Figure 2 depicts the narrow, low-density funnel that develops in the innermost regions. This can be characterized by the opening half-angle :
(6) |
Soft photons emitted within the funnel will scatter off its walls many times before escaping to infinity.
3 Coronal X-ray Emission
The X-ray emission of AGNs is thought to be produced via thermal Comptonization of UV photons by hot electrons in a corona (e.g., Haardt & Maraschi, 1991). The soft seed photons originates in a cold, geometrically thin accretion disk. The geometry of the hot X-ray corona is poorly known, but a variety of observational constraints point to a very compact X-ray source located within Schwarzschild radii of the black hole (e.g., Reis & Miller, 2013; Fabian et al., 2015). In general, the hard X-ray spectra of Seyfert galaxies can be well approximated by a power law with photon index and a high-energy cutoff. These depend on the physical properties of the scattering corona, i.e., the optical depth and the electron temperature (e.g., Tortosa et al., 2018).
To extend this two-phase disk-corona scenario to supercritical flows, we place an artificial, plane-parallel, thin corona just above the funnel inner surface, and denote with the fraction of the locally generated radiation power that emerges from this hot tenuous phase. Each area element within the funnel will “see” incoming radiation that must be included in the balance of forces. Equation (3) must then be modified as:
(7) |
where now is the sum of the flux locally generated (‘self-flux’) and the normal component of the incoming flux .
Consider first, for illustration, a geometrically thin disk scenario with , and let us write the total Compton luminosity per unit area in the corona as the sum of an upward-, , and downward-directed, , component. Ignoring the anisotropy of inverse Compton scattering, we can express the power radiated in all directions by the hot phase as , where is the mean energy gain per scattering, is the total soft photon flux emerging from the cool disk, and is the mean probability of scattering in the corona. The outgoing photon flux is . Hard X-ray photons emitted by the corona and directed downward toward the cool, optically thick disk layers are largely absorbed except for a small reflected component, (Lightman & White, 1988). Under the assumption of a nonreflective disk, energy balance for the cold and hot phases then gives
(8) | ||||
Solving for and using the previous definitions, we derive
(9) |
In the regime of unsaturated Comptonization, the mean energy gain per scattering is given by , where is the dimensionless electron temperature. Together with Equation (9), this implies a relation between the luminosity of the corona and its scattering optical depth and temperature, and therefore the photon index of the power-law Comptonized component (Haardt & Maraschi, 1993). The Comptonization radiative transfer models of Sunyaev & Titarchuk (1980) give
(10) |
According to the above expressions, spectral indeces and electron temperatures are in the range and keV for when one third of the soft photospheric emission of the disk results from the reprocessing of the hard X-ray photons emitted by the corona, i.e. when . Such hard X-ray photon indeces are seen above 2 keV in low-redshift radio-quiet AGNs (; Piconcelli et al., 2005). The highest luminosity AGNs at have X-ray spectra that are steeper than these values (; Zappacosta et al., 2023), and so do near-Eddington narrow-line Seyfert 1 (NLS1) galaxies (; Grupe et al., 2010). Note that the thin disk-corona scenario is highly anisotropic, as soft photons enter the corona only from below.
Consider instead the case of a thick supercritical disk. Near the bottom of the narrow funnel, where much of the luminosity is generated, the soft radiation field has a large isotropic component. The coronal plasma will Compton cool not just on the locally-generated soft photons coming from below, but also on those emitted by the surrounding funnel walls and entering the corona from above. The temperature of the hot electrons will drop, and the emitted X-ray spectrum will be much softer than derived previously. This can be seen by rewriting Equations (8) as
(11) | ||||
where the third term is the incoming external flux expressed as a multiple of the self-flux times a factor to account for the fraction that scatters off the corona. It has long been recognized that in narrow accretion funnels (small half-angle ) the self-flux is a much smaller quantity than either or because of the strong “reflection effect” (Sikora, 1981; Narayan et al., 1983; Madau, 1988). Solving for the amplification factor as before we now have
(12) |
Figure 3 demonstrates how, with increasing values of the paramater , AGN X-ray spectra become extremely soft as the corona cools down. Our scenario is schematically illustrated in Figure 4.
To estimate , we have subdivided the surface of the funnel into a finite number of -rings with quasi-uniform emission properties, and approximated the force balance Equation (7) as
(13) |
where
(14) |
is the fraction of radiation emitted from the th ring that reaches the th ring along the direction in the solid angle given by the term in square brackets. Here, is the outward unit normal to the surface at the element , and is the mesh size. Solving this set of linear equations for Models A and B gives and at , respectively. Models with higher supercritical rates (smaller ) yield values of near the bottom of the funnel.
Compare now a “standard” AGN corona in the open geometry of a thin accretion disk. Assuming and , the two-phase thin disk-corona model produces a cutoff power-law X-ray spectrum with and keV. The derived rest-frame keV bolometric correction, (for an assumed UV soft photon input at 100 eV), is comparable to the value observed in low-luminosity, “normal” type 1 AGNs (Duras et al., 2020). In the case of a corona deep in the funnel of a supercritical disk with and the same values of and , we predict instead , keV, and a bolometric correction of . We observe here that comparable values of are actually measured at low redshifts in near-Eddington NLS1 galaxies (Vasudevan & Fabian, 2007). Note also the steep dependence of the X-ray bolometric correction with the strength of the isotropic soft photon field: for narrower funnels (higher supercritical accretion rates) with , our model yields , keV, and . For comparison, the X-ray stacking of JWST-selected, faint type 1 AGNs at gives (Maiolino et al., 2024).
4 Discussion
The vast majority of AGNs identified by JWST, both type 1 and type 2, are not detected even with the deepest Chandra observations, not even in stacking (e.g., Maiolino et al., 2024; Yue et al., 2024). We have discussed the possibility that these AGNs may have intrinsically faint X-ray emission because they are fed at super-Eddington rates. Specifically, we have shown that in this case their inner hot coronas will be embedded in a funnel-like reflection geometry where most of the luminosity is generated. The nearly isotropic soft photon field Compton will cool the coronal plasma to much lower temperatures than in the standard open geometry of a thin accretion disk, and the emitted X-ray spectrum will be extremely soft.
In our calculations several approximations have been adopted. Radiation-pressure driven winds are expected to naturally arise from the innermost radii of supercritical disks, carrying away mechanical energy and modifying the emission properties of the funnel regions. The accretion flows discussed here, however, are only mildly super-Eddington, with in the range 7–12.5. In this regime, numerical simulations of accretion flows around non-spinning MBHs show that inner hot coronae form because of the dissipation of buoyantly rising magnetic fields above the photospheres of thick disks, funnels remain optically thin, and the kinetic luminosity of the outflows is only of the radiation luminosity (Jiang et al., 2014, 2019). It is these and similar simulation results that have inspired the analytical study presented here. Nevertheless, while thick disk models are constructed to have normal forces at the surface in balance (see Eq. 7), the tangential components of the absorbed radiation flux are left unbalanced, and one should account for the fact that the funnel surface layers and the corona above super-Eddington flows may not be static but moving upward (Sikora & Wilson, 1981; Narayan et al., 1983). Bulk Compton scattering could then produce an additional high-energy radiation term. Eventually, more physical Comptonization models that include an X-ray reflection component from neutral disk material should be used for detailed comparison with the observations.
In a recent work, Pacucci & Narayan (2024) have pointed out that super-Eddington accretion would also hide coronal X-ray emission from view at high inclination angles from the poles, and this effect may offer a complementary explanation for the X-ray weakness of many JWST-selected AGNs. Ultimately, our model of intrinsically faint X-ray emission in supercritical flows may not be able to explain objects like the “dormant” MBH GN-1146115 at (Juodžbalis et al., 2024), which seems to be accreting well below the Eddington limit and is undetected in X-rays with (Maiolino et al., 2024). Yet, the properties of this faint AGN with a high black hole-to-stellar mass ratio may support a scenario in which short bursts of super-Eddington accretion (Madau et al., 2014) have resulted in black hole overgrowth. During one of these bursts, the MBH in GN-1146115 was likely fed by a thick, radiation-supported supercritical torus like those discussed in this paper.
Acknowledgements
Support for this work was provided by NASA through grant TCAN 80NSSC21K027, by grant NSF PHY-2309135 to the Kavli Institute for Theoretical Physics (KITP), and by the Italian Ministry for Research and University (MUR) under Grant ‘Progetto Dipartimenti di Eccellenza 2023-2027’ (BiCoQ). One of us (PM) acknowledges useful discussions and inputs on this project with O. Blaes, A. Ferrara, R. Maiolino, and F. Pacucci. He is also grateful to V. Bromm, B. Robertson, R. Schneider, and R. Somerville for organizing the KITP workshop “Cosmic Origins: The First Billion Years”, which provided the initial motivation for this work.
References
- Abramowicz et al. (1988) Abramowicz, M. A., Czerny, B., Lasota, J. P., & Szuszkiewicz, E. 1988, ApJ, 332, 646
- Abramowicz & Fragile (2013) Abramowicz, M. A., & Fragile, P. C. 2013, Living Reviews in Relativity, 16, 1
- Ananna et al. (2024) Ananna, T. T., Bogdán, Á., Kovács, O. E., Natarajan, P., & Hickox, R. C. 2024, ApJ, 969, L18
- Duras et al. (2020) Duras, F., Bongiorno, A., Ricci, F., et al. 2020, A&A, 636, A73
- Fabian et al. (2015) Fabian, A. C., Lohfink, A., Kara, E., et al. 2015, MNRAS, 451, 4375
- Greene et al. (2024) Greene, J. E., Labbe, I., Goulding, A. D., et al. 2024, ApJ, 964, 39
- Grupe et al. (2010) Grupe, D., Komossa, S., Leighly, K. M., & Page, K. L. 2010, ApJS, 187, 64
- Haardt & Maraschi (1991) Haardt, F., & Maraschi, L. 1991, ApJ, 380, L51
- Haardt & Maraschi (1993) —. 1993, ApJ, 413, 507
- Harikane et al. (2023) Harikane, Y., Zhang, Y., Nakajima, K., et al. 2023, ApJ, 959, 39
- Jiang et al. (2014) Jiang, Y.-F., Stone, J. M., & Davis, S. W. 2014, ApJ, 796, 106
- Jiang et al. (2019) —. 2019, ApJ, 880, 67
- Juodžbalis et al. (2024) Juodžbalis, I., Maiolino, R., Baker, W. M., et al. 2024, arXiv e-prints, arXiv:2403.03872
- King (2024) King, A. 2024, MNRAS, 531, 550
- Kocevski et al. (2023) Kocevski, D. D., Onoue, M., Inayoshi, K., et al. 2023, ApJ, 954, L4
- Kokorev et al. (2024) Kokorev, V., Caputi, K. I., Greene, J. E., et al. 2024, ApJ, 968, 38
- Lasota et al. (2016) Lasota, J. P., Vieira, R. S. S., Sadowski, A., Narayan, R., & Abramowicz, M. A. 2016, A&A, 587, A13
- Lightman & White (1988) Lightman, A. P., & White, T. R. 1988, ApJ, 335, 57
- Lupi et al. (2024) Lupi, A., Trinca, A., Volonteri, M., Dotti, M., & Mazzucchelli, C. 2024, A&A, 689, A128
- Madau (1988) Madau, P. 1988, ApJ, 327, 116
- Madau et al. (2014) Madau, P., Haardt, F., & Dotti, M. 2014, ApJ, 784, L38
- Maiolino et al. (2023) Maiolino, R., Scholtz, J., Curtis-Lake, E., et al. 2023, arXiv e-prints, arXiv:2308.01230
- Maiolino et al. (2024) Maiolino, R., Risaliti, G., Signorini, M., et al. 2024, arXiv e-prints, arXiv:2405.00504
- Matthee et al. (2024) Matthee, J., Naidu, R. P., Brammer, G., et al. 2024, ApJ, 963, 129
- Narayan et al. (1983) Narayan, R., Nityananda, R., & Wiita, P. J. 1983, MNRAS, 205, 1103
- Ogawa et al. (2017) Ogawa, T., Mineshige, S., Kawashima, T., Ohsuga, K., & Hashizume, K. 2017, PASJ, 69, 33
- Pacucci & Narayan (2024) Pacucci, F., & Narayan, R. 2024, arXiv e-prints, arXiv:2407.15915
- Paczyńsky & Wiita (1980) Paczyńsky, B., & Wiita, P. J. 1980, A&A, 88, 23
- Pezzulli et al. (2016) Pezzulli, E., Valiante, R., & Schneider, R. 2016, MNRAS, 458, 3047
- Piconcelli et al. (2005) Piconcelli, E., Jimenez-Bailón, E., Guainazzi, M., et al. 2005, A&A, 432, 15
- Reis & Miller (2013) Reis, R. C., & Miller, J. M. 2013, ApJ, 769, L7
- Sadowski (2009) Sadowski, A. 2009, ApJS, 183, 171
- Sadowski et al. (2014) Sadowski, A., Narayan, R., McKinney, J. C., & Tchekhovskoy, A. 2014, MNRAS, 439, 503
- Sikora (1981) Sikora, M. 1981, MNRAS, 196, 257
- Sikora & Wilson (1981) Sikora, M., & Wilson, D. B. 1981, MNRAS, 197, 529
- Sunyaev & Titarchuk (1980) Sunyaev, R. A., & Titarchuk, L. G. 1980, A&A, 86, 121
- Tortosa et al. (2018) Tortosa, A., Bianchi, S., Marinucci, A., Matt, G., & Petrucci, P. O. 2018, A&A, 614, A37
- Vasudevan & Fabian (2007) Vasudevan, R. V., & Fabian, A. C. 2007, MNRAS, 381, 1235
- Volonteri et al. (2015) Volonteri, M., Silk, J., & Dubus, G. 2015, ApJ, 804, 148
- Wang et al. (2014) Wang, J.-M., Qiu, J., Du, P., & Ho, L. C. 2014, ApJ, 797, 65
- Wielgus et al. (2016) Wielgus, M., Yan, W., Lasota, J. P., & Abramowicz, M. A. 2016, A&A, 587, A38
- Wiita (1982) Wiita, P. J. 1982, ApJ, 256, 666
- Yue et al. (2024) Yue, M., Eilers, A.-C., Ananna, T. T., et al. 2024, arXiv e-prints, arXiv:2404.13290
- Zappacosta et al. (2023) Zappacosta, L., Piconcelli, E., Fiore, F., et al. 2023, A&A, 678, A201