Scholarly article on topic 'Equivalent off-diagonal cosmological models and ekpyrotic scenarios in $$f(R)$$ f ( R ) -modified, massive, and einstein gravity'

Equivalent off-diagonal cosmological models and ekpyrotic scenarios in $$f(R)$$ f ( R ) -modified, massive, and einstein gravity Academic research paper on "Physical sciences"

0
0
Share paper
Academic journal
Eur. Phys. J. C
OECD Field of science
Keywords
{""}

Academic research paper on topic "Equivalent off-diagonal cosmological models and ekpyrotic scenarios in $$f(R)$$ f ( R ) -modified, massive, and einstein gravity"

Eu, Phys. J. C (2015) 75:176 THE EUROPEAN , H , CrossMark

DOI 10.1140/epjc/s10052-015-3407-8 PHYSICAL JOURNAL C

Regular Article - Theoretical Physics

Equivalent off-diagonal cosmological models and ekpyrotic scenarios in f (R)-modified, massive, and einstein gravity

Sergiu I. Vacarua

Rector's Department, University "Al. I. Cuza" Iagi, 14 A. Lapugneanu street, Corpus R, UAIC, office 323, Ia§i 700057, Romania

Received: 9 January 2015 / Accepted: 15 April 2015 / Published online: 29 April 2015 © The Author(s) 2015. This article is published with open access at Springerlink.com

Abstract We reinvestigate how generic off-diagonal cosmological solutions depending, in general, on all spacetime coordinates can be constructed in massive and f -modified gravity using the anholonomic frame deformation method. New classes of locally anisotropic and (in-) homogeneous cosmological metrics are constructed with open and closed spatial geometries. By resorting to such solutions, we show that they describe the late time acceleration due to effective cosmological terms induced by nonlinear off-diagonal interactions, possible modifications of the gravitational action and graviton mass. The cosmolog-ical metrics and related Stuckelberg fields are constructed in explicit form up to nonholonomic frame transforms of the Friedmann-Lamaitre-Robertson-Walker (FLRW) coordinates. The solutions include matter, graviton mass, and other effective sources modeling nonlinear gravitational and matter field interactions with polarization of physical constants and deformations of metrics, which may explain dark energy and dark matter effects. However, we argue that it is not always necessary to modify gravity if we consider the effective generalized Einstein equations with nontrivial vacuum and/or non-minimal coupling with matter. Indeed, we state certain conditions when such configurations mimic interesting solutions in general relativity and modifications, for instance, when we can extract the general Painleve-Gullstrand and FLRW metrics. In a more general context, we elaborate on a reconstruction procedure for off-diagonal cosmological solutions which describe cyclic and ekpyrotic universes. Finally, open issues and further perspectives are discussed.

Associated visiting research affiliation of Sergiu I. Vacaru: Theory Division, CERN, CH-1211, Geneva 23, Switzerland ae-mails: sergiu.vacaru@uaic.ro; sergiu.vacaru@gmail.com

Contents

1 Introduction ..........................................1

2 Equivalent modeling of f -modified and massive gravity theories........................................2

3 Generating off-diagonal cosmological solutions ... 4

4 Examples of off-diagonal solutions with solitonic configurations ........................................6

4.1 One soliton solutions............................6

4.2 Three dimensional solitonic anistoropic waves . 7

4.3 Solitonic waves for a nontrivial vertical confor-mal f-factor......................................7

5 Reconstruction mechanism for off-cosmological solutions . 8

6 Conclusions ..........................................9

References........................10

1 Introduction

The "ekpyrotic" and "new ekpyrotic" mechanisms and cyclic models have been elaborated as alternatives to standard big bang inflationary cosmology [2-5]. Another alternative came from the idea that the graviton may have a nontrivial mass as was proposed by Fierz and Pauli [1]. For recent reviews and related f (R) modifications and cosmological models, in general, with non-minimal coupling and dilaton-brane cosmology, local anisotropies and/or effective modeling of massive gravity, see, respectively, [6-12] and [13-18]. The modern version of a ghost-free (bimetric) massive gravity theory was made in a series of papers (recent reviews can be found in [19]): developing generic nonlinear versions of the Fierz-Pauli theory, the so-called vDVZ discontinuity problem was solved via the Vainshtein mechanism [20-22] (avoiding the discontinuity by going beyond the linear theory). There are also more recent approaches based on the DGP model [2326]. During a long time, no solution was found for another problem with ghosts because at nonlinear order in massive

gravity there appears a sixth scalar degree of freedom as a ghost. This problem was considered in a paper by Boulware and Deser [27,28] together with similar issues related to the effective field theory approach.

A considerable amount of work has been made in order to understand the implications and find possible applications of such ghost-free models. The most substantial progress was made when de Rham et al. showed how to eliminate the scalar mode and Hassan and Rosen provided a complete proof for a class of bigravity/bimetric gravity theories; see [29-33]. In such approaches, the second metric describes an effective exotic matter related to massive gravitons and does not suffer from a ghost instability to all orders in a perturbation theory and away from the decoupling limit.

It is important to mention that the first black hole solutions in a nonlinear massive theory were found in the context of high-energy physics and for off-diagonal and/or higher dimensions generalizations [34,35]. A general decoupling of the generalized Einstein equations can be proven following the anholonomic frame deformation methods, AFDM, [3644]. Here we also note that the possibility that the graviton has a nonzero mass ¡1 results not only in fundamental theoretical implications but gives rise to straightforward phenomenolog-ical and cosmological consequences. For instance, a gravitational potential of Yukawa form, ~e—1 r/r, results in the decay of gravitational interactions at scales r > ¡1-1, and such an effect may result in the accelerated expansion of the Universe. In this way, a massive gravity theory, MGT, provides alternatives to dark energy and, via effective polarizations of the fundamental physical constants, may explain certain dark matter effects. We can treat this as a result of generic off-diagonal nonlinear interactions. Recently, various cosmological models derived for ghost-free (modified) massive gravity and bigravity theories have been elaborated and studied intensively [6-10,45-51].

The goal of this work is to construct generic off-diagonal cosmological solutions in MGT and state the conditions when such configurations are modeled equivalently in general relativity (GR). It extends the later variant of [52,53] with new results and more details of the proofs of the results, constructing exact solutions and providing a further analysis of results and speculation on the reconstruction mechanism.

As a first step, we consider off-diagonal deformations of a "prime" cosmological solution taken in a general Painleve-Gullstrand (PG) form, when theFriedman-Lamaitre-Robert-son-Walker (FLRW) can be recast for well-defined geometric conditions. Such constructions are performed in Sect. 2.

At the second step, the "target" metrics will be generated so as to possess one Killing symmetry (or other non-Killing symmetries) and depend on timelike and certain (or all) spacelike coordinates. In general, such off-diagonal solutions come with a local anisotropy and inhomogeneities for effective cosmological constants and polarizations of other

physical constants and coefficients of cosmological metrics, which can be modeled both in MGT and GR. We consider the method of constructing generic off-diagonal cosmological solutions in Sect. 3.

Then (the third step), we shall emphasize and speculate on the importance of off-diagonal nonlinear gravitational interactions for elaborating cosmological scenarios when dark matter and dark energy effects can be explained by anisotropic polarizations of vacuum and/or de Sitter like configurations. We provide examples of off-diagonal solutions with solitonic configurations; see Sect. 4.

A reconstruction mechanism for off-diagonal cosmologi-cal solutions with modified gravity and/or massive graviton effects is elaborated in Sect. 5. Finally, we conclude the paper in Sect. 6.

2 Equivalent modeling of f -modified and massive gravity theories

We shall work with MGTs modeled on a pseudo-Riemannian spacetime V with physical metric g = {g1v} and fiducial metrics. On massive gravity, see for a review [19] and on geometric methods in gravity and constructing exact solutions, see Refs. [36-44]. In addition to well-known approaches with diadic and tetradic (vierbein) variables, we shall work with nonholonomic manifolds where certain classes of frame transforms can be adapted to preserve a chosen splitting of the nonlinear and linear connection structures into some standard components (for instance, defining the Levi-Civita, LC, connection), and we shall deal with distortion tensors which can be fixed to be zero if additional constraints are imposed.

Using different nonholonomic frame variables, the action for our model can be written in two forms:

J &U4^\gap |

1(n/5M Vïg^ÏÏ f (R) + mL ].

f (R) - '—U(g/V, Kap) + mL

(1) (2)

In the above formulas, the physical/geometrical objects with a hat and/or written in boldface form are considered for a conventional 2 + 2 splitting when the two dimensional horizontal, h, and the two dimensional vertical, v, coordinates are labeled (respectively) in the form ua = (xi, ya), or u = (x, y), with indices i, j, k,... = 1, 2 and a, b,... = 3, 4. The scalar curvature R is for the LC-connection V.

We write R for the scalar curvature for an auxiliary (canonical) connection D uniquely determined by two conditions:

1. it is metric compatible, Dg = 0, and

2. the h-and v-torsions are zero (but there are nonzero h — v components of the torsion T completely determined by

g) for a conventional splitting N : TV = hV © vV,

with local coefficients of the so-called nonlinear connection structure, the N-connection structure, denoted in the form N = {Nb}.

General frame transforms can be parameterized in the form ea = Aaa(u)da, where the matrix Aaa is non-degenerate in a finite, or infinite region of V and da = d/du a . Using such Aaa, we can always re-define the geometric and physical object with respect to a class of N-adapted (dual) bases,

= (ei = di - Nfdb, ea = da = d/dya) and = (ej = dxl, eb = dyb + Nbcdyc)%

which are nonholonomic (equivalently, anholonomic) because, in general, relations of type

eaep - ep ea = W^p ey

are satisfied for certain nontrivial anholonomy coefficients WYp(u). The Einstein summation rule on repeated indices will be applied if the contrary is not stated.

The connection D allows us to decouple the field equations in various gravity theories and construct exact solutions in very general forms. Here we note that the distortion relation from the LC-connection,

D = V+Z[T ],

is uniquely determined by a distorting tensor Z completely defined by T and (as a consequence for such models) by (g, N). The main idea of the AFDM [36-44] is to use D as an auxiliary [for the (pseudo-) Riemannian spacetimes] and/or (with nonholonomically induced torsion) connection, which definitely allows us to decouple gravitational field equations for very general conditions with respect to N-adapted frames; see (3). We cannot perform such a decoupling in general form if we work from the very beginning with the data (V, да/), but there are proofs that this is possible for (D, ea; g, N). Having constructed integral varieties (for instance, locally anisotropic and/or inhomogeneous cosmological ones), we can impose additional nonholonomic (non-integrable constraints) when D^=0 ^ V and R ^ R, where R is the scalar curvature of V, and it is possible to extract exact solutions in GR.1

The MGTs with actions of type (1) generalize the so-called modified f ( R) gravity and the ghost-free massive gravity

[29-33]. We shall follow some conventions from [51]. Units will be used where h = c = 1 and the Planck mass MPl is defined via M2Pl = 1/8nG, with 4D Newton constant G. We write 5u4 instead of d4u because N-elongated differentials are used as in (3). A model can be specified by the corresponding constants, dynamical physical equations, and their solutions in the corresponding variables. It will be assumed that ¡1 = const is the mass of the graviton. For LC-configurations, we can fix conditions of type

f(R) - ^UKae) = f (R), or f(R) = f (R), or f(R) = R,

which depend on the type of models we elaborate and on what classes of solutions we want to construct. Such conditions can be very general ones for arbitrary frame transforms and N -connection deformations. We emphasize that it is possible to find solutions in explicit form if we choose the coefficients {Na} and the local frames for D when R = const in such forms that da f (R) = (dRf) x daR = 0, but, in general, da f (R) = 0.

The equations of motion for our nonholonomically modified massive gravity theory can be written

(dRf)R^v - 1 f(R)glxv + il2X^v = M-lTMV, (5)

where MPi is the Planck mass, R|v is the Einstein tensor for a pseudo-Riemannian metric g|v and D, T|v is the standard matter energy-momentum tensor. We note that in the above formulas, for D ^ V, we get R|v ^ Riv with a standard Ricci tensor Riv for V. Such limits give rise to the original modified massive gravity models, but it is important to find some generalized solutions and only at the end to consider additional assumptions on the limits and nonholo-nomic constraints. It should be emphasized that the effective energy-momentum tensor Xlv in (5) is defined by the potential of the graviton U = U2 + a3U3 + a4U4, where a3 and a4 are free parameters. The quantities U2, U3 and U4 are certain polynomials entailing the traces of some other polynomials

of a matrix KJ, = S^ -(Vi-1^)!

by the four Stuckelberg fields as

for a tensor determined

^IXV = д^ф^дуф-п^ь

when = (1, 1, 1, — 1). A series of arguments presented in [51] (geometrically, we can consider corresponding non-holonomic frame transforms and nonholonomic variables) prove that the parameter choice a3 = (a — 1)/3,a4 = (a2 — a + 1)/12 is useful for avoiding potential ghost instabilities. By frame transforms, we can fix

1 There will also be considered left and right upper/lower indices as labels for some geometric/physical objects.

X|v = a 1g|V.

By explicit computations, we can prove that for the configurations (7), de Sitter solutions with an effective cosmo-logical constant are possible, for instance, for an ansatz of PG type,

x h3(r, 0)[dp ■

■m (r, 0)dxi]2

ds2 = U2(r, t)[dr + ejf (r, t)dt] + a2r2(d0

+ sin2 0dp2) — V2(r, t)dt2.

In the above formula, spherical coordinates are used labeled in the form u& = (x1 = r, x2 = 0, y3 = p, y4 = t), where the function f takes non-negative values and the constant a = a/(a + 1) and e = ±1. We can consider bimetric configurations [determined as solutions of the system (5), (6), and (7)] with Stuckelberg fields parameterized in the unitary gauge as = t and = rn1, 02 = rn2,= rn3 where a three dimensional (3D) unit vector is defined as n = (n1 = sin 0 cos p, n2 = sin 0 sin p, n- = cos 0).

First we note that any PG metric of type (8) defines solutions both in GR and in MGT. For instance, we can extract the de Sitter solution, in the absence of matter, and obtain standard cosmological equations with a FLRW metric, for a perfect fluid source

TIV = [p(t) + p(t)] uIÀuv + p(t)g^v,

where uM = (0, 0, 0, -V ) can be reproduced for the effective cosmological constant eff À = I2/a. Secondly, it is also possible to express metrics of type (8) in a familiar cosmological FLRW form (see formulas (23), (24), and (27) in [51]).

Finally, in this section, we note that off-diagonal deformations of such solutions can be constructed for different classes of zero graviton mass and/or massive theories and various f -modifications; see examples in [36-44] and [34,35].

3 Generating off-diagonal cosmological solutions

We now make the crucial assumption that our Universe can be described by inhomogeneous cosmological metrics with Killing symmetry on d3 = dp and, in general, cannot be diag-onalized by coordinate transforms. As a matter of principle, we can consider dependencies on all spacetime coordinates, but this requires more cumbersome computations; see the examples in Refs. [36-44]. Inhomogeneities and anisotropies can be very small, but it is important to find certain classes of general solutions for generic nonlinear systems and impose at the end certain homogeneity and high symmetry conditions by selecting the corresponding subclasses of generation and integration functions. Up to general classes of frame transforms, we can consider the ansatz

ds2 = m(r, 0)gi(r)dr2 + n2(r, 0)g2(r)d02

x [dt + (wi (xk, t) + wi (xk))dxi]2

V4(r, 0, t)fl4(r, 0, t) (10)

The values na are called "polarization" functions, where rn is the so-called "vertical", v, conformal factor.

For metrics (10), we can consider off-diagonal N -coefficients labeled Nf (xk, y4), where (for parameteriza-tions corresponding to this class of ansatz)

N? = m (r, 0) and N4 = wi (xk, t) + Wi (xk).

The data for the "primary" metric are

gi(r) = U2 - h4(Wi)2, g2(r) = a2r2, h3 = a2r2 sin2 0, h4 = fU2 - V2|, = e\[f U2/h4, w2 = 0, ni = 0,

when the coordinate system is fixed in such a way that the values f, U, V in (8) result in a coefficient g1 depending only on r.

Using nonholonomic frame transforms, we can parameterize the energy momentum sources (9) and the effective equation (7) in the form

■y^a _

l(dRf ) 1

M2P lORf)

(Ta + a-1Xp

( mT + a-1)SaB = ( mT + aT)SaB, (12)

for constant values mT := M

M-f(dRÎ)-1

n —1 m

T and aT =

+ rn2(r, 0, p, t){n3(r, 0, t)

Pl (dRf)'

1, with respect to N-adapted frames (3). In general, such sources are not diagonal and may depend on all spacetime coordinates. Our assumption is that we prescribe a distribution with one Killing symmetry in a moment of time and then find the further evolution with respect to certain classes of nonholonomic frames. We emphasize that fixing such N-adapted parameterizations we can decouple the gravitational field equations in MGTs and construct exact solutions in explicit form.

Let us outline in brief the decoupling property of the gravitational and matter field equations in GR and various generalizations/modifications studied in detail in Refs. [36-44]. That anholonomic frame deformation method (AFDM) can be applied for the decoupling, and for constructing solutions of the MGT field equations (5) with any effective source parameterized in the form (12); see details in Refs. [36-44]).

We label the target off-diagonal metrics as g = (gi = nigi, ha = Vaha, Nj) with coefficients determined by the ansatz (10). In these formulas, there is no summation on the repeating indices. We shall use short-hand notations for the partial derivatives: = f •, 32f = f', d3f = f and d4f = f *. Computing the N-adapted coefficients of the Ricci tensors, when Rj = R|, R| = R^, R3k, and R4k

are not trivial, we write (5) as a system of nonlinear partial differential equations (PDEs):

f •• + f " = 2(m Y + aY),

$*h* = 2h3h4(mY + aY), n** + yn* = 0, Rwi - ai = 0, dk m = nk a? + Wk a*,

$ = ln

Vih3h4Í

y := i in^V,

Y 1 h4| / '

h * h* ai = — di$, R = -3-$*.

i 2h3 2h3

In the above formulas, we take the system of coordinates and the polarization functions to be fixed for configurations with g1 = g2 = e^(x ) and nonzero values and h*a.

We can extract solutions for the LC-configurations with zero torsion if the coefficients of metrics are subjected to additional conditions:

}*= ei \nj\h4\, ei in yThl = 0,

di wj = dj wi and n* = 0.

Step by step, the system of nonlinear PDEs (13)-(16) can be integrated in a general form for any œ constrained by a system of linear first order equations (14); see details in [3644]. The explicit solutions are given by quadratic elements:

ds2 = ef(xk )[(dx1)2 + (dx2)2 ] +

.2 ,.,2

h3[dç + (dk n) dxk]2 -

4 ( m y + aY) (3>*)2w2

( m y + aY)$,2

h4[dt + (di A)dxi]2,

for any

$ = $, (di $)* = di $*, wi + wi = di $/$* = di a.

We can construct exact solutions even if such conditions are not satisfied, i.e. the zero torsion conditions are not stated or there are given certain sources in non-explicit form. So, the AFDM can be applied to generate both off-diagonal metrics and nonholonomically induced torsions. There are various physical arguments for what type of generating/ integration functions and sources we have to choose in order to construct realistic scenarios for acceleration of the Universe and observable dark energy/dark matter effects.

Coming back to the properties of general solutions, we note that we can generate new classes of solutions for arbitrary nontrivial sources, m Y + aY = 0, and generating functions, $(xk, t) := e( and nk = dkn(xi). The resulting target metrics are generic off-diagonal and cannot be diago-nalized via coordinate transforms in a finite spacetime region because, in general, the anholonomy coefficients W^p for (3)

are not zero (we can check this by explicit computations). The polarization n-functions from (17) are

m = e^/g!, m = e+/§2, m = $2/4 (mY + aY), n4 = ($*)2/( mY + aY)$2. (18)

We conclude that prescribing any generating functions $(r, 0, t), n(r, 0), rn(r, 0, y, t) and sources mY, aY and then computing A(r,0,t), we can transform any PG (and, similarly, FLRW) metric g = (gi, ha ,Wi, ni) in MGT and/or GR into new classes of generic off-diagonal exact solutions depending on all spacetime coordinates. Such metrics define Einstein manifolds in GR with effective cosmological constants determined by mY + aY. With respect to N-adapted frames (3) the coefficients of the metric encode contributions from massive gravity, determined by aY, and matter fields, included in m Y.

We also note that it is possible to provide an "alternative" treatment of (17) as exact solutions in MGT. In such a case, we have to define and analyze the properties of fiducial Stuckelberg fields and the corresponding bimetric structure resulting in target solutions g = (gi, ha, Na). Let us analyze the primary configurations related to

(p! = ((pi- = a(r)pa~1nL, (p- = a(T)pa~ln-, (p4 = tk-1),

when the corresponding prime PG-metric g is taken in FLRW form,

ds2 = a2(dp2/(1 - Kp2) + p2(d62 + sin2 0dç2)) - dt2.

A corresponding fiducial tensor (6) is computed,

£ av du^du- = —r - a2

dp2 + p 2(d 0 2 + sin2 0 d ç2)

+ 2 Hp dp d t -

- H2p2 I dt

where the coefficients and coordinates are re-defined in the formr ^ p = ar/a(x) and t ^ t = kt, for K = 0, ±1; k is an integration constant; H := d ln a/dt and the local coordinates are parameterized in the form x1 = p, x2 = 0, y3 = y, y4 = T.

In the next step, for a target metric g = gap and frames

— ea eR naR

for ea- =

gij + NfNbgab Nfhab

i -Krh a

ef NS~

The values

gij ~i ~j -iay_

= e-e-nap = efSij = diag[nigi], hah

= e-e-nap = diag[naha], N3 = dk n, N4 = dk A

ea = e^da, we write

are related algebraically to the data (18) resulting in offdiagonal solutions (17). Then, to work out the "target" Stuckelberg fields, wecompute = e^ with e^ , these being inverse to ear, and the fiducial tensor

£ap = (ea<-)(e<)n^v = eore^ap-

bilities if the parameter choice is considered both for diagonal and off-diagonal cosmological solutions; see the details in [60,61]. This special choice also allows us to have a structure

gßV at least in N-adapted frames when the massive

gravity effects can be approximated by effective cosmological constants and exact solutions in MGT, which are also solutions in GR.

If the prime value £^v carries information as regards the two constants k and a, a target tensor is associated to off-diagonal solutions and encodes data about generating and integration functions and via superpositions on possible Killing symmetries, on various integration constants. Similar constructions were elaborated for holonomic and nonholo-nomic configurations in GR; see [54-56].

In the framework of MGT, the two cosmological solutions g and g related by nonholonomic deformations2 are characterized, respectively, by the two invariants

¡«1 = gaPda<adp<t and Iap = gapea<aep<p.

The tensor does not contain singularities because there are no coordinate singularities on the horizon for the PG metrics. It should be emphasized that the symmetry of is not the same as that of £and the singular behavior of Iap depends on the class of generating and integration functions we chose for constructing a target solution g.

In GR, MGTs, and/or Einstein-Finsler gravity theories [36-44], off-cosmological solutions of type (17) were found to generalize various models of Biachi, Kasner, Godel, and other universes. There are known locally anisotropic black holes and wormholes, in general, with solitonic background solutions; see [34,35,57-59]. For instance, Bianchi type anisotropic cosmological metrics are generated if we impose the corresponding Lie algebra symmetries on the metrics. It was emphasized in [51] that "any PG-type solution in general relativity (with a cosmological constant) is also a solution to massive gravity." Such a conclusion can be extended to a large class of generic off-diagonal cosmological solutions generated by effective cosmological constants but it is not true, for instance, if we consider nonholonomic deformations with nonholonomically induced torsion like in metric compatible Finsler theories.

Finally, we note that the analysis of cosmological perturbations around an off-diagonal cosmological background is not trivial because the fiducial and reference metrics do not respect the same symmetries. Nevertheless, fluctuations around de Sitter backgrounds seem to have a decoupling limit, which implies that one can avoid potential ghost insta-

4 Examples of off-diagonal solutions with solitonic configurations

We now consider three examples of off-diagonal cosmological solutions with solitonic modifications in MGT and (with alternative interpretation) GR. Two and three dimensional solitonic waves are typical nonlinear wave configurations which can be used for generating spacetime metrics with Killing, or non-Killing, symmetries and can be characterized by additional parametric dependencies and solitonic symmetries. Moving solitonic configurations can mimic various types of modified gravity dark energy and dark matter effects with a nontrivial gravitational vacuum, polarization of constants and additional nonlinear diagonal and off-diagonal interactions of the gravitational and matter fields.

4.1 One soliton solutions

We shall for simplicity work with a nonlinear radial (solitonic, with a left s-label) generating function,

$ = s $(r, t) = 4 arctan eqa(r-vt)+q0,

and with œ = 1, we construct a metric

ds2 = e+re){dr2 + d 62) + $

4 ( m y + a y)

h3(r,6)d y2

(dt s $)2 ~ 2

h4(r, t)[dt + (dr A)dr] . (20)

2 These solutions involve not only frame transforms but also deformation of the linear connection structure when at the end there are imposed additional constraints for a zero torsion.

( m y + ay) s $ 2

In this metric, for simplicity, we fixed n(r,0) = 0 and consider that A(r, t) is defined as a solution of s $'/s $* = dr A and ha are given by PG-data (11). The generating function (19), where a2 = (1 - v2)-1 for constants q, q0, v, is just a 1-soliton solution of the sine-Gordon equation

s $**- s + sin s $ = 0.

For any class of small polarizations with na ~ 1), we can assume that the source ( my + ay) is polarized by s$-2 when h3 ~ h3 and h4 ~ h4(s$*)2/ s$-4 with an offdiagonal term dr A, resulting in a stationary solitonic universe. If we consider (dRf )-1 = s $-2 in (12), we can model f -interactions of type (1) via off-diagonal interactions and "gravitational polarizations".

In the absence of matter, m Y = 0, the off-diagonal cosmology is completely determined by aY when s $ transforms /1 into an anisotropically polarized/variable mass of solitonic waves. Such configurations can be modeled if m Y ^ aY. If m Y ^ aY, we generate cosmological models determined by a distribution off matter fields when contributions from massive gravity come with a small anisotropic polarization. For a class of nonholonomic constraints on $ and f (which may be not of solitonic type), when the solutions (20) are of type (10) with na ~ 1 and ni ,wi ~ 0, we approximate PG-metrics of type (8).

Hence, by an appropriate choice of generating functions and sources, we can model equivalently modified gravity effects, massive gravity contributions or matter field configurations in GR and MGT interactions. For well-defined conditions, such configurations can be studied in the framework of some classes of off-diagonal solutions in Einstein gravity with effective cosmological constants.

4.2 Three dimensional solitonic anistoropic waves

More sophisticated nonlinear gravitational and matter field interactions can be modeled both in MGTs and GR if we consider more general classes of solutions of effective Einstein equations.

For instance, off-diagonal solitonic metrics with one Killing symmetry can be generated, for instance, if we take instead of (19) a generating function s$(r, 0, t) which is a solution of the Kadomtsev-Petivashvili, KdP, equations [62],

± s $" + (s $* + s $ s $' + e s$•••)• = 0,

when solutions induce a certain anisotropy on 0.3 In the dis-persionless limit, e ^ 0, we can assume that the solutions are independent on 0 and determined by Burgers' equation,

s $* + s $ s = 0.

Such solutions can be parameterized and treated similarly to (20) but with, in general, a nontrivial term (d0 A)d0 after h4, where s$' / s$* = A' and s/ s$* = A'. Similar metrics were constructed for Dirac spinor waves and solitons in anisotropic Taub-NUT spaces and in five dimensional brane gravity which can be encoded and classified by the corresponding solitonic hierarchies and geometric invariants; see [57-59]. Here we proved that AFDM can be extended to generate inhomogeneous off-diagonal cosmological solitonic solutions in various MGTs.

3 In a similar form, we can construct various types of vacuum gravitational 2D and 3D configurations characterized by solitonic hierarchies

and related bi-Hamilton structures, for instance, of KdP equations with possible mixtures with solutions for 2D and 3D sine-Gordon equations etc; see for the details Refs. [57-59].

4.3 Solitonic waves for a nontrivial vertical conformal f-factor

The cosmological solutions we look for also come with three dimensional solitons but for a f-factor as in (10). For instance, we consider solitons of KdP type, when m = m(r, y, t), when x1 = r, x2 = 0, y3 = y, y4 = t, for

± à™ + (dtœ + œœ' + aœ'")' = 0.

In the dispersionless limit, e ^ 0, we can take the solutions to be independent on the angle y and to be determined by Burgers' equation,

m* + rn rn' = 0.

The conditions (14) impose the additional constraint

e1m = m' + w\(r, 0, y)m* + H\(r, 0)of = 0.

In the system of coordinates when m' = 0, we can fix w2 = 0 and h2 = 0. For any arbitrary generating function with LC-configuration, $(r,0,t), we construct exact solutions,

2 _ 0f(r,e)

ds2 = e

(dr2 + d 0 2) +

4 ( m y + a X)

h3(r,0)d y2

(dt $)2œ

( m X + ay) $2

h4(r, t)[dt + (dr A)dr]2, (22)

which are generic off-diagonal and depend on all spacetime coordinates. Such stationary cosmological solutions come with polarizations on the two angles 0 and y. Nevertheless, the character of the anisotropies is different for metrics of type (20) and (22). In the third class of metrics, we obtain a Killing symmetry on only in the limit rn ^ 1, but in the first two such a symmetry exists generically. For (22), the value $ is not necessarily a solitonic one which can be used for additional off-diagonal modifications of solutions and various types of polarizations.

We can provide a physical interpretation of (22) which is similar to (20) if the generating and integration functions are chosen to satisfy the conditions na ~ 1 and Hi, wt ~ 0. We approximate the PG-metrics of type (8).

In a particular case, we can use a conformal f-factor which is a 1-solitonic one, i.e.

& ^ m(r, t) = 4 arctan eqa(r-vt)+q0,

where a2 = (1 - v2)-1 and constants q, q0, v, define a 1-soliton solution of the sine-Gordon equation,

of* - rn" + sin & = 0.

Such a soliton propagates in time along the radial coordinate.

We conclude that solitonic waves may mimic both particle type configurations as dark matter and encode certain hidden dark energy and off-diagonal gravitational and matter field interactions.

5 Reconstruction mechanism for off-cosmological

solutions

We now consider a reconstruction mechanism with distinguished off-diagonal cosmological effects [36-44] by generalizing some methods elaborated for f (R) gravity in [610]. The main idea was to present the MGT actions (with zero or nonzero gravitational mass) as sums of actions in GR and certain effective ideal fluid contributions with parameters defined by nontrivial ¡1 and f -deformations. The reconstruction method was developed for such formulations as lead to a cosmology with cyclic evolution. Then it was proven that the ekpyrotic scenario may also be realized for MGTs and that it is possible to reconstruct models of f (R) gravity which induce a little rip cosmology.

In our approach, we work with cosmological generic off-diagonal metrics and generalized connections. We can always choose such generating functions and parameters of solutions so that off-diagonal contributions are small, and the torsion is constrained to be zero, beginning at a certain fixed moment of time. Nevertheless, we cannot neglect for such classes of solutions possible nonlinear effects determined by generating functions and effective sources resulting in diag-onalized modifications. MGTs were elaborated as realistic alternatives for a unified description of inflation with dark energy when cosmological scenarios encode information on various massive and zero-mass gravitational modes. A crucial question is if and how such constructions have to be elaborated for generic off-diagonal cosmological spaces. This is not only a geometric problem for generalizing the reconstruction formalism for inhomogeneous and locally anisotropic cosmological theories. It is connected to a very important question of equivalent modeling of cosmological scenarios for different MGTs in the framework of GR with nonholo-nomic and off-diagonal nonlinear interactions.

Any cosmological solution in massive, MGT, and/or GR parameterized in the form (10) [in particular, as (20) and (22)] can be encoded into an effective functional (4) when

where t is the cosmological time. In order to extract a mono-tonically expanding and periodic cosmological scenario, we parameterize ln |a(t )| = H0t + a(t) for the periodic function a(t + t) = 1a cos(2nt/t), where 0 < 1a < H0. Our goal is to prove that such a behavior is encoded into off-diagonal solutions of type (20)-(22).

We write FLRW like equations with respect to N-adapted (moving) frames (3) for a generalized Hubble function H,

3 H2 = 8np and 3 H2 + 2e4 H = -8n p.

Using variables with da f (R)R=const = 0, we can consider a function H (t) when e4 H = dtH = H *. It should be noted that the approximation e4 ^ dt is considered "at the end" (after a class of off-diagonal solutions was found for a necessary type connection, which allowed one to decouple the equations) and the generating functions and effective sources are constrained to depend only on t. The energy-density and pressure of an effective perfect fluid are computed thus:

p = (8n)—1

2 f (R) + 3H e4(3 Rf )) - 3e4 H

= (8n)—1 dR ln/fj -3 H *

= (8n)~

dRlnj | ¡4-U + f (R)|—3H*

p = — (8n)—1 (dRf )-1 Q f ( R) + 2H e4(3Rf ) + e4e4(dRf )\ +e4 H

= (8n)—1

= —(8n)-

dR W | f \ + H

dR ln,/1 ^ U + f ( R)\ + H *

We emphasize that, in general, such values are defined with respect to nonholonomic frames. The effect of nonlinear deformations encoding physical data (¡ï, U, f, R) is preserved also in diagonal contributions for e4 ^ dt.

The equation of state, EoS, parameter for the effective dark fluid encoding MGTs parameters is defined by

f — ¡4u = f ( R),

R\D^V = R.

This allows us to work as in MGT; the conditions da f (R) = 0 if R = const simplify the computations substantially. Nevertheless, the contributions of the parameter ¡ï and effective potential U are included as functional dependencies of f and R. For D ^ V, we get nonlinear modifications both in diagonal and off-diagonal terms of the cosmological metrics.

The starting point of our approach is to consider a prime flat FLRW like metric

ds2 = a2(t )[(dx1)2 + (dx2 )2 + (dy3)2] — dt2,

f + 2H*d¡4U + f (R) + 2H*dRf

f — 6 H *d Rf ¡4- U + f ( R) — 6 H * dRf

The corresponding EoS is

p = —p — (2n)—1H *.

U(t) is computed, for simplicity, for a configuration of "target" Stuckelberg fields = e1 , and the solution found is, finally, modeled by generating functions with dependencies on t. The components of are computed with respect to N-adapted frames.

Taking a generating Hubble parameter H (t) = H0t + H1 sin mt, for m = 2n/x, we can recover the modified action for the oscillations of an off-diagonal (massive) universe (see similar details in [6-10]),

— 1 + sin Mt/3mH1 cos2 Mt ^ 1.

f (R(t)) = 6mHw dt[M sin Mt — 4 cos Mt (H0 + H1 sin Mt)]

x exp[H0t +--sin Mt].

Here we note that we cannot analytically make an inversion to find an explicit form R, or any nonholonomic deformation to R with respect to general N-adapted frames. Nevertheless, we can prescribe any values of constants H0 and H1 and of m and compute the effective dark energy and dark matter oscillating cosmology effects for any off-diagonal solution in massive gravity and/or effective MGT, GR. To extract the contributions of ¡1 we can fix, for instance, f (R) = R = R, and using (1) and (2) we can relate f (R(t)) and the respective constants to certain observable data in cosmology.

The MGT theories studied in this work encode, for respective nonholonomic constraints, the ekpyrotic scenario, which can be modeled similarly to f (R) gravity. A scalar field is introduced into the usual ekpyrotic models in order to reproduce a cyclic universe, and such a property exists if we consider off-diagonal solutions with massive gravity terms and/or f-modifications. The main idea is to develop the reconstruction techniques for the scalar-tensor theory using the AFDM, with nonholonomic off-diagonal metric and linear connection deformations. Working with general classes of off-diagonal cosmological solutions, the problem is to state the conditions for generating and integration functions when a corresponding ekpyrotic scenario will "survive" for certain constraints and in the respective limits.

Let us consider a prime configuration with energy-density for pressureless matter pm; for radiation and anisotropies we take, respectively, pr and pa. k is the spatial curvature of the universe and for the target effective energy-density p we have (23). A FLRW model can be described by

3 H2 =

Pm Pr Pa K

—T +--J +--T--Ô + P

a3 a4 a6 a2

We generate an off-diagonal/massive gravity cosmological cyclic scenario containing a contracting phase by solving the initial problems if w > 1; see (24). A homogeneous and isotropic spatially flat universe is obtained when the scale factor tends to zero and the effective f -terms (massive gravity and off-diagonal contributions) dominate over the rest. In such cases, the results are similar to those in the inflationary scenario. For recovering (25), the ekpyrotic scenario occurs and mimics the observable universe for t ~ k/2m in the effective EoS parameter

This allows us to conclude that in massive gravity and/or using off-diagonal interactions in GR, cyclic universes can be reconstructed in such forms that the initial, flatness, and/or horizon problems can be solved. We can compute possible locally anisotropic and inhomogeneous small contributions for self-consistent models with nonholonomic frames.

In the diversity of off-diagonal cosmological solutions which can be constructed using the above presented methods, there are cyclic ones with singularities of the type of a big bang/big crunch behavior. This is still a largely unexplored area both as regards the geometric methods of constructing exact solutions of PDEs and recovering procedures for certain "preferred" fundamental physical objects compatible with the experimental data. Choosing necessary types of generating and integration functions, we can avoid singularities and elaborate models with a smooth transition. Using the possibility to generate nonholonomically constrained f-models with equivalence to certain classes of solutions in massive gravity and/or off-diagonal configurations in GR, we can study in this context, following methods in [6-10,36-44], big rip and/or little rip cosmology models. In this case the phantom energy-density is modeled by off-diagonal interactions. We note that such nonsingular models for the dark energy were proposed as alternatives to ACDM cosmologies, [63-65]. Using the corresponding classes of the generating functions, we can reproduce the scenarios with a phantom scalar modeling a little rip in the framework of AFDM and nonholonomic MGTs. We omit such considerations in this article (see a summary of such a construction and further developments in [52,53]).

6 Conclusions

In this paper we deal with new classes of cosmological offdiagonal solutions in massive gravity with flat, open, and closed spatial geometries. These solutions can be systematically constructed for various types of modified gravity theories, MGTs, and in general relativity, GR. We applied an advanced geometric technique for decoupling the field equations and constructing exact solutions in massive and zero mass f (R) gravity, theories with nontrivial torsion, nonholo-nomic constraints to GR, and possible extensions on (co) tangent Lorentz bundles. The so-called anholonomic frame deformation method, AFDM, was elaborated during the last 15 years in our work [36-44,57-59], where a number of examples of off-diagonal solutions and new applications in gravity and modern cosmology were considered.

A very important property of such generalized classes of cosmological solutions is that they depend, in general, on all spacetime coordinates via generating and integration func-

tions and constants. They describe certain models of inho-mogeneous and locally anisotropic cosmology with a physical meaning and possible physical implications that are less clear [13-18]. After some classes of solutions were constructed in the most general form, we can impose at the end additional nonholonomic constraints, cosmological approximations, extract configurations with a prescribed spacetime symmetry, and/or a dependence on certain mass parameters; and we can consider asymptotic conditions, etc. Thus, our solutions can be used for elaborating homogeneous and isotropic cosmological models with an arbitrary spatial curvature, to study generalized Killing and non-Killing symmetries, with possible nonholonomically deformed (super-) symmetries [34,35,54-56] and to study the "non-spherical" collapse models of the formation of cosmic structure such as stars and galaxies (see also [51]).

Another aspect of the AFDM is that if we work only with cosmological scenarios for diagonalizable metrics, there are possibilities to discriminate the massive gravity theory from the f -gravity and/or GR. For diagonal metrics depending on a time like coordinate, we can formulate mathematical cosmology problems for certain nonlinear systems of ordinary differential equations with general solutions depending on integration constants. Identifying such a constant, for instance, with a graviton mass parameter, we do not have many possibilities to mimic a number of similar effects in GR for different MGTs. Following only such a "diagonal" approach, we positively have to modify the GR theory in order to explain the observational data in modern cosmology and elaborate realistic quantum models of massive gravity.

We now discuss a new and important feature of the off-diagonal anisotropic configurations, which allows us to model cosmic accelerations and massive gravity and/or dark energy and dark matter effects as certain effective Einstein spaces. Having integrated such a system of nonlinear partial differential equations, PDEs, for a large class of such solutions, we can ask such questions as: Maybe we do not need to modify radically the GR theory but only to extend the constructions to off-diagonal solutions and nonholonomic systems and try to apply this in modern cosmology? Could we explain observational data in modern cosmology via nonlinear diagonal and/or off-diagonal interactions with nonminimal coupling for matter and/or different phases of massive and zero mass gravity? This is a quite complicated theoretical and experimental problem and the main goal of this and our recent papers cited in Refs. [36-44,52,53] was to analyze such constructions from the viewpoint of massive gravity theory when off-diagonal effects can be alternatively explained by other types of gravity theories.

The reconstruction procedure for cosmological models with non-minimally coupled scalar fields evolving on a flat FLRW background and in different MGTs was studied in

[6-12]. In this work, we elaborated a reconstruction method for the massive gravity theory which admits an effective offdiagonal interpretation in GR and f -modified gravity with a cyclic and an ekpyrotic universe solution. We concluded that the expansion can be around the GR action even if we admit a nontrivial effective torsion. For zero torsion constraints, it is possible to construct off-diagonal cosmological models keeping the approach in the framework of the GR theory. We further investigated how our results indicate that theories with massive gravitons, with possible f -modified terms and off-diagonal interactions, may lead to more complicated scenarios of cyclic universes. Following such nonlinear (offdiagonal) approaches, the ekpyrotic (little rip) scenario can be realized without the need to introduce additional fields (or modifying gravity) but only in terms of massive gravity or GR. Another interesting construction can be related to the reconstruction of scenarios of f (R) and massive gravity theories leading to little rip universes both in locally anisotropic and isotropic variants. Finally, we note that the dark energy for little rip models present an example of non-singular phantom cosmology.

Acknowledgments The work is partially supported by the Program IDEI, PN-II-ID-PCE-2011-3-0256 and performed for a corresponding associate visiting program at CERN. The author thanks S. Capozziello, E. Guendelman, E. Elizalde, N. Mavromatos, M. Sami, D. Singleton, P. Stavrinos, and S. Rajpoot for important discussions, critical remarks, collaboration, and substantial support.

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP3.

References

1. M. Fierz, W. Pauli, Proc. R. Soc. Lond. A 173, 211 (1939)

2. J. Khoury, B.A. Ovrut, P.J. Steinhardt, N. Turok, Phys. Rev. D 64, 123522 (2001)

3. J.K. Erickson, D.H. Wesley, P.J. Steinhardt, N. Turok, Phys. Rev. D 69,063514 (2004)

4. E.I. Buchbinder, J. Khoury, B.A. Ovrut, Phys. Rev. D 76, 123503 (2007)

5. P.J. Steinhardt, N. Turok, Science 296, 1436 (2002)

6. K. Hinterbichler, Rev. Mod. Phys. 84, 671 (2012)

7. S. Nojiri, S. Odintsov, N. Shirai, JCAP 1305, 20 (2013)

8. S. Nojiri, S. Odintsov, Phys. Lett. B 716, 377 (2012)

9. S. Nojiri, S. Odintsov, D. Saez-Gomez, AIP Conf. Proc. 1458, 207 (2011). arXiv:1108.0767

10. Y.-F. Cai, E.N. Saridakis, Phys. Rev. D 90, 063528 (2014)

11. AYu. Kamenshchik, A. Tronconi, G. Ventruri, SYu. Vernov, Phys. Rev. D 87, 063503 (2013)

12. J. Ellis, N.E. Mavromatos, D.V. Nanopoulos, Phys. Lett. B 732, 380 (2014)

13. A.P. Kouretsis, Eur. Phys. J. C 74, 2079 (2014)

14. S. Basilakos, A.P. Kourtesis, E.N. Saridakis, P. Stavrinos, Phys. Rev. D 88,123510(2013)

15. P. Stavrinos, S. Vacaru, Class. Q. Grav. 30, 055012 (2013)

16. E. Guendelman, H. Nishino, S. Rajpoot, Phys. Lett. B 732, 156 (2014)

17. E.I. Guendelman, D. Singleton, JHEP 1405, 096 (2014)

18. E. Guendelman, A. Kaganovich, E. Nissimov, S. Pacheva, Phys. Lett. B 718, 1099 (2013)

19. C. de. Rham, Liv. Rev. Rel. 17, 7 (2014)

20. H. van Dam, J.J.G. Veltman, Nucl. Phys B 22, 397 (1970)

21. V. Zakharov, JETP Lett. (Sov. Phys.) 12, 312 (1970)

22. A.I. Vainshtein, Phys. Lett. B 39, 393 (1972)

23. G.R. Dvali, G. Gabadadze, M. Porrati, Phys. Lett. B 485, 2008 (2000)

24. C. Deffayet, G.R. Dvali, G. Babadaze, Phys. Rev. D 65, 044023 (2002)

25. M.A. Luty, M. Porrati, R. Rattazzi, JHEP 0309, 029 (2003)

26. A. Nicolis, R. Rattazzi, JHEP 0406, 059 (2004)

27. D.G. Boulware, S. Deser, Phys. Rev. D 6, 3368 (1972)

28. N. Arkani-Hamed, H. Gergi, M.D. Schwarz, Ann. Phys. 305, 96 (2003)

29. C. de Rham, Phys. Lett. B 681, 89 (2009)

30. C. de Rham, G. Gabadadze, Phys. Rev. D 82, 044020 (2010)

31. C. de Rham, G. Gabadadze, A.J. Tolley, Phys. Rev. Lett. 106, 231101 (2011)

32. S. Hassan, R.A. Rosen, Phys. Rev. Lett. 108, 041101 (2012)

33. S. Hassan, R.A. Rosen, JHEP 1204, 2012 (2012)

34. A. Salam, J.A. Stradhee, Phys. Rev. D 16, 2688 (1977)

35. T. Gheorghiu, O. Vacaru, S. Vacaru, Eur. Phys. J. C 74,3152 (2014)

36. S. Vacaru, Eur. Phys. J. C 74, 3132 (2014)

37. S. Vacaru, Int. J. Mod. Phys. D 21, 1250072 (2012)

38. S. Vacaru, Gen. Relat. Grav. 44, 1015 (2012)

39. S. Vacaru, Europhys. Lett. (EPL) 96, 5001 (2011)

40. S. Vacaru, Class. Q. Grav. 28, 215991 (2011)

41. S. Vacaru, Int. J. Geom. Meth. Mod. Phys. 8, 9 (2011)

42. S. Vacaru, J. Math. Phys. 49, 043504 (2008)

43. S. Vacaru, D. Singleton, Class. Q. Grav. 19, 3583-3601 (2002)

44. S. Vacaru, F.C. Popa, Class. Q. Grav. 18, 4921-4938 (2001)

45. G. D'Amico, C. de Rham, S. Dubovsky, G. Gabadadze, D. Pirt-skhalava, A.J. Tolley, Phys. Rev. D 84, 124046 (2011)

46. A.E. Gumrukcuoglu, C. Lin, S. Mukoyama, JCAP 1203, 006 (2012)

47. M. von Strauss, A. Schmidt-May, J. Enander, E. Mortsell, S.F. Hassan, JCAP 1203, 042 (2012)

48. M.S. Volkov, JHEP 1201, 035 (2012)

49. S. Koyama, G. Niz, G. Tasinato, JHEP 1112, 065 (2011)

50. P. Gratia, W. Hu, M. Wyman, Phys. Rev. D 86, 061504 (2012)

51. T. Kobayashi, S. Siino, M. Yamaguchi, D. Yoshida, Phys. Rev. D 86, 061505 (2012)

52. S. Vacaru, Off-diagonal ekpyrotic scenarios and equivalence of modified, massive and/or Einstein gravity. arXiv:1304.1080

53. E. Elizalde, S. Vacaru, Einstein spaces modeling non-minimal modified gravity. arXiv:1310.6868

54. R. Geroch, J. Math. Phys. 12, 918 (1971)

55. R. Geroch, J. Math. Phys. 12, 394 (1972)

56. S. Vacaru, Int. J. Geom. Meth. Mod. Phys. 4, 1285 (2008)

57. S. Vacaru, Acta Appl. Math. 110, 73 (2010)

58. S. Vacaru, J. Math. Phys. 50, 073503 (2009)

59. S. Vacaru, D. Singleton, Class. Q. Grav. 19, 2793 (2002)

60. L. Berezhiani, G. Chkareuli, C. de Rham, G. Gababadze, A.J. Tolley, Phys. Rev. D 85, 044024 (2012)

61. K. Koyama, G. Niz, G. Tasinato, JHEP 1112, 065 (2011)

62. B.B. Kadomtsev, V.I. Petviashvili, Doklady Akademii Nauk SSS 192, 753 (1970) [Sov. Phys. Dokl. 15, 539 (1970) (English translation)]

63. P.H. Frampton, K.J. Ludwick, R.J. Schrerrer, Phys. Rev. D 84, 063003 (2011)

64. P.H. Frampton, K.J. Ludwick, S. Nojiri, S.D. Odintsov, R.J. Schrerrer, Phys. Lett. B 708, 204-2011 (2012)

65. E. Elizalde, D. Saez-Gomez, Phys Rev. D 80, 044030 (2009)