---
title: Pseudomoment Cone and Hankel Spectrahedra
url: https://www.emergentmind.com/topics/pseudomoment-cone
type: topic
---

# Pseudomoment Cone and Hankel Spectrahedra

The homogeneous pseudo-moment cone $\Sigma_{n,2d}^*$ is the dual, in the sense of convex cones of linear functionals, of the sum-of-squares cone $\Sigma_{n,2d}$ of degree-$2d$ forms in $n$ variables. In matrix form, it is realized as a Hankel spectrahedral cone, and its facial geometry governs moment/SOS relaxations, the structure of truncated moment sequences, and the identifiability of atomic rank-one decompositions of moment matrices [2605.06854]. In SOS and moment hierarchies it also appears as a spectrahedral shadow, so questions about its extremal structure interact with the general geometry of linear images of the PSD cone, including possible non-closedness and associated boundary pathologies [2006.09956].

## 1. Definition and dual realizations

Let $H_{n,d}$ denote the real vector space of homogeneous polynomials of degree $d$ in $n$ variables. The cone of SOS forms of degree $2d$ is
$$
\Sigma_{n,2d}:=\{f\in H_{n,2d}\mid f=\sum_i q_i^2,\ q_i\in H_{n,d}\},
$$
and the cone of nonnegative forms is
$$
P_{n,2d}:=\{f\in H_{n,2d}\mid f(x)\ge 0,\ \forall x\in\mathbb{R}^n\}.
$$
Their dual cones of linear functionals $\ell:H_{n,2d}\to\mathbb{R}$ are
$$
P_{n,2d}^*:=\{\ell\in H_{n,2d}^*\mid \ell(p)\ge 0,\ \forall p\in P_{n,2d}\},
$$
and
$$
\Sigma_{n,2d}^*:=\{\ell\in H_{n,2d}^*\mid \ell(q^2)\ge 0,\ \forall q\in H_{n,d}\}.
$$
The dual $\Sigma_{n,2d}^*$ is called the homogeneous pseudo-moment cone [2605.06854].

By construction, $P_{n,2d}\supseteq \Sigma_{n,2d}$ and therefore $P_{n,2d}^*\subseteq \Sigma_{n,2d}^*$. This already separates the pseudo-moment cone from the cone coming from genuine measures: pseudo-moments are required to be nonnegative on squares, but not every such functional is induced by a representing measure [2605.06854].

A monomial basis is central to the standard matrix realization. For multi-indices $\alpha=(\alpha_1,\dots,\alpha_n)$ with $|\alpha|:=\sum_i \alpha_i=d$, the degree-$d$ monomial vector is
$$
v_d(x):=(x^\alpha)_{\alpha\in\mathbb{N}^n,\ |\alpha|=d}\in\mathbb{R}^N,
\qquad
N=\binom{n+d-1}{d}.
$$
The inhomogeneous vector of monomials up to degree $d$ is
$$
m_d(z):=(z^\alpha)_{\alpha\in\mathbb{N}^n,\ |\alpha|\le d}\in\mathbb{R}^{\binom{n+d}{d}},
$$
and homogenization gives $m_d(z)=v_d(\tilde z)$ with $\tilde z=[1;z]\in\mathbb{R}^{n+1}$ [2605.06854].

The SOS cone has a Gram representation:
$$
f(x)=v_d(x)^\top Q\,v_d(x),\qquad Q\succeq 0,
$$
so
$$
\Sigma_{n,2d}=\{v_d^\top Q v_d\mid Q\succeq 0\}.
$$
Dually, a linear functional $\ell$ is encoded by a homogeneous degree-$2d$ pseudo-moment sequence $y=(y_\gamma)$, with $\ell(x^\gamma)=y_\gamma$, and the associated homogeneous moment matrix is
$$
[M_d(y)]_{\alpha,\beta}=y_{\alpha+\beta},\qquad |\alpha|=|\beta|=d.
$$
Then
$$
M_d(y)\succeq 0\iff \ell(q^2)\ge 0\ \forall q\in H_{n,d},
$$
hence
$$
\Sigma_{n,2d}^*=\{\ell\mid M_d(y)\succeq 0\}
$$
[2605.06854].

## 2. Hankel spectrahedra, representable moments, and atoms

