Stability of Accelerated Expansion in Nonlinear Electrodynamics

This paper is devoted to study the phase space analysis of isotropic and homogenous universe model by taking a noninteracting mixture of electromagnetic and viscous radiating fluids whose viscous pressure satisfies a nonlinear version of the Israel-Stewart transport equation. We establish an autonomous system of equations by introducing normalized dimensionless variables. In order to analyze stability of the system, we find corresponding critical points for different values of the parameters. We also evaluate power-law scale factor whose behavior indicates different phases of the universe model. It is concluded that bulk viscosity as well as electromagnetic field enhances the stability of accelerated expansion of the isotropic and homogeneous universe model.


Introduction
Many astronomical observations (type Ia supernova, large scale structure and cosmic microwave background radiation) predict that our universe is expanding at an accelerating rate in its present stage [1]. These observations suggest two cosmos phases, i.e., the cosmic state before radiation (the primordial inflationary era) and ultimately the present cosmos phase after the matter dominated era. In the last couple of decades, it is speculated that the source for this observed cosmic acceleration with unusual anti-gravitational force may be an anonymous energy component dubbed as dark energy (DE). The existence of this energy with large negative pressure can be recognized by its distinctive nature from ordinary matter which may lead to cosmic expansion. The study of the dominant contents of matter distribution in the universe has remained one of the most challenging issues. Recent observations show that the visible part of our universe is made up of baryonic matter contributing only 5% of the total budget while the remaining ingredients yield the total energy density composed of non-baryonic fluids (68% DE and 27% dark matter) [2].
Several cosmological proposals have been introduced in literature to explore the ambiguous nature of DE. The cosmological constant (Λ) governed by a negative equation of state (EoS) parameter (γ = −1) is taken to be the simplest characterization of DE. However, this identification has two wellknown problems, i.e., fine-tuning and cosmic coincidence. In addition, there are several cosmological models which can be considered as an alternative to Λ like scalar field model [3], phantom model [4], tachyon field [5] and kessence [6] that also suggest expanding behavior of the universe. Another approach involves the generalization of simple barotropic EoS to more exotic forms such as Chaplygin gas [7] and its modification [8]. It has also been demonstrated that a fluid with bulk viscosity may cause accelerated expansion of the universe model without cosmological constant or scalar field [9]. Our main concern is to find another approach which can minimize exotic forms of matter by introducing dissipation through viscous effects of fluids.
During the last few years, cosmological models including nonlinear electromagnetic fields have attained remarkable interest. The application of this electrodynamics to different universe models may lead to many significant results. Nonlinear electrodynamics (NLED) is the generalization of Maxwell theory which is considered as the most viable theory to remove the initial singularities. Vollick [10] considered FRW universe model with NLED and found that the respective model will show a period of late-time acceleration for E 2 < 3B 2 . Kruglov [11] found that the universe tends to accelerate in magnetic background at the early era due to NLED model. Ovgun [12] formulated analytical nonsingular extension of isotropic and homogeneous solutions by presenting a new mathematical model in nonlinear magnetic monopole fields.
The study of possible stable late-time attractors has attained remarkable significance for different universe models. A phase space analysis manifests dynamical behavior of a cosmological model through a global view by reducing complexity of the equations (converting the system of equations to an autonomous system) which may help to understand different stages of evolution. Copeland et al. [13] studied a phase plane analysis of standard inflationary models and found that these models cannot solve density problem. Guo et al. [14] explored phase space analysis of FRW universe model filled with barotropic fluid and phantom scalar field in which phantom dominated solution is found to be a stable late-time attractor.
Garcia-Salcedo [15] examined the dynamics of FRW universe with NLED and found that the critical points have no effects. Yang and Gao [16] discussed phase space analysis of k-essence cosmology in which critical points play an important role for the final state of the universe. Xiao and Zhu [17] analyzed stability of FRW universe model in loop quantum gravity via phase space portraits by taking barotropic fluid as well as positive field potential. Acquaviva and Beesham [18] studied phase space analysis of FRW spacetime filled with noninteracting mixture of fluids (dust and viscous radiation) and found that nonlinear viscous model shows the possibility of current cosmic expansion.
This paper is devoted to study the phase space analysis of FRW universe model with nonlinear viscous fluid. The plan of the paper is as follows. In section 2, we provide basic formalism for NLED and general equations as well as a nonlinear model for bulk viscosity. An autonomous system of equations is established to analyze stability of the system by introducing normalized dimensionless variables in section 3. Section 4 provides the formulation of power-law scale factor. Finally, we conclude our results in the last section.

