---
title: 'Spectral Born Machines: Quantum Fourier Models'
url: https://www.emergentmind.com/topics/spectral-born-machines
type: topic
---

# Spectral Born Machines: Quantum Fourier Models

Spectral Born Machines (SBMs) are a class of quantum generative models derived by viewing IQP Born Machines through the lens of group Fourier analysis and then generalizing from bitstrings to integer-structured data over $\mathbb{Z}_d$. They define a model distribution by the Born rule, use a “Fourier sandwich” circuit $U(\theta) = (F^{\otimes n})^\dagger D(\theta) F^{\otimes n}$, and train through Fourier- or graph-spectral maximum mean discrepancy objectives that are classically estimable. In the binary setting, the same perspective identifies Quantum Circuit Born Machines (QCBMs) as quantum Fourier models whose output distributions are exactly Walsh–Hadamard transforms of correlators, making the IQP Born Machine the $d=2$ special case of the broader SBM framework [2607.06675, 2511.01845].

## 1. Formal definition and circuit structure

Given a parameterized quantum circuit $U(\theta)$ acting on an initial computational basis state $|0\rangle$, the output probability of integer-valued datum $x$ is
$$
p(x;\theta) = |\langle x|U(\theta)|0\rangle|^2,
$$
i.e., the Born rule. SBMs define a model distribution $q_\theta(x)$ as the computational-basis measurement distribution of the evolved state,
$$
q_\theta(x) = \operatorname{Tr}\!\left[(|x\rangle\langle x| \otimes I)\rho_\theta\right], \qquad \rho_\theta = U(\theta)|0\rangle\langle 0|U(\theta)^\dagger.
$$
The canonical SBM circuit is a “Fourier sandwich” over $\mathbb{Z}_d^n$,
$$
U(\theta) = (F^{\otimes n})^\dagger D(\theta) F^{\otimes n},
$$
where $F$ is the quantum Fourier transform over $\mathbb{Z}_d$ and $D(\theta)$ is a diagonal phase unitary in the computational basis. For $\mathbb{Z}_N$,
$$
F_N|x\rangle = \frac{1}{\sqrt{N}}\sum_{k=0}^{N-1} e^{2\pi i x k/N}|k\rangle,
$$
and for $\mathbb{Z}_d^n$,
$$
F^{\otimes n}|x\rangle = \frac{1}{\sqrt{d^n}}\sum_{k\in \mathbb{Z}_d^n}\omega^{k\cdot x}|k\rangle, \qquad \omega = e^{2\pi i/d}.
$$
The Fourier basis is spanned by characters $\chi_k(x)=\omega^{k\cdot x}$ [2607.06675].

In the computational basis, the diagonal layer satisfies
$$
D(\theta)|z\rangle = e^{i\Phi_\theta(z)}|z\rangle,
$$
with phase function $\Phi_\theta:\mathbb{Z}_d^n\to\mathbb{R}$ computable in $\operatorname{poly}(n,d)$. A particularly useful class is generated by Heisenberg–Weyl observables,
$$
D(\theta)=\prod_{g\in\mathcal{G}} e^{i\theta_g \mathcal{Q}_g}, \qquad \mathcal{Q}_g=\mathcal{Q}(g,0),
$$
where each generator contributes a separable phase eigenvalue
$$
\mathcal{Q}_g|z\rangle = \phi_g(z)|z\rangle, \qquad \phi_g(z)=\prod_{j=1}^n \sqrt{2}\cos(2\pi z_j g_j/d+\pi/4),
$$
so that
$$
\Phi_\theta(z)=\sum_{g\in\mathcal{G}}\theta_g \phi_g(z).
$$
Restricting $\mathcal{G}$ by weight and degree yields strong spectral bias and reduces parameter count. The binary limit $d=2$ recovers the IQP Born Machine, because SBMs reduce to IQP BM when $d=2$ and $F=H$ [2607.06675].

The same structure admits an equivalent Walsh–Hadamard description for qubit QCBMs. If $p(x;\theta)$ is the computational-basis measurement distribution, then
$$
\hat p(s;\theta)=\sum_{x\in\{0,1\}^n} p(x)\,(-1)^{s\cdot x}=\langle Z^s\rangle_{\rho(\theta)},
$$
and
$$
p(x;\theta)=2^{-n}\sum_{s\in\{0,1\}^n}\langle Z^s\rangle_{\rho(\theta)}\,(-1)^{s\cdot x}.
$$
Thus QCBMs are spectral models whose parameters control a trigonometric–Pauli polynomial whose $Z$-string coefficients are directly the Fourier/Walsh–Hadamard components of the output distribution [2511.01845].

## 2. Spectral representation and MMD training

The defining training perspective of SBMs is spectral. For distributions $P,Q$ on $\{0,1\}^n$ and a translation-invariant kernel $k(x,y)=\kappa(x\oplus y)$ with Walsh–Hadamard decomposition
$$
\kappa(z)=\sum_{s\in\{0,1\}^n}\lambda_s \chi_s(z),
$$
the maximum mean discrepancy diagonalizes as
$$
\operatorname{MMD}^2(P,Q;k)=\sum_{s\in\{0,1\}^n}\lambda_s |\hat P(s)-\hat Q(s)|^2.
$$
In the IQP-QCBM notation,
$$
\mathcal{L}(\theta)=\operatorname{MMD}_k^2(p,q_\theta)=\mathbb{E}_{a\sim\Lambda}\!\left[(C_p^a-C_\theta^a)^2\right]=\sum_a \Lambda(a)(C_p^a-C_\theta^a)^2,
$$
where $C_r^a=\mathbb{E}_{x\sim r}[(-1)^{a\cdot x}]$ is the characteristic function and $C_\theta^a=\hat q_\theta(a)$ [2602.11042].

For general SBMs on $\mathbb{Z}_d^n$, the kernel is chosen through graph spectral analysis. Each site $i$ is assigned a weighted graph $G_i=(V_i,E_i,w_i)$ encoding the appropriate notion of label closeness, with Laplacian $L_i=D_i-A_i$. The product graph $G=G_1 \square \cdots \square G_n$ has Laplacian
$$
L=\sum_i I\otimes\cdots\otimes L_i\otimes\cdots\otimes I,
$$
eigenvectors $u_k=\bigotimes_i u^{(i)}_{k_i}$, and eigenvalues $\Lambda_k=\sum_i \lambda^{(i)}_{k_i}$. With spectral kernel $K=f(L)$, typically the heat kernel $f(\lambda)=e^{-t\lambda}$, the MMD again diagonalizes:
$$
\operatorname{MMD}^2(p,q_\theta)=\sum_k f(\Lambda_k)\,|\tilde p(k)-\tilde q_\theta(k)|^2.
$$
For shift-invariant graphs such as the cycle $C_d$ for ordinal variables and the complete graph $K_d$ for categorical variables, the $u_k$ are $\mathbb{Z}_d^n$ characters and
$$
\tilde q_\theta(k)=d^{-n/2}\langle \mathcal{D}(k,0)\rangle_\theta.
$$
The MMD then reduces to a mixture of squared Heisenberg–Weyl moment differences,
$$
\operatorname{MMD}^2(p,q_\theta)=\mathbb{E}_{k\sim\mathcal{P}_t}\!\left[\,|\langle \mathcal{D}(k,0)\rangle_p-\langle \mathcal{D}(k,0)\rangle_\theta|^2\,\right],
$$
with
$$
\mathcal{P}_t(k)=\frac{e^{-t\Lambda_k}}{\sum_{k'} e^{-t\Lambda_{k'}}}.
$$
For $C_d$, $\mathcal{P}_t(k_i)\propto \exp(-4t\sin^2(\pi k_i/d))$, whereas for $K_d$, $\mathcal{P}_t(0)=\pi_0$ and $\mathcal{P}_t(k\neq 0)=(1-\pi_0)/(d-1)$ [2607.06675].

