Product state and speed distributions in photochemical triple fragmentations.

The clearest dynamical signature of a roaming reaction is a very cold distribution of energy into the rotational and translational degrees of freedom of the roaming donor fragment (e.g. CO) and an exceptionally hot vibrational distribution in the roaming acceptor fragment (e.g. H2, CH4). These signatures were initially identified in joint experimental/theoretical investigations of roaming in H2CO and CH3CHO and are now used to infer the presence of roaming mechanisms in other photodissociation reactions. In this paper we construct a phase space theory (PST) model of triple fragmentation (3F) and show that the dynamical signature of 3F is similar to that of the roaming donor fragment. The PST model starts with a calculation of two-body fragmentation (2F) of a generic molecule, ABC into AB + C. Every AB fragment with sufficient energy to undergo subsequence spontaneous dissociation is allowed to dissociate and the PST distribution of energy into A + B products is calculated for every initial AB state. Using CH3CHO --> HCO + CH3 --> H + CO + CH3 as an example, we calculate that the energy disposal into the rotational and translational degrees of freedom of the 3F products is very low, and is similar to the dynamical signature expected for production of CO via a roaming mechanism. We compare the 3F PST model with published experimental data for photodissociation of CH3CHO and CH3OCHO at energies above the 3F threshold.

Interpolated potential energy surfaces: How accurate do the second derivatives have to be?

A global potential energy surface for the water dimer is constructed using the modified Shepard interpolation scheme of Collins et al. According to this interpolation scheme, the energy at an arbitrary geometry is expressed as a weighted sum of Taylor series expansions from neighboring data points, where the energy and derivative data required are obtained from ab initio calculations. For some ab initio methods, errors are introduced into the second derivative matrix, either by numerical differencing of ab initio energies or numerical integration during the ab initio calculation. Therefore, we test the accuracy required of the second derivative data by truncation of the exact second derivatives to a series of approximate second derivatives, and assess the effect on the results of a quantum diffusion Monte Carlo (QDMC) simulation. Our results show that the calculated zero-point energy and wave function histograms converge to within the numerical uncertainty of the QDMC simulation by inclusion of either three significant figures or three decimal places in the second derivatives.