Hartree problem

Thus the Hartree problem is much more difficult to solve than the NICR... 

A highly readable account of early efforts to apply the independent-particle approximation to problems of organic chemistry. Although more accurate computational methods have since been developed for treating all of the problems discussed in the text, its discussion of approximate Hartree-Fock (semiempirical) methods and their accuracy is still useful. Moreover, the view supplied about what was understood and what was not understood in physical organic chemistry three decades ago is... 

Direct dynamics attempts to break this bottleneck in the study of MD, retaining the accuracy of the full electronic PES without the need for an analytic fit of data. The first studies in this field used semiclassical methods with semiempirical [66,67] or simple Hartree-Fock [68] wave functions to heat the electrons. These first studies used what is called BO dynamics, evaluating the PES at each step from the elech onic wave function obtained by solution of the electronic structure problem. An alternative, the Ehrenfest dynamics method, is to propagate the electronic wave function at the same time as the nuclei. Although early direct dynamics studies using this method [69-71] restricted themselves to adiabatic problems, the method can incorporate non-adiabatic effects directly in the electionic wave function. 

The problem with most quantum mechanical methods is that they scale badly. This means that, for instance, a calculation for twice as large a molecule does not require twice as much computer time and resources (this would be linear scaling), but rather 2" times as much, where n varies between about 3 for DFT calculations to 4 for Hartree-Fock and very large numbers for ab-initio techniques with explicit treatment of electron correlation. Thus, the size of the molecules that we can treat with conventional methods is limited. Linear scaling methods have been developed for ab-initio, DFT and semi-empirical methods, but only the latter are currently able to treat complete enzymes. There are two different approaches available. 

Within the periodic Hartree-Fock approach it is possible to incorporate many of the variants that we have discussed, such as LFHF or RHF. Density functional theory can also be used. I his makes it possible to compare the results obtained from these variants. Whilst density functional theory is more widely used for solid-state applications, there are certain types of problem that are currently more amenable to the Hartree-Fock method. Of particular ii. Icvance here are systems containing unpaired electrons, two recent examples being the clci tronic and magnetic properties of nickel oxide and alkaline earth oxides doped with alkali metal ions (Li in CaO) [Dovesi et al. 2000]. 

III. The Unrestricted Hartree-Fock Spin Impurity Problem... 

It should be noted that by moving to a matrix problem, one does not remove the need for an iterative solution the Fj y matrix elements depend on the Cy i LCAO-MO eoeffieients whieh are, in turn, solutions of the so-ealled Roothaan matrix Hartree-Foek equations- Zy Fj y Cy j = 8i Zy Sj y Cy j. One should also note that, just as F ( )i = 8i (l)j possesses a eomplete set of eigenfunetions, the matrix Fj y, whose dimension M is equal to the number of atomie basis orbitals used in the LCAO-MO expansion, has M eigenvalues 8i and M eigenveetors whose elements are the Cy j. Thus, there are oeeupied and virtual moleeular orbitals (mos) eaeh of whieh is deseribed in the LCAO-MO form with Cy j eoeffieients obtained via solution of... 

It is a well-known fact that the Hartree-Fock model does not describe bond dissociation correctly. For example, the H2 molecule will dissociate to an H+ and an atom rather than two H atoms as the bond length is increased. Other methods will dissociate to the correct products however, the difference in energy between the molecule and its dissociated parts will not be correct. There are several different reasons for these problems size-consistency, size-extensivity, wave function construction, and basis set superposition error. 

Hartree-Fock MO approach, the minimization of energy should provide the most accurate description of the electronic field. The mathematical problem is to define each of the terms, with being the most challenging. The formulation carmot be done exactly, but various approaches have been developed and calibrated by comparison with experimental data. The methods used most frequently by chemists were developed by A. D. Becke. " This approach is often called the B3LYP method. The computations can be done with... 

Reproducing the exact solution for the relevant n-electron problem a method ought to yield the same results as the exact solution to the Schrodinger equation to the greatest extent possible. What this means specifically depends on the theory underlying the method. Thus, Hartree-Fock theory should be (and is) able to reproduce the exact solution to the one electron problem, meaning it should be able to treat cases like HeH ... 

The predicted energy of this structure is approximately -215.76438 Hartrees, yielding a reaction barrier of about 86.6 kcal/mol. This value is more in line with expectations, although it is on the high side. Rotation of a double bond is a problem which often requires a higher level of theory than Hartree-Fock (for example, CASSCF) for accurate modeling. 

As we have seen throughout this book, the Hartree-Fock method provides a reasonable model for a wide range of problems and molecular systems. However, Hartree-Fock theory also has limitations. They arise principally from the fact that Hartree-Fock theory does not include a full treatment of the effects of electron correlation the energy contributions arising from electrons interacting with one another. For systems and situations where such effects are important, Hartree-Fock results may not be satisfactory. The theory and methodology underlying electron correlation is discussed in Appendix A. 

Prior to the widespread usage of methods based on Density Functional Theory, the MP2 method was one of the least expensive ways to improve on Hartree-Fock and it was thus often the first correlation method to be applied to new problems. It can successfully model a wide variety of systems, and MP2 geometries are usually quite accurate. Thus, MP2 remains a very useful tool in a computational chemist s toolbox. We ll see several examples of its utility in the exercises. 

AMI does not even get the correct sign for the energy difference. PM3 orders the isomers properly, but substantially underestimates the energy difference, and the Hartree-Fock value is only somewhat better. Electron correlation is needed to obtain even semi-quantitative results for this problem. ... 

Both the MP2 Onsager calculation and the IPCM calculaton are in good agreement with experiment. The SCI-PCM and Hartree-Fock Onsager SCRF calculations perform significantly less well for this problem. ... 

The problem has now become how to solve for the set of molecular orbital expansion coefficients, c. . Hartree-Fock theory takes advantage of the variational principle, which says that for the ground state of any antisymmetric normalized function of the electronic coordinates, which we will denote H, then the expectation value for the energy corresponding to E will always be greater than the energy for the exact wave function ... 

Any method which goes beyond SCF in attempting to treat this phenomenon properly is known as an electron correlation method (despite the fact that Hartree-Fock theory does include some correlation effects) or a post-SCT method. We will look briefly at two different approaches to the electron correlation problem in this section. 

A more general way to treat systems having an odd number of electrons, and certain electronically excited states of other systems, is to let the individual HF orbitals become singly occupied, as in Figure 6.3. In standard HF theory, we constrain the wavefunction so that every HF orbital is doubly occupied. The idea of unrestricted Hartree-Fock (UHF) theory is to allow the a and yS electrons to have different spatial wavefunctions. In the LCAO variant of UHF theory, we seek LCAO coefficients for the a spin and yS spin orbitals separately. These are determined from coupled matrix eigenvalue problems that are very similar to the closed-shell case. 

A Perturbation Theory is developed for treating a system of n electrons in which the Hartree-Fock solution appears as the zero-order approximation. It is shown by this development that the first order correction for the energy and the charge density of the system is zero. The expression for the second order correction for the energy greatly simplifies because of the special property of the zero order solution. It is pointed out that the development of the higher order approximation involves only calculations based on a definite one-body problem. 

I have dealt at length with the Hartree and the Hartree-Foek models. The father of this field, Sir Wilham Hartree, was eoneemed with the atomie problem where it is routinely possible to integrate numerieally the HF integro-differential equations in order to produee (numerieal) wavefunetions that eorrespond to the Hartree-Foek limit. For moleeular applieations the LCAO variant of HF theory assumes a dominant role beeause of the redueed symmetry of the problem. 

Notice that 1 haven t made any mention of the LCAO procedure Hartree produced numerical tables of radial functions. The atomic problem is quite different from the molecular one because of the high symmetry of atoms. The theory of atomic structure is simplified (or complicated, according to your viewpoint) by angular momentum considerations. The Hartree-Fock limit can be easily reached by numerical integration of the HF equations, and it is not necessary to invoke the LCAO method. 

The first two kinds of terms are called derivative integrals, they are the derivatives of integrals that are well known in molecular structure theory, and they are easy to evaluate. Terms of the third kind pose a problem, and we have to solve a set of equations called the coupled Hartree-Fock equations in order to find them. The coupled Hartree-Fock method is far from new one of the earliest papers is that of Gerratt and Mills. 

General expressions for the force constants and dipole derivatives of molecules are derived, and the problems arising from their practical application are reviewed. Great emphasis is placed on the use of the Hartree-Fock function as an approximate wavefunction, and a number... 

A problem with studies on inert gas is that the interactions are so weak. Alkali halides are important commercial compounds because of their role in extractive metallurgy. A deal of effort has gone into corresponding calculations on alkali halides such as LiCl, with a view to understanding the structure and properties of ionic melts. Experience suggests that calculations at the Hartree-Fock level of theory are adequate, provided that a reasonable basis set is chosen. Figure 17.7 shows the variation of the anisotropy and incremental mean pair polarizability as a function of distance. 

The ab initio methods used by most investigators include Hartree-Fock (FFF) and Density Functional Theory (DFT) [6, 7]. An ab initio method typically uses one of many basis sets for the solution of a particular problem. These basis sets are discussed in considerable detail in references [1] and [8]. DFT is based on the proof that the ground state electronic energy is determined completely by the electron density [9]. Thus, there is a direct relationship between electron density and the energy of a system. DFT calculations are extremely popular, as they provide reliable molecular structures and are considerably faster than FFF methods where correlation corrections (MP2) are included. Although intermolecular interactions in ion-pairs are dominated by dispersion interactions, DFT (B3LYP) theory lacks this term [10-14]. FFowever, DFT theory is quite successful in representing molecular structure, which is usually a primary concern. 

The problem of finding the best approximation of this type and the best one-electron set y2t. . ., y>N is handled in the Hartree-Fock scheme. Of course, a total wave function of the same type as Eq. 11.38 can never be an exact solution to the Schrodinger equation, and the error depends on the fact that the two-electron operator (Eq. 11.39) cannot be exactly replaced by a sum of one-particle operators. Physically we have neglected the effect of the "Coulomb hole" around each electron, but the results in Section II.C(2) show that the main error is connected with the neglect of the Coulomb correlation between electrons with opposite spins. 

The best wave function of the approximate form (Eq. 11.38) may then be determined by the variational principle (Eq. II.7), either by varying the quantity p as an entity, subject to the auxiliary conditions (Eq. 11.42), or by varying the basic set fv ip2,. . ., ipN subject to the orthonormality requirement. In both ways we are lead to Hartree-Fock functions pk satisfying the eigenvalue problem... 

---