We discuss the role of a particular combination of higher derivative terms in higher dimensional theories, particularly in the background of spontaneous compactification. Two classes of theories are proposed in this paper. The first model as a generalization of the critical gravity with the Maxwell field could have a de Sitter solution. We consider the Lanczos-Lovelock term and Horndeski term as well as the higher-order Maxwell term for the second model, which contains a possible longer expansion time for the inflationary phase. It is interesting that both models can be regarded as the generalization of the Randjbar-Daemi, Salam and Strathdee (RSS) model and give the well behavior for inflation stage under the specific assumptions.
1. Introduction
Over several years, the gravitational theory with higher derivative terms [1, 2] has been a highly interesting subject studied in physics. The treatment of Einstein gravity with perturbative field quantization naturally leads to the inclusion of such nonlinear terms in the spacetime curvatures for effective actions.
The generalization of Einstein gravity has been developed extensively by many authors. The Lanczos-Lovelock gravity [3–6] (for a review, see [7]), which is dictated by effective action with the dimensionally continued Euler forms, is defined in a generic dimensional spacetime, and it describes a massless spin-two fluctuation despite the existence of higher-order terms in its action. About a few decades ago, it was discovered that the effective field theory arising from string theory involves similar terms at least for the lowest-order correction to the Einstein-Hilbert term [8, 9].
The other specific extension is found in the study on the so-called critical gravity [10–12] (although the word “critical” originally implied the emergence of Log gravity [13–15]). In this gravitational theory, the behavior of spin-two modes is described by the Lagrangian of a special combination of higher-order terms in the Ricci tensor and the scalar curvature. There is a ghost field but no spin-zero mode in this theory. In the three-dimensional spacetime, the model governed by a similar Lagrangian expresses the new massive gravity [16, 17].
In this study, we consider incorporation of the Maxwell fields to the (next-to) critical gravity, which has been studied under the assumption of the maximally symmetric background spacetime. The higher derivative theory for the other type of background fields is worth studying. Inclusion of a nonvanishing flux field brings about partial compactification with a maximally symmetric extra space and might provide an interesting scenario for a very early stage of the universe.
The assumption of the partially symmetric solution could yield another possibility. It is feasible to attach the other type of higher-order terms to the Lagrangian which realizes the spontaneous compactification. Due to the symmetry, the use of the Riemann tensor is possible in the additional terms, which are expressed by special combinations of Gauss-Bonnet terms and Horndeski’s generalized Maxwell term [18]. These terms induce a second-order field equation displaying the cosmological evolution of the background geometry (the applications of the Horndeski’s vector-tensor theory have been discussed previously in four-dimensional [19, 20] and in six-dimensional cosmology [21]).
This paper is organized as follows. In Section 2, we build the Lagrangian of our first model, which involves the quadratic terms of the Ricci tensor, the Ricci scalar, and the two-form flux field, in a similar way as revealed before in the study of critical gravity. The spontaneous compactification in the first model is studied in Section 3 and cosmological solutions are examined in Section 4. In Section 5, we present the second model containing the quadratic term of the Riemann tensor. We also investigate the cosmological solution qualitatively in this model. The last section is devoted to the summary and further research prospects in this area of study.
It is to be noted that throughout this study, we deal with our models at the classical level only.
2. A Critical Modification of Einstein-Maxwell Theory
In the study on critical gravity, it is found that the Lagrangian of the specific combination of curvatures yields a spin-two massless graviton and a ghost mode; there appears no scalar mode as in Einstein gravity. In this section, we propose a model of higher-order extension of the Maxwell-Einstein theory. We will follow the method for constructing the Lagrangian examined in our previous study [22, 23]. In [22, 23], we constructed the action of higher derivative pure gravity which admits the (anti-)de Sitter solution. In similar sense, we manage to construct the model with higher-order terms in which background fields satisfy the Einstein equation in the presence of the Maxwell field in this paper.
First, we start with the action for the Einstein gravity and electromagnetism in D dimensions:(1)I0=∫dDx-gL0=∫dDx-gR-FMNFMN-Λ, where R and Λ denote the scalar curvature obtained from the metric gMN and the cosmological constant, respectively. The field strength FMN is defined by FMN≡∂MAN-∂NAM with a vector field AM and F2≡FMNFMN. The suffixes M,N,… vary over 0,1,…,D-1.
Apply the variational principle to the action results in the following equations of motion:(2)TMN=0 with(3)TMN≡1-gδI0δgMN=RMN-12RgMN-2FMN2-14F2gMN+12ΛgMN, where FMN2≡FMPFNP and(4)JN≡∇MFMN=0.
Next, we construct a Lagrangian for our model as follows:(5)L=αL0-βTMNΔMNPQTPQ, where(6)ΔMNPQ≡12gMPgNQ+gMQgNP-1D-1gMNgPQ. Note that the solution of TMN=0 is the solution for the equation of motion derived from L, provided that the field equation for the vector field is satisfied (though the stability is not necessarily guaranteed). (Generally speaking, we can add the term ∝JMJM to our Lagrangian, but we restrict ourselves to the case without such a term in the present paper.)
The Lagrangian L can then be written as(7)L=αR-F2-Λ-βRMNRMN-D4D-1R2-4FMN2RMN+D+22D-1F2R+4FMN2F2MN-D+84D-1F22+Λ2D-1D-2R-D-4F2-D2Λ.
Incidentally, we can pick up some “critical” cases: if α=D/4(D-1)βΛ, the contribution of the cosmological constant is zero in the D-dimensional spacetime; if α=D-2/2(D-1)βΛ, the Einstein-Hilbert term (R) is absent from the D-dimensional action; and if α=D-4/2(D-1)βΛ, the pure Maxwell term (F2) is absent from the D-dimensional action.
We should note that the auxiliary field method [14, 16, 17, 23], which can be used in critical higher-order gravity, also leads to the additional terms such as(8)L′=αL0-β2TMNSMN-SMNSMN+S2, where SMN is an auxiliary symmetric tensor field and S≡SMNgMN. Varying the auxiliary field SMN, we obtain(9)TMN=SMN-SgMN,or, solving this equation with respect to SMN, we get(10)SMN=TMN-1D-1TgMN=ΔMNPQTPQ. Therefore, one can see the equivalence between the Lagrangians, L and L′. We can derive the linearized field equation using the action I′=∫dDxL′.
Now, we investigate the behavior of the linear fluctuations around the background geometry, which satisfies TMN=0. We decompose the metric as follows:(11)gMN=g¯MN+hMN. The indices are raised and lowered by the background metric, g¯MN. Then, the trace of the fluctuation is expressed as(12)h≡g¯MNhMN. For convenience, we write down the curvature tensors up to the linear order in hMN here as [24](13)RMNPQ=R¯MNPQ-12R¯MLPQhLN-12R¯NLPQhLM+12-∇¯P∇¯MhQN+∇¯P∇¯NhQM+∇¯Q∇¯MhPN-∇¯Q∇¯NhPM,RNM=R¯NM-R¯MPNQhPQ+12R¯LMhNL-12R¯NLhLM+12-∇¯2hNM-∇¯N∇¯Mh+∇¯M∇¯LhNL+∇¯N∇¯LhLM,R=R¯-R¯MNhMN+∇¯M∇¯NhMN-∇¯2h,where the nabla stands for the covariant derivative and the barred symbols show that they are constructed from the background metric, g¯MN.
Assuming the background field satisfies T¯MN=R¯MN-1/2R¯g¯MN-2F¯MN2-1/4F¯2g¯MN+1/2Λg¯MN=0, we can use(14)δ2SPQTPQ-SPQSPQ+S2=2SPQδTPQ+higher-order terms in small fluctuations from the background, to obtain the field equation with a linearized order:(15)1-gδI′δgMN=0⟶-αTMN-β-∇¯2SMN-∇¯M∇¯NS+∇¯M∇¯LSLN+∇¯N∇¯LSLM-∇¯P∇¯QSPQg¯MN+∇¯2Sg¯MN-2R¯PMQNSPQ+R¯PQSPQg¯MN-R¯SMN+R¯LMSLN+R¯LNSLM+ΛSMN+4F¯MPF¯NQSPQ+F¯2SMN-2F¯PQ2SPQg¯MN=0.Noting that SMN is in the linear order (S¯MN=0), we have used the condition T¯MN=0 for the background field.
The generalized Maxwell equation in the linear order of fluctuations takes the form(16)1-gδI′δAN=0⟶2α∇¯MfMN-2β∇¯M2F¯MPSPN-2F¯NPSPM-F¯MNS=0, where fMN=FMN-F¯MN.
The linearized equation obtained from (9) is(17)TMN=SMN-gMNS⟶12-∇¯2hMN-∇¯M∇¯Nh+∇¯M∇¯LhLN+∇¯N∇¯LhLM-∇¯P∇¯QhPQg¯MN+∇¯2hg¯MN-2R¯PMQNhPQ+R¯PQhPQg¯MN-R¯hMN+R¯LMhLN+R¯LNhLM+ΛhMN+2F¯MPF¯NQhPQ+12F¯2hMN-F¯PQ2hPQg¯MN-2F¯PMfPN-2F¯PNfPM+F¯PQfPQg¯MN=SMN-g¯MNS.
The important observation is that, using (9), the field equation (15) becomes(18)-∇¯2SMN-∇¯M∇¯NS+∇¯M∇¯LSLN+∇¯N∇¯LSLM-∇¯P∇¯QSPQg¯MN+∇¯2Sg¯MN+αβSMN-g¯MNS-2R¯PMQNSPQ+R¯PQSPQg¯MN-R¯SMN+R¯LMSLN+R¯LNSLM+ΛSMN+4F¯MPF¯NQSPQ+F¯2SMN-2F¯PQ2SPQg¯MN=0.In the D-dimensional Minkowski vacuum, this equation is just the Fierz-Pauli equation [25] for the spin-two wave with mass-squared m2=α/β. Therefore, we recognize that SMN corresponds to the massive ghost field. The postulation of no tachyonic ghost requires α/β>0 in the present model.
The linear fluctuation from the background fields can be investigated by using the above equations. In the next section, we consider the solution for spontaneous compactification of spacetime as the background.
3. Spontaneous Compactification in Critically Modified Higher-Order Einstein-Maxwell Theory
In this section, we study the partial compactification of space with a nontrivial flux in our model. The extra dimensions are considered to be compactified into a sufficiently small size. A partially symmetric solution in our model describes compactification of the form MD-2⊗S2, where MD-2 is the (D-2) dimensional Minkowski spacetime and S2 is a two-dimensional sphere, as in the work of Randjbar-Daemi, Salam and Strathdee (RSS) [26], where D=6 is considered. We consider the solution for the flat D-2 dimensional spacetime and study its classical stability in this section, though other symmetric spacetimes are also interesting. Our notation is that the suffixes μ,ν,… run over 0,1,2,D-3 and the suffixes m,n,… are used for the extra dimensions of D-2 and D-1.
As stated in the previous section, the solution for TMN=0 and JM=0 is a solution in our model. A suitable ansatz for the flux field is(19)F=12FMNdxM∧dxN=12qɛmndxm∧dxn, where ɛmn is the antisymmetric symbol and the constant q is the strength of the two-form flux. (It can be checked that this configuration is the solution of our model.) This monopole flux yields the relation,(20)Fmn2=12F2gmn.
Then, the equation, TMN=0, can be decomposed to(21)Rμν=-1D-2F2gμν+1D-2Λgμν,Rmn=2Fmn2-1D-2F2gmn+1D-2Λgmn=D-3D-2F2gmn+1D-2Λgmn,for the spacetime of a direct product. As for the flat D-2 dimensional spacetime, where Rμν=0, we consider a fine-tuning of the parameters as(22)F¯2=Λ>0.And, therefore, the solution for the background fields reads(23)R¯=2F¯2,R¯mn=Λg¯mn=2F¯mn2,R¯mnpq=Λδpmδqn-δqmδpn=2F¯mnF¯pq.
Substituting this compactified background solution into (18), we obtain(24)-∇¯2SMN-∇¯M∇¯NS+∇¯M∇¯lSLN+∇¯N∇¯lSLM-∇¯P∇¯QSPQg¯MN+∇¯2Sg¯MN+R¯LMSLN+R¯LNSLM+αβSMN-g¯MNS=0.
To show the complete particle spectrum for the ghost field is rather complicated, but we can see that the lowest mass of the ghost is α/β for symmetric tensor fields (which does not couple to the Ricci curvature of the extra space) and vector bosons of the Kaluza-Klein origin (which come from the zero modes of the Lichnerowicz operator on S2, -∇¯2+R¯mn). It is also remarkable that the other graviton and vector field fluctuations are absent in the linearized equation (24).
Thus, taking the analysis by RSS for other fluctuation modes [26] into account, in the classical and linearized analysis, we find that the RSS background solution MD-2⊗S2 is stable if α/β>0, because there is no growing mode. (Of course, there remains the problem of nonlinear instability and the ghost problem in quantum physics.)
Before closing this section, we note that the effective D-2 dimensional Newton constant, which can be found by the coefficient of the Ricci scalar of the D-2 dimensional spacetime in the action, is not affected by the higher-order terms in our model and is independent of β, if the conditions on the background fields (22) and (23) are satisfied. (If we express the D-dimensional Ricci scalar as R=R¯+RD-2, where RD-2 is the Ricci scalar of the D-2 dimensional spacetime, the coefficient of RD-2 in the Lagrangian L (7) is found to be α-1/2(D-1)β-DR¯+(D+2)F¯2-(D-2)Λ.) It might be related to the absence of the graviton and the ghost modes in (24) at the linearized level mentioned above.
4. De Sitter Solutions and Cosmology in the Model
In this section, we perform the analysis on the stability of the compactification and the cosmological solution in our model. We consider an effective potential V for the radius of extra space S2 [27] in order to study them.
We can express the Riemann tensor of the S2 as(25)Rmnpq=1b2δpmδqn-δqmδpn, where b denotes the radius of the compactified sphere S2. Then, (21) with Rμν=0 impose [26](26)q2=12Λ,b2=1Λ, where the first equation represents the fine-tuning condition between the parameters.
The effective potential V as a function of b, the scale of the extra space, is obtained by replacing the background fields with functions of b in -I such that(27)Rnm=1b2δnm,F2nm=q2b4δnm=12Λb4δnm, and detg¯mn∝b2. The effective potential for the RSS model is calculated from -I0 and is given as(28)V0Λb2=-b22b2-2q2b4-Λ=-2+1Λb2+Λb2. Therefore, the effective potential for our model of which the Lagrangian is described by (7) is(29)Vy=α-2+1y+y+βΛ4D-14D-2-Dy+2Dy-12D-2y2+7D-16y3=y-124D-1y37D-16β+2D-4βy+4D-1α-Dβy2, with y≡Λb2. A minimum of the potential is found at y=1 and attains V(1)=0. The stability of the radius of the extra sphere is determined by the sign of the second derivative of the potential:(30)V′′1=2α+2D-3D-1βΛ. If D>3 and α,β>0, the scale b=1/Λ is found to be stable. We also find that it is impossible to create another potential minimum satisfying V=0 by any fine-tuning of β/α in the model.
At a large value of b, V can be indefinitely negative provided that D/4(D-1)βΛ-α>0. In this case, the D-dimensional cosmological constant becomes negative. This may not indicate instability of compactification against a large fluctuation, because the higher derivative terms induce a coupling between the compactification scale and the first and second derivatives of an expanding scale factor for the rest of the space.
Therefore, we investigate the cosmological equation including the dynamical scale factor and the compactification scale. The metric is usually assumed as(31)ds2=-dt2+a2tdΩD-32+b2tdΩ22, where b(t) denotes the radius of the compactified sphere S2 and dΩ22 is the line element in the unit sphere. For the maximally symmetric (D-3) dimensional space, the line element is denoted by dΩD-32 here, and its Riemann tensor is normalized to R~ijkl=k(δkiδlj-δliδkj) with k=0,±1, where i,j,…=1,2,…,D-3. The symbol a(t) is the scale factor for the D-3 dimensional homogeneous and isotropic space.
The Riemann curvature is computed with the metric as(32)R0i0j=a¨aδji,R0m0n=b¨bδnm,Rimjn=a˙ab˙bδjiδnm,Rijkl=a˙2+ka2δkiδlj-δliδkj,Rmnpq=b˙2+1b2δpmδqn-δqmδpn, where the dot (˙) denotes the derivative with respect to the cosmic time t.
Now we can obtain the effective action of the scale factor, a(t), and the compactification scale, b(t), for the RSS model as(33)I0⟶∫dtaD-3b22D-3a¨a+4b¨b+4D-3a˙ab˙b+D-3D-4a˙2+ka2+2b˙2+1b2-1Λb4-Λ~∫dtaD-3b2-4D-3a˙ab˙b+D-3D-4-a˙2+ka2+2-b˙2+1b2-1Λb4-Λ≡I~0, where the overall constant which comes from the volume of the space has been omitted. The action written after the symbol “~” has been integrated by parts.
The effective action for our model can be written, by adopting the auxiliary fields, as(34)I′⟶I~0-β∫dtaD-3b22s0-D-3D-42a˙2+ka2-2D-3a˙ab˙b-b˙2+1b2+12Λb4+12Λ+2D-3sa-D-4a¨a-2b¨b-2D-4a˙ab˙b-D-4D-52a˙2+ka2-b˙2+1b2+12Λb4+12Λ+4sb-D-3a¨a-b¨b-D-3a˙ab˙b-D-3D-42a˙2+ka2-12Λb4+12Λ-s02+D-3sa2+2sb2+s0+D-3sa+2sb2, where we set S00=s0, Sji=saδji, and Snm=sbδnm. Here, the equations of motion are omitted, but they can be obtained from the above action. (The Euler equations can be obtained from the action using the command EulerEquations in the add-on Calculus “VariationalMethods” for Mathematica [28].)
We seek for the de Sitter (⊗S2) solution, where k=0, a(t)=eH0t with constant H0, and b=b0= constant are taken as an ansatz. In this case, all the equations of motion become algebraic equations. With the ansatz, the auxiliary variables are expressed as(35)s0=sa=3-2Λb02+DD-3H02-ΛΛb042D-1Λb04,sb=5-2D+2D-2Λb02-D-2D-3H02+ΛΛb042D-1Λb04. Substituting these auxiliary relations to the other equations, we obtain the following two equations at last:(36)-7D-16β+12D-2βΛb02-22D-1α+D+2D-3D-4H02+DΛβΛb04+42D-1α+DD-3D-4H02-D-2ΛβΛ2b06+D-32D-6D2-6D+4H04β+2D-3D-4H022D-1α-D-2βΛ+-4D-1α+DβΛΛΛ2b08=0,(37)37D-16β-24D-2βΛb02+22D-1α+D+2D-2D-3H02+DΛβΛb04+D-32D-2D2-6D+4H04β+2D-2D-3H022D-1α-D-2βΛ+-4D-1α+DβΛΛΛ2b08=0.
The two sets of the solution for b0 and H0 can be found; one is(38)(i) H02=0,b02=1Λ,and the other is(39)(ii) H02=D-4D-33Λ,b02=D-3Λ.These values for solutions are the same as those of the RSS model [29], as expected (if D=6). In both cases, the values for s0, sa, and sb are equal to zero.
Since the simultaneous equations ((36), (37)) are fourth-order algebraic equations, there are two more sets of solutions (in general, the values for s0, sa, and sb are not zero for such solutions) at the most. We may infer that for sufficiently small value of β, only two original solutions ((38), (39)) are real solutions. For D=6, only two real solutions ((38), (39)) exist for βΛ/α<≈2.8. As a model of the universe, the multiple de Sitter phases would have some significance in cosmology.
We now study the behavior of the small fluctuations around the above exact solutions. In the neighborhood of (i), we set(40)bt=1Λ+Δbt,Ht=a˙a=ΔHt,Δsa=sa, and so on. Then, the linearized equations of motion can be expressed as(41)D-3Δsa+2Δsb=0,(42)D-1Δs¨a+2D-3Λ+D-1αβΔsa=0,(43)2D-1αD-2Δb¨+2D-4ΛΔb+D-3ΛD-1D-2α+4βΛΔsa=0,(44)D-1αD-2ΔH˙-4Λ3/2Δb-D-1D-2α-2D-3βΛΔsa=0,(45)D-1Δs0=-2D-3βΛαΔsa. Note that, for β=0, the above equations give Δsa=Δsb=Δs0=0, (D-2)Δb¨+2(D-4)ΛΔb=0, and (D-2)ΔH˙-4Λ3/2Δb=0, which are equivalent to the equations for the cosmology of the RSS model [29].
The oscillatory behavior of Δb as well as Δsa is already interesting enough. The variables Δs0, Δsa, and Δsb oscillate with a different frequency from that of Δb, as seen from (42) and (43). (Note that the angular frequency of the oscillation of Δs0,a,b is the square root of V′′(1)/2β in (30).) Thus, the tuning of the initial conditions provides enough time for accelerating expansion and a sufficient inflation can be achieved, though we cannot predict the choice of the initial state in the present classical analysis of the model.
The analysis of the perturbation from (ii) can be performed and it results in finding exponentially growing modes, as found in [29], though we omit the indication of equations here. A certain fine-tuning could give enough expansion time for a sufficient inflation as the de Sitter phase also in this case.
5. Spontaneous Compactification in Lovelock-Horndeski-Nonlinear Maxwell Theory
In the solution for the spontaneous compactification in previous sections, one can deduce the following relation in the background fields:(46)R¯MNPQ=2F¯MNF¯PQ.Similar to the consideration in the previous sections, we come up with an idea that the Minkowski compactification of the RSS model is attained in a higher derivative theory if the relation (46) is reflected in the higher-order correction term.
We then propose an additional higher-order term(47)L2=14δMNKLPQRSRPQMN-2FMNFPQRRSKL-2FKLFRS=RMNPQRMNPQ-4RMNRMN+R2-4RMNPQFMNFPQ-4RMNFMN2+RF2+4F2MNFMN2-8F22, where generalized Kronecker’s delta is defined by(48)δN1N2⋯NpM1M2⋯Mp≡δN1M1δN2M1⋯δNpM1δN1M2δN2M2⋯δNpM2⋮⋮⋱⋮δN1MpδN2Mp⋯δNpMp. Now, we write model 2 by the Lagrangian (therefore, we call the previous model model 1):(49)L=αL0+γL2. Note that L2 consists of the Lanczos-Lovelock Lagrangian (Gauss-Bonnet term) [3–6], the Horndeski’s vector-tensor term [18], and a nonlinear Maxwell term.
It is easy to see that adding L2 to the RSS Lagrangian L0 yields no modification of the linearized equation of motion on the RSS background fields after gauge fixing. Therefore, the Kaluza-Klein modes remain unchanged from those in the RSS model. It is worth noting that the effective D-2 dimensional Newton constant in model 2 is not affected by the value of γ.
Now, we study the cosmological evolution of the classical background in the model. Although the Kaluza-Klein cosmologies with the Gauss-Bonnet term or the generalized Euler invariant have been investigated in many studies including [30–35], the case with the RSS-type compactification with the Gauss-Bonnet term has been studied only in a few papers including [36]. (The compactification with the Horndeski Lagrangian has been considered by one of the present authors [21].)
As previously considered, we characterize the cosmological background by(50)ds2=-dt2+a2tdΩD-32+b2tdΩ22,F=12qɛmndxm∧dxn=1212Λɛmndxm∧dxn. Then, the Lagrangian L2 reduces to(51)L2=4D-3a¨a4D-4a˙ab˙b+D-4D-5a˙2+ka2+2b˙2+1b2-21Λb4+8b¨b2D-3a˙ab˙b+D-3D-4a˙2+ka2+D-3D-4D-5D-6a˙2+ka22+8D-3D-4D-5a˙ab˙ba˙2+ka2+8D-3D-4a˙2a2b˙2b2+4D-3D-4a˙2+ka2b˙2+1b2-1Λb4, and we obtain the following action for a and b:(52)I~2=∫dtaD-3b2D-3D-4D-5D-6-13a˙4a4-2ka2a˙2a2+k2a4+8D-3D-4D-5-13a˙3a3b˙b-ka2a˙ab˙b+4D-3D-4-a˙2a2b˙2+1b2+ka2-b˙2+1b2-4D-3D-4-a˙2+ka21Λb4-16D-3a˙ab˙b1Λb4. Again, we omit the exposition of equations of motion here.
We first explore the de Sitter phase or the maximally symmetric product spacetime. To find this, we set k=0, a(t)=eH0t, and b(t)=b0, where H0 and b0 are constants. We then obtain the simultaneous algebraic equations as follows:(53)-α+4D-3D-4H02γ+2α+2D-3D-4H02γΛb02+D-3D-4H02-Λα+D-3D-4D-5D-6H04γΛb04=0,(54)α+4D-2D-3H02γ+D-2D-3H02-Λα+D-2D-3D-4D-5H04γΛb04=0. Of course, the solution (i) (38) is a trivial solution of the above equations, while (ii) (39) is not a solution for γ≠0. One solution of the simultaneous equations ((53), (54)) with H0≠0 exists for all positive γ/α. The value of b0(>1/Λ) increases with the value of γ/α. The existence of the de Sitter phase with large values for H0 and b0 is suitable for the possible inflationary phase.
Next, we draw the phase diagram of the compactification scale b. We consider the case with k=0 and D=6. Figure 1 shows the phase space spanned by b and b˙ for γΛ/α=0,1,2. The small arrows in the figures indicate the flow in the phase space. One can see that b=1/Λ is an attractor and the area of the attracting region increases with the value of γ/α. (If ordinary matter is added, the dissipative motion becomes more apparent in an expanding era [37].) We conclude that an appropriate choice of the initial condition, which governs the flow passing through the branch point given by the de Sitter solution of (53) and (54), might result in a longer expansion time.
Phase space for model 2. We consider here a unit such that Λ=1. We show (a) for γ/α=0, (b) for γ/α=1, and (c) for γ/α=2, respectively.
6. Summary and Prospects
In this paper, we have considered special cases for extension of the Einstein-Maxwell theory with higher-order terms. The models proposed in the present paper have the same solution for spontaneous compactification MD-2⊗S2 utilizing the internal flux as the RSS model. The stability of the compactified vacuum has been confirmed at the classical level. Nevertheless, the models reveal the different cosmological behaviors because of the additional degrees of freedom or the existence of the additional de Sitter phases. We should investigate these evolutional equations by using numerical calculations in future work.
Our toy models seem to be quite special, but the analysis of the models will exhibit some typical behaviors including inflation, and we can regard general higher-order corrections as a modification to the critical models.
To obtain a realistic inflationary epoch which successfully explain the Planck result on [38, 39], we should analyze the models further in detail. Although the possible multiple de Sitter phases and the oscillations around the equilibrium scale certainly lead to an inflationary expansion of the universe in model 1, it is difficult to specify the initial condition because the dynamical compactification has uncertainty in the beginning of the multidimensional universe. Moreover, initial quantum fluctuations in higher dimensions should be reconsidered to obtain precise values of cosmological perturbations. Due to the same reasons, at the present stage, it is ambiguous whether the possible (tunable) slow-roll inflation in model 2 can explain Planck results. Extensive calculations for various cases are needed to conclude whether our model can completely explain the cosmological observables. (To try to explain the dark energy problem, we will face the fine-tuning problem in our models, as well as the RSS model.)
Recently, the black holes in higher derivative gravity have attracted renewed interest [40–42]. We consider that our model 1 is appropriate to examine magnetized black holes of a novel type.
The other subjects we should consider are, for instance, the coupling to matter fields; compactifications on (S2)N [43] in our models; the inclusion of a p-form flux field such as the Freund-Rubin compactification [44]; the treatment of our models in quantum cosmology [45–47]; the initial singularity problems [48] in the models; the possible supersymmetrization of the models; and the field content of the model in 3+2 dimensions.
The generalization of our model 1 to higher-order theory is considered a straightforward exercise. If we introduce an auxiliary field, SMN, the additional term can be replaced by(55)nδN1⋯Nn-1NM1⋯Mn-1MSM1N1⋯SMn-1Nn-1ΔMNPQTPQ-SMN+δN1⋯NnM1⋯MnSM1N1⋯SMnNn, which will be studied extensively elsewhere.
Finally, we must research the relation of our work to other theories. We suppose that the structure of model 1 owes itself to the matrix ΔMNPQ, which appears also in the spin-two propagator. It would be very interesting if our models could be regarded as an effective reformulation or revision of another theory.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
StelleK. S.Renormalization of higher-derivative quantum gravity197716495396910.1103/physrevd.16.953MR0452429StelleK. S.Classical gravity with higher derivatives197894353371LanczosC.Elektromagnetismus als natürliche Eigenschaft der Riemannschen Geometrie1932733-414716810.1007/BF013512102-s2.0-33748302419LanczosC.A remarkable property of the Riemann-Christoffel tensor in four dimensions193839484285010.2307/1968467LovelockD.The Einstein tensor and its generalizations1971123, article 49810.1063/1.1665613MR02758352-s2.0-0001748331LovelockD.The four-dimensionality of space and the Einstein tensor19721387487610.1063/1.1666069MR0297306ZBL0234.530202-s2.0-1642301328PadmanabhanT.KothawalaD.Lanczos-Lovelock models of gravity2013531311517110.1016/j.physrep.2013.05.007MR31440232-s2.0-84884595384ZuminoB.Gravity theories in more than four dimensions1986137110911410.1016/0370-1573(86)90076-1MR839214ZwiebachB.Curvature squared terms and string theories19851565-631531710.1016/0370-2693(85)91616-8LüH.PopeC. N.Critical gravity in four dimensions20111061818130210.1103/physrevlett.106.181302DeserS.LiuH.LüH.PopeC. N.ŞişmanT. Ç.TekinB.Critical points of D-dimensional extended gravities2011836061502(R)10.1103/physrevd.83.061502LüH.PangY.PopeC. N.Conformal gravity and extensions of critical gravity2011846906400110.1103/PhysRevD.84.064001AlishahihaM.FareghbalR.D-dimensional log gravity201183808405210.1103/physrevd.83.084052BergshoeffE. A.HohmO.RosseelJ.TownsendP. K.Modes of log gravity2011831010403810.1103/physrevd.83.104038JohanssonN.NasehA.ZojerT.Holographic two-point functions for 4d log-gravity201220129, article 11410.1007/JHEP09(2012)114BergshoeffE. A.HohmO.TownsendP. K.Massive gravity in three dimensions20091022020130110.1103/physrevlett.102.201301MR2507853BergshoeffE. A.HohmO.TownsendP. K.More on massive 3D gravity2009791212404210.1103/physrevd.79.124042HorndeskiG. W.Conservation of charge and the Einstein–Maxwell field equations197617111980198710.1063/1.522837MR0429007BarrowJ. D.ThorsrudM.YamamotoK.Cosmologies in Horndeski's second-order vector-tensor theory201320132, article 14610.1007/JHEP02(2013)146Beltrán JiménezJ.DurrerR.HeisenbergL.ThorsrudM.Stability of Horndeski vector-tensor interactions2013201310, article 064YoshidaK.ShiraishiK.A new vector-tensor theory and higher-dimensional cosmology199143212913210.1088/0031-8949/43/2/001MR1092694KanN.KobayashiK.ShiraishiK.‘Critical’ cosmology in higher order gravity20132013965168410.1155/2013/651684KanN.KobayashiK.ShiraishiK.Critical higher order gravities in higher dimensions201388404403510.1103/physrevd.88.044035BarthN. H.ChristensenS. M.Quantizing fourth-order gravity theories: the functional integral19832881876189310.1103/physrevd.28.1876MR720986FierzM.PauliW.On relativistic wave equations for particles of arbitrary spin in an electromagnetic field193917395321123210.1098/rspa.1939.0140Randjbar-DaemiS.SalamA.StrathdeeJ.Spontaneous compactification in six-dimensional Einstein-Maxwell theory1983214349151210.1016/0550-3213(83)90247-xMR6988942-s2.0-0001393994WetterichC.Spontaneous compactification in higher dimensional gravity1982113537738110.1016/0370-2693(82)90768-7Wolframhttp://www.wolfram.com/OkadaY.Inflation in Kaluza-Klein cosmology19851501–310310610.1016/0370-2693(85)90148-0MR775112MadoreJ.Kaluza-Klein theory with the Lanczos lagrangian1985110628929210.1016/0375-9601(85)90773-xMR799090Müller-HoissenF.Dimensionally continued Euler forms: Kaluza-Klein cosmology and dimensional reduction19863466510.1088/0264-9381/3/4/020HenriquesA. B.Higher-dimensional cosmological solutions with generalized Einstein equations198627762163310.1016/0550-3213(86)90460-8MR860892ShiraishiK.The Friedmann universe and compact internal spaces in higher-dimensional gravity theories198676132132410.1143/ptp.76.321MR864333ArikM.The analogue potential of Kaluza-Klein cosmology with Euler form actions1989328130831610.1016/0550-3213(89)90106-5MR1028024DeruelleN.Fariña-BustoL.Lovelock gravitational field equations in cosmology199041123696370810.1103/physrevd.41.3696MR10636292-s2.0-0000405070MignemiS.Spontaneous compactification in 6-dimensional einstein-lanczos-maxwell theory19861533734310.1142/s0217732386000427MR880825MaedaK.Stability and attractor in a higher-dimensional cosmology. I19863223324710.1088/0264-9381/3/2/017AdeP. A. R.AghanimN.ArnaudM.Planck 2015 results. XIII. Cosmological parametershttp://arxiv.org/abs/1502.01589Planck CollaborationPlanck 2015 results. XX. Constraints on inflationhttp://arxiv.org/abs/1502.02114LüH.PerkinsA.PopeC. N.StelleK. S.Black holes in higher derivative gravity20151141717160110.1103/physrevlett.114.171601LüH.PerkinsA.PopeC. N.StelleK. S.Black holes in D=4 higher-derivative gravity201530154501610.1142/s0217751x15450165KehagiasA.KounnasC.LüstD.RiottoA.Black hole solutions in R2 gravity201520155, article 14310.1007/JHEP05(2015)143BrownA. R.DahlenA.MasoumiA.Flux compactifications on (S2)N201490404501610.1103/physrevd.90.045016FreundP. G. O.RubinM. A.Dynamics of dimensional reduction198097223323510.1016/0370-2693(80)90590-0MR593842HalliwellJ. J.The quantum cosmology of Einstein-Maxwell theory in six dimensions1986266122824410.1016/0550-3213(86)90182-3MR828428ZhongY.LiX.-Z.Quantum cosmology of Einstein-Kalb-Ramond theory1990422712715ZhongY.LiX.-Z.The quantum cosmology of 4-form fields in d=11 spacetime1991151113118StarobinskyA. A.A new type of isotropic cosmological models without singularity19809119910210.1016/0370-2693(80)90670-x