{'fraction_non_alphanumeric': 0.05007213483936803, 'fraction_numerical': 0.01945881276610718, 'mean_word_length': 4.024752475247524, 'pattern_counts': {'":': 0, '<': 47, '<?xml version=': 0, '>': 23, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The Accretion of Dark Energy onto a Black Hole

The stationary, spherically symmetric accretion of dark energy onto a Schwarzschild black hole is considered in terms of relativistic hydrodynamics. The approximation of an ideal fluid is used to model the dark energy. General expressions are derived for the accretion rate of an ideal fluid with an arbitrary equation of state p=p(\rho) onto a black hole. The black hole mass was found to decrease for the accretion of phantom energy. The accretion process is studied in detail for two dark energy models that admit an analytical solution: a model with a linear equation of state, p=\alpha(\rho-\rho_0), and a Chaplygin gas. For one of the special cases of a linear equation of state, an analytical expression is derived for the accretion rate of dark energy onto a moving and rotating black hole. The masses of all black holes are shown to approach zero in cosmological models with phantom energy in which the Big Rip scenario is realized.


INTRODUCTION
In recent years, strong observational evidence that the Universe is currently expanding with acceleration has been obtained. In the Einstein theory of gravitation, this positive acceleration is explained by the dominance of dark energy with a negative pressure in the Universe [1][2][3][4]. Several theoretical models of dark energy have been suggested: the vacuum energy (the cosmological constant Λ) or such dynamical components as quintessence [5][6][7][8][9][10] and k-essence [11][12][13]. Models with dynamical dark energy seem more realistic, since tracker [14,15], or attractor, solutions are realized in them. Thus, the problem of fine tuning the parameters of the Universe is solved [11][12][13].
A peculiar property of cosmological models with dark energy is the possibility of a Big Rip [16,17]: an infinite increase in the scale factor of the Universe in a finite time. The Big Rip scenario is realized in the case of dark energy, the so-called phantom energy (for which ρ + p < 0). In the Big Rip scenario, the cosmological phantom energy density tends to infinity, and all of the bound objects are torn apart up to subnuclear scales. It should be noted, however, that the condition ρ + p < 0 alone is not enough for the Big Rip scenario to be realized [18]. Alam et al. [19] analyzed data on distant supernovas in a model-independent way and showed that the presence of phantom energy with −1.2 < w < −1 in the Universe at present is highly likely. The quantum properties of the phantom energy in curved spacetime were considered in [20]. The entropy of the Universe filled with phantom energy was discussed in [21]. Models with phantom energy are also used to construct mole burrows [22,23]. The accretion of a scalar field onto a black hole from special poten-tials V (φ) was considered in [24][25][26][27][28][29]. We use a different approach to describe the accretion of dark energy onto a black hole; more specifically, we model the dark energy by an ideal fluid with a negative pressure.
In our recent paper [30] (see also [31]), we showed that the masses of all black holes in the Universe with phantom energy gradually decrease, and the black holes disappear completely by the Big Rip. In this paper, we consider in detail the stationary spherical accretion of dynamical dark energy onto a black hole. The dark energy is modeled by an ideal fluid with a negative pressure. The history of research on the accre-tion of an ideal fluid onto a compact object begins with Bondi's classic paper [32]. A relativistic generalization was made by Michel [33] (see also [34][35][36][37][38][39][40][41] for further generalizations and supplements to Michel's solution). Carr and Hawking [42] considered the accretion of dust and radiation onto a black hole by solving the complete system of Einstein equations and taking into account the back reaction of the surrounding matter (see also [43] for a description of the progress made in this area and for a discussion of fundamental questions). Below, we obtain a solution for the stationary accretion of a test relativistic ideal fluid with an arbitrary equation of state p(ρ) onto a Schwarzschild black hole. Using this solution, we show that the black hole mass decreases during the accretion of phantom energy. The masses of black holes can decrease during accretion in the case of phan-tom energy due to the violation of the energy domi-nance condition (ρ + p ≥ 0) that underlies the theorem on the nondecreasing area of the event horizon of a classical black hole [44].
This paper is structured as follows. In Section 2, we derive general equations for the spherical accretion of an ideal fluid and describe basic parameters of the steady energy flux onto a black hole. We consider an arbitrary equation of state, w = p/ρ, where the pressure p can be positive (for ordinary matter) and negative (for dark energy, including phantom energy w < −1). Note that the parameter w of the equation of state need not be constant in our approach. Accretion causes the black hole mass to change: the mass increases for ρ + p > 0 and decreases for ρ + p < 0. The energy flux turns out to be completely determined by the black hole mass M , the dark energy density at infinity ρ ∞ , and the equation of state p = p(ρ) only if 0 < ∂p/∂ρ < 1. In this case, there is a critical point that fixes the flux just as for an ordi-nary fluid. When the condition 0 < ∂p/∂ρ < 1 is violated, the dark energy flux onto a black hole can formally be arbitrary. For 0 < ∂p/∂ρ < 1, we describe the method of calculating the fluid parameters at the critical point and the energy flux onto a black hole for given M , ρ ∞ , and p = p(ρ). In Section 3, we consider specific models of the equation of state for dark energy. In the first model, we use a simple equation of state with a linear density dependence of the pressure. We consider the special cases of accretion of several types of ideal fluid: ther-mal radiation, matter with an ultrahard equation of state, dark energy with ∂p/∂ρ ≶ 0, and linear phantom energy. The accretion rate of dark energy onto a moving black hole was calculated for the special case of ∂p/∂ρ = 1. As the second model, we investigate the accretion of a Chaplygin gas onto a black hole. The evolution of the black hole mass in the Universe with the Big Rip is considered in Section 4. The possibility that the presence of phantom energy will lead the Universe to the Big Rip in the future has been discussed in recent years. The problem of the fate of black holes in this Universe is solved in a rather unexpected way: black holes are not torn apart, but disappear by the Big Rip due to the accretion of phantom energy, irrespective of their initial masses. In Section 5, we discuss the cor-respondence between the accretion of dark energy modeled by an ideal fluid onto a black hole and the accretion of a scalar field. The results obtained are briefly discussed in Section 6.

GENERAL EQUATIONS FOR SPHERICAL ACCRE-TION
Let us consider the stationary, spherically symmetric accretion of an ideal fluid that models the dark energy in the special case of a negative pressure onto a black hole. The dark energy density is assumed to be low enough for the metric to be a Schwarzschild one with a high accuracy: where M is the mass of the black hole, r is the radial coordinate, and θ and φ are the angular spherical coordinates. We model the dark energy by an ideal fluid with the energy-momentum tensor where ρ is the density, p is the dark energy pressure, and u µ = dx µ /ds is the radial 4-velocity component. The pressure is assumed to be an arbitrary function of the density, p = p(ρ). Integrating the zeroth (time) compo-nent of the conservation law T µν ;ν = 0 yields the first integral of motion for stationary, spherically symmetric accretion (Bernoulli's relativistic equation or the energy equation): where x = r/M , u = dr/ds, and C 1 is the constant determined below. To find the second integral of motion, we use the equation for the component of the energy-momentum tensor conservation law along the 4-velocity In our case, this quation is [45] u µ ρ ,µ + (ρ + p)u µ ;µ = 0.
For the given equation of state p = p(ρ), the auxiliary function n can be defined by the relation The function n is identical to the particle concentration for an atomic gas, but it can also be used to describe a continuous medium that does not consist of any particles. In this case, the 'concentration' n is a formal auxiliary function. For an arbitrary equation of state p = p(ρ), we obtain a solution for n from Eq. (7): Using (8), we find the sought second integral of motion from Eq. (5): where n ∞ (the dark energy 'concentration' at infinity) was introduced for convenience. In the case of a fluid flow directed toward the black hole, u = (dr/ds) < 0, and therefore, the numerical constant A > 0. From (3) and (9), we can easily obtain where Let us now calculate the radial 4-velocity component and the fluid density on the event horizon of the black hole, r = 2M . Setting x = 2, we obtain from Eqs. (9), (10) and (11) where ρ H is the density on the x = 2 horizon. Thus, having specified the density at infinity ρ ∞ , the equation of state p = p(ρ), and the flux A and using definition (7) of the concentration, we can calculate the fluid density ρ H on the event horizon of the black hole from (12). Given the density on the horizon ρ H , we can easily determine the radial fluid 4-velocity component on the horizon from (9): Below, we will see that the constant A, which defines the energy flux onto the black hole, can be calculated for hydrodynamically stable ideal fluids with ∂p/∂ρ > 0. This can be done by determining the fluid parameters at the critical point. Following Michel [33], we find the relationship between the parameters at the critical point: where Together with (7), this yields where c 2 s = ∂p/∂ρ is the square of the effective speed of sound in the medium. We derive the following relation for the critical point from Eqs. (14), (16), (11), and (10): which fixes the fluid density at the critical point ρ * for an arbitrary equation of state p = p(ρ). Specifying ρ * and using (8), we can determine n(ρ * ). Accordingly, the quantities x * and u * can be calculated from (14) and (16). As a result, the numerical constant A can be calculated by substituting the derived quantities into (9). For c 2 s < 0 or c 2 s > 1, no critical point exists beyond the event horizon of the black hole (x * > 1), implying that the dark energy flux onto the black hole depends on the initial conditions for an unstable ideal fluid c 2 s < 0 or a 'superluminal' fluid c 2 s > 1. This result has a simple physical interpretation: the accreted fluid has a critical point if its speed increases from subsonic to supersonic values as it approaches the black hole. In contrast, for c 2 s or c 2 s > 1, the critical point either does not exist or is formally within the event horizon of the black hole. It should also be noted that fluids with c 2 s < 0 are hydrodynamically unstable (see [46,47] for a discussion).
Equation (10), together with (6), (8), and (9), defines the accretion rate onto a black hole. These equations are valid for an ideal fluid with an arbitrary equation of state p = p(ρ), in particular, for a gas of massless particles (thermalized radiation) and a gas of massive particles. For a gas of massive particles, the system of equations (9) and (10) reduces to a similar system of equations found by Michel [33]. It should be noted, however, that Eqs. (6), (8), (9), and (10) are also valid for dark energy, including phantom energy with ρ + p < 0. In these cases, the concentration n(ρ) is positive for any ρ, while the constant C 2 in Eq. (10) is negative.
The rate of change in the black hole mass (the energy flux onto the black hole) through accretion isṀ = −4πr 2 T r 0 . Using (9) and (10), this expression can be rewritten as [30] It follows from Eq. (18) than the mass of the black hole increases as it accretes the gas of particles when p > 0, but decreases as it accretes the phantom energy when p + ρ < 0. In particular, this implies that the black hole masses in the Universe filled with phantom energy must decrease. This result is general in nature. It does not depend on the specific form of the equation of state p = p(ρ); only the satisfaction of the condition p + ρ < 0 is important. The physical cause of the decrease in the black hole mass is as follows: the phantom energy falls to the black hole, but the energy flux associated with this fall is directed away from the black hole. If we ignore the cosmological evolution of the density ρ ∞ , then we find the law of change in the black hole mass from (18) to be where M i is the initial mass of the black hole, and τ is the evolution time scale:

Model of a Linear Equation of State
Let us consider the model of dark energy with a linear density dependence of pressure [30]: where α and ρ 0 0 are constants. Among the other cases, this model describes an ultrarelativistic gas (p = ρ/3), a gas with an ultrahard equation of state (p = ρ), and the simplest model of dark energy (ρ 0 = 0 and α < 0). The quantity α is related to the parameter w = p/ρ of the equation of state by w = α(ρ − ρ 0 )/ρ. An equation of state with w = const < 0 throughout the cosmological evolution is commonly used to analyze cosmological models. The matter with such an equation of state is hydrodynamically unstable and can exist only for a short period. Our equation of state (21) for α > 0 does not have this shortcoming. For α > 0, it also allows the case of hydrodynamically stable phantom energy to be described, which is not possible when using an equation of state with w = const < −1. In the real Universe, the equation of state changes with time (i.e., w depends on t). Therefore, Eq. (21) has the meaning of an approximation to the true equation of state only in a limited ρ range. From the physical point of view, the condition ρ > 0 must be satisfied for any equation of state in a comoving frame of reference. In particular, the state of matter with ρ = 0, but p = 0, is physically unacceptable. The corresponding constraints for the equation of state (21) are specified by conditions (29) and (30) given below.
For α < 0, there is no critical point for the accreted fluid flow. For α > 0, using (14) and (16), we obtain the parameters of the critical point Note that the parameters of the critical point (22) in the linear model (21) are determined only by ∂p/∂ρ = α and do not depend on ρ 0 , which fixes the physical nature of the fluid under consideration: a relativistic gas, dark energy, or phantom energy. Note also that no critical point exists beyond the event horizon of the black hole for α > 1 (this corresponds to a nonphysical situation with a superluminal speed of sound). Let us calculate the constant A, which defines the energy flux onto the black hole. We find from (8) where we introduced the effective density Using (17), we obtain ρ eff * ρ eff,∞ α/(1+α) where ρ eff * and ρ eff,∞ , are the effective densities at the critical point and at infinity, respectively. Substituting (24) into (23) and using (9), we obtain for the linear model It is easy to see that A ≥ 4 for 0 < α < 1. A = 4 for α = 1 (this corresponds to c s = 1); i.e., the constant A is on the order of 1 for relativistic speeds of sound. Using (25), we obtain from (20) To determine the fluid density on the event horizon of the black hole, we substitute (23) into (12) to yield where A is given by (25). For 0 < α < 1, the effective density on the horizon ρ eff,H cannot be lower than ρ eff,∞ . The radial 4-velocity component on the horizon can be found from (13) and (27): The value of u H changes from 1 to 1/2 for 0 < α < 1. The linear model (21) describes the phantom energy when ρ ∞ /ρ 0 < α/(1 + α). In this case, ρ + p < 0. However, the requirement that the density ρ be nonnegative should be taken into account. This parameter can formally be negative in the range 0 < α ≤ 1. Such a nonphysical situation imposes a constraint on the linear model (21) under consideration. For a physically proper description of the accretion process, we must require that the density ρ be nonnegative. We obtain the following constraint on the validity range of the linear model from (27) for hydrodynamically stable phantom energy: As follows from (29), at a given α, we can always choose the parameters ρ 0 and ρ ∞ in such a way that ρ > 0 for any r > 2M . On the other hand, model (21) describes the quintessence (not the phantom energy) for the entire r range only if p < 0. Consequently, a physically proper description of the quintessence can be obtained from (27) if For some of the specific choices of α (more specifically, for α = 1/3, 1/2, 2/3, and 1), ρ(x) and u(x) can be calculated analytically (see the Appendix for details on these calculations). In Figs. 1 and 2, the radial 4-velocity component, u, and the density normalized to the density at infinity, ρ/ρ ∞ are plotted against the coordinate x = r/2M . Figure 1 (left) shows the plots of u(x) for hydrodynamically stable fluids with c 2 s > 0 at α = 1/3, 1/2, 2/3, and 1 (the curves are arranged from top to bottom, respectively). Figure  1 (right) shows the plots of u(r) for c 2 s < 0 at α = −1.1, −2, and −1/2 (the curves are also arranged from top to bottom). For this case, we chose the boundary condition u H = 1 on the horizon. Figure 2 (left) shows the plots of ρ/ρ ∞ for a hydrodynamically stable fluid with α = 1 for various cases: ρ 0 = 0 (the model of neutron star matter); ρ 0 /ρ ∞ = 16/9 (the linear model of nonphantom dark energy); ρ 0 /ρ ∞ = 7/3 (the linear model of phantom energy); and ρ 0 /ρ ∞ = 7/3 (the linear model of phantom energy with ρ H = 0) (the curves are arranged from top to bottom, respectively). Figure 2 (right) shows the plots of (ρ/ρ ∞ ) for various cases: α = −2, ρ 0 = 0, and A = 4 (the linear model of phantom energy, the upper curve); and α = −1/2, ρ 0 = 0, and A = 4 (the linear model of nonphantom energy, the lower curve). For this case, we chose the velocity |u H | = 1 on the horizon.

Accretion onto a Moving and Rotating Black Hole
Let us consider the accretion onto a moving and rotating black hole in the special case of a linear equation of state with α = 1. The condition α = 1 allows an exact analytical expression to be derived for the accre-tion rate of dark energy onto a black hole.
For α = 1, we easily find from (23) that We obtain the following continuity equation for the particle concentration from (5): We can introduce the scalar field φ in terms of which the fluid velocity can be expressed as follows (there is no torsion in the fluid): We derive an equation for the auxiliary function φ by using Eqs. (31) and (32), Exactly the same equation arises in the problem of the accretion of a fluid with the equation of state p = ρ [41]. Thus, we reduced the problem of a black hole moving in dark energy with the equation of state p = ρ − ρ 0 to the problem of a fluid with an extremely hard equation of state, p = ρ. Using the method suggested in [41], we obtain the mass evolution law for a moving and rotating black hole immersed in dark energy with the equation of state p = ρ − ρ 0 : where is the radius of the event horizon for a rotating black hole, a = J/M is the specific angular momentum of the black hole (rotation parameter), and u 0 ∞ is the zeroth 4-velocity component of the black hole relative to the fluid. Expression (34) for u 0 BH = 0 reduces to (18) for a Schwarzschild (a = 0) black hole at rest.

Chaplygin Gas
Let us consider a Chaplygin gas with the following equation of state as another example of the solvable model: where α > 0. The range of parameters ρ 2 < α represents the phantom energy with a superluminal speed of sound, implying that the phantom energy flux onto the black hole is not fixed by the condition of its passage through the critical point. The case of ρ 2 > α corresponds to dark energy with ρ + p > 0 and 0 < c 2 s < 1. We can easily find from Eq. (8) that The density at the critical point can be calculated from (17) and (36): The velocity and the radial coordinate at the critical point are given by We then find the constant A from Eq. (9): For 0 < c 2 s < 1, the constant A cannot be smaller than 4, as in the case of the linear model. The evolution time scale of the black hole mass without any cosmological change in the dark energy density is given by Note that Eqs. (36)-(39) are applicable only for dark energy with ρ + p > 0 and are invalid for phantom energy. On the black hole horizon, For 0 < c 2 s < 1, the density on the horizon ρ H cannot be lower than ρ ∞ , and u H changes from 1 to 1/2. The Chaplygin gas density distribution can be determined from the general equations (9) and (10): The velocity distribution u(r) can be calculated by using Eqs. (9), (36), and (42). In Figs. 3 and 4, the velocity u and the density normalized to the density at infinity, ρ/ρ ∞ , are plotted against the coordinate x = r/2M . Figure 3 (left) shows the plots of u(r) for nonphantom dark energy at ρ 2 ∞ /α = 3, 2, and 1.1 (the curves are arranged from top to bottom, respectively). Figure 3 (right) shows the plots of u(r) for phantom energy at ρ 2 ∞ /α = 0.3, 0.5, and 0.9 (the curves are also arranged from top to bottom, respectively). In this case, the boundary condition u H = 1 is set on the horizon. Figure 4 (left) shows the plots of the normalized density, ρ/ρ ∞ , for nonphantom dark energy at ρ 2 ∞ /α = 3, 2, and 1.1 (the curves are arranged from top to bottom, respectively). Figure 4 (right) shows the plots of the normalized density, (ρ/ρ ∞ ), for phantom energy at ρ 2 ∞ /α = 0.9, 0.5, and 0.3 (the curves are also arranged from top to bottom, respectively). For this case, we chose the boundary condition u H = 1 on the horizon.

THE FATE OF BLACK HOLES DURING THE BIG RIP
Let us now consider the evolution of black holes in the cosmological Big Rip scenario, where the scale factor a(t) increases to infinity in a finite time [16,17]. For simplicity, we take into account only the dark energy and disregard the other forms of energy. In the linear model (21), the Big Rip takes place at ρ + p < 0 and α < −1. The following relation can be derived from Friedmann's equations in the case of a linear equation of state: |ρ + p| ∝ a −3(1+α) .
Setting, for simplicity, ρ 0 = 0, we find the evolution law of the phantom energy density in this Universe: where Here, ρ ∞,i is the initial cosmological phantom energy, and the initial time was chosen in such a way that the Big Rip occurs at time τ . We easily see from Eqs. (20) and (43) that the Big Rip takes place at α ≡ ∂p/∂ρ < −1. In general, the condition ρ + p < 0 alone is not enough for the cosmological evolution to be ended with the Big Rip [18]. Using Eq. (43), we find the evolution of the black hole mass in the cosmological Big Rip scenario from Eq. (18): and M i is the initial mass of the black hole. At α = −2 and a typical value of A = 4 (which correspond to u H = −1),Ṁ 0 = 3/8. In the limit t → τ (i.e., near the Big Rip), the t dependence of the black hole mass becomes linear, M ≃Ṁ 0 (τ − t). When t approaches τ , the rate of decrease in the black hole mass ceases to depend on the initial black hole mass and the phantom energy density: In other words, the masses of all black holes near the Big Rip are approximately equal and approach zero. This implies that the accretion of phantom energy dominates over the Hawking evaporation until the black hole mass decreases to the Planck mass. Formally, however, all black holes in the Universe completely evapo-rate during the Hawking radiation in the Planck time before the Big Rip occurs.

THE ACCRETION OF A SCALAR FIELD
In this section, we compare our calculations of the accretion of an ideal fluid with similar calculations of the accretion of a scalar (nonphantom) field onto a black hole [24][25][26][27][28][29]. The dark energy is commonly modeled by a scalar field with a potential V (φ). The approximation of an ideal fluid is rougher, since the scalar field φ and ∂ µ φ cannot be unambiguously reproduced for given ρ and p, which characterize an ideal fluid. Despite this difference between the scalar field and the ideal fluid, we will show that our results are in close agreement with the corresponding calculations of the accretion of a scalar field onto a black hole.
The Lagrangian of the scalar field is L = K − V , where K is the kinetic term and V is the potential. For the standard choice of the kinetic term K = φ ;µ φ ;µ /2, the corresponding energy flux onto the black hole is Jacobson [24] found a solution for the scalar field in the Schwarzschild metric for a zero potential, where φ ∞ is the scalar field at infinity. Frolov and Kofman [26] showed that this solution is also valid for many scalar fields with a nonzero potential V() under certain conditions. For this solution and, accordingly,Ṁ = 4π(2M ) 2φ2 ∞ . The energy-momentum tensor constructed using Jacobson's solution is identical to the energy-momentum tensor for an ideal fluid with an extremely hard equation of state, p = ρ, after the substitution p ∞ →φ 2 ∞ /2, ρ ∞ →φ 2 ∞ /2. This is not surprising, since the theory of a scalar field with a zero potential, V (φ), is identical to the model of an ideal fluid [48]. In view of this correspondence, we easily see agreement between our result (18) forṀ in the case of p = ρ and the corresponding results from [24,26].
The Lagrangian of the scalar field that describes the phantom energy must have a negative kinetic term [16,17], for example, K = −φ ;µ φ ;µ /2 (see [49] for more general cases). In this case, the phantom energy flux onto the black hole has the opposite sign, where φ is the solution of the same Klein-CGordon equation as that for the standard scalar field, but with the substitution V → −V . For a zero potential, this solution is identical to Jacobson's solution [24] obtained for a scalar field with a positive kinetic term. However, the Lagrangian with a negative kinetic term and V (φ) = 0 does not describe the phantom energy. At the same time, the solution for a scalar field with V (φ) = 0 is identical to the solution for a positive constant potential, V 0 = const, which can be chosen in such a way that ρ = −φ 2 /2 + V 0 > 0.
In this case, the scalar field describes the accreted phantom energy with ρ > 0 and p < −ρ, which leads to a decrease in the black hole mass at the ratė A simple example of phantom cosmology (but without the Big Rip) is realized by a scalar field with the potential V = m 2 φ 2 /2, where m ∼ 10 −33 eV [50]. After a short transition period, this cosmological model approaches an asymptotic state with H ≃ mφ/3 1/2 ,φ ≃ 2m/3 1/2 .

DISCUSSION AND CONCLUSIONS
In recent years, the concept of dark energy has been accepted and extensively discussed in cosmology. The possible existence of dark energy with a negative pressure leads to new cosmological scenarios, including the exotic model of the Universe in which all of the bound objects are destroyed and which dies itself as a result of the Big Rip. To determine the fate of black holes in this cosmological scenario, we considered the spherically symmetric, stationary accretion of dark energy mod-eled by an ideal fluid onto a black hole. We derived general equations for the accretion of an ideal fluid with the equation of state p = p(ρ) onto a Schwarzschild black hole. In particular, these equations can be used to describe the accretion of thermal radiation, dark energy, and phantom energy. We also considered the accretion onto a moving and rotating black hole in the special case of an extremely hard equation of state, p = ρ. We calculated the change in the black hole mass through accretion. The black hole masses for ρ+p > 0 were found to increase, as in the usual case. However, a qualitatively new result was obtained for phantom energy, i.e., for a medium with ρ + p < 0. We found that the black hole masses decrease in this situation. Using this result, we solved the problem of the fate of black holes in a universe that undergoes the Big Rip. It turns out that all black holes in this Universe must decrease their masses and disappear completely by the Big Rip. We also considered the correspondence between the accretion of dark energy in the model of an ideal fluid and the accretion of a scalar field.