A comparison of linear and nonlinear correlation factors for basis set limit Møller-Plesset second order binding energies and structures of He2, Be2, and Ne2

2006 ◽  
Vol 125 (9) ◽  
pp. 094302 ◽  
Author(s):  
David P. Tew ◽  
Wim Klopper
Keyword(s):  
Second Order ◽  
Binding Energies ◽  
Basis Set ◽  
Nonlinear Correlation ◽  
Linear And Nonlinear ◽  
Moller Plesset ◽  
Correlation Factors

Divergence of Many-Body Perturbation Theory for Noncovalent Interactions of Large Molecules

2019 ◽  
Author(s):  
Brian Nguyen ◽  
Guo P Chen ◽  
Matthew M. Agee ◽  
Asbjörn M. Burow ◽  
Matthew Tang ◽  
...  
Keyword(s):  
Perturbation Theory ◽  
System Size ◽  
State Density ◽  
Second Order ◽  
Binding Energies ◽  
Basis Set ◽  
Many Body ◽  
Many Body Perturbation Theory ◽  
Relative Errors ◽  
Ground State Density

Prompted by recent reports of large errors in noncovalent interaction (NI) energies obtained from many-body perturbation theory (MBPT), we compare the performance of second-order Møller–Plesset MBPT (MP2), spin-scaled MP2, dispersion-corrected semilocal density functional approximations (DFA), and the post-Kohn–Sham random phase approximation (RPA) for predicting binding energies of supramolecular complexes contained in the S66, L7, and S30L benchmarks. All binding energies are extrapolated to the basis set limit, corrected for basis set superposition errors, and compared to reference results of the domain-based local pair-natural orbital coupled-cluster (DLPNO-CCSD(T)) or better quality. Our results confirm that MP2 severely overestimates binding energies of large complexes, producing relative errors of over 100% for several benchmark compounds. RPA relative errors consistently range between 5-10%, significantly less than reported previously using smaller basis sets, whereas spin-scaled MP2 methods show limitations similar to MP2, albeit less pronounced, and empirically dispersion-corrected DFAs perform almost as well as RPA. Regression analysis reveals a systematic increase of relative MP2 binding energy errors with the system size at a rate of approximately 1‰ per valence electron, whereas the RPA and dispersion-corrected DFA relative errors are virtually independent of the system size. These observations are corroborated by a comparison of computed rotational constants of organic molecules to gas-phase spectroscopy data contained in the ROT34 benchmark. To analyze these results, an asymptotic adiabatic connection symmetry-adapted perturbation theory (AC-SAPT) is developed which uses monomers at full coupling whose ground-state density is constrained to the ground-state density of the complex. Using the fluctuation–dissipation theorem, we obtain a nonperturbative “screened second-order” expression for the dispersion energy in terms of monomer quantities which is exact for non-overlapping subsystems and free of induction terms; a first-order RPA-like approximation to the Hartree, exchange, and correlation kernel recovers the macroscopic Lifshitz limit. The AC-SAPT expansion of the interaction energy is obtained from Taylor expansion of the coupling strength integrand. Explicit expressions for the convergence radius of the AC-SAPT series are derived within RPA and MBPT and numerically evaluated. Whereas the AC-SAPT expansion is always convergent for nondegenerate monomers when RPA is used, it is found to spuriously diverge for second-order MBPT, except for the smallest and least polarizable monomers. The divergence of the AC-SAPT series within MBPT is numerically confirmed within RPA; prior numerical results on the convergence of the SAPT expansion for MBPT methods are revisited and support this conclusion once sufficiently high orders are included. The cause of the failure of MBPT methods for NIs of large systems is missing or incomplete “electrodynamic” screening of the Coulomb interaction due to induced particle–hole pairs between electrons in different monomers, leaving the effective interaction too strong for AC-SAPT to converge. Hence, MBPT cannot be considered reliable for quantitative predictions of NIs, even in moderately polarizable molecules with a few tens of atoms. The failure to accurately account for electrodynamic polarization makes MBPT qualitatively unsuitable for applications such as NIs of nanostructures, macromolecules, and soft materials; more robust non-perturbative approaches such as RPA or coupled cluster methods should be used instead whenever possible.<br>

Complete basis set limits of local second-order Møller–Plesset perturbation theory

2013 ◽  
Vol 111 (9-11) ◽  
pp. 1178-1189
Author(s):  
Kameron R. Jorgensen ◽  
Vinay V. Ramasesh ◽  
Sonja Hannibal ◽  
Nathan J. DeYonker ◽  
Angela K. Wilson
Keyword(s):  
Perturbation Theory ◽  
Second Order ◽  
Basis Set ◽  
Complete Basis ◽  
Complete Basis Set ◽  
Moller Plesset ◽  
Plesset Perturbation Theory

Divergence of Many-Body Perturbation Theory for Noncovalent Interactions of Large Molecules

2020 ◽  
Author(s):  
Brian Nguyen ◽  
Guo P Chen ◽  
Matthew M. Agee ◽  
Asbjörn M. Burow ◽  
Matthew Tang ◽  
...  
Keyword(s):  
Perturbation Theory ◽  
System Size ◽  
State Density ◽  
Second Order ◽  
Binding Energies ◽  
Basis Set ◽  
Many Body ◽  
Many Body Perturbation Theory ◽  
Relative Errors ◽  
Ground State Density

Prompted by recent reports of large errors in noncovalent interaction (NI) energies obtained from many-body perturbation theory (MBPT), we compare the performance of second-order Møller–Plesset MBPT (MP2), spin-scaled MP2, dispersion-corrected semilocal density functional approximations (DFA), and the post-Kohn–Sham random phase approximation (RPA) for predicting binding energies of supramolecular complexes contained in the S66, L7, and S30L benchmarks. All binding energies are extrapolated to the basis set limit, corrected for basis set superposition errors, and compared to reference results of the domain-based local pair-natural orbital coupled-cluster (DLPNO-CCSD(T)) or better quality. Our results confirm that MP2 severely overestimates binding energies of large complexes, producing relative errors of over 100% for several benchmark compounds. RPA relative errors consistently range between 5-10%, significantly less than reported previously using smaller basis sets, whereas spin-scaled MP2 methods show limitations similar to MP2, albeit less pronounced, and empirically dispersion-corrected DFAs perform almost as well as RPA. Regression analysis reveals a systematic increase of relative MP2 binding energy errors with the system size at a rate of approximately 0.1% per valence electron, whereas the RPA and dispersion-corrected DFA relative errors are virtually independent of the system size. These observations are corroborated by a comparison of computed rotational constants of organic molecules to gas-phase spectroscopy data contained in the ROT34 benchmark. To analyze these results, an asymptotic adiabatic connection symmetry-adapted perturbation theory (AC-SAPT) is developed which uses monomers at full coupling whose ground-state density is constrained to the ground-state density of the complex. Using the fluctuation–dissipation theorem, we obtain a nonperturbative “screened second-order” expression for the dispersion energy in terms of monomer quantities which is exact for non-overlapping subsystems and free of induction terms; a first-order RPA-like approximation to the Hartree, exchange, and correlation kernel recovers the macroscopic Lifshitz limit. The AC-SAPT expansion of the interaction energy is obtained from Taylor expansion of the coupling strength integrand. Explicit expressions for the convergence radius of the AC-SAPT series are derived within RPA and MBPT and numerically evaluated. Whereas the AC-SAPT expansion is always convergent for nondegenerate monomers when RPA is used, it is found to spuriously diverge for second-order MBPT, except for the smallest and least polarizable monomers. The divergence of the AC-SAPT series within MBPT is numerically confirmed within RPA; prior numerical results on the convergence of the SAPT expansion for MBPT methods are revisited and support this conclusion once sufficiently high orders are included. The cause of the failure of MBPT methods for NIs of large systems is missing or incomplete “electrodynamic” screening of the Coulomb interaction due to induced particle–hole pairs between electrons in different monomers, leaving the effective interaction too strong for AC-SAPT to converge. Hence, MBPT cannot be considered reliable for quantitative predictions of NIs, even in moderately polarizable molecules with a few tens of atoms. The failure to accurately account for electrodynamic polarization makes MBPT qualitatively unsuitable for applications such as NIs of nanostructures, macromolecules, and soft materials; more robust non-perturbative approaches such as RPA or coupled cluster methods should be used instead whenever possible.<br>

Divergence of Many-Body Perturbation Theory for Noncovalent Interactions of Large Molecules

2020 ◽  
Author(s):  
Brian Nguyen ◽  
Guo P Chen ◽  
Matthew M. Agee ◽  
Asbjörn M. Burow ◽  
Matthew Tang ◽  
...  
Keyword(s):  
Perturbation Theory ◽  
System Size ◽  
State Density ◽  
Second Order ◽  
Binding Energies ◽  
Basis Set ◽  
Many Body ◽  
Many Body Perturbation Theory ◽  
Relative Errors ◽  
Ground State Density

Prompted by recent reports of large errors in noncovalent interaction (NI) energies obtained from many-body perturbation theory (MBPT), we compare the performance of second-order Møller–Plesset MBPT (MP2), spin-scaled MP2, dispersion-corrected semilocal density functional approximations (DFA), and the post-Kohn–Sham random phase approximation (RPA) for predicting binding energies of supramolecular complexes contained in the S66, L7, and S30L benchmarks. All binding energies are extrapolated to the basis set limit, corrected for basis set superposition errors, and compared to reference results of the domain-based local pair-natural orbital coupled-cluster (DLPNO-CCSD(T)) or better quality. Our results confirm that MP2 severely overestimates binding energies of large complexes, producing relative errors of over 100% for several benchmark compounds. RPA relative errors consistently range between 5-10%, significantly less than reported previously using smaller basis sets, whereas spin-scaled MP2 methods show limitations similar to MP2, albeit less pronounced, and empirically dispersion-corrected DFAs perform almost as well as RPA. Regression analysis reveals a systematic increase of relative MP2 binding energy errors with the system size at a rate of approximately 0.1% per valence electron, whereas the RPA and dispersion-corrected DFA relative errors are virtually independent of the system size. These observations are corroborated by a comparison of computed rotational constants of organic molecules to gas-phase spectroscopy data contained in the ROT34 benchmark. To analyze these results, an asymptotic adiabatic connection symmetry-adapted perturbation theory (AC-SAPT) is developed which uses monomers at full coupling whose ground-state density is constrained to the ground-state density of the complex. Using the fluctuation–dissipation theorem, we obtain a nonperturbative “screened second-order” expression for the dispersion energy in terms of monomer quantities which is exact for non-overlapping subsystems and free of induction terms; a first-order RPA-like approximation to the Hartree, exchange, and correlation kernel recovers the macroscopic Lifshitz limit. The AC-SAPT expansion of the interaction energy is obtained from Taylor expansion of the coupling strength integrand. Explicit expressions for the convergence radius of the AC-SAPT series are derived within RPA and MBPT and numerically evaluated. Whereas the AC-SAPT expansion is always convergent for nondegenerate monomers when RPA is used, it is found to spuriously diverge for second-order MBPT, except for the smallest and least polarizable monomers. The divergence of the AC-SAPT series within MBPT is numerically confirmed within RPA; prior numerical results on the convergence of the SAPT expansion for MBPT methods are revisited and support this conclusion once sufficiently high orders are included. The cause of the failure of MBPT methods for NIs of large systems is missing or incomplete “electrodynamic” screening of the Coulomb interaction due to induced particle–hole pairs between electrons in different monomers, leaving the effective interaction too strong for AC-SAPT to converge. Hence, MBPT cannot be considered reliable for quantitative predictions of NIs, even in moderately polarizable molecules with a few tens of atoms. The failure to accurately account for electrodynamic polarization makes MBPT qualitatively unsuitable for applications such as NIs of nanostructures, macromolecules, and soft materials; more robust non-perturbative approaches such as RPA or coupled cluster methods should be used instead whenever possible.<br>

STACKING INTERACTION IN PYRAZINE DIMER

2006 ◽  
Vol 05 (03) ◽  
pp. 609-619 ◽  
Author(s):  
BRIJESH KUMAR MISHRA ◽  
N. SATHYAMURTHY
Keyword(s):  
Center Of Mass ◽  
Binding Energies ◽  
Mass Separation ◽  
Basis Set ◽  
Basis Sets ◽  
Mutual Orientation ◽  
Coupled Cluster ◽  
Basis Set Superposition Error ◽  
Moller Plesset ◽  
Plesset Perturbation Theory

The structure and stability of pyrazine dimer in different orientations have been investigated using second- and fourth-order Møller–Plesset perturbation theory (MP2, MP4) and coupled-cluster singles and doubles with non-iterative perturbative triples method [CCSD(T)] with various basis sets (6-31G*, 6-311G**, 6-311++G**, cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ). With the largest basis set (aug-cc-pVTZ) at the MP2 level of theory, the calculated binding energies (basis set superposition error corrected) for the sandwich, N–N axial-displaced, lateral-displaced, cross-displaced, T-shaped ( T N-ring and T H-ring ), and wedge-shaped ( W N–C and W C–C ) geometries are found to be 3.01, 4.18, 4.63, 5.92, 4.54, 2.78, 2.89, and 3.23 kcal/mol, respectively. Although MP2 calculations predict the cross-displaced geometry to be the most stable one, MP4 and CCSD(T) calculations show the most stable geometry to be T-shaped ( T N-ring ), with a center-of-mass separation of 4.2 Å. Dispersion interaction seems to be the major source of attraction in all the geometries considered, while the nature of electrostatic interaction depends on the mutual orientation of the two rings.

Sign in / Sign up

Export Citation Format

Share Document