SciELO - Scientific Electronic Library Online

vol.40 issue3Hot electron transport properties in characteristics of wurtzite GaN MESFETs using a five-valley modelTorsional vibration of carbon nanotubes under initial compression stress author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Brazilian Journal of Physics

Print version ISSN 0103-9733

Braz. J. Phys. vol.40 no.3 São Paulo Sept. 2010 


Effective cosmology a la Brans-Dicke with a non-minimally coupling massive inflaton field interacting with minimally coupling massless field



El-Nabulsi Ahmad Rami*

Department of Nuclear and Energy Engineering, Cheju National University, Ara-dong 1, Jeju 690-756, South Korea




We discuss an effective cosmology a la Brans-Dicke with two interacting scalar fields: a non-minimally coupling massive inflaton Higgs-like scalar field Φ interacting with a minimally coupling massless scalar field χ. Several features are observed and discussed in some details.

Keywords: Interacting scalar fields, non-minimal and minimal couplings, quintessence




The nature of the dark matter and dark energy component as one of the basic ingredients responsible for the accelerated expansion of the universe represents nowadays as one of the most profound and difficult problem in modern cosmology and theoretical physics.1,2,3,4,5,6 One simplest method to deal with this problem is by means of the Einstein's lambda or positive cosmological constant. This later, acts in fact, as an isotropic and homogenous de-Sitter inflationary phase-source in the sense that the cosmological equation of state is pΛ = -ρΛ = Λ and consequently, avoiding problems with fine-tunning of initial conditions in the early universe. Other options and alternatives exist in the literature: decaying cosmological constant (Λ), cold dark matter, dynamical Λ in the form of scalar field with self interacting exponential and inverse power-law potentials: quintessence and tracker field, k-essence, viscous fluid, Chaplygin gas, Brans-Dicke (BD) pressureless solutions, self-interacting BD cosmology with positive power-law potential, etc.7-19 Most of these theories are accompanied with problems and difficulties. For example, the theoretical predicted value of the cosmological constant surpasses by 120 order of magnitudes the observational value. As for the well-known quintessence theory with scalar field Φ and equation of state pΦ = ωρΦ acting as fluctuating dark energy, fine tuning parameters and several constraints are required. A lot of works have been devoted to the investigation of a realistic cosmological model with the non-minimal coupling between gravity and inflaton scalar field and to their connection with inflationary cosmology and phase transitions (spontaneous symmetry breaking (SSB)) in the early Universe.

For minimally coupled scalar field theories, ω > -1, but recent observations and phenomenological non-minimally coupled theories have showed that ω < -1 are also welcome.20,21,22,23 Minimally coupled theories with ω < -1 were proved to possess spatial gradient instabilities that would be ruled out by CMB observations.24 It is noteworthy that if the kinetic term has the wrong term, it is possible to achieve a minimally coupled theory with ω < -1. Due to such growing interests of scalar-tensor theories of gravitation, it is required to study the theory in some more detail. It is noteworthy that the inflaton scalar field Φ satisfies in a curved spacetime the Klein-Gordon equation and one need to set up, in general, a non-minimal coupling term that accommodates both the scalar field and the Ricci curvature of spacetime.

On the other hand, the non-minimal coupling represents an interesting deviation from the equivalence principle because it respects the geometrical nature of gravity and gauge symmetry of electromagnetism. It is in reality introduced by renormalization even if it is absent at the classical level. It is also required in classical general relativity by the Einstein equivalence principle.22,23,25,26,27 It is remarkable that scalar fields have been often used to explain the dark energy problem, and referred as quintessence scalar fields. The main parameter for these appealing cosmological models arises from the potential of the scalar field V(Φ), and one can note for example that scaling potentials have been highly investigated. It is also recognized that scalar fields can account either for dark matter or for dark energy. In both cases, the behaviour of the scalar field is determined by its scalar potential. If one desires to use the same scalar field to unify both problems, the difficulty is to determine which scalar potential can provide matter behaviour at local cosmological scales and satisfy the observational constraints. In reality, it is possible to build up a cosmological model which could advantageously restore a model containing two dark components. This fact may have interesting consequences in different aspects of gravity theories including supergravities.27,28

However, it has been recognized that reheating process and SSB may be performed through two scenarios using the Hartree mean-field approximation: the first one is that the inflaton scalar field Φ is converted through a self-interacting quartic potential λΦ4 (λ 10-12 from COBE observations) into many particles and the second one, is the presence of a new (real) source fieldχcreated through some coupling with the inflaton as g2Φ2 χ2, where g is a coupling constant acting as a free-parameter.29,30,31,32 The non-minimally coupled of the scalar field Φ to R for massive inflaton field were done and was proved to be successful to describe a geometric reheating model with large negative ξ. Surprisingly, the field χ was found to be amplified and that GUT scale gauge boson with mass 1016 GeV can be produced in the preheating stage, a scenario not successfully described with a massive-inflation quadratic potential alone with |ξ| 10-3.33 Besides, it was recently argued that large |ξ|( 104 ) allows using the Standard Model Higgs field for inflation. At the same time the case |ξ| 10-3 makes inflaton with quartic potential compatible with the CMB observations. For the coupling constant |ξ| 10-3 very small amount of tensor modes is expected for the quartic inflation. Interestingly, the non-minimally coupled inflation with quadratic potential and |ξ| 10-3 is no longer compatible with observations, generating too small spectral index.34

In this paper, we consider a more general dynamics of a non-minimally coupling inflaton scalar field Φ with the spacetime curvature R interacting with a minimally coupling massless and scalar field χ. Despite the fact that static massless scalar fields play a crucial role in Schwarzschild-like geometry, it was argued recently that they are connected to many cosmological solutions based on some-kind of duality transformations.35 They also play a crucial role in quantum fields in curved spacetime, in particular in the background of an Einstein-like universe36 and dilatonic electrodynamics in flat spacetime.37 On the other hand, we argue that massless static scalar field could contribute to the dark energy problem by supplying a substantial amount of energy density to the early universe, yet the astronomical data are not precise enough to pinpoint whether the dark energy is truly static or dynamical. Additionally, the term gΦ2 χ2 may have an electromagnetic origin38 as in the early universe, we strongly believe that matter was highly ionized and coupled to scalar field. Latter during cooling as a result of expansion of ions combined to form neutral matter. Furthermore, massless chaotic potentials play a leading role in inflationary scenarios and physics of the early universe and supergravity39,40.

The paper is organized as follows: In Sec.II, we introduce the Lagrangian of the theory and we analyze the cosmological FRW solutions around the background fields Φ and χ. In the same section, we derive the pressure, the density and the equation of state parameter and we analyze their dynamical evolutions in time. In Sec. III, we discuss briefly the role of the fields Φ and χ in the early evolution epoch of the universe. Sec. IV is devoted for conclusions and perspectives.



The Lagrangian of the theory is written in one of the following forms a la Brans-Dicke:

where κ2 = 8πG, G is the gravitational constant, g is the gauge coupling, ξΦ is a coupling constant which corresponds for the scalar field Φ, mΦ is the dynamical mass of the field Φ, hΦ and hc are free positive parameters, ΛΦ is the corresponding cosmological constant. Obviously, the effective mass of the scalar field Φ is = ξΦ R+g2 χ2 which in turns tends to zero for R = 0(flat spacetime) as well as in the absence of a the massless field χ. Amazingly, in the majority of inflationary cosmological scenarios, the coupling constant ξΦ is not a free parameter that could be tuned randomly, but its value is fixed by the gravitational theory and of the scalar field adopted. Moreover, the theoretical consistence of many inflationary scenarios is deeply affected by their proper values. Some phenomenological scenarios turn out to be theoretically contradictory, while others are viable according to the correct use of non-minimal coupling constants. Once their proper values are determined, one does not possess any longer the freedom to regulate their values and the fine-tuning problems that may outbreak the inflationary scenario re-emerge. It is noteworthy that in quantum field theory in curved spacetime, non-minimal coupling is to be expected when the spacetime curvature is large.

The field equations are obtained by varying the action with respect to the metric and to the scalar fields Φ and χ:



To analyze the cosmological solutions around the background fields Φ and χ, we adopt the Friedmann-Robertson-Walker (FRW) flat spacetime with metric d2 = -dt2 + a2 (t)d2 strongly favored by cosmological observations and where from theoretical point of views, the inflaton is distributed homogeneously. Here a(t) is the scale factor of the universe. In order to explore the time-variation of the equation of state parameter, we write the Higgs interacting-fields potential of the field Φ as V(Φ) = (1\2)µ2Φ2 +(h/4)Φ4 where µ2 = Rξ+g2 χ2 with h hΦ hχ > 0 for simplicity. In fact making use of the Bianchi identity, the total energy-momentum tensor is not a conserved quantity. However, since we have assumed that baryon-matter and radiation are considered as non-interacting fields from the decoupling age up to the present time, we presume through this work that the covariant divergence of the energy-momentum tensors of the baryonic matter and radiation vanishes. The resulting Klein-Gordon equation for the scalar fields Φ and χ, and the field equations are:

Equations (7)-(10) gives:

where H = /a is the Hubble parameter. The effective gravitational constant is = κ2 (1-κ2 ξΦ2 )-1 and the critical values of the scalar field are Φc = ±(κ2 ξ)-1/2 for ξ > 0. Hence, the effective Einstein equations are:

One interesting class of solution is obtained by setting Φ = Φc = ±(κ2 ξ)-1/2 and absence of the gauge coupling and the scalar field χ. This gives R = 6( +2H2 ) = κ2 (ρeff -3peff ) and therefore +2H2 = C where C is a constant. The solution is given by22:


for H < . Equation (15) describes an asymptotic contracting de Sitter spacetime as t, reaching a maximum at t = t0 and then superaccelerating with < 0 tending when t → ∞ to a de-Sitter expanding spacetime.




In fact, by setting Φ = ±Φc, one finds the following effective time dependent density and pressure:

and accordingly the corresponding state equation is

One amazing interesting class of solution is obtained if we conjecture that χ = exp(-Φ). As we expect that |Φc| = (κ2ξ)-1/2 >> 1, then we find:

Consequently for << H2 , w -1 which corresponds for quintessence whereas for h << 1, we obtain w < -1 which corresponds to phantom energy. The equation of state which may be written like:

suggests that the equation of state parameter may varies in time and this is quite appealing. Amazingly, for H = α/t, α > 0 (t is here the cosmological time), we find w -1+2/3α and accordingly for α = 3/2, we get w 0 . This special case corresponds for an accelerated universe dominated by pressureless matter and not dark energy. The case χ = exp(-Φ) results into the special exponential potential:

We plot in what follows the potential V(Φ) for h = 1/2 ,g = 1 for illustration purpose.



It is note worthy that exponential potential plays a crucial role in inflationary cosmology 31,32. The case χ = exp(-Φ) is not the only possible choice of solution. We may discuss the following independent interesting cases:

II-1: We discuss first the case of the scalar field χ = Φm, m is a real parameter. Accordingly, the state equation takes the special form:

Notice that for m < 0 and for |Φc| = (κ2ξ)-1/2 >> 1, equation (23) is approximated again by:

and the same previous arguments may hold. However, for m >> 1, we may approximate equation (24) by:

This case corresponds also for a quintessence dominated universe. The case where χ = Φm results on the following potential

We plot in what follows V(Φ) for h = 4,g = 1 and four different values for m:

II-2: Another interesting cases corresponds for H = β/χ,β is a real parameter augmented by χ = Φm, that is H = β /χ = βmΦ/Φ. Accordingly, equation (19) is reduced to:

As χ = Φm, then at |Φc| = (κ2ξ)-1/2 >> 1, equation (26) is reduced straightforwardly to w = -1.

II-3: At the end, we discuss the case H = eΦ and χ = Φm. Accordingly, equation (19) is reduced again to:

The previous case corresponds for certain modified gravity types having the following corresponding Lagrangians:

More generally, we will discuss the particular case ξ << 1. For this, we assume the following ansatzs Φ = Φ0Hq,χ = χ0Hp where (p,q) R and Φ0 ,χ0 are constants parameters chosen such that at t = 1,Φ = Φ0 = 1,χ = χ0 = 1 for mathematical simplicity. These ansatzs could be important building blocks in constructing the models of dark energy. For these particular solutions, as the universe expands in time, one naturally expects that the cosmological constant and accordingly the scalar curvature should decay, thereby leading to the tiny value observed presently. In fact, cold dark matter particles with time decreasing masses were proved to have an important measurable effect in the dynamical motion of the halo of spiral galaxies. At clusters scale this could have important consequences on dark matter halos (axions and mass varying neutrinos are good candidate for the cold dark matter of our universe).41.42,43 Accordingly, can easily check that for p = q:







The author may check that the quintessence solution is obtained one more if, for instance, h = -g2 , << H and << H. However, for H = α/t, α> 0, equation (28) is reduced to:

For h = -g2 and ξ = 1/6 which corresponds for the conformal coupling case, we obtain:

Notice that w > -1(dark energy dominance) for q2 > 1/6 whereas w < -1 (phantom energy dominance) for q2 < 1/6.

In order to analyze the dynamical equations in some details, we choose the scaling solutions which were proved to be useful to describe many interesting features of early and late-time dynamics of the cosmological evolutions. For this, we choose: Φ = tx, χ = ty and a = tz, (x,y,z) are real parameters. Equations (7)-(10) are now reduced for the case of conformal coupling ξ = 1/6 to:











It is an easy task to check that consistency is obtained if, for instance,R = 0, y = x = -1 and κ2 = t2 (increasing gravitational coupling constant). The case y = x = -1 indicates that χ Φ. Then:

After adding equation (38) to (39) we find:

which gives

which holds for h < 3.57. For h = -1, z 2.18 and for h << -1,z >> 1 which corresponds for an accelerated expansion of the universe. Besides, equation (38) or (39) gives g2 + h = 3z - 2. This fact corresponds for a chaotic quartic inflationary potential V(Φ) = (h+g24 /2. The effective pressure and density vary like:

and accordingly the equation of state parameter for z 2.18 and g2 + h = 3z - 2is:

whereas w -1 for z >> 1 as it is expected.

Finally, it is interesting to see that

and that the present day variation of it = 2H/z is too small for h << -1,z >> 1 in agreement with observational limit.44 In fact, time variation of the effective gravitational constant is tightly constrained by astronomical stellar observations with the Solar system and terrestrial laboratories. The time-variation of the effective gravitational coupling constant is not new as a variation of the gravitational constant is predicted in numerous scalar field tensor theories. Further, that there exit a lot of experimental evidence on the time variation of the 4D gravitational constant: radar ranging data to the Viking landers on Mars, lunar laser ranging experiments, measurements of the masses of young and old neutron stars in binary pulsars. It was recently observed that for late times, a modified cosmology with varying κeff is in accordance with the observed values of the cosmological parameters.44,45



In this section, we give a brief thought on the role of the inflaton dynamical massive field Φ interacting with the massless static scalar field χ in the early evolution epoch of the universe where the scalar field potential is:

The case where χ = Φm results on the following potential

We choose m = 1 for which

The partition function is written as:26

where ∂µ = (∂τ,∇), DΦ = ΠidΦi and:

whereβ is a constant having the dimension of time, is the wave number, ωn = (2π/β)n,n = 0,±1,±2,... is the Matsubara frequency, V is the volume in which we are considering the Fourier modes of the fields, ν here is assumed here to represent the = n = 0-mode of the scalar field assumed to be periodic to the imaginary Euclidean time τ and where the mass of the free part is M2 = 3(h+g2 ) < Φ2 > . All the = 0 modes are in the "ν-part" with d3 xψ = 0. The form of the effective Mexican potential will be responsible of the kind of symmetries we have in the theory. We can now expand the potential like:

The first term is the classical term, the second one is a 1-loop correction, the third one is the 2-loop-correction and so on.



In this work, we have discussed numerous new features of a cosmology a la Brans-Dicke with inflaton scalar field Φ with Higgs-like potential non-minimally coupled to the spacetime curvature R interacting with a minimally scalar field χ with a scalar fields potential of the form V(Φ,χ) = hΦ4/8 +hχ4/8 +g2 Φ2 χ2/2. The cosmological model described in this work corresponds for a homogeneous, isotropic and flat Friedmann-Robertson-Walker spacetime. It has interesting and promising features concerning the description of quintessence/ vacuum energy as non-minimally coupled scalar field as well as the important role of the interacting massive dynamics and massless static fields in the theory described. We have discussed many interesting ansatzs relating the inflaton field to a new (real) source fieldc created through some coupling with the inflaton as g2Φ2χ2, where g is a coupling constant acting as a free-parameter. The first case corresponds to χ = exp(-Φ) and the second case is χ = Φm and both cases where discussed at the critical value of the scalar field |Φc| = (κ2ξ)-1/2 >> 1. For the first case, it was found that for << H2, the equation of state parameter w -1 which corresponds for quintessence whereas for h << 1, we obtain w < -1 which corresponds to phantom energy. Furthermore, it was observed that the equation of state parameter may varies in time and this is quite appealing. One additional result obtained for this special case is obtained if the Hubble parameter decays like H =α/t, α > 0. We found that the equation of state parameter varies like w -1+2/3α and accordingly for α = 3/2, we get w 0 . This particular case is amazing as it corresponds for an accelerated universe dominated by pressureless matter and not dark energy. The case χ = exp(-Φ) results into the special exponential potential described by equation (22). It will be of interest to explore in the future a scalar tensor theory dominated by this type of potential. Whereas for the second case χ = Φm, it was as well observed that similar results occur just like the previous exponential behavior with the major difference of the form of the potential obtained. As it was observed in equation (26), the potential depends on m. For m = 1, the potential is of the form V(Φ) Φ4 which is chaotic-like.

Furthermore, we have made a detailed analysis of the dynamical equations obtained in some details where we have choose the scaling solutions. Other types of solutions may exist but we left this part for a future work, nevertheless, scaling laws were proved to give a good picture of the evolution of the universe. It was found that χ = Φ = t-1 and that the universe is in the stage of accelerated universe dominated by dark energy for z 2.18 and by quintessence for z >> 1. In addition, it was observed that the gravitational coupling constant increases in time whereas its present time-variation is too small for z >> 1, that is, in agreements with terrestrial and astrophysical observations.

We have also analyzed briefly the role of the inflaton dynamical massive field Φ interacting with the massless static scalar field χ in the early evolution epoch of the universe. One may extends all the previous arguments to include a dynamical gauge coupling parameter, e.g. a decrease of the gauge coupling scale due to a dilution of the scalar field. Work in this direction is under progress. It is interesting to have a decreasing gauge coupling constant during the matter evolution epoch of the universe. In fact, the physics responsible for a cosmic time variation of the coupling constants takes place at energies above the unification scale.46-51 In all cases, if g varies with time, this means that there exists some new physics beyond out ordinary physical theory and consequently the description of gravity from Einstein's General Relativity is incomplete. If it is the case, we need to introduce something to complete the theory and makes the scalar field vary.52,53,54,55 Further details and analysis are in progress.



The author would like to thanks the anonymous referees for their useful comments and suggestions.


[1] R R Caldwell, R Dave and P J Steinhardt (1998) Phys. Rev. Lett. 80, 1582.         [ Links ]

[2] L-M Yang et al (2000) Astrophys. J., 530, 17.         [ Links ]

[3] P J Steinhardt, L-M Yang and Zlatev I (1999) Phys. Rev. D59, 123509         [ Links ]

[4] A G Riess et al (1998) Astron. J. 116, 1009.         [ Links ]

[5] S Perlmutter et al (1999) Astrophys. J. 517, 565.         [ Links ]

[6] B R Schmidt et al (1998) Astrophys. J. 507, 46.         [ Links ]

[7] J C Fabris, S V B Goncalves and R de Sa Ribeiro (2006) Gen. Rel. Grav. 38, 495.         [ Links ]

[8] Ph Brax and J Martin (1999) Phys. Lett. B468, 45.         [ Links ]

[9] C Armendariz-Picon , V Mukhanov and P J Steinhardt (2001) Phys. Rev. D63, 103510.         [ Links ]

[10] A Yu Kamenshchik, V Moschella and V Pasquier (2001) Phys. Lett. B511, 265.         [ Links ]

[11] R R Caldwell , R Dave and P J Steinhardt (1998) Phys. Rev. Lett. 80, 1582.         [ Links ]

[12] A Starobinsky and V Sahni (2000) Int. J. Mod. Phys. D9, 373.         [ Links ]

[13] N A Bahcall, J P Ostriker, S Perlmutter and Steinhardt P J (1988) Science 284, 1481.         [ Links ]

[14] P J E Peebles and B Ratra (1988) Astrophys. J. Lett. 325, L17.         [ Links ]

[15] P J Steinhardt , L Wang and I Zlatev (1999) Phys. Rev. Lett. 82 , 896.         [ Links ]

[16] Zlatev I, Wang L and Steinhardt P J (1999) Phys. Rev. Lett. 59, 123504.         [ Links ]

[17] P G Ferreira and M Joyce (1998) Phys. Rev. D58, 023503.         [ Links ]

[18] V Sahni and L Wang (2000) Phys. Rev. D62, 103517.         [ Links ]

[19] J P Uzan (1999) Phys. Rev. D59, 123510.         [ Links ]

[20] F C Carvalho and A Saa , (2004) Phys. Rev. D70, 087302.         [ Links ]

[21] J S Alcaniz (2004) Phys. Rev. D69, 083521.         [ Links ]

[22] V Faraoni (2002) Int. J. Mod. Phys. D11, 471.         [ Links ]

[23] V Faraoni (2003) Phys. Rev. D68, 063508;         [ Links ] see also El-Nabulsi A R (2007) Rom. J. Phys. 52, 1-2, 155.         [ Links ]

[24] Z-K Guo, Y-S Piao and Y-S Zhang (2004) Phys. Lett. B594, 247.         [ Links ]

[25] El-Nabulsi A R (2005) Phys. Letts. B619, 26.         [ Links ]

[26] El-Nabulsi A R (2006) Chin. Phys. Lett. 23, 5, 1124.         [ Links ]

[27] P Brax and J Martin (1999) Phys. Letts. B468, 40.         [ Links ]

[28] G Dvali, G Gabadadze and M Shifman (2001) Mod. Phys. Lett. A16, 513.         [ Links ]

[29] S Tsujikawa, K Maeda and T Torii (1999), Phys. Rev. D60, 063515.         [ Links ]

[30] L Kofman, A Lind and A A Starobinsky (1994) Phys. Rev. Lett. 73, 3195.         [ Links ]

[31] E W Kolb and M S Turner (1990) The Early Universe, Addison-Wesley, Redwood, California.         [ Links ]

[32] A D Line (1990), Particle Physics and Inflationary Cosmology, Harwood, Chur, Switzerland.         [ Links ]

[33] D I Kaiser (1997) Phys. Rev. D56, 6175.         [ Links ]

[34] B A Bassett and S Liberati (1998) Phys. Rev. D58, 021302;         [ Links ] see also W-M Suen and C M Will(1988) Phys. Lett. B205, 447.         [ Links ]

[35] F L Bezrukov hep-ph/arXiv:0810.3165.

[36] M Gaudin et al (2006) Int. J. Mod. Phys. D15, 1387.         [ Links ]

[37] C A R Herdeiro and M Sampaio (2006) Class. Quantum Grav. 23, 473.         [ Links ]

[38] G Calcagni (2006) JHEP05, 012;         [ Links ] G Baym, G Bodeker and L McLerran (1996) Phys. Rev. D46, 662;         [ Links ] W D Garretson, Field G B and Caroll S M (1992) Phys. Rev. D46, 5346;         [ Links ] Barrow J D et al (2002) Class. Quantum Grav. 19, 155.         [ Links ]

[39] A Linde, Particle Physics and Inflationary Cosmology (Harwood, Chur, Switzerland, 1990).         [ Links ] A D Linde, Inflation and Quantum Cosmology (Academic Press, Boston, 1990).         [ Links ]

[40] A Linde (2005) New Astronomy Reviews 49, 35 and references therein.         [ Links ]

[41] R Kallosh , D Linde, S Prokushkin and M Shmakova (2002) Phys. Rev. D66, 123503.         [ Links ]

[42] C Rubano and P Scudellaro (2005) Gen. Rel. Gravit. 37 (3), 52.         [ Links ]

[43] A Bonanno and M Reuter (2002) Phys. Lett. B527, 9.         [ Links ]

[44] R Fardon, A E Nelson and N Weiner (2004) JCAP 0410, 00.         [ Links ]

[45] A Zee (1979) Phys. Rev. Lett. 42, 417;         [ Links ] L Smolin (1979) Nucl. Phys. B160, 253;         [ Links ] S L Adler (1982) Rev. Mod. Phys. 54, 729.         [ Links ]

[46] El-Nabulsi A R (2009) Int. J. Mod. Phys. D18, 5, 691.         [ Links ]

[47] El-Nabulsi A R (2009) Int. J. Mod. Phys. D18, 2, 289.         [ Links ]

[48] El-Nabulsi A R (2010) Comm. Theor. Phys. 53, 869.         [ Links ]

[49] J D Barrow et al (2003) Class. Quant. Grav., 20, 2945.         [ Links ]

[50] Witten E (1981) Nucl. Phys. B188, 513.         [ Links ]

[51] S Deser and B Zumino (1977) Phys. Rev. Lett. 38, 1433.         [ Links ]

[52] D J Gross (1989) Phys. Today 42, no. 12, 9.         [ Links ]

[53] J K Webb et al (1999) Phys. Rev. Lett. 82, 884.         [ Links ]

[54] A S Barabash (2003) Astrophys. Space Sci. 283, 607.         [ Links ]

[55] El-Nabulsi A R (2008) Chin. Phys. Lett. 25, 2785.         [ Links ]



(Received on 11 August, 2009)



* Electronic address: