Next Article in Journal
Temporal Fractal Nature of the Time-Fractional SPIDEs and Their Gradient
Next Article in Special Issue
Estimated Age of the Universe in Fractional Cosmology
Previous Article in Journal
Fixed-Time Distributed Time-Varying Optimization for Nonlinear Fractional-Order Multiagent Systems with Unbalanced Digraphs
Previous Article in Special Issue
Exact Solutions and Cosmological Constraints in Fractional Cosmology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Anisotropic Fractional Cosmology: K-Essence Theory

by
José Socorro
1,*,†,
J. Juan Rosales
2,† and
Leonel Toledo-Sesma
3,†
1
Departamento de Física, División de Ciencias e Ingenierías, Universidad de Guanajuato-Campus León, León 37150, Mexico
2
Departamento de Ingeniería Eléctrica, División de Ingenierías Campus Irapuato-Salamanca, Universidad de Guanajuato Carretera Salamanca-Valle de Santiago, km. 3.5 + 1.8 km, Comunidad de Palo Blanco, Salamanca 36885, Mexico
3
Unidad Profesional Interdisciplinaria de Ingeniería Campus Hidalgo, Instituto Politécnico Nacional, Carretera Pachuca—Actopan Kilómetro 1 + 500, San Agustín Tlaxiaca 42162, Mexico
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Fractal Fract. 2023, 7(11), 814; https://doi.org/10.3390/fractalfract7110814
Submission received: 15 September 2023 / Revised: 1 November 2023 / Accepted: 3 November 2023 / Published: 9 November 2023
(This article belongs to the Special Issue Fractional Gravity/Cosmology in Classical and Quantum Regimes)

Abstract

:
In the particular configuration of the scalar field k-essence in the Wheeler–DeWitt quantum equation, for some age in the Bianchi type I anisotropic cosmological model, a fractional differential equation for the scalar field arises naturally. The order of the fractional differential equation is β = 2 α 2 α 1 . This fractional equation belongs to different intervals depending on the value of the barotropic parameter; when ω X [ 0 , 1 ] , the order belongs to the interval 1 β 2 , and when ω X [ 1 , 0 ) , the order belongs to the interval 0 < β 1 . In the quantum scheme, we introduce the factor ordering problem in the variables ( Ω , ϕ ) and its corresponding momenta ( Π Ω , Π ϕ ) , obtaining a linear fractional differential equation with variable coefficients in the scalar field equation, then the solution is found using a fractional power series expansion. The corresponding quantum solutions are also given. We found the classical solution in the usual gauge N obtained in the Hamiltonian formalism and without a gauge. In the last case, the general solution is presented in a transformed time T ( τ ) ; however, in the dust era we found a closed solution in the gauge time τ .

1. Introduction

Fractional cosmology is a new line of research born approximately twenty years ago based on fractional calculus (FC). The FC is a non-local natural generalization to the arbitrary order of derivatives and integrals. Non-local effects occur in space and time. In the time domain, a non-local description becomes manifest as a memory effect, and in the space domain, it manifests as non-homogeneous similarity structures [1,2,3]. During the last decades, FC has been the subject of intense theoretical and applied research in almost all areas of the sciences and engineering from the point of view of classical and quantum systems [4,5,6,7,8,9,10,11,12,13,14]; recently, new studies on FC have been made [15,16,17,18]. This is because FC describes complex physical systems more accurately, and at the same time investigates more about simple dynamical systems [19,20]. General relativity is not an exception, in [21,22,23,24,25,26,27,28,29,30,31] the importance of FC and its potential applications in cosmology was introduced. In [32], the FRW universe was presented in the context of the variational principle of fractional action. In this new cosmological formulation, the accelerated expansion of the universe can be attributed to the fractional dissipative force without the need to introduce any kind of matter or scalar fields; similar results are obtained in [33,34]. In [22], the concept of fractional action cosmology was applied to massive gravity, where fractional graviton masses are introduced.
Unlike the previously described formalism to obtain fractional cosmology, in [35] it is mentioned that by quantifying different epochs of the k-essence theory, a fractional Wheeler–DeWitt equation in the scalar field component is naturally obtained. Recently, such an equation was solved for some epochs in the FRW model and communicated in [36]. In this work, we present the continuation of our previous investigation. In this case, we will analyze the Bianchi type I, which is the anisotropic generalization of the flat FRW cosmological model. In the quantum scheme, we introduce the factor-ordering problem in the variables ( Ω , ϕ ) and its corresponding momenta ( Π Ω , Π ϕ ), obtaining a fractional differential equation with variable coefficients in the scalar field equation. The solution is found using a fractional series expansion [37,38], generalizing our previous work [36].
This paper is organized in the following way: in Section 1, we give a brief review of fractional calculus and the main ideas of the k-essence formalism; in Section 2, we construct the Lagrangian and Hamiltonian densities for the anisotropic Bianchi type I cosmological model, considering a barotropic perfect fluid for the scale field in the variable X. We found the classical solution in the usual gauge N obtained in the Hamiltonian formalism, and without a gauge. In the last case, the general solution is presented in a transformed time T ( τ ) ; however, in the dust era we found a closed solution in the gauge time τ ; in Section 3, the quantization of the model for any era in our universe is performed and we present particular scenarios too. In this section, we introduce the factor ordering in both variables; finally, in Section 4, the conclusions are given.

2. Brief Review on Fractional Calculus and K-Essence Theory

2.1. Brief Review on Fractional Calculus

In the theory of fractional calculus, there are some definitions of fractional derivatives; Riemann–Liouville, Caputo, Caputo–Fabrizio, Atangana–Baleanu, to name a few, each with its advantages and disadvantages [39,40,41]. In this work, we use the Caputo fractional derivative of order γ , defined by using the Riemann–Liouville fractional integral [1]
I γ f ( t ) = 1 Γ ( γ ) 0 t f ( τ ) ( t τ ) 1 γ d τ , γ > 0 ,
recovering the ordinary integral when γ 1 . The Caputo fractional derivative of order γ 0 of a function f ( t ) , then, is defined as the fractional-order integral (1) of the integer-order derivative
0 C D t γ f ( t ) = I ( n γ ) 0 D t n f ( t ) = 1 Γ ( n γ ) 0 t f ( n ) ( τ ) ( t τ ) γ n + 1 d τ ,
with n 1 < γ n N = 1 , 2 , , γ R is the order of the fractional derivative, f ( n ) are the ordinary integer derivatives, and Γ ( x ) = 0 e t t x 1 d t is the gamma function. The Caputo derivative satisfies the following relations
0 C D γ [ f ( t ) + g ( t ) ] = 0 C D γ f ( t ) + 0 C D γ g ( t ) .
0 C D γ c = 0 , where c is a constant .
The Laplace transform of the function f ( t ) defined in the ordinary case is given by
L [ f ( t ) ] = 0 f ( t ) e s t d t F ( s ) ,
then, the Laplace transform of the Caputo fractional derivative (2) has the form
L [ 0 C D γ f ( t ) ] = s γ F ( s ) k = 0 n 1 s γ k 1 f ( k ) ( 0 ) ,
where f ( k ) is the ordinary derivative. Another definition which will be used is the Mittag–Leffler function [42,43,44],
E χ , σ ( z ) = n = 0 z n Γ ( n χ + σ ) ( χ , σ > 0 ) ,
for σ = 1 , we have a one-parameter Mittag–Leffler function:
E χ ( z ) = E χ , 1 ( z ) = n = 0 z n Γ ( n χ + 1 ) ( χ > 0 ) .
Other special cases are given in [43,44]
E 1 ( ± z ) = e ± z , E 2 ( z ) = cosh z , E 2 , 1 ( z 2 ) = cos z , E 2 , 2 ( z 2 ) = sinh z z , E 2 , 2 ( z 2 ) = sin z z .
A Laplace transform (5) of the Mittage–Leffler function is given by the formula
0 e s t t χ m + σ 1 E χ , σ ( m ) ( ± a t χ ) d t = m ! s χ σ ( s χ a ) m + 1 .
Consequently, the inverse Laplace transform is
L 1 m ! s χ σ ( s χ a ) m + 1 = t χ m + σ 1 E χ , σ ( m ) ( ± a t χ ) .
This expression will be very useful to obtain analytical solutions of fractional differential equations using the Laplace transform.

2.2. K-Essence Fractional in the Bianchi Type I Scenario

One of the fundamental problems of cosmology is to find an explanation consistent with experiments for the accelerated expansion of the universe. Many proposals to tackle this task suggest modifying the general relativity theory. A recent proposal suggests unifying the description of dark matter, dark energy, and inflation, employing a scalar field with a nonstandard kinetic term, known as k-essence theory. Usually, the action of the k-essence models [45,46,47,48,49,50] can be written as
S = d 4 x g 1 2 R + f ( ϕ ) G ( X ) + L m a t t e r ,
with g being the determinant of the metric, R the scalar curvature, f ( ϕ ) an arbitrary function of the dimensionless scalar field ϕ , X = 1 2 g μ ν μ ϕ ν ϕ the canonical kinetic energy, and L m a t t e r is the corresponding Lagrangian density of ordinary matter. So, performing the variation of the action (12) with respect to the metric g μ ν and X, the field equations are obtained:
  G μ ν f ( ϕ ) G X μ ϕ ν ϕ + G g μ ν = T μ ν ,
f ( ϕ ) G X μ μ ϕ + G X X μ X μ ϕ + d f d ϕ G 2 X G X = 0 ,
where we have assumed that 8 π G = 1 and a subscript X denotes differentiation with respect to X. K-essence was originally proposed as a model for inflation; and then, as a model for dark energy, along with explorations of unifying dark energy and dark matter [51,52].
The last set of field Equations (13) and (14) are the results of considering the scalar field X ( ϕ ) as part of the matter content, i.e., L X , ϕ = f ( ϕ ) G ( X ) , with the corresponding energy–momentum tensor
T μ ν ( ϕ ) = f ( ϕ ) G X μ ϕ ν ϕ + G ( X ) g μ ν .
Also, considering the energy–momentum tensor of a barotropic perfect fluid,
T μ ν ( ϕ ) = ( ρ ϕ + P ϕ ) u μ u ν + P ϕ g μ ν ,
with u μ being the four-velocity satisfying the relation u μ u μ = 1 , ρ ϕ the energy density, and P ϕ the pressure of the fluid. To simplify, we are going to consider a comoving perfect fluid, whose pressure and energy density corresponding to the energy–momentum tensor of the field X are
P ϕ ( X ) = f ( ϕ ) G , ρ ϕ ( X ) = f ( ϕ ) 2 X G X G ,
thus the barotropic parameter ω X = P ϕ ( X ) ρ ϕ ( X ) for the equivalent fluid is
ω X = G 2 X G X G .
Notice that the case of a constant barotropic index ω X (with the exception ω X = 0 ) can be obtained using the G function
G = X 1 + ω X 2 ω X .
At this point we can choose
G = X α , α = 1 + ω X 2 ω X ω X = 1 2 α 1 .
With this, we can write the states in the evolution of the universe summarized in the Table 1.
We are interested in the four-dimensional fractional cosmology in the scenario of k-essence within the anisotropic background, precisely, the Bianchi type I, whose metric has the line element g α β , which can be read as
d s 2 = N 2 ( t ) d t 2 + A 2 ( t ) d x 2 + B 2 ( t ) d y 2 + C 2 ( t ) d z 2 ,
where N ( t ) is the lapse function, and the functions A ( t ) , B ( t ) , and C ( t ) are the corresponding scale factors in the ( x , y , z ) directions, respectively. Moreover, in Misner’s parametrization, the radii for this anisotropic background have the explicit forms
A = e Ω + β + + 3 β , B = e Ω + β + 3 β , C = e Ω 2 β + ,
where the functions in the radii are dependent on time, Ω = Ω ( t ) , and β ± = β ± ( t ) . In this point, we notice that the line element (21) in the time d τ = N d t reads as
d s 2 = d τ 2 + e 2 ( Ω ( τ ) + β + ( τ ) + 3 β ( τ ) ) d x 2 + e 2 ( Ω ( τ ) + β + ( τ ) 3 β ( τ ) ) d y 2 + e 2 Ω ( τ ) 4 β + ( τ ) d z 2 ,
and employing the form of the functional G = X α , and the following quantities:
= d d τ = d N d t , g τ τ = g τ τ = 1 , G X = α X α 1 , G X X = α ( α 1 ) X α 2 , ϕ ; μ , μ = g μ ν ϕ , μ ν Γ μ ν θ ϕ , θ = ϕ + 3 Ω ϕ , X ; μ ϕ , μ = g μ ν ϕ , ν X ; μ = X ϕ , X = 1 2 ϕ 2 , ϕ 2 = 2 X , X = ϕ ϕ , ϕ = X ϕ ,
then Equation (14) is written as
α X α 1 ϕ + 3 Ω ϕ + α ( α 1 ) X α 2 X ϕ + ( 2 α 1 ) X α d d ϕ L n f = 0 ,
which can be transformed into
d d τ L n X + 6 Ω 2 α 1 + L n f 1 α = 0 ,
and in turn integrated, resulting in
f 1 2 α ( ϕ ) d ϕ = 2 λ e 3 Ω ( τ ) 2 α 1 d τ ,
where λ is an integration constant and has the same sign as f ( ϕ ) . In the gauge N = 24 e 3 Ω 2 α 1 , the right-hand side is
f 1 2 α ( ϕ ) d ϕ = 24 2 λ ( t t i ) ,
where t i is the initial time for the α scenario in the universe. At this point, we can introduce some structure for the function f ( ϕ ) and solve the integral.
When we consider the particular mathematical structure for the function f ( ϕ ) = p ϕ m or f ( ϕ ) = p e m ϕ with p and m constants, the classical solutions for the field ϕ in quadratures are
ϕ ( τ ) = ϕ ( τ i ) + ( 2 α + m ) 2 α p 1 2 α 2 λ e 3 Ω ( τ ) 2 α 1 d τ 2 α 2 α + m , f ( ϕ ) = p ϕ m , m 2 α E x p p 1 2 α 2 λ e 3 Ω ( τ ) 2 α 1 d τ , f ( ϕ ) = p ϕ 2 α , m = 2 α 2 α m L n m 2 α p 1 2 α 2 λ e 3 Ω ( τ ) 2 α 1 d τ , f ( ϕ ) = p e m ϕ , m 0 , p 1 2 α 2 λ e 3 Ω ( τ ) 2 α 1 d τ f ( ϕ ) = p , m = 0 .
The complete solution to the scalar field ϕ depends strongly on the mathematical structure of the scale factor Ω ( τ ) in the α scenario in our universe. In the gauge N = 24 e 3 Ω 2 α 1 , these solutions are
ϕ ( t ) = ϕ ( t i ) + 12 ( 2 α + m ) α p 1 2 α 2 λ t t i 2 α 2 α + m , f ( ϕ ) = p ϕ m , m 2 α E x p 24 p 1 2 α 2 λ t t i , f ( ϕ ) = p ϕ 2 α , m = 2 α 2 α m L n 12 m α p 1 2 α 2 λ t t i , f ( ϕ ) = p e m ϕ , m 0 , 24 p 1 2 α 2 λ t t i f ( ϕ ) = p , m = 0 ,
where t i and ϕ ( t i ) are the initial time and the scalar field in this time for the α scenario in the universe. In what follows, we perform the calculations to obtain the scale factor in some cases.

3. Lagrange and Hamilton Formalism

Introducing the line element (21) of the anisotropic Bianchi type I cosmological model into the Lagrangian (12), we have
L I = e 3 Ω 6 Ω ˙ 2 N 6 β ˙ + 2 N 6 β ˙ 2 N f ( ϕ ) 1 2 α ϕ ˙ 2 α N 2 α + 1 .
Using the standard definition of the momenta Π q μ = L q ˙ μ , where q μ are the coordinate fields q μ = ( Ω , β ± , ϕ ) , we obtain the momenta associated with each field
Π Ω = 12 N e 3 Ω Ω ˙ , Ω ˙ = e 3 Ω N Π Ω 12 , Π ± = 12 N e 3 Ω β ˙ ± , β ± ˙ = e 3 Ω N Π ± 12 , Π ϕ = f ( ϕ ) 1 2 α 2 α N 2 α 1 e 3 Ω ϕ ˙ 2 α 1 , ϕ ˙ = N e 3 Ω 2 α 1 2 α 2 α Π ϕ f ( ϕ ) 1 2 α 1 ,
and introducing them into the Lagrangian density, we obtain the canonical Lagrangian as L c a n o n i c a l = Π q μ q ˙ μ N H = Π q μ q ˙ μ H . When we perform the variation of this canonical Lagrangian with respect to N, δ L c a n o n i c a l δ N = 0 , we obtain the constraint H = 0 . In our model, this is the only constraint corresponding to the Hamiltonian density, which is weakly zero. So, the Hamiltonian is
H = N 24 e 3 2 α 1 Ω e 6 ( α 1 ) 2 α 1 Ω Π Ω 2 Π + 2 Π 2 12 ( 2 α 1 ) α 2 α 1 α f ( ϕ ) 1 2 α 1 Π ϕ 2 α 2 α 1 .

3.1. Exact Solution in the Gauge N = 24 e 3 2 α 1 Ω

Using the Hamilton equations for the momenta Π ˙ μ = H q μ and coordinates q ˙ μ = H Π μ , we have
Ω ˙ = 2 e 6 ( α 1 ) 2 α 1 Ω Π Ω
β + ˙ = 2 e 6 ( α 1 ) 2 α 1 Ω Π + ,
β ˙ = 2 e 6 ( α 1 ) 2 α 1 Ω Π ,
ϕ ˙ = 24 2 α 1 α f ( ϕ ) 1 2 α 1 Π ϕ 1 2 α 1 ,
Π ˙ Ω = 6 ( α 1 ) 2 α 1 e 6 ( α 1 ) 2 α 1 Ω Π Ω 2 Π + 2 Π 2 ,
Π ˙ ± = 0 , Π ± = p ± = c o n s t a n t ,
Π ˙ ϕ = 12 α 2 α 1 α 1 2 α 1 Π ϕ 2 α 2 α 1 f 2 α 2 α 1 f ˙ ϕ ˙ ,
solving Equation (39) using (36), we have Π ϕ = p ϕ f 1 2 α , with p ϕ an integration constant. With this result, and taking into account Equation (36), we obtain
f 1 2 α d ϕ = 24 p ϕ 2 α 1 α 1 2 α 1 t t i ,
that is, similar to (27), previously obtained, which was solved a Klein–Gordon-like equation directly. Using the Hamiltonian constraint and the solution to Equation (39) found previously, we have
e 6 ( α 1 ) 2 α 1 Ω Π Ω 2 Π + 2 Π 2 = 12 ( 2 α 1 ) α p ϕ 2 α 2 α 1 α 1 2 α 1 ,
then, the solution for the momenta becomes
Π Ω = η α t + p 0 ,
where the constant η α = 72 ( α 1 ) α p ϕ 2 α 2 α 1 α 1 2 α 1 and p 0 are constants of integration that, when introduced into the equation for Ω , give us the equation for the Ω function,
d Ω d t = 2 e 6 ( α 1 ) 2 α 1 Ω η α t + p 0 , e 6 ( α 1 ) 2 α 1 Ω = 6 ( α 1 ) 2 α 1 η α t 2 + 2 p 0 t + p 1 ,
whose solution becomes
Ω = 2 α 1 6 ( α 1 ) ln 6 ( α 1 ) 2 α 1 η α t 2 + 2 p 0 t + p 1 ,
and the solution for the scalar field is given by Equation (29). The solutions for the anisotropic function β ± are given by
β ± ( t ) = b ± ( 2 α 1 ) p ± 24 ( α 1 ) η α p 0 2 η α p 1 ln η α t + p 0 p 0 2 η α p 1 η α t + p 0 + p 0 2 η α p 1 ,
β ± = b ± 2 α 1 24 ( α 1 ) p ± η α λ α ln Σ ( t ) Σ + ( t ) ,
where
Σ ± ( t ) = η α t + p 0 ± λ α , and λ α = p 0 2 η α p 1 > 0 .
According to the last expressions, the radii associated with the Bianchi type I have the following behavior:
A ( t ) = e Ω + β + + 3 β = A 0 6 ( α 1 ) 2 α 1 Σ + ( t ) Σ ( t ) p + + 3 p 4 η α λ α ( η α t 2 + 2 p 0 t + p 1 ) 2 α 1 6 ( α 1 ) , B ( t ) = e Ω + β + 3 β = B 0 6 ( α 1 ) 2 α 1 Σ + ( t ) Σ ( t ) p + 3 p 4 η α λ α ( η α t 2 + 2 p 0 t + p 1 ) 2 α 1 6 ( α 1 ) , C ( t ) = e Ω 2 β + = C 0 [ 6 ( α 1 ) 2 α 1 ( Σ ( t ) Σ + ( t ) p + 2 η α λ α ( η α t 2 + 2 p 0 t + p 1 ) 2 α 1 6 ( α 1 ) ,
with the volume of this universe V ( t ) = A B C = e 3 Ω
V ( t ) = V 0 6 ( α 1 ) 2 α 1 ( η α t 2 + 2 p 0 t + p 1 ) 2 α 1 2 ( α 1 ) .
where we have graphed on different time scales in each scenario; in both cases the volume is increasing, as shown in Figure 1.

3.2. Exact Solution without Gauge N in the Time τ

For this case, the Hamilton procedure is not adequate, so we shall use the Hamilton–Jacobi procedure in order to find the solutions for the remaining mini-superspace variables, which arise by making the identification S ( Ω , β ± , ϕ ) q μ = Π μ in the Hamiltonian constraint (32), H = 0 , taking S ( Ω , β ± , ϕ ) = S Ω ( Ω ) + S + ( β + ) + S ( β ) + S ϕ ( ϕ ) , which results in
e 6 ( α 1 ) 2 α 1 Ω d S Ω d Ω 2 d S + d β + 2 d S d β 2 A α 1 f ( ϕ ) 1 2 α 1 d S ϕ d ϕ 2 α 2 α 1 = 0 .
Separating this equation, we have
e 6 ( α 1 ) 2 α 1 Ω d S Ω d Ω 2 d S + d β + 2 d S d β 2 = A α 1 f ( ϕ ) 1 2 α 1 d S ϕ d ϕ 2 α 2 α 1 = ϕ 2 .
with ϕ a separation constant. The solution in the variable ϕ is
Π ϕ = d S ϕ d ϕ = ϕ 2 A α 2 α 1 2 α f 1 2 α ( ϕ ) = p ϕ f 1 2 α ( ϕ ) ,
where S ( ϕ ) = p ϕ f 1 2 α ( ϕ ) d ϕ , obtaining similar results in the Hamilton procedure.
The specific values of the constants are p ϕ = ϕ 2 A α 2 α 1 2 α and A α = 12 ( 2 α 1 ) α 2 α 1 α 1 2 α 1 , in terms of the α parameter.
The other equations are read as
d S Ω d Ω 2 = + 2 + 2 + ϕ 2 e 6 ( α 1 ) 2 α 1 ,
d S + d β + 2 = + 2 , S + = s + ± + β + ,
d S d β 2 = 2 , S = s ± β ,
where i 2 are separation constants and s ± integration constants. On the other side, recalling the expressions for the momenta we can obtain solutions for Equations (52)–(54) in quadrature for the variable Ω and for α 1 .

3.3. Case for α 1

In this particular case, we have
d τ = 12 e 3 Ω d Ω 2 + ϕ 2 e 6 ( α 1 ) 2 α 1 Ω , 2 = + 2 + 2 ,
and for the anisotropic variables,
Δ β ± = ± 12 e 3 Ω ( τ ) d τ .
For solving Equation (55), we employ the transformation in the time variable d τ = e 3 2 α 1 Ω d T and U = 2 + ϕ 2 e 6 ( α 1 ) 2 α 1 Ω , so, d U = ϕ 2 6 ( α 1 ) 2 α 1 e 6 ( α 1 ) 2 α 1 Ω d Ω , resulting in
d T = 12 e 6 ( α 1 ) 2 α 1 Ω d Ω U = 2 ( 2 α 1 ) ( α 1 ) ϕ 2 d U U ,
and the solution is
T T 0 = 4 ( 2 α 1 ) ( α 1 ) ϕ 2 U U 0 ,
then, for the Ω variable, we have
Ω ( T ) = L n ϕ ( α 1 ) 4 ( 2 α 1 ) ( T T 0 ) + U 0 ϕ 2 2 ϕ 2 2 α 1 6 ( α 1 ) ,
and the time transformation becomes
d τ = ϕ ( α 1 ) 4 ( 2 α 1 ) ( T T 0 ) + U 0 ϕ 2 2 ϕ 2 1 2 ( α 1 ) d T .
To obtain the solutions in the time T for the anisotropic functions β ± ( T ) , we solve the integral
e 3 Ω ( τ ) d τ = e 6 ( α 1 ) 2 α 1 Ω ( T ) d T = d T ϕ ( α 1 ) 4 ( 2 α 1 ) ( T T 0 ) + U 0 ϕ 2 2 ϕ 2 = 4 ( 2 α 1 ) ( α 1 ) L n ϕ 2 ( α 1 ) 4 ( 2 α 1 ) ) ( T T 0 ) + μ 0 ϕ 2 ( α 1 ) 4 ( 2 α 1 ) ) ( T T 0 ) + μ 0 + ,
then, the anisotropic functions (56) become
β ± ( T ) = β ± ( T i ) ± 3 ( 2 α 1 ) ( α 1 ) L n ϕ 2 ( α 1 ) 4 ( 2 α 1 ) ) ( T T 0 ) + μ 0 ϕ 2 ( α 1 ) 4 ( 2 α 1 ) ) ( T T 0 ) + μ 0 + .
and the scalar field (28) takes the form
ϕ ( T ) = ϕ ( T i ) + ( 2 α + m ) 2 α p 1 2 α 2 λ ( T T i ) 2 α 2 α + m , f ( ϕ ) = p ϕ m , m 2 α , E x p p 1 2 α 2 λ ( T T i ) , f ( ϕ ) = p ϕ 2 α , m = 2 α , 2 α m L n m 2 α p 1 2 α 2 λ ( T T i ) , f ( ϕ ) = p e m ϕ , m 0 , p 1 2 α 2 λ ( T T i ) , f ( ϕ ) = p , m = 0 ,
On the other side, the only state when the time τ = T corresponds to the scenario α , which is calculated below.
  • Dust Scenario, α
For this particular case, we have
Ω ( τ ) = L n ϕ 8 ( τ τ 0 ) + U 0 ϕ 2 2 ϕ 2 1 3 ,
then, the volume function becomes
V ( τ ) = V 0 ϕ 8 ( τ τ 0 ) + U 0 ϕ 2 2 ϕ 2 ,
and the anisotropic functions are
β ± ( τ ) = β ± ( τ 0 ) ± 3 2 L n ϕ 2 8 ( τ τ 0 ) + μ 0 ϕ 2 8 ( τ τ 0 ) + μ 0 + .
We can see that the scalar field constant ϕ is huge, the anisotropic function tends to constant, and the anisotropic model can be an isotropic one. We rewrite the corresponding solutions in the scalar field (28) for this scenario:
ϕ ( τ ) = ϕ ( τ 0 ) + 2 λ ( τ τ 0 ) , f ( ϕ ) = p ϕ m , m 2 α , E x p 2 λ ( τ τ 0 ) , f ( ϕ ) = , m = 2 α , L n 2 λ ( τ τ i ) , f ( ϕ ) = p e 2 α ϕ , 2 λ ( τ τ 0 ) , f ( ϕ ) = p , m = 0 .

4. Quantum Regime

The WDW equation for these models is obtained by making the usual substitution Π q μ = i q μ in (32) and promoting the classical Hamiltonian density in the differential operator applied to the wave function Ψ ( Ω , β ± , ϕ ) , H ^ Ψ = 0 ; we have
2 e 6 ( α 1 ) 2 α 1 Ω 2 Ψ Ω 2 + 2 Ψ β + 2 + 2 Ψ β 2 12 ( 2 α 1 ) α 2 α 1 α f ( ϕ ) 1 2 α 1 2 α 2 α 1 2 α 2 α 1 ϕ 2 α 2 α 1 Ψ = 0 .
This fractional differential equation of degree β = 2 α 2 α 1 , belongs to different intervals depending on the value of the barotropic parameter [36]. We can write this equation in terms of the β parameter; we have
2 e 3 ( 2 β ) Ω 2 Ψ Ω 2 + 2 Ψ β + 2 + 2 Ψ β 2 24 β 2 α 1 α f ( ϕ ) 1 2 α 1 β β ϕ β Ψ = 0 .
For simplicity, the factor e 3 ( 2 β ) Ω may be the factor ordered with Π ^ Ω and f 1 2 α 1 ( ϕ ) may be the factor ordered with β ϕ β . In many ways, what we employ what might be called a semi-general factor ordering, which, in this case, would order the terms e 3 ( 2 β ) Ω Π ^ Ω 2 as e ( 3 ( 2 β ) Q ) Ω Ω e Q Ω Ω = e 3 ( 2 β ) Ω Ω 2 + Q e 3 ( 2 β ) Ω Ω , where Q is any real constant that measures the ambiguity in the factor ordering in the variables Ω and its corresponding momenta. For the other factor ordering, we make the following calculation which, in this case, would order the terms g ( ϕ ) f 1 2 α 1 ( ϕ ) β ϕ β , where in the particular case we choose g ( ϕ ) = ϕ s , similarly to f ( ϕ ) in the classical case, is
f 1 2 α 1 ( ϕ ) ϕ s β / 2 ϕ β / 2 ϕ s β / 2 ϕ β / 2 = f 1 2 α 1 ( ϕ ) β ϕ β + f 1 2 α 1 ( ϕ ) ϕ s β / 2 ϕ β / 2 ϕ s β / 2 ϕ β / 2 ,
where the Caputo fractional derivative of β / 2 ϕ β / 2 ϕ s becomes [1],
0 c D x β / 2 ϕ s = Γ ( s + 1 ) Γ ( s β / 2 + 1 ) ϕ s β / 2 .
Thus, Equation (68) is rewritten as
f 1 2 α 1 ( ϕ ) ϕ s β / 2 ϕ β / 2 ϕ s β / 2 ϕ β / 2 = f 1 2 α 1 ( ϕ ) β ϕ β + f 1 2 α 1 ( ϕ ) ϕ s Γ ( s + 1 ) Γ ( s + 1 β / 2 ) ϕ s β / 2 β / 2 ϕ β / 2 = f 1 2 α 1 ( ϕ ) β ϕ β + f 1 2 α 1 ( ϕ ) Γ ( s + 1 ) Γ ( s + 1 β / 2 ) ϕ β / 2 β / 2 ϕ β / 2
Assuming this factor ordering for the Wheeler–DeWitt equation, we obtain
e 3 ( 2 β ) Ω 2 Ψ Ω 2 + Q Ψ Ω + 2 Ψ β + 2 + 2 Ψ β 2 24 β 2 α 1 α f ( ϕ ) 1 2 α 1 β β ϕ β Ψ 24 β 2 α 1 α f ( ϕ ) 1 2 α 1 β ϕ s Γ ( s + 1 ) Γ ( s + 1 μ ) β 2 ϕ β 2 Ψ = 0 .
Using the ansatz for the wave function Ψ ( Ω , β + , β , ϕ ) = A ( Ω ) B + ( β + ) B ( β ) C ( ϕ ) , we obtain the following differential equations on the corresponding variables:
d 2 A d Ω 2 Q d A d Ω ± ρ 2 e 3 ( 2 β ) Ω + ρ 1 2 A = 0 ,
d 2 B + d β + 2 ρ 2 2 B + = 0 ,
d 2 B d β 2 ρ 3 2 B = 0 , ρ 3 2 = ρ 2 2 ρ 1 2
ϕ γ d 2 γ C d ϕ 2 γ + Γ ( s + 1 ) Γ ( s + 1 γ ) d γ C d ϕ γ ± p ϕ α 2 α 1 1 2 α 1 γ ρ 2 12 2 γ C = 0 , β = 2 γ , 0 < γ 1 .
We can see that the fractional differential Equation (75) has variable coefficients, so to solve it we can use the fractional power series [37,38], also f ( ϕ ) = p ϕ ϕ m , with p ϕ a constant and choosing m 2 α 1 = γ , as a particular case, which implies that the parameter m = α in the sense that the γ parameter is in accordance with its original definition (see Equation (75) and the definition of the β parameter, or the equation in the text after Equation (76)).
Following the book of Polyanin [53] (page 179.10), we find the solution for the first equation, considering different values in the factor-ordering parameter (we take the corresponding sign as minus in the constant ρ 2 ):
A ( Ω ) = e Q Ω 2 C 1 K ν ρ 3 ( 1 γ ) e 3 ( 1 γ ) Ω + C 2 I ν ρ 3 ( 1 γ ) e 3 ( 1 γ ) Ω ,
where K ν ( z ) and I ν ( z ) are the modified Bessel functions, and order ν = Q 2 + 4 ρ 1 2 6 ( 1 γ ) with γ = β 2 = α 2 α 1 the new order in the fractional derivative. However, for physical conditions we will only take the modified Bessel K ν .
The corresponding quantum solutions for Equations (73) and (74) are
B + = a 0 e ρ 2 β + + a 1 e ρ 2 β + ,
B = b 0 e ρ 3 β + b 1 e ρ 3 β .
with ρ 2 and ρ i 2 being separation constants.
The solution of Equation (75) with a positive sign and f ( ϕ ) = p ϕ = c o n s t a n t , with zero factor ordering, may be obtained by applying direct and inverse Laplace transforms [20,36], providing
C + ( ϕ , γ ) = E 2 γ z 2 , z = p ϕ α 2 α 1 1 2 ( 2 α ) γ ρ 2 3 γ ϕ γ , 0 < γ 1 ,
where E 2 γ is the Mittag–Leffler function (7), then, the probability density of the wave function for this particular case becomes
| Ψ | 2 = ψ 0 2 e Q Ω ± 2 ρ 2 β + ± 2 ρ 3 β K ν 2 ρ 3 ( 1 γ ) e 3 ( 1 γ ) Ω E 2 γ 2 ( z 2 ) ,
and its corresponding plot for two values in the ordering parameter Q is shown in the Figure 2.
Table 2 shows the differential equations obtained from (75), depending on the values of ω X , α , and γ .

Solution to FDE Associated with the Different State Evolutions

We write the fractional differential Equation (75) as follows:
ϕ γ d 2 γ C d ϕ 2 γ + A d γ C d ϕ γ + B C = 0 , 0 < γ 1 ,
where we have made the simplifications A = Γ ( s + 1 ) Γ ( s + 1 γ ) and B ( α , γ ) = ± p ϕ α 2 α 1 1 2 α 1 γ μ 2 12 2 γ . The last linear fractional differential Equation (81) will be solved using the fractional power series [37,38]
C = n = 0 a n ϕ n γ .
Then, the fractional derivatives are
d γ C d γ ϕ = n = 1 a n Γ [ n γ + 1 ] Γ [ ( n 1 ) γ + 1 ] ϕ ( n 1 ) γ , d 2 γ C d 2 γ ϕ = n = 2 a n Γ [ n γ + 1 ] Γ [ ( n 1 ) γ + 1 ] Γ [ ( n 1 ) γ + 1 ] Γ [ ( n 2 ) γ + 1 ] ϕ ( n 2 ) γ .
Substituting expressions (83) into (81), we obtain
n = 2 a n Γ [ n γ + 1 ] Γ [ ( n 2 ) γ + 1 ] ϕ ( n 1 ) γ + A n = 1 a n Γ [ n γ + 1 ] Γ [ ( n 1 ) γ + 1 ] ϕ ( n 1 ) γ + B ( α , γ ) n = 0 a n ϕ n γ = 0 .
Now, taking = n 1 into the first and second terms, and n = into the third term of (84), we have
= 1 a + 1 Γ [ ( + 1 ) γ + 1 ] Γ [ ( 1 ) γ + 1 ] ϕ γ + A = 0 a + 1 Γ [ ( + 1 ) γ + 1 ] Γ [ γ + 1 ] ϕ γ + B ( α , γ ) = 0 a ϕ γ = 0 .
Shifting one place in the second and third summations, we have
= 1 a + 1 Γ [ ( + 1 ) γ + 1 ] Γ [ ( 1 ) γ + 1 ] ϕ γ + A Γ [ γ + 1 ] Γ [ 1 ] a 1 + A = 1 a + 1 Γ [ ( + 1 ) γ + 1 ] Γ [ γ + 1 ] ϕ γ + B ( α , γ ) a 0 + B ( α , γ ) = 1 a ϕ γ = 0 .
From the last expression (86), we obtain ( s 0 )
a 1 = B ( α , γ ) A Γ ( γ + 1 ) = Γ [ s + 1 γ ] B ( α , γ ) Γ [ s + 1 ] Γ [ γ + 1 ] a 0 ,
and the recurrence relationship between the parameters a is
a + 1 = Γ [ s + 1 γ ] B ( α , γ ) Γ [ γ + 1 ] Γ [ ( 1 ) γ + 1 ] Γ [ ( + 1 ) γ + 1 ] Γ [ s + 1 γ ] Γ [ γ + 1 ] + Γ [ s + 1 ] Γ [ ( 1 ) γ + 1 ] a , 1 .
Some terms of this relation are
a 1 = Γ [ s + 1 γ ] B ( α , γ ) Γ [ s + 1 ] Γ [ γ + 1 ] a 0 , a 2 = Γ [ s + 1 γ ] B ( α , γ ) 2 Γ [ s + 1 ] Γ [ 2 γ + 1 ] Γ [ s + 1 ] + Γ [ s + 1 γ ] Γ [ γ + 1 ] a 0 , a 3 = Γ [ γ + 1 ] Γ [ s + 1 ] Γ [ 3 γ + 1 ] Γ [ s + 1 γ ] Γ [ 2 γ + 1 ] + Γ [ s + 1 ] Γ [ γ + 1 ] × Γ [ s + 1 γ ] B ( α , γ ) 3 Γ [ s + 1 ] + Γ [ s + 1 γ ] Γ [ γ + 1 ] a 0 , a 4 = Γ [ γ + 1 ] Γ [ 2 γ + 1 ] Γ [ s + 1 ] Γ [ 4 γ + 1 ] Γ [ s + 1 γ ] Γ [ 3 γ + 1 ] + Γ [ s + 1 ] Γ [ 2 γ + 1 ] × ( Γ [ s + 1 γ ] B ( α , γ ) ) 4 Γ [ s + 1 γ ] Γ [ 2 γ + 1 ] + Γ [ s + 1 ] Γ [ γ + 1 ] Γ [ s + 1 ] + Γ [ s + 1 γ ] Γ [ γ + 1 ] a 0 ,
Then, the solution of the fractional Equation (81) has the form
C s , α , γ = a 0 1 Γ [ s + 1 γ ] B ( α , γ ) Γ [ s + 1 ] Γ [ γ + 1 ] ϕ γ + Γ [ s + 1 γ ] B ( α , γ ) 2 Γ [ s + 1 ] Γ [ 2 γ + 1 ] Γ [ s + 1 ] + Γ [ s + 1 γ ] Γ [ γ + 1 ] ϕ 2 γ + .
For the dust-like scenario (see Table 2), α and γ = 1 2 , then B ( , 1 2 ) = μ 2 24 . The solution associated with this fractional differential equation is given in fractional series form by
C s , , 1 2 = a 0 1 2 Γ [ s + 1 2 ] B ( , 1 2 ) Γ [ s + 1 ] π ϕ 1 2 + 2 Γ [ s + 1 2 ] B ( , 1 2 ) 2 Γ [ s + 1 ] ( 2 Γ [ s + 1 ] + Γ [ s + 1 2 ] π ) ϕ 2 2 Γ [ s + 1 2 ] B ( , 1 2 ) 3 3 π Γ [ s + 1 ] ( 2 Γ [ s + 1 2 ] + Γ [ s + 1 ] π ) ( 2 Γ [ s + 1 ] + Γ [ s + 1 2 ] π ) ϕ 3 2 + ,
where we employed Γ [ 3 2 ] = π 2 and Γ [ 5 2 ] = 3 π 4 .
For the radiation stage, α = 2 and γ = 2 3 , then B ( 2 , 2 3 ) = p ϕ 3 μ 2 18 4 3 , the solution for the fractional differential equation is
C s , 2 , 2 3 = a 0 1 3 Γ [ s + 1 3 ] B ( 2 , 2 3 ) 2 Γ [ s + 1 ] Γ [ 2 3 ] ϕ 2 3 + 3 Γ [ s + 1 3 ] B ( 2 , 2 3 ) 2 4 Γ [ s + 1 ] Γ [ 4 3 ] 3 Γ [ s + 1 ] + 2 Γ [ s + 1 3 ] Γ [ 2 3 ] ϕ 4 3 2 3 Γ [ 2 3 ] 3 Γ [ s + 1 3 ] B ( 2 , 2 3 ) 3 2 Γ [ s + 1 ] 2 Γ [ s + 1 3 ] Γ [ 4 3 ] + Γ [ s + 1 ] Γ [ 2 3 ] 3 Γ [ s + 1 ] + 2 Γ [ s + 1 3 ] Γ [ 2 3 ] ϕ 8 3 + ,
with Γ [ 2 3 ] = 1.35412 and Γ [ 4 3 ] = 0.89298 .
We are going to present the graphical behavior of the wave function for the dust-like case, in which the solution is constrained to the variables Ω and ϕ . This means we are going to shrink the directions β + and β . So, the wave function takes the form
Ψ ( Ω , β + , β , ϕ ) = e Q Ω 2 C 1 K ν ρ 3 ( 1 γ ) e 3 ( 1 γ ) Ω × 1 a 1 ϕ 1 2 + a 2 ϕ a 3 ϕ 3 2 + a 4 ϕ 2 a 5 ϕ 5 2 + a 6 ϕ 3 a 7 ϕ 7 2 + a 8 ϕ 4 + .
By restraining ourselves to the values of ν = Q 2 + 4 ρ 1 2 6 ( 1 γ ) ) , Q = 10 , 6 , 2 , 1.6 , ρ = 1 2 , γ = 1 2 , s = 1 , ρ 1 = 0.1 , and C 1 = 1 , we see that the wave function ((93)) can be rewritten as
Ψ 2 ( Ω , ϕ ) = e Q Ω K ν 2 2 3 e 3 Ω 2 1 1 24 ϕ 1 / 2 + π 576 ( 4 + π ) ϕ 2 π 41472 ( 4 + π ) ϕ 3 / 2 + π 3 2 663552 ( 8 + 3 π ) ( 4 + π ) ϕ 2 π 3 2 69672960 ( 8 + 3 π ) ( 4 + π ) ϕ 5 2 + π 5 2 334430208 ( 8 + 3 π ) ( 4 + π ) ( 32 + 15 π ) ϕ 3 π 5 2 45649723392 ( 8 + 3 π ) ( 4 + π ) ( 32 + 15 π ) ϕ 7 2 + 5 π 7 2 208684449792 ( 8 + 3 π ) ( 4 + π ) ( 32 + 15 π ) ( 192 + 105 π ) ϕ 4 2 ,
where we have taken the cut-off order in ϕ 4 and the ’a’ parameters are read from (89). In the following, we present some plots of the probability density of the wave function, including the factor-ordering parameter Q and particular values in the parameter ρ , ρ 1 and the particular value to the ordering parameter s = 1. In all of them, we observe that for any value of Q, the probability density decays with respect to the scale factor, but has a different evolution in the scalar field. For small Q’s, the quantum universe has considerable existence in the evolution with respect to the scale factor and then decays. On the other hand, for large Q’s, this interval is small. What we can say about the evolution of the scalar field is that at small Q’s, the scalar field appears faster than for large Q’s, which enters late, but has existed forever (Figure 3).

5. Conclusions

Unlike the previous work [36], in the present paper we employed a barotropic equation with perfect fluid for the energy–momentum tensor in the k-essence scalar field into the Lagrangian and Hamiltonian formalism, obtaining the momentum of a scalar field with fractional numbers, while the momentum of the scale factor appears in the usual way. We obtained the classical solutions for different scenarios in the universe, employing different times ( t , T ( τ ) , τ ) . In the quantum scheme, we include the factor-ordering problem, and we find a fractional differential equation for the scalar field with variable coefficients, which was solved using the fractional series expansion. With this in mind, we visualize two alternatives in our analysis; the first one is within the traditional expectation over the behavior of the probability density, that the best candidates for quantum solutions are those that have a damping behavior with respect to the scale factor, appearing in all scenarios under our study, without saying anything about the scalar field. The other scenario is when we keep the scale factor, and we consider the values of the scalar field as significant in the quantum regime, appearing in various scenarios in the behavior of the universe; mainly in those where the universe has a huge behavior, for example, in the actual epoch, where the scalar field appears as background.
In other words, the interpretation of the probability density of the unnormalized wave function is given when we demand that Ψ does not diverge when the scale factor A (or Ω ) goes to infinity, and the scalar field is arbitrary. However, the evolution with the scalar field is important in this class of theory and others as it appears in some stages of the evolution of our universe.
In reference [54], the gravitational action integral is altered by hand, leading to a modified Friedmann-type equation. They employed the dynamical system approach in order to find the balance points providing a range for the order of the fractional derivatives in their investigation of the cosmological universe, and they mention that it can be confirmed that the solutions isotropize at a late time. In our approach, this occurs when the ϕ is huge (see Equation (64)); due to that, the anisotropic parameters become constants. On the other hand, as we use the volume of the universe V = A B C = e 3 Ω ( τ ) on any timescale, it depends only on the Ω function and not on the anisotropic parameters β ± . Whereas the scale factors A, B, and C depend on these anisotropic parameters. It remains to be studied whether the fractional derivatives alter the gravitational part and how the universe’s singularity can be avoided because, from our approach, this part needs to be revised since the gravitational part is not altered from the point of view of the equations of motion. It is modified in the scalar field part.
It would be interesting to extend the Bohm-type semi-classical formalism in this context, which we will explore in future work. Much work has been made in this direction [55,56,57,58,59,60,61], where the quantum potential emerges as the imaginary part in the Bohm formalism, appearing as a constraint equation. In this sense, in Reference [62] an approach appears that is based on the semi-classical limit of fractionary quantum cosmology using the Riesz derivative. It would be interesting to continue under our focus, where the corresponding Friedmann equation and the Hubble parameter depend on Levy’s fractional parameter, which is associated with the concept of the Lévy path in the corresponding quantum cosmology.
We briefly illustrate the main results of this work.
  • Using the k-essence formalism in a general way, applied to the anisotropic Bianchi type I cosmological model, we found the Hamiltonian density in the scalar field momenta raised to powers of non-integers, which produces in the quantum scheme a fractional differential equation in a natural way. We include the factor-ordering problem in both variables ( Ω , ϕ ) and its momenta ( Π Ω , Π ϕ ) , with the order β = 2 α 2 α 1 , where α ( 1 , ) , and it was solved in a general way, we include two particular scenarios of our universe.
  • We found the solution in the classical scheme employing two gauges, N = 24 e 3 Ω , for two forms of the function f ( ϕ ) in the time t; however, when we let the Lagrange multiplier N, we need to employ a transformed time T ( τ ) for solving the classical equation and, only in the dust era, we recover the gauge time τ .
  • In the quantum regime, when we include the factor-ordering problem, the fractional differential equation in the scalar field appears with variable coefficients, and it was necessary to use the fractional series expansion to solve it in a general way.
  • In one of our analyses presented on the probability density, we consider the values of the scalar field as significant in the quantum regime, appearing in various scenarios in the behavior of the universe, mainly in those where the universe has a huge behavior; for example, in the actual epoch, where the scalar field appears as a background, the quantum regime appears with big values, but it presents a moderate development in other scenarios with different ordering parameters Q and s.

Author Contributions

Conceptualization, J.S., J.J.R. and L.T.-S.; Methodology, J.S., J.J.R. and L.T.-S.; Formal analysis, J.S., J.J.R. and L.T.-S.; Investigation, J.S., J.J.R. and L.T.-S.; Writing—original draft, J.S. All authors have read and agreed to the published version of the manuscript.

Funding

J.S. was partially supported by PROMEP grants UGTO-CA-3. The authors were partially supported by SNI-CONACyT. J. Rosales is supported by PROMEP grants UGTO-CA-20 nonlinear photonics and the Department of Electrical Engineering. This work is part of the collaboration within the Instituto Avanzado de Cosmología and Red PROMEP: Gravitation and Mathematical Physics under project Quantum aspects of gravity in cosmological models, phenomenology and geometry of space-time. L.T.S. is supported by Secretaria de Investigación y Posgrado del Instituto Politécnico Nacional, grant SIP20211444. Many of the calculations were performed by the Symbolic Program REDUCE 3.8.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Podlubny, I. Fractional Differential Equations; Academic Press: New York, NY, USA, 1999. [Google Scholar]
  2. Uchaikin, V. Fractional Derivatives for Physicists and Engineers; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  3. Herrmann, R. Fractional Calculus: An Introduction For Physicists; World Scientific Publishing Co., Pte., Ltd.: Singapore, 2011. [Google Scholar]
  4. Caputo, M.; Mainarde, F.A. A new dissipation model based on memory mechanism. Pure Appl. Geophys. 1971, 91, 134–137. [Google Scholar] [CrossRef]
  5. Wyss, W. Fractional diffusion equation. J. Math. Phys. 1986, 27, 2782–2785. [Google Scholar] [CrossRef]
  6. Westerlund, S. Capacitor theory. IEEE Trans. Dielectr. Electr. Insul. 1994, 1, 826–839. [Google Scholar] [CrossRef]
  7. Rosu, H.C.; Madueno, A.L.; Socorro, J. Transform of Riccati equation of constant coefficients through fractional procedure. J. Phys. A Math. Gen. 2003, 36, 1087–1093. [Google Scholar] [CrossRef]
  8. Magin, R.L. Fractional Calculus in Bioengineering; Begell House Publisher: Rodding, Denmark, 2006. [Google Scholar]
  9. Tarasov, V.E. Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields and Media; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  10. Duarte Ortigueira, M.; Martynyuk, V.; Kosenkov, V.; Batista, G.A. A new look at the capacitor theory. Fractal Fract. 2023, 7, 86. [Google Scholar] [CrossRef]
  11. Laskin, N. Fractional Quantum Mechnics; World Scientific Publishing Co., Pte., Ltd.: Singapore, 2018. [Google Scholar]
  12. Moniz, P.V.; Jalalzadeh, S. From Fractional Quantum Mechanics to Quantum Cosmology: An Overture. Mathematics 2020, 8, 313. [Google Scholar] [CrossRef]
  13. Jalalzadeh, S.; da Silva, F.R.; Moniz, P.V. Prospecting black hole thermodynamics with fractional quantum mechanics. Eur. Phys. J. C 2021, 81, 632. [Google Scholar] [CrossRef]
  14. Jalalzadeh, S.; Moniz, P.V. Challenging Routes in Quantum Cosmology; World Scientific Publiching Company: Singapore, 2023. [Google Scholar]
  15. Valério, D.; Ortigueira, M.D.; Lopes, A.M. How Many Fractional Derivatives are There? Mathematics 2022, 10, 737. [Google Scholar] [CrossRef]
  16. Ortigueira, M.D. A New Look at the Initial Condition Problem. Mathematics 2022, 10, 1771. [Google Scholar] [CrossRef]
  17. Ortigueira, M.D.; Bohannan, G.W. Fractional Scale Calculus: Hadamard vs.Liouville. Fractal Fract. 2023, 7, 296. [Google Scholar] [CrossRef]
  18. Ortigueira, M.D.; Bengochea, G. A Simple Solution for the General Fractional Ambartsumian Equation. Appl. Sci. 2023, 13, 871. [Google Scholar] [CrossRef]
  19. Godínez, F.A.; Rosales, J.J.; Esparza-Posadas, M.F. Newton’s Law of Cooling with Caputo Derivative: Consistent Dimensionality to Compare with Experiments. Fractals 2022, 30, 2250187. [Google Scholar] [CrossRef]
  20. Rosales-García, J.; Andrade-Lucio, J.A.; Shulika, O. Conformable derivative applied to experimental Newton law of cooling. Rev. Mex. FÍsica 2020, 66, 224–227. [Google Scholar] [CrossRef]
  21. El-Nabulsi, A.R. Fractional Unstable Euclidean Universe. EJTP 2005, 8, 1–11. [Google Scholar]
  22. El-Nabulsi, R.A. Gravitons in fractional action cosmology. Int. J. Theor. Phys. 2012, 51, 3978–3992. [Google Scholar] [CrossRef]
  23. El-Nabulsi, R.A. Non-minimal coupling in fractional action cosmology. Indian J. Phys. 2013, 87, 835–840. [Google Scholar] [CrossRef]
  24. Rami, E.N.A. Fractional action oscillating phantom cosmology with conformal coupling. Eur. Phys. J. Plus 2015, 130, 102. [Google Scholar] [CrossRef]
  25. El-Nabulsi, R.A. A Cosmology Governed by a Fractional Differential Equation and the Generalized Kilbas-Saigo-Mittag-Leffler Function. Int. J. Theor. Phys. 2016, 55, 625–635. [Google Scholar] [CrossRef]
  26. El-Nabulsi, R.A. Implications of the Ornstein-Uhlenbeck-like fractional differential equation in cosmology. Rev. Mex. Fis. 2016, 62, 240. [Google Scholar]
  27. El-Nabulsi, R.A. Fractional Action Cosmology with Variable Order Parameter. Int. J. Theor. Phys. 2017, 56, 1159–1182. [Google Scholar] [CrossRef]
  28. El-Nabulsi, R.A. Wormholes in fractional action cosmology. Can. J. Phys. 2017, 95, 605–609. [Google Scholar] [CrossRef]
  29. Rasouli, S.M.M.; Jalalzadeh, S.; Moniz, P.V. Broadening quantum cosmology with a fractional whirl. Mod. Phys. Lett. A 2021, 36, 2140005. [Google Scholar] [CrossRef]
  30. Jalalzadeh, S.; Costa, E.W.O.; Moniz, P.V. de Sitter fractional quantum cosmology. Phys. Rev. D 2022, 105, L121901. [Google Scholar] [CrossRef]
  31. Rasouli, S.M.M.; Costa, E.W.O.; Moniz, P.V.; Jalalzadeh, S. Inflation and fractional quantum cosmology. Fractal Fract. 2022, 6, 655. [Google Scholar] [CrossRef]
  32. El-Nabulsi, A.R. Cosmology with Fractional Action Principle. Rom. Rep. Phys. 2007, 39, 763–771. [Google Scholar]
  33. García-Aspeitia, M.A.; Anaya, G.F.; Almada, A.H.; León, G.; Magaña, J. Cosmology under the fractional calculus approach. Mon. Not. 2022, 517, 4813–4826. [Google Scholar] [CrossRef]
  34. Leon, G.; García-Aspeitia, M.A.; Fernández-Anaya, G.; Hernández-Almada, A.; Magaña, J.; González, E. Cosmology under the fractional calculus approach: A possible H0 tension resolution. Mon. Not. R. Astron. Soc. 2022, 517, 4813–4826. [Google Scholar]
  35. Socorro, J.; Pimentel, L.O.; Espinoza-García, A. Classical Bianchi type I cosmology in K-essence theory. Adv. High Energy Phys. 2014, 2014, 805164. [Google Scholar] [CrossRef]
  36. Socorro, J.; Rosales, J.J. Quantum fraccionary cosmology: K-essence theory. Universe 2023, 9, 185. [Google Scholar] [CrossRef]
  37. El-Ajou, A.; Arqub, O.A.; Zhour, Z.A.; Momani, S. New Results on Fractional Power Series: Theories and Applications. Entropy 2013, 15, 5305–5323. [Google Scholar] [CrossRef]
  38. Rida, S.Z.; Arafa, A.A.M. New Method for Solving Linear Fractional Differential Equations. Int. J. Differ. Equations 2011, 2011, 814132. [Google Scholar] [CrossRef]
  39. Tarasov, V.E. No violation of the Leibniz rule. No fractional derivative. Commun. Nonlinear Sci. Numer. Simul. 2013, 18, 2945–2948. [Google Scholar] [CrossRef]
  40. Ortigueira, M.D.; Tenreiro Machado, J.A. What is a fractional derivative? J. Comput. Phys. 2015, 293, 4–13. [Google Scholar] [CrossRef]
  41. Teodoro, G.S.; Machado, J.T.; De Oliveira, E.C. A review of definitions of fractional derivatives and other operators. J. Comput. Phys. 2019, 388, 195–208. [Google Scholar] [CrossRef]
  42. Sandev, T.; Iomin, A. Special Functions of Fractional Calculus: Applications to Diffusion and Random Search Processes; World Scientific Publishing Co., Pte., Ltd.: Singapore, 2023. [Google Scholar]
  43. Erdelyi, A.; Magnus, W.; Oberthettinger, F.; Tricomi, F.G. Higher Transcendental Functions; McGraw-Hill: New York, NY, USA, 1955. [Google Scholar]
  44. Haubold, H.J.; Mathai, A.M.; Saxena, R.K. Mittag-Leffler functions and their applications. J. Appl. Math. 2011, 2011, 298628. [Google Scholar] [CrossRef]
  45. García, A.E.; Socorro, J.; Pimentel, L.O. Quantum Bianchi type IX cosmology in K-essence theory. Int. J. Theor. Phys. 2014, 53, 3066–3077. [Google Scholar] [CrossRef]
  46. de Putter, R.; Linder, E.V. Kinetic k-essence and Quintessence. Astropart. Phys. 2007, 28, 263. [Google Scholar] [CrossRef]
  47. Chiba, T.; Dutta, S.; Scherrer, R.J. Slow-roll k-essence. Phys. Rev. D 2009, 80, 043517. [Google Scholar] [CrossRef]
  48. Bose, N.; Majumdar, A.S. A k-essence model of inflation, dark matter and dark energy. Phys. Rev. D 2009, 79, 103517. [Google Scholar] [CrossRef]
  49. Arroja, F.; Sasaki, M. Note on the equivalence of a barotropic perfect fluid with a k-essence scalar field. Phys. Rev. D 2010, 81, 107301. [Google Scholar] [CrossRef]
  50. García, L.A.; Tejeiro, J.M.; Castaneda, L. K-essence scalar field as dynamical dark energy. arXiv 2012, arXiv:1210.5259. [Google Scholar]
  51. Bilic, N.; Tupper, G.; Viollier, R. Unification of dark matter and dark energy: The inhomogeneous Chaplygin gas. Phys. Lett. B 2002, 535, 17. [Google Scholar] [CrossRef]
  52. Bento, M.; Bertolami, O.; Sen, A. Dynamics of dark energy. Phys. Rev. D 2002, 66, 043507. [Google Scholar] [CrossRef]
  53. Polyanin, A.C.; Zaitsev, V.F. Handbook of Exact Solutions for Ordinary Differential Equations, 2nd ed.; Chapman Hall/CRC: Boca Raton, FL, USA, 2003. [Google Scholar]
  54. Micolta-Riascos, B.; Millano, A.D.; Genly, L.; Erices, C.; Paliathanasis, A. Revisiting Fractional Cosmology. Fractal Fract. 2023, 7, 149. [Google Scholar] [CrossRef]
  55. Socorro, J.; Rodríguez, P.A.; Nunez, O.E.; Hernández, R.; Espinoza-García, A. Quintom Potential from Quantum Anisotropic Cosmological Models. In Open Questions in Cosmology; Gonzalo, J.O., Ed.; Chapter 9; InTech Open: London, UK, 2012; pp. 219–240. ISBN 978-953-51-0880-1. [Google Scholar]
  56. Socorro, J.; Romero, P.; Pimentel, L.O.; Aguero, A. Quintom potentials from quantum cosmology using the FRW cosmological model. Int. J. Theor. Phys. 2013, 52, 2722–2734. [Google Scholar] [CrossRef]
  57. Zampeli, A.; Pailas, T.; Terzis, P.A.; Christodoulakis, T. Conditional symmetries in axisymmetric quantum cosmologies with scalar fields and the fate of the classical singularities. JCAP 2016, 2016, 66. [Google Scholar] [CrossRef]
  58. Socorro, J.; Nuñez, O.E. Scalar potentials with multi-scalar fields from quantum cosmology and supersymetric quantum mechanics. Eur. Phys. J. Plus 2017, 132, 168. [Google Scholar] [CrossRef]
  59. Paliathanasis, A. Dust fluid component from Lie symmetries in Scalar field Cosmology. Mod. Phys. Lett. A 2017, 32, 1750206. [Google Scholar] [CrossRef]
  60. Zampeli, A.; Paliathanasis, A. Quantization of inhomogeneous spacetimes with cosmological constant term. Class. Quantum Grav. 2021, 38, 165012. [Google Scholar] [CrossRef]
  61. Paliathanasis, A. Quantum potentiality in Inhomogeneous Cosmology. Universe 2021, 7, 52. [Google Scholar] [CrossRef]
  62. Oliveira Costa, E.W.; da Silva, P.F., Jr.; Rasouli, S.M.M.; Jalalzadeh, S. The estimation of the universe’s age as per fractional cosmology. arXiv 2023, arXiv:2310.09464. [Google Scholar]
Figure 1. Volume of the universe into the radiation and dust ages, respectively, according to Table 1, we choose p ϕ = 2 and p 1 = 1 .
Figure 1. Volume of the universe into the radiation and dust ages, respectively, according to Table 1, we choose p ϕ = 2 and p 1 = 1 .
Fractalfract 07 00814 g001
Figure 2. Probability density (80) of the universe dominated by dust era in this stage of the universe shows that the probability density has a decay, both in Ω and ϕ . Both are plotted with ρ = 0.5 , γ = 1 2 , and ρ 1 = 0.1 . (a) Q = 6 , ν = 2.00111 ; (b) Q = 2 , ν = 0.669992 .
Figure 2. Probability density (80) of the universe dominated by dust era in this stage of the universe shows that the probability density has a decay, both in Ω and ϕ . Both are plotted with ρ = 0.5 , γ = 1 2 , and ρ 1 = 0.1 . (a) Q = 6 , ν = 2.00111 ; (b) Q = 2 , ν = 0.669992 .
Fractalfract 07 00814 g002
Figure 3. Probability density (94) of the universe dominated by dust era in this stage of the universe shows that the probability density has a decay in Ω and exhibits considerable growth for certain values of ϕ . The plots have the parameters ρ = 0.5 , γ = 1 2 , and ρ 1 = 0.1 . (a) Q = 10 , ν = 3.334 ; (b) Q = 6 , ν = 2.00111 ; (c) Q = 2 , ν = 0.0669992 ; (d) Q = 1.6 , ν = 0.537489 .
Figure 3. Probability density (94) of the universe dominated by dust era in this stage of the universe shows that the probability density has a decay in Ω and exhibits considerable growth for certain values of ϕ . The plots have the parameters ρ = 0.5 , γ = 1 2 , and ρ 1 = 0.1 . (a) Q = 10 , ν = 3.334 ; (b) Q = 6 , ν = 2.00111 ; (c) Q = 2 , ν = 0.0669992 ; (d) Q = 1.6 , ν = 0.537489 .
Fractalfract 07 00814 g003
Table 1. States of the universe’s evolution according to the barotropic parameter ω X .
Table 1. States of the universe’s evolution according to the barotropic parameter ω X .
ω X α G ( X ) State of Evolution
11XStiff matter
1 3 2 X 2 Radiation
0 X m ,     m Dust-like
1 01,     f ( ϕ ) = Λ = c t e Inflation
1 3 −1 1 X Inflation-like
2 3 1 4 1 X 4 Inflation-like
Table 2. Fractionary equation in the field ϕ according to the barotropic parameter ω X .
Table 2. Fractionary equation in the field ϕ according to the barotropic parameter ω X .
ω X α γ Fractionary Equation
111 ϕ d 2 C d ϕ 2 + Γ ( s + 1 ) Γ ( s ) d C d ϕ ± p ϕ μ 2 24 2 C = 0
1 3 2 2 3 ϕ 2 3 d 4 3 C d ϕ 4 3 + Γ ( s + 1 ) Γ ( s + 1 3 ) d 2 3 C d ϕ 2 3 ± ( p ϕ ) 1 3 μ 2 18 4 3 C = 0
0 1 2 ϕ 1 2 d C d ϕ + Γ ( s + 1 ) Γ ( s + 1 2 ) d 1 2 C d ϕ 1 2 ± μ 2 24 C = 0
1 00 w i t h o u t e q u a t i o n
1 3 1 1 3 ϕ 1 3 d 2 3 C d ϕ 2 3 + Γ ( s + 1 ) Γ ( s + 2 3 ) d 1 3 C d ϕ 1 3 ± p ϕ α 2 α 1 1 2 α 1 μ 2 36 2 3 C = 0
2 3 1 4 1 6 ϕ 1 6 d 1 3 C d ϕ 1 3 + Γ ( s + 1 ) Γ ( s + 5 6 ) d 1 6 C d ϕ 1 6 ± p ϕ α 2 α 1 1 2 α 1 μ 2 72 1 3 C = 0
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Socorro, J.; Rosales, J.J.; Toledo-Sesma, L. Anisotropic Fractional Cosmology: K-Essence Theory. Fractal Fract. 2023, 7, 814. https://doi.org/10.3390/fractalfract7110814

AMA Style

Socorro J, Rosales JJ, Toledo-Sesma L. Anisotropic Fractional Cosmology: K-Essence Theory. Fractal and Fractional. 2023; 7(11):814. https://doi.org/10.3390/fractalfract7110814

Chicago/Turabian Style

Socorro, José, J. Juan Rosales, and Leonel Toledo-Sesma. 2023. "Anisotropic Fractional Cosmology: K-Essence Theory" Fractal and Fractional 7, no. 11: 814. https://doi.org/10.3390/fractalfract7110814

APA Style

Socorro, J., Rosales, J. J., & Toledo-Sesma, L. (2023). Anisotropic Fractional Cosmology: K-Essence Theory. Fractal and Fractional, 7(11), 814. https://doi.org/10.3390/fractalfract7110814

Article Metrics

Back to TopTop