Model for the dynamics of a water cluster in an x-ray free electron laser beam

A microscopic sample placed into a focused x-ray free electron laser beam will explode due to strong ionization on a femtosecond time scale. The dynamics of this Coulomb explosion has been modeled by Neutze et al. [Nature (London) 406, 752 (2000)] for a protein, using computer simulations. The results suggest that by using ultrashort exposures, structural information may be collected before the sample is destroyed due to radiation damage. In this paper a method is presented to include the effect of screening by free electrons in the sample in a molecular dynamics simulation. The electrons are approximated by a classical gas, and the electron distribution is calculated iteratively from the Poisson-Boltzmann equation. Test simulations of water clusters reveal the details of the explosion dynamics, as well as the evolution of the free electron gas during the beam exposure. We find that inclusion of the electron gas in the model slows down the Coulomb explosion. The hydrogen atoms leave the sample faster than the oxygen atoms, leading to a double layer of positive ions. A considerable electron density is located between these two layers. The fact that the hydrogens are found to explode much faster than the oxygens means that the diffracting part of the sample stays intact somewhat longer than the sample as a whole.