Correlation, electron Mpller-Plesset

So far, this method has only been applied to one-dimensional periodic systems, but, except for the involved coding, there is no perceived difficulty to extend it to three-dimensional crystals. As shown by studies on oligomers, these TDHF results suffer from the lack of electron correlation which can modify estimates by as much as one order of magnitude. Generalizing this approach by using traditional electron correlation methods (Mpller Plesset and coupled-cluster) or ad hoc DFT treatments is therefore a step to be pursued. 

One of us [1] reviewed the situation of electron correlation a quarter of a century ago in a paper with the title electron correlation in the seventies [2]. At that time most quantum chemists did not care about electron correlation, and standard methods for the large scale treatment of electron correlation, like Mpller-Plesset (MP) perturbation theory or coupled-cluster (CC) theory were not yet available. However precursors of these methods such as lEPA (independent electron pair approximation) and CEPA (coupled-electron-pair approximation) had already been developped and were being used, mainly in research groups in Germany [3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13]. 

HyperChem supports MP2 (second order Mpller-Plesset) correlation energy calculationsusing afe mi/io methods with anyavailable basis set. In order to save main memory and disk space, the HyperChem MP2 electron correlation calculation normally uses a so called frozen-core approximation, i.e. the inner shell (core) orbitals are omitted. A setting in CHEM.INI allows excitations from the core orbitals to be included if necessary (melted core). Only the single point calculation is available for this option. 

Specifies the calculation of electron correlation energy using the Mpller-Plesset second order perturbation theory (MP2). This option can only be applied to Single Point calculations. 

Things have moved on since the early papers given above. The development of Mpller-Plesset perturbation theory (Chapter 11) marked a turning point in treatments of electron correlation, and made such calculations feasible for molecules of moderate size. The Mpller-Plesset method is usually implemented up to MP4 but the convergence of the MPn series is sometimes unsatisfactory. The effect... 

How does a rigorously calculated electrostatic potential depend upon the computational level at which was obtained p(r) Most ab initio calculations of V(r) for reasonably sized molecules are based on self-consistent field (SCF) or near Hartree-Fock wavefunctions and therefore do not reflect electron correlation in the computation of p(r). It is true that the availability of supercomputers and high-powered work stations has made post-Hartree-Fock calculations of V(r) (which include electron correlation) a realistic possibility even for molecules with 5 to 10 first-row atoms however, there is reason to believe that such computational levels are usually not necessary and not warranted. The Mpller-Plesset theorem states that properties computed from Hartree-Fock wave functions using one-electron operators, as is T(r), are correct through first order (Mpller and Plesset 1934) any errors are no more than second-order effects. 

The difference between the Hartree-Fock energy and the exact solution of the Schrodinger equation (Figure 60), the so-called correlation energy, can be calculated approximately within the Hartree-Fock theory by the configuration interaction method (Cl) or by a perturbation theoretical approach (Mpller-Plesset perturbation calculation wth order, MPn). Within a Cl calculation the wave function is composed of a linear combination of different Slater determinants. Excited-state Slater determinants are then generated by exciting electrons from the filled SCF orbitals to the virtual ones ... 

As usual, the Hartree-Fock model can be corrected with perturbation theory (e.g., the Mpller-Plesset [MP] method29) and/or variational techniques (e.g., the configuration-interaction [Cl] method30) to account for electron-correlation effects. The electron density p(r) = N f P 2 d3 2... d3r can generally be expressed as... 

There is also a hierarchy of electron correlation procedures. The Hartree-Fock (HF) approximation neglects correlation of electrons with antiparallel spins. Increasing levels of accuracy of electron correlation treatment are achieved by Mpller-Plesset perturbation theory truncated at the second (MP2), third (MP3), or fourth (MP4) order. Further inclusion of electron correlation is achieved by methods such as quadratic configuration interaction with single, double, and (perturbatively calculated) triple excitations [QCISD(T)], and by the analogous coupled cluster theory [CCSD(T)] [8],... 

Prior to stretching C-S bond, we optimized the geometry of the anionic Me-S-Me molecule and parent neutral molecule at the unrestricted second-order Mpller-Plesset (UMP2) perturbation level of theory (in order to take into account the effect of electron correlation) with aug-cc-pVDZ basis sets [9]. We also... 

Subsequently, we performed more sophisticated calculations with a larger basis set and with inclusion of electron correlation at the level of Mpller-Plesset theory. From Figure 7 it is seen that the agreement between experiment and theory is even better than previously. Similar agreement was found for [1,2,3,4- H4]-1. ... 

In ab initio methods the HER approximation is used for build-up of initial estimate for and which have to be further improved by methods of configurational interaction in the complete active space (CAS) [39], or by Mpller-Plesset perturbation theory (MPn) of order n, or by the coupled clusters [40,41] methods. In fact, any reasonable result within the ab initio QC requires at least minimal involvement of electron correlation. All the technical tricks invented to go beyond the HFR calculation scheme in terms of different forms of the trial wave function or various perturbative procedures represent in fact attempts to estimate somehow the second term of Eq. (5) - the cumulant % of the two-particle density matrix. 

