XMDYNandXATOM: versatile simulation tools for quantitative modeling of X-ray free-electron laser induced dynamics of matter

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.

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.

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.

Multiple-core-hole resonance spectroscopy with ultraintense X-ray pulses

Understanding the interaction of intense, femtosecond X-ray pulses with heavy atoms is crucial for gaining insights into the structure and dynamics of matter. One key aspect of nonlinear light-matter interaction was, so far, not studied systematically at free-electron lasers-its dependence on the photon energy. Here, we use resonant ion spectroscopy to map out the transient electronic structures occurring during the complex charge-up pathways of xenon. Massively hollow atoms featuring up to six simultaneous core holes determine the spectra at specific photon energies and charge states. We also illustrate how different X-ray pulse parameters, which are usually intertwined, can be partially disentangled. The extraction of resonance spectra is facilitated by the possibility of working with a constant number of photons per X-ray pulse at all photon energies and the fact that the ion yields become independent of the peak fluence beyond a saturation point. Our study lays the groundwork for spectroscopic investigations of transient atomic species in exotic, multiple-core-hole states that have not been explored previously.

Theoretical study on time response of semiconductor photorefractive effects under subpicosecond ultra-fast X-rays

A theoretical model that can efficiently calculate the refractive index response of semiconductors under ultrafast X-ray radiation is established based on the photorefractive effect of semiconductors. The proposed model is used to interpret X-ray diagnostics experiments, and the results are in good agreement with experiments. In the proposed model, a rate equation model of free carrier density calculation is adopted with the X-ray absorption cross-sections calculated by atomic codes. The two-temperature model is used to describe the electron-lattice equilibration and the extended Drude model is applied to calculate the transient refractive index change. It is found that faster time response can be achieved for semiconductors with shorter carrier lifetime and sub-picosecond resolution can be obtained for InP and [Formula: see text]. The material response time is not sensitive to X-ray energy and the diagnostics can be used in the 1-10 keV energy range. This article is part of the theme issue 'Dynamic and transient processes in warm dense matter'.

Application of Boltzmann kinetic equations to model X-ray-created warm dense matter and plasma

In this review, we describe the application of Boltzmann kinetic equations for modelling warm dense matter and plasma formed after irradiation of solid materials with intense femtosecond X-ray pulses. Classical Boltzmann kinetic equations are derived from the reduced N-particle Liouville equations. They include only single-particle densities of ions and free electrons present in the sample. The first version of the Boltzmann kinetic equation solver was completed in 2006. It could model non-equilibrium evolution of X-ray-irradiated finite-size atomic systems. In 2016, the code was adapted to study plasma created from X-ray-irradiated materials. Additional extension of the code was then also performed, enabling simulations in the hard X-ray irradiation regime. In order to avoid treatment of a very high number of active atomic configurations involved in the excitation and relaxation of X-ray-irradiated materials, an approach called 'predominant excitation and relaxation path' (PERP) was introduced. It limited the number of active atomic configurations by following the sample evolution only along most PERPs. The performance of the Boltzmann code is illustrated in the examples of X-ray-heated solid carbon and gold. Actual model limitations and further model developments are discussed. This article is part of the theme issue 'Dynamic and transient processes in warm dense matter'.

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.

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.

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.

Nuclear Excitation by Electron Capture in Excited Ions

A nuclear excitation following the capture of an electron in an empty orbital has been recently observed for the first time. So far, the evaluation of the cross section of the process has been carried out widely using the assumption that the ion is in its electronic ground state prior to the capture. We show that by lifting this restriction new capture channels emerge resulting in a boost of more than three orders of magnitude to the electron capture resonance strength.

Approaching the Attosecond Frontier of Dynamics in Matter with the Concept of X-ray Chronoscopy

X-ray free electron lasers (XFELs) have provided scientists opportunities to study matter with unprecedented temporal and spatial resolutions. However, access to the attosecond domain (i.e., below 1 femtosecond) remains elusive. Herein, a time-dependent experimental concept is theorized, allowing us to track ultrafast processes in matter with sub-fs resolution. The proposed X-ray chronoscopy approach exploits the state-of-the-art developments in terahertz streaking to measure the time structure of X-ray pulses with ultrahigh temporal resolution. The sub-femtosecond dynamics of the saturable X-ray absorption process is simulated. The employed rate equation model confirms that the X-ray-induced mechanisms leading to X-ray transparency can be probed via measurement of an X-ray pulse time structure.

