format: text

text:
# Cosmology of Randall-Sundrum models11footnote 1Talk given at “International Workshop on Particle Physics and the Early Universe”, COSMO2000, Cheju Island, Korea, September 2000

## 1 Introduction

For past two years, many particle physicists and cosmologists were excited by the development of two ideas, the brane world and the warped extra dimensions, both of which are based on the existence of extra dimensions. The basic idea of the brane world is that standard model particles are localized on a (3+1)-dimensional brane (or a set of branes) embedded in a higher dimensional spacetime, while gravity propagates in the whole bulk . The warped extra dimension assumes that the background metric is curved along the extra dimensions, mainly due to the negative bulk cosmological constant . Why are these two ideas so exciting? They have brought us fresh views and perspectives in gravity, cosmology, particle physics and string theory. We have seen many interesting issues discussed so far, such as the localization of gravity, the gauge hierarchy problem, the cosmological constant problem and self-tuning models , the construction of supersymmetric RS models and the role of supersymmetry , the connection to string theory or Horava-Witten model and warped compactification, the interpretation in light of AdS/CFT holographic duality , the collider signatures of the KK modes and the radion, etc .

In this talk, we focus on the cosmological aspect of the two ideas, mostly in five dimensional models with one extra dimension. There has been much interest in this because the five dimensional nature of gravity and the brane setup might lead to the non-conventional cosmology even at low temperature as well as at high temperature above TeV scale. It was found that the Friedmann equation of the brane shows a $`H\rho `$ behavior and has an additional dark radiation term . There was also a difficulty concerning the negative tension brane, and it was not very clear what happens at temperatures above TeV, which can alter the early universe cosmology including inflation.

However, in the early models, a few important ingredients are not addressed, such as the mechanism for the localization of fields on the brane, the stabilization of the brane systems  and the way to achieve necessary fine tunings of parameters. They inevitably involve bulk matter and dynamics, and can change the whole picture, for example, changing the exponential warp factor to power law. For the cosmological consequences of brane world models, taking the stabilization into account is found to be crucial . The inclusion of the effect of stabilization recovers the conventional FRW cosmology at temperatures below TeV.

The talk is organized as follows. First, we derive the effective four-dimensional brane equations to see the generalities of brane dynamics, though its usefulness is limited by the lack of the knowledge of bulk effects. Then, we try to solve the five dimensional equations with an appropriate ansatz. The solution can be obtained in very restricted cases, but a framework can be found where we can study brane cosmology with the effect of stabilization taken into account. This is done through the $`\widehat{T}_5^5`$ component which is adjusted to stabilize the extra dimension in the presence of brane matter. Based on this, we analyze the background spacetime where we also discuss the mass scales and the hierarchy problem, and one and two brane models in turn.

## 2 Effective four-dimensional equations on the brane

### 2.1 Framework

We consider the five-dimensional spacetime with coordinates $`(\tau ,x^i,y)`$, and 3-branes embedded in it. The action describing our framework is

$`S`$ $`=`$ $`{\displaystyle _M}d^5x\sqrt{\widehat{g}}\left[{\displaystyle \frac{M^3}{2}}\widehat{R}\mathrm{\Lambda }_b+_{bM}\right]`$ (1)
$`+{\displaystyle \underset{i}{}}{\displaystyle d^4x\sqrt{g^{(i)}}\left[\mathrm{\Lambda }_i+_{iM}\right]}`$
where $`M`$ is the fundamental gravitational scale of the model, $`\mathrm{\Lambda }_b`$ and $`\mathrm{\Lambda }_i`$ represent the bulk cosmological constant and brane tensions, respectively. $`_{bM}`$ and $`_{iM}`$ are Lagrangians for the bulk fields and for the fields localized in the branes.

To investigate the role of bulk and bulk fields in brane dynamics  we introduce a bulk scalar $`\widehat{\mathrm{\Phi }}`$, with $`_{bM}=\frac{1}{2}(\widehat{}\widehat{\mathrm{\Phi }})^2`$. We also allow that $`\mathrm{\Lambda }_b`$, $`\mathrm{\Lambda }_i`$ and $`_i`$ can be functions of $`\widehat{\mathrm{\Phi }}`$. The bulk Einstein equations obtained from the action (1) are

$$\widehat{G}_{MN}=\widehat{T}_{MN}=\widehat{}_M\widehat{\mathrm{\Phi }}\widehat{}_N\widehat{\mathrm{\Phi }}\widehat{g}_{MN}\left[\frac{1}{2}(\widehat{}\widehat{\mathrm{\Phi }})^2+\mathrm{\Lambda }_b(\widehat{\mathrm{\Phi }})\right],$$
(2)
and the bulk scalar equation is

$$\widehat{\mathrm{}}\widehat{\mathrm{\Phi }}\frac{d\mathrm{\Lambda }_b(\widehat{\mathrm{\Phi }})}{d\widehat{\mathrm{\Phi }}}=0.$$
(3)
The existence of branes imposes the junction conditions on the above bulk equation, giving discontinuities in the first derivatives across the branes of metric and scalar field. First let us proceed the general setup. A brane can be described by the normal vector $`n_M`$. Then the induced metric on the brane is given by $`g_{MN}=\widehat{g}_{MN}n_Mn_N`$, while the extrinsic curvature by $`K_{MN}=_Mn_N`$. The junction conditions are

