Abstract
We show how the quantum to classical transition of the cosmological fluctuations produced during inflation can be described by means of the influence functional and the master equation. We split the inflaton field into the system-field (long-wavelength modes), and the environment, represented by its own short-wavelength modes. We compute the decoherence times for the system-field modes and compare them with the other time scales of the model.
Influence functional approach to decoherence during inflation
Fernando C. Lombardo
Departamento de Física Juan José Giambiagi, FCEyN UBA, Facultad de Ciencias Exactas y Naturales, Ciudad Universitaria, Pabellón I, 1428 Buenos Aires, Argentina
ABSTRACT
We show how the quantum to classical transition of the cosmological fluctuations produced during inflation can be described by means of the influence functional and the master equation. We split the inflaton field into the system-field (long-wavelength modes), and the environment, represented by its own short-wavelength modes. We compute the decoherence times for the system-field modes and compare them with the other time scales of the model.
1 Introduction
The emergence of classical physics from quantum behaviour is important for several physical phenomena in the early Universe. This is beyond the fundamental requirement that only after the Planck time can the metric of the Universe be assumed to be classical. For example, the inflationary era is assumed to have been induced by scalar inflaton fields, with simple potentials [1]. Such fields are typically assumed to have classical behaviour, although in principle a full quantum description should be used. The origin of large scale structure in the Universe can be traced back to quantum fluctuations that, after crossing the horizon, were frozen and became classical, stochastic, inhomogeneities [2].
It is generally assumed that several phase transitions have occurred during the expansion of the Universe [3]. As in the case for the inflaton fields, the (scalar) order parameter fields that describe these transitions are described classically. However, the description of early universe phase transitions from first principles is intrinsically quantum mechanical [4]. As a specific application [5] of the previous point, the very notion of topological defects (e.g. strings and monopoles) that characterize the domain structure after a finite-time transition, and whose presence has consequences for the early universe, is based on this assumption of classical behaviour for the order parameter [6], as it distributes itself between the several degenerate ground states of the ordered system.
In previous publications, we analysed the emergence of a classical order parameter during a second order phase transition and the role of decoherence in the process of topological defect formation [7,8,9,10].
In the present paper our concern is directly related with the first point above, the quantum to classical transition of the inflaton. Any approach must take into account both the quantum nature of the scalar field and the non-equilibrium aspects of the process [11]. The problem of the quantum to classical transition in the context of inflationary models was first addressed by Guth and Pi [12]. In that work, the authors used an inverted harmonic oscillator as a toy model to describe the early time evolution of the inflaton, starting from a Gaussian quantum state centered on the maximum of the potential. They subsequently showed that, according to Schrödinger's equation, the initial wave packet maintains its Gaussian shape (due to the linearity of the model). Since the wave function is Gaussian, the Wigner function is positive for all times. Moreover, it peaks on the classical trajectories in phase space as the wave function spreads. The Wigner function can then be interpreted as a classical probability distribution for coordinates and momenta, showing sharp classical correlations at long times. In other words, the initial Gaussian state becomes highly squeezed and indistinguishable from a classical stochastic process. In this sense, one recovers a classical evolution of the inflaton rolling down the hill.
A similar approach has been used by many authors to describe the appearance of classical inhomogeneities from quantum fluctuations in the inflationary era [13]. Indeed, the Fourier modes of a massless free field in an expanding universe satisfy the linear equation
For sufficiently long-wavelengths (k2 << a¢¢/a), this equation describes an unstable oscillator. If one considers an initial Gaussian wave function, it will remain Gaussian for all times, and it will spread with time. As with the toy model of Guth and Pi, one can show that classical correlations do appear, and that the Wigner function can again be interpreted as a classical probability distribution in phase space. (It is interesting to note that a similar mechanism can be invoked to explain the origin of a classical, cosmological magnetic field from amplification of quantum fluctuations.)
However, classical correlations are only one aspect of classical behaviour. It was subsequently recognized that, in order to have a complete classical limit, the role of the environment is crucial, since its interaction with the system distinguishes the field basis as the pointer basis [14]. (We are reminded that, even for the fundamental problem of the space-time metric becoming classical, simple arguments based on minisuperspace models suggest that the classical treatment is only correct because of the interaction of the metric with other quantum degrees of freedom [15].)
While these linear instabilities cited above characterise free fields, the approach fails when interactions are taken into account. Indeed, as shown again in simple quantum mechanical models (e.g. the anharmonic inverted oscillator), an initially Gaussian wave function becomes non-Gaussian when evolved numerically with the Schrödinger equation. The Wigner function now develops negative parts, and its interpretation as a classical probability breaks down [9]. One can always force the Gaussianity of the wave function by using a Gaussian variational wave function as an approximate solution of the Schrödinger equation, but this approximation deviates significantly from the exact solution as the wave function probes the non-linearities of the potential [9,16].
When interactions are taken into account, classical behaviour is recovered only for öpen systems", in which the unobservable degrees of freedom interact with their environment. When this interaction produces both a diagonalization of the reduced density matrix and a positive Wigner function, the quantum to classical transition is completed [17].
In Ref.[9] we have considered an anharmonic inverted oscillator coupled to a high temperature environment. We showed that it becomes classical very quickly, even before the wave function probes the non-linearities of the potential. Being an early time event, the quantum to classical transition can now be studied perturbatively. In general, recoherence effects are not expected [18]. Taking these facts into account, we have extended the approach to field theory models [7,8]. In field theory, one is usually interested in the long-wavelengths of the order parameter. Even the early universe is replete with fields of all sorts which comprise a rich environment, in the inflationary example, we considered a model in which the system-field interacts with the environment-field, including only its own short-wavelengths. This is enough during inflation. Assuming weak self-coupling constant (inflaton potential is flat) we have shown that decoherence is a short time event, shorter than the time tend, which is essentially the time by which different modes in the system sector cross the horizon. As a result, perturbative calculations are justified[8]. Subsequent dynamics can be described by a stochastic Langevin equation, the details of which are only known for early times [19].
In our approach, the quantum to classical transition is defined by the diagonalization of the reduced density matrix. In phase transitions the separation between long and short-wavelengths is determined by their stability, which depends on the parameters of the potential. During Inflation, this separation is set by the existence of the Hubble radius. Modes cross the aparent horizon during their evolution, and they are usually treated as classical. The main mativation of this talk is to present a formal way to understand this statement within the open quatum sytem approach. In the last sense, decoherence is the critical ingredient if we are to dynamically demonstrate the quantum-to-classical transition of the open system.
The paper is organized as follows. In Section 2 we introduce our model. This is a theory containing a real system field f, massless and minimally coupled to de Sitter background. We compute the influence functional by integrating out the environmental sector of the field, composed by the short-wavelength modes. Section 3 is dedicated to reviewing the evaluation of the master equation and the diffusion coefficients which are relevant in order to study decoherence. In Section 4 we evaluate upper bounds on the decoherence times. As we will see, decoherence takes place before the end of the inflationary period. Section 5 contains our final remarks.
2 The Influence Functional and the Density Matrix
Let us consider a massless quantum scalar field, minimally coupled to a de Sitter spacetime ds2 = a(h) [dh2 ] (where h is the conformal time dh = dt/a(t)), with a quartic self-interaction. The classical action is given by
where a(h) = 1/(Hh) and f¢ = df/dh (a(hi) = 1 [hi = H1], and H is the Hubble radius). Let us make a system-environment field splitting
where the system field contains the modes with wavelengths longer than the critical value L1 = 2p/lc, while the bath field contains wavelengths shorter than L1. As we set a(hi) = 1, a physical length lphys = a(h)l coincides with the corresponding comoving length at the initial time. Therefore, the splitting between system and environment gives a system sector constituted by all the modes with physical wavelengths shorter than the critical length lc at the initial time hi.
After splitting, the total action (2) can be written as
where S0 denotes the free field action and the interaction term is given by
The total density matrix (for the system and bath fields) is defined by
where and are the eigenstates of the field operators < and > , respectively. For simplicity, we will assume that the interaction is turned on at the initial time hi and that, at this time, the system and the environment are not correlated (we ignore, for the moment, the physical consequences of such a choice, it has been disussed in [8]). Therefore, the total density matrix can be written as the product of the density matrix operator for the system and for the bath
We will further assume that the initial state of the environment is the vacuum.
We are interested in the influence of the environment on the evolution of the system. Therefore the reduced density matrix is the object of relevance. It is defined by
The reduced density matrix evolves in time by means of
where Jr[h,hi] is the reduced evolution operator
The influence functional (or Feynman-Vernon functional) F[,] is defined as
This functional takes into account the effect of the environment on the system. The influence functional describes the averaged effect of the environmental degrees of freedom on the system degrees of freedom to which they are coupled. With this functional, one can indetify a noise and dissipation kernel related by some kind of fluctuation-dissipation relation. This relation is important when one is interested in possible stationary states where a balance in eventually reached. During inflation we have a very flat potential well away from its minimum, and we are, in general, only interested in the dynamics over some relatively small time. For example, we would neglect dissipation during the slow roll period; but it is not correct during the eventual reheating phase.
We define the influence action dA[,] and the coarse grained effective action (CGEA) A[,] as
We will calculate the influence action perturbatively in l and we will consider only terms up to order l2 and one loop in the expansion. The influence action has the following form
where á ñ0 is the quantum average, assuming the environment is initially in its vacuum state.
The influence functional can be computed, and the result is
where (x y) is the Feynmann propagator of the environment field, where the integration over momenta is restricted by the presence of the infrared cutoff L. We have also defined,
3 Master Equation and Diffusion Coefficients
In this Section we will obtain the evolution equation for the reduced density matrix (master equation), paying particular attention to the diffusion terms, which are responsible for decoherence. We will closely follow the quantum Brownian motion (QBM) example [21,22], translated into quantum field theory [7,24].
The first step in the evaluation of the master equation is the calculation of the density matrix propagator Jr from Eq.(9). In order to solve the functional integration which defines the reduced propagator, we perform a saddle point approximation
where is the solution of the equation of motion d ReA/ = 0 with boundary conditions (hi) = and (h) =. Since we are working up to l2 order, we can evaluate the influence functional using the solutions of the free field equations. This classical equation is + 2 Ñ2f < = 0, ( = a¢(h)/a(h)). A Fourier mode (x) = , satisfies
where we have used = a(h) and the fact that = . It is important to note that for longwavelength modes, k << 2/h2, Eq. (16) describes an unstable (upside-down) harmonic oscillator [12].
As a classical solution can be written as
where
We will assume that the system-field contains only one Fourier mode with = 0. This is a sort of "minisuperspace" approximation for the system-field that will greatly simplify the calculations, therefore we assume
where is given by (17).
In order to obtain the master equation we must compute the final time derivative of the propagator Jr. After that, all the dependence on the initial field configurations (coming from the classical solutions ) must be eliminated. In previous publicatiuons, we have shown that the free propagator satisfies [8]
These identities allow us to remove the initial field configurations , by expressing them in terms of the final fields and the derivatives , and obtain the master equation.
The full equation is very complicated and, as for quantum Brownian motion, it depends on the system-environment coupling. In what follows we will compute the diffusion coefficients for the different couplings described in the previous section. As we are solely interested in decoherence, it is sufficient to calculate the correction to the usual unitary evolution coming from the noise kernels (imaginary part of the influence action). The result reads
where we have defined G1 = and G2 = . V is the spatial volume inside which there are no coherent superpositions of macroscopically distinguishable states for the system field. The ellipsis denotes other terms coming from the time derivative that not contribute to the diffusive effects. This equation contains time-dependent diffusion coefficients Di(t). Up to one loop, only D1 and D2 survive. Coefficient D1 is related to the interaction term , while D2 to . These coefficients can be (formally) written as
and
where the function Fcl is
The explicit expression of these coefficient are complicated functions of conformal time, the particular mode k0, and the cutoff L, and we will show them in a separate publication [23]. For the porpuse of this talk, we will use some analytical approximations, which will allow to obtain an estimation of the scale of decoherence in a particular case.
4 Decoherence
The effect of the diffusion coefficient on the decoherence process can be seen considering the following approximate solution to the master equation
where is the solution of the unitary part of the master equation (i.e. without environment), and Gj includes the coefficients in front each diffusion term in Eq.(20). The system will decohere when the non-diagonal elements of the reduced density matrix are much smaller than the diagonal ones.
The decoherence time-scale sets the time after which we have a classical field configuration, and it can be defined as the solution to
We will solve this equation only in a very approximated way in order to find upper bounds for the decoherence times coming from each diffusion coefficients. A more refined evaluation of decoherence times will be shown in Ref. [23]. For example, when we consider a mode k0< H, we can probe that a good approximation to D1 is given by
This is valid for modes shorter or the same order than H; but it is an overestimation for k0 << H. As D1 is defined under the constraint L/3 < k0 < L, and considering we will set the critical length L < H, ours is a good estimation. Analysing the coefficient D2, it is possible to find an approximated expression for small values of k0 (respect to H). For k0/H > 1, that corresponds to those modes out of the aparent horizon at hi, the diffusion coefficient D2 is an oscilatory function and it has a maximun when k0 ~ L, this is the behavior noted for conformally coupled fields [24]. Finally, for long-wavelength modes, we can write
Both approximations above, are close to the exact coefficients when L < H.
In order to quantify decoherence time we have to fix the values of G1,2. For this, we have to assume values to l, V, D, and S. We will use the more conservative choice in order to have a lower bound to the decoherence time.
Assuming slow roll condition 1/2 (df/dt) << U = lf4, the classical equations (using << U¢) are
these equations are obtained under the following conditions
where = 1/G is the Plank mass.
Definig the end of the inflationary period setting ÎU ~ 1, one can set f(N) »mpl, where N = ln a(hf)/a(h) is the e-fold number. Thus, we assume the mean value of the system field at time of decoherence is f, and we set ( ) ~ 105f, and V ~ H3.
From previous considerations, we can show
and
Using N = H tend> 60 as an estimative scale to the end of inflationary period,
corresponding to the diffusion term D1 and values of l ~ 108. The scale coming from D2 is
From scales tD1, and tD2, we can see that for a given mode k0 < L < H, decoherence is effective by the time in which inflation is ending.
5 Final Remarks
Let us summarize the results contained in this paper. After the integration of the high frequency modes in Section 2, we obtained the CGEA for the low energy modes. From the imaginary part of the CGEA we obtained, in Section 3, the diffusion coefficients of the master equation. System and environment are two sectors of a single scalar field, and the results depend on the "size" of these sectors, which is fixed by the critical wavelength L1.
In Section 4 we analysed the decoherence times for those modes in the system whose wavelength is shorter than the critical value.
We have shown some analitycal approximations that allow to conclude that if we consider a critical length L ~ H, those modes with wavelength k0 << L are the more affected by diffusion throught coefficient D2. For these modes, we can show that the effect is not dependent of the critical L [23].
If one consider a cutoff L> H, and modes H < k0 < L, diffusive effects are larger for those modes in the system whose wavelength is close to the critical L1 [24].
In Ref. [23] we will present the complete expression of the diffusive terms, an also an extensive analysis of the evaluation of the timescale for decoherence in several cases of interest.
Acknowledgments
This work is supported by UBA, CONICET and Fundación Antorchas; Argentina. It is a pleasure to thank the organizers, specially to H.-Thomas Elze, for the invitation and hospitality during the meeting DICE '04 at Piombino. I also appreciate very interesting conversations with B.L. Hu, C. Kiefer, N. Mavromatos, and R. Rivers.
References
[1] A.D. Linde, Particle Physics and Inflationary Cosmology (Harwood, Chur, Switzerland, 1990).
[2] S.W. Hawking, Phys. Lett. B115, 295 (1982); A.A. Starobinsky, Phys. Lett. B117, 175 (1982); A.H. Guth and S.Y. Pi, Phys. Rev. Lett. 49, 1110 (1982).
[3] D.A. Kirzhnits and A.D. Linde, Phys. Lett. B42, 471 (1972).
[4] S.A. Ramsey and B.L. Hu, Phys. Rev. D56, 678 (1997); S.A. Ramsey, B.L. Hu, and A.M. Stylianopoulos, Phys. Rev. D57, 6003 (1998); D. Cormier and R. Holman, Phys. Rev. D62, 023520 (2000), and references therein.
[5] T.W.B. Kibble, Phys. Rep. 67, 183 (1980).
[6] A. Vilenkin, Phys. Rep. 121, 263 (1985); A. Rajantie, Int. J. Mod. Phys. A17, 1 (2002).
[7] F.C. Lombardo, F.D. Mazzitelli, and R.J. Rivers, Phys. Lett. B523, 317 (2001).,
[8] F.C. Lombardo, F.D. Mazzitelli, and R.J. Rivers, Nucl. Phys. B672, 462 (2003).
[9] F.C. Lombardo, F.D. Mazzitelli, and D. Monteoliva, Phys. Rev. D62, 045016 (2000).
[10] R.J. Rivers, F.C. Lombardo, and F.D. Mazzitelli, Phys. Lett. B539, 1 (2002).
[11] E. Calzetta and B.L. Hu, Phys. Rev. D52, 6770 (1995).
[12] A. Guth and S.Y. Pi, Phys. Rev. D32, 1899 (1991).
[13] D. Polarski and A.A. Starobinsky, Class. Quantum Grav. 13, 377 (1996); J. Lesgourgues, D. Polarski, and A.A. Starobinsky, Nucl. Phys. B497, 479 (1997).
[14] C. Kiefer, D. Polarski, and A.A. Starobinsky, Int. J. Mod. Phys. D7, 455 (1998).
[15] J.J. Halliwell, Phys. Rev. D 36, 3626 (1987).
[16] F. Cooper, S.Y. Pi and P. Stancioff, Phys. Rev. D34, 383 (1986).
[17] D. Giulini, C. Kiefer, E. Joos, J. Kupsch, I.O. Stamatescu, and H.D. Zeh, Decoherence and the appearance of a classical world in quantum theory (Springer, Berlin, Germany, 1996).
[18] N.D. Antunes, F.C. Lombardo, and D. Monteoliva, Phys. Rev. E64, 066118 (2001).
[19] A. Matacz, Phys. Rev. D55, 1860 (1997).
[20] G. Aarts and J. Berges, Phys. Rev. Lett. 88, 041603 (2002); W. Buchmüller and A. Jakovac, Nucl. Phys. B521, 219 (1998); G. Aarts and J. Smit, Nucl. Phys. B511, 451 (1998).
[21] W. G. Unruh and W. H. Zurek Phys. Rev. D40, 1071 (1989); A. Caldeira and A. Leggett, Phys. Rev. A31, 1059 (1985).
[22] B.L. Hu, J.P. Paz, and Y. Zhang, Phys. Rev. D45, 2843 (1993); D47, 1576 (1993); J.P. Paz, S. Habib, and W.H. Zurek, Phys. Rev. D47, 488 (1993).
[23] F.C. Lombardo and D. López Nacir, in preparation.
[24] F.C. Lombardo and F.D. Mazzitelli, Phys. Rev. D53, 2001 (1996).
Received on 19 December, 2004
- [1] A.D. Linde, Particle Physics and Inflationary Cosmology (Harwood, Chur, Switzerland, 1990).
- [2] S.W. Hawking, Phys. Lett. B115, 295 (1982);
- A.A. Starobinsky, Phys. Lett. B117, 175 (1982);
- A.H. Guth and S.Y. Pi, Phys. Rev. Lett. 49, 1110 (1982).
- [3] D.A. Kirzhnits and A.D. Linde, Phys. Lett. B42, 471 (1972).
- [4] S.A. Ramsey and B.L. Hu, Phys. Rev. D56, 678 (1997);
- S.A. Ramsey, B.L. Hu, and A.M. Stylianopoulos, Phys. Rev. D57, 6003 (1998);
- D. Cormier and R. Holman, Phys. Rev. D62, 023520 (2000),
- [5] T.W.B. Kibble, Phys. Rep. 67, 183 (1980).
- [6] A. Vilenkin, Phys. Rep. 121, 263 (1985);
- A. Rajantie, Int. J. Mod. Phys. A17, 1 (2002).
- [7] F.C. Lombardo, F.D. Mazzitelli, and R.J. Rivers, Phys. Lett. B523, 317 (2001).,
- [8] F.C. Lombardo, F.D. Mazzitelli, and R.J. Rivers, Nucl. Phys. B672, 462 (2003).
- [9] F.C. Lombardo, F.D. Mazzitelli, and D. Monteoliva, Phys. Rev. D62, 045016 (2000).
- [10] R.J. Rivers, F.C. Lombardo, and F.D. Mazzitelli, Phys. Lett. B539, 1 (2002).
- [11] E. Calzetta and B.L. Hu, Phys. Rev. D52, 6770 (1995).
- [12] A. Guth and S.Y. Pi, Phys. Rev. D32, 1899 (1991).
- [13] D. Polarski and A.A. Starobinsky, Class. Quantum Grav. 13, 377 (1996);
- J. Lesgourgues, D. Polarski, and A.A. Starobinsky, Nucl. Phys. B497, 479 (1997).
- [14] C. Kiefer, D. Polarski, and A.A. Starobinsky, Int. J. Mod. Phys. D7, 455 (1998).
- [15] J.J. Halliwell, Phys. Rev. D 36, 3626 (1987).
- [16] F. Cooper, S.Y. Pi and P. Stancioff, Phys. Rev. D34, 383 (1986).
- [17] D. Giulini, C. Kiefer, E. Joos, J. Kupsch, I.O. Stamatescu, and H.D. Zeh, Decoherence and the appearance of a classical world in quantum theory (Springer, Berlin, Germany, 1996).
- [18] N.D. Antunes, F.C. Lombardo, and D. Monteoliva, Phys. Rev. E64, 066118 (2001).
- [19] A. Matacz, Phys. Rev. D55, 1860 (1997).
- [20] G. Aarts and J. Berges, Phys. Rev. Lett. 88, 041603 (2002);
- W. Buchmüller and A. Jakovac, Nucl. Phys. B521, 219 (1998);
- G. Aarts and J. Smit, Nucl. Phys. B511, 451 (1998).
- [21] W. G. Unruh and W. H. Zurek Phys. Rev. D40, 1071 (1989);
- A. Caldeira and A. Leggett, Phys. Rev. A31, 1059 (1985).
- [22] B.L. Hu, J.P. Paz, and Y. Zhang, Phys. Rev. D45, 2843 (1993);
- D47, 1576 (1993); J.P. Paz, S. Habib, and W.H. Zurek, Phys. Rev. D47, 488 (1993).
- [24] F.C. Lombardo and F.D. Mazzitelli, Phys. Rev. D53, 2001 (1996).
Publication Dates
-
Publication in this collection
06 Sept 2005 -
Date of issue
June 2005
History
-
Received
19 Dec 2004