Issue 
A&A
Volume 587, March 2016



Article Number  A38  
Number of page(s)  7  
Section  Astrophysical processes  
DOI  https://doi.org/10.1051/00046361/201527878  
Published online  12 February 2016 
Limits on thickness and efficiency of Polish doughnuts in application to the ULX sources
^{1} N. Copernicus Astronomical Center, ul. Bartycka 18, 00716 Warszawa, Poland
email: maciek.wielgus@gmail.com
^{2} Institute of Micromechanics and Photonics, Warsaw University of Technology, ul. Św. A. Boboli 8, 02525 Warszawa, Poland
^{3} Institut d’Astrophysique de Paris, CNRS et Sorbonne Universités, UPMC Univ. Paris 06, UMR 7095, 98bis Bd Arago, 75014 Paris, France
^{4} Physics Department, Gothenburg University, 41296 Göteborg, Sweden
^{5} Institute of Physics, Silesian University in Opava, Bezručovo nám. 13, 746 01 Opava, Czech Republic
Received: 2 December 2015
Accepted: 10 January 2016
Polish doughnuts (PDs) are geometrically thick disks that rotate with superKeplerian velocities in their innermost parts, and whose long and narrow funnels along rotation axes collimate the emerging radiation into beams. In this paper we construct an extremal family of PDs that maximize both geometrical thickness and radiative efficiency. We then derive upper limits for these quantities and subsequently for the related ability to collimate radiation. PDs with such extreme properties may explain the observed properties of ultraluminous Xray sources without the need for the black hole masses to exceed ~10 M_{⊙}. However, we show that strong advective cooling, which is expected to be one of the dominant cooling mechanisms in accretion flows with superEddington accretion rates, tends to reduce the geometrical thickness and luminosity of PDs substantially. We also show that the beamed radiation emerging from the PD funnels corresponds to isotropic luminosities that obey L_{col} ≈ 0.1Ṁc^{2} for Ṁ ≫ Ṁ_{Edd}, and not the familiar and wellknown logarithmic relation, L ~ lnṀ.
Key words: accretion, accretion disks / stars: jets / stars: black holes / black hole physics / Xrays: bursts / stars: neutron
© ESO, 2016
1. Introduction
The research reported here was motivated by the question whether collimation of radiation in the funnels of Polish doughnuts (PDs) can explain superEddington luminosities of the ultraluminous Xray (ULX) sources, assuming that the ULXs are powered by accretion on the stellar mass compact objects.
1.1. SuperEddington accretion
The Eddington luminosity for an object with a mass M is given by the formula, (1)Here σ_{T} is the electron scattering crosssection and m_{p} is the proton mass. For objects powered by accretion, the corresponding Eddington accretion rate Ṁ_{Edd} is defined by^{1}(2)Observations provide several examples of objects that radiate at superEddington luminosities. In our Galaxy the bestknown examples are SS433 and GRS 1915+105 (see e.g. Fabrika et al. 2006; Fender & Belloni 2004, and references therein). Outside the Galaxy, superEddington luminosities are reached by ULX sources and tidal disruption events (e.g. Fabbiano 2006; van Velzen & Farrar 2014), as well as by several AGNs (see e.g. Du et al. 2015, and references therein). In the aspects that are relevant to our paper, the theory of superEddington accretion onto black holes was reviewed by Paczynski (1982, 1998), or Abramowicz (2005).
It is convenient to introduce the rescaled luminosity λ, rescaled accretion rate ṁ, maximal relative vertical thickness χ, and dimensionless advection strength ξ, Here H = H(R) is the vertical semithickness of the disk at the distance R from the black hole.
1.2. Polish doughnuts and the ULX sources
Fig. 1 Observed luminosities of a few ULX sources (horizontal dashed lines) compared with our extremal Polish doughnuts (blue lines) for different strength of advection ξ. We also indicate results of the global MHD simulations (red crosses) and various ad hoc models from recent literature (gray circles with letters). See text for a detailed explanation. 
The main results of our paper are summarized in Fig. 1, where luminosities inferred from different models are plotted against the corresponding mass accretion rates. Results related to the extremal PDs, discussed in Sect. 3, are indicated with blue lines. Different curves correspond to different advection strengths, expressed by the parameter ξ, as defined by Eq. (4). The solid blue line indicates ξ = 0.0 (no advection), the other three blue dashed lines are for ξ = 0.5,0.9, and 0.99. The collimated luminosity, L_{col} ≈ L/ (1−cosα), emerging from a narrow funnel with opening halfangle α, is an isotropic equivalent luminosity, that may be interpreted as the observed luminosity of sources in which radiation is strongly beamed (cf. Sect. 2.2). The theoretical limit λ_{col} = L_{col}/L_{Edd} = 0.625ṁ, derived in Sect. 4, is denoted with a black dashed line in Fig. 1. It approximates the luminosity of a nonadvective torus (thick blue continuous line) very accurately. Only this advectionless maximal configuration is close to points representing ULXs; accounting for advection leads to far too dim sources. The thin horizontal lines (with labels) correspond to the observed Xray luminosities of ULX sources. Two with known masses, a neutron star X2 in M 82 with mass M ≈ 1.4 M_{⊙} and luminosity L_{X} = 1.8 × 10^{40} erg/s, (Bachetti et al. 2014) and a black hole NGC 7793 with mass M< 15 M_{⊙} and L_{X} ≈ 5 × 10^{39} erg/s, (Motch et al. 2014), and two ULX sources with unknown or controversial masses, HLX1 in ESO 24349 with L_{X} = 1.2 × 10^{42} erg/s (Farrell et al. 2009; Godet et al. 2010), and NGC 5907 ULX1 with L_{X} ≈ 5 × 10^{40} erg/s (Walton et al. 2015). For the last source we assumed a mass of M = 10 M_{⊙}, for HLX1 we considered the cases of two proposed masses: 10^{4}M_{⊙} and 3 M_{⊙}. Circles on these lines show locations of the theoretical models, proposed by Godet et al. (2012; G), Kluźniak & Lasota (2015; K), and Lasota et al. (2015; L). The three crosses correspond to three models of black hole accretion flows from a recent magnetohydrodynamical (MHD) numerical simulation by Sa¸dowski & Narayan (2015). These simulations have been done assuming M = 3 × 10^{5}M_{⊙}, but in the L/L_{Edd} versus Ṁ/Ṁ_{Edd} relation, dependencies on the mass are scaled off.
Lasota et al. (2016) showed that slim accretion disks in which cooling is dominated by advection (Abramowicz et al. 1988; Sa¸dowski 2009) cannot be very geometrically thick, and that even at highly superEddington accretion rates, ṁ ≫ 1, the maximal relative disk thickness stays rather small, cotα = χ ≤ 5. Thus, the disks are indeed slim. On the other hand, models of PDs are wellknown to have arbitrarily large thickness, χ ≫ 5 (Abramowicz et al. 1978; Kozlowski et al. 1978; Paczyńsky & Wiita 1980; Jaroszynski et al. 1980). Here, we resolve this apparent contradiction by pointing out that the models of PDs constructed so far have been nonadvective. We reconsidered the PD models to include a strong, global advective cooling. We show that taking advection into account greatly reduces the PD thickness, deeming very thick tori construction impossible for realistic mass accretion rates. We proceed by evaluating the magnitude of radiation collimation in a funnel of a very thick PD (Sikora 1981), finding that it obeys the linear scaling L_{col} = 0.0625Ṁc^{2}, which agrees reasonably well with recent numerical simulations (Sa¸dowski & Narayan 2015) and observations of ULXs. However, when advection is accounted for, PDs cannot provide sufficient luminosity to explain ULXs.
Fig. 2 Left: shapes of tori with constant angular momentum. PDs with R_{in} → R_{mb} have an arbitrarily high H/R, see also Fig. 4. Shapes correspond to R_{in}/r_{S} = 2.025,2.01,2.005, and 2.0025. Black dots indicate the highest H/R. Right: Torus with R_{in}/r_{S} = 2.0025, corresponding to χ = cotα = 7.04. The dot indicates the highest H/R, crosses indicate the outer boundary of the region, where 75% and 90% of the total radiation is generated. We note that the scaling is different for vertical and horizontal axes. 
2. Polish doughnuts: assumptions and equations
Polish doughnuts are stationary and axially symmetric models of accretion structures around black holes. All properties of a PD are derived from a single assumed function ℓ(R): the specific angular momentum distribution at the PD surface. From an assumed ℓ(R) the PD shape H(R) is calculated together with the radiation flux at the surface f(R), the total luminosity L = ^{∫}fdS, and finally the accretion rate Ṁ. All these are given in terms of analytic (algebraic) formulae. No physical properties of the PD interior need to be considered. We stress that models of PDs do not assume anything specific about their interiors, not even about the equation of state, p = p(ρ,T). In particular, the pressure (gas and radiation) p, the density ρ, the temperature T, and the (nonazimuthal) velocity v do not appear in the model.
The PDs were constructed using the Einstein relativistic hydrodynamics equations in the Kerr geometry (, etc.), but they are often considered in the Newtonian model of the gravity of a nonrotating black hole introduced by Bohdan Paczyński. The Paczyński model assumes Newtonian hydrodynamics and the gravitational potential given by the Paczyńsky & Wiita (1980) formula, (5)where R, Z, φ are cylindrical coordinates, M is the black hole mass, and r_{S} is the gravitational (Schwarzschild) radius. In the Paczyński potential, the marginally stable orbit (ISCO) is located at R_{ms} = 3r_{S} and the marginally bound circular orbit at R_{mb} = 2r_{S}, exactly as in the case of the Schwarzschild (nonrotating) black hole.
The classic PD models assume that
 1.
 2.
