Next Article in Journal
Quantum Loop Corrections in the Modified Gravity Model of Starobinsky Inflation with Primordial Black Hole Production
Previous Article in Journal / Special Issue
Different Aspects of Entropic Cosmology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Combined Studies Approach to Rule Out Cosmological Models Which Are Based on Nonlinear Electrodynamics

by
Ricardo García-Salcedo
1,2,*,
Isidro Gómez-Vargas
3,
Tame González
4,
Vicent Martinez-Badenes
2 and
Israel Quiros
4
1
Centro de Investigación en Ciencia Aplicada y Tecnología Avanzada, Unidad Legaria del Instituto Politécnico Nacional, Legaria 694, Mexico City 11500, Mexico
2
Escuela Superior de Ingeniería, Ciencia y Tecnología, Universidad Internacional de Valencia—VIU, 46002 Valencia, Spain
3
Department of Astronomy, University of Geneva, 51 Chemin Pegasi, 1290 Versoix, Switzerland
4
Departamento de Ingeniería Civil, División de Ingeniería, Universidad de Guanajuato, Mexico City 36000, Mexico
*
Author to whom correspondence should be addressed.
Universe 2024, 10(9), 353; https://doi.org/10.3390/universe10090353
Submission received: 28 June 2024 / Revised: 12 August 2024 / Accepted: 27 August 2024 / Published: 4 September 2024
(This article belongs to the Special Issue Universe: Feature Papers 2024—'Cosmology')

Abstract

:
We apply a combined study in order to investigate the dynamics of cosmological models incorporating nonlinear electrodynamics (NLED). The study is based on the simultaneous investigation of such fundamental aspects as stability and causality, complemented with a dynamical systems investigation of the involved models, as well as Bayesian inference for parameter estimation. We explore two specific NLED models: the power-law and the rational Lagrangian. We present the theoretical framework of NLED coupled with general relativity, followed by an analysis of the stability and causality of the various NLED Lagrangians. We then perform a detailed dynamical analysis to identify the ranges where these models are stable and causal. Our results show that the power-law Lagrangian model transitions through various cosmological phases, evolving from a Maxwell radiation-dominated state to a matter-dominated state. For the rational Lagrangian model, including the Maxwell term, stable and causal behavior is observed within specific parameter ranges, with critical points indicating the evolutionary pathways of the universe. To validate our theoretical findings, we perform Bayesian parameter estimation using a comprehensive set of observational data, including cosmic chronometers, baryon acoustic oscillation (BAO) measurements, and supernovae type Ia (SNeIa). The estimated parameters for both models align with the expected values for the current universe, particularly the matter density Ω m and the Hubble parameter h. However, the parameters of the models are not tightly constrained within the prior ranges. Our combined studies approach rules out the mentioned models as an appropriate description of the cosmos. Our results highlight the need for further refinement and exploration of NLED-based cosmological models to fully integrate them into the standard cosmological framework.

1. Introduction

The present paradigm of cosmology includes an early period of inflationary expansion [1], a stage of decelerated matter dominance [2], and accelerated expansion at late times [3,4]. These distinct phases are well supported by observational evidence and are essential components of the standard cosmological model [5]. These dynamics can be achieved by introducing various material contents into Einstein’s theory of general relativity [6,7] or by appropriately modifying the theory itself [8].
Nonlinear electrodynamics (NLED) has been a subject of interest since the early 20th century, primarily due to its potential to address issues not explained by classical Maxwellian electrodynamics. One of the first significant models was introduced by Born and Infeld in 1934; it aimed to eliminate the singularities associated with the self-energy of point charges [9]. NLED theories have since been studied extensively for their ability to describe high-intensity electromagnetic fields and their applications in various areas of physics, including cosmology and astrophysics [10,11,12]. In cosmology, NLED can provide alternative mechanisms for inflation, structure formation, and the late-time accelerated expansion of the universe [13,14]. These models are particularly valuable for their potential to offer insights into the early universe’s dynamics and the nature of dark energy and dark matter.
In the context of cosmological models, NLED has been extensively studied, offering alternative mechanisms for early-universe inflation, structure formation, and late-time accelerated expansion. Various Lagrangians have been proposed to capture the nonlinear effects of electrodynamics in different cosmological scenarios. For instance, the Born–Infeld-type Lagrangian, L = β 2 ( 1 + F 2 β 2 G 2 16 β 4 1 ) , has been analyzed in early-universe models, FRW cosmologies, and Bianchi universes, with a focus on the squared sound speed c s 2 [13,15,16,17,18,19,20,21]. Another model, represented by L = β 2 [ ( 1 + β F σ β γ G 2 2 σ ) σ 1 ] , explores nonsingular universes and bouncing scenarios, demonstrating that the NLED behaves as a cosmological constant at early and late times [22]. Additionally, the Lagrangian L = 1 μ o 4 F + α F 2 + β G 2 has been studied extensively in the context of magnetic universes, inflation, and singularity-free scenarios within FRW cosmologies [23,24,25,26,27]. The Lagrangian L = 1 μ o 4 F + α F 2 has been analyzed primarily for early-time expansion, inflation, and singularity-free models in FRW cosmologies [19,21,23,24,28,29,30]. Other notable Lagrangians include L = F 2 β F + 1 , which has been studied in FRW and Bianchi I universes, showing early-time acceleration without singularities [31,32,33,34], and L = F e α F , which addresses inflation, singularity-free scenarios, and the squared sound speed c s 2 [29,35]. These studies highlight the diverse applications and significant potential of NLED in addressing key questions in cosmology, from early-universe dynamics to the late-time accelerated expansion.
The stability and causality analysis of NLED models ensures their physical viability, especially when considering cosmological applications. One of the critical indicators of stability is the behavior of the squared sound speed (SSS), defined as c s 2 = d p / d ρ [36,37]. A stable cosmological model requires a positive SSS ( c s 2 > 0 ) to prevent the uncontrolled growth of energy density perturbations, which leads to classical instabilities known as gradient instabilities [38]. Additionally, the SSS must not exceed the local speed of light ( c s 2 1 ) to avoid superluminal propagation, which would violate causality [39,40,41]. Early works on NLED, such as [25], laid the foundation for these analyses, and subsequent studies, including [21], have provided detailed discussions on the bounds of the SSS and their implications for the stability and causality of NLED models. Ensuring these conditions are met is important for the theoretical consistency and observational compatibility of NLED-based cosmological models. Similarly, the Lagrangian L = F e α F has been studied for its implications in inflationary scenarios and singularity-free models, with specific focus on ensuring that c s 2 remains within the acceptable range to avoid instabilities [29,35]. Ensuring the stability and causality of NLED models through c s 2 analysis is a theoretical exercise and a necessary step in validating these models against observational data and ensuring their applicability in describing the universe’s evolution.
The present paper explores several cosmological models based on NLED, so the background radiation goes beyond Maxwell’s radiation. Specifically, we investigate NLED models as potential candidates to replicate the observed dynamics of the universe. The novelty of the present investigation is in the combined study we perform through using different complementary tools: (1) the tools of dynamical systems theory, which are useful to expose the asymptotic cosmological dynamics of the models, i.e., to uncover those solutions which are generic and truly relevant; (2) the study of the instabilities which are due to small perturbations of the cosmic background filled with non-Maxwellian radiation, as well as of the tachyonic instabilities which arise due to violation of causality; and (3) a parameter fitting investigation through implementing a Bayesian analysis of the relevant cosmological parameters. The latter study allows us to determine the values of the model parameters that best match the empirical data, providing a robust test of the theoretical models against real-world observations. This procedure enables us to quantify the degree of agreement between the theoretical predictions and the observed data, thus assessing the viability of the models. In addition, by constraining the parameter values, we can reduce the uncertainty in the model predictions, leading to more precise cosmological insights. The dynamical systems tools, in the meantime, allow one to inspect the global (asymptotic) dynamics of the cosmological models in some phase space. The equilibrium asymptotic configurations amount to those generic solutions which are preferred by the system of differential equations on which the cosmological model is based. Attractor solutions represent those cosmological stages which the system naturally evolves to in the asymptotic future. Meanwhile, source critical points represent those asymptotic cosmological stages in the past from where the system evolves into the future. Saddle equilibrium points represent transient stages to which the system is temporarily attracted but from where it is repelled into the future attractor. The importance of the stability study has been already established.
The importance of the present study lies in the power of the combined study we perform here. There are separate parameter fitting studies as well as dynamical systems and stability investigations of several cosmological models which are based on nonlinear electromagnetic fields as the sources of gravitation. According to these studies some models can successfully explain the inflationary stage of early-time cosmic dynamics, while others can explain the present stage of accelerated expansion of the universe without the need for dark energy. As we shall show in this paper for two different NLED-based models, after the application of the combined study we are proposing here, we consistently show that none of these are acceptable models of our cosmos despite published work leading to contrary results. Hence, our approach serves as a consistent way to rule out cosmological models which otherwise are investigated as serious candidates for either the inflaton or dark energy.
We have organized the paper as follows: Section 2 presents the theoretical framework of NLED coupled to general relativity. Section 3 analyzes the stability and causality of various NLED Lagrangians. Section 4 performs a dynamical analysis of models with stable and causal Lagrangian densities. Section 5 provides Bayesian parameter estimation using observational data. Finally, Section 6 summarizes the conclusions and suggests directions for future research.

2. Nonlinear Electrodynamics Coupled to General Relativity

2.1. Modified Friedmann Equations

