Next Article in Journal
Optical Singly-Ionized Iron Emission in Radio-Quiet and Relativistically Jetted Active Galactic Nuclei
Previous Article in Journal
On the Dominant Lunisolar Perturbations for Long-Term Eccentricity Variation: The Case of Molniya Satellite Orbits
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Aspects of Gauss-Bonnet Scalarisation of Charged Black Holes

Departamento de Matemática da Universidade de Aveiro and Centre for Research and Development in Mathematics and Applications (CIDMA), Campus de Santiago, 3810-183 Aveiro, Portugal
*
Author to whom correspondence should be addressed.
Submission received: 18 November 2021 / Revised: 30 November 2021 / Accepted: 2 December 2021 / Published: 7 December 2021
(This article belongs to the Special Issue Gravitational Waves and Compact Objects: From Theory to Data Analysis)

Abstract

:
The general relativity vacuum black holes (BHs) can be scalarised in models where a scalar field non-minimally couples to the Gauss-Bonnet (GB) invariant. Such GB scalarisation comes in two flavours, depending on the GB sign that triggers the phenomenon. Hereafter these two cases are termed GB± scalarisation. For vacuum BHs, only GB+ scalarisation is possible in the static case, while GB scalarisation is spin induced. But for electrovacuum BHs, GB is also charged induced. We discuss the GB scalarisation of Reissner-Nordström and Kerr-Newman BHs, discussing zero modes and constructing fully non-linear solutions. Some comparisons with GB+ scalarisation are given. To assess the generality of the observed features, we also briefly consider the GB± scalarisation of stringy dilatonic BHs and coloured BHs which provide qualitative differences with respect to the electrovacuum case, namely on the distribution and existence of regions triggering GB scalarisation.

Contents
1Introduction     1
2The Einstein-Maxwell-Scalar-GB Model      3
    2.1 GBε Scalarisation of Electrovacuum Solutions..................................................................................................................................................................... 4
    2.2 Physical Quantities of Interest for Scalarised BHs............................................................................................................................................................... 4
3GBε Scalarisation of Reissner-Nordström BHs      5
    3.1 The Linear Scalar Clouds........................................................................................................................................................................................................6
    3.2 The Non-Linear Spherically Symmetric Scalarised BHs.....................................................................................................................................................7
4GBε Scalarisation of Kerr-Newman BHs     9
    4.1 Construction of the Scalarised Kerr-Newman BHs.............................................................................................................................................................10
    4.2 Numerical Results...................................................................................................................................................................................................................12
5Lessons from Alternative Charged BHs     12
    5.1 Einstein-Maxwell-Dilaton BHs..............................................................................................................................................................................................13
   5.2 Einstein-Yang-Mills BHs.........................................................................................................................................................................................................13
6Further Remarks   14
References   16

1. Introduction

It is conceivable that deviations from General Relativity (GR) occur only for sufficiently large curvatures. One explicit realization of this idea is the phenomenon of spontaneous scalarisation. The original idea proposed the scalarisation of neutron stars in scalar-tensor models [1]. More recently, it gained a new guise in which the scalarisation of the GR vacuum BH solutions becomes possible, in the context of extended scalar-tensor models that include the Gauss-Bonnet (GB) quadratic curvature invariant R GB 2 , as first pointed out in [2,3,4]. In the following we shall dub the latter GB scalarisation.
GB scalarisation circumvents well-known no-hair theorems (see [5] for a review) due to a certain class of non-minimal couplings between a real scalar field ϕ and the GB invariant. The phenomenon occurs for BHs in an appropriate mass range, defined by a dimensionful coupling in the model. Moreover, it can be triggered either if R GB 2 > 0 —hereafter dubbed GB + scalarisation [6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21]—or if R GB 2 < 0 —hereafter dubbed GB scalarisation. For the Kerr family of GR, the latter only occurs for sufficiently fast spinning BHs [22,23,24], which justifies the terminology spin-induced scalarisation [22]. By contrast, in the case of GB + scalarisation, Kerr BHs can also scalarise, but rotation actually suppresses the effects of scalarisation [17,25].
Enlarging the model to include charged BHs, however, GB scalarization ceases to rely solely on rotation. This can already be illustrated in electrovacuum GR. The Kerr-Newman solution develops a negative GB invariant for either sufficiently large dimensionless angular momentum j or sufficiently large dimensionless charge q. Thus, sufficiently near extremality, Kerr-Newman BHs develop regions with R GB 2 < 0 Figure 1. One may expect that the boundary of the region with R GB 2 < 0 marks the onset of the solutions prone to GB scalarisation, as for the Kerr case [26]. We shall confirm this expectation below, explicitly constructing some of the GB scalarised Kerr-Newman solutions and comparing them with the corresponding GB + scalarised solutions.
The discussion of the previous paragraph shows that, within electrovacuum BHs, GB scalarisation can be spin-induced or charge-induced (or both). Let us remark, however, that such GB charge-induced scalarisation is different from the scalarisation of charged BHs introduced in [27], where the non-minimal coupling occurs between the scalar field and the Maxwell field, with the GB term (or any curvature corrections) being absent.
Additionally, the aforementioned observations on the GB scalarisation of Kerr-Newman BHs show this process occurs even for spherically symmetric, non-spinning Reissner-Nordström (RN) BHs [28]. Moreover, the negative GB invariant always occurs in the immediate vicinity of the horizon. One may ask whether these features are generic. Is any charged BH model prone to GB scalarisation sufficiently close to the maximal allowed charge? And is the R GB 2 < 0 region supporting the scalarisation always occuring in the immediate vicinity of the horizon? Interestingly, neither of these features is generic, as we shall illustrate by considering two alternative models of charged (spherical) BHs.
This paper is organized as follows. In Section 2 we detail the Einstein-Maxwell-scalar-GB model and its equations of motion. A general discussion on the tachyonic instability associated with the scalarisation is given, followed by the relevant physical quantities for describing scalarised BHs. In Section 3 we consider the GB ± scalarisation of the RN BH, first discussing the sign of the GB invariant for the electrovacuum RN BH and then constructing both the linear scalar clouds and non-linear scalarised BHs. In Section 4 we consider the GB ± scalarisation of the Kerr-Newman BH. After a brief discussion on the sign of the GB invariant, we discuss the construction of the solutions and provide a sample of numerical results, focusing on the GB scalarisation case. In Section 5 we briefly consider Einstein-Maxwell-dilaton and Einstein-non-Abelian BHs, which yield two valuable lessons concerning GB ± scalarisation of charged BHs. We conclude with a brief discussion and final remarks in Section 6.

2. The Einstein-Maxwell-Scalar-GB Model

