Einstein-Maxwell-axion theory: Dyon solution with regular electric field

TL;DR

This paper explores static, spherically symmetric solutions in Einstein-Maxwell-axion theory, revealing axionic dyons with regular electric fields, including those with vanishing central electric fields and Coulombian asymptotes.

Contribution

It introduces and analyzes new solutions describing axionic dyons with regular electric fields, highlighting the effects of axion-photon coupling on electric charge structure.

Findings

Solutions with regular electric fields at the center are found.

Electric charge can be effectively induced by axion-photon interaction.

Constraints on electric and scalar charges are discussed.

Abstract

In the framework of 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 that solutions, which are characterized by the electric field regular at the center. Special attention is paid to the solutions with the electric field, which is vanishing at the center, has the Coulombian asymptote and thus display an extremum at some distant sphere. Constraints on the electric and effective scalar charges of such an object are discussed.

Equations

$$S_{EMa}=\int \mathcal{L}\sqrt{-g}{d}^{4}x,$$

$$\mathcal{L}=\frac{R}{2\kappa }+\frac{1}{4}{F}_{ik}{F}^{ik}+\frac{1}{4}\gamma F_{ik}\stackrel{\ast }{F}{}^{ik}\varphi$$

$$-\frac{1}{2}{\nabla }_i\varphi {\nabla }^i\varphi +\frac{1}{2}{m}_a^2{\varphi }^2.$$

$${\nabla }_i{F}^{ik}+\gamma \stackrel{\ast }{F}{}^{ik}{\nabla }_i\varphi =0,$$

$${\nabla }_i{\nabla }^i\varphi +m_a^2\varphi +\frac{\gamma }{4}{F}_{ik}\stackrel{\ast }{F}{}^{ik}=0,$$

$$R_{ik}-\frac{1}{2}R{g}_{ik}=\kappa \left(T_{ik}^{(M)}+T_{ik}^{(a)}\right).$$

$$T_{ik}^{(M)}=\frac{1}{4}{g}_{ik}{F}_{mn}{F}^{mn}-F_{in}{F_k}^n,$$

$$T_{ik}^{(a)}={\nabla }_i\varphi {\nabla }_k\varphi -\frac{1}{2}{g}_{ik}{\nabla }_n\varphi {\nabla }^n\varphi +\frac{1}{2}{m}_a^2{\varphi }^2{g}_{ik}.$$

$$d{s}^2={\mathrm{e}}^{-2\beta (u)}d{t}^2-{\mathrm{e}}^{2\beta (u)}\left[{\mathrm{e}}^{-4\rho (u)}d{u}^2+{\mathrm{e}}^{-2\rho (u)}d{\Omega }^2\right].$$

$$r={\mathrm{e}}^{\beta (u)-\rho (u)}.$$

$$\beta (0)=0,\rho (u){|}_{u\to 0}\sim \ln u.$$

$$F_{ut}=q(u){\mathrm{e}}^{-2\beta (u)},F_{\theta \phi }=\mu \sin \theta ,$$

$$E\equiv \frac{q}{r^2}=q{\mathrm{e}}^{2\rho -2\beta }.$$

$${\left(q-\gamma \mu \varphi \right)}_u^{\prime }=0,$$

$$q=Q+\gamma \mu \left(\varphi -{\varphi }_0\right).$$

$${\varphi }_{uu}^{\prime \prime }-\gamma \mu q{\mathrm{e}}^{-2\beta }-m_a^2\varphi {\mathrm{e}}^{2\beta -4\rho }=0.$$

$$q_{uu}^{\prime \prime }-{\gamma }^2{\mu }^{2}q{\mathrm{e}}^{-2\beta }-m_a^2\left(q+\gamma \mu {\varphi }_0-Q\right){\mathrm{e}}^{2\beta -4\rho }=0.$$

$$G_u^u={\mathrm{e}}^{-2\beta +4\rho }\left({{\beta }_u^{\prime }}^2-{{\rho }_u^{\prime }}^2+{\mathrm{e}}^{-2\rho }\right),$$

