Revealing the physical nature and the strength of charge‐inverted hydrogen bonds by SAPT(DFT), MP2, SCS‐MP2, MP2C,and CCSD(T) methods

The physical nature of charge‐inverted hydrogen bonds in H₃XH YH₃ (X = Si, Ge; Y = Al, Ga) dimer systems is studied by means of the SAPT(DFT)‐based decomposition of interaction energies and supermolecular interaction energies based on MP2, SCS‐MP2, MP2C, and CCSD(T) methods utilizing dimer‐centered aug‐cc‐pCVnZ (n = D, T, Q) basis sets as well as an extrapolation to the complete basis set limit. It is revealed that charge‐inverted hydrogen bonds are inductive in nature, although dispersion is also important. Computed interaction energies form the following relation: . It is confirmed that the aug‐cc‐pCVDZ basis set performs poorly and that very accurate values of interaction and dispersion energies require basis sets of at least quadrupole‐ζ quality. Considerably large binding energies suggest potential usefulness of charge‐inverted hydrogen bonds as an important structural motif in molecular binding. Terminology applying to σ‐ and π‐hole interactions as well as to triel and tetrel bonds is discussed. According to this new terminology the charge‐inverted hydrogen bond would become the first described case of a hydride‐triel bond. © 2017 Wiley Periodicals, Inc.     

DIESEL-MP2: a new program to perform large-scale multireference-MP2 computations

This article presents a new MR‐MP2 code (Multi‐Reference Møller–Plesset 2nd order) suitable for the computation MR‐MP2 energies of extended systems with strong near degeneracy effects (e.g., open shell systems). It is based on the DIESEL program package developed by Hanrath and Engels. Due to improved algorithms the new code is able to handle systems with 400–500 basis functions and more than 100 electrons. The code is made for parallel computers with distributed memory, but can also be run on local machines. It possesses two integral interfaces (MOLCAS, TURBOMOLE). The algorithms are briefly introduced and timings for the Neocarzinostatin chromophore are presented. The efficiencies of the codes obtained with Intel or GNU compilers are compared. © 2006 Wiley Periodicals, Inc. J Comput Chem 27: 1055–1062, 2006     

The GPU‐enabled divide‐expand‐consolidate RI‐MP2 method (DEC‐RI‐MP2)

We report porting of the Divide‐Expand‐Consolidate Resolution of the Identity second‐order Møller–Plesset perturbation (DEC‐RI‐MP2) method to the graphic processing units (GPUs) using OpenACC compiler directives. It is shown that the OpenACC compiler directives implementation efficiently accelerates the rate‐determining step of the DEC‐RI‐MP2 method with minor implementation effort. Moreover, the GPU acceleration results in a better load balance and thus in an overall scaling improvement of the DEC algorithm. The resulting cross‐platform hybrid MPI/OpenMP/OpenACC implementation has scalable and portable performance on heterogeneous HPC architectures. The GPU‐enabled code was benchmarked using a reduced version of the S12L test set of Stefan Grimme (Grimme, Chem. Eur. J. 2012, 18, 9955) consisting of supramolecular complexes up to 158 atoms and 4292 contracted basis functions (cc‐pVTZ). The test set results demonstrate the general applicability of the DEC‐RI‐MP2 method showing results consistent with the DEC‐RI‐MP2 introductory paper (Baudin et al., J. Chem. Phys. 2016, 144, 054102) on molecules of complicated electronic structures. © 2016 Wiley Periodicals, Inc.     

Fragment interaction analysis based on local MP2

We have developed a fragment interaction analysis based on local MP2 (FILM) in the context of the fragment molecular orbital (FMO) scheme. The primary purpose of this work is to provide a tool for analyzing inter-fragment interaction associated with dispersion interactions in a large molecule such as protein and DNA. Our implementation of local MP2 (LMP2) is based on the algorithm developed by Pulay and Werner. A potential of FILM was demonstrated using the human immunodeficiency virus type 1 protease (HIV-1 PR) complexed with lopinavir (LPV). The total energy, binding affinity, and inter-fragment interaction energy (IFIE) by the FMO method using LMP2 were compared with those obtained by canonical MP2 and the site-specific information in dispersion interaction was obtained. It turned out that the FILM is a useful tool for analyzing the dispersion interaction between an amino acid residue and a specific site of a ligand.     

Analytic energy gradients for orbital‐optimized MP3 and MP2.5 with the density‐fitting approximation: An efficient implementation

Efficient implementations of analytic gradients for the orbital‐optimized MP3 and MP2.5 and their standard versions with the density‐fitting approximation, which are denoted as DF‐MP3, DF‐MP2.5, DF‐OMP3, and DF‐OMP2.5, are presented. The DF‐MP3, DF‐MP2.5, DF‐OMP3, and DF‐OMP2.5 methods are applied to a set of alkanes and noncovalent interaction complexes to compare the computational cost with the conventional MP3, MP2.5, OMP3, and OMP2.5. Our results demonstrate that density‐fitted perturbation theory (DF‐MP) methods considered substantially reduce the computational cost compared to conventional MP methods. The efficiency of our DF‐MP methods arise from the reduced input/output (I/O) time and the acceleration of gradient related terms, such as computations of particle density and generalized Fock matrices (PDMs and GFM), solution of the Z‐vector equation, back‐transformations of PDMs and GFM, and evaluation of analytic gradients in the atomic orbital basis. Further, application results show that errors introduced by the DF approach are negligible. Mean absolute errors for bond lengths of a molecular set, with the cc‐pCVQZ basis set, is 0.0001–0.0002 Å. © 2017 Wiley Periodicals, Inc.     

