Inflation and fractional quantum cosmology

The Wheeler--DeWitt equation for a flat and compact Friedmann--Lema\^{i}tre--Robertson--Walker cosmology at the pre-inflation epoch is studied in the contexts of the standard and fractional quantum cosmology. Working within the semiclassical regime and applying the WKB approximation, we show that some fascinating consequences are obtained for our simple fractional scenario that are completely different from their corresponding standard counterparts: (i) The conventional de Sitter behavior of the inflationary universe for constant potential is replaced by a power-law inflation. (ii) The non-locality of the Riesz's fractional derivative produces a power-law inflation that depends on the fractal dimension of the compact spatial section of space-time, independent of the energy scale of the inflaton.


Introduction
Inflation is widely accepted as a formal solution to cosmological problems such as horizon and flatness [1,2,3]. The conventional approach to controlling inflation is that a scalar field with an appropriate potential, such as the Coleman-Weinberg potential [4], dominates the energy density of the universe from the start. In early theories of inflation, this energy density promotes a rapid expansion for the scale factor of the universe [5,6]. This accelerated expansion would be exponential, similar to de Sitter space, or power law, according to which the scalar field gradually decreases to the global minimum of its potential.
In addition to the well-known standard models of inflation, recent scenarios using noncommutative cosmological models, Generalized Uncertainty Principle (GUP) and other modified gravity models have shown that some problems with the standard models (such as graceful exit and the Hubble constant problem) can be solved, see, for instance [7,8,9,10,11,12].
It should be noted that quantum cosmology is a suitable method to study the essential initial conditions for the emergence of an inflationary phase. Tryon proposed in 1973 that a closed universe [13] could emerge spontaneously as a quantum fluctuation. He realized that in a spatially closed universe, all conserved charges are zero. Consequently, no conservation law prevents such a universe from forming spontaneously. Furthermore, according to general relativity, at the instant of creation of such a universe, not only matter fields, but also space-time itself is created, and there has not been anything earlier than that. In fact, if we assume that the universe is spatially homogeneous, isotropic, compact, and simply compact, the only feasible choice to have a space without a boundary is a closed space. The volume of open and flat spaces that are simply connected is infinite. By removing the restriction of simple connectivity, more sophisticated yet finite volume spaces for all possible flat, open, and close homogeneous and isotropic spaces can be obtained (please see [14] and references therein). In this work, we assume that the cosmic manifold is spatially compact since we are trying to understand the universe as a whole and it is implicitly inconceivable that the universe has a spatial boundary. This implies that the three-dimensional spacelike hypersurfaces are compact (in mathematical language, a compact manifold without a boundary is called a closed manifold).
Intuitively, a universe formed by a quantum fluctuation should be exceedingly small, having a Planck length linear size. This was initially a major challenge since it was unclear how to create the enormous universe we live in from a tiny compact quantum cosmological model. With the emergence of inflationary theories in which the universe goes through a de Sitter phase of exponential growth, the dilemma has vanished. All scales in the universe are expanded by a gigantic factor exp(Ht) as a result of inflation, where H is the constant expansion rate and t is the duration of the inflationary phase or t α where α must be greater than one and ideally at least of the order of ten.
Given the centrality of inflation in our present understanding of cosmology, it is reasonable and vital to try to grasp its details within the context of fractional quantum cosmology. The history of fractional derivatives is as long as that of classical calculus. A fractional derivative is a generalization of the integer-order derivative. It originated in the letter regarding the meaning of half-order derivative from L' Hôpital to Leibnitz in 1695 and is a promising tool for explaining various phenomena. Various definitions of fractional derivatives exist in the literature, including Riesz, Riemann-Liouville, Caputo, Hadamard, Marchand, and Griinwald-Letnikov, among others [15,16]. Recent quantum gravity conclusions have provided an essential push for the increased use of fractional calculus in quantum theory. Various approaches to quantum gravity, such as asymptotically safe quantum gravity [17,18,19,20], causal dynamical triangulations [21,22,23,24], loop quantum gravity, and spin foams [25,26,27], Hořava-Lifshitz gravity [28,29], non-local quantum gravity [30,31,32], and others, all lead to the same result: the dimension of space-time changes with scale. As a result of the abnormal scaling of the space-time dimension, all known theories of quantum gravity are multiscale. There are numerous applications for fractional calculus and active research in this field. Fractional quantum mechanics has been used to model fractional space-time in gravity and cosmology [33,34,35,36,37,38,39,40], and the fractional quantum field theory [41,42,43].
The main objective of this work is to reexamine the standard inflation in the context of fractional quantum cosmology (FQC) as a novel underlying scenario for studying the early state of the universe. Concretely, this study aims to employ the FQC with a minimally coupled scalar field and a compact 3-space as a probe model. We are aware of the limited breadth. However, fascinating aspects can be found. In the next section, we present a review of the model in the context of the canonical quantum cosmology and obtain a semi-classical wavefunction. Section 3 focuses on the fractional counterpart of the model mentioned above. Then, by analyzing the implications of fractional quantum cosmology, we examine why they may be more intriguing than those reached in Section 2.

Wheeler-DeWitt Equation in Slow Roll Regime
Let us consider the Friedmann-Lemaître-Robertson-Walker (FLRW) metric where N(t) is the lapse function, a(t) is the scale factor, and we will assume that the space-time manifold, M, is a spatially compact and globally hyperbolic Lorentzian manifold. We claimed that M is spatially closed since we want to investigate the cosmological model as a whole. It is basically impossible to imagine that the universe has a spatial boundary. Consequently, the spacelike sections of the space-time manifold, (Σ, h), have a finite volume with the normalized curvature index K taking values 0, ±1. Different ranges of variation for the coordinates (r, θ, ϕ) determine different topologies. Let us rewrite the metric as Then, we can define the volume of the spacelike hypersurfaces as Therefore, scalar curvature of the spacelike hypersurface (Σ, h) of the above metric is constant and equal to 6K/a 2 . One can show that any compact Riemannian three-manifold with constant curvature is homeomorphic toΣ/Γ, where Γ denotes the group of the covering transformations andΣ is the universal covering space, which is either R 3 (3-dimensional Euclidean space), S 3 (3-sphere), or H 3 (3-dimensional hyperbolic space) regarding the sign of K (K = 0, K = 1, or K = −1, respectively). Therefore, in three dimensions, an oriented compact and without boundary space is now a polyhedron whose faces are identified in pairs. For an overview, see Ref. [14]. As we would like to investigate a new scenario regarding an inflationary universe, let us consider the same action corresponding to standard inflationary theories: where φ is a scalar field minimally coupled to the Ricci scalar R, and we used the units where c = 1 = . As during inflation only, the inflaton φ dominates the dynamics, the action does not contain any other matter fields. By substituting the Ricci scalar corresponding to the metric (1), i.e., into action functional (4), we obtain Arnowitt-Deser-Misner (ADM) action of the model where an overdot denotes differentiation with respect to t. 3 The conjugate momenta associated with a and φ, respectively, are given by by which we can easily obtain the canonical ADM Hamiltonian: The above ADM Hamiltonian leads to the super-Hamiltonian constraint To obtain the Wheeler-DeWitt (WDW) equation, we apply the quantization map where the ordering parameter p is considered for the factor-ordering ambiguity [44]. By the action of the super-Hamiltonian operatorĤ on a wave function Ψ(a, φ), we can therefore write the WDW equation: where m P = 1/ √ G is the Planck mass. As mentioned, in this work we are interested in investigating a fractional counterpart of an inflationary model established during a 'slow roll' regime. For such a special but familiar scenario, we can ignore the φ-dependence of the wave function. Furthermore, in such a regime, we can assume that the scalar potential plays the role of an effective cosmological constant [45]. More precisely, we define Λ ≡ (8πV 0 )/m 2 P = constant. Additionally, let us define L ≡ √ 3/Λ, and assume p = 0. In the following, we are interested in semi-classical approximation, and in this case the factor ordering does not play a crucial role. According to Linde [46], a compact flat universe, such as the toroidal universe (K = 0, Σ = T 3 = S 1 × S 1 × S 1 ), is the simplest case to solve the problem of initial conditions for low-scale inflation. While the quantum creation of a closed or infinite open inflationary universe is exponentially inhibited, the toroidal universe is not.
Consequently, substituting p = 0 and K = 0 into WDW Equation (11), we obtain A solution for (12) in terms of the Bessel function is: where N is a normalization constant. The asymptotic form of the Bessel function, for large values of the scale factor, i.e., yields the asymptotic form of the wave function This wavefunction collects both the expanding exp(−i phases of the cosmological model simultaneously. We now discuss the classical limit of the wavefunction. First, we substitute the wavefunction written in form (Equation (16) is the single component WKB solutions to the WDW equation in the oscillatory regime [47]. To study the oscillatory behavior of the particular differential equations, see [48].) into the WDW Equation (12). Straightforward differentiation plus condition (14) (which is equivalent to WKB) (The semi-classical approximation and also the WKB wave function are important concepts in the context of quantum cosmology. It has been believed that the semi-classical approximation is employed to avoid the ambiguities arising from operator ordering issues in the Wheeler-DeWitt equation, as well as problems in the path integral formulation of the wave function. The WKB wave function is best suited to approximate the wave function in the semiclassical regime [47]), i.e., d 2 S /da 2 ≪ (dS /da) 2 ) lead to We recognize that (17) is the Hamilton-Jacobi equation in classical cosmology, which is equivalent to the Hamiltonian constraint. As we know, the momentum in the classical Hamilton-Jacobi theory is given by Π a = −dS /da. Therefore, with this definition and the corresponding relation in (7), we obtainȧ/a = 1/L or a(t) = a(t 0 ) exp((t − t 0 )/L). Thus, we have shown that the scale factor grows exponentially during the inflationary epoch, t i ≤ t ≤ t f , which is assumed to occur in a classical region, This solution implies the expansion rate during the de Sitter expansion. The so-called number of e-folds N e is often introduced as N e = ln(a(t f )/a(t i )). If we assume that the inflation occurred during a very small time interval t f − t i = 10 −37 s = 10 8 t P , the minimum number of e-foldings necessary for resolving the standard problems of the Big Bang is N e ≃ 60. Substituting these values into the above de Sitter solution, we obtain L ≃ 10 7 l p .

