arxiv: v2 [astro-ph.he] 7 Jan 2016

Transcription

1 Mon. Not. R. Aston. Soc.,?? (22) Pinted 26 July 28 (MN LATEX style file v2.2) Two-dimensional inflow-wind solution of black hole accetion with an evenly symmetic magnetic field axiv:5.2456v2 [asto-ph.he] 7 Jan 26 Amin Mosallanezhad,2, Defu Bu and Feng Yuan Shanghai Astonomical Obsevatoy, Chinese Academy of Sciences, 8 Nandan Road, Shanghai 23, China 2 Univesity of Chinese Academy of Sciences, 9A Yuquan Road, Beijing 49, China 26 July 28 INTRODUCTION Black hole accetion models can be dividedinto hot and cold types accoding to the tempeatue of the accetion flow. Since the pionee wok of Naayan & Yi (994) (see also Ichimau 977; Rees et al. 982), intensive studies have been pefomed on black hole hot accetion flows. It is now known thathot accetion flowis veycommon inthe univese,anging fom low-luminosity active galactic nuclei (LLAGNs), which is the majoity of neaby galaxies, to the quiescent and had states of black hole X-ay binaies. Yuan & Naayan (24) have ecently eviewed ou cuent undestanding of hot accetion flows, including its one-dimensional and multidimensional dynamics, jet fomation, adiation, and vaious astophysical applications. One of the most inteesting topics in ecent yeas in this field is wind (also often called outflow). Wind is not only a fundamental aspect of accetion flow, but also plays an impotant ole in AGN feedback, since wind can suppess the sta fomation and black hole gowth effectively (e.g., Ostike et al. 2). The study of wind was initiated by the hydodynamical and MHD numeical simulations of accetion flows (Stone, Pingle & Begelman 999; Igumenshchev & Abamowicz 999, 2; Stone & Pingle fyuan@shao.ac.cn (FY) ABSTRACT We solve the two-dimensional magnetohydodynamic (MHD) equations of black hole accetion with the pesence of magnetic field. The field includes a tubulent component, whose ole is epesented by the viscosity, and a lage-scale odeed component. The latte is futhe assumed to be evenly symmetic with the equatoial plane. The equations ae solved in the plane of a spheical coodinate by assuming timesteady and adially self-simila. An inflow-wind solution is found. Aound the equatoial plane, the gas is inflowing; while above and below the equatoial plane at a cetain citical angle, 47, the inflow changes its diection of adial motion and becomes wind. The diving foces ae analyzed and found to be the centifugal foce and the gadient of gas and magnetic pessue. The popeties of wind ae also calculated. The specific angula momentum of wind is found to be significantly lage than that of inflow, thus wind can tansfe angula momentum outwad. These analytical esults ae compaed to those obtained by the tajectoy analysis based on MHD numeical simulation data and good ageements ae found. Key wods: accetion, accetion discs black hole physics hydodynamics. 2). These simulations found that in contast to ou oiginal pictue, the mass accetion ate deceases with deceasing adius (see eview in Yuan et al. 22b). To explain such an astonishing esult, two models have been poposed, namely convection-dominated accetion flow (CDAF; Naayan et al. 2; Quataet & Guzinov 2) and adiabatic inflow-outflow solution (ADIOS; Blandfod & Begelman 999). By compaing the popeties of inflow and outflow and studying the convective stability of MHD accetion flow, Yuan et al. (22a; see also Li, Ostike & Sunyaev 23) show that stong wind must exist and it is the wind athe than convection that esults in the inwad decease of accetion ate. Most ecently, based on MHD simulation data, Yuan et al. (25) have convincingly show the existence of stong wind by obtaining the tajectoy of the vitual test paticles based on 3D GRMHD simulation data of black hole accetion 2. Moeove, the diving foces of wind have also been analyzed and the detailed popeties of winds have been calculated, Naayan et al. (22) also studied the outflow in hot accetion flow. While they also find the existence of wind, thei wind is much weake than that found by Yuan et al. 22a. The eason fo the discepancy was analyzed in Yuan et al. (25). 2 Bu et al. (23) show that the wind becomes much weake if the angula momentum of the accetion flow is vey low.

