This work studies the underlying nature of H‐bonds (HBs) of different types and strengths and tries to predict binding energies (BEs) based on the properties derived from wave function analysis. A total of 42 HB complexes constructed from 28 neutral and 14 charged monomers were considered. This set was designed to sample a wide range of HB strengths to obtain a complete view about HBs. BEs were derived with the accurate coupled cluster singles and doubles with perturbative triples correction (CCSD(T))(T) method and the physical components of the BE were investigated by symmetry‐adapted perturbation theory (SAPT). Quantum theory of atoms‐in‐molecules (QTAIM) descriptors and other HB indices were calculated based on high‐quality density functional theory wave functions. We propose a new and rigorous classification of H‐bonds (HBs) based on the SAPT decomposition. Neutral complexes are either classified as “very weak” HBs with a BE ≥ −2.5 kcal/mol that are mainly dominated by both dispersion and electrostatic interactions or as “weak‐to‐medium” HBs with a BE varying between −2.5 and −14.0 kcal/mol that are only dominated by electrostatic interactions. On the other hand, charged complexes are divided into “medium” HBs with a BE in the range of −11.0 to −15.0 kcal/mol, which are mainly dominated by electrostatic interactions, or into “strong” HBs whose BE is more negative than −15.0 kcal/mol, which are mainly dominated by electrostatic together with induction interactions. Among various explored correlations between BEs and wave function‐based HB descriptors, a fairly satisfactory correlation was found for the electron density at the bond critical point (BCP; ρBCP) of HBs. The fitted equation for neutral complexes is BE/kcal/mol =  − 223.08 × ρBCP/a. u. + 0.7423 with a mean absolute percentage error (MAPE) of 14.7%, while that for charged complexes is BE/kcal/mol =  − 332.34 × ρBCP/a. u. − 1.0661 with a MAPE of 10.0%. In practice, these equations may be used for a quick estimation of HB BEs, for example, for intramolecular HBs or large HB networks in biomolecules. © 2019 Wiley Periodicals, Inc.

A semi-empirical counterpoise-type correction for basis set superposition error (BSSE) in molecular systems is presented. An atom pair-wise potential corrects for the inter- and intra-molecular BSSE in supermolecular Hartree-Fock (HF) or density functional theory (DFT) calculations. This geometrical counterpoise (gCP) denoted scheme depends only on the molecular geometry, i.e., no input from the electronic wave-function is required and hence is applicable to molecules with ten thousands of atoms. The four necessary parameters have been determined by a fit to standard Boys and Bernadi counterpoise corrections for Hobza's S66×8 set of non-covalently bound complexes (528 data points). The method's target are small basis sets (e.g., minimal, split-valence, 6-31G*), but reliable results are also obtained for larger triple-ζ sets. The intermolecular BSSE is calculated by gCP within a typical error of 10%-30% that proves sufficient in many practical applications. The approach is suggested as a quantitative correction in production work and can also be routinely applied to estimate the magnitude of the BSSE beforehand. The applicability for biomolecules as the primary target is tested for the crambin protein, where gCP removes intramolecular BSSE effectively and yields conformational energies comparable to def2-TZVP basis results. Good mutual agreement is also found with Jensen's ACP(4) scheme, estimating the intramolecular BSSE in the phenylalanine-glycine-phenylalanine tripeptide, for which also a relaxed rotational energy profile is presented. A variety of minimal and double-ζ basis sets combined with gCP and the dispersion corrections DFT-D3 and DFT-NL are successfully benchmarked on the S22 and S66 sets of non-covalent interactions. Outstanding performance with a mean absolute deviation (MAD) of 0.51 kcal/mol (0.38 kcal/mol after D3-refit) is obtained at the gCP-corrected HF-D3/(minimal basis) level for the S66 benchmark. The gCP-corrected B3LYP-D3/6-31G* model chemistry yields MAD=0.68 kcal/mol, which represents a huge improvement over plain B3LYP/6-31G* (MAD=2.3 kcal/mol). Application of gCP-corrected B97-D3 and HF-D3 on a set of large protein-ligand complexes prove the robustness of the method. Analytical gCP gradients make optimizations of large systems feasible with small basis sets, as demonstrated for the inter-ring distances of 9-helicene and most of the complexes in Hobza's S22 test set. The method is implemented in a freely available FORTRAN program obtainable from the author's website.

