arxiv: v1 [astro-ph.ga] 2 Jun 2014

Transcription

1 Astronomy & Astrophysics manuscript no. komissarov c ESO 214 June 3, 214 Millimeter images and spectra of Sgr A* from a magnetized torus W. Yan 1, F. H. Vincent 1, M. A. Abramowicz 1,2,3, A. A. Zdziarski 1, and O. Straub 4 1 Nicolaus Copernicus Astronomical Center, ul. Bartycka 18, PL--716 Warszawa, Poland wyan@camk.edu.pl 2 Physics Department, Gothenburg University, SE Göteborg, Sweden marek.abramowicz@physics.gu.se 3 Institute of Physics, Faculty of Philosophy and Science, Silesian University in Opava, Bezručovo nám. 13, CZ-7461 Opava, Czech Republic 4 LUTH, Observatoire de Paris, CNRS, Université Paris Diderot, 5 place Jules Janssen, 9219 Meudon, France odele.straub@obspm.fr arxiv: v1 [astro-ph.ga] 2 Jun 214 Received ; accepted ABSTRACT Context. The supermassive black hole, Sagittarius Sgr) A*, in the centre of our Galaxy has the largest angular size in the sky among all astrophysical black holes. Its silhouette, assuming no rotation, spans 5 µas. Resolving such dimensions has long been out of reach for astronomical instruments until a new generation of interferometers being operational during this decade. Of particular interest is the Event Horizon Telescope EHT) with resolution 2 µas in the millimeter-wavelength range.87 mm 1.3 mm. Aims. We investigate the ability of the fully general relativistic Komissarov 26) analytical magnetized torus model to account for observable constraints at Sgr A* in the millimeter domain. The impact of the magnetic field geometry on the observables is also studied. Methods. We calculate ray-traced millimeter-wavelength images of a magnetized accretion torus surrounding the central black hole in Sgr A* for two different geometries of the magnetic field. We assume stationarity, axial symmetry, constant specific angular momentum, polytropic equation of state and small optical depth. Results. We show that the Komissarov model is capable of reproducing observable constraints in the millimeter domain, and in particular in the observable domain of the EHT. However, no signature of the magnetic field geometry can be retrieved from spectra or images. Conclusions. The ability of the Komissarov model to account for observations of Sgr A* in the millimeter domain is interesting in the perspective of the future EHT. Such an analytical model allows very fast computations. It will thus be a suitable test bed for investigating large domains of physical parameters, as well as non-black-hole compact object candidates. Key words. Galaxy: centre Accretion, accretion discs Black hole physics Relativistic processes 1. Introduction A black hole is characterised by its event horizon, a boundary that causally separates the external universe from its interior and makes the black hole appear black. The reason why there are so far no direct observations of their immediate environment is that astrophysical black holes observed from Earth have a very small apparent angular size in the sky. The measurable apparent size of a black hole refers to the size of its silhouette sometimes also referred to as black hole shadow, e.g., Falcke et al. 2), a very thin ring of light produced by photons that originate from just outside the photon orbit, i.e. the innermost circular orbit around the black hole. Photons coming from the region inside the photon orbit may either fall below the event horizon, or may still escape but with such high redshifts that the whole area inside the black hole silhouette appears dark. The black hole with largest silhouette of all is the supermassive black hole in the centre of our Galaxy, Sagittarius Sgr) A*. Given its mass M = M and distance D = 8.3 kpc Ghez et al. 28; Gillessen et al. 29a,b), the angular size of its silhouette projected on sky, taking into account the magnification due to gravitational lensing, is only 5 µas for a non-rotating black hole. This angular size is a decreasing function of the black hole spin, reaching 4 µas for a maximally rotating black hole. Sgr A* was first observed in the radio band Balick & Brown 1974), but its observed emission ranges from radio to X-ray energies. The most remarkable feature of Sgr A* is its complex variability at all observable wavelengths. The luminosity fluctuations increase with increasing wave energy, from a factor of a few at radio to a few orders of magnitude in the X-ray band see e.g. Genzel et al. 21, for a review). The spectral peak lies in the millimeter radio band and brings forth a peak luminosity of 1 36 erg s 1. The accretion structure around Sgr A* is thus extremely dim given its enormous mass. Therefore, adequate disc models describe a radiatively inefficient emitter like an advection dominated accretion flow ADAF, Narayan & Yi 1995). The term advection means here that a large part of the gravitationally liberated thermal energy is not converted into radiation but carried inward with the ionised, hot accretion flow. In the millimeter radio range, i.e. at wavelengths corresponding to the spectral peak of Sgr A*, the 1

