Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

X-Ray Weak AGNs from Super-Eddington Accretion onto Infant Black Holes

Piero Madau Department of Astronomy & Astrophysics, University of California, 1156 High Street, Santa Cruz, CA 95064 Dipartimento di Fisica “G. Occhialini”, Università degli Studi di Milano-Bicocca, Piazza della Scienza 3, I-20126 Milano, Italy Francesco Haardt Dipartimento di Scienza e Alta Tecnologia, Università degli Studi dell’Insubria, via Valleggio 11, I-22100 Como, Italy INAF, Osservatorio Astronomico di Brera, Via E. Bianchi 46, I-23807 Merate, Italy INFN, Sezione Milano-Bicocca, P.za della Scienza 3, I-20126 Milano, Italy
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 Γsimilar-to-or-equalsΓabsent\Gamma\simeq\,roman_Γ ≃ 2.8–4.0, making high-z𝑧zitalic_z super-Eddington AGNs largely undetectable by Chandra.

Accretion (14); Active galactic nuclei (16); James Webb Space Telescope (2291); Supermassive black holes (1663)
\journalinfo

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 z=410𝑧410z=4-10italic_z = 4 - 10 powered by accretion onto M=106108M𝑀superscript106superscript108subscriptMdirect-productM=10^{6}-10^{8}\,\,{\rm M_{\odot}}italic_M = 10 start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT - 10 start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT roman_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT 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, M˙M˙Eddmuch-greater-than˙𝑀subscript˙𝑀Edd\dot{M}\gg\dot{M}_{\rm Edd}over˙ start_ARG italic_M end_ARG ≫ over˙ start_ARG italic_M end_ARG start_POSTSUBSCRIPT roman_Edd end_POSTSUBSCRIPT, 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 M˙˙𝑀\dot{M}over˙ start_ARG italic_M end_ARG – slim disk solution.

Let us adopt a cylindrical coordinate system (r,φ,z)r,\varphi,z)italic_r , italic_φ , italic_z ) centered on a Schwarzschild black hole of mass M𝑀Mitalic_M, and use a pseudo-Newtonian potential to mimic general relativistic effects, Φ=GM/(RrS)Φ𝐺𝑀𝑅subscript𝑟𝑆\Phi=-GM/(R-r_{S})roman_Φ = - italic_G italic_M / ( italic_R - italic_r start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ), where R=(r2+z2)1/2𝑅superscriptsuperscript𝑟2superscript𝑧212R=(r^{2}+z^{2})^{1/2}italic_R = ( italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT is the spherical radius and rS=2GM/c2subscript𝑟𝑆2𝐺𝑀superscript𝑐2r_{S}=2GM/c^{2}italic_r start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT = 2 italic_G italic_M / italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT is the Schwarzschild radius. A physically realistic specific angular momentum distribution (r)𝑟\ell(r)roman_ℓ ( italic_r ) of the form (Paczyńsky & Wiita, 1980; Wiita, 1982)

(r)=K(rin)+𝒞(rrin)𝑟subscript𝐾subscript𝑟in𝒞𝑟subscript𝑟in\ell(r)=\ell_{K}(r_{\rm in})+{\cal C}(r-r_{\rm in})roman_ℓ ( italic_r ) = roman_ℓ start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT ) + caligraphic_C ( italic_r - italic_r start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT ) (1)

completely determines the shape of a thick disk (in the case of a polytropic gas \ellroman_ℓ depends only on the position coordinate r𝑟ritalic_r). The disk half-thickness z=h(r)𝑧𝑟z=h(r)italic_z = italic_h ( italic_r ) can be shown to be given by

h(r)={[GM(rinrS)GM(rinrS)J(r)+rS]2r2}1/2,𝑟superscriptsuperscriptdelimited-[]𝐺𝑀subscript𝑟insubscript𝑟𝑆𝐺𝑀subscript𝑟insubscript𝑟𝑆𝐽𝑟subscript𝑟𝑆2superscript𝑟212h(r)=\left\{\left[{GM(r_{\rm in}-r_{S})\over GM-(r_{\rm in}-r_{S})J(r)}+r_{S}% \right]^{2}-r^{2}\right\}^{1/2},italic_h ( italic_r ) = { [ divide start_ARG italic_G italic_M ( italic_r start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT - italic_r start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ) end_ARG start_ARG italic_G italic_M - ( italic_r start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT - italic_r start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ) italic_J ( italic_r ) end_ARG + italic_r start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT } start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT , (2)

where J(r)rinr2(r)𝑑r/r3𝐽𝑟superscriptsubscriptsubscript𝑟in𝑟superscript2superscript𝑟differential-dsuperscript𝑟superscript𝑟3J(r)\equiv\int_{r_{\rm in}}^{r}\ell^{2}(r^{\prime})dr^{\prime}/r^{\prime 3}italic_J ( italic_r ) ≡ ∫ start_POSTSUBSCRIPT italic_r start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_r end_POSTSUPERSCRIPT roman_ℓ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) italic_d italic_r start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT / italic_r start_POSTSUPERSCRIPT ′ 3 end_POSTSUPERSCRIPT. The intersections of (r)𝑟\ell(r)roman_ℓ ( italic_r ) with the Keplerian angular momentum distribution K(r)subscript𝐾𝑟\ell_{K}(r)roman_ℓ start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT ( italic_r ) provides the inner edge of the torus (rinsubscript𝑟inr_{\rm in}italic_r start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT), a pressure maximum (rcsubscript𝑟𝑐r_{c}italic_r start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT), and an outer transition radius (routsubscript𝑟outr_{\rm out}italic_r start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT) where the thick solution matches on to a thin accretion disk. The inner radius rinsubscript𝑟inr_{\rm in}italic_r start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT lies between the marginally bound and marginally stable orbits; its location determines the efficiency ε𝜀\varepsilonitalic_ε of the conversion of accreted matter into radiation.

