A class of solitons in Maxwell-scalar and Einstein-Maxwell-scalar models

Recently, no-go theorems for the existence of solitonic solutions in Einstein-Maxwellscalar (EMS) models have been established [1]. Here we discuss how these theorems can be circumvented by a specific class of non-minimal coupling functions between a real, canonical scalar field and the electromagnetic field. When the non-minimal coupling function diverges in a specific way near the location of a point charge, it regularises all physical quantities yielding an everywhere regular, localised lump of energy. Such solutions are possible even in flat spacetime Maxwell-scalar models, wherein the model is fully integrable in the spherical sector, and exact solutions can be obtained, yielding an explicit mechanism to de-singularise the Coulomb field. Considering their gravitational backreaction, the corresponding (numerical) EMS solitons provide a simple example of self-gravitating, localised energy lumps. ‡ herdeiro@ua.pt † jmiguel.oliveira@ua.pt ⋆ eugen.radu@ua.pt


Introduction
Einstein-Maxwell-scalar (EMS) models, described by the action wherein the real scalar field φ has a canonical kinetic term and it is non-minimally coupled to the Maxwell field strength F , via some function f (φ), emerge naturally in physics (R is the Ricci scalar). Well known contexts are Kaluza-Klein models [2][3][4] and supergravity/string theory [5]. In these cases the non-minimal coupling is provided by an exponential function of the sort f (φ) ∼ e −αφ , where α is a constant. But more general classes of coupling functions have been considered, for instance, in the context of cosmology [6,7].
Another interesting class of coupling functions emerged recently in the quest for models that deviate from general relativity only in the strong gravity regime. If f (φ) obeys some simple properties, the EMS model accommodates the phenomenon of spontaneous scalarisation of asymptotically flat charged black holes [8]. This means that even though the electrovacuum Reissner-Nordström (RN) black hole is a solution of the EMS model, for sufficiently high charge to mass ratio this black hole becomes unstable: it becomes energetically favourable for the RN black hole to scalarise. A new family of scalarised black holes bifurcates from the RN family, which contains the end states of this dynamical scalarisation mechanism -see also [9][10][11][12][13][14][15].
The existence of such a class of EMS models, that allow spontaneous scalarisation of charged black holes, and accommodates two distinct families of black holes, the bald RN and the "hairy" scalarised BHs, raised the question if solitons can also be found in EMS models, as this is often the case in theories that accommodate both bald and hairy BHs.
Building on this question, in [1] a set of theorems were established showing that, under various assumptions, no such solitonic solutions exist, similarly to the case of vacuum and electrovacuum. The purpose of this paper is to show that dropping one of the assumptions considered in [1], namely that the coupling function is everywhere finite, it becomes possible to circumvent the aforementioned theorems and obtain solitonic solutions, including in the flat spacetime (Maxwell-scalar theory) limit. In the latter case, moreover, the solitons can 2 Flat spacetime Maxwell-scalar models

A physical motivation
An awkward feature of classical electromagnetism is that the energy E of the Coulomb field of a point charge Q is divergent: Quantum considerations naturally introduce an ultraviolet cut-off to the validity of the classical Coulomb solution, regularising this integral. Quantum Electrodynamics (QED), however, is itself incomplete as a quantum field theory, due to the Landau pole [20]. But it yields the important lesson that the coupling constant g, which determines the strength of the electromagnetic interaction in the Maxwell Lagrangian runs with the energy scale.
Whatever fundamental theory turns out to complete QED, it may admit a covariant effective field theory description that captures the dynamics of the coupling. Then, g would emerge as a spacetime function with some dynamics. In a simple model, g would be a real scalar field with a standard kinetic term. Allowing a more general dynamics, one takes g as being an arbitrary function of the scalar field, keeping the latter with a standard kinetic term. This suggests considering the naive covariant effective field theory where φ is a real scalar field, F = dA is the covariant description of the electromagnetic dynamics and the background is Minkowski spacetime. The function f (φ) specifies the dynamics of the gauge coupling. This model ignores higher order corrections in F , so it is certainly incomplete. Nonetheless one may take the aforementioned reasoning as a motivation to consider this class of simple models. Can the Coulomb field of a point charge be de-singularised in this context?

