Parametric potential

The electronic wavefunctions depend — like the potentials — parametrically on all nuclear degrees of freedom (R, r). The stationary wavefunctions (including all electronic and all nuclear degrees of freedom) which represent the initial and the final molecular states for the situation sketched in Figure 15.1, are written in full detail as... 

A common alternative is to synthesize approximate state functions by linear combination of algebraic forms that resemble hydrogenic wave functions. Another strategy is to solve one-particle problems on assuming model potentials parametrically related to molecular size. This approach, known as free-electron simulation, is widely used in solid-state and semiconductor physics. It is the quantum-mechanical extension of the classic (1900) Drude model that pictures a metal as a regular array of cations, immersed in a sea of electrons. Another way to deal with problems of chemical interaction is to describe them as quantum effects, presumably too subtle for the ininitiated to ponder. Two prime examples are, the so-called dispersion interaction that explains van der Waals attraction, and Born repulsion, assumed to occur in ionic crystals. Most chemists are in fact sufficiently intimidated by such claims to consider the problem solved, although not understood. 

Since this approach does not account for long-range electrostatic potentials present in the extended material, the second approach chosen was the rigid-ion lattice energy minimization technique, widely used in solid-state chemistry. This method is based on the use of electrostatic potentials, as well as Born repulsion and bond-bending potentials parametrized such that computed atom—atom distances and angles and other material properties, such as, for instance, elastic constants, are well reproduced for related materials. In our case, parameters were chosen to fit a-quartz. 

Most calculations on CPs have used semiempirical methods. The Hiickel method yields useful results and many properties can be qualitatively understood [188,190], but numbers are not quite reliable. The valence effective Hamiltonian (VEH) method [191] has been applied successfully to CPs [187,192]. It uses atomic potentials parametrized on the results of ab initio HF-SCF calculations on small molecules, and not on experimental data in that sense, it is a purely theoretical method. 

The general implications of the electronic properties of the binding metal center (which are involved in redox potential parametrizations, see below) on the activation of unsaturated small... 

Starting point is QM calculation within the framework of density-functional theory (DFT) (Hohenberg and Kohn, 1964 Kohn and Sham, 1965 Payne et al., 1992). DFT-based energy calculations can be used to evaluate the parameters of classical interatomic interaction potentials, which can be used to perform MS, MC, and MD simulations such ab initio potential parametrization is a key to improving the transferability of the classical force field. In Fig. 1, an interatomic potential energy function for Si-H interactions is given as an example of such a parametrization (Ohira et al., 1995). 

The potential fiinctions for the mteractions between pairs of rare-gas atoms are known to a high degree of accuracy [125]. Flowever, many of them use ad hoc fiinctional fonns parametrized to give the best possible fit to a wide range of experimental data. They will not be considered because it is more instmctive to consider representations that are more finnly rooted in theory and could be used for a wide range of interactions with confidence. 

While simulations reach into larger time spans, the inaccuracies of force fields become more apparent on the one hand properties based on free energies, which were never used for parametrization, are computed more accurately and discrepancies show up on the other hand longer simulations, particularly of proteins, show more subtle discrepancies that only appear after nanoseconds. Thus force fields are under constant revision as far as their parameters are concerned, and this process will continue. Unfortunately the form of the potentials is hardly considered and the refinement leads to an increasing number of distinct atom types with a proliferating number of parameters and a severe detoriation of transferability. The increased use of quantum mechanics to derive potentials will not really improve this situation ab initio quantum mechanics is not reliable enough on the level of kT, and on-the-fly use of quantum methods to derive forces, as in the Car-Parrinello method, is not likely to be applicable to very large systems in the foreseeable future. 

The first point to remark is that methods that are to be incorporated in MD, and thus require frequent updates, must be both accurate and efficient. It is likely that only semi-empirical and density functional (DFT) methods are suitable for embedding. Semi-empirical methods include MO (molecular orbital) [90] and valence-bond methods [89], both being dependent on suitable parametrizations that can be validated by high-level ab initio QM. The quality of DFT has improved recently by refinements of the exchange density functional to such an extent that its accuracy rivals that of the best ab initio calculations [91]. DFT is quite suitable for embedding into a classical environment [92]. Therefore DFT is expected to have the best potential for future incorporation in embedded QM/MD. 

The semi-empirical methods have better MAD s than th Hartree-Fock-based methods, indicating that their parametrization ha accounted for some of the effects of electron correlation. However, thei maximum errors are very large. Semi-empirical methods are especiall poor at predicting ionization potentials and proton affinities. 

The angle bending in H9O occurs without breaking any bonds, and the electron correlation energy is therefore relatively constant over the whole curve. The HF, MP2 and MP4 bending potentials are shown in Figure 11.14, where the reference curve is taken from a parametric fit to a large number of spectroscopic data. ... 

We have used the basis set of the Linear-Muffin-Tin-Orbital (LMTO) method in the atomic sphere approximation (ASA). The LMTO-ASA is based on the work of Andersen and co-workers and the combined technique allows us to treat all phases on equal footing. To treat itinerant magnetism we have employed the Vosko-Wilk-Nusair parametrization for the exchange-correlation energy density and potential. In conjunction with this we have treated the alloying effects for random and partially ordered phases with a multisublattice generalization of the coherent potential approximation (CPA). 

Calculations were done with a full-potential version of the LMTO method with nonoverlapping spheres. The contributions from the interstitial region were accounted for by expanding the products of Hankel functions in a series of atom-ce- -ered Hankels of three different kinetic energies. The corrected tetrahedron method was used for Brillouin zone integration. Electronic exchange and correlation contributions to the total energy were obtained from the local-density functional calculated by Ceperley and Alder " and parametrized by Vosko, Wilk, and Nusair. ... 

Table 2 Elastic constants and bulk moduli for 4d cubic elements. Comparison is made between the results of our tight-binding parametrization (TB), first-principles full potential LAP., results (LAPW), where available, and experiment (Exp.). Calculations were performed at the experimental volume.

---