Next Article in Journal
A Prabhakar Fractional Approach for the Convection Flow of Casson Fluid across an Oscillating Surface Based on the Generalized Fourier Law
Next Article in Special Issue
Topological BF Description of 2D Accelerated Chiral Edge Modes
Previous Article in Journal
Investigation and Control of Damping in VSC–MTDC System with DC Circuit Breakers
Previous Article in Special Issue
Symmetries, Partners and Thresholds: The Case of the Xb
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Axion Electrodynamics: Axionically Induced Electric Flares in the Early Magnetized Universe

by
Alexander B. Balakin
1,*,
Vladimir V. Bochkarev
2 and
Albina F. Nizamieva
1
1
Department of General Relativity and Gravitation, Institute of Physics, Kazan Federal University, Kremlevskaya Str. 16a, 420008 Kazan, Russia
2
Department of Radiophysics, Institute of Physics, Kazan Federal University, Kremlevskaya Str. 16a, 420008 Kazan, Russia
*
Author to whom correspondence should be addressed.
Symmetry 2021, 13(11), 2038; https://doi.org/10.3390/sym13112038
Submission received: 6 October 2021 / Revised: 23 October 2021 / Accepted: 26 October 2021 / Published: 29 October 2021
(This article belongs to the Special Issue Physics and Symmetry Section: Feature Papers 2021)

Abstract

:
We consider the nonlinearly extended Einstein–Maxwell-axion theory, which is based on the account for two symmetries: first, the discrete symmetry associated with the properties of the axion field, and second, the Jackson’s symmetry, prescribing to the electrodynamics to be invariant with respect to the rotation in the plane coordinated by the electric and magnetic fields. We derive the master equations of the nonlinearly extended theory and apply them to the Bianchi-I model with magnetic field. The main result, describing the behavior of the nonlinearly coupled axion, electromagnetic, and gravitational fields is the anomalous growth of the axionically induced electric field in the early magnetized Universe. The character of behavior of this anomalous electric field can be indicated by the term flare. We expect, that these electric flares can produce the electron–positron pair creation, significant acceleration of the born charged particles, and the emission of the electromagnetic waves by these accelerated particles.

1. Introduction

The concept of symmetry has played and continues to play an outstanding role in the physics of fundamental interactions. Extending the theoretical models of the nonlinear cosmic axion electrodynamics, we happen to be inspired by beauty of the symmetries of two types. The first type of symmetry is connected with the idea that the Lagrangian of the axion electrodynamics is invariant with respect to the discrete symmetry, associated with the shift ϕ ˜ = ϕ + 2 π k (k is an integer) of the pseudoscalar field ϕ introduced by Peccei and Quinn [1] and indicated later as an axion field (see, e.g., [2,3,4,5,6] for details of this story). According to this, we need to have a Lagrangian, which contains the invariants based on the Maxwell tensor F i k , its dual F i k * , and the axion field ϕ introduced nonlinearly. Clearly, the well-known linear term ϕ F m n * F m n supports the discrete symmetry, since the rest term 2 π k g F m n * F m n in the action functional is the perfect divergence and thus can be eliminated. When one deals with the nonlinear modification of the theory, the nonlinear function of the true invariant f ϕ F m n * F m n does not satisfy the mentioned discrete symmetry, and in this sense, we have to search for new nonlinear terms. For this purpose, we attract for consideration the idea of Jackson [7] that the vacuum Faraday–Maxwell electrodynamics is symmetric with respect to linear transformation of the electric and magnetic fields E * E cos α + B sin α , B * B cos α E sin α , with the constant angle α . Regarding the nonlinear modification of the axion electrodynamics, we assume to use the Lagrangian, which satisfies both mentioned symmetries. The example of such Lagrangian has been constructed in the work [8]. The proposed nonlinear Lagrangian has the form L ( I ) , where L ( I ) is some nonlinear function of its argument, and the invariant I is
I = 1 4 cos ϕ F m n F m n + sin ϕ F m n * F m n .
We assume that for small I the Lagrangian can be transformed into L ( I ) I . In the linear theory, when ϕ 0 , the invariant I can be reduced to the standard invariant
I 1 4 F m n F m n + ϕ F m n * F m n
linear in ϕ and covering the basic term of the linear axion electrodynamics [4,5]. Clearly, the discrete symmetry ϕ ˜ = ϕ + 2 π k is supported identically by the choice of sin/cosine multipliers in front of the standard terms F m n F m n and F m n * F m n . As for the Jackson’s symmetry, if we use the transformation
F i k = F i k cos ϕ 2 + F i k sin ϕ 2 ,
we obtain the term 1 4 F i k F i k , which coincides with (1). This fact confirms the validity of the symmetry transformations
E i = E i cos ϕ 2 + B i sin ϕ 2 , B i = E i sin ϕ 2 + B i cos ϕ 2 ,
which include the electric field four-vector E i = F i k U k and the magnetic induction four-vector B i = F i k U k . The quantity U k is the velocity four-vector of the medium or of an observer.
The nonlinear models with such symmetries attract special attention, since they contain a latent instability which can be displayed by the anomalous growth of the axionically induced electric field. When the model is quasilinear, i.e., when the axion field enters the Lagrangian with sin/cosine nonlinearities but the electromagnetic field is described in the linear context, the problem of electric flares was discussed in the paper [8]. It was shown in [8] that the interaction with the dynamic aether stabilizes the evolution, and the axion field does not reach the catastrophic value ϕ = π 2 , which corresponds to the infinite value of the electric field proportional to the tan ϕ . Now, we exclude the dynamic aether from the model but consider the electromagnetic field in the nonlinear context. Our goal is to show that in the manifold of the guiding parameters of the model there exists an area, in which the presented model is stable. In particular, this means that the axionically induced electric field is finite, but its amplitude can reach very large values, thus producing the pair creation, particle acceleration, and electromagnetic waves emission in the early magnetized Universe.
The paper is organized as follows: In Section 2, we describe the mathematical formalism of the proposed nonlinear extension of the Einstein–Maxwell-axion theory. In Section 3, we reduce the basic equations taking into account the symmetry of the cosmological anisotropic Bianchi-I model with magnetic field and derive the pair of key equations for the axion field and for the derivative of the scale factor. In Section 4, we analyze the solutions to the key equations using qualitative and numerical methods and discuss the problem of the axionically induced electric field generation. Section 5 contains the discussion and conclusions.

2. Nonlinear Einstein–Maxwell-Axion Model

2.1. The Action Functional and Master Equations of the Model

The model is based on the action functional presented in the form
S = d 4 x g R + 2 Λ 2 κ + 1 2 Ψ 0 2 V ( ϕ ) k ϕ k ϕ + L ( I ) .
According to the standards, R is the Ricci scalar, Λ is the cosmological constant, and κ = 8 π G is the Einstein constant (we consider c = 1 ). The dimensionless quantity ϕ describes the pseudoscalar (axion) field, V ( ϕ ) is the potential of the axion field, and k is the covariant derivative. The constant Ψ 0 is reciprocal to the coupling constant of the axion–photon interaction g A γ γ , i.e., 1 Ψ 0 = g A γ γ < 1.47 × 10 10 GeV 1 . Many authors, keeping in mind physical motives, prefer to work with the axion field indicated as a; this quantity can be expressed in terms of ϕ , as follows: a = Ψ 0 ϕ , or equivalently ϕ = g A γ γ · a . In this representation, the term ϕ F m n * F m n in (2) converts into g A γ γ · a F m n * F m n , and the term Ψ 0 2 k ϕ k ϕ in (5) takes the form of k a k a . The new element of the theory is encoded in the last term L ( I ) , the Lagrangian of the electromagnetic field interacting with the axion field; this Lagrangian is the nonlinear function of the argument I , which has at least one continuous derivative and tends to I at I 0 , i.e., L ( I 0 ) I . The argument I is presented by (1).

2.1.1. Master Equations for the Electromagnetic Field

The variation of the action functional (5) with respect to the potential of the electromagnetic field A i yields
k L ( I ) cos ϕ F i k + sin ϕ F i k = 0 ,
where the prime denotes the derivative with respect to the argument I . In addition to the nonlinear electrodynamic Equation (6), we use the standard equations
k F i k = 0 ,
which are the consequence of the definition of the Maxwell tensor F m n = m A n n A m .
As examples, below we consider two model variants of the nonlinear function L ( I ) . The first one is of the power law type
L ( I ) = I + γ I ν , ν > 1 .
The second variant contains the Kohlrausch type function (stretched exponential [9])
L ( I ) = I 1 γ exp I γ 1 , γ > 0 .
When I is small and the parameter γ is positive, the Kohlrausch type function (9) can be linearized as follows: L ( I ) I 1 γ I γ = I , thus recovering the model of linear electrodynamics.

2.1.2. Master Equation for the Axion Field

Variation of the action functional (5) with respect to the axion field ϕ gives the equation
k k ϕ + 1 2 d V ( ϕ ) d ϕ = 1 4 Ψ 0 2 L ( I ) sin ϕ F i k F i k cos ϕ F i k F i k .
We consider the potential V ( ϕ ) to be of the periodic form
V ( ϕ ) = 2 m A 2 1 cos ϕ .
For small values ϕ , the potential (11) converts into the standard one V ( ϕ ) m A 2 ϕ 2 , where m A is the axion mass. Taking into account (11), one can rewrite Equation (10) as follows:
k k ϕ = m A 2 sin ϕ + 1 4 Ψ 0 2 L ( I ) sin ϕ F i k F i k cos ϕ F i k F i k .
This equation explicitly satisfies the requirement of the discrete symmetry ϕ ˜ = ϕ + 2 π k .
Remark 1.
Remark about the effective electromagnetic mass of the axion
For the linear version of electrodynamics and for small values of ϕ, the axion field equation takes the form
k k ϕ + m A 2 ϕ = 1 4 Ψ 0 2 ϕ F i k F i k F i k F i k .
Clearly, in the standard axion electrodynamics, the term ϕ 4 Ψ 0 2 F i k F i k , proportional to the axion field, did not appear. It is the residual term, which appears due to the sin/cosine extension of the Lagrangian and describes the idea that the axion field can mimic the effects attributed to the dilaton field (see [8] for details). Now, we can rewrite (13) in the following form:
k k ϕ + ϕ m A 2 1 4 Ψ 0 2 F i k F i k = 1 4 Ψ 0 2 F i k F i k .
The standard term in the right-hand side of this equation describes the electromagnetic source, which produces the axions. The new term, which contains the quantity
M ( eff ) = m A 2 1 4 Ψ 0 2 F i k F i k ,
can be indicated as the effective axion mass M ( eff ) , since the square of this term appears in front of the pseudoscalar field in the left-hand side of Equation (14). This effective mass depends on the value of the first invariant of the electromagnetic field F i k F i k . This finding is similar to the one known in the dilaton electrodynamics with the Lagrange term 1 8 Φ 2 F m n F m n , where the scalar field Φ describes the dilaton. However, the fact that the coupling constant Ψ 0 2 enters the effective mass M ( eff ) allows us to speak about the axionically induced dilaton-like contribution of the electromagnetic field into the effective axion mass.

2.1.3. Master Equations for the Gravitational Field

The variation of the action functional (5) with respect to the metric gives the gravity field equations
R i k 1 2 R g i k = Λ g i k + κ T i k ( EMA ) + κ T i k ( A ) .
Here, the term
T i k ( A ) = Ψ 0 2 i ϕ k ϕ + 1 2 g i k V ( ϕ ) p ϕ p ϕ
is the stress–energy tensor of the pseudoscalar (axion) field, and the tensor
T i k ( EMA ) = L ( I ) cos ϕ 1 4 g i k F m n F m n F i m F k m + g i k L ( I ) I · L ( I )
contains all the contributions of the electromagnetic field. It coincides with the standard stress–energy tensor of the electromagnetic field, when ϕ = 0 and L ( I ) = I . The trace of the tensor (18)
T i k ( EMA ) g i k = 4 L ( I ) I · L ( I )
is generally nonvanishing; it is equal to zero only in the linear version of the theory. The total stress–energy tensor satisfies the condition
k T i k ( A ) + T i k ( EMA ) = 0
as the differential consequence of (6), (7), and (10).

3. Magnetized Universe: Bianchi-I Model

We apply the established theory to the anisotropic homogeneous cosmological model with magnetic field; the appropriate spacetime platform for this model is the Bianchi-I model (see, e.g., [10,11,12,13,14,15,16] and the references therein). The metric
d s 2 = d t 2 a 2 ( t ) d x 2 b 2 ( t ) d y 2 c 2 ( t ) d z 2
depends on the cosmological time only. We assume that all physical quantities inherit this symmetry and also are a function of t only.

3.1. Exact Solution to the Electromagnetic Field Equations

With the fixed symmetry ansatz, Equation (7) standardly gives the solution for the magnetic field
F 12 = c o n s t B 0 ,
and hint that the electric field parallel to the magnetic field E ( t ) c ( t ) F 30 appears as the result of the axion–photon coupling. Indeed, four equations (6) reduce to the following one:
d d t L ( I ) a b cos ϕ E B 0 sin ϕ = 0 ,
and its solution is
E ( t ) = B 0 a b tan ϕ + const a b cos ϕ L ( I ) .
For the sake of simplicity we consider the model, in which the electric field is equal to zero at the moment t 0 ; moreover, we assume that the initial value of the axion field is equal to π m ( m = 0 , 1 , ), thus providing that the integration constant is equal to zero
E ( t 0 ) = 0 , ϕ ( t 0 ) = π m const = 0 .
With these initial conditions, we obtain the elegant formulas, first, for the electric field
E ( t ) = B 0 a b tan ϕ ,
second, for the invariants of the electromagnetic field
F p q * F p q = 4 B 0 2 a 2 b 2 tan ϕ , F p q F p q = 2 B 0 2 a 2 b 2 1 tan 2 ϕ ,
third, for the invariant I
I = B 0 2 2 a 2 b 2 cos ϕ .

3.2. Reduced Equation for the Axion Field

Equations (10) with (11) can be now transformed into
ϕ ¨ + ( a b c ) · a b c ϕ ˙ + m A 2 sin ϕ = B 0 2 L ( I ) sin ϕ 2 Ψ 0 2 a 2 b 2 cos 2 ϕ ,
where the dot denotes the derivative with respect to the cosmological time. There is also the following convenient form of this equation
ϕ ¨ + ( a b c ) · a b c ϕ ˙ = sin ϕ m A 2 + I L ( I ) Ψ 0 2 cos ϕ ,
which shows explicitly that the formula ϕ = π n with integer n gives the series of exact solutions to the axion field equation in the case of the arbitrary nonlinear term L ( I ) with a finite derivative.

3.3. Reduced Equations for the Gravity Field

Equations (16) with (17) and (18) convert into four nonlinear differential equations of the second order
a ˙ a b ˙ b + a ˙ a c ˙ c + b ˙ b c ˙ c = Λ + 1 2 κ Ψ 0 2 V ( ϕ ) + ϕ ˙ 2 + κ L ( I ) ,
b ¨ b + c ¨ c + b ˙ b c ˙ c = Λ + 1 2 κ Ψ 0 2 V ( ϕ ) ϕ ˙ 2 + κ L ( I ) 2 I L ( I ) ,
a ¨ a + c ¨ c + a ˙ a c ˙ c = Λ + 1 2 κ Ψ 0 2 V ( ϕ ) ϕ ˙ 2 + κ L ( I ) 2 I L ( I ) ,
b ¨ b + a ¨ a + b ˙ b a ˙ a = Λ + 1 2 κ Ψ 0 2 V ( ϕ ) ϕ ˙ 2 + κ L ( I ) .
Since only the magnetic and electric fields directed along the axis 0 z are the sources of the spacetime anisotropy, we can consider the truncated model with the so-called local isotropy, when a ( t ) = b ( t ) . In addition, we introduce the function H = a ˙ a and can rewrite the independent gravity field equations as the following set:
H 2 + 2 H c ˙ c = Λ + 1 2 κ Ψ 0 2 V ( ϕ ) + ϕ ˙ 2 + κ L ( I ) ,
a ¨ a + c ¨ c + a ˙ a c ˙ c = Λ + 1 2 κ Ψ 0 2 V ( ϕ ) ϕ ˙ 2 + κ L ( I ) 2 I L ( I ) ,
2 H ˙ + 3 H 2 = Λ + 1 2 κ Ψ 0 2 V ( ϕ ) ϕ ˙ 2 + κ L ( I ) .
As usual, we see that Equation (36) is the differential consequence of (35) and (37).

3.4. Key Equations of the Model

