A new view on vacuum stability in the MSSM

A consistent theoretical description of physics at high energies requires an assessment of vacuum stability in either the Standard Model or any extension of it. Especially supersymmetric extensions allow for several vacua and the choice of the desired electroweak one gives strong constraints on the parameter space. As the general parameter space in the Minimal Supersymmetric Standard Model is huge, any severe constraint on it unrelated to direct phenomenological observations enhances the predictability of the model. We perform an updated analysis of possible charge and color breaking minima without relying on fixed directions in field space that minimize certain terms in the potential (known as"D-flat"directions). Concerning the cosmological stability of false vacua, we argue that there are always directions in configuration space which lead to very short-lived vacua and therefore such exclusions are strict. In addition to existing strong constraints on the parameter space, we find even stronger constraints extending the field space compared to previous analyses and combine those constraints with predictions for the light CP-even Higgs mass in the Minimal Supersymmetric Standard Model. Low masses for supersymmetric partners are excluded from vacuum stability in combination with the 125 GeV Higgs and the allowed parameter space opens at a few TeV.


Introduction
The Standard Model (SM) of particle physics is completed with the final discovery of the Higgs boson (the SM scalar) [1,2] which shows the expected properties in the experiment [3] and only leaves small room for deviations from the SM predictions. However, this discovery finalized a set of problems within the SM from which one is the hierarchy problem of the Higgs mass [4][5][6][7] another one the discussion about the cosmological stability of the electroweak ground state [8][9][10][11][12][13][14][15][16]. Surprisingly, the most popular extension of the SM to solve the hierarchy problem simultaneously cures the stability problem, which is the Minimal Supersymmetric Standard Model (MSSM). Besides the well-known solution of the hierarchy problem by the existence of bosonic degrees of freedom that cancel loop contributions, similar contributions render the effective potential stable-besides the property of the MSSM having an intrinsically stable Higgs potential at the tree-level. This solution to all problems, however, comes along with a bunch of new problems from which a prominent one in connection to the stability of the vacuum state is the possible destabilization of the Higgs potential by additional scalar degrees of freedom. Finally, the true vacuum of the theory is related to the absolute ground state of the scalar potential which is not exclusively dedicated to vacuum expectation values (vevs) of Higgs scalars anymore but can be due to vevs of the additional scalars that break electric and/or color charge and/or additionally baryon and lepton number. While spontaneous breaking of lepton number may be a desired solution to the origin of neutrino masses [17][18][19][20], the spontaneous breakdown of good gauge symmetries in the SM should be avoided in a way that SU(3) c × U(1) em stays intact.
A sophisticated collection of codes checking for non-standard tree-level minima, improving with the one-loop effective potential and calculating tunneling rates in presence of finite temperatures by the help of COSMOTRANSITIONS [37] is given by the VEVACIOUS collaboration [38]. Recently, the old charge and color breaking (CCB) constraints have been analyzed and tested in the light of the Higgs discovery at the LHC [39,40] with an updated tunneling analysis [41]. An investigation of the one-loop Higgs potential in the MSSM [42] reveals an interrelation of one-loop stability constraints from the Higgs sector only and tree-level CCB constraints including colored directions [43]. Considerations of vacuum stability are a widely used ingredient in studies of MSSM-like scenarios [44][45][46][47][48].
A general paradigm is that charge and color breaking minima in the MSSM most probably appear in such directions in field space where the D-terms vanish. D-terms are the quadrilinear contributions to the full scalar potential proportional to squared gauge couplings and therefore always positive and always seen as to win over any negative contribution. A first more complete and rather exhaustive analysis taking basically all directions in field space into account was given about twenty years ago by [27], where a full list of many special cases had been discussed.
Still, a complete analysis of the problem that somehow resides in a satisfactory solution is not possible. We provide a possible way to handle the existence of non-standard vacua in the MSSM scalar potential that follows the spirit of [27] and goes beyond. The minimization procedure reduces then effectively to the optimization of the necessary condition for the existence on non-standard vacua. This optimization, however, is neither unique nor unambiguously to be determined. Moreover, once the vacuum tunneling probability is addressed, a new concern for the "optimized" field direction may arise: to give the strongest bound from the vacuum metastability, configurations are rather preferred that lead to the minimal tunneling time. Whether or not this requirement can be exploited in automated computer tools may be left to the programming skills of the developers. For the pedestrian, it appears sufficient to have a clear analytical cut although those rules are indeed not sufficient but necessary. This analytical cut, however, should only distinguish between a global CCB minimum and a strictly stable "desired" electroweak vacuum.
Why is a reassessment of this problem needed? Besides the complete analysis of [27] not so much has been done on the analytical level as it is quite hard and any access lacks generality. Since this great catalog of dangerous directions and associated bounds on the parameters has been worked out, the greatest further achievement is the discovery of the Higgs boson [1,2] that appears to be very SM-like and has (for MSSM purposes) a rather high mass of m h 0 = 125 GeV as follows from the combination of ATLAS and CMS data at 7 and 8 TeV [49]. This value requires sizeable radiative corrections, that are known to be large in the MSSM [50,51]. However, the available parameter space gets very much constraint imposing the correct Higgs mass, even if one allows for a generous theoretical error of about 3 GeV in the determination of this mass [52]. Especially, to achieve this shift a large stop mixing is needed which conversely requires large trilinear soft SUSY breaking couplings [53,54], assisted maybe by a large Higgsino mass parameter. These large trilinear scalar terms, however, unambiguously lead to CCB minima and render the desired vacuum unstable. It is therefore necessary and important to put severe constraints on those terms in order to assure theoretical consistency. As long as there persists to be no discovery of any sparticles at the LHC, inferring larger lower bounds on the sparticle masses will also lead to possibly more stable configurations as larger SUSY masses themselves lead to larger shifts in the Higgs mass [52] without the need for large left-right squark mixing. Anyhow, compressed scenarios that might be hidden in the collider searches are likely to be in trouble with the stability bounds; especially if they tuned [55] in such a way to reproduce weird signatures [56,57].
We proceed in this paper as follows: after introducing the four-field scalar potential, which is basically the necessary object to deal with in connection to the influence on the Higgs mass, we derive a generic exclusion bound in Section 2. The anatomy of the CCB states described by this bound is discussed in Section 3. Finally, we conclude.

