Falsifying cosmological models based on a non-linear electrodynamics

Recently, the nonlinear electrodynamics (NED) has been gaining attention to generate primordial magnetic fields in the Universe and also to resolve singularity problems. Moreover, recent works have shown the crucial role of the NED on the inflation. This paper provides a new approach based on a new model of NED as a source of gravitation to remove the cosmic singularity at the big bang and explain the cosmic acceleration during the inflation era on the background of stochastic magnetic field. Also, we found a realization of a cyclic Universe, free of initial singularity, due to the proposed NED energy density. In addition, we explore whether a NED field without or with matter can be the origin of the late-time acceleration. For this we obtain explicit equations for $H(z)$ and perform a MCMC analysis to constrain the NED parameters by using $31$ observational Hubble data (OHD) obtained from cosmic chronometers covering the redshift range $0<z<1.97$; and with the joint-light-analysis (JLA) SNIa compilation consisting in $740$ data points in the range $0.01<z<1.2$. All our constraints on the current magnetic field give $B_{0}\sim 10^{-31}\mathrm{cm^{-1}}$, which are larger than the upper limit $10^{-33}\mathrm{cm^{-1}}$ by the Planck satellite implying that NED cosmologies could not be suitable to explain the Universe late-time dynamics. However, the current data is able to falsify the scenario at late times. Indeed, one is able to reconstruct the deceleration parameter $q(z)$ using the best fit values of the parameters obtained from OHD and SNIa data sets. If the matter component is not included, the data sets predict an accelerated phase in the early Universe, but a non accelerated Universe is preferred in the current epoch...


I. INTRODUCTION
Universe started with a Big Bang, which had a singularity that all the laws of physics would have broken down [1,2]. Today, the first thing that is known that the two great theories of physics such as quantum mechanics and general relativity mostly work very well, except in some extreme conditions like the Big Bang [3]. There is something which is clearly still missing to explain this singularity [4]. Recently, cosmological models using non-linear electromagnetic fields (NEF) have been gain interest to remove singularity problem of the Universe at the Big Bang and also singularities of curvature invariants [5][6][7][8][9]. The Standard Cosmological Model (SCM) based on Friedmann-Robertson-Walker (FRW) geometry, also known as ΛCDM model, does not produce any solution for the singularity problem at the beginning of the Universe [10]. The SCM has a problem of singularities. If the Maxwell equations are intelligently modified, these singularities can be resolved. There are some cosmological models known as magnetic Universe that has no singularity because of the nonlinear modification of the Maxwell electrodynamics at strong fields such as early universe [11,12], because of the background of the conformally flat Robertson-Walker metric, but the cost is the break up of the conformal invariance of Maxwell theory [12][13][14][15][16][17][18][19][20][21].
Moreover, today it is widely accepted fact among the physicist that the Universe is accelerating. The idea of an accelerating Universe is supported and confirmed by type Ia supernovae and the cosmic microwave background (CMB) [1,10,22,23]. However the reason of the acceleration of the Universe is not entirely clear, nevertheless a number of solutions have been proposed [8,10,. One such solution is to introduce the cosmological constant in the Einstein's field equations. In this scenario, the acceleration of Universe is driven by dark energy (DE) which can be thought of as a kind of space-filling fluid with constant energy density through the Universe [10,51]. Another exotic form of matter proposed as a DE candidate is to consider a scalar field known as the quintessence [52]. On the other hand, the acceleration rate of the Universe in terms of the modified gravity theories continue to attract interest [51]. The simplest model which generalizes General Relativity is found by simply replacing the Ricci scalar (R) in the action by a function f (R). This idea led to many modified gravity models studied in the literature [53] or in the cosmological set up of a higher-order modified teleparallel theory (see [54][55][56][57][58][59][60][61][62][63] and references therein). Instead of doing modification of gravity, the nonlinear electrodynamics (NED) can be used to avoid singularities as well as resolve the horizon problem [12].
Recently, the idea of NED has been proposed as a solution to source the Universe acceleration [5][6][7]. In the early Universe the effect of the NED may have been very strong and, in principle, this may also explain the inflation. In this scenario the NEF can be considered as a source of the gravitational field and, as a consequence, nonlinear magnetic fields may be a driven mechanism of the inflation of the Universe. In this line of research, very recently many NED models have been investigated using a stochastic magnetic background, i. e. the cosmic background with the wavelength smaller than the curvature, with a non-vanishing < B 2 >, where matter should be identified with a primordial plasma [5, 6, 8, 12-14, 16, 64]. Thermal fluctuations in a dissipative plasma, i.e. plasma fluctuations, could source stochastic magnetic fields on a scale larger than the thermal wavelength [75]. Thus, there are the stochastic fluctuations of the electromagnetic field in a relativistic electron-positron plasma. Although homogeneous magnetic fields can affect the isotropy of the Universe, i.e. the energy-momentum tensor can become anisotropic which could cause an anisotropic expansion law and modify the CMB spectrum, the effect on the Universe geometry (isotropy) of magnetic fields tangled on scales much smaller than the Hubble radius are negligible [75]. Thus, averaging the magnetic fields, which are sources in general relativity [79], give the isotropy of the Friedman-Robertson-Walker (FRW) spacetime. On the other hand, bulk viscosity term is neglected in the electric conductivity of the primordial plasma by taking E 2 = 0 [65][66][67][68][69]. The NED is useful to remove singularities of Big Bang and try to explain inflation naturally. Note that using the scalar fields for the inflation and early Universe have a problem with the fine exit. The problem is that after the inflation is started, it goes forever which is known as eternal inflation. However, we will propose a model where there is not any eternal inflation problem. For this purpose, we use the magnetic Universe with Nonlinear Electromagnetic Field (NEF) in the stochastic background with a nonzero value of < B 2 > that supports the acceleration of the Universe. There are some cons and pros of NED from the Dirac-Born-Infeld (DBI) nonlinear electrodynamics such that for DBI theory, there is a duality symmetry, on the other hand, for the NED, it is broken but it is also violated for the QED. Other difference is the birefringence phenomenon that occur in NED and also in QED with quantum corrections, but not in DBI theory [14]. Moreover, the DBI model has a problem of causality.
It is believed that magnetic fields have important role on the evolution of our universe. However, we have known little information about the existence of magnetic fields at the early universe. Observations have manifested the existence of magnetic field in the universe, ranging from the stellar scale 10 −5 pc to the cosmological scale 10 4 Mpc [70,71]. In particular, the magnetic field on large scales (≤ 1 Mpc) is deemed to be formed in the early universe [72], namely, the primordial magnetic field. By the recent CMB observations, the strength of the magnetic field is less than a few nano-gauss at the 1 Mpc scale [73]. Additionally, the γ-ray detections of the distant blazars imply that the magnetic field might be larger than 10 −16 G on the scales 1 − 10 4 Mpc [74].
On the other hand, we have no direct observational evidence of primordial magnetic fields. The amplitude of primordial magnetic fields is also debatable. However, we believe that they existed because may have been needed to seed the large magnetic fields observed today. Nowadays, many theories are proposed to obtain the origin of cosmic magnetic fields for instance primordial vorticity plasma (vortical motion during the radiation era of the early Universe, vortical thermal background by macroscopic parity-violating currents), quantum-chromo-dynamics phase transition, first-order electroweak phase transition via a dynamo mechanism, etc ( [75] and references therein). In the last scenario, seed fields are provided by random magnetic field fluctuations which are always present on a scale of the order of a thermal wavelength. Primordial Nucleosynthesis limit the intensity of the magnetic seed fields to a current upper limit of 10 −9 G [76] and the lower limit B > 10 −19 G from the γ-ray observations [77]. Today, magnetic fields have been observed in different types of galaxies and also cluster of galaxies at wide range of redshifts. Furthermore, a lower bound, B ≥ 3 × 10 −16 G, has been obtained for intergalactic magnetic fields [78]. Constraints from the Planck satellite in 2015 show that the upper limit to be of the order of B < 10 −9 G [79].
In late time epochs, the reason to use NED may be different than the early universe: it can be implemented as a phenomenological approach, in which the cosmic substratum is modeled as a material media with electric permeability and magnetic susceptibility that depend in nonlinear way on the fields [126]. Another argument is based on the view that General Relativity is a low energy quantum effective field theory of gravity, provided that the Einstein-Hilbert classical action is augmented by the additional terms required by the trace anomaly characteristic of NED [127].
In this paper we use a new model of Lagrangian of NED, dubbed "NED with an exponential correction", which has a Maxwell limit at low energies and which is different than DBI nonlinear electrodynamics [5]. We assume that radiation of NED is dominated in the early Universe, to solve the initial singularity problem. We show that the NED with the gravitation field can create the negative pressure and cause the inflation and late cosmic acceleration of the Universe. In this regard, our manuscript is an extended version of the papers [5,8,14,20,21,80], at which are used Nonlinear magnetic fields, as a source of inflation. This model of NED is valid for the early and current regime of the Universe, on the other hand for the late Universe, the magnetic field is very weak. However, the question whether one can explain the late regime of the Universe in terms of only NED remains, at least theoretically, as an open possibility. Hence, we additionally investigate whether the exponential NED can provide the late-time accelerated expansion of the Universe.
We perform a phase-space analysis of Einstein-NED cosmology without a matter source and including it. The advantage of the using phase-space analysis [41][42][43][44][45][46][47][48][49][81][82][83][84][85][86][87][88][89][90][91] is that one can do more stability analysis with using visual plots using the trajectories in geometrical way so that it becomes easy to observe the property with the help of the attractors which are the most easily seen experimentally [92]. On the other hand, conceptually using NED has the advantage that no need to use some exotic fields such as scalar fields, branes or extra dimensions, it is just photon fields, and it is well known also in nonlinear optics which studying behavior of light in nonlinear media and also nonlinear collision of particles in quantum electrodynamics [93][94][95][96][97][98][99].
The paper is organized as follows. In section II we introduce the Lagrangian of NED. In section III we examine the acceleration and evolution of the Universe in terms of NED fields. In section IV we perform a detailed phase-space analysis of our Einstein-NED cosmology model. In section V we consider a more realistic scenario, namely we include a matter source in our setup. In section VI we put constraints on the NED parameters using observational Hubble data from cosmic chronometers and using the latest SNIa data. Finally we summarize and discuss our results in section VII.

II. GENERAL RELATIVITY COUPLED WITH NON-LINEAR ELECTRODYNAMICS WITH AN EXPONENTIAL CORRECTION
In highly nonlinear energy density situations such as in the early Universe, NED is expected to play a crucial role in the evolution of the Universe [11,19]. First it should be understand the contributions of nonlinear fields to inflation. For this purpose, we propose the following action GR coupled with the NED field as follows: where R is the Ricci scalar and L N ED is the Lagrangian of the NED fields and we have used units where 8πG = 1, c = 1. The new NED Lagrangian density is chosen as follows: where α is a constant with [B −2 0 ] units, and B 0 is the current value of the electromagnetic magnetic field, and β is a dimensionless parameter, After varying the action given in Eq. (1), one can find the Einstein field equations and the NED fields equation as follows: where The energy momentum tensor [7] T µν = K µλ F λ ν − g µν L N ED , K µλ = can be used to obtain the general form of the energy den-sity ρ N ED and the pressure p N ED of NED fields as and To find the solution of the Einstein field equations, we consider the homogeneous and isotropic cosmological metric of FRW with following line element: with scale factor a(t).
The key point of this study is that it can be supposed that the stochastic magnetic fields are the cosmic background with the wavelength smaller than the curvature so we can use the averaging of EM fields which are sources in GR and then we can obtain the isotropic FRW spacetime [101]. In general, the averaged EM fields have the properties: where the averaging brackets is used for a simplicity. The non-zero averaged magnetic field case is the most unexpected case [101], where the magnetic field of the Universe is frozen to occur the magnetic properties, it is necessary to screen the electric field of the charged primordial plasma. We use the Eqs. (6) and (7) (for E 2 = 0) and obtain: and (11) Through the analysis, it is imposed the energy condition ρ N ED ≥ 0, which implies 2β + αB 2 ≥ 0.
Afterwards, we use the FRW metric (8) and find the equation of the Friedmann: where we use the dot "." to denote the time derivative. Furthermore, we check the condition of the accelerated Universe ρ N ED + 3p N ED < 0 with the sources of NED fields using into Eqs. (6) and (7): FIG. 1: Are presented the regions which gives acceleration in the plane Y = αB vs. β.
The condition for the NED energy provide acceleration is (ρ N ED + 3p N ED < 0), which combined with the physical condition ρ N ED ≥ 0 gives β ≤ 0, α B 2 > −2β or β > 0, α B 2 > − 1 2 (1 + 2β) + 1 2 4β 2 + 12β + 1. As shown in Fig. (1) and the lower bound of αB that gives acceleration when β = 1 is at α B 2 = − 3 2 + 1 2 √ 17. Note that the source of the strong nonlinear electrodynamics field accelerates the Universe in the early stages, that is, for large αB 2 we enter the shadowed region in Fig. 1. Then we use the conservation of the energy-momentum tensor (∇ µ T µν = 0) for the FRW metric (8) and find the continuity equation: It is noted that the Hubble parameter is defined as H = a a which is the expansion rate of our Universe. This equation gives (15) One can integrate (one of the branches of) the above equation 1 using the ρ and p to obtain the evolution of the magnetic field respect to the scale factor where B 0 is for a(t) = 1. Afterward, we rewrite the energy density ρ and the pressure p using the evolution of the magnetic field: Using the equation of state (EoS) then we use Eqs. (16) and (17) and obtain the radiation or other relativistic fluid case: . Now, by choosing the value 1 We have also the trivial solution B = 0, the regime αB 2 → ∞ and the constant solutions B, such that −4β + α 2 B 4 + 2αβB 2 = 0. We submit the reader to section IV A for a more complete discussion of this special cases. , we have ω = 0 for non-relativistic matter (baryons, CDM), also for some values we have −1 < ω < −1/3 which is called "quintessence", a dynamical dark energy resulting in accelerating expansion. The models with ω < −1 has been termed "phantom" dark energy. We can also refer dark energy as vacuum energy, because one possible source of dark energy is quantum-mechanical fluctuations of the vacuum. In Fig. 2 is it plotted ω versus Y = αB 2 0 /a 4 . For small a (large Y ) the universe has large negative equation of state, it crosses the value 17 − 3 the universe changes from acceleration to deceleration. For large a (Y → 0) we see that ω approaches 1/3 (radiation dominated universe).
Now we show that the spacetime will be flat at t → ∞ (a → ∞) and that singularities are removed at the early/late phase of the Universe (a → 0). For this purpose we calculate the Ricci scalar (which represents the curvature of spacetime), the Ricci tensor squared, and the Kretschmann scalar. The Ricci scalar is calculated by using Einstein's field equation (3) and the energy-momentum tensor as follows: The Ricci tensor squared R µν R µν and the Kretschmann scalar R µναβ R µναβ are also obtained as and Furthermore, taking the limits a → 0 and at a → ∞ in (16) and in (17) Finally, when we check the limits of the energy density and pressure, and use the expressions (20), (21) and (22) we conclude that the Ricci scalar, the Ricci tensor squared, and the Kretschmann scalar are non singular at a(t) → 0 and at a(t) → ∞, and these show that the spacetime will be flat at t → ∞ and singularities are removed at the early/late phase of the Universe.

III. ACCELERATION AND EVOLUTION OF THE UNIVERSE
Now, we find the evolution of the Universe using the Einstein's equations and the energy density given in Eqs. (16) and (17). Without considering dust like matter, we study the evolution of the Universe. To calculate the scale factor as a function of time, first we use the second Friedmann's equation in flat Universe: and one can obtain the equation which shows conservation of energy for a particle moving in a effective potential V (a):ȧ We use equation (26), to get a qualitative feel for the evolution of the early Universe. Figure 3 shows the effective potential V eff (a) as function of the scale factor. The effective potential with a positive slope yields a force tending to slow down positive motion along the horizontal axis, while the portion of the effective potential with a negative slope yields a force tending to speed up positive motion along the horizontal axis. These two conditions occur to the right and the left, respectively, of the minimum scale factor at By analogy, then, a(t) accelerates to the left of a c and decelerates to the right of a c . This acceleration is due to nonlinear electrodynamics which behaves similarly to dark energy.
Using the Eq. (25) with the energy density in Eq. (16), we obtain the equation as follows: Expanding (27) 3 . Neglecting error terms, we calculate approximate cosmic time [102] as follows: where t 0 is a constant of integration which gives only the shift in time as shown in Fig. (4). Note that without losing generality, the integration constant has been shifted to t 0 = 0. Then, by simplicity we have imposed the condition β = 1. One can also obtain the evaluation of the scale factor for (α = B 0 = 1) as follows: Note that evolution of the scale factor has similar feature with radiation dominated Universe and for t = 0 it reduces to zero The function of a 0 is a radius of the Universe which shows that Universe begins from the zero point. One can also show them in terms of redshift where a = (1 + z) −1 . Now, looking at large a in figure 3 we see that V eff tends to zero as a → ∞ (that, is, when B = B 0 a −2 goes to zero); In figure 4, a goes to ∞.ȧ → B0 To complete the analysis we expand around large B, that is, small a, we haveȧ 2 = e − αB . Integrating we ob- that tends to (0, −∞) in the limit a → 0. This solution corresponds to the shallow region of the potential showed in (3) near the origin a = 0, whereȧ = −V eff → 0, and it is achieved in the early universe. The next step is to check the causality of the Universe using the speed of the sound. The speed of the sound should be smaller than the local light speed (c s ≤ 1) [103]. Moreover, we check the square sound speed to avoid the Laplacian instability, we require the conditions that must be positive (c 2 s > 0).Furthermore there exists a range in the parameter space where it is violated the condition for the absence of Laplace instability, "forbidden region". Now our model is free of Laplacian instability. Our all calculations are based on allowed region where there is no Laplacian instability.
If the Universe satisfied these conditions, there is a classical stability. The square of the sound speed is found from Eq.s (6) and (7); and the classical stability (c 2 s > 0) for the absence of Laplacian instabilities and causality of the Universe using the speed of sound c s ≤ 1, occurs as shown in Fig. 5 and Fig. 6, respectively. Finally, we deal with the configuration F = −β/α, which corresponds to 2β + αB 2 = 0. Using the parametrization B = B 0 a −2 , and the definition Y = αB 2 0 a −4 , this configuration appears when Y = −2β (in Fig 2. β has been chosen equal to 1), and it is attained as a → a s := (α negative). Expanding eq. (27) in a neighborhood of the value a = a s we haveȧ = The Hubble factor as t → 0 goes as H →  + O t 5/3 , such that we obtain correction terms to the previous asymptotic formula as t → 0. On the other hand, for large t we have H → 2 3 t −1 , which corresponds to matter domination. The scale factor around a s corresponds to the inflationary solution a(t) = exp[At f ] with f = 2 3 , that is followed with matter domination. A scale factor of the form a(t) = exp At f where A > 0 and 0 < f < 1 was introduced in [104][105][106] in the context of inflation. Since the expansion of the universe with this scale factor is slower than the de Sitter inflation (a(t) = exp(Ht) where H is constant), but faster than the power-law inflation (a(t) = t q where q > 1), it was called intermediate inflation. Intermediate inflationary models arise in the standard inflationary framework as exact cosmological solutions in the slowroll approximation to potentials that decay with inverse power-law of the inflaton field [107]. These models have been studied in some warm inflationary scenarios [108][109][110][111][112][113][114][115][116][117][118]. From a quantum mechanical point of view, this mechanism might provide an explanation for the large scale magnetic fields observed today. In particular the inflation period amplifies the quantum perturbations of the electromagnetic field leading to the current classical perturbations.

IV. PHASE SPACE ANALYSIS
The equation (26) represents the motion of a particle of the unit mass in the potential U eff (a) = − ρ N ED (a)a 2 6 .
We have that the equation (26) is satisfied on the zero energy level, where ρ N ED plays the role of effective energy density parameterized through the scale factor a(t). Therefore the standard cosmological model can be simply represented in the terms of a dynamical system of a Newtonian type: where the scale factor a plays the role of a positional variable of a fictitious particle of the unit mass, which mimics the expansion of the Universe. For simplicity, we use introduce the new constant b 0 = 1 2 B 0 2 α, and introduce the time rescaling τ = t B0 , B 0 > 0, and the variables The system is then equivalent to By definition x ≥ 0 (since we consider that the scale factor is non-negative). For this reason, it is convenient to define the new variables which takes values on the real line, and the time variable This system can be written in the form Thus, v 2 2 + W (u) = E, is the constant of energy. From the above system we see that, generically, the fixed points are situated on the axis u (v = 0). From the characteristic equation it follows that just three types of fixed points are admitted: 1. saddle if u 0 : ∂W ∂u | u=u0 = 0 and ∂ 2 W ∂u 2 | u=u0 < 0; 2. focus if u 0 : ∂W ∂u | u=u0 = 0 and ∂ 2 W ∂u 2 | u=u0 > 0; 3. degenerated critical point if u 0 : ∂W ∂u | u=u0 = 0 and In the concrete example we have Imposing the condition b 0 + βe 4u ≥ 0 we find that the fixed point at the finite region of the phase space is A : So that, according to the previous classification it is a focus. For the particular case b 0 = 0.5, β = 1 we find that the fixed point A has coordinates u = 0.144262 for which a = 1.15519 which corresponds to the minimum of V (a) in Fig 3. Additionally, in the limit u → +∞, we have W (u) → 0, W (u) → 0, hence, in this regime we can expect to have degenerate critical points. Furthermore, for b 0 β < 0 there are no fixed points at the finite region of the phase-plane.
Since the above system is in general unbounded, then we introduce the compactification we have the equivalent flow For the choice β = 2b 2 0 +b0 1−2b0 , 0 < b 0 < 1 2 which implies β > 0, b 0 > 0, the fixed point A exists. In Fig. 7 it is presented the dynamics of system (39)-(40) for the choice β = 2b 2 0 +b0 1−2b0 , b 0 = 0.06. In this case A is a focus according to the previous classification. Now, for =: k 2 > 0. Then, the equations near the origin of coordinates (fixed point A), can be written as where we have neglected the higher order terms, such that, close to the origin, the solutions can be approximated by V (η) = −kU 0 sin(ηk) + V 0 cos(ηk).
Thus, the orbits of the original system near the origin can be approximated by the ellipses given by which leads to periodic solutions: where U and V are given by (42), (43) The relation between η and the cosmic time t can be found by using which in general must be integrated numerically. In the figure 8 it is shown the behavior of the scale factor and its first derivative in terms of the parameter η. This is the realization of a cyclic Universe in our model supported by NED. The scale factor goes below and above the value a = 1, reaching a maximum and a minimum value of a, and a is bounded away zero (there is no initial singularity, as expected from our NED proposal). Now let us assume arbitrary b 0 > 0, β > 0 and obtain a second order expansion around A. For the sake of sim-plicity of notation we define µ = By neglecting higher orders terms we obtain Taking as initial condition , V (0) = δ V , and under the conditions b 0 > 0, β > 0, the solution is given by This solution approximates the exact solutions of the full system surrounding A. The orbits near A can be approximated by the ellipses with center A. Once we know the expressions of U (η), V (η), we can calculate the parametric expressions of a andȧ as functions of η through where U and V are defined by (49) and (50), respectively. Finally, the dynamics at the circle at infinity can be represented by the flow of The orbits lying on the circle at infinity can be parametrized as Thus, there are fixed points lying on the circumference at infinity (U, V ) = (±1, 0). Thus, apart of the singular points at infinity described before, and the point at the finite region (U, V ) =

A. Integrability and connection with the observables
In this section we comment on the integrability of the system at hand, and calculate some observables in terms of redshift.
As we commented before, from the conservation of the energy-momentum tensor (∇ µ T µν = 0) for the FRW metric we have and from (10) we have As we mentioned before, we have the trivial solution B = 0, the regime αB 2 → ∞ and the constant solutions B, such that −4β + α 2 B 4 + 2αβB 2 = 0. Concerning the former solutions we have the following results: β + β(β + 4) + 2 6α , but under the above conditions results in Λ < 0, which implies a ∝ e i √ |Λ|t . So, we discard these solutions. Furthermore, we have a second class of solutions given by 2. αB 2 = −β + β 2 + 4β is a real value for β > 0, α > 0 or β ≤ −4, α > 0. As before, equation (55) means that H is constant and equal to a cos- For β ≤ −4, α > 0 we have Λ < 0, which implies a ∝ e i √ |Λ|t . So, we discard these solutions. However, for β > 0, α > 0, we obtain Λ > 0 as required in, and then we arrive a a regime similar to de-Sitter type Universe, with a ∝ e √ Λt . Summarizing, our model supports the inflation similar to de-Sitter Universe, with provided β > 0, α > 0. Now, assuming that B is not a constant or trivial, it is easy to recast the field equations aṡ Using the previous restriction and introducing the logarithmic time variable N = ln a, we obtain the equations We use convention that a = 1, is the present value of the scale factor, so that N = 0 today. We define H(0) = H 0 and B(0) = B 0 the current values of the Hubble factor and the magnetic field respectively. It is convenient to define the dimensionless parameters Integrating the system (59), evaluating N = ln a and a = 1 1+z , we obtain where we have taken the reference values z = 0, a = 1 for today, such that, as z → −1, the contribution of magnetic fields at cosmological scales is negligible [119]. The current value of the deceleration parameter is Thus, to accommodate the current accelerated phase it is required h 0 > 1 6 e −b0 (1 − 2b 0 ). Summarizing, the function H(z) has the free parameters H 0 , h 0 , b 0 than can be contrasted against data and experiments. Then, we calculate the model parameters as: Notice that our analysis is performed in geometrized units where the dimension of B 0 and H 0 is [L −1 ]. In geometrized units, 8πG = 1 and c = 1, the conversion factors are 1Gauss = 1.44 × 10 −24 cm −1 , and H 0 = h 1.08 × 10 −30 cm −1 . Then we are dealing with magnetic fields of the order 10 −40 cm −1 B 0 10 −33 cm −1 in the present epoch [78,79].

V. GENERAL RELATIVITY COUPLED WITH NON-LINEAR ELECTRODYNAMICS WITH AN EXPONENTIAL CORRECTION IN THE PRESENCE OF MATTER
To introduce a more realistic model we include an additional matter source with constant equation of state parameter w m = p m /ρ m (with 0 ≤ w m ≤ 1 for standard matter) and with continuity equatioṅ Integrating out the last equation we obtain ρ m = ρ m,0 a −3(1+wm) . In this case, the second Friedmann's equation in flat Universe becomes Hence, the effective potential is now where, as in the previous case, we have introduced the new constants b 0 = 1 2 B 0 2 α and ρ m,0 = ρ 0 B 2 0 . Defining the variables which takes values on the real line, and the new time variable the system is then equivalent to where the effective potential is Now, the fixed points are found by solving numerically As before the above system is in general unbounded, so that we introduce the compactification Hence, we have the equivalent flow For the special choice of parameters e −b 0 (2b 2 0 +2b0β+b0−β) 6(b0+β) 2 − 1 6 ρ 0 (3w m + 1) = 0, the origin is a fixed point of the dynamical system. Given γ := U eff (0) = e −b 0 (b 2 0 (−(8b 2 0 +2b0+(6b0+3)wm+3))+β 2 (2(7−4b0)b0+(3−6b0)wm−1)−4b0β(b0(4b0+3wm−3)−3)) 6(b0+β) 3 .
Depending of whether sign it has, the origin is a saddle (U eff (0)<0) or a focus (U eff (0) > 0).
Assuming γ > 0, and by taking a linear expansion around the origin, we find the approximate system with solution Finally we have the parametric solution which is a periodic solution for γ > 0. For γ < 0, the solution is and has we commented before the origin is a saddle.
On the other hand, in the same way as for (39)-(40), the asymptotic system near the fixed point at infinity for The orbits lying on the circle at infinity can be parametrized as Thus, there are fixed points lying on the circumference at infinity (U, V ) = (±1, 0). Now, we proceed to the numerical integration of system (75)- (76) for some values of b 0 and β, and for a) dust and b) for stiff matter. There is a numerical evidence, as showed in in Fig. 9, of the realization of a cyclic Universe in our model supported by NED. That is, the scale factor goes below and above the value a = 1, reaching a maximum and a minimum value of a, and a is bounded away zero, as expected from our NED proposal. As shown in figure 9 (b) and (e), when A is a saddle, there are other cyclic solutions.