Fractional Quantum Cosmology for a Slow Roll Regime
In this section, we investigate our model in the context of the FQC. In this context, we first provide a brief overview of the FQC framework. In order to construct a fractional WDW equation, it is best to obtain the fractional Schrödinger equation (SE). 5 The classical Hamiltonian of a particle with mass m, space coordinates r, and momentum p is H = p 2 2m + V(r). Focusing on the quadratic dependence of p in the above equation might inspire to adopt a different function for the kinematic term so that it does not conflict with basic principles of mechanics. An acceptable approach from theoretical physics for this motivation is the Feynman path integral approach to quantum mechanics [49], first developed by Laskin [50]. Non-relativistic quantum mechanics was reconstructed by Feynmann and Hibbs as a path integral over the Brownian paths. Moreover, the space fractional quantum mechanics formulated by Laskin applying new fractional path integral on the basis of the Lévy flight [50,51]. In addition, Naber obtained the time fractional SE on the basis of fractional Brownian motion [52], see also [53]. Inspired by the above mentioned works, the generalized fractional SE was obtained where both the space and time derivatives were replaced by fractional counterparts.
In what follows, let us present a brief overview of the fractional quantum dynamics, which is generated by the Hamiltonian function H α (p, r) := D α |p| α + V(r) (for a detailed review, see [54]), where D α is a scale coefficient, and the Lévy's fractional parameter α relevant to the concept of Lévy path [55] is set to 1 < α ≤ 2. Then, in order to construct the corresponding quantum Hamiltonian of the system, we should apply the standard canonical quantization procedure, i.e., (r, p) → (r, −i ∇). Hence, we obtain H α = D α (− 2 ∆) α/2 + V(r). Now the usual SE, can be extended to its fractional counterpart by replacing the ordinary Laplace operator ∆ with the Riesz fraction (− 2 ∆) α/2 : Corresponding to a three-dimensional Euclidean space, the Riesz fractional Laplacian (− 2 ∆) α/2 in terms of the Fourier transformation is written as [56] (− 2 ∆) α/2 ψ(r, t) = F −1 |p| α F ψ(r, t) where p = p 2 1 + p 2 2 + p 2 3 . Finally, the fractional SE is proposed as Now, let us focus on our model. As for a homogeneous cosmological model, in addition to the supermomentum constraint for H i , the shift function N i , without neglecting any of Einstein equations, can be set to zero, the superspace WDW equation therefore reduces to the minisuperspace WDW equation: In Equation (23), is the d'Alembertian operator, q α and U(q ν ) are coordinates of N-dimensional minisuperspace and potential, respectively; f αβ stands for the corresponding minisuperspace metric with signature (−, +, +, . . . , +).
Regarding obtaining the fractional counterpart of the WDW Equation (23), the usual d'Alembertian operator should be replaced by [56,57,36,58] where |p| = −p 2 0 + p i p i , i = 1, 2, . . . , N, and F denotes Fourier transformation. Hence, the fractional counterpart of the WDW Equation (23) will be [14,35,34,59,33] As seen, the WDW Equation (12) is actually a one-dimensional SE with zero energy. Hence the final step of our procedure here is to obtain a particular fractional WDW equation similar to that proposed above for the SE case. Consequently, we can replace the ordinary derivative with the fractional Riesz derivative [14,35,34,59,33]: Therefore, we obtain the fractional version of the WDW Equation (12) for our herein model as Note that in a particular case where α = 2, Equation (27) reduces to its standard counterpart (12).
To obtain the semi-classical solution of this fractional WDW equation, we use the alternative equivalent representation of the Riesz fractional derivative (or fractional Laplacian) [60], given by where c 1,α = α2 α−1 √ π Γ((1+α)/2) Γ((2−α)/2) . Inserting the exponential form of the wavefunction, (16), into the above form of the fractional Laplacian, expanding Ψ(a ± v) in Taylor series about a, and using semi-classical approximation, we obtain where S ′ = dS /da. Therefore, in the classical regime, the fractional WDW equation yields the following fractional Hamiltonian constraint The corresponding fractional ADM Hamiltonian, counterpart to (8), is Hence, the WKB fractional wavefunction of the fractional WDW Equation (27), for expanding universe, disregarding the preexponential factor, is The Hamilton's equations,ȧ together with the Hamiltonian constraint (30), in the comoving frame, N = 1, lead to the following fractional solution: where according to Refs. [59,33], we defined D ≡ 2/α + 1.
In what follows, let us share some fascinating results for the fractional evolution of the universe during the fractional inflation epoch.
Inflation is frequently achieved by having a scalar field with a suitable potential that dominates the early energy density of the universe. As a result of this energy density, the scale factor of the universe expands rapidly. In the original theories of inflation, this expansion would be exponential and would resemble de Sitter space as the scalar field steadily rolls down to its global minimum. Therefore, during the inflation, the deceleration parameter is q = −1. However, in fractional modification of the same model, we obtain where we used (34). According to Equation (35), the early universe is accelerated if 2 < D < 2.5, and for these values of D, we have a power law inflationary model. Hence, the Lévy's fractional parameter α plays a crucial role in this model.
Assuming that the inflationary epoch occurred during a time interval t f − t i = 10 −37 s = 10 8 t P after the Planck's time, we obtain the number of e-foldings as where we used the fractional form of scale factor (34). For instance, assuming D = 2.019 we obtain N e ≃ 60, that is the minimum number of e-foldings necessary for resolving the standard problems of the Big Bang cosmology.

