Moller-Plesset perturbation theory coupled perturbed Hartree-Fock

Total energies E of PH have been calculated with the Moller-Plesset perturbation theory up to the fourth [4 to 7] or second order [8], with a Cl method involving single and double excitations (CISC) [8], with the perturbed Hartree-Fock method (coupled HF scheme) [9], with an MO-SCF method [2], and also with a united-atom (Ar) approximation [10]. 

One of the most dramatic changes in the standard theoretical model used most widely in quantum chemistry occurred in the early 1990s. Until then, ab initio quantum chemical applications [1] typically used a Hartree-Fock (HF) starting point, followed in many cases by second-order Moller-Plesset perturbation theory. For small molecules requiring more accuracy, additional calculations were performed with coupled-cluster theory, quadratic configuration interaction, or related methods. While these techniques are still used widely, a substantial majority of the papers being published today are based on applications of density functional theory (DFT) [2]. Almost universally, the researchers use a functional due to Becke, whose papers in 1992 and 1993 contributed to this remarkable transformation that changed the entire landscape of quantum chemistry. 

Figrire 8 An overview of quantum chemical methods for excited states. Bold-italic entries indicate methods that are currently applicable to large molecules. Important abbreviations used Cl (configuration Interaction), TD (time-dependent), CC (coupled-cluster), HF (Hartree-Fock), CAS (complete active space), RAS (restricted active space), MP (Moller-Plesset perturbation theory), S (singles excitation), SD (singles and doubles excitation), MR (multireference). Geometry optimizations of excited states for larger molecules are now possible with CIS, CASSCF, CC2, and TDDFT methods. 

Although a wide variety of theoretical methods is available to study weak noncovalent interactions such as hydrogen bonding or dispersion forces between molecules (and/or atoms), this chapter focuses on size consistent electronic structure techniques likely to be employed by researchers new to the field of computational chemistry. Not stuprisingly, the list of popular electronic structure techniques includes the self-consistent field (SCF) Hartree-Fock method as well as popular implementations of density functional theory (DFT). However, correlated wave function theory (WFT) methods are often required to obtain accmate structures and energetics for weakly bound clusters, and the most useful of these WFT techniques tend to be based on many-body perturbation theory (MBPT) (specifically, Moller-Plesset perturbation theory), quadratic configuration interaction (QCI) theory, and coupled-cluster (CC) theory. 

In a very recent study, Banerjee et al used different ab initio methods to study the properties of small Kjv clusters for even A with 2 < A < 20. They used both density-functional methods and post-Hartree-Fock approaches, where correlation is added either via Moller-Plesset perturbation theory or via the coupled-cluster approach (see, e.g., ref. 1). In order to determine the structures of the lowest total energy, they used as initial guesses structures from earlier studies on Nuat clusters that subsequently were allowed to relax to their closest total-energy minimum. Unfortunately, only the density-functional method was used in optimizing the structures, whereby a comparison between the two approaches is not made possible. 

Maroulis and Haskopoulos calculated the interaction electric dipole moment and polarizability for the C02-Rg systems, Rg = He, Ne, Ar, Kr and Xe. The potential minimum is very well defined for all these systems. In Fig. 19 is shown the potential energy surface for the C02-He interaction calculated at the MP2 level of theory. The most stable configuration corresponds to a T-shaped structure. The two local minima for the linear configuration of C02-He are also clearly visible. All interaction induced properties were extracted from finite-filed Moller-Plesset perturbation theory and coupled-cluster calculations with purpose-oriented basis sets. CCSD(T) values were calculated for the dipole moment pim of C02-He and C02-Ne the corresponding results are 0.0063 and 0.0107 eao, respectively. All post-Hartree-Fock methods yield stable values for this important property. For C02-He, = 0.0070 (SCF), 0.0063 (MP2), 0.0063 (MP4),... 

The calculation of molecular properties can be carried out at three distinct levels (i) ab initio, (ii) semi-empirical, (iii) empirical. Ab initio methods have increased enormously in accuracy and efficiency in the last two decades and are the focus of our discussion here. Ab initio methods have developed in two directions first, the level of approximation has become increasingly sophisticated and, hence, accurate. The earliest ab initio calculations used the Hartree-Fock/self-consistent field (HF/SCF) methodology, which is the simplest to implement. Subsequently, such methods as Moller-Plesset perturbation theory, multi-configuration self-consistent field theory (MCSCF) and coupled-cluster theory have been developed and implemented. Relatively recently, density functional theory (DFT) has become very popular, since it yields an accuracy much greater than that of HF/SCF while requiring relatively little additional computational effort. 

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. 

The set of atomic orbitals Xk is called a basis set, and the quality of the basis set will usually dictate the accuracy of the calculations. For example, the interaction energy between an active site and an adsorbate molecule might be seriously overestimated because of excessive basis set superposition error (BSSE) if the number of atomic orbitals taken in Eq. [4] is too small. Note that Hartree-Fock theory does not describe correlated electron motion. Models that go beyond the FiF approximation and take electron correlation into account are termed post-Flartree-Fock models. Extensive reviews of post-HF models based on configurational interaction (Cl) theory, Moller-Plesset (MP) perturbation theory, and coupled-cluster theory can be found in other chapters of this series. ... 

The question for a more systematic inclusion of electronic correlation brings us back to the realm of molecular quantum chemistry [51,182]. Recall that (see Section 2.11.3) the exact solution (configuration interaction. Cl) is found on the basis of the self-consistent Hartree-Fock wave function, namely by the excitation of the electrons into the virtual, unoccupied molecular orbitals. Unfortunately, the ultimate goal oi full Cl is obtainable for very small systems only, and restricted Cl is size-inconsistent the amount of electron correlation depends on the size of the system (Section 2.11.3). Thus, size-consistent but perturbative approaches (Moller-Plesset theory) are often used, and the simplest practical procedure (of second order, thus dubbed MP2 [129]) already scales with the fifth order of the system s size N, in contrast to Hartree-Fock theory ( N ). The accuracy of these methods may be systematically improved by going up to higher orders but this makes the calculations even more expensive and slow (MP3 N, MP4 N ). Fortunately, restricted Cl can be mathematically rephrased in the form of the so-called coupled clus-... 

---