$`[K_{\mu \nu }]`$ $`=`$ $`{\displaystyle \frac{1}{2M^3}}\left(T_{\mu \nu }{\displaystyle \frac{1}{3}}g_{\mu \nu }T\right),`$ (4)
$`[n_M^M\widehat{\mathrm{\Phi }}]`$ $`=`$ $`{\displaystyle \frac{1}{2M^3}}\left({\displaystyle \frac{d\mathrm{\Lambda }_i}{d\widehat{\mathrm{\Phi }}}}+{\displaystyle \frac{\delta _{iM}}{\delta \widehat{\mathrm{\Phi }}}}\right),`$ (5)
where the bracket in the left hand side means the difference across the brane and $`T_{\mu \nu }`$ is the energy momentum tensor of brane matter, i.e., $`\mathrm{\Lambda }_i`$ and $`_{iM}`$.

Suppose that we are localized on a certain brane and want to study the brane dynamics as observed by us. We may take two different approaches. The first approach is to derive the effective 3-brane equations localized on our brane. The second is to directly solve the whole bulk equations. In this section, we follow the first approach. The second will be dealt with in the next section.

### 2.2 The effective 3-brane Einstein equation

To derive the effective 3-brane equations, we need to know the extrinsic curvature and the intrinsic curvature of our brane in terms of the bulk metric $`\widehat{g}_{MN}`$ and the normal vector $`n_M`$, which are provided by the Codacci equation

$$_MK_\mu ^M_\mu K=g_\mu ^MG_{MN}n^N,$$
(6)
and the Gauss equation

$`G_{\mu \nu }`$ $`=`$ $`{\displaystyle \frac{2}{3}}\left[G_{MN}g_\mu ^Mg_\nu ^N+\left(G_{MN}n^Mn^N{\displaystyle \frac{1}{4}}G\right)g_{\mu \nu }\right]`$ (7)
$`+KK_{\mu \nu }K_\mu ^MK_{M\nu }{\displaystyle \frac{1}{2}}\left(K^2K_{MN}K^{MN}\right)g_{\mu \nu }E_{\mu \nu },`$
where $`E_{\mu \nu }=C_{MNOP}n^Mn^Og_\mu ^Ng_\nu ^P`$ and $`C_{MNOP}`$ is the bulk Weyl tensor.

Now we take $`y`$ as the Gaussian normal coordinates and impose $`Z_2`$ symmetry, $`yy`$. We expand $`\widehat{\mathrm{\Phi }}`$ around the brane,

$$\widehat{\mathrm{\Phi }}(x,y)=\varphi (x)+\mathrm{\Phi }_1(x)|y|+\frac{1}{2}\mathrm{\Phi }_2(x)y^2+𝒪(y^3).$$
(8)
Then from the bulk equations (2) and (3), and the junction conditions (4) and (5), we obtain the equation for $`\varphi `$

$$\mathrm{}\varphi \frac{d\mathrm{\Lambda }_b}{d\varphi }=\frac{T}{12M^3}\left(\frac{d\mathrm{\Lambda }}{d\varphi }\right)\mathrm{\Phi }_2,$$
(9)
and the four-dimensional effective Einstein equation for the brane

$$G_{\mu \nu }=\mathrm{\Lambda }_{\mathrm{eff}}(\varphi )g_{\mu \nu }+\frac{\mathrm{\Lambda }(\varphi )}{6M^3}T_{\mu \nu }+\frac{1}{M^6}\mathrm{\Pi }_{\mu \nu }E_{\mu \nu }+\frac{2}{3M^3}\widehat{T}_{\mu \nu }(\varphi ),$$
(10)
where

$`\mathrm{\Lambda }_{\mathrm{eff}}(\varphi )`$ $`=`$ $`{\displaystyle \frac{1}{2M^3}}\left[\mathrm{\Lambda }_b+{\displaystyle \frac{\mathrm{\Lambda }^2}{6M^3}}{\displaystyle \frac{1}{8}}\left({\displaystyle \frac{d\mathrm{\Lambda }}{d\varphi }}\right)^2\right],`$ (11)
$`\mathrm{\Pi }_{\mu \nu }`$ $`=`$ $`{\displaystyle \frac{1}{4}}T_{\mu \lambda }T_\nu ^\lambda +{\displaystyle \frac{1}{12}}TT_{\mu \nu }+{\displaystyle \frac{1}{8}}g_{\mu \nu }T_{\lambda \rho }T^{\lambda \rho }{\displaystyle \frac{1}{24}}g_{\mu \nu }T^2,`$ (12)
$`\widehat{T}_{\mu \nu }(\varphi )`$ $`=`$ $`_\mu \varphi _\nu \varphi {\displaystyle \frac{5}{8}}g_{\mu \nu }(\varphi )^2.`$ (13)
The first three terms in the right hand side of (10) deal with the sources on the brane. The first two terms are same as four-dimensional Einstein gravity, if we identify the four-dimensional Planck mass