2 2 Mosallanezhad, Bu & Yuan such as the mass flux, angula distibution, teminal velocity, and fluxes of enegy and momentum. The theoetical study on the existence of stong wind has been confimed by the detection of emission lines fom the accetion flow aound the supemassive black hole in the Galactic cente, Sg A* (Wang et al. 23). As an application of the accetion wind theoy, moe ecently, the fomation of the Femi bubbles detected by the Femi telescope in the Galaxy has been successfully explained by the inteaction of wind launched fom the accetion flow aound Sg A* and the intestella medium (Mou et al. 24; 25). Compaed to othe theoetical models of the Femi bubbles, this accetion-wind model has two advantages. One is that the paametes of wind adopted in the model ae taken fom the MHD simulation of accetion flow. The second one is that the model can explain the new obsevational esults obtained by Suzaku and Chanda which ae in conflict with the pedictions of othe models. Almost all of the above-mentioned theoetical woks on wind ae numeical. Compaed to numeical simulations, analytical study has its advantage of moe clealy evealing some undelying physics. In this pape, using analytical method, we solve the axisymmetic steady solution of hot accetion flow with the pesence of magnetic field. Self-simila assumption is adopted as in many pevious woks. Special attention is paid to whethe the winds exist in the solution and what ae thei diving mechanisms if they exist. The aim is to compae with the numeical esults pesented in Yuan et al. (25) thus to impove ou undestanding of numeical studies. Thee have been many analytical woks based on the self-simila assumption. In the one-dimensional woks, the existence of wind is often an assumption (e.g. Blandfod & Begelman 999; Akizuki & Fukue 26; Abbassi et al. 28; Zhang & Dai 28; Bu et al. 29). Winds have been found in some two-dimensional solutions (e.g. Xu & Chen 997; Blandfod & Begelman 24; Xue & Wang 25; Tanaka & Menou 26; Jiao & Wu 2; Mosallanezhad et al. 24; Gu 25; Samadi & Abbassi 26 ) but not all (e.g., Naayan & Yi 995 and efeences theein). In all these woks magnetic field is not included. It is well known that magnetic field must be pesent and even plays an impotant ole in the dynamics of accetion flow, such as the angula momentum tansfe by MRI (Balbus & Hawley 998), the convective stability of accetion flow, which is elated with the wind poduction (Yuan et al. 22a; Naayan et al. 22), and the diving mechanism of winds (Yuan et al. 25). So it is useful to solve the accetion solutions by including magnetic field. This is the aim of the pesent pape. The oganization of the pape is as follows. In 2, the basic MHD equations and assumptions ae given. The numeical solutions of the equations ae descibed in 3. The esults ae summaized in 4. 2 BASIC MHD EQUATIONS The basic MHD equations descibing accetion flows ead: d + v =, dt () dv dt = Ψ p+ T (J B), c (2) de dt p d dt = fq+, (3) B (v B t = ) ηj, (4) B =. (5) Intheaboveequations,,v,p,e,T,Bandη denotethedensity, velocity, gas pessue, gas intenal enegy, viscous stess tenso, magnetic field, and magnetic diffusivity, espectively. We adopt an equation of state of ideal gas, p = (γ )e, whee γ is the specific heat atio. A spheical coodinate (,, φ) is adopted in the pesent wok. Ψ(= GM/) is the Newtonian gavitational potential, whee M is the mass of the cental black hole, G is the gavitational constant. J c/4π( B) is the electic cuent. We must include this dissipation tem in eq. (4), because othewise we can t obtain the steady solution due to the accumulation of magnetic flux in the accetion. Q + is the heating ate. f is the advection facto descibing the faction of the heating enegy which is stoed in the gas and advected into the black hole. In this pape, we set f =. The heating ate Q + can be decomposed into two components, Q + = Q vis +Q es, (6) with Q vis and Q es denoting the viscous heating and the magnetic field dissipation heating, espectively. Following Stone et al. (999), we assume that the only non-zeo component of the viscous tenso T is the azimuthal component, T φ = ν (v φ/), (7) whee ν is the kinematic viscosity. The viscous heating and the magnetic field dissipation heating can be descibed as, Q vis = T φ ( vφ ), (8) Q es = 4π c 2 ηj2. (9) Weadopttheusualαdesciptionofviscosityν = αp/(ω K), whee Ω K = (GM/ 3 ) /2 is the Kepleian angula velocity. In ode to satisfy the adial self-simila condition we assume the magnetic diffusivity η = η p/(ω K). Following the suggestion fom numeical MHD simulations, we decompose the magnetic field into a lage-scale component and a tubulent component. Both of them can tansfe the angula momentum and can be dissipated and poduce heat. The magnetic field B in equations (2), (4) and (5) coesponds to the fome. We descibe the effects of the tubulent component of the magnetic field in tansfeing the angula momentum and in dissipating the enegy though the usual α desciption. Specifically, the viscous foce T in equation (2) epesents angula momentum tansfe by the tubulent magnetic field; the viscous heating ate Q vis in equation (6) is coespondingly associated with the tubulent component of the magnetic field. We assume that the lage-scale magnetic field is evenly

3 Wind fom hot accetion flow 3 symmetic about the equatoial plane (see Figue fo the poloidal configuation of the magnetic field). This kind of symmety is widely adopted in the studies of accetion disks (e.g. Blandfod & Payne 982; Lovelace et al. 994; Cao 2; Li & Begelman 24) 3. Thus we have, The thee components of induction equation (4) can be expessed as: B t = [sin(vb v B ηj φ )], (2) B (,) = B (,π ), () B (,) = +B (,π ), () B φ (,) = B φ (,π ), (2) B t = [sin(vb v B ηj φ )], (22) A tooidal component of the magnetic field is geneated due to the shea of the accetion flow. Theefoe the adial and the tooidal components of the magnetic field should have opposite sign. In this pape, we only focus on the egion above the equatoial plane, we set the adial component of themagnetic field to be positive (B > ). Theefoe, the tooidal component of the magnetic field geneated by shea is negative (B φ < ). We solve fo the steady state, axisymmetic ( / t = / φ = ) solutions of eqs. -5. The continuity equation can be ewitten as: ( ) 2 v 2 + sin (sinv ) =, (3) The momentum equation (2) eads: [ [ v v + v v v + v ( ) ] v v v2 φ = GM p 2 4π (J B φ J φ B ), (4) ( ) ] v +v v2 φ cot = p 4π (J φb J B φ ), (5) [ v φ v + v v φ + v ] φ (v +v cot) = ( 3 ) T 3 φ 4π (JB J B ), (6) whee the cuent (J) eads: J = sin (B φsin), (7) J = (B φ), (8) J φ = [ ] (B ) B. (9) The equation of enegy is expessed as: ( e v + v f ) e p ( v + v ( T φ ( vφ ) + η 4π ) = ( J 2 +J 2 +Jφ 2 ) ). (2) 3 Shaham Abassi s goup is woking on the case of an oddly symmetic magnetic field at the same time with us. B φ t = (v φb v B φ )+ (v φb v B φ )+ (ηj) (ηj ). (23) Fo a steady state, the ight hand side of equations (2), (22) and (23) has to be equal to zeo. So we have sin(v B v B ηj φ ) = o a constant. In this pape, we assume it is equal to zeo: sin(v B v B ηj φ ) =. (24) Finally, the divegence-fee equation (5) can be witten as, ( 2 ) B 2 + sin (sinb ) =. (25) Simila to pevious woks (e.g., Naayan & Yi 995; Xue & Wang 25; Jiao & Wu 2), we seek the adially self-simila solutions in the following foms: (,) = () n, (26) GM v (,) = v() = ΩK()v(), (27) v (,) = Ω K()v (), (28) v φ (,) = Ω K()v φ (), (29) p(,) = p()gm n, (3) B (,) = b () GM (n/2) (/2), (3) B (,) = b () GM (n/2) (/2), (32) B φ (,) = b φ () GM (n/2) (/2). (33) Using the above self-simila assumptions, equations (3)-(6), (2), (23)-(25) can be witten as, ( ) dln 3 v + d 2 n v +v cot + dv =, (34) d [ ] dv 2 v2 +v d v2 vφ 2 = +(n+)p 4π (j b φ j φ b ), (35) [ ] 2 vv dv +v d v2 φcot = dp d 4π (j φb j b φ ), (36) [ ] 2 vv dv φ φ +v d +v v φ cot = 3 2 α(n 2)pv φ 4π (jb j b ), (37)

