Femtosecond X-ray-induced explosion of C60 at extreme intensity

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.

The beam transport system for the Small Quantum Systems instrument at the European XFEL: optical layout and first commissioning results

The Small Quantum Systems instrument is one of the six operating instruments of the European XFEL, dedicated to the atomic, molecular and cluster physics communities. The instrument started its user operation at the end of 2018 after a commissioning phase. The design and characterization of the beam transport system are described here. The X-ray optical components of the beamline are detailed, and the beamline performances, transmission and focusing capabilities are reported. It is shown that the X-ray beam can be effectively focused as predicted by ray-tracing simulations. The impact of non-ideal X-ray source conditions on the focusing performances is discussed.

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.

Probing C60 Fullerenes from within Using Free Electron Lasers

Fullerenes, such as C60, are ideal systems to investigate energy redistribution following substantial excitation. Ultra-short and ultra-intense free electron lasers (FELs) have allowed molecular research in a new photon energy regime. FELs have allowed the study of the response of fullerenes to X-rays, which includes femtosecond multi-photon processes, as well as time-resolved ionization and fragmentation dynamics. This perspective: (1) provides a general introduction relevant to C60 research using photon sources, (2) reports on two specific X-ray FEL-based photoionization investigations of C60, at two different FEL fluences, one static and one time-resolved, and (3) offers a brief analysis and recommendations for future research.

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.

Complex energies and transition dipoles for shape-type resonances of uracil anion from stabilization curves via Padé

Absorption of slow moving electrons by neutral ground state nucleobases has been known to produce resonance metastable states. There are indications that such metastable states may play a key role in DNA/RNA damage. Therefore, herein, we present an ab initio non-Hermitian investigation of the resonance positions and decay rates for the low lying shape-type states of the uracil anion. In addition, we calculate the complex transition dipoles between these resonance states. We employ the resonance via Padé (RVP) method to calculate these complex properties from real stabilization curves by analytical dilation into the complex plane. This method has already been successfully applied to many small molecular systems, and herein, we present the first application of RVP to a medium-sized system. The presented resonance energies are optimized with respect to the size of the basis set and compared with previous theoretical studies and experimental findings. Complex transition dipoles between the shape-type resonances are computed using the optimal basis set. The ability to calculate ab initio energies and lifetimes of biologically relevant systems paves the way for studying reactions of such systems in which autoionization takes place, while the ability to also calculate their complex transition dipoles opens the door for studying photo-induced dynamics of such biological molecules.

Synchrotron X-ray-induced Synthesis of Copper Hydroxide Nitrate Nanoplates on Cu Thin Films in an Ambient Atmosphere

Synchrotron X-rays are widely used for material characterizations. However, they can also ionize atoms and molecules to damage and manipulate probed materials. We report here an X-ray-induced growth of copper hydroxide nitrate, Cu2(OH)3NO3, on copper thin films in the ambient atmosphere without solvents and thermal treatment. In situ synchrotron X-ray diffraction measurements showed that the time-dependent growth process of theCu2(OH)3NO3 is accompanied by the consumption of Cu metal and can be described by a sigmoidal model. The growth rate was reduced after the initial fast growth period. Scanning electron microscopy (SEM) images show that the isolated islands of Cu2(OH)3NO3 nanoplates formed in the beginning, which grew together with new nanoplates formed under continued X-ray irradiation. The result demonstrated that high-flux synchrotron X-rays may provide an unconventional approach to synthesizing and manipulating materials, which will inspire future investigation both experimentally and theoretically.

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.

Limitations of Structural Insight into Ultrafast Melting of Solid Materials with X-ray Diffraction Imaging

In this work, we analyze the application of X-ray diffraction imaging techniques to follow ultrafast structural transitions in solid materials using the example of an X-ray pump–X-ray probe experiment with a single-crystal silicon performed at a Linac Coherent Light Source. Due to the spatially non-uniform profile of the X-ray beam, the diffractive signal recorded in this experiment included contributions from crystal parts experiencing different fluences from the peak fluence down to zero. With our theoretical model, we could identify specific processes contributing to the silicon melting in those crystal regions, i.e., the non-thermal and thermal melting whose occurrences depended on the locally absorbed X-ray doses. We then constructed the total volume-integrated signal by summing up the coherent signal contributions (amplitudes) from the various crystal regions and found that this significantly differed from the signals obtained for a few selected uniform fluence values, including the peak fluence. This shows that the diffraction imaging signal obtained for a structurally damaged material after an impact of a non-uniform X-ray pump pulse cannot be always interpreted as the material’s response to a pulse of a specific (e.g., peak) fluence as it is sometimes believed. This observation has to be taken into account in planning and interpreting future experiments investigating structural changes in materials with X-ray diffraction imaging.

