---
title: Two-Dimensional Moment Problem
url: https://www.emergentmind.com/topics/two-dimensional-moment-problem
type: topic
---

# Two-Dimensional Moment Problem

The two-dimensional moment problem asks when a family of real numbers indexed by \(\mathbb{N}^2\) is representable as moments of a positive measure on \(\mathbb{R}^2\), either on the full plane or on a prescribed support set \(K\subseteq \mathbb{R}^2\). In its basic form, one seeks a positive Radon or Borel measure \(\mu\) such that
\[
s_{\alpha}=\int_{\mathbb{R}^2} x_1^{\alpha_1}x_2^{\alpha_2}\,d\mu(x_1,x_2),\qquad \alpha=(\alpha_1,\alpha_2)\in\mathbb{N}^2.
\]
The corresponding \(K\)-moment problem requires \(\operatorname{supp}\mu\subseteq K\). In two variables, the subject lies at the intersection of analysis, real algebraic geometry, operator theory, convex geometry, and semidefinite optimization, and it includes both full and truncated formulations [2604.12265].

## 1. Classical formulation and positivity

A two-dimensional moment sequence is a family \(s=(s_\alpha)_{\alpha\in\mathbb{N}^2}\), and it defines a linear functional \(L:\mathbb{R}[x_1,x_2]\to\mathbb{R}\) by
\[
L(p)=\sum_{\alpha} p_\alpha s_\alpha=\int_{\mathbb{R}^2} p(x_1,x_2)\,d\mu(x_1,x_2),
\]
for \(p(x_1,x_2)=\sum_\alpha p_\alpha x_1^{\alpha_1}x_2^{\alpha_2}\). The central existence statement is Haviland’s theorem in dimension two: if \(K\subseteq\mathbb{R}^2\) is closed, then \(L\) admits a representing positive Radon measure supported on \(K\) if and only if \(L(p)\ge 0\) for every polynomial \(p\) with \(p|_K\ge 0\). For \(K=\mathbb{R}^2\), this becomes positivity on the cone \(\operatorname{Pos}(\mathbb{R}^2)=\{p\in\mathbb{R}[x_1,x_2]:p\ge 0 \text{ on }\mathbb{R}^2\}\) [2604.12265].

In truncated form, only finitely many moments are prescribed. For degree \(2d\), one asks for \(\mu\) such that
\[
y_\alpha=\int x^\alpha\,d\mu(x),\qquad |\alpha|\le 2d,
\]
with \(x^\alpha=x_1^{\alpha_1}x_2^{\alpha_2}\). The truncated moment matrix \(M_d(y)\), indexed by \(A_d=\{\alpha\in\mathbb{N}^2:|\alpha|\le d\}\), is given by
\[
M_d(y)_{\alpha,\beta}=y_{\alpha+\beta}.
\]
Its positivity is equivalent to \(L_y(f^2)\ge 0\) for all polynomials \(f\) of degree \(\le d\) [1910.05251].

A persistent structural distinction from one dimension is that the positivity of a single Hankel or Toeplitz object does not generally solve the support-constrained problem in higher dimensions. One-dimensional Hamburger and circular moment problems admit such criteria, whereas in two dimensions “except some fortunate configurations, in higher dimensions there are no non-negativity criteria for the power moments of a measure to be supported by a prescribed closed set” [1910.05251].

## 2. Semialgebraic supports, quadratic modules, and positivity certificates

When the support is a basic closed semialgebraic set
\[
K=\{x\in\mathbb{R}^2:g_i(x)\ge 0,\ i=1,\dots,m\},
\]
the analytic cone of polynomials nonnegative on \(K\) is approximated algebraically by two standard cones. The quadratic module is
\[
M(g)=\Big\{\sigma_0+\sum_{i=1}^m \sigma_i g_i:\ \sigma_i\text{ are sums of squares in }\mathbb{R}[x_1,x_2]\Big\},
\]
and the preordering is
\[
T(g)=\Big\{\sum_{J\subseteq\{1,\dots,m\}} \sigma_J\,g^J:\ \sigma_J\text{ sums of squares},\ g^J:=\prod_{i\in J}g_i\Big\}.
\]
One always has \(M(g)\subseteq T(g)\subseteq \{p:p|_K\ge 0\}\) [2604.12265].

Two Positivstellensätze govern the compact case. Schmüdgen’s theorem states that if \(K\) is compact and \(p>0\) on \(K\), then \(p\in T(g)\). Putinar’s theorem states that if \(M(g)\) is Archimedean, meaning there exists \(N>0\) such that
\[
N-x_1^2-x_2^2\in M(g),
\]
then \(p>0\) on \(K\) implies \(p\in M(g)\). The Archimedean condition is the “algebraic shadow of compactness”: it implies \(K\subseteq\{x:x_1^2+x_2^2\le N\}\), hence compactness. The preordering \(T(g)\) is always Archimedean when \(K\) is compact, whereas \(M(g)\) may fail to be Archimedean unless the generators already include a coordinate bound such as a ball constraint [2604.12265].

The practical distinction is decisive. \(T(g)\) uses all products \(g^J\) and therefore has exponentially many generators in \(m\), while \(M(g)\) uses only linear combinations of the \(g_i\) with SOS weights. Under Archimedean \(M(g)\), semidefinite relaxations built from \(M(g)\) are typically cheaper and converge in Lasserre’s hierarchy; Schmüdgen’s route guarantees convergence under mere compactness at the cost of larger cones [2604.12265].

The standard two-variable examples are explicit. On the unit disk,
\[
K=\{1-x_1^2-x_2^2\ge 0\},
\]
\(M(g)\) is Archimedean with \(N=1\), so every polynomial strictly positive on the disk has a Putinar certificate
\[
p=\sigma_0+\sigma_1(1-x_1^2-x_2^2).
\]
On the box \(K=[-1,1]^2\), the generators \(g_1=1-x_1^2\), \(g_2=1-x_2^2\) satisfy
\[
2-(x_1^2+x_2^2)=g_1+g_2\in M(g),
\]
so again \(M(g)\) is Archimedean [2604.12265].

## 3. Operator-theoretic constructions and analytic parametrizations

A complementary route starts from the sesquilinear form induced by the moments. Given \(L:\mathbb{R}[x_1,x_2]\to\mathbb{R}\), define
\[
\langle p,q\rangle:=L(pq),\qquad N:=\{p:L(p^2)=0\}.
\]
The quotient \(\mathbb{R}[x_1,x_2]/N\) becomes a pre-Hilbert space, whose completion \(H_L\) carries densely defined symmetric multiplication operators
\[
X_i(p+N)=x_i p+N,\qquad i=1,2.
\]
Under compactness or Archimedean hypotheses ensuring boundedness, these extend to bounded self-adjoint commuting operators, and the joint spectral theorem yields a measure \(\mu\) on \(\mathbb{R}^2\) with
\[
L(x_1^{\alpha_1}x_2^{\alpha_2})=\langle X_1^{\alpha_1}X_2^{\alpha_2}v,v\rangle
=\int_{\mathbb{R}^2} x_1^{\alpha_1}x_2^{\alpha_2}\,d\mu(x),
\]
where \(v=1+N\). If \(L\) is nonnegative on \(T(g)\), or on \(M(g)\) under Putinar’s hypothesis, then \(\operatorname{supp}(\mu)\subseteq K\) [2604.12265].

