UNIVERSITÀ DEGLI STUDI DI PADOVA
DIPARTIMENTO DI FISICA “G. GALILEI”
Scuola di dottorato di ricerca in fisica, ciclo XXII
Cosmological correlation functions in scalar and vector inflationary models
Supervisori: Dottoranda:
Prof. Sabino Matarrese Emanuela Dimastrogiovanni
Dr. Nicola Bartolo
Contents:
 1 Introduction
 2 SchwingerKeldysh formalism
 3 Scalar loop corrections to
 4 Perturbative expansion of the Lagrangian in theories
 5 Tensor loop corrections to
 6 Complete expression for at one loop
 7 Study of perturbations in anisotropic cosmologies
 8 Inflation and primordial vector fields
 9 Correlation functions of in the model
 10 Amplitude of nonGaussianity: and
 11 Shape of nonGaussianity and statistical anisotropy features
 12 Trispectrum for models of gauge interactions
 13 Conclusions
 14 Acknowledgments

15 Appendices
 15.1 Computation of the interaction Hamiltonian
 15.2 Study of leading slow roll order vertices in the fourth order action
 15.3 Complete expressions of oneloop twovertex diagrams to leading order
 15.4 Background and first order perturbation equations for the gauge fields
 15.5 Calculation of the number of efoldings of single(scalar)field driven inflation in the presence of a vector multiplet
 15.6 Complete expressions for the functions appearing in the bispectrum from quartic interactions
 15.7 More details on computing the analytic expressions of vectorexchange diagrams
 15.8 Complete expressions for functions appearing in pointinteration diagrams
 16 References
to Matteo
1 Introduction
In the standard cosmological model, at very early times the Universe undergoes a quasi de Sitter exponential expansion driven by a scalar field, the inflaton, with an almost flat potential. The quantum fluctuations of this field are thought to be at the origin of both the Large Scale Structures and the Cosmic Microwave Background (CMB) fluctuations that we are able to observe at the present epoch [1]. CMB measurements indicate that the primordial density fluctuations are of order , have an almost scaleinvariant power spectrum and are fairly consistent with Gaussianity and statistical isotropy [2, 3, 4, 5, 6, 7, 8, 9]. All of these features find a convincing explanation within the inflationary paradigm. Nevertheless, deviations from the basic single(scalar)field slowroll model of inflation are allowed by experimental data. On one hand, it is then important to search for observational signatures that can help discriminate among all the possible scenarios; on the other hand, it is important to understand what the theoretical predictions are in this respect for the different models.
NonGaussianity and statistical anisotropy are two powerful signatures. A random field is defined “Gaussian” if it is entirely described by its twopoint function, higher order connected correlators being equal to zero. Primordial nonGaussianity [10, 11] is theoretically predicted by inflation: it arises from the interactions of the inflaton with gravity and from selfinteractions. However, it is observably too small in the singlefield slowroll scenario [12, 13, 152]. Alternatives to the latter have been proposed that predict higher levels of nonGaussianity such as multifield scenarios [15, 16, 17, 18, 19, 20, 21], curvaton models [22, 23, 24, 25, 26, 27] and models with noncanonical Lagrangians [28, 29, 30, 31, 32]. Many efforts have been directed to the study of higher order (three and fourpoint) cosmological correlators in these models [33, 34, 35, 36, 37, 31, 152, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48] and towards improving the prediction for the twopoint function, through quantum loop calculations [49, 50, 51, 13, 52, 53, 54, 55, 56]. From WMAP, the bounds on the bispectrum amplitude are given by [8] and by [9] at CL, respectively in the local and in the equilateral configurations. For the trispectrum, WMAP provides [57] ( is the “local” trispectrum amplitude from cubic contributions), whereas from LargeScaleStructures data [58], at CL. Planck [59] is expected to set further bounds on primordial nonGaussianity.
Statistical isotropy has always been considered one of the key features of the CMB fluctuations. The appearance of some “anomalies” [60, 61, 62] in the observations though, after numerous and careful data analysis, suggests a possible a breaking of this symmetry that might have occurred at some point of the Universe history, possibly at very early times. This encouraged a series of attempts to model this event, preferably by incorporating it in theories of inflation. Let us shortly describe the above mentioned “anomalies”. First of all, the large scale CMB quadrupole appears to be “too low” and the octupole “too planar”; in addition to that, there seems to exist a preferred direction along which quadrupole and octupole are aligned [63, 64, 60, 65, 66]. Also, a “cold spot”, i.e. a region of suppressed power, has been observed in the southern Galactic sky [61, 67]. Finally, an indication of asymmetry in the largescale power spectrum and in higherorder correlation functions between the northern and the southern ecliptic hemispheres was found [68, 62, 69]. Possible explanations for these anomalies have been suggested such as improper foreground subtraction, WMAP systematics, statistical flukes; the possibilities of topological or cosmological origins for them have been proposed as well. Moreover, considering a power spectrum anisotropy due to the existence of a preferred spatial direction and parametrized by a function as
(1.1) 
the fiveyear WMAP temperature data have been analyzed in order to find out what the magnitude and orientation of such an anisotropy could be. The magnitude has been found to be and the orientation aligned nearly along the ecliptic poles [70]. Similar results have been found in [71], where it is pointed out that the origin of such a signal is compatible with beam asymmetries (uncorrected in the maps) which should therefore be investigated before we can find out what the actual limits on the primordial are.
Several fairly recent works have taken the direction of analysing the consequences, in terms of dynamics of the Universe and of cosmological fluctuations, of an anisotropic preinflationary or inflationary era. A cosmic nohair conjecture exists according to which the presence of a cosmological constant at early times is expected to dilute any form of initial anisotropy [72]. This conjecture has been proven to be true for many (all Bianchi type cosmologies except for the the Bianchi typeIX, for which some restrictions are needed to ensure the applicability of the theorem), but not all kinds of metrics and counterexamples exist in the literature [73, 74, 75]. Moreover, even in the event isotropization should occur, there is a chance that signatures from anisotropic inflation or from an anisotropic preinflationary era might still be visible today [76, 77, 78, 79]. In the same context of searching for models of the early Universe that might produce some anisotropy signatures at late time, new theories have been proposed such as spinor models [80, 81, 82, 83], higher pforms [84, 85, 86, 87, 88, 89] and primordial vector field models.
Within vector field models, higher order correlators had been computed in [90, 91, 92, 93, 94] and, more recently, in [95, 96] for vector fields. We considered vector field models in [97, 98]. NonAbelian theories offer a richer amount of predictions compared to the Abelian case. Indeed, self interactions provide extra contributions to the bispectrum and trispectrum of curvature fluctuations that are naturally absent in the Abelian case. We verified that these extra contributions can be equally important in a large subset of the parameter space of the theory and, in some case, can even become the dominant ones.
The promising perspective of achieving more and more precise measurements for the cosmological observables thanks to Planck and future experiments and the search for signatures that may help identify the correct inflationary model, have also motivated studies of higher order corrections to cosmological correlation functions and to the power spectrum in particular. Indeed, loop corrections to the correlators arise from the interactions involving the fields during inflation and therefore carry some important information about the physics of the very early Universe.
Loop corrections may lead to interesting effects which scale like the power of the number of efolds between horizon exit of a given mode and the end of
inflation [99, 100, 101, 102, 103]. The interest in loop corrections to the correlators of cosmological perturbations generated during an early epoch of inflation has been recently stimulated by two papers of Weinberg [104, 105]. The reason is that oneloop corrections to the power spectrum of the
curvature perturbation seem to show some infrared divergences which scale like , where is some infrared comoving momentum cutoff [106, 107, 108, 109]. However, it has been discussed in [110, 111] (see also [112]) that such potentially large corrections do not appear in quantities that are directly observable.
As to the power spectrum of curvature perturbations, oneloop corrections have been computed in singlefield slowroll inflation by D. Seery [54, 55] and by N. Bartolo and myself [56], in singlefield slowroll inflation. In [56] we completed the analysis carried out in [54, 55], where the metric tensor fluctuations had been neglected for simplicity, by including them in the calculations and proving that their contribution is as important as the one from the scalar perturbations. In the context of loopcalculations, we have also been working on corrections to the power spectrum in theories with noncanonical Lagrangians, which allow for higher and possibly observable corrections [113].
It can be safely stated that in standard singlefield slowroll inflation, the perturbative expansion is wellbehaved, in the sense that the agreement with observations found at treelevel for the powerspectrum is not spoiled by the radiative corrections and, on a more general basis, higher order loop corrections introduce smaller and smaller corrections as the perturbation series expansion progresses. This is not generically true in more general theories, such as for instance models with noncanonical Lagrangians, for which bounds need to be requested on the parameters of the theory in order to preserve the validity of the perturbative approach [114, 115].
This thesis collects the main results of our work on loop corrections to the power spectrum in theories of scalar inflation (Secs. 2 to 6) [56, 113], on anisotropic preinflationary cosmologies (Sec. 7) [78] and on primordial nonGaussianity and anisotropy predictions from theories of inflation where vector fields can play a role in the production of the late time cosmological fluctuations (Secs. 8 to 12) [97, 98]. The N and the SchwingerKeldysh formalisms are some of the main tools of our computation and will be briefly reviewed.
2 SchwingerKeldysh formalism
The temperature fluctuations in the CMB are rather small, of order . Theoretical predictions for the powerspectrum of curvature perturbations during inflation provide a very good match at tree level: this suggests that it is correct to use perturbation theory to evaluate cosmological correlators. A formalism conveniently employed to implement the perturbative approach is the SchwingerKeldysh, also dubbed as “inin”, formalism. It was first formulated in [49, 50, 51], later applied by J. Maldacena in [13] to the calculation of the bispectrum of curvature fluctuations and revived by S. Weinberg in [52, 53]. In this formalism the expectation value of a field operator is given by
(2.1) 
where represents the vacuum of the interacting theory, and are timeordering and antiordering operators, the subscript indicates the fields in interaction picture and is the interaction Hamiltonian. The interaction picture has the advantage of allowing to deal with free fields only; the fields can be thus Fourier expanded in terms of quantum creation and annihilation operators
with commutation rules
The inin formula has many similarities with the Smatrix in quantum field theory in terms of mathematical structure and perturbative approach, but they also have fundamental differences: the Smatrix corresponds to a transition amplitude between an initial and a final state; a cosmological correlation function is instead the expectation value of a given observable at a given time; moreover, asymptotic states in cosmology are only defined at very early times, when the same initial conditions as in Minkowsky spacetime apply for the free fields.
Using the positive and negative path technique of the inin formalism [52, 53], the expectation value above can be recast in the form
(2.2) 
where the plus and minus signs indicate modified Feynman propagators, i.e. modified rules of contraction between interacting fields; schematically we have
(2.3) 
where the sum is taken over all of the possible sets of field contractions and
In momentum space we have
It is important to remember that, when we apply this formalism, the external fields are always supposed to be treated like fields.
3 Scalar loop corrections to
The power spectrum for the comoving curvature perturbation is defined by
(3.1) 
This and all other correlation functions presented in this thesis are computed using the formula. at a given time can be interpreted as a geometrical quantity indicating the fluctuations in the local expansion of the universe; in fact, if is the number of efoldings of expansion evaluated between times and , where the initial hypersurface is chosen to be flat and the final one is uniform density, we have
(3.2) 
The number of efoldings depends on all the fields and their perturbations on the initial slice. In principle, since the fields are governed by second order differential equations, it should also depend on their first time derivatives, but if we assume that slowroll conditions apply, then the time derivatives will not count as independent quantities.
Let us apply Eq. (3.2) to the computation of in singlefield slowroll inflation (the Lagrangian for the scalar field is given by )
(3.3) 
The sums can be expanded to the desired order. Up to one loop we have
(3.4)  
where a star indicates evaluation around the time of horizon crossing. Eq. (3.4) can finally be rewritten as [116, 109]
(3.5)  
is the tree level power spectrum (3.6)
(3.6) 
where is the Hubble parameter evaluated at horizon exit (when ). The variance per logarithmic interval in is given by . The one loop contribution to the power spectrum is given by
(3.7) 
where the first term on the righthand side, , accounts for the contributions coming from the inflaton selfinteractions and were computed by D. Seery in [108, 109]
(3.8) 
where and are numerical factors. Their diagrammatic representation is given in Fig. for the leading order and in Fig. for the nexttoleading order corrections. The loop corrections , arising from interactions between the tensor (graviton) modes and the scalar field, were ignored for simplicity in [108, 106], however they should be included since they are not slowroll suppressed compared to loops of scalar modes. Their computation was presented for the first time in our paper [56] and will be reviewed in Secs. 4 to 6 of this thesis.
Both and are evaluated at around the time of horizon crossing and as such they are due to genuine quantum effects.
The contributions in the third and fourth lines of Eq. (3.5), also dubbed as “classical oneloop”, can be considered as classical loop contributions arising after the perturbation modes leave the horizon. The distinction between classical and quantum loops is intended as for example in [109]: quantum loops find their origin in the Lagrangian interaction terms between the inflaton perturbations and the gravitational modes or from selfinteraction of ; classical loops are corrections merely coming from the expansion of using the formula and originate from zeroth order terms in the SchwingerKeldysh formula.
4 Perturbative expansion of the Lagrangian in theories
In this and in the next two sections, we will review the computation of the tensor loop corrections to . For our purposes, the exponentials in Eq. (2.2) need to be expanded up to second order in the interaction Hamiltonian
(4.1)  
where . Oneloop powerspectrum diagrams require an expansion of the interaction Hamiltonian to third and fourth order in the field fluctuations, i.e. . We provide in Figs. (3) and (4) the diagrammatic representation of the leading order corrections that we will find for the diagrams with tensor loops in singlefield slowroll inflation. The continuos lines represent scalar propagators, whereas the dotted lines indicate tensor propagators. In order to derive this result and the analytic expressions for these diagrams, we need to first calculate and expand up to fourth order in the field perturbations and . The starting point is the Lagrangian of the system.
We will begin with a more general Lagrangian for the scalar field than the usual , by introducing a nonconventional kinetic term, i.e.
(4.2) 
where is a generic function of the scalar field and of and is the Ricci scalar in four dimensions. Notice that the action (4.2) reduces to the standard case if , where is the potential for the scalar field.
Theories of inflation where the Lagrangian kinetic term is a generic function of the scalar field and its first derivatives, like in Eq. (4.2), are string theoryinspired. They represent interesting alternatives to the basic inflationary scenario because of their nonGaussianity predictions. The crucial quantity in this sense is represented by the speed of sound , which is allowed to vary between and . The perturbative expansion of the interaction Hamiltonian in this kind of models has coefficients proportional to inverse powers of the sound speed and therefore, for small values of , allows both for nonnegligible loop corrections to the power spectrum of the curvature fluctuations [113] and for large values for the amplitudes of three [31] and four [41, 42, 43, 44, 45, 46] point functions. In this thesis, we will carry out the calculations of the interaction Hamiltonian for these general theories up to a certain point and then, for simplicity in the presentation, focus on the canonical case (the remaining computations for more general Lagrangians will be found in [113]).
Let us list the background equations for the system
(4.3)  
(4.4)  
(4.5)  
where a dot indicates a derivative w.r.t. cosmic time and, to zeroth order, we have .
The so called flowparameters are defined as
(4.6)  
(4.7) 
These quantities reduce to the slowroll parameters in the standard case, so it is natural to assume and . It is not correct to talk about slowroll if is left as a generic function of and , since the smallness of and does not necessarily indicate that and . It can be convenient to decompose as the sum , where
(4.8)  
(4.9) 
The parameters that are expected to appear in the perturbative expansion of the Lagrangian are
(4.10)  
(4.11)  
(4.12)  
(4.13)  
(4.14)  
(4.15) 
where is the sound speed. is allowed to vary between and , so the quantity can freely range between and . The only assumption we make is , from being constant in the standard case.
4.1 ArnowittDeserMisner (ADM) decomposition for theories
The Lagrangian in Eq. (4.2) will now undergo a perturbative expansion in terms of the field fluctuations ( is the homogeneous background value for the field) and of the metric fluctuations.
It is convenient to adopt the ArnowittDeserMisner (ADM) splitting for the metric. In the spatially flat gauge the perturbed metric is
(4.16) 
(4.17) 
where is the scale factor, is a tensor perturbation with (traceless and divergenceless) and
det. Notice that repeated lower indices are summed up with a Kronecker delta, so stands for and .
In the ADM formalism, the action (4.2) becomes [13]
(4.18) 
where is the curvature scalar associated with the three dimensional metric and
A dot indicates derivatives w.r.t. time , all the spatial indices are raised and lowered with and units of will be from now on employed. To order we have
(4.19) 
The lapse and shift functions, and , can be written as
where , and are functions of time and space ( is divergenceless). We have exploited the gauge freedom to set two scalar and two vector modes to zero, thus leaving one scalar mode from , one scalar and two vector modes from and two tensor modes (the two independent polarizations of the graviton) from , together with the inflaton field perturbation . and are nondynamical degrees of freedom and can be expressed in terms of the other modes ( and ), once the Hamiltonian and the momentum constraints (we derive them in the next section) are solved.
Solving Hamiltonian and momentum constraint equations
Momentum and Hamiltonian constraints are derived from varying the action w.r.t. the shift and lapse functions respectively. It turns out that, in order to expand the action to a given order , it is only necessary to perturb and up to order [13, 31]. Therefore we will solve the constraints to second order in the metric and scalar field fluctuations.
Let us begin with the expansions
where and are respectively first and second order in the fields fluctuations (similarly for and , and for and ). Let us then expand to second order. is a generic function of and . We first need the expansion of
(4.20)  
where , is the zeroth order part, i.e. and is the perturbation to the desired order (). Notice that , but for simplicity we will suppress the subscript ’’ in the background value of the field.
The expressions for the perturbations become
(4.21)  
(4.22) 
and so on for and higher order terms. The expansion of up to second order becomes
where as usual the subscript ’’ indicated the zeroth order, , and similarly for the second order derivatives, and needs to be expanded up to the needed order.
We are now ready to write the momentum and Hamiltonian contraints
(4.24) 
(4.25) 
The momentum constraint to first order reads
(4.26) 
where is the Hubble parameter. Eq. (4.26) can be solved to derive . Taking the derivative of both sides of (4.26) and using the divergenceless condition for , we have
(4.27) 
Using the solution found for , we find , from which we can set . Here ,
which we will indicate in the rest of the thesis also as , and from now on we define for simplicity.
The momentum constraint to second order is
(4.28) 
The solutions are
(4.29)  
where , and
(4.30)  
Let us now move to the Hamiltonian constraint which provides
to first order and