This paper will consider the Einstein gravitational equations coupled with nonlinear radiation, represented by a nonlinear Lagrangian. In this context, the four-dimensional (4D) action of gravity coupled to nonlinear electrodynamics is given by ( 8 π G c 4 = 1 , geometrical units)
S = d 4 x g R + L m + L ( F , G ) ,
where R is the curvature scalar, L m is the Lagrangian density of the background perfect fluid, and L ( F , G ) is the gauge-invariant electromagnetic (EM) Lagrangian density, which is a function of the electromagnetic invariants
F = F μ ν F μ ν = 2 ( B 2 E 2 ) , G = 1 2 ϵ α β μ ν F α β F μ ν = 4 E · B ,
where B and E are the magnetic induction and electric fields, respectively, and F μ ν = μ A ν ν A μ , with A μ being the electromagnetic potential. In this work, we do not include a cosmological constant Λ in the action, since we want to investigate whether nonlinear electrodynamics can contribute to the current stage of accelerated expansion of the universe without relying on Λ . Our approach allows us to explore the possibility that the effects of nonlinear electrodynamics might provide an alternative explanation for dark energy.
In electrodynamics, the Maxwell–Lagrangian density can be formulated differently, affecting the resulting equations of motion and the interpretation of physical quantities. Specifically, the Lagrangian density for the electromagnetic field can be written as L = F , where F = 1 4 F μ ν F μ ν , or L = F / 4 , where F = F μ ν F μ ν . The first formulation, L = F , includes the normalization factor within the definition of the invariant F, simplifying some mathematical expressions. However, this approach can complicate the interpretation of the energy and momentum densities. On the other hand, the formulation L = 1 4 F explicitly includes the normalization factor in the Lagrangian, which aligns more closely with the standard conventions in classical electrodynamics, making the interpretation of physical quantities more straightforward. Therefore, we adopt this formulation in this work, as it is particularly useful in deriving and analyzing the equations of motion and the stability conditions for the nonlinear electrodynamics (NLED) models considered in this study.
The Born–Infeld theory was the first nonlinear extension of electromagnetism to avoid the singularities associated with infinite electric fields in Maxwell’s theory [9]. The Born–Infeld Lagrangian can be expressed as
L = b 2 1 1 1 2 b 2 F ,
where F = F μ ν F μ ν . Core concepts of NLED [10,42,43], their applications across various physics disciplines [44] and potential avenues for further research are explored. However, this work focuses solely on nonlinear electrodynamics as the material content for a cosmological model, excluding further research directions.
The corresponding gravitational field equations can be derived from the action (1) by performing variations with respect to the spacetime metric g μ ν to obtain
G μ ν = T μ ν m + T μ ν E M ,
where G μ ν is the Einstein tensor containing all the geometric information, and
T μ ν m = ρ m + p m u μ u ν p m g μ ν , T μ ν E M = g μ ν L ( F , G ) G L , G 4 F μ α F ν α L , F ,
are the energy–momentum tensors for ordinary matter (m) and the nonlinear electromagnetic field ( E M ), respectively. Here, ρ m = ρ m ( t ) and p m = p m ( t ) are the energy density and barotropic pressure of the background fluid, respectively, and u μ is the normalized ( u μ u μ = 1 ) velocity of the reference frame where the fields are measured. Also, L , X d L / d X and L , X X d 2 L / d X 2 , etc.
Variation of the action (1) with respect to the components of the electromagnetic potential A μ results in the electromagnetic field equations, referred to as modified Maxwell equations:
μ F μ ν L , F + 1 2 ϵ α β μ ν F α β L , G = 0 ,
where μ denotes the covariant derivative.
Observations indicate that the current universe is very close to a spatially flat geometry [45,46,47]. Therefore, in this paper, we consider a homogeneous and isotropic Friedman–Lemaitre–Robertson–Walker (FLRW) universe with flat spatial sections, described by the Robertson–Walker metric:
d s 2 = d t 2 a ( t ) 2 δ i j d x i d x j ,
where a ( t ) is the cosmological scale factor, and the Latin indices run over three-space.
In order to effectively incorporate NLED into a homogeneous and isotropic geometric framework, we employ an averaging technique that satisfies specific criteria [48]. These criteria include ensuring that the volumetric average of the electromagnetic field remains direction-independent [48], that field fluctuations are equally probable in all directions [49,50], and that there is no net energy flow as observed by comoving observers. Additionally, it is assumed that the electric and magnetic fields, as random fields, possess coherent lengths much shorter than cosmological horizon scales. This ensures that the NLED equations are compatible with the FLRW geometry, facilitating a consistent analysis of cosmological dynamics [14,20,21,23,51].
Under these conditions, the average energy–momentum tensor adopts the perfect fluid form:
T μ ν E M = ( ρ E M + p E M ) u μ u ν c 2 p E M g μ ν ,
where the density and pressure of the nonlinear radiation are given by
ρ E M = L + G L , G 4 L , F E 2 , p E M = L G L , G 4 3 2 B 2 E 2 L , F ,
where E 2 and B 2 are the averaged electric and magnetic fields squared.
In cosmological models coupled with NLED, particularly in the context of the early universe, the analysis often focuses solely on magnetic fields, setting the electric field to zero ( E = 0 ). This simplification stems from several compelling reasons. Firstly, the early universe’s hot, dense plasma environment favors magnetic fields due to magnetohydrodynamic effects and their energy density scaling with B 2 compared to E 2 [50]. Secondly, limiting the analysis to magnetic fields reduces the complexity of the dynamical equations. Thirdly, NLED can give rise to intriguing magnetic phenomena, such as magnetic monopoles or the self-organization of magnetic structures, making magnetic fields a central focus for understanding these phenomena and their potential cosmological implications. Finally, in many cosmological scenarios, the effects of the electric field are negligible compared to those of magnetic fields, particularly in the early universe. This approach, commonly referred to as studying magnetic universes (MUs), enables the investigation of the specific effects of NLED on magnetic phenomena in the early universe while maintaining a manageable analytical framework [14,20,21,23,51]. Thus,
ρ B = L , p B = L 4 3 F L , F ,
and F = 2 B 2 , G = 0 .
In the context of an MU, where NLED is incorporated, the Friedmann equations describe the dynamics of the universe. These equations can be written for the energy density and pressure components, considering the contribution of a magnetic field. The total energy density ρ t and pressure p t in such a universe include contributions from both matter and the nonlinear electromagnetic field.
We are now ready to formulate the dynamical equations for the FLRW model coupled with dark matter and NLED. Even this simplified framework can provide significant physical insights. We shall focus on MUs driven by electromagnetic Lagrangian densities that depend only on the invariant F = 2 B 2 . In this case, the cosmological equations take the following form:
3 H 2 = ρ m L , 2 H ˙ = ( ρ m + p m ) + 4 3 F L , F , ρ ˙ m = 3 H ρ m ( 1 + ω m ) , F ˙ = 4 H F ,
where H = a ˙ / a is the Hubble parameter (the overdot denotes a derivative with respect to cosmic time t), and ω m is the equation of state (EOS) parameter of the ordinary matter.
The solutions to the last equations in (8) are given by
ρ m = ρ 0 m a 3 ( 1 + ω m ) , F = F 0 a 4 ,
where F 0 and ρ 0 m are integration constants, and a = a ( t ) is the scale factor.

2.2. Cosmological Parameters

This section provides an in-depth introduction of cosmological parameters that play a fundamental role in our analysis [52]. These parameters allow us to characterize the different critical points that emerge from the analysis using dynamical systems. These critical points represent the initial conditions from which all trajectories of cosmological model evolution originate, the points toward which these trajectories inevitably converge during their evolution, and the transient stages through which the models pass on their evolution. Understanding these parameters provides insights into various features that can be tested against astronomical observations.

2.2.1. Energy Density Parameter

The energy density parameter, denoted as Ω , is a dimensionless quantity that plays an important role in cosmology. It is defined as the ratio of the energy density ρ x of a particular component x (such as matter, radiation, or dark energy) to the critical density ρ c r i t , which is the density required for the universe to be spatially flat. Mathematically, it is expressed as
Ω x = ρ x ρ c r i t , ρ c r i t = 3 H 2 ,
where H is the Hubble parameter.
The energy density parameter allows cosmologists to characterize the universe’s composition and predict its dynamics. It provides insights into whether the universe is open ( Ω < 1 ), closed ( Ω > 1 ), or flat ( Ω = 1 ).
The energy density parameter is essential for interpreting observational data, evaluating the stability of cosmological models, and understanding the universe’s ultimate fate. It helps identify the relative contributions of different components, such as dark matter, dark energy, and ordinary matter, to the overall energy budget of the universe.

2.2.2. EOS Parameter

The EOS parameter, ω , is a key relationship in cosmology and fluid physics that describes the connection between the pressure p and the energy density ρ of a perfect fluid:
ω = p ρ .
where p is the pressure of the fluid, and ρ is its energy density. This parameter is important because it determines how the fluid affects the universe’s expansion. Depending on the value of ω , the dynamic behavior of the fluid can vary significantly.
When considering a cosmological model with multiple material components, such as matter, radiation, and nonlinear electromagnetic fields, the effective barotropic parameter ω e f f provides a useful way to describe the overall relationship between the total pressure and the total energy density of the universe. The effective barotropic parameter is defined as
ω e f f = p t ρ t ,
where p t and ρ t are the total pressure and energy density, respectively, summed over all components. The effective barotropic parameter ω e f f can also be expressed in terms of the Hubble parameter H as follows:
ω e f f = 1 2 H ˙ 3 H 2 .
This expression provides a direct relationship between the effective barotropic parameter and the Hubble parameter, offering insights into the overall dynamics of the universe.
In a similar way, we can define the effective EOS parameter of the NLED electromagnetic field:
ω e f f B = p B ρ B = 1 + 4 L , F 3 L F .

2.2.3. Deceleration Parameter

The deceleration parameter, denoted by q, is a dimensionless quantity that describes the rate of change of the universe’s expansion rate. It indicates whether the expansion of the universe is accelerating or decelerating. The deceleration parameter is defined as
q = a ¨ a a ˙ 2 ,
where a is the scale factor and a ˙ and a ¨ are the first and the second time derivatives of a.
The deceleration parameter can be expressed in terms of the Hubble parameter H and its time derivative H ˙ as follows:
q = 1 H ˙ H 2 .
It is essential for understanding the universe’s expansion history and future evolution.

2.2.4. Squared Sound Speed

The SSS, denoted as c s 2 , is a fundamental parameter in cosmology and fluid dynamics that describes the propagation speed of pressure (acoustic) waves through a given medium. It is defined as the ratio of the change in pressure p to the change in energy density ρ and is mathematically expressed as [36,37]
c s 2 = d p d ρ .
The SSS is important for several reasons. The value of c s 2 bounds the stability and causality within a given cosmological model. In order to assure stability of the background against small perturbations, c s 2 must be non-negative ( c s 2 0 ) [21,38]. A negative SSS leads to exponential growth of perturbations, indicating a so-called gradient instability. Even if c s 2 is a positive quantity, a causality issue may arise whenever the squared sound speed exceeds the local speed of light1. Indeed, it is commonly assumed that c s 2 1 , while the complementary bound c s 2 > 1 is employed as a criterion for rejecting theories [21,39,40].
In a similar way, we can define the squared sound speed in terms of the NLED electromagnetic field as
c s 2 = d p B d ρ B = 1 3 + 4 3 F L , F F L , F ,
where p B and ρ B are given by Equations (7).

2.3. Statefinder Analysis

In this subsection, we present the statefinder analysis, an essential diagnostic tool in cosmology used to distinguish between different dark energy models and to compare them with the standard Λ CDM model. The statefinder parameters [53] provide a higher-order characterization of the expansion dynamics of the universe beyond the Hubble parameter H and the deceleration parameter q.
The statefinder analysis is based on two dimensionless parameters, r and s, derived from the universe’s scale factor a ( t ) . These parameters are defined as follows:
r = a a H 3 , s = r 1 3 ( q 1 / 2 ) ,
where, as before, H is the Hubble parameter and q is the deceleration parameter. The jerk parameter j is closely related to the statefinder parameter r and is defined as
j = 1 H 3 d 3 a d t 3 .
In a flat universe, the jerk parameter and the statefinder parameter r are identical, i.e., j = r [54]. The statefinder parameters r and s (19) can also be expressed in terms of the Hubble parameter and its derivatives as follows:
r = H ¨ H 3 + 3 H ˙ H 2 + 1 = 1 2 F L , F 3 H 2 1 + 4 L , F F L , F F = 1 2 c s 2 F L , F H 2 ,
while
s = 4 F L , F 1 + 4 L , F F L , F F 3 ( 4 F L , F 3 L ) = c s 2 4 F L , F 4 F L , F 3 L ,
where we have used (18).
The statefinder parameters, r and s, serve several critical purposes in cosmology. Providing a more nuanced view of the universe’s dynamics, they help differentiate between dark energy models, such as quintessence, Chaplygin gas, and braneworld models. The Λ CDM model, considered the standard model of cosmology, corresponds to the fixed point ( r , s ) = ( 1 , 0 ) . Any deviation from this point indicates a difference from the Λ CDM model. Additionally, while the Hubble parameter H and the deceleration parameter q offer information about the first and second derivatives of the scale factor, r and s incorporate the third derivative, enhancing the characterization of the universe’s evolution. In order to perform the statefinder analysis for a given cosmological model, one must first determine the form of the scale factor a ( t ) for the model, then derive H ( t ) and q ( t ) from a ( t ) , and finally, use the definitions of r and s to calculate these parameters.
For specific models, the Λ CDM model has statefinder parameters ( r , s ) = ( 1 , 0 ) , quintessence models generally exhibit r < 1 and s > 0 , and Chaplygin gas models typically fall in the region where r > 1 and s < 0 . As we shall see, for the NLED-based models under study the values of the statefinder parameters can even change sign during cosmic expansion. In particular, the sign of the parameter s is correlated with the sign of c s 2 , as seen from (21): for a definite sign of the ratio 4 F L , F / ( 4 F L , F 3 L ) , whenever the dynamics become unstable against small perturbations of the cosmic background, s changes sign. However, even for models that are stable against small perturbations, the parameter s changes sign whenever the above ratio changes sign (see Section 6).

3. Stability and Causality Analysis of Several NLED Lagrangians

We shall consider Lagrangians that only depend on the invariant F where E = 0 , so that only magnetic fields are assumed. The adiabatic speed of sound squared (SSS) for scalar perturbations is given by (18). Hence, the bounds on the SSS, 0 c s 2 1 , translate into the following bounds on the NLED Lagrangian and its F-derivatives [21]:
1 4 F L , F F L , F 1 2 .
These bounds ensure the given model is free of gradient and tachyon instabilities. Otherwise, following the causality principle, the group velocity of excitations over the background should be less than the speed of light, ensuring no tachyons are in the theory’s spectrum. The unitarity principle guarantees the absence of ghosts. Both principles lead to the inequalities [55]
L , F 0 , L , F F 0 , L , F + 2 F L , F F 0 ,
for a Lagrangian that only depends on the invariant F. As we can see, the stability and causality of the material content of the model corresponding to NLED radiation depend strongly on its Lagrangian dependence on the electromagnetic invariant F.

3.1. Power-Law NLED Lagrangian

Magnetic universes have been extensively explored within the framework of NLED theories, often characterized by simple Lagrangian densities. One of the simplest Lagrangians is given by
L = 1 4 F + α F 2 ,
where the nonlinear term F 2 [19,21,23,24,28,29,30] has been suggested to potentially induce a cosmic bounce, thus avoiding the initial singularity known as the Big Bang [23]. On the other hand, Lagrangians featuring inverse powers of the electromagnetic field F are intriguing due to the potential significance of nonlinear electromagnetic effects in both the early and late stages of cosmic evolution. Models employing Lagrangian densities like [14,21,51]
L = 1 4 F γ F ,
have been proposed to elucidate the late-time accelerated expansion of MUs [14]. Additionally, combinations of positive and negative powers of F have been investigated [21,51,56], as exemplified by
L = 1 4 F γ F + α F 2 .
This composite model captures key cosmic evolutionary stages and remarkably avoids the cosmological Big Bang singularity. Specifically, the quadratic term F 2 dominates during early epochs, facilitating a nonsingular bounce, while the Maxwell term F takes precedence in the radiation era. At late times, the term F 1 governs, driving cosmic acceleration [51].
In [21], it was demonstrated that numerous cosmological models rooted in NLED Lagrangians (24)–(26) are susceptible to curvature singularities of sudden and/or big rip varieties, or exhibit pronounced instability against minor perturbations of the cosmological background—often attributed to the negative sign of SSS. Moreover, concerns regarding causality may emerge due to the potential superluminal propagation of background perturbations.
The subsequent extension of the Lagrangian density involves incorporating a power-law dependence on the invariant F, expressed as
L = γ F α .
In this scenario, upon evaluating the SSS Equation (18), we find c s 2 = 1 + 4 3 α . It becomes evident that the Lagrangian density of NLED remains causal and stable within the range of α values spanning 3 / 4 α 3 / 2 .
Now, let us examine the Lagrangian density of NLED expressed as
L = 1 4 F γ F α .
This model characterizes a universe featuring linear Maxwell radiation alongside a nonlinear power-law term. A prior investigation [21] explored specific cases of α ( α = 1 and α = 2 ), revealing that these models fail to meet the necessary constraints for SSS. In [57], the authors examined a variant incorporating cold dark matter ( ρ m with p m = 0 ). Their analysis, anchored in astronomical observations and with γ < 0 , demonstrated that the Lagrangian (28) accurately replicates dark energy phenomena for α = 1 / 4 and α = 1 / 8 . More recently, Ref. [58] explored a cosmological framework involving a variable Newton’s constant, G ( t ) , alongside an NLED of the form L = F α 4 . They established its stability within the range 7 / 4 α 5 / 2 .
The expression for the SSS Equation (18), corresponding to the Lagrangian (28) in terms of the scale factor a and the parameter α , reads
c s 2 = 1 3 + 16 α ( α 1 ) γ F α 12 α γ F α + 3 F = 1 3 + 16 α ( α 1 ) γ F 0 α 1 a 4 ( 1 α ) 12 α γ F 0 α 1 a 4 ( 1 α ) + 3 ,
where on the right-hand side of this equation we have used the solution to the equation of the conservation of the electromagnetic field (9): F = F 0 a 4 .
The plot of the SSS (29) in the right panel of Figure 1 shows that the model (28) remains stable and causal within a narrow range of α values, specifically 3 / 4 α 3 / 2 . In [57], the NLED Lagrangian (28) was considered, and for the obtained values of α , namely, α = 1 / 4 and α = 1 / 8 ), the resulting SSSs are c s 2 = 4 / 3 and c s 2 = 7 / 6 , respectively, indicating gradient instability. Similarly, the cases with α = 1 and α = 2 , as shown in [21], yield c s 2 = 7 / 3 and c s 2 = 5 / 3 , respectively, reaffirming the gradient and also tachyon instabilities in these models.
We have already commented on the fact that for NLED-based models like the one given by the Lagrangian (28) a change in sign of the statefinder parameter s is correlated with a change in sign of the squared speed of sound c s 2 . In the present case,
s = c s 2 4 F L , F 4 F L , F 3 L = 4 c s 2 1 + 4 α γ F 0 α 1 a 4 ( 1 α ) 1 4 ( 3 4 α ) γ F 0 α 1 a 4 ( 1 α ) ,
where c s 2 is given by (29). It can be seen that for α > 3 / 4 , since c s 2 0 during the course of the cosmic evolution, and the ratio R 4 F L , F / ( 4 F L , F 3 L ) is a positive-valued function as well, the statefinder parameter s does not change sign. In general, for non-negative α -s, this parameter does not change sign. However, for negative α < 0 , we have the following picture. At very early times ( a 0 ), c s 2 = 1 / 3 > 0 and the ratio R = 4 > 0 , while at late times ( a ), c s 2 = 4 α / 3 1 < 0 and R = 4 α / ( 4 α 3 ) > 0 . This means that for α < 0 , the statefinder parameter s = c s 2 R changes sign from s > 0 to s < 0 during cosmic expansion. When s > 0 , the expansion is more similar to models like a quintessence, where the present accelerated expansion of the universe is explained as being due to a form of dark energy different from a cosmological constant. However, when s < 0 , the expansion suggests models that could be more similar to those with phantom energy or other exotic forms of energy.

3.2. Generalized Rational Nonlinear Electrodynamics

In this section, we explore a generalized model of NLED described by the following Lagrangian density:
L = F 4 b F 1 + ϵ ( 2 ϵ β F ) α ,
where b, ϵ = ± 1 and α are dimensionless parameters, and β is a parameter with dimensions [ L ] 4 . To ensure that the Lagrangian remains real, in [59] is stated that ϵ = 1 for B > E and ϵ = 1 for B < E . However, in our case, since E = 0 , we only consider the scenario where ϵ = 1 . We recover Maxwell’s electrodynamics when b = 0 . Furthermore, when we apply the weak field approximation ( β F 1 ), the power series expansion of L (Equation (31)) reveals that nonlinearity persists. This implies the existence of nonlinearity even at late times, which is important for exploring its implications for the model’s behavior during those periods, providing this Lagrangian density proves to be stable and causal.
Some particular cases of the Lagrangian density (31) have been previously analyzed in various studies [31,60,61,62,63,64], primarily investigating black hole solutions with different values of α , particularly when α = 1 / 2 , and for both values of ϵ . In [59], it was demonstrated that singularities of point electric charges are absent, and the electromagnetic energy is finite.
The expression for the SSS corresponding to the Lagrangian (31) in terms of the scale factor a and the parameters α , β , and b, is as follows:
c s 2 = 1 3 16 3 b α ( α 1 ) ( 2 β F ) 2 α ( α + 1 ) ( 2 β F ) α ( 1 + ( 2 β F ) α ) ( 2 4 b ( α 1 ) ) ( 2 β F ) α + ( 2 β F ) 2 α + 4 b + 1 = 1 3 + 16 b α v 1 α + ( 1 + α ) v 3 ( 1 + v ) 1 + 2 k v + ( 1 + 4 b ) v 2 ,
where we have defined the constant k 1 + 2 b ( 1 α ) and we have used the compact notation v ( a 4 / 2 β F 0 ) α . If we assume that the constants β and b are non-negative reals, a pure qualitative analysis complemented with the numerical study (see Figure 2) reveals that there is always a non-empty range of α values where the model is free of instability against small perturbations of the background and also of tachyon instability, i.e., where SSS meets the bounds 0 c s 2 1 , at all stages of the cosmic expansion. For instance, if we take b = 1 , for 2 / 3 α 1 , the model is both causal and stable. The width of the α —range depends on the constant b: For larger values of the constant b the range is wider.

Generalized Rational Nonlinear Electrodynamics without Maxwell Term

Let us examine in more detail a generalized model of rational nonlinear electrodynamics [59], excluding the linear Maxwell term and focusing on magnetic universes, described by the following Lagrangian density ( ϵ = 1 ):
L = b F 1 + ( 2 β F ) α ,
For the limit β 0 , the Lagrangian (33) reduces to Maxwell’s theory. Specific cases of this Lagrangian density have been studied in the cosmological context [59,60], detailed in [31,32,33,34,65]. A very special case of the Lagrangian density (33) where α = 1 has been studied in [31,32,33,60,66].
The squared sound speed for the Lagrangian density (33) is given by
c s 2 = 1 3 4 3 α [ ( α 1 ) ( 2 β F ) 2 α ( 2 β F ) α ( 1 + α ) ] [ 1 + ( 2 β F ) α ] [ ( α 1 ) ( 2 β F ) α 1 ] = 1 3 4 α 1 α + ( 1 + α ) v 3 ( 1 + v ) ( 1 α + v ) ,
where, as before, v ( a 4 / 2 β F 0 ) α . A straightforward inspection of (34) shows that, for 1 / 2 α 1 / 4 , the model (33) is stable against small perturbations of the background and is free of tachyonic instability as well (see Figure 3). As seen in Figure, in the region over the surface c s 2 = c s 2 ( v , α ) (34), lying between the vertical planes at α = 1 / 2 and α = 1 / 4 , the squared speed of sound meets the bounds 0 c s 2 1 , within which the model is both causal and stable.
The equations of motion of the model (33) are complex enough to warrant further investigation through a dynamical systems analysis, particularly within the parameter range 1 / 2 α 1 / 4 , where no issues with instability against small perturbations of the background or with tachyon instability arise.
Recently, in [65], the dynamics of this model have been analyzed. However, the stability and causality analysis in that work were verified only for values which fall outside of the range 1 / 2 α 1 / 4 . Therefore, in our work, we can perform a more comprehensive dynamical analysis within the parameter space we have identified as stable and causal. From the above, we can infer a specific region in the parameter space for the Lagrangian (33), where a dynamical analysis can be conducted in a homogeneous and isotropic framework to observe its effects during early and late times.

4. Dynamical Systems Analysis of Models with Stable and Causal Lagrangian Density

In this section, we perform a dynamical analysis of models where nonlinear radiation is stable and causal. Given the complexity of the equations of motion, here we utilize the tools provided by dynamical systems theory for this. These tools allow one to inspect the global (asymptotic) dynamics of cosmological models in some auxiliary phase space where the equilibrium asymptotic configurations amount to those solutions which are preferred by the system of differential equations of the cosmological model under study [67]. Attractor solutions represent those cosmological stages which the system naturally evolves to in the asymptotic future. Meanwhile, source critical points represent those asymptotic cosmological stages in the past from where the system evolves into the future. In particular, if past attractors exist, these mark the points in the phase space in whose neighborhood the initial conditions are to be given. Transient stages of the cosmic evolution such as, for instance, primordial inflation and matter dominance, are associated with saddle equilibrium points. Our task is to determine whether there are past and future attractors, as well as saddle critical points, in the phase space of the cosmological models under study.

4.1. Power-Law Lagrangian