We wish to consider the Einstein-Maxwell-scalar-GB (EMsGB) model, described by the following action
S EMsGB = d 4 x g 1 4 R 1 4 F μ ν F μ ν 1 2 μ ϕ μ ϕ + ϵ λ 2 4 f ( ϕ ) R GB 2 ,
where R is the Ricci scalar with respect to the spacetime metric g μ ν , R GB 2 is the GB invariant
R GB 2 R α β μ ν R α β μ ν 4 R α β R α β + 4 R 2 ,
R α β μ ν is the Riemann tensor, R α β is the Ricci tensor, F μ ν = μ A ν ν A μ is the Maxwell field strength tensor where A = A μ d x μ is the U ( 1 ) gauge potential, f ( ϕ ) is a coupling function of the real scalar field ϕ , λ is a constant of the theory with dimension of length and ϵ = ± 1 is chosen for GB ϵ scalarisation.
Varying the action (1) with respect to the metric tensor g μ ν , we obtain the Einstein equations,
R μ ν 1 2 g μ ν R = 2 T μ ν ( eff ) .
The effective energy-momentum tensor T μ ν ( eff ) has three distinct components:
T μ ν ( eff ) = T μ ν ( s ) + T μ ν ( M ) + T μ ν ( GB ) ,
consisting of the (pure) scalar and Maxwell parts, respectively,
T μ ν ( s ) = μ ϕ ν ϕ 1 2 g μ ν α ϕ α ϕ , T μ ν ( M ) = F μ α F ν α 1 4 F α β F α β ,
and a third contribution due to the scalar-GB term in (1)
T μ ν ( GB ) = 2 ϵ λ 2 P μ γ ν α α γ f ( ϕ ) ,
where
P α β μ ν = 1 4 ε α β ρ σ R ρ σ γ δ ε μ ν γ δ = R α β μ ν + g α ν R β μ g α μ R β ν + g β μ R α ν g β ν R α μ + 1 2 g α μ g β ν g α ν g β μ R ,
and ε α β ρ σ is the Levi-Civita tensor. The scalar field equation is
2 ϕ + ϵ λ 2 4 d f ( ϕ ) d ϕ R GB 2 = 0 ,
while the (source-free) Maxwell equations have the usual form
μ F μ ν = 0 .

2.1. GB ϵ Scalarisation of Electrovacuum Solutions

Spontaneous scalarisation manifests itself at the linear level as a tachyonic instability. Let us assume that ϕ = 0 solves (7), which will hold for a class of coupling functions. Then, the field equations reduce to those of electrovacuum GR, and the corresponding solutions provide solutions of the full model (1) as well. In particular, the Kerr-Newman geometry will be a solution of this model. Thus, for concreteness we shall refer to the scalarisation of the Kerr-Newman solution in the following; but a similar discussion would hold for any GR electrovacuum solution.
Next, we consider scalar perturbations of the Kerr-Newman solution within the full model (1). Assuming a small- ϕ expansion for the coupling function
f ( ϕ ) = f | ϕ = 0 + 1 2 d 2 f d ϕ 2 | ϕ = 0 ϕ 2 + O ( ϕ 3 ) ,
the linearised scalar field Equation (7) around the Kerr-Newman solution becomes
( μ eff 2 ) ϕ = 0 , where μ eff 2 = ϵ λ 2 4 d 2 f d ϕ 2 | ϕ = 0 R GB 2 ,
where □ and R GB 2 are computed for the scalar-free Kerr-Newman solution.
If μ eff 2 is not strictly positive, the scalar field possesses a (spacetime dependent) tachyonic mass. Wherever this tachyonic mass is supported, such region potentially supports a spacetime instability, which is precisely the GB scalarisation. To simplify the discussion, we assume without any loss of generality that d 2 f / d ϕ 2 | ϕ = 0 is strictly positive. Then the condition μ eff 2 < 0 is equivalent to
ϵ R GB 2 > 0 .
When this condition is obeyed in some region(s) outside the BH horizon, GB ϵ scalarisation is triggered.
If the Kerr-Newman BH reduces to a Schwarzschild BH of mass M, then
μ eff 2 = ϵ λ 2 4 d 2 f d ϕ 2 | ϕ = 0 48 M 2 r 6 ,
and only GB + scalarisation is possible.

2.2. Physical Quantities of Interest for Scalarised BHs

When the above instability is present there is also a different class of solutions for the model (1), besides the electrovacuum ones. These are the scalarised solutions. We are interested in the case of stationary BHs. Generically, these solutions posses three global charges: the mass M, the electric charge Q and the angular momentum J. For the solutions in this work, there is also a “scalar charge” Q s , which is not, however, associated with a conservation law. There are also a number of relevant horizon quantities: the Hawking temperature T H , the horizon area A H , the entropy S, and the horizon angular velocity Ω H . The BH entropy is the sum of two terms,
S = S E + S sGB , with S E = 1 4 A H , S sGB = ϵ 1 2 λ 2 H d 2 x h f ( ϕ ) R ( 2 ) ,
where R ( 2 ) is the Ricci scalar of the induced horizon metric h. The solutions satisfy the Smarr law
M = 2 Ω H J + 2 T H S + Φ Q + M s ,
where Φ is the electrostatic potential and M s is a contribution of the scalar field
M s = 1 2 d 3 g g ( a ϕ ) 2 .
Also, the solutions satisfy the first law of BH thermodynamics
d M = T H d S + Ω H d J + Φ d Q ,
in which there is no contribution from the scalar field.
In this work we shall focus on the quadratic coupling function,1
f ( ϕ ) = ϕ 2 2 ,
which is the simplest choice of f ( ϕ ) that guarantees that ϕ = 0 satisfies the scalar Equation (7).
It is useful to observe that the equations of the model are invariant under the transformation
r α r , λ α λ ,
with r the radial coordinate and α > 0 an arbitrary positive constant. Only quantities invariant under (18) (e.g., Q / M or Q / λ ) have a physical meaning. Following standard terminology, we define the following reduced quantities:
q Q M , a H A H 16 π M 2 , t H 8 π T H M , j J M 2 ,
which will be considered in what follows.

3. GB ϵ Scalarisation of Reissner-Nordström BHs

Let us start by considering the spinless limit of the Kerr-Newman family, the RN BH. The corresponding metric and gauge field can be written as (see e.g., [29])
d s 2 = N ( r ) σ ( r ) 2 d t 2 + d r 2 N ( r ) + r 2 ( d θ 2 + sin 2 θ d φ 2 ) , V = V ( r ) d t ,
where
σ ( r ) = 1 , N ( r ) = 1 2 M r + Q 2 r 2 , V ( r ) = Q r .
This BH possess an event horizon at
r = r h = M + M 2 Q 2 .
Thus, 0 q 1 and q = 1 for the extremal RN BH.
The GB invariant of the RN metric reads
R GB 2 = 8 r 8 6 M 2 r 2 12 Q 2 M r + 5 Q 4 .
For Q 0 this always becomes negative for some region with r > 0 . This region, however, is cloaked by the horizon unless the largest root of the quadratic equation in (the square brackets in) (23) exceeds r h . This condition is
q Q 1 + 1 6 > r h .
Using (22) one can easily show that is possible for
q > q c 0.957058 .
Thus, for q c < q 1 , a RN BH can undergo GB scalarisation.

3.1. The Linear Scalar Clouds

At the onset of the tachyonic instability, the linearised scalar field Equation (10) on the RN background allows solutions known as scalar clouds. These occur for a discrete set of RN solutions, each solution corresponding to a particular harmonic scalar field mode. To see this, we perform a harmonic decomposition of the scalar field,
ϕ = U ( r ) Y m ( θ , φ ) ,
where Y m are the real spherical harmonics and , m are the associated quantum numbers with the usual ranges, = 0 , 1 , and m . For a RN BH background (20) and (21), the linearised scalar Equation (10) becomes a radial equation
r 2 N ( r ) U ( r ) = ( + 1 ) U ( r ) ϵ λ 2 r 2 12 M 2 r 2 + 10 Q 4 r 4 24 M Q 2 r 3 ,
where a prime denotes a derivative w . r . t . the radial coordinate r.
This equation has the following asymptotic solutions: near the horizon
U ( r ) = u 0 + r h Q 2 r h 2 ( + 1 ) + ϵ λ 2 r h 2 3 6 Q 2 r h 2 + Q 4 r h 4 u 0 ( r r h ) + ,
where u 0 is an arbitrary nonzero constant which, in numerics, we set to 1; and near spatial infinity
U ( r ) = Q s r + 1 + ,
where Q s is the scalar charge for = 0 .
Solving (27) with the above asymptotic behaviours can be viewed as an eigenvalue problem. For a given , requiring the scalar field to smoothly interpolate between the asymptotics (28) and (29), crossing n times the r-axis, selects a discrete set of RN parameters. Thus, each scalar cloud is characterized by three ’quantum numbers’ ( n , , m ) , with the radial function being degenerate in terms of m as a consequence of the spherical symmetry of the RN background. In this work we shall report results on nodeless spherically symmetric fundamental solutions only2, i.e., with = m = n = 0 .
Having chosen the clouds quantum numbers, taking λ as fixed scale set in the action, and fixing the reduced charged q, the radial equation has a solution for a specific dimensionless ratio λ / M . For instance, for ϵ = + 1 , = m = n = 0 and q = 0 the selected value is λ / M 1.704 , corresponding to the initial point of the blue dotted curve in Figure 2 (left panel). This is the zero mode of the GB + instability of Schwarzschild. It selects a mass scale. Smaller masses (larger λ / M ) describe BHs unstable against scalarisation; larger masses (smaller λ / M ) correspond to stable BHs.
The variation of λ / M with increasing q can be interpreted as follows. GB scalarisation of Schwarzschild BHs may be attributed to a repulsive gravitational effect of the GB term, which only becomes dominant for sufficiently small BHs (in terms of λ ). Adding electric charge introduces two competing effects. On the one hand the electric charge provides a repulsive gravitational effect for RN BHs. This facilitates scalarisation, making it available for larger BHs (larger M, smaller λ / M ). On the other hand, the repulsive gravitational effect of the GB invariant, which is at the source of the scalarisation phenomenon, becomes supressed (and eventually the GB term even changes sign in some region), when increasing q. This supresses scalarisation, making it available only for smaller BHs (smaller M, larger λ / M ). The trend observed in the blue dotted curve in the main panel of Figure 2 (left panel) suggest that for small q the RN charge repulsion dominates and for large q the GB charge suppression becomes dominant. The inset, on the other hand, gives the behaviour for ϵ = 1 . In this case, as q increases, in the allowed (large q) interval, the GB behaviour dominates, and due to the opposite sign it provides an ever more significant repulsive contribution, thus faciliating GB scalarisation, which therefore is available for larger masses. We also remark that, for GB scalarisation, the ratio λ / M appears to diverge as q q c , while it stays finite as q 1 .
The profiles of typical scalar clouds are shown in Figure 2 (right panel). For ϵ = + 1 and moderate q we recover the picture found in the Schwarzschild case, Q = 0 : a monotonically decreasing profile starting with some finite value at the horizon (see the inset). This is also true for the ϵ = 1 case (red curve in main panel). But for q > q c and ϵ = + 1 , a new qualitative behaviour emerges: the maximal value of the scalar cloud can be attained outside the horizon3—see the blue curve in the main panel.

3.2. The Non-Linear Spherically Symmetric Scalarised BHs

The linear scalar clouds just discussed can be continued to the non-linear regime. Their backreaction originates scalarised BHs. We shall now discuss their construction for the case of the spherical, nodeless scalar clouds.
The ansatz to obtain the scalarised BH solutions is (20), together with a radial scalar field
ϕ ϕ ( r ) .
As in the usual electrovacuum case, the Maxwell Equation (8) yields a total derivative, leading to a first integral
V ( r ) = Q σ ( r ) r 2 .
This introduces the electric charge measured at infinity, Q. The scalar field satisfies the equation
ϕ + 2 r + N N + σ σ ϕ ϵ λ 2 r 2 N σ [ ( 3 5 N ) N σ + σ ( 1 N ) N N 2 ) + 2 ( 1 N ) N σ ϕ = 0 ,
while the equations for the metric functions N , σ are too involved and shall not be displayed here.
We are interested in BH solutions with an event horizon located at r = r h > 0 . The equations of the model are subject to the following boundary conditions.
N | r h = 0 , σ | r h = σ h , ϕ | r h = ϕ h , V | r h = 0 ; N | = 1 , σ | = 1 , ϕ | = 0 , V | = Φ ,
where σ h , ϕ h are constants fixed by numerics and Φ is the electrostatic potential at infinity. The horizon data fix the Hawking temperature, the horizon area and the entropy of the solutions,
T H = σ h N ( r h ) 4 π , A H = 4 π r h 2 , S = π r h 2 + ϵ λ 2 f ( ϕ h ) .
A local solution compatible with these asymptotics can be constructed both close to the horizon (as a power series in ( r r h ) ) and at infinity (as a power series in 1 / r ). For example, the first terms in the far field expression of the solutions read
N = 1 2 M r + Q 2 + Q s 2 r 2 + , ϕ ( r ) = Q s r + , V ( r ) = Φ Q r + , σ ( r ) = 1 Q s 2 2 r 2 + .
With the details just laid out we have numerically constructed the nonlinear continuation of the scalar clouds solving the full equations of the EMsGB model, for both signs of ϵ .
Technically, the construction of the scalarised BHs is a one parameter shooting problem in terms of the value of the scalar field at the horizon ϕ h . The input parameters are r h , Q and λ . Fixing the length scale λ , this leads to a two dimensional parameter space for the problem. The numerical results for several values of the ratio Q / λ are shown in Figure 3.
Figure 3 shows that for a given ratio Q / λ and both values of ϵ , one finds a continuum of solutions which bifurcate from the corresponding RN BH supporting a scalar cloud with these parameters. This line has a finite extent, ending in a critical configuration where the numerical process fails to converge. A general explanation for this behaviour can be traced back to the fact that the radicand of a square root in the horizon expansion of the scalar field vanishes as the critical set is approached. This is a generic feature of GB-scalar models. An exception here are the ϵ = + 1 solutions emerging from RN BHs with q > q c , in which case the critical configurations seem to possess a curvature singularity for some radius outside the event horizon (see [28] for a discussion).
From Figure 3 one may highlight two qualitatively different features when comparing ϵ = ± 1 . First, scalarisation reduces (increases) a H for ϵ = + 1 ( ϵ = 1 ). Secondly, “overcharged” solutions with q > 1 exist for ϵ = 1 only.

