Radion stabilization in higher curvature warped spacetime

TL;DR

This paper explores how higher curvature terms in a five-dimensional warped spacetime can naturally stabilize the modulus (radion) field, affecting its mass and interactions, with implications for phenomenology.

Contribution

It demonstrates that a $F(R) = R + eta R^2$ gravity model in warped spacetime can inherently stabilize the radion without additional fields, providing a geometric stabilization mechanism.

Findings

Radion stabilization is achievable in higher curvature warped spacetime.

Higher curvature terms modify radion mass and coupling.

The model identifies parameter regions for stable radion.

Abstract

We consider a five dimensional AdS spacetime in presence of higher curvature term like $F(R) = R + \alpha R^2$ in the bulk. In this model, we examine the possibility of modulus stabilization from the scalar degrees of freedom of higher curvature gravity free of ghosts. Our result reveals that the model stabilizes itself and the mechanism of modulus stabilization can be argued from a geometric point of view. We determine the region of the parametric space for which the modulus (or radion) can to be stabilized. We also show how the mass and coupling parameters of radion field are modified due to higher curvature term leading to modifications of its phenomenological implications on the visible 3-brane.

Equations

$$d{s}^2=e^{-2k{r}_c|\varphi |}{\eta }_{\mu \nu }d{x}^{\mu }d{x}^{\nu }-r_c^{2}d{\varphi }^2$$

$$M_{Pl}^2=\frac{M^3}{k}(1-e^{-2k\pi r_c})$$

$$m_H=m_0{e}^{-k\pi r_c}$$

$$S=-\frac{1}{2{\kappa }^2}\int d^{4}xd\varphi \sqrt{G}F(R)$$

$$S=-\frac{1}{2{\kappa }^2}\int d^{4}xd\varphi \sqrt{G}[F^{\prime }(A)(R-A)+F(A)]$$

