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

Complex Scalar fields in Scalar-Tensor and Scalar-Torsion theories

Andronikos Paliathanasis [email protected] Institute of Systems Science, Durban University of Technology, PO Box 1334, Durban 4000, South Africa Instituto de Ciencias Físicas y Matemáticas, Universidad Austral de Chile, Valdivia 5090000, Chile Mathematical Physics and Computational Statistics Research Laboratory, Department of Environment, Ionian University, Zakinthos 29100, Greece
Abstract

We investigate the cosmological dynamics in a spatially flat Friedmann–Lemaître–Robertson–Walker geometry in scalar-tensor and scalar-torsion theories where the nonminimally coupled scalar field is a complex field. We derive the cosmological field equations and we make use of dimensionless variables in order to determine the stationary points and determine their stability properties. The physical properties of the stationary points are discussed while we find that the two-different theories, scalar-tensor and scalar-torsion theories, share many common features in terms of the evolution of the physical variables in the background space.

Scalar field; Complex field; Scalar-tensor; Scalar-torsion; Dynamical analysis.
pacs:
98.80.-k, 95.35.+d, 95.36.+x

I Introduction

Scalar fields in gravitational theory is a simple mechanism to introduce new degrees of freedom which can play an important role in the description of observable cosmological phenomena ratra . The early acceleration phase of the universe is attributed to a scalar field known as inflaton in0 ; in1 ; in2 ; in3 . In addition, scalar fields have been used to describe also the late-time acceleration phase of the universe attributed to dark energy, or other matter components such as dark matter, see for instance sf1 ; sf2 ; sf3 ; sf4 ; sf5 ; sf6 ; sf7 ; sf8 ; sf9 ; sf10 and references therein.

Multi-scalar fields have been widely studied in the literature. Some well-known two-scalar field models are the quintom qq1 or the Chiral model which leads to hyperbolic inflation hy4 , while other proposed multi-scalar fields theories can be found for instance in cc1 ; mm2 ; mm3 ; mm4 ; mm5 ; mm6 and references therein. A simple mechanism to introduce a multi-scalar field theory is to consider the existence of a complex scalar field, the real and imaginary parts of which give the equivalent of a two scalar-field theory com1 ; com2 .

An inflationary model with a complex scalar field was proposed in com3 . Specifically it was found that, when inflation occurs, the imaginary component of the complex scalar field does not contribute in the cosmological fluid, that is, the phase of the complex scalar field is constant. The cosmological perturbations with a complex scalar field were investigated in com4 . Furthermore, in com5 the authors used cosmological observations to reconstruct the quintessence potential for a complex scalar field. A nonminimally coupled scalar field cosmological model has been studied in com6 , while for some recent studies of complex scalar field cosmological models we refer the reader to com7 ; com8 ; com9 ; com10 ; com11 ; com12 ; com13 ; com14 .

In this study we consider the scalar-tensor and the scalar-torsion theories with a complex scalar field sc1 . In scalar-tensor theory, the scalar field is minimally coupled to gravity. The scalar field interacts with the gravitational Action Integral of Einstein’s General Relativity, that is, the Ricci-scalar of the Levi-Civita connection, the scalar-tensor theory satisfies the Machian Principle Brans and the theory is defined in the so-called Jordan frame Jord . The Brans-Dicke Lagrangian Brans is the most common scalar-tensor theory. On the other hand, the scalar-torsion theory is the equivalent of scalar-tensor model in teleparallelism. In the latter, the fundamental geometric invariant is the torsion scalar determined by the curvatureless Weitzenböck connection Hayashi79 ; Tsamp . There are many important results in the literature on the cosmological studies of the scalar-tensor and scalar-torsion theories, for an extended discussion we refer the reader to omegaBDGR ; st1 ; st2 ; st3 ; st4 ; st5 ; st6 ; st7 ; st8 ; te2 ; te3 ; te4 and references therein.

The purpose of this study is to investigate the effects of a complex scalar field in the evolution of the cosmological dynamics for the background space for these two different gravitational theories and to compare the results. Such analysis provides us with important results in order to understand the differences between the use of the Ricci-scalar and of the torsion scalar in the background geometry. We use dimensionless variables to perform a detailed analysis of the dynamical systems which describe the evolution of the physical variables. Such an approach has been widely studied before with many interesting results about the viability of proposed gravitational theories, see for instance dyn1 ; dyn2 ; dyn3 ; dyn4 ; dyn6 ; dyn5 . The plan of the paper is as follows.

In Section II we consider a spatially flat Friedmann–Lemaître–Robertson–Walker geometry in scalar-tensor gravitational theory with complex scalar field. We derive the field equations and we write the point-like Lagrangian. Moreover, we determined the stationary points for the field equations and we investigate their stability properties. In Section III we perform a similar analysis but now in the context of scalar-tensor theory. Finally, in Section IV we summarize our results and compare the physical results provided by the two theories and draw our conclusions.

II Scalar-Tensor Cosmology

Consider the complex scalar field ψ\psi in the case of scalar-tensor theory sc1 , for which the gravitational Action Integral is defined as

STensor=d4xg[F(|ψ|)R+12gμνψ,μψ,νV(|ψ|)],S_{Tensor}=\int d^{4}x\sqrt{-g}\left[F\left(\left|\psi\right|\right)R+\frac{1}{2}g^{\mu\nu}\psi_{,\mu}\psi_{,\nu}^{\ast}-V\left(\left|\psi\right|\right)\right], (1)

where RR is the Ricci scalar related to the Levi-Civita connection for the metric tensor gμνg_{\mu\nu}; |ψ|\left|\psi\right| is the norm of the complex field ψ\psi, that is, |ψ|2=ψψ\left|\psi\right|^{2}=\psi\psi^{\ast}, F(|ψ|)F(\left|\psi\right|) is the coupling function between the gravitational and the scalar field ψ\psi and V(|ψ|)V\left(\left|\psi\right|\right) is the potential function which drives the dynamics. The Action Integral (1) admits the U(1)U\left(1\right) symmetry.

The Brans-Dicke theory with a complex field is recovered for F(|ψ|)=F0|ψ|2F(\left|\psi\right|)=F_{0}\left|\psi\right|^{2}, that is, the Action Integral (1) is sc1

SBD=d4xg[F0|ψ|2R+12gμνψ,μψ,νV(|ψ|)].S_{BD}=\int d^{4}x\sqrt{-g}\left[F_{0}\left|\psi\right|^{2}R+\frac{1}{2}g^{\mu\nu}\psi_{,\mu}\psi_{,\nu}^{\ast}-V\left(\left|\psi\right|\right)\right]. (2)

In the case of a spatially flat FLRW universe with line element,

ds2=N2(t)dt2+a2(t)(dx2+dy2+dz2),ds^{2}=-N^{2}\left(t\right)dt^{2}+a^{2}\left(t\right)\left(dx^{2}+dy^{2}+dz^{2}\right)~{}, (3)

we derive for the Ricciscalar

R=6(1NH˙+12H2),R=6\left(\frac{1}{N}\dot{H}+12H^{2}\right)~{},~{} (4)

where H=1Na˙aH=\frac{1}{N}\frac{\dot{a}}{a}, a˙=dadt\dot{a}=\frac{da}{dt}, is the Hubble function.

We substitute (4) into (1) which by integration by parts gives the point-like Lagrangian

LTensor(N,a,a˙,ψ,ψ˙)=1N(6F(|ψ|)aa˙2+6F˙(|ψ|)a2a˙+12a3(ψ˙ψ˙))a3NV(|ψ|),L_{Tensor}\left(N,a,\dot{a},\psi,\dot{\psi}\right)=\frac{1}{N}\left(6F\left(\left|\psi\right|\right)a\dot{a}^{2}+6\dot{F}\left(\left|\psi\right|\right)a^{2}\dot{a}+\frac{1}{2}a^{3}\left(\dot{\psi}\dot{\psi}^{\ast}\right)\right)-a^{3}NV\left(\left|\psi\right|\right), (5)

or in the case of Brans-Dicke

LBD(N,a,a˙,ψ,ψ˙)=1N(6F0|ψ|2aa˙2+6F0(|ψ|2)a2a˙+12a3(ψ˙ψ˙))a3NV(|ψ|).L_{BD}\left(N,a,\dot{a},\psi,\dot{\psi}\right)=\frac{1}{N}\left(6F_{0}\left|\psi\right|^{2}a\dot{a}^{2}+6F_{0}\left(\left|\psi\right|^{2}\right)^{\cdot}a^{2}\dot{a}+\frac{1}{2}a^{3}\left(\dot{\psi}\dot{\psi}^{\ast}\right)\right)-a^{3}NV\left(\left|\psi\right|\right). (6)

We focus now in the Brans-Dicke theory in which the field equations are described by the point-like Lagrangian (6). Moreover, the complex scalar field is written with the use of the polar form, ψ(t)=ϕ(t)eiθ(t)\psi\left(t\right)=\phi\left(t\right)e^{i\theta\left(t\right)}, such that the Lagrangian (6) becomes

LBD(N,a,a˙,ϕ,ϕ˙,θ,θ˙)=1N(6F0ϕ2aa˙2+12Fϕa2a˙ϕ˙+12a3(ϕ˙2+ϕ2θ˙2))a3NV(ϕ).L_{BD}\left(N,a,\dot{a},\phi,\dot{\phi},\theta,\dot{\theta}\right)=\frac{1}{N}\left(6F_{0}\phi^{2}a\dot{a}^{2}+12F\phi a^{2}\dot{a}\dot{\phi}+\frac{1}{2}a^{3}\left(\dot{\phi}^{2}+\phi^{2}\dot{\theta}^{2}\right)\right)-a^{3}NV\left(\phi\right). (7)

It is obvious that the Lagrangian function (7) describes a multi-scalar field cosmological model, where ϕ\phi is the Brans-Dicke field and θ\theta is a second-scalar field minimally coupled to gravity but coupled to the Brans-Dicke field ϕ\phi. The U(1)U\left(1\right) for the point-like Lagrangian (7) provides the invariant transformation θ=θ+ε\theta=\theta+\varepsilon, and the conservation law I0=1Na3ϕ2θ˙I_{0}=\frac{1}{N}a^{3}\phi^{2}\dot{\theta}.

Variation with respect to the dynamical variables {N,a,ϕ,θ}\left\{N,a,\phi,\theta\right\} of the Lagrangian (7) provides the cosmological field equations which are

0=6F0ϕ2H2+12FϕHϕ˙+12(ϕ˙2+ϕ2θ˙2)+V(ϕ),0=6F_{0}\phi^{2}H^{2}+12F\phi H\dot{\phi}+\frac{1}{2}\left(\dot{\phi}^{2}+\phi^{2}\dot{\theta}^{2}\right)+V\left(\phi\right)~{}, (8)
0=2F0ϕ2(2H˙+3H2)+8F0Hϕϕ˙12ϕ˙2+4F0ϕ˙212ϕ˙2θ˙2+4F0ϕϕ¨+V(ϕ),0=2F_{0}\phi^{2}\left(2\dot{H}+3H^{2}\right)+8F_{0}H\phi\dot{\phi}-\frac{1}{2}\dot{\phi}^{2}+4F_{0}\dot{\phi}^{2}-\frac{1}{2}\dot{\phi}^{2}\dot{\theta}^{2}+4F_{0}\phi\ddot{\phi}+V\left(\phi\right)~{}, (9)
0=ϕ¨+ϕ(12F0H˙θ˙2)+3ϕ˙H˙+12F0ϕH2+V,ϕ0=\ddot{\phi}+\phi\left(12F_{0}\dot{H}-\dot{\theta}^{2}\right)+3\dot{\phi}\dot{H}+12F_{0}\phi H^{2}+V_{,\phi}~{} (10)

and

0=ϕθ¨+(2ϕ˙+3Hϕ)θ˙or I0=a3ϕ2θ˙,0=\phi\ddot{\theta}+\left(2\dot{\phi}+3H\phi\right)\dot{\theta}~{}\text{or }I_{0}=a^{3}\phi^{2}\dot{\theta}\text{,} (11)

where without loss of generality we have assumed N(t)=1N\left(t\right)=1.

II.1 Cosmological dynamics

In order to reconstruct the cosmological history provided by this specific complex scalar-tensor theory we make use of dimensionless variables in the context of HH-normalization and we investigate the dynamical evolution of the field equations (8)-(11) by determining the stationary points and their stability properties.

We consider the new dimensionless dependent variables dyn1

ϕ˙=23Hϕx,V(ϕ)=6H2ϕ2y,ψ˙=23Hz,V,ϕ=λV(ϕ)ϕ.\dot{\phi}=2\sqrt{3}H\phi x~{},~{}V\left(\phi\right)=6H^{2}\phi^{2}y~{},~{}\dot{\psi}=2\sqrt{3}Hz~{},~{}V_{,\phi}=\lambda\frac{V\left(\phi\right)}{\phi}. (12)

Hence, for the new independent variable, τ=lna\tau=\ln a, the field equations (8)-(11) are written as the following dynamical system

2F0(12F01)dxdτ\displaystyle 2F_{0}\left(12F_{0}-1\right)\frac{dx}{d\tau} =23(748F0)x2+(324F0)x3\displaystyle=2\sqrt{3}\left(7-48F_{0}\right)x^{2}+\left(3-24F_{0}\right)x^{3} (13)
+23F0(F0+z2+(λ3)y)+3x(F0+(4λF01)y+(18F0)z2),\displaystyle+2\sqrt{3}F_{0}\left(F_{0}+z^{2}+\left(\lambda-3\right)y\right)+3x\left(F_{0}+\left(4\lambda F_{0}-1\right)y+\left(1-8F_{0}\right)z^{2}\right)~{},
F0(12F01)1ydydτ\displaystyle F_{0}\left(12F_{0}-1\right)\frac{1}{y}\frac{dy}{d\tau} =23F0x(4λ+12F0(λ2))+(324F0)x2\displaystyle=2\sqrt{3}F_{0}x\left(4-\lambda+12F_{0}\left(\lambda-2\right)\right)+\left(3-24F_{0}\right)x^{2} (14)
+3(F0(16F01)+(4λF01)y+(18F0)z2),\displaystyle+3\left(F_{0}\left(16F_{0}-1\right)+\left(4\lambda F_{0}-1\right)y+\left(1-8F_{0}\right)z^{2}\right)~{},
2F0(12F01)13zdzdτ\displaystyle-2F_{0}\left(12F_{0}-1\right)\frac{1}{3z}\frac{dz}{d\tau} =43F0(8F01)x+(8F01)x2\displaystyle=4\sqrt{3}F_{0}\left(8F_{0}-1\right)x+\left(8F_{0}-1\right)x^{2} (15)
+(14λF0)y+(8F01)(F0+z2),\displaystyle+\left(1-4\lambda F_{0}\right)y+\left(8F_{0}-1\right)\left(F_{0}+z^{2}\right)~{},
dλdτ=23λx(1λ+Γ(λ)),Γ(λ)=ϕV,ϕϕV,ϕ\frac{d\lambda}{d\tau}=2\sqrt{3}\lambda x\left(1-\lambda+\Gamma\left(\lambda\right)\right)~{},~{}\Gamma\left(\lambda\right)=\phi\frac{V_{,\phi\phi}}{V_{,\phi}} (16)

with algebraic constraint

x2+z2+y+F0(1+43x)=0.x^{2}+z^{2}+y+F_{0}\left(1+4\sqrt{3}x\right)=0~{}. (17)

For the scalar field potential we assume V(ϕ)=V0ϕλ0V\left(\phi\right)=V_{0}\phi^{\lambda_{0}}; such that λ=λ0\lambda=\lambda_{0} is always a constant and the dimension of the dynamical system is two, after we apply the constraint equation (17). From (17) we substitute the parameter yy and we end up with a two-dimensional system on the space of variables {x,z}\left\{x,z\right\}. Moreover, we observe that the dynamical system remains invariant under the change of variable zzz\rightarrow-z. Thus, without loss of generality, we can restrict our analysis to the region with z0z\geq 0. Furthermore, for completeness of our analysis the scalar field potential can be positive or negative such that yy\in\mathbb{R}.

The stationary points of the reduced dynamical system A=(x(A),z(A))A=\left(x\left(A\right),z\left(A\right)\right) are as follows

A1=(3F0(λ4)6F0(λ+2)3,0)A_{1}=\left(-\frac{\sqrt{3}F_{0}\left(\lambda-4\right)}{6F_{0}\left(\lambda+2\right)-3},0\right)~{}

and

A2=(x2,F0x2(43F0+x2)) .A_{2}=\left(x_{2},\sqrt{-F_{0}-x_{2}\left(4\sqrt{3}F_{0}+x_{2}\right)}\right)\text{\thinspace}.

Each stationary point describes an asymptotic solution for the field equations in which the effective fluid has the equation of state parameter

weff(x,y,z)=33F0x+(334F0)x2+3(4F02+(4F0λ1)y+(18F0)z2)3F0(12F01).w_{eff}\left(x,y,z\right)=\frac{3\sqrt{3}F_{0}x+\left(3-34F_{0}\right)x^{2}+3\left(4F_{0}^{2}+\left(4F_{0}\lambda-1\right)y+\left(1-8F_{0}\right)z^{2}\right)}{3F_{0}\left(12F_{0}-1\right)}. (18)

We proceed with the discussion of the physical properties and the stability properties for the admitted stationary points.

The stationary point A1A_{1} describes a universe dominated only by the scalar field ϕ\phi and its potential function V(ϕ)V\left(\phi\right), because y(A1)=F0(12F01)(3+F0(λ10)(2+λ))3(12F0(2+λ))2y\left(A_{1}\right)=\frac{F_{0}\left(12F_{0}-1\right)\left(3+F_{0}\left(\lambda-10\right)\left(2+\lambda\right)\right)}{3\left(1-2F_{0}\left(2+\lambda\right)\right)^{2}}. The effective equation of state parameter for the asymptotic solution is weff(A1)=3+2F0(2+λ(λ9))6F0(2+λ)3w_{eff}\left(A_{1}\right)=\frac{3+2F_{0}\left(2+\lambda\left(\lambda-9\right)\right)}{6F_{0}\left(2+\lambda\right)-3}, from which it follows that A1A_{1} describes an accelerated universe as is given in Fig. 1. Moreover, the eigenvalues of the linearized system around the stationary point are derived, e1(A1)=3+F0(λ10)(2+λ)2F0(2+λ)1e_{1}\left(A_{1}\right)=\frac{3+F_{0}\left(\lambda-10\right)\left(2+\lambda\right)}{2F_{0}\left(2+\lambda\right)-1}~{},e2(A1)=3+F0(λ10)(2+λ)2F0(2+λ)1~{}e_{2}\left(A_{1}\right)=\frac{3+F_{0}\left(\lambda-10\right)\left(2+\lambda\right)}{2F_{0}\left(2+\lambda\right)-1}.In Fig. 1 we present the region plot in the two-dimensional space of the free variables {F0,λ}\left\{F_{0},\lambda\right\}, where the eigenvalues have negative real parts, that is, point A1A_{1} is an attractor.

Refer to caption
Figure 1: Region in the space of the free variables {F0,λ}\left\{F_{0},\lambda\right\}, where weff(A1)<13w_{eff}\left(A_{1}\right)<-\frac{1}{3} (Left Fig.) and where A1A_{1} is an attractor (Right Fig.). It is easy to observe that in this specific region when A1A_{1} describes an accelerated universe is also an attractor for the dynamical system.

Points A2A_{2} describe a family of stationary points which are real and physical accepted when F0x2(43F0+x2)>0-F_{0}-x_{2}\left(4\sqrt{3}F_{0}+x_{2}\right)>0. For the asymptotic solution we calculate weff(A2)=1+83x2w_{eff}\left(A_{2}\right)=1+\frac{8}{\sqrt{3}}x_{2}, which means that the family of points A2A_{2} describe accelerated universes for x2<123x_{2}<-\frac{1}{2\sqrt{3}}. Furthermore, the eigenvalues are calculated e1(A2)=6+23(2+λ)x2e_{1}\left(A_{2}\right)=6+2\sqrt{3}\left(2+\lambda\right)x_{2} and e2(A2)=0e_{2}\left(A_{2}\right)=0. Then for 6+23(2+λ)x2<06+2\sqrt{3}\left(2+\lambda\right)x_{2}<0, where e1(A2)<0e_{1}\left(A_{2}\right)<0, the  Center Manifold Theorem (CMT) can be applied in order to investigate for a possible stable submanifold and to infer about the stability. However, such an analysis does not contribute to the physical discussion of the present theory and we select to work numerically. In Fig. 2 we present the phase space portrait in the two-dimensional place (x,z)\left(x,z\right) from which we observe that the stationary points are always saddle points for e1(A2)<0e_{1}\left(A_{2}\right)<0.

From the above analysis it is clear that various major eras of the cosmological evolution can be described by the stationary points, A1A_{1} and A2A_{2}. For instance for F0=32(29λ+λ2)1F_{0}=-\frac{3}{2}\left(2-9\lambda+\lambda^{2}\right)^{-1}, weff(A1)=0w_{eff}\left(A_{1}\right)=0 which means that A1A_{1} describes the matter dominated era, and points A2A_{2} can describe the early and late acceleration phases of the universe. That is not the unique case, since, for x2=38,x_{2}=-\frac{\sqrt{3}}{8}, weff(A2)=0,w_{eff}\left(A_{2}\right)=0, and A1A_{1} can describe the future acceleration phase of the universe.

Refer to caption
Figure 2: Phase-space portrait for the scalar-tensor theory on the two-dimensional plane (x,z)\left(x,z\right) for different values of the free parameters F0F_{0} and λ\lambda. From the figure it is clear that the surface of points described by A2A_{2} are always saddle points or source points. With solid line we present the family of points A2A_{2} for e1(A2)<0e_{1}\left(A_{2}\right)<0.

II.2 Analysis at Infinity

Because variables (x,z)\left(x,z\right) are not bounded, they can take values at all the range of the real numbers, which means that they can take values at infinity. Until now we have investigated the stationary points at the finite regime. Hence in order to search for stationary points at infinity we consider the Poincaré map

x=X1X2Z2,z=Z1X2Z2,dσ=1X2Z2dτ.x=\frac{X}{\sqrt{1-X^{2}-Z^{2}}}~{},~{}z=\frac{Z}{\sqrt{1-X^{2}-Z^{2}}}~{},~{}d\sigma=\sqrt{1-X^{2}-Z^{2}}d\tau. (19)

Therefore, at infinity, 1X2Z2=01-X^{2}-Z^{2}=0, i.e. Z=1X2Z=\sqrt{1-X^{2}}, the dynamical system is reduced to the single ordinary differential equation

dXdσ=3(λ4)12F01(X21).\frac{dX}{d\sigma}=\frac{\sqrt{3}\left(\lambda-4\right)}{12F_{0}-1}\left(X^{2}-1\right). (20)

Consequently, the stationary points at the infinity B=(X(B),Z(B))B=\left(X\left(B\right),Z\left(B\right)\right) are

B1=(1,0) and B2=(1,0).B_{1}=\left(1,0\right)\text{ and }B_{2}=\left(-1,0\right)\text{.} (21)

Hence, at infinity only the scalar field ϕ\phi contributes to the cosmological solution and the physical properties of the points are similar to those of A1A_{1}.

The eigenvalues of the two-dimensional dynamical system at the stationary points at infinity are

e1(B1)=3(λ4)12F01,Re(e2(B1))=0e_{1}\left(B_{1}\right)=\frac{\sqrt{3}\left(\lambda-4\right)}{12F_{0}-1}~{},~{}\operatorname{Re}\left(e_{2}\left(B_{1}\right)\right)=0 (22)

and

e1(B2)=3(λ4)12F01,Re(e2(B2))=0.e_{1}\left(B_{2}\right)=-\frac{\sqrt{3}\left(\lambda-4\right)}{12F_{0}-1}~{},~{}\operatorname{Re}\left(e_{2}\left(B_{2}\right)\right)=0. (23)

From the phase-space portraits of Fig. 3 we can easily observe that stationary point B1B_{1} is an attractor when e1(B1)<0e_{1}\left(B_{1}\right)<0, while, when e1(B1)>0,e_{1}\left(B_{1}\right)>0, point B2B_{2} is an attractor.

