Wigner distribution trajectory method for collision induced dissociation

Semiclassical Wigner theory of photodissociation in three dimensions: Shedding light on its basis

The semiclassical Wigner theory (SCWT) of photodissociation dynamics, initially proposed by Brown and Heller [J. Chem. Phys. 75, 186 (1981)] in order to describe state distributions in the products of direct collinear photodissociations, was recently extended to realistic three-dimensional triatomic processes of the same type [Arbelo-González et al., Phys. Chem. Chem. Phys. 15, 9994 (2013)]. The resulting approach, which takes into account rotational motions in addition to vibrational and translational ones, was applied to a triatomic-like model of methyl iodide photodissociation and its predictions were found to be in nearly quantitative agreement with rigorous quantum results, but at a much lower computational cost, making thereby SCWT a potential tool for the study of polyatomic reaction dynamics. Here, we analyse the main reasons for this agreement by means of an elementary model of fragmentation explicitly dealing with the rotational motion only. We show that our formulation of SCWT makes it a semiclassical approximation to an approximate planar quantum treatment of the dynamics, both of sufficient quality for the whole treatment to be satisfying.

New insights into the semiclassical Wigner treatment of photodissociation dynamics

The semiclassical Wigner treatment of Brown and Heller [J. Chem. Phys. 1981, 75, 186] is applied to direct triatomic (or triatomic-like polyatomic) photodissociations with the aim of accurately predicting final state distributions at relatively low computational cost, and having available a powerful interpretative tool. For the first time, the treatment takes rotational motions into account. The proposed formulation closely parallels the quantum description as far as possible. An approximate version is proposed, which is still accurate while numerically much more efficient. In addition to being weighted by usual vibrational Wigner distributions, final phase space states appear to be weighted by new rotational Wigner distributions. These densities have remarkable structures clearly showing that classical trajectories most contributing to rotational state j are those reaching the products with a rotational angular momentum close to [j(j + 1)](1/2) (in ℏ units). The previous methods involve running trajectories from the reagent molecule onto the products. The alternative backward approach [L. Bonnet, J. Chem. Phys., 2010, 133, 174108], in which trajectories are run in the reverse direction, is shown to strongly improve the numerical efficiency of the most rigorous method in addition to being state-selective, and thus, ideally suited to the description of state-correlated distributions measured in velocity imaging experiments. The results obtained by means of the previous methods are compared with rigorous quantum results in the case of Guo's triatomic-like model of methyl iodide photodissociation [J. Chem. Phys., 1992, 96, 6629] and close agreement is found. In comparison, the standard method of Goursaud et al. [J. Chem. Phys., 1976, 65, 5453] is only semi-quantitative.

Quasiclassical trajectory study of reactive and dissociative processes in H2+H2: Comparison with quantum-mechanical calculations

Quasiclassical trajectory calculations have been carried out for H(2)(v(1)=high)+H(2)(v(2)=low) collisions within a three degrees of freedom model where five different geometries of the colliding complex were considered. Within this approach, probabilities for different competitive processes are studied: four center reaction, collision induced dissociation, reactive dissociation, and three-body complex formation. The purpose is to compare in detail with equivalent quantum-mechanical wave packet calculations [Bartolomei et al., J. Chem. Phys 122, 064305 (2005)], especially the behavior of the probabilities near reaction thresholds. Quasiclassical calculations compare quite well with the quantum-mechanical ones for collision induced dissociation as well as for the four center reaction, although quantum effects become very important near thresholds, particularly for lower v(1)'s and for the four center process. Less quantitative agreement is found for reactive dissociation and three-body complex formation. It is found that most quantum effects are due to differences between quantum and classical vibrational distributions of H(2)(v(1)=high). Zero point energy violation has been found in the classical reactive-dissociative probabilities. Extension of these findings to full-dimensional treatments is examined.

Comparison of classical and quantum dynamics for collinear cluster scattering

The collinear dynamics of a cluster of four particles colliding with a fixed particle representing a surface is investigated using a four-dimensional wave packet approach. The properties of the system are chosen to resemble a water cluster interacting with graphite, but a deeper surface-particle potential is also considered causing significant dissociation of the cluster. Having four different product arrangement channels the system is quantum mechanically demanding but still manageable. The dynamical richness makes it a suitable benchmark system for evaluation of classical and quantum/classical schemes. The average energy transferred to the cluster and the three dissociation probabilities are presented as function of the initial state of the cluster. In addition to wave packet data, results obtained using quasiclassical as well as Wigner sampled classical trajectories are presented. The main conclusion is that classical mechanics can describe the dynamics of the system in a very satisfactory way. Including zero-point energy in the classical simulations is particularly important for a good description of dissociation but less important for energy transfer.