arxiv:gr-qc/ v1 4 Aug 1995

Transcription

1 QUANTUM EFFECTS IN DSF preprint 95/6, revised version FRIEDMANN-ROBERTSON-WALKER COSMOLOGIES Giampiero Esposito, Gennaro Miele, Luigi Rosa, Pietro Santorelli arxiv:gr-qc/ v1 4 Aug 1995 Istituto Nazionale di Fisica Nucleare, Sezione di Napoli Mostra d Oltremare Padiglione 20, Napoli, Italy; Dipartimento di Scienze Fisiche Mostra d Oltremare Padiglione 19, Napoli, Italy. Abstract. Electrodynamics for self-interacting scalar fields in spatially flat Friedmann- Robertson-Walker space-times is studied. The corresponding one-loop field equation for the expectation value of the complex scalar field in the conformal vacuum is derived. For exponentially expanding universes, the equations for the Bogoliubov coefficients describing the coupling of the scalar field to gravity are solved numerically. They yield a non-local correction to the Coleman-Weinberg effective potential which does not modify the pattern of minima found in static de Sitter space. Such a correction contains a dissipative term which, accounting for the decay of the classical configuration in scalar field quanta, may be relevant for the reheating stage. The physical meaning of the non-local term in the semiclassical field equation is investigated by evaluating this contribution for various background field configurations. PACS numbers: 0260, 0370, 0420, 1115,

2 1. Introduction Quantum effects in Friedmann-Robertson-Walker cosmologies In a recent series of papers [1-3], some of the authors have studied the one-loop effective potential for grand unified theories in de Sitter space. Our main results were a better understanding of the symmetry-breaking pattern first found in [4], a numerical approach to small perturbations of de Sitter cosmologies[2], and the analysis of SO(10) GUT theories in de Sitter cosmologies [3]. However, a constant Higgs field, with de Sitter four-space as a background in the corresponding one-loop effective potential, is only a mathematical idealization. A more realistic description of the early universe is instead obtained on considering a dynamical space-time such as the one occurring in Friedmann-Robertson- Walker (hereafter referred to as FRW) models. Indeed, our early work appearing in [2] tried to study the case of varying Higgs field by introducing a function of the Euclidean-time coordinate which reduces to the 4-sphere radius of de Sitter in the limit of constant Higgs field. Although the approximations made in [2] were legitimate for numerical purposes, the gravitational part of the action, and the one-loop effective potential, were not actually appropriate for studying a dynamical cosmological model. Thus, relying on the work in [5-6], this paper studies the first step towards the completion of our programme, i.e. scalar electrodynamics with a self-interaction term for the complex scalar field. The semiclassical field equations in time-variable backgrounds contain non-local terms, which describe the coupling of the scalar field to the gravitational background. The analysis of these equations is relevant for the reheating mechanism in inflationary cosmology and for the dynamics of dissipation via particle production [7]. Our 2

3 analysis deals with Lorentzian space-time manifolds with FRW symmetries, and does not rely on zeta-function regularization. Hence our geometric framework is substantially different from the Riemannian four-manifolds studied in [1-4], where the metric was positivedefinite. Section 2 derives the field equations for a class of cosmological models where scalar electrodynamics is studied in FRW universes. The Coleman-Weinberg potential, and its correction deriving from the non-local term involving the Bogoliubov coefficients for the coupling of the scalar field to gravity, are obtained numerically in section 3. The semiclassical field equation and the physical meaning of such non-local term are studied in detail in section 4. Concluding remarks are presented in section 5, and relevant details are given in the appendix. 2. Model and field equations For the reasons described in the introduction, we consider a complex scalar field φ with a mass and a self-interaction term, coupled to the electromagnetic potential A µ in curved space-time. Hence the action functional is I L det g d 4 x + boundary terms (2.1) where (cf [6]) [ L ( µ +iea µ )φ ][ ] ( µ iea µ )φ m 2 φ φ 1 4! λ(φ φ) 2 ξrφ φ 1 4 F µνf µν. (2.2) 3

4 With a standard notation, is the Levi-Civita connection on the background spacetime, ξ is a dimensionless parameter, R denotes the trace of the Ricci tensor, and F µν ν A µ µ A ν is the electromagnetic-field tensor. Boundary terms are necessary to obtain a well-posed variational problem, and their form is obtained after integration by parts in the volume integral in (2.1). Covariant derivatives of the scalar field are here used to achieve a uniform notation (cf [6] and [8]). Following [5], we now split the complex scalar field φ as the sum of a variable, realvalued background field φ c, and of a complex-valued fluctuation ϕ, i.e. φ = φ c +ϕ. (2.3) The conformal vacuum [9] is here chosen, and the quantum fluctuation ϕ has vanishing expectation value in such a state, < ϕ >= 0, so that < φ >= φ c. The field equations for A µ,φ c and ϕ are obtained by setting to zero the corresponding functional derivatives of the action. Asfar asthegaugepotentiala µ isconcerned, itisconvenient toimposethelorentz gauge µ A µ = 0. At this stage, to quantize the theory, one can follow the Gupta-Bleuler method, or the Faddeev-Popov procedure, or to eliminate the residual gauge freedom by imposing the relativistic gauge condition proposed by Ford in [6], i.e. X µ A µ = 0, where X is a timelike vector field. Such a field admits a natural form in FRW space-times (see below). Thus, on defining the operator g µν µ ν, the resulting form of the field equations is (cf [6]) ) ) (g µν R µν A ν = 2e 2 A µ φ φ ie (φ µ φ φ µ φ (2.4) 4

5 [ Quantum effects in Friedmann-Robertson-Walker cosmologies +m 2 +ξr λφ2 c e2 < A µ A µ > + 1 ] 24 λ < ϕ2 +(ϕ ) 2 +4ϕϕ > φ c = ie < A µ µ (ϕ ϕ ) > (2.5) [ +m 2 +ξr+ 1 6 λφ2 c ] ϕ λφ2 c ϕ = 0. (2.6) Note that equations (2.4)-(2.6) have been obtained by retaining in the action (2.1) only terms quadratic in the fluctuations ϕ and A µ, and setting to zero all terms involving ν A ν and its covariant derivatives in the field equations. The latter condition is sufficient to derive (2.4). Moreover, we require that µ φ c should be proportional to X µ [6]. Note also that the contribution of [ +m 2 +ξr+ 1 ] 12 λφ2 c φ c has been neglected in the course of deriving (2.6), since equation (2.5) implies that such a contribution is of second order in the quantum fluctuations. By taking the complex conjugate of equation (2.6), and defining ϕ (ϕ 1 +iϕ 2 )/ 2, λ 1 λ/2 and λ 2 λ/6, the addition and subtraction of the resulting equations leads to decoupled equations for the real and imaginary parts of ϕ, i.e. [ +m 2 +ξr+ 1 2 λ jφ 2 c ] ϕ j = 0 for all j = 1,2. (2.7) Moreover, by virtue of our particular gauge conditions, equation (2.5) takes the form (see the appendix) [ +m 2 +ξr λφ2 c e2 < A µ A µ > + 1 ] 4 λ 1 < ϕ 2 1 > +1 4 λ 2 < ϕ 2 2 > φ c = 0. (2.8) 5

6 Interestingly, the effects of quantum fluctuations in equation (2.8) reduce to a linear superposition of the self-interaction term studied in [5] and of the electromagnetic term studied in [6], without any coupling. The term < A µ A µ > is evaluated on considering the integral equation equivalent to (2.4), as shown in [6] and in the appendix. Such an analysis proves that, in the case of a spatially flat FRW background, the renormalized expectation value of A µ A µ in the conformal vacuum can be written as [6] < A µ A µ >= B X µ X ν R µν (2.9) where B = In particular, in a de Sitter universe, < A µ A µ >= 12BH 2, where H is the Hubble parameter [6]. Following [5], we can now write the set of equations leading to the solution of (2.7) and (2.8) in spatially flat FRW backgrounds. Without making any approximation, if one defines (a being the cosmic scale factor) τ t t 0 dy a(y) (2.10) and, for all j = 1,2 [ ( Ω 2 j,k k2 +a 2 (τ) m 2 + ξ 1 )R(τ)+ 1 ] 6 2 λ jφ 2 c (τ) (2.11) the one-loop field equations resulting from (2.7) and (2.8) are (cf (3.32) and (3.33) in [5]) 1 d 2 φ c a 2 dτ + 2 da dφ c 2 a 3 dτ dτ e2 < A µ A µ > φ c + V eff φ c + φ c 4π 2 a 2 2 j=1 1 2 λ j 0 ] dk k 2 Ω 1 j,k [s j,k +Re z j,k = 0 (2.12) 6

7 ( ) d d dτ s j,k = dτ log Ω j,k Re z j,k (2.13) ( ) d d ( dτ z j,k = dτ log Ω j,k s j,k + 1 ) 2i Ω j,k z j,k. (2.14) 2 The form of V eff / φ c is given in the appendix. Equations (2.13) and (2.14) are necessary to find solutions of (2.7) by using Fourier-transform techniques and a suitable change of coordinates, asshownin[5]. Theinitialconditionsfor s j,k andz j,k aretheonesappropriate for the choice of conformal vacuum, i.e. [5] s j,k (τ = 0) = 0 z j,k (τ = 0) = 0. (2.15) Of course, the values of m,ξ and λ j should be now regarded as the renormalized values of such parameters [5]. Note that the second line of equation (2.12) is a typical non-local term, resulting from the self-interaction of the scalar field (see (2.7)) and from its coupling to the geometric background. Moreover, a dissipative term exists which is part of the non-local correction to the Coleman-Weinberg potential, and it is due to the decay processes of φ c in scalar field quanta. By means of such decays, energy is transferred from φ c to the relativistic degrees of freedom ϕ (i.e. radiation). As it is well known, if this release of energy is sufficiently strong, radiation becomes dominant and hence the inflationary phase ends. This leads to the reheating stage, which is as important as the exponential expansion for the dynamics of the early universe. 7

8 3. Numerical evaluation of the one-loop effective potential In this section, for the physically relevant case of an exponentially expanding FRW universe, we compute the non-local term on the second line of (2.12). After integration of such aterm withrespect toφ c, theresulting expressioniscomparedwiththecoleman-weinberg potential V eff. The numerical analysis has been performed for a conformally invariant scalar field (henceξ = 1/6andm = 0),bysolvingequations(2.13)and(2.14)withthehelpoftheNAG routine D02BAF, withinitial conditions (2.15). We have fixed λ = 10 2, which ensures the reliability of the perturbative approach, and the Hubble parameter H = 10 1 M PL, which is an intermediate choice between a chaotic model (H M PL ) and a GUT inflationary phase (H 10 4 M PL ). In the case of exponential expansion in the time variable t, which implies (see (2.10)) a(τ) = a(0) (1 a(0)hτ) (3.1) and for a classical field configuration independent of τ, φ c = φ c0, we express the second line of (2.12) as a function of τ and φ c0. In figure 1 we plot the integral with respect to φ c0 of the above quantity when the conformal time τ varies between 0 and 1/H, which corresponds to a large e-fold number, φ c0 [0,10 M PL ], and a(0) has been set to 1. This configuration for a(τ) and fixed φ c0 corresponds to an exact de Sitter space, but unlike [1-4], with a Lorentzian signature. The non-local term describes an energy exchange between the gravitational field and the field φ (see (2.3)). Such an exchange may lead to dissipative 8

9 or non-dissipative processes in the early universe, depending on the sign of the non-local term, as it will be clear from the analysis in section 4. For a constant background field φ c, only a scale factor varying in time yields the above non-local corrections to the Coleman- Weinberg potential V eff. In figure 2, we plot the Coleman-Weinberg potential V eff (φ c ) (resulting from the integration of (A.7) for µ 1 = µ 2 = M PL and V eff (0) = V 0 = 3MPL 2 H2 /(8π)) for the same choice of parameters. As one can see from figures 1 and 2, the contribution of non-local terms to the one-loop effective potential is very small, and hence its effect on the semiclassical equation of motion is negligible in a first approximation (see section 4). 4. Semiclassical field equation and non-local effects The analysis of the previous section has shown that the effect of the Bogoliubov coefficients on the Coleman-Weinberg potential does not modify the pattern of minima found in static de Sitter space [8]. It is therefore legitimate to study the equation (2.12) when, in a first approximation, its second line is neglected, and then to use the resulting solution to evaluate the non-local correction. Here, the value of V eff (0) is set equal to V 0 as above, and it dominates the energy, as it occurs in a de Sitter universe, for φ c < M PL. On studying the limiting form of equation (2.12) when the second line is neglected, it can be easily seen that the self-interaction term (see (A.7)) plays a key role in obtaining a sensible physical model. In other words, for vanishing λ, equation(2.12) admits a run-away 9

10 exact solution for φ c in the form (denoting by D 1 and D 2 two integration constants) φ c = D 1 exp( β 1 Ht)+D 2 exp( β 2 Ht) (4.1) where β 1 (3+ ) 9+48Be 2 /2 and β 2 (3 ) 9+48Be 2 /2, and we have re-expressed τ in terms of t by means of (2.10). This behaviour results from the particular form taken by the potential in (2.12), i.e. 6BH 2 e 2 φ 2 c, which is unbounded from below. We thus study the equation given by the first line of (2.12) with non-vanishing λ, i.e. 1 d 2 φ c a 2 dτ dadφ c a 3 dτ dτ e2 < A µ A µ > φ c + V eff = 0 (4.2) φ c subject to the initial conditions φ c (0) = M PL and dφ c /dτ(τ = 0) = 0. In figures 3 and 4, φ c (τ) and dφ c /dτ are plotted. As one might have expected, φ c (τ) decreases until it reaches the zero value, by virtue of the nature of the potential (see figure 2), whilst the kinetic energy increases. We have then used the solution of equation (4.2) to evaluate the non-local effects in our model. For this purpose, starting from the definitions of energy density ρ φc and pressure p φc of the background scalar field in a de Sitter universe, i.e. ρ φc 1 2 φ 2 c +V eff 6Be 2 H 2 φ 2 c (4.3) p φc 1 2 φ 2 c V eff +6Be 2 H 2 φ 2 c (4.4) one gets from (2.12) the equation dρ φc dτ +3Ha(ρ φc +p φc ) = Θ(φ c (τ),τ) dφ c dτ. (4.5) 10

11 With our notation, Θ corresponds to the whole second line of (2.12), Θ(φ c (τ),τ) φ c 4π 2 a 2 2 j=1 1 2 λ j 0 ] dk k 2 Ω 1 j,k [s j,k +Re z j,k. (4.6) If the right hand side of (4.5) can be viewed as a dissipative term, which implies that it always takes negative values and depends quadratically on dφ c /dτ, the Bianchi identity leads to the following equation for the energy density ρ R of radiation: dρ R dτ +4Ha ρ R = Θ(φ c (τ),τ) dφ c dτ. (4.7) Infigure5weplottherighthandsideof(4.5)corresponding tothesolutionofequation (4.2). Note that its sign turns out to be positive for the largest part of the τ-interval, and hence cannot actually lead to dissipative effects. A naturally occurring question is how to interpret this lack of dissipation. Indeed, in the adiabatic approximation studied in [5], the coupling of the scalar field to a fermionic field by a Yukawa term produces a vacuum energy loss rate proportional to (Haφ c +dφ c /dτ)dφ c /dτ. In this last expression the second term is clearly a dissipative effect, since it does not depend on the sign of the velocity and in the formula reported in [5] it occurs with the correct sign to represent an energy loss. By contrast, the first term does depend on the sign of dφ c /dτ and in de Sitter, where H is constant, it only represents a further quantum correction to the energy and pressure of the φ c fluid. These considerations seem to suggest that also in our case, where one deals with decays of the classical field configuration into its quanta, the non-local term in the semiclassical field equation may be essentially a linear combination of φ c and dφ c /dτ. If this property 11

12 holds, one can expect that the evaluation of (Haφ c +dφ c /dτ)dφ c /dτ may indicate when dissipative effects are likely to occur, depending on whether the first or the second term of this linear combination is dominant. In other words, the sign of the right hand side of (4.5) shown in figure 5 can be understood by pointing out that for the solution of (4.2), shown in figures 3 and 4, the ratio (dφ c /dτ)/haφ c is smaller than 1 whilst dφ c /dτ is negative, when τ [0,9.2 M 1 PL ]. Note that, in the neighbourhood of τ = 9.2 M 1 PL, Θ dφ c/dτ vanishes, and this corresponds to the value of τ such that the linear combination (Haφ c +dφ c /dτ) vanishes. As a further check of the conjecture about the functional dependence of the non-local term of equation (2.12) on φ c and dφ c /dτ, one should also analyze other situations where, unlike before, the ratio (dφ c /dτ)/haφ c is larger than 1. This can be done for example by taking as a trial function for φ c (τ) the linear combination φ c (τ) = φ c0 +ατ, where α and φ c0 are two arbitrary parameters. The resulting analysis, performed by varying α and φ c0 so as to reproduce the conditions (dφ c /dτ)/haφ c > 1 or (dφ c /dτ)/haφ c < 1, seems to confirm that, in de Sitter, Θ(φ c (τ),τ) is indeed dominated by a term proportional to the combination (Haφ c +dφ c /dτ). 5. Concluding remarks This paper has studied scalar electrodynamics in a spatially flat FRW universe, by including a self-interaction term for a conformally invariant scalar field (cf [5,6,10]). On imposing the Lorentz gauge and the supplementary condition (A.4), the one-loop equation 12

13 for the expectation value of the complex scalar field in the conformal vacuum shows a linear superposition of the self-interaction term of [5] and of the electromagnetic term [6]. The results of our investigation are thus as follows. First, the numerical solution for the Bogoliubov coefficients in (2.13) and (2.14), subject to the initial conditions (2.15), has been obtained and used to find the correction to the Coleman-Weinberg effective potential by integrating with respect to φ c the second line of (2.12). The pattern of minima in the effective potential is not modified by the non-local term in (2.12). For the particular values of parameters considered in our investigation, the non-local corrections turn out to be several orders of magnitude smaller. Second, the limiting form of equation (2.12), i.e. equation (4.2), has been studied, and its numerical solution has been used to evaluate non-local effects in our cosmological model. Such a solution corresponds to a slow-roll dynamics. Interestingly, the approximate calculation of the function Θ defined in (4.6) shows that Θ does not lead necessarily to dissipative effects in the early universe. Nevertheless, in a de Sitter model, the right hand side of equation (4.5) is very well approximated by the same combination of φ c and dφ c /dτ which results from the adiabatic case studied in [5], where the coupling of a scalar field to a fermionic field was instead considered. More precisely, we have found that, for various forms of φ c (τ), the right hand side of (4.5) is essentially given by ( A Ha φ c + dφ c dτ ) dφc dτ where A is a positive function of τ. The first term in round brackets leads to a further quantum correction to the energy density of the scalar field, at least when the function A is slowly varying. The second term is purely dissipative (see the end of section 2). 13

14 One of the main motivations for studying non-local corrections to the Coleman- Weinberg potential was their possible relevance for the reheating of the early universe [7]. Our analysis confirms that they cannot be simply re-expressed by a term of the kind Γdφ c /dτ, as first found in [7], where a different approach to the effective action is used with respect to [5]. A non-trivial open problem is the numerical evaluation of the function A in our de Sitter model. Moreover, it appears necessary to understand how the form of A depends on the specific choice of the background field φ c (τ). Last, but not least, the whole analysis should be repeated for FRW models which are not in a de Sitter phase. Acknowledgments We are grateful to Renato Musto for useful discussions on scalar electrodynamics, and to the European Union for partial support under the Human Capital and Mobility Programme. Anonymous referees made comments which led to a substantial improvement of the original manuscript. Appendix Following[6], wedenotebyd µν R (x,x )thephotonretardedgreen s function, whichsatisfies the equation ( δ µ ρ µ ρ R µ ρ ) x D ρν R (x,x ) = g µν δ(x,x )/ det g. (A.1) 14

15 Thus, by writing A µ in (x) for the photon in-field, the solution of (2.4) to order e, in the Lorentz gauge µ A µ = 0 (A.2) is given by [6] A µ = A µ in ie d 4 x ) det g(x ) D µν R (x,x ) (φ ν φ φ ν φ. (A.3) in The residual gauge freedom of the problem is dealt with by imposing the additional condition [6] X µ A µ = 0 (A.4) where X is the same timelike vector field appearing in (2.9). Equations (A.2) (A.4) imply that, on inserting (A.3) into the right-hand side of (2.5), the only non-trivial contribution still vanishes after integration by parts, since D µν R νφ c = 0 µ D µν R = 0. (A.5) In the equation (2.12) for the expectation value φ c of the scalar field φ, the derivative with respect to φ c of the one-loop effective potential is given, in the case of non-vanishing renormalized mass, by (cf (3.15) of [5]) V eff = m 2 φ c +ξrφ c + 1 φ c 12 λφ3 c 1 (ξ 96π 2λ 1 ) Rφ c ( ( 384π 2λ m 2 + ξ 1 ) R+ 1 ) 6 12 λφ2 c φ c log + 1 ( ( 128π 2λ m 2 + ξ 1 ) R+ 1 ) 6 4 λφ2 c φ c log π 2λ2 φ 3 c m 2 +(ξ (1/6))R+(1/12)λφ 2 c m 2 m 2 +(ξ (1/6))R+(1/4)λφ 2 c m 2 (A.6)

16 and for m = 0 by (cf (3.16) of [5]) V eff = ξrφ c + 1 ( φ c 12 λφ3 c + 1 φ π 2λ2 c + 1 ( φ 3 512π 2λ2 c log + 1 (ξ 384π 2λ 1 ( )Rφ c (ξ 128π 2λ 1 ( )Rφ c 6 log (ξ (1/6))R+(1/4)λφ 2 c (λ/4)µ 2 1 log log (ξ (1/6))R+(1/12)λφ 2 c (λ/12)µ ) 3 (ξ (1/6))R+(1/12)λφ 2 c (ξ (1/6))µ 2 2 (ξ (1/6))R+(1/4)λφ 2 c (ξ (1/6))µ 2 2 ) 1 11 ) 3 ). 1 (A.7) In agreement with the notation of section 2, m,ξ and λ are the renormalized values of our parameters, and µ 1 and µ 2 are completely arbitrary renormalization points. Our equations (A.6) and (A.7) are obtained by imposing the renormalization conditions (3.12a) (3.12c) of [5] and bearing in mind that the numerical coefficients in our equation (2.8) differ from the ones occurring in equation (2.19) of [5], since we study a complex scalar field. References [1] Buccella F, Esposito G and Miele G 1992 Class. Quantum Grav [2] Esposito G, Miele G and Rosa G 1993 Class. Quantum Grav [3] Esposito G, Miele G and Rosa G 1994 Class. Quantum Grav [4] Allen B 1985 Ann. Phys., NY [5] Ringwald A 1987 Ann. Phys., NY [6] Ford L H 1985 Phys. Rev. D [7] Boyanovski D, de Vega H J, Holman R, Lee D S and Singh A 1995 Phys. Rev. D

17 [8] Allen B 1983 Nucl. Phys. B [9] Birrell N D and Davies P C W 1982 Quantum Fields in Curved Space (Cambridge: Cambridge University Press) [10] Shore G M 1980 Ann. Phys., NY Figure captions: Figure 1. In Planckunits (used also infigures 2,3,4and 5), thecorrection tothe Coleman- Weinberg potential resulting from the Bogoliubov coefficients in (2.12) (2.14) is plotted versus τ and φ c0. Figure 2. The Coleman-Weinberg one-loop effective potential is plotted versus φ c0, in the case of a de Sitter background. Figure 3. The solution of the semiclassical field equation (4.2) for the background field configuration is plotted versus τ. Figure 4. For the numerical solution plotted in figure 3, its derivative with respect to τ is here shown. Figure 5. The right hand side of equation (4.5) is plotted versus τ. 17

18

19

20

21

22