Surrounded Vaidya black holes: apparent horizon properties

We study the thermodynamical features and dynamical evolutions of various apparent horizons associated with the Vaidya evaporating black hole surrounded by the cosmological fields of dust, radiation, quintessence, cosmological constant-like and phantom. In this regard, we address in detail how do these surrounding fields contribute to the characteristic features of a surrounded dynamical black hole in comparison to a dynamical black hole in an empty background.


Introduction
The Vaidya solution [1,2] is one of the non-static solutions of Einstein field equations. This solution can be regarded as a dynamical generalization of the static Schwarzschild black hole solution. This solution is characterized by a dynamical mass function, i.e m = m(u) where u is the advanced time coordinate, and an outgoing null radiation flow. One can interpret this null radiation as high frequency electromagnetic or gravitational waves, massless scalar particles or neutrinos. Based on these properties, this solution has been used for studying the spherical symmetric gravitational collapse and for exploring the validity of the cosmic censorship conjecture [3,4,5,6], as a model for studying evaporating black holes as well as Hawking radiation [7,8,9,10,11], among the other applications. This solution has been generalized to the Vaidya-de Sitter solution [12], Bonnor-Vaidya solution [13], Bonnor-Vaidya-de Sitter solution [14,15,16,17,18], radiating dyon solution [19], and the surrounded Vaidya and Bonnor-Vaidya solutions [20,21], see also [22,23] for the other similar generalizations.
The connections between black holes and thermodynamics were initially discovered by Hawking [24] and Bekenstein [25]. Black holes behave as thermodynamic objects emitting radiation, the well-known Hawking radiation, from their event horizons which is described in the context of quantum field theory in curved spacetime. From thermodynamical point of view, there are (i) a characteristic temperature for a black hole which is proportional to the surface gravity associated with its event horizon (EH), and (ii) an entropy equal to one quarter of the area of the EH. Interestingly, the apparent horizon (AH) and the event horizon EH coincide for the static or stationary black holes. Then, one cannot distinguish exactly where is the actual emission surface for the Hawking radiation; AH or EH? For the case of Vaidya black hole, it is found in [26] that the Hawking radiation originates always in regions having geometry near to the geometry in a close neighborhood of AH. Then, EH appears as radiating surface only if it passes through such regions, similar to what happens for the case of static Schwarzschild black hole. Considering the collapse of a spherical shell leading to a black hole with the Vaidya solution [27], it is suggested that in order to maintain a local thermodynamics, the entropy of the black hole should be identified as A AH /4 where A AH is the area of the black hole's AH. Also, the fact that thermodynamics of a dynamical black hole can be built successfully on the AH is verified in [28,29].
In this paper, we aim to study the properties of the various horizons, mainly the AH, associated with the Vaidya evaporating black hole solution surrounded by the cosmological fields of dust, radiation, quintessence, cosmological constant-like and phantom. In this regard, we explore the surrounded Vaidya solutions introduced in [20] using the approach used in [30,31]. The organization of the paper is as follows. In section 2, we review briefly the surrounded Vaidya black hole solutions introduced in [20]. In section 3, we discuss on the proprieties of these solutions in a general form. In section 4, we study in detail the properties of these solutions for the surrounding fields of dust, radiation, quintessence, cosmological constant-like and phantom. In section 5, we give our concluding remarks with introducing our two new works.

The Surrounded Evaporating Vaidya Black Holes
The Vaidya black hole surrounded by generic fields in the advanced time coordinate u is obtained in [20] as where M (u), N s (u) and ω s are the time dependant decreasing mass of the evaporating Vaidya black hole, the normalization parameter and the equation of state parameter of the surrounding field, respectively. The corresponding total energy-momentum tensor T µ ν is given by [20,21,32,33] T µ ν = τ µ ν + T µ ν , where τ µ ν is the energy-momentum tensor associated to the Vaidya null radiation as such that σ = σ(u, r) is the density of the "outgoing radiation flow" and k µ is a null vector field directed radially outward, and T µ ν is the energy-momentum tensor of the surrounding dynamical field as T 0 0 = T 1 1 = −ρ s (u, r), where the energy density ρ s (u, r) of the surrounding field is given by The weak energy condition on the energy density (5) of the surrounding field requires implying that for the surrounding fields with ω s ≥ 0 and ω s ≤ 0, we need N s (u) ≤ 0 and N s (u) ≥ 0, respectively. Regarding (1) and (2), the original Vaidya evaporating black hole solution can be recovered by turning off the background field, i.e setting ρ s (u, r) = 0.
Following [30,35,36], we can rewrite the metric of our surrounded Vaidya black hole spacetime (1), using the null-vector decomposition, in the following form g µν = −n µ l ν − l µ n ν + h µν , where n µ and l µ are ingoing and outgoing null vectors, respectively, given by n µ = −δ u µ , and h µν = r 2 δ θ µ δ θ µ + r 2 sin 2 θδ φ µ δ φ µ , is the two dimensional metric for u = constant and r = constant surfaces such that n µ , l µ and h µν satisfy the following normalization conditions n µ n µ = l µ l µ = 0, l µ n µ = −1, As we will see in the following section, the null vector decomposition of the spacetime metric is more appropriate for studying its various geometric properties.

