Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
\ctexset

abstractname=Abstract \ctexsetfigurename=Figure \ctexsetbibname=References \ctexsettablename=Table

The shadow and gamma-ray bursts of a Schwarzschild black hole in asymptotic safety

Yuxuan Shi111E-mail address: shiyx2280771974@gmail.com, Hongbo Cheng222E-mail address: hbcheng@ecust.edu.cn
Department of Physics,
East China University of Science and Technology,
Shanghai 200237, China
The Shanghai Key Laboratory of Astrophysics,
Shanghai 200234, China
Abstract

We research on the neutrino pair annihilation ν+ν¯e+e+𝜈¯𝜈superscript𝑒superscript𝑒\nu+\overline{\nu}\longrightarrow e^{-}+e^{+}italic_ν + over¯ start_ARG italic_ν end_ARG ⟶ italic_e start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT + italic_e start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT around a massive source in asymptotic safety. Since neutrinos and photons have the same geodesic equation around black holes, we can estimate the radius where the neutrinos will be released by obtaining a series of trajectory curves with various correction values ξ𝜉\xiitalic_ξ. The black hole shadow radius is influenced by the correction parameter ξ𝜉\xiitalic_ξ. The black hole shadow radius decreases with increasing the ξ𝜉\xiitalic_ξ. The ratio Q˙Q˙Newt˙𝑄subscript˙𝑄𝑁𝑒𝑤𝑡\dfrac{\dot{Q}}{\dot{Q}_{Newt}}divide start_ARG over˙ start_ARG italic_Q end_ARG end_ARG start_ARG over˙ start_ARG italic_Q end_ARG start_POSTSUBSCRIPT italic_N italic_e italic_w italic_t end_POSTSUBSCRIPT end_ARG corresponding to the energy deposition per unit time over that in the Newtonian case is derived and calculated. We find that the quantum corrections to the black hole spacetime affect the emitted energy rate ratio for the annihilation. It is interesting that the more considerable quantum effect reduces the ratio value slightly. Although the energy conversion is damped because of the quantum correction, the energy deposition rate is enough during the neutrino-antineutrino annihilation. The corrected annihilation process can become a source of gamma ray burst. We also investigate the derivative dQ˙drd˙𝑄d𝑟\dfrac{\mathrm{d}\dot{Q}}{\mathrm{d}r}divide start_ARG roman_d over˙ start_ARG italic_Q end_ARG end_ARG start_ARG roman_d italic_r end_ARG relating to the star’s radius r𝑟ritalic_r to show that the quantum effect for the black hole will drop the ratio. The more manifest quantum gravity influence leads the weaker neutrino pair annihilation.

PACS number(s):
Keywords: quantum gravity; back hole; Gamma-ray burst;

I Introduction

A lot of efforts from the astrophysical community have been contributed to explain the gamma-ray bursts (GRBs) phenomenon with enough energy source. GRBs are abrupt increases in gamma-ray intensity coming from an identifiable point in the sky. After the Big Bang, it is the most powerful celestial explosion. Investigating the GRBs’ energy source has been a focus of the astrophysics community. One suitable theory is that the heated accretion disk can annihilate neutrinos and antineutrinos into electrons and positrons, νν¯ee+𝜈¯𝜈superscript𝑒superscript𝑒\nu\bar{\nu}\longrightarrow e^{-}e^{+}italic_ν over¯ start_ARG italic_ν end_ARG ⟶ italic_e start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT. This mechanism has the potential to produce GRBs energy [1-14]. A newborn stellar-mass black hole that is accreting matter at supercritical rates powers this system. Originating from the standpoint of the annihilation process, high-energy gamma rays are released when electron-positron couples close to the neutron star’s surface collapse, resulting in a supernova explosion. There is an obvious explosion like this one. The efficiency of annihilation of neutrino-antineutrino couples into electron-positron pairs is more than 4 times that of the Newtonian case in the background structure of a type II supernova, as demonstrated by Ref [15]. This rate of annihilation increases by 30 times on the surface of a collapsing neutron star. A variety of black holes that have been expanded upon by general relativity have been examined recently. For a Kerr black hole with a thin accretion disk, Ref [16, 17] examined the relativistic impact of the annihilation process νν¯ee+𝜈¯𝜈superscript𝑒superscript𝑒\nu\bar{\nu}\longrightarrow e^{-}e^{+}italic_ν over¯ start_ARG italic_ν end_ARG ⟶ italic_e start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT on the energy deposition rate close to the rotation axis; Ref [18, 19] addressed the implications of a strong off-axis. Additional generalizations of general relativity, including higher derivative gravity, Einstein dilation Gauss-Bonnet, charged Galileon, Brans-Dicke, Born-Infeld generalization of Reissner-Nordstrom solution, and Eddington-inspired Born-Infeld, all dramatically increase the energy deposition rate in neutron star and supernova envelopes [20]. When the quintessence field was taken into account, the authors in Ref [21] found that quintessence strengthened the energy of gamma ray bursts, which could perhaps explain how gamma ray bursts originate. The annihilation process of neutrinos around a massive source with a f(R) global magnetic monopole is investigated in this study, and it is found that this sort of black hole can be an even superior supplier of gamma ray bursts than the beforehand kind of explanation [22].

The black holes attract a lot of attentions from physical community for several decades. As solutions to Einstein’s field equation in the general relativity, the black holes are also the end-results of some stars during the process of gravitational collapse [23-25]. In the case of black holes, the classical description of spacetime can not be used to define the surrounding of singularities inside the event horizon [23-25]. In order to overcome the plague from the singular background, it is necessary to generalize the general relativity up to the high energy scale [26-35]. The generalization can lead the quantum gravity which could resolve the unphysical black hole singularities [26-35]. Several proposals such as loop quantum gravity [26-30], Horava-Lifshitz gravity [31-35], etc. have been put forward. Here we are going to focus on the powerful issue named as asymptotic safety scenario (AS) which is utilized to get rid of the singularities of the black holes [35]. The techniques of the functional renormalization group (FRG) was hired to limit the results of the theory at trans-Planckian energies, making the physical quantities to avoid the divergences at all scales [35]. The quantum effects are always shown as corrections to the classical metrics [36-44]. A series of corrected metrics include Schwarzschild [36], Kerr [37], Schwarzschild-AdS [38], Kerr-AdS [39] and Reissner-Nordstrom-AdS [40-44] without the classical singularity. More attentions have been contributed to the RGI-Schwarzschild black holes. The evaporation of the RGI-Schwarzschild black hole was investigated [45, 46], so was the mini-black hole production in colliders [47]. The radial accretion of matter onto a RGI-Schwarzschild black hole was also explored [48-50]. The shadows of black holes have been studied in the AS scenario for quantum gravity [51-54]. Recently, the quantum gravity effects on the radiation properties of a thin accretion disk surrounding a RGI-Schwarzschild black hole within the frame of IR-limit of the AS theory were discussed [55].

Here we will explore the energy deposition rate by the neutrino pair annihilation process around a Schwarzschild black hole in asymptotic safety. The metric of Schwarzschild black hole modified by the quantum effects needs to be studied in different directions. The motion of test particle around RGI-Schwarzschild black hole was described [56]. The electromagnetic ray reflection spectroscopy for the corrected black hole was also considered [25, 57-59]. The quantum gravity effects on radiation characteristics of thin accretion disks around a renormalization group improved Schwarzschild black hole were shown [55]. Besides the black hole features mentioned above, it is important to probe the annihilation energy deposition rate in the background of this kind of black holes to show the implications of asymptotic safety introducing the quantum corrections to the spacetime structures of the gravitational sources. However, to the best of our knowledge, little contributions have been paid for this topic. In this paper, we derive the integral form of the neutrino pair annihilation efficiency based on the Schwarzschild black hole with running Newton coupling G(r)𝐺𝑟G(r)italic_G ( italic_r ) corresponds to the quantum effects. Secondly, we calculate the ratio of total energy deposition to total Newtonian energy deposition for parameter ξ𝜉\xiitalic_ξ encoding the quantum effects on the spacetime geometry With the help of our numerical estimation, we wil discuss the influence from a RGI-Schwarzschild black hole in the IR-limit of the asymptotic safety scenario on the possibility that the astrophysical bodies attract the annihilation process generating the gamma-ray burst. The results and conclusions are listed in the end.

II The effective potential and photon orbits of the RGI-Schwarzschild black hole

We adopt the Schwarzschild metric modified by the quantum gravity in the infrared limit as follow [60-64],

ds2=f(r)dt2dr2f(r)r2(dθ2+sin2θdφ2)dsuperscript𝑠2𝑓𝑟dsuperscript𝑡2dsuperscript𝑟2𝑓𝑟superscript𝑟2dsuperscript𝜃2superscript2𝜃dsuperscript𝜑2\displaystyle\mathrm{d}s^{2}=f(r)\mathrm{d}t^{2}-\dfrac{\mathrm{d}r^{2}}{f(r)}% -r^{2}\left(\mathrm{d}\theta^{2}+\sin^{2}\theta\mathrm{d}\varphi^{2}\right)roman_d italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = italic_f ( italic_r ) roman_d italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - divide start_ARG roman_d italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_f ( italic_r ) end_ARG - italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( roman_d italic_θ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ roman_d italic_φ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) (1)

where

f(r)=12MG(r)r𝑓𝑟12𝑀𝐺𝑟𝑟\displaystyle f(r)=1-\dfrac{2MG(r)}{r}italic_f ( italic_r ) = 1 - divide start_ARG 2 italic_M italic_G ( italic_r ) end_ARG start_ARG italic_r end_ARG (2)

with mass M𝑀Mitalic_M. For the quantum-gravity-corrected Schwarzschild metric in the infrared limit, the running coupling G(r)𝐺𝑟G(r)italic_G ( italic_r ) takes the form [61],

G(r)=G0(1ξr2)𝐺𝑟subscript𝐺01𝜉superscript𝑟2\displaystyle G(r)=G_{0}\left(1-\dfrac{\xi}{r^{2}}\right)italic_G ( italic_r ) = italic_G start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( 1 - divide start_ARG italic_ξ end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) (3)

where ξ𝜉\xiitalic_ξ is a parameter with dimensions of length squared associated to the scale identification between the momentum scale and the radial coordinate. The general function of radial coordinate G(r)𝐺𝑟G(r)italic_G ( italic_r ) as the generalization of gravitational constant, comes as a result of quantum effects. It should be pointed out that the running coupling G(r)𝐺𝑟G(r)italic_G ( italic_r ) contains the parameter ξ𝜉\xiitalic_ξ standing for the renormalization group improvement and the parameter will change the spacetime structure [36, 54]. There are some upper bounds on ξ𝜉\xiitalic_ξ derived from a series of gravitational and nongravitational measurements [54]. We can choose c=G0=1𝑐subscript𝐺01c=G_{0}=1italic_c = italic_G start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1 in natural units. The radius of the event horizon is a root of f(r)=0𝑓𝑟0f(r)=0italic_f ( italic_r ) = 0 [55]. According to the calculation and discussion, the dimensionless parameter ξ¯=ξM2¯𝜉𝜉superscript𝑀2\overline{\xi}=\dfrac{\xi}{M^{2}}over¯ start_ARG italic_ξ end_ARG = divide start_ARG italic_ξ end_ARG start_ARG italic_M start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG and the mass M𝑀Mitalic_M should be limited as [55],

{0ξ¯ξ¯cM>Mccases0¯𝜉subscript¯𝜉𝑐otherwise𝑀subscript𝑀𝑐otherwise\displaystyle\begin{cases}0\leq\overline{\xi}\leq\overline{\xi}_{c}\\ M>M_{c}\\ \end{cases}{ start_ROW start_CELL 0 ≤ over¯ start_ARG italic_ξ end_ARG ≤ over¯ start_ARG italic_ξ end_ARG start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL italic_M > italic_M start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_CELL start_CELL end_CELL end_ROW (4)

where ξ¯c=1627subscript¯𝜉𝑐1627\overline{\xi}_{c}=\dfrac{16}{27}over¯ start_ARG italic_ξ end_ARG start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = divide start_ARG 16 end_ARG start_ARG 27 end_ARG and Mc=2716ξsubscript𝑀𝑐2716𝜉M_{c}=\sqrt{\dfrac{27}{16}\xi}italic_M start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = square-root start_ARG divide start_ARG 27 end_ARG start_ARG 16 end_ARG italic_ξ end_ARG, to keep the metric with two horizons [55]. It is significant to investigate the quantum gravity effects on the annihilation process within the region (4).

It is necessary to elaborate the accretion disk of Schwarzschild black holes in asymptotic safety because only sufficiently hot accretion disks can emit the neutrinos [65-68]. We plan to discuss the geodesics for a particle motion in the background specified by metric (1). We limit the particle motion to the equatorial plane with θ=π2𝜃𝜋2\theta=\dfrac{\pi}{2}italic_θ = divide start_ARG italic_π end_ARG start_ARG 2 end_ARG to write the Lagrangian as the description of the planar motion given by [15, 69],

\displaystyle\mathcal{L}caligraphic_L =12gμνdxμdτdxνdτabsent12subscript𝑔𝜇𝜈dsuperscript𝑥𝜇d𝜏dsuperscript𝑥𝜈d𝜏\displaystyle=\dfrac{1}{2}g_{\mu\nu}\dfrac{\mathrm{d}x^{\mu}}{\mathrm{d}\tau}% \dfrac{\mathrm{d}x^{\nu}}{\mathrm{d}\tau}\hskip 85.35826pt= divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT divide start_ARG roman_d italic_x start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT end_ARG start_ARG roman_d italic_τ end_ARG divide start_ARG roman_d italic_x start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT end_ARG start_ARG roman_d italic_τ end_ARG
=12[f(r)(dtdτ)21f(r)(drdτ)2r2(dφdτ)2]absent12delimited-[]𝑓𝑟superscriptd𝑡d𝜏21𝑓𝑟superscriptd𝑟d𝜏2superscript𝑟2superscriptd𝜑d𝜏2\displaystyle=\dfrac{1}{2}\left[f(r)\left(\dfrac{\mathrm{d}t}{\mathrm{d}\tau}% \right)^{2}-\dfrac{1}{f(r)}\left(\dfrac{\mathrm{d}r}{\mathrm{d}\tau}\right)^{2% }-r^{2}\left(\dfrac{\mathrm{d}\varphi}{\mathrm{d}\tau}\right)^{2}\right]= divide start_ARG 1 end_ARG start_ARG 2 end_ARG [ italic_f ( italic_r ) ( divide start_ARG roman_d italic_t end_ARG start_ARG roman_d italic_τ end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - divide start_ARG 1 end_ARG start_ARG italic_f ( italic_r ) end_ARG ( divide start_ARG roman_d italic_r end_ARG start_ARG roman_d italic_τ end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG roman_d italic_φ end_ARG start_ARG roman_d italic_τ end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] (5)

then the components of particle’s momentum are,

E𝐸\displaystyle Eitalic_E =f(r)dtdτabsent𝑓𝑟d𝑡d𝜏\displaystyle=f(r)\dfrac{\mathrm{d}t}{\mathrm{d}\tau}= italic_f ( italic_r ) divide start_ARG roman_d italic_t end_ARG start_ARG roman_d italic_τ end_ARG (6)
L𝐿\displaystyle Litalic_L =r2dφdτabsentsuperscript𝑟2d𝜑d𝜏\displaystyle=r^{2}\dfrac{\mathrm{d}\varphi}{\mathrm{d}\tau}= italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG roman_d italic_φ end_ARG start_ARG roman_d italic_τ end_ARG (7)

The geodesics introduce E𝐸Eitalic_E and L𝐿Litalic_L are their own constants and the momenta satisfy PμPμ=0subscript𝑃𝜇superscript𝑃𝜇0P_{\mu}P^{\mu}=0italic_P start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_P start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT = 0 for photons [15]. The equation for radial motion is [15],

(drdτ)2+Veff(r)=1b2superscriptd𝑟d𝜏2subscript𝑉𝑒𝑓𝑓𝑟1superscript𝑏2\displaystyle\left(\dfrac{\mathrm{d}r}{\mathrm{d}\tau}\right)^{2}+V_{eff}(r)=% \dfrac{1}{b^{2}}( divide start_ARG roman_d italic_r end_ARG start_ARG roman_d italic_τ end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_V start_POSTSUBSCRIPT italic_e italic_f italic_f end_POSTSUBSCRIPT ( italic_r ) = divide start_ARG 1 end_ARG start_ARG italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG (8)

where the effective potential is,

Veff(r)=1r2[12Mr(1ξr2)]subscript𝑉𝑒𝑓𝑓𝑟1superscript𝑟2delimited-[]12𝑀𝑟1𝜉superscript𝑟2\displaystyle V_{eff}(r)=\dfrac{1}{r^{2}}\left[1-\dfrac{2M}{r}\left(1-\dfrac{% \xi}{r^{2}}\right)\right]italic_V start_POSTSUBSCRIPT italic_e italic_f italic_f end_POSTSUBSCRIPT ( italic_r ) = divide start_ARG 1 end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ 1 - divide start_ARG 2 italic_M end_ARG start_ARG italic_r end_ARG ( 1 - divide start_ARG italic_ξ end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) ] (9)
Table 1: The event horion r+subscript𝑟r_{+}italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT, the photon sphere rpsubscript𝑟𝑝r_{p}italic_r start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT and the impact parameter bpsubscript𝑏𝑝b_{p}italic_b start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT of the photon sphere with the different ξ𝜉\xiitalic_ξ.
ξ=0𝜉0\xi=0italic_ξ = 0 ξ=1/9𝜉19\xi=1/9italic_ξ = 1 / 9 ξ=2/9𝜉29\xi=2/9italic_ξ = 2 / 9 ξ=1/3𝜉13\xi=1/3italic_ξ = 1 / 3 ξ=4/9𝜉49\xi=4/9italic_ξ = 4 / 9 ξ=5/9𝜉59\xi=5/9italic_ξ = 5 / 9
r+subscript𝑟r_{+}italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT 2.000002.000002.000002.00000 1.941021.941021.941021.94102 1.873361.873361.873361.87336 1.792521.792521.792521.79252 1.688061.688061.688061.68806 1.517491.517491.517491.51749
rpsubscript𝑟𝑝r_{p}italic_r start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT 3.000003.000003.000003.00000 2.935532.935532.935532.93553 2.864602.864602.864602.86460 2.785142.785142.785142.78514 2.693752.693752.693752.69375 2.583972.583972.583972.58397
bpsubscript𝑏𝑝b_{p}italic_b start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT 5.196155.196155.196155.19615 5.129755.129755.129755.12975 5.058185.058185.058185.05818 4.980134.980134.980134.98013 4.893564.893564.893564.89356 4.795024.795024.795024.79502
Refer to caption
(a) ξ=1/3𝜉13\xi=1/3italic_ξ = 1 / 3
Refer to caption
(b) ξ=5/9𝜉59\xi=5/9italic_ξ = 5 / 9
Figure 1: The effective potential profiles with M=1𝑀1M=1italic_M = 1 for ξ=13𝜉13\xi=\dfrac{1}{3}italic_ξ = divide start_ARG 1 end_ARG start_ARG 3 end_ARG (left panel) and ξ=59𝜉59\xi=\dfrac{5}{9}italic_ξ = divide start_ARG 5 end_ARG start_ARG 9 end_ARG (right panel). The photon sphere rpsubscript𝑟𝑝r_{p}italic_r start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT has radii indicated by dashed lines.

Here the impact factor is defined as b=LE𝑏𝐿𝐸b=\dfrac{L}{E}italic_b = divide start_ARG italic_L end_ARG start_ARG italic_E end_ARG. In view of the conditions of photon sphere like drdτ=0d𝑟d𝜏0\dfrac{\mathrm{d}r}{\mathrm{d}\tau}=0divide start_ARG roman_d italic_r end_ARG start_ARG roman_d italic_τ end_ARG = 0 and d2rdτ2=0superscriptd2𝑟dsuperscript𝜏20\dfrac{\mathrm{d}^{2}r}{\mathrm{d}\tau^{2}}=0divide start_ARG roman_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_r end_ARG start_ARG roman_d italic_τ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG = 0, Eq.(9) leads,

Veff(r)subscript𝑉𝑒𝑓𝑓𝑟\displaystyle V_{eff}(r)italic_V start_POSTSUBSCRIPT italic_e italic_f italic_f end_POSTSUBSCRIPT ( italic_r ) =1b2absent1superscript𝑏2\displaystyle=\dfrac{1}{b^{2}}= divide start_ARG 1 end_ARG start_ARG italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG (10)
dVeffdrdsubscript𝑉𝑒𝑓𝑓d𝑟\displaystyle\dfrac{\mathrm{d}V_{eff}}{\mathrm{d}r}divide start_ARG roman_d italic_V start_POSTSUBSCRIPT italic_e italic_f italic_f end_POSTSUBSCRIPT end_ARG start_ARG roman_d italic_r end_ARG =0absent0\displaystyle=0= 0 (11)

or equivalently,

rp2superscriptsubscript𝑟𝑝2\displaystyle r_{p}^{2}italic_r start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT =bp2f(rp)absentsuperscriptsubscript𝑏𝑝2𝑓subscript𝑟𝑝\displaystyle=b_{p}^{2}f(r_{p})= italic_b start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_f ( italic_r start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ) (12)
rp3df(r)dr|r=rp\displaystyle r_{p}^{3}\dfrac{\mathrm{d}f(r)}{\mathrm{d}r}\biggl{|}_{r=r_{p}}italic_r start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT divide start_ARG roman_d italic_f ( italic_r ) end_ARG start_ARG roman_d italic_r end_ARG | start_POSTSUBSCRIPT italic_r = italic_r start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT end_POSTSUBSCRIPT =2bp2f2(rp)absent2superscriptsubscript𝑏𝑝2superscript𝑓2subscript𝑟𝑝\displaystyle=2b_{p}^{2}f^{2}\left(r_{p}\right)= 2 italic_b start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_f start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ) (13)

where rpsubscript𝑟𝑝r_{p}italic_r start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT is the radius of photon sphere and bpsubscript𝑏𝑝b_{p}italic_b start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT is the relevant impact factor. According to the discussions on the metric (1), the outer horizon r+subscript𝑟r_{+}italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT, radius of photon sphere rpsubscript𝑟𝑝r_{p}italic_r start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT and impact factor bpsubscript𝑏𝑝b_{p}italic_b start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT relate to the dimensionless parameter ξ𝜉\xiitalic_ξ associated with the quantum gravity and the relations are listed in the Table 1. We plot the effective potential as function of radial coordinate in Figure 1. These figures show that the quantum gravity effect amends the features of the compact source although the shapes of curves corresponding to the different values of parameter ξ𝜉\xiitalic_ξ are similar.By combining the Eq.(7), Eq.(8) and Eq.(9), we have the trajectory equation like [15],

drdφ=±r21b21r2[12Mr(1ξr2)]d𝑟d𝜑plus-or-minussuperscript𝑟21superscript𝑏21superscript𝑟2delimited-[]12𝑀𝑟1𝜉superscript𝑟2\displaystyle\dfrac{\mathrm{d}r}{\mathrm{d}\varphi}=\pm r^{2}\sqrt{\dfrac{1}{b% ^{2}}-\dfrac{1}{r^{2}}\left[1-\dfrac{2M}{r}\left(1-\dfrac{\xi}{r^{2}}\right)% \right]}divide start_ARG roman_d italic_r end_ARG start_ARG roman_d italic_φ end_ARG = ± italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT square-root start_ARG divide start_ARG 1 end_ARG start_ARG italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - divide start_ARG 1 end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ 1 - divide start_ARG 2 italic_M end_ARG start_ARG italic_r end_ARG ( 1 - divide start_ARG italic_ξ end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) ] end_ARG (14)

By setting u1r𝑢1𝑟u\equiv\dfrac{1}{r}italic_u ≡ divide start_ARG 1 end_ARG start_ARG italic_r end_ARG, we can transform (14) into

dudφ=1b2u2[12Mu(1ξu2)]G(u)d𝑢d𝜑1superscript𝑏2superscript𝑢2delimited-[]12𝑀𝑢1𝜉superscript𝑢2𝐺𝑢\displaystyle\dfrac{\mathrm{d}u}{\mathrm{d}\varphi}=\sqrt{\dfrac{1}{b^{2}}-u^{% 2}\left[1-2Mu\left(1-\xi u^{2}\right)\right]}\equiv G(u)divide start_ARG roman_d italic_u end_ARG start_ARG roman_d italic_φ end_ARG = square-root start_ARG divide start_ARG 1 end_ARG start_ARG italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - italic_u start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [ 1 - 2 italic_M italic_u ( 1 - italic_ξ italic_u start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ] end_ARG ≡ italic_G ( italic_u ) (15)

For the sake of distinguishing among the light trajectories, it is useful to establish orbital fractions n=φ2π𝑛𝜑2𝜋n=\dfrac{\varphi}{2\pi}italic_n = divide start_ARG italic_φ end_ARG start_ARG 2 italic_π end_ARG [15, 69]. The direct emission of these rays corresponds to n<34𝑛34n<\dfrac{3}{4}italic_n < divide start_ARG 3 end_ARG start_ARG 4 end_ARG, the lensing ring to 34<n<5434𝑛54\dfrac{3}{4}<n<\dfrac{5}{4}divide start_ARG 3 end_ARG start_ARG 4 end_ARG < italic_n < divide start_ARG 5 end_ARG start_ARG 4 end_ARG, and the photon ring to n>54𝑛54n>\dfrac{5}{4}italic_n > divide start_ARG 5 end_ARG start_ARG 4 end_ARG. These rays intersect the equatorial plane once, twice, and more than twice. The range of the incident parameter b𝑏bitalic_b of the photon ring, lens ring, and direct image for the situations with ξ=13𝜉13\xi=\dfrac{1}{3}italic_ξ = divide start_ARG 1 end_ARG start_ARG 3 end_ARG and ξ=59𝜉59\xi=\dfrac{5}{9}italic_ξ = divide start_ARG 5 end_ARG start_ARG 9 end_ARG is given below,

ξ=13{Direct emission:n<34,b<4.74314andb>6.04424Lensing ring:34<n<54,4.74314<b<4.96569and5.02321<b<6.04424Photon ring:n>54,4.96569<b<5.02321𝜉13casesDirect emission:𝑛34𝑏expectation4.74314and𝑏6.04424Lensing ring:34𝑛544.74314𝑏4.96569and5.02321𝑏6.04424Photon ring:𝑛544.96569𝑏5.02321\displaystyle\xi=\dfrac{1}{3}\Rightarrow\begin{cases}\text{Direct emission:}n<% \dfrac{3}{4},\quad&b<4.74314\quad\text{and}\quad b>6.04424\\ \text{Lensing ring:}\dfrac{3}{4}<n<\dfrac{5}{4},\quad&4.74314<b<4.96569\quad% \text{and}\quad 5.02321<b<6.04424\\ \text{Photon ring:}n>\dfrac{5}{4},\quad&4.96569<b<5.02321\end{cases}italic_ξ = divide start_ARG 1 end_ARG start_ARG 3 end_ARG ⇒ { start_ROW start_CELL Direct emission: italic_n < divide start_ARG 3 end_ARG start_ARG 4 end_ARG , end_CELL start_CELL italic_b < 4.74314 and italic_b > 6.04424 end_CELL end_ROW start_ROW start_CELL Lensing ring: divide start_ARG 3 end_ARG start_ARG 4 end_ARG < italic_n < divide start_ARG 5 end_ARG start_ARG 4 end_ARG , end_CELL start_CELL 4.74314 < italic_b < 4.96569 and 5.02321 < italic_b < 6.04424 end_CELL end_ROW start_ROW start_CELL Photon ring: italic_n > divide start_ARG 5 end_ARG start_ARG 4 end_ARG , end_CELL start_CELL 4.96569 < italic_b < 5.02321 end_CELL end_ROW (16)
ξ=59{Direct emission:n<34,b<4.43019andb>5.95531Lensing ring:34<n<54,4.43019<b<4.76382and4.85521<b<5.95531Photon ring:n>54,4.76382<b<4.85521𝜉59casesDirect emission:𝑛34𝑏expectation4.43019and𝑏5.95531Lensing ring:34𝑛544.43019𝑏4.76382and4.85521𝑏5.95531Photon ring:𝑛544.76382𝑏4.85521\displaystyle\xi=\dfrac{5}{9}\Rightarrow\begin{cases}\text{Direct emission:}n<% \dfrac{3}{4},\quad&b<4.43019\quad\text{and}\quad b>5.95531\\ \text{Lensing ring:}\dfrac{3}{4}<n<\dfrac{5}{4},\quad&4.43019<b<4.76382\quad% \text{and}\quad 4.85521<b<5.95531\\ \text{Photon ring:}n>\dfrac{5}{4},\quad&4.76382<b<4.85521\end{cases}italic_ξ = divide start_ARG 5 end_ARG start_ARG 9 end_ARG ⇒ { start_ROW start_CELL Direct emission: italic_n < divide start_ARG 3 end_ARG start_ARG 4 end_ARG , end_CELL start_CELL italic_b < 4.43019 and italic_b > 5.95531 end_CELL end_ROW start_ROW start_CELL Lensing ring: divide start_ARG 3 end_ARG start_ARG 4 end_ARG < italic_n < divide start_ARG 5 end_ARG start_ARG 4 end_ARG , end_CELL start_CELL 4.43019 < italic_b < 4.76382 and 4.85521 < italic_b < 5.95531 end_CELL end_ROW start_ROW start_CELL Photon ring: italic_n > divide start_ARG 5 end_ARG start_ARG 4 end_ARG , end_CELL start_CELL 4.76382 < italic_b < 4.85521 end_CELL end_ROW (17)
Refer to caption
(a) ξ=1/3𝜉13\xi=1/3italic_ξ = 1 / 3
Refer to caption
(b) ξ=5/9𝜉59\xi=5/9italic_ξ = 5 / 9
Refer to caption
(c) ξ=1/3𝜉13\xi=1/3italic_ξ = 1 / 3
Refer to caption
(d) ξ=5/9𝜉59\xi=5/9italic_ξ = 5 / 9
Figure 2: The way photons behave in the RGI-Schwarzschild spacetime depends on the impact parameter b𝑏bitalic_b.

The Figure 2 describes a series of photon trajectories with coordinates (r,φ)𝑟𝜑(r,\varphi)( italic_r , italic_φ ) on the equatorial plane. In the Figure 2, The solid discs stand for black holes, and the photon orbit is represented by the dashed grey line. When b<bp𝑏subscript𝑏𝑝b<b_{p}italic_b < italic_b start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT, the electromagnetic rays directly enters the black hole. The photons orbit the gravitational source with the constraint b=bp𝑏subscript𝑏𝑝b=b_{p}italic_b = italic_b start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT. The potential shown in Figure 1 causes the light to deflect, and photons pass through and escape from the black hole if the impact parameter is larger, such as b>bp𝑏subscript𝑏𝑝b>b_{p}italic_b > italic_b start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT. Some photons escape the compact body and the remaining ones are either absorbed by the body or go around the source if the light is directed towards the black hole.

By comparing the aforementioned equations with those of the RGI-Schwarzschild black hole, it is clear that the parameter ξ𝜉\xiitalic_ξ decreases the regions of the direct emission, the lensing ring, and the photon ring. The black hole’s accretion disc is about to form and develop. Neutrinos are released from the hotter disc.

III The energy deposition rate by the neutrino annihilation process

We start to consider the energy deposition in the spacetime governed by Eq.(2) and (3). The energy deposition per unit time and per unit volume for the neutrino pair annihilation process is given by [15],

dE(𝐫)dtdV=2KGF2F(r)n(εν)n(εν¯)(εν+εν¯)εν3εν¯3dενdεν¯d𝐸𝐫d𝑡d𝑉2𝐾superscriptsubscript𝐺𝐹2𝐹𝑟double-integral𝑛subscript𝜀𝜈𝑛subscript𝜀¯𝜈subscript𝜀𝜈subscript𝜀¯𝜈superscriptsubscript𝜀𝜈3superscriptsubscript𝜀¯𝜈3differential-dsubscript𝜀𝜈differential-dsubscript𝜀¯𝜈\displaystyle\dfrac{\mathrm{d}E(\mathbf{r})}{\mathrm{d}t\mathrm{d}V}=2KG_{F}^{% 2}F(r)\iint n(\varepsilon_{\nu})n(\varepsilon_{\overline{\nu}})(\varepsilon_{% \nu}+\varepsilon_{\overline{\nu}})\varepsilon_{\nu}^{3}\varepsilon_{\overline{% \nu}}^{3}\mathrm{d}\varepsilon_{\nu}\mathrm{d}\varepsilon_{\overline{\nu}}divide start_ARG roman_d italic_E ( bold_r ) end_ARG start_ARG roman_d italic_t roman_d italic_V end_ARG = 2 italic_K italic_G start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_F ( italic_r ) ∬ italic_n ( italic_ε start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT ) italic_n ( italic_ε start_POSTSUBSCRIPT over¯ start_ARG italic_ν end_ARG end_POSTSUBSCRIPT ) ( italic_ε start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT + italic_ε start_POSTSUBSCRIPT over¯ start_ARG italic_ν end_ARG end_POSTSUBSCRIPT ) italic_ε start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_ε start_POSTSUBSCRIPT over¯ start_ARG italic_ν end_ARG end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT roman_d italic_ε start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT roman_d italic_ε start_POSTSUBSCRIPT over¯ start_ARG italic_ν end_ARG end_POSTSUBSCRIPT (18)

where

K=16π(1±4sin2θW+8sin4θW)𝐾16𝜋plus-or-minus14superscript2subscript𝜃𝑊8superscript4subscript𝜃𝑊\displaystyle K=\dfrac{1}{6\pi}(1\pm 4\sin^{2}\theta_{W}+8\sin^{4}\theta_{W})italic_K = divide start_ARG 1 end_ARG start_ARG 6 italic_π end_ARG ( 1 ± 4 roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT + 8 roman_sin start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT ) (19)

with the Weinberg angle sin2θW=0.23superscript2subscript𝜃𝑊0.23\sin^{2}\theta_{W}=0.23roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT = 0.23. For various neutrino pairs, the forms of Eq.(6) can be chosen as [15],

K(νμ,ν¯μ)=K(ντ,ν¯τ)=16π(14sin2θW+8sin4θW)𝐾subscript𝜈𝜇subscript¯𝜈𝜇𝐾subscript𝜈𝜏subscript¯𝜈𝜏16𝜋14superscript2subscript𝜃𝑊8superscript4subscript𝜃𝑊\displaystyle K(\nu_{\mu},\overline{\nu}_{\mu})=K(\nu_{\tau},\overline{\nu}_{% \tau})=\dfrac{1}{6\pi}\left(1-4\sin^{2}\theta_{W}+8\sin^{4}\theta_{W}\right)italic_K ( italic_ν start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT , over¯ start_ARG italic_ν end_ARG start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT ) = italic_K ( italic_ν start_POSTSUBSCRIPT italic_τ end_POSTSUBSCRIPT , over¯ start_ARG italic_ν end_ARG start_POSTSUBSCRIPT italic_τ end_POSTSUBSCRIPT ) = divide start_ARG 1 end_ARG start_ARG 6 italic_π end_ARG ( 1 - 4 roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT + 8 roman_sin start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT ) (20)

and

K(νe,ν¯e)=16π(1+4sin2θW+8sin4θW)𝐾subscript𝜈𝑒subscript¯𝜈𝑒16𝜋14superscript2subscript𝜃𝑊8superscript4subscript𝜃𝑊\displaystyle K(\nu_{e},\overline{\nu}_{e})=\dfrac{1}{6\pi}\left(1+4\sin^{2}% \theta_{W}+8\sin^{4}\theta_{W}\right)italic_K ( italic_ν start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT , over¯ start_ARG italic_ν end_ARG start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ) = divide start_ARG 1 end_ARG start_ARG 6 italic_π end_ARG ( 1 + 4 roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT + 8 roman_sin start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT ) (21)

respectively [15]. Here the Fermi constant GF=5.29×1044cm2MeV2subscript𝐺𝐹5.29superscript1044𝑐superscript𝑚2𝑀𝑒superscript𝑉2G_{F}=5.29\times 10^{-44}cm^{2}MeV^{-2}italic_G start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT = 5.29 × 10 start_POSTSUPERSCRIPT - 44 end_POSTSUPERSCRIPT italic_c italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_M italic_e italic_V start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT. The angular integration factor is represented by [15],

F(r)𝐹𝑟\displaystyle F(r)italic_F ( italic_r ) =(1𝛀𝝂𝛀𝝂¯)2dΩνdΩν¯absentdouble-integralsuperscript1subscript𝛀𝝂subscript𝛀bold-¯𝝂2differential-dsubscriptΩ𝜈differential-dsubscriptΩ¯𝜈\displaystyle=\iint\left(1-\bm{\Omega_{\nu}}\cdot\bm{\Omega_{\overline{\nu}}}% \right)^{2}\mathrm{d}\Omega_{\nu}\mathrm{d}\Omega_{\overline{\nu}}= ∬ ( 1 - bold_Ω start_POSTSUBSCRIPT bold_italic_ν end_POSTSUBSCRIPT ⋅ bold_Ω start_POSTSUBSCRIPT overbold_¯ start_ARG bold_italic_ν end_ARG end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_d roman_Ω start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT roman_d roman_Ω start_POSTSUBSCRIPT over¯ start_ARG italic_ν end_ARG end_POSTSUBSCRIPT
=2π23(1x)4(x2+4x+5)absent2superscript𝜋23superscript1𝑥4superscript𝑥24𝑥5\displaystyle=\dfrac{2\pi^{2}}{3}(1-x)^{4}\left(x^{2}+4x+5\right)\hskip 28.452% 74pt= divide start_ARG 2 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 3 end_ARG ( 1 - italic_x ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT ( italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 4 italic_x + 5 ) (22)

where

x=sinθr𝑥subscript𝜃𝑟\displaystyle x=\sin\theta_{r}italic_x = roman_sin italic_θ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT (23)

The angle θrsubscript𝜃𝑟\theta_{r}italic_θ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT is between the particle trajectory and the tangent vector to a circle orbit at radius r𝑟ritalic_r. For one kind of neutrino and antineutrino, Ων(Ων¯)subscriptΩ𝜈subscriptΩ¯𝜈\Omega_{\nu}(\Omega_{\overline{\nu}})roman_Ω start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT ( roman_Ω start_POSTSUBSCRIPT over¯ start_ARG italic_ν end_ARG end_POSTSUBSCRIPT ) is the unit direction vector and dΩν(dΩν¯)dsubscriptΩ𝜈dsubscriptΩ¯𝜈\mathrm{d}\Omega_{\nu}(\mathrm{d}\Omega_{\overline{\nu}})roman_d roman_Ω start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT ( roman_d roman_Ω start_POSTSUBSCRIPT over¯ start_ARG italic_ν end_ARG end_POSTSUBSCRIPT ) is a solid angle. At temperature T𝑇Titalic_T, n(εν)𝑛subscript𝜀𝜈n(\varepsilon_{\nu})italic_n ( italic_ε start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT ) and n(εν¯)𝑛subscript𝜀¯𝜈n(\varepsilon_{\overline{\nu}})italic_n ( italic_ε start_POSTSUBSCRIPT over¯ start_ARG italic_ν end_ARG end_POSTSUBSCRIPT ) are number densities for neutrino and antineutrino respectively in the phase space and satisfy the Fermi-Dirac distribution [15],

n(εν)=2h31exp(ενkT)+1𝑛subscript𝜀𝜈2superscript31subscript𝜀𝜈𝑘𝑇1\displaystyle n(\varepsilon_{\nu})=\dfrac{2}{h^{3}}\dfrac{1}{\exp\left({\dfrac% {\varepsilon_{\nu}}{kT}}\right)+1}italic_n ( italic_ε start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT ) = divide start_ARG 2 end_ARG start_ARG italic_h start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG divide start_ARG 1 end_ARG start_ARG roman_exp ( divide start_ARG italic_ε start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT end_ARG start_ARG italic_k italic_T end_ARG ) + 1 end_ARG (24)

where hhitalic_h is Planck constant and k𝑘kitalic_k is Boltzmann constant. With integrating the Eq.(5), the expression of deposition energy per unit time and unit volume is given by [15],

dEdtdV=21ζ(5)π4h6KGF2F(r)(kT)9d𝐸d𝑡d𝑉21𝜁5superscript𝜋4superscript6𝐾superscriptsubscript𝐺𝐹2𝐹𝑟superscript𝑘𝑇9\displaystyle\dfrac{\mathrm{d}E}{\mathrm{d}t\mathrm{d}V}=\dfrac{21\zeta(5)\pi^% {4}}{h^{6}}KG_{F}^{2}F(r)(kT)^{9}divide start_ARG roman_d italic_E end_ARG start_ARG roman_d italic_t roman_d italic_V end_ARG = divide start_ARG 21 italic_ζ ( 5 ) italic_π start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG italic_h start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT end_ARG italic_K italic_G start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_F ( italic_r ) ( italic_k italic_T ) start_POSTSUPERSCRIPT 9 end_POSTSUPERSCRIPT (25)

It is significant to derive the expression of dEdtdVd𝐸d𝑡d𝑉\dfrac{\mathrm{d}E}{\mathrm{d}t\mathrm{d}V}divide start_ARG roman_d italic_E end_ARG start_ARG roman_d italic_t roman_d italic_V end_ARG which is used to further the research on the converted energy rate on the different compact bodies [15]. The expression has something to do with the position, so has the temperature T=T(r)𝑇𝑇𝑟T=T(r)italic_T = italic_T ( italic_r ) called local temperature [15].

The local temperature T(𝐫)𝑇𝐫T(\mathbf{r})italic_T ( bold_r ) measured by a local observer is defined as T(𝐫)g00(𝐫)=constant𝑇𝐫subscript𝑔00𝐫constantT(\mathbf{r})\sqrt{g_{00}(\mathbf{r})}=\text{constant}italic_T ( bold_r ) square-root start_ARG italic_g start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT ( bold_r ) end_ARG = constant with g00subscript𝑔00g_{00}italic_g start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT, a component of spacetime metric [15]. The neutrino temperature at the neutrinosphere can be explained as [15],

T(r)g00(r)=T(R)g00(R)𝑇𝑟subscript𝑔00𝑟𝑇𝑅subscript𝑔00𝑅\displaystyle T(r)\sqrt{g_{00}(r)}=T(R)\sqrt{g_{00}(R)}italic_T ( italic_r ) square-root start_ARG italic_g start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT ( italic_r ) end_ARG = italic_T ( italic_R ) square-root start_ARG italic_g start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT ( italic_R ) end_ARG (26)

where R𝑅Ritalic_R is the radius of a gravitational source. It is elegant to replace the local temperature T(r)𝑇𝑟T(r)italic_T ( italic_r ) in the future calculation according to the identity (13). The luminosity relating to the redshift can be selected as [15],

L=g00(R0)L(R0)subscript𝐿subscript𝑔00subscript𝑅0𝐿subscript𝑅0\displaystyle L_{\infty}=g_{00}(R_{0})L(R_{0})italic_L start_POSTSUBSCRIPT ∞ end_POSTSUBSCRIPT = italic_g start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT ( italic_R start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) italic_L ( italic_R start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) (27)

where the luminosity for a single neutrino species at the neutrinosphere is [15],

L(R)=4πR0274ac4T4(R)𝐿𝑅4𝜋superscriptsubscript𝑅0274𝑎𝑐4superscript𝑇4𝑅\displaystyle L(R)=4\pi R_{0}^{2}\dfrac{7}{4}\dfrac{ac}{4}T^{4}(R)italic_L ( italic_R ) = 4 italic_π italic_R start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG 7 end_ARG start_ARG 4 end_ARG divide start_ARG italic_a italic_c end_ARG start_ARG 4 end_ARG italic_T start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT ( italic_R ) (28)

where a𝑎aitalic_a is the radiation constant and c𝑐citalic_c is the speed of light in vacuum. In order to replace the temperature associated with the observer’s position, we combine Eq.(26), Eq.(27) and Eq.(28) and substitute them into Eq.(25) to obtain [15],

dE(𝐫)dtdV=21ζ(5)π4h6KGF2k9(74πac)94L94F(r)[g00(R)]94[g00(r)]94R092d𝐸𝐫d𝑡d𝑉21𝜁5superscript𝜋4superscript6𝐾superscriptsubscript𝐺𝐹2superscript𝑘9superscript74𝜋𝑎𝑐94superscriptsubscript𝐿94𝐹𝑟superscriptdelimited-[]subscript𝑔00𝑅94superscriptdelimited-[]subscript𝑔00𝑟94superscriptsubscript𝑅092\displaystyle\dfrac{\mathrm{d}E(\mathbf{r})}{\mathrm{d}t\mathrm{d}V}=\dfrac{21% \zeta(5)\pi^{4}}{h^{6}}KG_{F}^{2}k^{9}\left(\dfrac{7}{4}\pi ac\right)^{-\frac{% 9}{4}}L_{\infty}^{\frac{9}{4}}F(r)\dfrac{\left[g_{00}(R)\right]^{\frac{9}{4}}}% {\left[g_{00}(r)\right]^{\frac{9}{4}}}R_{0}^{-\frac{9}{2}}divide start_ARG roman_d italic_E ( bold_r ) end_ARG start_ARG roman_d italic_t roman_d italic_V end_ARG = divide start_ARG 21 italic_ζ ( 5 ) italic_π start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG italic_h start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT end_ARG italic_K italic_G start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_k start_POSTSUPERSCRIPT 9 end_POSTSUPERSCRIPT ( divide start_ARG 7 end_ARG start_ARG 4 end_ARG italic_π italic_a italic_c ) start_POSTSUPERSCRIPT - divide start_ARG 9 end_ARG start_ARG 4 end_ARG end_POSTSUPERSCRIPT italic_L start_POSTSUBSCRIPT ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT divide start_ARG 9 end_ARG start_ARG 4 end_ARG end_POSTSUPERSCRIPT italic_F ( italic_r ) divide start_ARG [ italic_g start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT ( italic_R ) ] start_POSTSUPERSCRIPT divide start_ARG 9 end_ARG start_ARG 4 end_ARG end_POSTSUPERSCRIPT end_ARG start_ARG [ italic_g start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT ( italic_r ) ] start_POSTSUPERSCRIPT divide start_ARG 9 end_ARG start_ARG 4 end_ARG end_POSTSUPERSCRIPT end_ARG italic_R start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - divide start_ARG 9 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT (29)

In addition to the radial coordinate, the metric components for the massive source surface also appear in the expression of deposition energy per unit time and unit volume like Eq.(29). We can calculate the radiation energy power in the background of the gravitational source by means of the deposition energy density over time as Eq.(29). In order to compute the angular integration F(r)𝐹𝑟F(r)italic_F ( italic_r ) from Eq.(22), we should further the discussion on the variable x𝑥xitalic_x in Eq.(23). We follow the procedure of Ref.[15] and solve the null geodesic in the spacetime of a spherically symmetric gravitational object [15] to show [20],

x2superscript𝑥2\displaystyle x^{2}italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT =sin2θr|θR=0absentevaluated-atsuperscript2subscript𝜃𝑟subscript𝜃𝑅0\displaystyle=\sin^{2}\theta_{r}|_{\theta_{R}=0}= roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT | start_POSTSUBSCRIPT italic_θ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT = 0 end_POSTSUBSCRIPT
=1R2r2f(r)f(R)absent1superscript𝑅2superscript𝑟2𝑓𝑟𝑓𝑅\displaystyle=1-\dfrac{R^{2}}{r^{2}}\dfrac{f(r)}{f(R)}\hskip 14.22636pt= 1 - divide start_ARG italic_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG italic_f ( italic_r ) end_ARG start_ARG italic_f ( italic_R ) end_ARG (30)

Here f(r)=g00(r)𝑓𝑟subscript𝑔00𝑟f(r)=g_{00}(r)italic_f ( italic_r ) = italic_g start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT ( italic_r ), a component of metric (2). It is useful to relate the variable x=sinθ|θR=0𝑥evaluated-at𝜃subscript𝜃𝑅0x=\sin\theta|_{\theta_{R}}=0italic_x = roman_sin italic_θ | start_POSTSUBSCRIPT italic_θ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT end_POSTSUBSCRIPT = 0 to the environment structure around the gravitational object [20]. According to Eq.(22), the angular integration factor become the function of the metric. We can proceed the integration of rate per unit time and unit volume from Eq.(29) on the spherically symmetric volume around the gravitational source [21],

Q˙˙𝑄\displaystyle\dot{Q}over˙ start_ARG italic_Q end_ARG =dEg00dtabsentd𝐸subscript𝑔00d𝑡\displaystyle=\dfrac{\mathrm{d}E}{\sqrt{g_{00}}\mathrm{d}t}\hskip 284.52756pt= divide start_ARG roman_d italic_E end_ARG start_ARG square-root start_ARG italic_g start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT end_ARG roman_d italic_t end_ARG
=84ζ(5)π5h6KGF2k9(74πac)94L94[g00(R)]94R92R0r2g11(r)F(r)g00(r)𝑑rabsent84𝜁5superscript𝜋5superscript6𝐾superscriptsubscript𝐺𝐹2superscript𝑘9superscript74𝜋𝑎𝑐94superscriptsubscript𝐿94superscriptdelimited-[]subscript𝑔00𝑅94superscript𝑅92superscriptsubscriptsubscript𝑅0superscript𝑟2subscript𝑔11𝑟𝐹𝑟subscript𝑔00𝑟differential-d𝑟\displaystyle=\dfrac{84\zeta(5)\pi^{5}}{h^{6}}KG_{F}^{2}k^{9}\left(\dfrac{7}{4% }\pi ac\right)^{-\frac{9}{4}}L_{\infty}^{\frac{9}{4}}\left[g_{00}(R)\right]^{% \frac{9}{4}}R^{-\frac{9}{2}}\int_{R_{0}}^{\infty}\dfrac{r^{2}\sqrt{-g_{11}(r)}% F(r)}{g_{00}(r)}dr= divide start_ARG 84 italic_ζ ( 5 ) italic_π start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT end_ARG start_ARG italic_h start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT end_ARG italic_K italic_G start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_k start_POSTSUPERSCRIPT 9 end_POSTSUPERSCRIPT ( divide start_ARG 7 end_ARG start_ARG 4 end_ARG italic_π italic_a italic_c ) start_POSTSUPERSCRIPT - divide start_ARG 9 end_ARG start_ARG 4 end_ARG end_POSTSUPERSCRIPT italic_L start_POSTSUBSCRIPT ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT divide start_ARG 9 end_ARG start_ARG 4 end_ARG end_POSTSUPERSCRIPT [ italic_g start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT ( italic_R ) ] start_POSTSUPERSCRIPT divide start_ARG 9 end_ARG start_ARG 4 end_ARG end_POSTSUPERSCRIPT italic_R start_POSTSUPERSCRIPT - divide start_ARG 9 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ∫ start_POSTSUBSCRIPT italic_R start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT divide start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT square-root start_ARG - italic_g start_POSTSUBSCRIPT 11 end_POSTSUBSCRIPT ( italic_r ) end_ARG italic_F ( italic_r ) end_ARG start_ARG italic_g start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT ( italic_r ) end_ARG italic_d italic_r (31)

where g11(r)=1f(r)subscript𝑔11𝑟1𝑓𝑟g_{11}(r)=-\dfrac{1}{f(r)}italic_g start_POSTSUBSCRIPT 11 end_POSTSUBSCRIPT ( italic_r ) = - divide start_ARG 1 end_ARG start_ARG italic_f ( italic_r ) end_ARG, also a component of the renormaliztion group improved Schwarzschild metric (2). The metrics of the curved spacetime have a considerable influence on the deposition energy rate. Here Q˙˙𝑄\dot{Q}over˙ start_ARG italic_Q end_ARG can reflect the total amount of energy converted from neutrinos to electron-positron pairs per unit time at any radius [15]. The conversion may become an explosion with the extremely large values of Q˙˙𝑄\dot{Q}over˙ start_ARG italic_Q end_ARG. It is significant to compare the energy deposition rate (31) with the Newtonian quantities as [15, 20, 21],

Q˙Q˙Newt=3[g00(R)]941(x1)4(x2+4x+5)y2g11(Ry)g00(Ry)92𝑑y˙𝑄subscript˙𝑄𝑁𝑒𝑤𝑡3superscriptdelimited-[]subscript𝑔00𝑅94superscriptsubscript1superscript𝑥14superscript𝑥24𝑥5superscript𝑦2subscript𝑔11𝑅𝑦subscript𝑔00superscript𝑅𝑦92differential-d𝑦\displaystyle\dfrac{\dot{Q}}{\dot{Q}_{Newt}}=3\left[g_{00}(R)\right]^{\frac{9}% {4}}\int_{1}^{\infty}(x-1)^{4}\left(x^{2}+4x+5\right)\dfrac{y^{2}\sqrt{-g_{11}% (Ry)}}{g_{00}(Ry)^{\frac{9}{2}}}dydivide start_ARG over˙ start_ARG italic_Q end_ARG end_ARG start_ARG over˙ start_ARG italic_Q end_ARG start_POSTSUBSCRIPT italic_N italic_e italic_w italic_t end_POSTSUBSCRIPT end_ARG = 3 [ italic_g start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT ( italic_R ) ] start_POSTSUPERSCRIPT divide start_ARG 9 end_ARG start_ARG 4 end_ARG end_POSTSUPERSCRIPT ∫ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT ( italic_x - 1 ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT ( italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 4 italic_x + 5 ) divide start_ARG italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT square-root start_ARG - italic_g start_POSTSUBSCRIPT 11 end_POSTSUBSCRIPT ( italic_R italic_y ) end_ARG end_ARG start_ARG italic_g start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT ( italic_R italic_y ) start_POSTSUPERSCRIPT divide start_ARG 9 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT end_ARG italic_d italic_y (32)

with dimensionless variable y=rR𝑦𝑟𝑅y=\dfrac{r}{R}italic_y = divide start_ARG italic_r end_ARG start_ARG italic_R end_ARG and g00(r)subscript𝑔00𝑟g_{00}(r)italic_g start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT ( italic_r ), g11(r)subscript𝑔11𝑟g_{11}(r)italic_g start_POSTSUBSCRIPT 11 end_POSTSUBSCRIPT ( italic_r ) components of metric (1). According to the Ref.[15], we can present dQ˙drd˙𝑄d𝑟\dfrac{\mathrm{d}\dot{Q}}{\mathrm{d}r}divide start_ARG roman_d over˙ start_ARG italic_Q end_ARG end_ARG start_ARG roman_d italic_r end_ARG as the function of radial coordinate r𝑟ritalic_r to exhibit the enhancement,

dQ˙drd˙𝑄d𝑟\displaystyle\dfrac{\mathrm{d}\dot{Q}}{\mathrm{d}r}divide start_ARG roman_d over˙ start_ARG italic_Q end_ARG end_ARG start_ARG roman_d italic_r end_ARG =4π(dEdtdV)g11(r)r2absent4𝜋𝑑𝐸𝑑𝑡𝑑𝑉subscript𝑔11𝑟superscript𝑟2\displaystyle=4\pi\left(\dfrac{dE}{dtdV}\right)\sqrt{-g_{11}(r)}r^{2}\hskip 27% 0.30118pt= 4 italic_π ( divide start_ARG italic_d italic_E end_ARG start_ARG italic_d italic_t italic_d italic_V end_ARG ) square-root start_ARG - italic_g start_POSTSUBSCRIPT 11 end_POSTSUBSCRIPT ( italic_r ) end_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT
=168ζ(5)π73h6KGF2k9(74πac)94L94absent168𝜁5superscript𝜋73superscript6𝐾superscriptsubscript𝐺𝐹2superscript𝑘9superscript74𝜋𝑎𝑐94superscriptsubscript𝐿94\displaystyle=\dfrac{168\zeta(5)\pi^{7}}{3h^{6}}KG_{F}^{2}k^{9}\left(\dfrac{7}% {4}\pi ac\right)^{-\frac{9}{4}}L_{\infty}^{\frac{9}{4}}= divide start_ARG 168 italic_ζ ( 5 ) italic_π start_POSTSUPERSCRIPT 7 end_POSTSUPERSCRIPT end_ARG start_ARG 3 italic_h start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT end_ARG italic_K italic_G start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_k start_POSTSUPERSCRIPT 9 end_POSTSUPERSCRIPT ( divide start_ARG 7 end_ARG start_ARG 4 end_ARG italic_π italic_a italic_c ) start_POSTSUPERSCRIPT - divide start_ARG 9 end_ARG start_ARG 4 end_ARG end_POSTSUPERSCRIPT italic_L start_POSTSUBSCRIPT ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT divide start_ARG 9 end_ARG start_ARG 4 end_ARG end_POSTSUPERSCRIPT
×(x1)4(x2+4x+5)[g00(R)g00(r)]94R52g11(r)(rR)2absentsuperscript𝑥14superscript𝑥24𝑥5superscriptdelimited-[]subscript𝑔00𝑅subscript𝑔00𝑟94superscript𝑅52subscript𝑔11𝑟superscript𝑟𝑅2\displaystyle\quad\times(x-1)^{4}\left(x^{2}+4x+5\right)\left[\dfrac{g_{00}(R)% }{g_{00}(r)}\right]^{\frac{9}{4}}R^{-\frac{5}{2}}\sqrt{-g_{11}(r)}\left(\dfrac% {r}{R}\right)^{2}× ( italic_x - 1 ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT ( italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 4 italic_x + 5 ) [ divide start_ARG italic_g start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT ( italic_R ) end_ARG start_ARG italic_g start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT ( italic_r ) end_ARG ] start_POSTSUPERSCRIPT divide start_ARG 9 end_ARG start_ARG 4 end_ARG end_POSTSUPERSCRIPT italic_R start_POSTSUPERSCRIPT - divide start_ARG 5 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT square-root start_ARG - italic_g start_POSTSUBSCRIPT 11 end_POSTSUBSCRIPT ( italic_r ) end_ARG ( divide start_ARG italic_r end_ARG start_ARG italic_R end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT (33)

Here the derivative is a function of radial coordinate with the origin at the centre of the gravitational source while involving the metric functions. It is necessary to wonder how the structure of the compact body in asymptotic safety affect the neutrino annihilation around one. The derivative function could tell us which kind of stars attracting the annihilation could become source of gamma-ray burst.

It is important to discuss the ratio (32) in the quantum-gravity-corrected Schwarzschild spacetime in the infrared limit. We can relate that g00(r)=f(r)subscript𝑔00𝑟𝑓𝑟g_{00}(r)=f(r)italic_g start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT ( italic_r ) = italic_f ( italic_r ) and g11(r)=1f(r)subscript𝑔11𝑟1𝑓𝑟g_{11}(r)=-\dfrac{1}{f(r)}italic_g start_POSTSUBSCRIPT 11 end_POSTSUBSCRIPT ( italic_r ) = - divide start_ARG 1 end_ARG start_ARG italic_f ( italic_r ) end_ARG for metric (1). The ratio (32) can be reformed as [15, 60-64],

Q˙Q˙Newt=3[f(R)]941(x1)4(x2+4x+5)y2[f(Ry)]5dy˙𝑄subscript˙𝑄𝑁𝑒𝑤𝑡3superscriptdelimited-[]𝑓𝑅94superscriptsubscript1superscript𝑥14superscript𝑥24𝑥5superscript𝑦2superscriptdelimited-[]𝑓𝑅𝑦5differential-d𝑦\displaystyle\dfrac{\dot{Q}}{\dot{Q}_{Newt}}=3\left[f(R)\right]^{\frac{9}{4}}% \int_{1}^{\infty}(x-1)^{4}\left(x^{2}+4x+5\right)\dfrac{y^{2}}{\left[f(Ry)% \right]^{5}}\mathrm{d}ydivide start_ARG over˙ start_ARG italic_Q end_ARG end_ARG start_ARG over˙ start_ARG italic_Q end_ARG start_POSTSUBSCRIPT italic_N italic_e italic_w italic_t end_POSTSUBSCRIPT end_ARG = 3 [ italic_f ( italic_R ) ] start_POSTSUPERSCRIPT divide start_ARG 9 end_ARG start_ARG 4 end_ARG end_POSTSUPERSCRIPT ∫ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT ( italic_x - 1 ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT ( italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 4 italic_x + 5 ) divide start_ARG italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG [ italic_f ( italic_R italic_y ) ] start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT end_ARG roman_d italic_y (34)

where

f(R)𝑓𝑅\displaystyle f(R)italic_f ( italic_R ) =12MR+2MξR3absent12𝑀𝑅2𝑀𝜉superscript𝑅3\displaystyle=1-\dfrac{2M}{R}+\dfrac{2M\xi}{R^{3}}= 1 - divide start_ARG 2 italic_M end_ARG start_ARG italic_R end_ARG + divide start_ARG 2 italic_M italic_ξ end_ARG start_ARG italic_R start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG (35)
f(Ry)𝑓𝑅𝑦\displaystyle f(Ry)italic_f ( italic_R italic_y ) =12MR1y+2MξR31y3absent12𝑀𝑅1𝑦2𝑀𝜉superscript𝑅31superscript𝑦3\displaystyle=1-\dfrac{2M}{R}\dfrac{1}{y}+\dfrac{2M\xi}{R^{3}}\dfrac{1}{y^{3}}= 1 - divide start_ARG 2 italic_M end_ARG start_ARG italic_R end_ARG divide start_ARG 1 end_ARG start_ARG italic_y end_ARG + divide start_ARG 2 italic_M italic_ξ end_ARG start_ARG italic_R start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG divide start_ARG 1 end_ARG start_ARG italic_y start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG (36)

According to the metric function (2) and the running coupling (3), we rewrite the variable (30) as,

x2=11y212MR1y+2MξR31y312MR+2MξR3superscript𝑥211superscript𝑦212𝑀𝑅1𝑦2𝑀𝜉superscript𝑅31superscript𝑦312𝑀𝑅2𝑀𝜉superscript𝑅3\displaystyle x^{2}=1-\dfrac{1}{y^{2}}\dfrac{1-\dfrac{2M}{R}\dfrac{1}{y}+% \dfrac{2M\xi}{R^{3}}\dfrac{1}{y^{3}}}{1-\dfrac{2M}{R}+\dfrac{2M\xi}{R^{3}}}italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 1 - divide start_ARG 1 end_ARG start_ARG italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG 1 - divide start_ARG 2 italic_M end_ARG start_ARG italic_R end_ARG divide start_ARG 1 end_ARG start_ARG italic_y end_ARG + divide start_ARG 2 italic_M italic_ξ end_ARG start_ARG italic_R start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG divide start_ARG 1 end_ARG start_ARG italic_y start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG end_ARG start_ARG 1 - divide start_ARG 2 italic_M end_ARG start_ARG italic_R end_ARG + divide start_ARG 2 italic_M italic_ξ end_ARG start_ARG italic_R start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG end_ARG (37)
Refer to caption
Figure 3: The solid, dotted and dashed curves of the ratio Q˙QNewt˙˙𝑄˙subscript𝑄𝑁𝑒𝑤𝑡\dfrac{\dot{Q}}{\dot{Q_{Newt}}}divide start_ARG over˙ start_ARG italic_Q end_ARG end_ARG start_ARG over˙ start_ARG italic_Q start_POSTSUBSCRIPT italic_N italic_e italic_w italic_t end_POSTSUBSCRIPT end_ARG end_ARG as functions of the ratio RM𝑅𝑀\dfrac{R}{M}divide start_ARG italic_R end_ARG start_ARG italic_M end_ARG for quantum-gravity factors ξ=0,19,13,59𝜉0191359\xi=0,\dfrac{1}{9},\dfrac{1}{3},\dfrac{5}{9}italic_ξ = 0 , divide start_ARG 1 end_ARG start_ARG 9 end_ARG , divide start_ARG 1 end_ARG start_ARG 3 end_ARG , divide start_ARG 5 end_ARG start_ARG 9 end_ARG respectively.
Refer to caption
Figure 4: The solid, dashed and dot-dashed curves of the ratio Q˙QNewt˙˙𝑄˙subscript𝑄𝑁𝑒𝑤𝑡\dfrac{\dot{Q}}{\dot{Q_{Newt}}}divide start_ARG over˙ start_ARG italic_Q end_ARG end_ARG start_ARG over˙ start_ARG italic_Q start_POSTSUBSCRIPT italic_N italic_e italic_w italic_t end_POSTSUBSCRIPT end_ARG end_ARG as functions of the quantum-gravity factors ξ𝜉\xiitalic_ξ with R=4M,4.5M,5M,5.5M𝑅4𝑀4.5𝑀5𝑀5.5𝑀R=4M,4.5M,5M,5.5Mitalic_R = 4 italic_M , 4.5 italic_M , 5 italic_M , 5.5 italic_M respectively.

In order to show the influence from quantum gravity, we should quantify the integral expression of ratio Q˙Q˙Newt˙𝑄subscript˙𝑄𝑁𝑒𝑤𝑡\dfrac{\dot{Q}}{\dot{Q}_{Newt}}divide start_ARG over˙ start_ARG italic_Q end_ARG end_ARG start_ARG over˙ start_ARG italic_Q end_ARG start_POSTSUBSCRIPT italic_N italic_e italic_w italic_t end_POSTSUBSCRIPT end_ARG in Eq.(34) and depict the dependence on R0Msubscript𝑅0𝑀\dfrac{R_{0}}{M}divide start_ARG italic_R start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_M end_ARG with the effect parameter ξ𝜉\xiitalic_ξ in the Figures. Although the metric is corrected by the quantum gravity, it is difficult for us to estimate the influence on the energy deposition rate according to Eq.(34)-(36). We have to perform the complicated calculations to scrutinize the impact of neutrino annihilation associated with the spacetime structure around the compact bodies. It is necessary to gather data for the constraints over the values of ξ𝜉\xiitalic_ξ [54]. During our plotting, we have to magnify itsvalues more than huge times in order to reveal the differences among the curves of the ratio Q˙Q˙Newt˙𝑄subscript˙𝑄𝑁𝑒𝑤𝑡\dfrac{\dot{Q}}{\dot{Q}_{Newt}}divide start_ARG over˙ start_ARG italic_Q end_ARG end_ARG start_ARG over˙ start_ARG italic_Q end_ARG start_POSTSUBSCRIPT italic_N italic_e italic_w italic_t end_POSTSUBSCRIPT end_ARG for different magnitudes of factor ξ𝜉\xiitalic_ξ from quantum-gravity effects on the neutrino pair annihilation around the black hole in asymptotic safety gravity because of extremely tiny values of ξ𝜉\xiitalic_ξ from Ref.[54]. In Figure 1, the shape of ratio curves due to different values of parameter ξ𝜉\xiitalic_ξ are similar. The influence from this kind of sources consist of two parts, mass M𝑀Mitalic_M and gravitational coupling G(r)𝐺𝑟G(r)italic_G ( italic_r ) respectively. As a gravitational source, the values of ratio Q˙Q˙Newt˙𝑄subscript˙𝑄𝑁𝑒𝑤𝑡\dfrac{\dot{Q}}{\dot{Q}_{Newt}}divide start_ARG over˙ start_ARG italic_Q end_ARG end_ARG start_ARG over˙ start_ARG italic_Q end_ARG start_POSTSUBSCRIPT italic_N italic_e italic_w italic_t end_POSTSUBSCRIPT end_ARG are much more than one, so this kind of stellar objects attracting the neutrino pair annihilation process may provide with a source of gamma-ray burst. For the sake of showing the asymptotic behaviour of ratio Q˙Q˙Newt˙𝑄subscript˙𝑄𝑁𝑒𝑤𝑡\dfrac{\dot{Q}}{\dot{Q}_{Newt}}divide start_ARG over˙ start_ARG italic_Q end_ARG end_ARG start_ARG over˙ start_ARG italic_Q end_ARG start_POSTSUBSCRIPT italic_N italic_e italic_w italic_t end_POSTSUBSCRIPT end_ARG more clearly, we plot some curves to indicate that the ratio is decreasing function of ξ𝜉\xiitalic_ξ from quantum corrections in Figure 2. In contrast to the case of accretion disk around a stellar object in asymptotic safety that the quantum gravity effect can enhance the disk’s thermal properties [55], the more considerable influence will reduce the ratio Q˙Q˙Newt˙𝑄subscript˙𝑄𝑁𝑒𝑤𝑡\dfrac{\dot{Q}}{\dot{Q}_{Newt}}divide start_ARG over˙ start_ARG italic_Q end_ARG end_ARG start_ARG over˙ start_ARG italic_Q end_ARG start_POSTSUBSCRIPT italic_N italic_e italic_w italic_t end_POSTSUBSCRIPT end_ARG.

We also elaborate the derivative dQ˙drd˙𝑄d𝑟\dfrac{\mathrm{d}\dot{Q}}{\mathrm{d}r}divide start_ARG roman_d over˙ start_ARG italic_Q end_ARG end_ARG start_ARG roman_d italic_r end_ARG as a function of radius for several stellar masses denoted as rR𝑟𝑅\dfrac{r}{R}divide start_ARG italic_r end_ARG start_ARG italic_R end_ARG with R=3M𝑅3𝑀R=3Mitalic_R = 3 italic_M to exhibit the promotion of ee+superscript𝑒superscript𝑒e^{-}e^{+}italic_e start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT pair energy from the neutrino annihilation in Figure 3. It is similar that a slight larger variable ξ𝜉\xiitalic_ξ from quantum gravity decreases the derivative dQ˙drd˙𝑄d𝑟\dfrac{\mathrm{d}\dot{Q}}{\mathrm{d}r}divide start_ARG roman_d over˙ start_ARG italic_Q end_ARG end_ARG start_ARG roman_d italic_r end_ARG. The Figure 3 also indicate that the structure with smaller rR𝑟𝑅\dfrac{r}{R}divide start_ARG italic_r end_ARG start_ARG italic_R end_ARG lead to the greater dQ˙drd˙𝑄d𝑟\dfrac{\mathrm{d}\dot{Q}}{\mathrm{d}r}divide start_ARG roman_d over˙ start_ARG italic_Q end_ARG end_ARG start_ARG roman_d italic_r end_ARG. It is interesting that the increase is much larger near the surface of the neutron star similar to the case of Ref.[15] although the compact source surrounded by the annihilation includes the quantum gravity influence.

Refer to caption
Figure 5: The solid, dashed, dot-dashed and dotted curves of the derivative dQ˙drd˙𝑄d𝑟\dfrac{\mathrm{d}\dot{Q}}{\mathrm{d}r}divide start_ARG roman_d over˙ start_ARG italic_Q end_ARG end_ARG start_ARG roman_d italic_r end_ARG as functions of the ratio rR𝑟𝑅\dfrac{r}{R}divide start_ARG italic_r end_ARG start_ARG italic_R end_ARG with R=3M𝑅3𝑀R=3Mitalic_R = 3 italic_M for quantum-gravity factors ξ=0,19,13,59𝜉0191359\xi=0,\dfrac{1}{9},\dfrac{1}{3},\dfrac{5}{9}italic_ξ = 0 , divide start_ARG 1 end_ARG start_ARG 9 end_ARG , divide start_ARG 1 end_ARG start_ARG 3 end_ARG , divide start_ARG 5 end_ARG start_ARG 9 end_ARG respectively.

IV Conclusion

The neutrino pair annihilation νν¯ee+𝜈¯𝜈superscript𝑒superscript𝑒\nu\bar{\nu}\longrightarrow e^{-}e^{+}italic_ν over¯ start_ARG italic_ν end_ARG ⟶ italic_e start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT around a Schwarzschild black hole in asymptotic safety is discussed. The quantum gravity influence on the black holes can be explored in different directions. We derivative and calculate the emitted energy rate ratio as Q˙Q˙Newt˙𝑄subscript˙𝑄𝑁𝑒𝑤𝑡\dfrac{\dot{Q}}{\dot{Q}_{Newt}}divide start_ARG over˙ start_ARG italic_Q end_ARG end_ARG start_ARG over˙ start_ARG italic_Q end_ARG start_POSTSUBSCRIPT italic_N italic_e italic_w italic_t end_POSTSUBSCRIPT end_ARG and the emitted energy rate Q˙˙𝑄\dot{Q}over˙ start_ARG italic_Q end_ARG analyzed in the background of quantum-gravity corrected massive source. It is found that the quantum-gravity effect does not advance the ratio Q˙Q˙Newt˙𝑄subscript˙𝑄𝑁𝑒𝑤𝑡\dfrac{\dot{Q}}{\dot{Q}_{Newt}}divide start_ARG over˙ start_ARG italic_Q end_ARG end_ARG start_ARG over˙ start_ARG italic_Q end_ARG start_POSTSUBSCRIPT italic_N italic_e italic_w italic_t end_POSTSUBSCRIPT end_ARG, but damps Q˙˙𝑄\dot{Q}over˙ start_ARG italic_Q end_ARG, the energy released per unit time, in contrast to the case of accretion disk under the same conditions. The larger ξ𝜉\xiitalic_ξ, the parameter encoding the quantum gravity effects on the spacetime structure, will reduce the energy deposition rate. The ratio Q˙Q˙Newt˙𝑄subscript˙𝑄𝑁𝑒𝑤𝑡\dfrac{\dot{Q}}{\dot{Q}_{Newt}}divide start_ARG over˙ start_ARG italic_Q end_ARG end_ARG start_ARG over˙ start_ARG italic_Q end_ARG start_POSTSUBSCRIPT italic_N italic_e italic_w italic_t end_POSTSUBSCRIPT end_ARG shows numerically that the annihilation ν+ν¯e+e+𝜈¯𝜈superscript𝑒superscript𝑒\nu+\bar{\nu}\longrightarrow e^{-}+e^{+}italic_ν + over¯ start_ARG italic_ν end_ARG ⟶ italic_e start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT + italic_e start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT around the quantum-gravity corrected black hole can proceed as a source for the gamma-ray burst. This kind of annihilation can occur around the Type-II supernova or collapsing neutron stars. Within the region of ξ𝜉\xiitalic_ξ, this kind of annihilation process could be proposed as a well-qualified candidate of gamma-ray bursts.

Acknowledge

This work is partly supported by the Shanghai Key Laboratory of Astrophysics 18DZ2271600.

References

  • [1] D. Eichler, M. Livio, T. Piran, D. N. Schramm, Nature 340(1989)126
  • [2] B. Paczynski, Astrophys. J. 363(1990)218
  • [3] P. Meszaros, M. J. Rees, MNRAS 257(1992)29
  • [4] S. E. Woosley, Astrophys. J. 405(1993)273
  • [5] R. Mochkovitch, M. Hernanz, J. Isern, X. Martin, Nature 361(1993)236
  • [6] M. B. Davies, W. Benz, T. Piran, F. K. Thieleman, Astrophys. J. 431(1994)742
  • [7] M. Jaroszynski, A A, 305(1996)839
  • [8] H. Janka, M. Ruffert, A A, 307(1996)L33
  • [9] M. Ruffert, H. Janka, K. Takahashi, G. Schafer, A A, 319(1997)122
  • [10] M. Ruffert, H. Janka, A A, 338(1998)535
  • [11] R. Popham, S. E. Woosley, C. Fryer, Astrophys. J. 518(1999)356
  • [12] A. MacFadyen, S. E. Woosley, Astrophys. J. 524(1999)262
  • [13] M. Ruffert, H. Janka, A A, 344(1999)573
  • [14] M. A. Aloy, E. Muller, J. M. Ibanez, J. M. Martin, A. MacFadyen, Astrophys. J. 531(2000)L119
  • [15] J. D. Salmonson, J. R. Wilson, Astrophys. J. 517(1999)859
  • [16] K. Asano, T. Fukuyama, Astrophys. J. 531(2000)949
  • [17] K. Asano, T. Fukuyama, Astrophys. J. 546(2001)1019
  • [18] W. A. Miller, N. D. George, A. Kheyfets, J. M. McGhee, Astrophys. J. 583(2003)833
  • [19] R. Birkl, M. A. Aloy, H. Janka, E. Muller, A A 463(2007)51
  • [20] G. Lambiase, L. Mastrototaro, Astrophys. J. 904(2020)1
    A. Prasanna, S. Goswami, Phys. Lett. B526(2002)27
  • [21] G. Lambiase, L. Mastrototaro, Eur. Phys. J. C81(2021)932
  • [22] Y. Shi, H. Cheng, EPL 140(2022)49001
  • [23] T. K. Poddar, S. Goswai, A. K. Mishra, Eur. Phys. J. C83(2023)223
  • [24] S. Weinberg, Gravitation and Cosmology: Principles and Applications of the General Theory of Relativity, John Wiley Sons Inc., 1972
    M. Carmeli, Classical Fields: General Relativity and Gauge Theory, John Wiley Sons, Inc., 1982
  • [25] S. Chandrasekhar, The Mathematical Theory of Blck Holes, Oxford University Press, Inc., 1992
  • [26] A. Ashtekar, Phys. Rev. Lett. 57(1986)2244
  • [27] A. Ashtekar, Phys. Rev. D36(1987)1587
  • [28] C. Rovelli, Living Rev. Relativ. 1(1998)1
  • [29] A. Ashtekar, J. Baez, A. Corichi, K. Krasnov, Phys. Rev. Lett. 80(1998)904
  • [30] A. Perez, Class. Quantum Grav. 20(2003)R43
  • [31] P. Horava, Phys. Rev. D79(2009)084008
  • [32] P. Horava, Phys. Rev. Lett. 102(2009)161301
  • [33] E. Kiritsis, G. Kofinas, Nucl. Phys. B821(2009)467
  • [34] G. Cognola, R. Myrzakulov, L. Sebastiani, S. Vagnozzi, S. Zerbini, Class. Quantum Grav. 33(2016)225014
  • [35] M. Reuter, Phys. Rev. D57(1998)971
  • [36] A. Bonanno, M. Reuter, Phys. Rev. D62(2000)043008
  • [37] M. Reuter, E. Tuiran, Phys. Rev. D83(2011)044041
  • [38] B. Koch, F. Saueressig, Class. Quantum Grav. 31(2014)015006
  • [39] J. M. Pawlowski, D. Stock, Phys. Rev. D98(2018)106008
  • [40] C. Gonzalez, B. Koch, Int. J. Mod. Phys. A31(2016)1650141
  • [41] R. Casadio, S. D. H. Hsu, B. Mirza, Phys. Lett. B695(2011)317
  • [42] G. Kofinas, V. Zarikas, JCAP 1510(2015)069
  • [43] A. Bonano, B. Koch, A. Platania, Found. Phys. 48(2018)1393
  • [44] L. Bosma, B. Knorr, F. Sauressig, Phys. Rev. Lett. 123(2019)101301
  • [45] K. Falls, D. F. Litim, Phys. Rev. D89(2014)084002
  • [46] K. Falls, Int. J. Mod. Phys. A27(2012)1250019
  • [47] B. F. L. Ward, Acta Phys. Polon. B37(2006)1967
  • [48] R. Yang, Phys. Rev. D92(2015)084011
  • [49] M. U. Farooq, A. K. Ahmed, R. Yang, M. Jamil, Chin. Phys. C44(2020)065102
  • [50] F. H. Zuluaga, L. A. Sanchez, Chin. Phys. C45(2021)075102
  • [51] A. Held, R. Gold, A. Eichhorn, JCAP 1906(2019)029
  • [52] X. Lu, Y. Xie, Eur. Phys. J. C79(2019)1016
  • [53] R. Kumar, B. P. Singh, S. G. Ghosh, Ann. Phys. 420(2020)168252
    A. Eichhorn, A. Held, JCAP 2105(2021)073
  • [54] G. Lambiase, F. Scardigli, Phys. Rev. D105(2022)124054
  • [55] F. H. Zuluaga, L. A. Sanchez, Eur. Phys. J. C81(2021)840
  • [56] S. Haroon, M. P. Pavlovic, M. Sossich, A. Wang, Eur. Phys. J. C78(2018)519
  • [57] D. N. Page, K. S. Thorne, Astrophys. J. 191(1974)499
  • [58] B. Zhou, A. Abdikamalov, D. Ayzenberg, C. Bambi, S. Nampalliwar, A. Tripathi, JCAP 2101(2021)047
    S. Vagnozzi, R. Roy, Y. Tsai, L. Visinelli, M. Afrin, A. Allahyari, P. Bambhaniya, D. Dey, S. G. Ghosh, P. S. Joshi, K. Jusufi, M. Khodadi, R. K. Walia, A. Ovgun, C. Bambi, ”Horizon-scale tests of gravity theories and fundamental physics from the Event Horizon Telescope image of Sagittarius Asuperscript𝐴A^{\ast}italic_A start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT”, arXiv: 2205.07787
  • [59] Y. Zhang, M. Zhou, C. Bambi, Eur. Phys. J. C78(2018)376
  • [60] D. Benedetti, P. F. Machado, F. Saueressig, Mod. Phys. Lett. A24(2009)2233
  • [61] Y. Cai, D. A. Easson, JCAP 1009(2010)002
  • [62] K. Falls, C. R. King, D. F. Litim, K. Nikolakopoulos, C. Rahmede, Phys. Rev. D97(2018)086006
  • [63] K. Falls, D. F. Litim, J. Schroder, Phys. Rev. D99(2019)126015
  • [64] Y. Kluth, D. F. Litim, arXiv: 2008.09181
  • [65] V. Cardoso, V. F. Foit, M. Kleban, JCAP 1908(2019)006
  • [66] Z. Li, Y. Piao, Phys. Rev. D100(2019)044023
  • [67] K. A. Bronnikov, R. A. Konoplya, Phys. Rev. D101(2020)064004
  • [68] N. Tsukamoto, T. Kokubu, Phys. Rev. D101(2020)044030
  • [69] Y. Nambu, S. Nota, Y. Sakai, Phys. Rev. D100(2019)064037