Dynamics of Einstein-Aether Scalar field Cosmology

TL;DR

This paper investigates the evolution of Einstein-Aether cosmology with a scalar field in flat spacetime, analyzing stability, critical points, and solutions including inflationary and late-time attractors.

Contribution

It provides a comprehensive dynamical systems analysis of Einstein-Aether scalar field cosmology, including stability and critical points for various potentials.

Findings

Recovery of general relativity limit

Existence of de Sitter solutions for inflation and late-time acceleration

Numerical results for specific scalar potentials

Abstract

We study the cosmological evolution of the field equations in the context of Einstein-Aether cosmology by including a scalar field in a spatially flat Friedmann--Lema\^{\i}tre--Robertson--Walker spacetime. Our analysis is separated into two separate where a pressureless fluid source is included or absent. In particular, we determine the critical points of the field equations and we study the stability of the specific solutions. The limit of general relativity is fully recovered, while the dynamical system admits de Sitter solutions which can describe the past inflationary era and the future late-time attractor. Results for generic scalar field potentials are presented while some numerical behaviours are given for specific potential forms.

Equations

$$d{s}^2=-d{t}^2+a^2\left(t\right)\left(d{r}^2+r^{2}d{\theta }^2+r^2{\sin }^2\theta d{\varphi }^2\right),$$

$$S_{EA}=\int \sqrt{-g}d{x}^4\left(R+\frac{c_{\theta }}{3}{\theta }^2+c_{\sigma }{\sigma }^2+c_{\omega }{\omega }^2+c_{\alpha }{\alpha }^2\right)+S_{\varphi },$$

$$S_{\varphi }=\int \sqrt{-g}d{x}^4\left(\frac{1}{2}{g}^{\mu \nu }{\varphi }_{,\mu }{\varphi }_{,\nu }-V\left(\theta ,\varphi \right)\right).$$

$$\frac{1}{3}{\theta }^2=\frac{1}{2}{\dot{\varphi }}^2+V-\theta V_{\theta },$$

$$\frac{2}{3}\dot{\theta }=-{\dot{\varphi }}^2-\dot{\theta }{V}_{\theta \theta }-\dot{\varphi }{V}_{\theta \varphi },$$

$$\ddot{\varphi }+\theta \dot{\varphi }+V_{\varphi }=0.$$

$$V\left(\varphi ,\theta \right)=V_0{e}^{-\lambda \theta }+\sum _{r=0}^n{V}_r{\theta }^r{e}^{\frac{r-2}{2}\lambda \varphi },$$

$$V\left(\varphi ,\theta \right)=U\left(\varphi \right)+Y\left(\varphi \right)\theta ,$$

$$\frac{1}{3}{\theta }^2=\frac{1}{2}{\dot{\varphi }}^2+U\left(\varphi \right),$$

$$\frac{2}{3}\dot{\theta }+\frac{1}{3}{\theta }^2=-\frac{1}{2}{\dot{\varphi }}^2+U\left(\varphi \right)-\dot{\varphi }{Y}_{,\varphi },$$

$$\ddot{\varphi }+\theta \dot{\varphi }+U_{,\varphi }+Y_{,\varphi }\theta =0.$$

$$G_{ab}=T_{ab},$$

$$T_{ab}={\rho }_{\varphi }{u}_a{u}_b+p_{\varphi }{h}_{ab},$$

$${\rho }_{\varphi }=\frac{1}{2}{\dot{\varphi }}^2+U\left(\varphi \right),p_{\varphi }=\frac{1}{2}{\dot{\varphi }}^2-U\left(\varphi \right)+\dot{\varphi }{Y}_{,\varphi }.$$

$$\frac{1}{3}{\theta }^2=\frac{1}{2}{\dot{\varphi }}^2+U\left(\varphi \right)+{\rho }_m$$

$$\frac{2}{3}\dot{\theta }+\frac{1}{3}{\theta }^2=-\frac{1}{2}{\dot{\varphi }}^2+U\left(\varphi \right)-\dot{\varphi }{Y}_{,\varphi }-w_m{\rho }_m$$

$${\dot{\rho }}_m+\theta \left(1+w_m\right){\rho }_m=0.$$

$$x=\sqrt{\frac{3}{2}\frac{{\dot{\varphi }}^2}{{\theta }^2}},y=\sqrt{\frac{3U}{{\theta }^2}},\lambda =\frac{U_{,\varphi }}{U},\xi =\sqrt{2}\frac{Y_{,\varphi }}{\sqrt{U}},{\Omega }_m=\frac{3{\rho }_m}{{\theta }^2}.$$

$${\Omega }_m=1-x^2-y^2,$$

$$\frac{dx}{dN}=-3x(1+J\left(x,y,\xi ,{\Omega }_m\right))-\frac{1}{2}\left(\sqrt{6}\lambda y+3\xi \right)y,$$

$$\frac{dy}{dN}=\sqrt{\frac{3}{2}}\lambda xy-3yJ,$$

$$\frac{d\lambda}{dN}=\sqrt{6}x\lambda^{2}(\Gamma_{\lambda}\left(\lambda\right)-1),\$$

$$\frac{d\xi }{dN}=\sqrt{3}\xi x\left({\Gamma }_{\xi }\left(\xi \right)-\frac{\lambda }{\sqrt{2}}\right),$$

$$2J\left(x,y,\xi ,{\Omega }_m\right)=-1-x^2+y^2-\xi xy-w_m{\Omega }_m.$$

$${\Gamma }_{\lambda }\left(\lambda \right)=\frac{U_{,\varphi \varphi }U}{(U_{,\varphi }{)}^2},{\Gamma }_{\xi }\left(\xi \right)=\sqrt{2}\frac{Y_{,\varphi \varphi }}{Y_{,\varphi }}$$

$$U\left(\varphi \right)=U_0{e}^{\lambda \varphi },Y\left(\varphi \right)=Y_0+Y_1{e}^{-\frac{\lambda }{2}\varphi },$$

$$w_{\varphi }\left(x,y,\lambda ,\xi ,{\Omega }_m\right)=-1+\frac{2x^2+\xi xy}{x^2+y^2},$$

$$q\left(x,y,\lambda ,\xi ,{\Omega }_m\right)=-1-2J,$$

$$w_{\mathrm{tot}}\left(x,y,\lambda ,\xi ,{\Omega }_m\right)=\frac{1}{3}\left(q-1\right)\text{.}$$

$$w_{\varphi }\left(\lambda ,\xi \right)=-1+\lambda \frac{4\lambda +\xi \sqrt{6\left(4+{\xi }^2\right)-4{\lambda }^2}}{3\left(4+{\xi }^2\right)}.$$

Full text

Cosmology; Modified theories of gravity; Einstein-Aether; Scalar field; Critical points

I Introduction

Einstein-Aether theory is a Lorentz-violating theory in which a unitary timelike vector field, called the æther, is introduced into the Einstein-Hilbert action DJ ; DJ2 ; Carru ; carroll . The introduction of the timelike vector field in the action integral is also a specific selection of preferred frame at each point in the spacetime, and so this modification spontaneously breaks the Lorentz symmetry ea1 . The gravitational field equations are of second-order and correspond to variations of the action with respect to the metric tensor and the æther field. At this point we recall that the unitarity of the timelike vector field is guaranteed by introducing a lagrange multiplier. The Einstein-Aether theory can describe various cosmological phases, including those of early inflationary expansion and late dark-energy dominationdata1 ; data2 ; data3 ; data4 . It is important to mention here that the Einstein-Aether approach also describes the classical limit of Hořava gravity esf .

One of the ways to study a cosmological model is to perform a dynamical analysis by studying its critical points in order to connect them to the different observed eras, with their respective dynamical behaviours and characteristics DynSystemsWain ; DynSystemsColey ; Heinzle ; dyn1 ; dyn2 ; dyn3 ; dyn4 ; dyn5 ; dyn6 ; dyn7 ; dyn9 . For the Einstein-Aether cosmologies there have been several such studies dyn8 ; dyn8a ; dyn08b ; dyn08c ; dyn08e ; dyn08e ; coleysf ; KSEAPF .

For Einstein-Aether cosmologies Barrow provided exact solutions for specific forms of the scalar field potential in the framework of Friedmann–Lemaître–Robertson–Walker (FLRW) spacetime. There has been further study of the dynamical evolution and stability of those inflationary solutions in homogeneous and isotropic Einstein-Aether cosmologies containing a self interacting scalar field which interacts with the aether coleysf . Similar dynamical analysis can by found in Alhulaimi , where it was shown that for isotropic expansion the dynamics are independent of the aether parameters, but this is not the case for anisotropic expansion. In all cases there is a period of slow-roll inflation at intermediate times and, in some cases, accelerated expansion at late times.

Apart from the FLRW background scenario, there have been more wideranging studies. One such work, investigating the dynamical equations of the Einstein-Aether theory for the cases of FLRW as well as in an locally rotationally symmetric Bianchi Type III geometry LangRoum . There, it was found that the existence or the non-existence of the solutions to the reduced equations depends on the values of combinations of the initial parameters that enter the action integral. Results of this type have also been found elsewhere Foster ; Jacobson ; Zlosnik . For other dynamical studies in the context of the Einstein-Aether scenario we refer the reader the articles of KSEAPF ; ColeyLeon eda5 ; eda6 ; eda7 ; eda8 .

The plan of this paper is as follows. In Section II we present the model to be studied, which is an Einstein-Aether scalar field cosmology with spatially flat FLRW spacetime, where the scalar field lagrangian has been modified so that the scalar field potential is non-minimally coupled to the aether field, as proposed in DJ . The dimensionless dynamical analysis and the corresponding critical points are presented in Section III. The absence of the matter in the action integral implies that the dimension of the dynamical system can be either one or two, while adding a pressureless fluid raises the dimension of the system to two or three. Furthermore, the critical points are classified into three families. Sections IV and V include the main results of the current analysis, where we present the allowed critical points, It is interesting to mention that the case of general relativity (GR) with a minimally coupled scalar field is fully recovered, while new critical points are found which describe either power-law or de Sitter solutions. Finally, we draw our conclusions in Section VII.

II Einstein-aether cosmology

First, we consider a spatially flat FLRW spacetime with line element

[TABLE]

and as an aether field we consider the timelike vector field $u^{\mu}=\delta_{t}^{\mu}$. In Einstein-Aether Models the gravitational action is given by jacobo1

[TABLE]

where the action integral of the scalar field is assumed to be DJ

[TABLE]

Parameters $\theta,~{}\sigma,~{}\omega$ and $\alpha$ describe the kinematic quantities of the unitary vector field $u^{\mu},$ and correspond to the volume  expansion rate, shear, vorticity and fluid acceleration. Following the notations of DJ for the line element (1), the field equations are written as (with units $8\pi G=c=1$)

[TABLE]

[TABLE]

[TABLE]

Barrow Barrow previously found that for

[TABLE]

where $V_{0},~{}V_{r}$ and $\lambda$ are constants, the field equations (4)-(6) admit exact power-law solutions $\phi\left(t\right)=\frac{2}{\lambda}\ln t~{}$, $\theta\left(t\right)=3Bt^{-1}$, i.e. $a\left(t\right)\simeq t^{B}$ where $B=B\left(V_{0},V_{r},\lambda\right)$. The detailed dynamical analysis of (7) was performed in coleysf .

In this work we consider the potential function $V\left(\phi,\theta\right)$ to be

[TABLE]

where $U\left(\phi\right)$ now denotes the scalar field potential and $Y\left(\phi\right)$ is the coupling term between the scalar and aether fields. It is straightforward to observe that for $U\left(\phi\right)=U_{0}e^{-\lambda\theta}$ and $Y\left(\phi\right)\simeq\sqrt{U\left(\phi\right)}$ potential (7) is recovered for $n=1,V_{0}=0$.

With the aid of (8) the field equations (4)-(6) simplify to

[TABLE]

[TABLE]

[TABLE]

from which we can rewrite them as

[TABLE]

where $G_{ab}$ is the Einstein tensor, and $T_{ab}$ describes the effective Einstein-Aether fluid source with the scalar field, where in the $1+3$ decomposition is written as

[TABLE]

in which $h_{ab}=g_{ab}+u_{a}u_{b}$ is the projective tensor and $\rho_{\phi}~{}$and $p_{\phi}$ are given by

[TABLE]

We observe that for this specific potential, (8), the energy density $\rho_{\phi}$ is defined as in Einstein’s GR, while in the pressure term $p_{\phi}$ a new part is introduced due to the coupling between the scalar field and the aether field.

