-
Radiation and Heat Transport in Divergent Shock-Bubble Interactions
Authors:
Kelin Kurzer-Ogul,
Brian M. Haines,
David S. Montgomery,
Silvia Pandolfi,
Joshua P. Sauppe,
Andrew F. T. Leong,
Daniel Hodge,
Pawel M. Kozlowski,
Stefano Marchesini,
Eric Cunningham,
Eric Galtier,
Dimitri Khaghani,
Hae Ja Lee,
Bob Nagler,
Richard L. Sandberg,
Arianna E. Gleason,
Hussein Aluie,
Jessica K. Shang
Abstract:
Shock-bubble interactions (SBI) are important across a wide range of physical systems. In inertial confinement fusion, interactions between laser-driven shocks and micro-voids in both ablators and foam targets generate instabilities that are a major obstacle in achieving ignition. Experiments imaging the collapse of such voids at high energy densities (HED) are constrained by spatial and temporal…
▽ More
Shock-bubble interactions (SBI) are important across a wide range of physical systems. In inertial confinement fusion, interactions between laser-driven shocks and micro-voids in both ablators and foam targets generate instabilities that are a major obstacle in achieving ignition. Experiments imaging the collapse of such voids at high energy densities (HED) are constrained by spatial and temporal resolution, making simulations a vital tool in understanding these systems. In this study, we benchmark several radiation and thermal transport models in the xRAGE hydrodynamic code against experimental images of a collapsing mesoscale void during the passage of a 300 GPa shock. We also quantitatively examine the role of transport physics in the evolution of the SBI. This allows us to understand the dynamics of the interaction at timescales shorter than experimental imaging framerates. We find that all radiation models examined reproduce empirical shock velocities within experimental error. Radiation transport is found to reduce shock pressures by providing an additional energy pathway in the ablation region, but this effect is small ($\sim$1\% of total shock pressure). Employing a flux-limited Spitzer model for heat conduction, we find that flux limiters between 0.03 and 0.10 produce agreement with experimental velocities, suggesting that the system is well-within the Spitzer regime. Higher heat conduction is found to lower temperatures in the ablated plasma and to prevent secondary shocks at the ablation front, resulting in weaker primary shocks. Finally, we confirm that the SBI-driven instabilities observed in the HED regime are baroclinically driven, as in the low energy case.
△ Less
Submitted 5 March, 2024;
originally announced March 2024.
-
Fast Iterative Algorithms for Blind Phase Retrieval: A survey
Authors:
Huibin Chang,
Li Yang,
Stefano Marchesini
Abstract:
In nanoscale imaging technique and ultrafast laser, the reconstruction procedure is normally formulated as a blind phase retrieval (BPR) problem, where one has to recover both the sample and the probe (pupil) jointly from phaseless data. This survey first presents the mathematical formula of BPR, related nonlinear optimization problems and then gives a brief review of the recent iterative algorith…
▽ More
In nanoscale imaging technique and ultrafast laser, the reconstruction procedure is normally formulated as a blind phase retrieval (BPR) problem, where one has to recover both the sample and the probe (pupil) jointly from phaseless data. This survey first presents the mathematical formula of BPR, related nonlinear optimization problems and then gives a brief review of the recent iterative algorithms. It mainly consists of three types of algorithms, including the operator-splitting based first-order optimization methods, second order algorithm with Hessian,and subspace methods. The future research directions for experimental issues and theoretical analysis are further discussed.
△ Less
Submitted 12 November, 2022;
originally announced November 2022.
-
In-situ observation of the formation of laser-induced periodic surface structures with extreme spatial and temporal resolution
Authors:
K. Sokolowski-Tinten,
J. Bonse,
A. Barty,
H. N. Chapman,
S. Bajt,
M. J. Bogan,
S. Boutet,
A. Cavalleri,
S. Düsterer,
M. Frank,
J. Hajdu,
S. Hau-Riege,
S. Marchesini,
N. Stojanovic,
R. Treusch
Abstract:
Irradiation of solid surfaces with intense ultrashort laser pulses represents a unique way of depositing energy into materials. It allows to realize states of extreme electronic excitation and/or very high temperature and pressure, and to drive materials close to and beyond fundamental stability limits. As a consequence, structural changes and phase transitions often occur along unusual pathways a…
▽ More
Irradiation of solid surfaces with intense ultrashort laser pulses represents a unique way of depositing energy into materials. It allows to realize states of extreme electronic excitation and/or very high temperature and pressure, and to drive materials close to and beyond fundamental stability limits. As a consequence, structural changes and phase transitions often occur along unusual pathways and under strongly non-equilibrium conditions. Due to the inherent multiscale nature - both temporally and spatially - of these irreversible processes their direct experimental observation requires techniques that combine high temporal resolution with the appropriate spatial resolution and the capability to obtain good quality data on a single pulse/event basis. In this respect fourth generation light sources, namely short wavelength, short pulse free electron lasers (FELs) are offering new and fascinating possibilities. As an example, this chapter will discuss the results of scattering experiments carried at the FLASH free electron laser at DESY (Hamburg, Germany), which allowed us to resolve laser-induced structure formation at surfaces on the nanometer to sub-micron length scale and in temporal regimes ranging from picoseconds to several nanoseconds with sub-picosecond resolution.
△ Less
Submitted 9 June, 2022;
originally announced June 2022.
-
Overlapping Domain Decomposition Methods for Ptychographic Imaging
Authors:
Huibin Chang,
Roland Glowinski,
Stefano Marchesini,
Xue-cheng Tai,
Yang Wang,
Tieyong Zeng
Abstract:
In ptychography experiments, redundant scanning is usually required to guarantee the stable recovery, such that a huge amount of frames are generated, and thus it poses a great demand of parallel computing in order to solve this large-scale inverse problem. In this paper, we propose the overlapping Domain Decomposition Methods(DDMs) to solve the nonconvex optimization problem in ptychographic imag…
▽ More
In ptychography experiments, redundant scanning is usually required to guarantee the stable recovery, such that a huge amount of frames are generated, and thus it poses a great demand of parallel computing in order to solve this large-scale inverse problem. In this paper, we propose the overlapping Domain Decomposition Methods(DDMs) to solve the nonconvex optimization problem in ptychographic imaging. They decouple the problem defined on the whole domain into subproblems only defined on the subdomains with synchronizing information in the overlapping regions of these subdomains,thus leading to highly parallel algorithms with good load balance. More specifically, for the nonblind recovery (with known probe in advance), by enforcing the continuity of the overlapping regions for the image (sample), the nonlinear optimization model is established based on a novel smooth-truncated amplitude-Gaussian metric (ST-AGM). Such metric allows for fast calculation of the proximal mapping with closed form, and meanwhile provides the possibility for the convergence guarantee of the first-order nonconvex optimization algorithm due to its Lipschitz smoothness. Then the Alternating Direction Method of Multipliers (ADMM) is utilized to generate an efficient Overlapping Domain Decomposition based Ptychography algorithm(OD2P) for the two-subdomain domain decomposition (DD), where all subproblems can be computed with close-form solutions.Due to the Lipschitz continuity for the gradient of the objective function with ST-AGM, the convergence of the proposed OD2P is derived under mild conditions. Moreover, it is extended to more general case including multiple-subdomain DD and blind recovery. Numerical experiments are further conducted to show the performance of proposed algorithms, demonstrating good convergence speed and robustness to the noise.
△ Less
Submitted 4 February, 2021; v1 submitted 30 October, 2020;
originally announced November 2020.
-
Sparse Matrix-Based HPC Tomography
Authors:
Stefano Marchesini,
Anuradha Trivedi,
Pablo Enfedaque,
Talita Perciano,
Dilworth Parkinson
Abstract:
Tomographic imaging has benefited from advances in X-ray sources, detectors and optics to enable novel observations in science, engineering and medicine. These advances have come with a dramatic increase of input data in the form of faster frame rates, larger fields of view or higher resolution, so high performance solutions are currently widely used for analysis. Tomographic instruments can vary…
▽ More
Tomographic imaging has benefited from advances in X-ray sources, detectors and optics to enable novel observations in science, engineering and medicine. These advances have come with a dramatic increase of input data in the form of faster frame rates, larger fields of view or higher resolution, so high performance solutions are currently widely used for analysis. Tomographic instruments can vary significantly from one to another, including the hardware employed for reconstruction: from single CPU workstations to large scale hybrid CPU/GPU supercomputers. Flexibility on the software interfaces and reconstruction engines are also highly valued to allow for easy development and prototyping. This paper presents a novel software framework for tomographic analysis that tackles all aforementioned requirements. The proposed solution capitalizes on the increased performance of sparse matrix-vector multiplication and exploits multi-CPU and GPU reconstruction over MPI. The solution is implemented in Python and relies on CuPy for fast GPU operators and CUDA kernel integration, and on SciPy for CPU sparse matrix computation. As opposed to previous tomography solutions that are tailor-made for specific use cases or hardware, the proposed software is designed to provide flexible, portable and high-performance operators that can be used for continuous integration at different production environments, but also for prototyping new experimental settings or for algorithmic development. The experimental results demonstrate how our implementation can even outperform state-of-the-art software packages used at advanced X-ray sources worldwide.
△ Less
Submitted 22 April, 2020; v1 submitted 27 March, 2020;
originally announced March 2020.
-
Analyzer Free Linear Dichroic Ptychography
Authors:
Huibin Chang,
Matthew A Marcus,
Stefano Marchesini
Abstract:
Linear-dichroism is an important tool to characterize the transmission matrix and determine the crystal or orbital orientation in a material. In order to gain high resolution mapping of the transmission properties of such materials, we introduce the linear-dichroism scattering model in ptychographic imaging, and then develop an efficient two-stage reconstruction algorithm. Using proposed algorithm…
▽ More
Linear-dichroism is an important tool to characterize the transmission matrix and determine the crystal or orbital orientation in a material. In order to gain high resolution mapping of the transmission properties of such materials, we introduce the linear-dichroism scattering model in ptychographic imaging, and then develop an efficient two-stage reconstruction algorithm. Using proposed algorithm, the dichroic transmission matrix without an analyzer can be recovered by using ptychography measurements with as few as three different polarization angles, with the help of an empty region to remove phase ambiguities.
△ Less
Submitted 10 February, 2020;
originally announced February 2020.
-
Iterative X-ray Spectroscopic Ptychography
Authors:
Huibin Chang,
Ziqin Rong,
Pablo Enfedaque,
Stefano Marchesini
Abstract:
Spectroscopic ptychography is a powerful technique to determine the chemical composition of a sample with high spatial resolution. In spectro-ptychography, a sample is rastered through a focused x-ray beam with varying photon energy so that a series of phaseless diffraction data are recorded. Each chemical component in the material under investigation has a characteristic absorption and phase cont…
▽ More
Spectroscopic ptychography is a powerful technique to determine the chemical composition of a sample with high spatial resolution. In spectro-ptychography, a sample is rastered through a focused x-ray beam with varying photon energy so that a series of phaseless diffraction data are recorded. Each chemical component in the material under investigation has a characteristic absorption and phase contrast as a function of photon energy. Using a dictionary formed by the set of contrast functions of each energy for each chemical component, it is possible to obtain the chemical composition of the material from high resolution multi-spectral images. This paper presents SPA (Spectroscopic Ptychography with ADMM), a novel algorithm to iteratively solve the spectroscopic blind ptychography problem. We design first a nonlinear spectro-ptychography model based on Poisson maximum likelihood, and construct then the proposed method based on fast iterative splitting operators. SPA can be used to retrieve spectral contrast when considering both a known or an incomplete (partially known) dictionary of reference spectra. By coupling the redundancy across different spectral measurements, the proposed algorithm can achieve higher reconstruction quality when compared to standard state-of-the-art two-step methods. We demonstrate how SPA can recover accurate chemical maps from Poisson-noised measurements, and also show its enhanced robustness when reconstructing reduced redundancy ptychography data using large scanning stepsizes.
△ Less
Submitted 15 July, 2020; v1 submitted 31 May, 2019;
originally announced June 2019.
-
One-Dimensional Phase Retrieval: Regularization, Box Relaxation and Uniqueness
Authors:
Wing Hong Wong,
Yifei Lou,
Stefano Marchesini,
Tieyong Zeng
Abstract:
Recovering a signal from its Fourier magnitude is referred to as phase retrieval, which occurs in different fields of engineering and applied physics. This paper gives a new characterization of the phase retrieval problem. Particularly useful is the analysis revealing that the common gradient-based regularization does not restrict the set of solutions to a smaller set. Specifically focusing on bin…
▽ More
Recovering a signal from its Fourier magnitude is referred to as phase retrieval, which occurs in different fields of engineering and applied physics. This paper gives a new characterization of the phase retrieval problem. Particularly useful is the analysis revealing that the common gradient-based regularization does not restrict the set of solutions to a smaller set. Specifically focusing on binary signals, we show that a box relaxation is equivalent to the binary constraint for Fourier-types of phase retrieval. We further prove that binary signals can be recovered uniquely up to trivial ambiguities under certain conditions. Finally, we use the characterization theorem to develop an efficient denoising algorithm.
△ Less
Submitted 15 July, 2020; v1 submitted 23 April, 2019;
originally announced April 2019.
-
Iterative Joint Ptychography-Tomography with Total Variation Regularization
Authors:
Huibin Chang,
Pablo Enfedaque,
Stefano Marchesini
Abstract:
In order to determine the 3D structure of a thick sample, researchers have recently combined ptychography (for high resolution) and tomography (for 3D imaging) in a single experiment. 2-step methods are usually adopted for reconstruction, where the ptychography and tomography problems are often solved independently. In this paper, we provide a novel model and ADMM-based algorithm to jointly solve…
▽ More
In order to determine the 3D structure of a thick sample, researchers have recently combined ptychography (for high resolution) and tomography (for 3D imaging) in a single experiment. 2-step methods are usually adopted for reconstruction, where the ptychography and tomography problems are often solved independently. In this paper, we provide a novel model and ADMM-based algorithm to jointly solve the ptychography-tomography problem iteratively, also employing total variation regularization. The proposed method permits large scan stepsizes for the ptychography experiment, requiring less measurements and being more robust to noise with respect to other strategies, while achieving higher reconstruction quality results.
△ Less
Submitted 14 February, 2019;
originally announced February 2019.
-
Shaping Coherent X-rays with Binary Optics
Authors:
Stefano Marchesini,
Anne Sakdinawat
Abstract:
Diffractive lenses fabricated by lithographic methods are one of the most popular image forming optics in the x-ray regime. Most commonly, binary diffractive optics, such as Fresnel zone plates are used due to their ability to focus at high resolution and to manipulate the x-ray wavefront. We report here a binary zone plate design strategy to form arbitrary illuminations for coherent multiplexing,…
▽ More
Diffractive lenses fabricated by lithographic methods are one of the most popular image forming optics in the x-ray regime. Most commonly, binary diffractive optics, such as Fresnel zone plates are used due to their ability to focus at high resolution and to manipulate the x-ray wavefront. We report here a binary zone plate design strategy to form arbitrary illuminations for coherent multiplexing, structured illumination, and wavefront shaping experiments. Given a desired illumination, we adjust the duty cycle, harmonic order, and zone placement to vary both the amplitude and phase of the wavefront at the lens. This enables the binary lithographic pattern to generate arbitrary structured illumination optimized for a variety of applications such as holography, interferometry, ptychography, imaging, and others.
△ Less
Submitted 25 December, 2018;
originally announced December 2018.
-
Advanced Denoising for X-ray Ptychography
Authors:
Huibin Chang,
Pablo Enfedaque,
Jie Zhang,
Juliane Reinhardt,
Bjoern Enders,
Young-Sang Yu,
David Shapiro,
Christian G. Schroer,
Tieyong Zeng,
Stefano Marchesini
Abstract:
The success of ptychographic imaging experiments strongly depends on achieving high signal-to-noise ratio. This is particularly important in nanoscale imaging experiments when diffraction signals are very weak and the experiments are accompanied by significant parasitic scattering (background), outliers or correlated noise sources. It is also critical when rare events such as cosmic rays, or bad f…
▽ More
The success of ptychographic imaging experiments strongly depends on achieving high signal-to-noise ratio. This is particularly important in nanoscale imaging experiments when diffraction signals are very weak and the experiments are accompanied by significant parasitic scattering (background), outliers or correlated noise sources. It is also critical when rare events such as cosmic rays, or bad frames caused by electronic glitches or shutter timing malfunction take place.
In this paper, we propose a novel iterative algorithm with rigorous analysis that exploits the direct forward model for parasitic noise and sample smoothness to achieve a thorough characterization and removal of structured and random noise. We present a formal description of the proposed algorithm and prove its convergence under mild conditions. Numerical experiments from simulations and real data (both soft and hard X-ray beamlines) demonstrate that the proposed algorithms produce better results when compared to state-of-the-art methods.
△ Less
Submitted 14 February, 2019; v1 submitted 5 November, 2018;
originally announced November 2018.
-
Blind Ptychographic Phase Retrieval via Convergent Alternating Direction Method of Multipliers
Authors:
Huibin Chang,
Pable Enfedaque,
Stefano Marchesini
Abstract:
Ptychography has risen as a reference X-ray imaging technique: it achieves resolutions of one billionth of a meter, macroscopic field of view, or the capability to retrieve chemical or magnetic contrast, among other features. A ptychographyic reconstruction is normally formulated as a blind phase retrieval problem, where both the image (sample) and the probe (illumination) have to be recovered fro…
▽ More
Ptychography has risen as a reference X-ray imaging technique: it achieves resolutions of one billionth of a meter, macroscopic field of view, or the capability to retrieve chemical or magnetic contrast, among other features. A ptychographyic reconstruction is normally formulated as a blind phase retrieval problem, where both the image (sample) and the probe (illumination) have to be recovered from phaseless measured data. In this article we address a nonlinear least squares model for the blind ptychography problem with constraints on the image and the probe by maximum likelihood estimation of the Poisson noise model. We formulate a variant model that incorporates the information of phaseless measurements of the probe to eliminate possible artifacts. Next, we propose a generalized alternating direction method of multipliers designed for the proposed nonconvex models with convergence guarantee under mild conditions, where their subproblems can be solved by fast element-wise operations. Numerically, the proposed algorithm outperforms state-of-the-art algorithms in both speed and image quality.
△ Less
Submitted 24 September, 2018; v1 submitted 17 August, 2018;
originally announced August 2018.
-
Three dimensional localization of nanoscale battery reactions using soft X-ray tomography
Authors:
Young-Sang Yu,
Maryam Farmand,
Chunjoong Kim,
Yijin Liu,
Clare P. Grey,
Fiona C. Strobridge,
Tolek Tyliszczak,
Rich Celestre,
Peter Denes,
John Joseph,
Harinarayan Krishnan,
Filipe R. N. C. Maia,
A. L. David Kilcoyne,
Stefano Marchesini,
Talita Perciano Costa Leite,
Tony Warwick,
Howard Padmore,
Jordi Cabana,
David A. Shapiro
Abstract:
Battery function is determined by the efficiency and reversibility of the electrochemical phase transformations at solid electrodes. The microscopic tools available to study the chemical states of matter with the required spatial resolution and chemical specificity are intrinsically limited when studying complex architectures by their reliance on two dimensional projections of thick material. Here…
▽ More
Battery function is determined by the efficiency and reversibility of the electrochemical phase transformations at solid electrodes. The microscopic tools available to study the chemical states of matter with the required spatial resolution and chemical specificity are intrinsically limited when studying complex architectures by their reliance on two dimensional projections of thick material. Here, we report the development of soft X-ray ptychographic tomography, which resolves chemical states in three dimensions at 11-nm spatial resolution. We study an ensemble of nano-plates of lithium iron phosphate (LixFePO4) extracted from a battery electrode at 50% state of charge. Using a set of nanoscale tomograms, we quantify the electrochemical state and resolve phase boundaries throughout the volume of individual nano-particles. These observations reveal multiple reaction points and intra-particle heterogeneity that highlights the importance of electrical connectivity, providing novel insight to the design of the next generation of high-performance devices.
△ Less
Submitted 9 March, 2018; v1 submitted 4 November, 2017;
originally announced November 2017.
-
Partially Coherent Ptychography by Gradient Decomposition of the Probe
Authors:
Huibin Chang,
Pablo Enfedaque,
Yifei Lou,
Stefano Marchesini
Abstract:
Coherent ptychographic imaging experiments often discard over 99.9 % of the flux from a light source to define the coherence of an illumination. Even when coherent flux is sufficient, the stability required during an exposure is another important limiting factor. Partial coherence analysis can considerably reduce these limitations. A partially coherent illumination can often be written as the supe…
▽ More
Coherent ptychographic imaging experiments often discard over 99.9 % of the flux from a light source to define the coherence of an illumination. Even when coherent flux is sufficient, the stability required during an exposure is another important limiting factor. Partial coherence analysis can considerably reduce these limitations. A partially coherent illumination can often be written as the superposition of a single coherent illumination convolved with a separable translational kernel. In this paper we propose the Gradient Decomposition of the Probe (GDP), a model that exploits translational kernel separability, coupling the variances of the kernel with the transverse coherence. We describe an efficient first-order splitting algorithm GDP-ADMM to solve the proposed nonlinear optimization problem. Numerical experiments demonstrate the effectiveness of the proposed method with Gaussian and binary kernel functions in fly-scan measurements. Remarkably, GDP-ADMM produces satisfactory results even when the ratio between kernel width and beam size is more than one, or when the distance between successive acquisitions is twice as large as the beam width.
△ Less
Submitted 19 October, 2017;
originally announced October 2017.
-
Denoising Poisson Phaseless Measurements via Orthogonal Dictionary Learning
Authors:
Huibin Chang,
Stefano Marchesini
Abstract:
Phaseless diffraction measurements recorded by a CCD detector are often affected by Poisson noise. In this paper, we propose a dictionary learning model by employing patches based sparsity to denoise Poisson phaseless measurement. The model consists of three terms: (i) A representation term by an orthogonal dictionary, (ii) an $L^0$ pseudo norm of coefficient matrix, and (iii) a Kullback-Leibler d…
▽ More
Phaseless diffraction measurements recorded by a CCD detector are often affected by Poisson noise. In this paper, we propose a dictionary learning model by employing patches based sparsity to denoise Poisson phaseless measurement. The model consists of three terms: (i) A representation term by an orthogonal dictionary, (ii) an $L^0$ pseudo norm of coefficient matrix, and (iii) a Kullback-Leibler divergence to fit phaseless Poisson data. Fast Alternating Minimization Method (AMM) and Proximal Alternating Linearized Minimization (PALM) are adopted to solve the established model with convergence guarantee, and especially global convergence for PALM is derived. The subproblems for two algorithms have fast solvers, and indeed, the solutions for the sparse coding and dictionary updating both have closed forms due to the orthogonality of learned dictionaries. Numerical experiments for phase retrieval using coded diffraction and ptychographic patterns are performed to show the efficiency and robustness of proposed methods, which, by preserving texture features, produce visually and quantitatively improved denoised images compared with other phase retrieval algorithms without regularization and local sparsity promoting algorithms.
△ Less
Submitted 25 December, 2016;
originally announced December 2016.
-
A general framework for denoising phaseless diffraction measurements
Authors:
Huibin Chang,
Stefano Marchesini
Abstract:
We propose a general framework to recover underlying images from noisy phaseless diffraction measurements based on the alternating directional method of multipliers and the plug-and-play technique. The algorithm consists of three-step iterations: (i) Solving a generalized least square problem with the maximum a posteriori (MAP) estimate of the noise, (ii) Gaussian denoising and (iii) updating the…
▽ More
We propose a general framework to recover underlying images from noisy phaseless diffraction measurements based on the alternating directional method of multipliers and the plug-and-play technique. The algorithm consists of three-step iterations: (i) Solving a generalized least square problem with the maximum a posteriori (MAP) estimate of the noise, (ii) Gaussian denoising and (iii) updating the multipliers. The denoising step utilizes higher order filters such as total generalized variation and nonlocal sparsity based filters including nonlocal mean (NLM) and Block-matching and 3D filtering (BM3D) filters. The multipliers are updated by a symmetric technique to increase convergence speed. The proposed method with low computational complexity is provided with theoretical convergence guarantee, and it enables recovering images with sharp edges, clean background and repetitive features from noisy phaseless measurements. Numerous numerical experiments for Fourier phase retrieval (PR) as coded diffraction and ptychographic patterns are performed to verify the convergence and efficiency, showing that our proposed method outperforms the state-of-art PR algorithms without any regularization and those with total variational regularization.
△ Less
Submitted 2 November, 2016;
originally announced November 2016.
-
Nanosurveyor: a framework for real-time data processing
Authors:
Benedikt J. Daurer,
Hari Krishnan,
Talita Perciano,
Filipe R. N. C. Maia,
David A. Shapiro,
James A. Sethian,
Stefano Marchesini
Abstract:
Scientists are drawn to synchrotrons and accelerator based light sources because of their brightness, coherence and flux. The rate of improvement in brightness and detector technology has outpaced Moore's law growth seen for computers, networks, and storage, and is enabling novel observations and discoveries with faster frame rates, larger fields of view, higher resolution, and higher dimensionali…
▽ More
Scientists are drawn to synchrotrons and accelerator based light sources because of their brightness, coherence and flux. The rate of improvement in brightness and detector technology has outpaced Moore's law growth seen for computers, networks, and storage, and is enabling novel observations and discoveries with faster frame rates, larger fields of view, higher resolution, and higher dimensionality. Here we present an integrated software/algorithmic framework designed to capitalize on high throughput experiments, and describe the streamlined processing pipeline of ptychography data analysis. The pipeline provides throughput, compression, and resolution as well as rapid feedback to the microscope operators.
△ Less
Submitted 9 September, 2016;
originally announced September 2016.
-
SHARP: a distributed, GPU-based ptychographic solver
Authors:
Stefano Marchesini,
Hari Krishnan,
Benedikt J. Daurer,
David A. Shapiro,
Talita Perciano,
James A. Sethian,
Filipe R. N. C. Maia
Abstract:
Ever brighter light sources, fast parallel detectors, and advances in phase retrieval methods, have made ptychography a practical and popular imaging technique. Compared to previous techniques, ptychography provides superior robustness and resolution at the expense of more advanced and time consuming data analysis. By taking advantage of massively parallel architectures, high-throughput processing…
▽ More
Ever brighter light sources, fast parallel detectors, and advances in phase retrieval methods, have made ptychography a practical and popular imaging technique. Compared to previous techniques, ptychography provides superior robustness and resolution at the expense of more advanced and time consuming data analysis. By taking advantage of massively parallel architectures, high-throughput processing can expedite this analysis and provide microscopists with immediate feedback. These advances allow real-time imaging at wavelength limited resolution, coupled with a large field of view. Here, we introduce a set of algorithmic and computational methodologies used at the Advanced Light Source, and DOE light sources packaged as a CUDA based software environment named SHARP (http://camera.lbl.gov/sharp), aimed at providing state-of-the-art high-throughput ptychography reconstructions for the coming era of diffraction limited light sources.
△ Less
Submitted 20 June, 2016; v1 submitted 29 January, 2016;
originally announced February 2016.
-
Rank-1 accelerated illumination recovery in scanning diffractive imaging by transparency estimation
Authors:
Stefano Marchesini,
Hau-Tieng Wu
Abstract:
We consider the problem of blind ptychography, that is the joint estimation of an unknown object and an illumination function from diffraction intensity measurements. In ptychography, diffraction measurements from neighboring regions of the same object are related to each other by a pairwise relationship between overlapping frames. When the illumination is well known, the relationship among frames…
▽ More
We consider the problem of blind ptychography, that is the joint estimation of an unknown object and an illumination function from diffraction intensity measurements. In ptychography, diffraction measurements from neighboring regions of the same object are related to each other by a pairwise relationship between overlapping frames. When the illumination is well known, the relationship among frames is given by a linear projection operator. We propose a power iteration-projection algorithm that minimizes the global pairwise discrepancy among frames. We accelerate the convergence of power method by subtracting the estimated localized average transparency of the unknown object. The method is effective for weakly scattering and low contrast objects or piecewise smooth specimens.
△ Less
Submitted 2 February, 2016; v1 submitted 8 August, 2014;
originally announced August 2014.
-
Alternating Projection, Ptychographic Imaging and Phase Synchronization
Authors:
Stefano Marchesini,
Yu-Chao Tu,
Hau-tieng Wu
Abstract:
We demonstrate necessary and sufficient conditions of the local convergence of the alternating projection algorithm to a unique solution up to a global phase factor. Additionally, for the ptychography imaging problem, we discuss phase synchronization and graph connection Laplacian, and show how to construct an accurate initial guess to accelerate convergence speed to handle the big imaging data in…
▽ More
We demonstrate necessary and sufficient conditions of the local convergence of the alternating projection algorithm to a unique solution up to a global phase factor. Additionally, for the ptychography imaging problem, we discuss phase synchronization and graph connection Laplacian, and show how to construct an accurate initial guess to accelerate convergence speed to handle the big imaging data in the coming new light source era.
△ Less
Submitted 19 July, 2015; v1 submitted 3 February, 2014;
originally announced February 2014.
-
Phase retrieval by power iterations
Authors:
Stefano Marchesini
Abstract:
I show that the power iteration method applied to the phase retrieval problem converges under special conditions. One is given the relative phases between small non-overlapping groups of pixels of a recorded intensity pattern, but no information on the phase between the groups of pixels. Numerical tests show that the inverse block iteration recovers the solution in 1 iteration.
I show that the power iteration method applied to the phase retrieval problem converges under special conditions. One is given the relative phases between small non-overlapping groups of pixels of a recorded intensity pattern, but no information on the phase between the groups of pixels. Numerical tests show that the inverse block iteration recovers the solution in 1 iteration.
△ Less
Submitted 3 December, 2012;
originally announced December 2012.
-
Augmented projections for ptychographic imaging
Authors:
Stefano Marchesini,
Andre Schirotzek,
Chao Yang,
Hau-tieng Wu,
Filipe Maia
Abstract:
Ptychography is a popular technique to achieve diffraction limited resolution images of a two or three dimensional sample using high frame rate detectors. We introduce a relaxation of common projection algorithms to account for instabilities given by intensity and background fluctuations, position errors, or poor calibration using multiplexing illumination. This relaxation introduces an additional…
▽ More
Ptychography is a popular technique to achieve diffraction limited resolution images of a two or three dimensional sample using high frame rate detectors. We introduce a relaxation of common projection algorithms to account for instabilities given by intensity and background fluctuations, position errors, or poor calibration using multiplexing illumination. This relaxation introduces an additional phasing optimization at every step that enhances the convergence rate of common projection algorithms. Numerical tests exhibit the exact recovery of the object and the noise when there is high redundancy in the data.
△ Less
Submitted 29 August, 2013; v1 submitted 21 September, 2012;
originally announced September 2012.
-
Iterative Algorithms for Ptychographic Phase Retrieval
Authors:
C. Yang,
J. Qian,
A. Schirotzek,
F. Maia,
S. Marchesini
Abstract:
Ptychography promises diffraction limited resolution without the need for high resolution lenses. To achieve high resolution one has to solve the phase problem for many partially overlapping frames. Here we review some of the existing methods for solving ptychographic phase retrieval problem from a numerical analysis point of view, and propose alternative methods based on numerical optimization.
Ptychography promises diffraction limited resolution without the need for high resolution lenses. To achieve high resolution one has to solve the phase problem for many partially overlapping frames. Here we review some of the existing methods for solving ptychographic phase retrieval problem from a numerical analysis point of view, and propose alternative methods based on numerical optimization.
△ Less
Submitted 27 May, 2011;
originally announced May 2011.
-
Compressive auto-indexing in femtosecond nanocrystallography
Authors:
F. R. N. C. Maia,
C. Yang,
S. Marchesini
Abstract:
Ultrafast nanocrystallography has the potential to revolutionize biology by enabling structural elucidation of proteins for which it is possible to grow crystals with 10 or fewer unit cells on the side. The success of nanocrystallography depends on robust orientation-determination procedures that allow us to average diffraction data from multiple nanocrystals to produce a three dimensional (3D) di…
▽ More
Ultrafast nanocrystallography has the potential to revolutionize biology by enabling structural elucidation of proteins for which it is possible to grow crystals with 10 or fewer unit cells on the side. The success of nanocrystallography depends on robust orientation-determination procedures that allow us to average diffraction data from multiple nanocrystals to produce a three dimensional (3D) diffraction data volume with a high signal-to-noise ratio. Such a 3D diffraction volume can then be phased using standard crystallographic techniques. "Indexing" algorithms used in crystallography enable orientation determination of diffraction data from a single crystal when a relatively large number of reflections are recorded. Here we show that it is possible to obtain the exact lattice geometry from a smaller number of measurements than standard approaches using a basis pursuit solver.
△ Less
Submitted 12 November, 2010;
originally announced November 2010.
-
Compressive Phase Contrast Tomography
Authors:
F. R. N. C. Maia,
A. MacDowell,
S. Marchesini,
H. A. Padmore,
D. Y. Parkinson,
J. Pien,
A. Schirotzek,
C. Yang
Abstract:
When x-rays penetrate soft matter, their phase changes more rapidly than their amplitude. In- terference effects visible with high brightness sources creates higher contrast, edge enhanced images. When the object is piecewise smooth (made of big blocks of a few components), such higher con- trast datasets have a sparse solution. We apply basis pursuit solvers to improve SNR, remove ring artifacts,…
▽ More
When x-rays penetrate soft matter, their phase changes more rapidly than their amplitude. In- terference effects visible with high brightness sources creates higher contrast, edge enhanced images. When the object is piecewise smooth (made of big blocks of a few components), such higher con- trast datasets have a sparse solution. We apply basis pursuit solvers to improve SNR, remove ring artifacts, reduce the number of views and radiation dose from phase contrast datasets collected at the Hard X-Ray Micro Tomography Beamline at the Advanced Light Source. We report a GPU code for the most computationally intensive task, the gridding and inverse gridding algorithm (non uniform sampled Fourier transform).
△ Less
Submitted 7 September, 2010;
originally announced September 2010.
-
Cryptotomography: reconstructing 3D Fourier intensities from randomly oriented single-shot diffraction patterns
Authors:
N. D. Loh,
M. Bogan,
V. Elser,
A. Barty,
S. Boutet,
S. Bajt,
J. Hajdu,
T. Ekeberg,
F. R. N. C. Maia,
J. Schulz,
M. M. Seibert,
B. Iwan,
N. Timneanu,
S. Marchesini,
I. Schlichting,
R. L. Shoeman,
L. Lomb,
M. Frank,
M. Liang,
H. N. Chapman
Abstract:
We reconstructed the 3D Fourier intensity distribution of mono-disperse prolate nano-particles using single-shot 2D coherent diffraction patterns collected at DESY's FLASH facility when a bright, coherent, ultrafast X-ray pulse intercepted individual particles of random, unmeasured orientations. This first experimental demonstration of cryptotomography extended the Expansion-Maximization-Compres…
▽ More
We reconstructed the 3D Fourier intensity distribution of mono-disperse prolate nano-particles using single-shot 2D coherent diffraction patterns collected at DESY's FLASH facility when a bright, coherent, ultrafast X-ray pulse intercepted individual particles of random, unmeasured orientations. This first experimental demonstration of cryptotomography extended the Expansion-Maximization-Compression (EMC) framework to accommodate unmeasured fluctuations in photon fluence and loss of data due to saturation or background scatter. This work is an important step towards realizing single-shot diffraction imaging of single biomolecules.
△ Less
Submitted 5 April, 2010; v1 submitted 3 March, 2010;
originally announced March 2010.
-
Ab initio compressive phase retrieval
Authors:
S. Marchesini
Abstract:
Any object on earth has two fundamental properties: it is finite, and it is made of atoms. Structural information about an object can be obtained from diffraction amplitude measurements that account for either one of these traits. Nyquist-sampling of the Fourier amplitudes is sufficient to image single particles of finite size at any resolution. Atomic resolution data is routinely used to image…
▽ More
Any object on earth has two fundamental properties: it is finite, and it is made of atoms. Structural information about an object can be obtained from diffraction amplitude measurements that account for either one of these traits. Nyquist-sampling of the Fourier amplitudes is sufficient to image single particles of finite size at any resolution. Atomic resolution data is routinely used to image molecules replicated in a crystal structure. Here we report an algorithm that requires neither information, but uses the fact that an image of a natural object is compressible. Intended applications include tomographic diffractive imaging, crystallography, powder diffraction, small angle x-ray scattering and random Fourier amplitude measurements.
△ Less
Submitted 11 September, 2008;
originally announced September 2008.
-
X-Ray Powder Diffraction from Sub-Micron Crystals of Photosystem-1 Membrane Protein
Authors:
D. A. Shapiro,
D. DePonte,
R. B. Doak,
P. Fromme,
G. Hembree,
M. Hunter,
S. Marchesini,
K. Schmidt,
D. Starodub,
U. Weierstall,
H. Chapman,
J. Spence
Abstract:
We demonstrate that powder diffraction data can be collected from sub-micron crystals of a mbrane protein with nearly two orders of magnitude more atoms than the molecules commonly used for powder diffraction. The crystals of photosystem-1 protein were size-selected using a 500 nm pore- size filter and delivered to a soft x-ray beam with a photon energy of 1.5 keV using a dynamically focused mic…
▽ More
We demonstrate that powder diffraction data can be collected from sub-micron crystals of a mbrane protein with nearly two orders of magnitude more atoms than the molecules commonly used for powder diffraction. The crystals of photosystem-1 protein were size-selected using a 500 nm pore- size filter and delivered to a soft x-ray beam with a photon energy of 1.5 keV using a dynamically focused micro-jet developed for the serial crystallography experiment at beamline 9.0.1. The 10-micron jet places many such randomly oriented crystals in the x-ray beam simultaneously resulting in a powder diffraction pattern which extends to 28 angstrom resolution with just 200 seconds of x-ray exposure. The use of the jet for particle delivery allows for a thin sample, appropriate for the soft x-rays used, and continuously refreshes the crystals so that radiation damage is not possible. The small size of the crystals requires the use of lower energy photons for increased scattering strength and increased spacing between powder rings. The powder patterns obtained in this way, from abundant nano-crystals, could be used to provide low resolution molecular envelopes if phased using techniques such as compressive sensing which do not require atomic resolution data. The results also serve to test our aerojet injector system, with future application to femtosecond diffraction in Free Electron X-ray Laser schemes, and for Serial Crystallography using a single-file beam of aligned hydrated molecules.
△ Less
Submitted 25 September, 2008; v1 submitted 27 March, 2008;
originally announced March 2008.
-
Ultrafast, ultrabright, X-ray holography using a uniformly-redundant array
Authors:
S. Marchesini,
S. Boutet,
A. E. Sakdinawat,
M. J. Bogan,
S. Bajt,
A. Barty,
H. N. Chapman,
M. Frank,
S. P. Hau-Riege,
A. Szoke,
C. Cui,
M. R. Howells,
D. A. Shapiro,
J. C. H. Spence,
J. W. Shaevitz,
J. Y. Lee,
J. Hajdu,
M. M. Seibert
Abstract:
Advances in the development of free-electron lasers offer the realistic prospect of high-resolution imaging to study the nanoworld on the time-scale of atomic motions. We identify X-ray Fourier Transform holography, (FTH) as a promising but, so far, inefficient scheme to do this. We show that a uniformly redundant array (URA) placed next to the sample, multiplies the efficiency of X-ray FTH by m…
▽ More
Advances in the development of free-electron lasers offer the realistic prospect of high-resolution imaging to study the nanoworld on the time-scale of atomic motions. We identify X-ray Fourier Transform holography, (FTH) as a promising but, so far, inefficient scheme to do this. We show that a uniformly redundant array (URA) placed next to the sample, multiplies the efficiency of X-ray FTH by more than one thousand (approaching that of a perfect lens) and provides holographic images with both amplitude- and phase-contrast information. The experiments reported here demonstrate this concept by imaging a nano-fabricated object at a synchrotron source, and a bacterial cell at a soft X-ray free-electron-laser, where illumination by a single 15 fs pulse was successfully used in producing the holographic image. We expect with upcoming hard X-ray lasers to achieve considerably higher spatial resolution and to obtain ultrafast movies of excited states of matter.
△ Less
Submitted 8 February, 2008; v1 submitted 31 January, 2008;
originally announced January 2008.
-
Three-dimensional coherent X-ray diffraction imaging of a ceramic nanofoam: determination of structural deformation mechanisms
Authors:
A. Barty,
S. Marchesini,
H. N. Chapman,
C. Cui,
M. R. Howells,
D. A. Shapiro,
A. M. Minor,
J. C. H. Spence,
U. Weierstall,
J. Ilavsky,
A. Noy,
S. P. Hau-Riege,
A. B. Artyukhin,
T. Baumann,
T. Willey,
J. Stolken,
T. van Buuren,
J. H. Kinney
Abstract:
Ultra-low density polymers, metals, and ceramic nanofoams are valued for their high strength-to-weight ratio, high surface area and insulating properties ascribed to their structural geometry. We obtain the labrynthine internal structure of a tantalum oxide nanofoam by X-ray diffractive imaging. Finite element analysis from the structure reveals mechanical properties consistent with bulk samples…
▽ More
Ultra-low density polymers, metals, and ceramic nanofoams are valued for their high strength-to-weight ratio, high surface area and insulating properties ascribed to their structural geometry. We obtain the labrynthine internal structure of a tantalum oxide nanofoam by X-ray diffractive imaging. Finite element analysis from the structure reveals mechanical properties consistent with bulk samples and with a diffusion limited cluster aggregation model, while excess mass on the nodes discounts the dangling fragments hypothesis of percolation theory.
△ Less
Submitted 25 June, 2008; v1 submitted 30 August, 2007;
originally announced August 2007.
-
Phase retrieval and saddle-point optimization
Authors:
S. Marchesini
Abstract:
Iterative algorithms with feedback are amongst the most powerful and versatile optimization methods for phase retrieval. Among these, the hybrid input-output algorithm has demonstrated practical solutions to giga-element nonlinear phase retrieval problems, escaping local minima and producing images at resolutions beyond the capabilities of lens-based optical methods. Here, the input-output itera…
▽ More
Iterative algorithms with feedback are amongst the most powerful and versatile optimization methods for phase retrieval. Among these, the hybrid input-output algorithm has demonstrated practical solutions to giga-element nonlinear phase retrieval problems, escaping local minima and producing images at resolutions beyond the capabilities of lens-based optical methods. Here, the input-output iteration is improved by a lower dimensional subspace saddle-point optimization.
△ Less
Submitted 6 June, 2007; v1 submitted 24 November, 2006;
originally announced November 2006.
-
Femtosecond Diffractive Imaging with a Soft-X-ray Free-Electron Laser
Authors:
Henry N. Chapman,
Anton Barty,
Michael J. Bogan,
Sebastien Boutet,
Matthias Frank,
Stefan P. Hau-Riege,
Stefano Marchesini,
Bruce W. Woods,
Sasa Bajt,
W. Henry Benner,
Richard A. London,
Elke Plonjes,
Marion Kuhlmann,
Rolf Treusch,
Stefan Dusterer,
Thomas Tschentscher,
Jochen R. Schneider,
Eberhard Spiller,
Thomas Moller,
Christoph Bostedt,
Matthias Hoener,
David A. Shapiro,
Keith O. Hodgson,
David van der Spoel,
Florian Burmeister
, et al. (9 additional authors not shown)
Abstract:
Theory predicts that with an ultrashort and extremely bright coherent X-ray pulse, a single diffraction pattern may be recorded from a large macromolecule, a virus, or a cell before the sample explodes and turns into a plasma. Here we report the first experimental demonstration of this principle using the FLASH soft X-ray free-electron laser. An intense 25 fs, 4 10^13 W/cm^2 pulse, containing 10…
▽ More
Theory predicts that with an ultrashort and extremely bright coherent X-ray pulse, a single diffraction pattern may be recorded from a large macromolecule, a virus, or a cell before the sample explodes and turns into a plasma. Here we report the first experimental demonstration of this principle using the FLASH soft X-ray free-electron laser. An intense 25 fs, 4 10^13 W/cm^2 pulse, containing 10^12 photons at 32 nm wavelength, produced a coherent diffraction pattern from a nano-structured non-periodic object, before destroying it at 60,000 K. A novel X-ray camera assured single photon detection sensitivity by filtering out parasitic scattering and plasma radiation. The reconstructed image, obtained directly from the coherent pattern by phase retrieval through oversampling, shows no measurable damage, and extends to diffraction-limited resolution. A three-dimensional data set may be assembled from such images when copies of a reproducible sample are exposed to the beam one by one.
△ Less
Submitted 10 October, 2006; v1 submitted 7 October, 2006;
originally announced October 2006.
-
A unified evaluation of iterative projection algorithms for phase retrieval
Authors:
S. Marchesini
Abstract:
Iterative projection algorithms are successfully being used as a substitute of lenses to recombine, numerically rather than optically, light scattered by illuminated objects. Images obtained computationally allow aberration-free diffraction-limited imaging and the possibility of using radiation for which no lenses exist. The challenge of this imaging technique is transfered from the lenses to th…
▽ More
Iterative projection algorithms are successfully being used as a substitute of lenses to recombine, numerically rather than optically, light scattered by illuminated objects. Images obtained computationally allow aberration-free diffraction-limited imaging and the possibility of using radiation for which no lenses exist. The challenge of this imaging technique is transfered from the lenses to the algorithms. We evaluate these new computational ``instruments'' developed for the phase retrieval problem, and discuss acceleration strategies.
△ Less
Submitted 22 December, 2006; v1 submitted 24 March, 2006;
originally announced March 2006.
-
Phase Aberrations in Diffraction Microscopy
Authors:
S. Marchesini,
H. N. Chapman,
A. Barty,
C. Cui,
M. R. Howells,
J. C. H. Spence,
U. Weierstall,
A. M. Minor
Abstract:
In coherent X-ray diffraction microscopy the diffraction pattern generated by a sample illuminated with coherent x-rays is recorded, and a computer algorithm recovers the unmeasured phases to synthesize an image. By avoiding the use of a lens the resolution is limited, in principle, only by the largest scattering angles recorded. However, the imaging task is shifted from the experiment to the co…
▽ More
In coherent X-ray diffraction microscopy the diffraction pattern generated by a sample illuminated with coherent x-rays is recorded, and a computer algorithm recovers the unmeasured phases to synthesize an image. By avoiding the use of a lens the resolution is limited, in principle, only by the largest scattering angles recorded. However, the imaging task is shifted from the experiment to the computer, and the algorithm's ability to recover meaningful images in the presence of noise and limited prior knowledge may produce aberrations in the reconstructed image. We analyze the low order aberrations produced by our phase retrieval algorithms. We present two methods to improve the accuracy and stability of reconstructions.
△ Less
Submitted 6 October, 2005; v1 submitted 4 October, 2005;
originally announced October 2005.
-
Progress in Three-Dimensional Coherent X-Ray Diffraction Imaging
Authors:
S. Marchesini,
H. N. Chapman,
A. Barty,
A. Noy,
S. P. Hau-Riege,
J. M. Kinney,
C. Cui,
M. R. Howells,
R. Rosen,
J. C. H. Spence,
U. Weierstall,
D. Shapiro,
T. Beetz,
C. Jacobsen,
E. Lima,
A. M. Minor,
H. He
Abstract:
The Fourier inversion of phased coherent diffraction patterns offers images without the resolution and depth-of-focus limitations of lens-based tomographic systems. We report on our recent experimental images inverted using recent developments in phase retrieval algorithms, and summarize efforts that led to these accomplishments. These include ab-initio reconstruction of a two-dimensional test p…
▽ More
The Fourier inversion of phased coherent diffraction patterns offers images without the resolution and depth-of-focus limitations of lens-based tomographic systems. We report on our recent experimental images inverted using recent developments in phase retrieval algorithms, and summarize efforts that led to these accomplishments. These include ab-initio reconstruction of a two-dimensional test pattern, infinite depth of focus image of a thick object, and its high-resolution (~10 nm resolution) three-dimensional image. Developments on the structural imaging of low density aerogel samples are discussed.
△ Less
Submitted 6 October, 2005; v1 submitted 4 October, 2005;
originally announced October 2005.
-
High-resolution ab initio three-dimensional X-ray diffraction microscopy
Authors:
H. N. Chapman,
A. Barty,
S. Marchesini,
A. Noy,
C. Cui,
M. R. Howells,
R. Rosen,
H. He,
J. C. H. Spence,
U. Weierstall,
T. Beetz,
C. Jacobsen,
D. Shapiro
Abstract:
Coherent X-ray diffraction microscopy is a method of imaging non-periodic isolated objects at resolutions only limited, in principle, by the largest scattering angles recorded. We demonstrate X-ray diffraction imaging with high resolution in all three dimensions, as determined by a quantitative analysis of the reconstructed volume images. These images are retrieved from the 3D diffraction data u…
▽ More
Coherent X-ray diffraction microscopy is a method of imaging non-periodic isolated objects at resolutions only limited, in principle, by the largest scattering angles recorded. We demonstrate X-ray diffraction imaging with high resolution in all three dimensions, as determined by a quantitative analysis of the reconstructed volume images. These images are retrieved from the 3D diffraction data using no a priori knowledge about the shape or composition of the object, which has never before been demonstrated on a non-periodic object. We also construct 2D images of thick objects with infinite depth of focus (without loss of transverse spatial resolution). These methods can be used to image biological and materials science samples at high resolution using X-ray undulator radiation, and establishes the techniques to be used in atomic-resolution ultrafast imaging at X-ray free-electron laser sources.
△ Less
Submitted 8 September, 2005;
originally announced September 2005.
-
An assessment of the resolution limitation due to radiation-damage in x-ray diffraction microscopy
Authors:
M. R. Howells,
T. Beetz,
H. N. Chapman,
C. Cui,
J. M. Holton,
C. J. Jacobsen,
J. Kirz E. Lima,
S. Marchesini,
H. Miao,
D. Sayre,
D. A. Shapiro,
J. C. H. Spence
Abstract:
X-ray diffraction microscopy (XDM) is a new form of x-ray imaging that is being practiced at several third-generation synchrotron-radiation x-ray facilities. Although only five years have elapsed since the technique was first introduced, it has made rapid progress in demonstrating high-resolution threedimensional imaging and promises few-nm resolution with much larger samples than can be imaged…
▽ More
X-ray diffraction microscopy (XDM) is a new form of x-ray imaging that is being practiced at several third-generation synchrotron-radiation x-ray facilities. Although only five years have elapsed since the technique was first introduced, it has made rapid progress in demonstrating high-resolution threedimensional imaging and promises few-nm resolution with much larger samples than can be imaged in the transmission electron microscope. Both life- and materials-science applications of XDM are intended, and it is expected that the principal limitation to resolution will be radiation damage for life science and the coherent power of available x-ray sources for material science. In this paper we address the question of the role of radiation damage. We use a statistical analysis based on the so-called "dose fractionation theorem" of Hegerl and Hoppe to calculate the dose needed to make an image of a lifescience sample by XDM with a given resolution. We conclude that the needed dose scales with the inverse fourth power of the resolution and present experimental evidence to support this finding. To determine the maximum tolerable dose we have assembled a number of data taken from the literature plus some measurements of our own which cover ranges of resolution that are not well covered by reports in the literature. The tentative conclusion of this study is that XDM should be able to image frozen-hydrated protein samples at a resolution of about 10 nm with "Rose-criterion" image quality.
△ Less
Submitted 10 February, 2005;
originally announced February 2005.
-
Correction of non-linearity effects in detectors for electron spectroscopy
Authors:
N. Mannella,
S. Marchesini,
A. W. Kay,
A. Nambu,
T. Gresch,
S. -H. Yang,
B. S. Mun,
J. M. Bussat,
A. Rosenhahn,
C. S. Fadley
Abstract:
Using photoemission intensities and a detection system employed by many groups in the electron spectroscopy community as an example, we have quantitatively characterized and corrected detector non-linearity effects over the full dynamic range of the system. Non-linearity effects are found to be important whenever measuring relative peak intensities accurately is important, even in the low-countr…
▽ More
Using photoemission intensities and a detection system employed by many groups in the electron spectroscopy community as an example, we have quantitatively characterized and corrected detector non-linearity effects over the full dynamic range of the system. Non-linearity effects are found to be important whenever measuring relative peak intensities accurately is important, even in the low-countrate regime. This includes, for example, performing quantitative analyses for surface contaminants or sample bulk stoichiometries, where the peak intensities involved can differ by one or two orders of magnitude, and thus could occupy a significant portion of the detector dynamic range. Two successful procedures for correcting non-linearity effects are presented. The first one yields directly the detector efficiency by measuring a flat-background reference intensity as a function of incident x-ray flux, while the second one determines the detector response from a least-squares analysis of broad-scan survey spectra at different incident x-ray fluxes. Although we have used one spectrometer and detection system as an example, these methodologies should be useful for many other cases.
△ Less
Submitted 17 December, 2004;
originally announced December 2004.
-
Use of extended and prepared reference objects in experimental Fourier transform X-ray holography
Authors:
H. He,
M. Howells,
S. Marchesini,
H. N. Chapman,
U. Weierstall,
H. A. Padmore,
J. C. H. Spence
Abstract:
The use of one or more gold nanoballs as reference objects for Fourier Transform holography (FTH) is analysed using experimental soft X-ray diffraction from objects consisting of separated clusters of these balls. The holograms are deconvoluted against ball reference objects to invert to images, in combination with a Wiener filter to control noise. A resolution of ~30nm, smaller than one ball, i…
▽ More
The use of one or more gold nanoballs as reference objects for Fourier Transform holography (FTH) is analysed using experimental soft X-ray diffraction from objects consisting of separated clusters of these balls. The holograms are deconvoluted against ball reference objects to invert to images, in combination with a Wiener filter to control noise. A resolution of ~30nm, smaller than one ball, is obtained even if a large cluster of balls is used as the reference, giving the best resolution yet obtained by X-ray FTH. Methods of dealing with missing data due to a beamstop are discussed. Practical prepared objects which satisfy the FTH condition are suggested, and methods of forming them described.
△ Less
Submitted 7 May, 2004;
originally announced May 2004.
-
Benchmarking Iterative Projection Algorithms for Phase Retrieval
Authors:
S. Marchesini
Abstract:
Iterative projection algorithms for phase retrieval are tested on two simple toy models. The result provides useful insights in the behavior of these algorithms.
Iterative projection algorithms for phase retrieval are tested on two simple toy models. The result provides useful insights in the behavior of these algorithms.
△ Less
Submitted 19 April, 2004;
originally announced April 2004.
-
Phasing diffuse scattering. Application of the SIR2002 algorithm to the non-crystallographic phase problem
Authors:
B. Carrozzini,
G. L. Cascarano,
L. De Caro,
C. Giacovazzo,
S. Marchesini,
H. N. Chapman,
H. He,
M. R. Howells,
J. S. Wu,
U. Weierstall,
J. C. H. Spence
Abstract:
A new phasing algorithm has been used to determine the phases of diffuse elastic X-ray scattering from a non-periodic array of gold balls of 50 nm diameter. Two-dimensional real-space images, showing the charge-density distribution of the balls, have been reconstructed at 50 nm resolution from transmission diffraction patterns recorded at 550 eV energy. The reconstructed image fits well with sca…
▽ More
A new phasing algorithm has been used to determine the phases of diffuse elastic X-ray scattering from a non-periodic array of gold balls of 50 nm diameter. Two-dimensional real-space images, showing the charge-density distribution of the balls, have been reconstructed at 50 nm resolution from transmission diffraction patterns recorded at 550 eV energy. The reconstructed image fits well with scanning electron microscope (SEM) image of the same sample. The algorithm, which uses only the density modification portion of the SIR2002 program, is compared with the results obtained via the Gerchberg-Saxton-Fienup HIO algorithm. In this way the relationship between density modification in crystallography and the HiO algorithm used in signal and image processing is elucidated.
△ Less
Submitted 15 April, 2004;
originally announced April 2004.
-
SPEDEN: Reconstructing single particles from their diffraction patterns
Authors:
S. P. Hau-Riege,
H. Szoke,
H. N. Chapman,
A. Szoke,
S. Marchesini,
A. Noy,
H. He,
M. R. Howells,
U. Weierstall,
J. C. H. Spence
Abstract:
Speden is a computer program that reconstructs the electron density of single particles from their x-ray diffraction patterns, using a single-particle adaptation of the Holographic Method in crystallography. (Szoke, A., Szoke, H., and Somoza, J.R., 1997. Acta Cryst. A53, 291-313.) The method, like its parent, is unique that it does not rely on ``back'' transformation from the diffraction pattern…
▽ More
Speden is a computer program that reconstructs the electron density of single particles from their x-ray diffraction patterns, using a single-particle adaptation of the Holographic Method in crystallography. (Szoke, A., Szoke, H., and Somoza, J.R., 1997. Acta Cryst. A53, 291-313.) The method, like its parent, is unique that it does not rely on ``back'' transformation from the diffraction pattern into real space and on interpolation within measured data. It is designed to deal successfully with sparse, irregular, incomplete and noisy data. It is also designed to use prior information for ensuring sensible results and for reliable convergence. This article describes the theoretical basis for the reconstruction algorithm, its implementation and quantitative results of tests on synthetic and experimentally obtained data. The program could be used for determining the structure of radiation tolerant samples and, eventually, of large biological molecular structures without the need for crystallization.
△ Less
Submitted 18 March, 2004;
originally announced March 2004.
-
Spatial Coherence of Synchrotron Radiation
Authors:
R. Coisson,
S. Marchesini
Abstract:
Theory and measurement of spatial coherence of synchrotron radiation beams are briefly reviewed. Emphasis is given to simple relationships between electron beam characteristics and far field properties of the light beam.
Theory and measurement of spatial coherence of synchrotron radiation beams are briefly reviewed. Emphasis is given to simple relationships between electron beam characteristics and far field properties of the light beam.
△ Less
Submitted 7 January, 2004;
originally announced January 2004.
-
Direct Observation of High-Temperature Polaronic Behavior In Colossal Magnetoresistive Manganites
Authors:
N. Mannella,
A. Rosenhahn,
C. H. Booth,
S. Marchesini,
B. S. Mun,
S. -H. Yang,
K. Ibrahim,
Y. Tomioka,
C. S. Fadley
Abstract:
The temperature dependence of the electronic and atomic structure of the colossal magnetoresistive oxides $La_{1-x}Sr_{x}MnO_{3}$ (x = 0.3, 0.4) has been studied using core and valence level photoemission, x-ray absorption and emission, and extended x-ray absorption fine structure spectroscopy. A dramatic and reversible change of the electronic structure is observed on crossing the Curie tempera…
▽ More
The temperature dependence of the electronic and atomic structure of the colossal magnetoresistive oxides $La_{1-x}Sr_{x}MnO_{3}$ (x = 0.3, 0.4) has been studied using core and valence level photoemission, x-ray absorption and emission, and extended x-ray absorption fine structure spectroscopy. A dramatic and reversible change of the electronic structure is observed on crossing the Curie temperature, including charge localization and spin moment increase of Mn, together with Jahn-Teller distortions, both signatures of polaron formation. Our data are also consistent with a phase-separation scenario.
△ Less
Submitted 15 January, 2004; v1 submitted 15 January, 2004;
originally announced January 2004.
-
Coherent X-ray Diffractive Imaging; applications and limitations
Authors:
S. Marchesini,
H. N. Chapman,
S. P. Hau-Riege,
R. A. London,
A. Szoke,
H. He,
M. R. Howells,
H. Padmore,
R. Rosen,
J. C. H. Spence,
U. Weierstall
Abstract:
The inversion of a diffraction pattern offers aberration-free diffraction-limited 3D images without the resolution and depth-of-field limitations of lens-based tomographic systems, the only limitation being radiation damage. We review our experimental results, discuss the fundamental limits of this technique and future plans.
The inversion of a diffraction pattern offers aberration-free diffraction-limited 3D images without the resolution and depth-of-field limitations of lens-based tomographic systems, the only limitation being radiation damage. We review our experimental results, discuss the fundamental limits of this technique and future plans.
△ Less
Submitted 15 August, 2003;
originally announced August 2003.
-
Holographic Methods as Local Probes of the Atomic Order in Solids
Authors:
G. Faigel,
M Tegze,
G. Bortel,
Z. Jurek,
S. Marchesini,
M. Belakhovskyo,
A. Simionovici
Abstract:
In the last fifteen years several techniques based on the holographic principle have been developed for the study of the 3D local order in solids. These methods use various particles: electrons, hard x-ray photons, gamma photons, or neutrons to image the atoms. Although the practical realisation of the various imaging experiments is very different, there is a common thread; the use of inside ref…
▽ More
In the last fifteen years several techniques based on the holographic principle have been developed for the study of the 3D local order in solids. These methods use various particles: electrons, hard x-ray photons, gamma photons, or neutrons to image the atoms. Although the practical realisation of the various imaging experiments is very different, there is a common thread; the use of inside reference points for holographic imaging. In this paper we outline the basics of atomic resolution holography using inside reference points, especially concentrating to the hard x-ray case. Further, we outline the experimental requirements and what has been practically realized in the last decade. At last we give examples of applications and future perspectives.
△ Less
Submitted 19 August, 2003;
originally announced August 2003.
-
Differential Photoelectron Holography: A New Approach for Three-Dimensional Atomic Imaging
Authors:
S. Omori,
Y. Nihei,
E. Rotenberg,
J. D. Denlinger,
S. Marchesini,
S. D. Kevan,
B. P. Tonner,
M. A. Van Hove,
C. S. Fadley
Abstract:
We propose differential holography as a method to overcome the long-standing forward-scattering problem in photoelectron holography and related techniques for the three-dimensional imaging of atoms. Atomic images reconstructed from experimental and theoretical Cu 3p holograms from Cu(001) demonstrate that this method suppresses strong forward-scattering effects so as to yield more accurate three…
▽ More
We propose differential holography as a method to overcome the long-standing forward-scattering problem in photoelectron holography and related techniques for the three-dimensional imaging of atoms. Atomic images reconstructed from experimental and theoretical Cu 3p holograms from Cu(001) demonstrate that this method suppresses strong forward-scattering effects so as to yield more accurate three-dimensional images of side- and back-scattering atoms.
△ Less
Submitted 28 July, 2003;
originally announced July 2003.
-
X-ray image reconstruction from a diffraction pattern alone
Authors:
S. Marchesini,
H. He,
H. N. Chapman,
S. P. Hau-Riege,
A. Noy,
M. R. Howells,
U. Weierstall,
J. C. H. Spence
Abstract:
A solution to the inversion problem of scattering would offer aberration-free diffraction-limited 3D images without the resolution and depth-of-field limitations of lens-based tomographic systems. Powerful algorithms are increasingly being used to act as lenses to form such images. Current image reconstruction methods, however, require the knowledge of the shape of the object and the low spatial…
▽ More
A solution to the inversion problem of scattering would offer aberration-free diffraction-limited 3D images without the resolution and depth-of-field limitations of lens-based tomographic systems. Powerful algorithms are increasingly being used to act as lenses to form such images. Current image reconstruction methods, however, require the knowledge of the shape of the object and the low spatial frequencies unavoidably lost in experiments. Diffractive imaging has thus previously been used to increase the resolution of images obtained by other means. We demonstrate experimentally here a new inversion method, which reconstructs the image of the object without the need for any such prior knowledge.
△ Less
Submitted 12 September, 2003; v1 submitted 24 June, 2003;
originally announced June 2003.
-
X-ray fluorescence holography: going beyond the diffraction limit
Authors:
S. Marchesini,
C. S. Fadley,
F. J. Garcia de Abajo
Abstract:
X-ray fluorescence holography (XFH) is a method for obtaining diffraction-limited images of the local atomic structure around a given type of emitter. The reconstructed wave-field represents a distorted image of the scatterer electron density distribution; {\it i.e.} it is a convolution of the charge density distribution with a point spread function characteristic of the measurement. We here con…
▽ More
X-ray fluorescence holography (XFH) is a method for obtaining diffraction-limited images of the local atomic structure around a given type of emitter. The reconstructed wave-field represents a distorted image of the scatterer electron density distribution; {\it i.e.} it is a convolution of the charge density distribution with a point spread function characteristic of the measurement. We here consider several methods for the iterative deconvolution of such XFH holograms, and via theoretical simulations evaluate them from the point of view of going beyond the diffraction limit so as to image the electron charge density. Promising results for future applications are found for certain methods, and other possible image-enhancement techniques are also discussed.
△ Less
Submitted 4 June, 2002; v1 submitted 5 February, 2002;
originally announced February 2002.
-
Holographic analysis of diffraction structure factors
Authors:
S. Marchesini,
N. Mannella,
C. S. Fadley,
M. A. Van Hove,
J. J. Bucher,
D. K. Shuh,
L. Fabris,
M. H. Press,
M. W. West,
W. C. Stolte,
Z. Hussain
Abstract:
We combine the theory of inside-source/inside-detector x-ray fluorescence holography and Kossel lines/x ray standing waves in kinematic approximation to directly obtain the phases of the diffraction structure factors. The influence of Kossel lines and standing waves on holography is also discussed. We obtain partial phase determination from experimental data obtaining the sign of the real part o…
▽ More
We combine the theory of inside-source/inside-detector x-ray fluorescence holography and Kossel lines/x ray standing waves in kinematic approximation to directly obtain the phases of the diffraction structure factors. The influence of Kossel lines and standing waves on holography is also discussed. We obtain partial phase determination from experimental data obtaining the sign of the real part of the structure factor for several reciprocal lattice vectors of a vanadium crystal.
△ Less
Submitted 2 May, 2002; v1 submitted 4 February, 2002;
originally announced February 2002.