Resonance-Enhanced Multiphoton Ionization in the X-Ray Regime

Here, we report on the nonlinear ionization of argon atoms in the short wavelength regime using ultraintense x rays from the European XFEL. After sequential multiphoton ionization, high charge states are obtained. For photon energies that are insufficient to directly ionize a 1s electron, a different mechanism is required to obtain ionization to Ar^{17+}. We propose this occurs through a two-color process where the second harmonic of the FEL pulse resonantly excites the system via a 1s→2p transition followed by ionization by the fundamental FEL pulse, which is a type of x-ray resonance-enhanced multiphoton ionization (REMPI). This resonant phenomenon occurs not only for Ar^{16+}, but also through lower charge states, where multiple ionization competes with decay lifetimes, making x-ray REMPI distinctive from conventional REMPI. With the aid of state-of-the-art theoretical calculations, we explain the effects of x-ray REMPI on the relevant ion yields and spectral profile.

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.

Ultraintense, ultrashort pulse X-ray scattering in small molecules

We examine X-ray scattering from an isolated organic molecule from the linear to nonlinear absorptive regime. In the nonlinear regime, we explore the importance of both the coherent and incoherent channels and observe the onset of nonlinear behavior as a function of pulse duration and energy. In the linear regime, we test the sensitivity of the scattering signal to molecular bonding and electronic correlation via calculations using the independent atom model (IAM), Hartree-Fock (HF) and density functional theory (DFT). Finally, we describe how coherent X-ray scattering can be used to directly visualize femtosecond charge transfer and dissociation within a single molecule undergoing X-ray multiphoton absorption.

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.

Evaluation-oriented exploration of photo energy conversion systems: from fundamental optoelectronics and material screening to the combination with data science

Light is a form of energy that can be converted to electric and chemical energies. Thus, organic photovoltaics (OPVs), perovskite solar cells (PSCs), photocatalysts, and photodetectors have evolved as scientific and commercial enterprises. However, the complex photochemical reactions and multicomponent materials involved in these systems have hampered rapid progress in their fundamental understanding and material design. This review showcases the evaluation-oriented exploration of photo energy conversion materials by using electrodeless time-resolved microwave conductivity (TRMC) and materials informatics (MI). TRMC with its unique options (excitation sources, environmental control, frequency modulation, etc.) provides not only accelerated experimental screening of OPV and PSC materials but also a versatile route toward shedding light on their charge carrier dynamics. Furthermore, MI powered by machine learning is shown to allow extremely high-throughput exploration in the large molecular space, which is compatible with experimental screening and combinatorial synthesis.

Population inversion X-ray laser oscillator

Oscillators are at the heart of optical lasers, providing stable, transform-limited pulses. Until now, laser oscillators have been available only in the infrared to visible and near-ultraviolet (UV) spectral region. In this paper, we present a study of an oscillator operating in the 5- to 12-keV photon-energy range. We show that, using the [Formula: see text] line of transition metal compounds as the gain medium, an X-ray free-electron laser as a periodic pump, and a Bragg crystal optical cavity, it is possible to build X-ray oscillators producing intense, fully coherent, transform-limited pulses. As an example, we consider the case of a copper nitrate gain medium generating ∼ 5 × [Formula: see text] photons per pulse with 37-fs pulse length and 48-meV spectral resolution at 8-keV photon energy. Our theoretical study and simulation of this system show that, because of the very large gain per pass, the oscillator saturates and reaches full coherence in four to six optical-cavity transits. We discuss the feasibility and design of the X-ray optical cavity and other parts of the oscillator needed for its realization, opening the way to extend X-ray-based research beyond current capabilities.

A Perspective on Molecular Structure and Bond-Breaking in Radiation Damage in Serial Femtosecond Crystallography

X-ray free-electron lasers (XFELs) have a unique capability for time-resolved studies of protein dynamics and conformational changes on femto- and pico-second time scales. The extreme intensity of X-ray pulses can potentially cause significant modifications to the sample structure during exposure. Successful time-resolved XFEL crystallography depends on the unambiguous interpretation of the protein dynamics of interest from the effects of radiation damage. Proteins containing relatively heavy elements, such as sulfur or metals, have a higher risk for radiation damage. In metaloenzymes, for example, the dynamics of interest usually occur at the metal centers, which are also hotspots for damage due to the higher atomic number of the elements they contain. An ongoing challenge with such local damage is to understand the residual bonding in these locally ionized systems and bond-breaking dynamics. Here, we present a perspective on radiation damage in XFEL experiments with a particular focus on the impacts for time-resolved protein crystallography. We discuss recent experimental and modelling results of bond-breaking and ion motion at disulfide bonding sites in protein crystals.

Hybrid Plasma/Molecular-Dynamics Approach for Efficient XFEL Radiation Damage Simulations

X-ray free-electron laser pulses initiate a complex series of changes to the electronic and nuclear structure of matter on femtosecond timescales. These damage processes include widespread ionization, the formation of a quasi-plasma state and the ultimate explosion of the sample due to Coulomb forces. The accurate simulation of these dynamical effects is critical in designing feasible XFEL experiments and interpreting the results. Current molecular dynamics simulations are, however, computationally intensive, particularly when they treat unbound electrons as classical point particles. On the other hand, plasma simulations are computationally efficient but do not model atomic motion. Here we present a hybrid approach to XFEL damage simulation that combines molecular dynamics for the nuclear motion and plasma models to describe the evolution of the low-energy electron continuum. The plasma properties of the unbound electron gas are used to define modified inter-ionic potentials for the molecular dynamics, including Debye screening and drag forces. The hybrid approach is significantly faster than damage simulations that treat unbound electrons as classical particles, enabling simulations to be performed on large sample volumes.

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.

The local X-ray-induced dynamics that occur in protein crystals during serial femtosecond crystallography (SFX) measurements at XFELs are not well understood. Here the authors performed a time-resolved X-ray pump X-ray probe SFX experiment, and they observe distinct structural changes in the disulfide bridges and peptide backbone of proteins; complementing theoretical approaches allow them to further characterize the details of the X-ray induced ionization and local structural dynamics.

The role of transient resonances for ultra-fast imaging of single sucrose nanoclusters

Intense x-ray free-electron laser (XFEL) pulses hold great promise for imaging function in nanoscale and biological systems with atomic resolution. So far, however, the spatial resolution obtained from single shot experiments lags averaging static experiments. Here we report on a combined computational and experimental study about ultrafast diffractive imaging of sucrose clusters which are benchmark organic samples. Our theoretical model matches the experimental data from the water window to the keV x-ray regime. The large-scale dynamic scattering calculations reveal that transient phenomena driven by non-linear x-ray interaction are decisive for ultrafast imaging applications. Our study illuminates the complex interplay of the imaging process with the rapidly changing transient electronic structures in XFEL experiments and shows how computational models allow optimization of the parameters for ultrafast imaging experiments.

X-ray free electron lasers provide high photon flux to explore single particle diffraction imaging of biological samples. Here the authors present dynamic electronic structure calculations and benchmark them to single-particle XFEL diffraction data of sucrose clusters to predict optimal single-shot imaging conditions.

Core-level nonlinear spectroscopy triggered by stochastic X-ray pulses

Stochastic processes are highly relevant in research fields as different as neuroscience, economy, ecology, chemistry, and fundamental physics. However, due to their intrinsic unpredictability, stochastic mechanisms are very challenging for any kind of investigations and practical applications. Here we report the deliberate use of stochastic X-ray pulses in two-dimensional spectroscopy to the simultaneous mapping of unoccupied and occupied electronic states of atoms in a regime where the opacity and transparency properties of matter are subject to the incident intensity and photon energy. A readily transferable matrix formalism is presented to extract the electronic states from a dataset measured with the monitored input from a stochastic excitation source. The presented formalism enables investigations of the response of the electronic structure to irradiation with intense X-ray pulses while the time structure of the incident pulses is preserved.

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.

X-ray Emission Spectroscopy at X-ray Free Electron Lasers: Limits to Observation of the Classical Spectroscopic Response for Electronic Structure Analysis

X-ray free electron lasers (XFELs) provide ultrashort intense X-ray pulses suitable to probe electron dynamics but can also induce a multitude of nonlinear excitation processes. These affect spectroscopic measurements and interpretation, particularly for upcoming brighter XFELs. Here we identify and discuss the limits to observing classical spectroscopy, where only one photon is absorbed per atom for a Mn²⁺ in a light element (O, C, H) environment. X-ray emission spectroscopy (XES) with different incident photon energies, pulse intensities, and pulse durations is presented. A rate equation model based on sequential ionization and relaxation events is used to calculate populations of multiply ionized states during a single pulse and to explain the observed X-ray induced spectral lines shifts. This model provides easy estimation of spectral shifts, which is essential for experimental designs at XFELs and illustrates that shorter X-ray pulses will not overcome sequential ionization but can reduce electron cascade effects.

Simulated XUV photoelectron spectra of THz-pumped liquid water

Highly intense, sub-picosecond terahertz (THz) pulses can be used to induce ultrafast temperature jumps (T-jumps) in liquid water. A supercritical state of gas-like water with liquid density is established, and the accompanying structural changes are expected to give rise to time-dependent chemical shifts. We investigate the possibility of using extreme ultraviolet photoelectron spectroscopy as a probe for ultrafast dynamics induced by sub-picosecond THz pulses of varying intensities and frequencies. To this end, we use ab initio methods to calculate photoionization cross sections and photoelectron energies of (H₂O)₂₀ clusters embedded in an aqueous environment represented by point charges. The cluster geometries are sampled from ab initio molecular dynamics simulations modeling the THz-water interactions. We find that the peaks in the valence photoelectron spectrum are shifted by up to 0.4 eV after the pump pulse and that they are broadened with respect to unheated water. The shifts can be connected to structural changes caused by the heating, but due to saturation effects they are not sensitive enough to serve as a thermometer for T-jumped water.

Ultrafast x-ray-driven phenomena in nanocrystals: development and application of powerful simulation tools

We investigate the radiation damage dynamics of nanocrystals at high x-ray intensity, by using time-resolved scattering patterns. We present dynamics simulations for biologically relevant molecules using XMDYN extended to nanocrystals and scattering simulation with XSINC.

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.

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.

Modeling of the interaction of an x-ray free-electron laser with large finite samples

We describe a model for the study of the interaction of short x-ray free-electron laser (XFEL) pulses with large finite samples. Hydrodynamics is solved in one-dimensional planar geometry with consideration of the electron-ion energy exchange and of the possible elastoplastic behavior. From a time-dependent calculation of the complex refractive index and of the underlying atomic physics, XFEL energy deposition is modeled through a calculation of the radiation field in the material. In the case of hard x-ray irradiation, energetic electrons induced by the XFEL absorption can propagate and deposit their energy outside the interaction region. Simulations of the interaction of hard x-ray ultrashort pulses with solid materials Ru and Si at different grazing incidence angles are presented and discussed. The results obtained demonstrate the potential of this approach to predict damage dynamics for materials of interest for x-ray optics.

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.

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.

Interaction of molecular nitrogen with free-electron-laser radiation

Molecular double core hole contribution to the final atomic ion fragments of N₂ when driven by an FEL pulse.

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.

CCP-FEL: a collection of computer programs for free-electron laser research

The latest virtual special issue ofJournal of Applied Crystallography(http://journals.iucr.org/special_issues/2016/ccpfel) collects software for free-electron laser research and presents tools for a range of topics such as simulation of experiments, online monitoring of data collection, selection of hits, diagnostics of data quality, data management, data analysis and structure determination for both nanocrystallography and single-particle diffractive imaging. This article provides an introduction to the special issue.

Kinetic Boltzmann approach adapted for modeling highly ionized matter created by x-ray irradiation of a solid

We report on the kinetic Boltzmann approach adapted for simulations of highly ionized matter created from a solid by its x-ray irradiation. X rays can excite inner-shell electrons, which leads to the creation of deeply lying core holes. Their relaxation, especially in heavier elements, can take complicated paths, leading to a large number of active configurations. Their number can be so large that solving the set of respective evolution equations becomes computationally inefficient and another modeling approach should be used instead. To circumvent this complexity, the commonly used continuum models employ a superconfiguration scheme. Here, we propose an alternative approach which still uses "true" atomic configurations but limits their number by restricting the sample relaxation to the predominant relaxation paths. We test its reliability, performing respective calculations for a bulk material consisting of light atoms and comparing the results with a full calculation including all relaxation paths. Prospective application for heavy elements is discussed.