Potential energy surfaces, fitting

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. 

Lowest minima in the atom—atom potential energy surface, fitted to the ab initio results ... 

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. 

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. 

The same kind of isomerization can also be enforced in the parent compound, namely nitrosohydrogen HNO (150).198 Matrix irradiation of 150 leads to isonitroso hydrogen NOH (151). Calculations of the potential-energy surface of the system HNO/HON show a singlet ground state for HNO 150, which is also the global minimum.199 200 However, for HON 151 a triplet ground state is predicted. Indeed, the experimental IR spectrum of 151 fits much better the calculated spectrum of the triplet. Thus, the isomerization 150 151 involves... 

The various analytic functions and numerical interpolations depending on number of independent variables have been used for fitting the potential energy surfaces. 

Another major, future advance in the quantum chemical computation of potential energy surfaces for reaction dynamics will be the ability to routinely compute the energies of molecular systems on the fly . The tedious and time-consuming process of fitting computed quantum chemical values to functional forms could be avoided if it were possible to compute the PES as needed during a classical trajectory or quantum dynamics calculation. For many chemical reactions, it should be practical in the near future to prudently select a sufficiently rapid and accurate electronic structure method to facilitate dynamics computations on the fly. 

Even a molecularly smooth single-crystal face represents a potential energy surface that depends on the lateral position x, y) of the water molecule in addition to the dependence on the normal distance z. One simple way to introduce this surface corrugation is by adding the lattice periodicity. An example of this approach is given by Berkowitz and co-workers for the interaction between water and the 100 and 111 faces of the Pt crystal. In this case, the full (x, y, z) dependent potential was determined by a fit to the full atomistic model of Heinzinger and co-workers (see later discussion). 

Several researchers have used modified forms of the LEPS potential (see section 1.2 for a discussion of the modifications) to study the dynamics of H2 on the surfaces of various metals. Initial studies of this type were restricted to rigid surfaces, and the parameters in the LEPS surfaces were either determined by fitting to available experimental or theoretical data, or systematically varied to produce potential energy surfaces with specific properties. 

---