Next Article in Journal
Dynamic Analysis and FPGA Implementation of Fractional-Order Hopfield Networks with Memristive Synapse
Previous Article in Journal
Multi-Buoy Deployment Method Based on an Improved Tuna Swarm Optimizer Enhanced with Fractional-Order Calculus Method for Marine Observation
Previous Article in Special Issue
Estimated Age of the Universe in Fractional Cosmology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractional Einstein–Gauss–Bonnet Scalar Field Cosmology

1
Departamento de Física, Universidad Católica del Norte, Avenida Angamos 0610, Casilla 1280, Antofagasta 1270709, Chile
2
Departamento de Matemáticas, Universidad Católica del Norte, Avenida Angamos 0610, Casilla 1280, Antofagasta 1270709, Chile
3
Institute of Systems Science, Durban University of Technology, P.O. Box 1334, Durban 4000, South Africa
4
Department of Physics, Universidad de Antofagasta, Antofagasta 1240000, Chile
5
Escuela de Ingeniería, Universidad Central de Chile, Avenida Francisco de Aguirre 0405, La Serena 1710164, Chile
*
Author to whom correspondence should be addressed.
Fractal Fract. 2024, 8(11), 626; https://doi.org/10.3390/fractalfract8110626
Submission received: 1 October 2024 / Revised: 15 October 2024 / Accepted: 21 October 2024 / Published: 24 October 2024

Abstract

:
Our paper introduces a new theoretical framework called the Fractional Einstein–Gauss–Bonnet scalar field cosmology, which has important physical implications. Using fractional calculus to modify the gravitational action integral, we derived a modified Friedmann equation and a modified Klein–Gordon equation. Our research reveals non-trivial solutions associated with exponential potential, exponential couplings to the Gauss–Bonnet term, and a logarithmic scalar field, which are dependent on two cosmological parameters, m and α 0 = t 0 H 0 and the fractional derivative order μ . By employing linear stability theory, we reveal the phase space structure and analyze the dynamic effects of the Gauss–Bonnet couplings. The scaling behavior at some equilibrium points reveals that the geometric corrections in the coupling to the Gauss–Bonnet scalar can mimic the behavior of the dark sector in modified gravity. Using data from cosmic chronometers, type Ia supernovae, supermassive Black Hole Shadows, and strong gravitational lensing, we estimated the values of m and α 0 , indicating that the solution is consistent with an accelerated expansion at late times with the values α 0 = 1.38 ± 0.05 , m = 1.44 ± 0.05 , and μ = 1.48 ± 0.17 (consistent with Ω m , 0 = 0.311 ± 0.016 and h = 0.712 ± 0.007 ), resulting in an age of the Universe t 0 = 19.0 ± 0.7 [Gyr] at 1 σ CL. Ultimately, we obtained late-time accelerating power-law solutions supported by the most recent cosmological data, and we proposed an alternative explanation for the origin of cosmic acceleration other than Λ CDM. Our results generalize and significantly improve previous achievements in the literature, highlighting the practical implications of fractional calculus in cosmology.

1. Introduction

In general relativity, the Cosmological Principle is satisfied by the homogeneous and isotropic Friedmann–Lemaître–Robertson–Walker (FRLW) metric. It is believed that cold dark matter (CDM), which is not visible and only interacts through gravity, exists, and dark energy (in the form of a cosmological constant Λ ) explains the Universe’s accelerated expansion. This model is named Λ CDM and successfully explains the observed late-time acceleration of the Universe, as initially suggested by Type Ia supernovae (SNe Ia) observations [1] and later confirmed by Cosmic Microwave Background (CMB) measurements [2]. Moreover, it accurately describes the formation of the Universe’s large-scale structure. However, the model faces the long-standing Cosmological Constant problem [3,4], and the true mechanism driving the Universe’s accelerated expansion remains unknown [5]. Various alternative theories have been proposed to address these issues, including noncommutative theories, quantum cosmology, quantum deformed phase space models, and noncommutative minisuperspace approaches, in [6,7,8,9,10] and modified Brans–Dicke theory in [11,12]. Phantom fields introduce exotic physics through their negative kinetic energy and inherent quantum instabilities, though they are still supported by some observational evidence [13,14,15]. In previous studies, the cosmological evolution of phantom and quintom dark energy and their connection to late-time accelerated expansion were reviewed [16,17,18]. Constraints based on cosmic age and SNae I data and theoretical implications and observations can be found in [19,20]. The cosmological evolution of related hessence dark energy models was presented in [21,22]. The examination of crossing the phantom divide in the hybrid dark energy model context can be found in [23]. The reference [24] discusses cosmological perturbations of the quintom models of dark energy and their effects on observations. References [25,26] investigated the oscillatory behavior of quintom models. Additionally, ref. [27] studied interacting two-fluid scenarios for quintom dark energy, while ref. [28] presented a statefinder diagnostic for holographic dark energy model, and refs. [29,30] studied some holographic quintom dark energy models. Interestingly, ref. [31] proved that quintom models admit either tracking or phantom attractors, and refs. [32,33,34] analyzed the use of arbitrary potentials and classes of potential beyond the exponential potential. Furthermore, ref. [35] examined the transition from quintessence to phantom phase in the quintom model, and ref. [36] traced reconstructing attempts to relate these models with the Universe’s history, from inflation to acceleration, with phantom and canonical scalar fields. Additionally, refs. [37,38,39] extended the examination in [33,34]. The crossing of the phantom divide is also possible within the framework of scalar-tensor theories [40,41,42,43], as well as in various modified gravity theories [44,45,46,47,48,49]. Hořava–Lifshitz cosmology was investigated using phase-space analysis in [50] obtaining late-time bouncing oscillatory solutions under detailed balance. Recently, it was proven that the nonprojectable Hořava version is renormalizable [51,52]. In [53], the author investigated scaling solutions in the regime where the scalar field diverges, with an interest in describing the past asymptotic dynamics in nonminimally coupled dark energy.
Traditional mathematical models often need adjustments when dealing with power-law phenomena. Fractional calculus provides a suitable framework, accurately capturing non-local, frequency-dependent, and history-dependent properties inherent in power-law behavior. Fractional calculus is interdisciplinary and crucial for ongoing research. Recent studies show that fractional calculus is valuable in modeling fractional derivatives in quantum fields [54,55], as well as in exploring quantum gravity and cosmology [56,57,58,59,60], black holes [61,62], and fractional action cosmology [63,64,65,66,67]. We can develop fractional versions of Newtonian mechanics and Friedman–Robertson–Walker cosmology by substituting partial and fractional derivatives into familiar equations. For example, the study [68] investigated the emergence of a non-local-in-time fractional higher-order Newton’s second law of motion and the non-local-in-time fractional dynamics exhibiting disordered motions. Additionally, fractional variational problems have been presented in [69]. Fractional derivative methods have been developed through two primary approaches. The first, the last-step modification method, involves substituting the original cosmological field equations with fractional field equations tailored for a specific model. An illustrative example of this method is presented in [70]. The second approach, the first-step modification method, is more fundamental. In this method, fractional derivative geometry is established initially, followed by the application of the Fractional Action-like Variational Approach (FALVA) [71,72,73].
Using dynamical system methods, alongside observational data testing, provides a robust framework for analyzing the physical behavior of models, such as in fractional cosmology. This approach yields a cosmology with late-time acceleration without needing dark energy [74,75,76]. In the study by [74], researchers conducted a joint analysis using cosmic chronometers (CCs) and SNe Ia data to determine the best-fit values for the fractional order of the derivative. The observational tests for fractional cosmology in the study by [74] were improved in subsequent studies by [75,76]. In reference [77], the equation of state for a matter component was deduced based on compatibility conditions [77], which had not been previously analyzed in the study by [74].
In Gauss–Bonnet (GB) theory, the addition of the GB term changes the Einstein–Hilbert action as referenced in [78,79,80,81]. Studies such as [82,83] have demonstrated that scalar-coupled Gauss–Bonnet gravity in four dimensions can have significant effects on the early inflationary stage of our Universe. Various scenarios, including GB theory, F ( R ) , and Lorentz non-invariant models, are discussed in the reviews [48,49], providing unified descriptions of early-time inflation, bounce, and late-time cosmic acceleration. The mass of a scalar field is influenced when it interacts with the GB scalar through a non-constant coupling function, leading to non-zero effects on gravitational theory, as evidenced in [84,85,86,87]. In four dimensions, equilibrium points correspond to solutions with zero acceleration or de Sitter solutions in various coupling configurations, as found in [88,89]. A five-dimensional Einstein–Scalar–Gauss–Bonnet model was examined by [90], exploring different coupling functions and scalar field potentials while analyzing the dynamics and stability of equilibrium points. These models can explain the Universe’s early and late-time acceleration phases, making them suitable for studying inflation and dark energy. Reference [91] introduces a new theoretical framework for Einstein–Gauss–Bonnet theories, based on the requirement that the theory be compatible with the GW170817 event, ensuring that the gravitational wave speed is close to one in natural units.
In this paper, we introduce the Fractional Einstein–Gauss–Bonnet scalar field gravity framework, a novel model with significant cosmological implications. We present an evolution equation for an effective scalar field related to the coupling with the GB term. The research uncovers non-trivial solutions associated with exponential potentials, exponential couplings to the GB term, and logarithmic scalar fields. Stability conditions for the exact solutions are established using linear stability theory, and the scaling behavior of these solutions is investigated. Notably, the geometric corrections in the coupling to the GB scalar mimic the behavior of the dark sector in modified gravity. We estimate the free parameters using observational data from the SNe Ia, CC, gravitational lensing (GL), and Black Hole Shadow (BHS) data. Our solutions are consistent with current observational constraints and result in a late-time accelerating power-law solution for the scale factor. Additionally, we discuss the physical interpretation of the cosmological solutions, emphasizing the role of the fractional derivative order within our gravitational framework.
The manuscript is organized as follows. In Section 2, we will apply fractional differential calculus to compute specific physical quantities, highlighting their applications in cosmology. Rather than using covariant fractional derivatives as replacements for the standard covariant derivatives, we will adopt a point-like Lagrangian formulation of cosmology based on the flat FLRW metric and use FALVA to derive a modified cosmological model. Due to its relevance to cosmology, Section 3 explores scaling solutions and outlines a reconstruction procedure for both the potential and the coupling function. The stability analysis of the exact solution is conducted in Section 4. In Section 5, we will constrain the free parameters of the derived exact scaling solution. We will compute the best-fit parameters at the 1 σ of the confidence level for SNe Ia, CC, GL, and BHS data to achieve this. This approach offers the main advantages of combining dynamical systems analysis and observational testing, as it allows for using priors from the former analysis to inform the latter, aligning with the anticipated late-time behavior in a cosmological model. We present our conclusions in Section 6. For completeness, Appendix A contains our variational equations, while Appendix B outlines the formalism for the stability analysis of power-law solutions ψ s ( t ) of a generic ordinary differential equation given by F t , ψ ( t ) , ψ ˙ ( t ) , ψ ¨ ( t ) , 0 , where t > 0 is the independent variable and ψ ( t ) is the dependent variable [92,93,94]. Our proposal aims to construct regular differential equations on a bounded state space to describe dynamics globally, including the early- and late-time behavior of the intermediate stage of evolution, and alleviate cosmological problems like the H 0 tension and the Cosmic Coincidence problem.

2. Einstein–Gauss–Bonnet Scalar Field Gravity