In the previous section, we found that a power-law model of nonlinear electrodynamics (28) is stable and causal as long as 3 / 4 α 3 / 2 . We now analyze the evolution of a homogeneous and isotropic model whose matter content is given by this NLED and ordinary matter (cold dark matter). Thus, if we take the density for ordinary matter as ρ m (considering non-relativistic dust ( ω m = 0 ) and the nonlinear electromagnetic field given by the Lagrangian (28), the dynamical equations to study are the following (8):
3 H 2 = ρ m + F 4 + γ F α , 2 H ˙ = ρ m + 1 3 F + 4 3 γ α F α , ρ ˙ m = 3 H ρ m ( 1 + ω m ) , F ˙ = 4 H F .
Dynamical systems theory requires the transformation of this second-order system into a first-order system by going into an auxiliary phase space spanned by the variables [67]
x = F 12 H 2 = Ω rad , y = γ F α 3 H 2 = Ω rad nled ,
Thus, the first of Equations (35) can be written as Ω m + Ω rad + Ω rad nled = 1 , or equivalently,
Ω m = 1 x y ,
where we see that the variables are bounded within the phase space
Ψ = { ( x , y ) | 0 x + y 1 , 0 x 1 , 0 y 1 } .
The second Equation (35) takes the form
2 H ˙ H 2 = 3 Ω m + 4 x + 4 α y = 3 + x + ( 4 α 3 ) y .
The dynamical system reads
x = 2 x 2 + H ˙ H 2 = x [ x 1 + ( 4 α 3 ) y ] , y = 2 y 2 α + H ˙ H 2 = y [ ( 4 α 3 ) y + x + 3 4 α ] ,
where the prime denotes the derivative with respect to τ , which is the conformal time, τ = d t a ( t ) . In terms of the variables (36) of the phase space, the squared sound speed (29) can be written as
c s 2 = x + α ( 4 α 3 ) y 3 ( x + α y ) .
Finally, both the effective barotropic parameter, ω e f f , in Equation (13) and the deceleration parameter q in Equation (16) in the new variables are given by
ω e f f = 1 3 x + 4 α 3 3 y , q = 1 2 1 + x + ( 4 α 3 ) y
and the statefinder parameters are given by
r = 1 + 2 x + 2 α ( 4 α 3 ) y , s = 4 [ x + α ( 4 α 3 ) y ] 3 [ 3 x + ( 4 α 1 ) y ] .
The critical points of the dynamical system (40) within the phase space Ψ (38) are the following:
  • At the first critical point, P m ( 0 , 0 ) , ordinary matter dominates ( Ω m = 1 ). The behavior of this critical point, as determined by the eigenvalues λ 1 = 1 and λ 2 = 3 4 α , indicates that it is a decelerated ( q = 1 / 2 ) future attractor when α > 3 / 4 ( c s 2 -stable) and a saddle point when α < 3 / 4 ( c s 2 -unstable). The effective barotropic parameter ω eff = 0 indicates that the combined effect of ordinary matter and nonlinear radiation is such that the expansion of the universe behaves as if it is dominated by non-relativistic matter. This is consistent with a dust-dominated universe within the framework of standard cosmology. Finally, the statefinder parameters take the values r = 1 and s = undefined. This indicates that the cosmological model exhibits characteristics of the Λ CDM model with respect to the parameter r. However, the parameter s is undefined, so that there must be appreciable deviations from the Λ CDM model. This reflects the influence of the nonlinear electrodynamics on the universe’s dynamics, distinguishing it from a purely matter-dominated model.
  • The second point P rad ( 1 , 0 ) , where Maxwell radiation dominates with Ω rad = 1 and c s 2 = 1 3 , has eigenvalues λ 1 = 1 and λ 2 = 4 4 α . This indicates that it is a decelerated ( q = 1 ) source (past attractor) when α < 1 and a saddle point when α > 1 . An effective EOS parameter ω e f f = 1 / 3 indicates that the combined effect of ordinary matter and nonlinear radiation results in the expansion of the universe behaving as if it were dominated by radiation. The statefinder parameters take the values r = 3 and s = 4 / 9 . These values indicate significant deviations from the Λ CDM model. These deviations reflect the nonlinear electrodynamics’ impact on the dynamics of the universe.
  • For the point P rad nled ( 0 , 1 ) , where nonlinear radiation dominates Ω rad nled = 1 , the squared sound speed is given by c s 2 = ( 4 α 3 ) / 3 . The eigenvalues λ 1 = 4 α 3 and λ 2 = 4 ( α 1 ) indicate that this point is a future attractor if α < 3 / 4 ( c s 2 -unstable), a saddle point if 3 / 4 < α < 1 ( c s 2 -stable), and a source (past attractor) if α 1 ( c s 2 -stable if α 3 / 2 ). The deceleration parameter at this point is q = 2 α 1 , which indicates an accelerated expansion if α < 1 2 ; however, this value of the α -parameter falls outside of the range where the model is stable and causal. For α > 1 2 the expansion is always decelerated. The effective EOS parameter ω e f f = 1 + 4 α / 3 . The statefinder parameters for this critical point take the values r = 1 + 2 α ( 4 α 3 ) and s = 4 α ( 4 α 3 ) / 3 ( 4 α 1 ) . The dependence of the statefinder parameters on α reflects the impact of nonlinear electrodynamics on the dynamics of the universe, providing insights into appreciable deviations from the Λ CDM model.
The evolution of the model can be divided into three scenarios depending on the value of the parameter α :
  • For α < 3 / 4 the dynamics of the model start from a Maxwell radiation domain ( P rad ), transiently approach the matter-dominated phase ( P m ), and are finally attracted towards the nonlinear radiation domain ( P rad nled ). The transient matter-dominated stage could explain the correct amount of structure formation. In the case when α < 1 / 2 , the attractor solution depicts accelerated expansion, so that the nonlinear radiation could explain the accelerated expansion of the universe without the need for any dark energy. Unfortunately, this scenario is not physically acceptable because the values of the parameter α fall outside of the range 3 / 4 α 3 / 2 , where the model is free of gradient and tachyon instabilities.
  • For 3 / 4 α < 1 , all trajectories start from the Maxwell radiation-dominated point ( P rad ), then these approach to a transient stage where the nonlinear radiation dominates ( P rad nled ), and eventually, reach the matter-dominated attractor ( P m ). The values of the statefinder parameters vary from one critical point to another. This means that the model exhibits a dynamic behavior quite distinct from the standard cosmological model, indicating that the influence of nonlinear electrodynamics introduces substantial differences in the evolution of the universe compared to a Λ CDM framework, see the left panel of the Figure 4.
  • For 1 < α < 3 / 2 , the orbits in the phase space are sourced by the nonlinear radiation ( P rad nled ), evolve transiently near the Maxwell radiation-dominated point ( P rad ), and are finally attracted towards the matter-dominated point ( P m ). The values of the statefinder parameters change from critical point to critical point, for instance, at the past attractor
    r = 1 + 2 α ( 4 α 3 ) , s = 4 α ( 4 α 3 ) 3 ( 4 α 1 ) ,
    while at the future attractor r = 1 , s = undefined, and at the Maxwell point ( P rad ) r = 3 and s = 4 / 9 . These values do not closely align with the consensus Λ CDM model, signaling appreciable departures from it, see the right panel of the Figure 4.
In this section, we have analyzed the dynamic evolution of a nonlinear power-law model of electrodynamics characterized by the Lagrangian density (28). Previously, we have demonstrated that this model is stable and causal within the parameter range 3 / 4 α 3 / 2 . The dynamics of this model exhibit two relevant scenarios: In the first scenario, the model starts in a Maxwell radiation-dominated state, transitions through a transient phase dominated by nonlinear radiation, and ultimately culminates in a matter-dominated state. In the second scenario, the model begins in a nonlinear radiation-dominated state, transiently passes through a Maxwell radiation-dominated phase, and finally, reaches a matter-dominated state. These results provide a comprehensive understanding of the dynamic evolution and stability of the nonlinear power-law electrodynamics model. Furthermore, the statefinder analysis reveals that while the model’s parameters do not align closely with the Λ CDM model, they still offer valuable insights into the deviations introduced by nonlinear electrodynamics. Specifically, the statefinder parameters help quantify how much the model diverges from the standard cosmological model, highlighting the unique NLED effects on the evolution of the universe.

4.2. Generalized Rational Nonlinear Electrodynamics Dynamical Evolution without Maxwell Term

The next model we will analyze using dynamical systems corresponds to the Lagrangian density given by (33), which we have demonstrated to be causal and stable for the range 1 / 2 α 1 / 4 . The dynamical equations for this model are given by
3 H 2 = ρ m + b F 1 + ( 2 β F ) α , 2 H ˙ = ρ m 4 b F 3 ( α 1 ) ( 2 β F ) α 1 ( 1 + ( 2 β F ) α ) 2 ,
and F ˙ = 4 H F . Here, we consider ordinary (pressureless) matter ρ m ( p m = 0 ).
We consider the following new bounded phase plane variables to obtain an autonomous first-order dynamical system:
x = b F b F + 3 H 2 , y = ( 2 β F ) α 1 + ( 2 β F ) α .
Thus, Equations (46) in these new variables can be expressed as
Ω m = 1 Ω rad nled = 1 x ( 1 y ) 1 x , 2 H ˙ H 2 = 3 Ω m + 4 x ( 1 y ) ( 1 α y ) 1 x ,
where
Ω m = ρ m 3 H 2 , Ω rad nled = b F / 3 H 2 1 + ( 2 β F ) α .
The deceleration parameter, effective equation of state (EOS) parameter, and the SSS (34) are given by
q = 1 x + x ( 1 y ) ( 1 4 α y ) 2 ( 1 x ) , ω e f f = x ( 1 y ) ( 1 4 α y ) 3 ( 1 x ) , c s 2 = 1 3 4 α y ( 1 + α 2 α y ) 3 ( 1 α y ) ,
respectively, while the statefinder parameters read
r = 1 + 2 x ( 1 y ) [ 1 ( 5 + 4 α ) α y + 8 α 2 y 2 ] 1 x , s = 4 [ 1 ( 5 + 4 α ) α y + 8 α 2 y 2 ] 3 ( 1 4 α y ) .
The 2D dynamical system for this model reads
x = 2 x ( 1 x ) 2 + H ˙ H 2 , y = 4 α y ( 1 y ) ,
where the prime denotes differentiation with respect to τ = ln a . The phase space in which to look for equilibrium configurations of the dynamical system (47) is defined by
Ψ = ( x , y ) 0 x 1 , 0 y 1 , y 2 x 1 x ,
where the lower bound on the variable y comes from requiring that Ω m 0 .
Below, we show the critical points of the dynamical system (47). We include the values of several cosmological parameters of observational significance evaluated at the critical points:
  • The first critical point, P m ( 0 , 0 ) , dominated by ordinary matter ( Ω m = 1 ), is decelerated ( q = 1 2 ), and its effective barotropic parameter simulates dust ( ω e f f = 0 ). At this point, the SSS corresponds to Maxwell radiation, meaning c s 2 = 1 3 . The Jacobian matrix’s eigenvalues are λ 1 = 1 and λ 2 = 4 α . This indicates that P 1 is a future attractor if α > 0 and a saddle point if α < 0 . The values of the statefinder parameters at this point are both positive: r = 1 and s = 4 / 3 . The significant deviation of s from zero indicates a considerable divergence from the standard cosmological model. This substantial positive value of s reflects a significant influence of nonlinear electrodynamics on the universe’s dynamics.
  • The point P M ( 0 , 1 ) corresponds to decelerated expansion as well ( q = 1 / 2 ), and it is also dominated by ordinary matter. The effective EOS parameter ω e f f = 0 means that the background fluid mimics dust. The SSS at this point is given by c s 2 = ( 1 4 α ) / 3 . Therefore, the stability and causality depend on the value of α , which must satisfy the condition found in the previous section: 1 / 2 α 1 / 4 . The eigenvalues of the Jacobian matrix are λ 1 = 1 and λ 2 = 4 α , indicating that P 2 is a future attractor if α < 0 and a saddle point otherwise. The values of the statefinder parameters for this point are given by r = 1 and s = 4 ( 1 α ) / 3 . The parameter s varies with α , indicating how the nonlinear electrodynamics affects the model’s deviation from the standard cosmological model. Specifically, the expression for s highlights the sensitivity of the model’s behavior to the parameter α . Only when α = 1 does the model exhibit behavior similar to Λ CDM, where s = 0 .
  • The point P rad 1 / 2 , 0 is dominated by Maxwell radiation ( ω e f f = 1 / 3 ), with
    Ω rad = b F 12 H 2 = x 1 x = 1 ,
    and represents decelerated expansion ( q = 1 ). The SSS c s 2 = 1 / 3 . The eigenvalues of the Jacobian matrix at this point are λ 1 = 1 and λ 2 = 4 α , indicating that this point acts as a past attractor if α < 0 and as a saddle point for other values. For the statefinder parameters, we obtain that r = 3 and s = 4 / 3 . These values imply that the cosmological model at this critical point behaves very differently from the Λ CDM model.
  • For the critical point P rad nled = ( 1 , 1 ) , several of the relevant cosmological parameters such as Ω m , q , ω e f f , and the statefinder parameter r, are undefined. The only parameters which one can determine are c s 2 = ( 1 4 α ) / 3 and s = 4 ( 1 α ) / 3 . The eigenvalues of the Jacobian matrix at this point are λ 1 = 1 and λ 2 = 4 α , indicating that it acts as a past attractor if α > 0 , otherwise it is a saddle point. These features of the NLED-based equilibrium configuration suggests that the nonlinear effects of the magnetic field play an important role in the course of cosmic expansion.
The features of the cosmological model (33) are illustrated in Figure 5 for two distinct values of the parameter α which fall within the range 1 / 2 α 1 / 4 , where the model is free of gradient and tachyon instabilities: α = 0.1 and α = 0.1 , respectively. Regardless of the specific value of α , the phase space consistently displays a past attractor (blue point), a saddle point (black point), and a future attractor (red point). This structure indicates that the initial conditions should be set near those of the past attractor, from which the orbits of the dynamical system move towards the saddle point—representing a transient phase in cosmic evolution—before converging on the future attractor. The primary effect of the parameter α is to determine the roles of the different asymptotic states, interchanging their positions as past or future attractors and the saddle point configuration.

4.3. Generalized Rational Nonlinear Electrodynamics Dynamical Evolution with Maxwell Term

In this section, we analyze the dynamical evolution of a model corresponding to the Lagrangian density given by (31) with ϵ = 1 . Finding exact solutions of the equations of motion of this model is a complex task, and even if we find several exact cosmological solutions, these may belong in a set of measure zero. This is why in this and in the other cases in this paper we apply the tools of dynamical systems theory. Our task is to search for the critical points of the resulting dynamical system in some auxiliary phase space. These equilibrium configurations amount to those solutions of the original system of cosmological equations which are generic, and so, are the most relevant ones within the infinite set of possible exact solutions.
The dynamical equations for the model (31) are given by
3 H 2 = ρ m + F 4 + b F 1 + ( 2 β F ) α , 2 H ˙ = ρ m 4 F 3 1 4 b ( 1 + ( 2 β F ) α ) + b α ( 2 β F ) α ( 1 + ( 2 β F ) α ) 2 ,
where we are considering ordinary matter ρ m (non-relativistic dust ω m = 0 ) and F ˙ = 4 H F .
We consider the following new bounded variables:
x = F F + 12 H 2 , y = ( 2 β F ) α 1 + ( 2 β F ) α .
Thus, the dynamical Equation (49) in these new variables are traded by the following plane-autonomous system of ordinary differential equations:
x = 2 x ( 1 x ) 2 + H ˙ H 2 , y = 4 α y ( 1 y ) ,
where
Ω m = ρ m 3 H 2 = 1 x [ 1 + 4 b ( 1 y ) ] ( 1 x ) , H ˙ H 2 = 3 2 Ω m 2 x ( 1 x ) 1 + 4 b ( 1 y ) ( 1 α y ) .
The phase space where the physically meaningful critical points of the dynamical system (51) are located is defined by
Ψ = ( x , y ) 0 x 1 , 0 y 1 , 2 x ( 1 + 2 b ) 1 4 b x y 1 + 4 b 4 b .
There are four critical points associated with this dynamical system:
  • Matter-dominated critical point P m ( 0 , 0 ) ( Ω m = 1 ). At this point the electromagnetic effects are negligible, making this point representative of a purely matter-dominated phase. The deceleration parameter is q = 1 / 2 , and the effective EOS parameter ω e f f = 0 . The squared sound speed is c s 2 = 1 / 3 , while the statefinder parameters are r = 1 and s = ( 1 + 4 b ) / 3 ( 1 + b ) . The eigenvalues associated with the linearization matrix at this critical point are λ 1 = 1 and λ 2 = 4 α , indicating that this point is a decelerated future attractor if α > 0 and a saddle point otherwise.
  • Second matter-dominated critical point P M ( 0 , 1 ) ( Ω m = 1 ). This represents a matter-dominated phase with significant NLED effects, F . The deceleration parameter is q = 1 / 2 , the effective EOS parameter ω e f f = 0 and the squared sound speed c s 2 = 1 / 3 . The statefinder parameters are r = 1 and s = 1 / 3 . The eigenvalues associated with this critical point are λ 1 = 1 and λ 2 = 4 α , indicating that this point is a saddle point if α > 0 and decelerated future attractor if α < 0 .
  • The critical point P rad = ( 1 / 2 , 1 ) , with
    Ω rad = F 12 H 2 = x 1 x = 1 , Ω rad nled = b F / 3 H 2 1 + ( 2 β F ) α = 4 b x ( 1 y ) 1 x = 0 ,
    corresponds to Maxwell radiation domination (the effects of NLED radiation are negligible). This point is characterized by the following values of the relevant cosmological parameters: q = 1 , ω e f f = 1 / 3 , and the squared sound speed c s 2 = 1 3 . The values of the statefinder parameters are r = 3 and s = 1 / 3 . The eigenvalues of the linearization matrix evaluated at this point are λ 1 = 1 and λ 2 = 4 α , respectively, indicating that this point is a non-accelerated expanding past attractor if α > 0 and saddle point if not.
  • The critical point P rad nled = 1 2 ( 1 + 2 b ) , 0 is dominated by the combined effect of the nonlinear magnetic and Maxwell radiation fields,
    Ω rad = 1 1 + 4 b , Ω rad nled = 4 b 1 + 4 b Ω rad + Ω rad nled = 1 .
    The point exists for b 1 / 4 . The deceleration parameter is q = 1 , the effective EOS parameter ω e f f = 1 / 3 , and the statefinder parameters r = 4 and s = ( 1 + 4 b ) / 3 ( 1 + b ) . The eigenvalues of the linearization matrix are λ 1 = 1 and λ 2 = 4 α , respectively, indicating a non-accelerated expanding saddle point if α > 0 and a past attractor if α < 0 . When b = 0 , this point reduces to the point P = ( 1 / 2 , 0 ) , and when b , it reduces to the matter-dominant point P m = ( 0 , 0 ) . In the former case, there is not a nonlinear magnetic field, while in the latter case the nonlinear effects decouple from the spectrum of general relativity plus a Maxwell radiation field.
The evolution of the model is graphically depicted in Figure 6 for two different values of the parameter α . The left panel illustrates the case when α = 0.5 , while the right panel corresponds to α = 0.5 .

5. Parameter Fitting with Observational Data

Performing parameter fitting using observational data is essential to further validate and refine our models. Parameter fitting allows us to determine the values of the model parameters that best match the empirical data, providing a robust test of the theoretical models against real-world observations. This process is essential for several reasons. Firstly, it enables us to quantify the degree of agreement between the theoretical predictions and the observed data, thus assessing the viability of the models. Secondly, by constraining the parameter values we can reduce the uncertainty in the model predictions, leading to more precise and reliable cosmological insights.
In this study, we use 31 cosmic chronometers collected over several years [68,69,70,71,72,73,74,75]. Additionally, we use high-precision baryon acoustic oscillation (BAO) measurements at different redshifts up to z < 2.36 from BOSS DR14 quasars (eBOSS) [76], SDSS DR12 Galaxy Consensus [77], Ly- α DR14 cross-correlation [78], Ly- α DR14 auto-correlation [79], the Six-Degree Field Galaxy Survey (6dFGS) [80], and SDSS Main Galaxy Sample (MGS) [81]. In addition, we include SNeIa data from the Pantheon compilation [82] with 1048 supernovae. We also incorporate a compressed version of Planck-15 information, treating the CMB as a BAO experiment at redshift z = 1090 , measuring the angular scale of the sound horizon [83]. These datasets provide a comprehensive and high-precision set of observations that span a significant range of the universe’s history, making them ideal for testing cosmological models incorporating nonlinear electrodynamics (NLED). By fitting the parameters of our NLED models to these observations, we aim to derive constraints that either support or challenge the theoretical framework, ultimately contributing to a deeper understanding of the universe’s evolution.
We employ Bayesian inference to determine the values of the model parameters that best fit the observational data. In this study, we utilize a nested sampling algorithm, using 500 live points, within the library dynesty [84], and the SimpleMC code [85]. The resulting posterior distributions provide constraints on the parameters, enabling us to assess the viability and precision of the NLED models in describing the universe’s evolution. The results of the Bayesian inference are shown in Table 1 and in the plots of Figure 7. In Table 1, as a reference, we include the parameter estimation of the Λ CDM model using the same data and parameters for nested sampling. This allows us, using the logarithm of the Bayesian evidence log Z , to calculate the Bayes factor log B , and through the Bayes factor and Jeffrey’s scale [86,87], to make a model comparison.
Therefore, we employ Bayesian inference to investigate the parameters of the power-law and rational Lagrangian models, aiming to constrain them using observational data from the late universe. Specifically, we aim to assess whether the nonlinear electrodynamics models are compatible with the available observational evidence and whether the results of our dynamical analysis align with the estimated parameters.
By fitting the parameters of the power-law and rational Lagrangian models to the comprehensive observational datasets, we aim to derive constraints on the model parameters. This approach allows us to evaluate the compatibility of the NLED models with current observational evidence and to determine whether these models can provide a viable description of the universe’s evolution. In addition, we will quantify the degree of agreement between the theoretical predictions of the NLED models and the empirical data. The results of this analysis provide critical insights into the viability and precision of the NLED models in the context of modern cosmology.
In the case of the power-law Lagrangian, the Friedmann Equation (28), useful for performing Bayesian inference, is given by
H 2 H 0 2 = Ω m 0 ( 1 + z ) 3 ( 1 + ω m ) + Ω r 0 ( 1 + z ) 4 + Ω N L 0 ( 1 + z ) 4 α ,
where Ω m 0 = ρ m 0 3 H 0 2 , Ω r 0 = F 0 12 H 0 2 , and Ω N L 0 = γ F 0 α 3 H 0 2 . We use the following flat priors: Ω m [ 0.1 , 0.5 ] , h [ 0.4 , 0.9 ] , and α [ 2 , 2 ] .
For the rational Lagrangian,
L = F 4 b F 1 + 2 β F ,
which is a particular case of (33) for α = 1 ; the equation to be used is
H 2 H 0 2 = Ω m 0 ( 1 + z ) 3 + Ω r 0 ( 1 + z ) 4 + 4 b Ω r 0 ( 1 + z ) 4 24 β H 0 2 Ω r 0 ( 1 + z ) 4 + 1 ,
where Ω m 0 = ρ m 0 3 H 0 2 and Ω r 0 = F 0 12 H 0 2 . We assume β = 1 and use the flat priors Ω m [ 0.1 , 0.5 ] , h [ 0.4 , 0.9 ] , and b [ 27 / 125 , 108 / 343 ] for the Bayesian parameter estimation.
We observe that the datasets used cannot constrain the parameters α and b effectively; however, the values of the matter density Ω m and the Hubble parameter h for both models are consistent with the expected values for the current universe. The fact that the α and b parameters can take any value within the prior ranges aligns with the results found in the dynamical analysis.
In the triangle plots of Figure 7, the posterior sampling for both Lagrangians and their respective ω ( z ) values, using Equation (11) for both cases, can be seen. According to the Bayes factor in Table 1, the Λ CDM model is significantly preferred over the two NLED models when using late-universe observations. This is expected, as neither of the NLED models exhibits a cosmological constant behavior in the late universe, as shown in the EoS plots in Figure 7.
The Bayesian analysis results shown in Table 1 indicate that the Λ CDM model is favored significantly over the power-law and rational NLED models. The very strong evidence in favor of Λ CDM, indicated by the Bayes factor log B , aligns with the expectations based on cosmological constant behavior in the late universe, which neither NLED models can replicate. This suggests that while NLED models provide interesting theoretical insights, they may not be suitable for describing the current accelerated expansion of the universe.
Moreover, the inability to tightly constrain the parameters α and b within the prior ranges indicates that further theoretical refinement and possibly new observational data are needed to fully assess the viability of NLED models. The consistency of Ω m and h with expected values reinforces the reliability of the data and the robustness of the Λ CDM model as the standard cosmological model.

6. Discussion and Conclusions

This study explored the dynamics of cosmological models incorporating nonlinear electrodynamics (NLED), focusing on their stability and causality. Employing a combination of dynamical systems theory and Bayesian inference, we analyzed two specific NLED models: the power-law [57] and rational Lagrangian models [31,65].
We identified stable and causal parameter ranges for the power-law Lagrangian model. In the dynamical systems analysis, this model transitions through various cosmological phases, starting from a Maxwell radiation-dominated state and evolving to a matter-dominated state. For the rational Lagrangian model, including the Maxwell term, we observed stable and causal behavior for specific ranges of the parameter b. The phase space analysis revealed critical points indicating the evolutionary pathways of the universe, beginning from an early radiation-dominated state to an ordinary-matter-dominated state.
We performed Bayesian parameter estimation using a comprehensive set of observational data, including cosmic chronometers, baryon acoustic oscillation (BAO) measurements, and supernovae type Ia (SNeIa). The estimated parameters for both models were consistent with the expected values for the current universe, particularly the matter density Ω m and the Hubble parameter h. However, the parameters α and b could not be tightly constrained within the prior ranges, which aligns with the findings from the dynamical analysis.
Based on Bayesian evidence, our model comparison strongly favored the Λ CDM model over the NLED models for late-universe observations. This result was anticipated, as neither NLED model exhibited an accelerated stage of cosmic expansion for values of the parameters falling within the region in the parameter space where the models are stable and causal. Hence, based on our combined study, we may conclude that none of the analyzed models can satisfactorily explain the primordial inflation or the current stage of accelerated universe expansion. This is despite there being works in the bibliography where the models studied in the present paper are thoroughly investigated either as dark energy or as the inflaton field. For instance, in [57] the power-law model (28) is tested with observational data of type ia supernovae, long gamma-ray bursts, and Hubble parameter measurements. According to the authors, the statistical analysis shows that the inclusion of nonlinear electromagnetic matter is enough to produce the observed accelerated expansion and that there is no need to include a dark energy component. The α -parameter values favored by the observations are 1 , 1 / 4 and 1 / 8 . These values are far from the range 3 / 4 α 3 / 2 , where the model is free of the gradient and tachyon instabilities. As shown in our dynamical systems study of this model for 3 / 4 α 3 / 2 , there are not accelerated solutions in the phase portrait corresponding to this model. Another example is given by the work in [65], where the author analyses the inflation of the universe when the source of gravity is an electromagnetic field obeying nonlinear electrodynamics, with two parameters and without singularities, which is given by the Lagrangian (31). Even if there can be inflationary solutions, these amount to an empty set compared with the infinity of possible solutions of the equations of motion. As we have shown, the auxiliary phase space of this model does not have critical points (these correspond to the generic solutions which are ‘preferred’ by the model) which can be associated with inflation within the parameter range where the model is stable against small perturbations of the cosmic background and is causal.
We have shown, based on our combined study, that to account for early-time inflation and for the present stage of accelerated expansion, the models (28), (31), and (33) would require the inclusion of an inflaton field and a dark energy component, respectively. Consequently, these models fail to comprehensively explain inflation and the dark energy component of the cosmic background. This limits their interest as standalone models for describing the whole evolution of the universe.
This investigation comprehensively explains the stability, causality, and dynamical evolution of cosmological models driven by nonlinear electrodynamics. While the NLED models exhibit intriguing theoretical properties, their compatibility with observational data suggests that further refinement and exploration are needed to fully integrate them into the standard cosmological framework. Future work should focus on extending the analysis to other forms of NLED and exploring their implications for early-universe phenomena and high-energy astrophysical processes.

Author Contributions

Conceptualization, R.G.-S. and I.Q.; methodology, R.G.-S., V.M.-B. and I.Q.; software, I.G.-V.; validation, T.G., V.M.-B. and I.Q.; formal analysis, R.G.-S. and I.G.-V.; investigation, R.G.-S. and I.Q.; writing—original draft preparation, R.G.-S., I.G.-V. and V.M.-B.; writing—review and editing, R.G.-S., I.G.-V., V.M.-B. and I.Q.; visualization, R.G.-S., I.G.-V., T.G., V.M.-B. and I.Q.; funding acquisition, R.G.-S. and I.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Consejo Nacional de Ciencia y Tecnología: FORDECYT-PRONACES-CONACYT CF-MG-2558591. Secretaría de Investigación y Posgrado, Instituto Politécnico Nacional: SIP20230505-IPN and SIP20240638.

Data Availability Statement

This study generated no new data. All data used in this research have been previously published and are properly referenced in the manuscript. Further details can be found in the cited references.

Acknowledgments

We thank the anonymous reviewers for their valuable comments and suggestions, which have greatly improved our manuscript’s quality. R.G-S. thanks the institutional support by COFAA and EDI of the Instituto Politécnico Nacional, México.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the study’s design, in the collection, analyses, or interpretation of data, in the writing of the manuscript, or in the decision to publish the results.

Note

1
For a critical review on this issue, see [37].

References

  1. Guth, A.H. Inflationary universe: A possible solution to the horizon and flatness problems. Phys. Rev. 1981, 23, 347. [Google Scholar] [CrossRef]
  2. Peebles, P.J.E. Principles of Physical Cosmology; Princeton University Press: Princeton, NJ, USA, 1993; Volume 27. [Google Scholar]
  3. Riess, A.G.; Filippenko, A.V.; Challis, P.; Clocchiatti, A.; Diercks, A.; Garnavich, P.M.; Gilliland, R.L.; Hogan, C.J.; Jha, S.; Kirshner, R.P.; et al. Observational evidence from supernovae for an accelerating universe and a cosmological constant. Astron. J. 1998, 116, 1009. [Google Scholar] [CrossRef]
  4. Perlmutter, S.; Aldering, G.; Goldhaber, G.; Knop, R.; Nugent, P.; Castro, P.G.; Deustua, S.; Fabbro, S.; Goobar, A.; Groom, D.E.; et al. Measurements of Ω and Λ from 42 high-redshift supernovae. Astrophys. J. 1999, 517, 565. [Google Scholar] [CrossRef]
  5. Weinberg, S. Cosmology; OUP Oxford: Oxford, UK, 2008. [Google Scholar]
  6. Grøn, Ø.; Hervik, S. Einstein’s General Theory of Relativity: With Modern Applications in Cosmology; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2007. [Google Scholar]
  7. Ryder, L. Introduction to General Relativity; Cambridge University Press: Cambridge, UK, 2009. [Google Scholar]
  8. Clifton, T.; Ferreira, P.G.; Padilla, A.; Skordis, C. Modified Gravity and Cosmology. Phys. Rep. 2012, 513, 1–189. [Google Scholar]
  9. Born, M.; Infeld, L. Foundations of the new field theory. Proc. R. Soc. London. Ser. Contain. Pap. Math. Phys. Character 1934, 144, 425–451. [Google Scholar] [CrossRef]
  10. Plebański, J. Lectures on Non-Linear Electrodynamics: An Extended Version of Lectures Given at the Niels Bohr Institute and NORDITA, Copenhagen, in October 1968; Nordita: Stockholm, Sweden, 1970; Volume 25. [Google Scholar]
  11. Boillat, G. Nonlinear electrodynamics: Lagrangians and equations of motion. J. Math. Phys. 1970, 11, 941–951. [Google Scholar] [CrossRef]
  12. Novello, M.; Bergliaffa, S.P.; Salim, J. Singularities in general relativity coupled to nonlinear electrodynamics. Class. Quantum Gravity 2000, 17, 3821. [Google Scholar] [CrossRef]
  13. Garcia-Salcedo, R.; Breton, N. Born–Infeld cosmologies. Int. J. Mod. Phys. 2000, 15, 4341–4353. [Google Scholar] [CrossRef]
  14. Novello, M.; Bergliaffa, S.P.; Salim, J. Nonlinear electrodynamics and the acceleration of the universe. Phys. Rev. 2004, 69, 127301. [Google Scholar] [CrossRef]
  15. Dyadichev, V.; Gal’tsov, D.; Zorin, a.G.; Zotov, M.Y. Non-Abelian Born-Infeld cosmology. Phys. Rev. 2002, 65, 084007. [Google Scholar] [CrossRef]
  16. Moniz, P.V. Quintessence and Born-Infeld cosmology. Phys. Rev. 2002, 66, 103501. [Google Scholar] [CrossRef]
  17. Elizalde, E.; Lidsey, J.E.; Nojiri, S.; Odintsov, S.D. Born–Infeld quantum condensate as dark energy in the universe. Phys. Lett. 2003, 574, 1–7. [Google Scholar] [CrossRef]
  18. Dyadichev, V.V.; Gal’tsov, D.V.; Moniz, P.V. Chaos-order transition in Bianchi type I non-Abelian Born-Infeld cosmology. Phys. Rev. 2005, 72, 084021. [Google Scholar] [CrossRef]
  19. de Oliveira Costa, E.G.; Bergliaffa, S.E.P. A classification of the effective metric in nonlinear electrodynamics. Class. Quantum Gravity 2009, 26, 135015. [Google Scholar] [CrossRef]
  20. García-Salcedo, R.; Gonzalez, T.; Horta-Rangel, A.; Quiros, I. Comment on “Extended Born-Infeld theory and the bouncing magnetic universe”. Phys. Rev. 2014, 90, 128301. [Google Scholar] [CrossRef]
  21. Garcia-Salcedo, R.; Gonzalez, T.; Quiros, I. No compelling cosmological models come out of magnetic universes which are based on nonlinear electrodynamics. Phys. Rev. 2014, 89, 084047. [Google Scholar] [CrossRef]
  22. Kruglov, S. Notes on Born–Infeld-type electrodynamics. Mod. Phys. Lett. 2017, 32, 1750201. [Google Scholar] [CrossRef]
  23. De Lorenci, V.; Klippert, R.; Novello, M.; Salim, J. Nonlinear electrodynamics and FRW cosmology. Phys. Rev. 2002, 65, 063501. [Google Scholar] [CrossRef]
  24. Camara, C.; de Garcia Maia, M.; Carvalho, J.; Lima, J.A.S. Nonsingular FRW cosmology and nonlinear electrodynamics. Phys. Rev. 2004, 69, 123504. [Google Scholar] [CrossRef]
  25. Novello, M. Cosmological effects of nonlinear electrodynamics. Int. J. Mod. Phys. 2005, 20, 2421–2430. [Google Scholar] [CrossRef]
  26. Kruglov, S. Vacuum birefringence from the effective Lagrangian of the electromagnetic field. Phys. Rev. D—Part. Fields Gravit. Cosmol. 2007, 75, 117301. [Google Scholar] [CrossRef]
  27. Sharif, M.; Mumtaz, S. Stability of the accelerated expansion in nonlinear electrodynamics. Eur. Phys. J. 2017, 77, 1–8. [Google Scholar] [CrossRef]
  28. Camara, C.; Carvalho, J.; De Garcia Maia, M. Nonlinearity of electrodynamics as a source of matter creation in a flat FRW cosmology. Int. J. Mod. Phys. 2007, 16, 427–432. [Google Scholar] [CrossRef]
  29. Campanelli, L.; Cea, P.; Fogli, G.; Tedesco, L. Inflation-produced magnetic fields in nonlinear electrodynamics. Phys. Rev. 2008, 77, 043001. [Google Scholar] [CrossRef]
  30. Vollick, D.N. Homogeneous and isotropic cosmologies with nonlinear electromagnetic radiation. Phys. Rev. 2008, 78, 063524. [Google Scholar] [CrossRef]
  31. Kruglov, S. Universe acceleration and nonlinear electrodynamics. Phys. Rev. 2015, 92, 123523. [Google Scholar] [CrossRef]
  32. Otalora, G.; Övgün, A.; Saavedra, J.; Videla, N. Inflation from a nonlinear magnetic monopole field nonminimally coupled to curvature. J. Cosmol. Astropart. Phys. 2018, 2018, 3. [Google Scholar] [CrossRef]
  33. Singh, G.; Hulke, N.; Singh, A. Accelerating cosmologies with nonlinear electrodynamics. Can. J. Phys. 2018, 96, 992–998. [Google Scholar] [CrossRef]
  34. Kruglov, S. Inflation of universe due to nonlinear electrodynamics. Int. J. Mod. Phys. 2017, 32, 1750071. [Google Scholar] [CrossRef]
  35. Kruglov, S. Nonlinear electromagnetic fields as a source of universe acceleration. Int. J. Mod. Phys. 2016, 31, 1650058. [Google Scholar] [CrossRef]
  36. Peebles, P.J.E.; Ratra, B. The cosmological constant and dark energy. Rev. Mod. Phys. 2003, 75, 559. [Google Scholar] [CrossRef]
  37. Ellis, G.F.; Maartens, R.; MacCallum, M.A. Causality and the speed of sound. Gen. Relativ. Gravit. 2007, 39, 1651–1660. [Google Scholar] [CrossRef]
  38. Yang, R.; Qi, J. Dynamics of generalized tachyon field. Eur. Phys. J. 2012, 72, 1–7. [Google Scholar] [CrossRef]
  39. Hawking, S.W.; Ellis, G.F.R. Large Scale Struct.-Space-Time; Cambridge University Press: Cambridge, UK, 1973; Volume 1. [Google Scholar]
  40. Wald, R.M. General Relativity; University of Chicago Press: Chicago, IL, USA, 2010. [Google Scholar]
  41. Adams, A.; Arkani-Hamed, N.; Dubovsky, S.; Nicolis, A.; Rattazzi, R. Causality, analyticity and an IR obstruction to UV completion. J. High Energy Phys. 2006, 2006, 014. [Google Scholar] [CrossRef]
  42. Gaete, P.; Helayël-Neto, J. Remarks on nonlinear electrodynamics. Eur. Phys. J. 2014, 74, 1–8. [Google Scholar] [CrossRef]
  43. Gaete, P.; Helayël-Neto, J.A. A note on nonlinear electrodynamics. Europhys. Lett. 2017, 119, 51001. [Google Scholar] [CrossRef]
  44. Sorokin, D.P. Introductory Notes on Non-linear Electrodynamics and its Applications. Fortschritte Phys. 2022, 70, 2200092. [Google Scholar] [CrossRef]
  45. Spergel, D.N.; Bean, R.; Doré, O.; Nolta, M.; Bennett, C.; Dunkley, J.; Hinshaw, G.; Jarosik, N.E.; Komatsu, E.; Page, L.; et al. Three-year Wilkinson Microwave Anisotropy Probe (WMAP) observations: Implications for cosmology. Astrophys. J. Suppl. Ser. 2007, 170, 377. [Google Scholar] [CrossRef]
  46. Aghanim, N.; Akrami, Y.; Ashdown, M.; Aumont, J.; Baccigalupi, C.; Ballardini, M.; Banday, A.J.; Barreiro, R.; Bartolo, N.; Basak, S.; et al. Planck 2018 results-VI. Cosmological parameters. Astron. Astrophys. 2020, 641, A6. [Google Scholar]
  47. Alam, S.; Aubert, M.; Avila, S.; Balland, C.; Bautista, J.E.; Bershady, M.A.; Bizyaev, D.; Blanton, M.R.; Bolton, A.S.; Bovy, J.; et al. Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Cosmological implications from two decades of spectroscopic surveys at the Apache Point Observatory. Phys. Rev. 2021, 103, 083533. [Google Scholar] [CrossRef]
  48. Tolman, R.C.; Ehrenfest, P. Temperature equilibrium in a static gravitational field. Phys. Rev. 1930, 36, 1791. [Google Scholar] [CrossRef]
  49. Lemoine, D. Fluctuations in the relativistic plasma and primordial magnetic fields. Phys. Rev. 1995, 51, 2677. [Google Scholar] [CrossRef]
  50. Lemoine, D.; Lemoine, M. Primordial magnetic fields in string cosmology. Phys. Rev. 1995, 52, 1955. [Google Scholar] [CrossRef]
  51. Novello, M.; Goulart, E.; Salim, J.; Bergliaffa, S.P. Cosmological effects of nonlinear electrodynamics. Class. Quantum Gravity 2007, 24, 3021. [Google Scholar] [CrossRef]
  52. Dodelson, S.; Schmidt, F. Mod. Cosmol.; Academic Press: Cambridge, MA, USA, 2020. [Google Scholar]
  53. Sahni, V.; Saini, T.D.; Starobinsky, A.A.; Alam, U. Statefinder—A new geometrical diagnostic of dark energy. J. Exp. Theor. Phys. Lett. 2003, 77, 201–206. [Google Scholar] [CrossRef]
  54. Chaudhary, H.; Arora, D.; Debnath, U.; Mustafa, G.; Maurya, S.K. Addressing the rd Tension using Late-Time Observational Measurements in a Novel Deceleration Parametrization. arXiv 2023, arXiv:2308.07354v2. [Google Scholar]
  55. Shabad, A.E.; Usov, V.V. Effective Lagrangian in nonlinear electrodynamics and its properties of causality and unitarity. Phys. Rev. 2011, 83, 105006. [Google Scholar] [CrossRef]
  56. Maity, S.; Chakraborty, S.; Debnath, U. Correspondence Between Electromagnetic Field and Other Dark Energies in Nonlinear Electrodynamics. Int. J. Mod. Phys. 2011, 20, 2337–2350. [Google Scholar] [CrossRef]
  57. Montiel, A.; Bretón, N.; Salzano, V. Parameter estimation of a nonlinear magnetic universe from observations. Gen. Relativ. Gravit. 2014, 46, 1758. [Google Scholar] [CrossRef]
  58. Joseph, G.W.; Övgün, A. Cosmology with variable G and nonlinear electrodynamics. Indian J. Phys. 2021, 96, 1861–1866. [Google Scholar] [CrossRef]
  59. Kruglov, S.I. Nonlinear electrodynamics without singularities. Phys. Lett. A 2024, 493, 129248. [Google Scholar] [CrossRef]
  60. Kruglov, S. A model of nonlinear electrodynamics. Ann. Phys. 2015, 353, 299–306. [Google Scholar] [CrossRef]
  61. Kruglov, S. Nonlinear electrodynamics and magnetic black holes. Ann. Der Phys. 2017, 529, 1700073. [Google Scholar] [CrossRef]
  62. Mazharimousavi, S.H.; Halilsoy, M. Electric Black Holes in a Model of Nonlinear Electrodynamics. Ann. Phys. 2019, 531, 1900236. [Google Scholar] [CrossRef]
  63. Mazharimousavi, S.H.; Halilsoy, M. Note on regular magnetic black hole. Phys. Lett. 2019, 796, 123–125. [Google Scholar] [CrossRef]
  64. Kruglov, S. Nonlinearly charged AdS black holes, extended phase space thermodynamics and Joule–Thomson expansion. Ann. Phys. 2022, 441, 168894. [Google Scholar] [CrossRef]
  65. Kruglov, S. Universe inflation and nonlinear electrodynamics. Eur. Phys. J. 2024, 84, 205. [Google Scholar] [CrossRef]
  66. Kruglov, S. Rational nonlinear electrodynamics causes the inflation of the universe. Int. J. Mod. Phys. 2020, 35, 2050168. [Google Scholar] [CrossRef]
  67. García-Salcedo, R.; Gonzalez, T.; Horta-Rangel, F.A.; Quiros, I.; Sanchez-Guzmán, D. Introduction to the application of dynamical systems theory in the study of the dynamics of cosmological models of dark energy. Eur. J. Phys. 2015, 36, 025008. [Google Scholar] [CrossRef]
  68. Jimenez, R.; Verde, L.; Treu, T.; Stern, D. Constraints on the equation of state of dark energy and the Hubble constant from stellar ages and the cosmic microwave background. Astrophys. J. 2003, 593, 622–629. [Google Scholar] [CrossRef]
  69. Simon, J.; Verde, L.; Jimenez, R. Constraints on the redshift dependence of the dark energy potential. Phys. Rev. 2005, 71, 123001. [Google Scholar] [CrossRef]
  70. Stern, D.; Jimenez, R.; Verde, L.; Kamionkowski, M.; Stanford, S.A. Cosmic chronometers: Constraining the equation of state of dark energy. I: H(z) measurements. J. Cosmol. Astropart. Phys. 2010, 2010, 008. [Google Scholar] [CrossRef]
  71. Moresco, M.; Verde, L.; Pozzetti, L.; Jimenez, R.; Cimatti, A. New constraints on cosmological parameters and neutrino properties using the expansion rate of the Universe to z ∼ 1.75. J. Cosmol. Astropart. Phys. 2012, 2012, 53. [Google Scholar] [CrossRef]
  72. Zhang, C.; Zhang, H.; Yuan, S.; Liu, S.; Zhang, T.J.; Sun, Y.C. Four new observational H(z) data from luminous red galaxies in the Sloan Digital Sky Survey data release seven. Res. Astron. Astrophys. 2014, 14, 1221. [Google Scholar] [CrossRef]
  73. Moresco, M. Raising the bar: New constraints on the Hubble parameter with cosmic chronometers at z∼2. Mon. Not. R. Astron. Soc. Lett. 2015, 450, L16–L20. [Google Scholar] [CrossRef]
  74. Moresco, M.; Pozzetti, L.; Cimatti, A.; Jimenez, R.; Maraston, C.; Verde, L.; Thomas, D.; Citro, A.; Tojeiro, R.; Wilkinson, D. A 6% measurement of the Hubble parameter at z ∼ 0.45: Direct evidence of the epoch of cosmic re-acceleration. J. Cosmol. Astropart. Phys. 2016, 2016, 014. [Google Scholar] [CrossRef]
  75. Ratsimbazafy, A.; Loubser, S.; Crawford, S.; Cress, C.; Bassett, B.; Nichol, R.; Väisänen, P. Age-dating luminous red galaxies observed with the Southern African Large Telescope. Mon. Not. R. Astron. Soc. 2017, 467, 3239–3254. [Google Scholar] [CrossRef]
  76. Ata, M.; Baumgarten, F.; Bautista, J.; Beutler, F.; Bizyaev, D.; Blanton, M.R.; Blazek, J.A.; Bolton, A.S.; Brinkmann, J.; Brownstein, J.R.; et al. The clustering of the SDSS-IV extended Baryon Oscillation Spectroscopic Survey DR14 quasar sample: First measurement of baryon acoustic oscillations between redshift 0.8 and 2.2. Mon. Not. R. Astron. Soc. 2018, 473, 4773–4794. [Google Scholar] [CrossRef]
  77. Alam, S.; Ata, M.; Bailey, S.; Beutler, F.; Bizyaev, D.; Blazek, J.A.; Bolton, A.S.; Brownstein, J.R.; Burden, A.; Chuang, C.H.; et al. The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: Cosmological analysis of the DR12 galaxy sample. Mon. Not. R. Astron. Soc. 2017, 470, 2617–2652. [Google Scholar] [CrossRef]
  78. Blomqvist, M.; Des Bourboux, H.D.M.; de Sainte Agathe, V.; Rich, J.; Balland, C.; Bautista, J.E.; Dawson, K.; Font-Ribera, A.; Guy, J.; Le Goff, J.M.; et al. Baryon acoustic oscillations from the cross-correlation of Lyα absorption and quasars in eBOSS DR14. Astron. Astrophys. 2019, 629, A86. [Google Scholar] [CrossRef]
  79. de Sainte Agathe, V.; Balland, C.; Des Bourboux, H.D.M.; Blomqvist, M.; Guy, J.; Rich, J.; Font-Ribera, A.; Pieri, M.M.; Bautista, J.E.; Dawson, K.; et al. Baryon acoustic oscillations at z = 2.34 from the correlations of Lyα absorption in eBOSS DR14. Astron. Astrophys. 2019, 629, A85. [Google Scholar] [CrossRef]
  80. Beutler, F.; Blake, C.; Colless, M.; Jones, D.H.; Staveley-Smith, L.; Campbell, L.; Parker, Q.; Saunders, W.; Watson, F. The 6dF Galaxy Survey: Baryon acoustic oscillations and the local Hubble constant. Mon. Not. R. Astron. Soc. 2011, 416, 3017–3032. [Google Scholar] [CrossRef]
  81. Ross, A.J.; Samushia, L.; Howlett, C.; Percival, W.J.; Burden, A.; Manera, M. The clustering of the SDSS DR7 main Galaxy sample–I. A 4 per cent distance measure at z = 0.15. Mon. Not. R. Astron. Soc. 2015, 449, 835–847. [Google Scholar] [CrossRef]
  82. Scolnic, D.M.; Jones, D.; Rest, A.; Pan, Y.; Chornock, R.; Foley, R.; Huber, M.; Kessler, R.; Narayan, G.; Riess, A.; et al. The complete light-curve sample of spectroscopically confirmed SNe Ia from Pan-STARRS1 and cosmological constraints from the combined Pantheon sample. Astrophys. J. 2018, 859, 101. [Google Scholar] [CrossRef]
  83. Aubourg, É.; Bailey, S.; Bautista, J.E.; Beutler, F.; Bhardwaj, V.; Bizyaev, D.; Blanton, M.; Blomqvist, M.; Bolton, A.S.; Bovy, J.; et al. Cosmological implications of baryon acoustic oscillation measurements. Phys. Rev. 2015, 92, 123516. [Google Scholar] [CrossRef]
  84. Speagle, J.S. dynesty: A dynamic nested sampling package for estimating Bayesian posteriors and evidences. Mon. Not. R. Astron. Soc. 2020, 493, 3132–3158. [Google Scholar] [CrossRef]
  85. Vazquez, J.; Gomez-Vargas, I.; Slosar, A. Updated Version of a Simple MCMC Code for Cosmological Parameter Estimation Where Only Expansion History Matters. 2020. Available online: https://github.com/ja-vazquez/SimpleMC (accessed on 1 May 2024).
  86. Alberto Vazquez, J.; Bridges, M.; Hobson, M.; Lasenby, A. Reconstruction of the Dark Energy equation of state. J. Cosmol. Astropart. Phys. 2012, 9, 20. [Google Scholar] [CrossRef]
  87. Vázquez, J.A.; Tamayo, D.; Sen, A.A.; Quiros, I. Bayesian model selection on Scalar ϵ-Field Dark Energy. Phys. Rev. 2021, 103, 043506. [Google Scholar]
Figure 1. Plot of the surface c s 2 = c s 2 ( u , α ) given by Equation (29), where u γ ( a 4 / 2 β F 0 ) 1 α , is shown in the left panel, while the projection onto the plane ( u , α ) is shown in the right panel. The parallel vertical planes (red crosshatched diagonal line pattern) that cut the surface are located at the values α = 3 / 4 and α = 3 / 2 , respectively. In the part of the surface which is cut by these planes (red lines in the right panel) the SSS satisfies the bounds 0 c s 2 1 for any u (dark región between the red lines in the right panel).
Figure 1. Plot of the surface c s 2 = c s 2 ( u , α ) given by Equation (29), where u γ ( a 4 / 2 β F 0 ) 1 α , is shown in the left panel, while the projection onto the plane ( u , α ) is shown in the right panel. The parallel vertical planes (red crosshatched diagonal line pattern) that cut the surface are located at the values α = 3 / 4 and α = 3 / 2 , respectively. In the part of the surface which is cut by these planes (red lines in the right panel) the SSS satisfies the bounds 0 c s 2 1 for any u (dark región between the red lines in the right panel).
Universe 10 00353 g001
Figure 2. Plots of the squared speed of sound c s 2 in Equation (32) vs. v a 4 α / ( 2 β F 0 ) α , for different values of the parameter b, are drawn in the top figures, while in the bottom figures the corresponding projections onto the plane ( v , α ) are shown. From left to right, b = 5 , 1, and 1 / 2 , respectively. The region over the surface between the perpendicular planes (strip between red lines in the projection) is where the SSS meets the bounds 0 c s 2 1 . i.e., in these regions is where the model is free of instability against small perturbations of the background and of tachyon instability as well, for any value of v. As seen, for larger values of the constant b the strip is wider in the α parameter: for b = 5 5 / 2 α 5 / 4 , while for b = 1 2 / 3 α 1 and for b = 0.5 2 / 5 α 4 / 5 .
Figure 2. Plots of the squared speed of sound c s 2 in Equation (32) vs. v a 4 α / ( 2 β F 0 ) α , for different values of the parameter b, are drawn in the top figures, while in the bottom figures the corresponding projections onto the plane ( v , α ) are shown. From left to right, b = 5 , 1, and 1 / 2 , respectively. The region over the surface between the perpendicular planes (strip between red lines in the projection) is where the SSS meets the bounds 0 c s 2 1 . i.e., in these regions is where the model is free of instability against small perturbations of the background and of tachyon instability as well, for any value of v. As seen, for larger values of the constant b the strip is wider in the α parameter: for b = 5 5 / 2 α 5 / 4 , while for b = 1 2 / 3 α 1 and for b = 0.5 2 / 5 α 4 / 5 .
Universe 10 00353 g002
Figure 3. Plot of the surface c s 2 = c s 2 ( v , α ) given by Equation (34), where v ( a 4 / 2 β F 0 ) α , is shown in the left panel, while the projection onto the plane ( v , α ) is shown in the right panel. The parallel vertical planes (red crosshatched diagonal line pattern) that cut the surface are located at the values α = 1 / 4 and α = 1 / 2 , respectively. In the part of the surface which is cut by these planes, the SSS satisfies the bounds 0 c s 2 1 for any v.
Figure 3. Plot of the surface c s 2 = c s 2 ( v , α ) given by Equation (34), where v ( a 4 / 2 β F 0 ) α , is shown in the left panel, while the projection onto the plane ( v , α ) is shown in the right panel. The parallel vertical planes (red crosshatched diagonal line pattern) that cut the surface are located at the values α = 1 / 4 and α = 1 / 2 , respectively. In the part of the surface which is cut by these planes, the SSS satisfies the bounds 0 c s 2 1 for any v.
Universe 10 00353 g003
Figure 4. Phase portrait of the dynamical system (40) for power-law NLED. In the left panel, α = 0.85 , while in the right panel, α = 1.25 . These values of the parameter α fall within the interval 3 / 4 α < 3 / 2 , where the model is stable and causal, i.e., where the stability and causality conditions 0 c s 2 1 , are satisfied. The blue points mark the past attractors, while the red points label the future attractors. Saddle points are marked by the black points. Finally, the blue background region corresponds to the inequality where 0 Ω m 1 .
Figure 4. Phase portrait of the dynamical system (40) for power-law NLED. In the left panel, α = 0.85 , while in the right panel, α = 1.25 . These values of the parameter α fall within the interval 3 / 4 α < 3 / 2 , where the model is stable and causal, i.e., where the stability and causality conditions 0 c s 2 1 , are satisfied. The blue points mark the past attractors, while the red points label the future attractors. Saddle points are marked by the black points. Finally, the blue background region corresponds to the inequality where 0 Ω m 1 .
Universe 10 00353 g004
Figure 5. Phase portrait of the dynamical system (47) for different values of the parameter α . In the left panel, the orbits of the dynamical system with α = 0.1 are drawn. In the right panel, these are shown for α = 0.1 . These values of α fall within the interval 1 / 2 α 1 / 4 , where the model (33) is stable and causal. The blue background region corresponds to the inequality where 0 Ω m 1 .
Figure 5. Phase portrait of the dynamical system (47) for different values of the parameter α . In the left panel, the orbits of the dynamical system with α = 0.1 are drawn. In the right panel, these are shown for α = 0.1 . These values of α fall within the interval 1 / 2 α 1 / 4 , where the model (33) is stable and causal. The blue background region corresponds to the inequality where 0 Ω m 1 .
Universe 10 00353 g005
Figure 6. Phase portrait of the dynamical system (51) for b = 1 and different parameter values of the parameter α : α = 0.5 (left panel) and α = 0.5 (right panel). These values fall within the interval 2 / 3 α 1 , where the model is stable and causal for b = 1 . Notice that attractors change place with the saddle points. The blue background region corresponds to the inequality where 0 Ω m 1 .
Figure 6. Phase portrait of the dynamical system (51) for b = 1 and different parameter values of the parameter α : α = 0.5 (left panel) and α = 0.5 (right panel). These values fall within the interval 2 / 3 α 1 , where the model is stable and causal for b = 1 . Notice that attractors change place with the saddle points. The blue background region corresponds to the inequality where 0 Ω m 1 .
Universe 10 00353 g006
Figure 7. Posterior sampling for the parameters using Pantheon, BAO, cosmic chronometers, and the Planck point for the models given in Equations (53) and (55). The left panel shows the parameter estimation for the power-law Lagrangian, while the right panel depicts the rational Lagrangian model. The contour plots illustrate the constraints on Ω m , h, α , b, and Ω r , with the 1 σ and 2 σ confidence regions shaded. The insets display the equation of state parameter ω ( z ) as a function of redshift z for both models, indicating their deviation from the Λ CDM model at late times.
Figure 7. Posterior sampling for the parameters using Pantheon, BAO, cosmic chronometers, and the Planck point for the models given in Equations (53) and (55). The left panel shows the parameter estimation for the power-law Lagrangian, while the right panel depicts the rational Lagrangian model. The contour plots illustrate the constraints on Ω m , h, α , b, and Ω r , with the 1 σ and 2 σ confidence regions shaded. The insets display the equation of state parameter ω ( z ) as a function of redshift z for both models, indicating their deviation from the Λ CDM model at late times.
Universe 10 00353 g007
Table 1. Parameter estimation for the NLED power-law (Equation (53)) and NLED rational (Equation (55)) Lagrangians. The table presents the estimated values of each model’s matter density parameter Ω m , the Hubble parameter h, and the radiation density parameter Ω r . The parameter α is also estimated for the power-law model, while the parameter b is included for the rational model. The log Z values indicate the Bayesian evidence for each model, and the log B values provide the Bayes factor comparison with the Λ CDM model, showing a strong preference for the Λ CDM model. The last row qualitatively describes the strength of evidence favoring Λ CDM over the NLED models.
Table 1. Parameter estimation for the NLED power-law (Equation (53)) and NLED rational (Equation (55)) Lagrangians. The table presents the estimated values of each model’s matter density parameter Ω m , the Hubble parameter h, and the radiation density parameter Ω r . The parameter α is also estimated for the power-law model, while the parameter b is included for the rational model. The log Z values indicate the Bayesian evidence for each model, and the log B values provide the Bayes factor comparison with the Λ CDM model, showing a strong preference for the Λ CDM model. The last row qualitatively describes the strength of evidence favoring Λ CDM over the NLED models.
Power LawRational Λ CDM
Ω m 0.3481 ± 0.0866 0.36637 ± 0.0751 0.3005 ± 0.0064
h 0.7133 ± 0.09745 0.7413 ± 0.0829 0.6830 ± 0.0049
Ω r 5.16099 × 10 5 ± 1.5742 × 10 5 4.7019 × 10 5 ± 1.3266 × 10 5 5.3018 × 10 5 ± 7.6633 × 10 7
α 0.0679 ± 0.6644
b 0.04566 ± 0.1476
log Z 1784.16 ± 0.25 2355.57 ± 0.22 538.37 ± 0.215
log B with Λ CDM 1244.8 1817.2
Evidence over Λ CDMVery strong in favor of Λ CDMVery strong in favor of Λ CDM
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

García-Salcedo, R.; Gómez-Vargas, I.; González, T.; Martinez-Badenes, V.; Quiros, I. Combined Studies Approach to Rule Out Cosmological Models Which Are Based on Nonlinear Electrodynamics. Universe 2024, 10, 353. https://doi.org/10.3390/universe10090353

AMA Style

García-Salcedo R, Gómez-Vargas I, González T, Martinez-Badenes V, Quiros I. Combined Studies Approach to Rule Out Cosmological Models Which Are Based on Nonlinear Electrodynamics. Universe. 2024; 10(9):353. https://doi.org/10.3390/universe10090353

Chicago/Turabian Style

García-Salcedo, Ricardo, Isidro Gómez-Vargas, Tame González, Vicent Martinez-Badenes, and Israel Quiros. 2024. "Combined Studies Approach to Rule Out Cosmological Models Which Are Based on Nonlinear Electrodynamics" Universe 10, no. 9: 353. https://doi.org/10.3390/universe10090353

APA Style

García-Salcedo, R., Gómez-Vargas, I., González, T., Martinez-Badenes, V., & Quiros, I. (2024). Combined Studies Approach to Rule Out Cosmological Models Which Are Based on Nonlinear Electrodynamics. Universe, 10(9), 353. https://doi.org/10.3390/universe10090353

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