An improved algorithm for analytical gradient evaluation in resolution-of-the-identity second-order Møller-Plesset perturbation theory: application to alanine tetrapeptide conformational analysis

We present a new algorithm for analytical gradient evaluation in resolution‐of‐the‐identity second‐order Møller‐Plesset perturbation theory (RI‐MP2) and thoroughly assess its computational performance and chemical accuracy. This algorithm addresses the potential I/O bottlenecks associated with disk‐based storage and access of the RI‐MP2 t‐amplitudes by utilizing a semi‐direct batching approach and yields computational speed‐ups of approximately 2–3 over the best conventional MP2 analytical gradient algorithms. In addition, we attempt to provide a straightforward guide to performing reliable and cost‐efficient geometry optimizations at the RI‐MP2 level of theory. By computing relative atomization energies for the G3/99 set and optimizing a test set of 136 equilibrium molecular structures, we demonstrate that satisfactory relative accuracy and significant computational savings can be obtained using Pople‐style atomic orbital basis sets with the existing auxiliary basis expansions for RI‐MP2 computations. We also show that RI‐MP2 geometry optimizations reproduce molecular equilibrium structures with no significant deviations (>0.1 pm) from the predictions of conventional MP2 theory. As a chemical application, we computed the extended‐globular conformational energy gap in alanine tetrapeptide at the extrapolated RI‐MP2/cc‐pV(TQ)Z level as 2.884, 4.414, and 4.994 kcal/mol for structures optimized using the HF, DFT (B3LYP), and RI‐MP2 methodologies and the cc‐pVTZ basis set, respectively. These marked energetic discrepancies originate from differential intramolecular hydrogen bonding present in the globular conformation optimized at these levels of theory and clearly demonstrate the importance of long‐range correlation effects in polypeptide conformational analysis. © 2007 Wiley Periodicals, Inc. J Comput Chem, 2007     

Intermolecular interactions of formic acid with benzene: Energy decomposition analyses with ab initio MP2 and double‐hybrid density functional computations

The intermolecular interactions of formic acid (HCOOH) with benzene (C₆H₆) have been investigated using localized molecular orbital energy decomposition analyses (LMO‐EDA) with ab initio MP2 and several double‐hybrid density functionals. The molecular geometries of five HCOOH…C₆H₆ complexes and corresponding benchmark total interaction energies at the CCSD(T)/CBS level are taken from literature (Zhao et al., J. Chem. Theory Comput. 2009, 5, 2726). According to the results of LMO‐EDA with the MP2 method, the dispersion energies are found to be as important as the electrostatic energies for the total interaction energies of the five HCOOH…C₆H₆ complexes. Based on LMO‐EDA with the double‐hybrid density functionals of B2PLYP, B2K‐PLYP, B2T‐PLYP, and B2GP‐PLYP computations, two new parameters for the framework of B2PLYP are extrapolated. These two new parameters are tested with other 10 complexes involving C₆H₆ (Crittenden, J. Phys. Chem. A 2009, 113, 1663), and they perform well on predicting the corresponding total interaction energies. Interestingly, these two new parameters for the framework of B2PLYP also perform well on the noncovalent complexation energies database (NCCE31/05) developed by Truhlar's group (Zhao and Truhlar, J. Phys. Chem. A 2005, 109, 5656). Therefore, these two new parameters appear to be suitable for investigating the noncovalent interactions, and they are denoted as B2N‐PLYP, where N stands for the noncovalent interaction. This study is expected to provide new insight into the derivation of double‐hybrid density functionals for studying the noncovalent interactions. © 2013 Wiley Periodicals, Inc.     

The MC‐DFT approach including the SCS‐MP2 energies to the new minnesota‐type functionals

We have applied the multicoefficient density functional theory (MC‐DFT) to four recent Minnesota functionals, including M06‐2X, M08‐HX, M11, and MN12‐SX on the performance of thermochemical kinetics. The results indicated that the accuracy can be improved significantly using more than one basis set. We further included the SCS‐MP2 energies into MC‐DFT, and the resulting mean unsigned errors (MUEs) decreased by approximately 0.3 kcal/mol for the most accurate basis set combinations. The M06‐2X functional with the simple [6–311+G(d,p)/6–311+G(2d,2p)] combination gave the best performance/cost ratios for the MC‐DFT and MC‐SCS‐MP2|MC‐DFT methods with MUE of 1.58 and 1.22 kcal/mol, respectively. © 2014 Wiley Periodicals, Inc.