Wigner Function Moments Method

• Wigner Function Moments method is a family of phase-space techniques that extract information via moment hierarchies rather than full distribution reconstructions.
• It unifies diverse applications by encoding spatial-angular optics, quantum negativity detection, harmonic recovery on SO(3), and numerical discretizations in kinetic and electromagnetic theories.
• The approach enables efficient data compression and measurement through covariance matrices, permutation tests, and moment-preserving schemes, offering practical insights across physics and computation.

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 (Bekshaev et al., 2024, Mallick et al., 2024, Filbir et al., 2016, Chen et al., 2018, Perepelkin et al., 2023).

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 (Bekshaev et al., 2024). 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 (Mallick et al., 2024). 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 (Filbir et al., 2016).

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	$I(x,s;k_x,\Delta k)$0	Exact recovery of discrete measures
Kinetic and lattice Wigner theory	$I(x,s;k_x,\Delta k)$1	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 $I(x,s;k_x,\Delta k)$2 as

$I(x,s;k_x,\Delta k)$3

with marginals

$I(x,s;k_x,\Delta k)$4

Bekshaev, Angelsky and Zenkova generalize this construction to spatio-temporal light fields by introducing the retarded coordinate $I(x,s;k_x,\Delta k)$5 and the slowly-varying envelope $I(x,s;k_x,\Delta k)$6. In $I(x,s;k_x,\Delta k)$7 dimensions, the corresponding $I(x,s;k_x,\Delta k)$8D Wigner function is

$I(x,s;k_x,\Delta k)$9

and the total action is

$M_{k,\ell}=\int\!\!\int x^k p^\ell W(x,p)\,dx\,dp$0

This construction reduces to the standard spatial Wigner function if the $M_{k,\ell}=\int\!\!\int x^k p^\ell W(x,p)\,dx\,dp$1-dependence is ignored and to the time-frequency Wigner function if $M_{k,\ell}=\int\!\!\int x^k p^\ell W(x,p)\,dx\,dp$2 is ignored (Bekshaev et al., 2024).

The first and second moments are organized through the $M_{k,\ell}=\int\!\!\int x^k p^\ell W(x,p)\,dx\,dp$3D “ray” vector

$M_{k,\ell}=\int\!\!\int x^k p^\ell W(x,p)\,dx\,dp$4

with centroid and moment matrix

$M_{k,\ell}=\int\!\!\int x^k p^\ell W(x,p)\,dx\,dp$5

Assuming a centered packet, the real-space expressions include

$M_{k,\ell}=\int\!\!\int x^k p^\ell W(x,p)\,dx\,dp$6

$M_{k,\ell}=\int\!\!\int x^k p^\ell W(x,p)\,dx\,dp$7

and mixed terms such as

$M_{k,\ell}=\int\!\!\int x^k p^\ell W(x,p)\,dx\,dp$8

The diagonal entries characterize spatial size, temporal length, angular spread, and frequency spread, while the off-diagonal blocks $M_{k,\ell}=\int\!\!\int x^k p^\ell W(x,p)\,dx\,dp$9 encode mixed spatio-temporal correlations (Bekshaev et al., 2024).

A central physical application is the description of transverse orbital angular momentum in spatio-temporal optical vortices. The transverse OAM about the $w_n=\iint W^n(x,p)\,dx\,dp,$0 axis is written as

$w_n=\iint W^n(x,p)\,dx\,dp,$1

and in practice the $w_n=\iint W^n(x,p)\,dx\,dp,$2 term often vanishes by centroid definition, so that $w_n=\iint W^n(x,p)\,dx\,dp,$3. If the field carries a phase singularity in the $w_n=\iint W^n(x,p)\,dx\,dp,$4 plane, then $w_n=\iint W^n(x,p)\,dx\,dp,$5 and hence $w_n=\iint W^n(x,p)\,dx\,dp,$6. 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 (Bekshaev et al., 2024).

The transformation law is especially simple for first-order paraxial systems. If

$w_n=\iint W^n(x,p)\,dx\,dp,$7

with $w_n=\iint W^n(x,p)\,dx\,dp,$8 built from the usual $w_n=\iint W^n(x,p)\,dx\,dp,$9 ABCD matrix acting on the $SO(3)$0 subspace and the $SO(3)$1 subspace left untouched, then

$SO(3)$2

In particular,

$SO(3)$3

and the transverse-OAM entry obeys

$SO(3)$4

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

$SO(3)$5

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 (Bekshaev et al., 2024).

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

$SO(3)$6

with $SO(3)$7 by normalization. The key statement is that if a state has a nonnegative Wigner function, then the first three such moments satisfy

$SO(3)$8

Therefore,

$SO(3)$9

is a witness of Wigner negativity. The same condition can be written through a moment-matrix minor: positivity of $$y_{\ell,m,n}=\int_{SO(3)} D^\ell_{m,n}(R)\,d\mu(R),$$0 forces

$$y_{\ell,m,n}=\int_{SO(3)} D^\ell_{m,n}(R)\,d\mu(R),$$1

to be nonnegative (Mallick et al., 2024).

The paper gives explicit examples. For the single-photon Fock state,

$$y_{\ell,m,n}=\int_{SO(3)} D^\ell_{m,n}(R)\,d\mu(R),$$2

with

$$y_{\ell,m,n}=\int_{SO(3)} D^\ell_{m,n}(R)\,d\mu(R),$$3

so $$y_{\ell,m,n}=\int_{SO(3)} D^\ell_{m,n}(R)\,d\mu(R),$$4. For the single-photon-subtracted squeezed-vacuum state, one finds

$$y_{\ell,m,n}=\int_{SO(3)} D^\ell_{m,n}(R)\,d\mu(R),$$5

and again $$y_{\ell,m,n}=\int_{SO(3)} D^\ell_{m,n}(R)\,d\mu(R),$$6, while the two-mode squeezed vacuum satisfies $$y_{\ell,m,n}=\int_{SO(3)} D^\ell_{m,n}(R)\,d\mu(R),$$7, consistent with positivity (Mallick et al., 2024).

A major operational advantage is that $$y_{\ell,m,n}=\int_{SO(3)} D^\ell_{m,n}(R)\,d\mu(R),$$8 and $$y_{\ell,m,n}=\int_{SO(3)} D^\ell_{m,n}(R)\,d\mu(R),$$9 can be measured without full tomography or Wigner-function reconstruction. With two copies, the expectation value of a continuous-variable SWAP test equals $$M_{ij}=(1/F)\int (P_i-\bar P_i)(P_j-\bar P_j)I(P)\,dP$$0; with three copies, the expectation of a cyclic permutation equals $$M_{ij}=(1/F)\int (P_i-\bar P_i)(P_j-\bar P_j)I(P)\,dP$$1. The proposal therefore replaces reconstruction of $$M_{ij}=(1/F)\int (P_i-\bar P_i)(P_j-\bar P_j)I(P)\,dP$$2 by expectation values of permutation operators on a small number of copies (Mallick et al., 2024).

This approach is mathematically distinct from the moment method of Bednorz and Belzig. There the witness is based on quantities of the form

$$M_{ij}=(1/F)\int (P_i-\bar P_i)(P_j-\bar P_j)I(P)\,dP$$3

with $$M_{ij}=(1/F)\int (P_i-\bar P_i)(P_j-\bar P_j)I(P)\,dP$$4 a polynomial in $$M_{ij}=(1/F)\int (P_i-\bar P_i)(P_j-\bar P_j)I(P)\,dP$$5 and $$M_{ij}=(1/F)\int (P_i-\bar P_i)(P_j-\bar P_j)I(P)\,dP$$6. They show that no second-order choice of $$M_{ij}=(1/F)\int (P_i-\bar P_i)(P_j-\bar P_j)I(P)\,dP$$7 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

$$M_{ij}=(1/F)\int (P_i-\bar P_i)(P_j-\bar P_j)I(P)\,dP$$8

gives

$$M_{ij}=(1/F)\int (P_i-\bar P_i)(P_j-\bar P_j)I(P)\,dP$$9

A common misconception is that the statement “fourth moments are necessary” conflicts with the criterion $w_n=\iint W^n\,dx\,dp$0; it does not. The former concerns mixed moments of phase-space coordinates under $w_n=\iint W^n\,dx\,dp$1, whereas the latter concerns moments of the Wigner function itself (Bednorz et al., 2011, Mallick et al., 2024).

4. Wigner D-moments and exact recovery on $w_n=\iint W^n\,dx\,dp$2

On the rotation group $w_n=\iint W^n\,dx\,dp$3, the relevant Wigner objects are the Wigner D-functions $w_n=\iint W^n\,dx\,dp$4, which form a complete orthonormal basis of $w_n=\iint W^n\,dx\,dp$5. Filbir and Schröder study an unknown finite signed measure

$w_n=\iint W^n\,dx\,dp$6

and its moments up to degree $w_n=\iint W^n\,dx\,dp$7,

$w_n=\iint W^n\,dx\,dp$8

Recovery is posed as total-variation minimization over signed Borel measures with these moments fixed (Filbir et al., 2016).

The main theorem states that if $w_n=\iint W^n\,dx\,dp$9 and the support set obeys the minimal-separation condition

$M_{k,\ell}=\iint x^k p^\ell W\,dx\,dp$0

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 (Filbir et al., 2016).

The proof follows the now-standard dual-certificate strategy. One constructs a function

$M_{k,\ell}=\iint x^k p^\ell W\,dx\,dp$1

such that $M_{k,\ell}=\iint x^k p^\ell W\,dx\,dp$2 on the support and $M_{k,\ell}=\iint x^k p^\ell W\,dx\,dp$3 off the support. The ansatz uses a zonal interpolation kernel

$M_{k,\ell}=\iint x^k p^\ell W\,dx\,dp$4

together with left-invariant differential operators and Hermite interpolation conditions. The technical core is a collection of localization estimates for $M_{k,\ell}=\iint x^k p^\ell W\,dx\,dp$5 and its derivatives, with explicit constants, which yield invertibility of the interpolation system and strict off-support control (Filbir et al., 2016).

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 $M_{k,\ell}=\iint x^k p^\ell W\,dx\,dp$6 into a density matrix $M_{k,\ell}=\iint x^k p^\ell W\,dx\,dp$7, so that the Boltzmann equation becomes a Schrödinger-type evolution

$M_{k,\ell}=\iint x^k p^\ell W\,dx\,dp$8

The analysis is carried out in weighted Sobolev spaces $M_{k,\ell}=\iint x^k p^\ell W\,dx\,dp$9 that measure regularity in $SO(3)$0 and decay in $SO(3)$1 simultaneously. For $SO(3)$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

$SO(3)$3

remain in $SO(3)$4. 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 $SO(3)$5-theorem (Chen et al., 2018).

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

$SO(3)$6

the momentum density

$SO(3)$7

and the pressure tensor

$SO(3)$8

The exact hierarchy is infinite. Their Vlasov-Moyal approximation truncates the average acceleration field to leading order in $SO(3)$9, yielding closure at second order and recovering the classical Lorentz-force form at $I(x,s;k_x,\Delta k)$00. Within that approximation, the Boltzmann $I(x,s;k_x,\Delta k)$01-functional is entropy-conserving because $I(x,s;k_x,\Delta k)$02 for the pure Lorentz force; retaining higher $I(x,s;k_x,\Delta k)$03 terms gives nontrivial $I(x,s;k_x,\Delta k)$04-function evolution (Perepelkin et al., 2023).

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 $I(x,s;k_x,\Delta k)$05 with positive weights $I(x,s;k_x,\Delta k)$06 chosen to reproduce exactly all moments up to order $I(x,s;k_x,\Delta k)$07. The continuous moment

$I(x,s;k_x,\Delta k)$08

is approximated by

$I(x,s;k_x,\Delta k)$09

with $I(x,s;k_x,\Delta k)$10. 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 (Solorzano et al., 2017).

The discrete evolution gives, for all $I(x,s;k_x,\Delta k)$11,

$I(x,s;k_x,\Delta k)$12

and numerical stability is enforced by an artificial BGK-type collision operator

$I(x,s;k_x,\Delta k)$13

The equilibrium populations are constructed so as to preserve exactly the first $I(x,s;k_x,\Delta k)$14 moments:

$I(x,s;k_x,\Delta k)$15

Thus the collision operator damps only modes of order larger than $I(x,s;k_x,\Delta k)$16 and leaves the relevant moment dynamics exact up to $I(x,s;k_x,\Delta k)$17 (Solorzano et al., 2017).

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 $I(x,s;k_x,\Delta k)$18, and the method reproduces oscillations of $I(x,s;k_x,\Delta k)$19 and $I(x,s;k_x,\Delta k)$20 to less than $I(x,s;k_x,\Delta k)$21 accuracy. The authors also illustrate computational viability for three-dimensional open systems with a D3Q125 lattice and discuss the cost scaling in dimensions $I(x,s;k_x,\Delta k)$22 (Solorzano et al., 2017).

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 (Bekshaev et al., 2024). In quantum diagnostics it provides experimentally accessible witnesses of Wigner negativity without full tomography (Mallick et al., 2024). On $I(x,s;k_x,\Delta k)$23 it becomes a sparse recovery problem for discrete measures (Filbir et al., 2016). In kinetic and electromagnetic theory it supports regularity propagation, hierarchy closure, and $I(x,s;k_x,\Delta k)$24-function analysis (Chen et al., 2018, Perepelkin et al., 2023). In lattice formulations it becomes a moment-preserving discretization principle that ties numerical accuracy directly to the retained moment order (Solorzano et al., 2017).