Method For Intermolecular Interactions Research Articles

Benchmark comparison of dual-basis double-hybrid density functional theory and a neural-network-optimized method for intermolecular interactions

We present a computationally efficient implementation of double-hybrid density functional theory (DH-DFT) leveraging the dual basis methods of Head-Gordon and co-workers and the resolution-of-the-identity second-order Møller-Plesset (RI-MP2) theory. The B2PLYP, B2GP-PLYP, DSD-BLYP and DSD-PBEP86 density functionals are applied to assess the performance of dual-basis methods on several benchmark test cases, including the CONF set of conformational energy differences in C4-C7 alkanes, the S22 set of noncovalent interaction energies, and the RGC10 noble-gas dimer dissociation curves. The dual-basis DH-DFT approach is shown to give results in excellent agreement with conventional methods at a reduced computational cost. For noncovalent interaction energies, DH-DFT is compared against a leading neural-network-based approach, namely the SNS-MP2 method of McGibbon and coworkers (McGibbon et al., 2017). The DH-DFT and SNS-MP2 methods are shown to produce similar accuracies when compared to the established benchmark values.

Advances in molecular quantum chemistry contained in the Q-Chem 4 program package

A summary of the technical advances that are incorporated in the fourth major release of the Q-Chem quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and open-shell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Møller–Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly correlated Cr2 dimer, exploring zeolite-catalysed ethane dehydrogenation, energy decomposition analysis of a charged ter-molecular complex arising from glycerol photoionisation, and natural transition orbitals for a Frenkel exciton state in a nine-unit model of a self-assembling nanotube.

BSSE-free SCF methods for intermolecular interactions

Two versions of the intermolecular SCF theory based on the “chemical Hamiltonian approach” are presented. They permit avoiding the appearance of BSSE in the calculations, so that no a posteriori correction is necessary and one gets BSSE-free energies at the expense of a single (˜ N4) calculation. The results of the two versions are practically equal; they numerically are close to those obtained by the Boys–Bernardi correction scheme but do not exhibit the “overcorrection” behavior of the latter.

Scope and limitations of the SCS-MP2 method for stacking and hydrogen bonding interactions

Fluorobenzenes are pi-acceptor synthons that form pi-stacked structures in molecular crystals as well as in artificial DNAs. We investigate the competition between hydrogen bonding and pi-stacking in dimers consisting of the nucleobase mimic 2-pyridone (2PY) and all fluorobenzenes from 1-fluorobenzene to hexafluorobenzene (n-FB, with n = 1-6). We contrast the results of high level ab initio calculations with those obtained using ultraviolet (UV) and infrared (IR) laser spectroscopy of isolated and supersonically cooled dimers. The 2PY.n-FB complexes with n = 1-5 prefer double hydrogen bonding over pi-stacking, as diagnosed from the UV absorption and IR laser depletion spectra, which both show features characteristic of doubly H-bonded complexes. The 2-pyridone.hexafluorobenzene dimer is the only pi-stacked dimer, exhibiting a homogeneously broadened UV spectrum and no IR bands characteristic for H-bonded species. MP2 (second-order Møller-Plesset perturbation theory) calculations overestimate the pi-stacked dimer binding energies by about 10 kJ/mol and disagree with the experimental observations. In contrast, the MP2 treatment of the H-bonded dimers appears to be quite accurate. Grimme's spin-component-scaled MP2 approach (SCS-MP2) is an improvement over MP2 for the pi-stacked dimers, reducing the binding energy by approximately 10 kJ/mol. When applied to explicitly correlated MP2 theory (SCS-MP2-R12 approach), agreement with the corresponding coupled-cluster binding energies [at the CCSD(T) level] is very good for the pi-stacked dimers, within +/- 1 kJ/mol for the 2PY complexes with 1-fluorobenzene, 1,2-difluorobenzene, 1,2,4,5-tetrafluorobenzene, pentafluorobenzene and hexafluorobenzene. Unfortunately, the SCS-MP2 approach also reduces the binding energy of the H-bonded species, leading to disagreement with both coupled-cluster theory and experiment. The SCS-MP2-R12 binding energies follow the SCS-MP2 binding energies closely, being about 0.5 and 0.7 kJ/mol larger for the H-bonded and pi-stacked forms, respectively, in an augmented correlation-consistent polarized valence quadruple-zeta basis. It seems that the SCS-MP2 and SCS-MP2-R12 methods cannot provide sufficient accuracy to replace the CCSD(T) method for intermolecular interactions where H-bonding and pi-stacking are competitive.