We begin with an action with a scalar field with a non-zero coupling to the GB term,
S = d 4 x g R 2 κ 2 1 2 g μ ν μ ϕ ν ϕ V ( ϕ ) f ( ϕ ) G ,
where R is the curvature scalar and we have assumed units where κ 2 = 8 π G = 1 , and G is the GB term:
G = R 2 4 R α β R α β + R α β γ δ R α β γ δ .
Here, ϕ represents a scalar field that is minimally coupled to gravity and possesses a self-interacting potential V ( ϕ ) . Additionally, it is nonminimally coupled to the GB invariant G through the coupling function f ( ϕ ) . We are considering the FLRW
d s 2 = N 2 ( t ) d t 2 + a 2 ( t ) d x 2 + d y 2 + d z 2 .
Defining H = a ˙ / ( N a ) we obtain that the curvature scalar and the GB term are given by
R = 12 H 2 + 6 H ˙ / N , G = 24 H 2 H 2 + H ˙ / N .
By substituting R and G into the action, we can find the point-like Lagrangian:
L = N ˙ ( θ ) 24 a ˙ ( θ ) 3 f ( ϕ ( θ ) ) N ( θ ) 4 3 a ( θ ) 2 a ˙ ( θ ) N ( θ ) 2 + 3 a ( θ ) a ˙ ( θ ) 2 N ( θ ) + a ¨ ( θ ) 3 a ( θ ) 2 N ( θ ) 24 a ˙ ( θ ) 2 f ( ϕ ( θ ) ) N ( θ ) 3 + a ( θ ) 3 ϕ ˙ ( θ ) 2 2 N ( θ ) a ( θ ) 3 N ( θ ) V ( ϕ ( θ ) ) .
Extending this model using fractional calculus, we consider that the action is defined by
S ( τ ) = 1 Γ ( μ ) 0 τ L ( θ , q i ( θ ) , q ˙ i ( θ ) , q ¨ i ( θ ) ) ( τ θ ) μ 1 d θ .
We follow the procedures of reference [95], which start with the action (6). We make the steps in Appendix A to obtain the equations of motion of the fields. Implementing the derivatives in (A7), we have
L ( θ , q i ( θ ) , q ˙ i ( θ ) , q ¨ i ( θ ) ) q i d d θ L ( θ , q i ( θ ) , q ˙ i ( θ ) , q ¨ i ( θ ) ) q ˙ i + d 2 d θ 2 L ( θ , q i ( θ ) , q ˙ i ( θ ) , q ¨ i ( θ ) ) q ¨ 1 = 1 μ τ θ L ( θ , q i ( θ ) , q ˙ i ( θ ) , q ¨ i ( θ ) ) q ˙ i 2 d d θ L ( θ , q i ( θ ) , q ˙ i ( θ ) , q ¨ i ( θ ) ) q ¨ i ( 1 μ ) ( 2 μ ) ( τ θ ) 2 L ( θ , q i ( θ ) , q ˙ i ( θ ) , q ¨ i ( θ ) ) q ¨ i .
Then, by applying Hamilton’s principle, we obtain the Euler–Poisson equations modified by the fractional parameter μ given by (7) by taking the variations with respect to q i { N , a , ϕ } , and making N = 1 after the variation. Then, using the parameterization ( τ , θ ) = ( 2 t , t ) , we obtain respectively the Friedmann, Raychaudhuri and Klein–Gordon equations:
24 a ˙ 3 ϕ ˙ f ( ϕ ) a 3 24 ( μ 1 ) a ˙ 3 f ( ϕ ) t a 3 3 a ˙ t a ˙ μ a + a t a 2 + V ( ϕ ) + 1 2 ϕ ˙ 2 = 0 , 2 a ¨ a + 8 a ˙ 2 ϕ ˙ 2 f ( ϕ ) a 2 + a ˙ 2 ( μ 1 ) a t a ˙ t a 2 + f ( ϕ ) 16 ( μ 1 ) a ˙ 2 ϕ ˙ t a 2 + 8 a ˙ 2 ϕ ¨ a 2 + 16 a ˙ a ¨ ϕ ˙ a 2
+ f ( ϕ ) 8 ( μ 2 ) ( μ 1 ) a ˙ 2 t 2 a 2 16 ( μ 1 ) a ˙ a ¨ t a 2 ( μ 2 ) ( μ 1 ) t 2 + V ( ϕ ) 1 2 ϕ ˙ 2 = 0 ,
3 a ˙ ϕ ˙ a 24 a ˙ 2 a ¨ f ( ϕ ) a 3 + ( μ 1 ) ϕ ˙ t V ( ϕ ) ϕ ¨ = 0 .
Substituting H = a ˙ / a (with N 1 ) in (8)–(10), we have
24 H 3 ϕ ˙ f ( ϕ ) 24 ( μ 1 ) H 3 f ( ϕ ) t + 3 ( μ 1 ) H t 3 H 2 + V ( ϕ ) + 1 2 ϕ ˙ 2 = 0 , 8 H 2 ϕ ˙ 2 f ( ϕ ) + f ( ϕ ) 16 H H ˙ ϕ ˙ + H 2 8 ϕ ¨ 16 ( μ 1 ) ϕ ˙ t + 16 H 3 ϕ ˙ + f ( ϕ ) 16 ( μ 1 ) H H ˙ t + 8 ( μ 2 ) ( μ 1 ) H 2 t 2 16 ( μ 1 ) H 3 t
2 H ˙ + 2 ( μ 1 ) H t 3 H 2 ( μ 2 ) ( μ 1 ) t 2 + V ( ϕ ) 1 2 ϕ ˙ 2 = 0 ,
24 H 2 H ˙ f ( ϕ ) 24 H 4 f ( ϕ ) + ( 3 t H + μ 1 ) ϕ ˙ t V ( ϕ ) ϕ ¨ = 0 .
By solving the system (11)–(13) for V ( ϕ ) , ϕ ¨ and H ˙ , we obtain the following equations:
V ( ϕ ) = ( 3 3 μ ) H t + 3 H 2 + 24 ( 1 + μ ) f ( ϕ ) H 3 t 24 H 3 f ( ϕ ) ϕ ˙ 1 2 ϕ ˙ 2 , ϕ ¨ = 12 t ( 1 + μ ) H 3 f ( ϕ ) + t t V ( ϕ ) + ( 1 + μ ) ϕ ˙ + 288 t H 5 f ( ϕ ) ( 1 μ ) f ( ϕ ) + t f ( ϕ ) ϕ ˙ + 12 H 2 2 t ( 1 + μ ) f ( ϕ ) ϕ ˙ + f ( ϕ ) 2 3 μ + μ 2 + 3 t 2 ϕ ˙ 2 + H 8 ( 1 + μ ) f ( ϕ ) t V ( ϕ ) ( 1 + μ ) ϕ ˙ + t ϕ ˙ 3 t + 8 f ( ϕ ) t V ( ϕ ) ( 1 + μ ) ϕ ˙ 24 H 4 f ( ϕ ) 4 ( 2 + μ ) ( 1 + μ ) f ( ϕ ) + t t + 4 ϕ ˙ 2 ( 1 + μ ) f ( ϕ ) + t ϕ ˙ f ( ϕ )
t t + 8 H ( 1 + μ ) f ( ϕ ) + t f ( ϕ ) 12 H 3 f ( ϕ ) ϕ ˙ , H ˙ = 2 3 μ + μ 2 + t ( 1 + μ ) H + 192 t 2 H 6 f ( ϕ ) 2 + t 2 ϕ ˙ 2 + 8 t H 3 ( 1 μ ) f ( ϕ ) + 4 t f ( ϕ ) ϕ ˙ 8 H 2 ( 2 + μ ) ( 1 + μ ) f ( ϕ ) + t f ( ϕ ) t V ( ϕ ) ( 1 + μ ) ϕ ˙ + t ϕ ˙ 2 f ( ϕ )
2 t t + 8 H ( 1 + μ ) f ( ϕ ) + t f ( ϕ ) 12 H 3 f ( ϕ ) ϕ ˙ .
We now impose that the time derivative of Friedmann Equation (11) is zero under the hypothesis that μ 1 (for μ = 1 is trivial); that is,
d d t 24 H 3 ϕ ˙ f ( ϕ ) 24 ( μ 1 ) H 3 f ( ϕ ) t + 3 ( μ 1 ) H t 3 H 2 + V ( ϕ ) + 1 2 ϕ ˙ 2 = 0 .
Computing the time derivative and replacing V ( ϕ ) , ϕ ¨ and H ˙ given by (14)–(16) into (17), we obtain
0 = { 3 ( 2 + μ ) ( 1 + μ ) + t 2 ϕ ˙ 2 4608 t H 7 f ( ϕ ) 2 ( ( 3 + μ ) f ( ϕ ) ) + t f ( ϕ ) ϕ ˙ + 48 t H 3 ( ( 1 + μ ) f ( ϕ ) ) + 2 t f ( ϕ ) ϕ ˙ + t H 15 3 μ 16 ( 1 + μ ) f ( ϕ ) ϕ ˙ 2 + 16 t f ( ϕ ) ϕ ˙ 3 576 t H 5 ( 1 + μ ) f ( ϕ ) 2 + ( 3 + μ ) f ( ϕ ) 2 + 2 t f ( ϕ ) f ( ϕ ) ϕ ˙ 48 H 4 [ 4 ( 1 + μ ) ( 8 + 3 μ ) f ( ϕ ) 2 4 t 2 f ( ϕ ) 2 ϕ ˙ 2 + t f ( ϕ ) t + 4 f ( ϕ ) t V ( ϕ ) + ( 9 + 5 μ ) ϕ ˙ 4 t ϕ ˙ 2 f ( ϕ ) ] 6 H 2 4 f ( ϕ ) 10 14 μ + 4 μ 2 + t 2 ϕ ˙ 2 + t t 4 f ( ϕ ) t V ( ϕ ) + ( 7 + 3 μ ) ϕ ˙ + 4 t ϕ ˙ 2 f ( ϕ ) } / 2 t 2 t + 8 H ( 1 + μ ) f ( ϕ ) + t f ( ϕ ) 12 H 3 f ( ϕ ) ϕ ˙ .
We introduce the auxiliary field
ψ = f ( ϕ ) .
Moreover, we calculate the successive derivatives using the chain rule:
ψ ˙ = d d t f ( ϕ ) = f ( ϕ ) ϕ ˙ , ψ ¨ = f ( ϕ ) ϕ ˙ 2 + f ( ϕ ) ϕ ¨ = f ( ϕ ) ϕ ˙ 2 + ψ ˙ ϕ ˙ ϕ ¨ .
In this way, we obtain
f ( ϕ ) = ψ ˙ ϕ ˙ , f ( ϕ ) = 1 ϕ ˙ 2 ψ ¨ ψ ˙ ϕ ˙ ϕ ¨ .
Replacing the expression (21) into (14)–(16) and (18), then solving the resulting system for ϕ ¨ , ψ ¨ , H ˙ and V ( ϕ ) , it turns out that
ϕ ¨ = 24 ( μ 3 ) H 3 ψ ˙ ϕ ˙ + ( μ 1 ) ϕ ˙ t + 192 H 5 ψ ˙ 2 1 8 H 2 ψ ϕ ˙ + ϕ ˙ 8 H 2 ψ ˙ 1 8 H 2 ψ 3 H V ( ϕ ) , ψ ¨ = ( μ 1 ) μ 2 H 2 + 8 ( 8 3 μ ) ψ 8 t 2 + 2 ψ ˙ μ 1 8 H 2 ψ 1 + 3 μ 5 3 ( μ 1 ) H ψ + 2 ( μ 1 ) ψ ϕ ˙ 2 3 H 8 H 2 ψ 1 μ 5 8 H t
+ 16 H 2 ψ ˙ 2 1 8 H 2 ψ + H 2 8 H 2 ψ 1 + 3 ψ ˙ + ϕ ˙ 2 2 ψ ˙ 3 H 24 H 3 ψ + 2 8 H 2 ψ 1 + 3 24 H 2 1 4 ,
H ˙ = ( μ 3 ) H t + 8 H 3 ψ ˙ 8 H 2 ψ 1 + ϕ ˙ 2 24 H 2 ψ 3 H 2 ,
V ( ϕ ) = 3 ( μ 1 ) H 8 H 2 ψ 1 t 24 H 3 ψ ˙ + 3 H 2 1 2 ϕ ˙ 2 .
From Equation (25), we define the effective energy densities as follows:
ρ ϕ = 1 2 ϕ ˙ 2 + V ( ϕ ) , ρ fracc = 3 ( μ 1 ) t H , ρ GB = 24 H 3 ψ ˙ , ρ GB , fracc = 24 ( μ 1 ) t H 3 ψ ,
such that (11) is reduced to
3 H 2 = ρ ϕ + ρ fracc + ρ GB + ρ GB , fracc .

3. Exact Scaling Solutions

We can find an analytical solution to the system by considering
H = 2 m t p ,
ψ = F 0 t q ,
ϕ = 2 ln ( t ) λ , λ 0 ,
where m, F 0 , p, q and λ are constants.
Then, Equations (22)–(25) are written as
64 F 0 ( q μ ) t q + 2 μ m 2 t 2 p λ t 2 m 2 t 2 p 32 F 0 t q 96 F 0 λ q ( μ + q + 3 ) t q 3 p 1 m 3 + 3 λ m q 2 t p 1 m 2 t 2 p 32 F 0 t q
3 t p 1 λ 2 q 2 + 4 λ m V ln t 2 / λ = 0 , m 4 3 λ 2 ( μ 2 ) ( μ 1 ) + 4 t 4 p 2 96 λ 2 m 2 t 2 p 32 F 0 t q 1 4 + 32 F 0 2 11 μ + 3 μ 2 2 μ q + q ( q + 3 ) + 8 t 2 q 2 m 2 t 2 p 32 F 0 t q F 0 m 2 t 2 p + q 2 λ 2 ( 2 μ + q + 2 ) ( 2 μ + q + 5 ) + 4 λ 2 m 2 t 2 p 32 F 0 t q + 4 F 0 m 3 ( μ + q + 1 ) t 3 p + q 3 3 λ 2 m 2 t 2 p 32 F 0 t q
+ 1 16 m t p 1 64 F 0 q t q m 2 t 2 p 32 F 0 t q μ + 5 + 6 F 0 ( μ + q + 1 ) t q p 1 m = 0 ,
4 3 λ 2 t 2 32 F 0 t q 2 p m 2 1 + 2 t p 1 ( μ + p + q + 3 ) m 2 m q t p 1 m 2 t 2 p 32 F 0 t q 4 t 2 p m 2 = 0 ,
6 t 3 p 1 32 F 0 ( μ q 1 ) t q ( μ 1 ) m 2 t 2 p m 3 + 12 t 2 p m 2 V ln t 2 / λ 2 λ 2 t 2 = 0 .
By comparing the coefficients of t in the last equation, we have, by dimensional analysis, that
q = 2 p , 2 p = 2 p = 1 , q = 2 .
Then, for p = 1 and q = 2 , and assuming m 2 32 F 0 0 , m 0 and λ 0 we have
6144 F 0 2 λ 2 ( μ 5 ) + 64 F 0 m 2 3 λ 2 ( μ 3 ) + ( μ 2 ) m 6
+ λ m 3 t 2 m 2 32 F 0 V ln t 2 / λ 2 m 4 ( μ m 6 ) = 0 , 3072 F 0 2 λ 2 6 ( μ 3 ) + ( μ ( 3 μ 23 ) + 38 ) m 64 F 0 m 12 λ 2 + 2 ( μ 3 ) m 3 + m 2 3 λ 2 ( μ 2 ) ( 2 μ 7 ) + 6 + 6 λ 2 μ m
+ m 3 24 λ 2 + m 2 3 λ 2 ( μ 2 ) ( μ 1 ) + 4 6 λ 2 ( μ 5 ) m = 0 ,
192 F 0 λ 2 + 96 F 0 λ 2 ( μ 6 ) m 2 m 4 3 λ 2 ( μ 4 ) m 3 6 λ 2 m 2 = 0 ,
2 96 F 0 λ 2 ( μ 3 ) + m 3 + 3 λ 2 ( μ 1 ) m 2 6 λ 2 m λ 2 m 3 t 2 V ln t 2 / λ = 0 .
Equation (39) allows us to reconstruct the GB coupling and scalar field potentials, which replicate scaling behaviors.
In our scenario, say we acquire the exponential potential
V ϕ = 2 e λ ϕ 96 F 0 λ 2 ( μ 3 ) m 3 3 λ 2 ( μ 1 ) m 2 + 6 λ 2 m λ 2 m 3 .
Replacing the derivative of potential (40) in Equation (36), and assuming ( μ 6 ) m + 2 0 and m 0 , results in two independent equations:
F 0 = 1 96 m 2 2 m 3 λ 2 + m λ 2 ( ( μ 6 ) m + 2 ) + 3 ,
0 = 3 λ 2 + m m 3 λ 2 ( μ + 2 ) 4 μ + ( ( μ 5 ) μ + 2 ) m + 30 6 3 λ 2 + 2 .
The first condition is the definition of F 0 . Hence, we have the two-parametric family of solutions (recall λ = λ ( m , μ ) ):
R ( t ) = 12 ( m 4 ) m 2 t 2 Using ( 30 ) R ( ϕ ) = 12 ( m 4 ) e λ ϕ m 2 ,
G ( t ) = 192 ( m 2 ) m 4 t 4 G ( ϕ ) = 192 ( m 2 ) e 2 λ ϕ m 4 ,
H ( t ) = 2 m t H ( ϕ ) = 2 m e λ ϕ 2 , ψ ( t ) = 1 96 m 2 t 2 2 m 3 λ 2 + m λ 2 ( μ 6 ) m + 2 + 3
f ( ϕ ) = 1 96 m 2 e λ ϕ 2 m 3 λ 2 + m λ 2 ( μ 6 ) m + 2 + 3 ,
V ( t ) = 2 12 λ 2 + μ m 3 + 2 9 λ 2 1 m 2 + 6 λ 2 ( μ 8 ) m λ 2 m 2 t 2 ( μ 6 ) m + 2
V ( ϕ ) = 2 e λ ϕ 12 λ 2 + μ m 3 + 2 9 λ 2 1 m 2 + 6 λ 2 ( μ 8 ) m λ 2 m 2 ( μ 6 ) m + 2 .
Replacing the derivative of potential (40) in Equation (36), we obtain the conditions on the parameters
6144 F 0 2 λ 2 ( μ 4 ) + 32 F 0 m 6 λ 2 + m 3 λ 2 ( 7 3 μ ) + ( μ 3 ) m 6 + m 3 m 3 λ 2 ( μ 1 ) μ m + m + 6 6 λ 2 = 0 ,
together with (37) and (38). Reducing the conditions and assuming ( μ 6 ) m + 2 0 and m 0 , we obtain the cases
  • For 0 < μ 1 2 5 17 and 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 < m < 6 μ + 2 ,
    λ = 1 3 36 3 m 4 μ + ( ( μ 5 ) μ + 2 ) m + 30 ( μ + 2 ) m 6 and F 0 = m 2 ( m ( 4 μ + ( μ 2 ) ( μ 1 ) m + 18 ) 12 ) 32 ( m ( 4 μ + ( ( μ 5 ) μ + 2 ) m + 30 ) 12 ) . For μ = 1 2 5 17 and 1 83 60 6 17 < m < 1 16 27 + 3 17 , λ = 2 3 3 10 + 17 m 18 17 9 m + 12 and F 0 = m 2 17 5 m 2 2 4 + 17 m + 12 64 10 + 17 m 6 .
  • For 1 2 5 17 < μ 4 and m > 2 μ 15 μ 2 5 μ + 2 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 we have
    λ = 36 3 m 4 μ + ( ( μ 5 ) μ + 2 ) m + 30 3 ( μ + 2 ) m 6 and F 0 = m 2 ( m ( 4 μ + ( μ 2 ) ( μ 1 ) m + 18 ) 12 ) 32 ( m ( 4 μ + ( ( μ 5 ) μ + 2 ) m + 30 ) 12 ) .
  • For 4 < μ 1 2 5 + 17 and 6 μ + 2 < m < 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 we have λ = 36 3 m 4 μ + ( ( μ 5 ) μ + 2 ) m + 30 3 ( μ + 2 ) m 6 and F 0 = m 2 ( m ( 4 μ + ( μ 2 ) ( μ 1 ) m + 18 ) 12 ) 32 ( m ( 4 μ + ( ( μ 5 ) μ + 2 ) m + 30 ) 12 ) . For μ = 1 2 5 + 17 and 1 16 27 3 17 < m < 1 83 60 + 6 17 we have
    λ = 2 17 10 m + 6 3 9 + 17 m 36 and F 0 = m 2 12 m 5 + 17 m 2 17 + 8 64 17 10 m + 6 .
  • For 4 < μ < 1 2 5 + 17 and m > 2 μ 15 μ 2 5 μ + 2 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 we have
    λ = 36 3 m 4 μ + ( ( μ 5 ) μ + 2 ) m + 30 3 ( μ + 2 ) m 6 and F 0 = m 2 ( m ( 4 μ + ( μ 2 ) ( μ 1 ) m + 18 ) 12 ) 32 ( m ( 4 μ + ( ( μ 5 ) μ + 2 ) m + 30 ) 12 ) .
  • For μ > 1 2 5 + 17 and 6 μ + 2 < m < 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 we have λ = 1 3 36 3 m 4 μ + ( ( μ 5 ) μ + 2 ) m + 30 ( μ + 2 ) m 6 and F 0 = m 2 ( m ( 4 μ + ( μ 2 ) ( μ 1 ) m + 18 ) 12 ) 32 ( m ( 4 μ + ( ( μ 5 ) μ + 2 ) m + 30 ) 12 ) .
In Figure 1, we see the allowed regions for the parameters μ and m.
Replacing the ansatz
H = 2 m t , m 6 μ ,
in Equations (22)–(25), we obtain
ϕ ¨ = 6 ( m 2 ) ( μ 4 ) m + 2 m 2 t 2 32 ψ m 5 t 5 ϕ ˙ + ( μ 2 ) m 4 ϕ ˙ m t V ( ϕ ( t ) ) , ψ ¨ = m ( μ ( ( μ 7 ) m + 4 ) 6 ) + 4 ψ m 2 t 2
+ 1 96 m 2 t 2 ( 2 μ m ) ϕ ˙ 2 + 3 m ( 4 μ + ( 2 ( μ 5 ) μ ) m + 2 ) 12 ,
ψ ˙ = ( μ 4 ) m + 2 ψ m t 1 192 m t m 2 t 2 ϕ ˙ 2 + 6 ( μ 4 ) m + 12 ,
V ( ϕ ) = 6 ( 4 3 m ) m 2 t 2 + 192 ( 3 m 2 ) ψ m 4 t 4 + 1 2 ϕ ˙ 2 .
To investigate the stability of the solution (30) and (45)–(47), we specify the potential (48), obtaining the equations
ϕ ¨ = 192 λ ( 3 m 2 ) ψ m 4 t 4 + 6 λ ( 4 3 m ) m 2 t 2 + 192 ( m 2 ) ( μ 4 ) m + 2 ψ m 5 t 5 6 ( m 2 ) ( μ 4 ) m + 2 m 3 t 3 ϕ ˙
+ ( μ 2 ) m 4 ϕ ˙ m t + 1 2 λ ϕ ˙ 2 , ψ ¨ = m ( μ ( ( μ 7 ) m + 4 ) 6 ) + 4 ψ m 2 t 2
+ 1 96 m 2 t 2 ( 2 μ m ) ϕ ˙ 2 + 3 m ( 4 μ + ( 2 ( μ 5 ) μ ) m + 2 ) 12 ,
ψ ˙ = ( μ 4 ) m + 2 ψ m t 1 192 m t m 2 t 2 ϕ ˙ 2 + 6 ( μ 4 ) m + 12 .
After reducing the quadratic term ϕ ˙ 2 / 2 using the last equation, we obtain
ϕ ¨ = 96 λ ψ ˙ m 3 t 3 6 ( m 2 ) ( μ 4 ) m + 2 m 3 t 3 ϕ ˙ 3 λ ( μ + 2 ) m 6 m 2 t 2
+ ψ 192 ( m 2 ) ( μ 4 ) m + 2 m 5 t 5 ϕ ˙ + 96 λ ( μ + 2 ) m 2 m 4 t 4 + ( μ 2 ) m 4 ϕ ˙ m t , ψ ¨ = 1 32 m ( 4 μ + ( μ 2 ) ( μ 1 ) m + 18 ) 12
+ m ( μ ( μ m + m + 4 ) 22 ) + 12 ψ m 2 t 2 + 2 μ 2 m ψ ˙ t ,
ϕ ˙ 2 = 192 ( μ 4 ) m + 2 ψ m 4 t 4 192 ψ ˙ m 3 t 3 6 ( μ 4 ) m + 2 m 2 t 2 .
Equation (60) defines ϕ ˙ , and we analyze (59), our equation of interest. The general solution of (59) is
ψ ( t ) = m 2 t 2 m ( 4 μ + ( μ 2 ) ( μ 1 ) m + 18 ) 12 32 m ( 4 μ + ( ( μ 5 ) μ + 2 ) m + 30 ) 12 + c 1 t 2 μ m + m 4 m ( 8 μ m + m 96 ) + 64 2 m + c 2 t 2 μ m + m 4 + m ( 8 μ m + m 96 ) + 64 2 m .
Taking ψ ( t ) from (61), calculating its time derivative, ψ ˙ ( t ) , and substituting in (60), ϕ ( t ) is given through the quadrature
ϕ ( t ) = ± 1 t 6 m ( μ 4 ) + 2 m 2 η 2 32 ψ ( η ) 192 m η ψ ( η ) m 2 η 2 d η + c 3 .

4. Stability Analysis of the Exact Solution

Denoting
ψ c ( t ) = 1 96 m 2 t 2 2 m 3 λ 2 + m λ 2 ( μ 6 ) m + 2 + 3 ,
ϕ c ( t ) = 2 ln ( t ) λ .
In this case, the evolution equation for ψ is decoupled.
Specializing the procedure in Appendix B to the system (59) and (60), we define
ε ( τ ) = ψ ( τ ) ψ c ( τ ) 1 .
Hence,
ψ ¨ ( t ) = m 2 e 2 τ 6 λ 2 + 2 m 2 + 3 λ 2 ( μ 4 ) m ε ( τ ) + 3 ε ( τ ) + 2 ε ( τ ) + 2 96 λ 2 t 2 ( μ 6 ) m + 2 ,
ψ ˙ ( t ) = m 2 e 2 τ 6 λ 2 + 2 m 2 + 3 λ 2 ( μ 4 ) m ε ( τ ) + 2 ε ( τ ) + 2 96 λ 2 t ( μ 6 ) m + 2 ,
ψ ( t ) = m 2 e 2 τ ε ( τ ) + 1 6 λ 2 + 2 m 2 + 3 λ 2 ( μ 4 ) m 96 λ 2 ( μ 6 ) m + 2 .
Combining with Equation (59), we obtain
ε ( τ ) = 36 λ 2 2 m 3 λ 2 ( μ + 2 ) 4 μ + ( ( μ 5 ) μ + 2 ) m + 30 + 24 6 λ 2 + 2 m 2 + 3 λ 2 ( μ 4 ) m = 0 by ( 42 ) + 12 m ( 4 μ + ( ( μ 5 ) μ + 2 ) m + 30 ) m 2 ε ( τ ) + 2 μ 4 m 3 ε ( τ ) ,
where using the relation (42), the first term is zero. Defining v ( τ ) = ε ( τ ) , we obtain the linear dynamical system
ε = v ,
v = 12 m ( 4 μ + ( ( μ 5 ) μ + 2 ) m + 30 ) m 2 ε + 2 μ 4 m 3 v .

4.1. Stability of the Scaling Solution

The scaling solution corresponds to the fixed point ( ε , v ) = ( 0 , 0 ) . The matrix of the linear system is defined by
J ( ε , v ) = 0 1 12 m 4 μ + ( ( μ 5 ) μ + 2 ) m + 30 m 2 2 μ 4 m 3 .
Evaluating the linear matrix around the fixed point ( ε = 0 , v = 0 ) , we obtain the eigenvalues
λ 1 , 2 = 1 2 2 μ 4 m 3 ± m ( 8 μ m + m 96 ) + 64 m 2 .
The origin is a sink for
  • 0 < μ < 1 2 5 17 , 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 < m 48 8 μ + 1 8 8 μ 35 ( 8 μ + 1 ) 2 , or
  • 0 < μ < 1 2 5 17 , m 8 8 μ 35 ( 8 μ + 1 ) 2 + 48 8 μ + 1 , or
  • μ = 1 2 5 17 , 6 10 17 < m 48 1 + 4 5 17 8 35 4 5 17 1 + 4 5 17 , or
  • μ = 1 2 5 17 , m 48 1 + 4 5 17 + 8 35 4 5 17 1 + 4 5 17 , or
  • 1 2 5 17 < μ < 287 + 146 3 148 + 80 3 , 2 μ 15 μ 2 5 μ + 2 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 < m 48 8 μ + 1 8 8 μ 35 ( 8 μ + 1 ) 2 , or
  • 1 2 5 17 < μ < 287 + 146 3 148 + 80 3 , 8 8 μ 35 ( 8 μ + 1 ) 2 + 48 8 μ + 1 m < 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 , or
  • 287 + 146 3 148 + 80 3 μ < 146 3 287 80 3 148 , 2 μ 15 μ 2 5 μ + 2 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 < m 48 8 μ + 1 8 8 μ 35 ( 8 μ + 1 ) 2 .
The origin is a source for
  • 287 + 146 3 148 + 80 3 < μ 146 3 287 80 3 148 , 8 8 μ 35 ( 8 μ + 1 ) 2 + 48 8 μ + 1 m < 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 , or
  • 146 3 287 80 3 148 < μ < 35 8 , 2 μ 15 μ 2 5 μ + 2 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 < m 48 8 μ + 1 8 8 μ 35 ( 8 μ + 1 ) 2 , or
  • 146 3 287 80 3 148 < μ < 35 8 , 8 8 μ 35 ( 8 μ + 1 ) 2 + 48 8 μ + 1 m < 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 , or
  • 35 8 μ < 1 2 5 + 17 , 2 μ 15 μ 2 5 μ + 2 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 < m < 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 , or
  • μ = 1 2 5 + 17 , m > 6 17 10 , or
  • μ > 1 2 5 + 17 , m > 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 .
The origin is a saddle for
  • 0 < μ < 1 2 5 17 , 0 < m < 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 , or
  • μ = 1 2 5 17 , 0 < m < 6 10 17 , or
  • 1 2 5 17 < μ < 1 2 5 + 17 , 0 < m < 2 μ 15 μ 2 5 μ + 2 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 , or
  • 1 2 5 17 < μ < 1 2 5 + 17 , m > 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2
  • μ = 1 2 5 + 17 , 0 < m < 6 17 10 , or
  • μ > 1 2 5 + 17 , 0 < m < 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 .
In Figure 2, we can see the allowed regions for the parameters μ and m for the sink, source, and saddle cases.
Now, we will show the flow of system (70) and (71). In Figure 3, we show some phase space diagrams for sink cases, in Figure 4 the source cases, and in Figure 5, the saddle cases. We used specific numerical values of m and μ representing the regions depicted in Figure 2. These values are chosen to meet the specified constraints. These regions are then utilized to determine the priors for the Bayesian analysis carried out in Section 5, where we calculate the best-fit parameters at the 1 σ confidence level (CL) using the affine-invariant Markov chain Monte Carlo (MCMC) method [96], which is implemented in the Python code emcee version 3.1.6 [97].
This approach offers the main advantages of combining dynamical systems analysis and observational testing, as it allows for using priors from the former analysis to inform the latter, aligning with the anticipated late-time behavior in a cosmological model.

4.2. Dynamics at Infinity

To investigate the late-time dynamics of system (70) and (71) at infinity, we define the following variables:
X = ε 1 + ε 2 + v 2 , Y = v 1 + ε 2 + v 2 .
Hence, we obtain the system
d X d τ = Y m 2 m m X ( ( μ 5 ) μ X + X + ( 3 2 μ ) Y ) + m + 2 X ( ( 15 2 μ ) X + 2 Y ) 12 X 2 , d Y d τ = 1 m 2 [ X m 2 ( μ 5 ) μ ( ( μ 5 ) μ + 1 ) Y 2 + 2
+ 2 ( 2 μ 15 ) m Y 2 1 + 12 Y 2 1 + m Y Y 2 1 ( ( 2 μ 3 ) m 4 ) ] ,
defined in the compact set { ( X , Y ) R 2 : X 2 + Y 2 1 } .
There are five total equilibrium points for the system (75) and (76). However, the existence conditions for each point and their respective eigenvalues depend on both μ and m; therefore, the analysis is not trivial. The stationary points in the coordinates ( X , Y ) are
1.
A = ( 0 , 0 ) , with eigenvalues (73). The stability analysis is the same as in the previous finite regime.
2.
B , C = ± k 2 m ( k 1 + ( 2 μ 3 ) m 4 ) k 2 m ( 4 μ + ( ( μ 5 ) μ + 2 ) m + 30 ) 24 , ± k 2 k are antipodal points, where
k = 2 μ 4 10 μ 3 + 31 μ 2 22 μ + 10 m 4   8 2 μ 3 25 μ 2 + 81 μ 21 m 3 8 2 μ 2 + 30 μ 223 m 2 + 96 ( 2 μ 15 ) m + 288 ,
k 1 = m ( 8 μ m + m 96 ) + 64 ,
k 2 = m 3 ( 3 ( k 1 + 68 ) 2 μ ( k 1 + 4 μ ( 2 μ 25 ) + 320 ) ) + 4 m 2 ( k 1 4 μ ( μ + 15 ) + 436 ) + ( 2 ( μ 1 ) μ ( ( μ 9 ) μ + 21 ) + 13 ) m 4 + 96 ( 2 μ 15 ) m + 288 .
These points exist for the following intervals
(1)
μ = 1 2 5 17 , 0 m < 1 83 60 6 17 , or
(2)
μ = 1 2 5 17 , 1 83 60 6 17 < m < 8 15 + 4 17 48 4 17 21 , or
(3)
μ = 1 2 5 17 , m > 48 8 15 + 4 17 4 17 21 , or
(4)
35 8 < μ < 1 2 5 + 17 , m > 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 , or
(5)
1 4 11 2 3 < μ < 1 4 11 + 2 3 , m > 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 , or
(6)
1 2 5 17 < μ < 1 4 11 2 3 , m > 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 , or
(7)
μ > 1 2 5 + 17 , m > 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 , or
(8)
1 4 11 + 2 3 < μ 35 8 , m > 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 , or
(9)
m = 8 13 + 4 3 6 4 3 23 , μ = 1 4 11 2 3 , or
(10)
μ = 1 4 11 2 3 , 0 m < 152 + 16 3 32 13 + 4 3 55 + 4 3 , or
(11)
μ = 1 4 11 2 3 , 152 + 16 3 32 13 + 4 3 55 + 4 3 < m < 8 13 + 4 3 48 4 3 23 , or
(12)
μ = 1 4 11 2 3 , m > 48 8 13 + 4 3 4 3 23 , or
(13)
m = 8 6 + 13 4 3 23 + 4 3 , μ = 1 4 11 + 2 3 , or
(14)
μ = 1 4 11 + 2 3 , 0 m < 48 8 13 4 3 23 + 4 3 , or
(15)
μ = 1 4 11 + 2 3 , 48 + 8 13 4 3 23 + 4 3 < m < 152 + 16 3 32 13 4 3 4 3 55 , or
(16)
μ = 1 4 11 + 2 3 , m > 152 + 16 3 32 13 4 3 4 3 55 , or
(17)
μ = 1 2 5 + 17 , 0 m < 6 17 10 , or
(18)
μ = 1 2 5 + 17 , m > 6 17 10 , or
(19)
μ = 35 8 , m = 4 3 , or
(20)
μ > 1 2 5 + 17 , 0 m < 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 , or
(21)
0 μ < 1 2 5 17 , 0 m < 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 , or
(22)
1 4 11 + 2 3 < μ 35 8 , 0 m < 2 μ 15 μ 2 5 μ + 2 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 , or
(23)
35 8 < μ < 1 2 5 + 17 , 0 m < 2 μ 15 μ 2 5 μ + 2 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 , or
(24)
1 2 5 17 < μ < 1 4 11 2 3 , 0 m < 2 μ 15 μ 2 5 μ + 2 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 , or
(25)
1 4 11 2 3 < μ < 1 4 11 + 2 3 , 0 m < 2 μ 15 μ 2 5 μ + 2 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 , or
(26)
0 μ < 1 4 11 2 3 , m = 64 35 8 μ ( 8 μ + 1 ) 2 μ 8 35 8 μ ( 8 μ + 1 ) 2 + 48 8 μ + 1 , or
(27)
0 μ < 1 4 11 2 3 , m = 64 35 8 μ ( 8 μ + 1 ) 2 μ 8 35 8 μ ( 8 μ + 1 ) 2 48 8 μ 1 , or
(28)
0 μ < 1 2 5 17 , 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 < m < 48 8 μ + 1 8 8 μ 35 ( 8 μ + 1 ) 2 , or
(29)
0 μ < 1 2 5 17 , m > 8 8 μ 35 ( 8 μ + 1 ) 2 + 48 8 μ + 1 , or
(30)
1 4 11 + 2 3 < μ < 35 8 , m = 64 35 8 μ ( 8 μ + 1 ) 2 μ 8 35 8 μ ( 8 μ + 1 ) 2 + 48 8 μ + 1 , or
(31)
1 4 11 + 2 3 < μ < 35 8 , m = 64 35 8 μ ( 8 μ + 1 ) 2 μ 8 35 8 μ ( 8 μ + 1 ) 2 48 8 μ 1 , or
(32)
1 4 11 + 2 3 < μ 35 8 , 2 μ 15 μ 2 5 μ + 2 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 < m < 48 8 μ + 1 8 8 μ 35 ( 8 μ + 1 ) 2 , or
(33)
1 4 11 + 2 3 < μ < 35 8 , 8 8 μ 35 ( 8 μ + 1 ) 2 + 48 8 μ + 1 < m < 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 , or
(34)
1 4 11 2 3 < μ < 1 4 11 + 2 3 , m = 64 35 8 μ ( 8 μ + 1 ) 2 μ 8 35 8 μ ( 8 μ + 1 ) 2 + 48 8 μ + 1 , or
(35)
1 4 11 2 3 < μ < 1 4 11 + 2 3 , m = 64 35 8 μ ( 8 μ + 1 ) 2 μ 8 35 8 μ ( 8 μ + 1 ) 2 48 8 μ 1 , or
(36)
1 4 11 2 3 < μ < 1 4 11 + 2 3 , 2 μ 15 μ 2 5 μ + 2 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 < m < 48 8 μ + 1 8 8 μ 35 ( 8 μ + 1 ) 2 , or
(37)
1 4 11 2 3 < μ < 1 4 11 + 2 3 , 8 8 μ 35 ( 8 μ + 1 ) 2 + 48 8 μ + 1 < m < 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 , or
(38)
1 2 5 17 < μ < 1 4 11 2 3 , 2 μ 15 μ 2 5 μ + 2 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 < m < 48 8 μ + 1 8 8 μ 35 ( 8 μ + 1 ) 2 , or
(39)
1 2 5 17 < μ < 1 4 11 2 3 , 8 8 μ 35 ( 8 μ + 1 ) 2 + 48 8 μ + 1 < m < 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 , or
(40)
35 8 < μ < 1 2 5 + 17 ,
2 μ 15 μ 2 5 μ + 2 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 < m < 2 μ 15 μ 2 5 μ + 2 + 16 μ 2 120 μ + 249 μ 2 5 μ + 2 2 .
Recall that the flow in a neighborhood of antipodal points is topologically equivalent, and it may be reversed [98], but in this case, both points share the same eigenvalues. Therefore, the stability is the same for both B and C . We will write λ 1 ( B , C ) = f 1 ( μ , m ) and λ 2 ( B , C ) = f 2 ( μ , m ) to represent their eigenvalues. Given that many existence conditions and the eigenvalues depend on both free parameters, we will only consider some cases to analyze the stability of B and C . In particular, if we only consider the existence conditions for which the points are hyperbolic (most), we verify by numerical inspection that the points can never be attractors. However, they can be sources or saddles. When one parameter is fixed and the other free, the behavior is as depicted in Figure 6. On the other hand, when both parameters are free, the behaviour is as shown in Figure 7.
3.
D , E = ± k 2 m k 1 + ( 2 μ 3 ) m 4 k 2 m ( 4 μ + ( ( μ 5 ) μ + 2 ) m + 30 ) 24 , k 2 k , where k and k 1 are defined the same as in (77) and (78) but
k 2 = m 3 k 1 ( 2 μ 3 ) 8 μ ( μ ( 2 μ 25 ) + 80 ) + 204 4 m 2 k 1 + 4 μ ( μ + 15 ) 436 + 2 ( μ 1 ) μ ( ( μ 9 ) μ + 21 ) + 13 m 4 + 96 ( 2 μ 15 ) m + 288 .
As before, D and E are antipodal points, and once again, they share the same eigenvalues. We will write them as γ 1 = g 1 ( μ , m ) and γ 2 = g 2 ( μ , m ) . The conditions for D and E are the same as the forty cases. However, D and E can be attractors or saddles in those intervals. Additionally, numerical inspection in every interval shows that the points cannot be sources. See Figure 8 and Figure 9.
From Figure 10, we have that B, C, D, and E do not exist for the best-fit values of μ and m according to the analysis in Section 5.
Finally, we obtain the dynamics shown in Figure 11. Here, we observe that for ( μ = 1 , m = 3 ) or ( μ = 1 , m = 4 ) or ( μ = 2 , m = 3 ) , point A is an attractor. However, for ( μ = 2 , m = 4 ) , there are two complex eigenvalues with zero real parts, namely i 3 2 , i 3 2 which means A is a center. Finally, for ( μ = 3 , m = 3 ) , the point A is a source, B and C are saddles and D and E are sinks. Finally for ( μ = 3 , m = 4 ) , A is a saddle point, B and C are sources and D and E are sinks.

5. Observational Constraints

In this section, we constrain the free parameters of the exact scaling solution obtained for the exponential potential (see Section 3). To achieve this, we compute the best-fit parameters at the 1 σ confidence level (CL) using the affine-invariant Markov chain Monte Carlo (MCMC) method [96], implemented in the pure Python code emcee version 3.1.6 [97]. This analysis incorporates data from SNe Ia, CC, GL, and BHS.
In our procedure, we consider 100 chains, or “walkers”, and use the autocorrelation time τ corr provided by the emcee module as a convergence test. At every 50 steps, we compute the value of τ corr for each free parameter. If the current step exceeds 50 τ corr and the value of τ corr changes by less than 1 % , we conclude that the chains have converged and stop the constraints. We discard the first 5 τ corr steps as “burn-in” steps, thin the chains by τ corr / 2 , and flatten them. For this MCMC analysis, we consider the following Gaussian likelihood:
L I exp χ I 2 2 ,
where χ I 2 is the merit function. In particular, we consider the combinations of data sets χ SNe 2 , χ SNe + CC 2 = χ SNe 2 + χ CC 2 , and χ joint 2 = χ SNe 2 + χ CC 2 + χ GL 2 + χ BHS 2 . As initialization points, we consider a small vicinity around the values that maximize the likelihood.
We begin by constructing the Hubble parameters that describe the exact scaling solution in the following subsections. Next, we provide a brief overview of the merit function for each data set included in the cosmological constraints. Finally, we present the results along with a discussion of their implications.

5.1. Hubble Parameter for the Exponential Potential

To derive the theoretical Hubble parameter for the scaling solution associated with the exponential potential, we substitute Equations (30), (46) and (47) into Equation (24). This yields
H ˙ = 4 12 c 4 λ 2 t 3 H 3 + 1 3 λ 2 t 2 8 c 4 t 2 H 2 1 ( μ 3 ) H t H 2 ,
where
λ = 1 3 36 3 m 4 μ + ( ( μ 5 ) μ + 2 ) m + 30 ( μ + 2 ) m 6 , c 4 = 1 96 m 2 2 m 3 λ 2 + m λ 2 ( μ 6 ) m + 2 + 3 .
Hence, by introducing the logarithmic independent variable s = ln ( 1 + z ) , for which s as z , s 0 as z 0 , s as z 1 , and defining the age parameter as α = t H , we obtain the initial value problem
α ( s ) = 2 16 c 4 α ( s ) + 3 λ 2 3 λ 2 8 c 4 α ( s ) 2 1 μ 4 3 λ 2 α ( s ) α ( s ) + 6 ,
t ( s ) = t ( s ) α ( s ) ,
α ( 0 ) : = α 0 = t 0 H 0 , t ( 0 ) = t 0 .
The above system has the exact solution
α ( s ) = 2 m , m 6 μ ,
where α ( 0 ) : = α 0 = 2 m = t 0 H 0 implies m = 2 α 0 . So, integrating under the initial condition t ( 0 ) = t 0 , we obtain t ( s ) = t 0 e s α 0 , and, therefore
H ( s ) = α ( s ) t ( s ) 1 = H 0 e s α 0 H ( z ) = H 0 ( 1 + z ) 1 α 0 .
The deceleration parameter is given by q ( z ) = 1 + α 0 1 . Therefore, the Universe is accelerated without invoking dark energy for α 0 > 1 .
Finally, from Equation (25), we have the effective matter sources:
Ω ϕ : = 1 2 ϕ ˙ 2 + V ( ϕ ) 3 H 2 = 1 2 α 0 ( 3 α 0 μ 6 ) + 4 ( μ 1 ) α 0 2 α 0 μ + 3 α 0 ( α 0 5 ) μ 2 + 5 μ 2 ,
and
Ω X : = Ω fracc + Ω GB + Ω GB , fracc = 3 ( μ 1 ) H 8 H 2 ψ 1 t + 24 H 3 ψ ˙ 3 H 2 = 2 α 0 ( 3 α 0 μ 6 ) + 4 ( μ 1 ) α 0 2 μ α 0 + 3 α 0 ( α 0 5 ) μ 2 + 5 μ 2 ,
with
Ω ϕ + Ω X = 1 ,
where, for this scaling solution, the scalar field behaves as dark matter, say
Ω ϕ = Ω D M ,
and the extra terms Ω fracc + Ω GB + Ω GB , fracc mimic dark energy Ω X .
To compare with the fractional cosmological model, we also constrain the Λ CDM model, for which the corresponding theoretical Hubble parameter is given by
H = H 0 Ω m , 0 ( 1 + z ) 3 + 1 Ω m , 0 ,
whose corresponding parameter space is θ = ( H 0 , Ω m , 0 ) .
For this scaling solution (86), we have that Ω ϕ and Ω X are constants. Therefore,
1 2 α 0 ( 3 α 0 μ 6 ) + 4 ( μ 1 ) α 0 2 α 0 μ + 3 α 0 ( α 0 5 ) μ 2 + 5 μ 2 = Ω m , 0 .
Hence, we have two solutions
μ ± = 1 2 α 0 ( Ω m , 0 1 ) { α 0 ( 2 α 0 ( Ω m , 0 1 ) + 5 Ω m , 0 7 ) + 4 ± α 0 [ α 0 ( 16 α 0 2 ( Ω m , 0 1 ) 2 8 α 0 ( 5 Ω m , 0 7 ) ( Ω m , 0 1 ) + Ω m , 0 ( 17 Ω m , 0 86 ) + 73 ) + 8 ( 3 Ω m , 0 5 ) ] + 16 } .
Therefore, we use a Bayesian analysis to obtain the best-fit values of Ω m , 0 and α 0 and substitute in the previous expressions to obtain the physical values of μ .

5.2. Cosmic Chronometers

To constrain the model using cosmic chronometers, we utilize the data set from Ref. [99], which comprises 31 data points spanning the redshift range 0.0708 z 1.965 . These Hubble data points are derived using the differential age method, a model-independent approach [100]. Accordingly, the merit function for the CC data is formulated as
χ C C 2 = i = 1 31 H i H t h ( z i , θ ) σ H , i 2 ,
where H i represents the observational Hubble parameter at redshift z i , with an associated error σ H , i provided by the cosmic chronometer (CC) sample. The term H th denotes the theoretical Hubble parameter at the same redshift, while θ encompasses the free parameters of the model.
Note that the same theoretical Hubble parameter gives both analytical solutions of interest for the cosmological constraint, according to Equation (87), whose parameter space is θ = ( H 0 , α 0 ) . Therefore, we consider for our MCMC analysis the flat priors 0.55 < h < 0.85 and 0.5 < α 0 < 2.5 , where h is the reduced Hubble constant according to the expression H 0 = 100 k m / s M p c h . We also constrain the Λ CDM model as a further comparison, whose respective theoretical Hubble parameter is given by (92), whose respective parameter space is θ = ( H 0 , Ω m , 0 ) , for which we consider the same prior on h as in the fractional cosmology, plus the flat prior 0 < Ω m , 0 < 1 .

5.3. Type Ia Supernovae

For the SNe Ia data, we consider the Pantheon+ sample [101], which consists of 1701 data points in the redshift range 0.001 z 2.26 , whose respective merit function can be conveniently constructed in matrix notation (denoted by bold symbols) as
χ SNe 2 = D ( z , θ , M ) C 1 D ( z , θ , M ) ,
where D ( z , θ , M ) i = m B , i M μ t h ( z i , θ ) and C = C stat + C sys , with C the total uncertainty covariance matrix. The matrices C stat and C sys account for the statistical and systematic uncertainties, respectively. The quantity μ i = m B , i M corresponds to the observational distance modulus of the Pantheon+ sample, which is obtained by a modified version of Trip’s formula [102] and the BBC (BEAMS with Bias Corrections) approach [103]. In contrast, m B , i is the corrected apparent B-band magnitude of a fiducial SNe Ia at redshift z i , and M is the fiducial magnitude of an SNe Ia, which must be jointly estimated with the free parameters of the model under study.
The theoretical distance modulus for a spatially flat FLRW spacetime is given by
μ t h ( z i , θ ) = 5 log 10 d L ( z i , θ ) Mpc + 25 ,
with d L ( z i , θ ) as the luminosity distance given by
d L ( z i , θ ) = c ( 1 + z i ) 0 z i d z H t h ( z , θ ) ,
where c is the speed of light given in units of km/s.
In principle, there is a degeneration between M and H 0 . Hence, to constrain H 0 using SNe Ia data alone, it is necessary to include the SH0ES (Supernovae and H 0 for Equation of State of the dark energy program) Cepheid host distance anchors, with a merit function of the form
χ Cepheid 2 = D Cepheid M C 1 D Cepheid M ,
where D Cepheid M i = μ i M μ i Cepheid , where μ i Cepheid is the Cepheid calibrated host galaxy distance obtained by SH0ES [104]. So, we use the Cepheid distances as the “theory model” to calibrate M, considering that the difference μ i M μ i Cepheid is sensitive to M and largely insensitive to other parameters of the cosmological model. Taking into account that the total uncertainty covariance matrix for Cepheid is contained in the total uncertainty covariance matrix C , we define the merit function for the SNe Ia data as
χ SNe 2 = D ( z , θ , M ) C 1 D ( z , θ , M ) ,
where
D i = m B , i M μ i Cepheid i Cepheid host m B , i M μ t h ( z i , θ ) otherwise .
For the nuisance parameter M, we consider the flat prior 20 < M < 18 in our MCMC analysis.

5.4. Gravitational Lensing