Nonlinear Electrodynamics and General Equations
The standard cosmological model is successful in resolving many issues but still there are some issues which remain to be solved. One of them the initial singularity which leads to a troubling state of affairs because at this point, all known physical theories break down. If the early universe is governed by Maxwell's equations, then there will be a spacelike initial singularity in the past. However, if Maxwell's equations become modified in the early universe, when the electromagnetic field is large, it might help avoiding the occurrence of cosmic singularities [19]. For the situations where strong electromagnetic field occurs, it makes sense to couple gravitation with NLED. The coupling of Einstein gravity with NLED is defined by the action We consider nonlinear extension of Maxwell Lagrangian density up to second order terms in the field invariants F = F µν F µν and F * = F * µν F µν given by [10] where µ 0 denotes magnetic permeability, α, β > 0 are arbitrary constants which yield linear density for α, β → 0 and F * µν is the dual of electromagnetic field tensor. We do not consider the term F F * involving F * in order to preserve the parity [20,21]. The linear term of this Lagrangian dominates during radiation dominated era while the quadratic terms dominate in the early universe that corresponds to the bouncing behavior of the universe to avoid initial singularity [22]. The mechanisms behind the bounce have been demonstrated in [23,25]. The energy-momentum tensor associated with this Lagrangian has the following form In order to fulfill the requirement of isotropic and homogeneous universe, i.e., the electromagnetic field to act as its source, the energy density and the pressure corresponding to the electromagnetic field can be computed by averaging over volume [21,23]. It is assumed that electric and magnetic fields have coherent lengths that are much shorter than the cosmological horizon scales. After several conditions, the energy momentum tensor of the electromagnetic field associated with L(F, F * ) can be written as that of a perfect fluid such that where ∂ F represents partial derivative with respect to F = F µν F µν = 2(B 2 − E 2 ), E and B denote the averaged electric and magnetic fields, respectively. We consider isotropic and homogeneous universe model given by where a(t) is the scale factor. We assume the universe model to be filled with two cosmic fluids, i.e., a noninteracting electromagnetic fluid with energy density ρ EM as well as pressure p EM and a viscous fluid having energy density ρ v as well as pressure p = p v (ρ v ) + Ψ. Here p v represents the equilibrium part of viscous pressure whereas Ψ is the non-equilibrium part, i.e., bulk viscous pressure satisfying an evolution equation. Bulk viscosity plays an important role to stabilize the density evolution and overcomes the rapid changes in cosmos. It also promotes negative energy field in the fluid and hence can play the role of dark energy to describe the dynamics of cosmos. It has been suggested that a fluid with bulk viscosity may cause an accelerated expansion of the universe model without cosmological constant or scalar field [24]. The main contribution of bulk viscosity to the effective pressure is its dissipative effect. We obtain Raychaudhuri and constraint equations from the field equations given bẏ where dot means derivative with respect to time. The conservation of energymomentum tensor yields the following evolution equations for viscous and electromagnetic field componentṡ We consider a barotropic EoS for viscous fluid defined by where 1 ≤ γ ≤ 2. Using Eqs. (8) and (9), Raychaudhuri and conservation equations for viscous fluid turn out to bė We characterize the viscous pressure variable by the following evolution equation [26] where ζ, T , τ and τ * denote bulk viscosity, local equilibrium temperature, linear relaxation time and characteristic time in nonlinear background, respectively. This equation is derived by using a nonlinear model describing a relationship between thermodynamic flux "Ψ" and thermodynamic force "χ" in the form This is a nonlinear extension of Israel-Stewart equation which reduces to its linear form as τ * → 0. The nonlinear term in Eq.(15) must be positive for thermodynamic consistency and positivity of entropy production rate. The parameters involved in Eq. (15) can be defined by the relations Here k is a constant such that k = 0 gives linear (Israel-Steward) case while T 0 represents constant temperature. Also, ν corresponds to the dissipative effect of the speed of sound V such that V 2 = c 2 s + ν 2 , where c 2 s is its adiabatic contribution. By causality, V ≤ 1 and c 2 s = γ − 1 which yields The explicit form of evolution equation by using the above relations yieldṡ