Refer to caption
Figure 1: Meridional cross-sections (over one quadrant) for the supercritical thick disks described by Models A (magenta line) and B (green line). As rinsubscript𝑟inr_{\rm in}italic_r start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT decreases, routsubscript𝑟outr_{\rm out}italic_r start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT increases as does the ratio L/LEdd𝐿subscript𝐿EddL/L_{\rm Edd}italic_L / italic_L start_POSTSUBSCRIPT roman_Edd end_POSTSUBSCRIPT, while the efficiency of mass to energy conversion drops. Smaller values of rinsubscript𝑟inr_{\rm in}italic_r start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT imply steeper and deeper funnels, where pressure gradients are balanced by centrifugal forces rather than by gravity and luminosities exceed the Eddington limit. The square points mark the location on the surface inside which 90% of the disk luminosity Lradsubscript𝐿radL_{\rm rad}italic_L start_POSTSUBSCRIPT roman_rad end_POSTSUBSCRIPT is actually emitted. The inset shows the Keplerian specific angular momentum distribution for the adopted pseudo-Newtonian potential (black line), and the angular momentum distributions corresponding to our supercritical disks (Models A and B).

Radiation is emitted from the photosphere at the local Eddington rate

Frad=cκesgeff=cκes(Φ+2r3e^r),subscript𝐹rad𝑐subscript𝜅essubscript𝑔eff𝑐subscript𝜅esΦsuperscript2superscript𝑟3subscript^𝑒𝑟{\vec{F}}_{\rm rad}=-{c\over\kappa_{\rm es}}{\vec{g}}_{\rm eff}=-{c\over\kappa% _{\rm es}}\left(-{\vec{\nabla}}\Phi+{\ell^{2}\over r^{3}}{\hat{e}_{r}}\right),over→ start_ARG italic_F end_ARG start_POSTSUBSCRIPT roman_rad end_POSTSUBSCRIPT = - divide start_ARG italic_c end_ARG start_ARG italic_κ start_POSTSUBSCRIPT roman_es end_POSTSUBSCRIPT end_ARG over→ start_ARG italic_g end_ARG start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT = - divide start_ARG italic_c end_ARG start_ARG italic_κ start_POSTSUBSCRIPT roman_es end_POSTSUBSCRIPT end_ARG ( - over→ start_ARG ∇ end_ARG roman_Φ + divide start_ARG roman_ℓ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG over^ start_ARG italic_e end_ARG start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT ) , (3)

where geffsubscript𝑔eff{\vec{g}}_{\rm eff}over→ start_ARG italic_g end_ARG start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT is the effective gravity vector perpendicular to the surface of the disk, and κessubscript𝜅es\kappa_{\rm es}italic_κ start_POSTSUBSCRIPT roman_es end_POSTSUBSCRIPT is the electron scattering opacity. The total luminosity radiated by the torus, Lradsubscript𝐿radL_{\rm rad}italic_L start_POSTSUBSCRIPT roman_rad end_POSTSUBSCRIPT, is calculated by integrating the emitted flux over the two disk surfaces,

Lradsubscript𝐿rad\displaystyle L_{\rm rad}italic_L start_POSTSUBSCRIPT roman_rad end_POSTSUBSCRIPT =202πrinroutFrad𝑑Σ,absent2superscriptsubscript02𝜋superscriptsubscriptsubscript𝑟insubscript𝑟outsubscript𝐹raddifferential-dΣ\displaystyle=2\int_{0}^{2\pi}\int_{r_{\rm in}}^{r_{\rm out}}F_{\rm rad}\,d\Sigma,= 2 ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 italic_π end_POSTSUPERSCRIPT ∫ start_POSTSUBSCRIPT italic_r start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_F start_POSTSUBSCRIPT roman_rad end_POSTSUBSCRIPT italic_d roman_Σ , (4)