If a minimally coupled matter source is introduced with energy density $\rho_{m}$, pressure $p_{m}$, and constant parameter for the equation of state, $w_{m}=p_{m}/\rho_{m},$ the field equations (9)-(10) are modified as follows:

[TABLE]

[TABLE]

where, for the perfect fluid $\rho_{m}$, the conservation law is

[TABLE]

We carry out our analysis by writing the field equations (15)-(17) in dimensionless variables by using expansion-normalised variables.

III Dynamical system

In this section we present the main features of the dynamical analysis by using the method of critical points. This method is powerful because it provides information concerning the general evolution of the dynamical system. Hence, from such an analysis the overall cosmological viability of the model can be discussed.

The new dimensionless variables are defined as follows copeland

[TABLE]

After some calculations, the system of the field equations is written in these variables as

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

where $N=\ln\left(a\right)$, function $J\left(x,y,\xi,\Omega_{m}\right)$ is expressed as

[TABLE]

and

[TABLE]

In general, equations (19)-(23) form a three-dimensional system. Specifically, equations (22), (23) do not coexist, because parameters $\lambda$ and $\xi$ are not independent. Locally, the condition $\frac{\partial\lambda}{\partial\phi}\neq 0$ implies that the inverse function of $\lambda\left(\phi\right)$ exists so that $\phi=\phi\left(\lambda\right)$. Hence the function $\xi\left(\phi\right)$ depends on $\lambda$; indeed, $\xi=\xi\left(\phi\left(\lambda\right)\right)$, so $\xi=\xi\left(\lambda\right)$. In that case the only independent variables which survive are the $\left\{x,y,\lambda\right\}$.

On the other hand, if locally $\lambda=const$, then $\frac{d\lambda}{dN}\equiv 0$, and $\xi=\xi\left(\phi\right)$; hence, the only independent variables that survive are $\left\{x,y,\xi\right\},$ since $\phi=\phi\left(\xi\right)$. Consequently, there are three large families of potentials which we will study, that admit different dynamical systems:

Family (A) with $\lambda=const$. and $\xi=cons\not t,$ which corresponds to the potentials

[TABLE]

Family (B) where $U\left(\phi\right)=U_{0}e^{\lambda\phi}$, $\lambda=const$ with $\xi=\xi\left(\phi\right),$ and

Family (C) where the potential $U\left(\phi\right)$ is different from the exponential potential and so $\xi=\xi\left(\lambda\right).$

We observe from equation (19) that the variables $x,y$ obey the inequalities $0\leq x^{2}+y^{2}\leq 1$, where in the special case of $\Omega_{m}=0$, this reduces to $x^{2}+y^{2}=1.$ At this point we mention that the equation of state parameter for the perfect fluid is now

[TABLE]

while the deceleration parameter is

[TABLE]

and the equation of state parameter for the effective fluid is

[TABLE]

We continue our study with the analysis of the critical points of the system (19)-(23) for the aforementioned three families of potentials.

IV Scalar field without matter source

In this section we consider the case where the cosmic fluid does not include matter, namely $\Omega_{m}=0$. In this case our results are summarized as follows.

IV.1 Family A

For the first family of potentials the constraint (19) implies that the dynamical system (20)-(21) is reduced to a one-dimensional system. In this case we derive four critical points $\left(x_{P},y_{P}\right)$:

a. Points $A_{1\left(\pm\right)}$ with coordinates $A_{1\left(\pm\right)}=\left(\pm 1,0\right)$. These points describe a universe where $\rho_{\phi}=p_{\phi}=\dot{\phi}^{2}/2$, hence $w_{\phi}=\frac{p_{\phi}}{\rho_{\phi}}=1$. For the stability of the points we need to calculate the corresponding eigenvalues, which in this case are given by $e_{1}\left(A_{1\left(\pm\right)}\right)=3\pm\sqrt{\frac{3}{2}}\lambda.$ Therefore, point $A_{1\left(+\right)}$ is stable for $\lambda<-\sqrt{6}$ while point $A_{1\left(-\right)}$ is stable for $\lambda>\sqrt{6}$.

b. Point $A_{2}$ with coordinates $A_{2}=\left(-\frac{2\sqrt{2}\lambda+\xi\sqrt{3\left(4+\xi^{2}\right)-2\lambda^{2}}}{\sqrt{3}\left(4+\xi^{2}\right)},\frac{-\sqrt{2}\lambda\xi+2\sqrt{3\left(4+\xi^{2}\right)-2\lambda^{2}}}{\sqrt{3}\xi\left(4+\xi^{2}\right)}\right)~{}$exists for $\left\{\lambda<-\sqrt{6}~{},~{}\xi\geq\sqrt{\frac{2\left(\lambda^{2}-6\right)}{3}}\right\}~{}$or $\left\{\lambda>\sqrt{6}~{},~{}\xi\leq-\sqrt{\frac{2\left(\lambda^{2}-6\right)}{3}}\right\}$ or $\left\{-\sqrt{6}\leq\lambda\leq\sqrt{6}~{},~{}\xi\neq 0\right\}.$ The equation of state parameter is written as

[TABLE]

Point $A_{2}$ describes an accelerated universe when $w_{\phi}<-\frac{1}{3}$; that is, the parameters $\lambda,\xi$ are constrained by the following conditions  $\left\{\lambda<-\sqrt{2}~{},~{}\xi>\sqrt{2}\frac{\lambda^{2}-2}{\sqrt{3\lambda^{2}-2}}\right\}~{}$or  $\left\{\lambda>\frac{\sqrt{6}}{3}~{},~{}0<\xi<-\sqrt{2}\frac{\lambda^{2}-2}{\sqrt{3\lambda^{2}-2}}\right\}~{}$or $\left\{-\sqrt{2}\leq\lambda\leq\frac{\sqrt{6}}{3}~{},~{}\xi>0\right\}$ or $\left\{\lambda\geq\sqrt{6}~{},~{}\xi<-\sqrt{2}\frac{\lambda^{2}-2}{\sqrt{3\lambda^{2}-2}}\right\}.$ As far as the stability is concerned we find that $A_{2}$ is stable when $\left\{0<\lambda\leq\sqrt{6}~{},~{}\xi>0\right\}~{}$or $\left\{\lambda=\sqrt{6}~{},~{}\xi<0\right\}$ or $\left\{\lambda>\sqrt{6}~{},~{}\xi\leq-\sqrt{\frac{2}{3}\left(\lambda^{2}-6\right)}\right\}.~{}$ Moreover, if $A_{2}$ is an attractor describing cosmic acceleration, then the parameters $\lambda,\xi$ obey the inequalities $\left\{-\sqrt{2}<\lambda\leq\sqrt{\frac{2}{3}}~{},~{}\xi>0\right\}~{}$or $\left\{\lambda<-\sqrt{2}~{},~{}0<\xi<\sqrt{2}\frac{\lambda^{2}-2}{\sqrt{3\lambda-2}}\right\}$  or $\left\{\lambda\geq\frac{\sqrt{6}~{}}{3},0<~{}\xi<-\sqrt{2}\frac{\lambda^{2}-2}{\sqrt{3\lambda-2}}\right\}$ or $\left\{\lambda\geq\sqrt{6}~{},~{}\xi<-\sqrt{2}\frac{\lambda^{2}-2}{\sqrt{3\lambda-2}}~{}\right\}$

