Hartree-Fock first derivatives

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. 

Gaussian also predicts dipole moments and higher multipole moments (through hexadecapole). The dipole moment is the first derivative of the energy with respect to an applied electric field. It is a measure of the asymmetry in the molecular charge distribution, and is given as a vector in three dimensions. For Hartree-Fock calculations, this is equivalent to the expectation value of X, Y, and Z, which are the quantities reported in the output. 

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. 

The fii st term is zero because I and its derivatives are orthogonal. The fourth term involves second moments and we use the coupled Hartree-Fock procedure to find the terms requiring the first derivative of the wavefunction. 

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

As a consequence, field methods, which consist of computing the energy or dipole moment of the system for external electric field of different amplitudes and then evaluating their first, second derivatives with respect to the field amplitude numerically, cannot be applied. Similarly, procedures such as the coupled-perturbed Hartree-Fock (CPHF) or time-dependent Hartree-Fock (TDHF) approaches which determine the first-order response of the density matrix with respect to the perturbation cannot be applied due to the breakdown of periodicity. 

Further we can proceed similarly as in the case of adiabatic approximation. We shall not present here the details, these are presented in [21,22]. We just mention the most important features of our transformation (46-50). Firstly, when passing from crude adiabatic to adiabatic approximation the force constant changed from second derivative of electron-nuclei interaction ufcF to second derivative of Hartree-Fock energy Therefore when performing transformation (46-50) we expect change offeree constant and therefore change of the vibrational part of Hamiltonian... 

Hartree-Fock models are well defined and yield unique properties. They are both size consistent and variational. Not only may energies and wavefunctions be evaluated from purely analytical (as opposed to numerical) methods, but so too may first and second energy derivatives. This makes such important tasks as geometry optimization (which requires first derivatives) and determination of vibrational frequencies (which requires second derivatives) routine. Hartree-Fock models and are presently applicable to molecules comprising upwards of 50 to 100 atoms. 

Despite the fact that numerical integration is involved, pseudoanalytical procedures have been developed for calculation of first and second energy derivatives. This means that density functional models, like Hartree-Fock models are routinely applicable to determination of equilibrium and transition-state geometries and of vibrational frequencies. 

For quantum chemistry, first-row transition metal complexes are perhaps the most difficult systems to treat. First, complex open-shell states and spin couplings are much more difficult to deal with than closed-shell main group compounds. Second, the Hartree—Fock method, which underlies all accurate treatments in wavefunction-based theories, is a very poor starting point and is plagued by multiple instabilities that all represent different chemical resonance structures. On the other hand, density functional theory (DFT) often provides reasonably good structures and energies at an affordable computational cost. Properties, in particular magnetic properties, derived from DFT are often of somewhat more limited accuracy but are still useful for the interpretation of experimental data. 

Derive the detailed expression for the orbital Hessian for the special case of a closed shell single determinant wave function. Compare with equation (4 53) to check the result. The equation can be used to construct a second order optimization scheme in Hartree-Fock theory. What are the advantages and disadvantages of such a scheme compared to the conventional first order methods ... 

---