-
The universe is in a stage of accelerated expansion, which can be described by a negative pressure system [1]. Two sources of negative pressure are considered to drive this expansion: the cosmological constant and quintessence dark energy. Inspired by the thermodynamics of anti-de Sitter (AdS) black holes (BHs) against the background of the variable cosmological constant Λ, Kastor
etal. and Dolan proposed a scenario to investigate BH thermodynamics in the extended phase space. If Λ is given the meaning of thermodynamic pressure P, the BH mass is naturally explained by thermodynamic enthalpy [2–4].Based on this approach, various BH phase transitions have been advanced to AdS space-time. The phase transition of small-large BHs occurs in charged AdS BH systems, and it is similar to the van der Waals (vdW) fluid system in classical thermodynamics. In Ref. [5], Kubiznak
etal. calculated the universal constant and critical exponents of the Reissner-Nordström (RN) AdS BH and vdW systems, observing that these key parameters are consistent. Weietal. provided a universal formula of normalized scalar curvatureRN in the Ruppeiner geometry approach suitable for classical and BH systems and observed that the specific values of the critical exponent ofRN of RN AdS BHs and vdW fluids are identical [6, 7]. Hence, the difference in the critical behavior between the BH and vdW systems can be used to prove the deviation of the BH phase transition from the vdW phase transition. Some relevant thermodynamic studies on different charged BHs in extended phase space are reported in Refs. [8–14].Considering RN AdS BH coupling with quintessence dark energy, the phase transition characteristics of BH become interesting. Li showed the effect of quintessence dark energy on the critical behavior of RN AdS BHs [15]. Compared with the standard vdW phase transition, the universal constant has a significant difference and noticeably changes with the values of BH parameters, but the critical exponents retain the standard value. Toledo
etal. investigated the thermal behaviors of RN AdS BHs with cloud of strings and quintessence dark energy, observing that the critical point exhibited an analogous phenomenon. [16]. We believe that the reason for this phenomenon is that the thermal behavior of quintessence has not been fully considered.As a candidate to explain the accelerated expansion of the universe, quintessence dark energy can provide a pressure similar to the cosmological constant. The pressure
Pq is related to the state parameterωq and the energy densityρq of quintessence, satisfyingPq≡ωqρq [17, 18]. Based on the assumption that the energy tensor of the perfect fluid satisfies the linearity and additivity, Kiselev derived the solution of a spherically-symmetric BH in theωq range of quintessence dark energy (−1<ωq<−1/3 ) [19]. Because Kiselev's model is simple and feasible, most thermodynamic studies based on the coexistence of dark energy and BHs introduce quintessence using this approach, including reseach related to critical behavior [15, 16, 20–22].In contrast to the research considering Λ as the phase transition pressure of AdS BHs, we abandon the discussion of BH phase transitions in AdS space-time. We focus on the thermal characteristics of the interaction between quintessence dark energy and BHs. One of the reasons for this consideration is that combining the thermal behavior of quintessence with BH thermodynamics is easy. Using the practice of constructing the extended phase space with the
Λ−P relation, we can naturally construct the "quintessential" phase space using the state equation of quintessence.The paper is organized as follows. In Sec. II, we briefly review the solution of a BH surrounded by dark energy quintessence. In Sec. III, we construct the BH thermodynamics in the quintessential phase space and investigate the
Pq−υ critical behavior. Sec. IV provides our conclusions and discussion. -
With the condition of additivity and linearity in the energy momentum tensor, Kiselev derived the quintessential matter field equation and obtained a static spherically-symmetric charged BH solution [19]. We briefly review this work in this section. The general metric of the spherically symmetric static gravitational field is
ds2=eνdt2−eλdr2−r2(dθ2+sin2θdϕ2),
(1) where ν and λ are functions of r. The Einstein equations are expressed as
2Ttt=−e−λ(1r2−λ′r)+1r2,
(2) 2Trr=−e−λ(1r2+ν′r)+1r2,
(3) 2Tθθ=2Tϕϕ=−12e−λ(ν″+ν′22+ν′−λ′r−ν′λ′2).
(4) The energy-momentum tensor of the quintessence can be expressed as
Ttt=Trr=ρq,Tθθ=Tϕϕ=−12ρq(3ωq+1),
(5) where
ρq andωq are the energy density and state parameters of quintessence, respectively. By usingλ=−ln(1+f) , the energy-momentum tensor of the quintessence can be rewritten asTtt=Trr=−12r2(f+rf′),
(6) Tθθ=Tϕϕ=−14r(2f′+rf″).
(7) Using Eqs. (5) and (7), we can observe that f satisfies the linear equation
(3ωq+1)f+3(ωq+1)rf′+r2f″=0,
(8) with two solutions in the following form:
fq=cr3ωq+1,fBH=−rgr,
(9) where c and
rg are the normalization factors. Because the additivity and linearity are considered to be satisfied,fBH−q=fBH+fq is also considered a solution of the above equation. We can observe that there are two coexisting forms between the quintessence and BHs: non interacting form (fq andfBH ) and interacting form (fBH−q ). We can predict that the non interacting form contains a scenario similar to the cosmological constant, that is, whenfq has a larger scale evolution characteristic thanfBH . Our paper focuses on the interacting formfBH−q . Substituting the first solution of quintessence into Eq. (5), we obtainρq asρq=c23ωqr3(ωq+1).
(10) To ensure that
ρq is always greater than zero, c should be consistent withωq in sign. It is well known that dark energy is divided into three types according to the scope ofωq :● Case 1 asymptotically flat:
−1/3<ωq<0 .● Case 2 quintessence dark energy:
−1<ωq<−1/3 .● Case 3 phantom dark energy:
ωq<−1 .The metric of exact spherically-symmetric solutions for the Einstein equations describing BHs interacting with the quintessential matter is expressed as
ds2=f(r)dt2−f(r)−1dr2−r2(dθ2+sin2θdϕ2),
(11) where
f(r) is the metric potential, which can be expressed as [19]f(r)=1−2Mr−∑n(rnr)3ωn+1,
(12) where M is the BH mass,
rn is the dimensional normalization constant, andωn is the quintessential state parameter. From this expression, we can obtain the metric for RN BH interacting with quintessence or the Schwarzschild one if the BH chargeQ=0 ; we obtainf(r)q=1−2Mr+Q2r2−(rqr)3ωq+1.
(13) -
Assuming
r3ωq+1q≡kc (k>0) , Eq. (13) can be rewritten asf(r)q=1−2Mr+Q2r2−kcr3ωq+1.
(14) The mass M can be expressed with the horizon radius
r+ asM=r+2+Q22r+−kc2r3ωq+.
(15) The Hawking temperature can be obtained as
T=f′(r+)q4π=14πr+(1−Q2r2++3ωqkcr3ωq+1+).
(16) The entropy is
S=∫r+01T(∂M∂r+)dr+=πr2+.
(17) With the newly thermodynamic quantities defined from quintessence, the conventional first law of BH thermodynamics is
dM=TdS+ΦdQ+Cdc,
(18) where
C is the physical quantity conjugate to the parameter c. The quantity related to quintessence introduced in this manner can make the first law consistent with the Smarr relation. The physical quantities from Eq. (15) are given byΦ=(∂M∂Q)S,c=Qr+,
(19) C=(∂M∂c)S,Q=−k2r3ωq+.
(20) Additionally, the conventional Smarr relation is
M=2TS+ΦQ+(1+3ωq)Cc.
(21) It is possible that a significant component of the energy density of the universe has an equation of state different from those of matter, radiation, or cosmological constant (Λ), i.e., quintessence dark energy [23]. It is a dynamical, slowly-evolving, spatially inhomogeneous component with negative pressure [24]. Based on this concept, we assume that the interaction between the BH and quintessence causes the BH to obtain pressure, which is related to the pressure of quintessence. Using Eq. (10) and restoring the discarded
4π , the equation of state of quintessence with ideal fluid properties isPq=c8π3ω2qr3(ωq+1).
(22) Note that for theoretical analysis, we do not limit
ωq to the range of dark energy(−1,−1/3) ; therefore, the sign of c is not required, and the pressurePq in the quintessential phase space can be any real number. Hence, the following relationship can be easily derived:dc=8πr3(ωq+1)3ω2qdPq.
(23) Using the above equation, we can popularize the conventional first law of BH thermodynamics (18) and conventional Smarr relation (21) to the quintessential phase space:
dM=TdS+ΦdQ+VqdPq,
(24) M=2TS+ΦQ+(1+3ωq)VqPq,
(25) where
Vq is the conjugate quantity of the pressurePq , satisfyingVq=8πr3(ωq+1)3ω2qC=4πkr3+3ω2q.
(26) Note that
Vq can return to the general thermodynamic volumeV=43πr3+ if the free parameterk=ω2q . Furthermore, the Smarr relation in the quintessential phase space can return to the extended phase space ifωq=−1 . -
In the previous subsection, we use the equation of state of quintessence to construct the quintessential phase space. The first law of BH thermodynamics and Smarr relation in this scenario are derived. Next, we investigate
Pq−υ criticality under the framework of the quintessential phase space.Utilizing Eqs. (16) and (22), the equation of state of a BH can be expressed as
Pq=ωqk(T2r+−18πr2++Q28πr4+).
(27) According to the critical conditions
(∂P∂r+)=(∂2P∂r2+)=0 , the critical thermodynamic quantities of the RN BH in this scenario can be obtained:rqc=√6Q,Tqc=13√6πQ,Pqc=ωqk196πQ2.
(28) The sign consistency between
ωq and c makesPq andPqc have such characteristics, that is, the occurrence of this phase transition does not depend on whether the pressure provided by quintessence is positive or negative. For dark energy (−1<ωq<−1/3 ), this phase transition can occur during the accelerating expansion of the universe caused by negative pressure provided by quintessence. The universal constantε≡Pcυc/Tc can be expressed asεq=PqcυqcTqc=ωqk38.
(29) where
υqc is the specific volume, satisfyingυqc=2rqc [5]. We set the absolute value of universal constant of a standard vdW-type phase transition as3/8 . Judging from the universal constant, the ratioωq/k can make the phase transition far from the standard vdW-type. Using reduced parameter space, we investigate the phase transition characteristics near the critical point. The equation of state Eq. (27) is restructured aspq=38τν−2ν2+13ν4,
(30) where
pq=PqPqc,τ=TTqc,ν=υυqc.
(31) We observe that Eq. (30) and Ref. [7] (Eq. (93)) are identical, indicating that the phase transition of RN BH in the quintessential phase space has the same critical exponents and critical behavior of normalized scalar curvature
RN as the standard vdW-type phase transition, and more importantly, it is independent ofωq and k.Furthermore, the Gibbs free energy of an RN BH in the quintessential phase space can be expressed as
G=M−TS=3Q24r++r+4−(2+3ωq)pqr3+144Q2ωq.
(32) We can observe that the value of the free parameter k does not affect the behavior of the Gibbs free energy. For dark energy (
−1<ωq<−1/3 ), T as a function ofr+ and G as a function of T are depicted in Fig. 1. The typical swallowtail behavior of Gibbs function is retained only whenωq is less than−2/3 value, which corresponds to the vdW-type phase transition inT−r+ plane. The atypical swallowtail behavior (the swallowtail surrounded byT=0 ) always exists, which corresponds to a distinctive process in theT−r+ plane. The distinctive process is caused by the repulsive force provided by quintessence with negative pressure as the leading force in the later stage, which is exhibited as the temperature being depressed. It differs significantly from the scenario in which gravity dominates in the late stage of the extended phase space. -
When the free parameter k is fixed to
ω2q , the thermodynamic volume is consistent with that in the extended phase space, implying that the pressure of the BH system in this case can be directly composed of these two parts. In this scenario, the equation of state in the quintessential phase space isPq=1ωq(T2r+−18πr2++Q28πr4+).
(33) In Ref. [5], the equation of state in the extended phase space is expressed as
PΛ=T2r+−18πr2++Q28πr4+.
(34) Hence, we can obtain the equation of state after the combination:
Ptot=(1−1ωq)(T2r+−18πr2++Q28πr4+).
(35) Therefore, the critical thermodynamic quantity is
rc=√6Q,Tc=13√6πQ,Ptotc=196πQ2(1−1ωq).
(36) That is, we can obtain the critical universal constant as
εq=38(1−1ωq).
(37) If the dark energy case
−1<ωq<−1/3 ,εtot is within3/4<εtot<3/2 , which remains greater than the standard value3/8 of vdW-type phase transition. Based on Eqs. (26) and (35),Ptot as a function of V is depicted in Fig. 2. The increase inωq noticeably moves the isotherm upward, that is, the system pressure is increased. The features of vdW-type phase transition can be observed on thePtot−V plane with differentωq values, implying that, in the context of AdS space-time, gravity is the leading force in the later stage, which is different from the case only quintessence dark energy exists.Figure 2. (color online)
Pq as a function of V. The BH charge is taken asQ=1 and the quintessence state parameterωq=−1,−5/6, −2/3,−1/2,−1/3 .We can easily observe that the reduced equation of state has the same form as Eq. (27); therefore, this phase transition is still a second-order phase transition.
-
When the free parameter k is fixed to
1 , Eq. (16) becomes the general form of coupling quintessence, in which the direct coupling between non-relativistic matter and quintessence is not considered. The equation of state in the quintessential phase space can be expressed asPq=ωq(T2r+−18πr2++Q28πr4+).
(38) When the general thermodynamic volume V is used as the variable, we adjust the pressure in the quintessential phase when combining the pressure of the two phase spaces, which does not break the first law of thermodynamics and Smarr relation. Using Eqs. (26) and (38), we can obtain
VqPq=VP′q→P′q=VqPqV=kω2qPq,
(39) and the combined pressure
Ptot isPtot=PΛ−kω2qPq.
(40) The equation of state after combination can yield an expression similar to Eq. (35):
Ptot=PΛ−1ω2qPq=(1−1ωq)(T2r+−18πr2++Q28πr4+).
(41) We observe that the combination of the quintessential and extended phase spaces is not regulated by the free parameter k, which means that coupling between non-relativistic matter and quintessence does not affect the thermodynamic characteristics of the RN BH system.
-
The thermodynamic properties of an RN BH in the quintessential phase space are revealed in this analysis. When RN BH is surrounded by the thermal quintessence, the quintessence exerts a pressure
Pq the event horizon in spherically symmetric BH, which leads to the discussion of the BH thermodynamics in the quintessential phase space. For an RN BH interacting with quintessence, we use the equation of the state of quintessence to obtain the first law of BH thermodynamics and Smarr relation in the quintessential phase space and investigate the critical behavior of this thermodynamic system. We observe that the conjugate quantityVq of the pressurePq is often not equal to the thermodynamic volume V. A phase transition occurs during the evolution of an RN BH in the quintessential phase space, and its similarity to the standard vdW-type phase transition depends on the ratioωq/k in terms of the universal constant. From the perspective of critical exponents and critical behavior of normalized scalar curvatureRN , the phase transition is typically a second-order phase transition, which is the same as the standard vdW-type phase transition. We believe that the existence of this phase transition results in strange critical phenomena mentioned in the introduction when an RN AdS BH surrounded by quintessence dark energy is considered in the extended phase space.The
T−r+ andG−T diagrams show that, if−1<ωq<−2/3 , the features of vdW-type phase transition emerge, including typical equal area and typical swallowtail behaviors. If−2/3<ωq<−1/3 , the features of vdW-type phase transition completely lost, a distinctive phase transition occupied and an atypical swallowtail behavior (the swallowtail surrounded byT=0 ) replaces the typical swallowtail behavior, which is caused by the depression of quintessence with negative pressure. This effect increases with the increase in the state parameterωq .We further discuss the combination of RN BHs in the quintessential and extended phase spaces under two simple cases. Constructing the combination based on the general thermodynamic volume
V=(4/3)πr3 , the form of the equation of state has consistency and is not controlled by the free parameter k. We observe that the universal constantεtot is within3/4<εtot<3/2 ifωq is within−1<ωq<−2/3 . The interaction with quintessence dark energy causes the universal constant of RN AdS BHs to be larger than the standard vdW-type phase transition, but the phase transition grade is still of second-order. These results suggest that the quintessential phase space must be constructed with the pressure from quintessence to improve the thermodynamic research of AdS BHs involving quintessence dark energy, and it can be used as a new probe to explore BH thermodynamics.