Refer to caption
Figure 3: Phase-space portrait for the scalar-tensor theory on the two-dimensional plane (X,Z)\left(X,Z\right) for Z0,Z\geq 0, for different values of the free parameters F0F_{0} and λ\lambda. From the phase-space portaits we observe that points B1B_{1} and B2B_{2} can be attractors at the infinity regime for the dynamical system.

III Scalar-Torsion Cosmology

Scalar-torsion theory can be seen as the extension of scalar-tensor theory in teleparallelism. The gravitational Action Integral is te2

STorsion=d4xe[F(|ψ|)T+12gμνψ,μψ,νV(|ψ|)],S_{Torsion}=\int d^{4}x\sqrt{-e}\left[F\left(\left|\psi\right|\right)T+\frac{1}{2}g^{\mu\nu}\psi_{,\mu}\psi_{,\nu}^{\ast}-V\left(\left|\psi\right|\right)\right], (24)

where TT is the torsion scalar of the antisymmetric Weitzenböck connection, e=ge=\sqrt{-g}, eie_{i} describes the unholonomic frame, with g(ei,ej)=𝐞i𝐞i=ηijg(e_{i},e_{j})=\mathbf{e}_{i}\cdot\mathbf{e}_{i}=\eta_{ij} or in terms of coordinates, ei=hiμ(xκ)ie_{i}=h_{i}^{\mu}\left(x^{\kappa}\right)\partial_{i}, where now the Weitzenböck connection is expressed as revtel

Γ^λ=μνhaλμhνa\hat{\Gamma}^{\lambda}{}_{\mu\nu}=h_{a}^{\lambda}\partial_{\mu}h_{\nu}^{a}~{} (25)

and T=SβμνTβμνT={S_{\beta}}^{\mu\nu}{T^{\beta}}_{\mu\nu}, with TμνβT_{\mu\nu}^{\beta} to be the torsion tensor Tμνβ=Γ^νμβΓ^μνβT_{\mu\nu}^{\beta}=\hat{\Gamma}_{\nu\mu}^{\beta}-\hat{\Gamma}_{\mu\nu}^{\beta} and

Sβμν=12(Kμνβ+δβμTθνθδβνTθμθ),{S_{\beta}}^{\mu\nu}=\frac{1}{2}({K^{\mu\nu}}_{\beta}+\delta_{\beta}^{\mu}{T^{\theta\nu}}_{\theta}-\delta_{\beta}^{\nu}{T^{\theta\mu}}_{\theta}), (26)

where Kμνβ{K^{\mu\nu}}_{\beta} is the contorsion tensor

Kμνβ=12(TμνβTνμβTβμν).{K^{\mu\nu}}_{\beta}=-\frac{1}{2}({T^{\mu\nu}}_{\beta}-{T^{\nu\mu}}_{\beta}-{T_{\beta}}^{\mu\nu}). (27)

For the spatially flat FLRW spacetime with line element (3) we consider the diagonal frame

hμi(t)=diag(N(t),a(t),a(t),a(t)),h_{~{}\mu}^{i}(t)=\mathrm{diag}(N\left(t\right),a(t),a(t),a(t)), (28)

from which we calculate revtel

T=6H2.T=-6H^{2}. (29)

By replacing in the Action Integral (24) and substituting the complex scalar field as ψ=ϕeiθ\psi=\phi e^{i\theta}, we end up with the point-like Lagrangian function

LTorsion(N,a,a˙,ϕ,ϕ˙,ψ,ψ˙)=6NF(ϕ)aa˙2+a32N(ϕ˙2+ϕθ˙2)a3NV(ϕ).L_{Torsion}\left(N,a,\dot{a},\phi,\dot{\phi},\psi,\dot{\psi}\right)=-\frac{6}{N}F\left(\phi\right)a\dot{a}^{2}+\frac{a^{3}}{2N}\left(\dot{\phi}^{2}+\phi\dot{\theta}^{2}\right)-a^{3}NV\left(\phi\right). (30)

In an analogue to the Brans-Dicke model we assume now that F(ϕ)=F0ϕ2F\left(\phi\right)=F_{0}\phi^{2}. This is not a random choice, the Dilaton field, or the Brans-Dicke field for the specific potential function V(ϕ)=V0ϕ2V\left(\phi\right)=V_{0}\phi^{2} admits a discrete symmetry known as the Gasperini-Veneziano duality transformation. Similarly, the scalar-torsion theory for F(ϕ)=F0ϕ2F\left(\phi\right)=F_{0}\phi^{2} and potential function V(ϕ)=V0ϕ2V\left(\phi\right)=V_{0}\phi^{2} admits a discrete symmetry similar to the Gasperini-Veneziano duality transformation. Thus, the F(ϕ)=F0ϕ2F\left(\phi\right)=F_{0}\phi^{2} can be seen as the teleparallel Brans-Dicke equivalent model.

Variation with respect to the dynamical variables in the Lagrangian function (30) gives the field equations

0=6F0ϕ2H2+12(ϕ˙2+ϕ2θ˙2)+V(ϕ),0=-6F_{0}\phi^{2}H^{2}+\frac{1}{2}\left(\dot{\phi}^{2}+\phi^{2}\dot{\theta}^{2}\right)+V\left(\phi\right)~{}, (31)
4F0ϕ2(2H˙+3H2)+16F0ϕHϕ˙+(ϕ˙2+ϕ2θ˙2)V(ϕ)=0,4F_{0}\phi^{2}\left(2\dot{H}+3H^{2}\right)+16F_{0}\phi H\dot{\phi}+\left(\dot{\phi}^{2}+\phi^{2}\dot{\theta}^{2}\right)-V\left(\phi\right)=0~{}, (32)
ϕ¨ϕθ˙2+3H(4F0Hϕ+ϕ˙)+V,ϕ=0\ddot{\phi}-\phi\dot{\theta}^{2}+3H\left(4F_{0}H\phi+\dot{\phi}\right)+V_{,\phi}=0~{} (33)

and

ϕθ¨+θ˙(3Hϕ+2ϕ˙)=0.\phi\ddot{\theta}+\dot{\theta}\left(3H\phi+2\dot{\phi}\right)=0~{}. (34)

As in the case of the scalar-tensor theory we proceed with the analysis for the dynamics of the latter system of differential equations.

III.1 Cosmological dynamics

We make use of the dimensionless variables (12) and we write the field equations (31)-(34) in the equivalent dynamical system

F0dxdτ=(43F0+3x)x223F0(2F0+λy2z2)3x(F0+yz2),F_{0}\frac{dx}{d\tau}=\left(4\sqrt{3}F_{0}+3x\right)x^{2}-2\sqrt{3}F_{0}\left(2F_{0}+\lambda y-2z^{2}\right)-3x\left(F_{0}+y-z^{2}\right)~{},~{} (35)
F0dydτ=y(23F0(2+λ)x+3x2+3(F0y+z2)),F_{0}\frac{dy}{d\tau}=y\left(2\sqrt{3}F_{0}\left(2+\lambda\right)x+3x^{2}+3\left(F_{0}-y+z^{2}\right)\right)~{}, (36)
F0dzdτ=32z(x2y+z2F0),F_{0}\frac{dz}{d\tau}=\frac{3}{2}z\left(x^{2}-y+z^{2}-F_{0}\right)~{}, (37)
dλdτ=23λx(1λ+Γ(λ)),Γ(λ)=ϕV,ϕϕV,ϕ\frac{d\lambda}{d\tau}=2\sqrt{3}\lambda x\left(1-\lambda+\Gamma\left(\lambda\right)\right)~{},~{}~{}\Gamma\left(\lambda\right)=\phi\frac{V_{,\phi\phi}}{V_{,\phi}} (38)

with constraint equation

x2+y+z2F0=0,x^{2}+y+z^{2}-F_{0}=0, (39)

where now the parameter for the equation of state for the effective fluid is expressed as follows

weff(x,y,z)=83x+1F0(x2y+z2).w_{eff}\left(x,y,z\right)=\frac{8}{\sqrt{3}}x+\frac{1}{F_{0}}\left(x^{2}-y+z^{2}\right)\text{.} (40)

We observe that the field equations are invariant under the change of variable zzz\rightarrow-z. Thus without loss of generality we select to work in the region z0z\geq 0. Moreover, from the constraint equation (39) the dimension of the dynamical system is reduced by one, while, for the power-law potential function V(ϕ)=V0ϕλV\left(\phi\right)=V_{0}\phi^{\lambda}, the field dimension of the dynamical system is two.

The stationary points on the two-dimensional plane {x,z}\left\{x,z\right\}, are of the form C=(x(C),z(C))C=\left(x\left(C\right),z\left(C\right)\right), they are

C1=(F03(2+λ),0),C_{1}=\left(\frac{F_{0}}{\sqrt{3}}\left(2+\lambda\right),0\right)~{},
C2=(x2,F0(x2)2).C_{2}=\left(x_{2},\sqrt{F_{0}-\left(x_{2}\right)^{2}}\right)~{}.

Point C1C_{1} describes a universe in which the field ϕ\phi and the potential function contribute to the cosmological solution. The physical parameter weffw_{eff} is calculated to be weff(C1)=13(2F0(λ24)3)w_{eff}\left(C_{1}\right)=\frac{1}{3}\left(2F_{0}\left(\lambda^{2}-4\right)-3\right) from which it easily follows that C1C_{1} describes an accelerated asymptotic cosmological solution for {|λ|<2,F0>1λ24}\left\{\left|\lambda\right|<2,F_{0}>\frac{1}{\lambda^{2}-4}\right\} or {|λ|>2,F0<1λ24}\left\{\left|\lambda\right|>2,F_{0}<\frac{1}{\lambda^{2}-4}\right\} while, for |λ|=2\left|\lambda\right|=2, the de Sitter Universe is recovered. The eigenvalues of the linearized system near to the stationary point are

e1(C1)=3+F0(λ2+2),e2(C1)=3+F0(λ2+2) .e_{1}\left(C_{1}\right)=-3+F_{0}\left(\lambda^{2}+2\right)~{},~{}e_{2}\left(C_{1}\right)=-3+F_{0}\left(\lambda^{2}+2\right)\text{~{}}. (41)

Thus C1C_{1} is an attractor when 3+F0(λ2+2)<0-3+F_{0}\left(\lambda^{2}+2\right)<0.

The family of points C2C_{2} exists whenF0>0~{}F_{0}>0 and describes asymptotic solutions with weff(C2)=1+83x2w_{eff}\left(C_{2}\right)=1+\frac{8}{\sqrt{3}}x_{2} and eigenvalues

e1(C2)=6+23(2+λ)x2,e2(C2)=0.e_{1}\left(C_{2}\right)=6+2\sqrt{3}\left(2+\lambda\right)x_{2}~{},~{}e_{2}\left(C_{2}\right)=0.

Because one of the eigenvalues is zero, in Fig. 4 we present the phase-space portraits for the dynamical system in the two-dimensional plane (x,z)\left(x,z\right) for positive values of F0F_{0} and different parameters of λ\lambda. From the diagrams we observe that when e1(C2)<0e_{1}\left(C_{2}\right)<0 the family of points C2C_{2} describes attractors, otherwise the points are sources.

Refer to caption
Figure 4: Phase-space portrait for the scalar-torsion theory on the two-dimensional plane (x,z),\left(x,z\right), for different values of the free parameters F0F_{0} and λ\lambda in order to investigate the stability properties of the family of points C2C_{2} (solid lines). From the plots we observe that, when e1(C2)<0e_{1}\left(C_{2}\right)<0, the points ae attractors, otherwise they are source points

III.2 Analysis at Infinity

We consider now the Poincaré map (19) and we write the two-dimensional system in the equivalent form

F0dXdσ=(F0(1+F0)(X2+Z2))(3F0(2+λ)(X21)3X1X2Z2),F_{0}\frac{dX}{d\sigma}=\left(F_{0}-\left(1+F_{0}\right)\left(X^{2}+Z^{2}\right)\right)\left(\sqrt{3}F_{0}\left(2+\lambda\right)\left(X^{2}-1\right)-3X\sqrt{1-X^{2}-Z^{2}}\right)~{}, (42)
F0dZdσ=(F0(1+F0)(X2+Z2))(3F0(2+λ)X31X2Z2).F_{0}\frac{dZ}{d\sigma}=\left(F_{0}-\left(1+F_{0}\right)\left(X^{2}+Z^{2}\right)\right)\left(\sqrt{3}F_{0}\left(2+\lambda\right)X-3\sqrt{1-X^{2}-Z^{2}}\right). (43)

The stationary points at infinity are D=(X(D),Z(D))D=\left(X\left(D\right),Z\left(D\right)\right) with 1X(D)2Z(D)2=01-X\left(D\right)^{2}-Z\left(D\right)^{2}=0. We find the points

D1=(1,0) and D2=(1,0),D_{1}=\left(1,0\right)\text{ and }D_{2}=\left(-1,0\right)\text{,} (44)

from which it is clear that only the scalar field ϕ\phi contributes in the cosmological solution.

The eigenvalues of the linearized system are

e1(D1)=3(2+λ),Re(e2(D1))sign(F0)e_{1}\left(D_{1}\right)=-\sqrt{3}\left(2+\lambda\right)~{},~{}\operatorname{Re}\left(e_{2}\left(D_{1}\right)\right)\simeq-sign\left(F_{0}\right) (45)

and

e1(D2)=3(2+λ),Re(e2(D2))sign(F0).e_{1}\left(D_{2}\right)=\sqrt{3}\left(2+\lambda\right)~{},~{}\operatorname{Re}\left(e_{2}\left(D_{2}\right)\right)\simeq-sign\left(F_{0}\right). (46)

Therefore point D1D_{1} is an attractor when F0>0F_{0}>0 and 2+λ>02+\lambda>0, while D2D_{2} is attractor for F0>0F_{0}>0 and 2+λ<02+\lambda<0.

In Fig. 5 we present phase-space portraits for the dynamical system of scalar-torsion theory on the Poincaré variables.

Refer to caption
Figure 5: Phase-space portraits for the scalar-torsion theory on the two-dimensional plane (X,Z),\left(X,Z\right), for different values of the free parameters F0F_{0} and λ\lambda.

IV Conclusions

In this piece of work we considered a complex scalar field in the context of scalar-tensor and scalar-torsion theories in a spatially flat FLRW background geometry. For these two gravitational models we derived the field equations and we investigated the dynamical evolution of the physical quantities by using dimensional variables. With the use of the latter variables the field equations for both theories are reduced to a two-dimensional dynamical system of first-order ordinary differential equations. In order to construct the cosmological history as provided by the given models we determined the stationary points and we investigated their stability properties.

For the two different dynamical systems which correspond to the scalar-tensor and scalar-torsion theories respectively, we found that the stationary points for the two dynamical systems are four. They are two points at the finite regime and two points at infinity. The stationary points can describe accelerated asymptotic solutions which can describe the early and late-time acceleration phases of the universe. The stationary points and their stability properties are summarized in Table 1.

We can easily conclude that the two theories, that is, the scalar-torsion theory and the scalar-tensor theory have the same number of stationary points and similar stability properties. The physical variables of the background geometry at the stationary points have similar properties with a different functional dependence on the free variables F0F_{0} and λ\lambda. The dynamical equivalence of the two theories is an expected result. Torsion is not dynamical and there are no new degrees of freedom in scalar-torsion cosmology with respect to scalar-tensor cosmology.

From the above analysis it is clear that the two theories cannot be distinguished from the evolution of the background geometry and the analysis of the perturbations should be performed. However, that overpasses the scopus of this work and will be published elsewhere.

Finally, as far as the complex scalar field is concerned, we remark that there only at the (families) stationary points A2A_{2} and C2C_{2} the imaginary part of the complex scalar field contributes in the physical quantities of the asymptotic solutions. However, these families of points can describe important eras of the cosmological history, for instance, the early inflationary era, the radiation era or the matter dominated era. What will be of special interest is to investigate the existence of a scalar field potential which does not depend only on the norm of the complex scalar field.

Table 1: Stationary points and their stability properties for the cosmological models of our consideration.
Theory Point Finite/Infinity Stability
Scalar-Tensor
A1A_{1} Finite Attractor see Fig. 1
A2A_{2} Finite Saddle
B1B_{1} Infinity Attractor 3(λ4)12F01<0\frac{\sqrt{3}\left(\lambda-4\right)}{12F_{0}-1}<0
B2B_{2} Infinity Attractor 3(λ4)12F01<0-\frac{\sqrt{3}\left(\lambda-4\right)}{12F_{0}-1}<0
Scalar-Torsion
C1C_{1} Finite Attractor 3+F0(λ2+2)<0-3+F_{0}\left(\lambda^{2}+2\right)<0
C2C_{2} Finite Attractor 6+23(2+λ)x2<06+2\sqrt{3}\left(2+\lambda\right)x_{2}<0
D1D_{1} Infinity Attractor F0>0F_{0}>0, 2+λ>02+\lambda>0
D2D_{2} Infinity Attractor F0>0F_{0}>0, 2+λ>02+\lambda>0

Acknowledgements.
This work was partially supported by the National Research Foundation of South Africa (Grant Numbers 131604). The author thanks Dr. G. Anargirou for the hospitality provided while part of this work carried out.

References

  • (1) B. Ratra and P.J.E Peebles, Phys. Rev. D 37 3406 (1988)
  • (2) J. Rubio, Front. Astron. Space Sci. 5, 50 (2019)
  • (3) L.N. Granda and D.F. Jimenez, EPJC 79, 772 (2019)
  • (4) P. Parsons and J.D. Barrow, Phys. Rev. D 51, 6757 (1995)
  • (5) J.D. Barrow and A. Paliathanasis, Gen. Rel. Grav. 50, 82 (2018)
  • (6) E.V Linder, Phys. Rev. D. 70 023511 (2004)
  • (7) S. Barshay and G. Kreyerhoff, EPJC 5, 369 (1998)
  • (8) T. Matos, J.-R. Luevano, I. Quiros, L.A. Urena-Lopez and J.A. Vazquez, Phys. Rev. D 80, 123521 (2009)
  • (9) A.R. Liddle, C. Pahud and L.A. Urena-Lopez, Phys. Rev. D 77, 121301 (2008)
  • (10) L.P. Chimento, V. Mendez and N. Zuccala, Class. Quantum Grav. 16, 3749 (1999)
  • (11) C. Deffayet, S. Deser and G. Esposito-Farese, Phys. Rev. D 80, 064015 (2009)
  • (12) D. Bertacca, N. Bartolo and S. Matarrese, Adv. Astron. 2010, 904379 (2010)
  • (13) D. Bertacca, A. Raccanelli, O.F. Piatella, D. Pietrobon, N. Bartolo, S. Matarrese and T. Giannantonio, JCAP 03, 039 (2011)
  • (14) A. Suárez and T. Matos, MNRAS 416, 87 (2011)
  • (15) C. Cosme, J.G. Rosa and O. Bertolami, Phys. Lett. B 781, 639 (2018)
  • (16) W. Hu, Phys. Rev. D 71, 047301 (2005)
  • (17) D. Wands, Lect. Notes Phys. 738, 275 (2008)
  • (18) A.A. Coley and R.J. van den Hoogen, Phys. Rev. D 62, 023517 (2000)  
  • (19) A. Arbey, Phys. Rev. D 74, 043516 (2006)
  • (20) C.-J. Gao and Y.-G. Shen, Phys. Lett. B 541, 1 (2002)
  • (21) D.S.M. Alves and G.M. Kremer, JCAP 10, 009 (2004)
  • (22) V. Sivanesan, Phys. Rev. D 90, 104006 (2014)
  • (23) A. Paliathanasis, Universe 8, 325 (2022)
  • (24) I.M. Khalatnikov and A. Mezhlumian, Phys. Lett. A 169, 308 (1992)
  • (25) I.M. Khalatnikov, Lect. Notes Phys. 455, 343 (1995)
  • (26) D. Scialom, Helv. Phys. Acta 69, 190 (1996)
  • (27) P. Jetzer and D. Scialom, Phys. Rev. D 55, 7440 (1997)
  • (28) J.-A. Gu and W.-Y. Hwang, Phys. Lett. B 517, 1 (2001)
  • (29) A.Yu.Kamenshchik, I.M.Khalatnikov and A.V.Toporensky, Phys. Lett. B 357, 36 (1995)
  • (30) M. Arik, M. Calik and N. Katirci, Central Eur. J. Phys. 9, 1465 (2011)
  • (31) G. Rosen, EPL 89, 19002 (2010)
  • (32) Y.-G. Shen and X.-H. Ge, Int. J. Theor. Phys. 45, 17 (2006)
  • (33) A. Arbey, Phys. Rev. D 74, 043516 (2006)
  • (34) R.C.G. Landim, EPJC 76, 1 (2016)
  • (35) B. Li, P.R. Shapiro and T. Rindler-Daller, PoS BASH2015, 028 (2016)
  • (36) A. Bevilacqua, J. Kowalski-Glikman and W. Wislicki, Phys. Rev. D 105, 105004 (2022)
  • (37) H. Foidl and T. Rindler-Daller,  Phys. Rev. D 105, 123534 (2022)
  • (38) V. Faraoni, Cosmology in Scalar-Tensor Gravity, Fundamental Theories of Physics, Springer Dordrecht (2004)
  • (39) C. Brans and R.H. Dicke, Phys. Rev. 124, 195 (1961)
  • (40) P. Jordan, Nature 164, 637 (1937) and Schwerkraft und Weltfall, 2nd ed., Vieweg und Sohn, Braunschweig, (1955)
  • (41) K. Hayashi and T. Shirafuji, Phys. Rev. D 19, 3524 (1979)
  • (42) M. Tsamparlis, Phys. Lett. A 75, 27 (1979)
  • (43) V. Faraoni, Phys. Rev. D 59, 084021, (1999)
  • (44) J. Garcia-Bellido, D. Wands, Phys. Rev. D 52, 6739 (1995)
  • (45) S.Sen and A.A. Sen, Late time acceleration in Brans Dicke Cosmology, Phys. Rev. D 63, 124006 (2001)
  • (46) M. Artymowski, Z. Lalak and M. Lewicki, Inflation and dark energy from the Brans-Dicke theory, JCAP 06, 031 (2015)
  • (47) M. Gasperini and G. Veneziano, Astropart. Phys. 1, 317 (1993)
  • (48) L. Jarv and J. Lember, Universe 7, 178 (2021)
  • (49) M. Hohmann, Phys. Rev. D 98, 064002 (2018)
  • (50) M. Hohmann and C. Pfeifer, Phys. Rev. D 98, 064003 (2018)
  • (51) H. Wei, Phys. Dynamics of Teleparallel Dark Energy, Lett. B 712, 430 (2012)
  • (52) C.-Q. Geng, C.-C. Lee, E.N. Saridakis and Y.-P. Wu, ”Teleparallel” Dark Energy, Phys. Lett. B 704, 384 (2011)
  • (53) C. Xu, E.N. Saridakis and G. Leon, Phase-Space analysis of Teleparallel Dark Energy, JCAP 07, 005 (2012)
  • (54) A. Paliathanasis, Eur. Phys. J. Plus 136, 674 (2021)
  • (55) E.J. Copeland, A.R. Liddle and D. Wands, Phys. Rev. D 57, 4686 (1998)
  • (56) . Wainwright and G.F.R Ellis, Cambridge University Press, Cambridge (1997)
  • (57) A.A. Coley, Dynamical Systems and Cosmology, Springer, Dordrecht (2003)
  • (58) L. Amendola and S. Tsujikawa, Dark Energy, Cambridge University Press, Cambridge (2010)
  • (59) C.R. Fadragas and G. Leon, Class. Quantum Grav. 31, 195011 (2014)
  • (60) O. Hrycyna and M. Szydlowski, JCAP 1312, 016 (2013)
  • (61) S. Bahamonte, K.F. Dialektopoulos, C. Escamilla-Rivera, V. Gakis, M. Hendry, J.L. Said, J. Mifsud and E. Di Valentino, Teleparallel Gravity: From Theory to Cosmology, [arXiv:2106.13793] (2021)