Fourier Plane Tomographic Spectroscopy

• Fourier Plane Tomographic Spectroscopy is a technique that fuses angle-resolved spectroscopy with tomographic reconstruction to capture momentum and wavelength information in the Fourier domain.
• It employs optical configurations like back-focal-plane microscopy and rotational slitless spectroscopy to obtain four-dimensional scattering data with high angular and spectral resolution.
• The method addresses complex inverse problems using approaches such as filtered backpropagation and gradient-based optimization to recover structural and modal details from incomplete measurements.

Fourier Plane Tomographic Spectroscopy denotes a class of measurement and reconstruction strategies that combine spectroscopic information with tomography in a Fourier-domain representation. In back-focal-plane microscopy and spectroscopy, it is an angle-resolved, wavelength-resolved measurement carried out at the back focal plane (BFP, or Fourier plane) of a high-NA microscope objective, where the primary observable is the momentum-resolved light distribution $I(k_x,k_y,\lambda)$ (Vasista et al., 2018). In recent single-particle scattering work, the same concept has been extended to simultaneous four-dimensional characterization of scattering as a function of wavelength, incident direction, scattering direction, and polarization, yielding a scattering signature $S(\lambda,k_i,k_s,p)$ (Patzschke et al., 21 Jul 2025). Closely related diffraction-tomographic formulations use Fourier-domain mappings on Ewald-sphere manifolds and explicit filtered backpropagation formulas to reconstruct scattering potentials from wave measurements (Kirisits et al., 2024). A distinct, but terminologically adjacent, usage appears in standoff FPA-FTIR plume sensing, where the Fourier operation is the interferogram-to-spectrum transform and the tomographic output is assimilated as a concentration contour in a PDE-constrained inverse problem (Mattuschka et al., 10 Jun 2026). Rotational slitless spectroscopy extends the same tomography–spectroscopy logic to wide-field datacube recovery through repeated angular projections and matrix inversion (Zhang et al., 22 May 2026).

1. Terminology and conceptual scope

In the BFP-based literature, the Fourier plane is the objective pupil plane that encodes angular or wavevector information. Each point in the pupil corresponds to a unique emission or scattering direction, so a BFP image is a k-space image rather than a real-space image. On that basis, Fourier Plane Tomographic Spectroscopy measures the angular spectrum of light together with wavelength, and uses forward modeling or inversion to infer structural, orientational, or modal properties of the source or scatterer (Vasista et al., 2018).

The Janus-particle implementation makes this scope explicit by defining FPTS as a BFP-based method that resolves scattering simultaneously as a function of wavelength, incident direction, scattering direction, and polarization. In that setting, tomography is performed over illumination angles: the object is probed from many $k_i$, and the scattered field is recorded over the full $k_s$ manifold allowed by the objective NA (Patzschke et al., 21 Jul 2025).

The term is not uniform across all subfields. In dual FPA-FTIR standoff sensing, “FPA-FTIR” refers to focal plane array Fourier transform infrared spectroscopy, and the “Fourier transform” is in the spectral domain rather than the optical Fourier plane in spatial frequency. The associated inverse method consumes a tomographically reconstructed threshold contour rather than raw pupil-plane data (Mattuschka et al., 10 Jun 2026). This distinction matters because it separates two non-identical uses of “Fourier” in tomography: Fourier-plane angular spectroscopy in microscopy and Fourier-transform spectral reconstruction in standoff hyperspectral sensing.

A related extension appears in rotational slitless spectroscopy. There, multiple slitless projections acquired at different rotation angles are interpreted tomographically, and the reconstruction recovers a three-dimensional integral-field datacube. The governing language is that of the Radon transform and Fourier slice theorem rather than BFP microscopy, but the operational goal is similar: use spectrally encoded projections to reconstruct a higher-dimensional object (Zhang et al., 22 May 2026).

2. Fourier-plane encoding and k-space representations

The core optical principle in BFP methods is the one-to-one mapping between angle and pupil position. For a medium of refractive index $n$ and free-space wavelength $\lambda_0$, the wavevector magnitude is $k = 2\pi n/\lambda_0$, with components

$$k_x = k\sin\theta\cos\phi,\qquad k_y = k\sin\theta\sin\phi,\qquad k_z = k\cos\theta.$$

For an aplanatic high-NA objective, the pupil radius obeys

$$\rho_{\mathrm{BFP}} = f_{\mathrm{obj}}\, n \sin\theta,$$

and, after pupil magnification $M_p$, the camera-plane radius is $S(\lambda,k_i,k_s,p)$0. Using the pupil edge calibration $S(\lambda,k_i,k_s,p)$1, one obtains

$S(\lambda,k_i,k_s,p)$2

This gives a direct pixel-to-$S(\lambda,k_i,k_s,p)$3 conversion and defines the collection cone imposed by the objective NA (Vasista et al., 2018).

The Janus-particle formulation uses the same mapping in scattering notation. In that work, each BFP point $S(\lambda,k_i,k_s,p)$4 corresponds to a unique scattered direction $S(\lambda,k_i,k_s,p)$5, with

$S(\lambda,k_i,k_s,p)$6

and $S(\lambda,k_i,k_s,p)$7. The BFP image is spectrally dispersed along one camera axis, so that one detector coordinate carries wavelength and the orthogonal coordinate carries a fixed in-plane BFP coordinate. Scanning a slit across the BFP reconstructs the full angular distribution versus $S(\lambda,k_i,k_s,p)$8 (Patzschke et al., 21 Jul 2025).

Diffraction tomography expresses the same k-space logic through the generalized Fourier diffraction theorem. For plane-wave incidence $S(\lambda,k_i,k_s,p)$9 and detector plane $k_i$0 outside the object support, the partial Fourier transform of the scattered field satisfies

$k_i$1

with $k_i$2 and $k_i$3 for $k_i$4. The accessible object-space sample is therefore

$k_i$5

which lies on a hemisphere of the Ewald sphere centered at $k_i$6 (Kirisits et al., 2024).

These formulations differ in physical setting, but they share the same structural feature: the detector does not measure the object directly in real space. It measures either a momentum distribution or a wavefield whose Fourier-domain geometry determines which parts of object k-space are accessible.

3. Tomographic inverse models and reconstruction formalisms

The inverse problem in diffraction tomography is organized around Fourier coverage. With acquisition parameters collected in

$k_i$7

the k-space mapping is

$k_i$8

and the coverage set is

$k_i$9

The resulting filtered backpropagation reconstructs the $k_s$0-best approximation whose Fourier support lies in $k_s$1, with explicit Jacobian and multiplicity corrections through $k_s$2 and the Banach indicatrix $k_s$3 (Kirisits et al., 2024). In that framework, missing-cone structure, finite aperture, and limited angular diversity are properties of the coverage set itself rather than merely numerical artifacts.

Rotational slitless spectroscopy casts reconstruction as a linear inverse problem. After discretizing the datacube $k_s$4 and stacking all detector images into $k_s$5, the forward model is

$k_s$6

The least-squares formulation minimizes $k_s$7, with normal equations $k_s$8. The paper emphasizes on-the-fly forward and adjoint operators $k_s$9 and $n$0, gradient descent, and Lucy–Richardson updates,

$n$1

rather than explicit matrix storage. Tomographically, each rotated projection corresponds to a different family of plane integrals through the $n$2-cube, and the Fourier slice theorem links the 2D Fourier transform of each projection to a central plane in the 3D Fourier transform of the datacube (Zhang et al., 22 May 2026).

The standoff FTIR source-localization framework uses a different inverse model again. The contaminant concentration $n$3 evolves under the advection–diffusion PDE

$n$4

with initial condition $n$5. The tomographic measurement is a threshold contour

$n$6

and the inverse problem is written over nonnegative Radon measures,

$n$7

The optimization is solved by the Primal-Dual Active Point (PDAP) algorithm, with forward PDE solves, adjoint solves driven by a line delta source on $n$8, greedy activation by a dual certificate, and convex amplitude updates on the active set (Mattuschka et al., 10 Jun 2026).

A plausible implication is that “tomographic spectroscopy” in current usage is less a single inversion algorithm than a family of inverse problems distinguished by the geometry of the measurement operator: Ewald-sphere sampling in diffraction tomography, Radon-type projections in rotational slitless spectroscopy, and level-set constraints in standoff plume reconstruction.

4. Acquisition architectures and optical implementations

Classical BFP microscopy requires relay optics because the objective’s BFP lies inside the objective body. The review literature describes three common configurations for projecting a pupil conjugate to the detector: a Bertrand lens before the image plane, a Bertrand lens after the image plane, and a modified $n$9 relay. Telecentric pupil relays are emphasized because they stabilize the mapping against sample axial motion. For spectroscopy, a slit can be placed at a pupil conjugate and coupled to an imaging spectrograph, or the BFP can be projected directly onto the spectrograph entrance slit (Vasista et al., 2018).

The Janus-particle implementation realizes this architecture in a dark-field scattering microscope. Illumination uses a dark-field condenser with Fourier-plane apertures $\lambda_0$0 and $\lambda_0$1 to restrict incident wavevectors to a narrow cone of approximately $\lambda_0$2; detection uses an oil-immersion objective with adjustable NA up to $\lambda_0$3, a relay system that images either the sample plane or the BFP, a $\lambda_0$4 spatial pinhole $\lambda_0$5 to isolate a single particle, a blazed transmission grating in front of an sCMOS camera, and an adjustable slit $\lambda_0$6 that is translated to reconstruct the full BFP at each wavelength. Spectral response $\lambda_0$7 is obtained from the spectrally dispersed illumination pattern, and datasets with and without a long-pass filter are merged to cover approximately $\lambda_0$8 (Patzschke et al., 21 Jul 2025).

Rotational slitless spectroscopy uses a different acquisition geometry. Independent projections are generated either by rotating the dispersion direction relative to the detector or by rotating the entire telescope or field. The dispersion law is modeled as

$\lambda_0$9

and the rotation-dependent coordinate mappings determine how each datacube voxel contributes to detector pixels. The number of required independent rotations follows the counting relation

$k = 2\pi n/\lambda_0$0

which directly links instrument geometry, chosen pixelization, and observing strategy (Zhang et al., 22 May 2026).

The standoff FTIR configuration uses dual FPA-FTIR hyperspectral imaging systems operated with a suitable opening angle so that their fields of view overlap the suspected release area. Each system acquires interferograms that are Fourier transformed into calibrated radiance spectra; tomographic processing then yields spatial maps and, for the inverse model, a threshold contour $k = 2\pi n/\lambda_0$1. The paper emphasizes that representing the measurement as a contour rather than a volumetric field allows a mesh-independent formulation via a curve discretization and a trace operator, avoiding remeshing (Mattuschka et al., 10 Jun 2026).

5. Applications and demonstrated regimes

BFP imaging and spectroscopy have been used to study angular emission patterns of fluorescence and Raman signals from molecules, elastic scattering from nanostructures, energy–momentum dispersion of guided or leaky modes, metasurface emission, and polarization-resolved momentum-space observables. The review specifically highlights single-molecule dipole orientation retrieval, Raman radiation patterns in graphene, plasmonic nanorod and nanowire scattering, Yagi–Uda antenna emission, and E–k mapping in photonic crystals and microcavities (Vasista et al., 2018).

The Janus-particle study gives a detailed contemporary example of FPTS as four-dimensional scattering tomography. The particles are polystyrene spheres of diameter $k = 2\pi n/\lambda_0$2 with hemispherical Au caps, immersed in oil of refractive index approximately $k = 2\pi n/\lambda_0$3. The measurements and finite-element simulations identify three distinct multipolar families up to fifth order: axial-propagating transverse-electric, transverse-propagating transverse-electric, and transverse-propagating axial-electric. A shoulder at $k = 2\pi n/\lambda_0$4 is prominent for Au-side illumination, weaker for side-on illumination, and largely absent for PS-side illumination; side-on illumination shows a cluster of weaker peaks between $k = 2\pi n/\lambda_0$5. The measured angular patterns serve as orientation fingerprints, and the reported spectra show progressive red-shifts and linewidth narrowing of higher-order resonances (Patzschke et al., 21 Jul 2025).

In standoff plume reconstruction, the synthetic test case uses kinematic viscosity $k = 2\pi n/\lambda_0$6, inflow velocity $k = 2\pi n/\lambda_0$7 at the southern boundary, diffusion coefficient $k = 2\pi n/\lambda_0$8, measurement time $k = 2\pi n/\lambda_0$9, and threshold $$k_x = k\sin\theta\cos\phi,\qquad k_y = k\sin\theta\sin\phi,\qquad k_z = k\cos\theta.$$0. PDAP identifies a sparse set of 6 active Dirac components; the post-processed source achieves a localization error of approximately $$k_x = k\sin\theta\cos\phi,\qquad k_y = k\sin\theta\sin\phi,\qquad k_z = k\cos\theta.$$1, total reconstructed intensity $$k_x = k\sin\theta\cos\phi,\qquad k_y = k\sin\theta\sin\phi,\qquad k_z = k\cos\theta.$$2 for ground truth $$k_x = k\sin\theta\cos\phi,\qquad k_y = k\sin\theta\sin\phi,\qquad k_z = k\cos\theta.$$3, reconstruction error below $$k_x = k\sin\theta\cos\phi,\qquad k_y = k\sin\theta\sin\phi,\qquad k_z = k\cos\theta.$$4 across the domain at measurement and forecast times, and convergence in 16 forward and 16 adjoint solves (Mattuschka et al., 10 Jun 2026).

In rotational slitless spectroscopy, numerical experiments generate toy datacubes containing 100 stars with uniform spectra, sampled by $$k_x = k\sin\theta\cos\phi,\qquad k_y = k\sin\theta\sin\phi,\qquad k_z = k\cos\theta.$$5 photons per star. With datacube size $$k_x = k\sin\theta\cos\phi,\qquad k_y = k\sin\theta\sin\phi,\qquad k_z = k\cos\theta.$$6, detector size $$k_x = k\sin\theta\cos\phi,\qquad k_y = k\sin\theta\sin\phi,\qquad k_z = k\cos\theta.$$7, and $$k_x = k\sin\theta\cos\phi,\qquad k_y = k\sin\theta\sin\phi,\qquad k_z = k\cos\theta.$$8 rotations, both rotation modes reconstruct the datacube accurately. Flux profiles along $$k_x = k\sin\theta\cos\phi,\qquad k_y = k\sin\theta\sin\phi,\qquad k_z = k\cos\theta.$$9, $$\rho_{\mathrm{BFP}} = f_{\mathrm{obj}}\, n \sin\theta,$$0, and $$\rho_{\mathrm{BFP}} = f_{\mathrm{obj}}\, n \sin\theta,$$1 agree with ground truth to within approximately $$\rho_{\mathrm{BFP}} = f_{\mathrm{obj}}\, n \sin\theta,$$2 across the central region, with overall deviations $$\rho_{\mathrm{BFP}} = f_{\mathrm{obj}}\, n \sin\theta,$$3 at worst; 500 iterations run in less than $$\rho_{\mathrm{BFP}} = f_{\mathrm{obj}}\, n \sin\theta,$$4 on a 10-core M2 laptop with $$\rho_{\mathrm{BFP}} = f_{\mathrm{obj}}\, n \sin\theta,$$5 RAM (Zhang et al., 22 May 2026).

6. Limitations, ambiguities, and future development

A persistent ambiguity is terminological. In optical microscopy, the Fourier plane is the BFP of an objective and directly encodes wavevector content. In FPA-FTIR plume sensing, the Fourier transform is performed on interferograms to obtain spectra, and the tomographic inversion is driven by a geometric contour rather than by optical Fourier-plane imaging (Mattuschka et al., 10 Jun 2026). The two usages are compatible only at a high level of abstraction: both combine spectroscopy with tomography, but they do so through different physical measurement operators.

Coverage limitations are central in every branch of the subject. In diffraction tomography, the generalized Fourier diffraction theorem is derived in the scalar Helmholtz model under Born or Rytov approximations, and finite aperture or limited angle diversity produces missing-cone structure and reduced low-frequency coverage. The coverage remains bounded by radius $$\rho_{\mathrm{BFP}} = f_{\mathrm{obj}}\, n \sin\theta,$$6, and non-uniform sampling requires Jacobian compensation and multiplicity correction (Kirisits et al., 2024).

In BFP-based FPTS on single particles, the angular range is limited by the objective NA and the incident-angle range by the condenser NA. The reported Janus-particle measurements are intensity-only; no phase is recorded, and polarization channels are not separated in the reported experiments, even though the platform supports polarization-resolved operation. Signal-to-noise degrades beyond approximately $$\rho_{\mathrm{BFP}} = f_{\mathrm{obj}}\, n \sin\theta,$$7 because of detector quantum efficiency, and sample-to-sample variability, roughness, and thickness variation broaden or weaken higher-order resonances (Patzschke et al., 21 Jul 2025).

In rotational slitless spectroscopy, ill-conditioning arises from correlated projection angles, finite detector sampling, and edge truncation. The paper mitigates boundary artifacts by using only the central $$\rho_{\mathrm{BFP}} = f_{\mathrm{obj}}\, n \sin\theta,$$8 of the detector as an active region and leaving a $$\rho_{\mathrm{BFP}} = f_{\mathrm{obj}}\, n \sin\theta,$$9 border as a buffer. The authors identify regularization, adaptive voxelization, Bayesian inference, angle-set optimization, filtered backprojection adapted to realistic sampling, and hybrid rotation modes as natural extensions (Zhang et al., 22 May 2026).

In standoff contour-based spectroscopy, the inverse problem is severely ill-posed and underdetermined without sparsity. Accuracy depends on the quality of the extracted contour, the assumed advection–diffusion model, known wind and diffusion parameters, and sufficient multi-view coverage. The paper notes that strong winds, rapidly changing conditions, sensor noise, and occlusions can make the contour incomplete or ambiguous, and it proposes future work on simultaneous calibration of turbulence, reduced-order surrogates, and noise-aware formulations based on

$M_p$0

(Mattuschka et al., 10 Jun 2026).

Taken together, these developments show that Fourier Plane Tomographic Spectroscopy is best understood as a technically heterogeneous field organized around a common principle: spectroscopic measurements are acquired in a representation tied to angular, Fourier, or projection geometry, and tomography reconstructs latent spatial, modal, or dynamical structure from incomplete but physically structured data.