Standard ab initio techniques

Since the basis set does not have to be N-electron because we are solving for each electron separately, we just need atomic basis sets that will allow a linear combination of atomic orbitals (LCAO) for each MO. Naturally the atomic orbitals could consist of hydrogenic functions. For instance, to solve the He problem, the basis set could be the hydrogenic functions (Is, 2s, 2p, 

Gaussian functions are a very practical solution to the basis set problem. They differ from Slater orbitals but a fitted combination of gaussians [contraction, eq. (16)] can be used to simulate a Slater function. 

Gaussian type orbitals or functions are the most widely used in molecular calculations because the integrals Eire relatively simple so that they can be done analytically. They are also used, properly optimized, in density functional calculations. 

The HF method, based on MO-LCAO SCF theory, was a practical tool available for treating molecular systems. Molecular structure determinations were carried out for a large variety of systems, complementing experimental information or providing data that were not available. Many systems that were not known were calculated and their structures were succesfully predicted. However the energetics (bond 

The above correlated methods become very accurate as the order of the correction increases but the cost of the calculation (mainly CPU time) increases as powers of N. Therefore, the most accurate methods are only applicable to very small systems. Formally, the HF method scales as N4, MP2 and CISD as N5, MP3 as N6, MP4 and QCISD(T) as N7, and MP5 as N8. In all cases N is representative of the size of the calculation (for instance, number of basis functions or number of electrons). In practice, all the formal powers decrease as the size of the system increases. 

---

Further uses of pseudopotentials are numerous. The most obvious (and rather widely known ones) are to continue with the PP or MP Hamiltonians for a widely understood combination of the core and valence shells and to apply standard ab initio techniques to electrons in the valence subspace only. We do no elaborate further on this as the hybrid nature of the pseudopotential methods is rather obvious from the above and its more specific applications in a narrower QM/MM hybrid context will be described later. 

There are in general three common ways used today to go about examining the electronic structure of transition metal systems. All of these are concerned in one way or another in solving (or getting around) the correlation problem and we will give an overview of all three, underlining their respective strengths and weaknesses "standard" ab initio techniques (the surest and often... 

The intermolecular force field is calculated by DFT and in some cases by standard ab initio techniques. It is convenient to assign parameters to groups... 

In the earliest implementation applied to molecular problems, K. Johnson [39] used scattered-plane waves as a basis and the exchange-correlation energy was represented by (13). This SW-Xa method employed in addition an (muffin-tin) approximation to the Coulomb potential of (17) in which Vc is replaced by a sum of spherical potentials around each atom. This approximation is well suited for solids for which the SW-Xa method originally was developed [40]. However, it is less appropriate in molecules where the potential around each atom might be far from spherical. The SW-Xa method is computationally expedient compared to standard ab initio techniques and has been used with considerable success [41] to elucidate the electronic structure in complexes and clusters of transition metals. However, the use of the muffin-tin approximation precludes accurate calculations of total energies. The method has for this reason not been successful in studies involving molecular structures and bond energies [42]. 

It must be stressed that the use of GHOs in no way depends on this latter Gaussian expansion procedure since all molecular integrals reduce to standard STO forms. It has merely proved expedient to make use of the expansion method in calculations reported earlier. This ab initio technique provides the raw data with which to establish the patterns of behaviour of the GHOs. We can now address ourselves to one of our stated aims the development of approximation methods for a quantitative theory of valence. 

Several numerical tests and detailed comparisons of MEDLA electron densities to electron densities computed by traditional ab initio SCF technique using 3-21G and 6-31G basis sets have shown [67,71] that the MEDLA results are invariably of better quality than the standard 3-21G ab initio results, and the MEDLA results are virtually indistinguishable from the standard ab initio 6-31G basis set results obtained with the traditional Hartree-Fock method. 

Our own research is developed using molecular orbital calculations based on ab initio techniques these will be the primary focus of this review. Other excellent reviews of this technology with particular relation to experimental structure determination are given by Boggs [78,79] and is covered in standard textbooks [80]. The development of gradient techniques has been essential for the optimization and convergence problems in ab initio calculations [87] and has been carried over to semiempirical calculations as well [60],... 

Several excellent reviews and books have been written on the two main first principles techniques. The standard ab initio, which try to solve directly the Schrodinger equation using a multielectron wavefimction approach, and the DFT methods in which instead of the many-electron wavefimction a non-interactive wavefimction is calculated fi om which the electron density is... 

As far as standard ab initio calculations go, it appears to me that in the future there will be improvements to the computational techniques, though my guess is that at the moment they are ahead of the gas-phase experimental techniques. An example of such a computational improvement [94] is the frequency-dependent moment method in which hyperpolarizabilities are evaluated from configuration-interaction matrices without diagonalizing them. However, it should be remembered that calculations of high-order properties, such as the ones we have been discussing here, are a very tough test of computational quantum chemistry. 

Ab Initio Techniques. The use of ab initio techniques is possible with high resolution neutron and synchrotron data because more non-overlapping reflection data can be obtained than is possible from a standard diffractometer. As the resolving power of these diffractometers further increases, so will the successful use of ab initio calculations. 

This technique was first applied to H2 in a magnetic field by Wille [30], and is similar to the gauge-including (or independent) atomic orbital (GIAO) or London oibital often used in standard ab initio calculations. See references 7-12 and 64. 

Analytic gradient methods became widely used as a result of their implementation for closed-shell self-consistent field (SCF) wavefunctions by Pulay, who has reviewed the development of this topic. Since then, these methods have been extended to deal with all types of SCF wavefunctions, - as well as multi-configuration SCF (MC-SCF), - " configuration-interaction (Cl) wavefunctions, and various non-variational methods such as MoUer-Plesset (MP) perturbation theory - - and coupled-cluster (CC) techniques. - In short, it is possible to obtain analytic energy derivatives for virtually all the standard ab initio approaches. The main use of analytic gradient methods is, and will remain, the location of stationary points on a potential energy siuface, to obtain equilibrium and transition-state geometries. However, there is a specialized use in the calculation of quantities such as dipole derivatives. 

Characterization of amide vibrational modes as seen in IR and Raman spectra has developed from a series of theoretical analyses of empirical data. The designation of amide A, B, I, II, etc., modes stem from several early studies of the (V-methyl acetamide (NMA) molecule vibrational spectra which continues to be a target of theoretical analysis. 15 27,34 162 166,2391 Experimental frequencies were originally fitted to a valence force field using standard vibrational analysis techniques and subsequently were compared to ab initio quantum mechanical force field results. 

Figure 1.12 Determining the properties of drug molecules. Drug molecules may have their properties ascertained by either experimental or theoretical methods. Although experimental methods, especially X-ray crystallography, are the gold standard methods, calculational approaches tend to be faster and do provide high qnality information. Nonempirical techniques, such as ab initio quantum mechanics calcnlations, provide accnrate geometries and electron distribution properties for drng molecnles.

In our short survey of the computational techniques available for investigating TM compounds we first mention molecular mechanics (Chapter 3). It may seem humble by the standards of the quantum mechanical ab initio, semiempirical and DFT methods (Chapters 5, 6 and 8, respectively) but MM is useful for obtaining input structures for submission to one of those calculations, may even provide in itself useful information, and it is, of course, extremely fast. Indeed, a recent book on the modelling of inorganic compounds, mainly TM species, is devoted very largely to molecular mechanics and a program specially parameterized for TM compounds, Momec3 [105],... 

This technique utilizes an option, offered by most ab initio standard programs, to compute the energy of any guess function even if the latter is based on nonorthogonal orbitals. The technique orthogonalizes the orbitals without changing the Slater determinant, and then computes the expectation... 

We illustrate the MOVB method by the SN2 reaction between Cl- and CH3C1, and apply this technique to model substitution reactions. We show that the MOVB method can yield reasonable results for the ground state potential energy surface of the Sn2 reaction both in the gas phase and in solution in comparison with MO and ab initio VB calculations. In all calculations, the standard Gaussian 6-31G(d) basis function is used to construct the MOVB wave function. 

In the previous sections we have discussed the ligand-field theory from the point of view of quantum chemistry, and have presented an ab initio derivation of the ligand-field Hamiltonian (1-5). In principle this Hamiltonian can be constructed explicitly using the standard techniques of computational quantum chemistry, although in practice it is evident that this would be subject to the usual difficulties encountered with large molecules. Our concern in this section is with the use of Eq. (1-5) as the basis for a parameterisation scheme that permits the interpretation of the spectroscopic and magnetic properties of transition metal complexes in terms that are chemically intelligible. 

The first three of these are covered in Chapters 1,3, and 4. Computational methods for the study of molecular energetics, from semi-empirical approaches to sophisticated ab initio methods, although constantly under development and improvement, are now standard techniques and tools for most chemists (Burkert and Allinger 1982 Pople 1999). 

---