4 4 Mosallanezhad, Bu & Yuan Figue. Solid lines epesent the magnetic field line. Colos show the logaithm angula velocity (v φ /). It is clea that moving along the magnetic field lines stating fom the equatoial plane, the angula velocity deceases. Theefoe, if we assume that B above the equatoial plane is positive, shea of the flow will poduce a negative B φ. dp (nγ n )v p+v d γpv dln = ( d 9 f (γ ) 4 αpv2 φ + η ( j 2 +j 2 +j 2 ) ) φ, (38) 4π n 2 (vb db φ v φ b +ηj )+v φ d +b dv φ d v db φ d b dv φ d +η dj dη d +j =, (39) d whee db d 2 (n 3)b +b cot =. (4) η = η p, (4) j = db φ d +b φcot, (42) j = 2 (n )b φ, (43) j φ = η (vb v b ). (44) Equations (34)-(4) ae diffeential equations fo eight vaiables: v (), v (), v φ (), (), p(), b (), b () and b φ (). In this pape, we assume the accetion flow is evenly symmetic about the equatoial plane, then () = (π ), p() = p(π ), v () = v (π ), v φ () = v φ (π ) and v () = v (π ). At the equatoial plane, we have by symmety: = π 2 : b = v = dv d = dv φ d = d d = dp d = db d = db φ d =. (45) Wesetthedensityattheequatoial plane(π/2) =. Unde the above bounday conditions, equations (34)-(4) can be simplified into the following equations, dv d = ( n 3 ) v, (46) 2 2 v2 +vφ+(n+)p+ 2 ( ) 4π 2 (n )b2 φ + vb2 =, (47) η p v = 3α(n 2)p, (48) (nγ n )v = f (γ ) [ ( ( ) 2 ( ) )] αv2 φ + η vb 4π 2 (n )b φ +. (49) η p At the equatoial plane, we set the atio of the magnetic pessue to the gas pessue as, β = b2 8πp g, (5) π/2 β φ = b2 φ 8πp g. (5) In pinciple, to solve eqs. (34)-(4), in addition to the bounday condition at = 9, we should also supply the bounday condition at the otation axis, i.e., =. Mathematically, this is a two-point bounday poblem and it will assue that we have a well-behaved solution in the whole space. We have tied to obtain such solutions but failed. Theefoe, we only equie the solution satisfying the bounday condition at = 9. As we will explain in the next section, we think the solution obtained in this way should still be physically meaningful. π/2

5 Wind fom hot accetion flow 5 v/v K v/v K vφ/v K p v/cs, v/cs Mach v Mach v.3. Figue 2. Angula pofile of a vaiety of vaiables. The up-left plot shows that the adial velocity changes its sign at RESULTS We solve the above equations (34)-(4) numeically. We integate these equations fom the equatoial plane ( = 9 ) towads the otational axis ( = ). We adopt the values of β and β φ with efeence to the numeical simulation data of Yuan et al. (22a), but ou esults ae not sensitive to the exact values. As we have explained befoe, the tooidal component of the magnetic field, which is poduced by shea, should be negative. We do find that the tooidal component is negative at the beginning of integation. But at a cetain citical value of, denoted as B, we find that it begins to change its sign. At the egion of < B, the tooidal component of the magnetic field becomes positive thus the solution is no longe physical so we stop ou integation at B. We think ou solution in the egion of B < < 9 is still physical 4. This is because, on the one hand, the solution satisfies the equations and the bounday conditions at = 9. On the othe hand, the values of all physical quantities of the accetion flow at = B ae physical. Theefoe, we can easonably teat these values as bounday conditions at B. In fact, as we will descibe below, the main popeties of the solutions we obtain in this way ae in good ageement 4 The two-dimensional solution obtained by Jiao & Wu (2) is simila in the sense that they also have to stop thei integation at a cetain >. with those obtained in Yuan et al. (25) fom numeical simulations. Fom the bottom-ight plot of Figue 2, it seems that thee is a sonic point fo v /c s and the integation stops when the sonic point is appoached because of the singulaity thee. This is not tue. This is because, the value of v must be equal to zeo at the pole because of the axisymmetic bounday condition thee. Theefoe, v /c s should not pass though a sonic point. We note that Jiao & Wu (2) obtained a simila esult with us, i.e., thei solution also has to stop at a cetain angle. An inteesting explanation fo such a tuncation is given in thei wok. 3. Oveall inflow-wind stuctue of the solutions Figue shows the poloidal magnetic field configuation. We emphasize that while the even symmety of the field is ou assumption, the stength of the field is obtained by solving equations (34)-(4). Above the equatoial plane, the adial component of the magnetic field is positive (pointing to lage adii). So the tooidal component of the magnetic field, which is poduced by shea, is negative. We find that such kind of magnetic field configuation is simila to the time-aveaged field configuation in MHD numeical simulation in Yuan et al. (22a) and Yuan et al. (25). Figue 2 pesents the angula pofiles of fou quantities. The paametes we adopt ae α = η =.,γ = 5/3,β = 2 4,β φ = 2 3,n =.85 and f =. Fom the

6 6 Mosallanezhad, Bu & Yuan z/s / s Figue 3. Steam lines of the accetion flow, showing the inflow-wind stuctue. The length of aows indicates the absolute value of v + v. equatoial plane to the otation axis, the density and pessue decease, while v and v φ incease. This is consistent with that in Jiao & Wu (2). The value of v is negative and its absolute value inceases towads the pole. The Mach numbe v /c s inceases towads the pole, this is also consistent with that found in Jiao & Wu (2). Howeve, we note that beyond B, v have to decease towads the pole and eventually become zeo at pole due to the bounday condition thee. In the egion of 47 < < 9, the adial velocity is negative. This is the inflow egion. At the egion of < 47, the adial velocity becomes positive, i.e., wind is evident. Figue 3 shows the steam line of the flow. This esult is fully consistent with the esults obtained in Yuan et al. (25; see also Naayan et al. 22; Sadowski et al. 23). In those woks, they also found that inflow is pesent aound the equatoial plane of the accetion flow while wind is pesent in the pola egion. The bounday between inflow and wind obtained in the pesent pape is even quantitatively simila to that obtained in Yuan et al. (25). In both cases the bounday is located at about Angula momentum tansfe by wind The up-ight panel of Figue 2 shows the otation velocity of the accetion flow. We find that it significantly inceases fom the equatoial plane to the otation axis, which suggests that the specific angula momentum of wind may be lage than that of inflow. To examine this point in moe detail, we have calculated the mass flux-weighted specific angula momentum of both the wind and inflow (see Eqs. 8 & 9 in Yuan et al. (22a) fo the details of calculation). In unit of the Kepleian angula momentum l K, the mass-flux weighted angula momentum of wind and inflow is found to be l wind =.67l K, l inflow =.59l K, (52) espectively. This esult means that wind can tansfe the angula momentum of the accetion flow. This esult is again consistent with the numeical simulation esult pesented in Yuan et al. (22a). To compae the efficiency of angula momentum tansfe by outflow and by viscosity, we calculate the angula momentum flux tanspoted by outflow and viscosity as follows We find F Lout = 2π 2 π/2 B max(v,)(v φ )sind, (53) F Lvis = 2π 2 π/2 B T φ sind. (54) F Lout F Lvis =.26. (55) This means outflows do play an impotant ole in tansfeing angula momentum.

7 Wind fom hot accetion flow z/s 5 4 Outflow egion 3 2 Inflow egion F g F c -/ p gas -/(8π) B 2 F total / s Figue 4. The foce analysis at the wind egion to show the diving mechanism of wind. The length of the aows schematically denotes the magnitude of the foces while the diection fo the diection of foces. F g: gaviational foce; F c: centifugal foce; F total : the total foce. 3.3 Poloidal speed of wind We have also calculated the mass flux-weighted poloidal speed v p( v 2 +v 2 )of winds when theyae just launched. The esult is v p() =.8 v K(), (56) whee v K is the Kepleian velocity. This is consistent with that obtained in Yuan et al. (25), whee they found that v p.2v K(). 3.4 Benoulli paamete and teminal velocity of wind The Benoulli paamete is defined as Be() = v2 2 + γp (γ ) GM. (57) The fist, second and the thid tems on the ight-hand side of the equation coespond to the kinetic enegy, enthalpy and gavitational enegy, espectively. We have calculated the mass flux-weighted Benoulli paamete of the wind and obtained Be() =. v 2 K(). (58) The value of Be is positive, which means that the gas has enough enegy to ovecome the gavitational potential, to do wok to its suoundings, and to escape to infinity. As a compaison, the numeical simulation esult obtained in Yuanet al. (22a; eq. 9) is that Be()/v 2 K() (..2). So the pesent esult is again in good ageement with that obtained by numeical simulations 5. Now, as in Yuan et al. (22a), we ty to estimate the teminal poloidal velocity of wind based on Be(). To do this, we assume that, once launched, the viscous stess in the wind can be ignoed and Be() is conseved. When is lage enough, the last two tems in Equation (57) vanish. If the angula momentum is conseved (the wind is inviscid), the otational velocity is zeo at lage adius, so the teminal velocity is mainly the polodial velocity. Theefoe the teminal velocity is Combing with eq. (58), we have v p,tem v p 2Be() (59) v p,tem.47 v K(). (6) So the teminal poloidal velocity of wind oiginating fom adius is about half of the Kepleian velocity at thei oigin. As a compaison, Yuan et al. (22a) find v p,tem().5v K() while Yuan et al. (25) find v p,tem() (.2.4)v K(). 5 Yuan et al. (25) find that Be() is oughly constant of adius; but note that the constant Be is along the tajectoies of test paticles.