4. GB ϵ Scalarisation of Kerr-Newman BHs

Let us now address the GB ϵ scalarisation of the spinning, charged electrovacuum Kerr-Newman BH. This is a solution of the model (1), with the coupling (17), together with a vanishing scalar field, ϕ = 0 . This BH is described by its ADM mass M, total angular momentum per unit mass a = J / M and electric charge Q. In Boyer-Lindquist coordinates it reads (see e.g., [29])
d s 2 = Δ Σ d t a sin 2 θ d φ 2 + Σ d r 2 Δ + d θ 2 + sin 2 θ Σ a d t Σ + a 2 sin 2 θ d φ 2 ,
and
A μ d x μ = Q r Σ d t a sin 2 θ d φ ,
where
Δ r 2 2 M r + a 2 + Q 2 , Σ r 2 + a 2 cos 2 θ .
The event horizon of this solution is located at
r h = M 1 + 1 j 2 q 2 .
This implies the Kerr-Newman bound, j 2 + q 2 1 ; an extremal BH saturates this bound.
The Kerr-Newman metric has a GB invariant
R GB 2 = 48 M 2 Σ 3 1 2 a 2 Σ 3 3 r 2 a 2 cos 2 θ 2 cos 2 θ + 8 Q 2 Σ 6 r 4 ( 5 Q 2 12 M r ) + a 2 cos 2 θ 2 r 2 ( 19 Q 2 + 60 M r ) + 5 a 2 ( Q 2 12 M r ) cos 2 θ .
We have studied4 the sign of this quantity as a function of the parameters ( j , q ) Figure 1, observing that the qualitative picture found in the Kerr ( q = 0 ) or RN ( j = 0 ) cases is still valid for a Kerr-Newman BH. While R GB 2 is positive for large values of the radial coordinate, its sign close to the event horizon depends on the value of ( j , q ) . That is, for fixed j (or q), the GB invariant R GB 2 always becomes negative in a region outside the horizon, for large enough values of q (or j). In the presence of rotation, this region is located around the poles of the horizon, θ = 0 , π Figure 4. In Figure 1 we show the region in the ( j , q ) -domain where the GB invariant takes a negative sign at the poles of the horizon. Kerr-Newman BHs with R GB 2 < 0 have the potential to be scalarised for both signs of ϵ .

4.1. Construction of the Scalarised Kerr-Newman BHs

To construct the GB ϵ scalarised Kerr-Newman BHs, we shall use the ansatz in [23,25], supplemented with a nonzero gauge field. The metric ansatz reads:5
d s 2 = e 2 F 0 N d t 2 + e 2 F 1 d r 2 N + r 2 d θ 2 + e 2 F 2 r 2 sin 2 θ ( d φ W d t ) 2 , with N = 1 r h r ,
where the metric functions F i , W , as well as the scalar field ansatz ϕ = ϕ ( r , θ ) , depend on r , θ only. The ansatz for the U ( 1 ) potential is
A = A φ sin θ ( d φ W d t ) + V d t ,
with two potentials, a magnetic one A φ and an electric one V, both depending on r , θ .
Setting A φ = V = 0 results in the spontaneously scalarised Kerr BHs in [23,25], albeit with a different coupling function. The limit W = A φ = 0 results in the scalarised RN BHs [27,31,32] discussed above, albeit for a different radial coordinate.
The general problem is solved subject to the following boundary conditions. Asymptotic flatness requires
lim r F i = lim r W = 0 , and lim r ϕ = lim r A φ = 0 , lim r V = Φ .
Axial symmetry and regularity impose the following boundary conditions on the symmetry axis, i.e., at θ = 0 , π :
θ F i = θ W = θ ϕ = θ V = A φ = 0 .
Moreover, the absence of conical singularities implies also that F 1 = F 2 on the symmetry axis.
The event horizon is located at a constant value of r = r h > 0 . Only non-extremal solutions can be studied within the metric ansatz (41). We proceed by introducing a new radial coordinate x = r 2 r H 2 , which simplifies the boundary conditions at the horizon and also the numerical treatment of the problem. This results in the following boundary conditions at the horizon
x F i | r = r h = x ϕ | r = r h = 0 , W | r = r h = Ω H , x A φ | r = r h = V | r = r h = 0 ,
where the constant Ω H > 0 is the horizon angular velocity. Then, an approximate expansion of the solution compatible with these boundary conditions can easily be constructed.
Specializing some of the aforementioned physical quantities of interest for the ansatz in use, we obtain that the Hawking temperature and the event horizon area are determined by the following horizon data,
T H = 1 4 π r h e F 0 ( 0 ) ( θ ) F 1 ( 0 ) ( θ ) , A H = 2 π r h 2 0 π d θ sin θ e F 1 ( 0 ) ( θ ) + F 2 ( 0 ) ( θ ) ,
with the near horizon expansion F i = F i ( 0 ) ( θ ) + x 2 F i ( 2 ) ( θ ) + , and i = 0 , 1 , 2 .
The ADM mass M, the angular momentum J, the scalar ’charge’ Q s , together with the the magnetic dipole momentum q m , the electrostatic potential Φ and the electric charge Q are read off from the far field asymptotic of the metric and matter functions
g t t = 1 + 2 M r + , g φ t = 2 J r sin 2 θ + , ϕ = Q s r + , A φ = q m sin θ r + , V = Φ Q r + .
We remark that both the metric functions and the scalar field are invariant w . r . t . the transformation θ π θ .

4.2. Numerical Results

With the setup just described, we have employed a numerical approach similar to the one in [23,25]. The typical numerical errors for the solutions so obtained and reported below are of the order of 10 3 .
The scalarised Kerr-Newman solutions possess four independent charges: the three global charges shared with their electrovacuum counterparts, ( M , J , Q ) , plus the scalar charge Q s . In our approach, the input parameters are: the event horizon radius r h , the horizon angular velocity Ω H , the asymptotic value of the electrostatic potential Φ (or the electric charge Q), together with the coupling constant λ (which specifies the theory). Therefore, after fixing the scale λ , we are left with a three dimensional parameter space. A full scanning of such large parameter space of scalarised Kerr-Newman BHs is therefore a time consuming task. Here we focus on illustrative sets of solutions which already capture the generic behaviour. Moreover, although we have verified that spinning scalarised solutions exist for both signs of ϵ , we shall focus on the results for the more novel case of spin/charge scalarisation, ϵ = 1 . In practice we have scanned the parameter space by varying both r h and Q, for several different values of Ω H . Alternatively, we have varied both r h and Ω H , for fixed values of Q.
Our numerical results suggest that the solutions share most of the properties of the scalarised Kerr BHs discussed in [23,24]. In Figure 5 we display the reduced quantities ( q , j ) (left panel) and ( j , M / λ ) (right panel) for solutions with Ω H λ = 0.469 and illustrative values of the ratio Q / λ . These scalarised solutions emerge from a Kerr-Newman BH supporting a zero mode solution of the scalar equation—a scalar cloud—corresponding to the blue squares in Figure 5. The profiles of such clouds (with the maximal value of the scalar field normalised to one) are shown in Figure 6, for two different Kerr-Newman BHs. One can see that, while for small values of the reduced angular momentum j, the profile of the scalar field looks rather similar to that found in spherical case (see Figure 2 (right)), the large-j ones have a strong θ -dependence, the minimal value of ϕ being always found in the equatorial plane. The nonlinear continuation of the clouds (with constant Q / λ ) results in a sequence of solutions terminating at a critical configuration, as in the RN case reported before, corresponding to the red circles in Figure 5. The main trend observed in Figure 5 is that fixing Q and Ω H in units of λ , the scalarised BHs have more mass (and thus smaller q and λ / M ) and larger j.
Although our scanning of the full parameter space was limited, extrapolating the existing numerical data, we antecipate that the (three dimensional) domain of existence of ϵ = 1 spinning, charged scalarised BHs is bounded by four sets of solutions: (i) the existence surface, which corresponds to the set of Kerr-Newman solutions supporting scalar clouds; (ii) the set of critical solutions, which form again a two dimensional surface; (iii) the static configurations, J = 0 , which corresponds to the ϵ = 1 scalarised RN solutions discussed in the previous section; and (iv) the neutral configurations, Q = 0 , which were studied with the choice of the coupling function (17) in [24]. As for the Q = 0 case, the existence surface is universal for any expression of the coupling function allowing for scalarisation. Concerning the set (ii) (critical solutions) the numerical process fails to converge as it is approached, as in the static limit. The explanation for this behaviour can again be traced back to the fact that the radicand of a square root in the horizon expansion of the solutions vanishes as the critical set is approached. Note also that the sets (ii)–(iv) are not universal; they depend on the choice of the coupling function f ( ϕ ) .

5. Lessons from Alternative Charged BHs

It is interesting to test the generality of some of the results above concerning the interplay between the introduction of charge and the scalarisation phenomenon. For this purpose, in this section we shall be considering the simpler case of static, spherically symmetric, charged BHs in some alternative models, rather than electrovacuum.
As a first observation, we remark that for a finite mass, asymptotically flat solution, the GB invariant R GB 2 is strictly positive for large enough r. As with the RN BH, a sufficient condition for the occurrence of GB induced scalarization for b o t h signs of ϵ is that R GB 2 < 0 at the horizon. In general, however, the sign of the GB invariant at the horizon depends on the matter content. Indeed, for a generic spherically symmetric BH spacetime and using the metric ansatz (20), a straightforward computation leads to the simple relation
R GB 2 | r = r h = 12 r h 2 + 16 ρ ( H ) 2 16 r h 2 2 ρ ( H ) + p θ ( H ) ,
where r h is the horizon radius, ρ ( H ) = T t t ( r h ) and p θ ( H ) = T θ θ ( r h ) . One can easily see that, for a generic matter content, the above quantity has no definite sign.
One may then ask the following two questions: ( 1 ) is R GB 2 | r = r h < 0 close to the maximal charge for any charged BH model? ( 2 ) is R GB 2 | r = r h < 0 a necessary condition for GB scalarisation for any charged BH model?
We will now show, by concrete illustrations, that both these questions have a n e g a t i v e answer.

5.1. Einstein-Maxwell-Dilaton BHs

To answer question ( 1 ) above, we have investigated the sign of the GB invariant (together with its behaviour in the bulk) for the stringy generalisation of the RN BH—the Gibbons-Maeda-Garfinkle-Horowitz-Strominger (GMGHS) family of BHs [33,34]. In our context, these solutions are found for an action of the form
S EMdsGB = d 4 x g 1 4 R 1 2 μ ψ μ ψ 1 4 e 2 α ψ F μ ν F μ ν 1 2 μ ϕ μ ϕ + ϵ λ 2 4 f ( ϕ ) R GB 2 ,
describing an Einstein-Maxwell-dilaton-scalar-GB model, which includes an extra scalar field (the dilaton ψ ) with a non-minimally coupling with the Maxwell term, where α 0 is a constant of the theory. The GMGHS solution is found for ϕ = 0 , with f ( ϕ ) satisfying the condition (9). It is easy to prove that the behaviour found in the RN case ( α = 0 , ψ = 0 ) is recovered for small enough values of α . In that case, for large enough values of the electric charge R GB 2 becomes negative in a region between horizon and some maximal value of the radial coordinate. However, a direct computation shows that, for α > 0.9036 (for the conventions used in [34]), R GB 2 is strictly positive at the horizon and also in the bulk, irrespective of the value of the electric charge. Thus, as with Schwarzschild vacuum BHs, GB scalarisation of the GMGHS with large enough values of the dilaton coupling constant becomes impossible. This suggests that it would be interesting to check the status of GB scalarisation of the rotating counterpart of the GMGHS BH, the well known Kerr-Sen BH [35].

5.2. Einstein-Yang-Mills BHs

To answer question ( 2 ) above, we have investigated the sign of the GB invariant for the case of Einstein–Yang-Mills (EYM) BHs with SU(2) non-Abelian hair (nA) [36,37,38]. In our context, these solutions are found for an action of the form
S EYMsGB = d 4 x g 1 4 R 1 4 F μ ν ( a ) F μ ν ( a ) 1 2 μ ϕ μ ϕ + ϵ λ 2 4 f ( ϕ ) R GB 2 ,
with F μ ν ( a ) the nA field strength and a = 1 , 2 , 3 . These so called coloured BHs are asymptotically flat and possess a single global “charge”—the ADM mass, despite the presence of a local magnetic field (see [39,40] for reviews). Another striking difference with respect to their (magnetic) RN Abelian counterparts is the existence of a smooth solitonic limit [41], obtained as the horizon size shrinks to zero. At the same time, there is no upper bound on their horizon size. However, the large EYM BHs are essentially Schwarzschild solutions; the contribution of the YM fields to the total ADM mass becomes negligible, albeit these fields are still nontrivial.
Contact with question ( 2 ) above comes from observing that the GB invariant is always positive at the horizon for these solutions. However, R GB 2 may take negative values in a shell which does not touch the horizon, i.e., for some range of the radial coordinate r h < r 1 < r < r 2 —see Figure 7 (left panel). This feature occurs for small enough BHs: using the metric form (20) and the conventions in [39], we confirmed such shell is present for 0 < r h 0.71 , or, equivalently, 0 < a H 0.158 .
Given this qualitative difference with the RN case, one could ask whether GB ϵ scalarisation is still possible for both signs of ϵ . The answer is positive, and we have constructed the corresponding GB scalar clouds, i.e., solved Equation (10) for a large set of EYM BH backgrounds and both values of ϵ —see Figure 7 (right panel). As expected, ϵ = + 1 scalar clouds exist for all non-Abelian BHs. The value of the ratio λ / M 1.704 corresponding to a Schwarzschild BH is approach asymptotically, as a H 1 (i.e., large EYM BHs). Also, scalar clouds with ϵ = 1 exist for all BHs with R GB 2 < 0 in a shell outside the horizon. Nonlinear continuations of these scalar clouds should exist, but we did not attempt to construct them.
This example makes clear that GB ϵ scalarisation of BHs is not necessarily supported and triggered in the immediate vicinity of the event horizon.

