Bond length alternation in cyclic polyenes. V. Local finite-order perturbation theory approach

The finite-order many-body perturbation theory using the localized Wannier orbital basis is applied to the problem of bond length alternation in the Pariser–Parr–Pople model of cyclic polyenes C_N H_N, N = 4v + 2, which may be regarded as a simplified model of polyacetylene. Both the Møller–Plesset and the Epstein–Nesbet-type partitionings of the model Hamiltonian are employed. The localized orbital basis enables an efficient truncation of the perturbation theory summations over the intermediate states as well as an elimination of energetically unimportant diagrams, thus enabling one to obtain the fourth-order Møller–Plesset-perturbation energies with a relatively small computational effort even for large polyenes. The results obtained with the second-, third-, and fourth-order Møller–Plesset and with the third-order Epstein–Nesbet perturbation theories yield very similar bond length distortions (about 0.05 Å) and stabilization energies per site (about 0.04 eV) as obtained earlier with the RHF , one-parameter AMO , and delocalized orbital perturbation theories. The effects of truncation and diagram elimination in the fourth-order Møller–Plesset perturbation theory and the abnormal behavior of the second-order Epstein–Nesbet perturbation theory results in the localized Wannier basis near the instability threshold of the RHF solutions are discussed.     

Different forms of perturbation theory for the calculation of the correlation energy

Different forms of perturbation theory for the calculation of correlation energy in both closed-and open-shell systems are discussed. For closed-shell systems, Epstein–Nesbet perturbation theory is compared with Møller–Plesset (MP ) perturbation theory based on canonical Hartree–Fock orbitals and with MP theory based on internally consistent SCF orbitals. The traditional MP theory gives superior results despite its use of an inferior zeroth-order Hamiltonian. This behavior is rationalized in terms of the larger denominators present in the traditional MP theory. These conclusions are used to support the restricted open-shell perturbation methods proposed recently by Murray and Davidson, and these new methods are compared with spin-restricted Epstein–Nesbet theory and the unrestricted MP (UMP ) approach. © 1992 John Wiley & Sons, Inc.     

Many-body multireference Møller—Plesset and Epstein—Nesbet perturbation theory: Fast evaluation of second-order energy contributions

Multireference perturbation theory is examined in connection with the two partitions in the Møller—Plesset and Epstein—Nesbet schemes. The implementation of an efficient diagrammatic technique is described and two examples of application (diazene and the Cr₂ molecule), involving large variational spaces, are provided. © 1996 John Wiley & Sons, Inc.     

Quadratic Padé approximants and the intruder state problem of multireference perturbation methods

Simple and quadratic Padé resummation methods are applied to high‐order series from multireference many‐body perturbation theory (MR‐MBPT) calculations using various partitioning schemes (Møller–Plesset, Epstein–Nesbet, and forced degeneracy) to determine their efficacy in resumming slowly convergent or divergent series. The calculations are performed for the ground and low‐lying excited states of (i) CH₂, (ii) BeH₂ at three geometries, and (iii) Be, for which full configuration interaction (CI) calculations are available for comparison. The 49 perturbation series that are analyzed include those with oscillatory and monotonic divergence and convergence, including divergences that arise from either frontdoor or backdoor intruder states. Both the simple and quadratic Padé approximations are found to speed the convergence of slowly convergent or divergent series. However, the quadratic Padé method generally outperforms the simple Padé resummation. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Theory of relaxed density matrices: Application to second-order response properties

The evaluation of second derivatives of the electronic energy for nonvariational wave functions using an energy functional is discussed. It is shown that, in certain cases, the formation of the first-order relaxed density matrix leads to an efficient algorithm for the calculation of second-order response properties. Detailed formulas are given for second-order Møller–Plesset perturbation theory. © 1994 John Wiley & Sons, Inc.     

Multi‐state multi‐reference Møller–Plesset second‐order perturbation theory for molecular calculations

This work presents multi‐state multi‐reference Møller–Plesset second‐order perturbation theory as a variant of multi‐reference perturbation theory to treat electron correlation in molecules. An effective Hamiltonian is constructed from the first‐order wave operator to treat several strongly interacting electronic states simultaneously. The wave operator is obtained by solving the generalized Bloch equation within the first‐order interaction space using a multi‐partitioning of the Hamiltonian based on multi‐reference Møller–Plesset second‐order perturbation theory. The corresponding zeroth‐order Hamiltonians are nondiagonal. To reduce the computational effort that arises from the nondiagonal generalized Fock operator, a selection procedure is used that divides the configurations of the first‐order interaction space into two sets based on the strength of the interaction with the reference space. In the weaker interacting set, only the projected diagonal part of the zeroth‐order Hamiltonian is taken into account. The justification of the approach is demonstrated in two examples: the mixing of valence Rydberg states in ethylene, and the avoided crossing of neutral and ionic potential curves in LiF. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

BSSE‐free description of the formamide dimers

The different configurations (linear, zig‐zag, and cyclic) of formamide dimers have been studied at the level of both Hartree–Fock (HF) and second order Møller–Plesset perturbation theory (MP2). The widely used a posteriori Boys–Bernardi “counterpoise” (CP) correction scheme has been compared with our a priori methods utilizing the “chemical Hamiltonian approach” (CHA). The appropriate interaction energies have been calculated in six different basis sets (6‐31G, 6‐31G**, DZV, DZP, TZV, and cc‐pVDZ). © 2001 John Wiley & Sons, Inc. Int J Quant Chem, 2001     

Fully optimized implementation of the cluster-in-molecule local correlation approach for electron correlation calculations of large systems

A fully optimized implementation of the cluster-in-molecule (CIM) local correlation method for faster and more accurate electron correlation calculations of large systems is reported. The speedup comes from the new procedure of constructing virtual localized molecular orbitals of clusters. In the new procedure, Boughton–Pulay projection method is employed instead of a much more expensive Boys localization procedure. In addition, basis set superposition error correction for binding energy calculations and parallelized electron correlation calculations of clusters are now implemented. Benchmark calculations and illustrative applications at the Møller–Plesset perturbation theory, coupled cluster singles and doubles (CCSD), and CCSD with perturbative triples correction levels show that this newly optimized CIM approach is a reliable theoretical tool for electron correlation calculations of various large chemical systems. © 2018 Wiley Periodicals, Inc.     

Møller-Plesset perturbation energies and distances for HeC(20) extrapolated to the complete basis set limit

The potential energy surface for the C₂₀–He interaction is extrapolated for three representative cuts to the complete basis set limit using second‐order Møller–Plesset perturbation calculations with correlation consistent basis sets up to the doubly augmented variety. The results both with and without counterpoise correction show consistency with each other, supporting that extrapolation without such a correction provides a reliable scheme to elude the basis‐set‐superposition error. Converged attributes are obtained for the C₂₀– He interaction, which are used to predict the fullerene dimer ones. Time requirements show that the method can be drastically more economical than the counterpoise procedure and even competitive with Kohn‐Sham density functional theory for the title system. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2009     

Optimized accurate auxiliary basis sets for RI-MP2 and RI-CC2 calculations for the atoms Rb to Rn

The introduction of the resolution-of-the-identity (RI) approximation for electron repulsion integrals in quantum chemical calculations requires in addition to the orbital basis so-called auxiliary or fitting basis sets. We report here such auxiliary basis sets optimized for second-order Møller–Plesset perturbation theory for the recently published (Weigend and Ahlrichs Phys Chem Chem Phys, 2005, 7, 3297–3305) segmented contracted Gaussian basis sets of split, triple-ζ and quadruple-ζ valence quality for the atoms Rb–Rn (except lanthanides). These basis sets are designed for use in connection with small-core effective core potentials including scalar relativistic corrections. Hereby accurate resolution-of-the-identity calculations with second-order Møller–Plesset perturbation theory (MP2) and related methods can now be performed for molecules containing elements from H to Rn. The error of the RI approximation has been evaluated for a test set of 385 small and medium sized molecules, which represent the common oxidation states of each element, and is compared with the one-electron basis set error, estimated based on highly accurate explicitly correlated MP2–R12 calculations. With the reported auxiliary basis sets the RI error for MP2 correlation energies is typically two orders of magnitude smaller than the one-electron basis set error, independent on the position of the atoms in the periodic table.     

Approximate fourth-order perturbation theory of the electron correlation energy

An approximate fourth-order expression for the electron correlation energy in the Møller–Plesset perturbation scheme is proposed. It takes into account all the contributions to the fourthorder energy neglecting only those of the triple-substituted determinants. It is size consistent and correct to fourth order for an assembly of isolated two-electron systems. Illustrative calculations are reported for a series of small molecules.     

An approximate potential energy surface for HeI2 collisions

We have calculated the interaction potential for HeI₂ in T-shaped geometries using Hartree—Fock and Møller—Plesset third-order perturbation t     

Ab initioMO calculations on the acidities of water and methanol,and hydrogen bond energies of the conjugate ions with a water molecule

The acidities, deprotonation energies, of water and methanol were calculated by the use of the ab initio self-consistent-field (SCF ) molecular orbital (MO ) method with electron correlation computed by the thirdorder Møller–Plesset perturbation method and configuration interaction with double excitations. Zero-point vibrational energy correction translational energy change, and the PV work term were included to evaluate the accurate acidities. The calculated acidity difference including these corrections was 7 kcal/mol, which is somewhat smaller than the experimental ones (9.5–12.5 kcal/mol) recently determined. The hydrogen bond energies of the conjugate ions (OH^? and CH₃O^?) with a water molecule were calculated to be 2.3 kcal/mol near the Hartree–Fock limit; this energy only amounts to 25% of the (total) hydration energy difference between the two negative ions. The aqueous solvation effect on the acidity scale was discussed.     

Charge density by the use of the hiller–sucher–feinberg identity. Perturbative corrections and configuration interaction effects

