Spaces of sums of powers and real rank boundaries

We investigate properties of Waring decompositions of real homogeneous forms. We study the moduli of real decompositions, so-called Space of Sums of Powers, naturally included in the Variety of Sums of Powers. Explicit results are obtained for quaternary quadrics, relating the algebraic boundary of SSP\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{SSP}$$\end{document} to various loci in the Hilbert scheme of four points in P3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {P}^3$$\end{document}. Further, we study the locus of general real forms whose real rank coincides with the complex rank. In case of quaternary cubics the boundary of this locus is a degree forty hypersurface J(σ3(v3(P3)),τ(v3(P3)))\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J(\sigma _3(v_3(\mathbb {P}^3)),\tau (v_3(\mathbb {P}^3)))$$\end{document}.

The smallest possible r = rk K ( f ) is called the Waring rank of f .
Waring rank and decompositions have attracted attention of many mathematicians. On the one hand, their study is motivated by applications e.g. in computer sciences. On the other hand these problems are related to beautiful mathematics: geometry (through secant varieties) (Zak 2005), representation theory (through homogeneous varieties) (Landsberg and Manivel 2004), algebra (through apolar ideals and resolutions) (Ranestad and Schreyer 2000), moduli spaces (Mukai 1992a) and many more (Landsberg 2012).
The main questions that motivated us are: 1. How to find a Waring decomposition for a general form f ? 2. What is the geometry (moduli) of all decompositions? 3. How does the answer to first two questions depend on K = R or K = C? 4. When a real form admits a real Waring decomposition with r = rk C ( f )?
To investigate the geometry of the Waring decompositions one defines the Variety of Sums of Powers (Ranestad and Schreyer 2000): Here, the closure is taken in the (smoothable component of the) Hilbert scheme of r points in P(V * ) and is usually considered for K = C. Seminal works of Mukai, Ranestad, Schreyer Mukai (1992a, b, 2004; Schreyer (2000, 2013) and others provide descriptions of VSP's in special cases. Following Michałek et al. (2017), we investigate a semialgebraic subset of the real locus of VSP, Space of Sums of Powers, corresponding to those decompositions in which all i are real. We note that a real point of VSP does not have to correspond to fully real decomposition. The geometry and topology of the real locus of VSP is another very interesting subject, studied e.g. in Kollár and Schreyer (2004).

Definition 1.1
The Space of Sums of Powers, SSP( f ), is the semialgebraic subset of the real locus of VSP( f ), where all linear i forms are real. We also define its algebraic boundary ∂ alg SSP( f )-the Zariski closure of difference of the Euclidean closure of SSP( f ) minus the interior of the closure.
The first part of the paper is devoted to the study of the algebraic boundary ∂ alg SSP( f ). We start by defining the anti-polar form ( f ) for a general form of even degree in Definition 2.1. This is a generalization of the dual quadric. As we show, in several cases ( f ) governs the nonreduced structure of apolar schemes-cf. Proposition 2.6. This allows us to provide an explicit description of ∂ alg SSP( f ), when f is a quaternary quadric 2.11, extending the results obtained for ternary quadrics [Michałek et al. 2017, Section 2]. In the ternary case, due to results on Gorenstein resolutions, one could provide a description of ∂ alg SSP( f ) in terms of hyperdeterminants. We did not find a straightforward generalization to quaternary case, however our methods may be applied to ternary case. They provide a simpler, more explicit (but less intrinsic) description 2.1.1. Further, we study the geometry of ∂ alg SSP( f ) in detail, over Q and over C, relating it to the geometry of the Hilbert scheme. In particular, we discuss how the VSP intersects various loci of the Hilbert scheme.
We note that our approach towards ∂ alg SSP( f ) allows a description of intersection with subvarieties of VSP( f ). This is useful, when the VSP( f ) itself is a large variety. We apply it for quinary quadrics 2.1.3.
In the second part of the article we study the locus of real general forms for which the real rank equals the complex rank. In general, the complex rank is known by Alexander-Hrischowitz theorem (Alexander and Hirschowitz 1995). It follows that there exists a Zariski dense semialgebraic set R n,d of real forms with such real rank. The real rank boundary ∂ alg R n,d is defined as the Zariski closure of the topological boundary of R n,d .
We present an implemented, fast, deterministic algorithm, that for a general quaternary cubic f returns its unique Waring decomposition. This is an interesting example, one of the only two identifiable that is not binary (Galuppi and Mella 2016), studied both classically (Clebsch 1861) and in modern context (Oeding and Ottaviani 2013). Our algorithm can be used not only in applications, but working parametrically over the field K (t) it allows to provide a description of ∂ alg R 4,3 . We obtain Proposition 3.4; The join variety J (σ 3 (v 3 (P 3 )), τ (v 3 (P 3 ))) of the third secant variety of the third Veronese of P 3 and the tangential variety is an irreducible hypersurface of degree 40 in P 19 . It equals ∂ alg (R 4,3 ).
The geometry of ∂ alg R n,d is also related to the classical work of Hilbert on cones of sums of squares and nonnegative polynomials (Hilbert 1888). This allows us to provide a description of one component of ∂ alg R 4,4 as the dual of a variety of quartic symmetroids in Proposition 3.5.
In the Appendix 4 we present several algorithms used to: • prove that a given variety is irreducible (or reducible) over C, • compute the real locus of a variety, • describe a shape of the resolution of a generic member of a family of ideals.

Boundaries of spaces of sums of powers
Definition 2.1 (The anti-polar ( f )) Consider a homogeneous polynomial f ∈ S 2d (V * ) of degree 2d. It induces, through the middle catalecticant, a linear map Remark 2.2 Explicitly for a fixed base on V , to evaluate ( f ) on x one multiplies the inverse (or adjoint up to scalar) of the middle catalecticant from left and from right by the vector that evaluates all monomials of degree d on x.
The following proposition explains the name of ( f ), relating it to the anti-polar quartic defined in [Michałek et al. 2017, Theorem 4.1].

Proposition 2.3
For a homogeneous polynomial f ∈ S 2d (V * ) let us denote by C( f ) its middle catalecticant. If C( f ) is invertible, then (up to scalar) we have: Proof As the Definition 2.1 and the one in the proposition are intrinsic, we may fix a basis and assume l is a basis element. Then C(l 2d ) is simply given by a matrix with one nonzero entry on the diagonal. Hence, det C( f + l 2d ) − det C( f ) equals the complimentary minor to the nonzero entry. This exactly agrees with Definition 2.1 by Remark 2.2.

Example 2.4
If f is a quadric of full rank, than ( f ) is simply the dual quadric.
Remark 2.5 It is important to note that neither the inverse nor the adjoint of A f is the middle catalecticant of ( f )-cf. Dolgachev (2004) and [Dolgachev 2012, Remark 1.4.1], contrary to the case of quadrics. In fact, it is often the case that A −1 f is not a middle catalecticant of any equation of degree 2d-see e.g. [Michałek et al. 2016, Proposition 7.1].
The following proposition is based on results from Ranestad and Schreyer (2013).
Proposition 2.6 Let f ∈ S 2d (V * ) be such that the middle catalecticant A f is nondegenerate. Suppose S ⊂ P(V ) of length equal to the rank of A f is apolar to f . Then S has a nonreduced structure at a point l ∈ S if and only if ( f )(l) = 0.
Proof Consider a new schemeS defined by: where l ⊥ is the maximal ideal defining l. As A f is nondegenerate we know that (I S ) d = (0). On the other hand, by degree count we know that there exists g ∈ (IS) d . We may determine the g as follows: As A f is nondegenerate and l ⊥ d is a hypersurface in S d V * this determines g uniquely up to scalar and the above implication is an equivalence. We obtain that S is reduced at l if and only if where v d is the d-th Veronese (evaluating g as a polynomial on l is the same as evaluating g as a linear map on v d (l)). Hence, S is nonreduced at l if and only if < v d (l), A −1 f (v d (l)) >= 0, where < ·, · > is the usual dual pairing.
The following definitions will be useful in the study of SSP.
Definition 2.7 (V N S P, F) For a given form f ∈ S d (V * ) we define the Variety of Nonreduced Sums of Powers as a subscheme with the reduced structure of VSP( f ) corresponding to the locus of nonreduced schemes. In other words: or equivalently it parametrizes smoothable, nonsmooth apolar schemes. As a moduli space VSP comes with a universal family π : where the fiber over a given point of VSP( f ) equals the corresponding apolar scheme.
Let us describe how Proposition 2.6 may be used to find the boundary of SSP inside the VSP in special cases. Let f ∈ S 2d (V * ) have a nondegenerate middle catalecticant. Further assume that the rank of the catalecticant equals the generic rank in S 2d (V * ).
1. Assume we are given the universal family π : We note that set-theoretically π(B ∩ F) coincides with VNSP( f ). The construction above follows from the fact that a real decomposition can change to a complex one only by passing through a nonreduced scheme. As a consequence we obtain the following lemma that will allow us to find defining equations of ∂ alg SSP( f ).

Lemma 2.8 Suppose that SSP( f ) = ∅. If the top dimensional component of π(B∩F)
is irreducible, then its reduced structure coincides with ∂ alg SSP( f ).
In principal, this method could provide the description for the boundary of the SSP inside the VSP for quaternary quartics. Unfortunately, in this case the VSP, which is a 5-fold eludes an explicit description. Ranestad communicated to the authors that this is one of the most interesting outstanding problems on VSP's.

Quadrics
We now apply Proposition 2.6 to explicitly obtain the boundary of SSP for quadrics f in up to n ≤ 5 variables, i.e. in all cases when VSP is smooth. For n = 2, 3 this was achieved in Michałek et al. (2017). It that case, as the codimension of the apolar ideal f ⊥ was at most three, one could apply the classical results of Buchsbaum-Eisenbud on resolutions of Gorenstein schemes (Buchsbaum and Eisenbud 1977) and define the boundary by an appropriate hyperdeterminant. In case n > 3 we could still take the resolution of f ⊥ , however explicit results using this technique seem much harder. Instead we follow Proposition 2.6.

Lemma 2.9
Fix two real forms f 1 , f 2 , both with nonempty SSP's. Suppose f 1 can be obtained from f 2 by a complex change of coordinates. If the top dimensional Proof The isomorphism between f 1 and f 2 provides an isomorphism between their VSP's. We show that this isomorphism is also an isomorphisms of the algebraic boundaries. By the assumption the boundary is nonempty in both cases, hence of codimension one. However, it has to be contained in VNSP, which is irreducible and preserved by the isomorphism.
Remark 2.10 Note that we need the assumption that SSP's are nonempty. We have SSP(x 2 + y 2 ) = VSP(x 2 + y 2 ), hence the algebraic boundary is empty. On the other hand, a quadric with a different signature (but of course isomorphic over C) satisfies SSP(x 2 − y 2 ) VSP(x 2 − y 2 ) and the algebraic boundary consists of two points.
We will apply Lemma 2.9 to quadrics of different signature. Our aim is to describe the algebraic boundary of the SSP.

Ternary quadrics
The case of ternary quadrics is well-understood [Michałek et al. 2017, Section 2]. We use it as a warm-up. Fix f = x 1 x 3 + x 2 2 , which up to real isomorphism is the only indefinite quadric. The VSP( f ) is a smooth Fano 3-fold V 5 -quintic del Pezzo threefold-admitting a realization as an intersection G(3, 5) ∩ P 6 . The boundary ∂ alg SSP( f ) is given by a special hyperdeterminant, which turns out to be a degree 20 polynomial in 6 variables with 13956 terms.
We show how this polynomial simplifies, if we work in a local affine patch of VSP( f ) described in Ranestad and Schreyer (2013). Indeed, computing π(B ∩ F) in Appendix 4.1 we obtain a quartic surface: Its singular locus is a smooth curve-a complete intersection of a quadric and cubic surface. It corresponds to the locus where the apolar scheme is local, i.e. all three points coming together. As the quartic is irreducible it follows that it coincides with ∂ alg SSP( f ).
Three pictures of the affine patch of ∂ alg SSP( f = x 1 x 3 + x 2 2 ). The threedimensional ambient affine space represents the quintic Fano threefold-moduli of schemes apolar to f of length three. The surface corresponds to those schemes that are supported at two points: one smooth, one with nonreduced structure Spec C[x]/(x 2 ). The curve in the singular locus represents apolar schemes isomorphic to Spec C[x]/(x 3 ). The surface divides the threefold into two regions: one is the SSP( f ) corresponding to fully real decomposition, the other corresponds to decompositions where one linear form is real and the other two are conjugate.

Quaternary quadrics
The following theorem computes the boundary of the SSP relating its geometry to the Hilbert scheme and the geometry of the apolar schemes. Ranestad and Schreyer Ranestad and Schreyer (2000) provide an explicit local description of the variety VSP( f ). We will be working locally on such affine patches.
and equals π(B ∩ F). 2. The singular locus of Y has two four dimensional components over Q: Y 3,1 and Y 2,2 . The general point of Y 3,1 (resp. Y 2,2 ) corresponds to a nonreduced scheme with two support points; one of which has multiplicity 3 (resp. 2), the other has multiplicity 1 (resp. 2). 3. Over C, Y 2,2 however has two irreducible components; Y 2,2 and Y 2,2 intersecting along surface (identified later with Y 4,sing ). 4. The two components Y 3,1 , Y 2,2 intersect along a three dimensional threefold Y 4 , irreducible over Q. It is the singular locus of Y 3,1 . 5. Over C, Y 4 has two components, corresponding to intersections of Y 2,2 and Y 2,2 with Y 3,1 . By Lemma 2.9 the same result holds for all indefinite, full rank quadrics.
Proof By [Ranestad and Schreyer 2013, Theorem 1.1] VSP( f ) is a smooth, six dimensional variety. Ranestad and Schreyer provide an explicit local description of this variety and the universal family F in a Macaulay 2 package VarietyOfPolarSimplices.m2 (http://www.math.uni-sb.de/ag/schreyer/home/computeralgebra.htm). We extend their defining equations by ( f ), obtaining B ∩ F ⊂ VSP( f ) × P 3 . Eliminating the variables corresponding to P 3 automatically is not possible (due to a large ambient dimension of the VSP). However, one may find automorphisms of the (local) description of the VSP that reduce the number of variables. Afterwards we perform elimination, obtaining explicitly the defining equation of π(B ∩ F) (defined over Q)see Appendix 4.2. We check that it defines a prime ideal, over C, as follows-details of implementation are presented in Appendix 4.2.1. We fix four linear forms defined over Q and intersect them with the π(B ∩ F) obtaining a curve C. It is enough to show that C is irreducible over C, as π(B ∩ F) was a hypersurface-in particular equidimensional. We project, not changing the degree, until C becomes a plane curve, that is defined by g C . We prove it is irreducible as follows. We consider all possible factorizations of g C into a product of degrees 4 + 4, 3 + 5, 2 + 6, 1 + 7 with coefficients given by variables. Comparing all coefficients we obtain ideals, that equal the whole ring. By the previous discussion this proves the statement 1.
The computation of the singular locus of Y and its decomposition over Q is now straightforward. The two components Y 3,1 and Y 2,2 are distinguished by the dimension of their singular locus, which is respectively 3 and 2. There are several ways to prove that Y 3,1 corresponds to schemes of type (3, 1) and Y 2,2 of type (2, 2). One can restrict the family F to Y 3,1 and intersect it with ( f ). By Proposition 2.6 in case of schemes of type (3, 1) these yields a family of local (i.e. supported at one point) schemes (as only one of the points was nonreduced), while the schemes of type (2, 2) provide two distinct support points. Statements 2 and 4 follow.
To prove the statement 3 we project Y 2,2 obtaining a hypersurface H of degree 4. If Y 2,2 were irreducible, then H would have to be irreducible. However, the equation defining H decomposes as a product of two quadrics.
It may seem not clear why Y 4 -corresponding to all four points coming together-allows further degeneration. The reason is the geometry of the punctual Hilbert scheme of schemes of length four. The smoothable component is irreducible and consists of alignable schemes with a general point corresponding to the aligned scheme Spec C[x]/(x 4 ). However, not all schemes are aligned. The scheme Spec C[x, y]/(x 2 − y 2 , x y) is local, Gorenstein, but has a two dimensional tangent space. A very explicit degeneration of four points to this scheme is presented e.g. in [Landsberg and Michałek 2017, p. 11]. The main difference between C[x]/(x 4 ) and C[x, y]/(x 2 − y 2 , x y) is that in the first case the four points degenerate on a line, while in the second they come from distinct, linearly independent directions.
More explicitly, for a scheme corresponding to a point of Y 4 one can find an explicit projection to a line that preserves its degree-a feature possible only in case of Spec C[x]/(x 4 ). Further, one can find explicitly points of Y 4,sing and directly prove that the associated schemes are isomorphic to Spec C[x, y]/(x 2 − y 2 , x y).
The last statement follows by computation presented in Appendix 4.2.2. Precisely, we show that one equation in the ideal of Y 4 is of the form f 2 1 + f 2 2 and I (Y 4 ) + ( f 1 , f 2 ) = I (Y 4,sing ).

Remark 2.12
In the given embedding, Y 2,2 is of degree 4. Thus, Y 2,2 and Y 2,2 are of degree 2 each. As they are also of codimension 2, each of them should be defined by a quadric in a hyperplane section. Using arithmetics in a finite field we have checked that they are smooth.

Remark 2.13
As the locus in the Hilbert scheme corresponding to schemes isomorphic to C[x]/(x 4 ) is irreducible, it is surprising that Y 4 is reducible over C. It would be interesting to know if it could happen for generic form f in n + 1 variables of rank r that: • the codimension one subvariety of the smoothable component of the Hilbert scheme of r points in P n corresponding to nonreduced, Gorenstein schemes is irreducible, • the boundary ∂ alg SSP( f ) is reducible and positive dimensional.
There is always a divisor in the smoothable component of the Hilbert scheme, which is corresponding to r − 2 reduced points and a multiplicity 2 point. In some cases, however, there are other divisors corresponding to nonreduced Gorenstein schemes. For more information we refer the reader to [Iarrobino 1984, Theorem 1] and further consequences of this fact are discussed in [Buczyński et al. 2015, Example A.17], [Landsberg and Michałek (2017), p. 7]. However, even in such a case we do not know if ∂ alg SSP( f ) could be reducible (although we expect it).
Remark 2.14 All apolar schemes that we obtain must be Gorenstein by [Buczyńska and Buczyński 2014, Lemma 2.3]. Indeed, otherwise there would exist a Gorenstein scheme of smaller length that would be apolar to the form. In this range (e.g. in ambient dimension at most three or if length is at most ten) all Gorenstein schemes are smoothable. Hence, the form would have smaller border rank, i.e. would not be generic.

Quinary quadrics
In this case the VSP is a ten dimensional smooth variety. As before we may work locally in an affine patch described by Ranestad and Schreyer. However, in this case explicit computation of π(B∩F) is much harder, as the elimination of variables is both time and memory consuming. Instead we apply Lemma 2.8. We proceed as follows: 1. Eliminate by hand possibly many variables, using explicit isomorphisms of the ambient affine space of the affine patch of B ∩ F. 2. Consider the compactification C of the affine patch of B ∩ F in a projective space. 3. Fix a P 2 in a projective space that is a compactification of the affine patch of VSP. 4. Restrict the family C to the given P 2 and project obtaining C . 5. Check that C is a reduced, irreducible curve of degree 10.
As C is reduced, irreducible and of degree 10, it follows that the maximal dimensional component of π(B ∩ F) must be irreducible and defined by an unique polynomial of degree 10. A simplified, explicit computation is presented in Appendix 4.3. The code uses finite fields to make the computations fast (performing formal computation takes several hours, but is possible with the same code after changing the field). By Lemma 2.8 we obtain the following proposition.
We further note that the curve C representing ∂ alg SSP( f ) has 30 singular points.

Quaternary cubics
A general quaternary cubic f ∈ R[x 1 , x 2 , x 3 , x 4 ] has complex rank R(4, 3) = 5 and the decomposition is uniquely given by (1) which was claimed by Sylvester in 1851, and proved by Clebsch in 1861Clebsch (1861. It is known as Sylvester Pentahedral Theorem. Similarly to , Section 5], we propose the following algorithm to compute the five linear forms i .

Algorithm 3.1 I nput: A general quaternary cubic f . Out put: The decomposition
(1).

Compute a vector
To prove the correctness of this algorithm, we need the following lemma.

Lemma 3.2 For a general quaternary cubic f
, the minimal free resolution of its apolar ideal is given by Further, the map M can be represented by 6 × 10 matrix with two blocks where L is a 6 × 5 matrix with linear forms and Q is a 6 × 5 matrix with quadric forms. Moreover, there exists a row operation that makes the last row of L equal to 0. It distinguishes a unique five dimensional subspace of quadrics in f ⊥ with five linear syzygies.
Proof Since the apolar ideal of a quaternary cubic is Gorenstein, codimension four with six quadric generators, as in Reid (2015), the resolution has the form (2). So we only need to show that the matrix representation M has the form (3).

To prove this, consider a family of cubics in
where 16 of its coefficients are fixed. After applying G L R (4)-action on this family, it is dominant on the space of cubics, which we prove using an easy computation of a Jacobian. Hence, it is enough to prove the lemma holds for this family.
Hence, the first syzygy matrix M has five linear columns and five quadric columns, i.e.

M = [ L | Q ]
where L is a 6 × 5 matrix with linear forms and Q is a 6 × 5 matrix with quadric forms. Also, by multiplying M on left by the invertible matrix By applying the column operation given by the right multiplication of this matrix to the generators [ g 1 . . . g 6 ], we get the new generators [ c 6 g 1 − c 1 g 6 c 6 g 2 − c 2 g 6 . . . c 6 g 5 − c 5 g 6 , g 6 ].
The first five of them have five linear syzygies given by the five columns of U L. Moreover, choosing right basis as in Reid (2015) and applying Lemma 3.2, the second syzygy matrix N also has the form and the first row of L is 0.
From this, the ideal J generated by the first five generators c 6 g 1 − c 1 g 6 , c 6 g 2 − c 2 g 6 ,. . ., c 6 g 5 − c 5 g 6 has the minimal free resolution of the form where P is given by the 5 × 5 upper submatrix of L as we seen above.
Since the resolution has the form (5), the ideal J has dimension 0 and degree 5, and it is also contained in the apolar ideal f ⊥ . Further, by Lemma 3.2, the five chosen quadrics span the unique subspace contained in f ⊥ with this resolution.
Since there is the unique decomposition (1) for the given general cubic f by Clebsch (1861), the linear forms in the decomposition define an ideal K ⊂ f ⊥ of five points in P 3 . Since K is Gorenstein of codimension 3, by the Buchsbaum-Eisenbud structure theorem, it has the free resolution of the form It means that the ideal J we obtained by the Algorithm 3.1 is actually same as the ideal K , which defines the five linear forms in (1), as we wanted.
Using this algorithm, we can easily check whether the given quaternary cubic has real rank 5 or not. Namely, one computes the unique decomposition (1) and checks whether it is real. The real rank boundary can be obtained as in the following proposition. This proposition confirms [Michałek et al. 2017, Conjecture 5.5] for quaternary cubics.
Proof The parametrization defines a unirational variety Y in P 19 . The Jacobian of this parametrization is found to have corank 1. This means that Y has codimension 1 in P 19 . Hence Y is an irreducible hypersurface, defined by a unique (up to sign) irreducible homogeneous polynomial in 20 unknowns with rational coefficients. Let g be a real cubic with the form (7) that is a general point in Y . As goes to 0, the real cubics ( 4 + 5 ) 3 − 3 4 and (i 4 + 5 ) 3 + (−i 4 + 5 ) 3 converge to the cubic 2 4 5 in P 19 . It means that any small neighborhood of g in P 19 contains cubics of real rank 5 and cubics of real rank > 5. This implies that Y lies in the real rank boundary ∂ alg R 4,3 . Since Y is irreducible and codimension 1, it follows that ∂ alg R 4,3 exists and has Y as an irreducible component.
Using the algorithm 3.1, we can exactly compute the degree of the hypersurface ∂ alg R 4,3 which is 40. This is done as follows. First, fix the field K = Q(t) with a new variable t. We fix two cubic f 1 and f 2 in Q[x 1 , x 2 , x 3 , x 4 ] 3 , and run the algorithm 3.1 for Step 4 returns an ideal J in K [x 1 , x 2 , x 3 , x 4 ] that defines 5 points in P 3 over the algebraic closure of K . By eliminating each of two variables of {x 1 , x 2 , x 3 , x 4 }, we obtain six binary forms of degree 5 such that their coefficients are degree 15 polynomials in t. The discriminant of each binary form is a polynomial in Q[t] of degree 120 = 8 * 15. The greatest common divisor of these discriminants is a polynomial (t) of degree 40. We can check that By definition, is a homogeneous polynomial with integer coefficients, irreducible over Q, in the 20 coefficients of a general cubic f . Its specialization ( f 1 + t f 2 ) is a non-constant polynomial in Q[t], of degree deg(X ) in t. That polynomial divides (t) because Y lies in the real rank boundary. Since the latter is also irreducible over Q, we conclude that ( f 1 + t f 2 ) = γ · (t), where γ is a nonzero rational number. Hence has degree 40. We conclude that deg(Y ) = 40, and therefore Y = ∂ alg (R 4,3 ).

Quaternary quartics
Proposition 3.5 One component of the boundary ∂ alg R 4,4 equals the dual of a 24 dimensional variety of quartic symmetroids in P 3 , that is, the surfaces whose defining polynomial is the determinant of a symmetric 4 × 4-matrix of linear forms.
Proof By Hilbert's classification we know that for quaternary quartics the cone of positive forms P 4,4 strictly contains the cone of sums of squares 4,4 . Hence, the cone Q 4,4 of sums of (arbitrary many-c.f. Remark 3.7) 4-th powers of linear forms, which is the convex hull of the Veronese and equals the dual P * 4,4 , is strictly contained in the cone * 4,4 of forms with psd catalecticant. The boundary of Q 4,4 has two components. One (that is not interesting for our purposes) is the boundary of * 4,4 given by the determinant of the catalecticant. The other one, which we denote by B, is the dual of the Zariski closure of the extremal rays of P 4,4 \ 4,4 -the variety of quartic symmetroids in P 3 by [Blekherman et al. (2012), Theorem 3].
The forms in * 4,4 \Q 4,4 are obviously of real rank greater than 10, as they are not sums of powers by definition, and any other presentation would contradict the signature of the catalecticant. Further a generic point of B has real rank at most 10 by [Blekherman et al. (2012), Proposition 7]. Changing the linear forms of the Waring decomposition of quartics in B we obtain a Zariski dense set of forms of real rank at most 10, intersecting B in a relatively open set, which proves the proposition.
The following lemma is well-known to experts. We present a sketch of a proof based on [Blekherman et al. 2013, Lemma 4.18].
Lemma 3.6 For all positive integer k, the following set is a closed cone. In other words, the cone of forms whose positive rank is less than or equal to k is closed.
Proof Let = {(λ 1 , . . . , λ k ) ∈ R k ≥0 | k i=1 λ i = 1} be a simplex and let O be the unit sphere in V * . Consider the image C of the map: Clearly, C is compact and 0 / ∈ C. Thus, the cone over C is closed and it coincides with C k,2d .
Remark 3.7 There exists a nonempty open set inside Q 4,4 whose elements have real rank strictly greater than 10. Indeed, the quartic (x 2 1 + x 2 2 + x 2 3 + x 2 4 ) 2 has real rank 11 by [Reznick (1992), Proposition 9.26] and hence belongs to Q 4,4 \C 10,4 . The latter is open by Lemma 3.6. Further, every element of this set has real rank strictly greater than 10. Indeed, since every element of Q 4,4 has definite 10 × 10 catalecticant matrix, it cannot have decompositions whose coefficients have distinct signs.
One can easily prove that quaternary quartics f of signature (9, 1) have rank at least 11 if ( f ) has no real points. This and Remark 3.7 motivate the following conjecture.