Psi4 is a free and open-source ab initio electronic structure program providing Hartree-Fock, density functional theory, many-body perturbation theory, configuration interaction, density cumulant theory, symmetry-adapted perturbation theory, and coupled-cluster theory. Most of the methods are quite efficient thanks to density fitting and multi-core parallelism. The program is a hybrid of C++ and Python, and calculations may be run with very simple text files or using the Python API, facilitating post-processing and complex workflows; method developers also have access to most of Psi4's core functionality via Python. Job specification may be passed using The Molecular Sciences Software Institute (MolSSI) QCSchema data format, facilitating interoperability. A rewrite of our top-level computation driver, and concomitant adoption of the MolSSI QCArchive Infrastructure project, make the latest version of Psi4 well suited to distributed computation of large numbers of independent tasks. The project has fostered the development of independent software components that may be reused in other quantum chemistry programs. File list (2) download file view on ChemRxiv psi4.pdf (4.37 MiB) download file view on ChemRxiv supplementary_material.pdf (297.86 KiB)

Dispersion-corrected density functional theory is assessed on the new S66 and S66x8 benchmark sets for non-covalent interactions. In total, 17 different density functionals are evaluated. Two flavors of our latest additive London-dispersion correction DFT-D3 and DFT-D3(BJ), which differ in their short-range damping functions, are tested. In general, dispersion corrections are again shown to be crucial to obtain reliable non-covalent interaction energies and equilibrium distances. The corrections strongly diminish the performance differences between the functionals, and in summary most dispersion-corrected methods can be recommended. DFT-D3 and DFT-D3(BJ) also yield similar results but for most functionals and intermolecular distances, the rational Becke-Johnson scheme performs slightly better. Particularly, the statistical analysis for S66x8, which covers also non-equilibrium complex geometries, shows that the Minnesota class of functionals is also improved by the D3 scheme. The best methods on the (meta-)GGA or hybrid- (meta-)GGA level are B97-D3, BLYP-D3(BJ), PW6B95-D3, MPW1B95-D3 and LC-ωPBE-D3. Double-hybrid functionals are the most accurate and robust methods, and in particular PWPB95-D3 and B2-PLYP-D3(BJ) can be recommended. The best DFT-D3 and DFT-D3(BJ) approaches are competitive to specially adapted perturbation methods and clearly outperform standard MP2. Comparisons between S66, S22 and parts of the GMTKN30 database show that the S66 set provides statistically well-behaved data and can serve as a valuable tool for, for example, fitting purposes or cross-validation of other benchmark databases.

Activation of dihydrogen is typically a domain of transition metal chemistry.[1] Even nature uses metal-centered reactions to split the dihydrogen molecule in hydrogenase enzymes.[2] A recent development is the use of metal-free systems for H 2 activation: Stephan, Erker et al. have described frustrated Lewis pairs (FLP), that is, pairs of Lewis acids and bases that do not fully quench each other owing to the steric bulk of their substituents; these pairs heterolytically split the H 2 molecule (Scheme 1). [3][4][5][6] Phosphane/borane pairs, such as 1 or 3 (and an increasing number of related systems that have appeared in the literature), react rapidly and effectively with H 2 to yield the corresponding phosphonium cation/hydridoborate anion pairs (here 2 and 4, respectively). These systems have been used as active metal-free hydrogenation catalysts.Pµpai et al. [7] have presented the notion that the HÀH bond is cleaved by such systems in an almost linear P-H-H-B arrangement in the transition state (TS). We have now found that this proposal is probably a gross oversimplification of the mechanistic course taken, because the theoretical treatment used did not adequately take into account the interaction between the large substituents that are specifically used. We wish however to point out that the authors in Ref.[7a] pointed out correctly for the first time the importance of secondary, non-covalent C 6 F 5 ···tBu interactions. Herein, we present the results of our state-of-the-art calculations for this problem, which has resulted in a more realistic description of the TS involved, which features a non-linear P-H-H-B unit. Furthermore, we present an even simpler mechanistic picture of the basic activation step that emphasizes on the polarization of H 2 induced by the electric field of the FLP inside its cavity that can explain important (and hitherto unclear) experimental findings.One of the first and very basic questions involves the structure of the TS and in particular in how far the proposed linear P-H-H-B arrangement is required. For the intramolecular system 3 and the similar case of Sumerin et al.,[5c] a linear TS is geometrically not possible, although these systems also efficiently activate H 2 at ambient temperatures.For molecules 1-4, we performed high-level quantum chemical calculations at wavefunction (WF)-based levels (SCS-MP2[8a] and MP2, extrapolated to the complete basis set [CBS] limit [8b] ) and by state-of-the-art dispersion-corrected density functional theory (DFT-D) using the B97-D functional [8c] (for details, see the Supporting Information). For 1/H 2 , we first computed a relaxed two-dimensional potential energy surface (PES) with a fixed linear P-H-H-B unit with the most important HÀH and PÀB distances as variables. Full TS optimizations were then performed for both 1/H 2 and 3/H 2 . It should be noted that all our computations refer to isolated molecule conditions; although this makes direct comparisons with experimental observations difficult, [9] we think that it is very important to...

