Massless particle creation in a f(r) accelerating universe S. H. Pereira and J. C. Z. Aguilar Universidade Federal de Itajubá, Campus Itabira Rua São Paulo, 377 35900-373, Itabira, MG, Brazil arxiv:1108.3346v1 [gr-qc] 16 Aug 2011 In this paper we discuss the problem of massless particle creation in a flat, homogeneous and isotropic universe in the context of f(r) theories of gravity. We show that for a βr 4 theory the spectrum of created particles can be exactly calculated in a matter dominated universe and we find that particles of small values of the mode k are produced in the past and particles with large values of k are produced only in the future. The decay of these particles in massive particle-antiparticle pairs can be an alternative mechanism to generate particles of dark matter throughout the universe evolution. Moreover, this theory is naturally accelerated, giving an explanation for the current phase of acceleration of the universe. PACS numbers: I. INTRODUCTION The thermal history of the universe can be divided in three dynamically different stages [1 3]. The early universe is dominated by radiation, where the universe is dominated by ultra-relativistic matter. Then it goes through the matter dominated epoch, where the pressureless components determine the expansion rate. The late universe is believed to be dark energy dominated, where the component with negative pressure leads to an accelerated expansion, indicated by several kinds of complementary astronomical observations [4], and it represents about 70% of the total content in the universe. Another great problem of the actual cosmology is the dark matter, a non-luminous kind of matter that has zero pressure and is necessary to explain the gravitational dynamics of galaxy clusters, representing about 25% of the material content of the universe. What kind of particles can be candidates for the dark matter? It must be not interactive with radiation but gravitationally attractive, like ordinary matter. The origin and the nature of both dark energy and dark matter is still a mystery and its existence is beyond the domain of the standard model of particle physics and general relativity [5]. These are the two of the basic problems of the actual cosmology and several alternative theories has been developed in last decades in order to explain these intricate characteristics of the universe. Recently, non-standard gravity theories have been proposed as an alternative to explain the present accelerating stage of the universe only with cold dark matter, that is, with no appealing to the existence of dark energy. This is naturally obtained in f(r) gravity theories (for a review, see [6]). In such an approach, the curvature scalar R in the Einstein-Hilbert action is replaced by a general function f(r), so that the Einstein field equation is recovered as a particular case [7 15]. In a recent work [16] we have studied the quantum process of particle creation in a radiation dominated universe in view of a f(r) = R +βr n theory, which has the correct limit of the standard general relativity for β = 0. In that paper we have showed that both massive and massless scalar particles can be produced by purely expanding effects, contrary to Parker results [17], that forbids the creation of massless particles in radiation dominated universe. This occurs because the positive- and negative-frequency mode functions for massless fields are not mixed in course of the expansion due to the conformal invariance of the metric. In particular, this means that there is no creation of massless particles, either of zero or non-zero spin. Another way to understand this is that the dispersion relation in Fourier space is exactly the same one of the flat Minkowski space, and, as such, a wave propagation of a massless particle is not accompanied by particle Electronic address: shpereira@gmail.com
production. All this searches are in the framework of standard general relativity (see, for example, [18 20] and references therein, and more recently [21]. See also [22] for several interesting discussion on particle creation.) In the present work we investigate the particle production in a f(r) = βr n theory in a flat and matter dominated universe. We show that scalar massless particles can be produced as an possible alternative to create dark matter particles through the decay of these particles in particle-antiparticle pairs. Moreover, such f(r) naturally provide an accelerated stage for the evolution of the universe, without to appeal to a dark energy component. 2 II. GENERAL THEORY OF SCALAR PARTICLE CREATION The canonical quantization of a real minimally coupled scalar field in curved backgrounds follows in straight analogy with the quantization in a flat Minkowski background. The gravitational metric is treated as a classical external field which is generally non-homogeneous and non-stationary. The basic equation for the study of scalar particle creation in a spatially flat Friedmann-Robertson-Walker geometry is [21]: with χ k (η)+ω2 k (η)χ k(η) = 0, (1) ω 2 k(η) k 2 +m 2 eff and m 2 eff m 2 a(η) 2 a(η) a(η), (2) where k is the Fourier mode of the wave associated to the energy of the particle through the frequency ω k and m eff represents an effective mass of the particle. a(η) is the cosmological scale factor in terms of the conformal time η and the prime denotes derivatives with respect to it. The solutions of the above equation give the positive- and negative-frequency modes. Following standard lines, the quantization can be carried out by imposing equal-time commutation relations for the scalar field χ and its canonically conjugate momentum π χ, namely [χ(x,η),π(y,η)] = iδ(x y), and by implementing secondary quantization in the so-called Fock representation. After convenient Bogoliubov transformations, one obtains the transition amplitudes for the vacuum state and the associated spectrum of the produced particles in a non-stationary background [20, 21]. Usually, the calculations of particle production deal with comparing the particle number at asymptotically early and late times, or with respect to the vacuum states defined in two different frames and do not involve any loop calculation. However, the main problem one encounters when treating quantization in expanding backgrounds concerns the interpretation of the field theory in terms of particles. The absence of Poincaré group symmetry in curved space-time leads to the problem of the definition of particles and vacuum states. The problem may be solved by using the method of the diagonalization of instantaneous Hamiltonian by a Bogoliubov transformation, which leads to finite results for the number of created particles. To proceed, note that Eq. (1) is a second order differential equation with two independent solutions. Each solution χ k must be normalized for all times according to W k (η) χ k (η)χ k(η) χ k(η)χ k(η) = 2i. (3) The Bogoliubov coefficients can be calculated and a straightforward calculation leads to the final expression for the number of created particles in the k mode [20, 21]: N k (η) = 1 4 ω k (η) χ k (η) 2 + ω k(η) χ k (η) 2 1 4 2. (4) The density of particles per mode is given by n k (η) = N k(η) a(η) 3, (5)
3 and the total density of created particle is readily obtained by integrating over all the modes n(η) = n k (η)dk. (6) III. PARTICLE CREATION IN f(r) THEORY Now we will study the particle creation for a simple f(r) theory, discussed in detail by Allemandi et al [9], namely f(r) = βr n. (7) In the Palatini approach the generalized Einstein equations follows from the structural equation R df(r) dr 2f(R) = κτ, (8) where τ = 3p ρ for a perfect fluid and κ = 8πG. In the general case of a fluid equation of state in the form p = ωρ, we have ρ = ρ 0 a 3(1+ω) for the energy density, where ρ 0 is a constant that represents the present day value of the energy density of the fluid. Substituting (7) in (8) we are left with R n = κτ β(n 2), (9) which implies that the case n = 2 is singular. Taking into account the dominant energy condition ρ > 0, we should impose τ < 0, so that the condition ω < 1/3 must be satisfied. R n must be positive definite for even integer values of n and the dimensional constant β should be fixed as β > 0 for n < 2 and β < 0 for n > 2. The Hubble parameter in the flat case can be calculated to be [9] (ȧ ) 2 H 2 = = a ( ) 1 2nε κτ n, (10) 3(3ω 1)[3ω(n 1)+(n 3)] β(n 2) where ε = 1 for odd values of n and ε±1 for even values. The dot denotes derivatives with respect to physical time t. We see that this expression is undefined for the radiation case (ω = 1/3), since Eq. (9) implies R n = 0. This entails that our formalism fails to cover this case. A specific model to treat the particle production on the radiation dominated universe is given in [16]. It is also important to note that we do not need to assume that n is an integer. It can be any real or rational number, since that Eq. (9) be definite and Eq. (10) be positive. Now let us consider the case of matter dominated universe (ω = 0). The modified Friedmann equations can be obtained [9] (apart from integration constants) ( ) n/3 [ ] 1/3 3ε κρ 0 a(t) = t 2n/3. (11) 2n(3 n) β(2 n) Note that the case n = 3 is singular. With this expression we can calculate the deceleration parameter, (ä(t)a(t) ) q(t) ȧ(t) 2 = 3 2n 2n. (12) We see that accelerated solutions (q < 0) are obtained in the cases n < 0 or n > 3/2. In order to illustrate qualitatively the particle creation related to this theory we must to set particular values to the parameters n, β and ε. We choose n = 4, which corresponds to an accelerated a(t) according to (12), β < 0 and ε = 1. We see that the positivity of (9) and of the Hubble parameter (10) are satisfied. With such parameters, the scale factor is ( ) 4/3 ( ) 1/3 a(t) = At 8/3 3 κρ 0 ; A =, (13) 8 2β
and the deceleration parameter q = 5/8. The Hubble parameter H = ȧ/a = 8/3t determines the actual age of the universe as t 0 = 8/3H 0 18.3 10 9 years. It is interesting to note that it does not depends on the β parameter. To follow up we need the scale factor a(t) in terms of the conformal time η. This can be obtained by taking and after integration we are left with: a(η) = t dt η a(t ), (14) ( ) 8/5 B ( η) ; B = 3 A 3/5 ; < η < 0. (15) 8/5 5 In terms of the conformal time, the early universe (the past) corresponds to η and the late universe (the future) corresponds to η 0. Substituting this into the Eq. (1), we find that a /a = 104/25η 2 does not depends on the parameter B, and consequently does not depends on β too. The equation (2) for mode function in the case of a massless particle (m = 0) is [ χ k(η)+ k 2 104 ] 25η 2 χ k (η) = 0. (16) According to the normalization (3), the solution of this equation is given by π η [ χ k (η) = J21 2 (k η )], (17) 10 (k η ) iy21 10 where J ν and Y ν are the Bessel functions of first and second kind, respectively, and ν is its order. This function has the correct asymptotics at early times (η ) of the form χ k 1 k exp ( ikη +δ ), (18) with δ a phase. This corresponds to plane waves for the modes k in the past. Before we proceed with the calculation of the spectrum of non-massive particles created during the evolution of the universe, we must choose a value for the parameter β. We have already seen earlier that the differential equation (1) does not depends on β, so the functions χ k does not depends on this parameter too. However, to calculate the density of particles per mode (5), the parameter β is present on the denominator through the scale factor a(η). Let us choose β such that a(t 0 ) a 0 = 1, for t 0 = 8/3H 0, i.e., the scale factor is normalized to unity for the present time. To satisfy this condition we must have ( ) 8 81 8 β = κρ 0. (19) 8192 3H 0 Using (14) we get the value of the present conformal time as η 0 = 8/5H 0. With this the scale factor is ( ) 8/5 8 1 a(η) =. (20) 5H 0 ( η) 8/5 Now it is easy to calculate the evolution of particle density n k for each mode k, as given by (4) and (5). Fig. 1 shows the density of non-massive particles created during the evolution of the universe as a function of η/ η 0 for different values of the mode k. In the past (η ) the density is zero for all modes, however it grows throughout evolution. The value η/η 0 = 1 represents the present time. Analyzing the figure we can see that for every k there is a critical value of the time η for which the density of created particles grows abruptly and diverges. After that the creation process stops. Mathematically this occurs 4
5 FIG. 1: Spectrum of the density of created particles as function of the conformal time η/η 0 for different values of the mode k for the case n = 4. because in the expression (4) the frequency ω k (η) appears in the denominator, and when η 2 = 104/25k 2 we have ω k 0. Physically, the significance of this divergence is that for values of η 2 < 104/25k 2, the frequency ω 2 k becomes negative, consequently the state of minimum energy and the quantum vacuum are not well-defined in these cases, and the creation ceases for these values. Note that, in principle, an infinite number of particles can be created until the process stops. For k = 1 this occurs in the past, when η/ η 0 2.040. For k = 2 we have η/ η 0 1.019, very close to the present time. But for k = 3 this will only occur in the future, η/ η 0 0.680. We see that the greater the value of k, the divergence will occur in the future, meaning that particles of higher energies can only be created in the future. Until the present moment only particles with k < 2.0396 were created. Remember that k is continuous, because it represents the moment of free particles being created. Some values of k < 1 are also shown in the figure. IV. CONCLUDING REMARKS We haveaddressed the problem of real minimally coupled massless particle creation in a f(r) = βr 4 theory for a matter dominated universe. We have verified that this particular type of f(r) correctly reproduces a universe that evolves in accelerated manner (q = 5/8). In addition we have showed that non-massive particles with mode k < 2.0396 can be created in this particular model. Although they are non-massive particles, depending on its energy, on the framework of the quantum mechanics these particles could decay into massive particle-antiparticle pairs, as occur in pair production of electron-positron in the decay of a energetic photon. The vacuum quantum effects in gravitational field necessary to obtain an essential particle creation requires that the curvature of spacetime must be at least of the order of the inverse Compton length, thus one pair can be created on the Compton length. For a massless field such a threshold is absent and the spectrum of particles created is concentrated in the region of energies of the order of the characteristic curvature radius. Thus the effect will be essential only for a sufficiently large space-time curvature. This mechanism could be responsible, in principle, for the creation of massive particles of dark matter when the curvature radius was large, in the past. For example, analyzing the figure we see that for modes with k < 1 the amount of particles created is steadily increasing in the past, although the creation of such modes stop abruptly in the past. This shows that in the past a huge amount of particles with low k were created. In order to be a source of dark matter particles, the mechanism of decay of these particles in massive particle-antiparticle pairs needs to be studied more carefully.
6 Acknowledgments SHP is grateful to CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico, Brazilian research agency, for the financial support, process number 477872/2010-7. [1] S. Weinberg, Cosmology, Oxford University Press, New York, USA (2008). [2] E. W. Kolb and M. S. Turner, The Early Universe, Westview Press, USA, (1990). [3] V. Mukhanov, Physical Foundations of Cosmology, Cambridge University Press, (2005). [4] A. G. Riess et al., Astron. J. 116, 1009 (1998); S. Perlmutter et al., Astrophys. J. 517, 565 (1999); P. Astier et al., Astron. Astrophys. 447, 31 (2006); A. G. Riess et al., Astro. J. 659, 98 (2007). [5] T. Padmanabhan, Phys. Rept. 380, 235 (2003); P. J. E. Peebles and B. Ratra, Rev. Mod. Phys. 75, 559 (2003); J. A. S. Lima, Braz. J. Phys. 34, 194 (2004); V. Sahni and A. Starobinsky, Int. J. Mod. Phys. D 15, 2105 (2006). [6] T. P. Sotiriou and V. Faraoni, Rev. Mod. Phys. 82, 451 (2010). [7] D. N. Vollick, Phys. Rev. D 68, 063510 (2003); M. Amarzguioui, O. Elgaray, D. F. Mota and T. Multamaki, Astron. Astrophys. 454, 707, (2006); S. Nojiri and S. D. Odintsov, Phys. Lett. B 659, 821, (2008). [8] G. Allemandi, A. Borowiec, M. Francaviglia and S. D. Odintsov, Phys. Rev. D 72, 063505 (2005); L. Amendola, D. Polarski and S. Tsujikawa, Phys. Rev. Lett. 98, 131302 (2007); J. Santos, J. S. Alcaniz, F. C. Carvalho and N. Pires, Phys. Lett. B 669, 14 (2008); J. Santos and M. J. Reboucas, Phys. Rev. D 80, 063009 (2009). [9] G. Allemandi, A. Borowiec and M. Francaviglia, Phys. Rev. D 70, 043524 (2004); [10] T. P. Sotiriou, Phys. Rev. D 73 063515 (2006). [11] X. Meng and P. Wang, Class. Quant. Grav. 21, 2029 (2004). [12] A. D. Dolgov and M. Kawasaki, Phys. Lett. B 573, 1 (2003); T. Chiba, Phys. Lett. B 575, 1 (2003). [13] J. A. R. Cembranos, Phys. Rev. D 73, 064029 (2006); T. P. Sotiriou, Class. Quant. Grav. 23, 1253 (2006). [14] M. Gasperine and G. Veneziano, Phys. Lett. B 277, 256 (1992). [15] G. A. Vilkovisky, Class. Quant. Grav. 9, 895 (1992). [16] S. H. Pereira, C. H. G. Bessa and J. A. S. Lima, Phys. Lett. B 690, 103 (2010). [17] L. Parker, Phys. Rev. Lett. 21, 562 (1968); Phys. Rev. 183, 1057 (1969); Phys. Rev. D 3, 346 (1971); Phys. Rev. Lett. 28, 705 (1972); Phys. Rev. D 7, 976 (1973). [18] N. D. Birrell and P. C. W. Davies, Quantum Fields in Curved Space (Cambridge University Press, Cambridge, 1982). [19] S. A. Fulling, Aspects of Quantum Field Theory in Curved Spacetime (Cambridge University Press, Cambridge, 1989). [20] A. A. Grib, S. G. Mamayev and V. M. Mostepanenko, Vaccum Quantum effects in Strong Fields (Friedmann Laboratory Publishing, St. Petesburg, 1994). [21] V. F. Mukhanov and S. Winitzki, Introduction to Quantum Effects in Gravity, (Cambridge University Press, Cambridge, 2007). [22] Ya. B Zel dovich, Pisma Zh. Eksp. Teor. Fiz. 12, 443 (1970), (English transl. JETP 12, 307 (1970)); A. A. Grib, S. G. Mamayev and V. M. Mostepanenko, Gen. Rel. Grav., 7, 535 (1975); R. Brout, F.Englert and E. Gunzig, Ann. Phys. 115, 78 (1978); A. A. Grib and Yu. V. Pavlov, Grav. Cosmol. 11, 119 (2005); Grav. Cosmol. 12, 159 (2006); AIP Conf. Proc. 956, 96 (2007); M. Bordag, J. Lindig and V. M. Mostepanenko, Class. Quantum Grav. 15 581 (1998); M. Mijic, Phys. Rev. D 57, 2138 (1998); S. Biswas and I. Chowdhury, Int. J. Mod. Phys. D 15, 937 (2006); A. A. Grib and S. G. Mamayev, Yad. Fiz. 10, 1276 (1969)[English transl.: Sov. J. Nucl. Phys. 10, 722 (1970)]; Yu. V. Pavlov, Theor. Math. Phys. 126, 92 (2001).