Next Article in Journal
Dark Energy Is the Cosmological Quantum Vacuum Energy of Light Particles—The Axion and the Lightest Neutrino
Next Article in Special Issue
Running Vacuum in the Universe: Phenomenological Status in Light of the Latest Observations, and Its Impact on the σ8 and H0 Tensions
Previous Article in Journal
Decoupled Anisotropic Solutions Using Karmarkar Condition in f(G, T) Gravity
Previous Article in Special Issue
Generic Modification of Gravity, Late Time Acceleration and Hubble Tension
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamical Systems Analysis of f(Q) Gravity

1
Department of Mathematics, University College London, Gower Street, London WC1E 6BT, UK
2
Department of Physics and EHU Quantum Center, Faculty of Science and Technology, University of the Basque Country, 48080 Bilbao, Spain
*
Author to whom correspondence should be addressed.
Universe 2023, 9(4), 166; https://doi.org/10.3390/universe9040166
Submission received: 2 February 2023 / Revised: 27 February 2023 / Accepted: 7 March 2023 / Published: 29 March 2023
(This article belongs to the Special Issue Modified Gravity Approaches to the Tensions of ΛCDM)

Abstract

:
Modified gravity theories can be used for the description of homogeneous and isotropic cosmological models through the corresponding field equations. These can be cast into systems of autonomous differential equations because of their sole dependence on a well-chosen time variable, be it the cosmological time, or an alternative. For that reason, a dynamical systems approach offers a reliable route to study those equations. Through a model-independent set of variables, we are able to study all f ( Q ) modified gravity models. The drawback of the procedure is a more complicated constraint equation. However, it allows the dynamical system to be formulated in fewer dimensions than using other approaches. We focus on a recent model of interest, the power-exponential model, and generalize the fluid content of the model.

1. Introduction

We have an unprecedented understanding of the gravitational interaction as the main actor in the large-scale dynamics of the Universe, being responsible for the formation and evolution of structures on the largest scales. Einstein’s theory of General Relativity (GR) accounts successfully for a vast array of gravitational phenomena [1,2,3]. Unfortunately, the so-called dark sector, i.e., dark matter and dark energy, represents a true challenge to an otherwise successful paradigm. The aim of solving these problems has motivated the consideration of slight modifications of GR compatible with observations. The geometry of the Universe is assumed to be the spatially flat Friedmann–Lemaître—Robertson–Walker (FLRW) line element, in agreement with most current observational data [2,4,5]. The FLRW model appears to be the best cosmological model available at the moment [6,7,8,9,10,11,12,13,14,15,16]. It takes the form d s 2 = N ( t ) 2 d t 2 + a ( t ) 2 δ i j d x i d x j with scale factor a ( t ) and arbitrary lapse function N ( t ) . In what follows we can set N ( t ) = 1 without loss of generality, however, we note that some modified theories of gravity might not be compatible with this choice. Modifications of GR can help address observational tensions on the expansion rate of the universe as given by H 0 , and the S 8 parameter characterizing linear matter fluctuations on the scale of 8 h−1 [17,18,19].
Modified theories of gravity have been studied for almost as long as GR itself [20,21]. Possible modifications may bring extra geometrical structures, increase the number of dimensions, or introduce non-linearities into the Einstein–Hilbert action which is linear in curvature. Other options include non-minimal matter curvature couplings. The vast majority of modified gravity models being considered to fall somewhere in the above description [22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37].
Our work will deal with field equations that contain derivatives no higher than second order with respect to the independent variables. Due to the specific structure of the models, one can introduce the Hubble function so that all field equations take the form E i ( H , H ˙ , Ψ ) = 0 where we use Ψ to symbolize the presence of any matter fields. The rather simple description covers many (if not all) second-order modified gravity theories, provided non-minimal couplings are excluded. Let us also mention that we assume the usual matter conservation equations hold for each fluid individually. This is an assumption on the modified gravity theory in question and is linked to diffeomorphism invariance, and most models satisfy this assumption.
It is useful to cast the cosmological field equations into the form of a dynamical system, for a comprehensive summary of past work in the field see [38] and references therein. We briefly mention that the choice of variables can be problematic, something that is well-known in f ( R ) gravity, see for instance [39,40,41,42], see also [43] for f ( T ) gravity. Similar to these approaches we also tend to find complicated constraint equations when the most useful variables are chosen. However, it is quite remarkable that we can make a significant number of general statements about the system for arbitrary models. For example, all possible de Sitter points can be found by the introduction of convenient functions [44].
In the following, we deal with the class of modified gravity theories referred to generally as symmetric teleparallel gravity, or f ( Q ) gravity [45,46,47], which have similarities with the f ( G ) and f ( T ) theories, see [35,36,45,48,49,50] for details. In fact, the geometric scalars of these theories all coincide in cosmology Q = T = G = 6 H 2 . The cosmological field equations of these theories read
6 f H 2 1 2 f = ρ ,
( 12 H 2 f + f ) H ˙ = 1 2 ( ρ + p ) ,
and can be written as a dynamic system using standard techniques. Our formulation allows us a simple comparison with the successful Λ CDM model. Its early-time behavior is dominated by radiation which is a saddle (or repeller) from the dynamical systems viewpoint, whereas the late-time asymptotic regime represents a (de Sitter) cosmological constant-dominated attractor. The matter dominated epoch is a transient situation (saddle point). The specific f ( Q ) model we will study here displays properties compatible with the Λ CDM model.

2. Brief Review on the Standard Approach to the Construction of a Dynamical System

It is rarely easy to find appropriate variables to formulate modified gravity models as dynamical systems, again see [38]. Motivations to introduce different variables can vary, leading to distinctly different features [39,40,41]. Here we briefly run through the standard approach and show how it is equivalent to our formulation in one fewer dimensions.
A typical approach takes the Friedmann Equation (1) and divides by 6 f H 2 (where the assumption f 0 must be made). This excludes trivial constant functions, but also the case where f passes through zero dynamically, which may be important. The typical variables are
X i = ρ i 6 f H 2 , Y = f 12 f H 2 ,
where we are considering various matter-energy sources that may exist. As a reminder, General Relativity of course corresponds to f ( Q ) = Q + 2 λ 0 , where λ 0 stands for a cosmological constant term (with appropriate dimensions). With these variables, the Friedmann equation becomes the following constraint
i X i + Y = 1 .
The naturalness of the Friedmann constraint with this choice of variables is immediately apparent.
The next step is to introduce an additional variable to remove the explicit dependence on the Hubble function which would appear in the dynamical equations. A convenient choice for this extra variable is
Z = Q / Q 0 1 + Q / Q 0 = H 2 / H 0 2 1 + H 2 / H 0 2 ,
and clearly Z is positive and smaller than 1 for all times (that is Q = 6 H 2 is non-negative and finite).
The constraint equation allows us to remove one of the independent variables so that we are left with as many independent variables as one plus the number of matter sources. However, this is not the end of the story: since f is a function of Q it can always be rewritten in terms of Z for some given function. Therefore, Y can be removed altogether. This is exactly the approach we will take in the following section. One can think of Y as encapsulating information about the free parameters that can be found in f (such as the important case of a cosmological constant) but which are not depicted by Z. Further details on the meaning of this will be discussed with a specific example.
In the usual fashion, one then uses the time variable N = log a . We will not produce the explicit form of the dynamical equations d X i / d N and d Y / d N , which in the of for a single matter source can be found in our work [35]. It can be seen that defining the function
m = 12 H 2 ( log f ) ,
allows us to make very general statements about the stability of the system using Jacobian N equations.
We mention again the possibility of taking advantage of the dependence between variables Y and Z to remove one extra dimension from the system, see also [43]. To this end, we use the fact that Y = f / ( 2 f Q ) , which by virtue of the chain rule along with Equation (5) and the constraint i X i + Y = 1 allows us to write
i X ˙ i = Y Q ( Q 0 + Q ) 2 Q 0 Z ˙ .
We must remember throughout that Q can always be recast as a function of Z. The next step is to use (7) to replace X ˙ i in the equation where it appears. Next, the constraint i X i + Y = 1 must be used to remove X i and then it is necessary to remember that Y is a function of Z. This gives us a system with one fewer dimension and slightly different equations for the evolution of the dynamical variables. It is not difficult to see the agreement between the two approaches (with different dimensionalities), but again, we suggest the reader to check [44].
It may at first appear that this equivalent formulation poses no benefits, as a model needs to be specified in order to extract any useful information from the system. This, however, is incorrect, and we shall see that even in the completely general case we can analyze the dynamical system to some degree. Moreover, phase space and stability analysis become much simpler in fewer dimensions, despite representing the same physics.

3. Dynamical Systems Formulation

3.1. General Setup with Two Fluids

Generalizing the formulation of the problem with a reduced dimensionality we can now tackle a two-fluid case in just two dimensions. Their energy densities will be ρ 1 and ρ 2 with the equation of state parameters w 1 and w 2 . For this discussion we leave these equations of state parameters arbitrary, though we will later set w 1 = 0 and w = 1 / 3 . These choices render the cosmological equations as
f 6 H 2 ρ 1 3 H 2 ρ 2 3 H 2 + 2 f = 0 ,
( w 1 + 1 ) ρ 1 + ( w 2 + 1 ) ρ 2 + 2 H ˙ ( f + 12 H 2 f ) = 0 .
We now define the following dynamical variables
X i = ρ i 3 H 2 , for i = 1 , 2 ,
Z = H 2 / H 0 2 1 + H 2 / H 0 2 .
The first two variables are non-negative, have an easily recognizable form (being simply the standard matter density parameters Ω i ) and are different from their previous counterparts (3) as f does not appear in the denominator. There is now no dependence on f in any of the variables. The reason for this is the ability to write any function of H as a function of Z, with its dynamics being determined by the Friedmann constraint itself.
It is now possible to cast the Friedmann constraint as an expression that involves only the new variables
X 1 + X 2 = 1 1 Z f 6 H 0 2 + 2 f ,
where we are treating f = f ( 6 H 0 2 Z / ( 1 Z ) ) as an arbitrary function of Z.
A dual interpretation of f and f is always possible in the sense that they can be seen as functions of the scalar that governs the modified theory of gravity or as functions of H. Whatever the case, we will always present them as functions of Z. As a consequence the equation for X 1 + X 2 is also a function of Z until any particular form of f is specified. We can therefore use this equation to eliminate either X 1 or X 2 , making the phase space two-dimensional.
Let us choose to eliminate X 1 and consider the evolution of { X 2 , Z } ,
d X 2 d N = 3 X 2 m + 1 ( w 2 w 1 ) X 2 f ( w 2 + 1 ) m ( w 1 + 1 ) n + 2 w 1 w 2 + 1 ,
d Z d N = 3 ( Z 1 ) Z ( ( w 1 + 1 ) ( n 2 ) + ( w 1 w 2 ) X 2 / f ) m + 1 ,
where we have taken advantage of (12) and introduced the convenient functions
m ( Z ) : = 2 Q f f = 12 H 0 2 Z f ( 1 Z ) f ,
n ( Z ) : = f Q f = f ( 1 Z ) 6 H 0 2 Z f .
It is also worth keeping in mind that f is dimensionless whereas f has units of H 0 2 because f has dimensions H 0 2 . For this reason, Equations (13) and (14) and m ( Z ) and n ( Z ) are dimensionless as well.
Once a specific theoretical setting is chosen through f we are left with a closed system of equations ready to be studied. Fortunately, some of the key features of this set of equations do not depend on the chosen form of f, so a number of very broad conclusions may be drawn, which adds to the interest of our analysis and approach.

3.2. Fixed Points

For this discussion we assume w 1 w 2 and ( w 1 , w 2 ) 1 , that is, two different fluids and neither of them is a cosmological constant. There are then two families of fixed points for the system, which we will look at individually.
The first one corresponds to the points { X 2 , Z } = { 0 , Z } , where the second coordinate is specified through solutions of the algebraic equation
( n ( Z ) 2 ) ( Z 1 ) Z 1 + m ( Z ) = 0 .
Note that the locations of these points are independent of both fluid parameters w 1 and w 2 . The solutions where n ( Z ) = 2 and m ( Z * ) with n ( Z * ) finite will be particularly important. As such, we will name these points Pn and Pm, respectively. In the following section we show that these points possess very particular fixed properties that are of interest when assessing the validity of different cosmological models.
The second family is characterised by the two points B = { X 2 , 1 } and C = { X 2 , 0 } with
X 2 = f w 2 w 1 ( ( w 1 + 1 ) n ( Z ) + w 2 + ( w 2 + 1 ) m ( Z ) 1 2 w 1 ) ,
where X 2 * is evaluated at Z = 1 and Z = 0 . Here the location does depend on the particular values of w 1 and w 2 , and X 2 * should be treated as a function of Z. Note that the existence of the points and their belonging to the physical phase space is not guaranteed. To assess this one must study the Hubble constraint for that particular model, which can be written with the help of n ( Z ) as
X 1 = f ( 2 n ( Z ) ) X 2 .
We require that both variables X 1 and X 2 be positive in order to satisfy energy conditions, along with 0 Z 1 .
Lastly, we note the possibility that the dynamical equations diverge for some particular value of Z, which occurs if m ( Z ) = 1 or f = 0 . In these cases, trajectories cannot be extended beyond this Z coordinate and the phase space exhibits a ’critical line’ behavior, see [44]. For a full analysis including the existence criteria of the critical points, knowledge of the model f ( Q ) is needed.

3.3. Physical Parameters of the General System

The deceleration parameter q and the effective equation of state w eff can be expressed in terms of the dynamical variables
q : = a ¨ a a ˙ 2 = 1 3 2 X 1 + w 1 X 1 + X 2 + w 2 X 2 n ( Z ) 2 ( X 1 + X 2 ) ( m ( Z ) + 1 ) ,
w eff : = p tot ρ tot = 1 X 1 + w 1 X 1 + X 2 + w 2 X 2 n ( Z ) 2 ( X 1 + X 2 ) ( m ( Z ) + 1 ) ,
where the total energy density ρ tot is defined as ρ tot = ρ 1 + ρ 2 + ρ f with ρ f representing the additional non-GR terms
ρ f : = 3 H 2 + 1 2 f 6 H 2 f .
The total pressure is defined similarly p tot = p 1 + p 2 + p f . In Equations (20) and (21) the variable X 1 can equally be rewritten in terms of X 2 and Z using the Friedmann constraint, but a full analysis cannot be carried out until a function f is specified.
The density parameter of the additional non-GR terms will turn out to be useful later on, which we can express in terms of our dynamical variables as
Ω f : = ρ f 3 H 2 = 1 f ( 2 n ( Z ) ) ,
which satisfies Ω 1 + Ω 2 + Ω f = 1 from the Friedmann equation. Hence we have obtained a very neat expression representing the contributions of the modified theory beyond GR. As previously mentioned, a more standard dynamical systems formulation would introduce a dynamical variable for this Ω f (e.g., the Y in Section 2) but because it can be written totally in terms of Z this is not necessary. We therefore obtain a phase space with fewer dimensions at the expense of a more cumbersome Hubble constraint.
For the fixed points Pn satisfying n ( Z ) = 2 one immediately has Ω f = 1 . Similarly for points Pm, looking at the definition of m ( Z ) in (15), we see that if m ( Z ) whilst n ( Z ) stays finite, we also obtain Ω f = 1 . For these two solutions, the deceleration parameter and equation of state are fixed to be q = w eff = 1 , representing a de Sitter Universe. In fact, this is a necessary requirement for any de Sitter solution given that we have assumed w 1 w 2 and w 1 , w 2 1 . Models where n ( Z ) 2 and m ( Z ) for some Z in the range ( 0 , 1 ) cannot possess a de Sitter fixed point. This immediately rules out models such as f ( Q ) Q α for α 1 / 2 , or f ( Q ) = Q + β Q 2 for β 0 .
For the other fixed points B and C at { X 2 , 1 } and { X 2 , 0 } , we can also determine the deceleration parameter and effective equation of state. For both of these points one obtains w eff = w 2 and q = ( 1 + 3 w 2 ) / 2 , which is a remarkably general result independent from the model. These results are summarized in Table 1.
Note that there may exist other fixed points at { 0 , Z } , solutions to (17), which have not been included in the table. This is because their properties are more dependent on the specific model and do not lead to fixed values of q or w eff . Furthermore, note that we have not yet fully discussed the existence conditions for the fixed points, and their presence in the physical phase space depends on the Hubble constraint for that particular model. We now move on to studying a chosen model where a full analysis can be carried out.

4. Applications to f ( Q ) Models

4.1. Anagnostopoulos et al. Model

The power-exponential model proposed by Anagnostopoulos et al. in [51] displays a number of interesting features and was shown to pass a variety of observational tests [51,52]. In particular, the authors studied the model against Supernovae type Ia (SNIa), Baryonic Acoustic Oscillations (BAO), cosmic chronometers (CC), and Redshift Space Distortion (RSD) data and found that it is comparable, and for some datasets preferable, to the Λ CDM model. Moreover, it immediately passes early universe constraints. As such, it has been shown to be a genuine alternative to the Λ CDM concordance model and worthwhile studying from a dynamical systems perspective.
The model is given by the function
f ( Q ) = Q e λ Q 0 Q ,
with the single free parameter λ . A dynamical systems analysis was recently performed for this model in [53]. There the authors studied the background and perturbation equations of a universe with a single fluid matter component ( w = 0 ) , and the subsequent phase space was three-dimensional. It is interesting to then study this model in our reduced dimensionality formulation with an additional matter fluid component, which will turn out to be two-dimensional. Moreover, the reduced dimensions in our approach will turn out to make the stability analysis much simpler to compute.
In the limit that λ vanishes the model (24) reduces to GR without a cosmological constant. It does not however have a direct Λ CDM limit. When the parameter λ is small, to first order the function behaves like GR with a cosmological constant term Q 0 λ , and this behavior will be observed in the phase space analysis. The sign of the parameter λ leads to different phase spaces, and so we will investigate both cases. We will also assume that λ 0 , as this trivially leads back to GR.
The system is described by the dynamical Equations (13) and (14) along with Hubble constraint (19), with the functions n ( Z ) and m ( Z ) taking a remarkably simple form
n ( Z ) = Z Z + λ ( Z 1 ) ,
m ( Z ) = 2 ( Z 1 ) 2 λ 2 Z ( Z + λ ( Z 1 ) ) .
The f term written explicitly in terms of the variable Z is
f ( Q ) = ( Z + λ ( Z 1 ) ) e λ ( Z 1 ) / Z Z .
The fixed points are solutions to Equations (17) and (18), which can be easily solved using the exact forms of n ( Z ) , m ( Z ) and f given above. The first family of solutions along the X 2 = 0 line with Z = Z * are the points A = {0,1}, Pm = {0,0} and P n = { 0 , 2 λ / ( 1 + 2 λ ) } . Point A is an additional solution to the algebraic Equation (17) with properties that could not be determined in general, therefore it was left out of Table 1. The critical point P m satisfies m ( Z ) with n ( Z ) = finite , whilst the point P n is a solution to n ( Z ) = 2 . Hence these two points describe de Sitter attractors, as explained in the previous section and in Table 1.
The second set of solutions from Equation (18) include the point B = { 1 , 1 } and a conditional point C at { 0 , 0 } which requires λ < 0 . However, we will ignore this final point because it coincides with P m . Point B is the ρ 2 matter-dominated point. Note again that we have not yet assessed the validity of any of the fixed points; only those satisfying ( X 1 , X 2 ) 0 and 1 Z 0 are physically meaningful, for which we will need to use the Hubble constraint.
The Hubble constraint (19) can be written explicitly in terms of the variables as
X 1 + X 2 = e λ Z 1 / Z 1 + 2 λ ( Z 1 ) Z .
The requirement that our matter fluids have positive energy density leads to physical bounds on the phase space. In particular, one notes that for positive λ our fluid density parameters take the maximum value of one, whilst for negative λ we instead obtain ( X 1 , X 2 ) 2 e 1.21 . This situation, where the density parameters can be greater than one, can be understood in physical terms by considering the modified density parameter Ω f in Equation (23). For positive λ the density parameter is non-negative, and from the Hubble equation Ω 1 + Ω 2 + Ω f = 1 we can conclude that Ω 1 + Ω 2 1 . However, for negative λ we instead have the minimum of Ω f = 1 2 e , which leads to Ω 1 + Ω 2 2 e .
Using (28) we can determine X 1 at each of the fixed points, as well as the conditions for the point to be part of the physical phase space. The points A and B are always present irrespective of λ . The de Sitter points P m and P n require λ < 0 and λ > 0 respectiely. It is also interesting to note that all of the fixed points, their locations and their existence criteria are independent of the specific fluid equation of state. The deceleration parameter and effective equation of state can be evaluated at each of the fixed points using Equations (20) and (21). Lastly, linear stability theory has been applied to the fixed points, which can be found in Appendix A. In this formulation points A, B and P n are hyperbolic and P m is nonhypebolic. However, we show in Appendix A that P m acts as the late-time attractor within the physical phase space ( λ < 0 ). These results are collated in Table 2.

4.2. Phase Space Analysis

Next, we will fix the equations of state of our two fluid components to be w 1 = 0 and w 2 = 1 / 3 , representing matter and radiation. The phase portraits for this model with positive and negative values of the free parameter λ are shown in Figure 1. The absolute value of the free parameter is chosen to be | λ | = 0.371 , as this was shown in [51] to be the most consistent with observational constraints. The bordered region highlights the physical phase space and the red overlay represents regimes where the expansion of the Universe is accelerating q < 0 .
It is somewhat expected that for positive λ the qualitative results should be the same as GR with a positive cosmological constant, as the dynamics are similar to Λ CDM for this parameter value [51,53]. Comparing the phase portrait in Figure 1a to that of GR with a cosmological constant, which can be found in our previous work using the same dynamical systems formulation [44], reveals that the qualitative features are indeed the same. The phase space contains an early radiation-dominated repeller, point B, a matter dominate saddle, point A, and a late-time de Sitter attractor, point P n . Stability analysis indeed verifies that point B is unstable, A is a saddle and P n is stable.
For the case of a negative parameter λ , see Figure 1b, the phase space is similar but distinctly different. The fixed points of the system remain the same except point P m at { 0 , 0 } replaces P n at { 0 , 2 λ / ( 1 + 2 λ ) } . The new de Sitter point P m possess the same properties as P n , refer to Table 1 and Table 2. The obvious difference is that the physical phase space extends outwards beyond X 2 = 1 . As previously explained, this is due to the modified density parameter Ω f having a negative lower bound for λ < 0 . This leads to a noticeably different evolution of the density parameters and physical parameters q and w eff .
In Figure 2, the evolution of the matter and radiation density parameters Ω m = X 1 , Ω r = X 2 , the deceleration parameter q and the effective equation of state w eff are shown for both phase spaces. Figure 2a shows the evolution for λ > 0 of a trajectory following a heteroclinic orbit from points B → A → P n , whilst Figure 2b follows the heteroclinic orbit B → A → P m for λ < 0 . The dashed lines represent the evolution of these parameters for GR with a positive cosmological constant, f ( Q ) = Q + 2 Λ Q 0 = 6 H 2 + 12 H 0 2 Λ with Λ = | λ | . We have chosen Λ to equal λ such that the de Sitter point P n of GR and the de Sitter point P n of the Anagnostopoulos model (with positive λ ) take the same values. It is interesting to note that for f ( Q ) = Q + 2 Λ Q 0 the fixed points are exactly the A, B and P n given in Table 2 with λ replaced by Λ . For the case where the free parameter λ is negative, the same matching cannot be done because the point P m does not exist for GR with a cosmological constant. This can be seen from the function n ( Z ) = 1 2 Λ + 2 Λ Z and m ( Z ) = 0 . The point P m cannot exist because it requires m ( Z ) , see [44] for more details.
The background evolution of the model with positive λ can be seen to match very closely with its GR counterpart, given the same initial conditions. For the negative λ case in Figure 2b, a sharp spike in the matter density parameter can be noted before approaching the de Sitter point with w eff = 1 .
Overall, the dynamical system analysis gives a good understanding of the background dynamics of the Anagnostopoulos et al. model for both signs of the free parameter λ . In the positive λ case, the fixed points and stability of the system match what is found for GR with a positive cosmological constant. In fact, from a qualitative point of view, these models are identical. For the negative λ case, the fixed points of the system have the same properties but it is interesting to note the different physical phase space as well as the different evolutions for the density parameters.

5. Summary

We have presented a dynamical systems formulation that is well suited to the cosmological equations arising in various modified gravity theories, with f ( Q ) being the focus of this work. Once a model function f has been specified, the drawback of a more complicated constraint equation is indeed present but of little significance. This is especially true due to the clarity gained when dealing with a two-dimensional as opposed to a three-dimensional phase space (and similarly for higher-dimensional analogs). This same approach can easily be generalized to include additional degrees of freedom represented by additional dynamical variables. For example, the inclusion of extra matter sources, scalar fields, or non-zero spatial curvature could be easily realized within our formulation, see for instance [38]. The use of the functions m ( Z ) and n ( Z ) introduced in Equations (15) and (16) is particularly adept at assessing the validity of models, as the existence of late-time de Sitter points can be established almost immediately.
For the recently proposed f ( Q ) model by Anagnostopoulos et al. [51], we analyzed the phase space for a universe comprised of two fluid components, matter and radiation. This complements the previous dynamical systems analysis performed on this model for matter and matter perturbations for a positive value of the free parameter [53]. Indeed, for positive λ we reproduced the dynamics of Λ CDM. This is perhaps to be expected from the series expansion of the function f for small λ , with the leading order terms being Q + Q 0 λ , which is exactly the Lagrangian of GR plus a cosmological constant.
However, we also find the surprising result that a negative value of the parameter leads to a qualitatively similar dynamical system. In particular, it is interesting to note that a negative value of the parameter λ in fact still acts as a positive cosmological constant, leading to a late-time de Sitter point within the phase space. This was determined by studying the functions n ( Z ) and m ( Z ) , Equations (25) and (26), which took a remarkably simple form for this model.
The results of the dynamical systems analysis for the Anagnostopoulos et al. [51] model show that at the background level, it passes cosmological observational constraints, displaying the correct evolutionary behavior of the matter density parameters and the effective equation of state. Namely, for any non-zero value of the parameter λ , there exists an early-time radiation-dominated point, a matter saddle, and an accelerating de Sitter attractor. This study, the first to use both matter and radiation sources, gives more reason to continue to investigate this model in the future.
In summary, the approach taken leads to a number of model-independent results, which would be especially interesting to investigate in more detail. The moral behind the approach can ultimately be traced to the dynamical systems formulations of GR: for each matter source ρ i we introduce the corresponding density parameter Ω i as a variable. We then introduce one additional variable related to the remaining terms in the Hubble constraint in order to close the system, for which we chose the Hubble function H. In second order modifications, such as f ( T ) and f ( Q ) gravity, we have shown that this same prescription works for all models. This is in contrast to most of the dynamical systems formulations used in modified gravity. In the future it would be interesting to further study promising alternatives to the Λ CDM model using such a formulation. It would also be interesting to search for models that satisfy current observations yet exhibit a different and more complex fixed point behavior to GR, which could lead to qualitatively different predictions.

Author Contributions

Conceptualization, C.B. and R.L.; methodology, validation, formal analysis and investigation, writing—original draft preparation, writing—review and editing, C.B., R.L. and E.J.; visualization, C.B. and E.J.; funding acquisition, R.L. and E.J. All authors have read and agreed to the published version of the manuscript.

Funding

R.L. was supported by the Spanish Ministry of Science and Innovation through research projects PID2021-123226NB-I00 (comprising FEDER funds), and also by the Basque Government and Generalitat Valenciana through research projects IT1628-22 and PROMETEO/2020/079, respectively. Erik Jensko is supported by EPSRC Doctoral Training Programme (EP/R513143/1).

Data Availability Statement

No data were created.

Acknowledgments

R.L. thanks José Beltrán Jiménez for conversations.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A. Stability Analysis of Anagnostopoulos et al. Model

Here we apply linear stability theory to each of the fixed points in Table 2 for the model f ( Q ) = Q exp ( λ Q / Q 0 ) considered in Section 4. Where linear stability theory fails, we look to see what can be said about the nature of the fixed points by examining the autonomous equations and constraints directly.
For point A { 0 , 1 } we obtain the eigenvalues
λ 1 A = 3 ( w 1 w 2 ) , λ 2 A = 3 ( 1 + w 1 ) .
The point is therefore a saddle for w 1 < w 2 or unstable for w 1 > w 2 , as the second eigenvalue is always positive due to the assumption that w 1 > 1 . Point B { 1 , 1 } has eigenvalues
λ 1 B = 3 ( w 2 w 1 ) = λ 1 A , λ 2 B = 3 ( 1 + w 2 ) ,
which is unstable for w 1 < w 2 or a saddle point for w 1 > w 2 . Again, the second eigenvalue is always positive. Due to the freedom in the ordering of our matter fluids ρ 1 and ρ 2 , we choose w 1 < w 2 without loss of generality such that point A is a saddle and point B is unstable.
Point P n { 0 , 2 λ / ( 1 + 2 λ ) } has eigenvalues
λ 1 P n = 3 ( 1 + w 1 ) = λ 2 A , λ 2 P n = 3 ( 1 + w 2 ) = λ 2 B ,
and is therefore always stable. This is the stable late-time de Sitter point of the system.
Point P m has eigenvalues
λ 1 P m = 0 , λ 2 P m = 3 ( 1 + w 2 ) = λ 2 B ,
therefore methods beyond linear stability theory must be used to fully determine the stability.
A closer look at the nonhyperbolic point is shown in Figure A1, as well as the physically allowed values of X 2 and Z as determined from the Hubble constraint (28). Recall that we require λ < 0 for the existence of this point. Despite the fact that the point appears to be mathematically unstable, with trajectories moving away from the point in Figure A1, the Hubble constraint can be used to determine the fate of trajectories within the physical phase space. The boundary of the physical phase space is described by the equations
X 2 = e λ 1 Z 1 1 + 2 λ ( Z 1 ) Z , with X 2 0 , 0 Z 1 , λ < 0 .
As Z approaches zero (from above) X 2 goes to zero. One can also show that for all X 2 0 trajectories always travel in the negative Z direction whilst Z is between 0 and 1. This can be most easily seen by substituting the expression for X 2 on the physical boundary (A5) into the autonomous equation d Z / d N . This resulting equation is
d Z d N = 3 ( 1 + w 2 ) Z 2 ( 1 Z ) Z 2 λ ( 1 Z ) Z 2 + 2 λ 2 ( Z 1 ) 2 λ Z ( 1 Z ) ,
where the equation of state w 2 > 1 . All terms in the numerator and denominator are positive for λ < 0 , therefore d Z / d N is negative and trajectories on the boundary approach the origin.
Following the same logic, the same result can be shown for the general equation d Z / d N with λ < 0 . We can therefore conclude that all physical trajectories satisfying the Hubble constraint travel towards and terminate at the origin, point P m . This is because trajectories do not cross the boundary, and must end at Z = 0 which is only allowed at X 2 = 0 . This indeed matches what can be seen from the phase portraits, Figure 1b and Figure A1, and the numerical solutions in Figure 2b.
Figure A1. Nonhyperbolic fixed point P m .
Figure A1. Nonhyperbolic fixed point P m .
Universe 09 00166 g0a1

References

  1. Abbott, B.P.; Abbott, R.; Abbott, T.D.; Abernathy, M.R.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R.X.; et al. Observation of Gravitational Waves from a Binary Black Hole Merger. Phys. Rev. Lett. 2016, 116, 061102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Aghanim, N.; Akrami, Y.; Ashdown, M.; Aumont, J.; Baccigalupi, C.; Ballardini, M.; Banday, A.J.; Barreiro, R.B.; Bartolo, N.; Basak, S.; et al. Planck 2018 results. VI. Cosmological parameters. Astron. Astrophys. 2020, 641, A6, Erratum in: Astron. Astrophys. 2021, 652, C4. [Google Scholar] [CrossRef] [Green Version]
  3. Will, C.M. Theory and Experiment in Gravitational Physics; Cambridge University Press: Cambridge, UK, 2018; ISBN 978-1-108-67982-4/978-1-107-11744-0. [Google Scholar]
  4. Saadeh, D.; Feeney, S.M.; Pontzen, A.; Peiris, H.V.; McEwen, J.D. How isotropic is the Universe? Phys. Rev. Lett. 2016, 117, 131302. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Efstathiou, G.; Gratton, S. The evidence for a spatially flat Universe. Mon. Not. Roy. Astron. Soc. 2020, 496, L91–L95. [Google Scholar] [CrossRef]
  6. Copeland, E.J.; Sami, M.; Tsujikawa, S. Dynamics of dark energy. Int. J. Mod. Phys. D 2006, 15, 1753–1936. [Google Scholar] [CrossRef] [Green Version]
  7. Jain, B.; Zhang, P. Observational Tests of Modified Gravity. Phys. Rev. D 2008, 78, 063503. [Google Scholar] [CrossRef] [Green Version]
  8. Lombriser, L.; Lima, N.A. Challenges to Self-Acceleration in Modified Gravity from Gravitational Waves and Large-Scale Structure. Phys. Lett. B 2017, 765, 382–385. [Google Scholar] [CrossRef]
  9. Koyama, K. Cosmological Tests of Modified Gravity. Rep. Prog. Phys. 2016, 79, 046902. [Google Scholar] [CrossRef] [Green Version]
  10. Nunes, R.C.; Pan, S.; Saridakis, E.N. New observational constraints on f(T) gravity from cosmic chronometers. JCAP 2016, 08, 011. [Google Scholar] [CrossRef] [Green Version]
  11. Koyama, K. Gravity beyond general relativity. Int. J. Mod. Phys. D 2018, 27, 1848001. [Google Scholar] [CrossRef] [Green Version]
  12. Lombriser, L. Parametrizations for tests of gravity. Int. J. Mod. Phys. D 2018, 27, 1848002. [Google Scholar] [CrossRef] [Green Version]
  13. Lazkoz, R.; Lobo, F.S.N.; Nos, M.O.; Salzano, V. Observational constraints of f(Q) gravity. Phys. Rev. D 2019, 100, 104027. [Google Scholar] [CrossRef] [Green Version]
  14. Benetti, M.; Capozziello, S.; Lambiase, G. Updating constraints on f(T) teleparallel cosmology and the consistency with Big Bang Nucleosynthesis. Mon. Not. Roy. Astron. Soc. 2020, 500, 1795–1805. [Google Scholar] [CrossRef]
  15. Braglia, M.; Ballardini, M.; Finelli, F.; Koyama, K. Early modified gravity in light of the H0 tension and LSS data. Phys. Rev. D 2021, 103, 043528. [Google Scholar] [CrossRef]
  16. Valentino, E.D.; Mena, O.; Pan, S.; Visinelli, L.; Yang, W.; Melchiorri, A.; Mota, D.F.; Riess, A.G.; Silk, J. In the realm of the Hubble tension—a review of solutions. Class. Quant. Grav. 2021, 38, 153001. [Google Scholar] [CrossRef]
  17. Abdalla, E.; Abellán, G.F.; Aboubrahim, A.; Agnello, A.; Akarsu, O.; Akrami, Y.; Alestas, G.; Aloni, D.; Amendola, L.; Anchordoqui, L.A.; et al. Cosmology intertwined: A review of the particle physics, astrophysics, and cosmology associated with the cosmological tensions and anomalies. JHEAp 2022, 34, 49–211. [Google Scholar] [CrossRef]
  18. Valentino, E.D.; Anchordoqui, L.A.; Akarsu, O.; Ali-Haimoud, Y.; Amendola, L.; Arendse, N.; Asgari, M.; Ballardini, M.; Basilakos, S.; Battistelli, E.; et al. Snowmass2021 - Letter of interest cosmology intertwined II: The hubble constant tension. Astropart. Phys. 2021, 131, 102605. [Google Scholar] [CrossRef]
  19. Valentino, E.D.; Anchordoqui, L.A.; Akarsu, Ö.; Ali-Haimoud, Y.; Amendola, L.; Arendse, N.; Asgari, M.; Ballardini, M.; Basilakos, S.; Battistelli, E.; et al. Cosmology Intertwined III: 8 and S8. Astropart. Phys. 2021, 131, 102604. [Google Scholar] [CrossRef]
  20. Goenner, H.F.M. On the history of unified field theories. Living Rev. Rel. 2004, 7, 2. [Google Scholar]
  21. Goenner, H.F.M. On the History of Unified Field Theories. Part II. (ca. 1930–ca. 1965). Living Rev. Rel. 2014, 17, 5. [Google Scholar] [CrossRef] [Green Version]
  22. Capozziello, S. Curvature quintessence. Int. J. Mod. Phys. D 2002, 11, 483–492. [Google Scholar] [CrossRef] [Green Version]
  23. Ferraro, R.; Fiorini, F. Modified teleparallel gravity: Inflation without inflaton. Phys. Rev. D 2007, 75, 084031. [Google Scholar] [CrossRef] [Green Version]
  24. Sotiriou, T.P.; Faraoni, V. f(R) Theories Of Gravity. Rev. Mod. Phys. 2010, 82, 451–497. [Google Scholar] [CrossRef] [Green Version]
  25. Felice, A.D.; Tsujikawa, S. f(R) theories. Living Rev. Rel. 2010, 13, 3. [Google Scholar] [CrossRef] [Green Version]
  26. 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]
  27. Capozziello, S.; Laurentis, M.D. Extended Theories of Gravity. Phys. Rep. 2011, 509, 167–321. [Google Scholar] [CrossRef] [Green Version]
  28. 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]
  29. Clifton, T.; Ferreira, P.G.; Padilla, A.; Skordis, C. Modified Gravity and Cosmology. Phys. Rep. 2012, 513, 1–189. [Google Scholar]
  30. Bamba, K.; Capozziello, S.; Nojiri, S.; Odintsov, S.D. Dark energy cosmology: The equivalent description via different theoretical models and cosmography tests. Astrophys. Space Sci. 2012, 342, 155–228. [Google Scholar] [CrossRef] [Green Version]
  31. Nesseris, S.; Basilakos, S.; Saridakis, E.N.; Perivolaropoulos, L. Viable f(T) models are practically indistinguishable from ΛCDM. Phys. Rev. D 2013, 88, 103010. [Google Scholar] [CrossRef] [Green Version]
  32. Joyce, A.; Jain, B.; Khoury, J.; Trodden, M. Beyond the Cosmological Standard Model. Phys. Rep. 2015, 568, 1–98. [Google Scholar] [CrossRef] [Green Version]
  33. Cai, Y.F.; Capozziello, S.; Laurentis, M.D.; Saridakis, E.N. f(T) teleparallel gravity and cosmology. Rep. Prog. Phys. 2016, 79, 106901. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Nojiri, S.; Odintsov, S.D.; Oikonomou, V.K. Modified Gravity Theories on a Nutshell: Inflation, Bounce and Late-time Evolution. Phys. Rep. 2017, 692, 1–104. [Google Scholar] [CrossRef] [Green Version]
  35. Böhmer, C.G.; Jensko, E. Modified gravity: A unified approach. Phys. Rev. D 2021, 104, 024010. [Google Scholar] [CrossRef]
  36. Saridakis, E.N.; CANTATA. Modified Gravity and Cosmology: An Update by the CANTATA Network; Springer: Berlin/Heidelberg, Germany, 2021; ISBN 978-3-030-83715-0. [Google Scholar]
  37. Böhmer, C.G. Foundations of gravity—Modifications and extensions. In Modified Gravity and Cosmology: An Update by the CANTATA Network; Sari-dakis, E.N., Lazkoz, R., Salzano, V., Moniz, P.V., Capozziello, S., Beltrán Jiménez, J., De Laurentis, M., Olmo, G.J., Eds.; Springer International Publishing: Cham, Switzerland, 2021; pp. 27–38. [Google Scholar] [CrossRef]
  38. Bahamonde, S.; Böhmer, C.G.; Carloni, S.; Copeland, E.J.; Fang, W.; Tamanini, N. Dynamical systems applied to cosmology: Dark energy and modified gravity. Phys. Rep. 2018, 775–777, 1–122. [Google Scholar] [CrossRef] [Green Version]
  39. Amendola, L.; Gannouji, R.; Polarski, D.; Tsujikawa, S. Conditions for the cosmological viability of f(R) dark energy models. Phys. Rev. D 2007, 75, 083504. [Google Scholar] [CrossRef] [Green Version]
  40. Carloni, S. A new approach to the analysis of the phase space of f(R)-gravity. JCAP 2015, 9, 013. [Google Scholar] [CrossRef] [Green Version]
  41. Alho, A.; Carloni, S.; Uggla, C. On dynamical systems approaches and methods in f(R) cosmology. JCAP 2016, 8, 064. [Google Scholar] [CrossRef] [Green Version]
  42. Chakraborty, S.; Dunsby, P.K.S.; Macdevette, K. A note on the dynamical system formulations in f(R) gravity. Int. J. Geom. Meth. Mod. Phys. 2022, 19, 2230003. [Google Scholar] [CrossRef]
  43. Hohmann, M.; Jarv, L.; Ualikhanova, U. Dynamical systems approach and generic properties of f(T) cosmology. Phys. Rev. D 2017, 96, 043508. [Google Scholar] [CrossRef] [Green Version]
  44. Böhmer, C.G.; Jensko, E.; Lazkoz, R. Cosmological dynamical systems in modified gravity. Eur. Phys. J. C 2022, 82, 500. [Google Scholar] [CrossRef]
  45. Jiménez, J.B.; Heisenberg, L.; Koivisto, T.S.; Pekar, S. Cosmology in f(Q) geometry. Phys. Rev. D 2020, 101, 103507. [Google Scholar] [CrossRef]
  46. Jiménez, J.B.; Heisenberg, L.; Koivisto, T.S. The Geometrical Trinity of Gravity. Universe 2019, 5, 173. [Google Scholar] [CrossRef] [Green Version]
  47. Hohmann, M. General covariant symmetric teleparallel cosmology. Phys. Rev. D 2021, 104, 124077. [Google Scholar] [CrossRef]
  48. Jiménez, J.B.; Heisenberg, L.; Koivisto, T. Coincident General Relativity. Phys. Rev. D 2018, 98, 044048. [Google Scholar] [CrossRef] [Green Version]
  49. Harko, T.; Koivisto, T.S.; Lobo, F.S.N.; Olmo, G.J.; Rubiera-Garcia, D. Coupling matter in modified Q gravity. Phys. Rev. D 2018, 98, 084043. [Google Scholar] [CrossRef] [Green Version]
  50. Böhmer, C.G.; Jensko, E. Modified Gravity: A Unified approach to Metric-Affine Models. arXiv 2023. [Google Scholar] [CrossRef]
  51. 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]
  52. Anagnostopoulos, F.K.; Gakis, V.; Saridakis, E.N.; Basilakos, S. New models and big bang nucleosynthesis constraints in f(Q) gravity. Eur. Phys. J. C 2023, 83, 58. [Google Scholar] [CrossRef]
  53. Khyllep, W.; Dutta, J.; Saridakis, E.N.; Yesmakhanova, K. Cosmology in f(Q) gravity: A unified dynamical systems analysis of the background and perturbations. Phys. Rev. D 2023, 107, 044022. [Google Scholar] [CrossRef]
Figure 1. Phase portraits for the Anagnostopoulos et al. model (24). The physical phase space is within the bordered region whilst the red region represents accelerated expansion. (a) Phase space with positive λ . (b) Phase space with negative λ .
Figure 1. Phase portraits for the Anagnostopoulos et al. model (24). The physical phase space is within the bordered region whilst the red region represents accelerated expansion. (a) Phase space with positive λ . (b) Phase space with negative λ .
Universe 09 00166 g001
Figure 2. Evolution of density parameters Ω m , Ω r , effective equation of state w eff , and deceleration parameter q for GR with a positive cosmological constant (dashed) and Anagnostopoulos et al. model (solid line). (a) Evolution for positive λ . (b) Evolution for negative λ .
Figure 2. Evolution of density parameters Ω m , Ω r , effective equation of state w eff , and deceleration parameter q for GR with a positive cosmological constant (dashed) and Anagnostopoulos et al. model (solid line). (a) Evolution for positive λ . (b) Evolution for negative λ .
Universe 09 00166 g002
Table 1. Table of critical points with fixed values of deceleration parameter q and effective equation of state w eff .
Table 1. Table of critical points with fixed values of deceleration parameter q and effective equation of state w eff .
Point X 2 Zq w eff Requirement
P m 0 Z 1 1 n ( Z ) = 2
P n 0 Z 1 1 m ( Z )
B X 2 1 1 2 ( 1 + 3 w 2 ) w 2 X 2 evaluated at Z = 1
C X 2 0 1 2 ( 1 + 3 w 2 ) w 2 X 2 evaluated at Z = 0
Table 2. Table of fixed points for the Anagnostopoulos et al. model.
Table 2. Table of fixed points for the Anagnostopoulos et al. model.
Point X 1 X 2 Zq w eff Existence ConditionsStability
A101 1 2 ( 1 + 3 w 1 ) w 1 noneSaddle
B011 1 2 ( 1 + 3 w 2 ) w 2 noneUnstable
P m 000 1 1 λ < 0 Nonhyperbolic
P n 00 2 λ / ( 1 + 2 λ ) 1 1 λ > 0 Stable
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

Böhmer, C.; Jensko, E.; Lazkoz, R. Dynamical Systems Analysis of f(Q) Gravity. Universe 2023, 9, 166. https://doi.org/10.3390/universe9040166

AMA Style

Böhmer C, Jensko E, Lazkoz R. Dynamical Systems Analysis of f(Q) Gravity. Universe. 2023; 9(4):166. https://doi.org/10.3390/universe9040166

Chicago/Turabian Style

Böhmer, Christian, Erik Jensko, and Ruth Lazkoz. 2023. "Dynamical Systems Analysis of f(Q) Gravity" Universe 9, no. 4: 166. https://doi.org/10.3390/universe9040166

APA Style

Böhmer, C., Jensko, E., & Lazkoz, R. (2023). Dynamical Systems Analysis of f(Q) Gravity. Universe, 9(4), 166. https://doi.org/10.3390/universe9040166

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