Dissociative ionization of neon clusters Ne(n), n=3 to 14: a realistic multisurface dynamical study

The molecular dynamics with quantum transitions (MDQT) method is applied to study the fragmentation dynamics of neon clusters following vertical ionization of neutral clusters with 3 to 14 atoms. The motion of the neon atoms is treated classically, while transitions between the adiabatic electronic states of the ionic clusters are treated quantum mechanically. The potential energy surfaces are described by the diatomics-in-molecules model in a minimal basis set consisting of the effective 2p orbitals on each neon atom for the missing electron. The fragmentation mechanism is found to be rather explosive, with a large number of events where several atoms simultaneously dissociate. This is in contrast with evaporative atom by atom fragmentation. The dynamics are highly nonadiabatic, especially at shorter times and for the larger clusters. Initial excitation of the neutral clusters does not affect the fragmentation pattern. The influence of spin-orbit coupling is also examined and found to be small, except for the smaller size systems for which the proportion of the Ne+ fragment is increased up to 43%. From the methodological point of view, most of the usual momentum adjustment methods at hopping events are shown to induce nonconservation of the total nuclear angular momentum because of the nonzero electronic to rotation coupling in these systems. A new method for separating out this coupling and enforcing the conservation of the total nuclear momentum is proposed. It is applied here to the MDQT method of Tully but it is very general and can be applied to other surface hopping methods.