Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: soon after the invention of inflationary cosmology it was realized that inflation may divide our universe into many exponentially large domains corresponding to different metastable vacuum states @xcite . the total number of such vacuum states in string theory can be enormously large @xcite . a combination of these two facts with the kklt mechanism of vacuum stabilization @xcite gave rise to what is now called the string landscape scenario @xcite . the purpose of this paper is to clarify some features of the landscape . in particular , one may wonder what are the properties of the domain walls separating different vacua , how these domains may form , and how large will be the fraction of the volume of the universe inside different domains .    to this end , we will construct a broad class of static bps domain wall solutions interpolating between different ads or minkowski vacua with unbroken supersymmetry . quite generally , these bps domain walls may interpolate between various minima , maxima , or saddle points of the scalar field potential . bps domain wall solutions in 4-dimensional @xmath3 supergravity were first studied in @xcite . the methods used for that analysis have been later refined and generalized to a variety of models , for different number of dimensions and supersymmetries , see e.g. @xcite . our present aim is to further extend these methods to make them suitable for investigation of the string landscape scenario . following the logic of the recent developments , we will begin with an investigation of the supersymmetric vacua arising from flux compactifications prior to the kklt uplifting , and later introduce the more physically relevant metastable de sitter vacua . we start in section 2 with an overview of some known features of domain walls in effective d=4 , n=1 supergravity . the new results here include the derivation of bps domain wall gradient flow equations in generic models with f - term potentials with arbitrary khler potential , superpotential and number of moduli . we rewrite the gravitational action , including the boundary @xmath4-terms in a bps form : when the flow equations are satisfied , the integrand of the action defines the bps domain wall tension . afterward , in the spirit of the analogous work in five dimensions in the context of the ads / cft correspondence , where domain walls are viewed as dual to renormalization group flows @xcite , we present the bps flow equations in the form of rg flows and provide the relevant `` c - theorem '' .    in section 3 we give examples of bps domain walls where the potential has more than one critical point but there is no barrier between them . our examples include an ads minimum and an asymptotically minkowski vacuum . other examples have critical points at finite distance in the moduli space : an ads saddle point and an ads minimum , or an ads maximum and an ads saddle point . a particularly interesting case of flux vacua area codes in a calabi  yau hypersurface required the study of the flow of four moduli . they flow from one critical point of the potential , the landau  ginsburg fixed point , towards the second critical point with larger absolute value of the cosmological constant near the conifold . in section 4 we study a class of models in stringy landscape with a barrier between critical points . we have examples of a local minkowski or ads minimum and another absolute ads minimum of larger absolute value of the cosmological constant . in case of ads - ads wall , we display the domain walls with critical and super - critical tension . the last case has the superpotential crossing zero between two ads vacua . the results of our investigation in sections 2 - 4 directly apply to the transitions between supersymmetric ads and minkowski vacua . however , the most interesting parts of the string theory landscape correspond to ds vacua . these vacua can be obtained from supersymmetric ads or minkowski vacua by uplifting , which can be achieved , e.g. , by adding some anti - d3 branes @xcite . this procedure breaks supersymmetry and may destabilize some of the vacua , in particular those which corresponded to maxima and saddle points . as a result , some of the bps wall solutions may disappear after the uplifting , whereas some others will remain . however , instead of describing the static bps domain walls , they will describe domain walls moving with acceleration . such domain walls can be formed during the tunneling and bubble formation in the process of decay of metastable vacua , as discussed by coleman and de luccia ( cdl ) in @xcite . we will discuss this process in section 5 .    if the uplifting is relatively small , which may happen when we uplift a supersymmetric minkowski vacuum , or a shallow ads vacuum , as in the version of the kklt scenario developed in ref . @xcite , then the properties of the domain walls do not differ much form the properties of the bps walls prior to the uplifting . in this case we will be talking about near - bps , or near - extremal walls and bubbles . the concept of near - bps walls and bubbles is closely related to the concept of near - extremal black holes . whereas at the extremal limit the bps black holes are absolutely stable and have zero temperature , the near - extremal black holes have a non - vanishing temperature , they can evaporate and eventually approach a stable extremal limit . the physics of near - extremal black holes has been studied intensively over the years : a remarkable agreement with the physics of d - branes was established under condition of a small deviation from extremality @xcite .    in our new endeavor we will try to understand the rate of tunneling into a collapsing vacuum with a negative cosmological constant from uplifted slightly non - bps vacua using the proximity of the large size bubbles to the bps domain walls whose tension we can evaluate exactly . in the bps limit the corresponding bubbles have an infinite radius , but the near - bps ones have a large but finite radius @xcite . the description of these bubbles can be performed using the deviation from the bps limit as an expansion parameter . this simplifies investigation of bubble formation and vacuum decay in such models . in section 5 we will find , in particular , that the vacuum decay rate of the uplifted supersymmetric minkowski vacuum does not depend on the details of the model and is given by the universal expression @xmath5 , where @xmath6 is the entropy of the ds state . this decay rate is much higher than the typical rate of decay of a ds vacuum to a minkowski vacuum , which in many cases does not differ much from @xmath7 @xcite . in application to our vacuum with @xmath8 , we are talking about the difference in the decay rate by a factor @xmath9 . the situation is especially interesting in a more general case when we uplift the models with several different ads minima of the moduli potential . the ads vacuum with a smaller value of the volume modulus experiences a bigger uplifting . the scale of supersymmetry breaking in the uplifted vacuum is related to the magnitude of the uplifting . as a result , the decay of the uplifted ads vacuum to the collapsing space with a negative vacuum energy usually occurs much faster than the typical ds rate of the uplifted minkowski vacuum @xmath10 . under some conditions , the decay rate can be quite large , with the suppression factor @xmath11 being replaced by @xmath12 . here @xmath13 and @xmath14 is the depth of the ads minimum prior to the uplifting , which may be a hundred orders of magnitude greater than @xmath8 . in the kklt - type models , the depth of the ads vacuum prior to the uplifting is of the same order as @xmath15 , where @xmath16 is the gravitino mass after the uplifting @xcite . therefore in the class of models studied in section 5 , the vacuum decay rate due to the tunneling to a collapsing universe is related to the gravitino mass , @xmath17 . after the tunneling , the universe within each bubble with negative vacuum energy collapses within a microscopically small time @xmath18 . the condition that the lifetime of the uplifted vacuum is greater than the age of our part of the universe may serve as a vacuum superselection rule . an important feature of the decay of ds space to a minkowski vacuum or to a collapsing universe with a negative vacuum energy density is that the part of space that experiences such a decay never returns back to ds space . this is quite different from the back - and - forth transitions between various ds vacua , which may keep the system in a kind of thermal equilibrium . in this respect , the channels of irreversible decay of an uplifted ds vacuum to a minkowski vacuum or to a collapsing universe serve as sinks for the probability flow in the landscape . in section 6 we will show that the existence of such sinks has important implications for the calculations of the probabilities in the landscape . in particular , if one ignores the different rate of expansion of the universe in different vacua , i.e. uses a comoving volume distribution , then in the absence of the sinks , the fraction of the comoving volume of the universe in the state with a given vacuum energy density @xmath19 asymptotically approaches @xmath20 . this distribution has a sharp maximum at the smallest possible @xmath21 . meanwhile if the vacuum with the smallest possible @xmath21 is not the ground state but a sink , and the probability of tunneling to the sink is sufficiently large , then the fraction of the volume in the vacuum with the smallest @xmath21 will be exponentially suppressed , and the maximum of the probability distribution will be at larger values of @xmath21 . walls and bubbles considered in our paper interpolate between domains with different values of the scalar fields , but with the same values of fluxes . in section 7 we will briefly review our results and discuss a possibility of their generalization for the transitions between vacua with different fluxes .    1 cm in this paper we are going to analyze the domain wall solutions interpolating between different supersymmetric vacua of the landscape of string theory flux compactifications to 4 dimensions . we will focus on the case of @xmath3 supergravity . bps domain wall solutions with 1/2 of unbroken @xmath3 supersymmetry in the minimal 4-dimensional supergravity have been discussed in @xcite . the equations presented there refer to the case of a single scalar field . here we are going to extend this analysis to a more general setup allowing for the possibility of many moduli fields with generic metric in the moduli space . of course , in any sufficiently small patch , there is a possible reparametrization of the moduli space so that there is a single scalar field supporting the domain wall . in a rather generic setup , however , this is not possible globally , and therefore the equations we will propose in the following are more useful in this respect . the analysis presented here is mainly based on the approach presented in @xcite for the domain wall solutions of the minimal 5-dimensional supergravity . we briefly recall the main ingredients needed in order to discuss supersymmetric bosonic configurations of n=1 supergravity coupled to chiral multiplets . the bosonic sector of the action is @xmath22 where the potential can be expressed in terms of a superpotential @xmath23 as @xmath24 or , in terms of the covariantly holomorphic superpotential @xmath25 , as @xmath26 the covariant derivative is defined as @xmath27 on holomorphic quantities ( like @xmath28 ) and as @xmath29 on covariantly holomorphic quantities ( like @xmath30 ) . we will not consider vector multiplets and possible d - term extensions here . the relevant supersymmetry transformations ( of the gravitino and modulini ) are @xmath31 \delta\chi^{i}_{l } & = & \frac12 \gamma^{\mu }   \partial_{\mu } z^{i}\epsilon_{r } - \frac12 g^{i\bar\jmath}{\rm e}^{k/2 } \overline{d_{j}w}\epsilon_{l } , \label{deltachi } \end{aligned}\ ] ] where @xmath32 is the @xmath33 khler connection .      in the following we will consider flat space - like domain wall solutions of the supersymmetry conditions and of the supergravity equations of motion . for this reason we are going to use the following metric ansatz @xmath34 given the domain wall structure , we split the indices @xmath35 , @xmath36 . for the metric ( [ typei ] ) there is only one non - trivial component of the spin connection @xmath37 where @xmath38 , leading also to the ricci scalar ( which will become useful later ) @xmath39 . \label{scalari}\ ] ] plugging ( [ typei ] ) into the supersymmetry equations ( [ deltapsi])([deltachi ] ) and assuming that also the scalars @xmath40 depend only on the radial coordinate @xmath41 , one obtains @xmath42 \delta \psi_{i l } & = & \frac{a^\prime}{2 } \gamma_{\underline i}\gamma_{\underline r } \epsilon_l + \frac12 { \rm e}^{k/2 } w \gamma_i \epsilon_r \label{delpsiii } \\[2 mm ] \delta \chi_l & = & \frac12 \gamma^r z^{i\prime}(r ) \epsilon_r - \frac12 { \rm e}^{k/2 } g^{i\bar \jmath } \overline{d_j w } \epsilon_l , \label{delchii } \end{aligned}\ ] ] where @xmath43 was assumed . using also @xmath44 , one can start solving the various equations . from ( [ delpsiii ] ) , splitting the complex superpotential into norm and phase @xmath45 we get the projector condition @xmath46 from which the bps equation follows @xmath47 the condition ( [ proj0 ] ) can be assembled into a real projector @xmath48 ( satisfying @xmath49 ) acting on the full spinor as @xmath50 : @xmath51 . \label{projectori}\ ] ]    the radial equation on the gravitino gives a differential equation on the spinor parameter @xmath52 here , the projector condition was used . moreover , one can use the projector condition inside ( [ diffeps ] ) , or take the charge conjugate ( using the majorana condition ) . the two resulting equations are @xmath53 \epsilon_r^{\prime } & = & \frac{i}{2}a_r^b \epsilon_r \pm \frac12 { \rm e}^{k/2 } |w| \epsilon_r , \label{majorana } \end{aligned}\ ] ] whose consistency implies @xmath54    the solution is given in terms of an arbitrary constant left - handed spinor parameter @xmath55 ( hence the 1/2 supersymmetry preserved ) : @xmath56 the bps equation ( [ bpsi1 ] ) has been used to express the integration over the superpotential in terms of the warp factor and ( [ bpsextra ] ) for the integration over @xmath57 . in the last formula we have used the definition of the charged conjugate spinor @xmath58 that follows from choosing the charge conjugation matrix @xmath59 , with @xmath60 . in this way the spinor @xmath61 satisfies the relation ( [ proj0 ] ) .    finally , using again the projector condition in the modulini equations , one finds the bps equation for the scalar fields @xmath62 we can now also remove the dependence on the phase @xmath63 , by using the holomorphicity of the superpotential : @xmath64 this implies that the bps equations for the scalars can take the form : @xmath65 summarizing , the supersymmetry conditions that have to be satisfied in order to specify the domain wall solutions are @xmath66      z^{i \prime } & = & \mp 2 g^{i\bar \jmath }       \partial_{\bar \jmath } |{\cal z}| . \end{array } $ } \label{summary}\ ] ] the other consistency condition we derived from the solution of the supersymmetry equations , namely ( [ bpsextra ] ) , is identically satisfied . this can be seen by using the fact that @xmath63 is the phase of the superpotential , and by equating with ( [ bpsextra ] ) the condition coming from its definition , @xmath67 . the resulting equation is an identity upon using ( [ summary ] ) . domain walls in @xmath3 supergravity were studied before in @xcite in a more restrictive setting , so we compare our results . first there is a change in the signature . this implies that one should change the sign for every upper index and put an @xmath68 for each gamma matrix in the supersymmetry equations . this can be re - absorbed in sending @xmath69 . then there is a difference in the metric . on the solution they have @xmath70 , but also their radial coordinate @xmath71 is related to ours by @xmath72 . also , the decomposition of the superpotential in norm and phase implies that @xmath73 . altogether the equations in @xcite become @xmath74 t^{\prime } & = & -\zeta { \rm e}^{k/2+i \theta } g^{t\overline t } \overline{d_t w } , \label{cv2 } \\[2 mm ] \theta^{\prime } & = & \frac{i}{2 } \left(t^{\prime }   \partial_{t } k - \overline t^{\prime }   \partial_{\bar t } k\right ) . \label{cv3 } \end{aligned}\ ] ] it is easy to see that these equations are equivalent to the ones presented above , when restricted to a single scalar field flowing . we find it useful to present the domain wall equations in the form of renormalization group equations . for this purpose we will describe the flow not as a function of @xmath41 coordinate but as a function of the warp factor . we multiply eq . ( [ bpsi2a ] ) on @xmath75 and find , using the eq . ( [ bpsi1 ] ) @xmath76 and @xmath77 we will find below that in certain cases the domain wall bps equations in the form ( [ summary ] ) are useful . sometimes we will use equations in the form ( [ renorm1 ] ) , particularly in the case of many moduli participating in the flow and @xmath28 not crossing zero . a careful inspection of the bps conditions ( [ summary ] ) shows that we can associate to the domain wall solutions a monotonic function , which therefore specifies the direction of the flow . as it is known since the analysis of the 5-dimensional domain walls and the interpretation of such solutions as renormalization group flows in the context of the ads / cft correspondence , there is a c - function that is given by the first derivative of the logarithm of the warp factor @xmath78 @xcite . by using the bps equations we easily see that @xmath79 where the last inequality follows from the properties of the scalar metric .    there is , however , another interesting quantity that in most cases works as a c - function . this is the square of the covariantly holomorphic superpotential , which is related to the square of the mass of the gravitino at a generic point of the moduli space @xmath80 we can now look at the variations of @xmath81 along the flow . in order to compute it , we actually evaluate @xmath82 and use the bps equation ( [ bpsi2 ] ) . after some easy algebra we get that @xmath83 due to the positive definiteness of the metric of the scalar fields , this expression is clearly either always increasing , or always decreasing along the flow , depending on the sign of the bps equation , unless the flow crosses points where the superpotential vanishes . for these special cases the bps equations change sign as we cross the zeros of the superpotential and therefore @xmath81 is first decreasing until it reaches zero and then starts increasing again .      the supersymmetric domain walls studied above are static 1/2 bps configurations of @xmath84 supergravity . in an analogous situation for static 1/2 bps cosmic strings in n=1 supergravity @xcite it was possible to rewrite the gravitational action , including the boundary terms , in the bps form . the action therefore was given by some squares of the bps conditions and total derivative terms . a saturation of the bps bound is equivalent to solving first order differential equations for unbroken supersymmetry . the on shell action therefore is equal , up to a sign , to the total energy of the string . therefore the tension of the string , the energy per unit area , can easily be derived from the bps action . the metric of the cosmic string is @xmath85 and the relation between the @xmath86 on shell and a tension @xmath87 is @xmath88 where @xmath89 and @xmath90 is a fi term in the action .    for flat domain walls the metric is @xmath91 . therefore one would expect the following relation between the action and the tension ( the energy per unit of area ) @xmath92 in the case of cosmic strings @xmath93 and @xmath71 have flat metric , therefore the volume of the area is just @xmath94 . however , for our domain walls we have to distinguish the physical volume from coordinate volume and we have to ignore ( divide by ) the physical volume which is @xmath95 .    in @xcite the tension of the bps domain walls was calculated using a rather complicated nester construction . the formula for the tension / energy density of the wall between two ads critical points is given for two cases : in one case the superpotential does not vanish between critical points , in the other case it vanishes at some intermediate point :    * 2 ads vacua , @xmath96 in between @xmath97 * 2 ads vacua , @xmath98 in between @xmath99    here we show that for domain walls we can put the supergravity action in the standard bps form , namely as the integral of the square of the bps conditions plus boundary terms . we start with the gravitational action supplemented by the gaussian curvature at the boundaries . @xmath100 here @xmath101 is the metric at the boundaries . in our case the boundaries are at @xmath102 surface . therefore @xmath103 and @xmath104 where @xmath105 and @xmath106 and @xmath107 . one starts by plugging ( [ scalari ] ) in the einstein term and using that all the quantities depend only on the radial coordinate . then , we try to re - combine the various terms in the form of the bps equations ( [ bpsi1])([bpsi2a ] ) plus some boundary term . the action ( [ action1 ] ) becomes @xmath108\nonumber \\ & & + \int d^3 x \ , ( a^3)'| _ { r=+\infty } - ( a^3)'| _ { r=-\infty } ) \ .\end{aligned}\ ] ] this can now be recast in the form of bps equations squared plus additional boundary terms as @xmath109 \nonumber \\[2 mm ] & \mp & \int d^3 x \ , dr \left[2(a^3)^\prime |{\rm e}^{k/2}w| + a^3 z^{i\prime }   \partial_i ( e^{k/2}w){\rm e}^{-i\theta}+a^3 \bar z^{\bar \imath}\bar   \partial_i ( e^{k/2}\bar w){\rm e}^{i\theta}+ a^3 ( e^{k/2})^{\prime } |w|\right ] \label{bpsform0 } \end{aligned}\ ] ] finally , the second line becomes a total derivative by using ( [ holow ] ) , and the total action reads @xmath110 \nonumber \\[2 mm ] & \mp & \int d^3 x \ , dr \frac{d}{dr}\left[a^3\left(2 |{\cal z}|\right)\right ] , \label{bpsformi } \end{aligned}\ ] ] where we recovered the bps equations in the first line and we have a boundary term contribution in the second line . one can define the energy density functional @xmath111 , since we have a static solution , so that @xmath112 . if the superpotential does not vanish between the critical points we find @xmath113\nonumber\\ & & + \int d^3 x ( a^3 2 |{\cal z}| _ { r=+\infty}- a^3 2 |{\cal z}| _ { r=-\infty } ) . \label{bpsform2 } \end{aligned}\ ] ] in case that @xmath114 we have to make a choice of the sign in bps equations such that the tension is positive , as shown in eq . ( [ bpsform2 ] ) . thus when @xmath115 does not vanish between fixed points and the bps equation of motion are satisfied with the upper sign the energy density of the wall is given by @xmath116 note that @xmath117 . so our answer is exactly as in eqs . ( [ 1 ] ) , in agreement with @xcite . it was explained in @xcite that the bps domain walls with this tension are related to cdl bubbles with infinite radius . now we assume that there is a solution of the bps equations where the superpotential vanishes at some point @xmath118 between the ads critical points . the boundary integral in eq . ( [ bpsformi ] ) must be split into two terms since one part of the flow proceeds with one sign in the bps equations and for the other part , after @xmath119 point , the sign is flipped . therefore the answer for the energy density in this case is @xmath120 noting that @xmath117 , we have derived the correct tension as in eq . ( [ 2 ] ) in agreement with @xcite . this case , according to @xcite is not related to cdl tunneling . we would like to conclude this subsection by the observation that our new derivation of the bps domain wall tension follows the pattern discovered for the bps cosmic string in @xcite . we find the derivation using the bps form of the action significantly more elegant than using the nester construction . in this section we are going to give some examples of supersymmetric domain wall solutions for some of the generic models which capture the most important features of the landscape of flux vacua . we start with the simple scenario with just one khler modulus and all the complex structure moduli frozen : the kklt model @xcite . this model has one supersymmetric ads vacuum at a finite point of the volume modulus and a minkowski vacuum at infinite value of the volume modulus . we will solve the bps equations numerically and display the corresponding domain wall . after this simple example we show a more involved setup where also the complex structure moduli and the axion - dilaton are allowed to flow : the @xmath121 models . these include a wide variety of models ( and also the kklt model ) as a special case . the case which we will solve numerically has two different ads critical points at finite distance in the moduli space and three scalars are flowing between the critical points , two of them proportional to each other . finally we focus on a model where the effect of the khler deformations is neglected , but the form of the complex - structure deformations is more involved . the flux vacua in such models are described by a simple set of `` new attractor equations '' @xcite relating the points in the moduli space where the moduli are stabilized by fluxes . these flux vacua attractors have some resemblance as well as some differences with the black hole attractor equations @xcite . we will focus here on a model of a calabi  yau threefold obtained as a hypersurface in @xmath122 . this model has two different ads critical points @xcite at finite distance in the moduli space . we will find out how four scalars in this model are flowing from the landau  ginsburg ads critical point towards the near conifold region of the moduli space . the kklt model @xcite suggests a way to generate meta - stable de sitter vacua in string theory starting from supersymmetric anti de sitter minima with all moduli stabilized due to fluxes and using non - perturbative effects , such as gaugino condensation and/or d3-instantons . the uplifting of these stabilized ads vacua to ds vacua is possible under certain conditions by the addition of positive energy density sources , like anti - d3 branes . the simplest realization has a single volume modulus @xmath123 . after all the complex structure moduli have been fixed by the fluxes , the superpotential for the volume modulus ( whose khler potential is @xmath124 ) is @xmath125 with @xmath126 , @xmath127 and @xmath128 parameters depending on the details of the setup .    the corresponding potential , before uplifting , has one ads supersymmetric critical point at finite volume and an asymptotic minkowski vacuum at infinity of the volume . it is easy to check that we can consistently set the axion @xmath129 to zero , since this is a fixed point for any value of the volume @xmath130 . for the sake of simplicity we will also take the parameters @xmath127 , @xmath128 and @xmath126 to be real ( and positive the first two @xmath131 , @xmath132 , negative the last one @xmath133 ) . the potential for this case and the covariantly holomorphic superpotential @xmath30 are given by fig . [ kkltpot ] as functions of the volume modulus . [ kkltpot ]   ( blue upper curve ) and @xmath30 ( red curve below ) as functions of the modulus . one can see that the potential has an ads minimum at @xmath134 and tends to minkowski limit at large @xmath130 . @xmath30 is always negative in the interval between ads and minkowski vacua . potential energy density @xmath21 is shown in units @xmath135 , whereas the covariantly holomorphic superpotential @xmath30 is shown in units @xmath136.,title="fig : " ]    one finds that for these parameters during the total flow the covariantly holomorphic superpotential @xmath30 is negative and approaches zero at infinity . we can reduce the problem to a one - dimensional gradient - flow , where the bps equations follow by using @xmath137 we solve numerically the flow equations in the opposite direction ( for technical reasons ) , starting with the ads minimum towards the infinitely far minkowski . for this purpose we have to use ( [ summary ] ) : @xmath138 \phi^\prime ( r ) & = &   \partial_\phi |{\cal z}(r)|.\end{aligned}\ ] ] the ads critical point is at @xmath134 , and for values slightly above @xmath139 we get a domain wall solution that drives @xmath140 and @xmath141 . in this case we obviously have just one basin of attraction and the attractor point is the supersymmetric ads vacuum . asymptotically , at large @xmath41 , the value of @xmath130 grows , and @xmath115 becomes @xmath142 we can therefore solve the bps equations above in this limit by @xmath143 it should be noted that @xmath144 which means that the curvature of the space - time ( [ scalari ] ) vanishes asymptotically , as it should be in minkowski space .    , the covariantly holomorphic superpotential @xmath30 ( in units @xmath145 ) , the warp factor @xmath128 and the curvature @xmath146 ( in units @xmath147 ) , all as functions of @xmath41 . here @xmath41 is the coordinate of the domain wall configuration given in eq . ( [ typei ] ) ; we show @xmath41 in units @xmath148 . , title="fig:"]0.5 cm , the covariantly holomorphic superpotential @xmath30 ( in units @xmath145 ) , the warp factor @xmath128 and the curvature @xmath146 ( in units @xmath147 ) , all as functions of @xmath41 . here @xmath41 is the coordinate of the domain wall configuration given in eq . ( [ typei ] ) ; we show @xmath41 in units @xmath148 . , title="fig : " ] 0.5 cm , the covariantly holomorphic superpotential @xmath30 ( in units @xmath145 ) , the warp factor @xmath128 and the curvature @xmath146 ( in units @xmath147 ) , all as functions of @xmath41 . here @xmath41 is the coordinate of the domain wall configuration given in eq . ( [ typei ] ) ; we show @xmath41 in units @xmath148 . , title="fig:"]0.5 cm , the covariantly holomorphic superpotential @xmath30 ( in units @xmath145 ) , the warp factor @xmath128 and the curvature @xmath146 ( in units @xmath147 ) , all as functions of @xmath41 . here @xmath41 is the coordinate of the domain wall configuration given in eq . ( [ typei ] ) ; we show @xmath41 in units @xmath148 . , title="fig : " ]    in fig . [ kkltabcd ] we plot the wall : first on the top right we have a scalar interpolating between the ads value @xmath134 and infinity . on the top right we plot @xmath30 as the function of wall coordinate @xmath41 . it is always negative but here we see clearly that it takes a constant negative value at ads and starts interpolating towards zero at minkowski vacuum . the next plot on the bottom left shows the warp factor @xmath149 . finally on the right bottom we show the curvature of the domain wall : it is vanishing at large values of modulus in the flat space and grows towards the decreasing values of modulus and sharply drops down to its negative ads value upon reaching the ads vacuum .    from this example we learn that in the absence of a barrier between the ads and minkowski vacua there exists a static domain wall where the flow of the volume modulus from the infinite value at minkowski towards the finite value at ads is controlled by 1/2 of unbroken supersymmetry . the tension of this wall is equal to @xmath150 where @xmath151 is the cosmological constant of the ads minimum . the previous example had only one ads vacuum . however , in the landscape of flux compactifications there are also examples with multiple ads vacua and therefore it is interesting to display an example with multiple basins of attraction to different ads vacua and the supersymmetric domain wall solutions interpolating between them . in order to find a workable example , we need an analytic superpotential admitting multiple ads vacua . such an example is provided by various models of toroidal orbifolds analyzed in @xcite . the generic moduli content is given by the axio / dilaton @xmath152 and by the complex - structure moduli @xmath153 and 3 khler moduli @xmath154 , with khler potential @xmath155 for the sake of simplicity we focus on the case of one untwisted complex structure modulus , as for instance in the @xmath156 , @xmath157 orbifolds or in the @xmath158 orbifold with @xmath159 , @xmath160 , @xmath161 or @xmath162 lattice . the superpotential is a polynomial in the complex - structure and dilaton moduli that comes from the fluxes , added by some non - perturbative exponential terms involving the khler - structure moduli . if we consider the cases with just one complex structure modulus , the full untwisted superpotential reads @xcite @xmath163 here the parameters @xmath164 and @xmath165 are related to the fluxes ( they are the coefficients of the real sections of the 3-cycles ) , while the parameters @xmath166 and @xmath167 are to be associated with non - perturbative effects . this potential includes the one - modulus kklt scenario as a special example . we can then look for supersymmetric vacua , satisfying @xmath168 , @xmath169 and @xmath170 . for simplicity , let us consider the case that @xmath171 , @xmath172 and look for vacua where @xmath173 . then , the supersymmetry conditions become @xmath174 b+ d\ , s & = & \frac{1}{u+\bar u } w , \label{duw } \\[2 mm ] - gh\ , { \rm e}^{-h t } & = & \frac{1}{t+\bar t } w. \label{dtw } \end{aligned}\ ] ] a generic solution for @xmath175 of these equations requires that @xmath176 and if we further look for the case where im@xmath177 = im@xmath152 = 0 , then we also obtain that @xmath178 , \label{cond1 } \\[2 mm ] { \rm re}\,t & = & \frac{3c{\rm re}\,s(b + d\ , { \rm re}\,s)}{h(a b- cd \,{\rm re}\,s^2)}. \label{label3 } \end{aligned}\ ] ] when the values of the fluxes are chosen appropriately , one can obtain 2 distinct vacua at values of @xmath179 and @xmath180 . one instance is the choice @xmath181 . the critical points are at @xmath182 and @xmath183 . ( the value of the size of the manifold is tuned by @xmath184 . the smaller @xmath101 , the bigger @xmath93 . increasing @xmath185 on the other hand allows for smaller values of @xmath186 . ) we can then use the bps equations to obtain the interpolating solution between them . having seen how the critical points are obtained , we would like to reduce the problem to a two - dimensional system depending only on @xmath187 , @xmath188 . both critical points lie on the section of the moduli space where all the axions are set to zero . this condition can be kept consistently all along the flow . when the axions are vanishing , @xmath189 and therefore @xmath190 . this condition can be preserved along the flow as @xmath191 if the axions are vanishing . this also implies , since the metric is diagonal , that the real and imaginary parts decouple consistently . given the symmetry of the potential , we can also assume that @xmath173 all along the flow . the bps equations therefore reduce to @xmath192 t^{\prime } & = & \mp 8t^2   \partial_t |{\cal z}|(s , t , u ) , \label{boh } \\[2 mm ] u^{\prime } & = & \mp 8u^2   \partial_u |{\cal z}|(s , t , u ) , \nonumber \end{aligned}\ ] ] and , by using ( [ relus ] ) , one can prove that the differential equations for @xmath193 and @xmath186 are proportional to each other . this means that we can further reduce the problem to the description of the gradient flows for the flow potential @xmath194 , where @xmath195 and @xmath196 . in these variables , the bps equations become @xmath197      \displaystyle \frac{a^{\prime}}{a } & = & \pm |{\cal z}|(\phi_i )   \end{array } \right.\label{bpseff}\ ] ] where the absolute value of the central charge reads @xmath198 notice the double exponentials due to the non - perturbative terms in the khler moduli .    given this flow potential , we can look for the special case where we have two distinct critical points , for instance when the parameters are chosen as above : @xmath181 . for these parameter values ( [ z ] ) shows two supersymmetric ads critical points with different values of the cosmological constant and an asymptotically minkowski vacuum . although the bps equations can be written in a closed form , the integration can only be done numerically . in fig . [ fig : contplotw ] we show the contour plot of @xmath30 , where the @xmath199 modulus ( proportional to the complex structure and to the dilaton @xmath195 ) is on the horizontal axis and @xmath200 , proportional to the khler modulus @xmath196 is on the vertical one .     and the gradient flow between the two critical points . @xmath201 is everywhere positive between ads critical points , it vanishes only at the asymptotically minkowski vacuum . ] we plot the total potential on fig . [ fig : stupot ] .    . ]    we can clearly see from the pictures that there are two distinct critical points , and we have plotted the numerical solution interpolating between them . as it is clear from the picture , the two critical points have different behavior . using the nomenclature coming from 5d , one is approached by the flow in an ultra - violet direction ( i.e.  the warp factor goes as @xmath202 , for @xmath203 ) , the other in an infra - red ( i.e.  the warp factor goes as @xmath204 for @xmath205 ) . an ultra - violet critical point is a minimum of the absolute value of the superpotential ( minimum with positive @xmath28 or maximum with negative @xmath28 ) , an infra - red critical point is a maximum . this implies that there is no need for the superpotential to cross zero . the warp factor indeed shows no critical points , as it always increases along the flow . the superpotential will necessarily cross zero along the flow and therefore the bps equations will change sign only when the interpolation must be done between two infra - red ( two ultra - violet critical points are excluded as the wall should be singular to interpolate between them ) . this is the only way to obtain a smooth domain wall interpolating between them . a similar phenomenon occurred in the realization of the randall  sundrum scenario with a thick brane in @xcite . we will have examples of such domain walls in sec . [ sec : subcritical ] .    from fig . [ fig : contplotw ] we can also clearly see the two basins of attraction for the ads and dine - seiberg minkowski critical points and the ads repeller ( with a marginal line of attraction separating the two basins ) . this shows that there is one main basin of attraction for the ads vacuum with the larger absolute value of the cosmological constant . there is a line of attraction for the other critical point , which is really a saddle point of the flow potential . finally there is another basin of attraction for the values of the moduli that leads to the decompactification of the space . the size modulus @xmath93 runs to infinity at the same time as the dilaton goes to zero , so we get a decompactification of the space in the weak coupling regime . this example shows that the landscape attractors are different from the black - hole attractors . for the latter , supersymmetry vacua are always minima of the superpotential @xmath28 and therefore one always has clear distinct basins of attractions and there are no repeller - like vacua . moreover , if one is looking at the 5-dimensional case there is always only one basin of attraction for the vector - multiplet moduli .      the panorama of the landscape can be further refined by considering a representative example in the class of type ii compactifications on calabi  yau threefolds , defined as hypersurfaces in a weighted projective space . when fluxes are active for the rr and ns three - form field strengths on the orientifold of these hypersurfaces , a non - trivial superpotential is induced for the dilaton and complex structure moduli , leading to the many interesting physical properties that have driven a massive research activity ( for a review and bibliography see @xcite ) . these models are particularly relevant even for a single deformation parameter because of the results obtained in @xcite showing that , quite unexpectedly , the density of flux vacua is not uniform in moduli space , but rather it peaks around the conifold singularity , @xmath206 . indeed , upon extensive numerical computations , it has been shown that the conifold region plays the role of  attractor " basin and it serves as accumulation region for flux vacua . this supports the analytic predictions made in @xcite . on the other hand , also the landau  ginzburg singularity at @xmath207 enjoys very peculiar properties , deriving from it being a fixed point of the modular group . in general , it is the special point in moduli space where the number of vacua with vanishing superpotential @xmath208 and discrete symmetries can be of the same order as the total number of vacua ( when non - zero ) @xcite . for our purposes , it will be convenient for computations because of our need to use perturbative expansions for the cy periods , that are known around this point @xcite . thus , in spite of neglecting the khler structure moduli that in principle would be present in these systems , this case is particularly interesting as a representative of the large amount of possible flux vacua that are making the landscape so rich . another important related feature is that one may encounter multiple attraction basins @xcite . these consist of distinct supersymmetric vacua that can be found within a finite region in moduli space , or `` area code '' , for fixed values of the @xmath209 and @xmath210 fluxes . we thus proceed to identify such a situation and to construct the domain wall solution interpolating between different area codes . as a start , we quickly recall some necessary formulas , while referring for instance to @xcite for the relevant lore and more conventions . one examines in general the @xmath211 complex structure deformations . as usual , one chooses a symplectic basis @xmath212 for the @xmath213 homology cycles ( @xmath214 ) , with dual cohomology elements @xmath215 such that @xmath216 and one assembles into a symplectic vector the periods of the unique holomorphic three - form @xmath217 , @xmath218 the @xmath219 are projective coordinates on the complex structure moduli space and @xmath220 , where @xmath221 is the holomorphic prepotential . the khler potential for the complex - structure and axion - dilaton moduli space is given by @xmath222 where @xmath223 is the axion - dilaton and @xmath224 . the superpotential induced by non - trivial 3-form fluxes reads @xmath225 where @xmath226 and @xmath227 are the integer valued flux components @xmath228 we will take for convenience @xmath229 . let us consider in the following @xmath230 to be the calabi  yau locus defined by the equation @xmath231the supersymmetric vacua are found by solving the  f  flatness " condition , that is a ( khler covariantized ) extremization with respect to @xmath223 and the single complex - structure modulus @xmath232 , @xmath233 leading to the equivalent conditions @xmath234 alternatively , one can solve the flux vacua attractor equations @xcite . one finds that the moduli space of this model presents landau  ginzburg , conifold and large complex structure singular points respectively at @xmath207 , @xmath206 and @xmath235 @xcite . the monodromy group , @xmath236 , of the complex structure moduli space has two generators : @xmath127 , which generates phase rotations @xmath237 with @xmath238 around the landau  ginzburg point at @xmath239 , and @xmath177 which corresponds to the logarithmic monodromy @xmath240 around the conifold singularity @xmath241 . by itself , @xmath127 generates a @xmath242 subgroup , with an associated fixed point at @xmath239 ; @xmath177 , on the other hand , is of infinite order . the monodromy matrix @xmath127 generates rotations by a root of unity around @xmath243 : @xmath244 and is explicitly given by @xmath245 the periods admit , for @xmath246 , an expansion in a picard  fuchs basis @xmath247 where @xmath129 is the @xmath248 root of unity @xmath249 the usual choice for the periods in this basis appropriate for the polynomial m is @xmath250 . then , the periods in the symplectic basis @xmath251 are related to the picard  fuchs periods by a linear transformation @xmath252 the matrix of transformation from the picard  fuchs to the symplectic basis has rational elements that are given by @xmath253 quite generally , in the symplectic basis the periods have the expansion @xmath254 where the @xmath255 are complex constants and @xmath256 are constant four - vectors with components that are rational linear combinations of powers of @xmath129 .    as a further computational tools , we note in passing that for this model the meijer periods ( we find them following method developed in @xcite for quintic in the context of black hole attractors ) are given by the following expression in mathematica notations @xmath257,\\ & u^{-}_{1}(z ) & = { \rm meijerg}\left[\left\{\left\{\frac16,\frac26,\frac46,\frac56\right\},\{\}\right\},\{\{0,0\},\{0,0\}\},z\right],\\ & u^{-}_{2}(z ) & = { \rm meijerg}\left[\left\{\left\{\frac16,\frac26,\frac46,\frac56\right\},\{\}\right\},\{\{0,0,0\},\{0\}\},-z\right],\\ & u^{-}_{3}(z ) & = { \rm meijerg}\left[\left\{\left\{\frac16,\frac26,\frac46,\frac56\right\},\{\}\right\},\{\{0,0,0,0\},\{\}\},z\right ] , \end{aligned}\ ] ] where @xmath258 are valid for im@xmath259 .    the transformation to the picard  fuchs basis is given by the matrix @xmath260 so that @xmath261 with @xmath262 for im@xmath263 . we now proceed to specifying the vacua for given values of the fluxes . let us take , as a convenient choice among the solutions of the f - flatness conditions , the flux components to be @xmath264 and @xmath265 the virtue of this set of numbers is that the consequent vacua fall well within the convergence region of the picard  fuchs periods , thus making the numerical computations for the flow somewhat easier . it turns out that these fluxes give rise to two vacua : the first vacuum is a maximum at the landau  ginzburg point , with @xmath266 the critical value of the potential is @xmath267 . we find the following diagonalized mass matrix for this vacuum @xmath268 the second vacuum is a saddle point close to the conifold point , with @xmath269 and with potential equal to @xmath270 . in this case the diagonalized mass matrix is @xmath271    the interpolating solution going from one vacuum to the other is found by solving the flow equations ( [ renorm1 ] ) , that become @xmath272 where we have used for convenience the variable @xmath273 . the resulting flows are displayed below on the @xmath232 complex plane ( figure [ plot1a ] ) and on the @xmath223 complex plane ( figure [ plot2 ] ) . complex plane . ]     complex plane . ]    the corresponding flow of the potential @xmath274 as a function of @xmath275 is shown in fig . [ plot3 ] . it shows clearly that there is no barrier .    . ]    these domain walls are associated to flows that go from a maximum in the landau  ginzburg point , towards the conifold which is a saddle point . these are explicit examples of critical points which may be sometimes repellers ( for saddle points or maxima of the potential ) and sometimes attractors , in particular stable attractors in case of minima of the potential . in this case , the flow tends to escape from the lg point towards the conifold . we are confident that these results on interpolating domain walls are common to other one - parametric calabi  yau models that have been studied in @xcite and generically have similar distribution of vacua on moduli space . when the vacua are separated by a barrier , one can expect to find an absolute minimum as well as a meta - stable local minimum . this class of models in the stringy landscape was suggested in @xcite : in addition to the ads minimum and to the asymptotically minkowski vacuum , there is also a supersymmetric minkowski minimum at some finite value of the volume modulus . the main motivation for the new class of models with such a minimum was to have more flexibility with the gravitino mass , and to avoid certain cosmological problems specific to the standard kklt scenario . here we will show that the same class of models ( which we will call kl models ) allows two ads minima , both at finite value of the volume modulus . moreover , we will find two types of flows here : depending on the choice of parameters the superpotential may never vanish between the two ads vacua , or it may vanish at some point and therefore change the sign along the flow . this results in several different types of bps domain walls in the kl models . the model has one modulus as in the kklt case , but a more involved superpotential of the racetrack type @xcite . the superpotential is a modification of the kklt model by an additional exponential term @xmath276 this model is a simple instance of multiple basins of attraction . it has one supersymmetric minkowski vacuum at fixed volume , one supersymmetric ads vacuum also at fixed volume and one supersymmetric minkowski asymptotic vacuum at infinite volume as shown in an upper curve in figure [ klpot ] . the plot of the potential clearly shows the barrier between the left minkowski minimum and the ads minimum , the superpotential is negative everywhere away from minkowski vacua . [ h ! ] , blue upper curve , and the covariantly holomorphic superpotential @xmath30 , the red curve below , as the functions of the modulus @xmath130 for the kl model . the potential is shown in units @xmath277 , and @xmath30 is in units @xmath278 . both @xmath21 and @xmath30 have extrema at the same point , which is a reflection of the fact that supersymmetry is unbroken in the minima of @xmath21 . this plot is for @xmath131 , @xmath279 , @xmath280 , @xmath281 and @xmath282.,title="fig : " ]    for the choice @xmath131 , @xmath279 , @xmath280 , @xmath281 and @xmath282 , also using once more the definition @xmath283 , the vacua are at @xmath284 , @xmath285 and @xmath286 , with @xmath287 . therefore , one can construct two interpolating domain walls . one is between the ads minimum and the local minkowski minimum , which we will call m@xmath288 . another wall is between the ads minimum and the minkowski extremum at infinity , m@xmath289 . here we study only the domain walls interpolating between m@xmath288 and ads , see figure [ plotx ] . the features of the domain wall between m@xmath289 and ads minimum in this model will be the same as the ones studied in kklt case in section [ sec : kkltmodel ] .     and ads . we plot the scalar @xmath130 , @xmath30 ( in units @xmath290 ) , warp factor @xmath128 and the curvature scalar @xmath146 ( in units @xmath277 ) as functions of the coordinate @xmath41 . the curvature @xmath146 is equal to zero in the minkowski vacuum , has a peak at the domain wall , and approaches a constant negative value in the ads vacuum.,title="fig:"]0.5 cm   and ads . we plot the scalar @xmath130 , @xmath30 ( in units @xmath290 ) , warp factor @xmath128 and the curvature scalar @xmath146 ( in units @xmath277 ) as functions of the coordinate @xmath41 . the curvature @xmath146 is equal to zero in the minkowski vacuum , has a peak at the domain wall , and approaches a constant negative value in the ads vacuum.,title="fig : " ] 0.7 cm   and ads . we plot the scalar @xmath130 , @xmath30 ( in units @xmath290 ) , warp factor @xmath128 and the curvature scalar @xmath146 ( in units @xmath277 ) as functions of the coordinate @xmath41 . the curvature @xmath146 is equal to zero in the minkowski vacuum , has a peak at the domain wall , and approaches a constant negative value in the ads vacuum.,title="fig:"]0.5 cm   and ads . we plot the scalar @xmath130 , @xmath30 ( in units @xmath290 ) , warp factor @xmath128 and the curvature scalar @xmath146 ( in units @xmath277 ) as functions of the coordinate @xmath41 . the curvature @xmath146 is equal to zero in the minkowski vacuum , has a peak at the domain wall , and approaches a constant negative value in the ads vacuum.,title="fig : " ]    in figure [ plotx ] we plot on the top left the scalar potential @xmath21 , and on the right the covariantly holomorphic superpotential @xmath30 , which vanishes at m@xmath288 and stays negative all the way towards ads . the warp factor on the bottom left starts as a constant at m@xmath288 and interpolates into an exponent of @xmath41 near the ads . finally , the curvature jumps from zero at m@xmath288 to some positive value and falls down to a negative constant at ads . it is interesting that the curvature of the wall grows significantly between zero at m@xmath288 before becoming negative at ads . now we are going to consider the models with two ads minima . this can be done , e.g. , by taking the kl model with a supersymmetric minkowski minimum m@xmath291 , and change any of the parameters of this model . in particular , as one can see from fig.[fig : klppota ] , one can obtain an ads minimum from the minkowski minimum by a slight modification of a single parameter in the previous model : we changed @xmath292 to @xmath293 . using once more the definition @xmath283 , the ads vacua are at @xmath284 , @xmath294 and @xmath295 , with @xmath296 , and @xmath297 for @xmath284 . the potential @xmath21 ( upper curve ) and the covariantly holomorphic superpotential @xmath30 ( lower curve ) for the kl model as function of the modulus are shown in figure [ fig : klppota ] , for the choice @xmath298 , @xmath279 , @xmath280 , @xmath281 and @xmath282 . @xmath30 is always negative . ( upper curve , in units @xmath277 ) and the covariantly holomorphic superpotential @xmath30 ( lower curve , in units @xmath278 ) for the kl model , as function of the modulus . notice that we took @xmath298 instead of @xmath292 , all other parameters being the same as for the potential in fig . [ klpot].,title="fig : " ]    we show various features of the kl wall between the two ads vacua in fig . [ fig : plota ] . we plot the scalar and the covariantly holomorphic superpotential , which is a negative constant near the first ads vacuum and interpolates towards another negative value when it approaches another ads vacuum . the warp factor @xmath128 starts as a decreasing exponent and then changes into a faster decreasing exponent of @xmath41 near the second ads . the curvature @xmath146 jumps from some constant negative value near the first ads vacuum , becomes positive near the domain wall , and then falls down to a negative constant near the second ads vacuum .     and we plot the scalar @xmath130 and the covariantly holomorphic superpotential @xmath30 ( in units @xmath290 ) , which is always negative , warp factor @xmath128 , and the curvature scalar @xmath146 ( in units @xmath277 ) , as a function of the coordinate @xmath41 . the curvature @xmath146 has a peak at the domain wall , and approaches different constant negative values in two different ads minima.,title="fig:"]0.5 cm   and @xmath299 . we plot the scalar @xmath130 and the covariantly holomorphic superpotential @xmath30 ( in units @xmath290 ) , which is always negative , warp factor @xmath128 , and the curvature scalar @xmath146 ( in units @xmath277 ) , as a function of the coordinate @xmath41 . the curvature @xmath146 has a peak at the domain wall , and approaches different constant negative values in two different ads minima.,title="fig : " ] 0.7 cm   and @xmath299 . we plot the scalar @xmath130 and the covariantly holomorphic superpotential @xmath30 ( in units @xmath290 ) , which is always negative , warp factor @xmath128 , and the curvature scalar @xmath146 ( in units @xmath277 ) , as a function of the coordinate @xmath41 . the curvature @xmath146 has a peak at the domain wall , and approaches different constant negative values in two different ads minima.,title="fig:"]0.5 cm   and @xmath299 . we plot the scalar @xmath130 and the covariantly holomorphic superpotential @xmath30 ( in units @xmath290 ) , which is always negative , warp factor @xmath128 , and the curvature scalar @xmath146 ( in units @xmath277 ) , as a function of the coordinate @xmath41 . the curvature @xmath146 has a peak at the domain wall , and approaches different constant negative values in two different ads minima.,title="fig : " ]    the tension of this bps wall is given by eq . ( [ tension1 ] ) . @xmath300 it is critical in a sense explained in @xcite : the tension of such bps brane corresponds to a limit of the cdl bubble solution @xcite with infinite radius . when such bps wall exist there is no tunneling since these supersymmetric walls saturate the cdl bound . we will discuss this issue later when we will study tunneling , see section [ decay ] .      in the previous subsection we considered a kl model with two ads minima and the superpotential which did not change sign on the way from one minimum to another . now we will construct a kl model with two ads minima and superpotential vanishing at some point between them ( fig . this can be done by a small deviation of one of the parameters : instead of taking @xmath292 as in original kl model , or @xmath301 , as in the previous subsection with two ads vacua , it is sufficient to take @xmath302 . for example , for the choice @xmath303 , @xmath279 , @xmath280 , @xmath281 and @xmath282 , also using once more the definition @xmath283 , one finds that the ads vacua are at @xmath284 , @xmath304 , with @xmath305 , and at @xmath284 , @xmath306 , with @xmath307 . the covariantly holomorphic superpotential @xmath30 is positive near the ads minimum at @xmath308 . it crosses zero at @xmath309 and becomes negative when approaching the second ads minimum at @xmath306 . it approaches zero at infinite values of the volume modulus . ( upper curve , in units @xmath310 ) and the covariantly holomorphic superpotential @xmath30 ( lower curve , in units @xmath311 ) for the kl model as function of the modulus . @xmath30 is positive at @xmath308 and negative at @xmath312 . here we took @xmath303 instead of @xmath292 , all other parameters being the same as for the potential in fig . [ klpot].,title="fig : " ]    the corresponding bps solution is shown in fig . [ plotb ] . we plot on the top left the scalar field @xmath130 , and on the right the covariantly holomorphic superpotential , which changes the sign on the solution . the warp factor on the bottom left starts as a growing exponent near the first ads minimum and changes into a decreasing exponent of @xmath41 near the second minimum . finally , the curvature jumps from some constant negative value near the first ads vacuum , becomes positive near the domain wall , and then falls down to a negative constant near the second ads vacuum . as far as we know , this solution provides the first explicit example of a smooth bps domain walls in @xmath313 with the superpotential vanishing between the critical points .     changes the sign on the domain wall solution . we plot the scalar @xmath130 , the covariantly holomorphic superpotential @xmath30 ( in units @xmath290 ) , the warp factor @xmath128 , and the curvature @xmath146 ( in units @xmath277 ) . notice the characteristic peak of the scale factor at the point where @xmath30 changes sign.,title="fig:"]0.5 cm   changes the sign on the domain wall solution . we plot the scalar @xmath130 , the covariantly holomorphic superpotential @xmath30 ( in units @xmath290 ) , the warp factor @xmath128 , and the curvature @xmath146 ( in units @xmath277 ) . notice the characteristic peak of the scale factor at the point where @xmath30 changes sign.,title="fig : " ] 0.3 cm   changes the sign on the domain wall solution . we plot the scalar @xmath130 , the covariantly holomorphic superpotential @xmath30 ( in units @xmath290 ) , the warp factor @xmath128 , and the curvature @xmath146 ( in units @xmath277 ) . notice the characteristic peak of the scale factor at the point where @xmath30 changes sign.,title="fig : " ]   changes the sign on the domain wall solution . we plot the scalar @xmath130 , the covariantly holomorphic superpotential @xmath30 ( in units @xmath290 ) , the warp factor @xmath128 , and the curvature @xmath146 ( in units @xmath277 ) . notice the characteristic peak of the scale factor at the point where @xmath30 changes sign.,title="fig : " ]    the tension of this bps wall is given by eq . ( [ tension2 ] ) . @xmath314 it is super - critical in the sense explained in @xcite . the tension of such bps walls exceeds the cdl bound and therefore the tunneling is super - suppressed . in previous sections we studied domains with unbroken supersymmetry separated by static bps domain walls . in our world supersymmetry is broken , but if supersymmetry breaking is not very large as compared to other parameters , one may hope that our domain wall solutions will provide an approximate description of slowly moving near - bps domain walls . just as in the situation with the near - extremal black holes , one may hope the physics of the near - bps walls can be under a much better theoretical control than the physics of generic domain walls . it would be especially interesting to see whether the bps domain walls may play any role in the cosmological context , dividing the universe into different regions separated by static domain walls . however , from the point of view of the friedmann cosmology , the universe with a negative cosmological constant @xmath315 should collapse to a cosmological singularity within a very short time @xmath316 . therefore in the usual cosmological context one can hardly describe coexisting ads domains separated by static domain walls . one of the most interesting applications of the results obtained in the last section is the theory of vacuum decay and bubble formation . it is well known that ds vacua obtained in the kklt construction by uplifting of ads vacua are unstable with respect to the decay to the minkowski vacuum corresponding to infinitely large values of the volume modulus @xmath130 , i.e. to a supersymmetric 10d vacuum @xcite .    according to @xcite , the decay probability is given by @xmath317 where @xmath318 is the euclidean action for the tunneling trajectory , and @xmath319 is the euclidean action for the initial metastable configuration @xmath320 . the euclidean action calculated for the metastable vacuum ds solution @xmath321 is given by @xmath322 this action has a simple sign - reversal relation to the entropy of de sitter space @xmath323 : @xmath324 therefore the decay time of the metastable ds vacuum @xmath325 can be represented in the following way : @xmath326 here @xmath327 is the so - called recurrence time . it was shown in @xcite that if the decay of the metastable ds vacua occurs due to the tunneling through the barrier with a positive scalar potential , then @xmath318 is always negative , and therefore the decay always happens during the time shorter than the recurrence time @xmath328 . this result directly applied to the simplest kklt model where the tunneling occurs through the positive barrier separating the metastable ds vacuum and the supersymmetric 10d vacuum . however , the situation with the tunneling to the minima with a negative vacuum density remained less clear ; see refs . @xcite for the recent discussion . in what follows we will examine this issue using the results obtained in the previous section .    our consideration will be based on the observation of ref . @xcite that in all cases when @xmath329 does not vanish across the domain wall , the domain wall solutions can be represented as a limiting configuration describing the wall of the cdl bubbles of an infinitely large radius . for such bubbles , the tunneling action is also infinitely large , and the vacuum decay is impossible . in our context , this fact is related to the supersymmetry of the different vacua @xcite , and of the interpolating bps wall solution . however , in realistic situations supersymmetry is broken . for example , in the kklt construction the metastable ds vacuum state appears after the uplifting of a supersymmetric ads vacuum . uplifting can be achieved , e.g. , by adding a @xmath330 brane contribution , which results in a supersymmetry breaking in the ds vacuum . similarly , in the kl model @xcite one must uplift a minkowski vacuum or an ads vacuum to a ds vacuum , thus breaking the supersymmetry . as we will see , this may lead to a relatively rapid decay of the uplifted ds vacuum due to creation of the bubbles describing a collapsing universe with negative vacuum energy corresponding to the ads vacuum . note that in terms of a canonically normalized volume modulus @xmath130 , which we will use hereafter , the @xmath330 brane contribution looks as @xmath331 , where @xmath332 is some constant . ( the canonically normalized field , @xmath333 , is obtained by dividing the field @xmath130 used in the previous section by the factor @xmath334 . ) this contribution rapidly decreases at large @xmath130 . therefore this term may significantly uplift the ads ( or minkowski ) minimum at small values of @xmath130 , but it makes a much smaller effect on the ads minima at large @xmath130 . this suggests that the transitions should typically occur towards larger values of the volume of the compactified space . first of all , consider the kl potential with the supersymmetric minkowski vacuum , and uplift it by a tiny amount to make the vacuum energy there equal to the present cosmological constant , @xmath335 , in planck energy density units . this should make this ds vacuum metastable , and lead to a vacuum decay due to creation of bubbles with negative vacuum energy .    . ]    to simplify the investigation of its decay , we will consider potentials such that the difference between the depth of the two minima is much smaller than the height of the barrier . an example of such a potential is provided by the kl model with parameters @xmath336 , @xmath337 , @xmath338 , @xmath339 and @xmath340 . for these values of parameters the potential has a supersymmetric minkowski minimum at @xmath341 and an ads minimum at @xmath342 . if one slightly changes the parameters , e.g. takes @xmath343 , then the minkowski minimum becomes an ads minimum . if one uplifts the potential by adding a term @xmath331 , the potential acquires a ds minimum , as shown in fig . [ klthin ] . during all of these changes , the second ads minimum changes only slightly . one expects that the thin wall approximation should be valid for the description of tunneling in such potentials . one part of this condition is automatically satisfied for the near - bps bubbles . indeed , by considering sufficiently small uplifting one can always make such bubbles arbitrarily large without changing the thickness of the wall . however , in order to use the thin - wall approximation as formulated in ref . @xcite , one should also check that the difference between the depth of the two minima is much smaller than the height of the barrier , and that the width of the wall of the cdl bubble remains much smaller than the radius of the bubble even if one turns off the gravitational interaction . we checked numerically that this condition is also satisfied for our potential . having in mind potentials of this type , we will use the thin - wall approximation to investigate the decay of the ds vacuum with the field @xmath308 and vacuum energy density @xmath344 to the state with the field @xmath312 with a negative vacuum energy density @xmath345 . we will be interested in the realistic situation @xmath346 , as illustrated by fig . [ klthin ] . examples of the vacuum decay considered in @xcite described transitions between ds and minkowski space , and decay of the minkowski vacuum . we will use the results obtained by parke @xcite , who generalized the results of @xcite for the tunneling between the states with arbitrary values of @xmath347 in the thin - wall approximation .    according to @xcite , the bubble size is given by @xmath348 where @xmath349 is the bubble radius in the thin wall approximation ignoring gravity , @xmath111 is the bubble wall tension , and @xmath350 ( it can be positive or negative ; in our case it is negative ) , and @xmath351 in the thin - wall approximation , the tension of the wall can be approximated by the integral @xmath352^{1/2 } \ .\ ] ] here @xmath353 is the point near @xmath312 where @xmath354 . @xcite uses a slightly different prescription , which is equivalent to ours in the limit when the height of the barrier is much greater than the difference between @xmath347 . ) in the limit when the radius of the cdl bubble becomes infinitely large , the wall tension is known exactly : @xmath355 see @xcite , and also our derivation of this result in sect . 2 , eq .  ( [ tension1 ] ) .    using the variables @xmath356 one can represent eq . ( [ rad ] ) in the following convenient form : @xmath357 for a particular case of tunneling from minkowski space ( @xmath358 ) , the infinite bubble size regime ( the threshold regime for the tunneling ) corresponds to @xmath359 , @xmath360 , when the denominator vanishes .    the tunneling probability is given by @xmath361 , where @xmath362 here the first term is the no - gravity result , and the gravitational correction factor @xmath363 is given by @xcite @xmath364}{x^2(y^2 - 1)(1 + 2xy+x^2)^{\frac{1}{2 } } } \ . \label{gravfactor}\ ] ] this factor also blows up at @xmath359 , @xmath360 , i.e. the tunneling becomes forbidden . now let us uplift the minkowski vacuum , assuming that the ads vacuum at @xmath312 remains approximately at the same height @xmath365 . as we already argued , this is indeed the case if the ads vacuum is at large values of the volume modulus @xmath130 . uplifting changes not only @xmath366 , but also the bubble wall tension @xmath111 . typically uplifting decreases the wall tension , which increases the probability of the tunneling ; a sufficiently large uplifting may remove the wall altogether , see @xcite . any small uplifting makes the tunneling from @xmath308 to @xmath312 possible . in order to study the size of the bubbles and the probability of the tunneling , let us change @xmath367 and @xmath368 by @xmath369 and @xmath370 . in the linear approximation @xmath371 and both depend on the small uplifting @xmath372 . in general , one would expect that the result should depend both on @xmath369 and @xmath370 . however , we find that in the leading approximation in @xmath369 and @xmath370 , the results of the expansion near @xmath359 , @xmath360 depend only on @xmath370 , i.e. they are not sensitive to the small changes of the bubble wall tension . in particular , the denominator in ( [ size ] ) is equal to @xmath373 , where @xmath374 , for @xmath346 . this yields @xmath375 here @xmath376 is the hubble constant in ds space with vacuum energy density @xmath377 . thus , instead of an infinite size bubble , we have a bubble with the size equal to the size of the ds horizon !    now we will perform a similar investigation for the tunneling rate . by considering the function @xmath378 in eq . ( [ gravfactor ] ) and expanding it with respect to @xmath369 and @xmath370 , one finds that in the limit @xmath346 @xmath379 where @xmath380 is the entropy of ds space , @xmath381 this means that the tunneling probability is given by the universal model - independent equation which has a simple geometric interpretation in terms of ds entropy : @xmath382 several comments are in order . first of all , we have found that the decay rate of the vacuum obtained by uplifting of a supersymmetric minkowski minimum does not depend on the tension of the wall , or on the depth of the supersymmetric ads minimum . the only thing we needed to know is that prior to the uplifting these two minima were supersymmetric , which is a generic property of all known flux vacua . an unexpected feature of this result is that the tunneling is suppressed not by a factor @xmath383 , but by a factor @xmath384 . this means that the decay time of an uplifted minkowski vacuum will be shorter than the recurrence time @xmath385 by a huge factor @xmath386 . we will return to the discussion of this fact in section [ sink ] . as we have seen in the previous subsection , a tiny supersymmetry breaking achieved by uplifting of the supersymmetric minkowski vacuum to the ds vacuum with @xmath387 renders this vacuum unstable . now we will extend our results for a more general case . suppose that initially we have two supersymmetric ads vacua , with depth @xmath388 , such that @xmath389 . one can achieve it e.g. by considering the kl potential discussed in the previous section , but taking @xmath390 instead of @xmath391 . tunneling between the two supersymmetric ads minima is forbidden . let us also assume that the deeper minimum is at larger values of the volume modulus , so the kklt uplifting of the minimum @xmath366 to ds space with @xmath387 leaves the depth of the second ads minimum almost unaffected . this is indeed the case in the model discussed above . we will also assume that prior to the uplifting @xmath392 was much greater than the present value of the vacuum energy @xmath393 , so for all practical purposes it is not important whether we uplifted it to a ds state @xmath387 or to the minkowski state with @xmath358 . this subtle difference was important in our previous investigation because the degree of supersymmetry breaking was proportional to @xmath387 . now we are considering a much greater supersymmetry breaking , which is achieved due to a much more significant uplifting , from the initial value of @xmath366 to ( approximately ) zero . in order to simplify notation and avoid mixing the values of @xmath366 before and after the uplifting , we will assume that after the uplifting @xmath358 , and in all subsequent equations we will understand by @xmath366 its value in the ads minimum prior to the uplifting .    in order to study the decay rate of the ( nearly ) minkowski vacuum formed by uplifting , one should use , as before , eqs . ( [ coeffb ] ) and ( [ gravfactor ] ) . in minkowski vacuum , one has @xmath360 , so the expression ( [ gravfactor ] ) looks singular , but in fact this singularity is fictitious , which can be easily seen by expanding it near @xmath360 .    the outcome of the calculations can be represented in the following form : @xmath394 here @xmath111 is the bubble wall tension for the tunneling in the uplifted potential . the first term in this expression is the probability of the decay of the uplifted minkowski vacuum ignoring gravitational corrections . the second term represents the gravitational suppression of the probability of the tunneling due to the factor @xmath363 . as we see , the probability of the tunneling depends on the wall tension @xmath111 . however , the only part of the expression ( [ tunnads ] ) where this dependence is crucially important is the gravitational suppression term @xmath395 . if after uplifting we would still have two supersymmetric vacua ( minkowski and ads ) , then we would have @xmath396 , the last term in ( [ tunnads ] ) would blow up , and the tunneling would be forbidden . however , prior to the uplifting we had a smaller value of tension , @xmath397 . as we argued before , uplifting tends to decrease the barrier between the two vacua , and therefore it tends to decrease the tension . ( our numerical investigation of the tunneling in the kl model with the potential shown in fig . [ klthin ] is consistent with this expectation . ) if this is the case , the gravitational suppression gradually disappears , and the probability of the tunneling becomes very large .    to illustrate various options , let us assume , in the first approximation , that @xmath111 does not change during the uplifting . in this case we find that @xmath398 where @xmath399 .    for @xmath13 ( i.e. for @xmath365 of the same order as @xmath366 , prior to the uplifting ) , this expression coincides , up to a numerical coefficient @xmath400 , with the entropy of ds space with the vacuum energy @xmath392 , see eq . ( [ action2a ] ) . the main difference is that for the ds space in which we live now we have @xmath387 , and the factor @xmath401 is incredibly small , @xmath402 . meanwhile the depth of the ads vacuum prior to the uplifting can be a hundred orders of magnitude greater than @xmath393 , and therefore the decay rate of this vacuum after the uplifting can be very large ( see below for a numerical estimate ) . note that if , as we expect , @xmath111 slightly decreases during the uplifting in our model , then @xmath403 will be even smaller , and the decay rate will be even higher .    it would be nice to confirm our conclusions by other methods , or to calculate the tunneling probability directly , going beyond the thin - wall approximation . however , our main goal here was to show that generically there is no reason to expect a strong gravitational suppression of the decay of ds vacua if they are obtained by the kklt uplifting of the supersymmetric ads minima . qualitatively , the main argument is as follows : if there are many ads minima in the stringy landscape , and if the tunneling between them is suppressed by supersymmetry , then this suppression disappears after the uplifting and supersymmetry breaking . let us check whether these intuitive expectations match our estimates , eq . ( [ tunndsads ] ) . according to ( [ potential2 ] ) , the scalar potential prior to the uplifting is given by @xmath404 . in the supersymmetric minimum at @xmath308 we have @xmath405 , and @xmath406 . uplifting practically does not change the position of the minimum @xmath308 , and therefore it practically does not change @xmath407 . meanwhile the gravitino mass squared after the uplifting is given by @xmath408 . this means that the depth of the ads minimum before the uplifting is equal to @xmath409 , where @xmath16 is the present value of the gravitino mass @xcite . therefore our result for the tunneling rate of the uplifted ads state ( [ tunndsads ] ) can be represented in the form directly related to the gravitino mass : @xmath410    the scale of supersymmetry breaking in our world is much greater than the energy scale corresponding to the present value of the cosmological constant . for example , if the gravitino mass is @xmath411 gev @xmath412 in planck units , then @xmath413 . therefore eq .  ( [ tunndsadsgrav ] ) suggests that a typical decay rate of the ads vacua after their uplifting to the ds state with @xmath387 is greater than @xmath414 . most importantly , this decay rate is greater than @xmath402 by an enormously large factor @xmath415 . this conclusion may have important cosmological implications which we are going to discuss in the next section . but before doing so , let us point out a possible caveat in our discussion . as we already mentioned , if the tunneling between the two vacua should go through the region where @xmath30 changes the sign , then @xmath416 . in this case , if we ignore the decrease of the tension during the uplifting , the tunneling between the two vacua does not occur . one can clearly see it by looking at the gravitational suppression term @xmath395 . this term for @xmath416 is negative , which corresponds to super - suppression of the tunneling @xcite .    on the other hand , we expect that in the realistic theories with hundreds of moduli fields there should be many ads minima , some of which will remain ads after the uplifting . therefore it is natural to expect that the decay channels at least to some of these vacua will be open , and will be described by eqs . ( [ tunndsads ] ) , ( [ tunndsadsgrav ] ) . this is the only thing that matters to us , since the total rate of the vacuum decay is determined by the path of the least resistance , i.e. by the channels which maximize the tunneling probability . one of the main reasons to study vacuum transitions which occur on an incredibly large timescale is to find a probability distribution to live in one of the many vacuum states . this issue has a long history involving quantum cosmology , eternal inflation , and statistical properties of the stringy landscape .    for a long time , vacuum transitions were studied in the context of inflationary cosmology . the main subject was the possibility of transitions between ds vacuum and minkowski vacuum , or between two ds vacua . one of the most surprising and controversial results was obtained by hawking and moss @xcite . they found that a phase transition from a ds minimum with energy density @xmath366 may occur via a quantum jump to the top of the barrier separating the two vacua , with the height @xmath417 , and by a subsequent classical rolling to the second ds minimum with energy density @xmath418 . the probability of this jump , according to @xcite , was given by @xmath419 initial interpretation of this result was rather obscure because the corresponding instanton seemed to describe a homogeneous tunneling , which was certainly impossible in an infinite ( or exponentially large ) inflationary universe . moreover , from the derivation of this result it was not clear why the tunneling should occur to the top of the potential instead of going directly to the second ds minimum . a proper interpretation of the hawking  moss tunneling was achieved only after the development of the stochastic approach to inflation @xcite . one may consider quantum fluctuations of a light scalar field @xmath130 with @xmath420 . during each time interval @xmath421 this scalar field experiences quantum jumps with the wavelength @xmath422 and with a typical amplitude @xmath423 . then the wavelength of these fluctuations grows exponentially . as a result , quantum fluctuations lead to a local change of the amplitude of the field @xmath130 , which looks homogeneous on the horizon scale @xmath424 . from the point of view of a local observer , this process looks like a brownian motion of the homogeneous scalar field . if the potential has a ds minimum at @xmath425 with @xmath426 , then eventually the probability distribution to find the field with the value @xmath130 at a given point becomes time - independent , @xmath427 that is why the probability to gradually climb to the local maximum of the potential at @xmath428 and then fall to another ds minimum is given by eq . ( [ hm ] ) @xcite . note that the distribution @xmath429 is a fraction of the _ comoving _ volume of the universe corresponding to each of the ds vacua . this probability distribution does not take into account different rate of growth of different part of the universe . its interpretation can be understood as follows : at some initial moment one divides the universe into many domains of the same size , assigns one point to each domain , and follows the subsequent distribution @xmath429 of the points where the scalar field takes the value @xmath130 . a necessary condition for the derivation of this result in @xcite was the requirement that @xmath420 . this requirement is violated for all known scalar fields at the present ( post - inflationary ) stage of the evolution of the universe . however , one may try to generalize it for the situations with @xmath430 . following @xcite ( see also @xcite ) , we will look for a probability distribution @xmath431 to find a given _ point _ in a state with the vacuum energy @xmath347 . the main idea is to consider cdl tunneling between two ds vacua , with vacuum energies @xmath366 and @xmath365 , such that @xmath432 , and to study the possibility to tunneling in both directions , from @xmath366 to @xmath365 , or vice versa . the action on the tunneling trajectory , @xmath318 , does not depend on the direction in which the tunneling occurs , but the tunneling probability does depend on it . it is given by @xmath433 on the way up , and by @xmath434 on the way down @xcite . let us assume that the universe is in a stationary state , such that the comoving volume of the parts of the universe going upwards is balanced by the comoving volume of the parts going down . this can be expressed by the detailed balance equation @xmath435 which yields ( compare with eq . ( [ hm ] ) ) @xmath436 independently of the tunneling action @xmath318 .    equations ( [ e38a ] ) and ( [ weinb ] ) imply that the fraction of the comoving volume of the universe in a state @xmath130 ( or @xmath312 ) different from the ground state @xmath308 ( which is the state with the lowest , but positive , vacuum energy density ) is proportional to @xmath437 , with the normalization coefficient @xmath438 . the probability distribution @xmath439 coincides with the square of the hartle  hawking wave function describing the ground state of the universe @xcite . it has a simple physical meaning : the universe wants to be in the ground state @xmath308 with the lowest possible value of @xmath440 , and the probability of the deviations from the ground state is exponentially suppressed . this probability distribution also has a nice thermodynamic interpretation in terms of ds entropy @xmath441 @xcite : @xmath442 here , as before , @xmath443 . this result and its thermodynamic interpretation played a substantial role in the discussion of the string theory landscape @xcite .    unfortunately , there are some problems associated with this result . it suggests that the universe is similar to a pond filled with still water , where all things accumulate at the bottom : the fraction of the comoving volume with the smallest possible vacuum energy must be overwhelmingly large , @xmath444 . for a discussion of paradoxes associated with this conclusion and their possible resolution in the kklt scenario see , e.g. @xcite . investigation of the stationary probability distribution alone does not give us a full picture . for example , the probability distribution ( [ e38a ] ) tells us about the fraction of the comoving volume of the universe in a given state , but it tells us nothing about the evolution towards this state . a partial answer to this question can be given by investigation of the stochastic diffusion equations describing the evolution of the scalar field in the inflationary universe . but now , instead of looking for the most probable outcome of the evolution , one should follow the evolution backwards and look for the initial condition @xmath308 for the trajectories which bring the field to its final destination @xmath130 . in the stationary regime considered above , the corresponding solution looks very similar to ( [ e38a ] ) @xcite : @xmath445 in this equation , however , @xmath308 is not the position of the ground state , but a position of an arbitrary initial point for the diffusion process which eventually brings us to the point @xmath130 . as we see , the probability is maximized by the largest possible value of @xmath446 . interestingly , expression @xmath447 describing the probability of initial conditions coincides with the expression for the square of the tunneling wave function describing creation of a closed ds universe `` from nothing '' @xcite , whereas the second term looks like the square of the hartle - hawking wave function describing the ground state of the universe . in the stationary regime the squares of these two wave functions coexist in the same equation , but they provide answers to different questions .    here we should note that it is far from being clear whether the simplest probability distribution which we studied so far should be used in anthropic considerations . indeed , the probability that something happen at a given _ point _ , which is described by the probability distribution in the _ comoving volume _ as discussed above , does not take into account the exponential growth of the _ physical volume _ of the universe , which is the main property of inflationary models . for the discussion of the volume - weighted probability measure , see e.g. @xcite .    despite this reservation , it is quite important to study the probability distribution in the comoving volume since its investigation is an essential step in several different attempts to calculate the probabilities in the landscape . here we would like to discuss another problem with the stationary probability distribution described above : this stationary distribution is not very useful during inflation when one can ignore the existence of the lowest ds space with @xmath335 ; in order to obtain a stationary distribution one should take into account the growth of the physical volume of the universe @xcite . moreover , as we are going to show now , this distribution does not necessarily apply to the string theory landscape , simply because in the kklt scenario there are no stable ds vacua that could serve as a ground state of the universe . metastability of ds space in the kklt construction was emphasized in @xcite and in many subsequent papers . here we would like to look at this issue in a more detailed way . stationarity of the probability distribution ( [ weinb2 ] ) was achieved because the lowest ds state did not have any further way to fall . meanwhile , in string theory it is always possible for the metastable ds state to decay . it is important that if it decays by production of the bubbles of 10d minkowski space @xcite , or by production of bubbles containing a collapsing open universe with a negative cosmological constant @xcite , the standard mechanism of return back to the original ds state does not operate any longer . these processes work like sinks for the flow of probability in the landscape . the fraction of the comoving volume in ds vacua will decrease in time .    to get a visual understanding of this process , assume that we are looking at the universe in the comoving coordinates , and paint black all of its parts corresponding to one of the two ds states , and paint white the parts in another ds state . then the stationary regime which we discussed so far would imply that eventually the whole universe in average becomes gray , and the level of gray asymptotically will remain constant . suppose now that some parts of the universe may tunnel to a state with a negative cosmological constant . these parts will collapse , so they will not return to the initial ds vacua . if we paint such parts red , then the universe , instead of reaching the constant gray level , eventually will look completely red . this does not mean that eventually all parts of our universe will end up in the sink . due to the exponential expansion of the universe in various ds vacua , the universe eternally rejuvenates , and the total volume of the universe in different ds states continues to grow exponentially @xcite . however , if one wants to find the ratio of the comoving volume of the universe in different states , then instead of solving the detailed balance equation ( [ balance ] ) , one should solve the `` vacuum dynamics '' equations @xmath448 @xmath449 here @xmath450 is the probability current from the lower ds vacuum to the sink , i.e. to a collapsing universe , or to a minkowski vacuum , @xmath451 is the probability current from the upper ds vacuum to the sink , @xmath452 is the probability current from the lower ds vacuum to the upper ds vacuum , and @xmath453 is the probability current from the upper ds vacuum to the lower ds vacuum . combining this all together , gives us the following set of equations for the probability distributions : @xmath454 @xmath455 since the entropy is inversely proportional to the energy density , the entropy of the lower level is higher , @xmath456 . since the tunneling is exponentially suppressed , we have @xmath457 , so we have a hierarchy @xmath458 . we will associate the lower vacuum with our present vacuum state , with @xmath459 , and therefore in what follows we will often have in mind the regime @xmath460 .    for simplicity , we will study here the possibility that only the lower vacuum can tunnel to the sink , i.e. we will take the limit @xmath461 and drop the term @xmath462 in eq . . we will analyze more general solutions of equations ( [ v1 ] ) , ( [ v2 ] ) in a separate publication @xcite . on the other hand , we will keep in mind the results of the previous section , where we have found that typically the probability of the decay of a metastable ds vacuum to a sink can be quite high , @xmath463 . therefore we expect that @xmath464 . other mechanisms of a relatively fast vacuum decay have been discussed e.g. in @xcite . the main conclusions to be obtained in this section will be valid for any of these mechanisms .    by solving equations ( [ v1 ] ) , ( [ v2 ] ) , one can show that @xmath465 and @xmath466 exponentially decrease in the course of time due to the existence of the sink ( the term @xmath467 ) , but the corresponding physical volume of the universe exponentially grows , and the ratio @xmath468 approaches a stationary regime @xmath469 . in order to find @xmath470 , one can add to each other our equations ( without the term @xmath471 ) . this yields @xmath472 the solution is @xmath473    as for the ( asymptotically ) constant ratio @xmath474 , from eqs . ( [ v1 ] ) , ( [ v3 ] ) one finds @xmath475    one may consider two interesting regimes . suppose first that @xmath476 , i.e. the probability to fall to the sink from the lower vacuum is smaller than the probability of the decay of the upper vacuum . in this case one recovers the previous result , eq . ( [ weinb2 ] ) , which is related to the square of the hartle - hawking wave function : @xmath477    now let us consider an opposite regime , and assume that the decay rate of the uplifted ds vacuum to the sink is relatively large , @xmath478 ( which automatically means that @xmath479 ) . in this case the solution of eq . ( [ ttt ] ) is @xmath480 i.e. one has an inverted probability distribution . this result has a simple interpretation : if the `` thermal exchange '' between the two ds vacua occurs very slowly as compared to the rate of the decay of the lower ds vacuum , then the main fraction of the volume of the ds vacua will be in the state with higher energy density . note that the total fraction of the comoving volume remaining in both of the ds spaces exponentially decreases in time due to the existence of the sink in the lower ds vacuum : @xmath481 this result may have interesting methodological implications . in many recent discussions of the probabilities in the landscape the authors consider the world as seen by an imaginary eternal observer . in the context of the stringy landscape such considerations are problematic , because every would - be eternal observer living in the landscape ( unlike an observer in an eternally existing ds space ) begins his life in a cosmological singularity and ends up in a sink , within the typical time determined by eq . ( [ sinkpower ] ) . moreover , if the sinks occur in the higher ds vacua ( see below ) , then the `` eternal observers '' , or markers , will have a good chance to die well before reaching the thermalized regions produced after inflation . introduction of a finite lifetime of such observers may result in the modification of the probability measure as suggested in @xcite .    here we should make an important clarification . the probability distribution @xmath431 describes the fraction of the volume of the universe in a particular ds state . however , when the bubbles of a new phase expand , their interior becomes an empty ds space . if we are usual observers which are born after reheating of the universe , then one may argue that the probability _ to live _ in the bubble @xmath482 is proportional not to the relative volume of the bubbles @xmath483 but to the frequency of their production , i.e. to the ratio of the probability currents @xmath484 @xcite . in the absence of the sink , the fraction of the comoving volume which flows to the lower ds vacuum due to the tunneling from the upper ds vacuum is equal to the fraction of the volume jumping upwards from the lowest vacuum to the higher vacuum . in other words , the two probability currents are exactly equal to each other , j_01 = j_10  , which is the essence of the detailed balance equation ( [ balance ] ) . interestingly , our new results imply that this regime remains approximately valid even in the presence of the sink , under the condition @xmath485 .    on the other hand , in the regime described by eq . ( [ ooo ] ) , which occurs if the decay rate to the sink is large enough , @xmath486 , one has a completely different result : @xmath487 thus we have a crucial regime change at the moment when the decay rate of the lower vacuum to the sink starts competing with the decay rate of the upper ds vacuum . now we will compare the probability current to the sink , @xmath488 , and the probability current from the lower ds upwards , @xmath489 . in the absence of the sink , @xmath488 was zero , so every point falling to the lower ds vacuum eventually jumps up to the higher ds vacua . this leads to eternal recycling of ds vacua @xcite , and to the problems discussed in @xcite : in the absence of the sinks ( i.e. in the models with eternal ds space ) each point in the comoving volume would indefinitely wonder between different ds spaces . in this scenario , it would be more probable that our part of the universe was formed as a result of a large quantum fluctuation upwards from an empty ds space , instead of being formed at the stage of inflation . this would leave the homogeneity of our universe , and the origin of perturbations of metric unexplained @xcite . invention of the kklt mechanism demonstrated that all ds vacua in string theory are metastable @xcite , with the lifetime smaller than the recurrence time @xmath490 . this means that we do not have much time for the incredibly improbable jumps upwards to happen @xcite . however , the time necessary for the improbable fluctuations upwards , which could re - create conditions for the existence of life , is also somewhat smaller than the recurrence time @xcite . therefore in order to make sure that the universe does not recycle back from our ds state we would need to double - check that the probability to jump upwards is indeed strongly suppressed as compared to the probability to tunnel to the sink . the estimates of the decay probability of ds vacuum contained in our paper confirm that this is indeed the case . indeed , investigation of the tunneling in several models discussed in our paper suggests that the probability that any given part of ds space with @xmath491 will jump upwards and participate in the process of recycling of the universe is suppressed by a factor of @xmath492 . this can be expressed in terms of the ratio of the probability current from the lower ds vacuum to the upper vacuum @xmath489 with the probability current from the lower ds vacuum to the sink @xmath488 : @xmath493 the last part of this equation follows from the condition @xmath460 , and from the estimate of the decay rate obtained in the previous section , which suggests that in many cases this rate is quite high , @xmath464 . despite this fact , recycling may still be possible , since for @xmath485 one has @xmath494 . on the other hand , for @xmath486 one has @xmath495 . this suggests that in this regime the probability that our universe was created due to the jumps up @xcite becomes exponentially small as compared with the probability that it was formed as a result of inflation . using our results , one can easily check that this is indeed the case if the rate of spontaneous formation of the parts of the universe of our type in the vacuum @xmath496 is smaller than the decay rate of the vacuum @xmath496 to the sink . for wide sinks , this condition can be easily satisfied . this conclusion is closely related to the inversion of the probability distribution for @xmath486 ( wide sink regime ) , see eq . ( [ ooo ] ) . note that inverted probability distributions in the presence of the sink represented by minkowski vacuum naturally appear in the simplest chaotic inflation scenario with the volume - weighted probability measures proposed in @xcite ; see figs . 10 and 12 in @xcite , where the basic theory of the volume - weighted distributions in the inflationary cosmology was developed . with the first of the two volume - weighted probability measures proposed in @xcite the problem of non - inflationary jumps up does not appear at all ; the growth of the total volume of the universe ( and the corresponding fluxes ) is predominantly determined by the exponential expansion of domains with larger vacuum energy and by the subsequent inflation on the way down . in our investigation so far we considered the simplest model with two ds vacua , with the probability leak in the lower ds vacuum . on the other hand , if the probability leak occurs only in the upper ds vacuum , one may ignore the term @xmath497 in eq . ( [ v1 ] ) but keep the term @xmath498 in eq . ( [ v2 ] ) . in this case , for @xmath499 one again recovers the distribution ( [ hheq ] ) , whereas for @xmath500 one finds @xmath501 note that in this regime we have a flat probability distribution , which does not depend on @xmath502 . more generally , one should consider a dynamical equilibrium of a system of many ds , minkowski and ads vacua , where the probability leaks may occur at each level . the fraction of the comoving volume occupied by different metastable vacua in such a universe will be exponentially sensitive to the decay modes for each of these vacua . this makes the investigation of the probabilities in the landscape much more involved , but also much more interesting . previously we were trying to construct the wave function of the ground state of the universe . now we should learn how to study the universe with the holes in the ground . there are two possible approaches to string theory landscape . the first one is based on the calculation of the wave function of the universe , which may describe all branches of the wave function corresponding to all possible vacua . in fact , the first two proposals for the anthropic solution of the cosmological problem in the context of inflationary cosmology were based on the investigation of quantum creation of the universe with various values of scalar fields and fluxes @xcite , or with different types of compactification @xcite . another approach is to consider an inflationary universe consisting of many parts corresponding to all possible vacua @xcite , which allows , in particular , to consider a single universe divided into parts corresponding to all possible values of scalar fields , fluxes and compactifications , and look for the anthropic solution of the cosmological constant problem in this context , see e.g. @xcite . conceptually , this approach is much simpler than the approach based on quantum cosmology , but in order to make it sufficiently general one must find a way to describe the universe consisting of many different parts and describe possible transitions between them . this was the goal of our paper . one of the problems here is that a fully consistent `` landscape action '' with the bulk and brane sources is not known @xcite . however , one can hope that such action can be constructed at least in the case of unbroken supersymmetry following the `` supersymmetry in singular spaces '' proposal @xcite . this strategy was already applied in examples of fully consistent supersymmetric bulk and brane actions in @xcite which include the d8 domain walls with the piece - wise constant @xmath503 flux . in a more general situation with moduli stabilization when various fluxes may jump due to brane sources , the corresponding construction has still to be developed . therefore in the current project we restricted ourselves to studies of the regions of the landscape where stabilization of the moduli is under control and fluxes are fixed . even in this restricted and simplified situation we were able to find multiple vacua and bps domain walls separating them . in the first part of this paper we concentrated on the stringy landscape prior to the uplifting and found many bps domain wall solutions separating different parts of the stringy landscape from each other . we did it using the methods developed earlier in @xcite , and also more recent methods based on the `` new attractor equations '' @xcite .    in the second part of the paper we studied what happens with the domain walls after the uplifting . in general , potentials of many moduli before and after uplifting may have many different ads minima of the moduli potential , for the same values of fluxes . before the uplifting , the tunneling between these minima was forbidden by supersymmetry . one could expect , therefore , that after the uplifting the rate of tunneling may be related to the degree of supersymmetry breaking . indeed , in the particular examples that we have studied the probability of the decay of the uplifted vacuum to the remaining ads vacuum was suppressed by a factor @xmath504 , where @xmath505 is the gravitino mass . this decay probability is much greater than @xmath506 . one of the interesting features of the decay of ds vacua to minkowski or ads vacua is that this decay is irreversible . this process is substantially different from the stationary `` thermal exchange '' between various ds vacua . the existence of the channels of the irreversible vacuum decay , which serve as sinks for the flow of probability , is a distinguishing new feature of the stringy landscape which deserves further investigation . as we show in section 6 of our paper , the existence of the sinks strongly affects the probability distributions in string cosmology . it is a pleasure to thank l. alvarez  gaum , r. bousso , t. banks , m.  cveti , f. denef , b. freivogel , g. horowitz , s. kachru , m. mario , d. martelli , e. silverstein , s. shenker , l. susskind and w. taylor for stimulating conversations . the work of a.c . is supported in part by the european community s marie curie actions , project no . mrtn - ct-2004 - 005104 ( forces universe ) . the work of r.k . , a.l . and a.g . was supported by nsf grant phy-0244728 . the work of a.g . was also supported by nsf grant phy-0097915 , and by doe under contracts de - ac03 - 76sf00515 and rfbr 05 - 01 - 00758 . linde , `` nonsingular regenerating inflationary universe , '' print-82 - 0554 , cambridge university preprint , 1982 , see http://www.stanford.edu/@xmath507alinde/1982.pdf ; a.d . linde , `` the new inflationary universe scenario , '' in : _ the very early universe _ gibbons , s.w . hawking and s.siklos , cambridge university press ( 1983 ) , pp . 205 - 249 , see http://www.stanford.edu/@xmath507alinde/1983.pdf ; a. d. linde , `` eternally existing self - reproducing chaotic inflationary universe , '' phys .  lett . b * 175 * , 395 ( 1986 ) . w.  lerche , d.  lst and a.  n.  schellekens , `` chiral four - dimensional heterotic strings from selfdual lattices , '' nucl . b * 287 * , 477 ( 1987 ) . r.  bousso and j.  polchinski , `` quantization of four - form fluxes and dynamical neutralization of the cosmological constant , '' jhep * 0006 * , 006 ( 2000 ) [ arxiv : hep - th/0004134 ] . m.  r.  douglas , `` the statistics of string / m theory vacua , '' jhep * 0305 * 046 ( 2003 ) [ arxiv : hep - th/0303194 ] ; f.  denef and m.  r.  douglas , `` distributions of flux vacua , '' jhep * 0405 * , 072 ( 2004 ) [ arxiv : hep - th/0404116 ] . m.  cvetic and h.  h.  soleng , `` supergravity domain walls , '' phys . rept .   * 282 * , 159 ( 1997 ) [ arxiv : hep - th/9604090 ] ; m.  cvetic , s.  griffies and s.  j.  rey , `` nonperturbative stability of supergravity and superstring vacua , '' nucl . b * 389 * , 3 ( 1993 ) [ arxiv : hep - th/9206004 ] ; m.  cvetic , s.  griffies and s.  j.  rey , `` static domain walls in n=1 supergravity , '' nucl . b * 381 * , 301 ( 1992 ) [ arxiv : hep - th/9201007 ] ; m.  cvetic , s.  griffies and h.  h.  soleng , `` local and global gravitational aspects of domain wall space - times , '' phys .  rev . d * 48 * , 2613 ( 1993 ) [ arxiv : gr - qc/9306005 ] . a.  ceresole , g.  dallagata , r.  kallosh and a.  van proeyen , `` hypermultiplets , domain walls and supersymmetric attractors , '' phys . d * 64 * , 104006 ( 2001 ) [ arxiv : hep - th/0104056 ] . k.  behrndt , g.  lopes cardoso and d.  lst , `` curved bps domain wall solutions in four - dimensional n = 2 supergravity , '' nucl . b * 607 * , 391 ( 2001 ) [ arxiv : hep - th/0102128 ] ; g.  lopes cardoso , g.  dallagata and d.  lst , `` curved bps domain wall solutions in five - dimensional gauged supergravity , '' jhep * 0107 * ( 2001 ) 026 [ arxiv : hep - th/0104156 ] ; g.  lopes cardoso , g.  dallagata and d.  lst , `` curved bps domain walls and rg flow in five dimensions , '' jhep * 0203 * ( 2002 ) 044 [ arxiv : hep - th/0201270 ] ; k.  behrndt and m.  cvetic , `` bent bps domain walls of d = 5 n = 2 gauged supergravity coupled to hypermultiplets , '' phys . d * 65 * , 126007 ( 2002 ) [ arxiv : hep - th/0201272 ] . j.  louis and s.  vaula , `` n=1 domain wall solutions of massive type ii supergravity as generalized geometries , '' arxiv : hep - th/0605063 . d.  z.  freedman , s.  s.  gubser , k.  pilch and n.  p.  warner , `` renormalization group flows from holography supersymmetry and a c - theorem , '' adv . * 3 * , 363 ( 1999 ) [ arxiv : hep - th/9904017 ] . r.  kallosh and a.  linde , `` landscape , the scale of susy breaking , and inflation , '' jhep * 0412 * , 004 ( 2004 ) [ arxiv : hep - th/0411011 ] ; j.  j.  blanco - pillado , r.  kallosh and a.  linde , `` supersymmetry and stability of flux vacua , '' jhep * 0605 * , 053 ( 2006 ) [ arxiv : hep - th/0511042 ] . g.  t.  horowitz and a.  strominger , `` counting states of near - extremal black holes , '' phys . lett .   * 77 * , 2368 ( 1996 ) [ arxiv : hep - th/9602051 ] ; c.  g.  .  callan and j. m.  maldacena , `` d - brane approach to black hole quantum mechanics , '' nucl . b * 472 * , 591 ( 1996 ) [ arxiv : hep - th/9602043 ] ; a.  w.  peet , `` tasi lectures on black holes in string theory , '' arxiv : hep - th/0008241 . d.  z.  freedman , c.  nunez , m.  schnabl and k.  skenderis , `` fake supergravity and domain wall stability , '' phys . d * 69 * , 104027 ( 2004 ) [ arxiv : hep - th/0312055 ] ; a.  celi , a.  ceresole , g.  dallagata , a.  van proeyen and m.  zagermann , `` on the fakeness of fake supergravity , '' phys . d * 71 * , 045009 ( 2005 ) [ arxiv : hep - th/0410126 ] . r.  kallosh , `` new attractors , '' jhep * 0512 * , 022 ( 2005 ) [ arxiv : hep - th/0510024 ] . g.  dallagata , `` non - kaehler attracting manifolds , '' jhep * 0604 * ( 2006 ) 001 [ arxiv : hep - th/0602045 ] . j.  p.  hsu , a.  maloney and a.  tomasiello , `` black hole attractors and pure spinors , '' arxiv : hep - th/0602142 . a.  giryavets , `` new attractors and area codes , '' jhep * 0603 * , 020 ( 2006 ) [ arxiv : hep - th/0511215 ] . d.  lst , s.  reffert , w.  schulgin and s.  stieberger , `` moduli stabilization in type iib orientifolds . i : orbifold limits , '' arxiv : hep - th/0506090 . k.  behrndt and g.  dallagata , `` vacua of n = 2 gauged supergravity derived from non - homogenous quaternionic spaces , '' nucl . b * 627 * , 357 ( 2002 ) [ arxiv : hep - th/0112136 ] . m.  graa , `` flux compactifications in string theory : a comprehensive review , '' phys . rept .   * 423 * , 91 ( 2006 ) [ arxiv : hep - th/0509003 ] . a.  giryavets , s.  kachru and p.  k.  tripathy , `` on the taxonomy of flux vacua , '' jhep * 0408 * , 002 ( 2004 ) [ arxiv : hep - th/0404243 ] . s.  ashok and m.  r.  douglas , `` counting flux vacua , '' jhep * 0401 * , 060 ( 2004 ) [ arxiv : hep - th/0307049 ] . o.  dewolfe , a.  giryavets , s.  kachru and w.  taylor , `` enumerating flux vacua with enhanced symmetries , '' jhep * 0502 * , 037 ( 2005 ) [ arxiv : hep - th/0411061 ] . a.  klemm and s.  theisen , `` considerations of one modulus calabi - yau compactifications : picard - fuchs equations , kahler potentials and mirror maps , '' nucl .  phys . b * 389 * , 153 ( 1993 ) [ arxiv : hep - th/9205041 ] . f.  denef , b.  r.  greene and m.  raugas , `` split attractor flows and the spectrum of bps d - branes on the quintic , '' jhep * 0105 * , 012 ( 2001 ) [ arxiv : hep - th/0101135 ] . r.  bousso , b.  freivogel and m.  lippert , `` probabilities in the landscape : the decay of nearly flat space , '' arxiv : hep - th/0603105 . s.  weinberg , `` does gravitation resolve the ambiguity among supersymmetry vacua ? , '' phys . * 48 * ( 1982 ) 1776 . s.  j.  parke , `` gravity , the decay of the false vacuum and the new inflationary universe scenario , '' phys . b * 121 * , 313 ( 1983 ) . s.  w.  hawking and i.  g.  moss , `` supercooled phase transitions in the very early universe , '' phys . b * 110 * , 35 ( 1982 ) . a.  a.  starobinsky , `` stochastic de sitter ( inflationary ) stage in the early universe , '' in : _ current topics in field theory , quantum gravity and strings _ , lecture notes in physics , eds . de vega and n. sanchez ( springer , heidelberg 1986 ) * 206 * , p. 107 . a.  d.  linde , `` hard art of the universe creation ( stochastic approach to tunneling and baby universe formation ) , '' nucl . b * 372 * , 421 ( 1992 ) [ arxiv : hep - th/9110037 ] . a.  d.  linde , d.  a.  linde and a.  mezhlumian , `` from the big bang theory to the theory of a stationary universe , '' phys . d * 49 * , 1783 ( 1994 ) [ arxiv : gr - qc/9306035 ] . k.  m.  lee and e.  j.  weinberg , `` decay of the true vacuum in curved space - time , '' phys . d * 36 * , 1088 ( 1987 ) . j.  garriga and a.  vilenkin , `` recycling universe , '' phys . d * 57 * , 2230 ( 1998 ) [ arxiv : astro - ph/9707292 ] . n.  goheer , m.  kleban and l.  susskind , `` the trouble with de sitter space , '' jhep * 0307 * , 056 ( 2003 ) [ arxiv : hep - th/0212209 ] . a.  d.  linde , `` quantum creation of the inflationary universe , '' lett . nuovo cim .   * 39 * , 401 ( 1984 ) ; a.  vilenkin , `` quantum creation of universes , '' phys . d * 30 * , 509 ( 1984 ) . j.  garcia - bellido , a.  d.  linde and d.  a.  linde , `` fluctuations of the gravitational constant in the inflationary brans - dicke cosmology , '' phys . d * 50 * , 730 ( 1994 ) [ arxiv : astro - ph/9312039 ] ; a.  vilenkin , `` predictions from quantum cosmology , '' phys . lett .   * 74 * , 846 ( 1995 ) [ arxiv : gr - qc/9406010 ] ; j.  garcia - bellido and a.  d.  linde , `` stationarity of inflation and predictions of quantum cosmology , '' phys . d * 51 * , 429 ( 1995 ) [ arxiv : hep - th/9408023 ] ; j.  garcia - bellido and a.  d.  linde , `` stationary solutions in brans - dicke stochastic inflationary cosmology , '' phys .  rev . d * 52 * , 6730 ( 1995 ) [ arxiv : gr - qc/9504022 ] . a.  vilenkin , `` making predictions in eternally inflating universe , '' phys . d * 52 * , 3365 ( 1995 ) [ arxiv : gr - qc/9505031 ] ; s.  winitzki and a.  vilenkin , `` uncertainties of predictions in models of eternal inflation , '' phys . d * 53 * , 4298 ( 1996 ) [ arxiv : gr - qc/9510054 ] ; j.  garriga and a.  vilenkin , `` a prescription for probabilities in eternal inflation , '' phys . d * 64 * , 023507 ( 2001 ) [ arxiv : gr - qc/0102090 ] ; j.  garriga , d.  schwartz - perlov , a.  vilenkin and s.  winitzki , `` probabilities in the inflationary multiverse , '' jcap * 0601 * , 017 ( 2006 ) [ arxiv : hep - th/0509184 ] ; r.  easther , e.  a.  lim and m.  r.  martin , `` counting pockets with world lines in eternal inflation , '' jcap * 0603 * , 016 ( 2006 ) [ arxiv : astro - ph/0511233 ] ; a.  vilenkin , `` probabilities in the landscape , '' arxiv : hep - th/0602264 . a.  h.  guth , `` the inflationary universe : a possible solution to the horizon and flatness problems , '' phys . d * 23 * , 347 ( 1981 ) . p.  j.  steinhardt , `` natural inflation , '' in : _ the very early universe _ gibbons , s.w . hawking and s.siklos , cambridge university press , ( 1983 ) . a.  r.  frey , m.  lippert and b.  williams , `` the fall of stringy de sitter , '' phys . d * 68 * , 046008 ( 2003 ) [ arxiv : hep - th/0305018 ] . d.  green , e.  silverstein and d.  starr , `` attractor explosions and catalyzed vacuum decay , '' arxiv : hep - th/0605047 . a.  d.  linde , `` the inflationary universe , '' rept . phys .   * 47 * 925 ( 1984 ) . a. d. sakharov , `` cosmological transitions with a change in metric signature , '' sov . phys . jetp * 60 * , 214 ( 1984 ) [ zh . . fiz . * 87 * , 375 ( 1984 ) ] . j.  d.  brown and c.  teitelboim , `` neutralization of the cosmological constant by membrane creation , '' nucl .  phys . b * 297 * , 787 ( 1988 ) ; j.  d.  brown and c.  teitelboim , `` dynamical neutralization of the cosmological constant , '' phys . b * 195 * , 177 ( 1987 ) . j.  l.  feng , j.  march - russell , s.  sethi and f.  wilczek , `` saltatory relaxation of the cosmological constant , '' nucl . b * 602 * , 307 ( 2001 ) [ arxiv : hep - th/0005276 ] . t.  banks , m.  dine and e.  gorbatov , `` is there a string theory landscape ? , '' jhep * 0408 * , 058 ( 2004 ) [ arxiv : hep - th/0309170 ] . f.  denef , private discussion . de alwis , `` transitions between flux vacua , '' arxiv : hep - th/0605184 . e.  bergshoeff , r.  kallosh and a.  van proeyen , `` supersymmetry in singular spaces , '' jhep * 0010 * , 033 ( 2000 ) [ arxiv : hep - th/0007044 ] . e.  bergshoeff , r.  kallosh , t.  ortin , d.  roest and a.  van proeyen , `` new formulations of d = 10 supersymmetry and d8 - o8 domain walls , '' class . grav .   * 18 * , 3359 ( 2001 ) [ arxiv : hep - th/0103233 ] ; r.  kallosh , s.  prokushkin and m.  shmakova , `` domain walls with strings attached , '' jhep * 0107 * , 023 ( 2001 ) [ arxiv : hep - th/0107097 ] .
Abstract:

Answer: we develop a theory of static bps domain walls in stringy landscape and present a large family of bps walls interpolating between different supersymmetric vacua . examples include kklt models , stu models , type iib multiple flux vacua , and models with several minkowski and ads vacua . after the uplifting , some of the vacua become ds , whereas some others remain ads . the near - bps walls separating these vacua may be seen as bubble walls in the theory of vacuum decay . as an outcome of our investigation of the bps walls , we found that the decay rate of ds vacua to a collapsing space with a negative vacuum energy can be quite large . the parts of space that experience a decay to a collapsing space , or to a minkowski vacuum , never return back to ds space . the channels of irreversible vacuum decay serve as sinks for the probability flow . the existence of such sinks is a distinguishing feature of the landscape . we show that it strongly affects the probability distributions in string cosmology .     + su - itp-2006 - 17 +   + may 29 , 2006    _ @xmath0 infn @xmath1 dipartimento di fisica teorica + universita di torino , via pietro giuria 1 , 10125 torino , italy _    _ @xmath2 physics department , theory unit , cern , + ch 1211 , geneva 23 , switzerland _    5.5 mm