$$G_{\theta }^{\theta }=G_{\phi }^{\phi }={\mathrm{e}}^{-2\beta +4\rho }\left(-{{\beta }_u^{\prime }}^2+{{\rho }_u^{\prime }}^2+{\rho }_{uu}^{\prime \prime }\right),$$

$$G_t^t={\mathrm{e}}^{-2\beta +4\rho }\left(-{{\beta }_u^{\prime }}^2+{{\rho }_u^{\prime }}^2-2{\beta }_{uu}^{\prime \prime }+2{\rho }_{uu}^{\prime \prime }+{\mathrm{e}}^{-2\rho }\right).$$

$$T_u^u=\frac{1}{2}{\mathrm{e}}^{-2\beta +4\rho }\left[({\mu }^2+q^2){\mathrm{e}}^{-2\beta }-{{\varphi }_u^{\prime }}^2\right]+\frac{1}{2}{m}_a^2{\varphi }^2,$$

$$T_{\theta }^{\theta }=T_{\phi }^{\phi }=\frac{1}{2}{\mathrm{e}}^{-2\beta +4\rho }\left[-({\mu }^2+q^2){\mathrm{e}}^{-2\beta }+{{\varphi }_u^{\prime }}^2\right]+\frac{1}{2}{m}_a^2{\varphi }^2,$$

$$T_t^t=\frac{1}{2}{\mathrm{e}}^{-2\beta +4\rho }\left[({\mu }^2+q^2){\mathrm{e}}^{-2\beta }+{{\varphi }_u^{\prime }}^2\right]+\frac{1}{2}{m}_a^2{\varphi }^2.$$

$${\rho }_{uu}^{\prime \prime }+{\mathrm{e}}^{-2\rho }=0,$$

$${\beta }_{uu}^{\prime \prime }+\frac{\kappa }{2}\left({\mu }^2+q^2\right){\mathrm{e}}^{-2\beta }=0,$$

$${{\rho }_u^{\prime }}^2-{\mathrm{e}}^{-2\rho }={{\beta }_u^{\prime }}^2+\frac{\kappa }{2}{{\varphi }_u^{\prime }}^2-\frac{\kappa }{2}\left({\mu }^2+q^2\right){\mathrm{e}}^{-2\beta }.$$

$${{\rho }_u^{\prime }}^2-{\mathrm{e}}^{-2\rho }=C,$$

$$\rho =\ln \Pi (C,u),$$

$$\displaystyle\Pi(C,u)\equiv\left\{\begin{array}[]{l}\dfrac{\sinh\nu u}{\nu},\ \hbox{when $C=\nu^{2}>0$},\par\\ u,\ \hbox{when $C=0$},\par\\ \dfrac{\sin\nu u}{\nu},\ \hbox{when $C=-\nu^{2}<0$}.\end{array}\right.$$

$${{\beta }_u^{\prime }}^2+\frac{\kappa }{2{\gamma }^2{\mu }^2}{q_u^{\prime }}^2-\frac{\kappa }{2}\left({\mu }^2+q^2\right){\mathrm{e}}^{-2\beta }=C,$$

$$q_{uu}^{\prime \prime }-{\gamma }^2{\mu }^{2}q{\mathrm{e}}^{-2\beta }=0.$$

$$\varphi ={\varphi }_0+\frac{q-Q}{\gamma \mu },E(u)=q(u){\mathrm{e}}^{-2\beta }{\Pi }^2(C,u).$$

$${{\beta }_u^{\prime }}^2=C>0$$

$$\beta =Mu,$$

$$\rho =\ln \left(\frac{\sinh Mu}{M}\right).$$

$$d{s}^2={\mathrm{e}}^{-2Mu}d{t}^2-\frac{M^4{\mathrm{e}}^{2Mu}}{{\sinh }^{4}Mu}d{u}^2-\frac{M^2{\mathrm{e}}^{2Mu}}{{\sinh }^{2}Mu}d{\Omega }^2.$$

Full text

I 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 magnetic and effective electric charge Wilczek2 . This idea was based on the prediction of the axion electrodynamics that 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 one can say, that Wilczek presented in 1987 the first example of the so-called axionic dyon. The axion electrodynamics, on which this result was based, has been established and developed in the decade 1977-1987, being inspired by the theoretical discovery of Peccei and Quinn of the CP-invariance conservation PQ , and by discussions about a new light pseudo-Goldstone boson introduced by Weinberg Weinberg and Wilczek Wilczek . The model of coupling of the pseudoscalar and electromagnetic fields was formulated in covariant form by Ni in Ni77 ; the axion electrodynamics written in the 3-dimensional form was used by many authors (see, e.g., the work of Sikivie Sikivie83 ). Since the axions are considered to be candidates to the dark matter particles ADM1 ; ADM2 ; ADM3 ; ADM4 ; ADM5 ; ADM6 ; ADM7 ; ADM8 ; ADM9 , the physics of axions had become one of the key elements of numerous applications to cosmology and astrophysics. These applications involve into consideration various models of interaction of gravitational, electromagnetic, scalar and pseudoscalar fields, which are called nowadays the Einstein-Maxwell-axion, and Einstein-Maxwell-axion-dilaton models (see, e.g., EMAD1 ; EMAD2 ; EMAD3 ). Also, these applications attract the attention to the models, which belong to the class of theories indicated by term Extended Axion Electrodynamics ExtendedAE1 ; ExtendedAE2 ; ExtendedAE3 ; ExtendedAE4 ; ExtendedAE5 ; ExtendedAE6 ; ExtendedAE7 ; ExtendedAE8 .

In 1991 Lee and Weinberg LW1991 studied spherically symmetric solutions for static black holes with massless axionlike scalar field; in fact, it was a realization of the Wilczek idea in the framework of 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 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 LW1991 , we find the solution for the axionic dyon, which was obtained in the framework of the Einstein-Maxwell-axion model, and was predicted in Wilczek2 using the simple Maxwell-axion model. In this sense, one can say, that in Wilczek2 ; LW1991 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 pp1 , and in the context of search for fingerprints of relic axions in the terrestrial magnetosphere BG .

Now we are interested to find a regular solution for the axionic dyon. What does this means? In 1968 Bardeen Bardeen attracted the attention to the solutions of the field equations 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)=\frac{Q}{r^{2}+a^{2}}$ with $E(0)=\frac{Q}{a^{2}}$ and the Coulombian asymptote $E(r)\to\frac{Q}{r^{2}}$. In the framework of 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 search for regular solutions is worthy to be a 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., Reg1 ; Reg2 ; Reg3 ). Second, the nonminimal coupling can provide the electric field to be finite in the center (see, e.g., Reg4 ; Reg5 ; Reg6 ). 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, ${\vec{B}}^{2}-{\vec{E}}^{2}$, is not regular in the center, since the magnetic field, $B(r)=\frac{\nu}{r^{2}}$, in contrast to the electric one, cannot be finite there. As for the second (pseudo)invariant $(\vec{B}\cdot\vec{E})$, it is possible to be finite in the center, when the electric field is not only finite, but tends to zero not slowly than $r^{2}$. On the other hand, if the electric field strength $E$ is finite at the origin, but does not vanish, vector field $\vec{E}$ has a hedgehog-like singularity. Therefore if we expect to find a solution, which is characterized by the electric field 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 of the field equations, for which the electric component vanishes both at $r\to 0$ and $r\to\infty$. Below we intend to show that it is possible for 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 Section II, we remind basic details of the Einstein-Maxwell-axion theory, and recover well-known solutions with vanishing constant of the axion-photon coupling ($\gamma=0$), using the harmonic spacetime coordinates. In Section III, we analyze the solutions with nonvanishing $\gamma$; in Subsection III.1 we discuss an example of exact solution for the axionic dyon singular at the center; in Subsection III.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 Subsection III.3. Section IV contains conclusions.

II Einstein-Maxwell-axion model

II.1 Basic formalism

The action functional of the Einstein-Maxwell-axion model takes the form

[TABLE]

Here $R$ is the Ricci scalar; $g$ is the determinant of the metric tensor $g_{ik}$; $\kappa$ is the Einstein constant, $F_{ik}$ is the Maxwell tensor, ${\mathop{\rule{0.0pt}{0.0pt}{F}}\limits^{\,*}}\rule{0.0pt}{0.0pt}_{ik}$ denotes its dual tensor, $\phi$ stands for the pseudoscalar (axion) field; $\gamma$ 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 $\phi$, to the spacetime metric $g^{ik}$ gives, respectively, the equations of axion electrodynamics

[TABLE]

the equation for the axion field $\phi$

[TABLE]

and the equations for the gravitational field

[TABLE]

Here $T_{ik}^{(M)}$ and $T_{ik}^{(a)}$ are the energy-momentum tensors for the electromagnetic and axion fields, respectively, which are defined as follows

[TABLE]

The dual Maxwell tensor satisfies the equation $\nabla_{k}{\mathop{\rule{0.0pt}{0.0pt}{F}}\limits^{\,*}}\rule{0.0pt}{0.0pt}^{ik}=0$, which is free of information about the axion field.

II.2 Static spherically symmetric spacetime

Let us consider a static spherically symmetric spacetime with the metric

[TABLE]

We use the harmonic coordinate system Bron73 , 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., $\phi=\phi(u)$. This system is more convenient to analyze scalar field models, but when we will need to revive the usual spherical coordinate notation, we put

[TABLE]

At the spatial infinity, i.e, at $u=0$, asymptotic behavior of the spacetime metric is supposed to be Minkowskian. It means that

[TABLE]

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

[TABLE]

where the constant $\mu$ 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

[TABLE]

In fact, the scalar $E$ is the tetrad component of the Maxwell tensor $E=\sqrt{-F_{ut}F^{ut}}$. In these terms the equations of axion electrodynamics (2) reduce to one equation

[TABLE]

yielding the solution

[TABLE]

The constant of integration $\phi_{0}$ is the value of the axion field at the infinity, i.e., $\phi(u=0)=\phi_{0}$; similarly, we define $Q=q(0)$.

The axion field equation (3) takes now the form

[TABLE]

Using (13) we can rewrite this equation as follows:

[TABLE]

There are four nontrivial equations of the gravitational field. For the metric (7) four nonvanishing components of the Einstein tensor $G^{k}_{i}=R^{k}_{i}-\frac{1}{2}\delta_{i}^{k}R$ are

[TABLE]

The corresponding four nonvanishing components of the energy-momentum tensor $T^{k}_{i}=T^{k(M)}_{i}+T^{k(a)}_{i}$ take the form (see (5) and (6))

[TABLE]

If we assume (as in LW1991 ) that the axion field is massless, $m_{a}=0$, three independent equations for gravity field can be rewritten as

[TABLE]

Clearly, the first equation is decoupled from other ones, and can be immediately resolved. Indeed, the first integral of (22) is

[TABLE]

and the the solution satisfying the condition (9) takes the form

[TABLE]

Also, one can check directly, that with (29) Eq. (23) is a differential consequence of (24) with (13). Thus the key subsystem of master equations consists of the following pair of equations

[TABLE]

When the quantities $\beta(u)$ and $q(u)$ are found, the axion field and the electric field can be reconstructed as

[TABLE]

In other words, we have to find two functions $\beta$ and $q$, which satisfy the key system of equations (30). Since we use nonstandard coordinate $u$ instead of radial variable $r$, we would like to comment how the known solutions can be displayed in these terms.

II.3 Known solutions in the $u$-representation with vanishing constant of axion-photon coupling

In order to illustrate a behavior of the metric functions $\beta$, $\rho$ and the function $q$, let us give some examples of well-known spacetimes, for which the axion-photon coupling is supposed to be absent, $\gamma=0$.

II.3.1 Schwarzschild solution

In case when $\mu=0$ and $q=0$, the first equation from Eq. (30) reduces to the following form

[TABLE]

and the solution to it with the condition (9) can be found immediately

[TABLE]

where $C=M^{2}$. The formula (29) gives

[TABLE]

Thus, we obtain the Schwarzschild metric in the harmonic coordinates

[TABLE]

After transformation of the radial coordinate (see (8))

[TABLE]

this metric returns to its standard form

[TABLE]

The constant $M$ plays here the role of the mass. Mention should be made, that the $u$-coordinate system covers the Schwarzschild spacetime from the spatial infinity ($u=0$) till the horizon ($u\to\infty$) only. When $u\to\infty$ the metric component $g_{tt}={\rm e}^{-2Mu}$ tends to zero, i.e., $r\to 2M$.

II.3.2 Reissner-Nordström solution

Let the axion field and the function $q$ be constant, $\phi=\phi_{0}$, $\gamma=0$ and $q=Q$. Then the second equation from Eq. (30) is an identity, and the first one is simplified as

[TABLE]

The solution to this equation, which satisfies the condition $\beta(0)=0$, is

[TABLE]

where the value $u_{*}$ can be obtained from the condition

[TABLE]

In order to clarify the sense of the constant $C$ for the Reissner-Nordström solution, we consider the case $u\to 0$. At the origin the metric function $\beta$ behaves as

[TABLE]

and keeping in mind the Schwarzschild solution we can identify the factor in front of $u$ with the mass $M$, i.e.,

[TABLE]

Thus, the Reissner-Nordström solution in the harmonic coordinate system takes the form

[TABLE]

We have to compare this solution with the well-known one

[TABLE]

The solutions (41) can be identified with (42) keeping in mind the number of horizons.

(i) One horizon.

When $M^{2}>\kappa(Q^{2}+\mu^{2})/2$, i.e., when $C=\nu^{2}>0$, there is one horizon at $u=\infty$ and $r(u\to\infty)=M+\sqrt{C}$.

(ii) Naked singularity.

When $C=-\nu^{2}<0$, one obtains that

[TABLE]

[TABLE]

[TABLE]

If $u+u_{*}=\pi/\nu$ then $r\to 0$, therefore this point $u=u_{*}$ corresponds to the central naked singularity.

(iii) Double horizon. When $M^{2}=\kappa(Q^{2}+\mu^{2})/2$, i.e., when $C=0$, we obtain

[TABLE]

After transformation of the radial coordinate

[TABLE]

one can derive the standard form of the metric

[TABLE]

II.3.3 Penney and Fisher solutions

When the axion-photon coupling constant $\gamma$ is equal to zero and $m_{a}=0$, the Eq. (14) reduces to $\phi^{\prime\prime}_{uu}=0$, thus the axion field is linear in the variable $u$

[TABLE]

The integration constant $P$ can be indicated as an scalar (axion) “charge”. The constant $C$ is now a combination of the charges $P$, $Q$, and $\mu$

[TABLE]

The equations (24) give now

[TABLE]

and the solution to this equation takes the form

[TABLE]

Here the modified constant $\tilde{C}$ is of the form

[TABLE]

For the metric functions $\beta$ and $\rho$ given by (29)) the linear element (7) covers the Penney solution Penney

[TABLE]

Clearly, when $P=0$ the constants $C$ and $\tilde{C}$ coincide, and the Penney solution reduces to the Reissner-Nordström one.

When $\tilde{C}=0$, i.e.,

[TABLE]

we recover the “extremal” Penney solution

[TABLE]

In the particular case, if both electric and magnetic charges, $Q$ and $\mu$, vanish, the metric (53) turns into the Fisher metric Fisher

[TABLE]

where the constant

[TABLE]

can be positive, vanishing or negative depending on the relation between the mass $M$ and the scalar (axion) charge $P$.

III Solutions with nonvanishing constant of the axion-photon coupling, $\gamma\neq 0$

Let us consider the general case, for which the axion-photon coupling constant $\gamma$ does not vanish. We deal now with the key system of equations

[TABLE]

with the the boundary conditions

[TABLE]