Response of free-standing graphene monolayer exposed to ultrashort intense XUV pulse from free-electron laser

The response of a free-standing graphene monolayer exposed to a few tens of femtoseconds long extreme ultraviolet (XUV) pulse was studied theoretically in order to analyze and compare contributions of various mechanisms to the graphene damage, understood here as a global atomic disintegration. Our simulation results indicate that nonthermal disintegration of the atomic structure is the predominant damage mechanism for a free-standing graphene layer. Only at high absorbed doses, charge-induced disintegration of the graphene structure prevails. We also demonstrate that the progressing damage can be probed by femtosecond optical pulses in the soft UV regime (4 eV photon energy). The achieved quantitative understanding of the damage mechanisms may enable a better control of graphene-based devices when they are exposed to x-ray radiation, as well as an efficient processing of graphene layers with ultrashort intense XUV pulses.

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.

Uniform vs Partial Scaling within Resonances via Padé Based on the Similarities to Other Non-Hermitian Methods: Illustration for the Beryllium 1s22p3s State

Resonance via Padé (RVP) is an efficient method for calculating autoionization resonance states. It is based on the stabilization technique in which the basis set is scaled. The scaling can be uniform (i.e., all basis functions are scaled) or partial. Herein, we compare the two RVP scaling schemes for calculating an autoionization eigenvalue; moreover, the effect of freezing the core electrons is intertwined within this comparison. In order to study the different behavior of the RVP schemes, we associate each RVP scaling scheme with a complex contour of integration. Similarities between RVP and other non-Hermitian methods emerge from the generated contours, which suggest that RVP introduces similar outgoing boundary conditions as the complex scaling (CS), complex basis function (CBF), and reflection-free complex absorbing potential (RF-CAP) methods. A uniform-RVP contour, unlike a partial one, immediately penetrates the complex plane and influences the interaction region. Hence, uniform scaling within RVP destroys the description of the core electrons, as well as the description of the reference state, and yields less reliable results than partial scaling. The 1s²2p3s ¹P autoionization state of Be, at the equation-of-motion coupled-cluster level, is used as our case study model.

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.

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.

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.

Capturing the photo-induced dynamics of nano-molecules by X-ray free electron laser induced Coulomb explosion

We performed reaction dynamics simulations to demonstrate that the vibrational dynamics of C₆₀ induced by infrared (IR) pulses can be traced by triggering Coulomb explosion with intense femtosecond X-ray free electron laser (XFEL) probe pulses. The time series of the angular anisotropy β(t) of fast C⁺ and C²⁺ fragments of C₆₀ ⁶⁰⁺ produced by such an XFEL pulse reflects the instantaneous structure of C₆₀ vibrationally excited by IR pulses. The phases and amplitudes of excited vibrational modes and the coupling between excited modes can be successfully extracted from the expansion of β(t) in terms of vibrational modes. This proof-of-principle simulation clearly demonstrates that various information of the structures and reaction dynamics of large clusters or biomolecules can be retrieved by decomposing the experimentally determined β(t) into vibrational modes.

Photo-ionization and fragmentation of Sc3N@C80 following excitation above the Sc K-edge

We have investigated the ionization and fragmentation of a metallo-endohedral fullerene, Sc₃N@C₈₀, using ultrashort (10 fs) x-ray pulses. Following selective ionization of a Sc (1s) electron (hν = 4.55 keV), an Auger cascade leads predominantly to either a vibrationally cold multiply charged parent molecule or multifragmentation of the carbon cage following a phase transition. In contrast to previous studies, no intermediate regime of C₂ evaporation from the carbon cage is observed. A time-delayed, hard x-ray pulse (hν = 5.0 keV) was used to attempt to probe the electron transfer dynamics between the encapsulated Sc species and the carbon cage. A small but significant change in the intensity of Sc-containing fragment ions and coincidence counts for a delay of 100 fs compared to 0 fs, as well as an increase in the yield of small carbon fragment ions, may be indicative of incomplete charge transfer from the carbon cage on the sub-100 fs time scale.

Is radiation damage the limiting factor in high-resolution single particle imaging with X-ray free-electron lasers?

The prospect of single particle imaging with atomic resolution is one of the scientific drivers for the development of X-ray free-electron lasers. The assumption since the beginning has been that damage to the sample caused by intense X-ray pulses is one of the limiting factors for achieving subnanometer X-ray imaging of single particles and that X-ray pulses need to be as short as possible. Based on the molecular dynamics simulations of proteins in X-ray fields of various durations (5 fs, 25 fs, and 50 fs), we show that the noise in the diffracted signal caused by radiation damage is less than what can be expected from other sources, such as sample inhomogeneity and X-ray shot-to-shot variations. These findings show a different aspect of the feasibility of high-resolution single particle imaging using free-electron lasers, where employing X-ray pulses of longer durations could still provide a useful diffraction signal above the noise due to the Coulomb explosion.

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.

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.

Imaging electron-density fluctuations by multidimensional X-ray photon-coincidence diffraction

The ultrafast spontaneous electron-density fluctuation dynamics in molecules is studied theoretically by off-resonant multiple X-ray diffraction events. The time- and wavevector-resolved photon-coincidence signals give an image of electron-density fluctuations expressed through the four-point correlation function of the charge density in momentum space. A Fourier transform of the signal provides a real-space image of the multipoint charge-density correlation functions, which reveal snapshots of the evolving electron density in between the diffraction events. The proposed technique is illustrated by ab initio simulations of the momentum- and real-space inelastic scattering signals from a linear cyanotetracetylene molecule.

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.

The LAMP instrument at the Linac Coherent Light Source free-electron laser

The Laser Applications in Materials Processing (LAMP) instrument is a new end-station for soft X-ray imaging, high-field physics, and ultrafast X-ray science experiments that is available to users at the Linac Coherent Light Source (LCLS) free-electron laser. While the instrument resides in the Atomic, Molecular and Optical science hutch, its components can be used at any LCLS beamline. The end-station has a modular design that provides high flexibility in order to meet user-defined experimental requirements and specifications. The ultra-high-vacuum environment supports different sample delivery systems, including pulsed and continuous atomic, molecular, and cluster jets; liquid and aerosols jets; and effusive metal vapor beams. It also houses movable, large-format, high-speed pnCCD X-ray detectors for detecting scattered and fluorescent photons. Multiple charged-particle spectrometer options are compatible with the LAMP chamber, including a double-sided spectrometer for simultaneous and even coincident measurements of electrons, ions, and photons produced by the interaction of the high-intensity X-ray beam with the various samples. Here we describe the design and capabilities of the spectrometers along with some general aspects of the LAMP chamber and show some results from the initial instrument commissioning.

XUV/X-ray light and fast ions for ultrafast chemistry

The deposition of large amounts of energy in a molecule by XUV/X-ray photon absorption or fast-ion collision, triggers a set of complex ultrafast electronic and nuclear dynamics that allow a deep understanding and control of chemical reactivity. This themed issue showcases the research performed in the understanding, monitoring and control of these processes.

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.

A perspective for investigating photo-induced molecular dynamics from within with femtosecond free electron lasers

Photo-induced molecular dynamics can now be investigated using free electron lasers (FELs) whose attributes are unprecedented brightness, few femtosecond pulses duration and in the near future few hundreds of attosecond pulse duration.

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.

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.

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.

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

Rapid development of X-ray free-electron laser (XFEL) science has taken place in recent years owing to the consecutive launch of large-scale XFEL instruments around the world. Research areas such as warm dense matter physics and coherent X-ray imaging take advantage of the unprecedentedly high intensities of XFELs. A single XFEL pulse can induce very complex dynamics within matter initiated by core-hole photoionization. Owing to this complexity, theoretical modeling revealing details of the excitation and relaxation of irradiated matter is important for the correct interpretation of the measurements and for proposing new experiments.XMDYNis a computer simulation tool developed for modeling dynamics of matter induced by high-intensity X-rays. It utilizes atomic data calculated by theab initio XATOMtoolkit. Here these tools are discussed in detail.

Towards phasing using high X-ray intensity

X-ray free-electron lasers (XFELs) show great promise for macromolecular structure determination from sub-micrometre-sized crystals, using the emerging method of serial femtosecond crystallography. The extreme brightness of the XFEL radiation can multiply ionize most, if not all, atoms in a protein, causing their scattering factors to change during the pulse, with a preferential 'bleaching' of heavy atoms. This paper investigates the effects of electronic damage on experimental data collected from a Gd derivative of lysozyme microcrystals at different X-ray intensities, and the degree of ionization of Gd atoms is quantified from phased difference Fourier maps. A pattern sorting scheme is proposed to maximize the ionization contrast and the way in which the local electronic damage can be used for a new experimental phasing method is discussed.

Goniometer-based femtosecond X-ray diffraction of mutant 30S ribosomal subunit crystals

In this work, we collected radiation-damage-free data from a set of cryo-cooled crystals for a novel 30S ribosomal subunit mutant using goniometer-based femtosecond crystallography. Crystal quality assessment for these samples was conducted at the X-ray Pump Probe end-station of the Linac Coherent Light Source (LCLS) using recently introduced goniometer-based instrumentation. These 30S subunit crystals were genetically engineered to omit a 26-residue protein, Thx, which is present in the wild-type Thermus thermophilus 30S ribosomal subunit. We are primarily interested in elucidating the contribution of this ribosomal protein to the overall 30S subunit structure. To assess the viability of this study, femtosecond X-ray diffraction patterns from these crystals were recorded at the LCLS during a protein crystal screening beam time. During our data collection, we successfully observed diffraction from these difficult-to-grow 30S ribosomal subunit crystals. Most of our crystals were found to diffract to low resolution, while one crystal diffracted to 3.2 Å resolution. These data suggest the feasibility of pursuing high-resolution data collection as well as the need to improve sample preparation and handling in order to collect a complete radiation-damage-free data set using an X-ray Free Electron Laser.

Efficient electronic structure calculation for molecular ionization dynamics at high x-ray intensity

We present the implementation of an electronic-structure approach dedicated to ionization dynamics of molecules interacting with x-ray free-electron laser (XFEL) pulses. In our scheme, molecular orbitals for molecular core-hole states are represented by linear combination of numerical atomic orbitals that are solutions of corresponding atomic core-hole states. We demonstrate that our scheme efficiently calculates all possible multiple-hole configurations of molecules formed during XFEL pulses. The present method is suitable to investigate x-ray multiphoton multiple ionization dynamics and accompanying nuclear dynamics, providing essential information on the chemical dynamics relevant for high-intensity x-ray imaging.

Nanoplasma Formation by High Intensity Hard X-rays

Using electron spectroscopy, we have investigated nanoplasma formation from noble gas clusters exposed to high-intensity hard-x-ray pulses at ~5 keV. Our experiment was carried out at the SPring-8 Angstrom Compact free electron LAser (SACLA) facility in Japan. Dedicated theoretical simulations were performed with the molecular dynamics tool XMDYN. We found that in this unprecedented wavelength regime nanoplasma formation is a highly indirect process. In the argon clusters investigated, nanoplasma is mainly formed through secondary electron cascading initiated by slow Auger electrons. Energy is distributed within the sample entirely through Auger processes and secondary electron cascading following photoabsorption, as in the hard x-ray regime there is no direct energy transfer from the field to the plasma. This plasma formation mechanism is specific to the hard-x-ray regime and may, thus, also be important for XFEL-based molecular imaging studies. In xenon clusters, photo- and Auger electrons contribute more significantly to the nanoplasma formation. Good agreement between experiment and simulations validates our modelling approach. This has wide-ranging implications for our ability to quantitatively predict the behavior of complex molecular systems irradiated by high-intensity hard x-rays.

Efficient autoionization following intense laser-cluster interactions

Electron emission as a result of the interaction of clusters with intense laser pulses is commonly understood in terms of direct and evaporative ionization processes. In contrast, we provide evidence here of an important role played by autoionization in intense field ionization of molecular oxygen clusters. Superexcited states are populated during the cluster expansion, and their autoionization is observed on a ns time scale. Decay processes on fs to ps time scales are obscured by energy exchange of the emitted electrons with the environment.