6. Further Remarks

There is a well known analogy between the spinning vacuum Kerr BHs and the electrically charged static RN solutions.6 They possess many similar properties at the level of a thermodynamical description; in particular, both RN and Kerr BHs have an extremal limit, with a finite horizon size. In the context of this work, it is interesting to note that the GB invariant of a RN BH changes sign if q is large enough as the Kerr one changes sign for sufficiently large j. Thus, it is natural to conjecture that the qualitative picture found concerning the GB ϵ scalarisation of Kerr BHs [22,23,24,25] should be essentially recovered when replacing rotation by electric charge, with the existence of both ϵ = ± 1 scalarised solutions.
In this work we have confirmed that this conjecture is true, and construct the corresponding scalarised RN BHs. That is, we provide evidence for the following scenario: given an expression for the coupling function f ( ϕ ) , two classes of charged RN-scalarised solutions may exist for the same global charges. The first one has ϵ = 1 , and can be viewed as a generalization of the Q = 0 solutions in [2,3,4]. The second has ϵ = 1 , and in this case the condition μ eff 2 < 0 is supported by a large enough charge to mass ratio q > q c = 0.957058 , which implies R GB 2 < 0 for some region outside the horizon. We have also presented a preliminary investigation of the spinning generalisations of the above BHs, i.e., the scalarised Kerr-Newman BHs.
In the last part of this work, we have addressed the generality of these results. First, we established that the GB invariant of the stringy dilatonic generalisation of the RN BH becomes strictly positive for large enough values of the dilaton coupling constant. Also, we pointed out the possibility that the ϵ = 1 BH scalarisation may also appear in situations where R GB 2 is negative in a spherical s h e l l outside and disconnected from the horizon. This is the case of the coloured BHs in EYM theory.
In this paper, to simplify the picture, we have assumed the absence of a self-interaction term for the scalar field in the action (1). GB + scalarisation of spherical BHs including such self-interactions is discussed e.g., in [42,43]. Our results could be generalised to include such self-interactions.
The results in this paper emphasise that GB scalarisation is not only spin induced, but can also be induced by electromagnetic fields. Whereas electric charge is typically neglected for astrophysical BHs, the astrophysical environment surrounding astrophysical BHs has non-trivial electromagnetic fields, by virtue of the ionized matter in accretion disks, sourcing, for instance, magnetic fields. Thus, it would be interesting to consider in detail if and under which conditions such environments, if sufficiently “heavy”, trigger GB scalarisation.

Author Contributions

Investigation, C.A.R.H., A.M.P. and E.R.; Project administration, C.A.R.H.; Supervision, C.A.R.H. and E.R.; Writing—original draft, C.A.R.H.; Writing—review & editing, A.M.P. and E.R. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Center for Research and Development in Mathematics and Applications (CIDMA) through the Portuguese Foundation for Science and Technology (FCT—Fundacão para a Ciência e a Tecnologia), references UIDB/04106/2020 and UIDP/04106/2020, and by national funds (OE), through FCT, I.P., in the scope of the framework contract foreseen in the numbers 4, 5 and 6 of the article 23, of the Decree-Law 57/2016, of 29 August, changed by Law 57/2017, of 19 July. The authors acknowledge support from the projects PTDC/FIS-OUT/28407/2017, CERN/FIS-PAR/0027/2019 and PTDC/FIS-AST/3041/2020. A. Pombo is supported by the FCT grant PD/BD/142842/2018. This work has further been supported by the European Union’s Horizon 2020 research and innovation (RISE) programme H2020-MSCA-RISE-2017 Grant No. FunFiCO-777740.

Acknowledgments

The authors would like to acknowledge networking support by the COST Action CA16104. We would like to thank S. Hod for an useful comment.

Conflicts of Interest

The authors declare no conflict of interest.

Notes

1
For spherical symmetry, we have also explored the exponential coupling studied in [3,8], and observed that the behaviour is qualitatively similar.
2
Similar solutions are likely to exist for any other values of the quantum numbers, some preliminary results being found for the = 1 , m = n = 0 case. An investigation of Q = 0 , = 1 static solutions has been reported in [17].
3
This feature can be explained by studying Equation (27) [28].
4
An alternative expression for (40), in terms of P 1 ( 1 + 1 j 2 q 2 ) r / r h and P 2 j cos θ , is
R GB 2 = 48 M 4 1 ( P 1 2 + P 2 2 ) 3 1 2 ( P 1 2 + P 2 2 ) 3 P 2 2 ( 3 P 1 2 P 2 2 ) 2 + q 2 P 1 ( P 1 4 10 P 1 2 P 2 2 + 5 P 2 4 ) q 2 12 5 P 1 4 38 P 1 2 P 2 2 + 5 P 2 4 ,
a form which has been employed in our study.
5
The Kerr-Newman BH can also be written in this coordinate system. The corresponding expressions in the Kerr limit can be found in [30].
6
Since the electric-magnetic duality is still valid for the (Abelian) models in this work, the solutions possess a dual magnetic description.

