Application of Gaussian-type geminals in local second-order Møller-Plesset perturbation theory

In this work Gaussian-type Geminals (GTGs) are applied in local second-order Moller-Plesset perturbation theory to improve the basis set convergence. Our implementation is based on the weak orthogonality functional of Szalewicz et al., [Chem. Phys. Lett. 91, 169 (1982); J. Chem. Phys. 78, 1420 (1983)] and a newly developed program for calculating the necessary many-electron integrals. The local approximations together with GTGs in the treatment of the correlation energy are introduced and tested. First results for correlation energies of H(2)O, CH(4), CO, C(2)H(2), C(2)H(4), H(2)CO, and N(2)H(4) as well as some reaction and activation energies are presented. More than 97% of the valence-shell correlation energy is recovered using aug-cc-pVDZ basis sets and six GTGs per electron pair. The results are compared with conventional calculations using correlation-consistent basis sets as well as with MP2-R12 results.     

Improved third-order Møller-Plesset perturbation theory

Based on a partitioning of the total correlation energy into contributions from parallel‐ and antiparallel‐spin pairs of electrons, a modified third‐order Møller–Plesset (MP) perturbation theory is developed. The method, termed SCS–MP3 (SCS for spin‐component‐scaled) continues previous work on an improved version of MP2 (S. Grimme, J Chem Phys 2003, 118, 9095). A benchmark set of 32 isogyric reaction energies, 11 atomization energies, and 11 stretched geometries is used to assess to performance of the model in comparison to the standard quantum chemical approaches MP2, MP3, and QCISD(T). It is found, that the new method performs significantly better than usual MP2/MP3 and even outperforms the more costly QCISD method. Opposite to the usual MP series, the SCS third‐order correction uniformly improves the results. Dramatic enhancements are especially observed for the more difficult atomization energies, some of the stretched geometries, and reaction and ionization energies involving transition metal compounds where the method seems to be competitive or even superior to the widely used density functional approaches. Further tests performed for other complex systems (biradicals, C₂₀ isomers, transition states) demonstrate that the SCS–MP3 model yields often results of QCISD(T) accuracy. The uniformity with which the new approach improves for very different correlation problems indicates significant robustness, and suggests it as a valuable quantum chemical method of general use. © 2003 Wiley Periodicals, Inc. J Comput Chem 24: 1529–1537, 2003     

Tensor factorizations of local second-order Møller-Plesset theory

Efficient electronic structure methods can be built around efficient tensor representations of the wavefunction. Here we first describe a general view of tensor factorization for the compact representation of electronic wavefunctions. Next, we use this language to construct a low-complexity representation of the doubles amplitudes in local second-order M?ller-Plesset perturbation theory. We introduce two approximations--the direct orbital-specific virtual approximation and the full orbital-specific virtual approximation. In these approximations, each occupied orbital is associated with a small set of correlating virtual orbitals. Conceptually, the representation lies between the projected atomic orbital representation in Pulay-Saeb? local correlation theories and pair natural orbital correlation theories. We have tested the orbital-specific virtual approximations on a variety of systems and properties including total energies, reaction energies, and potential energy curves. Compared to the Pulay-Saeb? ansatz, we find that these approximations exhibit favorable accuracy and computational times while yielding smooth potential energy curves.     

Explicitly correlated second-order Møller-Plesset perturbation theory employing pseudospectral numerical quadratures

We implemented explicitly correlated second-order M?ller-Plesset perturbation theory with numerical quadratures using pseudospectral construction of grids. Introduction of pseudospectral approach for the calculation of many-electron integrals gives a possibility to use coarse grids without significant loss of precision in correlation energies, while the number of points in the grid is reduced about nine times. The use of complementary auxiliary basis sets as the sets of dealiasing functions is justified at both theoretical and computational levels. Benchmark calculations for a set of 16 molecules have shown the possibility to keep an error of second-order correlation energies within 1 milihartree (mH) with respect to MP2-F12 method with dense grids. Numerical tests for a set of 13 isogyric reactions are also performed.     

An atomic orbital-based reformulation of energy gradients in second-order Møller-Plesset perturbation theory

A fully atomic orbital (AO)-based reformulation of second-order M?ller-Plesset perturbation theory (MP2) energy gradients is introduced, which provides the basis for reducing the computational scaling with the molecular size from the fifth power to linear. Our formulation avoids any transformation between the AO and the molecular orbital (MO) basis and employs pseudodensity matrices similar to the AO-MP2 energy expressions within the Laplace scheme for energies. The explicit computation of perturbed one-particle density matrices emerging in the new AO-based gradient expression is avoided by reformulating the Z-vector method of Handy and Schaefer [J. Chem. Phys. 81, 5031 (1984)] within a density matrix-based scheme.     

Application of second-order Møller-Plesset perturbation theory with resolution-of-identity approximation to periodic systems

Efficient periodic boundary condition (PBC) calculations by the second-order M?ller-Plesset perturbation (MP2) method based on crystal orbital formalism are developed by introducing the resolution-of-identity (RI) approximation of four-center two-electron repulsion integrals (ERIs). The formulation and implementation of the PBC RI-MP2 method are presented. In this method, the mixed auxiliary basis functions of the combination of Poisson and Gaussian type functions are used to circumvent the slow convergence of the lattice sum of the long-range ERIs. Test calculations of one-dimensional periodic trans-polyacetylene show that the PBC RI-MP2 method greatly reduces the computational times as well as memory and disk sizes, without the loss of accuracy, compared to the conventional PBC MP2 method.     

New implementation of second-order Møller-Plesset perturbation theory with an analytic Slater-type geminal

The author introduces a new method for the exchange commutator integrals in explicitly correlated M?ller-Plesset second order perturbation theory. The method is well suited with an analytic Slater-type geminal correlation factor. He also explains the scheme for auxiliary integrals needed for the correlation factor. Based on different Ans?tze, he analyzes the performance of the method on correlation energies and reaction enthalpies in detail.     

Distributed memory parallel implementation of energies and gradients for second-order Møller-Plesset perturbation theory with the resolution-of-the-identity approximation

We present a parallel implementation of second-order M?ller-Plesset perturbation theory with the resolution-of-the-identity approximation (RI-MP2). The implementation is based on a recent improved sequential implementation of RI-MP2 within the Turbomole program package and employs the message passing interface (MPI) standard for communication between distributed memory nodes. The parallel implementation extends the applicability of canonical MP2 to considerably larger systems. Examples are presented for full geometry optimizations with up to 60 atoms and 3300 basis functions and MP2 energy calculations with more than 200 atoms and 7000 basis functions.     

Local explicitly correlated second-order Møller-Plesset perturbation theory with pair natural orbitals

We explore using a pair natural orbital analysis of approximate first-order pair functions as means to truncate the space of both virtual and complementary auxiliary orbitals in the context of explicitly correlated F12 methods using localised occupied orbitals. We demonstrate that this offers an attractive procedure and that only 10-40 virtual orbitals per significant pair are required to obtain second-order valence correlation energies to within 1-2% of the basis set limit. Moreover, for this level of virtual truncation, only 10-40 complementary auxiliary orbitals per pair are required for an accurate resolution of the identity in the computation of the three- and four-electron integrals that arise in explicitly correlated methods.     

General biorthogonal projected bases as applied to second-order Møller-Plesset perturbation theory

With low-order scaling correlated wave function theories in mind, we present second quantization formalism as well as biorthonormalization procedures for general--singular or nonsingular--bases. Of particular interest are the so-called projected atomic orbital bases, which are obtained from a set of atom-centered functions and feature a separation of occupied and virtual spaces. We demonstrate the formalism by deriving and implementing second-order M?ller-Plesset perturbation theory in it, and discuss the convergence and preconditioning of the iterative amplitude equations in detail.     

Intruder state avoidance multireference Møller-Plesset perturbation theory

A new perturbation approach is proposed that enhances the low‐order, perturbative convergence by modifying the zeroth‐order Hamiltonian in a manner that enlarges any small‐energy denominators that may otherwise appear in the perturbative expansion. This intruder state avoidance (ISA) method can be used in conjunction with any perturbative approach, but is most applicable to cases where small energy denominators arise from orthogonal‐space states—so‐called intruder states—that should, under normal circumstances, make a negligible contribution to the target state of interests. This ISA method is used with multireference Møller–Plesset (MRMP) perturbation theory on potential energy curves that are otherwise plagued by singularities when treated with (conventional) MRMP; calculation are performed on the 1³Σ state of O₂; and the 2¹Δ, 3¹Δ, 2³Δ, and 3³Δ states of AgH. This approach is also applied to other calculations where MRMP is influenced by intruder states; calculations are performed on the ³Π_u state of N₂, the ³Π state of CO, and the 2¹A′ state of formamide. A number of calculations are also performed to illustrate that this approach has little or no effect on MRMP when intruder states are not present in perturbative calculations; vertical excitation energies are computed for the low‐lying states of N₂, C₂, CO, formamide, and benzene; the adiabatic ¹A₁–³B₁ energy separation in CH₂, and the spectroscopic parameters of O₂ are also calculated. Vertical excitation energies are also performed on the Q and B bands states of free‐base, chlorin, and zinc–chlorin porphyrin, where somewhat larger couplings exists, and—as anticipated—a larger deviation is found between MRMP and ISA‐MRMP. © 2002 Wiley Periodicals, Inc. J Comput Chem 10: 957–965, 2002     

A kinetic energy fitting metric for resolution of the identity second-order Møller-Plesset perturbation theory

A kinetic-energy-based fitting metric for application in the context of resolution of the identity second-order M?ller-Plesset perturbation theory is presented, which is derived from the Poisson equation. Preliminary tests of the applicability include the evaluation of the error in the correlation energy, compared to standard M?ller-Plesset perturbation theory, with respect to the auxiliary basis set employed. We comment on the potential merits of this fitting metric, compared to standard resolution of the identity second-order M?ller-Plesset perturbation theory, and discuss its scaling behavior in the limit of large molecules.     

Optimization of orbital-specific virtuals in local Møller-Plesset perturbation theory

We present an orbital-optimized version of our orbital-specific-virtuals second-order M?ller-Plesset perturbation theory (OSV-MP2). The OSV model is a local correlation ansatz with a small basis of virtual functions for each occupied orbital. It is related to the Pulay-Saeb? approach, in which domains of virtual orbitals are drawn from a single set of projected atomic orbitals; but here the virtual functions associated with a particular occupied orbital are specifically tailored to the correlation effects in which that orbital participates. In this study, the shapes of the OSVs are optimized simultaneously with the OSV-MP2 amplitudes by minimizing the Hylleraas functional or approximations to it. It is found that optimized OSVs are considerably more accurate than the OSVs obtained through singular value decomposition of diagonal blocks of MP2 amplitudes, as used in our earlier work. Orbital-optimized OSV-MP2 recovers smooth potential energy surfaces regardless of the number of virtuals. Full optimization is still computationally demanding, but orbital optimization in a diagonal or Kapuy-type MP2 approximation provides an attractive scheme for determining accurate OSVs.     

Analytical energy gradients for local second-order Møller-Plesset perturbation theory using density fitting approximations

An efficient method to compute analytical energy derivatives for local second-order M?ller-Plesset perturbation energy is presented. Density fitting approximations are employed for all 4-index integrals and their derivatives. Using local fitting approximations, quadratic scaling with molecular size and cubic scaling with basis set size for a given molecule is achieved. The density fitting approximations have a negligible effect on the accuracy of optimized equilibrium structures or computed energy differences. The method can be applied to much larger molecules and basis sets than any previous second-order M?ller-Plesset gradient program. The efficiency and accuracy of the method is demonstrated for a number of organic molecules as well as for molecular clusters. Examples of geometry optimizations for molecules with 100 atoms and over 2000 basis functions without symmetry are presented.     

Two-level hierarchical parallelization of second-order Møller-Plesset perturbation calculations in divide-and-conquer method

A two‐level hierarchical parallelization scheme including the second‐order Møller–Plesset perturbation (MP2) theory in the divide‐and‐conquer method is presented. The scheme is a combination of coarse‐grain parallelization assigning each subsystem to a group of processors, with fine‐grain parallelization, where the computational tasks for evaluating MP2 correlation energy of the assigned subsystem are distributed among processors in the group. Test calculations demonstrate that the present scheme shows high parallel efficiency and makes MP2 calculations practical for very large molecules. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Analysis of self-consistency effects in range-separated density-functional theory with Møller-Plesset perturbation theory

Range-separated density-functional theory combines wave function theory for the long-range part of the two-electron interaction with density-functional theory for the short-range part. When describing the long-range interaction with non-variational methods, such as perturbation or coupled-cluster theories, self-consistency effects are introduced in the density functional part, which for an exact solution requires iterations. They are generally assumed to be small but no detailed study has been performed so far. Here, the authors analyze self-consistency when using M?ller-Plesset-type (MP) perturbation theory for the long range interaction. The lowest-order self-consistency corrections to the wave function and the energy, that enter the perturbation expansions at the second and fourth order, respectively, are both expressed in terms of the one-electron reduced density matrix. The computational implementation of the latter is based on a Neumann series which, interestingly, even though the effect is small, usually diverges. A convergence technique, which perhaps can be applied in other uses of Neumann series in perturbation theory, is proposed. The numerical results thus obtained show that, in weakly bound systems, self-consistency can be neglected since the long-range correlation does not affect the density significantly. Although MP is not adequate for multireference systems, it can still be used as a reliable analysis tool. Though the density change is not negligible anymore in such cases, self-consistency effects are found to be much smaller than long-range correlation effects (less than 10% for the systems considered). For that reason, a sensible approximation might be to update the short-range energy functional term while freezing its functional derivative, namely, the short-range local potential, in the wave function optimization. The accuracy of such an approximation still needs to be assessed.     

Analytic energy gradient in combined second-order Møller-Plesset perturbation theory and polarizable force field calculation

Second-order M?ller-Plesset perturbation theory (MP2) is used to describe electronic correlation on the basis of Hartree-Fock (HF) variational calculations that incorporate induced dipole polarizable force fields (i.e., QM/MMpol style HF and MP2). The Z-vector equations for regular closed shell and open shell MP2 methods (RMP2, ZAPT2, and UMP2) are extended to include induced dipole contributions to determine the MP2 response density so that nuclear gradient and other properties can be efficiently evaluated. A better estimation of the induced dipole polarization energy can be obtained using the MP2 relaxed density. QM/MMpol style MP2 molecular dynamics simulations are performed for the ground state and first triplet state of acetone solvated by 1024 polarizable water molecules. A switching function is used to ensure energy conservation in QM/MM simulation under periodic boundary condition.     

Calculation of weakly polar interaction energies in polypeptides using density functional and local Møller-Plesset perturbation theory

The interaction energies of ubiquitous weakly polar interactions in proteins are comparable with those of hydrogen bonds, consequently, they stabilize local, secondary, and tertiary structures. However, the most widely-used density functionals fail to describe the weakly polar interactions. Thus, it is important to find and test functionals which adequately describe and quantify the energetics of such interactions. For this purpose, interaction energies in the hydrophobic core of rubredoxin (PDB id: 1rb9) and in the S22 subset of the JSCH-2005 benchmark database were computed with the BHandHLYP and PWPW91 functionals and with the pseudospectral implementation of the local MP2 (PS-LMP2) method. The cc-pVDZ, cc-pVTZ(-f), cc-pVTZ, cc-pVQZ(-g), aug-cc-pVDZ, aug-cc-VTZ(-f), and aug-cc-pVTZ basis sets were used for the calculations. In the S22 subset the PS-LMP2 results were extrapolated to the complete basis set limit. Furthermore, the a posteriori counterpoise method of Boys and Bernardi was used to correct the basis set superposition errors in the calculation of interaction energies. Calculations using the BHandHLYP functional, both for the various weakly polar interactions in rubredoxin and for the dispersion interactions in the S22 subset, were in good agreement with those using the coupled cluster (CCSD(T)) and the resolution of identity MP2 (RIMP2) methods and clearly outperformed both the PWPW91 functional and the PS-LMP2 method. The results for the S22 hydrogen bonded subset, obtained with PWPW91 calculations, were closest to those of the reference high level calculations. For the "mixed" (hydrogen bonded and dispersive) interactions in the S22 subset, results obtained with the BHandHLYP and PS-LMP2 calculations agreed well with the reference calculations.     

Variational second-order Møller-Plesset theory based on the Luttinger-Ward functional

In recent years there have been some rather successful applications of a new variational technique for calculating the total energies of electronic systems. The new method is based on many-body perturbation theory and uses the one-electron Green function as the basic "variable" rather than the wave function of traditional variational calculations. It is the purpose of the present work to promote the new methods within the realm of traditional theoretical chemistry by demonstrating their utility for calculating the correlation energies of a number of atoms at a level corresponding to second-order M?ller-Plesset perturbation theory. The generalization to any desired order of perturbation theory is not hard to accomplish.     

Analytic energy gradient for second-order Møller-Plesset perturbation theory based on the fragment molecular orbital method

The first derivative of the total energy with respect to nuclear coordinates (the energy gradient) in the fragment molecular orbital (FMO) method is applied to second order M?ller-Plesset perturbation theory (MP2), resulting in the analytic derivative of the correlation energy in the external self-consistent electrostatic field. The completely analytic energy gradient equations are formulated at the FMO-MP2 level. Both for molecular clusters (H(2)O)(64) and a system with fragmentation across covalent bonds, a capped alanine decamer, the analytic FMO-MP2 energy gradients with the electrostatic dimer approximation are shown to be complete and accurate by comparing them with the corresponding numeric gradients. The developed gradient is parallelized with the parallel efficiency of about 97% on 32 Pentium4 nodes connected by Gigabit Ethernet.