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 d×dd\times d 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 t\mathfrak{t} denote the Lie algebra of a maximal torus of TT, it+∗i\mathfrak{t}^{*}_{+} a positive Weyl chamber, and π\pi the (complex) Lie algebra representation induced by Π\Pi (see section 2 below for precise definitions). We shall say that HH is a Ressayre element if 1) the hyperplane (−,H)=0(-,H)=0 is spanned by weights of M{M} and 2) there exists a vector ψ∈M\psi\in{M} annihilated by π(H)\pi(H) such that the “tangent map” at ψ\psi,

is an isomorphism; here, M(H<0){M}(H<0) denotes the direct sum of all negative eigenspaces of π(H)\pi(H) and n−(H<0)\mathfrak{n}_{-}(H<0) the sum of all root spaces for negative roots α\alpha such that (α,H)<0(\alpha,H)<0 (definition 3.14). Note that there are only finitely many Ressayre elements for any given representation π\pi. In section 3, we will prove the following result:

The moment cone for the KK-action on M{M} 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 HH is a Ressayre element then the domain and codomain of the tangent map (1.1) necessarily have the same dimension, i.e., dim⁡n−(H<0)=dim⁡M(H<0)\dim\mathfrak{n}_{-}(H<0)=\dim{M}(H<0). We call this the trace condition. Moreover, note that the determinant δH\delta_{H} of (1.1) (with respect to any fixed pair of bases) is a non-zero polynomial in ψ∈M(H=0)\psi\in{M}(H=0), the zero eigenspace of π(H)\pi(H). In fact, δH\delta_{H} is a canonical (up to scalar multiplication) lowest weight vector for the action of K(H=0)K(H=0), the centralizer of the torus generated by HH, on the space of polynomials on M(H=0){M}(H=0). This implies the following result in section 4, which we call the Horn condition:

is an element of the moment cone CK(H=0)(M(H=0))C_{K(H=0)}({M}(H=0)). In fact, δH\delta_{H} is a lowest weight vector of weight −κH-\kappa_{H}.

Here, Ω\Omega denotes the set of weights of M{M} and RG,−R_{G,-} denotes the set of negative roots of GG. 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 κH\kappa_{H}.

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 SU⁡(d)\operatorname{SU}(d), whose complexification is SL⁡(d)\operatorname{SL}(d), we will always use the maximal torus T(d)T(d) that consists of the diagonal unitary matrices of unit determinant. Its Lie algebra will be denoted by t(d)\mathfrak{t}(d). We will use as positive roots the αi,j(H):=Hii−Hjj\alpha_{i,j}(H):=H_{ii}-H_{jj} with i<ji<j, and abbreviate the (positive) roots by RdR_{d} and Rd,+R_{d,+}, respectively. Finally, we write Oλd{\mathcal{O}}^{d}_{\lambda} for the coadjoint SU⁡(d)\operatorname{SU}(d) orbit through a highest weight λ\lambda.

Now let Π ⁣:G→GL⁡(M)\Pi\colon G\rightarrow\operatorname{GL}({M}) a representation of GG on a finite-dimensional Hilbert space M{M} that is equipped with a KK-invariant Hermitian inner product ⟨−⟩−\braket{-}{-}, which we take to be antilinear in the first argument. We will oftentimes use Dirac’s notation ⟨ϕ⟩Aψ:=⟨ϕ⟩Aψ\braket{\phi}{A}{\psi}:=\braket{\phi}{A\psi} for ϕ,ψ∈M\phi,\psi\in{M} and A∈gl(M)A\in\mathfrak{gl}({M}). We denote by π ⁣:g→gl(M)\pi\colon\mathfrak{g}\rightarrow\mathfrak{gl}({M}) the induced Lie algebra representation, by Ω⊆PG\Omega\subseteq P_{G} the set of weights and write the weight space decomposition as M=⨁ω∈ΩMω{M}=\bigoplus_{\omega\in\Omega}{M}_{\omega}. The KK-action on M{M} admits a canonical (up to conventions) moment map, defined by

for all ψ∈M\psi\in{M} and X∈ikX\in i\mathfrak{k}. Here and in the following, we write (φ,X)=φ(X)(\varphi,X)=\varphi(X) for the duality pairing. The map μK\mu_{K} is indeed a moment map in the sense of symplectic geometry: it is KK-invariant and satisfies the basic identity

for all X∈kX\in\mathfrak{k}, where ωM(ϕ,ψ)=2Im⁡⟨ϕ⟩ψ\omega_{M}(\phi,\psi)=2\operatorname{Im}\braket{\phi}{\psi} denotes the symplectic form that we will use for M{M}. The moment cone then is defined as intersection of the moment map image with the positive Weyl chamber,

where R(M)=Sym⁡(M)∗R({M})=\operatorname{Sym}({M})^{*} denotes the space of polynomials on M{M}. A representation VG,λ∗V_{G,\lambda}^{*} occurs in R(M)R({M}) if and only if there exists a polynomial P∈R(M)P\in R({M}) of weight −λ-\lambda that is invariant under the action of the lower unipotent subgroup N−N_{-}. We shall call such a polynomial a lowest weight vector in R(M)R({M}).

(and CKC_{K} is a pointed cone with base ΔK\Delta_{K}). 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., dim⁡CK=dim⁡it∗\dim C_{K}=\dim i\mathfrak{t}^{*}. 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 (−,H)≥0(-,H)\geq 0. Since we have assumed that CKC_{K} is of maximal dimension, its facets are of codimension one in it∗i\mathfrak{t}^{*} and their inward-pointing normal vectors may be identified with the defining linear inequalities (−,H)≥0(-,H)\geq 0 of the moment cone. Since the moment cone is obtained by intersecting μK(M)\mu_{K}({M}) with the positive Weyl chamber, which itself is a maximal-dimensional polyhedral cone, some of the facets of CKC_{K} can be subsets of facets of it+∗i\mathfrak{t}^{*}_{+}, 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 (−,Hα)≥0(-,H_{\alpha})\geq 0 for some positive root α∈RG,+\alpha\in R_{G,+}. 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 CKC_{K} meets the relative interior it>0∗i\mathfrak{t}^{*}_{>0} of the positive Weyl chamber.

We first consider the moment map μT\mu_{T}, defined as in (2.2) for the action of the maximal torus T⊆KT\subseteq K. Let M=⨁ω∈ΩMω{M}=\bigoplus_{\omega\in\Omega}{M}_{\omega} be the decomposition of M{M} into weight spaces and let ψ∈M\psi\in{M} be a vector decomposed accordingly as ψ=∑ωψωvω\psi=\sum_{\omega}\psi_{\omega}v_{\omega}. Then μT\mu_{T} has the following concrete description:

