---
title: 'Spectral Sample: Theory and Applications'
url: https://www.emergentmind.com/topics/spectral-sample
type: topic
---

# Spectral Sample: Theory and Applications

“Spectral sample” is a field-dependent term used for constructions in which sampling is defined through spectral data rather than by direct coordinate selection. In algebraic-combinatorial statistics, it denotes the eigenvalue vector of a principal submatrix of a Haar-randomly conjugated diagonal matrix; in continuum percolation, it denotes a random finite subset of the plane whose law is given by squared Fourier coefficients averaged over random geometry; in graph signal processing, it denotes down- and upsampling performed directly on graph Fourier coefficients; and in optical turbulence simulation, it denotes randomized sampling of spatial frequencies in Fourier space to improve low-frequency phase statistics [1302.5892] [1906.12256] [1706.05147] [1905.07074].

## 1. Spectral sampling as an eigenvalue analogue of simple random sampling

Di Nardo, McCullagh, and Senato define spectral sampling from a finite population of real numbers
\[
x=(x_1,\dots,x_n), \qquad X=\operatorname{diag}(x),
\]
by first applying a Haar-distributed unitary matrix \(H\), so that \(HXH^\dagger\) is freely randomized, and then taking the leading \(m\times m\) principal submatrix
\[
Y=(HXH^\dagger)_{[m\times m]}.
\]
The spectral sample of size \(m\le n\) is the eigenvalue vector
\[
y=(y_1,\dots,y_m).
\]
For \(m=n\), \(y\) is a random permutation of \(x\); for \(m<n\), the \(y_i\) are typically nonatomic and do not coincide with the original \(x_i\) [1302.5892].

This construction is presented as the matrix or eigenvalue analogue of ordinary simple random sampling. If the unitary group is replaced by the permutation group \({\mathfrak S}_n\), then \(y\) becomes an ordinary simple random sample from \(x\). The associated notion of a natural statistic is likewise transferred from the classical setting: a statistic \(T=\{T_n\}\) is natural relative to spectral sampling if
\[
E\!\left(T_m(y)\mid x\right)=T_n(x),\qquad m\le n.
\]
In this sense, the statistic is inherited on the average. The paper places this viewpoint alongside Tukey’s classical natural symmetric polynomial statistics, including the normalized augmented symmetric functions
\[
\tilde{\mathfrak a}_{r,n}(x)=\frac1n\sum_i x_i^r,\qquad 
\tilde{\mathfrak a}_{rs,n}(x)=\frac{1}{(n)_2}\sum_{i\ne j}x_i^r x_j^s,
\]
with \((n)_r=n(n-1)\cdots(n-r+1)\) [1302.5892].

## 2. Natural statistics, spectral \(k\)-statistics, and free-cumulant asymptotics

The same paper develops the corresponding analogue of Fisher’s \(k\)-statistics and Tukey’s polykays. The central observables are trace powers of a random Hermitian matrix \(M\),
\[
\operatorname{Tr}(M^k)=\nu_1^k+\cdots+\nu_m^k,
\]
where \((\nu_1,\dots,\nu_m)\) are the eigenvalues of \(M\). Using umbral calculus, cumulants \(\{c_i\}\) are defined from moments \(\{a_i\}\) by
\[
1+\sum_{i\ge1} a_i\frac{z^i}{i!} =\exp\!\left\{\sum_{i\ge1} c_i\frac{z^i}{i!}\right\},
\]
and the cumulant umbra satisfies
\[
f(\kappa_\alpha,z)=1+\log f(\alpha,z), \qquad \kappa_\alpha \equiv \chi.\alpha.
\]
This symbolic machinery is extended to multivariate sequences and to eigenvalue data [1302.5892].

For a permutation \(\sigma\in{\mathfrak S}_i\), the paper defines
\[
\mu(Y)(\sigma)=\prod_{c\in C(\sigma)} \operatorname{Tr}\!\bigl(Y^{\ell(c)}\bigr),
\]
where \(C(\sigma)\) is the set of cycles of \(\sigma\) and \(\ell(c)\) is the cycle length. The inversion
\[
\tilde{\kappa}(y)=\mu(I_m)^{(-1)}\cdot \mu(Y)
\]
produces the spectral analogue of the classical moment–cumulant relation. The resulting matricial polykays, or spectral \(k\)-statistics, are unbiased estimators of products of cumulants of trace powers. For a partition \(\lambda=(1^{r_1},2^{r_2},\dots)\vdash i\), the corresponding quantity \(\kappa_\lambda(y)\) is constructed so that its expectation equals a product of cumulants of the trace process [1302.5892].

The asymptotic normalization connects spectral sampling to free probability. The normalized spectral \(k\)-statistics are defined by
\[
\tilde{\kappa}^{(N)}_\lambda(y):=m^{\,i-l(\lambda)}\tilde{\kappa}_\lambda(y),
\]
and satisfy
\[
E\!\left[\tilde{\kappa}^{(N)}_\lambda(y)\right] \longrightarrow c_1^{r_1}c_2^{r_2}\cdots,
\qquad \lambda=(1^{r_1},2^{r_2},\dots)\vdash i,
\]
where \(c_j\) are the free cumulants of the limiting noncommutative distribution. Generalized spectral polykays are then obtained by Möbius inversion on set partitions,
\[
l_\pi(y)=\sum_{\tau\ge \pi} {\mathfrak m}(\pi,\tau)\,\tilde{\kappa}_\tau(y),
\]
with normalized version
\[
l_\pi^{(N)}(y)=m^{i-|\pi|}l_\pi(y).
\]
A plausible implication is that spectral sampling provides a finite-\(m\) combinatorial bridge between symmetric-function statistics and free-cumulant asymptotics [1302.5892].

## 3. Annealed spectral sample in Voronoi percolation

In critical Voronoi percolation, the term refers to a different object. Given a bounded measurable Boolean function \(h:\Omega\to\mathbb R\) on colored Poisson–Voronoi configurations \(\omega\), and a fixed realization \(\eta\) of the underlying point process, the quenched Fourier expansion is
\[
h^\eta=\sum_{S\subseteq_f \eta}\widehat{h^\eta}(S)\,\chi_S^\eta,
\qquad 
\chi_S^\eta(\omega)=\prod_{x\in S}\omega(x).
\]
Rather than keeping the quenched spectral measure \(\widehat Q_{h^\eta}[S]=\widehat{h^\eta}(S)^2\), the paper defines the annealed version by averaging over \(\eta\):
\[
\widehat Q_h^{an}[A] = E\!\left[\sum_{S\subseteq_f \eta,\;S\in A}\widehat{h^\eta}(S)^2\right],
\]
and the normalized law
\[
\widehat P_h^{an}[A] = \frac{\widehat Q_h^{an}[A]}{E[h^2]}.
\]
A random variable with law \(\widehat P_h^{an}\) is the annealed spectral sample \(\mathcal S_h^{an}\). It is a random finite subset of \(\mathbb R^2\), i.e. a continuous point process [1906.12256].

A crucial identity is the covariance formula for frozen dynamical Voronoi percolation,
\[
E\big[h(\omega^{froz}(0))\,h(\omega^{froz}(t))\big]
=
\sum_{k\in\mathbb N}\widehat Q_h^{an}[|S|=k]\,e^{-kt},
\]
which diagonalizes temporal decorrelation in terms of spectral-sample cardinality. For the left-right crossing event \(g_n\) of \([0,n]^2\), the paper proves
\[
P\!\left[0<|\mathcal S_{g_n}^{an}|<r^2\alpha_4^{an}(r)\right]
\le
C\left(\frac nr\,\alpha_4^{an}(r,n)\right)^2+\frac{C}{n},
\qquad 1\le r\le n,
\]
and
\[
E|\mathcal S_{g_n}^{an}|\asymp n^2\alpha_4^{an}(n).
\]
For the one-arm event \(f_R=1_{\{0\leftrightarrow \partial B(0,R)\}}\), it proves
\[
P\!\left[0<|\mathcal S_{f_R}^{an}|\le r^2\alpha_4^{an}(r)\right]
\le C\,\alpha_1^{an}(r,R),
\]
together with a clustering estimate showing that if the spectral mass is small, it tends to be spatially concentrated near the origin [1906.12256].

The paper’s application is dynamical critical Voronoi percolation. The lower-tail and clustering bounds are used to show that there exist almost surely exceptional times with an unbounded monochromatic component, both for frozen color-resampling dynamics and for dynamics in which Voronoi centers move according to long-range stable Lévy processes [1906.12256].

## 4. Spectral-domain sampling of graph signals

In graph signal processing, spectral sampling is defined directly in the graph Fourier domain rather than in the vertex domain. For a graph Laplacian
\[
\mathbf{L}=\mathbf{D}-\mathbf{A},\qquad
\mathbf{L}=\mathbf{U}\mathbf{\Lambda}\mathbf{U}^*,
\]
the graph Fourier transform of a signal \(\mathbf f\) is
\[
\widetilde{f}[i]=\langle \mathbf{u}_i,\mathbf{f}\rangle
=\sum_{n=0}^{N-1}u_i^*[n]f[n].
\]
The paper argues that conventional vertex-domain sampling preserves retained values, but does not in general produce the classical frequency-domain effects of bandwidth expansion, aliasing, compression, or imaging on graphs [1706.05147].

For downsampling by a factor \(M\), the index-based spectral downsampling operator (GD2) is
\[
\widetilde{f}_{d}[k] = \sum_{p=0}^{M-1}\widetilde{f}\!\left[\frac{pN}{M}+k\right],
\]
or, in matrix form,
\[
\mathbf{f}_d=\mathbf{U}_1\mathbf{S}_d\mathbf{U}_0^*\mathbf{f},
\qquad
\mathbf{S}_d=
\begin{bmatrix}
\mathbf{I}_{N/M} & \mathbf{I}_{N/M} & \cdots
\end{bmatrix}.
\]
Because graph eigenvalues are generally nonuniformly spaced, the paper also defines the spectrum-based version (GD3),
\[
\widetilde{f}_{d}[k] = \sum_{p=0}^{M-1} 
\widetilde{f}_{\mathrm{int}\!\left(\frac{\rho}{M}\left(\lambda_{1,k}+p\lambda_{1,\max}\right)\right)},
\qquad
\rho=\frac{\lambda_{0,\max}}{\lambda_{1,\max}},
\]
together with symmetrized variants \((GD2′)\) and \((GD3′)\) designed to better mimic classical aliasing behavior [1706.05147].

Upsampling is defined analogously through spectral repetition. The index-based operator (GU2) repeats the original spectrum \(L\) times,
\[
\widetilde{f}_{u}[pN+k]=\widetilde{f}[k], \qquad p=0,\dots,L-1,
\]
and the spectrum-based operator (GU3) interpolates the repeated spectrum according to the enlarged graph. The key point is that these operators preserve frequency-domain behavior rather than vertex-domain identity. The paper shows that on ring graphs they reduce to classical time-domain sampling, while on path graphs and more general graphs the graph-dependent eigenstructure materially changes the sampling effect. Fractional sampling and Laplacian pyramid representation are presented as applications [1706.05147].

## 5. Randomized spectral sampling in optical turbulence simulation

In optical turbulence, randomized spectral sampling is an FFT-based phase-screen synthesis method that randomizes the sampled spatial frequencies instead of evaluating spectral coefficients only on the rigid DFT lattice. The standard synthesis uses
\[
\theta(j,l) = \sum_{n,m=-M/2}^{M/2-1} \tilde c(n\Delta\kappa_x,m\Delta\kappa_y)
\,\exp\!\left[2\pi i\frac{jn+lm}{M}\right],
\]
with
\[
\tilde c(\beta,\gamma)=(a+ib)\,k\,\sqrt{2\pi \Delta z\,\Delta\kappa_x\,\Delta\kappa_y\,\Phi_n(\beta \hat e_x+\gamma \hat e_y)}.
\]
The randomized method draws
\[
\delta \kappa_x \sim U\!\left[-\frac{\Delta \kappa_x}{2},\frac{\Delta \kappa_x}{2}\right],
\qquad
\delta \kappa_y \sim U\!\left[-\frac{\Delta \kappa_y}{2},\frac{\Delta \kappa_y}{2}\right],
\]
and replaces the fixed lattice by a randomly shifted one [1905.07074].

Using the Fourier transform shift theorem, the inverse transform can be implemented as a standard inverse FFT followed by a phase ramp,
\[
C(j,l) = M^2 \cdot \mathcal{F}^{-1}_2[\tilde C(n,m)],
\]
\[
\theta_R(j,l)=\exp\!\big[i(j\Delta x\,\delta\kappa_x+l\Delta y\,\delta\kappa_y)\big]\; C(j,l).
\]
Equivalently,
\[
\theta_R(j,l)= \sum_{n,m=-M/2}^{M/2-1} 
\tilde c\!\left(n\Delta\kappa_x+\delta\kappa_x,\;m\Delta\kappa_y+\delta\kappa_y\right)
\exp\!\left[i\big(j\Delta x(n\Delta\kappa_x+\delta\kappa_x)+l\Delta y(m\Delta\kappa_y+\delta\kappa_y)\big)\right].
\]
This produces phase screens that are no longer strictly periodic and that contain domain-wide low spatial frequency distortions, including tip/tilt-like structure, without introducing a separate subharmonic subgrid [1905.07074].

The principal validation metric is the phase structure function. For isotropic turbulence,
\[
D_\theta(\rho)=8\pi^2 k^2 \Delta z \int_0^\infty \kappa_\rho \Phi_n(\kappa_\rho)\big[1-J_0(\rho\kappa_\rho)\big]\,d\kappa_\rho,
\]
and the paper compares simulation to theory using RMS relative error
\[
\mathcal E=100\%\times \frac{\mathcal E_x+\mathcal E_y}{2}.
\]
For bounded spectra with finite outer scale, it reports that for
\[
1 \le L_0/(M\Delta x)\le 10^3,
\]
the randomized method kept the half-domain RMS error between about \(0.34\%\) and \(3.79\%\), whereas the traditional unrandomized FFT grid ranged up to about \(61.5\%\). For generalized non-Kolmogorov power laws with \(3<\alpha<4\), the paper proposes a hybrid method
\[
\theta_H(j,l)=\theta_R(j,l)+\sum_{p=1}^{N_p}\theta_{out}(j,l,p)+\theta_{in}(j,l),
\]
supplemented, when needed, by white-noise terms computed from unresolved spectral energy outside the sampled region [1905.07074].

## 6. Terminological scope, later extensions, and recurrent themes

The cited literature does not use “spectral sample” for a single universal mathematical object. Instead, the term is attached to several constructions whose common feature is that the sampling operation is defined on eigenvalues, Fourier coefficients, or spectral weights rather than on raw coordinates, vertices, or pixels.

| Setting | Spectral sample object | Defining mechanism |
|---|---|---|
| Finite-population matrix statistics | Eigenvalue vector \(y\) | Principal submatrix of \(HXH^\dagger\) |
| Voronoi percolation | Random finite subset \(\mathcal S_h^{an}\subset\mathbb R^2\) | Squared Fourier coefficients averaged over \(\eta\) |
| Graph signal processing | Down-/upsampled graph spectrum | Sampling directly on graph Fourier coefficients |
| Optical turbulence | Randomized Fourier-grid phase screen | Cell-wise random frequency shifts |

A common misconception is to treat these usages as interchangeable. They are not. The spectral sample of Di Nardo, McCullagh, and Senato is an eigenvalue sample associated with unitary conjugation; the annealed spectral sample of Voronoi percolation is a Fourier-analytic random point process; the graph-signal and turbulence literatures use spectral sampling to mean direct manipulation of spectral coefficients or sampled frequencies [1302.5892] [1906.12256] [1706.05147] [1905.07074].

Later work extends the phrase in yet other directions. In finite-sample spectral graph inference, one observed graph is treated as a spectral sample: a single adjacency matrix \(A\) obtained through a finite random sampling mechanism, with confidence regions constructed for latent eigenspaces and propagated to clustering, centrality, fairness-constrained post-processing, and topological summaries [2602.10566]. In large-language-model fine-tuning, GradSentry scores each training example independently by the normalized spectral entropy of its per-sample gradient,
\[
\bar{H}(G_i^\prime) = \frac{H(G_i^\prime)}{\log k},
\]
where
\[
H(G_i^\prime) = -\sum_{j=1}^{k} p_j \log p_j.
\]
This is a per-sample spectral diagnostic rather than a classical spectral sample, but it preserves the same underlying idea that sample-level structure can be exposed through the spectrum of an associated operator [2605.26574].

Taken together, these usages suggest a stable conceptual core: a spectral sample is a sample defined, weighted, or analyzed through spectral information, with the exact object determined by the ambient theory.

Source: https://www.emergentmind.com/topics/spectral-sample