Exploiting Symmetries in SDP-Relaxations for Polynomial Optimization 2: A Symmetric Quartic Is Nonnegative iff Each Restriction f^ω to a 2-Partition ω Is a Sum of SquaresResearch Paper
Motivation
Deciding whether a real polynomial is non-negative on is NP-hard already for quartics, and the standard computational substitute is to ask whether it is a sum of squares (SOS) of polynomials, which is a semidefinite feasibility problem. The two notions agree only in three cases, determined by Hilbert in 1888: univariate polynomials, quadratic forms, and ternary quartic forms (equivalently, polynomials of degree at most 4 in two variables). Outside those cases there are non-negative polynomials that are not SOS, among them symmetric ones: Choi and Lam exhibited a symmetric quartic form in four variables that is non-negative but not SOS.
Many polynomial optimization problems are symmetric, invariant under permutations of the variables. Riener, Theobald, Jansson Andrén and Lasserre (arXiv:1103.0486, Math. Oper. Res. 2013) study how symmetry reduces the size of SDP relaxations. Their §5 uses the degree principle, introduced by Timofte (2003) and refined by Riener (2012), to replace one symmetric problem in variables by finitely many problems in few variables. For symmetric quartics this yields an exact SOS certificate of non-negativity in every dimension: Theorem 5.5, the goal of this mission.
Timeline:
- 1888: Hilbert shows that non-negative ternary quartic forms are sums of squares (of quadratic forms).
- 1977: Choi and Lam give a non-negative symmetric quartic form in four variables that is not SOS.
- 2003: Timofte proves the half-degree principle for symmetric polynomials (J. Math. Anal. Appl.).
- 2012: Riener refines it to semialgebraic constraints, with (J. Pure Appl. Algebra).
- 2012–2013: Riener, Theobald, Jansson Andrén and Lasserre derive Theorem 5.5.
Setting
Let . A polynomial is symmetric if for every permutation ; its degree is the total degree. A polynomial is a sum of squares if for some polynomials in the same variables.
For symmetric put
and let be the set of points with at most distinct components.
An -partition of is a vector of positive, non-increasing integers with ; denotes the set of all of them. For ,
and . For the partitions are with , ; there are of them.
Formalization targets
Goal: Theorem 5.5 (p. 23)
For and a symmetric of degree ,
with the set of 2-partitions of .
Milestones
- Proposition 5.1 (p. 21, quoted from Riener [33, Theorem 4.5]), the degree principle:
- Display (5.1) (p. 22), for :
- Hilbert's theorem (p. 23): every non-negative ternary quartic form is a sum of squares.
- The step of the proof of Theorem 5.5 (p. 23): a polynomial in two variables of degree at most 4 is non-negative if and only if it is a sum of squares.
Significance
Theorem 5.5 turns non-negativity of a symmetric quartic in any number of variables into semidefinite feasibility problems of fixed size, one per 2-partition, each in two variables. This is a certificate scheme that is exact, in contrast with the SOS relaxation in variables, which fails for the Choi–Lam quartic. Proposition 5.1 and (5.1) are the general mechanism behind it, and they apply to any symmetric problem whose degrees are small relative to : they underlie the reduced relaxation hierarchy of Theorem 5.2 of the same paper.
All four milestones and the goal are proved results in the literature; the work is formalization. As far as a search of the platform shows, neither the degree principle nor Hilbert's ternary quartic theorem has a machine-checked statement on Prove2Me, and Mathlib has the vocabulary (MvPolynomial.IsSymmetric, IsSumSq) but neither theorem. Hilbert's theorem is a classical result of real algebraic geometry that is reusable far beyond this mission, for example in any SOS-exactness argument in low dimension.
Difficulty
The obvious route to the goal is to apply Hilbert's theorem to each , and that is the whole printed proof; it gives only the direction "non-negative implies every is SOS". The converse, "every SOS implies non-negative", is not a consequence of Hilbert's theorem: it says that a symmetric quartic attains its infimum, if anywhere, on points with at most two distinct coordinate values, which is the degree principle with . The paper quotes the degree principle from Riener's earlier work without proof, and it is not available in Mathlib.
Hilbert's theorem itself is the hardest item. It is a statement about all non-negative quartic forms in three variables, with no symmetry to exploit, and every known proof needs substantial real algebraic geometry beyond what Mathlib provides; it is not a computation on a Gram matrix.
Formalization scope
Polynomials are MvPolynomial (Fin n) ℝ; evaluation is MvPolynomial.eval; degree is totalDegree; symmetric is MvPolynomial.IsSymmetric; sum of squares is Mathlib's IsSumSq in the relevant polynomial ring (the empty sum is allowed, as on the page). An -partition is encoded by its block map , required to be monotone, surjective and with non-increasing fibre sizes; this is in bijection with the page's , and is the substitution . Infima in Proposition 5.1 and (5.1) are taken in EReal, so empty feasible sets give and unbounded problems give ; no extended-real arithmetic is used.
Disclosed additions and readings, each forced by the page:
- in Theorem 5.5. For there is no 2-partition and is a counterexample.
- in (5.1). For there is no -partition and the right-hand side is .
- "Ternary quartic polynomial" is a ternary quartic form in Hilbert's theorem; the inhomogeneous reading is false.
- "Of degree 4" is "degree at most 4" in the step of the proof (milestone 4), since substitution can lower the degree; the goal keeps "degree exactly 4" as printed.
Ruled out: replacing "sum of squares" by non-negativity of (which would make the goal Proposition 5.1 alone), real-valued infima (which are on empty or unbounded sets), or arbitrary maps in place of 2-partitions.
Welcome contributions: Hilbert's theorem on ternary quartics in any standard form; homogenization and dehomogenization lemmas for MvPolynomial and IsSumSq; the degree principle, first for (the case the goal needs) and then with constraints.
Selected references
- C. Riener, T. Theobald, L. Jansson Andrén, J. B. Lasserre, Exploiting symmetries in SDP-relaxations for polynomial optimization, Math. Oper. Res. 38(1):122–141, 2013; arXiv:1103.0486v3. https://arxiv.org/abs/1103.0486, https://doi.org/10.1287/moor.1120.0558
- C. Riener, On the degree and half-degree principle for symmetric polynomials, J. Pure Appl. Algebra 216(4):850–856, 2012. https://doi.org/10.1016/j.jpaa.2011.08.012
- V. Timofte, On the positivity of symmetric polynomial functions. Part I: General results, J. Math. Anal. Appl. 284(1):174–190, 2003. https://doi.org/10.1016/S0022-247X(03)00301-9
- D. Hilbert, Über die Darstellung definiter Formen als Summe von Formenquadraten, Math. Ann. 32:342–350, 1888. https://doi.org/10.1007/BF01443605
- M. D. Choi, T. Y. Lam, Extremal positive semidefinite forms, Math. Ann. 231:1–18, 1977. https://doi.org/10.1007/BF01360024