Inequalities for Moment Cones of Finite-Dimensional Representations
Michèle Vergne, Michael Walter
Introduction
The study of the convexity properties of the moment map and of its image has a long history in mathematics, starting from Schur and Horn’s observation that the diagonal entries of a Hermitian matrix are always contained in the convex hull of the spectrum Schur (1923); Horn (1954); cf. Kostant (1973). More generally, Atiyah and Guillemin-Sternberg have shown that, for any torus action on a compact, connected Hamiltonian manifold, the image of the moment map is a convex polytope, called the moment polytope Atiyah (1982); Guillemin and Sternberg (1982a). It can be explicitly computed as the convex hull of the images of torus fixed points. For non-abelian groups, the image of the moment map is no longer convex. Instead, Kirwan’s celebrated convexity theorem asserts that, for the action of an arbitrary compact, connected Lie group on a compact, connected Hamiltonian manifold, the intersection of the image of the moment map with a positive Weyl chamber is a convex polytope Kirwan (1984a). This is the correct generalization of the moment polytope to non-abelian group actions. Mumford has given a different proof in the case of projective subvarieties, which relies on a concrete description in terms of the decomposition of the homogeneous coordinate ring into irreducible representations Ness and Mumford (1984); cf. Guillemin and Sternberg (1982b); Brion (1987). However, no effective general methods are known for the computation of these polytopes (in contrast to the case of torus actions).
Our main contribution is a clean algebraic description of the moment cone in terms of finitely many linear inequalities. To state the result, let denote the Lie algebra of a maximal torus of , a positive Weyl chamber, and the (complex) Lie algebra representation induced by (see section 2 below for precise definitions). We shall say that is a Ressayre element if 1) the hyperplane is spanned by weights of and 2) there exists a vector annihilated by such that the “tangent map” at ,
is an isomorphism; here, denotes the direct sum of all negative eigenspaces of and the sum of all root spaces for negative roots such that (definition 3.14). Note that there are only finitely many Ressayre elements for any given representation . In section 3, we will prove the following result:
The moment cone for the -action on is given by
To prove theorem 1.1, we show that any facet corresponds to a Ressayre element by studying the moment map, which is quadratic, locally up to second order. To show that, conversely, any Ressayre element determines a valid inequality, we use Mumford’s description of the moment cone as in Ressayre (2010a). Indeed, our notion of a Ressayre element is closely related to Ressayre’s notion of a dominant pair. We note that the description in theorem 1.1 will typically contain redundancies. Thus our result differs from Brion (1999), where the non-trivial or “general” faces of the moment polytope are characterized precisely at the cost of requiring a recursive strategy for their computation.
If is a Ressayre element then the domain and codomain of the tangent map (1.1) necessarily have the same dimension, i.e., . We call this the trace condition. Moreover, note that the determinant of (1.1) (with respect to any fixed pair of bases) is a non-zero polynomial in , the zero eigenspace of . In fact, is a canonical (up to scalar multiplication) lowest weight vector for the action of , the centralizer of the torus generated by , on the space of polynomials on . This implies the following result in section 4, which we call the Horn condition:
is an element of the moment cone . In fact, is a lowest weight vector of weight .
Here, denotes the set of weights of and denotes the set of negative roots of . By applying theorem 1.1 to the lower-dimensional scenario, the Horn condition can be explicitly stated as a set of linear inequalities that have to be satisfied by .
Tangent maps and their determinants have been studied in great generality by Ressayre and Belkale from an algebro-geometric point of view Ressayre (2010a, b); Belkale (2010), and our theorems 1.1 and 1.2 can also be deduced from their results. In these works, the non-vanishing of the determinant has been in turn been translated into a cohomological condition. In contrast, we propose that, for the purposes of computing moment cones explicitly, it can be useful to instead test the non-vanishing of the determinant directly—either symbolically, which is easily possible in small dimensions, or numerically by using fast algorithms for polynomial identity testing, as we discuss in section 5 below. The challenge imposed by higher dimensions is rather in finding additional a priori constraints on the facets of the moment cone.
Moment Cones of Finite-Dimensional Representations
For , whose complexification is , we will always use the maximal torus that consists of the diagonal unitary matrices of unit determinant. Its Lie algebra will be denoted by . We will use as positive roots the with , and abbreviate the (positive) roots by and , respectively. Finally, we write for the coadjoint orbit through a highest weight .
Now let a representation of on a finite-dimensional Hilbert space that is equipped with a -invariant Hermitian inner product , which we take to be antilinear in the first argument. We will oftentimes use Dirac’s notation for and . We denote by the induced Lie algebra representation, by the set of weights and write the weight space decomposition as . The -action on admits a canonical (up to conventions) moment map, defined by
for all and . Here and in the following, we write for the duality pairing. The map is indeed a moment map in the sense of symplectic geometry: it is -invariant and satisfies the basic identity
for all , where denotes the symplectic form that we will use for . The moment cone then is defined as intersection of the moment map image with the positive Weyl chamber,
where denotes the space of polynomials on . A representation occurs in if and only if there exists a polynomial of weight that is invariant under the action of the lower unipotent subgroup . We shall call such a polynomial a lowest weight vector in .
(and is a pointed cone with base ). Thus we may equivalently study moment cones of representations or moment polytopes of the corresponding projective spaces.
Throughout this paper, we shall always work with moment cones (but see Walter (2014) for an exposition from the projective point of view). We shall moreover assume that the moment cone is of maximal dimension, i.e., . This is the case if and only if there exists a vector with finite stabilizer.
Facets of the Moment Cone
Like any polyhedral cone, the moment cone can be described by finitely many linear inequalities . Since we have assumed that is of maximal dimension, its facets are of codimension one in and their inward-pointing normal vectors may be identified with the defining linear inequalities of the moment cone. Since the moment cone is obtained by intersecting with the positive Weyl chamber, which itself is a maximal-dimensional polyhedral cone, some of the facets of can be subsets of facets of , and we shall call those the trivial facets of the moment cone:
A facet of the moment cone is trivial if it corresponds to an inequality of the form for some positive root . Otherwise, the facet is called non-trivial.
Non-trivial facets have also been called “general” in the literature Brion (1999). We record the following straightforward observation:
Any non-trivial facet of meets the relative interior of the positive Weyl chamber.
We first consider the moment map , defined as in (2.2) for the action of the maximal torus . Let be the decomposition of into weight spaces and let be a vector decomposed accordingly as . Then has the following concrete description:
Observe that is a conic combination of weights. It follows that the “abelian” moment cone of is precisely equal to the conical hull of the set of weights; it is maximal-dimensional since it contains . More generally, if is a subset of weights and then . For the next lemma recall that a critical point of a smooth map is a point where the differential df\big{|}_{m} is not surjective; a critical value is the image of a critical point. Then the following is well-known (e.g., (Christandl et al., 2014, Remark 4.14)):
The set of critical values of is equal to the union of the codimension-one conic hulls of subsets of weights.
Let with weight decomposition . By (3.1), is a conic combination of weights in . By the moment map property (2.3) and non-degeneracy of the symplectic form, is a critical point if and only if there exists such that (Guillemin and Sternberg, 1982a, Lemma 2.1), i.e., if and only if for all . It follows that is a critical point if and only if the conic hull of is of positive codimension.
In particular, any critical value is contained in a codimension-one conic hull of weights, since we may always add additional weights. Conversely, if is a subset of weights that spans a conic hull of codimension one then consists of critical values. ∎
We now derive a basic necessary condition that cuts down the defining inequalities of the moment cone to a finite set of candidates (cf. (Christandl et al., 2014, Remark 3.6)).
An element is called admissible if the linear hyperplane is spanned by a subset of weights in .
The notion of admissibility is invariant under the action of the Weyl group .
Let be an inequality corresponding to a non-trivial facet of the moment cone. Then is admissible.
By lemma 3.2, the intersection of with the interior of the positive Weyl chamber is non-empty. Each point in this intersection is a critical value for , hence of , and therefore according to lemma 3.3 contained in a linear hyperplane spanned by a subset of weights. Since this is true for all points in the intersection, which contains the relative interior of the facet, it follows that the facet is in fact contained in a single such hyperplane. ∎
2 Description of the Moment Cone by Ressayre Elements
Let us now fix an inequality corresponding to a non-trivial facet of the moment cone. Let be a preimage of a point on the facet . In the proof of lemma 3.5, we have used that is a critical point of (equivalently, that ) to gain information on the set of possible facets. To study the function in the vicinity of such a critical point it is natural to consider the Hessian, which is the quadratic form
For tangent vectors generated by the infinitesimal action of , we have the formula
where we have used that . We now decompose
where is the sum of the eigenspaces of the Hermitian operator with eigenvalue less than , etc. Then it is plain from (3.2) that the index of the Hessian , i.e., the dimension of a maximal subspace on which the quadratic form is negative definite, is equal to twice the complex dimension of .
We will now deduce a second formula for the index by observing that the Hessian is necessarily positive semidefinite on the subspace of those tangent vectors that are mapped to by the differential of the moment map. To see this, consider a curve with , and for all (such a curve can always be constructed by using the symplectic cross section (Guillemin and Sternberg, 1984b, Theorem 26.7)). Then, since and d(\mu_{K},H)\big{|}_{\psi}\equiv 0,
But , since and is a valid inequality for the moment cone. Together, this shows that, indeed, .
The subspace of all such can be computed in a different way. For this, let denote the sum of the root spaces of the compact Lie algebra. It follows from the moment map property (2.3) that
The Hessian is positive semidefinite on the symplectic complement
of .
In fact, it is well-known that is a symplectic subspace of (see, e.g., (Guillemin and Sternberg, 1982a, Lemma 6.7)). To see this, note that the restriction of the symplectic form to is given by
and can therefore be identified with the Kirillov-Kostant-Souriau symplectic form on the coadjoint orbit through . As a consequence, we have the decomposition .
The Hessian is block-diagonal with respect to the decomposition .
For all and we have that
since and therefore . ∎
The tangent map , is injective.
The stabilizer of the coadjoint action of at any is , while . As , the claim follows from this and -equivariance of the moment map. ∎
The index of the Hessian is equal to twice the number of positive roots such that .
Since is positive semidefinite on (lemma 3.6) and block-diagonal with respect to the decomposition (lemma 3.7), it suffices to compute the index of on . For this, recall from (3.3) that for all . Since the tangent map is injective (lemma 3.8), we may instead consider the form
on . Now observe that is block-diagonal with respect to , since for all and , , while . Therefore, it suffices to compute the index on a single root space . For this, define the “Pauli matrices” and , which satisfy the commutation relations etc. Then and form a basis of and
Likewise, , while . Since , we conclude that the index of is equal to twice the number of positive roots with . ∎
Since we have already seen above that the index of the Hessian is also equal to twice the complex dimension of , we obtain the following result, which can also be extracted from (Brion, 1999, Theorem 2):
If defines a non-trivial facet of the moment cone then
We now study the complexified group action. To this end, we consider the Lie algebra of the negative unipotent subgroup, which plays a role analogous to for vectors that are mapped into the positive Weyl chamber (compare the following with lemma 3.8).
The tangent map is injective.
Let be an arbitrary element in . Since , we find that with (the Cartan involution of ). Therefore,
so that by using we find that
In contrast to lemma 3.8, which continues to hold true if is mapped to the relative interior of a different Weyl chamber, it is important in lemma 3.11 to choose the negative unipotent subgroup (relative to the choice of positive Weyl chamber). For example, consider an irreducible -representation with highest weight and highest weight vector . Then and the “lowering operators” in indeed act injectively. On the other hand, the “raising operators” in the positive nilpotent Lie algebra annihilate the highest weight vector (by definition).
We now decompose the Lie algebra similarly to (3.4),
where is the sum of the complex root spaces with negative -weight , etc. We observe that corollary 3.10 can be equivalently stated as
Note that . Thus we obtain the following important result:
Let such that is a point on a non-trivial facet of the moment cone corresponding to the inequality . Then and the tangent map restricts to an isomorphism
The fact that is just a reformulation of . By the preceding discussion, the tangent map is well-defined as a map from to ; it is injective by lemma 3.11 and surjective since the dimensions agree according to (3.5). ∎
We now prove a partial converse to proposition 3.12, inspired by the argument of Ressayre Ressayre (2010a).
Suppose there exists such that the tangent map
is surjective. Then is a valid inequality for the moment cone.
Its differential at is the linear map
The assumption implies that this map is surjective. It follows that contains a small Euclidean ball around . In particular, any -invariant polynomial that is zero on is automatically zero everywhere on .
We now prove the inequality. By the description of the moment cone in (2.4), it suffices to show that for all highest weights such that . Recall that the highest weight of is , where is the longest Weyl group element that flips the positive and negative roots. Consider a lowest weight vector, i.e., a polynomial that is a weight vector of weight and invariant under the action of . Then and the restriction of to is non-zero by our discussion above. But this restriction is an element of , the space of polynomials on . Since all -weights in are non-positive, it follows that , as we set out to prove. ∎
We remark that proposition 3.13 holds unconditionally without any assumption on the dimension of the moment cone . We summarize our findings in the following definition and theorem that we had already advertised in the introduction.
An element is called a Ressayre element if
is admissible, i.e., the linear hyperplane is spanned by a subset of weights in , and
there exists such that the map
This follows directly from lemmas 3.5, 3.12 and 3.13. ∎
Theorem 1.1 gives a complete description of the moment cone of an arbitrary finite-dimensional -representation (under the assumption that is of maximal dimension). The set of inequalities thus obtained may still be redundant (i.e., not all inequalities necessarily correspond to facets of the moment cone). In contrast, Ressayre’s well-covering pairs Ressayre (2010a) characterize the facets of the moment polytope precisely. In our language, his condition amounts to requiring that that the generic fiber of the map is a point (as opposed to only requiring that the map be dominant). Our characterization is also related to (Brion, 1999, Theorem 2), which uses algebraic geometry to characterize non-trivial faces of arbitrary codimension. Unlike proposition 3.13, it relies on an assumption about lower-dimensional moment polytopes, which can in principle be obtained recursively.
Generalized Trace and Horn Condition
We now extract two useful necessary conditions that have to hold for any Ressayre element , and therefore for any non-trivial facet of the moment polytope. We will see in section 7.1 below that they are a generalization of the classical Horn inequalities, which justifies our terminology.
The first condition, which we call the trace condition is the observation that the domain and range of the tangent map (3.6) necessarily have to agree, i.e.,
We note that the right-hand side of (4.1) is invariant under the action of the Weyl group. This suggests that we first compute the dominant admissible and then determine those which satisfy the trace condition (4.1). For this, we will need the following well-known lemma.
Fix any such that for all positive roots , and set . Let be a Weyl group element such that . For all positive roots and small enough,
if and only if implies that . In other words, if and only if for all with . Since is regular, this immediately shows that such Weyl group elements exist and are unique. What is more, regularity also implies that
We obtain the following useful corollary:
The second condition, called the Horn condition, is based on the observation that for any Ressayre element the determinant polynomial
is non-zero (we take the determinant with respect to any fixed pair of bases). This can be understood to imply a statement about a smaller moment cone. To see this, let denote the identity component of the centralizer of the torus generated by , i.e., the connected subgroup of with Lie algebra ; denote by the corresponding negative unipotent subgroup and by the maximal compact subgroup. Then acts on and thus on the space of polynomial functions.
is an element of the moment cone . In fact, is a lowest weight vector of weight .
In view of (2.4), it suffices to argue that the determinant polynomial is a lowest weight vector in of weight .
To see this, fix a basis of such that each is a weight vector of weight , and denote by the negative roots with , so that is a basis of . Then the determinant polynomial with respect to this basis can be written as
where we write for the canonical homomorphism induced by a linear map . It follows that for any and ,
Computation of Moment Cones
Theorem 1.1 reduces the computation of the moment cone of an arbitrary finite-dimensional representation to an enumeration of all Ressayre elements, which in principle is straighforward: Since there are only finitely many weights, the admissibility condition cuts down the number of possible inequalities down to a finite list of candidates, and for each such candidate , the isomorphism condition can be easily checked. Indeed, we only need to verify the trace condition (4.1) and that the determinant polynomial (4.2) is non-zero. In this way, we obtain a deterministic algorithm to compute the moment cone for an arbitrary representation that can easily be implemented in a computer program onl . We remark that checking whether the determinant polynomial is non-zero can be sped up by using a fast probabilistic algorithm for polynomial identity testing, e.g., based on the Schwartz-Zippel lemma.
In practice, naively enumerating all admissible hyperplanes by considering all -element subsets of quickly becomes infeasible as one considers representations of larger dimensions. In this case, it is useful to first determine the admissible elements up to the Weyl group and to impose further necessary conditions by a more refined analysis of the representation at hand. Then corollary 4.2 can be used to obtain directly only those candidates that satisfy the trace condition. In this way, we may often cut down the number of candidates substantially for which we need to check that is non-zero. In section 6 below, we illustrate this analysis in the case of the one-body quantum marginal problem.
Quantum Marginal Problem and Kronecker Cone
The state of a quantum system is specified by a unit vector in a Hilbert space , which we will always assume to be finite-dimensional. Quantum systems composed of several distinguishable particles are described by the tensor product of the Hilbert spaces of their constituents, . The state of the -th subsystem can be described by the reduced density matrix , which is the unique positive semi-definite operator on such that
for all Hermitian operators on (it is also called the partial trace of ). The fundamental one-body quantum marginal problem asks which can arise as the reduced density matrices of a quantum state Klyachko (2004); Christandl and Mitchison (2006); Daftuar and Hayden (2004); Christandl et al. (2007); Walter (2014). Equivalently, it asks for the compatibility conditions that the have to satisfy in order for there to exist a global state . For two subsystems, it is a straightforward consequence of the singular value decomposition that and are compatible if and only if they have the same non-zero eigenvalues (including multiplicities). In general, however, the problem is much more involved; it can be shown that it is a strict generalization of the problem of computing the Horn cone of section 7 Klyachko (2004); Christandl et al. (2012); Walter (2014) and it has been solved in Klyachko (2004); Daftuar and Hayden (2004) by using similar methods Klyachko (1998); Berenstein and Sjamaar (2000). However, a concrete description akin to the Horn inequalities is still elusive.
By Schur-Weyl duality and Mumford’s description (2.4), the moment cone can equivalently be defined in terms of the representation theory of the symmetric group Christandl and Mitchison (2006); Klyachko (2004); Christandl et al. (2007, 2014): We have
where , , and vary over the set of Young diagrams with the same number of boxes and no more than , , and rows, respectively, and where denotes the Kronecker coefficient of the symmetric group, i.e., the multiplicity of the invariant subspace in the corresponding triple tensor product of irreducible representations. We will henceforth refer to as the Kronecker cone.
It will be convenient to assume without loss of generality that . In this case, is maximal-dimensional (see corollary A.2 in the appendix) and our method is directly applicable. In fact, we shall see below that all other cases can be reduced to the case . According to theorem 1.1, the moment cone is thus cut out by those which are Ressayre elements, i.e., which are admissible and satisfy the isomorphism condition. Naively determining the admissible by enumerating all subsets of with cardinality , determining whether they span a hyperplane and computing the normal vector amounts to considering subsets, which rapidly becomes infeasible (e.g., for , there are over 27 billion such subsets). We therefore need to derive additional constraints to make this approach computationally feasible.
By considering , the moment cone for the tripartite problem can now be written in terms of the bipartite cone and the moment polytopes of the action of on the coadjoint -orbits :
A first consequence is that the case can always be reduced to . Indeed, (6.2) and (6.3) imply that
if . Thus the moment cone for is isometric to , while for it is obtained as a projection of the latter. In the case where , it was observed in Manivel (1997); Klyachko (2004) that any normal vector of the moment cone necessarily has a rather special form. We state this result and give a succint alternative proof that does not rely on the results of Klyachko (2004); Berenstein and Sjamaar (2000):
Let be the normal vector of a non-trivial facet of the moment cone such that . Then:
is determined by a maximal number of equations of the form ,
the components of are precisely all possible partial sums , and
By lemma 3.2, there exists a regular dominant point in the interior of this facet. In a neighborhood of , the moment cone locally looks like a half-space, so that
is not only a valid inequality that holds for all , as follows from (6.3), but in fact a facet of the moment polytope , since we have assumed that .
While determining the moment polytopes is just as hard as determining the cone , this reformulation gives us an additional insight: Recall that the Duistermaat-Heckman measure for the action of the maximal torus of on the coadjoint -orbit is a measure on the abelian moment polytope. Now let denote the projection . The affine hyperplanes through some spanned by subsets of the restricted roots , , partition the abelian moment polytope into a finite number of polyhedral chambers, and it as an immediate consequence of the Heckman formula Harish-Chandra (1957); Heckman (1982) that the measure has a polynomial density function on each chamber (cf. Boysal and Vergne (2009)). The Duistermaat-Heckman measure for the action of can be recovered by applying a number of partial derivatives to the measure for (e.g., Christandl et al. (2014)). It follows that the moment polytope , which is the support of the latter measure, is equal to a finite union of chambers; in particular, its non-trivial facets are contained in hyperplanes of the form just described.
Applied to the facet (6.5), it follows that its normal vector is defined by a maximal number of equations of the form
for some indices . This shows the first assertion. Moreover, since the facet contains for some we obtain that
Following Klyachko, we shall call any that satisfies condition 1 of lemma 6.1 an extremal edge if it is in addition dominant and primitive (in the dual of the root lattice). There are only finitely many extremal edges and we shall denote them by . The extremal edges span the extreme rays of the cubicles, which are the full-dimensional convex cones of elements cut out by a maximal set of inequalities of the form . In other words, a cubicle is defined as a set of with fixed order of the and can therefore be encoded by a standard Young tableaux of rectangular shape Klyachko (2004). This gives a straightforward way of computationally determining all extremal edges (see table 1).
We remark that standard tableaux that correspond to cubicles have also been called additive in the literature Vallejo (2014); Manivel (2014). Manivel has shown that they can be associated with minimal regular faces of the moment cone for which the corresponding Kronecker coefficients stabilize Manivel (1997, 2014). The corresponding extremal edges determine non-trivial facets of incident to this face Manivel (1997); Klyachko (2004). However, not all non-trivial facets can be obtained in this way.
Let be the normal vector of a non-trivial facet of the moment cone , where . Then or is proportional to an element in the -orbit of (i.e., an extremal edge up to and rescaling).
In view of (6.4), we may obtain finite and complete set of inequalities for by taking any facet for and restricting to its first components (that is, set to be the traceless part of and define accordingly). Such an inequality will not necessarily define a facet of , but all facets arise in this way. If we started with a non-trivial inequality then the claim follows from lemma 6.1. If we started with a trivial inequality then we either obtain a trivial inequality or . ∎
Lemma 6.1 is efficient in finding candidates for facets of , but not necessarily so for with . Indeed, while we know from the proof of corollary 6.2 that any facet of the latter can be obtained by “restriction” of a facet of the former, for each given facet there are many possible restrictions, as we need to pick a subset of components out of the components of the given . If then it can be more efficient to apply corollary 6.2 to all three of , and (e.g., in the case of ). For this we will use the following lemma:
Let be an extremal edge. Then and are each either zero or primitive in the dual of the root lattice.
Case 2: . Since the matrix with columns the roots of is totally unimodular (e.g., (Schrijver, 1986, p. 274, (18))), it follows that spans the root lattice. To find a contradiction, suppose that is neither zero nor primitive. Then we can write for some and a non-zero element in the dual of the root lattice. But then,
for all . Since the spans the root lattice, we conclude that also is an element of the dual of the root lattice. But then is not primitive, which is the desired contradiction. ∎
Let be the normal vector of a non-trivial facet of the moment cone , where . Then is proportional to an element in the -orbit of
Since the following argument is equivariant under the Weyl group , we may without loss of generality assume that is dominant
Case 1: Only on component is non-zero, say, and therefore . We may rescale such that is primitive. Then and are also primitive and therefore extremal edges by corollary 6.2.
Case 2: At least two components are non-zero, say, and . We may rescale such that is primitive. Since and are non-zero, lemma 6.3 shows that both and are individually primitive. It follows that three of , and are primitive and therefore extremal edges by corollary 6.2. ∎
By theorem 1.1, any non-trivial facet is necessarily admissible. As admissibility is a Weyl group-invariant property, we immediately obtain the following corollary:
Let be the normal vector of a non-trivial facet of the moment cone , where . Then is proportional to an element in the -orbit of
For any given , there are in general several such that is admissible (or no such at all). The sets and consist of dominant, primitive elements, and they can be easily obtained algorithmically from the sets of extremal edges.
Let be the normal vector of a non-trivial facet of the moment cone , where . Then is proportional to an element in
where denotes the set of triples such that is a triple of shuffles with respect to the components of whose lengths sum up to .
We remark that it is straightforward to computationally generate all shuffles of a given length by adapting the algorithm of Effler and Ruskey Effler and Ruskey (2003). The upshot of corollary 6.6 then is that the moment cone is cut out by those candidates in that are also Ressayre elements (together with the trivial inequalities):
We conclude this section with an illustrative example of the method.
Consider the moment cone for any . It is not hard to verify that all pairs formed from and are extremal edges (formally, we should divide by to work with primitive elements, but we refrain from doing so in the interest of readability). Therefore,
There are two possible values of which extend the above to an admissible element . The first option is . However, the dimension of is , and therefore strictly larger than the dimension of . Thus the trace condition (4.1) can never be satisfied!
We now verify that is a Ressayre element. For this, observe that we have
for any . It follows at once from theorem 1.1 that is a valid inequality for the moment polytope. We have thus obtained the well-known polygonal inequality Higuchi et al. (2003)
in a completely mechanical fashion. By symmetry, the two other inequalities obtained by permuting the subsystems , , and are also valid. We remark that the moment cone is closely related to the Horn cone for discussed in section 7 below Christandl et al. (2012); Walter (2014); Littlewood (1958); Murnaghan (1955).
2 Computational Results
To verify that a given is a Ressayre element, we have implemented a computer program that works for arbitrary representations. To compute the set of candidates in the case of the one-body quantum marginal problem, we have used the strategy explained above. In table 2 we list some results obtained by our program in the symmetric scenario , which corresponds to three quantum particles with the same number of degrees of freedom. While the moment cones and had already been computed in Higuchi et al. (2003); Bravyi (2004); Franz (2002); Higuchi (2003) using different methods, the cone had been out of reach using current methods. In contrast, our method allows us the computation of in a few minutes, since it does not rely on an intermediate computation of the higher-dimensional cone .
There are further variants of the one-body quantum marginal problem, such as for fermionic systems Coleman (1963); Ruskai (1969), where the facets amount to strengthenings of the classical Pauli exclusion principle. Our theory is also applicable to these scenarios, and it would be interesting to undertake a similar analysis of the facets that would allow the computation of the corresponding moment cones beyond what has been possible in the literature Klyachko and Altunbulak (2008). We refer to (Walter, 2014, §3) for initial investigations, where a family of fermionic inequalities has been proved for all local dimensions by using our method. We remark that we have verified numerically that the pure-state fermionic inequalities list in Klyachko and Altunbulak (2008) are correct (but not their sufficiency).
Horn Cone and Howe-Lee-Tan-Willenbring Invariants
In this section, we consider the representation of on given by
Any non-negative Hermitian matrix can be written in the form ; since the spectra of and are equal, the moment cone is equal to
where denotes the eigenvalues of a Hermitian matrix , ordered non-increasingly. We will call the Horn cone in dimensions. As proved in Klyachko (1998); Knutson and Tao (1999), cf. Knutson and Tao (2001); Belkale (2006), is cut out by the Horn inequalities, which we will recall below.
In the following, we will show that the trace condition and the Horn condition assume their familiar form in the context of Horn’s problem (e.g., Knutson and Tao (2001)), thereby justifying our terminology. For any of cardinality , let denote the permutation that sends to and to the complement while preserving the order of each block. Then (we identify with a permutation matrix). However, does not satisfy the conditions of corollary 4.2. Instead, we shall write with , which reverses the order of each block. Then we find that satisfies the trace condition if and only if
where . For the left-hand side, we compute
where we have defined , which is a sequence of non-increasing non-negative integers. For the right-hand side, observe that consists of pairs of block matrices of the form , hence is isomorphic to , the vector space of pairs of complex -matrix. We conclude that the trace condition for amounts to
It is not hard to see that (7.1) is precisely equivalent to Horn’s trace condition.
where we have defined etc. A short computation then reveals that the element of proposition 1.2 identifies with the highest weight
of , where denotes the weight of the determinant representation of . On the other hand, by using the isomorphism
the moment cone on gets likewise identified with the direct product of the Horn cones . Thus the Horn condition (proposition 1.2) asserts that as well as . These two conditions are in fact equivalent by a well-known duality of the Littlewood-Richardson coefficients. Therefore, we arrive at a single condition
This condition is not only necessary for to be a facet, but it is also sufficient for to be a valid inequality, as is well-known (e.g., Knutson and Tao (2001)). We will show in the next section how in the context of our work this can be deduced from the saturation conjecture and Schubert calculus (in fact, we shall see that (7.5) implies that is a Ressayre element). In particular, we obtain the familiar recursive definition of Horn’s inequalities: Define to be the set of all triples of cardinality that satisfy the trace condition (7.1) as well as
for all and all triples . Then (7.5) implies that satisfies the Horn condition if and only if . Therefore, our trace and Horn conditions are indeed a generalization of the classical conditions due to Horn.
2 The Determinant Polynomial
In this section, we will show that any that satisfies the trace and Horn condition is automatically a Ressayre element, i.e., that the determinant polynomial is non-zero. To start, we note that by the saturation property of the Littlewood-Richardson coefficients Knutson and Tao (1999) (cf. Belkale (2006)), the Horn condition (7.5) implies that the following space of -invariants is non-zero:
A natural candidate is certainly the determinant polynomial itself (from which we had obtained the Horn condition), but it is not obvious that the Horn condition should imply that . We will give two alternative arguments that show that this is indeed the case.
The first proof follows an argument of Belkale Belkale (2006). We start with the observation that
Using that , we recognize that (7.6) is equivalent to the cohomological condition
a close variant of what we had analyzed in the proof of proposition 3.13. We will use the parametrizations
We show that the fibers of are generically non-empty. Thus let . Then is an element of the fiber if and only if
Any such and is automatically invertible, and therefore a general element in the stabilizer group of . Thus the fiber is non-empty if and only if we can find such that
Since etc., this is the case if and only if
By Kleiman’s transversality theorem, the cohomological condition (7.7) ensures that this is the case for generic . ∎
Now observe that using the identifications (7.2) and (7.4), the tangent map (3.6) at some base point reads
The determinant polynomial is non-zero, i.e., is a Ressayre element.
By Sard’s theorem and lemma 7.1, there exists a point where the differential of is surjective. By writing the differential of (7.9) in coordinates and comparing with (7.10), it is not hard to see that surjectivity of the former at some point implies surjectivity of the tangent map at , where and Thus is a Ressayre element. ∎
Let , , and , so that . Then and , and the associated Littlewood-Richardson coefficients are given by
In table 5 we list a set of homogeneous generators of the algebra of lowest weight vectors in . Note that and span the two-dimensional subspace of weight , , . An explicit calculation best left to a computer algebra system Stein et al. (2014) verifies that
Therefore, in agreement with (7.3) is indeed a lowest weight vector of weight
Note that if we consider and as coefficients rather than indeterminates, spans the two-dimensional subspace of lowest weight vectors of weight as we vary and over . Likewise, spans the subspace of lowest weight vectors of weight if we instead vary and .
where . Equation (7.11) can be shown by manual inspection, relating the matrix elements of the tangent map (3.6) with the matrix constructed by Howe et al. This gives a geometric interpretation of the invariants constructed in Howe et al. (2005) – namely, as the determinant of the tangent map (7.10) associated with a Ressayre element –, and it also serves as an alternative, second proof that is implied by the Horn condition via (7.6).
Acknowledgments. We would like to thank Velleda Baldoni, Soo Teck Lee, Nicolas Ressayre, and Jonathan Skowera for pleasant discussions. We acknowledge the National Technological University and the Institute for Mathematical Sciences at the National University of Singapore for their hospitality during the program on Inverse Moment Problems, where this work had been initiated. MW acknowledges financial support by the Swiss National Science Foundation (grants PP00P2-128455, 20CH21-138799 (CHIST-ERA project CQC)), the Swiss National Center of Competence in Research ‘Quantum Science and Technology (QSIT)’, the Swiss State Secretariat for Education and Research supporting COST action MP1006, the European Research Council under the European Union’s Seventh Framework Programme (FP/2007–2013)/ERC Grant Agreement no. 337603, the Simons Foundation, and FQXi.
Appendix A On the Dimension of the Kronecker Cone
implies that , , and are scalars (i.e., proportional to the identity matrix).
Now suppose that , , and are unitaries such that (A.1) holds. Then,
Let . Then the Kronecker cone is maximal-dimensional.
In the case where the following lemma strengthens corollary A.2:
for any triple of small perturbations . Therefore, the only constraints in the vicinity of are the Weyl chamber inequalities. ∎