Multiple images are produced when a background object (the source) is lensed due to the gravitational force of a massive body (the lens). Therefore, the light rays emitted from the source will take different paths through spacetime at different image positions and arrive at the observer at different times. In this sense, the time delay of two different images k and l depends on the mass distribution along the line of sight of the lensing object, which can be calculated as follows:
t k l = D t c ϕ k β 2 2 ψ ( ϕ k ) ϕ l β 2 2 + ψ ( ϕ l ) ,
where ϕ k and ϕ l are the angular position of the images, β is the angular position of the source, ψ ( ϕ k ) and ψ ( ϕ l ) are the lens potential at the image positions, and D t is the “time-delay distance”, which is theoretically given by the expression [105]
D t t h ( z , θ ) = 1 + z l d A , l ( z l , θ ) d A , s ( z s , θ ) d A , l s ( z l s , θ ) ,
where the subscripts l, s, and l s stand for the lens, the source, and between the lens and the source, respectively; z = ( z l , z s , z l s ) ; and d A , j is the angular diameter distance, which can be written in terms of the luminosity distance (98) as d L ( z j , θ ) = d A , j ( 1 + z j ) 2 or
d A , j ( z j , θ ) = c ( 1 + z j ) 0 z j d z H t h ( z , θ ) .
In this paper, we consider the gravitational lensing compilation provided by the H0LiCOW collaboration [106], which consists of six lensed quasars: B1608+656 [107], SDSS 1206+4332 [108], WFI2033-4723 [109], RXJ1131-1231, HE 0435-1223, and PG 1115-080 [110], whose respective merit function can be written as
χ GL 2 = i = 1 6 D t , i D t t h ( z i , θ ) σ D t , i 2 ,
where D t , i is the observational time-delay distance of the lensed quasar at redshift z i = ( z l , i ; z s , i ; z l s , i ) with an associated error σ D t , i (for more details, see Ref. [106]). It is important to note that, for z 0 , the angular diameter distance (104) tends to d A c z / H 0 and, therefore, the gravitational lensing data of the H0LiCOW collaboration is sensitive to H 0 , with a weak dependency on other cosmological parameters.

5.5. Black Hole Shadows

Black Hole Shadow (BHS) data are of interest to study our local Universe since their dynamic is quite simple and can be seen as standard rulers if the angular size redshift α , the relation between the size of the shadow, and the mass of the supermassive black hole that produces it are established [111]. In this paper, we are interested in two measures: the first one was made on the M87* supermassive black hole by The Event Horizon Telescope Collaboration [112] (the first detection of a BHS) and the second one corresponds to the detection of Sagittarius A* (Sgr A*) [113].
Light rays curve around its event horizon in a black hole (BH), creating a ring with a black spot center, the so-called shadow of the BH. The angular radius of the BHS for a Schwarzschild (SH) BH at redshift z i is given by
α S H z i , θ = 3 3 m d A ( z i , θ ) ,
where d A ( z i , θ ) is given by Equation (104) (note that the sub-index j is not necessary in this case) and m = G M B H / c 2 is the mass parameter of the BH, with M B H as the mass of the BH in solar masses units and G as the gravitational constant.
It is common to write Equation (106) in terms of the shadow radius α S H ( z i , θ ) = R S H / d A ( z i , θ ) , where R S H = 3 3 G M B H / c 2 (the speed of light is given in units of m/s in this case). Therefore, the merit function for the BHS data can be constructed as
χ B H S 2 = i = 1 2 α i α S H ( z i , θ ) σ α , i 2 ,
where α i is the observational angular radius of the BHS at redshift z i with an associated error σ α , i . It is important to note that for z 0 , the angular radius (106) tends to α S H R S H H 0 / c z , and, therefore, as well as in the gravitational lensing data, the BHS data are sensitive to H 0 , with a weak dependency on other cosmological parameters. On the other hand, we divide Equation (106) by a factor of 1.496 × 10 11 to obtain α S H in units of μ a s .

5.6. Results and Discussion

In Table 1, we present the total number of steps and the correlation time for the free parameter space of the Λ CDM and the fractional cosmology. In Table 2, we present their respective best-fit values at the 1 σ CL and χ min 2 criteria. In Figure 12 and Figure 13, we depict the posterior 1D distribution and joint marginalized regions of the free parameter space of the Λ CDM and fractional cosmologies at 1 σ ( 68.3 % ) , 2 σ ( 95.5 % ) , and 3 σ ( 99.7 % ) CL, respectively. These results were obtained by the MCMC analysis described in Section 5 for the SNe Ia, SNe Ia+CC, and SNe Ia+CC+GL+BHS (joint) data.
In Table 2, we can see that the Λ CDM model exhibits lower values of the χ min 2 criteria compared to the scaling solution explored in this paper. However, it is essential to note that this solution does not fully represent the behavior of the fractional cosmology studied in this paper. It only provides an approximate solution to the complete picture. Therefore, the constraint presented in the table offers insights into the ability of this fractional cosmology to describe the observed Universe. An explanation of this behavior can be seen in Figure 14, where we depict the Hubble parameter for the Λ CDM cosmology and the fractional cosmology for the exponential potential as a function of the redshift z. The figure shows that the solution can only mimic the Λ CDM model at a redshift z < 0.5 . The solutions are moving away from this point, with the differences being greater when the redshift increases.
To reconstruct the other parameters that characterized the solution, we use the chains obtained in our MCMC analysis described in Section 5 for the fractional cosmology in the case of the joint analysis. Following this line, considering that α = 2 / m , we obtain at 1 σ CL the value m = 1.44 ± 0.05 . Also, considering that α 0 = t 0 H 0 , we obtain at 1 σ CL the value t 0 = 19.0 ± 0.7 [Gyr]. On the other hand, to infer the value of μ for the scaling solution obtained for the exponential potential, we also consider the chains obtained in our MCMC analysis described in Section 5 for the Λ CDM model in the case of the joint analysis. Therefore, from Equation (88), we obtain the approximated values μ = 1.48 ± 0.17 and μ = 4.974 ± 0.10 .
Focusing on the observational suitable solution for the joint analysis, in Figure 15, we depict the deceleration parameter for the fractional cosmology obtained for the exponential potential as a function of the redshift z, with an error band at 1 σ CL. We also depict the deceleration parameter for the Λ CDM cosmology as a reference model. From this figure, we can conclude that this solution only represents an always expanding solution and does not exhibit a transition between a decelerated solution and an accelerated one such as the Λ CDM model. The deceleration parameter for this fractional solution is constant and given by q = 1 + α 0 1 , with a value at 1 σ CL of q 0 = 0.28 ± 0.03 at the current time; i.e., this solution also represents a less accelerated solution than the Λ CDM cosmology at the current time. Nevertheless, we must again emphasize that this solution is only a particular solution and does not represent the complete picture of fractional cosmology. This solution sheds some light on the capability of fractional cosmology in describing the cosmological background, where, as we can see, we can obtain an accelerated solution at the current time without invoking any dark energy.
Ultimately, we construct the H 0 ( z ) diagnostic [114] for the observational suitable fractional solution, which can give us insights into the possibility of alleviating the H 0 tension. Following this line, the H 0 ( z ) diagnostic is defined by
H 0 ( z ) = H ( z ) / Ω m , 0 ( 1 + z ) 3 + 1 Ω m , 0 ,
where the Hubble parameter is obtained numerically by H ( z ) = H t h ( z ) . So, in Figure 16, we illustrate the H 0 diagnostic for the Λ CDM cosmology and the fractional solution as a function of the redshift z, with an error band at 1 σ CL. This figure shows that at redshift z > 0.5 (the same redshift where both models are moving away according to Figure 14), the value of H 0 for the fractional cosmology presents a run to values lower than H 0 = 74.03 ± 1.42 k m / s M p c , obtained by model-independent measurements of cepheid [115]. Therefore, fractional cosmology can be a suitable framework to alleviate the H 0 tension.

6. Concluding Remarks

We have developed a new theoretical framework called the Fractional Einstein–Gauss–Bonnet scalar field cosmology, which has important physical implications. We used fractional calculus to modify the gravitational force with an arbitrary constant fractional-order derivative. This allowed us to derive a modified Friedmann equation and a modified Klein–Gordon equation using the Fractional Action-like Variational Approach.
One critical issue is whether the obtained numerical results will be the same if the model is formulated with another fractional-order derivative. Caputo’s definition of fractional and integral derivatives has several advantages compared to Riemann–Louville or Grünwald–Letnikov. First, it considers the values of the function and its derivatives at the origin (or, in general, at any lower or upper bound a), making it suitable for solving fractional-order initial value problems using Laplace transforms. Additionally, the Caputo fractional derivative of a constant is 0, making it more consistent with classical calculation. However, after applying different definitions of the fractional derivative to the same function, we found that they give different results, but all these definitions are equally well–founded. Physicists can test different definitions of the fractional derivative by comparing theoretical results with experimental data. Depending on the specific problem, one of the definitions will show the most significant agreement with the experiment.
In our research, we utilized fractional calculus to derive modified Friedmann and Klein–Gordon equations, resulting in changes to the gravitational action integral described in Section 2. We discovered non-trivial solutions related to exponential potential, exponential couplings to the GB term, and a logarithmic scalar field. Our findings show that we achieved late-time accelerating power-law solutions for the scale factor and Hubble parameter with inverse time power-law expansion, as discussed in Section 3. We also analyzed the phase space structure using linear stability theory and examined the dynamical effects of the Gauss–Bonnet couplings. The stability analysis of the solutions, which is a critical aspect of our research, was conducted in Section 4, where we presented a reconstruction procedure for the potential and the coupling functions. To reconstruct the other parameters characterizing the solution, we utilized the chains obtained in our MCMC analysis described in Section 5 for the fractional cosmology in the case of the joint analysis. Additionally, by considering that α = 2 / m , we derived specific approximate values for cosmological quantities, assuming they are definable, as the fractional modification of the derivative concept can yield results that closely align with observational cosmological measurements. Notably, the parameters m and α 0 were estimated using data from cosmic chronometers, type Ia supernovae, supermassive Black Hole Shadows, and strong gravitational lensing, indicating that the solution is consistent with an accelerated expansion ( q ( z ) < 0 ) at late times with the values α 0 = 1.38 ± 0.05 , and m = 1.44 ± 0.05 . The Universe’s age in the literature was reported as 26.7 [Gyr] for hybrid covarying coupling constants together with tired light models in [116] using Pantheon + data fit parameters with 5.8 [Gyr] at z = 10 and 3.5 [Gyr] at z = 20 , giving enough time to form massive galaxies. However, in previous works within the fractional cosmology context, the authors reported other Universe ages that are larger than expected under the standard paradigm; for example, 33 . 633 15.095 + 14.745 [Gyr] (CC), 33 . 837 10.788 + 27.833 [Gyr] (SNe Ia) and 33 . 617 4.511 + 3.411 [Gyr] (joint) were reported in [74]. Moreover, the value 25 . 62 4.46 + 6.89 [Gyr] (SNe Ia+OHD data) at 3 σ CL was obtained in [75]. Our research reports a value of t 0 = 19.0 ± 0.7 [Gyr] (SNe Ia+CC+GL+BHS) at 1 σ CL. Still, it could be larger than the age of the Universe expected under the standard paradigm. However, these values do not contradict the minimum bound expected for the Universe’s age imposed by globular clusters t 0 = 13 . 5 0.14 + 0.16 ± 0.23 [117], and, as far as we know, the maximum bound does not exist and is model-dependent.
We then used Bayesian analysis to determine the value for μ for the scaling solution found for the exponential potential; we also examined the chains generated in our MCMC analysis outlined in Section 5 for the Λ CDM model when performing joint analysis. Thus, using Equation (88), we derived the approximate values of μ = 1.48 ± 0.17 and μ = 4.974 ± 0.10 (consistent with Ω m , 0 = 0.311 ± 0.016 and h = 0.712 ± 0.007 ). We focused on finding an appropriate observational solution for joint analysis (with μ = 1.48 ± 0.17 ). Hence, we depicted the deceleration parameter for fractional cosmology with an exponential potential as a function of the redshift (z). We also included an error band at a confidence level of 1 σ . We used the Λ CDM cosmology as a reference model. This solution shows a constant deceleration parameter and does not transition from a decelerated phase to an accelerated one like the Λ CDM model. The deceleration parameter for this fractional solution is constant and given by q = 1 + α 0 1 , with a value at a confidence level of 1 σ of q 0 = 0.28 ± 0.03 at the current time. This means this solution represents a less accelerated expansion than the current Λ CDM cosmology. However, it is essential to note that this solution is only a specific case and does not represent the complete picture of fractional cosmology. This solution provides insight into the potential of fractional cosmology in describing the cosmological background, showing that we can achieve an accelerated expansion at the current time without invoking dark energy.
We have constructed the H 0 ( z ) diagnostic [114] to analyze the fractional solution for observational suitability. This diagnostic provided insights into the potential to alleviate the H 0 tension. Following this approach, the H 0 ( z ) diagnostic (108) involves obtaining the Hubble parameter numerically as H ( z ) = H t h ( z ) . We have plotted the H 0 diagnostic for both the Λ CDM cosmology and the fractional solution as a function of the redshift z, along with a 1 σ CL error band. Our analysis reveals that for redshifts z > 0.5 , the value of H 0 for the fractional cosmology consistently trended lower than H 0 = 74.03 ± 1.42 k m / s M p c , which was obtained from model-independent measurements of cepheids [115]. Therefore, fractional cosmology shows promise as a framework for addressing the H 0 tension.
Our model identifies various equilibrium points corresponding to cosmological scenarios, including accelerated and scaling solutions. The scaling behavior at specific equilibrium points suggests that the geometric corrections in the coupling to the GB scalar can replicate the behavior of the dark sector in modified gravity. Moreover, attracting scaling solutions solves the Coincidence Problem.
This cosmological model, derived from fractional calculus, is consistent with recent cosmological observations and provides an alternative explanation for cosmic acceleration, distinct from the Λ CDM model. Our results build upon and significantly enhance previous work in the field, emphasizing the practical implications of fractional calculus in cosmology. Additionally, our model shows promise in addressing the persistent tension in the Hubble constant, offering hope for future research [76].

Author Contributions

Conceptualization, B.M.-R., B.D., A.D.M. and G.L.; methodology, G.L.; software, B.M.-R., B.D., A.D.M., G.L., E.G. and J.M.; validation, B.M.-R., B.D., A.D.M., G.L., E.G. and J.M.; formal analysis, B.M.-R., B.D., A.D.M., G.L., E.G. and J.M.; investigation, B.M.-R., B.D., A.D.M., G.L., E.G. and J.M.; resources, G.L. and J.M.; writing—original draft preparation, G.L.; writing—review and editing, G.L.; visualization, B.M.-R., B.D., A.D.M., G.L., E.G. and J.M.; supervision, G.L., B.D. and A.D.M.; project administration, G.L., J.M., B.M.-R. and B.D.; funding acquisition, B.M.-R., B.D., A.D.M., G.L., E.G. and J.M. All authors have read and agreed to the published version of the manuscript.

Funding

B.M.-R., G.L., J.M., and A.D.M. acknowledge the financial support of Agencia Nacional de Investigación y Desarrollo (ANID), Chile, through Proyecto Fondecyt Regular 2024, Folio 1240514, Etapa 2024. G.L. was also funded by Vicerrectoría de Investigación y Desarrollo Tecnológico (VRIDT) at Universidad Católica del Norte (UCN) through Resolución VRIDT N°026/2023, Resolución VRIDT N°027/2023, Proyecto de Investigación Pro Fondecyt 2023 (Resolución VRIDT N°076/2023) and Resolución VRIDT N°09/2024. A.D.M. was also supported by Agencia Nacional de Investigación y Desarrollo - ANID Subdirección de Capital Humano/Doctorado Nacional/año 2020 folio 21200837, Gastos operacionales proyecto de Tesis/2022 folio 242220121 and VRIDT-UCN. E.G. was funded by Vicerrectoría de Investigación y Desarrollo Tecnológico (VRIDT) at Universidad Católica del Norte (UCN) through Proyecto de Investigación Pro Fondecyt 2023, Resolución VRIDT N°076/2023.

Data Availability Statement

The data supporting this article can be found in Section 5.

Acknowledgments

G.L. is thankful for the support of Núcleo de Investigación Geometría Diferencial y Aplicaciones, Resolución VRIDT N°096/2022. E.G. acknowledges the scientific support of Núcleo de Investigación No. 7 UCN-VRIDT 076/2020, Núcleo de Modelación y Simulación Científica (NMSC).

Conflicts of Interest

We declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A. Variational Equations

