The standard model on a domainwall brane?
Abstract
We propose a dimensional action that is a candidate for realising a standardmodellike effective theory for fields dynamically localised to a domainwall brane. Our construction is in part based on the conjecture that the DvaliShifman mechanism for dynamically localising gauge bosons to a domain wall works correctly in d. Assuming this to be so, we require the gauge symmetry to be SU(5) in the bulk, spontaneously breaking to SU(3)SU(2)U(1) inside the domain wall, thus dynamically localising the standard model gauge bosons provided that the SU(5) theory in the bulk exhibits confinement. The wall is created jointly by a real singletHiggs field configured as a kink, and an SU(5) adjointHiggs field that takes nonzero values inside the wall. Chiral dimensional quarks and leptons are confined and split along the bulk direction via their Yukawa couplings to and . The Higgs doublet and its colour triplet SU(5) partner are similarly localised and split. The splittings can suppress colouredHiggsinduced proton decay and, because of the different localisation profiles, the usual SU(5) mass relation does not arise. Localised gravity is generated via the RandallSundrum alternative to compactification.
pacs:
04.50.+h, 11.27.+d, 12.10.gI Introduction
There is no known fundamental principle requiring spacetime to be dimensional, so extra dimensions of space might exist. If so, then the effective dimensionality we observe in everyday life and in highenergy experiments has to be explained. It could be that the extra dimensions are topologically compact and small, as per the KaluzaKlein idea. Alternatively, the extra dimensions could be large but as yet unobserved because standard model (SM) fields are confined to a d brane. Large extra dimensions might be compact, as proposed by ArkaniHamed, Dimopoulos and Dvali ArkaniHamed et al. (1998), or infinite, as shown in the second of the RandallSundrum papers of 1999 Randall and Sundrum (1999a) (hereinafter RS2). See also Refs. Randall and Sundrum (1999b); Antoniadis (1990); Antoniadis et al. (1998); Akama (1983); Visser (1985).
The purpose of this paper is to propose a dimensional action that is a candidate for realising a standardmodellike effective theory for fields dynamically localised to a domainwall (DW) brane. Like RS2, there is one extra dimension and it is infinite. Unlike RS2, the brane is not a fundamental object but rather a solitonic solution of the theory, as per the Rubakov and Shaposhnikov Rubakov and Shaposhnikov (1983) proposal that we might live on a domain wall. Our construction assembles a number of dynamical localisation mechanisms into what we hope is a complete theory of a domainwalllocalised SM. These mechanisms are:

the localisation of d chiral fermion zero modes through the Yukawa coupling of d fermions to the background scalar fields;

the localisation of a SM Higgs doublet to the DW through its Higgs potential couplings to the DWforming scalar fields;

the localisation of SM gauge bosons via the DvaliShifman mechanism, instituted through a bulk that respects SU(5) gauge invariance Georgi and Glashow (1974);