where dΣ=[1+(dh/dr)2]1/2rdrdφ𝑑Σsuperscriptdelimited-[]1superscript𝑑𝑑𝑟212𝑟𝑑𝑟𝑑𝜑d\Sigma=[1+(dh/dr)^{2}]^{1/2}\,rdrd\varphiitalic_d roman_Σ = [ 1 + ( italic_d italic_h / italic_d italic_r ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT italic_r italic_d italic_r italic_d italic_φ is the disk area element. In the case of large tori, the total luminosity generated by viscosity is simply given by

Lgen=M˙ein=M˙c2rS(rin2rS)4(rinrS)2εM˙c2,subscript𝐿gen˙𝑀subscript𝑒in˙𝑀superscript𝑐2subscript𝑟𝑆subscript𝑟in2subscript𝑟𝑆4superscriptsubscript𝑟insubscript𝑟𝑆2𝜀˙𝑀superscript𝑐2L_{\rm gen}=\dot{M}e_{\rm in}=\dot{M}c^{2}\,{r_{S}(r_{\rm in}-2r_{S})\over 4(r% _{\rm in}-r_{S})^{2}}\equiv\varepsilon\dot{M}c^{2},italic_L start_POSTSUBSCRIPT roman_gen end_POSTSUBSCRIPT = over˙ start_ARG italic_M end_ARG italic_e start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT = over˙ start_ARG italic_M end_ARG italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG italic_r start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT - 2 italic_r start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ) end_ARG start_ARG 4 ( italic_r start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT - italic_r start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ≡ italic_ε over˙ start_ARG italic_M end_ARG italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (5)

where einsubscript𝑒ine_{\rm in}italic_e start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT 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. Lgen=Lradsubscript𝐿gensubscript𝐿radL_{\rm gen}=L_{\rm rad}italic_L start_POSTSUBSCRIPT roman_gen end_POSTSUBSCRIPT = italic_L start_POSTSUBSCRIPT roman_rad end_POSTSUBSCRIPT. By equating Equations (4) and (5) we then finally calculate M˙˙𝑀\dot{M}over˙ start_ARG italic_M end_ARG.

Table 1: Properties of Thick Disks.
\topruleParameter Values for Models A and B
rin/rSsubscript𝑟insubscript𝑟𝑆r_{\rm in}/r_{S}italic_r start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT / italic_r start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT 2.2, 2.3
𝒞𝒞{\cal C}caligraphic_C 0.04792, 0.063301
rout/rSsubscript𝑟outsubscript𝑟𝑆r_{\rm out}/r_{S}italic_r start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT / italic_r start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT 320, 163
ε𝜀\varepsilonitalic_ε 0.035, 0.044
Lrad/LEddsubscript𝐿radsubscript𝐿EddL_{\rm rad}/L_{\rm Edd}italic_L start_POSTSUBSCRIPT roman_rad end_POSTSUBSCRIPT / italic_L start_POSTSUBSCRIPT roman_Edd end_POSTSUBSCRIPT 4.3, 3.1
M˙/M˙Edd˙𝑀subscript˙𝑀Edd\dot{M}/\dot{M}_{\rm Edd}over˙ start_ARG italic_M end_ARG / over˙ start_ARG italic_M end_ARG start_POSTSUBSCRIPT roman_Edd end_POSTSUBSCRIPT 12.5, 7.1
ΘΘ\Thetaroman_Θ 35, 44
Refer to caption
Figure 2: The funnel regions of supercritical disks around black holes. The two configurations (Models A and B) have funnel opening half-angles of Θ=arccot(h/r)max=35Θarccotsubscript𝑟maxsuperscript35\Theta={\rm arccot}\,(h/r)_{\rm max}=35^{\circ}roman_Θ = roman_arccot ( italic_h / italic_r ) start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT = 35 start_POSTSUPERSCRIPT ∘ end_POSTSUPERSCRIPT and 44superscript4444^{\circ}44 start_POSTSUPERSCRIPT ∘ end_POSTSUPERSCRIPT, respectively. The dot on each funnel indicates the point where h/r𝑟h/ritalic_h / italic_r reaches a maximum, and the percentage next to it denotes the fraction of total disk luminosity emitted interior to that point.

The properties of two representative super-Eddington disk models are summarized in Table 1. Values for the disk inner edge rinsubscript𝑟inr_{\rm in}italic_r start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT and the constant 𝒞𝒞{\cal C}caligraphic_C in Equation (1) were chosen following Paczyńsky & Wiita (1980). As rinsubscript𝑟inr_{\rm in}italic_r start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT 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 r2similar-toabsentsuperscript𝑟2\sim r^{-2}∼ italic_r start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT in the inner funnel. The total isotropic luminosity radiated in Models A and B is Lrad/LEdd=4.3subscript𝐿radsubscript𝐿Edd4.3L_{\rm rad}/L_{\rm Edd}=4.3italic_L start_POSTSUBSCRIPT roman_rad end_POSTSUBSCRIPT / italic_L start_POSTSUBSCRIPT roman_Edd end_POSTSUBSCRIPT = 4.3 and 3.1, corresponding to M˙/M˙Edd=12.5˙𝑀subscript˙𝑀Edd12.5\dot{M}/\dot{M}_{\rm Edd}=12.5over˙ start_ARG italic_M end_ARG / over˙ start_ARG italic_M end_ARG start_POSTSUBSCRIPT roman_Edd end_POSTSUBSCRIPT = 12.5 and 7.1, respectively. Here, we have defined the critical accretion rate, M˙Edd10LEdd/c2subscript˙𝑀Edd10subscript𝐿Eddsuperscript𝑐2\dot{M}_{\rm Edd}\equiv 10\,L_{\rm Edd}/c^{2}over˙ start_ARG italic_M end_ARG start_POSTSUBSCRIPT roman_Edd end_POSTSUBSCRIPT ≡ 10 italic_L start_POSTSUBSCRIPT roman_Edd end_POSTSUBSCRIPT / italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, 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 ΘΘ\Thetaroman_Θ:

Θ=arccot(h/r)max.Θarccotsubscript𝑟max\Theta={\rm arccot}\,(h/r)_{\rm max}.roman_Θ = roman_arccot ( italic_h / italic_r ) start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT . (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 <10<absent10\mathrel{\hbox to0.0pt{\lower 3.0pt\hbox{$\mathchar 536\relax$}\hss}\raise 2.0% pt\hbox{$\mathchar 316\relax$}}10∼< 10 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 Γ=1.72.0Γ1.72.0\Gamma=1.7–2.0roman_Γ = 1.7 – 2.0 and a high-energy cutoff. These depend on the physical properties of the scattering corona, i.e., the optical depth τ𝜏\tauitalic_τ and the electron temperature kTe𝑘subscript𝑇𝑒kT_{e}italic_k italic_T start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT (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 fcsubscript𝑓𝑐f_{c}italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT 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:

FradFinn=cκesgeff,subscript𝐹radsubscriptsuperscript𝐹𝑛in𝑐subscript𝜅essubscript𝑔eff{\vec{F}}_{\rm rad}-{\vec{F}}^{n}_{\rm in}=-{c\over\kappa_{\rm es}}{\vec{g}}_{% \rm eff},over→ start_ARG italic_F end_ARG start_POSTSUBSCRIPT roman_rad end_POSTSUBSCRIPT - over→ start_ARG italic_F end_ARG start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT = - divide start_ARG italic_c end_ARG start_ARG italic_κ start_POSTSUBSCRIPT roman_es end_POSTSUBSCRIPT end_ARG over→ start_ARG italic_g end_ARG start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT , (7)

where now Fradsubscript𝐹rad\vec{F}_{\rm rad}over→ start_ARG italic_F end_ARG start_POSTSUBSCRIPT roman_rad end_POSTSUBSCRIPT is the sum of the flux locally generated (‘self-flux’) and the normal component of the incoming flux Finnsubscriptsuperscript𝐹𝑛in{\vec{F}}^{n}_{\rm in}over→ start_ARG italic_F end_ARG start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT.

Consider first, for illustration, a geometrically thin disk scenario with Fin=0subscript𝐹in0{F}_{\rm in}=0italic_F start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT = 0, and let us write the total Compton luminosity per unit area in the corona as the sum of an upward-, Fc+superscriptsubscript𝐹𝑐F_{c}^{+}italic_F start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT, and downward-directed, Fcsuperscriptsubscript𝐹𝑐F_{c}^{-}italic_F start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT, component. Ignoring the anisotropy of inverse Compton scattering, we can express the power radiated in all directions by the hot phase as 2Fc+=ApFs2superscriptsubscript𝐹𝑐𝐴𝑝subscript𝐹𝑠2F_{c}^{+}=ApF_{s}2 italic_F start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT = italic_A italic_p italic_F start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT, where A𝐴Aitalic_A is the mean energy gain per scattering, Fssubscript𝐹𝑠F_{s}italic_F start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT is the total soft photon flux emerging from the cool disk, and p=1exp(2τ)𝑝12𝜏p=1-\exp(-2\tau)italic_p = 1 - roman_exp ( - 2 italic_τ ) is the mean probability of scattering in the corona. The outgoing photon flux is Frad=Fs+Fc+subscript𝐹radsubscript𝐹𝑠superscriptsubscript𝐹𝑐F_{\rm rad}=F_{s}+F_{c}^{+}italic_F start_POSTSUBSCRIPT roman_rad end_POSTSUBSCRIPT = italic_F start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT + italic_F start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT. 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

Fs=subscript𝐹𝑠absent\displaystyle F_{s}=italic_F start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = (1fc)cκesgeff+12ApFs,1subscript𝑓𝑐𝑐subscript𝜅essubscript𝑔eff12𝐴𝑝subscript𝐹𝑠\displaystyle(1-f_{c})\,{c\over\kappa_{\rm es}}g_{\rm eff}+{1\over 2}ApF_{s},( 1 - italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ) divide start_ARG italic_c end_ARG start_ARG italic_κ start_POSTSUBSCRIPT roman_es end_POSTSUBSCRIPT end_ARG italic_g start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_A italic_p italic_F start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT , (8)
Fc+superscriptsubscript𝐹𝑐\displaystyle F_{c}^{+}italic_F start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT =fc2cκesgeff.absentsubscript𝑓𝑐2𝑐subscript𝜅essubscript𝑔eff\displaystyle={f_{c}\over 2}\,{c\over\kappa_{\rm es}}g_{\rm eff}.= divide start_ARG italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG divide start_ARG italic_c end_ARG start_ARG italic_κ start_POSTSUBSCRIPT roman_es end_POSTSUBSCRIPT end_ARG italic_g start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT .

Solving for A𝐴Aitalic_A and using the previous definitions, we derive

A=2fc(2fc)p.𝐴2subscript𝑓𝑐2subscript𝑓𝑐𝑝A={2f_{c}\over(2-f_{c})p}.italic_A = divide start_ARG 2 italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG start_ARG ( 2 - italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ) italic_p end_ARG . (9)

In the regime of unsaturated Comptonization, the mean energy gain per scattering is given by A=4θ+16θ2𝐴4𝜃16superscript𝜃2A=4\theta+16\theta^{2}italic_A = 4 italic_θ + 16 italic_θ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, where θ=kTe/mec2𝜃𝑘subscript𝑇𝑒subscript𝑚𝑒superscript𝑐2\theta=kT_{e}/m_{e}c^{2}italic_θ = italic_k italic_T start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT / italic_m start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT 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 ΓΓ\Gammaroman_Γ of the power-law Comptonized component FννΓ+1proportional-tosubscript𝐹𝜈superscript𝜈Γ1F_{\nu}\propto\nu^{-\Gamma+1}italic_F start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT ∝ italic_ν start_POSTSUPERSCRIPT - roman_Γ + 1 end_POSTSUPERSCRIPT (Haardt & Maraschi, 1993). The Comptonization radiative transfer models of Sunyaev & Titarchuk (1980) give

Γ=lnτ+2/(3+θ)ln(12θ2+25θ)+1.Γ𝜏23𝜃12superscript𝜃225𝜃1\Gamma={-\ln\tau+2/(3+\theta)\over\ln(12\theta^{2}+25\theta)}+1.roman_Γ = divide start_ARG - roman_ln italic_τ + 2 / ( 3 + italic_θ ) end_ARG start_ARG roman_ln ( 12 italic_θ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 25 italic_θ ) end_ARG + 1 . (10)

According to the above expressions, spectral indeces and electron temperatures are in the range Γ=1.702.0Γ1.702.0\Gamma=1.70-2.0roman_Γ = 1.70 - 2.0 and kTe>70>𝑘subscript𝑇𝑒70kT_{e}\mathrel{\hbox to0.0pt{\lower 3.0pt\hbox{$\mathchar 536\relax$}\hss}% \raise 2.0pt\hbox{$\mathchar 318\relax$}}70\,italic_k italic_T start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ∼> 70 keV for 0.3<τ<0.8<0.3𝜏<0.80.3\mathrel{\hbox to0.0pt{\lower 3.0pt\hbox{$\mathchar 536\relax$}\hss}\raise 2% .0pt\hbox{$\mathchar 316\relax$}}\tau\mathrel{\hbox to0.0pt{\lower 3.0pt\hbox{% $\mathchar 536\relax$}\hss}\raise 2.0pt\hbox{$\mathchar 316\relax$}}0.80.3 ∼< italic_τ ∼< 0.8 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 fc=0.5subscript𝑓𝑐0.5f_{c}=0.5italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = 0.5. Such hard X-ray photon indeces are seen above 2 keV in low-redshift radio-quiet AGNs (Γ=1.89±0.11Γplus-or-minus1.890.11\Gamma=1.89\pm 0.11roman_Γ = 1.89 ± 0.11; Piconcelli et al., 2005). The highest luminosity AGNs at z>6>𝑧6z\mathrel{\hbox to0.0pt{\lower 3.0pt\hbox{$\mathchar 536\relax$}\hss}\raise 2.% 0pt\hbox{$\mathchar 318\relax$}}6italic_z ∼> 6 have X-ray spectra that are steeper than these values (Γ=2.4±0.1Γplus-or-minus2.40.1\Gamma=2.4\pm 0.1roman_Γ = 2.4 ± 0.1; Zappacosta et al., 2023), and so do near-Eddington narrow-line Seyfert 1 (NLS1) galaxies (Γ=2.68±0.51Γplus-or-minus2.680.51\Gamma=2.68\pm 0.51roman_Γ = 2.68 ± 0.51; 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

Fs=subscript𝐹𝑠absent\displaystyle F_{s}=italic_F start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = (1fc)cκesgeff+12ApFs+mcκesgeff(1p),1subscript𝑓𝑐𝑐subscript𝜅essubscript𝑔eff12𝐴𝑝subscript𝐹𝑠𝑚𝑐subscript𝜅essubscript𝑔eff1𝑝\displaystyle(1-f_{c})\,{c\over\kappa_{\rm es}}g_{\rm eff}+{1\over 2}ApF_{s}+m% {c\over\kappa_{\rm es}}g_{\rm eff}(1-p),( 1 - italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ) divide start_ARG italic_c end_ARG start_ARG italic_κ start_POSTSUBSCRIPT roman_es end_POSTSUBSCRIPT end_ARG italic_g start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_A italic_p italic_F start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT + italic_m divide start_ARG italic_c end_ARG start_ARG italic_κ start_POSTSUBSCRIPT roman_es end_POSTSUBSCRIPT end_ARG italic_g start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( 1 - italic_p ) , (11)
Fc+superscriptsubscript𝐹𝑐\displaystyle F_{c}^{+}italic_F start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT =fc2cκesgeff,absentsubscript𝑓𝑐2𝑐subscript𝜅essubscript𝑔eff\displaystyle={f_{c}\over 2}\,{c\over\kappa_{\rm es}}g_{\rm eff},= divide start_ARG italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG divide start_ARG italic_c end_ARG start_ARG italic_κ start_POSTSUBSCRIPT roman_es end_POSTSUBSCRIPT end_ARG italic_g start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ,

where the third term is the incoming external flux expressed as a multiple m𝑚mitalic_m of the self-flux cgeff/κes𝑐subscript𝑔effsubscript𝜅es{cg_{\rm eff}/\kappa_{\rm es}}italic_c italic_g start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT / italic_κ start_POSTSUBSCRIPT roman_es end_POSTSUBSCRIPT times a factor 1p1𝑝1-p1 - italic_p to account for the fraction that scatters off the corona. It has long been recognized that in narrow accretion funnels (small half-angle ΘΘ\Thetaroman_Θ) the self-flux FradFinsubscript𝐹radsubscript𝐹in{F}_{\rm rad}-{F}_{\rm in}italic_F start_POSTSUBSCRIPT roman_rad end_POSTSUBSCRIPT - italic_F start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT is a much smaller quantity than either Fradsubscript𝐹rad{F}_{\rm rad}italic_F start_POSTSUBSCRIPT roman_rad end_POSTSUBSCRIPT or Finsubscript𝐹in{F}_{\rm in}italic_F start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT because of the strong “reflection effect” (Sikora, 1981; Narayan et al., 1983; Madau, 1988). Solving for the amplification factor A𝐴Aitalic_A as before we now have

A=2fc[(2fc)+2m(1p)]p.𝐴2subscript𝑓𝑐delimited-[]2subscript𝑓𝑐2𝑚1𝑝𝑝A={2f_{c}\over[(2-f_{c})+2m(1-p)]p}.italic_A = divide start_ARG 2 italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG start_ARG [ ( 2 - italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ) + 2 italic_m ( 1 - italic_p ) ] italic_p end_ARG . (12)

Figure 3 demonstrates how, with increasing values of the paramater m𝑚mitalic_m, AGN X-ray spectra become extremely soft as the corona cools down. Our scenario is schematically illustrated in Figure 4.

Refer to caption
Figure 3: Equilibrium electron temperature θ=kTe/mec2𝜃𝑘subscript𝑇𝑒subscript𝑚𝑒superscript𝑐2\theta=kT_{e}/m_{e}c^{2}italic_θ = italic_k italic_T start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT / italic_m start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT) of the hot coronal phase (lower curves) and the resulting photon index ΓΓ\Gammaroman_Γ of the X-ray Comptonized component (upper curves). The corona has an assumed electron scattering opacity of τ=0.5𝜏0.5\tau=0.5italic_τ = 0.5, comparable to the value inferred in Seyfert galaxies. Two values, fc=1/2subscript𝑓𝑐12f_{c}=1/2italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = 1 / 2 and fc=1/3subscript𝑓𝑐13f_{c}=1/3italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = 1 / 3, have been assumed for the fraction of the locally generated radiation power that emerges from the corona. The parameter m𝑚mitalic_m on the horizontal axis measures the strength of the incoming external radiation field in units of the local self-flux.

To estimate m𝑚mitalic_m, we have subdivided the surface of the funnel into a finite number of r𝑟ritalic_r-rings with quasi-uniform emission properties, and approximated the force balance Equation (7) as

Frad,i1πjBijFrad,j=cκesgeff,i,subscript𝐹radi1𝜋subscript𝑗subscript𝐵𝑖𝑗subscript𝐹radj𝑐subscript𝜅essubscript𝑔effiF_{\rm rad,i}-{1\over\pi}\sum_{j}B_{ij}\,F_{\rm rad,j}={c\over\kappa_{\rm es}}% {g}_{\rm eff,i},italic_F start_POSTSUBSCRIPT roman_rad , roman_i end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG italic_π end_ARG ∑ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT italic_F start_POSTSUBSCRIPT roman_rad , roman_j end_POSTSUBSCRIPT = divide start_ARG italic_c end_ARG start_ARG italic_κ start_POSTSUBSCRIPT roman_es end_POSTSUBSCRIPT end_ARG italic_g start_POSTSUBSCRIPT roman_eff , roman_i end_POSTSUBSCRIPT , (13)

where

Bij=n^iD|D|[(n^jD)|D|3dΣj]subscript𝐵𝑖𝑗subscript^𝑛𝑖𝐷𝐷delimited-[]subscript^𝑛𝑗𝐷superscript𝐷3𝑑subscriptΣ𝑗B_{ij}={{\hat{n}_{i}}\cdot\vec{D}\over|{\vec{D}}|}\left[{(\hat{n}_{j}\cdot\vec% {D})\over|\vec{D}|^{3}}d\Sigma_{j}\right]italic_B start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT = divide start_ARG over^ start_ARG italic_n end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_D end_ARG end_ARG start_ARG | over→ start_ARG italic_D end_ARG | end_ARG [ divide start_ARG ( over^ start_ARG italic_n end_ARG start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_D end_ARG ) end_ARG start_ARG | over→ start_ARG italic_D end_ARG | start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG italic_d roman_Σ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ] (14)

is the fraction of radiation emitted from the j𝑗jitalic_jth ring that reaches the i𝑖iitalic_ith ring along the direction D𝐷\vec{D}over→ start_ARG italic_D end_ARG in the solid angle given by the term in square brackets. Here, n^jsubscript^𝑛𝑗{\hat{n}}_{j}over^ start_ARG italic_n end_ARG start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT is the outward unit normal to the surface at the element dΣj=rj[1+(dh/dr)j2]1/2dφΔr𝑑subscriptΣ𝑗subscript𝑟𝑗superscriptdelimited-[]1superscriptsubscript𝑑𝑑𝑟𝑗212𝑑𝜑Δ𝑟d\Sigma_{j}=r_{j}[1+(dh/dr)_{j}^{2}]^{1/2}\,d\varphi\Delta ritalic_d roman_Σ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT = italic_r start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT [ 1 + ( italic_d italic_h / italic_d italic_r ) start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT italic_d italic_φ roman_Δ italic_r, and ΔrΔ𝑟\Delta rroman_Δ italic_r is the mesh size. Solving this set of linear equations for Models A and B gives m3similar-to-or-equals𝑚3m\simeq 3italic_m ≃ 3 and m2similar-to-or-equals𝑚2m\simeq 2italic_m ≃ 2 at ri=10rSsubscript𝑟𝑖10subscript𝑟𝑆r_{i}=10\,r_{S}italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = 10 italic_r start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT, respectively. Models with higher supercritical rates (smaller ΘΘ\Thetaroman_Θ) yield values of m>5>𝑚5m\mathrel{\hbox to0.0pt{\lower 3.0pt\hbox{$\mathchar 536\relax$}\hss}\raise 2.% 0pt\hbox{$\mathchar 318\relax$}}5italic_m ∼> 5 near the bottom of the funnel.

Refer to caption
Figure 4: Sketch illustrating the super-Eddington corona-funnel model discussed in the main text. Embedded in a funnel-like reflection geometry, the inner hot corona will cool down by Comptonizing a largely isotropic soft radiation field. The resulting X-ray spectra are predicted to be extremely soft, independently of the detailed geometry of the corona.

Compare now a “standard” AGN corona in the open geometry of a thin accretion disk. Assuming fc=0.5subscript𝑓𝑐0.5f_{c}=0.5italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = 0.5 and τ=0.5𝜏0.5\tau=0.5italic_τ = 0.5, the two-phase thin disk-corona model produces a cutoff power-law X-ray spectrum with Γ=1.9Γ1.9\Gamma=1.9roman_Γ = 1.9 and kTe=80𝑘subscript𝑇𝑒80kT_{e}=80\,italic_k italic_T start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT = 80keV. The derived rest-frame 2102102-102 - 10 keV bolometric correction, k210Lbol/L210=16subscript𝑘210subscript𝐿bolsubscript𝐿21016k_{2-10}\equiv L_{\rm bol}/L_{2-10}=16italic_k start_POSTSUBSCRIPT 2 - 10 end_POSTSUBSCRIPT ≡ italic_L start_POSTSUBSCRIPT roman_bol end_POSTSUBSCRIPT / italic_L start_POSTSUBSCRIPT 2 - 10 end_POSTSUBSCRIPT = 16 (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 m>3>𝑚3m\mathrel{\hbox to0.0pt{\lower 3.0pt\hbox{$\mathchar 536\relax$}\hss}\raise 2.% 0pt\hbox{$\mathchar 318\relax$}}3italic_m ∼> 3 and the same values of fcsubscript𝑓𝑐f_{c}italic_f start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT and τ𝜏\tauitalic_τ, we predict instead Γ>2.8>Γ2.8\Gamma\mathrel{\hbox to0.0pt{\lower 3.0pt\hbox{$\mathchar 536\relax$}\hss}% \raise 2.0pt\hbox{$\mathchar 318\relax$}}2.8roman_Γ ∼> 2.8, kTe<40<𝑘subscript𝑇𝑒40kT_{e}\mathrel{\hbox to0.0pt{\lower 3.0pt\hbox{$\mathchar 536\relax$}\hss}% \raise 2.0pt\hbox{$\mathchar 316\relax$}}40\,italic_k italic_T start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ∼< 40keV, and a bolometric correction of k210>68>subscript𝑘21068k_{2-10}\mathrel{\hbox to0.0pt{\lower 3.0pt\hbox{$\mathchar 536\relax$}\hss}% \raise 2.0pt\hbox{$\mathchar 318\relax$}}68italic_k start_POSTSUBSCRIPT 2 - 10 end_POSTSUBSCRIPT ∼> 68. We observe here that comparable values of k210=5080subscript𝑘2105080k_{2-10}=50-80italic_k start_POSTSUBSCRIPT 2 - 10 end_POSTSUBSCRIPT = 50 - 80 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 m>5>𝑚5m\mathrel{\hbox to0.0pt{\lower 3.0pt\hbox{$\mathchar 536\relax$}\hss}\raise 2.% 0pt\hbox{$\mathchar 318\relax$}}5italic_m ∼> 5, our model yields Γ>3.95>Γ3.95\Gamma\mathrel{\hbox to0.0pt{\lower 3.0pt\hbox{$\mathchar 536\relax$}\hss}% \raise 2.0pt\hbox{$\mathchar 318\relax$}}3.95roman_Γ ∼> 3.95, kTe<30<𝑘subscript𝑇𝑒30kT_{e}\mathrel{\hbox to0.0pt{\lower 3.0pt\hbox{$\mathchar 536\relax$}\hss}% \raise 2.0pt\hbox{$\mathchar 316\relax$}}30\,italic_k italic_T start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ∼< 30keV, and k210>1600>subscript𝑘2101600k_{2-10}\mathrel{\hbox to0.0pt{\lower 3.0pt\hbox{$\mathchar 536\relax$}\hss}% \raise 2.0pt\hbox{$\mathchar 318\relax$}}1600italic_k start_POSTSUBSCRIPT 2 - 10 end_POSTSUBSCRIPT ∼> 1600. For comparison, the X-ray stacking of JWST-selected, faint type 1 AGNs at z5similar-to-or-equals𝑧5z\simeq 5italic_z ≃ 5 gives k210>500subscript𝑘210500k_{2-10}>500italic_k start_POSTSUBSCRIPT 2 - 10 end_POSTSUBSCRIPT > 500 (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 M˙/M˙Edd˙𝑀subscript˙𝑀Edd\dot{M}/\dot{M}_{\rm Edd}over˙ start_ARG italic_M end_ARG / over˙ start_ARG italic_M end_ARG start_POSTSUBSCRIPT roman_Edd end_POSTSUBSCRIPT 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 10%similar-toabsentpercent10\sim 10\%∼ 10 % 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 z=6.67𝑧6.67z=6.67italic_z = 6.67 (Juodžbalis et al., 2024), which seems to be accreting well below the Eddington limit and is undetected in X-rays with k210>330subscript𝑘210330k_{2-10}>330italic_k start_POSTSUBSCRIPT 2 - 10 end_POSTSUBSCRIPT > 330 (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