Nuclear γray emission from very hot accretion flows
^{1}
MaxPlanckInstitut für Kernphysik, Saupfercheckweg 1, 69117 Heidelberg, Germany
email: ervin.kafexhiu@mpihd.mpg.de
^{2}
Dublin Institute for Advanced Studies, 31 Fitzwilliam Place, Dublin 2, UK
^{3}
National Research Nuclear University MEPhI, Kashirskoje Shosse, 31, 115409 Moscow, Russia
^{4}
Department of Physics and Astronomy, Purdue University, West Lafayette, IN 479072036, USA
^{5}
Astrophysical Big Bang Laboratory, RIKEN, 3510198 Saitama, Japan
^{6}
Space Research Institute of the Russian Academy of Sciences (IKI), 84/32 Profsoyuznaya Str, Moscow 117997, Russia
Received:
25
July
2018
Accepted:
5
February
2019
Optically thin accretion plasmas can reach ion temperatures T_{i} ≥ 10^{10} K and thus trigger nuclear reactions. Using a large nuclear interactions network, we studied the radial evolution of the chemical composition of the accretion flow toward the black hole and computed the emissivity in nuclear γray lines. In the advection dominated accretion flow (ADAF) regime, CNO and heavier nuclei are destroyed before reaching the last stable orbit. The overall luminosity in the deexcitation lines for a solar composition of plasma can be as high as few times 10^{−5} the accretion luminosity (Ṁc^{2}) and can be increased for heavier compositions up to 10^{−3}. The efficiency of transformation of the kinetic energy of the outflow into high energy (≥100 MeV) γrays through the production and decay of π^{0}mesons can be higher, up to 10^{−2} of the accretion luminosity. We show that in the ADAF model up to 15% of the mass of accretion matter can “evaporate” in the form of neutrons.
Key words: accretion, accretion disks / gamma rays: general / nuclear reactions, nucleosynthesis, abundances / stars: black holes
© E. Kafexhiu et al. 2019
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Open Access funding provided by Max Planck Society.
1. Introduction
Hard Xray emission from compact Xray binary systems is a distinct signature of a hot plasma accreting onto a black hole. In an attempt to explain the Xray emission from the stellar mass black hole Cygnus X1, Shapiro et al. (1976) developed an accretion disk model (the SLE model) whose inner region was considerably hotter and geometrically thicker than previous models. Accretion plasma is optically thin and is gaspressure dominated. This model naturally yields a twotemperature plasma in which the electron temperature is about kT_{e} ∼ 100 keV and the ion temperature is one or two orders of magnitude larger. However, later developments have shown that the SLE model is thermally unstable (see, e.g., Abramowicz et al. 1995, 2000) and that the advection of heat plays an important role in the viscous and thermal stability of the plasma (see, e.g., Ichimaru 1977; Paczynski & BisnovatyiKogan 1981; Abramowicz et al. 1988, 1995; Narayan & Popham 1993; Narayan & Yi 1994; Chen et al. 1995).
Another group of very hot and optically thin accretion disk solutions is linked to the socalled advectiondominated accretion flows (ADAF) Narayan & Yi (1994, Narayan & Yi (1995a,b). These models are radiative inefficient because most of the plasma heat is advected into the black hole. In ADAF models, plasma is geometrically thick with a quasispherical flow. The electron temperature for ADAF can be as high as kT_{e} ∼ 300 keV, whereas ion temperature is close to the viral temperature and can be on the order of kT_{i} ∼ 100 MeV.
Nuclear reactions play an important role in such high ion temperature accretion plasmas. These reactions change the plasma chemical composition and produce a characteristic γray emission, (Aharonian & Sunyaev 1984; Aharonian & Syunyaev 1987; Yi & Narayan 1997; Kafexhiu et al. 2012). Because of the high ion temperature, the nuclear spallation processes dominate. Therefore, plasma evolves toward a lighter composition. The γray emission of the accretion plasma is governed by prompt nuclear deexcitation γray lines. Their emissivities depend strongly on nuclear abundances. Therefore, their computation requires a detailed consideration of the chemical evolution of the accretion plasma (see, e.g., Kafexhiu et al. 2018).
When the ion temperature of plasma exceeds 10 MeV, it starts to produce π^{0}mesons. Their decay produces an additional γray emission component at higher energies (see, e.g., Dahlbacka et al. 1974; Kolykhalov & Syunyaev 1979; BisnovatyiKogan et al. 1980; Giovannelli et al. 1982; Dermer 1986; Mahadevan et al. 1997; Yi & Narayan 1997; Oka & Manmoto 2003; Kafexhiu et al. 2018). Pions have a higher threshold energy than nuclear reactions. Pions are produced by energetic particles that populate the high energy tail of the plasma distribution function. Therefore, the γray emission from π^{0}meson decay is sensitive on the high energy tail of the plasma distribution function. Nuclear γray lines, on the other hand, carry information about the plasma temperature because the most important part of their cross section is covered by the central part of the Maxwellian distribution function. Thus, the combined information about the prompt γray lines and the π^{0} → 2γ emissions should allow reconstruction of the plasma energy distribution function for a wide energy range.
In this paper, we present the results of our study of the γray emission related to twotemperature accretion plasma within the framework of the ADAF and SLE models. For this task, we utilize a massive nuclear reaction network and the recent parametrizations of πmeson production cross sections at low (close to the kinematic threshold) energies. We study the radial chemical evolution of the accretion flow and compute the resulting γray emission. We also discuss the efficiency of evaporation of free neutrons from the disk and the maximum amount of the secondary Li, Be, and B nuclei produced in the ADAF accretion plasma.
Fig. 1. Radial profiles of the temperature (top) and nucleon number density (bottom) profiles for the ADAF (left) and SLE (right) accretion disk models. The accretion disk parameters are α = 0.1 and ṁ = 0.1α^{2} (black lines) and α = 0.3 and ṁ = 0.17α^{2} (dashed red lines); the mass m = 10 and β = 0.5. 

Open with DEXTER 
2. Radial evolution of the chemical composition
2.1. Nuclear reaction network
Low density accretion plasma, unlike, for example the central core of a star, is dominated by binary nuclear collisions. For a binary reaction i + j → l + …, the production rate per unit volume of the nucleus l is written as (see, e.g., Fowler et al. 1967)
In this equation n_{i}, n_{j} and n_{l} are the number densities of species i, j, and l, respectively; δ_{ij} = 1 if nucleus i and j are identical, and zero otherwise. The is the thermal reaction rate, which is computed from
where, μ_{ij} = m_{i} m_{j}/(m_{i} + m_{j}) is the reduced mass of the two incoming nuclei, is the i + j → l + … reaction cross section, kT is the temperature in units of energy, and E is the collision energy.
Equation (1) is more convenient to solve in terms of the mass fraction abundance. For an element i, the mass fraction abundance is defined as X_{i} = A_{i} n_{i}/ρ, where A_{i} is the mass number, n_{i} is its number density, and ρ = ∑A_{i} n_{i} is the nucleon number density of the plasma. If the accretion plasma is initially composed of N nuclei, then we have to sum up over all possible channels that produce a given nucleus l. Moreover, nuclear reactions can produce nuclei that were not initially present in the plasma. Therefore, we have to consider all possible nuclear species that can be produced during the plasma lifetime and all their most important nuclear reaction channels. This system forms a socalled nuclear reaction network. For calculating the cross sections we use the nuclear code TALYS (Koning et al. 2005) and experimental data, for more details see Kafexhiu et al. (2012). We note that the recent version of the code TALYS can calculate cross sections for energies up to 1 GeV. This has allowed us to extend the computation of the reaction rates for temperatures on the order of 10^{12} K. In addition, the nuclear reaction network considers p + p → p + n + π^{+} and p + p → D + π^{+} channels, which can significantly change the abundances of neutrons and deuterium in case of a proton plasma with temperatures kT_{i} > 10 MeV; see Kafexhiu et al. (2018) for more details.
The solutions of the nuclear network are obtained by solving the following system of differential equations:
We converted the time derivative into the radial derivative using Ẋ^{l} = vdX^{l}/dR, thus solving X^{l} as a function of the disk radius.
In addition, our nuclear reaction network considers reactions i + j → l^{*} + …, which produce nucleus l in different excited states that lead to different prompt γray lines emission. The cross sections for the most important prompt nuclear γray lines are calculated using TALYS code and experimental data (see, e.g., Kozlovsky et al. 2002; Murphy et al. 2009). The number of γrays produced per unit volume ṅ_{γ} for a specific l^{*} transition, is calculated using Eq. (1). The reaction rate is calculated from the specific cross section. Since the accretion plasma temperatures are kT > 1 MeV, the Doppler broadening from thermal motion of l^{*} is larger than the natural width of the prompt γray line. Therefore, the γray line profile is well described by a Gaussian and its emissivity is given by
In this equation, is the central γray energy obtained by the energy difference of the two transition levels of the nucleus, is the Gaussian broadening, and m_{l} is the mass of the excited nucleus l^{*}.
With the exception of nuclear lines, nuclear collisions produce a γray continuum. For E_{γ} ≫ 10 MeV, the nuclear continuum is composed of mainly nuclear bremsstrahlung and the production and decay of the π^{0}mesons. Calculation of the emissivity ϕ_{γ} for these channels is described in Kafexhiu et al. (2018). To compute the differential luminosity of the disk, we have to integrate ϕ_{γ}(E_{γ}, R) over the entire volume of the disk, i.e., Φ_{γ}(E_{γ}) = ∫ ϕ_{γ}(E_{γ}, R) 2H 2πR dR.
2.2. Principal parameters of the accretion disk
Accretion disk is described by the following set of key parameters: (1) the viscosity parameter α; (2) the scaled black hole mass m = M/M_{⊙}, where M_{⊙} is the solar mass; (3) the scaled accretion rate ṁ = Ṁ/Ṁ_{edd}, where Ṁ_{edd} is the socalled Eddington accretion rate; (4) the ratio of the gas pressure to the total pressure β; (5) the adiabatic index γ; (6) the viscous heat parameter δ, which gives the ratio of electron heating to total heating; and (7) the advection parameter f, which gives the ratio of the advection heating to the total heating (see, e.g., Narayan & Yi 1995b). Apart from the system specific parameters such as the m and ṁ, most of the other parameters are constrained from numerical simulations or are limited by other parameters in a selfconsistent way.
By fixing the values of m, ṁ, α, and β, we can calculate the physical quantities of the plasma such as ion and electron temperatures (T_{i} and T_{e}), plasma nucleon number density (ρ), plasma radial speed (v), and the disk height (H), which are required to solve the nuclear reaction network. For the ADAF accretion, the above quantities are calculated from (see, e.g., Narayan & Yi 1995b)
Fig. 2. Radial chemical evolution (left) and the disk γray luminosity (right). The results are calculated for the ADAF model for parameters m = 10, α = 0.1, ṁ = 10^{−3}, β = 0.5, and the initial solar composition. Left panel: evolution of elements. The light elements and the initial nuclei with mass fraction X_{i} > 10^{−3} are explicitly shown. Right panel: disk differential luminosities (in units of photons/cm^{2} s MeV) of nuclear γrays consisting of resolved γray lines and the continuum composed of superposition of broad nuclear lines, neutron capture radiation, and nuclear bremsstrahlung emission). The curves denoted as (1), (2), and (3) show the disk differential π^{0}decay γray luminosities: curve (1) corresponds to the nominal Maxwellian distribution of particles, curve (2) corresponds to a Maxwellian distribution with a sharp cutoff at E_{k} = 4kT_{i}, whereas curve (3) corresponds to the Maxwellian distribution that has a powerlaw tail with an index p = 2.0 for E_{k} > 4 kT_{i}. The luminosity of the thermal electron bremsstrahlung (red dashed line) is also shown. The advection factor for this model is f = 0.86 and the accretion luminosity is L_{acc} = 1.4 × 10^{36} erg s^{−1}. The integrated luminosity that the electrons radiate away through bremsstrahlung is L_{e} = 1.9 × 10^{35} erg s^{−1}, the luminosity in the nuclear lines and continuum is L_{N} = 3.3 × 10^{31} erg s^{−1}, and the luminosity of radiation from decays of the secondary π^{0}mesons is L_{π} = 9.3 × 10^{33} erg s^{−1}. 

Open with DEXTER 
Fig. 3. Same as Fig. 2 but for the SLE disk model for the parameters m = 10, α = 0.1, ṁ = 10^{−3}, β = 0.5. The accretion luminosity for this model is L_{acc} = 1.4 × 10^{36} erg s^{−1}. The integrated luminosity that electrons radiate away through bremsstrahlung is L_{e} = 9.4 × 10^{35} erg s^{−1}, the luminosity due to nuclear interactions is L_{N} = 5.6 × 10^{30} erg s^{−1}, and the luminosity of the π^{0}decay γrays is negligible. 

Open with DEXTER 
In this case, r = R/R_{S}, where R is the disk radius and R_{S} is the Schwarzschild radius. c_{1} = g(α, ϵ′) × (5 + 2ϵ′)/3α^{2} and c_{3} = 2c_{1}/3, where, ϵ′=ϵ/f = f^{−1} (5/3 − γ)/(γ − 1) and g(α, ϵ′) = [1+18α^{2}/(5+2ϵ′)^{2}]^{1/2} − 1. The plasma ion and electron temperature profiles are computed from the local balance equations q^{+} = q^{adv} + q^{ie} for ions and q^{ie} = q^{−} for electrons. The quantity q^{+} is the plasma heating energy rate due to the viscous dissipation, q^{adv} = f q^{+} is the cooling rate due to advection, and q^{ie} is the cooling rate due to the Coulomb coupling with electrons. Electrons, on the other hand, are heated through Coulomb collisions with ions q^{ie} and are cooled by synchrotron, bremsstrahlung, and Comptonization, the sum of which gives q^{−}; for full details see Narayan & Yi (1995b). The solution of the local balance equations determines the value of f.
For the SLE accretion disk model the above quantities are computed from (see Shapiro et al. 1976)
In this equation, M_{*} = M/3, Ṁ_{∗} = Ṁ/10^{17} g s^{−1}, r_{*} = 2 r and . The radial velocity is obtained from the continuity equation and is given by υ = −Ṁ/(2H 2πRρm_{u}), where m_{u} is the nucleon mass. We note that H is half of the disk height.
2.3. Radial chemical evolution of the accreting flow
In this section, we apply the nuclear reaction network to the ADAF and SLE disk models. Throughout this article, we consider a 10 solar mass black hole, m = 10, that is typical for binary systems and assume an equipartition magnetic field with β = 0.5. The ADAF accretion disk exists below some critical values of the accretion rate and plasma density. Therefore, we choose α = 0.1 and ṁ = 0.1 × α^{2} = 10^{−3} far from the critical values and α = 0.3 and ṁ = 0.17 × α^{2} = 1.53 × 10^{−2} near the critical value. The ADAF and SLE temperature and density profiles for this set of parameters are shown in Fig. 1. In addition, we consider two initial chemical compositions: a solar composition realized when the donor is a normal star and a heavy composition expected when the donor is a Wolf–Rayet type star. More specifically we consider H, ^{4}He, ^{12}C, and ^{16}O with mass fractions X_{p} = X_{α} = 0.3 and X_{C12} = X_{O16} = 0.2.
The radial chemical evolution of the ADAF and SLE accretion disks are shown in Figs. 2–5. Figures 2 and 3 show the results for a solar initial composition and α = 0.1 and ṁ = 10^{−3} parameters. Figures 4 and 5 show the same results using ADAF accretion parameters near the critical values^{1} with α = 0.3 and ṁ = 1.53 × 10^{−2}. The results in Fig. 4 are calculated for a solar initial composition, whereas Fig. 5 corresponds to the heavier (HHeCO) initial composition.
Fig. 4. Same as Fig. 2, but for the parameters m = 10, α = 0.3, ṁ = 1.7 × α^{2} = 1.53 × 10^{−2}, β = 0.5. Left panel: evolution of some important nuclei. Right panel: differential photon luminosity of nuclear γrays (blue curve) and the π^{0}decay γrays (cyan line). The luminosity of the thermal electron bremsstrahlung is also shown (red dashed curve). The advection factor for this model is f = 0.58 and the accretion luminosity is L_{acc} = 2.1 × 10^{37} erg s^{−1}. The total luminosity in the nuclear γray lines and continuum is L_{N} = 1 × 10^{33} erg s^{−1} and the luminosity radiated through the thermal π^{0}meson production is L_{π} = 3.3 × 10^{35} erg s^{−1}. 