A. Integrability and connection with the observables
In this section we comment on the integrability of the system at hand, and calculate some observables in terms of redshift. We assume that B is not a constant or trivial.
By introducing the logarithmic time variable N = ln a, it is easy to recast the field equations as As before a = 1 corresponds to present value of the scale factor, so that N = 0 today and H(0) = H0, B(0) = B0 and ρm(0) = 3H 2 0 Ωm0 are the current values of the Hubble factor, the magnetic field, and the matter density, respectively, where we have defined the current normalized energy density Ωm(0) =: ρm 3H 2 | today = Ωm0. We use the dimensionless parameters h0 = αH 2 0 , b0 = αB 2 0 2 .
Integrating the above system and evaluating N = ln a and a = 1 1+z , we obtain Summarizing, the function H(z) has the free parameters H0, Ωm0, h0, b0 than can be contrasted against data and experiments. Then, we calculate the model parameters as: and the current value of q is

VI. OBSERVATIONAL CONSTRAINTS
In this section we test two NED models: without a fluid component and including it as dust matter (associated to a dark matter component), i.e. its equation of state is wm = 0. The free parameters H0, h0, b0, and Ωm0 (when dust matter is considered) are constrained using the observational Hubble data from cosmic chronometers and the latest SNIa data. Next we calculate the model parameters using Eqs. (65) and (92).
• Observational Hubble data (OHD). The differential age (DA) method measures H(z) between two passivelyevolving galaxies with similar metallicities and separated by a small redshift interval (cosmic chronometers)    [120,121]. The data provided by the DA method are cosmological-model-independent and then they can be used to probe alternative cosmological models. Here, we use the latest OHD obtained from DA technique, which contains 31 data points covering 0 < z < 1.97, compiled by [122] and references therein. The figure-ofmerit for the OHD is written as where H(zi) is the theoretical Hubble parameter (given by either Eq. 62 or Eq. 90, depending on whether we consider no matter or we include it, respectively), H da (zi) is the observational one at redshift zi, and σH i its error. Notice that in the last expression, we also consider the measurement of H0 = 73.24 ± 1.74Kms −1 Mpc −1 by [100] as a Gaussian prior.
• Type Ia Supernovae. The first evidence of the accelerating expansion of the Universe was provided by the observations of distant type Ia Supernovae (SNIa) [22,23]. Over the last years, the high-resolution SNIa observations have demonstrated be a key cosmological probe due the shape their light curves can be standardizable. Thus, any alternative cosmological model should be confronted with the latest SNIa data. Here, we use the joint-light-analysis (JLA) compilation by [123] consisting in 740 data points in the range 0.01 < z < 1.2. For the JLA sample, the observational distance modulus can be computed as where mB is the observed peak magnitude in rest-frame B band, X1 is the time stretching of the light-curve, C is the supernova color at maximum brightness, and MB 2 , a, and b are nuisance parameters in the distance estimate. On the other hand, the theoretical distance modulus is given by µ th = 5 log 10 (dL/10pc), where dL, the luminosity distance predicted by the NED cosmology (an interesting study of the SNIa luminosity distance in Born-Infeld NED cosmology is presented in [125]), reads as Therefore, the figure-of-merit for the SNIa data can be written as 2 Notice that if the host stellar mass M stellar < 10 10 M where Cov(a, b) is the covariance matrix 3 of µ obs provided by [123].
To constrain the free parameters, (Ωm0, h, h0, and b0), we perform a Bayesian Monte Carlo Markov Chain (MCMC) analysis employing the emcee python module using 800 walkers, 500 steps in the burn-in-phase, and 4000 MCMC steps to guaranty the convergence. We consider a Gaussian prior on Ωm0 as measured by [124] and uniform priors [0, 1] for both h0 and b0. Figures 10 and 11 show the 1D marginalized posterior distribution and the 2D confidence contours for the h, h0, and b0 parameters for the NED cosmology without matter and Ωm0, h, h0, and b0 when the matter component is included respectively. Notice that both OHD and SNIa data provide consistent constraints on the NED parameters for both models. Table I gives the mean values for the NED model parameters using different cosmological data for both cases: without and including a matter component. The chi-square values indicate a over-fitting to OHD but a good-fitting to SNIa data. All our constraints on the current magnetic field, B0 ∼ 10 −31 cm −1 , are larger than the upper limit 10 −33 cm −1 by the Planck satellite implying that NED cosmologies could not be suitable to explain the Universe dynamics at late times. Nevertheless, the Figure 12 illustrates a good fitting of the NED cosmology without (top panel) and with (bottom panel) matter to OHD. In addition, the Figures 13 and 14 show the reconstructed deceleration parameter q(z) in the range 0 < z < 2 for each data set without matter and including matter respectively. If the matter component is not included, although the data sets predict an accelerated phase in the early Universe, a non accelerated Universe is preferred in the current epoch. However, a late cosmic acceleration dynamics is allowed within the 2σ confidence levels. When a matter component is included in the NED cosmology, the data set predict a q(z) dynamics similar to that of the standard model. Indeed, q(z) → 1/2 when z → ∞. Moreover, both cosmological data favor up to 2σ confidence levels an accelerating expansion in the current epoch, i.e., the Universe passes of a decelerated phase to an accelerated stage at redshift ∼ 0.6 for the OHD (top panel) and SNIa (bottom panel) constraints. Therefore, although the NED cosmology including dust matter (wm = 0) predict a higher B0 value, it is able to drive a late-time cosmic acceleration which is consistent with the Y −β allowed regions of the Figure 1.

