text:
# hep-th/0103240 Exact anisotropic brane cosmologies

## I Introduction

In two recent papers Randall and Sundrum  have shown that a scenario with an infinite fifth dimension in the presence of a brane can generate a theory of gravity which mimics purely four-dimensional gravity, both with respect to the classical gravitational potential and with respect to gravitational radiation. The gravitational self-couplings are not significantly modified in this model. This result has been obtained from the study of a single $`3`$-brane embedded in five dimensions, with the $`5D`$ metric given by $`ds^2=e^{f(y)}\eta _{\mu \nu }dx^\mu dx^\nu +dy^2`$, which can produce a large hierarchy between the scale of particle physics and gravity due to the appearance of the warp factor. Even if the fifth dimension is uncompactified, standard $`4D`$ gravity is reproduced on the brane. In contrast to the compactified case, this follows because the near-brane geometry traps the massless graviton. Hence this model allows the presence of large or even infinite non-compact extra dimensions. Our brane is identified to a domain wall in a $`5`$-dimensional anti-de Sitter space-time.

The Randall-Sundrum model was inspired by superstring theory. The ten-dimensional $`E_8\times E_8`$ heterotic string theory, which contains the standard model of elementary particle, could be a promising candidate for the description of the real Universe. This theory is connected with an eleven-dimensional theory, M-theory, compactified on the orbifold $`R^{10}\times S^1/Z_2`$ . In this model we have two separated ten-dimensional manifolds.

The static Randall-Sundrum solution has been extended to time-dependent solutions and their cosmological properties have been extensively studied in papers too many to cite them all  (for a recent review of dynamics and geometry of brane Universes see ). In one of the first cosmological applications of this scenario it was pointed out that a model with a non-compact fifth dimension is potentially viable, while the scenario which might solve the hierarchy problem predicts a contracting Universe, leading to a variety of cosmological problems . By adding cosmological constants to the brane and bulk, the problem of the correct behavior of the Hubble parameter on the brane has been solved by Cline, Grojean and Servant . As a result one also obtains normal expansion during nucleosynthesis, but faster than normal expansion in the very early Universe. The creation of a spherically symmetric brane-world in AdS bulk has been considered, from a quantum cosmological point of view, with the use of the Wheeler-de Witt equation, by Anchordoqui, Nuñez and Olsen .

The effective gravitational field equations on the brane world, in which all the matter forces except gravity are confined on the $`3`$-brane in a $`5`$-dimensional space-time with $`Z_2`$-symmetry have been obtained, by using an elegant geometric approach, by Shiromizu, Maeda and Sasaki . The correct signature for gravity is provided by the brane with positive tension. If the bulk space-time is exactly anti-de Sitter, generically the matter on the brane is required to be spatially homogeneous. The electric part of the $`5`$-dimensional Weyl tensor $`E_{IJ}`$ gives the leading order corrections to the conventional Einstein equations on the brane. The four-dimensional field equations for the induced metric and scalar field on the world-volume of a $`3`$-brane in the five-dimensional bulk with Einstein gravity plus a self-interacting scalar field have been derived by Maeda and Wands . The effective four-dimensional Einstein equations include terms due to scalar fields and gravitational waves in the bulk.

From a general theoretical point of view brane cosmological models are two-measure theories, since they have two independent Lagrangians, one associated to the brane, the other associated to the bulk. An other two measure theory has been proposed by Guendelman , in which the energy-momentum tensor is also a non-linear function of the general-relativistic energy-momentum tensor. Due to the scale invariance, the true vacuum state has zero energy density, when the theory is analyzed in the conformal Einstein frame. Field theory models, constructed in a gravitational theory where a measure of integration $`\mathrm{\Phi }`$ in the action is not necessary $`\sqrt{g}`$, but is determined dynamically through additional degrees of freedom, have been presented in . In this type of theory it is possible to combine the solution of the cosmological constant problem with the possibility of inflation and spontaneously broken gauge unified theories. The models with global scale invariance allow for non-trivial scalar field potentials and masses for particles, so that the scale symmetry must be broken. The conserved quantities have been obtained by Guendelman , who showed that the infrared behavior of the conserved currents is singular, so there are no conserved charges associated with scale symmetry (which implies that in some high field region the potentials become flat). For closed strings and branes ( including the supersymmetric case) the modified measure formulation is possible and does not require the introduction of a particular scale (the string or brane tension) from the beginning, but rather these appears as integration constants .

The linearized perturbation equations in the generalized Randall-Sundrum model have been obtained, by using the covariant nonlinear dynamical equations for the gravitational and matter fields on the brane, by Maartens . The nonlocal energy density determines the tidal acceleration in the off-brane direction and can oppose singularity formation via the generalized Raychaudhuri equation. Isotropy of the cosmic microwave background may no longer guarantee a Friedmann-Robertson-Walker geometry. Vorticity on the brane decays as in general relativity, but nonlocal bulk effects can source the gravitomagnetic field, so that vector perturbations can also be generated in the absence of vorticity.

The behavior of an anisotropic Bianchi type I brane-world in the presence of inflationary scalar fields has been considered by Maartens, Sahni and Saini . According to their results the magnitude of the anisotropy parameter of the brane does not affect inflation, a large initial anisotropy introducing more damping into the scalar field equation of motion. Inflation at high energies proceeds at a higher rate than the corresponding rate in general relativity. By numerical integration of field equations it is also found that anisotropy always disappear within a fixed interval of time, no matter what its initial value. The shear dynamics in Bianchi type I cosmological model on a brane with a perfect fluid has been studied in . For $`1<\gamma <2`$ the shear has a maximum at some moment during a transition period from non-standard to standard cosmology, when the matter energy density is comparable to the brane tension.

A systematic analysis, using dynamical systems techniques, of the qualitative behavior of the Bianchi type I and V cosmological models in the Randall-Sundrum brane world scenario, with matter on the brane obeying a barotropic equation of state $`p=(\gamma 1)\rho `$, has been performed by Campos and Sopuerta . In particular, they constructed the state spaces for these models and discussed what new critical points appear, the occurrence of bifurcations and the dynamics of the anisotropy. Expanding Bianchi type I and V brane-worlds always isotropize, although there could be intermediate stages in which the anisotropy grows. Near the Big-Bang the anisotropy dominates for $`\gamma 1`$, in opposition to the general relativistic situation when it dominates for $`\gamma <2`$. Frolov  has shown that the geometrical construction of the Randall and Sundrum brane world can be carried out, without the assumption of spatial isotropy, by means of an homogeneous and anisotropic Kasner type solution of the Einstein-AdS equations in the bulk. The brane equations of motion in this anisotropic space-time are solved by a static brane configuration, with the geometry of the $`3`$-brane given by the $`4D`$ Kasner solution. The set of sufficient conditions, which must be satisfied by the brane matter and bulk metric so that a homogeneous and anisotropic brane asymptotically evolves to a de Sitter space-time in presence of a positive cosmological constant have been derived in . In the presence of a non-local energy density and/or a strong anisotropic stress (i.e. a magnetic field), a Bianchi type I brane, even initially expanding, is unstable and may collapse.

