This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Homogeneous Cosmological Models in Weyl’s Geometrical Scalar–Tensor Theory

A. Barros Centro de Desenvolvimento Sustentável do Semiárido, Universidade Federal de Campina Grande, 58540-000, Sumé, PB, Brazil. [email protected]    C. Romero Departamento de Física, Universidade Federal da Paraíba. C.P. 5008, 58059-970, João Pessoa, PB, Brazil. [email protected]
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.

Weyl geometry; Scalar-tensor theory; Cosmological models.
preprint:

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 gμνg_{\mu\nu} 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 gμνg_{\mu\nu} and a scalar field Φ\Phi. 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 ω\omega, 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 σμ\sigma_{\mu} in the geometry, which bears a great similarity with the electromagnetic four-potential. Weyl went on and introduced the second-order tensor Fμν=μσννσμF_{\mu\nu}=\partial_{\mu}\sigma_{\nu}-\partial_{\nu}\sigma_{\mu}, 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 αgμν=σαgμν\nabla_{\alpha}g_{\mu\nu}=\sigma_{\alpha}g_{\mu\nu}, and is easily verified that this condition is invariant under the conformal transformation gμνg¯μν=efgμνg_{\mu\nu}\rightarrow\bar{g}_{\mu\nu}=e^{f}g_{\mu\nu} carried out simultaneously with the gauge transformation σμσ¯μ=σμ+μf\sigma_{\mu}\rightarrow\bar{\sigma}_{\mu}=\sigma_{\mu}+\partial_{\mu}f, where ff 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 Fμν=0F_{\mu\nu}=0 (null second curvature), which is equivalent to say that the one-form σ\sigma is closed (dσ=0d\sigma=0), 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 ϕ\phi, such that σμ=μϕ\sigma_{\mu}=\partial_{\mu}\phi, and, instead of a vector field σ\sigma, we are left with a scalar field ϕ\phi, 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 ϕ\phi and is explicitly given by rom14

αgμν=gμνϕ,α.\triangledown_{\alpha}g_{\mu\nu}=g_{\mu\nu}\phi_{,\alpha}. (1)

Moreover, one can define the Weyl connection, whose coefficients in a local coordinate basis read

Γμνα={μνα}12gαβ(gβμϕ,ν+gβνϕ,μgμνϕ,β),\Gamma_{\mu\nu}^{\alpha}=\{_{\mu\nu}^{\alpha}\}-\frac{1}{2}g^{\alpha\beta}(g_{\beta\mu}\phi_{,\nu}+g_{\beta\nu}\phi_{,\mu}-g_{\mu\nu}\phi_{,\beta}), (2)

with {μνα}\{_{\mu\nu}^{\alpha}\} representing the usual Christoffel symbols.

In turn, the field equations of the Weyl geometrical scalar–tensor theory can be written as barros

Gμν\displaystyle G_{\mu\nu} =(ω32)Φ2(Φ,μΦ,νgμν2Φ,αΦ,α)\displaystyle=-\frac{(\omega-\frac{3}{2})\ }{\Phi^{2}\ }\left(\Phi_{,\mu}\Phi_{,\nu}-\frac{g_{\mu\nu}}{2}\Phi_{,\alpha}\Phi^{,\alpha}\right)
1Φ(Φ,μ;νgμνΦ)gμν2ΦV(Φ)8πTμν,\displaystyle-\frac{1}{\Phi}(\Phi_{,\mu;\nu}-g_{\mu\nu}\square\Phi\ )-\frac{g_{\mu\nu}}{2\Phi}V(\Phi)-8\pi T_{\mu\nu}, (3)
Φ=1ω(12dVdΦΦ+V(Φ)),\square\Phi\ =\frac{1}{\omega}\left(-\frac{1}{2}\frac{dV}{d\Phi}\Phi+V(\Phi)\right), (4)

where here, we are using the field variable Φ=eϕ\Phi=e^{-\phi}, ω=const\omega=const, V(ϕ)V(\phi) corresponds to the scalar field potential, and TμνT_{\mu\nu} represents the Weyl invariant energy–momentum tensor of the matter fields. We denote by GμνG_{\mu\nu} and \square the Einstein tensor and the d’Alembertian operator, respectively, defined with respect to the Christoffel symbols. If V(Φ)=2ΛΦV(\Phi)=2\Lambda\Phi, one can introduce the cosmological constant Λ\Lambda. However, let us take Λ=0\Lambda=0, and then the field equations are given by

Gμν\displaystyle G_{\mu\nu} =WΦ2(Φ,μΦ,ν12gμνΦ,αΦ,α)\displaystyle=-\frac{W\ }{\Phi^{2}\ }\left(\Phi_{,\mu}\Phi_{,\nu}-\frac{1}{2}g_{\mu\nu}\Phi_{,\alpha}\Phi^{,\alpha}\right)
1ΦΦ,μ;ν8πTμν,\displaystyle-\frac{1}{\Phi}\Phi_{,\mu;\nu}-8\pi T_{\mu\nu}, (5)
Φ=0,\square\Phi\ =0, (6)

where W=ω32W=\omega-\frac{3}{2}. Additionally, we can obtain from (5) and (6) that

Rμν=8πTμν+8πT2gμνWΦ2Φ,μΦ,νΦ,μ;νΦ,R_{\mu\nu}=-8\pi T_{\mu\nu}+\frac{8\pi T}{2}g_{\mu\nu}-\frac{W}{\Phi^{2}}\Phi_{,\mu}\Phi_{,\nu}-\frac{\Phi_{,\mu;\nu}}{\Phi}, (7)

with RμνR_{\mu\nu} denoting the Ricci tensor and T=gμνTμνT=g_{\mu\nu}T^{\mu\nu}. 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

ds2=dt2+R12dx2+R22dy2+R32dz2,ds^{2}=dt^{2}+R_{1}^{2}dx^{2}+R_{2}^{2}dy^{2}+R_{3}^{2}dz^{2}, (8)

with

Ri=ri(at+b)pi1+C, R_{i}=r_{i}(at+b)^{\frac{p_{i}}{1+C}},\text{ } (9)

(i=1,2,3i=1,2,3) and the Brans–Dicke scalar field

φ=φ0(at+b)C1+C,\varphi=\varphi_{0}(at+b)^{\frac{C}{1+C}}, (10)

where aa, bb, rir_{i}, and φ0\varphi_{0} are constants. The relations pi=1{\textstyle\sum}p_{i}=1 and

pi2=1C(ωC2){\textstyle\sum}p_{i}^{2}=1-C(\omega C-2) (11)

between the constants pi,Cp_{i},C and the scalar field coupling constant ω\omega 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 a=1a=1 and b=0b=0, Equations (9) and (10) may be written as