References

  1. Damour, T.; Esposito-Farese, G. Nonperturbative strong field effects in tensor—Scalar theories of gravitation. Phys. Rev. Lett. 1993, 70, 2220–2223. [Google Scholar] [CrossRef]
  2. Silva, H.O.; Sakstein, J.; Gualtieri, L.; Sotiriou, T.P.; Berti, E. Spontaneous scalarization of black holes and compact stars from a Gauss-Bonnet coupling. Phys. Rev. Lett. 2018, 120, 131104. [Google Scholar] [CrossRef] [Green Version]
  3. Doneva, D.D.; Yazadjiev, S.S. New Gauss-Bonnet black holes with curvature induced scalarization in the extended scalar-tensor theories. Phys. Rev. Lett. 2018, 120, 131103. [Google Scholar] [CrossRef] [Green Version]
  4. Antoniou, G.; Bakopoulos, A.; Kanti, P. Evasion of No-Hair Theorems and Novel Black-Hole Solutions in Gauss-Bonnet Theories. Phys. Rev. Lett. 2018, 120, 131102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Herdeiro, C.A.R.; Radu, E. Asymptotically flat black holes with scalar hair: A review. Int. J. Mod. Phys. D 2015, 24, 1542014. [Google Scholar] [CrossRef] [Green Version]
  6. Blázquez-Salcedo, J.L.; Doneva, D.D.; Kunz, J.; Yazadjiev, S.S. Radial perturbations of the scalarized Einstein-Gauss-Bonnet black holes. Phys. Rev. D 2018, 98, 084011. [Google Scholar] [CrossRef] [Green Version]
  7. Myung, Y.S.; Zou, D.C. Instability of Reissner–Nordström black hole in Einstein-Maxwell-scalar theory. Eur. Phys. J. C 2019, 79, 273. [Google Scholar] [CrossRef]
  8. Doneva, D.D.; Kiorpelidi, S.; Nedkova, P.G.; Papantonopoulos, E.; Yazadjiev, S.S. Charged Gauss-Bonnet black holes with curvature induced scalarization in the extended scalar-tensor theories. Phys. Rev. D 2018, 98, 104056. [Google Scholar] [CrossRef] [Green Version]
  9. Witek, H.; Gualtieri, L.; Pani, P.; Sotiriou, T.P. Black holes and binary mergers in scalar Gauss-Bonnet gravity: Scalar field dynamics. Phys. Rev. D 2019, 99, 064035. [Google Scholar] [CrossRef] [Green Version]
  10. Brihaye, Y.; Herdeiro, C.; Radu, E. The scalarised Schwarzschild-NUT spacetime. Phys. Lett. B 2019, 788, 295–301. [Google Scholar] [CrossRef]
  11. Myung, Y.S.; Zou, D.C. Quasinormal modes of scalarized black holes in the Einstein–Maxwell–Scalar theory. Phys. Lett. B 2019, 790, 400–407. [Google Scholar] [CrossRef]
  12. Minamitsuji, M.; Ikeda, T. Scalarized black holes in the presence of the coupling to Gauss-Bonnet gravity. Phys. Rev. D 2019, 99, 044017. [Google Scholar] [CrossRef] [Green Version]
  13. Silva, H.O.; Macedo, C.F.B.; Sotiriou, T.P.; Gualtieri, L.; Sakstein, J.; Berti, E. Stability of scalarized black hole solutions in scalar-Gauss-Bonnet gravity. Phys. Rev. D 2019, 99, 064011. [Google Scholar] [CrossRef] [Green Version]
  14. Brihaye, Y.; Ducobu, L. Hairy black holes, boson stars and non-minimal coupling to curvature invariants. Phys. Lett. B 2019, 795, 135–143. [Google Scholar] [CrossRef]
  15. Herdeiro, C.A.R.; Radu, E. Black hole scalarization from the breakdown of scale invariance. Phys. Rev. D 2019, 99, 084039. [Google Scholar] [CrossRef] [Green Version]
  16. Hod, S. Gauss-Bonnet black holes supporting massive scalar field configurations: The large-mass regime. Eur. Phys. J. C 2019, 79, 966. [Google Scholar] [CrossRef]
  17. Collodel, L.G.; Kleihaus, B.; Kunz, J.; Berti, E. Spinning and excited black holes in Einstein-scalar-Gauss–Bonnet theory. Class. Quantum Gravity 2020, 37, 075018. [Google Scholar] [CrossRef] [Green Version]
  18. Antoniou, G.; Lehébel, A.; Ventagli, G.; Sotiriou, T.P. Black hole scalarization with Gauss-Bonnet and Ricci scalar couplings. Phys. Rev. D 2021, 104, 044002. [Google Scholar] [CrossRef]
  19. East, W.E.; Ripley, J.L. Dynamics of Spontaneous Black Hole Scalarization and Mergers in Einstein-Scalar-Gauss-Bonnet Gravity. Phys. Rev. Lett. 2021, 127, 101102. [Google Scholar] [CrossRef]
  20. Annulli, L. CLAP for modified gravity: Scalar instabilities in binary black hole spacetimes. arXiv 2021, arXiv:2105.08728. [Google Scholar]
  21. Doneva, D.D.; Yazadjiev, S.S. Beyond the spontaneous scalarization: New fully nonlinear dynamical mechanism for formation of scalarized black holes. arXiv 2021, arXiv:2107.01738. [Google Scholar]
  22. Dima, A.; Barausse, E.; Franchini, N.; Sotiriou, T.P. Spin-induced black hole spontaneous scalarization. Phys. Rev. Lett. 2020, 125, 231101. [Google Scholar] [CrossRef] [PubMed]
  23. Herdeiro, C.A.R.; Radu, E.; Silva, H.O.; Sotiriou, T.P.; Yunes, N. Spin-induced scalarized black holes. Phys. Rev. Lett. 2021, 126, 011103. [Google Scholar] [CrossRef]
  24. Berti, E.; Collodel, L.G.; Kleihaus, B.; Kunz, J. Spin-induced black-hole scalarization in Einstein-scalar-Gauss-Bonnet theory. Phys. Rev. Lett. 2021, 126, 011104. [Google Scholar] [CrossRef] [PubMed]
  25. Cunha, P.V.; Herdeiro, C.A.; Radu, E. Spontaneously Scalarized Kerr Black Holes in Extended Scalar-Tensor–Gauss-Bonnet Gravity. Phys. Rev. Lett. 2019, 123, 011101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Hod, S. Onset of spontaneous scalarization in spinning Gauss-Bonnet black holes. Phys. Rev. D 2020, 102, 084060. [Google Scholar] [CrossRef]
  27. Herdeiro, C.A.R.; Radu, E.; Sanchis-Gual, N.; Font, J.A. Spontaneous Scalarization of Charged Black Holes. Phys. Rev. Lett. 2018, 121, 101102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Brihaye, Y.; Hartmann, B. Spontaneous scalarization of charged black holes at the approach to extremality. Phys. Lett. B 2019, 792, 244–250. [Google Scholar] [CrossRef]
  29. Townsend, P.K. Black holes: Lecture notes. arXiv 1997, arXiv:gr-qc/9707012. [Google Scholar]
  30. Herdeiro, C.; Radu, E. Construction and physical properties of Kerr black holes with scalar hair. Class. Quantum Gravity 2015, 32, 144001. [Google Scholar] [CrossRef] [Green Version]
  31. Fernandes, P.G.S.; Herdeiro, C.A.R.; Pombo, A.M.; Radu, E.; Sanchis-Gual, N. Spontaneous Scalarisation of Charged Black Holes: Coupling Dependence and Dynamical Features. Class. Quantum Gravity 2019, 36, 134002. [Google Scholar] [CrossRef] [Green Version]
  32. Astefanesei, D.; Herdeiro, C.; Pombo, A.; Radu, E. Einstein-Maxwell-scalar black holes: Classes of solutions, dyons and extremality. J. High Energy Phys. JHEP 2019, 10, 78. [Google Scholar] [CrossRef] [Green Version]
  33. Gibbons, G.W.; Maeda, K.I. Black Holes and Membranes in Higher Dimensional Theories with Dilaton Fields. Nucl. Phys. B 1988, 298, 741–775. [Google Scholar] [CrossRef]
  34. Garfinkle, D.; Horowitz, G.T.; Strominger, A. Charged black holes in string theory. Phys. Rev. D 1991, 43, 3140. [Google Scholar] [CrossRef]
  35. Sen, A. Rotating charged black hole solution in heterotic string theory. Phys. Rev. Lett. 1992, 69, 1006–1009. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Volkov, M.S.; Galtsov, D.V. Nonabelian Einstein Yang-Mills Black Holes. Sov. J. Exp. Theor. Phys. Lett. 1989, 50, 346–350. [Google Scholar]
  37. Kuenzle, H.P.; Masood-ul-Alam, A.K. Spherically Symmetric Static SU(2) Einstein Yang-Mills Fields. J. Math. Phys. 1990, 31, 928–935. [Google Scholar] [CrossRef]
  38. Bizon, P. Colored Black Holes. Phys. Rev. Lett. 1990, 64, 2844. [Google Scholar] [CrossRef] [PubMed]
  39. Volkov, M.S.; Gal’tsov, D.V. Gravitating non-Abelian solitons and black holes with Yang-Mills fields. Phys. Rep. 1999, 319, 1–83. [Google Scholar] [CrossRef] [Green Version]
  40. Volkov, M.S. Hairy black holes in the XX-th and XXI-st centuries. In Proceedings of the MG14 Meeting on General Relativity, Rome, Italy, 12–18 July 2015. [Google Scholar]
  41. Bartnik, R.; Mckinnon, J. Particlelike Solutions of the Einstein Yang-Mills Equations. Phys. Rev. Lett. 1988, 61, 141–144. [Google Scholar] [CrossRef] [PubMed]
  42. Doneva, D.D.; Staykov, K.V.; Yazadjiev, S.S. Gauss-Bonnet black holes with a massive scalar field. Phys. Rev. D 2019, 99, 104045. [Google Scholar] [CrossRef] [Green Version]
  43. Macedo, C.F.B.; Sakstein, J.; Berti, E.; Gualtieri, L.; Silva, H.O.; Sotiriou, T.P. Self-interactions and Spontaneous Black Hole Scalarization. Phys. Rev. D 2019, 99, 104041. [Google Scholar] [CrossRef] [Green Version]
