Toward true DNA base-stacking energies: MP2, CCSD(T), and complete basis set calculations

Stacking energies in low-energy geometries of pyrimidine, uracil, cytosine, and guanine homodimers were determined by the MP2 and CCSD(T) calculations utilizing a wide range of split-valence, correlation-consistent, and bond-functions basis sets. Complete basis set MP2 (CBS MP2) stacking energies extrapolated using aug-cc-pVXZ (X = D, T, and for pyrimidine dimer Q) basis sets equal to -5.3, -12.3, and -11.2 kcal/mol for the first three dimers, respectively. Higher-order correlation corrections estimated as the difference between MP2 and CCSD(T) stacking energies amount to 2.0, 0.7, and 0.9 kcal/mol and lead to final estimates of the genuine stacking energies for the three dimers of -3.4, -11.6, and -10.4 kcal/mol. The CBS MP2 stacking-energy estimate for guanine dimer (-14.8 kcal/mol) was based on the 6-31G(0.25) and aug-cc-pVDZ calculations. This simplified extrapolation can be routinely used with a meaningful accuracy around 1 kcal/mol for large aromatic stacking clusters. The final estimate of the guanine stacking energy after the CCSD(T) correction amounts to -12.9 kcal/mol. The MP2/6-31G(0.25) method previously used as the standard level to calculate aromatic stacking in hundreds of geometries of nucleobase dimers systematically underestimates the base stacking by ca. 1.0-2.5 kcal/mol per stacked dimer, covering 75-90% of the intermolecular correlation stabilization. We suggest that this correction is to be considered in calibration of force fields and other cheaper computational methods. The quality of the MP2/6-31G(0.25) predictions is nevertheless considerably better than suggested on the basis of monomer polarizability calculations. Fast and very accurate estimates of the MP2 aromatic stacking energies can be achieved using the RI-MP2 method. The CBS MP2 calculations and the CCSD(T) correction, when taken together, bring only marginal changes to the relative stability of H-bonded and stacked base pairs, with a slight shift of ca. 1 kcal/mol in favor of H-bonding. We suggest that the present values are very close to ultimate predictions of the strength of aromatic base stacking of DNA and RNA bases.     

Hydrogen-bonded nucleic acid base pairs containing unusual base tautomers: complete basis set calculations at the MP2 and CCSD(T) levels

The total interaction energies of altogether 15 hydrogen-bonded nucleic acid base pairs containing unusual base tautomers were calculated. The geometry properties of all selected adenine-thymine and guanine-cytosine hydrogen-bonded base pairs enable their incorporation into DNA. Unusual base pairing patterns were compared with Watson-Crick H-bonded structures of the adenine-thymine and guanine-cytosine pairs. The complete basis set (CBS) limit of the MP2 interaction energy and the CCSD(T) correction term, determined as the difference between the CCSD(T) and MP2 interaction energies, was evaluated. Extrapolation to the MP2 CBS limit was done using the aug-cc-pVDZ and aug-cc-pVTZ results, and the CCSD(T) correction term was determined with the 6-31G*(0.25) basis set. Final interaction energies were corrected while taking into account both tautomeric penalization determined at the CBS level and solvation/desolvation free energies. The situation for the adenine-thymine pairs is straightforward, and tautomeric pairs are significantly less stable than the Watson-Crick pair consisting of the canonical forms. In the case of the guanine-cytosine pair, the Watson-Crick structure made by canonical forms is again the most stable. The other two structures are, however, energetically rather similar (by 5 and 6 kcal/mol), which provides a very small but non-negligible chance of detecting these structures in the DNA double helix (1:5000). Due to the fact that DNA bases and base pairs incorporated into DNA are solvated less favorably than in isolated systems, this probability represents the very upper limit. The results clearly show how precisely the canonical building blocks of DNA molecules were chosen and how well their stability is maintained.     

Development and testing of a general amber force field

