text:
# 1 Introduction

## 1 Introduction

Extra dimensional theories where standard model fields are localized on a brane  have recently attracted a lot of attention, since such models have several distinct features from ordinary Kaluza-Klein (KK) theories. In particular, Randall and Sundrum  presented a simple model based on two branes and a single extra dimension, where the hierarchy problem could be solved due to the exponentially changing metric along the extra dimension. In order to obtain a phenomenologically acceptable model, the radion field (which corresponds to fluctuations in the distance of the two branes) has to get a mass, otherwise it would violate the equivalence principle , and also result in unconventional cosmological expansion equations . The simplest mechanism for radius stabilization has been suggested by Goldberger and Wise , who employed an additional bulk scalar which has a bulk mass term and also couples to both branes (for another issues related to the radion potential see ). Both the cosmology and collider phenomenology crucially depend on the mass and couplings of the radion. In particular, there is no radion moduli problem if the radion mass is $`O(\mathrm{TeV})`$ and its couplings to Standard Model (SM) fields is $`O(\mathrm{TeV}^1)`$. This is also the most favorable scenario for discovering the radion at a future collider.

In fact it was shown in , that the radion will have the above properties for the Goldberger–Wise scenario. However, the calculation of  were using a naïve ansatz for the radion field which ignores both the radion wavefunction and the backreaction of the stabilizing scalar field on the metric. The validity of this approximation has recently been questioned .

Therefore in this paper we analyze the coupled radion-scalar system in detail from the 5D point of view. We derive the coupled differential equations governing the dynamics of the system, and find the mass eigenvalues for some limiting cases. Due to the coupling between the radion and the bulk scalar, we find that there will be a single KK tower describing the system, with the metric perturbations non-vanishing for every KK mode. This implies that the Standard Model fields localized on the TeV brane will couple to every KK mode from the bulk scalar, and this could provide a means to directly probe the stabilizing physics.

Using the coupled equations for the radion–scalar system, we analyze the late–time behavior of the radion in an expanding universe, and find that the troubling 55–component of Einstein’s equation just determines the shift of the radion. This shift completely agrees with the shift obtained in  using the 4D effective theory.

Given that we have established that the radion mass is $`O(\mathrm{TeV})`$ and that its couplings to SM particles is $`O(\mathrm{TeV}^1)`$, it is reasonable to consider its effects on SM phenomenology. Some direct collider signatures for the radion and loop corrections have been discussed in . In the second half of this paper the effects of the radion on the oblique parameters***Loop effects for theories with large extra dimensions have been analyzed in . are calculated using an effective theory approach similar to . Since in the RS model the radion is the only new state well below the TeV scale, a low–energy effective theory including only the radion and SM fields is used. The effects of other heavy modes are accounted for by including non–renormalizable operators at the cutoff scale. In the absence of a curvature–scalar Higgs mixing operator, the corrections from the radion are small, but give a negative contribution to $`S`$. In the presence of such a mixing operator the corrections could be sizable due to the modified radion and Higgs couplings.

This paper is organized as follows: in Section 2 we review the Randall-Sundrum model and radius stabilization by bulk scalar fields. We also summarize the explicit example of  which we will be using for our explicit computation of the radion mass and couplings to SM fields. In Section 3 we present our ansatz for the coupled metric and scalar fluctuations based on the analysis of  of the radion without a stabilizing potential. We will derive a single ordinary differential equation, whose eigenmodes will yield the KK modes for the radion-scalar system. In Section 4 we analyze the generic properties of this equation. In the general case we find that the system is not described by a hermitian Schrödinger operator. However, we identify a convenient limit, in which the differential operator is in fact hermitian, and the eigenfunctions are manifestly orthogonal. In Section 5 we analyze the eigenfunctions in this limit, and find the approximate masses for the KK tower. In this analysis, the back reaction of the metric is neglected, which results in the lightest mode still being massless. The effect of the back reaction on the lightest mode is taken into account in Section 6, where we find the mass of the radion to be of the order (but slightly lighter) than the weak scale. In Section 7 we discuss the couplings of the radion and the KK tower to SM fields on the brane. We find that the radion coupling is exactly agrees with the results in , while the couplings of the other KK modes of the scalar field are suppressed by the mass of the given mode, and is proportional to the backreaction of the metric due to the scalar background. In Section 8 we demonstrate that in an expanding universe the shift in the radion at late times agrees with the 4D effective theory result obtained in . Having established the mass and coupling of the radion, we write an effective Lagrangian in Section 9 without any specific mechanism of radius stabilization and neglecting the contributions of the KK modes. In Section 10 we add a curvature-scalar Higgs coupling to the effective Lagrangian, and discuss how the couplings are modified. Then, in Section 11 we calculate the Feynman rules in a general gauge. These allow us to compute the oblique parameters via one-loop vacuum polarization diagrams with radions in the loop in Section 12. The radion correction is log divergent (unlike the Higgs), and so we also write the nonrenormalizable operators at the cutoff scale that provide the necessary counterterms. The size of the new contributions are shown for various cases in several figures in Section 12.1. We also estimate limits on the radion mass as a function of the cutoff scale in Section 13. Finally, we conclude in Section 14.

## 2 Review of the Randall-Sundrum Model and the
Goldberger-Wise Mechanism

Randall and Sundrum presented a very interesting proposal for solving the hierarchy problem . By introducing a fifth dimension where the bulk geometry is anti-de Sitter, a large hierarchy between the Planck scale and the TeV scale is obtained with only a mild finetuning. Two branes are introduced, located at the boundaries of the anti-de Sitter space. By tuning the bulk cosmological constant $`\mathrm{\Lambda }6k^2/\kappa ^2`$, the tensions $`V_P`$ and $`V_T`$ on the Planck and TeV branes, respectively, such that $`V_P=V_T=6k/\kappa ^2`$ (where $`\kappa ^2`$ is the 5D Newton constant related to the 5D Planck mass by $`\kappa ^2=1/2M^3`$) one obtains a 4-D Poincare invariant solution. The metric is then

$$ds^2=e^{2ky}\eta _{\mu \nu }dx^\mu dx^\nu dy^2,$$
(2.1)
where the Planck brane and TeV branes are located at $`y=0`$ and $`y=r_0`$. For a moderate choice of $`kr_0O(50)`$, a large hierarchy between the Planck scale and the weak scale is generated.

Since this solution is obtained for any value of $`r_0`$, some mechanism is required to fix $`r_050/k`$ as opposed to some other value of $`r_0`$. This must also be done without introducing any large finetuning. Further, small shifts in the separation between the two branes do not change the energy, and so are described in an effective theory by the fluctuations of a massless particle, the “radion”. This particle couples like a Brans-Dicke scalar and must be massive to recover ordinary 4-D Einstein gravity

One way to achieve these requirements is to introduce a bulk scalar field $`\varphi `$ that has a bulk potential $`V(\varphi )`$ . To stabilize the brane distance, potentials $`\lambda _{P,T}(\varphi )`$ on the Planck and TeV branes respectively are also included. The competition between the brane and bulk Lagrangians generates a vacuum expectation value (vev) for $`\varphi `$, which results in a 4-D vacuum energy that depends on $`r_0`$. For a simple choice of polynomial potentials a large hierarchy is then easily obtained with a mild finetuning , and the resulting mass for the radion is $`O(\mathrm{TeV})`$ .

The phenomenology of the radion depends on the strength of its coupling to the brane fields. Using the following naïve ansatz to describe the radion $`b(x)`$

$$ds^2=e^{2k|y|b(x)}ds_4^2b(x)^2dy^2,$$
(2.2)
 and  computed the normalization of the radion kinetic term to be

$$\frac{3}{4}e^{2kr_0}\frac{kr_0^2}{\kappa ^2}(b)^2.$$
(2.3)
Fields living on the TeV brane couple to the radion through the induced metric, with an interaction

$$\frac{kr_0}{2}b(x)\text{Tr}T_{\mu \nu }=\frac{r(x)}{\sqrt{6}\mathrm{\Lambda }_W}\text{Tr}T_{\mu \nu }$$
(2.4)
where $`r`$ is the canonically normalized radion, $`\mathrm{\Lambda }_W=M_{Pl}e^{kr_0}O(\mathrm{TeV})`$, $`M_{Pl}^2=(1e^{2kr_0})/(k\kappa ^2)1/(k\kappa ^2)`$, and $`T_{\mu \nu }`$ is the physical energy-momentum tensor of the TeV brane fields. It is then clear that the radion couples as $`1/`$TeV to the Standard Model fields. Obtaining an acceptable phenomenology then requires that the radion mass is $`O(`$TeV$`)`$, which is easily satisfied by the Goldberger–Wise mechanism.

The phenomenology of the radion is then crucially dependent on the normalization of the kinetic term. In fact, in the computation leading to the $`O(`$TeV$`)^2`$ prefactor in (2.3) there is a cancellation between two terms of $`O(M_{Pl}^2)`$. The origin of this cancellation remains somewhat mysterious, and the absence of this cancellation would clearly lead to different predictions. In  it was pointed out that there are additional contributions to the radion kinetic term not included in . In particular, the profile of the stabilizing field depends on $`r_0`$, and so a small change in $`r_0b(x)`$ distorts the background field. It was found that this results in an $`O(M_{Pl}^2)`$ correction to the radion kinetic term, thereby drastically changing the phenomenology of the radion.

We review the resolution of this issue in the first part of this paper. Some of the results presented in Sections 36 are already contained in the work by Tanaka and Montes , even though the results of this paper were obtained independently of Ref.. We explicitly determine the wavefunction of the radion when there is a stabilizing mechanism. We find that the radion mass is typically $`O(`$TeV$`)`$. In the limit that the backreaction of the stabilizing fields on the metric is small, we find that the correction of the stabilizing field to the radion kinetic term is subdominant to the gravitational contribution. We also find that once the stabilizing field has a non–zero vev, the Kaluza-Klein (KK) tower couples directly to the brane world fields, with $`1/`$TeV normalization, and amplitude depending on the size of the backreaction.

The action we consider is ***The action is integrated over the circle rather than the line segment.

$$M^3d^5x\sqrt{g}+d^5x\sqrt{g}\left(\frac{1}{2}\varphi \varphi V(\varphi )\right)d^4x\sqrt{g_4}\lambda _P(\varphi )d^4x\sqrt{g_4}\lambda _T(\varphi ),$$
(2.5)
where $`g_4`$ is the induced metric on the branes. The background vev for $`\varphi `$ and background metric that preserve 4-D Lorentz invariance is

$$\varphi (x,y)=\varphi _0(y),$$
(2.6)
$$ds^2=e^{2A}\eta _{\mu \nu }dx^\mu dx^\nu dy^2.$$
(2.7)
The Einstein equations are then

$$R_{ab}=\kappa ^2\stackrel{~}{T}_{ab}=\kappa ^2\left(T_{ab}\frac{1}{3}g_{ab}g^{cd}T_{cd}\right),$$
(2.8)
with $`\kappa ^2=1/(2M^3)`$. For this background the scalar and metric field equations are

$`4A^2A^{\prime \prime }`$ $`=`$ $`{\displaystyle \frac{2\kappa ^2}{3}}V(\varphi _0){\displaystyle \frac{\kappa ^2}{3}}{\displaystyle \underset{i}{}}\lambda _i(\varphi _0)\delta (yy_i),`$ (2.9)
$`A^2`$ $`=`$ $`{\displaystyle \frac{\kappa ^2\varphi _{0}^{}{}_{}{}^{2}}{12}}{\displaystyle \frac{\kappa ^2}{6}}V(\varphi _0)`$ (2.10)
$`\varphi _0^{\prime \prime }`$ $`=`$ $`4A^{}\varphi _0^{}+{\displaystyle \frac{V(\varphi _0)}{\varphi }}+{\displaystyle \underset{i}{}}{\displaystyle \frac{\lambda _i(\varphi _0)}{\varphi }}\delta (yy_i).`$ (2.11)
Here primes denote $`/y`$, and we reserve $`_\mu `$ to denote derivative with respect to the comoving 4-D spacetime coordinates $`x^\mu `$. The boundary equations for $`A`$ and $`\varphi _0`$ are obtained by matching the singular terms in the above equations. This gives

$$[A^{}]|_i=\frac{\kappa ^2}{3}\lambda _i(\varphi _0),$$
(2.12)
$$[\varphi _0^{}]|_i=\frac{\lambda _i(\varphi _0)}{\varphi }.$$
(2.13)
For analytical solutions we use an approach presented in . A particular class of potentials $`V`$ is considered which can be written in the form

$$V(\varphi )=\frac{1}{8}\left(\frac{W(\varphi )}{\varphi }\right)^2\frac{\kappa ^2}{6}W(\varphi )^2.$$
(2.14)
Then a solution to the following first order equations

$$\varphi _0^{}=\frac{1}{2}\frac{W}{\varphi }\text{ , }A^{}=\frac{\kappa ^2}{6}W(\varphi _0)$$
(2.15)
automatically solves both the Einstein and scalar field equations, once the appropriate boundary conditions are solved. The virtue of this method is that for simple choices of $`W`$ it is possible to also solve for the backreaction of $`\varphi _0`$ on the metric. This will be important for us, since we find that only after including the backreaction of the stabilizing field does one find that the radion acquires a mass.

In particular, to obtain some analytic results the following superpotential  will be used:

$$W(\varphi )=\frac{6k}{\kappa ^2}u\varphi ^2$$
(2.16)
with brane potentials

$$\lambda (\varphi )_\pm =\pm W(\varphi _\pm )\pm W^{}(\varphi _\pm )(\varphi \varphi _\pm )+\gamma _\pm ^2(\varphi \varphi _\pm )^2.$$
(2.17)
Here $`+/`$ refer to Planck/TeV brane. The solution is

$`\varphi _0(y)`$ $`=`$ $`\varphi _Pe^{uy}\text{ ,}`$ (2.18)
$`A(y)`$ $`=`$ $`ky+{\displaystyle \frac{\kappa ^2\varphi _P^2}{12}}e^{2uy}.`$ (2.19)
The separation distance $`r_0`$ is then fixed by matching $`\varphi _0`$ at $`0`$ and $`r_0`$ to $`\varphi _P`$ and $`\varphi _T`$ which gives $`ur_0=\mathrm{ln}\varphi _P/\varphi _T`$. So the quantity

$$e^{ur_0}=\frac{\varphi _T}{\varphi _P}$$
(2.20)
is not a (hierarchically) small number, since both $`\varphi _P`$ and $`\varphi _T`$ are $`O(M_{Pl}^{3/2})`$. This combination will appear later in the expression for the radion mass. Also for future reference, since the backreaction corresponds to the second term in $`A`$, the limit of a small backreaction is $`\kappa ^2\varphi _P^2,\kappa ^2\varphi _T^21`$, and $`u>0`$, but with $`\varphi _P/\varphi _T=`$constant, so that $`u`$ is kept constant.

## 3 The Coupled Field Equations

When $`\varphi _0=0`$ there is always a static solution independent of the value of $`r_0`$.***After two finetunings which are independent of $`r_0`$. But only one finetune remains after radius stabilization .  The small fluctuations in the relative position between the two branes then describe a massless particle (“the radion”), and its wavefunction is  $`G(x,y)=2F(x,y)=2ke^{2ky}R(x)`$ and where $`\text{ }\text{ }\text{ }\text{ }\text{ }R=0`$. Since the coupling of the radion to the Standard Model fields is $`1/`$TeV , obtaining an acceptable phenomenology requires that this radion acquires a mass.

We therefore consider the spectrum of perturbations about the above background which stabilizes the inter–brane separation. A general ansatz to describe the spin–0 fluctuations is

$$\varphi (x,y)=\varphi _0(y)+\phi (x,y)$$
(3.1)
$$ds^2=e^{2A2F(x,y)}\eta _{\mu \nu }dx^\mu dx^\nu (1+G(x,y))^2dy^2.$$
(3.2)
In order to describe all gravitational excitations of the model, one would need to add also the degrees of freedom in the graviton, by replacing $`\eta _{\mu \nu }\eta _{\mu \nu }+h_{\mu \nu }^{TT}`$. One can show that the Einstein equations with this replacement will have the radion and the graviton decoupled. This metric ansatz (3.2) (together with the two equations (3.12) and $`G=2F`$ which we will shortly derive) fixes our gauge choice. One can show, that the effect of the remaining gauge transformations that preserve the form of (3.2) and (3.12) just amount to a 4D gauge transformation on the graviton field $`h_{\mu \nu }`$ and can be used to impose a convenient 4D gauge for the graviton. In the following we will only concentrate on the radion field. Using this ansatz the Einstein and scalar field equations are linearized to obtain some coupled equations for $`F`$, $`G`$ and $`\phi `$. The linearized Einstein equations are

$$\delta R_{ab}=\kappa ^2\delta \stackrel{~}{T}_{ab}.$$
(3.3)
Inspecting the $`\delta R_{\mu \nu }`$ equation one immediately concludes that $`G=2F`$. For

$$\delta R_{\mu \nu }=\mathrm{}+2_\mu _\nu F_\mu _\nu G+\mathrm{}$$
(3.4)
where the ellipses all contain terms $`\eta _{\mu \nu }`$. Since to linear order in the perturbations the sources $`\delta \stackrel{~}{T}_{\mu \nu }`$ are also all $`\eta _{\mu \nu }`$, the $`_\mu _\nu `$ term in $`\delta R_{\mu \nu }`$ term must vanish. This gives $`G=2F+c`$. However in the limit $`F0`$, or $`G0`$ we should recover the background solution, so $`c=0`$. In what follows we set $`G=2F`$. Then the coupled field equations are

$`\delta R_{\mu \nu }`$ $`=`$ $`\eta _{\mu \nu }\text{ }\text{ }\text{ }\text{ }\text{ }F+e^{2A}\eta _{\mu \nu }\left(F^{\prime \prime }+10A^{}F^{}+6A^{\prime \prime }F24A_{}^{}{}_{}{}^{2}F\right),`$ (3.5)
$`\delta R_{\mu 5}`$ $`=`$ $`3_\mu F^{}6A^{}_\mu F,`$ (3.6)
$`\delta R_{55}`$ $`=`$ $`2e^{2A}\text{ }\text{ }\text{ }\text{ }\text{ }F+4F^{\prime \prime }16A^{}F^{}.`$ (3.7)
The source terms are

$`\delta \stackrel{~}{T}_{\mu \nu }`$ $`=`$ $`{\displaystyle \frac{2}{3}}e^{2A}\eta _{\mu \nu }\left(V^{}(\varphi _0)\phi 2V(\varphi _0)F\right)`$ (3.8)
$`{\displaystyle \frac{1}{3}}e^{2A}\eta _{\mu \nu }{\displaystyle \underset{i}{}}\left({\displaystyle \frac{\lambda _i(\varphi _0)}{\varphi }}\phi 4\lambda _i(\varphi _0)F\right)\delta (yy_i),`$
$`\delta \stackrel{~}{T}_{\mu 5}`$ $`=`$ $`\varphi _0^{}_\mu \phi ,`$ (3.9)
$`\delta \stackrel{~}{T}_{55}`$ $`=`$ $`2\varphi _0^{}\phi ^{}+{\displaystyle \frac{2}{3}}V^{}(\varphi _0)\phi +{\displaystyle \frac{8}{3}}V(\varphi _0)F`$ (3.10)
$`+{\displaystyle \frac{4}{3}}{\displaystyle \underset{i}{}}\left({\displaystyle \frac{\lambda _i(\varphi _0)}{\varphi }}\phi +2\lambda _i(\varphi _0)F\right)\delta (yy_i).`$
The linearized scalar field equation is

$`e^{2A}\text{ }\text{ }\text{ }\text{ }\text{ }\phi \phi ^{\prime \prime }+4A^{}\phi ^{}+{\displaystyle \frac{^2V}{\varphi ^2}}(\varphi _0)\phi `$ $`=`$ $`{\displaystyle \underset{i}{}}\left({\displaystyle \frac{^2\lambda _i(\varphi _0)}{\varphi ^2}}\phi +2{\displaystyle \frac{\lambda _i(\varphi _0)}{\varphi }}F\right)\delta (yy_i)`$ (3.11)
$`6\varphi _0^{}F^{}4{\displaystyle \frac{V}{\varphi }}F.`$
Notice that the $`R_{\mu 5}`$ may be integrated immediately to obtain

$$\varphi _0^{}\phi =\frac{3}{\kappa ^2}(F^{}2A^{}F).$$
(3.12)
An integration constant $`k(y)`$ has been sent to zero since we require that the fluctuations $`F`$ and $`\phi `$ are also localized in $`x`$. This equation (3.12) together with the metric ansatz (3.2) fixes our gauge choice. One can show, that the effect of the remaining gauge transformations that preserve the form of (3.2) and (3.12) just amount to a 4D gauge transformation on the graviton field $`h_{\mu \nu }`$ and can be used to impose a convenient 4D gauge for the graviton.

These equations must be supplemented by the boundary conditions for $`F`$ and $`\phi `$ on the two branes. These are obtained by identifying the singular terms in above equations. A priori the Einstein equations give two boundary conditions for each wall. It is however straightforward to show that one of them is trivially satisfied once $`A`$ satisfies the jump equation (2.12). The two remaining boundary equations are

$$[F^{}]=\frac{2\kappa ^2}{3}\lambda _i(\varphi _0)F+\frac{\kappa ^2}{3}\frac{\lambda _i}{\varphi }(\varphi _0)\phi .$$
(3.13)
$$[\phi ^{}]|_i=\frac{^2\lambda _i}{\varphi ^2}(\varphi _0)\phi +2\frac{\lambda _i}{\varphi }F$$
(3.14)
Upon using the jump equations for the background the first equation is seen to be equivalent to (3.12) and so provides no new constraints. Then only the second boundary condition must be implemented. A convenient limit will at times be considered in this paper. The second boundary condition simplifies in the limit of a stiff boundary potential. Namely, if $`^2\lambda _i/\varphi ^21`$ then the second boundary condition is just $`\phi |_i=0`$. Then in this limit the first boundary condition is just

$$(F^{}2A^{}F)|_i=0.$$
(3.15)
A single equation for $`F`$ is obtained as follows. One considers the combination $`e^{2A}\delta R_{\mu \nu }+\delta R_{55}`$ in the bulk. The point of this combination is to eliminate terms of the form $`V(\varphi _0)\phi `$. This leaves a bulk equation for $`F`$ and $`\phi ^{}`$ only:

$$e^{2A}\text{ }\text{ }\text{ }\text{ }\text{ }F+F^{\prime \prime }2A^{}F^{}=\frac{2\kappa ^2}{3}\varphi _0^{}\phi ^{}$$
(3.16)
One then eliminates $`\phi ^{}`$ in favor of $`F`$ using (3.12). This gives

$$F^{\prime \prime }2A^{}F^{}4A^{\prime \prime }F2\frac{\varphi _0^{\prime \prime }}{\varphi _0^{}}F^{}+4A^{}\frac{\varphi _0^{\prime \prime }}{\varphi _0^{}}F=e^{2A}\text{ }\text{ }\text{ }\text{ }\text{ }F,$$
(3.17)
to be solved in the bulk. This is the principle equation that will be studied and solved below. We note in passing that each eigenmode $`\text{ }\text{ }\text{ }\text{ }\text{ }F_n=m_n^2F_n`$ to this equation has two integration constants and one mass eigenvalue. One constant corresponds to the overall normalization. The remaining integration constant is fixed by the boundary condition at the Planck brane, and the mass is determined by the boundary condition on the TeV brane. In the stiff potential approximation we use the boundary condition given by (3.15).

It is then possible to show that a solution $`F`$ to the above equation automatically implies that the $`\phi `$ equation and the remaining Einstein equation are satisfied. In particular, starting with (3.16), one uses the derivative of (3.12) to eliminate $`F^{\prime \prime }`$. The resulting equation, call it $`E`$, is then differentiated and the combination $`0=E^{}2A^{}E`$ is constructed. Using the background field equations and (3.12) one arrives at the $`\phi `$ equation. Finally, the $`\delta R_{\mu \nu }`$ equation is obtained from the $`\delta R_{55}`$ equation after substituting for $`\phi ^{}`$.

## 4 General Properties of the Equation

First we show that the single ordinary differential equation for $`F(y)`$ given in (3.17) can always be brought into the Schrödinger form. For this we first transform the equations into the coordinate system where the background metric is conformally flat. This is achieved by the change of variables $`dze^{A(z)}=dy`$, where $`A(z)=A(y(z))`$. In these coordinates the equation simplifies to

$$F^{\prime \prime }3A^{}F^{}4A^{\prime \prime }F2\frac{\varphi _0^{\prime \prime }}{\varphi _0^{}}F^{}+4\frac{\varphi _0^{\prime \prime }}{\varphi _0^{}}A^{}F=m^2F.$$
(4.1)
After the rescaling of the field $`F`$ by $`F=e^{3/2A}\varphi _0^{}\stackrel{~}{F}`$ we obtain the Schrödinger-like equation

$$\stackrel{~}{F}^{\prime \prime }+\left(\frac{9}{4}A^2+\frac{5}{2}A^{\prime \prime }A^{}\frac{\varphi _0^{\prime \prime }}{\varphi _0^{}}+2\left(\frac{\varphi _0^{\prime \prime }}{\varphi _0^{}}\right)^2\frac{\varphi _0^{\prime \prime \prime }}{\varphi _0^{}}\right)\stackrel{~}{F}=m^2\stackrel{~}{F}.$$
(4.2)
However, this by itself does not guarantee hermiticity of the differential operator in (4.2). The reason is that this operator is defined only on a finite strip, and therefore in addition to writing the equation in a Schrödinger form one also has to ensure that one has hermitian boundary conditions for $`F`$. For the differential operator in (4.2) to actually be hermitian on the finite strip between the two branes, one also has to require that for any two functions $`F_1,F_2`$ on the strip $`F_1^{}(0)F_2(0)F_1(0)F_2^{}(0)F_1^{}(z_b)F_2(z_b)+F_1(z_b)F_2^{}(z_b)=0`$, where $`0`$ and $`z_b`$ denote the positions of the branes in the conformally flat $`z`$ coordinates. Once this condition is satisfied, it is automatically guaranteed by the usual theorems that all eigenvalues $`m_n^2`$ are real, that the eigenfunctions are orthogonal to each other and that they form a complete set. The actual boundary conditions that $`F`$ has to satisfy can be derived from the general boundary condition given in (3.14). In the particular model considered in this paper the boundary condition in the $`y`$ coordinates is given by

$$\pm \phi ^{}=\gamma _\pm ^2\phi \pm 2u\varphi _\pm F$$
(4.3)
In the special limit when $`\gamma _\pm \mathrm{}`$ this boundary condition reduces to $`\phi =0`$ on the two boundaries, which together with the constraint equation (3.12) between $`\phi `$ and $`F`$ just implies

$$(F^{}2A^{}F)|_i=0$$
(4.4)
at the two branes. Upon transforming to the Schrödinger basis and $`z`$ coordinates the boundary condition will be replaced by

$$\stackrel{~}{F}^{}=\stackrel{~}{F}(\frac{1}{2}A^{}\frac{\varphi _0^{\prime \prime }}{\varphi _0^{}})$$
(4.5)
at the branes. This boundary condition clearly satisfies the hermiticity properties and thus will ensure the appearance of only real mass eigenvalues of the coupled system. This will also be the case that we will analyze in full detail in the following sections. As for the general case, when $`\gamma _i`$ is finite, the boundary condition will not be hermitian. This can be easily seen from the fact that the general boundary condition involves $`\phi ^{}`$ at the branes, which should be expressed from (3.12) in terms of $`F^{\prime \prime }`$, $`F^{}`$ and $`F`$ at the brane. The appearance of $`F^{\prime \prime }`$ in the boundary condition will generically ruin the hermiticity of the operator. Nevertheless, one may eliminate $`F^{^{\prime \prime }}`$ in favor of the eigenvalue, and one can in principle solve for $`F`$. The non–hermiticity by itself however does not mean that the eigenvalues are not real. In fact, since $`\varphi `$ is a real scalar, and $`F`$ a component of the metric tensor, both of these functions have to be real to start with, which guarantees at least the appearance of only real eigenvalues. While for the model studied here (see Section 6) the radion is not tachyonic, it is unclear whether for a general potential this remains true. However, the orthogonality of the solutions is not guaranteed by anything, and will likely be violated in general for the non-hermitian boundary conditions. It would be interesting to understand the physics behind the non-orthogonality of these solutions in more detail.

## 5 Approximate solution for the KK tower

