{'fraction_non_alphanumeric': 0.03199076555221173, 'fraction_numerical': 0.01372799604237952, 'mean_word_length': 4.606193667668131, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 9, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Quadratic gravity and conformally coupled scalar fields

We construct black hole solutions in four-dimensional quadratic gravity, supported by a scalar field conformally coupled to quadratic terms in the curvature. The conformal matter Lagrangian is constructed with powers of traces of a conformally covariant tensor, which is defined in terms of the metric and a scalar field, and has the symmetries of the Riemann tensor. We find exact, neutral and charged, topological black hole solutions of this theory when the Weyl squared term is absent from the action functional. Including terms beyond quadratic order on the conformally covariant tensor, allows to have asymptotically de Sitter solutions, with a potential that is bounded from below. For generic values of the couplings we also show that static black hole solutions must have a constant Ricci scalar, and provide an analysis of the possible asymptotic behavior of both, the metric as well as the scalar field in the asymptotically AdS case, when the solutions match those of general relativity in vacuum at infinity. In this frame, the spacetime fulfils standard asymptotically AdS boundary conditions, and in spite of the non-standard couplings between the curvature and the scalar field, there is a family of black hole solutions in AdS that can be interpreted as localized objects. We also provide further comments on the extension of these results to higher dimensions.


I. INTRODUCTION
Higher curvature terms in addition to the Einstein-Hilbert Lagrangian appear in many different settings. The precise form of those correction depends on the specific details of the theory as well as on the frame used in the perturbative scheme that give rise to them. Higher curvature correction in four dimensions improve the renormalizability properties of the theory at the cost of spoiling its unitarity [1,2], a feature that may be absent when the whole series of corrections is considered. A particularly interesting theory with quadratic terms in the curvature in four dimensions is Critical Gravity whose Lagrangian is a precise combination of the Einstein-Hilbert term plus the Weyl square term. In this theory the coupling of the latter is fixed in terms of the negative, bared cosmological constant [3]. Restricting the fall-off of the asymptotically AdS solutions in this theory, it was argued that only Einstein manifolds may remain in the spectrum [4]. On the other hand, in asymptotically flat scenarios, the existence of black holes in General Relativity plus the most general quadratic theory in four dimensions, leads to spacetimes with vanishing Ricci scalar [5,6]. Such argument strongly relies on the tracelessness of the Bach tensor that is obtained as the Euler-Lagrange derivative of the conformal gravity action √ −gC µνλρ C µνλρ [5] (C µ νλρ being the conformally invariant Weyl tensor). Such theory, trivially admits Einstein manifolds as solutions since the Bach tensor vanishes on spacetimes that are conformally related to the former. The spherically symmetric solutions that non-trivially depart from the Schwarzschild black hole can be integrated numerically [5,6], and can also be obtained in terms of an infinite recurrence relation, which has a closed form [7][8][9]. Such construction can be extended to include a Maxwell field [10] due to the tracelessness of its energy-momentum tensor and a natural question therefore arises: is it possible to consider other matter sources in this construction?
Here we explore the consequences of introducing a conformally coupled scalar field in the action which is particularly relevant since conformally coupled scalars have been useful in the construction of black holes with secondary hair in GR [11][12][13][14][15][16][17][18][19][20][21][22][23]. Since the gravity Lagrangian contains operators up to mass dimension 4, we mimic the same structure in the matter sector, respecting a symmetry under local Weyl rescaling. A useful manner for doing so is to make use of the tensor that transforms under local Weyl rescalings as This tensor was introduced in [24] and was used to construct conformal couplings of a scalar field to Euler densities of a higher degree, leading to the conformally coupled version of Lovelock gravities [25]. These couplings have allowed to construct exact black holes with secondary hair in GR in higher dimensions [18], with interesting thermodynamic properties [19][20][21][26][27][28][29][30][31][32], extending to dimension d > 4 the known solutions in three and four dimensions of Einstein gravity and by-passing the no-go results [33,34].
Here we remain in dimension four, and introduce higher derivative terms constructed out from the tensor S µν λρ . The traces of the tensor S will be important, so we define S α β := S αγ βγ and S := S α α , which read Fixing the conformal weight of the scalar field in (2) to s = −1, one recovers, up to a global sign, the Lagrangian for the standard conformally coupled scalar field from the scalar quantity S defined in (4), with a canonical kinetic term [35]. In reference [24], the tensor S µν λρ defined in (1) was proven to be the unique, non-trivial tensorial combination of the scalar field and derivatives of the metric, that scales homogeneously under local Weyl rescalings and has the symmetries of the Riemann tensor. This tensor can also be obtained introducing a Weyl compensator for the metric which was proved useful in the construction of boundary terms in higher dimension [21] 1 . The inclusion of terms in the action that break the conformal invariance remove the pure-gauge nature of the field φ, making impossible to gauge it away. Similar higher curvature terms can be constructed using a Weyl-invariant curvature as explored at the quadratic level in arbitrary dimension in [37][38][39][40], where a gauge field A µ on top of the metric and the scalar have to be considered. In such scenario the whole 1 At the level of the Lagrangian for the standard conformally coupled scalar field, this was done in [36].
Lagrangian is Weyl invariant and therefore the scalar can be gauge fixed by a conformal transformation.
The paper is organized as follows: In Section II we present the theory consisting on the most general quadratic combination in the curvature and in the tensor S µν λρ , maintaining the conformal invariance of the matter sector. We show that due to the existence of some identities in four dimensions, the Lagrangian severely reduces. In Section III we show that the static black holes of the theory must have a constant Ricci scalar, forcing a second order constraint on the metric in a theory that has fourth-order field equations. Section IV is devoted to the construction of exact solutions describing spherically symmetric black holes in de Sitter (dS), or topological black holes with compact hyperbolic horizons in anti de Sitter (AdS). We also provide the dyonically charged version of these black holes. Section V contains a general discussion on the asymptotic behavior of the metric and the self-interacting scalar, with derivative self-interactions, admitting as a possible asymptotic behavior the standard Henneaux-Teitelboim asymptotic conditions [41]. Section VI contains an extension of these results when powers of the form S k>2 are included maintaining the conformal invariance in the matter sector. We show that asymptotically de Sitter black holes exist even when the potential for the scalar field is bounded from below. Partial generalizations of these results to higher, even dimensions, as well as further comments are given in Section VII.

II. QUADRATIC GRAVITY AND CONFORMALLY COUPLED SCALARS
We consider the theory where where C µνρσ is the Weyl tensor and W µνρσ (S) is the traceless part of the tensor S µνρσ , consequently defined as Some remarks are now in order: in the gravitational action (7), we have explicitly considered the most general quadratic theory in a particular basis of quadratic invariants. In four dimensions the term proportional to α 3 is the Euler density and therefore, its variation does not contribute to the field equations. Since we are interested only on the properties of the solutions of the theory (5), hereafter we set α 3 = 0. The matter action (8) is invariant where we have fixed the conformal weight of the scalar field to s = −1, which leads to a canonical kinetic term from the scalar field which comes from the term linear in S. We have also mimicked the same structure of the gravitational Lagrangian in the higher curvature terms. Since the scalar field has mass dimension −1, all the couplings of the quadratic terms, α i and β i in (7) and (8), are dimensionless. Therefore, these are all the terms that can appear up to quadratic order in an action, that preserve the conformal symmetry of the matter sector.
After explicitly computing the last two terms in (8), one obtains b.t. standing for boundary terms.
The first equality (10) is expected due to the fact mentioned above regarding the tensor S αβγδ as change of frame of R αβγδ . Consequently, from the point of view of the field equations it is enough to consider the following action principle where we redefined α 2 → α 2 − β 2 , and we have introduced Newton's constant. Here after we will work with this action.
The field equations coming from the action (12), taking variations with respect to the metric take a simple form The term proportional to α 2 is the Bach tensor, while the first three terms in the third line of (13) define the standard, quadratic, improved energy-momentum tensor [35]. On the other hand, the equation for the higher curvature, conformally coupled scalar reads As expected due to the conformal symmetry of the matter Lagrangian, the trace of the energy-momentum that can be read from equation (13), vanishes when one uses the equation for the scalar field (14). In the following sections we will be interested in deriving some properties of the black hole solutions of the system (13), (14). The trace of the field equations reduces to a wave equation for R − 4Λ, namely As it is well known for vanishing Λ in vacuum, equation (15) after linearisation around flat spacetime leads to a massive degree of freedom of gravitational origin, with mass given by m 2 eff = (6α 1 ) −1 , and therefore α 1 > 0 [1, 2] (see also [42]). It is interesting to notice that for negative cosmological constant Λ = −3/l 2 , the coupling α 1 could be negative, provided the Breitenlohner-Freedman (BF) [43,44] bound is fulfilled. Notice that this is not consistent with a perturbative approach since in such scenario |α 1 | << 1, and therefore the mass may be a negative large number, violating the BF bound. We therefore consider α 1 > 0.

III. RESTRICTING THE BLACK HOLE SOLUTIONS
An static black hole spacetime admits the following coordinates x µ = (t, x i ), where the metric takes the form Here N (x i ) is a function that vanishes on the would-be event horizon, while h ij is a regular, spacelike metric in the domain of outer communications for asymptotically flat and asymptotically AdS black holes. Now, the standard argument of no-hair theorems for asymptotically flat black holes can be extended to other asymptotic behavior. Let us define the quantity ξ = R − 4Λ, which by virtue of staticity depends only on the coordinates x i . We will assume that the conformal scalar φ decays at infinity such that R → −4Λ fast enough.
The wave equation in (15) can be integrated on the exterior region of the black hole, which after integration by parts, leads to The first term in the integrand contributes as a boundary term, which acquires two contributions at the would-be horizon and at infinity. Since N vanishes at the horizon, the regularity of ξ implies that both contributions vanish since we have also assumed that ξ approaches zero at infinity, fast enough. Since N is positive in the exterior region of asymptotically flat and AdS black holes and the metric h ij is positive definite in that region, we find that (17) implies that ξ has to vanish everywhere, and therefore, that the whole spacetime must have a constant Ricci scalar R = 4Λ. While the theory is of fourth order, this restriction provides a second order constraint that helps in the explicit integration of the field equation. Notice that for asymptotically de Sitter black hole, this argument can be directly extended leading to the same conclusion that R = 4Λ in the whole region between the event and cosmological horizons. Indeed, in that case, we must integrate in the region between the horizons and since the boundary term will vanish at both horizons we would reach the same conclusion 2 .
Assuming analyticity allows to extend this constraint to the whole spacetime.

A. Neutral case
If we allow for a single metric function in a Schwarzschild like ansatz, the integration of the constraint R = 4Λ that we have derived in the previous section, fixes the metric function allowing to construct an exact solution for the whole system of equation. In fact, the metric 2 Notice that a similar argument has been used in [45] to restrict black hole solutions in R 2 supergravity [46]. and scalar field define a solution of the equations (13) and (14), provided α 2 = 0, and λ = − 2πGΛ 9(1 + 128πGΛ(α 1 + β)) .
In the line element (18), dΩ γ denotes the metric of an Euclidean spacetime of constant curvature γ. Notice that while the coupling of the Weyl squared term in (5) has to be turned-off, this solution does receive a non-trivial contribution from the quadratic terms with couplings α 1 and β in (5).
The action (12) contains higher curvature terms as well as conformal couplings of the scalar with quadratic curvature terms, nevertheless the form of the metric function, and therefore the causal structure of this spacetime coincides with that reported in [14,15]. In the spherically symmetric case γ = 1 and therefore cosmic censorship implies Λ > 0. Provided the parameter µ is positive and below a certain critical value µ c , the causal structure is that of Reissner-Norstrom spacetime in de Sitter space, while for µ = µ c the event horizon coincides with the cosmological horizon. When dΩ γ is an smooth quotient of the hyperbolic space H 2 [47], one can have event horizons even in the asymptotically AdS case, containing a rich family of causal structures that include black holes with Cauchy horizons, as well as spacetimes that can be interpreted black holes inside a black hole [15]. Notice also that the reality of the scalar field φ(r) in (19) implies that the self-interaction potential λφ 4 will be bounded from below only in the asymptotically AdS case, i.e. when Λ < 0. We will see below that the inclusion of extra higher curvature couplings in the matter sector, may alleviate this situation even for positive cosmological constant.

B. Exact charged black hole solution
The exact solution of the previous section can be charged, by adding to the action principle (12), the Maxwell action Assuming a dyonic ansatz for the electromagnetic field, in the Schwarzschild-like coordinates of (18), one gets where vol(Ω γ ) represents the volume element of the two-dimensional, Euclidean manifold Ω γ in (18), while Q e and Q m stand for the electric and magnetic charge, respectively. For this charged configuration, the metric function f (r) and the scalar field φ(r) read provided the following constraint between the integration constants (µ, Q e , Q m ) is fulfilled As expected, this solution reduces to the one in the previous section when the charges Q e,m vanish. It is interesting to notice that this spacetime exists as a solution of the theory for arbitrary values of the couplings, and the strength of the self-interacting potential is not fixed in terms of the remaining couplings as in the uncharged case. Also notice that while the Weyl coupling α 2 has been set to zero, the remaining higher curvature couplings α 1 and β in (12) do contribute to the solution. Again, reality of the scalar field in (23) imply that the cosmological constant Λ and the algebraic, self-interacting coupling λ must have opposite sign.

V. ON THE ASYMPTOTIC BEHAVIOR
The black holes constructed in the previous sections are dressed by a secondary hair, and are not continuously connected with Schwazschild-(A)dS spacetime, since as we turn off the scalar field in (19), the spacetime approaches the vacuum. The general solution, for arbitrary values of the coupling with primary hair has to be integrated numerically, and here we provide some general properties of such solution, when it asymptotes to a solution of GR, in vacuum. Let us focus on the asymptotically AdS case, and therefore set Λ = −3/l 2 , which via the previous argument implies R = −12/l 2 , globally. Since the contributions to the energy-momentum tensor coming from the higher curvature conformal scalar contain negative powers of the scalar, one has to check whether such terms can indeed be considered as perturbations of the standard conformally coupled scalar in the asymptotic region where the scalar field tends to zero. As a matter of fact this is the case for the particular black holes solutions we found in the previous sections, as well as in higher dimensions d > 4 as can be seen from the explicit black hole solutions constructed in [18], but now in dimension four we have to check whether this is the case in a more generic setup. Since we are interested in the asymptotic behavior induced by the conformally coupled scalar field, we turn-off the purely gravitational higher curvature terms, i.e. we set α 1 = 0 and α 2 = 0 in (12). Let us introduce an effective curvature radius at infinity l eff . In global AdS f (r) = g(r) = r 2 /l 2 eff + 1 and in the absence of higher curvature terms and a self-interacting potential (i.e. λ = β = 0), the equation (14) reduces to that of a free, massive scalar with mass equal to the conformal mass m 2 = −2/l 2 eff , which is above the BF mass in four dimensions m 2 BF = −9/(4l 2 eff ). For this mass, the asymptotic behavior of the scalar field, which will be consistent with the AdS asymptotics in the backreacting case, read Dirichlet boundary conditions are such that a i = 0 while Neumann boundary conditions are defined by setting b i = 0. It is known that for the conformal mass, both sets of boundary conditions lead to a consistent quantization of the scalar in global AdS [43,44]. We want to see whether this asymptotic behavior is consistent in the fully backreacting case, which is not clear at a first sight due to the presence of negative powers of the scalar field in (14) and (13). We now consider in the gravitational field equations, the asymptotic expansion for the scalar field (25) as well as the following asymptotic expansion for the metric functions Let us first analyze the case when both branches in (25) are present, i.e. a 0 , b 0 = 0. For the metric equations (13) one gets the following expansions The contribution for the matter sector (25) appears explicitly in the O(r −4 ) terms of (27) and (28), and the vanishing of these terms, relate the arbitrary constants on the asymptotic expansions of the field and the metric functions. Notice that in all the equations, the leading terms leads to l eff = l, and therefore the curvature radius at infinity is completely controlled by the bared cosmological constant and does not receive contributions from the scalar.
Finally, the sub-leading term in the expansion at infinity of the field equations, relates the constants m f and m g , implying therefore that the spacetimes will not only be asymptotically AdS, but asymptotically Schwazschild-AdS.
When the slow branch is not present in (25), i.e. for Dirichlet boundary conditions for the scalar field (a 0 = 0), it is enough to consider the coefficients b i as non-vanishing. In such case, the asymptotic expansion of the components of the gravitational field equations (13) imply that b 1 has to vanish, and lead to This again implies that the curvature radius at infinity is determined only in terms of the bared cosmological constant. Even more, there is a single integration constant characterizing the scalar field at infinity since b 2 gets fixed in terms of b 0 . The remaining orders for the equations relate the b i>2 with b 0 and gravitational constant appearing in the expansion of the metric functions in (26). These asymptotic behavior are consistent with those of a localized matter distribution in AdS in four dimensions [41].
Is it interesting to notice that the for the exact black holes found in the previous section, the scalar field fulfils Neumann boundary conditions, since there is a single branch at infinity, controlled by the integration constant µ that falls-off with a leading term r −1 .

VI. BEYOND THE QUADRATIC LEVEL
Now, lets show that the neutral solution found in section IV.A, survives the inclusion of a certain type of higher curvature, conformal couplings. Let us consider the action It is simple here to substitute first a metric and scalar field of the form in the action (31) and then consider the field equations of the reduced action [48]. The solution in this case, takes the same form as before provided the following two constraints are fulfilled As before, these equations imply the existence of a single, independent integration constant that can be taken as µ or a. Nevertheless, the inclusion of the higher curvature terms, i.e. β i>2 = 0 allow to have spherically symmetric, asymptotically de Sitter black holes, maintaining the reality conditions of the scalar field. See Figure 1.  (33), but only the former leads to a potential that is bounded from below λ ∼ .043 (we have also set 16πG = 1).

VII. FURTHER COMMENTS
In this paper we have studied conformally coupled scalar fields containing conformal couplings with higher powers in the curvature, in four dimensions. We have shown that the solution found in [14] can be embedded in these theories, provided the relations between the couplings are updated due to the presence of the new terms in the action. In contrast with what occurs in such reference, the presence of the new terms allow to construct asymptotically de Sitter black holes with self-interactions of the form λφ 4 which are bounded from below. These solutions fulfil Neumann boundary conditions, which are still consistent in the backreacting case, even in the presence of the non-standard derivative, conformal, self-interactions of the scalar field.
We proved that even in the presence of quadratic terms in the curvature in the purely gravitational action, as well as a cosmological term, one has that the static, regular black holes must have a constant Ricci scalar. Such argument can be trivially extended to a certain family of theories in higher dimensions as follows. For simplicity we focus on the gravity part, but the argument extends to conformally invariant matter sources. Consider the following theory in dimension D = 4k where α 1,2 are are arbitrary, dimensionless, coupling constants. The field equations are of fourth order and or the form where B µν is a traceless, generalized Bach tensor, that comes from the variation of the Weyl invariant Lagrangian tr C 2k and E µν = 2kR 2k−1 R µν − 1 4k g µν R + 2kg µν R 2k−1 − 2k∇ µ ∇ ν R 2k−1 .
The trace of the field equations (37) reads Assuming the general ansatz for a static spacetime given in (16), leads to Where ∇ i is the covariant derivative constructed with the Levi-Civita connection associated to h ij . Then, after multiplying the trace of the field equations by R 2k−1 N one gets and then we integrating on the spacelike section and assuming that the boundary term vanishes as in Section III, we get Then, assuming α 1 > 0 we have that the vanishing of the integral, requires the vanishing of each of the terms in the integrand, therefore  (18), all these terms contribute with the same functional form to the field equations [49].
It is well know that solutions of GR with a standard conformally coupled scalar field can be embedded in eleven-dimensional supergravity (see e.g. [50]). It would be interesting to see whether the corrections to the conformal coupling we have introduced here can be mapped with some corrections of the eleven-dimensional theory.