Observe that μT(ψ)\mu_{T}(\psi) is a conic combination of weights. It follows that the “abelian” moment cone CTC_{T} of M{M} is precisely equal to the conical hull of the set of weights; it is maximal-dimensional since it contains CKC_{K}. More generally, if Ω′⊆Ω\Omega^{\prime}\subseteq\Omega is a subset of weights and MΩ′:=⨁ω∈Ω′Mω{M}_{\Omega^{\prime}}:=\bigoplus_{\omega\in\Omega^{\prime}}{M}_{\omega} then CT(MΩ′)=cone⁡Ω′C_{T}({M}_{\Omega^{\prime}})=\operatorname{cone}\Omega^{\prime}. For the next lemma recall that a critical point of a smooth map f ⁣:M→M′f\colon M\rightarrow M^{\prime} is a point m∈Mm\in M where the differential df\big{|}_{m} is not surjective; a critical value is the image f(m)f(m) 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 μT\mu_{T} is equal to the union of the codimension-one conic hulls of subsets of weights.

Let ψ∈M\psi\in{M} with weight decomposition ψ=∑ωψωvω\psi=\sum_{\omega}\psi_{\omega}v_{\omega}. By (3.1), μT(ψ)\mu_{T}(\psi) is a conic combination of weights in Ωψ:={ω∈Ω:ψω≠0}\Omega_{\psi}:=\{\omega\in\Omega:\psi_{\omega}\neq 0\}. By the moment map property (2.3) and non-degeneracy of the symplectic form, ψ\psi is a critical point if and only if there exists 0≠X∈t0\neq X\in\mathfrak{t} such that π(X)ψ=0\pi(X)\psi=0 (Guillemin and Sternberg, 1982a, Lemma 2.1), i.e., if and only if ω(X)=0\omega(X)=0 for all ω∈Ωψ\omega\in\Omega_{\psi}. It follows that ψ\psi is a critical point if and only if the conic hull of Ωψ\Omega_{\psi} 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 Ω′⊆Ω\Omega^{\prime}\subseteq\Omega is a subset of weights that spans a conic hull of codimension one then CT(MΩ′)=cone⁡Ω′C_{T}({M}_{\Omega^{\prime}})=\operatorname{cone}\Omega^{\prime} 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 H∈itH\in i\mathfrak{t} is called admissible if the linear hyperplane (−,H)=0(-,H)=0 is spanned by a subset of weights in Ω\Omega.

The notion of admissibility is invariant under the action of the Weyl group WKW_{K}.

Let (−,H)≥0(-,H)\geq 0 be an inequality corresponding to a non-trivial facet of the moment cone. Then HH is admissible.

By lemma 3.2, the intersection of (−,H)=0(-,H)=0 with the interior of the positive Weyl chamber it>0∗i\mathfrak{t}^{*}_{>0} is non-empty. Each point in this intersection is a critical value for (μK,H)=(μT,H)(\mu_{K},H)=(\mu_{T},H), hence of μT\mu_{T}, 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 (−,H)≥0(-,H)\geq 0 corresponding to a non-trivial facet (−,H)=0(-,H)=0 of the moment cone. Let ψ∈M\psi\in{M} be a preimage of a point μK(ψ)∈it>0∗\mu_{K}(\psi)\in i\mathfrak{t}^{*}_{>0} on the facet (−,H)=0(-,H)=0. In the proof of lemma 3.5, we have used that ψ\psi is a critical point of (μK,H)(\mu_{K},H) (equivalently, that π(H)ψ=0\pi(H)\psi=0) to gain information on the set of possible facets. To study the function (μK,H)(\mu_{K},H) in the vicinity of such a critical point ψ\psi it is natural to consider the Hessian, which is the quadratic form

For tangent vectors generated by the infinitesimal action of X,Y∈kX,Y\in\mathfrak{k}, we have the formula

where we have used that π(H)ψ=0\pi(H)\psi=0. We now decompose

where M(H<0)=⨁ω:(ω,H)<0Mω{M}(H<0)=\bigoplus_{\omega:(\omega,H)<0}{M}_{\omega} is the sum of the eigenspaces of the Hermitian operator π(H)\pi(H) with eigenvalue less than , etc. Then it is plain from (3.2) that the index of the Hessian QQ, i.e., the dimension of a maximal subspace on which the quadratic form is negative definite, is equal to twice the complex dimension of M(H<0){M}(H<0).

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 VV that are mapped to it∗i\mathfrak{t}^{*} by the differential of the moment map. To see this, consider a curve ψt\psi_{t} with ψ0=ψ\psi_{0}=\psi, ψ˙0=V\dot{\psi}_{0}=V and μK(ψt)∈it>0∗\mu_{K}(\psi_{t})\in i\mathfrak{t}^{*}_{>0} for all t∈(−ε,ε)t\in(-\varepsilon,\varepsilon) (such a curve can always be constructed by using the symplectic cross section (Guillemin and Sternberg, 1984b, Theorem 26.7)). Then, since (μK(ψ),H)=0(\mu_{K}(\psi),H)=0 and d(\mu_{K},H)\big{|}_{\psi}\equiv 0,

But (μK(ψt),H)≥0(\mu_{K}(\psi_{t}),H)\geq 0, since μK(ψt)∈CK\mu_{K}(\psi_{t})\in C_{K} and (−,H)≥0(-,H)\geq 0 is a valid inequality for the moment cone. Together, this shows that, indeed, Q(V,V)≥0Q(V,V)\geq 0.

The subspace of all such VV can be computed in a different way. For this, let r=⨁α∈RG,+kα\mathfrak{r}=\bigoplus_{\alpha\in R_{G,+}}\mathfrak{k}_{\alpha} 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 M0:=π(r)ψ⊆MM_{0}:=\pi(\mathfrak{r})\psi\subseteq{M}.

In fact, it is well-known that M0M_{0} is a symplectic subspace of M{M} (see, e.g., (Guillemin and Sternberg, 1982a, Lemma 6.7)). To see this, note that the restriction of the symplectic form to M0M_{0} is given by

and can therefore be identified with the Kirillov-Kostant-Souriau symplectic form on the coadjoint orbit through μK(ψ)∈it>0∗\mu_{K}(\psi)\in i\mathfrak{t}^{*}_{>0}. As a consequence, we have the decomposition M=M0⊕M0ωM{M}=M_{0}\oplus M_{0}^{\omega_{M}}.

The Hessian is block-diagonal with respect to the decomposition M=M0⊕M0ωM{M}=M_{0}\oplus M_{0}^{\omega_{M}}.

For all π(R)ψ∈M0\pi(R)\psi\in M_{0} and V∈M0ωMV\in M_{0}^{\omega_{M}} we have that

since [−iH,R]∈r[-iH,R]\in\mathfrak{r} and therefore π([−iH,R])ψ∈M0\pi([-iH,R])\psi\in M_{0}. ∎

The tangent map r→M\mathfrak{r}\rightarrow{M}, R↦π(R)ψR\mapsto\pi(R)\psi is injective.

The stabilizer of the coadjoint action of KK at any λ∈it>0∗\lambda\in i\mathfrak{t}^{*}_{>0} is TT, while k=t⊕r\mathfrak{k}=\mathfrak{t}\oplus\mathfrak{r}. As μK(ψ)∈it>0∗\mu_{K}(\psi)\in i\mathfrak{t}^{*}_{>0}, the claim follows from this and KK-equivariance of the moment map. ∎

The index of the Hessian is equal to twice the number of positive roots α∈RG,+\alpha\in R_{G,+} such that (α,H)>0(\alpha,H)>0.

Since QQ is positive semidefinite on M0ωMM_{0}^{\omega_{M}} (lemma 3.6) and block-diagonal with respect to the decomposition M=M0⊕M0ωM{M}=M_{0}\oplus M_{0}^{\omega_{M}} (lemma 3.7), it suffices to compute the index of QQ on M0={π(R)ψ:R∈r}M_{0}=\{\pi(R)\psi:R\in\mathfrak{r}\}. For this, recall from (3.3) that Q(π(R)ψ,π(S)ψ)=(μK(ψ),[[H,R],S])Q(\pi(R)\psi,\pi(S)\psi)=(\mu_{K}(\psi),[[H,R],S]) for all R,S∈rR,S\in\mathfrak{r}. Since the tangent map R↦π(R)ψR\mapsto\pi(R)\psi is injective (lemma 3.8), we may instead consider the form

on r\mathfrak{r}. Now observe that Q~\widetilde{Q} is block-diagonal with respect to r=⨁α∈RG,+kα\mathfrak{r}=\bigoplus_{\alpha\in R_{G,+}}\mathfrak{k}_{\alpha}, since for all R∈kαR\in\mathfrak{k}_{\alpha} and S∈kβS\in\mathfrak{k}_{\beta}, [[H,R],S]∈ikα±β[[H,R],S]\in i\mathfrak{k}_{\alpha\pm\beta}, while μK(ψ)∈it∗\mu_{K}(\psi)\in i\mathfrak{t}^{*}. Therefore, it suffices to compute the index on a single root space kα\mathfrak{k}_{\alpha}. For this, define the “Pauli matrices” Xα:=Eα+E−αX_{\alpha}:=E_{\alpha}+E_{-\alpha} and Yα:=i(E−α−Eα)Y_{\alpha}:=i(E_{-\alpha}-E_{\alpha}), which satisfy the commutation relations [Xα,Yα]=2iHα[X_{\alpha},Y_{\alpha}]=2iH_{\alpha} etc. Then iXαiX_{\alpha} and iYαiY_{\alpha} form a basis of kα\mathfrak{k}_{\alpha} and

Likewise, Q~(iYα,iYα)=−2(α,H)(μK(ψ),Hα)\widetilde{Q}(iY_{\alpha},iY_{\alpha})=-2(\alpha,H)(\mu_{K}(\psi),H_{\alpha}), while Q~(iXα,iYα)=0\widetilde{Q}(iX_{\alpha},iY_{\alpha})=0. Since (μK(ψ),Hα)>0(\mu_{K}(\psi),H_{\alpha})>0, we conclude that the index of QQ is equal to twice the number of positive roots α\alpha with (α,H)>0(\alpha,H)>0. ∎

Since we have already seen above that the index of the Hessian is also equal to twice the complex dimension of M(H<0){M}(H<0), we obtain the following result, which can also be extracted from (Brion, 1999, Theorem 2):

If (−,H)≥0(-,H)\geq 0 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 n−=⨁α∈RG,−gα\mathfrak{n}_{-}=\bigoplus_{\alpha\in R_{G,-}}\mathfrak{g}_{\alpha} of the negative unipotent subgroup, which plays a role analogous to r\mathfrak{r} for vectors ψ\psi that are mapped into the positive Weyl chamber (compare the following with lemma 3.8).

The tangent map n−→M,X→π(X)ψ\mathfrak{n}_{-}\rightarrow{M},X\to\pi(X)\psi is injective.

Let E−:=∑α∈RG,+zαE−αE_{-}:=\sum_{\alpha\in R_{G,+}}z_{\alpha}E_{-\alpha} be an arbitrary element in n−\mathfrak{n}_{-}. Since π(E±α)†=π(E∓α)\pi(E_{\pm\alpha})^{\dagger}=\pi(E_{\mp\alpha}), we find that π(E−)†=π(E+)\pi(E_{-})^{\dagger}=\pi(E_{+}) with E+:=∑α∈RG,+zˉαEαE_{+}:=\sum_{\alpha\in R_{G,+}}\bar{z}_{\alpha}E_{\alpha} (the Cartan involution of E−E_{-}). Therefore,

so that by using μK(ψ)∈it>0∗\mu_{K}(\psi)\in i\mathfrak{t}^{*}_{>0} we find that

In contrast to lemma 3.8, which continues to hold true if μK(ψ)\mu_{K}(\psi) 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 GG-representation M=VG,λ{M}=V_{G,\lambda} with highest weight λ∈it>0∗\lambda\in i\mathfrak{t}^{*}_{>0} and highest weight vector vλv_{\lambda}. Then μK(vλ)=λ∈it>0∗\mu_{K}(v_{\lambda})=\lambda\in i\mathfrak{t}^{*}_{>0} and the “lowering operators” in n−\mathfrak{n}_{-} indeed act injectively. On the other hand, the “raising operators” in the positive nilpotent Lie algebra n+\mathfrak{n}_{+} annihilate the highest weight vector (by definition).

We now decompose the Lie algebra n−\mathfrak{n}_{-} similarly to (3.4),

where n−(H<0)=⨁α∈RG,−:(α,H)<0gα\mathfrak{n}_{-}(H<0)=\bigoplus_{\alpha\in R_{G,-}:(\alpha,H)<0}\mathfrak{g}_{\alpha} is the sum of the complex root spaces with negative HH-weight (α,H)<0(\alpha,H)<0, etc. We observe that corollary 3.10 can be equivalently stated as

