Search

Signal-to-noise ratio and spatial resolution are quantitatively analysed in the context of in-line (propagation based) X-ray phase-contrast imaging. It is known that free-space propagation of a coherent X-ray beam from the imaged object to the detector plane, followed by phase retrieval in accordance with Paganin's method, can increase the signal-to-noise in the resultant images without deteriorating the spatial resolution. This results in violation of the noise-resolution uncertainty principle and demonstrates `unreasonable' effectiveness of the method. On the other hand, when the process of free-space propagation is performed in software, using the detected intensity distribution in the object plane, it cannot reproduce the same effectiveness, due to the amplification of photon shot noise. Here, it is shown that the performance of Paganin's method is determined by just two dimensionless parameters: the Fresnel number and the ratio of the real decrement to the imaginary part of the refractive index of the imaged object. The relevant theoretical analysis is performed first, followed by computer simulations and then by a brief test using experimental images collected at a synchrotron beamline. More extensive experimental tests will be presented in the second part of this paper., (open access.)

Double-slit interference experiments using monochromatic hard X-rays with the energy of 25 keV are presented. The experiments were performed at a synchrotron source with a distance of 110 m between the interferometer and the detector to produce an interference pattern with a sufficiently broad period that could be adequately sampled by a photon-counting detector with 75 micrometre pixels. In the single-particle version of the experiment, over one million image frames with a single registered photon in each one were collected. The sum of these frames showed a clear presence of the interference pattern with the expected period. Subsequent analysis provided an objective estimation of the minimal number of detected photons required to determine, in accordance with the Rose criterion, the presence of the photon interference. Apart from a general theoretical interest, these investigations were aimed at exploring the possibility of medical X-ray phase-contrast imaging in photon-counting regime at minimal radiation doses.

A unified method for three-dimensional reconstruction of objects from transmission images collected at multiple illumination directions is described. The method may be applicable to experimental conditions relevant to absorption-based, phase-contrast, or diffraction imaging using x rays, electrons, and other forms of penetrating radiation or matter waves. Both the phase retrieval (also known as contrast transfer function correction) and the effect of Ewald sphere curvature (in the cases with a shallow depth of field and significant in-object diffraction) are incorporated in the proposed algorithm and can be taken into account. Multiple scattering is not treated explicitly but can be mitigated as a result of angular averaging that constitutes an essential feature of the method. The corresponding numerical algorithm is based on three-dimensional gridding which allows for fast computational implementation, including a straightforward parallelization. The algorithm can be used with any scanning geometry involving plane-wave illumination. A software code implementing the proposed algorithm has been developed, tested on simulated and experimental image data, and made publicly available.

Breast microcalcifications are an important primary radiological indicator of breast cancer. However, microcalcification classification and diagnosis may be still challenging for radiologists due to limitations of the standard 2D mammography technique, including spatial and contrast resolution. In this study, we propose an approach to improve the detection of microcalcifications in propagation-based phase-contrast X-ray computed tomography of breast tissues. Five fresh mastectomies containing microcalcifications were scanned at different X-ray energies and radiation doses using synchrotron radiation. Both bright-field (i.e. conventional phase-retrieved images) and dark-field images were extracted from the same data sets using different image processing methods. A quantitative analysis was performed in terms of visibility and contrast-to-noise ratio of microcalcifications. The results show that while the signal-to-noise and the contrast-to-noise ratios are lower, the visibility of the microcalcifications is more than two times higher in the dark-field images compared to the bright-field images. Dark-field images have also provided more accurate information about the size and shape of the microcalcifications.

A method is proposed for high-resolution, three-dimensional reconstruction of internal structures of objects from planar transmission images. The described approach can be used with any form of radiation or matter waves, in principle, provided that the depth of field is smaller than the thickness of the sample. The physical optics basis for the method is elucidated, and the reconstruction algorithm is presented in detail. A simulated example demonstrates an application of the method to three-dimensional electron transmission imaging of a nanoparticle under realistic radiation dose and spatial resolution constraints. It is envisaged that the method can be applicable in high-resolution transmission electron microscopy, soft x-ray microscopy, ultrasound imaging, and other areas.

A method for three-dimensional reconstruction of objects from defocused images collected at multiple illumination directions in high-resolution transmission electron microscopy is presented. The method effectively corrects for the Ewald sphere curvature by taking into account the in-particle propagation of the electron beam. Numerical simulations demonstrate that the proposed method is capable of accurately reconstructing biological molecules or nanoparticles from high-resolution defocused images under conditions achievable in single-particle electron cryo-microscopy or electron tomography with realistic radiation doses, non-trivial aberrations, multiple scattering, and other experimentally relevant factors. The physics of the method is based on the well-known Diffraction Tomography formalism, but with the phase-retrieval step modified to include a conjugation of the phase (i.e., multiplication of the phase by a negative constant). At each illumination direction, numerically backpropagating the beam with the conjugated phase produces maximum contrast at the location of individual atoms in the molecule or nanoparticle. The resultant algorithm, Conjugated Holographic Reconstruction, can potentially be incorporated into established software tools for single-particle analysis, such as, for example, RELION or FREALIGN, in place of the conventional contrast transfer function correction procedure, in order to account for the Ewald sphere curvature and improve the spatial resolution of the three-dimensional reconstruction.

It has been argued that in atomic-resolution transmission electron microscopy (TEM) of sparse weakly scattering structures, such as small biological molecules, multiple electron scattering usually has only a small effect, while the in-molecule Fresnel diffraction can be significant due to the intrinsically shallow depth of focus. These facts suggest that the three-dimensional reconstruction of such structures from defocus image series collected at multiple rotational orientations of a molecule can be effectively performed for each atom separately, using the incoherent first Born approximation. The corresponding reconstruction method, termed here Differential Holographic Tomography, is developed theoretically and demonstrated computationally on several numerical models of biological molecules. It is shown that the method is capable of accurate reconstruction of the locations of atoms in a molecule from TEM data collected at a small number of random orientations of the molecule, with one or more defocus images per orientation. Possible applications to cryogenic electron microscopy and other areas are briefly discussed., (Copyright © 2021 Elsevier B.V. All rights reserved.)

Purpose: Breast cancer is the most common cancer in women in developing and developed countries and is responsible for 15% of women's cancer deaths worldwide. Conventional absorption-based breast imaging techniques lack sufficient contrast for comprehensive diagnosis. Propagation-based phase-contrast computed tomography (PB-CT) is a developing technique that exploits a more contrast-sensitive property of x-rays: x-ray refraction. X-ray absorption, refraction, and contrast-to-noise in the corresponding images depend on the x-ray energy used, for the same/fixed radiation dose. The aim of this paper is to explore the relationship between x-ray energy and radiological image quality in PB-CT imaging. Approach: Thirty-nine mastectomy samples were scanned at the imaging and medical beamline at the Australian Synchrotron. Samples were scanned at various x-ray energies of 26, 28, 30, 32, 34, and 60 keV using a Hamamatsu Flat Panel detector at the same object-to-detector distance of 6 m and mean glandular dose of 4 mGy. A total of 132 image sets were produced for analysis. Seven observers rated PB-CT images against absorption-based CT (AB-CT) images of the same samples on a five-point scale. A visual grading characteristics (VGC) study was used to determine the difference in image quality. Results: PB-CT images produced at 28, 30, 32, and 34 keV x-ray energies demonstrated statistically significant higher image quality than reference AB-CT images. The optimum x-ray energy, 30 keV, displayed the largest area under the curve ( AUC VGC ) of 0.754 ( p = 0.009 ). This was followed by 32 keV ( AUC VGC = 0.731 , p ≤ 0.001 ), 34 keV ( AUC VGC = 0.723 , p ≤ 0.001 ), and 28 keV ( AUC VGC = 0.654 , p = 0.015 ). Conclusions: An optimum energy range (around 30 keV) in the PB-CT technique allows for higher image quality at a dose comparable to conventional mammographic techniques. This results in improved radiological image quality compared with conventional techniques, which may ultimately lead to higher diagnostic efficacy and a reduction in breast cancer mortalities., (© 2021 Society of Photo-Optical Instrumentation Engineers (SPIE).)

Rationale and Objectives: Propagation-based phase-contrast CT (PB-CT) is an advanced X-ray imaging technology that exploits both refraction and absorption of the transmitted X-ray beam. This study was aimed at optimizing the experimental conditions of PB-CT for breast cancer imaging and examined its performance relative to conventional absorption-based CT (AB-CT) in terms of image quality and radiation dose., Materials and Methods: Surgically excised breast mastectomy specimens (n = 12) were scanned using both PB-CT and AB-CT techniques under varying imaging conditions. To evaluate the radiological image quality, visual grading characteristics (VGC) analysis was used in which 11 breast specialist radiologists compared the overall image quality of PB-CT images with respect to the corresponding AB-CT images. The area under the VGC curve was calculated to measure the differences between PB-CT and AB-CT images., Results: The highest radiological quality was obtained for PB-CT images using a 32 keV energy X-ray beam and by applying the Homogeneous Transport of Intensity Equation phase retrieval with the value of its parameter γ set to one-half of the theoretically optimal value for the given materials. Using these optimized conditions, the image quality of PB-CT images obtained at 4 mGy and 2 mGy mean glandular dose was significantly higher than AB-CT images at 4 mGy (AUC VGC  = 0.901, p = 0.001 and AUC VGC  = 0.819, p = 0.011, respectively)., Conclusion: PB-CT achieves a higher radiological image quality compared to AB-CT even at a considerably lower mean glandular dose. Successful translation of the PB-CT technique for breast cancer imaging can potentially result in improved breast cancer diagnosis., (Crown Copyright © 2020. Published by Elsevier Inc. All rights reserved.)

A method for extracting the dark-field signal in propagation-based phase-contrast imaging is proposed. In the case of objects consisting predominantly of a single material, or several different materials with similar ratios of the real decrement to the imaginary part of the complex refractive index, the proposed method requires a single image for extraction of the dark-field signal in two-dimensional projection imaging. In the case of three-dimensional tomographic imaging, the method needs only one image to be collected at each projection angle. Initial examples using simulated and experimental data indicate that this method can improve visualization of small sharp features inside a larger object, e.g. the visualization of microcalcifications in propagation-based x-ray breast cancer imaging. It is suggested that the proposed approach may be useful in other forms of biomedical imaging, where it can help one to obtain additional small-angle scattering information without increasing the radiation dose to the sample.

The relative roles of multiple electron scattering and in-molecule Fresnel diffraction in atomic-resolution transmission electron microscopy of small molecules are discussed. It is argued that, while multiple scattering tends to have only a moderate effect, the in-molecule Fresnel diffraction can be significant due to the shallow depth of focus in high-resolution electron microscopy. When the depth of focus is smaller than the size of molecule along the trajectory of the electrons, the projection approximation breaks down and pairs of images at orientations of the molecule rotated by 180 degrees are not the same. Put another way, it matters which way the electrons traverse the specimen. In such cases, diffraction tomography based on the first Born or first Rytov approximation represents a more suitable method for the reconstruction of the three-dimensional distribution of the electrostatic potential, compared to conventional computed tomography which intrinsically relies on the validity of the projection approximation. A simplified method for diffraction tomography based on the approximation of Fresnel diffraction by the transport of intensity equation is proposed and tested on numerically simulated examples., (Copyright © 2020 Elsevier B.V. All rights reserved.)

The deterioration of both the signal-to-noise ratio and the spatial resolution in the electron-density distribution reconstructed from diffraction intensities collected at different orientations of a sample is analysed theoretically with respect to the radiation damage to the sample and the variations in the X-ray intensities illuminating different copies of the sample. The simple analytical expressions and numerical estimates obtained for models of radiation damage and incident X-ray pulses may be helpful in planning X-ray free-electron laser (XFEL) imaging experiments and in analysis of experimental data. This approach to the analysis of partially coherent X-ray imaging configurations can potentially be used for analysis of other forms of imaging where the temporal behaviour of the sample and the incident intensity during exposure may affect the inverse problem of sample reconstruction.

Resolving the electronic structure of single biological molecules in their native state was among the primary motivations behind X-ray free-electron lasers. The ultra-short pulses they produce can outrun the atomic motion induced by radiation damage, but the electronic structure of the sample is still significantly modified from its original state. This paper explores the decoherence of the scattered signal induced by temporal evolution of the electronic structure in the sample molecule. It is shown that the undamaged electron density of a single-molecule sample can often be retrieved using only the two most occupied modes from the coherent mode decomposition of the partially coherent diffraction fluence., (© Alexander Kozlov et al. 2020.)

In recent years, X-ray speckle-tracking techniques have emerged as viable tools for wavefront metrology and sample imaging applications. These methods are based on the measurement of near-field images. Thanks to their simple experimental setup, high angular sensitivity and compatibility with low-coherence sources, these methods have been actively developed for use with synchrotron and laboratory light sources. Not only do speckle-tracking techniques give the potential for high-resolution imaging, but they also provide rapid and robust characterization of aberrations of X-ray optical elements, focal spot profiles, and sample position and transmission properties. In order to realize these capabilities, software implementations are required that are equally rapid and robust. To address this need, a software suite has been developed for the ptychographic X-ray speckle-tracking technique, an X-ray speckle-based method suitable for highly divergent wavefields. The software suite is written in Python 3, with an OpenCL back end for GPU and multi-CPU core processing. It is accessible as a Python module, through the command line or through a graphical user interface, and is available as source code under Version 3 or later of the GNU General Public License., (© Andrew J. Morgan et al. 2020.)

A method is presented for the measurement of the phase gradient of a wavefront by tracking the relative motion of speckles in projection holograms as a sample is scanned across the wavefront. By removing the need to obtain an undistorted reference image of the sample, this method is suitable for the metrology of highly divergent wavefields. Such wavefields allow for large magnification factors that, according to current imaging capabilities, will allow for nanoradian angular sensitivity and nanoscale sample projection imaging. Both the reconstruction algorithm and the imaging geometry are nearly identical to that of ptychography, except that the sample is placed downstream of the beam focus and that no coherent propagation is explicitly accounted for. Like other X-ray speckle tracking methods, it is robust to low-coherence X-ray sources, making it suitable for laboratory-based X-ray sources. Likewise, it is robust to errors in the registered sample positions, making it suitable for X-ray free-electron laser facilities, where beam-pointing fluctuations can be problematic for wavefront metrology. A modified form of the speckle tracking approximation is also presented, based on a second-order local expansion of the Fresnel integral. This result extends the validity of the speckle tracking approximation and may be useful for similar approaches in the field., (© Andrew J. Morgan et al. 2020.)

We show that the width of an arbitrary function and the width of the distribution of its values cannot be made arbitrarily small simultaneously. In the case of ergodic stochastic processes, an ensuing uncertainty relationship is then demonstrated for the product of correlation length and variance. A closely related uncertainty principle is also established for the average degree of fourth-order coherence and the spatial width of modes of bosonic quantum fields. However, it is shown that, in the case of stochastic and quantum observables, certain non-classical states with sub-Poissonian statistics, such as for example photon number squeezed states in quantum optics, can overcome the "classical" noise-resolution uncertainty limit. This uncertainty relationship, which is fundamentally different from the Heisenberg and related uncertainty principles, can define an upper limit for the information capacity of communication and imaging systems. It is expected to be useful in a variety of problems in classical and quantum optics and imaging.

Objectives: To evaluate and compare the image quality of propagation-based phase-contrast computed tomography (PB-CT) using synchrotron radiation and conventional cone-beam breast computed tomography (CBBCT) based on various radiological image quality criteria., Methods: Eight excised breast tissue samples of various sizes and containing different lesion types were scanned using PB-CT at a synchrotron facility and using CBBCT at a university-affiliated breast imaging centre. PB-CT scans were performed at two different mean glandular dose (MGD) levels: standard (5.8 mGy) and low (1.5 mGy), for comparison with CBBCT scans at the standard MGD (5.8 mGy). Image quality assessment was carried out using six quality criteria and six independent medical imaging experts in a reading room with mammography workstations. The interobserver agreement between readers was evaluated using intraclass correlation coefficient (ICC), and image quality was compared between the two breast imaging modalities using the area under the visual grading characteristic curve (AUC VGC )., Results: Interobserver agreement between the readers showed moderate reliability for five image criteria (ICC: ranging from 0.488 to 0.633) and low reliability for one criterion (image noise) (ICC 0.307). For five image quality criteria (overall quality, perceptible contrast, lesion sharpness, normal tissue interfaces, and calcification visibility), both standard-dose PB-CT images (AUC VGC 0.958 to 1, p ≤ .05) and low dose PB-CT images (AUC VGC 0.785 to 0.834, p ≤ .05) were of significantly higher image quality than standard-dose CBBCT images., Conclusions: Synchrotron-based PB-CT can achieve a significantly higher radiological image quality at a substantially lower radiation dose compared with conventional CBBCT., Key Points: • PB-CT using synchrotron radiation results in higher image quality than conventional CBBCT for breast imaging. • PB-CT using synchrotron radiation requires a lower radiation dose than conventional CBBCT for breast imaging. • PB-CT can help clinicians diagnose patients with breast cancer.

In this paper, we present and review the most recent computational advances in the BERTHA code. BERTHA can be regarded as the state of the art in fully relativistic four-component Dirac-Kohn-Sham (DKS) software. Thanks to the implementation of various parallelization and memory open-ended distribution schemes in combination with efficient "density fitting" algorithms, it greatly reduces the computational burden of four-component DKS calculations. We also report the newly developed OpenMP version of the code, that, together with the berthmod Python module, provides a significant leap forward in terms of usability and applicability of the BERTHA software. Some applications of the recently developed natural orbitals for chemical valence/charge displacement bonding analysis and the real-time time dependent DKS implementation are also reported.

We present a real-time time-dependent four-component Dirac-Kohn-Sham (RT-TDDKS) implementation based on the BERTHA code. This new implementation takes advantage of modern software engineering, including the prototyping techniques. The software design follows a three step approach: (i) the prototype implementation of a time-propagation algorithm in nonrelativistic real-time TDDFT within the Psi4NumPy framework, which provides a suitable environment for the creation of a clear, readable, and easy to test reference code in Python, (ii) the design of an original Python application programming interface for the relativistic four-component code BERTHA (PyBERTHA), which has an efficient computational kernel for relativistic integrals written in FORTRAN, and (iii) the porting of the time-propagation scheme enveloped within the Psi4NumPy framework to PyBERTHA. The propagation scheme consequently resides in a single readable Python computer code that is easy to maintain and in which the key quantities, such as the Dirac-Kohn-Sham and dipole matrices, can be accessed directly from the PyBERTHA module. For linear algebra operations (matrix-matrix multiplications and diagonalization) we use the highly optimized procedures implemented in the popular NumPy library. The overhead introduced by the Python interface to BERTHA is almost negligible (less than 1% evaluated on the SCF procedure), and the interoperability between different programming languages (FORTRAN, C, and Python) does not affect the numerical stability of the time-propagation scheme. Our new RT-TDDKS implementation has been employed to investigate the stability of the time-propagation procedure in combination with a density-fitting algorithm (both for the Coulomb and for the exchange-correlation matrix construction), which are employed in BERTHA to speed up the Dirac-Kohn-Sham matrix evaluation. On the basis of systematic calculations, employing several density-fitting basis sets of increasing accuracy, we showed that quantitative agreement can be achieved in combination with extended-fitting basis sets, with an error in the Coulomb energy below 1 μ-hartree. Convergence of the transition energies increasing of quality of the fitting basis sets has been also observed. Our data suggest that the error in the Coulomb energy may also represent a good estimate of the fitting basis set quality for real-time electron dynamics simulations. Further, we study the applicability of the RT-TDDKS method in combination with both weak- and extreme strong-field regime. Numerical results of excited-state transitions for the Group 12 atoms are reported and compared with a previous real-time Dirac-Kohn-Sham implementation (Repisky et al. J. Chem. Theory Comput. 2015, 11, 980-991). Finally, calculations of high harmonic generation in the hydrogen molecule and Au dimer have been also carried out. We were able to generate high harmonics with relatively well-defined peaks up to the 21st and 13th order in the case of H 2 and Au 2 , respectively. Our findings show that the four-component structure of the Dirac-Kohn-Sham Hamiltonian provides a suitable theoretical framework, with no intrinsic unfavorable features, to study molecules in the strong-field regime.

X-ray free-electron lasers (XFELs) enable crystallographic structure determination beyond the limitations imposed upon synchrotron measurements by radiation damage. The need for very short XFEL pulses is relieved through gating of Bragg diffraction by loss of crystalline order as damage progresses, but not if ionization events are spatially non-uniform due to underlying elemental distributions, as in biological samples. Indeed, correlated movements of iron and sulfur ions were observed in XFEL-irradiated ferredoxin microcrystals using unusually long pulses of 80 fs. Here, we report a femtosecond time-resolved X-ray pump/X-ray probe experiment on protein nanocrystals. We observe changes in the protein backbone and aromatic residues as well as disulfide bridges. Simulations show that the latter's correlated structural dynamics are much slower than expected for the predicted high atomic charge states due to significant impact of ion caging and plasma electron screening. This indicates that dense-environment effects can strongly affect local radiation damage-induced structural dynamics.

Purpose: Propagation-based phase-contrast computed tomography (PB-CT) is a method for three-dimensional x-ray imaging that utilizes refraction, as well as absorption, of x rays in the tissues to increase the signal-to-noise ratio (SNR) in the resultant images, in comparison with equivalent conventional absorption-only x-ray tomography (CT). Importantly, the higher SNR is achieved without sacrificing spatial resolution or increasing the radiation dose delivered to the imaged tissues. The present work has been carried out in the context of the current development of a breast CT imaging facility at the Australian Synchrotron., Methods: Seven unfixed complete mastectomy samples with and without breast cancer lesions have been imaged using absorption-only CT and PB-CT techniques under controlled experimental conditions. The radiation doses delivered to the mastectomy samples during the scans were comparable to those approved for mammographic screening. Physical characteristics of the reconstructed images, such as spatial resolution and SNR, have been measured and compared with the results of the radiological quality assessment of the complete absorption CT and PB-CT image stacks., Results: Despite the presence of some image artefacts, the PB-CT images have outperformed comparable absorption CT images collected at the same radiation dose, in terms of both the measured objective image characteristics and the radiological image scores. The outcomes of these experiments are shown to be consistent with predictions of the theory of PB-CT imaging and previous reported experimental studies of this imaging modality., Conclusions: The results presented in this paper demonstrate that PB-CT holds a high potential for improving on the quality and diagnostic value of images obtained using existing medical x-ray technologies, such as mammography and digital breast tomosynthesis (DBT). If implemented at suitable synchrotron imaging facilities, PB-CT can be used to complement existing imaging modalities, leading to more accurate breast cancer diagnosis., (© 2019 American Association of Physicists in Medicine.)

It is shown that the average signal-to-noise ratio (SNR) in the three-dimensional electron-density distribution of a sample reconstructed by coherent diffractive imaging cannot exceed twice the square root of the ratio of the mean total number of scattered photons detected during the scan and the number of spatially resolved voxels in the reconstructed volume. This result leads to an upper bound on Shannon's information capacity of this imaging method by specifying the maximum number of distinguishable density distributions within the reconstructed volume when the radiation dose delivered to the sample and the spatial resolution are both fixed. If the spatially averaged SNR in the reconstructed electron density is fixed instead, the radiation dose is shown to be proportional to the third or fourth power of the spatial resolution, depending on the sampling of the three-dimensional diffraction space and the scattering power of the sample.

We have recently introduced a simple yet powerful tool for analyzing quantitatively the coordination bond in terms of the donation and back-donation constituents of the Dewar-Chatt-Duncanson model. Our approach is based on the decomposition, via natural orbitals for chemical valence (NOCV), of the so-called charge-displacement (CD) function into additive chemically meaningful components (Bistoni et al. J. Chem. Phys. 2015, 142, 084112 ). The method, referred to as NOCV/CD, provides clear-cut measures of donation and back-donation charge flows following bond formation, and its robustness has been demonstrated by a tight correlation of the related charge-transfer estimates with experimental observables. In this paper we extend the NOCV/CD analysis scheme to the four-component relativistic framework, which includes spin-orbit coupling variationally. This formalism is incorporated into a recently developed, highly efficient parallel version of the relativistic Dirac-Kohn-Sham (DKS) program BERTHA (Rampino et al. J. Chem. Theory Comput. 2014, 10, 3766-3776 ). We test the accuracy and numerical stability of this new implementation through the analysis of the convergence properties of the basis sets employed to expand the DKS spinor solution and those used to linearize the electronic density in the density-fitting algorithm speeding up the evaluation of the DKS matrix. An illustration of NOCV/CD analysis in the relativistic framework is also given through a study of the metal-carbonyl coordination bond in a series of [M-CO] + (M = Cu, Ag, Au) complexes of group 11 metals, where relativistic effects, including spin-orbit coupling, are found to play an important role.

The effectiveness of reconstructive imaging using the homogeneous transport of intensity equation may be regarded as "unreasonable," because it has been shown to significantly increase signal-to-noise ratio while preserving spatial resolution, compared to equivalent conventional absorption-based imaging techniques at the same photon fluence. We reconcile this surprising behavior by analyzing the propagation of noise in typical in-line holography experiments. This analysis indicates that novel imaging techniques may be designed that produce high signal-to-noise images at low radiation doses without sacrificing spatial resolution.

The recent availability of extremely intense, femtosecond X-ray free-electron laser (XFEL) sources has spurred the development of serial femtosecond nanocrystallography (SFX). Here, SFX is used to analyze nanoscale crystals of β-hematin, the synthetic form of hemozoin which is a waste by-product of the malaria parasite. This analysis reveals significant differences in β-hematin data collected during SFX and synchrotron crystallography experiments. To interpret these differences two possibilities are considered: structural differences between the nanocrystal and larger crystalline forms of β-hematin, and radiation damage. Simulation studies show that structural inhomogeneity appears at present to provide a better fit to the experimental data. If confirmed, these observations will have implications for designing compounds that inhibit hemozoin formation and suggest that, for some systems at least, additional information may be gained by comparing structures obtained from nanocrystals and macroscopic crystals of the same molecule.

A reciprocal relationship between the autocovariance of the light intensity in the source plane and in the far-field detector plane is presented in a form analogous to the classical van Cittert-Zernike theorem, but involving intensity correlation functions. A "classical" version of the reciprocity relationship is considered first, based on the assumption of circular Gaussian statistics of the complex amplitudes in the source plane. The result is consistent with the theory of Hanbury Brown-Twiss interferometry, but it is shown to be also applicable to estimation of the source size or the spatial resolution of the detector from the noise power spectrum of flat-field images. An alternative version of the van Cittert-Zernike theorem for intensity correlations is then derived for a quantized electromagnetic beam in a coherent state, which leads to Poisson statistics for the intrinsic intensity of the beam.

The precise details of the interaction of intense X-ray pulses with matter are a topic of intense interest to researchers attempting to interpret the results of femtosecond X-ray free electron laser (XFEL) experiments. An increasing number of experimental observations have shown that although nuclear motion can be negligible, given a short enough incident pulse duration, electronic motion cannot be ignored. The current and widely accepted models assume that although electrons undergo dynamics driven by interaction with the pulse, their motion could largely be considered 'random'. This would then allow the supposedly incoherent contribution from the electronic motion to be treated as a continuous background signal and thus ignored. The original aim of our experiment was to precisely measure the change in intensity of individual Bragg peaks, due to X-ray induced electronic damage in a model system, crystalline C60. Contrary to this expectation, we observed that at the highest X-ray intensities, the electron dynamics in C60 were in fact highly correlated, and over sufficiently long distances that the positions of the Bragg reflections are significantly altered. This paper describes in detail the methods and protocols used for these experiments, which were conducted both at the Linac Coherent Light Source (LCLS) and the Australian Synchrotron (AS) as well as the crystallographic approaches used to analyse the data.

An uncertainty inequality is presented that establishes a lower limit for the product of the variance of the time-averaged intensity of a mode of a quantized electromagnetic field and the degree of its spatial localization. The lower limit is determined by the vacuum fluctuations within the volume corresponding to the width of the mode. This result also leads to a generalized form of the Heisenberg uncertainty principle for boson fields in which the lower limit for the product of uncertainties in the spatial and momentum localization of a mode is equal to the product of Planck's constant and a dimensionless functional which reflects the joint signal-to-noise ratio of the position and momentum of vacuum fluctuations in the region of the phase space occupied by the mode. Experimental X-ray synchrotron measurements provide an initial verification of the proposed theory in the case of Poisson statistics.

X-ray free-electron lasers (XFELs) deliver x-ray pulses with a coherent flux that is approximately eight orders of magnitude greater than that available from a modern third-generation synchrotron source. The power density of an XFEL pulse may be so high that it can modify the electronic properties of a sample on a femtosecond time scale. Exploration of the interaction of intense coherent x-ray pulses and matter is both of intrinsic scientific interest and of critical importance to the interpretation of experiments that probe the structures of materials using high-brightness femtosecond XFEL pulses. We report observations of the diffraction of extremely intense 32-fs nanofocused x-ray pulses by a powder sample of crystalline C60. We find that the diffraction pattern at the highest available incident power significantly differs from the one obtained using either third-generation synchrotron sources or XFEL sources operating at low output power and does not correspond to the diffraction pattern expected from any known phase of crystalline C60. We interpret these data as evidence of a long-range, coherent dynamic electronic distortion that is driven by the interaction of the periodic array of C60 molecular targets with intense x-ray pulses of femtosecond duration.

Ionization of atoms and molecules in strong laser fields is a fundamental process in many fields of research, especially in the emerging field of attosecond science. So far, demonstrably accurate data have only been acquired for atomic hydrogen (H), a species that is accessible to few investigators. Here, we present measurements of the ionization yield for argon, krypton, and xenon with percent-level accuracy, calibrated using H, in a laser regime widely used in attosecond science. We derive a transferable calibration standard for laser peak intensity, accurate to 1.3%, that is based on a simple reference curve. In addition, our measurements provide a much needed benchmark for testing models of ionization in noble-gas atoms, such as the widely employed single-active electron approximation.

Serial femtosecond X-ray crystallography (SFX) has created new opportunities in the field of structural analysis of protein nanocrystals. The intensity and timescale characteristics of the X-ray free-electron laser sources used in SFX experiments necessitate the analysis of a large collection of individual crystals of variable shape and quality to ultimately solve a single, average crystal structure. Ensembles of crystals are commonly encountered in powder diffraction, but serial crystallography is different because each crystal is measured individually and can be oriented via indexing and merged into a three-dimensional data set, as is done for conventional crystallography data. In this way, serial femtosecond crystallography data lie in between conventional crystallography data and powder diffraction data, sharing features of both. The extremely small sizes of nanocrystals, as well as the possible imperfections of their crystallite structure, significantly affect the diffraction pattern and raise the question of how best to extract accurate structure-factor moduli from serial crystallography data. Here it is demonstrated that whole-pattern fitting techniques established for one-dimensional powder diffraction analysis can be feasibly extended to higher dimensions for the analysis of merged SFX diffraction data. It is shown that for very small crystals, whole-pattern fitting methods are more accurate than Monte Carlo integration methods that are currently used.

During the last five years, serial femtosecond crystallography using X-ray laser pulses has been developed into a powerful technique for determining the atomic structures of protein molecules from micrometre- and sub-micrometre-sized crystals. One of the key reasons for this success is the 'self-gating' pulse effect, whereby the X-ray laser pulses do not need to outrun all radiation damage processes. Instead, X-ray-induced damage terminates the Bragg diffraction prior to the pulse completing its passage through the sample, as if the Bragg diffraction were generated by a shorter pulse of equal intensity. As a result, serial femtosecond crystallography does not need to be performed with pulses as short as 5-10 fs, but can succeed for pulses 50-100 fs in duration. It is shown here that a similar gating effect applies to single-molecule diffraction with respect to spatially uncorrelated damage processes like ionization and ion diffusion. The effect is clearly seen in calculations of the diffraction contrast, by calculating the diffraction of the average structure separately to the diffraction from statistical fluctuations of the structure due to damage ('damage noise'). The results suggest that sub-nanometre single-molecule imaging with 30-50 fs pulses, like those produced at currently operating facilities, should not yet be ruled out. The theory presented opens up new experimental avenues to measure the impact of damage on single-particle diffraction, which is needed to test damage models and to identify optimal imaging conditions.

We propose a new complete memory-distributed algorithm, which significantly improves the parallel implementation of the all-electron four-component Dirac-Kohn-Sham (DKS) module of BERTHA (J. Chem. Theory Comput. 2010, 6, 384). We devised an original procedure for mapping the DKS matrix between an efficient integral-driven distribution, guided by the structure of specific G-spinor basis sets and by density fitting algorithms, and the two-dimensional block-cyclic distribution scheme required by the ScaLAPACK library employed for the linear algebra operations. This implementation, because of the efficiency in the memory distribution, represents a leap forward in the applicability of the DKS procedure to arbitrarily large molecular systems and its porting on last-generation massively parallel systems. The performance of the code is illustrated by some test calculations on several gold clusters of increasing size. The DKS self-consistent procedure has been explicitly converged for two representative clusters, namely Au20 and Au34, for which the density of electronic states is reported and discussed. The largest gold cluster uses more than 39k basis functions and DKS matrices of the order of 23 GB.

Characterization of microscopic structural order and in particular medium range order (MRO) in amorphous materials is challenging. A new technique is demonstrated that allows analysis of MRO using X-rays. Diffraction data were collected from a sample consisting of densely packed polystyrene-latex micro-spheres. Ptychography is used to reconstruct the sample transmission function and fluctuation microscopy applied to characterize structural order producing a detailed `fluctuation map' allowing analysis of the sample at two distinct length scales. Independent verification is provided via X-ray diffractometry. Simulations of dense random packing of spheres have also been used to explore the origin of the structural order measured.

We present a gas jet array for use in high-order harmonic generation experiments. Precise control of the pressure in each individual gas jet has allowed a thorough investigation into mechanisms contributing to the selective enhancement observed in the harmonic spectra produced by dual-gas, multi-jet arrays. Our results reveal that in our case, the dominant enhancement mechanism is the result of a compression of the harmonic-producing gas jet due to the presence of other gas jets in the array. The individual control of the gas jets in the array also provides a promising method for enhancing the harmonic yield by precise tailoring of the length and pressure gradient of the interaction region.

The recent development of X-ray free-electron laser sources has created new opportunities for the structural analysis of protein nanocrystals. The extremely small sizes of the crystals, as well as imperfections of the crystal structure, result in an interference phenomenon in the diffraction pattern. With decreasing crystallite size the structural imperfections play a role in the formation of the diffraction pattern that is comparable in importance to the size effects and should be taken into account during the data analysis and structure reconstruction processes. There now exists a need to develop new methods of protein structure determination that do not depend on the availability of good-quality crystals and that can treat proteins under conditions close to the active form. This paper demonstrates an approach that is specifically tailored to nanocrystalline samples and offers a unique crystallographic solution.

We present a new interferometer technique whereby multiple extreme ultraviolet light pulses are generated at different positions within a single laser focus (i.e., from successive sources) with a highly controllable time delay. The interferometer technique is tested with two generating media to create two extreme ultraviolet light pulses with a time delay between them. The delay is found to be a consequence of the Gouy phase shift. Ultimately the apparatus is capable of accessing unprecedented time scales by allowing stable and repeatable delays as small as 100 zs.

We demonstrate Fresnel Coherent Diffractive Imaging (FCDI) tomography in the X-ray regime. The method uses an incident X-ray illumination with known curvature in combination with ptychography to overcome existing problems in diffraction imaging. The resulting tomographic reconstruction represents a 3D map of the specimen's complex refractive index at nano-scale resolution. We use this technique to image a lithographically fabricated glass capillary, in which features down to 70nm are clearly resolved.

Observations of long-lived coherences in photosynthetic light-harvesting complexes utilize short pulses with broad spectral bandwidths to coherently excite multiple transitions and coherent superpositions. In order to identify the role that such quantum effects might play in efficient energy transfer, however, an alternative approach is required. We have developed a technique for two-color photon echo spectroscopy to selectively excite the pathway of interest and measure its evolution in the absence of any other excitation. We use this technique to excite a coherence pathway in phycocyanin-645 from cryptophyte algae and measure the dynamics of this coherence. A decoherence time of 500 fs was measured, and clear signatures for strong coupling between the electronic states and phonon modes were observed, allowing coherent coupling between otherwise nonresonant transitions. This provides detailed experimental evidence of the long-lived coherences and the nature of the quantum mechanical interactions between electronic states and phonon modes in phycocyanin-645 from cryptophyte marine algae.

We present experimental data on strong-field ionization of atomic hydrogen by few-cycle laser pulses. We obtain quantitative agreement at the 10% level between the data and an ab initio simulation over a wide range of laser intensities and electron energies.

We demonstrate three-dimensional (3D) electronic spectroscopy of excitons in a double quantum well system using a three-dimensional phase retrieval algorithm to obtain the phase information that is lost in the measurement of intensities. By extending the analysis of two-dimensional spectroscopy to three dimensions, contributions from different quantum mechanical pathways can be further separated allowing greater insight into the mechanisms responsible for the observed peaks. By examining different slices of the complete three-dimensional spectrum, not only can the relative amplitudes be determined, but the peak shapes can also be analysed to reveal further details of the interactions with the environment and inhomogeneous broadening. We apply this technique to study the coupling between two coupled quantum wells, 5.7 nm and 8 nm wide, separated by a 4 nm barrier. Coupling between the heavy-hole excitons of each well results in a circular cross-peak indicating no correlation of the inhomogeneous broadening. An additional cross-peak is isolated in the 3D spectrum which is elongated in the diagonal direction indicating correlated inhomogeneous broadening. This is attributed to coupling of the excitons involving the two delocalised light-hole states and the electron state localised on the wide well. The attribution of this peak and the analysis of the peak shapes is supported by numerical simulations of the electron and hole wavefunctions and the three-dimensional spectrum based on a density matrix approach. An additional benefit of extending the phase retrieval algorithm from two to three dimensions is that it becomes substantially more reliable and less susceptible to noise as a result of the more extensive use of a priori information., (© 2011 American Institute of Physics)

We review recent theoretical and computational advances in the full relativistic four-component Dirac-Kohn-Sham (DKS) approach and its application to the calculation of the electronic structure of chemical systems containing many heavy atoms. We describe our implementation of an all-electron DKS approach based on the use of G-spinor basis sets, Hermite Gaussian functions, state-of-the-art density-fitting techniques and memory distributed parallelism. This approach has enormously extended the applicability of the DKS method, including for example large clusters of heavy atoms, and opens the way for future key developments. We examine the current limitations and future possible applications of the DKS approach, including the implementation of four-current density functionals and real-time propagation schemes. This would make possible to describe molecules in strong fields, accurately accounting for relativistic kinematic effects and spin-orbit coupling., (This journal is © the Owner Societies 2011)

The application of powder diffraction methods to problems in structural biology is generally regarded as intractable because of the large number of unresolved, overlapping X-ray reflections. Here, we use information about unit cell lattice parameters, space group transformations, and chemical composition as a priori information in a bootstrap process that resolves the ambiguities associated with overlapping reflections. The measured ratios of reflections that can be resolved experimentally are used to refine the position, the shape, and the orientation of low-resolution molecular structures within the unit cell, in leading to the resolution of the overlapping reflections. The molecular model is then made progressively more sophisticated as additional diffraction information is included in the analysis. We apply our method to the recovery of the structure of the bacteriorhodopsin molecule (bR) to a resolution of 7 Å using experimental data obtained from two-dimensional purple membrane crystals. The approach can be used to determine the structure factors directly or to provide reliable low-resolution phase information that can be refined further by the conventional methods of protein crystallography., (Copyright © 2011 The Protein Society.)

We present here the details of a phase retrieval technique that provides access to multidimensional modalities that are not currently available using existing interferometric techniques. The development of multidimensional optical spectroscopy has facilitated significant insights into electronic processes in physics, chemistry, and biology. The versatility and number of available techniques are, however, significantly limited by the requirement that the detection be interferometric. Many of these techniques are closely related to the vast range of multidimensional NMR spectroscopies, which revolutionized analytical chemistry more than 30 years ago. We focus here on the specific case of two-color multidimensional spectroscopy (analogous to heteronuclear NMR) and discuss the details of an iterative algorithm that recovers the relative phase relationships required to perform the Fourier transformation and find the unique solution for the 2D spectrum. A detailed guide is provided that describes the practical implementation of such algorithms. The effectiveness and accuracy of the phase retrieval process are assessed for simulated one- and two-color experiments. It is also compared with one-color experimental data for which the target phase information has been obtained independently by interferometry. In all the cases, the present algorithm yields results that compare well with the solutions obtained by other means. There are, however, some limitations and potential pitfalls that are identified and discussed. We conclude with a discussion of the potential applications and further advances that may be possible by adopting iterative phase retrieval algorithms of the type discussed here.

Tomographic coherent imaging requires the reconstruction of a series of two-dimensional projections of the object. We show that using the solution for the image of one projection as the starting point for the reconstruction of the next projection offers a reliable and rapid approach to the image reconstruction. The method is demonstrated on simulated and experimental data. This technique also simplifies reconstructions using data with curved incident wavefronts.

We show that all-electron relativistic four-component Dirac-Kohn-Sham (DKS) computations, using G-spinor basis sets and state-of-the-art density fitting algorithms, can be efficiently parallelized and applied to large molecular systems, including large clusters of heavy atoms. The performance of the parallel implementation of the DKS module of the program BERTHA is illustrated and analyzed by some test calculations on several gold clusters up to Au32, showing that calculations with more than 25 000 basis functions (i.e., DKS matrices on the order of 10 GB) are now feasible. As a first application of this novel implementation, we investigate the interaction of the atom Hg with the Au20 cluster.

We demonstrate a new imaging technique for cold atom clouds based on phase retrieval from a single diffraction measurement. Most single-shot diffractive imaging methods for cold atoms assume a monomorphic object to extract the column density. The method described here allows quantitative imaging of an inhomogeneous cloud, enabling recovery of either the atomic density or the refractive index, provided the other is known. Using ideas borrowed from density functional theory, we calculate the approximate paraxial diffracted intensity derivative from the measured diffracted intensity distribution and use it to solve the Transport of Intensity Equation (TIE) for the phase of the wave at the detector plane. Back-propagation to the object plane yields the object exit surface wave and then provides a quantitative measurement of either the atomic column density or refractive index. Images of homogeneous clouds showed good quantitative agreement with conventional techniques. An inhomogeneous cloud was created using a cascade electromagnetically induced transparency scheme and images of both phase and amplitude parts of refractive index across the cloud were separately retrieved, showing good agreement with theoretical results.

The measured spatial coherence characteristics of the illumination used in a diffractive imaging experiment are incorporated in an algorithm that reconstructs the complex transmission function of an object from experimental x-ray diffraction data using 1.4 keV x rays. Conventional coherent diffractive imaging, which assumes full spatial coherence, is a limiting case of our approach. Even in cases in which the deviation from full spatial coherence is small, we demonstrate a significant improvement in the quality of wave field reconstructions. Our formulation is applicable to x-ray and electron diffraction imaging techniques provided that the spatial coherence properties of the illumination are known or can be measured.

A method for numerically recovering the coherent modes and their occupancies from a known mutual optical intensity function is described. As an example, the technique is applied to previously published experimental data from an x-ray undulator source. The data are found to be described by three coherent modes, and the functional forms and relative occupancies of these modes are recovered.