Core-electron integrals

The core-electron integral, V, as in CNDO/INDO, is then equated to the corresponding two-electron integral ... 

It turns out, however, that extreme care must be taken in order to ensure that the surface energy, which is a small energy difference between two much larger total energies is computed reliably. All the issues that can affect the accuracy of the computational set-up, such as basis set, treatment of core electrons, integration... 

The amount of computation for MP2 is determined by the partial tran si ormatioii of the two-electron integrals, what can be done in a time proportionally to m (m is the u umber of basis functions), which IS comparable to computations involved m one step of(iID (doubly-excitcil eon figuration interaction) calculation. fo save some computer time and space, the core orbitals are frequently omitted from MP calculations. For more details on perturbation theory please see A. S/abo and N. Ostlund, Modem Quantum (. hern-isir > Macmillan, Xew York, 198.5. 

The two-center two-eleciroii. one-center two-electron, iwo-center one-electron, one-center one-electron, and core-core repulsion integrals involved in the above equations are discussed below. 

In addition to this term, account must be taken of the decreasing screen in g of then iicleus by th e electron s as the in teratom ic dis-tance becomes very small,. At very small distances the core-core term should approach the classical form. To account for this, an additional term is added to the basic core-core repulsion integral in MlXnO/3 to give ... 

In a closed-shell system, P = P) = P and the Fock matrix elements can be obtained by making this substitution. If a basis set containing s, p orbitals is used, then many of the one-centre integrals nominally included in INDO are equal to zero, as are the core elements Specifically, only the following one-centre, two-electron integrals are non-zero (/x/x /x/x), (pit w) and (fti/lfM/). The elements of the Fock matrix that are affected can then be written a." Uxllow s ... 

III fact, while this correction gives the desired behaviour at relatively long separations, it doLS not account for the fact that as two nuclei approach each other the screening by the core electrons decreases. As the separation approaches zero the core-core repulsion iimild be described by Coulomb s law. In MINDO/3 this is achieved by making the cure-core interaction a function of the electron-electron repulsion integrals as follows ... 

Semiempirical calculations are set up with the same general structure as a HF calculation in that they have a Hamiltonian and a wave function. Within this framework, certain pieces of information are approximated or completely omitted. Usually, the core electrons are not included in the calculation and only a minimal basis set is used. Also, some of the two-electron integrals are omitted. In order to correct for the errors introduced by omitting part of the calculation, the method is parameterized. Parameters to estimate the omitted values are obtained by fitting the results to experimental data or ah initio calculations. Often, these parameters replace some of the integrals that are excluded. 

In order to form the Fock matrix of an ab initio calculation, all the core-Hamiltonian matrix elements, H y, and two-electron integrals (pvIXa) have to be computed. If the total number of basis functions is m, the total number of the core Hamiltonian matrix elements is... 

By de ult, Gaussian will substitute the in-core method for direct SCF when there is enough memory because it is fester. When we ran these computations, we explicitly prevented Gaussian from using the in-core method. When you run your jobs, however, the in-core method will undoubtedly be used for some jobs, and so your values may differ. An in-core job is identified by the following line in the output Two-electron integrals will be kept in memory. 

The most elementary all valence electron NDO model is that known as Ippmplete neglect of differential overlap (CNDO). Segal and Pople introduced (his in 1966. Only valence electrons are explicitly treated, the inner shells being tijicen as part of the atomic core. The ZDO approximation is applied to the WO-electron integrals, so that... 

To construct the Fock matrix, eq. (3.51), integrals over all pairs of basis functions and the one-electron operator h are needed. For M basis functions there are of the order of of such one-electron integrals. These one-integrals are also known as core integrals, they describe the interaction of an electron with the whole frame of bare nuclei. The second part of the Fock matrix involves integrals over four basis functions and the g two-electron operator. There are of the order of of these two-electron integrals. In conventional HF methods the two-electron integrals are calculated and saved before the... 

Form the Fock matrix as the core (one-electron) integrals + the density matrix times the two-electron integrals. 

All calculations were performed on the Cray-2 computers at the Minnesota Supercomputer Center. In some cases the two-electron Integrals could be kept in the 256 megaword central memory of the Cray-2, and in these cases an "in-core" integral and SCF code(53) was used. The largest in-core calculations possible in... 

Thus, the current semi-empirical methods (MNDO, AMI and PM3) differ in the way in which core-repulsions are treated. Within the MNDO formalism the corerepulsion ( asmndo) is expressed in terms of two-centre, two-electron integrals (Eq. 5-4), where Za and Zb correspond to the core charges, Rab is the internuclear separation, and a a and aB are adjustable parameters in the exponential term [19]. 

Despite these modifications there remain a number of well-documented problems with the AM1/PM3 core-repulsion function [37] which has resulted in further refinements. For example, Jorgensen and co-workers have developed the PDDG (pair-wise distance directed Gaussian) PM3 and MNDO methods which display improved accuracy over standard NDDO parameterisations [38], However, for methods which include d-orbitals (e.g. MNDO/d [23,24], AMl/d [25] and AMI [39,40]) it has been found that to obtain the correct balance between attractive and repulsive Coulomb interactions requires an additional adjustable parameter p (previously evaluated using the one-centre two-electron integral Gss, Eq. 5-7), which is used in the evaluation of the two-centre two-electron integrals (Eq. 5-8). 

It is important to point out that recent results on density based overlap integrals [16] confirm the interest of the formulation of Erep as a sum of bond-bond, bond-lone pair and lone pair-lone pair repulsion indeed, core electrons do not contribute to the value of the overlap integrals. 

This integral is that of a core-electron interaction and therefore available through solution of the many-electron wavefunction using a variety of methods. 

Core-electron attraction integrals are calculated using the Goeppert-Mayer-Sklar approximation with neglect of penetration integrals. 

The one-electron core resonance integrals j3. are given by the Mulliken approximation ... 

---