Note that π(n−(H<0))M(H=0)⊆M(H<0)\pi(\mathfrak{n}_{-}(H<0)){M}(H=0)\subseteq{M}(H<0). Thus we obtain the following important result:

Let ψ∈M\psi\in{M} such that μK(ψ)∈it>0∗\mu_{K}(\psi)\in i\mathfrak{t}^{*}_{>0} is a point on a non-trivial facet of the moment cone corresponding to the inequality (−,H)≥0(-,H)\geq 0. Then ψ∈M(H=0)\psi\in{M}(H=0) and the tangent map restricts to an isomorphism

The fact that ψ∈M(H=0)\psi\in{M}(H=0) is just a reformulation of π(H)ψ=0\pi(H)\psi=0. By the preceding discussion, the tangent map is well-defined as a map from n−(H<0)\mathfrak{n}_{-}(H<0) to M(H<0){M}(H<0); 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 ψ∈M(H=0)\psi\in{M}(H=0) such that the tangent map

is surjective. Then (−,H)≥0(-,H)\geq 0 is a valid inequality for the moment cone.

Its differential at (1,ψ)(1,\psi) is the linear map

The assumption implies that this map is surjective. It follows that Π(N−)M(H≥0)⊆M\Pi(N_{-}){M}(H\geq 0)\subseteq{M} contains a small Euclidean ball around ψ\psi. In particular, any N−N_{-}-invariant polynomial that is zero on M(H≥0){M}(H\geq 0) is automatically zero everywhere on M{M}.

We now prove the inequality. By the description of the moment cone in (2.4), it suffices to show that (λ,H)≥0(\lambda,H)\geq 0 for all highest weights λ\lambda such that VG,λ∗⊆R(M)V^{*}_{G,\lambda}\subseteq R({M}). Recall that the highest weight of VG,λ∗V^{*}_{G,\lambda} is λ∗=−w0λ\lambda^{*}=-w_{0}\lambda, where w0w_{0} is the longest Weyl group element that flips the positive and negative roots. Consider a lowest weight vector, i.e., a polynomial P∈R(M)P\in R({M}) that is a weight vector of weight −λ-\lambda and invariant under the action of N−N_{-}. Then π(H)P=−(λ,H)P\pi(H)P=-(\lambda,H)P and the restriction of PP to M(H≥0){M}(H\geq 0) is non-zero by our discussion above. But this restriction is an element of R(M(H≥0))=Sym⁡(M(H≥0))∗R({M}(H\geq 0))=\operatorname{Sym}({M}(H\geq 0))^{*}, the space of polynomials on M(H≥0){M}(H\geq 0). Since all HH-weights in R(M(H≥0))R({M}(H\geq 0)) are non-positive, it follows that −(λ,H)≤0-(\lambda,H)\leq 0, as we set out to prove. ∎

We remark that proposition 3.13 holds unconditionally without any assumption on the dimension of the moment cone CKC_{K}. We summarize our findings in the following definition and theorem that we had already advertised in the introduction.

An element H∈itH\in i\mathfrak{t} is called a Ressayre element if

HH is admissible, i.e., the linear hyperplane (−,H)=0(-,H)=0 is spanned by a subset of weights in Ω\Omega, and

