Cosmological model due to dimensional reduction of higher-dimensional massive gravity theory

We investigate a cosmological model resulting from a dimensional reduction of the higher-dimensional dRGT massive gravity. By using the Kaluza-Klein dimensional reduction, we obtain an effective four-dimensional massive gravity theory with a scalar field. It is found that the resulting theory corresponds to a combined description of mass-varying massive gravity and quasi-dilaton massive gravity. By analyzing the cosmological solution, we found that it is possible to obtain the late-time expansion of the universe due to the graviton mass. By using a dynamical system approach, we found regions of model parameters for which the late-time expansion of the universe is a stable fixed point. Moreover, this also provides a mechanism to stabilize the extra dimensions.


I. INTRODUCTION
It is well-known that the universe is expanding with acceleration nowadays [1,2]. By using general relativity (GR) with cosmological constant, the theoretical predictions satisfy most of observational results. However, the value of the cosmological constant needs to be fine-tuned to the order of 10 −122 in Planck units. This does not provide a proper model in theoretical point of view. As a result, alternatively theoretical models have been intensively constructed in order to describe the accelerating expansion of the universe. One of possible ways to construct this theoretical model is to modify GR at large scale, driving the expansion at cosmological scale while recovering GR at local gravity scale. In this work, we focus on the modification of GR by adding suitable mass term into Einstein-Hilbert action proposed by de Rham, Gabadadze and Tolley [3,4] namely dRGT massive gravity theory.
In dRGT massive gravity theory, the mass terms are systematically constructed in order to eliminated the Boulware-Deser (BD) ghost from the theory. The BD ghost is an additional sixth degree of freedom which usually emerges in nonlinear massive gravity theory [5]. This healthy nonlinear massive gravity then provides a new window to various researches, not only in cosmological context but also in astronomical objects and black holes [6][7][8][9][10]. In cosmological context, it was found that the theory with Minkowski fiducial metric does not admit flat and closed Friedmann-Laîrmatre-Robertson-Walker (FLRW) solutions [11]. In order to obtain all kinds of FLRW solutions, it has been suggested to consider three options as follows [12]: 1) Use other forms of fiducial metric, for example, FLRW or de-Sitter [13][14][15][16], 2) Consider anisotropic universe [17][18][19] and 3) Add more degree of freedom into the theory. Although the self-accelerating expansion of the universe can be provided by using the FLRW fiducial metric, there are only two propagating degrees of freedom while there should be five for massive gravity theory in four-dimensional spacetime [13]. Adding more degrees of freedom is still one of possible ways for solving the lack of the number of propagating degrees of freedom. In this work, we will use the flat FLRW fiducial metric and specialize our attention to the third option where the additional degree of freedom is a scalar field.
The dRGT massive gravity models with additional degree of freedom have been investigated in various ways. By promoting the fiducial metric to be dynamical, the theory is known as "bi-gravity theory" [20]. Other class of this extension is obtained by adding a scalar field into the theory. One of simple arguments is that promoting the graviton mass to be varied according to a scalar field called mass-varying massive gravity (MVMG) [21][22][23][24][25].
One of drawback features of this model is that the accelerating expansion of the universe is driven by the scalar field instead of the graviton mass, the graviton mass shrinks to zero at late time expansion of the universe. Fortunately, the graviton mass can play a role of the cosmological constant by using the FLRW form of the fiducial metric [26]. Further investigation also shows that a part of graviton mass can be responsible for dark matter by promoting the scalar field as the k-essence field [26].
Another extension of dRGT massive gravity by using scalar field is quasi-dilaton massive gravity (QMG) [27,28]. The additional scalar field possesses the quasi-dilaton global symmetry. It was found that the self-accelerating solution to the model plagued by ghost instability [29][30][31]. Further extended models of the QMG are investigated by introducing the coupling between the scalar field and the Stuckelberg fields [32,33] and then found that it is possible to avoid the ghost instability in some classes of the solutions [34][35][36][37][38][39][40][41][42]. Another extension of the theory is investigated by introducing the Dirac-Born-Infeld (DBI) scalar field in such a way that the scalar field possesses the generalized Galileon shift symmetry [43][44][45]. Unfortunately, the model encounters problems such that the number degree of freedom is not correct or the ghost mode emerges in the linear perturbation level [46].
The theories with higher-dimensional spacetime have been intensively motivated from a more fundamental theory such as String theory (or M theory). Since the observations suggest that we live in four-dimensional spacetime, a mechanism to reduce the higher-dimensional spacetime to four-dimensional spacetime must be involved and one of well-known mechanisms is the Kaluza-Klein (KK) dimensional reduction [47,48]. Since gravity theories in higher-dimensional spacetime have more degree of freedom than ones in four-dimensional spacetime, it is generally to obtain the four-dimensional gravity theory with some additional fields associated with the structure of the extra dimensions. It is important to note that cosmological aspect of dRGT massive gravity with higher-dimensional spacetime has been investigated in the context of braneworld scenario [49] where the matter fields are confined in three-dimensional hypersurface [50,51]. In this work, we focus on other kinds of the higher-dimensional gravity theory where the matter fields can access to all dimensions called Universal Extra Dimension (UED) [52]. Note also that besides a dRGT theory, a dimensional reduction of a partially massless theory, which is closely related to a linear massive gravity, has been investigated [53]. Due to the simple structure of the extra dimensions, the additional field is a scalar field characterizing the radius of the extra dimensions. Therefore, it is worthwhile to consider the dRGT massive gravity theory in higher-dimensional   spacetime and choose a simple ansatz, assuming the maximally symmetric extra dimensions, to obtain a four-dimensional dRGT massive gravity theory with an additional scalar field.
According to this argument, we propose an alternative extension of dRGT massive gravity by using a scalar field which is obtained from the scalar degree of freedom characterizing the radius of the volume of the extra dimensions and then find a possibility to provide a cosmological solution such that the late-time expansion of the universe is obtained. This paper is organized as follows: in Section II, we consider dRGT massive gravity theory in (4+d)-dimensional spacetime and then perform the KK dimensional reduction.
As a result, we obtain an effective four-dimensional massive gravity theory with a scalar field. It is found that the theory admits a combined description of MVMG and QMG. In Section III, we investigate a cosmological model of this kind of massive gravity theory and found that it is possible to obtain the late-time expansion of the universe. In Section IV, we investigate this possibility by using the dynamical system approach to analyze dynamics of the universe. We found some regions of the model parameters that provide the standard evolution of the universe as well as the extra dimension can be stabilized. Finally, we summarize and discuss the results in Section V.

GRAVITY THEORY
One of possible ways to construct a consistent theory of massive gravity is to consider the Einstein gravity theory in (4 + d)-dimensional spacetime where d is the number of the extra dimensions. In fact it was shown that the dRGT massive gravity theory can be obtained by using discretization procedure of pure massless gravity in (4 + d)-dimensional spacetime [54][55][56]. As we have mentioned, one of possible ways to obtain the cosmological solution in dRGT massive gravity is adding the additional field into the theory. Since the scalar degree of freedom can be obtained from dimensional reduction, it is worthwhile to consider a reduced theory of the (4 + d)-dimensional dRGT massive gravity theory. In this section, we will perform a dimensional reduction of the (4 + d)-dimensional dRGT massive gravity theory. Let us consider the dRGT massive gravity theory in (4 + d)-dimensional spacetime as follows [57] where quantities with tilde,X, refer to the quantities in (4 + d)-dimensional spacetime.
R is the Ricci scalar, M (4+d) is the Planck mass in (4 + d)-dimensional spacetime, m g is the parameter associated to the graviton mass. Each potentialŨ i is a function of [K n ] The explicit form ofŨ i is presented in Appendix A. The metricf AB is the reference or fiducial metric and capital Latin indices refer to (4 + d)-dimensional spacetime indices. In order to perform the KK dimensional reduction which keeps only the scalar field left, we neglect the cross components between the extra dimension and the four-dimensional spacetime. As a result, the physical metric can be written asg where g µν is the metric in four-dimensional spacetime, γ ab is the metric in the extra dimensions, the small Latin indices refer to the indices in extra dimensions, and p(φ) is a function playing a role of the radius of the extra dimensions. In the same way as the physical metric, the fiducial metric in (4 + d)-dimensional spacetime can be decomposed as where f µν is the fiducial metric in four-dimensional spacetime and q(φ) is a function playing the role of the radius of the extra dimensions in the fiducial sector. For this ansatz, the action in Eq. (1) can be rewritten as where y a are coordinates on the extra-dimensional manifold which is assumed to be maximally symmetric with the metric γ ab . R[γ] is the Ricci scalar associated with the metric γ ab . Since we deal with the maximally symmetric (extra) space, the curvature parameter of γ ab for d > 1 is expressed as κ = R[γ] d(d−1) . Note that this curvature parameter for d = 1 is undefined and the geometry of the extra dimensions is hyperbolic, flat and spherical for κ < 0, κ = 0 and κ > 0 respectively. The potential in the mass termŨ can also be splitted into the usual form of the four-dimensional one, U, and the additional term, F . The usual potential, U is taken the same form asŨ but put For the additional potential, F , can be written as F = F 2 + 4+d j=3 α j F j where r = 1 − q/p. Each F i is written explicitly in Appendix B. Now we will perform the dimensional reduction of the (4 + d)-dimensional action in Eq. (4). As a result, the effective four-dimensional action becomes where the four-dimensional Planck mass, m pl , is related to the (4 From this scalar potential, one can see that it vanishes for the case of five-dimensional spacetime or if the spacetime geometry of the extra dimension is flat. The canonical form of the scalar field can be obtained by setting p = e 2 d(d+2) φ m pl . As we have seen, the graviton mass, M g is promoted as a function of the scalar field, Note that the radius function of the fiducial metric, q(φ), is arbitrary. For the simplest case, one can choose as the same form as the radius function of the physical metric. This corresponds to choosing r = 0 so that the potential reduces to one in four-dimensional massive gravity. Moreover, the potential U also respects a global symmetry where ψ a are four Stuckelberg fields defined via the fiducial metric as Therefore, the structure of the potential is exactly the same with QMG [27,28]. However, the mass function, M g (φ), breaks this global symmetry. Since the mass function depends on the scalar field, the theory is equivalent to MVMG [21][22][23][24][25]. As a result, the theory is somewhat a combination of QMG and MVMG. Up to our knowledge, this kind of extended dRGT massive gravity has not been investigated yet. Therefore, the main purpose of this paper is to investigate the cosmological models due to this kind of massive gravity.
Since the theory contains two dynamical fields, g µν and φ, equations of motion can be obtained by varying the action in Eq. (6) with respect to both fields. As a result, the equations of motion for g µν can be expressed as where X µν = δU δg µν − 1 2 g µν U and Y µν = δF δg µν − 1 2 g µν F . The explicit forms are lengthy, they will be presented in Appendix C. It is convenient to rewrite Eq. (9) in terms of the energy momentum tensors of the scalar field and the mass term as follows where By using the Bianchi identity, the conservation of the total energy momentum tensor can be obtained as follows Note that since both contents are coupled, they are not separately conserved. Actually, it can be viewed as a coupling model between a scalar field and the content corresponding to the graviton mass. We will show this argument explicitly in the next section. Note also that this conserved equation corresponds to the equation obtained from varying the action with respect to the fiducial metric. This will provide a constraint equation since the fiducial metric plays the role of Lagrange multiplier in the language of Lagrangian formulation.
The equation of motion corresponding to varying the action in Eq. (6) with respect to the scalar field can be written as where Φ = m pl In order to investigate the cosmological solutions, let us adopt the physical metric as the flat FLRW metric as follows where a(t) is a scale factor determining the scale of the spatial distance. Moreover, it is convenient to work in the unitary gauge where ψ a = x a and then choose the form of the fiducial metric as the flat FLRW metric, where b(t) is a scale factor for the spatial distance in fiducial metric. It is instructive and worthwhile to consider the simple case by choosing the radius function as q = p, corresponding to r = 0. This is instructive but useful for a more general case. Including Lagrangian density of the matter L m and radiation L r , the action of the model can be written as As a result, the modified Einstein field equation in Eq. (10) becomes By using ansatz we adopt in this section, the components (0, 0) and (i, j) can be written respectively as where H is the Hubble parameter, ρ m and p m are energy density and pressure of the matter respectively as well as ρ r and p r are ones of the radiation. The short-hand function A and B obtained from the tensor X µ ν can be expressed as For the scalar field, the equation of motion in Eq. (13) becomes In the same fashion, the constraint equation in Eq. (12) can be expressed as By using Eq. (21) and then rewrite the equation in a convenient form, we have where the over dot refers to the derivative with respect to t. From this constraint equation, an interesting branch of the solution is obtained A = B. For this branch, the equation of state parameter, w g , is exactly equal to −1 and then the graviton mass will play the role of cosmological constant. Note that the energy density and pressure contributed from graviton mass can be respectively written as Considering the A = B branch, called self-accelerating branch, one can see that the nontrivial solution isφ/H = 0. This means that the graviton mass, M g , is constant. This behaviour cannot be obtained for MVMG since the graviton mass always decrease at late time.
Moreover, it is important to note that the equations of motion can be viewed as a coupling equations between a scalar field and the graviton mass. By rewriting equations of motion, the coupling equations can be written as where and prime denotes the derivative with respect to N = ln a. From these equations, it is found that the interaction term, which actually corresponds to the energy transfer between two contents, vanishes when φ ′ = 0. This is consistent with the self-accelerating branch as we have discussed above. In the case of Φ = 0, the energy transfer between graviton mass and the scalar field also vanishes, but it is not necessary to obtain the solution which predicts the accelerating expansion of the universe because the solutions of y which satisfy the case it does not correspond to the fixed point, then the scalar field still evolve. We will see these behavior more clearly when we consider the dynamical system and this is a main issue for the next section.

IV. DYNAMICAL SYSTEM
To see the behavior of each cosmological contents for the model of this effective massive gravity theory in more detail, we will analyze the model by using a dynamical system approach. As a result in the previous section, (0,0) components of Einstein equation in Eq.
(18) can be written as where This is one of the constraint equations which contains six variables including y. The others come from Eq. (23) which is expressed in terms of dynamical variables as It is important to note that the dynamics of the universe with the additional scalar field The third constraint equation is the consequence of the definition of the potential of the Note that we have introduced a new parameter γ to characterize an effect of the potential term compared to the graviton mass. Actually, this parameter tells us how the curvature of the extra dimensions affects the dynamics of four-dimensional universe compared to the effect of the graviton mass. As a result, the geometry of the internal space can be characterized by this parameter as follows; the geometry of the internal space is hyperbolic, flat and spherical for γ < 0, γ = 0 and γ > 0 respectively. Note also that for d = 1, there is no internal curvature so that the potential term vanishes automatically.
Now we can choose to eliminate three dynamical variables, z, v and Ω m for convenience.
As a result, three dynamical equations for three variables can be written as where As we have mentioned, these fixed points correspond to self-accelerating expansion of the universe. Note also that x ′ = 0 in Eq. (33) is satisfied due toḢ = 0. Moreover, from Eq. (34), y ′ depends only on y so that one can examine the stability of the fixed points separately. As a result, the stability condition can be found by From this condition, one can see that the quantity in the left hand side depends on three parameters; α, β and s. One of them can be eliminated by using the initial conditions such as Ω m0 ∼ 0.25, Ω r0 ∼ 0 and Ω g0 = −Ax/3 + v ∼ 0.75. This can be done only in fivedimensional model where v = 0. For the model with higher than five dimensions, we have one more parameter, γ. Therefore, we will separate our analysis into two parts, five-dimensional model and higher than five-dimensional model.

A. Five-dimensional model
For five-dimensional model, the scalar potential vanishes then the constraint in Eq. (28) becomes −Ax/3 + Ω m0 + Ω r0 = 1. As we have mentioned above, one can find the stability condition in terms of parameter (α, β) by eliminating s due to the initial conditions such that Ω m0 ∼ 0.25, Ω r0 ∼ 0 and Ω g0 = −Ax/3 ∼ 0.75. By using a numerical method, the regions satisfying the stability condition for y + and y − can be illustrated as in Fig. 1.
It is important to note that the stability condition in Eq. (38) can be rewritten explicitly Since solutions y ± are solved from Y = 0, ∂ y Y | y + is positive and ∂ y Y | y − is negative. Therefore, the stability condition can be inferred from the sign of 2A − Φ; 2A − Φ < 0 for y + and  There is a special case where ∂ y Y | y ∓ = 0. In this case, the eigenvalue is equal to zero and the stability cannot be inferred from linear analysis. Since this point corresponds to the minimum of Y , this guarantees the stability of the nonlinear analysis. As a result, the condition for this point can be written as 3β = α 2 which turns out that y + = y − . Moreover, this point corresponds to the case z = 0 and z ′ = 0. This means that the scalar field is fixed all the time of the evolution of the universe. The radius of the extra dimension is also stabilized at this fixed point. By choosing α = 2, the evolution of the contents in the universe, Ω g , Ω m , Ω r and the w ef f can be numerically evaluated as shown in Fig. 2. From this figure, one can see that all contents evolve as in standard evolution such that there exist the radiation, matter and dark energy dominated periods.
Now we proceed our analysis to the case of ∂ y Y | y ∓ = 0. For this case, it is expected that the dynamics of the scalar field will affect the evolution of the universe since ∂ y z = 0 leading to z ′ = 0. For y − solution, the stability region in Fig 1 is a region below the line α 2 = 3β, α > 0.
Since w φ = 1, so that if z takes the value of order 0.1, w ef f will significantly deviate from the standard evolution. In order to obtain such a small effect of the kinetic of the scalar field, one has to take the parameter slightly from α 2 = 3β line. As a result, we choose α = 2, β = 1.32 and then the evolution of the contents in the universe is shown in Fig. 3.
From this figure, one found that there is a peak of z about z ∼ 0.001 at log a ∼ −2. This peak takes a small value due to choosing parameters α, β as discussed above. Note that the peak can be eliminated by choosing small value of α, β, for example α = 0.5, β = 0.0825, leading to small value of ∂z/∂y. Though it is not shown here, we demonstrate that we can obtain the standard evolution by setting as discussed.
For y + solution, the stability region is larger than y − solution. A behavior of the density parameters for each content in the universe is similar to ones in y − solution. In order to obtain the standard evolution of the density parameters − z does not significantly emerge during the whole predicted evolution − one has to choose the parameter α, β, for example α = −2.0, β = 0.5, such that ∂z/∂y −0.1. Since the evolution is similar to y − solution, we do not show it explicitly in this presentation.
It is important to note that z may dominate at early times since ∂z/∂y = 0. In this case, one has to tune parameters such that ∂z/∂y −0.01. For this condition, z still take the value less than 10 −4 at log a ∼ −50. This is enough to guarantee that the change of the radius of the extra dimension does not affect the Big Bang Nucleosynthesis.

B. Six-dimensional model
In this subsection, we will consider the dynamics of the universe including the potential term of the scalar field. As we have mentioned, the effect of the potential term can be characterized by parameter γ corresponding to the curvature of the internal space. One of important issues in this model is that the potential term is not involved in the dynamical system explicitly and then its dynamics is completely determined since it can be written in terms of y. Therefore, the effect of this term is concerned mostly in the region of stability of the fixed point.
As we have discussed in five dimensions case, there exists a special class of the solution such that y − = y + (or in terms of parameters β = α 2 /3). Therefore, it is instructive to consider this class of the solution to find out how the parameter affect the dynamics. As a result, the parameter s which determined from the initial value can be expressed as In order to obtain the real value of s, the parameter γ is restricted by From this condition, one finds that γ cannot take very large positive value without tuning parameter α. Actually, for order unity of α, γ can take order of 10 −2 . This will make the stability region narrow for positive value γ. On the other hand, γ can take very large negative value. This provide larger region of stability than one for γ = 0. The stability regions for y − solution with γ = 0.01 and γ = −0.1 are shown in Fig. 4. From this figure, one can see that for negative γ, the stability region is enlarged while the region is reduced for positive γ. Note that for y + solution, the change of the region due to the existence of γ is similar to one for y − . We do not show explicitly in this presentation. For the evolution of the density parameters, the behavior can be seen by using Eq. (40). From this equation, it is found that the existence of γ can make s decreased and ∂ y z ∝ (s − 1) 2 . Therefore the dynamics of z is reduced and then the peak reduces and moves to the earlier time. This behavior is confirmed by using numerical calculation as shown in the left panel of Fig. 5. By comparing this figure to the right panel of Fig. 3, one can see that the existence of potential term will reduce the peak of z and shift the peak to earlier time. We also found this behavior in both y + and y − solutions. However, it is not obvious in the case of y + since the stability region takes place in α < 0 region. One more feature of the dynamics in six-dimensions is that the contribution to drive the late-time expansion of the universe does not only come from the graviton mass but also from the potential term of scalar field. In the right panel of Fig. 5, the evolution of the contents in the universe is illustrated including the contribution from scalar potential, v (dashed-red line). In this case, we found that the expansion of the universe nowadays is driven by the graviton mass and the curvature of the extra dimension.
Note that for the case s = 1, y will be a non-dynamical variable then the scalar field is just a constant in the dynamical system. Therefore, this case gives us the same prediction as GR with the cosmological constant.
It is worthwhile to note that one of the problem from the cosmological model reduced from the higher dimensional space time is that it is difficult to find the mechanism to stabilize the extra dimensions. For example, in the Casimir dark energy models, one have to add other structure, Aether field [58] or Gauss-Bonnet terms [59], to stabilize the extra dimensions.
In our case, without the graviton mass term, the scalar potential is the exponentially decay so that the scalar field roll the potential down to the very large value. Ultimately, the extra dimensions cannot be compactified. However, by including the graviton mass terms, we found that the extra dimensions can be stabilized so that it is possible to compactified the extra dimensions to small value for hiding it from the observations. This is one of important result of the model. are also interpreted as the function of the scalar field. We found that this effective massive gravity theory in four-dimensional spacetime contain description of two extensions of dRGT theory, namely mass-varying massive gravity (MVMG) and quasidilaton massive gravity (QMG). Up to our knowledge, this kind of extended dRGT massive gravity has not been investigated yet. Therefore, we then investigate the cosmological model due to this effective massive gravity theory.
In order to investigate the cosmological solutions from our effective massive gravity theory, we adopt the form of physical metric and fiducial metric as the FLRW metric form. In order to complete our analysis, matter and radiation are introduced into our cosmological model.
It is found that the graviton mass and the scalar field are coupled. We also found that it is possible to obtain solutions providing the self-accelerating expansion of the universe (the case A = B). We pay attention to this possibility in detail by using the dynamical system approach.
To analyze dynamics of these cosmological models, the six dynamical variables x, y, z, v, Ω m and Ω r with four free parameters, α, β, s and γ are defined. For the five-dimensional model, in the special case where y − = y + , the kinetic term of the scalar field, or variable z, vanishes exactly. The free parameters α and β can be chosen without regarding the effect of z. However, the effect of z must be small for the case of y − = y + in order to obtain the expected evolution of the universe.
For the higher-dimensional model which is studied as six dimensions in this work, the potential of the scalar field characterized by the free parameter γ is introduced. We found that dynamics of z is affected by this free parameter through the ratio of two scale factors, s and the existence of the scalar potential provides us a larger stability region of α and β. The appropriated free parameters which give the standard evolution of the universe for each part of consideration have already been discussed (see in our Section IV). The crucial difference between the five-dimensional model and the higher one is the late-time expansion of the universe is obtained due to only the graviton mass in five dimensional spacetime while the both of graviton mass and the potential of the scalar field can drive the accelerating expansion of the universe. Moreover, by including the graviton mass term, we found that the extra dimensions can be stabilized. Therefore, it is possible to compactified the extra dimensions to small value for hiding it from the observations while it is not possible for the case such that the graviton mass is not introduced. This may shed light on the interplay between the consistent cosmological model and the fundamental theory with higher-dimensional spacetime.
Note that this effective massive gravity theory is not studied in the general case which is the case of r = 0 or p = q. Therefore, the cosmological implication in more complicated case is still the interesting challenge. The instability issue as well as the number of degree of freedom of the theory is one of interesting issues to investigate. We leave this investigation for further work.
In our investigation, the matter is not coupled to the scalar field. However, in principle, it is possible that they are coupled together since the conformal factor will influence matter Lagrangian. This may lead to coupling model between dark energy contributed from graviton mass and the matter field. Therefore it may be possible to use this coupling model to examine the possible way to solve the coincidence problem. We leave this investigation for further work.

ACKNOWLEDGMENTS
We would like to thank Lunchakorn Tannukij for value discussion and reading through the manuscript. This project is partially supported by the ICTP through grant No. OEA-NT-01.

Appendix A: The explicit form ofŨ
To construct the massive gravity theory which eliminates BD ghost, each potentialŨ i is implied from the characteristic equation given i × i square matrixK,K i×i [57]. In this work, we useK denotesK i×i in which their determinants can be constructedŨ i by detK i×i = detK = 1 2(i)!Ũ i . The form ofŨ i in Eq.
(1) will be expressed up to i = 6 as The additional potential obtained from the dimensional reduction was defined in Eq. (4) where with F 2 = 2[K] + r(d − 1), , , , .