Open with DEXTER 
Fig. 5. Same as Fig. 4 but for the initial composition with X_{p} = X_{α} = 0.3 and X_{C12} = X_{O16} = 0.2. The total luminosity in nuclear lines and continuum is L_{N} = 8.8 × 10^{33} erg s^{−1}, while the luminosity in π^{0}decay γrays is L_{π} = 2.9 × 10^{35} erg s^{−1}. 

Open with DEXTER 
3. γray emission of the accretion disk
3.1. γray luminosities
The differential γray production rates (“photon luminosity” in units of photons/(MeV s)) calculated for the ADAF and SLE accretion disk models are shown in Figs. 2–5. The key model parameters used in calculations are described in Sect. 2.3.
The γray luminosity of the accretion disk scales as L ∝ ρ^{2} × R^{3} ∝ mṁ^{2} α^{−2}, where ρ is the plasma density. We note, however, that this scaling is valid for accretion parameters far from their critical values. We let L_{acc} =Ṁc^{2} = ṁ × L_{Edd} be the accretion luminosity, where L_{Edd} is the socalled Eddington luminosity. Let us denote with L_{N} the γray luminosity linked to the nuclear reactions (including the prompt deexcitation lines and continuum) and with L_{π} the γray luminosity from the π^{0}meson decay. For an ADAF accretion plasma with a solar initial composition the calculations give (i) far from the critical regime (i.e., ṁ ≪ 0.17 α^{2} and α = 0.3),
(ii) near the critical regime, where the luminosities approach their maximum but do not exceed
For an initial heavy composition plasma with X_{p} = X_{α} = 0.3 and X_{C12} = X_{O16} = 0.2, these upper limits are changed to
For the SLE model, the luminosity radiated by nuclei for an initial solar composition and α = 0.3 and ṁ = 0.17α^{2}, is written as
while in the initially heavy composition plasma,
In the SLE model, the luminosity of π^{0}decay γrays is negligible because the relatively low ion temperature, kT_{i} ≤ 10 MeV.
3.2. Modified high energy tail of the Maxwellian distribution function
Generally, for the thermal plasma, the energy distribution of particles is silently assumed to be described by the Maxwellian distribution. However, this could be not a correct representation as long as it concerns the subrelativistic ion component of the accretion plasma. The particle distribution well above the temperature can be suppressed compared to the Maxwellian distribution if the relaxation time of the tail E ≫ kT supported by the binary elastic scatterings would be longer than the accretion time. In this case, the high energy tail would be underdeveloped and correspondingly, the γray luminosity suppressed, especially at energies well beyond 100 MeV (Aharonian & Atoyan 1983). On the other hand, even a slight acceleration of ions, for example due to plasma instabilities, can produce a suprathermal component well above the “nominal” tail of the Maxwellian distribution. Thus, the shape of the high energy tail of the accretion plasma distribution function can have dramatic effects on the pion production rate (Aharonian & Atoyan 1983; Kafexhiu et al. 2018.
In Fig. 2, we demonstrate the impact that this effect has on the spectrum of π^{0}decay γrays in the case of the ADAF accretion regime. Aside from the “nominal” Maxwellian distribution, two other proton distributions have been considered: a Maxwellian distribution with a sharp cutoff at 4 kT_{i0} and a Maxwellian distribution with a suprathermal tail with a powerlaw index p = 2.0 for energies above 4 kT_{i}. At energies above a few 100 MeV, the corresponding γray luminosities are significantly suppressed or enhanced compared to the luminosity corresponding to the Maxwellian distribution.
4. Evaporation of neutrons
Neutrons are not bounded by the electric and magnetic fields of the accretion plasma. If their energy – the sum of thermal and rotation energies of the accretion disk – is larger than the gravitational binding energy, neutrons can “evaporate” from the disk and reach the companion star before decaying (Aharonian & Sunyaev 1984; Jean & Guessoum 2001; Guessoum & Jean 2002). The ADAF predicts very hot plasma; therefore neutrons can easily escape from the disk almost isotropically. We calculated the neutron evaporation efficiency η in the framework of the ADAF model by adopting the calculations described in Aharonian & Sunyaev (1984), also considering the reduction of neutrons along the accretion disk. Figure 6 shows η as a function of the accretion radius. We find that the η = f(r) profile is universal and does not depend on the accretion parameters. This follows from the fact that the ADAF ion temperature is equal to the virial temperature independent of the accretion parameters.
Fig. 6. Neutron evaporation efficiency η as a function of the radius calculated for the ADAF model. This function is universal; it does not depend strongly on the details of accretion. 

Open with DEXTER 
The primary source of neutrons in the accretion plasma is the destruction of the CNO and heavier nuclei. In addition, ADAF can produce neutrons through the p + p → p + n + π^{+} channel converting protons into neutrons. The neutrons are primarily produced in the r ≲ 10 region of ADAF. Although this channel is less significant than that of nuclear destruction processes, the p + p → p + n + π^{+} channel ensures that ADAF can produce neutrons efficiently even if its accretion plasma is composed of merely hydrogen. Similar to the π^{0} → 2γ emissivity, the neutron production rate through the process p + p → p + n + π^{+} can vary depending on the shape of the high energy tail of the plasma distribution function.
The mass abundance scales of the neutrons such as X_{n} ∝ ρ/υ ∝ ṁα^{−2}. The proportionality constant is a function of radius r. Figure 7 shows the function X_{n} = f(ṁα^{−2}) for ADAF computed near the last stable orbit r = 3, which gives throughout the disk. The neutron abundance is estimated for two initial compositions: a solar composition and a heavier plasma with X_{p} = X_{α} = 0.3 and X_{C} = X_{O} = 0.2. For comparison, X_{n} is plotted with and without the p + p → p + n + π^{+} contribution.
Fig. 7. Neutron mass fraction abundance at r = 3 as a function of ṁ α^{−2}. The magenta dashdotted line and the cyan dashed line are calculated for the initial solar composition, obtained without and with the p + p → p + n + π^{+} contribution, respectively. The red line and the black dashed line are calculated for higher metallicity plasma with initial mass fraction for hydrogen, helium, carbon, and oxygen X_{p} = X_{α} = 0.3 and X_{C12} = X_{O16} = 0.2. The black dashed line takes into account the contribution of the p + p → p + n + π^{+}. 

Open with DEXTER 
5. Production of light elements
Within the ADAF model, the accretion flow has a positive Bernoulli integral. Thus, this model can produce outflows (e.g., Narayan & Yi 1994; Yi & Narayan 1997; Blandford & Begelman 1999). The outflows can carry some of the accretion material and inject it into the surrounding interstellar medium. The accretion material is richer with light elements such as Li and Be. Therefore, the expelled matter can enrich the interstellar medium with these elements in addition to the primordial (cosmological) source of Li and Be. Thus, the accretion disk outflows can, in principle, contribute to the solution of the socalled ^{7}Li problem (see, e.g., Fields 2011). Figure 8 shows the mass fraction abundance of ^{6}Li, ^{7}Li, and ^{7}Be as a function of , for an ADAF plasma close to the last stable orbit r = 3. The calculations are performed for the solar and for a heavy (X_{p} = X_{α} = 0.3 and X_{C} = X_{O} = 0.2) initial compositions.
Fig. 8. Mass fraction abundances of ^{6}Li (cyan), ^{7}Li (black), and ^{7}Be (red) at r = 3 for ADAF as a function of ṁ α^{−2}. The calculations are performed for two, solar, and “heavy” (X_{p} = X_{α} = 0.3 and X_{C} = X_{O} = 0.2) initial plasma compositions. 

Open with DEXTER 
6. Discussion
6.1. Radial chemical evolution
The radial chemical evolutions shown in Figs. 2–5, demonstrate that in the ADAF model the nuclear reactions play an important role starting from distances of 100 R_{S}. The SLE plasma is colder, and, consequently, the impact of nuclear reactions is limited by the region r ≲ 10. The initial composition of the plasma is significantly modified as the nuclei move toward the inner parts of the accretion disk. In the ADAF model, all heavy nuclei are destroyed before they reach the last stable orbit r = 3. As a consequence, the accretion matter that orbits near the black hole is composed of protons, neutrons, and light elements. In addition to ^{4}He, the destruction timescale of which is longer than the accretion falling time, the plasma contains traces of ^{3}He, D (at the level of 10^{−2}), somewhat less T, as well as nuclei of the LiBeB group (at the level of 10^{−4}).
6.2. γrays
Gamma rays produced in nuclear interactions can be observed between ≈1 MeV and 10 MeV. The region below 1 MeV is dominated by the thermal electron bremsstrahlung. The ADAF and SLE accretion models predict different gammaray spectra, basically because of differences in the temperature profiles. The plasma in the ADAF model is hotter. Therefore, γray lines are expected to be broad because of larger Doppler broadening. Within 1–10 MeV, ΔE_{γ}/E_{γ} ≈ 8%. The spectrum of the ADAF plasma continues above 10 MeV through the π^{0} → 2γ emission. Its emissivity peaks around E_{γ} ≈ 100 MeV and falls off exponentially at higher energies. The SLE plasma, on the other hand, produces narrower and more intensive nuclear lines and its γray spectrum has a cutoff at E_{γ} ≈ 10 MeV.
The ADAF and SLE models have about the same efficiency of converting accretion energy into nuclear γray lines, that is L_{N}/L_{acc} < 5 × 10^{−4}. This efficiency decreases if nuclei are destroyed before reaching the inner parts of the accretion disk, where the temperature and density reach their maximum. Therefore, by increasing the plasma accretion rate or density, the nuclear destruction rates become large and, correspondingly, the γray efficiency drops. The maximum possible efficiency depends on the initial plasma composition. It is larger if the plasma is composed of heavier nuclei. In the ADAF model, the maximum luminosity for a solar initial composition is ≈10^{33} × (m/10) erg s^{−1}. For a heavier composition, for example, with the mass abundances X_{p} = X_{α} = 0.3 and X_{C12} = X_{O16} = 0.2, the luminosity is larger, . We note that this estimate exceeds, by an order of magnitude, the result published by Yi & Narayan (1997).
The thermal π^{0}meson production efficiency in the ADAF model is on the order of L_{π}/L_{acc} < 10^{−2}. Unlike the nuclear γray lines, the π^{0}meson production efficiency has a strong dependence on the shape of the high energy tail of the plasma distribution function. The simultaneous observations of nuclear γray lines and π^{0}decay γrays can offer a unique opportunity to reconstruct the plasma energy distribution function, which can be used to discriminate between different accretion regimes. While, the nuclear γray lines contain information about the central part of the plasma distribution function, the π^{0} → 2γ spectra provide information about its high energy tail.
The detectability of nuclear radiation from hot accretion flows depends on the accretion model. The most favorable conditions can be realized for the ADAF regime when the expected energy fluxes in the 1–10 MeV and 50–200 MeV bands from a source located at the distance d can be as high as 10^{−11}(m/10)(d/1 kpc)^{−2} erg cm^{−2} s^{−1} and 3 × 10^{−9}(m/10)(d/1 kpc)^{−2} erg cm^{−2} s^{−1}, respectively (see Sect. 3.1). The flux sensitivities of the planned γray mission eASTROGAM (de Angelis et al. 2018) from the most promising galactic black hole binaries, such as Cyg X1 (m ≈ 15 and d ≈ 2 kpc) are significantly below the above fluxes, especially around 100 MeV. Therefore, eASTROGAM can provide meaningful probes of the ion temperature and accretion regimes through nuclear γrays from these objects.
6.3. Evaporation of neutrons
Figure 6 shows the neutron evaporation efficiency η in the ADAF model as a function of the accretion disk radius. The η = f(r) profile is universal for r < 100 because it depends mainly on the plasma temperature, which is equal to the virial temperature for a wide range of the accretion disk parameters. We find that about 70% of the neutrons can escape from the accretion disk.
The main source of neutrons are the CNO and heavier nuclei destruction. In addition, ADAF can effectively produce neutrons through the p + p → p + n + π^{+} channel that converts protons into neutrons. Like the p + p → π^{0} → 2γ emission, neutrons also originate from the central parts of the accretion disk (r ≲ 10) and their production rate depends on the shape of the high energy tail of the plasma distribution function. Figure 7 shows the maximum neutron fraction in the ADAF model that was obtained at r = 3 as a function of accretion parameters.
Using the evaporation efficiency and the neutron content at different distances, we find that the ADAF can expel up to 15 percent of the accretion mass flow in the form of neutrons. Some of these neutrons can reach the companion star before decaying. These neutrons are cooled down and then captured by protons in the atmosphere and produce the 2.22 MeV γrays that can serve as a signature of the discussed scenario (see, e.g., Aharonian & Sunyaev 1984; Jean & Guessoum 2001).
For consistency, we included the neutron evaporation process in our nuclear reaction network solution. We found no significant effect on the nuclear abundances and γray luminosity of the disk.
6.4. Production of light elements
Light elements are products of spallation of the CNO and heavier nuclei, as well as the synthesis of, for example, p + α and α + α reactions. Figure 8 shows the mass fraction abundances of ^{6}Li, ^{7}Li, and ^{7}Be as a function of ṁ/α^{2} near the last stable orbit for the ADAF accretion disk. For a fixed radius, these abundances are proportional to ṁ/α^{2}. However, for the accretion rates close to the critical rate, the destruction of ^{6}Li, ^{7}Li, and ^{7}Be becomes important and therefore their abundances are reduced; see Fig. 8. Yi & Narayan (1997) estimated lithium production in ADAF; their results are about an order of magnitude larger than the results we obtained via the nuclear reaction network. The reason for this discrepancy could be related to the neglection of Li destruction in their calculations.
Because of the positive Bernoulli integral, ADAF can expel matter into the surrounding interstellar medium providing an additional source of light elements. Detailed quantitative studies of this process could contribute to the solution of the socalled problem of primordial ^{7}Li (see, e.g., Fields 2011).
7. Summary and conclusions
Using a large nuclear reaction network and the recent studies of the πmeson production cross sections at low energies, we computed the radial evolution of the chemical composition of the accreting matter and the γray luminosities expected within the ADAF and SLE accretion disk models. Nuclear reactions are important within 100 R_{S}. The content of CNO and heavier nuclei decreases as the plasma moves toward the black hole event horizon. In the ADAF scenario, the plasma is sufficiently hot to destroy heavy nuclei before reaching the last stable orbit at r = 3. The efficiency of converting accretion energy into nuclear γray lines is L_{N}/L_{acc} < 5 × 10^{−5} for a solar composition and can be increased up to 10^{−3} for heavier compositions. The ADAF model can produce π^{0}decay γrays with an efficiency as high as L_{π}/L_{acc} < 10^{−2}. This efficiency is sensitive to the high energy tail of the plasma distribution function. The ADAF accretion disk can evaporate about 70% of its neutron content with a rate that can be as large as 15% of the accretion rate. The ADAF model can also produce a considerable amount of elements from the LiBeB group that can be injected into the surrounding interstellar medium and alter its chemical composition.
The ADAF solution is valid for accretion rates ≤0.17 × α^{2} (see, e.g., Abramowicz et al. 1995, 2000). For higher accretion rates the accretion regime switches to the relatively cold and thin solution of the Shakura & Sunyaev (1973) disk.
Acknowledgments
The authors are grateful to the members of the High Energy Astrophysics Theory Group for the very fruitful discussions. The authors are appreciative for the grant from the Russian Science Foundation 16121044, the NSF grant AST1306672, DoE grant DESC0016369, and NASA grant 80NSSC17K0757.
References
 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]
 Abramowicz, M. A., Lasota, J.P., & Igumenshchev, I. V. 2000, MNRAS, 314, 775 [NASA ADS] [CrossRef] [Google Scholar]
 Aharonian, F. A., & Atoyan, A. M. 1983, JETP, 58, 1079 [Google Scholar]
 Aharonian, F. A., & Sunyaev, R. A. 1984, MNRAS, 210, 257 [NASA ADS] [CrossRef] [Google Scholar]
 Aharonian, F. A., & Syunyaev, R. A. 1987, Astrophysics, 27, 413 [NASA ADS] [Google Scholar]
 BisnovatyiKogan, G. S., Khlopov, M. Y., Chechetkin, V. M., & Eramzhyan, R. A. 1980, Sov. Astron., 24, 716 [NASA ADS] [Google Scholar]
 Blandford, R. D., & Begelman, M. C. 1999, MNRAS, 303, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Chen, X., Abramowicz, M. A., Lasota, J.P., Narayan, R., & Yi, I. 1995, ApJ, 443, L61 [NASA ADS] [CrossRef] [Google Scholar]
 Dahlbacka, G. H., Chapline, G. F., & Weaver, T. A. 1974, Nature, 250, 36 [NASA ADS] [CrossRef] [Google Scholar]
 de Angelis, A., Tatischeff, V., Grenier, I. A., et al. 2018, J. High Energy Astrophys., 19, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Dermer, C. D. 1986, A&A, 157, 223 [NASA ADS] [Google Scholar]
 Fields, B. D. 2011, Annu. Rev. Nucl. Part. Sci., 61, 47 [NASA ADS] [CrossRef] [Google Scholar]
 Fowler, W. A., Caughlan, G. R., & Zimmerman, B. A. 1967, ARA&A, 5, 525 [NASA ADS] [CrossRef] [Google Scholar]
 Giovannelli, F., Karakula, S., & Tkaczyk, W. 1982, Acta Astron., 32, 121 [NASA ADS] [Google Scholar]
 Guessoum, N., & Jean, P. 2002, A&A, 396, 157 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ichimaru, S. 1977, ApJ, 214, 840 [NASA ADS] [CrossRef] [Google Scholar]
 Jean, P., & Guessoum, N. 2001, A&A, 378, 509 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kafexhiu, E., Aharonian, F., & Vila, G. 2012, Int. J. Mod. Phys. D, 21, 1250009 [NASA ADS] [CrossRef] [Google Scholar]
 Kafexhiu, E., Aharonian, F., & Barkov, M. 2018, ArXiv eprints [arXiv:1807.06079] [Google Scholar]
 Kolykhalov, P. I., & Syunyaev, R. A. 1979, Sov. Astron., 23, 189 [NASA ADS] [Google Scholar]
 Koning, A. J., Hilaire, S., & Duijvestijn, M. C. 2005, in International Conference on Nuclear Data for Science and Technology, eds. R. C. Haight, M. B. Chadwick, T. Kawano, & P. Talou, AIP Conf. Ser., 769, 1154 [NASA ADS] [CrossRef] [Google Scholar]
 Kozlovsky, B., Murphy, R. J., & Ramaty, R. 2002, ApJS, 141, 523 [NASA ADS] [CrossRef] [Google Scholar]
 Mahadevan, R., Narayan, R., & Krolik, J. 1997, ApJ, 486, 268 [NASA ADS] [CrossRef] [Google Scholar]
 Murphy, R. J., Kozlovsky, B., Kiener, J., & Share, G. H. 2009, ApJS, 183, 142 [NASA ADS] [CrossRef] [Google Scholar]
 Narayan, R., & Popham, R. 1993, Nature, 362, 820 [NASA ADS] [CrossRef] [Google Scholar]
 Narayan, R., & Yi, I.S. 1994, ApJ, 428, L13 [NASA ADS] [CrossRef] [Google Scholar]
 Narayan, R., & Yi, I. 1995a, ApJ, 444, 231 [NASA ADS] [CrossRef] [Google Scholar]
 Narayan, R., & Yi, I. 1995b, ApJ, 452, 710 [NASA ADS] [CrossRef] [Google Scholar]
 Oka, K., & Manmoto, T. 2003, MNRAS, 340, 543 [NASA ADS] [CrossRef] [Google Scholar]
 Paczynski, B., & BisnovatyiKogan, G. 1981, Acta Astron., 31, 283 [NASA ADS] [Google Scholar]
 Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 [NASA ADS] [Google Scholar]
 Shapiro, S. L., Lightman, A. P., & Eardley, D. M. 1976, ApJ, 204, 187 [NASA ADS] [CrossRef] [Google Scholar]
 Yi, I., & Narayan, R. 1997, ApJ, 486, 363 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
Fig. 1. Radial profiles of the temperature (top) and nucleon number density (bottom) profiles for the ADAF (left) and SLE (right) accretion disk models. The accretion disk parameters are α = 0.1 and ṁ = 0.1α^{2} (black lines) and α = 0.3 and ṁ = 0.17α^{2} (dashed red lines); the mass m = 10 and β = 0.5. 

Open with DEXTER  
In the text 
Fig. 2. Radial chemical evolution (left) and the disk γray luminosity (right). The results are calculated for the ADAF model for parameters m = 10, α = 0.1, ṁ = 10^{−3}, β = 0.5, and the initial solar composition. Left panel: evolution of elements. The light elements and the initial nuclei with mass fraction X_{i} > 10^{−3} are explicitly shown. Right panel: disk differential luminosities (in units of photons/cm^{2} s MeV) of nuclear γrays consisting of resolved γray lines and the continuum composed of superposition of broad nuclear lines, neutron capture radiation, and nuclear bremsstrahlung emission). The curves denoted as (1), (2), and (3) show the disk differential π^{0}decay γray luminosities: curve (1) corresponds to the nominal Maxwellian distribution of particles, curve (2) corresponds to a Maxwellian distribution with a sharp cutoff at E_{k} = 4kT_{i}, whereas curve (3) corresponds to the Maxwellian distribution that has a powerlaw tail with an index p = 2.0 for E_{k} > 4 kT_{i}. The luminosity of the thermal electron bremsstrahlung (red dashed line) is also shown. The advection factor for this model is f = 0.86 and the accretion luminosity is L_{acc} = 1.4 × 10^{36} erg s^{−1}. The integrated luminosity that the electrons radiate away through bremsstrahlung is L_{e} = 1.9 × 10^{35} erg s^{−1}, the luminosity in the nuclear lines and continuum is L_{N} = 3.3 × 10^{31} erg s^{−1}, and the luminosity of radiation from decays of the secondary π^{0}mesons is L_{π} = 9.3 × 10^{33} erg s^{−1}. 

Open with DEXTER  
In the text 
Fig. 3. Same as Fig. 2 but for the SLE disk model for the parameters m = 10, α = 0.1, ṁ = 10^{−3}, β = 0.5. The accretion luminosity for this model is L_{acc} = 1.4 × 10^{36} erg s^{−1}. The integrated luminosity that electrons radiate away through bremsstrahlung is L_{e} = 9.4 × 10^{35} erg s^{−1}, the luminosity due to nuclear interactions is L_{N} = 5.6 × 10^{30} erg s^{−1}, and the luminosity of the π^{0}decay γrays is negligible. 

Open with DEXTER  
In the text 
Fig. 4. Same as Fig. 2, but for the parameters m = 10, α = 0.3, ṁ = 1.7 × α^{2} = 1.53 × 10^{−2}, β = 0.5. Left panel: evolution of some important nuclei. Right panel: differential photon luminosity of nuclear γrays (blue curve) and the π^{0}decay γrays (cyan line). The luminosity of the thermal electron bremsstrahlung is also shown (red dashed curve). The advection factor for this model is f = 0.58 and the accretion luminosity is L_{acc} = 2.1 × 10^{37} erg s^{−1}. The total luminosity in the nuclear γray lines and continuum is L_{N} = 1 × 10^{33} erg s^{−1} and the luminosity radiated through the thermal π^{0}meson production is L_{π} = 3.3 × 10^{35} erg s^{−1}. 

