Correlated intermolecular interaction components from asymptotically corrected Kohn-Sham orbitals

In these years there was considerable interest inunderstanding of intermolecular forces in energetic(explosive) systems[1—3]. The supermolecular approach(SM) is widely adopted for calculating ab initio in-termolecular interactions. Nevertheless, it is unable toprovide physically meaningful interaction contribu-tions such as electrostatic, induction, repulsion anddispersion energies. In contrast, the symmetry-adaptedperturbation theory (SAPT)[4—8] has the ability to de-rive these correlated…     

Intermolecular potentials based on symmetry-adapted perturbation theory with dispersion energies from time-dependent density-functional calculations

Recently, three of us have proposed a method [Phys. Rev. Lett. 91, 33201 (2003)] for an accurate calculation of the dispersion energy utilizing frequency-dependent density susceptibilities of monomers obtained from time-dependent density-functional theory (DFT). In the present paper, we report numerical calculations for the helium, neon, water, and carbon dioxide dimers and show that for a wide range of intermonomer separations, including the van der Waals and short-range repulsion regions, the method provides dispersion energies with accuracies comparable to those that can be achieved using the current most sophisticated wave-function methods. If the dispersion energy is combined with (i) the electrostatic and first-order exchange interaction energies as defined in symmetry-adapted perturbation theory (SAPT) but computed using monomer Kohn-Sham (KS) determinants, and (ii) the induction energy computed using the coupled KS static response theory, (iii) the exchange-induction and exchange-dispersion energies computed using KS orbitals and orbital energies, the resulting method, denoted by SAPT(DFT), produces very accurate total interaction potentials. For the helium dimer, the only system with nearly exact benchmark values, SAPT(DFT) reproduces the interaction energy to within about 2% at the minimum and to a similar accuracy for all other distances ranging from the strongly repulsive to the asymptotic region. For the remaining systems investigated by us, the quality of the SAPT(DFT) interaction energies is so high that these energies may actually be more accurate than the best available results obtained with wave-function techniques. At the same time, SAPT(DFT) is much more computationally efficient than any method previously used for calculating the dispersion and other interaction energy components at this level of accuracy.     

Comparison of intermolecular interaction energies from SAPT and DFT including empirical dispersion contributions

The dispersion correction based on damped atom-atom long-range interaction contributions has been tested for an extended S22 database of intermolecular complexes using density functional theory (DFT) and symmetry adapted perturbation theory (SAPT) to account for the remaining interaction energy contributions. In the case of DFT, the dispersion correction of Grimme (J. Comput. Chem. 2006, 27, 1787) was used, while for SAPT, another damping function has been developed that has been optimized particularly for the database. It is found that both approaches yield about the same accuracy for the mixed-type complexes, while the DFT plus dispersion method performs better for the hydrogen-bridged systems and the SAPT plus dispersion approach is better for the dispersion-dominated complexes if compared with coupled cluster singles-doubles with perturbative triples interaction energies as a reference.     

Assessing the Accuracy of SAPT(DFT) Interaction Energies by Comparison with Experimentally Derived Noble Gas Potentials and Molecular Crystal Lattice Energies

The density functional version of symmetry‐adapted perturbation theory, SAPT(DFT), is a computationally efficient method for calculating intermolecular interaction energies. We evaluate its accuracy by comparison with experimentally determined noble gas interaction potentials and sublimation enthalpies, most of which have not been previously calculated using this method. In order to compare the results with wavefunction methods, we also calculate these quantities using MP2 and, for noble gas dimers, using CCSD(T). For the crystal lattice energy calculations, we include corrections to the dispersion, electrostatic, and induction energies that account for the finite interaction distance cutoff and higher‐order induction contributions. Overall, the energy values extrapolated to the complete basis set limit show that SAPT(DFT) achieves significantly better agreement with experiment than MP2.     

Density-functional theory-symmetry-adapted intermolecular perturbation theory with density fitting: a new efficient method to study intermolecular interaction energies

The previously developed DFT-SAPT approach, which combines symmetry-adapted intermolecular perturbation theory (SAPT) with a density-functional theory (DFT) representation of the monomers, has been implemented by using density fitting of two-electron objects. This approach, termed DF-DFT-SAPT, scales with the fifth power of the molecular size and with the third power upon increase of the basis set size for a given dimer, thus drastically reducing the cost of the conventional DFT-SAPT method. The accuracy of the density fitting approximation has been tested for the ethyne dimer. It has been found that the errors in the interaction energies due to density fitting are below 10(-3) kcal/mol with suitable auxiliary basis sets and thus one or two orders of magnitude smaller than the errors due to the use of a limited atomic orbital basis set. An investigation of three prominent structures of the benzene dimer, namely, the T shaped, parallel displaced, and sandwich geometries, employing basis sets of up to augmented quadruple-zeta quality shows that DF-DFT-SAPT outperforms second-order Moller-Plesset theory (MP2) and gives total interaction energies which are close to the best estimates inferred from combining the results of MP2 and coupled-cluster theory with single, double, and perturbative triple excitations.     

Potential energy surface and second virial coefficient of methane-water from ab initio calculations

Six-dimensional intermolecular potential energy surfaces (PESs) for the interaction of CH4 with H2O are presented, obtained from ab initio calculations using symmetry-adapted perturbation theory (SAPT) at two different levels of intramonomer correlation and the supermolecular approach at three different levels of electron correlation. Both CH4 and H2O are assumed to be rigid molecules with interatomic distances and angles fixed at the average values in the ground-state vibration. A physically motivated analytical expression for each PES has been developed as a sum of site-site functions. The PES of the CH4-H2O dimer has only two symmetry-distinct minima. From the SAPT calculations, the global minimum has an energy of -1.03 kcal/mol at a geometry where H2O is the proton donor, HO-H...CH4, with the O-H-C angle of 165 degrees, while the secondary minimum, with an energy of -0.72 kcal/mol, has CH4 in the role of the proton donor (H3C-H...OH2). We estimated the complete basis set limit of the SAPT interaction energy at the global minimum to be -1.06 kcal/mol. The classical cross second virial coefficient B12(T) has been calculated for the temperature range 298-653 K. Our best results agree well with some experiments, allowing an evaluation of the quality of experimental results.     

Highly Accurate CCSD(T) and DFT–SAPT Stabilization Energies of H‐Bonded and Stacked Structures of the Uracil Dimer

The CCSD(T) interaction energies for the H‐bonded and stacked structures of the uracil dimer are determined at the aug‐cc‐pVDZ and aug‐cc‐pVTZ levels. On the basis of these calculations we can construct the CCSD(T) interaction energies at the complete basis set (CBS) limit. The most accurate energies, based either on direct extrapolation of the CCSD(T) correlation energies obtained with the aug‐cc‐pVDZ and aug‐cc‐pVTZ basis sets or on the sum of extrapolated MP2 interaction energies (from aug‐cc‐pVTZ and aug‐cc‐pVQZ basis sets) and extrapolated ΔCCSD(T) correction terms [difference between CCSD(T) and MP2 interaction energies] differ only slightly, which demonstrates the reliability and robustness of both techniques. The latter values, which represent new standards for the H‐bonding and stacking structures of the uracil dimer, differ from the previously published data for the S22 set by a small amount. This suggests that interaction energies of the S22 set are generated with chemical accuracy. The most accurate CCSD(T)/CBS interaction energies are compared with interaction energies obtained from various computational procedures, namely the SCS–MP2 (SCS: spin‐component‐scaled), SCS(MI)–MP2 (MI: molecular interaction), MP3, dispersion‐augmented DFT (DFT–D), M06–2X, and DFT–SAPT (SAPT: symmetry‐adapted perturbation theory) methods. Among these techniques, the best results are obtained with the SCS(MI)–MP2 method. Remarkably good binding energies are also obtained with the DFT–SAPT method. Both DFT techniques tested yield similarly good interaction energies. The large magnitude of the stacking energy for the uracil dimer, compared to that of the benzene dimer, is explained by attractive electrostatic interactions present in the stacked uracil dimer. These interactions force both subsystems to approach each other and the dispersion energy benefits from a shorter intersystem separation.     

Intermolecular Interactions in Photodamaged DNA from Density Functional Theory Symmetry‐Adapted Perturbation Theory

The intermolecular interactions of the photodamaged cyclobutane pyrimidine dimer (CPD) lesion with adjacent nucleobases in the native intrahelical DNA double strand are investigated at the level of density functional theory symmetry‐adapted perturbation theory (DFT‐SAPT) and compared to the original (or repaired) case with pyrimidines (TpT) instead of CPD. The CPD aggregation is on average destabilized by about 6 kcal mol^?1 relative to that involving TpT. The effect of destabilization is asymmetric, that is, it involves a single H‐bonding (Watson–Crick (WC) type) base‐pair interaction.     

Ab initio calculations on low-energy conformers of alpha-cyclodextrin

In this article we carried out a comprehensive investigation of true minima on the potential energy surface (PES) for the alpha-cyclodextrin molecule using ab initio Hartree-Fock (HF) and density functional theory (DFT) quantum chemical methods, employing basis sets ranging from 6-31 G(d,p) to 6-311++G(2d,2p) triple-zeta quality. Thermodynamic quantities and the solvent effect were evaluated at the DFT level of theory. We believe that the most relevant conformers present on the multidimensional PES were sampled in our work, using an adequate treatment of electron correlation effects to describe the intramolecular hydrogen bonds that are present in cyclodextrin species. We present new structures not reported so far and discuss, in detail, the relevance of the DFT gas-phase equilibrium structures for the experimental and theoretical studies involving cyclodextrins and corresponding inclusion complexes, in the condensed phase. In addition, among the various true minimum energy structures located on the DFT PES, the preferred structures in the gas phase and aqueous media, needed to be used as representative minima on the PES in further studies involving the interaction of alpha-cyclodextrin with other species, were unambiguously identified.     

Theoretical investigation of the interactions between the π‐systems of molecular organic semiconductors and an analysis of the contributions of repulsion and electrostatics

A concept for the interactions between π‐systems is necessary to understand a number of phenomena in modern material sciences such as supramolecular properties and self‐assembly. In the present article, we investigate the intermolecular interaction energies between organic semiconductors with extended π‐systems using SAPT (symmetry‐adapted perturbation theory), LMO‐EDA (localized molecular orbital energy decomposition analysis), DFT‐D (density functional theory including dispersion corrections), and force‐field approaches. Both apolar organic molecules such as acenes and highly polarized π‐systems of merocyanines and squaraines were used to probe the influence of electrostatics on the shape of the potential energy surfaces (PES) governing the geometric structures of aggregates. Our results reveal that the shapes of the PESs result from variations in the short‐range, highly specific repulsion forces even for highly polar molecules. Using distributed quadrupoles, we show that it is nevertheless possible to mimic the intermolecular potentials with electrostatics. This is also possible with van‐der‐Waals potentials and a simple overlap‐based force‐field ansatz based on the overlap between p‐orbitals. © 2016 Wiley Periodicals, Inc.     

Pair potential for helium from symmetry-adapted perturbation theory calculations and from supermolecular data

Symmetry-adapted perturbation theory (SAPT) was applied to the helium dimer for interatomic separations R from 3 to 12 bohrs. The first-order interaction energy and the bulk of the second-order contribution were obtained using Gaussian geminal basis sets and are converged to about 0.1 mK near the minimum and for larger R. The remaining second-order contributions available in the SAPT suite of codes were computed using very large orbital basis sets, up to septuple-zeta quality, augmented by diffuse and midbond functions. The accuracy reached at this level was better than 1 mK in the same region. All the remaining components of the interaction energy were computed using the full configuration interaction method in bases up to sextuple-zeta quality. The latter components, although contributing only 1% near the minimum, have the largest uncertainty of about 10 mK in this region. The total interaction energy at R=5.6 bohrs is -11.000+/-0.011 K. For R< or =6.5 bohrs, the supermolecular (SM) interaction energies computed by us recently turned out to be slightly more accurate. Therefore, we have combined the SM results for R< or =6.5 bohrs with the SAPT results from 7.0 to 12 bohrs to fit analytic functions for the potential and for its error bars. The potential fit uses the best available van der Waals constants C(6) through C(16), including C(11), C(13), and C(15), and is believed to be the best current representation of the Born-Oppenheimer (BO) potential for helium. Using these fits, we found that the BO potential for the helium dimer exhibits the well depth D(e)=11.006+/-0.004 K, the equilibrium distance R(e)=5.608+/-0.012 bohrs, and supports one bound state for (4)He(2) with the dissociation energy D(0)=1.73+/-0.04 mK, and the average interatomic separation R=45.6+/-0.5 A.     

Density functional theory guided Monte Carlo simulations: application to melting of Na13

We present a density functional theory (DFT) based Monte Carlo simulation method in which a simple energy function gets fitted on-the-fly to DFT energies and gradients. The fitness of the energy function gets tested periodically using the classical importance function technique [R. Iftimie, D. Salahub, D. Wei, and J. Schofield, J. Chem. Phys. 113, 4852 (2000)]. The function is updated to fit the DFT energies and gradients of the most recent structures visited whenever it fails to achieve a preset accuracy. In this way, we effectively break down the problem of fitting the entire potential energy surface (PES) into many easier problems, which are to fit small local regions of the PES. We used the scaled Morse potential empirical function to guide a DFT Monte Carlo simulation of Na(13) at various temperatures. The use of empirical function guide produced a computational speed-up of about 7 in our test system without affecting the quality of the results.     

Six-dimensional potential energy surface of the dissociative chemisorption of HCl on Au(111) using neural networks

We constructed a six-dimensional potential energy surface(PES)for the dissociative chemisorption of HCl on Au(111)using the neural networks method based on roughly 70000 energies obtained from extensive density functional theory(DFT)calculations.The resulting PES is accurate and smooth,based on the small fitting errors and good agreement between the fitted PES and the direct DFT calculations.Time-dependent wave packet calculations show that the potential energy surface is very well converged with respect to the number of DFT data points,as well as to the fitting process.The dissociation probabilities of HCl initially in the ground rovibrational state from six-dimensional quantum dynamical calculations are quite diferent from the four-dimensional fixed-site calculations,indicating it is essential to perform full-dimensional quantum dynamical studies for the title molecule-surface interaction system.