Power spectrum oscillations from Planck-suppressed operators in effective field theory motivated monodromy inflation

Power spectrum oscillations from Planck-suppressed operators in effective field theory motivated monodromy inflation


We consider a phenomenological model of inflation where the inflaton is the phase of a complex scalar field . Planck-suppressed operators of modify the geometry of the vev at first order in the decay constant , which adds a first order periodic term to the definition of the canonically normalized inflaton . This correction to the inflaton induces a fixed number of extra oscillatory terms in the potential . We derive the same result in a toy scenario where the vacuum is an ellipse with an arbitrarily large eccentricity. These extra oscillations change the form of the power spectrum as a function of scale and provide a possible mechanism for differentiating EFT-motivated inflation from models where the angular shift symmetry is a gauge symmetry.

I Introduction

Single-field, slow-roll inflation models that generate detectable levels of primordial gravitational waves require a super-Planckian field displacement for the inflaton (1). The constraints on the primordial power spectrum from Planck (2); (3); (4) remain broadly consistent with the predictions of large classes of these models, particularly those with for . Knowledge of the detailed predictions of these models in the range has become increasingly important as estimates of the tensor-to-scalar ratio have improved (5); (6); (7).

However, Planck-scale effects should substantially alter the simple polynomial relationship over regions of field space that exceed . Monodromy is a mechanism that protects against these dangerous corrections by repeatedly wrapping the trajectory of the inflaton around a sub-Planckian region of a higher dimensional space, with an angular parameter . The mechanism relies on a corresponding change for each winding , with the global potential defined by patching together branch cuts. Large field excursions are then safe because the displacement in the “natural” parameter space remains small.

Importantly, monodromy inflation predicts detectable oscillatory features in the primordial power spectrum of curvature perturbations  (8); (9). If the angular field is coupled to gauge fields, then non-perturbative instanton effects modify the potential . In the absence of monodromy, this oscillatory contribution can drive natural inflation (10). Combined with a monodromy effect in , this contribution imprints oscillations in the power spectrum with a frequency that is related to the periodicity of the underlying theory through the axion decay constant . While the size of these oscillations is constrained by the WMAP (11) and Planck data (2); (12); (13); (14); (4), their detection would provide an unprecedented glimpse into the monodromy structure of the high energy theory. Therefore, it is important to understand the exact nature of these oscillations and to identify any mechanisms that might change this prediction.

In this paper we study an alternative source of oscillations in the primordial spectra in monodromy inflation, where the angular degree of freedom is the phase of a complex scalar with a vacuum expectation value (vev) of . Adding non-renormalizable terms to the action for alters the shape of the vacuum so that and the effective decay constant obtains an angular dependence. These higher order terms only contribute to the potential as , which is much smaller than unity for sub-Planckian decay constants . However, they change the form of the canonical kinetic energy term for in the vacuum at first order in . Consequently, the canonically normalized inflaton , which depends on the structure of the kinetic energy term for , obtains a quasi-periodic dependence at , where is the energy dimension of the non-renormalizable operator.

Therefore, the leading order gravitational operators with add first order oscillatory corrections to any term in the scalar action that depends on the angular coordinate . If the shift symmetry in the vacuum is broken, then the non-renormalizable operators add a set of sinusoidal functions to that have arbitrary relative phases. These new terms are independent of the scale of any instanton effects derived from gauge couplings and result in a superposition of high frequency power spectrum oscillations. These features must be taken into account when fitting the predictions of monodromy inflation to data.

Additionally, we consider a case where the structure of the vacuum is shaped like an ellipse with an arbitrarily large eccentricity. The coupling constant for the oscillatory terms are suppressed by a factor of at least compared to the angular symmetry breaking scale, even for a highly elliptical vev. A flat prior probability on the eccentricity , which corresponds to a naturally large correction to the vev, gives sinusoidal terms in that are suppressed at the order of compared to the contribution from the chaotic potential . This gives power spectrum features that are still safely inside the current bounds from Planck data (12).

Adding Planck-suppressed terms to the Lagrangian is motivated by an effective description of monodromy inflation at the energy scale . However, we do require the existence of global symmetries in the low-energy theory, which is a non-trivial assumption that makes obtaining large-field inflation simple. A more rigorous description of the low-energy theory might require significant fine-tuning in order to obtain a flat potential across super-Planckian field ranges. Consequently, we assume that an angular symmetry breaking term is technically natural in a field theory that contains a global angular shift symmetry (15) and not include operators that would spoil the chaotic inflation behavior of the potential . We take a phenomenological approach and remain agnostic about a specific high energy scenario that realizes this symmetry.

Monodromy arising from a gauge symmetry is a ubiquitous feature of high energy theories. In particular, string theory has many axions, which are angular parameters that are given a monodromy term by couplings to fluxes (16); (8); (17); (9); (18); (19); (20); (21). The monodromy potential for string axions can have a variety of exponents, with giving predictions that match the Planck posterior probabilities for the model within approximately - (21). While the primordial spectrum resulting from monodromy inflation with a broken global shift symmetry has an oscillatory contribution from each Planck-suppressed operator, the simplest models with a gauged axion symmetry give a coherent, sinusoidal signal in . In principle, these two signals can be distinguished via cosmic microwave background or large-scale structure data.

Furthermore, although a large number of angular parameters could contribute to inflation in principle, we only consider the single-field case. While many-field models can have simple predictions (22); (23); (24), the variance of the predicted values of observables for monodromy inflation can be undefined for  (25). Since small exponents are observationally favored, we will leave an analysis of many-field monodromy to future work.

This paper is organized as the following: Sect. II discusses how non-renormalizable operators alter the shape of the vacuum for a complex scalar with a Mexican hat type potential. Sect. III calculates the canonically normalized inflaton for an arbitrarily elliptical vev. Sect. IV estimates the shape of the features in from this mechanism both analytically and numerically. Sect. V compares the predicted spectra from Planck-corrected, EFT-motivated monodromy models to high energy axion monodromy scenarios. Sect. VI summarizes the findings.

Ii Modifying the geometry of the vacuum with irrelevant operators

We start by considering a Lagrangian for a complex scalar field,


where the renormalizable terms in the action satisfy the canonical Mexican hat shape with


We assume that the dimensionless real parameter is and that the decay constant is . The vev is independent of the phase and is sub-Planckian for small decay constants. In the vacuum, the radial component of has the mass . If is coupled to gauge fields, then non-perturbative instanton effects induce a sinusoidal correction to Eq. (2), which gives rise to natural inflation if  (10).

However, monodromy inflation can be generated in the more natural limit by explicitly breaking the global symmetry. We do this by adding a term to Eq. (2) like


where the subscript ASB denotes angular symmetry breaking. We assume that the form of this term remains stable over regions , which is satisfied in the vacuum. The domain of can be analytically extended to by piecing together each of the branches of the Mexican hat, which have a change in potential energy of for every winding. This symmetry breaking mechanism has also been invoked as a means to dynamically lower the value of the cosmological constant (26) and is well motivated from effective field theory independently of the monodromy mechanism we study here.

With Eq. (3) large-field inflation can be realized by the canonically normalized, pseudo–Nambu-Goldstone boson when , if the initial angular displacement is sufficiently large . For this gives an effective inflaton mass of , which is less than the Hubble scale independently of and is sub-Planckian if . To match the amplitude of scalar perturbations as estimated from Planck (7), the symmetry breaking parameter should be approximately , although the exact value depends on and .

To ensure that the mass of the radial direction is when the pivot scale leaves the horizon approximately -folds before the end of inflation we further require


where we have removed terms in the second approximation. A super-Hubble radial mass keeps the evolution close to the vacuum during inflation, which is an important requirement for this paper.1

For a concrete example, Refs (28); (29); (30) obtain monodromy in the phase of a complex scalar with a supergravity model where the global symmetry of the Kähler potential is broken by the superpotential in such a way that it can give a variety of inflationary potentials , including Eq. (3). Refs (31); (32); (33); (34) have similar helical potentials with a range of phenomenological predictions. Although realistic scenarios must account for analytical control over the symmetry breaking term at high energies, as well as the stability of a global minimum for reheating, we do not worry about these considerations here.

In addition to the Mexican hat shape of Eq (2) we can also add any operator that has the form


where is a dimensionless complex coupling, is the cutoff scale, and


The operator has energy dimension


and breaks the translational symmetry in Eq. (2), leaving a residual discrete symmetry . All operators of this type are generally expected from effective field theory considerations at energy dimension . With the cutoff scale set to , we interpret these operators as gravitational corrections to the renormalizable theory at sub-Planckian energies.2

The impact of operators like Eq. (6) has been well studied in the context of the Peccei-Quinn axion (35), where is coupled to the QCD gluon field and setting the vacuum dynamically solves the strong-CP problem. Experimental constraints on CP violation require and non-renormalizable terms like Eq. (6) force , giving a generally CP-violating vacuum in the absence of extreme fine-tuning in the relative phases of the operators (36); (37); (38); (39); (40).

Figure 1: Value of when corrected by Eq. (8) in comparison to . The parameters are , , and . We show the cases (red) , ; (blue) , ; and (yellow) , .

Typically, for inflation driven by the slowly rolling field , an operator is negligible since it is at most an correction to in the vacuum . However, in the absence of Planck-suppressed operators like Eq. (5), the local minimum is stabilized by the term in the Mexican hat potential (2). Consequently, the operators alter the minimum of in the radial direction of at first order in . Operators of higher energy dimension change the location of the vacuum as


where is the phase of the coupling and we have ignored terms and terms higher order in . Fig. 1 shows for different values of and , with the decay constant set to a large value that has been exaggerated in order to visualize the change in shape of the vev.

If we assume that , in addition to the negligible contributions from , the deviation of the vacuum geometry from in Eq. (8) affects the kinetic energy of the complex scalar as


We define a scalar field with a canonical kinetic energy in the usual fashion:


which is valid up to an arbitrary constant.

For simplicity we focus on a single leading order non-renormalizable operator with energy dimension . The derivative term in Eq. (10) for this operator is and can be effectively ignored compared to the contributions from . Up to the integrand reduces to


where we have defined the dimensionless constant


which we expect to be although it could be larger if . Assuming that the scale of does not disrupt the perturbative analysis with respect to , the inflaton can be approximated as


which we can invert order-by-order in to get


In the limit that we recover the standard definition .

Consequently, in the low-energy description we can fix and replace all instances of in the scalar sector of the action by the relationship given in Eq. (14). Defining the first and second terms in Eq. (14) as , the potential can be expanded like


For the combined Mexican hat action (2) and angular symmetry breaking potential (3), the effective potential at tree level with the Planck-corrected inflaton in Eq. (13) becomes


plus terms higher order in . Fig. 2 shows this potential for . Since it is the most conservative choice we use a Planckian cutoff scale for the rest of the paper.

The second term in Eq. (16) does not come directly from the Planck-scale operators (6), but rather from the substitution into Eq. (3), which is due to the modified shape of the vacuum in Eq. (8). The phase is important, since any of the possible oscillatory terms from Planck-suppressed operators and any oscillations coming from instanton effects generally have different relative phases, which we discuss in Sect. V.2. With and , the ratio fixes both the scale and frequency of the new oscillatory term compared to the angular symmetry breaking scale.

Figure 2: The portion corresponding to of the modified monodromy potential with one additional Planck-suppressed operator of energy dimension , with and .

Iii A non-perturbative elliptical vacuum

For sub-Planckian field values the corrections to the shape of the vacuum from non-renormalizable operators like Eq. (6) are naturally small. However, we show in this section that the corrections to the potential can be kept perturbatively small even when the vev is altered non-perturbatively.

For a simple example with a non-circular vev we consider the case where is an ellipse with semi-major axis and eccentricity , which are related by the semi-minor axis as


Instead of the Mexican hat potential (2), we assume that is structured so that keeps the potential energy constant and setting restores the symmetry. We leave the exact form of this potential implicit since we will only worry about the dynamics when . In analogy to the perturbative corrections to the vacuum in Sect. II, we keep a canonical kinetic energy term for , which breaks the -translational symmetry for the Lagrangian when .

We parametrize the complex scalar in terms of a radial coordinate and an angular coordinate as


Surfaces of constant are ellipses with foci at and we define the vev so that it follows one of these ellipses,


We assume that the effective mass in a direction locally orthogonal to is sufficiently high so that we can treat the system as if it were one dimensional at low energies, with non-dynamical. Near the vacuum we can therefore fix the location of the foci at .

Figure 3: The canonically normalized inflaton as a function of the elliptical angular variable , where the vev has semi-major axis and eccentricity . The dependence is invertible but non-elementary.

The canonical kinetic energy term for in the vacuum becomes


As in Eq. (10) we define the canonical inflaton as


which can be expressed in terms of the incomplete elliptic integral of the second kind as


Fig. 3 shows this relationship for a range of eccentricities . Although the function is bijective, its inverse is not elementary due to the special function . However, we can use the properties of to obtain an approximation to and calculate the potential energy of the canonical inflaton .

For , is quasi-periodic with period , satisfying


and where is an integer and is an arbitrary parameter. This relationship ensures that is composed of a linear and an oscillatory component, which Fig. 3 shows are approximately independent even as . Motivated by this, we decompose in Eq. (22) into a linear and oscillatory portion as


where Eqs (22) and (23) give the pre-factor


This coefficient is bounded in the range for . To determine the form of the oscillatory term in Eq. (24) we examine the limiting behavior of the elliptic integral.

For small eccentricities we can expand the elliptical integrals , which appear in and in Eq. (22), to get


We can compare this directly to the perturbative case of Eq. (13) and see that the functional form of these expressions matches in this limit, as expected. Similarly, for large eccentricities , Eq. (22) becomes


Therefore, Eq. (27) has an identical functional form for the oscillatory term as the small-eccentricity expression Eq. (26) at first order in . Since Fig. 3 shows that intermediate eccentricities have similar behavior, a close approximation to Eq. (22) for all eccentricities is


where we have defined the coefficient


in terms of the change in the elliptical integral, which is


The coefficient for the oscillatory term is bounded in the range for eccentricities .

Eq. (28) is valid for any ellipse and is not restricted to a perturbative expansion in . Although this expression could be inverted numerically to obtain , the inverse expression is transcendental. An analytical inversion, as in Eq. (14), can only be made once a suitably small parameter is identified. Rather than restricting our analysis to small eccentricities, which would replicate the results of Sect. II, we instead solve for perturbatively in the parameter


This ratio is naturally small since the oscillatory term in Eq. (28) is sub-dominant to the linear even in the large-eccentricity limit (27). Importantly, satisfies the exact bounds


The separation in scales between the linear and oscillatory contributions to justifies an order-by-order perturbative solution to the inverse function . At first order this yields


where the effective decay constant is


which is fixed for a given eccentricity. Eq. (33) is analogous to Eq. (14) in Sect. II. Fig. 4 compares this perturbative inverse function to the exact numerical inverse in the maximum eccentricity case . In this limit the ratio maximizes and the approximation is close to the numerical values even for this most extreme case. Fig. 4 only shows , but the periodicity of the elliptical integral in Eq. (23) ensures this approximation is valid for all values of .

Figure 4: The inflaton in terms of the angular variable when the expression (28) is inverted both numerically and approximately with Eq. (33). This shows the worst-case scenario with maximum eccentricity , with good agreement between the numeric and approximate solutions.

Given the relationship (33) we can replicate the analysis of Sect. II, substituting into the low-energy effective potential. In particular, since we constructed this potential so that it has the symmetry , we assume that the angular symmetry breaking mechanism acts as


rather than the expression in Eq. (3). Therefore, the elliptical vev yields a correction to the potential as


to first order in . This expression closely matches Eq. (16), which arises from modifying the geometry of the vev from adding non-renormalizable operators to the Lagrangian. However, the frequency of the oscillation term in Eq. (36) is not possible with an operator of energy dimension . Therefore the oscillatory term above is not Planck-suppressed. Instead, large elliptical deviations from a circular vev should come from relevant or marginal operator. We leave the exact construction of this type of scenario ambiguous.

This toy model demonstrates that even non-perturbatively modified vacuum geometries can still yield easily categorizable alterations to the basic shape of the potential. As long as we can identify an angular parameter that varies over some finite range in the absence of monodromy, then the relationship in the modified vacuum should also be periodic and invertible, since Fourier series like Eq. (6) are the natural basis with which to build operators over a compact subset of the real numbers. However, one might expect that being able to identify perturbative parameters that give a tractable perturbative inverse function is model dependent.

Iv Power spectrum features from non-circular monodromy

Changing the shape of the vev in Sects II and III has resulted in a modified monodromy potential of the form


where the constants , , and are related to the parameters in the potentials of Eqs. (16) and (36) in the obvious way. The ratio sets the amplitude of scalar perturbations and controls the oscillatory contribution to the primordial spectra. We estimate the power spectrum resulting from Eq. (37) using the formalism. This relates the curvature perturbation on equal density hypersurfaces to the field perturbations on spatially flat hypersurfaces via


to first order in , where is the number of -folds between the flat and equal density hypersurfaces and .

We approximate the power spectrum whenever a mode crosses the horizon by ignoring any non-standard sub-horizon evolution for the Fourier transformed field values . This gives the standard result for the power spectrum of scalar perturbations


where . From Eq. (38) the primordial curvature power spectrum is


Using the slow-roll approximation, the number of -folds of expansion can be calculated as


where is the value of the field at the end of inflation and is the field value when . The derivative of Eq. (41) is


where we have disregarded terms that are . The third term in the square brackets can be ignored for , which is the case for those modes that leave the horizon within a few -folds of the pivot scale , which has for .3

We estimate the horizon crossing value of the inflaton field in Eq. (42) using only the non-oscillatory portion of Eq. (37), i.e., with . This is only a rough approximation to the slow-roll dynamics for and results in a value for in Eq. (42) that artificially suppresses the oscillatory terms, since we are ignoring a component to that is in phase with these terms. A more accurate estimate of might be obtained by replicating the techniques4 of Refs (9); (41) or using generalized slow-roll techniques (42); (43); (44). However, this approximation suffices to understand the general behavior of the power spectrum as a function of scale.

Setting reduces Eq. (41) to


Since the power law term in Eq. (37) gives large-field inflation, we ignore any contributions to that result from throughout the rest of this section. For the pivot scale, Eq. (43) gives


We relate to the field value when the pivot scale leaves the horizon -folds before the end of inflation by


where we have used the slow-roll approximation and have related the scale factors at horizon crossing to the scale factor at the end of inflation by . Using Eq. (43), this can be expressed as


Modes that leave the horizon within a few -folds of each other have only a small change in . Consequently, we solve this equation for while ignoring in order to understand the behavior of only near . The solution to Eq. (46) in this limit is


This expression is valid only for and when depends only weakly on the oscillatory terms in .

Figure 5: Primordial power spectrum from the potential (37), normalized by the amplitude of the power spectrum at the pivot scale . The parameters are , , , , and , which describes an elliptical vacuum of eccentricity . The solid orange line is the approximation and the dashed blue line is the exact numerical value.

Substituting Eqs (42) and (47) into Eq. (40) gives an estimate of the primordial power spectrum. Fig. 5 compares this approximation to an exact numerical calculation for the potential (37) using the publicly available inflation solver MultiModeCode (45); (46); (47); (48). MultiModeCode numerically solves the Klein-Gordon and mode equations for on sub-horizon scales and is exact to first order in . The chosen parameters for Fig. 5 correspond to an elliptical vev with eccentricity , which results in sinusoidal oscillations in as a function of . The analytic solution is a close approximation, although the amplitude of the oscillations is suppressed relative to the exact calculation, as expected. For smaller amplitude oscillations with the approximation becomes a better fit.

Constraints on primordial tensor modes from Planck and BICEP2 (2); (3); (4); (6) restrict , although this is model dependent. To match recent studies that constrain oscillations in the primordial spectrum (11); (49); (13); (14), we generally expect and for the fiducial case of . When Eq. (37) is generated by Planck-suppressed operators, the upper limit on corresponds to with defined as in Eq. (12).5

V Differentiating Planck-corrected and gauged oscillations

v.1 The toy case

In Sects II and III we considered an effective description of monodromy inflation where a global angular symmetry is explicitly broken and irrelevant, Planck-suppressed operators. This differs from many axion monodromy scenarios that come from explicit high energy considerations, such as string axions (51) or extranatural inflation (52), where the translational symmetry in is a gauge symmetry.

Typically, the field in high energy theories is identified with an axion, which is obtained by integrating gauge fields over cycles in compact extra dimensions. The residual gauge symmetry for this scalar severely restricts the allowed structure of the Lagrangian and precludes operators like those in Eq. (6) from perturbatively altering the shape of the vacuum once the gauge is fixed. However, the gauge symmetry in these cases requires that is coupled to gauge fields, which induces an oscillatory term to the potential via non-perturbative instanton effects, which take the form


We treat this as the prototypical potential for gauged monodromy inflation, although the model might be different in more specific cases, particularly those that incorporate a natural variability in the oscillation frequency (31); (32); (34); (41).

In principle, signals in the primordial spectrum due to sinusoidal terms like Eq. (37) should be distinguishable from the features derived exclusively from non-perturbative instanton effects like in Eq. (48). The difference is most obvious for the case where the vacuum geometry is modified by a single Planck-suppressed term in the Lagrangian and the potential is closely approximated by Eq. (16). Most obviously, the coupling constant for the sinusoidal term varies as a function of for . For modes that leave the horizon closer to the end of inflation feel a relatively larger effect from the oscillatory portion of Eq. (16) than those that leave the horizon earlier.6 The resulting signal in the power spectrum for these smaller modes may be detectable in the correlation functions of large scale structure in large-volume surveys (53); (54); (55); (56), which are able to constrain the amplitude of small modes .

The non-renormalizable oscillations also have a larger frequency than the instanton terms by an integer multiple,


where for the leading Planck-scale correction. This gives features in that are of higher frequency than features from Eq. (48) for fixed values of .

Furthermore, if we assume that all dimensionless couplings in Eq. (16) are exactly unity, then the coupling for the new sinusoidal term from as compared to the leading order polynomial term is


With the angular symmetry breaking exponent fixed, this only has the frequency as a free parameter. The relative amplitude of oscillatory features in the primordial spectrum that result from the instanton and non-circular vev effects is then directly tied to their relative frequencies. However, the amplitude of power spectrum oscillations that result from Eq. (37) is independent of those coming from Eq. (48).

v.2 Power spectrum with multiple oscillatory components

In general, the analysis in Sect. V.1 is overly simplistic, since we should allow for all possible non-renormalizable operators to contribute to the modification of the vacuum, including derivative terms that we have ignored throughout this paper. However, adding only those Planck-suppressed operators of energy dimension that are described in Eq. (6) changes the relationship in Eq. (14) to


where the sum is over the three values of that are allowed from Eq. (7). The coefficients are defined as in Eq. (12) for each of the , where the phases are assumed to be uniformly distributed. Eq. (51) is a perturbative inverse function that is valid to and requires that the sum of the individual couplings is sub-dominant to the linear term . This marginally lowers the range of decay constants for which this expression is valid in comparison to Eq. (14). Potentials with multiple oscillatory terms have been studied in the absence of monodromy in Refs (57); (58); (59); (60) and with monodromy in Refs (61); (62). However, the tight relationships between the coupling constants and frequencies in Eq. (51), combined with the limited number of Planck-suppressed operators of energy dimension , severely reduce the effective dimensionality of parameter space for this scenario in comparison to previous work.

Figure 6: Power spectrum when all lowest-order Planck-suppressed operators of are included. The angular symmetry breaking exponent is and the decay constant is . The relative phases between the different oscillatory contributions and their overall couplings are chosen independently from the distributions and . Two different selections of and are plotted.

The general expression (51) modifies the potential (37) so that it contains a summation of out-of-phase oscillation terms. This similarly affects the derivative of the number of -folds as


We again assume that the total contribution of the oscillatory terms above only negligibly affects the relationship in Eq. (47). Therefore, the power spectrum can be estimated by simply substituting from Eq. (52) into Eq. (42). The power spectrum then roughly corresponds to a power-law with a superposition of sinusoidal terms, each of which has a form as in Fig. 5. Fig. 6 displays the power spectrum for the superimposed signal with two different choices of couplings and phases . This relatively complex primordial signal differs substantially from the easily distinguishable oscillations expected from the potential (48) alone.

Vi Summary

In Sect. II we studied the low-energy effective description of monodromy inflation resulting from the evolution of the angular component of a complex scalar field . We started with the only renormalizable potential that gives a non-zero vev and assume . Motivated by effective field theory we included higher order terms with energy dimension . Although these terms are suppressed at order , we showed that they modify the geometry of the vev so that it has an angular dependence that is affected more severely: .

When we restrict , the new geometry for the vacuum adds an angular dependence to the canonical kinetic energy term for . In the absence of irrelevant operators, the canonically normalized scalar field that drives slow-roll inflation is , but acquires a sinusoidal dependence when these terms are included. This correction is first order in and the inverse function is given schematically by .

In Sect. III we looked at more radical modifications to the shape of the vacuum. We parametrized the vev as an ellipse with arbitrary eccentricity and found that the elliptical angular coordinate variable is related perturbatively to the inflaton in terms of the small parameter , which is true for all possible eccentricities. This parameter describes the ratio of linear and oscillatory contributions to the elliptic integral , which is guaranteed to be small since for any ellipse. This gives a description of the effective potential in the vacuum that is analogous to Sect. II, without requiring a perturbatively small modification to the circular Mexican hat vev.

Since the relationship between the angular parameter and the inflaton has a first order correction, a term that breaks the symmetry is also corrected at first order in . The modified potential obtains one or more sinusoidal terms, which each contribute oscillations to the primordial power spectrum . In Sect. IV we calculated the power spectrum resulting from these types of corrections and found that they are oscillations in .

Without a reason to exclude Planck-suppressed operators from the theory, the existence of multiple oscillatory terms in the potential is completely general. As discussed in Sect. V, the from many such terms should be closely approximated by the linear superposition of signals with varying phases. However, the individual signals have a close relationship between their frequencies, which is constrained to a few possible integer values due to the energy dimension of the non-renormalizable operators. Furthermore, the overall coupling constant for the terms that determine the size of oscillations in can be related to the oscillation frequency, if the dimensionless parameters in the theory are .

The EFT-motivated description of the low-energy theory of a complex scalar contrasts with those cases where monodromy inflation results from the simplest non-perturbative breaking of a gauge symmetry. In principle, the effects of the modified vacuum mechanism presented here would distinguish a high energy axion monodromy theory from monodromy arising in an effective low-energy description. Consequently, searches for non-sinusoidal features in the primordial power spectrum might eventually allow us to extract a Planckian signal from the primordial spectrum. Conversely, the absence of complicated oscillations in would constrain the allowed form of Planck-suppressed operators in monodromy inflation, adding to existing constraints on high dimension operators from local non-Gaussianity (63). However, obtaining significant evidence in favor of a non-sinusoidal feature in the primordial spectrum remains a challenging problem.

I thank Richard Easther, Dan Green, Hiranya Peiris, Eva Silverstein, Ewan Stewart, and Jonathan White for helpful discussions. I acknowledge the use of the New Zealand eScience Infrastructure (NeSI) high-performance computing facilities, which are funded jointly by NeSI’s collaborator institutions and through the Ministry of Business, Innovation & Employment’s Research Infrastructure programme [http://www.nesi.org.nz]. This work has been facilitated by the Royal Society under their International Exchanges Scheme.


  1. Ref. (27) has shown that the dynamics of the heavy radial direction can reduce the speed of sound of the inflaton fluctuations and suppress the amount of gravitational waves. In order to isolate the predictions resulting from Planck-suppressed operators we do not include these effects in this paper.
  2. While the operators in Eq. (6) are built with Fourier series, we could use an alternative basis, e.g., the parametrization has operators that are also periodic in with a similar phenomenology to Eq. (5). However, we ignore those bases that have higher order terms that do not leave a residual angular symmetry, as well as any possible derivative terms.
  3. The derivatives vanish for adiabatically evolving scalar fields.
  4. This becomes more difficult when the angular symmetry breaking term does not have the exponent .
  5. The data can allow much wider ranges for the parameters, since the posterior probabilities for the model parameters can be highly non-trivial. Furthermore, large wavelength oscillations in can be confused for models that have large (positive or negative) running (50), due to the finite range of multipole moments that can be measured at high accuracy in the cosmic microwave background.
  6. For and , the sinusoidal term dominates over the diverging pre-factor so that here.


  1. D. H. Lyth, What would we learn by detecting a gravitational wave signal in the cosmic microwave background anisotropy?, Phys.Rev.Lett. 78 (1997) 1861–1863, [hep-ph/9606387].
  2. Planck Collaboration, P. Ade et al., Planck 2013 results. XXII. Constraints on inflation, Astron.Astrophys. 571 (2014) A22, [arXiv:1303.5082].
  3. Planck Collaboration, P. Ade et al., Planck 2013 Results. XXIV. Constraints on primordial non-Gaussianity, Astron.Astrophys. 571 (2014) A24, [arXiv:1303.5084].
  4. Planck Collaboration, P. Ade et al., Planck 2015 results. XX. Constraints on inflation, arXiv:1502.02114.
  5. Planck Collaboration, P. Ade et al., Planck 2013 results. XVI. Cosmological parameters, Astron.Astrophys. 571 (2014) A16, [arXiv:1303.5076].
  6. BICEP2, Planck Collaboration, P. Ade et al., Joint Analysis of BICEP2/ and Data, Phys.Rev.Lett. 114 (2015) 101301, [arXiv:1502.00612].
  7. Planck Collaboration, P. Ade et al., Planck 2015 results. XIII. Cosmological parameters, arXiv:1502.01589.
  8. L. McAllister, E. Silverstein, and A. Westphal, Gravity Waves and Linear Inflation from Axion Monodromy, Phys.Rev. D82 (2010) 046003, [arXiv:0808.0706].
  9. R. Flauger, L. McAllister, E. Pajer, A. Westphal, and G. Xu, Oscillations in the CMB from Axion Monodromy Inflation, JCAP 1006 (2010) 009, [arXiv:0907.2916].
  10. K. Freese, J. A. Frieman, and A. V. Olinto, Natural inflation with pseudo - Nambu-Goldstone bosons, Phys.Rev.Lett. 65 (1990) 3233–3236.
  11. H. Peiris, R. Easther, and R. Flauger, Constraining Monodromy Inflation, JCAP 1309 (2013) 018, [arXiv:1303.2616].
  12. R. Easther and R. Flauger, Planck Constraints on Monodromy Inflation, JCAP 1402 (2014) 037, [arXiv:1308.3736].
  13. P. D. Meerburg, D. N. Spergel, and B. D. Wandelt, Searching for Oscillations in the Primordial Power Spectrum: Perturbative Approach (Paper I), Phys.Rev. D89 (2014) 063536, [arXiv:1308.3704].
  14. P. D. Meerburg and D. N. Spergel, Searching for Oscillations in the Primordial Power Spectrum: Constraints from Planck (Paper II), Phys.Rev. D89 (2014) 063537, [arXiv:1308.3705].
  15. S. R. Behbahani, A. Dymarsky, M. Mirbabayi, and L. Senatore, (Small) Resonant non-Gaussianities: Signatures of a Discrete Shift Symmetry in the Effective Field Theory of Inflation, JCAP 1212 (2012) 036, [arXiv:1111.3373].
  16. E. Silverstein and A. Westphal, Monodromy in the CMB: Gravity Waves and String Inflation, Phys.Rev. D78 (2008) 106003, [arXiv:0803.3085].
  17. N. Kaloper and L. Sorbo, A Natural Framework for Chaotic Inflation, Phys.Rev.Lett. 102 (2009) 121301, [arXiv:0811.1989].
  18. X. Dong, B. Horn, E. Silverstein, and A. Westphal, Simple exercises to flatten your potential, Phys.Rev. D84 (2011) 026011, [arXiv:1011.4521].
  19. N. Kaloper, A. Lawrence, and L. Sorbo, An Ignoble Approach to Large Field Inflation, JCAP 1103 (2011) 023, [arXiv:1101.0026].
  20. N. Kaloper and A. Lawrence, Natural chaotic inflation and ultraviolet sensitivity, Phys.Rev. D90 (2014), no. 2 023506, [arXiv:1404.2912].
  21. L. McAllister, E. Silverstein, A. Westphal, and T. Wrase, The Powers of Monodromy, JHEP 1409 (2014) 123, [arXiv:1405.3652].
  22. J. Frazer, Predictions in multifield models of inflation, Journal of Cosmology and Astroparticle Physics 2014 (2013), no. 01 028, [arXiv:1303.3611].
  23. R. Easther, J. Frazer, H. V. Peiris, and L. C. Price, Simple predictions from multifield inflationary models, Phys.Rev.Lett. 112 (2014) 161302, [arXiv:1312.4035].
  24. D. Wenren, Tilt and Tensor-to-Scalar Ratio in Multifield Monodromy Inflation, arXiv:1405.1411.
  25. L. C. Price, H. V. Peiris, J. Frazer, and R. Easther, Gravitational wave consistency relations for multifield inflation, Phys.Rev.Lett. 114 (2015), no. 3 031301, [arXiv:1409.2498].
  26. L. F. Abbott, A Mechanism for Reducing the Value of the Cosmological Constant, Phys. Lett. B150 (1985) 427.
  27. A. Achúcarro, V. Atal, and Y. Welling, On the viability of and natural inflation, JCAP 1507 (2015) 008, [arXiv:1503.07486].
  28. T. Li, Z. Li, and D. V. Nanopoulos, Helical Phase Inflation, Phys. Rev. D91 (2015), no. 6 061303, [arXiv:1409.3267].
  29. T. Li, Z. Li, and D. V. Nanopoulos, Helical Phase Inflation and Monodromy in Supergravity Theory, arXiv:1412.5093.
  30. T. Li, Z. Li, and D. V. Nanopoulos, Helical Phase Inflation via Non-Geometric Flux Compactifications: from Natural to Starobinsky-like Inflation, arXiv:1507.04687.
  31. M. Berg, E. Pajer, and S. Sjors, Dante’s Inferno, Phys. Rev. D81 (2010) 103535, [arXiv:0912.1341].
  32. J. McDonald, Sub-Planckian Two-Field Inflation Consistent with the Lyth Bound, JCAP 1409 (2014), no. 09 027, [arXiv:1404.4620].
  33. J. McDonald, A Minimal Sub-Planckian Axion Inflation Model with Large Tensor-to-Scalar Ratio, JCAP 1501 (2015), no. 01 018, [arXiv:1407.7471].
  34. G. Barenboim and W.-I. Park, Spiral Inflation, Phys. Lett. B741 (2015) 252–255, [arXiv:1412.2724].
  35. R. D. Peccei and H. R. Quinn, CP Conservation in the Presence of Instantons, Phys. Rev. Lett. 38 (1977) 1440–1443.
  36. R. Holman, S. D. Hsu, T. W. Kephart, E. W. Kolb, R. Watkins, et al., Solutions to the strong CP problem in a world with gravity, Phys.Lett. B282 (1992) 132–136, [hep-ph/9203206].
  37. S. M. Barr and D. Seckel, Planck scale corrections to axion models, Phys.Rev. D46 (1992) 539–549.
  38. S. Ghigna, M. Lusignoli, and M. Roncadelli, Instability of the invisible axion, Phys.Lett. B283 (1992) 278–281.
  39. M. Kamionkowski and J. March-Russell, Planck scale physics and the Peccei-Quinn mechanism, Phys.Lett. B282 (1992) 137–141, [hep-th/9202003].
  40. R. Kallosh, A. D. Linde, D. A. Linde, and L. Susskind, Gravity and global symmetries, Phys.Rev. D52 (1995) 912–935, [hep-th/9502069].
  41. R. Flauger, L. McAllister, E. Silverstein, and A. Westphal, Drifting Oscillations in Axion Monodromy, arXiv:1412.1814.
  42. E. D. Stewart, The Spectrum of density perturbations produced during inflation to leading order in a general slow roll approximation, Phys. Rev. D65 (2002) 103508, [astro-ph/0110322].
  43. J.-O. Gong and E. D. Stewart, The Power spectrum for a multicomponent inflaton to second order corrections in the slow roll expansion, Phys. Lett. B538 (2002) 213–222, [astro-ph/0202098].
  44. H.-C. Lee, M. Sasaki, E. D. Stewart, T. Tanaka, and S. Yokoyama, A New delta N formalism for multi-component inflation, JCAP 0510 (2005) 004, [astro-ph/0506262].
  45. M. J. Mortonson, H. V. Peiris, and R. Easther, Bayesian Analysis of Inflation: Parameter Estimation for Single Field Models, Phys.Rev. D83 (2011) 043505, [arXiv:1007.4205].
  46. R. Easther and H. V. Peiris, Bayesian Analysis of Inflation II: Model Selection and Constraints on Reheating, Phys.Rev. D85 (2012) 103533, [arXiv:1112.0326].
  47. J. Norena, C. Wagner, L. Verde, H. V. Peiris, and R. Easther, Bayesian Analysis of Inflation III: Slow Roll Reconstruction Using Model Selection, Phys.Rev. D86 (2012) 023505, [arXiv:1202.0304].
  48. L. C. Price, J. Frazer, J. Xu, H. V. Peiris, and R. Easther, MultiModeCode: An efficient numerical solver for multifield inflation, JCAP 1503 (2015), no. 03 005, [arXiv:1410.0685].
  49. R. Easther and R. Flauger, Planck Constraints on Monodromy Inflation, JCAP 1402 (2014) 037, [arXiv:1308.3736].
  50. T. Kobayashi and F. Takahashi, Running Spectral Index from Inflation with Modulations, JCAP 1101 (2011) 026, [arXiv:1011.3988].
  51. D. Baumann and L. McAllister, Inflation and String Theory, arXiv:1404.2601.
  52. N. Arkani-Hamed, H.-C. Cheng, P. Creminelli, and L. Randall, Extra natural inflation, Phys.Rev.Lett. 90 (2003) 221302, [hep-th/0301218].
  53. D. J. Schlegel, C. Bebek, H. Heetderks, S. Ho, M. Lampton, et al., BigBOSS: The Ground-Based Stage IV Dark Energy Experiment, arXiv:0904.0468.
  54. LSST Science, LSST Project Collaboration, P. A. Abell et al., LSST Science Book, Version 2.0, arXiv:0912.0201.
  55. EUCLID Collaboration, R. Laureijs et al., Euclid Definition Study Report, arXiv:1110.3193.
  56. J. Chluba, J. Hamann, and S. P. Patil, Features and New Physical Scales in Primordial Observables: Theory and Observation, arXiv:1505.01834.
  57. M. Czerny and F. Takahashi, Multi-Natural Inflation, Phys. Lett. B733 (2014) 241–246, [arXiv:1401.5212].
  58. M. Czerny, T. Higaki, and F. Takahashi, Multi-Natural Inflation in Supergravity, JHEP 05 (2014) 144, [arXiv:1403.0410].
  59. M. Czerny, T. Kobayashi, and F. Takahashi, Running Spectral Index from Large-field Inflation with Modulations Revisited, Phys. Lett. B735 (2014) 176–180, [arXiv:1403.4589].
  60. M. Czerny, T. Higaki, and F. Takahashi, Multi-Natural Inflation in Supergravity and BICEP2, Phys. Lett. B734 (2014) 167–172, [arXiv:1403.5883].
  61. T. Kobayashi, O. Seto, and Y. Yamaguchi, Axion monodromy inflation with sinusoidal corrections, PTEP 2014 (2014), no. 10 103E01, [arXiv:1404.5518].
  62. T. Higaki, T. Kobayashi, O. Seto, and Y. Yamaguchi, Axion monodromy inflation with multi-natural modulations, JCAP 1410 (2014), no. 10 025, [arXiv:1405.0775].
  63. V. Assassi, D. Baumann, D. Green, and L. McAllister, Planck-Suppressed Operators, JCAP 1401 (2014), no. 01 033, [arXiv:1304.5226].
This is a comment super asjknd jkasnjk adsnkj
The feedback cannot be empty
Comments 0
The feedback cannot be empty
Add comment

You’re adding your first comment!
How to quickly get a good reply:
  • Offer a constructive comment on the author work.
  • Add helpful links to code implementation or project page.