#### Abstract

We consider three-dimensional BTZ black holes with three models of nonlinear electrodynamics as source. Calculating heat capacity, we study the stability and phase transitions of these black holes. We show that Maxwell, logarithmic, and exponential theories yield only type one phase transition which is related to the root(s) of heat capacity, whereas, for correction form of nonlinear electrodynamics, heat capacity contains two roots and one divergence point. Next, we use geometrical approach for studying classical thermodynamical behavior of the system. We show that Weinhold and Ruppeiner metrics fail to provide fruitful results and the consequences of the Quevedo approach are not completely matched to the heat capacity results. Then, we employ a new metric for solving this problem. We show that this approach is successful and all divergencies of its Ricci scalar and phase transition points coincide. We also show that there is no phase transition for uncharged BTZ black holes.

#### 1. Introduction

One of the interesting subjects for recent gravitational studies is the investigation of three-dimensional black holes [1–9]. Considering three-dimensional solutions helps us to find a profound insight in the black hole physics, quantum view of gravity, and also its relations to string theory [10–14]. Moreover, three-dimensional spacetimes play an essential role in improving our understanding of gravitational interaction in low dimensional manifolds [15]. Due to these facts, some of physicists have an interest in the -dimensional manifolds and their attractive properties [16–27].

The Maxwell theory is in agreement with experimental results, but it fails regarding some important issues such as self-energy of point-like charges which motivates one to regard nonlinear electrodynamics (NED). There are some evidences that motivate one to consider NED theories: solving the problem of point-like charge self-energy, understanding the nature of different complex systems, obtaining more information and insight regarding quantum gravity, compatible with AdS/CFT correspondence and string theory frames, description of pair creation for Hawking radiation, and the behavior of the compact astrophysical objects such as neutron stars and pulsars [28–30]. Therefore, many authors investigated the black hole solutions with nonlinear sources [31–49].

On the other hand, thermodynamical structure of the black holes has been of great interest. It is due to the fact that, according to AdS/CFT correspondence, black hole thermodynamics provides a machinery to map a solution in AdS spacetime to a conformal field on the boundary of this spacetime [13–15]. Also it was recently pointed out that considering cosmological constant as a thermodynamical variable leads to the behavior similar to the Van Der Waals liquid/gas system [50–57]. In addition, phase transition of the black holes plays an important role in exploring the critical behavior of the system near critical points. There are several approaches that one can employ to study the phase transition. One of these approaches is studying the behavior of the heat capacity. It is argued that roots and divergence points of the heat capacity are representing two types of phase transition [58–64]. In addition, studying the heat capacity and its behavior enables one to study the thermal stability of the black holes [65–70].

Another approach for studying phase transition of black holes is through thermodynamical geometry. The concept is to construct a spacetime by employing the thermodynamical properties of the system. Then, by studying the divergence points of thermodynamical Ricci scalar of the metric, one can investigate phase transition points. In other words, it is expected that divergencies of thermodynamical Ricci scalar (TRS) coincide with phase transition points of the black holes. Firstly, Weinhold introduced differential geometric concepts into ordinary thermodynamics [71, 72]. He considered a kind of metric defined as the second derivatives of internal energy with respect to entropy and other extensive quantities for a thermodynamical system. Later Ruppeiner [73, 74] introduced another metric and defined the minus second derivatives of entropy with respect to the internal energy and other extensive quantities. The Ruppeiner metric is conformal to the Weinhold metric with the inverse temperature as the conformal factor. It is notable that both metrics have been applied to study the thermodynamical geometry of ordinary systems [75–80]. In particular, it was found that the Ruppeiner geometry carries information of phase structure of thermodynamical system. For the systems with no statistical mechanical interactions (e.g., ideal gas), the scalar curvature is zero and the Ruppeiner metric is flat. Because of the success of their applications to ordinary thermodynamical systems, they have also been used to study black hole phase structures and lots of results have been obtained for different types of black holes [81–86].

It is notable that these two approaches fail in order to describe phase transition of several black holes [87]. In order to overcome this problem, Quevedo proposed new types of metrics for studying geometrical structure of the black hole thermodynamics [88, 89]. This method was employed to study the geometrical structure of the phase transition of the black holes [58–64, 90–94] and proved to be a strong machinery for describing phase transition of the black holes. But the results of this approach was not completely coincided with the results of classical thermodynamics arisen from the heat capacity [87]. In [87], a new metric was proposed in which the denominator of its Ricci scalar is only constructed of numerator and denominator of the heat capacity. Several phase transitions of the black holes have been studied in context of the HPEM (Hendi-Panahiyan-Eslam Panah-Momennia) metric [87].

In this paper, we study thermal stability and phase transition of the BTZ black holes in the presence of several NED models in context of heat capacity. Then, we employ Weinhold, Ruppeiner, and Quevedo methods for studying geometrothermodynamics of these black holes. We will see that Weinhold and Ruppeiner metrics fail to provide fruitful results and the consequences of the Quevedo approach are not completely matched to the heat capacity results. Then, we employ the HPEM metric and study the phase transition of these black holes in context of geometrothermodynamics. We end the paper with some closing remarks.

#### 2. Nonlinearly Charged BTZ Black Hole Solutions

The -dimensional action of Einstein gravity with NED field in the presence of cosmological constant is given by where is the Ricci scalar and the cosmological constant is in which is a scale factor. Also is the Lagrangian of NED, in which we consider three models. First model was proposed by Hendi (HNED) [95], second one is Soleng theory (SNED) [96], and third one is correction form of NED (CNED) [69]:where and are called the nonlinearity parameters, the Maxwell invariant in which is the electromagnetic field tensor, and is the gauge potential. We should note that for (HNED and SNED branches) and (CNED branch) the Maxwell Lagrangian can be recovered.

The nonlinearly charged static black hole solutions can be introduced with the following line element: where the metric function was obtained in [70, 95]where and and are integrations constant which are related to mass and the electric charge of the black hole, respectively. In addition, , , and the special function .

The entropy and the electric charge of the obtained NED black hole solutions can be calculated with the following forms [70, 95]: where is the event horizon of black hole. On the other hand, the quasilocal mass, which is related to geometrical mass, can be obtained as [70, 95]

Regarding (5) with (6) and obtaining by using , for these three cases of BTZ black holes, we can write where and . Having conserved and thermodynamic quantities at hand, it was shown that the first law of thermodynamics may be satisfied [95]. The main goal of this paper is investigating thermal stability and phase transition for these black holes.

##### 2.1. Heat Capacity

In order to investigate thermal stability and phase transition, one can usually adopt two different approaches to the matter at hand. In one method, the electric charge is considered as a fixed parameter and heat capacity of the black hole will be calculated. The positivity of the heat capacity is sufficient to ensure the local thermal stability of the solutions and its divergencies are corresponding to the phase transition points. This approach is known as canonical ensemble. Another approach for studying thermal stability of the black holes is grand canonical ensemble. In this approach, thermal stability is investigated by calculating the determinant of Hessian matrix of with respect to its extensive variables. The positivity of this determinant also represents the local stability of the solutions. Although these two approaches are different fundamentally, one expects the results to be the same for both ensembles; that is, ensemble independent. Here we use the first method for studying thermal stability. For this purpose, the system is considered to be in fixed charge and the heat capacity has the following form:

It is notable that when we study heat capacity for investigating the phase transition, we encounter with two different phenomena. In one, the changes in the signature of the heat capacity are representing a phase transition of the system. In other words, if the heat capacity is negative, then the system is in thermally unstable phase, whereas, for the case of the positive , the system is thermally stable. The roots of the heat capacity in this case are representing phase transition points which means one should solve the following equation: where hereafter we call this type of the phase transition as type one. It is a matter of calculation to show that

The other case of phase transition is the divergency of the heat capacity. In other words, the singular points of the heat capacity are representing places in which system goes under phase transition. This assumption leads to the fact that the roots of the denominator of the heat capacity are representing phase transitions. Therefore, we have the following relation for this type of phase transition: where we call this type of the phase transition as type two. Due to economical reasons, we did not write the explicit relations of for different BTZ solutions.

##### 2.2. Geometrical Thermodynamics

