Accurate potential energy surfaces

The systematic computational strategy outlined in this section of the review is necessary albeit demanding. The approach provides an accurate description of the entire spectrum of noncovalent interactions between fragments in a cluster. One can be confident in the calculated results regardless of cluster composition [i.e., whether examing the (HiOls, (CgH6)2, or a mixture of the two]. Less obviously but more importantly, one can also be confident in the calculated results across the entire (intermolecular) potential energy 

---

The above discussion represents a necessarily brief simnnary of the aspects of chemical reaction dynamics. The theoretical focus of tliis field is concerned with the development of accurate potential energy surfaces and the calculation of scattering dynamics on these surfaces. Experimentally, much effort has been devoted to developing complementary asymptotic techniques for product characterization and frequency- and time-resolved teclmiques to study transition-state spectroscopy and dynamics. It is instructive to see what can be accomplished with all of these capabilities. Of all the benclunark reactions mentioned in section A3.7.2. the reaction F + H2 —> HE + H represents the best example of how theory and experiment can converge to yield a fairly complete picture of the dynamics of a chemical reaction. Thus, the remainder of this chapter focuses on this reaction as a case study in reaction dynamics. 

Molecular dynamics studies can be done to examine how the path and orientation of approaching reactants lead to a chemical reaction. These studies require an accurate potential energy surface, which is most often an analytic... 

The case of water is particularly convenient because the required high Ka states may be detected in the solar absorption spectrum. However, it is difficult to observe the necessary high vibrational angular momentum states in molecules, which can only be probed by dispersed fluorescence or stimulated emission techniques. On the other hand, it is now possible to perform converged variational calculations on accurate potential energy surfaces, from which one could hope to verify the quantum monodromy and assess the extent to which it is disturbed by perturbations with other modes. Examples of such computed monodromy are seen for H2O in Fig. 2 and LiCN in Fig. 12. 

The Vibrational Energies of Ozone up to the Dissociation Threshold Dynamics Calculations on an Accurate Potential Energy Surface. 

Note that the conventional TST expression is simply the special case of VTST where evaluation is done exclusively for s = 0. As such, the VTST rate constant will always be less than or equal to the conventional TST rate constant (equal in the event that s = 0 minimizes Eq. (15.35)). Put differently, when very accurate potential energy surfaces are available, the conventional TST rate constant is typically an overestimate of the exact classical rate constant. (Note that it is possible, however, for a compensating or even offsetting error to arise from overestimation of the barrier height if the potential energy surface is not very accurate.)... 

Although the underlying approximations are too crude to obtain an accurate potential energy surface, another very important observation can be made when the London equation is compared to the energy expression for H2 the total energy is not equal to the sum of pairwise H-H interactions. Thus, E(Rab, Rac, Rbc) Z Eab + Eac + Ebg, where Eab corresponds to E+ of Eq. (3.31), and Eac and Ebc are given by similar expressions. The simple summation of pairwise H-H interactions only holds for the Coulomb integrals ... 

The harmonic approximation is unrealistic in a dynamical description of the dissociation dynamics, because anharmonic potential energy terms will play an important role in the large amplitude motion associated with dissociation. An accurate potential energy surface must be used in order to obtain a realistic dynamical description of the dissociation process and, as in the quasi-classical approach for bimolecular collisions, a numerical solution of the classical equations of motion is required [2]. 

A promising development in the latter direction was the implementation of a multicon-figurational DFT approach [52] with empirical parameters (in addition to those already contained in the mixed density functionals). This method can yield accurate potential energy surfaces of excited states, and has recently been adapted to perform spin-orbit Cl calculations [53], but the lack of analytic gradients has up to now prevented its use in the simulation of molecular dynamics and in geometry optimizations. A more classically... 

Thiophosgene is one of the simplest and best-studied prototype systems for which accurate potential energy surfaces can be obtained through experiment and high-level ab initio calculations. As such, thiophosgene is a molecule tailor-made for the fundamental understanding of electronic relaxation in polyatomic molecules. The results summarized in this chapter confirm the special role thiophosgene plays in the field of molecular photophysics and photochemistry. 

It seems, therefore, with the current renewal of theoretical interest in atomic and molecular collision problems, reactive scattering, and predissociation phenomena, that it is worthwhile to examine the VB theory as a useful model that is capable of yielding accurate potential energy surfaces. 

We can conclude that there are accurate potential energy surfaces to describe the reaction O ( D) + H2, which plays an important role in the ozone depletion cycle. The most recent PESs correctly reproduce the molecular beam experimental results, namely, the differential cross sections and energy distribution of the products, including the contribution of the abstraction mechanism in the first excited PES, within the present experimental resolution. 

If several electronically excited states are relevant for describing the photodissociation then one or more of the Rydberg orbitals of the molecule must be included in the (CAS) [13], As the number of orbitals and electrons increases in the CAS, the computational time increases dramatically. In order to obtain accurate potential energy surfaces for the excited electronic states, one must include diffuse functions in the basis set [4], For heavier atoms, a relativistic effective core potential (ECP) can be used to treat the scalar relativistic effects. The ECP basis sets have been developed by several research groups [15,16] and have been implemented in most of the standard electronic structure programs. 

The photodissociation dynamics of di- and triatomic molecules are well established. The great difficulty lies in obtaining highly accurate potential energy surfaces and their diabatic couplings. This can be an especially challenging task if many electronic excited states are available in the Franck-Condon region. 

Two methods are in common use for simulating molecular liquids the Monte Carlo method (MC) and molecular dynamics calculations (MD). Both depend on the availability of reasonably accurate potential energy surfaces and both are based on statistical classical mechanics, taking no account of quantum effects. In the past 10-15 years quantum Monte Carlo methods (QMC) have been developed that allow intramolecular degrees of freedom to be studied, but because of the computational complexity of this approach results have only been reported for water clusters. 

Mielke S. L., Lynch G. C., Truhlar D. G. and Schwenke D. W. (1993) A more accurate potential energy surface and quantiun mechanical cross section calculations for the F+ H2 reaction, Chem. Phys. Lett. 213, 10-16 Erratiun (1994) 217, 173. 

Accurate potential energy surfaces are essential if quantitative studies of molecular systems are to be performed. In principle, it is possible to solve equation 1 for the energy of a system to arbitrary accuracy. In practice, however, it proves difficult, if not impossible, to treat a complete system in this way quite simply because the methods that exist to solve... 

J.D. Kress, Z. Bacic Z, G.A Parker, and R.T Pack, Quantum reactive scattering in 3 dimensions using hyperspherical (aph) coordinates. 5. comparison between 2 accurate potential energy surfaces for H + H2 and D + H2. J. Phys. Chem., 94 8055-8058, 1990. 

---