The difference of Equations (35) and (37) yields
H c ˙ c H ˙ H 2 = 1 2 κ Ψ 0 2 ϕ ˙ 2 ,
and we immediately obtain the scale factor c ( t ) in quadratures
c ( t ) = c ( t 0 ) a ( t ) H ( t ) a ( t 0 ) H ( t 0 ) exp 1 2 κ Ψ 0 2 t 0 t d t H ( t ) ϕ ˙ 2 .
For the model with the local isotropy, it is convenient to use the new variable x introduced as follows:
x = a ( t ) a ( t 0 ) , d d t = x H d d x .
In these terms, the scale factor (39) takes the form
c ( x ) = c ( 1 ) x H ( x ) H ( 1 ) exp 1 2 κ Ψ 0 2 1 x z d z ϕ 2 ( z ) ,
and we obtain the pair of the key equations for two unknown functions H ( x ) and ϕ ( x ) :
1 x 2 d d x ( x 3 H 2 ) = Λ + 1 2 κ Ψ 0 2 2 m A 2 ( 1 cos ϕ ) x 2 H 2 ϕ 2 + κ L ( I ) .
1 x 2 d d x ( x 4 H 2 ϕ ) + 1 2 κ Ψ 0 2 x 3 H 2 ϕ 3 = sin ϕ m A 2 + I L ( I ) Ψ 0 2 cos ϕ .
Here and below depending on the context, we use two equivalent forms for the metric coefficient c ( 1 ) c ( x = 1 ) = c ( t 0 ) and for the Hubble function H ( 1 ) H ( x = 1 ) = H ( t 0 ) . We have to keep in mind now that the argument I of the nonlinear Lagrangian can be written as
I = ω 2 x 4 cos ϕ , ω 2 B 0 2 2 a 4 ( t 0 ) ,
and that the scale factor a ( t ) can be reconstructed using the formula
t t 0 = 1 a ( t ) a ( t 0 ) d x x H ( x ) .

4. Analysis of the Model: Exact Solutions and Numerical Study

The plan for the investigations of the presented model is the following: First of all, we consider the model with L ( I ) = I , as an intermediate case, i.e., the axion field is presented in the nonlinear context, but the electrodynamics is linear. The main question which we have to discuss here can be formulated as follows: are the solutions for the electric field finite or infinite, or equivalently, is the model stable or unstable? The second step is to consider three particular nonlinear models, for which the Hubble functions are found analytically, but the axion field and electric field are obtained using numerical calculations.

4.1. Quasilinear Model

4.1.1. The Test of Instability

For the model with L ( I ) = I , two key Equations (42) and (43) can be rewritten in the form:
1 x 2 d d x ( x 3 Y 2 ) = Λ * + κ Ψ 0 2 μ A 2 ( 1 cos ϕ ) 1 2 x 2 Y 2 ϕ 2 + 1 x 4 cos ϕ .
1 x 2 d d x ( x 4 Y 2 ϕ ) + 1 2 κ Ψ 0 2 x 3 Y 2 ϕ 3 = sin ϕ μ A 2 + 1 x 4 cos 2 ϕ .
In order to simplify the numerical simulation, we introduced here the auxiliary quantities Y ( x ) , μ A , Λ * given by
H ( x ) = Y ( x ) H * , m A = μ A H * , Λ = Λ * H * 2 , H * | B 0 | 2 a 2 ( t 0 ) Ψ 0 .
Similarly, we can introduce the dimensionless electric field e ( x ) as follows:
E ( x ) = e ( x ) E * , E * | B 0 | a 2 ( t 0 ) e ( x ) = tan ϕ ( x ) x 2 .
When ϕ π 2 , the electric field (26), the invariants of the electromagnetic fields (27), and the basic invariant (28) tend toward infinity. Respectively, the right-hand sides of the key Equations (46) and (47) also tend toward infinity in this limit. We have to test, whether the axion field ϕ can reach the catastrophic value ϕ = π 2 . For this purpose, we consider two submodels. In the first one, the initial value of the axion field is chosen as ϕ ( t 0 ) = 0 ; this value corresponds to the first minimum of the axion potential (11). In the second submodel, we assume that the axion field evolution starts with the value ϕ = π , which relates to the first maximum of the axion potential.

4.1.2. The First Submodel Is Stable

When ϕ ˙ ( t 0 ) = 0 , we see that ϕ ( x ) 0 is the exact solution to Equation (47). As for Equation (46), it can be directly integrated yielding
Y 2 ( x ) = 1 x 3 Y 2 ( 1 ) + 1 3 Λ * 1 1 x 3 + κ Ψ 0 2 1 x 3 1 x 4 .
Since E ( x ) = 0 for this solution, we recover the well-known solution attributed to the Bianchi-I model with the pure magnetic field. Asymptotically, this solution describes the de Sitter model with H Λ 3 and is characterized by isotropization c a c o n s t .
When ϕ ˙ ( t 0 ) 0 , the numerical analysis of the key Equations (46) and (47) shows that all the curves ϕ ( x ) are located below the horizontal line ϕ = π 2 . This means that the electric field E, Hubble function H and scale factors a, c take finite values, i.e., the catastrophic regime cannot be realized. We illustrate this fact using Figure 1, Figure 2 and Figure 3.
For the model under discussion, the oscillatory regime of the axion field evolution is predetermined by the structure of Equation (47). Indeed, in the vicinities of the extrema points, when ϕ 0 , the mentioned equation can be rewritten as
ϕ = 1 x 2 Y 2 μ A 2 + 1 x 4 cos 2 ϕ sin ϕ .
Clearly, the sign of the second derivative depends on the sign of the function sin ϕ : when ϕ > 0 , we see that ϕ < 0 , and we observe the maxima; when ϕ < 0 , we deal with the minima. If ϕ π 2 , the corresponding second derivative becomes large due to the term proportional to 1 cos 2 ϕ , this means that the degree of convexity/concavity of the graph increases, but the typical structure of the extrema is preserved. The damping of these oscillations is provided by the multipliers of the type 1 x 4 and 1 x 2 , which decrease when the reduced scale factor x = a ( t ) a ( t 0 ) grows.

4.1.3. The Second Submodel Is Unstable

When we put ϕ ( t 0 ) = π and ϕ ˙ ( t 0 ) = 0 , we again obtain the constant exact solution ϕ π to Equation (47). The Hubble function is now described by the formula
Y 2 ( x ) = 1 x 3 Y 2 ( 1 ) + 1 3 Λ ˜ 1 1 x 3 κ Ψ 0 2 1 x 3 1 x 4 ,
where the redefined cosmological constant is Λ ˜ = Λ * + 2 κ Ψ 0 2 μ A 2 . Again, at x , we obtain the de Sitter type behavior of the model.
When ϕ ˙ ( t 0 ) 0 , the numerical analysis of the key Equations (46) and (47) gives principally different results: the horizontal line ϕ = π 2 happens to be crossed for the arbitrary nonvanishing initial value of the derivative ϕ ˙ ( t 0 ) . This submodel is unstable. We illustrate this fact by Figure 4.

