Phase retrieval from images in the presence of first-order vortices

Estimation of dislocated phases and tunable orbital angular momentum using two cylindrical lenses

A first-order optical system consisting of two cylindrical lenses separated by a distance is considered. It is found to be non-conserving of orbital angular momentum of the incoming paraxial light field. The first-order optical system is effectively demonstrated to estimate phases with dislocations using a Gerchberg-Saxton-type phase retrieval algorithm by making use of measured intensities. Tunable orbital angular momentum in the outgoing light field is experimentally demonstrated using the considered first-order optical system by varying the distance of separation between the two cylindrical lenses.

Imaging Simulation of Charged Nanowires in TEM with Large Defocus Distance

In transmission electron microscope (TEM), both the amplitude and phase of electron beam change when electrons traverse a specimen. The amplitude is easily obtained by the square root of the intensity of a TEM image, while the phase affects defocused images. In order to obtain the phase map and verify the theoretical model of the interaction between electron beam and specimen, a lot of simulations have been performed by researchers. In this work, we have simulated defocus images of a SiC nanowire in TEM with the method of electron optics. Mean inner potential and charge distribution on the nanowire have been considered in the simulation. Besides, due to electron scattering coherence loss of the electron beam has been introduced. A dynamic process with Bayesian optimization was used in the simulation. With the infocus image as input and by adjusting fitting parameters the defocus image is determined with a reasonable charge distribution. The calculated defocus images are in a good agreement with the experimental ones. Here, we present a complete solution and verification method for solving nanoscale charge distribution in TEM.

Phase-retrieval Fourier microscopy of partially temporally coherent nanoantenna radiation patterns

We report an experimental technique for determining phase-resolved radiation patterns of single nanoantennas by phase-retrieval defocused imaging. A key property of nanoantennas is their ability to imprint spatial coherence, for instance, on fluorescent sources. Yet, measuring emitted wavefronts in absence of a reference field is difficult. We realize a defocused back focal plane microscope to measure phase even for partially temporally coherent light and benchmark the method using plasmonic bullseye antenna scattering. We outline the limitations of defocused imaging which are set by spectral bandwidth and antenna mode structure. This work is a first step to resolve wavefronts from fluorescence controlled by nanoantennas.

Estimation of dislocated phases in wavefronts through intensity measurements using a Gerchberg-Saxton type algorithm

Estimation of the phase of a singular paraxial light field from experimentally measured intensities using a Gerchberg-Saxton type algorithm is demonstrated. A combination of cylindrical lenses which does not conserve the orbital angular momentum of the light field is used in obtaining the measured intensities. Consistent extraction of the phases in regard of the orbital angular momentum is demonstrated both at the input and output transverse planes, using the measured intensities.

3D Imaging Based on Depth Measurement Technologies

Three-dimensional (3D) imaging has attracted more and more interest because of its widespread applications, especially in information and life science. These techniques can be broadly divided into two types: ray-based and wavefront-based 3D imaging. Issues such as imaging quality and system complexity of these techniques limit the applications significantly, and therefore many investigations have focused on 3D imaging from depth measurements. This paper presents an overview of 3D imaging from depth measurements, and provides a summary of the connection between the ray-based and wavefront-based 3D imaging techniques.

Probe reconstruction for holographic X-ray imaging

In X-ray holographic near-field imaging the resolution and image quality depend sensitively on the beam. Artifacts are often encountered due to the strong focusing required to reach high resolution. Here, two schemes for reconstructing the complex-valued and extended wavefront of X-ray nano-probes, primarily in the planes relevant for imaging (i.e. focus, sample and detection plane), are presented and compared. Firstly, near-field ptychography is used, based on scanning a test pattern laterally as well as longitudinally along the optical axis. Secondly, any test pattern is dispensed of and the wavefront reconstructed only from data recorded for different longitudinal translations of the detector. For this purpose, an optimized multi-plane projection algorithm is presented, which can cope with the numerically very challenging setting of a divergent wavefront emanating from a hard X-ray nanoprobe. The results of both schemes are in very good agreement. The probe retrieval can be used as a tool for optics alignment, in particular at X-ray nanoprobe beamlines. Combining probe retrieval and object reconstruction is also shown to improve the image quality of holographic near-field imaging.