A more elaborate operator model solves the unconstrained two-dimensional Hamburger problem through an extended moment problem indexed by
\[
(m,k,\ell;\,n,r,t)\in \mathcal{Q}
\]
and moments
\[
u_{m,k,\ell;\,n,r,t}
=
\int_{\mathbb{R}^2}
x_1^m (x_1+i)^k (x_1-i)^\ell
x_2^n (x_2+i)^r (x_2-i)^t\,d\mu.
\]
The extended problem has a solution if and only if the positivity condition (P-Ext) and the difference relations (9)–(12), equivalently (13)–(16), hold; in that case the solution is unique. The original two-dimensional moment problem is solvable if and only if the given \(s_{m,n}\) can be embedded as \(u_{m,0,0;\,n,0,0}\) in such an extended sequence, and all solutions are then obtained by the paper’s iterative Hilbert-space algorithm [1006.5803]. A second exposition of the same framework gives the same criterion and extends it to the complex moment problem through explicit conversion formulas between \(a_{m,n}\) and the real moments \(s_{p,q}\) [1009.4821].

Support restrictions can be incorporated operator-theoretically. For the strip
\[
\Pi(R)=\{(x_1,x_2)\in\mathbb{R}^2:|x_2|\le R\},
\]
solvability is equivalent to the positivity of the kernel
\[
K((m,n),(k,l))=s_{m+k,n+l}
\]
together with the boundedness inequality
\[
\sum a_{m,n}\overline{a_{k,l}}\,(R^2 s_{m+k,n+l}-s_{m+k,n+l+2})\ge 0,
\]
for all finitely supported coefficients. All solutions are parameterized by spectral measures of one multiplication operator commuting with the spectral measure of the other, via generalized resolvents and decomposable contractions [1005.1769].

Under the Carleman-type condition
\[
\sum_{k=1}^{\infty}\frac{1}{\sqrt[2k]{\,s_{2m,2k}+s_{2m+2,2k}\,}}=\infty,\qquad m=0,1,2,\dots,
\]
all solutions of the two-dimensional moment problem are parameterized by analytic contractive operator-valued functions
\[
\Phi\in \mathcal{S}_{V,U}(\mathbb{D};N_0(V),N_\infty(V))
\]
satisfying an intertwining constraint with the Cayley transform \(U\) of the second multiplication operator. In this setting, the determinate case is characterized by self-adjointness of the first multiplication operator \(A_1\), equivalently by the existence of a unique canonical solution [2508.10823].

## 4. Truncated problems, flat extensions, and canonical solutions

For truncated data \(s=(s_\alpha)\) with \(|\alpha|\le 2n\), the moment matrix is
\[
[M_n(s)]_{\alpha,\beta}=s_{\alpha+\beta},\qquad |\alpha|,|\beta|\le n.
\]
If \(K=\{g_i\ge 0\}\), each constraint polynomial \(g_i\) produces a localizing matrix
\[
[M_{n-d_i}(g_i s)]_{\alpha,\beta}
=
\sum_\gamma (g_i)_\gamma\, s_{\alpha+\beta+\gamma},
\qquad d_i=\left\lceil \frac{\deg(g_i)}{2}\right\rceil .
\]
The PSD conditions
\[
M_n(s)\succeq 0,\qquad M_{n-d_i}(g_i s)\succeq 0
\]
are necessary for a \(K\)-supported representing measure, and under compact or Archimedean hypotheses they become asymptotically sufficient in the moment-SOS hierarchy [2604.12265].

The decisive finite-atomic criterion is flatness:
\[
\operatorname{rank} M_n(s)=\operatorname{rank} M_{n-1}(s).
\]
If this holds together with \(M_n(s)\succeq 0\), then \(s\) admits a finitely atomic representing measure
\[
\mu=\sum_{k=1}^r w_k\,\delta_{\xi^{(k)}},\qquad r=\operatorname{rank}M_n(s).
\]
The kernel ideal \(\ker M_n(s)\) determines algebraic relations satisfied by the atoms, and multiplication matrices on the quotient recover the atom locations through joint eigenvalues; the weights follow from a Vandermonde or evaluation system [2604.12265].

