Next Article in Journal
On Rotationally Symmetrical Planar Networks and Their Local Fractional Metric Dimension
Next Article in Special Issue
Traversable Wormhole Solutions Admitting Noether Symmetry in f(R,T2) Theory
Previous Article in Journal
Video Anomaly Detection Based on Attention Mechanism
Previous Article in Special Issue
Observational Constraints on F(T,TG) Gravity with Hubble’s Parametrization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability Properties of Self-Similar Solutions in Symmetric Teleparallel f(Q)-Cosmology

by
Andronikos Paliathanasis
1,2
1
Institute of Systems Science, Durban University of Technology, P.O. Box 1334, Durban 4000, South Africa
2
Departamento de Matemáticas, Universidad Católica del Norte, Avda. Angamos 0610, Casilla, Antofagasta 1280, Chile
Symmetry 2023, 15(2), 529; https://doi.org/10.3390/sym15020529
Submission received: 1 February 2023 / Revised: 9 February 2023 / Accepted: 14 February 2023 / Published: 16 February 2023
(This article belongs to the Special Issue Symmetry and Asymmetry in Gravity Research)

Abstract

:
Self-similar cosmological solutions correspond to spacetimes that admit a homothetic symmetry. The physical properties of self-similar solutions can describe important eras of the cosmological evolution. Recently, self-similar cosmological solutions were derived for symmetric teleparallel f Q -theory with different types of connections. In this work, we study the stability properties of the self-similar cosmological solutions in order to investigate the effects of the different connections on the stability properties of the cosmic history. For the background geometry, we consider the isotropic Friedmann–Lemaître–Robertson–Walker space and the anisotropic and homogeneous Bianchi I space, for which we investigate the stability properties of Kasner-like universes.

1. Introduction

On the cosmological scale, General Relativity (GR) fails to provide an explanation for the observed phenomena [1,2,3,4,5]. At the present time, the universe is under in acceleration phase, which is attributed to a matter source called dark energy. Dark energy has not been directly observed until now, but new dynamical terms are attributed to it, which are necessarily introduced into the field equations of GR in order to provide anti-gravitational effects and explain the late-time acceleration. Because dark energy has been observed indirectly only, the nature and the origin of the dark energy is unknown to cosmologists. There is a specific approach adopted by cosmologists to explain the acceleration through the introduction of geometric invariants in the gravitational Action Integral [6], from which as a result the dark energy is a geometric effect of the gravitational theory. These new theories of gravity are known as modified theories of gravity [7,8,9,10]. The plethora of different models in the literature, which are based on the use of different geometric invariants for the modification of the Action Integral, form a “zoology” of gravitational theories, f R -theory [11,12], Gauss–Bonnet gravity [13,14], teleparallel theory of gravity [15,16,17,18] and others [19,20,21].
The fundamental invariant of GR is the curvature scalar R related to the Levi–Civita connection of the metric tensor which describes the physical space. However, an arbitrary connection can describe a manifold with a curvature scalar R, torsion scalar T and non-metricity scalar Q. These three scalars, R, T and Q, are also known as the geometrical trinity of gravity [22]. The Levi–Civita connection that is considered in GR is torsionless and with zero metricity parts. Thus, only the curvature scalar is used for the definition of the gravitational action. On the other hand, the existence of an unholonomic frame leads to a connection with a non-zero torsion scalar. The Weitzenböck connection [23] is a curvature-less connection, which leads to the Teleparallel Equivalent of General Relativity (TEGR), for which the torsion plays the role of the gravitational force [24]. Moreover, it has been found that a geometry with a torsion-free connection with non-zero scalar Q that describes a flat geometry, that is, the Riemann tensor is zero, is equivalent to GR and it is known as Symmetric Teleparallel General Relativity (STGR) [25].
When the Gravitational Action Integral is constructed by the scalars R, T and Q, the resulting three theories, GR, TEGR and STGR, are indeed equivalent. However, this equivalency is lost when additional invariants are introduced into the gravitational Action Integral in order to solve the dark energy problem. f R -gravity [9] is a fourth-order theory and includes various important models, which have been used to describe various epochs of the cosmological evolution such as the inflationary epoch [26,27] or the late-time acceleration phase related to the dark energy [28,29,30]. In a similar way, in teleparallelism a plethora of f T gravitational models have been introduced for the description of the observable phenomena. For a review of teleparallelism, we refer the reader to [31]. Furthermore, in [32], it is shown that in teleparallel gravity inflation can occur without the existence of an inflaton field, while f T theory can mimic the dark energy model [33]. See also [34,35,36] and references therein.
Recently, the f Q -gravity has drawn the attention of cosmologists. f Q -gravity can provide answers to various observation phenomena in cosmology, so as to provide dynamical terms in the field equations to explain the effects of dark energy, and to pass various tests by using the current observational data. In [37], the authors applied cosmological constraints for f Q models, which can reproduce the Λ CDM expansion history; it was found that the f Q theory is in agreement with the observations and it is supported by data over the Λ CDM model. A similar result was found dependently in [38]. The nonlinear f Q = α Q + β Q n theory studied in [39] as dark energy model and the cosmographic parameters derived. In [40], the authors considered f Q models, for which the equation of state parameter for the cosmological model can cross the phantom divine line. In [41,42,43] external and internal solutions of compact stars were studied. New cosmological asymptotic solutions in f Q theory with initial singularity were derived in [44]. See also the results presented in [45,46,47]. Bouncing cosmological scenarios in f Q -gravity were investigated in [48]. Moreover, in [49], the method of statefinder diagnostics was applied in symmetric teleparallel theory and it was found that for two nonlinear f Q models the field equations provide an evolution, which describes acceleration in the late time due to the dark energy component evolving from the non-metricity components. Finally, the ADM formulation and Hamiltonian analysis of f Q theory are presented in [50].
Unlike GR, in f Q -gravity it is possible to separate gravity from the inertial effects [51]. This is possible because we make use of a flat connection pertaining to the existence of affine coordinates in the coincident gauge in which the covariant derivative is reduced to the usual partial derivative. Indeed, because in STGR the curvature tensor is zero, there always exists a coordinate system in which the connection becomes zero, Γ μ ν λ = 0 [52].
The coincidence gauge can always be derived under a coordinate transformation. It has been found that the use of different connections for the definition of the non-metricity scalar, Q, affects the dynamics of f Q -gravity [53,54]. For the Friedmann–Lemaître–Robertson–Walker (FLRW) cosmological model four different families of connections have been found [55]. In [56], the existence of self-similar solutions was investigated in f Q -theory for FLRW cosmology. Self-similar solutions are of huge interest, because they can describe specific eras of the cosmological evolution. The self-similar solutions can be related to the asymptotic behaviour of the general dynamics [44]. In addition, the existence of the limit of solutions that are described by GR in f Q -theory is of special interest for the viability of the theory.
In this study, we are interested in the stability properties of the self-similar solutions in order to understand the effects of different connections during the different epochs of the cosmological history. Moreover, we are interested to investigate the existence of asymptotic scaling solutions that describe inflation to see if f Q symmetric teleparallel theory for different kinds of connections can reproduce the inflationary epoch. Such analysis gives constraints on the viability of the different kind of connections. The structure of the paper is as follows:
In Section 2, we present the basic properties and definitions for the f Q symmetric teleparallel theory. Section 3 includes the main results of this work in which we investigate the stability properties of self-similar cosmological solutions for the different connections in the coincidence gauge which generate the cosmological field equations of FLRW geometry in symmetric teleparallel theory. In Section 4, we discuss the existence of self-similar solutions in an anisotropic and homogeneous Bianchi I geometry in the context of f Q -theory. We find that a generalized Kasner solution is supported by f Q -theory in the coincidence gauge. Finally, in Section 5 we draw our conclusions.

2. Basic Properties and Definitions for the f ( Q ) Symmetric Teleparallel Theory

Assume a four-dimensional manifold V 4 with metric tensor g μ ν and generic connection Γ μ ν κ . From the connection, we can define the Riemann tensor
R λ μ ν κ = Γ λ ν κ x μ Γ λ μ κ x ν + Γ λ ν σ Γ μ σ κ Γ λ μ σ Γ μ σ κ ,
the torsion tensor
T μ ν λ = Γ μ ν λ Γ ν μ λ
and the non-metricity tensor
Q λ μ ν = g μ ν x λ Γ λ μ σ g σ ν Γ λ ν σ g μ σ .
In STGR, it follows that the torsion tensor vanishes, that is, T μ ν λ = 0 , and that the geometry is flat, i.e., R λ μ ν κ = 0 , while only the non-metricity tensor survives. The non-metricity scalar is defined as
Q = Q λ μ ν P λ μ ν ,
where P μ ν λ is the non-metricity conjugate tensor expressed as
P μ ν λ = 1 4 Q μ ν λ + 1 2 Q ( μ λ ν ) ( μ λ ν ) + 1 4 Q λ Q ¯ λ g μ ν 1 4 δ ( μ λ Q ν )
are Q λ = Q λ μ μ , Q ¯ λ = Q λ μ μ .
The gravitational Action Integral in In f Q -theory is as follows [57]
S = 1 2 d 4 x g f ( Q ) + d 4 x g L M + λ κ λ μ ν R λ μ ν κ + τ λ μ ν T μ ν λ ,
in which g = det ( g μ ν ) , L M is the Lagrangian density for the matter source and λ κ λ μ ν , τ λ μ ν are two Lagrange multipliers that impose the flatness of the connection, i.e., R λ μ ν κ = 0 , and that the connection is symmetric, that is, T μ ν λ = 0 .
The gravitational field equations are derived as
2 g λ g f ( Q ) P μ ν λ 1 2 f ( Q ) g μ ν + f ( Q ) P μ ρ σ Q ν ρ σ 2 Q ρ σ μ P ρ σ ν ρ σ = T μ ν ,
where f ( Q ) = d f Q d Q and T μ ν = 2 g g L M g μ ν is the energy-momentum tensor, which describes the matter components of the gravitational fluid.
Moreover, variation of the Action Integral (1) with respect to the connection gives the field equations
μ ν g f ( Q ) P μ ν σ μ ν = 0 .
The field equations (2) can be written in the equivalent form
f ( Q ) G μ ν + 1 2 g μ ν f ( Q ) Q f ( Q ) + 2 f ( Q ) λ Q P μ ν λ = T μ ν ,
in which G μ ν is the Einstein-tensor, where G μ ν = R ˜ μ ν 1 2 g μ ν R ˜ , with R ˜ μ ν and R ˜ are the Riemannian Ricci tensor and scalar, respectively, which are constructed by the Levi–Civita connection.
We can easily define the effective energy-momentum tensor
T μ ν = 1 f ( Q ) 1 2 g μ ν f ( Q ) Q f ( Q ) + 2 f ( Q ) λ Q P μ ν λ
which contributes the geometrodynamical degrees of freedom. Therefore, with the latter definition, the field equations can be written in the known form similar to that of general relativity
G μ ν = T μ ν + 1 f ( Q ) T μ ν .

3. FLRW Cosmology

According to the cosmological principle, at a large scale the universe is isotropic and homogeneous described by the FLRW line element
d s 2 = N ( t ) 2 d t 2 + a ( t ) 2 d r 2 1 k r 2 + r 2 d θ 2 + sin 2 θ d ϕ 2 ,
where a t is the scale factor and N t is the lapse function. The Hubble function is defined as H = 1 N a ˙ a , where a ˙ = d a d t . Parameter k is the spatial curvature and k = 0 corresponds to a spatially flat space, k = 1 is a closed universe and k = 1 is an open universe. Another important characteristic of the FLRW geometry is that the three-dimensional hypersurface is maximally symmetric and it admits a six-dimensional Killing algebra.
All the compatible connections for the FLRW for the symmetric teleparallel theory by enforcing on a generic connection the six Killing symmetries of the background geometry and the requirement to be flat have been derived before in [53,54]. In the following, without loss of generalitym, we assume N t = 1 .
For k = 0 , there are three compatible connections with common components
Γ θ θ r = r , Γ ϕ ϕ r = r sin 2 θ , Γ r θ θ = Γ θ r θ = Γ r ϕ ϕ = Γ ϕ r ϕ = 1 r , Γ ϕ ϕ θ = sin θ cos θ , Γ θ ϕ ϕ = Γ ϕ θ ϕ = cot θ .
The first connection, Γ 1 , has the additional non-zero component
Γ t t t = γ ( t ) ,
where γ ( t ) is a function of the time variable t.
The second connection, Γ 2 , has the extra non-zero components
Γ t t t = γ ˙ ( t ) γ ( t ) + γ ( t ) , Γ t r r = Γ r t r = Γ t θ θ = Γ θ t θ = Γ t ϕ ϕ = Γ ϕ t ϕ = γ ( t ) ,
while the third connection, Γ 3 , has the non-zero components for the connection
Γ t t t = γ ˙ ( t ) γ ( t ) , Γ r r t = γ ( t ) , Γ θ θ t = γ ( t ) r 2 , Γ ϕ ϕ t = γ ( t ) r 2 sin 2 θ .
For an FLRW geometry with non-zero spatial curvature, that is, k 0 , there exists only one compatible connection, namely Γ 4 k with non-zero coefficients
Γ t t t = k + γ ˙ ( t ) γ ( t ) , Γ r r t = γ ( t ) 1 k r 2 Γ θ θ t = γ ( t ) r 2 , Γ ϕ ϕ t = γ ( t ) r 2 sin 2 ( θ ) Γ t r r = Γ r t r = Γ t θ θ = Γ θ t θ = Γ t ϕ ϕ = Γ ϕ t ϕ = k γ ( t ) , Γ r r r = k r 1 k r 2 , Γ θ θ r = r 1 k r 2 , Γ ϕ ϕ r = r sin 2 ( θ ) 1 k r 2 , Γ r θ θ = Γ θ r θ = Γ r ϕ ϕ = Γ ϕ r ϕ = 1 r , Γ ϕ ϕ θ = sin θ cos θ , Γ θ ϕ ϕ = Γ ϕ θ ϕ = cot θ .
We remark that the latter connection in the limit k = 0 reduces to the third connection for the flat case. That is an important observation because these two connections may relate the field equations for the flat and non-flat spatial curvature cases.
Indeed, for each different connection the resulting cosmological field equations are different. The existence of self-similar solutions was the subject of study in [56]. For each connection, the functional form of the arbitrary function f Q is explicitly determined by the assumption that the background space and the connection admit a homothetic symmetry vector related to the existence of the self-similar solution, similarly to the case of GR.
Below, we present the self-similar solutions determined before and we investigate the stability of the solutions. Such an analysis provides important information about the evolution of the dynamical variables. Simultaneously, constraints can be constructed for the free parameters of the model and for the initial conditions of the theory.

3.1. First Connection

For the first connection, Γ 1 , we derive the non-metricity scalar
Q = 6 a ˙ 2 a 2 = 6 H 2 .
Moreover, the gravitational field equations are (4)
3 H 2 f ( Q ) + 1 2 f ( Q ) Q f ( Q ) = ρ , 2 N d d t f ( Q ) H 3 H 2 f ( Q ) 1 2 f ( Q ) Q f ( Q ) = p .
The parameters ρ and p correspond to the energy density and pressure components of the energy-momentum tensor for an external fluid. As it was found before in [56], there does not exist any self-similar solution for the first-connection without an external matter source.
Assume now that the energy-momentum tensor describes an ideal gas with a constant equation of state parameter w, that is p = w ρ , and the background geometry describes a self-similar universe with a t = a 0 t λ . Then, from the conservation law of energy for the energy-momentum tensor, it follows that ρ t = ρ 0 t 3 λ 1 + w , where ρ 0 is a constant.
By substituting in the field equation we find
f ( Q ) = f 0 Q + f 1 ( Q ) 3 2 λ ( w + 1 ) ,
where f 1 = 2 λ 3 λ ( 1 + w ) ρ 0 6 3 2 λ ( w + 1 ) 1 3 λ 1 + w .
The first coefficient of (7) does not contribute to the field equations. Indeed, the term g Q is a total derivative and when it is replaced by (5) in the Action Integral (1) it can be eliminated. Thus, the important coefficient is the second one, the power f Q = f 1 ( Q ) 3 2 λ ( w + 1 ) . The effective cosmological fluid for the self-similar solution has a constant effective equation of the state parameter w e f f = 1 + 2 3 λ , which does not depend upon the external fluid. That means that the gravitational effects of external fluid are overlapped or neutralized by the f Q terms.
With the use of (7) it is easy to see that the scaling solution a t = a 0 t λ is the analytic solution to the problem, which means that the solution is stable, because it is the unique solution. We remark that the stability analysis presented in [44] is for f Q = Q + f 1 Q n and the results there should not be confused with our conclusion for this work.

3.2. Second Connection

For the second connection we derive the non-metricity scalar
Q = 6 H 2 + 9 γ H + 3 γ ˙
in which the gauge function γ t is involved and it contributes to the cosmological dynamics.
The field equations in the case of a vacuum read
3 H 2 f ( Q ) + 1 2 f ( Q ) Q f ( Q ) + 3 γ Q ˙ f ( Q ) 2 = 0 ,
2 d d t f ( Q ) H 3 γ Q ˙ f ( Q ) = 0 ,
where there exists the constraint equation for the function f Q given by
Q ˙ 2 f ( Q ) + Q ¨ + 3 H Q ˙ f ( Q ) = 0 .
The self-similar solution [56] a t = a 0 t 1 3 m which describes a universe dominated by stiff fluid, i.e., w e f f = 1 , exists for f Q = f 1 Q + f 2 Q ln Q , and the gauge function γ t = C t 2 9 t ln t , where C = f 1 + f 2 3 + ln 4 3 9 f 2 . We remark that f Q = f 1 Q + f 2 Q ln Q for f 2 = f ¯ 2 f 1 and very small values of f ¯ 2 it can be seen as the first term of the Taylor expansion of the f Q = f 1 Q 1 + f ¯ 2 theory. This means that for small deviation from GR, the resulting solution describes a stiff fluid source and not the Minkowski spacetime. Recall that we do not consider the existence of an external matter source. This result differs from that found in the similar case of f R -cosmology [58].
A more general self-similar scaling solution, a t = a 0 t λ , λ 0 , 1 3 , is supported by the cosmological field equations for f ( Q ) = f 1 Q + f 2 Q 2 3 ( 1 λ ) and gauge function γ = q 0 t 3 λ + 12 λ 2 t 6 ( 3 λ 1 ) t 2 , where f 2 = 6 f 1 ( λ 1 ) λ q 0 2 3 ( λ 1 ) .

Stability Properties

We proceed with the investigation of the stability properties for the self-similar solution with λ 0 , 1 3 . We consider perturbations only on the scale factors and not on the gauge function γ t .
Thus, for the scale factor a t = a 0 t λ we derive the Hubble function H t = λ t . We assume now that H t = λ t + ε δ H t . We substitute into (9) and the first-order perturbations give the first-order differential equation
0 = 3 q 0 2 t 6 λ 8 39 λ + 36 λ 2 δ H t δ H ˙ + 16 t 2 λ 3 3 λ 2 2 3 λ δ H + t δ H ˙ 8 q 0 t 1 + 3 λ λ 6 λ 1 δ H + t δ H ˙ .
This has the analytic solution
δ H t = δ H 0 t 8 9 λ 3 q 0 t 3 λ + 4 t 2 3 λ λ 6 λ 7 q 0 t 3 λ + 4 λ 2 t 6 λ 3 .
Therefore, for λ > 1 3 and 8 + 36 λ 2 39 λ < 0 , the perturbations decay for large values of t, i.e., lim t + δ H t = 0 , which means that the self-similar solution a t = a 0 t λ is a stable solution and it is an attractor. For these values of the free parameter λ , it follows that 0.1754 < w e f f < 1 , that is, there is no stable self-similar solution which describes an accelerated universe.
However, we have not considered the equation of motion (11) in the stability analysis. We remark that we have not perturbed the gauge. In this case, the first order perturbations of (11) decay when 4 + 18 λ 2 21 λ . Thus, the field equations are stable when 1 3 < λ < 13 + 41 24 .

3.3. Third Connection

For the third connection the non-metricity scalar is calculated to be
Q = 6 H 2 + 3 γ a 2 H + 3 γ ˙ a 2 ,
while cosmological field equations are
3 H 2 f ( Q ) + 1 2 f ( Q ) Q f ( Q ) 3 γ Q ˙ f ( Q ) 2 a 2 = 0 ,
2 d d t f ( Q ) H + γ Q ˙ f ( Q ) a 2 = 0 ,
with the constraint
Q ˙ 2 f ( Q ) + Q ¨ + Q ˙ H + 2 γ ˙ γ f ( Q ) = 0 .
For the field equations of the third connection, we found that they support a self-similar universe with a t = a 0 t λ , for which the connection is also self-similar under a homothetic transformation, when γ t = γ 0 t 2 λ 1 , for the function [56]
f Q = f 1 Q 5 λ 1 2 ,
with the constraint equation γ 0 = 2 λ 2 5 λ 5 λ 3 and λ 1 5 , 3 5 or γ 0 = 2 λ 3 λ 1 and λ > 3 5 .
We proceed with the study of the perturbations for the two exact scaling solutions.

Stability Properties

Similarly to the previous connection, we performed a perturbation on the Hubble function, H t = λ t + ε δ H t and we substitute in (13) for the power-law function (16).
For γ 0 = 2 λ 2 5 λ 5 λ 3 , the first-order perturbations provide the ordinary differential equation
t 8 15 λ δ H ˙ + 32 + λ 30 λ 77 δ H = 0 ,
with analytic solution δ H t = δ H 0 t 32 77 λ + 30 λ 2 15 λ 8 , from which it follows that the perturbations decay and the self-similar solution is stable for 32 77 λ + 30 λ 2 15 λ 8 < 0 , i.e., λ < 0.52 and 0.53 < λ < 2.04 .
For the second case with γ 0 = 2 λ 3 λ 1 the first-order perturbations are always zero which means that the self-similar solution is always stable.
For the first case and the perturbations of the constraint equation, (15) we find that the field equations admit the self-similar solution a t = a 0 t λ as an attractor when 32 77 λ + 30 λ 2 15 λ 8 < 0 and 16 λ 7 + 45 λ 15 λ 8 < 0 , that is, 0.68 < λ < 0.52 and 0.53 < λ < 2.04 .

3.4. Fourth Connection

For a non-zero spatial curvature, k, the non-metricity scalar becomes
Q = 6 H 2 + 3 γ a 2 H + 3 γ ˙ a 2 + k 6 a 2 + 3 γ γ ˙ γ 3 H .
The gravitational field equations in a vacuum are
3 H 2 f ( Q ) + 1 2 f ( Q ) Q f ( Q ) 3 γ Q ˙ f ( Q ) 2 a 2 + 3 k f ( Q ) a 2 Q ˙ f ( Q ) 2 γ = 0 , 2 d d t f ( Q ) H + γ Q ˙ f ( Q ) a 2 + 4 k a 2 f ( Q ) = 0
and the equation of motion for the connection is written as
Q ˙ 2 f ( Q ) 1 + k a 2 γ 2 + Q ¨ 1 + k a 2 γ 2 + Q ˙ 1 + 3 k a 2 γ 2 H + 2 γ ˙ γ f ( Q ) = 0 .
For the last connection, which describes an FLRW universe with non-zero spatial curvature, for an arbitrary curvature the self-similar solution was found to be a t = a 0 t [56] for the gauge function γ ( t ) = γ 0 t , and
f Q = Q 2 ,
with the constraint equation γ 0 = 3 a 0 2 .
On the other hand, for the open universe, i.e., k = 1 , there exists the additional power theory
f ( Q ) = q 0 Q a 0 1 2 a 0 ,
where q 0 is a constant and it is required in order a 0 ± 1 and γ ( t ) = ± a 0 t . The case a 0 2 = 1 describes the Milne spacetime while the latter function, f Q , is reduced to that of STGR, i.e., f Q = Q .

Stability Properties

We now study the stability properties of the scaling-solution a t = a 0 t for the quadratic f Q = Q 2 model. In a similar approach to the previous one, we do not change the gauge γ t . Thus, we substitute H = 1 t + ε δ H t into (18) and from the first-order perturbations we define the ordinary differential equation
3 t δ H ˙ 5 δ H = 0 ,
that is, δ H t = δ H 0 t 5 3 . Hence, the scaling solution is always unstable.
For the open universe and the power-law theory f ( Q ) = q 0 Q a 0 1 2 a 0 it follows that the first-order perturbations of Equation (18) do not decay, that is, the Milne-line scaling solution is always unstable.

4. Existence of Self-Similar Solutions in an Anisotropic Bianchi I Cosmology in the Context of f Q -Gravity

Consider now the Bianchi I line element
d s 2 = N t 2 d t 2 + a 2 t d x 2 + b 2 t d y 2 + c 2 t d z 2 ,
in which a , b and c are the three scale factors. When a = b and a = c , the Bianchi I line-element reduces to that of the spatially flat FLRW cosmology. Without loss of generality we have selected a constant lapse function, i.e., N t = 1 .
For the Bianchi I geometry we follow [59] and in the coincidence gauge we assume the following metricity scalar
Q = 2 a ˙ b ˙ a b + a ˙ c ˙ a b + b ˙ c ˙ b c ,
where in the limit of isotropization the metricity scalar for the connection Γ 1 of FLRW spacetime is recovered. In [41] the isotropization of the Bianchi I spacetime in f Q theory was investigated; specifically, it was found that due to inflation which is related with the non-metricity components the Bianchi I spacetime can become isotropic. However, in this study, we are interested in a Kasner-like solution which always describes an anisotropic universe.
In the case in which there is no additional matter source, the gravitational field equations are
a ˙ b ˙ a b + a ˙ c ˙ a c + b ˙ c ˙ b c f ( Q ) + 1 2 f ( Q ) Q f ( Q ) = 0 ,
2 d d t f Q b ˙ c + c ˙ + b c f ( Q ) Q f ( Q ) = 0 ,
2 d d t f Q a ˙ c + c ˙ + a c f ( Q ) Q f ( Q ) = 0 ,
2 d d t f Q a ˙ b ˙ + a ˙ + a b f ( Q ) Q f ( Q ) = 0 .
Consider now the self-similar Kasner-like spacetime with scale-factors
a t = t p 1 , b t = t p 2 and c t = t p 3 ,
for the lapse function N t = 1 .
For these scale-factors the Bianchi I spacetime (20) admits the proper homothetic vector field [60]
X = t t + p 1 1 x x + p 2 1 y y + p 3 1 z z .
In the case of General Relativity, the exponents p 1 , p 2 and p 3 satisfy the Kasner relations. They are [61]
p 1 + p 2 + p 3 = 1
and
p 1 2 + p 2 2 + p 3 2 = 1
For a discussion and the importance of the Kasner spacetime we refer the reader to [62].
Hence, by replacing the scale factors (25) in the field Equations (21)–(24) and the power-law theory f Q = f 0 Q μ , it follows that the exponents p 1 , p 2 and p 3 satisfy the generalized Kasner relations
p 1 + p 2 + p 3 = 2 μ 1 ,
and
p 1 2 + p 2 2 + p 3 2 = 2 μ 1 2 .
However, in that case Q = 0 . Indeed, for any non-singular function f Q , where f Q 0 = 0 , this Kasner-like solution solves the field equations. For the power-law f Q = f 0 Q μ theory, the Kasner-like solution is the analytic solution. Hence there is no reason to study the stability of the self-similar solution. We remark that the generalized Kasner relations can be written in the form of the Kasner relations
P 1 + P 2 + P 3 = 1
P 1 2 + P 2 2 + P 3 3 = 1 ,
where P 1 = p 1 2 μ 1 , P 2 = p 2 2 μ 1 and P 3 = p 3 2 μ 1 .
It is important to mention that in terms of dynamics the field Equations (21)–(24) are of the same form as those of f T teleparallel theory of gravity which means that the algebraic results for the stability analysis presented in [63] are valid.
In a future work it will be of special interest to investigate the effects of different connections in Bianchi geometries.

5. Discussion

In this study, we investigated the stability properties of some self-similar cosmological solutions in symmetric teleparallel gravity. In particular, for the case of isotropic FLRW geometry we derived the field equations and self-similar solutions for the four different connections of the coincidence gauge. Three of the different corrections describe a spatially flat FLRW geometry, while the fourth connection describes an FLRW geometry with non-zero spatial curvature. The families of the connections introduce a gauge function, which affects the dynamics of the gravitational field equations for a non-linear f Q -theory.
We considered self-similar solutions that describe ideal gas solutions in the case of spatially flat FLRW geometry, and Milne-like geometries for non-zero spatially FLRW spacetime. The solutions exist for specific functions f Q and gauge functions for the corresponding connection. The self-similar solutions do not describe the full degrees of freedom of the field equations which means that they are special solutions and not the general solutions. In particular, they are asymptotic exact solutions of the generic solutions in each case.
Because of the complexity of the field equations, we were not able to study the phase-space by using dimensionless variables in the H-normalization approach. Instead, we took the standard approach to consider small perturbations of the scale factor, i.e., of the Hubble function, in the region of the similarity solutions. For each connection, we derived the differential equation for the perturbations and we solved each one of them. It was found that the behaviour of the perturbations depends upon the free parameters of the self-similar solution. Hence, the free parameter for the f Q -theory was a constraint for each connection with the requirement that the self-similar solution be an attractor for the field equations. The case for which the asymptotic self-similar solutions describe accelerated universes was studied.
Finally, we consider the anisotropic and homogeneous Bianchi I geometry and for the connection in the coincidence gauge we proved that Kasner-like solutions exist for a power-law theory f Q = Q μ , and for any non-singular function f Q with f Q 0 = 0 . Finally, the field equations for the background space are dynamically equivalent with the field equation with that of f T -theory in teleparallelism. Thus, the stability properties of the Kasner-like solutions are similar to those of the study [63].
In a future study, we plan to investigate further the self-similar solutions for anisotropic and homogeneous cosmological spacetimes.

Funding

This work was partially financially supported in part by the National Research Foundation of South Africa (Grant Numbers 131604). The author thanks the support of Vicerrectoría de Investigación y Desarrollo Tecnológico (Vridt) at Universidad Católica del Norte through Núcleo de Investigación Geometría Diferencial y Aplicaciones, Resolución Vridt No-096/2022.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Perlmutter, S.; Alderin, G.; Goldhabe, G.; Kno, R.A.; Nugen, P.; Castr, P.G.; Deustu, S.; Fabbr, S.; Gooba, A.; Groo, D.E.; et al. Measurements of Ω and Λ from 42 high-redshift supernovae. Astrophys. J. 1998, 517, 565. [Google Scholar] [CrossRef]
  2. Riess, A.G.; Filippenk, A.V.; Challi, P.; Clocchiatt, A.; Dierck, A.; Garnavic, P.M.; Gillilan, R.L.; Hoga, C.J.; Jh, S.; Kirshne, R.P.; et al. Observational Evidence from Supernovae for an Accelerating Universe and a Cosmological Constant. Astron. J. 1998, 116, 1009. [Google Scholar] [CrossRef] [Green Version]
  3. Suzuki, N.; Rubin, D.; Lidman, C.; Aldering, G.; Amanullah, R.; Barbary, K.; Barrientos, L.F.; Botyanszki, J.; Brodwin, M.; Connolly, N.; et al. The Hubble Space Telescope Cluster Supernova Survey. V. Improving the Dark-energy Constraints above z > 1 and Building an Early-type-hosted Supernova Sample. Astrophys. J. 2012, 746, 85. [Google Scholar] [CrossRef] [Green Version]
  4. Jarosik, N.; Bennett, C.L.; Dunkley, J.; Gold, B.; Greason, M.R.; Halpern, M.; Hill, R.S.; Hinshaw, G.; Kogut, A.; Komatsu, E.; et al. Seven-year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpretation. Astrophys. J. Suppl. 2011, 192, 18. [Google Scholar] [CrossRef] [Green Version]
  5. Ade, P.A.; Aghanim, N.; Alves, M.I.R.; Armitage-Caplan, C.; Arnaud, M.; Ashdown, M.; Atrio-Barandela, F.; Aumont, J.; Aussel, H.; Baccigalupi, C.; et al. Planck 2013 results. XVI. Cosmological Parameters. Astron. Astrophys. 2014, 571, A16. [Google Scholar]
  6. Tsujikawa, S. Modified gravity models of dark energy. Lect. Notes Phys. 2010, 800, 99. [Google Scholar]
  7. Nojiri, S.; Odintsov, S.D.; Oikonomou, V.K. Modified gravity theories on a nutshell: Inflation, bounce and late-time evolution. Phys. Rept. 2017, 692, 1. [Google Scholar] [CrossRef] [Green Version]
  8. Bamba, K.; Odintsov, S.D. Inflationary Cosmology in Modified Gravity Theories. Symmetry 2015, 7, 220. [Google Scholar] [CrossRef] [Green Version]
  9. Nojiri, S.; Odintsov, S.D. Unified cosmic history in modified gravity: From F(R) theory to Lorentz non-invariant models. Phys. Rep. 2011, 505, 59–144. [Google Scholar] [CrossRef] [Green Version]
  10. Clifton, T.; Ferreira, P.G.; Padilla, A.; Skordis, C. Modified gravity and cosmology. Phys. Rep. 2012, 513, 1–189. [Google Scholar] [CrossRef] [Green Version]
  11. Buchdahl, H.A. Non-Linear Lagrangians and Cosmological Theory. Mon. Not. Roy. Astron. Soc. 1970, 150, 1. [Google Scholar] [CrossRef] [Green Version]
  12. Barrow, J.D.; Cotsakis, S. Inflation and the conformal structure of higher-order gravity theories. Phys. Lett. B 1988, 214, 515. [Google Scholar] [CrossRef]
  13. Wang, D.; Mota, D. 4D Gauss-Bonnet gravity: Cosmological constraints, H0 tension and large scale structure. Phys. Dark Universe 2021, 32, 100813. [Google Scholar] [CrossRef]
  14. Li, B.; Barrow, J.D.; Mota, D.F. Cosmology of modified Gauss-Bonnet gravity. Phys. Rev. D 2007, 76, 044027. [Google Scholar] [CrossRef] [Green Version]
  15. Bahamonte, S.; Bohmer, C.G.; Wright, M. Modified teleparallel theories of gravity. Phys. Rev. D 2015, 92, 104042. [Google Scholar] [CrossRef]
  16. Paliathanasis, A. Cosmological evolution and exact solutions in a fourth-order theory of gravity. Phys. Rev. D 2017, 95, 064062. [Google Scholar] [CrossRef] [Green Version]
  17. Harko, T.; Lobo, F.S.N.; Nojiri, S.; Odintsov, S.D. f(R,T) gravity. Phys. Rev. D 2011, 84, 024020. [Google Scholar] [CrossRef] [Green Version]
  18. Nunes, R.C.; Bonilla, A.; Pan, S.; Saridakis, E.N. Observational constraints on f(T) gravity from varying fundamental constants. EPJC 2017, 77, 230. [Google Scholar] [CrossRef] [Green Version]
  19. Yousaf, Z.; Bamba, K.; Bhatti, M.Z.; Farwa, U. Quasi static evolution of compact objects in modified gravity. Gen. Rel. Grav. 2022, 54, 7. [Google Scholar] [CrossRef]
  20. Farwa, U.; Yousaf, Z.; Bhatti, M.Z. A measure of complexity for axial self-gravitating static fluids. Phys. Scr. 2022, 97, 105307. [Google Scholar] [CrossRef]
  21. Shankaranarayanan, S.; Johnson, J.P. Modified theories of gravity: Why, how and what? Gen. Rel. Grav. 2022, 54, 44. [Google Scholar] [CrossRef]
  22. Jimenez, J.B.; Heisenberg, L.; Koivisto, T.S. The Geometrical Trinity of Gravity. Universe 2019, 5, 173. [Google Scholar] [CrossRef] [Green Version]
  23. Weitzenböck, R. Invarianten Theorie; Nordhoff: Groningen, The Netherlands, 1923. [Google Scholar]
  24. Hayashi, K.; Shirafuji, T. New general relativity. Phys. Rev. D 1979, 19, 3524. [Google Scholar] [CrossRef]
  25. Nester, J.M.; Yo, H.-J. Symmetric teleparallel general relativity. Chin. J. Phys. 1999, 37, 113. [Google Scholar]
  26. Starobinsky, A.A. A new type of isotropic cosmological models without singularity. Phys. Lett. B 1980, 91, 99. [Google Scholar] [CrossRef]
  27. Barrow, J.D. The premature recollapse problem in closed inflationary universes. Nucl. Phys. B 1988, 296, 679. [Google Scholar] [CrossRef]
  28. Sotiriou, T.P.; Faraoni, V. f(R) theories of gravity. Rev. Mod. Phys. 2010, 82, 451. [Google Scholar] [CrossRef] [Green Version]
  29. Oikonomou, V.K.; Giannakoudi, I. A panorama of viable f(R) gravity dark energy models. IJMPD 2022, 31, 2250075. [Google Scholar] [CrossRef]
  30. Cognola, G.; Elizalde, E.; Nojiri, S.; Odintsov, S.D.; Sebastiani, L.; Zerbini, S. Class of viable modified f(R) gravities describing inflation and the onset of accelerated expansion. Phys. Rev. D 2007, 77, 046009. [Google Scholar] [CrossRef] [Green Version]
  31. Bahamonte, S.; Dialektopoulos, K.F.; Escamilla-Rivera, C.; Gakis, V.; Hendry, M.; Said, J.L.; Mifsud, J.; Valentino, E.D. Teleparallel Gravity: From Theory to Cosmology, Report on Progress in Physics. arXiv 2022, arXiv:2106.13793. [Google Scholar]
  32. Ferraro, R.; Fiorini, F. Modified teleparallel gravity: Inflation without an inflaton. Phys. Rev. D. 2007, 75, 084031. [Google Scholar] [CrossRef] [Green Version]
  33. Dent, J.B.; Dutta, S.; Saridakis, E.N. f(T) gravity mimicking dynamical dark energy. Background and perturbation analysis. JCAP 2011, 1, 9. [Google Scholar] [CrossRef] [Green Version]
  34. Krssak, M.; Hoogen, R.J.V.; Pereira, J.G.; Coley, C.G.B.A.A. Teleparallel Theories of Gravity: Illuminating a Fully Invariant Approach. Class. Quantum Grav. 2019, 36, 183001. [Google Scholar] [CrossRef] [Green Version]
  35. Ferraro, R.; Fiorini, F. Born-Infeld gravity in Weitzenböck spacetime. Phys. Rev. D 2008, 78, 124019. [Google Scholar] [CrossRef] [Green Version]
  36. Cai, Y.F.; Chen, S.H.; Dent, J.B.; Dutta, S.; Saridakis, E.N. Matter Bounce Cosmology with the f(T) Gravity. Class. Quantum Grav. 2011, 28, 215011. [Google Scholar] [CrossRef] [Green Version]
  37. Atayde, L.; Frusciante, N. Can f(Q) gravity challenge ΛCDM? Phys. Rev. D 2021, 104, 064052. [Google Scholar] [CrossRef]
  38. Anagnostopoulos, F.K.; Basilakos, S.; Saridakis, E.N. First evidence that non-metricity f(Q) gravity could challenge ΛCDM. Phys. Lett. B 2021, 822, 136634. [Google Scholar] [CrossRef]
  39. Solanki, R.; De, A.; Sahoo, P.K. Complete dark energy scenario in f(Q) gravity. Phys. Dark Universe 2022, 36, 100996. [Google Scholar] [CrossRef]
  40. Arora, S.; Sahoo, P.K. Crossing Phantom Divide in f(Q) gravity. Annalen. Phys. 2022, 534, 2200233. [Google Scholar] [CrossRef]
  41. De, A.; Mandal, S.; Beh, J.T.; Loo, T.-H.; Sahoo, P.K. Isotropization of locally rotationally symmetric Bianchi-I universe in f(Q)-gravity. EPJC 2022, 82, 72. [Google Scholar] [CrossRef]
  42. Lin, R.-H.; Zhai, X.-H. Spherically symmetric configuration in f(Q) gravity. Phys. Rev. D 2021, 103, 124001. [Google Scholar] [CrossRef]
  43. Ambrosio, F.D.; Fell, S.D.B.; Heisenberg, L.; Kuhn, S. Black holes in f(Q) gravity. Phys. Rev. D 2022, 105, 024042. [Google Scholar] [CrossRef]
  44. Khyllep, W.; Paliathanasis, A.; Dutta, J. Cosmological solutions and growth index of matter perturbations in f(Q) gravity. Phys. Rev. D 2021, 103, 103521. [Google Scholar] [CrossRef]
  45. Lymperis, A. Late-time cosmology with phantom dark-energy in f(Q) gravity. JCAP 2022, 11, 18. [Google Scholar] [CrossRef]
  46. Gadbail, G.N.; Mandal, S.; Sahoo, P.K. Parametrization of Deceleration Parameter in f(Q) Gravity. Physics 2022, 4, 1403. [Google Scholar] [CrossRef]
  47. Gadbail, G.N.; Mandal, S.; Sahoo, P.K. Reconstruction of ΛCDM universe in f(Q) gravity. Phys. Lett. B 2022, 835, 137509. [Google Scholar] [CrossRef]
  48. Bajardi, F.; Vernieri, D.; Capozziello, S. Bouncing cosmology in f(Q) symmetric teleparallel gravity. Eur. Phys. J. Plus 2020, 135, 912. [Google Scholar] [CrossRef]
  49. Solanki, R.; Sahoo, P.K. Statefinder Analysis of Symmetric Teleparallel Cosmology. Annalen der Physik 2022, 534, 2200076. [Google Scholar] [CrossRef]
  50. Hu, K.; Katsuragawa, T.; Qiu, T. ADM formulation and Hamiltonian analysis of f(Q) gravity. Phys. Rev. D 2022, 106, 044025. [Google Scholar] [CrossRef]
  51. Jimenez, J.B.; Heisenberg, L.; Koivisto, T.S. Cosmology in f(Q) geometry. Phys. Rev. D 2020, 101, 103507. [Google Scholar] [CrossRef]
  52. Eisenhart, L.P. Non-Riemannian Geometry; American Mathematical Society, Colloquium Publications: New York, NY, USA, 1927; Volume VIII. [Google Scholar]
  53. D’Ambrosio, F.; Heisenberg, L.; Kuhn, S. Revisiting cosmologies in teleparallelism. Class. Quantum Grav. 2022, 39, 025013. [Google Scholar] [CrossRef]
  54. Hohmann, M. General covariant symmetric teleparallel cosmology. Phys. Rev. D 2021, 104, 124077. [Google Scholar] [CrossRef]
  55. Dimakis, N.; Paliathanasis, A.; Roumeliotis, M.; Christodoulakis, T. FLRW solutions in f(Q) theory: The effect of using different connections. Phys. Rev. D 2022, 106, 043509. [Google Scholar] [CrossRef]
  56. Dimakis, N.; Roumeliotis, M.; Paliathanasis, A.; Apostolopoulos, P.S.; Christodoulakis, T. Self-similar cosmological solutions in symmetric teleparallel theory: Friedmann-Lemaître-Robertson-Walker spacetimes. Phys. Rev. D 2022, 106, 123516. [Google Scholar] [CrossRef]
  57. Zhao, D. Covariant formulation of f(Q) theory. Eur. Phys. J. C 2022, 82, 303. [Google Scholar] [CrossRef]
  58. Clifton, T.; Barrow, J.D. The power of general relativity. Phys. Rev. D 2005, 72, 103005. [Google Scholar] [CrossRef] [Green Version]
  59. Esposito, F.; Carloni, S.; Vignolo, S. Bianchi type-I cosmological dynamics in f(Q) gravity. Class. Quantum Grav. 2022, 39, 235014. [Google Scholar] [CrossRef]
  60. McIntosh, C.B.G.; Steel, J.D. All vacuum Bianchi I metrics with a homothety. Class. Quantum Grav. 1991, 8, 1173. [Google Scholar] [CrossRef] [Green Version]
  61. Kasner, E. Geometrical theorems on Einstein’s cosmological equations. Am. J. Math. 1921, 43, 217. [Google Scholar] [CrossRef]
  62. Gasperini, M. The Kasner Solution. In Theory of Gravitational Interactions; Undergraduate Lecture Notes in Physics; Springer: Milano, Italy, 2013. [Google Scholar]
  63. Paliathanasis, A.; Said, J.L.; Barrow, J.D. Stability of the Kasner Universe in f(T) Gravity. Phys. Rev. D 2018, 97, 044008. [Google Scholar] [CrossRef] [Green Version]
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

Paliathanasis, A. Stability Properties of Self-Similar Solutions in Symmetric Teleparallel f(Q)-Cosmology. Symmetry 2023, 15, 529. https://doi.org/10.3390/sym15020529

AMA Style

Paliathanasis A. Stability Properties of Self-Similar Solutions in Symmetric Teleparallel f(Q)-Cosmology. Symmetry. 2023; 15(2):529. https://doi.org/10.3390/sym15020529

Chicago/Turabian Style

Paliathanasis, Andronikos. 2023. "Stability Properties of Self-Similar Solutions in Symmetric Teleparallel f(Q)-Cosmology" Symmetry 15, no. 2: 529. https://doi.org/10.3390/sym15020529

APA Style

Paliathanasis, A. (2023). Stability Properties of Self-Similar Solutions in Symmetric Teleparallel f(Q)-Cosmology. Symmetry, 15(2), 529. https://doi.org/10.3390/sym15020529

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop