Vacuum stability of a general scalar potential of a few fields

We calculate analytical vacuum stability or bounded from below conditions for general scalar potentials of a few fields. After a brief review of copositivity, we show how to find positivity conditions for more complicated potentials. We discuss the vacuum stability conditions of the general potential of two real scalars, without and with the Higgs boson included in the potential. As further examples, we give explicit vacuum stability conditions for the two Higgs doublet model with no explicit CP breaking, and for the Z3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {Z}_{3}$$\end{document} scalar dark matter with an inert doublet and a complex singlet. We give a short overview of positivity conditions for tensors of quartic couplings via tensor eigenvalues.


Introduction
A scalar potential has to be bounded from below to make physical sense.In the Standard Model (SM), it simply means that the self-coupling of the Higgs boson has to be positive.In an extended model with more scalar fields, the potential has to be bounded from below -the vacuum has to be stable -in the limit of large field values in all possible directions of the field space.In this limit, any terms with dimensionful couplings -mass or cubic terms -can be neglected in comparison with the quartic part of the scalar potential. 1 In quantum field theories, scalar couplings change with energy due to the renormalisation group running.The vacuum stability conditions may be satisfied at some scales and not satisfied at others.Checking the vacuum stability of the tree level potential with running couplings can help to determine the scale of validity of a given model.On the other hand, in models with classical scale invariance [1], where tree-level mass terms are absent, violation of vacuum stability conditions at a finite field range can be used to produce minima and induce symmetry breaking via dimensional transmutation as in the Coleman-Weinberg mechanism [2].
The remarkable Tarski-Seidenberg theorem implies that the question of whether the vacuum is stable or not for given values of the scalar couplings is in principle always decidable.Nevertheless, the general problem of finding whether a given polynomial is non-negative is an NP-hard problem if the degree of the polynomial is at least four [3], which is the case for renormalisable scalar potentials in four dimensions.
The most general quartic potential of real scalars is, of course, where the coupling tensor λ ijkl can be always made completely symmetric under any exchange of the indices.Alas, relatively simple complete conditions of the positivity of the potential (1) can be given only in the case of two fields (Sect.3), or of three fields, if the potential is biquadratic in one of them (Sect.4).Since V (φ) is a homogenous quartic polynomial, scaling the fields by a positive real constant c gives V (c φ) = c 4 V (φ) and does not affect vacuum stability. 2Therefore, we can write the quartic potential as where φ2 = 1 and r 0. We see that if V ( φ) is negative for some φ, then the potential tends to negative infinity as r → ∞ and the vacuum is not stable.Thus, to determine whether a potential is bounded from below in the limit of large field values, we can minimise its quartic part on a unit hypersphere, enforced by a Lagrangian multiplier λ: which yields the stationary point equations 1 The requirement of strong stability means demanding that the quartic part of the potential V 4 > 0 as the fields ϕ i → ∞, whereas V 4 0 gives stability in the marginal sense and there can be flat directions (then the quadratic or mass squared part of the potential has to be positive and there cannot be any cubic terms).For simplicity, we give conditions for strong stability, which in practice means making inequalities strict. 2 We can even scale each field φ i separately by a different positive coefficient.In particular, we can scale it by φ i → φ i /λ 1 4 i , where λ i is its self-coupling, and make the coefficient of φ 4  i equal to unity for the purpose of calculating positivity conditions (this only holds at tree level or at a fixed energy scale).
Notice that λ = 4V min since if we write the constraint on φ as g(φ, c) = c 2 − φ 2 , then on one hand dV min /dc = λc, while on the other hand dV min /dc = 4c 3 V min with c = 1 in the end.
Global or gauge symmetries of the potential may help to simplify the problem.An important special case is given by quartic potentials that are biquadratic in fields and have the form The couplings λ ij can be written as a matrix on the basis of φ 2 i .Since the squares of real fields are nonnegative, the natural domain of such potentials is not R n but the non-negative orthant R n + .Positivity on R n + or copositivity (short for 'conditional positivity') was introduced in [4]. 3Copositivity has found wide use in the field of convex optimisation, and was first used to derive vacuum stability conditions in [6].The set of copositive matrices is larger than and includes the familiar set of positive semidefinite matrices. 4hile the positive definite part of the parameter space can be easily found via Sylvester's criterion [8], the criteria for copositivity are more involved, but definite analytic procedures exist to compute them.
Even if the fields are gauge multiplets, any potential can be written in terms of field magnitudes and orbit space variables [9,10,11,12,13].In many cases the potential is a monotonous function of orbit space parameters and its minimum occurs on the boundary of the orbit space.The symmetries of the potential may restrict the variables to a more complicated space such as the 'future light cone' orbit space of the two Higgs doublet model (2HDM) [14,15,16,17] or the similar orbit spaces of the 3HDM [18] and NHDM [19].The importance of taking the orbit space into account properly can be seen, for example, in the case of the type II seesaw: the vacuum stability conditions calculated in [20] were somewhat too strong, because the orbit space parameters do not vary independently [21].
In general, the conditions for a potential to be bounded from below can be expressed in many ways.It may be possible to produce conditions for vacuum stability that are analytical but of considerable length.For more complicated potentials with several fields, one has to resort to numerics, for which the methods we present can still be useful for reducing the parameter space to scan over.The purpose of the paper is to introduce into the 'toolbox' of calculating vacuum stability conditions some methods that are specific, but useable in many practical cases, and others that are more complicated but also more general.While recent mathematical literature is concerned with approximate methods of finding positivity for polynomials of many variables, particle physics models usually deal with a few scalar fields and analytical solutions may afford more insight.
On numerous occasions, the new addition to the scalar sector consists of just a couple of real scalar singlets, often in the guise of a complex singlet.In this case the vacuum stability reduces to the problem of positivity of a general quartic polynomial.For comparison, we also derive the vacuum stability conditions in another form, using the Sylvester criteria for the matrix of scalar couplings.And, of course, no low energy scalar potential is complete without the Standard Model Higgs doublet, which we learn to include as well.Similar conditions can be derived e.g. for the 2HDM, where the potential can be considered to be a quartic polynomial in magnitudes of fields, or for more complicated models, such as the Z 3 scalar dark matter [22,23] with an inert doublet and a complex singlet.To our knowledge, the results for the potentials of the two singlets (and the Higgs) and for the Z 3 scalar dark matter are new.As for the 2HDM with real couplings, our results are in shorter form than similar results in the literature [24].
In addition, we reconsider the notion of copositivity of matrices and discuss its relation to orbit space variables.In more complicated situations, tests of (co)positivity in terms of eigenvalues of the tensor of scalar couplings can help.Similarly to a positive matrix, a tensor is positive-definite if its eigenvalues associated with its real eigenvectors are positive.
In Sect. 2 we give a brief review of orbit spaces and copositivity.In Sect. 3 we give the conditions for the general potential of two real scalars to be positive.It is not too hard to include the SM Higgs doublet into the potential in Sect. 4. In Sect. 5 we derive vacuum stability conditions for the 2HDM with no explixit CP-breaking.Sect.6 provides another illustration in the vacuum stability conditions for Z 3 scalar dark matter.In Sect.7 we introduce tensor eigenvalues as a way to determine the vacuum stability conditions for a most general scalar potential.We conclude in Sect.8.A Mathematica notebook with the conditions for all examples is included with the L A T E X source of the paper.