It is the purpose of this paper to investigate some classes of exact solutions of the gravitational field equations in the brane world model for the anisotropic Bianchi type I and V geometries, for a conformally flat bulk (with vanishing Weyl tensor). For a perfect cosmological fluid obeying a linear barotropic equation of state, the general solution of the field equations can be obtained in an exact parametric form, for arbitrary $`\gamma `$ or in an exact analytic form for a stiff fluid. The inclusion of the quadratic terms in the energy-momentum tensor of the perfect cosmological fluids leads to major changes in the early dynamics of the anisotropic Universe, as compared to the standard general relativistic case. We also consider the dynamics of the scalar fields on the brane for different scalar field potentials of physical interest. The behavior of the observationally important physical quantities like anisotropy, shear and deceleration parameters is considered in detail for all these models.

The present paper is organized as follows. The field equations on the brane are written down in Section II. In Section III we present the general solution of the field equations for the Bianchi type I geometry with a perfect cosmological fluid. Evolution of the scalar fields in the brane-world is considered in Section IV. Bianchi type V space-times are investigated in Section V. In Section VI we discuss and conclude our results.

## II Brane geometry, field equations and consequences

In the $`5D`$ space-time the brane-world is located as $`Y(X^I)=0`$, where $`X^I,I=0,1,2,3,4`$ are $`5`$-dimensional coordinates. The effective action in five dimensions is

$$S=d^5X\sqrt{g_5}\left(\frac{1}{2k_5^2}R_5\mathrm{\Lambda }_5\right)+_{Y=0}d^4x\sqrt{g}\left(\frac{1}{k_5^2}K^\pm \lambda +L^{\text{matter}}\right),$$
(1)
with $`k_5^2=8\pi G_5`$ the $`5`$-dimensional gravitational coupling constant and where $`x^\mu ,\mu =0,1,2,3`$ are the induced $`4`$-dimensional brane world coordinates. $`R_5`$ is the $`5D`$ intrinsic curvature in the bulk and $`K^\pm `$ is the intrinsic curvature on either side of the brane.

On the $`5`$-dimensional space-time (the bulk), with the negative vacuum energy $`\mathrm{\Lambda }_5`$ as only source of the gravitational field the Einstein field equations are given by

$$G_{IJ}=k_5^2T_{IJ},T_{IJ}=\mathrm{\Lambda }_5g_{IJ}+\delta (Y)\left[\lambda g_{IJ}+T_{IJ}^{\text{matter}}\right],$$
(2)
In this space-time a brane is a fixed point of the $`Z_2`$ symmetry. In the following capital Latin indices run in the range $`0,\mathrm{},4`$ while Greek indices take the values $`0,\mathrm{},3`$.

Assuming a metric of the form $`ds^2=(n_In_J+g_{IJ})dx^Idx^J`$, with $`n_Idx^I=d\chi `$ the unit normal to the $`\chi =\text{const.}`$ hypersurfaces and $`g_{IJ}`$ the induced metric on $`\chi =\text{const.}`$ hypersurfaces, the effective four-dimensional gravitational equations on the brane take the form :

$$G_{\mu \nu }=\mathrm{\Lambda }g_{\mu \nu }+k_4^2T_{\mu \nu }+k_5^4S_{\mu \nu }E_{\mu \nu },$$
(3)
where

$$S_{\mu \nu }=\frac{1}{12}TT_{\mu \nu }\frac{1}{4}T_\mu {}_{}{}^{\alpha }T_{\nu \alpha }^{}+\frac{1}{24}g_{\mu \nu }\left(3T^{\alpha \beta }T_{\alpha \beta }T^2\right),$$
(4)
and $`\mathrm{\Lambda }=k_5^2(\mathrm{\Lambda }_5+k_5^2\lambda ^2/6)/2,k_4^2=k_5^4\lambda /6`$ and $`E_{IJ}=C_{IAJB}n^An^B`$. $`C_{IAJB}`$ is the $`5`$-dimensional Weyl tensor in the bulk and $`\lambda `$ is the vacuum energy on the brane. $`T_{\mu \nu }`$ is the matter energy-momentum tensor on the brane with components $`T^0{}_{0}{}^{}=\rho ,T^1{}_{1}{}^{}=T^2{}_{2}{}^{}=T^3{}_{3}{}^{}=p`$ and $`T=T^\mu _\mu `$ is the trace of the energy-momentum tensor.

In this paper we restrict our analysis to a conformally flat bulk geometry with $`C_{IAJB}0`$.

The Einstein equation in the bulk imply the conservation of the energy momentum tensor of the matter on the brane,

$$T_\mu {}_{}{}^{\nu }{}_{;\nu }{}^{}_{\chi =0}^{}=0.$$
(5)
In the following we investigate the Bianchi type Universes that are obtained from space-time geometries that admit a simple transitive three-dimensional group of isometries $`G_3`$ on space-like hypersurfaces. Hence the corresponding cosmological models are spatially homogeneous. There is only one essential dynamical coordinate, the time $`t`$, and the gravitational field equations reduce to ordinary differential equations. the metric on the brane is given in the general case by

$$ds^2=dt^2+\gamma _{ab}(t)e_i^a(x)dx^ie_j^b(x)dx^j,$$
(6)
where $`e_i^a(x)`$ are one-forms inverse to the spatial vector triad $`e_a^i(x)`$, which have the same commutators $`C^a{}_{bc}{}^{},a,b,c=1,2,3`$ as the structure constants of the group of isometries and commute with the unit normal vector $`e_0`$ to the surfaces of homogeneity, that is, $`e_a=e_a^i/x^i`$ and $`e_0=/t`$ obey the commutation relations $`[e_a,e_b]=C^c{}_{ab}{}^{}e_{c}^{},[e_0,e_a]=0`$.

In order to simplify the calculations we only consider geometries so that $`\gamma _{ab}(t)`$ can be taken to be a diagonal matrix. We denote by $`a_i(t),i=1,2,3`$ the components of the metric tensor on the anisotropic brane.

We define the following variable:

$`V`$ $`=`$ $`{\displaystyle \underset{i=1}{\overset{3}{}}}a_i,\text{(volume scale factor),}`$ (7)
$`H_i`$ $`=`$ $`{\displaystyle \frac{\dot{a}_i}{a_i}},i=1,2,3,\text{(directional Hubble parameters),}`$ (8)
$`H`$ $`=`$ $`{\displaystyle \frac{1}{3}}{\displaystyle \underset{i=1}{\overset{3}{}}}H_i,\text{(mean Hubble parameter),}`$ (9)
$`\mathrm{\Delta }H_i`$ $`=`$ $`H_iH,i=1,2,3.`$ (10)
From Eqs. (7), (8) and (9) we obtain $`H=\dot{V}/3V`$.

The physical quantities of observational interest in cosmology are the expansion scalar $`\theta `$, the mean anisotropy parameter $`A`$, the shear scalar $`\sigma ^2`$ and the deceleration parameter $`q`$, all are defined according to

$`\theta `$ $`=`$ $`3H,`$ (11)
$`A`$ $`=`$ $`{\displaystyle \frac{1}{3}}{\displaystyle \underset{i=1}{\overset{3}{}}}\left({\displaystyle \frac{\mathrm{\Delta }H_i}{H}}\right)^2,`$ (12)
$`\sigma ^2`$ $`=`$ $`{\displaystyle \frac{1}{2}}\sigma _{ik}\sigma ^{ik}={\displaystyle \frac{1}{2}}\left({\displaystyle \underset{i=1}{\overset{3}{}}}H_i^23H^2\right)={\displaystyle \frac{3}{2}}AH^2,`$ (13)
$`q`$ $`=`$ $`{\displaystyle \frac{d}{dt}}H^11=H^2\left(\dot{H}+H^2\right).`$ (14)
The sign of the deceleration parameter indicates whether the model inflates or not. The positive sign of $`q`$ corresponds to “standard” decelerating models whereas the negative sign indicates inflation.

## III Evolution of perfect cosmological fluids on Bianchi type I brane Universes

The line element of a Bianchi type I space-time, which generalizes the flat Robertson-Walker metric to the anisotropic case, is given by

$$ds^2=dt^2+a_1^2(t)dx^2+a_2^2(t)dy^2+a_3^2(t)dz^2.$$
(15)
We assume that the thermodynamic pressure $`p`$ of the cosmological fluid obeys a linear barotropic equation of state $`p=(\gamma 1)\rho ,\gamma =\text{const.}`$ and $`1\gamma 2`$.

Using the variables (7)-(9) the gravitational field equations and the Bianchi identity on the brane take the form

$`3\dot{H}+{\displaystyle \underset{i=1}{\overset{3}{}}}H_i^2`$ $`=`$ $`\mathrm{\Lambda }{\displaystyle \frac{3\gamma 2}{2}}k_4^2\rho {\displaystyle \frac{3\gamma 1}{12}}k_5^4\rho ^2,`$ (16)
$`{\displaystyle \frac{1}{V}}{\displaystyle \frac{d}{dt}}(VH_i)`$ $`=`$ $`\mathrm{\Lambda }{\displaystyle \frac{\gamma 2}{2}}k_4^2\rho {\displaystyle \frac{\gamma 1}{12}}k_5^4\rho ^2,i=1,2,3,`$ (17)
$`\dot{\rho }+3\gamma H\rho `$ $`=`$ $`0.`$ (18)
Thus, the time evolution of the energy density of the matter is given by

$$\rho =\rho _0V^\gamma ,\rho _0=\text{const.}>0.$$
(19)
By summing Eqs. (17) we find

$$\frac{1}{V}\frac{d}{dt}(VH)=\mathrm{\Lambda }\frac{\gamma 2}{2}k_4^2\rho \frac{\gamma 1}{12}k_5^4\rho ^2.$$
(20)
Subtracting Eq. (20) back to Eqs. (17) we obtain

$$H_i=H+\frac{K_i}{V},i=1,2,3,$$
(21)
with $`K_i,i=1,2,3`$ constants of integration satisfying the consistency condition $`_{i=1}^3K_i=0`$. By using the evolution law of the matter energy density the equation (20), the basic equation describing the dynamics of the anisotropic brane world for a conformally flat bulk, can be written as

$$\ddot{V}=3\mathrm{\Lambda }V\frac{3(\gamma 2)}{2}k_4^2\rho _0V^{1\gamma }\frac{\gamma 1}{4}k_5^4\rho _0^2V^{12\gamma },$$
(22)
and has the general solution

$$tt_0=\left(3\mathrm{\Lambda }V^2+3k_4^2\rho _0V^{2\gamma }+\frac{1}{4}k_5^4\rho _0^2V^{22\gamma }+C\right)^{1/2}𝑑V,$$
(23)
where $`C`$ is a constant of integration.

Therefore for a Bianchi type I induced brane geometry the general solution of the gravitational field equations can be expressed in the following exact parametric form, with $`V0`$ taken as parameter:

$`\theta `$ $`=`$ $`\left(3\mathrm{\Lambda }+3k_4^2\rho _0V^\gamma +{\displaystyle \frac{1}{4}}k_5^4\rho _0^2V^{2\gamma }+CV^2\right)^{1/2},`$ (24)
$`a_i`$ $`=`$ $`a_{0i}V^{1/3}\mathrm{exp}\left[K_i{\displaystyle \left(3\mathrm{\Lambda }V^4+3k_4^2\rho _0V^{4\gamma }+\frac{1}{4}k_5^4\rho _0^2V^{42\gamma }+CV^2\right)^{1/2}𝑑V}\right],i=1,2,3,`$ (25)
$`A`$ $`=`$ $`3K^2\left(3\mathrm{\Lambda }V^2+3k_4^2\rho _0V^{2\gamma }+{\displaystyle \frac{1}{4}}k_5^4\rho _0^2V^{22\gamma }+C\right)^1,`$ (26)
$`\sigma ^2`$ $`=`$ $`{\displaystyle \frac{K^2}{2}}V^2,`$ (27)
$`q`$ $`=`$ $`2{\displaystyle \frac{36\mathrm{\Lambda }V^2+18(2\gamma )k_4^2\rho _0V^{2\gamma }+3(1\gamma )k_5^4\rho _0^2V^{22\gamma }}{12\mathrm{\Lambda }V^2+12k_4^2\rho _0V^{2\gamma }+k_5^4\rho _0^2V^{22\gamma }+4C}},`$ (28)
where $`K^2=_{i=1}^3K_i^2`$.

With the use of Eqs. (19), (21) and (24), from Eq. (16) it follows that the arbitrary integration constants $`K_i,i=1,2,3`$ and $`C`$ must satisfy the consistency condition

