Correlation energy methods

The tautomerism of 2-hydroxypyridine in the gas phase has also been studied by microwave spectroscopy using both a conventional spectrometer and a jet-cooled millimeter-wave spectrometer. The spectra attributable to both (Z)-hydroxy tautomer and 2-pyridinone were observed in all cases, the relative abundance being 3 1 in favor of the hydroxy form. No ( )-hydroxy isomer has been detected (93JPC46). The increased stability of the <7,v-hydroxy tautomer compared to the tra/iv-isomer is confirmed by CNDO/2 calculations (AE— 0.64 kcal/mol), whereas the semiempirical effective pair correlation energy method favors the fra/iv-isomer by 0.69 kcal/mol (90BCJ2981). 

MO The atomic orbital integrals are generated once and stored externally, then transformed to the molecular orbital basis. The transformed (MO) integrals are also stored externally. The four index transformation, eq 11, requires work and disk storage (see the Integral Transformation section). This is the approach used until recently for all correlated energy methods. 

Although it is now somewhat dated, this book provides one of the best treatments of the Hartree-Fock approximation and the basic ideas involved in evaluating the correlation energy. An especially valuable feature of this book is that much attention is given to how these methods are actually implemented. 

Size-extensivity is of importance when one wishes to compare several similar systems with different numbers of atoms (i.e., methanol, ethanol, etc.). In all cases, the amount of correlation energy will increase as the number of atoms increases. However, methods that are not size-extensive will give less correlation energy for the larger system when considered in proportion to the number of electrons. Size-extensive methods should be used in order to compare the results of calculations on different-size systems. Methods can be approximately size-extensive. The size-extensivity and size-consistency of various methods are summarized in Table 26.1. 

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. 

The main difference between the G1/G2 and CBS methods is the way in which they try to extrapolate the correlation energy, as described below. Both tire G1/G2 and CBS methods come in different flavours, depending on the exact combinations of metliods used for obtaining the above four contributions. 

The Parameterized Configuration Interaction (PCI-X) method simply takes the correlation energy and scales it by a constant factor X (typical value 1.2), i.e. it is assumed that the given combination of method and basis set recovers a constant fraction of the correlation energy. 

G2/CBS models, essentially all die correlation energy of the bond being broken must be recovered. This in mm necessitates large basis sets and sophisticated correlation methods. This is also the reason why ab initio energies are not converted into heats of formation, as is normally done for semi-empirical methods (eq. (3.89)), since the resulting values are poor unless a very high level of theory is employed. 

For closed-shell systems LSDA is equal to LDA, and since this is the most common case, LDA is often used interchangeably with LSDA, although this is not true in the general case (eqs. (6.16) and (6.17)). The method proposed by Slater in 1951 can be considered as an LDA mediod where die correlation energy is neglected and the exchange term is given as... 

The correlation energy of a uniform electron gas has been determined by Monte Carlo methods for a number of different densities. In order to use these results in DFT calculations, it is desirable to have a suitable analytic interpolation formula. This has been constructed by Vosko, Wilk and Nusair (VWN) and is in general considered to be a very accurate fit. It interpolates between die unpolarized ( = 0) and spin polarized (C = 1) limits by the following functional. 

The LSDA approximation in general underestimates the exchange energy by 10%, thereby creating errors which are larger tlian the whole correlation energy. Electron correlation is furthermore overestimated, often by a factor close to 2, and bond strengths are as a consequence overestimated. Despite the simplicity of the fundamental assumptions, LSDA methods are often found to provide results with an accuracy similar to that obtained by wave mechanics HE methods. 

The magnitude of the core correlation can be evaluated by including the oxygen Is-electrons and using the cc-pCVXZ basis sets the results are shown in Table 11.9. The extrapolated CCSD(T) correlation energy is —0.370 a.u. Assuming that the CCSD(T) method provides 99.7% of the full Cl value, as indicated by Table 11.7, the extrapolated correlation energy becomes —0.371 a.u., well within the error limits on the estimated experimental value. The core (and core-valence) electron correlation is thus 0.063 a.u.. 

---