The matrix realization of the pseudo-moment cone is
$$
\Sigma_{n,2d}^{\mathrm{mat}}
=
\{M(\ell)\mid \ell\in \Sigma_{n,2d}^*\}
=
H_{n,d}^{\mathrm{Hankel}}\cap S_+^N,
$$
where
$$
H_{n,d}^{\mathrm{Hankel}}
:=
\{X\in S^N\mid X_{\alpha,\beta}\text{ depends only on }\alpha+\beta\}.
$$
Equivalently,
$$
X_{\alpha,\beta}=X_{\alpha',\beta'}\quad\text{whenever}\quad \alpha+\beta=\alpha'+\beta'.
$$
Thus the pseudo-moment cone is a spectrahedral cone obtained by intersecting the PSD cone with a Hankel linear subspace [2605.06854].

By contrast, the representable moment-matrix cone is
$$
M_{n,2d}
=
\{M(\ell)\mid \ell\text{ is a positive linear functional induced by a measure}\}
=
\operatorname{cone}\{v_d(z)v_d(z)^\top\mid z\in\mathbb{R}^n\},
$$
and one has
$$
M_{n,2d}\subseteq \Sigma_{n,2d}^{\mathrm{mat}}.
$$
A common misconception is that the pseudo-moment cone and the representable moment cone coincide. The inclusion above shows that representable moments form only a subcone of the full pseudo-moment cone [2605.06854].

Atoms arise from point evaluations. For any $a\in\mathbb{R}^n$, the functional $\ell_a(p):=p(a)$ belongs to $\Sigma_{n,2d}^*$, and
$$
M(\ell_a)=v_d(a)v_d(a)^\top,
$$
a rank-one atomic moment matrix. Weighted atoms are finite conic sums
$$
\sum_i w_i\,v_d(a_i)v_d(a_i)^\top,\qquad w_i>0.
$$
Every rank-one extreme ray of $\Sigma_{n,2d}^*$ is generated by a scaled point evaluation $c\cdot \ell_a$, $c\ge 0$ [2605.06854].

This rank-one description is only part of the global extremal structure. Higher-rank extreme rays do exist globally in $\Sigma_{n,2d}^*$, including rank $6$ for $(d=2,n\ge 4)$ and rank $3d-2$ for $d\ge 3,\ n\ge 3$ [2605.06854]. The significance of the recent local theory is not the absence of such rays globally, but the fact that they are excluded from certain minimal faces around generically generated atomic moment matrices.

## 3. Facial geometry and simplicial regularizability

The facial geometry of $\Sigma_{n,2d}^{\mathrm{mat}}$ is inherited from the PSD cone and then refined by Hankel constraints. If
$$
X=V\Lambda V^\top
$$
with $V\in\mathbb{R}^{N\times r}$ full column rank and $\Lambda\succeq 0$, then
$$
\minface(X,S_+^N)=V S_+^r V^\top.
$$
Intersecting with the Hankel subspace yields
$$
\minface(X,\Sigma_{n,2d}^{\mathrm{mat}})
=
(V S_+^r V^\top)\cap H_{n,d}^{\mathrm{Hankel}}
$$
[2605.06854].

The central geometric result of "Simplicial Regularizability of the Pseudo-Moment Cone and Carathéodory-Type Atomic Decomposition of Moment Matrices" establishes a local simplicial structure in the inhomogeneous-via-homogenization model [2605.06854]. Let
$$
\Sigma_{n+1,2d}^{\mathrm{mat}}
:=
S_+^{\binom{n+d}{d}}\cap H_{n+1,d}^{\mathrm{Hankel}},
$$
where $m_d(z)=v_d([1;z])$. If the number of atoms $s$ satisfies
$$
s\le \frac{1}{2}\left(\binom{n+d}{d}+1\right)-\binom{n+d}{d}^{-1}\binom{n+2d}{2d},
$$
then for fixed $d\ge 2$ this upper bound scales as $\Theta(n^d)$, and for generically chosen $\{(w_i,z_i)\}_{i=1}^s\in (\mathbb{R}_{>0}\times\mathbb{R}^n)^{\times s}$,
$$
\minface\!\left(
\sum_{i=1}^s w_i^2\,m_d(z_i)m_d(z_i)^\top,\,
\Sigma_{n+1,2d}^{\mathrm{mat}}
\right)
=
\operatorname{cone}\bigl(\{m_d(z_i)m_d(z_i)^\top\}_{i=1}^s\bigr),
$$
and the set $\{m_d(z_i)m_d(z_i)^\top\}_{i=1}^s$ is linearly independent [2605.06854].