$$K^2=\frac{2}{3}C.$$
(29)
At densities significantly greater than the nuclear one $`\rho _n`$, e.g. $`\rho >>\rho _n`$, we have $`p\rho `$, with the speed of sound $`c_s`$ tending to the speed of light, $`c_sc`$. A typical approach to the nuclear equation of state in the very high-density regime is to construct a relativistic Lagrangian that allows “bare” nucleons to interact attractively via scalar meson exchange and repulsively via the exchange of a more massive vector meson. But at the highest densities the vector meson exchange dominates and one still has $`p=\rho `$. Therefore the equation of state most appropriate to describe the high density regime of the early Universe is the stiff Zeldovich one, with $`\gamma =2`$.

For a stiff cosmological fluid, the dynamics of the matter on the brane is mainly determined by the correction terms from $`S_{\mu \nu }`$, quadratic in the matter variables and due to the form of the Gauss-Codazzi equations.

With $`\gamma =2`$ Eq. (23) becomes

$$tt_0=\frac{VdV}{\sqrt{3\mathrm{\Lambda }V^4+\alpha V^2+\beta }},$$
(30)
where $`\alpha =3k_4^2\rho _0+C`$ and $`\beta =\frac{1}{4}k_5^4\rho _0^2`$. The time dependence of the volume scale factor of the Bianchi type I Universe in the high-energy limit is therefore given by

$$V(t)=\frac{1}{2\sqrt{3\mathrm{\Lambda }}}e^{\sqrt{3\mathrm{\Lambda }}(tt_0)}\sqrt{\left(e^{2\sqrt{3\mathrm{\Lambda }}(tt_0)}\alpha \right)^212\mathrm{\Lambda }\beta }.$$
(31)
For $`t=t_S=t_0+\mathrm{ln}\left(2\sqrt{3\beta \mathrm{\Lambda }}+\alpha \right)/2\sqrt{3\mathrm{\Lambda }}`$ the volume scale factor is zero, $`V(t_S)=0`$. By reparameterizing the initial value of the cosmological time according to,

$$\mathrm{exp}(2\sqrt{3\mathrm{\Lambda }}t_0)=\alpha +2\sqrt{3\mathrm{\Lambda }\beta },$$
(32)
the evolution of the high-density brane Universe starts from $`t=0`$ from a singular state $`V_0=V(t=0)=0`$.

The time evolution of the expansion, scale factors, anisotropy, shear and deceleration parameter are given by

$`\theta (t)`$ $`=`$ $`\sqrt{3\mathrm{\Lambda }}{\displaystyle \frac{e^{4\sqrt{3\mathrm{\Lambda }}(tt_0)}\alpha ^2+12\mathrm{\Lambda }\beta }{\left(e^{2\sqrt{3\mathrm{\Lambda }}(tt_0)}\alpha \right)^212\mathrm{\Lambda }\beta }},`$ (33)
$`a_i(t)`$ $`=`$ $`a_{0i}e^{\sqrt{\frac{\mathrm{\Lambda }}{3}}(tt_0)}\left[\left(e^{2\sqrt{3\mathrm{\Lambda }}(tt_0)}\alpha \right)^212\mathrm{\Lambda }\beta \right]^{1/6}\mathrm{exp}\left[2K_iF\left(e^{\sqrt{3\mathrm{\Lambda }}(tt_0)}\right)\right],i=1,2,3,`$ (34)
$`A(t)`$ $`=`$ $`{\displaystyle \frac{12K^2e^{2\sqrt{3\mathrm{\Lambda }}(tt_0)}\left[\left(e^{2\sqrt{3\mathrm{\Lambda }}(tt_0)}\alpha \right)^212\mathrm{\Lambda }\beta \right]}{\left(e^{4\sqrt{3\mathrm{\Lambda }}(tt_0)}\alpha ^2+12\mathrm{\Lambda }\beta \right)^2}},`$ (35)
$`\sigma ^2(t)`$ $`=`$ $`{\displaystyle \frac{6K^2\mathrm{\Lambda }e^{2\sqrt{3\mathrm{\Lambda }}(tt_0)}}{\left(e^{2\sqrt{3\mathrm{\Lambda }}(tt_0)}\alpha \right)^212\mathrm{\Lambda }\beta }},`$ (36)
$`q(t)`$ $`=`$ $`12e^{2\sqrt{3\mathrm{\Lambda }}(tt_0)}{\displaystyle \frac{\alpha e^{4\sqrt{3\mathrm{\Lambda }}(tt_0)}+(\alpha ^212\mathrm{\Lambda }\beta )\left(\alpha 2e^{2\sqrt{3\mathrm{\Lambda }}(tt_0)}\right)}{\left[e^{4\sqrt{3\mathrm{\Lambda }}(tt_0)}\alpha ^2+12\mathrm{\Lambda }\beta \right]^2}}1,`$ (37)
where $`a_{0i},i=1,2,3`$ are arbitrary constants of integration and $`F(x)=\left[(x^2\alpha )^212\mathrm{\Lambda }\beta \right]^{1/2}𝑑x`$.

If the cosmological constant is zero, $`\mathrm{\Lambda }=0`$, the evolution of the brane Universe in the extreme limit of high densities is described by

$`V(t)`$ $`=`$ $`{\displaystyle \frac{1}{\sqrt{\alpha }}}\sqrt{\alpha ^2(tt_0)^2\beta },`$ (38)
$`\theta (t)`$ $`=`$ $`{\displaystyle \frac{\alpha ^2(tt_0)}{\alpha ^2(tt_0)^2\beta }},`$ (39)
$`a_i(t)`$ $`=`$ $`a_{0i}\left[\alpha ^2(tt_0)^2\beta \right]^{1/6}\left[\alpha (tt_0)+\sqrt{\alpha ^2(tt_0)^2\beta }\right]^{K_i/\sqrt{\alpha }},i=1,2,3,`$ (40)
$`A(t)`$ $`=`$ $`{\displaystyle \frac{3K^2}{\alpha ^3}}{\displaystyle \frac{\alpha ^2(tt_0)^2\beta }{(tt_0)^2}},`$ (41)
$`\sigma ^2(t)`$ $`=`$ $`{\displaystyle \frac{\alpha K^2}{2\left[\alpha ^2(tt_0)^2\beta \right]}},`$ (42)
$`q(t)`$ $`=`$ $`2+{\displaystyle \frac{3\beta }{\alpha ^2(tt_0)^2}}.`$ (43)
## IV Scalar Field Evolution on the Anisotropic Bianchi type I Brane

In the previous Sections we have considered the evolution of ordinary matter, obeying a barotropic equation of state, on the brane-like Universe. In this case the gravitational field equations can be solved exactly, with the general solution represented in an exact parametric form for arbitrary $`\gamma `$ and in an exact analytic form in the extreme limit of high densities, corresponding to $`\gamma =2`$.