Figure 1. GB invariant (in units of mass) of a Kerr-Newman BH with mass M, angular momentum J and charge Q, evaluated at the horizon’s poles ( r = r h and θ = 0 ), as a function of the dimensionless parameters j = J / M 2 and q = Q / M . Near extremality, R G B 2 < 0 .
Figure 1. GB invariant (in units of mass) of a Kerr-Newman BH with mass M, angular momentum J and charge Q, evaluated at the horizon’s poles ( r = r h and θ = 0 ), as a function of the dimensionless parameters j = J / M 2 and q = Q / M . Near extremality, R G B 2 < 0 .
Universe 07 00483 g001
Figure 2. (Left panel) Dimensionless ratio λ / M of the set of RN solutions supporting the = m = n = 0 scalar cloud vs. the reduced electric charge q for ϵ = ± 1 . (Right panel) Typical radial profiles of the spherical, nodeless scalar clouds on a RN BH background.
Figure 2. (Left panel) Dimensionless ratio λ / M of the set of RN solutions supporting the = m = n = 0 scalar cloud vs. the reduced electric charge q for ϵ = ± 1 . (Right panel) Typical radial profiles of the spherical, nodeless scalar clouds on a RN BH background.
Universe 07 00483 g002
Figure 3. Charge-temperature (main panels panel) and charge-horizon area (insets) diagrams, in units set by mass, for several illustrative families of GB ϵ scalarised RN solutions, for both ϵ = + 1 (left panel) and ϵ = 1 (right panel). The branches of scalarised solutions bifurcate from the electrovacuum RN BHs (blue line) and terminate in critical configurations (red circles).
Figure 3. Charge-temperature (main panels panel) and charge-horizon area (insets) diagrams, in units set by mass, for several illustrative families of GB ϵ scalarised RN solutions, for both ϵ = + 1 (left panel) and ϵ = 1 (right panel). The branches of scalarised solutions bifurcate from the electrovacuum RN BHs (blue line) and terminate in critical configurations (red circles).
Universe 07 00483 g003
Figure 4. The GB invariant at the horizon (in units of mass) as a function of the θ coordinate for several Kerr-Newman BHs.
Figure 4. The GB invariant at the horizon (in units of mass) as a function of the θ coordinate for several Kerr-Newman BHs.
Universe 07 00483 g004
Figure 5. Branches of GB scalarised Kerr-Newman BHs in a q vs. j plot (left panel) and j vs. λ / M plot (right panel). The branches are for a specific choice of Ω H λ = 0.469 and different choices of Q / λ . The blue squares correspond to Kerr-Newman BHs with a vanishing scalar field, while the red circles correspond to critical configurations.
Figure 5. Branches of GB scalarised Kerr-Newman BHs in a q vs. j plot (left panel) and j vs. λ / M plot (right panel). The branches are for a specific choice of Ω H λ = 0.469 and different choices of Q / λ . The blue squares correspond to Kerr-Newman BHs with a vanishing scalar field, while the red circles correspond to critical configurations.
Universe 07 00483 g005
Figure 6. The profile of typical scalar clouds are shown for two Kerr-Newman backgrounds with different values of the reduced charge and angular momentum ( q , j ) .
Figure 6. The profile of typical scalar clouds are shown for two Kerr-Newman backgrounds with different values of the reduced charge and angular momentum ( q , j ) .
Universe 07 00483 g006
Figure 7. (Left panel) GB invariant as a function of the radial coordinate for an EYM BH with horizon radius r h = 0.6 . One notices a shell with R G B 2 < 0 , located outside the horizon. (Right panel) λ / M of the set of EYM BHs supporting scalar clouds as a function of the reduced horizon area for ϵ = ± 1 .
Figure 7. (Left panel) GB invariant as a function of the radial coordinate for an EYM BH with horizon radius r h = 0.6 . One notices a shell with R G B 2 < 0 , located outside the horizon. (Right panel) λ / M of the set of EYM BHs supporting scalar clouds as a function of the reduced horizon area for ϵ = ± 1 .
Universe 07 00483 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Herdeiro, C.A.R.; Pombo, A.M.; Radu, E. Aspects of Gauss-Bonnet Scalarisation of Charged Black Holes. Universe 2021, 7, 483. https://0-doi-org.brum.beds.ac.uk/10.3390/universe7120483

AMA Style

Herdeiro CAR, Pombo AM, Radu E. Aspects of Gauss-Bonnet Scalarisation of Charged Black Holes. Universe. 2021; 7(12):483. https://0-doi-org.brum.beds.ac.uk/10.3390/universe7120483

Chicago/Turabian Style

Herdeiro, Carlos A. R., Alexandre M. Pombo, and Eugen Radu. 2021. "Aspects of Gauss-Bonnet Scalarisation of Charged Black Holes" Universe 7, no. 12: 483. https://0-doi-org.brum.beds.ac.uk/10.3390/universe7120483

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop