New algorithm enhances ptychographic image reconstruction
An international team of researchers that includes scientists from Berkeley Lab's Computational Research Division (CRD) and Center for Advanced Mathematics for Energy Research Applications (CAMERA) continues to find new ways to improve ptychographic image reconstruction.
In scientific experiments, ptychographic X-ray imaging is primarily used to characterize the structure and properties of matter and materials. While the method has been around for some 50 years, broad utilization has been hampered by the fact that the experimental process was slow and the computational processing of the data to produce a reconstructed image was expensive.
But in recent years advances in detectors and X-ray microscopes at light sources such as Berkeley Lab's Advanced Light Source (ALS) have made it possible to measure a ptychographic dataset in seconds. As a result, today ptychography is used in a range of scientific domains, including condensed-matter physics, cell biology and electronics.
In practice, X-ray ptychography works by focusing a beam of x-rays onto a spot in a sample. The scattering from the sample is recorded in the far field, and the recorded pattern is then phased to obtain the final image. The highest resolution attainable is not limited by the size of the focal spot, only by the numerical aperture and wavelength used. The phasing procedure in ptychography uses the overlap between consecutive exposures of the sample, plus the recorded far-field diffraction patterns, to reconstruct a high-resolution image of the sample.
As a result, reconstructing ptychographic datasets can be a data-intensive challenge that involves solving a difficult phase-retrieval problem, calibrating optical elements and dealing with experimental outliers and "noise." To address this challenge, Berkeley Lab scientists developed SHARP (scalable heterogeneous adaptive real-time ptychography), an algorithmic framework and computer software that enables the reconstruction of millions of phases of ptychographic image data per second. Since being introduce in 2016, SHARP has had a demonstrable impact on productivity for scientists working at the ALS and other light sources across the Department of Energy complex, with notable successes in the analysis of magnetic thin films, magnetozomes and 3-D battery materials.
Now researchers from CAMERA, the University of Texas and Tianjin Normal University—all members of the SHARP collaboration—have developed a model that further enhances SHARP's reconstruction capabilities. The new algorithm, GDP-ADMM (gradient decomposition of the probe/alternating direction method of multipliers), takes advantage of state-of-the-art mathematical aspects of phase retrieval, background noise optimization and detector "denoising" to enhance data acquisition and image resolution. With GDP-ADMM, SHARP is now able to handle more light than before, enabling faster acquisition and higher time resolution and ultimately more scientific discoveries.
A paper describing GDP-ADMM was the cover article in the May 2018 issue of Acta Crystallography Section A. GDP-ADMM allows more light to be used, opening the entrance slits of a ptychographic microscope and reducing the number of frames required to obtain enough data to reconstruct a meaningful image. The publication details how GDP-ADMM and partial coherence analysis help overcome stability issues inherent in coherent ptychographic imaging experiments, which often discard the majority of flux from a light source to define the coherence of an illumination (localized coherent X-ray probe). It also exploits translational kernel separability to speed analysis.
"The goal was to offer the ability to quickly discover interesting nanoparticles at full resolution by enabling rapid feedback from the microscopists at the beamlines," said Stefano Marchesini, a staff scientist in CRD and co-author on the Acta paper. "Even when the next-generation coherent light sources come online, we may be able to extend the X-ray energies that can be used in ptychography by using this model."