Abstract:Although covalent interactions determine the primary structure of a molecule, the noncovalent interactions are responsible for the tertiary and quaternary structure of a molecule and create the fascinating world of the 3D architectures of biomacromolecules. For example, the double helical structure of DNA is of fundamental importance for the function of DNA: it allows it to store and transfer genetic information. To fulfill this role, the structure is rigid to maintain the double helix with a proper positioning of the complementary base, and floppy to allow for its opening. Very strong covalent interactions cannot fulfill both of these criteria, but noncovalent interactions, which are about 2 orders of magnitude weaker, can. This Account highlights the recent advances in the field of the design of novel wave function theory (WFT) methods applicable to noncovalent complexes ranging in size from less than 100 atoms, for which highly accurate ab initio methods are available, up to extended ones (several thousands atoms), which are the domain of semiempirical QM (SQM) methods. Accurate interaction energies for noncovalent complexes are generated by the coupled-cluster technique, taking single- and double-electron excitations iteratively and triple-electron excitation perturbatively with a complete basis set description (CCSD(T)/CBS). The procedure provides interaction energies with high accuracy (error less than 1 kcal/mol). Because the method is computationally demanding, its application is limited to complexes smaller than 30 atoms. But researchers would also like to use computational methods to determine these interaction energies accurately for larger biological and nanoscale structures. Standard QM methods such as MP2, MP3, CCSD, or DFT fail to describe various types of noncovalent systems (H-bonded, stacked, dispersion-controlled, etc.) with comparable accuracy. Therefore, novel methods are needed that have been parametrized toward noncovalent interactions, and existing benchmark data sets represent an important tool for the development of new methods providing reliable characteristics of noncovalent clusters. Our laboratory developed the first suitable data set of CCSD(T)/CBS interaction energies and geometries of various noncovalent complexes, called S22. Since its publication in 2006, it has frequently been applied in parametrization and/or verification of various wave function and density functional techniques. During the intense use of this data set, several inconsistencies emerged, such as the insufficient accuracy of the CCSD(T) correction term or its unbalanced character, which has triggered the introduction of a new, broader, and more accurate data set called the S66 data set. It contains not only 66 CCSD(T)/CBS interaction energies determined in the equilibrium geometries but also 1056 interaction energies calculated at the same level for nonequilibrium geometries. The S22 and S66 data sets have been used for the verification of various WFT methods, and the lowest RMSE (S66, in kcal/mol) was found for the recently introduced SCS-MI-CCSD/CBS (0.08), MP2.5/CBS (0.16), MP2.X/6-31G* (0.27), and SCS-MI-MP2/CBS (0.38) methods. Because of their computational economy, the MP2.5 and MP2.X/6-31G* methods can be recommended for highly accurate calculations of large complexes with up to 100 atoms. The evaluation of SQM methods was based only on the S22 data set, and because some of these methods have been parametrized toward the same data set, the respective results should be taken with caution. For really extended complexes such as protein-ligand systems, only the SMQ methods are applicable. After adding the corrections to the dispersion energy and H-bonding, several methods exhibit surprisingly low RMSE (even below 0.5 kcal/mol). Among the various SMQ methods, the PM6-DH2 can be recommended because of its computational efficiency and it can be used for optimization (which is not the case for other SQM methods). The PM6-DH2 is the base of our novel scoring function used in in silico drug design.

Ab initio potential energy curves, transition dipole moments and spin–orbit coupling matrix elements for the first twenty states of the calcium dimer

Four-body singlet potential energy surface for reactions of calcium monofluoride

Rovibrational dynamics of the strontium molecule in the A(1)Σ(u)+, c(3)Π(u), and a(3)Σ(u)+ manifold from state-of-the-art ab initio calculations

Electric dipole transition amplitudes for atoms and ions with one valence electron

Improved Direct Diabatization and Coupled Potential Energy Surfaces for the Photodissociation of Ammonia

Changing the paradigm in f-containing cold molecules: the impact of spin-orbit coupling and f-d transitions on quasi-bound vibrational states

The Numerical Equivalence of Diabatic and Adiabatic Representations in Diatomic Molecules

Ab initio calculations on halogen-bonded complexes and comparison with density functional methods.

Re-examination of the Cs2 ground singlet X1Σg+ and triplet a3Σu+ states

Ab initio QED calculations in diatomic quasimolecules

Coupled 3D ( J ≥ 0) Time-Dependent Wave Packet Calculation for the F + H 2 Reaction on Accurate Ab Initio Multi-State Diabatic Potential Energy Surfaces

Quantum Mechanical Calculation of Noncovalent Interactions: A Large-Scale Evaluation of PMx, DFT, and SAPT Approaches

HAB79: A new molecular dataset for benchmarking DFT and DFTB electronic couplings against high-level ab initio calculations

Strictly correlated electrons approach to excitation energies of dissociating molecules

Excited states with selected CI-QMC: chemically accurate excitation energies and geometries

Ab initio investigation of potential energy curves of He2, He2+, and extrapolation by the machine learning method

Calculations on noncovalent interactions and databases of benchmark interaction energies

Extension of the B3LYP - Dispersion-Correcting Potential Approach to the Accurate Treatment of both Inter- and Intramolecular Interactions

Theoretical study of the LiNa molecule beyond the Born–Oppenheimer approximation: adiabatic and diabatic potential energy curves, radial coupling, adiabatic correction, dipole moments and vibrational levels

Benchmarking TD-DFT and Wave Function Methods for Oscillator Strengths and Excited-State Dipole Moments

Re-Examination Of The Cs-2 Ground Singlet X-1 Sigma(+)(G) And Triplet A(3)Sigma(+)(U) States