The Location and Properties of the Horizons
It is a well known fact that the black holes in general have not single horizons to characterize completely their local and global structures, see [30,34] as instances. They possess three important loci: (i) the future EH which is a null three-surface describing the locus of outgoing future-directed null geodesic rays that never reach arbitrarily large distances from the black hole, (ii) the AH which is the outermost marginally trapped three-surface for the outgoing photons. Classically, the AH can be either null or spacelike surfaces, and consequently it can move causally or non-causally. Considering the black hole quantum radiance, the AH can also be a timelike surface, and (iii) the timelike limit surface (TLS) which is a horizon-like locus. For the black holes possessing an exact timelike Killing vector field ∂/∂t in their exterior regions, TLS is the surface given by g(∂ t , ∂ t ) = g tt = 0. For black holes with a small dimensionless accretion or luminosity L = − dM du = −Ṁ (u), one can define the TLS (or quasi-static limit) as the locus where g(∂ u , ∂ u ) = g uu = 0, ∂/∂u being the timelike vector field for an observer sitting at rest at the far distance from the black hole. Generally, the TLS can be a null, timelike, or a spacelike three-surface. We note that when L refers to the quantum radiance, i.e given by the Hawking relation L = L H ∝ ℏ/M 2 , it is usually very small. In this case, it represents a black hole radiating an small amount of quanta in any reasonable time interval. However, the fluctuations in the radiation could be very large, at least comparable to L itself. Thus, we consider L as an average value such that the terms with O(L 2 ) are negligible. We also mention that for the spherical symmetric solutions of R µν = 0, there is a degeneracy for the horizons which means that all these three horizons coincide at r = 2M , known as the Schwarzschild radius. In general, such a coincidence of various horizons happens for the stationary situations. However, in the case of small accretion or luminosity, i.e for R µν = O(L), this degeneracy is partially lifted even for the spherical symmetric spacetimes. For the case of R µν = O(L), we have AH = T LS but EH = AH. For the case of R µν = 0 and without spherical symmetry but preserving stationarity, as for a Kerr hole, we have EH = AH but EH = T LS. So, we see that these three horizons do not coincide generally and they are very sensitive to small perturbations. Moreover, by including the generic background fields in our solution (1), we will considerṄ s (u) here also as an average value and the terms with O(Ṅ 2 s ) are negligible. In general, the structure and dynamics of various horizons for dynamical spacetimes can be understood using the behavior of null geodesic congruences, governed by the Raychaudhuri equation [36] where K, Θ, σ 2 = σ µν σ µν , ω 2 = ω µν ω µν and R µν are the surface gravity, expansion scalar, shear scalar, vorticity scalar, and Ricci tensor, respectively. The expansion scalar for the null rays parameterized by the parameter u is given by [30,37] where ∇ µ is the covariant derivative operator corresponding to the metric (1). The surface gravity is given by the non-affinely parameterized geodesic equation with the null tangent vector l µ as [37,34,38,39] l µ ∇ µ l ν = Kl ν .
Using the normalization conditions (10) on the null vectors n µ and l µ , one can obtain the surface gravity in (13) in the following form Spherically symmetric non-rotational spacetimes, such as our metric (1), are vorticity and shear free, i.e σ µν = 0 = ω µν . Then, for these spacetimes, the location, structure and dynamics of the various horizons depend on the surface gravity K, the expansion scalar Θ and the Ricci tensor R µν . One can follow York [30] and obtain the location of the AH and EH for a dynamical radiating black hole up to O(L) as • AH is defined at the surface such that Θ ≃ 0, • EH is the null surface such that dΘ du ≃ 0. Thus, one realizes that the EH can be obtained using the Raychaudhuri equation (11), and the AH can be determined by the equation (12).
Using the equations (8) and (14), one can obtain the general dynamical surface gravity associated to our metric (1) as Then, by using the equations (8), (12), and (15), the expansion scalar can be obtained as Regarding (1) and (16), we see that the solution to g(∂ u , ∂ u ) = g uu = 0 representing the TLS is the same as the solution to Θ = 0 representing the AH location except the trivial solution r = ∞ for the latter. Then, the AH and the TLS coincide for our spacetime (1). The Hawking temperature on the horizon (apparent or event) of black hole can be determined using its corresponding surface gravity as Then, regarding (15), we see that temperature is also a dynamical quantity with respect to the time coordinate u. Also, the semiclassical Bekenstein-Hawking entropy at the horizon is given by [25] where A H is the area of the horizon given by In the following of this section, we show that the condition dΘ du ≃ 0 determining the location of the EH up to O(L) is equivalent to the requirement that the acceleration of the congruences of null geodesics vanishes at the EH, i.er To prove this, using the metric (1), for a radial null geodesic parametrized by the parameter u, we havė and consequently, we find Then, using the equations (15), (16) and (22), we can rewrite the equation (20) as On the other hand, using the metric (1) and the null vector l µ defined in (8), we find Thus, by using the equations (23) and (24) and neglecting Θ 2 = O(L 2 ,Ṅ 2 s , LṄ s ) in the Raychaudhuri equation (11), we arrive at dΘ du | rEH ≃ 0.
Hence, the condition (25) determining the location of the event horizon becomes equivalent to solving the equation (23) in O(L,Ṅ s ) resulting from the vanishing acceleration of the congruences of null geodesics at EH. Then, the discussion in [30] for the EH location is also valid and verified here for the surrounded evaporating black holes. As mentioned previously, in contrast to the stationary spacetimes, the local definitions of the various horizons do not necessarily coincide with the location of the EH for dynamical black holes [34]. For such dynamical spacetimes, one is left with the question: "For which surface should one define the black hole area, surface gravity, temperature or entropy? " The canonical choice is to use the EH. However, as mentioned in the introduction, there are some evidences that it is the AH, and not the EH, that plays the key role in the Hawking radiation [26,27,28,29], see also [40,41,42]. This finding has became a key point in hopes to demonstrate the Hawking radiation in the laboratory using the models of analogue gravity [43]. Thus, regarding the importance of AH, we first introduce some specific subclasses of our general surrounded Vaidya black hole solution (1) and then, we investigate in detail the AH properties of these solutions including the AH location and dynamics as well as their associated Hawking temperature and Bekenstein-Hawking entropy. We consider the detailed study of EHs as our next work.

Apparent Horizon Properties for the Evaporating Vaidya Black
Hole in Generic Backgrounds

Evaporating Vaidya Black Hole in an Empty Background
In this section, we first explore the AH properties of the Vaidya black hole in empty background. This provides us with the possibility of comparing the properties of this black hole in the various backgrounds with its original form in an empty space. The case of Vaidya black hole in an empty space corresponds to N s (u) = 0. Then, the metric (1) takes the original Vaidya solution form as The surface gravity associated to the metric (26) is given by where the expansion scalar for the outgoing null geodesics reads as Neglecting the trivial solution of Θ = 0 at r = ∞, the location of the AH is given by Then, the dynamics of the single AH (29) is governed by the relatioṅ This shows that the single AH is moving inward during the evaporation, i.e it is shrinking. For this case, the nature of the AH (29) can be determined by looking at the induced metric on the horizon as showing that AH is a timelike surface for an evaporating BH. The Hawking temperature on the AH (29) can be determined using the surface gravity as This relation implies that the temperature associated with the AH (29) rises up by decreasing the mass of black hole during the evaporation. Finally, the semiclassical Bekenstein-Hawking entropy at AH (29) is given by

Evaporating Vaidya Black Hole Surrounded by the Dust Field
For the case of Vaidya black hole surrounded by the dust field with the equation of state parameter ω d = 0 [32,44], the metric (1) takes the following form where N d (u) is the normalization parameter for the dust field surrounding the back hole, with the dimension of [N d ] = l where l denotes the length. Regarding (34), the black hole in the dust background appears as a black hole with an effective dynamical mass 2M ef f (u) = 2M (u) + N d (u). In the following, we see that how the presence of this effective mass modifies the thermodynamics, causal structure and Penrose diagrams of the original Vaidya black hole. The surface gravity associated to the metric (34) is given by We see that the dust background appears with a positive surface gravity which increases the gravitational attraction of the black hole. The expansion scalar for the outgoing null geodesics reads as Then, we find the solution to Θ = 0 determining the location of the single physical AH as where r AH (u) corresponds to a modified larger AH compared with the case of black hole in an empty space given by (29). Then, the surrounding background dust field contributes positively to the black hole mass and size of AH . The dynamics of the AH (37) is governed by the relatioṅ Then, the AH behaviours for the cases of evaporating black hole in a static and dynamical dust backgrounds are different. For this case, we have the following possibilities: • ForṄ d (u) = 0, representing the static surrounding dust field, the AH (37) is shrinking. Regarding (29), this case also shows that the dynamics of AH for an evaporating black hole surrounded by a static dust field is the same as the black hole in the empty space.
• ForṄ d (u) = 2L, the AH (37) is frozen. This case shows that the black hole input by the surrounding field is equal to its output radiation. (37) is expanding. This case requiresṄ d (u) > 2L > 0 which shows an evaporating black hole in a collapsing dust field.
showing an evaporating black hole in a slowly collapsing or a decaying dust field, respectively. The AH (37) with 0 <Ṅ d (u) < 2L orṄ d (u) < 0, shrinks slower or faster than the AH of black hole in a static dust background.
The nature of AH (37) can be determined by looking at the induced metric on it as which implies that the AH (37) is a timelike surface for an evaporating black hole in the static dust field and for a dust field with −2L +Ṅ d (u) < 0, whereas it is a spacelike surface for a dust field with −2L +Ṅ d (u) > 0.
For the specific case ofṄ d (u) = 2L, AH (37) is a null surface. The Hawking temperature of the AH (37) of black hole is given by Then, we observe that the temperature associated with the AH (37) for the black hole in dust field is lower than the temperature of AH (29) of the black hole in an empty space. Finally, the semiclassical Bekenstein-Hawking entropy at the AH (37) is given by which implies that the entropy associated with the AH (37) in a dust field background is higher than the case of the black hole in an empty space.

Evaporating Vaidya Black Hole Surrounded by the Radiation Field
For the case of Vaidya black hole surrounded by a radiation field with the equation of state parameter ω r = 1 3 [32,44], the metric (1) takes the following form where N r (u) denotes the normalization parameter for the radiation field surrounding the black hole, with the dimension of [N r ] = l 2 . Regarding the positive energy condition on the surrounding radiation field, represented by the relation (6), it is required that N r (u) < 0. Then, by defining the positive parameter N r (u) = −N r (u), we arrive at the following metric This metric looks like a radiating charged Vaidya black hole with the dynamical charge Q 2 (u) = N r (u). Then, this result can be interpreted as the positive contribution of the characteristic feature of the surrounding radiation field as the effective charge-like term for the Vaidya black hole, see the Bonnor-Vaidya solution [13].
The appearance of such a charge-like term in the metric changes the causal structure and Penrose diagrams of this solution in comparison to the original Vaidya black hole. A similar effect in the causal structure of spacetime happens when one adds charge to the static Schwarzschild black hole leading to Reissner-Nordström black hole. Then, turning off the background radiation field which surrounds the dynamical Vaidya black hole is equal to turning off the charge in the static Reissner-Nordström case. The surface gravity associated to the metric (43) is given by Then, we see that the radiation background appears with a negative surface gravity which decreases the gravitational attraction of the black hole, exactly the same as the charge effect in Reissner-Nordström black hole. The expansion scalar for the outgoing null geodesics can be obtained as We can find the solutions to Θ = 0 determining the locations of the AH as and Then, depending on the sign of the discriminant ∆(u) = M 2 (u) − N r (u), we have the following three situations.
• For ∆(u) < 0, we have no real solutions. Then, there is no physical AH.
• For ∆(u) > 0, there are two real physical AHs as r AH − (u) = r 1 (u) and r AH + (u) = r 2 (u) for the black hole in the radiation background in contrast to the black hole in the empty space and in the dust field background having a single AH.
For the extremal case of ∆(u) = 0, there is only one single extremal AH given by which is half of the dynamical Schwarzschild radius. Then, the dynamics of the single extremal AH (48) is governed by the relationṙ This shows that the single extremal AH (48) for the black hole in the radiation background is moving inward during the evaporation, i.e it is shrinking. For this case, the nature of extremal AH (48) can be determined by looking at the induced metric on it as showing that the extremal AH (48) is a timelike surface for an evaporating BH in the radiation background. The Hawking temperature on the extremal AH (48) can be determined using the surface gravity as showing that the extremal AH (48) is a zero temperature surface. One may interpret this zero temperature as the suppressed or zero Hawking radiation for the extremal black hole in the radiation background. Finally, the semiclassical Bekenstein-Hawking entropy at AH (48) is given by Then, the entropy of the extremal AH (48) is a quarter of the entropy associated to the single AH of the black hole in an empty space. For the case of ∆(u) > 0, the location of the AHs (46) and (47) for the dilute radiation background, i.e. N r (u) ≪ M (u), are given by and It is seen that the contribution of the surrounding background radiation field leads to (i) the existence of a small inner AH (53) vanishing with N r (u) → 0, and (ii) an smaller outer AH (54) in comparison to the case of the AH (29) for the black hole in an empty background in which at the limit of N r (u) → 0 we recover the dynamical Schwarzschild radius r(u) = 2M (u). The dynamics of the AHs (53) and (54) are given bẏ andṙ Then, the AH behaviours for the cases of evaporating black hole in a static and dynamical radiation backgrounds are different. For this case, we have the following possibilities: • ForṄ r (u) = 0, representing the static surrounding radiation field, the inner AH (53) is moving outward (is expanding) while the outer AH (54) is moving inward (is shrinking) during the evaporation. By stopping the evaporation, if possible, both of the horizons will be frozen, i.
Another interesting result is that regarding (29) and (56), one realizes that the outer AH of the black hole in the static radiation field (54) is shrinking faster than the case of black hole in an empty background.
• ForṄ r (u) > 0, the inner and outer AHs (53) and (54) are moving outward and inward, respectively, but both moving faster than the AH of the black hole in the static radiation background. This case represents an evaporating black hole in a collapsing radiation field.
• ForṄ r (u) < 0, it is possible that the inner and outer AHs (53) and (54) move again outward and inward, respectively, but both will be slower than the AH of the black hole in the radiation backgrounds withṄ r (u) = 0 andṄ r (u) > 0. In the particular case of L N r (u) = −Ṅ r (u)M (u), the inner AH (53) will be frozen while the outer AH (54) will move inward exactly the same as AH (29) of the black hole in an empty background. Reversely, if the termṄ r (u) 2M(u) overcomes the other terms in (55) and (56), the inner and outer AHs (53) and (54) will move inward and outward, respectively.
The nature of AHs (53) and (54) can be determined by looking at the induced metrics on them up to O(N r ) as and Then, forṄ r (u) ≥ 0 the inner and outer AHs (53) and (54) are spacelike and timelike respectively while forṄ r (u) < 0, we see that the natures of AHs (53) and (54) depend on the value ofṄ r (u) 2M(u) term relative to the other terms in the above parentheses. For the totally positive parentheses, AHs (53) and (54) are spacelike whereas they are timelike for the totally negative parentheses . For the specific case of L N r (u) = −Ṅ r (u)M (u), the inner AH (53) will be null, whereas the outer AH (54) will be timelike.
The Hawking temperature on the AHs (53) and (54) of the black hole can be determined as and Regarding (59) and (60), one realizes that the temperature of the inner AH (53) is always negative whereas the temperature of the outer AH (54) is positive. This situation for the temperatures of the AHs is exactly similar to what happens for the Reissner-Nordström black hole where for the region r AH − < r < r AH + , the radial coordinate r is timelike and time coordinate t is spacelike. One may write the geodesic equation and see that the gravity force is repulsive at r = r − AH for the Reissner-Nordström black hole. That is, the inner horizon is repelling stuff while outer horizon is attracting them. In fact, the horizon temperature is related to surface gravity which is the gravitational acceleration of an object falling to the black hole (as measured by the asymptotic observer). Positivity of temperature then means that gravity force is attractive and vice versa. In our case, this point can be realized by the surface gravity given by (44). Here, the negative contribution of the background radiation field, playing the same role as the charge in the Reissner-Nordström black hole, leads to a negative surface gravity for small r values where the inner AH (53) lies there. This negative surface gravity yields the gravitational repulsion and the negative temperature (59) for the inner AH (53). One can also discuss about the gravitational repulsion of the inner AH (53), preventing anything to reach it, and its spacelike nature obtained in (57). Also regarding (32) and (60), we see that the temperature of the outer AH (54) of the black hole in the radiation background is lower than the temperature of the AH (29) in an empty background.
The semiclassical Bekenstein-Hawking entropy at the AHs (53) and (54) are given by and Comparing (61) and (62), the inner AH (53) has a negligible entropy relative to the entropy associated with the outer AH (54) in a dilute radiation background. Also regrading (62), we see that the background radiation field contributes negatively in the entropy of the outer AH leading to a lower entropy compared to the entropy (33) of the single AH (29) in an empty background. In the limit of N r (u) → 0, the entropy of the inner AH (53) vanishes as the horizon itself vanishes, while the entropy (62) of the outer AH (54) reduces to the entropy (33) of the evaporating black hole in an empty space.

Evaporating Vaidya Black Hole Surrounded by the Quintessence Field
For the case of Vaidya black hole surrounded by the quintessence field with the equation of state parameter ω q = − 2 3 [32,44], the metric (1) takes the following form where N q (u) is the normalization parameter for the quintessence field surrounding the black hole, with the dimension of [N q ] = l −1 . Regarding the positive energy condition on the surrounding quintessence field, represented by the relation (6), it is required that N q (u) > 0. Regarding the metric (63), one realizes the non-trivial contribution of the characteristic feature of the surrounding quintessence field to the metric of the Vaidya black hole. In this case, the background quintessence filed changes the causal structure of this black hole solution in comparison to the original Vaidya black hole in an empty background. A rather similar effect happens when one immerses a Vaidya black hole in a de Sitter background [12] with the difference that here our spacetime tends asymptotically to the quintessence rather than the de Sitter asymptotics.
The surface gravity associated to the metric (63) is given by Then, we see that the quintessence background appears with a negative surface gravity contribution which decreases the gravitational attraction of the black hole. This characteristic feature of the quintessence field, namely the gravitational repulsion, is favored by its cosmological application as one of the candidates for the dark energy which is responsible for the current accelerating expansion of the universe [44]. The expansion scalar for the outgoing null geodesics reads as Then, one can find the solutions to Θ = 0 determining the locations of the AH as and Then, depending on the sign of the discriminant ∆(u) = 1 − 8M (u)N q (u), we have the following three situations.
• For ∆(u) < 0, we have no real solutions. Then, there is no physical AH.
• For ∆(u) > 0, there are two real physical AHs as r AH − (u) = r 1 (u) and r AH + (u) = r 2 (u) for the black hole in the quintessence background in contrast to the black hole in the empty space and in the dust field background having a single AH.
For the extremal case of ∆(u) = 0, there is only one single extremal AH given by which is twice of the dynamical Schwarzschild radius. Then, the dynamics of the single extremal AH (68) is governed by the relationṙ This shows that the single extremal AH (68) for the black hole in the quintessence background is moving inward during the evaporation, i.e it is shrinking. For this case, the nature of the extremal AH (68) can be determined by looking at the induced metric on it as showing that the extremal AH (68) is a timelike surface for an evaporating BH in the quintessence background. The Hawking temperature on the extremal AH (68) can be determined using the surface gravity as indicating that the extremal AH (68) is a zero temperature surface. One may interpret this zero temperature as the suppressed or zero Hawking radiation for the extremal black hole in the quintessence background. Finally, the semiclassical Bekenstein-Hawking entropy at AH (68) is given by Then, the entropy of the extremal AH (68) is four times the entropy (33) associated to the single AH (29) of the black hole in an empty space. The location of two AHs (66) and (67) for the case ∆(u) > 0 with a small quintessence parameter, i.e N q (u) ≪ M (u), are given by and It is seen that the contribution of the surrounding background quintessence field leads to (i) the existence of an inner AH (73) larger than the dynamical Schwarzschild radius r(u) = 2M (u), and (ii) a large outer (cosmological) AH (74) tending to infinity for N q (u) ≪ 1. Regarding (73) and (74), the quintessence field has a positive contribution to the inner AH while the black hole mass has a negative contribution to the outer AH. Then, as the black hole evaporation continues, the inner AH (73) shrinks to vanish while the outer AH (74) expands toward its asymptotic value of N −1 q (u). For this case, the dynamics of the AHs (73) and (74) are given bẏ andṙ respectively. Then one realizes the following points: • ForṄ q (u) = 0 representing the static surrounding quintessence field, the inner and outer AHs (73) and (74) have equal velocities but with opposite signs. Then, the inner AH (73) shrinks while the outer (74) expands during the evaporation. Also, we observe that both of these AHs are moving faster than the case of single AH (29) for the black hole in an empty space.
• ForṄ q (u) < 0, the inner AH (73) shrinks while the outer (74) expands, and these are faster than the cases in static quintessence background and in empty spaces.
The nature of the AHs (73) and (74) can be determined by looking at the induced metrics on them up to O(N q ) as and These show that forṄ q (u) ≤ 0, the inner AH (73) is a timelike surface while the outer AH (73) is an spacelike surface. ForṄ q (u) > 0, the nature of these AHs depend on the M (u), L, N q (u) andṄ q (u) values. For the specific case of 4M 2 (u)Ṅ q (u) = 2L + 8LM (u)N q (u), the inner AH (73) is null while the outer (74) is timelike. The Hawking temperature on the AHs (73) and (74) of the black hole are and Then, regarding the condition 8M (u)N q (u) < 1, one realizes that the temperature of the inner AH (73) is always positive while the temperature of the outer AH (74) is negative. For the inner AH (73), the background quintessence field decreases the temperature of the AH relative to the case of black hole in empty space. Here, regarding (64), the negative temperature of the outer AH (74) is resulting from the fact that the quintessence background field effect dominates at large distances which leads to the negative surface gravity, namely the gravitational repulsion. Finally, the semiclassical Bekenstein-Hawking entropy at the horizons (73) and (74) are given by and (82) Then, for the quintessence background with N q (u) ≪ 1, the inner AH (73) has a negligible entropy in comparison to the entropy (74) associated with the outer AH. Also regrading (81), we see that the background quintessence field contributes positively in the entropy of the inner AH (73), which leads to a higher entropy relative to the single AH (29) of the black hole in an empty background. In the limit of N q (u) → 0, the entropy (81) of the inner AH (73) tends to the entropy (33) of AH (29) of the black hole in an empty background while the entropy of the outer AH blows up to infinity, as the horizon itself goes to infinity. Also regarding (82), one realizes that since the entropy is a non-negative quantity, then it is required that the first term overcomes to all of the next negative terms. This requires the value of quintessence normalizing parameter N q (u) to be very small. In fact, this entropy consideration on the cosmological AH (74) gives a direct restriction on the quintessence field normalization parameter N q (u).

Evaporating Vaidya Black Hole Surrounded by the Cosmological Constant-Like Field
For the case of Vaidya black hole surrounded by the cosmological constant-like field with the equation of state parameter ω c = −1 [32,44], which we call it in the following as the cosmological field in short, the metric (1) takes the following form where N c (u) is the normalization parameter for the cosmological field surrounding the black hole, with the dimension of [N c (u)] = l −2 . Regarding the positive energy condition on the surrounding cosmological field, represented by the relation (6), it is required that N c (u) > 0. Regarding the metric (83), one realizes the non-trivial contribution of the characteristic feature of the surrounding cosmological field to the metric of the Vaidya black hole. In this case, the background cosmological field changes the causal structure of the original Vaidya black hole in an empty background to the Vaidya-de Sitter space [12], similar to the case of Schwarzschild to Schwarzschild-de Sitter case. Notice to the difference between our more general solution here and that of obtained in [12] having a constant cosmological field. The surface gravity associated with the metric (83) is given by Then, we see that the cosmological background appears with a negative surface gravity contribution which decreases the gravitational attraction of the black hole. This characteristic feature of the cosmological field, namely the gravitational repulsion, is evidenced by its application as the most favored candidates for the dark energy as the field responsible for the current accelerating expansion of the universe [44]. The expansion scalar for the outgoing null geodesics reads as Then, one can find the solutions to Θ = 0 determining the locations of the AHs as where ∆(u) = 1 − 27M 2 (u)N c (u). Then, depending on the sign of the discriminant ∆(u) we have the following three situations: • For ∆(u) < 0, the first solution is negative while the second and third solutions are imaginary and complex conjugate. Then, there is no physical AH.
• For the extremal case ∆(u) = 0, equivalent to N c (u) = 1 27M 2 (u) , the first solution will be negative as and the second and third solutions coincide to a single extremal AH as r 2 (u) = r 3 (u) = 1 √ 3Nc(u) = 3M (u). Then, the first solution is not a physical horizon while the coinciding second and third solutions are representing a physical extremal AH.
• For ∆(u) > 0, the first solution will be negative, and gives no physical AH while the second and third solutions give two physical inner and outer AHs.
For the extremal case of ∆(u) = 0, there is only one single extremal AH given by which is one and a half times of the dynamical Schwarzschild radius. Then, the dynamics of the single extremal AH (87) is governed by the relationṙ This shows that the single extremal AH (87) for the black hole in the cosmological background is moving inward during the evaporation, i.e it is shrinking. For this case, the nature of the extremal AH (87) can be determined by looking at the induced metric on it as indicating that the extremal AH (87) is a timelike surface for an evaporating BH in the cosmological background. The Hawking temperature on the extremal AH (87) can be determined using the surface gravity as indicating that the extremal AH (87) is a zero temperature surface. One may interpret this zero temperature as the suppressed or zero Hawking radiation for the extremal black hole in the cosmological background. Finally, the semiclassical Bekenstein-Hawking entropy at extremal AH (87) is given by Then, the entropy of the extremal AH (87) is 9 4 the entropy (33) associated to the single AH (29) of the black hole in an empty space.
For the case of ∆(u) > 0, we have the following solutions representing the inner and outer AH, respectively. It is seen that the contribution of the surrounding background cosmological field leads to (i) the existence of an inner AH (92) larger than the dynamical Schwarzschild radius r(u) = 2M (u), and (ii) a large outer (cosmological) AH (93) tending to infinity for N c (u) ≪ 1. Regarding (92) and (93), similar to the quintessence case, the cosmological field has a positive contribution to the inner AH while the black hole mass has a negative contribution to the outer AH. Then, as the black hole evaporation continues, the inner AH (92) shrinks to vanish while the outer AH (93) expands toward its asymptotic value of N − 1 2 c (u). For this case, the dynamics of the AHs (92) and (93) are given bẏ anḋ respectively. Then, one realizes the following points: • ForṄ c (u) = 0, showing the static surrounding cosmological field, the inner AH (92) shrinks while the outer AH (93) expands. Also, we observe that the inner AH (92) is moving inward faster than the case of AH (29) for the black hole in an empty space.
• ForṄ c (u) < 0, the inner AH (92) shrinks while the outer AH (93) expands, and these are faster than the cases in static cosmological field background and in empty space.
• ForṄ c (u) > 0, the behaviour of these AHs (92) and (97) These show that forṄ c (u) ≤ 0, the inner AH (92) is a timelike surface while the outer AH (93) is spacelike surface. ForṄ c (u) > 0, the nature of these AHs depends on the M (u), L, N c (u) andṄ c (u) values. For the specific case of 8M 3 (u)Ṅ c (u) = 2L + 24LM 2 (u)N c (u), the inner AH (92) is null while the nature of the outer AH (93) depends on the total value of the remaining terms in the parenthesis in (97). The Hawking temperature on the inner and outer AHs (92) and (93) of black hole can be determined using the surface gravity as and respectively. One realizes that the temperature of the inner AH (92) is always positive while the temperature of the outer AH (93) is negative. For the inner AH (92), the background cosmological field decreases the temperature of the AH relative to the AH (29) of black hole in empty space. Here, regarding (84), the negative temperature of the outer AH (93) is resulting from the fact that the cosmological background field effect dominates at large distances which leads to the negative surface gravity, equivalently to the gravitational repulsion. Finally, the semiclassical Bekenstein-Hawking entropy at the AHs (92) and (93) are given by and Then, for the cosmological field background with N c (u) ≪ 1, the inner AH (92) has a negligible entropy (100) in comparison to the entropy (101) associated with the outer AH (93). Also regrading (100), we see that the background cosmological field contributes positively to the entropy of the inner AH (92), which leads to a higher entropy relative to the single AH (29) of the black hole in an empty background. In the limit of N c (u) → 0, the entropy (100) of the inner AH (92) tends to the entropy (33) of AH (29) of the black hole in an empty background while the entropy (100) of the outer AH (93) blows up to infinity. Also regarding (101), one realizes that since the entropy is a non-negative quantity, then this relation requires that the first term overcomes to all of the next negative terms. This requires the value of cosmological field normalizing parameter N c (u) to be very small. In fact, this entropy consideration on the cosmological AH (93) gives a restriction on the cosmological field normalization parameter N c (u).

Evaporating Vaidya Black Hole Surrounded by the Phantom Field
For the case of Vaidya black hole surrounded by the phantom field with the equation of state parameter ω p = − 4 3 [32,44], the metric (1) takes the following form where N p (u) is the normalization parameter for the phantom field surrounding the black hole, with the dimension of [N p (u)] = l −3 . Regarding the positive energy condition on the surrounding phantom field, represented by the relation (6), it is required that N p (u) > 0. Regarding the metric (102), one realizes the non-trivial contribution of the characteristic feature of the surrounding phantom field to the metric of the Vaidya black hole. In this case, the background phantom filed changes the causal structure of the original Vaidya black hole in an empty background to the case with a phantom asymptotic. The surface gravity associated with the metric (102) is given by Then, we see that the phantom background appears with a negative surface gravity contribution which decreases the gravitational attraction of the black hole. This characteristic feature of the phantom field, namely the gravitational repulsion, is favored by its cosmological application as one of the candidates for the dark energy [44]. The expansion scalar for the outgoing null geodesics reads as Then, we find the solutions of Θ = 0 determining the locations of the AH as , where Γ(u) = 8 2 and ∆(u) = 1 − 2048 27 M 3 (u)N p (u). Then, depending on the sign of the discriminant ∆(u) we have the following three situations.
• For ∆(u) < 0, the first and second solutions are negative while the third and fourth solutions are imaginary. Then, there is no physical AH.
• For the extremal case ∆(u) = 0, equivalent to N p (u) = • For ∆(u) > 0, the first and second solutions are imaginary and complex conjugate to each other and gives no physical AH, while the third and fourth solutions give two physical inner and outer AH.
For the extremal case of ∆(u) = 0, there is only one single extremal AH given by • ForṄ p (u) < 0, the inner AH (112) shrinks while the outer AH (112) expands, and these are faster than the cases in static phantom field background and in empty space.
• ForṄ p (u) > 0, the behaviour of the AHs (112) and (113) depend on the M (u), L, N p (u) andṄ p (u) values. For the specific case of 16M 4 (u)Ṅ p (u) = 2L + 64LM 2 (u)N p (u), the inner AH (112) is frozen while the outer's behaviour depends on the total value of the remaining terms in (115).
Then, for the phantom background with N p (u) ≪ 1, the inner AH (112) has a negligible entropy in comparison to the entropy associated with the outer AH (113). Also regrading (120), we see that the background phantom field contributes positively in the entropy of the inner AH (112), which leads to a higher entropy relative to the single AH (29) of the black hole in an empty background. In the limit of N p (u) → 0, the entropy (120) of the inner AH (112) tends to the entropy (33) of AH (29) of the black hole in an empty background while the entropy of the outer AH (113) blows up to infinity. Also regarding (121), one realizes that since the entropy is a non-negative quantity, then this relation requires that the first term overcomes to all of the next terms. This requires the value of phantom field normalizing parameter N p (u) to be very small. In fact, this entropy consideration on the cosmological AH (113) gives a restriction on the phantom field normalization parameter N p (u).

Conclusion
In this work, we have studied the thermodynamical features and dynamical evolutions of the various apparent horizons associated with the Vaidya evaporating black hole surrounded by the cosmological fields of dust, radiation, quintessence, cosmological constant-like and phantom. We have explored in detail how do these surrounding fields contribute to the characteristic features of a surrounded dynamical black hole. To this aim, using the null vector decomposition of the spacetime metric and following York [30], we have obtained the geometric properties of the solution in its general form. We have found that the solution to g(∂ u , ∂ u ) = g uu = 0 representing the TLS is the same as the solution to Θ = 0 representing the AH location, except the trivial solution r = ∞, for the latter. Then, the AH and the TLS coincide for our spacetime (1). Also, for the surrounded evaporating black hole solutions here, we have verified that the condition dΘ du ≃ 0 determining the location of EH is equal to the vanishing acceleration of the congruences of null geodesics at EH. In the following, we summarize some of our results for the AH properties of the evaporating black hole in various backgrounds.
• For the black hole with surrounding dust field, the dust background appears with a positive surface gravity increasing the gravitational attraction of the black hole. For this case, the size of AH is larger than that of black hole in an empty space. For the black hole surrounded by a static dust field, i.ė N d (u) = 0, AH is timelike surface and its dynamics is the same as the dynamics of AH in black hole in the empty space in which it is shrinking during the evaporation. For a dynamical background witḣ N d (u) = 2L, the AH is a frozen null surface. For the background with −2L+Ṅ d (u) > 0, AH is spacelike surface and expands during the evaporation while for −2L +Ṅ d (u) < 0, it is a timelike surface and shrinks during the evaporation. The temperature of AH for the black hole in dust field is lower than the temperature of AH of the black hole in an empty space. Finally, the entropy associated with the AH in a dust field background is higher than that of the black hole in an empty space.
• For the black hole with surrounding radiation field, the surrounding radiation field contributes as the effective charge-like term for the Vaidya black hole. For this case, the radiation background appears with a negative surface gravity decreasing the gravitational attraction of the black hole, exactly the same as the charge effect in Reissner-Nordström black hole. For the extremal case of ∆(u) = 0, there is only one single AH equal to half of the dynamical Schwarzschild radius which is a timelike and zero temperature surface, and shrinks during the evaporation. The entropy of the extremal AH is a quarter of the entropy associated to the single AH of the black hole in an empty space. For ∆(u) < 0, we have no real solutions while for ∆(u) > 0, there are two real physical AHs. For the latter case, the contribution of the background radiation field leads to (i) the existence of a small inner AH and (ii) an smaller outer AH in comparison to the case of the AH for the black hole in an empty background. For the static surrounding radiation field, i.eṄ r (u) = 0, the inner AH is expanding while the outer AH is shrinking during the evaporation, faster than the case of black hole in an empty background. ForṄ r (u) > 0, the inner and outer AH are expanding and shrinking, respectively, but both moving faster than the AH of the black hole in the static radiation background. ForṄ r (u) < 0, it is possible that the inner and outer AH move again outward and inward, respectively, but both will be slower than the AH of the black hole in the radiation backgrounds withṄ r (u) = 0 andṄ r (u) > 0. In the particular case of L N r (u) = −Ṅ r (u)M (u), the inner AH will be frozen while the outer AH will move inward exactly the same as AH of the black hole in an empty background. ForṄ r (u) ≥ 0 the inner and outer AH are spacelike and timelike respectively, while forṄ r (u) < 0, we see that the nature of AH depends on the value ofṄ r (u) 2M(u) . For the specific case of L N r (u) = −Ṅ r (u)M (u), the inner AH will be null while the outer AH will be timelike. The temperature of the inner AH is always negative while the temperature of the outer AH is positive, exactly the same as in the Reissner-Nordström black hole. The inner AH has a negligible entropy relative to the entropy associated with the outer AH in a dilute radiation background. The background radiation field contributes negatively in the entropy of the outer AH leading to a lower entropy relative to the entropy of the single AH in an empty background.
• For the black hole with surrounding quintessence field, the quintessence background contributes negatively to the surface gravity which decreases the gravitational attraction of the black hole. This characteristic feature of the quintessence field is favored by its cosmological application as one of the candidates for the dark energy. For the extremal case of ∆(u) = 0, there is only one single extremal AH twice of the dynamical Schwarzschild radius which is a timelike and zero temperature surface and shrinks during the black hole evaporation. The entropy of the extremal AH is four times the entropy associated to the single AH of the black hole in an empty space. For ∆(u) < 0, we have no real solutions while for ∆(u) > 0, there are two real physical AHs. For the latter case, the contribution of the background quintessence field leads to (i) the existence of an inner AH larger than the dynamical Schwarzschild radius r(u) = 2M (u), and (ii) a large outer cosmological AH tending to infinity for N q (u) ≪ 1. For the static surrounding quintessence field, i.eṄ q (u) = 0, the inner and outer AH have equal velocities but with opposite signs. Then, the inner AH is a timelike surface and shrinks while the outer is an spacelike surface and expands during the evaporation in which both of these AH are moving faster than the case of the single AH for the black hole in an empty space. ForṄ q (u) < 0, the inner AH is a timelike surface and shrinks while the outer is an spacelike surface and expands in which these are faster than the cases in static quintessence background and in empty spaces. Foṙ N q (u) > 0, the behaviour and nature of the AH depend on the M (u), L, N q (u) andṄ q (u) values. For the specific case of 4M 2 (u)Ṅ q (u) = 2L + 8LM (u)N q (u), the inner AH is a frozen null surface while the outer AH is shrinking timelike surface. The temperature of the inner AH is always positive while the temperature of the outer AH is negative. For the inner AH, the background quintessence field decreases the temperature of the AH relative to the case of black hole in empty space. Finally, the inner AH has a negligible entropy in comparison to the entropy associated with the outer AH. Also, the background quintessence field contributes positively in the entropy of the inner AH, which leads to a higher entropy relative to the single AH of the black hole in an empty background.
• For the black hole with surrounding cosmological field, the cosmological filed contributes negatively to the surface gravity decreasing the gravitational attraction of the black hole. This characteristic feature of the cosmological field is evidenced by its application as the most favored candidates for the dark energy. For the extremal case of ∆(u) = 0, there is only one single extremal AH, 3/2 times the dynamical Schwarzschild radius, which is a timelike and zero temperature surface and shrinks during the black hole evaporation. The entropy of the extremal AH is 9 4 times the entropy associated to the single AH of the black hole in an empty space. For ∆(u) < 0, we have no real solutions while for ∆(u) > 0, there are two real physical AHs. For the latter case, the contribution of the background cosmological field leads to (i) the existence of an inner AH larger than the dynamical Schwarzschild radius r(u) = 2M (u), and (ii) a large outer cosmological AH tending to infinity for N c (u) ≪ 1. For the static surrounding cosmological field, i.eṄ c (u) = 0, the inner AH is a timelike surface and shrinks faster than the case of the single AH for the black hole in an empty space while the outer is an spacelike surface and expands during the evaporation. ForṄ c (u) < 0, the inner AH is a timelike surface and shrinks while the outer is an spacelike surface and expands in which these are faster than the cases in static cosmological background and in empty spaces. ForṄ c (u) > 0, the behaviour and nature of the AH depend on the M (u), L, N c (u) andṄ c (u) values. For the specific case of 8M 3 (u)Ṅ c (u) = 2L + 24LM 2 (u)N c (u), the inner AH is a frozen null surface while the outer AH is shrinking timelike surface. The temperature of the inner AH is always positive while the temperature of the outer AH is negative. For the inner AH, the background cosmological field decreases the temperature of the AH relative to the case of black hole in empty space. Finally, the inner AH has a negligible entropy in comparison to the entropy associated with the outer AH. Also, we find that the background cosmological field contributes positively in the entropy of the inner AH, which leads to a higher entropy relative to the single AH of the black hole in an empty background.
• For the black hole with surrounding phantom field, the phantom filed contributes negatively to the surface gravity which decreases the gravitational attraction of the black hole. This characteristic feature of the phantom field is evidenced by its application as one of the candidates for the dark energy. For the extremal case of ∆(u) = 0, there is only one single extremal AH, 4/3 times the dynamical Schwarzschild radius, which is a timelike and zero temperature surface and shrinks during the black hole evaporation. The entropy of the extremal AH is 16 9 times the entropy associated to the single AH of the black hole in an empty space. For ∆(u) < 0, we have no real solutions while for ∆(u) > 0, there are two real physical AHs. For the latter case, the contribution of the background phantom field leads to (i) the existence of an inner AH larger than the dynamical Schwarzschild radius r(u) = 2M (u), and (ii) a large outer cosmological AH tending to infinity for N p (u) ≪ 1. For the static surrounding phantom field, i.eṄ p (u) = 0, the inner AH is a timelike surface and shrinks faster than the case of the single AH for the black hole in an empty space while the outer is an spacelike surface and expands during the evaporation. ForṄ p (u) < 0, the inner AH is a timelike surface and shrinks while the outer is an spacelike surface and expands in which these are faster than the cases in static cosmological background and in empty spaces. ForṄ p (u) > 0, the behaviour and nature of the AHs depend on the M (u), L, N p (u) andṄ p (u) values. For the specific case of 16M 4 (u)Ṅ p (u) = 2L + 64LM 2 (u)N p (u), the inner AH is a frozen null surface while the outer AH is shrinking timelike surface. The temperature of the inner AH is always positive while the temperature of the outer AH is negative. For the inner AH, the background phantom field decreases the temperature of the AH relative to the case of black hole in empty space. Finally, the inner AH has a negligible entropy in comparison to the entropy associated with the outer AH. Also, we find that the background phantom field contributes positively in the entropy of the inner AH, which leads to a higher entropy relative to the single AH of the black hole in an empty background.
We aim to report elsewhere the detailed study of EH properties together with the back-reaction and metric fluctuation for the surrounded Vaidya black holes.