QFT in curved space-time
1 Introduction
1.1 The Einstein-Hilbert action (31.01.2010) Before we turn to quantum field theory in curved space-times, we shortly review the classical theory of gravity. Einstein’s equations of classical general relativity can be derived from the Einstein-Hilbert action: 1 Z √ S = d4x g(R − 2Λ). 16πG with R the Ricci curvature and the cosmological constant Λ. Under a variation δgab, the variation of the action is given by: 1 Z √ √ √ δS = d4x δ ( g)(R − 2Λ) + gδ (R ) gab + gR δ gab . 16πG ab ab To analyze this expression, use:
h b b i δg = det (gab + δgab) − det gab = det gab det g c + δg c − 1 = b ab ab = det gab · Trδg c = gg δgab = −ggabδg ab ab c c ab c c g δRab = g (∂cδΓab − ∂bδΓac + ΓδΓ) = g (∇cδΓab − ∇bδΓac) ab c ab c = ∇c g δΓab − g δΓbc
For the last derivation, we used normal coordinates around some point p, i.e. Γ|p = 0 to get rid of terms quadratic in the Christoffel symbols. Next, we replaced the partial derivatives by covariant derivatives for the same reasons. Since the equation is now tensorial again, it holds bc everywhere. In the last step, we used the property of the Levi-Cività connection, ∇ag = 0. A slightly different derivation is given in [15], appendix E. The total derivative contributes only a boundary term, which can be eliminated adding an R a b additional piece 2 ∂M K where K is the extrinsic curvature, K = h b∇an , where hab = a gab ± nanb is the induced metric on the boundary, and n is a normal vector on the boundary. Ignoring these subtleties (for example for a manifold without boundary), we get:
1 Z √ 1 δS = d4x g R − Rg + Λg δgab 16πG ab 2 ab ab resulting in Einstein’s equations in vacuum after setting the variation to zero. To add matter simply add another term to the action, p.e. a field strength term of Yang-Mills gauge theory or electromagnetism: Z √ 1 Z 1 Z √ S = d4x gL = Tr F ∧ ?F = d4x g Tr F abF . M m 4g2 4g2 ab √ δ( gLm) √2 Define the energy-momentum tensor by Tab = g δgab and get Einstein’s equations 1 R − Rg + Λg = 8πGT . ab 2 ab ab ab
1 1.2 Attempts to quantize general relativity
General relativity is a classical field theory, with metric field gab, which plays a dual rôle, both being the dynamical variable in the game, and describing the background. Over the last decades there were various attempts to find a quantum theory of gravity. One of the major lines are the following:
(0) 1. Covariant perturbation method: gab = gab +hab, where the first part solves the classical Einstein equations. Write down Feynman rules. The problem of this attempt is that the perturbation theory is non-renormalizable. Let’s try to understand this. In natural units c = ~ = 1, the action is dimensionless. In dimension n, R has dimension −2 n n−2 L , while the measure contributes L , so Newton’s constant GN has dimension L . n−2 In fact, GN = lp in these units, where lp is the n dimensional Planck length. So our coupling parameter is:
n−2 1 1 1 Mp 2 = = n−2 = κ 16πGN 16πlp 16π
with Planck mass Mp. So at a distance X-times smaller (or at momentum X-times n−2 n−2 larger), GN seems X -times larger, so GN ∼ p and the critical dimension is 2. So in 4d, gravity seems non-renormalizable, i.e. gets stronger at shorter distances. In 4 2 dimensions GN = lp. Pure GR is finite to one-loop order (i.e. the divergences cancel due to an analogue of the Ward identity) but not in higher orders or if coupled to matter. Next step: supergravity. Make the theory supersymmetric, add a gravitino (spin 3/2), to cancel some divergences. Highest dimension (with one time dimension) is 11, so this is the most natural setup. Get solutions in 4 dimensions by dimensional reduction (like in Kaluza-Klein theories). Finally: Perturbative String theory.
2. Canonical quantization method: Wheeler-DeWitt equation: Hb |Ψi = 0. Background independance as a principle. This approach leads to Loop Quantum Gravity after introducing Ashtekar’s variables.
3. Path integral method: Try a generating functional: Z −S Z = D [gab] e
with the Einstein-Hilbert action S as above. It is not known how to define a measure on the moduli space of metrics. Another problem is again non-renormalizability. In LQG, the moduli space is reduced, one ends up with the spin foam formulation.
4. Alternate routes: Twistor theory, gravitational OR (Roger Penrose), non-commutative geometry (Alain Connes).
1.3 One step back: QFT in curved space-time Lacking a full theory of gravity, let’s see how far we get without it. We saw above, that ~ (reinserting ) the Einstein Hilbert action in four dimensions is proportional to 2 . Looking ~ lp
2 at the formal path integral
Z Z − ~ R R −S l2 D[gab] e = D[gab] e p
it seems that quantum effects of gravity become important only at the order of the Planck length. A first step towards quantum electrodynamics is to consider matter in a classical background of an electromagnetic field. In analogy, we might hope to get some insight to quantum effects of gravity by considering quantized fields in a classical background, i.e. in a classical curved space. Still there are some issues with this idea. The principle of general covariance states that matter and mass couple equally strongly to gravity. Nevertheless, if we consider for example a free field theory, only one-loop 1PI diagrams are possible (lacking interaction vertices). Up to order ~, we can therefore truncate gravity also to first-loop order. For an interacting theory of matter fields, higher loops of these fields are possible, so we should also better include higher graviton loops. Recall the very different nature of the coupling constant in, say, electromagnetism, and in Einstein gravity. In the first case, it is given by the 1 dimensionless fine-structure constant α ' 137 , while for gravity, the coupling is dimension- full, being proportional to the Planck length. While one single graviton loop is of zeroth order in GN , higher loops come from graviton interactions, contributing powers of GN . Hence for 2 2 GN lp e length scales l satisfying l2 = l2 α = 4π , quantum effects of matter dominate, and we can truncate gravity at one-loop order. It seems that this should give reasonable results, though there is still much uncertainty about the domain of validity. To sum up, there are several steps forward towards a quantum theory of gravity:
1. Quantum field theory in curved space-time: the background space-time is classical, meaning we work in zeroth order in ~. We ignore the back-reaction of the matter on space-time.
2. Semi-classical gravity: still treat the background as classical, but now take the back- reaction into account. One important element at this stage is the expectation value of the energy-momentum tensor (being quantized) in some state ψ. It contributes to Einstein’s equation: 1 D E Rab − Rgab = 8πG Tˆab . 2 ψ 3. Include also higher graviton loops.
We concentrate here on the first two stages, which already lead to interesting results, like particle creation, Unruh effect, and Hawking radiation.
2 Particle creation in a time-varying background (01.02.2010)
This subsection follows [15]. To do perturbative QFT (at least in the classical fashion), one needs a notion of particle states. Given such states, the S-matrix yields probabilities for a transition between a certain in-state as an element of some Fock space into an out-state in another Fock space. In particular, we want to define a vacuum state (a state without particle excitations). The notion of a particle is no longer automatic in a curved background. As we will see later on,
3 different observers might well disagree over the number of excitations. Only under certain assumptions, it is possible to define meaningful particle states. Possible setups are:
• asymptotically flat space-time, or
• non-vanishing curvature only in a compact region, or
• globally hyperbolic (space-time possesses a Cauchy hypersurface), or
• asymptotically stationary space-time (time-like Killing vector field gives rise to a split into positive and negative frequencies, see below), or
• Feynman propagator (i.e. Green’s function) exists (propagating positive frequency in- states into positive frequency out-states).
To define positive and negative frequencies, we want a positive-definite inner product. For concreteness, consider a Klein-Gordon field. There exists a conserved current, which gives rise to an inner product, the so-called Klein-Gordon product. We will come back to details later on. If we restrict to fields of positive frequency (the Fourier transform Φ(e ω, ~x) vanishes for ω < 0), the product is positive-definite. The Hilbert space of possible states is given by the symmetric Fock space FS(H) with vacuum state |0i = (1, 0, 0,...), where the i-th component in the infinite dimensional vector is in Ji−1 H. The field operator can now be given (in a distributional sense) by ˆ X ∗ − ∗ + Φ(x) = vi (x)ˆa + vi (x)ˆa i + − where a , a are creation and annihilation operators, and vi denotes an orthonormal basis of H. Specify Hilbert spaces Hin and Hout. The S-matrix describes, how a state evolves, p.e. for an in-state |0i ∈ FS(Hin), the out-state S |0i ∈ FS(Hout) tells about spontaneous creation of particle due to a time-dependant gravitational field. The latter state can be calculated. It vanishes for an odd tensor product of Hilbert spaces, hence particles can only be created in pairs. Particle creation only occurs, if some negative-frequency part is picked up.
3 Quantum fluctuations (13.02.2010)
We will from now on closely follow Mukhanov’s book [13]. p2 mω2x2 Let’s start with an harmonic oscillator in a flat background with Hamiltonian H = 2m + 2 . The classical equations of motion are x¨ + ω2x = 0. 2 2 ~ mω The Schrödinger equation is: − 2m ∆ + 2 |ψi = E |ψi. Solve for the wavefunction: r 1 mω 1/4 mωx2 rmω ψn(x) = exp − · Hn x , n = 0, 1, 2,.... 2nn! π~ 2~ ~ In particular, the vacuum state is given by1:
mω 1/4 mωx2 ψ0(x) = exp − . (3.1) π~ 2~
1 p mω i dψ0(x) alternative derivation: a |0i = 0, where a = x + p . This leads to: xψ0(x) + } = 0. 2~ mω mω dx
4 The vacuum state has expectation value hxi = 0. Nevertheless, there is a fluctuation δq ∼ q ~ mω . Now consider the field theory case. Here we have "one oscillator at each point". In momentum 2 2 2 space: Φ¨ k + (k + m )Φk = 0, leading to a wave functional : 1 Z Ψ[Φ] ∼ exp − d3k |Φ |2 ω . 2 k k In analogy, there occur vacuum fluctuations p 1 δΦk = h|Φk|i ∼ √ . ωk 1 If we consider fluctuations of a field in a box with length L, the fluctuations occur at k ∼ L and are given by: ( q L−1 for L 1 δΦ ∼ (δΦ )2 k3 ∼ m . L k −3/2 1 L for L m One consequence of these fluctuations is the Casimir effect (see below).
4 Driven harmonic oscillator
In QFTCS, a time-varying background geometry leads to particle excitations. A lower- dimensional analogy is a harmonic oscillator coupled to a time-varying source term. Take the Lagrangian: 1 1 L = q˙2 − ω2q2 + J(t)q. 2 2 Assume further that J(t) is non-vanishing for t ∈ [0,T ] only3. Now a− satisfies: i a˙ − = −iωa− + √ J(t), 2ω which can be integrated. If we define Z T i iωt J0 = √ dt e J(t), 2ω 0 a straight forward calculation shows, that the occupation number operator has expectation values4: 2 h0in|Nˆ(t)|0ini = |J0|
h0out|Nˆ(t)|0ini = 0. Hence the energy expectation value gets shifted due to the interactions: 1 h0 |Hˆ (t)|0 i = + |J |2 ω. in in 2 0
2 compare equation (3.1), with m = 1, ~ = 1, ω → ωk, x → Φk. 3This corresponds later on to a time-dependant gravitational field, which is stationary in the beginning and the end, to get notions of particles. 4 We are in the Heisenberg picture. If we start with the vacuum |0ini, the state stays in |0ini. So the first line indicates an exictation of the vacuum. However, after the interaction, the true vacuum is given by |0outi. This is reflected in the second equation.
5 5 Quantization of fields
Now we turn to quantum field theory. As stated above, there is "one harmonic oscillator at each point": 1 Z Z S[φ] = d4x φ˙2 − d4x d4y φ(x)M(x, y)φ(y). 2 2 M gets diagonalized to M(x, y) = −∆x + m δ(x−y). Due to the Laplacian, the oscillators are coupled. We can decouple them by a Fourier transformation. Next we perform a mode expansion in eigenfunctions in momentum space. The quantized solution for the field operator φˆ is given by: Z 3 ˆ d k 1 ∗ ik·x − −ik·x + φ(x) = √ v (t) e aˆ + vk(t) e aˆ , (5.1) (2π)3/2 2 k k k where the mode functions satisfy:
2 v¨k + ωk(t)vk = 0. In the case of a free field in flat space, one gets: 1 iωkt vk(t) = √ e , (5.2) ωk √ 2 2 with ωk = k + m . In particular, the mode functions are isotropic (independent of the direction of k). The vacuum energy is in general infinite. Each point is an oscillator by itself and therefore 1 contributes 2 ωk. In the free field case in flat space, one can renormalize the Hamiltonian by imposing a normal-ordering condition, thus substracting the vacuum energy. Another way to understand the appearence of mode functions goes as follows (see [6]): 1. Start with a Lagrangian, p.e. 1 L = ∂ φ∂µφ − m2φ2 − ξRφ2 2 µ ξ = 0 gives the standard Klein-Gordon equation and therefore easy equations of motion. λ For ξ = 0 the field is conformally coupled, so the action is invariant under gµν → e gµν. For the special case of conformally flat metrics (like FLRW metrics), the field decouples from gravity, since its action is equivalent to an action in flat space.
2. Derive the field equations5: 2 φ + m φ + ξRφ = 0. Define an inner product on solutions. To this end, pick a Cauchy surface Σ with normal unit vector nµ, and set:
Z ↔ Z ∗ µ ∗ ∗ µ (v1, v2) = i v2 ∂µ v1 dΣ n = i (v2∂µv1 − v1∂µv2) dΣ n ,,
where v1 and v2 are solutions to the equations of motion. By the equations of motion, the definition does not depend on the specific choice of
5 µ µ Here = ∇ ∇µ. For scalar fields: φ = ∇ ∂µφ.
6 Cauchy surface. Therefore, we get a well-defined inner product, the Klein-Gordon prod- uct. Note that the Klein-Gordon Product is a generalization of the Wronksian v0 v∗ − v0 ∗v W [v , v∗] =(v0 v∗) = k k k k = k k . k k 2i 2i 3. Now proceed with canonical quantization:
φ(x, t), π(x0, t) = iδ(x0, x), where δL π = , δφ˙ Z δ(x0, x) dΣ = 1.
If (vi) is a complete set of solutions of positive norm ("positive energy" in the flat case), then as in conventional QFT, φ gets operator valued, and can be expanded as:
ˆ X ∗ − + φ = vi aˆi + viaˆi , i where a+, a− are creation and annihilation operators. In the non-discrete case, we are back in (5.1).
The problem of curved space-times is the ambiguity in a choice of basis vi. While in the flat case, one picks positive frequency solutions (5.2), there is no canonical choice in the curved case.
5.1 Fields in FLRW models Now consider a field φ with Lagrangian (ξ = 0 above) 1 L = ∂ φ∂µφ − m2φ2 2 µ in a FLRW universe with metric
ds2 = dt2 − a2(t)dx2 = a2(η) dt2 − dx2 .
It is equivalent to a field χ = aφ in flat spacetime, where χ obeys the equation
00 2 χ − ∆χ + meff χ = 0, 2 2 2 a00 where meff = a m − a and the prime denotes taking the derivative with respect to the 2 2 2 1 conformal time η. For ξ non-zero, meff = a m + ξ − 6 R . Again, the oscillator decouple after a Fourier transformation. Now the mode functions satisfy
00 2 vk + ωk(η)vk, q 2 2 with ωk(η) = k + meff(η). This equation has two linear independent solutions, which we interpret as a complex solution vk. Making use of the Klein-Gordon product, we normalize:
(vk, vk0 ) = δkk0 .
7 5.2 Bogolyubov transformations For general manifolds, two isotropic mode functions are related via
∗ ∗ vk = αkuk + βkuk,
2 2 with Bogolyubov coefficients αk, βk. Normalization is kept, if |αk| − |βk| = 1. The operators are related via the Bogolyubov transformation:
ˆ− − ∗ + bk = αkaˆk + βkaˆ−k, ˆ+ ∗ + − bk = αkaˆk + βkaˆ−k. The b vacuum can be expressed in terms of the a vacuum:
" ∗ # Y 1 βk + + |0bi = p exp − aˆk aˆ−k |0ai . 2αk k |αk|
2 1 The product converges, if |βk| → 0 faster than k3 . This is the mean density of b-particles in 2 (3) the mode χk in the a-vacuum, since h0a|Nˆb|0ai = |βk| δ (0). The complete mean density of R 3 2 b-particles is d k |βk| which converges under the same condition as above.
5.3 Choice of vacua One possible vacuum is the instantaneous lowest energy vacuum: minimize the energy at a 2 given instant of time η0. A solution exists, if ωk(η0) > 0. The initial conditions for the mode functions in this case are6: 1 vk(η0) = p , ωk(η0) 0 p vk(η0) = i ωk(η0),
which turn out to be isotropic. The Hamiltonian is:
Z 1 Hˆ (η ) = d3k ω (η ) aˆ+aˆ− + δ(3)(0) . 0 k 0 k k 2
|0η0 i is the vacuum of instantaneous diagonalization, since the Hamiltonian is diagonal in its exited states. At later times, the vacuum usually turns into a superposition of excited states. In general, there is no unique prescription to define a vacuum. Several problems may occur:
• It makes sense to talk about a particle with momentum p, if ∆p p and hence the 1 spatial size λ p . But in this region, the spatial curvature may change. For high wave numbers and small curvature, this effect might become irrelevant, even on a cosmological scale.
• In general the 4d curvature counts, not only the spatial curvature. For example for a q 2 2 FLRW universe, the frequency ωk(η) = k + meff(η) can become imaginary. In this case, the modes become growing and decaying. No lowest energy state exists.
6 Compare also with the free field solution (5.2).
8 • Further the definition depends on the chosen reference frame.
• Even for small deviations in the curvature and on a small time interval, the particle
production from |0η1 i to |0η2 i can be infinite. Another possible vacuum choice which avoids the last difficulty is the adiabatic vacuum. If the variation of ωk is only small in some range of η one can use a WKB approximation. In this adiabatic regime, at fixed η0 one solves for the mode function of the adiabatic vacuum, 0 s.t. vk(η0) and vk(η0) coincide with the WKB approximation Z η WKB 1 vk (η) = p exp i dη˜ ωk(˜η) ωk(η) η0
at η0.
6 Particle creation
In the Heisenberg picture, the states are time-independant, while the operator change in time. ± Let us prepare our system in the state |0ai, the vacuum corresponding to operators aˆk . These operators can be chosen in a natural way, if the universe is static at initial time (see section 2). After an interaction with the gravitational field, let us assume, that our universe is again ˆ± in a static state. But now the natural operators are bk , with vacuum |0bi. Since our universe is still in the state |0ai, but particles are counted with Nˆb, the mean density of b-particles is Z 3 2 d k |βk| ,
as we have seen above. The energy density is Z 3 2 d k ωk |βk| .
6.1 Example Consider a massless scalar field in a FLRW universe which is static in the past and future. We want to solve
00 2 χk − ∆χk + meff χk = 0, 1 m2 = a2 ξ − R, eff 6 with initial condition (in-state) e−iωη χ(in)(η) = √ . k 2ω We can reformulate the problem as an integral equation Z η (in) 1 χk(η) = χk (η) + dη˜ χk(˜η) meff(˜η) sin ω(η − η˜), ω −∞ which can be solved perturbatively by iteration. One application is particle creation due to gravitational interactions during inflation. After a
9 deSitter expansion (p = −ρ, a ∼ eHt), the universe enters a matter dominated phase (p = 0). The energy contribution can be estimated by
2 ρVac ρ ' (1 − 6ξ) ρVac . ρPlanck
Here ρVac is the density driving inflation and ρPl is the Planckian density. For GUT scale inflation, this contribution is negligible, if reheating is efficient: