Surrounded Bonnor–Vaidya solution by cosmological fields

In the present work, we generalize our previous work (Heydarzade and Darabi in arXiv:1710.04485, 2018 on the surrounded Vaidya solution by cosmological fields to the case of Bonnor–Vaidya charged solution. In this regard, we construct a solution for the classical description of the evaporating-accreting charged Bonnor–Vaidya black holes in the generic dynamical backgrounds. We address some interesting features of these solutions and classify them according to their behaviors under imposing the positive energy condition. Also, we analyze the timelike geodesics associated with the obtained solutions and show that some new correction terms arise in comparison to the case of standard Schwarzschild black hole. Then, we explore all these features for each of the cosmological backgrounds of dust, radiation, quintessence and cosmological constant-like fields in more detail.


Introduction
In 1951, Vaidya introduced a new non-static solution, describing a spherical symmetric object possessing an outgoing null radiation, for the Einstein field equations [1,2]. This solution is characterized by a dynamical mass function, depending on the retarded time coordinate. Based on its dynamical nature, the Vaidya solution has been used for studying the process of spherical symmetric gravitational collapse and as a testing ground for the cosmic censorship conjecture [3][4][5][6][7], and as a dynamical generalisation of the Schwarzschild solution representing a spherically symmetric evaporating black hole, as well as studying the Hawking radiation [8][9][10][11][12][13][14][15]. This solution was generalized by Bonnor and Vaidya to the charged case, well known as the Bonnor-Vaidya solution [16]. This solution and its interesting features and applications are studied in [17][18][19][20] as instances. a e-mail: yheydarzade@bilkent.edu.tr b e-mail: f.darabi@azaruniv.edu Further generalization of the original Vaidya solution were introduced in [21] by Husain for a null fluid with a particular equation of state, and in [22] by wang and Wu using the fact that any linear superposition of particular solutions is also a solution to the Einstein field equations. Using this approach, one can find other general solutions such as the Vaidya-de Sitter [23], Bonnor-Vaidya-de Sitter [18,[24][25][26][27] and radiating dyon solutions [28]. The Vaidya solution and its generalizations are also studied in the context of modified theories of gravity, see for examples [4,[29][30][31][32][33].
Black holes have such an strong gravitational attraction that their nearby matter, even light, cannot escape from their gravitational field. Although, the black holes cannot be observed directly but there are some different ways to detect them in the binary systems as well as at the centers of their host galaxies. The most promising way for this detection is the accretion process. In the language of astrophysics, the accretion is defined as the inward flow of captured matter fields by a gravitating object towards its centre which leads to an increase of the mass and angular momentum of the accreting body. The observation of supermassive black holes at the center of galaxies represents that such massive black holes could have been gradually developed through the appropriate accretion processes. However, the accretion processes do not always increase the mass of the accreting bodies but they can also decrease their mass and lead them to shrink. It is shown that the accretion of phantom energy can decrease the black hole area [34][35][36][37][38][39]. For instance, in [34], it is shown that black holes will gradually vanish as the universe approaches to a cosmological big rip state. The shrink of the black hole area during the accretion of a potentially surrounding field is an interesting phenomena in the sense that it can be considered as an alternative for the black hole evaporation through the Hawking radiation or even as an auxiliary for speeding up the evaporation process. One physical explanation for diminishing the black hole mass through the accretion process is that the accreting particles of a phantom scalar field have a total negative energy [40]. Similar particles with negative energies are created through the Hawking radiation process as well as in the process of energy extraction from a black hole by the Penrose mechanism. Thus, the accretion process into the black holes is one of the most interesting research fields in relativistic astrophysics to answer how black holes affect their cosmological surrounding fields and what are the consequences or what are the influences of these surrounding fields on the features, dynamical behaviors and abundance of black holes [41][42][43][44][45][46][47][48][49]. See also [50] for the accretion of dark energy into black holes, and [38,39,51,52] for the accretion into the charged black holes.
In the present work, following the approach of [53,54] and [55,56], we construct a dynamical solution for the classical description of the evaporating-accreting Bonnor-Vaidya black holes in generic dynamical backgrounds. The organization of the paper is as follows. In Sect. 2, the surrounded Bonnor-Vaidya black hole solution and some of its general features are introduced. In Sects. 2.1-2.4, the special classes of this solution named as the Bonnor-Vaidya black hole surrounded by the dust, radiation, quintessence and cosmological constant fields, as well as their properties are studied in detail. Finally, the Sect. 3 is devoted to the summary and concluding remarks.

Surrounded evaporating-accreting Bonnor-Vaidya black hole solution
In this section, we generalize our previous solution [55,56] to the surrounded charged Bonnor-Vaidya black hole solution by following the approach of [53,54]. There are two main motivations for us for doing this generalization. The first one is that the existence of the charge can drastically change the global structure of the original spacetime [57]. For instance, we know the Reissner-Nordström black hole has a very distinct causal structure relative to the Schwarzschild case such that it predicts infinite series of parallel universes. The second reason is that a charged black hole possesses a spacetime structure almost similar to a rotating one, the Kerr black hole.
Regarding that the existing spherical symmetry in the charged case makes it more easily analyzable, then understanding the structure of a charged black hole may be a suitable ground to better understanding the structure of a more realistic rotating one. We consider the general spherical symmetric spacetime metric where d 2 = dθ 2 + sin 2 θ dφ 2 is the metric of two dimensional unit sphere and f (u, r ) is a generic metric function depending on both of the the radial coordinate r and advanced/retarded time coordinate u. The cases = ± 1 associated with the possible outgoing-ingoing flows corresponding to the effectively evaporating-accreting Bonnor-Vaidya black hole. For the metric (1), the nonvanishing components of the Einstein tensor are given by where dot and prime signs denote the derivatives with respect to the time coordinate u and the radial coordinate r , respectively. Thus, one can find that the total energy-momentum tensor supporting this spacetime must have the following non-diagonal form which must possess the same symmetries in the Einstein tensor G μ ν . Then, regarding the equations in (2), the equalities G 0 0 = G 1 1 and G 2 2 = G 3 3 in the Einstein tensor components demand the equalities T 0 0 = T 1 1 and T 2 2 = T 3 3 for the energy-momentum tensor components, respectively. Then, one may introduce an energy-momentum tensor obeying these properties as in our previous work giving the surrounded Vaidya black hole [55,56]. One possible generalization to [55,56] can be obtained by including the Maxwell electromagnetic energy-momentum tensor. In the following, we prove that the resulting total energy-momentum tensor obeys all the symmetries in G μ ν . Then, we show that this provides the possibility of finding the charged Bonnor-Vaidya black hole solutions [16] in a general dynamical background in the context of the Einstein-Maxwell theory. Thus, we consider the Einstein field equations, corresponding to the components of the Einstein tensor (2), with the total energymomentum tensor T μ ν given by where τ μ ν is the energy-momentum tensor associated to the Bonnor-Vaidya null radiation-accretion as such that σ = σ (u, r ) is the density of the "outgoing radiation-infalling accretion" flow and k μ = δ 0 μ is a null vector field and E μ ν is the trace-free Maxwell tensor given by where F μν is the antisymmetric Faraday tensor satisfying the vacuum Maxwell equations The spherical symmetry in the spacetime metric (1) dictates the only non-zero components of F μν tensor to be F 01 = −F 10 . Then, from the Eq. (7), one obtains where Q(u) is the dynamical electric charge and its associated null current is whereQ(u) = d Q(u) du . Using the Eqs. (1), (6) and (8), the only non-vanishing components of Maxwell tensor E μ ν will be Finally, T μ ν in (4) is the energy-momentum tensor of the surrounding perfect fluid defined as in [53] T 0 0 = − ρ s (u, r ), Here, the subscript "s stands for the surrounding field which generally can be a dust, radiation, quintessence and cosmological constant or even any complex field constructed by the combination of these fields. From (11), it is seen that the spatial profile of the surrounding energy-momentum tensor is proportional to its time component, representing the dynamical energy density ρ s (u, r ), with the arbitrary parameters α and β which depend on the internal structure of the corresponding surrounding fields. The isotropic averaging over the angles results in [53] T The last equality follows from the fact that r i r j = 1 3 δ i j r n r n which results in the barotropic equation of state for the surrounding field as where p s (u, r ) and ω s are the dynamical pressure and the constant equation of state parameter, respectively. Then, regarding the Einstein tensor components in (2) and the total energy-momentum tensor given by the Eqs. (3)- (5) and (11), we find that T 0 0 = T 1 1 and T 2 2 = T 3 3 . These exactly provide us the principle of additivity and linearity condition proposed in [53] for determining the free β parameter in the energy momentum-tensor (11) as Now, by substituting α and β parameters given in (13) and (14) into (11), one obtains the non-vanishing components of the surrounding energy-momentum T μ ν in the following forms Then, having the Einstein tensor components (2) and the corresponding general energy-momentum tensor T μ ν in (4), we have the corresponding field equations. The G 0 0 = T 0 0 and G 1 1 = T 1 1 components of the Einstein-Maxwell field equations give 1 Similarly, from and G 2 2 = T 2 2 and G 3 3 = T 3 3 components lead to By simultaneous solving the differential equations (16) and (18), one can find the following general solution for the metric function with the energy density of the surrounding filed in the form of in which M(u), Q(u) and N s (u) are integration coefficients representing the black hole dynamical mass and charge and dynamical surrounding field structure parameter, respectively. The weak energy condition on the energy density (20) of the surrounding field, i.e ρ s ≥ 0, requires implying that for the surrounding fields with ω s ≥ 0, it is needed to have N s (u) ≤ 0 and conversely for ω s ≤ 0 we have N s (u) ≥ 0. Regarding (19), the metric (1) takes the form of representing an effectively evaporating-accreting charged black hole in a dynamical background.
Here, it is worth to discuss about the stability of this black hole. The stability is achieved if the metric solution (22) be time independent, namely ∂ u g ab = 0 or However, because of different powers of r which yields an r -dependent differential equation, one cannot obtain a global stability condition. In other words, this metric solution cannot be stabilized unless in a local way. It seems this is the case for any other metric solution, studied throughout this paper, for which the corresponding differential equation of stability condition is r -dependent. Regarding (22), one may realize the following two distinct subclasses for this general solution for the field equations (16) and (18).
• The solution by setting f = f (u, r ) and ρ s = ρ s (r ) These considerations lead to M = M(u), Q = Q(u) and N s = constant in the metric function f (u, r ) and σ (u, r ) = 0 for the energy density. Then, there is no dynamics in the surrounding field and consequently the accretion of the surrounding field by the black hole cannot happen. Indeed, this case represents an evaporating charged black hole solution in a static background. The radiating charged black holes in an empty background (ρ s = 0) known as the original Bonnor-Vaidya solution [16], and in (anti)-de Sitter space (ρ s = ρ = constant) are special subclasses of this solution [18,58]. Some interesting features of these black holes can be found in [15,19,20,59]. • The solution by setting f = f (r ) and ρ s = ρ s (r ) These considerations lead to M = constant, Q = constant and N s = constant in the metric function f (u, r ) and consequently σ (u, r ) = 0 for the radiationaccretion density. This case represents a static charged back hole in a static background and consequently, there is no radiation-accretion. The Reissner-Nordström black hole as well as its generalization to (anti)-de Sitter background are special subclasses of this solution. For a general background, not just the (anti)-de Sitter background, it is interesting that for a constant mass and charge black hole in a static non-empty background, using the coordinate transformation one arrives at the solution of the Reissner-Nordström black hole surrounded by a surrounding field as This solution is a generalization of the Kiselev solution [53] to the charged case and its interesting properties are studied in [60][61][62][63]. Then, the generalized Kiselev solution is a subclass of our general dynamical solution (22) in the stationary limit.
Substituting f (u, r ) given by (19) in (17) gives the radiationaccretion density of the effectively evaporating-accreting Bonnor-Vaidya black hole as Then, we observe that the radiation-accretion density is resulted not only from the change in the black hole mass (σ M ) and surrounding field (σ N s ) but also from the change in the charge of the black hole (σ Q ), representing the electromagnetic energy. In this case, the black hole may have just the outgoing charged null radiation. Turning off the surrounding field dynamics, i.eṄ s (u) = 0, we recover the energy flux associated to the mass and charge changes of the central black hole corresponding to the Bonnor-Vaidya solution [16]. It is seen that ifṀ(u), Q(u)Q(u) andṄ s (u) have a same order of magnitude, the following distinct physical situations can be realized.
• For ω s < 0, the charge contribution is dominant near the black hole. For the far distances (r ), the black hole charge contribution falls down faster than the black hole mass and the surrounding field contributions, respectively, i.e |σ Q | < |σ M | < |σ N s |. Then, at large distances the surrounding field contribution is dominant.
• For 0 < ω s < 1/3, the charge contribution is dominant near the black hole. For the far distances (r ), the charge contributions falls down faster than the surrounding field and mass contributions, respectively, i.e |σ Q | < |σ N s | < |σ M |. Then, at large distances the black hole mass contribution is dominant.
• For ω s > 1/3, the surrounding field contribution is dominant near the black hole. For the far distances (r ), the the surrounding field contributions falls down faster than the charge and mass contributions, respectively, i.e |σ N s | < |σ Q | < |σ M |. Then, at large distances the black hole mass contribution is dominant again.
Considering the positive energy density condition on the total radiation-accretion density σ (u, r ) requires This inequality confines the dynamical behaviours of the charged Bonnor-Vaidya black hole and its background at arbitrary time and distance (u, r ). In the case of a static background and neutral black hole, as in the Vaidya's original solution, it is required that andṀ(u) have the same signs.
In the presence of the black hole charge and background field dynamics, it is not mandatory that andṀ(u) take the same signs, and the satisfaction of the positive energy density condition can be achieved even by their opposite signs depending on the dynamics of the black hole charge (Q(u)) and surrounding field parameters (Ṅ s (u) and ω s ). Then, the dynamical behaviour of the surrounding field is governed by Then, at an arbitrary distance r from the black hole, the surrounding field must obey the above conditions. Interestingly, for the special case ofṄ s (u) = 2r 3ω s −1 (Q(u)Q(u) − rṀ(u)), there is no pure radiation-accretion density, i.e σ (u, r ) = 0. This case is associated with two possible physical situations. The first one corresponds to the situation where for any particular distance r 0 , the background fielḋ N s (u) and black hole withṀ(u) andQ(u) behave such that their contributions cancel out each others, leading to σ (u, r 0 ) = 0. The second situation corresponds to the case where for the given dynamical behaviors of the black hole and its background, one can find the particular distance r * (u) = Roots o f [2Ṁ(u)r 3ω s − 2Q(u)Q(u)r 3ω s −1 +Ṅ s (u)], which is generally dynamical, possessing zero energy density σ (u, r * (u)). Generally, to have a particular distance at which the density σ (u, r * ) is zero, the reality and positivity of r * also requires thatṀ(u),Q(u) andṄ s (u) obey some specific conditions. Here, due to the fact that finding the location of r * (u) in its general form is very complicated, in comparison to our previous solution [55,56], we discuss in this regard by considering some specific surrounding fields through the following subsections. However, before we study the features of obtained solution for some specific cosmological surrounding fields, we would like to investigate the timelike geodesics corresponding to the metric (22) in its general form. Due to the spherical symmetry, the geodesics for this metric lie on a plane, in which one may choose θ = π/2 for the sake of simplicity. Considering the action where the star sign denotes the derivative with respect to the proper time τ , and using the variation, we arrive at the following equations for the ϕ, r and u variables, respectively, and where L is the conserved angular momentum per unit mass and dot and prime signs denote the derivative with respect to u and r coordinates, respectively. Substituting (30) in (31), we have Moreover, using the timelike geodesics condition, i.e g μνẋ μẋ ν = −1, we obtain where we have used also the Eq. (30). Now, substituting (33) and (34) in (32), we arrive at the following general equation of motion * * for the radial coordinate r . Substituting our metric function (19), our equation of motion (35) takes the form of * * Consequently, we realize the following interesting points.
• The terms in the first line are exactly the same as that of the standard Schwarzschild black hole solution except the time dependance in the mass of the black hole. Here, the terms represent the Newtonian gravitational force, the repulsive centrifugal force and the relativistic correction of general relativity (which accounts for the perihelion advance of planets), respectively. • The terms in the second line are new correction terms, in comparison to the standard Schwarzschild case, due to the charge of the central object. Here, the first term represents the Coulomb force while the second one represents a relativistic-like correction of GR through the coupling between the charge Q(u) and L angular momentum. These new correction terms may be small in general in comparison to their Schwarzschild counterparts. However, one can show that there are possibilities that these terms can be comparable or equal to them. Then, for finding the situations where these forces are comparable to the Newtonian gravitational force and the GR correction term in (36), we define the distances D q 1 and D q 2 corresponding to a q 1 a N 1 and a q 2 a L 1, respectively, in which a N , a L are the Newtonian and the relativistic correction accelerations, respectively, and a q 1 and a q 2 are defined as Accordingly, we obtain the distances D q 1 and D q 2 corresponding to a q 1 a N 1 and a q 2 a L 1, respectively, as • In the third line, we have two new correction terms due to the presence of the surrounding field. Here, the first term is similar to that of Newtonian gravitational term and the second term is similar to the relativistic correction of GR through the coupling between the background filed parameter N s (u) and angular momentum L. Then, we see that for the more realistic non-empty backgrounds, the geodesic equation of any object depends strictly not only on the mass of the central object of the system and the conserved angular momentum of the orbiting body, but also on the background field nature. Similar to the previous case, one can show that there are possibilities that the background correction terms can be comparable to their Schwarzschild counterparts. Thus, for this case, we define the distances D s 1 and D s 2 which correspond to a s 1 a N 1 and a s 2 a L 1, respectively, in which a s 1 and a s 2 are Then, we obtain the distances D s 1 and D s 2 as • The term in the fourth line is also a new non-Newtonian correction resulting from the dynamics of black hole and its surrounding field. It is associated with the radiationaccretion power of the black hole and its surrounding field. 1 Calling this acceleration as the induced acceleration a i by the dynamics, where the subscript i stands for "induced", we have in which, following Lindquist, Schwartz and Misner [64], we define the generalized "total apparent flux" as where L, Q and N are the apparent fluxes associated to the black hole mass, charge and its surrounding field, respectively.
Using these definitions, we can rewrite (41) as This new correction term may also be small in general in comparison to the Newtonian term [64]. However, one can show that there are also possibilities that these two terms can be comparable. Then, we define the distance D i which satisfies a i a N , and it will be given by the solutions of the following equation for different values of M, ω s and apparent fluxes L, Q and N as It is hard to find the general solutions to this equation in terms of its generic parameters L, Q, N, M and ω s . However, we will show that there are possible solutions for the various backgrounds of dust, radiation, quintessence and cosmological constant-like fields for some particular ranges of the parameters.
As we see from (40), the distances D q 1 and D q 1 depend only on the parameters of the black hole, and not on the background parameters, while D s 1 , D s 2 and D i depend on both the black hole and its background parameters. Then, in the following we give some plots denoting the possibility of having D q 1 and D q 1 representing a q 1 a N 1 and a q 2 a L 1, respectively, and postpone the studying of the remaining cases (D s 1 , D s 2 and D i ) till the following subsections. In Then, one can realize that there are possibilities for the phase space of our parameters such that the charge contributions can be comparable to their Schwarzschild counterparts.
In the following subsections, we consider the cosmological surrounding fields of dust, radiation, quintessence and cosmological constant-like fields as the special classes of the obtained general solution (22), and we will investigate some of their interesting features in more detail.

Evaporating-accreting Bonnor-Vaidya black hole surrounded by the dust field
For the dust surrounding field, we set ω d = 0 [53,65]. Then, the metric (22) appears in the following form It is seen that a charged black hole in the dust background appears as an effectively evaporating-accreting charged black hole with an effective mass 2M e f f = 2M(u) + N d (u). Then, the presence of effective mass term changes the thermodynamics, causal structure and Penrose diagrams of the original Bonnor-Vaidya black hole up to a mass re-scaling. The total radiation-accretion density in the dust background is given by and consequently dynamical behaviour of the background dust field at (u, r ) is governed by Then, at an arbitrary distance r from the black hole, the background dust field must obey the above conditions. Interestingly, for the special case ofṄ , there is no pure radiation-accretion density, i.e σ (u, r ) = 0, and then the total energy-momentum tensor (4) will be diagonalized. The case σ (u, r ) = 0 corresponds to two possible physical situations. The first one is related to the situation where for any particular distance r 0 , the background dust field (Ṅ d (u)) and black hole (Ṁ(u) andQ(u)) behave such that their contributions cancel out each others, leading to σ (u, r 0 ) = 0. The second situation is associated with the case where for the given dynamical behaviors of the Bonnor-Vaidya black hole and its surrounding dust field, one can find the particular distance possessing zero energy density, i.e σ (u, r * (u)) = 0. Then, regarding (45)-(47), the following points can be realized for a Bonnor-Vaidya black hole surrounded by the dust field.
• Regarding (45), , we find that the radiation-accretion density vanishes only for r * → ∞. This means that for the emission case, the outgoing charged radiation can penetrate through the dust background so far from the black hole horizon and for the accretion case by the black hole, the black hole affects the so far surrounding dust field.
• Regarding (47), for the case of constant rate forṄ d (u), M(u) andQ(u), the distance r * is fixed to a particular value. In general case whereṄ d (u) andṀ(u) andQ(u) have no constant rates, the r * is a dynamical position with respect to the time coordinate u, i.e r * = r * (u). • Regarding (47), to have a particular distance at which the energy density σ (u, r * ) is zero, the positivity of r * (u) also requires that Q(u)Q (u) and 2Ṁ e f f = 2Ṁ(u) +Ṅ d (u) have the same signs. For the cases in which r * (u) is not positive, the lack of a positive value radial coordinate is interpreted as follows: the total radiation-accretion density σ (u, r ) never and nowhere vanishes. • Regarding (47), demanding thatQ(u) andṀ(u) have the same signs for both of the radiation and accretion processes, the positivity condition of r * (u) requires the condition |2Ṁ(u)| ≥ |Ṅ d (u)| whenṄ d (u) takes opposite sign. • In the case of r * (u) being the positive radial distance, for the given radiation-accretion behaviors of the black hole and its surrounding dust field, i.eṀ(u),Q(u) anḋ N d (u), it is possible to find a distance at which we have no any radiation-accretion energy density contribution. In other words, it turns out that the rate of outgoing radiation energy density of the black hole is exactly balanced by the rate of ingoing absorption rate of surrounding field at the distance r * and vice versa. • Regarding (47), for the case of |Q(u)Q(u)| |Ṁ e f f |, we have r * → ∞. Considering the unit charge gauge, for the extremal caseQ(u) ≈Ṁ(u), for r * → ∞, we find that black hole evolves very slow relative to its background. Then, by satisfaction of these dynamical conditions to have r * → ∞, the positive energy density condition is respected everywhere in the spacetime. In other cases, the positive energy density is respected in some regions while it is violated beyond those regions.
• Another interesting situation happens whenṀ e f f = 0, i.e 2Ṁ(u) = −Ṅ d (u). In this case, regarding (45) and (47), the radiation-accretion density is only resulting and consequently r * → ∞. Also, in order to respect to the positive energy condition here, it is required that and Q(u)Q(u) have opposite signs. , ifṀ e f f = 0, we have r * → ∞.
In the following, we demonstrate the various general situations which can be realized for the Bonnor-Vaidya black hole surrounded by a dust field in the Tables 1 and 2. Here, we assume that the radiation case corresponds toṀ(u) < 0, Q(u) ≤ 0 and the accretion case corresponds toṀ(u) > 0, Q(u) ≥ 0.
Regarding Table 1, we see that for the cases I, II, VII and VIII, there are regions in spacetime that the positive energy condition is respected, while beyond these regions it is violated. The cases IV, V, VIIII and XII are not physical in the sense that the positive energy condition is violated in the whole spacetime. The cases III, VI and XI as well as X represent the situations that the positive energy condition is respected in the whole spacetime with and without a priory condition on the black hole and its surrounding dust field dynamics, respectively.
Regarding Table 2, we see that for the cases I, II, VII and VIII, there are regions in spacetime that the positive energy condition is respected, while beyond these regions it is violated. The cases III, VI, X and XI are not physical in the sense that the positive energy condition is violated in the Accretion/decay of SF by evaporating/vanishing BH Accretion/decay of SF by evaporating/vanishing BH Absorbtion of BH's radiation by SF whole spacetime. The cases IV as well as V, VIIII and XII represent the situations that the positive energy condition is respected in the whole spacetime without and with a priory condition on black hole and its surrounding dust filed dynamics, respectively. Regrading the conditions in the Tables 1 and 2 for = − 1 and = + 1, the behaviour of radiation-accretion density σ (u, r ) in (45)  Also, the Eq. (55) associated with a i a N takes the form of which has the solution Then, we see that how this particular distance depends on the parameters L, Q, N and M. In Fig. 5, we have plotted the solutions of (48) for some typical ranges of L, Q and N parameters. This figure indicates that in the dust background,

Evaporating-accreting Bonnor-Vaidya black hole surrounded by the radiation field
For the radiation surrounding field, we set ω r = 1 3 [53,65]. Then, the metric (22) appears in the following form The positive energy condition on the surrounding radiation field, represented by the relation (21), requires N r (u) 0. By defining the positive structure parameter N r (u) = − N r (u), the metric (50) reads as  N r (u). This result can be interpreted as the positive contribution of the characteristic feature of the surrounding radiation field to the effective charge of the black hole. As the consequence of arising the effective charge, the causal structure and Penrose diagrams for this black hole solution differs from the original Bonnor-Vaidya black hole up to a charge re-scaling.
The total radiation-accretion density is given by and consequently, the dynamical behaviour of the background radiation field is governed by the following conditions Then, at an arbitrary distance r from the black hole, the background field must obey the above conditions regarding the values. Interestingly, for the specific case ofṄ r (u) = 2 rṀ(u) − Q(u)Q(u) , there is no pure radiation-accretion density, i.e σ (u, r ) = 0, and the energy-momentum tensor (4) will be diagonalized. The case of zero energy density corresponds to two possible physical situations. The first one is related to the situation where the observer can be located at any distance r 0 such that the background radiation field (Ṅ r (u)), and black hole (Ṁ(u) andQ(u)) contributions cancel out each others, leading to σ (u, r 0 ) = 0 for a moment or even a period of time. The second situation is associated with the case where for the given dynamical behaviors of the charged black hole and its background, one can find the particular distance possessing zero energy density, i.e σ (u, r * (u)) = 0. Then, regarding (52)-(54), the following points can be realized for a Bonnor-Vaidya black hole surrounded by the radiation field.
• Regarding (52), forṄ r (u) = 2(rṀ(u) − Q(u)Q(u)), we find that the radiation-accretion density vanishes only for r * → ∞. This means that for the emission case, the outgoing charged radiation can penetrate through the radiation background so far from the black hole horizon and for the accretion case by the black hole, the black hole affects the so far surrounding radiation. • Regarding (54), for the case of constant rate forṄ r (u), M(u) andQ(u), the distance r * is fixed to a particular value. In general case whereṄ r (u) andṀ(u) andQ(u) have no constant rates, the r * is a dynamical position with respect to the time coordinate u, i.e r * = r * (u). • Regarding (54), for having a particular distance at which the density σ (u, r * ) is zero, the positivity of r * also requires thatṀ (u) and 2Q e f f (u)Q e f f (u) = 2QQ(u) + N r (u) have the same signs. For the cases in which r * is not positive, the lack of a positive value radial coordinate is interpreted as follows: the radiation-accretion density σ (u, r ) never and nowhere vanishes. • Regarding (54), demanding thatQ(u) andṀ(u) have the same signs for both of the radiation and accretion processes, the positivity condition of r * (u) requires the condition |2Q(u)Q(u)| ≥ |Ṅ r (u)| whenṄ r (u) takes opposite sign. • In the case of r * being the positive radial distance, for the given radiation-accretion behaviors of the black hole and its surrounding field, i.eṀ(u),Q(u) andṄ r (u), it is possible to find a distance at which we have no any radiationaccretion energy density contribution. In other words, it turns out that the rate of outgoing radiation energy density of the black hole is exactly balanced by the rate of ingoing absorption rate of surrounding field at the distance r * and vice versa. • Regarding (54), for the case of |Ṁ(u)| |Q e f f (u)Q e f f (u)|, we have r * → ∞. Considering the unit charge gauge, for the extremal caseQ(u) ≈Ṁ(u), we find that black hole evolves very slow relative to its radiation background. Then, by satisfaction of these dynamical conditions to have r * → ∞, the positive energy density is respected everywhere in the spacetime. In other cases, the positive energy density will be respected in some regions, while it is violated beyond those regions.
• Another interesting situation happens for two different cases as Q e f f (u) = 0 andQ e f f (u) = 0 corresponding to Q = N r = 0 and 2Q(u)Q(u) = −Ṅ r (u), respectively. In these cases, regarding (52) and (54), the radiationaccretion density is only resulting from the black hole mass contribution in the form of σ (u, r ) = 2Ṁ(u) r 2 and consequently must have the same sign asṀ(u) to have a positive energy density. In this case, the radiationaccretion density looks like the original neutral Vaidya solution in an empty space, while the black hole and its background here is completely different, and vanishes as r * → ∞. • Regarding (52), for both of the cases of neutral black hole (Q(u) = 0) and black hole with static charge (Q(u) = 0), we have Then, the positivity of r * demands thatṄ r (u) andṀ(u) have same signs, and for |2Ṁ(u)| |Ṅ r (u)|, we have r * → ∞.
In the following, we demonstrate the various general situations which can be realized for the Bonnor-Vaidya black hole surrounded by the radiation field in the Tables 3 and 4.
Regarding Table 1, we see that for the cases I, II, IV, VI, VII and VIIII, there are regions in spacetime where the positive energy condition is respected, while beyond these regions it is violated. The cases IV, V, VIIII and XII are not physical in the sense that the positive energy condition is violated in the whole spacetime. The cases VIII and X represent the situations where the positive energy condition is respected in the whole spacetime with and without a priory condition on the black hole and its surrounding radiation field dynamics, respectively.
Accretion/decay of SF by evaporating/vanishing BH Accretion/decay of SF by evaporating/vanishing BH

s o r b t i o no fB H ' sr a d i a t i o nb yS F
Accretion/decay of SF by evaporating/vanishing BH  Table 2, we see that for the cases I, II, IV, VI, VII and VIIII, there are regions in spacetime where the positive energy condition is respected, while beyond these regions it is violated. The cases VIII, X are not physical in the sense that the positive energy condition is violated in the whole spacetime. The cases IIII as well as V represent the situ-ations where the positive energy condition is respected in the whole spacetime with and without a priory condition on the black hole and its surrounding radiation filed dynamics, respectively.
Regrading the conditions in the Tables 3 and 4 for = − 1 and = + 1, the behaviour of radiation-accretion density Finally, considering the timelike geodesics, for this case, the distances D s 1 and D s 2 corresponding to a s 1 a N 1 and a s 2 a L 1, respectively, are given by In Fig. 8, the particular distances D s 1 and D s 2 versus some typical ranges of M(u) and N r (u) parameters are plotted. Then, we see that the possibility of the equality of Newtonian force and GR correction terms to the corresponding radiation background field contributions are provided.
Also, the Eq. (55) associated with a i a N takes the following form Then, it admits the following solution We see that how this particular distance depends on the values of parameters L, Q, N and M. In Fig. 9, we have plotted the solutions of (57) for some typical ranges of L, Q and N parameters. According to this figure, depending the parameter values, there are locations where the induced force, resulting from the radiation-accretion phenomena in the radiation background, can be equal to the Newtonian force.

Evaporating-accreting Bonnor-Vaidya black hole surrounded by the quintessence field
In the context of cosmology, the quintessence filed is known as the simplest scalar field dark energy model free of the theoretical problems such as Laplacian instabilities or ghosts. The energy density and the pressure profile of the quintessence field are generally supposed as time varying quantities and depend on the scalar field and its associated potential given by ρ = 1 2φ 2 + V (φ) and p = 1 2φ 2 − V (φ), respectively. Thus, the corresponding quintessence equation of state parameter lies in the range − 1 < ω q < − 1 3 . The static Schwarzschild black hole solution surrounded by a quintessence field was first introduced by Kiselev [53]. Then, this solution was generalized to the Reissner-Nordström case and investigated in [60][61][62]. For the quintessence surrounding field, we set ω q = − 2 3 [53,65]. Then, the metric (1) takes the following form This result is interpreted as the non-trivial contribution of the characteristic feature of the surrounding quintessence field to the metric of the Bonnor-Vaidya black hole. The presence of the background quintessence filed changes the causal struc-ture and Penrose diagrams of this black hole solution in comparison to the charged Vaidya black hole in an empty background. A rather similar effect happens when one immerses an static Schwarzschild in a (anti)-de Sitter background with the difference that here the spacetime tends asymptotically to quintessence rather than (anti)-de Sitter asymptotic state.
Regarding the positive energy density condition for this case, represented by the relation (21), it is required that N q (u) 0. In this case, the radiation density is given by Accretion/decay of SF by evaporating/vanishing BH Based on this relation, the dynamical behaviour of the surrounding quintessence field is governed by Then, at an arbitrary distance r from the black hole, the surrounding quintessence field must obey the above conditions. Interestingly, for the special case ofṄ q (u) = 2 r 3 (Q(u)Q(u)−rṀ(u)), there is no pure radiation-accretion density, i.e σ (u, r ) = 0, and the total energy-momentum tensor (4) will be diagonalized. This means that the black hole and its surrounding quintessence field completely cancel out the effects of each others. This case corresponds to two possible physical situations. The first one is related to the situation where the observer can be located at any distance r 0 such that the backgroundṄ q (u), and black holeṀ(u) andQ(u) contributions cancel out each others, leading to σ (u, r 0 ) = 0 for a moment or even a period of time. The second situation is related to the case where for the given dynamical behaviors of the charged black hole and its quintessence background, one can find the particular distance possessing zero energy density (σ (u, r * (u)) = 0), see the Appendix A for more details. Then, regarding (60)-(62), the following points can be realized for a Bonnor-Vaidya black hole surrounded by the quintessence field.
• Regarding (60), forṄ q (u) = 2 r 3 (Q(u)Q(u) − rṀ(u)), we find that in contrast to the cases of the Bonnor-Vaidya black hole surrounded by the dust and radiation fields, here the radiation-accretion density does not vanish at r * → ∞. This is due to the fact that the spacetime here has the quintessence asymptotic rather than an empty Minkowski.
• Regarding (62), for the case of constant rates forṄ q (u), M(u) andQ(u), the distance r * is fixed to a particular value. In general case whichṄ q (u) andṀ(u) andQ(u) have no constant rates, the r * has a dynamical position with respect to the time coordinate u, i.e r * = r * (u). • Regarding (62), to have a particular distance at which the density σ (u, r * ) is zero, the positivity of r * also requires thatṄ q (u) takes an opposite sign ofṀ(u) (andQ(u)), see (78). This is in agreement with our primary consideration for the signs of dynamical parameters (Ṅ q (u),Ṁ(u) anḋ Q(u)) for the radiation and accretion processes in the previous sections. For the cases in which r * is not positive, the lack of a positive value radial coordinate is interpreted as follows: the radiation-accretion density σ (u, r ) never and nowhere vanishes. • In the case of r * being the positive radial distance, for the given radiation-accretion behaviors of the black hole and its surrounding field, i.eṀ(u),Q(u) andṄ q (u), it is possible to find a distance at which we have no any radiationaccretion energy density contribution. In other words, it turns out that the rate of outgoing radiation energy density of the black hole is exactly balanced by the rate of ingoing absorption rate of surrounding field at the distance r * and vice versa. • Regarding (62) and (78), for both of the black holes with |Ṁ(u)| |Q(u)Q(u)| andṀ(u) → 0, we have r * → ∞ andṄ q (u) → 0. This means that for a black which is almost active only due to its dynamical charge, one can find that (i) there is a non-zero radiation density even at far distance from the black hole and (ii) positive energy condition is respected everywhere.
• Regarding (62) and (78), working in the unit charge gauge, for the extremal caseQ(u) ≈Ṁ(u), we find • Regarding (60), for both of the cases of neutral black hole (Q(u) = 0) and black hole with static charge (Q(u) = 0), we have Then, for |Ṅ q (u)| |2Ṁ(u)|, we have r * → ∞. This means that for an almost static background (the back-ground with negligible dynamics relative to the black hole mass), the zero of the radiation-accretion density lies at infinity and the positive energy density is respected everywhere in the spacetime. In other cases, one can find a finite value for r * representing the zero radiation-accretion density in which the positive energy density will be respected in some regions while it is violated beyond those regions, see our previous work on the neutral black hole case for more details [55,56].
• Regarding both the solutions (62) and (64), the signs ofṀ(u) andṄ q (u) should be opposite to have a zero radiation-accretion density for both of the radiationaccretion processes.
In the following, regarding the obtained solutions and the above discussions, we demonstrate the various situations which can be realized for the Bonnor-Vaidya black hole surrounded by the quintessence field in the Tables 5 and 6.
Regarding Table 5, we see that for the cases III and V, there are regions is spacetime that the positive energy condition is respected, while beyond these regions it is violated. The case II is not physical in the sense that the positive energy condition is violated in the whole spacetime. The case IV is also not physical in the sense that the positive energy condition is violated in whole spacetime except at the zero density point. The cases I and VI represent the situations that the positive energy condition is respected in the whole spacetime without a priory condition on the black hole and its surrounding quintessence filed dynamics.
Regarding Table 6, we see that for the cases III and V, there are regions in spacetime that the positive energy condition is respected, while beyond these regions it is violated. The cases I and VI are not physical in the sense that the positive energy condition is violated in whole spacetime. The cases II and IV represent the situations that the positive energy condition is respected in the whole spacetime without a priory condition on the black hole and its surrounding quintessence filed dynamics. Regrading the conditions in the Tables 5 and 6 for = −1 and = +1, respectively, the behaviour of radiation-accretion density σ (u, r ) in (60) is plotted for some typical values ofṀ(u),Q(u) andṄ q (u) in the Figs. 10 and 11, respectively. Using these plots, one can compare the radiation-accretion density values for the various situations.
Considering the timelike geodesics for this case, the distances D s 1 and D s 2 associated with a s 1 a N 1 and a s 2 a L 1, respectively, are given as In Fig. 12, we have plotted the location of these particular distances versus some typical ranges of the black hole mass Fig. 11 Radiation-accretion density σ versus the distance r for some typical constantṀ,Q andṄ q values for = − 1 in the quintessence background. Here, we have set Q = 1 for simplicity Fig. 12 The variation of D s1 (green plot) and D s2 (yellow plot) versus typical values of the M(u) and N r (u) parameters for the quintessence background M(u) and background quintessence field N q (u) parameters. Then, one finds that there are possibilities for the equality of the Newtonian force and GR correction terms to the corresponding quintessence background field contributions. Also, the Eq. (55) associated with a i a N for this case takes the following form The solutions to (66) in general are complicated and are given in the Appendix B. However, we can demonstrate those solutions in Fig. 13 for some typical values of our L, Q, N and M parameters. This figure represents that depending on the values of our parameter, we can find locations where the induced force, resulting from the radiation-accretion phenomena in the quintessence background, can be equal to the Newtonian force.

Evaporating-accreting Bonnor-Vaidya black hole surrounded by the cosmological field
For the cosmological surrounding field, we set ω c = − 1 [53,65]. Then, the metric (1) takes the following form This result is interpreted as the non-trivial contribution of the characteristic feature of the surrounding cosmological field to the metric of the charged Bonnor-Vaidya black hole. The presence of the background cosmological filed changes the causal structure and Penrose diagrams of this black hole solution in comparison to the black hole in an empty background. The similar effect happens when one immerse an static Schwarzschild black hole in a (anti)-de Sitter background. The positive energy density condition on the surrounding cosmological field, represented by the relation (21), requires N c (u) 0. Then, in this case, N c (u) plays the role of a positive dynamical cosmological field. This case may represents the dynamical black holes in more general cosmological models proposing a time varying cosmological term. The main purpose of these cosmological scenarios is to provide an explanation for the recent observed accelerating expansion of the universe, see [66][67][68][69][70][71][72] as some instances. For the case of N c = constant = > 0, we recover the Bonnor-Vaidya black hole embedded in a de Sitter space obtained by Patino and Rago [18]. The solution in [18] was generalized to the case of the rotating radiating charged black hole in a static de Sitter space in [25]. In [73], the causal structure of the solution obtained in [25] is studied.
In this case, the total radiation-accretion density is given by Then, the dynamical behaviour of the surrounding cosmological field is governed by the following conditions This represents that at an arbitrary distance r from the black hole, the surrounding cosmological field must obey the above conditions. Interestingly, for the special case ofṄ c (u) = 2 r 4 (Q(u)Q(u)−rṀ(u)), there is no pure radiation-accretion density, i.e σ (u, r ) = 0, and the total energy-momentum tensor (4) will be diagonalized. This case corresponds to two possible physical situations. The first one is associated with the situation where the observer can be located at any distance r 0 such that the background cosmological field (Ṅ c (u)), and the black hole (Ṁ(u) andQ(u)) contributions cancel out each others leading to σ (u, r 0 ) = 0 for a moment or even a period of time. The second situation is related to the case that for the given dynamical behaviors of the charged black hole and its cosmological background, one can find the particular distance possessing zero energy density, i.e σ (u, r * (u)) = 0, see the Appendix C for more details. Then, regarding (68)-(70), the following points can be realized for a Bonnor-Vaidya black hole surrounded by the cosmological field.
• Regarding (68), forṄ c (u) = 2 r 4 (Q(u)Q(u) − rṀ(u)), we find that in contrast to the cases of Bonnor-Vaidya black hole surrounded by dust and radiation fields, the radiation-accretion density does not vanish for r * → ∞. This is due to the fact that here the spacetime has the de Sitter asymptotic rather than an empty Minkowski.
• Regarding (70), for the case of constant rates forṄ c (u), M(u) andQ(u), the distance r * is fixed to a particular value. In general case whichṄ c (u) andṀ(u) andQ(u) have no constant rates, the r * has a dynamical position with respect to the time coordinate u, i.e r * = r * (u). • Regarding (70) and (84), to have a particular distance at which the density σ (u, r * ) is zero, the positivity of r * also requires thatṄ c (u) takes an opposite sign of QQ(u) (anḋ M(u)). This is in agreement with our primary consideration for the signs of dynamical parameters (Ṅ c (u),Ṁ(u) andQ(u)) for the radiation and accretion processes in the previous sections. Similarly, for the cases in which r * is not positive, the lack of a positive value radial coordinate is interpreted as follows: the radiation-accretion density σ (u, r ) never and nowhere vanishes. • In the case of r * being a real and positive radial distance, for the given radiation-accretion behaviors of the black hole and its surrounding field, i.eṀ(u),Q(u) anḋ N c (u), it is possible to find a distance at which we have no any radiation-accretion energy density contribution. In other words, it turns out that the rate of outgoing radiation energy density of the black hole is exactly balanced by the rate of ingoing absorption rate of surrounding cosmological field at the distance r * and vice versa. • Regarding (70) and (84), for both of the black holes with |Ṁ(u)| |Q(u)Q(u)| andṀ(u) → 0, we have r * → ∞ andṄ c (u) → 0. This means that for a black which is almost active only due to its dynamical charge, one can find (i) a non-zero radiation density even at far distance from the black hole and (ii) the respected positive energy condition at everywhere.
• Regarding (70) and (84), working in the unit charge gauge, for the extremal caseQ(u) ≈Ṁ(u), we find Then, regarding (63) and (71), comparing the quintessence and cosmological background fields for the extremal case, we have • Regarding (68), for both of the cases of neutral black hole (Q(u) = 0) and black hole with static charge (Q(u) = 0), we have Then, for |Ṅ c (u)| |2Ṁ(u)|, we have r * → ∞. This means that for an almost static background (the background with negligible dynamics relative to the black hole), the zero of the radiation-accretion density lies at infinity and the positive energy density is respected everywhere in the spacetime. In other cases, one can find a finite value for r * representing the total zero accretion density in which the positive energy density will be respected in some regions while violated beyond those regions, see our previous work on the neutral black hole case for more details [55,56].
• Regarding both the solutions (70) and (73), the signs ofṀ(u) andṄ c (u) should be opposite to have a zero radiation-accretion density for both of the radiationaccretion processes. • Regarding (64) and (73), in the case that the quintessence and cosmological backgrounds have a same behavior (|Ṅ q (u)| = |Ṅ c (u)|), for |Ṅ c,q (u)| < |2Ṁ(u)|, we have r * c (u) < r * q (u) while for |Ṅ c,q (u)| > |2Ṁ(u)|, we have r * c (u) > r * q (u). • Regarding (64) and (73) for both of the cases of neutral black hole (Q(u) = 0) and black hole with static charge (Q(u) = 0), whenṀ(u) andṄ q (u) as well asṄ c (u) have the same signs, r * (u) is imaginary and negative for the quintessence and cosmological fields, respectively. Then, for these cases, the radiation-accretion density never be zero and the positive energy condition is respected or violated in the whole spacetime.
In the following, regarding the obtained solutions and the above discussions, we demonstrate the various situations which can be realized for the Bonnor-Vaidya black hole surrounded by the cosmological field in the Tables 7 and 8. Regarding Table 7, we see that for the cases III and V, there are regions in spacetime that the positive energy condition is respected, while beyond these regions it is violated. The case II is not physical in the sense that the positive energy condition is violated in the whole spacetime. The case IV is also not physical in the sense that the positive energy condition is violated in the whole spacetime except at the zero Accretion/decay of SF by evaporating/vanishing BH Not physical density point. The cases I and VI represent the situations that the positive energy condition is respected in the whole spacetime without a priory condition on the black hole and its surrounding cosmological filed dynamics.
Regarding Table 8, we see that for the cases III and V, there are regions in spacetime that the positive energy condition is respected, while beyond these regions it is violated. The cases I and VI are not physical in the sense that the positive energy condition is violated in the whole spacetime. The cases II and IV represent the situations that the positive energy condition is respected in the whole spacetime without a priory condition on the black hole and its surrounding cosmological filed dynamics.
Regrading the conditions in the Tables 7 and 8 for = − 1 and = + 1, the behaviour of radiation-accretion density σ (u, r ) in (68) is plotted for some typical values ofṀ,Q andṄ c in the Figs. 14 and 15, respectively. Using these plots, one can compare the radiation-accretion density values for the various situations.
Concerning the timelike geodesics, the distances D s 1 and D s 2 associated with a s 1 a N 1 and a s 2 a L 1, respectively, read as According to (39), the case of D s 2 → ∞ is resulting from the fact that, in contrast to black hole itself, the cosmological constant-like field does not couple to angular momentum L, see a s 2 . Then, there is no similar effect to the GR correction term for the cosmological constant-like field. The location of the particular distance D s 1 for some typical ranges of the black hole mass M(u) and background cosmological constant-like field N c (u) parameters is plotted in Fig.  16. Thus, we realize that the possibilities for the equality of the Newtonian force to the corresponding cosmological constant-like background field contributions are provided Moreover, the Eq. (55) associated with a i a N takes the form of The solutions to (75) in general are complicated and are given in the Appendix D. However, we can demonstrate those solutions in Fig. 17 for some typical values of our L, Q, N and M parameters. From this figure, we see that, depending on the values of our parameter, we can find locations where the induced force of the radiation-accretion phenomena in the cosmological background can be equal to the Newtonian one.

Summary and concluding remarks
By generalizing our previous work [55,56], we have constructed a general solution for the classical description of the evaporating-accreting charged black holes in the generic dynamical backgrounds of dust, radiation, quintessence and cosmological constant, namely as the surrounded Bonnor-Vaidya black holes. We have shown that (i) the original Bonnor-Vaidya solution can be recovered by turning off the background field, and (ii) the charged Kiselev static solution can be obtained as another subclass of our general solution in • For ω s < 0, the charge contribution is dominant near the black hole while for the far distances, we have |σ Q | < |σ M | < |σ N s |, meaning that the surrounding field contribution is dominant at large distances. • For 0 < ω s < 1/3, the charge contribution is dominant near the black hole while for the far distances, we have |σ Q | < |σ N s | < |σ M |, meaning that the black hole mass contribution is dominant at large distances. • For ω s > 1/3, the surrounding field contribution is dominant near the black hole while for the far distances, we have |σ N s | < |σ Q | < |σ M | meaning that the black hole mass contribution is dominant again at large distances.
We have addressed some interesting features of these solutions and classified them according to their behaviors under imposing the positive energy condition. We have discussed that this condition gives some severe restrictions on the black hole and its background field dynamics. Using this condition, we have found a particular distance possessing zero effective radiation-accretion density. For all the Bonnor-Vaidya black holes surrounded by the dust, radiation, quintessence and cosmological constant fields, we have found a corre- sponding particular distance r * (u) possessing zero effective radiation-accretion density, i.e σ (u, r * ) = 0. This distance has a dynamical location in general, except for the case of constant rates for the black hole and its surrounding field parameters, i.eṀ(u),Q(u) andṄ s (u). For the cases in which there is no real and positive r * (u), the interpretation is that the total radiation-accretion density σ (u, r ) never and nowhere vanishes. Also, we have studied the timelike geodesics for the obtained solutions and have found that three new correction terms arise relative to the case of the standard Schwarzschild black hole. The first kind of corrections is due to the charge of the central black hole which includes two terms in which its first term is rather similar to the term of Newtonian gravitational force, while its second term is similar to the relativistic correction of GR due to the coupling of charge and angular momentum. The second kind of our correction terms is resulting from the presence of the background fields which surround the black hole. This corrections also include two terms in which the first term looks like to Newtonian gravitational term, while the second term is similar to the relativistic correction of GR due to the coupling of background field parameter to the angular momentum. We have discussed that for the various background fields, there are possibilities for the equality of Newtonian and GR correction terms to the corresponding charge and background field contributions. We have given some plots denoting these possibilities for each case. The third type of new corrections, as a non-Newtonian correction, is resulting from the dynamics of the central black hole and its surrounding field. For this case also, it is shown that depending on the dynamical features of black hole and its background, one can find possibilities that this new dynamical correction can be compared with the Newtonian case. We have given some plots denoting these particular situations. Accordingly, we realize that for more realistic cases which are non-static and possess non-empty backgrounds in nature, the geodesic equation of any object depends strictly not only on the (i) mass of the central object of the system, (ii) the angular momentum of the orbiting body, but also on (iii) background field features and (iv) black hole-background field dynamics. We have summarized our obtained results for the various surrounding fields as follows.
• For the dust background, we have the Bonnor-Vaidya black hole with the effective dynamical mass 2M e f f = 2M(u) + N d (u). To have a particular distance r * at which σ (u, r * ) = 0, the condition |2Ṁ(u)| ≥ |Ṅ d (u)| is required. For the case of |Q(u)Q(u)| |Ṁ e f f |, we have r * → ∞. In the unit charge gauge, for the extremal case (Q(u) ≈Ṁ(u)), for r * → ∞, we find that black hole evolves very slow relative to its dust background. Then, by satisfaction of these dynamical conditions to have r * → ∞, the positive energy density condition is respected everywhere in the spacetime. In other cases, the positive energy density is respected in some regions while it is violated beyond those regions. Another interesting situation happens whenṀ e f f = 0. In this case, the radiation-accretion density is only resulting from the charge contribution with r * → ∞.
• For the radiation background, we have the Bonnor-Vaidya black hole with the effective dynamical charge Q e f f (u) = Q 2 (u) + N r (u). To have a particular distance r * at which σ (u, r * ) = 0, the condition |2Q(u)Q(u)| ≥ |Ṅ r (u)| is required. For the case of |Ṁ(u)| |Q e f f (u)Q e f f (u)|, we have r * → ∞. In the unit charge gauge, for the extremal case (Q(u) ≈Ṁ(u)), we find that black hole evolves very slow relative to its radiation background. Then, by satisfaction of these dynamical conditions to have r * → ∞, the positive energy density is respected everywhere in the spacetime. In other cases, the positive energy density will be respected in some regions, while it is violated beyond those regions. Another interesting situation happens for two different cases as Q e f f (u) = 0 andQ e f f (u) = 0 corresponding to Q = N r = 0 and 2Q(u)Q(u) = −Ṅ r (u), respectively. In these cases, the radiation-accretion density is only resulting from the black hole mass contribution. For both of the cases of neutral black hole (Q(u) = 0) and black hole with static charge (Q(u) = 0), for |2Ṁ(u)| |Ṅ r (u)|, we have r * → ∞. • For the quintessence background, the positive energy density condition demands N q (u) 0. To have a particular distance r * at which σ (u, r * ) = 0, the positivity of r * requires thatṄ q (u) takes an opposite sign ofṀ(u) (andQ(u)). For both of the black holes with |Ṁ(u)| |Q(u)Q(u)| andṀ(u) → 0, we have r * → ∞ anḋ N q (u) → 0. This means that for a black which is almost active only due to its dynamical charge, one can find that (i) there is a non-zero radiation density even at far distance from the black hole and (ii) positive energy condition is respected everywhere. Here, in the unit charge gauge, for the extremal case (Q(u) ≈Ṁ(u)), we find r * → 3 2 andṄ q (u) → − 8 27Ṁ (u). For both of the cases of neutral black hole (Q(u) = 0) and black hole with static charge (Q(u) = 0), for |Ṅ q (u)| |2Ṁ(u)|, we have r * → ∞. This means that for an almost static background (the background with negligible dynamics relative to the black hole mass), the zero of the radiationaccretion density lies at infinity and the positive energy density is respected everywhere in the spacetime.
• For the cosmological background, the positive energy density condition demands N c (u) 0 representing (dynamical) de Sitter space. To have a particular distance r * at which σ (u, r * ) = 0, the positivity of r * requires thaṫ N c (u) with respect to QQ(u) (andṀ(u)) takes an opposite sign. For both of the black holes with |Ṁ(u)| |Q(u)Q(u)| andṀ(u) → 0, we have r * → ∞ anḋ N c (u) → 0. This means that for a black which is almost active only due to its dynamical charge, one can find (i) a non-zero radiation density even at far distance from the black hole and (ii) the respected positive energy condition at everywhere. In the unit charge gauge, for the extremal case (Q(u) ≈Ṁ(u)), we find r * → 4 3 andṄ c (u) → − 27 128Ṁ (u). Then, comparing the quintessence and cosmological background fields for the extremal case, we have r * c < r * q and |Ṅ c (u)| < |Ṅ q (u)|. Also, for both of the cases of neutral black hole (Q(u) = 0) and black hole with static charge (Q(u) = 0), for |Ṅ c (u)| |2Ṁ(u)|, we have r * → ∞. This means that for an almost static background (the background with negligible dynamics relative to the black hole), the zero of the radiationaccretion density lies at infinity and the positive energy density is respected everywhere in the spacetime. In the case that the quintessence and cosmological backgrounds have a same behavior (|Ṅ q (u)| = |Ṅ c (u)|), for |Ṅ c,q (u)| < |2Ṁ(u)|, we have r * c (u) < r * q (u) while for |Ṅ c,q (u)| > |2Ṁ(u)|, we have r * c (u) > r * q (u). Finally, for both of the cases of neutral black hole (Q(u) = 0) and black hole with static charge (Q(u) = 0), whenṀ(u) andṄ q (u) as well asṄ c (u) have the same signs, r * (u) is imaginary and negative for the quintessence and cosmological fields, respectively. Then, for these cases, the radiation-accretion density never be zero and the positive energy condition is respected or violated in the whole spacetime.
For all of the mentioned surrounding fields, by classifying the possible situations in Tables 1, 2, 3, 4, 5, 6, 7 and 8, we find that (i) there are situations possessing regions in spacetime that the positive energy condition is respected, while beyond these regions it is violated. (ii) there are some other situations that are not physical in the sense that the positive energy condition is violated in the whole spacetime, and (iii) there are also possible situations such that the positive energy condition is respected in the whole spacetime with or without a priory condition on black hole and its surrounding filed dynamics. To depict these situations, we have given some plots (Figs. 2 and 3 for dust, 6 and 7 for radiation, 10 and 11 for quintessence as well as 14 and 15 for the cosmological constant-like background) for the radiation-accretion density versus the distance r for some typical values ofṄ s (u),Ṁ(u) andQ(u) in the different backgrounds. Also, considering the timelike geodesics, the possibility of (i) the equality of the Newtonian force as well as GR correction terms to the corresponding contributions of various backgrounds, and (ii) the equality of the induced force resulting from the radiationaccretion phenomena to the Newtonian force are shown in the Figs. 4 and 5 for the dust background, in 8 and 9 for the radiation background, in 12 and 13 for the quintessence background, and in 16 and 17 for the cosmological constant-like background.
In our next work, we aim to report elsewhere on the causal structures, horizon and thermodynamical properties of our obtained solutions. where for σ (u, r ) = 0 where σ (u, r ) is given by (60). Due to the complexity of these solutions, one may consider the special case oḟ Then, Eq. (77) simplifies to and the solutions (76) reduce to Finally, regarding the fact thatQ(u) andṀ(u) have the same signs in both of the radiation and accretion processes, the only positive solution representing a zero radiation-accretion energy density will be r * (u) = r 2 (u) = r 3 (u). Note that here, the presence of negative solution r 1 (u) does not change the total physical properties of the solution, because there is a positive solution r * (u) = r 2 (u) = r 3 (u) for radial distance which possesses zero energy density.