Copositivity & Orbit Spaces
A scalar potential (5) biquadratic in fields is bounded from below if the matrix of couplings λ ij is copositive [6].Even if the fields are higher multiplets under a gauge group, any potential can be written in terms of squares of field magnitudes and a few dimensionless orbit space variables.

Copositivity
The criteria to determine whether a matrix is positive in the usual sense are well established.A symmetric matrix A is said to be positive semidefinite if the quadratic form x T Ax 0 for all vectors x in R n .A symmetric matrix A is positive definite if the inequality is strict, x T Ax > 0 for any non-zero vector x in R n .A matrix A is positive (semi)definite if and only if (i) the eigenvalues of A are positive (nonnegative), (ii) the principal minors of A are positive (non-negative), or (iii) the principal invariants of A are positive (non-negative).The principal minors of A are determinants of the principal submatrices, which are obtained by deleting k rows and columns from A in a symmetric way, i.e. both the i 1 , . . ., i k rows and the i 1 , . . ., i k columns are deleted.The largest principal submatrix of A is A itself.
On the other hand, copositive matrices are demanded to be positive not for all vectors in the R n , but only on positive vectors in R n + .A symmetric matrix A is copositive if the quadratic form x T Ax 0 for all vectors x 0 in the non-negative orthant R n + .(The notation x 0 means that x i 0 for each i = 0, . . ., n.) A symmetric matrix A is strictly copositive if the quadratic form x T Ax > 0 for all vectors x > 0 in the non-negative orthant R n + .For matrices of low order the copositivity conditions are relatively simple.A symmetric matrix A of order 2 is copositive if and only if [25] a 11 0, a 22 0, a 12 + √ a 11 a 22 0.
A symmetric matrix A of order 3 is copositive if and only if [26,27] a 11 0, The Cottle-Habetler-Lemke theorem [28] provides a practical way to find analytical copositivity conditions for matrices of low order.Let the order n − 1 submatrices of a real symmetric matrix A of order n be copositive.Then A is not copositive if and only if or A is copositive if and only if The adjugate of A is the transpose of the cofactor matrix of A: where M ij is the (i, j) minor of A, the determinant of the submatrix resulting from deleting the ith row and jth column of A.
Another general way to test copositivity is Kaplan's test [29]: A symmetric matrix A is copositive if and only if every principal submatrix of A has no eigenvector v 0 with associated eigenvalue λ 0.
We see that while the positivity of the matrix can be checked via its eigenvalues and the direction of its eigenvectors is irrelevant in this case, to check for copositivity at worst all the eigenvalues and eigenvectors of all principal submatrices of the matrix have to be calculated.
Copositivity with respect to a closed cone K or K-copositivity means If the cone K is polyhedral, it can be represented in the form where G is a real matrix whose columns {u 1 , . . ., u n } are positively linearly independent vectors in R n that map the extremal rays of R n + (basis vectors) into the extremal rays of the cone K.In such a case, the condition (11) takes the form [5] x Therefore to check whether A is copositive on K, we can check the usual copositivity of G T AG on R n + .Kaplan's test can also be generalised to copositivity on a closed cone K by requiring that every principal submatrix of A have no eigenvector v ∈ K with associated eigenvalue λ 0.