c. Point $A_{3}$  $=\left(\frac{-2\sqrt{6}\lambda+\sqrt{9\xi^{4}-6\left(\lambda^{2}-6\right)\xi^{2}}}{3\left(4+\xi^{2}\right)},\frac{-\sqrt{6}\xi^{2}\lambda+\sqrt{9\xi^{4}-6\left(\lambda^{2}-6\right)\xi^{2}}}{3\xi\left(4+\xi^{2}\right)}\right)$ exists when $\left\{\lambda<-\sqrt{6}~{},~{}\xi\geq\sqrt{6\left(\lambda^{2}-6\right)}\right\}$ or $\left\{\lambda=-\sqrt{6}~{},~{}\xi\neq 0\right\}~{}$ or $\left\{-\sqrt{6}<\lambda<\sqrt{6},~{}\xi<0\right\}$  or $\left\{\lambda>\sqrt{6}~{},~{}\xi<-\frac{1}{3}\sqrt{6\left(\lambda^{2}-6\right)}\right\}$. The parameter of the equation of state is

[TABLE]

where $w_{\phi}<-\frac{1}{3}$ when $\left\{-\sqrt{2}<\lambda<-\sqrt{\frac{2}{3}},-\sqrt{\frac{2\left(2-\lambda^{2}\right)^{2}}{3\lambda^{2}-2}}<\xi<0\right\}$ or $\ \left\{-\sqrt{\frac{2}{3}}\leq\lambda\leq\sqrt{2}~{},~{}\xi<0\right\}$ or $\left\{\lambda>\sqrt{2},~{}\xi<-\sqrt{\frac{2\left(2-\lambda^{2}\right)^{2}}{3\lambda^{2}-2}}\right\}$.  Point $A_{3}$ describes a stable solution only for $\xi<0$ and more specifically $\left\{0<\lambda<\sqrt{6}~{},~{}-\xi_{1}<\xi<0\right\}$ or $\left\{-\sqrt{6}\leq\lambda\leq\sqrt{6}~{},~{}\xi<0\right\}$ or $\left\{\lambda>\sqrt{6}~{},~{}\xi_{2}<\xi<0\right\}$, where $\xi_{1},\xi_{2}$ are the solutions of the algebraic equation

[TABLE]

In Fig. 1 we present the contour plots of the equation of state parameter $w_{\phi}\left(\lambda,\xi\right)$ for points $A_{2}$ and $A_{3}$. Notice that the stable critical points are represented by shaded regions. We observe that points with $\lambda>0$ and $\xi<0$ describe stable accelerated solutions, while it is possible for the EoS parameter to cross the phantom line, namely $w_{\phi}<-1$.

IV.2 Family B

We continue our analysis with the second family of critical points, namely Family B. Here the dynamical system is formed by equations (20), (21) and (23). By including the constraint, the dimension of the system is reduced by one, i.e. from three to two dimensions. We study the general evolution of the dynamical system by considering a general function $\Gamma_{\xi}\left(\xi\right)$ fm1 ; fm3 ; fm4 .

The critical points $\left(x_{p},y_{p},\xi_{p}\right)~{}$of the dynamical system are:

a. Points $B_{1\left(\pm\right)}=\left(A_{1\left(\pm\right)},\xi_{0}\right)$ for which $\xi_{0}$ is a solution of the algebraic equation $\Gamma_{\xi}\left(\xi_{0}\right)=\frac{\lambda}{\sqrt{2}}~{},~{}$or $\xi_{0}=0.$ The points describe the same physical solution as those for $A_{1\left(\pm\right)}$, hence the existence of the critical points is given in section 4.1, however the stability conditions change Specifically, the two eigenvalues are $e_{1}\left(B_{1\left(\pm\right)}\right)=6\pm\sqrt{6}\lambda$ and $e_{2}\left(B_{1\left(\pm\right)}\right)=\pm\sqrt{3}\xi_{0}\Gamma_{\xi}^{\prime}\left(\xi_{0}\right).~{}$ Therefore, $B_{1\left(+\right)}$ is stable when $\lambda<-\sqrt{6}~{},~{}\xi_{0}\Gamma_{\xi}^{\prime}\left(\xi_{0}\right)<0$, while $B_{1\left(-\right)}$ is stable as long as $\lambda>\sqrt{6}~{},~{}\xi_{0}\Gamma_{\xi}^{\prime}\left(\xi_{0}\right)>0$.

b. Point $B_{2}=\left(A_{2}\left(\xi_{0}\right),\xi_{0}\right)$ with $\Gamma_{\xi}\left(\xi_{0}\right)=\frac{\lambda}{\sqrt{2}}$. Again the properties of $B_{2}$ are similar with those of $A_{2}$ (see section 4.1). Concerning stability conditions $B_{2}$, describes an attractor solution when $\xi_{0}\Gamma_{\xi}^{\prime}\left(\xi_{0}\right)>0:$ $\left\{\lambda<0~{},\xi_{0}>\frac{\sqrt{6}}{3}\left|\lambda\right|\right\}$ or $\left\{0<\lambda\leq\sqrt{6}~{},~{}\xi_{0}>0\right\}$.

c. Point $B_{3}=\left(A_{3}\left(\xi_{0}\right),\xi_{0}\right)$ with $\Gamma_{\xi}\left(\xi_{0}\right)=\frac{\lambda}{\sqrt{2}}$. The physical properties of $B_{3}$ are those of point $A_{3}$ (see previous section). $B_{3}$ describes a stable solution for $\xi_{0}\Gamma_{\xi}^{\prime}\left(\xi_{0}\right)<0:\left\{-\sqrt{6}<\lambda\leq 0,~{}\xi_{0}<0\right\}$ or $\left\{\lambda>0,\xi<-\frac{\sqrt{6}}{3}\lambda\right\}$.