We describe here a general Amber force field (GAFF) for organic molecules. GAFF is designed to be compatible with existing Amber force fields for proteins and nucleic acids, and has parameters for most organic and pharmaceutical molecules that are composed of H, C, N, O, S, P, and halogens. It uses a simple functional form and a limited number of atom types, but incorporates both empirical and heuristic models to estimate force constants and partial atomic charges. The performance of GAFF in test cases is encouraging. In test I, 74 crystallographic structures were compared to GAFF minimized structures, with a root-mean-square displacement of 0.26 A, which is comparable to that of the Tripos 5.2 force field (0.25 A) and better than those of MMFF 94 and CHARMm (0.47 and 0.44 A, respectively). In test II, gas phase minimizations were performed on 22 nucleic acid base pairs, and the minimized structures and intermolecular energies were compared to MP2/6-31G* results. The RMS of displacements and relative energies were 0.25 A and 1.2 kcal/mol, respectively. These data are comparable to results from Parm99/RESP (0.16 A and 1.18 kcal/mol, respectively), which were parameterized to these base pairs. Test III looked at the relative energies of 71 conformational pairs that were used in development of the Parm99 force field. The RMS error in relative energies (compared to experiment) is about 0.5 kcal/mol. GAFF can be applied to wide range of molecules in an automatic fashion, making it suitable for rational drug design and database searching.     

Semi-empirical molecular orbital methods including dispersion corrections for the accurate prediction of the full range of intermolecular interactions in biomolecules

Semi-empirical calculations including an empirical dispersive correction are used to calculate intermolecular interaction energies and structures for a large database containing 156 biologically relevant molecules (hydrogen-bonded DNA base pairs, interstrand base pairs, stacked base pairs and amino acid base pairs) for which MP2 and CCSD(T) complete basis set (CBS) limit estimates of the interaction energies are available. The dispersion corrected semi-empirical methods are parameterised against a small training set of 22 complexes having a range of biologically important non-covalent interactions. For the full molecule set (156 complexes), compared to the high-level ab initio database, the mean unsigned errors of the interaction energies at the corrected semi-empirical level are 1.1 (AM1-D) and 1.2 (PM3-D) kcal mol(-1), being a significant improvement over existing AM1 and PM3 methods (8.6 and 8.2 kcal mol(-1)). Importantly, the new semi-empirical methods are capable of describing the diverse range of biological interactions, most notably stacking interactions, which are poorly described by both current AM1 and PM3 methods and by many DFT functionals. The new methods require no more computer time than existing semi-empirical methods and therefore represent an important advance in the study of important biological interactions.     

Probing phenylalanine/adenine pi-stacking interactions in protein complexes with explicitly correlated and CCSD(T) computations

To examine the effects of pi-stacking interactions between aromatic amino acid side chains and adenine bearing ligands in crystalline protein structures, 26 toluene/(N9-methyl)adenine model configurations have been constructed from protein/ligand crystal structures. Full geometry optimizations with the MP2 method cause the 26 crystal structures to collapse to six unique structures. The complete basis set (CBS) limit of the CCSD(T) interaction energies has been determined for all 32 structures by combining explicitly correlated MP2-R12 computations with a correction for higher-order correlation effects from CCSD(T) calculations. The CCSD(T) CBS limit interaction energies of the 26 crystal structures range from -3.19 to -6.77 kcal mol (-1) and average -5.01 kcal mol (-1). The CCSD(T) CBS limit interaction energies of the optimized complexes increase by roughly 1.5 kcal mol (-1) on average to -6.54 kcal mol (-1) (ranging from -5.93 to -7.05 kcal mol (-1)). Corrections for higher-order correlation effects are extremely important for both sets of structures and are responsible for the modest increase in the interaction energy after optimization. The MP2 method overbinds the crystal structures by 2.31 kcal mol (-1) on average compared to 4.50 kcal mol (-1) for the optimized structures.     

Polarizable model potential function for nucleic acid bases