Phase Space Analysis
In this section, we discuss the phase space analysis of isotropic and homogeneous universe model for radiation case. Due to many arbitrary parameters, it seems difficult to find analytical solution of the evolution equation. In this context, we define normalized dimensionless variables Ω = 3ρv Θ 2 andΨ = 3Ψ Θ 2 such that the corresponding dynamical system can be reduced to autonomous one. We also define a new variableτ for time through which the corresponding derivative is represented by prime such that dt dτ = 3 Θ . Here each term is associated with some physical explicit background since the chosen dimensionless variables Ω andΨ occur due to physical impact of viscous energy density and pressure, respectively. The system of Eqs. (13) and (14) in terms of these normalized variables takes the form Differentiation of the dimensionless variable for energy density gives Using Eqs. (19) and (20), this equation turns out to be Now we introduce the concept of a new evolution equation forΨ. The first derivative ofΨ with respect toτ through Eq. (19) leads to an evolution equation of the form It is mentioned here that Eqs. (22) and (23) play a remarkable role to describe the respective dynamical system for phase space analysis. In order to find the critical points {Ω c ,Ψ c }, we need to solve the dynamical system by imposing the condition Ω ′ =Ψ ′ = 0. The stability of FRW universe model can be analyzed according to the nature of critical points. Here we restrict the phase space region to a condition which is necessary for the positivity of entropy production rate given by [18,26] This condition tends the possible negative values ofΨ towards zero for k 2 ≫ ν 2 . Contrarily, the bulk pressure will be less restrictive if k 2 ≪ ν 2 . It is noted that finite values of k allow only positive values of bulk pressure in the limit ν → 0. It would be more convenient to consider k 2 ≤ ν 2 along with ν 2 ≤ 2 − γ and τ * = k 2 τ which leads to the fact that the characteristic time for nonlinear effects τ * does not exceed the characteristic time in linear background τ . We characterize the critical points by deceleration parameter q = −1 − Θ ′ Θ and effective EoS parameter γ ef f = − 2Θ ′ 3Θ which yield To examine a region of phase space undergoing accelerated expansion, we impose q < 0 in Eq.(25) which gives The possibility of accelerated expansion in the physical phase space is determined by comparing Eqs. (24) and (25) through q < 0 given by Substituting Ω ′ = 0 in Eq.(22), we find the following conditions We insert these conditions inΨ ′ to find the location of critical points. This analysis is carried out by characterizing the viscous fluid through the choice of its EoS parameter γ (radiation). We consider 0 < k 2 = ν 2 ≤ 2 − γ for which the case of stiff matter (γ = 2) is excluded from the analysis as it yields ν 2 = 0.

Radiation Case (γ = 4 3 )
We consider the radiation case for phase space analysis by taking γ = 4 3 . Imposing the condition (29) andΨ ′ = 0 in Eq.(23), we have This cubic equation yields three roots among which we retain only those roots that lie in the physical phase space. The general form of the dynamical system is given by The eigenvalues of the system can be determined by the Jacobian matrix The eigenvalues for the above stability matrix corresponding to the points P ± r are given by The fixed point is called a source (respectively, a sink) if both eigenvalues consist of positive (respectively, negative) real parts. In case of viscous radiating fluid, we can explore source and sink according to the sign of eigenvalues as well as direction of the trajectories. We investigate two critical points P + r = {1,Ψ + c } and P − r = {1,Ψ − c } corresponding to positive (Ψ + c ) and negative (Ψ − c ) roots, respectively. By taking Ω c = 0 and the second condition (30) withΨ c = − Ωc 3 − 3p EM , we obtain P 0 r = {0, −3p EM }.

Case I:
We are interested to analyze the impact of electromagnetic field on stability of the critical points in the presence of nonlinear bulk viscosity. The energy density (5) and pressure (6) are given by The dynamical behavior of critical points for different values of electric and magnetic fields as well as other parameters is shown in Figures 1-2. The green trajectory depicts a flow from the point P + d towards P − d . The white region corresponds to the negative entropy production rate that diverges on its boundary whereas the green region shows accelerated expansion of the universe (q < 0). Here the point P 0 d shows varying behavior, i.e., either it is a saddle point or a sink depending on the values of different parameters.
In these plots, we have taken ζ 0 = 0.2, 1 by varying ν, k, B and E. For ν = k = 1/5 and ζ 0 = 0.2, it is found that the global attractor P − d lies in green region showing accelerated expansion for the same values of B and E. This region tends to decrease by increasing E such that the point P − d lies in the deceleration region. By increasing ζ 0 , we find accelerated expansion with different values of B, E and larger values of the parameters ν and k. For ν = k = 1 and ζ 0 = 1, we find accelerated expansion of the universe model for all choices of electric and magnetic fields. The point P 0 d behaves as a sink for ζ 0 = 0.2 which becomes a saddle point for larger values of ζ 0 . We observe that the increasing value of bulk viscosity increases the region for accelerated expansion in the presence of NLED. In the following, we discuss two different cases for electric as well as magnetic universe.