The first condition for $\beta$ 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, it appears from the relationship $\phi=\phi_{0}+\frac{q-Q}{\gamma\mu}$ and the definition for the axion charge $\phi^{\prime}_{u}(0)=P$. As usual, we denote the asymptotic value of the pseudoscalar (axion) field as $\phi(0)=\phi_{0}$. and for this version of the key system of equations the constant $C$ is not arbitrary, it satisfies the condition (48) $C=M^{2}-\frac{\kappa}{2}\left(\mu^{2}+Q^{2}-P^{2}\right)$.

Clearly, the key system of equations (58) does not depend on $u$ explicitly, we see $u$ only as the argument of $\beta(u)$ and $q(u)$. This means that we can search for particular solutions of the form $\beta=\beta(q(u))$, and replace the derivative $\beta^{\prime}_{u}$ by $\frac{d\beta}{dq}q^{\prime}_{u}$ in the key system yielding the following equation:

[TABLE]

We will use this consequence in the next subsection to obtain particular exact solution to the key system.

III.1 Exact solution with the singularity at the center

In general case the key system of equations admits the numerical study only, that 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 (58) admits the solution quadratic in $q$:

[TABLE]

when five parameters $M$, $Q$, $\gamma$, $\kappa$, $\mu$ satisfy the following three relationships

[TABLE]

Since $C=0$, the second metric coefficient is of the form $\rho(u)=\ln u$. In order to find the function $q(u)$ we focus on the second equation (58). With the parameters given by (62) and boundary conditions (59) the first integral of that equation is

[TABLE]

so that its implicit solution

[TABLE]

is expressed in terms of the Gauss error function $\mathop{\rm erf}\nolimits(x)$, defined as

[TABLE]

When $|q|=\infty$, the first Gauss error function in (64) takes finite value; this means that there exists a finite value $u_{\infty}$, for which $|q(u_{\infty})|=\infty$. For instance, when $Q$ is positive, $q(u_{\infty})=-\infty$, and $u_{\infty}$ can be found as follows:

[TABLE]

The radial function $r$ (8) also can be presented in terms of Gauss error functions:

[TABLE]

According to this formula, $r(u_{\infty})=0$, thus, we obtain that $E(u_{\infty})=\infty$ and $\phi(u_{\infty})=\infty$. In other words, we deal with central singularity at $u=u_{\infty}$. On the other hand, $q=0$, when $u=u_{0}$, where

[TABLE]

which is valid for arbitrary signs of $Q$ and $\mu$. Similarly, we obtain

[TABLE]

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 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 $\phi(r)$ are presented in Fig. 1.

III.2 Axion electrostatics in the background field of the magnetic monopole

III.2.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 assume, first, that the background spacetime is flat, i.e., $\beta=0$ and $\rho=\ln u$. Then the last equation in (58) reduces to the form

[TABLE]

and we can obtain the solution discovered by Campbell, Kaloper, and Olive Campbell , which is regular at the center $u=\infty$ and satisfies the condition $q(u=0)=Q$

[TABLE]

Another boundary condition $q^{\prime}(0)=\gamma\mu P$ gives the constraint on the axion charge

[TABLE]

for which this regular solution exists.

III.2.2 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=-\nu^{2}<0$, and

[TABLE]

[TABLE]

[TABLE]

where $R_{\mu}=\sqrt{\frac{\kappa\mu^{2}}{2}}$ is the Reissner-Nordström radius. When $u+u_{*}\to\frac{\pi}{\nu}$, we obtain that ${\rm e}^{\beta}\to 0$ and $r\to 0$.

In this spacetime background the function $q(u)$, which determines the electric field induced by the axion-photon coupling, satisfies the equation

[TABLE]

The replacement $z=i\cot\nu(u+u_{*})$ transforms this equation into the Legendre equation

[TABLE]

where the parameter $\alpha$ is introduced as follows:

[TABLE]

The variable $z$ is complex; the quantity $|z|$ takes the value $|z|=\frac{M}{\nu}$ at $u=0$, and becomes infinite $|z|=\infty$, when $u=\frac{\pi}{\nu}-u_{*}$. We search for the solution $q(z)$, which is regular for the interval $\frac{M}{\nu}<|z|<\infty$, and we require especially, that the solution is regular at $u=\frac{\pi}{\nu}-u_{*}$. As usual, $q(z)$, the solution of the Legendre equation (77), is the linear combination of ${\sf P}_{\alpha}(z)$ and ${\sf Q}_{\alpha}(z)$, the Legendre functions of the first and second kinds, respectively (see, e.g., Erd for details). Keeping in mind the analytic properties of the Legendre functions, we can write the regular solution for $q(u)$ satisfying the condition $q(0)=Q$ in the following form:

[TABLE]

Here $\Theta(\cot\nu(u+u_{*}))$ is the Heaviside function; as was shown in ExtendedAE5 such a structure guarantees the regularity of the solution on the real axis of the complex plane $z$. Using (79) and (75) we can present the electric field $E(r)$ as a function of $r$ as follows:

[TABLE]

III.2.3 Integral representation of the solution

For the analysis of regularity of the electric field one can use also the convenient integral representations of the Legendre functions (see Erd )

[TABLE]

which yields, in particular,

[TABLE]

Using these representations, one can show that

[TABLE]

where the function $Z_{\alpha}(\xi)$ is defined as follows

[TABLE]

and $N$ is the standard Reissner-Nordström metric coefficient

[TABLE]

The function $Z_{\alpha}(\xi)$ satisfies the following relations:

[TABLE]

Using the formula (88), one can obtain that

[TABLE]

The electric field $E(r)=\frac{q}{r^{2}}$ is regular at the center $r=0$, when $\alpha\geq 1$. The value $E(0)$ is finite, when $\alpha=1$, and $E(0)=0$, when $\alpha>1$. The second invariant of the electromagnetic field $I_{(2)}\equiv\frac{1}{4}F_{mn}{\mathop{\rule{0.0pt}{0.0pt}{F}}\limits^{\,*}}\rule{0.0pt}{0.0pt}^{mn}$ is regular at the center, when $\alpha\geq 3$. Indeed,

[TABLE]

thus, at $r\to 0$ the invariant $\frac{1}{4}F_{mn}{\mathop{\rule{0.0pt}{0.0pt}{F}}\limits^{\,*}}\rule{0.0pt}{0.0pt}^{mn}\propto r^{\alpha-3}$. In Fig. 2 we present typical plots of the function $E(r)$ for three values of the parameter $\alpha$.

III.2.4 Behavior of the axion field

When the function $q(u)$ is found, the axion field $\phi=\phi_{0}+\frac{q-Q}{\gamma\mu}$ can be easily reconstructed. In particular, one can see 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\to 0$, the quantity $\frac{q-Q}{\gamma\mu}$ tends to $\frac{uq^{\prime}(0)}{\gamma\mu}$, so, in fact, we have to analyze the value of the quantity $\frac{P}{Q}$. This ratio can be calculated as

[TABLE]

Integral representation (81) of this quantity gives

[TABLE]

For two limiting cases, $M\to 0$ and $M\to R_{\mu}$, this expression takes the form

[TABLE]

[TABLE]

Using the formula (89), one can demonstrate the following detail: if $\alpha\to\infty$, i.e., if the axion-photon coupling constant $\gamma$ is much greater than $\sqrt{\kappa}$, the ratio $P/Q$ tends to a constant for any values of $M\in[0,R_{\mu})$:

[TABLE]

Thus, in this limit, $\alpha\to\infty$, we obtain the result coinciding with the flat spacetime case (see (72)).

III.2.5 Limiting case $M\to R_{\mu}$

Let us consider the extremal Reissner-Nordström case with $M\to R_{\mu}$. For this limit, we have

[TABLE]

The metric function $\beta$ takes the form

[TABLE]

while the electric field function $q$ according the formula (83) can be written as follows

[TABLE]

or, excluding the variable $r$,

[TABLE]

In contrast to the case $M<R_{\mu}$, the function $q$ vanishes at the double horizon $r=R_{\mu}\neq 0$ and ${\rm e}^{\beta}|_{q=0}\to\infty$.

III.3 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 (58) and (60). In contrast to the explicit example demonstrated in Subsection III.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 $\beta$ has to tend to $-\infty$, because the naked singularity associated with the magnetic monopole cannot be removed. Using Eq. (60), we obtain that $\beta$ behaves as

[TABLE]

where the parameter $\alpha>-1$ has to satisfy the condition

[TABLE]

Obviously, this relation does not differ from the corresponding expression (78) for the background solution. The formula (100) relates to the following asymptotic behavior of the functions $q(u)$ and $\beta(u)$:

[TABLE]

where the value $u_{0}$ corresponds to the value $r=0$ at the center. For instance, for the background solution considered above $u_{0}=\frac{\pi}{\nu}-u_{*}$. When ${\rm e}^{\rho}$ does not vanish at $u=u_{0}$, the standard radial coordinate $r={\rm e}^{\beta-\rho}$ behaves as follows:

[TABLE]

and we have

[TABLE]

The first formula coincides qualitatively with the corresponding expression for the background solution (see (90), and the electric field $E(r)=\frac{q}{r^{2}}$ vanishes at the center when $\alpha>1$ as well.

When $u=0$ the boundary conditions (59) give

[TABLE]

If we fix the electric and magnetic charges, $Q$ and $\mu$, the coupling constant $\gamma$, and the mass $M$, desired solution to Eq. (60) with conditions (100) and (103) exists only for a specific value of the axion charge $P$, and the inequality $\mu P/Q<0$ has to be valid. This latter constraint arises from the second equation of (58), because

[TABLE]

To illustrate dependence between $Q$, $\mu$, $P$, and $\alpha$ (or $\gamma$), we present Figs. 3 and 4. Each figure consists of three panels, which correspond to specific values of the coupling parameter $\alpha$, namely, for $\alpha=1$, 2, and 3, respectively.

On the other hand, if $C=0$ Eq. (60) admits a solution, which at $q=0$ behaves as follows (cf. (99))

[TABLE]

As it 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 $\mu$, and the mass $M$ for this limiting configuration are drawn using gray color.

Fig. 3 depicts dependence between the ratio $P/Q$ and the mass $M$ for fixed values of the electric charge $Q$ and the coupling parameter $\alpha$. The first (left) curve corresponds to the limiting case $Q\ll\mu$ described in Subsection III.2. Other curves correspond to $|Q|=n|\mu|$, where $n=1,\dots,5$.

Fig. 4 illustrates dependence between the ratio $Q\mu$ and the mass $M$ for fixed values of the axion scalar charge $P$ and the coupling parameter $\alpha$. The range $M/R_{\mu}\in[0,1)$ on the horizontal axis corresponds again to the limiting case $P\ll\mu$ described in Subsection III.2. Color curves correspond to $|P|=0.5n|\mu|$, where $n=1,\dots,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 $\gamma^{2}\gg\kappa$, or, equivalently, $\alpha\gg 1$, the color curves become horizontal, $|P|\to|Q|$ (see (72)). The solution of Campbell, Kaloper, and Olive (71) can be considered as a non-gravitational limit.

IV 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, Kaloper, and Olive Campbell , taking into account the gravitational field of the monopole.

In Subsection III.2 we present the four-parameter family of solutions (see Eqs. (80) and (83)) 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 $\mu$, the mass $M$, and the coupling parameter $\alpha$ (see Eq. (92)). Besides this relation, the parameters are bounded by two inequalities, which correspond to requirements of absence of horizons ($M^{2}<\kappa\mu^{2}/2$) and regularity at the origin ($\alpha>1$). In addition, when $\alpha\geq 3$ the invariant scalar $F_{ik}{\mathop{\rule{0.0pt}{0.0pt}{F}}\limits^{\,*}}\rule{0.0pt}{0.0pt}^{ik}$ appears to be regular in the center too.

In Subsection III.3, using numerical methods, we solve the total system of equations attributed to the Einstein-Maxwell-axion 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, and this background solution can be extracted from the general solution as an asymptotic case with $Q,\ P\ll\mu$.

Acknowledgements.

The work was supported by Russian Science Foundation (Project No. 16-12-10401), and, partially, by the Program of Competitive Growth of Kazan Federal University.