Phase transition and entropic force of de Sitter black hole in massive gravity

It is well known that de Sitter(dS) black holes generally have a black hole horizon and a cosmological horizon, both of which have Hawking radiation. But the radiation temperature of the two horizons is generally different, so dS black holes do not meet the requirements of thermal equilibrium stability, which brings certain difficulties to the study of the thermodynamic characteristics of black holes. In this paper, dS black hole is regarded as a thermodynamic system, and the effective thermodynamic quantities of the system are obtained. The influence of various state parameters on the effective thermodynamic quantities in the massive gravity space-time is discussed. The condition of the phase transition of the de Sitter black hole in massive gravity space-time is given. We consider that the total entropy of the dS black hole is the sum of the corresponding entropy of the two horizons plus an extra term from the correlation of the two horizons. By comparing the entropic force of interaction between black hole horizon and the cosmological horizon with Lennard-Jones force between two particles, we find that the change rule of entropic force between the two system is surprisingly the same. The research will help us to explore the real reason of accelerating expansion of the universe.


I. INTRODUCTION
In recent years, the study of thermodynamic properties of black holes has aroused great interest. Especially, after treating the cosmological constant as the thermodynamic pressure, the first law of black hole thermodynamics in extended phase space has been derived. Comparing the state parameters of black holes with that of the Van der Waals equation, the critical phenomena of various AdS space-time black holes have been studied. The influence of various parameters on the phase transition is discussed. Some progress has been made in the study of black hole phase transition . For de Sitter space-time, space-time has black hole horizon and cosmological horizon. Both horizons have Hawking radiation, and the radiation temperatures of the two horizons are generally different. Therefore dS black holes do not meet the requirements of thermal stability, which brings certain difficulties to the study of the thermodynamic properties of black hole. In recent years, the study of the thermodynamic properties of de Sitter space-time has attracted extensive attention [31][32][33][34][35][36][37][38][39][40][41][42][43][44][45][46][47][48][49]. Because in the early period of inflation, our universe is a quasi de Sitter space-time, and the constant term introduced in the study of de Sitter space-time is the contribution of vacuum energy, which is also a form of material energy. If the cosmological constant corresponds to dark energy, our universe will evolve into a new de Sitter phase. In order to construct the whole history of universe evolution, we should have a clear understanding of the classical and quantum properties of the de Sitter space-time.
Since the thermodynamic quantities corresponding to the two horizons in de Sitter space-time are all functions of mass M , charge Q and cosmological constant Λ. There must be certain relation between the thermodynamic quantities corresponding to the two horizons. Considering whether the corresponding thermodynamic quantity in the de Sitter space-time has the thermodynamic characteristics of AdS black hole after the correlation of black hole horizon and universe horizon; What is the relationship between the entropy of de Sitter space and the corresponding entropy of two horizons? What is the relationship between the entropic force caused by the interaction of two horizons and the ratio of the position of the two horizons? These problems are very important to study the stability and evolution of de Sitter space-time. Therefore, it is worth to establish a self-consistent relationship between the thermodynamic quantities of de Sitter space-time. In [50], the authors suggest a massive gravity which overcomes the traditional problems and yield an avenue for addressing the cosmological constant naturalness problem. Fierz and Pauli firstly provided the possibility of a massive graviton [51,52]. Further, van Dam and Veltman [53] and Zakharov [54] had found the linear theory which was coupled to a source, they found that even when the mass of the graviton approaches to zero, the prediction of the theory is different from that of the linear theory. Later, Boulware-Deser ghost theories had been studied, which were nonlinear and ghostlike instability massive gravity theories [55,56]. In the absence of such instability, significant progress has been made in establishing the theory of massive gravity [57][58][59][60][61][62]. The most straightforward way to construct the massive gravity theories is to simply add a mass term to the GR action, giving the graviton a mass in such a way that GR is recovered as mass vanishes. Recently, a charged BTZ black holes has been studied in [12,63], which consider a massive gravity. The massive BTZ black holes in the presence of Maxwell and Born-Infeld electrodynamics in asymptotically (A)dS spacetimes was studied in [6]. In [64][65][66][67], thermodynamics and entropy of AdS black hole in massive gravity were obtained, and the phase transition was discussed.
In this paper, on the basis of effective thermodynamic quantity of de Sitter black hole in massive gravity (DSBHMG), the thermodynamic characteristics of de Sitter space-time are discussed. The condition of the phase transition of DSBHMG space-time is found, and the influence of each parameter on the critical point is analyzed. It is noted that DSBHMG space-time has a second-order phase transition similar to that AdS black hole when the parameters meet certain conditions, except that the heat capacity at constant volume of DSBHMG is not zero. By calculating the entropic force of the interaction between the two horizons, the changing law of entropic force with the ratio of the position of two horizons is given. By comparing the curve with that of Lennard-Jones force [68] between two particles, it is found that the curve of entropic force changing with the position ratio of two horizons is similar to the curve of Lennard-Jones force changing with the distance between two particles, which has the same change rule. This discovery provides a new way to study the internal cause of accelerating expansion of the universe. This paper is arranged as follows: for the continuity of this paper, in the second part, we briefly introduce the corresponding thermodynamic quantities of black hole horizon and cosmological horizon in DSBHMG and the effective thermodynamic quantities of DSBHMG. In the third part, we discuss the critical phenomena of DSBHMG. In the fourth part, we discuss the entropic force of interaction between the two horizons, and obtain the entropic force and the Lennard-Jones force between two particles has a similar law of change. The fifth part of the paper discusses the conclusion. (We use the units G = h = k B = 1).

II. THERMODYNAMICS OF BLACK HOLES IN MASSIVE GRAVITY
We consider (3 + 1) -dimensional massive gravity with a Maxwell field [65-67, 69, 70], the action is as follows where Λ is the cosmological constant and k = 0,1, or −1 , correspond to a Ricci flat, sphere, or hyperbolic horizon for the black hole, respectively. Here F µν is the Maxwell field-strength tensor, c i are constants, and µ i are symmetric polynomials of the eigenvalues of the matrix g µν f µν where f µν is a fixed symmetric tensor. The action admits a static black hole solution with the space-time metric and reference metric as where h ij dx i dx j is the line element for an Einstein space with constant curvature, f (r) is the metric function, which is written by [71,72].
The positions of black hole horizon and universe horizon satisfy the equation f (r +,c ) = 0, thus, the mass m 0 can be expressed in terms of r +,c as where x = r+ rc and The radiation temperature of black hole horizon and cosmological horizon is The thermodynamic quantities corresponding to the two horizons satisfy the first law of thermodynamics where when the radiation temperature T + of the black hole horizon is equal to that of the cosmological horizon T c , the charge Q and Λ of the system satisfies the equation we can get the radiation temperature T when T + = T c .
Recently, through the study of the thermodynamic characteristics of dS space-time, one can get the thermodynamic volume of dS space-time is [33,48] if the energy M , charge Q and volume V in space-time are taken as the state parameters of a thermodynamic system, the first law of thermodynamics should be satisfied Considering the interaction between the black hole horizon and the cosmological horizon, we set the entropy of the system [47][48][49] here the undefined function f (x) represents the extra contribution from the correlations of the two horizons. When the radiation temperature of black hole horizon and cosmological horizon is equal, the effective temperature of space-time should also be equal to the radiation temperature of two horizons From equations (12), (13) and (14), the effective temperature T ef f , pressure P ef f and potential φ ef f of the system are obtained , where Considering the initial conditions f (0) = 0, we can get we set the initial parameters to m = 2.12, c 1 = 2,c 2 = 3.18,q = 1.7, r c = 1,k = 1, and then take different values for k, q, m, c 1 and c 2 respectively, to get the curve P ef f −

III. CRITICAL PHENOMENA
Heat capacity of the system at constant volume: , and B(x, q) is presented by Eq.(16). From FIG.3, we know that the heat capacity C V of the system in DSBHMG is not zero, which is different from the result for AdS black holes where C V = 0 .
The heat capacity at constant pressure is Here 0.01 0.02 .
The α − x curve is shown in FIG.5.  The isothermal compressibility is  It can be seen from the curve C P ef f − x, α − x, κ T ef f − x , C P ef f , α and κ T ef f are divergence at x = x c . The entropy and volume of the system are continuous. According to Ehrenfest's classification of phase transition, x c is the second-order phase transition point of the system. Gibbs function is given by the G − T ef f curve is following From Figures C P ef f − x, α − x, κ T ef f and G − T ef f , DSBHMG satisfies Ehrenfest's second-order phase transition conditions for phase transition classification at x = x c , so secondary phase transition occurs in DSBHMG space-time at x = x c .
From the expressions (11) and (13) of the total entropy and volume of space-time, it is known that in the range 0 < x < 1 the entropy and volume of DSBHMG are continuous. According to the curve G − T ef f , DSBHMG does not have a first-order phase transition, which is different from the first-order phase transition of AdS black holes.

IV. THE ENTROPIC FORCE OF INTERACTION BETWEEN TWO HORIZONS
For a general thermodynamic system, when two systems are independent of each other, the corresponding entropy of system A and system B are S A and S B . When two systems have interactions, the entropy of the total system is In the formula, the S AB is an extra entropy caused by the interaction of two systems. From formula (13), we know that the total entropy of the effective thermodynamic system in DSBHMG is divided into two parts, one is the entropy corresponding to the two horizons, and the other term is the increase of the system entropy after considering the interaction of the two horizons as a thermodynamic system, so the entropy is Entropic force in thermodynamic system is [73][74][75][76][77][78] F = −T ∂S ∂r (30) which T is the temperature of the system, r is the location of the boundary surface. The entropic force in DSBHMG is Where T ef f is the equivalent temperature of the system, r = r c − r + = r c (1 − x), from Eq.(29), we can obtain Eq.(32), The F (x) − x curve is shown in FIG.8. For different parameters, the intersection position of the F (x) − x curve and the axis is x 0 . The entropic force between the two horizons is positive(F > 0) when 1 > x > x 0 , which means that the two horizons are mutually exclusive; and the entropic force between the two horizons is negative (F < 0) when x 0 > x > 0, indicating that the two horizons are mutually attractive. Therefore, the position ratio of the two horizons x is different, and the entropic force between them is different. The two horizons are accelerated to separate under the action of entropic force, and the cosmological horizon expands faster than the black hole horizon when 1 > x > x 0 , and they are decelerating to be separated by the entropic force, and the cosmological horizon undergoes a decelerating expansion relative to the black hole horizon when x 0 > x > 0. So we get that the different values S 1 − s 2 between the area S 1 of the curve and the axis in interval 1 > x > x 0 and the area S 2 of the curve and the axis in interval x 0 > x > 0 determines whether the cosmological horizon accelerates or oscillates with respect to the black hole horizon. The cosmological horizon changes from accelerating expansion to decelerating expansion when S 1 − S 2 > 0. When it is equal to or less than zero, the cosmological horizon is changing from accelerating expansion to decelerating expansion, and then changing from accelerating contraction to decelerating contraction. One cycle ends and the next cycle begins. The cosmological horizon is oscillating relative to the horizon of the black hole.
Interactions between neutral molecules or atoms with a center of mass separation r are often approximated by the so-called Lennard-Jones potential energy φ L,J which is given by [68,[77][78][79] φ L,J = 4φ min r 0 r where the first term is a short-range repulsive interaction and the second term is a longer-range attractive interaction. By the definition of potential energy, the force between a molecule and a neighbor in the radial direction from the first molecule is F r = −dφL,J dr , which is positive (repulsive) for r < r min and negative (attractive) for r > r min . When the center of the first particle coincides with the coordinate center, let the radius of the particle be r0 2 . When the center of the second particle is at r, and the boundary of the second particle is at r 2 , then r = r 2 + r0 2 . Let's take y = r0 2r2 , 0 < y ≤ 1. Eq.(33) can be expressed as φ L,J (y) = 4φ min r 0 r 12 − r 0 r 6 = 4φ min 2 6 2 6 y 1 + y The interaction between the two particles is From FIG.8 and FIG.10, we can see the relation curve of entropic force between two horizons in DSBHMG with the ratio of the position of the two horizons is very similar to that of the Lenard-Jones interaction force with the position ratio of the two particle boundaries, which has the same change rule.
In order to show the similarity of the two curves more clearly, we show the two curves in the same coordinate as is shown in FIG.11. From FIG.11, we know that the trend of change in curve F (x) − x and F (y) − y is exactly the same, only the intersection position of the two curves and the abscissa is slightly different. It shows that the entropic force between the horizon of black hole and the cosmological horizon is very similar to the Lenard Jones interaction force between two particles in RN dssq space-time, which provides a new way to explore the interaction force between particles in black hole.

V. CONCLUSION AND DISCUSSION
When the black hole horizon and the cosmological horizon are viewed separately as independent thermodynamic systems without considering the correlation between them, the space-time does not meet the requirements of thermodynamic equilibrium stability because the radiation temperature of the two horizons are different. Thus the space-time is unstable. When the correlation of two horizons is considered, the effective temperature T ef f , effective pressure P ef f and effective potential φ ef f reflecting the thermal properties of DSBHMG are given by Eq. (15). From curve C p ef f − x, α − x and κ T ef f − x, when the ratio of the positions of the two horizons in DSBHMG is x = x c , phase transition occurs in the system. According to Ehrenfest's classification of phase transition, the phase transition of the system at this point is of second-order, which is similar to that of AdS black hole. In the study of the thermodynamic properties of spherically symmetric AdS black holes, people study the critical phenomena of AdS black holes by comparing the thermodynamic quantities of AdS black holes with the thermodynamic quantities of Van der Waals equation and obtain the critical exponents and the heat capacity of AdS black holes, which provides a basis for further study of the thermal effect of black holes and experimental observation. However, it is difficult to accept the fact that the heat capacity at constant volume of AdS black holes is zero. It is shown that the effective thermodynamic quantity of DSBHMG also has the phase transition characteristics similar to that of van der Waals system. From formula (19), the heat capacity of DSBHMG at constant volume is not zero, which is consistent with the result that this heat capacity of van der Waals system is not zero. From the curve of C P ef f − x, we know that the stability of DSBHMG depends on the value of x . When x < x c , space-time meets the requirements of thermodynamic stability. When 1 > x > x c , space-time does not meet the requirements of thermodynamic equilibrium stability, and space-time is unstable. Therefore, there is no DSBHMG steady-state space-time satisfying the position ratio 1 > x > x c of two horizons in the universe, which provides a theoretical basis to find black holes. The influence of parameters in DSBHMG on space-time stability is shown in C P ef f − x curve seeing FIG.4. In the framework of general relativity, the entropic force of interaction between the horizon of black hole and the cosmological horizon given by theory is very similar to the Lennard-Jones force between two particles confirmed by experiments. Since the space-time metric is derived from relativity, the thermodynamics of space-time is based on general relativity, and the thermal effect of black hole is obtained from quantum mechanics. The physical quantities obtained meet the first law of thermodynamics. Therefore, the conclusion we have given reveals the internal relationship among general relativity, quantum mechanics and thermodynamics, and provide a new way for us to study the interaction between particles in black holes and the micro state of particles in black holes, and the relationship between Lennard-Jones potential and micro state of particles in ordinary thermodynamic systems.
Through the analysis of the fourth part, it is known that under the action of entropic force, when the ratio of the positions of the two horizons satisfies x 0 < x < 1, the entropic force between the two horizons is mutually exclusive, and the two horizons in this region are accelerating expansion under the action of entropic force; when 0 < x < x 0 , the entropic force between the two horizons is mutually attractive, and the two horizons are decelerating expansion; until x → 0, the entropic force between the two horizons is zero, and the two horizons tend to be relatively static. Therefore, in the 0 < x < 1 range, the cosmological horizon accelerates to expand relative to the horizon of the black hole under the action of the entropic force. Our universe is a quasi de Sitter space-time, so our discussion provides a new way to explore the internal cause of the expansion of the universe.