$$G_{MN}(x,\varphi )\to {\tilde{G}}_{MN}=\exp (\sigma (x,\varphi )G_{MN}(x,\varphi )$$

$$\displaystyle R=e^{\sigma}\bigg{[}\tilde{R}-4\tilde{G}^{MN}\partial_{M}\partial_{N}\sigma-3\tilde{G}^{MN}\partial_{M}\sigma\partial_{N}\sigma\bigg{]}$$

$$S$$

$$\displaystyle S=(1/2\kappa^{2})\int d^{4}xd\phi\sqrt{\tilde{G}}\bigg{[}-\tilde{R}+3\tilde{G}^{MN}\partial_{M}\sigma\partial_{N}\sigma$$

$$\displaystyle-\frac{1}{2\kappa^{2}}\bigg{(}\frac{F(A)}{F^{\prime}(A)^{5/3}}-\frac{A}{F^{\prime}(A)^{2/3}}\bigg{)}\bigg{]}$$

$$S=\int d^{4}xd\phi\sqrt{\tilde{G}}\bigg{[}-\frac{\tilde{R}}{2\kappa^{2}}+\frac{1}{2}\tilde{G}^{MN}\partial_{M}\Phi\partial_{N}\Phi-V(\Phi)\bigg{]}\\$$

$$S$$

$$G_{MN}(x,\varphi )\to {\tilde{G}}_{MN}=\exp (\frac{1}{\sqrt{3}}\kappa \Phi (x,\varphi ))G_{MN}(x,\varphi )$$

$$S$$

$$V(\Phi )$$

$$\exp (\frac{3}{2\sqrt{3}}\kappa <\Phi >)=[\sqrt{9-40{\kappa }^2\alpha \Lambda }-2]$$

$$m_{\Phi }^2=-\frac{1}{8\alpha }[\sqrt{9-40{\kappa }^2\alpha \Lambda }][\sqrt{9-40{\kappa }^2\alpha \Lambda }-2{]}^{-\frac{2}{3}}$$

$$V(<\Phi >)$$

$$S$$

$$d{\tilde{s}}^2=e^{-2k{r}_c|\varphi |}{\eta }_{\mu \nu }d{x}^{\mu }d{x}^{\nu }-r_c^{2}d{\varphi }^2$$

$$-$$

$$\xi(\phi)=e^{2kr_{c}|\phi|}\big{[}Ae^{\nu kr_{c}|\phi|}+Be^{-\nu kr_{c}|\phi|}\big{]}$$

$$A=v_v{e}^{-(2+\nu )k{r}_c\pi }-v_h{e}^{-2\nu k{r}_c\pi }$$

$$B=v_h(1+e^{-2\nu k{r}_c\pi })-v_v{e}^{-(2+\nu )k{r}_c\pi }$$

$$\xi (0)$$

$$\xi (\pi )$$

$$V_{eff}(r_c)=k(\nu +2)A^2(e^{2\nu k{r}_c\pi }-1)+k(\nu -2)B^2(1-e^{-2\nu k{r}_c\pi })$$

$$k\pi r_c=\frac{4k^2}{m_{\Phi }^2}\ln [\frac{v_h}{v_v}]$$

$$\xi (\varphi )=\frac{2}{\sqrt{3}\kappa }\ln [1+2\alpha R]-<\Phi >$$

$$v_{h}=\frac{2}{\sqrt{3}\kappa}\ln\bigg{[}\frac{1+2\alpha R(0)}{\sqrt{9-40\kappa^{2}\alpha\Lambda}-2}\bigg{]}$$

$$v_{v}=\frac{2}{\sqrt{3}\kappa}\ln\bigg{[}\frac{1+2\alpha R(\pi)}{\sqrt{9-40\kappa^{2}\alpha\Lambda}-2}\bigg{]}$$

Full text

I Introduction

Till date, Standard Model (SM) of particle Physics is a widely accepted theory to describe the interactions of fundamental particles. Despite its enormous successes, the model is plagued with divergence of the Higgs mass due to radiative corrections which may run up to Planck scale. An unnatural fine tuning is needed to confine the Higgs mass within TeV scale.

Many attempts have been made to address this problem by considering the theories beyond SM of particle Physics. Few such candidates are - supersymmetry, technicolor and extra dimensions. Among many such attempts arkani ; horava ; RS1 ; kaloper ; cohen ; burgess ; chodos , Randall-Sundrum (RS) model RS1 of warped extra dimension draws special attention since it resolves the gauge hierarchy problem without choosing any intermediate scale in the theory.

In RS model, two 3-branes are embedded in a five dimensional anti de-Sitter spacetime compactified on a $M_{4}\times S_{1}/Z_{2}$ orbifold. The distance between these two branes is assumed to be $\sim$ of Planck length so that the required hierarchy between the two branes is generated. One of the crucial aspects of this braneworld scenario is to stabilize the distance between the branes (known as modulus or radion). For this purpose, it is necessary to generate an appropriate radion potential with a stable minimum consistent with the value proposed in RS model in order to solve the hierarchy problem. Goldberger and Wise (GW) proposed a mechanism GW1 to create such a radion potential by introducing a bulk scalar field with appropriate boundary values. Subsequently the phenomenology of radion field has also been studied extensively. This radion phenomenology GW_radion ; kribs ; julien ; wolfe along with the study of RS graviton dhr ; rizzo ; yong ; dhr1 ; thomas and RS black holes wiseman ; naresh ; dai are considered to be the testing ground of warped extra dimensional models in collider experiments atlas1 ; atlas2 . As the present experimental lower bound of the first graviton KK-mode mass climbs above 3 Tev, the RS-type resolution of the naturalness problem is undoubtedly under pressure. The question that how the Higgs is so much lighter than the 5-dimensional Planck scale needs to be settled properly. However, in any higher dimensional model with gravity in the bulk, the modulus must be stabilized to appropriate value to extract a meaningful low energy effective theory on the brane. The present work aims to address this issue specially in the context of a higher dimensional model where the fundamental curvature scale is of the order of Planck scale.

It is well known that Einstein-Hilbert action can be generalized by adding higher order curvature terms which naturally arise from diffeomorphism property of the action. Such terms also have their origin in String theory due to quantum corrections. $F(R)$ faraoni ; felice ; paliathanasis , Gauss-Bonnet (GB) nojiri2 ; nojiri3 ; cognola or more generally Lanczos-Lovelock gravity are some of the candidates in higher curvature gravitational theory.

Higher curvature terms become extremely relevant at the regime of large curvature. Thus for RS bulk geometry where the curvature is of the order of Planck scale, the higher curvature terms should play a crucial role. In general inclusion of higher curvature terms in the action leads to the appearance of ghost from higher derivative terms resulting into Ostragradsky instability. The Gauss-Bonnet model ( a special case of Lanczos-Lovelock model ) is however free of this instability due to appropriate choice of various quadratic combinations of Riemann tensor, Ricci tensor and curvature scalar . Some important modified solutions of the Randall-Sundrum model in presence of Gauss-Bonnet terms have been obtained by Kim et.al kim1 ; kim2 in the context of both static and inflationary scenario. A GB modified warped solution and it’s phenomenological implications was also discussed in ssgsc .

In contrast to GB model $F(R)$ gravity model however contains higher curvature terms consisting only of the scalar curvature $R$. Once again just as GB model, certain classes of $F(R)$ gravity models are free from ghost-like instability ( See section III for a detailed discussion ). In general $F(R)$ model can be mapped into a scalar-tensor theory at the action level by a conformal transformation of the metric faraoni ; felice ; barrow ; marino ; bahamonde ; catena ; ssg1 ; ssg2 . The issue of instability of the original $F(R)$ model is now reflected in the form of the kinetic and potential terms of the scalar field in the dual scalar-tensor model, where the potential will have a stable minimum and a kinetic term with proper signature if the original $F(R)$ model is free from ghosts ( see section III ). It is know that to achieve modulus stabilization in RS-like model, one needs a scalar field GW1 . We, in this work, therefore look for the possibility of having a geometric modulus stabilization mechanism due to the scalar degree of freedom originating from the higher curvature terms. The important questions in this context are :

• •

Is the RS braneworld modified by $F(R)$ gravity, stabilized even without introducing an external stabilizing field ?

• •

If the modulus can be stabilized in the dual scalar-tensor model, does it mean that it is also stabilized in the original $F(R)$ model ?

• •

Does the scalar kinetic and potential terms for the purpose of modulus stabilization correspond to a $F(R)$ model which is free of ghost-like instability ?

• •

If the braneworld scenario is stabilized consistently then how the radion mass and coupling parameters will change from that of RS scenario due to the presence of higher order curvature terms ?

We aim to address these questions in this work by using the equivalence between $F(R)$ and scalar-tensor theory.

The paper is organized as follows: Following two sections are devoted to brief reviews of RS scenario and conformal relationship between $F(R)$ and scalar-tensor (ST) theory. In section IV, we extend our analysis of section III for the specific $F(R)$ model considered in this work. Section V extensively describes the modulus stabilization, radion mass and coupling for the dual ST model while section VI addresses these for the original $F(R)$ model. After discussing the equivalence, the paper ends with some conclusive remarks in section VII.

II Brief description of RS scenario and its stabilization via GW mechanism

RS scenario is defined on a five dimensional AdS spacetime involving one warped and compact extra spacelike dimension. Two 3-branes known as TeV/visible and Planck/hidden brane are embedded in five dimensional spacetime where the intermediate region between the branes is termed as ’bulk’. If $\phi$ is the extra dimensional angular coordinate, then the branes are located at two fixed points $\phi=(0,\pi)$ while the latter one is identified with our known four dimensional universe. The opposite brane tensions along with the finely tuned five dimensional cosmological constant serve as energy-momentum tensor of RS scenario. The resulting spacetime metric RS1 is non-factorizable and expressed as,

[TABLE]

Here, $r_{c}$ is the compactification radius of the extra dimension. Due to $S^{1}/Z_{2}$ compactification along the extra dimension, $\phi$ ranges from $-\pi$ to $+\pi$. The quantity $k=\sqrt{\frac{-\Lambda}{12M^{3}}}$, is of the order of 5-dimensional Planck scale $M$. Thus $k$ relates the 5D Planck scale $M$ to the 5D cosmological constant $\Lambda$.

All the dimensionful parameters described above are related to the reduced 4-dimensional Planck scale ${M}_{Pl}$ as,

[TABLE]

In order to solve the hierarchy problem, it is assumed in RS scenario that the branes are separated by such a distance that $k\pi r_{c}\approx 36$. Then the exponential factor present in the metric, which is often called warp factor, produces a large suppression so that a mass scale of the order of Planck scale is reduced to TeV scale on the visible brane. A scalar mass say mass of Higgs is given as,

[TABLE]

where $m_{H}$ and $m_{0}$ are physical and bare Higgs mass respectively. But one of the crucial aspects of this braneworld scenario is to stabilize the distance between the branes (known as modulus or radion). For this purpose, Goldberger and Wise demonstrated that the modulus corresponding to the radius of the extra dimension in RS warped geometry model can be stabilized GW1 by invoking a massive scalar field in the bulk with appropriate vacuum expectation values (vev) at the two 3-branes that reside at the orbifold fixed points. Consequently the phenomenology of the radion field originating from 5D gravitational degrees of freedom has also been explored GW_radion .

III Transformation of a F(R) theory to scalar-tensor theory

In this section, we briefly describe how a higher curvature F(R) gravity model in five dimensional scenario can be recast into Einstein gravity with a scalar field. The F(R) action is expressed as,

[TABLE]

where $x^{\mu}=(x^{0},x^{1},x^{2},x^{3})$ are usual four dimensional coordinate and $\phi$ is the extra dimensional spatial angular coordinate. $R$ is the five dimensional Ricci curvature and $G$ is the determinant of the metric. Moreover $\frac{1}{2\kappa^{2}}$ as taken as $2M^{3}$ where $M$ is the five dimensional Planck scale. Introducing an auxiliary field $A(x,\phi)$, above action (4) can be equivalently written as,

[TABLE]

By the variation of the auxiliary field $A(x,\phi)$, one easily obtains $A=R$. Plugging back this solution $A=R$ into action (5), initial action (4) can be reproduced. At this stage, perform a conformal transformation of the metric as

[TABLE]

$M,N$ run form 0 to 5. $\sigma(x,\phi)$ is conformal factor and related to the auxiliary field as $\sigma=(2/3)\ln F^{\prime}(A)$. If $R$ and $\tilde{R}$ are the Ricci scalar formed by $G_{MN}$ and $\tilde{G}_{MN}$ respectively, then they are related as,

[TABLE]

Due to the above relation between $R$ and $\tilde{R}$, action (5) turns out to be,

[TABLE]

Considering $F^{\prime}(R)>0$ and using the aforementioned relation between $\sigma$ and $F^{\prime}(A)$, one lands up to the following scalar-tensor action

[TABLE]

where $\tilde{R}$ is the Ricci scalar formed by $\tilde{G}_{MN}$. $\sigma(x,\phi)$ is the scalar field, emerging from higher curvature degrees of freedom. At this stage, it is important to note that for $F^{\prime}(R)<0$, the kinetic term of the scalar field $\sigma$ as well as the Ricci scalar $\tilde{R}$ in the above action come with wrong sign, which indicates the existence of ghost field. Thus to avoid the ghost like structure, $F^{\prime}(R)$ must be greater than zero. Later we shall show that in the context of the present work, this condition is indeed satisfied.

The kinetic part of $\sigma(x,\phi)$ in eqn.(7), though correct in sign, is non-canonical. In order to make the scalar field canonical, transform $\sigma$ $\rightarrow$ $\Phi(x,\phi)=\sqrt{3}\frac{\sigma(x,\phi)}{\kappa}$. In terms of $\Phi(x,\phi)$, the above action takes the form,

[TABLE]

where $V(\Phi)=\frac{1}{2\kappa^{2}}[\frac{F(A)}{F^{\prime}(A)^{5/3}}-\frac{A}{F^{\prime}(A)^{2/3}}]$ is the scalar field potential which depends on the form of $F(R)$. Thus the action of $F(R)$ gravity in five dimension can be transformed into the action of a scalar-tensor theory by a conformal transformation of the metric.

IV RS like spacetime in F(R) model and corresponding scalar-tensor theory

In the present work, we consider a five dimensional AdS spacetime with two 3-brane scenario in F(R) model. To the leading order in higher curvature term, the form of $F(R)$ is taken as $F(R)=R+\alpha R^{2}$ where $\alpha$ is a constant with square of the inverse mass dimension. Considering $\phi$ as the extra dimensional angular coordinate, two branes are located at $\phi=0$ (hidden brane) and at $\phi=\pi$ (visible brane) respectively while the latter one is identified with the visible universe. Moreover the extra dimension is $S^{1}/Z_{2}$ orbifolded along the coordinate $\phi$. The action for this model is :

[TABLE]

where $\Lambda(<0)$ is the bulk cosmological constant and $V_{h}$, $V_{v}$ are the brane tensions on hidden, visible brane respectively.

This higher curvature like $F(R)$ model (in eqn.(8)) can be transformed into scalar-tensor theory by using the technique discussed in the previous section. Performing a conformal transformation of the metric as

[TABLE]

the above action (in eqn.(8)) can be expressed as a scalar-tensor theory with the action given by :

[TABLE]

where the quantities in tilde are reserved for ST theory. $\tilde{R}$ is the Ricci curvature formed by the transformed metric $\tilde{G}_{MN}$. $\Phi(x,\phi)$ is the scalar field corresponds to higher curvature degrees of freedom and $V(\Phi)$ is the scalar potential which for this specific choice form of $F(R)$ has the form,

[TABLE]

One can check that the above potential (in eqn.(11)) is stable for the parametric regime $\alpha<0$. This immediately ensures that for our model with negative bulk curvature, $F^{\prime}(R)=1+\alpha R$ is always greater than zero which in turn indicates that the original $F(R)$ model is free from Ostragradsky instability as discussed earlier.

The stable value ($<\Phi>$) as well as the mass squared ($m_{\Phi}^{2}$) of the scalar field ($\Phi$) are given by the following two equations

[TABLE]

and

[TABLE]

Furthermore, the minimum value of the potential i.e. $V(<\Phi>)$ is non zero and serves as a cosmological constant. Thus the effective cosmological constant in scalar-tensor theory is $\Lambda_{eff}=\Lambda+V(<\Phi>)$ where $V(<\Phi>)$ is,

[TABLE]

Above form of $V(<\Phi>)$ with $\Lambda<0$ clearly indicates that $\Lambda_{eff}$ is also negative or more explicitly,the corresponding scalar-tensor theory for the original $F(R)$ model has an AdS like spacetime. Considering $\xi$ as the fluctuation of the scalar field over its vev, the final form of action for the scalar-tensor theory in the bulk can be written as,

[TABLE]

where the terms up to quadratic order in $\xi$ are retained for $\kappa\xi<1$. A detailed justification of neglecting the higher order terms as well as their possible effects will be discussed in section VI-C.

V Modulus stabilization, radion mass and coupling in scalar-tensor (ST) theory

V.1 Modulus stabilization

In order to stabilize the modulus in ST theory, here we adopt the GW mechanism GW1 which requires a massive scalar field in the bulk. For the case of ST theory presented in eqn.(15), $\xi$ can act as a bulk scalar field with the mass given by eqn.(13). Considering a negligible backreaction of the scalar field ($\xi$) on the background spacetime, the solution of metric $\tilde{G}_{MN}$ is exactly same as RS model i.e.

[TABLE]

where $k=\sqrt{\frac{-\Lambda_{eff}}{24M^{3}}}$. With this metric, the scalar field equation of motion in the bulk is following,

[TABLE]

where the scalar field $\xi$ is taken as function of extra dimensional coordinate only. Considering non zero value of $\xi$ on branes, the above equation (17) has the general solution,

[TABLE]

with $\nu=\sqrt{4+m_{\Phi}^{2}/k^{2}}$. Moreover $A$ and $B$ are obtained from the boundary conditions, $\xi(0)=v_{h}$ and $\xi(\pi)=v_{v}$ as follows :

[TABLE]

and

[TABLE]

Upon substitution the form of $A$ and $B$ into eqn.(18), one finds that

[TABLE]

Above values of $\xi(0)$ and $\xi(\pi)$ matches with the boundary condition (i.e. $\xi(0)=v_{h}$ and $\xi(\pi)=v_{v}$) by neglecting the subleading powers of $e^{-kr_{c}\pi}$, as have been done earlier by the authors in GW1 .

Plugging back the solution of $\xi(\phi)$ (eqn.(18)) into scalar field action and integrating over $\phi$ yields an effective modulus potential having the following form,

[TABLE]

This potential has a minimum at

[TABLE]

According to GW mechanism GW1 , the branes are stabilized at that separation for which the modulus potential becomes minimum. Thus the above equation represents the stabilized value for inter-brane separation. Expression of $m_{\Phi}^{2}$ (eqn.(13)) clearly indicates that $r_{c}$ in eqn.(22) is positive only for $\alpha<0$. Thus the $F(R)$ model depicted in eqn.(8) with positive $\alpha$ can be transformed to a scalar-tensor theory where the scalar field has negative squared mass and the modulus of the ST theory can not be stabilized.

It may be observed that the scalar field degrees of freedom is related to the curvature as,

[TABLE]

Recall that $<\Phi>=\frac{2}{\sqrt{3}\kappa}\ln[\sqrt{9-40\kappa^{2}\alpha\Lambda}-2]$.

From the above expression, we can relate the boundary values of the scalar field (i.e $\xi(0)=v_{h}$ and $\xi(\pi)=v_{v}$) with the Ricci scalar as,

[TABLE]

and

[TABLE]

where $R(0)$ and $R(\pi)$ are the values of the curvature on Planck and TeV brane respectively. In section VI-A, we derive the expression of the bulk scalar curvature which in this scenario becomes dependent on the bulk coordinate $y$. Thus the parameters that are used in the scalar-tensor theory are actually related to the parameters of the original $F(R)$ theory.

Furthermore to derive the stabilization condition in scalar-tensor theory, the backreaction of the scalar field on spacetime geometry is neglected. It can be shown from GW1 , that this is valid as long as the stress energy tensor for the scalar field is less than the bulk cosmological constant which in turn implies that $v_{h}^{2}/M^{3}$ and $v_{v}^{2}/M^{3}$ are less than unity, where $v_{h}$ and $v_{v}$ are the boundary values of the scalar field. Now using eqn. (24) and eqn. (25), we can determine the conditions of negligible back-reactions in terms of the parameters appearing in the original $F(R)$ theory.

The effect of backreaction, though small, shall also incorporated in section VI-C. We will show that the backreaction modifies all the quantities described via eqn.(18, 19, 20, 22), though the modification is small in the limit $\kappa v_{h}<1$.

V.2 Radion Potential

In this section, we consider a fluctuation of branes around the stable separation ($r_{c}$). So the inter-brane separation can be considered as a field, and here, for simplicity we assume GW_radion that this new field depends only on the brane coordinates. The corresponding metric ansatz is,

[TABLE]

From the perspective of four dimensional effective theory, $\tilde{T}(x)$ is known as radion field. Recall that the quantities in tilde are reserved for ST theory. In order to find the radion mass, here we adopt the method proposed by Goldberger and Wise GW_radion . In the GW mechanism, the same bulk scalar field which stabilize the modulus can also generate the potential for radion field and in the present ST theory (eqn.15), $\xi(x,\phi)$ can fulfill the purpose.

With the metric in eqn.(26), a Kaluza-Klein reduction for the five dimensional Einstein-Hilbert action reduces to four dimensional effective action as,

[TABLE]

As we see that $T(x)$ is not canonical and thus we redefine the field by the following transformation,

[TABLE]

In terms of $\tilde{\Psi}$, the kinetic part of radion field becomes

[TABLE]

Correspondingly the radion potential is obtained from eqn.(21) by replacing $r_{c}$ by $\tilde{T}(x)$ i.e.

[TABLE]

where $A$ and $B$ are given by,

[TABLE]

Using the above potential in eqn.(28) with the transformation given in eqn.(27), one lands up with the following mass squared of radion field GW_radion in scalar-tensor theory as

[TABLE]

with $\epsilon=m_{\Phi}^{2}/4k^{2}$.

V.3 Coupling between radion and Standard Model fields

Being a gravitational degree of freedom, radion field interacts with brane energy-momentum tensor and the couplings of interaction are constrained by four dimensional general covariance. Following GW_radion , it can be shown that the coupling between radion and SM fields (Higgs for example) is given by

[TABLE]

where $\mu$ is mass of Higgs particle and $k\pi r_{c}$ is obtained from eqn.(22). Similarly the coupling between radion and $Z$ boson is,

[TABLE]

where $m_{Z}$ is the mass of $Z$ boson.

Now we turn our focus on modulus stabilization as well as on radion mass and coupling for the original $F(R)$ model (eqn.(8)) by using the stabilization condition of the corresponding scalar-tensor theory.

VI Modulus stabilization, radion mass and coupling in F(R) model

VI.1 Modulus stabilization

Recall that the original higher curvature $F(R)$ model is presented by action given in eqn.(8). Solutions of metric ($G_{MN}$) for this $F(R)$ model can be extracted from the solutions of corresponding scalar-tensor theory (eqn.(16) and eqn.(18)) with the help of eqn.(9). Thus the line element in $F(R)$ model turns out to be

[TABLE]

where $\Phi(\phi)=<\Phi>+\xi(\phi)$ and $\xi(\phi)$ is given by eqn.(18). This solution of $G_{MN}$ immediately leads to the separation between hidden ($\phi=0$) and visible ($\phi=\pi$) branes along the path of constant $x^{\mu}$ as follows :

[TABLE]

where $d$ is the inter-brane separation in $F(R)$ model. Using the explicit functional form of $\Phi(\phi)$ (eqn.(18)), above equation can be integrated and simplified to the following one,

[TABLE]

where the sub-leading terms of $\kappa\xi$ are neglected. $r_{c}$ is the modulus in the corresponding ST theory and it is stabilized which is shown in the previous section (eqn.(22)). So, it can be argued that due to the stabilization of ST theory, the modulus $d$ in $F(R)$ model is also stabilized with a value,

[TABLE]

Hence it can be concluded that $F(R)$ model where the only independent field is spacetime metric ($G_{MN}$), is a self stabilizing system. From the expression of $m_{\Phi}^{2}$ (eqn.(13)), it is clear that $d$ goes to zero at the limit $\alpha\rightarrow 0$. Moreover for $\alpha>0$, $m_{\Phi}^{2}$ becomes negative which in turn makes the modulus $d$ negative (see eqn.(35)), an unphysical situation. From the above two statements, it is clear that the self stabilization in $F(R)$ model arises entirely due to the presence of higher curvature term $\alpha R^{2}$ only when $\alpha<0$.

From the above relation (eqn.(35)) and using the expression of $m_{\Phi}^{2}$ (eqn.(13)), we obtain the figure 1 between $k\pi d$ and $|\alpha|M^{2}$.

The figure demonstrates that the brane separation ($d$) increases with the higher curvature parameter $\alpha$.

However, using the metric solution given in eqn.(33), one calculates the five dimensional Ricci scalar as follows:

[TABLE]

where $c$ is an integration constant. Recall that the boundary values of the curvature (i.e. $R(0)$ and $R(\pi)$, obtained from the eqn.(36)) are related with $v_{h}$ and $v_{v}$ by eqn.(24) and eqn.(25) respectively. Thus one can tune $R(0)$, $R(\pi)$ to fix $v_{h}$ and $v_{v}$ in a desirable value.

VI.2 Radion Potential

A fluctuation of branes around the stable configuration $d$ is now considered. This fluctuation can be taken as a field ($T(x)$) and for simplicity, this new field is assumed to be the function of brane coordinates only. The metric takes the following form,

[TABLE]

From the angle of four dimensional effective theory, $T(x)$ is known as radion field. Moreover $\Phi(x,\phi)$ is obtained from eqn.(18) by replacing $r_{c}$ to $T(x)$. Using the metric in eqn.(37), a Kaluza-Klein reduction for five dimensional Einstein-Hilbert action generates a kinetic as well as a potential part for the radion field $T(x)$. Kinetic part comes as

[TABLE]

It is evident that $T(x)$ is not canonical and in order to make it canonical, we redefine the field as

[TABLE]

where $f=\sqrt{\frac{24M^{3}}{k}}[1-\frac{20}{\sqrt{3}}\alpha k^{2}\kappa v_{h}]$. For $\alpha\rightarrow 0$, the action contains only the linear term in Ricci scalar and the factor $f$ goes to $\sqrt{24M^{3}/k}$ which agrees with GW_radion . With the redefined field, kinetic part of radion field becomes,

[TABLE]

Finally the Kaluza-Klein reduction leads to the potential part of radion field as follows,

[TABLE]

It may be observed that $V(\Psi)$ goes to zero as $\alpha$ tends to zero. This is expected because for $\alpha\rightarrow 0$, the action contains only the Einstein part which does not produce any potential term for the radion field GW_radion . Thus for five dimensional warped geometric model, the radion potential is generated from the higher order curvature term $\alpha R^{2}$. The potential in eqn.(40) has a vev at

[TABLE]

which leads to the inter-brane separation as,

[TABLE]

in the leading order of $\kappa\xi$. This above equation exactly resembles with the eqn.(35) which again indicates that the five dimensional warped spacetime we considered, is self stabilized by higher curvature gravitational theory. We emphasize that due to the presence of conformal factor connecting the two theories, the value of $k\pi d$ (in $F(R)$ model) is less than $k\pi r_{c}$ (in ST model). However we find that the stabilization of modulus remain intact in both the models. Finally the squared mass of radion field is as follows,

[TABLE]

It is noticed that the mass of radion field is enhanced by the higher curvature terms in five dimensional gravitational action.

VI.3 Radion potential with backreaction

It may be mentioned that the stabilized interbrane separation obtained in eqn. (42) is based on the conditions that the bulk scalar potential is retained up to quadratic term (see eqn.(15)) and the backreaction of the bulk scalar field is neglected on five dimensional spacetime. Both these conditions are followed from the assumption that $\kappa v_{h}<1$. Relaxation of this assumption is crucial to check the status of the stabilization condition in the presence of higher order self interaction terms in the bulk scalar potential. Here we examine whether the radion potential admits any stability when $V(\Phi)$ is retained up to cubic term in $\Phi$. In this scenario, the five dimensional action in ST theory turns out to be,

[TABLE]

where $g$ is the self cubic coupling of $\Phi(\phi)$ and can be easily determined from the form of $V(\Phi)$ presented in eqn.(11) as,

[TABLE]

Considering the metric ansatz as,

[TABLE]

the gravitational as well as the scalar field equations of motion take the following form,

[TABLE]

[TABLE]

[TABLE]

To determine the solutions of the above differential equations, we apply the iterative method by considering the form of metric determined in eqn. (16) as the zeroth order solution. In the leading order correction of $\kappa v_{h}$, $\xi(\phi)$ and $A(\phi)$ turn out to be

[TABLE]

and

[TABLE]

Thus due to the inclusion of the bulk scalar field backreaction, the warp factor gets modified and the correction term is proportional to $\kappa^{2}v_{h}^{2}$ which is indeed small for $\kappa v_{h}<1$. However if one includes these corrections in $\xi(\phi)$ and $A(\phi)$, one can extract the modified solution in F(R) model by a conformal transformation as indicated earlier:

[TABLE]

where $ds^{2}$ is the line element in $F(R)$ model and $\xi(\phi)$, $A(\phi)$ are given in eqn.(50) and eqn.(51) respectively. Plugging back the above solution of metric (in eqn.(52)) into the five dimensional $F(R)$ action and integrating over $\phi$, the radion potential is obtained as

[TABLE]

The above potential has a stable minimum at,

[TABLE]

which leads to the stabilized interbrane separation as:

[TABLE]

Comparing eqn.(42) and eqn.(55), it can be seen that the vev of the radion field and hence the stable modulus is shifted due to the presence of higher order self interaction terms in the bulk scalar potential or the inclusion of scalar field backreaction. However this shift is indeed small in the limit $\kappa v_{h}<1$

VI.4 Coupling between radion and Standard Model fields

The radion field arises as a scalar degree of freedom on the TeV brane and has interactions with the Standard Model (SM) fields. From the five dimensional metric (eqn.(37)), it is clear that the induced metric on visible brane is $(\frac{\Psi}{f})^{2}e^{-\kappa v_{v}\frac{1}{\sqrt{3}}}\eta_{\mu\nu}$ (where $f=\sqrt{\frac{24M^{3}}{k}}[1-\frac{20}{\sqrt{3}}\alpha k^{2}\kappa v_{h}]$) and consequently $\Psi(x)$ couples directly with SM fields.

For example, consider the Higgs sector of Standard Model,

[TABLE]

where $h(x)$ is the Higgs field. In order to get a canonical kinetic term, one needs to redefine $h(x)\longrightarrow H(x)=\frac{<\Psi>}{f}h(x)$. Therefore for $H(x)$, the above action can be written as,

[TABLE]

where $\mu=\mu_{0}\frac{<\Psi>}{f}e^{-\kappa v_{v}\frac{1}{\sqrt{3}}}$. Considering a fluctuation of $\Psi(x)$ about its vev as $\Psi(x)=<\Psi>+\delta\Psi$, one can obtain (from eqn.(57)) that $\delta\Psi$ couples to $H(x)$ through the trace of the energy-momentum tensor of the Higgs field:

[TABLE]

So, the coupling between radion and Higgs field become, $\lambda_{(H-\delta\Psi)}=\frac{\mu^{2}}{<\Psi>}$. Similar consideration holds for any other SM fields. For example for $Z$ boson, the corresponding coupling is $\lambda_{(Z-\delta\Psi)}=\frac{m^{2}_{Z}}{<\Psi>}$. Thus the inverse of $<\Psi>$ plays a crucial role in determining the coupling strength between radion and SM fields. In the present case, we obtain

[TABLE]

Hence finally we arrive at,

[TABLE]

and

[TABLE]

in the leading order of $\kappa$.

The coupling between radion and fermion field is similarly given by,

[TABLE]

If the fermion fields are allowed to propagate in the bulk, then the coupling with radion changes from that given in (eqn.(60)). This can be determined by Kaluza-Klein (KK) mode expansion of the bulk fermion grossman in a bulk governed by $F(R)$ gravity. Here we focus on the zeroth order KK mode and the solution of the zeroth order KK mode wave function in the bulk takes the form:

[TABLE]

where $\chi_{L}(\phi)$ is the left handed fermion wave function and $m$ is bulk fermionic mass. Recall that $d$ is the interbrane separation and $<\Phi>$ is given in eqn. (12). Similarly for the right handed mode ($\chi_{R}(\phi)$),

[TABLE]

Using these above solutions, we determine the coupling of radion with zeroth order fermionic KK mode, which yields,

[TABLE]

for left handed chiral mode and,

[TABLE]

for right handed mode.

The form of $v_{h}$ and $v_{v}$ in terms of five dimensional Ricci scalar can be extracted from eqn. (24) and eqn. (25). It is evident that the coupling $\lambda_{L}$, $\lambda_{R}$ is modified by the factor $(k\pm 2m)d$, in comparison to the coupling given in eqn. (60) for the fermion fields confined on the TeV brane.

From above analysis we note that the coupling between radion and SM fields is suppressed due to the presence of higher curvature parameter $\alpha$ which in turn modifies the phenomenology on visible 3-brane.

Before concluding, we mention about a recent work ssg3 , where a higher curvature gravity model with $R^{4}$ terms in the action is considered. The corresponding conformally transformed scalar action includes a quartic term which resembles closely to the scalar action considered in reference [11 and 13]. It has been shown that with such a specific choice of the scalar action one can estimate the exact modification of the warp factor due to the effects of the back-reaction of the scalar field on the background geometry and thus it enables us to address the role of back-reaction on the stabilization issue and also on various parameters of the low-energy effective action. However in such a model, the back-reacted geometry can be exactly determined if the scalar mass and the quartic term in the potential are inter-related and in the limit of the quartic term going to zero, the mass term also goes to zero. Therefore there is no smooth limit which takes this model to that considered by the GW where only quadratic mass term was present. The work reported in this article however has a different goal from that of ssg3 . Here we show that in the leading order quadratic curvature correction to Einstein action in the bulk we find a dual scalar tensor theory which under certain approximation is similar to the original GW scalar action which has a quadratic mass term only. We therefore explore and re-examine the originally proposed Goldberger-Wise modulus stabilization condition in the light of higher curvature gravity models where such a stabilizing field appears naturally from higher curvature degrees of freedom with a minimal curvature extension.

VII Conclusion

In this work, we consider a five dimensional AdS, compactified warped geometry model with two 3-branes embedded within the spacetime. Due to large curvature ($\sim$ Planck scale) in the bulk, the spacetime is assumed to be governed by a higher curvature gravity model , $F(R)=R+\alpha R^{2}$. In this scenario, we address how the modulus stabilization and radion phenomenology are affected by higher curvature term. The findings and implications of our results are as follows :

• •

The model comes as a self stabilizing system due to the presence of higher curvature term $\alpha R^{2}$. This is in sharp contrast to a model with only Einstein term in the bulk where the modulus can not be stabilized without incorporating any external degrees of freedom such as a scalar field. However for the higher curvature gravity model, this additional degree of freedom originates naturally from the higher curvature term and plays the role of a stabilizing field. It may also be noted that for $\alpha\rightarrow 0$, the stabilization condition (eqn.(35)) leads to zero brane separation.

• •

We scan the parametric space of $\alpha$ for which the modulus is going to be stabilized. Our result reveals that for $\alpha>0$, the inter-brane separation becomes negative which is an unphysical situation. Thus the braneworld we have considered is stabilized only when $\alpha<0$. This puts constraints on the $F(R)$ model itself. Moreover the distance between the branes increases with the value of the parameter $\alpha$ which is evident from figure (1). Thus the results obtained in this work clearly bring out the correlation between a geometrically stable warped solution resulting from negative bulk curvature and the stability of the higher curvature $F(R)$ model free from ghosts.

• •

Quadratic term in curvature also generates the radion potential with a stable minimum. We find the radion mass as well as radion coupling with SM fields. The expressions of mass (eqn.(43)) and coupling (eqn.(58) and eqn.(59)) clearly indicate that the radion mass is enhanced while the coupling is suppressed in comparison to the scenario where only Einstein gravity resides in the bulk GW_radion . Thus the cross section between radion and SM fields is overall suppressed due to the presence of higher order curvature terms in five dimensional gravitational action leading to a possible explanation of the invisibility of the radion field in the present experimental resolution.

• •

The effect of backreaction on the radion potential and its minima are studied. It is shown that the corrections due to the backreaction is indeed small in the limit $\kappa v_{h}<1$. The possible correction terms for the backreaction are determined.