In order to have an effective geometrical approach for studying phase transition of a system, one can build a suitable thermodynamical metric and investigate its Ricci scalar. Thermodynamical metrics were introduced based on the Hessian matrix of the mass (internal energy) with respect to the extensive variables. Therefore, although the electric charge is a fixed parameter for calculating the heat capacity in canonical ensemble, it may be an extensive variable for constructing thermodynamical metrics. In this method we expect that TRS diverges in both types of the mentioned phase transition points. In other words, the denominator of TRS must be constructed in a way that contains roots of the denominator and numerator of the heat capacity. In what follows, we will study the denominator of TRS of the several geometrical approaches and follow the recently proposed thermodynamical metric whose denominator only contains numerator and denominator of the heat capacity and, therefore, divergencies of TRS coincide with roots and divergences of the heat capacity.

In order to find the roots and divergence points of heat capacity, we should solve its numerator and denominator, separately. Solving the mentioned equations with respect to entropy leads to where

It is evident from obtained equation for HNED and SNED that there is only one real positive entropy in which heat capacity vanishes. Interestingly, in case of CNED theory, we find two roots for heat capacity. It is evident that the roots are increasing functions of the electric charge in these theories. As for nonlinearity parameter, in case of the HNED and SNED theories, the root is an increasing function of . For CNED theory, the smaller root is an increasing function of , whereas the larger root is a decreasing function of it.

Now we are in a position to study the existence of the type two phase transition point which is related to divergency of the heat capacity. Considering HNED branch of (13), one finds that is not real for all values of and . Therefore, there is no physical divergence point for the heat capacity of HNED model. Next, for the case of* SNED model*, the same behavior is seen. In other words, since for , we cannot obtain real for all values of and . Next, we should investigate* CNED model*. Regarding (13), we find that there is a divergence point for heat capacity (we should note that in this paper we consider positive ). In other words, this theory of nonlinear electromagnetic field enjoys the phase transition of type two. The divergence point is an increasing function of the nonlinearity and electric charge parameters. It is worthwhile to mention that, for the case of vanishing nonlinearity parameter in this theory, goes to zero. In other words, there is no divergence point for heat capacity of Maxwell theory. It is also clear that, in case of chargeless BTZ black holes, there is no root and divergence point for heat capacity. In other words, in case of chargeless BTZ black holes, there is not any kind of phase transition.

###### 2.2.1. Weinhold and Ruppeiner Metrics

The Weinhold metric was given in [71, 72] where and also , where denotes other extensive variables of the system. In case of Weinhold approach, one is considering the mass of the system as potential, other parameters such as entropy and electric charge as extensive parameters and related quantities such as temperature, and electric potential as intensive parameters.

The Ruppeiner metric was defined as [73, 74] where and . In this case, the thermodynamical potential is entropy. It is worthwhile to mention that, according to the proposal of the Quevedo, these two approaches are related to each other by a Legendre transformation [97].

Taking into account thermodynamical metrics of Weinhold and Ruppeiner, one can obtain their Ricci scalars. Since we would like to investigate divergence points of TRS, , we focus on its denominator (). One finds where and .

###### 2.2.2. The Quevedo Metrics

The Quevedo metrics have two kinds with the following forms [88, 89]: where iswithTaking into account Quevedo metrics, one can find that their related (denomerator of) Ricci scalars can be written as

###### 2.2.3. HPEM Metric

In order to avoid any extra divergencies in TRS which may not coincide with phase transitions of types one and two and also ensure that all the divergencies of the TRS coincide with phase transition points of the both types, HPEM metric was introduced [87]:

In this case, we have considered the total mass as thermodynamical potential and entropy and electric charge as extensive parameters. Calculations show that denominator of TRS leads to

#### 3. The Results of Various Approaches

Here, we investigate phase transitions of black holes using geometrothermodynamics. For this purpose, we used thermodynamical metrics introduced in previous section for the black holes solutions obtained in Section 2.

For Weinhold metric, none of divergencies of the Ricci scalar coincide with roots of the heat capacity in every theories of NED models that we have considered in this paper. On the other hand, one of the divergencies of TRS and divergence point of the heat capacity in CNED theory coincide with each other. It is notable that, in cases of the HNED and SNED theories, there is one divergence point for TRS (up panels of Figure 1) whereas, for CNED, there are two divergencies (down panels of Figure 1).

In case of Ruppeiner metric, for* HNED model* and* SNED model*, there is a root for heat capacity in which Ricci scalar of the Ruppeiner metric has a divergency. But there is also another divergence point for Ricci scalar which does not coincide with any phase transition point (up panels of Figure 2). Therefore, there is an extra divergence point. In case of the* CNED model*, two roots and one divergence point are observed for heat capacity in which Ricci scalar has related divergencies. In addition to these divergence points, one extra divergence point is also observed which is not related to any phase transition point (down panels of Figure 2).

As for Quevedo metrics, for case I, similar behavior as Weinhold is observed for all three theories of NED (Figure 3). On the other hand, in case of the other metrics of Quevedo, two divergence points for Ricci scalar were observed for HNED and SNED theories. One of these divergence points coincides with root of the heat capacity for these two nonlinear theories whereas the other one does not (up panels of Figure 4). In case of CNED theory, all the divergence points of TRS coincide with phase transition points except one (down panels of Figure 4). In other words, Quevedo’s metric predicts an extra divergence point, corresponding to the equation , which is not predicted in classical black hole thermodynamics.

**(a)**

**(b)**

**(c)**

**(d)**

**(a)**

**(b)**

**(c)**

**(d)**

It is evident that, in case of HPEM, all types of phase transition points of heat capacity coincide with divergencies of TRS of HPEM method (Figure 5). In other words, independent of the nonlinear theory under consideration, the HPEM method provides a machinery in which no extra divergency for TRS is observed and divergence points of TRS and phase transition points coincide. Another interesting and important property of the HPEM method is the behavior of TRS near divergence point for different types of phase transition. As one can see, the signature and behavior of TRS near divergence point for phase transition type one and two are different. Therefore, independent of studying the heat capacity, one can distinguish these types of phase transition from one another only by studying the behavior of TRS.

**(a)**

**(b)**

**(c)**

**(d)**

#### 4. Closing Remarks

In this paper, we have considered BTZ black holes, in the presence of three models of NED. We studied stability and phase transitions related to the heat capacity of the mentioned black holes. Next, we employed geometrical approach to study the thermodynamical behavior of the system. In other words, we have studied phase transitions of the system through Weinhold, Ruppeiner, and Quevedo methods. Also, we used the recently proposed approach to study geometrical thermodynamics.

We found that the Weinhold and Ruppeiner metrics for studying these BTZ solutions fail to provide a suitable result. In addition, the divergence points of the Quevedo TRS were not completely matched with the phase transition points of the heat capacity results. In other words, in these approaches, the existence of extra divergencies was observed which were not related to any phase transition point in the classical thermodynamics. In some of these approaches, no divergency of TRS coincided with phase transition points. In order to obtain a consistent results with the classical thermodynamic consequences (the heat capacity), we employed a new thermodynamical metric. In this approach, all the divergencies of TRS coincided with phase transition points. In other words, roots and divergence points of the heat capacity of the BTZ black holes in the presence of each nonlinear models matched with divergencies of TRS of this metric.

Also we found that, in case of HNED and SNED theories, there is no divergency for heat capacity. It means that, like Maxwell theory, these two theories have no second type phase transition. These two nonlinear theories of electrodynamics preserved the characteristic behavior of the Maxwell theory in case of heat capacity. On the other hand, for* CNED model*, two roots and one divergence point were found for the heat capacity. In other words, due to contribution of the nonlinear electromagnetic field, heat capacity enjoys the existence of one more phase transition point of type one and a phase transition of type two. In essence, this theory is a generalization of the Maxwell theory. But this generalization added another property to heat capacity that was not observed for the Maxwell theory.

Finally, it is worthwhile to mention a comment related to Legendre invariancy. It was shown that [98] the Legendre invariance alone is not sufficient to guarantee a unique description of thermodynamical metrics in terms of their curvatures. In addition to Legendre invariancy, one needs to demand curvature invariancy under a change of representation. Therefore, it will be worthwhile to investigate both Legendre and curvature invariancies. In addition, it will be interesting to think about the fundamental relation between the following two issues: (I) agreement of thermodynamical curvature results with usual thermodynamical approaches (such as the heat capacity) and (II) curvature invariancy in addition to the Legendre invariancy. It is also worthwhile to probe the fundamentality of cases (I) and (II) to find that considering which one may lead to satisfy the other one. Although first issue has been investigated for special cases [98], the second one has remained unanalyzed yet. We may address them in an independent work in the future.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The authors would like to thank the anonymous referee for valuable suggestions. They also thank the Shiraz University Research Council. This work has been supported financially by the Research Institute for Astronomy and Astrophysics of Maragha, Iran.