4.2. Analysis of the Behavior of Nonlinear Systems in the Vicinity of the First Minimum of the Axion Potential: Three Explicit Examples

4.2.1. Basic Equilibrium Solutions and Equation of Perturbation Dynamics

The first minimum of the axion potential corresponds to the value ϕ = 0 ; we add to the axion field a small spatially homogeneous perturbation ψ ( t ) . In the leading order approximation, the electric field (26) vanishes, and the key Equation (42) takes the form
1 x 2 d d x ( x 3 H 2 ) = Λ + κ L ( I ) .
Since now the invariant I is of the form I ω 2 x 4 z , we can rewrite the key equation as follows:
4 z d H 2 d z + 3 H 2 = Λ + κ L ( z ) .
The solution for H ( x ) can be presented in quadratures
H 2 ( x ) = H 2 ( 1 ) x 3 + Λ 3 ( 1 x 3 ) 1 4 κ ω 3 2 x 3 ω 2 ω 2 x 4 d ξ ξ 7 4 L ( ξ ) .
The scale factor c ( x ) also can be simplified:
c ( x ) = c ( 1 ) x H ( x ) H ( 1 ) .
In the first-order approximation with respect to ψ ( t ) , we obtain the equation
1 x 2 d d x x 4 H 2 ( x ) d d x ψ + ψ m A 2 + ω 2 L ( z ) x 4 Ψ 0 2 = 0 .
In the linear approximation, the electric field is proportional to ψ , i.e., E ( x ) = ψ x 2 E * .

4.2.2. The First Example: Power-Law Lagrange Function

Integration in (55) with the function (8) gives
H 2 ( x ) = H 2 ( 1 ) x 3 + Λ 3 ( 1 x 3 ) +
κ ω 2 x 4 ( x 1 ) + γ ( 4 ν 3 ) ω 2 x 4 ν x 4 ν 3 1 .
Asymptotically, we deal with the de Sitter law H ( x ) Λ 3 .
Figure 5 illustrates the behavior of H ( x ) , of the scale factor a ( t ) and of the ratio c ( t ) a ( t ) for this model. Figure 6 presents the corresponding graphs of the function ψ ( t ) .

4.2.3. The Second Example: Kohlrausch Type Function with γ = 1 8

Now, we consider the Lagrange function (9) with γ = 1 8 . The function H ( x ) can be presented in the explicit analytic form
H 2 ( x ) = H 2 ( 1 ) x 3 + Λ 3 ( 1 x 3 ) +
2 κ ω 3 x 6 ω 1 4 1 x 1 + exp ω 1 4 exp ω x 2 1 4 .
Again, we deal with the asymptotically de Sitter type behavior; however, in contrast to the results for the quasilinear model, one can see that H ( x ) and c ( t ) a ( t ) are nonmonotonic.
Figure 7 illustrates the behavior of H ( x ) , of the scale factor a ( t ) and of the ratio c ( t ) a ( t ) for this model. Figure 8 presents the corresponding graphs of the function ψ ( t ) .

4.2.4. The Third Example: Kohlrausch Type Function with γ = 2 (Anti-Gaussian Function)

When we consider the model Lagrangian of the form
L ( I ) = I + β I 11 4 e I 2 1 ,
we obtain explicit analytic formula for the Hubble function; its square is
H 2 ( x ) = H 2 ( 1 ) x 3 + Λ 3 ( 1 x 3 ) +
κ ω 2 x 4 ( x 1 ) + β 8 ω 2 x 4 11 4 ( 1 x 8 ) + ω 2 x 4 3 4 e ω 4 e ω 4 x 8 .
Figure 9 illustrates the behavior of H ( x ) , of the scale factor a ( t ) and of the ratio c ( t ) a ( t ) for this model. Figure 10 presents the corresponding graphs of the function ψ ( t ) .

5. Discussion and Conclusions

5.1. Axionically Induced Electric Field: Anomalies and Electric Flares

The mechanism of generation of the axionically induced electric field in the presence of the magnetic field is well known and well documented. This mechanism is based on the phenomenon of axion-photon coupling, and it has a lot of interesting applications (see, e.g., [5,6,17]). The new trend in the axion electrodynamics is connected with the elaboration of the nonlinear axion electrodynamics, which could have astrophysical and cosmological applications. In this sense, the early Universe with strong magnetic field [18,19,20], filled with the axionic dark matter [21,22,23,24], seems to be one of the most interesting objects of research. Indeed, it is well known that the dark matter axions interacting with the cosmic magnetic field produce the electric field; however, as the linear axion electrodynamics predicts, the contribution of such an electric field into the energetic balance of the Universe seems to be negligible, since the coupling constant of the axion–photon interactions is rather small, 1 Ψ 0 = g A γ γ < 1.47 × 10 10 GeV 1 . The situation changes principally when we work with nonlinear axion electrodynamics. As we have shown above, the electric field with zero initial value can grow proportionally to the function tan ϕ . This means that if the dimensionless function ϕ , describing the evolution of the axion field, tends to the value ϕ = π 2 , the axionically induced electric field grows abnormally. We have studied the special question: is the catastrophic value ϕ = π 2 reachable? The answer is negative; however, the maximum value of the electric field can be very large and even huge, depending on the initial value of the derivative ϕ ˙ ( t 0 ) . As was shown above, the anomalous behavior of the electric field can be indicated as flare, since this electric field rapidly decreases in amplitude; this phenomenon can occur only once in the early Universe. If the initial value of the axion field is equal to ϕ ( t 0 ) = π , i.e., if it corresponds to the first maximum of the axion field potential, the catastrophic value ϕ = π 2 is achieved inevitably; this model happens to be singular.
What are the physical consequences of such electric flares? From our point of view, one can expect, first, the creation of the electron–positron pairs; second, the acceleration of charged particles; and third, the emission of the electromagnetic waves by these accelerated particles.

5.2. On the Electron–Positron Pair Creation in the Axionically Induced Electric Field

According to Schwinger’s formula, the probability of the electron–positron pair creation in vacuum under the influence of a strong electric field E is estimated to be the following:
W = e 2 E 2 4 π 3 2 exp E ( c r i t ) E , E ( c r i t ) = π m e 2 e , ( c = 1 ) .
Clearly, the nonlinear axion electrodynamics of the early magnetized Universe does not prohibit the anomalous electric field energy to overcome the threshold 2 m e c 2 , thus providing the creation of the electron–positron pairs in vacuum.

5.3. Particle Acceleration

If the initial three-momenta of the electron and positron, born in the strong electric field, are equal to zero, these particles begin to move along the direction 0 z indicated by the parallel electric and magnetic fields. For this field configuration, the particles have only two components of the velocity four-vector: U z and U 0 = 1 + c 2 ( t ) U z 2 , and the longitudinal component of the acceleration four-vector. From the equation of the particle dynamics
d P j d s = Γ j l k P k P l m e + e m e F j l P l ,
where P i = m e U i is the particle four-momentum, we obtain for the longitudinal acceleration and velocity, respectively, the following relationships:
d U z d t = e c ( t ) B 0 m e a 2 ( t ) tan ϕ ( t ) , U z ( t ) = e B 0 m e t 0 t d t c ( t ) a 2 ( t ) tan ϕ ( t ) .
These kinematic characteristics depend on time, and thus, we deal with the nonuniform accelerated motion. Clearly, electrons and positrons move in opposite directions and can take part in the process of anomalous acceleration, if ϕ ( t ) π 2 .

5.4. On the Radiation of Electromagnetic Waves by the Charged Particles Accelerated by the Axionically Induced Electric Field

The nonuniformly moving particles emit the electromagnetic wave. We consider the radiation of particles, which have the parallel three-velocity v and three-acceleration w . In this case, the standard formulas for the radiation intensity per single solid angle has the form
d I d Ω = e 2 4 π w 2 sin 2 θ 1 v cos θ 6 , v = | v | , w = | w | , ( c = 1 ) ,
where θ is the angle between the axis 0 z and the direction to the observation point (see, e.g., [25]). In our case, for the very sharp peak of the electric flare, we can use the formulas v U z U 0 and w d U z d t (see (64)). In the case of the anomalous axionically induced electric field, the charged particle born in this field can reach the ultrarelativistic state of motion, and the radiation is focused in the small angle θ . The intensity of radiation is proportional to B 0 2 tan 2 ϕ . In other words, one can expect that the directions θ = 0 and θ = π are marked by the extra portions of light.

5.5. Conclusions

We presented the new version of the nonlinear extension of the Einstein–Maxwell-axion theory, and based on this theory, we can propose the following scenario of the evolution of the spatially homogeneous early magnetized Universe.
  • Let the axion field evolution start with zero value ϕ ( t 0 ) = 0 and with the nonvanishing derivative ϕ ˙ ( t 0 ) 0 ; in the presence of the cosmological magnetic field, the axion–photon coupling provides the creation of the axionically induced electric field parallel to the magnetic field; in its turn, such electromagnetic field configuration produces a new portion of axions, thus unwinding the spiral of the nonlinear axion–photon interactions.
  • In the process of evolution, the axion field and its derivative remain finite; the axion field itself happens to be trapped in the zone π 2 < ϕ < π 2 and thus does not play an essential role in the Universe expansion.
  • The axionically induced electric field, in contrast to the axion field, can reach anomalously large values, when ϕ π 2 ; the so-called electric flares can appear, and the electromagnetic contribution to the sources of the gravity field can become essential.
  • The geometric characteristics of the Universe, such as scale factors and their derivatives, remain nonsingular at t > t 0 and inherit the oscillatory behavior of the electromagnetic sources; when the cosmological constant is nonvanishing, the final stage of the Universe evolution is of the de Sitter type and thus corresponds to the instruction about the late-time accelerated expansion.
  • The axionically induced anomalous electric field produces the electron–positron pairs and accelerates the born charged particles; the accelerated charged particles emit the electromagnetic waves, which are focused in the direction pointed by the parallel magnetic and electric fields; these flares can take part in the formation of fluctuations of the cosmic microwave background.
In the next work, we hope to revive this scenario by adding estimations and other necessary details.

Author Contributions

Conceptualization, A.B.B.; methodology, A.B.B. and V.V.B.; formal analysis, A.B.B. and V.V.B.; investigation, A.B.B., V.V.B., and A.F.N.; writing—original draft preparation, A.B.B.; writing—review and editing, A.B.B. and A.F.N.; visualization, V.V.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Russian Foundation for Basic Research (Grant N 20-52-05009).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The study did not report any data.

Acknowledgments

