[Paper Review] Feasibility of 3D reconstructions from a single 2D diffraction measurement

Pierre Thibault|arXiv (Cornell University)|Sep 9, 2009

Advanced X-ray Imaging Techniques3 citations

TL;DR

This paper challenges the feasibility of reconstructing 3D object densities from a single 2D diffraction measurement using ankylography, arguing that the method relies on the Born approximation, which breaks down for strongly absorbing samples. The authors demonstrate through theoretical analysis and multislice simulations that the experimental data used in the original claim do not support 3D reconstruction, rendering the reported reconstruction invalid due to invalid physical assumptions and flawed algorithmic interpretation of the Ewald sphere sampling.

ABSTRACT

We comment on the recent manuscript by Raines et al. [arXiv:0905.0269v2] (now published in Nature, vol. 463, p. 214-217, 2010), which suggests that in certain conditions a single diffraction measurement may be sufficient to reconstruct the full three-dimensional density of a scatterer. We show that past literature contains the tools to assess rigorously the feasibility of this approach. We question the formulation of the reconstruction algorithm used by the authors and we argue that the experimental data used as a demonstration is not suitable for this method, and thus that the reconstruction is not valid. This second version was produced for documentation purposes. In addition to the minimally modified original comment, it includes in appendix a subsequent reply to one of the authors (J. Miao).

Motivation & Objective

• To rigorously assess the feasibility of 3D reconstruction from a single 2D diffraction measurement using established sampling theory.
• To challenge the validity of the ankylography method proposed by Raines et al., particularly its reliance on the Born approximation.
• To demonstrate that the experimental data used in the original study are inconsistent with the assumptions required for 3D reconstruction.
• To argue that the reconstruction algorithm misinterprets the Ewald sphere sampling and fails to account for strong absorption effects.
• To provide a theoretical and simulation-based rebuttal showing that the observed data are better explained as a 2D projection with artificial 3D extrusion.

Proposed method

• Theoretical analysis of non-uniform Fourier sampling on the Ewald sphere using tools from non-uniform sampling theory and autocorrelation function support.
• Application of the Born approximation to relate measured diffraction intensities to the 3D Fourier transform of the object density.
• Use of multislice simulations to model wave propagation through a highly absorbing, tilted membrane with 3×10⁻⁴ transmittance.
• Comparison of simulated diffraction patterns with experimental data to assess whether the observed rings and intensity modulations are consistent with the Born approximation.
• Reconstruction of 3D density using phase retrieval algorithms (e.g., hybrid input-output), but with critical evaluation of their validity under non-Born conditions.
• Derivation of the physical interpretation of the Ewald sphere sampling, showing that it only represents a 2D amplitude distribution when absorption is strong.

Experimental results

Research questions

• RQ1 Can a single 2D diffraction measurement on the Ewald sphere uniquely reconstruct the 3D density of a scatterer under realistic scattering conditions?
• RQ2 Is the ankylography method valid when the Born approximation breaks down due to strong absorption in the sample?
• RQ3 Do the experimental data used in the original study support a 3D reconstruction, or are they better explained as a 2D projection with artificial depth extension?
• RQ4 Can the observed diffraction rings and intensity modulations be explained by dynamical scattering effects rather than 3D Fourier transform sampling?
• RQ5 What is the physical meaning of the reconstructed 3D array if the Born approximation is invalid and the sample is highly absorbing?

Key findings

• The reconstruction presented in the original ankylography paper is not valid because the sample's strong absorption invalidates the Born approximation, which is a prerequisite for interpreting the Ewald sphere data as a 3D Fourier transform.
• The measured intensities on the Ewald sphere do not represent a valid sampling of the 3D Fourier transform of the object's density when absorption is high, rendering the 3D reconstruction physically meaningless.
• Multislice simulations show that the observed diffraction pattern with dark rings is consistent with a 2D object with strong absorption and edge effects, not a 3D density distribution.
• The 3D array produced by the reconstruction is not the Fourier transform of the 3D electron density; instead, it is an artificial extrusion of a 2D amplitude distribution.
• The experimental data used in the original study do not support the claim of 3D reconstruction, as they are inconsistent with the theoretical framework required by ankylography.
• The use of phase retrieval algorithms like hybrid input-output does not resolve the fundamental issue of incorrect physical modeling when the Born approximation fails.