there exists ψ∈M(H=0)\psi\in{M}(H=0) 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 KK-representation M{M} (under the assumption that CKC_{K} 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 N−×N−(H≥0)M(H≥0)→MN_{-}\times_{N_{-}(H\geq 0)}{M}(H\geq 0)\rightarrow{M} 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 HH, 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 H0H_{0} and then determine those H∈WK⋅H0H\in W_{K}\cdot H_{0} which satisfy the trace condition (4.1). For this, we will need the following well-known lemma.

Fix any ρK∈it+\rho_{K}\in i\mathfrak{t}_{+} such that (α,ρK)>0(\alpha,\rho_{K})>0 for all positive roots α∈RG,+\alpha\in R_{G,+}, and set Hε:=H−ερKH^{\varepsilon}:=H-\varepsilon\rho_{K}. Let w∈WKw\in W_{K} be a Weyl group element such that H=w⋅H0H=w\cdot H_{0}. For all positive roots α∈RG,+\alpha\in R_{G,+} and ε>0\varepsilon>0 small enough,

if and only if (α,H0)=0(\alpha,H_{0})=0 implies that w⋅α∈RG,−w\cdot\alpha\in R_{G,-}. In other words, H0ε:=w−1⋅Hε∈it+H^{\varepsilon}_{0}:=w^{-1}\cdot H^{\varepsilon}\in i\mathfrak{t}_{+} if and only if w0w⋅α∈RG,+w_{0}w\cdot\alpha\in R_{G,+} for all α∈RG,+\alpha\in R_{G,+} with (α,H0)=0(\alpha,H_{0})=0. Since HεH^{\varepsilon} is regular, this immediately shows that such Weyl group elements ww 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 HH 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 G(H=0)G(H=0) denote the identity component of the centralizer of the torus generated by HH, i.e., the connected subgroup of GG with Lie algebra g(H=0)=h⊕⨁α:(α,H)=0gα\mathfrak{g}(H=0)=\mathfrak{h}\oplus\bigoplus_{\alpha:(\alpha,H)=0}\mathfrak{g}_{\alpha}; denote by N−(H=0)N_{-}(H=0) the corresponding negative unipotent subgroup and by K(H=0)=G(H=0)∩KK(H=0)=G(H=0)\cap K the maximal compact subgroup. Then G(H=0)G(H=0) acts on M(H=0){M}(H=0) and thus on the space R(M(H=0))=Sym⁡(M(H=0))∗R({M}(H=0))=\operatorname{Sym}({M}(H=0))^{*} of polynomial functions.

is an element of the moment cone CK(H=0)(M(H=0))C_{K(H=0)}({M}(H=0)). In fact, δH\delta_{H} is a lowest weight vector of weight −κH-\kappa_{H}.

In view of (2.4), it suffices to argue that the determinant polynomial δH\delta_{H} is a lowest weight vector in R(M(H=0))R({M}(H=0)) of weight −κH-\kappa_{H}.

To see this, fix a basis ψ1,…,ψk\psi_{1},\dots,\psi_{k} of M(H<0){M}(H<0) such that each ψj\psi_{j} is a weight vector of weight ωj\omega_{j}, and denote by α1,…,αk\alpha_{1},\dots,\alpha_{k} the negative roots with (α,H)<0(\alpha,H)<0, so that Ek:=EαkE_{k}:=E_{\alpha_{k}} is a basis of n−(H<0)\mathfrak{n}_{-}(H<0). Then the determinant polynomial δH\delta_{H} with respect to this basis can be written as

where we write ΛkA\Lambda^{k}A for the canonical homomorphism ΛkV→ΛkW\Lambda^{k}V\rightarrow\Lambda^{k}W induced by a linear map A ⁣:V→WA\colon V\rightarrow W. It follows that for any g∈Gg\in G and ψ∈M(H=0)\psi\in{M}(H=0),

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 (−,H)≥0(-,H)\geq 0, 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 (rK−1)(r_{K}-1)-element subsets of Ω\Omega 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 δH\delta_{H} 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 ψ\psi in a Hilbert space M{M}, 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, M=⨂k=1nMk{M}=\bigotimes_{k=1}^{n}{M}_{k}. The state of the kk-th subsystem can be described by the reduced density matrix ρk\rho_{k}, which is the unique positive semi-definite operator on Mk{M}_{k} such that

for all Hermitian operators XkX_{k} on Mk{M}_{k} (it is also called the partial trace of ψ\psi). The fundamental one-body quantum marginal problem asks which ρ1,…,ρn\rho_{1},\dots,\rho_{n} can arise as the reduced density matrices of a quantum state ψ∈M\psi\in{M} 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 ρk\rho_{k} have to satisfy in order for there to exist a global state ψ\psi. For two subsystems, it is a straightforward consequence of the singular value decomposition that ρ1\rho_{1} and ρ2\rho_{2} 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 C(a,b,c)C(a,b,c) 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 α\alpha, β\beta, and γ\gamma vary over the set of Young diagrams with the same number kk of boxes and no more than aa, bb, and cc rows, respectively, and where gα,β,γg_{\alpha,\beta,\gamma} denotes the Kronecker coefficient of the symmetric group, i.e., the multiplicity of the invariant subspace in the corresponding triple tensor product of irreducible Sk−S_{k}-representations. We will henceforth refer to C(a,b,c)C(a,b,c) as the Kronecker cone.

It will be convenient to assume without loss of generality that 1<a≤b≤c≤ab1<a\leq b\leq c\leq ab. In this case, C(a,b,c)C(a,b,c) 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 c=abc=ab. According to theorem 1.1, the moment cone C(a,b,c)C(a,b,c) is thus cut out by those HH which are Ressayre elements, i.e., which are admissible and satisfy the isomorphism condition. Naively determining the admissible HH by enumerating all subsets of Ω\Omega with cardinality r−1r-1, determining whether they span a hyperplane and computing the normal vector amounts to considering (∣Ω∣r−1)=(abca+b+c−3)\binom{\lvert\Omega\rvert}{r-1}=\binom{abc}{a+b+c-3} subsets, which rapidly becomes infeasible (e.g., for a=b=c=4a=b=c=4, there are over 27 billion such subsets). We therefore need to derive additional constraints to make this approach computationally feasible.

By considering SU⁡(a)×SU⁡(b)×SU⁡(c)⊆SU⁡(ab)×SU⁡(c)\operatorname{SU}(a)\times\operatorname{SU}(b)\times\operatorname{SU}(c)\subseteq\operatorname{SU}(ab)\times\operatorname{SU}(c), the moment cone for the tripartite problem can now be written in terms of the bipartite cone C(ab,c)C(ab,c) and the moment polytopes Δ(a,b∣λAB)\Delta(a,b|\lambda_{AB}) of the action of SU⁡(a)×SU⁡(b)\operatorname{SU}(a)\times\operatorname{SU}(b) on the coadjoint SU⁡(ab)\operatorname{SU}(ab)-orbits OλABab{\mathcal{O}}^{ab}_{\lambda_{AB}}:

A first consequence is that the case c≠abc\neq ab can always be reduced to c=abc=ab. Indeed, (6.2) and (6.3) imply that

if c<abc<ab. Thus the moment cone for c>abc>ab is isometric to C(a,b,ab)C(a,b,ab), while for c<abc<ab it is obtained as a projection of the latter. In the case where c=abc=ab, it was observed in Manivel (1997); Klyachko (2004) that any normal vector of the moment cone C(a,b,ab)C(a,b,ab) 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 H=(HA,HB,HAB,z)H=(H_{A},H_{B},H_{AB},z) be the normal vector of a non-trivial facet of the moment cone C(a,b,ab)C(a,b,ab) such that (HA,HB)≠0(H_{A},H_{B})\neq 0. Then:

(HA,HB)(H_{A},H_{B}) is determined by a maximal number of equations of the form HA,i+HB,j=HA,k+HB,lH_{A,i}+H_{B,j}=H_{A,k}+H_{B,l},

the components of HABH_{AB} are precisely all possible partial sums −HA,i−HB,j-H_{A,i}-H_{B,j}, and

By lemma 3.2, there exists a regular dominant point λ=(λA,λB,λAB)∈it>0∗\lambda=(\lambda_{A},\lambda_{B},\lambda_{AB})\in i\mathfrak{t}^{*}_{>0} in the interior of this facet. In a neighborhood of λ\lambda, the moment cone locally looks like a half-space, so that

is not only a valid inequality that holds for all (μA,μB)∈Δ(a,b∣λAB)(\mu_{A},\mu_{B})\in\Delta(a,b|\lambda_{AB}), as follows from (6.3), but in fact a facet of the moment polytope Δ(a,b∣λAB)\Delta(a,b|\lambda_{AB}), since we have assumed that (HA,HB)≠0(H_{A},H_{B})\neq 0.

While determining the moment polytopes Δ(a,b∣λAB)\Delta(a,b|\lambda_{AB}) is just as hard as determining the cone C(a,b,ab)C(a,b,ab), this reformulation gives us an additional insight: Recall that the Duistermaat-Heckman measure for the action of the maximal torus T(a)×T(b)T(a)\times T(b) of SU⁡(a)×SU⁡(b)\operatorname{SU}(a)\times\operatorname{SU}(b) on the coadjoint SU⁡(ab)\operatorname{SU}(ab)-orbit OλABab{\mathcal{O}}^{ab}_{\lambda_{AB}} is a measure on the abelian moment polytope. Now let π\pi denote the projection it(ab)∗→it(a)∗⊕it(b)∗i\mathfrak{t}(ab)^{*}\rightarrow i\mathfrak{t}(a)^{*}\oplus i\mathfrak{t}(b)^{*}. The affine hyperplanes through some π(wλAB)\pi(w\lambda_{AB}) spanned by subsets of the restricted roots π(α)\pi(\alpha), α∈Rab,+\alpha\in R_{ab,+}, 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 SU⁡(a)×SU⁡(b)\operatorname{SU}(a)\times\operatorname{SU}(b) can be recovered by applying a number of partial derivatives to the measure for T(a)×T(b)T(a)\times T(b) (e.g., Christandl et al. (2014)). It follows that the moment polytope Δ(a,b∣λAB)\Delta(a,b|\lambda_{AB}), 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 (HA,HB)(H_{A},H_{B}) is defined by a maximal number of equations of the form

for some indices i,j,k,li,j,k,l. This shows the first assertion. Moreover, since the facet contains π(wλAB)\pi(w\lambda_{AB}) for some w∈Sabw\in S_{ab} we obtain that

Following Klyachko, we shall call any (HA,HB)(H_{A},H_{B}) 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 E+(a,b)\mathcal{E}_{+}(a,b). The extremal edges span the extreme rays of the cubicles, which are the full-dimensional convex cones of elements (HA,HB)(H_{A},H_{B}) cut out by a maximal set of inequalities of the form HA,i+HB,j≥HA,k+HB,lH_{A,i}+H_{B,j}\geq H_{A,k}+H_{B,l}. In other words, a cubicle is defined as a set of (HA,HB)(H_{A},H_{B}) with fixed order of the HA,i+HB,jH_{A,i}+H_{B,j} and can therefore be encoded by a standard Young tableaux of rectangular shape a×ba\times b 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 C(a,b,ab)C(a,b,ab) incident to this face Manivel (1997); Klyachko (2004). However, not all non-trivial facets can be obtained in this way.

Let H=(HA,HB,HC,z)H=(H_{A},H_{B},H_{C},z) be the normal vector of a non-trivial facet of the moment cone C(a,b,c)C(a,b,c), where c≤abc\leq ab. Then (HA,HB)=0(H_{A},H_{B})=0 or (HA,HB)(H_{A},H_{B}) is proportional to an element in the Sa×SbS_{a}\times S_{b}-orbit of E+(a,b)\mathcal{E}_{+}(a,b) (i.e., an extremal edge up to Sa×SbS_{a}\times S_{b} and rescaling).

In view of (6.4), we may obtain finite and complete set of inequalities for C(a,b,c)C(a,b,c) by taking any facet H=(HA,HB,HAB,0)H=(H_{A},H_{B},H_{AB},0) for C(a,b,ab)C(a,b,ab) and restricting HABH_{AB} to its first cc components (that is, set HCH_{C} to be the traceless part of (HAB,1,…,HAB,c)(H_{AB,1},\dots,H_{AB,c}) and define zz accordingly). Such an inequality will not necessarily define a facet of C(a,b,c)C(a,b,c), 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 (HA,HB)=0(H_{A},H_{B})=0. ∎

Lemma 6.1 is efficient in finding candidates for facets of C(a,b,ab)C(a,b,ab), but not necessarily so for C(a,b,c)C(a,b,c)with c<abc<ab. 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 cc components out of the abab components of the given HABH_{AB}. If c≪abc\ll ab then it can be more efficient to apply corollary 6.2 to all three of (HA,HB)(H_{A},H_{B}), (HA,HC)(H_{A},H_{C}) and (HB,HC)(H_{B},H_{C}) (e.g., in the case of a=b=c=4a=b=c=4). For this we will use the following lemma:

Let (HA,HB)(H_{A},H_{B}) be an extremal edge. Then HAH_{A} and HBH_{B} are each either zero or primitive in the dual of the root lattice.

Case 2: m=it(b)\mathfrak{m}=i\mathfrak{t}(b). Since the matrix with columns the roots of Ab−1A_{b-1} is totally unimodular (e.g., (Schrijver, 1986, p. 274, (18))), it follows that SBS_{B} spans the root lattice. To find a contradiction, suppose that HAH_{A} is neither zero nor primitive. Then we can write HA=nHA′H_{A}=nH^{\prime}_{A} for some n>1n>1 and a non-zero element HA′H^{\prime}_{A} in the dual of the root lattice. But then,

for all (α,β)∈S(\alpha,\beta)\in S. Since the SBS_{B} spans the root lattice, we conclude that also HB/nH_{B}/n is an element of the dual of the root lattice. But then H=(HA,HB)H=(H_{A},H_{B}) is not primitive, which is the desired contradiction. ∎

Let H=(HA,HB,HC,z)H=(H_{A},H_{B},H_{C},z) be the normal vector of a non-trivial facet of the moment cone C(a,b,c)C(a,b,c), where c≤abc\leq ab. Then (HA,HB,HC)(H_{A},H_{B},H_{C}) is proportional to an element in the Sa×Sb×ScS_{a}\times S_{b}\times S_{c}-orbit of

Since the following argument is equivariant under the Weyl group Sa×Sb×ScS_{a}\times S_{b}\times S_{c}, we may without loss of generality assume that (HA,HB,HC)(H_{A},H_{B},H_{C}) is dominant

Case 1: Only on component is non-zero, say, HC≠0H_{C}\neq 0 and therefore (HA,HB)=0(H_{A},H_{B})=0. We may rescale HH such that HCH_{C} is primitive. Then (HA,HC)(H_{A},H_{C}) and (HB,HC)(H_{B},H_{C}) are also primitive and therefore extremal edges by corollary 6.2.

Case 2: At least two components are non-zero, say, HAH_{A} and HBH_{B}. We may rescale HH such that (HA,HB)(H_{A},H_{B}) is primitive. Since HAH_{A} and HBH_{B} are non-zero, lemma 6.3 shows that both HAH_{A} and HBH_{B} are individually primitive. It follows that three of (HA,HB)(H_{A},H_{B}), (HA,HC)(H_{A},H_{C}) and (HB,HC)(H_{B},H_{C}) 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 HH be the normal vector of a non-trivial facet of the moment cone C(a,b,c)C(a,b,c), where c≤abc\leq ab. Then HH is proportional to an element in the Sa×Sb×ScS_{a}\times S_{b}\times S_{c}-orbit of

For any given (HA,HB,HC)∈E+(a,b,c)(H_{A},H_{B},H_{C})\in\mathcal{E}_{+}(a,b,c), there are in general several zz such that H=(HA,HB,HC,z)H=(H_{A},H_{B},H_{C},z) is admissible (or no such zz at all). The sets E+(a,b,c)\mathcal{E}_{+}(a,b,c) and E+,adm\mathcal{E}_{+,\text{adm}} consist of dominant, primitive elements, and they can be easily obtained algorithmically from the sets of extremal edges.

Let HH be the normal vector of a non-trivial facet of the moment cone C(a,b,c)C(a,b,c), where c≤abc\leq ab. Then HH is proportional to an element in

where W(H0)W(H_{0}) denotes the set of triples w=(wA,wB,wC)∈Sa×Sb×Scw=(w_{A},w_{B},w_{C})\in S_{a}\times S_{b}\times S_{c} such that (w0wA,w0wB,w0wC)(w_{0}w_{A},w_{0}w_{B},w_{0}w_{C}) is a triple of shuffles with respect to the components of H0H_{0} whose lengths sum up to dim⁡M(H0<0)\dim{M}(H_{0}<0).

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 C(a,b,c)C(a,b,c) is cut out by those candidates in E(a,b,c)\mathcal{E}(a,b,c) 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 C(d,d,d)C(d,d,d) for any d>1d>1. It is not hard to verify that all pairs formed from (1,…,1,1−d)(1,\dots,1,1-d) and (d−1,−1,…,−1)(d-1,-1,\dots,-1) are extremal edges (formally, we should divide by dd to work with primitive elements, but we refrain from doing so in the interest of readability). Therefore,

There are two possible values of zz which extend the above to an admissible element H0∈E+,adm(d,d,d)H_{0}\in\mathcal{E}_{+,\text{adm}}(d,d,d). The first option is z=−1z=-1. However, the dimension of M(H0<0){M}(H_{0}<0) is 2(d−1)2+d2(d-1)^{2}+d, and therefore strictly larger than the dimension of n−\mathfrak{n}_{-}. Thus the trace condition (4.1) can never be satisfied!

We now verify that HH is a Ressayre element. For this, observe that we have

for any ψ=∑ijkψijkeijk∈M(H=0)\psi=\sum_{ijk}\psi_{ijk}e_{ijk}\in{M}(H=0). It follows at once from theorem 1.1 that (H,λ)≥0(H,\lambda)\geq 0 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 AA, BB, and CC are also valid. We remark that the moment cone CK(H=0)(M(H=0))C_{K(H=0)}({M}(H=0)) is closely related to the Horn cone for U⁡(d−1)\operatorname{U}(d-1) discussed in section 7 below Christandl et al. (2012); Walter (2014); Littlewood (1958); Murnaghan (1955).

2 Computational Results

To verify that a given HH is a Ressayre element, we have implemented a computer program that works for arbitrary representations. To compute the set of candidates E(a,b,c)\mathcal{E}(a,b,c) 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 a=b=ca=b=c, which corresponds to three quantum particles with the same number of degrees of freedom. While the moment cones C(2,2,2)C(2,2,2) and C(3,3,3)C(3,3,3) had already been computed in Higuchi et al. (2003); Bravyi (2004); Franz (2002); Higuchi (2003) using different methods, the cone C(4,4,4)C(4,4,4) had been out of reach using current methods. In contrast, our method allows us the computation of C(4,4,4)C(4,4,4) in a few minutes, since it does not rely on an intermediate computation of the higher-dimensional cone C(4,4,16)C(4,4,16).

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 G=GL⁡(d)×GL⁡(d)×GL⁡(d)G=\operatorname{GL}(d)\times\operatorname{GL}(d)\times\operatorname{GL}(d) on M=gl(d)⊕gl(d){M}=\mathfrak{gl}(d)\oplus\mathfrak{gl}(d) given by

Any non-negative Hermitian matrix can be written in the form a†aa^{\dagger}a; since the spectra of aa†aa^{\dagger} and a†aa^{\dagger}a are equal, the moment cone is equal to

where spec⁡X\operatorname{spec}X denotes the eigenvalues of a Hermitian matrix XX, ordered non-increasingly. We will call C(d)C(d) the Horn cone in dd dimensions. As proved in Klyachko (1998); Knutson and Tao (1999), cf. Knutson and Tao (2001); Belkale (2006), C(d)C(d) 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 I={i1<⋯<ir}⊆[d]I=\{i_{1}<\dots<i_{r}\}\subseteq[d] of cardinality rr, let wIw_{I} denote the permutation that sends [r][r] to II and [r]c[r]^{c} to the complement Ic={i1c<⋯<id−rc}I^{c}=\{i^{c}_{1}<\dots<i^{c}_{d-r}\} while preserving the order of each block. Then EI=wI⋅Er=wIErwI−1E_{I}=w_{I}\cdot E_{r}=w_{I}E_{r}w_{I}^{-1} (we identify wIw_{I} with a permutation matrix). However, wIw_{I} does not satisfy the conditions of corollary 4.2. Instead, we shall write EI=w~I⋅ErE_{I}=\widetilde{w}_{I}\cdot E_{r} with w~I:=wIcw0\widetilde{w}_{I}:=w_{I^{c}}w_{0}, which reverses the order of each block. Then we find that HIJKH_{IJK} satisfies the trace condition if and only if

where Hr:=(Er,Er,Er)H_{r}:=(E_{r},E_{r},E_{r}). For the left-hand side, we compute

where we have defined λI:=(d−r−(ia−a))a=1r\lambda_{I}:=(d-r-(i_{a}-a))_{a=1}^{r}, which is a sequence of non-increasing non-negative integers. For the right-hand side, observe that M(Hr<0){M}(H_{r}<0) consists of pairs of block matrices of the form (00∗0)\left(\begin{smallmatrix}0&0\\ *&0\end{smallmatrix}\right), hence is isomorphic to Md−r,r2M_{d-r,r}^{2}, the vector space of pairs of complex (d−r)×r(d-r)\times r-matrix. We conclude that the trace condition for HIJKH_{IJK} amounts to

It is not hard to see that (7.1) is precisely equivalent to Horn’s trace condition.

where we have defined TI:={X∈Md−r,r:⟨b⟩Xa=0 if ia>ibc}T_{I}:=\{X\in M_{d-r,r}:\braket{b}{X}{a}=0\text{ if }i_{a}>i^{c}_{b}\} etc. A short computation then reveals that the element κHIJK\kappa_{H_{IJK}} of proposition 1.2 identifies with the highest weight

of U⁡(r)3×U⁡(d−r)3\operatorname{U}(r)^{3}\times\operatorname{U}(d-r)^{3}, where χr\chi_{r} denotes the weight of the determinant representation of U⁡(r)\operatorname{U}(r). On the other hand, by using the isomorphism

the moment cone K(HIJK=0)K(H_{IJK}=0) on M(HIJK=0){M}(H_{IJK}=0) gets likewise identified with the direct product of the Horn cones C(r)×C(d−r)C(r)\times C(d-r). Thus the Horn condition (proposition 1.2) asserts that (λI,λJ,λK−2(d−r)χr)∈C(r)(\lambda_{I},\lambda_{J},\lambda_{K}-2(d-r)\chi_{r})\in C(r) as well as (λIc,λJc,λKc−rχd−r)∈C(d−r)(\lambda_{I^{c}},\lambda_{J^{c}},\lambda_{K^{c}}-r\chi_{d-r})\in C(d-r). 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 HIJKH_{IJK} to be a facet, but it is also sufficient for HIJKH_{IJK} 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 HIJKH_{IJK} is a Ressayre element). In particular, we obtain the familiar recursive definition of Horn’s inequalities: Define Horn⁡(d,r)\operatorname{Horn}(d,r) to be the set of all triples I,J,K⊆[d]I,J,K\subseteq[d] of cardinality r<dr<d that satisfy the trace condition (7.1) as well as

for all s<rs<r and all triples (A,B,C)∈Horn⁡(r,s)(A,B,C)\in\operatorname{Horn}(r,s). Then (7.5) implies that HIJKH_{IJK} satisfies the Horn condition if and only if (I,J,K)∈Horn⁡(d,r)(I,J,K)\in\operatorname{Horn}(d,r). 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 HIJKH_{IJK} that satisfies the trace and Horn condition is automatically a Ressayre element, i.e., that the determinant polynomial δIJK:=δHIJK\delta_{IJK}:=\delta_{H_{IJK}} 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 GL⁡(r)\operatorname{GL}(r)-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 δIJK≠0\delta_{IJK}\neq 0. 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 Ω(wI−)=Ω(w0wI)\Omega(w_{I_{-}})=\Omega(w_{0}w_{I}), 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 RIJKR_{IJK} are generically non-empty. Thus let (g,h)∈GL⁡(d)2⊆gl(d)2(g,h)\in\operatorname{GL}(d)^{2}\subseteq\mathfrak{gl}(d)^{2}. Then (x,y,z,p,q)(x,y,z,p,q) is an element of the fiber RIJK−1(g,h)R_{IJK}^{-1}(g,h) if and only if

Any such pp and qq is automatically invertible, and therefore a general element in the stabilizer group of Vr∈Gr⁡(r,d)V_{r}\in\operatorname{Gr}(r,d). Thus the fiber is non-empty if and only if we can find (x,y,z)(x,y,z) such that

Since {x~⋅Vr:x∈TI}=(w0wI)−1Ω(w0wI)\{\widetilde{x}\cdot V_{r}:x\in T_{I}\}=(w_{0}w_{I})^{-1}\Omega(w_{0}w_{I}) 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 (g,h)(g,h). ∎

Now observe that using the identifications (7.2) and (7.4), the tangent map (3.6) at some base point (a,b,a′,b′)∈gl(r)2⊕gl(d−r)2(a,b,a^{\prime},b^{\prime})\in\mathfrak{gl}(r)^{2}\oplus\mathfrak{gl}(d-r)^{2} reads

The determinant polynomial δIJK=det⁡VIJK\delta_{IJK}=\det V_{IJK} is non-zero, i.e., HIJKH_{IJK} is a Ressayre element.

By Sard’s theorem and lemma 7.1, there exists a point where the differential of RIJKR_{IJK} 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 (x,y,z,p,q)(x,y,z,p,q) implies surjectivity of the tangent map at (a,b,a′,b′)(a,b,a^{\prime},b^{\prime}), where p=(aa′′0a′)p=\left(\begin{smallmatrix}a&a^{\prime\prime}\\ 0&a^{\prime}\end{smallmatrix}\right) and q=(bb′′0b′).q=\left(\begin{smallmatrix}b&b^{\prime\prime}\\ 0&b^{\prime}\end{smallmatrix}\right). Thus HIJKH_{IJK} is a Ressayre element. ∎

Let d=6d=6, r=3r=3, and I=J=K={1<3<5}I=J=K=\{1<3<5\}, so that Ic=Jc=Kc={2<4<6}I^{c}=J^{c}=K^{c}=\{2<4<6\}. Then λI=λJ=λK=(3,2,1)\lambda_{I}=\lambda_{J}=\lambda_{K}=(3,2,1) and λIc=λJc=λKc=(2,1,0)\lambda_{I^{c}}=\lambda_{J^{c}}=\lambda_{K^{c}}=(2,1,0), 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 R(gl(3)2)R(\mathfrak{gl}(3)^{2}). Note that F1(x,y):=f1(x)f2(y)h(y,x)F_{1}(x,y):=f_{1}(x)f_{2}(y)h(y,x) and F2(y,x):=F1(x,y)F_{2}(y,x):=F_{1}(x,y) span the two-dimensional subspace of weight ((−2,−1,0)((-2,-1,0), (−2,−1,0)(-2,-1,0), (1,2,3))(1,2,3)). An explicit calculation best left to a computer algebra system Stein et al. (2014) verifies that

Therefore, in agreement with (7.3) δIJK\delta_{IJK} is indeed a lowest weight vector of weight

Note that if we consider a′a^{\prime} and b′b^{\prime} as coefficients rather than indeterminates, δIJK(−,a′,b′)\delta_{IJK}(-,a^{\prime},b^{\prime}) spans the two-dimensional subspace of lowest weight vectors of weight ((−3,−2,−1),(−3,−2,−1),(3,4,5))((-3,-2,-1),(-3,-2,-1),(3,4,5)) as we vary a′a^{\prime} and b′b^{\prime} over gl(3)\mathfrak{gl}(3). Likewise, δIJK(a,b,−)\delta_{IJK}(a,b,-) spans the subspace of lowest weight vectors of weight ((−2,−1,0),(−2,−1,0),(1,2,3))((-2,-1,0),(-2,-1,0),(1,2,3)) if we instead vary aa and bb.

where (D,E,F)=(λIc,λJc,λKc∗+rχd−r)(D,E,F)=(\lambda_{I^{c}},\lambda_{J^{c}},\lambda_{K^{c}}^{*}+r\chi_{d-r}). 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 δIJK≠0\delta_{IJK}\neq 0 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 UAU_{A}, UBU_{B}, and UCU_{C} are scalars (i.e., proportional to the identity matrix).

Now suppose that UAU_{A}, UBU_{B}, and UCU_{C} are unitaries such that (A.1) holds. Then,

Let 1<a≤b≤c≤ab1<a\leq b\leq c\leq ab. Then the Kronecker cone C(a,b,c)C(a,b,c) is maximal-dimensional.

In the case where a=b=ca=b=c the following lemma strengthens corollary A.2:

for any triple of small perturbations εA,εB,εC\varepsilon_{A},\varepsilon_{B},\varepsilon_{C}. Therefore, the only constraints in the vicinity of (τd,τd,τd)(\tau_{d},\tau_{d},\tau_{d}) are the Weyl chamber inequalities. ∎

References