Estimated Age of the Universe in Fractional Cosmology

Our proposed cosmological framework, which is based on fractional quantum cosmology, aims to address the issue of synchronicity in the age of the universe. To achieve this, we have developed a new fractional Λ CDM cosmological model. We obtained the necessary formalism by obtaining the fractional Hamiltonian constraint in a general minisuperspace. This formalism has allowed us to derive the fractional Friedmann and Raychaudhuri equations for a homogeneous and isotropic cosmology. Unlike the traditional de Sitter phase, our model exhibits a power-law accelerated expansion in the late-time universe, when vacuum energy becomes dominant. By fitting the model’s parameters to cosmological observations, we determined that the fractional parameter of L´evy equals α = 1 . 986. Additionally, we have calculated the age of the universe to be 13.8196 Gyr. Furthermore, we have found that the ratio of the age to Hubble time from the present epoch to the distant future is finite and confined within the interval 0 . 9858 ≤ Ht < 95 . 238.


Introduction
Estimates for the age of the universe have been made over time.Einstein and de Sitter's model [1] predicted an age of the universe that was roughly 9 Gyr, which is quite different from the age of the ancient stars.The age of ancient stars has been estimated to be larger than 12 Gyr, as per the research by Chaboyer [2].This discrepancy led to the development of new models that could account for the universe's observed age and expansion rate.Despite the universe experiencing deceleration for the first 9 billion years and then transitioning to acceleration through cosmic jerk for the past 5 billion years [3], the present value of the dimensionless age, H 0 t 0 , where H 0 and t 0 are the current value of the Hubble parameter and the universe's age, respectively, is limited by supernova distances to be very close to 1.For example, Tonry et al. [4] estimated this to be H 0 t 0 = 0.96 ± 0.04.We seem to be experiencing a period of privilege at the moment, wherein the age of the universe in terms of Hubble time, especially in the case of a ΛCDM Universe, is extremely close to unity.This coincidence is known as the synchronicity problem of the universe's age [3,5].
Various other methods have verified this unexpected closeness to unity, which has been the subject of intense speculation and debate.The synchronicity problem has been investigated by cosmological models with the inclusion of scalar fields [6,7], effective dependency between dark energy and dark matter [8,9,10,11,12], interacting holographic dark energy [13,14,15], and other interpretations at the level of standard cosmology [16,17].From another perspective, quantum cosmology in the paradigm of the Wheeler-DeWitt (WDW) equation in the minisuperspace of a homogeneous and isotropic universe was used to propose a solution to the question in a unified way to other gaps in the ΛCDM model [18].Proposals that modify the Friedmann-Lamaître-Robertson-Walker (FLRW) cosmology were adopted to solve cosmological problems.They have a linear or quasi-linear evolution of the universe in common and present synchronicity Additionally, some proposals in cosmology, including the widely accepted Kolb model (1989) [19], explain the presence of K-matter, a unique form of exotic matter that dominates at low redshifts.K-matter has an equation of state with a value of ω = −1/3, and its density decreases as the universe expands, following an inverse square relation with the scale factor.Another proposal worth mentioning is the Allen model [20], which incorporates a S U(2) cosmological instanton in addition to the Einstein field equations.The instanton dominates during the later stages of the Universe's evolution, resulting in an approximately proportional relationship between the scale factor and cosmic time t.Moving on, Pimentel and Diaz-Rivera [21] presented a solution to the synchronicity problem by considering the effect of a time-dependent cosmological constant Λ in a family of scalar-tensor theories.Interestingly, the cosmological constant in these models exhibits an inverse square relationship with time, reaching its minimum value in the present era.However, a point of concern arises when considering the gravitational constant G, which decreases inversely with time and contradicts the experimental constraints on time variation in G.Moreover, Melia and Shevchuk [22] raised an intriguing question regarding the coincidence between the cosmic time t 0 = H −1 0 , describing the distance that light has traveled since the Big Bang and the gravitational horizon R h (t 0 ).The query was whether these two quantities were indeed equal or not.However, Lewis [23] criticized this model by highlighting that if the dark energy component of the universe possesses an equation of state with ω = −1/3, the addition of any amount of matter will inevitably lead to a deviation from the strictly linear evolution of the R h = ct universe.This deviation would manifest the universe transitioning to a matter-dominated phase during its earliest epochs.
The primary aim of this study is to establish an appropriate fractional cosmological framework to examine the synchronicity issue, among other enigmas linked to the ΛCDM model.Our investigation endeavors to tackle the synchronicity problem through the fractional quantum cosmology (FQC) proposal [24,25,26,27,28,18,29].The basic and preliminary methodology of the FQC is to combine fractional quantum mechanics (FQM) [30,31,32] and quantum cosmology (QC) [33] to achieve the corresponding cosmological consequences (see [29] for a detailed discussion).However, it should be mentioned that in this study, we only consider the (semi-)classical limit of this theory; to see other semiclassical models established applying various gravitational theories as underlying frameworks, see, for instance, [34,35,36,37,38,39,40,41,42].Several studies have put forth the idea of fractional classical cosmology (FCC) [43,44,45], since the application of fractional calculus in cosmology is now producing encouraging and promising results [46].
There have been numerous instances in the fields of gravity and cosmology where fractional calculus has been employed, showcasing the broad and significant applications of this mathematical tool.The current research in this domain is thriving, with scientists actively exploring the potential of fractional calculus.Its effectiveness as a valuable instrument becomes apparent when confronting a wide range of problems associated with gravitational forces and cosmolog-2 ical models, as evidenced by the wealth of studies conducted [46,45,47,48,49,50,51,52,53,54,55,56,57].Moreover, fractional calculus finds application in various other areas, such as the study of the stochastic gravitational-wave background in quantum gravity [58], the investigation of gravitational-wave luminosity distance [59], the exploration of inflation and CMB spectrum [28,60,61], the development of fractional action cosmology [62,63,64], the analysis of the fractional geodesic equation and discrete gravity [65], and the examination of nonminimal coupling and chaotic inflation [66].Additionally, fractional calculus has been utilized in investigating phantom cosmology with conformal coupling [67], Ornstein-Uhlenbeck-like fractional differential equations in cosmology [68], fractional action cosmology with a variableorder parameter [69], and wormholes in fractional action cosmology [70].Notably, new metrics have been considered [71], while the application of fractional calculus has also extended to the study of dark energy models in emergent, logamediate, and intermediate scenarios of the universe [72,73].As an example, the work of [50,52] yielded a value of α = 0.926 for the order of the Riemann-Liouville fractional integral, providing further insights into the intricacies of the field.
In [43,47,48], several exact solutions were obtained for cosmological models, which, due to the fractal nature of space-time, exhibit a significant deviation from the standard model [74,75].The investigation carried out in ref. [27] delves into the interval 1 ≤ α < 2 and employs Riesz's fractional derivative to derive the nonboundary and tunneling wave functions for a closed de Sitter geometry.Furthermore, ref. [28] explores the preinflation epoch within the framework of fractional quantum cosmology.With regard to the thermodynamics of fractional BHs, this subject has been thoroughly examined in ref. [27], and the findings from this research have been effectively utilized to obtain the Friedmann equations in the context of emergent gravity in ref. [76].It is worth mentioning that fractional calculus has also been employed to modify the Friedmann and Raychaudhuri equations, thereby enabling an investigation into the dynamics of the universe without the need for cold dark matter (CDM) or dark energy [77].In addition, an alternative approach involves the utilization of fractional calculus to determine the value of the cosmological constant, which necessitates a restructuring due to the well-recognized ultraviolet divergence in conventional quantum field theory [78,79,80].To expand the scope of this fractional approach, [81,82] explore modified Newtonian dynamics (MOND) and quantum cosmology.Lastly, it is crucial to note that there exists a multitude of definitions for fractional derivatives and fractional integrals, including those proposed by Riemann-Liouville, Caputo, Riesz, Hadamard, Marchand, and Griinwald-Letnikov, as well as more recent formulations (as comprehensively documented in [83] and [84], along with their respective references).Despite the extensive research conducted on these operators, they do not always exhibit the conventional properties associated with function differentiation, such as Leibniz's rule, the chain rule, and the semigroup property [84,83].
This paper is organized as follows: We begin in Section 2 by introducing fractional quantum cosmology, which allows us to define the fractional extension of the ADM Hamiltonian and Hamiltonian constraint within the minisuperspace approximation of gravity.Next, in Section 3, we apply the general formalism obtained in Section 3 to the FLRW cosmology.This leads to deriving the fractional extensions of the Friedmann and Raychaudhuri equations.We then proceed with a thorough discussion of various fractional cosmological parameters, including the Hubble and deceleration parameters.In Section 4, we briefly review the available cosmological datasets.Afterwards, we consider the observational datasets, such as SNIa, CMB, BAO, BBN, and OHD, to estimate both the universe's age and the model's free parameters.Section 5 focuses on the perspective of fractional cosmology regarding the synchronicity problem.Finally, in 3 Section 6, we summarize the results and offer additional discussions.