We start with an action defined as [95]
S = 1 Γ ( μ ) 0 τ L θ , q i ( θ ) , q ˙ i ( θ ) , q ¨ i ( θ ) ( τ θ ) μ 1 d θ ,
where q ˙ i ( θ ) , q ¨ i ( θ ) means integer order derivatives with respect to the argument θ .
To obtain the equations of motion of the fields, we do the infinitesimal transformation of the fields
q i q i + δ q i , δ L = L ( q i + δ q i ) L ( q i ) .
Then, we write the variation of the action:
δ S = 1 Γ ( μ ) 0 τ ( δ L ) ( τ θ ) μ 1 + L δ ( ( τ θ ) μ 1 ) d θ = 1 Γ ( μ ) 0 τ L q i δ q i + L q ˙ i δ q ˙ i + L q ¨ i δ q ¨ i ( τ θ ) μ 1 d θ .
We can integrate by parts:
0 τ L q ˙ i ( τ θ ) μ 1 δ q ˙ i d θ = L q ˙ i ( τ θ ) μ 1 δ q i | 0 τ 0 τ d d θ L q ˙ i ( τ θ ) μ 1 δ q i d θ ,
0 τ L q ¨ i ( τ θ ) μ 1 δ q ¨ i d θ = L q ¨ i ( τ θ ) μ 1 δ q ˙ i | 0 τ 0 τ d d θ L q ¨ i ( τ θ ) μ 1 δ q ˙ i d θ = d d θ L q ¨ ( τ θ ) μ 1 δ q i | 0 τ + L q ¨ i ( τ θ ) μ 1 δ q ˙ i | 0 τ + 0 τ d 2 d θ 2 L q ¨ i ( τ θ ) μ 1 δ q i d θ .
Therefore, the variation of the action can be written as
δ S = 1 Γ ( μ ) 0 τ L q i ( τ θ ) μ 1 d d θ L q ˙ i ( τ θ ) μ 1 + d 2 d θ 2 L q ¨ i ( τ θ ) μ 1 δ q i d θ + L q ˙ i ( τ θ ) μ 1 d d θ L q ¨ ( τ θ ) μ 1 δ q i | 0 τ + L q ¨ i ( τ θ ) μ 1 δ q ˙ i | 0 τ .
Then, assuming that δ q i ( 0 ) = 0 at the lower end, we have the equations of motion:
L q i ( τ θ ) μ 1 d d θ L q ˙ i ( τ θ ) μ 1 + d 2 d θ 2 L q ¨ i ( τ θ ) μ 1 = 0 .
For the more general case where the Lagrangian depends on derivatives of up to order n > 2 of the generalized coordinates q i , the same procedure can be used to show that, through successive integration by parts, alternating signs emerge, and the equations of motion are
k = 0 n ( 1 ) k d k d θ k L q i ( k ) τ θ μ 1 = 0 .

Appendix B. Stability Analysis of Power-Law Solutions

In general, when we have an ordinary differential equation
F t , ψ t , ψ ˙ t , ψ ¨ t , . . . 0 ,
where t is the independent variable and ψ t is the dependent variable; we can provide the analysis of the stability of the solution ψ s ( t ) in the interval 0 < t < using similar methods as in [92,93,94]. Defining the new time variable
t = e τ , < τ < ,
such that t 0 as τ and t as τ , as well as the ratio
u ( τ ) = ψ ( τ ) ψ s ( τ ) ,
where ψ ( τ ) = ψ ( e τ ) and ψ s ( τ ) = ψ s ( e τ ) are given. Notice that, evaluating at ψ s ( τ ) , we have u = 1 ; therefore, defining ε = ψ ( τ ) ψ s ( τ ) 1 , the solution is shifted to ε = 0 . Let
ψ d ψ d τ ,
then,
ψ ˙ = d τ d t ψ = e τ ψ , ψ ¨ = e 2 τ ψ ψ .
Therefore, Equation (A9) becomes
G τ , u τ , u τ , u τ , 0 .
According to (59), ψ ( t ) satisfies
ψ ¨ = c 1 + c 2 ψ t 2 + c 3 ψ ˙ t ,
where c 1 , c 2 , and c 3 are the constants. Additionally, given the solution found for the exponential potential, (46), the critical solution (63) can be expressed as
ψ c ( t ) = c 4 t 2 ,
where c 4 is also a constant.
Thus, replacing ψ ¨ , ψ ˙ , and t, (A15) becomes
e 2 τ ψ ( τ ) ψ ( τ ) = c 1 + c 2 e 2 τ ψ ( τ ) + c 3 e 2 τ ψ ( τ ) ,
That is,
ψ ( τ ) = c 1 e 2 τ + c 2 ψ ( τ ) + c 3 + 1 ψ ( τ ) .
Using the definition of ε , we see that
ψ ( τ ) = ( ε ( τ ) + 1 ) ψ c ( τ ) , ψ ( τ ) = ε ( τ ) ψ c ( τ ) + ( ε ( τ ) + 1 ) ψ c ( τ ) , ψ ( τ ) = ε ( τ ) ψ c ( τ ) + 2 ε ( τ ) ψ c ( τ ) + ( ε ( τ ) + 1 ) ψ c ( τ ) .
Therefore,
ε ψ c + 2 ε ψ c + ( ε + 1 ) ψ c = c 1 e 2 τ + c 2 ( ε + 1 ) ψ c + ( c 3 + 1 ) ε ψ c + ( ε + 1 ) ψ c .
Equation (A20) can be written as
ε ψ c + 2 ε ψ c + ε ψ c + ψ c c 1 e 2 τ c 2 ψ c ( c 3 + 1 ) ψ c = 0 , ψ c ( τ ) satisfies ( A 18 ) . = c 2 ε ψ c + ( c 3 + 1 ) ε ψ c + ε ψ c
Due to ψ c ( τ ) being a non-trivial particular solution of Equation (A18), the term inside the square brackets on the left-hand side of (A21) is zero. After simplification, we acquire
ε + 2 ε ψ c ψ c + ε ψ c ψ c = c 2 ε + ( c 3 + 1 ) ε + ε ψ c ψ c .
Note that
ψ c ( τ ) = c 4 e 2 τ , ψ c ( τ ) = 2 c 4 e 2 τ , ψ c ( τ ) = 4 c 4 e 2 τ .
Therefore,
ψ c ψ c = 2 , ψ c ψ c = 4 .
Hence, by reducing terms, we have
ε = ε ( c 3 3 ) + ε ( c 2 + 2 c 3 2 ) .
Defining v = ε , we have the linear system
ε = v , v = ε ( c 2 + 2 c 3 2 ) + v ( c 3 3 ) ,
that was investigated using dynamical systems tools in Section 4 for the choice of the constants:
c 1 = 1 32 m ( 4 μ + ( μ 2 ) ( μ 1 ) m + 18 ) 12 , c 2 = m μ ( μ m + m + 4 ) 22 + 12 m 2 , c 3 = 2 μ 2 m , c 4 = 1 96 m 2 2 m 3 λ 2 + m λ 2 ( μ 6 ) m + 2 + 3 .

References

  1. Riess, A.G.; Filippenko, A.V.; Challis, P.; Clocchiatti, A.; Diercks, A.; Garnavich, P.M.; Gillil, 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]
  2. Aghanim, N.; Akrami, Y.; Ashdown, M.; Aumont, J.; Baccigalupi, C.; Ballardini, M.; Banday, A.J.; Barreiro, R.B.; Bartolo, N.; Basak, S.; et al. Planck 2018 results. VI. Cosmological parameters. Astron. Astrophys. 2020, 641, A6, Erratum in Astron. Astrophys. 2021, 652, C4. [Google Scholar] [CrossRef]
  3. Zeldovich, Y.B. The cosmological constant and the theory of elementary particles. Sov. Phys. Uspekhi 1968, 11, 381. [Google Scholar] [CrossRef]
  4. Weinberg, S. The cosmological constant problem. Rev. Mod. Phys. 1989, 61, 1–23. [Google Scholar] [CrossRef]
  5. Carroll, S.M. The Cosmological constant. Living Rev. Rel. 2001, 4, 1. [Google Scholar] [CrossRef]
  6. Rasouli, S.M.M.; Ziaie, A.H.; Marto, J.; Moniz, P.V. Gravitational Collapse of a Homogeneous Scalar Field in Deformed Phase Space. Phys. Rev. D 2014, 89, 044028. [Google Scholar] [CrossRef]
  7. Jalalzadeh, S.; Rasouli, S.M.M.; Moniz, P.V. Quantum cosmology, minimal length and holography. Phys. Rev. D 2014, 90, 023541. [Google Scholar] [CrossRef]
  8. Rasouli, S.M.M.; Vargas Moniz, P. Noncommutative minisuperspace, gravity-driven acceleration, and kinetic inflation. Phys. Rev. D 2014, 90, 083533. [Google Scholar] [CrossRef]
  9. Rasouli, S.M.M.; Vargas Moniz, P. Gravity-Driven Acceleration and Kinetic Inflation in Noncommutative Brans-Dicke Setting. Odessa Astron. Pub. 2016, 29, 19. [Google Scholar] [CrossRef]
  10. Jalalzadeh, S.; Capistrano, A.J.S.; Moniz, P.V. Quantum deformation of quantum cosmology: A framework to discuss the cosmological constant problem. Phys. Dark Univ. 2017, 18, 55–66. [Google Scholar] [CrossRef]
  11. Rasouli, S.M.M.; Ziaie, A.H.; Jalalzadeh, S.; Moniz, P.V. Non-singular Brans–Dicke collapse in deformed phase space. Ann. Phys. 2016, 375, 154–178. [Google Scholar] [CrossRef]
  12. Rasouli, S.M.M.; Farhoudi, M.; Vargas Moniz, P. Modified Brans–Dicke theory in arbitrary dimensions. Class. Quant. Grav. 2014, 31, 115002. [Google Scholar] [CrossRef]
  13. Melchiorri, A.; Mersini-Houghton, L.; Odman, C.J.; Trodden, M. The State of the dark energy equation of state. Phys. Rev. D 2003, 68, 043509. [Google Scholar] [CrossRef]
  14. Vikman, A. Can dark energy evolve to the phantom? Phys. Rev. D 2005, 71, 023515. [Google Scholar] [CrossRef]
  15. Nesseris, S.; Perivolaropoulos, L. Crossing the Phantom Divide: Theoretical Implications and Observational Status. J. Cosmol. Astropart. Phys. 2007, 2007, 018. [Google Scholar] [CrossRef]
  16. Elizalde, E.; Nojiri, S.; Odintsov, S.D. Late-time cosmology in (phantom) scalar-tensor theory: Dark energy and the cosmic speed-up. Phys. Rev. D 2004, 70, 043539. [Google Scholar] [CrossRef]
  17. Guo, Z.K.; Piao, Y.S.; Zhang, X.M.; Zhang, Y.Z. Cosmological evolution of a quintom model of dark energy. Phys. Lett. B 2005, 608, 177–182. [Google Scholar] [CrossRef]
  18. Zhang, X.F.; Li, H.; Piao, Y.S.; Zhang, X.M. Two-field models of dark energy with equation of state across-1. Mod. Phys. Lett. A 2006, 21, 231–242. [Google Scholar] [CrossRef]
  19. Feng, B.; Wang, X.L.; Zhang, X.M. Dark energy constraints from the cosmic age and supernova. Phys. Lett. B 2005, 607, 35–41. [Google Scholar] [CrossRef]
  20. Cai, Y.F.; Saridakis, E.N.; Setare, M.R.; Xia, J.Q. Quintom Cosmology: Theoretical implications and observations. Phys. Rept. 2010, 493, 1–60. [Google Scholar] [CrossRef]
  21. Wei, H.; Cai, R.G. Cosmological evolution of hessence dark energy and avoidance of big rip. Phys. Rev. D 2005, 72, 123507. [Google Scholar] [CrossRef]
  22. Wei, H.; Cai, R.G.; Zeng, D.F. Hessence: A New view of quintom dark energy. Class. Quant. Grav. 2005, 22, 3189–3202. [Google Scholar] [CrossRef]
  23. Wei, H.; Cai, R.G. A note on crossing the phantom divide in hybrid dark energy model. Phys. Lett. B 2006, 634, 9–13. [Google Scholar] [CrossRef]
  24. Zhao, G.B.; Xia, J.Q.; Li, M.; Feng, B.; Zhang, X. Perturbations of the quintom models of dark energy and the effects on observations. Phys. Rev. D 2005, 72, 123515. [Google Scholar] [CrossRef]
  25. Feng, B.; Li, M.; Piao, Y.S.; Zhang, X. Oscillating quintom and the recurrent Universe. Phys. Lett. B 2006, 634, 101–105. [Google Scholar] [CrossRef]
  26. Xia, J.Q.; Feng, B.; Zhang, X.M. Constraints on oscillating quintom from supernova, microwave background and galaxy clustering. Mod. Phys. Lett. A 2005, 20, 2409–2416. [Google Scholar] [CrossRef]
  27. Zhang, X. An interacting two-fluid scenario for quintom dark energy. Commun. Theor. Phys. 2005, 44, 762–768. [Google Scholar] [CrossRef]
  28. Zhang, X. Statefinder diagnostic for holographic dark energy model. Int. J. Mod. Phys. D 2005, 14, 1597–1606. [Google Scholar] [CrossRef]
  29. Zhang, X.; Wu, F.Q. Constraints on holographic dark energy from Type Ia supernova observations. Phys. Rev. D 2005, 72, 043524. [Google Scholar] [CrossRef]
  30. Zhang, X. Dynamical vacuum energy, holographic quintom, and the reconstruction of scalar-field dark energy. Phys. Rev. D 2006, 74, 103505. [Google Scholar] [CrossRef]
  31. Lazkoz, R.; Leon, G. Quintom cosmologies admitting either tracking or phantom attractors. Phys. Lett. B 2006, 638, 303–309. [Google Scholar] [CrossRef]
  32. Alimohammadi, M.; Sadjadi, H.M. The w = -1 crossing of the quintom model with arbitrary potential. Phys. Lett. B 2007, 648, 113–118. [Google Scholar] [CrossRef]
  33. Lazkoz, R.; Leon, G.; Quiros, I. Quintom cosmologies with arbitrary potentials. Phys. Lett. B 2007, 649, 103–110. [Google Scholar] [CrossRef]
  34. Leon, G.; Leyva, Y.; Socorro, J. Quintom phase-space: Beyond the exponential potential. Phys. Lett. B 2014, 732, 285–297. [Google Scholar] [CrossRef]
  35. Mohseni Sadjadi, H.; Alimohammadi, M. Transition from quintessence to phantom phase in quintom model. Phys. Rev. D 2006, 74, 043506. [Google Scholar] [CrossRef]
  36. Elizalde, E.; Nojiri, S.; Odintsov, S.D.; Saez-Gomez, D.; Faraoni, V. Reconstructing the Universe history, from inflation to acceleration, with phantom and canonical scalar fields. Phys. Rev. D 2008, 77, 106005. [Google Scholar] [CrossRef]
  37. Leon, G.; Paliathanasis, A.; Morales-Martínez, J.L. The past and future dynamics of quintom dark energy models. Eur. Phys. J. C 2018, 78, 753. [Google Scholar] [CrossRef]
  38. Tot, J.; Yildirim, B.; Coley, A.; Leon, G. The dynamics of scalar-field Quintom cosmological models. Phys. Dark Univ. 2023, 39, 101155. [Google Scholar] [CrossRef]
  39. Leon, G.; Coley, A.; Paliathanasis, A.; Tot, J.; Yildirim, B. Global dynamics of two models for Quintom Friedman–Lemaître– Robertson–Walker Universes. Phys. Dark Univ. 2024, 45, 101503. [Google Scholar] [CrossRef]
  40. Apostolopoulos, P.S.; Tetradis, N. Late acceleration and w = −1 crossing in induced gravity. Phys. Rev. D 2006, 74, 064021. [Google Scholar] [CrossRef]
  41. Bamba, K.; Nojiri, S.; Odintsov, S.D. Inflationary cosmology and the late-time accelerated expansion of the Universe in non-minimal Yang-Mills-F(R) gravity and non-minimal vector-F(R) gravity. Phys. Rev. D 2008, 77, 123532. [Google Scholar] [CrossRef]
  42. Bamba, K.; Geng, C.Q.; Nojiri, S.; Odintsov, S.D. Crossing of the phantom divide in modified gravity. Phys. Rev. D 2009, 79, 083014. [Google Scholar] [CrossRef]
  43. Setare, M.R.; Saridakis, E.N. Braneworld models with a non-minimally coupled phantom bulk field: A Simple way to obtain the -1-crossing at late times. J. Cosmol. Astropart. Phys. 2009, 2009, 002. [Google Scholar] [CrossRef]
  44. Nojiri, S.; Odintsov, S.D. Modified gravity with negative and positive powers of the curvature: Unification of the inflation and of the cosmic acceleration. Phys. Rev. D 2003, 68, 123512. [Google Scholar] [CrossRef]
  45. Nojiri, S.; Odintsov, S.D. Modified f(R) gravity consistent with realistic cosmology: From matter dominated epoch to dark energy Universe. Phys. Rev. D 2006, 74, 086005. [Google Scholar] [CrossRef]
  46. Cognola, G.; Elizalde, E.; Nojiri, S.; Odintsov, S.D.; Sebastiani, L.; Zerbini, S. A Class of viable modified f(R) gravities describing inflation and the onset of accelerated expansion. Phys. Rev. D 2008, 77, 046009. [Google Scholar] [CrossRef]
  47. Capozziello, S.; Nojiri, S.; Odintsov, S.D.; Troisi, A. Cosmological viability of f(R)-gravity as an ideal fluid and its compatibility with a matter dominated phase. Phys. Lett. B 2006, 639, 135–143. [Google Scholar] [CrossRef]
  48. Nojiri, S.; Odintsov, S.D. Unified cosmic history in modified gravity: From F(R) theory to Lorentz non-invariant models. Phys. Rep. 2011, 505, 59–144. [Google Scholar] [CrossRef]
  49. Nojiri, S.; Odintsov, S.D.; Oikonomou, V.K. Modified Gravity Theories on a Nutshell: Inflation, Bounce and Late-time Evolution. Phys. Rept. 2017, 692, 1–104. [Google Scholar] [CrossRef]
  50. Leon, G.; Saridakis, E.N. Phase-space analysis of Horava-Lifshitz cosmology. J. Cosmol. Astropart. Phys. 2009, 11, 006. [Google Scholar] [CrossRef]
  51. Bellorin, J.; Borquez, C.; Droguett, B. Effective action of the Hořava theory: Cancellation of divergences. Phys. Rev. D 2024, 109, 084007. [Google Scholar] [CrossRef]
  52. Bellorin, J.; Borquez, C.; Droguett, B. Renormalization of the nonprojectable Hořava theory. Phys. Rev. D 2024, 110, 084020. [Google Scholar] [CrossRef]
  53. Leon, G. On the Past Asymptotic Dynamics of Non-minimally Coupled Dark Energy. Class. Quant. Grav. 2009, 26, 035008. [Google Scholar] [CrossRef]
  54. Lim, S.C. Fractional derivative quantum fields at positive temperature. Physica A 2006, 363, 269–281. [Google Scholar] [CrossRef]
  55. Lim, S.C.; Eab, C.H. (Eds.) Fractional quantum fields. In Volume 5 Applications in Physics, Part B; De Gruyter: Berlin, Germany, 2019; pp. 237–256. [Google Scholar] [CrossRef]
  56. El-Nabulsi, A.R. Fractional derivatives generalization of Einstein’s field equations. Indian J. Phys. 2013, 87, 195–200. [Google Scholar] [CrossRef]
  57. El-Nabulsi, A.R. Non-minimal coupling in fractional action cosmology. Indian J. Phys. 2013, 87, 835–840. [Google Scholar] [CrossRef]
  58. Moniz, P.V.; Jalalzadeh, S. From Fractional Quantum Mechanics to Quantum Cosmology: An Overture. Mathematics 2020, 8, 313. [Google Scholar] [CrossRef]
  59. Rasouli, S.M.M.; Jalalzadeh, S.; Moniz, P.V. Broadening quantum cosmology with a fractional whirl. Mod. Phys. Lett. A 2021, 36, 2140005. [Google Scholar] [CrossRef]
  60. Moniz, P.V.; Jalalzadeh, S. Challenging Routes in Quantum Cosmology; World Scientific Publishing: Singapore, 2020. [Google Scholar] [CrossRef]
  61. Vacaru, S.I. Fractional Dynamics from Einstein Gravity, General Solutions, and Black Holes. Int. J. Theor. Phys. 2012, 51, 1338–1359. [Google Scholar] [CrossRef]
  62. Jalalzadeh, S.; da Silva, F.R.; Moniz, P.V. Prospecting black hole thermodynamics with fractional quantum mechanics. Eur. Phys. J. C 2021, 81, 632. [Google Scholar] [CrossRef]
  63. El-Nabulsi, R.A. Cosmology with a fractional action principle. Rom. Rep. Phys. 2007, 59, 763–771. [Google Scholar]
  64. El-Nabulsi, A.R. Fractional Lagrangian formulation of general relativity and emergence of complex, spinorial and noncommutative gravity. Int. J. Geom. Meth. Mod. Phys. 2009, 6, 25–76. [Google Scholar] [CrossRef]
  65. EL-Nabulsi, R.A. Fractional dynamics, fractional weak bosons masses and physics beyond the standard model. Chaos Solitons Fractals 2009, 41, 2262–2270. [Google Scholar] [CrossRef]
  66. Jamil, M.; Momeni, D.; Rashid, M.A. Fractional Action Cosmology with Power Law Weight Function. J. Phys. Conf. Ser. 2012, 354, 012008. [Google Scholar] [CrossRef]
  67. Shchigolev, V.K. Testing Fractional Action Cosmology. Eur. Phys. J. Plus 2016, 131, 256. [Google Scholar] [CrossRef]
  68. El-Nabulsi, R.A. non-local-in-time kinetic energy in nonconservative fractional systems, disordered dynamics, jerk and snap and oscillatory motions in the rotating fluid tube. Int. J.-Non-Linear Mech. 2017, 93, 65–81. [Google Scholar] [CrossRef]
  69. Rami, E.N.A. Fractional variational problems from extended exponentially fractional integral. Appl. Math. Comput. 2011, 217, 9492–9496. [Google Scholar] [CrossRef]
  70. Barrientos, E.; Mendoza, S.; Padilla, P. Extending Friedmann equations using fractional derivatives using a Last Step Modification technique: The case of a matter dominated accelerated expanding Universe. Symmetry 2021, 13, 174. [Google Scholar] [CrossRef]
  71. El-Nabulsi, R.A.; Torres, D.F. Fractional Action-Like Variational Problems. J. Math. Phys. 2008, 49, 053521. [Google Scholar] [CrossRef]
  72. Baleanu, D.; Trujillo, J. A new method of finding the fractional Euler-Lagrange and Hamilton equations within Caputo fractional derivatives. Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 1111–1115. [Google Scholar] [CrossRef]
  73. Odzijewicz, T.; Malinowska, A.; Torres, D. A generalized fractional calculus of variations. Comput. Math. Appl. 2012, 64, 3351–3366. [Google Scholar] [CrossRef]
  74. García-Aspeitia, M.A.; Fernandez-Anaya, G.; Hernández-Almada, A.; Leon, G.; Magaña, J. Cosmology under the fractional calculus approach. Mon. Not. Roy. Astron. Soc. 2022, 517, 4813–4826. [Google Scholar] [CrossRef]
  75. González, E.; Leon, G.; Fernandez-Anaya, G. Exact solutions and cosmological constraints in fractional cosmology. Fractal Fract. 2023, 7, 368. [Google Scholar] [CrossRef]
  76. Leon Torres, G.; García-Aspeitia, M.A.; Fernandez-Anaya, G.; Hernández-Almada, A.; Magaña, J.; González, E. Cosmology under the fractional calculus approach: A possible H0 tension resolution? PoS 2023, 248, 1–15. [Google Scholar] [CrossRef]
  77. Micolta-Riascos, B.; Millano, A.D.; Leon, G.; Erices, C.; Paliathanasis, A. Revisiting Fractional Cosmology. Fractal Fract. 2023, 7, 149. [Google Scholar] [CrossRef]
  78. Nojiri, S.; Odintsov, S.D.; Sasaki, M. Gauss–Bonnet dark energy. Phys. Rev. D 2005, 71, 123509. [Google Scholar] [CrossRef]
  79. Nojiri, S.; Odintsov, S.D.; Sami, M. Dark energy cosmology from higher-order, string-inspired gravity and its reconstruction. Phys. Rev. D 2006, 74, 046004. [Google Scholar] [CrossRef]
  80. Cognola, G.; Elizalde, E.; Nojiri, S.; Odintsov, S.; Zerbini, S. String-inspired Gauss–Bonnet gravity reconstructed from the Universe expansion history and yielding the transition from matter dominance to dark energy. Phys. Rev. D 2007, 75, 086002. [Google Scholar] [CrossRef]
  81. Nojiri, S.; Odintsov, S.D.; Tretyakov, P.V. Dark energy from modified F(R)-scalar-Gauss–Bonnet gravity. Phys. Lett. B 2007, 651, 224–231. [Google Scholar] [CrossRef]
  82. Chakraborty, S.; Paul, T.; SenGupta, S. Inflation driven by Einstein–Gauss–Bonnet gravity. Phys. Rev. D 2018, 98, 083539. [Google Scholar] [CrossRef]
  83. Fomin, I.V. Cosmological Inflation with Einstein–Gauss–Bonnet Gravity. Phys. Part. Nucl. 2018, 49, 525–529. [Google Scholar] [CrossRef]
  84. Kanti, P.; Gannouji, R.; Dadhich, N. Gauss–Bonnet Inflation. Phys. Rev. D 2015, 92, 041302. [Google Scholar] [CrossRef]
  85. Hikmawan, G.; Soda, J.; Suroso, A.; Zen, F.P. Comment on “Gauss–Bonnet inflation”. Phys. Rev. D 2016, 93, 068301. [Google Scholar] [CrossRef]
  86. Motaharfar, M.; Sepangi, H.R. Warm-tachyon Gauss–Bonnet inflation in the light of Planck 2015 data. Eur. Phys. J. C 2016, 76, 646. [Google Scholar] [CrossRef]
  87. Rashidi, N.; Nozari, K. Gauss–Bonnet Inflation after Planck2018. Astrophys. J. 2020, 890, 58. [Google Scholar] [CrossRef]
  88. Millano, A.D.; Leon, G.; Paliathanasis, A. Phase-Space Analysis of an Einstein–Gauss–Bonnet Scalar Field Cosmology. Mathematics 2023, 11, 1408. [Google Scholar] [CrossRef]
  89. Millano, A.D.; Leon, G.; Paliathanasis, A. Global dynamics in Einstein–Gauss–Bonnet scalar field cosmology with matter. Phys. Rev. D 2023, 108, 023519. [Google Scholar] [CrossRef]
  90. Millano, A.D.; Michea, C.; Leon, G.; Paliathanasis, A. Dynamics of a higher-dimensional Einstein–Scalar–Gauss–Bonnet cosmology. Phys. Dark Univ. 2024, 46, 101589. [Google Scholar] [CrossRef]
  91. Odintsov, S.; Oikonomou, V.; Fronimos, F. Rectifying Einstein–Gauss–Bonnet inflation in view of GW170817. Nucl. Phys. B 2020, 958, 115135. [Google Scholar] [CrossRef]
  92. Ratra, B.; Peebles, P.J.E. Cosmological Consequences of a Rolling Homogeneous Scalar Field. Phys. Rev. D 1988, 37, 3406. [Google Scholar] [CrossRef]
  93. Liddle, A.R.; Scherrer, R.J. A Classification of scalar field potentials with cosmological scaling solutions. Phys. Rev. D 1999, 59, 023509. [Google Scholar] [CrossRef]
  94. Uzan, J.P. Cosmological scaling solutions of nonminimally coupled scalar fields. Phys. Rev. D 1999, 59, 123510. [Google Scholar] [CrossRef]
  95. Frederico, G.A.S.F.; Torres, D. Necessary Optimality Conditions for Fractional Action-like Problems with Intrinsic and Observer Times. WSEAS Trans. Math. 2008, 7, 6–11. [Google Scholar]
  96. Goodman, J.; Weare, J. Ensemble samplers with affine invariance. Commun. Appl. Math. Comput. Sci. 2010, 5, 65–80. [Google Scholar] [CrossRef]
  97. Foreman-Mackey, D.; Hogg, D.W.; Lang, D.; Goodman, J. emcee: The MCMC Hammer. Publ. Astron. Soc. Pac. 2013, 125, 306–312. [Google Scholar] [CrossRef]
  98. Perko, L. Differential Equations and Dynamical Systems, 3rd ed.; Springer: Berlin/Heidelberg, Germany, 2001. [Google Scholar]
  99. Capozziello, S.; D’Agostino, R.; Luongo, O. Cosmographic analysis with Chebyshev polynomials. Mon. Not. R. Astron. Soc. 2018, 476, 3924–3938. [Google Scholar] [CrossRef]
  100. Jimenez, R.; Loeb, A. Constraining cosmological parameters based on relative galaxy ages. Astrophys. J. 2002, 573, 37–42. [Google Scholar] [CrossRef]
  101. Brout, D.; Scolnic, D.; Popovic, B.; Riess, A.G.; Carr, A.; Zuntz, J.; Kessler, R.; Davis, T.M.; Hinton, S.; Jones, D.; et al. The Pantheon+ Analysis: Cosmological Constraints. Astrophys. J. 2022, 938, 110. [Google Scholar] [CrossRef]
  102. Tripp, R. A Two-parameter luminosity correction for type Ia supernovae. Astron. Astrophys. 1998, 331, 815–820. [Google Scholar]
  103. Kessler, R.; Scolnic, D. Correcting Type Ia Supernova Distances for Selection Biases and Contamination in Photometrically Identified Samples. Astrophys. J. 2017, 836, 56. [Google Scholar] [CrossRef]
  104. Riess, A.G.; Yuan, W.; Macri, L.M.; Scolnic, D.; Brout, D.; Casertano, S.; Jones, D.O.; Murakami, Y.; An, G.S.; Breuval, L.; et al. A Comprehensive Measurement of the Local Value of the Hubble Constant with 1 km s−1 Mpc−1 Uncertainty from the Hubble Space Telescope and the SH0ES Team. Astrophys. J. Lett. 2022, 934, L7. [Google Scholar] [CrossRef]
  105. Treu, T.; Marshall, P.J. Time Delay Cosmography. Astron. Astrophys. Rev. 2016, 24, 11. [Google Scholar] [CrossRef]
  106. Wong, K.C.; Suyu, S.H.; Chen, G.C.-F.; Rusu, C.E.; Millon, M.; Sluse, D.; Bonvin, V.; Fassnacht, C.D.; Taubenberger, S.; Auger, M.W.; et al. H0LiCOW – XIII. A 2.4 per cent measurement of H0 from lensed quasars: 5.3σ tension between early- and late-Universe probes. Mon. Not. R. Astron. Soc. 2020, 498, 1420–1439. [Google Scholar] [CrossRef]
  107. Jee, I.; Suyu, S.H.; Komatsu, E.; Fassnacht, C.D.; Hilbert, S.; Koopmans, L.V.E. A measurement of the Hubble constant from angular diameter distances to two gravitational lenses. Science 2019, 365, 1134–1138. [Google Scholar] [CrossRef] [PubMed]
  108. Birrer, S.; Treu, T.; Rusu, C.E.; Bonvin, V.; Fassnacht, C.D.; Chan, J.H.H.; Agnello, A.; Shajib, A.J.; Chen, G.C.-F.; Auger, M.; et al. H0LiCOW - IX. Cosmographic analysis of the doubly imaged quasar SDSS 1206+4332 and a new measurement of the Hubble constant. Mon. Not. R. Astron. Soc. 2019, 484, 4726. [Google Scholar] [CrossRef]
  109. Rusu, C.E.; Wong, K.C.; Bonvin, V.; Sluse, D.; Suyu, S.H.; Fassnacht, C.D.; Chan, J.H.; Hilbert, S.; Auger, M.W.; Sonnenfeld, A.; et al. H0LiCOW XII. Lens mass model of WFI2033 - 4723 and blind measurement of its time-delay distance and H0. Mon. Not. R. Astron. Soc. 2020, 498, 1440–1468. [Google Scholar] [CrossRef]
  110. Chen, G.C.; Fassnacht, C.D.; Suyu, S.H.; Rusu, C.E.; Chan, J.H.; Wong, K.C.; Auger, M.W.; Hilbert, S.; Bonvin, V.; Birrer, S.; et al. A SHARP view of H0LiCOW: H0 from three time-delay gravitational lens systems with adaptive optics imaging. Mon. Not. R. Astron. Soc. 2019, 490, 1743–1773. [Google Scholar] [CrossRef]
  111. Escamilla-Rivera, C.; Torres Castillejos, R. H0 Tension on the Light of Supermassive Black Hole Shadows Data. Universe 2023, 9, 14. [Google Scholar] [CrossRef]
  112. Akiyama, K.; Alberdi, A.; Alef, W.; Asada, K.; Azuly, R. First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole. Astrophys. J. Lett. 2019, 875, L1. [Google Scholar] [CrossRef]
  113. Akiyama, K.; Alberdi, A.; Alef, W.; Algaba, J.C.; Anantua, R.; Asada, K.; Azulay, R.; Bach, U.; Baczko, A.K.; Ball, D.; et al. First Sagittarius A* Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole in the Center of the Milky Way. Astrophys. J. Lett. 2022, 930, L12. [Google Scholar] [CrossRef]
  114. Krishnan, C.; Ó Colgáin, E.; Sheikh-Jabbari, M.M.; Yang, T. Running Hubble tension and a H0 diagnostic. Phys. Rev. D 2021, 103, 103509. [Google Scholar] [CrossRef]
  115. Riess, A.G.; Casertano, S.; Yuan, W.; Macri, L.M.; Scolnic, D. Large Magellanic Cloud Cepheid Standards Provide a 1% Foundation for the Determination of the Hubble Constant and Stronger Evidence for Physics beyond ΛCDM. Astrophys. J. 2019, 876, 85. [Google Scholar] [CrossRef]
  116. Gupta, R.P. JWST early Universe observations and ΛCDM cosmology. Mon. Not. R. Astron. Soc. 2023, 524, 3385–3395. [Google Scholar] [CrossRef]
  117. Valcin, D.; Jimenez, R.; Verde, L.; Bernal, J.L.; Wandelt, B.D. The age of the Universe with globular clusters: Reducing systematic uncertainties. J. Cosmol. Astropart. Phys. 2021, 2021, 17. [Google Scholar] [CrossRef]
