The fragmentation upon electron impact ionization of Ar4He1000 is investigated by means of mixed quantum–classical dynamics simulations. The dopant dynamics is described by a surface hopping method coupled with a diatomics-in-molecules model to properly take into account the multiple electronic surfaces and possible transitions between them. Helium atoms are treated individually using zero-point averaged dynamics, a method based on the building of an effective He–He potential. Fast electronic relaxation is observed from less than 2 ps to ∼30 ps, depending on initial conditions. The main fragments observed are and (q ≤ 1000), with a strong contribution of the bare ion, and neither Ar+ nor Ar+Heq fragments are found. The smaller fragments (q ≤ 50) are found to mostly come from ion ejection, whereas larger fragments (q > 500) originate from long-term ion trapping. Although the structure of the trapped ions is the same as in the gas phase, trapped and are rather slightly bound and structures (i.e., an core with one or two argon atoms roaming within the droplet). These loose structures can undergo geminate recombination and release or (q ≤ 50) in the gas phase and/or induce strong helium droplet evaporation. Finally, the translational energy of the fragment center of mass was found to be suitable to provide a clear signature of the broad variety of processes at play in our simulations.
REFERENCES
We call the dynamics “mixed quantum-classical” rather than “classical” since the latter usually refers to a specific method based on the JWKB approximation of the wave function, resulting from a development in powers of Planck’s constant.
The “temperature” is a kinetic temperature since the dynamics is microcanonical. The trajectory is periodically interrupted, and velocities are rescaled if needed in order to maintain the temperature.
An error has been noticed in Ref. 23, p. 5 where tanh[a(R + b)] is indicated instead of tanh[a(R − b)] for the switching function T(R).