arxiv: v2 [hep-th] 24 Apr 2018

Similar documents
d 2 x 0a d d =0. Relative to an arbitrary (accelerating frame) specified by x a = x a (x 0b ), the latter becomes: d 2 x a d 2 + a dx b dx c

Geometry of the homogeneous and isotropic spaces

EM Boundary Value Problems

Gravitational Memory?

I. CONSTRUCTION OF THE GREEN S FUNCTION

( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.

arxiv: v1 [physics.pop-ph] 3 Jun 2013

BLACK HOLES IN STRING THEORY

Lecture 8 - Gauss s Law

Problems with Mannheim s conformal gravity program

2. Electrostatics. Dr. Rakhesh Singh Kshetrimayum 8/11/ Electromagnetic Field Theory by R. S. Kshetrimayum

Spherical Solutions due to the Exterior Geometry of a Charged Weyl Black Hole

Scattering in Three Dimensions

The Strain Compatibility Equations in Polar Coordinates RAWB, Last Update 27/12/07

On the integration of the equations of hydrodynamics

Compactly Supported Radial Basis Functions

The Schwartzchild Geometry

A New Approach to General Relativity

Vectors, Vector Calculus, and Coordinate Systems

An Exact Solution of Navier Stokes Equation

The Schwarzschild Solution

ASTR415: Problem Set #6

PROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.

Chapter 13 Gravitation

Liquid gas interface under hydrostatic pressure

Exceptional regular singular points of second-order ODEs. 1. Solving second-order ODEs

Pressure Calculation of a Constant Density Star in the Dynamic Theory of Gravity

arxiv: v2 [gr-qc] 18 Aug 2014

Problems with Mannheim s conformal gravity program

Chem 453/544 Fall /08/03. Exam #1 Solutions

A thermodynamic degree of freedom solution to the galaxy cluster problem of MOND. Abstract

Homework # 3 Solution Key

Fresnel Diffraction. monchromatic light source

ON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION. 1. Introduction. 1 r r. r k for every set E A, E \ {0},

Physics 2A Chapter 10 - Moment of Inertia Fall 2018

Perturbation theory and stability analysis for string-corrected black holes in arbitrary dimensions

working pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50

AST 121S: The origin and evolution of the Universe. Introduction to Mathematical Handout 1

Conformal transformations + Schwarzschild

Absorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere

arxiv:gr-qc/ v1 1 Sep 2005

4. Kruskal Coordinates and Penrose Diagrams.

From Gravitational Collapse to Black Holes

A NEW VARIABLE STIFFNESS SPRING USING A PRESTRESSED MECHANISM

MATH 415, WEEK 3: Parameter-Dependence and Bifurcations

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Problem Set 10 Solutions. r s

The Substring Search Problem

Chapter 3 Optical Systems with Annular Pupils

Contact impedance of grounded and capacitive electrodes

Does a black hole rotate in Chern-Simons modified gravity?

Nuclear size corrections to the energy levels of single-electron atoms

Math 2263 Solutions for Spring 2003 Final Exam

The Millikan Experiment: Determining the Elementary Charge

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Electromagnetism II September 15, 2012 Prof. Alan Guth PROBLEM SET 2

DOING PHYSICS WITH MATLAB COMPUTATIONAL OPTICS

Homework 7 Solutions

The Precession of Mercury s Perihelion

Supplementary Figure 1. Circular parallel lamellae grain size as a function of annealing time at 250 C. Error bars represent the 2σ uncertainty in

arxiv:hep-th/ v2 11 Nov 2004

arxiv: v1 [physics.gen-ph] 18 Aug 2018

Right-handed screw dislocation in an isotropic solid

Is there a magnification paradox in gravitational lensing?

6 PROBABILITY GENERATING FUNCTIONS

To Feel a Force Chapter 7 Static equilibrium - torque and friction

Physics 235 Chapter 5. Chapter 5 Gravitation

Analytical Solutions for Confined Aquifers with non constant Pumping using Computer Algebra

As is natural, our Aerospace Structures will be described in a Euclidean three-dimensional space R 3.

ON THE TWO-BODY PROBLEM IN QUANTUM MECHANICS

Classical Mechanics Homework set 7, due Nov 8th: Solutions

you of a spring. The potential energy for a spring is given by the parabola U( x)

8 Separation of Variables in Other Coordinate Systems

Physics 161 Fall 2011 Extra Credit 2 Investigating Black Holes - Solutions The Following is Worth 50 Points!!!

Physics 107 TUTORIAL ASSIGNMENT #8

arxiv:gr-qc/ v1 5 Jun 1996

Introduction to Nuclear Forces

A Relativistic Electron in a Coulomb Potential

Geodesic motion in Kerr spacetime

Question 1: The dipole

Review: Electrostatics and Magnetostatics

Single Particle State AB AB

Quantum Mechanics II

Exploration of the three-person duel

Vectors, Vector Calculus, and Coordinate Systems

On the ratio of maximum and minimum degree in maximal intersecting families

Encapsulation theory: the transformation equations of absolute information hiding.

Math 124B February 02, 2012

On the ratio of maximum and minimum degree in maximal intersecting families

= e2. = 2e2. = 3e2. V = Ze2. where Z is the atomic numnber. Thus, we take as the Hamiltonian for a hydrogenic. H = p2 r. (19.4)

3.1 Random variables

Higher Derivative Corrected Black Holes: Perturbative Stability and Absorption Cross Section in Heterotic String Theory

7.2. Coulomb s Law. The Electric Force

Abelian Symmetry Breaking and Modified Gravity

Hawking Radiation. Alex Fontana

1 Spherical multipole moments

Euclidean Figures and Solids without Incircles or Inspheres

! E da = 4πkQ enc, has E under the integral sign, so it is not ordinarily an

KOEBE DOMAINS FOR THE CLASSES OF FUNCTIONS WITH RANGES INCLUDED IN GIVEN SETS

Solution of a Spherically Symmetric Static Problem of General Relativity for an Elastic Solid Sphere

MATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE. We consider second order constant coefficient scalar linear PDEs on R n. These have the form

A Valid Finite Bounded Expanding Carmelian Universe Without Dark Matter

Transcription:

Black-Hole Solutions with Scala Hai in Einstein-Scala-Gauss-Bonnet Theoies axiv:1711.07431v2 [hep-th] 24 Ap 2018 G. Antoniou (a,b) 1, A. Bakopoulos (b) 2 and P. Kanti (b) 3 (a) School of Physics and Astonomy, Univesity of Minnesota, Minneapolis, MN 55455, USA (b) Division of Theoetical Physics, Depatment of Physics, Univesity of Ioannina, Ioannina GR-45110, Geece Abstact In the context of the Einstein-scala-Gauss-Bonnet theoy, with a geneal coupling function between the scala field and the quadatic Gauss-Bonnet tem, we investigate the existence of egula black-hole solutions with scala hai. Based on a pevious theoetical analysis, that studied the evasion of the old and novel no-hai theoems, we conside a vaiety of foms fo the coupling function (exponential, even and odd polynomial, invese polynomial, and logaithmic) that, in conjunction with the pofile of the scala field, satisfy a basic constaint. Ou numeical analysis then always leads to families of egula, asymptotically-flat black-hole solutions with non-tivial scala hai. The solution fo the scala field and the pofile of the coesponding enegy-momentum tenso, depending on the value of the coupling constant, may exhibit a non-monotonic behaviou, an unusual featue that highlights the limitations of the existing no-hai theoems. We also detemine and study in detail the scala chage, hoizon aea and entopy of ou solutions. 1 Email: anton296@umn.edu 2 Email: abakop@cc.uoi.g 3 Email: pkanti@cc.uoi.g

1 Intoduction The constuction of genealised gavitational theoies, with the inclusion of exta fields o highe-cuvatue tems in the action, has attacted an enomous inteest duing the last decades [1, 2]. The main eason is that these theoies may povide the famewok fo the ultimate theoy of gavity in the context of which seveal poblems of the taditional Geneal Relativity may be esolved. Theefoe, in the context of these modified gavitational theoies, diffeent aspects of gavity, fom black-hole solutions to cosmological solutions, have been e-addessed and, in seveal occasions, shown to lead to novel, inteesting solutions. Such a class of solutions was the one descibing egula black holes with a non-tivial scala field in the exteio egion, a type of solutions fobidden by Geneal Relativity. In an attempt to constuct a citeion whose fulfillment o violation could, espectively, fobid o allow the emegence of such solutions, no-hai theoems wee fomed. In its fist fom, the old no-hai theoem [3] excluded static black holes with a scala field; howeve, this was soon outdated by the discovey of black holes with Yang-Mills [4], Skyme fields [5] o with a confomal coupling to gavity [6]. These led to the fomulation of the novel no-hai theoem [7] that was ecently extended to cove the case of standad scala-tenso theoies [8]; a new fom was also poposed, that coves the case of Galileon fields [9]. Howeve, in a limited numbe of theoies, even the novel no-hai theoem may be evaded. The fist counte-example appeaed vey soon afte its fomulation, and demonstated the existence of egula black holes with a scala hai in the context of the Einstein-Dilaton-Gauss-Bonnet theoy [10] (fo some ealie studies that paved the way, see [11 15]). Additional solutions in the context of the same theoy wee deived in the pesence of a Yang-Mills field [16,17], in an abitay numbe of dimensions [18] o in the case of otation [19 22] (fo a numbe of inteesting eviews on the topic, see [23 25]). The moe ecent fom of the novel no-hai theoem [9] was also shown to be evaded [26] and concete solutions wee constucted [27, 28]. A common featue of the theoies that evaded the no-hai theoems was the pesence of highe-cuvatue tems, such as the quadatic Gauss-Bonnet (GB) tem inspied by the sting theoy [29] o Hondeski theoy [30]. It is the pesence of such tems that invalidate basic equiements of the no-hai theoems, and open the way fo the constuction of black-hole solutions with scala hai. In a ecent wok of ous [31], we consideed a geneal class of Einstein-scala-GB theoies, of which the cases [10, 26] constitute paticula examples. We demonstated that, unde cetain constaints on the fom of the coupling function between the scala field and the Gauss-Bonnet tem, and in conjunction with the pofile of the scala field itself, a egula black-hole hoizon egime and an asymptoticallyflat egime may be smoothly connected, and thus the no-hai theoems may be evaded. A numbe of novel black-hole solutions with scala hai was thus detemined and biefly pesented [31]. In the pesent wok, we povide additional suppot to the aguments pesented in [31]. We conside a geneal gavitational theoy containing the Ricci scala, a scala field and the GB tem, with the latte two quantities being coupled togethe though a coupling 1

function f(φ). Guided by the findings of ou pevious wok [31], we impose the afoementioned constaints on the coupling function f and the scala field φ, and investigate the existence of egula, black-hole solutions with a non-tivial scala hai. We find a lage numbe of such solutions fo a vaiety of foms fo the coupling function: exponential, polynomial (even and odd), invese polynomial (even and odd) and logaithmic. In all cases, the solutions fo the metic components, scala field, cuvatue invaiant quantities and components of the enegy-momentum tenso ae deived and discussed. Futhe chaacteistics of the poduced solutions, such as the scala chage, hoizon aea and entopy, ae also detemined, studied in detail and compaed to the coesponding Schwazschild values. 2 The Theoetical Famewok We stat with the following action functional that descibes a geneal class of highecuvatue gavitational theoies [31]: S = 1 d 4 x [ g R 1 ] 16π 2 µφ µ φ + f(φ)rgb 2. (1) In this, the Einstein-Hilbet tem given by the Ricci scala cuvatue R is accompanied by the quadatic Gauss-Bonnet (GB) tem RGB 2 defined as R 2 GB = R µνρσ R µνρσ 4R µν R µν + R 2, (2) in tems of the Riemann tenso R µνρσ, Ricci tenso R µν and the Ricci scala R. A scala field φ also appeas in the action (1), and couples to the GB tem though a geneal coupling function f(φ). This is a necessay equiement in ode fo the GB tem, that is a total deivative in fou dimensions, to contibute to the field equations. The gavitational field equations and the equation fo the scala field may be deived by vaying the action (1) with espect to the metic tenso g µν and the scala field φ, espectively. These have the fom: G µν = T µν, (3) 2 φ + f(φ)r 2 GB = 0, (4) whee G µν is the Einstein tenso and T µν the enegy-momentum tenso. The latte eceives contibutions fom both the scala field and the Gauss-Bonnet tem and is given by T µν = 1 4 g µν ρ φ ρ φ + 1 2 µφ ν φ 1 2 (g κλαβ ρµg λν + g λµ g ρν ) η Rργ αβ γ κ f(φ), (5) with R ργ αβ = ηργστ R σταβ = ɛργστ g R σταβ. (6) Thoughout this wok, the dot ove the coupling function denotes its deivative with espect to the scala field (i.e. f = df/dφ), and we employ units in which G = c = 1. 2

Ou aim is to find solutions of the set of Eqs. (3)-(4) that descibe egula, static, asymptotically-flat black-hole solutions with a non-tivial scala field. In paticula, we will assume that the line-element takes the following spheically-symmetic fom ds 2 = e A() dt 2 + e B() d 2 + 2 (dθ 2 + sin 2 θ dϕ 2 ), (7) and that the scala field is also static and spheically-symmetic, φ = φ(). In ou quest fo the afoementioned solutions, we will conside a vaiety of coupling functions f(φ), that will howeve need to obey cetain constaints [31]. By employing the line-element (7), we may obtain the explicit foms of Eqs. (3)-(4): the (tt), () and (θθ) components of Einstein s equations then, espectively, ead 4e B (e B + B 1) = φ 2[ 2 e B + 16 f(e B 1) ] 8 f [ B φ (e B 3) 2φ (e B 1) ], (8) 4e B (e B A 1) = φ 2 2 e B + 8 ( e B 3 ) fa φ, (9) e B[ A 2 2B + A (2 B ) + 2A ] = φ 2( e B 8 fa ) + while the equation fo the scala field takes the fom: 2φ + (4 + A B ) φ + 4 fe B +4 f [ A 2 φ + 2φ A + A (2φ 3B φ ) ], (10) [ (e B 3)A B (e B 1)(2A + A 2 ) ] = 0. (11) In the above, the pime in the metic functions A and B and the scala field φ denotes thei diffeentiation with espect to the adial coodinate. Befoe being able to numeically integate the system of equations (8)-(11), we need to educe the numbe of independent vaiables. We obseve that the ()-component (9) may take the fom of a second-ode polynomial with espect to e B, i.e. e 2B +βe B +γ = 0. Theefoe, this can be solved to give: whee β = 2 φ 2 e B = β ± β 2 4γ 2 (2 fφ + )A 1,, (12) γ = 6 fφ A. (13) 4 Accoding to the above, the solution fo the metic function B() may be easily found once the solutions fo the scala field φ() and the metic function A() ae detemined. In the emaining system of equations (8), (10) and (11), e B may theefoe be eliminated by using Eq. (12) while B may be eplaced by the expession B = γ + β e B 2e 2B + βe B, (14) 3

that follows by diffeentiating Eq. (12) with espect to the adial coodinate. Subsequently, the emaining thee equations (8), (10) and (11) fom a system of only two independent, odinay diffeential equations of second ode fo the functions A and φ: A = P S, (15) φ = Q S. (16) In the above equations, P, Q and S ae complicated expessions of (, e B, φ, A, f, f) that ae given in the Appendix. Note, that in these expessions we have eliminated, via Eq. (14), B that involves A and φ, but etained e B fo notational simplicity. 2.1 Asymptotic Solution at Black-Hole Hoizon We will stat ou quest fo black-hole solutions with a non tivial scala hai by detemining fist the asymptotic solutions of the set of Eqs. (8)-(11) nea the black-hole hoizon and at asymptotic infinity. These solutions will seve as bounday conditions fo ou numeical integation but will also povide impotant constaints on ou theoy (1). Nea the black-hole hoizon h, it is usually assumed that the metic functions and the scala field may be expanded as e A = e B = φ = a n ( h ) n, (17) n=1 b n ( h ) n, (18) n=1 n=0 φ (n) ( h ) n! ( h ) n, (19) whee (a n, b n ) ae constant coefficients and φ (n) ( h ) denotes the (nth)-deivative of the scala field evaluated at the black-hole hoizon. Equations (17)-(18) eflect the expected behaviou of the metic tenso nea the hoizon of a spheically-symmetic black hole with the solution being egula if the scala coefficients φ (n) ( h ) in Eq. (19) emain finite at the same egime. In a ecent wok of ous [31], we have followed instead the altenative appoach of assuming meely that, nea the hoizon, the metic function A() diveges, in accodance to Eq. (17). Then, the system of diffeential equations was evaluated in the limit h, and the finiteness of the quantity φ was demanded. This appoach was followed in [10] whee an exponential coupling function was assumed between the scala field and the GB tem. In [31], the fom of the coupling function f(φ) was left abitay, and the equiement of the finiteness of φ h was shown to be satisfied only unde the constaint hφ 3 h + 12f h + 2hφ 2 2 f h = 0, (20) h 4

whee all quantities have been evaluated at the hoizon h. The above is a second-ode polynomial with espect to φ h, and may be easily solved to yield the solutions: φ h = h 4f 1 ± 1 96 f h 2. (21) h h 4 The above ensues that an asymptotic black-hole solution with a egula scala field exists fo a geneal class of theoies of the fom (1). The only constaint on the fom of the coupling function aises fom the demand that the fist deivative of the scala field on the hoizon must be eal, which tanslates to the inequality f 2 h < 4 h 96. (22) Assuming the validity of the constaint (20), Eq. (15) then uniquely detemines the fom of the metic function A in the nea-hoizon egime; though Eq. (12), the metic function B is also detemined. Theefoe, the asymptotic solution of the Eqs. (12), (15) and (16), in the limit h, is given by the expessions e A = a 1 ( h ) +..., (23) e B = b 1 ( h ) +..., (24) φ = φ h + φ h( h ) + φ h( h ) 2 +.... (25) and descibes, by constuction, a black-hole hoizon with a egula scala field povided that φ obeys the constaint (21) and the coupling function f satisfies Eq. (22). We note that the desied fom of the asymptotic solution was deived only fo the choice of the (+)-sign in Eq. (12) as the ( )-sign fails to lead to a black-hole solution [31]. The afoementioned egulaity of the nea-hoizon solution should be eflected to the components of the enegy-momentum tenso T µν as well as to the scala invaiant quantities of the theoy. The non-vanishing components of the enegy momentum tenso (5) ae: [ T t t = e 2B 4 2 [ T = e B φ 4 T θ θ =T ϕ ϕ = e 2B 4 φ ( 2 2 e B + 16 f(e ) B 1) ] ( φ 8e B e B 3 ) fa 2 [φ 2 (e B 8 fa ) 8 f ( B φ (e B 3) 2φ (e B 1) )], (26), (27) 4 f ( A 2 φ + 2φ A + A (2φ 3B φ ) )]. (28) Employing the asymptotic behaviou given in Eqs. (23)-(25), we eadily deive the following appoximate behaviou: T t t = + 2e B B φ f + O( 2 h ), (29) T = 2e B A φ f + O( 2 h ), (30) T θ θ =T ϕ ϕ = e 2B φ f ( 2A + A 2 3A B ) + O( h ). (31) 5

The above expessions, in the limit h, lead to constant values fo all components of the enegy-momentum tenso. Similaly, one may see that the scala invaiant quantities R, R µν R µν and R µνρσ R µνρσ, whose exact expessions ae listed in the Appendix, educe to the appoximate foms R = + 2e B 2 ( e B 2A ) + O( h ), (32) R µν R µν = + 2e 2B 4 ( e B A ) 2 + O( h ), (33) R µνρσ R µνρσ = + 4e 2B 4 ( e 2B + 2 A 2) + O( h ). (34) In the above, we have used that, nea the hoizon, A B 1/( h ) and A 2 A, as dictated by Eqs. (23)-(24). Again, the dominant tem in each cuvatue invaiant adopts a constant, finite value in the limit h. Subsequently, the GB tem also comes out to be finite, in the same limit, and given by 2.2 Asymptotic Solution at Infinity R 2 GB = + 12e 2B 4 A 2 + O( h ). (35) At the othe asymptotic egime, that of adial infinity ( ), the metic functions and the scala field may be again expanded in powe seies, this time in tems of 1/. Demanding that the metic components educe to those of the asymptotically-flat Minkowski space-time while the scala field assumes a constant value, we wite: e A = 1 + e B = 1 + φ = φ + p n, (36) q n, (37) n=1 n=1 n=1 d n. (38) The abitay coefficients (p n, q n, d n ) ae in pinciple detemined upon substitution of Eqs. (36)-(38) in the field equations of the theoy. Howeve, two of the coefficients, namely p 1 and d 1, emain as fee paametes and ae associated with the Anowitt- Dese-Misne (ADM) mass and scala chage, espectively: p 1 2M and d 1 = D. The emaining coefficients may be calculated to an abitay ode: we have pefomed this 6

calculation up to ode O(1/ 6 ), and deived the following expessions: e A = 1 2M + MD2 12 3 96M 3 D 3MD 3 + 512M 2 + 24MD f + M 2 D 2 6 4 f 64D 2 f + 128MD 2 f + O(1/ 6 ), (39) 90 5 + 768M 4 208M 2D2 384MD f + 3D 4 48 4 e B = 1 + 2M + 16M 2 D 2 + 32M 3 5MD 2 4 2 4 3 + 6144M 5 2464M 3 D 2 + 97MD 4 6144M 2 Df + 192D 3 f 384MD 2 f 192 5 + O(1/ 6 ), φ = φ + D + MD + 32M 2 D D 3 2 24 + 12M 3D 24M 2 f MD 3 3 6 4 + 6144M 4 D 928M 2 D 3 + 9D 5 12288M 3 f 1536MD 2 f + O(1/ 6 ). (41) 1920 5 We obseve that the scala chage D modifies significantly the expansion of the metic functions at ode O(1/ 2 ) and highe. The existence itself of D, and thus of a non-tivial fom fo the scala field, is caused by the pesence of the GB tem in the theoy. The exact fom, howeve, of the coupling function does not ente in the above expansions ealie than the ode O(1/ 4 ). This shows that an asymptotically flat solution of Eqs. (8)-(11), with a constant scala field does not equie a specific coupling function and, in fact, aises fo an abitay fom of this function. The asymptotic solution at infinity, given by Eqs. (36)-(38), is also chaacteised by egula components of T µν and cuvatue invaiants. Employing the facts that, as, (e A, e B, φ) O(1) while (A, B, φ ) O (1/ 2 ), we find fo the components of the enegy-momentum tenso the asymptotic behaviou: T t t T T θ θ 1 ( ) 1 4 φ 2 + O. (42) Clealy, all of the above components go to zeo, as expected. A simila behaviou is exhibited by all cuvatue invaiants and the GB tem, in accodance to the asymptotically-flat limit deived above. In paticula, fo the GB tem, we obtain 6 (40) RGB 2 48M 2. (43) 6 2.3 Connecting the two asymptotic solutions In the pevious two subsections, we have constucted a nea-hoizon solution with a egula scala field and an asymptotically-flat solution with a constant scala field that was achieved unde mild constaints on the fom of the coupling function f(φ). Howeve, given the complexity of the equations of the theoy, it is the numeical integation of the system (15)-(16) that will eveal whethe these two asymptotic solutions may be 7

smoothly matched to ceate a black-hole solution with scala hai valid ove the entie adial domain. In fact, theoetical aguments developed decades ago, the so-called no-hai theoems, excluded in the past the emegence of such solutions in a vaiety of scala-tenso theoies. The olde vesion of the no-hai theoem [3] was applied in theoies with minimallycoupled scala fields: it employs the scala equation of motion and elies on the sign of V (φ), whee V (φ) is the potential of the scala field. In most theoies studied, the quantity V (φ) had the opposite sign fom the demanded one, and this excluded the emegence of the desied black-hole solutions. In [31], the same agument was applied in the case of the theoy (1), appopiately alteed to yield a constaint on the effective potential V eff (φ) = f(φ)rgb 2 of the scala field. Although this constaint was of an integal fom ove the entie adial egime, in special cases it meely demanded that f(φ)rgb 2 > 0. In Eqs. (35) and (43), the asymptotic values of the GB tem nea the hoizon and at infinity wee deived: they ae both positive, with the latte deceasing as inceases. These esults point to a monotonic deceasing behaviou of the GB tem fom an initial positive value nea the hoizon to a vanishing value at adial infinity - as we will shotly demonstate in the next section, this is indeed the behaviou of the GB tem. In that case, the only equiement fo the evasion of the old no-hai theoem is appaently the positivity of the coupling function f(φ). As the old no-hai theoem imposes in geneal mild constaints on a theoy, in [31] we consideed in addition the novel no-hai theoem [7] that applies also in theoies with a confomal coupling of the scala fields to gavity (see, also, [8]). This agument elies on the pofile of the T component of the enegy-momentum tenso of the theoy in tems of the adial coodinate. In a lage class of theoies, unde the assumptions of positivity and consevation of enegy, it is in geneal extemely difficult to smoothly match the neahoizon and asymptotic-infinity values of T. That pevents the emegence of black-hole solutions with scala hai. Howeve, the novel no-hai theoem was shown to be evaded in the context of the Einstein-Scala-Gauss-Bonnet theoy with an exponential coupling function [10] o a linea coupling function [26]. In both theoies, the pesence of the GB tem played a catalytic ole fo the emegence of the solutions. Theefoe, in [31], we kept again the fom of the coupling function f(φ) abitay, and econsideed the agument of the novel no-hai theoem. Ou analysis evealed that the evasion of this theoem holds fo a geneal class of theoies involving the GB tem: the pofile of T in tems of may be easily made smooth and monotonic unde the assumptions that, nea the hoizon, fφ < 0, fφ + fφ 2 > 0. (44) The fist constaint, accoding to Eq. (30) ensues the positivity of T in the nea-hoizon egime; the second that (T ) is negative in the same egime. Then, in conjunction with the behaviou descibed by Eq. (42) at adial infinity, T is positive and deceasing ove the whole adial egime. This behaviou invalidates the equiements set by the novel no-hai theoem and thus causes its evasion. In fact, the fist of the constaints listed in Eq. (44) is aleady satisfied: Eq. (21) dictates that the combination fφ at the hoizon is always negative to ensue the egulaity of the black-hole hoizon. Theefoe, if f > 0, then φ h must be necessaily negative, o 8

vice vesa. The emaining constaint fφ + fφ 2 > 0 may be altenatively witten as ( fφ ) h > 0; this meely demands that the afoementioned negative value of the quantity ( fφ ) h should be constained away fom the hoizon so that the two asymptotic solutions (23)-(25) and (36)-(38) can smoothly match. As we will see in the next section, this second constaint is automatically satisfied fo all the solutions found and does not demand any fine-tuning of ou paametes. 3 Numeical Solutions The deivation of exact solutions, valid ove the entie adial domain, demands the numeical integation of the system (15)-(16). Ou integation stats at a distance vey close to the hoizon of the black hole, i.e. at h + O(10 5 ) (fo simplicity, we set h = 1). Thee, we use as bounday conditions the asymptotic solution (23)-(25) togethe with Eq. (21) fo φ h upon choosing a paticula coupling function f(φ). The integation poceeds towads lage values of the adial coodinate until the fom of the deived solution matches the asymptotic solution (36)-(38). In the next sub-sections, we pesent a vaiety of egula black-hole solutions with scala hai fo diffeent choices of the coupling function f(φ). 3.1 Exponential Coupling Function Fist, we conside the case whee f(φ) = αe κφ. Accoding to the aguments pesented in Section 2.3, the evasion of the old no-hai theoem is ensued fo f(φ) > 0; theefoe, we focus on the case with α > 0. As the exponential function is always positive-definite, the sign of f = ακe κφ is then detemined by the sign of κ. In ode to evade also the novel no-hai theoem and allow fo egula black-hole solutions to emege, we should satisfy the constaint fφ < 0, o equivalently κ φ < 0, nea the hoizon. Theefoe, fo κ > 0, we should have φ h < 0, which causes the deceasing of the scala field as we move away ϕ() 1.6 1.4 1.2 1.0 α=0.1, a 1 =1 ϕ h =1.5 ϕ h =1 ϕ h =0.5 0.10 0.05 0.00 0.00002 0.00000-0.00002-0.00004-0.00006 2.0 2.2 2.4 2.6 2.8 3.0 α=0.1, a 1 =1, ϕ h =0.5 0.8-0.05 T 0.6-0.10 T t t T θ θ 1 5 10 50 100 500 1000 1.0 1.5 2.0 2.5 3.0 Figue 1: The scala field φ (left plot), and the enegy-momentum tenso T µν (ight plot) in tems of the adial coodinate, fo f(φ) = αe φ. 9

6 5 4 α=0.01 g tt g 12 10 8 α=0.01 R 2 GB 3 R 2 GB 6 2 4 1 2 0 1 5 10 50 100 500 0 1.0 1.5 2.0 2.5 3.0 Figue 2: The metic components g tt and g (left plot), and the Gauss-Bonnet tem R 2 GB (ight plot) in tems of the adial coodinate, fo f(φ) = αe φ. fom the black-hole hoizon. The situation is evesed fo κ < 0 when φ h > 0 and the scala field inceases with. The case of f(φ) = αe φ, with α > 0, was studied in [10] and led to the well-known family of Dilatonic Black Holes. The solutions wee indeed egula and asymptoticallyflat with the scala field deceasing away fom the hoizon, in ageement with the above discussion. Hee, we pesent the complimentay case with f(φ) = αe φ (the exact value of κ does not alte the physical chaacteistics of the solution and, hee, we set it to κ = 1). At the left plot of Fig. 1, we pesent a family of solutions fo the scala field φ fo diffeent initial values φ h : fo κ = 1 < 0, the scala field must necessaily have φ h > 0, and theefoe inceases as inceases. The value of the coupling constant α, once the fom of the coupling function and the asymptotic value φ h ae chosen, is esticted by the inequality (22) - hee, we pesent solutions fo indicative allowed values of α. At the ight plot of Fig. 1, we pesent the enegy-momentum tenso components fo an indicative solution of this family: clealy, all components emain finite ove the whole adial egime. In paticula, T emains positive and monotonically deceases towads infinity, exactly the behaviou that ensues the evasion of the novel no-hai theoem. As discussed in the pevious section, apat fom choosing the input value φ h fo ou numeical integation in accodance with Eq. (21), no othe fine-tuning is necessay; the second constaint fo the evasion of the novel no-hai theoem, i.e. fφ + fφ 2 > 0, is automatically satisfied without any futhe action, and this is eflected in the deceasing behaviou of T component nea the hoizon. At the left and ight plots of Fig. 2, we also pesent the solution fo the two metic components ( g tt, g ) and the GB tem RGB 2, espectively. The metic components exhibit the expected behaviou nea the black-hole hoizon with g tt vanishing and g diveging at h = 1. In ode to ensue asymptotic flatness at adial infinity, the fee paamete a 1 appeaing in the nea-hoizon solution (23) is appopiately chosen. On the othe hand, the GB tem emains finite and positive-definite ove the entie adial domain - in fact it displays the monotonic behaviou, hinted by its two asymptotic limits 10

0.00 0.0 D/α 1/2-0.05-0.10-0.15-0.20-0.25 α=0.1, h=1 D/α 1/2-0.2-0.4-0.6 0.0002<α<0.113 φ h =0.1, h =1-0.30-0.8-0.35 0 2 4 6 8 10 φ h 2 5 10 20 M/α 1/2 Figue 3: The scala chage D as a function of the nea-hoizon value φ h (left plot), and of its mass M (ight plot), fo f(φ) = αe φ. (35) and (43), that causes the evasion of the old no-hai theoem. As expected, it contibutes significantly nea the hoizon, whee the cuvatue is lage, and quickly fades away as we move towads lage distances. The pofile of the metic components and GB tem exhibit the same qualitative behaviou in all families of black-hole solutions pesented in this wok, so we efain fom giving additional plots of these two quantities in the next sub-sections. The pofile of the scala chage D as a function of the nea-hoizon value φ h and of the mass M is given at the left and ight plot, espectively, of Fig. 3 (each dot in these, and subsequent, plots stand fo a diffeent black-hole solution). Fo the exponential coupling function f(φ) = αe φ, and fo fixed α and h, the scala field nea the hoizon φ h may ange fom a minimum value, dictated by Eq. (22), up to infinity. As the left plot eveals, as φ h, the coupling of the scala field to the GB tem vanishes, and we ecove the Schwazschild case with a tivial scala field and a vanishing chage. On the othe hand, in the ight plot, we obseve that as the mass of the black-hole inceases, the scala chage deceases in absolute value, and thus lage black holes tend to have smalle chages. It is also inteesting to study the pofiles of the aea of the black-hole hoizon, A h = 4πh 2, and of the entopy S h of this class of solutions. The entopy is defined though the elation [32] [ ] (βf ) S h = β F, (45) β whee F = I E /β is the Helmholtz fee-enegy of the system given in tems of the Euclidean vesion of the action I E. Also, β = 1/(k B T ) with the tempeatue following easily fom the definition [33, 34] T = k 2π = 1 4π ( 1 gtt g dg tt d ) h = a1 b 1 4π. (46) The calculation of the tempeatue and entopy of the dilatonic black hole, with an exponential coupling function of the fom f(φ) = αe φ, was pefomed in detail in [17]. 11

2500 1.00 2000 1500 0.94< h<4.71 α=0.1, φ h=0.1 A h/α S h/α 0.95 0.001<α<0.113 φ h=0.1, h=1 A h/a sch 1000 0.90 S h/s sch 500 0 2 3 4 5 6 7 0.85 2 3 4 5 6 M/α 1/2 M/α 1/2 Figue 4: The hoizon aea A h and entopy S h of the black hole (left plot) and thei atios to the coesponding Schwazschild values (ight plot) in tems of the mass M fo f(φ) = αe φ. By closely epeating the analysis, we find the following expessions fo the tempeatue and entopy T = 1 4π (2M + D) 2 h + 4f(φ h), (47) S h = A h 4 + 4πf(φ h) (48) of a GB black-hole aising in the context of ou theoy (1) with a geneal coupling function f(φ) between the scala field and the GB tem. We easily confim that, in the absence of the coupling function, the above quantities educe to the coesponding Schwazschild ones, T = 1/(4π h ) and S h = A h /4, espectively. By employing the expessions fo A h and S h as given above, we depict the hoizon aea and entopy, in tems of the mass of the black hole, at the left plot of Fig. 4: we obseve that both quantities incease fast as the mass inceases. The ight plot allows us to compae moe effectively ou solutions to the Schwazschild one. The lowe cuve depicts the atio of A h to the aea of the Schwazschild solution A Sch = 16πM 2, as a function of M; we obseve that, fo lage black-hole masses, the atio A h /A Sch appoaches unity, theefoe, lage GB black holes ae not expected to deviate in thei chaacteistics fom the Schwazschild solution of the same mass. On the othe hand, in the small-mass limit, the atio A h /A Sch significantly deviates fom unity; in addition, a lowe bound appeas fo the black-hole adius, and thus of the mass of the black hole, due to the constaint (22) not pesent in the Schwazschild case this featue has been noted befoe in the case of the dilatonic black holes [10,25]. It is woth noting that the GB tem, as an exta gavitational effect, causes the shinking of the size of the black hole as the atio A h /A Sch emains fo all solutions below unity. Tuning finally to the entopy S h of ou black-hole solutions, we find a simila patten: vey small GB black holes diffe themselves fom the Schwazschild solution by having a lowe entopy wheeas lage GB black holes tend to acquie, among othe chaacteistics, and the entopy of the Schwazschild solution. We obseve howeve that, apat fom the 12

vey small-mass egime close to the minimum value, the exponential GB black holes have in geneal a highe entopy than the Schwazschild solution, a chaacteistic that points to the themodynamical stability of these solutions. In fact, the dilatonic GB black holes [10], which compise a subclass of this family of solutions with a coupling function of the fom f(φ) = αe φ, have an identical entopy patten, and wee shown to be linealy stable unde small petubations moe than twenty yeas ago. 3.2 Even Polynomial Function Next, we conside the case whee f(φ) = αφ 2n, with n 1. Since the coupling function must be positive-definite, we assume again that α > 0. The fist constaint fo the evasion of the novel no-hai theoem, fφ < 0 nea the hoizon, now tanslates to φ h φ h < 0. Theefoe, two classes of solutions appea fo each value of n: one fo φ h > 0, whee φ h < 0 and the solution fo the scala field deceases with, and one fo φ h < 0, whee φ h > 0 and the scala field inceases away fom the black-hole hoizon. At the left plot of Fig. 5, we depict the fist family of solutions with φ h < 0 and φ h > 0 fo the choice f(φ) = αφ 2 while, at the left plot of Fig. 6, we depict the second class with φ h > 0 and φ h < 0 fo the choice f(φ) = αφ4. The complimentay classes of solutions may be easily deived in each case by evesing the signs of φ h and φ h. The fom of the enegy-momentum tenso components fo the two choices f(φ) = αφ 2 and f(φ) = αφ 4, and fo two indicative solutions, ae depicted at the ight plots of Figs. 5 and 6, espectively. We obseve that the qualitative behaviou of the thee components lagely emain the same despite the change in the fom of the coupling function f(φ) (note also the esemblance with the behaviou depicted at the ight plot of Fig. 1). In fact, the asymptotic behaviou of T µν nea the black-hole hoizon and adial infinity is fixed, accoding to Eqs. (29)-(31) and (42), espectively. Independently of the fom of the coupling function f(φ), at asymptotic infinity T appoaches zeo fom the positive side, while T t t and T θ θ do the same fom the negative side. In the nea-hoizon egime, -1 0.010 α=0.01, a 1 =1, ϕ h =-1 T ϕ() -2-3 -4-5 α=0.01, a 1 =1 ϕ h =-1 ϕ h =-2 ϕ h =-5 1 5 10 50 100 500 1000 0.005 0.000-0.005-0.010 T t t T θ θ 0.0001 0.0000-0.0001-0.0002-0.0003-0.0004 2.0 2.2 2.4 2.6 2.8 3.0 1.0 1.5 2.0 2.5 3.0 Figue 5: The scala field φ (left plot), and the enegy-momentum tenso T µν (ight plot) in tems of the adial coodinate, fo f(φ) = αφ 2. 13

ϕ() 1.0 0.9 0.8 0.7 0.6 α=0.01, a 1 =1 ϕ h =1 ϕ h =0.8 ϕ h =0.6 1 5 10 50 100 500 1000 0.0020 0.00002 0.0015 0.00000 0.0010-0.00002-0.00004 0.0005-0.00006 2.0 2.2 2.4 2.6 2.8 3.0 0.0000 α=0.01, a 1 =1, ϕ h =0.6-0.0005-0.0010-0.0015 T T t t T θ θ -0.0020 1.0 1.5 2.0 2.5 3.0 Figue 6: The scala field φ (left plot), and the enegy-momentum tenso T µν (ight plot) in tems of the adial coodinate, fo f(φ) = αφ 4. Eqs. (29)-(31) dictate that sign ( ) ( T t t, sign (T h ) h sign φ ) hf h, (49) sign ( ( T θ θ )h + sign φ ) hf h. (50) Using that f h φ h < 0 and the scaling behaviou of the metic functions nea the hoizon, we may easily deive that (T t t) h and (T ) h always assume positive values while (T θ θ ) h assumes a negative one. The fom of f(φ) meely changes the magnitude of these asymptotic values: in the case of a polynomial coupling function, the highe the degee is, the lage the asymptotic values nea the hoizon ae. One could assume that the intemediate behaviou of the scala field and the enegymomentum tenso always emains qualitatively the same. In fact, this is not so. Let us fix fo simplicity the values of h and φ h, and gadually incease the value of the coupling paamete α; this has the effect of inceasing the magnitude of the GB souce-tem appeaing in the equation of motion (4) fo φ. As an indicative case, at the left plot of Fig. 7, we depict the behaviou of φ fo fou lage values of α, and, at the ight plot, the ϕ() 0.010 0.008 0.006 0.004 0.002 0.000-0.002-0.004 α=2 α=3 α=4 α=5 1 2 5 10 20 0.10 0.05 α=3, a 1 =0.9, ϕ h =0.01 T T t t T θ θ 0.00 2. 10-9 -0.05 0-2. 10-9 -4. 10-9 -0.10 8.0 8.5 9.0 9.5 10.0 1.0 1.1 1.2 1.3 1.4 1.5 Figue 7: The scala field φ (left plot), and the enegy-momentum tenso T µν (ight plot) in tems of the adial coodinate, fo f(φ) = αφ 2 and vaious values of the coupling constant α. 14

D/α 1/2 0.05 0.00-0.05 0.0005<α<0.51 φ h =-0.1, h =1 0.995 0.990 0.985 0.005<α<0.51 φ h=-0.1, h=1 A h/a sch S h/s sch -0.10 0.980-0.15 1 2 5 10 20 M/α 1/2 1.0 1.5 2.0 2.5 3.0 M/α 1/2 Figue 8: The scala chage D (left plot), and the atios A h /A Sch and S h /S Sch (ight plot) in tems of the mass M, fo f(φ) = αφ 2. T µν fo one of these solutions. We obseve that all of these quantities ae not monotonic any moe; they go though a numbe of maxima o minima with that numbe inceasing with the value of α befoe eaching thei asymptotic values at infinity. Note that the nea-hoizon behaviou of both φ and T µν is still the one that guaantees the evasion of the no-hai theoem. We may thus conclude that the pesence of the GB tem in the theoy not only ensues that the asymptotic solutions (23)-(25) and (36)-(38) may be smoothly connected to ceate a egula black hole but it allows fo this to happen even in a non-monotonic way. Let us also study the chaacteistics of this class of black-hole solutions aising fo an even polynomial coupling function. In Fig. 8 (left plot), we depict the scala chage D in tems of the mass M of the black hole, fo the quadatic coupling function f(φ) = αφ 2 : we obseve that, in this case, the function D(M) is not monotonic in the small-mass egime but it tends again to zeo fo lage values of its mass. In tems of the nea-hoizon value φ h, the scala chage exhibits the expected behaviou: fo lage values of φ h, the effect of the GB tem becomes impotant and D inceases; on the othe hand, fo vanishing φ h, i.e. a vanishing coupling function, the scala chage also vanishes in ode to minimise the numbe of figues, we efain fom showing plots depicting the anticipated behaviou; in the same spiit, we pesent no new plots fo the quatic coupling function as it leads to exactly the same qualitative behaviou. Tuning to the hoizon aea A h and entopy S h of these black-hole solutions, we find that, in tems of the mass M, they both quickly incease, showing a pofile simila to that of Fig. 4 (left plot) fo the exponential case. The atio A h /A Sch emains again below unity ove the whole mass egime, and intepolates between a value coesponding to the lowest allowed value of the mass, accoding to Eq. (22), and the asymptotic Schwazschild value at the lage-mass limit. The entopy atio S h /S Sch, on the othe hand, is found to have a diffeent pofile by emaining now always below unity - this featue points pehaps towads a themodynamic instability of the even polynomial GB black holes compaed to the Schwazschild solution. 15

ϕ() 5 4 3 2 1 0 α=0.1, a 1 =0.6 ϕ h =5 ϕ h =2 ϕ h =1 1 5 10 50 100 500 1000 0.5 0.0-0.5-1.0 0.002 0.000-0.002-0.004-0.006-0.008-0.010 2.0 2.2 2.4 2.6 2.8 3.0 α=0.1, a 1 =0.6, ϕ h =5 T T t t T θ θ 1.0 1.5 2.0 2.5 3.0 Figue 9: The scala field φ (left plot), and the enegy-momentum tenso T µν (ight plot) in tems of the adial coodinate, fo f(φ) = αφ. 3.3 Odd Polynomial Function We now conside the case of f(φ) = α φ 2n+1, with n 0 and α > 0. Hee, the constaint fφ < 0 tanslates to φ 2n φ < 0, o simply to φ h < 0 fo all solutions. In Fig. 9, we have chosen the linea case, i.e. f(φ) = αφ, and pesented an indicative family of solutions fo the scala field (left plot) and the components of the enegy-momentum tenso fo one of them (ight plot). The deceasing pofile of φ fo all solutions, as we move away fom the black-hole hoizon, is evident and in ageement with the above constaint. The enegy-momentum tenso clealy satisfies the analytically pedicted behaviou at the two asymptotic egimes, that once again ensues the evasion of the novel no-hai theoem. The two plots in Fig. 10, left and ight, depict the same quantities but fo the case f(φ) = αφ 3. Thei pofile agees with that expected fo a egula, black-hole solution with a scala hai. The aleted eade may notice that, hee, we have chosen to pesent solutions with φ h < 0; fo an odd polynomial coupling function, these should have been pohibited unde the constaint f(φ) > 0, that follows fom the old no-hai theoem [3,31]. Nevetheless, egula black-hole solutions with a non-tivial scala field do emege, that do not seem to satisfy f(φ) > 0. A set of such solutions ae shown at the left plot ϕ() -0.8-1.0-1.2-1.4-1.6-1.8-2.0 α=0.01, a 1 =1 ϕ h =-0.7 ϕ h =-1 ϕ h =-1.50 1 5 10 50 100 500 1000 0.004 0.002 0.000-0.002-0.004 0.00005 0.00000-0.00005-0.00010-0.00015-0.00020-0.00025 α=0.01, a 1 =1, ϕ h =-0.7 T T t t T θ θ 1.0 1.5 2.0 2.5 3.0 2.0 2.2 2.4 2.6 2.8 3.0 Figue 10: The scala field φ (left plot), and the enegy-momentum tenso T µν plot) in tems of the adial coodinate, fo f(φ) = αφ 3. (ight 16

ϕ() 0.00-0.05-0.10 α=50 α=60 α=70-0.15 α=80-0.20 1 5 10 50 100 500 1000 0.015 0.010 0.005 0.000-0.005-0.010-0.015 2. 10-9 1. 10-9 0-1. 10-9 -2. 10-9 α=80, a 1 =1, ϕ h =0.01 T T t t T θ θ 1 2 3 4 70 80 90 100 Figue 11: The scala field φ (left plot), and the enegy-momentum tenso T µν (ight plot) in tems of the adial coodinate, fo f(φ) = αφ 3 and a vaiety of values of the coupling constant α. of Fig. 10 (we efain fom showing the set of solutions with φ h > 0 as these have simila chaacteistics). As we obseve, all of them obey the φ h < 0 constaint imposed by the evasion of the novel no-hai theoem, and lead to the expected behaviou of T µν ; the latte may be clealy seen at the ight plot of Fig. 10 whee such a pohibited solution is plotted. The behaviou of the metic components and the GB tem continue to be given by plots simila to the ones in Fig. 2. The emegence of solutions that violate the constaint f(φ) > 0 is in fact a geneal featue of ou analysis and not an isolated finding in the case of odd polynomial coupling function. Appaently, the pesence of the coupling of the scala field to the GB tem not only opens the way fo black-hole solutions to emege but endes the old no-hai theoem incapable of dictating when this may happen. Looking moe caefully at the agument on which the old no-hai theoem was based on [3,31], one eadily ealizes that this involves the integal of the scala equation ove the entie exteio egime, and thus the global solution of the field equations whose chaacteistics cannot be pedicted befoehand. In contast, the evasion of the novel no-hai theoem is based on local quantities, such as the enegy-momentum components and thei deivatives at paticula adial egimes, which may be easily computed. In addition, it is indiffeent to the behaviou of the solution in the intemediate egime, which may indeed exhibit an abitay pofile as the one pesented in the plots of Fig. 7. In fact, all solutions found in the context of ou analysis, with no exception, satisfy the constaints that ensue the evasion of the novel no-hai theoem. In this case, too, one may deive solutions fo a vaiety of values of the coupling constant α, as long as these obey the constaint (22). In Fig. 11 (left plot), we depict a family of solutions with φ h = 0.01 and a vaiety of values of α. This family of solutions pesent a less monotonic pofile compaed to the one exhibited by the solutions in Fig. 10. The components of the enegy-momentum tenso fo one of these solutions is depicted at the ight plot of Fig. 11, and they pesent a moe evolved pofile with the emegence of minima and maxima between the black-hole hoizon and adial infinity. We also notice that the solutions fo the scala field, although they stat fom the positive-value egime (φ h = 0.01), they coss to negative values fo faily small values of the adial coodinate. 17

1.2 1.0 0.0003<α<0.102 φ h=0.2, h=1 1.00 0.95 D/α 1/2 0.8 0.6 0.90 0.0003<α<0.102 φ h=0.2, h=1 A h/a sch 0.4 0.85 S h/s sch 0.2 0.0 2 5 10 20 0.80 2 3 4 5 6 M/α 1/2 M/α 1/2 Figue 12: The scala chage D (left plot), and the atios A h /A Sch and S h /S Sch (ight plot) in tems of the mass M, fo f(φ) = αφ. This behaviou causes the odd coupling function to change sign along the adial egime, a featue that makes any application of the old no-hai theoem even moe challenging. Tuning again to the chaacteistics of the odd polynomial black-hole solutions, in Fig. 12 (left plot), we depict the scala chage D in tems of the mass M of the black hole, fo the linea coupling function f(φ) = αφ: hee, the function D(M) is monotonic and appoaches a vanishing asymptotic value as M inceases. The scala chage D has no dependence on the initial scala-field value φ h since it is the fist deivative f (φ) that appeas in the scala equation (4) and, fo a linea function, this is meely a constant. The hoizon aea A h and entopy S h exhibit again an inceasing pofile in tems of M simila to that of Fig. 4 (left plot) fo the exponential case. The moe infomative atios A h /A Sch and S h /S Sch ae given in the ight plot of Fig. 12: as in the case of quadatic and quatic coupling functions, both quantities emain smalle than unity and intepolate between a lowest value coesponding to the black-hole solution with the lowest mass and the Schwazschild limit acquied at the lage mass limit. We now addess sepaately the chaacteistics of the black-hole solutions aising in the case of the cubic coupling function f(φ) = αφ 3 since hee we find a distinctly diffeent behaviou. As mentioned above, also in this case, as the coupling constant α inceases, fom zeo to its maximum value (fo h and φ h fixed), solutions with no monotonic pofile in tems of the adial coodinate aise (see left plot of Fig. 11). We depict the scala chage D in tems of the mass M of the black hole, fo the whole α-egime, in the left plot of Fig. 13: we may easily obseve the emegence of two diffeent banches of solutions (with a thid, shot one appeaing at the end of the uppe banch) coesponding to the same mass M. These banches appea at the small-mass limit of the solutions wheeas fo lage masses only one banch suvives with a vey small scala chage. In the ight plot of Fig. 13, we show the atio S h /S Sch in tems of the mass M: this quantity also displays the existence of thee banches with the one that is smoothly connected to the Schwazschild limit having the highe entopy. The two additional banches with the lage values of scala chage, compaed to the one of the Schwazschild banch, have a lowe entopy and they ae pobably less themodynamically stable. This behaviou was not obseved in the case of the quadatic coupling function whee moe evolved solutions fo the scala field also appeaed (see left plot of Fig. 5): thee, the function D(M) was not monotonic 18

1.0 0.15 D/α 1/2 0.10 3.402<α<340.207 φ h=0.01, h=1 0.8 0.6 0.4 3.402<α<340.207 φ h=0.01, h=1 A h/a sch S h/s sch 0.05 0.2 0.00 0.06 0.07 0.08 0.09 0.10 0.11 0.12 M/α 1/2 0.0 0.06 0.07 0.08 0.09 0.10 0.11 0.12 M/α 1/2 Figue 13: The scala chage D (left plot), and the atios A h /A Sch and S h /S Sch (ight plot) in tems of the mass M, fo f(φ) = αφ 3. but was always single-valued; that ceated shot, disconnected banches of solutions with slightly diffeent values of entopy atio S h /S Sch but all lying below unity. Let us finally note that the hoizon aea atio A h /A Sch, not shown hee fo bevity, has the same pofile as the one displayed in Fig. 12 fo the linea function while the D(φ h ) function shows the anticipated inceasing pofile as φ h, and thus the GB coupling, inceases. 3.4 Invese Polynomial Function The next case to conside is the one whee f(φ) = αφ k, whee k > 0, and α is also assumed to be positive, fo simplicity. Let us conside diectly some indicative cases: k = 1: In this case, the constaint fo the evasion of the novel no-hai theoem becomes: fφ = 2αφ /φ 2 < 0 which demands φ h > 0 fo all solutions. At the left plot of Fig. 14, we pesent a family of solutions fo the scala field emeging fo this coupling function. All solutions ae inceasing away fom the black-hole hoizon in accodance to the above comment. The components of the enegy-momentum ϕ() 1.20 1.15 1.10 1.05 1.00 α=0.01, a 1 =1 ϕ h =+1.00 ϕ h =+1.05 ϕ h =+1.10 ϕ h =+1.15 5 10 50 100 500 0.003 0.002 0.001 0.000-0.001-0.002 α=0.01, a 1 =1, φh=+1.00 T T t t T θ θ -1.0-0.003 1.0 1.5 2.0 2.5 3.0 1.0 1.5 2.0 2.5 3.0 10-4 1.0 0.5 0.0-0.5 Figue 14: The scala field φ (left plot), and the enegy-momentum tenso T µν plot) in tems of the adial coodinate, fo f(φ) = α/φ. (ight 19

1.20 0.010 α=0.01, a 1 =1, φh=+1.00 T ϕ() 1.15 1.10 1.05 1.00 α=0.01, a 1 =1 ϕ h =+1.00 ϕ h =+1.05 ϕ h =+1.10 ϕ h =+1.15 5 10 50 100 500 0.005 0.000-0.005-0.010 10-4 T t t T θ θ 1.0 1.5 2.0 2.5 3.0 1.0 0.5 0.0-0.5-1.0 1.0 1.5 2.0 2.5 3.0 Figue 15: The scala field φ (left plot), and the enegy-momentum tenso T µν plot) in tems of the adial coodinate, fo f(φ) = α/φ 2. (ight tenso ae also well behaved, as may be seen fom the ight plot of Fig. 14. As in the case of the odd polynomial function, an additional set of solutions aises with φ h < 0 with simila chaacteistics. k = 2: In this case, the constaint becomes: fφ = αφ /φ 3 < 0 which demands φ h φ h > 0. A family of solutions fo the scala field emeging fo this coupling function, with φ h > 0 and inceasing with, ae pesented at the left plot of Fig. 15 a complementay family of solutions with φ h < 0 and deceasing away fom the black-hole hoizon wee also found. The components of the enegy-momentum tenso fo an indicative solution ae depicted at the ight plot of Fig. 15, and clealy emain finite ove the whole exteio egime. Concening the chaacteistics of the invese polynomial GB black holes, we find again an inteesting behaviou in Fig. 16, we depict the invese linea case, as an indicative one. The scala chage D exhibits a monotonic deceasing behaviou in tems of the mass M, as one may see in the left plot of the figue, appoaching zeo at the lagemass limit. In tems of the input paamete φ h, the scala chage pesents the anticipated 0.00 1.10 D/α 1/2-0.05-0.10-0.15-0.20-0.25-0.30-0.35 0.0009<α<0.919 φ h =3, h =1 1.05 1.00 0.95 0.90 0.85 0.009<α<0.919 φ h=3, h=1 A h/a sch S h/s sch 0.5 1 5 10 M/α 1/2 1.0 1.5 2.0 M/α 1/2 Figue 16: The scala chage D (left plot), and the atios A h /A Sch and S h /S Sch (ight plot) in tems of the mass M, fo f(φ) = α/φ. 20

behaviou (and thus is not shown hee): fo an invese coupling function, D inceases as the value of φ h, and thus of the GB coupling, deceases. The quantities A h and S h incease once again quickly with the mass M, as in Fig. 4. The atio A h /A Sch is shown in the ight plot of Fig. 16, and eveals again the constantly smalle size of the GB black holes compaed to the asymptotic Schwazschild solution, as well as the existence of a lowest-mass solution. The atio S h /S Sch, depicted also in Fig. 16, eveals that the entie class of these black-hole solutions independently of thei mass have a highe entopy compaed to the asymptotic Schwazschild solution. This featue is in fact unique fo the invese linea function; a simila analysis fo the invese quadatic coupling function has poduced simila esults fo the quantities D(M), D(φ h ), A h /A Sch and S h /S Sch with the only diffeence being that the vey-low mass egime of the invese quadatic GB black holes have a lowe entopy than the Schwazschild solution, i.e. the situation esembles moe the one depicted in Fig. 4. 3.5 Logaithmic Function We finally addess the case whee f(φ) = α ln(φ). Hee, black-hole solutions emege fo fφ = αφ /φ < 0, nea the black-hole hoizon; fo α > 0, this tanslates to φ h < 0 (since the agument of the logaithm must be a positive numbe, i.e. φ > 0). As a esult, the solutions fo the scala field ae esticted to have a deceasing behaviou as we move away fom the black-hole hoizon - this is indeed the behaviou obseved in the class of solutions depicted at the left plot of Fig. 17. One may also obseve that the plot includes solutions with eithe φ h > 1 o φ h < 1, o equivalently with f > 0 o f < 0. Once again, the old no-hai theoem is poven to be inadequate to exclude the pesence of egula black holes with scala hai even in subclasses of the theoy (1). In contast, the deived solutions continue to satisfy the constaints fo the evasion of the novel no-hai theoem. The components of the enegy-momentum tenso ae pesented in the ight plot of Fig. 17: they exhibit the same chaacteistics as in the cases pesented in the pevious subsections with the most impotant being the monotonic, deceasing pofile of the T ϕ() 1.10 1.05 1.00 0.95 0.90 0.85 0.80 0.75 α=0.01, a 1 =1 ϕ h =+0.8 ϕ h =+0.9 ϕ h =+1.0 ϕ h =+1.1 5 10 50 100 500 0.002 0.001 0.000-0.001-0.002 10-4 α=0.01, a 1 =1, ϕ h =+1.00 T T t t T θ θ 1.0 1.5 2.0 2.5 3.0 1.0 0.5 0.0-0.5-1.0 1.0 1.5 2.0 2.5 3.0 Figue 17: The scala field φ (left plot), and the enegy-momentum tenso T µν plot) in tems of the adial coodinate, fo f(φ) = α ln φ. (ight 21

1.0 1.00 D/α 1/2 0.8 0.6 0.4 0.0002<α<0.204 φ h =2, h =1 0.95 0.90 0.85 0.002<α<0.204 φ h=2, h=1 A h/a sch S h/s sch 0.2 0.80 0.0 2 5 10 20 M/α 1/2 0.75 1.5 2.0 2.5 3.0 3.5 4.0 4.5 M/α 1/2 Figue 18: The scala chage D (left plot), and the atios A h /A Sch and S h /S Sch (ight plot) in tems of the mass M, fo f(φ) = α ln φ. component. The coupling constant α can also take a vaiety of values as long as it satisfies Eq. (22); in this case, the monotonic behaviou of φ ove the whole exteio space of the black hole is peseved independently of the value of α. The enegy-momentum tenso also assumes the same fom as in Fig. 17, and thus we efain fom pesenting any new plots. The scala chage D of the logaithmic GB black holes, in tems of the mass M is shown in the left plot of Fig. 18. Once again, we obseve that, as the mass of the blackhole solution inceases, D deceases towads a vanishing value ou pevious analysis has shown that this featue is often connected with the themodynamical stability of the solutions. Indeed, as it may be seen fom the ight plot of Fig. 18, the atio S h /S Sch is above unity fo a vey lage pat of the mass-egime; it is again the small-mass egime that is excluded fom the themodynamical stable solutions. The scala chage D has again the anticipated behaviou in tems of the paamete φ h : since it is the f (φ) = α/φ that entes the scala equation (4), the scala chage inceases as φ h deceases. The atio A h /A Sch is also shown in the ight plot of Fig. 18: as in the pevious cases, it coesponds to a class of smalle black-hole solutions, compaed to the Schwazschild solution of the same mass, with a minimum-allowed mass that appoaches the asymptotic Schwazschild solution as M eaches lage values. 4 Conclusions The emegence of egula black-hole solutions in theoies with scala fields has always attacted the inteest of eseaches especially in conjunction with the no-hai theoems. The latte wee fomed in an attempt to limit the emegence of black-hole solutions in the context of theoies beyond Geneal Relativity. Howeve, a small numbe of theoies, in which the evasion of those theoems was ealised, wee found in the liteatue ove the yeas leading to novel solutions with scala hai. A main chaacteistic of those theoies was the coupling of the scala field with higheode gavitational tems. When one focuses on theoies whee the highe-cuvatue tem 22