Under some simplifications associated with the symmetry of fullerenes, it has been possible to perform calculations of type Hartree-Fock in which the interelec-tronic correlation has been included up to second order Mpller-Plesset (Moller et al. 1934 Purcell 1979 Cioslowski 1995), and calculations based on the density functional (Pople et al. 1976). However, given the difficulties faced by ab initio computations when all the electrons of these large molecules are taken into account, other semiempirical methods of the Hiickel type or tight-binding (Haddon 1992) models have been developed to determine the electronic structure of C60 (Cioslowski 1995 Lin and Nori 1996) and associated properties like polarizabilities (Bonin and Kresin 1997 Rubio et al. 1993) hyperpolarizabilities (Fanti et al. 1995) plasmon excitations (Bertsch et al. 1991) etc. These semiempirical models reproduce the order of monoelectronic levels close to the Fermi level. Other more sophisticated semiempirical models, like the PPP (Pariser-Parr-Pople) (Pariser and Parr 1953 Pople 1953) obtain better quantitative results when compared with photoemission experiments (Savage 1975). 

Instead, practical methods involve a subset of possible Slater determinants, especially those in which two electrons are moved from the orbitals they occupy in the HF wavefunction into empty orbitals. These doubly excited determinants provide a description of the physical effect missing in HF theory, correlation between the motions of different electrons. Single and triple excitations are also included in some correlated ab initio methods. Different methods use different techniques to decide which determinants to include, and all these methods are computationally more expensive than HF theory, in some cases considerably more. Single-reference correlated methods start from the HF wavefunction and include various excited determinants. Important methods in inorganic chemistry include Mpller-Plesset perturbation theory (MP2), coupled cluster theory with single and double excitations (CCSD), and a modified form of CCSD that also accounts approximately for triple excitations, CCSD(T). 

The Mpller-Plesset (MP) treatment of electron correlation [84] is based on perturbation theory, a very general approach used in physics to treat complex systems [85] this particular approach was described by M0ller and Plesset in 1934 [86] and developed into a practical molecular computational method by Binkley and Pople [87] in 1975. The basic idea behind perturbation theory is that if we know how to treat a simple (often idealized) system then a more complex (and often more realistic) version of this system, if it is not too different, can be treated mathematically as an altered (perturbed) version of the simple one. Mpller-Plesset calculations are denoted as MP, MPPT (M0ller-Plesset perturbation theory) or MBPT (many-body perturbation theory) calculations. The derivation of the Mpller-Plesset method [88] is somewhat involved, and only the flavor of the approach will be given here. There is a hierarchy of MP energy levels MPO, MP1 (these first two designations are not actually used), MP2, etc., which successively account more thoroughly for interelectronic repulsion. 

The configuration interaction (Cl) treatment of electron correlation [83,95] is based on the simple idea that one can improve on the HF wavefunction, and hence energy, by adding on to the HF wavefunction terms that represent promotion of electrons from occupied to virtual MOs. The HF term and the additional terms each represent a particular electronic configuration, and the actual wavefunction and electronic structure of the system can be conceptualized as the result of the interaction of these configurations. This electron promotion, which makes it easier for electrons to avoid one another, is as we saw (Section 5.4.2) also the physical idea behind the Mpller-Plesset method the MP and Cl methods differ in their mathematical approaches. 

The coupled cluster (CC) method is actually related to both the perturbation (Section 5.4.2) and the Cl approaches (Section 5.4.3). Like perturbation theory, CC theory is connected to the linked cluster theorem (linked diagram theorem) [101], which proves that MP calculations are size-consistent (see below). Like standard Cl it expresses the correlated wavefunction as a sum of the HF ground state determinant and determinants representing the promotion of electrons from this into virtual MOs. As with the Mpller-Plesset equations, the derivation of the CC equations is complicated. The basic idea is to express the correlated wave-function Tasa sum of determinants by allowing a series of operators 7), 73,... to act on the HF wavefunction ... 

First let us review static and dynamic electron correlation. Dynamic (dynamical) electron correlation is easy to grasp, if not so easy to treat exhaustively. It is simply the adjustment by each electron, at each moment, of its motion in accordance with its interaction with each other electron in the system. Dynamic correlation and its treatment with perturbation (Mpller-Plesset), configuration interaction, and coupled cluster methods was covered in Section 5.4. 

An ab initio version of the Mpller-Plesset perturbation theory within the DPCM solvation approach was introduced years ago by Olivares et al. [26] following the above intuitive considerations based on the fact that the electron correlation which modifies both the HF solute charge distribution and the solvent reaction potential depending on it can be back-modified by the solvent. To decouple these combined effects the authors introduced three alternative schemes ... 

In this section, we briefly discuss some of the electronic structure methods which have been used in the calculations of the PE functions which are discussed in the following sections. There are variety of ab initio electronic structure methods which can be used for the calculation of the PE surface of the electronic ground state. Most widely used are Hartree-Fock (HF) based methods. In this approach, the electronic wavefunction of a closed-shell system is described by a determinant composed of restricted one-electron spin orbitals. The unrestricted HF (UHF) method can handle also open-shell electronic systems. The limitation of HF based methods is that they do not account for electron correlation effects. For the electronic ground state of closed-shell systems, electron correlation effects can be accounted for relatively easily by second-order Mpller-Plesset perturbation theory (MP2). In modern implementations of MP2, linear scaling with the size of the system has been achieved. It is thus possible to treat quite large molecules and clusters at this level of theory. 

---