Use of effective core potentials

The use of Effective Core Potential operators reduces the computational problem in three ways the primitive basis set can be reduced, the contracted basis set can be reduced and the occupied orbital space can be reduced. The reduction of the occupied orbital space is almost inconsequential in molecular calculations, since it neither affects the number of integrals nor the size of the matrices which has to be diagonalized. The reduction of the primitive basis set is of course more important, but since the integral evaluation time is in general not the bottleneck in molecular calculations, this reduction is still of limited importance. There are some cases where the size of the primitive basis set indeed is important, e.g. in direct SCF procedures. The size of the contracted basis set is very important, however. The bottleneck in normal SCF or Cl calculations is the disc storage and/or the iteration time. Both the disc storage and the iteration time depend strongly on the number of contracted functions. 

The process is terminated when the energy converges satisfactorily. Computers presently in use permit a maximum of about 1000 linear parameters and about 10,000 exponential parameters. Up to six particles have been treated explicitly by this method, which can be extended to larger systems by use of effective core potentials. The accuracy of calculated binding energies can be substantially less than 1 meV. 

Calculations are performed for both nickel and cobalt clusters. When the metal atoms are in the atomic d s state (see below) the atoms are treated computationally as one-electron systems by the use of Effective Core Potentials (ECP s). This is a crude but qualitatively correct description and treats the core including the 3d shell bv potentials and projection operators. Similarly for atoms in the d s state a two-electron ECP is used. 

Glukhovtsev MN, Pross A, McGrath MP, Radom L. Extension of Gaussian-2 (G2) theory to bromine-containing and iodine-containing molecules—use of effective core potentials. J Chem Phys 1995 103 1878-1885. 

Hamiltonian have been incorporated in a straightforward manner through the use of effective core potentials. 

The most important modifications of the MNDO method are the use of effective core potentials for the inner orbitals and the inclusion of orthogonalization corrections in a way as was suggested and implemented a long time ago in the SINDOl method [249] at first developed for organic compounds of first-row elements and later extended to the elements of the second and third row [250,251]. 

The results of calculations using effective core potentials of the several types may be compared with experimental measurements, but more useful comparisons can be made with all-electron calculations for the same systems. For example, in studying the use of effective core potentials in QMC calculations, Lao and Christiansen calculated the valence correlation energy for Ne and found excellent agreement with previous full-CI benchmark calculations. They recovered 98-100% of the valence correlation energy and could detect no significant error due to the effective potential approximation. 

The basis sets used for transition metal systems should be of at least double-zeta quality, and commonly larger sets are used it is especially common to have the d functions more flexibly contracted, The use of effective core potentials (ECPs) can reduce the computational effort significantly and allow inclusion of relativistic effects. Experience has shown that if the system contains significant M" " character (where M represents a transition metal atom) it can be important to include the semicore electrons (i.e., the 3s and 3p orbitals for the first transition row) in the valence treatment and not include them in the ECP. While relativistic effects such as the core contraction can be included in the ECP, it is perhaps not surprising that for the third transition row, one must go to treatments that include additional relativistic effects, such as spin-orbit coupling. 

The use of effective core potentials (ECPs) may introduce some errors into computational results. For example, when the ECP fits the DFT functional properly, the associated error may be reduced to 3 kcal mol" in case of 3d transition... 

The intrinsically relativistic nature of the electronic structure of the halogen monoxides is dictated by the presence of two spin-orbit components in the ground electronic state. Any accurate characterization of the XO Xj and X2 rii/2 potentials must therefore treat relativistic effects explicitly. However, diis requirement greatly increases the cost and complexity of the ab initio effort(26), resulting in few relativistic potential surface calculations such as the 10 study by Roszak et al.(21) One more frequently finds relativistic effects incorporated as corrections to non-relativistic energies or treated through the use of effective core potentials (23). 

Hay, P. J., Martin, R. L., 1998, Theoretical Studies of the Structures and Vibrational Frequencies of Actinide Compounds Using Relativistic Effective Core Potentials With Hartree-Fock and Density Functional Methods ... 

Initially, the level of theory that provides accurate geometries and bond energies of TM compounds, yet allows calculations on medium-sized molecules to be performed with reasonable time and CPU resources, had to be determined. Systematic investigations of effective core potentials (ECPs) with different valence basis sets led us to propose a standard level of theory for calculations on TM elements, namely ECPs with valence basis sets of a DZP quality [9, 10]. The small-core ECPs by Hay and Wadt [11] has been chosen, where the original valence basis sets (55/5/N) were decontracted to (441/2111/N-11) withN = 5,4, and 3, for the first-, second-, and third-row TM elements, respectively. The ECPs of the second and third TM rows include scalar relativistic effects while the first-row ECPs are nonrelativistic [11], For main-group elements, either 6-31G(d) [12-16] all electron basis set or, for the heavier elements, ECPs with equivalent (31/31/1) valence basis sets [17] have been employed. This combination has become our standard basis set II, which is used in a majority of our calculations [18]. 

Before any computational study on molecular properties can be carried out, a molecular model needs to be established. It can be based on an appropriate crystal structure or derived using any technique that can produce a valid model for a given compound, whether or not it has been prepared. Molecular mechanics is one such technique and, primarily for reasons of computational simplicity and efficiency, it is one of the most widely used technique. Quantum-mechanical modeling is far more computationally intensive and until recently has been used only rarely for metal complexes. However, the development of effective-core potentials (ECP) and density-functional-theory methods (DFT) has made the use of quantum mechanics a practical alternative. This is particularly so when the electronic structures of a small number of compounds or isomers are required or when transition states or excited states, which are not usually available in molecular mechanics, are to be investigated. However, molecular mechanics is still orders of magnitude faster than ab-initio quantum mechanics and therefore, when large numbers of... 

Figure 1. Calculated geometries of HAIOH, HAINH2, HAIF, HAISH, HBOH, HBeOH, and HMgOH. The three geometries given for HAIOH are obtained using an effective core potential on Al[CP(sp)], an all-electron calculation employing a double- sp basis set [AE(sp)], and an allelectron calculation employing d polarization functions on the Al and O atoms in addition to the...

Two levels of theory are commonly used in the design of the nickel-based catalysts shown in Figure 11 Density Functional Theory (B3LYP functional used with effective core potentials for Ni and 6-3IG for everything else in the complex) and molecular mechanics (both the UFF (4) and reaction force field, RFF (85,86) are used) (87). All these methods are complementary, and the experiments are guided from the results of several calculations using different molecular modeling techniques. 

In Table 6.3, the values of De for RfCU are compared with those obtained within various approximations using relativistic effective core potentials (RECP) Kramers-restricted Hartree-Fock (KRHF) (Han et al 1999), averaged RECP including second-order M0ller-Plesset perturbation theory (AREP-MP2) for the correlation part (Han et al. 1999), RECP coupled-cluster single double (triple) [CCSD(T)] excitations (Han et al. 1999), and a Dirac-Fock-Breit (DFB) method (Malli and Styszynski 1998). The AREP-MP2 calculation of De gives 20.4 eV, while the RECP-CCSD(T) method with correlation leads to 18.8 eV. Our value of De of 19.5 eV is just between these calculated values. 

Rohifing C M and Raghavachari K1990 A theoretical study of small silicon clusters using an effective core potential Chem. Phys. Lett. 167 559... 

---