---
title: SPADE Optical Demultiplexing
url: https://www.emergentmind.com/topics/spatial-mode-demultiplexing-spade
type: topic
---

# SPADE Optical Demultiplexing

Searching arXiv for recent and foundational papers on SPADE to ground the encyclopedia entry.
arXiv search query: spatial-mode demultiplexing SPADE incoherent imaging superresolution
Spatial-mode demultiplexing (SPADE) is a class of optical measurement schemes that replaces position-resolved image-plane detection by projection of the optical field onto an orthonormal spatial-mode basis followed by mode-resolved photon counting. In far-field incoherent imaging, SPADE is typically constructed from modes matched to the point-spread function (PSF), most often Hermite–Gaussian (HG) modes for a Gaussian PSF. Its central role is to concentrate subdiffraction information—especially separation, scale, and higher spatial moments—into a small number of informative channels while suppressing the bright, parameter-insensitive background that dominates direct imaging. In the canonical problem of discriminating or estimating properties of closely spaced incoherent sources, SPADE avoids the vanishing Fisher information of direct imaging at small separations, attains quantum-optimal scaling under ideal conditions, and has become a unifying framework across superresolution imaging, asymmetric hypothesis testing, spectroscopy, microscopy, and related mode-sorting architectures [1608.03211], [1703.08833], [2605.14432].

## 1. Measurement principle and optical construction

SPADE operates by projecting the image-plane field onto an orthonormal set of spatial modes tailored to the imaging system’s PSF, and then detecting the discrete mode label or mode-resolved intensity. For a Gaussian PSF, the natural choice is the HG basis. In one dimension, the Gaussian amplitude PSF may be written as
\[
\psi(x)=\left(\frac{1}{2\pi\sigma^2}\right)^{1/4}\exp\!\left(-\frac{x^2}{4\sigma^2}\right),
\]
with \(\int dx\,|\psi(x)|^2=1\). In the corresponding HG-SPADE receiver, each detected photon is assigned to a mode \(q\in\{0,1,2,\dots\}\), with projectors \(\Pi_q=|\phi_q\rangle\langle\phi_q|\) and detection probabilities
\[
P_a(q)=\operatorname{Tr}(\Pi_q \rho_a),
\]
where \(\rho_a\) denotes the one-photon state under hypothesis \(a\) [2605.14432].

A general construction beyond the Gaussian case uses PSF-adapted derivative-like modes. In the semiclassical formulation, one begins from the Fourier transform \(\Psi(k)\) of the PSF and constructs orthogonal polynomials \(g_q(k)\) with respect to the weight \(|\Psi(k)|^2\). The resulting point-spread-function-adapted (PAD) modes are
\[
\Phi_q(k)=(-i)^{|q|}g_q(k)\Psi(k), \qquad
\phi_q(x)=(-i)^{|q|}g_q(-i\partial)\psi(x),
\]
which remain orthonormal in real space [1703.08833]. For Gaussian \(\Psi(k)\), the \(g_q\) reduce to scaled Hermite polynomials, and the real-space PAD modes become the familiar TEM/HG modes [1703.08833].

The physical significance of this basis choice is that spatial perturbations of a centered source couple preferentially into higher-order modes. Under perfect alignment and a single on-axis source, all photons ideally populate the lowest mode \(q=0\). Small displacements cause “leakage” out of \(q=0\), with the first nontrivial contribution appearing in \(q=1\). This leakage is the core signal exploited by SPADE in sub-Rayleigh problems [2605.14432]. A widely used coarse graining is binary-SPADE, which records only whether a photon is in the fundamental mode \(q=0\) or in the complement \(q\ge 1\). Near alignment, this retains the full leading \(O(s^2)\) leakage contrast, while redistribution into higher modes enters only at \(O(s^4)\) [2605.14432].

Closely related architectures include interferometric PAD (iPAD), where superpositions of two spatial modes are interfered to access off-diagonal coherences in the modal mutual coherence matrix. If \(\varphi_{qq'}^\pm=(\phi_q\pm\phi_{q'})/\sqrt{2}\), then the measured powers encode \(\operatorname{Re}\Gamma_{qq'}\), extending SPADE from diagonal moment access to full moment tomography [1703.08833].

## 2. Statistical model and relation to direct imaging

In direct imaging, the measurement outcome is the photon position \(x\), with image-plane distribution
\[
p_0(x)=|\psi(x)|^2,\qquad
p_1(x)=(1-\epsilon)|\psi(x)|^2+\epsilon|\psi(x-s)|^2
\]
for the asymmetric one-versus-two-source problem, where \(s\) is the separation, \(\sigma\) the PSF width, and \(\epsilon\ll 1\) the brightness ratio of the secondary source [2605.14432]. The data are i.i.d. samples \(X_i\sim p_a\), and the likelihood over \(n\) detections is \(\prod_i p_a(X_i)\) [2605.14432].

SPADE changes only the measurement, not the underlying optical state. In the same problem, the quantum states per detected photon are
\[
|\psi_x\rangle:=\int dx'\,\psi(x'-x)\,a^\dagger(x')|0\rangle,
\]
\[
\rho_0=|\psi_0\rangle\langle\psi_0|,\qquad
\rho_1=(1-\epsilon)|\psi_0\rangle\langle\psi_0|+\epsilon|\psi_s\rangle\langle\psi_s|.
\]
A SPADE receiver measures the induced classical distribution over mode labels \(q\), producing i.i.d. outcomes \(Y_i=q_i\), or equivalently multinomial counts after \(n\) detections [2605.14432].

In a broader semiclassical treatment of incoherent imaging, one writes the mutual coherence function as
\[
\Gamma(x,x'|\theta)=\int dX\,\psi(x-X)\psi^*(x'-X)\,F(X|\theta),
\]
and the average image intensity as
\[
f(x|\theta)=\Gamma(x,x|\theta)=\int dX\,|\psi(x-X)|^2\,F(X|\theta),
\]
with \(F(X|\theta)\) the object intensity function. After linear-optical projection onto modes \(\phi_j(x)\), the mean power in channel \(j\) is
\[
p_j(\theta)=\int dX\left|\int dx\,\phi_j^*(x)\psi(x-X)\right|^2F(X|\theta),
\]
and photon counts are independent Poisson variables with means \(\tau p_j(\theta)\), where \(\tau=\eta T/(\hbar\omega)\) [1703.08833]. This Poisson channel model underlies both moment estimation and Fisher-information calculations.

The statistical contrast between direct imaging and SPADE is therefore not a change of source model but a change of sufficient statistics. Direct imaging samples a PSF-dominated intensity profile, whereas SPADE samples a discrete mode distribution in which the background-dominated fundamental mode can be rejected or separated from the informative leakage channels [1703.08833].

## 3. Superresolution, Fisher information, and quantum limits

SPADE is most closely associated with the removal of Rayleigh’s curse. In direct imaging of two equally bright incoherent point sources with Gaussian PSF, the Fisher information for separation vanishes as the separation tends to zero. In one formulation, for a Gaussian PSF with width \(\sigma\),
\[
F_{\mathrm{DI}}(d)\propto \frac{N d^2}{\sigma^4}\qquad (d\to 0),
\]
while HG-SPADE yields a separation-independent value
\[
F_{\mathrm{SPADE}}(d)=\frac{N}{4\sigma^2},
\]
matching the quantum Fisher information over the full separation range in the ideal Gaussian case [2509.17115]. An earlier moment-based treatment expresses the same phenomenon as
\[
\mathcal{J}^{(\mathrm{direct})}(d)\approx \frac{N d^2}{8},\qquad
\mathcal{J}^{(\mathrm{TEM})}(d)\approx \frac{N}{4},
\]
for small \(d\), again showing direct-imaging collapse and constant SPADE information [1608.03211].

The mechanism is transparent in the small-displacement expansion of the modal probabilities. For a displaced Gaussian source, the HG-mode probabilities are Poisson in the mode index:
\[
P_s(q)=e^{-\tau}\frac{\tau^q}{q!},\qquad \tau=\frac{s^2}{4\sigma^2}.
\]
Thus,
\[
P_s(0)=1-\tau+O(\tau^2),\qquad
P_s(1)=\tau+O(\tau^2),\qquad
\sum_{q\ge 2}P_s(q)=O(\tau^2)=O(s^4),
\]
so the first-order leakage channel alone already contains the leading \(O(s^2)\) separation signal [2605.14432]. This explains why binary-SPADE and first-order HG detection can perform nearly optimally in the sub-Rayleigh regime.

In the language of object moments, SPADE also linearizes access to higher-order structure. For a subdiffraction object with support width \(\Delta\ll 1\), the PAD-mode mutual coherence matrix satisfies
\[
\Gamma_{qq'}=\sum_{r,r'} H_{qr}H_{q'r'}\theta_{r+r'},
\]
so \(\Gamma_{qq'}\) is chiefly sensitive to the moment \(\theta_{q+q'}\) at leading order. The resulting estimator variances scale as
\[
\mathrm{Var}[\widehat{\theta}_\mu']\sim \frac{\theta_0^2}{N_s}O(\Delta^{|\mu|}),
\]
whereas direct-imaging Cramér–Rao bounds do not improve with shrinking object size in the same way [1703.08833]. For Gaussian PSF and two equal incoherent point sources, SPADE yields a finite separation Fisher information \(N_s/4\) as \(d\to 0\), while direct imaging gives \(F^{\mathrm{DI}}(d)\to 0\) [1703.08833].

The same logic extends to other parameters. For location and scale parameters of subdiffraction objects, SPADE approaches the optimal precision allowed by quantum mechanics in the small-object limit [1608.03211]. This suggests that SPADE is not merely a special-purpose separation estimator but a measurement paradigm aligned with the quantum geometry of incoherent imaging.

## 4. Singular asymptotics and finite-photon hypothesis testing

Beyond Fisher information and large-sample exponents, recent work has analyzed SPADE near singular model boundaries, where regular asymptotics break down. In the asymmetric one-versus-two-source discrimination problem, the null hypothesis is realized at the singular point \((\epsilon,s)=(0,0)\), so the Kullback–Leibler function vanishes on a nonregular boundary. Singular learning theory then describes the Bayes free energy through the real log canonical threshold (RLCT) \(\lambda\) and multiplicity \(m\) [2605.14432].

For the aligned Gaussian problem, direct imaging and SPADE have different local divergences,
\[
D_{\mathrm{DI}}(\epsilon,s)\asymp \epsilon^2 s^2,\qquad
D_{\mathrm{SPADE}}(\epsilon,s)\asymp \epsilon s^2,
\]
leading to zeta functions
\[
\zeta_{\mathrm{DI}}(z)\sim \frac{C_{\mathrm{DI}}}{(2z+1)^2},\qquad
\zeta_{\mathrm{SPADE}}(z)\sim \frac{C_{\mathrm{SPADE}}}{(z+1)(2z+1)}.
\]
Both schemes share the same rightmost pole \(z=-1/2\), hence the same RLCT \(\lambda=1/2\), but their multiplicities differ:
\[
(\lambda_{\mathrm{DI}},m_{\mathrm{DI}})=\left(\frac12,2\right),\qquad
(\lambda_{\mathrm{SPADE}},m_{\mathrm{SPADE}})=\left(\frac12,1\right).
\]
This yields Bayes free-energy asymptotics
\[
F_n^{\mathrm{DI}}=\frac{1}{2}\log n-\log\log n+O_p(1),\qquad
F_n^{\mathrm{SPADE}}=\frac{1}{2}\log n+O_p(1),
\]
so aligned SPADE has a universal subleading advantage in the local prior-weighted singular regime [2605.14432].

In the same setting, the ideal aligned SPADE classical Stein exponent matches the quantum benchmark to leading order. Direct imaging has
\[
D(p_0\|p_1)=\frac{\epsilon^2 s^2}{2\sigma^2}+O(\epsilon^2 s^4,\epsilon^3),
\]
whereas the quantum Stein exponent is
\[
D(\rho_0\|\rho_1)=\frac{\epsilon s^2}{4\sigma^2}+O(\epsilon s^4,\epsilon^2),
\]
and aligned SPADE achieves
\[
D_{\mathrm{SPADE}}(\epsilon,s)=\frac{\epsilon s^2}{4\sigma^2}+O(\epsilon s^4,\epsilon^2),
\]
matching the quantum limit to leading order [2605.14432].

However, finite-photon performance is not determined solely by asymptotic exponents. In misaligned binary-SPADE, the local Kullback–Leibler divergence scales as \(\epsilon^2 s^4\), and nontrivial local power first appears on the intrinsic scale \(s=O(n^{-1/4})\), whereas direct imaging appears on \(s=O(n^{-1/2})\) [2605.14432]. This suggests a formal asymptotic advantage for misaligned binary-SPADE, but numerical Neyman–Pearson comparisons under common physical conditions in that analysis show direct imaging stronger on the plotted grids and reveal an exact blind separation \(s^\ast=2\theta\), at which the binary-SPADE power collapses to the test size \(\alpha\) [2605.14432]. A plausible implication is that singular asymptotics organize the local theory, but experimental utility depends sensitively on which asymptotic regime is actually entered.

## 5. Imperfections, crosstalk, misalignment, and robustness

SPADE’s theoretical advantage is contingent on mode fidelity. Several distinct imperfections recur across the literature: misalignment between the sorter basis and the source centroid, crosstalk between modal channels, background and dark noise, finite bandwidth, and dynamic centroid motion.

### Misalignment

Misalignment replaces the ideal support mismatch of aligned SPADE by a genuine fluctuating null. In the binary-SPADE reduction with detector offset \(\theta\), the \(q=0\) success probabilities are
\[
p_0(\theta)=e^{-\gamma^2},\qquad
p_s(\theta)=e^{-(\gamma-g_s)^2},\qquad
\gamma=\frac{\theta}{2\sigma},\qquad g_s=\frac{s}{2\sigma}.
\]
Near alignment, the Bernoulli contrast obeys
\[
\Delta_{\epsilon,s}(\theta)=a(\theta)\,u(\epsilon)\,s^2+O(s^3),
\qquad
u(\epsilon)=\epsilon(1-\epsilon),
\]
with
\[
a(\theta)=\frac{p_0(\theta)[2\gamma^2-1]}{4\sigma^2}.
\]
This preserves the leading \(O(s^2)\) leakage contrast but removes the ideal null singularity [2605.14432]. In adaptive or Brownian settings, centroid fluctuations can reintroduce Rayleigh-type scaling if the measurement interval is too long. For a centroid undergoing Brownian motion with diffusion constant \(D\), and segment duration \(T\), the dimensionless parameter \(\tau=DT/w^2\) controls the regime. For short times \(x\gg \sqrt{\tau}\), where \(x=s/(2w)\), the per-photon FI remains approximately constant, but for long times \(x\ll \sqrt{\tau}\), the FI scales as \(x^2\), so Rayleigh’s curse reappears [2407.13723]. Adaptive realignment with \(\sqrt{\tau}=\kappa x\) and small \(\kappa\) preserves near-optimal performance [2407.13723].

### Crosstalk

Crosstalk has been analyzed both phenomenologically and experimentally. In an asymmetric one-versus-two-source test using the combined odd channel, balanced crosstalk with probability \(C\) modifies the Bernoulli distributions via
\[
\begin{pmatrix}P'(0|H)\\P'(1|H)\end{pmatrix}
=
\begin{pmatrix}C_{00} & C_{01}\\ C_{10} & C_{11}\end{pmatrix}
\begin{pmatrix}P(0|H)\\P(1|H)\end{pmatrix}.
\]
The resulting leading SPADE error exponent becomes
\[
D_{\mathrm{SD}}(p_0'\|p_1')\simeq
\frac{\epsilon^2 s^4}{32}\frac{(C_{11}-C_{10})^2}{C_{10}(1-C_{10})}
\simeq \frac{\epsilon^2 s^4}{32C_{10}}
\]
for small \(C_{10}\), while direct imaging has \(D_{\mathrm{DI}}\simeq \epsilon^2 s^4/4\). Balanced SPADE remains exponentially superior provided
\[
C\le C_{\mathrm{th}}=\frac{3-\sqrt{6}}{6}\simeq 0.1.
\]
At \(C_{10}\approx 10^{-2}\), SPADE requires about \(0.08N_{\mathrm{DI}}\), or roughly \(12.5\times\) fewer photons, to achieve the same error rate [2605.15929]. This threshold result is notable because it identifies a concrete imperfection level below which SPADE’s practical advantage survives.

A different crosstalk analysis for resolving separation between unbalanced sources finds that, for any imbalance, crosstalk-affected SPADE performs worse than ideal direct imaging in the strict vanishing-separation limit, although it remains superior over several orders of magnitude of sub-Rayleigh separations above an imbalance-dependent threshold \(x_c(\nu)\approx c(\nu)\sqrt{p_c}\) [2211.09157]. This is not a contradiction: one result addresses asymmetric hypothesis testing with calibrated odd-channel statistics, while the other studies asymptotic estimation FI for unbalanced two-source separation. Together they indicate that the effect of crosstalk depends materially on the inference task.

### Detector and background noise

With noisy detection, SPADE’s resolution becomes signal-to-noise-ratio limited. In a binary-SPADE or first-order mode detector with background \(b\), the small-separation FI is suppressed by the factor \(1/[1+4\sigma^2\beta/d^2]\), where \(\beta=1/\mathrm{SNR}\), so the minimum resolvable separation scales as
\[
d_{\min}\propto \mathrm{SNR}^{-1/2}
\]
rather than vanishing arbitrarily [1911.04744]. Photon counting retains the largest advantage; homodyne and heterodyne variants also beat direct imaging in an SNR-dependent window but are limited to a maximum FI of \(Q/4\) [1911.04744].

### Noise suppression protocols

Noise need not only be mitigated passively. For random unitary noise generated by polynomials of creation and annihilation operators, repeated demultiplexers interleaved with rotation operators can decouple the noise in the large-repetition, weak-noise limit. For displacement noise, a simpler two-pass protocol with an interleaved parity operator can achieve perfect decoupling if the noise configuration is frozen between passes, restoring SPADE’s constant small-separation FI [2407.01995]. This suggests a direct connection between SPADE robustness and dynamical-decoupling ideas from quantum control.

## 6. Implementations, applications, and extensions

SPADE has moved well beyond the canonical two-source thought experiment. Its implementations now span mode sorters, interferometric schemes, multi-plane light conversion (MPLC), digital holography, double-clad fibers, and spectroscopy platforms.

### Multi-plane light conversion and mode sorters

MPLC is a prominent implementation technology because it realizes arbitrary unitary spatial-mode transforms with cascaded phase profiles and free-space propagation. Reflective folded MPLC architectures have been used to sort low-order HG modes with experimental average fidelity \(90.6\%\) for a three-mode sorter and simulated diagonal fidelities \(0.94\)–\(0.96\) for a six-mode design [2509.17115]. These devices support direct SPADE measurements for superresolution experiments and make the ideal HG-basis picture experimentally concrete.

A larger 45-mode MPLC multiplexer/demultiplexer based on a separable HG basis achieves average insertion loss of about \(4\) dB and average crosstalk of about \(-28\) dB across the C band [1803.07907]. Although developed for space-division multiplexing rather than imaging, it implements precisely the kind of high-dimensional orthonormal projection required by SPADE, and therefore functions as a scalable mode-sorting backend for large-mode-count SPADE architectures.

### Distance and intensity metrology

In the bright-source regime, SPADE can estimate transverse separation and relative intensity of two incoherent point sources with large dynamic range. In one implementation using a commercial HG demultiplexer at \(1.55\,\mu\mathrm{m}\), the measured odd-mode signals \(HG_{01}\) and \(HG_{10}\) scaled quadratically with sub-Rayleigh displacement, yielding a resolving power \(r\simeq 0.023\) of the beam waist and a relative-intensity detection limit \(\epsilon_\ell\simeq 2.7\times 10^{-5}\), about \(21\times\) better than direct imaging in the same apparatus [2206.05246].

A simultaneous nine-mode HG demultiplexing experiment likewise demonstrated 2D distance estimation beyond the Rayleigh limit with a large dynamic range, while also showing that calibrated crosstalk acts as an additive offset in first-order channels and must be explicitly modeled [2008.02157].

### Asymmetric source discrimination and exoplanet-inspired tests

In asymmetric one-versus-two-source detection, SPADE is especially attractive because odd modes or leakage channels suppress the bright on-axis source. A 2026 tabletop experiment implemented a parameter-independent test based on odd-channel counts and showed that SPADE retains an advantage over direct imaging at measured crosstalk \(C\approx 10^{-2}\), with false-negative rates well below direct imaging in the small-separation, low-brightness-ratio regime [2605.15929].

This theme extends to spectroscopy. A proof-of-principle experiment sorted light from overlapping incoherent “star” and “planet” sources into HG modes and measured spectra in each mode. The star coupled almost entirely to \(HG_0\), while the off-axis planet preferentially populated \(HG_1\). In the overlapping-PSF regime, the Fisher information for planet spectroscopy scaled as \(\epsilon\) under HG-SPADE but as \(\epsilon^2\) under direct detection, a quadratic gain in information [2409.01190]. The implementation used an MPLC sorter and SNSPDs, with measured first-order cross-talk \(\chi\simeq 0.0035\) [2409.01190].

An on-sky binary-SPADE demonstration has now been reported using a double-clad fiber coupler that separates a PSF-matched fundamental mode from its complement. In the photon-starved regime, with fewer than \(5000\) photons per measurement, the instrument detected a binary star below the diffraction limit and achieved type-II error rates consistently lower than perfect direct imaging, though unbalanced loss in the coupler substantially limited the scaling relative to ideal SPADE [2606.18025]. This constitutes the first on-sky demonstration of binary-source hypothesis testing using SPADE-based detection.

### Dynamic sources and motion estimation

SPADE remains useful for dynamic scenes. For rotating or oscillating incoherent source pairs with Gaussian PSF, the small-separation SPADE FI remains strictly positive under broad orientation dynamics, while direct imaging often reverts to \(d^2\)-scaling [2407.10507]. Under random orientation on the sphere,
\[
w^2F(d)\approx \frac{2}{3}-\frac{8}{9}x^2,\qquad
w^2F_{\mathrm{DI}}(d)\approx \frac{16}{9}x^2,
\]
with \(x=d/(2w)\) [2407.10507]. The same work proposes a rotation-averaging algorithm that removes the need to estimate the source orientation angle while preserving near-quantum-limited precision [2407.10507].

SPADE has also been used for frequency estimation of a micro-oscillating point source. A two-mode “plus-minus” SPADE, using the superpositions \(|\phi_\pm\rangle=(|\phi_0\rangle\pm|\phi_1\rangle)/\sqrt{2}\), concentrates position information into only two channels. In the presence of background, its Fisher-information penalty scales as \(1/(1+2b/\nu)\), making it substantially more robust than direct imaging and standard HG-SPADE for estimating oscillation frequency [2504.04350].

### Fluctuation-enhanced SPADE

Temporal fluctuations can simplify as well as improve SPADE. With blinking emitters, temporal cumulants of low-order SPADE outputs encode higher spatial moments. In the subdiffraction regime, for SPADE mode \(j\),
\[
\kappa(I_{j_1}^S,\dots,I_{j_r}^S)
=
\frac{\tilde\kappa_r}{\prod_{l=1}^r 4^{j_l} j_l!}\,
\theta_{2(j_1+\cdots+j_r)}
+O(\Delta^{2(j_1+\cdots+j_r)+2}),
\]
so higher-order cumulants improve higher-moment estimation [2511.20790]. Remarkably, in the presence of fluctuations, full even-moment recovery can be obtained with the much simpler image inversion interferometer, because
\[
\kappa^{(r)}(I_-^{\mathrm{III}})=\frac{\tilde\kappa_r}{4^r}\theta_{2r}+O(\Delta^{2r+2}),
\]
allowing fluctuation-enhanced parity sorting to replace full SPADE for even moments [2511.20790]. This suggests that blinking can act as an informational resource rather than merely a nuisance.

### Quantum and communication extensions

SPADE has also entered explicitly quantum regimes. For spatially entangled bi-photons, coincidence detection after HG projection yields a Fisher-information enhancement by \(\sqrt{K}\), where \(K\) is the Schmidt number of the two-photon state. In the 2D case,
\[
F_{2\mathrm{ph}}^{(2D)}=\frac{\sqrt{K}}{2},
\]
showing that entanglement enhances superresolution beyond ordinary SPADE [2212.10468].

Outside imaging, the same projective mode-sorting principle appears in communications-oriented demultiplexers. A tri-degree-of-freedom multi-vortex geometric beam basis supports a SPADE-like matched-filter receiver with high divergence degeneracy and low bit-error rate, illustrating that the abstract notion of SPADE extends naturally to large-scale mode demultiplexing beyond superresolution imaging [2110.08815].

---

SPADE is therefore best understood not as a single instrument but as a measurement principle: project onto a basis adapted to the PSF and the inference task, then count in the resulting mode channels. In ideal aligned Gaussian settings, this principle yields quantum-optimal discrimination and estimation, eliminates Rayleigh’s curse, and gives direct access to low-order object moments [1608.03211], [1703.08833], [2605.14432]. In realistic systems, performance depends critically on alignment, crosstalk, noise, and temporal dynamics; under sufficiently good calibration, those imperfections can be tolerated, modeled, or partially corrected [2605.15929], [2407.01995]. The present research trajectory points toward robust mode sorters, adaptive and on-sky receivers, joint spectral-spatial schemes, and task-specific reduced measurements that preserve SPADE’s core informational advantage while relaxing its experimental demands [2409.01190], [2606.18025], [2511.20790].

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