Fitted potential energy surface

The computed and fitted potential energy surfaces and their symmetry and resonance assignment are shown in Fig. 3a. Figure 3b shows the parameter geometry dependence that accomplishes the fit of the singlet states. Their semiquantitative relevance is validated by the same reasons discussed in the previous section. 

Spline functions have been used successfully to fit potential energy surfaces in one and two dimensions [52,53]. However for higher dimensional surfaces im-physical oscillations [53] can easily occur, and densely distributed ab initio points maybe required [53,54]. 

G.C. Maisuradze, et al.. Interpolating moving least-squares methods for fitting potential energy surfaces Detailed analysis of one-dimensional applications, /. Chem. Phys. 119 (19) (2003) 10002-10014. 

H.G. Yu, S. Andersson, G. Nyman, A generalized discrete variable representation approach to interpolating or fitting potential energy surfaces, Chem. Phys. Lett. 321 (3-4) (2000) 275-280. 

In this section, we will discuss how we have used CVPT to follow the arrows up through Fig. 1—that is, how we have used CVPT to fit potential energy surfaces using CVPT. After presenting the general approach, we will discuss potentials for HCN and... 

There are three reasons for this improvement. First the 12 state was misassigned (20). Second, much of the earlier work used an approximation to the force field in order to simplify the calculation of the vibrational states. In particular, the internal coordinate force field was reexpressed as a low-order expansion in terms of the rectilinear normal coordinates (96,74). The final reason for the improvement is a result of the quality of the fit potential energy surface which we have obtained. The earlier adjusted force fields were obtained by changing only the diagonal quadratic force constants. 

Hu C H and Thakkar A J 1996 Potential energy surface for interactions between N2 and He ab initio calculations, analytic fits, and second virial coefficients J. Chem. Phys. 104 2541... 

The fitting parameters in the transfomi method are properties related to the two potential energy surfaces that define die electronic resonance. These curves are obtained when the two hypersurfaces are cut along theyth nomial mode coordinate. In order of increasing theoretical sophistication these properties are (i) the relative position of their minima (often called the displacement parameters), (ii) the force constant of the vibration (its frequency), (iii) nuclear coordinate dependence of the electronic transition moment and (iv) the issue of mode mixing upon excitation—known as the Duschinsky effect—requiring a multidimensional approach. 

The Ar-HCl and Ar-HF complexes became prototypes for the study of intennolecular forces. Holmgren et al [30] produced an empirical potential energy surface for Ar-HCl fitted to the microwave and radiofrequency spectra,... 

The full dynamical treatment of electrons and nuclei together in a laboratory system of coordinates is computationally intensive and difficult. However, the availability of multiprocessor computers and detailed attention to the development of efficient software, such as ENDyne, which can be maintained and debugged continually when new features are added, make END a viable alternative among methods for the study of molecular processes. Eurthemiore, when the application of END is compared to the total effort of accurate determination of relevant potential energy surfaces and nonadiabatic coupling terms, faithful analytical fitting and interpolation of the common pointwise representation of surfaces and coupling terms, and the solution of the coupled dynamical equations in a suitable internal coordinates, the computational effort of END is competitive. 

In Chapter IX, Liang et al. present an approach, termed as the crude Bom-Oppenheimer approximation, which is based on the Born-Oppen-heimer approximation but employs the straightforward perturbation method. Within their chapter they develop this approximation to become a practical method for computing potential energy surfaces. They show that to carry out different orders of perturbation, the ability to calculate the matrix elements of the derivatives of the Coulomb interaction with respect to nuclear coordinates is essential. For this purpose, they study a diatomic molecule, and by doing that demonstrate the basic skill to compute the relevant matrix elements for the Gaussian basis sets. Finally, they apply this approach to the H2 molecule and show that the calculated equilibrium position and foree constant fit reasonable well those obtained by other approaches. 

Molecular mechanics methods are not generally applicable to structures very far from equilibrium, such as transition structures. Calculations that use algebraic expressions to describe the reaction path and transition structure are usually semiclassical algorithms. These calculations use an energy expression fitted to an ah initio potential energy surface for that exact reaction, rather than using the same parameters for every molecule. Semiclassical calculations are discussed further in Chapter 19. 

Ah initio trajectory calculations have now been performed. However, these calculations require such an enormous amount of computer time that they have only been done on the simplest systems. At the present time, these calculations are too expensive to be used for computing rate constants, which require many trajectories to be computed. Semiempirical methods have been designed specifically for dynamics calculations, which have given insight into vibrational motion, but they have not been the methods of choice for computing rate constants since they are generally inferior to analytic potential energy surfaces fitted from ah initio results. 

At first sight, the easiest approach is to fit a set of points near the saddle point to some analytical expression. Derivatives of the fitted function can then be used to locate the saddle point. This method has been well used for small molecules (see Sana, 1981). An accurate fit to a large portion of the potential energy surface is also needed for the study of reaction dynamics by classical or semi-classical trajectory methods. 

---