arxiv: v2 [gr-qc] 9 Nov 2018

Transcription

1 Cosmological Models in the Einstein-Aether Theory arxiv: v2 [gr-qc] 9 Nov 2018 M. Campista 1,2, R. Chan 1, M. F. A. da Silva 3, O. Goldoni 3, V. H. Satheeshkumar 1,3,4 and Jaime F. Villas da Rocha 5 1 Coordenação de Astronomia e Astrofísica, Observatório Nacional (ON), Rua General José Cristino, 77, São Cristóvão, CEP , Rio de Janeiro, RJ, Brazil. 2 Departamento de Física Teórica, Instituto de Física, Universidade Federal do Rio de Janeiro (UFRJ), Rio de Janeiro, Brazil. 3 Departamento de Física Teórica, Instituto de Física, Universidade do Estado do Rio de Janeiro (UERJ), Rua São Francisco Xavier 524, Maracanã, CEP , Rio de Janeiro, RJ, Brazil. 4 Centro Brasileiro de Pesquisas Físicas (CBPF), Rua Dr. Xavier Sigaud, 150, Urca, CEP , Rio de Janeiro, RJ, Brazil. 5 Universidade Federal do Estado do Rio de Janeiro (UFERJ), Instituto de Biociências, Departamento de Ciências Naturais, Avenida Pasteur 458, Urca, CEP , Rio de Janeiro, RJ, Brazil. marcelacampista@hotmail.com, chan@on.br, mfasnic@gmail.com, otaviosama@gmail.com, vhsatheeshkumar@gmail.com, jfvroch@pq.cnpq.br Abstract. We study cosmological solutions in Einstein-Aether (EA) theory with an isotropic fluid where the energy density is constant. Assuming a Friedmann-Lemaître- Robertson-Walker (FLRW) metric for the spacetime, we have obtained a contracting, expanding, bouncing and static universes, depending on the parameters of EA theory and the metric. In the case of null energy density, we have obtained a vacuum solution for the spacetime for which the Kretschmann scalar is not null but singular at the initial time. We show that the universe can be only accelerated. This accelerated behavior is due to the existence of dark energy during the evolution of the universe. PACS numbers: h, Nx, Cc, Dm, Bw

2 Cosmological Models in the Einstein-Aether Theory 2 1. Introduction A generally covariant theory in which local Lorentz Invariance (LI) is broken by a dynamical unit timelike vector field u a, often referred to as aether, was newly introduced in 2001 aiming to study the effects of a dynamical preferred frame in gravitation and cosmology [1]. The presence of a preferred frame defined by the timelike unit vector field or aether field breaks LI, which makes Einstein-aether (EA) theory a lowenergy effective theory. Infact, the actionofeatheory isthemost general generally covariant functional of the spacetime metric g ab and aether field u a involving no more than two derivatives. The parameters of theory have been severely constrained using many observational tests such as the primordial nucleosynthesis [2], ultra-high energy cosmic rays [3], the solar system tests [4, 5], binary pulsars [6, 7], and more recently gravitational waves [8, 9]. These results confirm that there exists a family of Einstein-aether theories with small-enough couplings that passes all current observational tests [10]. There is an extensive literature about the EA field equation solutions [11]-[17]. In this article, we analyze cosmological solutions of EA theory with an isotropic fluid where the energy density is constant. We show that for a Friedmann-Lemaître- Robertson-Walker (FLRW) metric, there exist universes that are either contracting or expanding or bouncing, depending on several parameters of the model. In the particular case of null energy density, we have obtained a vacuum solution for which the Kretschmann scalar is not null but singular at the initial time. The paper is organized as follows. In Section II, we present the EA field equations. In Section III, we show the cosmological model field equations. In Section IV, we show the solution of a cosmological model with constant density. In Section V, we obtain a particular vacuum solution in EA theory. In Section VI, we present the energy conditions for the effective aether fluid in GR where we analyze the energy conditions for several different combinations of the values of the parameters β, k and α. Finally, in Section VII, we present the final remarks. 2. The Einstein-Aether Field Equations In this section, we outline the essential ingredients of the EA theory. The metric signature is ( +++) and the units are chosen so that the speed of light defined by the metric g ab is unity. The action of EA is given by, g S = (LEinstein +L aether +L matter) d 4 x (1) and where L Einstein = 1 16πG R (2) L aether = 1 16πG [ Kab mn a u m b u n +λ(g ab u a u b +1)] (3)

3 Cosmological Models in the Einstein-Aether Theory 3 and L matter is the matter Lagrangian, which depends on the metric tensor and the matter field. Here, R is the usual Ricci scalar and K ab mn is defined as K ab mn = c 1 g ab g mn +c 2 δ a m δb n +c 3δ a n δb m c 4u a u b g mn, (4) where the c i are dimensionless coupling constants, λ is a Lagrange multiplier enforcing the unit timelike constraint on the aether, and δ a m δb n = gaα g αm g bβ g βn, (5) In the weak-field, slow-motion limit EA theory reduces to Newtonian gravity with a value of Newton s constant G N related to the parameter G in the action (1) by, ( G = G N 1 c ) 1 +c 4. (6) 2 Note that if c 1 = c 4 the EA coupling constant G becomes the Newtonian coupling constant G N, without imposing c 1 = c 4 = 0. Now, we extract the equations of motion by extremizing the action with respect to the variables involved, u a, g mn. The variation with respect to the Lagrange multiplier λ simply imposes the condition that u a is a unit timelike vector, g ab u a u b = 1. (7) The variation of the action with respect u a leads to [18] where, and a J a b +c 4 a a b u a +λu b = 0, (8) J a m = Kab mn bu n, (9) a a = u b b u a. (10) Thevariationoftheactionwithrespect tothemetricg mn gives thedynamical equations, where G Einstein ab G Einstein ab T aether ab T matter ab = T aether ab = R ab 1 2 g abr +8πGT matter ab, (11) = c [J c (au b) +u c J (ab) J (a c u b) ] 1 2 g abj c d c u d +λu a u b +c 1 [ a u c b u c c u a c u b ]+c 4 a a a b, = 2 g δ( gl matter ) δg ab. (12) Sticking to the convention, we have considered aether on the matter-side of the field equations. Itmight aswell beonthegeometry-sideasthemanifolditselfisnowendowed with a vector field which is reflected in the fact that coupling constant of the EA theory is different from that of GR. So, it is unfortunate if one gets an impression that EA theory is simply GR coupled to vector field because of Eq.(12).

4 Cosmological Models in the Einstein-Aether Theory 4 3. Cosmology in Einstein-Aether Theory In order to study cosmology, we consider a spacetime described by the metric of the FLRW form, given by [ ] dr ds 2 = dt 2 +B(t) kr 2 +r2 dθ 2 +r 2 sin 2 θdφ 2, (13) where k is a constant representing the Gaussian curvature of the space at a given time. It can be 0,+1 or 1. In accordance with Eq.(7), the aether field is assumed unitary, timelike and constant, chosen as u a = (1,0,0,0). (14) We now work out different terms in the field equations Eq.(12). Firstly, the usual geometric part of GR is represented by, G Einstein tt = 3 [ ] k B 2 +Ḃ2, (15) G Einstein rr = 1 [ k 1 kr 2 +Ḃ2 +2B B ], (16) G Einstein θθ = (kr 2 +r 2 Ḃ 2 +2r 2 B B), (17) G Einstein φφ = G Einstein θθ sin 2 θ. (18) The contribution of the aether field to the stress-energy tensor is represented by, T aether tt = 3β 2 T aether rr = T aether θθ T aether φφ Ḃ 2 B β 2(1 kr 2 ) = r2 β 2 Here, the constant β is defined as 2, (19) [Ḃ2 +2B B], (20) [Ḃ2 +2B B], (21) = T aether θθ sin 2 θ. (22) β = c 1 +3c 2 +c 3. (23) It can be proved from the equations (15)-(22) that we can not have an anisotropic pressure fluid. We notice also that the field equations do not depend on the parameter c 4. This result comes naturally without imposing any restriction or any parameter redefinition made by previous authors [2][21]. We call attention to the fact that if β = 0 then c 2 = (c 1 + c 3 )/2, although the field equations be reduced to Einstein s field equations since the EA terms of the field equations are nulled, they still preserve the EA coupling constant. On the other hand, from equations (4) and (6), we can see that the total elimination of the EA parts of the field equations would require the more restrictive condition c 1 = c 2 = c 3 = c 4 = 0.

5 Cosmological Models in the Einstein-Aether Theory 5 Now, as per the matter part, we assume that the universe is filled with an isotropic fluid of a constant energy density ρ = ρ 0 and pressure p. Then, the energy-momentum tensor for the matter is given by T matter ab = (ρ+p)v a v b +pg ab, (24) where v a = δt a is a unit timelike vector representing the aether fluid four-velocity. We are also assuming a comoving frame of reference. The component tt of the field equations can be written as 1 [ ] 6k 2B 2 6Ḃ2 3Ḃ2 β +16πGρ 0 B 2 = 0. (25) The component rr of the field equations can be written 1 [2k 2(1 kr 2 ) +2Ḃ2 +4B B +Ḃ2 β +2B Bβ ] +16πGpB 2 = 0. (26) Using equations (25) and (26) we can obtain the second equation of Friedmann B +2) = 3(β B 8πG (ρ 0 +3p). (27) The Kretschmann scalar for the metric (13) is written as K = 12 [ k 2 B 4 +2kḂ2 +Ḃ4 + B ] 2 B 2. (28) The geometrical radius for the metric (13) is defined as R ge = r B 2 = r B. (29) 4. Solution for an Isotropic Fluid in Einstein-Aether with Constant Density The equation (25) gives us the solutions 6k ( B(t) = ) 3β 6+e 2αwt e αwt, k < 0, (30) 6w(β +2) B(t) = e αwt, k = 0, (31) 6k ( B(t) = 3β +6+e 2αwt e αwt, k > 0, 6w(β +2) (32) where α = ±1 and w = 4 πgρ 0 3(β +2), (33) and β +2 > 0 in order to have w real and not divergent necessarily, to Ḃ 0. An additional static solution can be obtained by the assumption that Ḃ = 0 in equation (25), 6k B(t) = 4, (34) πgρ 0 where we can notice that we must have k > 0 only (since ρ 0 > 0).

6 Cosmological Models in the Einstein-Aether Theory 6 We assume the values k = 1,0,1 in the equations (30)-(32), respectively for the sake of simplicity. Substituting the equations (30)-(32) into (26) we get 3w 2 p = 16( 3β 6+e 2αwt ) 2 Gπ ( 24 24β 6β e 2αwt β +e 2αwt β +72e 2αwt +2e 2αwt + 9e 2αwt β 3 +54e 2αwt β 2) e 2αwt, (35) and p = 3 16πG w2 (β +2), (36) 3w 2 p = 16(3β +6+e 2αwt ) 2 Gπ ( 24+24β +6β e 2αwt β +e 2αwt β +72e 2αwt +2e 2αwt + 9e 2αwt β 3 +54e 2αwt β 2) e 2αwt. (37) For the static solution (34) we get p = 1 3 ρ 0. (38) From the covariant divergence of the EA momentum-energy tensor, we obtain 3(ρ 0 +p)ḃ(t) + ρ = 0. (39) B(t) Applying this additional condition to the solutions (30)-(32) we can conclude that the solutions (30) and (32) are valid only for β = 0 and the solution (31) is the unique one that is valid for any value of β +2 > 0. Thus, we can write 6 ( B(t) = ) 6+e 2αw 0 t e αw0t, k = 1 (40) 12w 0 B(t) = e αwt, k = 0 (41) where B(t) = 6 12w 0 ( 6+e 2αw 0 t ) e αw 0t, k = 1 (42) πgρ0 w 0 = 4, (43) 6 then the pressures (35)-(37) reduce to and 3w 2 ( 0 p = 8( 6+e 2αw 0t 12+36e 2αw 0 t +e ) 2αw 0t e 2αw0t, (44) ) 2 Gπ p = 3 16πG w2 (β +2), (45)

7 Cosmological Models in the Einstein-Aether Theory 7 3w 2 ( 0 p = 8(6+e 2αw 0t 12+36e 2αw 0 t +e ) 2αw 0t e 2αw0t, (46) ) 2 Gπ for k = 1,0,1, respectively. Substituting the equations (40), (42), (43) (for k = ±1) into (26), or substituting the equation (41), (33) (for k = 0) into (26) we always get p = ρ 0, in order to fulfill the continuity equation. Substituting equations (40)-(42) into the metric (13) we get for k = 1,0,1, respectively ds 2 = dt ( 6+e 2αw0t ) 2 e 2αw0t 24w0 [ 2 ] dr 2 1+r 2 +r2 dθ 2 +r 2 sin(θ) 2 dφ 2, (48) (47) ds 2 = dt 2 +e 2αwt[ dr 2 +r 2 dθ 2 +r 2 sin(θ) 2 dφ 2], (49) ds 2 = dt (6+e 2αw0t ) 2 e 2αw0t 24w0 [ 2 ] dr 2 1 r 2 +r2 dθ 2 +r 2 sin(θ) 2 dφ 2. (50) We will assume the values k = 1 in the equation (51), for the sake of simplicity. For the static solution (34) we get ds 2 = dt [ ] dr 2 8πGρ 0 1 r 2 +r2 dθ 2 +r 2 sin(θ) 2 dφ 2. (51) Substituting equations (40)-(42) into the Kretschmann scalar (28) we have 24w0 4 K = ( 6+e 2αw 0t ) ( e 2αw 0 t 864e 2αw0t +e 4αw0t +1296e ) 4αw 0t e 4αw0t, (52) K = 2048π2 G 2 ρ 2 0 3(β +2) 2, (53) 24w0 4 K = (6+e 2αw 0t ) ( e 2αw 0 t +864e 2αw0t +e 4αw0t +1296e ) 4αw 0t e 4αw0t, (54) for k = 1,0,1, respectively. For the static solution (34) we get K = π2 G 2 ρ 2 0. (55) We can notice that K is singular only for the case k = 1 at t sing = 1 2αw 0 ln(6). (56)

8 Cosmological Models in the Einstein-Aether Theory 8 Substituting equations (40)-(42) into the geometric radius (29) we get for k = 1, 0, 1, respectively 6 ( R ge = r ) 3(β +2)+e 2αw 0 t e αw 0 t, (57) 12w 0 R ge = re 2αwt, (58) 6 ( R ge = r 3(β +2)+e 2αw 0 t e αw0t, 12w 0 (59) since β +2 > 0 and w > 0. For the static solution (34) we get R ge = r 6 4 πgρ 0. (60) We plot in Figure 1 the geometrical radii for β = 3,0, 3/2, k = 1,0,1, α = ±1 and πgρ 0 = 1. We can notice that the universe can contract after/or expand depending of parameters β, k and α. We can see also the the universe has a Big Bang event for k = 1. The acceleration is defined in terms of the derivative of R ge as a(t) = R ge r = B B B = 8πG 3(β +2) (ρ 0 +3p) B (61) We plot in Figure 2 the acceleration for β = 3,0, 3/2, k = 1,0,1, α = ±1 and πgρ 0 = 1. We can notice that the universe can contract after/or expand depending of parameters β, k and α. We can see also the the universe has a Big Bang event for k = 1. Calculating the acceleration parameter using (61) for any value of β, α = ±1 and πgρ 0 = 1 we have that the universe is always in accelerating (a(t) > 0). 5. Vacuum Solution in Einstein-Aether Theory A vacuum solution is obtained substituting ρ 0 = 0 into equation (25). Solving for B(t), we get 2k B(t) = β +2 t+b 0, (62) where b 0 is a constant of integration and k < 0. We assume the value k = 1 in the equation (62), for the sake of simplicity. We can rewrite this equation considering a rescaling of initial time, B(t) = where t 0 = b 0 β β +2 (t t 0), (63)

9 Cosmological Models in the Einstein-Aether Theory 9 Figure 1. Geometrical radius R ge /r for β = 3,0, 3/2, k = 1,0,1, α = ±1 and πgρ 0 = 1. The geometrical radii for β = 0 and different values of k = 1,0,1 correspond to the solutions (48) and (50), respectively (the two top figures). The geometrical radii for k = 0 and different values of β = 3,0 3/2 correspond only to the solution (49) (the two bottom figures). and Substituting equation (63) into equations (26) and (28) we have p = ρ 0 = 0, (64) K = 3β2 (t t 0 ) 4. (65) Notice that the Kretschmann scalar is not null unless β = 0. Moreover it presents a temporal singularity at t = t 0, although we do not have any matter in the spacetime. This is an unexpected result because we would expect a flat spacetime instead of a curved spacetime. However, this seems an important result because of the lack of a

10 Cosmological Models in the Einstein-Aether Theory 10 Figure 2. Acceleration a(t)/r for β = 3,0, 3/2, k = 1,0,1, α = ±1 and πgρ 0 = 1. The acceleration for β = 0 and different values of k = 1, 0, 1 correspond to the solutions (48) and (50), respectively (the two top figures). The acceleration for k = 0 and different values of β = 3,0 3/2 correspond only to the solution (49) (the two bottom figures). general exterior vacuum solution in EA theory, even with a static aether field, up to now [22]. Substituting equation(63) into the metric(13) and into the geometrical radius (29), assuming t 0 = 0, we get and ds 2 = dt 2 + 2t2 β +2 [ ] dr 2 1+r 2 +r2 dθ 2 +r 2 sin(θ) 2 dφ 2, (66) 2 R ge = r t, (67) β +2 (see Figure 3). Thus, we have a cosmological solution for a spacetime that apparently

11 Cosmological Models in the Einstein-Aether Theory 11 Figure 3. The geometrical radius R ge /r of the vacuum solution for β = 3,0, 3 2, assuming k = 1 and t 0 = 0, for the equation (67). Table 1. Summary of the Analysis of All Possible Solutions Equation k K ρ 0 Spacetime EA-Static Universe (p = ρ 0 ) EA-FRLW (Steady State Universe) (β 0) 0 (p = ρ 0 ) EA-FRLW (Steady State Universe) 0 (β = 0) 0 (p = ρ 0 ) EA-FRLW (Steady State Universe) (p = ρ 0 ) EA-FRLW (Steady State Universe) (β 0) = 0 EA-Schwarzschild-type = 0 (β = 0) = 0 Minkowski (Riemann Tensor 0) there is vacuum, but with a constant aether field velocity. It seems that we can interpret this result as the evolution of the spacetime itself, like de Sitter spacetime. In the following Sections, we show an interpretation of the effective aether field contribution in terms of GR. 6. Energy Conditions in GR for the Effective Aether Fluid with ρ 0 > 0 The energy conditions have not been completely established in EA theory yet, except only for the the energy density that must be non negative [20]. Thus, next we will analyze the energy conditions of GR theory assuming the geometrical part of the EA field equations (G aether ab ) and the EA energy-momentum tensor of a perfect fluid as an effective energy-momentum tensor. In this way we can interpret our solutions (48)-(50) in the context of GR theory. Then, let us now consider Tab aether as part of an effective

12 Cosmological Models in the Einstein-Aether Theory 12 energy-momentum tensor Tab e given by T e ab = Gaether ab +8πGT matter ab, (68) and calculate the energy conditions as in GR. In order to simplify our analysis, we will assume that GR coupling constant 8πG N = 1, for the sake of simplicity. The energy condition for ρ 0 > 0 must be fulfilled [20]. Let us assume the energy-momentum tensor for the aether field as an isotropic fluid with an energy density ρ e, pressure p e, given by T e ab = (ρ e +p e )w a w b +p e g ab, (69) where w a = δ a t is a unit timelike vector representing the fluid velocity. We are also assuming a comoving reference. Thus, and Ttt e = Tmatter +Ttt aether = 8πGρ 0 3β 2 T e rr = T matter rr = T e θθ = Tmatter θθ 1 1 kr 2 +Trr aether = [ 8πGpB 2 + β 2 +T aether θθ Ḃ 2 B 2 = ρ e, (70) (Ḃ2 +2B B) ] = p e B 2 = 8πGpr 2 B 2 + r2 β 2 1 kr2, (71) [Ḃ2 B] +2B = p e r 2 B 2, (72) T e φφ = Te θθ sin2 θ. (73) Using the equations (70) and (71) we get an equation of state giving by ) (Ḃ ρ e +p e = 8πG(ρ 0 +p)+β d, (74) dt B ρ e +3p e = 8πG(ρ 0 +3p)+3β ( ) B. (75) B Now, we will studied the energy conditions in GR of this fluid, using the same spacetime obtained in Section III (see Figures 4, 5, 6 and 7). We will analyze the energy conditions ρ e 0, EC 1 = ρ e +p e 0, EC 2 = ρ e p e 0 and EC 3 = ρ e + 3p e 0 for different combinations of the parameters. The matter classification follows the reference [19]. The static solution (51) fulfills all the energy conditions. We show in Figure 4 the energy conditions EC 1, EC 2 and EC 3, assuming k = 0, β = 3,0, 3/2, α = 1 and πgρ 0 = 1, for the solution (49). We show in Figure 5 the energy conditions EC 1, EC 2 and EC 3, assuming k = 0,β = 3,0, 3/2, α = 1 and πgρ 0 = 1, for the solution (49). From the Figures 4 and 5 we can conclude that we have only dark energy and accelerated universe, since EC 3 is always negative, ρ e is always positive and EC 1 and EC 2 are always non negative.

13 Cosmological Models in the Einstein-Aether Theory 13

14 Cosmological Models in the Einstein-Aether Theory 14

15 Cosmological Models in the Einstein-Aether Theory 15 Figure 6. This is a plot of the energy conditions EC 1 = ρ e +p e, EC 2 = ρ e p e and EC 3 = ρ e +3p e, assuming k = ±1, β = 3,0, 3/2, α = ±1, πgρ 0 = 1, for the solutions (48) and (50). We show in Figure 6 the energy conditions EC 1, EC 2 and EC 3, assuming k = ±1,β = 3,0, 3/2, α = ±1 and πgρ 0 = 1, for the solutions (48) and (50). We show in Figure 7 the energy conditions EC 1, EC 2 and EC 3, assuming k = 1,β = 3,0, 3/2 and πgρ 0 = 1, for the vacuum solution (66). From the Figures 7 we can conclude that we have standard energy for β = 3/2 and β = 0. On the other hand, for β = 3, we have attractive phantom energy in the time interval t 1 t t 3, we have nonphysical energy in the time interval t 3 t t 4 and we have standard energy in the time intervals t < t 1 and t > t 4. Where t 1 = 3/4 and t 2 = 3/4 are the roots of ρ e = 0 and where t 3 = 6/4 and t 4 = 6/4 are the roots of EC 1 = Conclusions We have studied cosmological solutions in EA theory with an isotropic fluid where the energydensityisconstant. WehaveassumedaFLRWmetricforthespacetime. Wehave obtained a contracting, expanding or bouncing universe, depending on the parameters of the EA and the metric. The Kretschmann scalar has singularities at t sing = 1 2αw 0 ln(6) πgρ only for k = 1 where w 0 = 4 0, meaning that this is unique model with a Big 6 Bang event in its evolution. In the particular case of null energy density, we have obtained a vacuum solution for the spacetime, for which the Kretschmann scalar is not null but singular at the initial time. We have also studied the contribution of the aether field equations, considering it as an effective source in GR. We have shown that the universe can be only accelerated. This accelerated behavior can be interpreted as due

16 Cosmological Models in the Einstein-Aether Theory 16

17 Cosmological Models in the Einstein-Aether Theory 17 Table 2. Summary of Energy Conditions in GR of the Effective Aether Fluid for the Dynamical Solutions. Here, EC 1 = ρ e + p e, EC 2 = ρ e p e and EC 3 = ρ e + 3p e. We have assumed πgρ 0 = 1 in order to calculate t 1 and t 2, solving ρ e = 0, where t 1 < t 2 and in order to calculate t 3 and t 4, solving EC 1 = 0, where t 3 < t 4. Figure k β α ρ e EC 1 EC 2 EC 3 Matter > 0 = 0 > 0 < 0 dark energy > 0 = 0 > 0 < 0 dark energy 0-3/2 1 > 0 = 0 > 0 < 0 dark energy > 0 = 0 > 0 < 0 dark energy > 0 = 0 > 0 < 0 dark energy 0-3/2-1 > 0 = 0 > 0 < 0 dark energy 6 ±1 0 ±1 > 0 = 0 > 0 < 0 dark energy -1 3 > 0 > 0 > 0 < 0 dark energy (t 1 < t < t 2 ) -1 3 > 0 < 0 > 0 > 0 attractive phantom energy (t 2 < t < t 3 ) < 0 < 0 > 0 > 0 nonphysical energy (t 3 < t < t 4 ) -1 3 > 0 > 0 > 0 > 0 standard energy (t < t 1 or t > t 4 ) -1 0 > 0 > 0 > 0 > 0 standard energy -1-3/2 > 0 > 0 > 0 > 0 standard energy to the existence of dark energy during the evolution of the universe. This result is expected in the cases described by the metrics (48)-(50), since the fluid corresponds to dark energy. For the static solution (51) the matter fluid satisfies marginally the strong energy conditions and therefore we can attribute a repulsive characteristic of the EA contribution. Acknowledgments The financial assistance from FAPERJ/UERJ (MFAdaS) is gratefully acknowledged. The author (RC) acknowledges the financial support from FAPERJ (no.e- 26/ /2000, E-26/ /2002 and E-26/ /2006). MFAdaS and RC also acknowledge the financial support from Conselho Nacional de Desenvolvimento Científico e Tecnológico - CNPq - Brazil. The author (MFAdaS) also acknowledges the financial support from Financiadora de Estudos e Projetos - FINEP - Brazil (Ref. 2399/03). VHS gratefully acknowledges the financial support from FAPERJ though Programa Pós-doutorado Nota 10. VHS also thanks Jailson Alcaniz for the hospitality at Observatório Nacional (ON), and José Abdalla Helayäl-Neto and Sofiane Faci at Centro Brasileiro de Pesquisas Físicas (CBPF). (MC) thanks ON and UFRJ for their hospitality. References [1] T. Jacobson and D. Mattingly, Phys. Rev. D 64, (2001) [arxiv: ].

18 Cosmological Models in the Einstein-Aether Theory 18 [2] S. M. Carroll and E. A. Lim, Phys. Rev. D 70, (2004)the [arxiv: ]. [3] J. W. Elliott, G. D. Moore and H. Stoica, JHEP 0508, 066 (2005) [arxiv: ]. [4] C. Eling and T. Jacobson, Phys. Rev. D 69, (2004) [arxiv: ]. [5] M. L. Graesser, A. Jenkins and M. B. Wise, Phys. Lett. B 613, 5 (2005) [arxiv: ]. [6] B. Z. Foster, Phys. Rev. D 76, (2007) [arxiv: ]. [7] K. Yagi, D. Blas, E. Barausse and N. Yunes, Phys. Rev. D 89, no. 8, (2014) Erratum: [Phys. Rev. D 90, no. 6, (2014)] Erratum: [Phys. Rev. D 90, no. 6, (2014)] [arxiv: ]. [8] Y. Gong, S. Hou, D. Liang and E. Papantonopoulos, Phys. Rev. D 97, (2018) [arxiv: ]. [9] J. Oost, S. Mukohyama and A. Wang, Phys. Rev. D 97, (2018) [arxiv: ]. [10] B. Z. Foster and T. Jacobson, Phys. Rev. D 73, (2006) [arxiv: ]. [11] Barrow, J.D., Phys. Rev. D 85, (2012) [12] Solomon, A.R. and Barrow, J.D., Phys. Rev. D 89, (2014) [13] Ding, C., Wang, A., Wang, X., Phys. Rev. D 92, (2015) [14] Ding, C., Wang, A., Wang, X., Zhu, T., Phys Rev D (2016) [15] Ho, F.H. Zhang, S.J. Liu, H.S. Wang, A., Phys. Let. B 782, 723 (2018) [16] Oost, J., Bhattacharjee, M, Wang, A., [arxiv: ] [17] Bhattacharjee, M. Mukohyama, S., Wan, M.B., Wang, A., Phys. Rev. D 98, (2018) [18] D. Garfinkle, C. Eling and T. Jacobson, Phys. Rev. D 76, (2007) [19] R. Chan, M. F. A. da Silva and J. F. Villas da Rocha, Mod. Phys. Lett. A, 24, 1137 (2009) [20] D. Garfinkle and T. Jacobson, Phys. Rev. Lett. 107, (2011) [21] A. A. Coley, G. Leon, P. Sandin and J. Latta, JCAP 12, 010 (2015) [22] D. Vernieri and S. Carloni, (2017) [arxiv:1706:06608v1] [23] T. Jacobson, Phys. Rev. D 89, (2014)