![](https://csdnimg.cn/release/download_crawler_static/12262127/bg3.jpg)
Eur. Phys. J. C (2018) 78:413 Page 3 of 15 413
presented in the form of an algebraic algorithm [4] and later
written in compact closed form [13] as inequalities imposed
not on the parameters
a
themselves but on four eigenval-
ues
ˆ
i
of a real symmetric 4 × 4matrix
ij
which encodes
all quartic interaction terms. The form of these conditions is
very simple and basis-invariant,
ˆ
0
> 0,
ˆ
0
>
ˆ
1,2,3
, (7)
but checking them within a specific 2HDM requires first find-
ing these eigenvalues, though this step can be easily imple-
mented in numerical scans of the parameter space.
A somewhat similar systematic method of deriving the
exact BFB conditions exists for models, in which the Higgs
potential can be written in terms of independent positive-
definite field bilinears r
i
. In this case, the quartic potential can
again be rewritten as a quadratic form V =
ij
r
i
r
j
, but its
positive definiteness must be insured only in the first orthant
r
i
≥ 0. These conditions are known as copositivity (condi-
tional positivity) criteria. They were developed in [14–16]
and applied to such cases as some 2HDMs, singlet-doublet
models, models with Z
3
symmetric scalar dark matter, and
left-right symmetric models.
Beyond two Higgs doublets, in the general N -Higgs-
doublet model (NHDM), the exact BFB conditions in closed
form are still not known. Several attempts to attack the prob-
lem with the bilinear space formalism [4,17–20] did not cul-
minate in a closed set of inequalities. The technical chal-
lenge is that, with N Higgs doublets, the space of bilinears
r
a
, a = 1,...,N
2
−1, does not span the entire R
N
2
−1
space
but only a lower-dimensional algebraic manifold, which is
described with a series of polynomial constraints. Positive-
definiteness of a quadratic form on a complicated algebraic
manifold cannot be decided with linear algebra and requires
algebraic-geometric tools, that have not been found yet.
For larger gauge symmetries and for scalars in higher-
dimensional representations, it is appropriate to analyze the
scalar potential not in the scalar fields space but in the space
of gauge orbits. This approach flourished in 1980’s with the
advent of Grand Unification models, see, for example [21–
24], and a short historical overview in [16].
In specific multi-Higgs models, in which large continu-
ous or discrete symmetry groups dramatically simplify the
potential, the exact conditions can be established [25–30].
We mention, in particular, the method developed in [28,30]
to rewrite the Higgs potential as a linear combination of new
variables, the group-invariant quartic field combinations, and
to determine the exact shape of the space spanned by these
variables. This method is similar to the so-called linear pro-
graming, and it gives the BFB constraints directly from the
description of the shape of the space available.
In certain cases, when the exact necessary and sufficient
conditions are not known but a parameter scan still needs
to be performed, it may be enough to write down a set of
sufficient conditions. They may be overly restrictive, but if a
point satisfies them, the potential is guaranteed to be positive-
definite. An example of such conditions was given for a spe-
cific 3HDM in [31]. The idea is to pick up all terms with
“angular” dependence in the scalar field space and find a
lower bound for each term separately. For example, if the
potential contains a term (φ
†
1
φ
2
)(φ
†
1
φ
3
) with real coefficient
, one can place the following lower boundary on it in the
r
i
≡|φ
i
|
2
space:
(φ
†
1
φ
2
)(φ
†
1
φ
3
) ≥−||r
1
√
r
2
r
3
≥−||r
1
(r
2
+r
3
)/2. (8)
In this way, the original potential V can be limited from below
by another potential
˜
V , which is a quadratic form in terms of
r
i
and for which the copositivity criteria are applicable.
In this work we present an algorithm, which in principle
solves the problem in a generic setting. The algorithm uses
elements of the theory of resultants and of the recently devel-
oped spectral theory of tensors. However, solving the prob-
lem in principle is quite different from solving it in practice.
To our best knowledge, the approach was only briefly men-
tioned in [16] but was not developed any further nor imple-
mented in any code. We have implemented the method in a
computer-algebra code, which is available at GitHub [32],
and tested it in cases, in which analytical solutions already
exist. The complexity of the algorithm implementation grows
so fast that, with limited computer resources, we could not
apply it to cases where the results are not yet known.
This does not imply, of course, that this direction is a
dead-end. The method itself is innovative but the specific
algorithm we propose is clearly not optimal. We believe that
with additional efforts, it can be seriously improved and may
eventually produce a ready-to-use solution in various popular
classes of multi-scalar models, such as the general 3HDM.
The structure of the paper is the following. In the next
Section, we present our strategy and formulate the algorithm.
Section 3 contains an introduction to the spectral theory of
tensors, its application to the BFB problem, and describes a
practical algorithm to calculate the characteristic polynomial
of a symmetric tensor. In Sect. 4, we show how this method
works. We first do it with two elementary examples, in which
all calculations can be performed manually, and then apply
the computer-algebra package to the case of a Z
2
-symmetric
2HDM, where the BFB conditions are known. We find agree-
ment of the results, which serves as a check of the validity
of our algorithm. We end with a discussion in Sect. 5 of how
the algorithm can be improved in the future and draw con-
clusions. The Appendix contains a pedagogical introduction
to polynomial rings and polynomial division with an appli-
cation to the theory of resultants.
123