$$M_P^2=\frac{\mathrm{\Lambda }(\varphi )}{6M^6}.$$
(14)
The third term gives a correction quadratic in energy-momentum tensor. The last two terms in Eq. (10) are bulk effect terms, which reflect the existence of bulk in brane dynamics. They are inputs from bulk dynamics and not determined in the brane equations. The brane equations (9) and (10) are not closed equations. In this sense, they would not be very useful without the knowledge of the bulk effect terms. However, they reveal a most general structure of the brane equations. In this regard, we note that the Planck mass (14) seems at first to be determined solely by the brane tension. This seems strange because the graviton (and its zero mode) comes from the bulk fields. Therefore there must be some correlation between the brane tension and bulk dynamics, which are incorporated in the brane equations through these bulk effect terms. We will see this in the following sections where the bulk solutions are treated.

## 3 Five-dimensional cosmological solutions

### 3.1 Framework

In this section, we investigate the cosmological bulk solutions of five-dimensional models. For the two brane models, the fifth dimension $`y`$ is assumed to be an orbifold $`S^1/Z_2`$ with $`yy+1`$ and $`yy`$ identified. Two 3-branes reside at two fixed points (boundaries) $`y=0,\frac{1}{2}`$.

Since we are interested in the cosmological solution, we consider the metric where the 3-dimensional spatial section is homogeneous and isotropic.

$$ds^2=n^2(\tau ,y)d\tau ^2+a^2(\tau ,y)\gamma _{ij}dx^idx^j+b^2(\tau ,y)dy^2,$$
(15)
where $`\gamma _{ij}`$ is the 3-dimensional homogeneous and isotropic metric, and we will use $`K=1,0,+1`$ to represent its spatial curvature. Einstein equations are given by $`\widehat{G}_{MN}=(1/M^3)\widehat{T}_{MN}`$ where

$`\widehat{G}_{00}`$ $`=`$ $`3\left\{{\displaystyle \frac{\dot{a}}{a}}\left({\displaystyle \frac{\dot{a}}{a}}+{\displaystyle \frac{\dot{b}}{b}}\right){\displaystyle \frac{n^2}{b^2}}\left[{\displaystyle \frac{a^{\prime \prime }}{a}}+{\displaystyle \frac{a^{}}{a}}\left({\displaystyle \frac{a^{}}{a}}{\displaystyle \frac{b^{}}{b}}\right)\right]+K{\displaystyle \frac{n^2}{a^2}}\right\}`$ (16)
$`\widehat{G}_{ij}`$ $`=`$ $`{\displaystyle \frac{a^2}{n^2}}\gamma _{ij}\left\{2{\displaystyle \frac{\ddot{a}}{a}}+{\displaystyle \frac{\dot{a}}{a}}\left({\displaystyle \frac{\dot{a}}{a}}2{\displaystyle \frac{\dot{n}}{n}}\right)+{\displaystyle \frac{\ddot{b}}{b}}+{\displaystyle \frac{\dot{b}}{b}}\left(2{\displaystyle \frac{\dot{a}}{a}}{\displaystyle \frac{\dot{n}}{n}}\right)\right\}K\gamma _{ij}`$ (17)
$`+{\displaystyle \frac{a^2}{b^2}}\gamma _{ij}\left\{2{\displaystyle \frac{a^{\prime \prime }}{a}}+{\displaystyle \frac{n^{\prime \prime }}{n}}+{\displaystyle \frac{a^{}}{a}}\left({\displaystyle \frac{a^{}}{a}}+2{\displaystyle \frac{n^{}}{n}}\right){\displaystyle \frac{b^{}}{b}}\left(2{\displaystyle \frac{a^{}}{a}}+{\displaystyle \frac{n^{}}{n}}\right)\right\}`$
$`\widehat{G}_{55}`$ $`=`$ $`3\left\{{\displaystyle \frac{b^2}{n^2}}\left[{\displaystyle \frac{\dot{a}}{a}}\left({\displaystyle \frac{\dot{a}}{a}}{\displaystyle \frac{\dot{n}}{n}}\right)+{\displaystyle \frac{\ddot{a}}{a}}\right]+{\displaystyle \frac{a^{}}{a}}\left({\displaystyle \frac{a^{}}{a}}+{\displaystyle \frac{n^{}}{n}}\right)\right\}`$ (18)
$`\widehat{G}_{05}`$ $`=`$ $`3\left({\displaystyle \frac{n^{}}{n}}{\displaystyle \frac{\dot{a}}{a}}+{\displaystyle \frac{a^{}}{a}}{\displaystyle \frac{\dot{b}}{b}}{\displaystyle \frac{\dot{a}^{}}{a}}\right)`$ (19)
$`\widehat{T}_N^M`$ $`=`$ $`\mathrm{diag}[\widehat{\rho },\widehat{p},\widehat{p},\widehat{p},\widehat{p}_5]+{\displaystyle \underset{i=0,\frac{1}{2}}{}}{\displaystyle \frac{\delta (y_i)}{b}}\mathrm{diag}[\rho _i,p_i,p_i,p_i,0]`$ (20)
In addition to Einstein equations, we use the energy-momentum conservation equation, $`_M\widehat{T}_N^M=0`$