d. Point $B_{4}$ with coordinates $B_{4}=\left(-\frac{\lambda}{\sqrt{6}},\sqrt{1-\frac{\lambda^{2}}{6}},0\right)$  exists for $\left|\lambda\right|<\sqrt{6}$. This situation describes a tracking solution of the exponential potential with $\xi=0$. The equation of state parameter reads $w_{\phi}\left(\lambda,\xi\right)=-1+\frac{\lambda^{2}}{3}$, hence we have acceleration when $\left|\lambda\right|<\sqrt{2}$. The eigenvalues of the linearized system are determined to be $e_{1}\left(B_{4}\right)=-3+\lambda^{2}~{},~{}e_{2}\left(B_{4}\right)=\frac{\lambda}{2}\left(\lambda-\sqrt{2}\Gamma_{\xi}\left(0\right)\right),~{}$where $e_{1}\left(B_{4}\right)<0,\ e_{2}\left(B_{4}\right)<0$ for $\left|\lambda\right|<\sqrt{3}$, and $\left|\Gamma_{\xi}\left(0\right)\right|>\frac{\sqrt{2}}{2}\lambda$.

e. Finally, point $B_{5}$ with coordinates $B_{5}=\left(0,1,\xi_{0}\right),~{}\xi_{0}=-\sqrt{\frac{2}{3}}\lambda$ describes a de Sitter solution for which $w_{\phi}=-1$. The eigenvalues of the linearized system are $e_{1}\left(B_{5}\right)=-\frac{3+\sqrt{3\left(3-2\lambda^{2}+2\sqrt{2}\lambda\Gamma_{\xi}\left(\xi_{1}\right)\right)}}{2}~{},~{}e_{1}\left(B_{5}\right)=-\frac{3-\sqrt{3\left(3-2\lambda^{2}+2\sqrt{2}\lambda\Gamma_{\xi}\left(\xi_{1}\right)\right)}}{2}~{}$and thus point $B_{5}$ is an attractor when $\left|\Gamma_{\xi}^{\prime}\left(\xi_{0}\right)\right|<\frac{\lambda}{\sqrt{2}}$. Notice that the de Sitter solution does not exist for the exponential case in the context of the scalar field cosmology which reduces to GR.

IV.2.1 Application

Consider $Y\left(\phi\right)=Y_{1}e^{\left(\nu+\frac{\lambda}{2}\right)\phi},~{}\nu\neq 0$; we calculate that $\phi=\frac{1}{2}\ln\left(\frac{2\xi^{2}}{Y_{1}^{2}\left(2\nu+\lambda\right)^{2}}\right)$ and $\Gamma\left(\xi\right)=\left(\sqrt{2}v+\frac{\lambda}{\sqrt{2}}\right)$.

Therefore, equation (23) is simplified to

[TABLE]

while the possible critical points are now only points $B_{1\left(\pm\right)}~{}$with $\xi_{0}=0$, $B_{4}$ and $B_{5}$. Points $B_{4}$ and $B_{5}$ are attractors when $\left\{\left|\lambda\right|<\sqrt{6},\nu\lambda>0\right\}$ and $\left\{\nu<0~{},~{}0<\lambda\leq-\frac{3}{2\nu}\right\}\cup\left\{\nu>0~{},~{}-\frac{3}{2\nu}<\lambda\leq 0\right\},$ respectively.

In Fig. 2 we present the phase space diagram for the dynamical system in the variables $\left\{x,\xi\right\}$ for two sets of the variables $\lambda~{}$and $\nu$. For $\lambda=-2$ and $\nu=1$ it is clear that the de Sitter universe $B_{5}~{}$is an attractor while for $\lambda=-\sqrt{3}$, $\nu=-\frac{1}{2}$ the unique attractor is the scaling solution $B_{4}$.

In a similar way we continue with the third family of critical points that we considered.

IV.3 Family C

The third family of critical points correspond to the dynamical system (20), (21) and (22) with the constraint condition (19). Recall that in this case $\xi=\xi\left(\lambda\right)$.

The critical points are:

a. Points $C_{1\left(\pm\right)}=\left(A_{1\left(\pm\right)},\lambda_{0}\right)$ with $\Gamma\left(\lambda_{0}\right)=1$ or $\lambda_{0}=0.$ The physical descriptions of the points are those of $A_{1\left(\pm\right)}$. The eigenvalues of the linearized system are calculated to be $e_{1}\left(C_{1\left(\pm\right)}\right)=6-\sqrt{6}\lambda_{0}~{},~{}e_{2}\left(C_{1\left(\pm\right)}\right)=-\sqrt{6}\lambda_{0}^{2}\Gamma^{\prime}\left(\lambda_{0}\right),$. We observe that at least of one of the points $C_{1\left(\pm\right)}$ is stable only for $\Gamma^{\prime}\left(\lambda_{0}\right)>0$ and $\left|\lambda_{0}\right|>\sqrt{6}$.

b. Point $C_{2,3}=\left(A_{2,3},\lambda_{0}\right)$ with $\Gamma\left(\lambda_{0}\right)=1$. The existence conditions and the physical description are the same as those of $A_{2,3}$ (see section 4.1). Because of the nonlinearity of the eigenvalues, the map of $\lambda_{0},\xi\left(\lambda_{0}\right)$ in which the points are stable is presented in Fig. 3.

c. Point $C_{4}=\left(-\frac{\xi}{\sqrt{4+\xi^{2}}},\frac{2}{\sqrt{4+\xi^{2}}},\lambda_{0}\right)$ with $\lambda_{0}=0,$ describes a de Sitter solution, $w_{\phi}\left(C_{4}\right)=-1$, and actually reduces to points $C_{2}$, $C_{3}$, with $\lambda_{0}=0$ respectively. The eigenvalues of the linearized system are given by $e_{1}\left(C_{4}\right)=0~{},~{}e_{2}\left(C_{4}\right)=-3.~{}$ Since $e_{1}=0$ we apply the central manifold theorem in order to decide the stability and we find that $C_{4}$ is always an attractor.

d. Point $C_{5}=\left(B_{5},\lambda_{0}\right)$ with $\lambda_{0}=-\sqrt{\frac{3}{2}}\xi$ is found to be stable when the following condition holds $\left(\Gamma_{\lambda}\left(0\right)-1\right)\left(2+\sqrt{6}\xi^{\prime}\left(0\right)\right)>0\,$. The physical description of $C_{5}$ is the same as that of $B_{5}$.

IV.3.1 Application

Let us consider $U\left(\phi\right)=U_{0}\phi^{n}$ and $Y\left(\phi\right)=Y_{0}\phi^{1+\frac{2}{n}}$ from where we calculate $\phi=\frac{n}{\lambda}$, $\xi=\frac{Y_{0}}{\sqrt{2U_{0}}}\left(2+n\right)=const$ and $\Gamma_{\lambda}\left(\lambda\right)=\frac{n-1}{n}$. Therefore, equation (22) reduces to

[TABLE]

Therefore, the critical points are $C_{1\left(\pm\right)},~{}C_{4}$ and $C_{5}$. As far as stability is concerned we find that points $C_{1\left(\pm\right)}$ are always unstable and point $C_{5}~{}$is stable only when $n\leq-2\xi^{2}$. For the latter case using various values of the free parameters $n,~{}\xi$ we plot in Figs. 4 and 5 the phase space diagram $\left\{x,\lambda\right\}$.

V Scalar field in the presence of matter

In this section we include in our analysis a pressureless matter component with $w_{m}=\frac{p_{m}}{\rho_{m}}=0$. In this case since $0\leq\Omega_{m}\leq 1$, the constraint equation (19) yields $x^{2}+y^{2}\leq 1$. Following the lines of the previous section we study the same family of potentials, namely A, B and C.

V.1 Family A

a. The first Point is $\bar{A}_{0}=\left(0,0\right)$ and the arbitrary parameters $\lambda,~{}\xi$ describe a universe for which $\Omega_{m}=1$, $w_{\mathrm{tot}}=0$. The eigenvalues of the linearized system are determined to be $e_{1}\left(\bar{A}_{0}\right)=-\frac{3}{2}$, $e_{2}\left(\bar{A}_{0}\right)=\frac{3}{2}$, hence the point is always unstable.

b. Points $\bar{A}_{1\left(\pm\right)}$ with eigenvalues $e_{1}\left(\bar{A}_{1\left(\pm\right)}\right)=3\pm\sqrt{\frac{3}{2}}\lambda~{},~{}e_{2}\left(\bar{A}_{1\left(\pm\right)}\right)=3$, from which we infer that points $\bar{A}_{1\left(\pm\right)}$ are unstable points.

c. Point $\bar{A}_{2}$ with eigenvalues $e_{1}\left(\bar{A}_{2}\right)=-3+\frac{4\lambda^{2}+\lambda\sqrt{6\xi^{4}-4\left(\lambda^{2}-6\right)}}{4+\xi^{2}}~{},~{}e_{2}\left(\bar{A}_{2}\right)=-3+\frac{2\lambda^{2}+\lambda\sqrt{6\xi^{4}-4\left(\lambda^{2}-6\right)}}{2\left(4+\xi^{2}\right)}$ describes a stable solution when $\left\{\lambda<-\sqrt{3}~{},~{}\xi>\frac{2\left(\lambda^{2}-3\right)}{\sqrt{6\lambda^{2}-9}}\right\}~{}$or $\left\{-\sqrt{3}\leq\lambda\leq\frac{\sqrt{6}}{2}~{},~{}\xi>0\right\}$ or $\left\{\lambda>\frac{\sqrt{6}}{2}~{},~{}\xi<\frac{6-2\lambda^{2}}{\sqrt{6\lambda^{2}-9}}\right\}$ or $\left\{\lambda=-\sqrt{6},~{}\xi>\frac{2\sqrt{3}}{3}\right\}$.

d. Point $\bar{A}_{3}$ with eigenvalues $e_{1}\left(\bar{A}_{3}\right)=-3+\frac{4\lambda^{2}-\lambda\sqrt{6\xi^{4}-4\left(\lambda^{2}-6\right)}}{4+\xi^{2}}~{},~{}e_{2}\left(\bar{A}_{3}\right)=-3+\frac{2\lambda^{2}-\lambda\sqrt{6\xi^{4}-4\left(\lambda^{2}-6\right)}}{2\left(4+\xi^{2}\right)}$ describes a stable solution when $\left\{\lambda<-\frac{\sqrt{6}}{2}~{},~{}\frac{2\left(\lambda^{2}-3\right)}{\sqrt{6\lambda^{2}-9}}<\xi<0\right\}~{}$or $\left\{-\frac{\sqrt{6}}{2}\leq\lambda\leq\sqrt{3}~{},~{}\xi<0\right\}$ or $\left\{\lambda>\sqrt{3}~{},~{}\xi<-\frac{2\left(\lambda^{2}-3\right)}{\sqrt{6\lambda^{2}-9}}\right\}$ or $\left\{\lambda=\sqrt{6}~{},~{}\xi<-\frac{2\sqrt{3}}{3}\right\}$.

e. Point $\bar{A}_{4}$ with coordinates $\bar{A}_{4}=\left(-\sqrt{\frac{3}{2}}\frac{1}{\lambda},\sqrt{\frac{3}{2}}\frac{\sqrt{\lambda^{2}\left(4+\xi^{2}\right)}-\lambda\xi}{2\lambda^{2}}\right)$ describes a universe where the scalar field mimics the pressureless fluid, i.e. $w_{\phi}=0$. The effective parameter is $w_{\mathrm{tot}}=0$, where $\Omega_{m}=1-\frac{3}{\lambda^{2}}-\frac{3\xi^{2}}{4\lambda^{2}}+\frac{3}{4}\frac{\xi}{\lambda^{3}}\sqrt{\lambda^{2}\left(4+\xi^{2}\right)}$. The point exists when $\left\{\lambda<-\frac{\sqrt{6}}{2}~{},~{}\xi\leq\frac{2\left(\lambda^{2}-3\right)}{\sqrt{6\lambda^{2}-9}}\right\}$ or$~{}\left\{\sqrt{\frac{3}{2}}<\lambda<\sqrt{3}~{},~{}\xi\geq\frac{2\left(\lambda^{2}-3\right)}{\sqrt{6\lambda^{2}-9}}\right\}$ or $\left\{\lambda\geq\sqrt{3}~{},~{}\xi\geq-\frac{2\left(\lambda^{2}-3\right)}{\sqrt{6\lambda^{2}-9}}\right\}$.

In Fig. 6 we show $\left\{\lambda,\xi\right\}$ diagrams for where the critical points $\bar{A}_{2}~{},~{}\bar{A}_{3}$ and $\bar{A}_{4}$ exist and have negative eigenvalues, namely we have stable solutions. Moreover, from Fig. 6 we observe that only one of the critical points $\bar{A}_{2}$, $\bar{A}_{3}$ and $\bar{A}_{4}$ can be stable points of the system.

In Fig. 7 we present the phase space diagrams for the dynamical variables $\left\{x,y\right\}$, using different values of the free parameters $\lambda,\xi$. In particular, we provide three diagrams where in each case one of the critical points $\bar{A}_{2},~{}\bar{A}_{3}$ and $\bar{A}_{4}$ is stable. The dynamical evolution of the cosmological parameters $\Omega_{m}\left(a\right)$ and $w_{tot}\left(a\right)$ are demonstrated in Fig. 8 for the real trajectories presented in Fig. 7.

V.2 Family B

The critical points which correspond to family B are the following:

a. Point $\bar{B}_{0}=\left(\bar{A}_{0},\xi\right)$ where $\xi$ is arbitrary: since $\xi$ is arbitrary, point $B_{0}$ describes a line in the space $\left\{x,y,\xi\right\}$. The physical description is that of point $A_{0}$. The eigenvalues of the linearized system are determined to be $e_{1}\left(\bar{B}_{0}\right)=-\frac{3}{2}~{},~{}e_{2}\left(\bar{B}_{0}\right)=\frac{3}{2}$ and $e_{3}\left(\bar{B}_{0}\right)=0$, from which we conclude that $B_{0}$ describes an unstable solution and $B_{0}$ is a source point.

b. Points $\bar{B}_{1\left(\pm\right)}$ with eigenvalues $e_{1}\left(\bar{B}_{1\left(\pm\right)}\right)=3~{}$, $e_{2}\left(\bar{B}_{1\left(\pm\right)}\right)=3\pm\sqrt{6}\lambda~{},$ $e_{3}\left(\bar{B}_{1\left(\pm\right)}\right)=\pm\sqrt{3}\xi_{0}\Gamma_{\xi}^{\prime}\left(\xi_{0}\right)$, hence the solutions at points $B_{1\left(\pm\right)}$ are unstable.

c. Point $\bar{B}_{2}$ which is found to describe a stable solution if and only if $\xi_{0}\Gamma_{\xi}^{\prime}\left(\xi_{0}\right)>0$. Notice that the parameters $\left\{\lambda,\xi_{0}\right\}$ can be viewed in Fig. 9.

d. Point $\bar{B}_{3}$ which is found to describe a stable solution if and only if $\xi_{0}\Gamma_{\xi}^{\prime}\left(\xi_{0}\right)<0$, while the parameters $\left\{\lambda,\xi_{0}\right\}$ are given in Fig. 9.

e. Point $\bar{B}_{4}$ with eigenvalues $e_{1}\left(\bar{B}_{4}\right)=-3+\frac{\lambda^{2}}{2}$, $e_{2}\left(\bar{B}_{4}\right)=-3+\lambda^{2}$ and $e_{3}\left(\bar{B}_{3}\right)=\frac{\lambda^{2}}{2}-\frac{\sqrt{2}\lambda}{2}\Gamma_{\xi}\left(0\right).$It describes a stable power-law solution when $\left\{-\sqrt{3}<\lambda<0~{},~{}\Gamma_{\xi}\left(0\right)<\frac{\lambda}{\sqrt{2}}\right\}$ or $\left\{0<\lambda<\sqrt{3}~{},~{}\Gamma_{\xi}\left(0\right)>\frac{\lambda}{\sqrt{2}}\right\}$.

f. Point $\bar{B}_{5}$ describes a stable de Sitter universe when $\left|\Gamma_{\xi}^{\prime}\left(\xi_{0}\right)\right|<\frac{\lambda}{\sqrt{2}}.$

g. Point $\bar{B}_{6}$ with coordinates $\bar{B}_{6}=\left(\bar{A}_{4},\xi_{0}\right)$ with $\xi_{0}=0$ or $\Gamma_{\xi}\left(\xi_{0}\right)=0$. Because of the nonlinearity of the eigenvalues, the regions of the parameter space, $\left\{\lambda,\xi_{0}\right\}$, for which point $\bar{B}_{6}$ is an attractor, are shown in Fig. 10.

V.3 Family C

We complete our analysis with the third family of critical points which correspond to the case where $\lambda_{,\phi}\neq const$ and $\xi=\xi\left(\lambda\right)$. Using the dynamical system (20), (21), (22) we find the following critical points $\left(x_{p},y_{p},\lambda_{p}\right)$:

a. Point $\bar{C}_{0}=\left(\bar{A}_{0},\lambda\right)$ where $\lambda$ is arbitrary. This point describes a line of points in the space $\left\{x,y,\lambda\right\}$. The eigenvalues of the linearized system are $e_{1}\left(\bar{C}_{0}\right)=-\frac{3}{2}~{},~{}e_{2}\left(\bar{C}_{0}\right)=\frac{3}{2}$ and $e_{3}\left(\bar{C}_{0}\right)=0$, from which we infer that the current point is a source “line”.

b. Points $\bar{C}_{1\left(\pm\right)}=\left(\bar{A}_{1\left(\pm\right)},\lambda_{0}\right)$ with $\lambda_{0}=0$ or $\Gamma_{\lambda}\left(\lambda_{0}\right)=1$. These points are always sources, because they always have a positive eigenvalue. Indeed, the corresponding eigenvalues are $e_{1}\left(\bar{C}_{1\left(\pm\right)}\right)=3~{},~{}e_{2}\left(\bar{C}_{1\left(\pm\right)}\right)=3+\sqrt{\frac{3}{2}}\lambda_{0}$ and $e_{3}\left(\bar{C}_{1\left(\pm\right)}\right)=\sqrt{6}\lambda^{2}\Gamma_{\lambda}^{\prime}\left(\lambda_{0}\right)$.

c. Points $\bar{C}_{2,3}$ describes a power-law solution with $\Omega_{m}\left(\bar{C}_{2,3}\right)=0$. Due to the nonlinearity of the eigenvalues in Fig. 11 we present the region of the parameters $\left\{\lambda_{0},\xi\left(\lambda_{0}\right)\right\}$ for which $\bar{C}_{2,3}$ are stable.

