#### Abstract

We construct the higher order terms of curvatures in Lagrangians of the scale factor in the *D*-dimensional Robertson-Walker metric, which are linear in the second derivative of the scale factor with respect to cosmic time. It is shown that they are composed of the Lovelock tensors at the first step; iterative construction yields arbitrarily high order terms. The relation to the former work on higher order gravity is discussed. Despite the absence of scalar degrees of freedom in cosmological models which come from our Lagrangian, it is shown that an inflationary behavior of the scale factor can be found. The application to the thick brane solutions is also studied.

#### 1. Introduction

It is well known that higher-derivative gravity has a scalar degree of freedom in general [1–4]. In cosmological models of higher-derivative gravity, the scalar mode is expected to play an important role [5–7]. On the other hand, some cases are also known that higher order terms in curvatures for a gravitational action do not affect cosmological development of a scale factor. For example, it is known that terms which consist of contraction of Weyl tensors in a gravitational Lagrangian do not change evolutional equations for a scale factor in a model with homogeneous and isotropic space. The other special combinations of curvatures are known. In the specific dimension, the Euler form as a Lagrangian does not produce the dynamics of gravity at all, because the action becomes a topological quantity in such a case.

The dimensionally continued Euler densities have also been studied [8–17], because of their relation to the effective Lagrangian of string theory, and are found to give no scalar mode since the second derivative of the metric disappears in the action if we perform integration by parts. The absence of scalar modes is interesting for studying black holes in the theory, because the scalar modes lead to singularities, in general, which avoid expected horizons.

In recent years, it turned out that there is a special case where a scalar mode disappears in higher-derivative gravity. Originally, this fact was found in research of a three-dimensional theory of massive gravity [18], and an extended version in four dimensions was proposed [19]. The authors of those papers intended to study the renormalizability and unitarity of gravitation theory in a maximally symmetric spacetime. Thus, the absence of a massive scalar mode is at least a necessary condition of such theories referred as critical gravity (in our analysis, we do not care for the other critical values for the cosmological constant). Until now, however, only the cases with curvature tensors of a limited number have been investigated in higher dimensions [20, 21]. We are interested in the higher order theory of gravitation in which a scalar mode does not appear in a general higher dimensions.

In the present paper, we generalize the structure of the Lagrangian of critical gravity, to models with higher order terms in curvature tensors in higher dimensions. We show that such extensions can be attained by use of the Lovelock tensors. In order to offer a systematic way to construct the required higher order term, we take an explanatory approach by assuming the Robertson-Walker (RW) metric. In this approach, the absence of second derivatives of the scale factor from the Lagrangian with appropriate total derivatives is considered as a necessary condition for disappearing scalar modes.

It should be noted that the combination defined in -dimensional spacetime is used in critical gravities in three and four dimensions [18, 19]. This term can be considered as a trace of multiplication of the Einstein tensor and a linear combination of the Einstein tensor and its trace part: where is the Einstein tensor. Incidentally, is known as the Schouten tensor up to the factor . Now, if the RW metric is assumed, the time-time component of the Einstein tensor does not include the second time-derivative of the scale factor. The Schouten tensor appearing in the pervious term is made for its spatial component to have no second time-derivative of the scale factor. Therefore, the trace of the product of two tensors is linear in the second time-derivative of the scale factor, and if a surface term is suitably assigned, the Lagrangian is expressed only with the scale factor and its first time-derivative. Thus, additional scalar modes do not appear. From this observation, we find that it can be extended by using the Lovelock tensor instead of the Einstein tensor when the dimension of spacetime is higher. In the present paper, we do not analyze the massive tensor modes in our models. Thus, the genuine criticality as quantum gravity is left for future works.

The RW geometry is known to be conformally flat [22]; that is, the Weyl tensor for the RW cosmological metrics vanishes. An extension of Lovelock gravity for conformally flat geometry was considered by Meissner and Olechowski [23, 24]. They showed that the extension is possible for a term in dimensions provided , whereas the Lovelock gravity has the term at most (where denotes a general curvature tensor). Oliva and Ray also construct higher-derivative gravity with the second order equation utilizing the Weyl and Riemann tensors [25]. Their Lagrangian involves up to term with because of the use of the generalized Kronecker delta as in the case of the Lovelock gravity.

In the present paper, we show that it is possible to continue the higher-curvature terms beyond the number of dimensions.

The outline of this paper is as follows. In Section 2, we construct the candidate Lagrangian for cosmological models without scalar modes written in curvature tensors as the first step. The confirmation of the property of the model Lagrangian is performed in Section 3, substituting the RW metric. In Section 4, we show that an extension to more higher order terms in curvatures can be obtained. Using the higher order Lagrangian so far obtained, we propose (toy) models for the scale factor with inflationary behavior Section 5. In Section 6, the application of our Lagrangian to constructing domain wall solutions is studied. In the last section, we offer some concluding remarks.

#### 2. Lovelock Tensors and Generalization of the Higher Order Term in Critical Gravity

In this section, we construct the higher order term in curvatures by generalizing that of critical gravity. We will verify the absence of scalar modes in cosmological models with the terms in the next section. First, we introduce the dimensionally continued Euler density where the generalized Kronecker delta is defined as The dimensionally continued Euler density consists of th order in the curvature tensors (). For example, for , we find the Einstein-Hilbert term and for , we find the Gauss-Bonnet term as is well known.

Next, we consider the Lovelock tensor [8, 9]. The Lovelock tensor is a generalization of the Einstein tensor and defined as This is a symmetric tensor of th order in the curvature tensors. For example, for , we find and this is known as the Einstein tensor especially. For , we find It should be noted that where denotes the dimension of spacetime.

Here, we construct the new combination of the Lovelock tensor and the metric multiplied by the trace of the Lovelock tensor. That is, It is worth noting that For example, for , we obtain which is proportional to the Schouten tensor. Therefore, we obtain the following combination: which appears in critical gravities [18, 19].

Now, we find that the natural generalization of this is given by It is worth pointing out that the expression is symmetric against the exchange of and .

In the next section, we confirm that this combination is suitable for an extension of critical gravity in higher dimensions, by utilizing the RW metric.

#### 3. Higher Order Term for RW Metric

We consider the following RW metric in dimensions: where is the scale factor, and denotes the line element of a maximally symmetric space of -dimensions, whose scalar curvature is normalized to with .

First in this section, we examine the Lovelock tensors of . By explicit calculation of curvatures, we find, for , where and . Here and hereafter, we use the suffixes denoting the spatial dimensions, . Also, for , we obtain

By the combinatorial property of the generalized Kronecker delta, we can find the Lovelock tensor for a general as follows: and . It should be noted that the component of the Lovelock tensor does not include .

Next, we calculate the generalized Schouten tensor for the RW metric. Because the trace of the Lovelock tensor is given as we find By construction, is absent in .

Now, we consider the combined term . We obtain Because this combination is apparently linear in , the action including this term can be expressed as the functional of and , by means of integrations by part. Therefore, we realize that there is no scalar mode in a cosmological setting.

The combination can exist for limited numbers , which depends on . For the exchange symmetry, we assume that . We find that is nonzero when for , for , and for , and so forth. Therefore, we have new higher-derivative terms including with for odd and for even .

Up to now, we become aware of a similarity to the Lovelock Lagrangian. It takes the form for the RW metric as follows: The equivalence up to the overall constant is obvious; that is, . Incidentally, we can consider as the Lovelock Lagrangian . We note, however, that the Lovelock Lagrangian has its meaning as a part of Lagrangian when for , for , and for , and so forth, because becomes a total derivative for .

Before closing this section, discussion on relation to the work of Meissner and Olechowski [23, 24] is in order. Their approach is equivalent to considering the Lagrangian constructed from the Lovelock Lagrangian in which the Riemann tensor is expressed by the Schouten tensor under the assumption of vanishing Weyl tensor (strictly speaking, they replace Riemann tensor by the Kulkarni-Nomizu product of the Schouten tensor with the metric in the Lovelock tensor). They used the combinatorial property of the generalized Kronecker delta for extension to higher order in curvatures. If the RW metric is substituted, we find that the Lagrangian of Meissner and Olechowski of order of coincides with our Lagrangian , where , up to the constant factor. The allowed spacetime dimension is the same for odd , . Therefore, our Lagrangian and theirs are almost equivalent. The variety with respect to two integers in is due to the use of Riemann tensors as well as scalar curvatures and Ricci tensors in our approach. It is notable that differences may occur if we consider the black hole or nonconformally flat solutions in the theory governed by the Lagrangians of higher order terms.

Later, we show that the restriction by the dimensions can be overcome. To exhibit the discussion on the subject, we examine the cosmological action in the present model again in the next section.

#### 4. Further Extension to Higher Order in Curvature (Especially for )

We consider the action for : with arbitrary coefficients . Here, we regard as .

Then, the action for the scale factor can be read as where

Using the expansion we rewrite the action, after partial integration, as Especially for , we find a simple action

Now, we define the Lagrangian for the scale factor by Then, we find the equation of motion for the scale factor

The Hamiltonian constraint is regarded as the variation equation , where is the lapse function defined as . We set after the manipulation. We now find the Hamiltonian constraint equation

It is known that the variation of the lapse function corresponds to the variation of , and the variation of the scale factor corresponds to the variation of , up to certain factors. We can find, from (32) and (31), the following generalized Lovelock tensor as the variation of the action with respect to the metric:

The generalized Lovelock tensor can be written as variation of the action, which is proportional to as seen from the construction, but with arbitrary coefficients in the definition of the action . In spite of the arbitrariness, the functional form of the Lagrangian is unambiguous if the conformally flat metric is substituted.

In the similar manner, the corresponding generalized Schouten tensor is defined as in the previous section, with . Then, we find Then, the trace of the product of these tensors is found to be

Again, we obtain the same functional form as those constructed from the Lovelock tensors in the previous section. More iterative operations yield the more higher order terms. For this time, however, there is no limitation on the relation between the number of dimensions and the order of curvatures , except for the case with (the Lagrangian becomes a total derivative in this case for the conformally flat spacetime). Although we do not exhibit an explicit form of the higher order term in mass dimensions, obviously it can be expressed as the combination of curvature tensors and their covariant derivatives.

Unfortunately, only for , our iterative method cannot give term, since the term becomes a total derivative for the conformally flat metric. We define the expression for by the manner followed by Meissner and Olechowski as Incidentally, it is obvious that is at most linear in if the RW metric is substituted, because only is the component of the Schouten tensor which is linear in . Then, the desired terms of all the order in curvatures and derivatives in any dimensions can be created by our iterative method.

Another way to cross the dimensional limitation, which is inspired by the work of Meissner and Olechowski [23, 24], is also found. In the manner of their paper [23, 24], the following Lagrangian (modulo volume form) is proposed: where it is expressed by our present notation. We become aware of an extension of it as where are arbitrary integers. For the RW metric, this term turns out to be proportional to where . This is also the same form of the candidate Lagrangian. Many equivalent combinations exist for higher order terms under the vanishing-Weyl-tensor condition.

In the next section and after, we will turn our attention to apply the higher order Lagrangian obtained here to cosmological models.

#### 5. Cosmology

In this section, we investigate the possible inflationary stage in evolution of the universe in the model with higher order terms discussed in the previous sections.

Let us consider the general cosmological action with an appropriate total derivative term with respect to the time, which removes the second derivative of the scale factor with respect to time.

As usual, the energy momentum tensor of matter is taken as with the action for matter . The nonzero components of the energy momentum tensor of matter are assumed as and , where is the energy density and is the pressure, for the RW universe. As seen in the previous sections, the equation of motion derived from the action (42) coincides with the linear combination of the components of the Lovelock tensors in functional form with the assumption of the RW metric. Thus, the energy conservation holds generally. This fact is due to the absence of the scalar mode and the fact that we did not need to rescale the metric.

We here give a few models in four-dimensional spacetime with focus on the possibility of an inflationary growth of the scale factor. Furthermore, we take , that is, assume the flat space. Then, the action (42) is equivalent to The 00 component of the equation of motion reads

If we can choose the coefficient freely, almost arbitrary equations including and can be made. That is, where is a function which can be expressed by a series and does not include the term. Here, the factor given in the right hand side is chosen as for the similarity with the standard cosmology. We wish to call the cosmology of the model “ cosmology” (if , we obtain “ cosmology”).

Now, we specify the function . For example, let us consider with a typical mass scale . Solving the equation for , we obtain

If the energy density is sufficiently low, such as , the relation in the standard cosmology is valid: On the other hand, in the era of the high energy density such that , we find and we obtain an approximate de Sitter inflation (although the analysis here holds for any equation of state of matter, physical constraints may arise from the possible ranges of the energy density and the pressure of specific matter).

Another model can be chosen, in which the correction is a monomial of the higher order term. That is, In the rapid expansion phase, , the correction would be dominant. Further, if we assume the dust matter, that is, , we obtain This indicates the phase of the power-law inflation.

Even though we show special toy models here, we find that the existence of higher order terms yields the inflationary growth of the scale factor with ordinary matter, and with no scalar mode and no redefinition of the metric.

#### 6. Domain Wall

We will apply the previous discussion on a conformally flat metric to solutions for domain walls (“thick branes”). The -dimensional metric suitable for a codimension one domain wall is given by where . In addition, we consider a neutral, minimally coupled, self-interacting scalar field as a classical matter field. It is known that superpotential method [26–30] is available in this case to find BPS kink equations. Let us investigate the case with higher order terms in a similar manner.

Substituting the metric (54), the gravitational Lagrangian becomes where the prime () denotes a derivative with respect to . The term in can be removed by discarding a total divergence. Thus, we obtain The action for the real scalar field for static branes is written by Since the coefficients can be arbitrarily chosen, we use an arbitrary function to represent the general action. Then, we can write the total action as

The field equations can now be derived in a usual manner. The differential equation for the scalar is where . The equation for is written as where and . The reparametrization invariance of leads to the first integral given by

Define . Then, one can find .

Note that as well as does not possess the term for even . Using these, the equations can be rewritten as

Now, we take a BPS ansatz [26–31] Then, the equations reduce to the first order equation and the potential in terms of .

If the function is a polynomial up to the quadratic order, a domain wall solution exists, for the potential has two minima [26–31]. The solution can be expressed as a kink solution in the theory. We find here that the potential with many vacua naturally corresponds to the general higher order gravity.

For example, we try to express the sine-Gordon equation. We take , with a positive constant . Further, we choose with constants and . This choice is possible if the spacetime dimension is odd. Then, the scalar obeys the sine-Gordon equation and an exact static solution is known as Then, the potential takes the form: The minima of the potential are located at . Although it is difficult to obtain exact solutions of other types, we can suppose multiple domain walls with distinct topological numbers in this and similar models. The model with many vacua may also serve an interesting mechanism to realize naturally a small cosmological constant in (thick) brane world. This possibility will be studied in future.

#### 7. Summary and Conclusion

In the present paper, we have attempted to show the possible quasilinear second-order theory of gravity in conformally flat spacetimes. Models with arbitrary higher order of curvatures have been obtained. As long as we adopt an isotropic and homogeneous cosmological setting, the energy conservation holds in the models because there is no scalar mode and no requirement of frame rescaling. In spite of them, inflationary expansion can be found in the models. We have also found that the domain wall solution in the present type of the higher order gravity can be obtained naturally with the potential having many minima.

Our work corresponds to the extension of the Lovelock higher-curvature gravity in arbitrarily higher order terms in higher dimensions. Our analysis, however, has been limited for the case with the conformally flat spacetime, and the explicit form of the action written in curvature tensors and derivatives has been still ambiguous. The stability of the solutions obtained here is problematic for anisotropic perturbations or tensor modes. To study it, we should classify the possible form of the Lagrangian when the Weyl tensor vanishes. It is known that the dimensionally continued Euler forms have the property of factorization in terms of those of lower orders when the spacetime is represented by the direct product of spaces [12, 13]. Thus, compactification should be worth studying to seek some special combination of curvature tensors. It is expected that the hopeful “critical” relation among the coefficients of different orders of curvatures will be selected by consideration on various background spacetimes.

Nevertheless, we emphasize that our arbitrarily high order gravity can be applied to many models in various contexts. The problem of initial singularity can be reconsidered by studying classical bouncing universes in our model. On the other hand, the Wheeler-DeWitt equation of the Lagrangian should lead to higher-derivative quantum cosmology. The equation must be difficult to treat with, but the study on it may shed new light on quantum gravity. The possible black hole solutions are interesting in both cases of asymptotically flat and asymptotically AdS spacetime. We will return to some of the problems in future.