For triangular truncations
\[
K_r=\{(k,\ell)\in\mathbb{Z}_+^2:k+\ell\le r\},
\]
canonical solutions are defined as those generated by commuting self-adjoint extensions inside the associated Hilbert space. There is a one-to-one correspondence between canonical solutions and dimensionally stable close extensions, and in the two-dimensional triangular case this is a one-to-one correspondence between canonical solutions and flat extensions of the given data. The paper introduces the index of nonself-adjointness
\[
i_s=\dim(H\ominus H_0),
\]
with \(i_s=0\) corresponding to flatness. For \(i_s=1\), explicit necessary and sufficient conditions are obtained; for \(i_s=2\), the existence question is reduced either to explicit conditions or to a single quadratic equation with several unknowns [2406.19916].

Rectangular truncations admit different constructive criteria. For the truncated two-dimensional problem with data \(\{s_{m,n}\}_{0\le m\le M,\ 0\le n\le N}\), one reduction places the representing measure on a finite union of parallel lines
\[
K_N(a_0,\dots,a_N)=\bigcup_{j=0}^N \{(x_1,x_2):x_2=a_j\},
\]
and reduces solvability to \(N+1\) one-dimensional truncated Hamburger problems via a Vandermonde system. This yields explicit necessary and sufficient conditions for \(M=N=1\) and \(M=1,N=2\), and explicit sufficient conditions for \(M=N=2\), \(M=2,N=3\), \(M=3,N=2\), and \(M=N=3\) [1707.09501].

## 5. Determinacy, algebraic varieties, and special support geometries

Determinacy means uniqueness of the representing measure. For compact \(K\), determinacy is automatic. On \(\mathbb{R}^2\), one sufficient condition is the two-coordinate Carleman divergence
\[
\sum_{n=1}^{\infty} s_{(2n,0)}^{-1/(2n)}=\infty,
\qquad
\sum_{n=1}^{\infty} s_{(0,2n)}^{-1/(2n)}=\infty.
\]
More generally, if the one-dimensional marginals are determinate, then the measure is determinate; this is Petersen’s principle. Operator-theoretic quasi-analyticity and essential self-adjointness criteria also imply determinacy [2604.12265].

For the Stieltjes problem on \(\mathbb{R}_+^2\), positivity
\[
L(p^2)\ge 0,\qquad L(x_1 p^2)\ge 0,\qquad L(x_2 p^2)\ge 0
\]
together with the multivariate Stieltjes divergence
\[
\sum_{n=1}^{\infty} \big(L(x_1^n)\big)^{-1/(2n)}=\infty,
\qquad
\sum_{n=1}^{\infty} \big(L(x_2^n)\big)^{-1/(2n)}=\infty
\]
implies existence and uniqueness of a representing measure on \(\mathbb{R}_+^2\). The same paper extends this to unbounded semialgebraic sets \(K(f)=\{x:f_j(x)\ge 0\}\), assuming positivity on \(p^2\) and \(f_j p^2\) plus divergence along the generators \(f_j\) [2011.04025].

Algebraic varieties provide a particularly transparent support condition. If \(\operatorname{supp}\mu\subseteq V=\{h(x_1,x_2)=0\}\), then one can encode \(V\) as the basic closed set \(K=\{h\ge 0,-h\ge 0\}\), and
\[
T(h,-h)=\Sigma+I(V).
\]
Consequently, if \(p\ge 0\) on \(V\), then for every \(\varepsilon>0\) there exist \(\sigma\in\Sigma\) and \(q\in I(V)\) such that \(p+\varepsilon=\sigma+q\). In the moment problem on \(V\), positivity on the variety is thus captured algebraically up to an arbitrarily small perturbation [2604.12265].

A substantial extension of this principle is the class of “curves with bumps.” If \((q)\subset\mathbb{R}[x_1,x_2]\) is a real principal ideal whose zero set \(V(q)\) satisfies the Scheiderer–Plaumann conditions, and \(Q\) is an Archimedean quadratic module, then
\[
\Sigma^2+qQ
\]
has the strong moment property, and its positivity set is
\[
\operatorname{Pos}(\Sigma^2+qQ)=V(q)\cup [K_Q\cap\{q>0\}].
\]
This yields non-negativity certificates and truncated SDP conditions for measures supported on a curve together with compact semialgebraic “bumps” [1910.05251].

Certain plane curves admit exact reduction to one-dimensional strong Hamburger problems. For the variety \(xy=1\), the two-dimensional truncated problem reduces to a strong truncated Hamburger problem with Laurent moments, and existence of a representing measure is equivalent to existence of a flat extension of the moment matrix for all degrees \(k\). For \(x^2y=1\), the analogous reduction leads to a strong truncated Hamburger problem with one missing moment, which is solved through partial PSD conditions and Schur-complement formulas [2101.00486].

## 6. Convex geometry, atom counts, and computational frameworks

In finite-dimensional truncated settings, the set of moment functionals forms a convex cone. For a finite-dimensional real vector space \(V\subset C(K;\mathbb{R})\), the moment cone is
\[
\mathcal{M}(V,K)=\{L\in V^*:\exists\,\mu\ge 0,\ L(f)=\int_K f\,d\mu\ \forall f\in V\}.
\]
Richter–Tchakaloff implies that every moment functional has an atomic representing measure with at most \(\dim V\) atoms [1703.01497]. The moment cone is convex, and for differentiable \(V\) one can analyze faces, exposed faces, inner points, and singularities through the differential of the atomic moment map
\[
S_k(C,X)=\sum_{j=1}^k C_j S_F(x_j)
\]
and through the sets \(W(L)\), \(V_+(L)\), and the core variety \(\mathrm{V}(L)\), with
\[
W(L)=\mathrm{V}(L)
\]
for nonzero \(L\) [1703.01497].

For \(A=\mathbb{R}[x_1,x_2]_{\le d}\), the Carathéodory number \(\mathcal{C}_A\) measures the maximal minimal number of atoms needed to represent an \(A\)-moment sequence. If
\[
m=\dim \mathbb{R}[x_1,x_2]_{\le d}=\binom{2+d}{d},
\]
then the differential-geometric lower bound is
\[
N_A=\left\lceil \frac{m}{3}\right\rceil,
\]
except for \(d=2\), where \(N_A=3\), and \(d=4\), where \(N_A=6\). On connected \(X\subseteq\mathbb{R}^2\), one has \(\mathcal{C}_A\le m-1\) [1703.01494, 1809.00584].

Several computational frameworks coexist with moment-SOS methods. In one analytic approach, the two-variable Cauchy or Stieltjes transform
\[
F(z,\zeta)= -\int\!\!\int \frac{1}{(z-x)(\zeta-y)}\,d\mu(x,y)
\]
is expanded at infinity, and a two-variable Schur algorithm produces continued-fraction representations for all solutions of certain truncated problems. In the non-symmetric and symmetric settings, every solution is described by finitely many Schur atoms together with a free remainder \(T\); when \(T\equiv 0\), the fraction terminates and the measure is atomic [2404.03342].

A different computational viewpoint is maximum entropy on bounded domains. On \(\Omega=[-1,1]^2\), the moment-constrained entropy maximizer has density
\[
p_\lambda(x,y)=\frac{1}{Z(\lambda)}\exp\!\left(\sum_{k=1}^M \lambda_k \phi_k(x,y)\right),
\]
and the unknown multipliers solve nonlinear moment-matching equations. The equation-by-equation method solves these constraints sequentially through a homotopy/Newton scheme, can discard infeasible constraints, and in reported two-dimensional tests was more accurate than classical Newton’s method, MATLAB’s generic solver, and a BFGS-based method [1702.02242].

Taken together, these viewpoints show that the two-dimensional moment problem is not a single theorem but a family of closely related existence, uniqueness, extension, and reconstruction problems. Haviland-type positivity gives the abstract criterion; quadratic modules and preorderings supply semialgebraic certificates; GNS and spectral theory reconstruct measures from functionals; flat extensions and multiplication matrices recover atoms from truncated data; and convex, analytic, and optimization frameworks describe the geometry of solutions and the algorithms used to compute them [2604.12265].

Source: https://www.emergentmind.com/topics/two-dimensional-moment-problem