Using Diffusion Monte Carlo to Evaluate the Initial Conditions for Classical Studies of the Photodissociation Dynamics of HCl Dimer

Efficiency considerations in the construction of interpolated potential energy surfaces for the calculation of quantum observables by diffusion Monte Carlo

A modified Shepard interpolation scheme is used to construct global potential energy surfaces (PES) in order to calculate quantum observables--vibrationally averaged internal coordinates, fully anharmonic zero-point energies and nuclear radial distribution functions--for a prototypical loosely bound molecular system, the water dimer. The efficiency of PES construction is examined with respect to (a) the method used to sample configurational space, (b) the method used to choose which points to add to the PES data set, and (c) the use of either a one- or two-part weight function. The most efficient method for constructing the PES is found to require a quantum sampling regime, a combination of both h-weight and rms methods for choosing data points and use of the two-part weight function in the interpolation. Using this regime, the quantum diffusion Monte Carlo zero-point energy converges to the exact result within addition of 50 data points. The vibrationally averaged O-O distance and O-O radial distribution function, however, converge more slowly and require addition of over 500 data points. The methods presented here are expected to be applicable to both other loosely bound complexes as well as tightly bound molecular species. When combined with high quality ab initio calculations, these methods should be able to accurately characterize the PES of such species.

Photodissociation of hydrogen iodide on the surface of large argon clusters: The orientation of the librational wave function and the scattering from the cluster cage

A set of photodissociation experiments and simulations of hydrogen iodide (HI) on Arn clusters, with an average size n = 139, has been carried out for different laser polarizations. The doped clusters are prepared by a pick-up process. The HI molecule is then photodissociated by a UV laser pulse and the outgoing H fragment is ionized by resonance enhanced multiphoton ionization in a (2 + 1) excitation scheme within the same laser pulse at the wavelength of 243 nm. The measured time-of-flight spectra are transformed into hydrogen kinetic energy distributions. They exhibit a strong fraction of caged H atoms at zero-kinetic energy and peaks at the unperturbed cage exit for both spin-orbit channels nearly independent of the polarization. At this dissociation wavelength, the bare HI molecule exhibits a strict state separation, with a parallel transition to the spin-orbit excited state and perpendicular transitions to the ground state. The experimental results have been reproduced using molecular simulation techniques. Classical molecular dynamics was used to estimate the HI dopant distribution after the pick-up procedure. Subsequently, quasi-classical molecular dynamics (Wigner trajectories approach) has been applied for the photodissociation dynamics. The following main results have been obtained: (i) The HI dopant lands on the surface of the argon cluster during the pick-up process, (ii) zero-point energy plays a dominant role for the hydrogen orientation in the ground state of HI-Arn surface clusters, qualitatively changing the result of the photodissociation experiment upon increasing the number of argon atoms, and, finally, (iii) the scattering of hydrogen atoms from the cage which originate from different dissociation states seriously affects the experimentally measured kinetic energy distributions.