Bled Workshops in Physics Vol. 14, No. 2 1 Mass Hierarchy and Physics Beyond the Standard Model I. Antoniadis* Department of Physics, CERN - Theory Division CH-1211 Geneva 23, Switzerland** Abstract. I discuss the status of the mass hierarchy problem and prospects for beyond the Standard Model physics in the light of the Higgs scalar discovery at the LHC and the experimental searches for new physics. In particular, I will discuss in this context low energy supersymmetry and large extra dimensions with low string scale. Povzetek. Izmerjeno Higgsovo skalarno polje, ki ga je napovedal standardni model, klice po oceni stanja teorij v fiziki osnovnih delcev. Predstavim problem hierarhije mas doslej poznanih osnovnih delcev in napovedi za meritve na LHC, ki jih ponujata supersimetricna teorija - napoveduje supersimetricne partnerje doslej poznanim delcem in poljem - in teorija strun - ki napoveduje, da bodo meritve potrdile obstoj vec kot stirih dimenzij. 1.1 Introduction During the last few decades, physics beyond the Standard Model (SM) was guided from the problem of mass hierarchy. This can be formulated as the question of why gravity appears to us so weak compared to the other three known fundamental interactions corresponding to the electromagnetic, weak and strong nuclear forces. Indeed, gravitational interactions are suppressed by a very high energy scale, the Planck mass MP - 1019 GeV, associated to a length lP - 10-35 m, where they are expected to become important. In a quantum theory, the hierarchy implies a severe fine tuning of the fundamental parameters in more than 30 decimal places in order to keep the masses of elementary particles at their observed values. The reason is that quantum radiative corrections to all masses generated by the Higgs vacuum expectation value (VEV) are proportional to the ultraviolet cutoff which in the presence of gravity is fixed by the Planck mass. As a result, all masses are "attracted" to about 1016 times heavier than their observed values. Besides compositeness, there are two main theories that have been proposed and studied extensively during the last years, corresponding to different approaches of dealing with the mass hierarchy problem. (1) Low energy supersym-metry with all superparticle masses in the TeV region. Indeed, in the limit of exact supersymmetry, quadratically divergent corrections to the Higgs self-energy * On leave from CPHT (UMR CNRS 7644) Ecole Polytechnique, F-91128 Palaiseau ** e-mail: ignatios.antoniadis@cern.ch are exactly cancelled, while in the softly broken case, they are cutoff by the su-persymmetry breaking mass splittings. (2) TeV scale strings, in which quadratic divergences are cutoff by the string scale and low energy supersymmetry is not needed. Both ideas are experimentally testable at high-energy particle colliders and in particular at LHC. On the other hand, the recent major discovery of the Higgs boson at the LHC with a mass around 126 GeV is so far compatible with the Standard Model within 2ct and its precision tests. It is also compatible with low energy supersymmetry, although with some degree of fine-tuning in its minimal version. Indeed, in the minimal supersymmetric Standard Model (MSSM), the lightest Higgs scalar mass mh satisfies the following inequality: mh < mz cos2 2| + 3 m4 (4n)2 "V i m2 A2 ln—2 + —2 22 1 m2 m2 AL 12mt2, < (130GeV)2, (1.1) where the first term in the r.h.s. corresponds to the tree-level prediction and the second term includes the one loop corrections due to the top and stop loops. Here, mZ,mt,mt are the Z-boson, the top and stop quark masses, respectively, v = = \Jv2 + v2 with vt the VEVs of the two higgses, tan | = v2/vi, and At the trilinear stop scalar coupling. Thus, a Higgs mass around 126 GeV requires a heavy stop mt — 3 TeV for vanishing At, or At — 3mt — 1.5 TeV in the 'best' case. These values are obviously consistent with the present LHC bounds on supersymmetry searches, but they will certainly be probed in the next run at double energy. Theoretically, they imply a fine-tuning of the electroweak (EW) scale at the percent to per mille level. This fine-tuning can be alleviated in supersymmetric models beyond the MSSM, as discussed in the next session. 1.2 MSSM Higgs sector with dimension-five and dimension-six operators Although extremely successful, the Standard Model or its supersymmetric version (MSSM) is not a fundamental theory, and this motivated the theoretical efforts to understand the nature of new physics beyond it. This search can be done using an effective field theory approach, in which the "new physics" is parametrised by effective operators. The power of this approach resides in arranging these operators in powers of 1/M* where M* is the scale of new physics that generated them. To improve the predictive power, one considers additional organising principles, such as: (i) symmetry constraints that these operators should respect, often inspired by phenomenology (for example: R-parity, lepton or baryon number conservation, etc). (ii) a truncation of the series of operators to a given order in the power of the inverse scale 1/M*. The effective low-energy Lagrangian then takes the form L = + L Mn (^2) i,n * where is the SM or the MSSM Lagrangian; On>i is an operator of dimension d = n + 4 with the index i running over the set of operators of a given dimension; cn,i are some coefficients of order O(1 ). This description is appropriate for scales E which satisfy E C M*. Constraints from phenomenology can then be used to set bounds on the scale of new physics M*. Regarding the origin of operators On,i, they can be generated classically or at the quantum level. At the classical level, this can happen by integration of some new massive states, via the equations of motion and one then generates an infinite series. This can happen even in 4D renormalisable theories; indeed, even though the low energy interaction looks nonrenormalisable, it may actually point to a renormalisable theory valid up to a much higher scale (a familiar example is the Fermi interaction). Such effective operators are also generated at the quantum level, for example following compactification of a higher dimensional theory, by the radiative corrections associated with momentum and winding modes of the compactification [1-5]. The effects of these operators on the low energy observables can be comparable to the radiative effects of light states in the SM/MSSM [6] and this shows the importance of their study. In the following we shall study these effects to the case of the MSSM Higgs sector with additional operators of dimensions d = 5 and d = 6 [6,7]. In particular we show that the mass of lightest SM-like Higgs can easily be increased close to the observed value by new physics in the region of few TeV. We then discuss the nature of the "new physics" behind the effective operators. 1.2.1 MSSM Higgs sector with d=5 and d=6 operators In the leading order, new physics beyond the MSSM Higgs sector can manifest itself as operators of either d = 5 [6,12-14] or d = 6 [7,15] or both. If generated by the same new physics, by comparing O (1/M*) and O(1/M^) terms one can estimate when the series expansion in 1/M* breaks down. There is only one operator in the Higgs sector of dimension d = 5: r 1 Ll = M d2e AH (S) (H2.H1)2 + h.c. = 2 Z10 (h2.hi )(h2.Fi + F2.H1 ) + Zii mo (h2.hi )2 + h.c, (1.3) where AH (S)/M* = Z10 + C11 m0 00 and Z10, C11 ~ 1/M*. It can be generated by integrating out a massive gauge singlet or SU(2) triplet. Indeed, in the MSSM with a massive gauge singlet, with an F-term of type M*I2 +1H1 .H2, when integrating out I generates L1. With the standard notation, here S = mo02 is the spurion superfield and m0 the supersymmetry breaking scale. We assume that m0 C M*, so that the effective approach is reliable. If this is not respected and the "new physics" is represented by "light" states (like the MSSM states), the 1/M* expansion is not reliable and one should work in a setup where these are not integrated out. The list of d = 6 operators is longer [16]: oi = 1 m2 . d4e zi (Hi eVl Hi)2, o2 = i M2 . d4e z2 (h2 eV2 H2)2, 03 = 1 M2 . d4e z3 (Hi eVl Hi) (h2 eV2 H2), 04 = 1 M2 . d4e Z4 (H2.Hi)(H2.Hi)t, 05 = 1 M2 . d4e Z5 (h| eVl Hi ) H2. Hi + h.c. 06 = 1 M2 . d4e Z6 (H2 eV2 H2) H2. Hi + h.c. 07 = 1 M2 . d2e Z7 Tr Wa Wa (H2 Hi) + h.c., 08 = 1 m2 . d4eZ8 (H2 Hi)2 + b.c., (1.4) where Wa = (-1/4) D"e-VDa eV is the chiral field strength of SU(2)L or U(1 )Y vector superfields Vw and Vy respectively. Also V^2 = VW(^a/2) + (t 1/2) Vy with the upper (minus) sign for Vi. Finally, the wavefunction coefficients are spurion dependent and have the structure (1/M") Zi(S, Sf) = Oi0 + Oil mo 66 + mo 66 + ai2 m" 6666, - 1/M". (1.5) Regarding the origin of these operators: O1,2,3 can be generated in the MSSM by an additional, massive U(1)' gauge boson or SU(2) triplets, when integrated out [12]. O4 can be generated by a massive gauge singlet or SU(2) triplet, while O5,6 can be generated by a combination of SU(2) doublets and massive gauge singlet. 07 is essentially a threshold correction to the gauge coupling, with a moduli field replaced by the Higgs. O8 exists only in non-susy case, but is generated when removing a d = 5 derivative operator by field redefinitions [6], so we keep it. 1.2.2 Higgs mass corrections from d = 5 and d = 6 operators. With the above set of independent, effective operators, one finds the scalar potential V and its EW minimum which is perturbed by O(1/M^) corrections from that of the MSSM. From V one computes the mass of CP-odd/even Higgs fields: mA = (mA)MSSM 2 Zio Ho V2 sin 2p + 2mo Zii v2 + 6mA, 6mA = O(1/M2) (1.6) for the pseudoscalar Higgs, with (m^)MSSM the MSSM value, with 6m^ due to O(1/M") corrections from d = 5 and d = 6 operators. For the CP-even Higgs one has [6,12,14] mh,H = (mh,H )MSSM + (2 Zio Po) v2 sin 2p ! ± mA + mz" + (—2 Zii mo) v2 ^ _ (mA - m|) cos2 2p - 1 T + 6m2,H, where óm2^ = O(1/M2) (1.7) The upper (lower) signs correspond to h (H), and w = (mA + m|)2 - 4 mA m| cos2 2p. With this result one can show that the mass mh can be increased near the observed value, also with the help of quantum corrections [6,12-14]. Regarding the O(1/M2) corrections of óm^ H, 6mA and óm^ H of (1.6), (1.7), in the general case of including all operators and their associated supersymmetry breaking, they have a complicated form. Exact expressions can be found in [7,15]. For most purposes, an expansion of these in 1/ tan p is accurate enough. At large tan p, d = 6 operators bring corrections comparable to those of d = 5 operators. The relative tan p enhancement of O(1/M2) corrections compensates for the extra suppression that these have relative to O(1/M*) operators (which involve both hi, h2 and are not enhanced in this limit). Note however that in some models only d = 6 operators may be present, depending on the details of the "new physics" generating the effective operators. Let us present the correction O(1/M2) to m^ for the case mA is kept fixed to an appropriate value. The result is, assuming mA >mZ, (otherwise 6mh and ómH are exchanged): 6mh = —2v2 «22 m2 + («30 + a4o)^2 + 2«6i mo Po — «20 m; - (2 Zio Ho)2 v4 (mA - m2)-1 + v2cot P (mA —m|)-^4mA ( (2a2i + «3i + «4i + 2«si ) mo Ho + (2a5o +«6o) H2 + «62 m0) — (2a6o —3a/o) mA m| — (2a.6o +«zo) mZ) + 8 (mA + m|) (ho mo Zio Zii) v2/(mA—mZ)2 + 0(1/tan2 p) (1.8) The mass corrections in (1.8) must be added to the rhs of eq.(1.7) to obtain the full value of m^. Together with (1.4), (1.5), these corrections identify the operators of d = 6 with the largest contributions, which is important for model building beyond the MSSM Higgs sector. These operators are 02,3,4 in the absence of supersymmetry breaking and 02,6 when this is broken. It is preferable, however, to increase mh by supersymmetric rather than supersymmetry-breaking effects of the effective operators, because the latter are less under control in the effective approach; also, one would favour a supersymmetric solution to the fine-tuning problem associated with increasing the MSSM Higgs mass. Therefore 02,3,4 are the leading operators, with the remark that O2 has a smaller effect, of order (mZ/^o)2 relative to O3,4 (for similar aj0, j = 2,3,4). At smaller tan O5,6 can also give significant contributions, while O7 has a relative suppression factor (mZ/^0)2. Note that we kept all operators Oi independent. By doing so, one can easily single out the individual contribution of each operator, which helps in model building, since not all operators are present in a specific model. One limit to consider is that where the operators of d = 6 have coefficients such that their contributions add up to maximise Sm^. Since a^ are not known, one can choose: -a22 = -«61 = -«30 = -«40 = «20 > 0 (1.9) In this case, at large tan | : Smh « 2v2a20 [m0 + 2m0^0 + 2 ^ + mZ] (1.10) A simple numerical example is illustrative. For m0 = 1 TeV, |x0 = 350 GeV, and with v « 246 GeV, one has Smh « 2.36 a20 x 1011 (GeV)2. Assuming M* = 10 TeV and ignoring d = 5 operators, with a20 - 1/M2 and the MSSM value of mh taken to be its upper classical limit mZ (reached for large tan |), we obtain an increase of mh from d = 6 operators alone of about Amh = 12.15 GeV to mh « 103 GeV. An increase of a20 by a factor of 2.5 to a20 - 2.5/M2 would give Amh « 28 GeV to mh « 119.2 GeV, which is already above the LEP bound. Note that this increase is realised even for a scale M* of new physics beyond the LHC reach. The above choice of M* = 10 TeV was partly motivated by the fine-tuning results [13] and on convergence grounds: the expansion parameter of our effective analysis is mq/M* where mq is any scale of the theory, in particular it can be m0. For a susy breaking scale m0 - O(1) TeV (say m0 = 3 TeV) and c1,2 or Oij of Zi(S,St) of order unity (say c1,2 = 2.5) one has for M* = 10 TeV that ci,2 m0/M* = 0.75 which is already close to unity, and at the limit of validity of the effective expansion in powers of 1/M*. To conclude, even for a scale of "new physics" above the LHC reach, one can still classically increase mh to near the LHC measured value. 1.3 Strings and extra dimensions The appropriate and most convenient framework for low energy supersymmetry and grand unification is the perturbative heterotic string. Indeed, in this theory, gravity and gauge interactions have the same origin, as massless modes of the closed heterotic string, and they are unified at the string scale Ms. As a result, the Planck mass is predicted to be proportional to Ms: Mp = Ms/g, (1.11) where g is the gauge coupling. In the simplest constructions all gauge couplings are the same at the string scale, given by the four-dimensional (4d) string coupling, and thus no grand unified group is needed for unification. In our conventions aGUT = g2 ~ 0.04, leading to a discrepancy between the string and grand unification scale Mgut by almost two orders of magnitude. Explaining this gap introduces in general new parameters or a new scale, and the predictive power is essentially lost. This is the main defect of this framework, which remains though an open and interesting possibility. The other other perturbative framework that has been studied extensively in the more recent years is type I string theory with D-branes. Unlike in the heterotic string, gauge and gravitational interactions have now different origin. The latter are described again by closed strings, while the former emerge as excitations of open strings with endpoints confined on D-branes [17]. This leads to a braneworld description of our universe, which should be localized on a hypersurface, i.e. a membrane extended in p spatial dimensions, called p-brane (see Fig. 1.1). Closed strings propagate in all nine dimensions of string theory: in those extended along the p-brane, called parallel, as well as in the transverse ones. On the contrary, open strings are attached on the p-brane. Obviously, our p-brane world must Fig. 1.1. D-brane world universe in type I string framework. have at least the three known dimensions of space. But it may contain more: the extra d|| = p — 3 parallel dimensions must have a finite size, in order to be unobservable at present energies, and can be as large as TeV-1 - 10-18 m [18]. On the other hand, transverse dimensions interact with us only gravitationally and experimental bounds are much weaker: their size should be less than about 0.1 mm [19]. In the following, I review the main properties and experimental signatures of low string scale models [20]. 1.3.1 Framework of low scale strings In type I theory, the different origin of gauge and gravitational interactions implies that the relation between the Planck and string scales is not linear as (1.11) of the heterotic string. The requirement that string theory should be weakly coupled, constrain the size of all parallel dimensions to be of order of the string length, while transverse dimensions remain unrestricted. Assuming an isotropic transverse space of n = 9 — p compact dimensions of common radius R^, one finds: MP = 92 M2+nRT , 9s ~ g2 . (1.12) 9s where gs is the string coupling. It follows that the type I string scale can be chosen hierarchically smaller than the Planck mass [21,20] at the expense of introducing extra large transverse dimensions felt only by gravity, while keeping the string coupling small [20]. The weakness of 4d gravity compared to gauge interactions (ratio MW/MP) is then attributed to the largeness of the transverse space R^ compared to the string length ls = M-1. An important property of these models is that gravity becomes effectively (4 + n)-dimensional with a strength comparable to those of gauge interactions at the string scale. The first relation of Eq. (1.12) can be understood as a consequence of the (4 + n)-dimensional Gauss law for gravity, with MÍ4+n) = M2+n/g4 the effective scale of gravity in 4 + n dimensions. Taking Ms ~ 1 TeV, one finds a size for the extra dimensions R^ varying from 108 km, .1 mm, down to a Fermi for n = 1,2, or 6 large dimensions, respectively. This shows that while n = 1 is excluded, n > 2 is allowed by present experimental bounds on gravitational forces [19,22]. Thus, in these models, gravity appears to us very weak at macroscopic scales because its intensity is spread in the "hidden" extra dimensions. At distances shorter than R^, it should deviate from Newton's law, which may be possible to explore in laboratory experiments. 1.3.2 Experimental implications in accelerators We now turn to the experimental predictions of TeV scale strings. Their main implications in particle accelerators are of four types, in correspondence with the four different sectors that are generally present: 1. New compactified parallel dimensions; In this case RMs > 1, and the associated compactification scale R-1 would be the first scale of new physics that should be found increasing the beam energy [18,23]. The main consequence is the existence of KK excitations for all SM particles that propagate along the extra parallel dimensions. These can be produced on-shell at LHC as new resonances [24]. 2. New extra large transverse dimensions and low scale quantum gravity,. The main experimental signal is gravitational radiation in the bulk from any physical process on the world-brane [25]. 3. Genuine string and quantum gravity effects. Direct production of string res- onances in hadron colliders leads generically to a universal deviation from Standard Model in jet distribution [26]. In particular, the first Regge excitation of the gluon has spin 2 and a width an order of magnitude lower than the string scale, leading to a characteristic peak in dijet production; similarly, the first excitations of quarks have spin 3/2. 4. Extra U(1 )'s arising generically in D-brane models as part of unitary gauge group factors. They obtain in general masses due to four- or higher-dimensional anomalies, via the so-called Green-Schwarz anomaly cancellation mechanism involving axionic fields from the closed string sector. The resulting masses are therefore suppressed by a loop factor compared to the string scale. From the low energy point of view, they gauge global symmetries of the Standard Model, such as the baryon and lepton number. An important property of the anomaly cancellation mechanism is that the anomalous U(1) gauge bosons acquire masses leaving behind the corresponding global symmetries unbroken in perturbation theory. Thus, this is a way to guarantee proton stability (from unbroken baryon number) and avoid large Majorana neutrino masses (from unbroken lepton number) due to dimension-5 operators involving two higgses and two leptons that are suppressed only by the TeV string scale. Such extra U(1 )s have interesting properties and distinct experimental signatures [27-29]. 5. Concerning possible micro-black hole production, note that a string size black hole has a horizon radius rH - 1 in string units, while the Newton's constant behaves as GN - g^ .It follows that the mass of a d-dimensional black hole is [30]: Mbh - rH/2-1 / Gn ~ 1/g2. Using the value of the SM gauge couplings gs ~ g2 - 0.1, one finds that the energy threshold MBH of micro-black hole production is about four orders of magnitude higher than the string scale, implying that one would produce 104 string states before reaching MBH. 1.3.3 Electroweak symmetry breaking Non-supersymmetric TeV strings offer also a framework to realize gauge symmetry breaking radiatively. Indeed, from the effective field theory point of view, one expects quadratically divergent one-loop contributions to the masses of scalar fields. The divergences are cut off by Ms and if the corrections are negative, they can induce electroweak symmetry breaking and explain the mild hierarchy between the weak and a string scale at a few TeV, in terms of a loop factor [31]. More precisely, in the minimal case of one Higgs doublet H, the scalar potential is: V = A(HH)2 + M2(HtH), (1.13) where A arises at tree-level. Moreover, in any model where the Higgs field comes from an open string with both ends fixed on the same brane stack, it is given by an appropriate truncation of a supersymmetric theory. On the other hand, m2 is generated at one loop: M2 = —e2 g2 M2 , (1.14) where e is a loop factor that can be estimated from a toy model computation and varies in the region e - 10-1 — 10-3. Indeed, consider for illustration a simple case where the whole one-loop effective potential of a scalar field can be computed. We assume for instance one extra dimension compactified on a circle of radius R > 1 (in string units). An interesting situation is provided by a class of models where a non-vanishing VEV for a scalar (Higgs) field ^ results in shifting the mass of each KK excitation by a constant a($): Mm=(m+Ra^) , (1.15) with m the KK integer momentum number. Such mass shifts arise for instance in the presence of a Wilson line, a = q | ^y A, where A is the internal component of a gauge field and q the charge of a given state under the corresponding generator. A straightforward computation shows that the ^-dependent part of the one-loop effective potential is given by [32]: Veff = —Tr(-)F -—2 e2nina dl l3/2fs(l) e-n n R l (1.16) where F = 0,1 for bosons and fermions, respectively. We have included a regulating function fs (l) which contains for example the effects of string oscillators. To understand its role we will consider the two limits R >> 1 and R << 1. In the first case only the l —» 0 region contributes to the integral. This means that the effective potential receives sizable contributions only from the infrared (field theory) degrees of freedom. In this limit we would have fs (l) —» 1. For example, in the string model considered in [31]: 4 fs(l) = 4$ Cl +1 ) ->1 for l->0, (1.17) and the field theory result is finite and can be explicitly computed. As a result of the Taylor expansion around a = 0, we are able to extract the one-loop contribution to the coefficient of the term of the potential quadratic in the Higgs field. It is given by a loop factor times the compactification scale [32]. One thus obtains - g2/R2 up to a proportionality constant which is calculable in the effective field theory. On the other hand, if we consider R —» 0, which by T-duality corresponds to taking the extra dimension as transverse and very large, the one-loop effective potential receives contributions from the whole tower of string oscillators as appearing in fs(l), leading to squared masses given by a loop factor times M:2, according to eq. (1.14). More precisely, from the expression (1.16), one finds: 14 i(R) = 2^2 dl 04 r-1 , M d3 V n2e-2nn2R2l (21)5/2 V2 V 2 il +1 R3X n2e-2n . (1.18) For the asymptotic value R 0 (corresponding upon T-duality to a large transverse dimension of radius 1/R), e(0) ~ 0.14, and the effective cut-off for the mass term is Ms, as can be seen from Eq. (1.14). At large R, |j.2 (R) falls off as 1/R2, which is the effective cut-off in the limit R —» oo, as we argued above, in agreement with field theory results in the presence of a compactified extra dimension.In fact, in the limit R —» oo, an analytic approximation to e(R) gives: £(R) - M^R , = " °.°°8. (L19) The potential (1.13) has the usual minimum, given by the VEV of the neutral component of the Higgs doublet v = \J— ^2/A. Furthermore, from (1.14), one can compute Ms in terms of the Higgs mass mH = —2|j.2: Ms = , yielding naturally values in the TeV range. Acknowledgments Work supported in part by the European Commission under the ERC Advanced Grant 226371. References 1. D. M. Ghilencea and H. M. Lee, JHEP 0509 (2005) 024. 2. D. M. Ghilencea, Phys. Rev. D 70 (2004) 045018. 3. S. Groot Nibbelink, M. Hillenbach, Phys. Lett. B 616 (2005) 125. 4. S. Groot Nibbelink, M. Hillenbach, Nucl.Phys. B 748 (2006) 60. 5. J. F. Oliver, J. Papavassiliou and A. Santamaria, Phys. Rev. D 67 (2003) 125004. 6. I. Antoniadis, E. Dudas, D. Ghilencea, P. Tziveloglou, Nucl. Phys. B 808 (2009) 155. 7. I. Antoniadis, E. Dudas, D. Ghilencea, P. Tziveloglou, Nucl. Phys. B 831 (2010) 133. 8. G. F. Giudice, C. Grojean, A. Pomarol and R. Rattazzi, JHEP 0706 (2007) 045. 9. S. Cassel, D. M. Ghilencea, G. G. Ross, Nucl. Phys. B 835 (2010) 110. 10. S. Cassel, D. M. Ghilencea and G. G. Ross, Phys. Lett. B 687 (2010) 214. 11. R. Barate et al. [LEP Working Group for Higgs boson searches and ALEPH, DELPHI, L3, OPAL Collaborations, Phys. Lett. B 565 (2003) 61; S. Schael et al. ALEPH, DELPHI, L3, OPAL and LEP Working Group for Higgs Boson Searches, Eur. Phys. J. C47(2006)547 12. M. Dine, N. Seiberg, S. Thomas, Phys. Rev. D 76 (2007) 095004. 13. S. Cassel, D. Ghilencea, G. G. Ross, Nucl. Phys. B 825 (2010) 203. 14. K. Blum, C. Delaunay and Y. Hochberg, arXiv:0905.1701 [hep-ph]. 15. M. Carena, K. Kong, E. Ponton and J. Zurita, arXiv:0909.5434 [hep-ph]. M. Carena, E. Ponton and J. Zurita, arXiv:1005.4887 [hep-ph]. 16. D. Piriz and J. Wudka, Phys. Rev. D 56 (1997) 4170. A. Brignole, J. A. Casas, J. R. Espinosa and I. Navarro, 17. C. Angelantonj and A. Sagnotti, Phys. Rept. 371 (2002) 1 [Erratum-ibid. 376 (2003) 339] [arXiv:hep-th/0204089]. 18. I. Antoniadis, Phys. Lett. B 246 (1990) 377. 19. D. J. Kapner, T. S. Cook, E. G. Adelberger, J. H. Gundlach, B. R. Heckel, C. D. Hoyle and H. E. Swanson, Phys. Rev. Lett. 98 (2007) 021101. 20. N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 429 (1998) 263 [arXiv:hep-ph/9803315]; I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 436 (1998) 257 [arXiv:hep-ph/9804398]. 21. J. D. Lykken, Phys. Rev. D 54 (1996) 3693 [arXiv:hep-th/9603133]. 22. J.C. Long and J.C. Price, Comptes Rendus Physique 4 (2003) 337; R.S. Decca, D. Lopez, H.B. Chan, E. Fischbach, D.E. Krause and C.R. Jamell, Phys. Rev. Lett. 94 (2005) 240401; R.S. Decca et al., arXiv:0706.3283 [hep-ph]; S.J. Smullin, A.A. Geraci, D.M. Weld, J. Chi-averini, S. Holmes and A. Kapitulnik, arXiv:hep-ph/0508204; H. Abele, S. Haeßler and A. Westphal, in 271th WE-Heraeus-Seminar, Bad Honnef (2002). 23. I. Antoniadis and K. Benakli, Phys. Lett. B 326 (1994) 69. 24. I. Antoniadis, K. Benakli and M. Quiros, Phys. Lett. B 331 (1994) 313 and Phys. Lett. B 460 (1999) 176; P. Nath, Y. Yamada and M. Yamaguchi, Phys. Lett. B 466 (1999) 100 T. G. Rizzo and J. D. Wells, Phys. Rev. D 61 (2000) 016007; T. G. Rizzo, Phys. Rev. D 61 (2000) 055005; A. De Rujula, A. Donini, M. B. Gavela and S. Rigolin, Phys. Lett. B 482 (2000) 195. 25. G.F. Giudice, R. Rattazzi and J.D. Wells, Nucl. Phys. B 544 (1999) 3; E.A. Mirabelli, M. Perelstein and M.E. Peskin, Phys. Rev. Lett. 82 (1999) 2236; T. Han, J.D. Lykken and R. Zhang, Phys. Rev. D 59 (1999) 105006; K. Cheung, W.-Y. Keung, Phys. Rev. D 60 (1999) 112003; C. Balazs et al., Phys. Rev. Lett. 83 (1999) 2112; J.L. Hewett, Phys. Rev. Lett. 82 (1999) 4765. 26. L.A. Anchordoqui, H. Goldberg, D. Lust, S. Nawata, S. Stieberger and T.R. Taylor, Phys. Rev. Lett. 101 (2008) 241803 [arXiv:0808.0497 [hep-ph]]. 27. I. Antoniadis, E. Kiritsis and T. N. Tomaras, Phys. Lett. B 486,186 (2000) [arXiv:hep-ph/0004214]; I. Antoniadis, E. Kiritsis and J. Rizos, Nucl. Phys. B 637, 92 (2002) [arXiv:hep-th/0204153]; I. Antoniadis, E. Kiritsis, J. Rizos and T. N. Tomaras, Nucl. Phys. B 660, 81 (2003) [arXiv:hep-th/0210263]. 28. G. Shiu and S.-H. H. Tye, Phys. Rev. D 58 (1998) 106007; Z. Kakushadze and S.-H. H. Tye, Nucl. Phys. B 548 (1999) 180; L. E. Ibanez, C. Munoz and S. Rigolin, Nucl. Phys. B 553 (1999) 43. 29. L. A. Anchordoqui, I. Antoniadis, H. Goldberg, X. Huang, D. Lust and T. R. Taylor, Phys. Rev. D 85 (2012) 086003 [arXiv:1107.4309 [hep-ph]]. 30. G.T. Horowitz and J. Polchinski, Phys. Rev. D 55 (1997) 6189. 31. I. Antoniadis, K. Benakli and M. Quiros, Nucl. Phys. B 583 (2000) 35. 32. I. Antoniadis, K. Benakli and M. Quiros, New Jour. Phys. 3 (2001) 20.