the DW generalisation of the RS2 mechanism for localising gravitons.
Three of these four mechanisms involve wellestablished phenomena. The DvaliShifman (DS) gauge boson localisation idea remains an interesting conjecture in the d context, not as yet proven to work. What we shall show in this paper is that if one takes the DS mechanism to work in d, then the construction of a DWlocalised standard model follows readily, and even elegantly. We hope that our model spurs rigorous studies of the DS mechanism in d, to either confirm it or disprove it. Were it to be confirmed, then our modelbuilding setup would provide a clear pathway to the construction of phenomenologicallyrealistic effective theories of DWlocalised fields. We shall review the DS mechanism below.
The main aesthetic motivation for our model is to treat all spatial dimensions on an equal footing in the action. In particular, all these dimensions are infinite, as in the RS2 setup. But “dimensional democracy” is taken further than in RS2, because that theory has translational invariance along the extra dimension explicitly broken through the introduction of an infinitelythin fundamental brane into the action. To achieve dimensional democracy we must have no explicit brane terms in the action, but replace the RS2 fundamental brane with a finitethickness stable domainwall configuration of scalar fields.
We shall argue that our theory is likely to be the minimal way to get a purely fieldtheoretic realisation of a DWconfined SM. It is interesting that in order to achieve this the DvaliShifman mechanism immediately motivates an extension to SU(5). We are also encouraged by the fact that some of the usual problems of SU(5) grand unification have solutions automatically provided by the minimal theory, without “epicyclic” adhoc fixes. As we shall explain below, the usual style SU(5) mass relations are simply absent, because the fermion localisation realises a modified version of the splitfermion idea of ArkaniHamed and Schmaltz ArkaniHamed and Schmaltz (2000) (see also Coulthurst et al. (2006)). The downtype quarks necessarily have different bulk profile functions from the charged leptons, and because the d masses are computed from overlap integrals of profile functions, the quarklepton mass degeneracy just does not arise. The fermion splitting can also suppress colouredHiggsinduced proton decay. An important loose end is that we are not yet able to analyse gauge coupling constant unification in our unusual version of SU(5). We shall explain below why a full unification study is premature.
Our focus in this paper is on modelbuilding rather than detailed phenomenology. We wish to explain the logic of our construction, and provide evidence that it has good phenomenology without supplying absolute proof.
We review the DvaliShifman mechanism in the next section, describe our model in the following section, and conclude in the last section.
Ii DvaliShifman mechanism
The most plausible mechanism for localising gauge bosons to a DW in such a way as to preserve gauge invariance is that proposed by Dvali and Shifman Dvali and Shifman (1997). This requires a confining nonAbelian gauge theory in the bulk, with the symmetry G broken to a subgroup group H inside the DW. Massless gauge bosons corresponding to H are then localised to the wall. As we wish to localise the standard model fields, the minimal choice is to take G=SU(5) and H=SU(3)SU(2)U(1).
The truth of the DS mechanism rests on quite a firm foundation for DWs residing in a background d spacetime Rubakov (2001); Dubovsky and Rubakov (2001); ArkaniHamed and Schmaltz (1999). Following DS, let us consider the simple toy example of G=SU(2) and H=U(1). Place a U(1) source charge inside the wall. Because the SU(2)respecting bulk is in confinement phase, the electric field lines of the source charge cannot penetrate into the bulk. Instead, the field lines are repelled from the DWbulk interface thus reducing the effective dimensionality of the Coulomb field by one. Adopting the ’t HooftMandelstam proposal that confinement arises from the magnetic dual of superconductivity, the repulsion of field lines from the interface is readily understood from the dual Meissner effect Dubovsky and Rubakov (2001); ArkaniHamed and Schmaltz (1999).
Now place the source charge in the bulk. By confinement, which is tantamount to the expulsion of electric fields, the electric flux from the source must form a flux tube that ends on the domain wall Dubovsky and Rubakov (2001); ArkaniHamed and Schmaltz (1999). Once inside the wall the field lines are able to spread out in the plane of the wall. It is as if the charge was actually inside the wall: the electric field configuration is the same at large distances inside the wall irrespective of the position of the source. In the quantal situation where the position of a source charge is indefinite, it follows that the long range Coulomb field is independent of how the wavefunction depends on the coordinate perpendicular to the wall (the “extra” dimension). We shall be using this result below when we assume that gaugeuniversality for H holds independently of the bulk profiles of the trapped fields.
If H is nonAbelian, then these arguments generalise to the case of chromoelectric field line expulsion from the bulk.
Another perspective on the localisation physics is provided by the mass gap Dvali and Shifman (1997). In the bulk, becaue of confinement, the gauge bosons of H cannot themselves propagate but instead form constituents of propagating G glueballs. But the glueballs of G are massive. In the G=SU(2) and H=U(1) example, the U(1) gauge boson which is both massless and free inside the wall, must somehow incorporate itself into a massive SU(2) glueball if it propagates into the bulk. But the mass gap implies an energy cost in doing so, thus any U(1) gauge boson inside the wall is dynamically constrained to remain there. If H has nonAbelian factors that are themselves in confinement phase inside the wall, then the mass gap suppression corresponds to the H glueballs inside the wall being less massive than the G glueballs in the bulk.
These arguments are rather convincing because they rest on the wellestablished confinement property for asymptoticallyfree nonAbelian gauge theories in d. In the d case, the DS mechanism is a conjecture, because d confinement (or lack thereof) is not properly understood. The main issue is that pure YangMills theory is not renormalisable in d (or larger). At the level of lattice gauge calculations, this corresponds to the lack of a physical limit when taking the lattice spacing to zero. To expand on this point, it is known that d SU(2) has a first order phase transition for finite lattice spacing Creutz (1979). We have verified this conclusion for d SU(5) and so presumably SU(5) has a confining phase for sufficiently large values of the gauge coupling constant. This analysis cannot be extended to the continuum limit, and so we must be content with d SU(5) exhibiting confinement below a relevant cutoff of the theory. Thus we consider d DS to be an effective mechanism, valid below this cutoff, which does the job of confining gauge fields to the DW. As we remark below, any field theoretic braneworld model is nonrenormalisable and hence must be defined with an ultraviolet cutoff, so in our context we do not need to take the continuum limit.
To the best of our knowledge, the DS mechanism has not been directly checked in d, which would require more than just an analysis of the phase structure of pure YangMills theory. But we are encouraged by lattice gauge calculations in d Laine et al. (2004) which do verify the mechanism. We shall assume that it works also in d, and show that realistic model building is then quite possible.
Iii The model
We now describe our model. As stated above, the DS mechanism immediately motivates that the bulk should respect at least an SU(5) gauge symmetry. By one definition of “minimal”, the bulk symmetry should be exactly SU(5) and it should also be the symmetry of the action; the model presented below has these features. (It is also interesting to consider models not adhering to these strictures. For example, Ref. Davidson et al. (2007) describes a theory where the symmetry of the action is larger than the symmetry of the bulk.)
The SU(5) d field content is:
(1) 
plus gauge fields. The field is real, is conveniently represented as a Hermitian traceless matrix, while is a fivefold column vector of complex fields. Chirality does not exist in d, so both the ’s are Dirac fields, with being a antisymmetric matrix. The SU(5) transformations are: , , and . We shall, for simplicity, consider only one quarklepton family here, though the generalisation to three families is straightforward. The neutrino mass question is also deferred to later work.
Let us begin by ignoring gravity, to focus on the purely particlephysics aspects of the model. Later we discuss what remains the same, and what changes, when RS2style warped gravity is added. The action is
(2) 
where contains the SU(5) gaugecovariant kineticenergy terms, has the Yukawa couplings of the fermions to and ,
(3)  
and is the SU(5) Yukawa Lagrangian used to generate quark and lepton masses:
(4) 
The last term can only be written in SU(5) index notation: .
The Higgs potential is where
(5)  
(6)  
The action is invariant under the reflection discrete symmetry , , and . The d Dirac matrices are , where , and .
The theory is nonrenormalisable in d. As usual in these kinds of models, there is an implicitly assumed ultraviolet cutoff and an ultraviolet completion above that scale. We shall adopt the agnostic stance for both the existence and nature of this UV completion. Our action is perhaps best considered as the set of lowestdimensional operators, consistent with the stated symmetries, of a nonrenormalisable effective theory that is putatively to be derived from the UV completion.
The background DW is found by solving the EulerLagrange equations for an independent configuration obeying the boundary conditions
(7) 
corresponding to degenerate global minima of . The spontaneously broken reflection symmetry ensures topological stability for the DW. Numerical solutions exist for a significant region of parameter space. Purely for the sake of giving a concrete example, we can impose the parameter conditions
(8) 
with , permitting the analytic solution
(9) 
where , , and is the adjoint component associated with the weakhypercharge generator . All other components vanish. The configuration is the usual kink, while induces SU(5) SU(3)SU(2)U(1) within the DW which has width . This background solution creates the brane, and simultaneously confines SM gauge fields to it provided the DvaliShifman mechanism works with this kind of an SU(5) bulk in d.
We have checked numerically that configurations such as Eq. (9) are perturbatively stable against the formation of additional nonzero components. The other fields, , and , propagate in this background. Within the wall, the SU(5) confinement dynamics are suppressed so we can analyse classical localisation solutions for fermions and scalars in the usual way. Outside the wall, the nonperturbative SU(5) physics makes calculating impossible absent a dedicated lattice program. Since the localisation takes place with a characteristic distance scale of , ignoring the nonperturbative corrections is approximately valid.
It may be worthwhile to expand on this point. The computation of localised lowestenergy modes, such as fermion zero modes, is but the start of a systematic mode analysis whereby d fields are reinterpreted as infinite towers of d fields (generalisation of a KaluzaKlein decomposition); see, for example, Refs.George and Volkas (2007); Davies and George (2007) for an introduction to this procedure. Schematically, one employs a mode decomposition of the form () where the sum is over suitable modes and includes an integration if the modes contain a continuum. The are the d fields, and the are mode functions. The mode functions are usually chosen to obey certain suitable differential equations so that the fields are those of definite mass in the effective d theory (in the familiar KaluzaKlein case of a circular extra dimension, the mode functions are chosen to be sinusoidal for precisely this reason). However, from a mathematical point of view the set of mode functions is just some complete set of functions that permits the decomposition of without loss of generality, and so one has the usual freedom to change basis by changing the modefunction set. This is a pertinent observation for theories that employ the nonperturbative quantumfieldtheoretic DvaliShifman mechanism. In the bulk, the component fields are subject to these dynamics unless is a gauge singlet, and thus the physical meaning ascribed to the mode functions has to take this into account. There is no problem in using the same modedecomposition one would use in the absence of the nonperturbative bulk, because that is simply a mathematicallyvalid recasting of as an infinite set of components. If the bulk is indeed in confinement phase, then the gauge nonsinglet fields will not propagate as free particles, so their physical interpretation will be as constituent particles. This is conceptually no different from expressing the QCD Lagrangian in terms of quarks and gluons even though the propagating states are hadrons. Fortunately we are mainly interested in the lowest modes, whose mode functions are sharply peaked inside the domain wall, and so to a first approximation we need not be concerned with interpretive complications because of the nonperturbative bulk.
The d fermions couple to the background dependent scalar fields as per . A full mode decomposition analysis would involve writing each d fermion field as
(10) 
substituting this into d Dirac equation
(11) 
where is given by the relevant background domainwall scalar field configuration (see below), and requiring that the components satisfy the d Dirac equations,
(12) 
The mode functions then obey the Schrödingerlike equations
(13) 
with effective potentials
(14) 
The nature of the mode functions is then readily deduced from the analogy with the equivalent quantummechanical problem. In particular, note that as the potentials tend to the positive constant . We shall return to this observation later on when we consider gravity.
To analyse the localisation of the lowest mode (the zero mode) for each fermion, the full mode analysis above is unnecessary.
Instead, it suffices to
solve the Dirac equations with separated variable configurations
where the are d
zeromode leftchiral fields.
(15) 
where and is the weakhypercharge of the SM components denoted and . The SU(5) structure automatically gives different localisation points and profiles to the different SM components – splitting ArkaniHamed and Schmaltz (2000); Coulthurst et al. (2006) – depending on hypercharge and whether they are in the or the . The zeroes of
(16) 
are the localisation centres, with the bulk profiles
(17) 
To localise d leftchiral fields, all the must pass through zero with positive slope. Examples of these split profiles are given in Fig. 1.
The Higgsdoublet and colouredscalar contained in are similarly localised George and Volkas (2007) and split by their interaction with the background fields, as given by the terms in Eq. (6). Writing , where are required to satisfy a massive d KleinGordon equation with masssquared parameters , the profiles obey the Schrödingerlike equation,
(18) 
with a weakhyperchargedependent effective potential,
(19) 
The full spectrum of localised and delocalised modes is obtained by solving these eigenvalue equations, but we are interested here in only the lowest mass eigenstates. There is sufficient parameter freedom to allow while , thus setting the stage for an effective Mexicanhat potential for and hence electroweak symmetry breakdown inside the wall. An example of the effective potentials are given in Fig. 2. (The scalar spectrum also contains the kink translational zero mode; Ref. George and Volkas (2007) explains how this mode can be frozen out.)
We can now see how natural resolutions arise to some of the usual problems with an SU(5) GUT. The mass relation is not obtained, because the d Yukawa couplings depend on overlap integrals in the extra dimension, which will be different because of the different fermion localisation profiles. The coloured scalar induces proton decay through the Yukawa terms and , but this effect can be suppressed by making the relevant profile overlaps very small. For example, splitting and so that that they overlap exponentially little would suffice Coulthurst et al. (2006).
For the onefamily standard model, it is obvious that we have enough parameters to fit the quark and lepton masses. For the threefamily case, it is plausible that sufficient parameter freedom exists, though this has not been rigorously proven as yet. It is a complicated problem, because the physical observables depend on profile functions which depend in complicated ways on the Lagrangian parameters (and corrections to the classical calculations due to the effect of the nonperturbative bulk will also exist at some level).
Gauge coupling constant evolution cannot be examined until a proper phenomenological parameter fitting is done, because the higher mass modes both depend on these parameters and affect the coupling constant evolution. Since the higher mass modes are split multiplets, the running will be different from standard d nonsupersymmetric , and successful unification may be possible. Note that coupling constants run logarithmically, not through a powerlaw, in the effective d theory of localised fields.
We now turn on gravity, with Eq. (2) modified to
(20) 
where is the determinant of the metric, the 5D gravitational mass scale, the scalar curvature, and the bulk cosmological constant. The other terms now include minimal coupling to gravity. We first seek a background metric configuration that will simultaneously localise gauge bosons and gravitons. For a significant parameterspace region, the EinsteinKleinGordon equations admit numerical solutions where is a kink, is an even function that asymptotes to zero at , and the metric assumes the Minkowskibrane warped form,
(21) 
with asymptotically. The usual RandallSundrum finetuning condition involving the bulk cosmological constant must be imposed to ensure a Minkowski brane. For the special parameter choices,
(22) 
an analytic solution exists; this is useful because it serves as a concrete example:
(23)  
(24)  
(25) 
where . As in RS2, the linearised graviton fluctuation equation has a confined zeromode that is identified as the usual graviton Csaki et al. (2000).
The d fermion spectrum still contains a localised zero mode for each species. However, far from the brane the effective potentials that replace those in Eq. 14 are now driven asymptotically to zero by the exponentially decreasing warp factor Dubovsky et al. (2000), whereas in the gravityfree case they tended to the strictly positive constants . The gravitycase effective potentials are thus volcanolike and consequently support modes of arbitrarily small energy: continua starting at zero mass. This feature is quite analogous to the wellknown gravitonmode situation in the RS2 model: there is no mass gap, but rather a continuum a modes starting immediately above the localised massless graviton mode Randall and Sundrum (1999a). The spectrum similarly contains a localised standardmodel Higgs doublet plus a continuum starting at zero mass Davies and George (2007). We need to make sure that the absence of a mass gap does not spoil the existence of a lowenergy effective theory displaying dimensional reduction down to d.
Let us first for simplicity ignore the DvaliShifmanlike bulk physics, or focus, if you like, on the modes of the gaugesinglet scalar field . Except for discrete resonant modes, corresponding to quasilocalised states Dubovsky et al. (2000); Davies and George (2007), the lowestmass continuum modes are suppressed on the brane because they have to tunnel through the potential barrier of the volcanolike effective potential. As such, their integrated effects at low energies will be dominated by the zero modes, just as in the wellknown graviton case Randall and Sundrum (1999a). Because of this phenomenon, the localised modes do indeed form a low energy effective d theory. A detailed discussion of these matters can be found in Ref. Davies and George (2007).
The analysis of the gauge nonsinglet fermion and scalar modes is affected by the DS phenomenon. Since the continuum modes penetrate into the bulk, they feel the full effects of the confinementphase physics we assume holds there. We therefore expect the lowmass continuum modes to manifest physically as the constituents of massive “hadrons” in the bulk. But the lowest mass hadrons still have to tunnel through the volcanolike potential barriers to get inside the domain wall, and since the nonperturbative effects switch off near the wall, the situation analysed in the previous paragraph is regained and with some plausibility the same conclusions follow.
Having described the construction of the model, it is now worth surveying the various scales it contains and how they should relate to each other. Of the many scales in the model, four need careful consideration: the ultraviolet cutoff , the SU(5) breaking scale on the brane , the bulk SU(5) confinement scale and the DW inverse width . All of these scales must be well above the electroweak scale. Within the four, the required hierarchy is
(26) 
For obvious reasons, the UV cutoff must be the highest scale in the theory. The SU(5) breaking scale on the brane must be higher than the SU(5) bulk confinement scale, because we need to suppress the SU(5) confinement dynamics on the brane. If the opposite were the case, then the dynamics of the field would be everywhere dominated by the strong SU(5) interactions and our classical background scalar field configuration would have no physical relevance. Finally, the SU(5) bulk glueball radius scale must be smaller than the width of the DW in order for the DvaliShifman effect to work, as discussed in the lattice gauge analysis of Ref.Laine et al. (2004). This translates into the confinement scale being higher than the inverse wall width. The UV, DWwidth and SU(5) breaking scales are governed by free parameters, so the required hierarchy amongst those three can always be achieved. The SU(5) confinement scale is in principle to be calculated from the UVcutoff bulk SU(5) gauge theory, and will depend on and the dimensionful gauge coupling constant . If the qualitative behaviour of the pure YangMills theory discussed in Sec. II also holds for the complete theory, then we expect there to be a critical coupling above which the theory is confining. The hypothetical lattice gauge theory calculation would have to allow values of to furnish a that obeyed Eq. (26). This calculation has not been performed.
Iv Conclusion
In summary, we have proposed a candidate d action for realising a SMlike theory plus gravity dynamically localised to a domain wall. The dynamical localisation mechanisms for fermions, scalars and gravitons are well understood, whereas gauge boson localisation is postulated by way of the DvaliShifman mechanism. The DS mechanism is at this stage a conjecture in the d context because of an incomplete understanding of confinement. What we have shown is that it is quite straightforward to construct a DWlocalised SM if confinement exists for an SU(5) gauge theory bulk.
The proposed model – a d SU(5) gauge theory minimally coupled to gravity – enjoys some interesting qualitative features. Notably, the usual treelevel SU(5) relation is automatically absent and colouredHiggsinduced proton decay can be suppressed.
There are a number of open problems, including the following:

the veracity of the DvaliShifman mechanism in d, as discussed above;

to understand the phenomenological implications, including for proton decay, of the gauge bosons that are massive inside the domain wall;

to see whether there is enough parameter freedom to fit the threefamily standardmodel masses and mixing angles while obeying experimental bounds on proton decay;

to study how the effective d SU(3)SU(2)U(1) gauge coupling constants unify into a d SU(5) gauge coupling constant;

to generate nonzero neutrino masses;

to understand the phenomenology of the kink translational zero mode in the gravity case Shaposhnikov et al. (2005).
The DvaliShifman gauge boson localisation mechanism appears to be a keystone. If it can work in d, then a whole world of domainwall brane model building is opened up, of which the theory presented above is but an example. If it does not work, then it is not at all clear that realistic fieldtheoretic domainwall brane models exist when the extra dimension is noncompact. We hope that our efforts lead to renewed interest in the issue of confinement in higherdimensional gauge theories.
Acknowledgements.
We thank A. Kobakhidze and M. Trodden for useful discussions. DPG would like to thank Simon Catterall for discussions about confinement in dimensions and for making available some relevant lattice gauge theory software. DPG would also like to thank Martin Schmaltz for a useful discussion about the scales in the model. This work was supported by the Australian Research Council, the Puzey Bequest and the Henry & Louisa Williams Bequest to the University of Melbourne.Footnotes
 As pointed out by Dvali and Shifman Dvali and Shifman (1997), as well as localising gauge bosons the confining bulk can localise gauge nonsinglet fermions and scalar fields. However, for our application, we have to retain the seemingly redundant localisationtoakink mechanism. The DS mechanism on its own will not suffice, because it will localise vectorlike fermions not massless chiral fermions. The kink configuration is necessary for the spontaneous generation of chirality in the d effective theory.
References
 N. ArkaniHamed, S. Dimopoulos, and G. R. Dvali, Phys. Lett. B429, 263 (1998).
 L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 (1999a).
 L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999b).
 I. Antoniadis, Phys. Lett. B246, 377 (1990).
 I. Antoniadis, N. ArkaniHamed, S. Dimopoulos, and G. Dvali, Phys. Lett. B436, 257 (1998).
 K. Akama, Lecture Notes in Physics, Berlin SpringerVerlag 176, 267 (1983).
 M. Visser, Phys. Lett. B159, 22 (1985).
 V. A. Rubakov and M. E. Shaposhnikov, Phys. Lett. B125, 136 (1983).
 H. Georgi and S. L. Glashow, Phys. Rev. Lett. 32, 438 (1974).
 N. ArkaniHamed and M. Schmaltz, Phys. Rev. D61, 033005 (2000).
 A. Coulthurst, K. L. McDonald, and B. H. J. McKellar, Phys. Rev. D74, 127701 (2006).
 G. R. Dvali and M. A. Shifman, Phys. Lett. B396, 64 (1997).
 V. A. Rubakov, Phys. Usp. 44, 871 (2001).
 S. L. Dubovsky and V. A. Rubakov, Int. J. Mod. Phys. A16, 4331 (2001).
 N. ArkaniHamed and M. Schmaltz, Phys. Lett. B450, 92 (1999).
 M. Creutz, Phys. Rev. Lett. 43, 553 (1979).
 M. Laine, H. B. Meyer, K. Rummukainen, and M. Shaposhnikov, JHEP 04, 027 (2004).
 A. Davidson, D. P. George, A. Kobakhidze, R. R. Volkas, and K. C. Wali (2007), eprint arXiv:0710.3432[hepph].
 D. P. George and R. R. Volkas, Phys. Rev. D75, 105007 (2007).
 R. Davies and D. P. George, Phys. Rev. D76, 104010 (2007).
 C. Csaki, J. Erlich, T. J. Hollowood, and Y. Shirman, Nucl. Phys. B581, 309 (2000).
 S. L. Dubovsky, V. A. Rubakov, and P. G. Tinyakov, Phys. Rev. D62, 105011 (2000).
 M. Shaposhnikov, P. Tinyakov, and K. Zuleta, JHEP 0509, 062 (2005).