arxiv: v4 [hep-th] 8 Sep 2014

Transcription

1 CTPU Natural inflation with multiple sub-planckian axions Kiwoon Choi, a, Hyungjin Kim, a,b, and Seokhoon Yun a,b, a Center for Theoretical Physics of the Universe, Institute for Basic Science (IBS), Daejeon , Korea b Department of Phyiscs, KAIST, Daejeon , Korea arxiv: v4 [hep-th] 8 Sep 014 Abstract We extend the Kim-Nilles-Peloso (KNP) alignment mechanism for natural inflation to models with N > axions, which obtains super-planckian effective axion decay constant f eff M Pl through an alignment of the anomaly coefficients of multiple axions having sub-planckian fundamental decay constants f 0 M Pl. The original version of the KNP mechanism realized with two axions requires that some of the anomaly coefficients should be of the order of f eff /f 0, which would be uncomfortably large if f eff /f 0 O(100) as suggested by the recent BICEP results. We note that the KNP mechanism can be realized with the anomaly coefficients of O(1) if the number of axions N is large as N lnn ln(f eff /f 0 ), in which case the effective decay constant can be enhanced as f eff /f 0 N!n N 1 for n denoting the typical size of the integer-valued anomaly coefficients. Comparing to the other multiple axion scenario, the N-flation scenario which requires N feff /f 0, the KNP mechanism has a virtue of not invoking to a too large number of axions, although it requires a specific alignment of the anomaly coefficients, which can be achieved with a probability of O(f 0 /f eff ) under a random choice of the anomaly coefficients. We also present a simple model realizing a multiple axion monodromy along the inflaton direction. kchoi@kaist.ac.kr hjkim06@kaist.ac.kr yunsuhak@kaist.ac.kr 1

2 I. INTRODUCTION Recent detection of tensor modes in the cosmic microwave background by BICEP suggests that the inflation scale is rather high, H I GeV [1]. In the context of slow roll inflation, such a high inflation scale implies that the inflaton field has experienced a super-planckian excursion over the inflation period []. This poses a question if the inflaton field can be decoupled from high scale physics above the scale of quantum gravity, so that an effective field theory description of inflation does make sense. An attractive solution to this puzzle is to introduce an approximate continuous shift symmetry along the inflaton direction. To implement this idea, in natural inflation [3] the inflaton is assumed to be a pseudo-nambu-goldstone boson φ with a periodic potential ( )] φ V(φ) = Λ [1 cos 4, (1) f where f is the axion decay constant which defines the fundamental domain of the axionlike inflaton field: φ φ+πf. () This inflaton potential is stable against perturbative quantum corrections, which can be assured through the approximate shift symmetry φ φ + constant. Also, under a reasonable assumption on the nonperturbative dynamics generating the axion potential, one can justify that the above simple potential provides a good approximation to the full inflaton potential. However, there is a difficulty in this simple setup. For successful slow roll inflation, the axion decay constant is required to have a super-planckian value. The recent BICEP result makes the problem even more severe as it suggests [4] f 10M Pl, (3) where M Pl GeV is the reduced Planck mass. It appears to be difficult to get such a large axion decay constant from a sensible fundamental theory, particularly from string theory. In the limit where a controllable approximation is available, string theory predicts that the axion scale is typically around g M Pl /8π [5]. During the past years, there have been several proposals to circumvent this problem. They include the two-axion model of Kim et al. [6], which obtains a super-planckian effective

3 axion decay constant through an alignment of the anomaly coefficients of two axions having sub-planckian fundamental decay constants, a five-dimensional (5D) gauge-axion unification model in which the axionlike inflaton obtains a super-planckian decay constant in the limit that 5D gauge coupling becomes weaker than the 5D gravitational coupling [7], a model with non-minimal axion-gravity coupling where gravitationally enhanced Hubble friction makes natural inflation operative with sub-planckian axion decay constant [8], the N-flation scenario [9] based on the idea of assisted inflation [10] with many axions, and the axion monodromy based on either a string theoretic [11, 1] or field theoretic [13, 14] scheme for multiple windings in the axion field space. In this paper, we revisit the Kim-Nilles-Peloso (KNP) alignment mechanism to extend the scheme to models with N > axions. The original version of the KNP mechanism realized with two axions requires that some of the anomaly coefficients should be of the order of f eff /f i, where f eff M Pl is the super- Planckian effective decay constant of the axionic inflaton, while f i M Pl are the sub- Planckian fundamental axion decay constants in the model. In case that f i g M Pl /8π [5] as suggested by string theory, this would require that some anomaly coefficients should be uncomfortably large as f eff /f i = O( ). We note that the KNP mechanism can be realized with the anomaly coefficients of O(1) if the number of axions is large as N lnn ln(f eff /f i ), inwhich casetheeffective decayconstant canbeenhanced asf eff /f i N!n N 1 for n denoting the typical size of the integer-valued anomaly coefficients. We examine also the probability for the KNP alignment to be achieved under a random choice of the anomaly coefficients. In regard to enhancing the effective axion decay constant, a relevant question is how many fields do we need to get super-planckian f eff M Pl. As the Planck scale receives a quadratically divergent radiative correction from each light field, schematically we have δm Pl N l Λ, where N l denotes the number of light fields and Λ is the cutoff scale of loop momenta. Then the scheme would be in trouble if it requires a too large number of light fields as N l f eff /f i. In our multiple axion scenario, f eff/f i grows exponentially as a function of N for a fixed value of n > 1, so the number of required axions is of the order of ln(f eff /f i ). On the other hand, for the original two-axion KNP model [6], one needs n = O(f eff /f i ), where the anomaly coefficient n can be identified as the number of gaugecharged fermions generating the axion coupling to instantons. In the N-flation scenario [9], the number of required axions is O(f eff /f i ). So our scheme can enhance f eff by introducing 3

4 a parametrically smaller number of fields, as compared to the KNP two-axion model and the N-flation scenario. In certain cases, the KNP mechanism can be interpreted as enhancing the effective axion decay constant as f eff /f i n 1 through the Z n monodromy structure of a light axion, which is induced along the inflaton direction by the mixing with heavy axions. In this context, we present a simple model yielding f eff /f i N i= n i through a multiple axion monodromy described by N i= Z n i (n i > 1). We present also a model yielding f eff /f i N 1 even when all the integer-valued anomaly coefficients are restricted as n ij 1. The organization of this paper is as follows. In Sec. II, we review the original two-axion model of KNP to illustrate the basic idea and set the notations. In Sec. III, we extend the KNP mechanism to models with N > axions. Sec. IV is the conclusion. II. KIM-NILLES-PELOSO MECHANISM WITH TWO AXIONS We begin with a brief review of the original Kim-Nilles-Peloso mechanism realized with two axions [6]. In the field basis of periodic axions: φ i φ i +πf i, (4) the axion potential consistent with the axion periodicity is generically given by [ ( V(φ i ) = Λ 4 n1 φ cos + n )] [ ( φ +Λ 4 m1 φ 1 1 cos + m )] φ, (5) f 1 f f 1 f where n = (n 1,n ) and m = (m 1,m ) are linearly independent integer-valued coefficients, and f i (i = 1,) denote the fundamental axion decay constants which are presumed to be comparable to each other, while being significantly lower than the reduced Planck scale: f 1 f M Pl. Here we include only the leading nonperturbative effects generating the axion potential, under the assumption that the next order nonperturbative effects are small enough. The integer-valued coefficients n i,m i parametrize the discrete degrees of freedom in the underlying nonperturbative dynamics generating the axion potential. A simple possibility is that the axion potential is generated by hidden sector gauge field instantons through the symmetry breaking by anomalies. In such case, the model involves two non-abelian hidden 4

5 sector gauge groups G a (a = 1,), together with the gauge-charged fermions having the following couplings to axions: λ ii f i e iq iiφ i /f i ψil ψ IR +H.c., (6) I i where λ ii denote dimensionless Yukawa couplings and ψ I are assumed to be charged Dirac fermions for simplicity. Then the Noether current of the nonlinearly realized Peccei-Quinn symmetries have the U(1) i -G a -G a anomalies as U(1) i : φ i φ i +α i f i, ψil ψ IR e iq iiα i ψil ψ IR (7) µ J µ i = n i 16π F 1 F 1 + m i 16π F F, (8) where F a are the gauge field strength of the gauge group G a, and the anomaly coefficients are given by n i = I q ii Tr(T 1 (ψ I)), m i = I q ii Tr(T (ψ I)) (9) for T a (ψ I ) (a = 1,) denoting the G a -charge matrix of ψ I normalized as Tr(T a) = 1/ for the fundamental representation of G a. With this symmetry breaking by anomalies, the gauge field instantons of G a generate the axion potential of the form (5). Based on this observation, in the following we will call n i,m i the anomaly coefficients. However it should be noted that the axion potential (5) can be generated by different kinds of nonperturbative effects, for instance string theoretic instantons or hidden gaugino condensations. In such case, the integer coefficients n i,m i can be determined by a variety of different discrete UV quantum numbers, e.g. the quantized fluxes, the number of stacked D-branes, and/or the number of windings for stringy instantons. To discuss the KNP mechanism, it is convenient to consider the limit Λ = 0, in which the axion potential is given by V(φ i ) = Λ 4 1 [ ( n1 φ 1 1 cos + n )] φ. (10) f 1 f Obviously this axion potential has a one-dimensional periodic flat direction in the twodimensional fundamental domain of the axion fields spanned by φ i = [0,πf i ], φ flat n φ 1 f n 1φ f 1, (11) 5

6 φ φ 1 FIG. 1: Flat direction in the fundamental domain of axion fields in the limit Λ = 0. Even though the fundamental domain is sub-planckian with f i M Pl, the flat direction can have a super-planckian length if one (or both) of n i /gcd(n 1,n ) is large enough. The right panel depicts the flat direction in the fundamental domain for which the axion periodicity is manifest. which can be identified as the inflaton direction. One easily finds that the length of this periodic flat direction is given by l flat = π n 1 f +n f 1 gcd(n 1,n ), (1) where gcd(n 1,n ) denotes the greatest common divisor of n 1 and n. This shows that a super-planckian flat direction with l flat M Pl f i can be developed on the twodimensional sub-planckian domain if n 1 gcd(n 1,n ) or n gcd(n 1,n ) M Pl f i 1. (13) In Fig. 1, we depict the flat direction in the fundamental domain of axion fields, which has a length given by (1). Since the axionic inflaton of natural inflation rolls down along this periodic flat direction, its effective decay constant is bounded as f eff l flat π, which means that at least one of n i should be as large as gcd(n 1,n )f eff /f i. Turning on the second axion potential [ ( V = Λ 4 m1 φ 1 1 cos + m )] φ, (14) f 1 f a nontrivial potential is developed along the periodic flat direction having a length (1). Even when l flat M Pl, natural inflation is not guaranteed as the inflaton potential induced 6

7 by V generically has multiple modulations along the flat direction. We find that the number of modulations over the full range of the periodic flat direction is given by which results in the effective decay constant N mod = n 1m n m 1, (15) gcd(n 1,n ) f eff = l flat πn mod = n 1 f +n f 1 n 1 m n m 1. (16) It is straightforward to derive the above form of the effective axion decay constant [6]. Taking the rotated axion field basis: ψ = n 1f φ 1 +n f 1 φ n 1 f +n f1, φ flat = n f 1 φ 1 n 1 f φ, (17) n 1 f +n f1 the axion potential (5) can be written as V(ψ,φ flat ) = Λ 4 1 [ ( )] ψ 1 cos +Λ 4 f ψ [ ( ψ 1 cos f ψ + φ flat f eff )], (18) where f ψ = f 1 f, n 1 f +n f1 f ψ = f 1f n 1 f +n f1, n 1 m 1 f +n m f1 f eff = n 1 f +n f 1 n 1 m n m 1. (19) Taking the limit Λ 4 1 Λ4, the heavy field component ψ can be integrated out, yielding the effective potential of the light inflaton direction φ flat as [ V eff (φ flat ) = Λ 4 1 cos ( φflat f eff )]. (0) From the above expression of f eff, it is clear that f eff f i is not a generic feature of the model, but requires a special alignment of the anomaly coefficients n i,m i. Note that in the absence of any alignment, generically f eff f i n, (1) where n represents the typical size of the anomaly coefficients n i,m i. Obviously there can be many different choices (or alignments) of the anomaly coefficients yielding f eff /f i 1. A 7

8 simple possibility is that one or both of n i are large, while the denominator (n 1 m n m 1 ) is tuned to be O(1). As a specific example, KNP considered the case that n 1 and the other three anomaly coefficients are given by m 1 = n 1 = 1 and m = n + O(1) [6]. Of course, one can consider different examples as we do below, for instance n 1 with n 1 = m = 1, m 1 = 0. To understand the geometric meaning of the required alignment, one can introduce an alignment angle δθ as sinδθ n 1 m m 1 n () (n 1 +n )(m 1 +m ), and rewrite the effective decay constant as f eff = 1 ( ) n 1f +n f1 1/ sinδθ (n 1 +n )(m 1 +m ). (3) Note that δθ corresponds to the angle between the heavy field direction ψ and the other field direction φ of the second axion potential V: ψ n 1φ 1 f 1 + n φ f, φ m 1φ 1 f 1 + m φ f. (4) For the case of two axions, one of the economic ways to get f eff /f i 1 is to have n 1 n = O(f eff /f i ), m 1 m = O(1), (5) for which δθ = O(f i /f eff ). (6) This implies that the probability for achieving f eff f i under a random choice of the anomaly coefficients, but within the specific region of (5), is given by P(f eff /f i ) = O(f i /f eff ). (7) On the other hand, if one extends the random choice to the generic region where all the anomaly coefficients can be of the order of f eff /f i, one needs an alignment angle δθ = O(f i /f eff ), and then the probability for achieving f eff f i is reduced as P(f eff /f i ) = O(fi /f eff ). (8) 8

9 A particularly interesting choice [15, 16] of the anomaly coefficients realizing the KNP alignment is n 1 = m = 1, m 1 = 0, n 1, (9) for which the light inflaton direction is identified as φ flat = n f 1 φ 1 f φ f eff φ 1 (30) with the effective decay constant f eff = n f 1 +f n f 1 (31) where we assumed f 1 f. In this case, we have e iφ 1/f 1 = e in φ /f (3) along the inflaton direction. Then the enhanced effective axion decay constant can be considered as a consequence of the Z n monodromy structure of φ 1, which is induced by the mixing with the heavy axion component φ. III. GENERALIZATION TO N > AXIONS It is in fact straightforward to generalize the KNP mechanism to the case with N > axions. 1 In the presence of N axions, again in the periodic axion field basis φ i φ i +πf i (i = 1,,...,N), (33) the axion potential takes the form V = N i=1 Λ 4 i [ 1 cos ( N j=1 n ij φ j f j )], (34) where n i = (n i1,n i,..,n in ) are linearly independent integer-valued anomaly coefficients, and f i M Pl are the fundamental decay constants which are presumed to be comparable to each other. To proceed, let us take the limit Λ N = 0, 1 AgeneralizationtothecasewithN > axionshasbeendiscussedin[17]toaccommodatetheintermediate scale QCD axion as well as a quintessence axion with Planck scale decay constant in the model. 9

10 for which the one-dimensional periodic flat direction is parametrized as where φ flat i X i φ i Det φ 1 φ φ N n 11 f 1 n 1 f. n N 1,1 f 1. n N 1, f n 1N f N. n N 1,N f N, (35) for X i = C if i j f j n 11 n 1... n 1,(i 1) n 1,(i+1)... n 1N C i = ( 1) i+1 n Det 1 n... n,(i 1) n,(i+1)... n N..... n N 1,1 n N 1,... n N 1,(i 1) n N 1,(i+1)... n N 1,N (36). (37) The length of this periodic flat direction is determined by the minimal discrete shift φ i along the flat direction, under which the axion field configuration comes back to the original configuration. One then finds φ i = yielding the length of the flat direction: πc i f i gcd(c 1,C,,C N ), (38) N π i=1 l flat = C i f i gcd (C 1,C,,C N ), (39) where gcd(c 1,C,,C N ) denotes the greatest common divisor of all C i. For generic anomaly coefficients n ij, the magnitude of C i is quite sensitive to the number of axions, as well as to the typical size of n ij. Here we are interested in the limit N 1 with n ij n max = O(few). (40) To proceed, we can regard n ij as a random integer-valued variable with flat probability distribution: P(n ij ) = 10 1 n max +1.

11 We then have where n max n ij = n ij P(n ij ) = 0, n ij = n max n max n ij = n ij = n max n ij P(n ij) = 1 3 n max(1+n max ) n, n = n max (1+n max )/3 denotes the typical size of the anomaly coefficients in the range (40). One can similarly compute the expectation value of C i. For C N = σ sgn(σ)n 1σ(1) n (N 1)σ(N 1), where the summation is over all possible permutations of {1,,...,N 1}, one easily finds C N = σ (n 1σ(1) n (N 1)σ(N 1) ) = (N 1)!n (N 1), and therefore N i=1 C i = N! n (N 1). As implied by this expectation value, in most cases we have ( N 1/ Ci) N!n N 1, gcd(c 1,C,...,C N ) = O(1), (41) i=1 and therefore a flat direction enhanced as l flat πf i N!n N 1, (4) where the sub-planckian axion decays constants f i are assumed to be comparable to each other. Here we are interested in the case that the sub-planckian axions φ i originate from higher-dimensional tensor gauge fields in compactified string theory, in which case all f i are comparable to M st /8π for the string scale M st [5]. Note that the flat direction length is exponentially enhanced in the limit N 1 when the typical anomaly coefficient n > 1. As We have in fact an upper bound i C i < N N/ n N 1 max. We found through a numerical analysis that i C i 0. N!n N 1 in most cases of our interest. It is also known that a probability for gcd(c i ) = 1 under a random choice of C i within the range C i Q is given by P(gcd(C i ) = 1) = 1/ζ(N)+O(1/Q), which is close to the unity in the limit N 1 and Q 1. Although in our case C i is not a randomly chosen integer, but a determinant of the randomly chosen anomaly coefficients n ij, we confirmed again through a numerical analysis that gcd(c i ) = 1 in most cases of our interest. 11

12 a result, an exponentially long flat direction can be obtained with the anomaly coefficients n ij n max = O(few), with a moderately large number of axions N lnn ln(l flat /πf i ). (43) This can be understood by noting that the flat direction generically winds each of the additional axion dimensions by about n times, which would explain the factor n N 1, and there are also combinatoric degrees of freedom for the windings of the N-dimensional torus of axion fields, which would explain the factor N!. Introducing the last part of the axion potential [ ( V = Λ 4 N 1 cos j n Nj φ j f j a periodic potential is developed along the flat direction. Again super-planckian flat direction does not guarantee a super-planckian effective decay constant. Instead we have )], f eff = l flat πn mod, (44) where N mod is the number of modulations of the axion potential over the full period of the flat direction. Taking the limit Λ i Λ N (i = 1,,...,N 1) to integrate out the (N 1) heavy axions, we find that the effective potential of the flat direction is given by [ ( )] V eff (φ flat ) = Λ 4 φflat N 1 cos, (45) where N f eff f eff = i=1 C i f i Det N, (46) for n N1 n N n NN n DetN = Det 11 n 1 n 1N =... i n (N 1)1 n (N 1) n (N 1)N C i n Ni. (47) This tells that the number of modulations over the full range of the flat direction is N mod = i C in Ni gcd(c 1,C,,C N ). 1

13 To justify our approach to integrate out the (N 1) heavy axions, let us briefly examine the axion masses in our scheme. For the axion potential (34), the N N axion mass-square matrix is given by m kl = i Λ4 in ik n il /f k f l, yielding Detm = (DetN) N In the presence of light flat direction, its mass square is given by m flat Λ4 N f eff Then the mass-square determinant can be written as Detm = (DetN) N = m flat i C j f j f N Λ 4 i fi N 1 i = Λ4 N f eff Λ 4 i f i i Λ 4 i fi. (48) = Λ4 N (DetN). (49) j C j f j f eff f N m flat N 1 (DetN) N 1 i i Nn Λ 4 i ef i Λ 4 i fi m flat N 1 i m heavy (i), (50) where m heavy (i) is the ith heavy axion mass, and we have used j C j f j /f N N!n(N 1) N N n (N 1) /e N under the assumption that all f i are comparable to each other. We then find m heavy (i) m flat Nn e f eff f i Λ 4 i Λ 4 N, (51) which shows that the heavy axions are heavy enough compared to the flat direction if the anomaly coefficients are aligned to yield f eff f i, even when Λ N is comparable to Λ i (i = 1,...,N 1). As in the case with two axions, it is clear that f eff f i is not a generic feature, but requires a specific alignment of the anomaly coefficients. Yet, compared to the case of two axions, a notable difference is that the mechanism does not require large anomaly coefficients of O(f eff /f i ), but a moderately large number of axions N lnn ln(f eff /f i ) together with the anomaly coefficients n ij = O(1). To quantify the required degree of alignment, let us introduce an alignment angle as in the case of two axions: sinδθ C n N C n N Then the effective decay constant reads as = DetN C n N. (5) f eff = 1 sinδθ N i=1 C i f i C n N 13 f i n N sinδθ, (53)

14 implying that we need to align δθ to be small as δθ = O(f i /f eff ). (54) This also suggests that the probability for having f eff f i under a random choice of the anomaly coefficients in the range n ij n max = O(1) is given by P(f eff /f i ) = O(f i /f eff ). (55) Note that in the case of two axions we have a similar probability only when the random selection is limited to a specific (economic) region of the anomaly coefficients, for instance the region of (5), while P(f eff /f i ) = O(fi /f eff ) if one extends the random selection to the general region where all the anomaly coefficients can be comparable to each other. In view of that the fundamental axion decay constants in string theory are typically in the range f i GeV [5], while a successful natural inflation compatible with the recent BICEP results [1] requires f eff 10M Pl, we are particularly interested in the minimal number of axions which can yield f eff /f i = O( ) for n ij n max (n max = 1,,3). (56) According to our discussion above, the corresponding range of N is roughly given by N lnn ln(f eff /f i ). (57) We have performed a numerical analysis to evaluate P(f eff /f i ) for the three different values of f eff /f i : f eff /f i = 10, 5 10, 10 3, when N = 8 15 (n max = 1), N = 7 14 (n max = ), N = 5 1 (n max = 3). The results are depicted in Table I, which confirms that the probability for the necessary alignment to be achieved under a random choice of the anomaly coefficients is indeed of the order of f i /f eff. Before closing this section, let us present a couple of explicit models which achieve an exponentially enhanced effective axion decay constant within the framework discussed above. 14

15 N f eff /f i N f eff /f i N f eff /f i TABLE I: The probability to f i /f eff ratio, R = P/(f i /f eff ), for the necessary alignment under 10 5 random choices of the anomaly coefficients. We have considered N = 8 15 for n max = 1, N = 7 14 for n max =, and N = 5 1 for n max = 3. Our first model is [ V = Λ cos + Λ 4 N 1 ( φ1 [ 1 cos + n )] φ f 1 f ( φn 1 + n Nφ N f N 1 f N +Λ 4 [ 1 cos )] +Λ 4 N ( φ + n )] 3φ 3 + (58) f f ( 3 )] φn, [ 1 cos which is designed to realize a multiple axion monodromy along the inflaton direction. The anomaly coefficient matrix of the model takes the form 1 n 1 n 3 N = (n ij ) = 1..., (59)... nn 1 for which DetN = 1, C i = N j=i+1 f N n j. (60) The resulting effective axion decay constant is given by N ( N ) 1/ N f eff = = n j f i n n 3 n N f 1 (61) i=1 f i C i DetN i=1 j=i+1 15

16 φ 3 φ πf 3 /n n 3 φ 1 πf /n 3 FIG. : Multiple monodromy structure for the three-axion model with n = n 3 =. The solid red line represents the inflaton direction in the fundamental domain of three axions. Note that φ = πf along the inflaton direction requires φ 1 = πn f 1, and φ 3 = πf 3 requires φ = πn 3 f. As a result, φ 3 = πf 3 along the inflaton direction yields φ 1 = πn n 3 f 1. if n i > 1 and f 1 f f N. In the limit Λ N Λ i (i = 1,...,N 1), the inflaton direction is determined to be φ flat i C i φ i = ( N ( N n i )φ 1 n i )φ + n N φ N 1 +φ N (6) i= i=3 This model can be considered as a generalization of the two axion models of [15, 16], and realizes a multiple axion monodromy N i= Z n i along the inflaton direction. As a consequence, in order for the Nth axion φ N to travel one period along the inflaton direction, i.e. φ N = πf N, the other axions φ i (i = 1,,...,N 1) should experience a multiple winding as φ i πf i = N j i+1 n j. (63) In Fig., we depict such multiple monodromy structure for the case of N = 3 and n = n 3 =. 16

17 Our second model is [ ( V = Λ 4 φ1 1 1 cos + φ φ 3 + +( ) Nφ )] N f 1 f f 3 f [ ( N + Λ 4 φ 1 cos + φ 3 φ 4 + +( ) N 1φ )] N + f f 3 f 4 f [ ( N + Λ 4 φn N 1 cos + φ N 1 φ )] N f N f N 1 f [ ( N + Λ 4 φn 1 N 1 1 cos + φ )] [ ( )] N +Λ 4 φn N 1 cos, (64) f N 1 f N f N which is designed to achieve an exponentially enhanced effective decay constant even when all the integer-valued anomaly coefficients are restricted as n ij 1. In this case, we have DetN = 1, C i = 4 (N 1 i) (i = 1,,, N 1), C N = 1. (65) Then, assuming that all fundamental decay constants are comparable to each other, the effective decay constant is exponentially enhanced as f eff 1 3 (4N 1 +)f i N 1 3 f i (66) although we have n ij 1. IV. CONCLUSION Natural inflation provides an attractive framework for large field chaotic inflation which can explain the recent detection of primordial gravitational waves by BICEP. The KNP alignment mechanism offers an interesting scheme to get a super-planckian effective axion decay constant f eff M Pl, which is a necessary component of natural inflation, starting from sub-planckian fundamental axion decay constants f i M Pl of multiple axions. In this paper, we extended the original KNP model with two axions to models with N > axions. Compared to the original KNP model, a notable difference is that large anomaly coefficients ofo(f eff /f i )arenotrequiredanymore ifthenumber ofaxionsismoderatelylarge as N lnn ln(f eff /f i ). With such N, the effective decay constant can be exponentially enhanced asf eff /f i N!n N 1 forndenoting thetypical size oftheinteger-valued anomaly coefficients, and the probability for the necessary alignment to be achieved under a random choice of the anomaly coefficients is of the order of f i /f eff. 17

18 The structure of our setup is rich enough to realize a variety of different possibilities. For instance, it can realize a multiple axion monodromy N i= Z n i yielding f eff /f i N i= n i. The setup can also give rise to an exponentially enhanced effective axion decay constant as f eff /f i N 1, even when all the integer-valued anomaly coefficients are restricted as n ij 1. Acknowledgements We thank Chang Sub Shin, Kwang Sik Jeong and Jinn Ouk Gong for helpful discussions. [1] P. A. R. Ade et al. [BICEP Collaboration], Phys. Rev. Lett. 11, (014) [arxiv: [astro-ph.co]]. [] D. H. Lyth, Phys. Rev. Lett. 78, 1861 (1997) [hep-ph/ ]. [3] K. Freese, J. A. Frieman and A. V. Olinto, Phys. Rev. Lett. 65, 333 (1990); F. C. Adams, J. R. Bond, K. Freese, J. A. Frieman and A. V. Olinto, Phys. Rev. D 47, 46 (1993) [hep-ph/90745]. [4] K. Freese and W. H. Kinney, arxiv: [astro-ph.co]. [5] K. Choi and J. E. Kim, Phys. Lett. B 154, 393 (1985) [Erratum-ibid. 156B, 45 (1985)]; P. Svrcek and E. Witten, JHEP 0606, 051 (006) [hep-th/060506]; T. Banks, M. Dine, P. J. Fox and E. Gorbatov, JCAP 0306, 001 (003) [hep-th/03035]. [6] J. E. Kim, H. P. Nilles and M. Peloso, JCAP 0501, 005 (005) [hep-ph/ ]. [7] N. Arkani-Hamed, H. -C. Cheng, P. Creminelli and L. Randall, Phys. Rev. Lett. 90, 130 (003) [hep-th/030118]. [8] C. Germani and A. Kehagias, Phys. Rev. Lett. 106, (011) [arxiv: [hep-ph]]; C. Germani and Y. Watanabe, JCAP 1107, 031 (011) [Addendum-ibid. 1107, A01 (011)] [arxiv: [astro-ph.co]]. [9] S. Dimopoulos, S. Kachru, J. McGreevy and J. G. Wacker, JCAP 0808, 003 (008) [hep-th/050705]. [10] A. R. Liddle, A. Mazumdar and F. E. Schunck, Phys. Rev. D 58, (1998) [astro-ph/ ]; E. J. Copeland, A. Mazumdar and N. J. Nunes, Phys. Rev. D 60, 18

19 (1999) [astro-ph/ ]; P. Kanti and K. A. Olive, Phys. Lett. B 464, 19 (1999) [hep-ph/ ]. [11] E. Silverstein and A. Westphal, Phys. Rev. D 78, (008) [arxiv: [hepth]]; L. McAllister, E. Silverstein and A. Westphal, Phys. Rev. D 8, (010) [arxiv: [hep-th]]. [1] F. Marchesano, G. Shiu and A. M. Uranga, arxiv: [hep-th]; A. Hebecker, S. C. Kraus and L. T. Witkowski, arxiv: [hep-th]. [13] N. Kaloper and L. Sorbo, Phys. Rev. Lett. 10, (009) [arxiv: [hep-th]]; N. Kaloper, A. Lawrence and L. Sorbo, JCAP 1103, 03 (011) [arxiv: [hep-th]]. [14] K. Harigaya and M. Ibe, arxiv: [hep-ph]. [15] S. -H. H. Tye and S. S. C. Wong, arxiv: [astro-ph.co] [16] I. Ben-Dayan, F. G. Pedro and A. Westphal, arxiv: [hep-th]. [17] A. Chatzistavrakidis, E. Erfani, H. P. Nilles and I. Zavala, JCAP 109, 006 (01) [arxiv: [hep-ph]]. 19