A density-fitting implementation of the density-based basis-set correction method.

This work reports an efficient density-fitting implementation of the density-based basis-set correction (DBBSC) method in the MOLPRO software. This method consists in correcting the energy calculated by a wave-function method with a given basis set by an adapted basis-set correction density functional incorporating the short-range electron correlation effects missing in the basis set, resulting in an accelerated convergence to the complete-basis-set limit. Different basis-set correction density-functional approximations are explored and the complementary-auxiliary-basis-set single-excitation correction is added. The method is tested on a benchmark set of reaction energies at the second-order Møller-Plesset (MP2) level and a comparison with the explicitly correlated MP2-F12 method is provided. The results show that the DBBSC method greatly accelerates the basis convergence of MP2 reaction energies, without reaching the accuracy of the MP2-F12 method but with a lower computational cost.

Thermochemical evaluation of adaptive and fixed density functional theory quadrature schemes.

A systematic study is made of the accuracy and efficiency of a number of existing quadrature schemes for molecular Kohn-Sham Density-Functional Theory (DFT) using 408 molecules and 254 chemical reactions. Included are the fixed SG-x (x = 0-3) grids of Gill et al., Dasgupta, and Herbert, the 3-zone grids of Treutler and Ahlrichs, a fixed five-zone grid implemented in Molpro, and a new adaptive grid scheme. While all methods provide a systematic reduction of errors upon extension of the grid sizes, significant differences are observed in the accuracies for similar grid sizes with various approaches. For the tests in this work, the SG-x fixed grids are less suitable to achieve high accuracies in the DFT integration, while our new adaptive grid performed best among the schemes studied in this work. The extra computational time to generate the adaptive grid scales linearly with molecular size and is negligible compared with the time needed for the self-consistent field iterations for large molecules. A comparison of the grid accuracies using various density functionals shows that meta-GGA functionals need larger integration grids than GGA functionals to reach the same degree of accuracy, confirming previous investigations of the numerical stability of meta-GGA functionals. On the other hand, the grid integration errors are almost independent of the basis set, and the basis set errors are mostly much larger than the errors caused by the numerical integrations, even when using the smallest grids tested in this work.

The Molpro quantum chemistry package.

Molpro is a general purpose quantum chemistry software package with a long development history. It was originally focused on accurate wavefunction calculations for small molecules but now has many additional distinctive capabilities that include, inter alia, local correlation approximations combined with explicit correlation, highly efficient implementations of single-reference correlation methods, robust and efficient multireference methods for large molecules, projection embedding, and anharmonic vibrational spectra. In addition to conventional input-file specification of calculations, Molpro calculations can now be specified and analyzed via a new graphical user interface and through a Python framework.

Intermolecular interaction energies from fourth order many-body perturbation theory. Impact of individual electron correlation contributions.

The performance of Møller-Plesset perturbation theory methods for describing intermolecular interaction energies has been investigated with the focus on illuminating the impact of individual electron correlation energy contributions in fourth order. It is shown that a physically meaningful decomposition of the fourth order correlation energy can be obtained by grouping individual correlation energy terms that share the same diagrammatic loop structure. This decomposition of the fourth order singles (S), doubles (D), triples (T), and quadruples (Q) terms revealed that individual terms from each excitation class can have a huge impact on the energy that is much larger than the total fourth order correlation contribution. A partial summation of S, D, and Q terms has been derived that can reproduce the full fourth order interaction energies with a good accuracy and which does not include the computationally expensive triples energy term.

The coulombic σ-hole model describes bonding in CX3IY- complexes completely.

Contrary to recent reports, the σ-hole interaction energies of complexes between the carbon tetrahalides CX3I (X = F, Cl, Br, I) and halide anions Y- (Y = F, Cl, Br, I) are described very well by the simple Coulombic σ-hole concept if it is applied properly. There is no need to invoke charge transfer, which in any case is not uniquely distinguishable from polarization.

Geometry optimisations with a nonlocal density-functional theory method based on a double Hirshfeld partitioning.

Energy gradients have been derived for the nonlocal density-functional theory (NLDFT) method from Heßelmann [J. Chem. Theory Comput. 9, 273 (2013)]. It is shown that the derivative of the NLDFT functional can easily be obtained analytically due to the fact that the inherent Hirshfeld weights are described in terms of analytic expressions of the atomic densities determined by Slater's rules. The accuracy of the NLDFT gradient has been tested by performing geometry optimisations for a range of 76 tripeptide molecules and a number of small noncovalently bonded dimer systems for which high level coupled cluster reference structures are accessible. It was found that the resulting optimised structures are in good agreement with corresponding structures optimised using second-order Møller-Plesset or coupled cluster wave function electron correlation methods. Moreover, conformer energies as well as intermolecular interaction energies are shown to be in fair agreement with corresponding density-functional theory methods employing pairwise atom-atom dispersion models.

Correlation effects and many-body interactions in water clusters.

Background: The quantum-chemical description of the interactions in water clusters is an essential basis for deriving accurate and physically sound models of the interaction potential for water to be used in molecular simulations. In particular, the role of many-body interactions beyond the two-body interactions, which are often not explicitly taken into account by empirical force fields, can be accurately described by quantum chemistry methods on an adequate level, e.g., random-phase approximation electron correlation methods. The relative magnitudes of the different interaction energy contributions obtained by accurate ab initio calculations can therefore provide useful insights that can be exploited to develop enhanced force field methods. Results: In line with earlier theoretical studies of the interactions in water clusters, it has been found that the main contribution to the many-body interactions in clusters with a size of up to N = 13 molecules are higher-order polarisation interaction terms. Compared to this, many-body dispersion interactions are practically negligible for all studied sytems. The two-body dispersion interaction, however, plays a significant role in the formation of the structures of the water clusters and their stability, since it leads to a distinct compression of the cluster sizes compared to the structures optimized on an uncorrelated level. Overall, the many-body interactions amount to about 13% of the total interaction energy, irrespective of the cluster size. The electron correlation contribution to these, however, amounts to only about 30% to the total many-body interactions for the largest clusters studied and is repulsive for all structures considered in this work. Conclusion: While this shows that three- and higher-body interactions can not be neglected in the description of water complexes, the electron correlation contributions to these are much smaller in comparison to the two-body electron correlation effects. Efficient quantum chemistry approaches for describing intermolecular interactions between water molecules may therefore describe higher-body interactions on an uncorrelated Hartree-Fock level without a serious loss in accuracy.

Intramolecular interactions in sterically crowded hydrocarbon molecules.

A molecular fragmentation method has been used to analyze the intramolecular interactions in the three molecules coupled diamantane, hexaphenylethane, and all-meta-tert-butyl substituted hexaphenylethane. The significance of these systems lies in the fact, that steric crowding effects enable a stabilization of the central carbon bond that possesses an extended length (1.6 to 1.7 Å) beyond conventional carbon-carbon bonds due to the steric repulsion of the attached hydrocarbon groups. The total stability of these molecules therefore depends on a delicate balance between attractive interaction forces on the one hand and on repulsive forces on the other hand. We have quantified the different interaction energy contributions using symmetry-adapted perturbation theory based on a density functional theory description of the monomers. It has been found that the attractive dispersion interactions increase more strongly with the level of crowding in the systems than the counteracting exchange interactions. This shows that steric crowding effects can have a significant impact on the structure and stability of large and branched molecules. © 2017 Wiley Periodicals, Inc.

Trifluoromethyl: An Amphiphilic Noncovalent Bonding Partner.

The traditional "Fδ- " picture of fluorine suggests that it can only interact with electrophilic centers such as backbone-carbonyl carbon atoms or hydrogen-bond donors in proteins. We show that this view, which neglects polarization, is incomplete and the trifluoromethyl groups can act both as electrophiles and nucleophiles to form noncovalent interactions. The underlying polarization mechanism is based on the anomeric effect and is only fully operative if the geometry is allowed to relax. MP2/aug-cc-pVDZ calculations on model systems demonstrate the effect of the unusual group polarizability of trifluoromethyl. A survey of the Protein Databank reveals more than 600 weak interactions involving a trifluorotoluene moiety. The unique combination of the anomeric effect and the group-polarization process associated with it in CF3 allows its most negative molecular electrostatic potential (MEP) on the surface in contact with a nucleophile to become zero, so that the area of positive MEP on the backside of the carbon atom becomes dominant. However, the reverse polarization is also facile, so that CF3 can also act as an H-bond acceptor for cations such as the guanidinium group of arginine.

Low scaling random-phase approximation electron correlation method including exchange interactions using localised orbitals.

A random-phase approximation electron correlation method including exchange interactions has been developed which reduces the scaling behaviour of the standard approach by two to four orders of magnitude, effectively leading to a linear scaling performance if the local structures of the underlying quantities are fully exploited in the calculations. This has been achieved by a transformation of the integrals and amplitudes from the canonical orbital basis into a local orbital basis and a subsequent dyadic screening approach. The performance of the method is demonstrated for a range of tripeptide molecules as well as for two conformers of the polyglycine molecule using up to 40 glycine units. While a reasonable agreement with the corresponding canonical method is obtained if long-range Coulomb interactions are not screened by the local method, a significant improvement in the performance is achieved for larger systems beyond 20 glycine units. Furthermore, the control of the Coulomb screening threshold allows for a quantification of intramolecular dispersion interactions, as will be exemplified for the polyglycine conformers as well as a highly branched hexaphenylethane derivate which is stabilised by steric crowding effects.

On the Stability of Cyclophane Derivates Using a Molecular Fragmentation Method.

A molecular fragmentation method is used to study the stability of cyclophane derivates by decomposing the molecular energy into the molecular strain and intramolecular interaction energies. The molecular strain energies obtained by utilising the fragmentation method are in good agreement with existing experimental data. The intramolecular interaction energies calculated as the difference between the supermolecular energy and the bonded fragment energies are repulsive in the cyclophanes studied. The nature of this interaction is studied for groups of systematically extended doubled layered paracyclophane systems using the random-phase approximation (RPA), two recently developed extensions to the RPA and standard density functional theory (DFT) methods including dispersion corrections. Upon a systematic increase in conjugation the strongly repulsive intramolecular interaction energy reduces and thus leads to an increase in the stability. Finally, existing experimental and theoretical estimates of the molecular strain are compared with the results of this work.

Molecular energies from an incremental fragmentation method.

The systematic molecular fragmentation method by Collins and Deev [J. Chem. Phys. 125, 104104 (2006)] has been used to calculate total energies and relative conformational energies for a number of small and extended molecular systems. In contrast to the original approach by Collins, we have tested the accuracy of the fragmentation method by utilising an incremental scheme in which the energies at the lowest level of the fragmentation are calculated on an accurate quantum chemistry level while lower-cost methods are used to correct the low-level energies through a high-level fragmentation. In this work, the fragment energies at the lowest level of fragmentation were calculated using the random-phase approximation (RPA) and two recently developed extensions to the RPA while the incremental corrections at higher levels of the fragmentation were calculated using standard density functional theory (DFT) methods. The complete incremental fragmentation method has been shown to reproduce the supermolecule results with a very good accuracy, almost independent on the molecular type, size, or type of decomposition. The fragmentation method has also been used in conjunction with the DFT-SAPT (symmetry-adapted perturbation theory) method which enables a breakdown of the total nonbonding energy contributions into individual interaction energy terms. Finally, the potential problems of the method connected with the use of capping hydrogen atoms are analysed and two possible solutions are supplied.

Polarisabilities of long conjugated chain molecules with density functional response methods: The role of coupled and uncoupled response.

The longitudinal component of the dipole-dipole polarisability of polyacetylene molecules containing 4 to 20 carbon atoms has been calculated with density-functional theory (DFT) response methods. In order to analyse the effect of the uncoupled and coupled contributions to the response matrix, a number of different sets of orbitals were combined with different approximations for the Hessian matrix. This revealed a surprising result: a qualitatively correct increase of the polarisability with the chain length can already be reproduced on the uncoupled level if the response matrix is constructed from Hartree-Fock (HF) or exact-exchange (EXX) DFT orbitals. The nonlocal HF and the local EXX exchange potentials both produce a displacement of charge from the chain ends to the centre of the polyacetylene molecule compared to DFT methods using standard exchange-correlation potentials. In this way, the reduced increase of the transition dipole moments along the molecular axis counteracts the decrease of the occupied-virtual orbital energy gaps and leads to a linear dependence of the polarisabilities (normalised by the number of carbon atoms) on the chain length. A new DFT response approach is tested which utilises unitary transformed Hartree-Fock orbitals as input and which resolves the failure of standard DFT response methods.

Intermolecular symmetry-adapted perturbation theory study of large organic complexes.

Binding energies for the complexes of the S12L database by Grimme [Chem. Eur. J. 18, 9955 (2012)] were calculated using intermolecular symmetry-adapted perturbation theory combined with a density-functional theory description of the interacting molecules. The individual interaction energy decompositions revealed no particular change in the stabilisation pattern as compared to smaller dimer systems at equilibrium structures. This demonstrates that, to some extent, the qualitative description of the interaction of small dimer systems may be extrapolated to larger systems, a method that is widely used in force-fields in which the total interaction energy is decomposed into atom-atom contributions. A comparison of the binding energies with accurate experimental reference values from Grimme, the latter including thermodynamic corrections from semiempirical calculations, has shown a fairly good agreement to within the error range of the reference binding energies.

Efficient self-consistent treatment of electron correlation within the random phase approximation.

A self-consistent Kohn-Sham (KS) method is presented that treats correlation on the basis of the adiabatic-connection dissipation-fluctuation theorem employing the direct random phase approximation (dRPA), i.e., taking into account only the Coulomb kernel while neglecting the exchange-correlation kernel in the calculation of the Kohn-Sham correlation energy and potential. The method, denoted self-consistent dRPA method, furthermore treats exactly the exchange energy and the local multiplicative KS exchange potential. It uses Gaussian basis sets, is reasonably efficient, exhibiting a scaling of the computational effort with the forth power of the system size, and thus is generally applicable to molecules. The resulting dRPA correlation potentials in contrast to common approximate correlation potentials are in good agreement with exact reference potentials. The negatives of the eigenvalues of the highest occupied molecular orbitals are found to be in good agreement with experimental ionization potentials. Total energies from self-consistent dRPA calculations, as expected, are even poorer than non-self-consistent dRPA total energies and dRPA reaction and non-covalent binding energies do not significantly benefit from self-consistency. On the other hand, energies obtained with a recently introduced adiabatic-connection dissipation-fluctuation approach (EXXRPA+, exact-exchange random phase approximation) that takes into account, besides the Coulomb kernel, also the exact frequency-dependent exchange kernel are significantly improved if evaluated with orbitals obtained from a self-consistent dRPA calculation instead of an exact exchange-only calculation. Total energies, reaction energies, and noncovalent binding energies obtained in this way are of the same quality as those of high-level quantum chemistry methods, like the coupled cluster singles doubles method which is computationally more demanding.

Long-range correlation energies from frequency-dependent weighted exchange-hole dipole polarisabilities.

Long-range correlation energies are calculated using an approximation of the single-particle density-density response function of the system that leads to an expression requiring only occupied orbitals and eigenvalues. Dipole-dipole polarisabilities and isotropic leading-order dispersion coefficients obtained from this approximation are shown to be in a reasonable agreement with corresponding values from the experiment or dipole oscillator strength distributions. The localised polarisabilities were used to calculate a long-range correlation correction to a hybrid-generalised gradient approximation functional using a proper damping function at short ranges. It was found that the hybrid density-functional theory+dispersion method obtained in this way has a comparable accuracy than high-level ab initio wave function methods at a much lower computational cost. This has been analysed for a number of systems from the GMTKN30 database including subsets for noncovalently bound complexes, relative energies for sugar conformers and reaction energies and barrier heights of pericyclic reactions of some medium sized organic molecules.

Resolution of identity approach for the Kohn-Sham correlation energy within the exact-exchange random-phase approximation.

Two related methods to calculate the Kohn-Sham correlation energy within the framework of the adiabatic-connection fluctuation-dissipation theorem are presented. The required coupling-strength-dependent density-density response functions are calculated within exact-exchange time-dependent density-functional theory, i.e., within time-dependent density-functional response theory using the full frequency-dependent exchange kernel in addition to the Coulomb kernel. The resulting resolution-of-identity exact-exchange random-phase approximation (RI-EXXRPA) methods in contrast to previous EXXRPA methods employ an auxiliary basis set (RI basis set) to improve the computational efficiency, in particular, to reduce the formal scaling of the computational effort with respect to the system size N from N(6) to N(5). Moreover, the presented RI-EXXRPA methods, in contrast to previous ones, do not treat products of occupied times unoccupied orbitals as if they were linearly independent. Finally, terms neglected in previous EXXRPA methods can be included, which leads to a method designated RI-EXXRPA+, while the method without these extra terms is simply referred to as RI-EXXRPA. Both EXXRPA methods are shown to yield total energies, reaction energies of small molecules, and binding energies of noncovalently bonded dimers of a quality that is similar and in some cases even better than that obtained with quantum chemistry methods such as Mo̸ller-Plesset perturbation theory of second order (MP2) or with the coupled cluster singles doubles method. In contrast to MP2 and to conventional density-functional methods, the presented RI-EXXRPA methods are able to treat static correlation.

Correct description of the bond dissociation limit without breaking spin symmetry by a random-phase-approximation correlation functional.

A correlation functional that is termed exact-exchange random phase approximation (EXX-RPA) functional and is obtained with the exact frequency-dependent exchange kernel via the fluctuation-dissipation theorem is shown to correctly describe electron pair bonds in the dissociation limit without the need to resort to symmetry breaking in spin space. Because the functional also yields more accurate electronic energies for molecules in their equilibrium geometry than standard correlation functionals, it combines accuracy at equilibrium bond distances and in dissociation processes with a correct description of spin, something all commonly employed correlation functionals fail to do. The reason why the EXX-RPA correlation functional yields distinctively and qualitatively better results than RPA approaches based on Hartree-Fock and time-dependent Hartree-Fock is explained.

On the accuracy of DFT-SAPT, MP2, SCS-MP2, MP2C, and DFT+Disp methods for the interaction energies of endohedral complexes of the C(60) fullerene with a rare gas atom.

Selected points on the potential energy surface for the complexes Rg@C(60) (Rg = He, Ne, Ar, Kr) are calculated with various theoretical methods, like symmetry-adapted perturbation theory with monomers described by density functional theory (DFT-SAPT), supermolecular Møller-Plesset theory truncated on the second order (MP2), spin-component-scaled MP2 (SCS-MP2), supermolecular density functional theory with empirical dispersion correction (DFT+Disp), and the recently developed MP2C method that improves the MP2 method for long-range electron correlation effects. A stabilization of the endohedral complex is predicted by all methods, but the depth of the potential energy well is overestimated by the DFT+Disp and MP2 approaches. On the other hand, the MP2C model agrees well with DFT-SAPT, which serves as the reference. The performance of SCS-MP2 is mixed: it produces too low interaction energies for the two heavier guests, while its accuracy for He@C(60) and Ne@C(60) is similar to that of MP2C. Fitting formulas for the main interaction energy components, i.e. the dispersion and first-order repulsion energies are proposed, which are applicable for both endo- and exohedral cases. For all examined methods density fitting is used to evaluate two-electron repulsion integrals, which is indispensable to allow studies of noncovalent complexes of this size. It has been found that density-fitting auxiliary basis sets cannot be used in a black-box fashion for the calculation of the first-order SAPT electrostatic energy, and that the quality of these basis sets should be always carefully examined in order to avoid an unphysical long-range behavior.

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.