We analyze the error compensations that are responsible for the relatively good performance of the popular B3LYP/6-31G* model chemistry for molecular thermochemistry. We present the B3LYP-gCP-D3/6-31G* scheme, which corrects for missing London dispersion and basis set superposition error (BSSE) in a physically sound manner. Benchmark results for the general main group thermochemistry, kinetics, and noncovalent interactions set (GMTKN30) are presented. A detailed look is cast on organic reactions of several arenes with C(60), Diels-Alder reactions, and barriers to [4 + 3] cycloadditions. We demonstrate the practical advantages of the new B3LYP-gCP-D3/6-31G* scheme and show its higher robustness over standard B3LYP/6-31G*. B3LYP-gCP-D3/6-31G* is meant to fully substitute standard B3LYP/6-31G* calculations in the same black-box sense at essentially no increase in computational cost. The energy corrections are made available by a Web service ( http://www.thch.uni-bonn.de/tc/gcpd3 ) and by freely available software.

Explicit solvent atomistic molecular dynamics (MD) simulations represent an established technique to study structural dynamics of RNA molecules and an important complement for diverse experimental methods. However, performance of molecular mechanical (MM) force fields (ff's) remains far from satisfactory even after decades of development, as apparent from a problematic structural description of some important RNA motifs. Actually, some of the smallest RNA molecules belong to the most challenging systems for MD simulations and, among them, the UUCG tetraloop is saliently difficult. We report a detailed analysis of UUCG MD simulations, depicting the sequence of events leading to the loss of the UUCG native state during MD simulations. The total amount of MD simulation data analyzed in this work is close to 1.3 ms. We identify molecular interactions, backbone conformations, and substates that are involved in the process. Then, we unravel specific ff deficiencies using diverse quantum mechanical/molecular mechanical (QM/MM) and QM calculations. Comparison between the MM and QM methods shows discrepancies in the description of the 5′-flanking phosphate moiety and both signature sugar−base interactions. Our work indicates that poor behavior of the UUCG tetraloop in simulations is a complex issue that cannot be attributed to one dominant and straightforwardly correctable factor. Instead, there is a concerted effect of multiple ff inaccuracies that are coupled and amplifying each other. We attempted to improve the simulation behavior by some carefully tailored interventions, but the results were still far from satisfactory, underlying the difficulties in development of accurate nucleic acid ff's.

Molecular mechanical (MM) force fields are commonly employed for biomolecular simulations. Despite their success, the nonpolarizable nature of contemporary additive force fields limits their performance, especially in long simulations and when strong polarization effects are present. Guanine quadruplex D(R)NA molecules have been successfully studied by MM simulations in the past. However, the G-stems are stabilized by a chain of monovalent cations that create sizable polarization effects. Indeed, simulation studies revealed several problems that have been tentatively attributed to the lack of polarization. Here, we provide a detailed comparison between quantum chemical (QM) DFT-D3 and MM potential energy surfaces of ion binding to G-stems and assess differences that may affect MM simulations. We suggest that MM describes binding of a single ion to the G-stem rather well. However, polarization effects become very significant when a second ion is present. We suggest that the MM approximation substantially limits accuracy of description of energy and dynamics of multiple ions inside the G-stems and binding of ions at the stem-loop junctions. The difference between QM and MM descriptions is also explored using symmetry-adapted perturbation theory and quantum theory of atoms in molecules analyses, which reveal a delicate balance of electrostatic and induction effects.