An integrable model
The naive model (4), which is the decoupling limit of the EMS model (1) wherein back reaction is neglected, is integrable in the spherical sector. Taking the following ansatz for the fields in spherical coordinates in Minkowski spacetime (t, r, θ, ϕ): the Maxwell equations yield a first integral: where Q is interpreted as the electric charge. Using this first integral, the Klein-Gordon equation reads which, introducing the coordinate x ≡ 1/r, yields another first integral It is a simple application of the virial theorem, or a Derrick-type scaling theorem [21], to show that solutions must have E = 0. For instance, this can be seen from the condition [22] that holds for time-independent, finite energy field configuration in Minkowski spacetime, where i, j are spatial indices in Cartesian coordinates. Relation (9) is a simple consequence of energy-momentum conservation and can be interpreted as the balancing of the total stresses in an extended object. There are regions where matter is in tension and regions where it is in compression, for any static balanced soliton. Thus, the problem of finding solutions is reduced to solving, from (8), and then inverting x(φ) → φ(x) → φ(r). Fixing the coupling function f (φ) one can obtain φ(r) and, from (6), the electrostatic potential, both as line integrals. Due to the two first integrals the system is fully integrable.

Everywhere regular solutions
To assess if the solutions have finite energy one must consider the energy-momentum of the model, This yields the energy density ρ, after using (8): and the total energy, E, obtained by integrating the energy density on a spacelike slice Σ In order to obtain regular solutions at the origin we assume the scalar field admits a power series expansion near the origin: We do not constrain the constant coefficient φ 0 , which may or may not vanish. Apart from φ 0 , let φ N , where N ∈ N 1 be the first non-vanishing coefficient in this expansion. Then, from (8), Thus, as r → 0, Regularity of the scalar field at the origin then requires the coupling to diverge as ∼ 1/r 2N +2 . From (12) this implies that the energy density is finite therein and from (6), close to the origin. Thus, all physical quantities are finite close to the origin, for this class of behaviours of the coupling.

A class of examples
There is still, of course, some freedom in choosing the coupling function, within the class with the correct divergent behaviour at the origin. Let us consider examples.

A simple coupling yielding regular solutions
As an explicit example, consider where α is a non-zero constant. Then (10) immediately yields, taking the integration constant such that φ → 0 as r → ∞: One observes that φ(r) is regular and smooth as r → 0, φ(r) ≃ 1/α − r/(Qα 2 ); thus we expect, from (16), that the coupling to diverge as 1/r 4 . Asymptotically , on the other hand, φ(r) ≃ Q/r. Thus the scalar "charge" coincides with the electric charge. Plugging (19) into (18) yields: The coupling diverges as 1/r 4 at the origin, as anticipated. This divergence precisely cancels the divergence of the Maxwell field at the origin, cf. (12), making it finite and non-zero. In fact, the energy density, from (12), is It follows that the total energy (13) is Now, using (6) we obtain for the electrostatic potential: All the quantities (21), (22), (23) manifestly reduce to the usual Coulombic ones upon taking α → 0. In such case (19) reduces to the profile of a scalar charge Q at the origin.
The expressions make manifest how α regularises the solution.

A family of couplings yielding regular solutions
As further examples, with slightly different features, we generalise the coupling (18) as where n is an integer. Using this coupling, equation (10) gives which has a different indefinite integral for n = 2 and n = 2.
For n = 2, imposing φ(r → ∞) = 0 to fix the integration constant, one obtains which reduces to (18) for n = 4. For regular solutions at the origin we require lim r→0 φ to be finite. This implies n > 2, in which case which is finite, as required. For n = 3 we see that the second term goes as r 2 ; but for n > 4, the second term has a non-integer power. In the former case we anticipate, from (16), that the coupling diverges as 1/r 6 . In the latter case, φ is not analytic at the origin. It will, nonetheless yield a regular solution, when analysing the usual physical quantities.
The coupling f (φ) as a function of r then reads: which diverges as ∼ 1/r 2n n−2 at the origin, for n > 2, but respects the condition lim r→∞ f (r) = 1. We confirm, in particular, the 1/r 6 divergence, for n = 3 and a divergence with (generically) a non-integer inverse power for n 5. The electric field E µ = −∂ µ V (r) has only one non-zero component which reads, from (6) which behaves as r 4 n−2 near the origin, and it is thus regular for n > 2. The total energy now reads Thus, the family of cases with n > 2 illustrate how regular solutions can be obtained, with a different analytic behaviour of the scalar field near the origin (the cases n = 3 and n = 4) and non-analytic behaviour (n > 4).
With n = 2, following a similar reasoning one obtains which is a regular solution at r = 0 with lim r→0 φ(r) = 1/α. The coupling function f (φ) becomes f (r) = e 2αQ/r , which, as before, also diverges at r = 0 but respects lim r→∞ f (r) = 1. Observe, however, it does not diverge as an inverse power of r, which was the conclusion in Section 2.3. This is because, again, φ in this case does not admit a power series expansion near the origin.
This illustrates yet a different example of divergent coupling that yields regular solutions.
The electric field is now and the total energy is In these considerations αQ was assumed to be positive. Otherwise the total energy (30)-(34) would be negative, which would violate the weak energy condition. Interestingly enough, despite the seemingly different solution for n = 2, the total energy E is a smooth function of the power n, as (34) coincides with setting n = 2 in (30).

Dilatonic coupling: a spherically symmetric solution
As mentioned in the Introduction, a dilatonic coupling where α is a constant, emerges in relevant scenarios. Let us thus briefly mention the existence of a spherically symmetric, exact solution for this coupling.
Considering (35) in (10), and taking the integration constant so that the scalar field vanishes asymptotically we immediately get Thus, the coupling, as a function of r is Thus, the coupling diverges at the origin and, if αQ > 0 it is regular elsewhere. Moreover, using now (6) we get One finds the following small-r expansion of the solution thus, the scalar field diverges at the origin. Asymptotically, on the other hand, both fields are well behaved The energy density of this solution diverges at the origin: the total mass, however, is finite This solution is interesting in that it shows a divergent coupling can source a finite mass configuration which, nonetheless, is not fully regular, as the scalar field and the energy density diverge at the origin.

The gravitating solitons
The above flat spacetime solutions can be made to self-gravitate by coupling (4) to Einstein's general relativity. For the case of the regular solutions described in the previous section, this yields, perhaps, the simplest models of charged soliton.
One now considers the EMS model (1), where c = G = 1. In addition to the ansatz (5) we consider the metric form and m(r) is the Misner-Sharp mass; r is thus the areal radius, a geometrically meaningful coordinate.
The ansatz (5) and (43) yield the following effective Lagrangian: As in the flat spacetime case, the equation of the electric potential possesses a first integral, which generalises (6), and reads where again the integration constant Q is the electric charge, which we shall assume to be strictly positive, without any loss of generality. Using this integral, the remaining equations of motion become 1 where we have defined .
The smooth of a spacetime configurations can be assessed by analysing the Ricci scalar and the Kretschmann scalar

Near the origin
A small r analysis of the field equations confirms the conclusion observed in the flat space analysis: for a scalar field admitting a power series expansion near the origin φ = φ 0 + φ 1 r + . . . and φ 1 = 0, if the coupling diverges as 1/r 4 , finite energy, everywhere regular solutions are possible. To see this, we again start by assuming the existence of a power series expansion of solutions, with the scalar field approaching a finite nonzero value where φ 0 is arbitrary. Then, the equations of motion, together with the assumption of regularity, impose, for the n th derivative of S(φ) computed at the origin, denoted S (n) (φ 0 ), This implies the advertised behaviour: the coupling function f (φ) diverges as 1/r 4 as r → 0. This behaviour cancels the divergence associated with the presence of an electric charge, providing a smooth configuration as r → 0.
The small r expansion of the matter functions reads while for the metric functions we find where φ 2 and δ 0 are constants that are fixed by the numerics when integrating the field equations from the origin to infinity and requiring the correct asymptotic behaviour. With this expansion, both the Kretschmann curvature scalar and Ricci scalar are finite as r → 0, taking the form and R = 48 The small r expansion of S(φ) reads which implies the following generic approximate form of the coupling function Of course, we could have assumed that in the scalar field expansion φ 1 = 0 and the power series starts at a higher order term. This would impact in the way the coupling diverges at the origin, similarly to the flat spacetime analysis of section 2.3. For concreteness, here we focus on the case with φ 1 = 0.

Near infinity
A large r analysis of the field equations, on the other hand, imposing yields the following approximate solutions: Here M is the ADM mass and Q is the electric charge; Φ is the electrostatic potential at infinity and Q s is the scalar 'charge' which in general needs not equal the electric charge (it did in the flat spacetime illustration above). In fact, the equations of motion possess again two first integrals implying that the gravitating solutions satisfy the following relation This last equation, in particular, shows manifestly the curved background breaks the equality between Q and Q s .
Interestingly, one can show that there is a Smarr relation in terms of these asymptotic quantities, which is not affected by the scalar field, Moreover, a first law of thermodynamics can be obtained in the form We emphasise the absence of a scalar field contribution in these relations.

The full solutions
The gravitating version of the exact solution in Minkowski spacetime described in subsection 2.4.1, with coupling (18), and whose asymptotic limits have been described in subsection 3.1, can be constructed numerically. The set of four ordinary differential equations obtained from the above setup was solved numerically by using a standard Runge-Kutta ordinary differential equation solver and appropriate boundary conditions. Fixing α, gravitating solitons exist for arbitrary large values of Q. The profile of a typical solution is shown in Fig. 1. As one can see, the profiles of the various functions, and in particular that of the Kretschmann scalar K, are smooth as r → 0. In Fig. 2 we show the ADM mass vs. electric charge diagram for families of solutions with different values of α. One can see that for all families, the solutions trivialise as α → 0. Moreover, a smaller α implies that the same value of the electric charge can support a more massive soliton. Obviously, the solutions also trivialise as Q → 0. The electric charge supports the soliton. This is also manifest from the following virial identity that can be derived for these solutions:  For Q = 0 the right hand side vanishes, and so must the left hand side, which implies φ ′ = 0 and hence no non-trivial scalar profile exists.
Self-gravitating solitons with the coupling (24) and n = 3 where also obtained. They follow the same pattern as the n = 3 case, which is therefore illustrative.

Perturbative solutions
The existence of a flat spacetime solution, whose total mass-energy is proportional to 1/α, suggests that the self-gravitating solitons may be expressed as a perturbative series in 1/α. Let us indeed show that the numerical solutions of the previous subsection can be approximated by such perturbative solutions. This approximation, as we will show and as one may anticipate, is accurate for sufficiently large α.
The perturbative solutions are obtained by performing a power series expansion for all relevant functions As for the numerical solutions of the previous subsection, we focus on the quartic coupling function (18). Solving iteratively the field equations order by order in 1/α, we arrive at the following expressions 2 m 1 (r) = 0 = m 3 (r) , m 2 (r) = qr 3 3(q + r) 3 , m 4 (r) = − qr 4 (10q 2 + qr + r 2 ) 90(q + r) 6 , δ 1 (r) = 0 = δ 3 (r) , δ 2 (r) = q 2 (q + 3r) 6(q + r) 3 , where q is a free parameter, whose physical significance becomes transparent by computing the far field asymptotics of the electric potential. One finds it is related to the electric charge measured at infinity Q, as The perturbative solution yields the following ADM mass and scalar charge, valid to fourth order in perturbation theory: Observe that the first terms in (71) reproduce the flat spacetime limit, eq. (22) and the fact that the electric and scalar charge coincide in that limit.
In Fig. 2 the perturbative solutions (71) are compared with the numerical solutions.
One can observe that the former provide a good approximation for large values of α; for instance, for α = 10 the relative difference between the numerical result for M (Q) and the theory one is around 10 −4 . However, the differences start to increase for smaller α. This is illustrated by the results for α = 0.2 in Fig. 2.
Finally let us mention that a similar solution has been derived for the self-gravitating solitons with the coupling (24) and n = 3. In this case one finds, e.g.

Discussion
Recently, a set of theorems were shown [1] establishing that the model (1) does not allow self-gravitating solitons. One of the observations therein is that if the coupling would diverge the theorems could, potentially, be circumvented. The purpose of this paper was to provide the mechanism how this can happen providing a simple construction of flat spacetime and gravitating solitons.
Preliminary analysis shows the solitons we have described herein are stable against spherical perturbations. If this is the case for generic perturbations, these solitons can be used for dynamical studies in many setups, as, for instance, boson stars [23]. Moreover, this construction reveals how to de-singularise the Coulomb field in a classical field theory, without resorting to non-linear electrodynamics, as in the Born-Infeld model [24], or invoking a manifestly extended object, such as in the Dirac model of the electron as a spherical membrane [25].
Finally, let us remark that there is a well known similarity between the EMS model and the extended scalar-tensor-Gauss-Bonnet model, where black hole scalarisation was first pointed out in [26][27][28]. Very recently, a family of particle-like solutions in the latter model were discussed [29]. These particle like solutions are also supported by a divergent coupling making them the counterparts of the solutions described herein. But in the cases reported in [29] the scalar field also diverges at the origin, in contrast with our fully regular solutions.