d. Point $\bar{C}_{4}$ describes a de Sitter universe and the eigenvalues for the linearized system are $e_{1}\left(\bar{C}_{4}\right)=-3~{},~{}e_{2}\left(\bar{C}_{4}\right)=-3$ and $e_{3}\left(\bar{C}_{4}\right)=0$. Since $e_{3}=0$ means that we need to use central manifold theorem: it implies that $\bar{C}_{4}$ is always a future attractor of the dynamical system.

e. Point $\bar{C}_{5}$ with eigenvalues $e_{1}\left(\bar{C}_{5}\right)=-3$, $e_{2,3}\left(\bar{C}_{5}\right)=-\frac{3}{2}\left(1\pm\sqrt{1+\left(\Gamma_{\lambda}\left(0\right)-1\right)\left(\sqrt{6}\right)\xi^{\prime}-2}\right)$ describes a stable de Sitter solution when $\left\{\Gamma_{\lambda}\left(0\right)<1,\xi^{\prime}>\sqrt{\frac{2}{3}}\right\}$ or $\left\{\Gamma_{\lambda}\left(0\right)>1,\xi^{\prime}<\sqrt{\frac{2}{3}}\right\}$.

f. Point \bar{C}_{6}=\left(\bar{A}_{4},\lambda_{0}\right)\,\with $\lambda_{0}=0$ or $\Gamma\left(\lambda_{0}\right)=0$, is stable when $\left\{\sqrt{\frac{3}{2}}<\lambda~{},~{}\xi>\frac{6-2\lambda^{2}}{\sqrt{6\lambda^{2}-9}}\right\}$ or $\left\{\lambda<\frac{\sqrt{6}}{2}~{},~{}\xi<\frac{6-2\lambda^{2}}{\sqrt{6\lambda^{2}-9}}\right\}$.

Overall, from the aforementioned analysis, we conclude that in the case of matter there are two possible de Sitter solutions which can act as attractors for the expansion dynamics.

VI Exact solutions

In the above analysis we have shown that the evolution of the dimensionless dynamical system is always in a three-dimensional phase space, and there is an extra free function which must satisfy constraints in order for the critical point, that is, the solution at the critical point, to exist. In order to understand this more fully let us consider the field equations (15)-(16) and assume that there is no matter source other than the scalar field, i.e. $\rho_{m}=p_{m}=0$.

From system (15)-(16) we find

[TABLE]

Assume now that $\theta=\theta_{0}t^{-1},~{}$which describes a perfect fluid solution, and for $\phi=t^{-1}$ we find that

[TABLE]

Hence, this solution is described by a point in family C. In order to study the stability of the solution we substitute in (15)-(16) $\theta=\theta_{0}t^{-1}+\varepsilon\Theta\left(t\right),~{}\phi=t^{-1}+\varepsilon\Phi\left(t\right)$ for the latter potentials and the solution of the linearized system reveals that functions $\Theta$ and $\Phi$ decay when$~{}\theta_{0}>3$.

However, this it is not the unique power-law solution as is the case in GR. Indeed, by selecting the same expansion rate $\theta=\theta t^{-1}$, but now $\phi=\ln t^{\frac{1}{\alpha}}$, we find

[TABLE]

and this is a solution that is now described by a critical point of family B.

Moreover, when $\theta=\frac{2}{3}\coth\left(\theta_{0}\tau\right)$ and $\phi=\frac{1}{\alpha t}$ we find

[TABLE]

In a similar way we can construct other kinds of solutions for the scalar field with the aether field and other kinds of matter sources. For instance, the latter solution is that of GR with cosmological constant term and a perfect fluid source.

VII Conclusions

In this paper we have performed a detailed analysis of the dynamical evolution of an Einstein-Aether scalar field cosmology in the framework of a spatially-flat FLRW background universe, where the scalar field is coupled to the aether field. The model that we analysed depends on two unknown functions, the first, $U\left(\phi\right),$ corresponds to the scalar field potential, while the second function $Y\left(\phi\right)$ defines the coupling between the scalar field and the aether field. The Friedmann equation is the same with that of Einstein’s GR, while the Raychaudhuri acceleration equation is modified, since the effective pressure term of the scalar field fluid differs from that in GR.

We use expansion-normalised variables to rewrite the field equations as a system of algebraic-differential equations, which contains at most three independent variables. The possible critical points correspond to three different families of solutions. In family A the scalar field potential is exponential $U\left(\phi\right)=U_{0}e^{\lambda\phi}~{}$while the coupling function is given by $~{}Y\left(\phi\right)=Y_{0}+Y_{1}e^{-\frac{\lambda}{2}\phi}$. Family B corresponds to the exponential potential for $U\left(\phi\right),$ while $Y\left(\phi\right)$ is arbitrary, while in family $C$ the scalar field potential $U\left(\phi\right)$ is arbitrary.

When we assume there is no other fluid source in the universe, we find that family A admits four-critical points, while families B and C admit six critical points. On the other hand, when a pressureless matter source is introduced, the maximum number of possible critical points is increased by two for the three families. At this point it is important to mention that family A corresponds to a specific case of the function $V\left(\phi,\theta\right)$ which is determined in Barrow and for which the the field equations admit exact power-law solutions. Moreover, we found that in family B, power-law solutions exist only when $Y\left(\phi\right)$ is approximated locally by the exponential function $Y\left(\phi\right)\simeq e^{-\frac{\lambda}{2}\phi}$. On the other hand, we found that for family $C$ there is a critical point which describes a de Sitter universe as a future attractor. Additionally, the critical points of families A and C reduce to solutions of GR when Y\left(\phi\right)\is constant.

Finally, in order to demonstrate the evolution of the dynamical system, we presented some phase-space diagrams for specific cases as well as a qualitative evolution for the respective cosmological parameters. From the latter it is clear that this specific scalar field type of cosmology approach can describe some key epochs in cosmological evolution.

Acknowledgements.

GP is supported by the scholarship of the Hellenic Foundation for Research and Innovation (ELIDEK grant No. 633). SB acknowledges support by the Research Center for Astronomy of the Academy of Athens in the context of the program “Tracing the Cosmic Acceleration. JDB is supported by the Science and Technology Facilities Council (STFC) of the United Kingdom.