Open with DEXTER  
In the text 
Fig. 5. Same as Fig. 4 but for the initial composition with X_{p} = X_{α} = 0.3 and X_{C12} = X_{O16} = 0.2. The total luminosity in nuclear lines and continuum is L_{N} = 8.8 × 10^{33} erg s^{−1}, while the luminosity in π^{0}decay γrays is L_{π} = 2.9 × 10^{35} erg s^{−1}. 

Open with DEXTER  
In the text 
Fig. 6. Neutron evaporation efficiency η as a function of the radius calculated for the ADAF model. This function is universal; it does not depend strongly on the details of accretion. 

Open with DEXTER  
In the text 
Fig. 7. Neutron mass fraction abundance at r = 3 as a function of ṁ α^{−2}. The magenta dashdotted line and the cyan dashed line are calculated for the initial solar composition, obtained without and with the p + p → p + n + π^{+} contribution, respectively. The red line and the black dashed line are calculated for higher metallicity plasma with initial mass fraction for hydrogen, helium, carbon, and oxygen X_{p} = X_{α} = 0.3 and X_{C12} = X_{O16} = 0.2. The black dashed line takes into account the contribution of the p + p → p + n + π^{+}. 

Open with DEXTER  
In the text 
Fig. 8. Mass fraction abundances of ^{6}Li (cyan), ^{7}Li (black), and ^{7}Be (red) at r = 3 for ADAF as a function of ṁ α^{−2}. The calculations are performed for two, solar, and “heavy” (X_{p} = X_{α} = 0.3 and X_{C} = X_{O} = 0.2) initial plasma compositions. 

Open with DEXTER  
In the text 