Theoretical Studies of Anisotropic Melting of Ice Induced by Ultrafast Nonthermal Heating.

Water and ice are routinely studied with X-rays to reveal their diverse structures and anomalous properties. We employ a hybrid collisional-radiative/molecular-dynamics method to explore how femtosecond X-ray pulses interact with hexagonal ice. We find that ice makes a phase transition into a crystalline plasma where its initial structure is maintained up to tens of femtoseconds. The ultrafast melting process occurs anisotropically, where different geometric configurations of the structure melt on different time scales. The transient state and anisotropic melting of crystals can be captured by X-ray diffraction, which impacts any study of crystalline structures probed by femtosecond X-ray lasers.

MolDStruct: Modeling the dynamics and structure of matter exposed to ultrafast x-ray lasers with hybrid collisional-radiative/molecular dynamics.

We describe a method to compute photon-matter interaction and atomic dynamics with x-ray lasers using a hybrid code based on classical molecular dynamics and collisional-radiative calculations. The forces between the atoms are dynamically determined based on changes to their electronic occupations and the formation of a free electron cloud created from the irradiation of photons in the x-ray spectrum. The rapid transition from neutral solid matter to dense plasma phase allows the use of screened potentials, reducing the number of non-bonded interactions. In combination with parallelization through domain decomposition, the hybrid code handles large-scale molecular dynamics and ionization. This method is applicable for large enough samples (solids, liquids, proteins, viruses, atomic clusters, and crystals) that, when exposed to an x-ray laser pulse, turn into a plasma in the first few femtoseconds of the interaction. We present four examples demonstrating the applicability of the method. We investigate the non-thermal heating and scattering of bulk water and damage-induced dynamics of a protein crystal using an x-ray pump-probe scheme. In both cases, we compare to the experimental data. For single particle imaging, we simulate the ultrafast dynamics of a methane cluster exposed to a femtosecond x-ray laser. In the context of coherent diffractive imaging, we study the fragmentation as given by an x-ray pump-probe setup to understand the evolution of radiation damage in the time range of hundreds of femtoseconds.

Coherent diffractive imaging of proteins and viral capsids: simulating MS SPIDOC.

MS SPIDOC is a novel sample delivery system designed for single (isolated) particle imaging at X-ray Free-Electron Lasers that is adaptable towards most large-scale facility beamlines. Biological samples can range from small proteins to MDa particles. Following nano-electrospray ionization, ionic samples can be m/z-filtered and structurally separated before being oriented at the interaction zone. Here, we present the simulation package developed alongside this prototype. The first part describes how the front-to-end ion trajectory simulations have been conducted. Highlighted is a quadrant lens; a simple but efficient device that steers the ion beam within the vicinity of the strong DC orientation field in the interaction zone to ensure spatial overlap with the X-rays. The second part focuses on protein orientation and discusses its potential with respect to diffractive imaging methods. Last, coherent diffractive imaging of prototypical T = 1 and T = 3 norovirus capsids is shown. We use realistic experimental parameters from the SPB/SFX instrument at the European XFEL to demonstrate that low-resolution diffractive imaging data (q < 0.3 nm-1) can be collected with only a few X-ray pulses. Such low-resolution data are sufficient to distinguish between both symmetries of the capsids, allowing to probe low abundant species in a beam if MS SPIDOC is used as sample delivery.

Imaging of femtosecond bond breaking and charge dynamics in ultracharged peptides.

X-ray free-electrons lasers have revolutionized the method of imaging biological macromolecules such as proteins, viruses and cells by opening the door to structural determination of both single particles and crystals at room temperature. By utilizing high intensity X-ray pulses on femtosecond timescales, the effects of radiation damage can be reduced. Achieving high resolution structures will likely require knowledge of how radiation damage affects the structure on an atomic scale, since the experimentally obtained electron densities will be reconstructed in the presence of radiation damage. Detailed understanding of the expected damage scenarios provides further information, in addition to guiding possible corrections that may need to be made to obtain a damage free reconstruction. In this work, we have quantified the effects of ionizing photon-matter interactions using first principles molecular dynamics. We utilize density functional theory to calculate bond breaking and charge dynamics in three ultracharged molecules and two different structural conformations that are important to the structural integrity of biological macromolecules, comparing to our previous studies on amino acids. The effects of the ultracharged states and subsequent bond breaking in real space are studied in reciprocal space using coherent diffractive imaging of an ensemble of aligned biomolecules in the gas phase.

Ultrafast dynamics and scattering of protic ionic liquids induced by XFEL pulses.

X-rays are routinely used for structural studies through scattering, and femtosecond X-ray lasers can probe ultrafast dynamics. We aim to capture the femtosecond dynamics of liquid samples using simulations and deconstruct the interplay of ionization and atomic motion within the X-ray laser pulse. This deconstruction is resolution dependent, as ionization influences the low momentum transfers through changes in scattering form factors, while atomic motion has a greater effect at high momentum transfers through loss of coherence. Our methodology uses a combination of classical molecular dynamics and plasma simulation on a protic ionic liquid to quantify the contributions to the scattering signal and how these evolve with time during the X-ray laser pulse. Our method is relevant for studies of organic liquids, biomolecules in solution or any low-Z materials at liquid densities that quickly turn into a plasma while probed with X-rays.

Structural Heterogeneity in Single Particle Imaging Using X-ray Lasers.

One of the challenges facing single particle imaging with ultrafast X-ray pulses is the structural heterogeneity of the sample to be imaged. For the method to succeed with weakly scattering samples, the diffracted images from a large number of individual proteins need to be averaged. The more the individual proteins differ in structure, the lower the achievable resolution in the final reconstructed image. We use molecular dynamics to simulate two globular proteins in vacuum, fully desolvated as well as with two different solvation layers, at various temperatures. We calculate the diffraction patterns based on the simulations and evaluate the noise in the averaged patterns arising from the structural differences and the surrounding water. Our simulations show that the presence of a minimal water coverage with an average 3 Å thickness will stabilize the protein, reducing the noise associated with structural heterogeneity, whereas additional water will generate more background noise.

Femtosecond bond breaking and charge dynamics in ultracharged amino acids.

Historically, structure determination of nanocrystals, proteins, and macromolecules required the growth of high-quality crystals sufficiently large to diffract X-rays efficiently while withstanding radiation damage. The development of the X-ray free-electron laser has opened the path toward high resolution single particle imaging, and the extreme intensity of the X-rays ensures that enough diffraction statistics are collected before the sample is destroyed by radiation damage. Still, recovery of the structure is a challenge, in part due to the partial fragmentation of the sample during the diffraction event. In this study, we use first-principles based methods to study the impact of radiation induced ionization of six amino acids on the reconstruction process. In particular, we study the fragmentation and charge rearrangement to elucidate the time scales involved and the characteristic fragments occurring.