We have seen that the coupled radion-scalar system leads to a single ordinary second order differential equation. From now on we will always assume that we can use the limit $`\gamma _i\mathrm{}`$, and be able to use the hermitian boundary conditions (4.4). In the following we will present an approximate solution to these equations. For this, we will first neglect the backreaction of the non-vanishing scalar background on the metric. This will lead us to a simple Bessel-type equation, which will give a very good approximation for the masses of the KK tower of the fields. However, surprisingly, in this approximation the radion (which we identify as the lowest lying solution of (3.17) remains massless. Therefore, after presenting this approximation, we will give a perturbative analysis for the effect of the backreaction of the metric on the radion mass. We will find that as expected, the radion mass will be of order TeV, but somewhat lighter just as predicted in .

To find the actual wave functions and masses for the radion-scalar system, we will use the particular model put forward by de Wolfe et al.  and summarized in Section 2. First we neglect the backreaction of the scalar field background on the metric, which seems to be a good approximation as long as $`\kappa \varphi _{P,T}1`$. In this case the equation for the radion field $`F`$ reduces in the Schrödinger frame to:

$$F^{\prime \prime }+\frac{\alpha (\alpha +1)k^2}{(kz+1)^2}F=m^2F,$$
(5.1)
where $`\alpha `$ is given by

$$\alpha =\frac{3}{2}\frac{u}{k}.$$
(5.2)
In these coordinates the boundary conditions at the brane simplify to

$$F^{}+\frac{\alpha k}{kz+1}F=0$$
(5.3)
at the locations of the branes at $`z=0`$ and $`z_b\frac{1}{k}(e^{kr_0}1)`$. The solutions of these equations are given by linear combinations of the Bessel functions $`J_{\alpha +\frac{1}{2}}`$ and Neumann functions $`N_{\alpha +\frac{1}{2}}`$:

$$F_n(z)=a_n(z+\frac{1}{k})^{\frac{1}{2}}N_{\alpha +\frac{1}{2}}(m_n(z+\frac{1}{k}))+b_n(z+\frac{1}{k})^{\frac{1}{2}}J_{\alpha +\frac{1}{2}}(m_n(z+\frac{1}{k})).$$
(5.4)
The mass eigenvalues $`m_n`$ can be determined from the boundary condition (5.3). Using the relation for Bessel functions

$$Z_n^{}(x)=Z_{n1}(x)\frac{n}{x}Z_n(x)$$
(5.5)
the boundary conditions at the two branes simply reduce to

$`a_nN_{\alpha \frac{1}{2}}({\displaystyle \frac{m_n}{k}})+b_nJ_{\alpha \frac{1}{2}}({\displaystyle \frac{m_n}{k}})=0,`$
$`a_nN_{\alpha \frac{1}{2}}({\displaystyle \frac{m_ne^{kr_0}}{k}})+b_nJ_{\alpha \frac{1}{2}}({\displaystyle \frac{m_ne^{kr_0}}{k}})=0,`$ (5.6)
which yields the simple equation

$$b(m_n)=J_{\alpha \frac{1}{2}}(\frac{m_n}{k})\frac{N_{\alpha \frac{1}{2}}(\frac{m_ne^{kr_0}}{k})}{N_{\alpha \frac{1}{2}}(\frac{m_n}{k})}J_{\alpha \frac{1}{2}}(\frac{m_ne^{kr_0}}{k})=0,$$
(5.7)
which can be used to determine the mass eigenvalues $`m_n`$. This can be done numerically. In Fig. 1 we show the lowest mass eigenvalues for $`\alpha =2.5`$, which corresponds to the somewhat unrealistic value $`\frac{u}{k}=1`$. In Fig. 2 we show the dependence of the first non-vanishing mass eigenvalue on the value of $`\alpha =3/2u/k`$. One can easily see from (5.7), that $`m=0`$ is always a solution to (5.7), therefore in the approximation we are using the radion is still massless. For the higher states of the KK tower it is a good approximation to use the mass eigenvalues obtained from (5.7), because the masses are of the order (and even larger) than the TeV scale, thus in the limit of small backreaction that we are considering throughout the paper these masses will be only slightly modified. The radion (which appeared as the zero mode above) however needs special treatment, since the shift in the mass (which is usually negligible for the higher KK modes) coming from the backreaction of the metric background due to the scalar field is the leading order contribution to the mass for the radion. Below we will estimate the size of the radion mass in perturbation theory.

## 6 The radion mass

In the previous section we have seen what the approximate wave-functions and masses are for the KK tower of the coupled radion-scalar system. In this approximation of neglecting the backreaction, however, we have still found a vanishing radion mass. This is in fact easy to show for a general stabilizing potential. From (3.17), $`F=e^{2A}`$ is always a solution with zero mass if $`A^{\prime \prime }`$ is neglected in the bulk. Thus the radion mass is always proportional to the backreaction of the metric independently of the details of the potential of the stabilizing scalar field. In the following, we will show how the backreaction generates a non-vanishing mass for the radion field. For this, we start with the equation describing the radion wave function in the $`y`$ coordinates:

$$F^{\prime \prime }2A^{}F^{}4A^{\prime \prime }F+2uF^{}4uA^{}F+m^2e^{2A}F=0,$$
(6.1)
where $`A(y)`$ is given in (2.19). The appropriate boundary condition is $`F^{}2A^{}F=0`$ at the branes. In the special limit $`\gamma _\pm \mathrm{}`$ the other boundary condition is $`\phi =0`$. Thus we will treat the backreaction as a perturbation, and look for the solution in terms of a perturbative series in $`l\kappa \varphi _P/\sqrt{2}`$. Then we write the solution as

$$F_0=e^{2k|y|}(1+l^2f_0(y))\text{ , }m_r^2=l^2\stackrel{~}{m}^2\text{ , }A(y)=k|y|+\frac{l^2}{6}e^{2u|y|}.$$
(6.2)
Expanding the solution as above and keeping only the leading terms in $`l^2`$ we obtain the equation

$$f_0^{\prime \prime }+2(k+u)f_0^{}=\stackrel{~}{m}^2e^{2k|y|}\frac{4}{3}(ku)ue^{2u|y|}$$
(6.3)
along with the boundary conditions

$$f_0^{}+\frac{2}{3}ue^{2u|y|}=0$$
(6.4)
at the location of the branes. One can easily find the most general solution for $`f`$ from the equation in the bulk, which is given by

$$f_0^{}(y)=Ce^{2(k+u)|y|}\frac{\stackrel{~}{m}^2}{2(2k+u)}e^{2k|y|}\frac{2(ku)u}{3k}e^{2u|y|},$$
(6.5)
where the integration constant $`C`$ along with the radion mass $`\stackrel{~}{m}`$ is determined by the boundary conditions at the brane. This way we obtain the radion mass to be

$$m_{radion}^2=\frac{4l^2(2k+u)u^2}{3k}e^{2(u+k)r_0},$$
(6.6)
where $`r_0`$ denotes the location of the brane. Note that this result is very similar to the answer obtained from the effective theory computation using the naïve ansatz, except for the important difference in the power of $`u/k`$. The exact result obtained here scales as $`(u/k)^2`$, whereas the effective theory result would scale as $`(u/k)^{3/2}`$.***We thank Jim Cline for these observations. It would be very interesting to understand the origin of this different scaling. For this model to give the correct value of the weak scale without reintroducing a large fine-tuning again one needs $`u/k1/37`$, thus the radion mass turns out to be somewhat lighter than the TeV scale. It is suppressed by the factor $`l\frac{u}{k}e^{ur_0}`$ compared to the TeV scale. Thus in this approximation $`m_{radion}\frac{l}{40}`$ TeV, which could be at least in the range of a few GeV’s. Of course, we need to emphasize that $`l`$ is not necessarily small for the stabilization mechanism to work, we took this limit only for calculational convenience.

## 7 Coupling to SM fields

In this section the coupling of the radion and KK tower of $`\varphi `$ to the TeV brane are obtained. In particular we demonstrate that the bulk scalar field gives a small correction to the radion kinetic term, and thus the kinetic terms obtained from the Einstein–Hilbert part of the action dominate, justifying the results obtained using the naïve ansatz .

In the previous section it was seen that by including the backreaction an $`𝒪(`$TeV$`{}_{}{}^{2})`$ for the radion is obtained. The wavefunction is then

$$F_0(x,y)=e^{2k|y|}\left(1+l^2f_0(y)\right)R(x)$$
(7.1)
where $`f_0(y)`$ is given by the integral of (6.5). Since the radion mass is $`𝒪(`$TeV$`{}_{}{}^{2})`$, and by assumption $`l^21`$, we see by inspection that the backreaction induces a small correction to the unperturbed wavefunction. So for the purposes of determining the coupling of the radion to the TeV brane it is sufficient to include only the unperturbed wavefunction, namely $`F(x,y)=e^{2k|y|}R(x)`$. Then a straightforward calculation gives

$$M^3𝑑y\sqrt{g}6M^3(R)^2e^{2A}e^{4k|y|}=\frac{6M^3}{k}(e^{2kr_0}1)(R)^2$$
(7.2)
So the normalized radion $`r(x)`$ is $`R(x)=r(x)e^{kr_0}/\sqrt{6}M_{Pl}`$, since $`M^3/k=M_{Pl}^2/2`$. This implies a coupling to the TeV brane fields which is

$$R(x)e^{2kr_0}\left(1+𝒪(l^2)\right)\text{Tr}T_{\mu \nu }=\frac{r(x)}{\sqrt{6}M_{Pl}e^{kr_0}}\text{Tr}T_{\mu \nu }\left(1+𝒪(l^2)\right),$$
(7.3)
where the left hand side of this equation is a consequence of the fact that the induced metric on the TeV brane is $`e^{2A(r_0)}(1e^{2kr_0}R(x)\eta _{\mu \nu })`$. The coupling obtained this way agrees precisely with (2.4). This is perhaps surprising, since the latter computation used an ansatz which did not satisfy the equations of motion. This makes us suspect that results which depend only on the leading unperturbed form of the radion wave function will be correctly captured by the naïve ansatz.

Now we address the issue that was originally raised by . Is the radion kinetic term dominated by the kinetic term of the bulk scalar field or the bulk gravity action, and in particular is the former hierarchically larger? To answer this, we need the change in $`\phi `$ caused by a fluctuation in the radion. Since $`\phi =0`$ when the backreaction is not included, we must include the leading backreaction correction to the radion wavefunction, given by the integral of (6.5). From (3.12) we compute that the change in $`\phi `$ to $`𝒪(l^2)`$ due to the radion is

$$\phi =\frac{3}{\kappa ^2\varphi _0^{}}(F^{}2A^{}F)=\frac{3l^2}{\kappa ^2\varphi _0^{}}(F_0^{}+\frac{2u}{3}e^{2u|y|})=\frac{3l^2R(x)}{\kappa ^2\varphi _0^{}}e^{2k|y|}f_3(y)$$
(7.4)
where

$`f_3(y)`$ $``$ $`f_0^{}(y)+{\displaystyle \frac{2}{3}}ue^{2u|y|}`$ (7.5)
$`=`$ $`Ce^{2(k+u)|y|}{\displaystyle \frac{\stackrel{~}{m}^2}{2(2k+u)}}e^{2k|y|}+{\displaystyle \frac{2u^2}{3k}}e^{2u|y|}.`$
This fluctuation in $`\varphi `$ then contributes to the radion kinetic term at $`O(l^2)`$ an amount

$$𝑑ye^{4A}g^{\mu \nu }_\mu \phi _\nu \phi =\frac{9l^2}{2\kappa ^2u^2}(R)^2𝑑ye^{2(ku)|y|}f_3(y)^2.$$
(7.6)
From (7.2) the unnormalized contribution from bulk gravity to the kinetic term is $`e^{2kr_0}`$. So we only need to consider those contributions from $`\varphi `$ which are comparable or larger to this. Recalling that $`m^2e^{2kr_0}`$, it is seen that the largest terms in (7.6) are at best $`e^{2kr_0}`$. Explicitly performing the integral one finds that it is

$$\delta =2l^2\frac{u^2}{\kappa ^2k^2}(R)^2e^{2kr_06ur_0}\left(\frac{1}{3ku}+\frac{1}{k3u}\frac{1}{ku}\right)\left(1+O(l^2)\right).$$
(7.7)
This is typically $`l^2u^2e^{2kr_0}M^3/k^3`$, which is smaller than (7.2) since we assuming that the backreaction is small, $`l1`$, and also that $`uk`$ to obtain a realistic hierarchy. So the radion kinetic term is dominated by the contribution from the bulk gravity, and receives a small correction from the stabilizing bulk scalar field.

In Section 4 it was found that for the simple boundary conditions $`\phi =0`$ (corresponding to the limit $`^2\lambda _\pm /\varphi ^21`$) a self–adjoint equation for $`F`$ was obtained. The general solution to this is

$$F(x,y)=\underset{n}{}\alpha _nF_n(x,y)$$
(7.8)
where $`F_n`$ is a mass eigenstate, and the $`\alpha _n^{}s`$ are some numbers. We expect that $`F`$ includes the massive radion, but where did all the other states come from? It is helpful to reconsider what happens when the backreaction is neglected. In this limit the KK tower in $`F`$ completely disappears and only the (massless) radion remains. This may be observed from (3.12), since neglecting the backreaction corresponds to $`\kappa ^2\varphi _01`$, and this amounts to setting $`F^{}2A^{}F=0`$. The only solution for $`F`$ in this case is the radion zero mode $`F=e^{2k|y|}`$. Once the backreaction is included, however, the fluctuating modes in $`\varphi `$ and $`F`$ are correlated through (3.12). In particular, a general fluctuation $`\phi `$ induces a change in $`F`$. The sum over KK states appearing above is then just the decomposition of $`F`$ into these KK eigenstates. It is then expected that the coefficients $`\alpha _n`$ for the non–radion states to be suppressed by the backreaction.

The preceding remarks imply that the TeV brane fields, which couple to the induced metric $`F`$, also directly couple to the KK tower, by an amount suppressed by the backreaction.***The KK modes of the scalar field do not mix with the KK modes of the graviton. The reason is that the only way they could mix is by a coupling of the 5D trace of the metric to the scalar KK modes. However, the graviton is traceless, and the trace is basically identified with the radion, therefore no additional graviton-scalar mixing could be introduced. Since $`F\varphi _0^{}\phi `$ is already suppressed by the backreaction, therefore in order to compute the induced metric to lowest order in the backreaction we can use the zeroth order wavefunctions for $`\phi `$. The normalized KK fields are given by

$$\phi _n(x,z)=\frac{\psi _n(x)}{N_n}(kz+1)^2J_{2+u/k}(m_nz+1/k)\left(1+l^2f_n(y)\right).$$
(7.9)
Here $`\psi _n(x)`$ are the normalized 4–D fields satisfying $`\text{ }\text{ }\text{ }\text{ }\text{ }\psi _n=m_n^2\psi _n`$. The orthogonality of these solutions when the backreaction vanishes $`(l=0)`$ follows from the boundary condition $`\phi _n=0`$ and the properties of the Bessel functions. Also $`\psi _n(x)`$ are the normalized 4–D fields satisfying $`\text{ }\text{ }\text{ }\text{ }\text{ }\psi _n=m_n^2\psi _n`$. The normalization constant is

$$N_n=\frac{1}{\sqrt{k}}e^{kr_0}J_{3+u/k}\left(\frac{m_n}{k}e^{kr_0}\right).$$
(7.10)
As discussed previously, the lowest order masses $`m_n`$ are determined by $`J_{2+u/k}(e^{2kr_0}m_n/k)=0`$ and are real since the operator equation with these boundary conditions is self-adjoint. The coupling of these fields to the TeV brane is given by

$$F(x,y=r_0)\text{Tr}T_{\mu \nu }$$
(7.11)
where $`F`$ is the solution of (3.12) for the solutions $`\phi _n`$ given above. One finds the coupling

$$\frac{l}{m_n}\lambda _sx_n^u\psi _n(x)\text{Tr}T_{\mu \nu }.$$
(7.12)
The model-dependent couplings that appear are

$$\lambda _S=\frac{\sqrt{2}}{3}u\kappa \sqrt{k}e^{ur_0}𝒪(\frac{u}{k})$$
(7.13)
and

$$x_n^u=\frac{J_{1+u/k}(\frac{m_n}{k}e^{kr_0})}{J_{3+u/k}(\frac{m_n}{k}e^{kr_0})}$$
(7.14)
is a numerical constant of $`O(1)`$. While the inclusion of the backreaction leads to a TeV suppressed coupling for the KK modes, the size decreases rather rapidly due to the $`1/m_n1/`$TeV suppression, as may be observed from inspecting Fig. 1.

The coupling discussed here implies that the KK modes of $`\varphi `$ can be directly produced at future colliders, and they also decay directly to Standard Model fields. This may be puzzling at first, since the stabilizing potential may have a global discrete symmetry, such as $`Z_2`$, which would naïvely imply that some of these KK modes are stable. The background vev for $`\varphi `$ explicitly breaks this symmetry, however, and this allows for all the KK modes to decay into the brane world fields.

The direct coupling of the KK modes from the stabilizing fields may have interesting implications for search strategies and current limits on the Randall-Sundrum framework. In particular, it may be important to not neglect the stabilizing potential when discussing these issues. However, when the backreaction is small, the size of their couplings is suppressed by $`u/k1/40`$ compared to the that of the radion. Therefore, in what follows, we neglect these states in the loop computations.

## 8 Cosmological Implications

The subject of brane cosmology has recently attracted lots of interest . Most of this was due to the realization, that the expansion of a brane universe could be significantly different from the ordinary Friedmann–Robertson–Walker (FRW) cosmology . However, it did not take very long to realize, that this is simply due to the fact, that a generic brane model (like the one presented in ) can not give the ordinary cosmological evolution, since gravity is in general manifestly higher dimensional. This means that in these models the 4D effective theory is usually not described by ordinary Einstein gravity, but generically a complicated scalar-tensor theory of gravity. However, observations show that our Universe is described by Einstein’s theory of relativity to a high precision, therefore one has to require from the outset that a brane model reproduces ordinary Einstein gravity, at least at long enough distances. Once this is achieved, the cosmological expansion will be automatically described by the ordinary FRW Universe, which simply follows from the fact that the effective theory is ordinary Einstein gravity. Thus one can see that the issue of unconventional cosmologies is nothing else but the issue of whether one recovers 4D gravity. This issue manifests itself in the case of the Randall-Sundrum two-brane model due to the presence of the radion field. Without a stabilizing potential, the radion field will be massless, and yield additional long range forces, and also contribute to the expansion of the Universe, yielding an unconventional cosmology, which is presumably excluded by the requirement for a successful nucleosynthesis . Thus the radion field has to obtain a mass. Once it is massive, gravity on both branes will be ordinary 4D gravity, and thus the cosmology will be conventional below temperatures comparable to the radion mass. This has been explained in great detail in , and also in . In  a simplified calculation for the cosmological expansion has been presented, where the wave function of the radion has been neglected, and also the effects of the stabilizing scalar field were included by adding a five dimensional potential for the radion field $`V(b)`$. Assuming that that the potential $`V(b)`$ is very steep, it was shown from a perturbative solution of the bulk equations that the ordinary FRW Universe is recovered. It was also argued, that the $`55`$ component of the Einstein equation, which in the absence of a stabilizing potential usually leads to the unconventional expansion equations will only determine the shift in the radion field due to matter on the wall, and does not result in unconventional cosmologies once the radius is stabilized. Below we demonstrate, that the results of  which were neglecting the radion wave function, and also did not include the fluctuations of the scalar field at the brane remain valid in the more precise framework of radion stabilization explained in the previous sections. In particular, we will show, that the result obtained for the shift in the radion field due to matter on the walls in Eq. (4.15) of  is exactly reproduced in the full calculation.

To compute $`G_{55}`$ we use the ansatz

$$ds^2=n(t,y)^2dt^2a^2(t,y)d^2xb(t,y)^2dy^2$$
(8.1)
for which

$$G_{55}=3\left[\frac{a^{}}{a}\left(\frac{a^{}}{a}+\frac{n^{}}{n}\right)\frac{b^2}{n^2}\left(\frac{\dot{a}}{a}\left(\frac{\dot{a}}{a}\frac{\dot{n}}{n}\right)\frac{\ddot{a}}{a}\right)\right].$$
(8.2)
The jump equations for $`a`$ and $`n`$ on the TeV brane imply

$$\frac{[a^{}]}{a}=\frac{\kappa ^2}{3}(\lambda _{}(\varphi )+\rho )b\text{ , }\frac{[n^{}]}{n}=\frac{\kappa ^2}{3}(\lambda _{}(\varphi )+3p+2\rho )b$$
(8.3)
Here $`\rho `$ and $`p`$ are the bare energy matter density and pressure on the TeV brane, which are related to the physically measured quantities on the TeV brane by $`\rho _0=\rho e^{4A_0}`$, etc. where $`e^{A_0}`$ is the scale factor on the TeV brane. Then averaging the $`G_{55}`$ equation about the TeV brane and linearizing to $`𝒪(\rho ,F,\phi )`$ gives

$`<G_{55}>`$ $`=`$ $`\kappa ^4{\displaystyle \frac{\lambda _{}^2(\varphi _0)}{6}}3e^{2A}\left(\left({\displaystyle \frac{\dot{a}}{a}}\right)^2+{\displaystyle \frac{\ddot{a}}{a}}\right){\displaystyle \frac{\kappa ^4\lambda _{}(\varphi _0)}{12}}(3p\rho )`$ (8.4)
$`+{\displaystyle \frac{\kappa ^4\lambda _{}(\varphi _0)}{3}}{\displaystyle \frac{\lambda _{}(\varphi _0)}{\varphi }}\phi +{\displaystyle \frac{2\kappa ^4\lambda _{}^2(\varphi _0)}{3}}F.`$
Terms with $`\dot{n}`$ are higher order in $`\rho `$ and are dropped, and $`b=1+2F`$ has been used. For late–time cosmology in the presence of radion stabilization it is reasonable to use the FRW equation

$$\left(\frac{\dot{a_0}}{a_0}\right)^2+\frac{\ddot{a_0}}{a_0}=\frac{1}{6M_{Pl}^2}\left(3p_{}\rho _{}+3p_0\rho _0\right).$$
(8.5)
Note this includes a contribution from matter $`(\rho _{})`$ on the Planck brane. Also implicit in the use of this equation is the assumption that the time variation of the radion is negligible, which is justified a posteriori. Then using the relation $`\kappa ^4\lambda _{}(\varphi _0)=6(1e^{2A_0})/M_{Pl}^2`$ gives

$`<G_{55}>`$ $`=`$ $`{\displaystyle \frac{\kappa ^4\lambda _{}^2(\varphi _0)}{6}}+{\displaystyle \frac{e^{4A_0}}{2M_{Pl}^2}}\left(3p_0\rho _0+e^{2A_0}(3p_{}\rho _{})\right)`$ (8.6)
$`+{\displaystyle \frac{\kappa ^4\lambda _{}(\varphi _0)}{3}}{\displaystyle \frac{\lambda _{}(\varphi _0)}{\varphi }}\phi +{\displaystyle \frac{2\kappa ^4\lambda _{}^2(\varphi _0)}{3}}F.`$
The $`G_{55}`$ equation is

$$G_{55}=\kappa ^2T_{55}=\kappa ^2\left(\frac{1}{2}\varphi _{}^{}{}_{}{}^{2}g_{55}(\frac{1}{2}(\varphi )^2V(\varphi ))\right).$$
(8.7)
Then the averaging of $`T_{55}`$ and linearizing using (3.2) gives

$$\kappa ^2<T_{55}>=\kappa ^2(\frac{1}{2}\varphi _{}^{}{}_{0}{}^{2}V(\varphi _0))+\kappa ^2(\varphi _0^{}\phi ^{}4FV\frac{V}{\varphi }\phi )$$
(8.8)
with all quantities are evaluated on the TeV brane. Using the background bulk equation (2.10) and the jump equation (2.12) the leading terms are seen to cancel. Then after using the background equations, the jump equations for the background fields, and some algebra gives

$$\frac{e^{4A_0}}{2M_{Pl}^2}\left(3p_0\rho _0+e^{2A_0}(3p_{}\rho _{})\right)=\kappa ^2(\varphi _{}^{}{}_{0}{}^{}\phi ^{}2\varphi _{}^{^{}}{}_{0}{}^{2}F\frac{V}{\varphi }\phi 4A^{}\varphi _0^{}\phi ).$$
(8.9)
Using (3.16) to eliminate $`\phi ^{}`$, (3.17) to eliminate $`F^{\prime \prime }`$ in favor of the mass eigenvalue, and (2.11) to eliminate $`\varphi _0^{\prime \prime }`$ finally gives

$$\frac{e^{4A_0}}{2M_{Pl}^2}\left(3p_0\rho _0+e^{2A_0}(3p_{}\rho _{})\right)=3e^{2A_0}m_r^2F.$$
(8.10)
But the shift in the distance between the two branes is obtained by integrating the line element, which gives $`\delta r_0=R(e^{2kr_0}1)/k`$. Then since $`F=Re^{2kr_0}`$, one obtains

$$\frac{\delta r_0}{r_0}=\frac{1}{6kr_0}\frac{(1e^{2A_0})}{m_r^2M_{Pl}^2e^{2A_0}}\left(\rho _03p_0+e^{2A_0}(\rho _{}3p_{})\right)$$
(8.11)
which is precisely the result found in  obtained by using a 4D effective theory.***A translation dictionary between two different notations is required: $`kr_0=m_0b_0/2`$. This is perhaps not surprising, since for constant radion field the naïve ansatz and the full metric including the wavefunction of the radion are equivalent up to a coordinate transformation. So in an adiabatic approximation the leading order result using the naïve ansatz should agree with that obtained from using the correct radion wavefunction, if the fluctuations in the scalar field are ignored. It is less clear why the full answer including the contribution from the scalar field turns out to be exactly equal to the calculation using the naïve ansatz. Note that matter on the Planck brane causes a smaller shift in the radion compared to an equal amount of matter on the TeV brane. This is because the radion wavefunction is peaked at the TeV brane, and it couples more weakly to the Planck brane relative to the TeV brane by precisely the amount $`e^{2A_0}`$. Thus one finds the very general result that in the presence of matter on the branes and a stabilizing mechanism, the $`G_{55}`$ equation determines the shift in the radion.

## 9 The Effective 4D Lagrangian

In the previous sections we have argued that in the presence of a stabilizing potential the linear couplings of the radion and bulk scalars is given by

$`{\displaystyle \frac{1}{2}}(r)^2{\displaystyle \frac{1}{2}}m^2r^2+{\displaystyle \underset{n}{}}{\displaystyle \frac{1}{2}}\left((\psi _n)^2m_n^2\psi _n^2\right)+DH^{}DH`$
$`+\left({\displaystyle \frac{r(x)}{\sqrt{6}\mathrm{\Lambda }}}+{\displaystyle \underset{n}{}}\alpha _n{\displaystyle \frac{\psi _n(x)}{\mathrm{\Lambda }_n}}\right)\text{Tr}T_{\mu \nu }+\xi H^{}HV(H).`$ (9.1)
The masses appearing here are $`O(`$TeV$`)`$, and their particular value depends on the details of the stabilizing mechanism. The scale $`\mathrm{\Lambda }=e^{kr_0}M_{Pl}`$ in the Randall-Sundrum model, but here we have left it general. The other scales are $`\mathrm{\Lambda }_nm_n`$, and the $`\alpha _n`$ are also model–dependent, and vanish in the limit of small backreaction. In the remaining sections we restrict ourselves to the above Lagrangian, and do not commit ourselves to any specific mechanism of radius stabilization. For the electroweak analysis we neglect the contributions of the KK modes from $`\varphi `$.

Note that in the above Lagrangian we have also included a curvature Higgs scalar operator $`H^{}H`$. The presence of this operator leads to interesting signals for discovering the Higgs and radion at future colliders . In particular, the branching fractions of the Higgs and radion to $`gg`$ and $`\overline{b}b`$ can be substantial different from that of the SM Higgs.

As discussed in , the presence of the conformal term $`H^{}HR`$ leads to both kinetic and mass mixing between the neutral Higgs and radion. Below we summarize the relevant formulae for mixing and couplings. The interested reader is referred to the next section for details.

One finds that the “gauge” $`h`$ and $`r`$ are related to the mass eigenstates $`h_m`$ and $`r_m`$ by

$$h=(\mathrm{cos}\theta \frac{6\xi \gamma }{Z}\mathrm{sin}\theta )h_m+(\mathrm{sin}\theta +\frac{6\xi \gamma }{Z}\mathrm{cos}\theta )r_m\text{ , }$$
(9.2)
$$r=\mathrm{cos}\theta \frac{r_m}{Z}\mathrm{sin}\theta \frac{h_m}{Z}\text{ ,}$$
(9.3)
where

$$\mathrm{tan}2\theta =12\xi \gamma Z\frac{m_h^2}{m_r^2m_h^26\xi \gamma ^2(112\xi )}\text{ , }$$
(9.4)
$$\gamma =\frac{v}{\sqrt{6}\mathrm{\Lambda }}\text{ , }Z^2=1+6\xi \gamma (16\xi ),$$
(9.5)
where $`v246`$ GeV is the electroweak vev. Requiring that the quantity $`Z^2`$ be positive (in order to avoid ghost-like states) places an upper bound on the value of $`\xi `$, for a given $`\gamma `$. Physically this requirement comes from maintaining positive definite kinetic terms for $`h`$ and $`\varphi `$.

In this basis, the couplings of the physical radion and Higgs appropriate for tree–level studies are

$$\left((\mathrm{cos}\theta (6\xi 1)\frac{\gamma \mathrm{sin}\theta }{Z})h_m+(\mathrm{sin}\theta +(6\xi 1)\frac{\gamma \mathrm{cos}\theta }{Z})r_m\right)\text{Tr}T_{\mu \nu }.$$
(9.6)
In the $`\xi 0`$ limit one recovers

$$\left(h\gamma r\right)\text{Tr}T_{\mu \nu }\text{ , }$$
(9.7)
obtained in . Note that $`\text{Tr}T_{\mu \nu }`$ includes SM Higgs contributions.

## 10 Curvature–Scalar Mixing

In this Section the effect of introducing a curvature–scalar interaction are reviewed. The discussion parallels , however some of the resulting formulae are slightly different because here terms of $`𝒪(\gamma ^2)`$ and $`𝒪(\gamma ^2\xi ^2)`$ are kept.

We begin with the couplings of the radion and Higgs to the SM fields before electroweak symmetry breaking. The induced metric on the TeV wall is

$$g_{\mu \nu }^{ind}(x)=e^{2A(r_0)2e^{2kr_0}R(x)}g_{\mu \nu }(x),$$
(10.1)
where the warp factor includes the backreaction, although its inclusion is not necessary for our purposes. The canonically normalized radion $`r`$ is

$$R(x)=e^{kr_0}\frac{r(x)}{\sqrt{6}M_{Pl}}.$$
(10.2)
So we express the induced metric expanded about a Minkowski metric as

$$g_{\mu \nu }^{ind}(x)=e^{2A(r_0)2\frac{\gamma }{v}r(x)}\eta _{\mu \nu }e^{2A(r_0)}\mathrm{\Omega }^2(r)\eta _{\mu \nu }$$
(10.3)
with

$$\gamma =\frac{v}{\sqrt{6}\mathrm{\Lambda }}\text{ , }\mathrm{\Lambda }=M_{Pl}e^{kr_0}.$$
(10.4)
The four dimensional effective action we consider is

$`S_{TeV}`$ $`=`$ $`{\displaystyle d^4x\sqrt{g_{ind}}\left(g_{ind}^{\mu \nu }D_\mu H^{}D_\nu HV(H)\right)}`$ (10.5)
$`+{\displaystyle d^4x\sqrt{g}\frac{1}{2}\left((r)^2m_r^2r^2\right)}+{\displaystyle d^4x\sqrt{g_{ind}}\xi R(g_{ind})H^{}H}+S_{SM}.`$
To canonically normalize the Higgs and other SM fields, we perform the field–independent redefinition

$$He^{A(r_0)}H\text{ , }\psi e^{3A(r_0)/2}\psi $$
(10.6)
In this basis the Higgs–radion potential is

$$V(H,r)=\mathrm{\Omega }^4(r)V(H).$$
(10.7)
Note that $`V`$ also includes the effective 4–dimensional cosmological constant, which we assume to vanish. Clearly this potential has a minimum at the same location as $`V(H)`$, so that the EWSB vacuum is $`r=0`$ and $`H^0=v/\sqrt{2}`$.

We consider the presence of the curvature mixing term

$$_\xi =\sqrt{g_{ind}}\xi (g_{ind})H^{}H.$$
(10.8)
Our choice of signs for $`\xi `$ is such that the Higgs potential receives a positive mass–squared correction in a de Sitter phase when $`\xi `$ is positive. Since this is a renormalizable interaction, there is no reason for it not to be present, or to be suppressed. What makes this operator important in this case is that $``$ contains the induced metric, rather than just the ordinary 4-dimensional metric. In particular,

$$(\mathrm{\Omega }^2(r)\eta _{\mu \nu })=6\mathrm{\Omega }^2\left(\text{ }\text{ }\text{ }\text{ }\text{ }\mathrm{ln}\mathrm{\Omega }+(\mathrm{ln}\mathrm{\Omega })^2\right).$$
(10.9)
So the curvature–scalar interaction is

$$_\xi =6\xi \mathrm{\Omega }^2\left(\text{ }\text{ }\text{ }\text{ }\text{ }\mathrm{ln}\mathrm{\Omega }+(\mathrm{ln}\mathrm{\Omega })^2\right)H^{}H$$
(10.10)
To see the effect of the curvature scalar interaction we expand $`H^0=(v+h)/\sqrt{2}`$ and $`\mathrm{\Omega }(r)=1\gamma r/v+\mathrm{}`$. We need to only expand $`\mathrm{\Omega }`$ to linear order since the derivative terms are already of $`O(r)`$. This gives at quadratic order

$$_\xi =6\xi \gamma h\text{ }\text{ }\text{ }\text{ }\text{ }r+3\xi \gamma ^2(r)^2.$$
(10.11)
where a total derivative has been dropped. The $`\xi `$ terms clearly introduce kinetic mixing. The full radion–Higgs Lagrangian to be diagonalized is

$$=\frac{1}{2}h\text{ }\text{ }\text{ }\text{ }\text{ }h\frac{1}{2}m_h^2h^2\frac{1}{2}(1+6\xi \gamma ^2)r\text{ }\text{ }\text{ }\text{ }\text{ }r\frac{1}{2}m_r^2r^2+6\xi \gamma h\text{ }\text{ }\text{ }\text{ }\text{ }r.$$
(10.12)
The mass parameters $`m_r`$, $`m_h`$ are the masses of the radion and Higgs, respectively, in the limit $`\xi =0`$. The kinetic terms are diagonalized by the shift $`h=h^{}+6\xi \gamma r^{}/Z`$ , and $`r=r^{}/Z`$. Here

$$Z^2=1+6\xi \gamma ^2(16\xi )$$
(10.13)
is the coefficient of the radion kinetic term after undoing the kinetic mixing, and is therefore required to be positive in order to keep the radion kinetic term positive definite. For a fixed cutoff $`\mathrm{\Lambda }`$ this restricts the size of the mixing parameter $`\xi `$. It must lie in the range

$$\frac{1}{12}(1\sqrt{1+\frac{4}{\gamma ^2}})\xi \frac{1}{12}(1+\sqrt{1+\frac{4}{\gamma ^2}})$$
(10.14)
for non-zero values of $`\gamma `$. Otherwise one has a ghost-like radion field, which presumably signals an instability of the theory.

This rescaling diagonalizes the kinetic terms, but introduces mixing in the mass matrix. A final rotation $`h^{}=\mathrm{cos}\theta h_m+\mathrm{sin}\theta r_m`$ and $`r^{}=\mathrm{cos}\theta r_m\mathrm{sin}\theta h_m`$ brings the Lagrangian to canonical form. With the above definition of the sign of the rotation, the rotation angle is

$$\mathrm{tan}2\theta =12\xi \gamma Z\frac{m_h^2}{m_r^2m_h^2(Z^236\xi ^2\gamma ^2)}.$$
(10.15)
We note that for moderate values of $`\xi `$ and $`\gamma `$ (i.e., $`Z^2>36\xi ^2\gamma ^2`$) the mixing angle $`\mathrm{tan}2\theta `$ is negative when $`m_h>m_r`$. For small $`\gamma `$ we can expand

$$\mathrm{tan}2\theta =12\xi \gamma \frac{m_h^2}{m_r^2m_h^2}+O(\gamma ^2).$$
(10.16)
Putting everything together, the relation between the gauge and mass eigenstates is

$`h`$ $`=`$ $`(\mathrm{cos}\theta {\displaystyle \frac{6\xi \gamma }{Z}}\mathrm{sin}\theta )h_m+(\mathrm{sin}\theta +{\displaystyle \frac{6\xi \gamma }{Z}}\mathrm{cos}\theta )r_m\text{ ,}`$ (10.17)
$`r`$ $`=`$ $`\mathrm{cos}\theta {\displaystyle \frac{r_m}{Z}}\mathrm{sin}\theta {\displaystyle \frac{h_m}{Z}}.`$ (10.18)
The mass eigenvalues are easily obtained

$$m_\pm ^2=\frac{1}{2Z^2}\left(m_r^2+(1+6\xi \gamma ^2)m_h^2\pm ((m_r^2m_h^2(1+6\xi \gamma ^2))^2+144\gamma ^2\xi ^2m_r^2m_h^2)^{1/2}\right)$$
(10.19)
The heavier state $`(+)`$ is identified with the state with the larger of $`(m_h^2,m_r^2)`$.

## 11 Radion Couplings and Feynman Rules

In this section we derive the Feynman rules relevant to the computation of the oblique parameters $`S`$, $`T`$, $`U`$.

Before proceeding, we pause to ask whether higher-order couplings such as

$$\varphi ^2\text{Tr}T_{\mu \nu }$$
(11.1)
also affect in particular the electroweak precision measurements $`S`$, $`T`$ and $`U`$. This operator could either be directly present, or generated from the above linear coupling due to a non-trivial kinetic term for the radion. Although this operator contributes at one-loop to the gauge boson two point functions, it is easy to see that they do not contribute to $`T`$, since the $`m_V^2`$ contained in $`\text{Tr}T_{\mu \nu }`$ is canceled by the $`1/m_V^2`$ appearing in the expression for $`T`$, nor to $`S`$ or $`U`$ since the contribution of this operator to the vacuum polarizations is momentum–independent. Thus we need to only consider the linear coupling

$$\frac{\gamma r}{v}\text{Tr}T.$$
(11.2)
This operator will have a contribution to the oblique corrections similar to that of the standard model Higgs. First we discuss the Feynman rules for the interactions from the $`\frac{\gamma r}{v}\mathrm{Tr}T`$ operator. The interaction Lagrangian term relevant for the gauge-boson propagator corrections is just given by

$$_{int}=\frac{\gamma }{v}r(2M_W^2W_\mu ^+W^\mu +M_Z^2Z_\mu Z^\mu ).$$
(11.3)
In addition, to ensure gauge invariance of the results, one also has to examine the gauge fixing terms carefully. The gauge fixing Lagrangian for the $`W`$ and the $`Z`$ are given by

$$_{gf}=\sqrt{g}[\frac{1}{\alpha }(D_\mu W^{\mu +}+i\alpha M_W\mathrm{\Psi }^+)(D_\mu W^\mu i\alpha M_W\mathrm{\Psi }^{})\frac{1}{2\alpha }(D_\mu Z^\mu +\alpha M_Z\mathrm{\Psi })^2,]$$
(11.4)
where the $`\mathrm{\Psi }`$’s are the would-be-Goldstone bosons, and $`\alpha `$ is the gauge fixing parameter in the Rα gauge. Note, that since the gravitational background is non-trivial, we have replaced the ordinary derivatives by covariant derivatives. The background metric is given by $`g_{\mu \nu }=\mathrm{\Omega }^2(r)\eta _{\mu \nu }=e^{\frac{2\gamma r}{v}}\eta _{\mu \nu }`$, therefore the covariant derivative of a vector will take the form

$$D_\mu V^\mu =\mathrm{\Omega }^2(_\mu V^\mu \frac{2\gamma }{v}_\mu rV^\mu ).$$
(11.5)
Thus, from the gauge fixing terms one also obtains three-point interaction vertices of the form

$$_{gf}=\frac{2\gamma }{v\alpha }(_\mu Z^\mu )(_\nu rZ^\nu )+\frac{2\gamma }{v\alpha }(_\mu W^{\mu +})(_\nu rW^\nu )+\frac{2\gamma }{v\alpha }(_\mu W^\mu )(_\nu rW^{\nu +}).$$
(11.6)
With these operators added, the Feynman rule for the three-point function is given by

$$\frac{i2M_V^2\gamma }{v}\eta _{\mu \nu }\frac{i2\gamma }{\alpha v}(p_{2\mu }p_{1\nu }+p_{3\nu }p_{1\mu })$$
(11.7)
In addition to the cubic vertices evaluated above, there are also 4-point couplings of the radion and the gauge bosons. These terms will not contribute to the oblique electroweak corrections, however they will be important to obtain a gauge invariant answer for the gauge boson vacuum polarization diagrams. These terms arise from two different sources. The first source is the conformal coupling $`e^{\frac{\gamma r}{v}}\mathrm{Tr}T`$ to the trace of the energy momentum tensor. In the formalism of  this can be obtained by a very careful expansion of the interaction terms $`|D_\mu H|^2+\frac{1}{\sqrt{6}\mathrm{\Lambda }}r(H^{}DH+h.c.)`$. Either way one finds the additional operators

$$\frac{\gamma ^2}{v^2}r^2(\frac{M_W^2}{2}W_\mu ^+W^\mu +M_Z^2Z_\mu Z^\mu ).$$
(11.8)
The other source of quartic interaction terms are again the gauge fixing terms. One simply expands these to higher order to obtain the interaction terms

$$\frac{2\gamma ^2}{\alpha v^2}_\mu r_\nu r(Z^\mu Z^\nu +2W^{\mu +}W^\nu ).$$
(11.9)
The interaction terms involving the would-be Goldstone bosons just contribute a total derivative, and thus they can be omitted. The above operators give rise to the following Feynman rule for the 4 point function:

$$\frac{i2M_V^2\gamma ^2}{v^2}\eta _{\mu \nu }+\frac{i4\gamma ^2}{\alpha v^2}(p_{1\mu }p_{2\nu }+p_{1\nu }p_{2\mu })$$
(11.10)
The presence of the terms proportional to $`\frac{1}{\alpha }`$ in (11.7) and (11.10) are in fact crucially important to obtain gauge invariant amplitudes, however for calculations in the unitary gauge $`\alpha \mathrm{}`$ their effect vanishes.

This is however not the complete story for the Feynman rules. The reason is that the radion couples conformally to the metric and so is not the same as the Higgs. Thus for example one finds that at one–loop the radion has the anomalous coupling

$$\frac{r}{\mathrm{\Lambda }}b_G\frac{\alpha _G}{8\pi }G_{\mu \nu }G^{\mu \nu }$$
(11.11)
in addition to the usual momentum–dependent coupling obtained from one-loop diagrams with internal fermions. Here $`b_G`$ is the beta–function. This may be understood as due to the scaling anomaly together with $`r`$ as a generator of scale transformations. Diagrammatically this result is obtained by preserving the conformal coupling of $`r`$ when the theory is regulated. For dimensional regularization this means that the radion must couple conformally to the $`D`$-dimensional metric. Since the linear coupling of the radion is obtained from varying the induced metric, for loop computations the radion should couple instead to $`\text{Tr}_DT_{\mu \nu }`$, but where now the trace is evaluated in $`D`$-dimensions. This differs from the above coupling by some operators whose coefficient vanishes when $`D4`$. Since this $`ϵ=4D`$ dependence can be offset by poles appearing in the loops, the appearance of these additional operators can result in finite non-zero results in the $`D4`$ limit. These operators will indeed have a non-vanishing contribution to the S and U parameter. Next we calculate the Feynman rules for these “anomalous” couplings. The interactions we should therefore study are

$$=\frac{h}{v}\text{Tr}_4T+\frac{\gamma r}{v}\text{Tr}_DT.$$
(11.12)
We point out that it is the original ‘gauge’ radion, $`r`$, that has the conformal coupling to the metric, and consequently it is this field which appears in the above interactions. But a result of the curvature–scalar term interaction is to introduce mixing between the radion and Higgs, which implies that after transforming to the mass basis the physical Higgs, $`h_m`$, will have couplings similar to those above.

The radion will thus have the interaction terms $`(D=4ϵ)`$

$$=\gamma \frac{r}{v}\text{Tr}_DT_{\mu \nu }=\gamma \frac{r}{v}\text{Tr}_4T\gamma \frac{ϵ}{4}\frac{r}{v}F_{\mu \nu }F^{\mu \nu }+\gamma \frac{ϵ}{2}\frac{r}{v}(2M_W^2W^+W^{}+M_Z^2Z^2),$$
(11.13)
where the first term is the one we have already discussed above. We will later add in the Higgs and mixing coefficients. The terms relevant to the gauge boson propagators from the last two terms are:

$`_{int}^{anom}=`$ $`{\displaystyle \frac{ϵ\gamma }{2v}}r[(_\mu Z_\nu )^2(_\mu Z_\nu )(_\nu Z_\mu )+2(_\mu W_\nu ^+)(^\mu W^\nu )2(_\mu W_\nu ^+)(^\nu W^\mu )`$ (11.14)
$`2M_W^2W_\mu ^+W^\nu M_Z^2Z_\mu Z^\nu ].`$
In $`D`$ dimensions, one also has to modify the gauge fixing terms. The covariant derivative of a vector field will be modified to

$$D_\mu V^\mu =\mathrm{\Omega }^2(_\mu V^\mu \frac{(D2)\gamma }{v}_\mu rV^\mu ).$$
(11.15)
In addition, the $`\sqrt{g}`$ factor in front of the gauge fixing terms have to be modified to $`\mathrm{\Omega }^D`$. Thus these interaction terms modify the Feynman rules for the interaction vertex given in (11.7) to

$$\frac{2iM_V^2\gamma }{v}(1\frac{ϵ}{2})\eta _{\mu \nu }+\frac{i\gamma ϵ}{v}(p_2p_3\eta _{\mu \nu }p_{2\nu }p_{3\mu })\frac{i(2ϵ)\gamma }{\alpha v}(p_{2\mu }p_{1\nu }+p_{3\nu }p_{1\mu })\frac{i\gamma ϵ}{\alpha v}p_{2\mu }p_{3\nu }.$$
(11.16)
In addition, the four-point vertex in (11.10) is also modified by terms proportional to $`ϵ`$, which however do not contribute to a calculation in the unitary gauge.

We close this section by discussing how to take the mixing between the Higgs and the radion due to the possible presence of the curvature-scalar mixing operator into account. The interaction in the gauge basis is

$$=\frac{h}{v}\text{Tr}_4T+\frac{\gamma r}{v}\text{Tr}_DT.$$
(11.17)
Using

$`h`$ $`=`$ $`ah_m+br_m\text{ ,}`$ (11.18)
$`r`$ $`=`$ $`ch_m+dr_m,`$ (11.19)
where the coefficients $`a`$, $`b`$, $`c`$ and $`d`$ can be read off from (10.17),(10.18), then

$``$ $`=`$ $`\left((a\gamma c){\displaystyle \frac{h_m}{v}}+(\gamma db){\displaystyle \frac{r_m}{v}}\right)\text{Tr}_4T`$ (11.20)
$`+(ch_m+dr_m){\displaystyle \frac{\gamma }{v}}\left({\displaystyle \frac{ϵ}{4}}F_{\mu \nu }F^{\mu \nu }+{\displaystyle \frac{ϵ}{2}}M_V^2V^2\right).`$
Thus the Feynman rules for the radion (and also for the Higgs) have to be modified such, that the above mixing terms are taken properly into account, for example for the mass eigenstate radion the Feynman rule will be

$$\frac{2iM_V^2}{v}(\gamma db\frac{\gamma dϵ}{2})\eta _{\mu \nu }+\frac{i\gamma dϵ}{v}(p_2p_3\eta _{\mu \nu }p_{2\nu }p_{3\mu })\frac{i(2ϵ)\gamma d}{\alpha v}(p_{2\mu }p_{1\nu }+p_{3\nu }p_{1\mu })\frac{i\gamma ϵd}{\alpha v}p_{2\mu }p_{3\nu }.$$
(11.21)
## 12 Electroweak Precision Measurements

In this Section we consider the corrections of the Randall–Sundrum model to the oblique parameters. Our analysis also applies more generally to a model with a light scalar coupled conformally to the SM metric but with a typical coupling of $`𝒪(`$TeV$`)^1`$.

Our approach is to use an effective theory with cutoff $`𝒪(\mathrm{\Lambda })`$, similar to the approach taken in . Below this scale the only light fields are the radion and Higgs whose contributions we are going to calculate explicitly. In our approach the effect of any modes heavier than the cutoff are included by introducing higher dimension operators that directly contribute to the oblique parameters. This in principle includes the effects of the heavy spin–2 KK states, for example, which are typically heavier than the radion. A direct computation of the effect of the heavy spin-2 states using a momentum–dependent regulator has been presented in .

In the previous Section the radion coupling to the gauge bosons was obtained and found to be similar to that of the Higgs. The contribution of the Higgs to the oblique parameters is by itself divergent, but these divergences are canceled by the contribution from the pseudo-Goldstone bosons (or the longitudinal states of the massive gauge bosons). Thus for the radion one expects a divergent contribution, but in contrast to the Higgs there is no additional source to cancel this. This is perhaps not surprising since the radion interactions are non-renormalizable.

A set of operators that provide the necessary counterterms for the wavefunction renormalization is

$$𝒪_X=\frac{g^2Z_X}{\mathrm{\Lambda }^2}\left(H^{}H\text{Tr}W_{\mu \nu }W^{\mu \nu }+\frac{1}{2}\mathrm{tan}^2\theta _WH^{}HB_{\mu \nu }B^{\mu \nu }+\mathrm{tan}\theta _WH^{}W_{\mu \nu }B^{\mu \nu }H\right),$$
(12.1)
where $`W_{\mu \nu }=W^a\tau ^a`$, with the generators normalized to $`1/2`$. Note that the last operator is gauge invariant, since for a gauge transformation $`U`$, $`W_{\mu \nu }UW_{\mu \nu }U^{}`$. Setting the Higgs to its vev in the above operator gives

$$𝒪_X\frac{g^2Z_Xv^2}{2\mathrm{\Lambda }^2}\left(W_{\mu \nu }^+W_{}^{}{}_{}{}^{\mu \nu }+\frac{1}{2\mathrm{cos}^2\theta _W}Z_{\mu \nu }Z^{\mu \nu }\right).$$
(12.2)
In this model the radion does not contribute to $`\gamma \gamma `$ and $`\gamma Z`$ wavefunction renormalization, and the absence of these counterterms uniquely fixes the coefficients in (12.1). The explicit computation of the $`ZZ`$ and $`WW`$ wavefunction renormalizations demonstrates that the above operator has the correct relative factor between the two gauge bosons. This is perhaps non–trivial, since there is no additional degree of freedom to fix this relative factor. We also note that an identical operator to (12.2) is also required in technicolor theories in order to cancel the divergent contribution of pseudo–Goldstone bosons to the oblique parameters .

The operator which provides the counterterms for the mass renormalization is

$$𝒪_M=\frac{Z_M}{2\mathrm{\Lambda }^2}\left(g_{}^{}{}_{}{}^{2}(D_\mu H^{}H)(H^{}D^\mu H)+g^2H^{}H(D_\mu H^{}D^\mu H)\right).$$
(12.3)
The first operator that appears here violates the custodial symmetry, and in particular contributes to the $`Z`$ but not the $`W`$ mass. After electroweak symmetry breaking they together reduce to

$$𝒪_M\frac{Z_M}{\mathrm{\Lambda }^2}\left(M_W^4W_\mu ^+W^\mu +\frac{M_Z^4}{2}Z^\mu Z_\mu \right).$$
(12.4)
Note that it is $`m_V^4`$ which appears, so we explicitly see that this operator contributes to the $`\rho `$ parameter. In this case it is trivial to obtain the correct mass renormalization from (12.3), since here there are two coefficients to be determined from only two constraints. So in addition to the usual Standard Model renormalizations, these wavefunction and mass counterterms are also required to renormalize the model.

The model with the radion represents an effective theory valid for $`E\begin{array}{c}<\hfill \\ \hfill \end{array}\mathrm{\Lambda }`$. The dimension–6 operators discussed above are obtained by integrating some unknown degrees of freedom in the full theory, and in the effective theory they appear with some unknown coefficients $`Z_i(\mathrm{\Lambda })`$. These for example could include the effects of integrating out the heavy spin-2 KK modes. Since the divergences for which the two above operators act as counterterms arise at one loop order, they are proportional to $`\gamma ^2/(16\pi ^2)`$. Thus it is reasonable to expect that the finite part of the operator is also of the same order, and in order to match the form of the explicitly calculated one-loop corrections we will write the (finite) coefficients as

$$Z_M=\frac{\gamma ^2\mathrm{\Lambda }^2}{16\pi ^2v^2}a_M,Z_X=\frac{\gamma ^2\mathrm{\Lambda }^2}{16\pi ^2v^2}a_X,$$
(12.5)
where we expect that the dimensionless parameters $`a_M`$ and $`a_X`$ are at most of order one. Note that since $`\gamma ^2\mathrm{\Lambda }^2`$, in this parameterization the dimension–6 operators are still suppressed by $`\mathrm{\Lambda }^2`$.

These coefficients parameterize the unknown physics integrated out at the scale $`E\mathrm{\Lambda }`$. Since, however, the radion at one loop contributes to the anomalous dimension of these operators, when comparing to the experimental results evaluated at the $`Z`$ mass large logarithms of $`O(\mathrm{ln}\mathrm{\Lambda }/M_Z)`$ appear from this anomalous scaling and this effect should be included. Following Wilson, a one-loop Wilsonian renormalization group equation is obtained for the operator coefficients. In the leading logarithm approximation the value of these coefficients at the weak scale is determined to be

$$Z_i(M_Z)=Z_i(\mathrm{\Lambda })+\frac{\beta _i}{16\pi ^2}\mathrm{ln}\frac{\mathrm{\Lambda }^2}{M_Z^2},$$
(12.6)
with the $`\beta _i`$ determined from the coefficient of the $`\mathrm{ln}\mu `$ term (or equivalently, from the $`1/ϵ`$ poles) in an explicit one–loop computation. To compute the oblique parameters, one then adds the contribution of these renormalized operators to the finite parts of the one-loop diagrams. In the leading logarithm approximation this amounts to simply replacing the $`1/ϵ`$ poles in the gauge–boson self-energies with $`\mathrm{ln}\mathrm{\Lambda }/M_Z`$.

To compute the oblique parameters one uses the Feynman rules in the previous Section to compute the two Feynman diagrams that contribute to the vacuum polarizations. Both diagrams have one internal radion, and one uses the three point function and the other uses the four point function. Our convention for the sign of the vacuum polarizations $`\mathrm{\Pi }_{VV}`$ is that
$$=i\mathrm{\Pi }_{VV}^{\mu \nu }(p^2)=i\eta ^{\mu \nu }\mathrm{\Pi }_{VV}(p^2)+ip^\mu p^\nu \stackrel{~}{\mathrm{\Pi }}_{VV}(p^2)$$
(12.7)

and only the first term is computed. The generic form of the radion contribution is

$$\mathrm{\Pi }_{VV}(p^2)=\mathrm{\Pi }_{VV}^S(p^2)+\mathrm{\Pi }_{VV}^A(p^2).$$
(12.8)
Here ’$`A`$’ denotes the anomalous contribution due to the conformal coupling of the radion, and ’$`S`$’ denotes the standard contribution which is also similar to the Higgs contribution (when $`\xi =0`$). The anomalous couplings are discussed in the previous Section. By an appropriate rescaling of the coupling we can also use $`\mathrm{\Pi }^S`$ for the Higgs. When $`\xi 0`$, the results given below can also be used to compute the oblique parameters after an appropriate redefinition of the couplings and masses. The modification to the expression for the oblique parameters is summarized in (12.21).

Inspecting the Feynman diagrams for the vacuum polarizations one finds that the quadratic divergences cancel between the two diagrams, leaving only a logarithmic divergence. Therefore this justifies the use of dimensional regularization. An explicit computation of the Feynman diagrams in unitary gauge and for vanishing curvature scalar parameter $`\xi `$ gives

$`\mathrm{\Pi }_{VV}^S(0)`$ $`=`$ $`{\displaystyle \frac{\gamma ^2}{16\pi ^2}}{\displaystyle \frac{m_V^4}{v^2}}\left({\displaystyle \frac{6}{ϵ}}+{\displaystyle \frac{5}{2}}{\displaystyle \frac{m_r^2}{2m_V^2}}+3{\displaystyle \frac{m_V^2\mathrm{ln}m_V^2/\mu ^2m_r^2\mathrm{ln}m_r^2/\mu ^2}{m_r^2m_V^2}}\right)\text{ ,}`$ (12.9)
$`\mathrm{\Pi }_{VV}^S(m_V^2)`$ $`=`$ $`{\displaystyle \frac{\gamma ^2}{16\pi ^2}}{\displaystyle \frac{m_V^4}{v^2}}\left({\displaystyle \frac{20}{3ϵ}}{\displaystyle \frac{2}{3}}{\displaystyle \frac{m_r^2}{m_V^2}}+{\displaystyle \frac{1}{3}}{\displaystyle \frac{m_r^4}{m_V^4}}+{\displaystyle \frac{10}{9}}\right)`$ (12.10)
$`+{\displaystyle \frac{\gamma ^2}{16\pi ^2}}{\displaystyle \frac{m_V^4}{v^2}}\left(4{\displaystyle \frac{4}{3}}{\displaystyle \frac{m_r^2}{m_V^2}}+{\displaystyle \frac{m_r^4}{m_V^4}}\right){\displaystyle _0^1}𝑑x\mathrm{ln}(x^2m_V^2+(1x)m_r^2)/\mu ^2`$
$`+{\displaystyle \frac{\gamma ^2}{16\pi ^2}}{\displaystyle \frac{m_V^4}{v^2}}\left({\displaystyle \frac{m_r^2}{m_V^2}}({\displaystyle \frac{m_r^2}{3m_V^2}}1)\mathrm{ln}{\displaystyle \frac{m_r^2}{\mu ^2}}+{\displaystyle \frac{1}{3}}({\displaystyle \frac{m_r^2}{m_V^2}}2)\mathrm{ln}{\displaystyle \frac{m_V^2}{\mu ^2}}\right)\text{ ,}`$
where to avoid confusion with defining too many $`\gamma ^{}s`$, the renormalization scale $`\mu `$ appearing includes the usual factors of $`4\pi `$ and Euler’s constant $`\gamma _E`$. An analytic expression for the Feynman parameter integral may be obtained, but it is not very illuminating. A powerful check on these expressions is gauge invariance. We have explicitly checked that in the general Rα gauge, the gauge parameter $`\alpha `$ cancels from the expression and reproduces the above results. We note that the divergences appearing here have the form as given by the operators in (12.1) and (12.3). For $`\mathrm{\Pi }(0)_{VV}m_V^4`$ as required by (12.3). The difference between $`\mathrm{\Pi }(m_V^2)`$ and $`\mathrm{\Pi }(0)`$ gives the divergence proportional to $`p^2`$, but in the above equations $`p^2`$ has already been set to $`m_V^2`$. With this in mind, by inspection the coefficient of $`p^2`$ is $`m_V^2`$, as required by (12.1). So the set of operators given by (12.1) and (12.3) provide an appropriate set of counterterms. We renormalize using $`\overline{\text{MS}}`$. Then using the Wilsonian approach outlined above, in the leading logarithm approximation we just replace $`1/ϵ`$ with $`\mathrm{ln}\mathrm{\Lambda }/M_Z`$, that is

$$\frac{2}{ϵ}\gamma _E+\mathrm{ln}4\pi \mathrm{ln}\frac{\mathrm{\Lambda }^2}{M_Z^2},$$
(12.11)
and set $`\mu =M_Z`$.

The anomalous contribution is finite and is

$$\mathrm{\Pi }_{VV}^A(p^2)=\frac{\gamma ^2}{16\pi ^2}\frac{m_V^4}{v^2}\left(\frac{10}{3}\frac{p^2}{m_V^2}+62\frac{m_r^2}{m_V^2}\right).$$
(12.12)
We note that the $`m_r^2`$ term does not contribute to the oblique parameters since it is $`p^2`$–independent and also $`m_V^2`$.

The PDG convention for the oblique parameters that we use here is

$`T`$ $`=`$ $`{\displaystyle \frac{1}{\alpha }}\left({\displaystyle \frac{\mathrm{\Pi }_{WW}(0)}{M_W^2}}{\displaystyle \frac{\mathrm{\Pi }_{ZZ}(0)}{M_Z^2}}\right)`$ (12.13)
$`S`$ $`=`$ $`{\displaystyle \frac{4\mathrm{sin}^2\theta _W\mathrm{cos}^2\theta _W}{\alpha }}\left({\displaystyle \frac{\mathrm{\Pi }_{ZZ}(M_Z^2)}{M_Z^2}}{\displaystyle \frac{\mathrm{\Pi }_{ZZ}(0)}{M_Z^2}}\right)`$ (12.14)
$`S+U`$ $`=`$ $`{\displaystyle \frac{4\mathrm{sin}^2\theta _W}{\alpha }}\left({\displaystyle \frac{\mathrm{\Pi }_{WW}(M_W^2)}{M_W^2}}{\displaystyle \frac{\mathrm{\Pi }_{WW}(0)}{M_W^2}}\right)`$ (12.15)
More generally one must also include the $`\mathrm{\Pi }_{Z\gamma }`$ and $`\mathrm{\Pi }_{\gamma \gamma }`$ self-energies. They have been dropped here since they do not receive contributions from either the Higgs or the radion.

Using the above expressions one can evaluate the contribution of the radion to $`S`$ and $`T`$ in the limit of a large radion mass. One obtains for $`\xi =0`$,

$`S`$ $`=`$ $`{\displaystyle \frac{\gamma ^2}{\pi }}\left(a_X{\displaystyle \frac{1}{12}}\mathrm{ln}{\displaystyle \frac{\mathrm{\Lambda }^2}{m_r^2}}{\displaystyle \frac{5}{72}}\right)`$ (12.16)
$`T`$ $`=`$ $`{\displaystyle \frac{3\gamma ^2}{16\pi \mathrm{cos}^2\theta _W}}\left({\displaystyle \frac{a_M}{3}}+\mathrm{ln}{\displaystyle \frac{\mathrm{\Lambda }^2}{m_r^2}}+{\displaystyle \frac{5}{6}}\right).`$ (12.17)
Inspecting the above expressions for the $`\mathrm{\Pi }`$’s one finds that there is no divergent contribution to $`U`$, and this is consistent with the fact that $`𝒪_X`$ and $`𝒪_M`$ provide no counterterms for $`U`$.

It is interesting that the radion contribution to $`S`$ is negative and to $`T`$ is positive. This is easy to understand by comparing this result to the contribution of the Higgs in the SM. In fact, for $`\gamma =1`$, $`m_r=m_h`$ and $`Z_i(\mathrm{\Lambda })=0`$, the radion result is identical to the (logarithmic) contribution of the SM physical Higgs. But there the total correction to $`S`$ and $`T`$ is finite, so that for the Higgs the $`\mathrm{\Lambda }`$ dependence in the above expression is canceled by loops of $`W`$’s and $`Z`$, leaving a $`\mathrm{ln}m_h/M_Z`$ dependence. Then one obtains the usual positive (negative) correction to $`S`$ $`(T)`$. For the radion though the large logarithms are present, and since $`\mathrm{\Lambda }>m_r`$, the radion contribution to $`S`$ is negative, and to $`T`$ it is positive. Recalling that $`\gamma ^2=v^2/6\mathrm{\Lambda }^2`$, the size of the above correction is only significant for small values of $`\mathrm{\Lambda }`$.

We conclude with a comment on the decoupling behavior of the radion. Inspecting the above expressions we see that for large $`m_r`$ the radion contribution scales as

$$\frac{1}{\mathrm{\Lambda }^2}\mathrm{ln}\frac{\mathrm{\Lambda }_C}{m_r}.$$
(12.18)
For the purposes of this paragraph we distinguish the cutoff of the effective theory, $`\mathrm{\Lambda }_C`$, from the mass scale $`\mathrm{\Lambda }=M_{Pl}e^{kr_0}`$ appearing in the radion coupling. In our analysis we have approximated $`\mathrm{\Lambda }_C\mathrm{\Lambda }`$. If the radion mass is much larger then $`\mathrm{\Lambda }`$, then the cutoff of the effective theory is the much higher scale $`Em_r`$ or larger. Then it is more natural to express all higher dimension operators as suppressed by $`m_r^1`$ or $`\mathrm{\Lambda }_C^1`$. But then the coupling of the radion to the gauge bosons contains the very large coefficient $`m_r/\mathrm{\Lambda }`$. It is then inappropriate to use the one–loop approximation. So to remain within the validity of the approximations used here one must also increase $`\mathrm{\Lambda }`$ for large $`m_r`$. Then the radion contribution decouples.

For large radion mass it is conceivable that the coupling of the radion to Tr$`T`$ decreases. This cannot be seen in our computations here because we have made the approximation of using only the zero mode wavefunction to determine the coupling, but have included the backreaction perturbatively to compute the radion mass. For very large radion mass these approximations are invalid, and then one must exactly solve the equations. It is then conceivable that for large radion mass the radion wavefunction on the TeV brane decreases in such a manner that the coupling to Tr$`T`$ remains natural, i.e. $`O(m_r^1)`$.

Following standard practice we define a reference model in which one computes the oblique parameters within the Standard Model, which means for some specific value for the Higgs mass. Since the curvature scalar operator mixes the radion and Higgs, the two physical scalars are some mixture of the gauge Higgs and radion. This mixture is not a unitary rotation due to the kinetic mixing between the states. The couplings of the ‘Higgs’ in this case is somewhat different than in the Standard Model, and for the purposes of computing the oblique parameters it is easiest to think of this as a new model, rather than as a perturbation to the Standard Model. So in computing the oblique parameters the Standard Model Higgs contribution (for the reference Higgs mass) should be subtracted out, and the contribution of the physical states in this model added back in. That is,

$$XX_{SM}^{ref}=X_{new}(m_h,m_r,\mathrm{\Lambda },\xi )X_H^{ref}(m_h=m_h^{ref}).$$
(12.19)
As mentioned, $`X_H^{ref}(m_h=m_h^{ref})`$ is the contribution from only the Higgs, with mass $`m_h^{ref}`$ and with Standard Model couplings, and it is independent of $`\mathrm{\Lambda }`$, $`m_r`$ and the curvature-mixing parameter $`\xi `$. The quantity $`X_{SM}^{ref}`$ is the full SM contribution, with the Higgs set at the reference mass $`m_h^{ref}`$. The new physics contribution contains two pieces,

$$X_{new}=X_R+X_H\text{ ,}$$
(12.20)
which describe the contribution of the physical radion $`(X_R)`$, and the physical Higgs $`(X_H)`$. In the limit of no curvature–scalar mixing, this last contribution is just that of a Standard Model Higgs with mass $`m_h`$. For a general curvature scalar mixing one just needs to include the effect of the mixing coefficients $`a`$, $`b`$, $`c`$ and $`d`$ given in the previous section. Then

$`X_{new}`$ $`=`$ $`\left(\mathrm{cos}\theta {\displaystyle \frac{\gamma }{Z}}(6\xi 1)\mathrm{sin}\theta \right)^2X(m_h^{phys},\gamma =1)`$ (12.21)
$`+\left(\mathrm{sin}\theta +{\displaystyle \frac{\gamma }{Z}}(6\xi 1)\mathrm{cos}\theta \right)^2X(m_r^{phys},\gamma =1)`$
$`{\displaystyle \frac{\gamma ^2}{Z^2}}(6\xi 1)X^A.`$
Here $`X^A`$ is the anomalous contribution, obtained from the vacuum polarization (12.12), but dropping the $`m_r^2`$ term since this does not contribute to any of the oblique parameters. The other $`X`$’s are obtained from using the vacuum polarizations (12.9), (12.10), and inserting the physical mass of the appropriate state.

We note that, for example, the full anomalous contribution to $`S`$ from the above formula is

$$S_{new}^A=\frac{5\gamma ^2}{6Z^2\pi }(6\xi 1)$$
(12.22)
and is negligible for reasonable values of $`\gamma `$ and $`\xi `$. The anomalous contribution to $`T`$ is even smaller.

### 12.1 Numerical Results

The “new” contribution, $`X`$, is constrained to lie within the measured values (extracted assuming $`m_h^{\mathrm{SM}\mathrm{ref}}=100`$ GeV)

$`S_{\mathrm{meas}}`$ $`=`$ $`0.07\pm 0.11`$ (12.23)
$`T_{\mathrm{meas}}`$ $`=`$ $`0.10\pm 0.14`$ (12.24)
$`U_{\mathrm{meas}}`$ $`=`$ $`0.11\pm 0.15.`$ (12.25)
(The errors are for one sigma.) $`X`$ can be easily calculated once $`m_h`$, $`m_r`$, $`\mathrm{\Lambda }`$, and $`\xi `$ are specified. Notice that since the current best-fit values for the electroweak parameters are nonzero, the new contributions can be more weakly or more strongly constrained depending on whether they add destructively or constructively with the Higgs, respectively. As a first example, we show the contribution to $`S`$ and $`T`$ in Fig. 3 as a function of $`m_h=m_r`$ (the “gauge” masses), fixing $`\mathrm{\Lambda }=1`$ TeV. Each contour corresponds to a different value of $`\xi `$, and the contours end when a physical mass exceeds the cutoff. The unshaded region corresponds to the $`1\sigma `$ allowed region. Notice that $`T`$ is a strong constraint on small (gauge) masses, while $`S`$ is a strong constraint for large masses. Also, the the case with $`\xi =0`$ is nearly identical to the ordinary SM Higgs contribution, since the radion contribution that can be separated out is strongly suppressed by the coupling $`\gamma ^2`$. In Fig. 4 we show the contribution to $`S`$ and $`T`$ as a function of $`m_r`$ with $`m_h=300`$ GeV. Notice that the contributions are nearly independent of $`m_r`$ for small curvature scalar mixing.

The above results illustrate a general trend that with small or absent curvature-scalar mixing, the bound on the Higgs mass is not significantly affected by the presence of the radion. This is not true, however, if we allow $`\xi `$ to take larger values. It is easiest to first illustrate that radion physics with large curvature scalar mixing can significantly relax the bound on the Higgs mass, by scanning through the parameter space (choosing $`m_h=m_r`$ for simplicity) for values that satisfy one or two sigma limits on the electroweak parameters. We find sets of parameters that are not minor perturbations on the SM limit allow the physical masses of the Higgs and radion to be several hundred GeV, and perhaps even TeV scale. In Fig. 5 we show the range of physical masses and the range of $`\xi `$ as a function of the cutoff scale. In general there is not a unique mapping between the figures, however the “shark fin” structure for the one sigma region in Fig. 5(a) does correspond to the “inverse fin” in Fig. 5(b). Notice that at two sigma the physical Higgs mass can be much larger than the SM bound throughout the parameter space shown, and even at one sigma there exists a narrow range of large, negative curvature-scalar mixing where the physical Higgs mass could be of order a TeV. The latter result arises from a cancellation between the physical Higgs and radion contributions with the SM reference contribution. This can be seen in the limit of a large radion and Higgs mass. Since the dependence on the masses is only logarithmic, we can approximate the masses as being equal. Then, for example,

$$S_{\mathrm{new}}=\frac{1}{\pi }(\frac{1}{12}\mathrm{ln}\frac{m^2}{M_Z^2}\frac{5}{72})(6\xi 1)^2\frac{\gamma ^2}{Z^2\pi }(\frac{1}{12}\mathrm{ln}\frac{\mathrm{\Lambda }^2}{m^2}+\frac{5}{72}).$$
(12.26)
The first contribution is just the usual correction from the Higgs. But the second correction can be potentially large and negative due to the dependence on the $`\xi `$ parameter. It should be emphasized that the large correction is due to the (non–unitary) kinetic mixing between the radion and Higgs, or equivalently, due to the non–standard couplings of the radion and Higgs in the mass basis. Hence, while this region is provocative, it nonetheless requires fine-tuning.

These results have assumed that the contribution from the nonrenormalizable counterterms is small, meaning $`a_X`$ and $`a_M`$ are less than order $`1`$. For larger coefficients the allowed regions of parameter space, albeit only at moderately large $`\gamma `$. In Fig. 6 we show the shift in the contours, for the two sigma region, resulting from taking $`a_X=\pm 10`$. ($`a_M`$ was also taken to be $`10`$, but the effect on the contours was negligible.)

## 13 Limits on Radion Mass

As we found in Section 6, the mass of the radion is expected to be significantly below the the cutoff scale, placing it in a region that can be directly probed by experiments. The previous section has shown that the radion couples much like a Higgs boson, and in the limit $`\xi 0`$, the tree-level couplings of the radion are simply scaled by $`\gamma `$. Let us first consider the bounds in this case.

In the SM, the current bound on the Higgs mass comes primarily from the LEP processes $`e^+e^{}Z^{}Zh`$, with the value $`m_h^{\mathrm{SM}}<108`$ GeV . For the radion, an exactly analogous production process occurs $`e^+e^{}Z^{}Zr`$, except that the $`ZZr`$ coupling has a factor of $`\gamma `$. To a good approximation, we can therefore estimate the production cross section of radions at LEP by simply scaling the Higgs cross section by $`\gamma ^2`$.

The decay of the radion is somewhat more complicated, however. As we discussed in Section 11, the radion couples directly to gauge bosons through the conformal anomaly. Although this coupling is one-loop suppressed, it competes with Yukawa suppressed interactions and, for the case $`rgg`$, can be comparable or even dominate . In the radion mass range well below the $`t\overline{t}`$ threshold, the ratio of the two largest widths can be expressed as

$$\frac{\mathrm{\Gamma }(rgg)}{\mathrm{\Gamma }(rb\overline{b})}=\frac{\alpha _s^2c_3^2}{12\pi ^2\beta ^3}\left(\frac{m_r}{m_b}\right)^2,$$
(13.1)
where $`\beta ^2=14m_b^2/m_r^2`$ and $`c_323/3`$ is roughly the one-loop QCD $`\beta `$-function coefficient (approximately including the smaller contribution from the one-loop triangle diagram with top quarks.) Notice that the coupling $`\gamma ^2`$ cancels in this ratio. This can be written in the suggestive form

$$\frac{\mathrm{\Gamma }(rgg)}{\mathrm{\Gamma }(rb\overline{b})}1/\beta ^3\left(\frac{m_r}{12m_b}\right)^2.$$
(13.2)
Hence, $`rgg`$ dominates for the region $`12m_b<m_r<2M_W`$. The search strategy for the radion is therefore significantly different from the Higgs in this mass window, namely searching for a pair of gluon jets instead of a pair of $`b`$-jets. Similarly, the radion has a different production cross section at hadron colliders via gluon fusion, proportional to the conformal anomaly-enhanced width into gluons but suppressed by the usual $`\gamma ^2`$ .

Determining an accurate bound on the radion mass in the region that can be probed by LEP requires a detailed analysis of detecting a two gluon plus Z signal. We will not attempt this here. Instead, the expected bound on the radion mass can be roughly estimated as a function of the coupling if we assume that some number of production events $`N`$ at LEP could not have escaped detection (or be lumped into SM backgrounds). Near the kinematical limit, the best bound will always come from the highest energy data. For lower mass radions, a lower center-of-mass energy results in a slightly higher cross section.***For instance, $`\sigma _{\sqrt{s}=189\mathrm{GeV}}(e^+e^{}rZ)/\sigma _{\sqrt{s}=202\mathrm{GeV}}(e^+e^{}rZ)1.3`$ for small $`m_r`$. Since the bound for a given radion mass is limited only by luminosity, we can combine the multitude of LEP runs at various energies by weighting by the integrated luminosity accumulated. The bound is then simply

$$\gamma ^2<\frac{N}{\sigma _\sqrt{s}(e^+e^{}Zh;m_h=m_r)\times _\sqrt{s}}$$
(13.3)
where the sum is over the various recent LEP runs with center-of-mass energy $`\sqrt{s}`$. $`N`$ encodes all of the detailed analyses of backgrounds, signal efficiencies, etc., and is in general not independent of energy or radion mass. In Fig. 7 we simply show the bound obtained if $`N=20`$ or $`100`$ (the integrated luminosity for each energy was also summed) corresponding to producing $`20`$ or $`100`$ events summed over all four LEP experiments. These numbers were chosen since searches for Higgs bosons typically need a few to tens of events (per detector) for a statistically significant signal-to-background ratio. Notice that no bound on the radion mass is expected from the recent LEP runs once $`\gamma `$ is less than about $`0.1`$.

One could also search for light radions, $`m_r<60`$ GeV, at LEP I via the decay $`Zf\overline{f}+r`$, through the same coupling discussed above. However, the expected bound obtained by using this procedure is no better than that found above for $`m_r`$ larger than about $`10`$ GeV.

We have not attempted to estimate a bound on $`\gamma `$ for a radion mass less than about $`10`$ GeV. We really do not expect the radion to be several orders of magnitude below the cutoff scale, and so at the outset it seems this mass region is unnatural. But, the presence of several low energy production processes (and rare decays) could be important, so a considerably more careful analysis than what we have attempted here is needed.

When curvature scalar mixing is included, the coupling $`ZZr`$ is modified as shown in Eq. (11.20). The above analysis can be translated into this more general case, but now the coupling is not simply $`\gamma `$ but a function of the curvature-scalar mixing as well. In addition, the SM Higgs couplings are also modified, and so its production and decay are also affected. In particular, the production cross section could be either enhanced or suppressed. (This is similar to what happens in two Higgs doublet models, such as the MSSM.) An interesting signal for RS with curvature-scalar mixing could be observing a nonstandard cross section or decay rate for a SM-like Higgs.

## 14 Conclusions

In this paper we have analyzed the coupled radion-scalar system in detail, including the backreaction of the bulk stabilizing scalar on the metric. We derived the coupled differential equations governing the dynamics of the system, and found the mass eigenvalues for some limiting cases. We find that due to the coupling between the radion and the bulk scalar, there will be a single KK tower describing the system, with the metric perturbations non-vanishing for every KK mode. This implies that the Standard Model fields localized on the TeV brane will couple to every KK mode from the bulk scalar, and this could provide a means to directly probe the stabilizing physics. We also found that in an expanding universe the shift in the radion at late times completely agrees with the effective theory result of .

We also calculated the contributions of the radion to the oblique parameters using an effective theory approach. Since the radion is the only new state well below the TeV scale, we argued that a low–energy effective theory including only the radion and SM fields is sufficient, as long as appropriate nonrenormalizable counterterms at the cutoff scale are added. In the absence of a curvature–scalar Higgs mixing operator, the size of the contribution to the oblique parameters due to the radion is small. In the presence of such a mixing operator, the corrections can be much larger due to the modified radion and Higgs couplings. In particular, including only the mixed radion and Higgs fields as “new physics”, we calculated the range of curvature-scalar mixing for a given cutoff scale that allows the physical Higgs mass to be up to of order the cutoff scale, while $`S_{\mathrm{new}}`$ and $`T_{\mathrm{new}}`$ were within the experimental limits. However, the parameters must be increasingly fine-tuned to achieve a Higgs mass that exceeds a few hundred GeV.

## Acknowledgements

After Sections 3-6 were completed, we were informed that many of the results of these sections are also contained in Ref. . We thank Riccardo Rattazzi for pointing us to this reference. We also thank Tanmoy Bhattacharya, Kiwoon Choi, Jim Cline, Gábor Cynolter, Josh Erlich, Christophe Grojean, Howie Haber, Tao Han, György Pócsik, Lisa Randall, John Terning and James Wells for useful discussions, Jim Cline for reading an earlier version of the manuscript, as well as Tao Han and John Terning for comments on the manuscript. C.C. is an Oppenheimer fellow at the Los Alamos National Laboratory. The research of C.C. is supported by the Department of Energy under contract W-7405-ENG-36. The work of M.G. was supported in part by the Department of Energy. The research of G.K. is supported in part by the U.S. Department of Energy under grant number DOE-ER-40682-143.