2 Event Horizon Telescope EHT, Doeleman et al. 29), operational in , and the orbital telescope RadioAstron Sokolovsky et al. 213), launched in 211, will be able to perform high resolution Very Long Baseline Interferometry VLBI) observations. Like this, images of the close environment of the black hole will be obtained, in particular images of the accretion flow. Recently, the emitting zone around Sgr A* was constrained to merely 37 µas by millimeter VLBI Doeleman et al. 28). Since this is smaller than the silhouette size of the presumed black hole, the observed emission should originate from the surrounding accretion flow. These new observational possibilities on Sgr A* have stimulated a lot of recent research dedicated to modeling the accretion flow surrounding this black hole see e.g. Goldston et al. 25; Noble et al. 27; Chan et al. 29; Mościbrodzka et al. 29; Dexter et al. 21; Shcherbakov et al. 212; Mościbrodzka & Falcke 213), as well as works specifically related to the future EHT reviewed by Broderick et al. 214). The hope is that a detailed knowledge of theoretically predicted observational appearance of the structure of the accretion flow in Sgr A* will provide powerful and reliable tools to test Einstein s general relativity at its strong-field limit. While eventually sophisticated numerical models like Global General Relativistic Radiation Magnetohydrodynamics - GGRRMHD, McKinney et al. 213) of Sgr A* will be used to make a meaningful comparison between theory and observations, in the foreseeable future simple analytic models will be invaluable as a secure guide in the vast parameter space that needs to be explored. Following this idea, we have recently constructed an analytic optically thin Polish Doughnut model of Sgr A* Straub et al. 212). It assumed that the magnetic field in Sgr A* had no global structure, but instead was chaotic i.e. locally isotropic. In this paper we make the next logical step by considering a model with a globally ordered toroidal) magnetic field. For this, we explicitly calculate images and spectra of Sgr A* relevant for the EHT observations, i.e. in the wavelength range.87 mm 1.3 mm. We use the Komissarov 26) analytic model of a magnetized optically thin Polish Doughnut and follow all its assumptions. In the Komissarov model, all general relativistic effects, and influence of the toroidal) magnetic field are fully and exactly taken into account. They are calculated from the first principles with no approximation. The presence of a magnetic field is important in calculations of the synchrotron radiation emissivity, which we also derive following Wardziński & Zdziarski 2). We consider here a torus-shaped, barotropic, and stationary disk with axisymmetry and constant angular momentum around a Kerr black hole. The disk is fully ionized. These assumptions reflect the basic physics of the real object. For these two scenarios - chaotic and toroidal magnetic fields - we explicitly calculate images and spectra of Sgr A* relevant for EHT observations, i.e. in the wavelength range.87 mm 1.3 mm. Our main goal is to determine whether the Komissarov model can account for millimeter observational constraints for Sgr A*, as well as whether the geometry of the magnetic field chaotic or toroidal) leads to observable differences. We summarize the basic features of the magnetized torus model and its synchrotron radiation in Sections 2 and 3, respectively. Section 4 describes the ability of the Komissarov model to reproduce observable constraints and Section 5 presents conclusions and perspectives. 2. Magnetized accretion torus We want to investigate the observable appearance of an accretion torus for two distinct magnetic field configurations: toroidal and isotropic i.e. chaotic) Toroidal magnetic field the Komissarov model) We constructed a magnetized accretion torus at the Galactic centre using the model developed by Komissarov 26), which describes analytically a polytropic accretion torus with toroidal magnetic field in the Kerr spacetime. The fluid 4-velocity is assumed to be u = u t,,, u ϕ ) 1) using Boyer-Lindquist coordinates. We assume a constant specific angular momentum l u ϕ /u t. 2) This quantity is expressed in terms of the dimensionless specific angular momentum λ = l l ms l mb l ms, λ 1 3) where l ms and l mb are the specific angular momentum at the marginally stable and bound orbits respectively. These assumptions fully determine the 4-velocity. The gas and magnetic pressures are assumed to follow the polytropic prescription p = κh k ; p m = κ m L k 1 h k 4) where p is the gas pressure, p m is the magnetic pressure, κ and κ m are polytropic constants, k is the polytropic index assumed identical for gas and magnetic pressures), h = p + ρ is the particle enthalpy where ρ is the gas energy density, and L g 2 tϕ g tt g ϕϕ where g µν is the Kerr metric. The conservation of stress-energy leads to W s W + k κ + κm L k 1) h k 1 = 5) k 1 where the potential W = ln u t is used. We assume that the torus fills its Roche lobe, which fixes the central radius of the torus and its surface, thus the values of the potential at the centre, W c, and at the surface, W s, of the torus. This immediately gives h = h c ω κ + κ ml k 1 ) 1/k 1) c κ + κ m L k 1 6) where ω = W W s )/W c W s ) and L c is the value of L at the center of the torus. The polytropic constants κ and κ m can be expressed according to κ = h 1 k c W c W s ) k 1 β c, 7) k 1 + β c κ m = L1 k c κ β c 2

3 where β c is the central magnetic pressure ratio, β c p c /p m,c. The electron number density n e = h κhk µ e m u 8) is then known analytically, as well as the magnetic pressure, the expression of which is p m = h c β c β c ) ) k 1 k L L c 1/k 1) ) k 1 L W c W s ) ω k/k 1), L c k 1 k where the first line contains fluid parameters whereas the second line is dictated by spacetime geometry. We note that the magnetic pressure is directly proportional to the central enthalpy, but has a non-trivial dependency on β c. The magnetic pressure has a different expression depending on the magnetic field geometry. For a chaotic, isotropic field it is 1/3 of the magnetic energy density, whereas it is equal to the magnetic energy density for a toroidal field Hughes 1991). Thus, here 9) p m = B2, toroidal magnetic field 1) 8π where B is the magnetic field 3-vector magnitude in the fluid frame. The magnetic field in the Boyer-Lindquist frame is assumed to be toroidal, b = b t,,, b ϕ ), and to be orthogonal to the fluid 4-velocity, b u =. In the fluid frame, i.e. in an orthonormal tetrad u, e i ), its components are b =,,, B). The norm of b is thus equal to the quantity B, the magnitude of the magnetic field 3-vector in the fluid frame, which is known analytically. It is then easy to get b ϕ = b t = l b ϕ. B gϕϕ + 2l g tϕ + l 2 g, 11) tt which is fully known analytically. Let us consider one synchrotron photon emitted at a given point of the torus. Let p be the 4-vector tangent to the photon geodesic and l be its projection orthogonal to u. The angle ϑ between the magnetic field b and the direction of emission is given by l b = l b cos ϑ. It is known analytically as well. We now note that such an accretion torus cannot be made of a perfect gas. If it were, then pm u /ρk B ) = T/µ e where T is the electron temperature and k B is the Boltzmann constant. However, it is easy to see that p/ρ is independent of the central value of the enthalpy h c. Thus the temperature would be independent of h c as well, and would be purely determined by the geometry of spacetime, which does not make sense. We will still assume that there exists a relation T = Cp/ρ where C is a constant, but does not take its perfect-gas value. Rather, we choose T c at the center of the torus and define the constant C by T c = Cp c /ρ c. Then T = T c ρ ρ c ) k 1 12) depends on the choice of T c, and no longer only on spacetime geometry Isotropic magnetic field An accretion torus with isotropic i.e. chaotic) magnetic field has already been studied around Sgr A* by Straub et al. 212). This model is simply the limit of the Komissarov 26) model with κ m =. The section above thus directly applies to this simpler case. The magnetic field strength is obtained by assuming that the magnetic pressure is everywhere related to the gas pressure through p m = 1 β p 13) thus, the β parameter is valid in the whole torus, not only at its center. Then the magnetic field magnitude is known from its link to the magnetic pressure, which for a chaotic magnetic field is given by Hughes 1991) p m = B2, chaotic magnetic field. 14) 24π 3. Synchrotron radiation Millimeter-wavelength emission coming from Sgr A* has been attributed to a region close to the event horizon of the supermassive black hole and can be explained by thermal synchrotron radiation see e.g. Genzel et al. 21) Toroidal magnetic field Wardziński & Zdziarski 2) show that for a mildly relativistic Maxwellian electron distribution, n e γ) = n e γγ 2 1) 1/2 exp[ γ ], 15) θ e K 2 1/θ e ) θ e where θ e = k B T/m e c 2 ), m e being the electron mass, γ = 1 v 2 /c 2 ) 1/2 is the Lorenz factor and K 2 is a modified Bessel function, the emission coefficient is j dir ν ) = πe2 2c νν ) 1/2 X γ )n e γ ) cot2 ϑ γ 2 [ 1 1 γ 2 ) cos 2 ϑ ] 1/4 Zϑ, γ ) 16) where ν eb/2πm e c) is the cyclotron frequency. The superscript dir means that this emission coefficient depends on the angle between the magnetic field and the direction of emission, no angle averaging has been performed. Then, 3

