Thermodynamics, phase transitions and Ruppeiner geometry for Einstein–dilaton–Lifshitz black holes in the presence of Maxwell and Born–Infeld electrodynamics

In this paper, we first obtain the higher-dimen-sional dilaton–Lifshitz black hole solutions in the presence of Born–Infeld (BI) electrodynamics. We find that there are two different solutions for the cases of z=n+1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$z=n+1$$\end{document} and z≠n+1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$z\ne n+1$$\end{document} where z is the dynamical critical exponent and n is the number of spatial dimensions. Calculating the conserved and thermodynamical quantities, we show that the first law of thermodynamics is satisfied for both cases. Then we turn to the study of different phase transitions for our Lifshitz black holes. We start with the Hawking–Page phase transition and explore the effects of different parameters of our model on it for both linearly and BI charged cases. After that, we discuss the phase transitions inside the black holes. We present the improved Davies quantities and prove that the phase transition points shown by them are coincident with the Ruppeiner ones. We show that the zero temperature phase transitions are transitions in the radiance properties of black holes by using the Landau–Lifshitz theory of thermodynamic fluctuations. Next, we turn to the study of the Ruppeiner geometry (thermodynamic geometry) for our solutions. We investigate thermal stability, interaction type of possible black hole molecules and phase transitions of our solutions for linearly and BI charged cases separately. For the linearly charged case, we show that there are no phase transitions at finite temperature for the case z≥2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ z\ge 2$$\end{document}. For z<2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$z<2$$\end{document}, it is found that the number of finite temperature phase transition points depends on the value of the black hole charge and there are not more than two. When we have two finite temperature phase transition points, there is no thermally stable black hole between these two points and we have discontinuous small/large black hole phase transitions. As expected, for small black holes, we observe finite magnitude for the Ruppeiner invariant, which shows the finite correlation between possible black hole molecules, while for large black holes, the correlation is very small. Finally, we study the Ruppeiner geometry and thermal stability of BI charged Lifshtiz black holes for different values of z. We observe that small black holes are thermally unstable in some situations. Also, the behavior of the correlation between possible black hole molecules for large black holes is the same as for the linearly charged case. In both the linearly and the BI charged cases, for some choices of the parameters, the black hole system behaves like a Van der Waals gas near the transition point.


Abstract
In this paper, we first obtain the higher-dimensional dilaton-Lifshitz black hole solutions in the presence of Born-Infeld (BI) electrodynamics. We find that there are two different solutions for the cases of z = n+1 and z = n+1 where z is the dynamical critical exponent and n is the number of spatial dimensions. Calculating the conserved and thermodynamical quantities, we show that the first law of thermodynamics is satisfied for both cases. Then we turn to the study of different phase transitions for our Lifshitz black holes. We start with the Hawking-Page phase transition and explore the effects of different parameters of our model on it for both linearly and BI charged cases. After that, we discuss the phase transitions inside the black holes. We present the improved Davies quantities and prove that the phase transition points shown by them are coincident with the Ruppeiner ones. We show that the zero temperature phase transitions are transitions in the radiance properties of black holes by using the Landau-Lifshitz theory of thermodynamic fluctuations. Next, we turn to the study of the Ruppeiner geometry (thermodynamic geometry) for our solutions. We investigate thermal stability, interaction type of possible black hole molecules and phase transitions of our solutions for linearly and BI charged cases separately. For the linearly charged case, we show that there are no phase transitions at finite temperature for the case z ≥ 2. For z < 2, it is found that the number of finite temperature phase transition points depends on the value of the black hole charge and there are not more than two. When we have two finite temperature phase trana e-mail: mkzangeneh@shirazu.ac.ir b e-mail: mehdizadeh.mr@uk.ac.ir c e-mail: wang_b@sjtu.edu.cn d e-mail: asheykhi@shirazu.ac.ir sition points, there is no thermally stable black hole between these two points and we have discontinuous small/large black hole phase transitions. As expected, for small black holes, we observe finite magnitude for the Ruppeiner invariant, which shows the finite correlation between possible black hole molecules, while for large black holes, the correlation is very small. Finally, we study the Ruppeiner geometry and thermal stability of BI charged Lifshtiz black holes for different values of z. We observe that small black holes are thermally unstable in some situations. Also, the behavior of the correlation between possible black hole molecules for large black holes is the same as for the linearly charged case. In both the linearly and the BI charged cases, for some choices of the parameters, the black hole system behaves like a Van der Waals gas near the transition point.

