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.     

Second order Møller-Plesset perturbation theory based upon the fragment molecular orbital method

The fragment molecular orbital (FMO) method was combined with the second order M?ller-Plesset (MP2) perturbation theory. The accuracy of the method using the 6-31G(*) basis set was tested on (H(2)O)(n), n=16,32,64; alpha-helices and beta-strands of alanine n-mers, n=10,20,40; as well as on (H(2)O)(n), n=16,32,64 using the 6-31 + + G(**) basis set. Relative to the regular MP2 results that could be afforded, the FMO2-MP2 error in the correlation energy did not exceed 0.003 a.u., the error in the correlation energy gradient did not exceed 0.000 05 a.u./bohr and the error in the correlation contribution to dipole moment did not exceed 0.03 debye. An approximation reducing computational load based on fragment separation was introduced and tested. The FMO2-MP2 method demonstrated nearly linear scaling and drastically reduced the memory requirements of the regular MP2, making possible calculations with several thousands basis functions using small Pentium clusters. As an example, (H(2)O)(64) with the 6-31 + + G(**) basis set (1920 basis functions) can be run in 1 Gbyte RAM and it took 136 s on a 40-node Pentium4 cluster.     

Analytic energy gradients in combined second order Møller-Plesset perturbation theory and conductorlike polarizable continuum model calculation

The analytic energy gradients in combined second order M?ller-Plesset perturbation theory and conductorlike polarizable continuum model calculations are derived and implemented for spin-restricted closed shell (RMP2), Z-averaged spin-restricted open shell (ZAPT2), and spin-unrestricted open shell (UMP2) cases. Using these methods, the geometries of the S(0) ground state and the T(1) state of three nucleobase pairs (guanine-cytosine, adenine-thymine, and adenine-uracil) in the gas phase and aqueous solution phase are optimized. It is found that in both the gas phase and the aqueous solution phase the hydrogen bonds in the T(1) state pairs are weakened by ~1 kcal/mol as compared to those in the S(0) state pairs.     

Accuracy of the three-body fragment molecular orbital method applied to Møller-Plesset perturbation theory

The three‐body energy expansion in the fragment molecular orbital method (FMO) was applied to the 2nd order Møller–Plesset theory (MP2). The accuracy of both the two and three‐body expansions was determined for water clusters, alanine n‐mers (α‐helices and β‐strands) and one synthetic protein, using the 6‐31G* and 6‐311G* basis sets. At the best level of theory (three‐body, two molecules/residues per fragment), the absolute errors in energy relative to ab initio MP2 were at most 1.2 and 5.0 mhartree, for the 6‐31G* and 6‐311G* basis sets, respectively. The relative accuracy was at worst 99.996% and 99.96%, for 6‐31G* and 6‐311G*, respectively. A three‐body approximation was introduced and the optimum threshold value was determined. The protein calculation (6‐31G*) at the production level (FMO2/2) took 3 h on 36 3.2‐GHz Pentium 4 nodes and had the absolute error in the MP2 correlation energy of only 2 kcal/mol. © 2007 Wiley Periodicals, Inc. J Comput Chem, 2007     

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.     

A parallelized integral-direct second-order Møller–Plesset perturbation theory method with a fragment molecular orbital scheme

We propose a parallelized integral-direct algorithm of the second-order Møller–Plesset perturbation theory (MP2) as a size-consistent correlated method. The algorithm is a modification of the recipe by Mochizuki et al. [(1996) Theor Chim Acta 93:211]. There is no need to communicate the bulky data of integrals across worker processes, keeping the formal fifth-power dependence on the number of basis functions. A multiple integral screening procedure is incorporated to reduce the operation costs effectively. An approximate MP2 density matrix can also be directly calculated through the integral contraction with orbital energies. We implement the MP2 code by accepting Kitauras fragment molecular orbital (FMO) scheme as in the program ABINIT-MP developed by Nakano et al. [(2002) Chem Phys Lett 351:475]. The error in the FMO–MP2 energies is found to be within the order of the chemical accuracy. Timing and parallel acceleration results are shown for test molecules.     

Alternative linear-scaling methodology for the second-order Møller-Plesset perturbation calculation based on the divide-and-conquer method

A new scheme for obtaining the approximate correlation energy in the divide-and-conquer (DC) method of Yang [Phys. Rev. Lett. 66, 1438 (1991)] is presented. In this method, the correlation energy of the total system is evaluated by summing up subsystem contributions, which are calculated from subsystem orbitals based on a scheme for partitioning the correlation energy. We applied this method to the second-order Moller-Plesset perturbation theory (MP2), which we call DC-MP2. Numerical assessment revealed that this scheme provides a reliable correlation energy with significantly less computational cost than the conventional MP2 calculation.     

Resolution of the identity atomic orbital Laplace transformed second order Møller-Plesset theory for nonconducting periodic systems

An improvement in performance of the atomic orbital Laplace transformed second-order M?ller-Plesset (AO-LT-MP2) method for periodic systems is reported using the resolution of identity (RI) technique. Transformation of the two-electron integrals constitutes the main computational bottleneck of the AO-LT-MP2 method. A substitution of regular four-center integrals by their three center counterparts in the RI approximation naturally reduces the computational cost of the integral transformation step. The RI divergence problem in the presence of periodic boundary conditions is solved in our implementation by restricting the fitting domain. Accuracy and computational efficiency of the RI-AO-LT-MP2 approach are assessed on a set of one-dimensional test systems: trans-polyacetylene and anti-transoid polymethineimine.     

An efficient atomic orbital based second-order Møller-Plesset gradient program

Based on the orbital-invariant atomic orbital formulation of the MP2 (M?ller-Plesset second-order perturbation theory) energy and gradient [P. Pulay and S. Saeb?, Theor. Chim. Acta 69, 357 (1986)], we have derived and programmed detailed working equations for closed-shell MP2 gradients. The orbital-invariant form avoids the difficulties of other formulations with frozen orbitals, and allows the use of arbitrary occupied orbitals, an important consideration for local correlation theories, although the present program uses canonical molecular orbitals. The atomic orbital formulation offers savings both in storage and computer time. Test calculations on systems containing up to approximately 100 atoms and approximately 1000 basis functions, performed on a single personal computer, are reported. Parallelization of the code is underway.     

Theoretical interpretation of Grimme's spin-component-scaled second order Møller-Plesset theory

It is shown that spin-component-scaled second order M?ller-Plesset theory proposed by Grimme [J. Chem. Phys. 118, 9095 (2003)] can be interpreted as a two-parameter scaling of the zero order Hamiltonian, a generalization of the approach reported by Feenberg [Phys. Rev. 103, 1116 (1956)].     

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.     

Quadratically convergent algorithm for orbital optimization in the orbital-optimized coupled-cluster doubles method and in orbital-optimized second-order Møller-Plesset perturbation theory

Using a Lagrangian-based approach, we present a more elegant derivation of the equations necessary for the variational optimization of the molecular orbitals (MOs) for the coupled-cluster doubles (CCD) method and second-order M?ller-Plesset perturbation theory (MP2). These orbital-optimized theories are referred to as OO-CCD and OO-MP2 (or simply "OD" and "OMP2" for short), respectively. We also present an improved algorithm for orbital optimization in these methods. Explicit equations for response density matrices, the MO gradient, and the MO Hessian are reported both in spin-orbital and closed-shell spin-adapted forms. The Newton-Raphson algorithm is used for the optimization procedure using the MO gradient and Hessian. Further, orbital stability analyses are also carried out at correlated levels. The OD and OMP2 approaches are compared with the standard MP2, CCD, CCSD, and CCSD(T) methods. All these methods are applied to H(2)O, three diatomics, and the O(4)(+) molecule. Results demonstrate that the CCSD and OD methods give nearly identical results for H(2)O and diatomics; however, in symmetry-breaking problems as exemplified by O(4)(+), the OD method provides better results for vibrational frequencies. The OD method has further advantages over CCSD: its analytic gradients are easier to compute since there is no need to solve the coupled-perturbed equations for the orbital response, the computation of one-electron properties are easier because there is no response contribution to the particle density matrices, the variational optimized orbitals can be readily extended to allow inactive orbitals, it avoids spurious second-order poles in its response function, and its transition dipole moments are gauge invariant. The OMP2 has these same advantages over canonical MP2, making it promising for excited state properties via linear response theory. The quadratically convergent orbital-optimization procedure converges quickly for OMP2, and provides molecular properties that are somewhat different than those of MP2 for most of the test cases considered (although they are similar for H(2)O). Bond lengths are somewhat longer, and vibrational frequencies somewhat smaller, for OMP2 compared to MP2. In the difficult case of O(4)(+), results for several vibrational frequencies are significantly improved in going from MP2 to OMP2.     

The divide-expand-consolidate family of coupled cluster methods: numerical illustrations using second order Møller-Plesset perturbation theory

Previously, we have introduced the linear scaling coupled cluster (CC) divide-expand-consolidate (DEC) method, using an occupied space partitioning of the standard correlation energy. In this article, we show that the correlation energy may alternatively be expressed using a virtual space partitioning, and that the Lagrangian correlation energy may be partitioned using elements from both the occupied and virtual partitioning schemes. The partitionings of the correlation energy leads to atomic site and pair interaction energies which are term-wise invariant with respect to an orthogonal transformation among the occupied or the virtual orbitals. Evaluating the atomic site and pair interaction energies using local orbitals leads to a linear scaling algorithm and a distinction between Coulomb hole and dispersion energy contributions to the correlation energy. Further, a detailed error analysis is performed illustrating the error control imposed on all components of the energy by the chosen energy threshold. This error control is ultimately used to show how to reduce the computational cost for evaluating dispersion energy contributions in DEC.     

Influence of the water molecule on cation-pi interaction: ab initio second order Møller-Plesset perturbation theory (MP2) calculations

The influence of introducing water molecules into a cation-pi complex on the interaction between the cation and the pi system was investigated using the MP2/6-311++G method to explore how a cation-pi complex changes in terms of both its geometry and its binding strength during the hydration. The calculation on the methylammonium-benzene complex showed that the cation-pi interaction is weakened by introducing H(2)O molecules into the system. For example, the optimized interaction distance between the cation and the benzene becomes longer and longer, the transferred charge between them becomes less and less, and the cation-pi binding strength becomes weaker and weaker as the water molecule is introduced one by one. Furthermore, the introduction of the third water molecule leads to a dramatic change in both the complex geometry and the binding energy, resulting in the destruction of the cation-pi interaction. The decomposition on the binding energy shows that the influence is mostly brought out through the electrostatic and induction interactions. This study also demonstrated that the basis set superposition error, thermal energy, and zero-point vibrational energy are significant and needed to be corrected for accurately predicting the binding strength in a hydrated cation-pi complex at the MP2/6-311++G level. Therefore, the results are helpful to better understand the role of water molecules in some biological processes involving cation-pi interactions.     

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.     

An explicitly correlated second order Møller-Plesset theory using a frozen Gaussian geminal

A variant of the MP2-R12 class of theories is introduced using an arbitrary geminal function in the place of r12. Integrals are derived for the case where the geminal is expanded in a basis of Gaussian functions in the interelectronic distance. Recurrence relations are derived that do not depend on the exponents of the Gaussian geminals, allowing much of the integration work to be performed after summations over the geminal expansion. Sample calculations at various levels of explicitly correlated MP2 theory are presented for He, Ne, and water.     

Intermolecular potentials of the methane dimer calculated with Møller-Plesset perturbation theory and density functional theory

We have calculated the intermolecular interaction potentials of the methane dimer at the minimum-energy D(3d) conformation using the Hartree-Fock (HF) self-consistent theory, the correlation-corrected second-order M?ller-Plesset (MP2) perturbation theory, and the density functional theory (DFT) with the Perdew-Wang (PW91) functional as the exchange or the correlation part. The HF calculations yield unbound potentials largely due to the exchange-repulsion interaction. In the MP2 calculations, the basis set effects on the repulsion exponent, the equilibrium bond length, the binding energy, and the asymptotic behavior of the calculated intermolecular potentials have been thoroughly studied. We have employed basis sets from the Slater-type orbitals fitted with Gaussian functions (STO-nG) (n=3-6) [Quantum Theory of Molecular and Solids: The Self-Consistent Field for Molecular and Solids (McGraw-Hill, New York, 1974), Vol. 4], Pople's medium size basis sets of Krishnan et al. [J. Chem. Phys. 72, 650 (1980)] [up to 6-311++G(3df,3pd)] to Dunning's correlation consistent basis sets [J. Chem. Phys. 90, 1007 (1989)] (cc-pVXZ and aug-cc-pVXZ) (X=D, T, and Q). With increasing basis size, the repulsion exponent and the equilibrium bond length converge at the 6-31G** basis set and the 6-311++G(2d,2p) basis set, respectively, while a large basis set (aug-cc-pVTZ) is required to converge the binding energy at a chemical accuracy (approximately 0.01 kcal/mol). Up to the largest basis set used, the asymptotic dispersion coefficient has not converged to the destined C6 value from molecular polarizability calculations. The slow convergence could indicate the inefficacy of using the MP2 calculations with Gaussian-type functions to model the asymptotic behavior. Both the basis set superposition error (BSSE) corrected and uncorrected results are presented to emphasize the importance of including such corrections. Only the BSSE corrected results systematically converge to the destined potential curve with increasing basis size. The DFT calculations generate a wide range of interaction patterns, from purely unbound to strongly bound, underestimating or overestimating the binding energy. The binding energy calculated using the PW91PW91 functional and the equilibrium bond length calculated using the PW91VP86 functional are close to the MP2 results at the basis set limit.     

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.     

Scaled opposite spin second order Møller-Plesset theory with improved physical description of long-range dispersion interactions

Separate scaling of the same-spin and opposite spin contributions to the second-order M?ller-Plesset energy can yield statistically improved performance for a variety of chemical problems. If only the opposite spin contribution is scaled, it is also possible to reduce the computational complexity from fifth order to fourth order in system size, with very little degradation of the results. However neither of these scaled MP2 energies recovers the full MP2 result for the dispersion energy of nonoverlapping systems. This deficiency is addressed in this work by using a distance-dependent scaling of the opposite spin correlation energy. The resulting method is compared against the previously proposed scaled MP2 methods on a range of problems involving both short and long-range interactions.     

Calculating stacking interactions in nucleic acid base-pair steps using spin-component scaling and local second order Møller-Plesset perturbation theory

Stacking interaction energies for ten B-DNA base-pair steps are computed with density fitted local second-order M?ller-Plesset perturbation theory (DF-LMP2), and with the spin-component scaled (SCS) and spin-component scaled for nucleobases (SCSN) variants of DF-LMP2. Comparison with existing CBS(T) reference data indicates larger than expected energy differences for both SCS variants. After an analysis of the errors involved, an alternative method of producing reference data is proposed where DF-LMP2/aug-cc-pVTZ and DF-LMP2/aug-cc-pVQZ energies for the whole complex are extrapolated to produce interaction energies that do not require many-body correction and show reduced error in estimation of the basis set limit. A literature correction term from coupled cluster theory with perturbative triples is then added to the DF-LMP2 estimated basis set limit. These new reference data are consistently around 1 kcal mol(-1) less than previous literature data. DF-SCSN-LMP2/aug-cc-pVTZ is found to reproduce the new reference interaction energies with a root mean square error (RMSE) of 0.71 kcal mol(-1), while SCS consistently underestimates the binding energy.