Capturing and visualizing transient X-ray wavefront topological features by single-grid phase imaging

The detection, localisation and characterisation of stationary and singular points in the phase of an X-ray wavefield is a challenge, particularly given a time-evolving field. In this paper, the associated difficulties are met by the single-grid, single-exposure X-ray phase contrast imaging technique, enabling direct measurement of phase maxima, minima, saddle points and vortices, in both slowly varying fields and as a means to visualise weakly-attenuating samples that introduce detailed phase variations to the X-ray wavefield. We examine how these high-resolution vector measurements can be visualised, using branch cuts in the phase gradient angle to characterise phase features. The phase gradient angle is proposed as a useful modality for the localisation and tracking of sample features and the magnitude of the phase gradient for improved visualization of samples in projection, capturing edges and bulk structure while avoiding a directional bias. In addition, we describe an advanced two-stage approach to single-grid phase retrieval.

Transport of intensity phase imaging in the presence of curl effects induced by strongly absorbing photomasks

We report theoretical and experimental results for imaging of electromagnetic phase edge effects in lithography photomasks. Our method starts from the transport of intensity equation (TIE), which solves for phase from through-focus intensity images. Traditional TIE algorithms make an implicit assumption that the underlying in-plane power flow is curl-free. Motivated by our current study, we describe a practical situation in which this assumption breaks down. Strong absorption gradients in mask features interact with phase edges to contribute a curl to the in-plane Poynting vector, causing severe artifacts in the phase recovered. We derive how curl effects are coupled into intensity measurements and propose an iterative algorithm that not only corrects the artifacts, but also recovers missing curl components.

Phase imaging for absorptive phase objects using hybrid uniform and structured illumination transport of intensity equation

Transport of intensity equation (TIE) has been a popular and convenient phase imaging method that retrieves phase profile from the measurement of intensity differentials. Conventional 2-shot uniform illumination TIE can give reliable inversion of the phase from intensity in many situations of practical interest; however, it has a null space consisting of fields with non-zero circulation of the Poynting vector. Here, we propose the hybrid illumination TIE method to disambiguate such objects. By comparing the diffraction signals using uniform and structured (sinusoidal) illumination patterns, we obtain a modulation-induced signal that depends solely on the phase gradient. In this way, we also increase signal sensitivity in the low spatial frequency region.

Phase discrepancy analysis and compensation for fast Fourier transform based solution of the transport of intensity equation

The transport of intensity equation (TIE) has long been recognized as a quantitative method for phase retrieval and phase contrast imaging. However, it is shown that the most widely accepted fast Fourier transform (FFT) based solutions do not provide an exact solution to the TIE in general. The root of the problem lies in the so-called "Teague's assumption" that the transverse flux is considered to be a conservative field, which cannot be satisfied for a general object. In this work, we present the theoretical analysis of the phase discrepancy owing to the Teague's assumption, and derive the necessary and sufficient conditions for the FFT-based solution to coincide with the exact phase. An iterative algorithm is then proposed aiming to compensate such phase discrepancy in a simple yet effective manner.

Reconstruction of wave front and object for inline holography from a set of detection planes

We illustrate the errors inherent in the conventional empty beam correction of full field X-ray propagation imaging, i.e. the division of intensities in the detection plane measured with an object in the beam by the intensity pattern measured without the object, i.e. the empty beam intensity pattern. The error of this conventional approximation is controlled by the ratio of the source size to the smallest feature in the object, as is shown by numerical simulation. In a second step, we investigate how to overcome the flawed empty beam division by simultaneous reconstruction of the probing wavefront (probe) and of the object, based on measurements in several detection planes (multi-projection approach). The algorithmic scheme is demonstrated numerically and experimentally, using the defocus wavefront of the hard X-ray nanoprobe setup at the European Synchrotron Radiation Facility (ESRF).

Argand-plane vorticity singularities in complex scalar optical fields: an experimental study using optical speckle

The Cornu spiral is, in essence, the image resulting from an Argand-plane map associated with monochromatic complex scalar plane waves diffracting from an infinite edge. Argand-plane maps can be useful in the analysis of more general optical fields. We experimentally study particular features of Argand-plane mappings known as "vorticity singularities" that are associated with mapping continuous single-valued complex scalar speckle fields to the Argand plane. Vorticity singularities possess a hierarchy of Argand-plane catastrophes including the fold, cusp and elliptic umbilic. We also confirm their connection to vortices in two-dimensional complex scalar waves. The study of vorticity singularities may also have implications for higher-dimensional fields such as coherence functions and multi-component fields such as vector and spinor fields.

Electron vortex production and control using aberration induced diffraction catastrophes

An aberration corrected electron microscope is used to create electron diffraction catastrophes, containing arrays of intensity zeros threading vortex cores. Vortices are ascribed to these arrays using catastrophe theory, scalar diffraction integrals, and experimentally retrieved phase maps. From measured wave function phases, obtained using focal-series phase retrieval, the orbital angular momentum density is mapped for highly astigmatic electron probes. We observe vortex rings and topological reconnections of nodal lines by tracking the vortex cores using the retrieved phases.

Phase measurement using an optical vortex lattice produced with a three-beam interferometer

A new phase-measurement technique is proposed, which utilizes a three-beam interferometer. Three-wave interference in the interferometer generates a uniform lattice of optical vortices, which is distorted by the presence of an object inserted in one arm of the interferometer. The transverse displacement of the vortices is proportional to the phase shift in the object wave. Tracking the vortices permits the phase of the object to be reconstructed. We demonstrate the method experimentally using a simple lens and a more complex object, namely the wing of a common house fly. Since the technique is implemented in real space, it is capable of reconstructing the phase locally.

Wavefield imaging via iterative retrieval based on phase modulation diversity

We present a fast and robust non-interferomentric wavefield retrieval approach suitable for imaging of both amplitude and phase distributions of scalar coherent beams. It is based on the diversity of the intensity measurements obtained under controlled astigmatism and it can be easily implemented in standard imaging systems. Its application for imaging in microscopy is experimentally studied. Relevant examples illustrate the approach capabilities for image super-resolution, numerical refocusing, quantitative imaging and phase mapping.

Transmission and emission x-ray microscopy: operation modes, contrast mechanisms and applications

Advances in microscopy techniques based on x-rays have opened unprecedented opportunities in terms of spatial resolution, combined with chemical and morphology sensitivity, to analyze solid, soft and liquid matter. The advent of ultrabright third and fourth generation photon sources and the continuous development of x-ray optics and detectors has pushed the limits of imaging and spectroscopic analysis to structures as small as a few tens of nanometers. Specific interactions of x-rays with matter provide elemental and chemical sensitivity that have made x-ray spectromicroscopy techniques a very attractive tool, complementary to other microscopies, for characterization in all actual research fields. The x-ray penetration power meets the demand to examine samples too thick for electron microscopes implementing 3D imaging and recently also 4D imaging which adds time resolution as well. Implementation of a variety of phase contrast techniques enhances the structural sensitivity, especially for the hard x-ray regime. Implementation of lensless or diffraction imaging helps to enhance the lateral resolution of x-ray imaging to the wavelength dependent diffraction limit.

Phase and amplitude imaging from noisy images by Kalman filtering

We propose and demonstrate a computational method for complex-field imaging from many noisy intensity images with varying defocus, using an extended complex Kalman filter. The technique offers dynamic smoothing of noisy measurements and is recursive rather than iterative, so is suitable for adaptive measurements. The Kalman filter provides near-optimal results in very low-light situations and may be adapted to propagation through turbulent, scattering, or nonlinear media.

Astigmatic phase retrieval: an experimental demonstration

We present the first experimental demonstration of the astigmatic phase retrieval technique, in which the diffracted wavefield is distorted by cylindrical curvature in two orthogonal directions. A charge-one vortex, a charge-two vortex, and a simple test image are all correctly reconstructed.

Astigmatic intensity equation for electron microscopy based phase retrieval

Phase retrieval, in principle, can be performed in a transmission electron microscope (TEM) using arbitrary aberrations of electron waves; provided that the aberrations are well-characterised and known. For example, the transport of intensity equation (TIE) can be used to infer the phase from a through-focus series of images. In this work an "astigmatic intensity equation" (AIE) is considered, which relates phase gradients to intensity variations caused by TEM objective lens focus and astigmatism variations. Within the paraxial approximation, it is shown that an exact solution of the AIE for the phase can be obtained using efficient Fourier transform methods. Experimental requirements for using the AIE are the measurement of a through-focus derivative and another intensity derivative, which is taken with respect to objective lens astigmatism variation. Two quasi-experimental investigations are conducted to test the validity of the solution.

Projected potential profiles across interfaces obtained by reconstructing the exit face wave function from through focal series

An iterative method for reconstructing the exit face wave function from a through focal series of transmission electron microscopy image line profiles across an interface is presented. Apart from high-resolution images recorded with small changes in defocus, this method works also well for a large defocus range as used for Fresnel imaging. Using the phase-object approximation the projected electrostatic as well as the absorptive potential profiles across an interface are determined from this exit face wave function. A new experimental image alignment procedure was developed in order to align images with large relative defocus shift. The performance of this procedure is shown to be superior to other image alignment procedures existing in the literature. The reconstruction method is applied to both simulated and experimental images.

Complex signal recovery from multiple fractional Fourier-transform intensities

The problem of recovering a complex signal from the magnitudes of any number of its fractional Fourier transforms at any set of fractional orders is addressed. This problem corresponds to the problem of phase retrieval from the transverse intensity profiles of an optical field at arbitrary locations in an optical system involving arbitrary concatenations of lenses and sections of free space. The dependence of the results on the number of orders, their spread, and the noise is investigated. Generally, increasing the number of orders improves the results, but with diminishing return beyond a certain point. Selecting the measurement planes such that their fractional orders are well separated or spread as much as possible also leads to better results.

Movable aperture lensless transmission microscopy: a novel phase retrieval algorithm

We propose an iterative phase retrieval method that uses a series of diffraction patterns, measured only in intensity, to solve for both amplitude and phase of the image wave function over a wide field of view and at wavelength-limited resolution. The new technique requires an aperture that is scanned to two or more positions over the object wave function. A simple implementation of the method is modeled and demonstrated, showing how the algorithm uses overlapping data in real space to resolve ambiguities in the solution. The technique opens up the possibility of practical transmission lensless microscopy at subatomic resolution using electrons, x rays, or nuclear particles.

Imaging a dense nanodot assembly by phase retrieval from TEM images

Phase retrieval is a classical inverse problem in many fields dealing with waves that is becoming of increasing interest in transmission electron microscopy (TEM). A non-interferometric approach is here applied to TEM images. Phase retrieval possibilities given by the transport intensity equation are compared to the ones deriving from the weak phase object approximation. In the limit of small angles, both methods lead to a similar equation between the phase and a set of defocus images. This equation can be solved by an image processing equivalent to using a specific filter in Fourier space. This processing leads to phase images with a spatial resolution here essentially limited by the defocus amount between images. A dense assembly of silicon nanodots is used as a model case to illustrate the interest of this approximate phase retrieval method which can be carried out on standard equipment. The dot heights estimated using the phase images are found to be in good agreement with ones measured by atomic force microscopy. Since image noise and large defocus values may strongly affect the solution given by the approximate method, an iterative phase retrieval method is also used as a test for working conditions.

Exit wave reconstruction at atomic resolution

An iterative method for exit wave function reconstruction based on wave function propagation in free space is presented. The method, which has the potential for application to many forms of microscopy, has been tailored to work with a through focal series of images measured in a high-resolution transmission electron microscope. Practical difficulties for exit wave reconstruction which are pertinent in this experimental environment are the slight incoherence of the electron beam, sample drift and its effect upon the defocus step size that can be utilised, and the number of image measurements that need to be made. To gauge the effectiveness of the method it is applied to experimental data that has been analysed previously using a maximum likelihood formalism (the MAL method).

Lensless imaging from diffraction intensity measurements by use of a noniterative phase-retrieval method

A method of reconstructing the complex amplitude of an object that is illuminated by a coherent wave from its Fresnel diffraction patterns is proposed for high-frequency wave phenomena such as x-rays and electron waves. A noniterative phase-retrieval method that uses a Gaussian filter is employed here, and it is shown that the object's illumination with amplitude distribution in the Fraunhofer diffraction pattern of a circular aperture can be used as a substitute for the Gaussian filter. This method has an advantage over other noniterative phase-retrieval methods in that it can retrieve phase vortices.

Resolution extension and exit wave reconstruction in complex HREM

Direct methods in real and reciprocal space are developed for structural reversion. The direct method in real space involves the use of a novel method to retrieve the phase in the image plane using transport of intensity equation/maximum entropy method (TIE/MEM) and exit wave reconstruction by self-consistent propagation. Since the exit wave is restored from the complex signal in the image planes, no image model between the exit wave and image is assumed. The structural information in the reconstructed exit wave is then further extended by a "complex" maximum entropy method as a direct method in reciprocal space to extrapolate the phase to higher frequencies.

Wave-function reconstruction of complex fields obeying nonlinear parabolic equations

We present a generalized Gerchberg-Saxton (GS) algorithm for reconstructing a (2+1)-dimensional complex scalar wave field which obeys a known nonlinear nondissipative parabolic differential equation, given knowledge of the wave-field modulus at three or more values of an evolution parameter such as time. This algorithm is used to recover the complex wave function of a (2+1)-dimensional Bose-Einstein condensate (BEC) from simulated modulus data. The Gross-Pitaveskii equation is used to model the dynamics of the BEC, with the modulus information being provided by a temporal sequence of simulated absorption images of the condensate. The efficacy of the generalized GS algorithm is examined for a wide range of simulation conditions, including strong nonlinearities, vortex states and Poisson noise. The general form of this algorithm, which allows one to reconstruct a time-dependent wave function, will be useful for studying the phase dynamics of topological defects in coherent quantum systems.

Phase recovery and lensless imaging by iterative methods in optical, X-ray and electron diffraction

Thomas Young's quantitative analysis of interference effects provided the confidence needed to revive the wave theory of light, and firmly established the concept of phase in optics. Phase plays a similarly fundamental role in matter-wave interferometry, for which the field-emission electron microscope provides ideal instrumentation. The wave-particle duality is vividly demonstrated by experimental 'Young's fringes' using coherent electron beams under conditions in which the flight time is less than the time between particle emission. A brief historical review is given of electron interferometry and holography, including the Aharonov-Bohm effect and the electron Sagnac interferometer. The simultaneous development of phase-contrast imaging at subnanometre spatial resolution has greatly deepened our understanding of atomic processes in biology, materials science and condensed-matter physics, while electron holography has become a routine tool for the mapping of electrostatic and magnetic fields in materials on a nanometre scale. The encoding of phase information in scattered farfield intensities is discussed, and non-interferometric, non-crystallographic methods for phase retrieval are reviewed in relationship to electron holography. Examples of phase measurement and diffraction-limited imaging using the hybrid input-output iterative algorithm are given, including simulations for soft X-ray imaging, and new experimental results for coherent electron and visible-light scattering. Image reconstruction is demonstrated from experimental electron and visible-light Fraunhofer diffraction patterns. The prospects this provides for lensless imaging using particles for which no lenses exist (such as neutrons, condensates, coherent atom beams and X-rays) are discussed. These new interactions can be expected to provide new information, perhaps, for example, in biology, with the advantage of less damage to samples.

Phase measurement of waves that obey nonlinear equations

We consider the problem of phase retrieval for classical and quantum wave fields that obey a wide class of nonlinear wave equations. This problem is tackled by means of a suitable generalization of existing methods based on the linear transport-of-intensity equation. The extension of these ideas to systems of coupled nonlinear wave equations is also given.

Observation of an x-ray vortex

Phase singularities are a ubiquitous feature of waves of all forms and represent a fundamental aspect of wave topology. An optical vortex phase singularity occurs when there is a spiral phase ramp about a point phase singularity. We report an experimental observation of an optical vortex in a field consisting of 9-keV x-ray photons. The vortex is created with an x-ray optical structure that imparts a spiral phase distribution to the incident wave field and is observed by use of diffraction about a wire to create a division-of-wave-front interferometer.

Computational aberration correction for an arbitrary linear imaging system

We show that aberration corrections can be made in any arbitrary linear imaging system provided the aberrations are well characterized and at least one of these aberrations can be independently varied in a well-controlled manner. We derive a generalization of the Schrödinger equation for wave propagation in aberration space assuming forward scattering. Transport equations in aberration space are derived. A general iterative algorithm which can retrieve the phase, and is robust in the presence of noise, is also derived. This is demonstrated using simulated data pertinent to electron microscopy, from a series of images with differing spherical aberration.