Homogeneous Cosmological Models in Weyl’s Geometrical Scalar–Tensor Theory
Abstract
In this paper, we consider homogeneous cosmological solutions in the context of the Weyl geometrical scalar–tensor theory. Firstly, we exhibit an anisotropic Kasner type solution taking advantage of some similarities between this theory and the Brans–Dicke theory. Next, we consider an isotropic model with a flat spatial section sourced by matter configurations described by a perfect fluid. In this model, we obtain an analytical solution for the stiff matter case. For other cases, we carry out a complete qualitative analysis theory to investigate the general behaviour of the solutions, presenting some possible scenarios. In this work, we do not consider the presence of the cosmological constant nor do we take any potential of the scalar field into account. Because of this, we do not find any solution describing the acceleration of the universe.
I Introduction
As is well known, scalar–tensor theories of gravity were proposed some years ago by Jordan jor59 , and Brans and Dicke Bra61 . Later, they were extended in a more general framework ber68 ; wag70 ; nor70 . In fact, they represent a generalization of the simplest scalar–tensor theory of gravity which is the Brans–Dicke theory far04 ; wil14 . In general scalar–tensor theories of gravity, the gravitational field is not described only by the usual tensor field of general relativity. In addition to this, we have one or several long-range scalar fields which also mediate gravitational interaction.
Scalar–tensor theories of gravity have been a subject of renewed interest. Certainly, one motivation for this is the belief that, at least at sufficiently high energy scales, gravity becomes scalar–tensorial in nature dam94 and, therefore, these theories are important in the very early Universe. On the other hand, two important theoretical developments have been achieved such as, for example, unification models based on superstrings, which naturally associate long-range scalar partners to the usual tensor gravity of Einstein theory gre87 . Another motivation for the investigation of scalar–tensor theories is that inflationary cosmology in this framework seems to solve the fine-tuning problem and, in this way, give us a mechanism of terminating inflationary eras la89 . Apart from the solution of this problem, the scalar–tensor theories by themselves have direct implications for cosmology and for experimental tests of the gravitational interaction dam00 and for this reason, they are relevant in the investigation of the early Universe.
Among alternative theories of gravity, scalar–tensor theories are perhaps the most popular ones. As we have pointed out before, in these theories, gravitational effects are described by both a metric field and a scalar field . A well-known example is the Brans–Dicke theory Bra61 ; tah21 , in which the geometry of the underlying space-time manifold is assumed to be Riemannian, and the scalar field replaces the gravitational constant being interpreted as the inverse of a varying gravitational coupling parameter. In addition to the reasons mentioned above, the scalar–tensor theories are studied because they admit key ingredients of string theories, such as a dilaton-like gravitational scalar field that has a non-minimal coupling to the curvature mor11 . On the other hand, a different approach, in which the scalar field appears as part of the space-time geometry, namely, the Weyl geometrical scalar–tensor theory, has been discussed recently in the literature rom14 . Indeed, in this new approach, one considers the space-time structure as a very special case of the framework adopted in the original Weyl unified field theory wey18 ; wey52 , the geometrical space-time structure being that of a Weyl integrable space-time (WIST) rom08 ; rom09 ; rom10 ; rom11 ; rom12 . It is important to remark that other gravity theories in which a scalar field plays a geometrical role have also been proposed Fonseca ; Fonseca1 ; Fonseca2 .
Recently, some theoretical aspects concerning the Weyl geometrical scalar–tensor theory have been studied, in particular the behaviour of the solutions when , the scalar field’s coupling constant, goes to infinity barros . The investigation of cosmological vacuum models for different scalar potentials has also been carried out rom16 . In the present article, we extend this research to include anisotropic models of Kasner type. Here, we take advantage of some similarities between vacuum solutions of the Weyl geometrical scalar–tensor theory and those coming from the Brans–Dicke theory. We also examine cosmological solutions in the presence of matter, a scenario that has not yet been investigated in Weyl geometrical scalar–tensor theory, and at the same time, we compare the results obtained with similar solutions already known from general relativity and the Brans–Dicke theory.
The paper is organized as follows. In Section II, we briefly review Weyl’s original theory, which inspired the geometrical scalar–tensor approach. In Section III, the field equations of the Weyl geometrical scalar–tensor theory are obtained. Then, a Kasner type solution is exhibited in Section IV, while in Section V we work with a homogeneous and isotropic cosmological model having a perfect fluid as a source, such that we find an analytical solution for the stiff matter case and we study the other cases using the qualitative analysis of dynamical systems. Finally, Section VI is devoted to our conclusions.
II Weyl’s Theory
In the first scalar–tensor theories, the so-called Jordan–Brans–Dicke theories, it is assumed, as in general relativity, that the space-time geometry is purely Riemannian. On the other hand, if the Palatini variational method is applied to deduce the field equations from the action, then in a large class of scalar–tensor theories, a non-Riemannian compatibility condition between the metric and the affine connection appears naturally (for a more general result, see bur98 ). In this way, we have a theory that establishes the space-time geometry from first principles, that is, the space-time manifold is dynamically generated by the choice of the particular coupling of the scalar field in the gravitational sector. In the case where the action is that of the Brans–Dicke theory, this procedure leads to the so-called Weyl integrable space-times, a particular version of the geometry conceived by H. Weyl in his attempt to unify gravity and electromagnetism wey18 . Note, however, that here, it is the scalar field that is being geometrized.
It is true that the Weyl geometry is one of the simplest generalizations of Riemannian geometry, in which the Riemannian compatibility condition between the metric and the affine connection is weakened. This was an ingenious way that Weyl devised to introduce a covariant vector field in the geometry, which bears a great similarity with the electromagnetic four-potential. Weyl went on and introduced the second-order tensor , which he interpreted as representing another kind of curvature, namely, the length curvature. As a consequence of this modification in the Riemannian compatibility condition, the covariant derivative of the metric tensor does not vanish, as in Riemannian geometry, and the length of vectors when parallel transported along a curve may change. However, such theory suffered from a severe criticism by Einstein, who objected that the nonintegrability of length implies that the rate at which a clock measures time, i.e., its clock rate, in this case would depend on the past history of the clock. As a consequence of this fact, spectral lines with sharp frequencies would not appear and the spectrum of neighbouring elements of the same kind would be different goe04 . This became known in the literature as the second clock effect (incidentally, the first clock effect refers to the well-known effect corresponding to the “twin paradox”, which is predicted by special and general relativity theories).
Weyl’s new compatibility condition is given by , and is easily verified that this condition is invariant under the conformal transformation carried out simultaneously with the gauge transformation , where is an arbitrary scalar function. The discovery of this new symmetry is now considered by some authors as the birth of modern gauge theories ale13 . Now, if (null second curvature), which is equivalent to say that the one-form is closed (), then there is no electromagnetic field. In this case, we know that, from Poincaré’s lemma tu11 , it follows that there exists a scalar field , such that , and, instead of a vector field , we are left with a scalar field , which, in addition to the metric, is the fundamental object that characterizes the geometry. A space-time endowed with this particular version of Weyl’s geometry came to be known as aWeyl integrable space-time.
III The Field Equations
As we have already mentioned, in the Weyl geometrical scalar–tensor theory, the underlying space-time manifold is that of a Weyl integrable space-time rom08 . In this sense, the Weyl nonmetricity condition involves a purely geometrical scalar field and is explicitly given by rom14
(1) |
Moreover, one can define the Weyl connection, whose coefficients in a local coordinate basis read
(2) |
with representing the usual Christoffel symbols.
In turn, the field equations of the Weyl geometrical scalar–tensor theory can be written as barros
(3) |
(4) |
where here, we are using the field variable , , corresponds to the scalar field potential, and represents the Weyl invariant energy–momentum tensor of the matter fields. We denote by and the Einstein tensor and the d’Alembertian operator, respectively, defined with respect to the Christoffel symbols. If , one can introduce the cosmological constant . However, let us take , and then the field equations are given by
(5) |
(6) |
where . Additionally, we can obtain from (5) and (6) that
(7) |
with denoting the Ricci tensor and . Equations (6) and (7) constitute the field equations we use in the following.
IV Kasner Type Solution
As is well known, the Kasner metric was obtained by the mathematician E. Kasner in 1921 and represents an exact solution to Einstein’s field equations. It describes an anisotropic universe without matter, that is, it is a vacuum solution. Historically, interest in the Kasner solution came from the fact that, although it may have a singularity (“big bang” or a “big crunch”), an isotropic expansion or contraction of space is not allowed, and this led to the generic singularity studies, the so-called BKL singularities tip79 .
The Kasner type solution in the Brans–Dicke theory of gravity is given by Bra61 ; rub72
(8) |
with
(9) |
() and the Brans–Dicke scalar field
(10) |
where , , , and are constants. The relations and
(11) |
between the constants and the scalar field coupling constant are also satisfied.
The space-time given by (8) corresponds to a homogeneous universe, without matter and rotation, with distinct expansions along the three orthogonal axes, which reflects anisotropy. Note that if and , Equations (9) and (10) may be written as
(12) |
(13) |
In order to obtain a solution in the Weyl geometrical scalar–tensor theory, let us consider the following result: a vacuum solution of the Weyl geometrical scalar–tensor theory can be found if we make the change in the correspondent vacuum solution of the Brans–Dicke theory. In fact, the two theories are not physically equivalent given that in Weyl’s geometrical scalar–tensor theory test particles follow affine Weyl geodesics (autoparallels) and not metric geodesics as in the case of the Brans–Dicke theory. Nonetheless, there is a formal equality between the vacuum field equations of the two theories rom14 .
Thus, the Kasner type solution in the Weyl geometrical scalar–tensor theory is given by Equation (12) and
(14) |
where and
(15) |
Now, if we choose , it follows that
(16) |
V A Perfect Fluid Cosmological Model
The Friedmann–Robertson–Walker metric with a flat spatial section is given by
(24) |
where denotes the scale factor. In this cosmological model, the matter content is a perfect fluid represented by the energy–momentum tensor
(25) |
with , , being the thermodynamic pressure, the energy density, and the four-velocity vector field. Then, field Equations (6) and (7) reduce to
(26) |
(27) |
(28) |
The dot means differentiation with respect to time. Moreover, due to the assumption of spatial homogeneity, the scalar field is supposed to be a function of only. Additionally, with the definitions and , one can express (26)–(28) in the form
(29) |
(30) |
(31) |
After some calculations and by using Equations (5) and (6), it is easy to show that
(33) |
which reduces to
(34) |
in the context of the cosmological model considered.
V.1 Stiff Matter Solution
Next, we obtain the equations of a dynamic system which lead us to carry out a rich analysis of the solutions. For this purpose, let us consider the following equation, which results from (29)–(31):
(35) |
This equation, together with (31), constitutes a homogeneous autonomous planar dynamic system. It is important to note that the solutions of this system, and , must necessarily satisfy the constraint imposed by Equation (32).
Cosmological scenarios modelled by stiff matter have been investigated recently, particularly in connection with the problem of dark matter pie15 . Now, let us consider the stiff matter case in the geometrical scalar–tensor theory. Then, it follows from (35), in the case known as stiff matter (), that
(36) |
Clearly, an immediate solution of the system of Equations (31) and (36) is given by , which leads to the particular solution , , and being constants. Hence, we have , , which then leads to
(37) |
(38) |
where and are constants, which we recognize as a de Sitter type solution, with the scalar field also having an exponential behaviour. Furthermore, from (32), we can find
(39) |
Now, by defining , let us find the general solution for the stiff matter case. To do this, one can add Equations (31) and (36) to obtain
(40) |
whose solution is
(41) |
where is a constant. In turn, from Equation (34) with , it follows that
(42) |
whose solution is
(43) |
with constant.
To obtain the expression of , let us consider (32) and (43) and use that . In this way we are led to
(44) |
(45) |
where
(46) |
while the condition is required to be satisfied.
The solutions for the scale factor and the scalar field can be obtained by integrating the expressions and , giving the following:
(48) |
(49) |
with and being constants.
It should be noted that the constant can also be written as
(50) |
Thus, for a large , we obtain
(51) |
Let us now consider that , being a function which tends to one when is large. Therefore, takes the form
(52) |
Under these conditions, one can obtain from (49) that
(53) |
When behaves as in (53) for a large , it has been verified that any vacuum solution of the Weyl geometrical scalar–tensor theory reduces to the corresponding general relativistic solution in the limit barros . This fact also occurs here, since Equations (48) and (49) become equal to the Einstein solution
(54) |
(55) |
for (we take ). Naturally, the geometry of the space-time becomes Riemannian, according to Equations (21)–(23).
V.2 Qualitative Analysis for
For values of the parameter in the interval , we use the qualitative analysis theory andronov , by which many of the general characteristics of the integral solutions of the system can be studied without working out explicit solutions and . For this purpose, let us start by writing Equations (31) and (35) of the dynamic system as
(56) |
(57) |
An equilibrium point of the system, i.e., a solution that occurs when , is the origin of the phase plane, the point . This solution represents Minkowski’s space-time, being the only finite equilibrium point that is significant in the system.
In the qualitative analysis of solutions of Equations (56) and (57), one must construct the phase diagrams. For this, we make use of the Poincaré compactification method, which projects the phase plane into a sphere. A second mapping, in turn, projects this sphere orthogonally onto a disk, whose circumference represents the infinity of the initial phase plane andronov .
V.3 Invariant Rays and Regions of Negative Energy Density
Initially, in our analysis, we obtain the invariant rays of the dynamic system defined above. For this, let us make the change of variables and , and being polar coordinates of the plane. In this way, we find
(58) |
(59) |
Now, from the relations between the variables , , , , and Equations (58) and (59), it can be shown that
(60) |
Next, we obtain the invariant rays, which, by definition, consist of solutions where the ratio . Thus, putting in expression (60) leads to
(61) |
For , the roots of (62) are and . The solutions representing these invariant rays appear in phase diagrams such as curves and , respectively (see Figure 1, for example). When , in addition to the roots and already mentioned, there are four more:
(63) |
(64) |
which correspond to the curves , , , and , respectively (see Figure 2, for instance). These invariant rays depend on and, as its value increases, the following behaviour is observed: the line rotates anticlockwise approaching the -axis, while the line moves clockwise tending to make an angle of with the positive direction of the -axis. It should also be noted that if , the lines and coincide, making an angle of with the -axis.


To continue, let us check if there are regions of the phase diagrams in which . In these regions, the solutions should not be admitted as physical solutions, at least classically.
We start by replacing in (32). We thus obtain
(65) |
It is easy to verify, taking into account (62), that the invariant rays lying on the lines and represent vacuum solutions. Moreover, we have no region with a negative energy density if . On the other hand, when , we find regions where that are delimited by the invariant rays that lie on the lines and . In the next section, these regions are represented as dotted regions in the phase diagrams, which widen as the value of increases, tending to leave the classically allowed solutions localized in a narrow region that includes the -axis.
V.4 Phase Diagrams
Now, one can obtain the basic representation of Weyl’s cosmological solutions on the Poincaré sphere (the phase diagrams). This allows us to make a qualitative analysis of the solutions at infinity. First, let us make some comments about the diagrams (Figures 1–3), which are valid for and are separated into intervals of (or )111The cases and were not analysed because they contain multiple equilibrium points or singularities..

Initially, for (see Figure 1), the closed curves appearing in the diagram represent nonsingular cosmological models, which start in the infinitely distant past from Minkowski’s space-time (the point ) and tend to it again in the infinitely distant future; these universes present an initial phase of contraction, and then move into an expansive phase. For some of these solutions, the scalar field is increasing (if ), while for the others, it is decreasing, in which case . On the other hand, it is possible to have singular solutions with a constant scalar field (): they are represented by the curves, which correspond to solutions that start with a “big bang”, and then undergo an expansive phase, finally tending to Minkowski’s space-time, and the curves, solutions that start from Minkowski’s space-time (in the infinitely distant past, with the cosmic time ), and follow a contraction regime until the final collapse.
In fact, the curves and also correspond to solutions of general relativity, since from (56) with , it follows that
(66) |
whose solution is
(67) |
where is an arbitrary constant. Therefore, by setting , we obtain the known scale factor
(68) |
In Figure 2, we consider the interval . In this diagram, there are six invariant rays: , , , , , and . It is interesting to recall that the dotted regions in the diagram contain solutions with , so that the curves restricted to these regions do not correspond to physical models. Furthermore, solutions lying on the lines and are vacuum solutions (), possessing singularities in their geometries, i.e., they are “big bang” models ( and ) or models that collapse (and ), but with the scalar field varying. In the region where , one finds solutions similar to those in the previous diagram and also expanding universes with decreasing () and collapsing universes with increasing ().
For , it turns out that there are no nonsingular solutions in the diagrams. In Figure 3 (), in addition to solutions that appeared in Figure 2 when , we now observe the existence of expanding universes with increasing () and collapsing universes with decreasing (). As mentioned before, if increases, the line moves anticlockwise approaching the line , while the line moves clockwise, also approaching ; as a consequence, the “forbidden” regions (sectors and ), where , become wider. In the limit , the line remains in the region where the energy density is positive, representing the solutions of general relativity given by (68). Actually, for each value of , the line contains the solutions (68) because (which implies ) is a solution to Equation (6).
In most of the diagrams, the equilibrium points do not appear as isolated points. In these cases, they correspond to multiple equilibrium points, constituting the invariant rays. In the other cases they appear on the Poincaré sphere as points at the infinity, whose nature are indicating in the table below
Intervals | |||
---|---|---|---|
saddle points | - | - | |
saddle points | two-tangent nodes | saddle points | |
() | saddle points | two-tangent nodes´ | two-tangent nodes |
Table 1: Behaviour of the equilibrium points on the Poincaré sphere.
VI Conclusions
In this paper, we sought to find cosmological solutions in the context of the Weyl geometrical scalar–tensor theory. The vacuum field equations of this theory are formally identical to those of the Brans–Dicke theory, so we were able to obtain a Kasner type solution from the corresponding solution in the Brans–Dicke theory. We also found that, in the limit , the Kasner solution of general relativity was recovered. On the other hand, we investigated the existence of solutions for homogeneous and isotropic models sourced by a perfect fluid. In this case, we found an analytic solution for stiff matter and also showed that the corresponding solution of general relativity could be obtained in the limit . For values of the parameter , no analytical solution was possible, and we used dynamical systems theory to display the phase diagrams of the solutions in intervals of (or ). When , we highlighted solutions representing universes with and an increasing geometric scalar field, which started with a “big bang” and expanded to a final phase that tended toward Minkowski’s space-time (the curves ).
An interesting fact regarding the phase diagrams examined here is that there was no difference between the cosmological models when different values of the parameter were considered. In that sense, it can be seen that Equation (62), which determines the invariant rays, did not depend on . Moreover, it should be noted that in the present context, matter was not a source of the geometric scalar field in Equation (6). By contrast, in the Brans–Dicke theory, the scalar field equation is
(69) |
where, as is well known, denotes the trace of the energy–momentum tensor. For the case of a perfect fluid source, and only when . As a consequence, in the present scalar–tensor theory, cosmological models differed according to the value of the parameter rom89 .
In this work, we did not consider the presence of the cosmological constant,
nor did we take any potential of the scalar field into account. Because of
this, we did not find any solution describing the acceleration of the
universe. Incidentally, models describing cosmological scenarios in which the
acceleration of the cosmos is driven by a scalar field, quintessence
models ste99 ; toh23 and Chaplygin gas models fab17 ; ben02 among
others sen01 , have been investigated with interest. Two lines of
research that we leave for further work are (i) an investigation of the role
the geometric scalar field could play in approaching the problem of dark
matter and (ii) considering scenarios where the cosmological constant is
present with the hope that they can give some light to the problem of
dark energy.
Acknowledgements. C. Romero thanks CNPq (Brazil) for financial
support.
References
- (1) Jordan, P. Zum gegenwärtigen Stand der Diracschen kosmologischen Hypothesen. Z. Phys. 1959, 157, 112.
- (2) Brans, C.; Dicke, R.H. Mach’s Principle and a Relativistic Theory of Gravitation. Phys. Rev. 1961, 124, 925.
- (3) Bergmann, P.G. Comments on the scalar-tensor theory. Int. J. Theor. Phys. 1968, 1, 25.
- (4) Wagoner, R.V. Scalar-tensor theory and gravitational waves. Phys. Rev. D 1970, 1, 3209.
- (5) Nordtvedt, K. Post-Newtonian metric for a general class of scalar-tensor gravitational theories and observational consequences. Astrophys. J. 1970, 161, 1059.
- (6) Faraoni, V. Cosmology in Scalar-Tensor Gravity, 1st ed.; Springer: Dordrecht, The Netherlands, 2004.
- (7) Will, C.M. The Confrontation between General Relativity and Experiment. Living Rev. Relativ. 2014, 17, 4.
- (8) Damour, T.; Polyakov, A. The String Dilaton and a Least Coupling Principle. Nucl. Phys. B 1994, 423, 532.
- (9) Green, M.B.; Schwarz, J.H.; Witten, E. Superstring Theory, 1st ed.; Cambridge University Press: Cambridge, UK, 1987.
- (10) La, D.; Steinhardt, P.J. Extended Inflationary Cosmology. Phys. Rev. Lett. 1989, 62, 376.
- (11) Damour, T. Experimental tests of relativistic gravity. Nucl. Phys. B (Proc. Suppl.) 2000, 80, 41.
- (12) Tahura, S.; Nichols, D.A.; Yagi, K. Gravitational-Wave Memory Effects in Brans-Dicke Theory: Waveforms and Effects in the Post-Newtonian Approximation. Phys. Rev. D 2021, 104, 104010.
- (13) Morris, J.R. Dilatonic Effects on a Falling Test Mass in Scalar-Tensor Theory. Gen. Rel. Grav. 2011, 43, 2821.
- (14) Almeida, T.S.; Pucheu, M.L.; Romero, C.; Formiga, J.B. From Brans-Dicke Gravity to a Geometrical Scalar-Tensor Theory. Phys. Rev. D 2014, 89, 064047.
- (15) Weyl, H. Gravitation und Elektrizität. Sitzungsber. Preuss. Akad. Wiss. Phys. Math. Kl. 1918, 1918, 465.
- (16) Weyl, H. Space, Time, Matter, 1st ed.; Dover: New York, NY, USA, 1952.
- (17) Novello, M.; Heintzmann, H. Weyl Integrable Space-Time: A Model of our Cosmos? Phys. Lett. A 1983, 98, 10.
- (18) Bronnikov, K.A; Konstantinov, M.Y.; Melnikov, V. Integral Weyl Geometry in Cosmology. Grav. Cosmol. 1995, 1, 60.
- (19) Poulis, F.P.; Salim, J.M. Weyl Geometry and Gauge-invariant Gravitation. Int. J. Mod. Phys. D 2014, 23, 1450091.
- (20) Moon, T.; Lee, J.; Oh, P. Conformal Invariance in Einstein-Cartan-Weyl Space. Mod. Phys. Lett. A 2010, 25, 3129.
- (21) Scholz, E. The Unexpected Resurgence of Weyl Geometry in late 20th-Century Physics. In Beyond Einstein; Rowe, D., Sauer, T., Walter, S., Eds.; Birkhäuser: New York, NY, USA, 2018; pp. 261–360.
- (22) Paliathanasis, A. Dynamics in Interacting Scalar-Torsion Cosmology. Universe 2021, 7, 244.
- (23) Novello, M.; Bittencourt, E.; Moschella, U.; Goulart, E.; Salim, J.M.; Toniato, J.D. Geometric Scalar Theory of Gravity. JCAP 2013, 6, 14.
- (24) Sen, D.K.; Dunn, K.A. A Scalar-Tensor Theory of Gravitation in a Modified Riemannian Manifold. J. Math. Phys. 1971, 12, 578.
- (25) Barros, A.; Romero, C. A Brief Note on the Limit in Weyl Geometrical Scalar-Tensor Theory. Int. J. Mod. Phys. D 2021, 30, 2150113.
- (26) Pucheu, M.L.; Alves Junior, F.A.P.; Barreto, A.B.; Romero, C. Cosmological Models in Weyl Geometrical Scalar-Tensor Theory. Phys. Rev. D 2016, 94, 064010.
- (27) Burton, H.; Mann, R.B. Palatini variational principle for an extended Einstein-Hilbert action. Phys. Rev. D 1998, 57, 4754.
- (28) Goenner, H.F.M. On the History of Unified Field Theories. Living Rev. Relativ. 2004, 7, 2.
- (29) Afriat, A. Weyl’s Gauge Argument. Found. Phys. 2013, 43, 699.
- (30) Tu, L.W. An Introduction to Manifolds, 2nd ed.; Springer: New York, NY, USA, 2011.
- (31) Barrow, J.D.; Tipler, F.J. Analysis of the generic singularity studies by Belinskii, Khalatnikov, and Lifschitz. Phys. Rep. 1979, 56, 371.
- (32) Ruban, V.A.; Finkelstein, A.M. Generalization of the Taub-Kazner Cosmological Metric in the Scalar-Tensor Gravitation Theory. Nuovo Cimento Lett. 1972, 5, 289.
- (33) Chavanis, P.H. Cosmology with a stiff matter era. Phys. Rev. D 2015, 92, 103004.
- (34) Andronov, A.A.; Leontovich, E.A.; Gordon, I.I.; Maier, A.G. Qualitative Theory of Second Order Dynamic Systems, 1st ed.; John Wiley & Sons: New York, NY, USA, 1973.
- (35) Romero, C.; Oliveira, H.P.; de Mello Neto, J.R.T. Qualitative Analysis of Cosmological Models in Brans-Dicke Theory. Astrophys. Spce Sci. 1989, 158, 229.
- (36) Zlatev, I.; Wang, L.; Steinhardt, P.J. Quintessence, Cosmic Coincidence, and the Cosmological Constant. Phys. Rev. Lett. 1999, 82, 896.
- (37) Thompson, R.I. Evolution of Cosmological Parameters and Fundamental Constants in a Flat Quintessence Cosmology: A Dynamical Alternative to CDM. Universe 2023, 9, 172.
- (38) Del Campo, S.; Fabris, J.C.; Herrera, R.; Zimdahl, W. Is the cosmological dark sector better modeled by a generalized Chaplygin gas or by a scalar field? Eur. Phys. J. C 2017, 77, 479.
- (39) Bento, M.C.; Bertolami, O.; Sen, A.A. Generalized Chaplygin gas, accelerated expansion, and dark-energy-matter unification. Phys. Rev. D 2002, 66, 043507.
- (40) Sen, S.; Sen, A.A. Late time acceleration in Brans-Dicke cosmology. Phys. Rev. D 2001, 63, 124006.