Post Hartree-Fock Schemes

The average treatment of electron correlation is inadequate in single determinant HF theory. Several schemes exist to remedy this, ranging from the explicit inclusion of Configuration Interaction to the use of perturbation theory. The formalism for improving on the HF approximation is well defined. However, the practical implementation of Cl is problematic. 

The two main concerns are (1) Cl scales badly (N5 to N8) and exacts too large a penalty in term of computational resources and (2) Cl is not always well-behaved in as much as an apparently systematic improvement of the Cl expansion does not necessarily lead to a systematic increase in the accuracy of the computed results. 

However, while post-HF methods can, in principle, provide a systematic way of improving the calculation until any desired level of accuracy is achieved, the 

---

DFT and post-Hartree-Fock ab initio studies on the different tautomers and rotamers of 2-hydroxy- (94) and 2,3-dihydroxypyrazine (95) indicate that the former species is stabilized by about -3 kcal/mol with respect to its keto tautomer in the gas phase [99JST229]. In solution the opposite appears to be true. For 95, the hydroxyoxo 95b and diketo tautomers 95c are most stable in the gas phase (Scheme 61). 

Trioxane 210 has been used as a model system by Gu and coworkers to study the antimalarial drug artemisinin 211 (Scheme 137) [97CPL234, 99JST103]. It is the boat/twist form rather than the chair conformer of 210 that describes the subunit in 211. Moreover, geometric parameters and vibrational frequencies can only reliably be computed at the DFT level and by post-Hartree-Fock methods. B3-LYP/6-31G calculations on the conformers of 3,3,6,6-tetramethyl-1,2,4,5-tetroxane show that the chair conformer is stabilized with respect to the twisted conformer by about -2.8 kcal/mol [00JST85]. No corresponding boat conformer was found. 

Structure and reactivity of (thionitroso)ethylene 9 have been examined by ab initio methods at the post-Hartree-Fock (HF) level of theory. This compound is expected to undergo electrocyclic ring closure to the more stable 4/7-1,2-thiazete 3. Geometrical bond length and dihedral angle parameters at the stationary points of the electrocyclic reactions of thionitrosoethylene were calculated at the RMP2/6-31G level of theory (Scheme 1) <1996LA1615>. 

This description of quantum mechanical methods for computing (hyper)polarizabilities demonstrates why, nowada, the determination of hyperpolarizabilities of systems containing hundreds of atoms can, at best, be achieved by adopting, for obvious computational reasons, semi-empirical schemes. In this study, the evaluation of the static and dynamic polarizabilities and first hyperpolaiizabilities was carried out at die Time-Dependent Hartree-Fock (TDOT) [39] level with the AMI [50] Hamiltonian. The dipole moments were also evaluated using the AMI scheme. The reliability of the semi-empirical AMI calculations was addressed in two ways. For small and medium-size push-pull polyenes, the TDHF/AMl approach was compared to Hartree-Fock and post Hartree-Fock [51] calculations of die static and dynamic longitudinal first hyperpolarizability. Except near resonance, the TDHF/AMl scheme was shown to perform appreciably better than the ab initio TDHF scheme. Then, the static electronic first hyperpolaiizabilities of the MNA molecule and dimer have been calculated [15] with various ab initio schemes and compared to the AMI results. In particular, the inclusion of electron correlation at the MP2 level leads to an increase of Paaa by about 50% with respect to the CPHF approach, similar to the effect calculated by Sim et al. [52] for the longitudinal p tensor component of p-nitroaniline. The use of AMI Hamiltonian predicts a p aa value that is smaller than the correlated MP2/6-31G result but larger than any of the CPHF ones, which results fi-om the implicit treatment of correlation effects, characteristic of die semi-empirical methods. This comparison confirms that a part of die electron... 

This is not the place for a full overview of the wave function based post Hartree-Fock methods currently applied for the calculation of intermolecular interactions and in particular molecule/surface interactions. Table 3 contains a brief characterization of the most widely applied schemes. The two most popular methods are MP2 (second order Moller-Plesset perturbation theory), because it covers large part of electronic correlation at comparably low ex-... 

The majority of routine calculations in quantum chemistry are done with variational methods (mainly the Hartree-Fock scheme). If we consider post-Hartree-Fock calculations then non-variational [CCSD, CCSD(T)] as well as perturbational (among them MBPT) approaches take the lead. The perturbational methods are based on the simple idea that the system, in a slightly modified condition, is similar to that before the perturbation is applied (cf. p. 203). 

There is no fundamental change in the concept of correlation between relativistic and nonrelativistic quantum chemistry in both cases, correlation describes the difference between a mean-field description, which forms the reference state for the correlation method, and the exact description. We can also define dynamical and non-dynamical correlation in both cases. There is in fact no formal difference between a nonrelativistic spin-orbital-based formalism and a relativistic spinor-based formalism. Thus we should be able to transfer most of the schemes for post-Hartree-Fock calculations to a relativistic post-Dirac-Hartree-Fock model. Several such schemes have been implemented and applied in a range of calculations. The main technical differences to consider are those arising from having to deal with integrals that are complex, and the need to replace algorithms that exploit the nonrelativistic spin symmetry by schemes that use time-reversal and double-group symmetry. 

---