Case II (E = 0):
It is well-known that NLED helps to diminish the initial singularity in the early universe where only the primordial plasma identifies matter [27]. Some recent results indicate that a magnetic universe is appropriate to avoid the initial singularity and ultimately shows late time accelerated expansion [28].
Here we assume the squared electric field < E 2 > to be zero such that the magnetic field (F = 2B 2 ) rules over the universe known as magnetized universe. Thus the energy density (5) and pressure (6) take the form The respective evolution plots are given in Figure 3. For ν = k = 2/3 and ζ 0 = 1, we find that sink lies in green region showing the stability of accelerated expansion for the magnetized universe. This region tends to decrease by increasing the value of magnetic field B. The point P 0 d behaves as saddle for small values of magnetic field. It is mentioned here that increasing values of bulk the viscosity and the parameters ν as well as k with different values of B give rise to the stability of accelerated expansion of the universe for different choices of B. We also find that a smaller value of bulk viscosity show decelerated expansion with increasing values of B. Here, we deal with the electric universe by setting < B 2 >= 0. The corresponding energy density and pressure are given by The plots corresponding to different choices of electric field E are shown in Figure 4. For ν = k = 2/3 and ζ 0 = 1, we analyze the sink P − d in green region showing accelerated expansion of the universe for different values of E. The point P 0 d behaves as a saddle for small values of magnetic field. We find that the region for accelerated expansion tend to decrease by increasing electric field E. It is observed that accelerated expanding region exists for increasing values of bulk viscosity and parameters ν as well as k with all choices of E. It supports the fact that the role of bulk viscosity and electric field is to increase the stability of accelerated expansion of the universe model. The summary of our results filled with viscous radiating fluid is given in Table 1.

Power-Law Scale Factor
In this section, we discuss the power-law behavior of scale factor corresponding to the critical points. For this purpose, we integrate Eq. (19) For Θ = 0, we formulate power-law scale factor whenever 1 + 3p EM + (γ − 1)Ω +Ψ = 0. Solving Θ = 3ȧ a for a(t), we obtain the generic critical point as The following condition must hold for exponentially expanding models (identified by the condition 1 + 3p EM + (γ − 1)Ω +Ψ = 0) to be present in the physical phase space region (bounded by Eq. (24)) This condition is not satisfied in the physical phase space for ν 2 = k 2 . If ν 2 > k 2 , the above inequality must be satisfied in the following physical phase space region It is mentioned here that sign of the term 1 + 3p EM + (γ − 1)Ω +Ψ is quite important to evaluate different cosmological stages. If 1 + 3p EM + (γ − 1)Ω + Ψ = 0, it corresponds to the exponential expansion of the universe model. Also, 1+3p EM +(γ −1)Ω+Ψ ≷ 0 yields accelerated expansion or contraction of the cosmological model, respectively. If ν 2 < k 2 , the possibility of having accelerated expansion will narrow down. Figure 5 shows the physical phase space region (excluding the white region with negative entropy production rate) whereas yellow and dark gray regions correspond to accelerated and exponential expansion of the universe model for v 2 > k 2 , respectively. Table  2 provides the polynomial behavior of power-law scale factor for different critical points with 1 + 3p EM + (γ − 1)Ω +Ψ = 0.

Outlook
In this paper, we have discussed the impact of NLED on the phase space analysis of isotropic and homogeneous universe model by taking noninteracting mixture of electromagnetic and viscous radiating fluids. This analysis has been proved to be a remarkable technique for the stability of dynamical system. An autonomous system of equations has been developed by defining normalized dimensionless variables. We have evaluated the corresponding critical points for different values of the parameters to discuss stability of the system. We have also calculated eigenvalues which characterize these critical points. We summarize our results as follows. Firstly, we have discussed stability of critical points through their eigenvalues corresponding to different values of E and B for viscous radiation dominated universe model. It is found that the critical points P + d and P − d correspond to source (unstable) and sink (stable), respectively (Figures 1-2). It is mentioned here that the green region corresponds to accelerated expansion of the universe. The point P − d is a global attractor in the physical phase space region which leads to an expanding model dominated by viscous matter for various choices of cosmological parameters. In the presence of both electric and magnetic fields, we find that bulk viscosity increases the region for accelerated expansion while the increasing values of E shows deceleration region for smaller values of bulk viscosity. It is mentioned here that large values of bulk viscosity as well as other parameters correspond to accelerated expansion of the respective universe model for all choices of electric and magnetic fields.
We have also studied electric and magnetic universe cases separately. It is found that sink lies in the green region showing accelerated expansion of the magnetized universe for smaller values of bulk viscosity and other parameters while increasing value of magnetic field decreases this region. For B = 0, we have analyzed accelerated expansion of the universe model corresponding to large values of the parameters which tends to decrease by increasing E. It is worth mentioning here that the role of bulk viscosity is to increase the green region for accelerated expansion with different choices of E and B for both electric as well as magnetic universe. Moreover, we have also studied the behavior of power-law scale factor corresponding to the critical points. It is found that the power-law scale factor indicates various phases of evolution (accelerated or exponential expansion) of the respective universe model.