Empirical PESs

So far, the only approximation in our description of the FMS method has been the use of a finite basis set. When we test for numerical convergence (small model systems and empirical PESs), we often do not make any other approximations but for large systems and/or ab //i/Y/o-determined PESs (AIMS), additional approximations have to be made. These approximations are discussed in this subsection in chronological order (i.e., we begin with the initial basis set and proceed with propagation and analysis of the results). 

Ab initio force fields Ab initio (semi)global PESs Empirical PESs... 

The energy levels determined with the help of empirical PESs cannot match the extreme accuracy of the MARVEL levels but they also provide a complete set. 

Empirical PESs interpolate very well but their extrapolation potential is inferior fo fhose of fhe ab initio PESs. Therefore, even if a highly accurafe empirical PES is available, fhe ab initio surface must be retained as it might prove to be a better choice for finding new fransifions in a new region of fhe specfrum and a better starting point for further refinement of the surface if more defailed experimental information became available. 

The polarizability radial functions obtained in this work can be used for the construction of semi-empirical PES necessary for molecular dynamics studies. 

Figure 3.12. Inelastic scattering of Ar from Pt(lll) at the various input energies listed in the figure and for an initial angle of incidence 0, = 45° and Ts = 800 K. Results are plotted as EfIE vs. the final scattered angle . Points are the experimental results and the lines marked adjacently in the label are results of molecular dynamics simulations on an empirical PES. The long dot-dashed curve is the prediction of a cube model of energy transfer, while the dashed curve is the prediction from hard sphere scattering. From Ref. [135].

Numerous PES s have been proposed for the H + CH4 reaetion. The early ones were empirical form fit to experimental results and neglected the symmetry of the methane molecule.[48] Later, semiempirical and ab initio surfaees, with partially correct permutation symmetry, were developed. Notable among them is the semi-empirical PES by Jordan and Gilbert (JG), symmetric with respect to all four hydrogens on methane.[49] It was built on earlier PESs given by Joseph et a/.[50] and Duchovic et al.[51 ]... 

In part, Simmons s attempt [60] to detail an empirical PES (see Section 8.5.1 and Figure 8.15) for [Cu(bipy)2(0N0)](N03) from variable-temperature crystallographic data on one structure reflects the philosophy of the alternate approach to the one discussed above. A study of metallacyclopentene ring inversion in (s-d5-/7 -buta-... 

Here the bracketed term represents the ensemble average of the mean-sqnare displacement [91Loml]. Clearly the accuracy of the adatom potential energy surface (PES) (Fig. 5) governs the qnality of MD simulations. Both classical (empirical) and semi-empirical PES s can be employed in MD simnlations. 

In empirical PES calculations of interacting molecular species the energy appears directly in this form. 

The "flexibility" of the approximation is an important requirement inq)osed during dynamic calculations to empirical PES. 

We have alluded to the comrection between the molecular PES and the spectroscopic Hamiltonian. These are two very different representations of the molecular Hamiltonian, yet both are supposed to describe the same molecular dynamics. Furthemrore, the PES often is obtained via ab initio quairtum mechanical calculations while the spectroscopic Hamiltonian is most often obtained by an empirical fit to an experimental spectrum. Is there a direct link between these two seemingly very different ways of apprehending the molecular Hamiltonian and dynamics And if so, how consistent are these two distinct ways of viewing the molecule ... 

In the original study (61), NTU /NTU was calculated for thousands of hypothetical design cases as a function of both Pe and Pep -pu- The results were correlated and empirical expressions were given that can be evaluated on a handheld calculator, just as equations 76 and 77, but which approximate the computer calculation much better, to within about 5%. 

A number of empirical tunneling paths have been proposed in order to simplify the two-dimensional problem. Among those are MEP [Kato et al. 1977], sudden straight line [Makri and Miller 1989], and the so-called expectation-value path [Shida et al. 1989]. The results of these papers are hard to compare because slightly different PES were used. As to the expectation-value path, it was constructed as a parametric line q(Q) on which the vibration coordinate q takes its expectation value when Q is fixed. Clearly, for the PES at hand this path coincides with MEP, since is a harmonic oscillator. 

In addition, the frequency cooo, as well as the tunneUng distance can also be extracted from the same empirical data. Thus all the information needed to construct a PES is available. Of course, this PES is a rather crude approximation, since all the skeleton vibrations are replaced by a single mode with effective frequency cooo and coupling parameter C. From the experimental data it is known that the strong hydrogen bond (roo < 2.6 A) is usually typical of intramolecular hydrogen transfer. 

Our studies on the three enzymes have involved the use of semi-empirical methods, using published and also SRP parameter sets. For both LADH and MADH (Figures 5-3a and b) hybrid QM/MM models were employed [8, 9, 88-90], In LADH the PES surface was calculated at the AMI level [20] but scaled by data from the HF/3-21G surface [91]. The results of the CVT calculation with the SCT correction show quite modest yet contributory degrees of tunnelling, an RTE... 

For the discrete bubble model described in Section V.C, future work will be focused on implementation of closure equations in the force balance, like empirical relations for bubble-rise velocities and the interaction between bubbles. Clearly, a more refined model for the bubble-bubble interaction, including coalescence and breakup, is required along with a more realistic description of the rheology of fluidized suspensions. Finally, the adapted model should be augmented with a thermal energy balance, and associated closures for the thermophysical properties, to study heat transport in large-scale fluidized beds, such as FCC-regenerators and PE and PP gas-phase polymerization reactors. 

In classical molecular dynamics, on the other hand, particles move according to the laws of classical mechanics over a PES that has been empirically parameterized. By means of their kinetic energy they can overcome energetic barriers and visit a much more extended portion of phase space. Tools from statistical mechanics can, moreover, be used to determine thermodynamic (e.g. relative free energies) and dynamic properties of the system from its temporal evolution. The quality of the results is, however, limited to the accuracy and reliability of the (empirically) parameterized PES. 

---