Configuration interaction (CI ) calculations are carried out for He atom and H₂, LiH, and BH molecules in order to obtain the value of the charge density at the nucleus by the use of the Hiller–Sucher–Feinberg (HSF ) identity. The HSF density also can be calculated with the double perturbation theory based on the Møller–Plesset-type theory and also on the Epstein–Nesbet type. It is found that each value of the HSF density for these perturbative corrections is very similar to the corresponding value by the CI with single- and double-electron excitations. Especially for the He atom, the HSF density value by the full-CI wave function coincides with the exact value of the charge density in three significant figures. The electron correlation effects are found to be small but cannot be ignored in the HSF density. It can be shown that the superiority of the HSF density over the usual delta-function-type density is excellently confirmed.     

Multireference perturbation theory with optimized partitioning. II. Applications to molecular systems

The second‐order multireference perturbation theory using an optimized partitioning, denoted as MROPT(2), is applied to calculations of various molecular properties—excitation energies, spectroscopic parameters, and potential energy curves—for five molecules: ethylene, butadiene, benzene, N₂, and O₂. The calculated results are compared with those obtained with second‐ and third‐order multireference perturbation theory using the traditional partitioning techniques. We also give results from computations using the multireference configuration interaction (MRCI) method. The presented results show very close resemblance between the new method and MRCI with renormalized Davidson correction. The accuracy of the new method is good and is comparable to that of second‐order multireference perturbation theory using Møller‐Plesset partitioning. © 2003 Wiley Periodicals, Inc. J Comput Chem 24: 1390–1400, 2003     

Computerized implementation of higher‐order electron‐correlation methods and their linear‐scaling divide‐and‐conquer extensions

We have implemented a linear‐scaling divide‐and‐conquer (DC)‐based higher‐order coupled‐cluster (CC) and Møller–Plesset perturbation theories (MPPT) as well as their combinations automatically by means of the tensor contraction engine, which is a computerized symbolic algebra system. The DC‐based energy expressions of the standard CC and MPPT methods and the CC methods augmented with a perturbation correction were proposed for up to high excitation orders [e.g., CCSDTQ, MP4, and CCSD(2)_TQ]. The numerical assessment for hydrogen halide chains, polyene chains, and first coordination sphere (C1) model of photoactive yellow protein has revealed that the DC‐based correlation methods provide reliable correlation energies with significantly less computational cost than that of the conventional implementations. © 2017 Wiley Periodicals, Inc.     

Treatment of aperiodicity and correlation in polymers

The negative factor counting method (in its simple and matrix block form) for the determination of the density of states of disordered polymer chains and its applications to different aperiodic organic polymers are reviewed. The problems of the calculation of the correlation energy in large systems like polymers will be discussed. Different ways for the partitioning of an energy band into regions make it possible to perform Møller–Plesset perturbation theoretical calculations on polymers. Applications to hydrogen chains are presented. Finally, possible applications to disordered polymers are also discussed.     

Effect of electron correlation on theoretical vibrational frequencies

Theoretical HF /6-31G * (Hartree–Fock, 6-31G * basis set) and MP 2/6-31G * (second-order Møller–Plesset, 6-31G * basis set) vibrational frequencies based on complete quadratic force fields have been obtained for a set of 36 one- and two-heavy-atom molecules comprising first-row elements for which experimental spectroscopic data are available. Frequencies calculated at the HF /6-31G * level are an average of 12.6% higher than experimental values. Partial treatment of electron correlation via the perturbation method of Møller and Plesset, terminated at second order, leads to a significant reduction in this error, although theoretical MP 2/6-31G * frequencies are still larger than the experimental quantities by 7.3%. Part of the difference may be traced to the restriction of quadratic force fields, as comparison with experimental harmonic frequencies shows deviations of only 9.5% and 4.7% for the two levels, respectively. The calculated frequencies are used in conjunction with the corresponding theoretical equilibrium structures to obtain absolute molecular entropies, which may in turn be used to yield entropies of reaction. These latter quantities are generally in good accord with entropies derived using experimental structures and frequencies.     

Implementation of dynamical nucleation theory with quantum potentials

A method is implemented within the context of dynamical nucleation theory in order to efficiently determine the ab initio water dimer evaporation rate constant. The drive for increased efficiency in a Monte Carlo methodology is established by the need to use relatively expensive quantum mechanical interaction potentials. A discussion is presented illustrating the theory, algorithm, and implementation of this method to the water dimer. Hartree–Fock and second order Møller–Plesset perturbation theories along with the Dang–Chang polarizable classical potential are utilized to determine the ab initio water dimer evaporation rate constant. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2009     

A general efficient implementation of the BSSE-free SCF and MP2 methods based on the chemical Hamiltonian approach

We describe some details related to a new, general, and efficient implementation of the BSSE‐free SCF and second‐order Møller–Plesset perturbation theories of intermolecular interactions, based on the “Chemical Hamiltonian Approach” (CHA). The program is applicable for both open‐shell and closed‐shell systems and for an arbitrary number of interacting subsystems. With the new program the CHA method is faster than the usual “counterpoise correction” scheme for single point calculations, especially for clusters consisting of several molecules. The numerical results provided by these conceptually different schemes, however, have again found to be very close to each other. The CHA scheme is particularly good for providing truly BSSE‐free MP2 data for intermolecular potentials. © 2006 Wiley Periodicals, Inc. J Comput Chem 27: 1505–1516, 2006