SciELO - Scientific Electronic Library Online

 
vol.53 número1The Iras 08589−4714 Star-Forming regionLight curve and orbital period analysis of VX LAC índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista mexicana de astronomía y astrofísica

versión impresa ISSN 0185-1101

Rev. mex. astron. astrofis vol.53 no.1 Ciudad de México abr. 2017

 

Articles

New cosmological solutions in massive gravity theory

S. S. A. Pinho1 

S. H. Pereira1 

E. L. Mendonça1 

1Faculdade de Engenharia de Guaratinguetá, UNESP, Brazil.


Abstract

In this paper we present some new cosmological solutions in massive gravity theory. Some homogeneous and isotropic solutions correctly describe accelerated evolutions for the universe. The study was realized considering a specific form to the fiducial metric and found different functions and constant parameters of the theory that guarantee the conservation of the energy momentum tensor. Several accelerating cosmologies were found, all of them reproducing a cosmological constant term proportional to the graviton mass, with a de Sitter type solution for the scale factor. We have also verified that when the fiducial metric is close to the physical metric the solutions are absent, except for some specific open cases.

Key Words: cosmology: theory - dark energy - dark matter

Resumen

En este trabajo presentamos algunas nuevas soluciones en la teoría de gravedad masiva. Algunas soluciones homogéneas e isotrópicas describen correctamente evoluciones aceleradas de nuestro universo. Hicimos la investigación con formas específicas de la métrica fiducial y sin embargo encontramos en la teoría distintas funciones y para ́metros constantes que garantizan la conservación de energía y momento. Se encontraron varias cosmologías y en todas se constata la reproducción de un término de constante cosmológica que es proporcional a la masa del gravitón con una solución del tipo Sitter para el factor de escala. También hemos comprobado que no hay soluciones cuando la métrica de referencia esta ́ cerca de la métrica física, con excepción de algunos casos específicos, aunque abiertos.

1. Introduction

The idea to give a non-null mass to the graviton is an old one in the history of physics (Fierz & Pauli 1939; Boulware & Deser 1972). The proof of the existence of a non-linear generalization of the so-called massive gravity theory is a problem that has stimulated several studies in last years (Hassan & Rosen 2012, 2011; Hassan et. al. 2012; Rham & Gabadadze 2010; Hinterbichler 2012; Arkani-Hamed et. al. 2003; Volkov 2012a,b; Kobayashi 2014; Gümrükçuöglu et. al. 2012, 2011; de Rham et. al. 2014; Dubovsky 2004; de Rham et. al. 2014; Creminelli et. al. 2005), since the current observations of supernovae of type Ia (SNIa) (Riess 1998; Amanullah et. al. 2010), cosmic microwave background (CMB) radiation (Komatsu et. al. 2011; Ade et. al. 2014) and Hubble parameter data (Farooq & Ratra 2013; Sharov & Vorontsova 2014) indicate an accelerated expansion of the universe. Massive gravitons could perfectly mimic the effect of a cosmological constant term, rendering the theory a good alternative to the ΛCDM model of cosmology, which is plagued with several fundamental issues related to the cosmological constant in a Friedmann-Lemaıtre-Robertson-Walker (FLRW) background (Weinberg 1989).

In general, massive gravity theories have ghost instabilities at the non-linear order, known as BoulwareDeser (BD) ghosts (Boulware & Deser 1972). A procedure recently outlined in Rham & Gabadadze (2010) successfully obtained actions that are ghost-free for the fully non-linear expansion, resulting in theories sometimes called dRGT (de Rham-Gabadadze-Tolley) models, which were first shown to be BD ghost free in Hassan & Rosen (2012, 2011); Hassan et. al. (2012). Massive gravity is constructed with an additional metric f μν (called fiducial metric) together with the physical metric g μν . When the fiducial metric is also endowed with a dynamic, the theory is called bimetric gravity or just bigravity. In this theory the problems of ghosts and instabilities are absent, but a new mass scale related to the fiducial metric is also present in addition to the graviton mass. For this reason we prefer to study the massive gravity version where just the graviton mass is present, and the reference metric is assumed to be of non-FLRW type, admitting non-isotropic and non-homogeneous forms. In some sense we study cases similar to those already presented in Volkov (2012a), but several other solutions are also analysed, including an approach based just on the conservation of the energy-momentum tensor. In some sense we have also generalized the cases of Volkov (2012a), since we have considered a physical metric of Friedmann type.

It is well known that when the fiducial metric is assumed to be flat, isotropic flat and closed FLRW cosmologies do not exist even at the background level. On the other hand, isotropic open cosmologies exist as classical solutions, but there are also perturbations that are unstable (Gümrükҫuög˘lu et. al. 2012, 2011). Quantum cosmological models in the framework of nonlinear massive gravity have also been studied in Vakili & Khosravi (2012), but some constraints present in the theory prevent a flat and closed FLRW solution, although this is not the case for the open universe. Cosmological constraints with observational data for massive gravity theories have also been studied recently (Gannouji et. al. 2013; Pereira et. al. 2016), allowing the determination of free parameters of the theory.

The paper is organized as follows. In § 2 we present the general massive gravity theory. In § 3 we present several cases with accelerating solutions based on the knowing of the fiducial metric. In § 4 we present solutions based on the conservation of the energy momentum tensor. We conclude in § 5.

2. Massive gravity theory without ghosts

The covariant action in the massive gravity theory according to Rham & Gabadadze (2010) can be written as:

S=18πGd4x-g-R2+m2UK+Smat (1)

with S(mat) representing the ordinary matter content (radiation, baryons, dust, etc) and

UK=U2+α3U3+α4U4 (2)

acting like a potential due to the graviton mass term, m. In (2) α3 and α4 are constant parameters and the functions U2, U3 and U4 are given by:

U2=12!k2-k2 (3)

U3=13!k3-3kk2+2k3 (4)

U4=14!k4-6k2k2+3k22+8kk3-6k4 (5)

where traces are represented by [K] and superior orders of K n (n = 2,3,4) are Kμ αKα ν, Kμ αKαβKβ ν and Kμ αKα βKβ γKγ ν, respectively.

We introduce K as:

kυμ=δνμ-Υνμ (6)

with γ given by:

γνμ=gμαfαν (7)

where g stands for the physical metric and f represents the fiducial metric. The matrix represented by kνμ, and consequently the mass term, breaks the general covariance due to the presence of the fiducial metric. However thanks to the well-known Stückelberg trick (Arkani-Hamed et. al. 2003) one can think f μν as a covariant tensor field which can be constructed as follows:

fμν=ηABδμXAδνXB (8)

where XA are a set of four fields which transform as scalars under a general coordinate transformation of spacetime and are called Stückelberg scalars. At the same time ηAB = diag[1, −1, −1, −1]. Notice that when the physical metric and the fiducial metric coincide we have kνμ identically zero, which cancel the contribution of massive gravity. In order to have a non-trivial contribution from the graviton mass term we must consider a fiducial metric slightly different from the physical one.

The variation of the action with respect to gμν results in the Einstein equations:

Gνμ=m2Tνμ+8πGT(mat)νμ (9)

where T(mat)νμ corresponds to the energy-momentum tensor of the ordinary matter,

T(mat)νμ= 18πGdiagρt,-ρt,-ρt,-ρt (10)

and Tνμ represents the contribution to the energy-momentum tensor coming from the graviton mass term,

Tμν=2δUδgμν-Ugμν (11)

Substituting (2) in (11) one finds:

Tνμ=Υαμkνα-kδνα-α3ΥαμU2δνα-kkνα+k2να+α4ΥαμU3δνα-U2kνα-kk2να+k3να-Uδνμ (12)

By considering that there is no interaction between ordinary matter and the massive graviton we have that the energy-momentum tensor must be conserved separately,

μTμν(mat)=0 μTμν (13)

which will impose conditions on the Stückelberg fields and also on the α parameters.

In which follows we will considerer just the graviton mass term contribution to the energy momentum tensor, present in the fiducial metric. The matter part can be added in the standard way.

3. Stückelberg’s analysis: fiducial metric

When we are looking for proper Stückelberg fields the starting point is to define the physical and the fiducial metric. In order to reproduce cosmological results we posit that the physical metric must be homogeneous and isotropic, i.e.:

dsg2=N(t)2dt2-αt2dr2-αt2fr2(dθ2+sin2θdϕ2) (14)

where N(t) is known as the lapse function and f(r) can be equal to r, sinh(r) or sin(r) depending on what kind of universe we are dealing with: flat (k = 0), open (k = −1) or closed (k = +1), respectively. As in Volkov (2012a), we assume a generalized form for the fiducial metric, spherically symmetric:

dsf2=T˙t,r2-U˙t,r2dt2+2T˙t,rT't,r-U˙t,rU't.rdtdr+T't,r2-U't,r2dr2-Ut,r2(dθ2+sin2θdϕ2) (15)

Here the functions T (t, r) and U (t, r) are closely related to the Stückelberg fields. Besides, dots and primes represent derivatives with respect to t and r, respectively. The ghost-free massive gravity introduced in the last section can be rewritten in terms of such functions through function (7). We bind both metrics (14) and (15) according to (7) and then we are able to describe the dRGT theory using the functions T (t, r) and U (t, r).

From now on we assume the condition:

T˙t,rT't,r=U˙(t,r)U'(t,r) (16)

which guarantees that fiducial metric (15) is diagonal and consequently homogeneous and symmetric, depending on the functions T(t,r) and U(t,r). With such a constraint one can check, following (7), that:

(Υνμ)2=A(t,r)N(t)20000B(t,r)α(t)20000U(t,r)2α(t)2f(r)20000U(t.r)2α(t)2f(r)2 (17)

where we have defined A(t,r)= T˙(t,r)2U˙(t,r)2 and B(t,r)=U′(t,r)2T′(t,r)2. We assume that thereis a square root of (17) given by the ansatz:

Υνμ=A(t,r)±N(t)2Ft,r00Gt,rB(t,r)±α(t)0000U(t,r)α(t)f(r)0000U(t,r)α(t)f(t,r) (18)

where G(t,r) and F(t,r) are a priori undefined functions. But by squaring (18) we see that such functions cannot assume any form. Therefore we have three different cases to analyse. The first one concerns a diagonal fiducial metric (F(r,t) = G(t,r) = 0) and the other two are non-diagonal (F(r,t) ̸= 0 and G(t,r) = 0 or F(r,t) = 0 and G(t,r) ̸= 0).

In order to check the viability of such choices for the Stückelberg fields we will consider the conservation of the energy-momentum tensor, ∇μTνμ= 0, for each component. Due to the specific form of (18) we have ∇μT2μ= ∇μT3μ= 0. The condition of conservation of the T0μ component will be the same for all the following cases. It is given by:

μT0μ=-3αtAt,rNt3+3α3+α4±2αt2-2Ut,r+Bt,rfrAt,rfrNt1+2α3+α4+1α(t)3(-Ut,r+2B(t,r)f(r))At,rU(t,r)ft2N(t)(α3+α4)da(t)dt±12at-4δUt,rδtBt,r+δBt,rδtfrfrBt,r3+3α3+α4+1at2-2Ut,rδUt,rδtBt,r+Ut,rδBt,rδtfr+2Bt,rδUt,rδtfrfr2Bt,r1+2α3+α4±12a(t)3Ut,rδBt,rδt+4Bt,rδUt,rδtU(t,r)f(t)2B(t,r)(α3+α4) (19)

In the following sections we present the analysis for the conservation of the energy momentum tensor for several cases, according to the values of α3 and α4 or A(t,r), B(t,r), F(t,r) and G(t,r), depending on which case we are analysing.

3.1. Diagonal Metric

Here we consider that γ, given in (18), is diagonal (G(t,r) = F(r,t) = 0) and the functions A(t,r) and B(t,r) are independent. However, we still have the constraint (16). There are several functions which respect this constraint, but we are interested in cases for which ∇μT0μ = 0 [given by (19)], and ∇μT1μ = 0, which is given by:

μT1μ=12δA(t,r)δrNtAt,r3+3α3+α4+2αtdfrdrBt,r+δUt,rδtfr3α2+3α3+α4±1αt2At,rδUt,rδr+2dfrdrBt,rAt,r+Ut,rδAt,rδrfrNtAt,r1+2α3+α4+2at2dfrdrBt,r+δUt,rδrfr21+2α3+α4+12at24At,rδUt,rδr+Ut,rδAt,rδr+4dfrdrBt,rAt,rU(t,r)fr2NtAt,r(α3+α4) (20)

We now analyse some different cases.

3.1.1. A(t,r) = 0, B(t,r) = 4/f(r) 6 and U(t,r) = 1/f(r) 2

The conservation condition (19) gives an expression for T(t,r), which depends on what kind of universe we are dealing with. One can check that T(t,r) is a constant for f(r) = r (flat universe), T(t,r) = ±2coth(r)2 for f (r) = sinh(r) (open universe) and T(t,r) = ±2i cot(r)2 for f(r) = sin(r) (closed universe). Therefore, the equations of motion obtained by varying the action with respect to a(t) and N(t) are, respectively:

4a¨(t)a(t)+2a˙(t)2a(t)2+m2a(t)2f(r)61+2α3+α4+2ka(t)2-m26+4α3+α4=0 (21)

and

6a¨(t)2a(t)2+3m2a(t)2f(r)6+1+2α3+α4+6ka(t)2-2m2at3fr9α3+α4-m26+4α3+α4=0 (22)

These are the FLRW equations obtained for the particular choice of A(t,r), B(t,r) and U(t,r). The presence of the anisotropic terms with f(r) means that the solutions are of little interest for cosmology. However if we choose α3 = −1 and α4 = 1, the solutions turn out to be isotropic, with the presence of a cosmological constant-like term, related to the massive graviton mass,

2a(t)¨a(t)+a˙(t)2a(t)2+ka(t)2-32m2=0 (23)

a˙(t)2a(t)2+ka(t)2-m22=0 (24)

Such solutions are accelerating and could perfectly reproduce the cosmological observations.

3.1.2. U(t,r)=0

In this case we have T(t,r) = constant, which leads to A(t,r) = B(t,r) = 0. The energy-momentum tensor is always conserved and both the fiducial metric and γ are null. Although this seems to be a pathological case, it just reproduces a cosmology with a cosmological constant, represented by the m 2 term, which comes from U(t,r) in the last term of (12). The equations of motion are:

2a¨(t)a(t)+a˙(t)2a(t)2+ka(t)2-m226+4α3+α4=0 (25)

a˙(t)2a(t)2+ka(t)2-m226+4α3+α4=0 (26)

3.1.3. A(t,r) = W(t) 2 , B(t,r) = K(t) 2 and U(t,r) = −K(t)f(r)

When one tries to impose conservation of the energy-momentum tensor it is found that the universe must be open (k = −1), and the relation between W(t) and K(t) is described by

Kt=tWt¨Nt~a˙t¨a˙(t¨)dt~ (27)

However, the relation (27) together with the conditions for T(t,r) results in a a˙(t) = −N(t), which has no cosmological interest. As can be seen from equation (27), when W(t) = N(t), K(t) would be equal to −a(t) and this results in a standard Friedmann metric. This is due the fact that there is no cosmology when the fiducial metric mimics the Friedmann one (Rham & Gabadadze 2010).

3.1.4. At,r=N(t)2a(t)2B(t,r)

With this condition we intend to verify the compatibility of cosmology and massive gravity for Friedmannlike metrics. Such compatibility has been already checked in de Rham et. al. (2014) and states that only an open universe could afford both theories. In some sense this analysis is also present in Volkov (2012a), but the scale factor a(t)2 present in the time component of the physical metric used in Volkov (2012a) does not represents a true FLRW metric. Here we have addressed a similar analysis but we have taken N(t)2 in (14), and after doing calculations the function N(t) is assumed equal to one in order to recover the FLRW background.

Here we present another approach to analyse a solution close to the Friedmann metric. Our choice is to find solutions of T(t,r) and U(t,r) and then try to manipulate them in order to get conservative cases. Besides, we also find results that conserve the energy-momentum tensor like in the previous subsection. By confronting both results we conclude that there is really something special in the theory that excludes solutions for an isotropic and homogeneous fiducial metric.

There are two constraints on the functions T(t,r) and U(t,r) which guarantee diagonality of γ2. By making Υ002Υ112 in equation (17) we get At,r=N(t)2a(t)2B(t,r), which gives:

T˙(t,r)2-U˙t,r2=-N(t)2a(t)2T't,r2-U'(t,r)2 (28)

After substituting (16) in equation (28) one can find

U˙t,r=±N(t)a(t)T'(t,r) (29)

and

T˙t,r=±N(t)a(t)U'(t,r) (30)

By deriving (29) with respect to t and (30) with respect to r one can combine the results to get:

U¨t,r=a(t)N(t)ddtN(t)a(t)U˙t,r+N(t)2a(t)2U''(t,r) (31)

It is possible to check that T(t,r) has the same solutions of U(t,r) since that function satisfies an equation equivalent to (31). We have chosen U(t,r) = R(r)Z(t) in order to separate the functions of r and t in (31). For this solutions C will stand for the constant of the Fourier method.

When C = 0 the solution is:

Ut,r=a1tN(t~)a(t~)dt~+a2a3+a4r (32)

where a1, a2, a3 and a4 are constants.

For a positive constant (C > 0) one obtains:

Ut,r=b1sinhCtN(t~)a(t~)dt~+b2coshCtN(t~)a(t~)×b3sinhCr+b4coshCr (33)

where b1, b2, b3 and b4 are constants.

For a negative constant (C < 0):

Ut,r=c1sin|C|tN(t~)a(t~)dt~+c2cos|C|tN(t~)a(t~)×c3sin|C|r+c4cos|C|r (34)

with c1, c2, c3, and c4 constants.

By checking equations (32), (33) and (34) it is possible to conclude that both functions U(t,r) and T(t,r) can be modeled in many ways according to the values of the constants. It is worth to pay attention to constant C, which is directly related to the curvature k of the Friedmann metric. That is, after adjusting the constants ai,bi and ci (i=1,2,3,4) to be zero, we can make R(r) mimic f(r).

For instance, when a3 = 0 in (32) we have R(r) = r (a 4 can be absorbed by the temporal part) which is in accord to f(r) = r for k = 0. In equation (33), which corresponds to C > 0, it is possible to obtain f(r) by taking C = 1. This choice recovers f(r) = sinh(r), that represents k = −1. Last but not least, equation (34) corresponds to the closed universe (k = +1) when C = −1 and c4 = 0. For this case R(r) and f(r) are equal to sin(r).

Therefore, one can recover a Friedmann behavior in the r dependent part, which can match a quasi-isotropic metric. Such behavior is desirable once an isotropic and homogeneous metric has already been proved to be compatible only with an open universe.

Some specific cases related to this are:

4.1±At,r=N(t)u12,±B(t,r)=atu12,Ut,r=a(t)f(r)u1

For this case we have that u1 is an arbitrary constant. The motion equations are given by

2a¨(t)a(t)+a˙(t)2a(t)2-m226+4α3+α4+m26α4u14-m21+2α3+α4u12+4m233+3α3+α4u1+kat2=0 (35)

and

a˙(t)2a(t)2+ka(t)2-m266+4α3+α4+2m233+3α3+α4±m21+2α3+α4u12+2m23α3+α4u13-m26α4u14=0 (36)

4.2±A(t,r)=|Ntu22,±B(t,r)=atu22,Ut,r=-a(t)f(r)u2

In this case u2 is a constant given by:

u2=±3+3α3+α4α3+α4 (37)

The Friedmann equations are:

2a¨(t)a(t)+a˙(t)2a(t)2-m226+4α3+α4+m231+2α3+α4u22+kat2+m26(α4)u24=0 (38)

and

a˙(t)2a(t)2-m266+4α3+α4+2m231+2α3+α4u22+kat2-m26α4u24=0 (39)

This case is more interesting due its resemblance with the open universe case, which allows compatibility with massivegravity. In this case (Gümrükcuöglu et. al. 2011) the fiducial metricis defined as fuv=ηABδμXAδvXB with X=f(t)1+kx2+y2+z2,kf(tx,kfty,kf(t)z). For such a definition we have that ft=a(t)|k|u3, where u3 has the same value of u from equation (47). Then we obtain the following pair of equations of motion:

12a¨(t)a(t)+6a˙(t)2a(t)2+6kat2-3m26+4α3+α4+2m21+2α3+α4u32+m2(α4)u34=0 (40)

and

6a˙(t)2a(t)2-m26+4α3+α4+2m21+2α3+α4u32-m2α4u34+6kat2=0 (41)

In the next two sections we analyze massive gravity by considering off-diagonal elements in the function γ given by (18). In Section B we consider that G(t,r) is zero while F(t,r) is non zero. Section D presents the opposite situation (G(t, r) ̸= 0 and F(t, r) = 0). For both cases it is imposed that At,r=-N(t)a(t)B(t,r) in order to have γ From (17) diagonal.

3.2. Non-diagonal metric: Case 1

As described previously, this case consists of making F(t,r)0 and G(t,r) = 0. After substituting the condition for keeping γ2 (17) diagonal into equation (19) we get the condition to have uT0μ null. The other term necessary to ensure conservation is given by the conservation of T1μ, given by:

(42)

Two different cases are analysed next.

3.2.1. A(t,r) = B(t,r) = 0,F(t,r) = 1/U(t,r) 2 , G(t,r) = 0, α 3 = −2 and α 4 = 3

With this choice there are some restrictions on the functions T(t,r) and U (t,r). By having A(t,r) and B(t,r) null, the only condition is that both are constants or equal. So, it is possible to make U(t,r) equal to a(t)f(r), which would turn this metric closer to the Friedmann one. In any case, the equations of motion are given by

2a¨ta(t)+a˙(t)2a(t)2+ka(t)2-m22=0 (43)

and

a˙(t)2a(t)2+ka(t)2-m26 =0 (44)

Once again, the massive gravity contribution is present as a cosmological constant in the Friedmann equations.

3.2.2. A(t,r) = u 2 N(t) 2 , B(t,r) = u 2 a(t) 2 , U(t,r) = −ua(t)f(r), G(t,r) = 0 and F(t,r) = ua(t)

Here u is a constant to be determined. After some calculation we have the equations of motion:

12a¨(t)a(t)+6a˙(t)2a(t)2-3m26+4α3+α4+5m2u3+3α3+α4+6kat2±m2u21+2α3+α4-m2u3α3+α4=0 (45)

And

6a˙t2a(t)2-m26+4α3+α4+6kat2+m2u3+3α3+α4+m2u21+2α3+α4-m2u3α3+α4=0 (46)

where u is given by:

u=-21+2α3+α4α3+α4±12α3+α4×41+2α3+α42-4α3+α43+3α3+α412 (47)

With a particular choice of α3 = 0 and α4 = −3 we have u = −2 or u = +2/3. Both cases furnish a positive cosmological constant to the mass term, which can also lead to an accelerated expansion.

3.3. Non-diagonal metric: Case 2

Here we follow the same idea of the last subsection, but we have F(t,r) = 0 and G(t,r) 0. Therefore, anisotropy is still kept. Conditions for conservation are obtained after plugging At,r=-N(t)a(t)B(t,r) into (19) and for the T1μ we have:

(48)

We analyse just one case:

3.3.1. A(t,r)=B(t,r)=U(t,r)=F(t,r)=0 and α 4 =−3(α 3 +1)

This case is dictated by the fact that function (18) has only one non-zero element that is off-diagonal. After squaring it, we see that equation (17) is identically zero. This means that K becomes an identity matrix in the Lagrangian. G(t,r) is an arbitrary function with α4 = −3(α3 + 1) being the only necessary condition for conservation of the energy-momentum tensor. From this we obtain the equations of motion,

2a¨(t)a(t)+a˙(t)2a(t)2+ka(t)2-123+α3m2=0 (49)

and

a˙(t)2a(t)2+ka(t)2-163+α3m2=0 (50)

Both equations (49) and (50) show that massive gravity causes a cosmological constant to appear in the Friedmann equations. For this case it is possible to eliminate the gravitation contribution by making α3= −3.

4. CONCLUSION

In this paper we have obtained new cosmological solutions in the massive gravity theory, constructed by means of to set specific forms for the fiducial metric (15). We found different functions A(t,r), B(t,r), U(t,r) and T(t,r) that guarantee the conservation of the energy-momentum tensor. The equations for the energy momentum tensor conservation were presented in the general form for three different cases, one being diagonal in the fiducial metric and two being non-diagonal. The expressions for the energy-momentum conservation are (19), (20), (42) and (48), the equation for the T0μ component being valid for all the cases treated here. In some cases the parameters α were also specified in order to give accelerated solutions. A particular choice for the α parameters also permits to turn an anisotropic solution into an isotropic one, as occurs in the first case, equations (21-22), where the choices α3 = −1 and α4 = 1 lead to the isotropic equations (23-24). We are interested in homogeneous and isotropic solutions in order to respect the cosmological principle. Anisotropic solutions may be interesting in other specific cases, but for the current evolution of the universe we consider homogeneous and isotropic cases.

Several homogeneous and isotropic solutions that correctly describe an accelerated evolution for the universe were also found, all of them reproducing a cosmological constant term proportional to the graviton mass. From the pairs of equations (23-24), (25-26), (35-36), (38-39), (40-41), (43-44), (45-46) and (49-50) it is easy to see that they can all be reduced to an single equation of the form:

a¨a=Cα,um2 (51)

where c(α, u) is a constant function depending on the free parameters α and on the constants u. It is easy to see that the solution to the above equation is:

at=a0exp±C(α,u)mt (52)

where the positive solution is a de Sitter type solution, with the property of providing an accelerated expansion to the universe, here sourced by the graviton mass m.

The free parameters α3 and α4 are the contributions coming from the lagrangians counter terms (2) needed to eliminate the Bouware-Deser ghost present in the theory. From a classical point of view, the constants α are as yet unknown, but maybe their values will follow from a perturbative quantum construction of the theory. Although a theoretical approach in order to obtain the values of such constants is not yet available, observational data can provide constraints. (Pereira et. al. 2016).

We have also verified that when the fiducial metric is close to the physical metric the solutions are absent. As already shown in different works, this is a characteristic of the massive gravity theory, which forbids the fiducial metric to be of the same type as the physical one.

REFERENCES

Ade, P. A. R. et al. 2014, A&A, 571, A16 [ Links ]

Amanullah, R., Lidman, C., Rubin, D., et al. 2010, ApJ, 716, 712 [ Links ]

Arkani-Hamed, N., Georgi, H., & Schwartz, M. D. 2003, AnPhy, 305, 96 [ Links ]

Boulware, D. G. & Deser, S. 1972, PhRvD, 6, 3368 [ Links ]

Creminelli, P., Nicolis, A., Papucci, M., & Trincherini, E. 2005, JHEP, 09, 003 [ Links ]

de Felice, A., Gümrükçuöğlu, A. E., Lin, C., & Mukohyama, S. 2012, PhLB, 717, 295 [ Links ]

de Felice, A., Gümrükcuöglu, A. E., & Mukohyama, S. 2011, JCAP, 11, 030 [ Links ]

de Rham, C., Fasiello, M. & Tolley, A. J., 2014, JMPD, 23, 1443006 [ Links ]

de Rham, C. & Gabadadze, G. 2010, PhRvD, 82, 044020 [ Links ]

de Rham, C., Matas, A. & Tolley, A. J. 2014, CQGra, 31, 165004 [ Links ]

Dubovsky, S. L. 2004, JHEP, 10, 076 [ Links ]

Farooq, O. & Ratra, B. 2013, ApJ, 766, L7 [ Links ]

Fierz, M. & Pauli, W. 1939, PRIAA, 173, 211 [ Links ]

Gannouji, R., Hossain, Md. W., Sami, M. & Saridakis, N. 2013, PhRvD, 87, 3536 [ Links ]

Hassan, S. F. & Rosen, R. A. 2011, JHEP, 07, 009 [ Links ]

______. 2012, PhRvD, 108, 041101 [ Links ]

Hassan, S. F., Rosen, R. A., Schmidt-May, A., & von Strauss, M. 2012, PhLB, 715, 335 [ Links ]

Hinterbichler, K. 2012, RvMP, 84, 671 [ Links ]

Kobayashi, T., Siino, M., Yamaguchi, M., & Yoshida, D. 2014, NuPhS, 246, 76 [ Links ]

Komatsu, E., et al. 2011, ApJS, 192, 18 [ Links ]

Pereira, S. H., Mendonça, E. L., Pinho A. S. S. & Jesus, J. F. 2016, RMxAA, 52, 125 [ Links ]

Riess, A. G. 1998, AJ, 166, 1009 [ Links ]

Sharov, G. S. & Vorontsova, E. G. 2014, JCAP, 10, 57 [ Links ]

Vakili, B. & Khosravi, N. 2012, PhRvD, 85, 083529 [ Links ]

Volkov, M. S. 2012, PhRvD, 86, 061502 [ Links ]

______. 2012, PhRvD, 86, 104022 [ Links ]

Weinberg, S. 1989, RVMP, 61, 1 [ Links ]

SHP is grateful to CNPq Conselho Nacional de Desenvolvimento Científico e Tecnológico, Brazilian research agency, for financial support, grants number 304297/2015-1. E.L.M thanks CNPq 449806/2014-6. APSS thanks CAPES Coordenacao de Aperfeicoamento de Pessoal de Nível Superior, for financial support.

Received: January 26, 2016; Accepted: September 05, 2016

E. L. Mendonça, S. H. Pereira, & A. Pinho S. S.: Faculdade de Engenharia de Guaratinguetá, UNESP Univ. Estadual Paulista “Júlio de Mesquita Filho”, Departamento de Física e Química, Av. Dr. Ariberto Pereira da Cunha 333 Pedregulho, 12516-410-Guaratinguetá, SP, Brazil (shpereira@gmail.com).

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License