DRAFTVERSION APRIL 13, 2021 Typeset using LATEX default style in AASTeX62
Accretion-modified Stars in Accretion Disks of Active Galactic Nuclei: Slowly Transient Appearance
JIAN-MIN WANG,1, 2, 3 JUN-RONG LIU,1, 2 LUIS C.HO,4, 5 AND PU DU1
1Key Laboratory for Particle Astrophysics, Institute of High Energy Physics, Chinese Academy of Sciences, 19B Yuquan Road, Beijing 100049, China 2School of Astronomy and Space Sciences, University of Chinese Academy of Sciences, 19A Yuquan road, Beijing 100049, China 3National Astronomical Observatory of China, 20A Datun Road, Beijing 100020, China 4Kavli Institute for Astronomy and Astrophysics, Peking University, Beijing 100871, China 5Department of Astronomy, School of Physics, Peking University, Beijing 100871, China
(Received 2020 Feburary 11; Revised 2021 March 10, 2021; Accepted 2021 March 12, 2021)
ABSTRACT Compact objects are expected to exist in the accretion disks of supermassive black holes (SMBHs) in active galactic nuclei (AGNs), and in the presence of such a dense environment ( 1014 cm−3), they will form a new ∼ kind of stellar population denoted as Accretion-Modified Stars (AMSs). This hypothesis is supported by recent LIGO/Virgo detection of the mergers of very high-mass stellar binary black holes (BHs). We show that the AMSs will be trapped by the SMBH-disk within a typical AGN lifetime. In the context of SMBH-disks, the 9 2 rates of Bondi accretion onto BHs are 10 LEdd/c , where LEdd is the Eddington luminosity and c is the speed ∼ of light. Outflows developed from the hyper-Eddington accretion strongly impact the Bondi sphere and induce episodic accretion. We show that the hyper-Eddington accretion will be halted after an accretion interval of 5 ta 10 m1 s, where m1 = m•/10M is the BH mass. The kinetic energy of the outflows accumulated during ∼ ta is equivalent to 10 supernovae driving an explosion of the Bondi sphere and developing blast waves. We demonstrate that a synchrotron flare from relativistic electrons accelerated by the blast waves peaks in the soft X-ray band ( 0.1 keV), significantly contributing to the radio, optical, UV, and soft X-ray emission of typical ∼ radio-quiet quasars. External inverse Compton scattering of the electrons peaks around 40 GeV and is detectable through Fermi-LAT. The flare, decaying with t−6/5 with a few months, will appear as a slowly varying transient. The flares, occurring at a rate of a few per year in radio-quiet quasars, provide a new mechanism for explaining AGN variability.
Keywords: Active galactic nuclei (16); Galaxy accretion disks (562); Supermassive black holes (1663)
1. INTRODUCTION Compact objects (neutron stars and stellar black holes) may exist in accretion disks of active galactic nuclei (AGNs) and quasars. On the one hand, pioneering ideas of self-gravitating accretion disks in AGNs (Paczynski´ 1977; Kolykhalov & Sunyaev 1980; Shlosman & Begelman 1989) have received increasing attention (Collin & Zahn 1999, 2008; Goodman 2003; Goodman & Tan 2004; Wang et al. 2010, 2011, 2012), for they offer a possible explanation for the super-solar metallicities inferred in AGNs and quasars (Hamann & Ferland 1999; Warner et al. 2003; Nagao et al. 2006; Shin et al. 2013; Du & Wang 2014). The high metallicities can be naturally linked to star formation in self-gravitating disks, which inevitably produces compact objects from arXiv:2103.07708v2 [astro-ph.HE] 10 Apr 2021 supernova explosions. On the other hand, the high metallicity of quasars can also be explained by stars from nuclear clusters captured by SMBH-disks (Artymowicz et al. 1993), a process that will also introduce compact objects in the disks (Cheng & Wang 1999; Cantiello et al. 2020). Enveloped by the very dense gaseous medium of the accretion disk, compact objects inevitably will form a special kind of object, whose envelope is very massive, even more massive that the central compact object itself (see §2.2). We use the terminology of Accretion-Modified Stars (AMSs) to denote these special objects in the SMBH-disks1. The fate of AMSs in such an environment, however, is unknown.
1 The original concept of the Thorne-Zytkow˙ objects (TZOs: Thorne & Zytkow˙ 1975, 1977), which only pertained to neutron stars enshrouded by envelopes are in both hydrostatic and thermal equilibrium, with some accretion occurring just in the very inner region close to the neutron stars. They are expected to have specific surface abundances due to exotic nuclear burning processes happening deeper into the star and close to the compact object. The present AMSs are neither in hydrostatic nor thermal equilibrium. These compact objects are episodically accreting from the Bondi sphere. 2 WANGETAL.
Bondi sphere
outflows BH SMBH-disk �a SMBH RBon
2 Figure 1. Left: Stellar mass black holes (1M − 10 M ) form Accretion-modified Stars (AMSs) in AGN accretion disks. Right: The Bondi sphere hosts a BH undergoing hyper-Eddington accretion. Powerful outflows developed from the slim disk of the BH produce strong feedback 4 to the Bondi sphere. The hyper-accretion is halted by the mechanical feedback at the radius ra(∼ 10 gravitational radii) within interval ta, generating flares from radio to γ-rays (see text for details). In this paper, we neglect the spins of the Bondi sphere (see a brief discussion in §2.2).
The recent detection of GW190521, consisting of (85+66)M binary BHs, by the Advanced LIGO/Virgo consortium (Abbott et al. 2020) has received much attention because it is harbored by the quasar SDSS J1249+3449 monitored by the Zwicky Transient Facility (Graham et al. 2020). BH binaries with such high masses far exceeds the limit of pair instability of massive stars (e.g., Woosley 2017). As the BH masses of GW190521 are much higher than the upper limit of supernova explosions from an isolated star, this event has inspired the idea that mergers of compact objects can occur in AGN disks (Cheng & Wang 1999; Bartos et al. 2017; McKernan et al. 2019; Yang et al. 2019, 2020; Tanaga et al. 2020; Samsing et al. 2020). Such a scenario would involve AMSs in SMBH-disks. In this paper, we suggest that BHs are undergoing episodic accretion governed by powerful outflows developed from hyper- Eddington accretion. The outflows drive a Bondi explosion and result in non-thermal flares from the blast waves that are predicted to appear as slowly varying transients in the radio, optical, UV, and soft X-ray bands. AMSs offer a new mechanism for explaining AGN variability but also for origins of massive stellar black holes.
2. A SIMPLE MODEL Stars, which originate either from capture from nuclear star clusters or direct formation from self-gravitating disks, accrete gas from SMBH-disks to form more massive objects that evolve quickly into compact objects (Artymowicz et al. 1993; Cheng & Wang 1999; Cantiello et al. 2020). The massive envelope around a compact object forms an AMS through accretion of gas from the local SMBH-disk. In this paper, we only focus on AMSs of black holes. Other more complicated options can arise. −2 The dimensionless accretion rate of the central SMBH is defined by M˙ = M˙ •/M˙ Edd, where M˙ Edd = LEddc is the 46 −1 Eddington limit rate, M˙ • is the accretion rate of the SMBH, LEdd = 4πGM•mpc/σT = 1.3 10 M8 erg s , G is the 8 8 × gravitational constant, M8 = M•/10 M is the SMBH mass in units of 10 M , c is the speed of light, mp is the proton mass, and σT is the Thomson cross section. The half-thickness, density, mid-plane temperature, and radial velocity of the SMBH-disk are (e.g., Kato et al. 2008) 14 −1/10 9/10 ˙3/20 9/8 H = 4.3 10 α0.1 M8 M r4 cm, × −11 −7/10 ˙11/20 −15/8 −3 ρd = 6.9 10 (α0.1M8) M r4 g cm , × (1) 3 −1/5 ˙3/10 −3/4 Tc = 4.6 10 (α0.1M8) M r4 K, × 2 4/5 −1/5 3/10 −1/4 −1 vr = 2.6 10 α M M˙ r cm s , × 0.1 8 4 4 2 where α0.1 = α/0.1 is the viscosity parameter, r4 = R/10 Rg is the radius of the disk from the SMBH, and Rg = GM•/c is the gravitational radius. The self-gravity of the disk can be described by the Toomre parameter, defined as Q = ΩKcs/πGρdH, p 1/2 −1 p 3 where cs 3kTc/mp 15.7 T km s is the local sound speed of SMBH-disk and ΩK = GM•/R . The disk becomes ≈ ≈ 4 self-gravitating (SG) beyond a critical radius where Q = 1, which is given by
3 28/45 −52/45 −22/45 RSG/Rg = 1.2 10 α M M˙ . (2) × 0.1 8
We consider the regions beyond RSG, where massive stars are formed either through capture by the SMBH-disk or they are formed in situ, likely with a top-heavy initial mass function because of the high temperatures in the nuclear environment. ACCRETION-MODIFIED STARS IN AGNS 3
We would like to point out that the self-gravitating regions are undergoing more complicated physics beyond the solutions of Eqn. (1). For example, star formation is unavoidable (Collin & Zahn 1999, 2008; Sirko & Goodman 2003; Thompson et al. 2005), gravito turbulence dominates (Rafikov 2007, 2015), and the region has a clumpy distribution. The actual density and temperature could be different from that used here so that the fate of AMSs could be more complicated than currently assumed. The main characteristics of their evolution, however, can be predicted qualitatively. Future studies will be carried out about compact objects colliding with gaseous clumps in the self-gravitating regions.
2.1. AMSs Considering that black holes formed from massive stars are kicked off with random velocities of v0 102 103 km s−1 with • ≈ − respect to their remnant (the black holes are still tightly bound by the SMBHs), accretion onto the black holes in SMBH-disks can be described by the Hoyle-Lyttleton-Bondi (HLB) formulation. Neglecting the weak dependence on adiabatic index and self-gravity of the accreted gas (Wandel 1984), we take the simplest form
2 2 4πG m•ρd M˙ Bon = , (3) 2 2 3/2 (cs + v•) where m• is the black hole mass from supernova explosion and v• is the kick-off velocity relative to the motion of the SMBH- disk. The kick-off velocity of the black holes is uncertain, but it could be up to v0 103 km s−1 if comparable with the • ∼ case of neutron stars (e.g., Nakamura et al. 2019). We consider two extreme scenarios, in which the black hole is co-rotating 0 −1/2 −1 and counterrotating with the SMBH-disk. If v• VK = 3000 r4 km s , with VK the Keplerian velocity, we expect two ∼ 0 populations of accreting black holes in the disk: co-rotating black holes with v• VK v• cs and counterrotating ones with 0 ≈ | − | v• VK + v 2VK cs (see §2.2). It should be stressed that the current choice of kick-off velocities are just for the two ≈ • ≈ possible regimes, and is not necessarily representative of the bulk population. The corresponding accretion rates are 8 −2/5 ˙1/10 −3/4 M˙ Bon 8.9 10 m1(α0.1M8) M r4 (for corotation), m˙ Bon = × (4) −3/8 M˙ Edd ≈ −7/10 ˙11/20 5.1 m1(α0.1M8) M r4 (for counter rotation),