VII. CONCLUSION
In this paper, we have considered a new field of NED as a source of gravity to shed light on the dynamics behind the accelerating Universe and solve the singularity problem of the Big Bang and the Universe curvature.
We have argued that the NED on cosmological scales could reveal the acceleration of the Universe during the inflationary era. We have also shown that, after this period of cosmic inflation the Universe undergoes decelerated expansion and asymptotically approaches the Minkowski spacetime. From the dynamical system approach we have found for parameters β = 2b 2 0 +b 0 1−2b 0 , b0 = 0.06, B0 = 0.5, U0 = V0 = 0.2, and assuming no background matter, the realization of a cyclic Universe in our model supported by NED near the value a0 = 1. The scale factor goes below and above the value a = 1, reaching a maximum and a minimum value of a, and a is bounded away zero (there is no initial singularity, as expected from our NED proposal). For our two models (not including matter, and for prefect fluid), We have found approximated solutions in parametric form for the scalar factor and its time derivative, which are valid in the neighborhood of the fixed point at the finite region. For some values of the parameters we found also the realization of a cyclic Universe.
For the model without matter, we have found solutions with B constant (αB 2 = −β + β 2 + 4β), which supports the inflation similar to de-Sitter Universe, with a ∝ e The top and bottom panels show the reconstruction of the deceleration parameter q(z) using the h, h 0 , and b 0 mean values from OHD and SNIa data respectively into the theoretical q(z) given by (63). The shadowed areas represent the 1σ and 2σ confidence levels. The inset shows the q(z) behavior in the range 0 < z < 0.2.
In the case of evolving H, it was provided explicit expressions for H(z) by direct integration of the equations of motion in the two models: with and without including a perfect fluid. We have tested these NED models by performing a Bayesian Monte Carlo Markov Chain (MCMC) analysis using OHD and SNIa data. We have found that the first case predicts B0 values higher limits than the Planck satellites bounds. That is B0 ∼ 10 −31 cm −1 , are larger than the upper limit 10 −33 cm −1 . In addition, although NED without matter is able to fit the OHD, it prefers no late cosmic acceleration. For the model including a fluid, we consider dust matter with an equation of state wm = 0 which is associated to a dark matter component. In this particular example we conclude that this is The top and bottom panels show the reconstruction of the deceleration parameter q(z) using using the Ω m0 , h, h 0 , and b 0 mean values from OHD and SNIa data respectively into the theoretical q(z) given by (91). The shadowed areas represent the 1σ and 2σ confidence levels.
a good model for the early time Universe, but there are not statistical differences with the usual model for the radiation epoch. We have reconstructed the deceleration parameter q(z) in the range 0 < z < 2 as shown in Figure 14 for OHD (top panel) and SNIa data (bottom panel). Notice that q(z) → 1/2 when z → ∞. In addition, both cosmological data predict an accelerating expansion, i.e., the Universe passes of a decelerated phase to an accelerated stage at redshift ∼ 0.6. Nevertheless, the B0 estimation is higher than the upper limit by Planck measurements. That is, first we have proved that a model based on NED alone does not passes the late cosmic acceleration test. Second, we have proved that dust matter plus NED passed this test, and could provide late-time acceleration (the combined effect of both matter contributions). Summarizing, using in a combined way the powerful of phase-space analysis, and the observational fit, one is able to falsify (in the sense of test) a theoretical cosmological model based on NED. From one side, dynamical systems tools allows to identify regions in the parameter space to provide stability conditions for fixed points with physical meaning. Furthermore, it was possible to use 2D phase spaces, showing trajectories in a geometrical way, so that it becomes easy to observe the property with the help of the attractors which are the most easily seen experimentally. On the other hand, comparison with data allows to refine more the region of parameters. In principle, although the NED cosmology including dust matter predict a value of B0 greater than the upper bound found by Planck satellite, it is able to drive a late-time cosmic acceleration which is consistent with our dynamical systems analysis and it is preferred by OHD and SNIa data sets. In our opinion this is a very interesting result and deserves discussion.