$`{\displaystyle \frac{d\widehat{\rho }}{d\tau }}+3(\widehat{\rho }+\widehat{p}){\displaystyle \frac{\dot{a}}{a}}+(\widehat{\rho }+\widehat{p}_5){\displaystyle \frac{\dot{b}}{b}}=0,`$ (21)
$`\widehat{p}_5^{}+\widehat{p}_5\left({\displaystyle \frac{n^{}}{n}}+3{\displaystyle \frac{a^{}}{a}}\right)+{\displaystyle \frac{n^{}}{n}}\widehat{\rho }3{\displaystyle \frac{a^{}}{a}}\widehat{\rho }=0.`$ (22)
Brane sources in Eq. (20) can be converted to boundary conditions (Junction conditions): $`n`$, $`a`$ and $`b`$ must be continuous and $`n^{}`$, $`a^{}`$ are discontinuous due to boundary sources by

$$\frac{1}{b}\frac{n^{}}{n}|_{y_i}^{y_i+}=\frac{2\rho _i+3p_i}{3M^3},\frac{1}{b}\frac{a^{}}{a}|_{y_i}^{y_i+}=\frac{\rho _i}{3M^3}.$$
(23)
### 3.2 Five-dimensional spacetime with the bulk cosmological constant and brane tensions

First, we consider the five-dimensional spacetime supported by the negative bulk cosmological constant $`\widehat{\rho }=\widehat{p}=\mathrm{\Lambda }_b<0`$ and brane tensions $`\rho _i=p_i=\mathrm{\Lambda }_i`$. We define the parameters

$$k=\left(\frac{\mathrm{\Lambda }_b}{6M^3}\right)^{1/2},k_i=\frac{\mathrm{\Lambda }_i}{6M^3}$$
(24)
In general cases with $`k_0,k_{1/2}k`$, the metric is given by

$$ds^2=\frac{d\tau ^2+\delta _{ij}dx^idx^j+(k\tau b_0)^2dy^2}{\left[k\tau \mathrm{sinh}(kb_0|y|+c_0)+g_0\right]^2}$$
(25)
where $`b_0`$ and $`c_0`$ are determined by

$$k\mathrm{cosh}(c_0)=k_0,k\mathrm{cosh}(\frac{1}{2}kb_0+c_0)=k_{\frac{1}{2}}.$$
(26)
The metric (25) describes a slice of AdS5 space inflating both in extra dimension and in spatial dimensions. We note two special cases. If we have a fine tuning $`(k_2k)/(k_1k)=e^{kb_o}(k_2+k)/(k_1+k)`$, we obtain a inflationary solution with static extra dimension

$$ds^2=\frac{d\tau ^2+e^{2k\tau }\delta _{ij}dx^idx^j+b_0^2dy^2}{\mathrm{sinh}^2(kb_0|y|+c_0)}$$
(27)
With two fine tunings $`k=k_0=k_{\frac{1}{2}}`$, we can get a static solution, the Randall-Sundrum model

$$ds^2=e^{2kb_0|y|}\eta _{\mu \nu }dx^\mu dx^\nu +b_0^2dy^2$$
(28)
This model attracted much attention because it convert the gauge hierarchy problem into a geometric problem. The four-dimensional Planck scale in this model is given by $`M_P^2=(M^3/k)[1e^{kb_0}]`$. Any mass parameter $`m_0`$ on the visible brane corresponds to a physical mass $`m=m_0e^{\frac{1}{2}kb_0}`$, which we identify as the weak scale. Hence the large hierarchy between the Planck scale and the weak scale $`M_P/M_W10^{16}`$ can be explained by a warp factor $`e^{\frac{1}{2}kb_0}`$ with $`\frac{1}{2}kb_037`$. The property is not spoiled by $`R^2`$ corrections, if $`R^2`$ corrections are given by Gauss-Bonnet interactions .

Note that the space described by the metric (25) is locally AdS5, because it shares the same bulk equation. The metric (25) can be transformed to (28) by a coordinate transformation

$$e^{kb_{\mathrm{RS}}y_{\mathrm{RS}}}=k\tau \mathrm{sinh}(kb_0y+c_0)+g_0,\tau _{\mathrm{RS}}=\tau \mathrm{cosh}(kb_0y+c_0).$$
(29)
Therefore, we can say that above metrics are the slices of AdS5 with different boundary geometries.

### 3.3 Cosmological solutions with static extra dimension

Let us turn to the more interesting situation where the matter is added on the brane(s) and possibly in the bulk. But we require the extra dimension to be static, that is $`\dot{b}=0`$. This requires a fine tuning between matter densities. For the bulk energy-momentum, we assume

$$\widehat{\rho }=\mathrm{\Lambda }_b,\widehat{p}=\mathrm{\Lambda }_b,\widehat{p}_5=\mathrm{\Lambda }_b+p_5(\tau ,y).$$
(30)
where $`\mathrm{\Lambda }_b<0`$. The form of $`p_5(\tau ,y)`$ is constrained by Eq. (22) to be

$$p_5(\tau ,y)=\frac{\stackrel{~}{p}_5(\tau )}{n(\tau ,y)a^3(\tau ,y)}$$
(31)
With the gauge fixing $`n(\tau ,y=0)=1`$, we obtain the solution

$`a^2(\tau ,y)`$ $`=`$ $`a_0^2[(1+{\displaystyle \frac{\dot{a}_0^2+K}{2k^2a_0^2}})\mathrm{cosh}(2kby){\displaystyle \frac{\dot{a}_0^2+K}{2k^2a_0^2}}`$ (32)
$`\pm (1+{\displaystyle \frac{1}{k^2}}[{\displaystyle \frac{\dot{a}_0^2+K}{a_0^2}}+{\displaystyle \frac{C}{a_0^4}}])^{1/2}\mathrm{sinh}(2kby)],`$
and $`n(\tau ,y)=\dot{a}(\tau ,y)/\dot{a}_0(\tau )`$. Here $`a_0(\tau )a(\tau ,y=0)`$ and $`C(\tau )`$ is determined by $`\stackrel{~}{p}_5(\tau )`$ up to a constant through the equation

$$\dot{C}=\frac{2\dot{a}_0(\tau )}{3M^3}\stackrel{~}{p}_5(\tau ).$$
(33)
$`a_0(\tau )`$ is fixed by the boundary condition. Suppose that a brane with the energy density $`\rho _0`$ is placed at $`y=0`$, and assume $`Z_2`$ symmetry $`yy`$. The junction condition (23) at the brane gives the evolution equation for $`a_0(\tau )`$

$$\left(\frac{\dot{a}_0}{a_0}\right)^2+\frac{K}{a_0^2}=k^2+\left(\frac{\rho _0}{6M^3}\right)^2\frac{C}{a_0^4}$$
(34)
If we further assume that we are living on the brane, this equation is nothing but the Friedmann equation of our universe. However, in the simple case where $`\mathrm{\Lambda }_b=0`$ and $`p_5=0`$, this equation differs from the four-dimensional Friedmann equation in two aspects. First, the Hubble parameter $`H=\dot{a}_0/a_0`$ is proportional to $`\rho _0`$ instead of $`\sqrt{\rho _0}`$. Second, there is a $`C/a_0^4`$ term which looks like a radiation term. This means that we can have dynamic solution without any matter on the brane or in the bulk. This term seems to have an interesting interpretation in view of AdS/CFT . Both of these alter the big bang nucleosynthesis result: $`H\rho _0`$ is not compatible and $`C`$ is required to be small enough.

If we consider the negative bulk cosmological constant together with the positive brane tension, $`\rho _0=\mathrm{\Lambda }_0+\rho _{0M}`$, the above equation is written as

$$\left(\frac{\dot{a}_0}{a_0}\right)^2+\frac{K}{a_0^2}=(k_0^2k^2)+\frac{\mathrm{\Lambda }_0}{18M^6}\rho _{0M}+\frac{1}{36M^6}\rho _{0M}^2\frac{C}{a_0^4}.$$
(35)
The four-dimensional Friedmann equation is obtained for $`k_0^2k^2=0`$, $`C=0`$, $`\rho _{0M}\mathrm{\Lambda }_0`$. We can obtain a viable cosmology for the positive tension brane attached to the infinite size extra dimension. The effective cosmological constant is given by $`\mathrm{\Lambda }_{\mathrm{eff}}(k_0^2k^2)`$ At high energy/temperature, that is, in the very early universe $`\rho _{0M}^2`$ term dominates and results in very interesting consequences .

## 4 The stabilized RS models

### 4.1 Stabilization by balanced bulk matter

In the previous section, we saw that the brane at $`y=0`$ fixes the whole bulk solution. If we have the second brane at $`y=\frac{1}{2}`$, the junction condition at $`y=\frac{1}{2}`$ gives

$$\overline{\rho }_{\frac{1}{2}}=\frac{\mathrm{sinh}(kb)+\frac{1}{2}(\overline{\rho }_0^2\overline{C}1)\mathrm{sinh}(kb)\overline{\rho }_0\mathrm{cosh}(kb)}{\mathrm{cosh}(kb)+\frac{1}{2}(\overline{\rho }_0^2\overline{C}1)(\mathrm{cosh}(kb)1)\overline{\rho }_0\mathrm{sinh}(kb)},$$
(36)
where

$$\overline{\rho }_i=\frac{\rho _i}{6M^3k}(i=0,\frac{1}{2}),\overline{C}=\frac{C}{k^2a_0^4}.$$
(37)
This is a constraint between energy densities $`\rho _i`$ on two branes. The reason why this constraint is necessary is obvious. We used the ansatz with the static extra dimension, which is not the general case for the two brane model. But the (almost) static extra dimension is required from the phenomenological view point. We need a stabilization mechanism to make the extra dimension static.

Here we consider a simple modeling of how the stabilization mechanism works . It must involve some bulk dynamics. We suppose that it is through the role of $`C`$. The basic idea is that the back reaction by the stabilization mechanism to the brane energy densities which, if left alone, would destabilize the extra dimension, induces the bulk energy momentum $`\widehat{T}_5^5`$ which forces $`C`$ fitted to the constraint (36) through (33) and keep $`\dot{b}=0`$. The constraint (36) can be solved to give the necessary $`C`$

$$\overline{C}=(\overline{\rho }_0^21)2\frac{\overline{\rho }_0+\overline{\rho }_{\frac{1}{2}}(1+\overline{\rho }_0\overline{\rho }_{\frac{1}{2}})\mathrm{tanh}(kb)}{\mathrm{tanh}(kb)\{1\overline{\rho }_{\frac{1}{2}}\mathrm{tanh}(kb/2)\}},$$
(38)
and the bulk energy-momentum component $`p_5`$ through

$$p_5(\tau ,y)=\frac{6M^3}{4a^3(\tau ,y)\dot{a}(\tau ,y)}\dot{C}.$$
(39)
Inserting (38) into (34), we obtain the 3-brane Friedmann equation for the stabilized two brane system

$$\left(\frac{\dot{a}_0}{a_0}\right)^2+\frac{K}{a_0^2}=2k^2\frac{\overline{\rho }_0+\overline{\rho }_{\frac{1}{2}}(1+\overline{\rho }_0\overline{\rho }_{\frac{1}{2}})\mathrm{tanh}(kb)}{\mathrm{tanh}(kb)\{1\overline{\rho }_{\frac{1}{2}}\mathrm{tanh}(kb/2)\}}.$$
(40)
We have two comments here. First, in general it is expected that the stabilization mechanism induces $`\widehat{T}_i^i`$ as well as $`\widehat{T}_5^5`$, and $`b`$ may not be strictly static but be shifted somewhat as $`\rho _i`$ changes in time. Considering $`\widehat{T}_5^5`$ only and requiring strictly static $`b`$ seems to give a limit of infinitely steep stabilization potential. It is of course unrealistic, but gives the leading behaviors of RS models with a certain, unknown stabilization mechanism. Second, we can see from (39) that the induced $`p_5`$ does not fix the additive constant of $`C`$. This cannot be controlled by $`\widehat{T}_5^5`$, but required to vanish to satisfy the constraint. This may require another mechanism behind. Actually non-vanishing $`C`$ implies the breakdown of conformal symmetry in the bulk, and there might be a connection between the stabilization mechanism which requires non-trivial $`C`$ and the conformal symmetry breaking.

Now we rephrase the metric for the stabilized RS model

$$ds^2=n^2(\tau ,y)d\tau ^2+a^2(\tau ,y)\delta _{ij}dx^idx^j+b^2(\tau ,y)dy^2$$
(41)
where

$`b(\tau ,y)`$ $`=`$ $`b=\text{constant},`$
$`n(\tau ,y)`$ $`=`$ $`{\displaystyle \frac{\dot{a}(\tau ,y)}{\dot{a}(\tau ,0)}},`$
$`a(\tau ,y)`$ $`=`$ $`a_0(\tau )[\mathrm{cosh}(2kby)\overline{\rho }_0\mathrm{sinh}(2kby)`$ (42)
$`+{\displaystyle \frac{(\overline{\rho }_0+\overline{\rho }_{\frac{1}{2}})\mathrm{cosh}(kb)(1+\overline{\rho }_0\overline{\rho }_{\frac{1}{2}})}{\mathrm{sinh}(kb)\overline{\rho }_{\frac{1}{2}}\{\mathrm{cosh}(kb)1\}}}\{\mathrm{cosh}(2kby)1\}]^{1/2}`$
and $`a_0(\tau )`$ satisfies (40).

### 4.2 The static five-dimensional spacetime: The mass hierarchy and the cosmological constant

With the metric in (41) and (4.1), we first consider the case where there is no matter and only the bulk cosmological constant and the brane tensions are involved. The metric is now given by $`a(\tau ,y)=a_0(\tau )n(y)`$ and

$`n(y)`$ $`=`$ $`[\mathrm{cosh}(2kby)\overline{k}_0\mathrm{sinh}(2kby)`$ (43)
$`+{\displaystyle \frac{(\overline{k}_0+\overline{k}_{\frac{1}{2}})\mathrm{cosh}(kb)(1+\overline{k}_0\overline{k}_{\frac{1}{2}})}{\mathrm{sinh}(kb)\overline{k}_{\frac{1}{2}}\{\mathrm{cosh}(kb)1\}}}\{\mathrm{cosh}(2kby)1\}]^{1/2}.`$
From (38) and (39), the balanced $`p_5`$ is found to be

$$p_5(y)=\frac{6M^3k^2}{n(y)^4}\left[(\overline{k}_0^21)2\frac{(\overline{k}_0+\overline{k}_{\frac{1}{2}})(1+\overline{k}_0\overline{k}_{\frac{1}{2}})\mathrm{tanh}(kb)}{\mathrm{tanh}(kb)\{1\overline{k}_{\frac{1}{2}}\mathrm{tanh}(kb/2)\}}\right].$$
(44)
The scale factor $`a(\tau ,y)`$ undergoes inflation, and the Hubble parameter can be defined independently of $`y`$ because $`H=\dot{a}(\tau ,y)/a(\tau ,y)=\dot{a}_0(\tau )/a_0(\tau )`$.

The static background spacetime is obtained when we make a fine tuning to satisfy the condition for the vanishing cosmological constant

$$(\overline{k}_0+\overline{k}_{\frac{1}{2}})(1+\overline{k}_0\overline{k}_{\frac{1}{2}})\mathrm{tanh}(kb)=0.$$
(45)
With this condition, (43) and (44) are simplified to

$`n(y)`$ $`=`$ $`\left[\mathrm{cosh}(2kby)\overline{k}_0\mathrm{sinh}(2kby)\right]^{1/2},`$ (46)
$`p_5(y)`$ $`=`$ $`{\displaystyle \frac{6M^3(k_0^2k^2)}{n(y)^4}}.`$ (47)
We can identify the four-dimensional Planck scale in two ways. Firstly, we can get it from the 4-dimensional effective theory which is obtained by integrating the action over the extra dimension,

$$M_P^2=M^3_{\frac{1}{2}}^{\frac{1}{2}}b𝑑yn(y)^2=\frac{M^3}{k}\left[\mathrm{sinh}(kb)\overline{k}_0\{\mathrm{cosh}(kb)1\}\right].$$
(48)
Secondly, we can deduce it from the 3-brane Friedmann equation (40), which leads to

$$M_P^2=\frac{M^3}{k}\frac{\mathrm{tanh}(kb)\{1\overline{k}_{\frac{1}{2}}\mathrm{tanh}(kb/2)\}}{1\overline{k}_{\frac{1}{2}}\mathrm{tanh}(kb)}$$
(49)
The two derived Planck scales (48) and (49) coincide under the condition (45), that is, when the cosmological constant vanishes.

Let us consider the gauge hierarchy problem in this background spacetime. Physical mass scale at two branes are given by

$`M_W`$ $`=`$ $`Mn(\tau ,0)=M,`$ (50)
$`M_H`$ $`=`$ $`Mn(\tau ,{\displaystyle \frac{1}{2}})=M\left[\mathrm{cosh}(kb)\overline{k}_0\mathrm{sinh}(kb)\right]^{1/2}.`$ (51)
Note that we placed the visible brane at $`y=0`$ and identified the physical mass scale of the visible brane as the weak scale. The physical mass scale on the hidden brane, called the hidden scale here, is in general different from the Planck scale. Now the ratio of the electroweak scale and the Planck scale is

$$\frac{M_P^2}{M_W^2}=\left(\frac{M}{k}\right)\left[\mathrm{sinh}(kb)\overline{k}_0\left\{\mathrm{cosh}(kb)1\right\}\right]10^{32}.$$
(52)
The condition for the vanishing cosmological constant and the solution to the gauge hierarchy problem impose two conditions among four parameters $`k`$, $`k_0`$, $`k_{\frac{1}{2}}`$ and $`b`$. Hence, in this model, there are continuous set of static background spacetimes which solve the cosmological constant problem and the gauge hierarchy problem together, specified by, for example, two parameters $`(k/M,kb)`$ and $`k_0`$ and $`k_{\frac{1}{2}}`$ can be expressed in terms of them from (45) and (52)

$`\overline{k}_0`$ $`=`$ $`{\displaystyle \frac{\mathrm{sinh}(kb)10^{32}\left(k/M\right)}{\mathrm{cosh}(kb)1}},`$ (53)
$`\overline{k}_{\frac{1}{2}}`$ $`=`$ $`{\displaystyle \frac{\overline{k}_0+\mathrm{tanh}(kb)}{1\overline{k}_0\mathrm{tanh}(kb)}}.`$ (54)
The original RS model with $`k=k_0=k_{\frac{1}{2}}`$ is a special case where $`p_5`$ in (47) vanishes.

Note that while the hidden scale is larger than the weak scale by the ratio of warp factors at two branes, $`n(\tau ,\frac{1}{2})/n(\tau ,0)`$, it is in general different from the four-dimensional Planck scale by a factor $`M/k`$ and a different warp factor combination. In the original RS model where $`k=k_0=k_{\frac{1}{2}}`$, $`e^{kb}1`$ and $`M/k1`$, the difference disappears. If we introduce another hierarchy $`M/k1`$ (without any feasible reason at the moment), this will show up in the hierarchy of the hidden scale and the Planck scale.

### 4.3 One brane model

Let us turn to the cases where matters are added on the brane. First, we look at the case $`\rho _{\frac{1}{2}}=0`$. This corresponds to compactifying the extra dimension without the second brane. This is made possible by the tuned $`\widehat{p}_5`$ distribution along the extra dimension

$$\widehat{p}_5=\frac{a_0^3}{na^3}\left[6M^3k^2\frac{1}{2}k\mathrm{coth}(kb)(\rho _03p_0)\frac{\rho _0(\rho _0+3p_0)}{12M^3}\right].$$
(55)
The 3-brane Friedmann equation (40) becomes

$$\left(\frac{\dot{a}_0}{a_0}\right)^2+\frac{K}{a_0^2}=2k^2\left[1+\overline{\rho }_0\mathrm{coth}(kb)\right].$$
(56)
If we split out the brane tension from the brane energy density, the above equation takes the form of four-dimensional Friedmann equation

$$\left(\frac{\dot{a}_0}{a_0}\right)^2+\frac{K}{a_0^2}=2k^2\left[\left(\frac{k_0}{k}\mathrm{coth}(kb)1\right)+\mathrm{coth}(kb)\overline{\rho }_{0M}\right],$$
(57)
with the identifications

$`M_p^2`$ $`=`$ $`{\displaystyle \frac{M^3}{k}}\mathrm{tanh}(kb),`$ (58)
$`\mathrm{\Lambda }_{\mathrm{eff}}`$ $`=`$ $`\mathrm{\Lambda }_0(6M^3\mathrm{\Lambda }_b)^{\frac{1}{2}}\mathrm{tanh}(kb).`$ (59)
An interesting characteristic of this model is its implication for the cosmological constant problem. Suppose that we have the relation $`k=k_0`$ in some way, that is, assume a solution to the ‘big’ cosmological constant problem. Then the size of the cosmological constant has an exponential dependence on the size of extra dimension. This is just the RS-type solution to the ‘small’ cosmological constant problem. The currently observed cosmological constant $`\mathrm{\Omega }_\mathrm{\Lambda }1`$ can be fitted with $`kb140`$.

### 4.4 Two brane model

Now we turn to the two brane case. We split the brane energy densities into the brane tensions and the brane matter energy densities as we did in the one brane model, and impose the vanishing cosmological constant condition (45). The the 3-brane Friedmann equation (40) can be written as

$$\left(\frac{\dot{a}_0}{a_0}\right)^2+\frac{K}{a_0^2}=\frac{\rho _{0M}+\stackrel{~}{\rho }_{\frac{1}{2}M}+{\displaystyle \frac{\rho _{0M}\stackrel{~}{\rho }_{\frac{1}{2}M}}{6M^3k}}{\displaystyle \frac{\mathrm{tanh}(kb)}{1\overline{k}_0\mathrm{tanh}(kb)}}}{3M_P^2\left[1{\displaystyle \frac{\stackrel{~}{\rho }_{\frac{1}{2}M}}{6M_P^2k^2}}{\displaystyle \frac{\mathrm{tanh}(kb)\mathrm{tanh}(kb/2)}{1\overline{k}_0\mathrm{tanh}(kb)}}\right]},$$
(60)
where $`\stackrel{~}{\rho }_{\frac{1}{2}M}=\rho _{\frac{1}{2}M}n(\frac{1}{2})^4=\rho _{\frac{1}{2}M}\left[\mathrm{cosh}(kb)\overline{k}_0\mathrm{sinh}(kb)\right]^2`$, which is the physically observed energy density on the hidden brane. Up to small correction, this equation is nothing but the four-dimensional Friedmann equation. The energy densities on both branes contribute equally. So the matter on the hidden brane acts as dark matter for our brane.

The inclusion of the effect of stabilization mechanism, through the balanced $`\widehat{T}_5^5`$ component in this simplified model, gives a ordinary FRW cosmology at least up to TeV scale, resolving the peculiarities caused by the existence of the negative tension brane. Above the TeV scale, we meet a complicated situation where in addition that the quadratic terms become important, we need to consider the excitation of other dynamical variables which may spoil the stabilization. Further study is required to clarify it.

## 5 Conclusion

There are many interesting issues in the brane world models and large/warped extra dimensions, such as the mass hierarchy, the cosmological constant, localization of gravity, confinement of fields on the brane, fine tunings of the bulk cosmological constant and brane tensions, stabilization of the extra dimension, the role of supersymmetry, the role of AdS/CFT duality, connection to string theory, and collider signals, etc.

In this talk, we focused on the cosmological implications, together with the cosmological constant and the gauge hierarchy. In cosmological side, the model with one positive tension brane and the infinite warped extra dimension (RS2) has a viable cosmology without the need of stabilization. But for models with two or more branes or with the compact extra dimension, taking the effects of stabilization into account is very crucial in studying the cosmology of the models. We have shown that, through the simple method using the balaced $`\widehat{T}_5^5`$ component, the inclusion of the effects of stabilization recovers the conventional four-dimensional FRW cosmology. Therefore, the stabilized RS models have viable cosmology below TeV scale, with interesting new perspective on the mass hierarchy and the cosmological constant. Further study is required to clarify the cosmology of these models at temperatures above TeV scale.

## Acknowledgments

I would like to express my gratitute to the organizers of the “International Workshop on Particle Physics and the Early Universe”, COSMO2000 and acknowledge the hospitality of KIAS. I have benefited from discussions with K. Choi and H. D. Kim.