# Black Hole Solution of Einstein-Born-Infeld-Yang-Mills Theory

Authors: Kun Meng; Da-Bao Yang; Zhan-Ning Hu
Journal: Advances in High Energy Physics
 (2017)

---

## Abstract

A new four-dimensional black hole solution of Einstein-Born-Infeld-Yang-Mills theory is constructed; several degenerated forms of the black hole solution are presented. The related thermodynamical quantities are calculated, with which the first law of thermodynamics is checked to be satisfied. Identifying the cosmological constant as pressure of the system, the phase transition behaviors of the black hole in the extended phase space are studied.


---

## Body

## 1. Introduction

When considering electromagnetic field as the matter source to construct black hole solution in Einstein gravity, one usually use the standard Maxwell’s theory ofU
(
1
) gauge field. However, when the field is strong enough the linear Maxwell’s theory becomes invalid; the nonlinearities of electromagnetic field should be introduced as the correct electrodynamics. In 1930s, Born and Infeld proposed nonlinear electrodynamics with the aim of obtaining a finite value of the self-energy of electron [1]. Fradkin and Tseytlin reproduced Born-Infeld (BI) action in the framework of string theory [2]. The D3-brane dynamics was also noticed to be governed by BI action [3]. Hoffmann first found a solution of Einstein gravity coupled to BI electromagnetic field [4], which is devoid of essential singularity at the origin. Many black hole solutions of Einstein-Born-Infeld (EBI) theory with or without a cosmological constant have been found [5–8]. Black hole solutions of Gauss-Bonnet and three-order Lovelock gravity coupled to BI electromagnetic field have also been found in [9, 10].Since the condition for no-hair theorem in asymptotic (A)dS spacetime is much relaxed, besides the usually used matter sources such as scalar and Maxwell fields, it is natural to couple Yang-Mills (YM) field to gravity and construct black hole solutions. However, the YM equations are so complicated that early attempts of searching for black hole solutions of Einstein-Yang-Mills (EYM) theory were performed numerically. The first analytic black hole solution was found by Yasskin applying Wu-Yang ansatz [11]. This kind of solutions have also been generalized to higher dimensions and higher derivative gravity [12, 13]. Regular nonminimal magnetic black hole solutions have been investigated recently in [14]. To our knowledge, black hole solutions of Einstein gravity coupled to both BI electromagnetic field and YM field have not been studied; in this paper, we want to construct black hole solution of Einstein-Born-Infeld-Yang-Mills (EBIYM) gravity and study some properties of the black hole.Among the properties of black holes, thermodynamical properties have gained much attention in the past several decades. Since the early work [15] studying thermal phase transition between Schwarzchild-AdS black hole and AdS vacuum, the behaviors of phase transition of black holes have been studied intensively. In [16, 17], the thermodynamics of charged AdS black hole have been studied and found that there exist phase transitions between large/small black holes at the sacrifice of viewing the electric charge as intensive variable; the phase transition behavior resembles the one of van der Waals liquid-gas system. In order to solve the problem encountered in [16, 17], Kubizňák et al. extended the phase space by identifying the cosmological constant as the pressure of the gravitational system [18, 19]; then the thermodynamic volume conjugate to pressure can be defined; in this framework the phase transitions between large/small black holes are founded too. Subsequently, the study was extended to various other kinds of black holes [20–23]. In this paper, we will also explore the phase transitions of the black hole given below in extended phase space.The paper is organized as follows. In Section2, we present the black hole solution of four-dimensional Einstein-BI-YM gravity and discuss its several degenerated forms. In Section 3, we calculate the thermodynamical quantities, check the first law of thermodynamics, and study the phase transitions of the black hole in extended phase space. Finally, some concluding remarks are given in Section 4.


## 2. Black Hole Solution of Einstein-BI-YM Theory