A.B. is grateful for the Russian Foundation for Basic Research for financial support (Grant N 20-52-05009). This paper was supported by the Kazan Federal University Strategic Academic Leadership Program.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Peccei, R.D.; Quinn, H.R. CP conservation in the presence of instantons. Phys. Rev. Lett. 1977, 38, 1440–1443. [Google Scholar] [CrossRef] [Green Version]
  2. Weinberg, S. A new light boson? Phys. Rev. Lett. 1978, 40, 223–226. [Google Scholar] [CrossRef]
  3. Wilczek, F. Problem of strong P and T invariance in the presence of instantons. Phys. Rev. Lett. 1978, 40, 279–282. [Google Scholar] [CrossRef]
  4. Ni, W.-T. Equivalence principles and electromagnetism. Phys. Rev. Lett. 1977, 38, 301–304. [Google Scholar] [CrossRef]
  5. Sikivie, P. Experimental tests of the “invisible” axion. Phys. Rev. Lett. 1983, 51, 1415–1417. [Google Scholar] [CrossRef]
  6. Wilczek, F. Two applications of axion electrodynamics. Phys. Rev. Lett. 1987, 58, 1799–1802. [Google Scholar] [CrossRef] [PubMed]
  7. Jackson, J.D. Classical Electrodynamics; John Wiley and Sons: New York, NY, USA, 1999. [Google Scholar]
  8. Balakin, A.B.; Galimova, A.A. Towards nonlinear axion-dilaton electrodynamics: How can axionic dark matter mimic dilaton-photon interactions? Phys. Rev. D 2021, 104, 044059. [Google Scholar] [CrossRef]
  9. Kohlrausch, R. Theorie des elektrischen Ruckstandes in der Leidner Flasche. Ann. Phys. Chem. 1854, 91, 56–82. [Google Scholar] [CrossRef] [Green Version]
  10. Turner, M.S.; Widrow, L.M. Inflation-produced, large-scale magnetic fields. Phys. Rev. D 1988, 37, 2743–2769. [Google Scholar] [CrossRef]
  11. Grasso, D.; Rubinstein, H.R. Magnetic fields in the early Universe. Phys. Rept. 2001, 348, 163–266. [Google Scholar] [CrossRef] [Green Version]
  12. Lee, D.-S.; Lee, W.; Ng, K.-W. Primordial magnetic fields from dark energy. Phys. Lett. B 2002, 542, 1–7. [Google Scholar] [CrossRef] [Green Version]
  13. Giovannini, M. The magnetized Universe. Int. J. Mod. Phys. D 2004, 13, 391–502. [Google Scholar] [CrossRef]
  14. Bronnikov, K.A.; Chudaeva, E.N.; Shikin, G.N. Magneto-dilatonic Bianchi-I cosmology: Isotropization and singularity problems. Class. Quantum Gravity 2004, 21, 3389–3403. [Google Scholar] [CrossRef] [Green Version]
  15. Bamba, K.; Odintsov, S.D. Inflation and late-time cosmic acceleration in non-minimal Maxwell-F(R) gravity and the generation of large-scale magnetic fields. J. Cosmol. Astropart. Phys. 2008, 4, 24. [Google Scholar] [CrossRef] [Green Version]
  16. Balakin, A.B. Magnetic relaxation in the Bianchi-I universe. Class. Quantum Gravity 2007, 24, 5221–5245. [Google Scholar] [CrossRef]
  17. Balakin, A.B.; Muharlyamov, R.K.; Zayats, A.E. Axion-induced oscillations of cooperative electric field in a cosmic magneto-active plasma. Eur. Phys. J. D 2014, 68, 159. [Google Scholar] [CrossRef] [Green Version]
  18. Novello, M.; Goulart, E.; Salim, J.M.; Perez Bergliaffa, S.E. Cosmological effects of nonlinear electrodynamics. Class. Quantum Gravity 2007, 24, 3021–3036. [Google Scholar] [CrossRef] [Green Version]
  19. Kruglov, S.I. Universe acceleration and nonlinear electrodynamics. Phys. Rev. D 2015, 92, 123523. [Google Scholar] [CrossRef] [Green Version]
  20. Övgün, A.; Leon, G.; Magaña, J.; Jusufi, K. Falsifying cosmological models based on a non-linear electrodynamics. Eur. Phys. J. C 2018, 78, 462. [Google Scholar] [CrossRef] [Green Version]
  21. Duffy, L.D.; van Bibber, K. Axions as dark matter particles. New J. Phys. 2009, 11, 1050088. [Google Scholar] [CrossRef]
  22. Steffen, F.D. Dark Matter candidates—Axions, neutralinos, gravitinos, and axinos. Eur. Phys. J. C 2009, 59, 557–588. [Google Scholar] [CrossRef] [Green Version]
  23. Marsh, D.J.E. Axion Cosmology. Phys. Rept. 2016, 643, 1–79. [Google Scholar] [CrossRef] [Green Version]
  24. Sikivie, P. Invisible axion search methods. Rev. Mod. Phys. 2021, 93, 15004. [Google Scholar] [CrossRef]
  25. Landau, L.D.; Lifshitz, E.M. The Classical Theory of Field; Pergamon Press: Oxford, UK, 1971. [Google Scholar]