Introduction
It has been over 40 years since Bekenstein and Hawking first disclosed that black hole can be considered as a thermodynamic system, with characteristic temperature and entropy [1][2][3][4]. Taking into account the fact that black holes have no hair, there are no classical degrees of freedom to account for such thermodynamic properties. It is a general belief that thermodynamic properties of a system may reflect the statistical mechanics of underlying relevant microscopic degrees of freedom. But the detailed nature of these microscopic gravitational states has remained a mystery. The Bekenstein-Hawking entropy, S = A/(4hG), depends on both Planck's constant and Newtonian gravitational constant, implying that the thermodynamics of black holes may relate quantum mechanics and gravity. Recently, there has been some progress in understanding the microscopic degrees of free-dom of the black hole entropy, for example in string theory [5][6][7] as well as loop quantum gravity [8][9][10]. But the accounts of the black hole entropy are not complete and they only work within some particular models and some special domains where string theory and loop quantum gravity can apply. Besides, despite counting very different states, many inequivalent approaches to quantum gravity obtain identical results and it is not clear why any counting of microstates should reproduce the same Bekenstein-Hawking entropy [11]. The statistical mechanical description of the black hole entropy is still not elegant.
On the other side, black hole can be heated or cooled through absorption and evaporation processes. According to Boltzmann's insight, if a system can be heated, it must have microscopic structures. Recently, in [12], possible microscopic structures of a charged anti-de Sitter black hole have been studied and some kind of interactions between possible micromolecules have been investigated by an interesting physical tool, the Ruppeiner geometry. Derived from the thermodynamic fluctuation theory, the Ruppeiner geometry [13,14] is considered a powerful tool in exploring the possible interactions between black hole microscopic structures. The sign of the Ruppeiner invariant R (the Ricci scalar of the Ruppeiner geometry) was argued to be useful for identifying the physical systems similar to the Fermi (Bose) ideal gas when R > 0 (R < 0) or the classical ideal gas when R = 0 [15]. Besides, the sign of the Ruppeiner invariant R can further be used to interpret the type of dominant interaction between molecules of a thermodynamic system. When R > 0, there is a repulsive interaction between molecules, when R < 0 the interaction is attractive, and for R = 0 there is no interaction in the microstates [16][17][18]. Moreover, the magnitude of the Ruppeiner invariant |R| measures the average number of correlated Planck areas on the event horizon for a black hole system [19]. For a review of the description of the Ruppeiner geometry in black hole systems, we refer to [20,21] and the references therein. Further studies of molecular interactions of black holes, based on the Ruppeiner geometry, have been carried out in [12,22,23].
The phase transition is another interesting topic in black holes thermodynamics. Davies discussed the thermodynamic phase transition of the black holes by looking at the behavior of the heat capacity [24][25][26]. He claimed that the discontinuity of the heat capacity marks a second order phase transition in black holes. However, it was argued that physical properties do not show any particularity at this discontinuity point if compared with other heat capacity values; for example the regularity of the event horizon is not lost and the black hole internal state remains uninfluenced [27]. Thus, it is hard to accept the discontinuity point of the heat capacity as a true physical point of the phase transition. Employing the Landau-Lifshitz theory of thermodynamic fluctuations [28,29], Pavon and Rubi gave a deep under-standing of the black hole phase transition [30,31]. They found that some second moments in the fluctuation of relevant thermodynamic quantities diverge when the black hole becomes extreme. This divergence shows that the thermodynamic fluctuation is tremendous and the rigorous meaning of the thermodynamical quantities is broken down. This is exactly the characteristic of the thermodynamic phase transition point. At this phase transition point, the Hawking temperature is zero which indicates that for the extreme black hole there is only super-radiation but no Hawking radiation, which is in sharp difference from that of the non-extreme black holes. Black holes phase transition in the context of Landau-Lifshitz theory have been investigated in [32,33]. Recently, further differences in dynamical properties before and after the black hole thermodynamical phase transition has been disclosed in [34][35][36]. Also, in [37][38][39][40][41][42], black hole phase transitions have been studied from holographic point of view. A question now arises: how we can further understand this macroscopic thermodynamic phase transition in black hole physics? for example whether there is a microscopic explanation of this thermodynamic phase transition. The Ruppeiner geometry is a possible tool we can use to investigate the thermodynamic phase transitions from microscopic point of view. This method is safer to determine true phase transitions than other methods since, regardless of the microscopic model, R has a unique status in identifying microscopic order (which is at the foundation of phase transitions at microscopic level) from thermodynamics [20,21]. Some attempts in this direction have been reported in [43][44][45][46][47][48][49][50][51][52][53][54]. In a recent work [46], it was found that the divergence of the Ruppeiner invariant coincides with the critical point in the phase transition in a holographic superconductor model. It is interesting to investigate whether the Ruppeiner geometry [20,21] can present us further reason to determine which of the thermodynamical discussions mentioned above is valid for describing the thermodynamical phase transition. In particular, we would like to explore whether the Davies phase transition conjecture can reflect some special properties in microstructures and be in consistence with the Ruppeiner geometry description. If the Davies conjecture does not have the microscopic explanation, we will further think about how to improve the Davies conjecture to describe the black hole phase transition.
We will employ the black hole in Lifshitz spacetime as a configuration to study our physical problems mentioned above. This spacetime was first introduced in [55], which respects the anisotropic conformal transformation t → λ z t, x → λ x, where z is dynamical critical exponent. For the Lifshitz spacetime, it is necessary to include some matter sources such as massive gauge fields [56][57][58][59][60] or higher-curvature corrections [61] to guarantee the asymptotic behavior of the Lifshitz black hole. It is difficult to find an analytic Lifshitz black hole solution for arbitrary z, although some attempts have been performed [62]. This makes the discussion of ther-modynamics for such a black hole difficult. Fortunately, in Einstein-dilaton gravity with a massless gauge field, it is possible to find an exact Lifshitz black solution for arbitrary z [63,64]. This model is suggested in the low energy limit of string theory [65]. While thermodynamical behaviors of uncharged and charged Einstein-dilaton-Lifshitz black holes have been revealed in [63,66] and [64], respectively, thermodynamics of uncharged Gauss-Bonnet-dilaton-Lifshitz solution has been studied in [67]. It is also interesting to study Lifshitz black hole solutions in the presence of other gauge fields such as the power-law Maxwell field [68], the logarithmic [69] and exponential [70] nonlinear electrodynamics. For example, thermodynamics and thermal and dynamical stabilities of Einstein-dilaton-Lifshitz solutions in the presence of power-law Maxwell field have been studied in [71]. In the context of AdS/CFT [72-74] application, the electrical conductivity were explored for exponentially [75] and logarithmic [76] charged Lifshitz solutions. In the present work, we shall consider the Born-Infeld (BI) nonlinear electrodynamics in the context of Einstein-dilaton-Lifshitz black holes. The motivation for considering the BI-dilaton action comes from the fact that the dynamics of D-branes and some soliton solutions of supergravity are governed by the BI action [77][78][79][80][81][82]. Besides, the low energy limit of open superstring theory suggest the BI electrodynamic action to be coupled to a dilaton field [77][78][79]. It is surprising that many years before the appearance of the BI action in superstring theory, in the 1930s, this nonlinear electrodynamics was introduced for the first time, with the aim of solving the infinite selfenergy problem of a point-like charged particle by imposing a maximum strength for the electromagnetic field [83].
In this paper, we will first look for a general (n + 1)dimensional Lifshitz black hole solution in the context of Einstein-dilaton gravity in the presence of BI electrodynamics. We will show that the general metric function has different solutions for z = n + 1 and z = n + 1 cases. It is important to note that the difference in the metric function has not been observed in the previous studies on Lifshitzdilaton black holes [71,75,76]. Based on this general solution, we will study the thermodynamics of Lifshitz-dilaton black holes coupled to a linear Maxwell field and BI nonlinear electrodynamics. We will show that the Hawking-Page phase transition [84] exists both in the presence of linear and nonlinear electrodynamics. There are some attempts to study phase transitions of uncharged Lifshitz solutions for fixed z [85] or in three [86] and four [87] dimensions. The Hawking-Page phase transition revealed in this paper is interesting, since it depends on different values of z in different spacetime dimensions in the presence of linear Maxwell and nonlinear BI electrodynamic fields. We will further concentrate our attention on understanding the thermodynamic phase transition from microstructures. We shall examine the relation between the Ruppeiner geometry and thermodynamical descriptions of the phase transition such as the Davies conjecture and the Landau-Lifshitz method. We try to give more microscopic understanding of the thermodynamical phase transitions in the black hole system. We explore the thermodynamic geometry (Ruppeiner geometry) for linearly and nonlinearly charged Lifshitz solutions separately and show the properties of interactions between possible black hole molecules. To the best of our knowledge, there is no study of thermodynamic geometry on Lifshitz solutions in the literature. Interestingly enough, by studying Ruppeiner geometry, we have found that our solutions show the Van der Waals like behavior near critical point in some cases.
The layout of the paper is as follows. In the next section, we give the basic field equations and obtain the BI charged Lifshitz-dilaton black hole solutions. In Sect. 3, we first explore the satisfaction of the first law of thermodynamics for Lifshitz-dilaton black holes in the presence of BI electrodynamics. Then we study different phase transitions including the Hawking-Page phase transition and phase transition at zero temperature for linearly and BI charged cases. In Sect. 4, we investigate thermodynamic geometry of the obtained solutions for linearly and nonlinearly BI charged cases by adopting the Ruppeiner approach. We finish with a summary and closing remarks in Sect. 5.