the specific angular momentum ℓ at the photosphere is a known (assumed) function, (7)The angular momentum is Keplerian at the inner and outer radii of a PD, ℓ(R_{in}) = ℓ_{K}(R_{in}) and ℓ(R_{out}) = ℓ_{K}(R_{out}), where R_{mb}<R_{in}<R_{ms}.
Radiation is emitted from the photosphere at the local Eddington rate, that is, the local flux f^{rad} is given by (8)where κ is the mass absorption coefficient and g^{eff} is the effective gravity, given by
2.1. Polish doughnut shape H(R)
By integrating the differential equation (12)we derive an explicit analytic formula for the PD shape, H = H(R), (13)The vertical thickness H(R) is zero at the inner edge, R_{in}. It is easy to see that the thickness is also zero at the outer edge R_{out} given by the integral condition, (14)In the particular case of a constant angular momentum distribution, (15)from Eq. (14) we obtain (16)Then the shape of the constant angular momentum torus is (17)Two limiting shapes are the infinite, unbounded torus for R_{in} = 2r_{S} (ρ = 0) and the ring whose cross section is reduced to a point for R_{in} = 3r_{S} (ρ = 1). Examples of H(R) profiles for the constant angular momentum distribution, calculated assuming different R_{in}, are shown in Fig. 2 (left).
2.2. Polish doughnut luminosity L
Integrating f^{rad} over the PD photosphere, whose location is given by Eq. (13), gives the PD total luminosity L, (18)For thick tori most contribution to the integral (18) comes from the inner region, that is, the funnel, L_{fun} ≈ L, Fig. 2 (right). When radiation collimated by the funnel of opening halfangle α is observed, which yields a measured flux f^{r}, the collimated equivalent isotropic luminosity is calculated to be (19)where we have introduced the collimating factor β.
Fig. 3 Conservation of energy and angular momentum. 
2.3. Polish doughnut efficiency ϵ
Figure 3 illustrates the conservation of energy and angular momentum in a PD. Energy and angular momentum flow in at the outer edge and flow out at the inner edge. The rate at which energy is deposited in the PD interior is L_{0}, Here Ṁ is the accretion rate, (22)is the specific mechanical energy, Ω = ℓ/R^{2} is the angular velocity, and is the torque. Assuming the usual notorque inner boundary condition, , we derive from Eqs. (20), (21) (23)where ϵ(R_{in},R_{out}) is the efficiency of energy generation, (24)For very large PDs, R_{out} → ∞, we have e_{out} = 0 = Ω_{out}, and using the inner boundary condition ℓ_{in} = ℓ_{K}(R_{in}), we may write (25)The efficiency ϵ_{∞}(R_{in}) is the upper limit for the efficiency of a PD with an inner radius R_{in}. Nevertheless, we note that for R_{in} → R_{mb} = 2r_{S} the efficiency of a large PD tends to zero, (26)In thermal equilibrium the energy gain must be compensated for by the radiative and advective losses (27)
2.4. Polish doughnut advective cooling
Since in the PD formalism no interior physics is considered, the advective losses cannot be calculated directly. We parameterize the advective losses by L_{adv} = ξL_{0}, with some 0 ≤ ξ ≤ 1. With this parametrization we write an explicit formula for the accretion rate in terms of luminosity, (28)The total radiative efficiency of large PDs, that is, the upper limit of PD efficiency for given R_{in} and ξ, may be defined as (29)Since we are interested in placing constraints on the thickness and luminosity for a given ṁ of PDs, we use the highest efficiency ϵ_{rad} in the calculations below. Radiatively inefficient accretion flows, or RIAFs, have ϵ_{rad}(R_{in},ξ) ≪ 1. The RIAFtype large PDs have either R_{in} ≈ R_{mb}, or ξ ≈ 1, or both.
3. Extremal family of PDs
Models of PDs are determined by the specific angular momentum distribution. In this section we determine a particular distribution that for a given location of the inner edge R_{in} gives the greatest possible relative thickness, h = H/R, and the highest possible efficiency of a PD.
Let R_{max} denote the (radial) location of the maximal relative thickness h, for some unspecified angular momentum distribution (30)where χ = h(R_{max}) is the actual value of the greatest relative thickness, or in other words, the quantity that we wish to determine at its greatest extent. By “funnel” we understand the inner region of the disk, for which R<R_{max}, see Fig. 2 (right).
From Eqs. (12) and (30) we derive (31)We introduce a parameter y = ℓ(R_{max}) /ℓ_{K}(R_{max}) that indicates how close the angular momentum at R_{max} is to the local Keplerian value, and rewrite Eq. (30) as (32)where a^{2} = 1 + χ^{2} and Δ = r_{S}/R_{max}. We note that (33)and therefore in the limit Δ → 0, that is, (34)it is clear that the lower the angular momentum, the thicker the torus. Admissible specific angular momentum distributions are nondecreasing in radius as a consequence of the Rayleigh stability condition. Hence, tori with angular momentum ℓ = ℓ_{K}(R_{in}) = const. (at least for R<R_{max}) represent a family of extremely thick tori: the relative thickness h is maximal for them. For these tori we find an analytic formula for R_{max}, (35)and, given Eqs. (17) and (35), for the relative thickness, (36)The angular momentum distribution for R>R_{max} is irrelevant for the maximization of h – it may be constant, but it may as well not be. In particular, any distribution that is constant for R_{in} ≤ R ≤ R_{max} and has a tail for R_{max} ≤ R ≤ R_{out} chosen to satisfy Eq. (14) for R_{out} → ∞, maximizes both geometrical thickness and radiative efficiency. Such tori constitute our family of extremal PDs, with the relative thickness given by Eq. (36) and efficiency given by Eq. (29).
Fig. 4 Highest value of the relative thickness χ = (H/R)_{max} for PDs. The straight dashdotted line represents the asymptotic relation . The cusp opening halfangle ψ is also indicated. 
4. Properties of very thick PDs
In this section we derive the asymptotic expansion of the PD properties in the limit of very geometrically thick disks, χ ≫ 1. We expand the characteristic quantities up to a leading term in the dimensionless parameter ρ, as defined by Eq. (11), around ρ = 0. Immediately, (37)From Eq. (35) we find (38)To find the expression for the geometric thickness we use Eq. (17) to find, after some algebra, (39)from which (40)follows. Hence, very thick tori indeed correspond to ρ → 0, R_{in} → R_{mb}, see Fig. 4. At the inner edge of the torus a cusp is formed, through which matter may be dynamically advected into the black hole. From Eq. (17) we find that the cusp opening halfangle ψ (see Fig. 2) obeys (41)In the limit case of ρ → 0 we find ψ → π/ 4. As can be seen in Fig. 4, the dependence of the cusp opening halfangle on ρ is rather weak, and ψ = π/ 4 is a good approximation for all very thick tori.
Fig. 5 Radiative efficiency ϵ_{rad}. The continues line shows ξ = 0, the dashed lines represent advection strength ξ = 0.5,0.9, and 0.99. The higher the advection, the lower the efficiency. The dashdotted line represents the asymptotic formula for an efficiency η = ρ/ 4. 
The limit efficiency, Eq. (29), can be expanded as (42)This is shown in Fig. 5. Since , the funnel collimating factor is (43)Expanding the luminosity is much more involved because we need to integrate over a nontrivial PD surface. However, assuming a conical shape of the funnel and integrating the radiation flux over such a surface (we recall that for very thick tori L_{fun} ≈ L), we find the following result (44)which can be expanded for ρ ≈ 0 as (45)Formula (45) is not expected to yield an accurate result. We find that it estimates the luminosity with a relative error of about 50%. We proceed by assuming that Eq. (45) at least provides a proper functional form of the asymptotic relation between ρ and L in the form (46)for some constants A and B. An accurate fit can be found for A = 2 and B = 0.225, see Fig. 6. Combining the results for luminosity and efficiency, we find the mass accretion rate (47)We see now that because ϵ_{rad} ~ ρ ~ 1 /β, the collimated luminosity is (48)or, in dimensionless units, (49)Clearly, for very thick PDs, the luminosity collimated by the narrow funnel scales linearly with the mass accretion rate. From plotting the relative thickness χ against the mass accretion rate ṁ, Fig. 7, we see that while formally the PDs can be arbitrarily thick, even a huge accretion rate of 10^{5}Ṁ_{Edd} only provides χ ≈ 40. Furthermore, strong advection reduces the thickness by a factor on the order of (1−ξ)^{1/2}, and for an advection parameter ξ = 0.99 even Ṁ = 10^{5}Ṁ_{Edd} cannot produce a disk with a relative thickness larger than χ = 6.
Fig. 6 Luminosity λ and accretion rate ṁ for large, nonadvective PDs with constant angular momentum. The luminosity integrated over whole torus is indicated with a thick line, while the thin line corresponds to the funnel luminosity (i.e., integrated over R < R_{max}). The dashdotted lines correspond to the asymptotic formulae for very thick tori, Eqs. (46), (47). 
Fig. 7 PD relative thickness χ as a function of ṁ. The continuous line shows ξ = 0, the dashed lines represent ξ = 0.5,0.9, and 0.99. 
4.1. Accretion “branches”
It is customary to display analytic models of accretion disks in the ṁ versus Σ(R) parameter space, where (50)is the surface density at a given cylindrical location R. An example is shown in Fig. 8 which locates the main types of thin (or semithin) disk models: ShakuraSunyaev, slim, and ADAFs. The question of where in this figure the PDs are located, cannot be answered precisely and unambiguously. There is no unique model for the PD interiors, therefore the function Σ(R) is not known.
Fig. 8 Accretion disk equilibria at a particular radial location R = 5r_{S} and for a black hole mass 10 M_{⊙}. Flows above the dotted lines τ = 1 are optically thin. The solid Sshaped line on the right represents the optically thick disks, the solid line on the left represents the optically thin disks. The upper branches of both lines represent advectiondominated disks (slim and ADAFs). The ShakuraSunyaev viscosity parameter is assumed to be α = 0.1. We note the location of the PDs. This figure is adopted (with some simplifications) from Lasota (2015), the original appeared in Abramowicz et al. (1995). 
In the “classic” zeroadvection case considered in the 1980s , the region of the thick PDs is above the H/R = 1 line and to the right of the τ = 1 line, as for example in the review by Abramowicz & Fragile (2013). However, when strong advection is added, as we did here, the PDs move into the region occupied by slim disks.
5. Conclusions
We have demonstrated that at high accretion rates, when advection probably is a dominant source of cooling, the relative thickness of PDs is significantly reduced. Doughnuts with strong advection may be considered as approximate models of slim disks, which is consistent with the conclusion of Sadowski et al. (2016). We considered advection as an additional (to radiation) cooling process, using only very general, global, conservation laws for mass, energy, and angular momentum. In particular, we did not specify whether advection is radial or vertical. Thus it follows from our general considerations that any type of strong advection would keep the thickness of the PDs at relatively low values. Our results cast doubts on whether collimation by a thickdisk funnel is an adequate model for the ULXs. When advection is taken into account, even very high mass accretion rates cannot produce sufficient collimated luminosity, cf. Fig. 1. Moreover, while a nonadvective disk seems to agree with the numerical GRMHD results of Sądowski & Narayan (2015), the latter does not report any increase of thickness with the mass accretion rate.
To conclude, PDs are very minimalist models of accretion flows at superEddigton accretion rates. They only give the photospheric properties of these flows and do not refer to their interiors. This is both a deficiency and a virtue.
Acknowledgments
Discussions on superEddington flows with Andrew King and Olek Sa¸dowski were of great help when writing this paper. This work was supported by the Polish NCN grants 2013/09/B/ST9/00060, DEC2012/04/A/ST9/00083 and UMO2013/08/A/ST9/00795, the Czech “Synergy” (Opava) grant ASCRM100031242 CZ.1.07/2.3.00/20.0071. M.W. acknowledges support of the Foundation for Polish Science within the START Programme. J.P.L. was supported in part by the French Space Agency CNES.
References
 Abramowicz, M. A. 2005, in Growing Black Holes: Accretion in a Cosmological Context, eds. A. Merloni, S. Nayakshin, & R. A. Sunyaev, 257 [Google Scholar]
 Abramowicz, M. A., & Fragile, P. C. 2013, Liv. Rev. Relativ., 16, 1 [Google Scholar]
 Abramowicz, M., Jaroszynski, M., & Sikora, M. 1978, A&A, 63, 221 [NASA ADS] [Google Scholar]
 Abramowicz, M. A., Czerny, B., Lasota, J. P., & Szuszkiewicz, E. 1988, ApJ, 332, 646 [NASA ADS] [CrossRef] [Google Scholar]
 Abramowicz, M. A., Chen, X., Kato, S., Lasota, J.P., & Regev, O. 1995, ApJ, 438, L37 [NASA ADS] [CrossRef] [Google Scholar]
 Bachetti, M., Harrison, F. A., Walton, D. J., et al. 2014, Nature, 514, 202 [NASA ADS] [CrossRef] [Google Scholar]
 Du, P., Hu, C., Lu, K.X., et al. 2015, ApJ, 806, 22 [NASA ADS] [CrossRef] [Google Scholar]
 Fabbiano, G. 2006, ARA&A, 44, 323 [NASA ADS] [CrossRef] [Google Scholar]
 Fabrika, S., Karpov, S., Abolmasov, P., & Sholukhova, O. 2006, in Populations of High Energy Sources in Galaxies, eds. E. J. A. Meurs, & G. Fabbiano, IAU Symp., 230, 278 [Google Scholar]
 Farrell, S. A., Webb, N. A., Barret, D., Godet, O., & Rodrigues, J. M. 2009, Nature, 460, 73 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Fender, R., & Belloni, T. 2004, ARA&A, 42, 317 [NASA ADS] [CrossRef] [Google Scholar]
 Godet, O., Barret, D., Webb, N., & Farrell, S. 2010, ATel, 2821, 1 [NASA ADS] [Google Scholar]
 Godet, O., Plazolles, B., Kawaguchi, T., et al. 2012, ApJ, 752, 34 [NASA ADS] [CrossRef] [Google Scholar]
 Jaroszynski, M., Abramowicz, M. A., & Paczynski, B. 1980, Acta Astron., 30, 1 [NASA ADS] [Google Scholar]
 Kluźniak, W., & Lasota, J.P. 2015, MNRAS, 448, L43 [NASA ADS] [CrossRef] [Google Scholar]
 Kozlowski, M., Jaroszynski, M., & Abramowicz, M. A. 1978, A&A, 63, 209 [NASA ADS] [Google Scholar]
 Lasota, J.P. 2015, in Astrophysics of Black Holes – From fundamental aspects to latest developments, ed. C. Bambi (Springer: Astrophysics and Space Science Library), in press [arXiv:1505.02172] [Google Scholar]
 Lasota, J.P., King, A. R., & Dubus, G. 2015, ApJ, 801, L4 [NASA ADS] [CrossRef] [Google Scholar]
 Lasota, J.P., Vieira, R. S. S., Sadowski, A., Narayan, R., & Abramowicz, M. A. 2016, A&A, 587, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Motch, C., Pakull, M. W., Soria, R., Grisé, F., & Pietrzyński, G. 2014, Nature, 514, 198 [NASA ADS] [CrossRef] [Google Scholar]
 Paczynski, B. 1982, Mitteilungen der Astronomischen Gesellschaft Hamburg, 57, 27 [Google Scholar]
 Paczynski, B. 1998, Acta Astron., 48, 667 [NASA ADS] [Google Scholar]
 Paczyńsky, B., & Wiita, P. J. 1980, A&A, 88, 23 [NASA ADS] [Google Scholar]
 Sądowski, A. 2009, ApJS, 183, 171 [NASA ADS] [CrossRef] [Google Scholar]
 Sądowski, A., & Narayan, R. 2015, MNRAS, 453, 3213 [NASA ADS] [CrossRef] [Google Scholar]
 Sadowski, A., Lasota, J.P., Abramowicz, M. A., & Narayan, R. 2016, MNRAS, 456, 3915 [NASA ADS] [CrossRef] [Google Scholar]
 Sikora, M. 1981, MNRAS, 196, 257 [Google Scholar]
 van Velzen, S., & Farrar, G. R. 2014, ApJ, 792, 53 [NASA ADS] [CrossRef] [Google Scholar]
 Walton, D. J., Harrison, F. A., Bachetti, M., et al. 2015, ApJ, 799, 122 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
Fig. 1 Observed luminosities of a few ULX sources (horizontal dashed lines) compared with our extremal Polish doughnuts (blue lines) for different strength of advection ξ. We also indicate results of the global MHD simulations (red crosses) and various ad hoc models from recent literature (gray circles with letters). See text for a detailed explanation. 

In the text 
Fig. 2 Left: shapes of tori with constant angular momentum. PDs with R_{in} → R_{mb} have an arbitrarily high H/R, see also Fig. 4. Shapes correspond to R_{in}/r_{S} = 2.025,2.01,2.005, and 2.0025. Black dots indicate the highest H/R. Right: Torus with R_{in}/r_{S} = 2.0025, corresponding to χ = cotα = 7.04. The dot indicates the highest H/R, crosses indicate the outer boundary of the region, where 75% and 90% of the total radiation is generated. We note that the scaling is different for vertical and horizontal axes. 

In the text 
Fig. 3 Conservation of energy and angular momentum. 

In the text 
Fig. 4 Highest value of the relative thickness χ = (H/R)_{max} for PDs. The straight dashdotted line represents the asymptotic relation . The cusp opening halfangle ψ is also indicated. 

In the text 
Fig. 5 Radiative efficiency ϵ_{rad}. The continues line shows ξ = 0, the dashed lines represent advection strength ξ = 0.5,0.9, and 0.99. The higher the advection, the lower the efficiency. The dashdotted line represents the asymptotic formula for an efficiency η = ρ/ 4. 

In the text 
Fig. 6 Luminosity λ and accretion rate ṁ for large, nonadvective PDs with constant angular momentum. The luminosity integrated over whole torus is indicated with a thick line, while the thin line corresponds to the funnel luminosity (i.e., integrated over R < R_{max}). The dashdotted lines correspond to the asymptotic formulae for very thick tori, Eqs. (46), (47). 

In the text 
Fig. 7 PD relative thickness χ as a function of ṁ. The continuous line shows ξ = 0, the dashed lines represent ξ = 0.5,0.9, and 0.99. 

In the text 
Fig. 8 Accretion disk equilibria at a particular radial location R = 5r_{S} and for a black hole mass 10 M_{⊙}. Flows above the dotted lines τ = 1 are optically thin. The solid Sshaped line on the right represents the optically thick disks, the solid line on the left represents the optically thin disks. The upper branches of both lines represent advectiondominated disks (slim and ADAFs). The ShakuraSunyaev viscosity parameter is assumed to be α = 0.1. We note the location of the PDs. This figure is adopted (with some simplifications) from Lasota (2015), the original appeared in Abramowicz et al. (1995). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.