A polarizable model potential (PMP) function for adenine (A), cytosine (C), guanine (G), thymine (T), and uracil (U) is developed on the basis of ab initio molecular orbital calculations at the MP2/6-31+G* level. The PMP function consists of Coulomb, van der Waals, and polarization terms. The permanent atomic charges of the Coulomb term are determined by using electrostatic potential (ESP) optimization. The multicenter polarizabilities of the polarization term are determined by using polarized one-electron potential (POP) optimization in which the electron density changes induced by a test charge are target. Isotropic and anisotropic polarizabilities are adopted as the multicenter polarizabilities. In the PMP calculations using the optimized parameters, the interaction energies of Watson-Crick type A-T and C-G base pairs were -15.6 and -29.4 kcal/mol, respectively. The interaction energy of Hoogsteen type A-T base pair was -17.8 kcal/mol. These results reproduce well the quantum chemistry calculations at the MP2/6-311++G(3df,2pd) level within the differences of 0.6 kcal/mol. The stacking energies of A-T and C-G were -9.7 and -10.9 kcal/mol. These reproduce well the calculation results at the MP2/6-311++G (2d,2p) level within the differences of 1.3 kcal/mol. The potential energy surfaces of the system in which a sodium ion or a chloride ion is adjacent to the nucleic acid base are calculated. The interaction energies of the PMP function reproduced well the calculation results at the MP2/6-31+G* or MP2/6-311++G(2d,2p) level. The reason why the PMP function reproduces well the high-level quantum mechanical interaction energies is addressed from the viewpoint of each energy terms.     

Quantum dynamics of vibrationally activated OH-CO reactant complexes

The MP2 complete basis set (CBS) limit for the binding energy of the two low-lying water octamer isomers of D2d and S4 symmetry is estimated at -72.7+/-0.4 kcal/mol using the family of augmented correlation-consistent orbital basis sets of double through quintuple zeta quality. The largest MP2 calculation with the augmented quintuple zeta (aug-cc-pV5Z) basis set produced binding energies of -73.70 (D2d) and -73.67 kcal/mol (S4). The effects of higher correlation, computed at the CCSD(T) level of theory, are estimated at <0.1 kcal/mol. The newly established MP2/CBS limit for the water octamer is reproduced quite accurately by the newly developed all atom polarizable, flexible interaction potential (TTM2-F). The TTM2-F binding energies of -73.21 (D2d) and -73.24 kcal/mol (S4) for the two isomers are just 0.5 kcal/mol (or 0.7%) larger than the MP2/CBS limit.     

Stabilization energies of the hydrogen-bonded and stacked structures of nucleic acid base pairs in the crystal geometries of CG, AT, and AC DNA steps and in the NMR geometry of the 5'-d(GCGAAGC)-3' hairpin: Complete basis set calculations at the MP2 and CCSD(T) levels

Stabilization energies of the H-bonded and stacked structures of a DNA base pair were studied in the crystal structures of adenine-thymine, cytosine-guanine, and adenine-cytosine steps as well as in the 5'-d(GCGAAGC)-3' hairpin (utilizing the NMR geometry). Stabilization energies were determined as the sum of the complete basis set (CBS) limit of MP2 stabilization energies and the Delta E(CCSD(T)) - Delta E(MP2) correction term evaluated with the 6-31G*(0.25) basis set. The CBS limit was determined by a two-point extrapolation using the aug-cc-pVXZ basis sets for X = D and T. While the H-bonding energies are comparable to those of base pairs in a crystal and a vacuum, the stacking energies are considerably smaller in a crystal. Despite this, the stacking is still important and accounts for a significant part of the overall stabilization. It contributes equally to the stability of DNA as does H-bonding for AT-rich DNAs, while in the case of GC-rich DNAs it forms about one-third of the total stabilization. Interstrand stacking reaches surprisingly large values, well comparable to the intrastrand ones, and thus contributes significantly to the overall stabilization. The hairpin structure is characterized by significant stacking, and both guanine...cytosine pairs possess stacking energies larger than 11.5 kcal/mol. A high portion of stabilization in the studied hairpin comes from stacking (similar to that found for AT-rich DNAs) despite the fact that it contains two GC Watson-Crick pairs having very large H-bonding stabilization. The DFT/B3LYP/6-31G** method yields satisfactory values of interaction energies for H-bonded structures, while it fails completely for stacking.     

Probing the effects of heterogeneity on delocalized pi...pi interaction energies

Dimers composed of benzene (Bz), 1,3,5-triazine (Tz), cyanogen (Cy) and diacetylene (Di) are used to examine the effects of heterogeneity at the molecular level and at the cluster level on pi...pi stacking energies. The MP2 complete basis set (CBS) limits for the interaction energies (E(int)) of these model systems were determined with extrapolation techniques designed for correlation consistent basis sets. CCSD(T) calculations were used to correct for higher-order correlation effects (deltaE(CCSD)(T)(MP2)) which were as large as +2.81 kcal mol(-1). The introduction of nitrogen atoms into the parallel-slipped dimers of the aforementioned molecules causes significant changes to E(int). The CCSD(T)/CBS E(int) for Di-Cy is -2.47 kcal mol(-1) which is substantially larger than either Cy-Cy (-1.69 kcal mol(-1)) or Di-Di (-1.42 kcal mol(-1)). Similarly, the heteroaromatic Bz-Tz dimer has an E(int) of -3.75 kcal mol(-1) which is much larger than either Tz-Tz (-3.03 kcal mol(-1)) or Bz-Bz (-2.78 kcal mol(-1)). Symmetry-adapted perturbation theory calculations reveal a correlation between the electrostatic component of E(int) and the large increase in the interaction energy for the mixed dimers. However, all components (exchange, induction, dispersion) must be considered to rationalize the observed trend. Another significant conclusion of this work is that basis-set superposition error has a negligible impact on the popular deltaE(CCSD)(T)(MP2) correction, which indicates that counterpoise corrections are not necessary when computing higher-order correlation effects on E(int). Spin-component-scaled MP2 (SCS-MP2 and SCSN-MP2) calculations with a correlation-consistent triple-zeta basis set reproduce the trends in the interaction energies despite overestimating the CCSD(T)/CBS E(int) of Bz-Tz by 20-30%.     

Physical nature of silane⋯carbene dimers revealed by state-of-the-art ab initio calculations

Using the SAPT2 + 3(CCD)δMP2 method in complete basis set (CBS) limit, it is shown that the interactions in the recently studied silane⋯carbene dimers are mainly dispersive in nature. Consequently, slow convergence of dispersion energy also forces slow convergence of the interaction energy. Therefore, obtaining very accurate values requires extrapolation of the correlation part to the CBS limit. The most accurate values obtained at the CCSD(T)/CBS level of theory show that the studied silane⋯carbene dimers are rather weakly bound, with interaction energies ranging from about −1.9 to −1.3 kcal/mol. Comparing to CCSD(T)/CBS, it will be shown that SCS-MP2 and MP2C methods clearly underestimate and methods based on SAPT2+ and having some third-order corrections, as well as the MP2 method, overestimate values of interaction energies. Popular SAPT(DFT) method performs better than SCS-MP2 and MP2C; nevertheless, underestimation is still considerable. The underestimation is slightly quenched if third-order dispersion energy and its exchange counterpart is added to the SAPT(DFT). The closest value of CCSD(T)/CBS has been given by the SAPT2 + (3)(CCD)δMP2 method in quadruple-ζ basis set. © 2019 Wiley Periodicals, Inc.     

Ab initio study of hydrogen-bond formation between aliphatic and phenolic hydroxy groups and selected amino acid side chains

Hydrogen bonding was studied in 24 pairs of isopropyl alcohol and phenol as one partner, and water and amino-acid mimics (methanol, acetamide, neutral and protonated imidazole, protonated methylalamine, methyl-guanidium cation, and acetate anion) as the other partner. MP2/6-31+G* and MP2/aug-cc-pvtz calculations were conducted in the gas phase and in a model continuum dielectric environment with dielectric constant of 15.0. Structures were optimized in the gas phase with both basis sets, and zero-point energies were calculated at the MP2/6-31+G* level. At the MP2/aug-cc-pvtz level, the BSSE values from the Boys-Bernardi counterpoise calculations amount to 10-20 and 5-10% of the uncorrected binding energies of the neutral and ionic complexes, respectively. The geometry distortion energy upon hydrogen-bond formation is up to 2 kcal/mol, with the exception of the most strongly bound complexes. The BSSE-corrected MP2/aug-cc-pvtz binding energy of -27.56 kcal/mol for the gas-phase acetate...phenol system has been classified as a short and strong hydrogen bond (SSHB). The CH3NH3+...isopropyl alcohol complex with binding energy of -22.54 kcal/mol approaches this classification. The complete basis set limit (CBS) for the binding energy was calculated for twelve and six complexes on the basis of standard and counterpoise-corrected geometry optimizations, respectively. The X...Y distances of the X-H...Y bridges differ by up to 0.03 A as calculated by the two methods, whereas the corresponding CBS energy values differ by up to 0.03 kcal/mol. Uncorrected MP2/aug-cc-pvtz hydrogen-bonding energies are more negative by up to 0.35 kcal/mol than the MP2/CBS values, and overestimate the CCSD(T)/CBS binding energies generally by up to 5% for the eight studied complexes in the gas phase. The uncorrected MP2/aug-cc-pvtz binding energies decreased (in absolute value) by 11-18 kcal/mol for the ionic species and by up to 5 kcal/mol for the neutral complexes when the electrostatic effect of a polarizable model environment was considered. The DeltaECCSD(T) - DeltaEMP2 corrections still remained close to their gas-phase values for four complexes with 0, +/-1 net charges. Good correlations (R2 = 0.918-0.958) for the in-environment MP2/aug-cc-pvtz and MP2/6-31+G* hydrogen-bonding energies facilitate the high-level prediction of these energies on the basis of relatively simple MP2/6-31+G* calculations.     

