Femtosecond Reduction of Atomic Scattering Factors Triggered by Intense X-Ray Pulse.

X-ray diffraction of silicon irradiated with tightly focused femtosecond x-ray pulses (photon energy, 11.5 keV; pulse duration, 6 fs) was measured at various x-ray intensities up to 4.6×10^{19} W/cm^{2}. The measurement reveals that the diffraction intensity is highly suppressed when the x-ray intensity reaches of the order of 10^{19} W/cm^{2}. With a dedicated simulation, we confirm that the observed reduction of the diffraction intensity can be attributed to the femtosecond change in individual atomic scattering factors due to the ultrafast creation of highly ionized atoms through photoionization, Auger decay, and subsequent collisional ionization. We anticipate that this ultrafast reduction of atomic scattering factor will be a basis for new x-ray nonlinear techniques, such as pulse shortening and contrast variation x-ray scattering.

Tree-Code Based Improvement of Computational Performance of the X-ray-Matter-Interaction Simulation Tool XMDYN.

In this work, we report on incorporating for the first time tree-algorithm based solvers into the molecular dynamics code, XMDYN. XMDYN was developed to describe the interaction of ultrafast X-ray pulses with atomic assemblies. It is also a part of the simulation platform, SIMEX, developed for computational single-particle imaging studies at the SPB/SFX instrument of the European XFEL facility. In order to improve the XMDYN performance, we incorporated the existing tree-algorithm based Coulomb solver, PEPC, into the code, and developed a dedicated tree-algorithm based secondary ionization solver, now also included in the XMDYN code. These extensions enable computationally efficient simulations of X-ray irradiated large atomic assemblies, e.g., large protein systems or viruses that are of strong interest for ultrafast X-ray science. The XMDYN-based preparatory simulations can now guide future single-particle-imaging experiments at the free-electron-laser facility, EuXFEL.

Atomic-Scale Visualization of Ultrafast Bond Breaking in X-Ray-Excited Diamond.

Ultrafast changes of charge density distribution in diamond after irradiation with an intense x-ray pulse (photon energy, 7.8 keV; pulse duration, 6 fs; intensity, 3×10^{19} W/cm^{2}) have been visualized with the x-ray pump-x-ray probe technique. The measurement reveals that covalent bonds in diamond are broken and the electron distribution around each atom becomes almost isotropic within ∼5 fs after the intensity maximum of the x-ray pump pulse. The 15 fs time delay observed between the bond breaking and atomic disordering indicates nonisothermality of electron and lattice subsystems on this timescale. From these observations and simulation results, we interpret that the x-ray-induced change of the interatomic potential drives the ultrafast atomic disordering underway to the following nonthermal melting.

xcalib: a focal spot calibrator for intense X-ray free-electron laser pulses based on the charge state distributions of light atoms.

The xcalib toolkit has been developed to calibrate the beam profile of an X-ray free-electron laser (XFEL) at the focal spot based on the experimental charge state distributions (CSDs) of light atoms. Characterization of the fluence distribution at the focal spot is essential to perform the volume integrations of physical quantities for a quantitative comparison between theoretical and experimental results, especially for fluence-dependent quantities. The use of the CSDs of light atoms is advantageous because CSDs directly reflect experimental conditions at the focal spot, and the properties of light atoms have been well established in both theory and experiment. Theoretical CSDs are obtained using xatom, a toolkit to calculate atomic electronic structure and to simulate ionization dynamics of atoms exposed to intense XFEL pulses, which involves highly excited multiple core-hole states. Employing a simple function with a few parameters, the spatial profile of an XFEL beam is determined by minimizing the difference between theoretical and experimental results. The optimization procedure employing the reinforcement learning technique can automatize and organize calibration procedures which, before, had been performed manually. xcalib has high flexibility, simultaneously combining different optimization methods, sets of charge states, and a wide range of parameter space. Hence, in combination with xatom, xcalib serves as a comprehensive tool to calibrate the fluence profile of a tightly focused XFEL beam in the interaction region.

Radiation-Induced Chemical Dynamics in Ar Clusters Exposed to Strong X-Ray Pulses.

We show that electron and ion spectroscopy reveals the details of the oligomer formation in Ar clusters exposed to an x-ray free electron laser (XFEL) pulse, i.e., chemical dynamics triggered by x rays. With guidance from a dedicated molecular dynamics simulation tool, we find that van der Waals bonding, the oligomer formation mechanism, and charge transfer among the cluster constituents significantly affect ionization dynamics induced by an XFEL pulse of moderate fluence. Our results clearly demonstrate that XFEL pulses can be used not only to "damage and destroy" molecular assemblies but also to modify and transform their molecular structure. The accuracy of the predictions obtained makes it possible to apply the cluster spectroscopy, in connection with the respective simulations, for estimation of the XFEL pulse fluence in the fluence regime below single-atom multiple-photon absorption, which is hardly accessible with other diagnostic tools.

