Interpolation of potential

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. 

The method which has been developed to construct a potential energy surface (PES) by interpolation of ab initio data has been reviewed in some... 

Certain transition metal ions such as Co2, Ti3 are known to form chelates with trimethylphosphate, i.e., dimethoxyphosphato complexes84 8S>. The chelate effect is responsible for the high stabilities of such complexes, which is expressed in the more negative values for the half-wave potentials. All ions producing such complexes are expected to undergo reduction in TMP at more negative potentials than would be expected from interpolation of the curves. 

Figure 4. Selected iso-n and and iso-h contours superimposed on the potential energy surface for maltose. The iso-energy contours are drawn by interpolation of 1 kcal/mol with respect to the relative energy minimum ( ). The iso-h = 0 contour divides the map into two regions, corresponding to right-handed (h>0) and left-handed (h<0) chirality.

Some of the potential energy functions used to calculate the total strain energy of a molecule are similar to the functions used in the analysis of vibrational spectra. Because the parameters used to derive the strain energies from these functions are fitted quantities, which are based on experimental data (for example X-ray structures), molecular mechanics may be referred to as empirical force field calculations (more often the simplification force field calculations is used). The quality of such calculations is strongly dependent on the reliability of potential energy functions and the corresponding parameters (the force field). Thus, the selection of experimental data to fit the force field is one of the most important steps in a molecular mechanics study. An empirical force field calculation is in essence a method where the structure and the strain energy of an unknown molecule are interpolated from a series of similar molecules with known structures and properties. 

Before a detailed presentation of the ab initio dynamics simulations, first the fundamental difference between atomic and molecular adsorption on the one hand and dissociative adsorption on the other hand has to be addressed. Then I will briefly discuss the question whether quantum or classical methods are appropriate for the simulation of the adsorption dynamics. This section will be followed by a short introduction into the determination of potential energy surfaces from first principles and their continuous representation by some analytical or numerical interpolation schemes. Then the dissociative adsorption and associative desorption of hydrogen at metal and semiconductor surfaces and the molecular trapping of oxygen on platinum will be discussed in some detail. 

An alternative approach is the interpolation of the ab initio PES by some suitable analytical or numerical scheme. For the six-dimensional quantum dynamical studies of hydrogen dissociation on Pd(l 0 0) and Cu(l 0 0) discussed in the next section, ab initio potential energy surfaces have been fitted to an analytical representations [5, 10, 13, 15, 38]. 

CVT approach is particularly attractive due to the limited amount of potential energy and Hessian information that is required to perform the calculations. Direct dynamics with CVT thus offers an efficient and cost-effective methodology. Furthermore, several theoretical reviews60,61 have indicated that CVT plus multidimensional semi-classical tunneling approximations yield accurate rate constants not only for gas-phase reactions but also for chemisorption and diffusion on metals. Computationally, it is expensive if these Hessians are to be calculated at an accurate level of ab initio molecular orbital theory. Several approaches have been proposed to reduce this computational demand. One approach is to estimate rate constants and tunneling contributions by using Interpolated CVT when the available accurate ab initio electronic structure information is very limited.62 Another way is to carry out CVT calculations with multidimensional semi-classical tunneling approximations. 

By means of the Pohl commutator P the potentiometer was now connected alternately to the potential leads, p and p, coming from the ends of the cathode wire, and to those of the ten-ohm comparison resistance, readings of the fall of potential across the two resistances being taken several times in rapid succession, while the time was in each case noted, in order that the value for the ten-ohm might be calculated by interpolation for the exact moment at which that for the cathode wire was observed. This made it possible to eliminate the error which would otherwise have resulted from the gradual drift in the resistance of the cathode, and hence in the measuring current. 

In the table below, column 1 shows the measuring current in tnilli-amperes column 2 the time elapsed from the first reading in seconds and column 3 the interval between successive observations. The unenclosed numbers in the column headed Ei0 give the observed fall of potential across the ten-ohm resistance, and the numbers in parentheses are interpolated from the preceding for the time at which Ex was read. Ex is the observed fall of potential across the cathode wire and Rx is the resistance of the latter calculated from the relation Rx = 10mEx/(Eio). 

---