Intercalators. 1. Nature of stacking interactions between intercalators (ethidium, daunomycin, ellipticine, and 4',6-diaminide-2-phenylindole) and DNA base pairs. Ab initio quantum chemical, density functional theory, and empirical potential study

Properties of isolated intercalators (ethidium (E), daunomycin (D), ellipticine (EL), and 4,6'-diaminide-2-phenylindole (DAPI)) and their stacking interactions with adenine...thymine (AT) and guanine...cytosine (GC) nucleic acid base pairs were investigated by means of a nonempirical correlated ab initio method. All intercalators exhibit large charge delocalization, and none of them (including the DAPI dication) exhibits a site with dominant charge. All intercalators have large polarizability and are good electron acceptors, while base pairs are good electron donors. MP2/6-31G*(0.25) stabilization energies of intercalator...base pair complexes are large (E...AT, 22.4 kcal/mol; D...GC, 17.8 kcal/mol; EL...GC, 18.2 kcal/mol; DAPI...GC, 21.1 kcal/mol) and are well reproduced by modified AMBER potential (van der Waals radii of intercalator atoms are enlarged and their energy depths are increased). Standard AMBER potential underestimates binding, especially for DAPI-containing complexes. Because the DAPI dication is the best electron acceptor (among all intercalators studied), this difference is explained by the importance of the charge-transfer term, which is not included in the AMBER potential. For the neutral EL molecule, the standard AMBER force field provides correct results. The Hartree-Fock and DFT/B3LYP methods, not covering the dispersion energy, fail completely to reveal any energy minimum at the potential energy curve of the E...AT complex, and these methods thus cannot be recommended for a study of intercalation process. On the other hand, an approximate version of the DFT method, which was extended to cover London dispersion energy, yields for all complexes very good stabilization energies that are well comparable with referenced ab initio data. Besides the vertical dependence of the interaction, an energy twist dependence of the interaction energy was also investigated by a reference correlated ab initio method and empirical potentials. It is concluded that, despite the cationic (E +1, D +1, DAPI +2) or polar (EL) character of the intercalators investigated, it is the dispersion energy which predominantly contributes to the stability of intercalator...base pair complexes. Any procedure which does not cover dispersion energy is thus not suitable for studying the process of intercalation.     

Reliable predictions of the thermochemistry of boron-nitrogen hydrogen storage compounds: BxNxHy, x = 2, 3

Thermochemical data calculated using ab initio molecular orbital theory are reported for 16 BxNxHy compounds with x = 2, 3 and y > or = 2x. Accurate gas-phase heats of formation were obtained using coupled cluster with single and double excitations and perturbative triples (CCSD(T)) valence electron calculations extrapolated to the complete basis set (CBS) limit with additional corrections including core/valence, scalar relativistic, and spin-orbit corrections to predict the atomization energies and scaled harmonic frequencies to correct for zero point and thermal energies and estimate entropies. Computationally cheaper calculations were also performed using the G3MP2 and G3B3 variants of the Gaussian 03 method, as well as density functional theory (DFT) using the B3LYP functional. The G3MP2 heats of formation are too positive by up to approximately 6 kcal/mol as compared with CCSD(T)/CBS values. The more expensive G3B3 method predicts heats of formation that are too negative as compared with the CCSD(T)/CBS values by up to 3-4 kcal/mol. DFT using the B3LYP functional and 6-311+G** basis set predict isodesmic reaction energies to within a few kcal/mol compared with the CCSD(T)/CBS method so isodesmic reactions involving BN compounds and the analogous hydrocarbons can be used to estimate heats of formation. Heats of formation of c-B3N3H12 and c-B3N3H6 are -95.5 and -115.5 kcal/mol at 298 K, respectively, using our best calculated CCSD(T)/CBS approach. The experimental value for c-B3N3H6 appears to be approximately 7 kcal/mol too negative. Enthalpies, entropies, and free energies are calculated for many dehydrocoupling and dehydrogenation reactions that convert BNH6 to alicyclic and cyclic oligomers and H2(g). Generally, the reactions are highly exothermic and exergonic as well because of the release of 1 or more equivalents of H2(g). For c-B3N3H12 and c-B3N3H6, available experimental data for sublimation and vaporization lead to estimates of their condensed phase 298 K heats of formation: DeltaHf degrees [c-B3N3H12(s)] = -124 kcal/mol and DeltaHf degrees [c-B3N3H6(l)] = -123 kcal/mol. The reaction thermochemistries for the dehydrocoupling of BNH6(s) to c-B3N3H12(s) and the dehydrogenation of c-B3N3H12(s) to c-B3N3H6(l) are much less exothermic compared with the gas-phase reactions due to intermolecular forces which decrease in the order BNH6 > cyclo-B3N3H12 > cyclo-B3N3H6. The condensed phase reaction free energies are less negative compared with the gas-phase reactions but are still too favorable for BNH6 to be regenerated from either c-B3N3H12 or c-B3N3H6 by just an overpressure of H2.     

Performance of spin-component-scaled Møller-Plesset theory (SCS-MP2) for potential energy curves of noncovalent interactions

An examination of the performance of density-fitted, spin-component-scaled, second-order M?ller-Plesset theory (SCS-MP2), SCS-MP2 with parameters optimized for nucleic acids (SCSN-MP2), and their local-correlation variants, SCS-LMP2 and SCSN-LMP2, is presented for the sandwich and T-shaped benzene dimers, the methane-benzene and H(2)S-benzene complexes, and the methane dimer over entire potential energy curves. These are compared to benchmark-quality estimates of the complete-basis-set limit for coupled-cluster theory through perturbative triple excitations, CCSD(T)/CBS. With the exception of the methane dimer, SCSN-LMP2/CBS tends to outperform SCS-LMP2/CBS with maximum relative errors of 6 and 18%, respectively, at the optimal CCSD(T)/CBS intermolecular distances. For the methane dimer, errors for SCS(N)-(L)MP2/CBS remain in the 0.2-0.3 kcal mol(-1) range, corresponding to a larger relative error of 40-50%. Although the local MP2 methods perform very similarly to their conventional counterparts when aug-cc-pVTZ or larger basis sets are used, in the absence of counterpoise correction the local approximation becomes significantly worse for the aug-cc-pVDZ basis set. The changes due to local correlation approximations for the aug-cc-pVDZ basis are reduced when diffuse functions are neglected for hydrogen atoms.     

Stabilisation energy of C(6)H(6)...C(6)X(6) (X = F, Cl, Br, I, CN) complexes: complete basis set limit calculations at MP2 and CCSD(T) levels

Stabilisation energies of stacked structures of C(6)H(6)...C(6)X(6) (X = F, Cl, Br, CN) complexes were determined at the CCSD(T) complete basis set (CBS) limit level. These energies were constructed from MP2/CBS stabilisation energies and a CCSD(T) correction term determined with a medium basis set (6-31G**). The former energies were extrapolated using the two-point formula of Helgaker et al. from aug-cc-pVDZ and aug-cc-pVTZ Hartree-Fock energies and MP2 correlation energies. The CCSD(T) correction term is systematically repulsive. The final CCSD(T)/CBS stabilisation energies are large, considerably larger than previously calculated and increase in the series as follows: hexafluorobenzene (6.3 kcal mol(-1)), hexachlorobenzene (8.8 kcal mol(-1)), hexabromobenzene (8.1 kcal mol(-1)) and hexacyanobenzene (11.0 kcal mol(-1)). MP2/SDD** relativistic calculations performed for all complexes mentioned and also for benzene[dot dot dot]hexaiodobenzene have clearly shown that due to relativistic effects the stabilisation energy of the hexaiodobenzene complex is lower than that of hexabromobenzene complex. The decomposition of the total interaction energy to physically defined energy components was made by using the symmetry adapted perturbation treatment (SAPT). The main stabilisation contribution for all complexes investigated is due to London dispersion energy, with the induction term being smaller. Electrostatic and induction terms which are attractive are compensated by their exchange counterparts. The stacked motif in the complexes studied is very stable and might thus be valuable as a supramolecular synthon.     

An ab initio benchmark study of hydrogen bonded formamide dimers

The five singly and doubly hydrogen bonded dimers of formamide are calculated at the correlated level by using resolution of identity M?ller-Plesset second-order perturbation theory (RIMP2) and the coupled cluster with singles, doubles, and perturbative triples [CCSD(T)] method. All structures are optimized with the Dunning aug-cc-pVTZ and aug-cc-pVQZ basis sets. The binding energies are extrapolated to the complete basis set (CBS) limit by using the aug-cc-pVXZ (X = D, T, Q) basis set series. The effect of extending the basis set to aug-cc-pV5Z on the geometries and binding energies is studied for the centrosymmetric doubly N-H...O bonded dimer FA1 and the doubly C-H...O bonded dimer FA5. The MP2 CBS limits range from -5.19 kcal/mol for FA5 to -14.80 kcal/mol for the FA1 dimer. The DeltaCCSD(T) corrections to the MP2 CBS limit binding energies calculated with the 6-31+G(d,p), aug-cc-pVDZ, and aug-cc-pVTZ basis sets are mutually consistent to within < or =0.03 kcal/mol. The DeltaCCSD(T) correction increases the binding energy of the C-H...O bonded FA5 dimer by 0.4 kcal/mol or approximately 9% over the distance range +/-0.5 Angstrom relative to the potential minimum. This implies that the ubiquitous long-range C-H...O interactions in proteins are stronger than hitherto calculated.     

Structure and binding energy of anion-pi and cation-pi complexes: a comparison of MP2, RI-MP2, DFT, and DF-DFT methods

Several complexes of benzene with cations, hexafluorobenzene with anions, 1,3,5-trifluorobenzene with cations and anions, and s-triazine with cations and anions have been evaluated and compared at the MP2 and resolution of the identity MP2 (RI-MP2) levels. The RI-MP2 method is considerably faster than the MP2 and the interaction energies and equilibrium distances are almost identical for both methods. A similar result is found when comparing DFT and density fitting DFT (DF-DFT) levels. Therefore RI-MP2 and DF-DFT methods are well suited for the study of ion-pi interactions.     

Toward a less costly but accurate calculation of the CCSD(T)/CBS noncovalent interaction energy

The popular method of calculating the noncovalent interaction energies at the coupled-cluster single-, double-, and perturbative triple-excitations [CCSD(T)] theory level in the complete basis set (CBS) limit was to add a CCSD(T) correction term to the CBS second-order Møller-Plesset perturbation theory (MP2). The CCSD(T) correction term is the difference between the CCSD(T) and MP2 interaction energies evaluated in a medium basis set. However, the CCSD(T) calculations with the medium basis sets are still very expensive for systems with more than 30 atoms. Comparatively, the domain-based local pair natural orbital coupled-cluster method [DLPNO-CCSD(T)] can be applied to large systems with over 1,000 atoms. Considering both the computational accuracy and efficiency, in this work, we propose a new scheme to calculate the CCSD(T)/CBS interaction energies. In this scheme, the MP2/CBS term keeps intact and the CCSD(T) correction term is replaced by a DLPNO-CCSD(T) correction term which is the difference between the DLPNO-CCSD(T) and DLPNO-MP2 interaction energies evaluated in a medium basis set. The interaction energies of the noncovalent systems in the S22, HSG, HBC6, NBC10, and S66 databases were recalculated employing this new scheme. The consistent and tight settings of the truncation parameters for DLPNO-CCSD(T) and DLPNO-MP2 in this noncanonical CCSD(T)/CBS calculations lead to the maximum absolute deviation and root-mean-square deviation from the canonical CCSD(T)/CBS interaction energies of less than or equal to 0.28 kcal/mol and 0.09 kcal/mol, respectively. The high accuracy and low cost of this new computational scheme make it an excellent candidate for the study of large noncovalent systems.     

Ab initio prediction of the gas- and solution-phase acidities of strong Brønsted acids: the calculation of pKa values less than -10

The intrinsic gas-phase acidities of a series of 21 Br?nsted acids have been predicted with G3(MP2) theory. The G3(MP2) results agree with high level CCSD(T)/CBS acidities for H(2)SO(4), FSO(3)H, CH(3)SO(3)H, and CF(3)SO(3)H to within 1 kcal/mol. The G3(MP2) results are in excellent agreement with experimental gas-phase acidities in the range 342-302 kcal/mol to within <1 kcal/mol for 14 out of 15 acids. Five of the six acids in the range of 302-289 kcal/mol had an average deviation of 5.5 kcal/mol and the strongest acid, (CF(3)SO(2))(3)CH, deviated by 15.0 kcal/mol. These high-level calculations strongly suggest that the experimental acidities in this very acidic part of the scale need to be remeasured. The CCSD(T)/CBS (mixed exponential Gaussian) additive approach for CH(3)CO(2)H, HNO(3), H(2)SO(4), CH(3)SO(3)H, FSO(3)H, and CF(3)SO(3)H gives excellent agreement (+/-1 kcal/mol) with experiment for the DeltaH(f)(0)'s of non-sulfur containing species, and supports the low end of the experimental values for H(2)SO(4) and FSO(3)H. Use of a larger basis set (aug-cc-pV5Z) in the CBS extrapolation improves the agreement with experiment for both H(2)SO(4) and FSO(3)H. The G3(MP2) heats of formation for RSO(3)H molecules tend to be underestimated as compared to the CCSD(T)/CBS approach by 2.5-7.0 kcal/mol. COSMO solvation calculations were used to predict solution free energies and pK(a) values with pK(a)'s up to -17.4. Including the solvation of the proton gives good agreement with experimental pK(a) values in the very acidic regime, whereas it is less reliable for weaker acids. The use of CH(3)CO(2)H and HNO(3) as reference acids in the less acidic and more acidic regions of the scale, respectively, provided improved results to within +/-2 pK(a) units in nearly all cases (+/-3 kcal/mol accuracy).