As suggested by particle physics models, scalar fields, $`\varphi `$, with energy density and pressure given by $`\rho _\varphi =\dot{\varphi }^2/2+U(\varphi )`$ and $`p_\varphi =\dot{\varphi }^2/2U(\varphi )`$, respectively, where $`U(\varphi )`$ is the scalar field potential, are supposed to play a fundamental role in the evolution of the early Universe. For a Bianchi type I brane-Universe filled with a scalar field the gravitational field equations take the form:

$`3\dot{H}+{\displaystyle \underset{i=1}{\overset{3}{}}}H_i^2`$ $`=`$ $`k_4^2\left[U(\varphi )\dot{\varphi }^2\right]k_5^4\left[{\displaystyle \frac{1}{6}}\left({\displaystyle \frac{\dot{\varphi }^2}{2}}+U(\varphi )\right)^2+{\displaystyle \frac{1}{4}}\left({\displaystyle \frac{\dot{\varphi }^4}{4}}U^2(\varphi )\right)\right],`$ (44)
$`{\displaystyle \frac{1}{V}}{\displaystyle \frac{d}{dt}}(VH_i)`$ $`=`$ $`k_4^2U(\varphi ){\displaystyle \frac{1}{12}}k_5^4\left({\displaystyle \frac{\dot{\varphi }^4}{4}}U^2(\varphi )\right),i=1,2,3.`$ (45)
The scalar field also obeys the following evolution equation:

$$\ddot{\varphi }+3H\dot{\varphi }+\frac{dU(\varphi )}{d\varphi }=0.$$
(46)
Depending on the functional form of the potential $`U(\varphi )`$, several cosmological scenarios of scalar field evolution on the brane can be obtained.

As a first example we consider that the scalar field potential is a constant, $`U(\varphi )=\mathrm{\Lambda }=\text{const.}>0`$. The potential acts like a cosmological constant. Hence Eq. (46) gives $`\dot{\varphi }=2\varphi _0/V`$, with $`\varphi _0>0`$ a constant of integration. Adding Eqs. (45) we find the basic equation describing the dynamics of the constant potential scalar field on the brane:

$$\ddot{V}=\lambda V\frac{\lambda _0}{V^3},$$
(47)
where we denoted $`\lambda =3k_4^2\mathrm{\Lambda }+k_5^4\mathrm{\Lambda }^2/4`$ and $`\lambda _0=k_5^4\varphi _0^4>0`$.

The general solution of Eq. (47) is given by

$$V(t)=\frac{1}{2\sqrt{\lambda }}e^{\sqrt{\lambda }(tt_0)}\sqrt{\left(e^{2\sqrt{\lambda }(tt_0)}C\right)^24\lambda _0\lambda },$$
(48)
with $`C`$ a constant of integration. This solution is similar to the solution corresponding to the stiff matter cosmological fluid and all the other physical quantities can be obtained from Eqs. (33-37) by means of the formal substitutions $`3\mathrm{\Lambda }\lambda ,\alpha C`$ and $`\beta \lambda _0`$. Therefore, by choosing the value of $`t_0`$ such that $`\mathrm{exp}(2\sqrt{\lambda }t_0)=C+\sqrt{4\lambda _0\lambda }`$, then the brane Universe evolves from $`t=0`$ with a singular initial state.

As a second example of scalar field evolution on the brane we consider a scalar field with potential energy proportional to the kinetic one, so that $`U(\varphi )=(\gamma _\varphi 1)\dot{\varphi }^2/2,\gamma _\varphi =\text{const.}>1`$ . The energy density and pressure of the scalar field are given by $`\rho _\varphi =\gamma _\varphi \dot{\varphi }^2/2`$ and $`p_\varphi =(2\gamma _\varphi )\dot{\varphi }^2/2`$, respectively. The evolution equation for the scalar field gives $`\dot{\varphi }=2\varphi _0V^{1/\gamma _\varphi }`$. The evolution equation for the volume scale factor becomes:

$$\ddot{V}=6k_4^2\varphi _0^2(\gamma _\varphi 1)V^{12/\gamma _\varphi }k_5^4\varphi _0^4\gamma _\varphi (2\gamma _\varphi )V^{14/\gamma _\varphi },$$
(49)
and has the general solution

$$tt_0=\left[6k_4^2\varphi _0^2\gamma _\varphi V^{22/\gamma _\varphi }+k_5^4\varphi _0^4\gamma _\varphi ^2V^{24/\gamma _\varphi }+C\right]^{1/2}𝑑V,$$
(50)
with $`t_0`$ and $`C`$ arbitrary constants of integration. The time variation of the physically important parameters is given in a parametric form, with $`V`$ taken as parameter, by

$`\theta `$ $`=`$ $`\sqrt{6k_4^2\varphi _0^2\gamma _\varphi V^{2/\gamma _\varphi }+k_5^4\varphi _0^4\gamma _\varphi ^2V^{4/\gamma _\varphi }+CV^2},`$ (51)
$`a_i`$ $`=`$ $`a_{i0}V^{1/3}\mathrm{exp}\left[K_i{\displaystyle \left(6k_4^2\varphi _0^2\gamma _\varphi V^{42/\gamma _\varphi }+k_5^4\varphi _0^4\gamma _\varphi ^2V^{44/\gamma _\varphi }+CV^2\right)^{1/2}𝑑V}\right],`$ (52)
$`A`$ $`=`$ $`{\displaystyle \frac{3K^2}{6k_4^2\varphi _0^2\gamma _\varphi V^{22/\gamma _\varphi }+k_5^4\varphi _0^4\gamma _\varphi ^2V^{24/\gamma _\varphi }+C}},`$ (53)
$`\sigma ^2`$ $`=`$ $`{\displaystyle \frac{K^2}{2}}V^2,`$ (54)
$`q`$ $`=`$ $`2{\displaystyle \frac{18k_4^2\varphi _0^2(\gamma _\varphi 1)V^{22/\gamma _\varphi }+3k_5^4\varphi _0^4\gamma _\varphi (\gamma _\varphi 2)V^{24/\gamma _\varphi }}{6k_4^2\varphi _0^2\gamma _\varphi V^{22/\gamma _\varphi }+k_5^4\varphi _0^4\gamma _\varphi ^2V^{24/\gamma _\varphi }+C}},`$ (55)
$`\varphi \varphi _0^{}`$ $`=`$ $`2\varphi _0{\displaystyle \left[6k_4^2\varphi _0^2\gamma _\varphi V^{22/\gamma _\varphi }+k_5^4\varphi _0^4\gamma _\varphi ^2V^{24/\gamma _\varphi }+C\right]^{1/2}V^{1/\gamma _\varphi }𝑑V},`$ (56)
with $`a_{i0},i=1,2,3`$ and $`\varphi _0^{}`$ constants of integration. In this model the scalar field satisfies an equation of state of the form $`p_\varphi =(2/\gamma _\varphi 1)\rho _\varphi `$.

The cosmological behavior of Universes filled with a scalar field, $`\varphi `$, as well as a Liouville type exponential potential $`U(\varphi )=\mathrm{\Lambda }e^{k\varphi },\mathrm{\Lambda },k=\text{const.}>0`$, has been extensively investigated in the physical literature for both homogeneous and inhomogeneous scalar fields (for a discussion of the physical relevance of the exponential potential see  and references therein). The exponential type potential arises in the four-dimensional effective Kaluza-Klein type theories from compactification of the higher-dimensional supergravity or superstring theories. In string or Kaluza-Klein theories the moduli fields associated with the geometry of the extra-dimensions may have effective exponential potentials due to the curvature of the internal spaces or to the interaction of the moduli with form fields on the internal spaces. Exponential potentials can also arise due to non-perturbative effects such as gaugino condensation.

For the flat isotropic Friedmann-Robertson-Walker geometry a particular solution of the field equations with $`\varphi =\delta \mathrm{ln}t`$, has been obtained by Barrow . In the framework of general relativity there is no anisotropic Bianchi type I solution of this type.

We try to generalize this solution for the case of the scalar fields confined on the brane. By assuming that the constants $`k`$ and $`\delta `$ satisfy the condition $`k\delta =2`$, it is easy to check by direct calculations that the following expressions satisfy the field equations (44)-(46):

$`V(t)`$ $`=`$ $`V_0t^2,H(t)={\displaystyle \frac{2}{3}}t^1,a_i(t)=a_{i0}t^{\frac{2}{3}}\mathrm{exp}\left(K_iV_0^1t^1\right),i=1,2,3,`$ (57)
$`A(t)`$ $`=`$ $`{\displaystyle \frac{3}{4}}K^2V_0^2t^2,\sigma ^2(t)={\displaystyle \frac{1}{2}}K^2V_0^2t^4,q(t)={\displaystyle \frac{1}{2}},`$ (58)
$`\varphi (t)`$ $`=`$ $`{\displaystyle \frac{2}{\sqrt{3}k_4}}\mathrm{ln}t,U(t)={\displaystyle \frac{2}{3k_4^2}}t^2.`$ (59)
The scalar field and the potential are determined by the four-dimensional coupling constant, with $`\mathrm{\Lambda }=2/3k_4^2,k=\sqrt{3}k_4`$ and $`\delta =2/\sqrt{3}k_4`$. The integrations constants $`K_i,i=1,2,3`$ satisfy the consistency conditions $`_{i=1}^3K_i=0`$.

However, if we substitute above results to the Eq. (44), then we will find a constraint on the constants of integration $`K^2=8k_5^4V_0^2/27k_4^4`$. From the expression of $`K^2`$ we can see that for the general relativistic case, $`k_5=0`$, only isotropic solution can exist, namely $`K^2=_{i=1}^3K_i^2=0`$, leading to $`K_i0,i`$ and to a zero pressure scalar field, $`p_\varphi 0`$. For the brane world model, the constraint condition for $`K^2`$ indicates that the scalar field driven brane Universe with exponential type potential and with a simple logarithmic time dependence of the scalar field cannot exist, in general.

In the case of this simple model the scalar field potential can also be given, in the isotropic case, as a function of the volume scale factor, $`UV^1`$. It is interesting to find the general solution of the field equations for this choice of $`U`$ for the anisotropic brane, without imposing any supplementary restrictions on the physical variables. For the sake of generality we consider a more general problem, with a general functional dependence of the scalar field potential on the volume scale factor, $`U=U(V)`$. With this assumption the Klein-Gordon equation for the scalar field (46) can be integrated to give

$$\dot{\varphi }^2=\frac{2}{V^2}\left(BV^2\frac{dU}{dV}𝑑V\right),$$
(60)
and

$$\rho _\varphi =\frac{\dot{\varphi }^2}{2}+U=\frac{1}{V^2}\left[B+Y(V)\right],$$
(61)
where $`B`$ is a constant of integration and $`Y(V)=2VU(V)𝑑V`$. The general solution of the gravitational field equations on the brane can be expressed again in an exact parametric form, with $`V`$ as parameter:

$`tt_0`$ $`=`$ $`{\displaystyle \left\{\frac{3}{2}K^2+3k_4^2[B+Y(V)]+\frac{k_5^4}{4}V^2[B+Y(V)]^2\right\}^{1/2}𝑑V},`$ (62)
$`\theta `$ $`=`$ $`\sqrt{{\displaystyle \frac{3}{2}}K^2V^2+3k_4^2V^2[B+Y(V)]+{\displaystyle \frac{k_5^4}{4}}V^4[B+Y(V)]^2},`$ (63)
$`a_i`$ $`=`$ $`a_{i0}V^{1/3}\mathrm{exp}\left\{K_i{\displaystyle \left[\frac{3}{2}K^2V^2+3k_4^2V^2[B+Y(V)]+\frac{k_5^4}{4}[B+Y(V)]^2\right]^{1/2}𝑑V}\right\},i=1,2,3,`$ (64)
$`A`$ $`=`$ $`{\displaystyle \frac{3K^2}{\frac{3}{2}K^2+3k_4^2[B+Y(V)]+\frac{k_5^4}{4}V^2[B+Y(V)]^2}},`$ (65)
$`\sigma ^2`$ $`=`$ $`{\displaystyle \frac{K^2}{2V^2}},`$ (66)
$`q`$ $`=`$ $`2{\displaystyle \frac{18k_4^2VY^{}(V)+3k_5^4V^1[B+Y(V)]Y^{}(V)3k_5^4V^2[B+Y(V)]^2}{6K^2+12k_4^2[B+Y(V)]+k_5^4V^2[B+Y(V)]^2}},`$ (67)
$`\varphi \varphi _0`$ $`=`$ $`{\displaystyle \frac{\sqrt{2B+2Y(V)VY^{}(V)}dV}{\sqrt{\frac{3}{2}K^2V^2+3k_4^2[B+Y(V)]V^2+\frac{k_5^4}{4}[B+Y(V)]^2}}},`$ (68)
with $`t_0,a_{i0},i=1,2,3`$ and $`\varphi _0`$ constants of integration and $`Y^{}(V)=dY(V)/dV`$.

Suppose that the scalar field potential is given, as a function of the volume scale factor, by the simple relation $`U=u_0/V`$, with $`u_0>0`$ a constant. Then, by denoting $`\omega =3k_4^2u_0/2`$ and $`\tau =(3K^2+2k_5^4u_0^2)/12k_4^2u_0`$ and taking the integration constant $`B=0`$, the general solution of the gravitational field equations for the scalar field confined on the anisotropic brane take the form:

$`V(t)`$ $`=`$ $`\omega (tt_0)^2\tau ,\theta (t)={\displaystyle \frac{2\omega (tt_0)}{\omega (tt_0)^2\tau }},`$ (69)
$`a_i(t)`$ $`=`$ $`a_{i0}{\displaystyle \underset{j=1,2}{}}\left[\sqrt{\omega }(tt_0)+\epsilon _j\sqrt{\tau }\right]^{1/3\epsilon _jK_i/2\sqrt{\omega \tau }},i=1,2,3,`$ (70)
$`A(t)`$ $`=`$ $`{\displaystyle \frac{3K^2}{4\omega ^2(tt_0)^2}},\sigma ^2(t)={\displaystyle \frac{K^2}{2[\omega (tt_0)^2\tau ]^2}},`$ (71)
$`q(t)`$ $`=`$ $`{\displaystyle \frac{1}{2}}+{\displaystyle \frac{3\tau }{2\omega (tt_0)^2}},\varphi (t)={\displaystyle \frac{2}{\sqrt{3}k_4}}\mathrm{sinh}^1\left[\sqrt{{\displaystyle \frac{\omega }{\tau }}(tt_0)^21}\right],`$ (72)
$`U(t)`$ $`=`$ $`{\displaystyle \frac{u_0}{\omega (tt_0)^2\tau }},U(\varphi )={\displaystyle \frac{u_0}{\tau }}\mathrm{sinh}^2\left({\displaystyle \frac{\sqrt{3}k_4}{2}}\varphi \right),`$ (73)
where we have also taken $`\varphi _0=0`$ and denoted $`\epsilon _1=+1,\epsilon _2=1`$.

## V Perfect fluid Bianchi type V brane cosmologies

The line element of a Bianchi type V space-time, representing the anisotropic generalization of the open $`k=1`$ Robertson-Walker geometry, is given by

$$ds^2=dt^2+a_1^2(t)dx^2+a_2^2(t)e^{2x}dy^2+a_3^2(t)e^{2x}dz^2.$$
(74)
The effective Einstein field equations on the brane are given by

$`3\dot{H}+{\displaystyle \underset{i=1}{\overset{3}{}}}H_i^2`$ $`=`$ $`\mathrm{\Lambda }{\displaystyle \frac{3\gamma 2}{2}}k_4^2\rho {\displaystyle \frac{3\gamma 1}{12}}k_5^4\rho ^2,`$ (75)
$`{\displaystyle \frac{1}{V}}{\displaystyle \frac{d}{dt}}(VH_i)2a_1^2`$ $`=`$ $`\mathrm{\Lambda }{\displaystyle \frac{\gamma 2}{2}}k_4^2\rho {\displaystyle \frac{\gamma 1}{12}}k_5^4\rho ^2,i=1,2,3,`$ (76)
$`2H_1H_2H_3`$ $`=`$ $`0,`$ (77)
$`\dot{\rho }+3\gamma H\rho `$ $`=`$ $`0.`$ (78)
From Eq. (77) we obtain $`a_2a_3=a_1^2`$, leading to $`V=a_1^3`$. The distribution of the matter energy-density on the brane is given again by $`\rho =\rho _0V^\gamma `$. Adding Eqs. (76) we obtain

$$\frac{1}{V}\frac{d}{dt}(VH)=\mathrm{\Lambda }+2V^{2/3}\frac{\gamma 2}{2}k_4^2\rho \frac{\gamma 1}{12}k_5^4\rho ^2.$$
(79)
Subtracting Eq. (79) from Eqs. (76) leads again to

$$H_i=H+\frac{K_i}{V},i=1,2,3,$$
(80)
with $`K_i,i=1,2,3`$ constants of integration satisfying the consistency condition $`_{i=1}^3K_i=0`$.

Substitution of the Hubble parameters $`H_i,i=1,2,3`$ into Eq. (77) gives $`2K_1=K_2+K_3`$ and, consequently, we obtain

$$K_1=0,K_2=K_3.$$
(81)
The general dynamics of the Bianchi type V space-time is described by the equation

$$\ddot{V}=6V^{1/3}+3\mathrm{\Lambda }V\frac{3(\gamma 2)}{2}k_4^2\rho _0V^{1\gamma }\frac{\gamma 1}{4}k_5^4\rho _0^2V^{12\gamma },$$
(82)
with the general solution

$$tt_0=\left(9V^{4/3}+3\mathrm{\Lambda }V^2+3k_4^2\rho _0V^{2\gamma }+\frac{1}{4}k_5^4\rho _0^2V^{22\gamma }+C\right)^{1/2}𝑑V,$$
(83)
where $`C`$ is a constant of integration.

Therefore for a Bianchi type V induced brane geometry the general solution of the gravitational field equations can be expressed in the following exact parametric form, with $`V0`$ taken as parameter:

$`\theta `$ $`=`$ $`\left(9V^{2/3}+3\mathrm{\Lambda }+3k_4^2\rho _0V^\gamma +{\displaystyle \frac{1}{4}}k_5^4\rho _0^2V^{2\gamma }+CV^2\right)^{1/2},`$ (84)
$`a_i`$ $`=`$ $`a_{0i}V^{1/3}\mathrm{exp}\left[K_i{\displaystyle \left(9V^{10/3}+3\mathrm{\Lambda }V^4+3k_4^2\rho _0V^{4\gamma }+\frac{1}{4}k_5^4\rho _0^2V^{42\gamma }+CV^2\right)^{1/2}𝑑V}\right],i=1,2,3,`$ (85)
$`A`$ $`=`$ $`3K^2\left(9V^{4/3}+3\mathrm{\Lambda }V^2+3k_4^2\rho _0V^{2\gamma }+{\displaystyle \frac{1}{4}}k_5^4\rho _0^2V^{22\gamma }+C\right)^1,`$ (86)
$`\sigma ^2`$ $`=`$ $`{\displaystyle \frac{K^2}{2}}V^2,`$ (87)
$`q`$ $`=`$ $`2{\displaystyle \frac{72V^{4/3}+36\mathrm{\Lambda }V^2+18(2\gamma )k_4^2\rho _0V^{2\gamma }+3(1\gamma )k_5^4\rho _0^2V^{22\gamma }}{36V^{4/3}+12\mathrm{\Lambda }V^2+12k_4^2\rho _0V^{2\gamma }+k_5^4\rho _0^2V^{22\gamma }+4C}},`$ (88)
where $`K^2=_{i=1}^3K_i^2`$.

Substitution of the Hubble parameter and of the density into Eq. (75) gives the second consistency condition satisfied by the constants $`K_i,i=2,3`$:

$$K^2=\frac{2}{3}C.$$
(89)
Hence we can express the integration constants $`K_2,K_3`$ as functions of the constant $`C`$:

$$K_2=\pm \sqrt{\frac{C}{3}},K_3=\sqrt{\frac{C}{3}}.$$
(90)
In the limit of high matter densities ($`\gamma =2`$) the general solution of the gravitational field equations for a Bianchi type V geometry cannot be expressed in an exact analytic form.

## VI Discussions and final remarks

In the present paper we have considered some classes of exact solutions of the gravitational field equations on the brane world in the framework of the Randall-Sundrum scenario. The Randall-Sundrum mechanism was originally motivated as a possible mechanism for evading Kaluza-Klein compactification by localizing gravity in the presence of an uncompactified extra dimension. In this model, the 5-dimensional bulk space-time is assumed to be vacuum except for the presence of a cosmological constant. Matter fields on the brane are regarded as responsible for the dynamics of the brane. The quadratic terms in the energy-momentum tensor, due to the form of the Gauss-Codazzi equations, may be important in the very early stages of the Universe, leading to major changes in the dynamics of the Universe.

For a perfect linear barotropic cosmological fluid confined on the brane, the early cosmological evolution in anisotropic Bianchi type I and V geometries is fundamentally changed by the inclusion of the terms proportional to the squre of the energy density. The time variation of the mean anisotropy parameter of the Bianchi type I space-time is represented, for different values of $`\gamma `$, in Fig.1. At high densities the brane Universe starts its evolution from an isotropic state, with $`A(t_0)=0`$. The anisotropy increases and reaches a maximum value after a finite time interval $`t_{\mathrm{max}}`$, and for $`t>t_{\mathrm{max}}`$, the mean anisotropy is a monotonically decreasing function, tending to zero in the large time limit. This behavior is in sharp contrast to the usual general relativistic evolution, illustrated by the case $`\gamma =1`$ in Fig.1 (for this choice the quadratic contribution to the energy-momentum tensor vanishes), in which the Universe is born in a state of maximum anisotropy. In the initial stage the evolution is non-inflationary, but in the large time limit the brane Universe ends in an inflationary stage. In Fig.2 we present the dynamics of the deceleration parameter for different values of $`\gamma `$. In the limit of small $`t`$ (and small $`V`$, also), and with $`C=0`$, from Eq. (23) we obtain $`Vt^{1/\gamma }`$ (this result is valid for both Bianchi type I and V geometries). Therefore in the early stages of evolution of the Universe the mean anisotropy varies as $`At^{22/\gamma }`$. Near the singular point the scale factors behave like $`a_{i0}t^{1/3\gamma }\mathrm{exp}\left(\frac{\gamma }{1+\gamma }K_i^{}t^{1+1/\gamma }\right)`$. The deceleration parameter is given by $`q=3\gamma 1`$. Figs.3 and 4 present the time evolution of the mean anisotropy and deceleration parameter for the brane with Bianchi type V geometry, the general features of the cosmological dynamics being similar to the Bianchi type I case.

The general evolution of scalars fields confined on the anisotropic brane is also modified by the presence of the quadratic terms in the energy-momentum tensor. For a constant potential scalar field, the solution of the gravitational field equations is very similar to the stiff cosmological fluid case, with the cosmic evolution starting from an isotropic state, with zero anisotropy and with a maximum anisotropic phase reached after a finite time $`t_{\mathrm{max}}`$. The dynamics of the scalar field satisfying the equation of state $`p_\varphi =(2/\gamma _\varphi 1)\rho _\varphi `$, $`\gamma _\varphi >1`$, depends on the values of $`\gamma _\varphi `$. For $`1<\gamma _\varphi <3/2`$, the Universe is born in an isotropic state, but for $`\gamma _\varphi >3/2`$ the behavior is similar to the general relativistic case, with an infinite initial anisotropy, which rapidly decays to zero. For small $`\gamma _\varphi `$ the evolution is non-inflationary. In the limit of large $`\gamma _\varphi `$ the scalar field behaves like a cosmological fluid obeying an equation of state of the form $`p_\varphi =\rho _\varphi `$. In this case and for large times the behavior of the Universe is described by

$`V`$ $`=`$ $`V_0\mathrm{exp}(H_0t),H={\displaystyle \frac{1}{3}}H_0,a_i=a_{i0}\mathrm{exp}\left({\displaystyle \frac{H_0}{3}}t{\displaystyle \frac{K_i}{V_0H_0}}e^{H_0t}\right),i=1,2,3,`$ (91)
$`A`$ $`=`$ $`{\displaystyle \frac{3K^2}{V_0^2H_0^2}}\mathrm{exp}(2H_0t),q=1,\varphi =2\varphi _0t,\rho _\varphi =2\varphi _0^2\gamma _\varphi =\text{const.},`$ (92)
where $`H_0=\varphi _0\sqrt{6k_4^2\gamma _\varphi +k_5^4\varphi _0^2\gamma _\varphi ^2}`$. In the large time limit the Universe ends in an isotropic de Sitter inflationary phase. The time variations of the anisotropy and deceleration parameters are represented, for different values of $`\gamma _\varphi `$, in Figs. 5 and 6.

Scalar fields with potentials of the form $`U(\varphi )=U_0\mathrm{sinh}^\beta (\alpha \varphi ),U_0,\beta ,\alpha =\text{consts.}`$, have been recently considered in quintessence models as tracker solutions of the gravitational field equations, which can drive the Universe in its present non-decelerating phase . For such a type of potential we have obtained the general solution of the field equations in the case of the anisotropic brane world for $`\beta =2.`$ The asymptotic behavior of the potential corresponds to an inverse power-law like form for small $`\varphi `$, $`U(\varphi )\varphi ^2`$ and to an exponential one for large values of the scalar field, $`U(\varphi )e^{\sqrt{3}k_4\varphi }`$. In this model the Universe starts from a state with infinite anisotropy and ends in an isotropic phase, but the evolution is generally non-inflationary, with $`q>0`$ for all times.

It is interesting to note that in two measure gravitational theories with scale invariance a dilaton field with exponential type potentials has to be also introduced, leading to non-trivial mass generation. Interpolating models for natural transition from inflation to a slowly accelerated Universe at late times appear naturally . For cosmological models without the cosmological constant problem and theories of extended objects (strings, branes) without a fundamental scale see Guendelman .

We have not presented the solutions corresponding to scalar fields confined to the Bianchi type V brane world. Their mathematical form is very similar to the Bianchi type I case, with one correction term, proportional to $`V^{4/3}`$ added to the parametric time equation. Due to the presence of this term, exact analytic solutions cannot be generally obtained for scalar field models in Bianchi type V space-times, but the general physical behavior of scalar fields in both geometries has the same qualitative features.

## Acknowledgments

The work of CMC is supported by the Taiwan CosPA project and, in part, by the Center of Theoretical Physics at NTU and National Center for Theoretical Science.