---
title: Spatial Demultiplexing (SPADE) in Optical Imaging
url: https://www.emergentmind.com/topics/spatial-demultiplexing-spade
type: topic
---

# Spatial Demultiplexing (SPADE) in Optical Imaging

Spatial-mode demultiplexing, commonly abbreviated **SPADE**, is a far-field optical measurement scheme in which the image-plane field is not measured in the position basis but is instead projected onto an orthonormal basis of spatial modes matched to the point-spread function (PSF), followed by photon counting or intensity readout. In the subdiffraction regime, this re-encodes separation, size, and higher-order structural information into mode occupancies rather than into small distortions of a blurred intensity profile. Foundational analyses showed that, for incoherent imaging, SPADE can estimate second and higher moments of an object much more precisely than direct imaging under photon shot noise, and can approach the quantum-optimal precision for location and scale parameters [1608.03211, 1703.08833].

## 1. Measurement principle

In the standard incoherent-imaging model, the one-photon state at the image plane is written as
\[
\rho_1 = \int d^2\mathbf R\, F(\mathbf R)\, |\psi_{\mathbf R}\rangle\langle \psi_{\mathbf R}|,
\]
with \(F(\mathbf R)\) the normalized source intensity distribution and \(|\psi_{\mathbf R}\rangle\) the shifted PSF state. Direct imaging measures the position-basis intensity
\[
f(\mathbf r)=\langle \mathbf r|\rho_1|\mathbf r\rangle,
\]
whereas SPADE measures mode probabilities
\[
p(j)=\langle \phi_j|\rho_1|\phi_j\rangle = \int d^2\mathbf R\,F(\mathbf R)\,|\langle \phi_j|\psi_{\mathbf R}\rangle|^2
\]
for a chosen orthonormal mode set \(\{|\phi_j\rangle\}\) [1608.03211].

For a Gaussian PSF, the natural basis is the Hermite–Gaussian (HG), or equivalently TEM, basis. In that setting, the fundamental mode is matched to the on-axis PSF, while small displacements or subdiffraction structure populate first-order and higher-order modes. This is the operational reason SPADE avoids the usual image-plane collapse of information: the parameter of interest is mapped into power in modes that are orthogonal to the dominant background [1703.08833].

A particularly transparent formulation appears for two incoherent point sources with known centroid. Expanding the field in the HG basis, the first-order modes \(HG_{10}\) and \(HG_{01}\) encode the \(x\)- and \(y\)-components of separation. In practical terms, SPADE is therefore a **mode-resolved intensity measurement** rather than a pixel-resolved image measurement [2008.02157].

This mode-based viewpoint generalizes beyond separation estimation. For a one-photon state from two incoherent sources, one formulation writes
\[
P(q\mid \theta)=\operatorname{Tr}\bigl[\ket{1,u_q}\bra{1,u_q}\rho_1\bigr]
\]
with \(u_q\) an orthonormal mode basis and \(\theta\) the source separation. The observable is the mode index and its count, not the image-plane coordinate [2509.17115].

## 2. Information-theoretic structure and Rayleigh’s curse

The central statistical contrast between direct imaging and SPADE is expressed through Fisher information. In direct imaging, when two incoherent point sources become very close, the image intensity changes only weakly with separation, and the Fisher information tends to zero. This is the standard **Rayleigh curse**. For example, in a direct-imaging model of two equally bright sources separated by \(2d\), the Fisher information behaves as \(F(d)\propto s\,d^2\) for \(d\ll \sigma\), so the Cramér–Rao bound diverges as \(d\to 0\) [1911.04744].

By contrast, for Gaussian PSFs and ideal HG-mode demultiplexing, the SPADE Fisher information remains finite in the sub-Rayleigh regime. One formulation gives
\[
\mathcal{F}^{\mathrm{SPADE}}(\theta)=\frac{N}{4\sigma^2},
\]
independent of separation and equal to the quantum Fisher information in that model [2509.17115]. In a related derivation for two incoherent beams with known centroid,
\[
I_F^q \simeq \frac{N}{w_0^2}, \qquad I_F^c = \frac{N}{w_0^2},
\]
showing that ideal SPADE reaches the quantum-limited scaling in the small-displacement regime [2008.02157].

This advantage extends from separation to moment estimation. In the subdiffraction limit, direct imaging estimates high-order spatial moments poorly, with relative-error scaling
\[
\frac{\delta \hat\theta_\mu}{\theta_\mu}=\mathcal O(\Delta^{-\mu}),
\]
whereas SPADE and iSPADE achieve
\[
\frac{\delta \hat\theta_\mu}{\theta_\mu}=\mathcal O(\Delta^{-\lceil \mu/2\rceil}),
\]
which is the quantum-optimal scaling stated for that regime [2511.20790]. Foundational analyses therefore treat diffraction not as a fundamental limit on information extraction, but as a consequence of measuring in the position basis rather than in a PSF-adapted mode basis [1608.03211].

## 3. Bases, receiver architectures, and practical implementations

The canonical SPADE receiver uses HG modes, but several experimentally relevant variants appear in the literature. Full SPADE measures many orthogonal modes; **binary SPADE** coarse-grains this to a null mode and its orthogonal complement, as in
\[
\mathcal M_1=\{R_0=\ketbra{0}{0},\; R_1=\mathbb I-R_0\},
\]
which is often sufficient for subdiffraction detection and hypothesis testing [2407.01995].

Interferometric extensions are also important. Tsang’s original framework introduced interferometric TEM measurements to access odd and mixed moments through two-mode superpositions [1608.03211]. More recent work describes **iSPADE**, where neighboring HG modes are interfered so that odd moments can be recovered in addition to the even moments available from plain HG sorting [2511.20790].

For noise-robust motion estimation, a two-mode implementation called **PM-SPADE** uses the “plus-minus” modes
\[
\ket{\phi_\pm} \equiv \frac{1}{\sqrt{2}}(\ket{\phi_0}\pm \ket{\phi_1}),
\]
designed so that each output has occupancy near \(1/2\) for small displacements, improving robustness to background and dark counts [2504.04350].

Experimentally, SPADE has been realized with several mode-sorting technologies. A **Multi-Plane Light Conversion (MPLC)** device from Cailabs was used as a 9-mode demultiplexer converting free-space HG modes into spatially separated outputs [2008.02157]. A later MPLC implementation based on a folded SLM architecture used wavefront matching to sort HG modes and reported simulated average overlap around **98%**, experimental average overlap around **88%**, experimental average sorting fidelity around **90.6%**, and simulated fidelity around **98.3%** [2509.17115]. Commercial HG demultiplexers such as **Cailabs PROTEUS-C** were used in spectroscopy and bright-source experiments [2409.01190].

A simplified two-channel architecture is the **double-clad fiber (DCF)** receiver, which separates the fundamental mode into a single-mode core and higher-order content into a multimode cladding. This implements a binary-SPADE-type measurement and was used in the first on-sky binary-source hypothesis-testing demonstration [2606.18025].

## 4. Imperfections, crosstalk, and nonstationary scenes

The ideal SPADE advantage is not unconditional. Detector noise, crosstalk, misalignment, and source dynamics each modify the asymptotic picture.

For noisy detectors, the informative antisymmetric or first-order mode acquires a background floor. In a Poissonian model with background mean \(b\) and signal mean \(s\), the smallest resolvable separation scales as
\[
d_{\min}\sim \sigma\,\mathrm{SNR}^{-1/2},
\]
more precisely, under the paper’s half-resolution convention,
\[
d_{1/2}\equiv \frac{2\sigma}{\sqrt{\mathrm{SNR}}}.
\]
Thus, ideal arbitrarily-small-separation superresolution is lost once detector noise is included [1911.04744].

Mode-sorting **crosstalk** is another dominant limitation. For asymmetric source discrimination with balanced crosstalk
\[
C_{10}=C_{01}=C,\qquad C_{11}=C_{00}=1-C,
\]
one analysis derived the threshold
\[
C_{\mathrm{th}}=\frac{3-\sqrt 6}{6}\simeq 0.1,
\]
below which SPADE remains asymptotically superior to direct imaging in the relevant small-separation, low-intensity-ratio regime [2605.15929]. In separation estimation with arbitrary brightness imbalance, the picture is more nuanced: ideal SPADE has constant Fisher information \(w^2F_{\rm HG}(d)=1\), but for crosstalk-affected SPADE and any source imbalance, the paper states that SPADE performs worse than ideal direct imaging in the asymptotic limit \(x\to 0\); nevertheless, for practical sub-Rayleigh separations with \(x \gtrsim 3\sqrt{p_c}\), SPADE is effectively optimal [2211.09157].

The literature also contains explicit noise-mitigation proposals. For random-unitary mode-mixing noise generated by polynomials of creation and annihilation operators, repeated demultiplexers interlaced with phase-space rotations can decouple the noise in the limit of large repetition number and small noise strength. For displacement noise, a two-step sequence with parity,
\[
\mathcal G_1=\{1,\Pi\}, \qquad \Pi=e^{i\pi a^\dagger a},
\]
cancels identical displacement errors exactly when the noise is frozen between the two steps [2407.01995].

Dynamic scenes introduce a different limitation. Under Brownian motion of the source centroid, adaptive SPADE retains a short-timescale regime with \(d_{\min}/w \sim N^{-1/2}\), but in the long-timescale regime the Fisher information scales as \(x^2\), restoring Rayleigh-like behavior and giving \(d_{\min}/w \sim N^{-1/4}\) [2407.13723]. By contrast, for rotational and oscillatory dynamics of the source pair, HG-SPADE can remain free from Rayleigh’s curse: for isotropic random orientation one reported small-separation behavior is
\[
w^2 F(d)\approx \frac23 - \frac{8x^2}{9},
\]
which stays finite as \(d\to 0\) [2407.10507].

## 5. Separation estimation, discrimination, spectroscopy, and motion sensing

A substantial experimental literature treats SPADE as a practical superresolution instrument rather than only an asymptotic construct. Simultaneous multimode demultiplexing was used to estimate two-dimensional transverse separation with sensitivity around \(2\times 10^{-3}\) in \(d/w_0\), useful performance up to \(d/w_0\approx 1.2\), and raw scans from \(0\) to \(3w_0\) [2008.02157].

In the bright-source regime, SPADE has also been used for estimating both separation and relative intensity. For two bright incoherent point sources, one experiment reported a resolving power of about
\[
r \simeq 0.023
\]
in units of beam waist, and for relative-intensity estimation found a detection limit
\[
\epsilon_\ell \simeq 2.7\times 10^{-5},
\]
compared with
\[
\epsilon_\ell^{\mathrm{DI}} \simeq 5.7\times 10^{-4}
\]
for direct imaging, i.e. about **21 times worse** [2206.05246].

A major recent application is **asymmetric source discrimination** motivated by exoplanet detection. In one universal, parameter-independent test, photons in \(\mathrm{HG}_{01}+\mathrm{HG}_{10}\) are counted and compared with a threshold set from the \(H_0\) distribution at \(\alpha=0.05\). With crosstalk \(C\simeq 10^{-2}\), the paper reports that SPADE requires roughly
\[
\sim 1/(8C_{10}) \approx 12.5
\]
times fewer photons than direct imaging in the asymptotic regime [2605.15929]. A related imperfection analysis found that noisy SPADE and direct imaging can share the same realistic \(\nu^2 s^4\) scaling, but SPADE retains a superior prefactor and remains the most efficient method under the stated sub-Rayleigh noise conditions [2505.00064]. The first on-sky binary-SPADE instrument, based on a DCF receiver, demonstrated binary-star detection of **Alpha Centauri** below the diffraction limit in a photon-starved regime and reported type-II error always lower than perfect direct imaging, although heavily limited by unbalanced loss in the coupler [2606.18025].

SPADE has also been extended from imaging to **spectroscopy**. In a proof-of-principle star–planet experiment, light was sorted into HG modes before spectral analysis, so that the star was concentrated in \(\mathrm{HG}_{00}\) while the off-axis planet coupled preferentially to first-order modes. The paper reports that spectral discrimination becomes clear for \(d_a \gtrsim 0.2\) in its setup, whereas simulated direct detection needs roughly \(d_a \gtrsim 2\) for comparable separation [2409.01190].

Another application is **micro-oscillation frequency estimation**. In that setting, PM-SPADE outperformed direct imaging in the presence of background noise because it used only two outputs and maintained substantial occupancy near small displacement. The paper reports that, without background noise, both PM-SPADE and direct imaging approach the QCRB, but with background noise PM-SPADE remains much more stable and significantly outperforms direct imaging [2504.04350].

## 6. Extensions, finite-sample subtleties, and acronym disambiguation

Several extensions push SPADE beyond the standard single-photon, static, Gaussian scenario. In a biphoton setting based on SPDC, coincidence measurements after projection of both photons onto HG modes yielded a Fisher-information enhancement proportional to \(\sqrt{K}\), where \(K\) is the Schmidt number. In the ideal full-basis theory,
\[
\mathcal{FI}_{\mathrm{coinc}}^{\mathrm{tot}}=\frac{1}{2}\sqrt{K},
\]
reducing to the ordinary SPADE benchmark \(1/2\) when \(K=1\) [2212.10468].

Temporal fluctuations can also be exploited as a resource. **SOFSPADE** combines SPADE with emitter blinking and uses temporal cumulants of detector outputs. The paper reports that temporal fluctuations both improve precision and simplify hardware, because in the presence of blinking full even-moment information can be recovered with **image inversion interferometry (III)** instead of a full HG sorter. For the 8th moment, SOFSPADE reaches the same precision with roughly **\(10^3\) fewer frames** than mean-signal SPADE in the reported simulation [2511.20790].

The finite-sample literature adds an important caution. Under ideal alignment, SPADE attains the quantum-optimal large-sample Stein exponent for one-versus-two source discrimination, but a singular-learning analysis showed that ideal aligned benchmarks do not automatically survive misalignment. In a binary-SPADE reduction with offset \(\theta\), the paper identifies an exact blind separation
\[
s^\ast=2\theta,
\]
at which the coarse-grained null and alternative become identical and the test power collapses to \(\alpha\). On the plotted Neyman–Pearson grids, direct imaging was stronger than misaligned binary-SPADE despite the latter’s distinctive local scaling structure [2605.14432]. This suggests that ideal aligned asymptotics should not be conflated with finite-\(n\), imperfectly aligned performance.

The acronym **SPADE** is also overloaded outside this imaging literature. The temporal-to-spatial multiphoton-routing device in “Single-active-element demultiplexed multi-photon source” is explicitly *not* the spatial-mode demultiplexing scheme used in optical imaging, even though both involve demultiplexing into spatial outputs [2304.12956]. Likewise, “SPADE: Spatial Transcriptomics and Pathology Alignment Using a Mixture of Data Experts for an Expressive Latent Space” uses the same acronym for a pathology foundation model rather than for optical mode demultiplexing [2506.21857].

Across these variants, the stable core of SPADE remains the same: project the field onto a PSF-adapted orthonormal mode basis, measure the resulting mode occupancies, and use the fact that subdiffraction information is often concentrated in a small set of orthogonal modes rather than in the intensity profile of a blurred image.

Source: https://www.emergentmind.com/topics/spatial-demultiplexing-spade