Electron repulsion functional

A recent account of the thermodynamic analogy of the present derivatives given in eqns (61), (63), (64) and (65) can be found in ref. 72. In these equations, E is the total energy of the system and F the Hohenberg-Kohn functional, the sum of the kinetic energy and electron-electron repulsion functional. 

VV e now wish to establish the general functional form of possible wavefunctions for the two electrons in this pseudo helium atom. We will do so by considering first the spatial part of the u a efunction. We will show how to derive functional forms for the wavefunction in which the i change of electrons is independent of the electron labels and does not affect the electron density. The simplest approach is to assume that each wavefunction for the helium atom is the product of the individual one-electron solutions. As we have just seen, this implies that the total energy is equal to the sum of the one-electron orbital energies, which is not correct as ii ignores electron-electron repulsion. Nevertheless, it is a useful illustrative model. The wavefunction of the lowest energy state then has each of the two electrons in a Is orbital ... 

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 ... 

Set this threshold to a small positive constant (the default value is 10" ° Hartree). This threshold is used by HyperChem to ignore all two-electron repulsion integrals with an absolute value less than this value. This option controls the performance of the SCF iterations and the accuracy of the wave function and energies since it can decrease the number of calculated two-electron integrals. 

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 uses 16 bytes (two double-precision words) of storage for each electron repulsion integral. The first 8 bytes save the compressed four indices and the second 8 bytes store the value of the integral. Each index takes 16 bits. Thus the maximum number of basis functions is 65,535. This should satisfy all users of HyperChem for the foreseeable future. 

From electronic structure theory it is known that the repulsion is due to overlap of the electronic wave functions, and furthermore that the electron density falls off approximately exponentially with the distance from the nucleus (the exact wave function for the hydrogen atom is an exponential function). There is therefore some justification for choosing the repulsive part as an exponential function. The general form of the Exponential - R Ey w function, also known as a ""Buckingham " or ""Hill" type potential is... 

The electron-electron repulsion operator has a singularity for r 12 = 0 which results in the exact wave function having a cusp (discontinuous derivative).- ... 

The correlation of electron motion in molecular systems is responsible for many important effects, but its theoretical treatment has proved to be very difficult. Thus many quantum valence calculations use wave functions which are adjusted to optimize kinetic energy effects and the potential energy of interaction of nuclei and electrons but which do not adequately allow for electron correlation and hence yield excessive electron repulsion energy. This problem may be subdivided into cases of overlapping and nonoverlapping electron distributions. Both are very important but we shall concern ourselves here with only the nonoverlapping case. 

Electron-electron repulsion integrals, 28 Electrons bonding, 14, 18-19 electron-electron repulsion, 8 inner-shell core, 4 ionization energy of, 10 localization of, 16 polarization of, 75 Schroedinger equation for, 2 triplet spin states, 15-16 valence, core-valence separation, 4 wave functions of, 4,15-16 Electrostatic fields, of proteins, 122 Electrostatic interactions, 13, 87 in enzymatic reactions, 209-211,225-228 in lysozyme, 158-161,167-169 in metalloenzymes, 200-207 in proteins ... 

Table 1.2 The excess inter-electronic repulsion for T configurations (column 4), relative to a smoothly varying baseline function, g(n), drawn through the formulae for f, f and f. ...

The existence of the first HK theorem is quite surprising since electron-electron repulsion is a two-electron phenomenon and the electron density depends only on one set of electronic coordinates. Unfortunately, the universal functional is unknown and a plethora of different forms have been suggested that have been inspired by model systems such as the uniform or weakly inhomogeneous electron gas, the helium atom, or simply in an ad hoc way. A recent review describes the major classes of presently used density functionals [10]. 

Consider a crude approximation to the ground state of the lithium atom in which the electron-electron repulsions are neglected. Construct the ground-state wave function in terms of the hydrogen-like atomic orbitals. 

The solution of the Schrodinger equation with Helec is the electronic wave function xTelec and the electronic energy Eelec. depends on the electron coordinates, while the nuclear coordinates enter only parametrically and do not explicitly appear in Pelec. The total energy Etot is then the sum of Eelec and the constant nuclear repulsion term, M M... 

An important consequence of the only approximate treatment of the electron-electron repulsion is that the true wave function of a many electron system is never a single Slater determinant We may ask now if SD is not the exact wave function of N interacting electrons, is there any other (necessarily artificial model) system of which it is the correct wave function The answer is Yes it can easily be shown that a Slater determinant is indeed an eigenfunction of a Hamilton operator defined as the sum of the Fock operators of equation (1-25)... 

The energy due to the external potential is determined simply by the density and is therefore independent of the wave function generating that density. Hence, it is the same for all wave functions integrating to a particular density and we can separate it from the kinetic and electron-electron repulsion contributions... 

To understand how Kohn and Sham tackled this problem, we go back to the discussion of the Hartree-Fock scheme in Chapter 1. There, our wave function was a single Slater determinant SD constructed from N spin orbitals. While the Slater determinant enters the HF method as the approximation to the true N-electron wave function, we showed in Section 1.3 that 4>sd can also be looked upon as the exact wave function of a fictitious system of N non-interacting electrons (that is electrons which behave as uncharged fermions and therefore do not interact with each other via Coulomb repulsion), moving in the elfective potential VHF. For this type of wave function the kinetic energy can be exactly expressed as... 

Redress can be obtained by the electron localization function (ELF). It decomposes the electron density spatially into regions that correspond to the notion of electron pairs, and its results are compatible with the valence shell electron-pair repulsion theory. An electron has a certain electron density p, (x, y, z) at a site x, y, z this can be calculated with quantum mechanics. Take a small, spherical volume element AV around this site. The product nY(x, y, z) = p, (x, y, z)AV corresponds to the number of electrons in this volume element. For a given number of electrons the size of the sphere AV adapts itself to the electron density. For this given number of electrons one can calculate the probability w(x, y, z) of finding a second electron with the same spin within this very volume element. According to the Pauli principle this electron must belong to another electron pair. The electron localization function is defined with the aid of this probability ... 

The first two tenns in Eq. (2) represent the kinetic energy of the nuclei and the electrons, respectively. The remaining three terms specify the potential energy as a function of the interaction between the particles. Equation (3) expresses the potential function for the interaction of each pair of nuclei. In general, this sum is composed of terms that are given by Coulomb s law for the repulsion between particles of like charge. Similarly, Eq. (4) corresponds to the electron-electron repulsion. Finally, Eq. (5) is the potential function for the attraction between a given electron (<) and a nucleus (j). 

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). 

---