Fractional Quantum Cosmology
This section briefly explains the FQC framework.To develop a fractional WDW equation, we first briefly describe how to construct a fractional Schrödinger equation (SE).
In two seminal papers, Laskin [30,85] introduced a groundbreaking concept known as space fractional quantum mechanics (SFQM) through the application of a purely physics-based approach.Laskin's revolutionary work was built upon the solid foundation of nonrelativistic quantum mechanics as formulated by Feynman and Hibbs, which extensively utilized a path integral representation employing Brownian paths.However, Laskin's innovative contribution to the field was to replace these conventional Brownian paths with paths characterized by Lévy flights, thus giving rise to the fascinating realm of SFQM and expanding our understanding of quantum phenomena in space.
In the framework of SFQM, one can initiate the analysis by considering the conventional time-dependent SE, which is expressed as follows: This equation serves as a starting point for further exploration, wherein one can delve into its fractional counterpart, as proposed by Laskin [31].The fractional SE is obtained by replacing the traditional Laplace operator ∆ with the Riesz fractional derivative (−ℏ 2 ∆) α/2 , which can be expressed as This replacement paves the way for the exploration of the space-fractional SE: where α ∈ (1, 2] is the Lévy index, and D α is, according to Laskin, a generalized "fractional quantum diffusion coefficient" of dimension J 1−α m α s −α .As mentioned earlier, (−ℏ 2 ∆) α/2 is a generalization of the fractional Riesz derivative, known as the fractional Laplacian [86]: which, by means of Fourier transforms, relates ψ(r, t) and φ(p, t).For the special case when α = 2, D α reduces to D 2 = 1/2m, where m is the particle's mass.The Riesz fractional derivative (the fractional Laplacian), in general, characterizes a unique diffusion process that arises from the random displacements of individuals who can move to nearby or neighboring locations and even venture to distant sites by means of Lévy flights.Both literal and conceptual flights have been observed or claimed to occur in various scenarios, such as turbulent fluid motion and the transportation of materials in fractured media.In the field of mechanics, the fractional Laplacian describes the motion of a chain or group of particles interconnected by elastic springs that connect not only to their immediate neighbors but also to all other particles.The strength of the springs decreases as the particles move apart, while the particle arrangement can take on either a regular or fractal pattern.In a more abstract sense, the fractional Laplacian represents the influence of a nonlocal process on a conservation law, which is not only influenced by local conditions but also by the overall state of a field of interest at a specific moment in time.We suggest that the reader consults Pozrikidis's exceptionally well-written and comprehensive book [86] in order to obtain detailed explanations and develop a deeper understanding of fractional Laplacian.
On the other hand, as demonstrated by Laskin in his work [87], the evolution of the fractional wave function (3) can be effectively described by a kernel.This kernel is expressed as an integral equation, which can be written as follows: In this equation, ψ(r a , t a ) represents the wave function of the initial state, while ψ(r b , t b ) represents the wave function of the final state.The kernel itself, denoted as K L (r b , t b |r a , t a ), is given by the following expression: Here, it is important to note that H α (r, p, t) serves as the fractional extension for the Hamiltonian of the system under consideration.This Hamiltonian can be mathematically expressed as It is worth mentioning that in the particular scenario where α takes a value of 2 and D α is equal to 1/2m, the equation ( 7) simplifies to the standard Hamiltonian, which is defined as We can further expand the aforementioned procedure in order to acquire fractional quantum cosmology.In cosmology, we assume that the space-time metric and matter field(s) exhibit homogeneity and isotropy.Consequently, the lapse function is considered to be homogeneous, denoted as N = N(t), while the shift is assigned a value of zero, N i = 0.As a result, the space-time line element reduces to Inserting the restricted form of the line element (8) into the ADM action functional, one generally obtains the following reduced form of the action [29]: In this equation, f µν denotes the reduced DeWitt metric of n-dimensional minisuperspace with signature (−, +, ..., +).Also, U(q), the potential, includes the Ricci scalar of the threedimensional t = const.hypersurfaces and the potential of matter fields.Generally, the minisuperspace coordinates q µ may include matter fields and three-metric components.Variation in the action with respect to q A yields the Euler-Lagrange equations where µ αβ are the connection coefficients (Christoffel symbols) determined from the metric of the minisuperspace f µν .Also, in a variation with respect to the lapse function, N, one obtains the constraint 1 2N To ensure consistency, it is necessary for Equations ( 10) and ( 11) to be equal to the 00 and i j components of the complete Einstein equations, respectively.
To define the Hamiltonian, we first need to define canonical momenta in a usual way, which is This leads us to the canonical ADM Hamiltonian in the minisuperspace where f µν is the inverse metric of the minisuperspace.The lapse function N serves as a Lagrange multiplier with the purpose of ensuring the satisfaction of the Hamiltonian constraint In the coordinate representation of canonical quantization, the process involves the incorporation of a wave function, denoted as Ψ(q), which is independent of time and necessitates the condition that it is rendered null by the operator associated with the aforementioned Hamiltonian constraint.As a result, the WDW equation is obtained: where is the curvature of the minisuperspace metric f µν , and ζ is an arbitrary constant [29].It is important to note that the minisuperspace metric's dependence on the local coordinates of the minisuperspace leads to a nontrivial operator ordering issue in the resulting WDW equation.However, this issue can be somewhat addressed by ensuring that the quantization process is covariant in minisuperspace.This means that it should not be affected by general coordinate transformations in the minisuperspace.The factor ordering procedure is the origin of the minisuperspace curvature term, denoted by R.
To acquire the fractional equivalent of the WDW equation (15), it is necessary to substitute the typical d'Alembertian operator with the fractional Riesz-d'Alembertian operator given by [88,89] (−∇) where |Π| = f µν Π µ Π ν , and F represents a Fourier transformation.Therefore, the fractional equivalent of the WDW equation ( 15) can be expressed as [24,25,26,27,29] where D α (m P ) is the coefficient, depending on the Planck mass m P = 1/ √ G, which reduces to unity at α = 2 limit, and the value of Lévy's fractional parameter, denoted as α, is assigned the range of values 1 < α ≤ 2.
By following the methodology described by Laskin [87], it becomes feasible to illustrate the comprehensive procedure for tje fractional path integral quantization of constrained systems, as expounded upon by Balatin, Fradkin, and Vilkovisky (for further elaboration, please refer to [90]).This approach culminates in the derivation of the resulting propagation amplitude which arises when this technique is applied to the minisuperspace model that is precisely defined by the Hamiltonian ( 14) where H α (q µ , Π µ ) is the fractional Hamiltonian constraint, counterpart to (14), defined by This equation gives us the fractional extension of Equation ( 11).Also, one can easily obtain the fractional Euler-Lagrange counterpart of ( 10) by Hamilton's equations where H (α) ADM is the fractional ADM Hamiltonian defined by In the following section, we derive the fractional ΛCDM cosmology using Equations ( 19) and (20).

