Shadow, absorption and Hawking radiation of a Schwarzschild black hole surrounded by a cloud of strings in Rastall gravity

This paper studies the black hole shadow, absorption cross section, and Hawking radiation of a massless scalar field in the background of a static spherically symmetric black hole spacetime that is surrounded by a cloud of strings in Rastall gravity. Specifically, the effects of the parameters $a$ and $\beta$ on the photon sphere and shadow radii are investigated. The results show that as the negative parameter $\beta$ decreases, the photon sphere and shadow radii change in an N-shape. In addition, the absorption cross section obtained after solving the massless Klein-Gordon equation is calculated using the sinc approximation and the partial waves method. We compare the absorption cross section obtained by the sinc approximation and the partial waves method, and find it to be exceptionally consistent in the mid-to-high frequency region. Furthermore, the effects of parameters $a$ and $\beta$ on absorption are examined in detail. Finally, we study in detail the effects of the parameters $a$, $\beta$ and $l$ on the Hawking radiation power emission spectrum of the considered black hole. It turns out that the string parameter $a$ always suppresses the power emission spectrum, indicating that such black holes live longer when the string parameter $a$ is increased while other parameters are fixed.

or strong field conditions. Particularly, black holes, as one of the predictions, are arguably the most interesting and mysterious celestial bodies in our universe. The mystery of a black hole is that nothing, including light, can escape its event horizon. For the past few decades, the existence of black holes was only studied through indirect methods, until the first images of black holes appeared in 2019 [3]. This discovery provides many inspiring answers for our exploration of Einstein's theory of general relativity and for testing other revised theories of gravity, taking our understanding of black hole physics a major step forward. However, the basic theory proposed by Einstein cannot explain some phenomena or solve certain fundamental problems, e.g., the singularity problem and the conjecture that the covariant divergence of the energy-momentum tensor may be non-zero.
To account for the special case where the covariant divergence of the energy-momentum tensor does not vanish, Rastall [4] proposed a special modification of general relativity where the field equation is T µν ;µ = λR ,ν and λ = 0 corresponds to the Einstein equation. An important feature of Rastall's gravity is that the field equation T µν ;µ = λR ,ν is obtained directly by violating the normal conservation law, which does not rely on the metric or palatini formalism [5]. Also, it is important to note that Rastall's gravity appears to be consistent with experimental observations in the context of cosmology [6]. Specifically, the observational data include, but are not limited to, the age of the universe, helium nucleosynthesis, and Hubble parameters. What is more interesting is that the modified gravity gives us a lot of novel and interesting results at the cosmological level. Besides, some attention has been focused on a debate, namely, whether Rastall gravity is equivalent to Einstein gravity. Visser [7] thought that the modified gravity proposed by Rastall is a rearrangement of the matter sector of Einstein gravity. In other words, the geometrical part of the field equation is the same in both theories, so we just need to construct a new energy-momentum tensor to fulfill the ordinary conservation law. So the author claimed there is nothing new, such as gravity description, in Rastall proposal. Das et al. [8] had a conclusion that in the framework of non-equilibrium thermodynamics (for homogeneous and isotropic FLRW black hole model), generalized Rastall gravity theory is equivalent to Einstein gravity theory. However, other researchers disagree with Visser's ones, see for example the research of Darabi and his colleagues [9] who support that Rastall theory is not equivalent to Einstein gravity theory and give a simple example to prove that the claim proposed by Visser is incorrect. Moreover, they indicated that Rastall gravity theory is an "open" gravity theory in comparison to basic general relativity and has more compatible with observational cosmology. Hansraj et al. [10] also discussed this dispute and their results are consistent with Darabi et al. [9].
In this work, they showed that Rastall gravity can satisfy the fundamental conditions for physically viable model whereas Einstein gravity doesn't fulfill its requirements (see [10] for more detailed discussion). Some works [11][12][13][14] have shown the difference between Rastall gravity and Einstein gravity from theoretical or cosmological perspectives. Finally, regardless of whether the Rastall gravity is equivalent to Einstein gravity, Rastall gravity theory is worth studying or discussing because it faces a challenge from cosmological (astrophysical) observations.
String theory, on the other hand, holds that the fundamental unit of nature is not the point particle in particle physics, but an extended one-dimensional string. Letelier [15] proposed for the first time that the source of the gravitational field could be a cloud of strings, and gave an exact solution for a Schwarzschild black hole surrounded by a collection of strings in the context of Einstein's general relativity. In addition, black holes that treat a cloud of strings as the source of the gravitational field in the modified gravity have been studied [16][17][18][19][20]. For instance, Cai and Miao proposed a black hole solution in which a cloud of strings is the source of the gravitational field of a Schwarzschild black hole in the context of Rastall gravity [21]. The authors also analyzed fundamental thermal properties, quasinormal modes of gravitational perturbations, area spectra [22,23], and entropy spectra.
The experimental results reported by the Event Horizon Telescope Collaboration [3] not only directly prove the existence of black holes, but also allow us to directly observe the shadows of black holes. The theoretical analysis of black hole shadows has a long history.
For example, Synge [24] discussed the shadows of Schwarzschild spacetime, and Bardeen et al. [25] analyzed the shadows of Kerr black holes. In addition to shadow analysis performed in basic general relativity, it extends to various modified forms of gravity or arbitrarydimensional spacetime. Abbas and Sabiullah [26] studied the structure of timelike as well as null geodesics of regular Hayward black hole and found the massive particles, which move along timelike geodesics path, are dragged toward the black hole. To the best of our knowledge, numerous studies [27][28][29][30][31][32][33][34][35][36][37] have been devoted to studying the shadows of black holes under various modified gravity. More concretely, Gyulchev et al. [38][39][40] analyzed the shadows cast by different rotating traversable wormholes. Interestingly, the near horizon geometry is determined by the shadow cast by the black hole. Instead, the trajectory of light is affected by the plasma surrounding the black hole. This causes the geometric size and shape of the shadow on the Kerr spacetime to change [41]. In general, gravitational light deflection causes black hole shadows, and the trajectory of a photon in a vacuum depends on its impact parameter [42]. Therefore, we cannot ignore the role of the impact parameter in shadow formation.
Due to the special properties of black holes, we cannot directly study their internal structure. However, black hole is not an isolated system because it interacts with its surrounding environment, such as absorption, scattering and Hawking radiation. These interactions can convey information about the interior of the event horizon. In particular, as one of the interactions, the absorption cross section of black holes has received extensive attention from researchers. That's because one of the most useful and efficient ways to understand the properties of a black hole is to analyze the absorption of matter waves and the test field around the black hole. This series of studies began in the 1970s [43][44][45][46][47][48][49]. During that period, Sanchez found that the absorption cross section of Schwarzschild spacetime for scalar waves oscillates around the geometric capture cross section. About twenty years later, Das et al. [50] presented a key result that, in the low-energy regime, the absorption cross section of a coupled massless scalar field is equal to its event horizon area. Consequently, the literature on this particular topic has proliferated over the past few decades, covering various fields of research and several revision theories [51][52][53][54][55][56][57][58][59][60][61][62][63][64][65].
Furthermore, Hawking predicted that black holes are thermal systems, like black bodies, and then have associated temperature and entropy. Based on the analysis of quantum field dynamics in the context of curved space-time, Hawking pointed out that black holes emit radiation, known as Hawking radiation, from their event horizons [66,67]. Intriguingly, Hawking radiation depends on the type of particle and the geometry of the black hole. This is because the Hawking temperature T BH = f (r + ) 4π is one of the influencing factors. Moreover, Yale [68] has analyzed the Hawking radiation of particle scalars, fermions and bosons spin-1 using the tunneling method. In recent years, a large body of literature [69][70][71][72][73][74][75][76][77][78][79] has emerged on Hawking radiation on various modified gravity, including high-dimensional black holes. This paper investigates the black hole shadow, absorption cross section and Hawking radiation of the test scalar field of a Schwarzschild black hole surrounded by a cloud of strings in Rastall gravity. Specifically, Cai and Miao [21] presented the corresponding quasinormal modes of odd parity gravitational field by the WKB approximation. On this basis, our re-search contributes to further understanding of this black hole and its physical characteristics. This paper is organized as follows. The second section outlines the basic information of the black hole solution, that is, a Schwarzschild black hole surrounded by a cloud of strings in the context of Rastall gravity, and also gives the meaning of the influencing parameters.
The third part is devoted to the derivation of massless scalar equations and the analysis of related effective potentials. Section 4 analyzes the radius of the photon sphere and the shadow radius of the black hole in detail. Next, the absorption cross section of the scalar field is calculated using the sinc approximation and the partial wave method, and the effects of the parameters are also investigated. Section 6 gives the expression of Hawking radiation and the corresponding results for the Hawking radiation power emission spectra. The last section contains the summary and conclusions. Besides, we use the natural unit that c = G = = 1 in this paper.

II. THE SOLUTION OF A SCHWARZSCHILD BLACK HOLE SURROUNDED BY A CLOUD OF STRINGS IN RASTALL GRAVITY
The field equations of the Rastall gravity [4] are as follows, where κ and λ represent the Rastall gravitational coupling constant and the Rastall parameter, respectively. Moreover, β is defined as the product of these two parameters, i.e., β ≡ κλ. From the above equations we have that where R, T denote the Ricci scalar and the trace of the energy-momentum tensor, respectively. Besides, κ = 4β−1 6β−1 8π under the Newtonian limit [80]. It can be seen from the above equations that, the Einstein gravity is recovered and the energy-momentum tensor is conserved when the Rastall gravity parameter λ vanishes, i.e., β = 0.
We consider the case where the metric is static and spherically symmetric, with the metric [21] f (r) It is worth noting that the Rastall theory should satisfy the Newtonian limit [80]. Therefore, the cases β = 1 6 and β = 1 4 are not allowed. The parameter a needs to satisfy a specific constraint, namely a ≡ κb where b is a constant of integration associated with a cloud of strings. Specifically, β and a represent the influence of the Rastall gravity and the string, respectively. Consequently, the Rastall gravity is converted to Einstein gravity when β = 0.
Meanwhile, when a equals to 0, the Schwarzschild spacetime is restored.

III. SCALAR WAVE EQUATION
The massless scalar field Ψ governed by the massless Klein-Gordon equation in curved spacetime can be formulated as and then the massless scalar field Ψ can be decomposed as follows where P l (cos θ) denotes the Legendre polynomial, l and m represent the corresponding angular quantum number and magnetic quantum number, respectively. In addition, the function Ψ ωl satisfies the following ordinary differential equation, where V ef f (r) stands for the corresponding effective potential that is defined as Moreover, by substituting the metric in the effective potential, the specific potential is reformulated as Additionally, we define the following tortoise coordinate change Consequently, the equation (9) is equivalent to Note that both the metric f (r) and the effective potential V ef f (r) are divergent when β is set to -0.5. Besides, to satisfy the condition that the effective potential V (r) → 0 when r → ∞, we have β < 1 6 . Hence, the domain of β should be in (−0.5, 1 6 ). Moreover, due to the condition a ≡ κb, the domain of a strictly relies on the positivity (negativity) of the parameter β. For β < 0, the barrier of the effective potential V ef f (r) disappears as the parameter a approaches to 1. Accordingly, the domain of the parameter a is set to [0, 1). In contrast, for β > 0, when the parameter a is larger, a black hole surrounded by a cloud of strings in Rastall gravity has no event horizon. For instance, when β = 0.1, the domain of a is set to [0, 0.3]. Fig.1 shows the behaviour of the effective potential V ef f (r) with respect to r for different angular quantum numbers l when a = 0.1, β = 1 10 . We find that the peak value of the effective potential increases when the angular quantum number l is increased. Furthermore, the potential V ef f (r) first increases, then decreases, and finally tends to zero at r → ∞.
As shown in Fig.2, to compare the effects of parameters a and β on the effective potential V ef f (r), we depict the behaviour of V ef f (r) with respect to a and β when β < 0 and β > 0, respectively. Specifically, for β > 0, i.e., when the parameter β is fixed to 1 10 , the barrier height of the effective potential decreases as the string parameter a increases. It is clear that the peak of the effective potential becomes smaller and shifts to the right side as a increases.
Next, we vary the Rastall parameter β and fix a to 0.1. It can be seen that with the increase of β, the peak value of the effective potential decreases, and the position of the peak value does not change much compared with the case where the string parameter a changes.
Meanwhile, when β < 0, one can see that for the same value of a, the barrier height of the potential first increases and then decreases with decreasing β. Also, the peak position firstly shifts to the left and then to the right. Furthermore, when the parameter a is varied, at the same value of the Rastall parameter β, the barrier height decreases and the peak position shifts to the right as the parameter a increases. 10 . The plot in the upper right corner is that the parameter β is variable when a = 0.1. For β < 0, the picture in the lower left corner is that the parameter β is variable with a = 0.6. the figure in the lower right corner is that the parameter a is variable for fixed β = − 1 3 .
We add some boundary conditions for the Schörding-like equation (13) because we are interested in the absorption cross section and Hawking radiation. Near the horizon regime and at infinity, one can find that ψ ωl (r * ) need to satisfy the following boundary conditions where R ωl and T ωl in the action denote reflection and transmission coefficients, respectively.
Due to the conservation of flux, R ωl and T ωl satisfy the following constraint Furthermore, the phase shift δ l can be defined as Next, we will discuss the black hole shadows, absorption cross section and Hawking radiation based on the last two sections.

IV. SHADOWS
In this section, we investigate the role of the Rastall parameter β and the string parameter a on the shadow radius of a black hole enclosed by a cloud of strings in Rastall gravity.
The photon trajectories of a black hole surrounded by a cloud of strings in Rastall gravity can be represented by null geodesics [26,81]. The Lagrangian of geodesic equations for the curve spacetime have the following form where the overdot symbol denotes the differentiation with respect to the affine parameter τ . Without loss of generality, we consider an analysis restricted to the equatorial plane, i.e., θ = π 2 . By using the Euler-Lagrange equation, the t and φ coordinates are expressed as, where E, L are motion constants, representing the energy and angular momentum of the massless test particle, respectively.
Hence, by substituting Eq. (18) and Eq. (19) in the Lagrangian equation (17), the Lagrangian expression can be written aṡ furthermore, we define where V stands for the effective potential of the massless test particle. Besides, the nulllike geodesics of the equatorial circular motion in static spherically symmetric spacetime should satisfy the conditionsṙ = 0 andr = 0. Consequently, we have V = E 2 and dV dr = 0, indicating the stability of circular null geodesics. The equations V (r p ) = 0 and V (r) |r=rp = 0 [82] represent the circular orbit of the photon, that is, the photon sphere radius r p .
Moreover, the critical impact parameter b c can be expressed as On the other hand, the black hole shadow radius r s is represented by the celestial coordinates (x, y) as follows Specifically, the effects of the parameters a and β on the photon sphere and shadow radii are shown in Table I. For β > 0, from Table I one can see that for fixed a = 0.1 (β = 1 10 ), the photon sphere and shadow radii increase as the parameter β (a) increases. Furthermore, as the string parameter a tends to 0.3 when β = 1 10 , the black hole shadow radius increases rapidly. For β < 0, when we set a = 0.3, we observe that the photon sphere and shadow radii first increase, then decrease and finally increases as the parameter β decreases. A possible reason is that the metric f (r) is not a monotonic function of the Rastall parameter β in the range −0.5 < β < 0. Therefore, when the parameter β is set to − 1 3 , the photon sphere and shadow radii increase as a approaches to its parameter maximum. In this section, we calculate the absorption cross section using two methods, viz., the sinc approximation method and the partial waves method where the gray-body factor is calculated by sixth-order WKB method. Besides, we view the capture cross section as a reference.
It is known that the absorption cross section at the low-frequency and high-frequency limits can be calculated by different analytical approximations. The total absorption cross section of massless scalar waves in an arbitrary-dimensional general spherically symmetric black hole inclines to its area [50] in the low-frequency regime, which is the event horizon of the black hole. In the high-frequency regime, the total absorption cross section of the massless scalar field converges to the geometric capture cross section, described by the following null where b c denotes the above critical impact parameter.

A. sinc approximation
Sanchez [49] proposed that in the high-frequency regime, the total absorption cross-section oscillates near the above-mentioned capture cross section (27/4)πr 2 s , where r 2 s = 2M , and has an interval of oscillation peaks, ∆ = 2 √ 27M . In addition, Sanchez also presented the following analytical approximation of the absorption cross section which has the best fit when A = 1.14 ∼ √ 2 and B < 10 −4 .
Furthermore, the Sanchez approximation was generalized by Décanini et al. to static spherically symmetric spacetimes of arbitrary dimensions. Décanini et al. [83] showed that in the eikonal state, the fluctuation of the absorption cross section was completely and very simply described by the properties of the null unstable geodesics located on the photon sphere. Important characteristics are the orbital period and the Lyapunov exponent.
Specifically, the sinc approximation of the absorption cross section in a d-dimensional static and spherically symmetric black hole is given by where the oscillation part of the absorption, i.e., σ osc abs , is expressed as with sinc(x) denoting the sine cardinal and d representing the dimension of the black hole. Besides 2π rc √ f (rc) = 2πb c indicates the orbital period of the black hole on the photon sphere [84]. The parameter η c for measuring the instability of the circular orbit on the photon sphere is defined as for instance, the sinc approximation of the absorption cross section of a Schwarzschild black hole at the high-frequency limit is written as

B. Partial wave approach
We consider that the field Φ, which is purely ingoing waves at the event horizon, is the sum of the monochromatic incident plane wave Φ I and outgoing scattered wave Φ S in the far-field, that is, Without loss of generality, we assume that the direction of wave propagation is along the z-axis. Accordingly, the monochromatic incident plane wave Φ I and the outgoing scattered wave Φ S are respectively defined as wheref (θ) denotes the scattering amplitude. Moreover, e iωz can be decomposed as [85] e iωz = ∞ l=0 (2l + 1)i l j l (ωr)P l (cos θ), with j l (.) representing the spherical Bessel function.
Hence, Eq. (33) in the far-field can be rewritten as follows, C ωl e −iωr + e −iπ(l+1)e iωr P l (cos θ), where C ωl is given by The field solution Φ depends on the boundary conditions (14). This means that the ingoing part of Φ should match the incident plane wave Φ I . Therefore we obtain The absorption cross section depends on the flux of particles that enter the black hole through the effective potential. Hence, we can introduce the four-current density vector as follows and the above equation satisfies the conservation law, that is By substituting Eq.(38) into Eq.(39) under the boundary condition Eq. (14), we obtain the four-current density vector by surface integral as where N (r) is the flux that passes the surface Σ with a constant radius r and dΩ = sin θdθdϕ.
The flux is a constant, and when we consider the stationary scenarios, N (minus) represents the particles passing through the potential and entering the black hole [86]. Besides, we have used the orthogonality of Legendre polynomials, i.e., where x = cos θ. Furthermore, the absorption cross section σ abs is defined as the ratio of the particle flux |N | to the plane wave incident current ω. Hence, the absorption cross section can be written as and the partial absorption cross section can be expressed as In order to study the effects of Rastall and string parameters on the absorption cross section of the scalar field, we need to calculate the phase shift δ l , that is, the transmission coefficient. In this paper, we use the WKB approximation to obtain the transmission coefficient T ω . Assuming that the probability of the incident plane wave is equal to 1, Eq. (15) can be expressed as The transmission probability of different multipole numbers l can be obtained with the help of the sixth-order WKB method, with where α is obtained by In Eq.(48), V 0 represents the maximum value of the potential at r = r 0 , and the prime denotes the derivative of the potential at r = r 0 with respect to r * . Moreover, Λ i (K) indicates a higher-order correction of the WKB method, which depends on K and the 2i order derivative of the potential at its maximum position [87,88].
Specifically, we express the third-order method as follows, In Eqs. (49), the superscripts (2,3,4,5,6) of the effective potential represent the corresponding differentials with respect to the tortoise coordinate r * , and b = n + 1 2 . Besides, since the specific expressions of Λ 4 (K), Λ 5 (K) and Λ 6 (K) are overly cumbersome (see Ref. [88]), they will not be described in detail here. In addition, during the calculation, we find that when the Rastall parameter β is set as a fraction, the results and figures of the WKB approximation calculation are more accurate than when β is set as a decimal [21]. This phenomenon can be attributed to the term r value, and finally decreases to almost the same value with increasing ω. Furthermore, it is easy to see that as the string parameter a increases, the partial absorption increases and its peak position shifts to the left. When we fix the string parameter and change the Rastall parameter, one can get that the peak value of the partial absorption cross section increases as the Ratall parameter β increases.
In Fig. 4 we present the total absorption cross section of a Schwarzschild black hole surrounded by a cloud of strings in Rastall gravity for different values of the string parameter, where l goes from 0 to 10 and β = 1 10 . Specifically, the horizontal solid line represents the geometric capture section. As shown in Fig. 4, the dashed curve is the sinc approximation result, and the solid curve is the partial wave result using the sixth-order WKB approximation. We show that increasing the parameter a results in incrementing the absorption cross section. We also notice that the two curves are significantly different at small values of frequency. Moreover, as the string parameter increases, the difference is more pronounced and the range of oscillation amplitudes is significantly wider. However, in the high frequency regime, the total absorption cross sections obtained by these two methods are in good agreement and converge to the geometric capture cross section. In Fig. 5, our results show that when we fix the value of a and increase β, the absorption cross section increases. The difference between the two curves also increases significantly in the low frequency regime due to the Rastall parameter.
As shown in Fig. 6, we describe the behavior of the partial absorption cross section for different parameters obtained by the sixth-order WKB method when β < 0. From the left figure, where the Rastall parameter is treated as a variable and the parameter a is fixed, we observe that the partial absorption cross section does not monotonically increase as the Rastall parameter β decreases. The partial cross sections intersect in the range of 0.2 < ω < 0.3 due to the effective potential. Therefore, the variation trend of the partial cross section is 'N' type. We also observe in the right plot that when we increase the string parameter, the partial cross section increases monotonically. Moreover, the peak position of the partial cross section is evidently shifted to the left. Finally, we observe that as the multipole number l increases, the partial cross section decreases and its peak position shifts to the right.
In Fig. 7 we give the total absorption cross sections of the massless scalar field by varying the Rastall parameter β (β < 0) and fixing the string parameter a = 0.6. We can observe that the change of the total absorption cross section as a function of β is similar to that of the partial absorption cross section. This is because the higher the potential barrier, the more particles are scattered back to the black hole by the potential barrier. In addition, we can see that when we reduce the Rastall parameter to −0.5, the difference between the solid curve and the dashed curve gradually decreases.
In Fig. 8 we present the total absorption cross section of the massless scalar field when β < 0, changing the string parameters and fixing β = − 1 3 . It can be observed that the difference between the two curves is the smallest at the low-frequency limit compared to the above three cases. Furthermore, when ω is large, the total absorption cross section as a function of ω goes into the capture cross section. We also notice that the total absorption cross section, as well as the oscillation amplitude, increases with increasing string parameters.

VI. HAWKING RADIATION
In this section, we employ the sixth-order WKB method to calculate the Hawking radiation for massless scalar fields. Furthermore, we analyze the effects of the string and the Rastall parameters on Hawking radiation in the background of a Schwarzschild black hole surrounded by a cloud of strings in Rastall gravity.
A black hole behaves almost in the same way as a black body, emitting particles when its temperature is proportional to the surface gravity [66]. Hawking also presented that black holes can radiate particles in the form of thermal. This is due to the quantum tunneling effect created by the vacuum fluctuations near the event horizon of the black hole. Therefore, if we consider quantum effects and the laws of thermodynamics are satisfied, black holes can produce radiation. This phenomenon is known as Hawking radiation. The Hawking radiation calculated by the gray-body factor has the following expression [89,90] where N l is the multiplicity that depends only on the black hole dimension. Moreover, for the massless scalar field in a four-dimensional black hole, l and N l satisfy the condition N l = 2l + 1. T ωl denotes the above gray-body factor and T BH represents the Hawking temperature. Specifically, the Hawking temperature of static spherically symmetric spacetime can be written as where f (r h ) = 0 and r h is the radius of the event horizon. Besides, the string and Rastall parameters need to satisfy the previous parameter range, i.e., −0.5 < β < 1 6 and 0 ≤ a < 1. By substituting Eq. (52) and N l = 2l + 1 into Eq.(50), we can further obtain the Hawking power emission spectrum Fig.9 compares the effects of parameters a and β on the Hawking power emission spectrum of the massless scalar wave when β is non-negative. We can clearly observe in the upper panel that for a given l and β, increasing the parameter a depresses the power emission spectrum. Moreover, the peak power emission spectrum gradually shifts to low frequencies as a increases. It is clear from the middle panel that when we fix l and a, but increase the parameter β, the peak power emission spectrum gradually decreases and moves to low frequencies. As the multipole number l increases, we can get from the lower panel that for a massless scalar field, the power emission spectrum decreases and the peak position shifts towards high frequencies. In conclusion, the parameters a, β and l suppress the power emission spectrum. Besides, it is easy to see that if the values of parameters a and b are chosen larger, the lifespan of the black hole will be longer.
This trait is more easily observed in Fig.10, which plots the effects of parameters a, β and l on the power emission rate (as a function of ω) for the scalar wave in the range β < 0. From the upper figure we can see that when we increase the parameter a, the power emission spectrum decreases. That is, under the condition that β is constant, the increase of the string parameter a leads to a decrease in the energy emission rate, thus making the lifetime of the black hole longer. Furthermore, we also observe in the center panel that with decreasing Rastall parameter, for fixed l and a, the peak value of the power emission rate increases and then decreases, and the peak position first shifts to high frequency and then moves to low frequency. Finally, we fix the two parameters a = 0.6 and β = − 1 3 and analyze the effects of the multipole number l in the lower panel. It is clear that a larger multipole number results in a lower power emission spectrum. Besides, it is worth noting that the low multipole number l dominates the energy emission rate, while the contribution of the high multipole number l is extremely small and thus negligible.

VII. CONCLUSION AND DISCUSSION
In the previous sections, we have comprehensively studied the black hole shadow, absorption cross section and power emission spectrum of Hawking radiation for the massless scalar shadow and photon sphere radii caused by the curved light ray. Because we consider the black hole to be static spherically symmetric, the radii of the photon sphere and shadow are constant. In other words, the black hole shadow has spherical symmetry. Besides, the radius of the photon sphere increases as the parameter a increases. However, when we consider β as a variable, the photon sphere and shadow radii fluctuate abnormally. The reason is that when the Rastall parameter is less than zero, the metric f (r) changes abnormally. Second, with the help of the sixth-order WKB method, we calculated the absorption cross section of the scalar field in detail. To compare the accuracy of the sixth-order WKB, we also presented the results of the sinc approximation with the geometric capture cross section as a reference. From Figs.3, 4 and 5, we can clearly observe that larger values of the parameters a and β enhance the partial or total absorption cross section when β > 0. However, in the low frequency range, when a or β is set to a larger value, the results calculated by the two methods are quite different. Furthermore, in Figs.6, 7 and 8, we plotted the partial and total absorption cross sections when β<0. Unlike the case where β is positive, the absorption cross section does not always grow as the Rastall parameter decreases. Since the potential barrier reflects waves, the change in the absorption cross section is exactly the opposite of the change in the potential barrier. Hence, as β decreases, the total absorption cross section first increases, then decreases and finally increases again. It is worth mentioning that the smaller the value of β, the smaller the difference between the two approximations. Very importantly, in the mid-high frequency region, the total absorption cross section and the sinc approximation are in good agreement and in all cases oscillate around the geometric capture cross section σ geo .
Finally, we investigated the energy emission rate of Hawking radiation. Specifically, the power emission rate is affected by the string parameter, the Rastall parameter as well as the multipole number. In Fig.9, we found that both a and β suppress the power emission spectrum, and the peak position shifts to a lower energy region. Moreover, the multipole number l also significantly depresses the power emission spectra whereas the peak position shifts to the higher frequency regime. The case of β < 0 is also similar to the case of β > 0 above, except the case where β varies and a is fixed. As the Rastall parameter decreases, the power emission spectrum first increases and then decreases, at the same time, the peak position first moves to the higher frequency region and then enters the lower energy region.