Arxiv:2001.11032V2 [Astro-Ph.GA] 25 Apr 2020 Plain the Nature of Low-Luminosity Accreting Smbhs (Ichimaru 1977; Narayan & Yi 1994; Stone Et Al

Draft version April 28, 2020 Typeset using LATEX twocolumn style in AASTeX63

Universal transition diagram from dormant to actively accreting supermassive black holes

Kohei Inayoshi,1 Kohei Ichikawa,2, 3 and Luis C. Ho1, 4

1Kavli Institute for Astronomy and Astrophysics, Peking University, Beijing 100871, China 2Astronomical Institute, Graduate School of Science Tohoku University, 6-3 Aramaki, Aoba-ku, Sendai 980-8578, Japan 3Frontier Research Institute for Interdisciplinary Sciences, Tohoku University, Sendai 980-8578, Japan 4Department of Astronomy, School of Physics, Peking University, Beijing 100871, China

(Received April 28, 2020)

ABSTRACT The vast majority of supermassive black holes (SMBHs) in the local universe exhibit levels of activity much lower than those expected from gas supplying rates onto the galactic nuclei, and only a small fraction of silent SMBHs can turn into active galactic nuclei. Revisiting observational data of very nearby SMBHs whose gravitational spheres of influence are spatially reached by the Chandra X-ray satellite, we find that the level of BH activity drastically increases from the quiescent phase when the inflow rate outside of the BH influence radius is higher than 0.1% of the Eddington accretion rate. We also show that the relation between the nuclear luminosity and gas accretion rate from the BH influence radius measured from X-ray observations is well described by the universal state transition of accreting SMBHs, as predicted by recent hydrodynamical simulations with radiative cooling and BH feedback. After the state transition, young massive stars should form naturally in the nucleus, as observed in the case of the nearest SMBH, Sagittarius A∗, which is currently quiescent but was recently active.

Keywords: galaxies: nuclei — galaxies: Seyfert — quasars: supermassive black holes.

1. INTRODUCTION vatories around the globe and succeeded in imaging the Supermassive black holes (SMBHs) are ubiquitously accretion flow onto the M87 SMBH and presumably its harbored at the centers of massive nearby galaxies, and black hole shadow (Event Horizon Telescope Collabora- are believed to coevolve with their host galaxies through tion et al. 2019a,b). By comparing the observed image BH feeding and energetic feedback over cosmic time (Ko- to those obtained by ray-traced general relativistic mag- rmendy & Ho 2013, hereafter KH13). Most SMBHs in netohydrodynamic simulations of BH accretion, our un- the local universe are nearly quiescent, low-luminosity derstanding of gas dynamics in the strong-field regime active galactic nuclei (AGNs) with bolometric luminosi- of general relativity will be improved significantly (e.g., −3 Dexter et al. 2012; Mo´scibrodzka et al. 2016; Ryan et al. ties Lbol/LEdd 10 , where LEdd is the Eddington luminosity (Ho 2008, 2009). The radiatively inefficient 2018). accretion flow (RIAF) model has been proposed to ex- Observations of nearby SMBHs also provide a great opportunity to explore BH accretion dynamics at the arXiv:2001.11032v2 [astro-ph.GA] 25 Apr 2020 plain the nature of low-luminosity accreting SMBHs (Ichimaru 1977; Narayan & Yi 1994; Stone et al. 1999; BH gravitational influence radius (hereafter, Bondi ra- Narayan et al. 2000; Quataert & Gruzinov 2000; Yuan dius), which is on the order of ∼ 1 − 100 parsecs ( & Narayan 2014). Two notable examples are Sagittar- the event horizon scales), from which mass inflows oc- ˙ ius A∗ (Sgr A∗) and the SMBH at the center of the cur at the so-called Bondi accretion rate MB. This giant elliptical galaxy NGC 4486 (M87). Recently, the Bondi accretion rate is supposed to be an upper limit of Event Horizon Telescope project pulled together obser- the BH feeding rate and lead to radiative output with −3 Lbol/LEdd ≈ 10 . However, this level of energy release is far more luminous than actually observed, assum- Corresponding author: Kohei Inayoshi, Kohei Ichikawa ing a canonical value of 10% for the radiative efficiency [email protected], [email protected] (Shakura & Sunyaev 1973) or an even lower value (∼ a 2 Inayoshi et al. few %) for RIAFs (Xie & Yuan 2012; Ryan et al. 2017). et al. 2009), NGC 3115 (Wong et al. 2014), and NGC The deficit of radiative luminosity suggests that only a 1407 (Humphrey et al. 2006). Moreover, Russell et al. fraction of inflowing matter from the Bondi scale ac- (2013, hereafter R13) compiled Chandra observations to tually feeds the BH and produces radiative/mechanical discuss the Bondi accretion rates for the SMBHs in 13 output that we observe in the nuclei of low-luminosity nearby galaxies, providing us with eight new sources1. SMBHs. Therefore, drawing the relation between the Finally, we have added four galaxies (NGC 315, NGC mass inflow rate and radiative output is essential to un- 2681, NGC 4278, and NGC 5005) with clear hot diffuse derstand the mechanisms governing the fate of the fuel gas emission (Younes et al. 2011). In total, we have reservoir, and to obtain a consistent picture to link quan- added 16 more sources, resulting in the total number of titatively dormant SMBHs and more luminous AGNs. 31 sources. In this paper, we revisit observational data of very In what follows, we describe the manner of observa- nearby SMBHs for which the spatial resolution of the tional data assembly for the 31 low-luminosity SMBHs. Chandra X-ray satellite reaches the Bondi scales, and The estimated physical quantities for those objects are report that the observational data clearly show the ex- summarized in Table1. istence of a universal state transition for BH activity. The level of radiative output before the state transition 2.1. Distance (D) and critical mass inflow rate for the transition can be It is crucial to obtain a homogeneous set of reliable described well by recent numerical simulations taking distances for our sample. Since most of our sample are into account radiative processes and BH feedback. very nearby sources, redshift-independent distances are used whenever possible. We follow the priorities for the 2. SAMPLE ASSEMBLY choices of distance from KH13 (see their Tables 2 and 3). For 23 out of the 31 sources, we use distances based on High-spatial resolution X-ray observations using the surface brightness fluctuation measurements for individ- Chandra satellite give unique information on the prop- ual galaxies in the Virgo and Fornax clusters (Blakeslee erties of inflowing gas around the Bondi radius, where et al. 2009, 2010). For NGC 1291, we adopt the distance the gas is turned into hot and diffuse plasma with based on the tip of the red giant branch (McQuinn et al. temperature T ≈ 0.2 − 1 keV and electron density 2017). For Sgr A∗, we adopt the distance measured from n ≈ 0.1 − 10 cm−3. Even with the excellent spatial e resolved stellar dynamics (Genzel et al. 2010). Distances resolution of ∼ 0.500, the Bondi scale, at most on the for the remaining six sources are taken from the mean order of ∼ 100 pc, is resolvable only for a small number values of several distance determinations listed in NED, of sources at distances up to 50 Mpc. mainly based on Cepheid variables, surface brightness In this paper, we adopt the compilation by Pellegrini fluctuations, tip of the red giant branch, and RR Lyrae (2005), which contains a total of 50 X-ray-faint galax- 42 −1 stars. ies with X-ray luminosities L2−10 . 10 erg s in the 2–10 keV band, for which Chandra data are publicly 2.2. Black Hole Mass (M•) available as of 2005. For 15 of these sources, the phys- We collect the BH masses of 16 sources based on spa- ical parameters of diffuse hot gas surrounding the nu- tially resolved stellar or ionized gas kinematic observa- clear regions, such as thermal temperature and electron tions, as compiled in KH13 and the references therein. density, have been directly measured with Chandra ob- Since the other 15 sources do not have reliable dynam- servations. We note that NGC 821, NGC 1553, and ical mass measurements, we estimate their BH masses NGC 5128 are excluded from our sample because the using the empirical relation given in Eq. (7) of KH13, electron densities for the first two sources in the literature are not observed values but just assumed to be 4.38±0.29 +0.37 8 σ −3 M• = 3.09 × 10 M , (1) ne ' 0.1 cm , and the density for the third source −0.33 200 km s−1 is not measured by Chandra observations but XMM- where σ is the central velocity dispersion of the bulge Newton, which has a lower spatial resolution of ' 6000. In addition, we collect some interesting sources whose stars, taken from the Hyperleda database (Paturel et al. Bondi scales are reported to be resolved. One is NGC 2003). The mean error of the BH masses for the sample 224 (M31), which is the nearest major galaxy to us and is ∆M•/M• = 0.26. contains an SMBH in the center, and detailed Chan- dra observations are available (Dosaj et al. 2002; Gar- 1 Note that NGC 4778, the brightest galaxy in HCG 62, is excluded cia et al. 2005). The others are nearby galaxies that because the BH mass measurement is not based on reliable meth- have been extensively studied: NGC 1332 (Humphrey ods (see §2.2). Universal phase transition of AGNs 3

2.3. Bolometric Luminosity (Lbol) diffuse emission that has a peak at . 1 keV. Therefore, Ideally, the AGN bolometric luminosities should be contamination from the nuclear component is not a se- measured directly from their broadband spectral energy rious issue for obtaining the temperature of diffuse hot distributions (SEDs). In the case of Sgr A∗, an almost gas surrounding the SMBH. complete nuclear SED is available from the radio to X- The spatially resolved region containing diffuse gas is ray band (Narayan et al. 1998), integration of which divided into radial annuli with sufficient counts for spec- 36 −1 tral fitting (see a typical example in Wong et al. 2014). yields Lbol = 1.2 × 10 erg s . Unfortunately, complete broadband nuclear SEDs are The MEKAL (Mewe et al. 1985) and APEC (Smith et al. not available for most of our sources. One of the most 2001) models are commonly used for spectral fitting, and both are implemented in the standard X-ray fitting secured indirect methods is to estimate Lbol from the 2– tool XSPEC (Arnaud 1996). The model fit yields the 10 keV absorption-corrected X-ray luminosity (L2−10) temperature T from the peak of the thermal bump at using a bolometric correction of CX ≡ Lbol/L2−10 = 15.8 (Ho 2009). Since our targets are very faint, the X- E = 0.3 − 1 keV, and the electron density ne from the ray data must be obtained with spatial resolution high spectral normalization, producing radial profiles of the enough to distinguish the nuclear source from the host two quantities. The literature does not always provide galaxy. We collect the X-ray luminosities from Chandra the uncertainties for estimating T and ne. In such cases, observations achieving the best spatial resolution in the we adopt the mean error of the sources with robust error X-ray band. Note that there are six sources that have no estimates; ∆T/T = 0.10 and ∆ne/ne = 0.43. X-ray core coincident with the assumed optical nucleus, In Table1, we list the temperatures and electron den- have no power-law component in their nuclear spectra, sities measured either at r = RB or at the inner-most or have a significant level of contamination from X-ray radii that the Chandra observations can reach if the binaries (see Table1). We treat their nuclear X-ray Bondi radii are not resolved. As exceptions, however, luminosities as upper limits in Figure1. we present the electron densities extrapolated to the It is quite difficult to evaluate the uncertainties of the Bondi radii of seven unresolved sources studied in R13 estimated bolometric luminosities. A major concern is (NGC 507, NGC 1316, NGC 4696, NGC 5044, NGC the possible underestimation of intrinsic absorption of 5813, NGC 5846, and NGC 6166), for which the density the emergent X-ray radiation. Low-luminosity AGNs, profiles are fitted with three different models (a power- however, generally suffer from very minimal intrinsic ab- law model continuing a steep density gradient to RB, a sorption (Ho 2008), and this effect typically introduces β-model flattening to a constant, and a shallow Sérsic an uncertainty only at the level of a factor of 2 to our profile with n = 4) with errors arising from the different luminosity estimates (Ho 2009). Accordingly, we assign model assumptions. errors of 50% to Lbol. 2.5. Bondi Radius (RB) and Accretion Rate (M˙ B) The Eddington ratio is defined by λEdd ≡ Lbol/LEdd, 38 −1 where LEdd = 1.26×10 (M•/M ) erg s . We estimate The Bondi radius is the characteristic radius within the uncertainty of λEdd by propagating the errors on which the gravitational force by the central BH domi- Lbol and M•. nates the thermal pressure gradient force of the gas surrounding the BH (Bondi & Hoyle 1944; Bondi 1952),

2.4. Temperature (T ) and Electron Density (ne) defined by Both the temperature and electron density are key pa- GM R ≡ • ' 1.68 M T −1 pc, (2) rameters needed to estimate the Bondi radius and the B 2 8 keV cs Bondi accretion rate. For each object in our sample, 8 we select the deepest Chandra observations from the where M8 = M•/(10 M ) and TkeV = kT/keV. We literature to obtain T and ne of the diffuse hot gas in emphasize that our definition of RB is smaller by a fac- the nuclear region. Although the details of the analysis tor of 2 compared to the one in Pellegrini(2005), which differ from study to study, the general methodology to uses the diameter instead of the radius. Note that the estimate these two quantities follows a similar manner. uncertainties in estimating the Bondi radius are mainly Whenever possible, the diffuse component is extracted due to the uncertainties from the BH mass measure- after removing point sources arising from X-ray binaries ments, and the mean uncertainty is ∆RB/RB ≈ 0.3. and the central nuclear component. Since the nuclear We divide the sample into two groups based on component usually has a characteristic power-law spec- whether or not the angular size of the Bondi radius trum in the hard X-ray band above ∼ 1 keV, this fact (θB ≡ RB/D) within the estimated uncertainty is larger enables us to distinguish the nuclear emission from the than 0.500, which is the best achievable resolution by 4 Inayoshi et al.

Chandra. Based on this criterion, we define sources hav- information is integrated. To derive the density from 00 ing sufficiently large Bondi radii (θB ≥ 0.5 ) as resolved the emission measure, we assume that the volume con- (see the first 11 sources labeled as “Bondi = Y” in Ta- tribution to the emission measure is inside a radius of ble1). Note that the Bondi radius of the BH in NGC 100 and that the slope of the density profile is β = −1. 4594 might be marginally resolved within the error of The choice of radius is quite arbitrary, but the estimated −3 the BH mass measurement, but the spatial resolution of value for NGC 315 (ne = 0.36 cm ) agrees well with the its Chandra observation reaches only ' 2.500 (Pellegrini density reported by other studies (Worrall et al. 2003). et al. 2003). Thus, we consider this source to be unresolved (see the remaining 20 sources labeled as “Bondi 3. RESULTS = N” in Table1). Figure1 presents the relation between the luminosity The Bondi accretion rate is defined as the mass inflow and Bondi accretion rate for the samples whose Bondi rate through the Bondi radius, scales are resolved (blue) and unresolved (green). The luminosity and accretion rate are both normalized by G2M 2 ˙ • their Eddington values, as BH accretion systems are MB ≡ 4πq(γ)ρ 3 cs characterized by these two dimensionless quantities, not −6 2 −3/2 −1 ' 7.0 × 10 ne,0.1M8 TkeV M yr , (3) by the actual physical scales. Assuming that the BH feeding rate is equal to the Bondi accretion rate, as in the where q(γ) = 1/4 for γ = 5/3 and ne,0.1 = original advection-dominated accretion flow model (gray −3 ne/(0.1 cm ). The Bondi accretion rate normalized thin dashed line), the radiative luminosity is expected to by the Eddington accretion rate is given as be much higher than the observed values by several orders of magnitude. While the luminosity discrepancy is M˙ B −6 −3/2 large at lower Bondi rates (M˙ 10−3 M˙ ), the lumi- m˙ B ≡ ' 3.1 × 10 ne,0.1M8TkeV . (4) B Edd M˙ −3 Edd nosity increases dramatically at M˙ B ' 10 M˙ Edd and 00 catches up to the lowest levels of the activity observed in For the unresolved sources (θB < 0.5 ), their Bondi ac- low-luminosity AGNs (L ≈ 10−4 L ; orange sym- cretion rates are estimated by extrapolating the elec- bol Edd bols). Even in the face of uncertainties of electron den- tron densities at RB from the measured values, namely 00 sity measurements for several unresolved objects, the by multiplying by a factor of f (≡ 0.5 /θB > 1) with −1 overall trend covering many orders of magnitude in the an assumed density profile of ne ∝ r . The choice of phase diagram of Figure1 is seen robustly. the slope β is consistent with the profiles for two well- resolved objects: β = −1.0 ± 0.2 for NGC 4486 (M87) Figure2 shows the two-dimensional distribution of the (Russell et al. 2015) and β = −1.05±0.25 for NGC 3115 density distribution of rotating accretion flows onto an SMBH, obtained from two-dimensional hydrodynamical (Wong et al. 2014). Note that the extrapolation has been taken into account for the seven unresolved sources simulations with radiative processes and BH feedback studied by R13. (Inayoshi et al. 2019), where feedback is modeled with the prescription obtained by Sazonov et al.(2005) with- For NGC 1291 and NGC 4594, the extrapolated val- out injecting disk winds (see also Yuan et al. 2018)2. ues of electron density are treated as lower limits because their nuclear regions are poorly resolved (θ ' 2.900 Unlike previous studies that have investigated accretion flows on smaller scales assuming a compact and gravita- and 2.500, respectively; Irwin et al. 2002; Pellegrini et al. tionally bound disk as the initial state, these numerical 2003). For NGC 221 (M 32), the electron density is the volume-averaged value within 3000 and is thus treated simulations focus on the accretion dynamics at larger as a lower limit (Ho et al. 2003). For convenience, we scales covering the Bondi radius (dashed circles), where the gas is weakly bound initially. This type of simu- list the “uncorrected” Bondi accretion rate and the cor- lation allows us to set plausible initial and boundary rection factor separately in Table1, but show the “corrected” Bondi accretion rate in Figure1. It is also worth conditions, which can be directly measured in the nuclear regions surrounding quiescent SMBHs. For the noting that removing these three galaxies and five ob- lower density case with M˙ /M˙ 10−3 (left panel), jects with no detection of nuclear X-ray emission from B Edd the unresolved sample, all of the other 12 unresolved radiative cooling is inefficient and viscous energy dis- sipation heats the gas, leading to a hot and turbulent sources are located in the transition region in Figure1 within the errors. Younes et al.(2011) did not provide ne for the four 2 Several previous works proposed the importance of mass loss via low-luminosity sources included in their study, only es- winds to explain the nature of low-luminosity AGNs (e.g., Bland- timates of the emission measure for which the volume ford & Begelman 1999; Li et al. 2013; Bu et al. 2016) Universal phase transition of AGNs 5 1

10-2

Thin Disk Cold 10-4 ・ B

Edd ・ IC 1459 M・ = M IC 4296

NGC 4486

bol -6 10 (M87) L / NGC NGC NGC 4374 4552 1407 NGC 3115 -8 SgrA* NGC 10 4472 RB - resolved NGC 224 NGC (M31) 4649 RB - unresolved -10 Luminous AGNs 10 Adiabatic Flow Transition

10-6 10-5 10-4 10-3 10-2 10-1 ・・ M B / MEdd Figure 1. Radiative luminosity and Bondi accretion rate for nearby quiescent and low-luminosity accreting SMBHs whose Bondi scales are resolved (blue) and unresolved (green) by Chandra observations. Numerical simulation results are shown by −3 the red solid and dashed curves. For lower values of M˙ B/M˙ Edd (< 10 ), the radiative luminosity is significantly reduced from that for M˙ • = M˙ B due to the suppression of the BH feeding by turbulent gas motion in the disk (red region). The transition −4 −3 from RIAFs to cold accretion disks occurs at 4 × 10 < M˙ B/M˙ Edd < 2 × 10 (yellow shaded region). After the transition, the cold thin disk tends to be unstable to its self-gravity (magenta region). These observational and theoretical results bridge the gap between the dormant SMBH population and the faint end of the Seyfert galaxy population in the local universe. Orange circles and squares represent AGNs from the Swift/BAT AGN and PG quasar catalog, respectively. accretion flow through a geometrically thick disk. Im- where the radiative efficiency model for a hot accretion portantly, the inflow rate decreases toward the center disk is adopted (Ryan et al. 2017; see also Figure 2 in as M˙ inf (r) ' M˙ B · (r/RB) due to turbulent motion and Inayoshi et al. 2019). The level of radiative luminosity results in a low value onto the central BH (Quataert & agrees well with the observational results for SMBHs ac- −3 Gruzinov 2000; Inayoshi et al. 2018). In such a hot accreting at M˙ B/M˙ Edd 10 , as shown in Figure1 (red cretion flow, thermal conductivity of electrons, instead region). Note that the width of the red region reflects of turbulence, transports energy outward at the inner the uncertainties of the conductivity suppression factor, −2 region of r < 10 RB. Since suppression of the accre- which has a range of 0.03 ≤ fc ≤ 0.3. ˙ ˙ −3 tion by turbulence ceases in the inner region, the final By contrast, for the highest rate of MB/MEdd & 10 BH feeding rate is the accreting gas collapses to the midplane and forms !3/5 a geometrically thin disk because of radiative cooling M˙ M˙ • ' 1.5 × 10−6 T −4/5α0.37f 2/5 B , (right panel). Since the entropy generated by viscos- 7 −2 −1 −3 M˙ Edd 10 M˙ Edd ity is radiated away, turbulent gas motion ceases, and (5) the gas accretion rate through the disk becomes as high 7 where T7 ≡ T/10 K is the gas temperature, α−2 ≡ as M˙ • ' M˙ B. This cooling transition dramatically in- α/0.01 is the strength of viscosity, and f−1 ≡ fc/0.1 creases the luminosity by several orders of magnitude is the conductivity suppression factor due to magnetic (Figure1, yellow region). We note that this critical rate fields (Narayan & Medvedev 2001). In other words, the required for the cooling transition around the BH in- ˙ −3 ˙ BH feeding rate is reduced by 2 − 3 orders of magnitude fluence radius, MB & 10 MEdd (Gaspari et al. 2015; from the Bondi accretion rate. As a result, the radiative Inayoshi et al. 2019) is 1 − 2 orders of magnitude lower luminosity reduces to than the critical rate for a compact accretion disk near ! the BH (Yuan 2001, 2003), the latter of which would L M˙ bol ' 3 × 10−8 T −4/3α0.63f 0.68 B , (6) 7 −2 −1 −3 LEdd 10 M˙ Edd 6 Inayoshi et al.

・・ 10-3 ・・ ~ 10-3 and host galaxy components, i.e., star-formation activ- M B / MEdd << M B / MEdd - ity (Ichikawa et al. 2019). The UV/optical-selected sample is based on the Palomar-Green (PG) quasar sur- vey (Boroson & Green 1992), and its star formation activity has also been investigated with decomposition of its IR SEDs (Shangguan et al. 2018). The rich multi- wavelength data sets of these objects enable us to estimate the bolometric luminosity and BH mass (i.e., Eddington luminosity and accretion rate), and also to study the statistical properties of their star formation rates (SFRs) on the host galaxy scales with the empirical relation (Kennicutt 1998), −1 LFIR SFRglobal ' 0.45 M yr . (7) 1043 erg s−1 We restrict the sample to bright objects with far-infrared 10 (FIR) luminosities LFIR > 10 L because the empirical relation is valid for star-forming galaxies and the IR 0.3 5.5 100 0.3 1700 107 SED fitting model is based on SED templates of galaxies density (n e /ne,∞ ) with FIR luminosities above this threshold. We note that because of their large distances, the Figure 2. Density distribution of accretion flows onto an Bondi radii for the Swift/BAT and PG quasar sam- SMBH accreting at two different rates, from the edge-on view ples are not resolved even with Chandra, and thus (i.e., the horizontal and vertical axes are parallel and per- pendicular to the equator of the accretion flow). For lower their Bondi accretion rates are not estimated properly. −3 For the purpose of illustration in Figure1, we replace accretion rates (M˙ B 10 M˙ Edd; left), the RIAF is kept adiabatic and highly turbulent, forming a geometrically thick the Bondi accretion rates on the horizontal axis with disk inside the Bondi radius (dashed line). For the higher 0.01 SFRglobal for all the 426 luminous AGNs. This fac- ˙ −3 ˙ accretion rates (MB ' 10 MEdd; right), the accreting gas tor of 0.01 is arbitrary but may be of a similar order of begins to cool, collapses toward the equator, and forms a magnitude as the Bondi accretion rate for the following cold, geometrically thin disk. reasons. First, local Seyfert galaxies with high-angular resolution (∼ 0.4 − 0.800) mid-infrared spectroscopy in- dicate that the SFRs in the nuclear regions (< 70 pc) be more relevant to the disk state transition for X-ray are, on average, 5 times lower than those measured in binaries. the circumnuclear regions on kiloparsec scales (Esquej 4. DISCUSSION et al. 2014). Second, numerical simulations studying the structure of circumnuclear disks suggest that strong 4.1. Relationship with local bright AGNs star formation activity leads to turbulence in the disk, Among the sample, 15 objects (five resolved and 10 and turbulent viscosity efficiently induces mass accre- unresolved ones) are located in the transition region of tion onto the Bondi scale at a rate of 10% of the nu- −7 Figure1 and are as luminous as 10 < Lbol/LEdd < clear SFR (Wada & Norman 2002; Inayoshi et al. 2019). −4 10 . Those SMBHs bridge the gap between the adia- With these simple assumptions, the bright AGN popu- batic accretion flows and the faint end of Seyfert galaxies lation seems to follow the red curve in Figure1, which in the local universe (Ichikawa et al. 2019). assumes M˙ • ' M˙ B with the radiative efficiency of a ge- In order to extend the state transition diagram to the ometrically thin accretion disk. brighter end, we consider two well-defined and comple- ∗ mentary samples of local AGNs selected by X-rays and 4.2. Past Activity of Sgr A UV/optical surveys (orange circles and squares, respec- Sgr A∗ is the SMBH in the center of the Milky Way 6 tively). The X-ray-selected sample was derived from with a mass of M• ' 4.4 × 10 M , whose activity is the all-sky 70-month Swift/Burst Alert Telescope (BAT; known to be very quiescent at present (Lbol/LEdd . Baumgartner et al. 2013), whose 3 − 500 µm infrared 10−8). However, several lines of observational evidence (IR) properties have been investigated (Meléndezet al. suggest that its past AGN activity was higher (Kaifu 2014; Ichikawa et al. 2017) and quantified systemati- et al. 1972; Koyama et al. 1996; Bland-Hawthorn & Co- cally from decomposition of the IR SED into the AGN hen 2003; Totani 2006; Ryu et al. 2013). One of the Universal phase transition of AGNs 7 most striking clues is the discovery of the Fermi bub- tationally unstable after the transition is indicated with bles, which are expanding above and below the Galactic the magenta region in Figure1(Inayoshi et al. 2019, see plane with an age of roughly a few Myr (Su et al. 2010). also Menou & Quataert 2001). This process naturally A short episode of AGN activity lasting ∼ O(Myr) and explains the existence of young massive stars in a thin injecting a total energy of order ∼ 1055 erg are required stellar disk in the Galactic center (Levin & Beloborodov to create the bubbles. This level of energy injection 2003). Major episodes of BH accretion with & 0.01 LEdd is achieved by assuming that the past AGN luminos- and star formation would likely blow away the accreting 40 −1 −4 ∗ ity was LX ' 10 erg s (Lbol/LEdd ' 3 × 10 , gas and quench the activity of Sgr A . adopting a typical value of the bolometric correction for low-luminosity AGNs). The Bondi accretion rate required to explain the energy output is on the order of M˙ ≈ 10−3 M˙ (≈ 10−4 M yr−1), which is pre- B Edd ACKNOWLEDGMENTS sumably the accretion rate through the cold disk (i.e., M˙ d ∼ M˙ B). In fact, the accretion rate is high enough We greatly thank Jeremiah P. Ostriker for use- for the disk to fragment into clumps and form stars by ful discussions. This work is partially supported by a spiral-mode gravitational instability, which is charac- the National Science Foundation of China (11721303, terized by the Toomre parameter (Toomre 1964): 11991052, 11950410493), the National Key R&D Pro- gram of China (2016YFA0400702), the Program for Es- 3α c3 Q ' eff s tablishing a Consortium for the Development of Human GM˙ d Resources in Science and Technology, Japan Science !−1 and Technology Agency, and the Japan Society for the α T 3/2 M˙ ' 0.83 eff c d 1, (8) Promotion of Science (JSPS) KAKENHI (18K13584). 0.03 300 K ˙ . MB Numerical computations were carried out with High- performance Computing Platform of Peking University where αeff (& 0.01) is the effective viscous parameter 1/2 and Cray XC50 at the Center for Computational As- caused by the spiral arms and cs (∝ Tc ) is the sound trophysics of the National Astronomical Observatory of speed of cold gas with a temperature of Tc. The pa- Japan. rameter range where the accretion disk would be gravi-

REFERENCES

Arnaud, K. A. 1996, in Astronomical Society of the Pacific Boroson, T. A., & Green, R. F. 1992, ApJS, 80, 109 Conference Series, Vol. 101, Astronomical Data Analysis Bu, D.-F., Yuan, F., Gan, Z.-M., & Yang X.-H., 2016, ApJ, Software and Systems V, ed. G. H. Jacoby & J. Barnes, 818, 83 17 Cappellari, M., Verolme, E. K., van der Marel, R. P., et al. Baganoff, F. K., Maeda, Y., Morris, M., et al. 2003, ApJ, 2002, ApJ, 578, 787 591, 891 Dexter, J., McKinney, J. C., & Agol, E. 2012, MNRAS, Baumgartner, W. H., Tueller, J., Markwardt, C. B., et al. 421, 1517 2013, ApJS, 207, 19 Dosaj, A., Garcia, M., Forman, W., et al. 2002, in Bender, R., Kormendy, J., Bower, G., et al. 2005, ApJ, 631, Astronomical Society of the Pacific Conference Series, 280 Vol. 262, The High Energy Universe at Sharp Focus: Blakeslee, J. P., Cantiello, M., Mei, S., et al. 2010, ApJ, Chandra Science, ed. E. M. Schlegel & S. D. Vrtilek, 147 724, 657 Emsellem, E., Dejonghe, H., & Bacon, R. 1999, MNRAS, Blakeslee, J. P., Jordán,A., Mei, S., et al. 2009, ApJ, 694, 303, 495 556 Esquej, P., Alonso-Herrero, A., González-Mart´ın,O., et al. Blakeslee, J. P., Vazdekis, A., & Ajhar, E. A. 2001, 2014, ApJ, 780, 86 MNRAS, 320, 193 Event Horizon Telescope Collaboration, Akiyama, K., Blandford, R. D., & Begelman, M. C., 1999, MNRAS, 303, Alberdi, A., et al. 2019a, ApJL, 875, L1 L1 Event Horizon Telescope Collaboration, Akiyama, K., Bland-Hawthorn, J., & Cohen, M. 2003, ApJ, 582, 246 Alberdi, A., et al. 2019b, ApJL, 875, L4 Bondi, H. 1952, MNRAS, 112, 195 Fabbiano, G., Elvis, M., Markoff, S., et al. 2003, ApJ, 588, Bondi, H., & Hoyle, F. 1944, MNRAS, 104, 273 175 8 Inayoshi et al.

Ferrarese, L., Ford, H. C., & Jaffe, W. 1996, ApJ, 470, 444 Meléndez,M., Mushotzky, R. F., Shimizu, T. T., Barger, Garcia, M. R., Williams, B. F., Yuan, F., et al. 2005, ApJ, A. J., & Cowie, L. L. 2014, ApJ, 794, 152 632, 1042 Menou, K., & Quataert, E. 2001, ApJ, 552, 204 Gaspari, M., Brighenti, F., & Temi, P. 2015, A&A, 579, A62 Mewe, R., Gronenschild, E. H. B. M., & van den Oord, Gebhardt, K., Adams, J., Richstone, D., et al. 2011, ApJ, G. H. J. 1985, A&AS, 62, 197 729, 119 Mo´scibrodzka, M., Falcke, H., & Shiokawa, H. 2016, A&A, Genzel, R., Eisenhauer, F., & Gillessen, S. 2010, Reviews of 586, A38 Modern Physics, 82, 3121 Narayan, R., Igumenshchev, I. V., & Abramowicz, M. A. González-Mart´ın,O., Masegosa, J., Márquez,I., Guerrero, 2000, ApJ, 539, 798 M. A., & Dultzin-Hacyan, D. 2006, A&A, 460, 45 Narayan, R., Mahadevan, R., Grindlay, J. E., Popham, Ho, L. C. 2008, ARA&A, 46, 475 R. G., & Gammie, C. 1998, ApJ, 492, 554 Ho, L. C. 2009, ApJ, 699, 626 Narayan, R., Medvedev, M. V., 2001, MNRAS, 562, L129 Ho, L. C., Terashima, Y., & Ulvestad, J. S. 2003, ApJ, 589, Narayan, R., & Yi, I. 1994, ApJL, 428, L13 783 Nowak, N., Saglia, R. P., Thomas, J., et al. 2008, MNRAS, Houghton, R. C. W., Magorrian, J., Sarzi, M., et al. 2006, 391, 1629 MNRAS, 367, 2 Paggi, A., Fabbiano, G., Kim, D.-W., et al. 2014, ApJ, 787, Humphrey, P. J., & Buote, D. A. 2004, ApJ, 612, 848 134 Humphrey, P. J., Buote, D. A., Brighenti, F., Gebhardt, Paturel, G., Petit, C., Prugniel, P., et al. 2003, A&A, 412, K., & Mathews, W. G. 2009, ApJ, 703, 1257 45 Humphrey, P. J., Buote, D. A., Gastaldello, F., et al. 2006, Pellegrini, S. 2005, ApJ, 624, 155 ApJ, 646, 899 Pellegrini, S., Baldi, A., Fabbiano, G., & Kim, D.-W. 2003, Ichikawa, K., Ricci, C., Ueda, Y., et al. 2017, ApJ, 835, 74 ApJ, 597, 175 Ichikawa, K., Ricci, C., Ueda, Y., et al. 2019, ApJ, 870, 31 Quataert, E., & Gruzinov, A. 2000, ApJ, 539, 809 Ichimaru, S. 1977, ApJ, 214, 840 Rusli, S. P., Thomas, J., Erwin, P., et al. 2011, MNRAS, Inayoshi, K., Ichikawa, K., Ostriker, J. P., & Kuiper, R. 410, 1223 2019, MNRAS, 486, 5377 Rusli, S. P., Thomas, J., Saglia, R. P., et al. 2013, AJ, 146, Inayoshi, K., Ostriker, J. P., Haiman, Z., & Kuiper, R. 45 2018, MNRAS, 476, 1412 Russell, H. R., Fabian, A. C., McNamara, B. R., & Irwin, J. A., Sarazin, C. L., & Bregman, J. N. 2002, ApJ, Broderick, A. E. 2015, MNRAS, 451, 588 570, 152 Russell, H. R., McNamara, B. R., Edge, A. C., et al. 2013, Jardel, J. R., Gebhardt, K., Shen, J., et al. 2011, ApJ, 739, MNRAS, 432, 530 21 Ryan, B. R., Ressler, S. M., Dolence, J. C., et al. 2017, Kaifu, N., Kato, T., & Iguchi, T. 1972, Nature Physical ApJL, 844, L24 Science, 238, 105 Kennicutt, Robert C., J. 1998, ARA&A, 36, 189 Ryan, B. R., Ressler, S. M., Dolence, J. C., Gammie, C., & Kormendy, J., & Ho, L. C. 2013, ARA&A, 51, 511 Quataert, E. 2018, ApJ, 864, 126 Koyama, K., Maeda, Y., Sonobe, T., et al. 1996, PASJ, 48, Ryu, S. G., Nobukawa, M., Nakashima, S., et al. 2013, 249 PASJ, 65, 33 Lanz, L., Jones, C., Forman, W. R., et al. 2010, ApJ, 721, Sazonov S. Y., Ostriker J. P., Ciotti L., Sunyaev R. A., 1702 2005, MNRAS, 358, 168 Levin, Y., & Beloborodov, A. M. 2003, ApJL, 590, L33 Schulze, A., & Gebhardt, K. 2011, ApJ, 729, 21 Li, J., Ostriker, J. P., & Sunyaev, R., 2013, MNRAS, 767, Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 105 Shangguan, J., Ho, L. C., & Xie, Y. 2018, ApJ, 854, 158 Loewenstein, M., Mushotzky, R. F., Angelini, L., Arnaud, Shen, J., & Gebhardt, K. 2010, ApJ, 711, 484 K. A., & Quataert, E. 2001, ApJL, 555, L21 Smith, R. K., Brickhouse, N. S., Liedahl, D. A., & Mathews, W. G., & Brighenti, F. 2003, ARA&A, 41, 191 Raymond, J. C. 2001, ApJL, 556, L91 McCourt, M., Parrish, I. J., Sharma, P., & Quataert, E. Soria, R., Fabbiano, G., Graham, A. W., et al. 2006, ApJ, 2011, MNRAS, 413, 1295 640, 126 McQuinn, K. B. W., Skillman, E. D., Dolphin, A. E., Berg, Stone, J. M., Pringle, J. E., & Begelman, M. C. 1999, D., & Kennicutt, R. 2017, AJ, 154, 51 MNRAS, 310, 1002 Universal phase transition of AGNs 9

Su, M., Slatyer, T. R., & Finkbeiner, D. P. 2010, ApJ, 724, Wong, K.-W., Irwin, J. A., Shcherbakov, R. V., et al. 2014, 1044 ApJ, 780, 9 Worrall, D. M., Birkinshaw, M., & Hardcastle, M. J. 2003, Tonry, J. L., Dressler, A., Blakeslee, J. P., et al. 2001, ApJ, MNRAS, 343, L73 546, 681 Xie, F.-G., & Yuan, F. 2012, MNRAS, 427, 1580 Toomre, A. 1964, ApJ, 139, 1217 Younes, G., Porquet, D., Sabra, B., & Reeves, J. N. 2011, A&A, 530, A149 Totani, T. 2006, PASJ, 58, 965 Yuan, F. 2001, MNRAS, 324, 119 van den Bosch, R. C. E., & de Zeeuw, P. T. 2010, MNRAS, Yuan, F. 2003, ApJL, 594, L99 Yuan, F., & Narayan, R. 2014, ARA&A, 52, 529 401, 1770 Yuan F., Yoon D., Li Y.-P., Gan Z.-M., Ho L. C., Guo F., Wada, K., & Norman, C. A. 2002, ApJL, 566, L21 2018, MNRAS, 857, 121 Walsh, J. L., Barth, A. J., & Sarzi, M. 2010, ApJ, 721, 762 Zhang, Z.-L., & Xu, H.-G. 2004, ChJA&A, 4, 221 10 Inayoshi et al. 5, 42, 14, 23 ∗ Bondi Reference 4 N 28, 29, 30, 31 1 N 5, 18, 32, 32 2 N 1, 18, 14, 25 6 N 1, 18, 32, 32 . . . . 04 N 1, 18, 32, 32 35 N 33, 18, 34, 34 32 N 5, 37, 38, 39 89 N 12, 40, 16, 16 07 N 5, 41, 14, 39 23 N 24, 18, 32, 32 38 N 11 N 5, 18, 14, 16 93 N 12, 43, 44, 44 f ...... 8114 –06 – – Y – Y Y 1, 2, Y 3, 5, 4 6, 7, 8 5, 9, 12, 10, 13, 11 14, 15 025 17 171914 –18 –12 – Y14 – Y 12,24 – Y 6, 16, 16 12, 17, – Y 14, 11 5,08 – Y 18, 14, 15 12, 19, Y 20, 20 3 5, Y 21, 22, 22 24, 18, 23, 25 26, 26, 27, 27 06 – N 1, 18, 10, 15 06 3 015 – N 12, 35, 36, 15 07 1 14 1 0020 48 020 5 07 2 006 17 19 1 027 4 06 – N 5, 18, 15, 15 05 2 009 11 070506 –013 – – N – N 5, N 18, 15, 15 N 5, 18, 25, 15 5, 18, 14, 1, 15 18, 15, 15 ...... ) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 ( ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± B 3 9 . . θ 72 51 66 42 54 80 54 36 40 16 16 15 38 26 22 36 17 17 19 21 24 ...... 029 082 099 029 122 043 047 ...... 0104 . 24 0 23 0 31 0 23 0 26 0 24 0 ...... 352326 1 22 0 0 0 28 0 222424 1 23 0 28 0 22 0 0 1 30 0 27 0 27 0 23 0 23 0 23 0 25 0 22 0 22 0 23 0 282424 0 25 0 0 0 0 0 0 0 0 0 ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ± ± ± ± ± ± Edd ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± λ 62 18 83 13 58 55 ...... 11 32 90 14 28 92 41 93 39 01 43 44 45 54 07 36 11 25 76 65 90 40 09 18 94 7 7 6 6 7 6 ...... log 9 6 7 6 7 5 6 7 5 5 8 4 5 5 6 5 5 5 5 5 4 6 7 6 5 − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − < < < < < < 11 23 12 21 33 31 08 4 10 20 19 22 26 4 24 4 19 29 18 20 20 20 28 20 31 21 18 22 23 4 30 ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 B ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± m 1 9 5 5 . . . . 82 51 23 47 91 64 87 34 25 72 64 34 24 03 62 08 45 32 43 23 67 47 21 93 56 54 27 ...... 3 . . . . . 2 . 3 ...... 3 . 3 3 3 3 4 6 2 3 3 2 4 3 4 3 3 5 3 3 4 4 4 3 5 3 3 3 3 − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − continued 2 8 1 10 010 6 5 17 35 9 9 ...... 16 5 15 7 3 5 9 2 8 8 8 3 0 25 17 13 0 10 14 17 32 21 2 1 14 0 0 1 11 0 11 (pc) log ˙ ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± 3 3 B 8 5 3 6 1 ...... 11 86 62 R . . . 056 . Table 1 ) 3 04014 4 10 100 1013 41 46 53 028 0 055012 114 13 40 63 16 74 81 13 51 7 56 11 6 2 8 4 41 18 26 15 0 4 15 4 18 4 1 06 17 04 8 25 30 008 10 4 4 10 28 13357 33 29 28 ...... − 0 0 0 0 0 0 11 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± (cm 6 8 3 0 0 2 0 2 1 ...... 10 25 41 31 62 31 97 37 28 44 43 15 11 40 30 28 40 e ...... 072 069 439 019 n . . . . . Physical quantities of low-luminosity SMBHs in nearby galaxies. 04020 0 0 11010 0 09 0 0 28 0 20 26 11010 0 030 0 0 07030 0 0 08 0 020 0 05 0 007 1 06 1 09 0 04 0 020 1 04 1 10 1 34 0 06 0 04 0 030 0 030 1 010 0 010009 0 07 0 1 ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 (K) ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± T 50 37 37 80 30 91 56 55 34 60 80 67 62 58 65 60 41 10 ...... 780 340 350 250 560 490 343 580 330 640 310 330 378 ...... Table 1 )

28 0 18 0 0510 0 0 0713 0 0 04 1 02711 0 10 0 08 1 17 0 0 20 0 16 0 16 0 08 0 06 0 22 0 08 0 09 0 12 0 06 0 027 0 08 0 15 0 11 0 09 0 17 0 1111 0 12 0 1 ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 /M • ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± M 16 39 67 22 97 40 95 63 92 67 39 43 21 12 23 17 10 53 72 82 75 49 86 31 20 71 81 82 27 ...... 789 823 . . 3 9 0 9 3 8 0 9 4 9 8 9 2 7 5 8 7 8 5 8 0 8 2 8 29 8 8 . 5172 8 9 68 9 4692 9 9 95 8 66 9 85 9 36 8 61 7 54 8 ...... 54 8 08 8 87 8 ...... 784 8 . . . 8057 6 . (Mpc) log( 00828 6 . . 0 D ∗ Name NGC 224NGC 1407NGC 0 3115 29 NGC 4374 9 NGC 4472 18 16 NGC 221 0 NGC 4486NGC 4552 16 NGC 4649 15 IC 1459 16 IC 4296Sgr A 28 46 NGC 315 64 NGC 507 70 NGC 1291 9 NGC 1316 20 NGC 1332 22 NGC 1399 20 NGC 2681 17 NGC 4261 32 NGC 4278 16 NGC 4438 11 NGC 4594 9 NGC 4636 14 NGC 4696 35 NGC 4697 12 NGC 5005 20 NGC 5044 31 NGC 5813NGC 5846 32 NGC 6166 24 125 Universal phase transition of AGNs 11 T , 10 − 2 Chandra L Bondi Reference , • M f , D ) 00 ( B ); Column 8: Eddington ratio θ Edd ˙ M / B ˙ M ≡ Edd B λ m log B m (continued) (pc) log ˙ B R Table 1 ) 3 − (cm e n (K) . T 00 5 . 2 ' )

/M • M ); Column 9: angular size of the Bondi radius; Column 10: correction factor of density for unresolved sources; Column 11: Bondi Edd /L ): (1) NED, (2) Bender et al. 2005 , (3) Garcia et al. 2005 , (4) Dosaj et al. ,2002 (5) Tonry et al. 2001 , (6) Rusli et al. 2013 , (7) Zhang & e — Column 1: galaxy name; Column 2: distance; Column 3: BH mass; Column 4 and 5: electron temperature and density at the Bondi radius; n bol (Mpc) log( L Column 6: Bondi radius; Column 7: Bondi accretion rate normalized by the Eddington rate ( ˙ ( resolved source = Y (the first 11 sources) and Bondi unresolved source = N (others); Column 12: references for the data of ( & Xu 2004 , (8) 2010 ,Humphrey et (14) al. al. González-Mart´ınet 2006 ,2013 ,2006 , (15) (19) (9) RussellShenEmsellemBlakeslee et et & et al. al. Gebhardt al. ,1999 2013 ,2010 , Zeeuw 2001 , (10) (20) (16) 2010 , (25) Ho PaggiLoewenstein (30) Pellegrini 2009 , et etNowak2005 ,Ho (11) al. et al. (26) etRussell2014 , al. Genzel2001 , al. et (21) 2006 ,2008 , et2003 , (17) al. Cappellari (36) (41) al. (31) Gebhardt,2015 et2010 ,LanzFerrareseMathews (12) etof et al. (27) et &Blakeslee al. the al. Baganoff2002 , al. Brighenti et2011 , Bondi et1996 , (22) 2010 ,,2003 al. (18) al. radius (42) Fabbiano (37) (32) ,2009 Kormendy et2003 , ofJardelRusliYounes (13) & (28) et al. the et etWalsh Ho Blakeslee et al. 2003 , central al. al. et al. ,2011 (23) SMBH,2011 2011 , al. (33) Pellegrini in (38) (43) 2010 ,McCourt etHumphrey NGCSchulze (29) et al. & & 4594van al. 2003 , Buote Gebhardt might den2011 , (24) ,2004 be,2011 Bosch (34) (39) and marginally &IrwinHumphrey de (44) resolved et et within ( Soria al. al. the et,2002 2009 , error, al. (35) (40) 2006 ). butHoughton the Note et spatial that al. resolution the of angular its size observation reaches only D Note Name