Orbit Spaces
Because the scalar potential is a homogenous polynomial of fields, the question whether the potential is bounded from below does not depend on the norm of the field for a single scalar field (e.g. the Higgs in the SM) and depends only on ratios of the norms for multiple fields.Orbit spaces for different gauge groups and potentials were actively studied in the 1980s, especially in the context of Grand Unified Theories.We give a short review with references to some works.The space of magnitudes of fields is the non-negative orthant R n + and hence the problem of positivity of the potential can be put in terms of copositivity of the matrix (or tensor) of quartic couplings on the basis of norms of fields, with a few orbit space parameters to minimise over.
Orbit spaces in the context of spontaneous symmetry breaking are described in [30,31].The method of orbit spaces to minimise potentials is detailed in [11].The case of the adjoint representation of SU (N ), in particular of the 24 of SU (5), is detailed in [10].An illustration for a potential of two scalars in different representations is given in [9] for the orbit space for the quartic potential of the 5 and the 24 representation of SU (5).In [12] the case of SO(N ) with adjoint + vector representation is analysed.For several groups, the orbit spaces for simple potentials are described in [13].
We present a short outline of the method for a potential of one scalar multiplet.In our cursory review we will heed to the presentation of Jai Sam Kim [11].For a theory with the non-Abelian gauge group G, the quartic potential for a scalar φ in an n-dimensional irreducible representation R of G can be written as where f and g specify different gauge invariant contractions of indices.The potential V (φ) is invariant under a group transformation where T (θ) is an n-dimensional matrix representing a group element.In general where X L are the generators of the group G and θ L are parameters that specify the group element.
Because in general the multiplet φ has many components, it is hard to solve the minimisation equations (4) for them.Moreover, there is a degeneracy of the components of φ that give the same minimum of V .A gauge transformation rotates the components of φ while leaving the value of the potential unchanged.
More formally, the orbit of a particular φ with constant components (such as a vacuum expectation value) is the set of states φ θ = T (θ)φ with T (θ) an element of G.It can be shown that all the states φ θ respect the same group, the little group of the orbit, as φ does.If the group is unitary, then all the states φ θ have the same norm φ † φ.The set of orbits that respect the same little group is called the stratum of the little group.Therefore we have to look for the orbit -and its little group -that minimises the potential.
Orbits of φ are specified by invariant polynomials P (φ) [32,33,34].There is a basis set of invariant polynomials I a (φ) such that every invariant polynomial P (φ) can be expressed as a polynomial in the polynomial basis: P (φ) = P [I a (φ)].Each representation R has a different number of basic invariants.An orbit can be pictured as a point in the -dimensional space of I a .
The magnitude of φ is irrelevant to minimising the potential to find the vacuum stability conditions.The strata can be specified by dimensionless ratios of invariants, e.g.
These dimensionless ratios are called orbit parameters and can be thought of as a set of angles.The potential ( 13) can then be written as where The potential V (φ) is bounded from below if Therefore, we have to minimise the potential with respect to the orbit space parameters.It is obvious that for any φ, the range of α i is bounded from below and above: α imin α i α imax .We have to calculate the orbit space -the physical region in the orbit space parameters α i .Because λ φ + λ φ α 1 + λ φ α 2 = C describes a line in the orbit space, the minimum of the potential is on the boundary of the orbit space, in particular it can be on a cusp of the boundary curve.Similar considerations apply for more than two orbit space parameters and for several scalars in different representations.For many more details, we refer the interested reader to the works cited.