4 [ γ = [ 1 + 2νθ e ν νθ e sin 2 ) 1 ] ϑ 2ν 4νθe ) 2 ] ν sin ϑ θ e.8 θ e.8 17) is the Lorenz factor of those thermal electrons that contribute most to the emission at ν, and [ 2θe γ 2 ] 1/2 1) γ3γ 2, θ e.8 1) X γ) = 18) ) 1/2 2θe, θ e.8 3γ t γ 2 1) 1 ν1 + t 2 ) 2 sin ϑ, n, ν γ { } 2n t exp[1 + t 2 ) 1 2 ] Zϑ, γ) = t 2 ) ) Synchrotron radiation becomes self-absorbed below a critical frequency ν c where the plasma becomes optically thick. The emitted spectrum below this frequency will follow the Rayleigh-Jeans emission law. Following Narayan & Yi 1995) we determine the self-absorption critical frequency by asking that at this frequency and at the current radius r, the synchrotron emission over the volume of the sphere with radius r equates the Rayleigh-Jeans emission from the surface of this sphere, i.e. 4 3 πr3 j dir ν ν c ) = πb ν ν c ) 4πr 2 2) where B ν = 2k B T ν 2 /c 2 is the Rayleigh-Jeans emission law. The synchrotron emitted spectrum is then smoothly connected with a Rayleigh-Jeans spectrum below ν c Isotropic magnetic field Wardziński & Zdziarski 2) give the angle-averaged limit of the synchrotron emission [ jν avg = 21/6 π 3/2 e 2 ) ] 1/3 n e ν 9v exp 21) 3 5/6 ck 2 1/θ e )v1/6 2 where v = ν/ν θ 2 e). The superscript avg now refers to an emission coefficient after angle averaging. There is a multiplicative correction factor a in this expression in Wardziński & Zdziarski 2), which is close to unity, and which we do not take into in this work as we are not interested in very precise values of synchrotron fluxes. We note that the equation above is only valid for a mildly relativistic plasma, θ e 1. This condition is satisfied in our model. The self-absorbed synchrotron is treated in the same way as for the toroidal magnetic field case. 4. Ray-traced images and spectra We use the open-source 1 ray-tracing code GYOTO Vincent et al. 211) to compute images and spectra of the magnetized accretion torus. We have studied the impact of the models parameters on the observables in order to produce simulations that are in agreement with current constraints. Restricting ourselves to the millimeter domain, we want our model to be able to reproduce the observed millimeter spectrum of Sgr A*. The corresponding images should predict a correct size of the emitting zone which was constrained to µas at 1.3 mm by Doeleman et al. 28). In order to find a correct set of parameters, we have proceeded partly by trial and error and partly using the minimization of a χ 2 criterium on a grid of parameters. However, the bestbet parameters that we obtained should not be considered as a proper fit to the data. For instance, we never varied the polytropic index. We only consider these parameters as being able to reasonably fit all observable constraints. Note that this fitting procedure was performed for the toroidal magnetic field model. The same parameters are used for the isotropic model except for the central enthalpy h c which is divided by a factor 3 in order to compensate the difference of magnetic pressure between the two models and allow to obtain the same value of magnetic field magnitude for both cases see Eqs. 9, 1 and 14). A more general and robust fit of our model to observed data goes beyond the scope of this article and will be investigated in a future work. Table 1 gives the best-bet parameters of our model. We note that they are in good agreement with independent GRMHD results see e.g. Mościbrodzka et al. 29). In particular we find the same tendency towards high spin values as many other authors using different models Mościbrodzka et al. 29; Dexter et al. 21; Shcherbakov et al. 212; Drappeau et al. 213). However, Broderick et al. 29) find a most probable spin consistent with a slowly-rotating black hole. The mean number density and electron temperature of our model reach n e = cm 3 and θ e =.25. The mean magnitude of the magnetic field is of 2 G. We note that the value of λ has a very important impact on the size of the accretion torus Straub et al. 212): the torus inflates as λ increases towards 1. The value of λ =.9 used here corresponds to a torus with radial and vertical extensions of the order of 3 GM/c 2. Figure 1 illustrates the dependency of the torus cross-section with λ and spin. Table 1. Best-bet Komissarov model parameters m u is the atomic mass unit). parameter value spin a.95 inclination i 85 angular momentum λ.9 gas/magnetic pressure ratio β c, β 1 polytropic index k 5/3 central density cm 3 ) h c/m uc 2 ) central electron temperature K) T c Figure 2 shows the synchrotron spectra associated to the toroidal and isotropic configuration of the accretion torus, for the same values of parameters given in Table 1. The 1 Freely available at the URL 4

5 6 6 4 a=.5, λ=.9 4 a=.5, λ=.9 2 a=.5, λ=.7 2 a=.95, λ=.9 best-bet) z z r r Fig. 1. Cross-section of the torus surface when varying the dimensionless angular momentum λ left) and the spin a right). The torus inflates for increasing λ and decreasing spin. The best-bet torus is highlighted in the right panel. The axes are in units of the gravitational radius GM/c 2. ν Fν [erg s 1 cm 2 ] Best-bet Komissarov model Corresponding isotropic model ν [Hz] Fig. 2. Best-bet synchrotron spectrum emitted by an accretion torus surrounding Sgr A* with a toroidal solid) or isotropic dashed) magnetic field. The observed error bars are taken from Falcke et al. 1998); Marrone et al. 28). Note that the low-frequency spectrum requires non-thermal synchrotron see e.g. Yuan et al. 23) which we do not model. The EHT will observe at two wavelengths marked with vertical dashed black lines at.87 mm and 1.3 mm). isotropic spectrum is higher than the toroidal spectrum because synchrotron radiation is stronger perpendicularly to the magnetic field direction in the toroidal case: most of the radiation is emitted orthogonally to the emitter s 4-velocity the magnetic field being toroidal), where it is not boosted by the beaming effect. Our best-bet Komissarov model allows to account for the observed spectral data in the range [1 GHz 7 GHz] which includes the observation domain of the EHT 23 and 345 GHz). The corresponding isotropic magnetic field model leads to a much worse fit, but this cannot be used as a way to constrain the magnetic field geometry because a different set of parameters for the isotropic field would lead to a fit as satisfactory as the one obtained here for the Komissarov model. Figure 3 shows the dependency of the millimeter spectrum for the Komissarov model as a function of the black hole spin and inclination parameters. It is interesting to note that, although the spectrum is weakly dependent on the inclination, it changes significantly with the spin parameter. In order to compute images that may be compared to observed data, we model the smearing effect due to scattering by interstellar electrons. Following Mościbrodzka et al. 29) we use a Gaussian profile with FWHM 1.39 λ 2 Bower et al. 26). At.87 mm this boils down to smearing the image with a Gaussian of FWHM 1 µas, thus the effect is not very significant. Figure 4 shows the.87 mm where the EHT will reach its highest resolution) images of the accretion torus for both magnetic field orientations, using the parameters listed in Table 1. Both magnetic field geometry allow to account for the observable constraint related to the size of the emitting 5

6 1 1 a= i=85 a=.7 i=6 ν Fν [erg s 1 cm 2 ] 1 12 a=.5 ν Fν [erg s 1 cm 2 ] 1 12 i= ν [Hz] ν [Hz] Fig. 3. Dependency of the spectrum of the Komissarov model with spin parameter a left) and inclination i right). All other parameters are fixed at their best-bet values Fig. 4. Scatter-broadened best-bet synchrotron images at.87 mm of an accretion torus surrounding Sgr A* with an isotropic left) or toroidal right) magnetic field. The dotted circles show the 1σ confidence contours for the size of the emitting zone obtained at 1.3 mm by Doeleman et al. 28). The linear color code is the same for both panels. The maximum value of specific intensity is of erg s 1 cm 2 str 1 Hz 1. region. We note that the size of these images is of only a few tens of µas, corresponding to less than 1 GM/c 2, whereas the total size of the torus is of the order of 3 GM/c 2 as illustrated in Figure 1. This is due to two reasons. Firstly, the relativistic beaming effect concentrates the light in the region where the emitter s 4-velocity is oriented towards the observer on the left of our images). Secondly, the electron temperature and number density decreases from the center of the torus outwards. As the synchrotron radiation is very sensitive to these two quantities, it quickly drops when going away from the torus center. The flux distribution is very similar for the two panels isotropic, toroidal) which makes it unlikely to be able to constrain the magnetic field geometry only from imaging the surroundings of Sgr A*. We note that a thinner torus e.g. with λ.3) leads to much more different images that may be distinguished with the EHT) because the difference in relativistic beaming effect between the two models is more pronounced. However, such a thin torus leads to a predicted emitting area much too small to account for the constraint of Doeleman et al. 28). It may seem strange that the two panels of Figure 4 do not show any silhouette, i.e. any obscure region in their center. This is due to the very high inclination of our best-bet model: the relativistic beaming effect is then so strong that the flux from the right side of our images is negligible with respect to the left side. We note also that, even for lower inclination, the silhouette may be erased depending on the index of the spectrum F ν ν α ). Figure 5 shows images and spectra of a Komissarov torus as seen from an inclination i = 45. Two different electron central densities are considered, the best-bet value, and twice this value. The corresponding spectra are illustrated on the top panel: they both fit reasonably well the data, but lead to different spectral index in the range of EHT observation. This behavior has a direct impact on the image. Indeed, the frequency at which the synchrotron radiation is emitted at one point of the torus, corresponding to one observed wavelength of.87 mm, varies by a factor of a few depending on the red- 6

