Observational appearances of a f(R) global monopole black hole illuminated by various accretions

In this paper, we explore three simple models of accretions on a global monopole black hole in f(R) theory, and numerically study the corresponding observational appearances as seen by an observer located at the asymptotic infinity and the certain region out of black hole. For the thin-disk accretion, the results here show that the brighter lensing ring and the darker photon ring that around black hole shadow, always make a small contribution and a negligible contribution to total observed intensity respectively. While, the direct emission of disk contributes a dominant part, and the size of shadow always depends on the disk’s location. For the static and infalling spherical accretions, it turns out that the radiuses of the shadows and photon spheres are always same for both accretions, which implies that the boundary of shadow represents the signature of the spacetime geometry in this case. However, we also find that the brightness of shadow in infalling accretion is darker than that in static case since the Doppler effect is taken into account. In addition, the effect of the global monopole parameter η\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta $$\end{document} and f(R) parameter ψ0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\psi _0$$\end{document} on observational appearances of black hole are clearly emphasized throughout of this paper. Finally, we conclude that black hole shadows and the related rings with some different observable features can be used for us to distinguish black holes from different gravity theories and set the upper limits to the f(R) parameter ψ0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\psi _0$$\end{document}.


Introduction
As a most fascinating prediction of General relativity, black hole has always attracted much interests for a long time [1][2][3][4][5]. In particular, an international collaboration of Event Horizon Telescope (EHT) has recently observed that the first shadow image of a supermassive black hole in the centre of M87 has a e-mail: gpliphys@yeah.net (corresponding author) b e-mail: kjhe94@163.com been captured through the very long baseline interferometer experiment [6][7][8][9][10][11]. And, the second image with inclusion of the magnetic fields and the plasma properties has also been imaged in this year [12,13]. This fact receives a wide cheer of scientific community since it can give rise to some significant breakthroughs on fundamental properties of black hole. From the image of M87, it shows that there are a dark interior at the core of image, and outside of which is a compact asymmetric ring. In general, the dark interior is now called as black hole shadow, while the outside bright ring is known as photon ring.
The light ray emitted from the infinity will be deflected as a result of the strong gravitational field when it passes through the near region of black hole. The photon will asymptotically moved to the bound photon orbit if the curvature of the deflected light approach to that of the critical curve. The bound photon orbit is the so-called photon sphere, which can cast a shadow of black hole. The curvature of the deflected light is referred to the impact parameter, i.e., the critical curve corresponds to the critical impact parameter b p . For Schwarzschild case, the photon sphere and the critical impact parameter are found to be 3M and 3 √ 3M, respectively. So far, the study of shadow of black holes and wormholes has received more and more interests in recent year since the first calculations of shadow . The shadow of the Schwarzschild black hole was first studied by Synge and Luminet, and then Bardeen investigated the D-shape shadow for Kerr black hole in 1973 [73][74][75][76]. In modified gravity theories, the shadows cast by the non-rotating and rotating black holes are studied for different values of the parameter [20,21]. Different from the single shadow, double shadows of black holes and wormholes have also been presented in [16,43]. And, the effect of the Melvin magnetic field and the plasma on the structure of shadow have been carefully addressed, and references therein [52,53]. However, it should be pointed out that there are always various kinds of matters existed in real universe, thereby the near region of black hole would arise various accretions naturally, instead of remain the vacuum. Combined with this facts, when Schwarzschild black hole surrounded by a simple model of spherically symmetric accretion, the properties of image has been carefully studied in 2019, and the results show that the accretion has no influence on the shadow size [77]. In this case, the optical appearance of shadow with two photon spheres in a Hairy Black Hole has been presented in [78], where the smaller photon sphere is considered to be the shadow boundary. For another model of accretion, namely the optically thin and geometrically thin disk, the later studies [79][80][81][82][83][84][85][86][87] show that the size of shadow is closely related to the position of the accretion, and the bright region outside of shadow are always composed of direct emission, lensing rings and photon rings. With this consideration, an optically observational signature of the asymmetric thin-shell wormhole has been thoroughly investigated, which shows that this signature can be used to distinguish wormholes from black holes [88]. And, one also finds that the polarized image of black hole produced by photon coupled to Weyl tensor is very important to understand the thin accretion disk [89]. Hence one can see that, it is still very necessary and important to study the shadow of black holes with inclusion of various models of accretions and the other possible factors.
To explain the accelerated-inflation problem without of dark energy and dark matter, Buchdahl in 1970 proposed a type of modified gravity theory, namely the f (R) gravity theory [90][91][92][93]. Since then, many efforts have devoted to study the f (R) gravity theories [94][95][96][97][98][99][100][101][102][103]. In 2012, one has obtained a global monopole black hole solution in f (R) theory [94]. Effects of the f (R) parameter on the thermodynamic properties and strong gravitational lensing of black hole have been discussed by Man and Cheng [95,96]. And in the context of f (R) theory, one has used the standard procedure of the reconstructive approach to derive a d-dimensional black hole solution with a spherically symmetric string cloud configuration [97]. The space-time geometry generated by a f (R) global monopole black hole has also been investigated in the metric-affine formalism [98]. In addition, the Hawking-Page phase transition in this system are inspected with particular attention in 2019 [99]. Nevertheless, the shadows and rings of the f (R) black hole surrounded by various accretions remains unclear. So, our aim in this paper is to focus on this issue. In particular, under the illumination of virous accretions, we will numerically study the observational appearances as seen by an observer located at the asymptotic infinity and the certain region out of black hole, and reveal the effects of the global monopole parameter and f (R) parameter on observational appearances of black hole.
The organization of the paper is as follows: In Sect. 2, we will first study the effective potential and photon orbits of a f (R) global monopole black hole. Section 3 is devoted to investigate shadows and rings of thin disk emission in this spacetime. In Sect. 4, observational appearances of black hole surrounded by the static and infalling spherical accretions will be carefully discussed. Finally, Sect. 5 ends up with some conclusions and discussions.

The effective potential and photon orbits of f (R) global monopole black hole
At first, we in this section will first to study the effective potential of a f (R) black hole with a global monopole. The Lagrangian with inclusion of the global monopole can be expressed as and the self-coupling triplet of scalar fields presents a global monopole, which is where, η and λ are all the model parameters, h(r ) is a dimensionless function, and x a x a = r 2 . In f (R) gravity theory, the action coupled to a matter field can written as f (R) represents the analytical function of the Ricci scalar and κ = 8π G with G is the Newton gravitational constant. S m is the action associated with the matter fields, which is S m = d 4 x √ −gL. In this case, a spherically symmetric solution has been obtained, it is with The factor ψ 0 originated from the f (R) theory of gravity represents the deviation of standard general relativity. And, the linear correction ψ 0 r 1, which is very different from black holes in de Sitter spacetime. The part with a global monopole 8π Gη 2 is also very small, which is approximately equal to 10 −5 . As the conditions ψ 0 → 0 and η → 0 satisfied, this black hole will naturally reduce to the schwarzschild black hole. By solving the equation B(r ) = 0, we have For simplicity, here the Newton constant G has been set to 1. It is clear that there are two horizons(r + and r − ) of f (R) global monopole black hole, the smaller one is the event horizon r − = r h , another is supposed to be the cosmological horizon r + = r c . In the next, our aim is to study the null geodesic of f (R) global monopole black hole. Generally, here we consider the situation that the motion of a particle located at the equatorial plane(θ = π/2). So, it is easy to obtain the Lagrangian L for a particle in this spacetime, which is where,ẋ μ = ∂ x μ ∂λ is the four-velocity of the photon and λ is the affine parameter. In this spacetime, there exists two Killing fields ∂ t and ∂ φ , which correspond to two conserved quantities, namely the energy E and the orbital angular momentum L, respectively. They are 1 Meanwhile, by considering g μνẋ μẋ ν = 0 for null geodesic, then solving Eq. (7) with the aid of Eq. (8), it yieldṡ Here, b c = |L| E is the impact parameter, the affine parameter λ is redefined as λ/|L|, and the sign ± represents the counterclockwise and clockwise direction of the light ray. By introducing the effective potential V e f f into (11), then we havė where , Considering the photon sphere conditionṙ = 0 andr = 0, V e f f should satisfy the condition , V e f f = 0.
The sign represents the derivative with respect to r . In a four-dimensional symmetric black hole, this condition can be rewrote as In Eq. (15), b p and r p are the impact parameter and radius of the photon sphere, respectively. By choosing different values of η and ψ 0 , the numerical results of the event horizon r h , the cosmological horizon r c , the radius r p and impact parameter b p of photon sphere have been shown in Table 1. From it, one can see that the cosmological horizon r c decreases when the f (R) parameter ψ 0 increased, and it will vanish as ψ 0 = 0. However, the event horizon r h , the radius r p and impact parameter b p of photon sphere are all increases with increase of the η and ψ 0 . This implies that the f (R) global monopole black hole does not only have the cosmological horizon, but also the event horizon is larger than that of Schwarzschild black hole. By taking η = 0.07 and ψ 0 = 0.003 as an example, the effective potential V e f f has been intuitively shown in Fig. 1.
From Fig. 1, one can see that the effective potential V e f f has its value existed only when r is greater than or equal to the value of event horizon r h . At the event horizon, the value of V e f f is zero. From this point to the photon sphere r p , it increases monotonously and gets to the maximum value at the photon sphere. Then, it will gradually decrease to zero with the increase of r from photon sphere to infinity. In this sense, the potential V e f f will naturally produce the potential barrier for the entering light rays. When light rays at the position of observer emitted in the radially inward direction, some of them (namely, b c > b p ) will encounter with the potential barrier and will be reflected by it. This region at where those light rays located corresponds to region 1 in Fig. 1. And, some of them (namely, b c < b p ) will not encounter the potential barrier and will drop into black hole finally. This region is related to region 3 in Fig. 1. Moreover, there still exists an interesting region, namely region 2, in which light rays (namely, b c = b p ) will asymptotically approach to the orbit of the photon sphere, and then circle around black hole infinitely many times. On the other hand, it is obvious that the effective potential V e f f decreases in f (R) global monopole black hole by comparing with Schwarzschild black hole.
To carefully address the trajectories of light ray, here we use the Eqs. (10) and (11) to reexpress the motion equation of photon, it is Table 1 The numerical results of the event horizon r h , the cosmological horizon r c , the radius r p and impact parameter b p of photon sphere, where the value of M is set to 1  With a rearrangement of r , i.e., u = 1/r , the above equation can be expressed as Using the ray-tracing code, the trajectories of light ray are shown in Fig. 2 with the aid of Eq. (17). The spacing of impact parameter b c in Fig. 2 are all chosen as 1 5 . The solid disk represents the black hole, while the dashed grey circle is the photon orbit. When the impact parameter b c < b p (region 3 in Fig. 1), the photon will be captured by black hole. The trajectories of those photons correspond to the black line that shown in Fig. 2. On the contrary, if b c > b p (region 1 in Fig. 1), the photon will be reflected back and move to infinity. And, its trajectories correspond to the green line in Fig. 2. While b c = b p (region 2 in Fig. 1), the photon will revolve around black hole infinitely many times. The trajectories of those photons approximately correspond to the red line in Fig. 2. More importantly, it is clear from Fig. 2 that the incident light rays are almost parallel to the abscissa axis for subfigure (a), while are more or less sloping for subfigures (b) and (c). The reason is that there is no cosmological horizon existed in (a), thereby the observer is assumed to be located at the infinity. But for (b) and (c), due to the existence of cosmological horizon, the location of the observer is fixed to r c 6 in Fig. 2. 2 In addition, the global monopole parameter η and the f (R) parameter ψ 0 produce an obvious influence of trajectories of light rays that cannot be ignored.

Shadows and rings of thin disk emission
For a real universe, it is general considered that black hole surrounded by various matters. In this case, the study of shadows and rings of black hole will be more interesting. Therefore, we in this section will investigate the shadows and rings of f (R) global monopole black hole when it surrounded by an optically and geometrically thin disk accretion.

Direct emission, Lensed ring and photon ring
As described by Gralla et al. [79], when one traced the light ray from the observer backwards towards the black hole, the emissions with different value of b c will give rise to different appearances of the black hole to a distant observer. According to the orbital plane azimuthal angle φ, the total number of orbits n = φ/2π as a function of the impact parameter b c has been defined to distinguish the trajectories of light rays emitted from the north pole direction. There exists three cases for different values of n. They are direct emission, lensing ring and photon ring respectively, which correspond to n < 3/4, 3/4 < n < 5/4 and n > 5/4. And, the corresponding trajectories of light rays will intersect with the equatorial plane once, twice and more than twice, respectively. For a f (R) black hole with the global monopole, we have taken the following choice of η and ψ 0 as the examples to present the regions of direct emission, lensing ring and photon ring in Eqs. (18)- (20).
Photon ring: n > 5/4, 6.42633 < b c < 6.49637 (20) It is obvious from the above equations that two parameters η and ψ 0 all enlarge the regions of direct emission, lensing ring and photon ring by comparing with that of Schwarzschild black hole. In the following Fig. 3, we take a more careful attention to show the trajectories of light ray, where the regions of direct emissions and rings are emphasized.
The first line directly shows the fractional number of orbits n versus the impact parameter b c , where the singularities always occur at b c = b p . And, the corresponding trajectories of light rays in the f (R) global monopole black hole are presented in the second line. In this line, the direct, lensing ring and photon ring trajectories correspond to the red line, blue line and green line, for which the related spacings of b c reads 1/5, 1/50 and 1/500, respectively. Similarly to Fig. 2, the photon orbit is the dashed grey line and black hole is presented as a solid disk.

Observed specific intensities and transfer functions
After discussing direct emissions and rings in the f (R) global monopole black hole, our aim in this subsection is to investigate the observed specific intensity of the thin-disk accretion, which viewed in a face-on orientation. Here, the thin disk is assumed to be lied in the equatorial plane of f (R) global monopole black hole. And, the emissions from it are isotropically in the rest frame of static worldlines. In this case, when the static observer located at the north pole, the observed specific intensity can be wrote as the following form. It is Here, we have defined I obs (r ) as the observed specific intensity and its frequency as ν, while I em (r ) as the emitted specific intensity and its frequency as ν e . The total observed specific intensity can be reexpressed as The relationship ν = A(r ) 1/2 ν e has been used to simplify Eq. (22). The total emitted specific intensity I emi (r ) is an integral of I em (r ) with respect to ν e , i.e., I emi (r ) = I em (r )dν e . Since the thin disk lies in the equatorial plane, the light rays emitted from the north pole will intersect with the disk. So, the brightness will be picked up by the light rays for each intersection. When the photons emitted from thin disk, the light rays located at the regions of lensed rings (which are the blue lines in Fig. 3) will intersect with the opposite side of the disk once. Those rays will pick up the additional brightness once. For the light rays located at the regions of photon rings (which are the green lines in Fig. 3), it does not only intersect with the opposite side of the disk, but also with the front side. Those rays will pick up the additional brightness twice. In view of this, the total observed intensity should be regarded as a sum of those intensities from each intersection. So, The function r n (b c ) is the radial location at which the light rays intersect with the disk. It is referred as the transfer function, where n = 1, 2, 3 . . . . At each point of b c , there will produce a slope of the transfer function dr/db c , which represents the demagnification factor. Note that, the absorption of brightness which may abate the observation intensity has been neglected here. In Fig. 4, we have plotted the first three transfer functions r n (b) with respect to the impact parameter b c for the f (R) global monopole black hole. In Fig. 4, the red, blue and green lines correspond to the first(n = 1), second(n = 2) and third(n = 3) transfer function, respectively. The slope of those lines are very different. For the red line, the value of its slope is in the neighbourhood of 1. So, this line is related to the direct image of the disk, which represents the redshifted source profile. While, the blue line's slope increase quickly and is much greater  Fig. 4 The first three transfer functions for the f (R) global monopole black hole than 1, thereby this line is related to the lensing ring (with the photon ring included in), which represents the image of the back side of the disk that has been largely demagnified. For the green line, its slope is close to infinity. This line is related to the photon ring. So, it implies that the image of the front side of the disk is extremely demagnified. Similar to previous figures, it is easy to see that the transfer function depends on the global monopole parameter η and the f (R) parameter ψ 0 closely.

Observational appearances of direct emissions and rings
Based on above discussions, our aim in this subsection is continue to study the observed specific intensity with some typical toy-model functions. Firstly, we assume that the emissions of the thin disk occur at the innermost stable circular orbit and the emitted function I emi (r ) of it is a decay function suppressed by the second power. In this case, the emitted function I emi (r ) is assumed to be With the help of Eq. (23), the corresponding observed intensity of direct emissions and rings for the f (R) global monopole black hole have been intuitively shown in Fig. 5. In Fig. 5, the left, middle and right column correspond to the emitted function I emi (r ), the observed specific intensity I (r ) and the corresponding observational appearances of the thin disk under various choices of parameters η and ψ 0 . From the left column, it is obvious that the emissions peaked sharply at r ∼ 6.8M for subfigure (a)(r ∼ 6.4M for subfigure (b) and r ∼ 7.7M for subfigure (c)), and then decayed quickly with the increase of r . From middle column, one can see that the observed intensities of direct image are very similar to the emitted intensities, and the peaks of it move to the position b c ∼ 8.2M for the case (a)(b c ∼ 7.7M for the case (b) and b c ∼ 9.6M for the case (c)) due to the gravitational lensing. The lensing ring as an image of the disk's back side is found to in the range of 6.8M ∼ 7.6M for (a)(5.6M ∼ 6.1M for (b) and 6.9M ∼ 7.7M for (c)). This thin region makes a small contribution to the total observed intensity. Also, the photon ring as a more narrow region located at the position b c ∼ 6.3M for (a)(b c ∼ 5.3M for (b) and r ∼ 6.4M for (c)). The intensity of this ring is so small that it is hard to see it in the right column. Also, we find that the parameter η decreased the observed intensity of direct emission, while ψ 0 increased it. For the lensing ring, two parameters all cut down the observed intensity, even if they hardly affect the intensity of photon ring. From the right column, the corresponding shadows and rings can be seen in Fig. 5, where the hardly saw photon rings will appear by enlarging the plot.
Secondly, when the thin disk located at the photon sphere, we assume the emitted function I emi (r ) is an exponential decay function. In this situation, the emitted function I emi (r ) is In a similar way, the corresponding observed intensity of direct emissions and rings can be easily obtained, which have been plotted in Fig. 6. From left column of Fig. 6, the peak of emission occurs at r ∼ 3.4M for case (a)(r ∼ 3.1M for case (b) and r ∼ 3.4M for case (c)). The observed intensity I (r ) of direct emissions peaks near the position b c ∼ 4.5M for case (a)(b c ∼ 4M for case (b) and b c ∼ 4.7M for case (c)). Interestingly, the direct emissions, the lensing ring as well as photon ring in the middle column overlapped at the region 6.3M ∼ 8M for case (a)(5.3M ∼ 6M for case (b) and 6.4M ∼ 7.2M for case (c)). However, it shows that the lensing ring and photon ring are indistinguishable for case (a) and (b) while distinguishable for case (c). Similar to the function I emi (r ), the lensing ring also contributes a small part to the total observed intensity, and photon ring remains an ignorable part. Meanwhile, we find that the parameters η and ψ 0 all decreased observed intensities of direct emission and rings. The related shadows and rings are clearly presented at the right column of Fig. 6.
Thirdly, by considering the thin disk lies in the position of the event horizon, we assume I emi (r ) is a more or less moderate decay function. In this consideration, the emitted function I emi (r ) can be read as Similarly, we in Fig. 7 have presented the corresponding shadows and rings of this geometrically and optically thin disk.
In Fig. 7, the emission of the disk decayed slowly and spread to the position r ∼ 2.3M for case (a)(r ∼ 2.1M for case (b) and r ∼ 2.3M for case (c)). In this sense, the observed intensity of direct emission occurs at b c ∼ 3.4M for case (a)(b c ∼ 3M for case (b) and b c ∼ 3.5M for case (c)). It shows in the middle column that the case (a) is very similar to that of the Schwarzschild black hole, where the lensing ring and photon ring are indistinguishable and they superimposed on the direct image. But unlike Schwarzschild black hole, there is only a peak for case (b) since the peak of direct emission and rings are all superimposed and very close. And interestingly, we also found in case (c) that the peaks of direct emission, lensing and photon ring are all distinguishable, which locate at the b c ∼ 6.6M, 6.4M and 6.8M respec- tively. Similarly, η and ψ 0 all decreased observed intensities of direct emission and rings as before. The right column of Fig. 7 shows the corresponding shadows and rings for different value of η and ψ 0 .

Shadows and photon spheres with spherical accretions
Besides the thin disk, there also may exist various kinds of spherical accretions around black hole. In view of this, we will continue to investigate the shadows of f (R) global monopole black hole surrounded by the static and infalling spherical accretions in this section.

The static spherical accretion
When the static accretion surrounded the f (R) global monopole black hole, the observed intensity for a distant observer can be expressed as And in current context, we have The symbol l, g, j (ν e ) and d prop represent the trajectory of emitted photons, the redshift factor, the emissivity per unit volume in the rest frame, and the infinitesimal proper length, respectively. And, ν o , ν e are the observed and emitted photon frequency, where ν o is assumed to be the frequency when the emitter radiate photons monochromatically. With the aid of Eqs. (16), (28), (31) and (32), the observed intensity for a static observer can be obtained as In the context of the static accretion, it is obvious that the obtained luminosity of the f (R) global monopole black hole is a function of r . Since the radius r connected with the impact parameter b c with in trajectory equation Eq. (17), we can plot the observed luminosity with respect to b c in Fig. 8 to intuitively show the shadows and photon spheres of the f (R) global monopole black hole. From Fig. 8, it is easy to see from the first row that, the observed intensity quickly increase to the maximum value firstly and then gradually decrease to a lower value with the increase of the impact parameter b c . Since the photon will revolve many times around the black hole, the maximum value of I obs (ν o ) occurs at b c = b p . When b c < b p , the

The infalling spherical accretion
Considering the fact that the accretion matters are always dynamical, it is also interesting to investigate the shadows and photon spheres when the f (R) global monopole black hole surrounded by the infalling spherical accretion. Similar to the static spherical accretion, the expression of the observed intensity for a distant observer Eq. (27) still holds true for the infalling spherical accretion. While in this case, the redshift factor should be rewritten as 3 Due to the existence of radiation field, there is always a small value of I obs (ν o ).
where, κ μ ≡ẋ μ is the four-velocities of the photon, which can be found in Eqs. (9)- (11). And, v μ obs ≡ (1, 0, 0, 0) represents the four-velocities of the static observer. While, the four-velocities of the infalling accretion is v μ emi , which is given by The proper distance in this case reads In addition, here the emissivity j (ν e ) in Eq. (29) remains unchanged. By using Eqs. (32)-(34), the observed intensity for a static observer in the context of the infalling spherical accretion can be obtained, which is In a similar way, the shadows and photon spheres of the f (R) global monopole black hole cast by the infalling accretion have been presented in Fig. 9. Similar to the case of the static accretion, the infalling accretion emission is also sharply peaked at b c = b p , and then decrease slowly with the increase of b c , which are shown in the first row of Fig. 9. The parameter η more or The specific intensity I obs (ν o ) seen by a distant observer for an infalling accretion less decreases the intensity of I obs (ν o ), while parameter ψ 0 increases it obviously. By comparing the second row of Fig. 9 with that of Fig. 8, it can be saw that the central region inside the photon sphere of the corresponding shadows for the infalling accretion are darker than that for the static accretion, which is a result of the Doppler effect. In addition, we can also find that the region outside the photon sphere is brighter for infalling accretion, which caused by the reason that I obs (ν o ) decreased slower than that of the static accretion as b c > b p .

Conclusions and discussions
In this paper, by considering three simple models of accretions, we have numerically investigated observational appearances of a f (R) global monopole black hole when the observer located at the asymptotic infinity and the certain region out of black hole. At first, we studied the effective potential and photon orbits of a f (R) global monopole black hole, and analysed the trajectories of light rays emitted from the north pole direction. Then, when emissions emitted from an optically and geometrically thin disk in the near region of black hole, the observational appearances of Direct emissions and rings have been carefully addressed with the aid of some typical toy-model functions. Finally, we also discussed shadows and photon spheres of black hole surrounded by the static and infalling spherical accretions. Throughout of this paper, the effect of the global monopole parameter η and f (R) parameter ψ 0 on shadows, rings and photon spheres of this black hole are clearly emphasized.
The result shows that the event horizon r h , the radius r p and impact parameter b p of photon sphere are all increases with increase of η and ψ 0 , while the cosmological horizon r c decreases. And, the effective potential decreases in the f (R) global monopole black hole by comparing with that of Schwarzschild black hole. For the thin-disk accretion, it turns out that the location of direct emission occurs at b c < 6.16715 and b c > 7.73130M for η = 0.07 and ψ 0 = 0.003, and the ranges of lensing ring are 6.16715M < b c < 6.42633M and 6.49637M < b c < 7.73130M, while photon ring is within the scope of 6.42633M < b c < 6.49637M. The related results for other two cases of η and ψ 0 have been presented in Eqs. (18) and (19). Those regions as well as the corresponding trajectories of light ray are clearly presented in Fig. 3. Based on the first three transfer functions, when the emission function is I emi (r ), one can see from Fig. 5 that the observed intensities of direct image, lensing ring and photon ring peaked at b c ∼ 9.6M, 6.9M and 6.4M, for subfigure (c) respectively. For the emission function I emi (r ), it shows that lensing ring and photon ring are indistinguishable for case (a) and (b) while distinguishable for case (c) in Fig. 6. Those rings are all superimposed on the direct emission at the region 6.4M ∼ 7.2M for case (c). For I emi (r ), as the emissions extended to the region of event horizon in Fig. 7, it shows in case (c) that the peaks of direct emission, lensing and photon ring are distinguishable, which locate at the bc ∼ 6.6M, 6.4Mand 6.8M, respectively. And, they are more or less indistinguishable for case (a) and (b). More importantly, we can conclude that the direct image is always similar to the emission profile, and contributes a dominant part of the total observed intensity. While, lensing ring and photon ring that around black hole shadow, always make a small contribution and a negligible contribution, for all three toy-model functions. For the static and infalling spherical accretions, we find that shadows and photon spheres cast by a f (R) global monopole black hole in the (x,y) plane are always similar. But, the brightness of shadow in the infalling accretion is darker than that in the static accretion, which caused by the Doppler effect. And, the brightness of photon sphere with infalling spherical accretion decreased slower than that of the static one. In addition, it shows that parameter ψ 0 increases the observed intensity for infalling accretion, although two parameters(η and ψ 0 ) all decrease the observed intensity in static case. In a word, observational appearances of a f (R) global monopole black hole surrounded by various accretions present some very different features so that it can be used for us to distinguish black holes from different gravity theories.
On the other hand, we in principle are also possible to constrain the f (R) parameter ψ 0 by using the shadow of M87 detected by the Event Horizon Telescope. Based on the discussions [55,104,105], the diameter of the shadow in units of mass for M87 has been obtained, which is d M87 ≡ D·δ/M ≈ 11±1.5, where D and δ are the angular size of the shadow and the distance to M87, respectively. Within 1δ and 2δ uncertainties, it is easy to see that the ranges of the diameter locate at the regions 9.5 ∼ 12.5 and 8 ∼ 14. So, this range may give rise to some upper or lower limits to black hole parameters. Combined with this fact, we can expect to set an upper limit on the f (R) parameter ψ 0 by using the shadow of a f (R) global monopole black hole. By computing the diameter of the resulting shadow, when the parameter η = 0.01, it is easy for us to find that the upper limits of the parameter ψ 0 are ψ 0 ≈ 0.034M for 1δ uncertainties and ψ 0 ≈ 0.051M for 2δ uncertainties. The detailed analysis of the observational constraints on black hole parameters will be carefully addressed in our next work. In a conclusion, it is true that black hole shadow can indeed present some observational constraints on the related parameters of black hole. In addition, there may also exists an optically thin but geometrically thick accretion around the black hole. Therefore, it is interesting for us to further discuss the observational appearances of black hole in the near future.
Data Availability Statement This manuscript has no associated data or the data will not be deposited. [Authors' comment: All the data are shown as the figures and formulae in this paper. No other associated movie or animation data.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .