Recent advances in electronic structure theory and their influence on the accuracy of ab initio potential energy surfaces

Recent advances in electronic structure theory and the availability of high speed vector processors have substantially increased the accuracy of ab initio potential energy surfaces. The recently developed atomic natural orbital approach for basis set contraction has reduced both the basis set incompleteness and superposition errors in molecular calculations. Furthermore, full CI calculations can often be used to calibrate a CASSCF/MRCI approach that quantitatively accounts for the valence correlation energy. These computational advances also provide a vehicle for systematically improving the calculations and for estimating the residual error in the calculations. Calculations on selected diatomic and triatomic systems will be used to illustrate the accuracy that currently can be achieved for molecular systems. In particular, the F+H2 yields HF+H potential energy hypersurface is used to illustrate the impact of these computational advances on the calculation of potential energy surfaces

A complete quantum mechanical study of chlorine photodissociation.

A fully quantum mechanical dynamical calculation on the photodissociation of molecular chlorine is presented. The magnitudes and phases of all the relevant photofragment T-matrices have been calculated, making this study the computational equivalent of a "complete experiment," where all the possible parameters defining an experiment have been determined. The results are used to simulate cross-sections and angular momentum polarization information which may be compared with experimental data. The calculations rigorously confirm the currently accepted mechanism for the UV photodissociation of Cl(2), in which the majority of the products exit on the C(1)Π(1u) state, with non-adiabatic couplings to the A(3)Π(1u) and several other Ω = 1 states, and a small contribution from the B(3)Π state present at longer wavelengths

The vibrationally mediated photodissociation of Cl2.

The photodissociation of vibrationally excited Cl(2)(v = 1) has been investigated experimentally using the velocity mapped ion imaging technique. The experimental measurements presented here are compared with the results of time-dependent wavepacket calculations performed on a set of ab initio potential energy curves. The high level calculations allow prediction of all the dynamical information regarding the dissociation, including electronic polarization effects. Using a combination of theory and experiment it was found that there was negligible cooling of the vibrational degree of freedom of the parent molecule in the molecular beam. The results presented are compared with those following the photodissociation of Cl(2)(v = 0). Although the same electronic states are found to be important for Cl(2)(v = 1) as for Cl(2)(v = 0), significant differences were found regarding many of the observables. The overall level of agreement between theory and experiment was found to be reasonable and confirms previous assignments of the photodissociation mechanism

Electronic polarization effects in the photodissociation of Cl2.

Velocity mapped ion imaging and resonantly enhanced multiphoton ionization time-of-flight methods have been used to investigate the photodissociation dynamics of the diatomic molecule Cl(2) following excitation to the first UV absorption band. The experimental results presented here are compared with high level time dependent wavepacket calculations performed on a set of ab initio potential energy curves [D. B. Kokh, A. B. Alekseyev, and R. J. Buenker, J. Chem. Phys. 120, 11549 (2004)]. The theoretical calculations provide the first determination of all dynamical information regarding the dissociation of a system of this complexity, including angular momentum polarization. Both low rank K = 1, 2 and high rank K = 3 electronic polarization are predicted to be important for dissociation into both asymptotic product channels and, in general, good agreement is found between the recent theory and the measurements made here, which include the first experimental determination of high rank K = 3 orientation

Accurate DMBE Potential Energy Surface For the N(2D) + H2(1S g+ ) Reaction Using an Improved Switching Function Formalism

Abstract A single-sheeted double many-body expansion (DMBE) potential energy surface is reported for the 1 2 A'' state of NH2. To approximate its true multi-sheeted nature, a novel switching function that imposes the correct behavior at the H2(X 1S g +)+ N(2 D) and NH(X 3S-) + H(2 S) dissociation limits has been suggested. The new DMBE form is shown to fit with high accuracy an extensive set of new ab initio points (calculated at the multi-reference configuration interaction level using the full valence complete active space as reference and aug-cc-pVQZ and aug-cc-pV5Z basis sets) that have been semiempirically corrected at the valence regions by scaling the n-body dynamical correlation terms such as to account for the finite basis set size and truncated configuration interaction expansion. A detailed study of the N(2 D) ... H2(X 1S g +) van der Waals region has also been carried out. These calculations predict a nearly free rigid-rotor with two shallow van der Waals wells of C 2v and C 8 v symmetries. Such a result contrasts with previous cc-pVTZ calculations which predict a single T-shaped van der Waals structure. Except in the vicinity of the crossing seam, which is replaced by an avoided intersection, the fit shows the correct physical behavior over the entire configurational space. The topographical features of the new DMBE potential energy surface are examined in detail and compared with those of other potential functions available in the literature. Amongst such features, we highlight the barrier for linearization (11,802 cm-1) which is found to overestimate the most recent empirical spectroscopic estimate by only 28 cm-1. Additionally, the T-shaped N(2 D) ... H2 van der Waals minimum is predicted to have a well depth of 90 cm-1, being 11 cm-1 deeper than the C 8 v minimum. The title DMBE form is therefore recommendable for dynamics studies of both non-reactive and reactive N(2 D)+H2 collisions