No Arabic abstract
Coherent diffraction imaging (CDI) is high-resolution lensless microscopy that has been applied to image a wide range of specimens using synchrotron radiation, X-ray free electron lasers, high harmonic generation, soft X-ray laser and electrons. Despite these rapid advances, it remains a challenge to reconstruct fine features in weakly scattering objects such as biological specimens from noisy data. Here we present an effective iterative algorithm, termed oversampling smoothness (OSS), for phase retrieval of noisy diffraction intensities. OSS exploits the correlation information among the pixels or voxels in the region outside of a support in real space. By properly applying spatial frequency filters to the pixels or voxels outside the support at different stage of the iterative process (i.e. a smoothness constraint), OSS finds a balance between the hybrid input-output (HIO) and error reduction (ER) algorithms to search for a global minimum in solution space, while reducing the oscillations in the reconstruction. Both our numerical simulations with Poisson noise and experimental data from a biological cell indicate that OSS consistently outperforms the HIO, ER-HIO and noise robust (NR)-HIO algorithms at all noise levels in terms of accuracy and consistency of the reconstructions. We expect OSS to find application in the rapidly growing CDI field as well as other disciplines where phase retrieval from noisy Fourier magnitudes is needed.
In both light optics and electron optics, the amplitude of a wave scattered by an object is an observable that is usually recorded in the form of an intensity distribution in a real space image or a diffraction image. In contrast, retrieval of the phase of a scattered wave is a well-known challenge, which is usually approached by interferometric or numerical methods. In electron microscopy, as a result of constraints in the lens setup, it is particularly difficult to retrieve the phase of a diffraction image. Here, we use a defocused beam generated by a nanofabricated hologram to form a reference wave that can be interfered with a diffracted beam. This setup provides an extended interference region with the sample wavefunction in the Fraunhofer plane. As a case study, we retrieve the phase of an electron vortex beam. Beyond this specific example, the approach can be used to retrieve the wavefronts of diffracted beams from a wide range of samples.
This paper proposes a novel algorithm for image phase retrieval, i.e., for recovering complex-valued images from the amplitudes of noisy linear combinations (often the Fourier transform) of the sought complex images. The algorithm is developed using the alternating projection framework and is aimed to obtain high performance for heavily noisy (Poissonian or Gaussian) observations. The estimation of the target images is reformulated as a sparse regression, often termed sparse coding, in the complex domain. This is accomplished by learning a complex domain dictionary from the data it represents via matrix factorization with sparsity constraints on the code (i.e., the regression coefficients). Our algorithm, termed dictionary learning phase retrieval (DLPR), jointly learns the referred to dictionary and reconstructs the unknown target image. The effectiveness of DLPR is illustrated through experiments conducted on complex images, simulated and real, where it shows noticeable advantages over the state-of-the-art competitors.
The first step to gain optical control over the ultrafast processes initiated by light in solids is a correct identification of the physical mechanisms at play. Among them, exciton formation has been identified as a crucial phenomenon which deeply affects the electro-optical properties of most semiconductors and insulators of technological interest. While recent experiments based on attosecond spectroscopy techniques have demonstrated the possibility to observe the early-stage exciton dynamics, the description of the underlying exciton properties remains non-trivial. In this work we propose a new method called extended Ptychographic Iterative engine for eXcitons (ePIX), capable of reconstructing the main physical properties which determine the evolution of the quasi-particle with no prior knowledge of the exact relaxation dynamics or the pump temporal characteristics. By demonstrating its accuracy even when the exciton dynamics is comparable to the pump pulse duration, ePIX is established as a powerful approach to widen our knowledge of solid-state physics.
In this work we develop an algorithm for signal reconstruction from the magnitude of its Fourier transform in a situation where some (non-zero) parts of the sought signal are known. Although our method does not assume that the known part comprises the boundary of the sought signal, this is often the case in microscopy: a specimen is placed inside a known mask, which can be thought of as a known light source that surrounds the unknown signal. Therefore, in the past, several algorithms were suggested that solve the phase retrieval problem assuming known boundary values. Unlike our method, these methods do rely on the fact that the known part is on the boundary. Besides the reconstruction method we give an explanation of the phenomena observed in previous work: the reconstruction is much faster when there is more energy concentrated in the known part. Quite surprisingly, this can be explained using our previous results on phase retrieval with approximately known Fourier phase.
Single-particle imaging with X-ray free-electron lasers depends crucially on algorithms that merge large numbers of weak diffraction patterns despite missing measurements of parameters such as particle orientations. The Expand-Maximize-Compress (EMC) algorithm is highly effective at merging single-particle diffraction patterns with missing orientation values, but most implementations exhaustively sample the space of missing parameters and may become computationally prohibitive as the number of degrees of freedom extend beyond orientation angles. Here we describe how the EMC algorithm can be modified to employ Metropolis Monte Carlo sampling rather than grid sampling, which may be favorable for cases with more than three missing parameters. Using simulated data, this variant is compared to the standard EMC algorithm. Higher dimensional cases of mixed target species and variable x-ray fluence are also explored.