Hartree-Fock orbitals, restricted

If we except the Density Functional Theory and Coupled Clusters treatments (see, for example, reference [1] and references therein), the Configuration Interaction (Cl) and the Many-Body-Perturbation-Theory (MBPT) [2] approaches are the most widely-used methods to deal with the correlation problem in computational chemistry. The MBPT approach based on an HF-SCF (Hartree-Fock Self-Consistent Field) single reference taking RHF (Restricted Hartree-Fock) [3] or UHF (Unrestricted Hartree-Fock ) orbitals [4-6] has been particularly developed, at various order of perturbation n, leading to the widespread MPw or UMPw treatments when a Moller-Plesset (MP) partition of the electronic Hamiltonian is considered [7]. The implementation of such methods in various codes and the large distribution of some of them as black boxes make the MPn theories a common way for the non-specialist to tentatively include, with more or less relevancy, correlation effects in the calculations. 

Both and can involve double replacements with respect to the single determinantal reference function constructed from canonical Hartree-Fock orbitals. The intermediate state can involve single, double, triple and quadruple replacements with respect to the reference function. The renormalization terms in (44) involve only double replacement. Brueckner91 showed that diagonal terms in the principal term which have a non linear dependence of N are completely cancelled by the renormalization terms. This cancellation is incomplete if the level of replacement employed in generating the intermediate states is restricted. 

For the lighter elements we have available the VSIP data defined for orbitals of definite occupancy (11, 12, 13). For the heavier elements we have only the Hartree-Fock orbital energies (20, 22). Thus, for the lighter combinations we used either source of data while if either element is heavy, our choice is restricted. It seems preferable to use data... 

Let us first consider the calculations on the H2 system. In Table 12.2, the MP2 natural orbitals are used in the first iteration in Table 12.3, we use the canonical Hartree-Fock orbitals. Because of the different choices of orbitals, the optimizations proceed rather differently. For the optimization based on the MP2 natural orbitals, the optimization begins in the local region each step corresporxls to a Newton step with no step-length restrictions. Quadratic convergence is therefore observed in all outer iterations - see the reduction in the gradient and step norms in Table 12.2. The ratio parameter r (12.3.21), which probes the quadratic dominance of the energy function, is close to 1 in all iterations. 

In addition, using restricted Hartree-Fock orbitals, jp2 m l is the matrix element over spatial orbitals assumed to the independent of spin. The matrix M describes multipair excitations. In M = 0, Eqs. (2.33) are separable. Then their diagonal elements yield... 

The structure of /,p ( ) may be determined using many-body perturbation theory. If one chooses self-consistence field-restricted Hartree-Fock orbitals, all diagrams contributing to canonical Hartree-Fock are omitted. The second-order expression for (w ) is... 

A variation on the HF procedure is the way that orbitals are constructed to reflect paired or unpaired electrons. If the molecule has a singlet spin, then the same orbital spatial function can be used for both the a and P spin electrons in each pair. This is called the restricted Hartree-Fock method (RHF). 

Another way of constructing wave functions for open-shell molecules is the restricted open shell Hartree-Fock method (ROHF). In this method, the paired electrons share the same spatial orbital thus, there is no spin contamination. The ROHF technique is more difficult to implement than UHF and may require slightly more CPU time to execute. ROHF is primarily used for cases where spin contamination is large using UHF. 

Introductory descriptions of Hartree-Fock calculations [often using Rootaan s self-consistent field (SCF) method] focus on singlet systems for which all electron spins are paired. By assuming that the calculation is restricted to two electrons per occupied orbital, the computation can be done more efficiently. This is often referred to as a spin-restricted Hartree-Fock calculation or RHF. 

You can order the molecular orbitals that are a solution to equation (47) according to their energy. Electrons populate the orbitals, with the lowest energy orbitals first. Anormal, closed-shell, Restricted Hartree Fock (RHF) description has a maximum of two electrons in each molecular orbital, one with electron spin up and one with electron spin down, as shown ... 

The Roothaan equations just described are strictly the equations for a closed-shell Restricted Hartree-Fock (RHF) description only, as illustrated by the orbital energy level diagram shown earlier. To be more specific ... 

A restricted Hartree-Fock description means that spin-up and spin-down electrons occupy the same spatial orbitals /j—there is no allowance for different spatial orbitals for different electron spins. 

You will need to decide whether or not to request Restricted (RHF) or Unrestricted (UHF) Hartree-Fock calculations. This question embodies a certain amount of controversy and there is no simple answer. The answer often depends simply on which you prefer or what set of scientific prejudices you have. Ask yourself whether you prefer orbital energy diagrams with one or two electrons per orbital. 

A UHF wave function may also be a necessary description when the effects of spin polarization are required. As discussed in Differences Between INDO and UNDO, a Restricted Hartree-Fock description will not properly describe a situation such as the methyl radical. The unpaired electron in this molecule occupies a p-orbital with a node in the plane of the molecule. When an RHF description is used (all the s orbitals have paired electrons), then no spin density exists anywhere in the s system. With a UHF description, however, the spin-up electron in the p-orbital interacts differently with spin-up and spin-down electrons in the s system and the s-orbitals become spatially separate for spin-up and spin-down electrons with resultant spin density in the s system. 

Here we give the molecule specification in Cartesian coordinates. The route section specifies a single point energy calculation at the Hartree-Fock level, using the 6-31G(d) basis set. We ve specified a restricted Hartree-Fock calculation (via the R prepended to the HF procedure keyword) because this is a closed shell system. We ve also requested that information about the molecular orbitals be included in the output with Pop=Reg. 

So far, we have considered only the restricted Hartree-Fock method. For open shell systems, an unrestricted method, capable of treating unpaired electrons, is needed. For this case, the alpha and beta electrons are in different orbitals, resulting in two sets of molecular orbital expansion coefficients ... 

Here, occ means occupied and virt means virtual. In the restricted Hartree-Fock model, each orbital can be occupied by at most one a spin and one (i spin electron. That is the meaning of the (redundant) Alpha in the output. In the unrestricted Hartree-Fock model, the a spin electrons have a different spatial part to the spin electrons and the output consists of the HF-LCAO coefficients for both the a spin and the spin electrons. 

So far there have not been any restrictions on the MOs used to build the determinantal trial wave function. The Slater determinant has been written in terms of spinorbitals, eq. (3.20), being products of a spatial orbital times a spin function (a or /3). If there are no restrictions on the form of the spatial orbitals, the trial function is an Unrestricted Hartree-Fock (UHF) wave function. The term Different Orbitals for Different Spins (DODS) is also sometimes used. If the interest is in systems with an even number of electrons and a singlet type of wave function (a closed shell system), the restriction that each spatial orbital should have two electrons, one with a and one with /3 spin, is normally made. Such wave functions are known as Restricted Hartree-Fock (RHF). Open-shell systems may also be described by restricted type wave functions, where the spatial part of the doubly occupied orbitals is forced to be the same this is known as Restricted Open-shell Hartree-Fock (ROHF). For open-shell species a UHF treatment leads to well-defined orbital energies, which may be interpreted as ionization potentials. Section 3.4. For an ROHF wave function it is not possible to chose a unitary transformation which makes the matrix of Lagrange multipliers in eq. (3.40) diagonal, and orbital energies from an ROHF wave function are consequently not uniquely defined, and cannot be equated to ionization potentials by a Koopman type argument. 

In the ordinary Hartree-Fock scheme, the total wave function is approximated by a single Slater determinant and, if the system possesses certain symmetry properties, they may impose rather severe restrictions on the occupied spin orbitals see, e.g., Eq. 11.61. These restrictions may be removed and the total energy correspondingly decreased, if instead we approximate the total wave function by means of the first term in the symmetry adapted set, i.e., by the projection of a single determinant. Since in both cases,... 

---