one of the most important aspects of string theory with potential applications to beyond - the - standard - model physics and cosmology is the existence of additional extended objects .
a particularly interesting class of these , @xmath2-branes , occur as @xmath3-dimensional surfaces arising as supersymmetric solutions of the low - energy supergravity field equations but have also a high - energy interpretation as surfaces confining the endpoints of open strings @xcite . by now many different configurations of @xmath2-branes
have been studied , including parallel and intersecting sets of branes @xcite , the combination of branes with different dimensionalities @xcite and also several possible embeddings in the compactified extra - dimensions @xcite .
these models have an enormous potential in the construction of gauge theories with matter and several new extensions of the standard model have been obtained in this way @xcite .
branes have also been intensively studied in their application to the problems of early universe cosmology , such as inflation and reheating . in the context of effective field theory treatments of brane - world model building , an immense amount of work has been performed on the astrophysical and cosmological implications of these theories for both flat and warped extra dimensions @xcite , while many studies in string theory have also been undertaken @xcite .
one important property of @xmath2-branes is the fact that they interact at large distances via the exchange of closed string modes , which allows one to study them as massive point - like charged objects moving along the dimensions transverse to their world - volumes
. an interesting application of this was proposed in @xcite , where it was shown that branes with opposite charges may form bound states similar to electron - positron states , despite the different origin of the central potentials involved .
this analogy has motivated the name of _ branonium _ to designate this type of systems .
several properties of these bound states have been analyzed in @xcite , including the classical orbits , the associated quantum dynamics and possible compactification schemes .
this analysis assumed a probe ( anti)brane moving in the background spacetime created by a fixed central stack of branes , thus neglecting the gravitational backreaction of the probe . in the case of @xmath1-branes ,
the probe s trajectory follows closed elliptical orbits similar to those found in planetary systems .
it was also suggested that the probe s motion would lead to radiation of particles into both the bulk space and the probe s world - volume , making otherwise stable orbits decay .
the stability of this system was further studied in @xcite , where it was shown that orbital decay necessarily occurs when geometric moduli of the background spacetime are allowed to vary .    in this work ,
we explore some properties of branonium systems , namely the mechanism leading to particle production in the world - volume of the probe brane . for the particular case of a scalar field living in the world - volume of a probe @xmath1-brane
, we show that particle production may occur in a regime of parametric resonance similar to the one behind the phenomenon of preheating after inflation @xcite , leading to an exponential growth of the particle number in the probe s world - volume .
we determine how the orbital parameters of the system constrain this mechanism and analyze the possible properties of the resonantly produced particles .
also , we study how the probe s trajectory is modified by compactifying three of the dimensions parallel to the branes , so that a 4-dimensional spacetime with a flat friedmann - robertson - walker geometry is obtained . in this case
, we analyze the effects of hubble expansion on the orbital parameters and how it modifies the parametric resonance regime .
other instability sources such as radiation into bulk modes are discussed as well .
we also consider the effects of transverse space compactification on the probe s motion in the case of a compact 3-torus and analyze the associated production of angular momentum .
finally , we explore some of the cosmological applications of the resonant particle production mechanism and relate them to possible final configurations of the system , including stable states or brane - antibrane annihilation .    in the next section ,
we begin by introducing some of the general properties of @xmath4 branonia , with particular emphasis on the @xmath5 case , and discuss the effects of hubble expansion on the orbital parameters . in section iii , we study the evolution of scalar particle - modes confined to the probe brane and analyze the properties of the associated parametric resonance mechanism in both the non - expanding and expanding universe cases .
we distinguish between the production of massless and massive particles and determine the significance of energy damping into these modes on the probe s motion .
we also consider , in this section , the effects of radiation into bulk modes on the resonant production of brane - bound particles . in section iv , we analyze how compactifying the directions transverse to the branes modifies the probe s trajectory and dynamically generates the required angular momentum .
possible final states of the system are described in section v , where we also analyze the relevance of the particle production mechanism to cosmology in some particular scenarios .
we summarize the main results and conclusions of this work in section vi .
branonium corresponds , in its simplest form proposed in @xcite , to a system whose components are a fixed stack of @xmath6 parallel @xmath7-branes and a probe @xmath7-brane , parallel to the stack , whose gravitational backreaction is neglected so that it has no influence on the geometry of the background spacetime .
let the coordinates @xmath8 , @xmath9 , correspond to the dimensions parallel to the stack and the coordinates @xmath10 , @xmath11 , correspond to those transverse to the branes .
then , the general solution of the supergravity equations of motion in @xmath12 dimensions for @xmath6 parallel @xmath7-branes is given , in the einstein frame , by @xcite : [ eq1 ] ds^2&=&h^-dx_p^2+h^dy^2  , + e^&=&h^  , c_01p=(1-h^-1 )  , which give the @xmath12-dimensional metric , the dilaton @xmath13 and the non - vanishing components of the bulk @xmath3-form field @xmath14}$ ] which couples to the @xmath3-dimensional world - volume of the stack of branes .
the harmonic function @xmath15 is a function of the radial coordinate in transverse space , @xmath16 , and has the form : [ eq2 ] h(r)=1+q_pr^  , where @xmath17 and the factor of @xmath18 ensures that the solution is asymptotically the @xmath12-dimensional minkowski space .
for the cases we will be interested in studying , the transverse space is at least 2-dimensional , so that @xmath19 .
the stack of branes is , thus , located at the origin in the transverse space .
the exponents @xmath20 , @xmath21 , @xmath22 and the constant @xmath23 that characterize the solution are given by : [ eq3 ] & = & 4d(d-2 )  , = 4(d-2 )  , + & = & 2  ,  = 2  , where @xmath24 and [ eq4 ] = _
n^2 + 2dd-2  .
the constant @xmath25 defines the coupling of the @xmath26-form field strength to the dilaton in the supergravity action
. we will be interested in the particular case of @xmath2-branes in type ii supergravities , which arise as the low energy limit of type ii string theories , so that the @xmath27-form field strength corresponds to a closed string state in the ramond - ramond sector , and @xmath28 .
taking @xmath29 , we have @xmath30 , @xmath31 , @xmath32 , @xmath33 , @xmath34 and @xmath35 .
it is clear that the condition @xmath19 implies @xmath36 .
the constant @xmath37 is related to the charge carried by the stack of @xmath6 branes , and for the case of @xmath2-branes we have [ eq5 ] q_p = c_pg_snl_s^  , where @xmath38 , @xmath39 is the string length and @xmath40 is the string coupling constant .
this solution holds in the limit of validity of string perturbation theory and for small curvature , which allows a low energy description of the system in terms of supergravity fields .
this approximation is valid as long as the local string coupling , @xmath41 , is everywhere small .
also , we need to require the radial coordinate @xmath42 , which will give the distance between the probe and the stack branes , to be large compared to the string length @xmath39 , so that mediation by supergravity bulk fields ( graviton , dilaton and rr - forms ) is the dominant source of brane interactions .    before analyzing the motion of the probe brane , we need to consider possible compactifications of the 10 dimensional spacetime that yield our 4-dimensional world at low energies .
compactification of dimensions parallel and transverse to the world volume of the branes involve different procedures . in the case of dimensions parallel to the branes
, it suffices to consider the kaluza - klein _
ansatz_. if we start with the @xmath12-dimensional spacetime eq .
( [ eq1 ] ) and compactify along a periodic coordinate @xmath43 with radius @xmath44 , so that the branes are wrapped around this dimension , we obtain a @xmath45-dimensional spacetime with metric given by @xcite : [ eq6 ] ds_d^2&=&e^2ds_d^2+e^2(dz+b_dx^)^2  , + ds_d^2&=&h^-dx_p^2+h^dy^2  , + e^&=&h^  , where [ eq7 ] & = & 12(d-1)(d-2 )  , = -(d-2 )  , + & = & 2(d-2 )  , + &=&-2  , =+2  .
the dilaton field @xmath13 is not altered by this compactification , while the @xmath46-form rr field can be decomposed as follows : [ eq8 ] c_[n-1]=b_[n-1]+b_[n-2]dz  .
we take the following truncation of these fields : [ eq9 ] b_[n-1]=0  , ( b_[n-2])_01p-1=(1-h^-1 )  .
this truncation leads to a @xmath47-brane supergravity solution in @xmath48 dimensions , characterized by the same harmonic function of the transverse space radial coordinate @xmath15 .
applying this _ toroidal compactification _
procedure @xmath49 times we obtain @xmath50-brane solutions living in a @xmath51-dimensional spacetime , with the number of transverse dimensions remaining the same .
the process of compactifying the dimensions transverse to the branes is not as straightforward due to the lack of translational invariance which results from the presence of the branes .
in particular , to make a given transverse coordinate @xmath52 periodic , so that a similar kaluza - klein mechanism can be applied , one needs to include the appropriate  image branes " of the source branes which will ensure the invariance of the configuration under @xmath53 , @xmath54 being the radius of the compact dimension .
for example , if a source brane is located at @xmath55 , an image brane has to be considered at @xmath56 .
this implies a generalization of the harmonic function in eq .
( [ eq2 ] ) to include branes located at different points in transverse space @xcite : [ eq10 ] h()=1+_iq_p|-|^  , where the vectors @xmath57 denote the positions of the different image branes and @xmath58 .
these two compactification schemes allow us to reduce a 10-dimensional solution of the supergravity field equations to a 4-dimensional spacetime which resembles our world .
we will be interested in studying a particular case of the general branonium configuration described so far where we have a source stack of d6-branes and also a probe @xmath1 brane . as was shown in @xcite ,
when the probe carries an antibrane charge , i.e. opposite to the charge of the source branes , the system yields closed elliptical orbits which can be solved analytically , while for other values of @xmath36 the orbits fail to close .
this feature has , as we will describe later on , important consequences for particle production in the world - volume of the probe .
compactifying 3 of the dimensions parallel to the branes , we obtain effective d3-branes moving in a 3-dimensional transverse space , the shape of the orbits not being affected by this toroidal compactification scheme as mentioned before .
if the transverse dimensions are finite but their typical size is much larger than that of the compact parallel dimensions and the interbrane distance , we may consider them to be infinite and neglect the brane images described before .
although this would correspond to a 7-dimensional spacetime , from the point of view of fields confined to the world - volume of the source and probe branes it is effectively 4-dimensional .
the dynamics in the transverse space will , nevertheless , affect the 4-dimensional dynamics of the fields .    in order to make the system more similar to the observable universe
, we consider a modification of the solution eq .
( 1 ) so that , after compactification of 3 of the parallel dimensions , the 4-dimensional world - volume of the branes has a flat friedmann - robertson - walker ( frw ) geometry .
this is a particular case of the the system studied in @xcite , where the moduli governing the size of both parallel and transverse dimensions were allowed to vary . in this work , we will assume that some dynamical mechanism at the string scale ( see , e.g. @xcite ) or at some other high - energy scale ( such as susy breaking ) fixes all moduli except the scale factor associated with the 3 non - compact dimensions parallel to the branes .
hence , in the einstein frame , the complete 7-dimensional line element is given by : [ eq11 ] ds_7 ^ 2=h^-_4[-dt^2+a^2_ijdx^idx^j]+h^_4_mndy^mdy^n  , where @xmath59 , @xmath60 and @xmath61 .
the harmonic function is in this case @xmath62 .
the exponents @xmath63 and @xmath64 can be obtained using eq .
( [ eq3 ] ) and applying eq .
( [ eq7 ] ) for the compactification on a 3-torus , but their precise values will not be necessary in the subsequent calculations .
the motion of the probe brane is described by its action , which is divided in two parts .
the born - infeld contribution , which concerns the geometry of the world volume and possible gauge fields living on it , is in the string frame given by , for a generic @xmath7-brane probe @xcite : [ eq12 ] s_bi =- t_pd^p+1  e^-  , where @xmath65 are world - volume coordinates , @xmath66 gives the @xmath7-brane tension and @xmath67 is the induced metric on the probe .
the hatted metric tensor is defined in the string frame , and is related to the einstein frame one via @xmath68 with @xmath69 in the general case .
the antisymmetric field strength tensor @xmath70 refers to gauge fields describing open string modes with endpoints on the probe ( in the case of dirichlet branes ) , but in what follows we will set these fields to zero .
the wess - zumino part of the action describes the minimal coupling of the ramond - ramond @xmath3-form to the @xmath3-world volume of the probe brane and can be written as @xcite : [ eq13 ] s_wz =- qt_pc_[p+1 ]  .
the constant @xmath71 gives the sign of the charge of the probe brane , which being a bps saturated state is equal or opposite to the brane tension @xmath66 , so that @xmath72 corresponds to a brane and @xmath73 corresponds to an antibrane @xcite .
if we now consider the particular case of d6 branes in 10 dimensions , compactifying 3 of the parallel dimensions and including the scale factor describing the expansion ( or contraction ) of the 3 remaining non - compact dimensions along the brane , and choose the world - volume coordinates such that @xmath74 , we arrive at the following expression for the total action of the probe brane : [ eq14 ] s_4=-t_6v_3d^4x  a^3(t )  , where @xmath75 gives the volume of the compact 3-torus around which the d6-branes are wrapped to give the effective 4-dimensional space .
we have assumed that the transverse space coordinates are exclusively time - dependent , @xmath76 and defined the velocity of the probe as @xmath77 , with dots denoting time derivatives .
as mentioned before , the bound orbits of the probe brane can be exactly solved for the case where the scale factor is constant and it was shown in @xcite that they correspond to closed ellipses .
we now wish to analyze the effects of the varying scale factor on these orbits . in this case , a full analytical solution is hard to obtain and it is simpler to consider the evolution of the system when the probe brane moves at larges distances from the stack with small velocities .
these two conditions can be expressed as @xmath78 .
one must recall that we considered initially that the transverse dimensions can be considered infinite in extent .
hence , our analysis will hold as long as the radial coordinate is large compared to the typical length scale of the harmonic function , given by @xmath79 ( assumed larger than the string scale and the compactification scale ) , but still much smaller than the typical size of the transverse dimensions ( assumed much larger than @xmath79 ) . with this hierarchy in mind
, we may expand the probe brane s action to obtain , to lowest order , [ eq15 ] s_4-t_6v_3d^4x  a^3(t )  .
our subsequent analysis will not be affected by the constant term and we may drop it . for a probe antibrane , [ eq16 ] s_4t_6v_3d^4x  a^3(t )  .
the rotational symmetry of the problem implies that the probe s trajectory will be confined to a plane , which we may choose to be @xmath80 .
then , we may define the canonically normalized complex scalar field [ eq17 ] ( y^1+iy^2 )  , so that its action is , in the non - relativistic and large distance limit , [ eq18 ] s_=d^4x  a^3(t)(-g^__^*-v ( ) )  , where @xmath81 is the 4-dimensional flat frw metric of the spacetime where the field is defined and the potential is given by : [ eq19 ] v()=2t_6v_3h^-12t_6v_3(1- )  , so that it only depends on @xmath82 .
note that in writing the action in this form one must take into account that the spatial derivatives of the field , @xmath83 , @xmath84 , are assumed to vanish , i.e. the probe brane should remain parallel to the source branes during its motion .
the analysis of the classical motion of the probe brane is in this way reexpressed as a classical field theory problem of a complex scalar field in the background of a flat frw spacetime . before deriving the equations of motion for the field @xmath85 ,
it is useful to consider some of the symmetries of the action .
first , the rotational symmetry of the problem in the plane transverse to the branes is associated with a conserved angular momentum . in terms of the field theory approach
, this corresponds to the global u(1 ) invariance of the action under @xmath86 , where @xmath87 is a constant .
the associated noether current induces a conserved charge which gives , in the quantum theory , the particle number operator associated with the field .
writing the field as @xmath88 and defining the comoving angular momentum as @xmath89 , the comoving particle number density is given by : [ eq20 ] n = i(^*-^*)=2 ^ 2=2l  , while the total particle number , which corresponds to the conserved noether charge , is given by @xmath90 .
the energy density and pressure of the field can be obtained from its energy - momentum tensor , yielding : [ eq21 ] & & t_00=||^2+v()=^2+l^2
^ 2+v ( )  , + p & & t_iia^2=||^2-v()=^2+l^2 ^ 2-v ( )  .    the covariant conservation of the energy - momentum tensor can then be expressed as @xmath91 , where the energy of the field is defined as @xmath92 and @xmath93 is the volume of the expanding ( or contracting ) flat frw spacetime . defining the hubble parameter @xmath94 , we can write this in the form : [ eq22 ] ddt=-3h(+p )  .    from eqs .
( [ eq21 ] ) and ( [ eq22 ] ) we conclude that @xmath95 , the equality holding only for @xmath96 .
thus , in an expanding universe ( @xmath97 ) , the system evolves with a strictly decreasing energy density .
the angular momentum evolves according to : [ eq23 ] dldt+3hl=0  , which corresponds to conservation of the total particle number @xmath6 .
hence , in an expanding universe , the absolute value of the angular momentum will necessarily decrease , vanishing asymptotically .
we , thus , conclude from this simple analysis of conservation laws that the expansion of the universe will make the system reduce both its energy density and its angular momentum . to analyze the details of this evolution , we need to determine the equations of motion for the field @xmath85 . varying the action eq .
( [ eq18 ] ) with respect to @xmath98 , and taking this variation to vanish , we obtain : [ eq24 ] + 3h+v(||)^*=0  . in terms of the transverse space coordinates @xmath10 , @xmath99 , this can be written as : [ eq25 ] y^m+3hy^m+2q_6y^mr^3=0  .
this implies that the polar variables @xmath100 and @xmath101 satisfy : [ eq26 ] & + 3h - l^2 ^ 3 + 12 ^ 2=0  , + & + 2 + 3h=0  , where we have defined the constant [ eq27 ] 4q_6(t_6v_32)^3/2  , so that , apart from constant factors , the large distances potential is @xmath102 .
it is easy to check that the equation for @xmath101 simply gives the evolution of the angular momentum that we obtained previously in eq .
( [ eq23 ] ) .
the evolution of the radial field @xmath100 is determined by the form of the hubble parameter , given by the energy density content of the universe via the friedmann equation in the usual way .
this may include not only the energy density of the field @xmath85 but also all other matter , radiation or vacuum energy components .
the possibility that this scalar field may drive ( slow - roll ) inflation if , in the early universe , it dominates the energy density has been analyzed in @xcite , for the case where the probe brane moves in a linear trajectory , i.e. @xmath103 .
it was shown that the interbrane potential arising from type ii supergravity / string theories is not flat enough to produce the required number of e - foldings , exhibiting an  @xmath104-problem " .
this is intrinsically related to the assumption that @xmath42 is much smaller than the size of the transverse dimensions .
an inflationary period driven by the angular field variable @xmath101 was proposed in @xcite , where it is argued that , despite having a flat potential , this field may provide the constant energy density necessary for inflation if the probe brane moves in a circular orbit with @xmath105 . however , we need to take into account the fact that , during inflation , the angular momentum decays exponentially as @xmath106 , from eq .
( [ eq23 ] ) .
the system will hence quickly tend to the @xmath103 case , where sufficient inflation is difficult to obtain .
furthermore , the decay of the angular momentum will necessarily alter the evolution of @xmath100 and the value of @xmath107 , so that a slow - roll inflationary mechanism with the probe brane moving in a circular orbit would be hard to construct - branes moving in warped throats @xcite .
in particular , the important role of the angular variables in providing accelerated periods of expansion was discussed in @xcite . ] .    although we do not wish to completely discard such inflationary mechanisms
, we will from now on assume that some other field is responsible for inflation and analyze how the compactified @xmath108 branonium system evolves in the post - inflationary eras .    before analyzing in detail the effects of the expansion on the trajectory of the probe ,
let us recall the general bound orbits of the system in a non - expanding universe . in this case
, the system reduces to the well - known problem of a particle in a central @xmath109 attractive potential , admitting closed orbit solutions of the form : [ eq28 ] ( ) = r(1-e^2)1+e  .
these are closed ellipses with eccentricity and semi - major axis given by : [ eq29 ] e&=&(1 + 4l^2 ^ 2)^12  , + r&=&-2  .
the properties of such orbits may be obtained by studying the effective potential @xmath110 .
this potential has a minimum at @xmath111 , tends to @xmath112 at the origin and to zero at large distances .
the condition @xmath113 then implies that orbits with @xmath114 will be bounded , with @xmath115 . in this case
, @xmath100 will oscillate between its minimum and maximum values , @xmath116 and @xmath117 , which can be obtained from eq .
( [ eq28 ] ) by setting @xmath118 and @xmath119 , respectively . at the minimum of @xmath120 ,
the orbits will be circular with @xmath121 , and @xmath122 . for @xmath123 , the orbits will only be bounded from below with e.g. @xmath124 for @xmath125 , which corresponds to an @xmath126 parabolic orbit .
the linear trajectory , with @xmath103 , is a particular case of the latter with no bounds on the value of @xmath100 except for the trivial @xmath127 . for @xmath128 ,
the orbits are hyperbolic .
consider now the expanding case with the scale factor evolving as @xmath129 .
this power law behavior is typical of the post - inflationary stages of the universe s evolution , where one may consider a single fluid to give the dominant contribution to the total energy density .
for example , in the radiation era we have @xmath130 and in the matter era @xmath131 .
we can also take @xmath131 at the end of inflation , when the oscillations of the inflaton about the minimum of its potential , with an equation of state corresponding to that of non - relativistic matter , dominate the energy density .
this model only holds away from the transitions between these periods as in these cases at least two of the components give similar contributions to the energy density .
we will consider the evolution of the branonium system for arbitrary @xmath87 , starting at some instant @xmath132 when the scale factor has a value @xmath133 .
the hubble parameter is then of the form @xmath134 , which implies from eq .
( [ eq23 ] ) that the angular momentum evolves according to : [ eq30 ] l = l_0(tt_0)^-3  , with @xmath135 being its initial value .
it is then clear that @xmath136 as @xmath137 , as we have mentioned earlier .
the equation for the radial field @xmath100 can now be written as : [ eq31 ] + 3t - l_0 ^ 2 ^ 3(tt_0)^-6 + 12 ^ 2=0  .
let us start by analyzing how the system evolves when placed initially in a would - be circular orbit , so that [ eq32 ] ( t_0)=2l_0 ^ 2_c0  , as @xmath138 .
it is easy to see that @xmath100 can not remain constant for @xmath139 due to the decay of the angular momentum .
we may , however , look for solutions where the condition for circular orbits is maintained during the motion of the probe brane , i.e. [ eq33 ] ( t)=2l^2(t)=_c0(tt_0)^-6_c(t )  .    substituting into eq .
( [ eq31 ] ) we obtain [ eq34 ] _
c0t_0 ^ 2(tt_0)^-6 - 26(3 + 1)=0  , which is satisfied only in the non - expanding case , @xmath140 , and for @xmath141 .
as we are interested in cases where @xmath142 , we conclude that @xmath143 is not an exact solution of the equations of motion .
it is , however , an approximate solution at late times and , as we will check later , gives the global evolution of the system in an expanding universe , so that it is useful to consider some of its properties .
in particular , the angular frequency evolves as : [ eq35 ] _ c = l_c^2(tt_0)^9  , so that the orbital period decreases as @xmath144 .
such a solution eventually fails to satisfy the non - relativistic and large distance approximation as the radius of the orbit decreases and its angular velocity grows . in particular ,
the parameter @xmath145 which controls this approximation grows , for @xmath146 , as [ eq36 ] hv^2(t_6v_3)^4(q_6l_0)^4(tt_0)^12  .
hence , given the values of @xmath147 and @xmath75 , it is the ratio @xmath148 that controls how long the approximations remain valid .
consider now trajectories with a non - vanishing initial eccentricity . in an expanding universe ,
the effective potential becomes time - dependent , although its asymptotic properties at the origin and at infinity remain the same . as the angular momentum redshifts away , the minimum of the effective potential decreases , the same happening with the value of @xmath100 at this minimum , according to the circular solution defined in eq .
( [ eq33 ] ) .
as discussed earlier , the energy density of the field always decreases during the motion of the probe brane in an expanding universe . hence ,
if the system is initially in one of the closed orbits with @xmath149 , the latter condition will be satisfied for all @xmath139 .
this means that , at all times , the probe brane will be in one of the closed elliptical orbits defined earlier . from eq .
( [ eq29 ] ) , we also conclude that the semi - major axis of the orbit will decrease as the energy density evolves to more negative values .
the eccentricity of the orbits may , however , not remain constant , as it depends on the variation of both the energy density and the angular momentum .
in fact , using eqs .
( [ eq22 ] ) and ( [ eq23 ] ) , one can show that the eccentricity varies according to [ eq37 ] de^2dt=-12hl^2 ^ 2(3+p)=-24hl^2 ^ 2(2t+v )  .
where we defined the kinetic energy of the field as @xmath150 , so that @xmath151 and @xmath152 .
thus , in an expanding universe , for finite @xmath135 , one expects the eccentricity of the orbit to vary unless @xmath153 .
this last condition is satisfied on average in the non - expanding case , corresponding to the virial theorem .
it does not necessarily hold , however , in an expanding universe , with the semi - major axis of the orbits decaying as discussed above .
it is nevertheless clear that , at late times , as @xmath107 and @xmath154 decrease , the variation of the eccentricity of the orbits should be smaller .
the sign of @xmath155 is also not definite , so that the eccentricity may either increase or decrease during the motion of the probe brane .    from eqs .
( [ eq29 ] ) and ( [ eq33 ] ) , we obtain : [ eq38 ] r(t)=_c(t)1-e^2(t )  .
hence , the variation of the semi - major axis follows the decrease of the circular solution obtained earlier , being also affected by the variation of the eccentricity .    in order to have a better understanding of the evolution of the probe brane s motion in an expanding universe ,
we have solved the equations of motion numerically .
measuring all quantities in terms of the string length , i.e. setting @xmath156 , we choose , as an example , the values @xmath157 and @xmath158 we follow the example analyzed in @xcite . ] .
we set the initial conditions at @xmath159 to be those of a non - expanding elliptical orbit with angular momentum @xmath135 and eccentricity @xmath160 , which we take to be the only free orbital parameters .
the field is initially at its maximum value @xmath161 , with @xmath162 and @xmath119 . after determining the numerical solution for @xmath163 , we computed the corresponding energy density and eccentricity evolution , according to eqs .
( [ eq21 ] ) and ( [ eq29 ] ) , respectively .
we used the numerical solution for the eccentricity to compute @xmath164 , according to eq .
( [ eq38 ] ) , and also the maximum and minimum values of @xmath163 at each orbit , given by : [ eq39 ] _
min(t)&=&r(t)(1-e(t ) )  , + _
max(t)&=&r(t)(1+e(t ) )  .
we have also determined the evolution of the angular field variable @xmath101 using @xmath165 .    in figure 1
we have plotted the results obtained for the values @xmath166 and @xmath167 in a matter - dominated universe , @xmath131 , with @xmath168
. seems quite large , one must take into account that it is measured in units of the string time , which is of order @xmath169 seconds if the string scale is of the order of the planck scale .
thus , we are considering the probe brane to be moving quite early in the history of the universe . ]     and @xmath167 in a matter - dominated universe , @xmath131 , with @xmath168 .
the plots show ( a ) the radial field and associated quantities , ( b ) the energy density , ( c ) the eccentricity and ( d ) the motion of the probe brane in the complex plane of the field @xmath85.,title="fig : " ] [ fig : figure_1 ]    observing the plots shown in this figure , we conclude that , as expected , the probe brane evolves continually through elliptical orbits of decreasing semi - major axis , so that they fail to close .
the radial field then oscillates between minimum and maximum values which decrease in time .
the global decrease follows that of the circular solution , eq .
( [ eq33 ] ) , but exhibits an oscillating behavior associated with the variation of the energy density .
although the latter strictly decreases , as expected , it also oscillates with decreasing period and amplitude , so that at later stages its evolution tends to be smooth . consequently
, one observes oscillations in the evolution of the eccentricity of the orbit . in the example shown in figure 1
, these oscillations have a small amplitude .
however , at earlier times , when the effects of hubble expansion are more significant , numerical simulations show that this amplitude can be quite large .
all simulations show that @xmath170 starts increasing and , as the period and amplitude of the oscillations decay , it tends to a constant value .
a constant eccentricity at late times implies that , in this limit , the energy density should , on average , vary as [ eq40 ] l^-2(t)(tt_0)^6  , tending smoothly to more negative values as observed in figure 1 .
then , from eq .
( [ eq22 ] ) , [ eq41 ] = -3h=6h  , so that @xmath171 or , equivalently , @xmath172 .
thus , although at early times the effect of the expansion makes the eccentricity vary significantly , at late times the system _ virializes _ and starts evolving smoothly between orbits of constant eccentricity and decaying radius .
this result will be useful later for studying particle production in the probe brane at late times .
another interesting consequence of the expansion of the universe is a rotation of the axis of the elliptical orbits with time , although this is a small effect in the example shown in figure 1 .
this is simply a result of the growing eccentricity of the orbits which deviates the maximum and minimum values of @xmath100 from @xmath119 and @xmath118 , respectively .
the orbital axis will , however , stabilize at late times , as @xmath170 tends to a constant value .
numerically , one also observes a decrease in the amplitude of the oscillations of @xmath163 .
this is given by [ eq42 ] ( t)_max(t)-_min(t)=r(t)e(t )  , so that the observed decay of the amplitude is mainly due to the decrease of @xmath164 as discussed above .
the asymptotic value of the eccentricity depends not only on the initial time @xmath132 at which the probe starts its motion but also on the initial values of the orbital angular momentum and eccentricity .
numerical simulations show that larger values of @xmath135 lead to a more pronounced growth of the eccentricity , in agreement with eq .
( [ eq37 ] ) .
the dependence on @xmath160 is less trivial and numerically one finds a larger eccentricity variation for @xmath160 close to 0 and to @xmath173 .
this is , however , very small for @xmath160 close to 1 , which is expected , as the energy density of the field can not increase to produce hyperbolic orbits with @xmath174 if initially @xmath149 .
the eccentricity always grows if initially @xmath105 , so that @xmath175 and @xmath176 , as most of the eccentricity variation occurs initially when hubble expansion is more significant .
the opposite behavior should be observed if initially @xmath177 .
although the example we have shown refers to a matter - dominated universe , @xmath131 , our qualitative discussion and analytical results hold for all @xmath142 , in particular to the @xmath130 radiation - dominated universe .    to summarize the results of this section ,
we conclude that , in an expanding universe , the probe brane follows elliptical orbits with decreasing radius and increasing frequency and whose eccentricity exhibits an oscillating behavior and asymptotically tends to a constant value larger than the initial one .
if the motion of the probe brane begins at late times with a small angular momentum , the latter effect is negligible .
the probe brane will asymptotically collide with the branes in the stack , as @xmath178 , but our approximations will break down before this happens , so that we expect the motion of the probe brane to deviate significantly from our previous results as its angular velocity grows and it spirals towards the central stack . the collision will , however , occur if nothing else prevents the probe from losing energy and angular momentum as the universe expands , and may lead to the annihilation of the probe antibrane with one of the source branes . for now , we will use the results obtained in this section to study the production of particles in the probe brane and we will return to this issue in section v of this work .
in the original discussion of the branonium system @xcite , it was argued that the fields confined to the probe brane , such as the gauge fields mentioned earlier , become time - dependent due to the motion of the probe through the background spacetime created by the central stack .
consequently , an observer bound to the probe interprets this variation as corresponding to the production of particles associated with these fields .
such radiation of energy into brane particle - modes arises only if the distance between the probe and the stack varies in time , thus creating a time - dependent background from the point of view of brane - bound observers .
although the power radiated into these modes was estimated in @xcite , many aspects of this particle production mechanism remain unclear . in this section
, we will analyze this mechanism in more detail , revealing some new properties of particle production in branonium , and discuss how the expansion of the universe modifies this process .
the fields confined to the brane arise , from the string theory point of view , as states associated with open strings whose endpoints are attached to the brane , as discussed earlier . for a single @xmath7-brane , which is a @xmath179 bps supersymmetric state ,
these include bosonic fields such as gauge bosons and scalar fields , as well as their fermionic superpartners .
if the probe corresponds itself to a stack of @xmath180 parallel @xmath7-branes , other degrees of freedom associated with a @xmath181 vector supermultiplet may arise @xcite . here
, we will study the simplest case of a real scalar field @xmath104 confined to the probe brane and follow an effective field theory approach
.    for the @xmath108 system described earlier , we write the effective action for the scalar field @xmath104 in the form : [ eq43 ] s_=-t_6d^7   e^-(-12^__-12m^2 ^ 2 )  , where @xmath182 denotes the mass of the scalar field , which includes its bare mass as well as quantum corrections , and all other quantities are those defined earlier in this work .
when compactifying 3 of the dimensions parallel to the brane as before , the scalar field will be decomposed in its kaluza - klein ( kk ) modes , leading to a tower of 4-dimensional massive modes , or kk states .
we will focus on the evolution of the zero mode of the field @xmath104 , whose mass will simply be given by @xmath182 .
to simplify the notation , we will denote this mode as @xmath104 , although one must bear in mind that it is not the original 7-dimensional field .
the global factor of @xmath147 will not affect our discussion and we may absorb it into the definition of the field .    first , let us consider the case @xmath183 .
then , the classical equations of motion arising from the action eq .
( [ eq43 ] ) are given by : [ eq44 ] _
( e^-^_)=0  .    using the results obtained in section ii , and defining the function @xmath184^{1\over2}$ ] , we can write this as [ eq45 ] + ( 3h - f)-1a^2(1-hv^2)^2=0  , where @xmath185 and @xmath186 is the flat 3-dimensional laplacian of the field .    to construct the associated quantum description ,
we follow the semi - classical approach to the quantization of scalar fields in curved space @xcite . in a curved background spacetime with a time - varying geometry , as is the case of the world - volume of the probe brane , the induced time - variation of the fields modifies the usual canonical quantization procedure , as the quantum operators become themselves time - dependent .
in particular , the creation and annihilation operators associated with the field will now evolve in time as the background changes , the same happening with the associated multi - particle states .
this is the main reason behind the production of particles in a dynamical background .
let us start by expanding the field in fourier modes of the form : [ eq46 ] ( , t)=(a__(t)e^i+ a^_^*_(t)e^-i )  .    in the quantum theory , @xmath187 and
@xmath188 become the annihilation and creation operators associated with the fourier mode @xmath189 of the field .
expanding the field in this way , we include all the time dependence of the field in the functions @xmath190 and @xmath191 , while the operators remain time - independent .
the conjugate momentum to the field @xmath104 , obtained by computing @xmath192 , can then be written as [ eq47 ] ( , t)=a^3f(a__e^i+ a^_^*_e^-i ) .
if the creation and annihilation operators satisfy the canonical commutation relations [ eq48 ] [ a_,a_]=[a^_,a^_]=0  ,  [ a_,a^_]=^3(- )  , then the canonical commutators for the field and its conjugate momentum , @xmath193=i\delta^3(\mathbf{x}-\mathbf{y})$ ] , can only be obtained if the mode functions satisfy the following wronskian normalization condition : [ eq49 ] _
^*_-^*__=ifa^3  .
this condition is essential for consistency of the quantization procedure .
we now follow a semi - classical approach and consider the evolution of the field modes to be given by the classical equations of motion , assuming that any quantum corrections to their propagation can be neglected .
substituting eq .
( [ eq46 ] ) into eq .
( [ eq45 ] ) , we see that each of the fourier modes evolves independently according to [ eq50 ] _
+ ( 3h - f)_+k^2a^2(1-hv^2)_=0  . and similarly for the complex conjugate modes @xmath194 .
notice that this equation depends only on @xmath195 but not the direction of the momenta @xmath196 .
this is a consequence of the isotropy of our metric _ ansatz _ , and we may write without loss of generality @xmath197 .
it is also clear , from eq .
( [ eq50 ] ) , that the physical momentum of the modes is related to their comoving momentum @xmath198 via [ eq51 ] k_phys = ka  .
this expression includes both the expected redshift of the modes in an expanding universe and the effect of the probe s motion through the dynamical background .
let us now rescale the field modes by defining the mode functions @xmath199 .
these satisfy [ eq52 ] _
k+_k^2x_k=0  , which corresponds to the equation for a harmonic oscillator with a variable frequency given by : [ eq53 ] _ k^2&=&k^2a^2(1-hv^2)-14(3h - f)^2 - 12(3h - f)= + & = & k_phys^2+^2  , where we have defined [ eq54 ] ^2(t)-14(3h - f)^2
- 12(3h - f )  .    as
each mode behaves as a harmonic oscillator , we can define the associated particle number via [ eq55 ] e_k&=&12|x_k|^2 + 12_k^2|x_k|^2=_k(n_k+12 ) + n_k&=&_k2(|x_k|^2_k^2+|x_k|^2)-12  .
it is well known that this quantity is an adiabatic invariant for a harmonic oscillator with variable frequency , as happens , for example , in an oscillating pendulum whose length is decreased infinitely slowly @xcite .
thus , the number of _ quanta _ in a given fourier mode @xmath198 can only change if its frequency varies in a non - adiabatic way , which can be expressed by the following condition @xcite : [ eq56 ] |d_kdt|_k^2  .    for our particular case , this gives , in the limit @xmath78 , [ eq57 ] |hk_phys^2-|(k_phys^2+^2)^32  .
this condition will then constrain the number of particles produced in each physical momentum mode by the motion of the probe brane .    to quantify this number of particles , we need to look for solutions of the mode equation , eq . ( [ eq52 ] ) .
we will start by discussing what happens in a non - expanding universe and then we will analyze the effects of hubble expansion on particle production .
consider the motion of the probe brane in a non - expanding universe , @xmath96 , and assume that the non - relativistic and large distance approximations discussed earlier are valid . as we have seen in section ii , the probe follows closed elliptical trajectories if its energy ( density ) is negative , @xmath149 .
if @xmath78 , we have to leading order : [ eq58 ] 1-hv^21  , k_physk  , fh^1/2
. then , we may take the following approximations : [ eq59 ] f&=&ff  , + f&&12  , so that we can write : [ eq60 ] ^2=-14f^2 + 12f-516(hh)^2 + 14hh  .
an exact analytical expression for the time variation of the radial field @xmath163 can be obtained in a non - expanding universe , being given implicitly by [ eq61 ] ( t)=r  .
the angular variable @xmath200 satisfies [ eq62 ] t - t_0=^-1  , where @xmath201 and @xmath202 is the angular frequency of the orbit . the angular variables @xmath203 and @xmath101 are related by : [ eq63 ] =  .    for small eccentricity orbits , @xmath204 , eq .
( [ eq62 ] ) gives , to leading order in @xmath205 : [ eq64 ] _ 0+(t - t_0 )  .
choosing @xmath206 we have , in this approximation , [ eq65 ] ( t)r  .
this allows us to determine the non - relativistic and large distance expression for @xmath207 , in the case of small eccentricity orbits .
to lowest order , we then obtain : [ eq66 ] ^2e^2((t - t_0 ) )  . defining [ eq67 ] ^2e^2  , the variable frequency of the mode with comoving momentum @xmath198 can ,
hence , be approximated by : [ eq68 ] _ k^2(t)k^2+^2((t - t_0 ) )  .    if we now rewrite eq .
( [ eq52 ] ) in terms of the rescaled time variable @xmath208 , we obtain : [ eq69 ] x^_k+(a_k-2q(2z))x_k=0  , where [ eq70 ] a_k  , q-2 ^ 2 ^ 2  .
( [ eq69 ] ) has the form of the well - known mathieu equation @xcite with parameters @xmath209 and @xmath71 determined by the orbital parameters of the probe brane and by the value of the comoving momentum @xmath198 of each mode . in the parameter space @xmath210 , the mathieu equation exhibits both stable and unstable solutions , the latter being closely associated with the phenomenon of parametric resonance @xcite . as we are working under the assumption that the probe brane s orbit has a large radius and a small eccentricity , it is easy to check that @xmath211 ( the analysis is independent of the sign of @xmath71 ) . in this region of parameter space , the mathieu equation exhibits instabilities in a series of narrow resonance bands near @xmath212 , with @xmath27 being a positive integer , and with a width in comoving momentum space approximately given by @xmath213 . in these resonance bands ,
the solution evolves as @xmath214 , with a real exponent @xmath215 , according to floquet s theorem .
the most important of these is the first resonance band , which occurs when [ eq71 ] 1-|q|-18q^2a_k1+|q|-18q^2  .
the exponent @xmath216 can then be approximately written as [ eq72 ] _ k  .
this is , as claimed , real for @xmath217 , having a maximum value of @xmath218 for @xmath219 .
this implies that particle - modes lying inside this resonance band will be exponentially amplified , leading to a resonant production of particles .
thus , from eq .
( [ eq70 ] ) , the center of the resonance band occurs for comoving momentum @xmath220 , while its lower and upper bounds are given by [ eq73 ] k_min&&k_c  , + k_max&&k_c  , so that the resonance band has a width in comoving momentum space of @xmath221 .
this is quite small compared to the value of the center of the resonance , as @xmath222 , which means that particle production occurs in a regime of narrow parametric resonance .
the particles produced by this mechanism have typical energies of order @xmath223 , i.e. with half the typical energy of the field @xmath85 .
this agrees with the discussion in @xcite , where it is stated that particles are produced in pairs with opposite momenta , as if resulting directly from the decay of the interbrane distance field - particles ensures momentum conservation . ] .
we have found in this work , however , that these particles are produced resonantly , the associated field modes being exponentially amplified .    after the modes inside the resonance band have been amplified sufficiently , we may write , from eq .
( [ eq55 ] ) , [ eq74 ] n_ke^2_kz = e^_k(t - t_0 )  , where we used that @xmath224 , so that the term involving @xmath225 in eq .
( [ eq55 ] ) can be neglected .
this also implies that , although the particle number in each mode grows exponentially , the resonance takes a long time to develop .
in particular , the typical resonance time , @xmath226 , is much larger than the orbital period of the probe brane , which is of order @xmath227 .
this is characteristic of the narrow resonance regime and , hence , a significant particle number is only produced after the probe brane has completed a large number of orbits , given by @xmath228 .
to have an idea of the order of magnitude of these quantities , consider the scenario studied in section ii , where @xmath229 with @xmath158
. then , for an orbit with eccentricity @xmath230 and an angular momentum @xmath231 , we have @xmath232 , which satisfies the large distances approximation , and @xmath233 .
this implies that the probe brane will have to complete approximately @xmath234 orbits for the particle number to be amplified by one e - folding .    to illustrate the results obtained earlier , we have solved eq .
( [ eq52 ] ) for the example described above .
an initial vacuum state with zero particle number corresponds to the following initial conditions : [ eq75 ] x_k(t_0)=1  , _
k(t_0)=-i  .    in figure 2
we have plotted the numerical results obtained for the comoving momentum at the center of the resonance band , @xmath235 .
the linear evolution observed for the logarithm of the particle number shows that , as expected , this quantity is being amplified exponentially with time , but that a significant number of particles is produced only after the probe brane has completed a few thousands of orbits .     and
@xmath231 , giving @xmath236 .
the time coordinate @xmath237 is given in units of the orbital period @xmath238.,title="fig : " ] [ fig : figure_2 ]    we note that , although this solution was obtained numerically , a formal analytical solution of the mathieu equation can be written in terms of the mathieu sine and cosine functions .    in figure 3
we also show the evolution of the particle number for a mode outside the limits of the resonance band . in this case
, @xmath239 takes imaginary values and @xmath240 exhibits an oscillating behavior with a small amplitude , confirming that only for comoving momenta @xmath241 a significant particle number is produced by the elliptical motion of the probe brane .
it is also worth mentioning that , in the limit of circular orbits , @xmath242 , the parameter @xmath71 vanishes and there is no amplification of any of the particle momentum modes .
this confirms our early assumption that the interbrane brane distance needs to vary for particle production to occur .     and @xmath231 , giving @xmath236 . the comoving momentum of the mode is @xmath243 .
the time coordinate @xmath237 is given in units of the orbital period @xmath238.,title="fig : " ] [ fig : figure_3 ]    so far we have not included the backreaction effects of the produced particles on the probe s motion .
we will return to this later in section iii c , where we estimate the energy density damped into brane - modes by the resonance and discuss how this affects the probe brane s orbital decay .
consider now the case where the field has a non - vanishing mass , @xmath244 . in this case , the classical equations of motion eq .
( [ eq44 ] ) can be generalized to [ eq76 ] _
( e^-^_)-e^-m^2=0  , or , explicitly , [ eq77 ] + ( 3h - f)-1a^2(1-hv^2)^2+(1-hv^2)h^12m^2=0  .
each of the fourier modes then evolves independently according to [ eq78 ] _ k+(3h - f)_k+  .
this implies that all modes have , in this case , an effective physical mass which is different from the mass parameter @xmath182 due to the motion of the probe brane : [ eq79 ] m_eff^2=(1-hv^2)h^12m^2  .
notice that the mass is not , however , affected by the expansion of the universe .
when writing the equations of motion in terms of the rescaled mode function @xmath245 defined earlier , we again obtain the equation for a harmonic oscillator with a varying frequency , which is now given by : [ eq80 ] _
k^2=k_phys^2+m_eff^2+^2  , where all the quantities are time - dependent in the general case .
let us now set @xmath96 and compare these results with the ones obtained for the particle production mechanism in the massless case .
first , notice that , in the non - relativistic and large distance approximation , we have @xmath246 .
thus , the leading order modification introduced by the mass of the field can be obtained by simply replacing @xmath247 .
following the same procedure as before , one can reduce the particle - mode equation to the mathieu equation , with the parameter @xmath209 being now given by : [ eq81 ] a_k=4(k^2+m^2)^2  .
the parameter @xmath71 , which quantifies the strength and width of the resonance phenomenon , is not modified by the introduction of the mass parameter , so that we expect the particle production mechanism to occur as before . in this case , however , the center of the resonance band is shifted to a lower momentum , @xmath248 , which gives @xmath249 .
the upper and lower limits of the resonance band are also modified , being now given by : [ eq82 ] k_min=  , + k_max=  , which in turn alters the width of the resonance band in the obvious way . if , as in the massless case , we have @xmath250 , we may write @xmath251 , which is larger than the corresponding value in the massless case , satisfying nevertheless the narrow resonance condition .
notice that this implies an upper bound for the mass of the particles which can be produced by this mechanism , [ eq83 ] m  . in the particle description discussed earlier
, this simply means that each  particle " in the interbrane distance field @xmath85 needs to have enough energy to decay into two @xmath104-particles at rest .
a residual kinetic energy is also necessary if all the modes inside the resonance band are to be excited . in the non - relativistic and large distance approximation , the orbital frequency of the probe brane
is assumed to be small , so that only particles with small masses may be excited during its motion .
one must recall , however , that we are measuring all quantities in terms of the typical string parameters .
in particular , the orbital frequency should be small compared to the string energy scale .
if the latter corresponds to the planck scale , then the produced particles can still be quite massive .
for the example studied earlier in this section , with orbital parameters @xmath252 , particles up to masses of order @xmath253 gev can be resonantly created .
recalling the definition of @xmath85 in terms of the physical interbrane distance , eq .
( [ eq17 ] ) , we may write the upper bound on the mass of the produced particles in terms of the physical semi - major axis of the orbit , @xmath254 . inserting back the missing @xmath39 factors , we obtain : [ eq84 ] m_max = m_s  , where @xmath255 is the string energy scale . as we assumed in the beginning that @xmath256
, it is clear from this expression that @xmath257 , which allows nevertheless for the production of very massive particles . in general ,
particles living in the world - volume of the probe brane will be massive , so that we need to use the results obtained above to describe the parametric resonance .
for masses which are small compared to the typical orbital energies of the probe , the simpler treatment of the @xmath183 case will be sufficient .
let us now consider the case of an expanding universe , with @xmath258 .
the inclusion of the scale factor will modify our previous analysis of particle production in several ways .
first , it will change the frequency of each harmonic oscillator particle - mode , namely modifying the factor @xmath259 , given in eq .
( [ eq54 ] ) .
also , the physical momenta of the particle - modes will be redshifted as the universe expands , so that , in the non - relativistic and large distance approximation , we have @xmath260 .
finally , the motion of the probe brane will be altered , as we have seen in section ii , so that the semi - major axis of elliptical orbits will decrease as @xmath261 , making the angular frequency consequently increase as @xmath262 .
furthermore , if the hubble parameter is sufficiently large , the orbital eccentricity may exhibit a significant growth .
all the modifications make a complete analytical study of particle production in the probe brane rather complex . in order to determine the leading effects of the expansion
, we will consider the motion of the probe brane at sufficiently late times , so that some of the effects of the expansion can be discarded . as we have concluded in section ii , the variation of the eccentricity of the orbits is , in this regime , very small , so that it may be neglected .
the explicit effects of the expansion in modifying the factor @xmath207 may also be discarded in this limit , provided that @xmath263 and @xmath264 , or explicitly [ eq85 ] t(r_physq_6 )  , t^2(r_physq_6 )  .
the first of these conditions will , _ a priori _ , be more constraining . for the example we have been following in this section , the first condition gives @xmath265 while the second only implies the constraint @xmath266 .    with these constraints in mind
, we now have to analyze how hubble expansion modifies the factors @xmath267 and @xmath268 in the variable harmonic oscillator frequency .
the radial interbrane distance field @xmath100 can now be written , in the low eccentricity orbit approximation , as : [ eq86 ] ( t)r(t )  , with [ eq87 ] r(t)=r_0(tt_0)^-6  , ( t)=_0(tt_0)^9  .
computing its time derivatives , we find that they involve terms arising from hubble expansion which are suppressed by at least one power of @xmath237 , being negligible at late times .
there are also terms involving the quantity @xmath269 .
if the time necessary for particle production to occur is smaller than the age of the universe at that time , we may also discard these terms .
after some algebra , the leading order modification to @xmath207 is then given by [ eq88 ] ^2(t)^2(t)((t)(t - t_0 ) )  .
the quantity @xmath270 is now time - dependent and can be obtained trivially by including the appropriate time variation of the orbital parameters in its non - expanding expression , eq .
( [ eq67 ] ) .
this gives [ eq89 ] ^2(t)=_0 ^ 2(tt_0)^24  , _ 0 ^ 2e_0 ^ 2  .
then , taking into account the momentum redshift , we can write the particle - mode equation , to lowest order , as : [ eq90 ] _ k+x_k=0  .    following the same reasoning as in the non - expanding case , we rescale the time coordinate by defining the variable @xmath271 . as the angular frequency of the orbit is now time - dependent , this change of variables will introduce new terms in the equation .
these will , however , be suppressed by powers of @xmath237 and involve first and second derivatives of the mode function @xmath272 . in the non - expanding case , these were quite small compared to @xmath272 , and one expects them to be negligible in this case as well . using that , to lowest order , @xmath273 , we may write the leading order equation for the mode function in terms of the variable @xmath274 as [ eq91 ] x^_k+[a_k(z)-2q(z)(2z)]x_k=0  , with [ eq92 ] a_k(z)a_k0(1-z )  ,
q(z)q_0(1+z )  , where @xmath275 and @xmath276 correspond to the initial values of the parameters of the mathieu equation , given by their non - expanding case expressions , and the coefficients @xmath20 and @xmath277 are given by [ eq93 ]  , = 310  .
we , hence , conclude that the leading order effect of the expansion is to make the coefficients of the mathieu equation vary in time .
in particular , their variation is , to a first approximation , linear in the variable @xmath274 and controlled by the single parameter @xmath20 , which will be small at late times .
one must take into account , however , that eq .
( [ eq91 ] ) is not the mathieu equation and , namely , does not satisfy the conditions of floquet s theorem which allowed us to define the coefficient @xmath239 .
nevertheless , we can use the results obtained in the non - expanding case from the properties of the mathieu equation as a guide to describe the behavior of the modes when @xmath97 .
first , notice that , as @xmath209 decreases with @xmath274 , the position of the center of the resonance band , at @xmath249 , will be shifted to higher momentum modes : [ eq94 ] k_c(z)=k_c01-z  , k_c0  .
thus , we expect modes with comoving momentum larger than @xmath278 to be excited as the probe s orbit decays . to better understand the evolution of each mode , let us analyze the evolution of the exponent @xmath279 , which is given by : [ eq95 ] _
the zeros of @xmath279 occur for [ eq96 ] z_1=-|q_0|+a_k0 - 1a_k0+|q_0|  , z_2=|q_0|+a_k0 - 1a_k0-|q_0|  , with @xmath280 .
this means that @xmath279 will be real for @xmath281 and pure imaginary for @xmath282 and @xmath283 .
figure 4 illustrates the typical evolution of @xmath279 for a generic fourier mode .     for a generic mode .
the shaded area corresponds to the period the mode spends inside the resonance band.,title="fig : " ] [ fig : figure_4 ]    thus , each mode may only experience the resonant regime during a finite amount of time , when @xmath279 is real and the mode is inside the resonance band . during the periods where @xmath279 is pure imaginary , we expect the particle number to exhibit an oscillating behavior . in the non - expanding case , the frequency and amplitude of these oscillations for modes outside the resonance band
is determined by the value of @xmath284 , so that the larger the latter the smaller the period and the amplitude of the oscillations .
as , in an expanding universe , @xmath285 changes in time , we expect the frequency and amplitude of the oscillations of the particle number to vary as the probe brane s orbit decays . similarly , as @xmath286 is not constant when the mode is inside the resonance band , the strength of the resonance is also expected to vary .
particle modes will then have different behaviours according to their position relative to the initial resonance band , which can be parametrized by : [ eq97 ] a_k01+|q_0|  .
recalling that the centre of the original resonance band corresponds to @xmath287 and that the probe s motion begins at @xmath288 , we need to distinguish three different types of modes :    1 .
@xmath289 + as both @xmath290 and @xmath291 are negative in this case , these modes will always be outside the resonance band , with an imaginary exponent @xmath279 for all @xmath292 .
2 .   @xmath293 + this case corresponds to the particle - modes inside the resonance band in a non - expanding universe , with @xmath294 but @xmath295 .
hence , they will start inside the resonance band but stop experiencing the resonant regime for @xmath283 .
3 .   @xmath296 + although these modes are outside the initial resonance band , they will experience the resonant behavior at some later stage , as @xmath297 .    for both type ii and type iii modes ,
the maximum value of @xmath298 occurs for @xmath299 , giving , for @xmath300 , [ eq98 ] _
k(z_max)(1 + 310|q_0| )  .
this gives the same result as in the non - expanding limit for @xmath287 .
it is easy to check that @xmath301 as @xmath302 , the same happening for @xmath290 and @xmath291 .
this implies that only modes for which @xmath294 ( @xmath303 ) and @xmath295 ( @xmath304 ) will be excited in this limit , in agreement with the fact that only type ii modes are excited in the non - expanding case ( @xmath305 ) and that these experience the parametric resonance during an infinite period of time .
it is also clear from eq .
( [ eq98 ] ) that higher momentum modes are excited by the resonance at later times , as expected from the evolution of the center of the resonance band .
the total time each type ii or type iii mode spends inside the resonance band is given by : [ eq99 ] z_bandz_2-z_1(1 - 710|q_0| )  .
it is clear that this quantity decreases with @xmath20 , so that the smaller the hubble parameter the more time each mode is excited by the resonance . in the limit @xmath302 , @xmath306 , but as previously discussed this limit only applies to type ii modes , as all the others are outside the resonance band in this limit .
one also concludes that higher momentum modes will spend less time inside the resonance band , which suggests these will be less excited .
the analysis of the properties of the exponent @xmath279 thus provides a very clear insight on the qualitative evolution of all momentum modes .
we would like , however , to be able to quantify the particle number produced by the resonance in each mode . in the non - expanding case , @xmath307 , which implied @xmath308 in the limit of large particle number for type ii modes .
despite the variation of @xmath239 in the @xmath97 case , we may take this exponent to be approximately constant during an infinitesimal interval between @xmath274 and @xmath309 .
then , during this interval , @xmath272 will be amplified by a factor @xmath310 if the mode is inside the resonance band . integrating this result
, we expect that , after a significant number of particles has been produced , [ eq100 ] n_k(z)(2_z_i^z_k(z)dz )  .
the initial time @xmath311 refers to the time the mode enters the resonance band , i.e. @xmath312 and @xmath313 for type ii and type iii modes , respectively .
one can use eq .
( [ eq100 ] ) to estimate the total particle number produced in each mode by the parametric resonance .
the details of this calculation are given in appendix a. for type ii and type iii modes , one obtains : [ eq101 ] n_k^res ( ii)(z_2)&=&q_0 ^ 2211+|q_0|(+2 + + & + & )  , + n_k^res ( iii)(z_2)&=&q_0 ^ 221+|q_0|  .
( [ eq101 ] ) constitutes the main result of this section , giving the leading order expressions for the particle number density produced by the resonance in an expanding universe . from these approximated expressions
, we can see that the parameter @xmath314 controls the strength of the resonance , a significant number of particles being produced only if @xmath315 .
recall that , in the non - expanding case , the typical resonance time was given by @xmath316 .
the typical time a mode spends inside the resonance band in the expanding case is , from eq .
( [ eq99 ] ) , of order @xmath317 .
hence , apart from numerical factors , the condition for a significant number of particles to be produced is simply stating that the time a given mode spends inside the resonance band should be greater than the typical time the resonance takes to develop , @xmath318 .    from the results in eq .
( [ eq101 ] ) , we can also conclude that the total particle number produced increases with @xmath319 for @xmath320 . on the other hand , for @xmath296
, it exhibits a slow decrease with @xmath319 , showing that the production of high momentum modes is suppressed .
this agrees with the fact that high momentum modes spend less time inside the resonance band , as obtained above .      to check the results obtained so far for the resonant particle production in an expanding universe , we have solved eq .
( [ eq91 ] ) numerically . in order to simplify the computation
, we have rescaled the mode functions via @xmath321 , so that initially we have approximately [ eq102 ] _
k(0)1  , ^_k(0)-i  .
we then computed the function [ eq103 ] _
k(z)=|_k(z)|^2+|^_k(z)|^2  .
this is related to the physical particle number density of each mode , given by eq .
( [ eq55 ] ) , approximately by @xmath322 , taking the limit of large particle number .
note that @xmath323 initially , as we neglected the @xmath179 factor in eq .
( [ eq55 ] ) , so that the effects of this normalization have to be taken into account before the resonance produces a significant number of particles .    in figure 5
, we illustrate the results obtained for modes of types i , ii and iii , including also the analytical prediction obtained from eq .
( [ eq100 ] ) .
these simulations correspond to an original resonance strength @xmath324 , the effects of the expansion being quantified by @xmath325 ( recall eq .
( [ eq93 ] ) ) .
these values ensure the validity of the approximations made and also that a significant particle number is produced for modes of types ii and iii .     and @xmath325 .
the plots correspond to ( a ) a type i mode , with @xmath326 , ( b ) a type ii mode , with @xmath287 and ( c ) a type iii mode with @xmath327 . the solid line corresponds to the numerical solution in all three cases , while in ( b ) and ( c ) the dashed line gives the corresponding analytical prediction .
the shaded area corresponds to the period that each mode spends inside the resonance band . ]    observing the plots shown in this figure , we see that all three particle - modes follow the expected behavior . the type i mode exhibits an oscillating particle number , with oscillations of decreasing amplitude and period , confirming that , in fact , no net particle number is produced reflects the normalization chosen in eq .
( [ eq103 ] ) and does not correspond to a physical particle number . ] .
both the type ii and the type iii modes are exponentially amplified during a finite period , after which their particle number oscillates with decreasing amplitude and period , tending to an adiabatically constant value . also , as expected , the type iii mode exhibits an oscillating behavior before entering the resonance band
.    we can also conclude from figure 5 that the particle number follows the predicted evolution inside the resonance band , the main differences between the numerical and the analytical solutions occurring near the endpoints of the resonance band , where the transitions between resonant and oscillating regimes take place .
the discrepancies for low particle number are also due to the normalization of @xmath328 , as discussed earlier .
it is nevertheless clear that eq .
( [ eq100 ] ) gives a quite good description of the resonant particle production regime in an expanding universe , although it underestimates the particle number density , and that the contribution of the non - resonant periods are subdominant , as expected .
the total particle number produced by the end of the resonant regime , @xmath329 , was computed for several values of @xmath319 , keeping @xmath330 and @xmath20 fixed at the values chosen above .
the results we have obtained are shown in figure 6 .    .
plot ( a ) shows the results for small values of @xmath319 , including type ii and iii modes , and plot ( b ) shows the results for large values of @xmath319 , where only type iii modes are present .
the solid line corresponds to the numerical results while the dashed line gives the corresponding analytical prediction .
all results correspond to @xmath324 and @xmath325.,title="fig : " ] [ fig : figure_6 ]    one observes that @xmath329 increases with @xmath319 for type ii modes and slowly decreases with @xmath319 for type iii modes , in agreement with eq .
( [ eq101 ] ) .
the main discrepancies between the numerical results and the analytical prediction are again due to the smooth transition between the resonant and non - resonant regimes at @xmath291 , with the latter giving a subdominant contribution to the total particle number , which is more significant for type iii modes .
the different number of oscillations the modes undergo before entering the resonance band give the oscillations that can be observed in figure 6 .
these are suppressed for large @xmath319 as the pre - resonance oscillations of the particle number give a negligible contribution for high momentum modes .
numerical simulations also show that , as expected , the particle number in all modes decreases with @xmath20 , so that hubble expansion may inhibit the resonant amplification of particle modes at early times .
thus , we conclude that our analysis of the resonance regime gives a good description of the mechanism of particle production in branonium for an expanding universe and that the corrections arising from the non - resonant regimes are subdominant .
extrapolation of this discussion to high momentum modes is , however , difficult , as these will only be excited at late times , where the approximations we have considered may no longer hold .
the production of high momentum particles is nevertheless suppressed , as we have concluded earlier .
furthermore , at each given time @xmath331 , only the modes with @xmath332 have entered the resonance band , which means that only those modes for which [ eq104 ] a_0k <
a_0k^max(z ) have started being excited .
this gives the following momentum cut - off for modes which have already been amplified at @xmath139 : [ eq105 ] k_max(t)=12 ^ 12(1 - 20)^-12  .
notice that this gives the expected endpoint of the resonance band in the limit @xmath333 , given by eq .
( [ eq73 ] ) .
the generalization of these results for massive particles follows the same reasoning used in the non - expanding case , with the substitution @xmath334 within our approximations .
note that , in the expanding case , the physical momentum of the particles is redshifted while their mass remains the same as the universe expands .
this modification will then produce a more complex evolution of the parameter @xmath335 in the modified mathieu equation , eq .
( [ eq91 ] ) , given to lowest order by : [ eq106 ] a_k(z)=a_k0(1-z)+4m^2_0 ^ 2(1-z )  , where @xmath336 and all other quantities are defined as before .
although this will alter the way the resonance band moves into higher momentum modes as the probe s orbit decays , the evolution of the modes is expected to follow the same qualitative behavior as in the massless case , so that we will not analyze this in more detail .
we , thus , conclude that , to lowest order , the parametric resonance survives in an expanding universe , although modes are excited only during a finite period .
these results hold only at late times but we note once again that all quantities are measured with respect to the typical string values .
we need , however , to bear in mind that eq .
( [ eq91 ] ) is not the mathieu equation and that several other terms in eq .
( [ eq52 ] ) may become important at early times .
these may alter the excitation time and the amplification of each mode , or even prevent any particle production .
it is nevertheless clear that the effects of the expansion are suppressed at the typical energies involved in the probe s motion , so that one still expects a significant number of particles to be produced .
as the resonance develops , the probe brane loses energy to excited @xmath104-particle modes and one may wonder whether the energy radiated by this mechanism is sufficient to affect the probe brane s motion .
let us start by computing the energy density of massless particles produced by the resonance in a non - expanding universe .
we have seen that each rescaled mode function @xmath272 has an associated harmonic oscillator energy function @xmath337 . for the mode functions @xmath338 ,
the associated energy function is approximately the same in the non - relativistic and large distance limit , where @xmath339 .
hence , the total energy density radiated into brane particle - modes is given by : [ eq107 ] _ p = e_k=12 ^ 2_0^+dk  k^2_k(n_k+12 )  , where we have used the momentum space isotropy discussed earlier in this work .
after a significant number of particles has been produced , we may take @xmath340 and neglect the @xmath179 factor in eq .
( [ eq107 ] ) .
also , we know that particles are produced in a narrow resonance band centered at @xmath341 with a band width @xmath342 . the particle number distribution in momentum space can then be well approximated by a gaussian distribution of the form : [ eq108 ] n_k(t)=n_k_c(t)(-12(k - k_c)^2k(t)^2 )  .
the gaussian distribution width @xmath343 should correspond to a fraction of the resonance band width , so that we may write @xmath344 .
the coefficient @xmath345 is expected to be time - dependent as the central mode @xmath235 is more amplified than all the other modes in the resonance band , so that the distribution should become more sharply peaked about this value with time .
the advantage of writing the particle number distribution in this form is that it allows us to further approximate the result by a @xmath346-function distribution , taking into account that @xmath347 , as discussed earlier . introducing the correct normalization
, we may then write : [ eq109 ] n_k(t)n_k_c(t)k(t)(k - k_c )  . using @xmath348
, we obtain : [ eq110 ] _ p & & ( 2)^3k_resn_k_c + & & |q|(2)^4 n_k_c  , where the time dependence of @xmath349 , @xmath345 and consequently @xmath350 is implicit . notice that , when @xmath351 , we have @xmath352 , as the resonance vanishes in this limit . to have an estimate of the order of magnitude of this quantity ,
take the example we have considered before in the discussion of particle production in a non - expanding universe , for which @xmath353 and @xmath354 , in units of the string length . by solving the equations of motion numerically and fitting the obtained particle number distribution to the gaussian distribution given in eq .
( [ eq108 ] ) , we may determine the values of @xmath345 and @xmath349 .
for an initial time @xmath355 , we obtained the following results : [ eq111 ] n_k_c(5t_0)&&7.510 ^ 2  , ( 5t_0)0.20  , + n_k_c(10t_0)&&1.610 ^ 7  , ( 10t_0)0.12  .
these two examples confirm that , indeed , the resonance width decreases in time . from eq .
( [ eq110 ] ) , we obtain for these two cases @xmath356 and @xmath357 , in units of the string energy scale .    the energy density radiated into brane particle - modes should then be compared with the energy density of the interbrane distance field @xmath85 , which is given by @xmath358 .
for the example considered above , @xmath359 , so that the value obtained for the energy density radiated into particles up to @xmath360 is negligible when compared to the probe brane s energy density . in general , we obtain : [ eq112 ] _ p||(2)^103 n_k_c  , which gives the main result of this subsection .
this means that , in the previous example , if we take @xmath361 , up to a @xmath104-particle number density of @xmath362 can be produced without affecting significantly the motion of the probe brane .
these results can be easily generalized for the massive case , by taking into account the changes in @xmath235 and @xmath363 discussed earlier . for @xmath250
, we obtain : [ eq113 ] _
p||k_c(2)^73 n_k_c  , with @xmath364 and which , as expected , reduces to eq .
( [ eq112 ] ) in the limit @xmath365 . as the momentum of the produced particles decreases with @xmath182 , it is easy to conclude that the energy damped into brane particle - modes is smaller for more massive particles ( recall that the strength of the resonance is not affected by the mass of the particles ) .
the analysis of the energy radiated into brane particle - modes is more difficult to perform in the expanding universe case , as more effects have to be taken into account .
first , we need to recall the approximate relation between @xmath328 and the physical particle number .
next , we need to notice that the mode functions @xmath366 are redshifted by a factor of @xmath367 with respect to the harmonic oscillator mode functions @xmath272 . finally , for the purposes of determining the energy density of the produced particles , we may take @xmath368 .
then , we have [ eq114 ] _ p(t)_k  .
we see that the energy density of the produced particles is redshifted by the usual power of @xmath369 , but is further reduced due to the redshift of the physical momentum and energy of the modes .
we can use the expressions obtained for the particle number produced by the resonance in eq .
( [ eq101 ] ) to determine @xmath328 , as we have concluded that the effects of the non - resonant periods are subdominant . writing them in terms of the comoving momentum of the modes ,
we conclude that , for type iii modes , [ eq115 ] _
p(t)dk  k^3(q_0 ^ 2_0 ^ 281k^2 )  .    to determine the limits of this integral
, one must recall that , at a given time , only a finite number of modes has entered the resonance band . in particular , one may use the value computed in eq .
( [ eq105 ] ) as the upper limit of this integral , although some of these modes have not yet been completely excited by the resonance . the lower limit of the integral will be the first type iii mode , with @xmath370 .
this gives a finite energy density radiated into brane modes at each finite time @xmath139 .
this energy density is redshifted by hubble expansion , but at the same time modes with increasing energy enter the resonance band .
this process will , however , stop at some point , as our approximations will break down when the motion of the probe becomes relativistic and gets too close to the central stack . even if the resonance mechanism persists in this limit , although certainly in a different form , the probe will eventually modify its trajectory , either by colliding with the central stack or via some other mechanism that stabilizes its motion .
the contribution of type ii modes is obviously finite and can be computed by integrating eq .
( [ eq114 ] ) using the corresponding expression for @xmath328 given in eq .
( [ eq101 ] ) , with a lower limit @xmath278 and an upper limit @xmath371 for the integral .
no particular insight is , however , gained by computing the exact expressions for the energy density for both type ii and type iii modes .
instead , we note that the dominant contribution to the energy density will be given by high momentum modes which enter the resonance band at late times . for these , the effects of the expansion of the universe will be less significant and one may use the results obtained in the non - expanding case to compute the energy density produced by the resonance . in this case
one should use the values of @xmath71 and @xmath372 at the end of the resonance mechanism .      as discussed in @xcite , the probe s trajectory in the background created by the central stack
may induce not only the production of brane - bound particles but also radiation into bulk modes , namely gravitational , rr - form and dilatonic fields .
this effect is due to the accelerated motion of the probe brane and most of the power is radiated into lower - spin fields , in particular scalar fields which couple to the brane s orbital monopole moment .
the probe will thus lose energy through this process , contributing to the decay of its orbit .
this will certainly modify the mechanism of resonant production of particle - modes in the brane , namely by varying the orbital frequency , which determines the comoving momentum of the modes that are excited by the parametric resonance .
one then expects the effects of radiation into bulk modes to be quite similar to those induced by the universe s expansion , so that we will not analyze this process in detail .
nevertheless , it is important to estimate its contribution to the decay of the probe s orbit .
the power radiated into bulk scalar modes was estimated in @xcite , where it was shown that the number of orbits the probe can complete before the interbrane distance becomes of order @xmath39 is approximately given by : [ eq116 ] ( r_il_s)^321(g_s^3n)^12  , where @xmath373 is the initial value of the physical radius of the orbit , assuming it is circular .
( this gives a good estimate for the decay time , even though we are interested in small eccentricity orbits for the resonance process .
@xmath374 corresponds to the ratio of the decay time to the initial orbital period . )
this expression shows that , for the probe to complete a large number of orbits before decay , it needs to be at a large distance from the central stack , in units of the string length .
also , one needs @xmath375 to obtain a sufficiently large value of @xmath374 .
this value should be compared to the number of orbits of the probe necessary for the resonance to be effective , which in the non - expanding case is of order @xmath376 .
hence , the production of a significant number of particles requires @xmath377 , which gives for the initial radius of the orbit : [ eq117 ] r_il_s(32|q|)^23n^13g_s  .
for example , if @xmath378 and there are 10 branes in the central stack , we need the initial interbrane distance to be larger than about @xmath379 for the resonance to be effective , which is not too large a number taking into account that @xmath40 is parametrically small .    as , in the analysis of the parametric resonance mechanism , we have assumed the interbrane distance to be large compared to the string length , we expect the effects of radiation into bulk modes to be initially negligible .
this process should , however , become more important as the orbit decays and increases its acceleration .
it is even possible that , at late times , it overcomes hubble expansion as the main energy loss process .
nevertheless , the end of the resonance should still be determined by the break down of the non - relativistic and large distance approximation , as this also controls the amount of energy damped into bulk modes .
so far we have considered the probe to move at distances from the central stack which are small compared to the typical size of the transverse space directions , so that we may neglect the effects of compactification on the probe s motion . in this section , we will analyze the leading order effects introduced by the finite size of these directions , and show that they may lead to the creation of orbital angular momentum , a necessary condition for resonant particle production to take place .
we will consider the simplest case of compactifying the three transverse directions on an isotropic 3-torus of size @xmath380 , which is defined by the identifications @xmath381 .
as mentioned earlier , the harmonic function associated with the central stack configuration needs to be modified in this case by including the appropriate
brane images " , according to eq .
( [ eq10 ] ) .
these images are placed at points with coordinates @xmath382 , with integer @xmath383 , defining a hypercubic lattice corresponding to the covering space of the 3-torus , @xmath384 .
although toroidal compactifications are too simple to give realistic particle physics phenomenology , this will be sufficient to illustrate the main effects of compactification on the probe brane s trajectory .    the interbrane potential associated with the generalized harmonic function eq .
( [ eq10 ] ) is obtained by computing the propagator for a massless field on a torus , which satisfies : [ eq118 ] ^2g(,)=(-)-1v _  , where @xmath385 is the volume of the transverse 3-torus .
the term @xmath386 is included for consistency , so that the integral over the compact manifold of both sides of eq .
( [ eq118 ] ) vanishes . if one expands the green s function @xmath387 in terms of eigenfunctions of the laplacian operator @xmath388 , one concludes that the @xmath386 term removes the unphysical zero - mode which makes the green s function diverge .
this term arises naturally from the curvature of the non - compact 4-dimensional spacetime @xcite .
it has been shown that one can write the massless propagator on the 3-torus in terms of an integral involving the jacobi - theta function @xmath389 @xcite : [ eq119 ] g()g(,)=1r__0^s  .
although this expression gives the full compact space propagator , we are mainly interested in the leading order modifications introduced by compactification when the probe is moving closer to the central stack than to any of its images ( orbits around any of the brane images are equivalent due to the symmetries of the hypercubic lattice ) . in @xcite ,
an expansion of this propagator about the origin was computed using ewald s method for calculating potentials in hypercubic lattices in the context of solid state physics .
the leading corrections to the central @xmath390 potential are given by : [ eq120 ] g()=-14r - r^26v_-c_sv_^13-a_4h_4()-a_6h_6()-  the values for the coefficients of the various terms were found numerically to be @xmath391 , @xmath392 and @xmath393 .
the functions @xmath394 and @xmath395 are harmonic functions of order @xmath396 and @xmath397 , respectively , and the series continues with harmonic terms of higher ( even ) orders . the leading harmonic correction to the propagator
is given by : [ eq121 ] h_4()=1r_^5  .
recalling that the function @xmath398 contributes to the interbrane potential through the graviton - dilaton and rr - form interactions , we may write this potential for @xmath399 , to order @xmath400 and discarding constant terms , as [ eq122 ] v()&=&-2q_6  , where we defined the function [ eq123 ] f(y^ir )  .
we conclude that , to this order , the interbrane potential is modified by two terms .
both terms correspond to repulsive contributions to the potential arising from the overall attraction of the image branes and their coefficients are suppressed by a factor @xmath401 , so that their effects are negligible for large transverse volume .
jellium " term , of order @xmath402 and whose name arises in the context of solid state physics , gives an isotropic contribution while the  asymmetry " term , of order @xmath400 , gives an anisotropic contribution to the potential
.    we will be more interested in the effects of the asymmetry term , as its breaks the rotational symmetry of the transverse space and leads to the generation of orbital angular momentum .
the effects of the jellium term will , however , be more significant close to the central stack .
it is easy to show that , for large transverse volume , this term will , to leading order , make the probe brane s orbits precess .
this is due to its repulsive nature and , for small eccentricity orbits , the deficit angle is approximately given by @xmath403 .
this effect will , however , be suppressed at late times as the probe s angular momentum is redshifted away by hubble expansion .
the inclusion of the asymmetry term makes the analysis quite difficult , as the potential is no longer central and depends on both angular coordinates in the transverse space .
this may lead , in particular , to non - planar orbits of the probe brane .
however , if we set @xmath404 initially , the probe will feel no force along this direction and its trajectory will be confined to the @xmath405 plane as before .
we will focus on this particular case , bearing in mind that in general non - planar trajectories may arise .    within these assumptions ,
the problem reduces , as before , to the evolution of the canonically normalized complex scalar field @xmath85 , defined in eq .
( [ eq17 ] ) , in an expanding universe .
its potential can be written as [ eq124 ] v()=-||-_j||^2-_a(_r^4+_i^4 - 6_r^2_i^2 )  , where @xmath406 and @xmath407 are , respectively , the real and imaginary parts of the field , @xmath408 was defined in eq .
( [ eq27 ] ) and [ eq125 ] _
j  , _ a=32a_4t_6v_3q_6r_^5 give the strength of the jellium and asymmetry terms .
this potential is no longer invariant under the global u(1 ) symmetry of the non - compact case , signaling that the probe s angular momentum is no longer conserved .
this is quite similar to the u(1)-violating potential governing the evolution of the scalar field in the affleck - dine mechanism @xcite , where non - conservation of angular momentum plays a crucial role in generating the u(1 ) baryon number asymmetry in our universe .
the potential remains , however , invariant under the hypercubic group symmetries @xmath409 and @xmath410 . in terms of the polar angle @xmath101 ,
this means that @xmath411 rotations as well as reflections about the @xmath412 axis are preserved by the asymmetry term , as one can easily conclude by writing eq .
( [ eq123 ] ) for @xmath80 as [ eq126 ] f()=1 - 2 ^ 2(2 )  .    hence , we only need to consider initial conditions such that @xmath413 or equivalently @xmath414 .
the u(1 ) symmetry violation precludes a complete analytical description of the orbits .
however , it is not difficult to obtain a qualitative insight on the main features of the probe s trajectories . if the probe is placed significantly far from the central stack ( but still closer to it than to any of the image branes ) , the asymmetry term will not be negligible and some angular momentum will be created or destroyed . if the probe has no angular momentum initially , as we would expect immediately after inflation , then it will necessarily acquire some and be placed in an orbit around the central stack , instead of just falling towards it along the radial direction . as both the jellium and the asymmetry interactions are repulsive , it is possible that this trajectory does not remain bound to the central stack , becoming connected to one or more branes in the hypercubic lattice .
hubble expansion will , however , redshift the probe s energy and angular momentum , so that one expects the probe to become bound to only one brane stack at late times .
as the orbital radius decreases , the u(1)-violating term becomes less significant and the particle number associated with @xmath85 should asymptotically become constant , i.e. the angular momentum should vary only due to the universe s expansion .
one also expects the jellium term , as well as the asymmetry term , to induce some precession of the orbital axis , which should stop at late times when these terms become negligible .
due to the hypercubic symmetries , along @xmath415 the force acting on the probe will be in the radial direction , as @xmath416 , and no angular momentum will be created in these particular directions . for @xmath417 , angular momentum creation
should , however , be a generic feature .    to have a better understanding of how the initial conditions affect the amount of angular momentum created or destroyed , we have computed the force acting on the probe brane and the associated torque , given by : [ eq127 ] = -v =- v_r - v_i  , + = ||=-_rv_i+_iv_r  , where we have considered the field @xmath85 as a vector in the @xmath418 plane .
recall that the torque gives the variation of the angular momentum , which in terms of the particle number density @xmath27 defined in eq .
( [ eq20 ] ) can be written as : [ eq128 ] dndt+3hn=  .
we have plotted the isotorque contours for some particular parameter values in figure 7 .
we have restricted the position of the probe to the region @xmath419 , where the expansion of the potential about @xmath420 is sufficiently accurate and the probe s motion is bound to the central stack ( near the limits of this region , the orbits may not be bound to the central stack , as the actual region where the force points towards the central stack has a more complicated shape ) .    ) plane .
coordinates are normalized to the size of the transverse torus , with the value @xmath421 in this example , in units of the string length .
we considered @xmath158 and @xmath422 in these units .
note that the torque changes sign when crossing one of the @xmath423 axis , with integer @xmath27 , being negative for @xmath417 .
also , the absolute value of the torque increases as one moves away from the origin.,title="fig : " ] [ fig : figure_7 ]    observing this figure we see that the torque is maximized close to the boundary of the region mentioned above ( darker areas for @xmath414 ) . also , the longer the probe s trajectory remains close to this boundary the more angular momentum it is likely to gain .
thus , maximum angular momentum creation should occur if the probe is placed below the centre of the hypercubic cell .
the initial value of the hubble parameter @xmath107 will also affect the amount of angular momentum created .
if it is of the same order of magnitude as the orbital frequency , one expects the probe to be driven into lower torque regions within a few periods .
lower values of the hubble parameter should allow the probe s angular momentum to oscillate significantly as the probe moves through alternate regions of positive and negative torque , but no net angular momentum will be gained after completing the first few orbits .
the asymptotic value of the comoving particle number @xmath424 will , hence , be determined by the initial path of the probe .
we have simulated the evolution of the probe brane numerically to illustrate this discussion .
the initial values for the real and imaginary parts of the field were defined as follows : [ eq129 ] _ r , i(t_0)=_r , ir _  , _
r , i=0  , with @xmath425 , giving @xmath417 , and @xmath426 , which ensures the probe is placed initially within the region where it is bound to the central stack .
we have also computed the eccentricity of the orbit , as defined in eq .
( [ eq29 ] ) , with the probe s energy density including both the jellium and asymmetry terms .
this quantity should only become meaningful at late times , when the latter terms become negligible , but its evolution tracks the creation of angular momentum along the probe s trajectory .
figure 8 shows an example of the results obtained for the probe s orbit . in this case
, the probe is placed below the centre of the hypercubic cell and , as expected , a significant angular momentum is produced during its motion .
one observes , as expected , a small precession of the orbital axis and the decay of the orbit into regions of lower torque , which drives the orbital eccentricity towards a constant value .
initial conditions were chosen in this example so that the inverse hubble parameter is of the same order of magnitude as the orbital period and , hence , the orbit is quickly stabilized .
the eccentricity tends in this example to a constant value of @xmath427 .
notice the similarities with the example illustrated in figure 1 , the main difference residing in the fact that , in this case , the probe has initially no angular momentum .    ,
@xmath421 and @xmath428 .
the probe is initially stationary at the point @xmath429 for @xmath430 and we consider a matter - dominated universe with @xmath131 . the asymptotic value obtained for
the eccentricity is @xmath427.,title="fig : " ] [ fig : figure_8 ]    in figure 9 we plot the asymptotic values of the eccentricity obtained numerically for different initial positions in the plane @xmath418 , considering two distinct values of the initial hubble parameter in a matter - dominated universe . in both cases , one observes that smaller eccentricity orbits are obtained near the larger torque regions below the centre of the hypercubic cell , as expected .
we see that , in general , some angular momentum is produced , although the majority of initial conditions leads to highly eccentric orbits , with @xmath431 .
we also conclude that the probe s orbit is generically less eccentric if its motion starts at later times , so that the initial hubble damping of the asymmetry term is smaller . as one may observe in figure 9 , in particular in plot ( b ) , the combination of the asymmetry term and hubble expansion makes the final eccentricity vary in a non - trivial way with the initial position of the probe , namely near @xmath432 , where the orbits may in some cases become unbound .
note that these results correspond to particular values of the parameters of the problem and that the values of the asymptotic eccentricity are highly dependent on these parameters , namely the relevant distance scales @xmath79 and @xmath380 .
, with @xmath433 ( with similar results for @xmath434 ) , assuming @xmath158 , @xmath421 and @xmath428 .
both plots were obtained for a matter - dominated universe with @xmath131 and the motion of the probe starts at ( a ) @xmath435 and ( b ) @xmath436 .
all parameters are given in units of the string length .
contours were plotted for @xmath437 , with @xmath438 , and darker regions correspond to lower eccentricities.,title="fig : " ] [ fig : figure_9 ]    these results show that a significant amount of angular momentum can be created by u(1)-violating terms in the probe s potential arising from compactification of the transverse dimensions , which is a key ingredient for the development of the parametric resonance we have studied in this work .
the simple 3-torus compactification thus ilustrates how low eccentricity orbits can be created for a range of initial conditions , even if the probe has initially no angular momentum .
potential applications of the resonant particle production mechanism to cosmology will be largely determined by the system s final configuration . as we have discussed before ,
the @xmath108 branonium we have analyzed in this work is unstable , as the probe ( anti)brane loses energy and angular momentum due to hubble expansion . at some stage ,
the distance between the probe and the central stack of branes will consequently become of the order of the string length ( although our approximations will break down long before this happens ) .
when this occurs , the scalar mode associated with open strings stretching between the probe and the source branes becomes tachyonic , signaling an instability of the system @xcite .
condensation of this tachyonic mode will then lead to the annihilation of the probe brane with one of the branes in the central stack .
the non - supersymmetric nature of the @xmath108 configuration , which is a general property of all systems mixing brane and antibrane states , is behind this instability .
namely , it makes the system evolve into a @xmath179 bps state with @xmath439 parallel @xmath1 branes , with the same total charge as the original configuration but with a lower energy .
interesting scenarios may arise if some additional mechanism stabilizes the probe brane at a finite distance from the central stack , in which case the produced particles may survive after the resonance mechanism ends .
an attractive possibility for stabilizing multiple brane systems was suggested in @xcite .
when supersymmetry breaking occurs , one or more fields within the supermultiplets of the theory typically acquire masses , their value being set by the energy scale at which susy is broken .
massive gauge potentials will then induce short - range interactions which , if repulsive , may balance the gravitational atraction between different branes .
notice that the presence of a probe antibrane itself breaks susy but that this does not induce bulk field masses , so that other sources of susy breaking need to be considered in this case .
suppose , for example , that the @xmath440 rr - form and consequently its magnetic dual form @xmath441 become massive after susy breaking , the same happening to the dilaton field , @xmath13 .
both the rr - form and dilaton - mediated interactions will then be described by yukawa potentials whose range is determined by the masses @xmath442 and @xmath443 , respectively . if @xmath444 , the dilaton - mediated interaction will be exponentially suppressed for interbrane distances of order @xmath445 .
assuming @xmath446 , we may write the relevant terms in the potential for @xmath447 as    [ eq130 ] v(r)=m^4(-1r )  , where @xmath448 gives the overall constant factor . the constant @xmath345 parametrizes the effective strength of the rr - interaction compared to the gravitational atraction and one expects @xmath449 as the dilaton - mediated interaction is negligible . if , as we have assumed so far , the probe is a @xmath0-brane , the rr - interaction , although short - ranged , is attractive and can not balance the @xmath390 gravitational part . if , however , the probe is a @xmath1-brane , we have @xmath450 and stabilization of the system may be possible . in this case , for @xmath451 , the rr - interaction is exponentially suppressed and the probe s motion is governed only by the gravitational term
. the evolution of the probe brane will , in this case , be very similar to that of a probe antibrane in a massless rr - form background which we have studied in this work , the potential having a smaller strength than in the latter case .
this difference does not modify the qualitative features of the probe s trajectory and all the results we have derived in this work remain valid using an effective central stack charge @xmath452 .
in particular , one expects the parametric resonance to develop in the probe s world - volume under the same conditions as before and the effects of transverse space compactification to provide the necessary source of angular momentum for similar initial conditions
.    the fate of a probe @xmath1-brane will , however , be quite different than that of a probe @xmath0-brane .
as the probe s orbit decays due to the universe s expansion , the interbrane distance will eventually become of order @xmath445 . at these distances ,
the repulsive rr - interaction becomes relevant to the probe s motion and stabilization may be possible .
the potential given in eq .
( [ eq130 ] ) has a local minimum formally given by : [ eq131 ] r_0=-m_rr^-1  , where @xmath453 is the lambert w - function , defined as the inverse of @xmath454 .
after its angular momentum has been redshifted away by hubble expansion , one expects the probe to settle at this local minimum .
this may not be the absolute minimum of the full potential , which includes the dilaton - mediated interaction and other supergravity and string theory corrections closer to the central stack .
hence , absolute stabilization of the probe can not be guaranteed , but it is reasonable to expect at least a long - lived metastable state .    from the properties of the lambert w - function ,
one concludes that @xmath455 for @xmath456 , being positive for @xmath457 , and that @xmath458 strictly grows with @xmath345 .
hence , as expected , stabilization at large distances from the central stack is only possible if the rr - repulsion is stronger than the gravitational attraction . in this case
, one expects the probe brane to stabilize away from the central stack at a distance @xmath459 .
the scale of supersymmetry breaking will then determine the cosmological implications of the particle production mechanism analyzed in this work .
if susy is broken at very high energies , inducing a large mass for the rr - potential ( but still smaller than the dilaton s mass ) , the probe will stabilize very close to the central stack and no tachyonic modes are present to induce annihilation . the stable probe and the central stack
will then be characterized by a general broken gauge symmetry group @xmath460 , where @xmath461 is the number of branes in the probe , with an exponentially large number of particles charged under the @xmath181 gauge group if the resonance develops before stabilization occurs .
this may be relevant , for example , to the generation of the baryon asymmetry in our universe .
on the other hand , a parametrically small soft susy breaking mass will stabilize the probe at large distances .
the parametric resonance may then be relevant for dark matter production , as particles in the probe s world - volume will necessarily interact weakly with the visible sector if the latter is embedded in the central stack .
as we have shown , these particles may be supermassive ( although parametrically small compared to the string scale ) , so that this scenario could provide an effective mechanism for a non - thermal production of wimpzillas @xcite .
one also expects the parametric resonance to amplify other bosonic modes living in the probe s world - volume , such as yang - mills fields if the probe has more than one brane .
if these mediate baryon number - violating forces , their interactions with central stack fields could then provide a mechanism for baryogenesis as the one discussed in @xcite , in which case the probe would need to be stabilized before nucleosynthesis .
it is also possible for fermionic particles living in the probe s world - volume to be produced in resonance , although pauli blocking makes this process significantly different from the bosonic case , according to the discussion given in @xcite .
thus , if the probe can be stabilized , a resonant production of fermionic dark matter particles may also be achieved through this process .
one also expects the interbrane distance field to oscillate about the local minimum of the potential , possibly generating a second stage of resonant amplification .
this may be induced from coupling to other fields , as in the standard preheating mechanism , or gravitationally , as in the case described in this work .
although this second stage may contribute significantly to the final particle number density in the probe brane s world - volume , we will not discuss them in further detail but rather refer the reader to the discussions given in this work and in the literature @xcite .
we emphasize that all these scenarios depend on the particular embedding of the standard model fields in this setup , an issue for which there is still no complete answer and that is closely related to the nature of susy breaking and its communication to the standard model .
in this work we have analyzed the mechanism of particle production in @xmath108 branonium systems , as well as in @xmath462 systems with susy breaking .
we have shown that , in the limit of large distances and small velocities and for small eccentricity orbits , a parametric resonance will develop , producing scalar particles confined to the probe brane s world - volume .
massless particles are produced in a narrow resonance band centered at the comoving momentum corresponding to one half of the probe s orbital frequency @xmath372 , as if resulting from the direct decay of the particles associated with the interbrane distance field .
the strength of the resonance is given by the exponent @xmath239 defined in eq .
( [ eq72 ] ) .
massive particles are produced in the same way , with energies of order @xmath463 .
these energies are small compared to the string energy scale but , as the latter should be close to the planck scale , very massive particles can be produced by the resonance mechanism , according to eq .
( [ eq84 ] ) . for both massless and massive particles , the associated particle number grows exponentially in time and a large number of particles can be produced by this mechanism , just as in the case of preheating after inflation @xcite .
however , the resonance surely requires the probe to complete a large number of orbits before the particle - modes can be significantly amplified , due to the small value of the exponent @xmath464 .
hubble expansion of the 4-dimensional effective flat frw spacetime makes the probe s orbit decay , redshifting its energy and angular momentum .
this alters the development of the parametric resonance significantly if the probe s motion does not occur at sufficiently late times .
if , however , its motion begins late in the history of the universe , at least compared to typical string times , the resonance will still develop with the associated resonance band being shifted towards higher momentum modes .
each mode will be excited during the finite period it spends inside the resonance band , after which the associated particle number becomes an adiabatic invariant . a significant number of particles , given appproximately by eq .
( [ eq101 ] ) , can be produced if this period is sufficiently long for the resonance to develop , according to the condition @xmath465 , where @xmath20 quantifies the effects of the expansion and was defined in eq .
( [ eq93 ] ) .    we have also concluded that a large number of brane - bound particles can be produced without affecting the probe s orbital motion , as can be seen in eq .
( [ eq112 ] ) .
radiation into bulk closed string modes will also have a negligible effect if the probe moves at sufficiently large distances , as assumed for the validity of our study .
a realistic implementation of this mechanism requires particular initial configurations with large interbrane distances , non - relativistic velocities and almost circular orbits .
we have shown that angular momentum creation may result from the effects of compactification of the directions transverse to the branes , illustrating this for the case of a compact 3-torus , where low eccentricity orbits can be produced for some range of initial conditions .    in this work ,
we have also discussed how , for a probe @xmath1-brane , an interbrane potential may arise from supersymmetry breaking , if e.g. the relevant rr - form and the dilaton , @xmath13 , gain a mass in the process .
stabilization of the probe is then possible at late times and at distances of order @xmath445 , if @xmath466 .
such a stabilization avoids annihilation with the central stack , as necessarily happens for a probe antibrane , making the parametric resonance regime potentially relevant for baryon number or dark matter generation , for example .    in our analysis ,
the parametric resonance arises only for @xmath108 or @xmath462 systems , where the @xmath390 potential allows the probe to move in closed elliptical orbits in a non - expanding universe .
it is crucial for the development of the resonance that the interbrane distance field @xmath100 exhibits an oscillating behavior , although the orbits do not need to close , as we have seen for the case of an expanding universe .
thus , other @xmath4 or @xmath467 branonium systems with @xmath468 , where the probe trajectories are necessarily unbounded from below , can not exhibit such resonant particle production .
additional stabilization potentials may , however , provide the required oscillatory behavior , so that it may be possible to find a resonant particle production mechanism in other cases .
it is nevertheless clear that resonant particle production in branonium systems may play an important role in our universe s evolution and we hope with this work to motivate future exploitation of their properties and applications .    0.2 cm    * acknowledgements *    0.2 cm the authors would like to thank graham ross , martin schvellinger , andre lukas and yang - hui he for useful discussions .
jgr is supported by fct ( portugal ) under the grant sfrh / bd/23036/2005 .
jmr is partially supported by the ec network @xmath469 framework programme research and training network
quest for unification " ( mrtn - ct-2004 - 503369 ) .
this work was partially supported by the eu fp6 marie curie research and training network
 universenet " ( mrtn - ct-2006 - 035863 ) .    0.5 cm
the particle number produced by the parametric resonance can be estimated by computing the integral in eq .
( [ eq100 ] ) .
this gives the physical particle number apart from constant factors , subdominant non - resonant contributions and also an @xmath369 redshift factor due to hubble expansion .
we have ,          in order to compute this function , let us write the exponent @xmath279 in the form : [ eqa2 ] _
k(z)=12  , where [ eqa3 ] a&=&q_0 ^ 2 ^ 2-a_k0 ^ 2 ^ 2  , + b&=&2(q_0 ^ 2+a_k0(a_k0 - 1 ) )  , + c&=&q_0 ^ 2-(a_k0 - 1)^2  . for all modes of interest , @xmath470 , while we have @xmath471 for type ii modes and @xmath472 for type iii modes .
the sign of @xmath473 will not affect our results .
it is also useful to define @xmath474 , which is explicitly given by [ eqa4 ] d^2=4q_0 ^ 2(2a_k0(a_k0 - 1)+(a_k0 - 1)^2 ^ 2+a_k0 ^ 2 ^ 2 )  , and is positive for all type ii and type iii momentum modes . with these considerations ,
we find , for @xmath475 , [ eqa5 ] n_k^res(z)&=12a_z_i^z  .    taking into account that @xmath476 , @xmath477 and also that @xmath478 , we can write the total particle number produced by the resonance for type ii modes as : [ eqa6 ] n_k^res(ii)(z_2)=-14a  . similarly , for type iii modes , we obtain [ eqa7 ] n_k^res(iii)(z_2)=8d^2(-a)^32  .
it is easy to check that eqs .
( [ eqa6 ] ) and ( [ eqa7 ] ) give the same result for the mode in the transition between types ii and iii , at @xmath479 , as expected . recalling that @xmath300
, these expressions can be approximated by those given in eq .
( [ eq101 ] ) .
j.  dai , r.  g.  leigh and j.  polchinski , mod .
lett .   a * 4 * ( 1989 ) 2073 ; j.  polchinski , phys .
lett .   *
75 * , 4724 ( 1995 ) [ arxiv : hep - th/9510017 ] .
m.  berkooz , m.  r.  douglas and r.  g.  leigh , nucl .
b * 480 * , 265 ( 1996 ) [ arxiv : hep - th/9606139 ] .
g. aldazabal _
et al . _ , j.  math .
phys .   * 42 * , 3103 ( 2001 ) [ arxiv : hep - th/0011073 ] ; l.  e.  ibanez , f.  marchesano and r.  rabadan , jhep * 0111 * , 002 ( 2001 ) [ arxiv : hep - th/0105155 ] ; d.  bailin , g.  v.  kraniotis and a.  love , phys .
b * 530 * , 202 ( 2002 ) [ arxiv : hep - th/0108131 ] ; d.  cremades , l.  e.  ibanez and f.  marchesano , jhep * 0207 * , 022 ( 2002 ) [ arxiv : hep - th/0203160 ] ; a.  m.  uranga , class .
quant .  grav .
* 20 * , s373 ( 2003 ) [ arxiv : hep - th/0301032 ] ; c.  kokorelis , jhep * 0209 * ( 2002 ) 029 [ arxiv : hep - th/0205147 ] ; c.  kokorelis , jhep * 0208 * , 036 ( 2002 ) [ arxiv : hep - th/0206108 ] ; n.  ohta and p.  k.  townsend , phys .
b * 418 * , 77 ( 1998 ) [ arxiv : hep - th/9710129 ] .
+ for recent reviews see , e.g. , r.  blumenhagen , m.  cvetic , p.  langacker and g.  shiu , ann .  rev .
nucl .  part .
* 55 * , 71 ( 2005 ) [ arxiv : hep - th/0502005 ] ; r.  blumenhagen , b.  kors , d.  lust and s.  stieberger , phys .
* 445 * , 1 ( 2007 ) [ arxiv : hep - th/0610327 ] .
m. r. douglas , in _ cargese 1997 , strings , branes and dualities _ [ arxiv : hep - th/9512077 ] ; a.  a.  tseytlin , nucl .
b * 475 * , 149 ( 1996 ) [ arxiv : hep - th/9604035 ] ; a.  a.  tseytlin , class .
* 14 * , 2085 ( 1997 ) [ arxiv : hep - th/9702163 ] ; m.  s.  costa , nucl .
b * 490 * , 202 ( 1997 ) [ arxiv : hep - th/9609181 ] ; z.  guralnik and s.  ramgoolam , nucl .
b * 499 * , 241 ( 1997 ) [ arxiv : hep - th/9702099 ] .
n.  arkani - hamed , s.  dimopoulos and g.  r.  dvali , phys .
d * 59 * , 086004 ( 1999 ) [ arxiv : hep - ph/9807344 ] ; g.  r.  dvali and s.  h.  h.  tye , phys .
b * 450 * , 72 ( 1999 ) [ arxiv : hep - ph/9812483 ] ; n.  arkani - hamed , _
et al . _ ,
b * 567 * , 189 ( 2000 ) [ arxiv : hep - ph/9903224 ] ; p.
binetruy , c.  deffayet and d.  langlois , nucl .
b * 565 * , 269 ( 2000 ) [ arxiv : hep - th/9905012 ] ; j.  m.  cline , c.  grojean and g.  servant , phys .
rev .  lett .
* 83 * , 4245 ( 1999 ) [ arxiv : hep - ph/9906523 ] ; c.  csaki , _ et al . _ , phys .  lett .
b * 462 * ( 1999 ) 34 [ arxiv : hep - ph/9906513 ] ; s.  s.  gubser , phys .
d * 63 * , 084017 ( 2001 ) [ arxiv : hep - th/9912001 ] ; a.  hebecker and j.  march - russell , nucl .
b * 608 * , 375 ( 2001 ) [ arxiv : hep - ph/0103214 ] ; h.  a.  bridgman , k.  a.  malik and d.  wands , phys .
d * 65 * , 043502 ( 2002 ) [ arxiv : astro - ph/0107245 ] ; a.  kehagias and e.  kiritsis , jhep * 9911 * , 022 ( 1999 ) [ arxiv : hep - th/9910174 ] ; p.  kraus , jhep * 9912 * ( 1999 ) 011 [ arxiv : hep - th/9910149 ] .
+ for reviews of brane cosmology from the effective field theory viewpoint see , eg , d.  langlois , prog .
suppl .   * 148 * , 181 ( 2003 ) [ arxiv : hep - th/0209261 ] ; r.  maartens , living rev .
* 7 * , 7 ( 2004 ) [ arxiv : gr - qc/0312059 ] .
c.  p.  burgess _ et al .
_ , jhep * 0107 * , 047 ( 2001 ) [ arxiv : hep - th/0105204 ] .
c. p. burgess _
et al . _ ,
jhep * 0203 * , 052 ( 2002 ) [ arxiv : hep - th/0111025 ] . g.  shiu and s.  h.  h.  tye , phys
.  lett .
b * 516 * , 421 ( 2001 ) [ arxiv : hep - th/0106274 ] ; n.  t.  jones , h.  stoica and s.  h.  h.  tye , jhep * 0207 * , 051 ( 2002 ) [ arxiv : hep - th/0203163 ] ; j.  h.  brodie and d.  a.  easson , jcap * 0312 * , 004 ( 2003 ) [ arxiv : hep - th/0301138 ] ; h.  firouzjahi and s.  h.  h.  tye , phys .
b * 584 * , 147 ( 2004 ) [ arxiv : hep - th/0312020 ] .
+ for reviews of brane inflation in string theory see , eg , e.  kiritsis , fortsch .
phys .   * 52 * , 200 ( 2004 ) [ phys .
* 421 * , 105 ( 2005  errat,429,121 - 122.2006 ) ] [ arxiv : hep - th/0310001 ] ; s.  h.  henry tye , arxiv : hep - th/0610221 .
j.  h.  traschen and r.  h.  brandenberger , phys .
d * 42 * , 2491 ( 1990 ) ; l.  kofman , a.  d.  linde and a.  a.  starobinsky , phys .
lett .   * 73 * , 3195 ( 1994 ) [ arxiv : hep - th/9405187 ] ; y.  shtanov , j.  h.  traschen and r.  h.  brandenberger , phys .
d * 51 * , 5438 ( 1995 ) [ arxiv : hep - ph/9407247 ] ; d.  boyanovsky , h.  j.  de vega , r.  holman , d.  s.  lee and a.  singh , phys .
d * 51 * , 4419 ( 1995 ) [ arxiv : hep - ph/9408214 ] ; m.  yoshimura , prog .
phys .   *
94 * , 873 ( 1995 ) [ arxiv : hep - th/9506176 ] ; d.  i.  kaiser , phys .
d * 53 * , 1776 ( 1996 ) [ arxiv : astro - ph/9507108 ] ; s.  y.  khlebnikov and i.  i.  tkachev , phys .  rev .
lett .   * 77 * , 219 ( 1996 ) [ arxiv : hep - ph/9603378 ] .
l.  kofman , a.  d.  linde and a.  a.  starobinsky , phys .
d * 56 * , 3258 ( 1997 ) [ arxiv : hep - ph/9704452 ] ; p.  b.  greene , l.  kofman , a.  d.  linde and a.  a.  starobinsky , phys .
d * 56 * , 6175 ( 1997 ) [ arxiv : hep - ph/9705347 ] ; g.  n.  felder , l.  kofman and a.  d.  linde , phys .
d * 59 * ( 1999 ) 123523 [ arxiv : hep - ph/9812289 ] .
s.  kachru , r.  kallosh , a.  linde , j.  m.  maldacena , l.  p.  mcallister and s.  p.  trivedi , jcap * 0310 * , 013 ( 2003 ) [ arxiv : hep - th/0308055 ] .
s.  shandera , b.  shlaer , h.  stoica and s.  h.  h.  tye , jcap * 0402 * , 013 ( 2004 ) [ arxiv : hep - th/0311207 ] .
a.  sen , jhep * 9808 * , 012 ( 1998 ) [ arxiv : hep - th/9805170 ] ; a. sen , in _
cargese 1999 , progress in string theory and m - theory _ , arxiv : hep - th/9904207 .
e. kolb , d. j. h. chung and a. riotto , in _ buenos aires 1998 , trends in theoretical physics 2 _ , arxiv : hep - ph/9810361 .
e. w. kolb , a. linde and a. riotto , phys .
lett .   * 77 * , 4290 ( 1996 ) [ arxiv : hep - ph/9606260 ] . p. b. greene and l. kofman , phys
.  rev .
d * 62 * , 123516 ( 2000 ) [ arxiv : hep - ph/0003018 ] .