Thermodynamic Geometry and Phase Transitions of Dyonic Charged AdS Black Holes

We investigate phase transitions and critical phenomena of four dimensional dyonic charged AdS black holes in the framework of thermodynamic geometry. In a mixed canonical grand canonical ensemble with a fixed electric charge and varying magnetic charge these black holes exhibit liquid gas like first order phase transition culminating in a second order critical point similar to the Van der Waals gas. We show that the thermodynamic scalar curvature R for these black holes follow our proposed geometrical characterization of the R-crossing Method for the first order liquid gas like phase transition and exhibits a divergence at the second order critical point. The pattern of R crossing and divergence exactly corresponds to those of a Van der Waals gas described by us in an earlier work.


Introduction
In the last two decades the investigation of the thermodynamics of black holes has evolved into one of the most crucial issues in the study of quantum theories of gravity [1][2][3][4][5][6][7][8] 1 . The thermodynamics of asymptotically Anti de Sitter (AdS) black holes has assumed critical importance in this context owing to the famed AdS-CFT duality. Although a complete understanding of the microscopic statistical basis of black hole thermodynamics has eluded a clear explanation the study of the thermodynamics of AdS black holes has revealed a rich variety of phase structures and critical phenomena [9][10][11][12]. Remarkably in [10][11][12][13] an interesting similarity between the phase diagrams of electrically charged Reisner Nordstrom (RN) AdS black holes and that of a Van der Waals fluid could be established.
In a completely different context a geometrical approach to thermodynamics and phase transitions has been developed through the work of Weinhold [14,15] and Ruppeiner [16]. In this framework a Riemannian geometry of classical thermodynamic fluctuations, with a Euclidean signature, may be defined for the equilibrium thermodynamic state space of any thermodynamic system. The probability distribution of such fluctuations in a Gaussian approximation could then be related to the positive definite invariant interval defined by this geometry. It was found that the thermodynamic scalar curvature arising from this equilibrium state space geometry encodes the microscopic interactions underlying the thermodynamic system. From standard scaling and hyperscaling arguments it was possible to show that the thermodynamic scalar curvature was proportional to the correlation volume of the system and this diverged at a critical point of second order phase transition. This direct connection between the thermodynamic scalar curvature and the microscopic correlation length make this geometrical framework extremely suitable for the description of black hole thermodynamics where a complete microscopic structure is till an elusive issue. In fact the application of this framework has provided interesting insights into the thermodynamics of both non extremal and extremal black holes [17][18][19][20][21][22][23][24][25][26][27][28].
Although the geometrical framework described characterized second order phase transitions through the divergence of the thermodynamic scalar curvature, a similar characterization based on the thermodynamic curvature, for first order phase transitions was lacking. One of the authors (GS) in the collaborations [29][30][31] showed that the thermodynamic scalar curvature exhibited a multiple valued branched structure similar to the free energy at a first order phase transition which could be used for their characterization. Application of this to the study of liquid gas like first order phase transitions for a large class of asymptotically AdS black holes provided a confirmation of this alternative characterization. Following this in [32] a complete alternative geometrical characterization for first order phase transition, phase coexistence and supercritical phenomena could be established in the framework of thermodynamic geometries. This characterization extended Widoms [33] microscopic approach based on the correlation length to describe first order phase transitions, to propose the equality of the correlation length in the coexisting phases. In the framework of thermodynamic geometries the above characterization translates to the equality of the thermodynamic scalar curvature R for the coexisting phases. This essentially amounted to the crossing of the branches of the multiple valued curvature R for the coexisting phases. This R-Crossing Formula [32] applied to conventional simple fluids as an alternative to the Maxwell equal area construction or the equality of the free energy, showed remarkable correspondence with experimental data . Additionally this framework [32] also provided a first theoretical method for the construction of the Widom line, the locus of the maxima of the correlation length in the supercritical regime. The Widom line serves as a line of dynamical crossover for fluid properties that seems to retain the memory of the distinct subcritical phases [33][34][35][36]. Thus our construction led to a complete unified geometrical framework for the characterization of subcritical, critical and supercritical phenomena based on the thermodynamic scalar curvature and has led to interesting further applications [37][38][39][40].
Following earlier studies of the thermodynamics of electrically charged AdS black holes [10,11,13] in [41] the authors have investigated the thermodynamics of dyonic charged AdS black holes [42] in four dimensions which involves both electric and magnetic charges. These dyonic charged black holes are solutions to the equations of motion of a Einstein-Maxwell theory with a negative cosmological constant. In [41] the cosmological constant Λ is considered as the thermodynamic pressure [43] and it was shown that in a mixed canonical-grand canonical ensemble with a fixed magnetic charge and a varying electric charge the black hole undergoes a liquid gas like first order phase transition culminating in a critical point that resembles the phase diagram of a Van der Waals fluid [44]. A similar result is also expected to follow for the fixed electric charge and varying magnetic charge scenario due to the symmetrical fashion in which the charges occur in the Smarr formula for the black hole mass. As these black holes exhibit first order phase transitions in mixed ensembles it is naturally interesting to study the phase structure using the thermodynamics scalar curvature R referred earlier in the framework of thermodynamic geometries. In this article we investigate the thermodynamic geometry of such dyonic charged AdS black holes in four dimensions using mixed canonical grand canonical ensembles and study their phase transition and critical phenomena. However in our analysis we adopt a more conventional thermodynamic approach without considering the cosmological constant as a thermodynamic pressure. A study of the thermodynamics of this simpler case reproduces the first order phase transition and critical phenomena in a mixed ensemble observed in [41] where the cosmological constant is considered as the thermodynamic pressure. In particular we implement our R-Crossing Formula for the thermodynamic scalar curvature R to characterize the phase coexistence in the first order liquid gas like phase transition and show that the curvature R diverges at the critical point in accordance with the conclusions in [29][30][31][32].
This article is organized as follows, in Section 2 we briefly review the four dimnsional dyonic charged AdS black holes as solutions to the equation of motion of a Einsteic Maxwell AdS action. In Section 3 we discuss the thermodynamics of these dyonic charged AdS black holes in a more conventional scenario without considering the cosmological constant as a thermodynamic pressure and reproduce the first order phase transition and critical phenomena described in [41]. In Section 4 we briefly review the essential elements of thermodynamic geometries and employ this geometrical framework to investigate the thermodynamics and phase transitions of the four dimensional dyonic charged AdS black holes and implement the R crossing Method to characterize the first order liquid gas like phase transition for these black holes. In Section 5 we present a summary of our results and discussions including future open issues.

Dyonic black holes
In this section we briefly review the dyonic charged black hole solutions to the equations of motion of the Einstein Maxwell theory in a four dimensional AdS space time as described in [41]. The action for this theory may be expressed as where R is the scalar curvature, F µν is the electromagnetic field strength tensor and l is the AdS length scale and the gravitational constant G has been set equal to one. The equations of motions may be written down as A static spherically symmetric solution to the equations of motion are where, A µ , q e , q m and M are identified as the electromagnetic four-potential, electric charge, magnetic charge and mass of the black hole respectively. The quantity r + is the outer horizon radius of the black hole given by the zeroes of the lapse function f (r) described above.
A dyonic black hole is associated with a dyonic charge that involves both an electric (q e ) and a magnetic (q m ) charge [13,41]. The presence of magnetic charge leads to a corresponding magnetic potential (φ m ) in addition to the electric potential (φ e ). The expression for electric potential (φ e ) may be given as, here, φ E may be defined as the asymptotic value of the electric potential measured at infinity.

Thermodynamics of Dyonic charged black holes
It is now well established that in a semi classical framework, black holes may be considered as thermodynamic systems characterized by a Hawking temperature T and an entropy S. For the case of the four dimensional dyonic charged AdS black hole these may be expressed as The first law of black hole thermodynamics in this case may be written down as where, T and S are the Hawking temperature and the entropy of the dyonic charged AdS black hole respectively. Using this the Smarr formula for the mass M of the black hole in terms of s, q e and q m using Eq.(5) may be obtained as in [41] Following [29] we have scaled out the AdS length scale l which do not figure in our subsequent analysis. To analyze the thermodynamics of the dyonic charged AdS black hole we consider a mixed ensemble that is grand canonical with respect to the varying electric charge φ e but canonical with respect to the magentic charge q m which is fixed. Since both the magnetic as well as the electric charges enter in the Smarr relation given by Eq.(11) in a symmetric fashion, hence the alternative mixed ensemble with a varying magnetic charge q m and fixed electric charge q e is expected to have an identical thermodynamic behavior and does not need to be considered seprately. Using the Smarr formula the Hawking temperature T may be expressed as The Gibbs free energy specific to the mixed ensemble involving a varying electric charge q e and a fixed magnetic charge q m is as follows where φ e is the fixed electric potential.
It is now possible to express the Gibbs free energy G given above in terms of the variables φ e , q m and S as Similarly the temperature of the dyonic charged AdS black hole may now be expressed in terms of the variables S, φ e and q m as Now using these equations (14) and (15) for the Gibb' s free energy and the Hawking temperature it is straightforward to obtain the variation of the Gibb's free energy G with respect to the temperature T as shown in Fig.(1). The plot for the free energy G shows a typical multiple valued branched swallowtail behavior characteristic of a first order phase transition at a certain temperature T . This behavior disappears as the value of the magnetic charge q m is increased demonstrating the culmination of the subcritical phase coexistence to a critical point of a second order phase transition. This may be observed from the Fig.(1) where the blue curve corresponds to the free energy near the critical point. The red curve in Fig.(1) corresponds to the value q m = 0.095 for the magnetic charge which is below the critical value q m = 0.107. The subcritical red curve possesses three distinct branches that form the typical swallowtail structure for the fixed value of the magnetic charge q m = 0.095 whereas above the critical value q m = 0.107 this structure disappears. It may be noted that the qualitative nature of the plots do not change with the variation of electric potential φ e .  The first order phase transition described here occurs between the two distinct phases typified by small and large black holes characterized by their (outer) horizon radius r + , as shown in the Fig.(2), where the temperature T is plotted against the entropy S. Branch I or the 'small black hole' branch here corresponds to the 'liquid-like' phase and branch III or the 'large black hole' branch corresponds to the 'gas-like' phase in relation to a the usual liquid-gas phase transition. The three branches of the red curve in the T − S plot labeled as I, II and III correspond to the three respective branches of the red curve in the G − T plot displayed in Fig.(1). As one increases the temperature T the system moves along the branch-I, reaches the point of crossing of the branch-I and branch-III, jumps to the branch-III and continues along that branch. Branch-II is least favored as it corresponds to a higher free energy G than the other two branches at any given temperature T . As is well known this branch corresponds to a meta-stable state where the specific heat becomes negative. As one proceeds beyond the critical point, branch-II disappears and this corresponds to the usual uniform supercritical phase. Notice that in the limit q m → 0, the resulting black hole is essentially a conventional electrically charged RN-AdS black hole in grand-canonical (fixed φ e ensemble). In this limit the only possible phase transition is the usual Hawking-Page transition [5,9] between a thermal AdS space time and a black hole solution. In this limit the small black hole branch coincides with the T -axis in G − T plot.

Thermodynamic Geometry of a Dyonic Charged Black Hole
In this section we briefly collect the relevant results of thermodynamic geometry and then procced to obtain the thermodynamic geometry and the scalar curvature of the dyonic charged black hole in the mixed canonical-grand canonical ensemble elucidated earlier. As mentioned in the introduction it is possible to associate a curved Riemannian geometry with an Euclidean signature based on the fluctuations of the thermodynamic variables in the equlibrium state space of the system. The positive definite invariant line element for this geometry is then related to the probability distribution of the fluctuations connecting such equlibrium thermodynamic states in a Gaussian approximation. The Riemannian metric describing this geometry may be obtained from the thermodynamic potential appropriate to the representation being considered. In the entropy representation [32], the metric is given by the Hessian of the total entropy S = S(U, V, N ) with respect to the extensive variables with the volume V held fixed and U, N being the the internal energy and the number density respectively. The metric may be represented as where x µ are the extensive variables which serve as the coordinates in the state space subject to the condition dS ≥ 0 in accordance with the second law of thermodynamics and s is the combined entropy of the system and its surroundings. The volume V is considered to be infinite in the thermodynamic limit and in this case s may be considered to be just the entropy of the system only. Alternatively in the energy representation [14,15], the thermodynamic metric is obtained from the Hessian of the internal energy U with respect to the extensive variables and may be expressed as Using scaling and hyperscaling arguments the thermodynamic scalar curvature R obtained from this metric may be shown to be proportional to the correlation volume of the system as where ξ is the correlation length and d is the physical dimension of the system. Thus the thermodynamic scalar curvature R actually relates to the underlying microscopic interactions with R = 0 describing a non interacting microscopic basis. For our case of the four dimensional dyonic charged AdS black hole in the mixed canonical-grand canonical ensemble it is convenient to use the energy representation due to Weinhold [15] to obtain the thermodynamic metric. From the Smarr formula Eq.(11) the extensive variables which serve as the coordinates for the geometry are x µ = {s, q e } the electric charge and the entropy respectively. Here the magnetic charge q m is held fixed. The thermodynamic metric may then be expressed as It is to be noted that the internal energy of the black hole is determined to be its ADM mass M in the formulation for black hole thermodynamics which we are using. The corresponding metric components in this representation with the coordinates x µ = {s, q e } may be expressed as follows g ss = 3π 2 q 2 e + q 2 m − πs + 3s 2 2s (−π 2 (q 2 e + q 2 e ) + πs + 3s 2 ) , g sqe = g qes = − 2π 2 q e −π 2 (q 2 e + q 2 m ) + πs + 3s 2 , From the metric in the energy representation in Eq. (20) using the preceeding formula the thermodynamic scalar curvature R in terms of the variables s, q e and q m may be expressed from the standard Riemannian geometry formulation as, N = (π 2 q 2 m + 3s 2 )(3π 4 (q 2 e + q 2 m ) 2 + π 3 (−3q 2 e + q 2 m )s + 12π 2 (q 2 e + 3q 2 m )s 2 − 9πs 3 + 9s 4 ), D = s(π 2 (q 2 e + 3q 2 m ) − πs + 3s 2 ) 2 (−π 2 (q 2 e + q 2 m ) + πs + 3s 2 ). (21) For our purpose of analyzing the phase structure of the four dimensional dyonic charged black hole it would be convenient to express the thermodynamic scalar curvature R in terms of the alternative variables s, φ e and q m as A = (π 2 q 2 m + 3s 2 )(3π 4 q 4 m + 9s 4 + 3πs 3 (−3 + 4φ 2 e ) + π 3 q 2 m s(1 + 6φ 2 e ) + 3π 2 s 2 (12q 2 m − φ 2 e + φ 4 e )), B = s(3π 2 q 2 m + 3s 2 + πs(−1 + φ 2 e )) 2 (−π 2 q 2 m + 3s 2 + π(s − sφ 2 e )). (22) The phase structure of the dyonic charged black hole may now be investigated by studying the variation of the thermodynamic scalar curvature R with the temperature T using the equations (22) and (15) 2 . A comparison of the plots for the free energy of the dyonic charged AdS black hole in Fig.(1) and those for the thermodynamic scalar curvature R with the temperature T in Fig.(3) clearly shows the correspondence between the two. As mentioned earlier the thermodynamic scalar curvature R is a multivalued function of the thermodynamic parameters in the neighbourhood of a first order phase transition simmilar to the free energy G. The red curve in the Fig.(1) for the value of the magnetic charge q m = 0.095 corresponds to the subcritical regime accordingly has three distinct branches Moving in the direction of the increasing temperature T we first encounter the branch on the left that corresponds to the 'small black hole' phase (I).
As the temperature T is increased, the system moves towards the right in the graph to reach the first order phase transition temperature at T = 0.2137 where the 'small black hole' branch (I) and the 'large black hole' branch (III) intersect each other forming a 'swallowtail'. At this temperature, the black hole makes a transition from the 'small black hole' phase to the 'large black hole' phase. As the temperature is further increased the system continues along the branch corresponding to the 'large black hole' phase. Now from the red curve for the value of the magnetic charge q m = 0.095 shown in the Fig.(3a) for the R vs. T plot a behaviour simmilar to that of the red curve in free energy G vs the temperature T plot described above may be observed. In this plot also, in the Fig.(3a) moving in the direction of the increasing temperature T the branch on the left that corresponds to the 'small black hole' pahse (I) is first encountered. As the temperature T is further increased, the system moves towards the right to reach the first order phase transition temperature at T = 0.2174 where the 'small black hole' branch (I) and the 'large black hole' branch (III) intersect each other. For higher temperatures the system continues along the branch corresponding to the'large black hole' phase. In table (1) we summarize the values of the first order phase transition temperatures T f obtained from the free energy G vs. the temperature T plots in Fig.(1) and the thermodynamic scalar curvature R vs T plots in the Fig.(3). It may be noted from the preceeding table that the difference in the values of the temperature T f for the first order phase transition obtained from both the G − T plots and the R − T plots respectively matches upto the second decimal place. At this point it must be emphasized that the first order phase transition temperature obtained from the free energy G vs the temperature T plot essentially corresponds to the Maxwell construction signifying equal free energy in the coexisting phases. Whereas the corresponding temperature obtained from the thermodynamic scalar curvature R vs the temperature T plot through the crossing of the distinct branches is characterized by the R Crossing Method proposed by us as an alternative to the Maxwell construction outlined in [29,30,32]. This essentially corresponds to the equality of the correlation length ξ in the coexisting phases proposed as an extension of Widom's microscopic approach to phase transitions. It also must be emphasized here that unlike the subject of black hole thermodynamics where no experimental results are available for comparison our alternative R Crossing Method for characterization of first order phase transition exhibits a better fit with experimental data for simple fluids as described in [32].  Further observe now from the plots for the thermodynamic scalar curvature R vs temperature T that upon increasing the value of the magnetic charge q m , the 'small black hole' branch (I) and the 'large black hole' branch (III) recede away from each other as one approaches criticality. At the critical point the crossing of the two branches of the thermodynamic scalar curvature R completely disappear resulting in a divergence of R as predicted by Eq. (18). The blue curve in the Fig.(3) corresponds to the critical point of phase transition which clearly show that both the branches of the thermodynamic scalar curvature extend to the infinity. It may be observed from the Table. (1) that the thermodynamic scalar curvature diverges at a critical temperature which matches well with the critical temperature obtained from the free energy G vs temperature T plot (1). In other words the critical value of the temperature obtained from the thermodynamic scalar curvature R vs the temperature T plot in Fig.(3c) almost coincides with the value of the critical temperature obtained from the free energy G vs the temperature T plot ( the blue curve in Fig.(1)). Beyond the critical point the distinct subcritical phases of the small and large black holes disappear resulting in a uniform supercritical phase exactly similar to the case of a liquid gas phase transition. Notice however that for the case of fluids as described in [32], the fluid retains the memory of the subcritical phases across a locus of the maxima of the correlation length which is termed the Widom line in the literature, which serves as a crossover line for dynamical fluid properties and divides the supercritical regime into a liquid like and gas like region. The significance of such a Widom line in black hole thermodynamics however is not quite clear although just as in [32] it is possible to construct the Widom line through the locus of the maxima of the thermodynamic scalar curvature R in the supercritical regime.

Summary and Discussions
In summary we have investigated the thermodynamics, phase transition and critical phenomena of four dimensional dyonic charged AdS black holes using the framework of thermodynamic geometries. Unlike [41] where the cosmological constant is considered as the thermodynamic pressure in our study we have adopted a more conventional approach where the cosmological constant is a fixed parameter for the black hole solution. In this context we have reproduced the first order liquid gas like phase transition with coexisting phases in this more simple scenario using a mixed canonical grand canonical ensemble where the magnetic charge is held fixed and the electric charge is allowed to vary. Due to the symmetrical appearence of the electric and the magnetic charge in the Smarr formula, in the alternative setting of a fixed electric charge and varying magnetic charge the thermodynamic behaviour is expected to be identical and does not merit a separate analysis. Our analysis in this simpler scenario for the Gibb' s free energy as a function of the temperature clearly reproduces the typical swallowtail feature characterizing the first order liquid gas like phase transition and coexistence between a small and a large black hole phase in the subcritical regime.
The thermodynamic metric for the dyonic charged black hole is then obtained in the more convenient Weinhold picture as a Hessian of the internal energy with respect to the entropy S and the electric charge φ e . The thermodynamic scalar curvature R for this black hole is then obtained using standard Riemannian geometric techniques for the Ricci scalar. Analysis of the variation of the thermodynamic scalar curvature R as a function of the temperature T clearly demonstrates the expected branching behaviour signifying a first order phase transition. The plot of R vs the temperature T exhibits the crossing of the branches at the first order phase transition temperature and the fixed value of the magnetic charge φ m in accordance with our R Crossing Method that serves as an alternative to the conventional Maxwell construction. The value of the magnetic charge for the first order phase transition obtained through the R Crossing Method compares favorably and is almost identical to that obtained from the usual Maxwell construction and the equality of the Gibbs free energy G for the coexisting phases. For a further increase of the magnetic charge away from the subcritical regime for the first order phase transition values finally leads to a critical point signifying a second order phase transition where the thermodynamic scalar curvature R exhibits an expected divergence. The branching and subsequent crossing behaviour leading to a divergence at the critical point of the thermodynamic scalar curvature R is completely simmilar to that for the van der Waals gas described in [29,45].
Our analysis clearly establishes and further corroborates our approach for studying the phase structure of AdS black holes in the framework of thermodynamic geometry and proves our proposal for the R Crossing Method as an alternative to the conventional Maxwell construction and equality of the Gibbs free energy required for phase coexistence. An important open problem for future investigation in this scenario is to implement the framework of thermodynamic geometry to describe extended black hole thermodynamics which includes the cosmological constant as a thermodynamic pressure. In addition it would be interesting to construct the Widom line in the supercritical regime through the locus of the maxima of the thermodynamic scalar curvature R max as suggested by us and investigate the significance of this Widom line in the context of black hole phase transitions.