---
title: Wigner Function Moments Method
url: https://www.emergentmind.com/topics/wigner-function-moments-method
type: topic
---

# Wigner Function Moments Method

The expression **Wigner Function Moments method** denotes a family of phase-space techniques in which information is extracted from moments associated with a Wigner object rather than from full pointwise reconstruction. In the literature represented here, this includes spatial-angular and spatio-temporal moment matrices of optical Wigner distributions, moments of powers of a quantum Wigner function used as nonclassicality witnesses, Wigner D-moments on $SO(3)$ used for exact measure recovery, and moment hierarchies derived from Wigner transforms in kinetic and electromagnetic transport theory [2404.06708][2407.12116][1606.05306][1804.04019][2308.12613].

## 1. Definitions and conceptual scope

A common source of confusion is that **“Wigner moments”** are not a single construction. In the optical setting of Bekshaev, Angelsky and Zenkova, moments are taken with respect to a spatio-temporal Wigner distribution $I(x,s;k_x,\Delta k)$ and assembled into centroid and covariance data [2404.06708]. In continuous-variable quantum information, Mallick et al. distinguish ordinary phase-space moments
$$
M_{k,\ell}=\int\!\!\int x^k p^\ell W(x,p)\,dx\,dp
$$
from **Wigner-function moments**
$$
w_n=\iint W^n(x,p)\,dx\,dp,
$$
which are moments of the quasiprobability itself [2407.12116]. In harmonic analysis on $SO(3)$, Filbir and Schröder define moments against Wigner D-functions,
$$
y_{\ell,m,n}=\int_{SO(3)} D^\ell_{m,n}(R)\,d\mu(R),
$$
for recovery of discrete measures [1606.05306].

| Context | Moment object | Role |
|---|---|---|
| Spatio-temporal optics | $M_{ij}=(1/F)\int (P_i-\bar P_i)(P_j-\bar P_j)I(P)\,dP$ | Packet size, spread, correlations, transverse OAM |
| Continuous-variable states | $w_n=\iint W^n\,dx\,dp$ or $M_{k,\ell}=\iint x^k p^\ell W\,dx\,dp$ | Negativity witnessing |
| $SO(3)$ inverse problems | $y_{\ell,m,n}=\int D^\ell_{m,n}(R)\,d\mu(R)$ | Exact recovery of discrete measures |
| Kinetic and lattice Wigner theory | $M^{(n)}=\int p^n W\,dp$ | Hierarchy, closure, and numerical preservation |

This plurality is not merely terminological. The moment objects encode different structures: covariance geometry in optics, quasiprobability negativity in quantum-state diagnostics, low-pass harmonic data on a compact Lie group, or transport closures in kinetic equations. A plausible implication is that the unifying idea is algebraic compression of phase-space information, whereas the analytic content depends strongly on the ambient problem.

## 2. Spatio-temporal optical moment formalism

For quasimonochromatic paraxial beams, the Wigner distribution is introduced from a scalar complex amplitude $u(x)$ as
$$
W(x,k)=\int_{-\infty}^{\infty} u^*(x-\xi/2)\,u(x+\xi/2)\,e^{i k \xi}\,d\xi,
$$
with marginals
$$
\int W(x,k)\,dk=|u(x)|^2,\qquad \int W(x,k)\,dx=|U(k)|^2.
$$
Bekshaev, Angelsky and Zenkova generalize this construction to spatio-temporal light fields by introducing the retarded coordinate $s=z-c t$ and the slowly-varying envelope $u(x,y,z,s)e^{-i\omega t}$. In $1+1$ dimensions, the corresponding $4$D Wigner function is
$$
I(x,s;k_x,\Delta k)=\frac{1}{2\pi}\iint u^*(x-\xi/2,s-\sigma/2)\,u(x+\xi/2,s+\sigma/2)e^{i(k_x\xi+\Delta k\sigma)}\,d\xi\,d\sigma,
$$
and the total action is
$$
F=\int I\,dx\,ds\,dk_x\,d\Delta k.
$$
This construction reduces to the standard spatial Wigner function if the $s$-dependence is ignored and to the time-frequency Wigner function if $x$ is ignored [2404.06708].

The first and second moments are organized through the $4$D “ray” vector
$$
P=(x,s;k_x,\Delta k)^T,
$$
with centroid and moment matrix
$$
\bar P_i=\frac{1}{F}\int P_i\,I(P)\,dP,\qquad
M_{ij}=\frac{1}{F}\int (P_i-\bar P_i)(P_j-\bar P_j)\,I(P)\,dP.
$$
Assuming a centered packet, the real-space expressions include
$$
M_{11}=\frac{1}{F}\int x^2|u|^2\,dx\,ds,\qquad
M_{22}=\frac{1}{F}\int s^2|u|^2\,dx\,ds,
$$
$$
M_{33}=\frac{1}{F}\int |\partial_x u|^2\,dx\,ds,\qquad
M_{44}=\frac{1}{F}\int |\partial_s u|^2\,dx\,ds,
$$
and mixed terms such as
$$
M_{13}=\frac{1}{F}\,\mathrm{Im}\int x\,u^*\,\partial_x u\,dx\,ds.
$$
The diagonal entries characterize spatial size, temporal length, angular spread, and frequency spread, while the off-diagonal blocks $M_{13},M_{14},M_{23},M_{24}$ encode mixed spatio-temporal correlations [2404.06708].

A central physical application is the description of transverse orbital angular momentum in spatio-temporal optical vortices. The transverse OAM about the $(x=0,s=0)$ axis is written as
$$
\Lambda_y=\int (s\,p_x-x\,p_s)\,dx\,ds \simeq \frac{1}{8\pi\omega}\,2\,M_{23},
$$
and in practice the $p_s$ term often vanishes by centroid definition, so that $\Lambda_y\propto M_{23}$. If the field carries a phase singularity in the $(x,s)$ plane, then $M_{23}\neq 0$ and hence $\Lambda_y\neq 0$. The paper also presents non-vortex spatio-temporal wave packets with transverse OAM, so the presence of transverse OAM is not restricted to phase-singular configurations [2404.06708].

The transformation law is especially simple for first-order paraxial systems. If
$$
P_{\mathrm{out}}=H\,P_{\mathrm{in}},
$$
with $H$ built from the usual $2\times 2$ ABCD matrix acting on the $(x,k_x)$ subspace and the $(s,\Delta k)$ subspace left untouched, then
$$
M_{\mathrm{out}}=H\,M_{\mathrm{in}}\,H^T.
$$
In particular,
$$
M_{11}^{\mathrm{out}}=A^2M_{11}^{\mathrm{in}}+2ABM_{13}^{\mathrm{in}}+B^2M_{33}^{\mathrm{in}},
$$
and the transverse-OAM entry obeys
$$
M_{23}^{\mathrm{out}}=k_0 C M_{12}^{\mathrm{in}}+D M_{23}^{\mathrm{in}}.
$$
This rule underlies the explicit schemes in the paper for generating a spatio-temporal optical vortex from a non-vortex Hermite-Gaussian-like packet by lens action, generating non-vortex fields with transverse OAM, and analyzing visible rotation of the intensity ellipse with principal-axis angle
$$
\theta=\frac{1}{2}\atan\!\left[\frac{2M_{12}}{M_{11}-M_{22}}\right].
$$
The authors further state that the same regular and unified formalism can be generalized to inhomogeneous and random media, and that partially coherent or vector spatio-temporal fields can be treated by matrix-valued Wigner functions and generalized moment matrices along the same lines [2404.06708].

## 3. Nonclassicality detection by Wigner-function moments

In continuous-variable quantum theory, Mallick et al. formulate a distinct **Wigner-Function-Moments** method in which the relevant quantities are not monomial phase-space moments but
$$
w_n=\iint W^n(x,p)\,dx\,dp,\qquad n=1,2,3,\dots
$$
with $w_1=1$ by normalization. The key statement is that if a state has a nonnegative Wigner function, then the first three such moments satisfy
$$
w_2^2\le w_3.
$$
Therefore,
$$
\Delta=w_2^2-w_3>0
$$
is a witness of Wigner negativity. The same condition can be written through a moment-matrix minor: positivity of $W$ forces
$$
\det\!\begin{pmatrix}1 & w_2\\ w_2 & w_3\end{pmatrix}=w_3-w_2^2
$$
to be nonnegative [2407.12116].

The paper gives explicit examples. For the single-photon Fock state,
$$
W_{|1\rangle}(x,p)=\frac{-1}{\pi}e^{-(x^2+p^2)}\bigl[1-2(x^2+p^2)\bigr],
$$
with
$$
w_2=0.159155,\qquad w_3\approx 0.00375,
$$
so $\Delta>0.02>0$. For the single-photon-subtracted squeezed-vacuum state, one finds
$$
w_2=\frac{1}{16\pi^4},\qquad w_3=\frac{1}{81\pi^4},
$$
and again $\Delta>0$, while the two-mode squeezed vacuum satisfies $w_2^2<w_3$, consistent with positivity [2407.12116].

A major operational advantage is that $w_2$ and $w_3$ can be measured without full tomography or Wigner-function reconstruction. With two copies, the expectation value of a continuous-variable SWAP test equals $w_2$; with three copies, the expectation of a cyclic permutation equals $w_3$. The proposal therefore replaces reconstruction of $W(x,p)$ by expectation values of permutation operators on a small number of copies [2407.12116].

This approach is mathematically distinct from the moment method of Bednorz and Belzig. There the witness is based on quantities of the form
$$
\langle f^2(q,p)\rangle_W=\int dq\,dp\,f^2(q,p)\,W(q,p),
$$
with $f$ a polynomial in $q$ and $p$. They show that no second-order choice of $f$ can reveal negativity, that fourth-order moments are necessary and sufficient in the generic non-rotationally-invariant case, and that rotationally invariant states require eighth-order moments. For the single-photon state, the polynomial
$$
f(q,p)=(q^2+p^2)^2-12(q^2+p^2)+26
$$
gives
$$
\langle f^2\rangle_W=-28<0.
$$
A common misconception is that the statement “fourth moments are necessary” conflicts with the criterion $w_2^2\le w_3$; it does not. The former concerns mixed moments of phase-space coordinates under $W$, whereas the latter concerns moments of the Wigner function itself [1103.1245][2407.12116].

## 4. Wigner D-moments and exact recovery on $SO(3)$

On the rotation group $SO(3)$, the relevant Wigner objects are the Wigner D-functions $D^\ell_{m,n}(R)$, which form a complete orthonormal basis of $L^2(SO(3))$. Filbir and Schröder study an unknown finite signed measure
$$
\mu=\sum_{k=1}^K a_k\,\delta_{R_k}
$$
and its moments up to degree $N$,
$$
y_{\ell,m,n}=\int_{SO(3)} D^\ell_{m,n}(R)\,d\mu(R)
=\sum_{k=1}^K a_k D^\ell_{m,n}(R_k).
$$
Recovery is posed as total-variation minimization over signed Borel measures with these moments fixed [1606.05306].

The main theorem states that if $N\ge 20$ and the support set obeys the minimal-separation condition
$$
\rho(C)\ge \frac{36}{N+1},
$$
then the true discrete measure is the unique solution of the total-variation problem. In that sense, low-degree Wigner D-moments suffice for exact recovery provided the support is sufficiently spread out [1606.05306].

The proof follows the now-standard dual-certificate strategy. One constructs a function
$$
q\in \mathrm{span}\{D^\ell_{m,n}:\ell\le N\}
$$
such that $q(R_k)=\mathrm{sign}(a_k)$ on the support and $|q(R)|<1$ off the support. The ansatz uses a zonal interpolation kernel
$$
\sigma_N(R,S)=\sum_{\ell=0}^N h_N(\ell)\sum_{m,n=-\ell}^{\ell}
D^\ell_{m,n}(R)\,\overline{D^\ell_{m,n}(S)}
=\sum_{\ell=0}^N h_N(\ell)\,U_{2\ell}(\cos(\tfrac12 d(S,R))),
$$
together with left-invariant differential operators and Hermite interpolation conditions. The technical core is a collection of localization estimates for $\sigma_N$ and its derivatives, with explicit constants, which yield invertibility of the interpolation system and strict off-support control [1606.05306].

This use of Wigner moments is structurally different from phase-space covariance methods or nonclassicality witnesses. Here the moments are harmonic coefficients on a compact Lie group, and their role is sparse measure recovery rather than diagnostics of a quantum state. The paper explicitly relates the result to exact super-resolution from low-pass data.

## 5. Moment hierarchies from Wigner transforms in kinetic and electromagnetic theory

In the Boltzmann setting studied by Chen, Denlinger and Pavlović, the Wigner transform is used to convert a phase-space density $f(t,x,v)$ into a density matrix $\gamma(x,x')$, so that the Boltzmann equation becomes a Schrödinger-type evolution
$$
(i\partial_t+\tfrac12(\Delta_x-\Delta_{x'}))\gamma=B(\gamma,\gamma).
$$
The analysis is carried out in weighted Sobolev spaces $H^{\alpha,\beta}$ that measure regularity in $x$ and decay in $v$ simultaneously. For $\alpha,\beta>(d-1)/2$ and bounded collision kernel, the authors prove local well-posedness, persistence of regularity, and continuity of the solution map. They also show propagation of spatial moments and velocity derivatives: if the initial data are sufficiently regular, then quantities such as
$$
\langle x\rangle^k\langle \nabla_v\rangle^\ell f
$$
remain in $L^\infty_TL^2_{x,v}$. The mechanism is an exchange of regularity in return for moments of the inverse Wigner transform, and the resulting regular solutions satisfy non-negativity, conservation of energy, and the $H$-theorem [1804.04019].

Perepelkin et al. develop a gauge-invariant Wigner formalism for a charged scalar particle in an electromagnetic field using the Weyl-Stratonovich transform. The Wigner function is defined with a phase factor containing the line integral of the vector potential, which guarantees gauge invariance. From the exact Wigner-Vlasov or Moyal equation, they derive moment equations for the density
$$
n(x,t)=\int f(x,P,t)\,d^3P,
$$
the momentum density
$$
j(x,t)=\int P\,f(x,P,t)\,d^3P,
$$
and the pressure tensor
$$
P_{ij}(x,t)=\int P_iP_j\,f\,d^3P.
$$
The exact hierarchy is infinite. Their Vlasov-Moyal approximation truncates the average acceleration field to leading order in $\hbar^2$, yielding closure at second order and recovering the classical Lorentz-force form at $\hbar^0$. Within that approximation, the Boltzmann $H$-functional is entropy-conserving because $\nabla_P\!\cdot a=0$ for the pure Lorentz force; retaining higher $\hbar^2$ terms gives nontrivial $H$-function evolution [2308.12613].

These two lines of work use Wigner moments in a PDE sense rather than as finite-dimensional descriptors alone. In one case the moments propagate regularity and decay properties of Boltzmann solutions; in the other they generate a gauge-invariant fluid hierarchy with a systematically improvable quantum correction structure.

## 6. Discrete moment recovery and numerical schemes

The lattice Wigner equation of Sol et al. realizes a numerical **Wigner-function-moments** method by discretizing momentum space into a finite set of lattice momenta $p_i$ with positive weights $w_i$ chosen to reproduce exactly all moments up to order $N_\Pi$. The continuous moment
$$
M^{(n)}(x,t)=\int_{-\infty}^{\infty} p^n W(x,p,t)\,dp
$$
is approximated by
$$
M^{(n)}(x,t)\approx \sum_{i=1}^{N_q} w_i\,p_i^n\,f_i(x,t),
$$
with $f_i(x,t)=W(x,p_i,t)$. The method is designed so that the moments of the Wigner function are recovered exactly up to the desired order determined by the number of discrete momenta retained [1709.05934].

The discrete evolution gives, for all $n\le N_\Pi$,
$$
\partial_t[M^{(n)}]+\partial_x[M^{(n+1)}]
=\sum_{i=1}^{N_q} v_i^n S_i + O(\delta t^2),
$$
and numerical stability is enforced by an artificial BGK-type collision operator
$$
\Omega_i=-\frac{1}{\tau_w}\bigl[\bar W_i-\bar W_i^{\mathrm{eq}}\bigr].
$$
The equilibrium populations are constructed so as to preserve exactly the first $N_\Pi$ moments:
$$
\sum_i v_i^n\bar W_i^{\mathrm{eq}}=\sum_i v_i^n\bar W_i,\qquad n=0,1,\dots,N_\Pi.
$$
Thus the collision operator damps only modes of order larger than $N_\Pi$ and leaves the relevant moment dynamics exact up to $O(\delta t^2)$ [1709.05934].

The paper validates the scheme for harmonic and anharmonic oscillators and applies it to one- and two-dimensional open systems with potential barriers. For the harmonic oscillator, the root-mean-square density error after one period satisfies $\Delta\propto (\delta x)^2$, and the method reproduces oscillations of $\rho(x,t)=M^{(0)}$ and $j(x,t)=M^{(1)}$ to less than $10^{-3}$ accuracy. The authors also illustrate computational viability for three-dimensional open systems with a D3Q125 lattice and discuss the cost scaling in dimensions $d=1,2,3$ [1709.05934].

Taken together, these developments show that the Wigner Function Moments method is best understood as a methodological class rather than a single algorithm. In optics it yields a compact algebra for size, spread, correlation, transverse OAM, and ABCD propagation [2404.06708]. In quantum diagnostics it provides experimentally accessible witnesses of Wigner negativity without full tomography [2407.12116]. On $SO(3)$ it becomes a sparse recovery problem for discrete measures [1606.05306]. In kinetic and electromagnetic theory it supports regularity propagation, hierarchy closure, and $H$-function analysis [1804.04019][2308.12613]. In lattice formulations it becomes a moment-preserving discretization principle that ties numerical accuracy directly to the retained moment order [1709.05934].

Source: https://www.emergentmind.com/topics/wigner-function-moments-method