A particle-particle grand canonical Monte Carlo model with Coulomb pair potential interaction was used to simulate modification of ionization potentials by electrostatic microfields. The Barnes-Hut tree algorithm [J. Barnes and P. Hut, Nature 324, 446 (1986)] was used to speed up calculations of electric potential. Atomic levels were approximated to be independent of the microfields as was assumed in the original paper by Ecker and Kröll [Phys. Fluids 6, 62 (1963)]; however, the available levels were limited by the corresponding mean inter-particle distance. The code was tested on hydrogen and dense aluminum plasmas. The amount of depression was up to 50% higher in the Debye-Hückel regime for hydrogen plasmas, in the high density limit, reasonable agreement was found with the Ecker-Kröll model for hydrogen plasmas and with the Stewart-Pyatt model [J. Stewart and K. Pyatt, Jr., Astrophys. J. 144, 1203 (1966)] for aluminum plasmas. Our 3D code is an improvement over the spherically symmetric simplifications of the Ecker-Kröll and Stewart-Pyatt models and is also not limited to high atomic numbers as is the underlying Thomas-Fermi model used in the Stewart-Pyatt model.

Warm dense matter is at the boundary between a plasma and a condensed phase and plays a role in astrophysics, planetary science and inertial confinement fusion research. However, its electronic structure and ionic structure upon irradiation with strong laser pulses remain poorly understood. Here, we use an intense and ultrafast X-ray free-electron laser pulse to simultaneously create and characterize warm dense copper using L-edge X-ray absorption spectroscopy over a large irradiation intensity range. Below a pulse intensity of 10¹⁵ W cm⁻², an absorption peak below the L edge appears, originating from transient depletion of the 3d band. This peak shifts to lower energy with increasing intensity, indicating the movement of the 3d band upon strong X-ray excitation. At higher intensities, substantial ionization and collisions lead to the transition from reverse saturable absorption to saturable absorption of the X-ray free-electron laser pulse, two nonlinear effects that hold promise for X-ray pulse-shaping. We employ theoretical calculations that combine a model based on kinetic Boltzmann equations with finite-temperature real-space density-functional theory to interpret these observations. The results can be used to benchmark non-equilibrium models of electronic structure in warm dense matter.

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.

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.

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.

In the warm dense matter (WDM) regime, where condensed, gas, and plasma phases coexist, matter frequently exhibits unusual properties that cannot be described by contemporary theory. Experiments reporting phenomena in WDM are therefore of interest to advance our physical understanding of this regime, which is found in dwarf stars, giant planets, and fusion ignition experiments. Using 7.1 keV X-ray free electron laser radiation (nominally 5×10⁵J/cm²), we produced and probed transient WDM in liquid water. Wide-angle X-ray scattering (WAXS) from the probe reveals a new ~9 Å structure that forms within 75 fs. By 100 fs, the WAXS peak corresponding to this new structure is of comparable magnitude to the ambient water peak, which is attenuated. Simulations suggest that the experiment probes a superposition of two regimes. In the first, fluences expected at the focus severely ionize the water, which becomes effectively transparent to the probe. In the second, out-of-focus pump radiation produces O¹⁺ and O²⁺ ions, which rearrange due to Coulombic repulsion over 10 s of fs. Our simulations account for a decrease in ambient water signal and an increase in low-angle X-ray scattering but not the experimentally observed 9 Å feature, presenting a new challenge for theory.