Ri=ritpi1+C, R_{i}=r_{i}t^{\frac{p_{i}}{1+C}},\text{ } (12)
φ=φ0tC1+C. \varphi=\varphi_{0}t^{\frac{C}{1+C}}.\text{ } (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 ωW=ω3/2\omega\rightarrow W=\omega-3/2 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

Φ=Φ0tC1+C, \Phi=\Phi_{0}t^{\frac{C}{1+C}},\text{ } (14)

where pi=1{\textstyle\sum}p_{i}=1 and

pi2=1C(WC2)=1C[(ω32)C2].{\textstyle\sum}p_{i}^{2}=1-C(WC-2)=1-C\left[\left(\omega-\frac{3}{2}\right)C-2\right]. (15)

Now, if we choose C=2WC=\dfrac{2}{W}, it follows that

pi2=1.{\textstyle\sum}p_{i}^{2}=1. (16)

Furthermore, (12) and (14) become

Ri=ritWpiW+2=rit[(ω3/2)/(ω+1/2)]pi,R_{i}=r_{i}t^{\frac{Wp_{i}}{W+2}}=r_{i}t^{\left[\left(\omega-3/2\right)/\left(\omega+1/2\right)\right]p_{i}}, (17)
Φ=Φ0t2W+2=Φ0t[2/(ω+1/2)].\Phi=\Phi_{0}t^{\frac{2}{W+2}}=\Phi_{0}t^{\left[2/\left(\omega+1/2\right)\right]}. (18)

In the limit ω\omega\rightarrow\infty, (17) and (18) tend to

Ri=tpi, R_{i}=t^{p_{i}},\text{ } (19)
Φ=Φ0, \Phi=\Phi_{0},\text{ } (20)

where we have taken ri=1r_{i}=1. On the other hand, from (1) and (2) we find that

αgμν=gμν(Φ,αΦ),\nabla_{\alpha}g_{\mu\nu}=-g_{\mu\nu}\left(\frac{\Phi_{,\alpha}}{\Phi}\right), (21)
Γμνα={αμν}+12Φgαβ(gβμΦ,ν+gβνΦ,μgμνΦ,β),\Gamma_{\mu\nu}^{\alpha}=\genfrac{\{}{\}}{0.0pt}{1}{\alpha}{\mu\nu}+\frac{1}{2\Phi}g^{\alpha\beta}\left(g_{\beta\mu}\Phi_{,\nu}+g_{\beta\nu}\Phi_{,\mu}-g_{\mu\nu}\Phi_{,\beta}\right), (22)

by considering the scalar field in the form Φ=eϕ\Phi=e^{-\phi}. Thus, when ω\omega\rightarrow\infty, the space-time geometry becomes Riemannian as we have

αgμν=0, and Γμνα={αμν}.\nabla_{\alpha}g_{\mu\nu}=0,\text{ \ and\ \ }\Gamma_{\mu\nu}^{\alpha}=\genfrac{\{}{\}}{0.0pt}{1}{\alpha}{\mu\nu}. (23)

Therefore, also taking into account (19) and (20), the Kasner solution of general relativity is recovered in this limit.

V A Perfect Fluid Cosmological Model

The Friedmann–Robertson–Walker metric with a flat spatial section is given by

ds2=dt2R2(t)[dr2+r2(dϑ2+sin2ϑdχ2)],ds^{2}=dt^{2}-R^{2}(t)\left[dr^{2}+r^{2}\left(d\vartheta^{2}+\sin^{2}\vartheta d\chi^{2}\right)\right], (24)

where R(t)R(t) denotes the scale factor. In this cosmological model, the matter content is a perfect fluid represented by the energy–momentum tensor

Tμν=(p+ρ)uμuνpgμν,T_{\mu\nu}=\left(p+\rho\right)u_{\mu}u_{\nu}-pg_{\mu\nu}, (25)

with p=λρp=\lambda\rho, 0λ10\leq\lambda\leq 1, pp being the thermodynamic pressure, ρ\rho the energy density, and uμ=(1,0,0,0)u_{\mu}=\left(1,0,0,0\right) the four-velocity vector field. Then, field Equations (6) and (7) reduce to

3R¨R=4πρ(1+3λ)WΦ˙2Φ2Φ¨Φ,\frac{3\ddot{R}}{R}=-4\pi\rho\left(1+3\lambda\right)-W\frac{\dot{\Phi}^{2}}{\Phi^{2}}-\frac{\ddot{\Phi}}{\Phi}, (26)
R¨R+2R˙2R=4πρ(1λ)R˙Φ˙RΦ,\frac{\ddot{R}}{R}+\frac{2\dot{R}^{2}}{R}=4\pi\rho\left(1-\lambda\right)-\frac{\dot{R}\dot{\Phi}}{R\Phi}, (27)
Φ¨Φ+3R˙Φ˙RΦ=0.\frac{\ddot{\Phi}}{\Phi}+\frac{3\dot{R}\dot{\Phi}}{R\Phi}=0. (28)

The dot means differentiation with respect to time. Moreover, due to the assumption of spatial homogeneity, the scalar field Φ\Phi is supposed to be a function of tt only. Additionally, with the definitions θ=3R˙R\theta=\frac{3\dot{R}}{R} and Ψ=Φ˙Φ\Psi=\frac{\dot{\Phi}}{\Phi}, one can express (26)–(28) in the form

θ˙=θ234πρ(1+3λ)(W+1)Ψ2Ψ˙,\dot{\theta}=-\frac{\theta^{2}}{3}-4\pi\rho\left(1+3\lambda\right)-\left(W+1\right)\Psi^{2}-\dot{\Psi}, (29)
θ˙=θ2+12πρ(1λ)θΨ,\dot{\theta}=-\theta^{2}+12\pi\rho\left(1-\lambda\right)-\theta\Psi, (30)
Ψ˙=Ψ2θΨ.\dot{\Psi}=-\Psi^{2}-\theta\Psi. (31)

By combining (29)–(31), we can derive the equation

θ23WΨ22+θΨ=8πρ.\frac{\theta^{2}}{3}-\frac{W\Psi^{2}}{2}+\theta\Psi=8\pi\rho. (32)

After some calculations and by using Equations (5) and (6), it is easy to show that

T ;νμν=T2Φ,μΦΦ,νΦT μν,T_{\text{ \ \ };\nu}^{\mu\nu}=\frac{T}{2}\frac{\Phi^{,\mu}}{\Phi}-\frac{\Phi_{,\nu}}{\Phi}T_{\text{ \ \ }}^{\mu\nu}, (33)

which reduces to

ρ˙=[(1+λ)θ+(1+3λ2)Ψ]ρ\dot{\rho}=-\left[\left(1+\lambda\right)\theta+\left(\frac{1+3\lambda}{2}\right)\Psi\right]\rho (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):

θ˙=(1+λ)2θ2+(13λ)2θΨ3W(1λ)4Ψ2.\dot{\theta}=-\frac{\left(1+\lambda\right)}{2}\theta^{2}+\frac{\left(1-3\lambda\right)}{2}\theta\Psi-\frac{3W\left(1-\lambda\right)}{4}\Psi^{2}. (35)

This equation, together with (31), constitutes a homogeneous autonomous planar dynamic system. It is important to note that the solutions of this system, θ(t)\theta(t) and Ψ(t)\Psi(t), 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 (λ=1\lambda=1), that

θ˙=θ2θΨ.\dot{\theta}=-\theta^{2}-\theta\Psi. (36)

Clearly, an immediate solution of the system of Equations (31) and (36) is given by Ψ=θ\Psi=-\theta, which leads to the particular solution (θ=θ0(\theta=\theta_{0} , Ψ=Ψ0)\Psi=\Psi_{0}), θ0\theta_{0} and Ψ0\Psi_{0} being constants. Hence, we have 3R˙R=θ0\frac{3\dot{R}}{R}=\theta_{0}, Φ˙Φ=Ψ0\frac{\dot{\Phi}}{\Phi}=\Psi_{0}, which then leads to

R(t)=R0exp(θ0t3),R(t)=R_{0}\exp\left(\frac{\theta_{0}t}{3}\right), (37)
Φ(t)=Φ0exp(Ψ0t),\Phi(t)=\Phi_{0}\exp\left(\Psi_{0}t\right), (38)

where R0R_{0} and Φ0\Phi_{0} 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

ρ=(3W+4)48πθ02.\rho=-\frac{\left(3W+4\right)}{48\pi}\theta_{0}^{2}. (39)

Now, by defining α=θ+Ψ0\alpha=\theta+\Psi\neq 0, let us find the general solution for the stiff matter case. To do this, one can add Equations (31) and (36) to obtain

α˙+α2=0,\dot{\alpha}+\alpha^{2}=0, (40)

whose solution is

α=1t+D,\alpha=\frac{1}{t+D}, (41)

where DD is a constant. In turn, from Equation (34) with λ=1\lambda=1, it follows that

ρ˙+2αρ=0,\dot{\rho}+2\alpha\rho=0, (42)

whose solution is

ρ=ρ0(t+D)2,\rho=\frac{\rho_{0}}{\left(t+D\right)^{2}}, (43)

with ρ0\rho_{0} constant.

To obtain the expression of θ\theta, let us consider (32) and (43) and use that Ψ=αθ=1t+Dθ\Psi=\alpha-\theta=\frac{1}{t+D}-\theta. In this way we are led to

θ=Bt+D,\theta=\frac{B}{t+D}, (44)
Ψ=1Bt+D,\Psi=\frac{1-B}{t+D}, (45)

where

B=3(W+1)±3[(3+2W)16πρ0(4+3W)]4+3W,B=\frac{3(W+1)\pm\sqrt{3\left[\left(3+2W\right)-16\pi\rho_{0}\left(4+3W\right)\right]}}{4+3W}, (46)

while the condition (3+2W)16πρ0(4+3W)0\left(3+2W\right)-16\pi\rho_{0}\left(4+3W\right)\geq 0 is required to be satisfied.

On the other hand, if we replace (44) and (45) in (32), we obtain

ρ=18π(t+D)2[B23W2(1B)2+B(1B)].\rho=\frac{1}{8\pi\left(t+D\right)^{2}}\left[\frac{B^{2}}{3}-\frac{W}{2}\left(1-B\right)^{2}+B\left(1-B\right)\right]. (47)

The solutions for the scale factor and the scalar field can be obtained by integrating the expressions θ=3R˙R\theta=\frac{3\dot{R}}{R} and Ψ=Φ˙Φ\Psi=\frac{\dot{\Phi}}{\Phi}, giving the following:

R(t)=R0(t+D)B/3,R(t)=R_{0}\left(t+D\right)^{B/3}, (48)
Φ(t)=Φ0(t+D)1B,\Phi(t)=\Phi_{0}\left(t+D\right)^{1-B}, (49)

with R0R_{0} and Φ0\Phi_{0} being constants.

It should be noted that the constant BB can also be written as

B=114+3W±3[(3+2W)16πρ0(4+3W)]4+3W.B=1-\frac{1}{4+3W}\pm\frac{\sqrt{3\left[\left(3+2W\right)-16\pi\rho_{0}\left(4+3W\right)\right]}}{4+3W}. (50)

Thus, for a large WW, we obtain

B=113W±13W(248πρ0) .B=1-\frac{1}{3W}\pm\sqrt{\frac{1}{3W}\left(2-48\pi\rho_{0}\right)}\text{ }. (51)

Let us now consider that ρ0=f(W)24π\rho_{0}=\frac{f(W)}{24\pi}, f(W)f(W) being a function which tends to one when WW is large. Therefore, BB takes the form

B=113W=1+O(1W).B=1-\frac{1}{3W}=1+O\left(\frac{1}{W}\right). (52)

Under these conditions, one can obtain from (49) that

Φ(t)=Φ0+O(1W).\Phi(t)=\Phi_{0}+O\left(\frac{1}{W}\right). (53)

When Φ\Phi behaves as in (53) for a large WW, 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 WW\rightarrow\infty barros . This fact also occurs here, since Equations (48) and (49) become equal to the Einstein solution

R(t)=R0t1/3,R(t)=R_{0}t^{1/3}, (54)
Φ=Φ0,\Phi=\Phi_{0}, (55)

for WW\rightarrow\infty (we take D=0D=0). Naturally, the geometry of the space-time becomes Riemannian, according to Equations (21)–(23).

V.2 Qualitative Analysis for λ1\lambda\neq 1

For values of the parameter λ\lambda in the interval 0λ<10\leq\lambda<1, 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 θ(t)\theta(t) and Ψ(t)\Psi(t). For this purpose, let us start by writing Equations (31) and (35) of the dynamic system as

θ˙=F(θ,Ψ)=(1+λ)2θ2+(13λ)2θΨ3W(1λ)4Ψ2,\dot{\theta}=F(\theta,\Psi)=-\frac{\left(1+\lambda\right)}{2}\theta^{2}+\frac{\left(1-3\lambda\right)}{2}\theta\Psi-\frac{3W\left(1-\lambda\right)}{4}\Psi^{2}, (56)
Ψ˙=H(θ,Ψ)=Ψ2θΨ.\dot{\Psi}=H(\theta,\Psi)=-\Psi^{2}-\theta\Psi. (57)

An equilibrium point of the system, i.e., a solution that occurs when F(θ,Ψ)=H(θ,Ψ)=0F(\theta,\Psi)=H(\theta,\Psi)=0, is the origin of the phase plane, the point MM (θ=0,Ψ=0)(\theta=0,\Psi=0). 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 θ=rcosβ\theta=r\cos\beta and Ψ=rsinβ\Psi=r\sin\beta, rr and β\beta being polar coordinates of the plane. In this way, we find

θ˙=r2[(1+λ)2cos2β+(13λ)2cosβsinβ3W(1λ)4sin2β]=r2F¯(β),\dot{\theta}=r^{2}\left[-\frac{\left(1+\lambda\right)}{2}\cos^{2}\beta+\frac{\left(1-3\lambda\right)}{2}\cos\beta\sin\beta-\frac{3W\left(1-\lambda\right)}{4}\sin^{2}\beta\right]=r^{2}\overline{F}(\beta), (58)
Ψ˙=r2[sin2βcosβsinβ]=r2H¯(β).\dot{\Psi}=r^{2}\left[-\sin^{2}\beta-\cos\beta\sin\beta\right]=r^{2}\overline{H}(\beta). (59)

Now, from the relations between the variables θ\theta, Ψ\Psi, rr, β\beta, and Equations (58) and (59), it can be shown that

β˙=r(F¯(β)sinβ+H¯(β)cosβ).\dot{\beta}=r\left(-\overline{F}(\beta)\sin\beta+\overline{H}(\beta)\cos\beta\right). (60)

Next, we obtain the invariant rays, which, by definition, consist of solutions where the ratio Ψθ=tanβ=const\frac{\Psi}{\theta}=\tan\beta=const. Thus, putting β˙=0\dot{\beta}=0 in expression (60) leads to

tanβ=H¯(β)F¯(β).\tan\beta=\frac{\overline{H}(\beta)}{\overline{F}(\beta)}. (61)

Again, with the help of Equations (58) and (59), it follows from (61) that

tanβ(W2tan2βtanβ13)=0.\tan\beta\left(\frac{W}{2}\tan^{2}\beta-\tan\beta-\frac{1}{3}\right)=0. (62)

For W<32W<-\frac{3}{2}, the roots of (62) are β1=0\beta_{1}=0 and β2=π\beta_{2}=\pi. The solutions representing these invariant rays appear in phase diagrams such as curves AMAM and MAMA^{\prime}, respectively (see Figure 1, for example). When W>32W>-\frac{3}{2}, in addition to the roots β1\beta_{1} and β2\beta_{2} already mentioned, there are four more:

β3=tan1[32(1+1+2W3)]1, β4=β3+π,\beta_{3}=\tan^{-1}\left[-\frac{3}{2}\left(1+\sqrt{1+\frac{2W}{3}}\right)\right]^{-1},\text{ \ \ \ }\beta_{4}=\beta_{3}+\pi, (63)
β5=tan1[32(11+2W3)]1, β6=β5+π,\beta_{5}=\tan^{-1}\left[-\frac{3}{2}\left(1-\sqrt{1+\frac{2W}{3}}\right)\right]^{-1},\text{ \ \ \ }\beta_{6}=\beta_{5}+\pi, (64)

which correspond to the curves BMBM, MBMB^{\prime}, CMCM, and MCMC^{\prime}, respectively (see Figure 2, for instance). These invariant rays depend on WW and, as its value increases, the following behaviour is observed: the line BBBB^{\prime} rotates anticlockwise approaching the θ\theta-axis, while the line CCCC^{\prime} moves clockwise tending to make an angle of 180o-180{{}^{o}} with the positive direction of the θ\theta-axis. It should also be noted that if W=32W=-\frac{3}{2}, the lines BBBB^{\prime} and CCCC^{\prime} coincide, making an angle of 33.69o-33.69{{}^{o}} with the θ\theta-axis.

Refer to caption
Figure 1: W<32\ \ W<-\frac{3}{2} (ω<0\omega<0).
Refer to caption
Figure 2: 32<W<43\ \ -\frac{3}{2}<W<-\frac{4}{3} (0<ω<160<\omega<\frac{1}{6}).

To continue, let us check if there are regions of the phase diagrams in which ρ<0\rho<0. In these regions, the solutions should not be admitted as physical solutions, at least classically.

We start by replacing Ψ=θtanβ\Psi=\theta\tan\beta in (32). We thus obtain

θ2(W2tan2βtanβ13)=8πρ.-\theta^{2}\left(\frac{W}{2}\tan^{2}\beta-\tan\beta-\frac{1}{3}\right)=8\pi\rho. (65)

It is easy to verify, taking into account (62), that the invariant rays lying on the lines BBBB^{\prime} and CCCC^{\prime} represent vacuum solutions. Moreover, we have no region with a negative energy density if W<32W<-\frac{3}{2}. On the other hand, when W>32W>-\frac{3}{2}, we find regions where ρ<0\rho<0 that are delimited by the invariant rays that lie on the lines BBBB^{\prime} and CCCC^{\prime}. In the next section, these regions are represented as dotted regions in the phase diagrams, which widen as the value of WW increases, tending to leave the classically allowed solutions localized in a narrow region that includes the θ\theta-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 13), which are valid for λ1\lambda\neq 1 and are separated into intervals of WW (or ω\omega)111The cases W=32,W=-\frac{3}{2}, W43W-\frac{4}{3} and W=0W=0 were not analysed because they contain multiple equilibrium points or singularities..

Refer to caption
Figure 3: W>0\ \ W>0 (ω>32\omega>\frac{3}{2}).

Initially, for W<3/2W<-3/2 (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 MM (0,0)(0,0)) 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 Φ\Phi is increasing (if Ψ>0\Psi>0), while for the others, it is decreasing, in which case Ψ<0\Psi<0. On the other hand, it is possible to have singular solutions with a constant scalar field (Ψ=0\Psi=0): they are represented by the AMAM 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 MAMA^{\prime} curves, solutions that start from Minkowski’s space-time (in the infinitely distant past, with the cosmic time tt\rightarrow-\infty), and follow a contraction regime until the final collapse.

In fact, the curves AMAM and MAMA^{\prime} also correspond to solutions of general relativity, since from (56) with Ψ=0\Psi=0, it follows that

θ˙=(1+λ)2θ2,\dot{\theta}=-\frac{\left(1+\lambda\right)}{2}\theta^{2}, (66)

whose solution is

1θ=(1+λ)2t+δ,\frac{1}{\theta}=\frac{(1+\lambda)}{2}t+\delta, (67)

where δ\delta is an arbitrary constant. Therefore, by setting δ=0\delta=0, we obtain the known scale factor

R(t)=R0t2/3(1+λ).R(t)=R_{0}t^{2/3(1+\lambda)}. (68)

In Figure 2, we consider the interval 3/2<W<4/3-3/2<W<-4/3. In this diagram, there are six invariant rays: AMAM, MAMA^{\prime}, BMBM, MBMB^{\prime}, CMCM, and MCMC^{\prime}. It is interesting to recall that the dotted regions in the diagram contain solutions with ρ<0\rho<0, so that the curves restricted to these regions do not correspond to physical models. Furthermore, solutions lying on the lines BBBB^{\prime} and CCCC^{\prime} are vacuum solutions (ρ=0\rho=0), possessing singularities in their geometries, i.e., they are “big bang” models (BMBM and CMCM) or models that collapse (MBMB^{\prime}and MCMC^{\prime}), but with the scalar field varying. In the region where ρ>0\rho>0, one finds solutions similar to those in the previous diagram and also expanding universes with decreasing Φ\Phi (BMBM) and collapsing universes with increasing Φ\Phi (MBMB^{\prime}).

For W>43W>-\frac{4}{3}, it turns out that there are no nonsingular solutions in the diagrams. In Figure 3 (W>0W>0), in addition to solutions that appeared in Figure 2 when ρ0\rho\geq 0, we now observe the existence of expanding universes with increasing Φ\Phi (CMC^{\prime}M) and collapsing universes with decreasing Φ\Phi (MCMC). As mentioned before, if WW increases, the line BBBB^{\prime} moves anticlockwise approaching the line AAAA^{\prime}, while the line CCCC^{\prime} moves clockwise, also approaching AAAA^{\prime}; as a consequence, the “forbidden” regions (sectors MBCMB^{\prime}C^{\prime} and MBCMBC), where ρ<0\rho<0, become wider. In the limit WW\rightarrow\infty, the line AAAA^{\prime} remains in the region where the energy density ρ\rho is positive, representing the solutions of general relativity given by (68). Actually, for each value of WW, the line AAAA^{\prime} contains the solutions (68) because Φ=const\Phi=const (which implies Ψ=0\Psi=0) 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 A,AA,A^{\prime} B,BB,B^{\prime} C,CC,C^{\prime}
W<3/2W<-3/2 saddle points - -
3/2<W<4/3-3/2<W<-4/3 saddle points two-tangent nodes saddle points
W>4/3W>-4/3 (W0W\neq 0) 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 ω\omega\rightarrow\infty, 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 WW\rightarrow\infty. For values of the parameter λ1\lambda\neq 1, no analytical solution was possible, and we used dynamical systems theory to display the phase diagrams of the solutions in intervals of WW (or ω\omega). When W>0W>0, we highlighted solutions representing universes with ρ>0\rho>0 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 CMC^{\prime}M).

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 λ\lambda were considered. In that sense, it can be seen that Equation (62), which determines the invariant rays, did not depend on λ\lambda. Moreover, it should be noted that in the present context, matter was not a source of the geometric scalar field Φ\Phi in Equation (6). By contrast, in the Brans–Dicke theory, the scalar field equation is

φ=8πT2ω+3,\square\varphi\ =\frac{8\pi T}{2\omega+3}, (69)

where, as is well known, TT denotes the trace of the energy–momentum tensor. For the case of a perfect fluid source, T=T(λ)T=T(\lambda) and T=0T=0 only when λ=13\lambda=\frac{1}{3}. As a consequence, in the present scalar–tensor theory, cosmological models differed according to the value of the parameter λ\lambda 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 ω\omega\rightarrow\infty 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 Λ\LambdaCDM. 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.