Topological black hole in the theory with nonminimal derivative coupling with power-law Maxwell field and its thermodynamics
M. M. Stetsko

Nanorobots-assisted Detection of Multifocal Cancer: A Multimodal Optimization Perspective
    Shaolong Shi, Student Member, IEEE,
        Yifan Chen, Senior Member, IEEE,
        and Xin Yao, Fellow, IEEE
S. Shi is with the Harbin Institute of Technology, Harbin 150001, China,
and also with the Department of Electrical and Electronic Engineering,
Southern University of Science and Technology, Shenzhen 518055,
China (e-mail: shishaolong@hit.edu.cn).
Y. Chen is with the Faculty of Science and Engineering and the Faculty
of Computing and Mathematical Sciences, University of Waikato, Hamilton
3240, New Zealand, and also with the Department of Electrical and
Electronic Engineering, Southern University of Science and Technology,
Shenzhen 518055, China (e-mail: yifan.chen@waikato.ac.nz).
X. Yao is with the Department of Computer Science and Engineering,
Southern University of Science and Technology, Shenzhen 518055,
China, and also with the School of Computer Science, University of
Birmingham, Birmingham B15 2TT, U.K. (e-mail: yaox@sustc.edu.cn).
    December 30, 2023
=====================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================



We obtain topological black hole solutions in scalar-tensor gravity with nonminimal derivative coupling between scalar and tensor components of gravity and power-law Maxwell field minimally coupled to gravity. The obtained solutions can be treated as a generalization of previously derived charged solutions with standard Maxwell action <cit.>. We examine the behaviour of obtained metric functions for some asymptotic values of distance and coupling. To obtain information about singularities of the metrics we calculate Kretschmann scalar. We also examine the behaviour of gauge potential and show that it is necessary to impose some constraints on parameter of nonlinearity in order to obtain reasonable behaviour of the filed. The next part of our work is devoted to the examination of black hole's thermodynamics. Namely we obtain black hole's temperature and investigate it in general as well as in some particular cases. To introduce entropy we use well known Wald procedure which can be applied to quite general diffeomorphism-invariant theories. We also extend thermodynamic phase space by introducing thermodynamic pressure related to cosmological constant and as a result we derive generalized first law and Smarr relation. The extended thermodynamic variables also allow us to construct Gibbs free energy and its examination gives important information about thermodynamic stability and  phase transitions. We also calculate heat capacity of the black holes which demonstrates variety of behaviour for different values of allowed parameters.


§ INTRODUCTION

General Relativity is very renowned theory that gives explanations for numerous facts and is in agreement with lots of observations  starting from planetary scales and up to the cosmological ones <cit.>. One of the biggest prediction of General Relativity that was made in early years of this theory, namely gravitational waves, finally have been discovered experimentally <cit.>  after the century of great ideas but unsuccessful attempts to detect them. Despite the great success of General Relativity there are still several open issues which are extremely important in order to have more complete picture of our Universe. Among those puzzles which remain unsolved yet we distinguish the problems about the origin of cosmological singularities, cosmological constant problem and closely related to it dark energy and dark matter issues and also the question about the presence of terms of higher order in curvature than the standard Einstein-Hilbert term. All the mentioned questions bring the idea of modification of General Relativity <cit.>, but there are a lot of possible ways of this generalization and some of them have received great interest and have been applied to numerous problems of Gravity and Cosmology, the others have not become so popular and known. Systematic description of open issues and the present day challenges of General Relativity and the approaches which face these challenges and give reasonable solutions for them is given in quite recent review <cit.>. 

Amidst these generalizations of standard General Relativity that have got broad interest are the so called scalar-tensor theories, which can be classified by the type of coupling between gravitational tensor and scalar degrees of freedom. We would also like to distinguish Horndeski gravity which appears to be the most general scalar-tensor theory where the equations of motion for tensor as well as scalar components are of the second order <cit.>.  Being almost unknown for more than three decades it has been studied intensively since the early years of the last decade. This interest to some extent was motivated by the examination of the so called Galileon theories which are a sort of the scalar-tensor theories with a specific Galilean invariance regarding the scalar field <cit.>.  It was shown that Galileon theory and Horndeski gravity are equivalent <cit.>. Galileon action was also obtained with help of the standard Kaluza-Klein procedure of compactification of higher dimensional Lovelock theory <cit.>. It should be pointed out that here dilaton fields which appear due to the use of the Kaluza-Klein procedure takes the role of the Galileons. We also note that some particular cases of general Horndeski theory, namely the theory with nonminimal derivative coupling, were derived in low energy String Theory <cit.>. Galileon theory also appears in the framework of massive gravity <cit.>. It should be pointed out that more general multiscalar versions of Horndeski theory were formulated <cit.> and theories beyond Horndeski were considered <cit.>.

Horndeski gravity and its particular cases have been applied to numerous problems of cosmology and black holes' physics. Namely, exact cosmological solutions within the nonominimal derivative coupling theory were examined <cit.>. Cosmological scenarios including quintessence and phantom fields were investigated and transition between different types of de Sitter solutions was shown <cit.>. Models of dark energy/dark matter in the framework of nonminimally coupled gravity and scalar field were considered and their dynamics was examined <cit.>. Influence of power-law scalar potential on cosmological dynamics was investigated <cit.>. New inflationary mechanisms which were supposed to be more realistic were examined <cit.>. Nonminimally coupled scalar-tensor theory can successfully generate slow rolling inflation without violation of unitarity bounds <cit.>. Reheating process during rapid oscillations of the inflation and curvaton scenario in nonminimal derivative coupling models were studied <cit.>. Reconstruction of Horndeski gravity from effective field theory of dark energy was also investigated <cit.>. Constraints on Horndeski theory and its predictions in the light of modern gravitational observation were also discussed <cit.>. 

Black holes and other compact astrophysical objects are also very interesting area of investigation. Taking into account quite recent observation of  gravitational waves and the fact that the sources of the waves are the compact objects they become of extreme importance for the verification of predictions of different generalizations of standard of General Relativity, including Horndeski theory, its generalizations or particular cases. Due to complicated structure of general Horndeski theory the special interest was devoted to the examination of black holes in some particular cases of Horndeski theory. Namely, a static black hole solution in the model with nonminimal derivative coupling but without cosmological constant was obtained <cit.>. Four and five-dimensional solutions in the the theory with cosmological constant was found and some aspects of its thermodynamics were investigated <cit.>. More general multidimensional and topological solutions were obtained <cit.>. Stealth Schwarzschild and de Sitter self-tuned solutions were found <cit.>. Slowly rotating neutron stars and black holes in the nonminimal derivative coupling sector of Horndeski gravity were examined in <cit.>. Three dimensional solutions were considered in <cit.>. Some aspects of black hole thermodynamics were studied <cit.>.  Black holes in more general Lovelock Galileon theory were investigated <cit.>. Stability of black holes in Horndeski gravity was paid attention to in the papers <cit.>. Causal structure of black holes in shift-symmetric Horndeski theories was examined <cit.>. Some aspects of neutron and boson stars in Horndeski gravity were considered <cit.>. The existence of black hole's hair in Horndeski gravity was studied <cit.>. Hairy black holes in the theory with nonminimal coupling and Gauss-Bonnet term were investigated <cit.>. 

As it has been noted above charged topological black holes in Einstein-Horndeski gravity with standard (linear) Maxwell action term were considered and their thermodynamics was studied in the paper <cit.>.  In our work we examine topological black holes in the theory with nonminimal derivative coupling between tensor and scalar components and nonlinear Maxwell field with Lagrangian of power-law type. The obtained in our work black holes' solutions can be considered as generalization of previously mentioned solutions <cit.> and when the factor (power) of nonlinearity p=1  those previously obtained solutions <cit.> are recovered. It is known that standard Maxwell Lagrangian in a space-time of dimension other than four is not conformally invariant and its generalization to power-law nonlinear expression allows to recover the conformal invariance <cit.>. We point out here that Maxwell term with the power-law nonlinearity in Lagrangian was considered in different works related to black holes <cit.>. Having obtained the black hole solutions we investigate their thermodynamics, namely we calculate temperature, entropy, Gibbs free energy and heat capacity and study their behaviour. We note that some keypoints of the present work were discussed in the work <cit.>.

This work is organized as follows. In the second section we derive topological black hole solutions and investigate the behaviour of metric functions for different distances and regimes of coupling, we also examine Kretschmann scalar in order to make some conclusions about singular points of the obtained solution, we also obtain gauge potential and derive some constraints on parameter of nonlinearity in order to have this potential reasonable. The third section which is divided on several subsections is devoted to the study of thermodynamic properties of the obtained black  holes. Namely using standard geometric relation we calculate black hole's temperature, then with help of Wald procedure we obtain a relation for entropy and it allows us to construct the first law of black hole's thermodynamics. Using extended thermodynamics approach we calculate Gibbs free energy and derive the generalized first law and Smarr relation, we also investigate heat capacity of the black holes. Finally, the forth section contains some conclusions.



§ EQUATIONS OF MOTION FOR THE SYSTEM WITH NONLINEAR MAXWELL FIELD

We consider an action for a system which can be described by a standard Einstein-Hilbert term with cosmological constant and some scalar field with nonminimal derivative coupling and we also suppose that the additional Maxwell field we take into account is represented by a nonlinear term. So, the action integral can be written in the form:

    S=∫ d^n+1x√(-g)( R-2Λ-1/2(α g^μν-η G^μν)∂_μ∂_ν +(-F_μνF^μν)^p)

where g_μν  denotes the metric tensor, g=detg_μν is the determinant of the metric, G_μν and R are the Einstein tensor and Ricci scalar correspondingly, Ł is the cosmological constant and  is the scalar field, nonminimally coupled to gravity and finally F_μν denotes the Maxwell field tensor and p is the factor (power) of nonlinearity and when p=1 the standard Maxwell action term is recovered.

Having used the principle of the least action we derive equations of motion which can be written in the form:

    G_μν+Λ g_μν=1/2(α T^(1)_μν+η T^(2)_μν)+T^(3)_μν

and here we use the following notations

    T^(1)_μν=_μ_ν-1/2g_μν^λ_λ,


    T^(2)_μν=1/2_μ_ν R-2^λ_ν R_λμ+1/2^λ_λ G_μν-g_μν(-1/2_λ_κ^λ^κ.
    .+1/2(^2)^2-R_λκ^λ^κ)
    -_μ^λ_ν_λ+
    _μ_ν^2-R_λμκν^λ^κ


    T^(3)_μν=g_μν/2(-F_łF^ł)^p+2p(-F_łF^ł)^p-1F_μρF_ν^ρ

It should be noted that the term T^(1)_μν is the standard energy-momentum tensor for a scalar field minimally coupled to gravity whereas the collection of terms denoted by T^(2)_μν is caused by the nonminimal coupling to gravity. We also remark that the terms denoted by T^(3)_μν comprise the energy-momentum tensor for nonlinear Maxwell field and it is easy to check that when p=1 one arrives at the standard Maxwell field energy-momentum tensor. Variation of the action (<ref>) with respect to the scalar field  gives rise to a corresponding equation of motion of the form:

    (α g_μν-η G_μν)^μ^ν=0.

And for the Maxwell field one obtains:

    _μ((-F_łF^ł)^p-1F^μν))=0.

One can verify that when p=1 the standard Maxwell equations can be recovered.

We are to obtain so called topological solutions and we assume that the metric we find takes the form:

    ds^2=-U(r)dt^2+W(r)dr^2+r^2dΩ^2()_(n-1),

where dΩ^2()_n-1 is the line element of a n-1-dimensional hypersurface of a constant curvature and it takes the form. 

    dΩ^2()_(n-1)=
    
    dþ^2+sin^2þdΩ^2_(n-2),   =1,
    
    dþ^2+þ^2 dΩ^2_(n-2),  =0,
    
    dþ^2+sinh^2þdΩ^2_(n-2),  =-1,

and here dΩ^2_(n-2) is the line element of a n-2–dimensional sphere. The given expression for dΩ^2()_(n-1) represents surfaces of positive, null or negative curvature for given values of parameter . 

Having supposed that the electromagnetic field has a scalar component only (A=A_0(r)dt) and taking into consideration the chosen form of the metric (<ref>) one can solve the Maxwell equations (<ref>) and obtain electromagnetic field tensor in the form:

    F_rt=-q/r^(n-1)/(2p-1)√(UW)

and here q is an integration constant related to the charge of the black hole.
Now we consider the equation for the scalar field (<ref>) and having integrated it for a once we derive:

    √(U/W)r^n-1[-(n-1)/2rW(U/U-(n-2)/r( W-1))]'=C

We assume that the constant C is equal to zero and this condition is equivalent to the following relation:

    g_rr- G_rr=0

Having used the metric (<ref>) we can write the equations (<ref>) in the following form:

    (n-1)/2rW(W/W+(n-2)/r( W-1))(1+3/4()^2/W)-Ł=/4W()^2+
    /2((n-1)/rW^2”+(n-1)(n-2)/r^2W^2()^2( W-1/2))+2p-1/2(2q^2/r^2(n-1)/(2p-1))^p


    (n-1)/2rW(U/U-(n-2)/r( W-1))(1+3/4()^2/W)+Ł=
    /4W()^2-/2(n-1)(n-2)/2r^2W()^2-2p-1/2(2q^2/r^2(n-1)/(2p-1))^p



    [1/2UW(U-(U)^2/2U-UW/2W)+n-2/2rW(U/U-W/W)-(n-2)(n-3)/2r^2W( W-1)]×
    (1+/4()^2/W)+Ł=-/4W()^2-/2W^2(U/2U+n-1/r)+
    /2()^2/W(UW/4UW^2+(n-2)W/2rW^2-(n-2)(n-3)/2r^2)+1/2(2q^2/r^2(n-1)/(2p-1))^p

Having combined  the equations (<ref>) and (<ref>) we arrive at the relations:

    ()^2=-4r^2W/(2 r^2+(n-1)(n-2))(+Ł+2^p-1(2p-1) q^2p/r^2p(n-1)/(2p-1))


    UW=((-Ł)r^2+(n-1)(n-2)-2^p-1(2p-1)q^2pr^2(1-p(n-1)/(2p-1)))^2/(2 r^2+(n-1)(n-2))^2.

It is worth being emphasized that the square of derivative of the scalar field should be positive definite in the domain outside black hole's horizon (since we are to obtain black hole's solution). To provide the positive definiteness of (')^2 a condition on the parameters of coupling ,  and cosmological constant Ł should be imposed, for example for positive  and  one should take the cosmological constant Ł to be a negative one since the metric function W(r) is supposed to be positive definite in the outer domain (if the parameter  is negative it can be shown that the negative value of the cosmological constant also should be taken negative).  In the inner domain (inside the black hole) the right hand side of the relation (<ref>) might change its sign and it means that the scalar field might turn to be a phantom-like field, but the behaviour of the field  inside the black hole does not affect on its behaviour in the outer domain. We would also like to point out that the right hand side of the relation (<ref>) is always nonnegative and it means that the sign of the metric functions U(r) and W(r) is always the same (when one of them changes its sign, so does the other function), as it has to be for a black hole. It should be noted that when p=1 the given above relations (<ref>) and (<ref>) reduce to the corresponding relations, obtained for linear Maxwell field <cit.> .
 The metric function U(r) can be represented in the form:
 
    U(r)=-μ/r^n-2-2Ł/n(n-1)r^2-2^p(2p-1)^2q^2p/(n-1)(2p-n)r^2(1-p(n-1)/2p-1)+
    (+Ł)^2/2(n-1)r^n-2∫r^n+1/r^2+d^2dr+2^p-1(2p-1)(+Ł)q^2p/(n-1)r^n-2×
    ∫r^n+1-2p(n-1)/2p-1/r^2+d^2dr+2^2p-3(2p-1)^2 q^4p/(n-1)r^n-2∫r^n+1-4p(n-1)/2p-1/r^2+d^2dr

 and here d^2=(n-1)(n-2)/2.
 The integrals in the written above formula gives rise to different expressions in case of the exponents of the r in the numerators of the integrals are integer or not. The form of the final expression for the metric function U(r) also depends on the sign of the introduced parameter d^2 and the relation between integration variable r and this parameter. The first of the mentioned integrals always has positive integer exponent, so when d^2>0 it can be represented in the form:

    ∫r^n+1/r^2+d^2dr=(-1)^(n+1)/2d^narctan(r/d)+∑^(n-1)/2_j=0(-1)^jd^2jr^n-2j)/n-2j

for odd values of n and

    ∫r^n+1/r^2+d^2dr=(-1)^n/2d^n/2ln(r^2/d^2+1)+∑^n/2-1_j=0(-1)^jd^2jr^n-2j)/n-2j

for the case of even n. One can conclude that the first integral in the formula (<ref>) can be expressed in terms of elementary functions consisting of power terms and a transcendental function which depends on the fact whether the dimension of space n is odd or even. 

The second and the third integrals in the right hand side of the relation (<ref>) might have nonpositive integer or noninteger exponents in the numerators. It should be pointed out here that the second integral might have positive exponent when p>(n+1)/4, but the exponent of the third integral always takes negative values.  In the case when 2p(n-1)/2p-1 is integer both exponents are integer and the integrals can be represented in the form:

    ∫r^[N_i]/r^2+d^2dr=(-1)^-[N_i]/2/d^1-[N_i]arctan(r/d)+∑^-[N_i]/2-1_j=0(-1)^jr^[N_i]+2j+1/([N_i]+2j+1)d^2(j+1)
 
when [N_i] is even and

    ∫r^[N_i]/r^2+d^2dr=(-1)^-([N_i]-1)/2/2d^1-[N_i]ln(1+d^2/r^2)+∑^-([N_i]+3)/2_j=0(-1)^jr^[N_i]+2j+1/([N_i]+2j+1)d^2(j+1)
 
for odd [N_i], and here we denoted N_1=n+1-2p(n-1)/2p-1 and N_2=n+1-4p(n-1)/2p-1 and the brackets [N_i] mean that the parameter N_i (i=1,2) takes an integer value. In general nor N_1 neither N_2 are integers, so the corresponding integrals can be expressed only through the special functions. In particular, for large r (large r means that r>b) we can write:

    ∫r^N_i/r^2+d^2dr=r^N_i-1∑^+∞_j=01/N_i-1-2j(-d^2/r^2)^j=r^N_i-1/(N_i-1)_2F_1(1,1-N_i/2;3-N_i/2;-d^2/r^2)

The series written above is convergent, when d^2/r^2<1. In contrast for small distances, namely, when r<d  another representation of the hypergeometric function can be used, so the latter integral can be rewritten in the form:

    ∫r^N_i/r^2+d^2dr=r^N_i+1/(1+N_i)d^2_2F_1(1,1+N_i/2;3+N_i/2; -r^2/d^2)

It is clear that when N_i is integer the given hypergeometric functions can be rewritten through the elementary logarithmic or inverse tangent functions.
Using the derived forms for integrals we can write the general relation for the metric function. Namely, for odd n and noninteger N_1 and N_2 we obtain:
 
    U(r)=-μ/r^n-2-2Ł/n(n-1)r^2+(+Ł)^2/2(n-1)[(-1)^(n+1)/2d^n/r^n-2arctan(r/d)+∑^(n-1)/2_j=0(-1)^jd^2jr^2(1-j)/n-2j]
    -2^p(2p-1)^2q^2p/(2p-n)(n-1)r^2(3p-pn-1)/2p-1+2^p-1(2p-1)^2(+Ł)q^2p/(n-1)(2p-n)r^2(1-p(n-1)/2p-1)×
    _2F_1(1,p(n-1)/2p-1-n/2;p(n-1)/2p-1-n/2+1;-d^2/r^2)+2^2p-3(2p-1)^3 q^4p/(n-1)(4p-2pn-n)r^2(1-2p(n-1)/2p-1)×
    _2F_1(1,2p(n-1)/2p-1-n/2;2p(n-1)/2p-1-n/2+1;-d^2/r^2).

For even n we arrive at the following expression:

    U(r)=-μ/r^n-2-2Ł/n(n-1)r^2+(+Ł)^2/2(n-1)[(-1)^n/2d^n/2r^n-2ln(r^2/d^2+1)+∑^n/2-1_j=0(-1)^jd^2jr^2(1-j)/n-2j]
    -2^p(2p-1)^2q^2p/(2p-n)(n-1)r^2(3p-pn-1)/2p-1+2^p-1(2p-1)^2(+Ł)q^2p/(n-1)(2p-n)r^2(1-p(n-1)/2p-1)×
    _2F_1(1,p(n-1)/2p-1-n/2;p(n-1)/2p-1-n/2+1;-d^2/r^2)+2^2p-3(2p-1)^3 q^4p/(n-1)(4p-2pn-n)r^2(1-2p(n-1)/2p-1)×
    _2F_1(1,2p(n-1)/2p-1-n/2;2p(n-1)/2p-1-n/2+1;-d^2/r^2).

 It should be noted that in the written above expressions for the metric function U(r) the relation (<ref>) is used and it allows to describe the behaviour of the solution for large r (when r>b), whereas for small r (r<b), the relation (<ref>) should be used. 

For linear electromagnetic field (p=1) the relations for the metric function U(r) obtained in the paper <cit.> are restored, namely for the odd n it can be represented in the form:

    U(r)=-μ/r^n-2-2Ł/n(n-1)r^2+(+Ł)^2/2(n-1)[(-1)^n+1/2d^n/r^n-2arctan(r/d)+∑^(n-1)/2_j=0(-1)^jd^2jr^2(1-j)/n-2j]
    +2q^2/(n-1)(n-2)r^2(2-n)+(+Ł)q^2/(n-1)r^n-2[(-1)^n-3/2/d^n-2arctan(r/d)+∑^(n-5)/2_j=0(-1)^jr^4-n+2j/(4-n+2j)d^2(j+1)]
    + q^4/2(n-1)r^n-2[(-1)^3n-5/2/d^3n-4arctan(r/d)+∑^(3n-7)/2_j=0(-1)^jr^6-3n+2j/(6-3n+2j)d^2(j+1)]

 and for even n we arrive at:

    U(r)=-μ/r^n-2-2Ł/n(n-1)r^2+(+Ł)^2/2(n-1)[(-1)^n/2d^n/2r^n-2ln(r^2/d^2+1)+∑^n/2-1_j=0(-1)^jd^2jr^2(1-j)/n-2j]
    +2q^2/(n-1)(n-2)r^2(2-n)+(+Ł)q^2/(n-1)r^n-2[(-1)^n-2/2/2d^n-2ln(1+d^2/r^2)+∑^n/2-3_j=0(-1)^jr^4-n+2j/(4-n+2j)d^2(j+1)]
    + q^4/2(n-1)r^n-2[(-1)^3n/2/2d^3n-4ln(1+d^2/r^2)+∑^3n/2-4_j=0(-1)^jr^6-3n+2j/(6-3n+2j)d^2(j+1)].

It should be pointed out here that for the linear case p=1 as well as for the cases when 2p(n-1)/2p-1 is integer then the relations (<ref>) and (<ref>) are used and the resulting form of the metric function U(r) contains only power-law and elementary transcendental functions.
 
When =0 (flat horizon surface)  the metric function can be written in the following way:

    U(r)=-μ/r^n-2+(-Ł)^2/2 n(n-1)r^2-2^p-1(-Ł)(2p-1)^2q^2p/(n-1)(2p-n)r^2(1-p(n-1)/2p-1)
    +2^2p-3(2p-1)^3 q^4p/(n-1)(4p-n-2pn)r^2(1-2p(n-1)/2p-1)


We investigate the behaviour of the metrics (<ref>) and (<ref>) for some asymptotic values of r and η. Firstly, we examine asymptotic  behaviour of metric function U(r) when r→ 0. To perform this task we decompose  all the functions comprising the mentioned function U(r) into a series for small distances. As a result, for small r→ 0 the function can be written in the form:

    U(r)=-μ/r^n-2-2Ł/n(n-1)r^2-2^p(2p-1)^2q^2p/(2p-n)(n-1)r^2(1-p(n-1)/2p-1)+(+Ł)^2/^2(n-1)^2(n^2-4)r^4
    +2^p(2p-1)^2(+Ł)q^2p/(n-1)^2(n-2)(6p-n-2)r^2(2-p(n-1)/2p-1)+2^2(p-1)(2p-1)^3/(n-1)^2(n-2)×
    q^4p/2p(4-n)-n-2r^4(1-p(n-1)/2p-1)+ O_1(r^6)+ O_2(r^2(3-p(n-1)/2p-1))+ O_3(r^2(3-2p(n-1)/2p-1))

and here O_i(r^γ_i) denotes the terms that appear as the following terms in the decomposition of corresponding functions in U(r). It should be noted that the first four terms in the expansion (<ref>) completely recover the standard solution for nonlinearly charged black hole with a cosmological constant, but the relation (<ref>) does not give rise to the condition U(r)W(r)≃ 1 which takes place in the standard Einsteinian theory.

Having used asymptotic expansions of the functions in the relations (<ref>) and (<ref>) for when r→∞ we can write the expressions for the metric function U(r) when the distance r tends to infinity, namely for odd n we arrive at:

    U(r)=(1-(+Ł)^2/4^2)+(-Ł)^2/2n(n-1)r^2-μ/r^n-2+(+Ł)^2/2(n-1)((-1)^n+1/2π d^n/2 r^n-2+.
    .∑^(n-1)/2_j=2(-1)^jd^2jr^2(1-j)/n-2j)-2^p-1(2p-1)^2(-Ł)q^2p/(n-1)(2p-n)r^2(1-p(n-1)/2p-1)+
    2^2p-3(2p-1)^3 q^4p/(n-1)(4p-n-2pn)r^2(1-2p(n-1)/2p-1)+ O_1(1/r^n-1)+ O_2(1/r^2p(n-1)/2p-1)+ O_3(1/r^4p(n-1)/2p-1),

and for even n we obtain:

    U(r)=(1-(+Ł)^2/4^2)+(-Ł)^2/2n(n-1)r^2-μ/r^n-2+(+Ł)^2/2(n-1)((-1)^n/2d^n/r^n-2ln(r/d)+.
    .∑^n/2-1_j=2(-1)^jd^2jr^2(1-j)/n-2j)-2^p-1(2p-1)^2(-Ł)q^2p/(n-1)(2p-n)r^2(1-p(n-1)/2p-1)+
    2^2p-3(2p-1)^3 q^4p/(n-1)(4p-n-2pn)r^2(1-2p(n-1)/2p-1)+ O_1(1/r^n)+ O_2(1/r^2p(n-1)/2p-1)+ O_3(1/r^4p(n-1)/2p-1),

and here O_i denotes the subleading terms in the expansions. It should be noted that the behaviour of the leading term of the metric function U(r) at infinity does not depend on choice of the parameter  and is completely the same for three types of geometry and this fact has complete analogy with topological black holes in the framework of the standard General Relativity. So, we can conclude that the behaviour of the metric functions at the infinity is deSitterian or anti-deSitterian which depends on the sign of .  The other moment should be paid attention to is the fact that the leading terms related to the charge (those which contain the parameter q) is the same for all three types of geometry. It should also be stressed that the asymptotic behaviour of the metric function when r→ 0 and r→∞ is very important for the calculation of Kretschmann scalar which gives information about true and coordinate singularities of the metric.

We also examine the regime of large nonminimal derivative  coupling (large ), this situation means that the term with minimal coupling (the term with factor ) is considerably smaller than the term which represents nonminimal coupling (the term with factor ). For this situation the metric function U(r) can be represented as follows:

    U(r)=-μ/r^n-2-2Ł/n(n-1)r^2-2^p(2p-1)^2q^2p/(2p-n)(n-1)r^2(1-p(n-1)/2p-1)+Ł^2/(n-1)^2(n^2-4)r^4+
    2^p(2p-1)^2Ł q^2p/(n-1)^2(n-2)(6p-n-2)r^2(2-p(n-1)/2p-1)+2^2(p-1)(2p-1)^3 q^4p/(n-1)^2(n-2)(8p-2pn-n-2)r^4(1-p(n-1)/2p-1)+ O(1/)

The latter relation can be compared with the relation (<ref>) obtained for the case of flat horizon geometry and as it is easy to verify that the flat geometry metric contains the terms with the dependence ∼ whereas for the case of nonflat geometry (=± 1) such kind of terms does not appear.
The product of the metric functions (<ref>) for the regime of large  takes the form:

    UW=(1-Ł r^2/(n-1)(n-2)-2^p-1(2p-1)q^2pr^2(1-p(n-1)/2p-1)/(n-1)(n-2))^2+ O(1/).

The written above relation is again valid when =± 1 and for the flat-type geometry (=0) the product of metric functions in the regime of large  is as follows:

    UW=(Ł^2+2^p(2p-1)Ł q^2pr^-2p(n-1)/2p-1+2^2(p-1)(2p-1)^2q^4pr^-4p(n-1)/2p-1)^2/4^2+ O()

which goes in agreement with the fact that the function U(r) behaves as ∼ in this regime.


The figure Fig.[<ref>] shows the behaviour of the metric function U(r) for linear (p=1) and conformal cases (p=n+1/4). In general the behaviour of these functions is very similar. Both of them go up to infinity with increasing of r and the difference between them becomes negligibly small because of the fact that at the large distances the dominant terms in both functions are  of AdS-type which have the same fixed parameters. It should be pointed out that this conclusion is valid for any allowed value of p. For small distances (r→ 0) the metric function for linear case goes faster to minus infinity than in the conformal case. The quite general conclusion we can make here is the fact that with increasing of the parameter of nonlinearity p the behaviour of the metric function for r→ 0 becomes weaker than for smaller value of p. 

In order to find out where the metric has physical singularities we should analyze the behaviour of Kretschmann scalar. For our types of metrics it can be written in the form:

    R_μνł̨R^μνł̨=1/UW(d/dr[U'/√(UW)])^2+(n-1)/r^2W^2((U')^2/U^2+(W')^2/W^2)+2(n-1)(n-2)/r^4W^2( W-1)^2.

The represented form of the Kretschmann scalar is suitable for all three types of geometry and certainly the evident form of this scalar would depend on it but for simplicity we have not substituted the final forms of the functions U and W into the written expression of the scalar (<ref>).

Having written the Kretschmann scalar (<ref>) one can analyze its behaviour for different values of r. First of all it is necessary to investigate its behaviour at the horizon points, namely at the point r_+ where  U(r_+)=0 and the function W(r) is singular. Taking into account the relation (<ref>) it is easy to persuade oneself that the horizon points are the points of coordinate singularity, where the Kretschmann scalar has finite values that is usually takes place for black hole's solution. When r→ 0 we arrive at:

    R_μνł̨R^μνł̨∼16(n-2)^2(3p-pn-1)^2/(8p-2pn-n-2)^2((4p-2pn-1)^2/(2p-1)^2+n(n-1)/n-2)1/r^4

and as a consequence a physical singularity appears at the origin. It is worth being noted that for all of dimensions n the character of singularity is identical what is in contrast for example with the case of chargeless black hole, where the character of this singularity depends on the dimension of space <cit.> or with the case of a charged black hole in the standard General Relativity where this dependence also takes place. This peculiarity can be explained by the fact that in the case we consider now the product of the metric functions UW is also singular what does not happen nor in the chargeless case <cit.> neither in the standard GR, where we have UW=1. When r→∞ we find:

    R_μνł̨R^μνł̨∼8(n+1)^2/n(n-1)^2^2.

It should be noted that the Kretschmann scalar at the infinity behaves identically to corresponding value for the chargeless black hole <cit.> and this fact is obvious because both types of black holes have asymptotically AdS-like behaviour at the infinity. To sum up the investigation of Kretschmann scalar we can conclude that its behaviour is rather typical for a black hole with the only singularity at the origin, coordinate singularities at the horizons and some asymptotic value at the infinity.

Now we can proceed to the investigation of the solutions of equations for  gravitational field in the case of negative d^2 (d^2<0). This situation takes place when the parameters  and  have opposite signs.  The general formula (<ref>) is also valid in this case, but due to the negativity of d^2 the results of integration would be different. Not to introduce new notations here we denote d^2=-|d|^2 and |·| means the absolute value of the parameter. The results of integration again depend on the parity of n and the fact whether N_1 and N_2 are integer or not. For the first integral in the relation (<ref>) we write:

    ∫r^n+1/r^2-|d|^2dr=|d|^n/2ln|r-|d|/r+|d||+∑^(n-1)/2_j=0|d|^2jr^n-2j/n-2j

for odd values of n and

    ∫r^n+1/r^2-|d|^2dr=|d|^n/2ln|r^2/|d|^2-1|+∑^n/2-1_j=0|d|^2jr^n-2j/n-2j

for even n. For noninteger values of N_1 and N_2 we can utilize the relation (<ref>) making the replacement -d^2→|d|^2 and for integer values of N_i it can calculated directly and written in terms of elementary functions or having the representation through the special functions (<ref>) and using some limit procedures, when parameters of the hypergeometric functions become integer.

Now we can write the expression for the metric function U(r), namely for odd n and r>|d| we arrive at the following expression:
 
    U(r)=-μ/r^n-2-2Ł/n(n-1)r^2+(+Ł)^2/2(n-1)[|d|^n/2r^n-2ln|r-|d|/r+|d||+∑^(n-1)/2_j=0|d|^2jr^2(1-j)/n-2j]
    -2^p(2p-1)^2q^2p/(2p-n)(n-1)r^2(3p-pn-1)/2p-1+2^p-1(2p-1)^2(+Ł)q^2p/(n-1)(2p-n)r^2(1-p(n-1)/2p-1)×
    _2F_1(1,p(n-1)/2p-1-n/2;p(n-1)/2p-1-n/2+1;|d|^2/r^2)+2^2p-3(2p-1)^3 q^4p/(n-1)(4p-2pn-n)r^2(1-2p(n-1)/2p-1)×
    _2F_1(1,2p(n-1)/2p-1-n/2;2p(n-1)/2p-1-n/2+1;|d|^2/r^2),

 whereas for the even n and r>|d| we obtain:
 
    U(r)=-μ/r^n-2-2Ł/n(n-1)r^2+(+Ł)^2/2(n-1)[|d|^n/2r^n-2ln|r^2/|d|^2-1|+∑^n/2-1_j=0|d|^2jr^2(1-j)/n-2j]
    -2^p(2p-1)^2q^2p/(2p-n)(n-1)r^2(3p-pn-1)/2p-1+2^p-1(2p-1)^2(+Ł)q^2p/(n-1)(2p-n)r^2(1-p(n-1)/2p-1)×
    _2F_1(1,p(n-1)/2p-1-n/2;p(n-1)/2p-1-n/2+1;|d|^2/r^2)+2^2p-3(2p-1)^3 q^4p/(n-1)(4p-2pn-n)r^2(1-2p(n-1)/2p-1)×
    _2F_1(1,2p(n-1)/2p-1-n/2;2p(n-1)/2p-1-n/2+1;|d|^2/r^2).

 We remark that similarly to the previously considered case of positive d^2  another representation of hypergeometric function should be utilized when r<|d|. It should be noted that the given relations for the metric function have coordinate singularity at the point r=d which is also the point  where the function (')^2 becomes singular. It is easy to check that for the distances larger than |d| the function U(r) for both parities of n ((<ref>) and (<ref>)) has a cosmological horizon in complete analogy with the chargeless case <cit.>, thus it is difficult to treat the obtained solution as a black hole. 
 
 In the following we will consider thermodynamics of black holes represented by the metric functions (<ref>) and (<ref>). Since the black hole's mass is extremely important for the thermodynamic we use the relations (<ref>) and (<ref>) to express the mass parameter μ as a function of horizon radius r_+. As a result we obtain:

    μ=(1-(+Ł)^2/4^2)r^n-2_++(-Ł)^2/2n(n-1)r^n_++(+Ł)^2/2(n-1)[(-1)^n+1/2d^narctan(r_+/d)+.
    .∑^(n-1)/2_j=2(-1)^jd^2jr^n-2j_+/n-2j]-2^p(2p-1)^2q^2p/(2p-n)(n-1)r^2p-n/2p-1_+(1-(+Ł)/2×.
    ._2F_1(1,p(n-1)/2p-1-n/2;p(n-1)/2p-1-n/2+1;-d^2/r^2_+))+2^2p-3(2p-1)^3 q^4p/(n-1)(4p-2pn-n)×
    r^4p-2pn-n/2p-1_+_2F_1(1,2p(n-1)/2p-1-n/2;2p(n-1)/2p-1-n/2+1;-d^2/r^2_+)

for odd n and for even n we arrive at:

    μ=(1-(+Ł)^2/4^2)r^n-2_++(-Ł)^2/2n(n-1)r^n_++(+Ł)^2/2(n-1)[(-1)^n/2d^n/2ln(r^2_+/d^2+1)+.
    .∑^n/2-1_j=2(-1)^jd^2jr^n-2j_+/n-2j]-2^p(2p-1)^2q^2p/(2p-n)(n-1)r^2p-n/2p-1(1-(+Ł)/2×.
    ._2F_1(1,p(n-1)/2p-1-n/2;p(n-1)/2p-1-n/2+1;-d^2/r^2_+))+2^2p-3(2p-1)^3 q^4p/(n-1)(4p-2pn-n)×
    r^4p-2pn-n/2p-1_+_2F_1(1,2p(n-1)/2p-1-n/2;2p(n-1)/2p-1-n/2+1;-d^2/r^2_+)


The dependence μ=μ(r_+) is illustrated on the Fig.[<ref>]. One can see that the crucial difference between them is in range of small r_+, whereas for large r_+ the difference between them becomes negligibly small and this fact can be explained by the dominance of the ∼ r^n_+ term, which has the same factor for all values of the parameter p.
 


 §.§ Field potential and black hole's charge

Since we have obtained a charged black hole and we are going to obtain the first law of black hole thermodynamics in the following it is very important to have the relation for the total charge of the black hole and the evident form of the vector potential. As it was noted above, the potential has the only component A_0(r) which can be recovered by integration of the relation (<ref>), namely we can write:

    A_t≡ψ=ψ_0-∫q/r^n-1/2p-1√(UW)dr

where ψ_0 denotes a constant which is taken to impose some condition on the value of potential, namely in our case we suppose that ψ_0 is nonzero, taken in a such way that at the event horizon the the potential ψ is equal to zero. It should be noted that we can choose the constant in other way, namely we can impose that at the infinity ψ→ 0, but our choice is a bit convenient when we consider black hole thermodynamics. Having substituted the relation (<ref>) into (<ref>) and performing integration we can write:

    ψ(r)=ψ_0+2p-1/n-2pqr^2p-n/2p-1+(2p-1)(+Ł)q/2(2p-n)r^2p-n/2p-1_2F_1(1,n-1/2(2p-1)-1/2;n-1/2(2p-1)+1/2;-d^2/r^2)
    
    +2^p-2(2p-1) q^2p+1/(4p-2pn-n)r^1-(2p+1)(n-1)/2p-1_2F_1(1,(2p+1)(n-1)/2(2p-1)-1/2;(2p+1)(n-1)/2(2p-1)+1/2;-d^2/r^2)

and here it is assumed that 2p(n-1)/2p-1 is noninteger. It should also be remarked that the written hypergeometric functions are valid for large distances (r>d) and they tend to one when r→∞, so at large distances the potential ψ demonstrate power-law dependence. Similarly to the metric function U(r)  for small distances (r<d) another representation for the hypergeometric functions should be used, but again in the limit r→ 0 the potential shows power-law dependence. To obtain finite value of the potential at the infinity one should impose that 2p-n/2p-1<0 which lead to restriction on the value of parameter of nonlinearity p, namely we obtain:

    1/2<p<n/2.

It should be noted that written above condition for the parameter p is only related to electromagnetic part of the action and there are no imprints of gravity part. We also remark that some conditions on the parameter of nonlinearity was obtained in the paper <cit.> where a black hole in Einstein-dilaton-nonlinear Maxwell theory was investigated. The quite general conclusion one can make here is that the parameter of nonlinearity cannot be chosen completely arbitrary, in any kind of theory one has some conditions which are imposed on the parameter of nonlinearity and those conditions are taken to obtain reasonable  behaviour of electric potential and possibly some other functions.

In case of linear field p=1 (it satisfies the condition (<ref>)) the gauge potential (<ref>) can be written in terms of elementary functions. We write them in the separate form for odd and even n respectively. Namely, for odd n one arrives at:

    ψ(r)=ψ_0+q/n-2r^2-n+(+Ł)q/2[∑^(n-5)/2_j=0(-1)^jr^2j-n+4/d^2(j+1)(2j-n+4)+(-1)^n-3/2/d^n-2arctan(r/d)]
    + q^3/2[∑^(3n-7)/2_j=0(-1)^jr^2j-3n+6/d^2(j+1)(2j-3n+6)+(-1)^3n-5/2/d^3n-4arctan(r/d)],
 
and for even n the potential is as follows:

    ψ(r)=ψ_0+q/n-2r^2-n+(+Ł)q/2[∑^(n-6)/2_j=0(-1)^jr^2j-n+4/d^2(j+1)(2j-n+4)+(-1)^n/2+1/2d^n-2ln(1+d^2/r^2)]
    + q^3/2[∑^(3n-8)/2_j=0(-1)^jr^2j-3n+6/d^2(j+1)(2j-3n+6)+(-1)^3n/2/2d^3n-4ln(1+d^2/r^2)].
 
It should be noted that for large r the potential for both parities of n (relations (<ref>) and (<ref>)) has reasonable behavior, because it contains the terms of inverse proportionality (∼1/r^m, m>0), and it means that at the infinity the potential is finite.

We also note that for the case when p=n+1/4 (conformal case) the gauge potential can also be written in terms of elementary functions. So for this type of nonlinearity we can write for odd n:

    ψ(r)=ψ_0+q/r+(+Ł)q/2 darctan(r/d)+2^n-3/4(n-1) q^n+3/2/4×
    [∑^(n+1)/2_j=0(-1)^jr^2j-2-n/d^2(j+1)(2j-2-n)+(-1)^n+3/2/d^n+4arctan(r/d)];

and for even n we arrive at:

    ψ(r)=ψ_0+q/r+(+Ł)q/2 darctan(r/d)+2^n-3/4(n-1) q^n+3/2/4×
    [∑^n/2_j=0(-1)^jr^2j-2-n/d^2(j+1)(2j-2-n)+(-1)^n/2/2d^n+4ln(1+d^2/r^2)];

Having compared the written above relations for potential in case of conformal nonlinearity with the case of linear field one can conclude that    the conformal case has the terms that go slower to zero at the infinity than in linear case, so in the considered above linear case the potential goes faster to its limit values than in conformal case.

Apart of the gauge potential we have just considered, the other important quantity is the charge of the black hole. To obtain it we utilize Gauss law which in nonlinear case can be represented in the form:

    Q=1/4π∫_Σ(-F_μνF^μν)^p-1*F,

where *F denotes a Hodge dual of electromagnetic field form and integral is taken over n-1-dimensional closed spacelike hypersurface Σ. Having calculated the latter integral we arrive at:

    Q=ω_n-1/4π2^p-1q^2p-1

and here ω_n-1 denotes the volume (or surface area) of “unit” hypersurface of constant curvature (it reduces to the surface area of unit sphere in case of spherical symmetry). Now we can define total electric charge of the black hole per unit volume ω_n-1 as follows:

    Q̅=1/ω_n-1Q=1/4π2^p-1q^2p-1

The electric potential that an observer can measure at infinity with respect to the horizon can be defined in the form:

    Φ_q=A_μχ^μ|_+∞-A_μχ^μ|_r_+

where χ^μ is a Killing vector, null on the horizon. If one considers linear gauge theory (p=1) the Killing vector might be the time translation vector χ^μ=∂/∂ t, but for p≠ 1 to obtain consistent theory one should make other choice for the Killing vector, namely it was shown that for power-law field theory one should take : χ^μ=p∂/∂ t <cit.> and as a result the electric potential of the black hole takes the form:

    Φ_q=pψ_0.

One can see that the electric potential the observer can measure is completely defined by the integration constant ψ_0, which provides the  condition that the gauge potential A_μ  is equal to zero at the horizon and the parameter of nonlinearity p.  It should be pointed out that one can impose that A_μ=0 at the infinity but the electric potential Φ_q defined by the relation (<ref>) would be completely identical to the expression (<ref>).



§ THERMODYNAMICS OF THE BLACK HOLE



 §.§ Black hole temperature

To obtain black holes thermodynamics we start from calculation of its temperature. The temperature is supposed to be well defined quantity in quite general setting due to its strict geometrical basis <cit.>. To calculate the temperature we use the so called surface gravity which is defined as follows:

    κ^2=-1/2∇_aχ_b∇^aχ^b,

where χ^μ is a Killing vector field which is null on the horizon. Because we consider the static metric we can chose χ^ν=∂/∂ t. As a result we can write relation for black holes temperature:

    T=/2π=1/4πU(r_+)/√(U(r_+)W(r_+))
  
and here r_+ similarly as in the previous section denotes the   event horizons radius. It is well understood that the evident form for black hole's temperature will depend on the parity of n and the fact whether the mentioned before numbers N_1 and N_2 are integers or not. Under the assumption that the parameters N_1 and N_2 are noninteger we can write the relation for temperature as follows:

    T=1/4π2 r^2_++(n-1)(n-2)/(-Ł)r^2_++(n-1)(n-2)-2^p-1(2p-1) q^2pr^2(1-p(n-1)/2p-1)_+((-Ł)^2/2(n-1)r_++.
    .(1-(+Ł)^2/4^2)(n-2)/r_++(+Ł)^2/2(n-1)[(-1)^(n+1-σ)/2d^n+1-σ/r^n-2-σ_+(r^2_++d^2)+∑^(n-1-σ)/2_j=2(-1)^jd^2jr^1-2j_+].
    .-2^p(2p-1)q^2p/n-1r^1-2p(n-1)/2p-1_+[1-(+Ł)/2(_2F_1(1,p(n-1)/2p-1-n/2;p(n-1)/2p-1-n/2+1;-d^2/r^2_+)-...
    ...2(2p-1)d^2/(n+2p-2)r^2_+_2F_1(2,p(n-1)/2p-1-n/2+1;p(n-1)/2p-1-n/2+2;-d^2/r^2_+))]+2^2p-3(2p-1)^2 q^4p/(n-1)×.
    .r^1-4p(n-1)/2p-1_+[_2F_1(1,2p(n-1)/2p-1-n/2;2p(n-1)/2p-1-n/2+1;-d^2/r^2_+)-2(2p-1)d^2/(2pn+n-2)r^2_+×..
    .._2F_1(2,2p(n-1)/2p-1-n/2+1;2p(n-1)/2p-1-n/2+2;-d^2/r^2_+)])
 
and here σ=0 when n is odd and σ=1 for even n. It is worth noting that latter relation represents the temperature as a function of the horizon radius r_+, for r_+>d, if r_+<d the representation of integrals given by the relation (<ref>) should be utilized in the metric function (<ref>) instead of the relation (<ref>) which has brought  us to the written above relation (<ref>). The behaviour of the function T=T(r_+) is rather complicated, but nevertheless looking at the relation (<ref>) some conclusions about the behaviour of this function can be made. The most important point here that should be paid attention to is the fact that for large r_+ the temperature grows up almost linearly due to corresponding linear term and the fact that all the other terms are inversely proportional to some power of r_+. The hypergeometric functions in the relation (<ref>) go to one when r_+→∞ and due to the presence of factors r^1-2p(n-1)/(2p-1)_+, r^1-4p(n-1)/(2p-1)_+ and r^-2_+ the contribution of the terms related to electromagnetic field in the temperature (<ref>) becomes negligibly small when r_+→∞. As it has just been noted, to obtain general relation for the temperature when r_+<d one should use the relation (<ref>) and the resulting hypergeometric functions one arrives at will go to one when r_+→ 0, but the terms with hypergeometric functions contain the factors of inverse proportionality ∼ 1/r^γ_+  (γ>0), thus the black hole temperature would have power-law singularity when r_+→ 0. So having relied on the mentioned above facts we can conclude that the temperature might have nonmonotonous behaviour which gives rise to the consequence that in this case we have a phase transition. It is also worth being remarked that a Hawking-Page phase transition usually takes place in AdS-like cases, due to this fact  the phase transition for our black hole supposedly would be of that type. 

For the particular case of linear field (p=1) we can write the temperature in the following form:

    T=1/4π2 r^2_++(n-1)(n-2)/(-Ł)r^2_++(n-1)(n-2)- q^2r^2(2-n)_+((-Ł)^2/2(n-1)r_++(1-(+Ł)^2/4^2)×.
    .(n-2)/r_++(+Ł)^2/2(n-1)[(-1)^(n+1-σ)/2d^n+1-σ/r^n-2-σ_+(r^2_++d^2)+∑^n-1-σ/2_j=2(-1)^jd^2jr^1-2j_+]-2q^2/n-1r^3-2n_++(+Ł)q^2/(n-1)×.
    .[(-1)^(n-3-σ)/2d^3+σ-n/r^n-2+σ_+(r^2_++d^2)+∑^n-5-σ/2_j=0(-1)^jr^2j+5-2n_+/d^2(j+1)]+ q^4/2(n-1)[(-1)^(3n-5-σ)/2d^5+σ-3n/r^n-2+σ_+(r^2_++d^2)+∑^3n-7-σ/2_j=0(-1)^jr^2j+7-4n_+/d^2(j+1)])

The written expression for the temperature in linear case (p=1) inherits the main features that we described for the nonlinear case (p≠ 1). In particular, it has almost linear increase for large radii, distorted by a weak nonlinearities which become negligibly small for large r_+. If the raduis of horizon r_+ decreases and for some value of r_+ it might turn negative, but this negativity would mean that the black holes with the other fixed parameters does not exist. When r_+→ 0 the temperature remains negative and goes down to infinitely low values and this decrease is mainly driven by the term ∼ -1/r^2n-3.

For the other interesting particular case, namely when p=(n+1)/4 we have the following expression for temperature:

    T=1/4π2 r^2_++(n-1)(n-2)/(-Ł)r^2_++(n-1)(n-2)-2^n-7/4(n-1) q^n+1/2r^1-n_+((-Ł)^2/2(n-1)r_++.
    .(1-(+Ł)^2/4^2)(n-2)/r_+-2^n-3/4q^n+1/2/r^n_++2^n-7/4(+Ł)q^n+1/2/(r^2_++d^2)r^n-2_++(+Ł)^2/2(n-1)[∑^n-1-σ/2_j=2(-1)^jd^2jr^1-2j_+..
    ..+(-1)^n+1-σ/2d^n+1-σ/r^n-2-σ_+(r^2_++d^2)]+2^n-9/2(n-1) q^n+1/[∑^n-1-σ/2_j=0(-1)^jr^2(j-n)+1_+/d^2(j+1)+(-1)^n+1-σ/2/d^n+1-σr^n-2+σ_+(r^2_++d^2)]).
 
The given expression shows that similarly to the previously considered case for large r_+ the temperature (<ref>) grows up almost linearly with increasing of r_+, while for infinitely low values of the horizon radius r_+ (r_+→ 0) the behaviour of the temperature is defined by the leading term ∼ -1/r^n_+, and in comparison with the previously examined linear case in the conformal case (p=(n+1)/4) the temperature decreases slowly. It should be stressed here that the mentioned above ranges of horizon radii r_+ when the temperature of the black holes turns to be negative are unphysical and we refer them just to have more complete information about the function T=T(r_+).  


The Fig.[<ref>] shows the dependence of the temperature T on the horizon radius r_+ for linear as well as conformal cases. Those dependences have some common features, namely for some values of cosmological constant greater in absolute value than some fixed (critical) one the temperature goes up monotonically with increasing of the radius of horizon r_+. As it will be shown below that the thermodynamic pressure we are going to introduce is proportional to the absolute value of the cosmological constant <cit.> and as a consequence the increasing of the module of Ł means the increasing of the pressure. Under the assumption about thermodynamic pressure the Fig.[<ref>] has direct analogy with the behaviour of TV-diagram of a substance described by Van der Waals equation of state which is applicable to real systems in condensed matter physics such as liquids or gases <cit.>.



 §.§ Wald procedure and entropy

Consistent approach which allows one to obtain conserved charges of a black hole was developed by Wald <cit.>. Having used the derived conserved quantities one can develop black hole's thermodynamics. It should be noted that Wald formalism is applicable to general diffeomorphism-invariant theories even those containing the derivatives of higher order. Wald formalism was applied to numerous  examples of black holes, in particular, Einstein-scalar theory <cit.>, Einstein-Proca theory <cit.>, Einstein-Yang-Mills theory <cit.>, gravity with quadratic-curvature invariants <cit.>, Lifshitz black hole <cit.>, black holes in Horndeski theory <cit.>.

It was argued <cit.> that the Wald relation for entropy should be corrected due to singular behaviour of the scalar potential at the horizon. To obtain the first law of black hole and relation for entropy the general Wald procedure was used <cit.>. The resulting expression for entropy differs from the well-known relation in General Relativity. In our work we also use the Wald procedure and before its application we briefly review the keypoints of this formalism.  Having supposed that we have some diffeomorphism-invariant theory described by a Lagrangian function L, we perform the variation of the Lagrangian which can be written in the following form:

    δ L= e.o.m.+√(-g)∇_νJ^ν

and here e.o.m. represents the terms which give the equations of motion and J^ν is the surface “current” and the corresponding “current” 1-form can be constructed J_(1)=J_νdx^ν. This surface term allows to construct Hodge dual form Θ_(n)=*J_(1). Having chosen a vector of infinitesimal diffeomorphism δ x^μ=ξ^μ one constructs a form:

    J_(n)=Θ_(n)-i_ξ* L=e.o.m.-d*J_(2)

where i_ξ* L is the contraction of the vector field ξ^μ with the form * L (here L denotes the Lagrangian of the considered system). It should be noted that when the equations  of motion are satisfied (on-shell condition) one can consider n-1-form  Q_(n-1)=*J_(2) for which J_(n)=dQ_(n-1). To obtain some thermodynamic relations the vector field ξ^μ should be taken to be a Killing vector field that is null at the horizon. According to Wald approach the variation of the Hamiltonian can be represented in the form: 

    δ H=1/16π(δ∫_cJ_(n)-∫_cd(i_ξΘ_(n)))=1/16π∫_Σ_(n-1)(δ Q-i_ξΘ_(n))

and here c is a n–dimensional Cauchy surface and Σ_(n-1) denotes its n-1-dimensional boundary with two components, namely one on the infinity and the other is on the event horizon. The first law of black hole thermodynamics follows from the relation:

    δ H_∞=δ H_+

where indices of the variation of the Hamiltonian, namely ∞ and + mean that these variations are taken at the infinity and at the event horizon respectively. It is clear that the first law of thermodynamics can be constructed if the variation δ H is a well behaved function in the domain in between two mentioned above boundaries.

Now we can proceed to the direct calculations of the variation of the Hamiltonian. It should be noted that the Wald method in Horndeski theory was examined in <cit.> for a chargeless black hole and in charged case it was considered in <cit.>.  We can write:

    (δ Q-i_ξΘ)_min=r√(UW)((n-1)/rW^2δ W+2p(2p-1)2^p-1/(UW)^p(ψ')^2(p-1)×.
    .ψ(ψ'(δ W/W+δ U/U)-2δψ')-'/Wδ)Ω_(n-1),


    (δ Q-i_ξΘ)_nm=η(n-1)/2r^n-2√(U/W)((')^2/2W^2δ W-δ((')^2/W)+2 r/(n-1)'δ)Ω_(n-1),

where the relations (<ref>) and (<ref>) are related to the minimally and nonminimally coupled parts of the action respectively. Their difference takes the form:

    (δ Q-i_ξΘ)_tot=r^n-1√(UW)((n-1)/rW^2(1+/4(')^2/W)δ W+2p(2p-1)2^p-1/(UW)^p(ψ')^2(p-1)×.
    .ψ(ψ'(δ W/W+δ U/U)-2δψ')-(n-1)/2rWδ((')^2/W))Ω_(n-1).

It should be pointed out that if p=1 the latter relation gets transformed into the corresponding relation obtained for a linear field <cit.>. The important point we should pay attention to is the fact that for a chargeless black hole the latter term in the total variation (<ref>), namely δ((')^2/W) is equal to zero identically for all r, whereas for the charged linear case p=1 it equals to zero at infinity because of the decay of electromagnetic filed at infinity, at the same time the variation might not be equal to zero at  the horizon but it has the metric function 1/W as a factor and as a consequence the last term in the total variation also disappears on the horizon. In our case of nonlinear field we can write the this last term of the variation (<ref>) at infinity in the following form:

    r^n-2√(UW)/Wδ((')^2/W)|_∞=r^n-2U/√(UW)δ((')^2/W)|_∞⟶∼ r^n-2p(n-1)/2p-1.

As it is easy to see that for providing the disappearance of the given variation one has to impose some conditions on the parameter p, namely we get completely the same condition as the previously derived (<ref>). This fact is an additional evidence that the parameter p can not be taken arbitrary, but it should be chosen in the  interval (<ref>) to obtain reasonable behaviour of the electric potential at infinity and as we have just seen to derive the  first law of black hole thermodynamics. At the horizon this variation is equal to zero identically as it was for linear case. Finally, we can write total variation at the infinity in the form:

    (δ Q-i_ξΘ)_tot=(n-1)δμ-4p(2p-1)2^p-1q^2(p-1)ψ_0δ q.

Taking into account the relation for the total charge (<ref>) and potential Φ_q (<ref>) and performing integration over angular variables we can write the variation of total Hamiltonian at the infinity in the form:

    H_∞=δ M-Φ_qδ Q,

where δ M is the variation of black hole's mass which can be defined as follows:

    M=(n-1)ω_n-1/16πμ.

The obtained expression for the black hole's mass (<ref>) does not depend anyhow on the nonminimal coupling of scalar field  with gravity and has the same form form as for Schwarzchild-AdS (or RN-AdS) black hole. Variating the Hamiltonian at the horizon we can write:

    δ H_+=(n-1)ω_n-1/16πU(r_+)r^n-2_+δ r_+=√(U(r_+)W(r_+))Tδ( A/4)=(1+/4()^2/W|_r_+)Tδ( A/4).

where A=ω_n-1r^n-1_+ is the horizon area of the black hole. One can see that it is not possible to represent the obtained variation in the form Tδ S which is common in the standard General Relativity.  Now following the procedure proposed in <cit.> and developed in the work <cit.> we can rewrite the latter relation as follows:

    H_+=Tδ S+Φ^+_δ Q^+_,

where S is entropy of the black hole, Q^+_ is the so called scalar “charge” and corresponding to it potential Φ^+_  can be defined in the form:

    S=(1+/4()^2/W|_r_+) A/4,


    Q^+_=ω_n-1√(1+/4()^2/W|_r_+),   Φ^+_=- AT/2ω_n-1√(1+/4()^2/W|_r_+).

We note that the chosen form for the scalar “charge” Q^+_ and conjugate potential Ψ^+_ is not unique, in particular in the work <cit.> different expressions were utilized. It can be verified easily that the chosen above forms for the scalar potential and its conjugate value (<ref>) allows one to obtain the following important relation:  

    Φ^+_Q^+_=- AT/2(1+/4()^2/W|_r_+)=-2TS


Having equated the right hand sides of the relations (<ref>) and (<ref>) one arrives at:

    δ M=Tδ S+Φ^+_δ Q^+_+Φ_qδ Q



 §.§ Extended thermodynamics
 
The phase space of thermodynamic values used in the first law (<ref>) and the very same relation (<ref>) can be generalized if one assumes that the cosmological constant Ł is not a fundamental one but a parameter which can be varied. The corresponding variation of gravitational Lagrangian over parameter Ł shows that the respective variation has a dimension of volume and might be called as a thermodynamic volume and it means that the parameter Ł might be identified with the thermodynamic pressure <cit.>. The exact definitions of this pressure and conjugate thermodynamic volume are as follows:

    P=-Λ/8π,    V=(∂ M/∂ P)_S, Q^+_.

Having used the definition of the thermodynamic volume (<ref>) we can write the evident form for it:

    V=ω_n-1(r^n_+/n-+Ł/2[∑^n-1/2_j=0(-1)^jd^2jr^n-2j_+/n-2j+(-1)^n+1/2d^narctan(r_+/d)]+.
    .2^p-2(2p-1)^2 q^2p/(n-2p)r^n-2p(n-1)/2p-1_+_2F_1(1,p(n-1)/2p-1-n/2;p(n-1)/2p-1-n/2+1;-d^2/r^2_+))

for odd n and

    V=ω_n-1(r^n_+/n-+Ł/2[∑^n/2-1_j=0(-1)^jd^2jr^n-2j_+/n-2j+(-1)^n/2d^nln(r^2_+/d^2+1)]+.
    .2^p-2(2p-1)^2 q^2p/(n-2p)r^n-2p(n-1)/2p-1_+_2F_1(1,p(n-1)/2p-1-n/2;p(n-1)/2p-1-n/2+1;-d^2/r^2_+))

for even n correspondingly. We point out that the only substantial difference in the given relations is the ∼arctan(r_+/d) term for odd n and the term ∼ln(r^2_+/d^2+1) for even n.

We also consider an additional intensive variable, which, as was shown <cit.> is necessary when one tries to obtain Smarr relation. This additional variable and its conjugate value can be written as follows:

    Π=/8π,   Ψ=(∂ M/∂Π)_S,Q^+_,P

Having used the latter relation we obtain the evident form for the value Ψ, namely for odd n we arrive at:

    Ψ=ω_n-1/2[(^2-Ł^2^2)/2^2(∑^n-1/2_j=0(-1)^jd^2jr^n-2j_+/n-2j+(-1)^n+1/2d^narctan(r_+/d))-(+Ł)^2/4^2×.
    .(∑^n-1/2_j=1(-1)^j2jd^2j/n-2jr^n-2j_++(-1)^n+1/2d^n(narctan(r_+/d)-dr_+/r^2_++d^2))-Ψ_1],

whereas for even n we derive:

    Ψ=ω_n-1/2[(^2-Ł^2^2)/2^2(∑^n/2-1_j=0(-1)^jd^2jr^n-2j_+/n-2j+(-1)^n/2d^n/2ln(r^2_+/d^2+1))-(+Ł)^2/4^2×.
    .(∑^n/2-1_j=1(-1)^j2jd^2j/n-2jr^n-2j_++(-1)^n/2d^n/2(nln(r^2_+/d^2+1)-2r^2_+/r^2_++d^2))-Ψ_1],

where we have used the following notation in order to simplify the previous two relations:
 
    Ψ_1=2^p-1(2p-1)^2q^2p/^2(2p-n)r^n-2p(n-1)/2p-1_+[Ł^2_2F_1(1,p(n-1)/2p-1-n/2;p(n-1)/2p-1-n/2+1;-d^2/r^2_+)+.
    .(2p-n)/(n+2p-2)(+Ł) d^2/r^2_+_2F_1(2,p(n-1)/2p-1-n/2+1;p(n-1)/2p-1-n/2+2;-d^2/r^2_+)]+
    2^2p-3(2p-1)^3^2 q^4p/^2(4p-2pn-n)r^n-4p(n-1)/2p-1_+[_2F_1(1,2p(n-1)/2p-1-n/2;2p(n-1)/2p-1-n/2+1;-d^2/r^2_+)+.
    .(4p-2pn-n)d^2/(2pn+n-2)r^2_+_2F_1(2,2p(n-1)/2p-1-n/2+1;2p(n-1)/2p-1-n/2+2;-d^2/r^2_+)].

It should be pointed out that the main difference in the relations (<ref>) and (<ref>) is of similar character as it was for the thermodynamic volume, related to the fact that odd and even dimensions have inverse tangential and logarithmic terms respectively.

Having extended thermodynamic phase space by additional variables such as P and Π we can write the extended first law in the form:

    δ M=Tδ S+Φ^+_δ Q^+_+Φ_qδ Q+Vδ P+ΨδΠ.

Now, having utilized all the new thermodynamic variables and taking into account the equality (<ref>) we obtain the Smarr relation, which can be represented as follows:

    (n-2)M=(n-1)TS-2VP-2ΨΠ+1/p(3p-pn-1)Φ_qQ,
 
and as a consequence for linear case (p=1) we have:

    (n-2)M=(n-1)TS-2VP-2ΨΠ+(2-n)Φ_qQ,

which is typical for a black hole with linear Maxwell field. It should be emphasized that the Smarr relation (<ref>) we have just obtained to great extent is grounded on the chosen form of scalar potential (<ref>) and this form seems to be preferable among the other possible variants.  


 §.§ Gibbs free energy and heat capacity

Since we have introduced thermodynamic pressure and volume by the relations (<ref>) the mass function M would play the role of enthalpy <cit.>. It is known that more convenient function to analyze thermodynamic behaviour of a system, especially in case there is some critical behaviour, is the Gibbs free energy which can be introduced in the following way:

    G=M-TS

Utilizing the latter relation we can write the explicit relation for the Gibbs free energy, namely for odd n we can write:

    G=ω_n-1/16π( r^n-2_++2Ł/n(n-1)r^n_++(+Ł)^2/2(n-1)[∑^n-1/2_j=0(-1)^j(2j-1)d^2j/n-2jr^n-2j_++(-1)^n+1/2d^n×..
    ..((n-1)arctan(r_+/d)-dr_+/r^2_++d^2)]+2^p(2p-1)(4p-2pn-1)q^2p/(n-1)(2p-n)r^n-2p(n-1)/2p-1_++G_1),

whereas for even n we obtain:

    G=ω_n-1/16π( r^n-2_++2Ł/n(n-1)r^n_++(+Ł)^2/2(n-1)[∑^n/2-1_j=0(-1)^j(2j-1)d^2j/n-2jr^n-2j_++(-1)^n/2d^n/2×..
    ..((n-1)ln(r^2_+/d^2+1)-2r^2_+/r^2_++d^2)]+2^p(2p-1)(4p-2pn-1)q^2p/(n-1)(2p-n)r^n-2p(n-1)/2p-1_++G_1),

and here we have used the following notation:

    G_1=-2^p-1(2p-1)(4p-2pn-1)(+Ł)q^2p/(n-1)(2p-n)r^n-2p(n-1)/2p-1_+_2F_1(1,p(n-1)/2p-1-n/2;p(n-1)/2p-1-n/2+1;-d^2/r^2_+)
    +2^p(2p-1)^2(+Ł)d^2q^2p/(n-1)(2p+n-2)r^n-2-2p(n-1)/2p-1_+_2F_1(2,p(n-1)/2p-1-n/2+1;p(n-1)/2p-1-n/2+2;-d^2/r^2_+)
    +2^2p-3(2p-1)^2(4pn-6p+1) q^4p/(n-1)(4p-2pn-n)r^n-4p(n-1)/2p-1_+_2F_1(1,2p(n-1)/2p-1-n/2;2p(n-1)/2p-1-n/2+1;-d^2/r^2_+)
    +2^2(p-1)(2p-1)^3 d^2q^2p/(n-1)(2pn+n-2)r^n-2-4p(n-1)/2p-1_+_2F_1(2,2p(n-1)/2p-1-n/2+1;2p(n-1)/2p-1-n/2+2;-d^2/r^2_+).

For the linear case (p=1) the Gibbs free energy can be written in a bit simpler from, in particular for odd n we can write:

    G=ω_n-1/16π( r^n-2_++2Ł/n(n-1)r^n_++2(2n-3)/(n-1)(n-2)q^2r^2-n_++(+Ł)^2/2(n-1)[∑^n-1/2_j=0(-1)^j(2j-1)d^2j/n-2jr^n-2j_++..
    ..(-1)^n+1/2d^n((n-1)arctan(r_+/d)-dr_+/r^2_++d^2)]+(+Ł)q^2/(n-1)[∑^n-5/2_j=0(-1)^j(2n-2j-5)/(2j+4-n)d^2(j+1)r^2j+4-n_++..
    ..(-1)^n-3/2/d^n-2((n-1)arctan(r_+/d)-dr_+/r^2_++d^2)]+ q^4/2(n-1)[∑^3n-7/2_j=0(-1)^j(4n-2j-7)/(2j+6-3n)d^2(j+1)r^2j+6-3n_++..
    ..(-1)^3n-5/2/d^3n-4((n-1)arctan(r_+/d)-dr_+/r^2_++d^2)])
 
and for even n it takes the following form:

    G=ω_n-1/16π( r^n-2_++2Ł/n(n-1)r^n_++2(2n-3)/(n-1)(n-2)q^2r^2-n_++(+Ł)^2/2(n-1)[∑^n/2-1_j=0(-1)^j(2j-1)d^2j/n-2jr^n-2j_++..
    ..(-1)^n/2d^n((n-1)/2ln(r^2_+/d^2+1)-r^2_+/r^2_++d^2)]+(+Ł)q^2/(n-1)[∑^n-6/2_j=0(-1)^j(2n-2j-5)/(2j+4-n)d^2(j+1)r^2j+4-n_++..
    ..(-1)^n-2/2/d^n-2((n-1)/2ln(d^2/r^2_++1)+d^2/r^2_++d^2)]+ q^4/2(n-1)[∑^3n-8/2_j=0(-1)^j(4n-2j-7)/(2j+6-3n)d^2(j+1)r^2j+6-3n_++..
    ..(-1)^3n/2/d^3n-4((n-1)/2ln(d^2/r^2_++1)+d^2/r^2_++d^2)])
 
It should be emphasized that the second sums in the relations (<ref>) and (<ref>) (namely the sums ∑^n-5/2_j=0 and ∑^n-6/2_j=0) do not appear in the given above relations when n=3 and n=4 respectively.

 
We can also write the expressions for the Gibbs free energy in conformal case, namely, for odd n we have:

    G=ω_n-1/16π( r^n-2_++2Ł/n(n-1)r^n_++2^n-3/4nq^n+1/2/r_++2^n-7/4(+Ł)/ dq^n+1/2((n-1)arctan(r_+/d)-r_+d/r^2_++d^2).
    .(+Ł)^2/2(n-1)[∑^n-1/2_j=0(-1)^j(2j-1)d^2j/n-2jr^n-2j_++(-1)^n+1/2d^n((n-1)arctan(r_+/d)-r_+d /r^2_++d^2)]+.
    .2^n-9/2(n-1)/q^n+1[∑^n-1/2_j=0(-1)^j(2(n-j)-1)/2j-nr^2j-n_+/d^2(j+1)+(-1)^n+1/2/d^n+2((n-1)arctan(r_+/d)-r_+d/r^2_++d^2)]).

It should be stressed that when n=3 the latter relation coincides with corresponding relation obtained from the relation for linear field (<ref>) as it should be. For even n the Gibbs free energy takes the form:

    G=ω_n-1/16π( r^n-2_++Ł/n(n-1)r^n_++2^n-3/4nq^n+1/2/r_++2^n-7/4(+Ł)/ dq^n+1/2((n-1)arctan(r_+/d)-r_+d/r^2_++d^2)+.
    .(+Ł)^2/2(n-1)[∑^n/2-1_j=2(-1)^j(2j-1)d^2j/n-2jr^n-2j_++(-1)^n/2d^n((n-1)/2ln(r^2_+/d^2+1)-r^2_+/r^2_++d^2)]+.
    .2^n-9/2(n-1)/q^n+1[∑^n/2-1_j=0(-1)^j(2(n-j)-1)/2j-nr^2j-n_+/d^2(j+1)+(-1)^n/2+1/d^n+2((n-1)/2ln(d^2/r^2_++1)+d^2/r^2_++d^2)]).

The obtained above relations have some similarity with the linear case (p=1). The Fig.[<ref>]  represents the behaviour of the Gibbs free energy for linear p=1 and conformal  p=n+1/4 cases.  They show some similarity namely both of them have smooth behaviour for the pressures above the critical one (P_c=-Ł_c/8π). The critical value of pressure P_c (or cosmological constant Ł_c) can be defined as the value when nonmonotonous  behvaviour in the dependence T=T(r_+) disappears (see the Fig.[<ref>]). For the pressures above this critical one the function T=T(r_+) is monotonous, whereas for the pressures slightly below the critical values a specific maximum and minimum appear (Fig.[<ref>]). It should be pointed out that in linear case (p=1) there  are some similarities with the behaviour of the function G=G(T,P) for real gases and liquids, namely for the pressures above the critical value we have bigger values of the Gibbs free energy for the same temperatures, whereas for the conformal case the situation is opposite. Another important point we should emphasize is the fact that in both cases (p=1 and p=(n+1)/4) for the pressures substantially below the critical one a specific loop on the graphs of the function G=G(T,P) appears, this situation in general is similar to the case of standard General Relativity or real condensed systems which are described by a Van der Waals-like equations of state. We also mention that for Van der Waals-systems (both gravitational and condensed matter) a first order phase transition appears for the pressures below the critical one and specific swallow-tail behaviour of the Gibbs free energy takes place. This swallow-tail loop appears immediately below the critical pressure and increase with corresponding decreasing of the pressure. We also point out that for some gravitational systems, namely for black holes with dilaton field additional phase transition of the zeroth order takes place. This zeroth order phase transition takes place when the pressure is below the critical one but the swallow-tail loop has not formed yet. In quite general features we have similar situation here, but to make a conclusion about type of phase transition deeper investigation of thermodynamic functions should be made and it will be performed elsewhere.

The other important function which brings considerable information about thermodynamic behaviour and especially about the phase transitions that might occur for a given system is the heat capacity:

    C_P=T(∂ S/∂ T)_P,Π,Q^+_,Q=T(∂ S/∂ r_+)_P(∂ r_+/∂ T)_P

It  should be pointed out that in the given above relation for the heat capacity we keep the parameters Q^+_, Q and Π fixed and this fact allows us to write the evident form for the heat capacity as follows:

    C_p=(n-1)ω_n-1r^n-2_+/4√(U(r_+)W(r_+))U'(r_+)((2-n)/r^2_+(1-(+Ł)^2/4^2)+(-Ł)^2/2(n-1)-2^p(4p-2pn-1)/n-1×.
    .r^-2p(n-1)/2p-1_++(+Ł)^2/2(n-1)[∑^n-1-σ/2_j=2(-1)^j(1-2j)d^2j/r^2j_++(-1)^n+1-σ/2d^n+1-σ((2+σ-n)d^2+(σ-n)r^2_+)/r^n-1-σ_+(r^2_++d^2)^2]+.
    .2^p-1(+Ł)/(n-1)q^2pr^-2p(n-1)/2p-1_+[(4p-2pn-1)_2F_1(1,p(n-1)/2p-1-n/2;p(n-1)/2p-1-n/2+1;-d^2/r^2_+)+..
    ..2(4p^2-1)(n-1)d^2/(n+2p-2)r^2_+_2F_1(2,p(n-1)/2p-1-n/2+1;p(n-1)/2p-1-n/2+2;-d^2/r^2_+)-8(2p-1)^2d^4/(n+6p-4)r^4_+×..
    .._2F_1(3,p(n-1)/2p-1-n/2+2;p(n-1)/2p-1-n/2+3;-d^2/r^2_+)]+2^2p-3(2p-1)^2/(n-1)q^4pr^-4p(n-1)/2p-1_+[(6p-4pn-1)/2p-1×..
    .._2F_1(1,2p(n-1)/2p-1-n/2;2p(n-1)/2p-1-n/2+1;-d^2/r^2_+)+2(6p+1)(n-1)d^2/(2pn+n-2)r^2_+_2F_1(2,2p(n-1)/2p-1-n/2+1;...
    ...2p(n-1)/2p-1-n/2+2;-d^2/r^2_+)-8(2p-1)d^4r^-4_+/(2pn+n+4p-4)_2F_1(3,2p(n-1)/2p-1-n/2+2;2p(n-1)/2p-1-n/2+3;-d^2/r^2_+)])^-1.

We note that for simplification of the latter relation we have not substituted the explicit relations for the derivative U'(r_+) and for the square root of the product of the metric functions  √(U(r_+)W(r_+)). We also write the expression for the heat capacity in case of linear gauge field:

    C_p=(n-1)ω_n-1r^n-2_+/4√(U(r_+)W(r_+))U'(r_+)((2-n)/r^2_+(1-(+Ł)^2/4^2)+(-Ł)^2/2(n-1)-2(3-2n)/n-1×.
    . q^2r^2(1-n)_++(+Ł)^2/2(n-1)[∑^n-1-σ/2_j=2(-1)^j(1-2j)d^2j/r^2j_++(-1)^n+1-σ/2d^n+1-σ/r^n-1-σ_+((2+σ-n)d^2+(σ-n)r^2_+)/(r^2_++d^2)^2]+.
    .(+Ł)q^2/(n-1)[∑^n-5-σ/2_j=0(-1)^j(2j+5-2n)/d^2(j+1)r^2(j+2-n)_++(-1)^n-3-σ/2d^3+σ-n/r^n-1+σ_+((2-σ-n)d^2-(n+σ)r^2_+)/(r^2_++d^2)^2]+.
    . q^4/2(n-1)[∑^3n-7-σ/2_j=0(-1)^j(2j+7-4n)/d^2(j+1)r^2(j+3-2n)_++(-1)^3n-5-σ/2d^5-3n+σ/r^n-1+σ_+((2-σ-n)d^2-(n+σ)r^2_+)/(r^2_++d^2)^2])^-1.


For the conformal case the heat capacity can be written in the form:

    C_P=(n-1)ω_n-1r^n-2_+/4√(U(r_+)W(r_+))U'(r_+)((2-n)/r^2_+(1-(+Ł)^2/4^2)+(-Ł)^2/2(n-1)+.
    .(+Ł)^2/2(n-1)[∑^n-1-σ/2_j=2(-1)^j(1-2j)d^2j/r^2j_++(-1)^n+1-σ/2d^n+1-σ/r^n-1-σ_+((2+σ-n)d^2+(σ-n)r^2_+)/(r^2_++d^2)^2]+.
    .2^n-3/4n q^n+1/2/r^n+1_++2^n-7/4(+Ł)q^n+1/2/ r^n-1_+(r^2_++d^2)^2((2-n)d^2-nr^2_+)+2^n-9/2(n-1) q^n+1/×.
    .[∑^n-1-σ/2_j=0(-1)^j(2(j-n)+1)r^2(j-n)_+/d^2(j+1)+(-1)^n+1-σ/2/d^n+1-σ((2-σ-n)d^2-(σ+n)r^2_+)/r^n-1+σ_+(r^2_++d^2)^2])^-1.

The comparison of the given above heat capacities is shown on the Fig.[<ref>]. The behaviour of the heat capacities is very similar, namely both of them have discontinuity which separate stable and unstable domains, for the horizon radii larger than these specific ones where discontinuity takes place the heat capacities are positive and it means that black holes are stable, whereas the black holes with radii smaller than these specific values are unstable because of the negative values of the heat capacities. It should be noted that for both cases p=1 and p=(n+1)/4 numerical values of heat capacities are very close in case all the fixed parameters take the same value. Since the function T(r_+) might have two extrema points it means that the heat capacity possesses two points of discontinuity and one of those two points is depicted on the Fig.[<ref>], the second point of discontinuity takes place for the smaller values of the horizon radius r_+ and again it separates stable and unstable domains. We point out that the domain of instability is located in between those two points of discontinuity. When the pressure P (or absolute value of Ł) goes up the points of discontinuity become closer and for some specific value of pressure they merge and get transformed into a point of a local maximum of the heat capacity (Fig.[<ref>]). Further increasing of the pressure gives rise to the diminishing with the following disappearance  of the maximum, so for high pressures when the temperature of the black hole becomes monotonous function T=T(r_+) the heat capacity gets transformed into a monotonous function also. 




§ CONCLUSIONS

In our work we have considered black hole solutions in a scalar-tensor theory with nonminimal derivative coupling between gravity and scalar field  and nonlinear electromagnetic field minimally coupled to gravity. We have taken the power law nonlinearity for the electromagnetic field Lagrangian in order to have possibility to obtain conformally invariant electrodynamics, but in our work we do not fix the parameter p which characterizes degree of nonlinearity. As a result we have obtained solutions valid for different values of the parameter p, including p=1 which was previously derived in the paper <cit.> and p=(n+1)/4 as particular cases. Thus, the solutions we have obtained in our work are a generalization of the solutions derived in <cit.>. It should be pointed out that  the parameter p is not completely arbitrary, to have reasonable behaviour of the electric field at the infinity and well-defined thermodynamic relations the parameter p  should satisfy the following condition 1/2<p<n/2 and it is easy to see that linear p=1 and conformal p=(n+1)/4 cases fulfill the mentioned above condition.

 In this work we focus on static solutions and since we take into account cosmological constant Ł it gives us an opportunity to obtain topological black holes with different types of horizon surfaces, namely spherical (=1), flat (=0) and hyperbolic (=-1) horizons. The obtained three types of solution have some common features, whereas the other ones might be distinct, namely at large distances (r→∞) all of them show anti-de Sitterian behaviour, at the same time for small distances (r→ 0) all three solutions represent different types of geometry. It should be pointed out that for the flat type geometry the obtained black hole's solution is the simplest one and is expressed in terms of elementary power-law functions of r for any allowed value of p, whereas for other two types of geometry the general solution can be represented by hypergeometric functions which reduce to elementary ones only when the parameters of the hypergeometric functions are integer, in particular they are integer for linear p=1 and conformal cases p=(n+1)/4. Another important moment we should mention is the fact that in the regime of strong nonminimal coupling (when the term with factor  is relatively large in comparison with the term with the factor ) the behaviour of metric functions for flat (=0) and nonflat (=± 1) geometries differs substantially. To investigate the question about  singularities of the obtained metrics we have considered the Kretschmann scalar (<ref>) which shows that there is just a coordinate singularity at the black hole's horizon r_+ and the only physical singularity is at the origin. It should be emphasized here that in the vicinity of the origin of coordinates the behaviour of the Kretschmann scalar is ∼1/r^4 for all types of solutions and arbitrary n and all allowed values of p, and the leading term of this scalar (<ref>) does not depend on the integration constants μ and q which are related to black hole's mass and charge respectively, this fact is specific peculiarity that takes place for the obtained black hole solutions and does not happen in the framework of standard General Relativity. At the infinity the Kretschmann scalar takes finite value which is completely defined by the parameters of coupling  and  and the dimension of space n, whereas the cosmological constant is not present here and this fact is the major difference from the situation in standard General Relativity where the presence of the constant Ł means that it appears in the Kretschmann scalar at the infinity.
 
 We have also obtained the potential of the gauge field and using Gauss law we have calculated black hole's charge. Here we point out that the electric potential, namely its change between the horizon and the infinity and the charge are extremely important for constructing of the first law of black hole's thermodynamics. It is worth being stressed that the mentioned above conditions imposed on the parameter p are derived after careful analysis of the behaviour of the electric potential which should not diverge at the infinity. We also note that the same constraints on the parameter p are obtained when we use Wald procedure to derive the first law of black hole's thermodynamics.
 
 Using standard geometrical relation we calculate black hole's temperature in general case and we also derive temperature for linear and conformal cases. The quite general trait of the temperature is the fact that it goes up almost linearly when the radius of horizon r_+ becomes large enough, similar behaviour of the temperature takes place in standard General Relativity when one takes into account the cosmological constant Ł, but in our case instead of “pure” cosmological constant Ł an “effective” cosmological constant appears which depends on the “pure” value and the parameters of coupling  and . For small radii of horizon the temperature becomes negative. The Fig.[<ref>] demonstrates the dependences T=T(r_+) for some intermediate values of radii of horizon r_+, these graphs show that the temperature might be a nonmonotonous function of r_+ if the cosmological constant Ł (or the pressure we identify with) is lower than the specific value Ł_c, whereas for the cosmological constant is larger than Ł_c the temperature becomes monotonic. The described behaviour of the temperature is similar to the situation that takes place for real systems in condensed matter physics, namely real gases and liquids. The nonmonotonic behaviour of temperature gives rise to the appearance of  phase transitions which separate stable and unstable phases.
 
 To obtain entropy of the black hole we have applied Wald's approach <cit.> which has shown its efficiency in scalar-tensor theories of gravity <cit.>. It should be pointed out that the relation for the entropy we have obtained differs from the well-known relation of standard General Relativity where the entropy is defined as a quarter of horizon area surface, namely we have obtained the relation (<ref>) which has this quarter of the area multiplied by some specific factor, caused by the nonminimal coupling and the obtained general relation for the entropy is the same for all types of the black holes we investigate. To derive the relation for entropy we have also introduced additional specific scalar “charge” Q^+_ as it was proposed by Feng et al. <cit.> but this scalar “charge” can be chosen in different ways and similarly as it was done in our earlier work <cit.> we take it in a different form than it was used in the original paper <cit.>. The evident form of this scalar “charge” we have taken allows us to derive Smarr relation (<ref>)  easily. Having introduced the entropy and the scalar charge Q^+_ together with conjugate potential we have obtained  the first law of black hole thermodynamics. 
 
Having used the concept of thermodynamic pressure, which is related to the cosmological constant Ł <cit.> we have developed the so called extended thermodynamics. Namely, we have derived the generalized first law (<ref>) and the mentioned above Smarr relation (<ref>). Having made Legendre transformation we have obtained the Gibbs free energy which is extremely useful for analysis of thermodynamic behaviour of the black hole. The Gibbs free energy has monotonous behaviour for the pressures above some specific value, which depends on the dimension of space, parameters of coupling  and , parameter of nonlinearity p, black hole's charge q and type of geometry of the horizon. When the pressure decreases the Gibbs free energy becomes nonmonotonous and finally a closed loop appears that has some similarities with closed loop of a swallow-tail form which is typical for a Van der Waals systems or charged black holes in standard General Relativity <cit.>. The appearance of the loops gives rise to the conclusion about the existence of phase transitions and rich phase behaviour, but these issues need more detailed examination and will be considered elsewhere. Finally, we have calculated heat capacity for the black holes, namely we have written the evident form for it in general case as well as in particular cases linear p=1 and conformal p=(n+1)/4 electrodynamics. The obtained relations and as a consequence the corresponding graphs for them show some common features, namely for small pressures the heat capacity has two discontinuity points which separates stable and unstable phases, the increase of the pressure leads to merging of these discontinuities and its transformation into a peak which vanishes with corresponding increase of the pressure.  
 

§ ACKNOWLEDGEMENTS

This work was partly supported by Project FF-30F (No. 0116U001539) from the Ministry of Education and Science of Ukraine.
99
Will_LRR2014 C. M. Will, Liv. Rev. Rel. 17, 4 (2014).
Berti_CQG2015 E. Berti, E. Barausse, V. Cardoso, et al. Class. Quant. Grav. 32, 243001 (2015). 
Abbott_PRL16 B. P. Abbott, Phys. Rev. Lett. 116, 061102 (2016).
Abbott_PRL16_2 B. P. Abbott, Phys. Rev. Lett. 116, 241103 (2016).
Abbott_PRL17 B. P. Abbott, Phys. Rev. Lett. 118, 221101 (2017).
Clifton_PhysRept2012 T. Clifton, P. G. Ferreira, A. Padilla, C. Skordis, Phys. Rept. 513, 1 (2012).
Heisenberg_2018 L. Heisenberg, arxiv:1807.01725
Horndeski_IJTP74 G. W. Horndeski, Int. Journ. Theor. Phys. 10, 363, (1974).
Deffayet_PRD09 C. Deffayet, S. Deser, G. Esposito-Farese, Phys. Rev. D 80, 064015 (2009).
Kobayashi_PTEP11 T. Kobayashi, M. Yamaguchi, J. Yokoyama, PTEP 126, 511 (2011).
Acoleyen_PRD11 K. Van Acoleyen, J. Van Doorsselaere, Phys. Rev. D 83, 084025 (2011).
Charmousis_LNP15 C. Charmousis, Lect. Notes Phys. 892, 25 (2015).
Metsaev_NPB87 R. R. Metsaev, A. A. Tseytlin, Nucl. Phys. B 293, 385 (1987).
Meissner_PLB97 K. A. Meissner, Phys. Lett. B 392, 298 (1997).
Cartier_PRD01 C. Cartier, J.-C. Hwang, E. J. Copeland, Phys. Rev. D 64, 103504 (2001).
deRham_PRD11 C. de Rham, L. Heisenberg, Phys. Rev. D 84, 043503 (2011). 
Deffayet_PRD10 C. Deffayet, S. Deser, G. Esposito-Farese, Phys. Rev. D 82, 061501 (2010).
Padilla_JHEP13 A. Padilla, V. Sivanesan, JHEP 04, 032 (2013).
Charmousis_JHEP14 C. Charmousis, T. Kolyvaris, E. Papantonopoulos, M. Tsoukalas, JHEP 07, 085 (2014).
Ohashi_JHEP15 S. Ohashi, N. Tanahashi, T. Kobayashi, M. Yamaguchi, JHEP 07 008 (2015).
Zumala_PRD14 M. Zumalacárregui, J. García-Bellido, Phys. Rev. D 89, 064046 (2014).
Gleyzes_PRL15 J. Gleyzes, D. Langlois, F. Piazza, F. Vernizzi, Phys. Rev. Lett. 114, 211101 (2015).
Gleyzes_JCAP15 J. Gleyzes, D. Langlois, F. Piazza, F. Vernizzi, JCAP 1502, 018 (2016).
Langlois_JCAP16 D. Langlois, K. Noui, JCAP 1602, 034 (2016).
Chrisostomi_JCAP16 M. Crisostomi, K. Koyama, G. Tasinato, JCAP 1604, 044 (2016).
BenAchour_PRD16 J. Ben Achour, D. Langlois, K. Noui, Phys. Rev. D 93, 124005 (2016).
BenAchour_JHEP16 J. Ben Achour, M. Crisostomi, K. Koyama, et al., JHEP 1612, 100 (2016).
Sushkov_PRD09 S. V. Sushkov, Phys. Rev. D 80, 103505 (2009).
Saridakis_PRD10 E. N. Saridakis, S. V. Sushkov, Phys. Rev. D 81, 083510 (2010).
Granda_JCAP10 L. N. Granda, W. Cadona, JCAP 07, 021 (2010).
Gao_JCAP10 C. Gao, JCAP 06, 023 (2010).
Sadjadi_PRD11 H. Mohseni Sadjadi, Phys. Rev. D 83, 107301 (2011).
Skugoreva_PRD13 M. A. Skugoreva, S. V. Sushkov, A. V. Toporensky, Phys. Rev. D 88, 083539 (2013); M. A. Skugoreva, S. V. Sushkov, A. V. Toporensky, Phys. Rev. D 88, 109906(E) (2013).
Sushkov_PRD12 S. V. Sushkov, Phys. Rev. D 85, 123520 (2012).
Germani_PRL10 C. Germani, A. Kehagias, Phys. Rev. Lett. 105, 011302 (2010).
Germani_PRL11 C. Germani, A. Kehagias, Phys. Rev. Lett. 106, 161302 (2011).
Sadjadi_JCAP13 H. M. Sadjadi, P. Goodarzi, JCAP 1302, 038 (2013); H. M. Sadjadi, P. Goodarzi, JCAP 1307, 038 (2013); 
Dalianis_JCAP17 I. Dalianis, G. Koutsoumbas, K. Nterkis, E. Papantonopoulos, JCAP 1702, 027 (2017).
Feng_PLB14 K. Feng, T. Qiu, Y.-S. Piao, Phys. Lett. B 729, 99 (2014).
Feng_PRD14 K. Feng, T. Qiu, Phys. Rev. D 90, 123508 (2014).
Qiu_EPJC17 T. Qiu, K. Feng, EPJC 77, 687 (2017).
Kennedy_PRD17 J. Kennedy, L. Lombriser, A. Taylor, Phys. Rev. D 96, 084051 (2017).
Hou_EPJC17 S.  Hou, Y. Gong, Eur. Phys. J. C 78, 247 (2018).
Crisostomi_PRD18 M. Crisostomi, K. Koyama, Phys. Rev. D 97, 021301 (2018); M. Crisostomi, K. Koyama, Phys. Rev. D 97, 084004 (2018).
Langlois_PRD18 D. Langlois, R. Saito, D. Yamauchi, K. Noui, Phys. Rev. D 97, 061501 (2018).
Rinaldi_PRD12 M. Rinaldi, Phys. Rev. D 86, 084048 (2012).
Minamitsuji_PRD14 M. Minamitsuji, Phys. Rev. D 89, 064017, (2014).
Anabalon_PRD14 A. Anabalon, A. Cisterna, J. Oliva, Phys. Rev. D 89, 084050 (2014).
Cisterna_PRD14 A. Cisterna, C. Erices, Phys. Rev. D 89, 084038 (2014).
Kobayashi_PTEP14 T. Kobayashi, N. Tanahashi, Prog. Theor. Exp. Phys. 2014, 073E02 (2014).
Babichev_JHEP14 E. Babichev, C. Charmousis, JHEP 08, 106 (2014).
Cisterna_PRD16 A. Cisterna, T. Delsate, L. Ducobu, M. Rinaldi, Phys. Rev. D 93, 084046 (2016).
Maselli_PRD15A. Maselli, H. O. Silva, M. Minamitsuji, E. Berti, Phys. Rev. D 92, 104049 (2015).
Stetsko M. M. Stetsko, arXiv:1811.05030.
Bravo-Gaete_PRD14 M. Bravo-Gaete, M. Hassaine, Phys. Rev. D 90, 024008 (2014).
Giribet_PRD15 G. Giribet, M. Tsoukalas, Phys. Rev. D 92, 064027 (2015).
Clement_arxiv G. Clément, Kh. Nouicer, arXiv:1805.09796.
Feng_JHEP15 X.-H. Feng, H.-S. Liu, H. Lu, C. N. Pope, JHEP 1511, 176 (2015).
Feng_PRD16 X.-H. Feng, H.-S. Liu, H. Lu, C. N. Pope, Phys. Rev. D 93 044030 (2016). 
Charmousis_PRD15 Ch. Charmousis, M. Tsoukalas, Phys. Rev. D 92, 104050 (2015).
Cisterna_PRD15 A. Cisterna, M. Cruz, T. Delsate, J. Saavedra, Phys. Rev. D 92, 104018 (2015).
Takahashi_PRD17 K. Takahashi, T. Suyama, Phys. Rev. D 95, 024034 (2017).
Tretyakova_CQG17 D. A. Tretyakova, K. Takahashi, Class. Quant. Grav. 34, 175007 (2017).
Ganguly_CQG18 A. Ganguly, R. Gannouji, M. Gonzalez-Espinoza, C. Pizarro-Moya, Class. Quant. Grav. 35, 145008  (2018).
Babichev_PRL18 E. Babichev, C. Charmousis, G. Esposito-Farese, A. Lehébel, Phys. Rev. Lett. 120, 241101 (2018).
Benkel_arxiv R. Benkel, N. Franchini, M. Saravani, T. P. Sotiriou, arXiv:1806.08214.
Cisterna_PRD15 A. Cisterna, T. Delsate, M. Rinaldi, Phys. Rev. D 92, 044050 (2015).
Brihaye_PRD16 Y. Brihaye, A. Cisterna, C. Erices, Phys. Rev. D 93, 124057 (2016).
Verbin_PRD18 Y. Verbin, Y. Brihaye, Phys. Rev. D 97, 044046 (2018). 
Sotiriou_PRL14 T. P. Sotiriou, S.-Y. Zhou, Phys. Rev. Lett. 112, 251102 (2014);  T. P. Sotiriou, S.-Y. Zhou, Phys. Rev. D 90, 124063 (2014).
Antoniou_PRL18 G. Antoniou, A. Bakopoulos, P. Kanti, Phys. Rev. Lett. 120, 131102 (2018); G. Antoniou, A. Bakopoulos, P. Kanti, Phys. Rev. D 97, 084037 (2018).
Hassaine_PRD07 M. Hassaine, C. Martinez, Phys. Rev. D 75, 027502 (2007); M. Hassaine, C. Martinez, Class. Quant. Grav. 25, 195023, (2008).
Maeda_PRD09 H. Maeda, M. Hassaine, C. Martinez, Phys. Rev. D 79, 044012 (2009);
Gurtug_PRD12 O. Gurtug, S. H. Mazharimousavi, M. Halilsoy, Phys. Rev. D 85, 104004 (2012).
Kord_PRD15 M. K. Zangeneh, A. Sheykhi, M. H. Dehghani, Phys. Rev. D 91, 044035 (2015).
Kord_PRD15_2 M. K. Zangeneh, A. Sheykhi, M. H. Dehghani, Phys. Rev. D 92, 024050 (2015); M. K. Zangeneh, M. H. Dehghani, A. Sheykhi, Phys. Rev. D 92, 104035 (2015).
Stetsko_2 M. M. Stetsko, arXiv:1810.00061.
bateman H. Bateman, A. Erdelyi, Higher transcendental functions, vol.1, N.-Y., McGraw-Hill (1953).
Gibbons_PRD77 G. W. Gibbons, S. W. Hawking, Phys. Rev. D 15, 2752 (1977).
Wald_PRD93 R. M. Wald, Phys. Rev. D 48, 3427 (1993).
Iyer_PRD94 V. Iyer, R. M. Wald, Phys. Rev. D 50, 846 (1994).
Liu_PLB14 H. S. Liu, H. Lu, Phys. Lett. B 730, 267 (2014).
Lu_JHEP15 H. Lu, C. N. Pope, Q. Wen, JHEP 1503, 165 (2015).
Liu_JHEP14 H. S. Liu, H. Lu, C.N. Pope, JHEP 1406, 109 (2014).
Fan_JHEP15 Z. Y. Fan, H. Lu, JHEP 1502, 013 (2015).
Fan_PRD15 Z. Y. Fan, H. Lu, Phys. Rev. D 91, 064009 (2015).
Liu_JHEP14_2 H. S. Liu, H. Lu, JHEP 1412, 071 (2014). 

Kastor_CQG09 D. Kastor, S. Ray , J. Traschen, Class. Quant. Grav. 26, 195011, (2009). 
Miao_EPJC16 Y.-G. Miao, Z.-M. Xu, EPJC 76, 638 (2016).
Kastor_CQG09 D. Kastor, S. Ray , J. Traschen, Class. Quant. Grav. 26, 195011, (2009). 
Kubiznak_JHEP12 D. Kubiznak, R. B. Mann, JHEP 07, 033 (2012).
Gunasekaran_JHEP12 S. Gunasekaran, R. B. Mann, D. Kubiznak, JHEP 11, 110 (2012).
Kubiznak_CQG2017 D. Kubiznak, R. B. Mann, M. Teo, Class. Quant. Grav. 34 063001 (2017).