Conclusions and Discussions
Our simple fractional quantum cosmological model with characteristic features effectively influenced the early inflationary phase of the universe. Observations of CMB fluctuations and the distribution of matter in the universe constrain the properties of the inflaton field. A heavy 8 scalar field with a mass of 10 13 GeV, close to the grand unified theory scale, is preferred by wellknown arguments even though the mass of such an inflaton and its interaction with matter fields are unknown. This field is widely cited as evidence that there are new physics at the intersection of the electroweak and Planck scales. Our results imply that power-law inflation within fractional quantum cosmology may emerge from a new angle by accounting for quantum gravity effects via the fractional derivative. In this paper, by considering a FLRW universe with all three curvature indices, a canonical kinetic energy associated with a scalar field minimally coupled to gravity, and a scalar potential, we obtained the corresponding Hamiltonian. Then, considering a non-zero ordering parameter, we constructed the standard WDW equation. To explore the goals of our paper, we focused on studying the early state of the universe. With this in mind, we confined ourselves to a semiclassical case, where the slow roll regime and a vanishing ordering parameter have been assumed for a compact flat universe, which represents the simplest case study to the problem of initial conditions [46].
The above assumptions motivated us to investigate our main objectives in two different contexts as follows.
Standard quantum cosmology: For this case, the subsequent steps have been followed.
• According to the arguments presented above, we applied the slow roll regime, and substituted p = 0 and K = 0 into generalized WDW Equation (11).
• We then showed that the solution of the simplified differential Equation (12) can be expressed by the Bessel function.
• The model was then examined according to the arguments set out previously; please see the expressions below Equation (15). Concretely, the model was studied assuming the WKB condition. For such a semiclassical regime, we have shown that the scale factor of the universe evolves exponentially at early times, which corresponds to the accelerated de Sitter expansion.
• We have additionally shown that L ≡ √ 3/Λ ought to take a value 10 7 l p to obtain a sufficient e-folding number required to solve the problems of the standard cosmology.
Fractional quantum cosmology: The fractional extensions of the above-referred outcomes have been the principle goal of this paper. In what follows, allow us to summarize the relevant consequences and analyze them.
• We gave a brief overview of fractional quantum mechanics and then explained how a suitable model is constructed for the fractional quantum cosmology. Based on the FQC procedure, we derived the fractional WDW equation and fractional ADM Hamiltonian from those retrieved for the corresponding standard model.
• To solve the fractional WDW equation, we applied the WKB approximation again and set N = 1.
• Despite the standard case, we have shown that the scale factor of the universe takes the form of a power-law function of time. It is seen that it not only depends on the Lévy's fractional parameter α, but also, as can be seen from the relation (34), on the volume of the compact 3-space Σ, V 0 .
Based on the points summarized above, let us reiterate the main differences in the consequences of the standard quantum cosmology and FQC for our simple model. 1. It is worth mentioning that the standard model and its fractional counterpart yield completely different results for the evolution of the scale factor: for the standard case, not only does the scale factor a(t) accelerate with the exponential-law of the time, but also such evolution is only affected by matter (i.e., the constant potential; cf. Equation (18)), and therefore the number of e-folding depends only on the constant L. In the fractional case, but the scale factor is a power-law function of time, which is completely independent of the matter; but instead, only the Lévy's fractional parameter α and the volume of the compact 3-space Σ, V 0 , determine how the scale factor evolves (cf. relation (34)). This gives us another dependency for the number of e-folding. Such consequences are interpreted as distinguishing features of the fractional quantum gravity. 2. As mentioned, unlike the usual de Sitter solution (18), the evolution of the scale factor associated with the fractional case, i.e., (34), depends on the volume of the compact 3space Σ, V 0 . More concretely, the evolution of the scale factor depends on the global geometry that affects the topology of the universe as a whole. The global geometry of the universe, i.e., the spatial curvature and topology and consequently the shape of the entire universe, are not determined by the Einstein gravitational field equations since they are differential equations that only determine the local features of space-time. Therefore, the 3-volume dependence of the scale factor is entirely a fractional quantum gravity effect. As we see in (34), the value of the scale factor at early time t = L is proportional to the inverse root square of V 0 , such that the smaller the value of the 3-volume, the bigger the initial scale factor. Interestingly, the authors of Ref. [61] (for a recent analysis, see [62,63,64]) have shown that the WMAP 3 year data are consistent with the possibility that we live in a 'small universe' shaped like a flat 3-torus T 3 whose fundamental domain is a cube with side length l = 4 corresponding to a volume of 5 × 10 3 Gpc 3 . It turns out that the torus model describes the data much better than the best-fit ΛCDM model, as it shows the suppression of the CMB anisotropy at large scales first observed by COBE. 3. As seen in Equation (28), the Riesz fractional derivative is a non-local operator that can explain non-local processes in the minisuperspace. In fact, non-locality is a general behavior of all fractional derivatives and integrals and hence they describe processes with non-locality in time (memory) and space (large jumps) [65]. The fractional derivative has lagged far behind the integer-order calculus due to the unclear physical meaning, which is a major obstacle. In 1974, it was proposed as an open problem to ask "What are the physical interpretations of fractional calculus?" Using analogous reasoning, a physical explanation for a fractional time derivative was put forth in 2002 [66], although no experiments were performed to support the new time scale. In Ref. [67], the authors have shown that the time fractional order is a memory index. Time disappears in quantum gravity and cosmology, and the space-fractional derivative may play a crucial role. In Refs. [59,33], the authors show that Lévy's fractional parameter α represents the fractal dimension (denoted by D in Equations (34), (35) and (36)) of the black hole horizon or the cosmological horizon. In addition, as we saw in (34), as a result of big jumps in minisuperspace, the initial value of emerged classical scale factor depends on the global geometry (or topology) of the universe. Furthermore, Equation (36) shows that the acceleration of the universe and e-folding of the inflation epoch are direct consequences of its topology. 10