Figure 1. This figure illustrates the behavior of the axion field with the initial value ϕ ( t 0 ) = 0 as a function of the initial value of the axion field derivative ϕ ˙ ( t 0 ) ; three values of this parameter are fixed in the lower-right corner of the panel. All the values of the axion field ϕ belong to the interval π 2 < ϕ < π 2 . The character of evolution of the axion field can be characterized as damping oscillations. The sharpness of the first line in the zones of the first maximum and first minimum is the result of the choice of a special unified scale, which allowed us to visualize the oscillation pattern simultaneously for both: small and large values of the scale factor a.
Figure 1. This figure illustrates the behavior of the axion field with the initial value ϕ ( t 0 ) = 0 as a function of the initial value of the axion field derivative ϕ ˙ ( t 0 ) ; three values of this parameter are fixed in the lower-right corner of the panel. All the values of the axion field ϕ belong to the interval π 2 < ϕ < π 2 . The character of evolution of the axion field can be characterized as damping oscillations. The sharpness of the first line in the zones of the first maximum and first minimum is the result of the choice of a special unified scale, which allowed us to visualize the oscillation pattern simultaneously for both: small and large values of the scale factor a.
Symmetry 13 02038 g001
Figure 2. This figure presents three graphs of the reduced electric field e ( x ) = E ( x ) E * , which relate to the graphs of ϕ depicted on Figure 1. ( E * = | B 0 | a 2 ( t 0 ) ). The electric field reaches very big values for the time moments, when the values of the axion field are close to ± π 2 . The behavior of the electric field inherits the oscillatory regime of the axion field; however, this property is hidden on the graph because of the incommensurability of the amplitudes of the first and second maxima of the corresponding graphs.
Figure 2. This figure presents three graphs of the reduced electric field e ( x ) = E ( x ) E * , which relate to the graphs of ϕ depicted on Figure 1. ( E * = | B 0 | a 2 ( t 0 ) ). The electric field reaches very big values for the time moments, when the values of the axion field are close to ± π 2 . The behavior of the electric field inherits the oscillatory regime of the axion field; however, this property is hidden on the graph because of the incommensurability of the amplitudes of the first and second maxima of the corresponding graphs.
Symmetry 13 02038 g002
Figure 3. Illustrations of the geometric properties of the model. The behavior of the reduced Hubble function Y ( x ) = H ( x ) H * is presented on the left panel; the central panel contains the graph of the scale factor a ( t ) ; and the right panel presents the evolution of the ratio c ( t ) a ( t ) . Oscillations of the axion field declare themselves in the form of ripples on the graph of the Hubble function. Asymptotically, H ( x ) tends toward the de Sitter constant Λ 3 , and the ratio c ( t ) a ( t ) tends toward one, thus confirming the fact of the Universe isotropization.
Figure 3. Illustrations of the geometric properties of the model. The behavior of the reduced Hubble function Y ( x ) = H ( x ) H * is presented on the left panel; the central panel contains the graph of the scale factor a ( t ) ; and the right panel presents the evolution of the ratio c ( t ) a ( t ) . Oscillations of the axion field declare themselves in the form of ripples on the graph of the Hubble function. Asymptotically, H ( x ) tends toward the de Sitter constant Λ 3 , and the ratio c ( t ) a ( t ) tends toward one, thus confirming the fact of the Universe isotropization.
Symmetry 13 02038 g003
Figure 4. Illustrations of the behavior of the axion field with initial value ϕ ( t 0 ) = π . The graphs are monotonic, and the axion field reaches the value ϕ = π 2 for the finite value of the reduced scale factor x. The corresponding values of the electric field are infinite. The geometric characteristics of the Universe reveal singularity.
Figure 4. Illustrations of the behavior of the axion field with initial value ϕ ( t 0 ) = π . The graphs are monotonic, and the axion field reaches the value ϕ = π 2 for the finite value of the reduced scale factor x. The corresponding values of the electric field are infinite. The geometric characteristics of the Universe reveal singularity.
Symmetry 13 02038 g004
Figure 5. Illustration of the behavior of the Hubble function H ( x ) , of the scale factor a ( t ) and of the ratio c ( t ) a ( t ) for a few values of the parameters ν and ω attributed to the Power-Law Lagrange function.
Figure 5. Illustration of the behavior of the Hubble function H ( x ) , of the scale factor a ( t ) and of the ratio c ( t ) a ( t ) for a few values of the parameters ν and ω attributed to the Power-Law Lagrange function.
Symmetry 13 02038 g005
Figure 6. Illustration of the behavior of the function ψ for the model with Power-Law Lagrange function. Perturbations remain finite, and the model is stable.
Figure 6. Illustration of the behavior of the function ψ for the model with Power-Law Lagrange function. Perturbations remain finite, and the model is stable.
Symmetry 13 02038 g006
Figure 7. Illustration of the behavior of the function H ( x ) , a ( t ) and of the ratio c ( t ) a ( t ) for the Kohlrausch model with γ = 1 8 .
Figure 7. Illustration of the behavior of the function H ( x ) , a ( t ) and of the ratio c ( t ) a ( t ) for the Kohlrausch model with γ = 1 8 .
Symmetry 13 02038 g007
Figure 8. Illustration of the behavior of the function ψ for the Kohlrausch model with γ = 1 8 . Perturbations of the axion field remain finite and their moduli do not exceed the value π 2 .
Figure 8. Illustration of the behavior of the function ψ for the Kohlrausch model with γ = 1 8 . Perturbations of the axion field remain finite and their moduli do not exceed the value π 2 .
Symmetry 13 02038 g008
Figure 9. Illustration of the behavior of the functions H ( x ) , a ( t ) , and c ( t ) a ( t ) for the anti-Gaussian Lagrange function (60).
Figure 9. Illustration of the behavior of the functions H ( x ) , a ( t ) , and c ( t ) a ( t ) for the anti-Gaussian Lagrange function (60).
Symmetry 13 02038 g009
Figure 10. Illustration of the behavior of the function ψ ( t ) for the anti-Gaussian Lagrange function (60).
Figure 10. Illustration of the behavior of the function ψ ( t ) for the anti-Gaussian Lagrange function (60).
Symmetry 13 02038 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Balakin, A.B.; Bochkarev, V.V.; Nizamieva, A.F. Nonlinear Axion Electrodynamics: Axionically Induced Electric Flares in the Early Magnetized Universe. Symmetry 2021, 13, 2038. https://doi.org/10.3390/sym13112038

AMA Style

Balakin AB, Bochkarev VV, Nizamieva AF. Nonlinear Axion Electrodynamics: Axionically Induced Electric Flares in the Early Magnetized Universe. Symmetry. 2021; 13(11):2038. https://doi.org/10.3390/sym13112038

Chicago/Turabian Style

Balakin, Alexander B., Vladimir V. Bochkarev, and Albina F. Nizamieva. 2021. "Nonlinear Axion Electrodynamics: Axionically Induced Electric Flares in the Early Magnetized Universe" Symmetry 13, no. 11: 2038. https://doi.org/10.3390/sym13112038

APA Style

Balakin, A. B., Bochkarev, V. V., & Nizamieva, A. F. (2021). Nonlinear Axion Electrodynamics: Axionically Induced Electric Flares in the Early Magnetized Universe. Symmetry, 13(11), 2038. https://doi.org/10.3390/sym13112038

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop