Supersymmetry preserving and breaking degenerate vacua, and radiative moduli stabilization
Tatsuo Kobayashi, Naoya Omoto, Hajime Otsuka, and Takuya H. Tatsuishi

TL;DR
This paper introduces a novel moduli stabilization method where supersymmetric and breaking vacua are initially degenerate, with loop corrections resolving this degeneracy and predicting light axions.
Contribution
It presents a new scenario for moduli stabilization with degenerate vacua that are distinguished by loop corrections, leading to potential light axions.
Findings
Degenerate vacua are stabilized by loop corrections.
Light axions are predicted in the model.
Loop effects resolve initial vacuum degeneracy.
Abstract
We propose a new type of moduli stabilization scenario where the supersymmetric and supersymmetry-breaking minima are degenerate at the leading level. The inclusion of the loop-corrections originating from the matter fields resolves this degeneracy of vacua. Light axions are predicted in one of our models.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
EPHOU-17-006
WU-HEP-17-08
**Supersymmetry preserving and breaking degenerate vacua, and radiative moduli stabilization
**
Tatsuo Kobayashi,1,*** E-mail address: [email protected] Naoya Omoto,1,††† E-mail address: [email protected] Hajime Otsuka,2,‡‡‡ E-mail address: [email protected]
and Takuya H. Tatsuishi1§§§ E-mail address: [email protected]
1Department of Physics, Hokkaido University, Sapporo 060-0810, Japan
2*Department of Physics, Waseda University, Tokyo 169-8555, Japan *
( Abstract
We propose a new type of moduli stabilization scenario where the supersymmetric and supersymmetry-breaking minima are degenerate at the leading level. The inclusion of the loop-corrections originating from the matter fields resolves this degeneracy of vacua. Light axions are predicted in one of our models.
)
1 Introduction
The superstring theory predicts six-dimensional (6D) compact space in addition to four-dimensional (4D) spacetime. The size and shape of the 6D compact space are determined by moduli. Thus, moduli are a characteristic feature in the superstring theory on compact space. Indeed, moduli fields play important roles in the superstring theory and its 4D low-energy effective field theory, in particular in particle phenomenology and cosmology.(See for a review, e.g. Refs. [1, 2].) Studies on physics relevant to moduli would provide a remnant of the superstring theory on compact space.
Gauge couplings, Yukawa couplings and other couplings in the 4D low-energy effective field theory are given by vacuum expectation values (VEVs) of moduli fields. In the phenomenological point of view, the spectrum of supersymmetric particles is sensitive to the supersymmetry (SUSY) breaking of moduli fields through the gravitational interactions between the moduli fields and matter fields [3, 4, 5, 6]. This aspect would be relevant to dark matter physics, since the lightest superpartner is a candidate of dark matter. Moreover, the thermal history of the Universe highly depends on the dynamics of moduli fields as well as axion fields, which are imaginary parts of moduli fields. From the theoretical point of view, these moduli fields are originating from the vector and tensor fields in low-energy effective action of the superstring theory as well as the higher-dimensional supergravity. The stabilization of the moduli field is one of the most important issues to realize a consistent low-energy effective action of the superstring theory.
The moduli potential is prohibited by the higher-dimensional gauge and Lorentz symmetries at the perturbative level. On the other hand, the nontrivial background fields and nonperturbative effects generate the moduli potential. Then, one can stabilize the moduli fields and also study the dynamics relevant to moduli fields. The vacuum structure of the moduli potential is of particular importance. For example, the flat direction of the moduli potential can drive cosmological inflation111See for the detail of moduli inflations as well as axion inflations, e.g., Ref. [7]. and the lifetime of our Universe depends on the (meta)stability of the vacuum.
So far, there are several mechanisms to stabilize the moduli, in particular closed string moduli fields, e.g., the Kachru-Kallosh-Linde-Trivedi (KKLT) scenario [8] and the LARGE volume scenario [9]. (See, e.g., Ref. [1] and references therein.) In this paper, we propose a new type of moduli stabilization scenario by using the string-derived four-dimensional supergravity action. We find that the supersymmetric and SUSY-breaking vacua are degenerate at the tree level and they are independent of the term of certain moduli fields. The loop effects originating from the matter fields generate the moduli potential and resolve this degeneracy of vacua.
This paper is organized as follows. In section 2, we study a new type of potential for a complex structure modulus within the framework of the type IIB superstring theory. Its minimum has the SUSY preserving and breaking degenerate vacua, which can be resolved by loop effects. Then, the modulus field including the axion is stabilized. In section 3, we study similar aspects for the Kähler modulus. In this case, the real part of the Kähler modulus is stabilized, but the axion remains massless at this stage. Finally, Sec. 4 is devoted to the conclusion. In the Appendix, we give an illustrating model for the moduli stabilization.
2 Complex structure moduli
In sections 2 and 3, we consider two illustrative supergravity models where the moduli fields correspond to the complex structure modulus and Kähler modulus within the framework of the type IIB superstring theory. In both scenarios, we show that the supersymmetric and SUSY-breaking vacua are degenerate at the leading level. The inclusion of the one-loop corrections from the matter fields resolves this degeneracy.
In the type IIB superstring theory on the Calabi-Yau (CY) orientifold, the Kähler potential of moduli fields is described by
[TABLE]
where is the holomorphic three-form of the CY manifold and is the volume of the CY manifold. Here, , and denote the dilaton, complex structure moduli and Kähler moduli, respectively.
The three-form flux can induce the superpotential [10]
[TABLE]
where is an imaginary self-dual three-form. Also nonperturbative effects such as D-brane instantons and gaugino condensations can generate the superpotential of and , e.g.,
[TABLE]
where represent the -dependent one-loop corrections and and are the numerical constants.
2.1 The degenerate scalar potential
First, let us consider the following Kähler potential and superpotential based on the four-dimensional supergravity:
[TABLE]
where are the complex constants. Throughout this paper, we use the reduced Planck unit . In the type IIB superstring theory on the CY orientifold, the modulus field could be identified with one of the complex structure moduli of the CY manifold. We now assume that the other complex structure moduli and dilaton are stabilized at the minimum by the three-form fluxes [11]. When the parameters are determined only by the three-form fluxes as well as VEVs of the other complex structure moduli, these would be of . On the other hand, when the above superpotential appears from the instanton effects, are characterized as , with being a certain Kähler modulus of the CY manifold, and could be suppressed. In the following analysis, we treat as complex constants by further assuming that all the Kähler moduli are stabilized at the minimum by the other nonperturbative effects. When are of the order of unity, it is a challenging issue to stabilize the Kähler modulus at the scale above the mass of . Conversely, when both are exponentially suppressed, one can achieve the above assumption as shown later. In the following, we use the parametrization as and .
To see the degenerate supersymmetric and SUSY-breaking vacua, we calculate the scalar potential in the notation of and ,
[TABLE]
where is the inverse of Kähler metric and
[TABLE]
As a result, the scalar potential (5) remains flat in the direction of . On the other hand, from the extremal condition of ,
[TABLE]
is stabilized at the minimum
[TABLE]
where should be negative to justify our low-energy effective action. Note that is typically of the order of unity in both scenarios: and . If is smaller than unity, the moduli space of deviates from the large complex structure regime and our discussing logarithmic the Kähler potential is not reliable. It turns out that the mass squared of canonically normalized ,
[TABLE]
is positive and the vacuum energy,
[TABLE]
is negative at this minimum. Note that the mass squared of canonically normalized at this vacuum is taken smaller than the other moduli fields, in particular, the Kähler moduli. For example, when we consider the nonperturbative superpotential for the Kähler moduli irrelevant to our focusing , the Kähler moduli can be stabilized at certain minima by them and the mass squared of the overall Kähler modulus is given by ) for the KKLT scenario. It is much heavier than that of our considered as discussed in the Appendix A.
Although the vacuum energy is independent of , the -term of the modulus depends on ,
[TABLE]
where
[TABLE]
Thus, we find that the SUSY is preserved at and broken at , respectively. From the fact that the scalar potential is independent of , supersymmetric and SUSY-breaking vacua are degenerate.
2.2 Loop corrections
The analysis in Sec. 2.1 shows that the vacuum energy is negative, , and at the same time, the scalar potential remains flat in the direction of . First, we introduce the uplifting sector to obtain the tiny cosmological constant. In particular, we assume that the -independent potential induced by the anti-D-branes uplifts the anti-de Sitter vacuum to the Minkowski one such as the KKLT scenario.222Uplifting by spontaneous -term SUSY breaking is also possible [12, 13]. Even in that case, the direction would remain flat. Next, the nonvanishing term of gives rise to the soft terms of matter fields. These massive supersymmetric particles induce the -dependent scalar potential through one-loop corrections [14]:
[TABLE]
In the following analysis, we illustrate how we can stabilize by one-loop effects. For such a purpose, we focus on the situation that the gauginos and supersymmetric scalar fields contribute to the one-loop corrections. The gaugino masses are provided by
[TABLE]
where with representing the gauge kinetic functions for the standard model gauge groups, respectively. On the other hand, the soft scalar masses are given by
[TABLE]
where
[TABLE]
with being the Kähler metric of the matter fields.
To simplify our illustrating analysis, we assume that the typical soft scalar mass and gaugino mass mainly contribute to the one-loop potential. Then, those are characterized as
[TABLE]
where and are real constants and denotes the soft scalar mass induced by the -independent -term contributions. The gaugino mass may also have another contribution such as . Even in such a case, the following discussion is similar. For simple illustration, we restrict ourselves to the above spectrum of superpartners.
By rescaling the term of , one can set , and the corresponding one-loop potential can be written by
[TABLE]
where , and correspond to the multiplicities of the scalars and gauginos, respectively. Now, we include the constant coming from the terms of Kähler moduli and anti-D-brane effects to achieve the tiny cosmological constant at the vacuum. By using , , and , the above scalar potential is simplified as
[TABLE]
The first derivative of the one-loop potential with respect to is given by
[TABLE]
from which there are two possible minima leading to and . To see the nonvanishing , we draw the one-loop scalar potential as a function of by setting the following illustrative parameters:
[TABLE]
in the left (right) panel in Fig. 1 and
[TABLE]
in Fig. 2. It turns out that the nonvanishing of depends on the sign of and the value of .
To discuss the vacuum structure of the one-loop scalar potential, we analytically derive the condition of vanishing . From the second derivative of the one-loop potential with respect to
[TABLE]
becomes negative at the origin for
[TABLE]
It implies that when and
[TABLE]
the minimum of is taken as because of the instability at . However, we assume that the other moduli fields are decoupled from our system and such a high-scale SUSY breaking is not reliable. Thus, when , supersymmetric minimum is favorable. In this case, the soft terms are determined by terms of the Kähler moduli. Such a vanishing term of is also interesting from the aspects of the flavor structure of the matter fields. Indeed, Yukawa couplings among the standard model particles depend on the complex structure moduli through the compactification of an extra dimension as derived in the type IIB superstring theory with magnetized D-branes [15]. Thus, a sizable term of a complex structure modulus is dangerous for the flavor-changing processes among the supersymmetric particles, which are severely constrained in the low-scale SUSY-breaking scenario.
Let us take a closer look at the vacuum structure of the one-loop scalar potential. First, we further simplify the scalar potential by setting
[TABLE]
Then, the total scalar potential consist of the radiative corrections from bosonic and fermionic parts:
[TABLE]
where
[TABLE]
Note that gives rise to the tachyonic soft scalar masses from the definition of soft scalar mass in Eq. (17) and . Thus, is constrained to the range, .
First of all, let us discuss the bosonic contribution. From the first derivative of with respect to ,
[TABLE]
we find that there are two possible SUSY and SUSY-breaking vacua in the physical domain . When , the SUSY is preserved as
[TABLE]
whereas the SUSY is broken under ,
[TABLE]
It can be confirmed in Fig. 3, where we show the parameter dependence of the bosonic part .
Next, we focus on the fermionic contribution. From the first derivative of with respect to ,
[TABLE]
we find that there are two minima
[TABLE]
and
[TABLE]
However, the latter one leads to the tachyonic mass squared, . Thus the fermionic contribution leads to the SUSY-preserving minimum without depending on the value of . In Fig. 4, the fermionic part is drawn as a function of . It turns out that larger than causes the instability of the vacuum.
Now, we combine the bosonic and fermionic contributions. When , the bosonic part is the monotonically increasing function with respect to , whereas the fermionic one depends on the magnitude of . As shown in the left and right panels in Fig. 5, the vacuum instability arises by large values of included in the fermionic contribution. Since the vacuum energy is provided by
[TABLE]
one can achieve the tiny cosmological constant by setting .
On the other hand, in the case of , the vacuum structure is determined by . Indeed, when , the vacuum of is shifted to the origin by the fermionic contribution as shown in the left panel in Fig. 6. However, the large value of triggers the instability of the vacuum as can be seen in the right panel in Fig. 6.
So far, we have discussed the case with . It is found that the bosonic contribution gives the source of SUSY breaking, whereas the fermionic contribution preserves the SUSY irrelevant to the parameters. However, a large fermionic contribution causes the instability of the vacuum. In a similar way, we can discuss the case of , which changes the bosonic part of the radiative corrections as .
Finally, we discuss the possibility of nonvanishing by adding the nonperturbative effects of into the one-loop scalar potential. The nonperturbative effects of are expected to appear through e.g. the gaugino condensation on hidden D-branes where the gauge kinetic function involves the -dependent one-loop corrections [16]. Since the real part of is already stabilized by the superpotential in Eq. (4), the potential of can be extracted as
[TABLE]
where is the typical decay constant and is a real constant. In the following, we set for simplicity. Then, the term of ,
[TABLE]
leads to the following total scalar potential:
[TABLE]
Here, the constant is inserted to realize the tiny cosmological constant at the vacuum in a way similar to the previous scenario. By rescaling the parameters as
[TABLE]
we analyze the following potential
[TABLE]
Figure. 7 shows that the nonvanishing is achieved even when is positive. Since the origin of and are the nonperturbative effects, one can realize the low-scale SUSY-breaking scenario in this model.
We have assumed that loop corrections are dominant in the potential of . When other nonperturbative effects are dominant, obviously is stabilized by such nonperturbative effects and loop effects provide subdominant corrections.
3 Kähler moduli
In this section, we consider another example where the supersymmetric and SUSY-breaking minima are degenerate at the leading level.
By use of the flux-induced superpotential (2), the -term scalar potential is calculated as
[TABLE]
where is canceled by the no-scale structure of the Kähler moduli,
[TABLE]
Note that the above no-scale structure is valid only at the tree level.
Then, the dilaton and complex structure moduli are stabilized at the minimum,
[TABLE]
which lead to the Minkowski minimum . When , the supersymmetry is broken by the term of the Kähler moduli. In contrast to the previous section, we now assume that all the complex structure moduli and dilaton are stabilized by the flux-induced superpotential. Although the terms of and vanish at this Minkowski minimum, the terms of the Kähler moduli are nonvanishing, in general:
[TABLE]
when .
For simplicity, we study the model with the overall Kähler modulus with the CY volume . Then, the term of the Kähler modulus is simplified as
[TABLE]
Thus, supersymmetric and SUSY-breaking minima are also degenerate in a way similar to the previous section, since the scalar potential is independent of and . However, the supersymmetric vacuum corresponds to , that is, the decompactification limit.
When the leading corrections are involved, the Kähler potential of the Kähler modulus is corrected as [17]
[TABLE]
where with and being the Euler characteristic of CY and string coupling. These -corrections break the no-scale structure, and the scalar potential is generated as
[TABLE]
The sign of depends on the number of complex structure moduli and Kähler moduli. When the number of Kähler moduli is smaller than that of complex structure moduli, is positive. In the case of single Kähler modulus, the -term potential reduces to
[TABLE]
Along the same step outlined in Sec. 2, we take into account the loop corrections originating from the supersymmetric particles whose soft terms are dominated by the term of the Kähler modulus. It is remarkable that the loop corrections give rise to the stabilization of unlike the case in Sec. 2. Then, by assuming that the typical gaugino and supersymmetric scalar fields mainly contribute to the loop effects, the total scalar potential becomes
[TABLE]
where . Here, we employ the same notation of Sec. 2 and is chosen as a real constant, for simplicity.
By setting the illustrative parameters
[TABLE]
the scalar potential is drawn as in Fig. 8. As a result, the degeneracy of vacua is resolved by the loop corrections. In contrast to the discussion in Sec. 2, the vanishing corresponds to the unphysical domain . Thus, the SUSY-breaking vacuum is selected. Indeed, the above illustrative parameters give rise to the high-scale SUSY-breaking minimum, where the vacuum expectation value of ,
[TABLE]
resides in a reliable range of the supergravity approximation. After canonically normalizing the modulus,
[TABLE]
with , its mass squared is evaluated as
[TABLE]
in the reduced Planck unit, i.e. , which should be smaller than the other complex structure moduli and dilaton to justify our low-energy effective action. Since those complex structure moduli and dilaton fields have been stabilized at the SUSY-breaking minimum, their masses are typically greater than or equal to the gravitino mass in our numerical example. Our situation is thus justified.
Interestingly, the tuning of allows us to consider the tiny cosmological constant. In the above scenario, can be stabilized at a fine value, but its imaginary part, i.e. the axion, remains massless.
Let us take a closer look at the scalar potential. To simplify our analysis, we redefine the scalar potential as a function of which is constrained to be to avoid the tachyonic soft scalar masses,
[TABLE]
where
[TABLE]
and the corresponding parameters are defined as
[TABLE]
The simplified scalar potential is similar to the scalar potential in Eq. (27). However, in the current case, there exists a cubic term of and the cosmological constant is not included. As outlined in Sec. 2, we concentrate on the case with . When , the scalar potential is drawn as in the left panel in Fig. 9, where we take the same parameters as in Fig. 5. The vacuum is settled into . On the other hand, in the case of , the cubic term of lifts the potential compared with the case in Sec. 2. In the right panel in Fig. 9, we plot the scalar potential by taking the same parameters in Fig. 6.
Thus, the bosonic contribution gives the source of SUSY breaking, and the fermionic contribution gives the instability of the vacuum.
So far, we have studied the model with the overall Kähler modulus. However, we can discuss the model with many Kähler moduli fields . Even in such a model, the tree-level scalar potential is flat along all of the Kähler moduli directions because of the no-scale structure (42). Also, the tree-level potential is independent of terms of , and these themselves depend only on but not . In such multi-moduli models, gaugino masses and soft scalar masses would be written by
[TABLE]
Using these, we obtain the one-loop potential . Then, one can stabilize and in a similar way. However, all the axionic parts of remain massless at this stage.
Let us consider the axion potentials by adding the nonperturbative effects such as the gaugino condensation on hidden D-branes or world-sheet instanton effects. Those nonperturbative effects appear in the superpotential as shown in Eq. (3), and the potentials of the canonically normalized axions are extracted as
[TABLE]
where , , and are the real constants determined by the stabilization of the complex structure moduli, dilaton, and the real parts of the Kähler moduli, respectively. If the axions appear in the above nonperturbative effects, their masses are much smaller than the other moduli fields. Light axions would be interesting. For example, they can be candidates for the QCD axion and dark matter. In such a case, the QCD axion appears only in the gauge kinetic functions of the visible sector and not in the nonperturbative effects (3). To identify the Kähler axion with the QCD axion and dark matter, the axion decay constant is relatively larger than the desirable one. For that reason, we require certain dilution effects below the QCD scale. If the axion potentials are generated as in Eq. (58), the axion may derive inflationary expansion of the Universe by using the alignment mechanism [18]. Furthermore, the string axiverse scenario may be realized if they are light [19].333See for a realization of light axions, e.g., Refs. [20, 21]. At any rate, these aspects of axion phenomenology and cosmology highly depend on details of nonperturbative effects (58) below the energy scale, where the other moduli are stabilized. Such an analysis is beyond our scope of this paper. The detailed discussion of axion phenomenology will be left for a future work.
4 Conclusion
We have studied a new type of moduli potential and stabilization. In the model with the complex structure modulus , the supersymmetric and SUSY-breaking minima are degenerate at the leading order. That is, the tree-level potential is independent of the term of , but depends on . The term itself depends on . Loop effects due to -dependent gaugino and sfermion masses resolve the degeneracy of vacua and stabilize the axion . The SUSY vacuum or SUSY-breaking vacuum is selected depending on parameters in the potential. Low-scale SUSY breaking is also possible when additional proper nonperturbative effects are involved.
We have also studied the model with the Kähler modulus . This model has the flat direction along both and at the leading level. The SUSY vacuum and SUSY-breaking vacuum are degenerate, but the SUSY vacuum corresponds to the decompactification limit . In this model, the modulus term depends only on . The real part can be stabilized by inclusion of corrections and loop effects due to -dependent gaugino and sfermion masses. However, the axion remains massless at this stage.
We can extend the model with the single Kähler modulus to the models with many Kähler moduli. Their real parts can be stabilized in a similar mechanism, but many axionic parts would remain light. Such axions would be interesting, e.g., for candidates of dark matter and the QCD axion. Also, one of the light axions could derive the cosmological inflation if a proper potential is generated. Moreover, these axions would be interesting from the viewpoint of a string axiverse. Such axion phenomenology would be studied elsewhere.
Acknowledgments
T.K. was supported in part by the Grant-in-Aid for Scientific Research No. 26247042 and No. 17H05395 from the Ministry of Education, Culture, Sports, Science and Technology in Japan. H.O. was supported in part by a Grant-in-Aid for Young Scientists (B) (No. 17K14303) from Japan Society for the Promotion of Science.
Appendix A Moduli stabilization in the KKLT scenario
In this Appendix, we present the stabilization of complex structure moduli and Kähler moduli in the setup of Sec. 2. To see the typical mass scale of complex structure moduli and Kähler moduli, we consider one of the complex structure moduli and overall Kähler modulus . In D effective supergravity, those are stabilized by the nonperturbative superpotential of the form
[TABLE]
and the Kähler potential is given by
[TABLE]
Here, in Eq. (4) corresponds to , and we assume that is a real constant for simplicity, although it is a function of , in general.444When is a function of , the overall Kähler modulus could be stabilized at the racetrack minimum. Then, the total scalar potential becomes
[TABLE]
where , , are the Kähler covariant derivatives and is the inverse of Kähler metric .
We first stabilize the overall Kähler modulus by the superpotential
[TABLE]
which is the same as the KKLT scenario. It can be stabilized at the minimum satisfying
[TABLE]
which leads to the stabilization of as
[TABLE]
in the limit of . At this minimum, the mass squared of canonically normalized is evaluated as
[TABLE]
For example, when we obtain .
Let us next stabilize at the fixed minimum of . In this case, by setting and , the superpotential is the same as in Eq. (4) except for the factor . It turns out that the mass squared of is given by
[TABLE]
which is much lighter than in the limit of . To obtain the correct vacua of both the moduli fields, we have to stabilize them simultaneously. However, our analysis is justified because of the mass hierarchy between them.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] L. E. Ibanez and A. M. Uranga, “String theory and particle physics: An introduction to string phenomenology,” Cambridge University Press (2012).
- 2[2] R. Blumenhagen, B. Kors, D. Lust and S. Stieberger, Phys. Rept. 445 , 1 (2007) [hep-th/0610327].
- 3[3] V. S. Kaplunovsky and J. Louis, Phys. Lett. B 306 , 269 (1993) [hep-th/9303040].
- 4[4] A. Brignole, L. E. Ibanez and C. Munoz, Nucl. Phys. B 422 , 125 (1994) Erratum: [Nucl. Phys. B 436 , 747 (1995)] [hep-ph/9308271]; T. Kobayashi, D. Suematsu, K. Yamada and Y. Yamagishi, Phys. Lett. B 348 , 402 (1995) [hep-ph/9408322]; A. Brignole, L. E. Ibanez, C. Munoz and C. Scheich, Z. Phys. C 74 , 157 (1997) [hep-ph/9508258]; L. E. Ibanez, C. Munoz and S. Rigolin, Nucl. Phys. B 553 , 43 (1999) [hep-ph/9812397].
- 5[5] K. Choi, A. Falkowski, H. P. Nilles and M. Olechowski, Nucl. Phys. B 718 , 113 (2005) [hep-th/0503216].
- 6[6] K. Choi, K. S. Jeong, T. Kobayashi and K. i. Okumura, Phys. Lett. B 633 , 355 (2006) [hep-ph/0508029]; Phys. Rev. D 75 , 095012 (2007) [hep-ph/0612258].
- 7[7] D. Baumann and L. Mc Allister, ar Xiv:1404.2601 [hep-th].
- 8[8] S. Kachru, R. Kallosh, A. D. Linde and S. P. Trivedi, Phys. Rev. D 68 , 046005 (2003) doi:10.1103/Phys Rev D.68.046005 [hep-th/0301240].