The four-field scalar potential
The MSSM in fact is a multi-scalar theory and its scalar potential is a complicated object potentially leading to undesired configurations. The configuration space depends on the vacuum expectation values of each field that are the field values at the minima of the potential. The potential in general has multiple minima where only the global one is considered to be the true ground state of the theory. If in any case the current electroweak vacuum we are believing to be sitting in is not the true one, this configuration will only be stable for a certain amount of time and due to quantum tunneling the global minimum will be reached.
Moreover, we have to take care that the potential is not unbounded from below (constraints known as UFB, i. e. unbounded from below bounds in the literature). Taking quantum corrections (and at least the one-loop effective potential) into account, those will always be rescued and the quantum potential will be bounded from below [27,42,58], whereas a new deep minimum will appear at very large field values. Contrary to large field-valued minima that usually come along with low tunneling rates into the true ground state, the minima discussed in this paper are close-by roughly with vevs around the SUSY scale (few TeV).
We are especially interested in the cross-relations of current analyses in the MSSM Higgs sector with the formation of non-standard vacua. The missing observational evidence for SUSY partners at all paired with a relatively heavy SM-like Higgs requests extreme parameter configurations. Existing analytical and semi-analytical bounds on the parameter space from the stability of the standard electroweak vacuum still are in agreement with what is needed to cope with the current situation. However, as we will see, most scenarios in the phenomenological MSSM (pMSSM) where all parameters are defined as input values at the SUSY scale suffer from charge and color breaking minima already at the SUSY scale (or slightly above). Moreover, the usual argument that tunneling rates to the deeper minimum are sufficiently small does not hold as there can be always a path in field space found where a closer vacuum shows up and fast tunneling proceeds to rolling down towards the final true vacuum. We shall explain this further.
Knowing the ground state of the theory means knowing the origin of spontaneous symmetry breaking means knowing the structure of the scalar potential. Each non-vanishing vev of fermionic or vector component fields would in addition break Lorentz symmetry and destroy the structure of space-time. Only the scalar part can break inner symmetries spontaneously and in a way which keeps external symmetries intact (not to speak about supersymmetry, but to break it we rely on soft breaking and stay ignorant about its deeper origin).
The ground state of the theory is given by the state which minimizes the potential energy density; therefore the relevant object is actually the effective potential, which at tree-level is equivalent to the classical scalar potential. In principle, quantum (one and higher loop) effects are calculable [59,60] and allow for spontaneous breaking radiatively. While the SM effective potential can be trivially made stable at the tree-level by choosing the Higgs self-coupling positive, the same coupling runs negative at higher energies and renders the electroweak vacuum metastable on cosmological scales [9,61]. In multi-scalar theories as in the MSSM, the situation is more involved already at the tree-level; a tree-level analysis of the scalar potential will result in regions of allowed parameters. Loop corrections are not expected to make unstable regions more stable around the scale of the relevant vev, although purely loop-induced minima may be missed.
The MSSM scalar potential is calculated according to some simple rules and consists of three basic contributions to which we will refer as the soft breaking, the F -term and the D-term contribution: The soft breaking part breaks supersymmetry softly and mimics the couplings of the superpotential plus additional scalar mass terms, where the F -terms basically follow from the superpotential as derivatives with respect to the scalar components where the sum over all scalar degrees of freedom is implicitly assumed to keep a plain notation. In our discussion and analysis, we consider only the chiral supermultiplets of third generation quarks as they couple with comparably large Yukawa couplings (as superpotential parameters) to the Higgs sector and also their corresponding trilinear soft SUSY breaking couplings are assumed to be large. For cleanliness and a first understanding of the "new" phenomena hidden in an old setup, we leave leptons and their superpartners out of the game as we are primarily interested in the appearance of color breaking minima. The inclusion of third generation (s)leptons is, however, trivial and follows the same procedure. We then define the (reduced) superpotential of "our" version of the MSSM by where we denote the left-handed quark doublet as Q L = (T L , B L ) and the two Higgs doublets , respectively, and the SU(2) L -invariant multiplication by the dot product. The SU(2) L singlets are put into the left-chiral supermultipletsT R = {t * R , t c R } andB R = {b * R , b c R }, respectively. Additionally, we have to break SUSY softly which is done in the usual way with scalar mass terms and trilinear couplings: The D-term part, finally, gives additional quadrilinear terms for the scalar potential associated with gauge couplings, with the corresponding hypercharges Υ φ , weak charges σ (Pauli matrices for SU(2) L -doublets φ) and color charge matrices T . Again, summation over all gauge multiplets φ is implicitly understood.
The charged Higgs directions play no role in the forthcoming discussion since on one hand the potential is SU(2)-invariant and may be always rotated into the desired shape-on the other hand, the soft SUSY breaking terms also break SU (2) in the squark sector (as top and bottom squarks are treated differently and additional left-right mixing is introduced by the A-terms). Any charge breaking Higgs vev will then be related to a color breaking squark vev anyway and we shall be able to express everything in neutral Higgs vevs, h 0 d and h 0 u (for simplicity, we drop the superscript " 0 " in the following), as well as stop and sbottom vevst andb. 1 Finally, we have the combined top/bottom-squark-Higgs scalar potential Some remarks are necessary on the structure of the scalar potential given above and how to treat the field values and their possible phases. In the previous honorable and groundbreaking works introducing charge and color breaking solutions for the first time [21,22] it is correctly stated that for potentials considered in these cases, the trilinear couplings as well as the corresponding field vevs can always be chosen real and positive. This obvious observation, however, might be used to overconstrain the field space and therefore underconstrain the constraints on the involved parameters. Indeed, the potential of Eq. (6) has some freedom in the field redefinitions; especially, it is rephasing invariant apart from the trilinear terms and the Higgs bilinear ∼ B µ h d h u . The last term is real by construction, all the others (besides the trilinears) are absolute squares of field values. Still, we do not have the freedom to rephase the fields in such a way, that the trilinear terms behave in a well defined way. In particular, the choice of all fields real and positive is not possible! We can, for sure, find a convention for the scalar quarks but not anymore for the Higgs fields. We therefore allow both h u and h d to vary in the positive and negative regime and only constrain |t| = α|φ| as well as |b| = β|φ| with a certain scalar field value φ (where we choose h u = φ).
Moreover, we set h d = ηφ with η any real number and α, β real and positive. In case, we are considering real parameters only (not the complex MSSM), the potential is symmetric in L ↔ R exchange of left-and right-handed field labels. Setting all squark fieldsq L =q R (with q = t, b) simplifies also the D-terms in the sense, that the g 2 3 contribution vanishes and the g 2 2 and g 2 1 are the same in terms of the fields. The commitment to real parameters (and fields!) nevertheless is also a severe constraint, that may be, however, compassed by imposing global CP-invariance of the theory (i. e. CP-invariance of SUSY breaking if one refers to the A-terms). It is therefore a good assumption to consider real fields only and just constrain the colored scalars to be positive (as the colored potential is invariant underq → −q). 2 Applying the considerations from above, we now have where we applied in Eq. (6) Rewriting finally the potential, we get with Each of the effective parameters in the potential Eq. (9) depends implicitly on the scaling The minimization of the one field potential is done trivially and also the requirement for the global minimum at Knowing about the dependence on the actual field direction, this bound can be improved as or rather Note, that we easily recover the famous "traditional" CCB bound by Frère et al. [21] setting η = β = 0 and α = 1 which corresponds to the ray |t L | = |t R | = |h u | in field space: Similar expressions can be easily achieved for different field directions. The specific choices have been made to make all gauge coupling contributions in Eq. (10c) vanish-though the quartics from the Yukawa couplings, which are numerically much larger, remain. There exists no real solution for η with non-vanishing α and/or β to have λ = 0. So there will be (large) quartics anyway, which on the other hand means that we do not necessarily need to restrict to the D-flat condition which explicitly forces all g 2 i -terms in the scalar potential to be absent.
Similarly, by employing other alignments, we also recover the recently proposed [43] µ y b bound and the corresponding bound from the h u -b D-flat direction with either η = 0 or h d = ± 1 + α 2 |h u |, corresponding to |h d | 2 = |h u | 2 + |b| 2 , and leading to with α > 0; a reasonable fit to the numerically derived exclusion limits can be found for α ≈ 0.8.
In the past, many attempts have been exercised to significantly improve the stability bound on the trilinear A-term according to Uneq. (12). Possible replacements range from which was given (actually for t ↔ u on the first generation A-term) by [27] and improved considering the cosmological stability of the potential through tunneling effects by [34] to 4 and recently updated by [41] in the light of the Higgs discovery as which is more in agreement (numerically) with Uneq. (15) but shall only be applied to smaller values of µ and larger pseudoscalar masses m A , whereas moderate tan β. How exactly this "small", "large" and "moderate" is defined may be left to the gusto of the user.
All in all, the bounds (12)-(17) leave an undecided feeling behind and remain open the question for a robust, roughly unique and unambiguous constraint (which we also fail to provide).
We insist on the smaller sign (<) in Uneq. (12) and later on because the smaller or equal (≤) includes a degenerate vacuum with 〈φ〉 = 0 which also leads to undesired phenomenology (where we do not want to speculate about multiple degenerate vacua as done for the SM Higgs case [62]). To be on the safe side, the < is always preferred. The optimized class of conditions given in Uneq. (11) lead in general to a more involved interplay of different field directions that cannot be displayed in such a nice expression like Uneq. (12).
The meaning of such bounds stayed controversial in the literature and history. One significant improvement has been achieved by the discussion about the stability on cosmological grounds, the question whether or not the desired vacuum has had the possibility to decay to the true vacuum within the life-time of the universe. However, any (semi-)analytical constraint suffers from a distinct choice of the field configurations as any such choice influences the tunneling rate, as well.
The main task is now to find the "optimized" directions, meaning certain combinations of η, α and β that give rise to the most severe bounds à la Uneq. (11) leading to the deepest CCB minimum (and therefore the true vacuum of the theory). Numerical minimization (and maximization) can be efficiently done with many available tools. However, as we will see, the optimized direction is not necessarily the most dangerous direction as the former one is in certain cases related to very large field vevs accompanied with a rather high barrier between the trivial (local) minimum at 〈φ〉 = 0 and the true vacuum. Those configurations are related to very large tunneling times for the vacuum-to-vacuum transition and thus considered to be less dangerous. There are nevertheless slightly tilted or shifted directions in field space where the non-standard minimum lies closer and also the barrier is more complanate and therewith easier to be reached by quantum tunneling. Once the barrier is overcome, the true vacuum can be approached directly.
Before we continue with the actual analysis of the (reduced) MSSM incarnated in the full scalar potential of Eq. (6), we make a brief but necessary digression and discuss the issue of vacuum tunneling.

Instability vs. metastability
The process of finding the global minimum of a complicated potential is hazardous, even more the interpretation of the newly found configuration. Is the standard (local) vacuum stable against quantum tunneling towards this preferred true vacuum-or may there even be a path to gently roll down into the desired state? The estimate of the tunneling rate via the so-called bounce action itself is a tricky business, however, for a wide class of potentials a very pictorial approximation can be used where only the position (i.e. vev) of the deeper minimum and the maximum in between and the height of the wall is needed. For a thick wall separating the false from the true vacuum, a very convenient approximation formula was provided by [63] which is an exact solution for  Fig. 1 Eq. (18) is very convenient to check the stability of a given configuration in the reduced one-field potential without going into the details of the non-perturbative calculation. In comparison to the life-time of the universe, one finds metastable vacua for B 400, see [34].
It is not necessarily the global minimum that determines the tunneling rate to a nonstandard minimum. Numerical procedures may overlook the vacuum on one hand, but on the other hand the decay time to a local minimum may be much smaller and the transition to the deeper one does not play a role anymore. 5 We want to illustrate at two sample points with different phenomenology that both show deeper charge and color breaking minima.
As we only check for CCB minima, we do not impose this constraint in addition; a parameter To check for metastability, one may be tempted to define the field configuration and the specific ray that shows the deepest non-standard vacuum as the ideal or optimal one. However, as the new vev appears at say O(10 TeV) and the barrier in between gets sufficiently high, say O(few TeV 4 ), B is 400 in that specific direction as for the point in Fig. 2.
However, there are other directions via which the global minimum can be accessed with a much smaller tunneling time. For the sample point of Fig. 1 from above, we show a tomographic view of the scalar potential in theb-h u plane for increasing η = h d /h u in Fig. 3 and the same potential sliced differently for increasing β =b/h u in the h d -h u plane in Fig. 4. This is to illustrate that there is no unique choice for some fixed values of η and β that exclusively show a non-standard vacuum. There are wide regions in field space and all paths should be treated equal to estimate the tunneling rate. The "optimal" direction for the determination of the bounds on the potential parameters (masses, trilinear and quadrilinear couplings) should be rather given by the shortest tunneling time. As recommendation how to deal with any CCB exclusion, we declare each point that fails the condition A(η, α, β) 2 < 4λ(η, α, β)M 2 (η, α, β)  for any specific combination of η, α and β as clearly unstable. An easy (but maybe CPU intensive) way to check this is to scan over a reasonable range, e. g. η ∈ [−3, 3] and α, β ∈ [0, 2]; with a binning of 0.1 this procedure should find CCB configurations (since the field space regions are quite extended, even coarser binnings should lead to a trustable result).

Anatomy of Charge and Color Breaking
The main stability condition is given by the unequation depending on the field misalignments, as well as on all relevant model parameters. We distinguish between parameters in configuration space that independently of the model parameters can lead to instable configurations, such as α, β and η-and the model parameters that change the shape of the scalar potential as whole object (such as the soft masses, the trilinear couplings A t and A b as well as the µ-parameter). Furthermore, we request the parameters of the one-loop Higgs potential (i. e. the genuine type-II 2HDM of the MSSM) to allow for spontaneous electroweak breaking with the correct vevs. The ratio of the standard v u /v d is what we call tan β, obeying v 2 d + v 2 u = v 2 = (246 GeV) 2 . Note that finally, the "true" tan β will be given by 1/η. We neither keep tan β fixed in a sense that the true vacuum has to respect this relation, nor do we infer that from an original tan β > 1 the ratio 〈h u 〉/〈h d 〉 has to have the same property. What we actually find is that for most configurations the true vacuum seems to have 〈h d 〉 > 〈h u 〉 and the CCB vev typically shows 〈q〉 0.7〈h u 〉. Unfortunately, for the bounds deduced numerically by attempting to find the global minimum of parameter points violating (19), no expression of the scaling parameters α, β and η can be found in terms of the relevant potential parameters although they crucially depend on the MSSM parameter point. The ideal solution would be an exclusion of the form (19) with η, α and β given in terms of the model parameters. Similar attempts have been achieved by [27] where one trilinear operator at a time was considered only. For four operators this appears to be impossible.
For the numerical analysis in the following, we consider a very phenomenological version of the MSSM with all SUSY breaking parameters defined at the low (SUSY) scale without referring to any high-scale unified scenario. As the developing CCB minima typically also   show up around the same low scale, we ignore any effects from the renormalization group as the corresponding logarithms are small and only have a mild impact on the shape of the potential (see e. g. [34] and the reference therein to [72]). For the qualitative discussion this point is irrelevant anyway. Quantitatively, if desired, parameters at the relevant scale can be employed as input values for the analytical bounds.
We determine the soft SUSY breaking Higgs masses m 2 H u and m 2 H d requiring electroweak symmetry breaking via the conditions with the one-loop Higgs potential V 1 [42]. The bilinear soft breaking term is related to the pseudoscalar mass at tree-level via B µ = m 2 A sin β cos β. Our free parameters are the ratio of the two Higgs vevs at tree-level, tan β = v u /v d , the soft squark masses, which we for Including bottom Yukawa effects in the analysis of CCB minima has not been done to great extend in the literature, as y b usually is neglected because of its smallness. However, for large tan β and certain other regions in parameter space this cannot be done anymore.
Especially the ∆ b resummation for the bottom quark mass effectively changes the bottom Yukawa coupling dramatically for such regions. While y b gets lowered compared to m b /v d for large tan β, it grows severely for negative µ as can be seen from the expressions and even runs into a non-perturbative region (what is a well-known behavior). The reduction at large tan β and small but positive µ keeps this window open in the following analysis. We include the dominant contributions to ∆ b from the gluino and the higgsino loop [73][74][75][76] and get the corrected bottom Yukawa coupling with Another remark is inevitable on the relevance of the parameters in the discussion. Usually, when MSSM effects on the Higgs mass are discussed, the "stop mixing parameter"  [42]. An additional minimum seems to appear at a larger field value h u which is driven by the µterm and therefore the requirement is that this non-standard (apparently charge and color conserving!) vev does not lead to a minimal value of the potential that is lower than at the electroweak vev. Actually, this behavior is an artifact of neglecting colored directions in the potential already at the tree-level leading to an imaginary part in the example of Ref. [42] that was not understood (and therefore just ignored). As this imaginary part is related to If one commits to the genuine D-flat direction only (say |h d | 2 = |h u | 2 + |b| 2 ), similarly wrong exclusions (conclusions?) can be drawn. The comparison of these two choices has been elaborated in [43] together with the corresponding analytic bound Uneqs. (13) and (14).
The combined exclusion limit interpolates between the two and is shown in the upper right plot of Fig So far, we only analyzed the very generic potential of Eq. (6), rewritten as single-field potential (9), without any reference to current phenomenology of the MSSM. It appears that the pure theoretical consideration to have a self-consistent theory (especially having no deeper minimum than the electroweak ground state) already excludes wide regions of the available parameter space. The constraints are even stronger than the well-known strong constraints of [27]. Reasons are that we do not insist on tan β > 1 for the new vacuum and  [42,43] with no stop and h d vevs; top right: including also h d = 0, where the exclusion is now interpolating between the two scenarios of [43] and is a bit stronger (as not necessarily |h d | 2 = |h u | 2 + |b| 2 is fixed). Down left: including all field directions discussed in this paper, as in the upper row we kept A b = 0; down right: now switching on A b = A t , nearly the complete area seems to be excluded. To compare with the usual (p)MSSM literature, we have to rescale the A-terms A t → A t / y t and A b → A b / y b . In this area, y b ranges from ∼ 0.12 to ∼ 0.8 and gets large in the upper left corner of the µ-tan β plane including the ∆ b resummation but less large than m b /v d (which seems to rescue this corner once A b is switched on).   anyway not yet excluded by experiment in the simplified analyses M SUSY ≥ 1500 GeV (for a small µ-term and rather large tan β), we can enter the correct range. Increasing µ shifts the allowed regime to even larger M SUSY . It is therefore with hindsight not surprising at all, that there have been no signals of SUSY found so far in combination to the measured light Higgs mass. Without this additional crucial ingredient one might get depressed seeing the parameter space being closed, especially when the trilinear soft SUSY breaking couplings are taken equally large, A b = A t , as usually done. One the other hand, this is exactly what is observed by the non-observation of light stops so far. Very light squarks (below say 1 TeV) in connection with large squark mixing are inconsistent with a stable electroweak vacuum. In that sense, SUSY awaits her discovery in the very near future.

Conclusions
We have reported on a new view of charge and color breaking minima in the Minimal Supersymmetric Standard Model and consequently derived novel bounds on the parameter space from the self-consistency of the theory. In order to avoid any configurations that lead to an unstable electroweak ground state, large portions of the available parameter space are excluded. We have argued that the exclusions cannot be treated as metastability bounds requiring only a life-time of the false vacuum of about the age of the universe and any CCB exclusion an MSSM parameters is to be seen strict. We have extended the exhaustive work of Ref. [27] mostly by relaxing the constraint on h d to be strictly smaller than h u but lacking simple analytic expressions to cover the numerical exclusions. By analyzing both Higgs and third generation squark directions simultaneously (four fields), we cannot fix the signs of the trilinear terms to make them positive. This in addition opens a new window to exclude larger parameter regions as h d = −h u is allowed and particularly enhances the effect for certain sign combinations. A generic analytic bound on the four-field level is rather impossible; the remaining freedom, however, allows not to be too restrictive and especially allow for short-lived vacua in formerly metastable parameter regions.  Fig. 6 for A t ≈ A b ≈ ±1500 GeV).