Action and asymptotic Lifshitz solutions
In this section, we intend to obtain exact (n +1)-dimensional dilaton-Lifshitz black holes in the presence of BI nonlinear electrodynamics. Our ansatz for the line elements of the spacetime is [64,88] where z(≥ 1) is the dynamical critical exponent and is an (n − 1)-dimensional hypersurface with constant curvature (n − 1)(n − 2) and volume ω n−1 . As r → ∞, the line elements ( 1) reduce asymptotically to the Lifshitz spacetime, On the other side, as pointed out above, we would like to consider BI nonlinear electrodynamics. In the absence of the dilaton field, the BI Lagrangian density is written as [83] where β is the Born-Infeld parameter related to the Regge slope where A μ is electromagnetic potential. One of the effects of the presence of a dilaton field is its coupling with the electromagnetic field. Thus, in the presence of the dilaton field we deal with a modified form for the BI Lagrangian density including its coupling with the dilaton scalar field [89,90], where λ is a constant. The Lagrangian density of the stringgenerated Einstein-dilaton model [65] with two Maxwell gauge fields [64] in the presence of BI electrodynamics can be written in an Einstein frame as where R is the Ricci scalar and and the λ i are some constants. In Lagrangian (5), and (B i ) μ is gauge potential. In the large β limit, L recovers the Einstein-dilaton-Maxwell Lagrangian in its leading order [64,71] Varying the action S = M d n+1 x √ −gL with respect to the metric g μν , the dilaton field and electromagnetic potentials A μ and (B i ) μ leads to the following field equations: where we use the convention X Y = ∂ X/∂Y . Using the metric ansatz (1), electromagnetic field equations (9) and (10) can be solved immediately as where ϒ = 1 + q 2 l 2z−2 /(β 2 r 2n−2 ), and (r ) can be obtained by subtracting (tt) and (rr) components of Eq. (7 ) and solving the resulting equation. We find Substituting Eqs. (11), (12) and (13) in the field equations (7) and (8), one can solve the equations for f (r ) to obtain where we should set so that the field equations are fully satisfied. In the solution (14), m is a constant which is related to the total mass of black brane as we will see in next section. The integral of the last term of f (r ) for z = n + 1 can be done in terms of hypergeometric function. Thus, f (r ) can be written as Note that solution (16) obviously satisfies the fact that which reproduces the result of [71] for every z in linear Maxwell case. The behaviors of the metric function for z = n + 1 and z = n + 1 have been depicted in Fig. 1a and b, respectively. It is notable to mention that in the case of z = n + 1, there is no Schwartzschild-like black hole since in this case f (r ) goes to positive infinity as r goes to zero. However, for z = n + 1, we may have Schwartzschild-like black hole (dash-dotted line in Fig. 1a) in addition to nonextreme (solid line) and extreme (dotted line) black holes and naked singularity (dashed line). For non-extreme case, there are two inner (Cauchy) and outer (event) horizons. In both Fig. 1a and b, we see that the larger the nonlinearity parameter β is, the smaller the distance between two inner and outer horizons is so that, for large enough β, we have just one horizon (extreme case) or naked singularities. The Schwartzschild-like case occurs for lower β's in the case of z = n + 1 as Fig. 1a shows.
As one can see in (17), the fourth term in expansions for both z = n + 1 and z = n + 1 cases reproduce the charge term of [71] in linear Maxwell case as one expects. The temperature of the black hole horizon can be obtained via where χ = ∂ t is the Killing vector and r + is the radius of event horizon. Using (18), one can calculate the Hawking temperature as where prime denotes the derivative with respect to r and ϒ + = ϒ (r = r + ). For the temperature one has the same formula as (19) for the two cases z = n + 1 and z = n + 1. One can check that, for large β, (19) reduces to the temperature of Einstein-Maxwell-dilaton-Lifshitz black holes [71], namely The entropy of the black holes can be calculated by using the area law of the entropy [2,91,92] which is applied to almost all kinds of black holes in Einstein gravity including dilaton black holes [93][94][95][96]. Therefore, the entropy of the black brane per unit volume ω n−1 becomes Having Eqs. (11), (13) and (15) at hand, we can find electromagnetic gauge potential A t = F rt dr in terms of hypergeometric function as The large β behavior of gauge potential is in agreement with [71] In next section, we will study thermodynamics of dilaton Lifshitz black holes in the presence of BI electrodynamics by seeking for satisfaction of thermodynamics first law through calculation of conserved and thermodynamic quantities. We also show that our Lifshitz solutions can exhibit the Hawking-Page phase transition. Then we discuss the inside phase transitions of our Lifshitz black holes.

First law of thermodynamics
This subsection is devoted to a study of the first law of thermodynamics for Lifshitz-dilaton black hole solutions in the presence of BI nonlinear electrodynamics. As the first step, we calculate the fundamental quantity for thermodynamics discussions namely mass. For this purpose, we apply the modified subtraction method of Brown and York (BY) [97][98][99]. In order to use this method, the metric should be written in the form For our case, it is clear that R = r and thus The metric of background is chosen to be the Lifshitz metric (24) i.e.
The quasilocal conserved mass can be obtained through where σ is the determinant of the boundary B metric, K 0 ab is the background extrinsic curvature, n a is the timelike unit normal vector to the boundary B and ξ b is a timelike Killing vector field on the boundary surface. Applying the above modified BY formalism, the mass of the space-time per unit volume ω n−1 can be calculated as where the mass parameter m can be obtained from the fact that f (r + ) = 0 as Now, we turn to a calculation of the electric charge of the solution. Using the Gauss law, we can calculate the electric charge via where are, respectively, the unit spacelike and timelike normals to the hypersurface of radius r . Using (30), the charge per unit volume ω n−1 can be computed as The electrostatic potential difference (U ) between the horizon and infinity is defined as Using Eqs. (22) and (32), one can obtain the electric potential which is the same for the cases z = n + 1 and z = n + 1. In order to investigate the first law of black hole thermodynamics, we should obtain the Smarr-type formula for mass (28). With Eqs. (29), (31) and (21) at hand, the mass can be written as a function of the extensive thermodynamic quantities S and Q in the form of where = 1 + π 2 Q 2 /(β 2 S 2 ). Calculations show that intensive quantities coincide with those computed by Eqs. (19) and (33). Therefore, the thermodynamics quantities satisfy the first law of thermodynamics, for the two solutions z = n + 1 and z = n + 1.
In the next part of this section, we will discuss the Hawking-Page and inside black hole phase transitions for our Lifshitz solutions.

Hawking-Page phase transition
As it is clear from Fig. 1, there are some parameter choices for which we have extreme black holes and therefore zero tem-perature. In addition, as one can see from Fig. 2, there are some other choices of the parameters that show a non-zero positive minimum for temperature T min . The influences of different parameters on T min can be seen from Fig. 2. When we increase the dimension n, T min increases too, while it decreases with increasing z. Comparing Fig. 2a and b, one finds that the effect of nonlinearity implies increasing in T min . The behaviors illustrated in Fig. 2 show a Hawking-Page phase transition for the obtained solutions. Let us have a closer look on Fig. 2. In the first part of the T -S curves where we have small black holes (note that S = r n−1 + /4), ∂ T /∂ S < 0, which implies negative heat capacity and therefore small black holes are thermally unstable. But in the large black hole part of the curves we have a positive heat capacity and therefore large black holes are thermally stable. In addition to small and large black holes, we have a thermal Lifshitz or radiation solution too. Since small black holes are thermally unstable, the system has a choice between large black hole and thermal Lifshitz solutions and chooses one of them according to the Gibbs free energy. The Gibbs free energy, can be obtained by using (19), (21) Fig. 4 The behavior of G versus T for nonlinearly charged case with l = 1, b = 1 and n = 3 and the critical exponent z makes T HP lower. Also, the lower the nonlinearity parameter β is, the lower the Hawking-Page temperature T HP is. Note that a lower β makes the electrodynamics more affected by nonlinearity.

Phase transitions inside the black hole
There are at least three well-known ways to discuss the phase transitions inside the black hole. Two of these ways are based on the macroscopic point of view and one of them is based on the microscopic viewpoint. The two macroscopic ways are Davies [24] and Landau-Lifshitz [28,29] methods, which discuss, respectively, the behavior of heat capacities and thermodynamic fluctuations. Thermodynamic geometry or Ruppeiner geometry [16,20,21] is the microscopic way; it discusses the phase transitions in addition to type and strength of interactions. In what follows, we discuss the relation between the phase transitions predicted by Ruppeiner geometry and the Davies method. Next, we will turn to the Landau-Lifshitz theory of thermodynamic fluctuations.

Ruppeiner and Davies phase transitions
In order to discuss thermodynamic geometry, one should study the divergences, sign and magnitude of Ricci scalar corresponding to Ruppeiner metric (usually called Ruppeiner invariant) to determine phase transitions and strength and type of dominated interaction between possible black hole molecules [16,20,21]. To do that, we define the Ruppeiner metric in (M, Q) space where the entropy S is the thermodynamic potential, The above metric can also be rewritten in the Weinhold form, The Ruppeiner invariant corresponding to (39) can be expressed in a general form as It is remarkable that at the point where M SS vanishes or equivalently heat capacity at constant charge C Q diverges, we have a thermally unstable system due to negativity of H M S,Q if M S Q = 0 (which occurs in many black hole systems). Thus, the heat capacity at constant charge C Q cannot be a suitable thermodynamic quantity to show the phase transition of such systems when we have two changing thermodynamic parameters, for instance S and Q. There is some work in the literature (for instance [102]) in which the correctness of the Ruppeiner method for recognizing the phase transitions has been judged by comparing the Ruppeiner and C Q transition points. This procedure of course seems to be incorrect according to what we pointed out above. Also, as we discussed in the introduction, the divergences of R are safer if we determine phase transitions. On the other hand, in [43,48], the authors have suggested some suitable thermodynamic quantities to show the phase transitions pre-dicted by the Ruppeiner invariant. These quantities are the specific heat at constant electrical potential, C U , the analog of the volume expansion coefficient, α, and the analog of the isothermal compressibility coefficient κ T defined as As one can see in the Appendix, these thermodynamic quantities have the forms It is obvious that these quantities show the same phase transitions as the Ruppeiner geometry because all of them diverge at roots of H M S,Q and C U vanishes at zero temperature where R diverges. To show the coincidence of the Ruppeiner phase transitions and C U divergences, some proofs have been presented in [50,54]. The above quantities can be considered as improved Davies quantities [24] which show the phase transitions to coincide with the Ruppeiner ones. In the next part, we study the Landau-Lifshitz theory of thermodynamic fluctuations to explore the possible signature of black hole phase transitions and the properties of black hole radiance.

Landau-Lifshitz theory (non-extreme/extreme phase transition)
Here, we seek any possible effect of a transition on black hole radiance by using the Landau-Lifshitz theory of thermodynamic fluctuations [28,29]. We focus on the (3 + 1)dimensional linearly charged case. The extension to higherdimensional or nonlinearly charged cases is trivial and gives no novel result. Based on Landau-Lifshitz theory [28,29], in a fluctuation-dissipative process, the fluxẊ i of a given thermodynamic quantity X i is given bẏ where a dot shows the temporal derivative and χ i and i j are, respectively, the thermodynamic force conjugate to the fluxẊ i and the phenomenological transport coefficients. In addition, the rate of entropy production is expressed bẏ where "+" ("−") holds for the entropy rate contributions which come from the non-concave (concave) parts of S. The second moments corresponding to the flux fluctuations are (we set k B = 1) where the mean value with respect to the steady state is denoted by the angular brackets and the fluctuations δẊ i are the spontaneous deviations from the value of steady state Ẋ i . To guarantee that the correlations are zero when two fluxes are independent, the Kronecker δ i j is put in Eq. (45). According to [71], the mass M, electric potential energy U and temperature T can be obtained for (3+1)-dimensional linearly charged case as We know that in extreme black hole case, the Hawking temperature on the event horizon vanishes and therefore in this case we have = 0. Using Eq. (46), we can obtain the entropy production rate aṡ where The mass loss rate is given by [103] dM dt The first term on the right side of Eq. (49) is the thermal mass loss corresponding to Hawking radiation, which is just the Stefan-Boltzmann law, with b = π 2 /15 (we seth = 1) as the radiation constant. The constant α depends on the number of species of massless particles and the quantity σ is the cross-section of geometrical optics. The second term on the right side of Eq. (49) is responsible for the loss of mass corresponding to charged particles. In fact, it is the U dQ term which arises in the first law of black hole mechanics.
With references to what was explained and computed above, one can calculate the second moments or correlation functions of the thermodynamical quantities, It is clear that the second moments δṠδṠ , δṪ δṪ and δṠδṪ diverge for the extreme black hole case where vanishes (see Eq. 47). It means that there is a phase transition in this case. This phase transition is between extreme and non-extreme black holes for which we have a sudden change in emission properties. In the non-extreme case, the black hole can give off particles and radiation through both spontaneous Hawking emission and superradiant scattering, whereas in the extreme case, the black hole can just radiate via superradiant scattering.
As one can see from Eq. (49),Ṁ andQ are related. Therefore, all of the above second moments can be reexpressed in terms ofQ. Let us calculateQ for our case. The rate of charge loss can be stated as where is the rate of electron-positron pair creation per fourvolume and e is charge of electron. According to Schwinger's theory [104] for (3 + 1) dimensions, the rate of electronpositron pair creation in a constant electric field E is where Q 0 = 4π eb 2z−2 /π m 2 l z−1 and m is the mass of electron. In the presence of linear Maxwell electrodynamics, the electric field is E = Q/r z+1 and therefore Combining Eqs. (55) and (56), we arrive at where [a, b] is incomplete gamma function. When r + Q, Eq. (57) reduces to where we have used in which x −1 1. In the following section, we turn to the study of thermodynamic geometry of our black hole solutions to figure out the behavior of the black hole's possible molecules and phase transitions.

Ruppeiner geometry
In this section we study thermodynamic geometry of the Lifshitz-dilaton black holes for linearly Maxwell and nonlinearly BI gauge fields, separately. We have introduced this method in Sect. 3.2.2 with focus on the study of the phase transitions which occur at divergence of the Ruppeiner invariant R. In addition to divergences, R has other properties which give us information about thermodynamic of the system. The sign of R gives us the information about the dominated interaction between possible black hole molecules, while its magnitude measures the average number of correlated Planck areas on the event horizon [16,[19][20][21]. R > 0 means the domination of repulsive interaction, R < 0 shows the attraction dominated regime and when R vanishes the system behaves like an ideal gas i.e. there is no interaction. In the following, we first study thermodynamic geometry in the presence of linear Maxwell electrodynamics. Then we extend our study to nonlinearly charged black holes where BI electrodynamics has been employed. There is just a necessary comment. As we stated before in Sect. 3 As we mentioned above, for investigating thermal stability we need to check the signs of M Q Q and H M S,Q . In our case Thus, in order to disclose the thermal stability of system, we need to study the sign of the determinant of the Hessian matrix. We find where The numerator N of (40) is a complicated finite function of S and Q in this case, including long terms that we do not express explicitly for brevity. However, as mentioned in Sect. 3.2.2, one can find the denominator D in the form of where T L M and H M S,Q have been give in (60) and (63), respectively.
Having Eqs. (63) and (65) at hand, we are in a position to investigate the divergences of R, which play the central role in thermodynamic geometry discussions, and also thermal stability of system. As one can see from Eqs. (63) and (65), for z = 2, the divergences occur just in the case of the extremal black holes where T L M = 0. For z = 2, in addition to extremal black hole case, R diverges in zeros of (63). In the latter case, we can calculate the corresponding temperature by solving F = 0 for Q and then putting this Q in Eq. (60) to arrive at The above temperature is negative for z > 2 i.e. there is no black hole at this diverging point and therefore the divergences of R occur just for extremal black hole case when z > 2. However, for z < 2 when T > 0, we can see an upper limit in entropy and charge of system. The largest entropy S for which F = 0 (which we call it critical entropy S c ) can be calculated by finding the extremum point where ∂F/∂ S = 0 as at which One should note that the absolute value of Q c is also the largest charge value which satisfies F = 0. Another remark to be mentioned is that (67) imposes an upper limit on the size of black hole too (see (21)). At this point, the corresponding temperature can be obtained: For charges greater than Q c , the Ruppeiner invariant diverges only in the case of extremal black holes. For Q = Q c , in addition to T L M = 0, we have one other divergence in R specified by (67) and (70). For Q < Q c , in addition to T L M = 0, we have at most two other divergences, since the order of the polynomial in terms of S is always lower than 3 for n ≥ 3 and z < 2. One should note that, in the latter case, the temperature region between two divergences is not allowed since H M S,Q < 0 (Fig. 5).
We have summarized the above discussion in Figs. 6 and 7. These figures also show the sign of the Ruppeiner invariant for different choices of the parameters that determines the type of interaction between black hole molecules [16,20,21]. Figure 6a is depicted for RN-AdS case (n = 3, z = 1). In this case, it can be seen that, for Q > Q c , the Ruppeiner invariant diverges only for extremal black holes. As Fig. 6a shows, there is also a range of T for which R < 0, namely the dominated interaction between black hole molecules is attractive. Furthermore, the interactions near zero temperature is the same as interactions of Fermi gas molecules near zero temperature [16]. According to Fig. 5a, for Q > Q c , H M S,Q is positive (also M Q Q > 0 (see (62))), and therefore the system is stable for all T region. For Q = Q c , in addition to the zero temperature, we have another temperature (T c ), where a divergence of R occurs (see Fig. 6a). At zero temperature, the Ruppeiner invariant goes to +∞, while at T c it goes to −∞. The latter case is similar to the Van der Waals gas phase transition at the critical point in the sense that at the phase transition temperature, R goes to −∞ [14,21]. For Q = Q c , R becomes positive when we go away from the second divergence (T c ) on the temperature axis. In Q = Q c case, H M S,Q is positive and just vanishes at T c (Fig. 5a), so, (a) (b) Fig. 7 The behavior of the Ruppeiner invariant R versus T for the linear Maxwell case with l = b = 1 the system is always thermally stable. For Q < Q c , there are three divergences; one at T = 0, one at T < T c and one at T > T c . In this case, according to Fig. 5a, H M S,Q is negative in the temperature region between two roots and show instability. This not-allowed region is equivalent to the temperature region between two divergences of the Ruppeiner invariant for Q < Q c (Fig. 6a). Figures 5b and 6b show the same properties for black holes with different parameters. In this case, T c is greater than one of previous case while Q c is lower. Figure 7  Thus, since H M S,Q nowhere vanishes for z > 2 (see discussions below (66)) and is positive at T = 0 according to above equation, it is positive throughout the temperature region and therefore the system is always thermally stable for z > 2.
Regarding the nature of the phase transition occurring at zero temperature where the Ruppeiner invariant diverges, we discussed in previous section via Landau-Lifshitz theory of thermodynamic fluctuations. However, regarding the phase transitions occurred at divergences of R at finite temperatures, we can give some comments here. We have seen two kinds of phase transitions here for z < 2 (see Fig. 6) namely continuous (for Q = Q c where R diverges at just one finite temperature or entropy) and discontinuous (for Q < Q c where R diverges at two finite temperatures or entropies and we have a jump between these two points since there is no thermally stable black hole between them). Both of these phase transitions can be considered as small/large black holes phase transitions. The first reason for this argument is that as temperature increases, entropy or equivalently size of black hole increases (note that ∂ S/∂ T = M −1 SS > 0). Therefore,on the left side of phase transition points where the temperature is lower, we have small size black holes and the right side where the temperature is higher we have large size ones. This fact can also be seen from the behavior of the Ruppeiner invariant magnitude at the two sides of the phase transitions. For small black holes, we expect a finite correlation between possible black hole molecules (of course far from phase transition points) because those are close to each other. For large black holes, we expect the correlation between possible molecules to tend to a small value near zero since molecules become approximately free. These expected behaviors can be seen in Fig. 6.

Born-Infeld case
For the Born-Infeld case, we can calculate the Ruppeiner invariant by using Eqs. (19), (34) and (39). The Ruppeiner invariant in this case is very complicated due to the presence of hypergeometric functions. Therefore, in this case we discuss the thermodynamic geometry non-analytically and by looking at plots. We study the cases z < 2, z > 2, z = 2 and z = n + 1 separately. First, we study the case z < 2. Figure  8a shows that changing β can cause a change in the dominant interaction. For instance, in a range of T , we have negative R (attraction) for β = 1 (note that in this range the system is thermally stable as one can see from Figs. 9a and 10a). For β = 1, the system behaves like a Fermi gas at zero temperature, namely R goes to positive infinity at zero temperature [16]. For β = 0.82, the Ruppeiner invariant diverges at two points; one of them is the zero temperature. According to Fig. 9a, for temperatures lower than the second divergence point, H M S,Q is negative and therefore the system is thermally unstable. Since M Q Q > 0 (Fig. 10a), the system is thermally stable just for temperatures greater than the temperature of second divergence for β = 0.82. Figure 8a shows Fig. 8b shows that the behavior of system looks like Van der Waals gas at phase transition temperature i.e. R goes to negative infinity at this point [14,21]. For z > 2, the behavior of R is depicted in Fig. 11a. It can be seen that the type of dominated interaction changes for different β and we have negative R for some cases. In this case, we have a behavior like a Fermi gas at zero temperature for extremal black holes. For β = 0.13, there is a divergence at non-zero temperature; for lower temperatures, the system is unstable (Fig. 11b, c) and which are always positive and therefore the system is always stable and R experiences no divergence (Fig. 12) which is positive for all temperatures. The behavior of the Ruppeiner invariant and H M S,Q are depicted in Fig. 13a and b for this case, respectively. As one can see the type of interaction is β-dependent for some temperatures. For β = 0.04, H M S,Q is positive just for temperatures greater than the finite temperature of divergence (Fig. 13b) and therefore the system is thermally stable for this range of temperatures.
Most of the phase transitions discovered above in the presence of BI electrodynamics at finite temperatures cannot be interpreted as small/large black hole phase transitions because in these cases small size black holes are unstable. Further studies to disclose the nature of these phase transitions are called for.

Summary and closing remarks
In many condensed matter systems, fixed points governing the phase transitions respect dynamical scaling, t → λ z t, x → λ x where z is the dynamical critical exponent. The gravity duals of such systems are Lifshitz black holes. In this paper, we first sought for the (n+1)-dimensional Born-Infeld (BI) charged Lifshitz black hole solutions in the context of dilaton gravity. We found that these solutions are different for the cases z = n + 1 and z = n + 1. We found that both solutions and showed that the solution for the case z = n + 1 can never be Schwartzschild-like. Then we studied the thermodynamics of both cases by calculating conserved and thermodynamical quantities and checking the satisfaction of the first law of thermodynamics. After that, we looked for the Hawking-Page phase transition for our solutions, both in the cases of linearly and BI charged black holes. We studied this phenomenon and the effects of different parameters on it by presenting the behaviors of temperature T with respect to entropy S at fixed electrical potential energy U and also Gibbs free energy G with respect to T . Then we turned to discuss the phase transitions inside the black holes. In this part, we first presented the improved Davies quantities that show the phase transition points coincided with the ones of the Ruppeiner geometry. This coincidence has been proved directly in the Appendix. All of our solutions, provided that those are thermally stable at zero temperature, show the divergence at this point both from the Ruppeiner and the Davies points of view. Using the Landau-Lifshitz theory of thermodynamic fluctuations, we showed that this phase transition is a transition of the radiance properties of black holes. At zero temperature, an extreme black hole can just radiate through superradiant scattering whereas a non-extreme black hole at finite temperature can give off particles and radiation via both spontaneous Hawking radiation and superradiant scattering.
Next, we turned to study Ruppeiner geometry for our solutions. We investigated thermal stability, interaction type of possible black hole molecules and phase transitions of our solutions for linearly and nonlinearly BI charged cases separately. For the linearly charged case, we showed that there are no diverging points for Ricci scalar of the Ruppeiner geometry (Ruppeiner invariant) at finite temperature for the case z ≥ 2. For z < 2, it was found that the number of divergences (which show the phase transitions) at finite temperatures depend on the value of the charge Q. We introduced a critical value for the charge, Q c ; for greater values there is no divergence at finite temperature, for values lower than it there are at most two divergences and for Q = Q c , there is just one diverging point for the Ruppeiner invariant. For the case of Q < Q c , there is a thermally unstable region for systems between two divergences at finite temperatures. So, this phase transition can be claimed as a discontinuous phase transition between small and large black holes. For small black holes not close to the transition point, we observed a finite magnitude for the Ruppeiner invariant R. This is reasonable since the magnitude of R shows the correlation of possible black hole molecules. Also, for large black holes the magnitude of the Ruppeiner invariant tends to a very small value as expected. For Q = Q c , the solutions show a continuous small/large black holes phase transition at finite temperature. In the case of BI charged solutions, we investigated the Ruppeiner geometry and thermal stability for z < 2, z > 2, z = 2 and z = n + 1 separately. In some of these cases, small black holes were thermally unstable. So, more studies are called for to discover the nature of phase transitions at diverging points of R. In both the linearly and nonlinearly charged cases, for some choices of the parameters, the black hole system behaves like a Van der Waals gas near the transition point.
Finally, we would like to suggest some related interesting issues which can be considered for future studies. It is interesting to repeat the studies such as regards the Hawking-Page phase transition, Ruppeiner geometry and Landau-Lifshitz theory for black branes to discover the effect of different constant curvatures of the (n − 1)-dimensional hypersurface on those phenomena. One can also seek for any signature of different phase transitions discovered here such as the Hawking-Page phase transition, and phase transitions determined by the Ruppeiner geometry, in dynamical properties of solutions by investigating quasi-normal modes. Some of this work is in progress by the authors.