We intend to search for a four-dimensional black hole solution of EBIYM theory; the action of the theory is given by(1)
I
=
∫
d
4
x
-
g
R
-
2
Λ
16
π
G
+
L
F
-
Tr
⁡
F
μ
ν
a
F
a
μ
ν
,where G is the Newton constant; we take the convention 16
π
G
=
1 in the following for convenience, Tr
⁡
(
·
) means sum over the gauge group indices, and(2)
L
F
=
4
β
2
1
-
1
+
F
μ
ν
F
μ
ν
2
β
2
.Here β is the BI parameter with dimension of mass. In the limit β
→
0, L
(
F
)
→
0. In the limit β
→
∞, L
(
F
) reduces to the standard Maxwell form:(3)
L
F
=
-
F
μ
ν
F
μ
ν
+
O
F
4
.Take variation with respect to the metric, the Maxwell fieldA
μ, and the YM field A
μ
(
a
), respectively, we have the equations of motion(4)
G
μ
ν
+
Λ
g
μ
ν
=
T
μ
ν
M
+
T
μ
ν
YM
,
(5)
∇
μ
F
μ
ν
1
+
F
ρ
σ
F
ρ
σ
/
2
β
2
=
0
,
(6)
∇
μ
F
a
μ
ν
+
f
b
c
a
A
μ
b
F
c
μ
ν
=
0
,where f
(
b
)
(
c
)
(
a
) are the real structure constants of the gauge group and(7)
T
μ
ν
EM
=
1
2
g
μ
ν
L
F
+
2
F
μ
ρ
F
ν
ρ
1
+
F
ρ
σ
F
ρ
σ
/
2
β
2
,
T
μ
ν
YM
=
-
1
2
g
μ
ν
F
ρ
σ
a
F
a
ρ
σ
+
2
F
μ
ρ
a
F
ν
a
ρ
.In order to find a static solution of the field equations, we take a general metric ansatz(8)
d
s
2
=
-
f
r
d
t
2
+
d
r
2
f
r
+
r
2
d
Ω
2
2
,where d
Ω
2
2 is the line element of 2-sphere. Before solving (4), we need to solve equations of motion of the matter fields (5) and (6) first. Under the electrostatic potential assumption, all other components of the strength tensor F
μ
ν vanish except F
t
r; solving (5) one obtains(9)
F
t
r
=
β
q
β
2
r
4
+
q
2
,where q is an integration constant. Note that unlike the field strength in Maxwell electrodynamics, the one in BI electrodynamics is finite at r
=
0. To solve equations of motion of the YM field (6), we take the magnetic Wu-Yang ansatz of the gauge potential [14]. Since we want to find a black hole solution with global YM charge, the gauge group is supposed to be S
U
(
2
); then the real structure constants f
(
a
)
(
b
)
(
c
) are the complete antisymmetric symbols ε
(
a
)
(
b
)
(
c
). We use the position dependent generators t
(
r
)
,
t
(
θ
)
,
t
(
φ
) of the gauge group; the relations between the generators t
(
r
)
,
t
(
θ
)
,
t
(
φ
) and the standard S
U
(
2
) generators are(10)
t
r
=
cos
⁡
ν
φ
sin
⁡
θ
t
1
+
sin
⁡
ν
φ
sin
⁡
θ
t
2
+
cos
⁡
θ
t
3
,
t
θ
=
cos
⁡
ν
φ
cos
⁡
θ
t
1
+
sin
⁡
ν
φ
cos
⁡
θ
t
2
-
sin
⁡
θ
t
3
,
t
φ
=
-
sin
⁡
ν
φ
t
1
+
cos
⁡
ν
φ
t
2
,where ν is a nonvanishing integer; as will become clear later, it relates to global YM charge. Using the commutation relations [
t
(
a
)
,
t
(
b
)
]
=
ε
(
a
)
(
b
)
(
c
)
t
(
c
), where a, b, and c take values from 1 to 3 and ε
(
1
)
(
2
)
(
3
)
=
1, it is straightforward to check the following commutation relations:(11)
t
r
,
t
θ
=
t
φ
,
t
θ
,
t
φ
=
t
r
,
t
φ
,
t
r
=
t
θ
,which are satisfied. The gauge field characterized by Wu-Yang ansatz is of the form (12)
A
0
a
=
0
,
A
r
a
=
0
,
A
θ
a
=
-
δ
φ
a
,
A
φ
a
=
ν
sin
⁡
θ
δ
θ
a
.One can check that the equations of motion of YM field (6) are satisfied under the choice of gauge (12). Under the gauge (12), the only nonvanishing component of the YM field strength tensor is(13)
F
θ
φ
=
-
ν
sin
⁡
θ
t
r
.With the field strength tensors (9), (13), and the metric ansatz (8), one can work out the expression of energy momentum tensor straightforwardly. Substituting (7) and (8) into (4) and solving (4) we find the following new black hole solution:(14)
f
r
=
1
-
m
0
r
+
1
3
2
β
2
-
Λ
r
2
+
ν
2
r
2
-
2
β
3
q
2
+
r
4
β
2
+
4
q
2
3
r
2
F
1
2
1
2
,
1
4
,
5
4
,
-
q
2
r
4
β
2
.Here m
0 is the mass parameter and F
1
2
[
1
/
2
,
1
/
4
,
5
/
4
,
-
q
2
/
r
4
β
2
] is the Gaussian hypergeometric function. Before further analysis of solution (14), we discuss its several limiting forms. When the parameter ν in (14) vanishes, the Einstein-BI black hole in (A)dS space is reproduced [8]. When β
→
∞, BI electrodynamics reduces to the standard Maxwell form; solution (14) reduces to(15)
f
r
=
1
-
m
0
r
-
1
3
Λ
r
2
+
q
2
r
2
+
ν
2
r
2
,which is a new black hole with both electric charge and magnetic YM charge. When β
→
0, the electric charge vanishes automatically; one obtains a magnetic black hole with only YM charge(16)
f
r
=
1
-
m
0
r
-
1
3
Λ
r
2
+
ν
2
r
2
,which is of the same form as the 4-dimensional (A)dS-Reinssner-Nordström (RN) black hole. When ν
=
0 and β
→
∞, solution (14) degenerates to RN black hole in (A)dS space.


