Einstein–Maxwell-axion theory: dyon solution with regular electric field

In the framework of the Einstein–Maxwell-axion theory we consider static spherically symmetric solutions which describe a magnetic monopole in the axionic environment. These solutions are interpreted as the solutions for an axionic dyon, the electric charge of which is composite, i.e. in addition to the standard central electric charge it includes an effective electric charge induced by the axion–photon coupling. We focus on the analysis of those solutions which are characterized by the electric field regular at the center. Special attention is paid to the solutions with the electric field that is vanishing at the center, and that has the Coulombian asymptote, and thus displays an extremum at some distant sphere. Constraints on the electric and effective scalar charges of such an object are discussed.


Introduction
In 1987 Wilczek has formulated the idea that for a distant observer the magnetic monopole in an axionic environment looks like a dyon with a magnetic and effective electric charge [1]. This idea was based on the prediction of the axion electrodynamics where the interaction between the radial magnetic field attributed to the monopole and the surrounding pseudoscalar (axion) field produces the radial electric field without real electric charge at the center. That is why it is said that Wilczek in 1987 presented the first example of the so-called axionic dyon. The axion electrodynamics on which this result was based have been established and developed in the decade 1977-1987, being inspired by the theoretical discovery of Peccei and Quinn of the CP-invariance conservation [2] and by discussions about a new light pseudo-Goldstone boson introduced by Weinberg [3] and Wilczek [4]. The model of a e-mail: Alexander.Balakin@kpfu.ru b e-mail: Alexei.Zayats@kpfu.ru the coupling of the pseudoscalar and electromagnetic fields was formulated in a covariant form by Ni in [5]; the axion electrodynamics written in the 3-dimensional form was used by many authors (see, e.g., the work of Sikivie [6]). Since the axions are considered to be candidates to the dark matter particles [7][8][9][10][11][12][13][14][15] the physics of axions had become one of the key elements of numerous applications to cosmology and astrophysics. These applications take into consideration various models of interaction of gravitational, electromagnetic, scalar and pseudoscalar fields which are nowadays called the Einstein-Maxwell-axion and Einstein-Maxwell-axiondilaton models (see, e.g., [16][17][18]). Also, these applications focus the attention on the models which belong to the class of theories associated with extended axion electrodynamics [19][20][21][22][23][24][25][26].
In 1991 Lee and Weinberg [27] studied spherically symmetric solutions for static black holes with a massless axionlike scalar field; in fact it was a realization of the Wilczek idea in the framework of the Einstein-Maxwell-axion theory. Lee and Weinberg have obtained self-consistent master equations for the axion field and metric coefficients, analyzed the asymptotic properties of the solutions and studied the analytic and numeric solutions for the cases of large and small values of the constant of the axion-photon coupling. If we omit the initial electric charge at the center of the object described in [27] we find the solution for the axionic dyon, which was obtained in the framework of the Einstein-Maxwell-axion model and was predicted in [1] using the simple Maxwell-axion model. In this sense it can be said that in [1,27] the authors presented the first (static) example of the so-called longitudinal magneto-electric cluster in which the magnetic and axionically induced electric fields are parallel to one another. Later the solutions describing the Longitudinal Clusters were found in the systems with the pp-wave symmetry [28] and in the context of the search for fingerprints of relic axions in the terrestrial magnetosphere [29]. Now we are interested to find a regular solution for the axionic dyon. What does this mean? In 1968 Bardeen [30] attracted the attention to solutions of the field equations which are regular in the center. The first idea was to modify the equations for the electric field so that it will be finite at r = 0; for instance it might be the function E(r ) = Q r 2 +a 2 with E(0) = Q a 2 and the Coulombian asymptote E(r ) → Q r 2 . In the framework of the Einstein-Maxwell theory the regularity in the center assumes that not only the electric field is finite but the metric coefficients and all curvature invariants are finite as well. The story of the search for regular solutions is worthy to be subject of special review; we would like to mention only three details in this context. First, the nonminimal coupling of electromagnetic and gauge fields can provide the gravitational field to be regular (see, e.g., [31][32][33]). Second, the nonminimal coupling can provide the electric field to be finite in the center (see, e.g., [34][35][36]). Third, for solutions with a magnetic monopole field the situation is not perfect; for the mentioned solutions the first invariant of the electromagnetic (or gauge) field, B 2 − E 2 is not regular in the center since the magnetic field, B(r ) = ν r 2 , in contrast to the electric one, cannot be finite there. As for the second (pseudo)invariant (B · E) it is possible for it to be finite in the center, when the electric field is not only finite, but tends to zero not more slowly than r 2 . On the other hand, if the electric field strength E is finite at the origin but does not vanish, the vector field E has a hedgehog-like singularity. Therefore if we expect to find a solution which is characterized by the electric field being regular in the center in the strict sense of the word, we should require the condition E(0) = 0 to be satisfied.
Thus, searching for the regular axionic dyons we are faced with the problem to find an exact solution for the field equations for which the electric component vanishes both at r → 0 and r → ∞. Below we intend to show that this is possible for a magnetic monopole surrounded by the pseudoscalar (axion) field when the guiding parameters of the model are specifically coupled.
The paper is organized as follows. In Sect. 2 basic details of the Einstein-Maxwell-axion theory are revived and wellknown solutions with the vanishing constant of the axionphoton coupling (γ = 0) are recovered, using the harmonic spacetime coordinates. In Sect. 3 we analyze the solutions with nonvanishing γ ; in Sect. 3.1 we discuss an example of an exact solution for the axionic dyon singular at the center; in Sect. 3.2 we study (analytically) the regular solutions of the axion electrostatics in the background of magnetic monopole; the results of the numerical study are presented in Sect. 3.3. Section 4 contains conclusions.

Basic formalism
The action functional of the Einstein-Maxwell-axion model takes the form Here R is the Ricci scalar; g is the determinant of the metric tensor g ik ; κ is the Einstein constant, F ik is the Maxwell tensor, * Fik denotes its dual tensor, φ stands for the pseudoscalar (axion) field; γ is the constant of the axion-photon coupling; and m a is the axion mass.
The variation of the action (1) with respect to potentials of the electromagnetic field A i , to the axion field φ, to the spacetime metric g ik gives, respectively, the equations of axion electrodynamics the equation for the axion field φ and the equations for the gravitational field Here T (M) ik and T (a) ik are the energy-momentum tensors for the electromagnetic and axion fields, respectively, which are defined as follows: The dual Maxwell tensor satisfies the equation ∇ k * F ik = 0 which is free from information on the axion field.

Static spherically symmetric spacetime
Let us consider a static spherically symmetric spacetime with the metric We use the harmonic coordinate system [37] in which the variable u plays the role of a radial coordinate; the spatial infinity corresponds to u = 0. We assume that the axion field depends on the radial coordinate only, i.e. φ = φ(u). This system is more convenient to analyze scalar field models. However, when we will need to revive the usual spherical coordinate notation we put At the spatial infinity, i.e at u = 0, asymptotic behavior of the spacetime metric is supposed to be Minkowskian. It means that In this paper we focus on the study of configurations with a magnetic monopole located at the center; the Maxwell tensor components are chosen to be equal to where the constant μ relates to the magnetic charge, q(u) is a function to be found. To characterize the electric field it is convenient to introduce also the scalar quantity E defined as follows: In fact, the scalar E is the tetrad component of the Maxwell tensor E = √ −F ut F ut . In these terms the equations of axion electrodynamics (2) reduce to one equation yielding the solution The constant of integration φ 0 is the value of the axion field at the infinity, i.e. φ(u = 0) = φ 0 ; similarly we define Q = q(0). The axion field equation (3) takes now the form Using (13) this equation can be rewritten as follows: There are four nontrivial equations of the gravitational field. For the metric (7) four nonvanishing components of The corresponding four nonvanishing components of the energy-momentum tensor take the form (see (5) and (6)) If we assume (as in [27]) that the axion field is massless, m a = 0, three independent equations for gravity field can be rewritten Clearly the first equation is decoupled from other ones and can be immediately resolved. Indeed, the first integral of (22) is and the solution satisfying the condition (9) takes the form Also one can check directly that, with (26), Eq. (23) is a differential consequence of (24) with (13). Thus, the key subsystem of master equations consists of the following pair of equations: When the quantities β(u) and q(u) are found the axion field and the electric field can be reconstructed, In other words, we have to find two functions β and q which satisfy the key system of equations (27). Since we use the nonstandard coordinate u instead of the radial variable r we would like to comment how the known solutions can be displayed in these terms.

Known solutions in the u-representation with vanishing constant of axion-photon coupling
In order to illustrate a behavior of the metric functions β, ρ and the function q, some examples of well-known spacetimes will be given, for which the axion-photon coupling is supposed to be absent, γ = 0.

Schwarzschild solution
When μ = 0 and q = 0 the first equation from Eq. (27) reduces to the following form: and the solution to it with the condition (9) can be found immediately where C = M 2 . Equation (26) gives Thus, the Schwarzschild metric is obtained in the harmonic coordinates After transformation of the radial coordinate (see (8)) this metric returns to its standard form The constant M plays here the role of the mass. It should be mentioned that the u-coordinate system covers the Schwarzschild spacetime from the spatial infinity (u = 0) till the horizon (u → ∞) only. When u → ∞ the metric component g tt = e −2Mu tends to zero, i.e., r → 2M.

Reissner-Nordström solution
Let the axion field and the function q be constant, φ = φ 0 , γ = 0 and q = Q. Then the second equation from Eq. (27) is an identity and the first one is simplified as The solution to this equation, which satisfies the condition where the value u * can be obtained from the condition In order to clarify the sense of the constant C for the Reissner-Nordström solution, we consider the case u → 0. At the origin the metric function β behaves as and keeping in mind the Schwarzschild solution it is possible to identify the factor in front of u with the mass M, i.e., Thus, the Reissner-Nordström solution in the harmonic coordinate system takes the form We have to compare this solution with the well-known one The solutions (38) can be identified with (39) keeping in mind the number of horizons.
When C = −ν 2 < 0, one obtains If u +u * = π/ν then r → 0, therefore this point u = u * corresponds to the central naked singularity. (iii) Double horizon. When M 2 = κ(Q 2 +μ 2 )/2, i.e., when C = 0, we obtain After transformation of the radial coordinate one can derive the standard form of the metric

Penney and Fisher solutions
When the axion-photon coupling constant γ is equal to zero and m a = 0, the Eq. (14) reduces to φ uu = 0, thus the axion field is linear in the variable u The integration constant P can be indicated as a scalar (axion) "charge". The constant C is now a combination of the charges P, Q and μ Equation (24) gives now and the solution to this equation takes the form Here the modified constantC is of the form For the metric functions β and ρ given by (26) the linear element (7) covers the Penney solution [38] Clearly, when P = 0 the constants C andC coincide and the Penney solution reduces to the Reissner-Nordström one. WhenC = 0, i.e., we recover the "extremal" Penney solution In the particular case that both electric and magnetic charges, Q and μ, vanish, the metric (50) turns into the Fisher metric [39] where the constant can be positive, vanishing or negative depending on the relation between the mass M and the scalar (axion) charge P.

Solutions with nonvanishing constant of the axion-photon coupling, γ = 0
Let us consider the general case for which the axion-photon coupling constant γ does not vanish. We deal now with the key system of equations with the boundary conditions The first condition for β is the requirement that the spacetime is asymptotically Minkowskian; the second one introduces the asymptotic Schwarzschild mass M. The first condition for q means that Q is the asymptotic electric charge. As for the last condition, this appears from the relationship φ = φ 0 + q−Q γ μ and the definition for the axion charge φ u (0) = P. As usual we denote the asymptotic value of the pseudoscalar (axion) field as φ(0) = φ 0 . For this version of the key system of equations the constant C is not arbitrary, it satisfies the condition (45) C = M 2 − κ 2 (μ 2 + Q 2 − P 2 ). Clearly, the key system of Eq. (55) does not depend on u explicitly, we see u only as the argument of β(u) and q(u). This means that particular solutions of the form β = β(q(u)) can be searched and replace the derivative β u by dβ dq q u in the key system yielding the following equation: We will use this consequence in the next subsection to obtain a particular exact solution to the key system.

Exact solution with the singularity at the center
In general the key system of equations admits the numerical study only, which is why we would like to start our discussion with a particular but explicit example of a solution when the constant C is vanishing, C = 0. Then the first equation (55) admits the solution quadratic in q: when the five parameters M, Q, γ , κ, μ satisfy the following three relationships: Since C = 0, the second metric coefficient is of the form ρ(u) = ln u. In order to find the function q(u) we focus on the second equation (55). With the parameters given by (59) and boundary conditions (56) the first integral of that equation is so that its implicit solution is expressed in terms of the Gauss error function erf(x), defined as When |q| = ∞, the first Gauss error function in (61) takes finite value; this means that there exists a finite value u ∞ , for which |q(u ∞ )| = ∞. For instance, when Q is positive, q(u ∞ ) = −∞ and u ∞ can be found as follows: The radial function r (8) also can be presented in terms of Gauss error functions: According to this formula, r (u ∞ ) = 0, we obtain E(u ∞ ) = ∞ and φ(u ∞ ) = ∞. In other words, we deal with central singularity at u = u ∞ . On the other hand, q = 0 when u = u 0 , where which is valid for arbitrary signs of Q and μ. Similarly, we obtain Thus the electric field takes zero value, when u = u 0 and r = r (u 0 ). Since E(u = u 0 ) = E(u = 0) = 0, the function E(u) reaches its extremum at the finite value of the variable u (the type of extremum, minimum or maximum, is predetermined by the sign of the electric charge Q). Typical plots of E(r ) and φ(r ) are presented in Fig. 1.

Preamble: the regular solution to the equation of axion electrostatics in the flat spacetime
In order to simplify further interpretation of solutions, let us first assume that the background spacetime is flat, i.e., β = 0 and ρ = ln u. Then the last equation in (55) reduces to the form and the solution can be obtained that is discovered by Campbell et al. [40] which is regular at the center u = ∞ and satisfies the condition q( Another boundary condition q (0) = γ μP gives the constraint on the axion charge for which this regular solution exists.

Exact solution to the equation of axion electrostatics
We assume now that the background gravitational field is formed by the magnetic monopole without horizons and with the naked singularity at the center. In fact, the background metric relates to the Reissner-Nordström solution with a magnetic charge. This means that C = −ν 2 < 0 and where R μ = κμ 2 2 is the Reissner-Nordström radius. When u + u * → π ν we obtain e β → 0 and r → 0. In this spacetime background the function q(u) which determines the electric field induced by the axion-photon coupling satisfies the equation The replacement z = i cot ν(u +u * ) transforms this equation into the Legendre equation, where the parameter α is introduced as follows: The variable z is complex; the quantity |z| takes the value |z| = M ν at u = 0 and becomes infinite |z| = ∞, when u = π ν − u * . We search for the solution q(z), which is regular for the interval M ν < |z| < ∞, and we especially require that the solution is regular at u = π ν − u * . As usual, q(z), the solution of the Legendre equation (74) is the linear combination of P α (z) and Q α (z), the Legendre functions of the first and second kinds, respectively (see, e.g., [41] for details). Keeping in mind the analytic properties of the Legendre functions, the regular solution for q(u) can be written, satisfying the condition q(0) = Q, in the following form: Here (cot ν(u + u * )) is the Heaviside function; as was shown in [23] such a structure guarantees the regularity of the solution on the real axis of the complex plane z. Using (76) and (72) the electric field E(r ) can be presented as a function of r as follows:

Integral representation of the solution
For the analysis of the regularity of the electric field one can use also the convenient integral representations of the Legendre functions (see [41]) which yields in particular Using these representations, it is shown that where the function Z α (ξ ) is defined as follows: and N is the standard Reissner-Nordström metric coefficient The function Z α (ξ ) satisfies the following relations: Using Eq. (85), we obtain The electric field E(r ) = q r 2 is regular at the center r = 0 when α ≥ 1. The value E(0) is finite when α = 1, and E(0) = 0 when α > 1. The second invariant of the electromagnetic field I (2) ≡ 1 4 F mn * F mn is regular at the center, when α ≥ 3. Indeed, thus, at r → 0 the invariant 1 4 F mn * F mn ∝ r α−3 . In Fig. 2 we present typical plots of the function E(r ) for three values of the parameter α.

Behavior of the axion field
When the function q(u) is found, the axion field φ = φ 0 + q−Q γ μ can easily be reconstructed. In particular one sees that the axion field is regular at the center when the function q(u) takes finite value at r = 0. We focus now on the following detail: when u → 0 the quantity q−Q γ μ tends to uq (0) γ μ . Thus in fact we have to analyze the value of the quantity P Q . This ratio can be calculated as E/E 0 r/R The integral representation (78) of this quantity gives For two limiting cases, M → 0 and M → R μ , this expression takes the form Thus, in this limit, α → ∞, we obtain the result coinciding with the flat spacetime case (see (69)).

Limiting case M → R μ
Let us consider the extremal Reissner-Nordström case with M → R μ . For this limit, we have The metric function β takes the form while the electric field function q according the Eq. (80) can be written as follows: or, excluding the variable r , In contrast to the case M < R μ , the function q vanishes at the double horizon r = R μ = 0 and e β | q=0 → ∞.

Qualitative and numerical studies of the regular solutions
When the spacetime background is not fixed, i.e. the model is self-consistent, we have to solve the general system of the key equations (55) and (57). In contrast to the explicit example demonstrated in Sect. 3.1 regular solutions to this system can be presented in a numerical form only. In this subsection we will study solutions with the electric field, regular at the center, when the function q vanishes at r = 0. The metric function β has to tend to −∞, because the naked singularity associated with the magnetic monopole cannot be removed. Using Eq. (57), we see that β behaves as where the parameter α > −1 has to satisfy the condition Obviously, this relation does not differ from the corresponding expression (75) for the background solution. Equation where the value u 0 corresponds to the value r = 0 at the center. For instance, for the background solution considered above u 0 = π ν − u * . When e ρ does not vanish at u = u 0 , the standard radial coordinate r = e β−ρ behaves as follows: and we have q ∝ r α+1 , β ∝ ln r.
The first formula coincides qualitatively with the corresponding expression for the background solution (see (87) and the electric field E(r ) = q r 2 vanishes at the center when α > 1 as well.
When u = 0 the boundary conditions (56) give If we fix the electric and magnetic charges, Q and μ, the coupling constant γ and the mass M, the desired solution to Eq. (57) with conditions (97) and (100) exists only for a specific value of the axion charge P and the inequality μP/Q < 0 has to be valid. The latter constraint arises from the second equation of (55), because To illustrate the dependence between Q, μ, P and α (or γ ) we present Figs. 3 and 4. Each figure consists of three panels, which correspond to specific values of the coupling parameter α, namely α = 1, 2 and 3, respectively.
On the other hand, if C = 0, then Eq. (57) admits a solution which at q = 0 behaves as follows (cf. (96)): As was mentioned above, such a solution corresponds to the metric, which possesses a double horizon. Curves describing in Figs. 3 and 4 the relationship between the charges P, Q and μ, and the mass M for this limiting configuration are drawn using gray color. Figure 3 depicts dependence between the ratio P/Q and the mass M for fixed values of the electric charge Q and the coupling parameter α. The first (left) curve corresponds to the limiting case Q μ described in Sect. 3.2. Other curves correspond to |Q| = n|μ|, where n = 1, . . . , 5. Figure 4 illustrates dependence between the ratio Qμ and the mass M for fixed values of the axion scalar charge P and the coupling parameter α. The range M/R μ ∈ [0, 1) on the horizontal axis corresponds again to the limiting case P μ described in Sect. 3.2. Color curves correspond to |P| = 0.5n|μ| where n = 1, . . . , 5. The gray line defines the mass-charge relation for the spacetime metric with double (extremal) horizon. If the gravitational interaction is much weaker than the axion-photon coupling, i.e. when γ 2 κ, or, equivalently, α 1, the color curves become horizontal, |P| → |Q| (see (69)). The solution of Campbell, Kaloper and Olive (68) can be considered as a non-gravitational limit.

Conclusions
In the present paper we realize Wilczek's idea about a magnetic monopole surrounded by an axion-induced radial electric field in the framework of the Einstein-Maxwell model with the massless axion field. Since this electric field is created by interaction between the magnetic field of the monopole and the axion field and is not related to any real electric charge, the electric field has to be regular in the center in the strict sense, i.e. E(0) = 0. In this sense our solution is a generalization of the result of Campbell et al. [40], taking into account the gravitational field of the monopole.
In Sect. 3.2 we present the four-parameter family of solutions (see Eqs. (77), (80)) in the framework of the axion electrodynamics on the background of the magnetic monopole gravitational field with the metric of the Reissner-Nordström type. The fifth parameter, the axion field charge P, is determined by other parameters, namely the electric and magnetic charges Q and μ, the mass M and the coupling parameter α (see Eq. (89)). Besides this relation, the parameters are bounded by two inequalities which correspond to the requirements of the absence of horizons (M 2 < κμ 2 /2) and regularity at the origin (α > 1). In addition, when α ≥ 3 the invariant scalar F ik * F ik appears to be regular in the center too.
In Sect. 3.3, using numerical methods, we solve the total system of equations attributed to the Einstein-Maxwellaxion model, in which the gravitational field is self-consistent, not the background one. We demonstrate that the behavior of the solutions to the self-consistent system qualitatively coincides with the background solution. This background solution can be extracted from the general solution as an asymptotic case with Q, P μ.