Hartree-Fock equation derivation

These are the glorious Hartree-Fock equations derived in general in the spin orbital basis. But wait - there s a problem. These are coupled integro-differential equations, and while they are not strictly unsolvable, they re a pain. It would be nice to at least uncouple them, so let s do that. 

We shall initially consider a closed-shell system with N electroris in N/2 orbitals. The derivation of the Hartree-Fock equations for such a system was first proposed by Roothaan [Roothaan 1951] and (independently) by Hall [Hall 1951]. The resulting equations are known as the Roothaan equations or the Roothaan-Hall equations. Unlike the integro-differential form of the Hartree-Fock equations. Equation (2.124), Roothaan and Hall recast the equations in matrix form, which can be solved using standard techniques and can be applied to systems of any geometry. We shall identify the major steps in the Roothaan approach. 

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. 

When deriving the Hartree-Fock equations it was only required that the variation of the energy with respect to an orbital variation should be zero. This is equivalent to the first derivatives of the energy with respect to the MO expansion coefficients being equal to zero. The Hartree-Fock equations can be solved by an iterative SCF method, and... 

Although a calculation of the wave function response can be avoided for the first derivative, it is necessary for second (and higher) derivatives. Eq. (10.29) gives directly an equation for determining the (first-order) response, which is structurally the same as eq. (10.36). For an HF wave function, an equation of the change in the MO coefficients may also be formulated from the Hartree-Fock equation, eq. (3.50). 

In the Hartree-Fock approach, the many-body wave function in form of a Slater determinant plays the key role in the theory. For instance, the Hartree-Fock equations are derived by minimization of the total energy expressed in terms of this determinantal wave function. In density functional theory (3,4), the fundamental role is taken over by an observable quantity, the electron density. An important theorem of density functional theory states that the correct ground state density, n(r), determines rigorously all electronic properties of the system, in particular its total energy. The totd energy of a system can be expressed as a functional of the density n (r) and this functional, E[n (r)], is minimized by the ground state density. 

Ehf from equation (1-20) is obviously a functional of the spin orbitals, EHF = E[ XJ]. Thus, the variational freedom in this expression is in the choice of the orbitals. In addition, the constraint that the % remain orthonormal must be satisfied throughout the minimization, which introduces the Lagrangian multipliers e in the resulting equations. These equations (1-24) represent the Hartree-Fock equations, which determine the best spin orbitals, i. e., those (xj for which EHF attains its lowest value (for a detailed derivation see Szabo and Ostlund, 1982)... 

The Hartree-Fock approach derives from the application of a series of well defined approaches to the time independent Schrodinger equation (equation 3), which derives from the postulates of quantum mechanics [27]. The result of these approaches is the iterative resolution of equation 2, presented in the previous subsection, which in this case is solved in an exact way, without the approximations of semiempirical methods. Although this involves a significant increase in computational cost, it has the advantage of not requiring any additional parametrization, and because of this the FIF method can be directly applied to transition metal systems. The lack of electron correlation associated to this method, and its importance in transition metal systems, limits however the validity of the numerical results. 

Before deriving the Hartree-Fock equations for the excited state orbitals, we consider the orthogonality constraints imposed on these orbitals. 

VARIATIONAL DERIVATION OF HARTREE-FOCK EQUATIONS FOR EXCITED STATES... 

Derive the closed shell Hartree-Fock equations from the condition that the gradient (4 9) is zero for stationaiy orbitals. The wave function I0> is assumed to be a closed shell Slater determinant. Start by deriving equation (4 41), and from there (4 46). 

The relativistic or non-relativistic random-phase approximation (RRPA or RPA)t is a generalized self-consistent field procedure which may be derived making the Dirac/Hartree-Fock equations time-dependent. Therefore, the approach is often called time-dependent Dirac/Hartree-Fock. The name random phase comes from the original application of this method to very large systems where it was argued that terms due to interactions between many alternative pairs of excited particles, so-called two-particle-two-hole interactions ((2p-2h) see below) tend to... 

As Dewar points out in ref. [30a], this derivation is not really satisfactory. A rigorous approach is a simplified version of the derivation of the Hartree-Fock equations (Chapter 5, Section 5.2.3). It starts with the total molecular wavefunction expressed as a determinant, writes the energy in terms of this wavefunction and the Hamiltonian and finds the condition for minimum energy subject to the molecular orbitals being orthonormal (cf. orthogonal matrices, Section 4.3.3). The procedure is explained in some detail in Chapter 5, Section 5.2.3)... 

The next step in deriving the Hartree-Fock equations is to express the energy of the molecule or atom in terms of the total wavefunction T the energy will then be minimized with respect to each of the component molecular (or atomic an atom is a... 

The KS equations are obtained by differentiating the energy with respect to the KS molecular orbitals, analogously to the derivation of the Hartree-Fock equations, where differentiation is with respect to wavefunction molecular orbitals (Section 5.2.3.4). We use the fact that the electron density distribution of the reference system, which is by decree exactly the same as that of the ground state of our real system (see the definition at the beginning of the discussion of the Kohn-Sham energy), is given by (reference [9])... 

Exact second derivatives methods require the solution of the coupled perturbed Hartree-Fock equations, CPHF [11,34,35]. At the Hartree-Fock level this requires several steps in addition to the usual SCF procedure and the evaluation of the first derivatives. 

It is important that the orbitals r ) that satisfy (5.17) are orthogonal for different eigenvalues We can use them to construct an orthonormal set with which to express the many-electron problem. We have assumed orthogonality in deriving (5.17). We now show that this is consistent by considering (5.17) as an equation for the matrix elements formed for the bra orbital ( from a Schrodinger equation for rj). The set of such equations for N different rj) is called the Hartree—Fock equations. 

The Density Functional theorem states that the total ground state energy is a unique functional of the electron density, p [40]. This simple but enormously powerful result means that it is possible, in principle, to provide an exact description of all electron correlation effects within a one-electron (i.e. orbital-based) scheme. Khon and Sham (KS) [41] have derived a set of equations which embody this result. They have an identical form to the one-electron Hartree Fock equations. The difference is that the exchange-correlation term, Vxc, is not the same. 

---