## 3. Thermodynamics

In this section, we will first calculate the thermodynamic quantities of the black hole obtained in the last section, then check the first law of thermodynamics, and finally extend the phase space to explore the phase transition behaviors of the black hole.


### 3.1. Calculations of Thermodynamical Quantities

Since the black hole (14) is asymptotic to (A)dS, in order to calculate the mass of the black hole, one should adopt the method proposed by Abbott and Deser [24], which gives(17)
M
=
8
π
m
0
.Note that the area of unit 2-sphere ω
2
=
4
π is used. In order to calculate the temperature, one should take the Euclidean continuation (t
→
-
i
τ) and demand the absence of conical singularity at the horizon; the period of the Euclidean time τ is the inverse of Hawking temperature 1
/
T; one obtains(18)
T
=
1
4
π
r
+
1
+
2
β
2
r
+
2
-
Λ
r
+
2
-
ν
2
r
+
2
-
2
β
q
2
+
r
+
4
β
2
,where r
+ is the radius of the outmost horizon of the black hole. In model (1) we considered that entropy of the black hole (14) obeys the area law, which is(19)
S
=
16
π
2
r
+
2
;note that the convention 16
π
G
=
1 is adopted. The electric charge is defined as(20)
Q
e
=
1
4
π
∫
∗
F
d
Ω
2
=
q
,where ∗
F represents Hodge dual of the strength tensor and the integration is performed on the t
=
const and r
→
∞ hypersurface. The electrostatic potential Φ, measured at infinity with respect to the horizon, is defined as(21)
Φ
=
A
μ
χ
μ
r
→
∞
-
A
μ
χ
μ
r
=
r
+
.Our calculation gives Φ
=
q
/
r
+
F
1
2
[
1
/
2
,
1
/
4
,
5
/
4
,
-
q
2
/
r
+
4
β
2
]. Since the gauge group is S
U
(
2
), we can define a global YM charge, for which we take the definition [25–27](22)
Q
YM
=
1
4
π
∫
F
θ
φ
a
F
θ
φ
a
d
θ
d
φ
.Just as the definition of electric charge, the integration in (22) is performed on t
=
const and r
→
∞ hypersurface too. Our calculation gives(23)
Q
YM
=
-
ν
. In order to preserve the YM charge being positive, ν is taken to be negative. Now with all the thermodynamic quantities above in hand, one can check that the first law of thermodynamics(24)
d
M
=
T
d
S
+
Φ
d
Q
e
+
U
d
Q
YM is satisfied, where U defined by U
≡
∂
M
/
∂
Q
YM
S
,
Q
e is the thermodynamic potential conjugate to the YM charge Q
YM.


### 3.2. Phase Transition in Extended Phase Space

In extended phase space, cosmological constant is identified as pressure of the gravitational system:(25)
P
=
-
Λ
8
π
. There are some reasons to do so. First, one can suppose there exist more fundamental theories, where the cosmological constant arises as vacuum expectation value. Second, the Smarr relation becomes inconsistent with the first law of thermodynamics unless the variation of Λ is included [18]. Now the mass of the black hole is viewed as enthalpy H
≡
M rather than internal energy of the system [28]; then the Gibbs free energy of the system is G
=
H
-
T
S. If the BI parameter β is viewed as a free thermodynamic variable too, one can check the first law of thermodynamics in extended phase space(26)
d
H
=
T
d
S
+
Φ
d
Q
e
+
U
d
Q
YM
+
V
d
P
+
B
d
β is satisfied, where(27)
V
=
∂
H
∂
P
S
,
β
,
Q
e
,
Q
YM
=
64
π
2
r
+
3
3is the thermodynamic volume conjugate to the pressure, and(28)
B
=
∂
H
∂
β
S
,
P
,
Q
e
,
Q
YM
=
4
π
3
8
β
r
+
3
-
8
r
+
r
+
4
β
2
+
Q
e
2
+
4
Q
e
2
β
r
+
F
1
2
1
2
,
1
4
,
5
4
,
-
Q
e
2
r
+
4
β
2 is the thermodynamic quantity conjugate to β.Before studying the phase transitions of the black hole (14) for a general β, let us first discuss two limiting cases β
→
0 and β
→
∞. When β
→
0, in terms of the temperature and horizon radius, the pressure of the system can be written as(29)
P
=
T
2
r
+
-
1
8
π
r
+
2
+
Q
YM
2
8
π
r
+
4
.The critical points of the thermal phase transition are determined by the following critical equations:(30)
∂
P
∂
r
+
r
+
=
r
c
,
T
=
T
c
=
∂
2
P
∂
r
+
2
r
+
=
r
c
,
T
=
T
c
=
0
.The critical equations (30) can be solved analytically, which gives the critical horizon radius(31)
r
c
=
6
Q
YM and the critical temperature(32)
T
c
=
1
3
6
π
Q
YM
. It is easy to see that the critical horizon radius and temperature and consequently the critical pressure are solely determined by the YM charge. This is because, as previously mentioned, when β
→
0, the electric charge vanishes automatically and the black hole carries only YM charge as shown in (16).Whenβ
→
∞, this case corresponds to the black hole described by the metric (8) with f
(
r
) given in (15). The pressure of the system is given by(33)
P
=
T
2
r
+
-
1
8
π
r
+
2
+
Q
e
2
+
Q
YM
2
8
π
r
+
4
.Solving the critical equations (30), we obtain the critical horizon radius (34)
r
c
=
6
Q
e
2
+
Q
YM
2and the critical temperature(35)
T
c
=
1
3
π
6
Q
e
2
+
Q
Y
M
2
.The critical horizon radius, temperature, and pressure are determined by both electric and YM charge. Now, for the black hole (15) by taking proper values of the related parameters we can plot the isobars as displayed in Figure 1. On T
-
r
+ plot of Figure 1, the dotted line is the isobar of P
>
P
c; in this case there is only a single phase which is in analogy to the thermal behavior of “ideal gas.” The dashed line is the isobar of P
=
P
c. The solid line is the isobar of P
<
P
c; one can see that there are two branches of black holes whose temperatures increase with the increase of horizon radius; one branch is in small radius region and the other one is in large radius region. Between the two branches, there is a region whose temperature decreases with the increase of horizon radius, which means the thermal expansion coefficient is negative; this region corresponds to an unstable phase. For appropriate value of temperature, two horizon radii are permitted; thus large/small black holes phase transition occurs. This kind of thermal behavior resembles the one of van der Waals gas/liquid system. The phase transition can also be seen in G
-
T plot, where “swallow tail” appears when P
<
P
c, just as the solid line shows us. The dashed line on G
-
T plot is the isobar of P
=
P
c, and the dotted line corresponding to P
>
P
c describes “ideal gas” phase.Figure 1
Isobaric plots forQ
YM
=
2, Q
e
=
1. The isobars on T
-
r
+ plot correspond one-to-one with the G
-
T plot. On both (a) and (b), from top to bottom, the pressures are, respectively, P
=
1.2
P
c (dotted), P
=
P
c (dashed), and P
=
0.8
P
c (solid).
(a)
(b)Now, let us consider phase transition of the black hole (14) for a general value of β. For general β, the pressure of the system is (36)
P
=
-
1
8
π
r
+
2
+
T
2
r
+
+
β
2
4
π
1
+
Q
e
2
β
2
r
+
4
-
β
2
4
π
+
Q
YM
2
8
π
r
+
4
. In this case, the critical equations (30) are too complicated to be solved analytically, we solve them numerically, the results are collected in Table 1. From the data in the table, one learns that, with the increase of β or Q
YM, T
c and P
c decrease while r
c and the ratio P
c
r
c
/
T
c increase. Figure 2 gives the isobaric plots near the critical point of P
c
=
0.0006631 and Figure 3 displays the effects of YM charge and BI parameter on the critical pressures. On T
-
r
+ plot of Figure 2, just as the case of β
→
∞, P
>
P
c isobar describes the “ideal gas” phase. The isobar of P
<
P
c shows that the “large black hole” region and “small black hole” region are thermodynamical stable while the medium region is thermodynamical unstable, which indicates the existence of phase transition. On G
-
T plot of Figure 2, the isobars show that when P
<
P
c phase transition characterized by the “swallow tail” occurs. On Figure 3(b) one can see that, although corresponding to different orders of magnitude of β, the dotted line (β
=
1) and the dashed line (β
=
10) are almost coincident, and the deviation of the solid line (β
=
0.1) from the dotted or dashed one is small; it seems that the influence of the values of β to critical pressure is not that large just as the authors found in [23]. This may be attributed to the parameter region we choose.Table 1
Critical values forQ
e
=
1.
β
Q
YM
r
c
T
c
P
c
P
c
r
c
/
T
c
10
3
7.7459654
0.0136979
0.0003316
0.1875150
10
2
5.4772185
0.0193717
0.0006631
0.1875024
10
1
3.4640315
0.0306296
0.0016579
0.1874989
1
2
5.4765161
0.0193728
0.0006632
0.1874807
0.1
2
5.4113483
0.0194752
0.0006714
0.1865541Figure 2
Isobaric plots forQ
YM
=
2, Q
e
=
1, and β
=
10. The isobars on T
-
r
+ plot and G
-
T plot are in one-to-one correspondence with each other. On both (a) and (b), from top to bottom, the pressures are, respectively, P
=
1.2
P
c (dotted), P
=
P
c (dashed), and P
=
0.8
P
c (solid).
(a)
(b)Figure 3
Isobaric plots that reflect the effects of YM charge (a) and BI parameter (b). On (a), the parameters of all the isobars are taken asQ
e
=
1, β
=
10, and P
=
0.7
P
c; the difference is that the solid, dashed, and dotted lines correspond to Q
YM
=
1,2
,
3, respectively. On (b), the parameters of all the isobars are taken as Q
e
=
1, Q
YM
=
2, and P
=
0.7
P
c; the solid, dotted, and dashed lines correspond to β
=
0.1,1
,
10, respectively.
(a)
(b)


## 3.1. Calculations of Thermodynamical Quantities

Since the black hole (14) is asymptotic to (A)dS, in order to calculate the mass of the black hole, one should adopt the method proposed by Abbott and Deser [24], which gives(17)
M
=
8
π
m
0
.Note that the area of unit 2-sphere ω
2
=
4
π is used. In order to calculate the temperature, one should take the Euclidean continuation (t
→
-
i
τ) and demand the absence of conical singularity at the horizon; the period of the Euclidean time τ is the inverse of Hawking temperature 1
/
T; one obtains(18)
T
=
1
4
π
r
+
1
+
2
β
2
r
+
2
-
Λ
r
+
2
-
ν
2
r
+
2
-
2
β
q
2
+
r
+
4
β
2
,where r
+ is the radius of the outmost horizon of the black hole. In model (1) we considered that entropy of the black hole (14) obeys the area law, which is(19)
S
=
16
π
2
r
+
2
;note that the convention 16
π
G
=
1 is adopted. The electric charge is defined as(20)
Q
e
=
1
4
π
∫
∗
F
d
Ω
2
=
q
,where ∗
F represents Hodge dual of the strength tensor and the integration is performed on the t
=
const and r
→
∞ hypersurface. The electrostatic potential Φ, measured at infinity with respect to the horizon, is defined as(21)
Φ
=
A
μ
χ
μ
r
→
∞
-
A
μ
χ
μ
r
=
r
+
.Our calculation gives Φ
=
q
/
r
+
F
1
2
[
1
/
2
,
1
/
4
,
5
/
4
,
-
q
2
/
r
+
4
β
2
]. Since the gauge group is S
U
(
2
), we can define a global YM charge, for which we take the definition [25–27](22)
Q
YM
=
1
4
π
∫
F
θ
φ
a
F
θ
φ
a
d
θ
d
φ
.Just as the definition of electric charge, the integration in (22) is performed on t
=
const and r
→
∞ hypersurface too. Our calculation gives(23)
Q
YM
=
-
ν
. In order to preserve the YM charge being positive, ν is taken to be negative. Now with all the thermodynamic quantities above in hand, one can check that the first law of thermodynamics(24)
d
M
=
T
d
S
+
Φ
d
Q
e
+
U
d
Q
YM is satisfied, where U defined by U
≡
∂
M
/
∂
Q
YM
S
,
Q
e is the thermodynamic potential conjugate to the YM charge Q
YM.


## 3.2. Phase Transition in Extended Phase Space

In extended phase space, cosmological constant is identified as pressure of the gravitational system:(25)
P
=
-
Λ
8
π
. There are some reasons to do so. First, one can suppose there exist more fundamental theories, where the cosmological constant arises as vacuum expectation value. Second, the Smarr relation becomes inconsistent with the first law of thermodynamics unless the variation of Λ is included [18]. Now the mass of the black hole is viewed as enthalpy H
≡
M rather than internal energy of the system [28]; then the Gibbs free energy of the system is G
=
H
-
T
S. If the BI parameter β is viewed as a free thermodynamic variable too, one can check the first law of thermodynamics in extended phase space(26)
d
H
=
T
d
S
+
Φ
d
Q
e
+
U
d
Q
YM
+
V
d
P
+
B
d
β is satisfied, where(27)
V
=
∂
H
∂
P
S
,
β
,
Q
e
,
Q
YM
=
64
π
2
r
+
3
3is the thermodynamic volume conjugate to the pressure, and(28)
B
=
∂
H
∂
β
S
,
P
,
Q
e
,
Q
YM
=
4
π
3
8
β
r
+
3
-
8
r
+
r
+
4
β
2
+
Q
e
2
+
4
Q
e
2
β
r
+
F
1
2
1
2
,
1
4
,
5
4
,
-
Q
e
2
r
+
4
β
2 is the thermodynamic quantity conjugate to β.Before studying the phase transitions of the black hole (14) for a general β, let us first discuss two limiting cases β
→
0 and β
→
∞. When β
→
0, in terms of the temperature and horizon radius, the pressure of the system can be written as(29)
P
=
T
2
r
+
-
1
8
π
r
+
2
+
Q
YM
2
8
π
r
+
4
.The critical points of the thermal phase transition are determined by the following critical equations:(30)
∂
P
∂
r
+
r
+
=
r
c
,
T
=
T
c
=
∂
2
P
∂
r
+
2
r
+
=
r
c
,
T
=
T
c
=
0
.The critical equations (30) can be solved analytically, which gives the critical horizon radius(31)
r
c
=
6
Q
YM and the critical temperature(32)
T
c
=
1
3
6
π
Q
YM
. It is easy to see that the critical horizon radius and temperature and consequently the critical pressure are solely determined by the YM charge. This is because, as previously mentioned, when β
→
0, the electric charge vanishes automatically and the black hole carries only YM charge as shown in (16).Whenβ
→
∞, this case corresponds to the black hole described by the metric (8) with f
(
r
) given in (15). The pressure of the system is given by(33)
P
=
T
2
r
+
-
1
8
π
r
+
2
+
Q
e
2
+
Q
YM
2
8
π
r
+
4
.Solving the critical equations (30), we obtain the critical horizon radius (34)
r
c
=
6
Q
e
2
+
Q
YM
2and the critical temperature(35)
T
c
=
1
3
π
6
Q
e
2
+
Q
Y
M
2
.The critical horizon radius, temperature, and pressure are determined by both electric and YM charge. Now, for the black hole (15) by taking proper values of the related parameters we can plot the isobars as displayed in Figure 1. On T
-
r
+ plot of Figure 1, the dotted line is the isobar of P
>
P
c; in this case there is only a single phase which is in analogy to the thermal behavior of “ideal gas.” The dashed line is the isobar of P
=
P
c. The solid line is the isobar of P
<
P
c; one can see that there are two branches of black holes whose temperatures increase with the increase of horizon radius; one branch is in small radius region and the other one is in large radius region. Between the two branches, there is a region whose temperature decreases with the increase of horizon radius, which means the thermal expansion coefficient is negative; this region corresponds to an unstable phase. For appropriate value of temperature, two horizon radii are permitted; thus large/small black holes phase transition occurs. This kind of thermal behavior resembles the one of van der Waals gas/liquid system. The phase transition can also be seen in G
-
T plot, where “swallow tail” appears when P
<
P
c, just as the solid line shows us. The dashed line on G
-
T plot is the isobar of P
=
P
c, and the dotted line corresponding to P
>
P
c describes “ideal gas” phase.Figure 1
Isobaric plots forQ
YM
=
2, Q
e
=
1. The isobars on T
-
r
+ plot correspond one-to-one with the G
-
T plot. On both (a) and (b), from top to bottom, the pressures are, respectively, P
=
1.2
P
c (dotted), P
=
P
c (dashed), and P
=
0.8
P
c (solid).
(a)
(b)Now, let us consider phase transition of the black hole (14) for a general value of β. For general β, the pressure of the system is (36)
P
=
-
1
8
π
r
+
2
+
T
2
r
+
+
β
2
4
π
1
+
Q
e
2
β
2
r
+
4
-
β
2
4
π
+
Q
YM
2
8
π
r
+
4
. In this case, the critical equations (30) are too complicated to be solved analytically, we solve them numerically, the results are collected in Table 1. From the data in the table, one learns that, with the increase of β or Q
YM, T
c and P
c decrease while r
c and the ratio P
c
r
c
/
T
c increase. Figure 2 gives the isobaric plots near the critical point of P
c
=
0.0006631 and Figure 3 displays the effects of YM charge and BI parameter on the critical pressures. On T
-
r
+ plot of Figure 2, just as the case of β
→
∞, P
>
P
c isobar describes the “ideal gas” phase. The isobar of P
<
P
c shows that the “large black hole” region and “small black hole” region are thermodynamical stable while the medium region is thermodynamical unstable, which indicates the existence of phase transition. On G
-
T plot of Figure 2, the isobars show that when P
<
P
c phase transition characterized by the “swallow tail” occurs. On Figure 3(b) one can see that, although corresponding to different orders of magnitude of β, the dotted line (β
=
1) and the dashed line (β
=
10) are almost coincident, and the deviation of the solid line (β
=
0.1) from the dotted or dashed one is small; it seems that the influence of the values of β to critical pressure is not that large just as the authors found in [23]. This may be attributed to the parameter region we choose.Table 1
Critical values forQ
e
=
1.
β
Q
YM
r
c
T
c
P
c
P
c
r
c
/
T
c
10
3
7.7459654
0.0136979
0.0003316
0.1875150
10
2
5.4772185
0.0193717
0.0006631
0.1875024
10
1
3.4640315
0.0306296
0.0016579
0.1874989
1
2
5.4765161
0.0193728
0.0006632
0.1874807
0.1
2
5.4113483
0.0194752
0.0006714
0.1865541Figure 2
Isobaric plots forQ
YM
=
2, Q
e
=
1, and β
=
10. The isobars on T
-
r
+ plot and G
-
T plot are in one-to-one correspondence with each other. On both (a) and (b), from top to bottom, the pressures are, respectively, P
=
1.2
P
c (dotted), P
=
P
c (dashed), and P
=
0.8
P
c (solid).
(a)
(b)Figure 3
Isobaric plots that reflect the effects of YM charge (a) and BI parameter (b). On (a), the parameters of all the isobars are taken asQ
e
=
1, β
=
10, and P
=
0.7
P
c; the difference is that the solid, dashed, and dotted lines correspond to Q
YM
=
1,2
,
3, respectively. On (b), the parameters of all the isobars are taken as Q
e
=
1, Q
YM
=
2, and P
=
0.7
P
c; the solid, dotted, and dashed lines correspond to β
=
0.1,1
,
10, respectively.
(a)
(b)


## 4. Conclusion

In this paper, we construct a black hole solution of Einstein gravity coupled to BI electromagnetic field andS
U
(
2
) YM field. When the parameter ν
=
0, one reproduces the Einstein-BI balck hole in (A)dS space. When the BI parameter β
→
∞, the solution degenerates to a new black hole with electrostatic charge and magnetic YM charge. When β
→
0, the solution reduces to a black hole with only YM charge. When the parameters ν
=
0 and β
→
∞, one obtains the (A)dS-RN black hole. We calculate the thermodynamical quantities such as mass, temperature, entropy, electric charge, and YM charge and check that the fist law of thermodynamics is satisfied.We study the phase transition behaviors of the black hole in extended phase space by identifying the cosmological constant as pressure of the system. The critical equations can be solved analytically when takingβ
→
∞ limit while they have to be solved numerically for general β. After plotting the isobars one sees that when above the critical pressure there is only an “ideal gas” phase, when below the critical pressure, there are two stable regions corresponding, respectively, to large/small black hole phases with an unstable medium region, which implies the existence of the large/small black hole phase transition; this thermal behavior of black hole is in analogy to the one of van der Waals gas. The effects of β and Q
YM on critical pressures are also considered and displayed in the isobaric plots.



---

*Source: 2038202-2017-01-31.xml*