This spectral diagonalization is the mechanism that makes SBMs classically trainable at scale. In the IQP setting, the same principle appears as a Walsh-diagonal positive-definite kernel
$$
k(x,y)=\sum_{s\in\{0,1\}^n}\lambda_s \chi_s(x)\chi_s(y), \qquad \lambda_s\ge 0,
$$
for which
$$
\operatorname{MMD}^2(p_\theta,p_{\mathrm{data}})=\sum_{s\neq 0}\lambda_s \big(\langle \chi_s\rangle_{p_\theta}-\langle \chi_s\rangle_{p_{\mathrm{data}}}\big)^2
=z(\theta)^\top \operatorname{diag}(\lambda)\, z(\theta),
$$
with residual vector $z_s(\theta)\equiv \langle Z_s\rangle_{p_\theta}-\langle Z_s\rangle_{p_{\mathrm{data}}}$ [2605.27735].

## 3. IQP realizations and trainability

An $n$-qubit IQP-QCBM prepares $U_\theta|+\rangle^{\otimes n}$ and measures all qubits in the $X$ basis, with
$$
U_\theta=\prod_{j=1}^D \exp(i\theta_j Z^{s^j}),
$$
where the generators are Pauli-$Z$ strings indexed by $s^j\in\mathbb{F}_2^n$. For $a\in\mathbb{F}_2^n\setminus\{0\}$, the characteristic function is
$$
C_\theta^a=\operatorname{Tr}[U_\theta^\dagger X^a U_\theta \rho],
$$
and depends only on the generators that anticommute with $X^a$. Defining
$$
S^a=\{s^j\mid a\cdot s^j=1\},
$$
the paper derives a subset-sum expansion over
$$
\Xi^a=\{J\subseteq S^a:\bigoplus_{j\in J} s^j = 0\},
$$
which yields exact formulas for the variances of $C_\theta^a$ and its partial derivatives under symmetric i.i.d. initialization [2602.11042].

The central trainability result is the uniform-initialization theorem. With i.i.d. uniform $\theta_j\sim \mathcal{U}[0,2\pi)$, define the critical rank
$$
r^a := \dim \operatorname{span}_{\mathbb{F}_2}(S^a).
$$
Then
$$
\operatorname{Var}_\theta[C_\theta^a]=2^{-r^a},
$$
and, for $s^\ell\in S^a$,
$$
\operatorname{Var}_\theta\!\left[\frac{\partial C_\theta^a}{\partial \theta_\ell}\right]=2^{2-r^a},
$$
while the derivative is zero for $s^\ell\notin S^a$. These variances monotonically decrease as more generators are added. Hence, the gradient magnitude is explicitly controlled by the linear-algebraic property $r^a$ [2602.11042].

This result makes the kernel spectrum and the generator topology jointly decisive. If the kernel spectrum is sufficiently flat or the architecture forces $r^a\ge c n$ for all $a\neq 0$, then $\operatorname{Var}_\theta[\partial \mathcal{L}/\partial \theta_\ell]\in \exp(-\Omega(n))$. Conversely, if $\Lambda$ concentrates an inverse-polynomial fraction on frequencies with $r^a=O(\log n)$, then
$$
\operatorname{Var}_\theta\!\left[\frac{\partial C_\theta^a}{\partial \theta_\ell}\right]=2^{2-O(\log n)}=\operatorname{poly}(1/n),
$$
avoiding exponential suppression. Hamming-weight-dependent kernels are “low-weight-biased” if $\lambda_s=f(|s|)$ concentrates mass on small $|s|$, and in structured IQP topologies such kernels can avoid exponential gradient suppression at lower-weight frequencies [2602.11042].

The architecture dependence is explicit in the examples analyzed. For product-state generators, $r^a=|a|$. For a 2D lattice with single-qubit and nearest-neighbor $Z\otimes Z$, $r^a=|N[A]|$ with $|N[A]|\in[|a|,5|a|]$. For sparse Erdős–Rényi graphs $G\sim G(n,p)$ with $p=c\ln n/n$, one has $r^a=|N[A]|$ and low-weight modes remain subexponential or polylogarithmic in rank. For the complete graph with all single- and two-qubit generators, $r^a=n$ for any $a\neq 0$, implying $\operatorname{Var}_\theta[\partial C^a]=2^{2-n}$ and exponential suppression everywhere; in this case no kernel can prevent a barren plateau [2602.11042].

The same analysis identifies an alternative trainable regime. If $\theta_j$ are i.i.d. Gaussian with $\theta_j\sim \mathcal{N}(0,\gamma^2)$ and $\gamma^2=c/D$ for constant $c$ and $D=\operatorname{poly}(n)$, then for any frequency $a$ and any generator $s^\ell\in S^a$,
$$
\operatorname{Var}_\theta\!\left[\frac{\partial C_\theta^a}{\partial \theta_\ell}\right]=\Omega(\operatorname{poly}(n^{-1})).
$$
Under mild conditions, the lower bound for the MMD gradient is polynomial, avoiding barren plateaus. The paper notes that this concentrates $U_\theta$ near identity, which can ease classical simulability; use carefully when preserving potential advantage matters [2602.11042].

The complexity-theoretic picture does not coincide with trainability. Product and 2D lattice families are not $L^2$ anti-concentrated. Sparse ER with $p=c\ln n/n$, $c>2$, satisfies $\mathbb{E}_G[\sum_a 2^{-r^a(G)}]=2+o(1)$ and is anti-concentrated. The complete graph is also anti-concentrated but untrainable at initialization. The existence theorem in the sparse ER case therefore connects a trainable low-weight regime with hardness arguments based on anti-concentration and the Bremner-style polynomial-hierarchy consequences for additive-error classical sampling [2602.11042].

## 4. Classical surrogation, truncation, and deployment discrepancy

The spectral viewpoint is independent of the loss function. QCBMs are identified as a quantum Fourier model independently of the loss function, because the Born distribution over computational-basis bitstrings is exactly the Walsh–Hadamard transform of a set of $Z$-string correlators. This permits direct spectrum truncation. The paper studies $k$-order truncation,
$$
p^{(k)}(x;\theta)=2^{-n}\sum_{s:|s|\le k}\langle Z^s\rangle_{\rho(\theta)}\,(-1)^{s\cdot x},
$$
and Random Fourier Correlator truncation,
$$
p^\Omega(x;\theta)=2^{-n}\sum_{s\in\Omega}\langle Z^s\rangle_{\rho(\theta)}\,(-1)^{s\cdot x}.
$$
These approximations preserve normalization if $s=0$ is included, but not positivity; thus, training must use losses defined on $\mathbb{R}$ [2511.01845].

The omitted spectrum quantitatively controls the approximation error. If $D_k(x):=p(x)-p^{(k)}(x)$, then the deterministic upper bound is
$$
|D_k(x)|\le 1-\frac{N_k}{2^n}, \qquad N_k=\sum_{p=0}^k \binom{n}{p},
$$
and for Haar-random unitaries,
$$
\mathbb{E}[D_k(x)^2]=\frac{2^n-N_k}{2^{3n}}.
$$
The paper recasts training as regression in an RKHS induced by the parity kernel
$$
s(x,x')=2^{-n}\sum_{s\in\Omega} (-1)^{s\cdot (x\oplus x')},
$$
which yields random-feature approximations over frequencies $s$ with a sampling distribution proportional to $|\langle Z^s\rangle|$ [2511.01845].

The resulting dequantization theorem states that if the sampling distribution over frequencies is efficiently samplable, polynomially concentrated, and aligned with the optimal quantum spectrum, then RFC-based classical training achieves, with high probability and $M,D\in \operatorname{poly}(n,\epsilon^{-1})$, an excess risk within $\epsilon$ of the optimal quantum risk. The paper also gives a discrepancy theorem for train-classical, deploy-quantum pipelines:
$$
|R[p^*_{D,c_{C\ell}}]-R[p_Q(\theta^*_{C\ell})]|
\le C\Big(\|c_{C\ell}(\theta^*_{C\ell})-c_Q(\theta^*_Q)\|+\|c_Q(\theta^*_Q)-c_Q(\theta^*_{C\ell})\|\Big).
$$
The first term is the feature-cardinality gap; the second is surrogate mismatch. Together, they quantify deployment-induced distribution shifts [2511.01845].

Two surrogate families are analyzed in detail. Tensor-network surrogates use Matrix Product States and RMPS averages, with closed-form variance formulas for marginals and truncated probabilities. Pauli-propagation surrogates propagate observables in the Heisenberg picture and, for IQP circuits, give a closed-form surrogate for truncated probabilities with flip-budget $h$. The paper studies IQP circuits, matchcircuits, Heisenberg-chain circuits, and Haldane-chain circuits, derives analytic variance formulas for matchcircuits, and characterizes a new dynamical Lie algebra for the Haldane chain. The overarching conclusion is that not all correlators are trainable; high-order terms may vanish or offer little gradient signal individually, which motivates spectral truncation and surrogate training on a controllable subset [2511.01845].

## 5. Implementations, software, and empirical demonstrations

The software and hardware-facing realization of SBMs is now split between fully general qudit implementations and specialized IQP pipelines. The paper “Spectral Born machines: classically trainable quantum generative models for discrete data” makes the SBM training machinery available in a new *tcdq* module of the PennyLane software platform. The module provides batched HW moment estimators implementing the $J/E$ factoring, graph spectral kernel utilities for cycle and complete graphs, unbiased MMD U-statistic estimators with JAX autodiff, gate-set construction helpers for weight/degree restrictions and sparse $\mathcal{G}$ selection by empirical Fourier magnitude, and GPU-ready pipelines for large $|\theta|$, $|\mathcal{Z}|$, and $|\mathcal{K}|$ [2607.06675].

For tractable $\Phi_\theta$, expectation values of HW displacements can be estimated to inverse polynomial additive precision by uniform sampling over $z\in \mathbb{Z}_d^n$:
$$
\langle \mathcal{D}(k,m)\rangle_\theta
=
\mathbb{E}_{z\sim\mathrm{Unif}}
\left[
\exp\!\left(i\pi m\cdot(2z-k)/d + i\Delta_\theta^k(z)\right)
\right],
$$
where $\Delta_\theta^k(z)=\Phi_\theta(z)-\Phi_\theta(z\ominus k)$. This is unbiased, with $O(1/\epsilon^2)$ samples for additive error $\epsilon$. Batching many moments and $z$-samples can be reduced to matrix multiplications, with overall complexity
$$
O(|\mathcal{O}||\mathcal{Z}|n+\sum_w 2^w |\mathcal{G}_w| (|\mathcal{O}|+|\mathcal{Z}|+|\mathcal{O}||\mathcal{Z}|)),
$$
which is practical with low weights such as one- and two-site gates [2607.06675].

A distinct implementation line appears in the 64-qubit calorimeter study, which trains an IQP Born machine on real high-energy-physics calorimeter shower images and compiles the trained model into a single sampling-hard IQP circuit for quantum deployment. The pipeline has three components: a Mixture-of-IQP (\moiqp{}) architecture, a Pearson-Stabilized Correlation Kernel (\psck{}), and an exact deferred-measurement compilation of \moiqp{} into a single IQP circuit on $\nfeat + \lceil \log_2 \Lcomp \rceil$ qubits (\ciqp{}). With $D=n_{\mathrm{feat}}=8$, $B=8$, and $L_{\mathrm{comp}}=8$, the compiled circuit acts on $67$ qubits. Across five seeds at $\Lcomp = 8$, $1500$ epochs, the model reaches $\maerho = 0.069 \pm 0.008$ against a $0.052$ encoding-fidelity floor on the training split and $0.071 \pm 0.008$ on a held-out test split, versus a Liu--Wang baseline at $\maerho = 0.100$. The compiled \ciqp{} reproduces the \moiqp{} marginal to $0.591 \pm 0.012$ times the Monte Carlo noise floor [2605.27735].

The reported large-scale SBM demonstrations extend beyond qubit IQP models:

| Study | System and model | Reported outcome |
|---|---|---|
| Potts/clock model | $36$ qudits, $d=16$ | All models reach similar train/test $\operatorname{MMD}^2 \approx 0.0016$ |
| Ribosomal RNA | $n=95$ qudits, $d=4$; $190$-qubit circuit when encoded to qubits | Model 1 final $\operatorname{MMD}^2 \approx 0.0020$; Model 2 $\approx 0.0035$ |
| Calorimeter images | IQP Born machine at $64$ qubits, compiled to $67$ qubits | $\operatorname{MAE}_\rho = 0.069 \pm 0.008$ train, $0.071 \pm 0.008$ test |

In the synthetic cyclic-data experiment, four SBMs on a $6\times 6$ periodic lattice with $d=16$ and train/test size $10^4$ each were trained for $10000$ iterations; the Degree 1, Degree 2, Degree 3, and All degrees models used $2592$, $10224$, $22896$, and $142290$ parameters, respectively, and all reached similar train/test $\operatorname{MMD}^2 \approx 0.0016$. The largest converges fastest, consistent with overparameterization benefits without overfitting. In the ribosomal RNA experiment, positions were filtered to retain $93$ positions out of $120$, with train size $428$ and test size $53$; the two-qudit model used $40{,}470$ parameters, the two-qudit plus selected three-qudit model used $1{,}297{,}661$ parameters, and the reported observation was “no obvious overfitting despite $>1$M parameters” [2607.06675].

## 6. Limitations, misconceptions, and open questions

Several limitations are structural rather than implementation-specific. In the IQP trainability analysis, the average-case lower bound uses Assumption 2, namely an unstructured target distribution ensemble with mean-zero and pairwise-uncorrelated characteristic values; highly structured targets can violate it, and tighter analysis of cross-frequency covariances is then needed. Assumption 1, symmetric i.i.d. initialization, is crucial for the exact subset-sum variance formulas; asymmetric or correlated initializations require modified analysis. The small-variance Gaussian bound assumes $D=\operatorname{poly}(n)$ so that $1/D=\operatorname{poly}(n^{-1})$ [2602.11042].

A common misconception is that classical hardness of sampling and trainability necessarily coincide. The architecture comparison in the IQP setting shows the opposite. Product and 2D lattice families are trainable at low weights but are not anti-concentrated. The complete graph is anti-concentrated but not trainable at initialization. Sparse Erdős–Rényi architectures with $p=c\ln n/n$ occupy the regime in which anti-concentration coexists with partial-spectrum trainability at low Hamming weight. This suggests that the feasible operating regime is architecture- and kernel-dependent rather than generic [2602.11042].

A second misconception is that spectral training guarantees accurate distribution learning in stronger distances. Minimizing MMD does not guarantee small total variation distance $\operatorname{TV}(p,q_\theta)$; MMD is controlled by TV but not vice versa. In the truncated-spectrum setting, pseudo-probabilities can be negative, so losses must be defined on $\mathbb{R}$. The dequantization conditions in RFC-based learning hinge on access to the optimal spectrum and its sparsity or low-degree structure; if the learned quantum distribution’s mass is spread across high-order frequencies, RFC dequantization fails [2511.01845].

At the model-design level, kernel choice sensitivity is substantial. A mismatch, such as using complete graphs for ordinal data, can slow training or reduce accuracy. The group structure itself can also be mismatched to the data; if the data do not exhibit cyclic or categorical structure, $\mathbb{Z}_d$ assumptions may not help. On hardware, approximate QFT compilation errors or noise degrade sampling fidelity, even though training is classical. The papers identify several open lines: broader dequantization criteria beyond shift-invariant kernels and RFC alignment, robustness-to-noise analyses in the spectral picture, non-Abelian SBMs aligned with domain symmetries, and extensions to broader architectures including IQP universality with ancilla, boson sampling or continuous-variable Spectral Born Machines, and photonic implementations [2607.06675, 2511.01845].

Within these constraints, the published SBM literature now delineates a coherent picture. SBMs are QCBMs viewed through their Fourier spectrum, with Fourier coefficients equal to directly trainable correlators or Heisenberg–Weyl moments; they admit graph-spectral MMD objectives, classically tractable estimators, and train-classical/deploy-quantum workflows; and their practical performance depends on the interaction between kernel spectrum, generator topology, initialization, and the spectral structure of the data [2607.06675].

Source: https://www.emergentmind.com/topics/spectral-born-machines