8 8 Mosallanezhad, Bu & Yuan 3.5 Analysis of foces diving the wind To study the diving mechanisms of wind, we have calculated the foces at the wind egion. Figue 4 shows the esult. We can see fom the figue that the dominant diving foces ae the centifugal foce, the gadient of the gas pessue and the magnetic pessue. The stength of the thee foces ae compaable. This esult is again fully consistent with that found in Yuan et al. (25), although hee we have an odeed lage-scale magnetic field while the field in Yuan et al. (25) is tangled. Jiao & Wu (2) found that the diving foces ae the gas pessue gadient and centifugal foce because they do not include magnetic field in thei model. Compaed to the Blandfod & Payne(982) model, it is simila in the sense that the centifugal foce is impotant; the diffeence is that hee the gas pessue and magnetic pessue gadient foces play an equally impotant ole in diving the winds. 4 SUMMARY AND DISCUSSION In this pape, we have solved the two-dimensional MHD solution of black hole accetion, with the inclusion of an evenly symmetic magnetic field (Fig. ). The steady state and selfsimila assumption in the adial diection ae adopted to simplify the calculation. One caveat in ou solution is that we fail to obtain the solution in the full space. Ou integation stats fom the equatoial plane but stops at a cetain B angle since the solution becomes unphysical beyond that angle. But we believe that the solution obtained is still physical since all the physical quantities fom = 9 to B ae physical. An inflow-outflow solution is found. Aound the equatoial plane, it is the inflow. Above and below the equatoial plane beyond a cetain angle, it is outflow (wind) (Fig. 3). This stuctue is same as that obtained in pevious numeical simulation woks (e.g., Yuan et al. 22; Naayan et al. 22; Sadowski et al. 23; Yuan et al. 25). We have also calculated some popeties of winds and found good consistency with those obtained in Yuan et al. (25) and Yuan et al. (22), which ae based on the MHD numeical simulation esults. Fo example, we find that the specific angula momentum of wind is lage than that of inflow thus wind can effectively tansfe angula momentum; the poloidal velocity of wind at both thei oigin and infinity calculated in the pesent analytical wok is also quite simila to those obtained by numeical simulations. At last, the diving foce of wind is analyzed and found to be the sum of centifugal foce, gadient of gas and magnetic pessue. This is again consistent with the esults obtained in Yuan et al. (25). ACKNOWLEDGMENTS We thank Zhao-Ming Gan and Mao-chun Wu fo useful discussions. This poject was suppoted in pat by the National Basic Reseach Pogam of China (973 Pogam, gant 24CB8458), the Stategic Pioity Reseach Pogam The Emegence of Cosmological Stuctues of CAS (gant XDB9), the Natual Science Foundation of China (gants 335 and 5735), and the CAS/SAFEA Intenational Patneship Pogam fo Ceative Reseach Teams. REFERENCES Abbassi, S., Ghanbai, J., Najja, S., 28, MNRAS, 388, 663 Akizuki, C., Fukue, J., 26, PASJ, 58, 469 Balbus, S. A., & Hawley, J. F., 998, Rev. Mod. Phys., 7, Blandfod, R. D., Begelman, M. C. 999, MNRAS, 33, L Blandfod, R. D., Begelman, M. C. 24, MNRAS, 349, 68 Blandfod, R., Payne D. G., 982, MNRAS, 99, 883 Bu, D., Yuan, F., Xie, F., 29, MNRAS, 392, 325 Bu, D., Yuan, F., Wu, M., Cuada, J., 23, MNRAS, 434, 692 Cao, Xinwu., 2, ApJ, 737, 94 Gu, W. M., 25, ApJ, 799, 7 Ichimau, S., 977, ApJ, 24, 84 Igumenshchev, I. V., Abamowicz M. A., 999, MNRAS, 33, 39 Igumenshchev, I. V., Abamowicz M. A., 2, ApJS, 3, 463 Jiao, C. L., Wu, X. B., 2, ApJ, 733, 2 Li, S. L., Begelman, M. C., 24, ApJ, 786, 6 Li, J., Ostike J., Sunyaev R., 23, ApJ, 767, 5 Lovelace, R. V. E., Romanova, M. M., Newman, W. I., 994, ApJ, 437, 36 Mosallanezhad A., Abbassi S., Beianvand N., 24, MN- RAS, 437, 32 Mou, G., Yuan, F., Bu, D., Sun, M., Su, M., 24, ApJ, 79, 9 Mou, G., Yuan, F., Gan, Z., Sun, M. 25, ApJ, axiv: Naayan, R., Igumenshchev, I. V., & Abamowicz, M. A. 2, ApJ, 593, 798 Naayan, R., Sadowski, A., Penna, R. F., et al. 22, MN- RAS, 426, 324 Naayan, R., & Yi, I., 994, ApJ,428, L3 Naayan, R., & Yi, I., 995, ApJ, 444, 238 Ostike, J. P., Choi E., Ciotti L., Novak G. S., Poga D., 2, ApJ, 722, 642 Quataet, E., & Guzinov, A. 2, ApJ, 539, 89 Rees M. J., Begelman M. C., Blandfod R. D., Phinney E. S., 982, Natue, 295, 7 Sadowski, A., Naayan, R., Penna, R., & Zhu, Y. 23, MNRAS, 436, 3856 Samadi M., Abbassi S., 26, MNRAS, 455, 338 Stone,J.M., Pingle J.E., Begelman M. C., 999, MNRAS, 3, 2 Stone, J. M., Pingle J. E., 2, MNRAS, 322, 46 Tanaka, T., & Menou, K. 26, ApJ, 649, 345 Wang, Q. D., Nowak, M. A, Makoff, S. B, et al., 23, Science, 34, 98 Xu, G., Chen, X., 997, ApJ, 489, L29 Xue, L., Wang, J.-C., 25, ApJ,623, 372 Yuan, F., Bu, D., Wu, M., 22a, ApJ, 76, 3 Yuan, F., Naayan, R., 24, ARA&A, 52, 529 Yuan, F., Wu M., Bu D., 22b, ApJ, 76, 29 Yuan, F., Gan, Z., Naayan, R., Sadowski, A., Bu, D., Bai, X. 25, ApJ, 84,

9 Zhang, D., Dai, Z. G., 28, MNRAS, 388, 49 Wind fom hot accetion flow 9