---
title: Transverse Phase Space Tomography
url: https://www.emergentmind.com/topics/transverse-phase-space-tomographic-reconstruction
type: topic
---

# Transverse Phase Space Tomography

Searching arXiv for the specified paper and closely related transverse phase-space tomography work.
Transverse phase space tomographic reconstruction is a measurement and inversion framework for recovering a phase-space distribution from a set of lower-dimensional projections acquired under controlled phase-space rotations or equivalent linear transforms. In accelerator and cold-atom settings, the reconstructed object is typically a distribution such as $f(x,x')$, $f(y,y')$, $f(x,p)$, or, in coupled problems, a higher-dimensional density over $(x,x',y,y')$ or beyond. The method is used because rms-based emittance or Twiss analyses can miss non-Gaussian structure, coupling, halo, and weak nonlinear distortions, whereas the reconstructed density can reveal the actual geometry of the occupied phase-space volume, including dephasing, crescent formation, X-shaped superpositions, and slice-dependent structure [1403.3432].

## 1. Definition, scope, and relation to conventional beam characterization

Transverse phase space describes the distribution of particle coordinates in a transverse plane, for example $(x,x')$ with $x$ the transverse position and $x' = dx/ds$ the angle with respect to the beam axis, and similarly for $(y,y')$ [2508.06013]. In cold-atom work, the analogous pair is $(x,p)$, where $x$ is position and $p$ momentum; the same tomographic logic applies because the measured density is a projection of a two-dimensional phase-space distribution after a controlled rotation [1403.3432].

The basic objective is to reconstruct the full distribution rather than only its second moments. This distinction is central in systems where the beam or cloud “does not produce a strictly elliptical phase space,” where profiles show non-Gaussian traits, or where coupling introduces interplane correlations that are not represented by uncoupled rms fits [1301.6765]. At FACET-II, this is stated explicitly: emittance measurement and Twiss characterization traditionally rely on fitting rms beam sizes versus quadrupole strength under a Gaussian-beam assumption, and this can be inaccurate for non-Gaussian beams; tomography instead recovers the actual phase-space density and reveals non-Gaussian features such as X-shaped superpositions [2508.06013].

The same formalism extends beyond uncoupled 2D reconstructions. FAST reported a fully coupled 4D transverse reconstruction using mixed projections generated by seven quadrupoles and 2D screen images [1811.04114]. SwissFEL ACHIP used nanofabricated wire scanners to reconstruct $f(x,x')$ and $f(y,y')$ for micrometer-scale beams [2102.08028]. FLASHForward experimentally demonstrated a 5D reconstruction of $f(x,x',y,y',t)$ by combining quadrupole-based phase-space rotation with a polarizable X-band transverse deflection structure [2505.13724]. PAL-XFEL implemented a full 6D generative reconstruction with standard accelerator elements, within which the transverse tomograms appear as marginals of the reconstructed six-dimensional density [2502.20654].

A common misconception is that tomography is merely a more elaborate quadrupole scan. The published work indicates a stricter claim: tomography is an inverse problem whose fidelity depends on how projections are generated, how the forward operator is modeled, and how incomplete or noisy measurements are regularized. This is especially evident when weak nonlinearities, dispersion, space charge, or limited-angle coverage alter the effective projection geometry [1403.3432].

## 2. Phase-space rotation, transport, and projection formation

The operational core of transverse tomography is the generation of a family of distinct projections. In the ideal harmonic case, phase-space points rotate uniformly. For cold atoms in a nearly harmonic trap, the exact harmonic mapping is
$$
x \to x\cos\theta - \frac{p}{m\omega_0}\sin\theta,\qquad
p \to p\cos\theta + m\omega_0 x\sin\theta,
$$
with $\theta(t)=\omega_0 t$, which permits tomography by interpreting time evolution as a rotation of the phase-space frame [1403.3432]. In accelerator beamlines, the corresponding mechanism is linear transport through drifts, quadrupoles, solenoids, or TDS streaking sections. For a single-quadrupole scan at FACET-II, the linear optics between the quadrupole and the screen is represented by
$$
M(k)=
\begin{pmatrix}
R_{11}(k) & R_{12}(k)\\
R_{21}(k) & R_{22}(k)
\end{pmatrix},
$$
and each setting produces a distinct projection labeled by
$$
\theta = \tan^{-1}\!\left(\frac{M_{1,1}(k)}{M_{1,2}(k)}\right),\qquad
s = \sqrt{M_{1,1}^2(k)+M_{1,2}^2(k)}.
$$
The measured line profiles are then rescaled and labeled by $\theta$ and $s$ to form a sinogram [2508.06013].

The projection itself is a Radon-type marginal. In the cold-atom formulation,
$$
P_{\theta}(x) = \int_{-\infty}^{\infty} f\big(x\cos\theta + p\sin\theta,\; -x\sin\theta + p\cos\theta\big)\, dp,
$$
while in FACET-II the idealized projection is written as
$$
\mathcal{R}[f](\theta,s)=\int f(x,x')\,\delta\!\big(s-x\cos\theta-x'\sin\theta\big)\,dx\,dx'.
$$
FAST expresses the same idea in 4D: for optics setting $k$, the 2D screen image is a projection of the upstream 4D density along the unmeasured angles $(x',y')$ under a transport matrix $R_k$ [1403.3432; 2508.06013; 1811.04114].

Several measurement schemes instantiate this principle differently. Fermilab’s 400 MeV Linac used an upstream quadrupole focusing triplet and a 10 m, dispersion-free, magnet-free straight with three multiwire profile monitors; different triplet settings altered the projection direction through the effective transport matrix $M(s)$, with the measured coordinate
$$
x = m_{11}(s)x_0 + m_{12}(s)x_0',
$$
and the corresponding normal angle
$$
\theta(s)=\arctan\!\big(m_{12}(s)/m_{11}(s)\big).
$$
Tomography was described as “rotating the phase space distribution using different waist focusing conditions of the upstream triplet and performing a de-convolution of the profile data” [1301.6765].

SwissFEL ACHIP replaced screens with nanofabricated wire projections. Nine free-standing gold stripes, each $1\,\mu\mathrm{m}$ wide and arranged radially and uniformly in angle, were scanned across the beam at six longitudinal positions around the waist. The measured beam loss monitor signal is proportional to the line integral of the beam density along the wire, yielding 54 projections in total. Near the waist, small changes in longitudinal position produce comparatively large changes in phase advance, which makes that region especially informative for tomography [2102.08028].

In higher-dimensional tomography, projection generation combines more than one rotation mechanism. The 5D method proposed for ARES and later demonstrated at FLASHForward used transverse phase advances $(\mu_x,\mu_y)$ from quadrupole scans together with adjustable streaking angle $\psi$ or $\theta$ from a polarizable X-band TDS. Near RF zero-crossing, the TDS maps time into an angular kick, and after downstream drift this becomes a screen shear proportional to time, enabling reconstruction of 3D $(x,y,t)$ charge density before slice-wise 4D inversion [2305.03538; 2505.13724].

## 3. Inversion algorithms and reconstruction representations

The inversion stage converts a set of measured projections into a phase-space density. Different experiments use different representations and regularization strategies, but all solve an ill-conditioned inverse problem.

In the cold-atom work, the reconstruction uses the Leonhardt kernel inversion rather than standard filtered back-projection:
$$
P(\bar{q},\bar{p}) = \frac{1}{2\pi^2}\int_0^{\pi} d\theta \int_{-\infty}^{\infty} dx\; K(\bar{q}\cos\theta + \bar{p}\sin\theta - x)\; {\rm pr}(x,\theta),
$$
with a regularized, nondivergent kernel $K(x)$ and cutoff $k_c = 0.43\,\mu\mathrm{m}^{-1}$, which sets the phase-space resolution and stabilizes inversion against noise and finite sampling [1403.3432]. The same paper notes that a maximum-likelihood approach was also tested and yielded similar reconstructions.

FACET-II uses Maximum Likelihood Expectation Maximization (MLEM) in a discretized non-negative grid model,
$$
g_i \approx \sum_j A_{ij}f_j,\qquad
f_j^{(n+1)} = f_j^{(n)} \; \frac{ \sum_i A_{ij} \, \frac{g_i}{\sum_{l} A_{il} f_l^{(n)} } }{ \sum_i A_{ij} },
$$
with the measured projections stacked into a sinogram. The paper states that MLEM is used rather than filtered back-projection and that regularization such as early stopping or smoothing can be applied, though the details are not specified there [2508.06013].

FAST formulates the coupled 4D reconstruction as a large linear system $Af=b$ over a 6×6×6×6 grid, giving $6^4=1296$ voxel amplitudes as unknowns. To suppress ringing and keep projections analytically tractable, each voxel basis function is a multivariate normal in a normalized 4D phase-space basis constructed via Williamson normalization from an initial second-moment estimate. The inversion uses truncated singular value decomposition, retaining 200 singular values, with a 1296×40659 linear system whose SVD took approximately one hour on a single Intel Xeon E5-1650 (3.5 GHz) core [1811.04114].

SwissFEL ACHIP departs from grid-based inversion and represents the distribution by an ensemble of macro-particles. Each particle carries $(x_i,x'_i,y_i,y'_i)$, and a smooth density is generated by separable Gaussian kernels with typical choices $\rho_{x,y}\approx 80\,\mathrm{nm}$ and $\rho_{x',y'}=\rho_{x,y}/z_{\max}$. The reconstruction iterates a stochastic birth–death update based on the normalized projection residuals,
$$
\Delta_{z,\theta}(\xi)=\frac{P^m_{z,\theta}(\xi)-P^r_{z,\theta}(\xi)}{\max_\xi P^r_{z,\theta}(\xi)},
$$
and stops when
$$
\frac{|p_k-p_{k-1}|}{|p_k|}<\tau,\qquad \tau\approx 0.005,
$$
which produced stable convergence in approximately 100–110 iterations for the reported data [2102.08028].

The 5D ARES and FLASHForward reconstructions use Simultaneous Algebraic Reconstruction Technique (SART), implemented slice-wise. In the 5D method, the discrete forward model is
$$
g=Af+\epsilon,
$$
with positivity and charge-related constraints. The authors report that SART converges to good results in 1–2 iterations; two iterations were used in both the simulation study and the later FLASHForward experiment [2305.03538; 2505.13724].

The 2026 resolution framework introduces a different level of analysis. It states that the experimentally accessible degrees of freedom are determined by the Gram or sampling operator $G$, with eigenvalues that quantify which Hilbert-space modes are transmitted reliably and which are weakly supported. In that framework, reconstruction in the Gram eigenbasis is a measurement-adapted compression of the tomographic problem, and features dominated by small-$\lambda_k$ modes should be treated as artefacts or as hypotheses requiring additional measurements [2605.29784]. This suggests a unifying interpretation of regularization across kernel, SVD, MLEM, and algebraic methods: each, implicitly or explicitly, suppresses modes that are poorly supported by the measurement geometry.

## 4. Experimental realizations across cold atoms and accelerators

The published implementations span cold atoms, proton linacs, photoinjectors, FEL beamlines, and GeV-class electron transport. The measurement geometries differ, but the recurring structure is projection acquisition under known transport plus inversion under explicit regularization.

| Platform | Projection generation | Reconstruction |
|---|---|---|
| Atom-chip cold atoms | Harmonic rotation with 13 images over $[0,\pi]$ | Leonhardt kernel [1403.3432] |
| Fermilab FAST | Seven quadrupoles, 22 screen projections | Truncated SVD in 4D [1811.04114] |
| SwissFEL ACHIP | Nine wire angles at six $z$ positions, 54 projections | Particle-based birth–death inversion [2102.08028] |
| FACET-II | Single-quadrupole scan upstream of a screen | MLEM sinogram inversion [2508.06013] |
| FLASHForward | Phase-advance scans plus ten TDS angles | 3D + slice-wise 4D SART [2505.13724] |

The cold-atom demonstration remains notable because it isolates deformation by weak anharmonicity with unusually clear control. The experiment used approximately 3,000 ${}^{87}\mathrm{Rb}$ atoms at approximately $160\,\mathrm{nK}$ in a magnetic atom-chip trap with longitudinal harmonic frequency $\omega_0 = 2\pi\times 38\,\mathrm{Hz}$ and transverse frequencies approximately $2\pi\times 110\,\mathrm{Hz}$. For an initial oscillation amplitude of $85\,\mu\mathrm{m}$ and corrugation wire current of $5\,\mathrm{mA}$, slight deformation appears at $t_1=100\,\mathrm{ms}$, becomes distinct at $300\,\mathrm{ms}$, and becomes a pronounced crescent-shaped distribution at $500\,\mathrm{ms}$. Without corrugation, the distribution remains approximately isotropic Gaussian up to $500\,\mathrm{ms}$ [1403.3432].

Fermilab’s 400 MeV proton Linac provided an earlier beamline realization motivated by the statement that the beam “does not produce a strictly elliptical phase space.” The 10 m straight used three multiwire profile monitors, with MW5 at the waist having $0.5\,\mathrm{mm}$ pitch, and preliminary rms analyses reported $\epsilon_y \approx 2.00\,\pi\,\mathrm{mm\cdot mrad}$ by the analytical three-screen method and $\epsilon_x \approx 1.45\,\pi\,\mathrm{mm\cdot mrad}$, with notably larger approximate 95% emittances of approximately $18\pi\,\mathrm{mm\cdot mrad}$ in the vertical plane and approximately $13\pi\,\mathrm{mm\cdot mrad}$ in the horizontal plane [1301.6765]. The tomography program there was explicitly intended to recover the non-elliptical distribution hidden behind those rms summaries.

FAST moved to a coupled 4D reconstruction at 32 MeV and 100 nC bunch charge. Seven of eight quadrupoles were varied so as to rotate the beam’s transverse phase space by $\pi$ in each plane while keeping the Twiss functions at the screen approximately constant, producing 22 distinct projections and 40,659 informative pixel values after image preparation. The reconstructed 4D density revealed small but nonzero interplane correlations consistent with solenoidal coupling and possible laser spot asymmetries or space-charge distortions [1811.04114].

SwissFEL ACHIP addressed the opposite scale regime: micrometer beams requiring submicrometer diagnostics. The reported experiment used a beam energy of approximately 3.2 GeV and bunch charge of approximately 1 pC, with scans at six $z$-locations spanning approximately 8 cm around the waist. Typical SNR ranged from 25 to 45, charge fluctuations were approximately 1.3% rms, and PMT readout noise was less than 1%. The reconstructed waist beam sizes were $\sigma_x^* \approx 1.04(6)\,\mu\mathrm{m}$ and $\sigma_y^* \approx 1.26(5)\,\mu\mathrm{m}$, with $\beta^* \approx 3.7(2)\,\mathrm{cm}$ in both planes; normalized emittances from Gaussian fits to the reconstructed phase spaces were $\epsilon_{n,x}\approx 186(15)\,\mathrm{nm\cdot rad}$ and $\epsilon_{n,y}\approx 278(18)\,\mathrm{nm\cdot rad}$ [2102.08028].

FACET-II demonstrated that a single quadrupole and a downstream screen can suffice for useful transverse tomography. In nominal two-bunch operation, with a 1.6 nC driver and 0.5 nC witness at about 7 ps separation, injector tomography revealed a distinctly non-Gaussian X-shaped structure in the horizontal phase space when the two bunches were combined. The interpretation was verified by characterizing each bunch individually. With zero pulse separation and charge varied from approximately 0.5 nC to 2 nC, both planes showed increasingly stretched distributions with higher charge, and in the horizontal plane the phase space was observed to “flip” when the bunch charge was halved [2508.06013].

FLASHForward demonstrated a 5D experimental reconstruction at 1.09 GeV and 297 pC using a PolariX TDS and a GAGG:Ce screen. Ten steps over approximately $180^\circ$ in each of $\mu_x$ and $\mu_y$ and ten TDS angles over approximately $180^\circ$ yielded 960 successfully recorded projections over 28 hours. The tomography gave an rms bunch duration of $199^{+8}_{-7}\,\mathrm{fs}$, in excellent agreement with $194\pm 16\,\mathrm{fs}$ from independent streaked-screen measurements, and projected normalized emittances of $\epsilon_{n,x}\approx 6.63\pm 0.10\,\mu\mathrm{m}$ and $\epsilon_{n,y}\approx 8.25\pm 0.20\,\mu\mathrm{m}$ [2505.13724].

## 5. Nonlinearity, coupling, multi-bunch structure, and higher-dimensional generalizations

One of the clearest lessons from the literature is that tomography is unusually sensitive to weak distortions that are difficult to detect in single projections. In the atom-chip experiment, weak static corrugation $\delta V(x)$ produced energy-dependent angular velocity,
$$
\dot{\theta} = -\omega_0\left[1 - \frac{x\, \delta F(x)/2}{E - \delta V(x)}\right],
$$
with $\delta F(x)=-\partial_x\delta V(x)$, making the effective rotation frequency $\Omega(E)$ energy dependent and causing angular-velocity dispersion across the ensemble. For $x_{\rm shift}=85\,\mu\mathrm{m}$ and $T=160\,\mathrm{nK}$, the ensemble spans $\delta\omega/\omega_0 \approx -0.015$ to $+0.009$, giving an angular dispersion after $t_1=500\,\mathrm{ms}$ of $\Delta\theta\sim 2.85\,\mathrm{rad}$, which explains the observed crescent-shaped reconstructions [1403.3432].

The accelerator analogue is amplitude- or energy-dependent distortion from field errors, fringe fields, misalignments, dispersion, or space charge. FAST frames the coupled 4D problem in terms of linear transport and symplectic invariants such as eigenemittances, showing that mixed projections can reveal correlations that uncoupled fits miss [1811.04114]. FACET-II shows a different class of departure from rms optics: superimposed bunches with different transverse phase spaces can form an X-shaped composite distribution, and charge tuning can reveal space-charge-dominated dynamics that are not evident in standard Gaussian-based emittance extraction [2508.06013].

The higher-dimensional extensions make these couplings explicit rather than treating them as nuisances. The 5D method reconstructs
$$
z=(x,x',y,y',t)^{\mathrm T},\qquad f(z)\ge 0,
$$
by combining quadrupole-based transverse rotations with a TDS whose polarization angle $\psi$ is adjustable over approximately $180^\circ$. In normalized coordinates, the transport in each plane is a pure rotation,
$$
\begin{pmatrix}x_N\\x_N'\end{pmatrix}_f=
\begin{pmatrix}\cos\mu_x & \sin\mu_x\\-\sin\mu_x & \cos\mu_x\end{pmatrix}
\begin{pmatrix}x_N\\x_N'\end{pmatrix}_i,
$$
and similarly in $y$, which is why slice-wise 4D inversion becomes possible [2305.03538]. FLASHForward then used the reconstructed 5D density to generate a particle distribution for simulations and to extract the transverse 4D slice emittance $\epsilon_{4D}(z)=\sqrt{\det\Sigma_\perp(z)}$, finding an average relative discrepancy of approximately 5% between $\epsilon_{4D}(t)$ and $\epsilon_x(t)\epsilon_y(t)$, with a minimum of approximately 1% near $t\approx -170\,\mathrm{fs}$, indicating minor residual cross-plane correlations in the core [2505.13724].

The 6D extensions pursue a still broader goal: using transverse images under multiple machine settings to infer the full source distribution over $(x,x',y,y',z,\delta)$. The two-stage CNN trained on KEK-ATF ASTRA data reconstructs 15 pairwise 2D histograms from only sixteen $x$–$y$ screen images measured at a dispersive location under different RF phases and solenoid fields, with inference in under a minute on an NVIDIA RTX A400 GPU [2603.02733]. PAL-XFEL’s generative phase space reconstruction instead uses a transformer-based generative neural network combined with differentiable Bmad-X tracking and standard accelerator elements, including a quadrupole, bunch compressor, X-band linearizer, and screens, to reconstruct a near-unique 6D state whose downstream predictions match withheld measurements [2502.20654].

These higher-dimensional works do not replace transverse tomography so much as subsume it. The published descriptions state directly that within a full 6D reconstruction, transverse tomography corresponds to the marginal distributions over $(x,x')$ and $(y,y')$ extracted from the reconstructed 6D density [2502.20654]. This suggests that the boundary between “transverse tomography” and “full phase-space diagnostics” is increasingly methodological rather than conceptual.

## 6. Resolution limits, artefacts, validation, and methodological outlook

Tomography is limited not only by detector resolution and shot noise but by the measurement geometry itself. The most explicit formulation of this point is given by the Gram-operator framework, which states that the measurement induces a sampling operator $G$ whose eigenvalues determine the experimentally accessible modes and the reconstruction bandwidth. Large eigenvalues correspond to reliably resolved modes; small eigenvalues correspond to weakly accessible modes that are highly sensitive to noise. Features supported predominantly by modes below a threshold $\lambda_{\mathrm{thr}}$ should be treated as artefacts or as hypotheses requiring additional measurements [2605.29784].

In practical accelerator tomography, the same limitation appears as finite-angle conditioning, sparse projection sets, or insufficiently diverse optics. FAST notes that 22 images, corresponding to approximately $\pi$ coverage in each plane, were sufficient for a stable 4D inversion on a $6^4$ grid, but also remarks that more angles or mixed-coupling optics can improve conditioning, especially for finer grids [1811.04114]. The 5D ARES study quantified this dependence: to reach less than or equal to 5% discrepancies in all reconstructed transverse planes, approximately 50 transverse rotation angles and 50 streaking angles were needed; for less than or equal to 10% discrepancies, approximately 30 rotation angles and approximately 25 streaking angles sufficed [2305.03538].

Detector and instrument models matter comparably. Fermilab’s multiwire data require deconvolution of a point-spread function set primarily by wire pitch, wire diameter, mechanical alignment, and electronic response [1301.6765]. ACHIP explicitly convolves the forward model with the measured wire profile of width $w=1\,\mu\mathrm{m}$ and reports that spatial resolution is ultimately limited by wire roughness at approximately 100 nm [2102.08028]. FLASHForward propagates uncertainty from the measured shear parameters by performing 100 reconstructions with $\eta_\theta$ drawn from Gaussian distributions and notes that a 1% beam-energy error leads to approximately 5% emittance error [2505.13724].

Validation strategies across the literature are correspondingly stringent. The atom-chip work compares reconstructions against simulations using both energy-dependent phase angles and full numerical integration of Newton’s equations including rare collisions, finding quantitative agreement and phase-space area conservation when dimensions are separable and collisions are infrequent [1403.3432]. FAST compares measured profiles for the $x$-scan and $y$-scan with corresponding projections of the reconstructed model and propagates the reconstructed second moments to the screen using $\Sigma_{\mathrm{out},k}=R_k\Sigma_{\mathrm{in}}R_k^{\mathrm T}$ [1811.04114]. FLASHForward tracks a five-million-particle distribution sampled from the reconstructed 5D density to the screen for all measured settings and reports a centroid-based discrepancy of $32\%\pm 15\%$ across projections, with qualitative agreement of the main features [2505.13724]. PAL-XFEL validates its 6D reconstruction by predicting downstream longitudinal phase space and $(y,p_z)$ correlations that were withheld from training, supporting the claim of near-unique reconstruction [2502.20654].

The field also contains an objective tension between model-based and learned reconstructions. Kernel inversions, SVD, MLEM, and SART make the forward model explicit and their regularization comparatively transparent [1403.3432; 1811.04114; 2508.06013; 2505.13724]. CNN and generative methods reduce measurement time and can incorporate rich coupled physics, but they depend on the representativeness of simulation data, the fidelity of the transport model, and calibration stability [2603.02733; 2502.20654]. A plausible implication is that future transverse phase-space tomography will increasingly combine both approaches: measurement-adapted subspaces or Gram-eigenmode criteria for what is genuinely resolved, together with flexible learned priors for interpolation within that accessible bandwidth [2605.29784].

Across cold atoms, linacs, FEL injectors, and GeV-class beamlines, the topic is now defined less by a single inversion formula than by a shared architecture: generate controlled phase-space projections, model the transport operator faithfully, regularize according to measurement support, and validate by forward prediction. Within that architecture, transverse phase space tomographic reconstruction has become a diagnostic for weak anharmonicity, 4D coupling, multi-bunch superposition, slice emittance, and full 5D or 6D phase-space structure [1403.3432].

Source: https://www.emergentmind.com/topics/transverse-phase-space-tomographic-reconstruction