Fractional FLRW Cosmology
Let us commence our discussion by considering the FLRW metric, which serves as the underlying geometric structure In this metric, we can observe various components contributing to the overall structure.Firstly, the term N(t) represents the lapse function, which characterizes the rate at which time elapses within the cosmic framework.Secondly, the variable t represents the cosmic time, enabling us to measure and understand the progression of events within the universe.Furthermore, the term dΩ 2  2 stands for the standard line element for S 2 , encompassing both the angular coordinates θ and ϕ.Lastly, the parameter k assumes values of −1, 0, and 1, signifying the sign of spatial curvature.These values correspond to an open, flat, and closed universe, respectively.
The functional, which encompasses both gravitational and matter components, with the matter field being treated as a perfect fluid [91,92,93,94,95,36,96], can be expressed as follows: where K denotes the trace of the extrinsic curvature of the spacelike submanifolds, (Σ, h), of the space-time, (M, g), in which M = R × Σ and M represents the manifold of the space-time with boundary ∂M.In addition, ρ is the total energy density of the cosmic fluid with components where ρ (r) , ρ (m) , and ρ (Λ) = Λ 8πG are the energy density of the radiation, matter density of the cosmic dust (constitute of baryonic and cold matter), and the energy density of cosmological constant, respectively.
The act of introducing the homogeneous and isotropic metric (22) into the action functional (23) results in a simplified Arnowitt-Deser-Misner (ADM) form as follows [29]: where is the spatial volume of the spacelike hypersurfaces (Σ, h).
To ensure that the spacelike sections, represented by Σ, have a finite volume, it is assumed that they are both compact and without boundary.In the scientific community, it is widely accepted that any compact Riemannian three-manifold with constant curvature is homeomorphic to Σ/Γ, where Γ represents the group of covering transformations, and Σ represents the universal covering space.The shape of Σ can take on different forms depending on the sign of the spatial curvature k (where k can be 0, 1, or -1), such as R 3 (representing three-dimensional Euclidean space), S 3 (representing a 3-sphere), or H 3 (representing three-dimensional hyperbolic space) [29].
Using this particular approach, the action mentioned above is clearly defined and has a finite 3-volume, which is an important characteristic.The formula for accurately determining the volume of spherical 3-spaces is given by V k=1 = 2π 2 /|Γ|, where |Γ| represents the order of the group Γ.It is worth noting that when the 3-manifold is topologically complex, the order of the group Γ tends to be large, resulting in a relatively small corresponding volume.Additionally, it is important to mention that there is no lower bound on the volume since the group Γ can have an arbitrarily large number of elements according to [29].
Ongoing research and exploration within the scientific community are still being carried out to determine the allowed topologies of H 3 /Γ.However, Thurston and Jørgensen have made significant contributions in this field by demonstrating the existence of the smallest volume for compact orientable hyperbolic 3-manifolds, as discussed in [97].This achievement can be attained through a finite number of manifolds, with the prime candidate being the Weeks-Matveev-Fomenko manifold.Its volume is approximately 0.942707.The research conducted by Thurston and Jørgensen has greatly enhanced our understanding of the intricate and fascinating properties of hyperbolic 3-manifolds and their volumes.
The Lagrangian that corresponds to the aforementioned action (25) is As a result, the corresponding ADM Hamiltonian is given by where is the conjugate momenta of the scale factor.The ability to select the lapse function provides us with gauge freedom, resulting in the super-Hamiltonian constraint The minisuperspace of the FLRW cosmology, corresponding to the Hamiltonian constraint (30), is unidimensional with 3  and [H] = 1/(length), one can easily verify that D α (m P ) = m 3(2−α) P .Therefore, the fractional extension of the ADM Hamiltonian (28) and Hamiltonian constraints (30) are given by where l P = √ G is the Planck length.In the comoving gauge N = 1, using the Hamilton's equation, ȧ = ∂H (α)  ADM /∂Π, and rescaling the cosmic time by t → α lead us to the following fractional extension of the Friedmann equation where ρ is given by relation (24), and the parameter D is related to the Lévy fractional parameter α by Substituting α = 2 (or equivalently D = 2) into Equation (34), we obtain the standard counterpart, i.e., (48).We should note that there are several variants of the fractional Friedmann Equation in the literature; see, for instance, [43,44].
Similar to the standard model of cosmology, one can define the fractional extension of the density parameters It should be noted that the density parameters mentioned above follow the standard definitions in ΛCDM when D = 2.As a result, the fractional Friedmann (34) takes the following standard form Ω(m) + Ω(rad) which is also known as the closure relation.The Raychaudhuri equation, also known as the second Friedmann equation, can be derived by taking the time derivative of the Friedmann Equation ( 34) and combining it with the continuity equation.In fractional cosmology, the energy-momentum tensor of a perfect fluid is subject to the same covariant conservation equation as in standard cosmology.This equation can be expressed as follows: where These two lead us to the fractional Raychaudhuri equation In addition, the fractional Friedmann Equation (34), utilizing the continuity Equation (38), can be rewritten as where Ω(i) 0 and Ω(k) 0 denote the density parameters at the present epoch, (a = 1).In the same method, one can obtain the deceleration parameter q.By taking a time derivative of the Friedmann equation (34) and employing the continuity equation and fractional density parameters (36), we find where, in the second equality, we used the definition of the redshift 1 + z = 1/a.To observe the overall impact of the Lévy fractional parameter on cosmic evolution, it is instructive to examine a simple flat (k = 0) model of the universe, consisting of cosmic dust and the cosmological constant.The behavior of the q in relation to the redshift z is depicted in Figures 1 and 2. In Figure 1, we maintain a constant value of D = 2.2 and examine the impact of Ω(m) 0 on the variable q in terms of z.Our observation reveals that q does not respond significantly to the current values of Ω(m) 0 .However, in Figure 2, it becomes apparent that q is highly sensitive to the Lévy fractional parameter.By analyzing Figure 2, we can deduce that as α increases, the transition from a deceleration phase to an acceleration phase occurs at lower redshifts.

Fractional Cosmology Versus Data
It is important to mention that the analysis in this section is limited to the background level and does not take any disturbances into account.In the following section, we compare our theoretical findings with observations.We examine four significant datasets, namely SNIa, CMB, BAO, BBN, and OHD.Some articles suggest that OHD, which is obtained based on redshift, can be compared with the SNIa dataset.If one is interested, one can refer to reference [98] and other relevant sources.The largest combined sample of SNIa, named the Pantheon Sample, consists of 1048 SNIa in the redshift range of [0.01, 2.3], as presented by Scolnic et al. [99].In this analysis, we use the binned Pantheon SNIa dataset, which consists of 40 data points within the 0.014 < z < 1.61 range.We perform an analytical marginalization of the nuisance parameter in the distance estimates M. In ref. [100], the introduction of χ 2 S N as in which the covariance matrix, denoted by S , is calculated using the binned Pantheon sample, taking into account both statistical and systematic errors.Additionally, V is a row vector consisting of unitary elements, while W i represents the difference between µ b,i and µ(z i ).The distance modulus, µ, is predicted theoretically using the luminosity distance d L , which is given by the following equation: Here, E(x) is the dimensionless Hubble parameter, and h is the Hubble constant.The normalization constants independent of the cosmological model can be eliminated for this particular dataset.Afterwards, we incorporate the data obtained by observing acoustic signatures within the large-scale clustering of galaxies.By utilizing the BAO data, it is possible to minimize the χ 2 BAO that is defined as [101] Since the SNIa and BAO data contain valuable information regarding the universe at lower redshifts, we incorporate the CMB shift information by considering the compressed CMB likelihood.This likelihood is discussed in Table 1 of ref. [102]) and focuses on the angular scale of the sound horizon at the last scattering, denoted as l a , as well as the baryon density parameter.By including this information, we aim to explore the entire history of expansion leading up to the last scattering surface.To achieve this, we refer to Equations ( 22)-( 33) in ref. [103].Additionally, the Big Bang Nucleosynthesis (BBN) offers a data point that helps constrain mainly Ω (b) 0 , as mentioned in the work by Serra et al. [104].
The χ 2 BBN is given by Finally, we additionally incorporate the data obtained from the observational Hubble parameter.In this particular investigation, we utilize the 31 data points derived from the recent and precise estimates of H(z) within the redshift range of 0.07 ≤ z ≤ 1.965.These particular data points are independent of the Baryon Acoustic Oscillation (BAO) data points and were presented initially in [105].In this case, we can write where σ i is the Gaussian error on the measured value of H i .
The ΛCDM model, also known as the Concordance Model, is widely accepted in cosmology as it accurately fits current observations.It has been extensively tested using various cosmological measurements.Our statistical analysis results are presented in Table 1.The fractional ΛCDM model, on the other hand, introduces an additional independent parameter (α or D) that represents the attribute of nonlocality due to the fractional derivative.The statistical analysis of this model is presented in Table 2. Figure 3 shows the evolution of the deceleration parameter q(z) against redshift for the ΛCDM and fractional ΛCDM models.As we see from this figure, the deceleration parameter of the standard model of cosmology is relatively greater than the fractional deceleration parameter in all redshifts.Figure 4 illustrates the 1σ and 2σ confidence regions that are obtained from fitting the fractional ΛCDM model to a comprehensive set of datasets, including Baryon Acoustic Oscillations (BAO), Supernovae (SN), Cosmic Microwave Background (CMB), Big Bang Nucleosynthesis (BBN), and Observational Hubble Data (OHD).This fitting process determines the best possible agreement between the model and the observed data.In order to further evaluate the accuracy of the fractional ΛCDM model, Table 2 was prepared to showcase the consistency between the derived values of Ωm and H 0 obtained from the BAO+SN+CMB+BBN+OHD datasets and the corresponding values reported by the Planck 2018 Collaboration.By comparing the two, we can gain valuable insights into the reliability of our fits.Additionally, Figure 5 was included to provide a visual representation of the performance of our fits as compared with the well-established ΛCDM model.It is worth mentioning that both the fractional ΛCDM model and the ΛCDM model exhibit agreement with the most up-to-date cosmological observations when it comes to their respective values of the deceleration parameter q 0 and the associated transition redshift z t = 0.72.This further solidifies the credibility of our findings.Observations of Type SNIa provide crucial data for understanding the expansion of the universe.These observations serve as primary evidence for the universe's accelerated expansion.To obtain the best possible results from teh SNIa data, we compare the observed distance modulus of SNIa detection with the theoretical value.In this regard, we utilized the Pantheon sample,   which is an updated dataset of SNIa with 1048 distance moduli (µ) at various redshifts within the range of 0.01 < z < 2.26.In Figure 6, we illustrate the Pantheon Survey as the standard Hubble diagram of SNIa with an absolute magnitude of M 0 = −19.37.This figure also compares our model with the ΛCDM model, showing significant similarities.The Hubble parameter plotted against the redshift z is depicted in the graph illustrated in Figure 7, utilizing the data extracted from Table 1.Upon perusing the figure, it becomes evident that our model demonstrates a remarkable degree of concurrence with the conventional model of cosmology, thereby substantiating its validity and reliability., baryon density parameter Ω (b) 0 , matter density parameter Ω (m) 0 , curvature density parameter Ω (k) 0 , and Hubble parameter H 0 when using data from the BAO+SN+CMB+BBN+OHD.

The Synchronicity Problem
Let us review the critical property of the ΛCDM model related to the synchronicity problem.The limiting case of D = 2 in the Friedmann Equation (40) gives us the Friedmann equation in the standard (ΛCDM) model of cosmology    2 and 1, respectively. where . The solutions of the aforementioned equation for a flat universe, when radiation, cosmic dust, or the cosmological constant dominate, are given by where, in the last equality, C is a constant of integration.It is easy to verify that we have Ht = 1/2 and Ht = 2/3 in the first two cases, respectively.In the context of power-law expansion, the ratio of the universe's age to the Hubble time is a constant.This is a general property.However, when the cosmological constant dominates in late-time cosmology, we have the equation tH = √ Λ/3t.This means that at the present epoch, concerning cosmological observations, we have t 0 ∝ √ 3/Λ.This equation, in fact, is an expression of the synchronicity problem of the universe's age.
Let us analyze the effect of the Lévy fractional parameter on the synchronicity problem of the age of the universe.We obtain explicit solutions for the scale factor in a flat universe with a single-component perfect fluid.Thus, let us consider the i-th component is dominated in the Friedmann equation (40).Regarding Equations ( 40) and ( 41), the scale factor and the deceleration parameter for a perfect fluid with the equation of state This demonstrates that within the radiation-dominated universe (ω = 1/3), the expansion rate is a(t) ∝ t 3 6−D .Furthermore, it is noteworthy that the nature of the scale factor in a matterdominated universe (ω = 0) remains unaffected by the Lévy fractional parameter and is given by a ∝ t 2 3 .On the other hand, for a de Sitter universe (ω = −1) and D 2, the above form of scale factor reduces to a(t) ∝ t 2 3(D−2) .Furthermore, according to Equation ( 51), it can be inferred that in the context of fractional cosmology, our universe has experienced a period of deceleration during the phases dominated by radiation and pressureless matter.To summarize, when considering the presence of radiation and pressureless matter, fractional cosmology is unable to elucidate the cosmic phase transition from deceleration to acceleration that occurred throughout the history of the universe unless the inclusion of dark energy (in our model, a cosmological constant) is taken into consideration.Also, the deceleration parameter (51) shows that we have an accelerated universe for D < 8/3 restricted values of D.
Solution (50) shows that in the fractional ΛCDM model, the universe follows a power-law accelerated expansion even in recent times, unlike the standard model, where the late-time acceleration is explained by exponential de Sitter expansion.In the fractional ΛCDM model, for the i-th component of cosmic fluid, we have Thus, for a de Sitter universe, we find Unlike the ΛCDM model, the present epoch in the history of the universe is not a particular moment; instead, the age is always proportional to the Hubble time at that moment.Figure 8 displays the evolution of Ht in relation to the scale factor for both the standard model of cosmology and its fractional extension.It is evident that at the current epoch, both models yield nearly identical values for H 0 t 0 .However, in the future, their behaviors will diverge.In the  and 1, respectively.
ΛCDM model, Ht increases for a > 1, whereas in the case of fractional ΛCDM, it remains constant, as we observe in the aforementioned single-component fluid.
The content matter of the universe has been predominantly composed of cosmic dust throughout most of its entire history.Regarding this fact and ignoring the radiation and curvature terms (cold matter plus the cosmological constant), one can easily show that Utilizing the above equation, we find that when 2 ≤ D < 3, the current age of the universe, in units of Hubble time, is restricted to 0.6667 < H 0 t 0 ≤ 0.9887.
Figure 9 illustrates how the universe's current age is restricted for different fractional parameter values.Based on inequality (55), we can observe that when the fractional parameter D approaches 2, the value of H 0 t 0 becomes almost equal to one.However, in the standard model of cosmology, where D = 2, Ht tends towards infinity in the distant future.For other fractional parameter values, Ht remains finite.Based on the information presented in Table 2, the fractional parameter has a fitted value of D = 2.007.For this value of D, we obtain In addition, in the upper limit of the fractional parameter, D → 3, we have lim a→∞ Ht = 0.6667.
The analysis conducted shows that in the fractional ΛCDM model, the value of Ht is always finite, and at the present epoch, its value is consistently close to unity.Therefore, it can be confidently concluded that the value of Ht is finite in the fractional ΛCDM model.It is worth mentioning that fractional cosmology may have a considerable impact on an inflationary phase.Utilizing a scalar field with an appropriate and dominant potential compared with the universe's initial energy density is the typical method employed to achieve inflation.This inflationary approach triggers a rapid expansion of the scale factor of the universe.In the initial theories of inflation, this expansion would manifest as exponential growth, resembling de Sitter space, as the scalar field gradually moves towards its global minimum.By considering that the potential of the scalar field primarily affects inflation, we can deduce, based on the fractional deceleration parameter, that during the early universe, acceleration was experienced when the value of D lay between 2 and 2.667.This specific range of D values establishes a power-law inflationary model, and within this model, the fractional parameter α assumes a significant role.If we assume that the inflationary epoch took place within a time interval of t i = 10 −36 s, and t f = 10 −33 s, we can calculate the number of e-foldings: We obtain D ≤ 2.072 or α ≥ 1.866 if we assume a minimum of 64 e-foldings for addressing common Big Bang issues.This conclusion is similar to what we discovered in 4.

Conclusions
The concept of synchronicity is perplexing because the dimensionless cosmic time, represented as H 0 t 0 , can vary significantly within the standard model of cosmology.This model explains the late-time accelerated expansion as a de Sitter universe.In this model, the current coincidence of cosmic time and Hubble time is only valid for the present moment.However, the age of the universe will not be as near to the Hubble time in the past or future, and tH might even be infinite in the future.
In this paper, we present a comprehensive demonstration indicating that in the fractional ΛCDM model, during the late times, when the cosmological constant is dominated, the universe exhibits a remarkable and intriguing phenomenon characterized by a power-law accelerated expansion.This stands in sharp contrast to the standard model, in which an exponential de Sitter expansion is obtained to explain the late-time acceleration.Therefore, the age of the universe is continuously linked to the Hubble time at any given moment, creating a proportional relationship between the two phenomena.This contrasts the traditional cosmology framework, where the present epoch in the vast expanse of the universe is perceived as a specific instant in time.
Our statistical analysis revealed that the fractional parameter proposed by Lévy, which characterizes the degree of nonlocality in a system, is determined to be precisely α = 1.986.We have also determined the age of the universe through the meticulous study of various cosmological phenomena and the use of sophisticated techniques.Our calculations show that the age of the universe is t 0 = 13.8196 billion years.Furthermore, by calculating the dimensionless age of the universe, we find that H 0 t 0 = 0.9858, and Ht = 95.238 for the distant future.These findings indicate that our model is a plausible and valid cosmological model.It successfully encompasses a crucial early stage of matter domination, essential for forming and emerging intricate structures within the universe.Furthermore, our model also incorporates a more recent acceleration phase that corresponds harmoniously with the observations made in the field of cosmology.
In conclusion, it is essential to emphasize that our fractional model, which serves as a modified gravity theory, is currently in its nascent stage and necessitates further exploration and development.By delving into numerous issues within modern cosmology, such as the Hubble tension, gravitational lensing, gravitational-wave propagation, and density perturbation generation, we not only facilitate the advancement of fractional models but also potentially uncover fascinating findings at both the quantum and classical levels.However, it is prudent to defer the aforementioned evaluations to future endeavors.

Figure 3 :
Figure 3: The evolution of the deceleration parameter q(z) is shown against redshift.The dotted line represents the ΛCDM model, while the dashed line represents the fractional ΛCDM model.

DFigure 4 :
Figure 4: The two-dimensional contours and one-dimensional marginalized distributions represent the fractional ΛCDM model parameters with 68% CL and 95% CL, respectively.

Figure 6 :
Figure 6: The distance modulus of supernovae (data points) is compared with the theoretically predicted distance modulus (red line) in a fractional ΛCDM model.

Figure 7 :
Figure 7: The graph displays the evolution of the Hubble parameter, H(z), in units of km • sec −1 • Mpc −1 versus redshift z, with error bars.The red line represents the dynamics of the Hubble parameter obtained from Equation (40).The values of the density parameters, fractional density parameters, and Hubble constant at the present epoch are the best-fit values from Tables2 and 1, respectively.

Figure 8 :
Figure 8: The evolution of Ht with the scale factor is depicted for D = 2 (dashed) and D = 2.3 (line).The best-fit values for the density parameters and fractional density parameters are substituted from Tables2 and 1, respectively.

Table 1 :
The best-fit parameters with 1σ and 2σ confidence levels (CLs) for the FLRW model.

Table 2 :
The best-fit parameters with 1σ and 2σ confidence levels (CLs) for fractional FLRW model.
Figure 9:The dependence of H 0 t 0 to the fractional parameter, D. The best-fit value of the fractional density parameter of dust is substituted from Table2.