Issue 
A&A
Volume 543, July 2012



Article Number  A83  
Number of page(s)  11  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201219209  
Published online  03 July 2012 
Modelling the black hole silhouette in Sagittarius A* with ion tori
^{1} Nicolaus Copernicus Astronomical Center, Bartycka 18, 00716 Warsaw, Poland
email: odele@camk.edu.pl
^{2} LESIA, CNRS UMR 8109, Observatoire de Paris, Université Pierre et Marie Curie, Université Paris Diderot, 92190 Meudon, France
^{3} LUTH, CNRS UMR 8102, Observatoire de Paris, Université Paris Diderot, 92190 Meudon, France
^{4} Department of Physics, Göteborg University, 41296 Göteborg, Sweden
^{5} Institute of Physics, Silesian University, 74601 Opava, Czech Republic
Received: 12 March 2012
Accepted: 16 May 2012
We calculate the “observed at infinity” image and spectrum of the accretion structure in Sgr A*, by modelling it as an optically thin, constant angular momentum ion torus in hydrodynamic equilibrium. The physics we consider includes a twotemperature plasma, a toroidal magnetic field, as well as radiative cooling by bremsstrahlung, synchrotron, and inverse Compton processes. Our relativistic model has the virtue of being fully analytic and very simple, depending only on eight tunable parameters: the black hole spin and the inclination of the spin axis to our line of sight, the torus angular momentum, the polytropic index, the magnetic to total pressure ratio, the central values of density and electron temperature, and the ratio of electron to ion temperatures. The observed image and spectrum are calculated numerically using the raytracing code GYOTO. Our results demonstrate that the ion torus model is able to account for the main features of the accretion structure surrounding Sgr A*.
Key words: black hole physics / accretion, accretion disks / Galaxy: center
© ESO, 2012
1. Introduction
The Galactic centre is one of the most interesting regions for scientific investigation because it is the closest available galactic nucleus, hence can be studied with resolutions that are impossible to achieve with other galaxies. The radio source Sagittarius A* (henceforth Sgr A*) is associated with a supermassive black hole at the centre of our Galaxy. The black hole mass, M = 4.3 ± 0.5 × 10^{6} M_{⊙}, has been determined in particular from the complete and highly eccentric orbit of the star S2, which passes Sgr A* at its pericentre in a mere 17 lighthours, i.e. 1400r_{S}, where r_{S} ≡ 2GM/c^{2} is the Schwarzschild (or gravitational) radius (Schödel et al. 2002; Bower et al. 2004; Gillessen et al. 2009). Given the Sgr A* distance of 8.3 ± 0.35 kpc, the Schwarzschild radius corresponds to an angular size of around 10 μas in the sky, making the Galactic centre black hole an ideal candidate for nearfuture microarcsecond interferometric technologies (Paumard et al. 2008; Eisenhauer et al. 2008; Eisner et al. 2010; Broderick et al. 2011).
In this context, several authors have calculated images of various theoretical models of accretion structures around Sgr A*. Numerical simulations extend from threedimensional nonrelativistic magnetohydrodynamics (MHD, Goldston et al. 2005; Huang et al. 2007) to twodimensional general relativistic magnetohydrodynamics (GRMHD, Mościbrodzka et al. 2009; Hilburn et al. 2010) and threedimensional GRMHD (Dexter et al. 2010; Shcherbakov et al. 2012; Dibi et al. 2012). The size and shape of the “black hole silhouette” cast by the black hole onto the accretion structure, is determined by the photon orbit and depends only on the black hole mass and spin. Therefore, with knowledge of the black hole mass at the Galactic centre the spin can be estimated by fitting the calculated size and shape of the shadow to the size and shape observed. The appearance (spatial extent, shape, brightness distribution, etc.) of the entire image, however, depends not only on black hole mass and spin but also on the details of the accretion structure around Sgr A*. Most of them are uncertain, e.g., the chemical abundance, the radiative processes, and the inclination. They have not yet been sufficiently examined, and a wide range of the relevant parameter space remains unexplored.
Fig. 1 Meridional cut through a λ = 0.3 (black curve) and a λ = 0.7 torus (grey curve) around a black hole (the black semicircle is the event horizon). Left panel: Schwarzschild case. Right panel: Kerr case with a = 0.9 M, the longdashed semicircle marking the ergosphere. 
The accretion structure around Sgr A* is most probably a radiatively inefficient advection dominated accretion flow (ADAF), as Narayan et al. (1995) rather convincingly demonstrated by means of spectral fitting. In the following, an approximate analytic torus model associated with the family of Polish doughnuts introduced by Abramowicz et al. (1978) and Jaroszynski et al. (1980) is constructed for a radiatively inefficient accretion flow to describe the accretion structure in Sgr A*. The model depends on a few tunable parameters that can be changed and adjusted so that they can cover the whole parameter space of the problem. The general model of the source and the raytracing tailored especially for it, can be applied to spectral calculations of other black hole candidates that are both optically thin (ion tori that can be used to model e.g. some spectral states in Xray binaries) and optically thick (Polish doughnuts, to model e.g. some luminous AGNs and ULXs).
Shapiro et al. (1976) designed a hot, twotemperature disc model to describe the strong Xray emission observed in Cygnus X1, which was too hot to be understood in terms of the standard model. They found an optically and geometrically thin solution branch in which ions and electrons are weakly coupled, having different temperatures and being in energy balance. This solution, however, is thermally unstable.
Ion tori were proposed independently by Ichimaru (1977) (though using different terminology) and Rees et al. (1982) to model a gas flow that emits little detectable radiation but is at the same time able to power radio jets, both in galactic sources and microquasars. Ion tori are geometrically thick, gas (ion) pressuresupported spheroidal structures, located in the inner regions of accretion flows. They are extremely optically thin and very radiatively inefficient. These ion tori consist of a fully ionised plasma, hence their name, in which protons and electrons are thermally decoupled (each following its own temperature distribution) and not in local energy balance. They are assumed to be lowα accretion flows threaded with magnetic fields, so that the ordinary molecular viscosity is suppressed in favour of magnetic/turbulent viscosity^{1}. Owing to their vertically extended shape, H/r ≲ 1, ion tori naturally create a pair of funnels through which magnetic flux can collimate and escape. Ion tori only exist in the subEddington accretion regime, Ṁ < Ṁ_{crit} < Ṁ_{Edd} (see e.g. Ichimaru 1977; Rees et al. 1982), where radiative cooling is inefficient enough for ions to remain sufficiently hot (near virial) during the whole inflow time, i.e., t_{diff} ≫ t_{accr}. Otherwise, at higher massaccretion rates, cooling causes the puffedup inner spheroidal region to gradually become opaque and deflate to a standard thin disc. As a consequence of t_{diff} ≫ t_{accr}, the ions carry most of the energy with them into the black hole. Instead of “the ions are not cooling down fast enough”, one may also and equivalently say that the disc is locally “advectively cooled” because energy is carried away by ions. Therefore, this branch of solutions was renamed to “advection dominated accretion flow” (ADAF, see e.g. Narayan & Yi 1994, 1995; Abramowicz et al. 1995; Esin et al. 1997).
In this work, we model Sgr A* by ion tori very much similar to these described by Rees et al. (1982). They are based on the mathematical description introduced by Paczyński and collaborators in their studies of Polish doughnuts Abramowicz et al. 1978; Jaroszynski et al. 1980. Although Polish doughnuts share some of their important characteristics that are completely the opposite of these of ion tori (they are optically thick, radiation pressure supported, and correspond to superEddington accretion rates, Ṁ ≫ Ṁ_{Edd}), both classes of tori have the same equipotential structure and very similar dynamical properties. In particular, they both have a “Roche lobe”, i.e. a critical equipotential that crosses itself along the “cusp” at r = r_{cusp} (see Fig. 1). Roche lobe overflow causes dynamical mass loss from the torus to the black hole, with no need of help from viscosity. Thus, the accretion flow at radii r ≤ r_{cusp} is regulated by the black hole strong gravity and not by viscous processes. The cusp should be regarded as the inner edge of the disc; for low accretion rates it corresponds to ISCO, whereas for higher accretion rates the cusp is closer to the black hole. Both ion tori and Polish doughnuts are dynamically unstable in terms of the Papaloizou & Pringle (1984) instability. However, as proven by Blaes (1987), this instability is suppressed by the Roche lobe overflow (see also Narayan & Popham 1993).
These properties, which are shared by ion tori and Polish doughnuts, are also genuine and robust for much more general toroidal structures around black holes. In particular, they do not depend on the angular momentum distribution inside the torus. It is customary to assume that angular momentum is constant, ℓ = ℓ_{0} = const., as this leads to remarkably simple final analytic formulae. This is an assumption adopted here. With this, we do not pretend to provide a realistic accretion scenario, but an easily accessible yet powerful analytical analogy. In followup papers, we will relax this assumption and calculate the observed properties of ion tori with ℓ ≠ const. Dynamical models for these tori have been calculated e.g. by Qian et al. (2009).
The paper is structured as follows. In Sect. 2, we construct the hydrodynamical torus and in Sects. 3 and 4 we discuss the radiative properties. Section 5 shows the resulting spectra in relation to broadband data of Sgr A* and a series of images. The presented spectra are neither chisquared fits nor best guesses, they merely illustrate the performance of the ion torus model. Section 6 presents the conclusions.
2. The geometry of fat tori
2.1. Kerr metric
In BoyerLindquist spherical coordinates (t,r,θ,φ), geometrical units c = 1 = G, and signature (−, +, +, +), the Kerr metric line element has the form (1)where M is the black hole mass, a ≡ J/M its reduced angular momentum (J being the black hole angular momentum), Ξ ≡ r^{2} + a^{2}cos^{2}θ, and Δ ≡ r^{2} − 2Mr + a^{2}.
Circular orbits in Kerr spacetime are obtained from the timelike geodesics of the metric in Eq. (1) for which the fourvelocity takes the form u^{μ} = (u^{t},0,0,u^{φ}). The two constants of geodesic motion, which are associated respectively with stationarity and axisymmetry, are then expressible as where (4)is the angular velocity with respect to a distant observer. The specific angular momentum is defined by (5)At circular geodesic orbits, the angular momentum has its “Keplerian” form (Bardeen et al. 1972): (6)Here and in the remainder of the article, we consider M as a fixed parameter, so that we can assume ℓ_{K} to be a function only of (r,a).
Important circular orbits in the Kerr metric are the marginally stable and marginally bound ones, corresponding to values of r (Bardeen et al. 1972) where z_{1} ≡ 1 + (1 − a^{2}/M^{2})^{1/3}^{ [}(1 + a/M)^{1/3} + (1 − a/M)^{1/3}^{] } and . On these special orbits, the specific angular momentum takes the values (9)
2.2. Fluid torus of constant specific angular momentum
We consider a fluid torus of negligible selfgravitation around a Kerr black hole. Modelling the fluid as a perfect one, the stressenergy tensor is (10)where P is the fluid pressure and ϵ the fluid proper energy density. For a purely circular motion (u^{μ} = u^{t}(1,0,0,Ω)), it can be shown that the energymomentum conservation equation, , takes the form (11)with ℓ related to the fluid fourvelocity components by Eq. (5). Assuming that the fluid is barotropic, P = P(ϵ), we introduce the enthalpy function (12)so that the lefthand side of Eq. (11) becomes simply ∇_{μ}H. Following Abramowicz et al. (1978), we consider models for which ℓ is constant within the entire torus (see also Sect. 1): (13)Accordingly, the equation of motion in Eq. (11) reduces to ∇_{μ}H = −∇_{μ}ln( − u_{t}) yielding (14)where we have introduced the potential (15)From the normalisation relation u_{μ}u^{μ} = −1, we get (16)where Ω should be considered as the function of (r,θ) (17)This last relation is easily derived by combining Eqs. (2), (3) and (5).
Abramowicz et al. (1978) proved that the cusp location should be between the mariginaly stable and the mariginaly bound orbit, which implies that ℓ_{0} must obey (18)with ℓ_{ms}(a) and ℓ_{mb}(a) given by Eq. (9). Recast in terms of the dimensionless parameter (19)the condition in Eq. (18) is equivalent to (20)Given the radial dependence in Eq. (6) of the Keplerian specific angular momentum, which takes its minimum at r = r_{ms}(a), and the constraint of Eq. (18), there are two values of r, r_{in}, and r_{c} say, for which ℓ_{K}(r,a) = ℓ_{0}. One has necessarily (21)Since for r = r_{in} or r = r_{c}, the actual angular momentum ℓ_{0} is equal to the Keplerian one, the gravity and centrifugal forces balance each other, implying that ∇_{μ}W = 0. From the first integral in Eq. (14), we then have ∇_{μ}H = 0, or equivalently, via Eq. (12), ∇_{μ}P = 0. A point where the gradient vanishes can be called either a singular point (selfcrossing of an equipotential) or an extremum. In the present case, r = r_{in} is the location where one of the equipotentials of W, the socalled Roche lobe, selfcrosses, as shown in Fig. 1. This gives rise to a cusp at the torus surface (Abramowicz et al. 1978). On the other side, r = r_{c} corresponds to the maximum of P. The circle r = r_{c} in the equatorial plane is called the centre of the torus.
2.3. Torus solution for a polytropic equation of state
To go further, we assume a polytropic equation of state (22)where K and n are two constants, n being the polytropic index and K the polytropic constant. The total energy density ϵ = ρ + Π is the sum of energy density ρ and internal energy Π, which in the nonrelativistic limit, Π ≪ ρ, reduces to ϵ ≃ ρ. Equation (12) is then readily integrated, yielding (23)The surface of the torus is defined by P = 0. From Eq. (22), this corresponds to ϵ = 0, and from Eq. (23) to H = 0. Therefore, we may rewrite the first integral of motion in Eq. (14) as (24)where the constant W_{s} is the value of the potential W at the torus surface. Denoting by H_{c} and W_{c} the values of H and W at the torus centre, Eq. (24) implies that (25)We introduce the dimensionless potential (26)From Eqs. (24) and (25), we have (27)Substituting Eq. (23) for H, we get
where ϵ_{c} is the energy density at the torus centre. Solving for ϵ, we obtain (28)At this stage, our torus model is determined by five parameters (in addition to the black hole mass M): the Kerr spin parameter a, the dimensionless specific angular momentum λ, the polytropic index n, the polytropic constant K, and the central density ϵ_{c}. From λ and a, we evaluate ℓ_{0} via Eq. (19). The values of ℓ_{0} and a fully determine the potential W(r,θ) according to Eqs. (16), (17). Since we are seeking a Rochelobe filling torus, the value W_{s} of the potential W at the torus surface must be set to the value at the Roche lobe (cf. Sect. 2.2) (29)Given W_{s} and H_{c} (deduced from ϵ_{c} in Eq. (23)), we determine W_{c} using Eq. (25). We know then entirely the dimensionless potential ω(r,θ) as given by Eq. (26). We can thus compute the energy density everywhere in the torus using Eq. (28).
We note that, by construction (cf. Eq. (26)), ω is zero at the surface of the torus and one at the centre. Between 0 ≤ ω(r,θ) ≤ 1, the toroidal equipotentials extend from the largest possible area, the torus surface, down to a single line, a circle at r_{c}. Open equipotential surfaces have ω(r,θ) < 0. The size of the Roche lobe tori depends on the spin. The higher the spin, the narrower the angular momentum distribution that sets the location of the critical radii, hence the smaller the tori. For a fixed value of a, the Roche torus is the largest, having an infinite outer radius r_{out}, when λ = 1, i.e., ℓ_{0} = ℓ_{mb}(a).
3. Adding physics
In what follows, we set n = 3/2, which corresponds to the adiabatic index γ = 1 + 1/n = 5/3 of a nonrelativistic ideal gas with no radiation pressure. This is consistent with the existence of a very optically thin medium.
3.1. Thermodynamic quantities
In optically thin gaspressuresupported ion tori, the radiation pressure can be neglected. The total pressure P is expressed as the sum of the magnetic and gas contributions, P_{mag} and P_{gas}(30)The magnetic and gas pressures are assumed to be some fixed parts of the total pressure (an assumption that is often made in analytic models of accretion structures) (31)were B is the intensity of magnetic field. The gas is assumed to be a twotemperature plasma, with μ_{i}, μ_{e}, T_{i}, and T_{e} being the mean molecular weights and temperatures of ions and electrons, respectively. The gas pressure is then expressed as the sum of the ion and electron contributions (32)where k_{B} is the Boltzmann constant and m_{u} is the atomic mass unit. We write the ion and electron temperatures as (33)where f and g are linear functions of the equipotential function ω.
We note that the model does not consider heating and cooling of the accretion structure, but assumes thermal equilibrium with a given temperature distribution. While this may not represent the real situation, it is still an adequate approximation supported by MHD simulations that selfconsistently include cooling and heating in the dynamics. Dibi et al. (2012) found that these processes are not dynamically important in the parameter range corresponding to Sgr A*. We assume here that the ion and electron temperatures are equal at the surface: T_{i} = T_{e} when ω = 0, and that they are at some particular ratio ξ at the centre ξ T_{i} = T_{e} when ω = 1. This leads to and since Eq. (32) is fulfilled when f(ω) + g(ω) = 1, one easily finds values of g(0) and g(1) that allow us to write explicit expressions for the two temperatures It is often practical to express equations in terms of the dimensionless temperatures (38)were m_{e} and m_{i} are the electron and ion masses.
Fig. 2 Ion (dashed/dotted) and electron (solid/dotdashed) temperature distribution throughout a λ = 0.3 torus (black) and a λ = 0.7 torus (grey) around a Schwarzschild black hole. We assume a central ion to electron temperature ratio of ξ = 0.1. The thin black line marks T_{vir}. 
Figure 2 shows the temperature distribution in tori. It is different from the ADAF selfsimilar solution (e.g., Narayan & Yi 1995), which being selfsimilar must be a monotonic function of radius. Here, the nonmonotonicity of the temperature curves has its roots in the existence of boundaries, the outer occurring at the outer torus radius and the inner at the cusp. The ion and electron temperatures have a maximal separation (given by ξ) at the torus centre and converge towards either end of the torus. The temperature in our model is obviously zero at these boundaries. We note that in the torus interior our curves and those of Narayan & Yi (1995) are clearly consistent. Using Eq. (32), one can now derive the relation for the polytropic constant at the torus centre (39)where T_{e,c} is the central electron temperature. The total pressure is (40)and the magnetic field (41)The thermodynamic relations between T, P, and ϵ require the knowledge of the mean molecular weight μ, which is linked to the ion and electron molecular weights by (42)with (43)where X and Y are the hydrogen and helium abundances, which are assumed to be equal to X = 0.75 and Y = 0.25 as in Narayan & Yi (1995). On the one hand, this implies that the effective molecular weight for ions is μ_{i} = 1.23 and for electrons μ_{e} = 1.14, and on the other hand it defines the electron and ion number densities, n_{e} and n_{i}, with respect to the energy density. In the following, we provide general formulae in which the factor n_{i} is replaced by the sum over ion species, (44)where Z_{j} and n_{j} are respectively the charge and number density of the jth ion species.
3.2. Radiative processes
We have to complete the set of equations given in Sect. 3 by specifying the physical processes that are involved in the radiative cooling. We consider a twotemperature plasma cooled by synchrotron radiation, inverse Compton process, and bremsstrahlung emission. A very convenient and general description of these cooling processes was presented by Narayan & Yi (1995). We closely follow their approximations and procedures, and also use their formulae (correcting one typo). For the reader’s convenience and completeness, we summarize these author’s description here. Electronpositron pair creation and annihilation is neglected but as shown by Björnsson et al. (1996) and Kusunose & Mineshige (1996) this is justified in most cases of interest.
In the following subsections, we determine the emission and absorption coefficients inside the torus corresponding to bremsstrahlung, synchrotron radiation, and Compton processes. The emission coefficient respectively (resp.) absorption coefficient, j_{ν} (resp. α_{ν}), allows us to compute the increment (resp. decrement) of specific intensity I_{ν} when progressing through the emitting medium by a small distance dl(45)The cgs unit of j_{ν} is erg cm^{3} s^{1} ster^{1} Hz^{1}, while α_{ν} is expressed in cm^{1}.
3.2.1. Bremsstrahlung
The rate at which energy is lost due to bremsstrahlung, , includes emission from both ionelectron and electronelectron collisions (Svensson 1982; Stepney & Guilbert 1983) (46)The ionelectron bremsstrahlung cooling rate is given by (47)where α_{f} = 1/137 is the fine structure constant and the function F_{ei}(θ_{e}) has the approximate form in units of [erg s^{1}cm^{3}] (48)In the original formula quoted by Stepney & Guilbert (1983), there is a number 0.42 instead of 0.48 (see Narayan & Yi 1995).
For the electronelectron bremsstrahlung cooling rate, Svensson (1982) gives the following formula (49)where the function F_{ee}(θ_{e}) is given in units of [erg s^{1} cm^{3}] and has the approximate form (50)where r_{e} = e^{2}/m_{e}c^{2} is the classical radius of electron and the Euler number η = exp( − γ_{E}) = 0.5616 is given by the EulerMascheroni constant, γ_{E}. And again, we replace in the original formula 5/4 with 1.28 (see Narayan & Yi 1995).
With the bremsstrahlung cooling rate given above, one can express the bremsstrahlung emission coefficient as (51)where h is the Planck constant, the 1/4π factor assumes isotropic emission in the emitter’s frame and is the velocityaveraged Gaunt factor given by Rybicki & Lightman (1986)(52)
3.2.2. Synchrotron cooling
The emission coefficient for synchrotron emission by a relativistic Maxwellian distribution of electrons is given by (see Pacholczyk 1970) (53)with a factor 1/4π again for isotropic emission in the emitter’s frame and the fitting function, (54)where (55)and the parameters α, β, and γ are tabulated for a range of temperatures in Mahadevan et al. (1996). The fitting formula is valid only for θ_{e} ≳ 1, i.e. T_{e} ≳ 6 × 10^{9} K, which is satisfied for applications to ion tori because the synchrotron emission that dominates the central regions of a torus comes, as in other advective flows, from relativistic electrons in the tail of the Maxwellian distribution. At lower temperatures or in the outer torus regions, the emission is dominated by bremsstrahlung. Below a critical frequency, ν_{c}, the synchrotron spectrum becomes selfabsorbed. As the flow in ion tori is very much akin to a spherical flow, this frequency can be obtained by equating the synchrotron emission in a sphere of some radius, R [cm], to the RayleighJeans blackbody emission from the surface of that sphere and solving for x_{M}, which is then substituted into (56)We assume that at low frequencies, ν < ν_{c}, the absorption is locally given by Kirchhoff’s law, i.e., the lowfrequency part of the synchrotron emission behaves like a blackbody.
3.2.3. Compton cooling
The soft bremsstrahlung and synchrotron photons in an ion torus filled with a thermal distribution of electrons are (inverse) Compton scattered to higher energies. In the central regions of the flow in particular, this can be an important cooling mechanism. There is a probability that a seed photon of some initial energy, E_{in} = hν, is in optically thin material scattered to an amplified energy , where (57)Thus, the energy exchange between electrons and photons depends only on the electron temperature θ_{e} and the probability that a photon will interact with an electron, which is given by the electron scattering optical depth τ_{es} = ^{∫}n_{e} σ_{T} dl. Dermer et al. (1991) and Esin et al. (1996) derived an approximate prescription for the energy enhancement factor due to Compton scattering, which is defined as the average energy change of a seed photon (58)where (59)The dimensionless energy is given by . As the emerging photons cannot gain more energy than the electrons they collide with have, there is an upper limit to the energy gained, x ≲ 3 θ_{e}. We note that in our simple ion torus scenario we apply, without loss of generality, the Thomson crosssection rather than the KleinNishina crosssection.
Comptonised emission is η − 1 times the seed photon distribution. The part of the spectrum that can be Comptonised lies between the critical synchrotron selfabsorption edge, x = x_{c} = hν_{c}/m_{e}c^{2}, and x = 3 θ_{e}.
Comptonisation of bremsstrahlung emission is then given by (60)We note that the corresponding expression in Narayan & Yi (1995) has a typo that we have corrected here.
Comptonisation of synchrotron radiation that is emitted mostly at the selfabsorption frequency, ν_{c}, is given by (61)We calculate Comptonisation to second order and assume that photons that are upscattered to θ_{e} form a Wien tail.
Fig. 3 Radiative cooling rates of a λ = 0.7 torus around a spin a = 0.9M black hole. The contributions from bremsstrahlung (black, dashed), synchrotron (black, solid), Comptonised bremsstrahlung (grey, dotdashed), and Comptonised synchrotron emission (grey, dotted) are shown as a function of radius. 
3.2.4. Total cooling
The total emission coefficient is the sum of all radiative contributions (62)For a medium in local thermodynamic equilibrium (LTE) at temperature T, the emission coefficient given in Eq. (62) and the absorption coefficient are related by means of Kirchhoff’s law (63)where B_{ν} is Planck’s law of blackbody radiation. For the typical optical depth of ion tori, α_{ν} is negligible, thus the absorption can be safely ignored.
Integration of j_{ν} over the whole frequency range gives the total cooling rate q^{ − } at each point (r,θ) in the torus. Figure 3 shows the different contributions to the total cooling rate from the equatorial plane, θ = π/2, of a λ = 0.3 torus.
Radiation that originates close to a black hole is influenced by various relativistic effects, such as the gravitational bending of light rays, gravitational redshift, and Doppler beaming. We raytrace the emission derived in Eq. (62) from each point (r, θ) inside the Roche lobe equipotential back to the observer explicitly including all relativistic effects in light propagation. This is described in the next two sections.
Parameters for our reference ion torus.
4. Raytracing
We use the General relativitY Orbit Tracer of the Observatory of Paris^{2} (GYOTO), presented in Vincent et al. (2011), which permits the integration of null and timelike geodesics in any analytical or numerical metric. It is optimised to compute geodesics in the Kerr spacetime efficiently, and is able to integrate the radiative transfer equation in optically thin media.
The required inputs for GYOTO are the values of the emission coefficient, j_{ν}, and absorption coefficient, α_{ν}, at any point within the torus. These quantities, defined in Sect. 3.2, are known analytically for the ion torus, so the integration is straightforward. In particular, the absorption coefficient is zero everywhere.
Photons are launched from the observer’s screen, which is assumed to be spatially at rest at r = 100 M, in a given specified solid angle. The null geodesics are integrated until the torus is reached (or until the photon escapes too far from the torus) and the specific intensity I_{ν} is updated at each step inside the torus by using Eq. (45). A map of specific intensity is thus obtained for a given set of directions on sky. It is then straightforward to compute the observed specific flux F_{ν} by summing these values of specific intensities over all directions of incidence.
Two kinds of GYOTOcomputed quantities are presented below. The spectra (Figs. 4–9) show the quantity ν F_{ν} for different values of frequencies. The images (Fig. 10) are maps of specific intensity I_{ν}.
We stress that the whole C++ code for the ion torus is now included in the open source GYOTO code, which is available at the abovementioned URL.
5. The spectra and images of ion tori
This section presents the spectra and images of an ion torus surrounding a Kerr black hole. Table 1 gives the reference values of the parameters used for the computations. The value chosen for the central density ϵ_{c} is standard for Sgr A* (see e.g. the values in Yuan et al. 2003; Liu et al. 2004; Mościbrodzka et al. 2009; DoddsEden et al. 2010). Various parameters are varied in the computations: the spin, a, the dimensionless angular momentum, λ, the magnetic to total pressure ratio, β, the electron to ion temperature ratio, ξ, the ratio of the central electron temperature to the virial temperature, and the inclination, i.
Sections 5.1 and 5.2 describe the impact of these parameters on the spectrum and image of the ion torus.
Fig. 4 Impact of spin on an ion torus spectrum: a = 0 (solid blue), 0.5M (dashed red), or 0.9M (dotted green). All other parameters are set to their reference values listed in Table 1. The black observed data are taken from Zhao et al. (2001), Zylka et al. (1995), Marrone et al. (2008) for radio and submm data, Telesco et al. (1996), Cotera et al. (1999), Eckart et al. (2006), Schödel et al. (2007) for far and midinfrared data, Genzel et al. (2003) for nearinfrared data, and Baganoff et al. (2003) for the Xray bow tie. The downward pointing arrows refer to upper values. Note that we do not present the fitted spectra here. 
Fig. 5 Impact of central temperature: T_{0}/T_{vir} = 0.01 (solid blue), 0.015 (dashed red), or 0.02 (dotted green). 
Fig. 6 Impact of angular momentum: λ = 0.3 (solid blue), 0.45 (dashed red), or 0.6 (dotted green). 
Fig. 7 Impact of pressure ratio: β = 0.01 (solid blue), 0.1 (dashed red), or 0.2 (dotted green). 
5.1. Torus spectra
Figures 4 to 9 show the impact of various parameters on the observed spectrum of an ion torus. All parameters that have not explicitly specified values are set according to Table 1. These figures show that the ion torus model is able to account for the general features of the observed data, by appropriately tuning the various parameters. This includes the Xray flare “bow tie” which in the ion torus model may originate from soft photons that are inverse Compton scattered by the same population of hot electrons that is responsible for the synchrotron emission. Only the flattening of the spectrum at low frequencies is never matched, which is due to the absence of a nonthermal electron distribution in our model, as first pointed out by Mahadevan (1998) and demonstrated by Özel et al. (2000) and Yuan et al. (2003). A more evolved model taking into account this effect will thus be developed in future work.
Figure 4 shows that varying the spin significantly alters the emission. This is because at high black hole spins the accretion torus shrinks and moves very close to the black hole. Therefore, the flux is shifted to slightly higher energies. Increasing the spin enhances the synchrotron flux, but softens bremsstrahlung.
Figure 5 shows that increasing the central temperature implies that there is a displacement of the whole spectrum towards higher fluxes. The same effect is obtained by increasing the dimensionless angular momentum λ (Fig. 6). This is because the torus puffs up like a balloon as λ increases (see Fig. 1), thus providing more flux.
Figures 7 and 8 show that increasing either the pressure ratio β or temperature ratio ξ has an opposite impact: a higher β leads to higher synchrotron flux, while a higher ξ gives a smaller flux. This is because increasing β increases the magnetic field strength, while increasing ξ decreases B (see Eq. (41)).
Figure 9 shows that the inclination has little effect on the spectrum. Although the number of illuminated pixels is smaller at higher inclination, the beaming effect is more important. These two phenomena have opposite effects, and the resulting spectrum does not change much.
Future work will be devoted to constraining the various parameters, in particular Sgr A*’s spin, by fitting these spectra to observed data.
Fig. 8 Impact of temperature ratio: ξ = 0.08 (solid blue), 0.1 (dashed red), or 0.2 (dotted green). 
5.2. Torus images
Each image is a superimposition of a firstorder image (the thick distorted ring) and higherorder images (the circles centred on the black hole) produced by photons that swirl around the black hole before reaching the observer. There is also a very fine circle of light in each image, that consists of photons originating from a location just outside the photon orbit (see e.g. Bardeen et al. 1972, for a definition). Photons that escape from the region inside the photon orbit are so severely redshifted that they create an area of reduced intensity on the observer’s screen: this is the socalled “black hole silhouette” (which is also known as the black hole shadow; Falcke et al. 2000).
Fig. 9 Impact of inclination: i = 40° (solid blue), 60° (dashed red), or 80° (dotted green). 
Fig. 10 Upper left: image of the reference ion torus corresponding to the parameters listed in Table 1, as observed by an observer on Earth, in inverted colours (the darker, the more luminous). The display is intentionally saturated for a better rendering of the less luminous parts of the image. Upper right: same as in the upper left but with the spin parameter increased to a = 0.9M. Lower left: same as in the upper left but with dimensionless angular momentum increased to λ = 0.6. Lower right: same as in the upper left but with the inclination angle decreased to i = 40°. 
Our interest in these images is the following. Knowing the particular extent and shape of the accreting region from observations allows us to constrain various parameters of the flow. In particular, the location and shape of the photon orbit contains information on the geometry of spacetime. Increasing spin displaces the thin photon ring in the image plane off centre and a ≳ 0.9M introduces in addition distortion. Black hole spin, torus dimensionless angular momentum and inclination have each a huge and characteristic impact on the observed image of the torus as depicted in Fig. 10. This is also revealed in the spectra Figs. 4 and 6. The inclination angle, however, does not influence the spectra in a commensurate way (Fig. 9).
The required technique of measurement of the black hole silhouette is already available thanks to recent progress in millimetre Very Long Baseline Interferometry (mmVLBI) (see Doeleman et al. 2008). With the nearfuture submm VLBI instruments direct observation of processes in the range between 5–40 Schwarzschild radii will become feasible and with it also the potential to image radiation from an accretion structure and, possibly, to deduce the existence of an event horizon.
6. Conclusions and perspectives
We have calculated electromagnetic spectra and images of accretion structures around the central black hole in Sgr A* using a simple analytic model of “ion tori” very similar to these described in the wellknown paper by Rees et al. (1982). Our results depend on observationally unknown tunable parameters of the model, particularly the black hole spin. The hope is that fitting these analytic models to observations could practically restrict the allowed parameter range, so that future sophisticated MHD simulations of Sgr A* (which are similar to the more advanced ones of Dibi et al. 2012) will be more tightly constrained.
One interesting question here is how accurate the simple analytic models are compared with advanced MHD simulations. We are working on a comparison of our spectra and images with these calculated by Dibi et al. (2012). If the preliminary impression that these images and spectra are indeed very similar were to be confirmed by the test, then it would be much easier to fit theory to observations, and in particular to eventually measure the black hole spin in Sgr A*.
We plan to expand our models of ion tori to include a more complete coverage of the parameter space in connection with a “general” angular momentum distribution (as in Qian et al. 2009) and to include, directly and analytically, largescale magnetic fields (as in Komissarov 2006). We also plan to calculate models of ion tori in a nonKerr background (as advocated by Psaltis and collaborators, e.g., Psaltis et al. 2008; Johannsen & Psaltis 2010) to possibly constrain gravity in the strongfield regime of alternative gravity theories to Einstein’s general relativity. In the context of flares it may also be interesting to examine nonstationary ion tori (as in Abramowicz et al. 1983).
Freely available at the URL http://gyoto.obspm.fr
Acknowledgments
We gratefully acknowledge support from the Polish grants NN203 383136 (O.S.) and NCN 2011/01/B/ST9/05439 (M.A.A.), the stipends from the Swiss National Science Foundation (O.S.) and the Région IledeFrance (F.H.V.), as well as the Swedish VR and the Czech CZ.1.07/2.3.00/20.0071 “Synergy” grants (M.A.A.) and support from the FrenchPolish LEA AstroPF programme.
References
 Abramowicz, M., Jaroszynski, M., & Sikora, M. 1978, A&A, 63, 221 [NASA ADS] [Google Scholar]
 Abramowicz, M. A., Henderson, P. F., & Ghosh, P. 1983, MNRAS, 203, 323 [NASA ADS] [Google Scholar]
 Abramowicz, M. A., Chen, X., Kato, S., Lasota, J.P., & Regev, O. 1995, ApJ, 438, L37 [NASA ADS] [CrossRef] [Google Scholar]
 Baganoff, F. K., Maeda, Y., Morris, M., et al. 2003, ApJ, 591, 891 [NASA ADS] [CrossRef] [Google Scholar]
 Bardeen, J. M., Press, W. H., & Teukolsky, S. A. 1972, ApJ, 178, 347 [NASA ADS] [CrossRef] [Google Scholar]
 Björnsson, G., Abramowicz, M. A., Chen, X., & Lasota, J. 1996, ApJ, 467, 99 [NASA ADS] [CrossRef] [Google Scholar]
 Blaes, O. M. 1987, MNRAS, 227, 975 [NASA ADS] [Google Scholar]
 Bower, G. C., Falcke, H., Herrnstein, R. M., et al. 2004, Science, 304, 704 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Broderick, A. E., Loeb, A., & Reid, M. J. 2011, ApJ, 735, 57 [NASA ADS] [CrossRef] [Google Scholar]
 Cotera, A., Morris, M., Ghez, A. M., et al. 1999, in The Central Parsecs of the Galaxy, ed. H. Falcke, A. Cotera, W. J. Duschl, F. Melia, & M. J. Rieke, ASP Conf. Ser., 186, 240 [Google Scholar]
 Dermer, C. D., Liang, E. P., & Canfield, E. 1991, ApJ, 369, 410 [NASA ADS] [CrossRef] [Google Scholar]
 Dexter, J., Agol, E., Fragile, P. C., & McKinney, J. C. 2010, ApJ, 717, 1092 [NASA ADS] [CrossRef] [Google Scholar]
 Dibi, S., Drappeau, S., Fragile, P. C., Markoff, S., & Dexter, J. 2012 [arXiv:1206.3976] [Google Scholar]
 DoddsEden, K., Sharma, P., Quataert, E., et al. 2010, ApJ, 725, 450 [NASA ADS] [CrossRef] [Google Scholar]
 Doeleman, S. S., Weintroub, J., Rogers, A. E. E., et al. 2008, Nature, 455, 78 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Eckart, A., Schödel, R., Meyer, L., et al. 2006, A&A, 455, 1 [Google Scholar]
 Eisenhauer, F., Perrin, G., Brandner, W., et al. 2008, SPIE Conf. Ser., 7013, 70132A [Google Scholar]
 Eisner, J. A., Akeson, R., Colavita, M., et al. 2010, SPIE Conf. Ser., 7734, 773411 [CrossRef] [Google Scholar]
 Esin, A. A., Narayan, R., Ostriker, E., & Yi, I. 1996, ApJ, 465, 312 [NASA ADS] [CrossRef] [Google Scholar]
 Esin, A. A., McClintock, J. E., & Narayan, R. 1997, ApJ, 489, 865 [NASA ADS] [CrossRef] [Google Scholar]
 Falcke, H., Melia, F., & Agol, E. 2000, ApJ, 528, L13 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Genzel, R., Schödel, R., Ott, T., et al. 2003, Nature, 425, 934 [CrossRef] [PubMed] [Google Scholar]
 Gillessen, S., Eisenhauer, F., Trippe, S., et al. 2009, ApJ, 692, 1075 [NASA ADS] [CrossRef] [Google Scholar]
 Goldston, J. E., Quataert, E., & Igumenshchev, I. V. 2005, ApJ, 621, 785 [NASA ADS] [CrossRef] [Google Scholar]
 Hilburn, G., Liang, E., Liu, S., & Li, H. 2010, MNRAS, 401, 1620 [NASA ADS] [CrossRef] [Google Scholar]
 Huang, L., Cai, M., Shen, Z., & Yuan, F. 2007, MNRAS, 379, 833 [NASA ADS] [CrossRef] [Google Scholar]
 Ichimaru, S. 1977, ApJ, 214, 840 [NASA ADS] [CrossRef] [Google Scholar]
 Jaroszynski, M., Abramowicz, M. A., & Paczynski, B. 1980, Acta Astron., 30, 1 [NASA ADS] [Google Scholar]
 Johannsen, T., & Psaltis, D. 2010, ApJ, 718, 446 [NASA ADS] [CrossRef] [Google Scholar]
 Komissarov, S. S. 2006, MNRAS, 368, 993 [NASA ADS] [CrossRef] [Google Scholar]
 Kusunose, M., & Mineshige, S. 1996, ApJ, 468, 330 [NASA ADS] [CrossRef] [Google Scholar]
 Liu, S., Petrosian, V., & Melia, F. 2004, ApJ, 611, L101 [NASA ADS] [CrossRef] [Google Scholar]
 Mahadevan, R. 1998, Nature, 394, 651 [NASA ADS] [CrossRef] [Google Scholar]
 Mahadevan, R., Narayan, R., & Yi, I. 1996, ApJ, 465, 327 [NASA ADS] [CrossRef] [Google Scholar]
 Marrone, D. P., Baganoff, F. K., Morris, M. R., et al. 2008, ApJ, 682, 373 [NASA ADS] [CrossRef] [Google Scholar]
 Mościbrodzka, M., Gammie, C. F., Dolence, J. C., Shiokawa, H., & Leung, P. K. 2009, ApJ, 706, 497 [NASA ADS] [CrossRef] [Google Scholar]
 Narayan, R., & Popham, R. 1993, Nature, 362, 820 [NASA ADS] [CrossRef] [Google Scholar]
 Narayan, R., & Yi, I. 1994, ApJ, 428, L13 [NASA ADS] [CrossRef] [Google Scholar]
 Narayan, R., & Yi, I. 1995, ApJ, 452, 710 [NASA ADS] [CrossRef] [Google Scholar]
 Narayan, R., Yi, I., & Mahadevan, R. 1995, Nature, 374, 623 [NASA ADS] [CrossRef] [Google Scholar]
 Özel, F., Psaltis, D., & Narayan, R. 2000, ApJ, 541, 234 [NASA ADS] [CrossRef] [Google Scholar]
 Pacholczyk, A. G. 1970, Radio astrophysics. Nonthermal processes in galactic and extragalactic sources (San Francisco: Freeman) [Google Scholar]
 Papaloizou, J. C. B., & Pringle, J. E. 1984, MNRAS, 208, 721 [NASA ADS] [CrossRef] [Google Scholar]
 Paumard, T., Perrin, G., Eckart, A., et al. 2008, in The Power of Optical/ IR Interferometry: Recent Scientific Results and 2nd Generation, ed. A. Richichi, F. Delplancke, F. Paresce, & A. Chelli, 313 [Google Scholar]
 Psaltis, D., Perrodin, D., Dienes, K. R., & Mocioiu, I. 2008, Phys. Rev. Lett., 100, 119902 [NASA ADS] [CrossRef] [Google Scholar]
 Qian, L., Abramowicz, M. A., Fragile, P. C., et al. 2009, A&A, 498, 471 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rees, M. J., Begelman, M. C., Blandford, R. D., & Phinney, E. S. 1982, Nature, 295, 17 [NASA ADS] [CrossRef] [Google Scholar]
 Rybicki, G. B., & Lightman, A. P. 1986, Radiative Processes in Astrophysics (WileyVCH) [Google Scholar]
 Schödel, R., Ott, T., Genzel, R., et al. 2002, Nature, 419, 694 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Schödel, R., Eckart, A., Mužić, K., et al. 2007, A&A, 462, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Shapiro, S. L., Lightman, A. P., & Eardley, D. M. 1976, ApJ, 204, 187 [NASA ADS] [CrossRef] [Google Scholar]
 Shcherbakov, R. V., Penna, R. F., & McKinney, J. C. 2012, ApJ, accepted [arXiv:1007.4832] [Google Scholar]
 Stepney, S., & Guilbert, P. W. 1983, MNRAS, 204, 1269 [NASA ADS] [Google Scholar]
 Svensson, R. 1982, ApJ, 258, 335 [NASA ADS] [CrossRef] [Google Scholar]
 Telesco, C. M., Davidson, J. A., & Werner, M. W. 1996, ApJ, 456, 541 [NASA ADS] [CrossRef] [Google Scholar]
 Vincent, F. H., Paumard, T., Gourgoulhon, E., & Perrin, G. 2011, Class. Quant. Grav., 28, 225011 [Google Scholar]
 Yuan, F., Quataert, E., & Narayan, R. 2003, ApJ, 598, 301 [NASA ADS] [CrossRef] [Google Scholar]
 Zhao, J.H., Bower, G. C., & Goss, W. M. 2001, ApJ, 547, L29 [NASA ADS] [CrossRef] [Google Scholar]
 Zylka, R., Mezger, P. G., WardThompson, D., Duschl, W. J., & Lesch, H. 1995, A&A, 297, 83 [NASA ADS] [Google Scholar]
All Tables
All Figures
Fig. 1 Meridional cut through a λ = 0.3 (black curve) and a λ = 0.7 torus (grey curve) around a black hole (the black semicircle is the event horizon). Left panel: Schwarzschild case. Right panel: Kerr case with a = 0.9 M, the longdashed semicircle marking the ergosphere. 

In the text 
Fig. 2 Ion (dashed/dotted) and electron (solid/dotdashed) temperature distribution throughout a λ = 0.3 torus (black) and a λ = 0.7 torus (grey) around a Schwarzschild black hole. We assume a central ion to electron temperature ratio of ξ = 0.1. The thin black line marks T_{vir}. 

In the text 
Fig. 3 Radiative cooling rates of a λ = 0.7 torus around a spin a = 0.9M black hole. The contributions from bremsstrahlung (black, dashed), synchrotron (black, solid), Comptonised bremsstrahlung (grey, dotdashed), and Comptonised synchrotron emission (grey, dotted) are shown as a function of radius. 

In the text 
Fig. 4 Impact of spin on an ion torus spectrum: a = 0 (solid blue), 0.5M (dashed red), or 0.9M (dotted green). All other parameters are set to their reference values listed in Table 1. The black observed data are taken from Zhao et al. (2001), Zylka et al. (1995), Marrone et al. (2008) for radio and submm data, Telesco et al. (1996), Cotera et al. (1999), Eckart et al. (2006), Schödel et al. (2007) for far and midinfrared data, Genzel et al. (2003) for nearinfrared data, and Baganoff et al. (2003) for the Xray bow tie. The downward pointing arrows refer to upper values. Note that we do not present the fitted spectra here. 

In the text 
Fig. 5 Impact of central temperature: T_{0}/T_{vir} = 0.01 (solid blue), 0.015 (dashed red), or 0.02 (dotted green). 

In the text 
Fig. 6 Impact of angular momentum: λ = 0.3 (solid blue), 0.45 (dashed red), or 0.6 (dotted green). 

In the text 
Fig. 7 Impact of pressure ratio: β = 0.01 (solid blue), 0.1 (dashed red), or 0.2 (dotted green). 

In the text 
Fig. 8 Impact of temperature ratio: ξ = 0.08 (solid blue), 0.1 (dashed red), or 0.2 (dotted green). 

In the text 
Fig. 9 Impact of inclination: i = 40° (solid blue), 60° (dashed red), or 80° (dotted green). 

In the text 
Fig. 10 Upper left: image of the reference ion torus corresponding to the parameters listed in Table 1, as observed by an observer on Earth, in inverted colours (the darker, the more luminous). The display is intentionally saturated for a better rendering of the less luminous parts of the image. Upper right: same as in the upper left but with the spin parameter increased to a = 0.9M. Lower left: same as in the upper left but with dimensionless angular momentum increased to λ = 0.6. Lower right: same as in the upper left but with the inclination angle decreased to i = 40°. 

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.