1. 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,
, where
and
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
. 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
, 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
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
, describing the distance that light has traveled since the Big Bang and the gravitational horizon
. 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
, the addition of any amount of matter will inevitably lead to a deviation from the strictly linear evolution of the
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 [
18,
24,
25,
26,
27,
28,
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 cosmological models, as evidenced by the wealth of studies conducted [
45,
46,
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 variable-order 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
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
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,
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 [
83,
84].
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
Section 6, we summarize the results and offer additional discussions.
2. 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
, which can be expressed as
This replacement paves the way for the exploration of the space-fractional SE:
where
is the Lévy index, and
is, according to Laskin, a generalized “fractional quantum diffusion coefficient” of dimension
. As mentioned earlier,
is a generalization of the fractional Riesz derivative, known as the fractional Laplacian [
86]:
which, by means of Fourier transforms, relates
and
. For the special case when
,
reduces to
, 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,
represents the wave function of the initial state, while
represents the wave function of the final state. The kernel itself, denoted as
, is given by the following expression:
Here, it is important to note that
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
is equal to
, 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
, while the shift is assigned a value of zero,
. 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,
denotes the reduced DeWitt metric of
n-dimensional minisuperspace with signature
. Also,
, the potential, includes the Ricci scalar of the three-dimensional
hypersurfaces and the potential of matter fields. Generally, the minisuperspace coordinates
may include matter fields and three-metric components. Variation in the action with respect to
yields the Euler–Lagrange equations
where
are the connection coefficients (Christoffel symbols) determined from the metric of the minisuperspace
. Also, in a variation with respect to the lapse function,
N, one obtains the constraint
To ensure consistency, it is necessary for Equations (
10) and (
11) to be equal to the 00 and
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
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
, 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 d’Alembertian,
is the curvature of the minisuperspace metric
, 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
.
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
, and
represents a Fourier transformation. Therefore, the fractional equivalent of the WDW Equation (
15) can be expressed as [
24,
25,
26,
27,
29]
where
is the coefficient, depending on the Planck mass
, which reduces to unity at
limit, and the value of Lévy’s fractional parameter, denoted as
, is assigned the range of values
.
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
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
is the fractional ADM Hamiltonian defined by
In the following section, we derive the fractional
CDM cosmology using Equations (
19) and (
20).
3. 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 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 stands for the standard line element for , encompassing both the angular coordinates and . Lastly, the parameter k assumes values of 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 [
36,
91,
92,
93,
94,
95,
96], can be expressed as follows:
where
K denotes the trace of the extrinsic curvature of the spacelike submanifolds,
, of the space–time,
, in which
and
represents the manifold of the space–time with boundary
. In addition,
is the total energy density of the cosmic fluid with components
where
,
, and
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
.
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
(representing three-dimensional Euclidean space),
(representing a 3-sphere), or
(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
, 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
. 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
. Regarding dimensions of
and
, one can easily verify that
. Therefore, the fractional extension of the ADM Hamiltonian (
28) and Hamiltonian constraints (
30) are given by
where
is the Planck length.
In the comoving gauge
, using the Hamilton’s equation,
, and rescaling the cosmic time by
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
(or equivalently
) 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
by
It should be noted that the density parameters mentioned above follow the standard definitions in
CDM when
. As a result, the fractional Friedmann (
34) takes the following standard form
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
and
denote the density parameters at the present epoch,
.
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
.
To observe the overall impact of the Lévy fractional parameter on cosmic evolution, it is instructive to examine a simple flat (
) 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
Figure 1 and
Figure 2. In
Figure 1, we maintain a constant value of
and examine the impact of
on the variable
q in terms of
z. Our observation reveals that
q does not respond significantly to the current values of
. 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.
4. 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
, 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
range. We perform an analytical marginalization of the nuisance parameter in the distance estimates
M. In ref. [
100], the introduction of
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
represents the difference between
and
. The distance modulus,
, is predicted theoretically using the luminosity distance
, which is given by the following equation:
Here, 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
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
, 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
, as mentioned in the work by Serra et al. [
104].
The
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
within the redshift range of
. 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
is the Gaussian error on the measured value of
.
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
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
and
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
and
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
and the associated transition redshift
. 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
. In
Figure 6, we illustrate the Pantheon Survey as the standard Hubble diagram of SNIa with an absolute magnitude of
. 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.
5. The Synchronicity Problem
Let us review the critical property of the
CDM model related to the synchronicity problem. The limiting case of
in the Friedmann Equation (
40) gives us the Friedmann equation in the standard (
CDM) model of cosmology
where
, and
. 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
and
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
. This means that at the present epoch, concerning cosmological observations, we have
. 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
are
This demonstrates that within the radiation-dominated universe (
), the expansion rate is
. Furthermore, it is noteworthy that the nature of the scale factor in a matter-dominated universe (
) remains unaffected by the Lévy fractional parameter and is given by
. On the other hand, for a de Sitter universe (
) and
, the above form of scale factor reduces to
. 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
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
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
. However, in the future, their behaviors will diverge. In the
CDM model,
increases for
, 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
, the current age of the universe, in units of Hubble time, is restricted to
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
becomes almost equal to one. However, in the standard model of cosmology, where
,
tends towards infinity in the distant future. For other fractional parameter values,
remains finite. Based on the information presented in
Table 2, the fractional parameter has a fitted value of
. For this value of
D, we obtain
In addition, in the upper limit of the fractional parameter,
, we have
The analysis conducted shows that in the fractional CDM model, the value of 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 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
s, and
s, we can calculate the number of e-foldings:
We obtain
or
if we assume a minimum of 64 e-foldings for addressing common Big Bang issues. This conclusion is similar to what we discovered in
Section 4.
6. Conclusions
The concept of synchronicity is perplexing because the dimensionless cosmic time, represented as , 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 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 . 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 billion years. Furthermore, by calculating the dimensionless age of the universe, we find that , and 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.