Dilatonic black holes in gravity's rainbow with a nonlinear source: the effects of thermal fluctuations

This paper is devoted to investigate nonlinearly charged dilatonic black holes in the context of gravity's rainbow with two cases: I) by considering the usual entropy, II) in the presence of first order logarithmic correction of entropy. First, exact solutions of dilatonic Born-Infeld gravity with an energy dependent Liouville-type potential with black hole interpretation are obtained. Then, thermodynamic properties of the mentioned cases are studied, separately. It will be shown that although mass, entropy and heat capacity are modified due to the presence of first order correction, the temperature remains independent of it. Furthermore, it will be shown that divergencies of the heat capacity, hence phase transition points are also independent of first order correction whereas the stability conditions are highly sensitive to variation of the correction parameter. Except for the effects of first order correction, we will also draw a limit on values of the dilatonic parameter and show that it is possible to recognize AdS and dS thermodynamical behaviors for two specific branches of the dilatonic parameter. In addition, the effects of nonlinear electromagnetic field and energy functions on thermodynamical behavior of the solutions will be highlighted and dependency of critical behavior, on these generalizations will be investigated.


I. INTRODUCTION
Albeit extremely successful consequences of Einstein's theory of gravity, there are various reasons which signal that this theory should be modified, at least at high energy regime. Regarding the physical processes at energies of order of the Planck scale, one is expected to modify Einstein gravity to a quantum theory of gravity which is related to Quantum Field Theory and could be the correct framework to describe processes at ultraviolet (UV) energies. Considering the UV-completion regime of black holes, the gravitational interaction at energies exceeding the Planck mass implies that the UV-completion must be achieved by new quantum degrees of freedom of wavelength much shorter than the Planck length. However, it is expected to recover the semi-classical (thermodynamical) behavior of black holes in the mean-field approximation [1,2]. Horava-Lifshitz gravity [3] is one of the interesting UV modifications of general relativity (GR), in which reduces to GR at infrared (IR) limit. Based on Lifshitz scaling for space and time, different attractive subjects have been considered such as string theory of types IIA and IIB [4,5], AdS/CFT correspondence [6][7][8], dilatonic black holes/branes [9][10][11][12][13] and phase transitions/geometrical thermodynamics of black holes [14][15][16].
On the other hand, another point of view regarding the UV completion of GR is related to gravity's rainbow * email address: hendi@shirazu.ac.ir † email address: behzad.slampanah@gmail.com ‡ email address: shahram.panahiyan@uni-jena.de § email address: m.momennia@shirazu.ac.ir [17,18]. In order to introduce the gravity's rainbow, we can start from its historical origin; the doubly special relativity. It is known that the special theory of relativity is based on two postulates [19]; (i) the laws of physics are invariant in all inertial frames, and (ii) the constant invariant speed of light in the vacuum for all observers, regardless of the relative motion of the light source. This results into the fact that the speed of a massive particle cannot be equal or larger than the speed of light. Based on high energy point of view, if one consider an upper bound for the energy of a test particle and add this assumption as a postulate to the others mentioned above, the doubly special relativity will be constructed [20][21][22]. So, in this manner, the energy of test particle cannot be greater than the Planck energy. Following the gravity supplementing to the special relativity which leads to GR, one can generalize doubly special relativity in the presence of gravity to obtain the so-called doubly general relativity or gravity's rainbow [18].
Another way for constructing such theories is considering a modification of the standard energy-momentum relation. The modification of usual energy-momentum relation in gravity's rainbow is given as where E and E p are the energy of a test particle and the Planck energy, respectively and they are related to together with ε = E/E p in which ε ≤ 1, because the energy of a test particle E cannot be greater than E p [23]. The functions f (ε) and g (ε) are called rainbow functions. They are phenomenologically motivated and could be extracted by experimental data [24], see Refs. [17,25,26] for more details. It is worthwhile to mention that the rainbow functions are required to satisfy lim ε→0 f (ε) = 1 and lim ε→0 g (ε) = 1, where these conditions ensure that the modified energy-momentum relation reduces to its usual form in the IR limit. It is worth mentioning that such justification is based on the standard model of particle physics. Since quantum theories of gravity are different from classical one in high energy regime, a full quantum theory has to be valid in UV regime. It is notable that both of the Horava-Lifshitz gravity and gravity's rainbow are valid in the UV regime. On the other hand, considering suitable choice of the energy functions, the Horava-Lifshitz gravity can be related to the gravity's rainbow [27]. For the reasons mentioned, we consider the gravity's rainbow in this paper.
Recently, there has been a growing interest in gravity's rainbow because of its interesting achievements in the context of theoretical physics, such as providing possible solutions for information paradox [28,29], admitting usual uncertainty principle [30,31], existence of remnants for black holes after evaporation [32,33], and absence of black hole production at LHC [34]. In addition, from cosmological point of view, it is possible to remove the big bang singularity by using gravity's rainbow [35][36][37][38]. Black hole thermodynamics in the presence of gravity's rainbow coupled to (non)linear electrodynamics has been discussed in literature [39][40][41]. Moreover, the effects of rainbow functions on the thermodynamics properties and phase transition of black holes have been studied in Gauss-Bonnet gravity [42], Kaluza-Klein theory [43], massive gravity [44], and F (R) theories of gravity [45,46]. In the view point of astrophysics, the structure of neutron stars have been investigated and it is found that by increasing the rainbow function more than one, the maximum mass of neutron stars increases [47,48].
Regarding various aspects of black holes in the context of gravity's rainbow, it will be interesting to look for possible coupling of rainbow functions with scalar field and its self interacting potential. Astronomical observations indicated that our Universe has an accelerated expansion [49,50]. Einstein gravity cannot explain such acceleration. In this regards, various modifications of Einstein gravity are proposed in literature. Among them, scalar tensor gravity and string theory are most acceptable ones with interesting properties. It is worth mentioning that the low energy limit of effective string theory could lead to appearance of a dilaton scalar field [51,52] which motivates one to investigate GR coupled with such a scalar field. In addition, dilaton field coupled with other gauge fields has significant effects on the solutions [53][54][55][56]. In Ref. [57], it was shown that dilaton field can be an appropriate candidate for dark matter. In addition, in order to have deep picture regarding the nature of dark energy, a new scalar field can be added to the field content of original theory [58,59]. Also, GR coupled with a dilaton field can explain the accelerated expansion of our Universe, properly. Moreover, it was shown that in the presence of dilaton scalar field, the asymptotic behavior of black hole spacetimes is neither flat nor (a)dS [60][61][62][63][64]. This is due to the fact that in the limit r → ∞, the ef-fects of dilaton field is still present. Besides, it was shown that the dilatonic black hole solutions can be constructed in the (a)dS spacetime background [65]. Black objects in the presence of dilaton gravity and their thermodynamics have been studied in Refs. [66][67][68][69][70][71][72][73][74][75]. Recently, considering the dilaton field, the hydrostatic equilibrium equation of compact objects has been obtained and the properties of neutron stars have been analyzed in Refs. [76,77]. The evaporation of quantum black holes has also been investigated using two dimensional dilaton gravity [78,79]. In addition, there is a strong motivation to study rainbow deformation of geometries that occur in the string theory.
Nonlinear field theory is one of the most interesting branches in physical sciences because most physical systems presented in the nature are nonlinear. On the other hand, existence of some limitations in the linear Maxwell theory motivates one to consider nonlinear electrodynamics (NED) [80][81][82][83]. One of the main advantages of considering NED theories comes from the fact that these theories are richer than the linear Maxwell theory and in some special cases they reduce to the Maxwell field. Besides, it was shown that NED can remove the black hole and big bang singularities [84][85][86]. In addition, the effects of NED become quite important in superstrongly magnetized compact objects [87,88]. Studying general relativity coupled to NED is an attractive subject because of its specific properties in gauge/gravity duality. One of the most important NED theories is so-called Born-Infeld NED (BI NED) which has been introduced by Born and Infeld in 1934 [89]. This interesting type of NED removes the divergence point of electric field of a point-like charge and such a property was the main motivation for introducing this kind of NED. Other strong motivation for considering BI NED comes from the fact that it naturally arises in the low-energy limit of open string theory [90][91][92][93]. The following incomplete list shows the investigation of gravitational fields coupled to BI NED for static black holes [94][95][96][97][98][99][100][101], wormholes [102,103], rotating black objects [104,105], and superconductors [106].
Motivated by these applications, in this paper, we analyze dilaton field using the formalism of gravity's rainbow in the presence of BI NED. So in the present paper, we are going to study the dilatonic black hole solutions in Einstein gravity's rainbow coupled to BI NED with considering the effects of thermal fluctuations.

II. BLACK HOLE SOLUTIONS OF DILATON GRAVITY'S RAINBOW IN THE PRESENCE OF BI NED
Here, we are interested in dilatonic charged black holes in the context of gravity's rainbow with BI source. Since we are interested in minimal coupling, the action can be written as [107] where R is the Ricci scalar, and Φ and V (Φ) are, respectively, dilaton field and its corresponding potential. L(F, Φ) is the Lagrangian of the BI electromagnetic field under consideration with the following explicit form where the Maxwell invariant is denoted by . Also, β is the nonlinearity parameter and dilatonic constant is identified by α which is a parameter for determining the strength of coupling between the scalar field and electrodynamics.
For the sake of simplification in calculation, we use the following redefinition where L(Y ) = 1 − √ 1 + Y and Y = e −4αΦ F 2β 2 . Variation of the action (1) leads to the following field equations Here, we consider the following energy dependent line element in which R(r) and Ψ(r) are metric functions which should be determined. We recall that these two functions g(ε) and f (ε) are energy (or rainbow) functions that are chosen phenomenologically. In addition, dΩ 2 k is the metric of two dimensional subspace which depends on topology of the boundary of spacetime and could have positive (elliptic k = 1), zero (flat k = 0) or negative (hyperbolic k = −1) curvature such as Since we are interested in black holes with radial electric field, the suitable choice of gauge potential is which by using Eq. (6), we can obtain the following field equation where F tr = A ′ t is the tr-component of the electromagnetic tensor field and prime denotes d/dr. Solving Eq. (10), we obtain where q is an integration constant which is related to the total electric charge of the solutions. In order to find the metric functions, we consider the following modified version of the Liouville-type dilaton potential where Λ is a free parameter which plays the role of cosmological constant. It is worthwhile to mention that in IR limit (g(ε) = f (ε) = 1), Eq. (12) reduces to the known Liouville-type dilaton potential that is employed for finding Friedman-Robertson-Walker scalar field cosmologies [108] and Einstein-Maxwell-dilaton black holes [107,109].
Here, we consider R(r) = e αΦ(r) as an ansatz for finding the metric function. This ansatz is supported by studies conducted in Einstein-Maxwell-dilaton gravity [110]. Now, by using Eqs. (4), (5) and the metric (7) with obtained electromagnetic field tensor, one can find the following differential equation for calculating the metric function, analytically, where Θ i,j = i + jα 2 , η = q 2 f 2 (ε)g 2 (ε) 2γ and γ = 2α 2 /Θ 1,1 . Using obtained field equations, one can find both the metric function and dilatonic field in following forms where 1 4 , −η is the hypergeometric function. Also, b is an arbitrary constant related to the scalar field and m is an integration constant which is related to total mass of the solutions. It is worthwhile to mention that for limiting case of β −→ ∞, obtained metric function will lead to which is charged black hole in gravity's rainbow [111]. On the other hand, for limiting case of α = γ = 0 and β −→ ∞, our solutions reduce to which is the metric function of 4-dimensional asymptotically AdS topological charged black hole in gravity's rainbow [46]. In order to study the effects of matter fields as well as dilatonic gravity, we will investigate the behavior of the Kretschmann scalar for small and large values of radial coordinate. The existence of divergency for this scalar denotes that our solution contains a curvature singularity. If this singularity is covered with a horizon, obtained solution is interpreted as black hole. It is a matter of calculation to show that for this black hole, we have It is evident from Eq. (18) that there is an essential singularity located at the origin for this solution. Therefore, the first condition for having black hole is satisfied. On the other hand, the asymptotical behavior of system is modified due to dilatonic gravity and it is not (A)dS. It is notable that, the existence of horizon for the solution is investigated through following diagrams (Fig. 1).
Plotted diagram shows that depending on the choices of different parameters, obtained solutions may present black holes with two horizons, extremal black holes and naked singularity. For exampe, in the case of BI theory, for small values of nonlinearity parameter observed behavior is Schwarzschild like (continues line of Fig. 1). Increasing nonlinearity parameter will change the behavior of the metric function into Rissner-Nordstrom like which may yield two horizons (dotted line of Fig. 1), extremal black holes (dashed line of Fig. 1) or naked singularity (dashed-dotted line of Fig. 1). Here, we saw that generalization to nonlinear electrodynamics as well as gravity's rainbow, provided a richer phenomenologies for black holes. We see that geometrical structure of the black holes have been modified due to these generalizations. Depending on the choices of different parameters, the type of singularity, the general behavior of metric function and number and type of horizons have been modified. These modifications highlight the contributions of this matter field and gravities. We will continue to study the effects of these generalizations on conserved and thermodynamical quantities in the next section.

A. Usual thermodynamics (non-correction)
This section is devoted to the calculation of conserved quantities without thermal fluctuations. Here, we investigate the effects of gravity's rainbow, dilaton scalar field and nonlinear BI electrodynamics on thermodynamical quantities and check the validity of the first law of black hole thermodynamics. Using the concept of surface gravity, one can show that the Hawking temperature for obtained solutions will be where η + = η |r=r + . Since we are working in Einstein framework, the entropy of these black holes could be ob-tained by the area law in which by setting α = 0 and g(ε) = 1, the entropy of black holes in Einstein gravity is recovered. In order to find the total electric charge of these black holes, one can use the Gauss's law which leads to following result The electric potential of the black holes at the horizon radius with respect to spacial infinity as a reference could be calculated with the following relation where for nonlinear BI theory can be obtained as In order to obtain the total mass of black holes, one can use the definition of Abbott and Deser [112][113][114], which leads to It is worthwhile to mention that for IR limit (g(ε) = f (ε) = 1) and by setting α = 0, Eq. (25) reduces to the mass of Einstein-Maxwell black holes. Now, we are in a position to examine the validity of the first law of thermodynamics. To do so, first, we should obtain geometrical mass of the solutions, m, as a function of other parameters, by using Ψ (r = r + ) = 0. Then, by replacing it in Eq. (25), we obtain a relation for the total mass of the black holes versus entropy and electric charge. In this case, the extensive parameters will be the entropy and total electric charge and their corresponding conjugating quantities are the temperature and the electric potential, respectively. Therefore, the validation of the first law of black holes thermodynamics is done by It is a matter of the calculation to show that obtained quantities in Eq. (26) are exactly the same with those obtained in Eqs. (20) and (24). Therefore, we find that the first law is valid as Our final subject of the interest in this section is heat capacity. The information that are provided by this quantity could be used to render the thermodynamical structure of black holes in the context of their thermal stability/instability. In addition, the existence of discontinuity for this quantity signals the presence of thermal phase transition. It is worthwhile to mention that the discontinuities in the heat capacity usually are observed. For black holes in canonical ensemble, the heat capacity is calculated by following relation in which by using obtained temperature (20) and entropy (21), one can find it for this specific thermodynamical quantities as The thermodynamical stability and possible phase transition of the black holes in this case will be discussed in following sections. In the next section, we will obtain the thermodynamical quantities for the case where entropy is corrected to include a logarithmic correction.

B. Thermal fluctuations: correction
Here, we want to investigate the effects of thermal fluctuations on obtained solutions. Consideration of thermal fluctuations results into modifications of different ther-modynamical quantities, though some of the quantities remain fixed. To the leading order, the entropy gets logarithmic correction and given by in which S 0 is the uncorrected entropy which is given in the equation (21). Also, ζ is the thermal fluctuation parameter which we will call it correction parameter through the paper. Using Eqs. (20) and (21), we can obtain the corrected entropy as Our investigation regarding surface gravity confirms that the temperature is one of the thermodynamical quantities which is not affected by the presence of first order correction. In other words, for this case, the ob-tained temperature will be same as that was previously obtained (20). Using the corrected entropy (31) and temperature (20), we are able to compute the Helmholtz free energy as where Now, by using the obtained Helmholtz free energy and employing its relation with mass, M = F + T S, one can calculate the corrected mass in the following form We obtain the electric charge (Q) and the electric (chemical) potential (U ), which are similar to Eqs. (22) and (24), respectively. Therefore, by using the obtained thermodynamic quantities, the first law of black hole thermodynamics satisfy as Finally, the heat capacity of this case is of interest. It is a matter of calculation to obtain this quantity by using Eqs. (20) and (31) with (28) as Let us highlight some important properties of the ob-tained heat capacity in the presence of first order correc-tion. First of all, the effects of first order correction are observed only in the numerator of heat capacity which indicates that divergencies of the heat capacity, hence phase transition points are independent of first order correction. On the other hand, since the first order correction is present in numerator, it is expected that roots of the heat capacity and stability conditions are first order correction dependent. The contributions of first order correction and differences between correction and noncorrection cases will be discussed in details in next section.
C. Thermal structure: comparison between correction and non-correction In this section, the main goal is to understand the possible scenarios regarding thermal structure of the black hole solutions. Since we have obtained thermodynamical properties in the context of both first order correction and non-corrected solutions, we will give details in the context of both of them. Furthermore, we also investigate the details of contribution of the first order correction.
The obtained temperatures for correction and noncorrection cases are the same. On the other hand, from the first law of thermodynamics, we can see that temperature is calculated as a function of fluctuation in internal energy with changes of entropy. This shows that change in entropy due to variation of the correction parameter results into modifications in internal energy on a level that the temperature of this system remains fixed. In other words, modifications of entropy is realized by internal energy in a manner that temperature remains fixed. This provides us with a tool to increase/decrease the entropy and internal energy without any concern regarding the possible changes in temperature. This is in a manner, isothermal like behavior, although further tests are required to establish the isothermal nature of this property.
Returning to the obtained temperature, one can recognizes specific properties for it: I) obtained temperature provides us with imposing a limit on values that dilatonic parameter could acquire (α = 1). Such limitation prevents temperature to have divergent value, II) using this limit, one can see that essentially two branches exist for dilatonic parameter, hence temperature: α > 1 and α < 1. This is due to fact that the signature of different terms in temperature would be opposite for these two cases. Only exceptions are the cosmological constant and purely nonlinear terms. Later, we will show how these two branches give different pictures regarding the thermodynamics of these black holes, III) dilatonic parameter is coupled with all presented terms in the temperature. This shows that dilatonic gravity has profound effects on the behavior of temperature, hence thermodynamical structure of the black holes. Now, in order to investigate the behavior of temperature, we have plotted various diagrams for variation of dilatonic parameter, α. As it was pointed out, we have divided the effects of dilatonic parameter into two branches of α < 1 (left panel of Fig. 2) and α > 1 (left panel of Fig. 3).
First, we investigate α < 1 case. Evidently, depending on the choices of dilatonic parameter, temperature could have different properties such as: I) existence of one root and being increasing function of the horizon radius (continuous and dotted lines in left panel of Fig. 2). II) existence of one root and one extremum (dashed line in left panel of Fig. 2). III) existence of one root and two extrema (dashed-dotted line in left panel of Fig. 2). Evidently, for small values of the dilatonic parameter, the contribution of scalar field is not significant and general behavior of temperature is similar to the absence of this parameter. Increasing the dilatonic parameter leads to formation of an extremum. Further increment leads to existence of two extrema: a maximum and then a minimum in which maximum is formed in smaller horizon radius comparing to minimum. The extrema in temperature are matched with divergencies in the heat capacity. Therefore, one can conclude that in the case of α < 1, increasing the dilatonic parameter leads to introduction of critical behavior in thermodynamical structure of the black holes. Later in studying the heat capacity, we will give more details regarding the types of critical behavior that could exist for these black holes.
Let us now turn our attention to the case of α > 1. Interestingly, for small values of the dilatonic parameter, the temperature is negative valued everywhere (continuous line in left panel of Fig. 3). This indicates that our solutions are thermodynamically non-physical. It is worthwhile to mention that for this case, temperature enjoys a maximum in its structure. The maximum value of temperature is an increasing function of dilatonic parameter. For specific value of dilatonic parameter, one can find a root for temperature (dotted line in left panel of Fig. 3). The root is an extreme point, but then again, we should point it out that except for the root, temperature is negative everywhere. Increasing dilatonic parameter more than this specific value leads to formation of two roots for the temperature (dashed and dashed-dotted lines in left panel of Fig. 3). The maximum is located between these roots, which shows that the physical black holes could only be observed for medium black holes. Whereas, for small and large black holes, temperature is negative valued and solutions are thermally non-physical.
The behavior that we observed in plotted diagrams for α < 1 case actually shows the existence of subcritical isobars. Presence of subcritical isobar, so far has been reported only for AdS black holes. On the other hand, the behavior that we observed for temperature in α > 1 case was similar to the one that previously was reported for dS black holes. These two specific properties confirms a very important result regarding the dilatonic parameter: for α < 1 case, thermodynamical behavior or temperature of black holes is AdS like while for α > 1 case, this behavior is dS like. This indicates that in general, for α < 1, the black holes have AdS like behavior while for the α > 1 the behavior is dS like. If one is interested to study the AdS/CFT duality in the context of these black holes, the valid branch is where α < 1. Now, we turn our attention to the heat capacity which contains information regarding thermal stability and possible phase transition points. Unlike temperature, the heat capacity could be correction dependent. In other words, considering the first order correction to the entropy, the heat capacity would be modified. Interestingly, the effects of correction could be observed in the numerator of heat capacity while its denominator is independent of it. This confirms two important things regarding the effects of first order correction: I) the phase transition points that are realized through divergencies in the heat capacity are not affected by the first order correction. In other words, the first order correction has no effect on phase transition points, II) in usual black holes without the first order correction, the heat capacity and temperature share same roots. But since the numerator of heat capacity in the presence of first order correction is modified, temperature and heat capacity no longer share the same roots. This modifies the conditions regarding thermal stability of the solutions.
We recall that in the case of α < 1, the behavior is AdS like. For this case, there exists a critical value for dilatonic parameter, say α c (c stands for "critical") which could be used to divide possible scenarios available for the heat capacity. For 0 < α < α c , the effects of dilatonic gravity is significant on the place of root of heat capacity. Before the root, the heat capacity is negative and solutions are thermally unstable while after it, the opposite is seen and solutions are thermally stable (see continuous and dotted lines in middle panel of Fig. 2). This shows that for small values of the dilatonic parameter, region of stability is modified. At the critical value (α = α c ), the heat capacity acquires a divergency which is interpreted as critical behavior (dashed line in middle panel of Fig. 2). The sign of heat capacity around this divergence point is positive which indicates that the crit- ical characteristic takes place between two stable phases, as expected. In this case, the heat capacity enjoys a root as well, and it is located before the divergency. Finally, for α c < α < 1, the heat capacity has two divergencies in its structure (dashed-dotted line in middle panel of Fig.  2). Between the divergencies, the sign of heat capacity is negative which shows that solutions are thermally unstable. This indicates that there is a phase transition taking place between two divergencies. Therefore, the only stable phases provided for the black holes in this case are small and large black holes while the medium black holes suffers thermal instability. It is worthwhile to mention that before the smaller divergency, there also exists a root for the heat capacity. Before this root, the temperature is negative, therefore solutions are non-physical. The behavior that we observed for the heat capacity here completely matches to the one that was observed for the heat capacity of AdS black holes (see appendix for more details). This confirms that thermal stability structure of these black holes in case of α < 1 is same as AdS black holes. This provides us with further proof to recognize the branch α < 1 as AdS spacetime.
In case of α > 1, the general behavior in temperature was dS like. Here, it is possible to divide the general behavior of the heat capacity into three groups by a specific value of α, say α er (er stands for extreme root). For 1 < α < α er , the heat capacity has only one divergency in which the heat capacity switches from negative to positive (continuous line in middle panel of Fig. 3). This case has negative temperature everywhere. Therefore, although the heat capacity signals the existence of stable state, the negative temperature shows that no physical solution exists in this case. This highlights the importance of studying the temperature alongside of the heat capacity to separate physical solutions from non-physical ones. For α = α er , interestingly no divergency is observed for the heat capacity, although the temperature gives us the detail of its existence by having a maximum (dotted line in left and middle panels of Fig. 3). In this case, the temperature is negative valued everywhere except at its root which is an extreme one. But the heat capacity in this case shows the existence of only one root which after it, the heat capacity is positive valued. The absence of divergency in heat capacity is due to the fact that extremum and root of the temperature are identical. Considering the relation for obtaining heat capacity (28), it is obvious that no divergency could be observed for the heat capacity in this case. Finally, for α er < α, the heat capacity enjoys two roots with a divergency located between them (dashed and dashed-dotted lines in middle panel of Fig. 3). Between smaller root and divergency, and after larger root, the heat capacity is positive. Whereas, before smaller root, and between divergency and larger root, the heat capacity is negative and solutions are thermally unstable. But we should remind the reader that the only physical region for this case is between two roots (positive T ). Therefore, there is one physical stable phase and a physical unstable one around divergency and small black holes are stable. The behaviors that we have observed for the heat capacity in this case (α > 1) is exactly the same as the one that was observed for dS case (see appendix for more details). Therefore, this confirms the analogy of consideration of α > 1 as dS case. Now, we give more details regarding the effects of first order correction on thermodynamical behavior of the solutions. To do so, we have considered different behaviors that were reported for the heat capacity before and study the effects of variation of the first order parameter, ζ.
In brief, we should note that for α < 1, we had three cases for the heat capacity: I) existence of only one root (0 < α < α c ); for this case, the variation of correction parameter leads to modification in place of the root which is an increasing function of ζ (up-left panel of Fig. 4), II) existence of one root and one divergency (α = α c ); for small values of the first order correction parameter, ζ, the place of root is modified to larger values, but the sign of heat capacity around divergence point remains positive. Interestingly, for sufficiently large values of the correction parameter, the place of root will be shifted to after the divergency and the sign of heat capacity before it and also around divergence point will be negative (up-middle panel of Fig. 4). Therefore, the only thermally stable phase is after the root and around the divergence point, phase transition is between two unstable black holes, III) Existence of one root and two divergencies (α c < α); in this case, for small values of the correction parameter, the place of root is modified to larger values of the horizon radius but before the divergencies. The stability conditions remain unchanged. But for medium values of the correction parameter, the root will be located between the divergencies (up-right panel of Fig. 4). Interestingly, before smaller divergency, and between the root and larger divergency, the heat capacity is negative. Whereas, between the root and smaller divergency, and also after the larger divergency, the heat capacity is positive (dashed and dashed-dotted lines in up-right panel of Fig. 4). Increasing the value of correction parameter leads to root being place after the larger divergency. In this case, the only stable phases (positive heat capacity) are between the divergencies and after the root (bold continuous and bold dotted lines in up-right panel of Fig. 4).
For α > 1, we only consider the physical case in which the temperature has two roots and between them, temperature has positive values. Evidently, contribution of the first order correction results into modification of the place of root for the heat capacity to higher values of the horizon radius. The root will be located after divergency (up panels of Fig. 5). Before root (around divergence point as well) the heat capacity is negative. The positive heat capacity could be observed after its root. But the root of heat capacity is located after larger root of the temperature, therefore, even though the heat capacity changes its sign to positive after its root, but it is within non-physical region. This results into a conclusion that existence of the first order correction for this case leads to instability of the black holes within physical region. We see that contribution of the first order correction destabilizes the solutions.
Our next study here is measuring the modifications in entropy by studying the fluctuation in temperature. To do so, we have plotted diagrams for both cases of α < 1 and α > 1 corresponding to those plotted for temper-ature and heat capacity. In addition, we have plotted other diagrams to understand the effects of first order correction on thermodynamical behavior of the entropy as a function of temperature.
For α < 1 which is interpreted as AdS case, evidently, three distinctive behaviors could be observed for the entropy. These behaviors could be divided and addressed by α c which was introduced before in the context of heat capacity. For 0 < α < α c , the general behavior of the entropy is not modified on a significant level. Here, entropy is an increasing function of temperature (continuous and dotted lines in right panel of Fig. 2). By setting α = α c , an extremum is formed for entropy versus temperature diagrams (dashed line in right panel of Fig. 2). This case corresponds to the cases where the temperature and the heat capacity enjoy an extremum and divergency, respectively, in their structures. Therefore, one can conclude that the entropy of extremum point is the critical entropy where system goes under thermal phase transition. Increasing the dilatonic parameter to reach the region of α c < α results into formation of two extrema which could be recognized by T 1 and T 2 (dashed-dotted line in right panel of Fig. 2). Between these two indications, the temperature is a decreasing function of entropy and for every temperature, there exists three different entropies. For T = T 1 , T 2 every temperature has two specific entropies. This case corresponds to the existence of two divergencies in the heat capacity. Thermodynamical principle informs us that system is in favor of increasing its entropy. This indicates that for the case where three (two) entropies are available for temperature, system moves to the case, which has the highest entropy. This is indeed the characteristic behavior of the phase transition. Therefore, we see that measurement of the entropy as a function of the fluctuation of temperature provides us specific characteristic that enables us to recognize critical behavior.
For α > 1 which is interpreted as dS case, we see that for non-physical case where system has the negative temperature, the entropy is positive valued (continuous line in right panel of Fig. 3). Increasing dilatonic parameter leads to temperature acquires an extreme root where, interestingly, entropy has positive and non-zero value (dotted line in right panel of Fig. 3). Increasing the dilatonic parameter furthermore leads to formation of a region of positive temperature with a maximum provided for temperature (dashed and dashed-dotted lines in right panel of Fig. 3). In this case, except for the maximum, every temperature has two different entropies. The maximum of temperature is where the system acquires divergency in its heat capacity. Now remembering that system de-sires higher values of the entropy, we can see that characteristic phase transition behavior is observed by system jumping from smaller entropy to larger one at the same temperature. Interestingly, even for vanishing temperature, there are two values provided for the entropy of system. Now, let us focus on the effects of first order correction. As it was pointed out, the case α < 1 (AdS case) enjoys larger number of the possibilities in its structure depending on the value of dilatonic parameter. To have more conclusive discussion regarding the results of first order correction, we study the effects of variation of the correction parameter for different cases provided α < 1 and for the physical case observed in α > 1 separately; I) for 0 < α < α c in case of α < 1, we see that adding first order correction results into formation of a minimum for entropy, S min (down-left panel of Fig. 4). In addition, for vanishing temperature, entropy is positive and non-zero valued, S 0 . In the region S min < S < S 0 , for every entropy, there exists two different temperatures. The minimum and S 0 are increasing function of the correction parameter, II) for α = α c in case of α < 1, where an extremum was observed in the absence of correction term, adding the correction results into formation of a minimum for the entropy alongside the extremum that was observed before (down-middle panel of Fig. 4). But interestingly, the entropies of minimum and extremum points are increasing functions of the correction parameter. Here too, similar to previous case, for specific region of entropy, every entropy has two different temperatures. Now, remembering that extremum point is where system goes under phase transition, one can conclude that entropy of the critical point depends on value of the correction parameter. This shows that critical behavior in the presence of first order correction is reached in higher values of the disorder in system, hence entropy, III) for α c < α in case of α < 1, as it was pointed out, two distinctive temperatures where available, T 1 and T 2 , in which between them, every temperature has three different entropies. In addition, the entropy related to T 1 is larger than entropy corresponding to T 2 (continuous line in down-right panel of Fig. 4). Now, the effects of first order correction in this case could be divided into four groups which are recognized with ζ 1 , ζ 2 and ζ 3 in which ζ 1 < ζ 2 < ζ 3 . For 0 < ζ < ζ 1 , the general behavior of the entropy versus temperature is same as that observed for vanishing T with one difference that entropies corresponding to T 1 and T 2 are higher than vanishing temperature (dotted line in down-right panel of Fig. 4). Increasing the correction parameter to reach the region of ζ 1 < ζ < ζ 2 , interestingly, results into formation of cycle for entropy versus temperature (dashed line in downright panel of Fig. 4). This means that between T 1 and T 2 , there is yet another temperature, T ′ that has two different entropies corresponding to it. In this case, the entropy of T 1 is still bigger than entropy of T 2 . If we choose the correction parameter from the region of ζ 2 < ζ < ζ 3 , the cyclic behavior could be observed, but interestingly, the entropy corresponding to T 2 is larger than the one related to T 1 (dashed-dotted line in down-right panel of Fig. 4). This means that although the entropies related to these two temperatures are increasing functions of the correction parameter, but the entropy related to T 2 grows faster comparing to the one corresponding to T 1 . Finally, for the region of ζ 3 < ζ, the cyclic behavior is vanished, but since the entropy of T 2 is now bigger than the one related to T 1 , the diagrams has opposite curve comparing to small values of the correction parameter (bold continuous and dotted lines in down-right panel of Fig. 4), IV) finally, for the physical case observed in α > 1, in the absence of correction, we observed a maximum for temperature. Interestingly, the existence of first order correction results into a cyclic (closed) diagram for the entropy versus temperature (down panels of Fig. 5). In addition, there is a minimum obtained for the entropy which has maximum temperature as its correspondence. In this case, for every temperature (entropy) there exists two entropies (temperatures). The minimum of entropy is an increasing function of the correction parameter.
In order to complete our study, we finally investigate the effects of nonlinearity parameter and energy functions on thermodynamical behavior of these black holes for two branches of α < 1 (AdS case) and α > 1 (dS case).
First, we turn our attention to the nonlinear nature of solutions. For case of α < 1, evidently, the effects of nonlinear electrodynamics could be divided into three categories with a critical value for nonlinearity parameter, β c (in which c stands for critical). In the absence of nonlinearity parameter, temperature has a minimum and it is positive valued everywhere (continuous line in up-left panel of Fig. 6). In the presence of nonlinearity parameter and for the region of β < β c , temperature enjoys a root, a maximum and a minimum in its structure (dotted line in up-left panel of Fig. 6). As for β = β c , the number of extremum is reduced to one which is located after the temperature (dashed line in up-left panel of Fig. 6). Increasing the nonlinearity parameter beyond β c < β results into vanishing the extremum in temperature, and temperature will be an increasing function of the horizon radius with a root (dashed-dotted line in upleft panel of Fig. 6). The corresponding heat capacity diagram shows that temperature and heat capacity share same roots and the extrema of temperature are matched with divergencies in the heat capacity (compare up-left and up-right panels of Fig. 6). Considering this fact, one can conclude two important points for α < 1 branch; I) for small values of the nonlinearity parameter, thermodynamical structure of the black holes enjoys the existence of phase transition, II) increasing the nonlinearity parameter results into absence of the phase transition for these black holes. It is worthwhile to mention that for vanishing the nonlinearity parameter, since temperature has a minimum, heat capacity has a divergency, which signals for possible critical behavior. The sign of heat capacity changes from negative to positive around divergency which indicates a phase transition between small unstable black holes to large stable ones. For case of α > 1, in the absence of the nonlinearity parameter, temperature is a decreasing function of the horizon radius with a root (continuous line in down-left panel of Fig. 6). In the presence of nonlinearity parameter, temperature acquires a maximum with two roots. The maximum and number of the roots are decreasing functions of the nonlinearity parameter. This indicates that for sufficiently large values of the nonlinearity parameter, the maximum will be relocated to negative values and roots will be vanished (dashed-dotted line in down-left panel of Fig. 6). Therefore, one can conclude that for α > 1 branch, the effects of increasing nonlinearity parameter result into elimination of physical solutions. As for the stability, interestingly, in the absence of nonlinearity parameter and in the region where temperature is positive, the heat capacity is negative, hence solutions are thermally unstable. After the root, although heat capacity is positive, temperature is negative valued. Therefore, for this case, only unstable solutions exist. In the presence of nonlinearity parameter, if the maximum of the temperature is located at positive values (leads to presence of two roots for temperature), heat capacity enjoys a phase transition between large unstable to small stable black holes (dotted and dashed lines in down-right panel of Fig. 6).
The situation for the effects of the energy functions is investigated in Fig. 7. Here, we have investigated the effects of gravity's rainbow for both branches of α < 1 (up panels of Fig. 7) and α > 1 (down panels of Fig. 7). Evidently, for small values of the rainbow functions, temperature enjoys a root and a maximum and a minimum in its structure which correspondingly, heat capacity would have the same and two divergencies matching extrema points in temperature (continuous line in up panel of Fig. 7). This shows that these black holes, for small values of the energy functions enjoy a phase transition over a region which is located between two divergencies of heat capacity. The number of the extremum in temperature, hence divergency in heat capacity is a decreasing function of the energy function. For certain value of the energy function, temperature will have an extremum and heat capacity enjoys a divergency (dashed lines in up panel of Fig. 7). In this case, a phase transition between two stable black holes takes place. Increasing the energy function beyond this value results into the absence of extremum in temperature, hence divergency in the heat capacity. Therefore, existence of the critical behavior depends on values of the energy functions for α < 1 branch. As for α > 1, it is evident that maximum exists for temperature which corresponds to a divergency in heat capacity (down panels of Fig. 7). For small values of the energy function, this maximum is within positive valued region and temperature also enjoys existence of two roots. The maximum is a decreasing function of f (ε) and for sufficiently large values of this energy function, it will be relocated into negative region indicating negative temperature without root, hence absence of physical solutions. Therefore, one can conclude that existence of physical solutions for α > 1, depends on the values that energy functions can acquire.

IV. CONCLUSION
Regarding the fact that we should consider high energy (UV) regime near the black holes motivates us to consider an energy dependent spacetime with a minimal coupling between gravity, dilaton scalar field as well as a nonlinear U (1) gauge field. In this paper, we have studied the black hole solutions in dilaton gravity's rainbow in the presence of BI source. It was shown that the nonlinearity parameter, energy functions and dilatonic parameter modify the type of the singularity, the places and number of the horizon radius and their type as well. These modifications highlighted the importance of matter fields and gravities.
Next, we obtained conserved and thermodynamical quantities and we proved that despite the mortifications of the gravity's rainbow, dilatonic gravity and BI field, the first law of the black holes thermodynamics is valid.
Our study in the context of temperature and heat capacity provided us with interesting properties for dila-tonic parameter. First of all, we were able to impose an specific limit on dilatonic parameter in order to avoid divergent temperature. In addition, it was shown that for α < 1, the characteristic behaviors of temperature and heat capacity are exactly the same as those that were observed for black holes in AdS spacetime. Whereas for α > 1, the properties extracted for the temperature and the heat capacity matches those observed for black holes in dS case. These two specifications provide us with the possibility of conducting studies that are specified for the branch AdS/dS such as AdS/CFT correspondence.
We have investigated the effects of thermal fluctuations on thermodynamical quantities. Although the entropy and mass were affected by thermal fluctuation, the first law remained valid in this case as well. In studying the effects of first order correction, it was shown that although the phase transition points are independent of first order correction, the stability conditions, which determine thermal structure of the black holes, are highly sensitive to variation of the first order correction parameter. Here, we observed that depending on the choice for correction parameter, the stability regions for black holes are modified on high level depending on the possible scenarios provided for α < 1. On the other hand, the contribution of first order correction resulted into destabilization of the solutions for α > 1 case. Remembering that α < 1 could be interpreted as AdS spacetime and α > 1 as dS case, one can conclude that for the AdS case, the stability conditions depend on the correction parameter while for dS case, first order correction makes the solutions unstable within physical region. Furthermore, one can observe that generalization to include the first order correction, results into larger families of thermal structures provided for these black holes in AdS case.
The measurement of entropy as a function of the fluctuation in temperature revealed profound deviation from the non-correction case. In α < 1, it was shown that depending on choices of the correction parameter, it was possible to introduce phenomena such as cyclic behavior and transition from one specific curve for diagrams to opposite curve. As for α > 1, existence of first order correction resulted into cyclic behavior as well for entropy versus temperature diagrams and a minimum for entropy which has maximum of temperature as its correspondence. Generally speaking, the possibilities of such behaviors were provided due to the fact that entropy is correction dependent while the temperature is independent of it. In addition to mentioned important effects of the first order correction, there is another effect that is of importance; the temperatures of the critical point through entropy versus temperature diagrams remained fixed but the critical entropies where shifted to higher values. This confirms that generically, the black holes with first order correction included, go under phase transition in higher level of disorder in their system.
To complete our thermodynamical investigation, we studied the effects of the nonlinearity parameter and energy functions on temperature and heat capacity. In α < 1 case, interestingly, it was observed that for small values of the energy function and nonlinearity parameter, system enjoys a phase transition over a region. Increasing these quantities leads to formation of an extremum for temperature, hence a divergency for heat capacity. In this case, the phase transition was taking place over a single point. As for α > 1, the maximum for temperature was a decreasing function of the energy functions and nonlinearity parameter, and for sufficiency large values of these quantities, maximum will be in negative region of temperature which is interpreted as absence of physical solutions. It is worthwhile to mention that for vanishing nonlinearity parameter, the general behavior of temperature and type of divergencies where completely different indicating the high contribution of the nonlinear electromagnetic field generalization.
Our study in this paper confirmed that consideration of the first order correction, results into larger classes of thermodynamical structure provided for black holes. But there are two questions that should be answered: I) are we free to choose any value for the correction parameter? II) are all possibilities that are provided for the thermal structure of the black holes, physical ones? One method to regard these questions is through the extended phase space. The concept of extended phase space provides the possibility of studying specific properties for black holes such as compressibility coefficient and speed of sound. These properties along with the concept that speed of sound must not exceed the speed of light, provide us with the possibility of finding upper/lower limit on values that correction parameter can acquire and distinguish physical thermodynamical structures from non-physical ones that were obtained in this paper. The only shortcoming of this method is the fact that extended phase space could only be introduced for AdS case. In any case, this issue is now under investigation and we hope to address these questions in another paper.