Semi functions

Functional fonns based on the above ideas are used in the FIFD [127] and Tang-Toeimies models [129], where the repulsion tenn is obtained by fitting to Flartree-Fock calculations, and in the XC model [92] where the repulsion is modelled by an ab initio Coulomb tenn and a semi-empirical exchange-repulsion tenn Cunent versions of all these models employ an individually damped dispersion series for the attractive... 

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. 

DFT calculations offer a good compromise between speed and accuracy. They are well suited for problem molecules such as transition metal complexes. This feature has revolutionized computational inorganic chemistry. DFT often underestimates activation energies and many functionals reproduce hydrogen bonds poorly. Weak van der Waals interactions (dispersion) are not reproduced by DFT a weakness that is shared with current semi-empirical MO techniques. 

I hcre arc two types of Cl calculations im piemen ted in Hyper-Ch ern sin gly exciled Cl an d in icroslate Cl. I hc sin gly excited C which is available for both ah initio and sem i-etn pirical calculations may be used to generate CV spectra and the microstate Cl available only for the semi-empirical methods in HyperChern is used to improve the wave function and energies including the electron ic correlation. On ly sin gle point calculation s can he perform cd in HyperChetn using Cl. 

TTyperChem s semi-einpirical calculations solve (approximately) the Schriidinger equation for this electronic Tlainiltonian leading to an electronic wave function ] j,, , (r) for the electrons ... 

In ub initio calculations all elements of the Fock matrix are calculated using Equation (2.226), ii re peifive of whether the basis functions ip, cp, formally bonded. To discuss the semi-empirical melh ids it is useful to consider the Fock matrix elements in three groups (the diagonal... 

Chapter 2 we worked through the two most commonly used quantum mechanical models r performing calculations on ground-state organic -like molecules, the ab initio and semi-ipirical approaches. We also considered some of the properties that can be calculated ing these techniques. In this chapter we will consider various advanced features of the ab Itio approach and also examine the use of density functional methods. Finally, we will amine the important topic of how quantum mechanics can be used to study the solid state. 

HyperChem can plot orbital wave functions resulting from semi-empirical and ab initio quantum mechanical calculations. It is interesting to view both the nodal properties and the relative sizes of the wave functions. Orbital wave functions can provide chemical insights. 

Parameters for elements (basis functions in ab initio methods usually derived from experimental data and empirical parameters in semi-empirical methods usually obtained from empirical data or ab initio calculations) are independent of the chemical environment. In contrast, parameters used in molecular mechanics methods often depend on the chemical environment. 

The simplest many-electron wave function that satisfies the Exclusion Principle is a product of N different one-electron functions that have been antisymmetrized, or written as a determinant. Here, N is the number of electrons (or valence electrons) in the molecule. HyperChem uses this form of the wave function for most semi-empirical and ab initio calculations. Exceptions involve using the Configuration Interaction option (see page 119). HyperChem computes one-electron functions, termed molecular spin orbitals, by relatively simple integration and summation calculations. The many-electron wave function, which has N terms (the number of terms in the determinant), never needs to be evaluated. 

You can extend the calculation of the Hartree-Eock semi-empirical wave function by choosing Configuration Interaction (Cl) in the... 

Semi-empirical Options dialog box and you can also extend the calculation of the Hartree-Fock ab initio wave function by choosing Cl in the Ab Initio Options dialog box. Use Cl for these electron configurations ... 

This equation is important in interpreting the results of calculations. In ab initio and semi-empirical calculations, atomic orbitals are functions of the x, y, and z coordinates of the electron that closely resemble the valence orbitals of the isolated atoms. 

The setup of these calculations is very similar for both quantum and molecular mechanics. In practice, molecular dynamics calculation s using the nl) initio and semi-empirical quantum mechanical SCFmethods are limited to relatively small systems. Each time step requires a complete calculation of the wave function and the forces. 

In order to conserve the total energy in molecular dynamics calculations using semi-empirical methods, the gradient needs to be very accurate. Although the gradient is calculated analytically, it is a function of wavefunction, so its accuracy depends on that of the wavefunction. Tests for CH4 show that the convergence limit needs to be at most le-6 for CNDO and INDO and le-7 for MINDO/3, MNDO, AMI, and PM3 for accurate energy conservation. ZINDO/S is not suitable for molecular dynamics calculations. 

Once you have calculated an ab initio or a semi-empirical wave function via a single point calculation, geometry optimization, molecular dynamics or vibrations, you can plot the electrostatic potential surrounding the molecule, the total electronic density, the spin density, one or more molecular orbitals /i, and the electron densities of individual orbitals You can examine orbital energies and select orbitals for plotting from an orbital energy level diagram. 

Many problems with MNDO involve cases where the NDO approximation electron-electron repulsion is most important. AMI is an improvement over MNDO, even though it uses the same basic approximation. It is generally the most accurate semi-empirical method in HyperChem and is the method of choice for most problems. Altering part of the theoretical framework (the function describing repulsion between atomic cores) and assigning new parameters improves the performance of AMI. It deals with hydrogen bonds properly, produces accurate predictions of activation barriers for many reactions, and predicts heats of formation of molecules with an error that is about 40 percent smaller than with MNDO. 

HyperChem quantum mechanical calculations are ab initio and semi-empirical. Ab initio calculations use parameters (contracted basis functions) associated with shells, such as an s shell, sp shell, etc., or atomic numbers (atoms). Semi-empirical calculations use parameters associated with specific atomic numbers. The concept of atom types is not used in the conventional quantum mechanics methods. Semi-empirical quantum mechanics methods use a rigorous quantum mechanical formulation combined with the use of empirical parameters obtained from comparison with experiment. If parameters are available for the atoms of a given molecule, the ab initio and semi-empirical calculations have an a priori aspect when compared with a molecular mechanics calculation, letting... 

When you perform a single point semi-empirical or ab initio calculation, you obtain the energy and the first derivatives of the energy with respect to Cartesian displacement of the atoms. Since the wave function for the molecule is computed in the process, there are a number of other molecular properties that could be available to you. Molecular properties are basically an average over the wave function of certain operators describing the property. For example, the electronic dipole operator is basically just the operator for the position of an electron and the electronic contribution to the dipole moment is... 

If a molecule is rotated by changing the position of the viewer (left mouse button rotation) then the molecule s position in the molecular coordinate system has not changed and another contour plot can be requested without recomputing the wave function. That is, many orbitals can be plotted after a single point ab initio or semi-empirical calculation. Any contour map is available without recomputation of the wave function. 

A recent victim of the killer micros was Evans and Sutherland s parallel computer development effort, halted ia 1990. Their architecture combiaed a small number of approximately 1-MFLOPS processors iato semi-iadependent functional units. Several of these units could, ia turn, be combiaed to form a processor hierarchy, building up to systems that were expected to cost between 1 and 8 million dollars. With the advent of lO-MFLOPS uniprocessor killer micros, such an architecture became irrelevant and the project was halted. The RISC killer micro could deUver the same level of performance as could the combiaed efforts of 10 of the 1-MFLOPS processors, evea with the unlikely assumptioa that the problem could be perfectiy parallelized across 10 processors. 

---