7 shift factor. Thus, the behavior of the spectrum away from the observing frequency has a direct impact on the image. In particular, it affects the degree of visibility of the silhouette. Finally, we stress that the scatter-broadening, even restricted to 1 µas at.87 mm, has an important blurring effect on the silhouette. 5. Conclusion and perspectives We have constructed a millimeter-wavelength synchrotron radiative model for Sgr A* based on the fully general relativistic, analytical magnetized torus model of Komissarov 26), who assumes a purely toroidal magnetic field. We show that such a model is able to account for the observable spectral and interferometric constraints in the millimeter domain. We investigated the difference between the observable predictions of the Komissarov model with those of a magnetized torus with a chaotic magnetic field in order to investigate the impact of the magnetic field geometry on the spectra and images. We show that it is unlikely that future observations of these observables only will be able to constrain the magnetic field geometry in the vicinity of Sgr A*. Our finding that the Komissarov model is able to reproduce observable constraints in the millimeter domain is interesting in the perspective of the future millimeter very-long-baseline interferometer EHT. Such an analytical model allows very fast computations as opposed to GRMHD simulations. For example, one spectrum or one image at a relatively low resolution, sufficient to fit the data, requires less than one minute of computation on a standard laptop for our model. This allows us to investigate large domains of physical parameters, which is not doable with GRMHD simulations because of the computational time limitation. In this perspective, we believe that the Komissarov model for Sgr A* as developed in this work may be of interest for the future data analysis linked with the EHT project. In particular, this model may be a suitable test bed for investigating the observational counterparts of compact objects alternative to the Kerr black hole. In order to further test the possibility to constrain the magnetic field geometry in the vicinity of Sgr A*, future work will be dedicated to investigating the polarization predictions of the Komissarov and isotropic models. Future work will also be dedicated to developing a full analysis of the parameter space of the Komissarov model in order to provide a robust fit to observed data. Broderick, A. E., Fish, V. L., Doeleman, S. S., & Loeb, A. 29, ApJ, 697, 45 Broderick, A. E., Johannsen, T., Loeb, A., & Psaltis, D. 214, ApJ, 784, 7 Chan, C.-k., Liu, S., Fryer, C. L., et al. 29, ApJ, 71, 521 Dexter, J., Agol, E., Fragile, P. C., & McKinney, J. C. 21, ApJ, 717, 192 Doeleman, S., Agol, E., Backer, D., et al. 29, in Astronomy, Vol. 21, astro21: The Astronomy and Astrophysics Decadal Survey, 68 Doeleman, S. S., Weintroub, J., Rogers, A. E. E., et al. 28, Nature, 455, 78 Drappeau, S., Dibi, S., Dexter, J., Markoff, S., & Fragile, P. C. 213, MNRAS, 431, 2872 Falcke, H., Goss, W. M., Matsuo, H., et al. 1998, ApJ, 499, 731 Falcke, H., Melia, F., & Agol, E. 2, ApJ, 528, L13 Genzel, R., Eisenhauer, F., & Gillessen, S. 21, Reviews of Modern Physics, 82, 3121 Ghez, A. M., Salim, S., Weinberg, N. N., et al. 28, ApJ, 689, 144 Gillessen, S., Eisenhauer, F., Fritz, T. K., et al. 29a, ApJ, 77, L114 Gillessen, S., Eisenhauer, F., Trippe, S., et al. 29b, ApJ, 692, 175 Goldston, J. E., Quataert, E., & Igumenshchev, I. V. 25, ApJ, 621, 785 Hughes, P. A. 1991, Beams and jets in astrophysics Komissarov, S. S. 26, MNRAS, 368, 993 Marrone, D. P., Baganoff, F. K., Morris, M. R., et al. 28, ApJ, 682, 373 McKinney, J. C., Tchekhovskoy, A., Sadowski, A., & Narayan, R. 213, ArXiv e-prints Mościbrodzka, M. & Falcke, H. 213, A&A, 559, L3 Mościbrodzka, M., Gammie, C. F., Dolence, J. C., Shiokawa, H., & Leung, P. K. 29, ApJ, 76, 497 Narayan, R. & Yi, I. 1995, ApJ, 452, 71 Noble, S. C., Leung, P. K., Gammie, C. F., & Book, L. G. 27, Classical and Quantum Gravity, 24, 259 Shcherbakov, R. V., Penna, R. F., & McKinney, J. C. 212, ApJ, 755, 133 Sokolovsky et al. 213, ArXiv e-prints astro-ph/ Straub, O., Vincent, F. H., Abramowicz, M. A., Gourgoulhon, E., & Paumard, T. 212, A&A, 543, A83 Vincent, F. H., Paumard, T., Gourgoulhon, E., & Perrin, G. 211, Classical and Quantum Gravity, 28, Wardziński, G. & Zdziarski, A. A. 2, MNRAS, 314, 183 Yuan, F., Quataert, E., & Narayan, R. 23, ApJ, 598, 31 Acknowledgements FHV acknowledges interesting discussions with Monika Mościbrodzka at the first Black Hole Cam workshop in Effelsberg. This work was supported by three Polish NCN grants: 211/1/B/ST9/5439, 212/4/M/ST9/78, and 213/8/A/ST9/795, and the European Synergy grant CZ.1.7/2.3./2.71 aimed to support international collaboration at the Institute of Physics of the Silesian University in Opava. References Balick, B. & Brown, R. L. 1974, ApJ, 194, 265 Bower, G. C., Goss, W. M., Falcke, H., Backer, D. C., & Lithwick, Y. 26, ApJ, 648, L127 7

8 1 1 i=45,ne=ne,bb i=45, ne=2 ne,bb ν Fν [erg s 1 cm 2 ] ν [Hz] 5 i=45, ne=ne,bb 5 i=45, ne=ne,bb scatter-broadened i=45, ne=2 ne,bb 5 i=45, ne=2 ne,bb scatter-broadened Fig. 5. Spectra upper panel) and images lower panels) of a Komissarov torus seen with an inclination i = 45. The other parameters have their best-bet values, except the electron density which is varied n e,bb being its best-bet value). The images of the right column are the scatter-broadened equivalents of the pure images of the left column. Images have been normalized to their maximum value of intensity, which differs from panel to panel: only the relative distribution of flux is meaningful here. Depending on the spectral behavior, the silhouette the external limit of the non-emitting central zone, clearly visible in the bottom left panel) is more or less visible. Scatter-broadening in all cases has an important blurring effect on the silhouette. 8