The theorem is a local facial statement. Its interpretation is that the minimal face containing
$$
X=\sum_i w_i^2\,m_d(z_i)m_d(z_i)^\top
$$
is a simplicial cone generated exactly by the planted rank-one atoms. Near such generically generated $X$, the pseudo-moment cone locally regularizes to a cone isomorphic to $\mathbb{R}_+^s$, enabling unique atomic recovery [2605.06854].

This local regularity does not contradict the existence of higher-rank extreme rays. The same source states explicitly that higher-rank extreme rays exist globally, but they do not appear in the minimal face of $X$ under the stated genericity and $O(n^d)$ bound [2605.06854]. This distinction is crucial for interpreting the result: it concerns a generic atomic regime, not the entirety of the cone.

## 4. Carathéodory-type extreme-ray decomposition

The same work develops a general decomposition procedure for spectrahedral cones, denoted RayDecomp [2605.06854]. Given a nonempty spectrahedral cone $C\subseteq S^N$ and a nonzero point $X\in C$, the algorithm returns coefficients and extreme rays
$$
X=\sum_{k=1}^r t_k M_k,\qquad t_k>0,
$$
with each $M_k$ generating an extreme ray of $C$.

Its procedure is:

1. Initialize $k:=1$, $X_1:=X$.
2. While $X_k\neq 0$:
   1. $F_k:=\minface(X_k,C)$.
   2. Draw $B_k$ from any absolutely continuous distribution on $S^N$.
   3. Solve
      $$
      M_k\in \arg\min_{Y\in F_k,\ \operatorname{tr}(Y)=1}\langle B_k,Y\rangle,
      $$
      and set $M_k$ to the minimizer.
   4. Compute
      $$
      t_k:=\sup\{t\ge 0\mid X_k-tM_k\in F_k\}.
      $$
   5. Update $X_{k+1}:=X_k-t_k M_k$.
3. Return the decomposition [2605.06854].

If each SDP is solved exactly, then each step is well-defined, and with probability $1$ the algorithm terminates in finitely many steps and returns an extreme-ray decomposition of $X$ [2605.06854]. The probabilistic mechanism is geometric: minimizing a random linear functional over the normalized face almost surely exposes a generator of an extreme ray.

When specialized to the pseudo-moment cone in the generic simplicial regime, the algorithm provably recovers the planted atoms and weights uniquely up to permutation. If
$$
X=\sum_{i=1}^s w_i^2\,m_d(z_i)m_d(z_i)^\top
$$
with $s$ in the theorem’s bound and $\{(w_i,z_i)\}$ generic, then with probability $1$
$$
r=s,\qquad
t_k=w_k^2\|m_d(z_k)\|_2^2,\qquad
M_k=\frac{1}{\|m_d(z_k)\|_2^2}m_d(z_k)m_d(z_k)^\top
$$
[2605.06854].

For a rank-one output $M_k$, if $\lambda_k$ is its unique nonzero eigenvalue and $h_k$ is the associated unit eigenvector, then
$$
m_d(z_k)=\frac{1}{h_{k,1}}\,h_k,
$$
and $z_k$ is read off from the degree-$1$ entries of $m_d(z_k)$. The planted weight recovers as
$$
w_k^2=t_k h_{k,1}^2
$$
[2605.06854].

Compared to flatness-based extraction on degree-$2d$ truncated moment matrices, which can certify at most $\Theta(n^{d-1})$ atoms from degree-$2d$ data, this method recovers up to $\Theta(n^d)$ atoms without solving higher-degree relaxations [2605.06854]. This suggests a different route to identifiability: not via higher-order flat extension, but via local simplicial facial geometry.

## 5. Numerical stabilization, complexity, and empirical behavior

A finite-precision robust implementation, RayDecompRestart, is described for the pseudo-moment setting [2605.06854]. It includes facial reduction at each iteration: the algorithm restricts to the PSD face determined by the numerical range space of $X_k$, truncates small eigenvalues, removes redundant linear constraints via pivoted QR, and solves the reduced SDP in the face. The stated purpose is to restore strict feasibility and numerical stability.

The implementation also replaces the one-dimensional feasibility search by the closed-form step size
$$
t_k=\lambda_{\max}^{-1}\!\bigl(X_k^{-1/2}M_kX_k^{-1/2}\bigr),
$$
uses alternating projections onto the linear subspace $\{X\mid \mathcal{A}(X)=0\}$ and onto $S_+^N$, and employs a restart mechanism that detects SDP failures, reduces rank-thresholds, reprojects, and restarts [2605.06854].

For $\Sigma_{n+1,2d}^{\mathrm{mat}}$, the ambient matrix size is
$$
N=\binom{n+d}{d}=\Theta(n^d)
$$
for fixed $d$. One inner-loop iteration costs, in the worst case,
$$
O(r_{\max}^6+r_{\max}^4N^2+N^3)
$$
time and
$$
O(r_{\max}^4+r_{\max}^2N^2)
$$
memory, where $r_{\max}\le N$ is the numerical rank of the input $X$. The total cost over $O(r_{\max})$ iterations is
$$
O\bigl(r_{\max}(r_{\max}^6+r_{\max}^4N^2+N^3)\bigr).
$$
In the worst case $r_{\max}=N$, this becomes
$$
O(N^7)\sim O(n^{7d})
$$
time and
$$
O(N^4)\sim O(n^{4d})
$$
memory [2605.06854].

The reported numerical behavior is twofold. First, a stabilized implementation using the MOSEK SDP solver, facial reduction, and alternating projections exhibits strong recovery across a range of $(n,d)$, with near-perfect recovery below a data-driven phase-transition threshold in $s$, often substantially larger than the conservative theoretical bound [2605.06854]. Second, outside the guaranteed regime, the algorithm returns non-unique extreme rays, and empirically high-rank extreme rays appear once $s$ crosses the phase-transition threshold; one cited example is rank-$6$ rays for $(d=2,n=3)$ in the homogenized model [2605.06854]. This suggests that the same procedure can function as a practical sampler of high-rank extreme rays, although that role is empirical rather than part of the proved guarantee.

## 6. Spectrahedral shadows, closedness, and broader convex-algebraic context

Pseudo-moment cones also arise as linear images of PSD cones. In SOS and moment hierarchies one works with moment or pseudo-moment matrices $Y\in S_+^N$ subject to linear constraints, and the truncated moment sequences extracted from $Y$ form a spectrahedral shadow
$$
C_d=M_d(S_+^N)=\pi_{\mathcal{L}}(S_+^N)
$$
for an appropriate subspace $\mathcal{L}\subset S^N$ [2006.09956]. Consequently, the geometry of pseudo-moment cones is linked to the general theory of projections of the PSD cone.

"Bad Projections of the PSD Cone" studies when such linear images fail to be closed [2006.09956]. For a subspace $\mathcal{L}\subset S^n$ with associated map $\pi_{\mathcal{L}}$, the image cone
$$
C=\pi_{\mathcal{L}}(S_+^n)
$$
need not be closed, even though $S_+^n$ itself is closed. The closure satisfies
$$
\overline{\pi_{\mathcal{L}}(S_+^n)}
\simeq
(\mathcal{L}\cap S_+^n)^\vee
$$
[2006.09956]. In the context of pseudo-moment cones, this means that some boundary points of the closure may not be attainable by actual PSD moment matrices after projection.

The paper provides several equivalent diagnostics for badness, including Pataki’s block characterization, an intrinsic criterion in terms of spectrahedral ranks and ideals, and a normal-cycle criterion based on complementarity [2006.09956]. It also identifies the Zariski closure of the bad locus as a hypersurface in the Grassmannian whose irreducible components are coisotropic hypersurfaces of symmetric determinantal varieties [2006.09956]. These results apply verbatim to pseudomoment cones viewed as spectrahedral shadows.

The practical implications are explicit. Non-closed pseudomoment cones lead to weak infeasibility and failure of strong duality in relaxations, resulting in duality gaps or numerical instability [2006.09956]. This is a different phenomenon from the local simplicial regularizability established in the Hankel spectrahedral realization of $\Sigma_{n,2d}^*$ [2605.06854]. A plausible implication is that the pseudo-moment cone should be understood through two complementary geometric lenses: as an intersection $H_{n,d}^{\mathrm{Hankel}}\cap S_+^N$, where local facial structure can become simplicial around generic atomic points, and as a spectrahedral shadow, where projection-induced non-closedness can create global boundary pathologies.

Within SOS optimization and moment hierarchies, these two perspectives address different technical issues. The simplicial regularizability result explains identifiability and decomposition in a generic $O(n^d)$ atomic regime [2605.06854]. The bad-projection framework explains when projected pseudomoment sets may violate Slater-type regularity and require facial reduction or extended dual constructions [2006.09956]. Taken together, they place the pseudomoment cone at the intersection of convex algebraic geometry, spectrahedral facial theory, and algorithmic moment decomposition.

Source: https://www.emergentmind.com/topics/pseudomoment-cone