Charged Gauss–Bonnet black holes supporting non-minimally coupled scalar clouds: analytic treatment in the near-critical regime

Recent numerical studies have revealed the physically intriguing fact that charged black holes whose charge-to-mass ratios are larger than the critical value (Q/M)crit=2(9+6)/5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(Q/M)_{\text {crit}}=\sqrt{2(9+\sqrt{6})}/5$$\end{document} can support hairy matter configurations which are made of scalar fields with a non-minimal negative coupling to the Gauss–Bonnet invariant of the curved spacetime. Using analytical techniques, we explore the physical and mathematical properties of the composed charged-black-hole-nonminimally-coupled-linearized-massless-scalar-field configurations in the near-critical Q/M≳(Q/M)crit\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Q/M\gtrsim (Q/M)_{\text {crit}}$$\end{document} regime. In particular, we derive an analytical resonance formula that describes the charge-dependence of the dimensionless coupling parameter η¯crit=η¯crit(Q/M)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{\eta }_{\text {crit}}=\bar{\eta }_{\text {crit}}(Q/M)$$\end{document} of the composed Einstein–Maxwell-nonminimally-coupled-scalar-field system along the existence-line of the theory, a critical border that separates bald Reissner–Nordström black holes from hairy charged-black-hole-scalar-field configurations. In addition, it is explicitly shown that the large-coupling -η¯crit(Q/M)≫1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-\bar{\eta }_{\text {crit}}(Q/M)\gg 1$$\end{document} analytical results derived in the present paper for the composed Einstein–Maxwell-scalar theory agree remarkably well with direct numerical computations of the corresponding black-hole-field resonance spectrum.


I. INTRODUCTION
Early no-hair theorems [1][2][3][4][5][6], which were mainly motivated by Wheeler's no-hair conjecture for black holes [7,8], have revealed the physically interesting fact that black holes with spatially regular horizons cannot support static matter configurations which are made of scalar fields.
The physically intriguing phenomenon of spontaneous scalarization of black holes in generalized Einstein-Gauss-Bonnet-scalar field theories [9][10][11][12][13][14][15][16] owe its existence to the presence of an effective spatially-dependent mass term of the linearized form −ηG, which reflects the direct non-trivial coupling between the scalar field and the Gauss-Bonnet curvature invariant, in the Klein-Gordon wave equation of the supported scalar configurations [see Eq. (10) below].The dimensionless physical parameter η of the composed field theory controls the strength of the direct non-minimal interaction between the scalar field and the spatially-dependent Gauss-Bonnet invariant of the curved spacetime.
The spontaneous scalarization phenomenon of charged black holes in composed Einstein-Maxwell-Gauss-Bonnetscalar-field theories has been explored, using numerical techniques, in the physically interesting works [25,26] (see [13,14,[27][28][29][30][31][32] for recent studies of the spontaneous scalarization phenomenon of asymptotically flat spinning black holes).Intriguingly, it has been revealed in [25,26] that the composed charged-black-hole-nonminimally-coupledscalar-field system is characterized by a charge-dependent existence-line η = η( Q) [33] which marks the onset of the spontaneous scalarization phenomenon in the generalized Einstein-Maxwell-scalar field theory.In particular, for a given value of the black-hole electric charge Q, the critical existence line marks the boundary between bald Reissner-Nordström black holes and composed charged-black-hole-nonminimally-coupled-scalar-field hairy configurations.
The charge-dependent critical existence-line η = η( Q) of the generalized Einstein-Maxwell-scalar theory is composed of charged Reissner-Nordström black holes that support scalar 'clouds' [34,35], spatially regular matter configurations which are made of the non-trivially coupled linearized scalar fields.The characteristic existence-line of the physical system is universal in the sense that different non-trivially coupled Einstein-Maxwell-scalar field theories that share the same weak-field functional behavior f (φ) = 1 + ηφ 2 /2 + O(φ 4 ) [25,26] of the coupling function are characterized by the same critical boundary between bald and hairy black-hole spacetimes.
Interestingly, the numerical results presented in [25,26] have revealed that the spontaneous scalarization phenomenon of black holes in composed Einstein-Maxwell-scalar field theories with negative values of the non-minimal coupling parameter η may be induced by the electric charge of the supporting black hole.In particular, one finds that, in the η < 0 regime, the onset of the spontaneous scalarization phenomenon is marked by the dimensionless black-hole electric charge [25,26] Qcrit = 2(9 + √ 6)

5
. ( Only charged black holes with Q ≥ Qc can support non-minimally coupled spatially regular scalar clouds in the negative coupling η < 0 regime.Intriguingly, it has been demonstrated numerically in [25,26] that, in the nearcritical regime Q/ Qcrit → 1 + , the hairy charged black holes are characterized by the large-coupling asymptotic relation The main goal of the present paper is to study, using analytical techniques, the physical and mathematical properties of the composed charged-Reissner-Nordström-black-hole-nonminimally-coupled-scalar-field cloudy configurations in the near-critical regime Q > ∼ Qcrit .In particular, a remarkably compact WKB resonance formula, which describes the charge-dependence η = η( Q) of the characteristic critical existence-line of the composed Einstein-Maxwell-Gauss-Bonnet-nonminimally-coupled-massless-scalar-field theory in the dimensionless large-coupling −η ≫ 1 regime, will be derived.
Interestingly, we shall explicitly show below that the analytically derived near-critical resonance formula of the present paper [see Eq. ( 40) below] provides a simple analytical explanation for the numerically observed [25,26] monotonically decreasing functional behavior |η| = |η( Q)| of the critical existence-line that characterizes, in the η < 0 regime, the composed charged-Reissner-Nordström-black-hole-nonminimally-coupled-linearized-scalar-field configurations.

II. DESCRIPTION OF THE SYSTEM
We explore the physical and mathematical properties of spatially regular scalar 'clouds' (linearized scalar field configurations) which are supported by charged Reissner-Nordström black holes.The supported massless scalar fields are characterized by a direct non-trivial (non-minimal) coupling to the Gauss-Bonnet invariant of the curved charged spacetime.The action of the composed Einstein-Maxwell-Gauss-Bonnet-nonminimally-coupled-massless-scalar-field system is given by the expression [26,36] where is the Gauss-Bonnet curvature invariant.Following [11,12,18], we assume that the scalar function f (φ), which controls the non-trivial direct coupling of the scalar field to the Gauss-Bonnet curvature invariant, is characterized by the leading order universal functional behavior in the weak-field regime.As discussed in [11,12,18], the functional behavior (5) of the scalar coupling function in the weak scalar field regime guarantees that the Einstein-matter field equations are satisfied by the familiar scalarless black-hole solutions of general relativity (in our case, the Reissner-Nordström black-hole spacetime) in the φ → 0 limit.The physical parameter η [37], which may take either positive or negative values [25,26], determines the strength of the direct (non-minimal) interaction between the spatially regular massless scalar field configurations and the Gauss-Bonnet invariant (4) of the charged curved spacetime.
The supporting charged Reissner-Nordström black-hole spacetime is characterized by the line element [38][39][40]] where the metric function ∆ is given by the functional expression The physical parameters {M, Q} are respectively the mass and electric charge of the central supporting Reissner-Nordström black hole.The horizon radii of the curved black-hole spacetime (6) are determined by the zeros of the metric function ∆: The radially-dependent Gauss-Bonnet invariant, which characterizes the curved Reissner-Nordström black-hole spacetime (6), is given by the expression [26] A variation of the action (3) with respect to the scalar field yields the generalized Klein-Gordon equation [26] for the non-minimally coupled scalar field with the radially-dependent effective mass term This effective mass term reflects the direct (non-minimal) coupling of the scalar field φ to the Gauss-Bonnet invariant G [see Eq. ( 9)] of the charged curved spacetime.Interestingly, one finds that, in the regime [26,41] the charged-dependent Gauss-Bonnet curvature invariant becomes negative in the exterior region of the charged black-hole spacetime.Thus, the spatially-dependent effective mass term (11) with η < 0 may become negative in the interval (13).As we shall explicitly prove below, this intriguing property of the coupled Einstein-Maxwell-nonminimally-coupled-scalar-field system (3) allows the central charged Reissner-Nordström black hole (6) to support bound-state linearized cloudy configurations of the non-minimally coupled scalar field.
Using the functional field decomposition [42] one obtains from Eq. ( 10) the differential equation [26,43,44] for the radial part of the linearized non-minimally coupled scalar field in the supporting curved black-hole spacetime (6).
The ordinary differential equation ( 15) determines the spatial behavior of the static non-minimally coupled linearized massless scalar field configurations in the supporting charged Reissner-Nordström black-hole spacetime (6).Following [25,26], we shall assume that the scalar eigenfunction ψ(r) of the bound-state field configurations is spatially well behaved with the physically motivated boundary conditions [25,26] at the black-hole horizon and at spatial infinity.
In the next section we shall study, using analytical techniques, the physical and mathematical properties of the negatively coupled charged-black-hole-nonminimally-coupled-linearized-massless-scalar-field configurations that characterize the composed Einstein-Maxwell-scalar field theory (3).In particular, we shall explicitly prove that the composed Reissner-Nordström-black-hole-scalar-field cloudy configurations in the Q > Q crit regime [see Eq. ( 12)] with η < 0 owe their existence to the presence of an effective near-horizon binding potential well [see Eq. ( 31) below].

III. COMPOSED REISSNER-NORDSTR ÖM-BLACK-HOLE-NONMINIMALLY-COUPLED-SCALAR-FIELD CONFIGURATIONS: A WKB ANALYSIS
In the present section we shall determine the charge-dependent functional behavior η = η(Q/M ) of the critical existence-line, which characterizes the composed charged-black-hole-nonminimally-coupled-scalar-field cloudy configurations, in the dimensionless large-coupling regime Defining the new radial eigenfunction and using the coordinate y(r), which is defined by the radial differential relation [45] one obtains from equation ( 15) the Schrödinger-like radial equation where the radially-dependent potential V [r(y)] of the composed charged-Reissner-Nordström-black-holenonminimally-coupled-scalar-field system is given by the functional expression with [see Eq. ( 7)] The presence of the non-trivial term in the effective interaction potential ( 21) is a direct consequence of the non-trivial (non-minimal) coupling between the Gauss-Bonnet curvature invariant (9) of the charged spacetime (6) and the supported scalar field.
We shall now prove that, in the dimensionless large-coupling regime (17), the ordinary differential equation (20), which characterizes the composed Einstein-Maxwell-scalar field theory (3), is amenable to an analytical treatment.In particular, a standard second-order WKB analysis of the Schrödinger-like ordinary differential equation ( 20) yields the familiar quantization condition [46][47][48] for the discrete spectrum {η(M, Q; n} n=∞ n=0 of the dimensionless coupling parameter which characterizes the composed charged-Reissner-Nordström-black-hole-nonminimally-coupled-scalar-field cloudy configurations.Here n ∈ {0, 1, 2, ...} is the discrete resonant parameter of the physical system.The critical existence-line of the field theory is determined by the fundamental n = 0 mode.The integration limits {y t− , y t+ } in the WKB integral relation (24) are the classical turning points of the effective binding potential (21).Using Eq. ( 19), one finds that the WKB relation (24) can be expressed in the integral form Defining the dimensionless physical parameters and one obtains the near-critical (ǫ ≪ 1) near-horizon (x ≪ 1) relations [see Eqs. ( 8), (12), and ( 22)] and Substituting the relations ( 28), (29), and (30) into Eqs.( 21) and ( 23), one finds the (rather cumbersome) nearcritical near-horizon expression for the effective interaction potential of the composed black-hole-massless-scalar-field cloudy configurations.We shall henceforth consider composed charged-Reissner-Nordström-black-hole-massless-scalar-field cloudy configurations in the dimensionless large-coupling regime [see Eqs. ( 17), ( 26), and Eq. ( 38) below] in which case one finds from Eqs. ( 30) and ( 31) the remarkably compact leading order functional expression for the effective binding potential of the composed charged-black-hole-scalar-field system.Taking cognizance of Eqs. ( 25), (27), and (33), one finds the resonance condition for the composed black-hole-field system.The WKB integral relation (34) can be expressed in the compact mathematical form where Using the integral relation one finds from (35) the remarkably simple discrete resonance formula [49] η(ǫ; which characterizes the composed charged-Reissner-Nordström-black-hole-nonminimally-coupled-linearized-masslessscalar-field cloudy configurations in the dimensionless large-coupling −η ≫ 1 regime (or equivalently, in the nearcritical ǫ ≪ 1 regime).The discrete resonance spectrum (38) of the composed Einstein-Maxwell-Gauss-Bonnet-nonminimally-coupledscalar-field system (3) can be expressed in the dimensionless form [see Eqs. ( 1) and ( 26 where Q ≡ Q/M .

IV. NUMERICAL CONFIRMATION
In the present section we shall test the accuracy of the analytically derived large-coupling resonance spectrum (38), which characterizes the composed charged-Reissner-Nordström-black-hole-nonminimally-coupled-linearized-scalarfield cloudy configurations.The charge-dependent resonance spectrum of the black-hole-field system has recently been computed numerically in [26].
In Table I we present, for various values of the dimensionless coupling parameter λ ≡ 2 |η| [44] used in [26], the charge-dependent ratio R(λ) ≡ ǫ analytical /ǫ numerical between the analytically calculated value of the dimensionless critical parameter ǫ [as calculated directly from the analytically derived large-coupling resonance formula (38)] and the corresponding exact (numerically computed [26]) values of the critical parameter.The data presented in Table I for the cloudy Reissner-Nordström-black-hole-scalar-field configurations reveals the fact that the agreement between the analytically derived resonance formula (38) and the corresponding numerically computed resonance spectrum of [26] is remarkably good in the large-coupling λ ≫ 1 regime [50] of the Einstein-Maxwell-scalar field theory (3).[26,44] of the theory, the dimensionless ratio R(λ) ≡ ǫ analytical /ǫ numerical between the analytically calculated value of the dimensionless critical parameter ǫ [as calculated directly from the resonance formula (38) for the fundamental n = 0 mode] and the corresponding exact values of the critical parameter as computed numerically in [26].One finds that, in the large-coupling λ ≫ 1 regime of the composed Reissner-Nordström-black-hole-scalar-field cloudy configurations, the agreement between the analytically derived resonance formula (38) and the corresponding numerically computed resonance spectrum of [26] is remarkably good [50].
The spontaneous scalarization phenomenon of charged black holes in composed Einstein-Maxwell-Gauss-Bonnetscalar field theories has recently been studied numerically in the physically important works [25,26].In particular, it has been revealed in [25,26] that a charge-dependent critical existence-line ηcrit = ηcrit (Q/M ) separates bald Reissner-Nordström black-hole spacetimes from the composed charged-black-hole-nonminimally-coupled-massless-scalar-field hairy configurations of the Einstein-Maxwell-scalar theory (3) [51], where the dimensionless physical parameter η quantifies the strength of the direct non-trivial interaction between the supported scalar field and the Gauss-Bonnet curvature invariant.
Interestingly, it has been demonstrated [25,26] that, in the negative coupling η < 0 regime, the composed chargedblack-hole-linearized-scalar-field cloudy configurations that sit on the critical existence-line of the Einstein-Maxwellscalar field theory (3) are restricted to the dimensionless charge regime Q > Qcrit = 2(9 + √ 6)/5 [see Eq. ( 1)].In particular, the numerical results presented in [25,26] provide important evidence for an intriguing divergent functional behavior −η( Q) → ∞ of the non-minimal coupling parameter of the theory along the existence-line of the system in the near-critical limit Q/ Qcrit → 1 + .
In the present paper we have used analytical techniques in order to explore the physical and mathematical properties of the composed charged-Reissner-Nordström-black-hole-nonminimally-coupled-massless-scalar-field hairy configurations of the Einstein-Maxwell-scalar theory (3) in the large-coupling −η( Q) ≫ 1 regime.In particular, using a WKB procedure we have derived the analytical formula (38) for the discrete resonant spectrum of the dimensionless coupling parameter η( Q) which characterizes the composed charged-black-hole-scalar-field cloudy configurations in the near-critical Q > ∼ Q crit regime.The analytically derived discrete WKB resonance formula (38) yields the remarkably compact charge-dependent functional relation [52] η for the critical existence-line of the composed Einstein-Maxwell-Gauss-Bonnet-nonminimally-coupled-massless-scalar field theory (3), where 0 ≤ ǫ = ( Q − Qcrit )/ Qcrit ≪ 1 [see Eqs. ( 1) and ( 26)] is the dimensionless distance of the system from the exact critical (−η → ∞) configuration.
It is worth stressing the fact that the physical significance of the critical existence-line (40) in the Einstein-Maxwellscalar field theory (3) stems from the fact that it separates, in the large-coupling −η ≫ 1 regime, bald Reissner-Nordström black holes from hairy charged-black-hole-nonminimally-coupled-massless-scalar-field configurations.From the analytically derived resonance formula (40) one finds, in accord with the physically interesting numerical results presented in [25,26], that the dimensionless coupling parameter |η( Q)| of the theory (3) is a monotonically decreasing function of the black-hole electric charge Q in the negative coupling η < 0 regime.
Finally, it is interesting to point out that our results provide a simple analytical explanation for the physically intriguing numerical observation originally made in [25,26], according to which the charge-dependent coupling parameter η( Q) of the Einstein-Maxwell-scalar field theory (3) diverges along the critical existence line of the system in the Q → Qcrit limit.In particular, the analytically derived resonance formula (40) reveals the fact that the coupling parameter η( Q), which characterizes the cloudy black-hole-field configurations, diverges quadratically in the near-critical ǫ → 0 limit.