Computational study of diffraction image formation from XFEL irradiated single ribosome molecule.

Single particle imaging at atomic resolution is perhaps one of the most desired goals for ultrafast X-ray science with X-ray free-electron lasers. Such a capability would create great opportunity within the biological sciences, as high-resolution structural information of biosamples that may not crystallize is essential for many research areas therein. In this paper, we report on a comprehensive computational study of diffraction image formation during single particle imaging of a macromolecule, containing over one hundred thousand non-hydrogen atoms. For this study, we use a dedicated simulation framework, SIMEX, available at the European XFEL facility. Our results demonstrate the full feasibility of computational single-particle imaging studies for biological samples of realistic size. This finding is important as it shows that the SIMEX platform can be used for simulations to inform relevant single-particle-imaging experiments and help to establish optimal parameters for these experiments. This will enable more focused and more efficient single-particle-imaging experiments at XFEL facilities, making the best use of the resource-intensive XFEL operation.

Water layer and radiation damage effects on the orientation recovery of proteins in single-particle imaging at an X-ray free-electron laser.

The noise caused by sample heterogeneity (including sample solvent) has been identified as one of the determinant factors for a successful X-ray single-particle imaging experiment. It influences both the radiation damage process that occurs during illumination as well as the scattering patterns captured by the detector. Here, we investigate the impact of water layer thickness and radiation damage on orientation recovery from diffraction patterns of the nitrogenase iron protein. Orientation recovery is a critical step for single-particle imaging. It enables to sort a set of diffraction patterns scattered by identical particles placed at unknown orientations and assemble them into a 3D reciprocal space volume. The recovery quality is characterized by a "disconcurrence" metric. Our results show that while a water layer mitigates protein damage, the noise generated by the scattering from it can introduce challenges for orientation recovery and is anticipated to cause problems in the phase retrieval process to extract the desired protein structure. Compared to these disadvantageous effects due to the thick water layer, the effects of radiation damage on the orientation recovery are relatively small. Therefore, minimizing the amount of residual sample solvent should be considered a crucial step in improving the fidelity and resolution of X-ray single-particle imaging experiments.

Plasma environmental effects in the atomic structure for simulating x-ray free-electron-laser-heated solid-density matter.

High energy density (HED) matter exists extensively in the Universe, and it can be created with extreme conditions in laboratory facilities such as x-ray free-electron lasers (XFEL). In HED matter, the electronic structure of individual atomic ions is influenced by a dense plasma environment, and one of the most significant phenomena is the ionization potential depression (IPD). Incorporation of the IPD effects is of great importance in accurate modeling of dense plasmas. All theoretical treatments of IPD so far have been based on the assumption of local thermodynamic equilibrium, but its validity is questionable in ultrafast formation dynamics of dense plasmas, particularly when interacting with intense XFEL pulses. A treatment of transient IPD, based on an electronic-structure calculation of an atom in the presence of a plasma environment described by classical particles, has recently been proposed [Phys. Rev. E 103, 023203 (2021)2470-004510.1103/PhysRevE.103.023203], but its application to and impact on plasma dynamics simulations have not been investigated yet. In this work, we extend XMDYN, a hybrid quantum-classical approach combining Monte Carlo and molecular dynamics, by incorporating the proposed IPD treatment into plasma dynamics simulations. We demonstrate the importance of the IPD effects in theoretical modeling of aluminum dense plasmas by comparing two XMDYN simulations: one with electronic-structure calculations of isolated atoms (without IPD) and the other with those of atoms embedded in a plasma (with IPD). At equilibrium, the mean charge obtained in the plasma simulation with IPD is in good agreement with the full quantum-mechanical average-atom model. The present approach promises to be a reliable tool to simulate the creation and nonequilibrium evolution of dense plasmas induced by ultraintense and ultrashort XFEL pulses.

Chemical effects on the dynamics of organic molecules irradiated with high intensity x rays.

The interaction of a high intensity x-ray pulse with matter causes ionization of the constituent atoms through various atomic processes, and the system eventually goes through a complex structural dynamics. Understanding this whole process is important from the perspective of structure determination of molecules using single particle imaging. XMDYN, which is a classical molecular dynamics-Monte Carlo based hybrid approach, has been successful in simulating the dynamical evolution of various systems under intense irradiation over the past years. The present study aims for extending the XMDYN toolkit to treat chemical bonds using the reactive force field. In order to study its impact, a highly intense x-ray pulse was made to interact with the simplest amino acid, glycine. Different model variants were used to highlight the consequences of charge rearrangement and chemical bonds on the time evolution. The charge-rearrangement-enhanced x-ray ionization of molecules effect is also discussed to address the capability of a classical MD based approach, i.e., XMDYN, to capture such a molecular phenomenon.

Erratum: Transient ionization potential depression in nonthermal dense plasmas at high x-ray intensity [Phys. Rev. E 103, 023203 (2021)].

This corrects the article DOI: 10.1103/PhysRevE.103.023203.

Erratum: Molecular-dynamics approach for studying the nonequilibrium behavior of x-ray-heated solid-density matter [Phys. Rev. E 96, 023205 (2017)].

This corrects the article DOI: 10.1103/PhysRevE.96.023205.

Transient ionization potential depression in nonthermal dense plasmas at high x-ray intensity.

The advent of x-ray free-electron lasers (XFELs), which provide intense ultrashort x-ray pulses, has brought a new way of creating and analyzing hot and warm dense plasmas in the laboratory. Because of the ultrashort pulse duration, the XFEL-produced plasma will be out of equilibrium at the beginning, and even the electronic subsystem may not reach thermal equilibrium while interacting with a femtosecond timescale pulse. In the dense plasma, the ionization potential depression (IPD) induced by the plasma environment plays a crucial role for understanding and modeling microscopic dynamical processes. However, all theoretical approaches for IPD have been based on local thermal equilibrium (LTE), and it has been controversial to use LTE IPD models for the nonthermal situation. In this work, we propose a non-LTE (NLTE) approach to calculate the IPD effect by combining a quantum-mechanical electronic-structure calculation and a classical molecular dynamics simulation. This hybrid approach enables us to investigate the time evolution of ionization potentials and IPDs during and after the interaction with XFEL pulses, without the limitation of the LTE assumption. In our NLTE approach, the transient IPD values are presented as distributions evolving with time, which cannot be captured by conventional LTE-based models. The time-integrated ionization potential values are in good agreement with benchmark experimental data on solid-density aluminum plasma and other theoretical predictions based on LTE. The present work is promising to provide critical insights into nonequilibrium dynamics of dense plasma formation and thermalization induced by XFEL pulses.

Effects of radiation damage and inelastic scattering on single-particle imaging of hydrated proteins with an X-ray Free-Electron Laser.

We present a computational case study of X-ray single-particle imaging of hydrated proteins on an example of 2-Nitrogenase-Iron protein covered with water layers of various thickness, using a start-to-end simulation platform and experimental parameters of the SPB/SFX instrument at the European X-ray Free-Electron Laser facility. The simulations identify an optimal thickness of the water layer at which the effective resolution for imaging the hydrated sample becomes significantly higher than for the non-hydrated sample. This effect is lost when the water layer becomes too thick. Even though the detailed results presented pertain to the specific sample studied, the trends which we identify should also hold in a general case. We expect these findings will guide future single-particle imaging experiments using hydrated proteins.

Structural dynamics in proteins induced by and probed with X-ray free-electron laser pulses.

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.

Towards the theoretical limitations of X-ray nanocrystallography at high intensity: the validity of the effective-form-factor description.

X-ray free-electron lasers (XFELs) broaden horizons in X-ray crystallography. Facilitated by the unprecedented high intensity and ultrashort duration of the XFEL pulses, they enable us to investigate the structure and dynamics of macromolecules with nano-sized crystals. A limitation is the extent of radiation damage in the nanocrystal target. A large degree of ionization initiated by the incident high-intensity XFEL pulse alters the scattering properties of the atoms leading to perturbed measured patterns. In this article, the effective-form-factor approximation applied to capture this phenomenon is discussed. Additionally, the importance of temporal configurational fluctuations at high intensities, shaping these quantities besides the average electron loss, is shown. An analysis regarding the applicability of the approach to targets consisting of several atomic species is made, both theoretically and via realistic radiation-damage simulations. It is concluded that, up to intensities relevant for XFEL-based nanocrystallography, the effective-form-factor description is sufficiently accurate. This work justifies treating measured scattering patterns using conventional structure-reconstruction algorithms.

Relativistic and resonant effects in the ionization of heavy atoms by ultra-intense hard X-rays.

An accurate description of the interaction of intense hard X-ray pulses with heavy atoms, which is crucial for many applications of free-electron lasers, represents a hitherto unresolved challenge for theory because of the enormous number of electronic configurations and relativistic effects, which need to be taken into account. Here we report results on multiple ionization of xenon atoms by ultra-intense (about 1019 W/cm2) femtosecond X-ray pulses at photon energies from 5.5 to 8.3 keV and present a theoretical model capable of reproducing the experimental data in the entire energy range. Our analysis shows that the interplay of resonant and relativistic effects results in strongly structured charge state distributions, which reflect resonant positions of relativistically shifted electronic levels of highly charged ions created during the X-ray pulse. The theoretical approach described here provides a basis for accurate modeling of radiation damage in hard X-ray imaging experiments on targets with high-Z constituents.

Molecular-dynamics approach for studying the nonequilibrium behavior of x-ray-heated solid-density matter.

When matter is exposed to a high-intensity x-ray free-electron-laser pulse, the x rays excite inner-shell electrons leading to the ionization of the electrons through various atomic processes and creating high-energy-density plasma, i.e., warm or hot dense matter. The resulting system consists of atoms in various electronic configurations, thermalizing on subpicosecond to picosecond timescales after photoexcitation. We present a simulation study of x-ray-heated solid-density matter. For this we use XMDYN, a Monte Carlo molecular-dynamics-based code with periodic boundary conditions, which allows one to investigate nonequilibrium dynamics. XMDYN is capable of treating systems containing light and heavy atomic species with full electronic configuration space and three-dimensional spatial inhomogeneity. For the validation of our approach we compare for a model system the electron temperatures and the ion charge-state distribution from XMDYN to results for the thermalized system based on the average-atom model implemented in XATOM, an ab initio x-ray atomic physics toolkit extended to include a plasma environment. Further, we also compare the average charge evolution of diamond with the predictions of a Boltzmann continuum approach. We demonstrate that XMDYN results are in good quantitative agreement with the above-mentioned approaches, suggesting that the current implementation of XMDYN is a viable approach to simulate the dynamics of x-ray-driven nonequilibrium dynamics in solids. To illustrate the potential of XMDYN for treating complex systems, we present calculations on the triiodo benzene derivative 5-amino-2,4,6-triiodoisophthalic acid (I3C), a compound of relevance of biomolecular imaging, consisting of heavy and light atomic species.

Start-to-end simulation of single-particle imaging using ultra-short pulses at the European X-ray Free-Electron Laser.

Single-particle imaging with X-ray free-electron lasers (XFELs) has the potential to provide structural information at atomic resolution for non-crystalline biomolecules. This potential exists because ultra-short intense pulses can produce interpretable diffraction data notwithstanding radiation damage. This paper explores the impact of pulse duration on the interpretability of diffraction data using comprehensive and realistic simulations of an imaging experiment at the European X-ray Free-Electron Laser. It is found that the optimal pulse duration for molecules with a few thousand atoms at 5 keV lies between 3 and 9 fs.

Calculation of x-ray scattering patterns from nanocrystals at high x-ray intensity.

We present a generalized method to describe the x-ray scattering intensity of the Bragg spots in a diffraction pattern from nanocrystals exposed to intense x-ray pulses. Our method involves the subdivision of a crystal into smaller units. In order to calculate the dynamics within every unit, we employ a Monte-Carlo-molecular dynamics-ab-initio hybrid framework using real space periodic boundary conditions. By combining all the units, we simulate the diffraction pattern of a crystal larger than the transverse x-ray beam profile, a situation commonly encountered in femtosecond nanocrystallography experiments with focused x-ray free-electron laser radiation. Radiation damage is not spatially uniform and depends on the fluence associated with each specific region inside the crystal. To investigate the effects of uniform and non-uniform fluence distribution, we have used two different spatial beam profiles, Gaussian and flattop.

A comprehensive simulation framework for imaging single particles and biomolecules at the European X-ray Free-Electron Laser.

The advent of newer, brighter, and more coherent X-ray sources, such as X-ray Free-Electron Lasers (XFELs), represents a tremendous growth in the potential to apply coherent X-rays to determine the structure of materials from the micron-scale down to the Angstrom-scale. There is a significant need for a multi-physics simulation framework to perform source-to-detector simulations for a single particle imaging experiment, including (i) the multidimensional simulation of the X-ray source; (ii) simulation of the wave-optics propagation of the coherent XFEL beams; (iii) atomistic modelling of photon-material interactions; (iv) simulation of the time-dependent diffraction process, including incoherent scattering; (v) assembling noisy and incomplete diffraction intensities into a three-dimensional data set using the Expansion-Maximisation-Compression (EMC) algorithm and (vi) phase retrieval to obtain structural information. We demonstrate the framework by simulating a single-particle experiment for a nitrogenase iron protein using parameters of the SPB/SFX instrument of the European XFEL. This exercise demonstrably yields interpretable consequences for structure determination that are crucial yet currently unavailable for experiment design.