Figure 1. Allowed regions for the parameters μ and m. The red star represents the point in the parameter space with the best-fit values of μ and m according to the analysis in Section 5.
Figure 1. Allowed regions for the parameters μ and m. The red star represents the point in the parameter space with the best-fit values of μ and m according to the analysis in Section 5.
Fractalfract 08 00626 g001
Figure 2. Allowed regions by the parameters μ and m, for the sink, source, and saddle cases. The red star represents the point in the parameter space with the best-fit values of μ and m according to the analysis in Section 5.
Figure 2. Allowed regions by the parameters μ and m, for the sink, source, and saddle cases. The red star represents the point in the parameter space with the best-fit values of μ and m according to the analysis in Section 5.
Fractalfract 08 00626 g002
Figure 3. Flow of system (70) and (71) for the sink cases. The streamlines are colored by default according to the magnitude of the vector field, with the arrow pointing in the time-increasing direction.
Figure 3. Flow of system (70) and (71) for the sink cases. The streamlines are colored by default according to the magnitude of the vector field, with the arrow pointing in the time-increasing direction.
Fractalfract 08 00626 g003
Figure 4. Flow of system (70) and (71) for the source cases. The streamlines are colored by default according to the magnitude of the vector field, with the arrow pointing in the time-increasing direction.
Figure 4. Flow of system (70) and (71) for the source cases. The streamlines are colored by default according to the magnitude of the vector field, with the arrow pointing in the time-increasing direction.
Fractalfract 08 00626 g004
Figure 5. Flow of system (70) and (71) for the saddle cases. The streamlines are colored by default according to the magnitude of the vector field, with the arrow pointing in the time-increasing direction.
Figure 5. Flow of system (70) and (71) for the saddle cases. The streamlines are colored by default according to the magnitude of the vector field, with the arrow pointing in the time-increasing direction.
Fractalfract 08 00626 g005
Figure 6. Real part of the eigenvalues for B and C for a fixed value of one of the parameters according to the respective case. In cases (1) and (2) μ is fixed as μ = 1 2 5 17 while m moves in two different intervals. In case (34), μ = 35 8 and 4 5 < m < 16 . Both points behave as sources or saddles whenever they exist and are hyperbolic.
Figure 6. Real part of the eigenvalues for B and C for a fixed value of one of the parameters according to the respective case. In cases (1) and (2) μ is fixed as μ = 1 2 5 17 while m moves in two different intervals. In case (34), μ = 35 8 and 4 5 < m < 16 . Both points behave as sources or saddles whenever they exist and are hyperbolic.
Fractalfract 08 00626 g006
Figure 7. Real part of the eigenvalues for B and C where both of the parameters μ , m remain free in some intervals. The cases depicted are cases (5), (37) and (38). Once again, as in Figure 6, both points behave either as sources or saddles for every existence condition in which they are hyperbolic.
Figure 7. Real part of the eigenvalues for B and C where both of the parameters μ , m remain free in some intervals. The cases depicted are cases (5), (37) and (38). Once again, as in Figure 6, both points behave either as sources or saddles for every existence condition in which they are hyperbolic.
Fractalfract 08 00626 g007
Figure 8. Real part of the eigenvalues for D and E for a fixed value of one of the parameters according to the respective case. In cases (1) and (2), μ is fixed as μ = 1 2 5 17 while m moves in two different intervals. In case (34), μ = 35 8 and 4 5 < m < 16 . Both points behave as sinks or saddles whenever they exist and are hyperbolic.
Figure 8. Real part of the eigenvalues for D and E for a fixed value of one of the parameters according to the respective case. In cases (1) and (2), μ is fixed as μ = 1 2 5 17 while m moves in two different intervals. In case (34), μ = 35 8 and 4 5 < m < 16 . Both points behave as sinks or saddles whenever they exist and are hyperbolic.
Fractalfract 08 00626 g008
Figure 9. Real part of the eigenvalues for D and E where both of the parameters μ , m remain free in some intervals. The cases depicted are cases (5), (37) and (38). Once again, as in Figure 8, both points behave as sinks or saddles whenever they exist and are hyperbolic.
Figure 9. Real part of the eigenvalues for D and E where both of the parameters μ , m remain free in some intervals. The cases depicted are cases (5), (37) and (38). Once again, as in Figure 8, both points behave as sinks or saddles whenever they exist and are hyperbolic.
Fractalfract 08 00626 g009
Figure 10. Existence region for points B, C, D and E. The red star represents the point in the parameter space with the best-fit values of μ and m according to the analysis in Section 5.
Figure 10. Existence region for points B, C, D and E. The red star represents the point in the parameter space with the best-fit values of μ and m according to the analysis in Section 5.
Fractalfract 08 00626 g010
Figure 11. Dynamics of system (75) and (76) in the compactified region X 2 + Y 2 1 for different values of the parameters μ and m. We observe that for ( μ = 1 , m = 3 ) or ( μ = 1 , m = 4 ) or ( μ = 2 , m = 3 ) , point A is an attractor. However, for ( μ = 2 , m = 4 ) , there are two complex eigenvalues with zero real parts, namely i 3 2 , i 3 2 which means A is a center. For ( μ = 3 , m = 3 ) , the point A is a source, B and C are saddles and D and E are sinks. Finally for ( μ = 3 , m = 4 ) , A is a saddle point, B and C are sources and D and E are sinks. The streamlines are colored by default according to the magnitude of the vector field, with the arrow pointing in the time-increasing direction.
Figure 11. Dynamics of system (75) and (76) in the compactified region X 2 + Y 2 1 for different values of the parameters μ and m. We observe that for ( μ = 1 , m = 3 ) or ( μ = 1 , m = 4 ) or ( μ = 2 , m = 3 ) , point A is an attractor. However, for ( μ = 2 , m = 4 ) , there are two complex eigenvalues with zero real parts, namely i 3 2 , i 3 2 which means A is a center. For ( μ = 3 , m = 3 ) , the point A is a source, B and C are saddles and D and E are sinks. Finally for ( μ = 3 , m = 4 ) , A is a saddle point, B and C are sources and D and E are sinks. The streamlines are colored by default according to the magnitude of the vector field, with the arrow pointing in the time-increasing direction.
Fractalfract 08 00626 g011aFractalfract 08 00626 g011b
Figure 12. Posterior 1D distribution and joint marginalized regions of the free parameter space of the Λ CDM cosmology for the SNe Ia, SNe Ia+CC, and SNe Ia+CC+GL+BHS (joint) data. The admissible joint regions correspond to 1 σ , 2 σ , and 3 σ CL, respectively. The best-fit values for each model’s free parameter are shown in Table 2.
Figure 12. Posterior 1D distribution and joint marginalized regions of the free parameter space of the Λ CDM cosmology for the SNe Ia, SNe Ia+CC, and SNe Ia+CC+GL+BHS (joint) data. The admissible joint regions correspond to 1 σ , 2 σ , and 3 σ CL, respectively. The best-fit values for each model’s free parameter are shown in Table 2.
Fractalfract 08 00626 g012
Figure 13. Posterior 1D distribution and joint marginalized regions of the free parameter space of the fractional cosmology for the SNe Ia, SNe Ia+CC, and SNe Ia+CC+GL+BHS (joint) data. The admissible joint regions correspond to 1 σ , 2 σ , and 3 σ CL, respectively. The best-fit values for each model free parameter are shown in Table 2.
Figure 13. Posterior 1D distribution and joint marginalized regions of the free parameter space of the fractional cosmology for the SNe Ia, SNe Ia+CC, and SNe Ia+CC+GL+BHS (joint) data. The admissible joint regions correspond to 1 σ , 2 σ , and 3 σ CL, respectively. The best-fit values for each model free parameter are shown in Table 2.
Fractalfract 08 00626 g013
Figure 14. Theoretical Hubble parameter for the Λ CDM cosmology (red dashed line) and the fractional cosmology for the exponential potential (solid blue line) as a function of the redshift z, contrasted with the CC sample. The shaded curve represents the confidence regions of the Hubble parameter at a 1 σ CL. The figure was obtained using the chains of the MCMC procedure described in Section 5 for the joint analysis.
Figure 14. Theoretical Hubble parameter for the Λ CDM cosmology (red dashed line) and the fractional cosmology for the exponential potential (solid blue line) as a function of the redshift z, contrasted with the CC sample. The shaded curve represents the confidence regions of the Hubble parameter at a 1 σ CL. The figure was obtained using the chains of the MCMC procedure described in Section 5 for the joint analysis.
Fractalfract 08 00626 g014
Figure 15. Deceleration parameters for the Λ CDM cosmology (red dashed line) and the fractional cosmology for the exponential potential (solid blue line) as a function of the redshift z. The shaded curves represent the confidence region of the deceleration parameter at a 1 σ CL. The figure was obtained using the chains of the MCMC procedure described in Section 5 for the joint analysis.
Figure 15. Deceleration parameters for the Λ CDM cosmology (red dashed line) and the fractional cosmology for the exponential potential (solid blue line) as a function of the redshift z. The shaded curves represent the confidence region of the deceleration parameter at a 1 σ CL. The figure was obtained using the chains of the MCMC procedure described in Section 5 for the joint analysis.
Fractalfract 08 00626 g015
Figure 16. H 0 diagnostics for the Λ CDM cosmology (red dashed line) and the fractional cosmology for the exponential potential (solid blue line) as a function of the redshift z. The shaded curve represents the confidence regions of the H 0 diagnostic at a 1 σ CL. The figure was obtained using the chains of the MCMC procedure described in Section 5 for the joint analysis.
Figure 16. H 0 diagnostics for the Λ CDM cosmology (red dashed line) and the fractional cosmology for the exponential potential (solid blue line) as a function of the redshift z. The shaded curve represents the confidence regions of the H 0 diagnostic at a 1 σ CL. The figure was obtained using the chains of the MCMC procedure described in Section 5 for the joint analysis.
Fractalfract 08 00626 g016
Table 1. Total number of steps and autocorrelation time τ corr for the free parameter space of the Λ CDM and the fractional cosmology, respectively. These values were obtained when the convergence test described in Section 5 is fulfilled for an MCMC analysis with 100 chains and the flat priors 0.55 < h < 0.85 , 0 < Ω m , 0 < 1 , 0.5 < α 0 < 2.5 , and 20 < M < 18 . This information is provided so that our results are replicable.
Table 1. Total number of steps and autocorrelation time τ corr for the free parameter space of the Λ CDM and the fractional cosmology, respectively. These values were obtained when the convergence test described in Section 5 is fulfilled for an MCMC analysis with 100 chains and the flat priors 0.55 < h < 0.85 , 0 < Ω m , 0 < 1 , 0.5 < α 0 < 2.5 , and 20 < M < 18 . This information is provided so that our results are replicable.
DataTotal Steps τ corr
h Ω m , 0 α 0 M
Λ CDM cosmology
SNe Ia1250 24.3 22.7 24.3
SNe Ia+CC1300 25.2 24.3 25.0
joint1350 23.4 24.6 23.6
Fractional cosmology
SNe Ia1300 24.5 25.0 24.6
SNe Ia+CC1600 25.0 26.3 25.5
joint1350 24.8 23.0 24.3
Table 2. Best-fit values and χ min 2 criteria of the Λ CDM and the fractional cosmology respectively, for the SNe Ia, SNe Ia+CC, and SNe Ia+CC+GL+BHS (joint) data. The uncertainties presented correspond to 1 σ CL.
Table 2. Best-fit values and χ min 2 criteria of the Λ CDM and the fractional cosmology respectively, for the SNe Ia, SNe Ia+CC, and SNe Ia+CC+GL+BHS (joint) data. The uncertainties presented correspond to 1 σ CL.
DataBest-Fit Values χ min 2
h Ω m , 0 α 0 M
Λ CDM cosmology
SNe Ia 0.734 ± 0.010 0.333 ± 0.018 19.25 ± 0.03 1523
SNe Ia+CC 0.719 ± 0.009 0.314 ± 0.016 19.30 ± 0.03 1547
Joint 0.724 ± 0.008 0.311 ± 0.016 19.28 ± 0.02 1681
Fractional cosmology
SNe Ia 0.729 ± 0.010 1.41 ± 0.06 19.25 ± 0.03 1532
SNe Ia+CC 0.717 ± 0.009 1.39 ± 0.05 19.28 ± 0.02 1564
Joint 0.712 ± 0.007 1.38 ± 0.05 19.29 ± 0.02 1706
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

Micolta-Riascos, B.; Millano, A.D.; Leon, G.; Droguett, B.; González, E.; Magaña, J. Fractional Einstein–Gauss–Bonnet Scalar Field Cosmology. Fractal Fract. 2024, 8, 626. https://doi.org/10.3390/fractalfract8110626

AMA Style

Micolta-Riascos B, Millano AD, Leon G, Droguett B, González E, Magaña J. Fractional Einstein–Gauss–Bonnet Scalar Field Cosmology. Fractal and Fractional. 2024; 8(11):626. https://doi.org/10.3390/fractalfract8110626

Chicago/Turabian Style

Micolta-Riascos, Bayron, Alfredo D. Millano, Genly Leon, Byron Droguett, Esteban González, and Juan Magaña. 2024. "Fractional Einstein–Gauss–Bonnet Scalar Field Cosmology" Fractal and Fractional 8, no. 11: 626. https://doi.org/10.3390/fractalfract8110626

APA Style

Micolta-Riascos, B., Millano, A. D., Leon, G., Droguett, B., González, E., & Magaña, J. (2024). Fractional Einstein–Gauss–Bonnet Scalar Field Cosmology. Fractal and Fractional, 8(11), 626. https://doi.org/10.3390/fractalfract8110626

Article Metrics

Back to TopTop