Vacuum Stability Conditions from Positivity of a Quartic Polynomial
The most general scalar potential of two real scalar fields φ 1 and φ 2 is5 Here we use the same notation λ ij in a different way than in eq. ( 5) and copositivity, but the meaning should be clear from the context.Set either φ 1 or φ 2 to zero, and it follows that the self-couplings λ 40 and λ 04 have to be positive in order for V to be positive. 6When both fields are non-zero, we can divide the potential V by φ 4  2 and choose the ratio x = φ 1 /φ 2 as the new variable (equivalently, set φ 2 = 1 to dehomogenise V ), reducing the question of vacuum stability of the potential to that of the positivity of a general quartic polynomial of one variable, which obviously is positive if it has no real roots and a 4 > 0 and a 0 > 0.
The nature of the roots of P can be determined by considering its discriminant and two additional polynomials of its coefficients, The condition for P to have only complex roots is [35,36] In the marginal case D = 0 (which we can usually ignore), the conditions must hold.For V (φ 1 , φ 2 )/φ 4 2 , the discriminant D and the polynomials Q and R are given by and the vacuum stability conditions for V (φ 1 , φ 2 ) are given by For the record, the polynomial S = λ 3 13 + 4λ 40 (2λ 40 λ 31 − λ 22 λ 13 ).Of course, we could as well divide by φ 4  1 and choose φ 2 /φ 1 as the variable. 7The discriminant D is invariant under the exchange, but Q and R are not.The allowed parameter space, of course, stays the same.If either λ 31 or λ 13 is zero due to some symmetry, this freedom permits us to simplify the expressions for Q and R. For example, if λ 31 = 0, the conditions are simpler if we choose φ 1 /φ 2 as the variable, since the remaining λ 13 term is only linear in φ 1 .In this case the condition Q > 0 ∨ R > 0 takes the form λ 22 + 2 √ λ 40 λ 04 > 0. As a cross-check, we can set both λ 31 and λ 13 to zero.Then the vacuum stability conditions become which can be simplified to the usual copositivity conditions for V with φ 2 1 and φ 2 2 as the non-negative variables.An illustration of the vacuum stability conditions for some values of parameters is given in Figure 1.The parts of the parameter space that are allowed are shown in dark green.In addition, the light green area is allowed if the scalars can have only non-negative values (see Section 3.2).If both λ 31 = λ 13 = 0, the vacuum stability conditions reduce to the usual copositivity conditions (31) (dashed line in the right panel).

Bounded From Below Conditions for a Quartic Polynomial on R +
The domain of the general quartic polynomial (21) can be restricted to non-negative real numbers R + .In this case, the allowed range of parameters is somewhat larger.The positivity conditions for the polynomial (21) for x 0 are given by [37] 8 where The main difference with the case of two real scalars is that a range of positive and opposite-sign λ 31 and λ 13 is allowed, as seen in the left hand panel of figure 1.Also note that it is trivial to restrict the domain to the negative numbers instead by taking x → −x, equivalent to changing a 1 → −a 1 and a 3 → −a 3 in the above conditions.
If the non-negative variables are magnitudes of scalar fields, then the coefficients a i may depend on additional orbit space parameters, notably phases.These may allow, in effect, to always choose the λ 31 and/or λ 13 terms to be negative.For that reason, as we will see below for the 2HDM, the conditions (32) for a positive variable can often be eschewed in favour of the simpler conditions (24) for a real variable.

Vacuum Stability Conditions from Positivity with an Affine Space
We will derive another, different form of vacuum stability conditions for the potential (20).The matrix of quartic couplings of the potential (20) in the monomial basis (φ 2  1 , φ where c is an arbitrary constant due to the ambiguity (φ The matrices Λ(c) form an affine space.If for some value of c the matrix Λ is positive-definite, then the potential (20) is bounded from below.Note that since φ 1 φ 2 ∈ R, one has to demand the usual positivity, not copositivity.
The Sylvester criterion for the positivity of Λ is given by First of all, the last inequality of (36) implies that c is real and has the same sign as λ 22 .The first inequality of ( 37) is trivially satisfied.To satisfy the last two inequalities of (37) one must have |c| |c 0 |, where The extrema of the left hand side (LHS) of the inequality (38) with respect to c are Which of c ± gives the maximum of the LHS of (38)?The second derivative of the LHS of (38) with respect to c is (1 − 3 2 c)λ 3  22 .Inserting here the solutions c ± , we see that for λ 22 > 0, the maximum is given by c + and likewise for λ 22 < 0 the maximum is c − .Of course, if for λ 22 > 0, for example, c + < c 0 , one must take c = c 0 .
Therefore, the optimal value for c is The conditions ( 36), ( 37) and (38) together with (41) define the same region of the parameter space as (29).We have traded the relative complexity of the latter for the apparent simplicity of the former at the cost of introducing the optimal coefficient c as a piece-wise function of the couplings that has a discontinuity at λ 22 = 0.
In some regions of the parameter space it is easy to find the conditions.If the coefficient matrix of M 2 , given by is positive-definite, that is then M 2 > 0 for any values of the fields.If, on the other hand, then M 2 0 for any values of the fields, and V |H| 2 =|H| 2 min > 0 must hold for all values of the fields.The intermediate situation, where for some values of the fields M 2 < 0 and for others M 2 > 0, occurs if the eigenvalues of M 2 have opposite sign, that is, the determinant of M 2 is negative: Then we can always make a transformation of the singlet fields to bring M 2 into the anti-diagonal form It is evident that such an M 2 is negative in two opposite quadrants in the φ 1 φ 2 -plane.If we dehomogenise V |H| 2 =|H| 2 min by taking e.g.φ 1 = 1, we have to restrict the other field φ 2 to a half-axis to respect the constraint M 2 < 0. But we already have conditions for a quartic polynomial to be positive (or negative) on a half-axis: the conditions (32) for a quartic in a positive variable.
To begin to bring M 2 into the form (50), first we diagonalise M 2 by the orthogonal matrix where sin θ = λ H11 sgn(λ H02 − λ H20 ) For our purposes, sgn 0 = 1.The eigenvalues of M 2 are given by In the intermediate case, we have λ H− < 0 and λ H+ > 0. After the diagonalisation Next we scale M 2 by to make it proportional to the unit matrix.Finally we rotate M 2 by into the anti-diagonal form.Altogether, we transform the fields in yielding the transformed potential V |H| 2 =|H| 2 min .After the transformation, the coefficient matrix has the form hence λ H11 = sgn(λ H02 − λ H20 ).
For comparison, if we set all the terms containing odd powers of φ 1 and φ 2 to zero, λ H11 = λ 31 = λ 13 = 0, then the positivity of the potential is given by the much simpler (strict) copositivity constraints on the matrix of couplings in the (φ 1 , φ 2 , |H| 2 ) basis, which are still necessary, but not sufficient conditions for positivity in the general case.The transformed potential V |H| 2 =|H| 2 min is rather complicated, so let us look at special cases.If λ H20 = λ H02 = 0, then M 2 already has the anti-diagonal form, min and λ H11 = λ H11 .If, on the other hand, λ H11 = 0, then M 2 is diagonal and U θ is the unit matrix.Then λ H20 and λ H02 must have opposite signs for the determinant of M 2 to be negative and sgn(λ H02 −λ H20 ) = sgn λ H02 .The transformed potential (multiplied by the irrelevant positive coefficient 4λ 2 H20 λ 2 H02 to make it simpler) is where In several particle physics models, the equivalent of M 2 already is diagonal.An example is given by the classically scale invariant Z 3 symmetric dark matter model [38].

Vacuum Stability of the 2HDM with Real Couplings
The scalar potential of two Higgs doublets H 1 and H 2 in the 2HDM with no explicit CP-violation is where we have taken the potentially complex couplings λ 5 , λ 6 and λ 7 real and parameterised the field bilinears as [39] |H The orbit space parameter ρ ∈ [0, 1] as implied by the Cauchy inequality 0 While the general form of the vacuum stability conditions for the most general 2HDM potential has been given [15,16,17] in the elegant 'light cone' formalism, the conditions can be given in a simple explicit and analytical form in terms of the potential couplings only in special cases.
If λ 6 = λ 7 = 0, we recover the vacuum stability conditions [40,41,42,43,39] for the inert doublet model If λ 6 = 0 or λ 7 = 0, then the condition λ 3 + λ 4 − λ 5 + 2 √ λ 1 λ 2 > 0 is a necessary condition [24].Finding the minimum of the general potential (65) is complicated because of its non-linear dependence on orbit parameters ρ and φ [12].It is practically impossible to minimise the polynomials D, Λ 1 and Λ 2 in the conditions (32) with respect to these parameters.Instead, we minimise the potential with respect to φ, ρ, h 1 , and h 2 , with the fields lying on the circle h as in eq. ( 3).The minimisation equations are Eq. (69) reduces to sin φ = 0 which yields φ = 0 and φ = π or cos φ = ±1.(There is another solution for φ that holds in the special case λ 4 = λ 5 , but in this limit it gives the same solution for V min as we will obtain below.)The solutions for h 1 = 0, h 2 = 0 and ρ = 0 reproduce the conditions (67).While each of the two solutions gives a region of parameter space corresponding to the allowed region (32) for a positive variable, the allowed region is their intersection which is given by the allowed region (24) of parameter space for a real variable for V φ=0 or V φ=π , illustrated in figure 3.
The extremum solutions to the equations (70), (71), ( 72) and (73) are given by Note that the the solution for ρ 2 in (74) is non-negative if h 2 1 and h 2 2 are non-negative.Because the extremum may not be a minimum or may lie outside the rectangular orbit space ρ ∈ [0, 1], cos φ ∈ [−1, 1], the potential must be separately minimised on the edges and vertices of that square.The potential is required to be positive in these parameter regions in any case as a necessary condition.
The conditions for ρ = 0 (and any φ) are given by (67).The extremum solutions (75), ( 76) and (74) are already calculated for cos φ = ±1.By continuity, the minimum value of cos φ at ρ = 1 must be cos φ = ±1 as well.Thus the conditions for ρ = 1 can be found by applying the conditions (24) to V cos φ=±1,ρ=1 /h 4  2 with x = h 1 /h 2 .The discriminant D and the polynomials Q and R are the same for φ = 0 and for φ = π and given by D cos φ=±1,ρ=1 = 16[16λ where λ 345 ≡ λ 3 + λ 4 + λ 5 .Altogether, the conditions for the 2HDM potential with real couplings to be bounded from below are where the conditions for V ρ=0 > 0 are given by (67) and p =⇒ q is equivalent to ¬p ∨ q.In fact, it is enough to check that either h 2 1 or h 2 2 is within bounds, since they are related by h 2 1 + h 2 2 = 1.The conditions (81) are illustrated in Figure 4.
The approach we use gives a simpler result than directly minimising V cos φ=±1,ρ=1 with respect to h 1 and h 2 as we did for V min , especially as the conditions (24) automatically take into account the two different values cos φ = ±1.For a potential that depends on three or more moduli of fields, e.g. the 3HDM, where we would have a h 3 besides h 1 and h 2 , this is not possible and one has to minimise the potential on a hypersphere or use tensor eigenvalues (see Section 7).

Vacuum Stability for Z 3 Scalar Dark Matter
Another physical example is given by scalar dark matter stable under a Z 3 discrete group.The most general scalar quartic potential of the SM Higgs H 1 , an inert doublet H 2 and a complex singlet S which is symmetric under a Z 3 group is [22,23] V where we have used the parametrisation (66) for the doublet bilinears and S = se iφ S , and we have minimised cos(φ+2φ S +φ λ S12 ) = −1 so λ S12 = −|λ S12 | without loss of generality.We define ).The situation is similar to the case of two real scalars and the Higgs boson in Section 4. First of all, λ S > 0 and V (h 1 , h 2 ) > 0. The conditions for V (h 1 , h 2 ) > 0 are the same as in the inert doublet model with λ 5 = 0: We minimise the potential with respect to h 1 , h 2 , s and ρ with the fields lying on a sphere, enforced by a Lagrange multiplier λ.The minimisation equations are The solution with all fields non-zero is Note that h 2 1 , h 2 2 , s 2 and V min share the same denominator.The solution s = 0 will repeat (83) and (84) and the solutions h 1 = 0 or h 2 = 0 will be made redundant by eq. ( 95).
Since ρ = 0 sets the λ S12 term to zero, V ρ=0 > 0 is biquadratic in the fields and we can calculate the positivity conditions for V ρ=0 > 0 via (strict) copositivity of the matrix of couplings in the (h 2 1 , h 2 2 , s 2 ) basis, part of which repeats conditions (83).For ρ = 1, instead of direct minimisation of the potential in all variables, we find it easier to calculate conditions very similar to the case (61) of two real scalars and the Higgs boson.To reduce V ρ=1 to a polynomial of two variables, we minimise it with respect to s 2 : Again, either M 2 (h 1 , h 2 ) > 0 and the solution for s 2 min is unphysical, or else The coefficient matrix of M 2 is given by Repeating the procedure of Section 4 step by step, the conditions for V ρ=1 > 0, in addition to λ S > 0 and (83) and (84), are given by where we take λ 31 → − sgn(λ S2 − λ S1 )λ 31 and λ 13 → − sgn(λ S2 − λ S1 )λ 13 .As in Section 4, λ ij denote the coefficients of the transformed potential V s 2 =s 2 min , where M 2 is anti-diagonal. 9Note that we are justified to use the conditions (24) for the positivity of the general quartic on reals in the second line of (99), because in (97) the coefficients of the h 3 1 h 2 and h 1 h 3 2 terms are non-positive in that case.Altogether, the conditions for the potential (82) symmetric under a Z 3 to be bounded from below are where the conditions for V (h 1 , h 2 ) > 0 are given by ( 83) and (84), the conditions for V ρ=0 > 0 are given by (95), the conditions for V ρ=1 > 0 are given by (99), and the extremum solutions h 2 1 , h 2 2 , s 2 , ρ 2 and V min are given by (90), (91), (92), ( 93) and (94).The conditions (100) are illustrated in Figure 5.

Positive Tensors
The most general scalar potential of n real singlet scalar fields φ i can be written as where Λ is the tensor of scalar couplings and φ = (φ 1 , . . ., φ n ).Clearly, the tensor Λ can always made fully symmetric under permutations of the indices of its elements.A symmetric matrix is positive if its eigenvalues are greater than zero.Can one generalise matrix eigenvalues and eigenvectors to tensors in such a way that they have similar properties?Indeed, tensor eigenvalues and eigenvectors have been defined independently by Qi in [44] (which we follow in our exposition) and Lim [45].
A real mth-order n-dimensional tensor A has n m elements A i1...im , where i j = 1, . . ., n for j = 1, . . ., m.A homogenous polynomial f (x) of n variables and degree m can be written as the tensor product In analogy with non-negative matrices, a n-dimensional tensor A is called non-negative if Ax m 0 for all x ∈ R n .For a vector x ∈ R n we define (x [m] ) i = x m i .The number λ is an eigenvalue of A if it is a solution to the equation where x is an eigenvector of A. For m = 2, the tensor A is a matrix and the equation ( 103) coincides with the usual matrix eigenvalue equation.When m > 2, the eigenvalues and eigenvectors of a tensor can be complex, but for even m, there always exist real eigenvalues and eigenvectors, so-called H-eigenvalues and H-eigenvectors.
An n-dimensional symmetric tensor of order m has n(m − 1) n−1 eigenvalues, so a coupling tensor of a renormalisable potential of n real scalar fields in 4-dimensional spacetime has n 3 n−1 eigenvalues: for example, with two real fields, there are six eigenvalues.The product of all the eigenvalues of A is the hyperdeterminant det A, that is the resultant of Ax m−1 = 0.When m = 2, the hyperdeterminant reduces to the usual matrix determinant.The sum of all eigenvalues is (m − 1) n−1 times the sum of diagonal elements or trace tr A = i A iiii .
For even m, an eigenvalue equation similar to the matrix eigenvalue equation can be given, where I is the unit tensor with I i1...im = δ i1...im .The degree of the eigenvalue equation is d = n(m−1) n−1 .
Obviously only real eigenvectors have relevance to whether the tensor A is non-negative.The tensor A is non-negative if all its eigenvalues with real eigenvectors are non-negative.All the principal subtensors of A, that is the tensors obtained by setting one or more of the variables x i in f (x) to zero, must be positive as well.
Qi also defines tensor E-eigenvalues and E-eigenvectors, which are exactly the same as solutions to the equations (4) with fields constrained to a hypersphere.For even m, there always exist real E-eigenvalues and E-eigenvectors.
Copositive tensors are defined in obvious analogy to copositive matrices [46]: A real symmetric tensor A of order m and dimension n is copositive if Ax m 0 for all x ∈ R n + .It has been shown that Kaplan's test of matrix co-positivity [29] directly generalises to copositive tensors [47]: A symmetric tensor A is copositive if and only if every principal subtensor of A has no eigenvector v > 0 with associated H-eigenvalue λ < 0.

Vacuum Stability of the Potential of Two Real Scalars
As an example, we consider again the general potential of two real scalars given by (20):

Figure 2 :
Figure 2: An example of parameter space allowed by the vacuum stability conditions (61) for the potential (42) of two real scalars and the Higgs.The values of parameters are λ H = 0.125, λ 40 = 0.125, λ 04 = 0.25, λ 31 = λ 13 = λ H20 = λ H02 = 0, except for plots where these couplings vary.Both light and dark green regions are allowed for for λ H11 = 0, while only the dark green region is allowed for λ H11 = 0.5.

Figure 5 :
Figure 5: An example of parameter space allowed by the vacuum stability conditions (100) for the potential (82) of the SM Higgs, an inert doublet and a complex singlet.The values of parameters are λ H = 0.125, λ S = 0.125, λ 2 = 0.25, λ 3 = λ 4 = λ S1 = λ S2 = 0, except for plots where these couplings vary.Both light and dark green regions are allowed for for |λ S12 | = 0, while |λ S12 | = 0.5 only allows the dark green region.In the last two panels, where |λ S12 | varies, only the dark blue region is allowed for λ S2 = 0, while λ S2 = 0.5 allows both the dark and light blue regions.