Basis Set Limit Coupled Cluster Study of H-Bonded Systems and Assessment of More Approximate Methods

High-Order Quantum-Mechanical Analysis of Hydrogen Bonding in Hachimoji and Natural DNA Base Pairs

High-order quantum chemistry is applied to hydrogen-bonded natural DNA nucleobase pairs [adenine:thymine (A:T) and guanine:cytosine (G:C)] and non-natural Hachimoji nucleobase pairs [isoguanine:1-methylcytosine (B:S) and 2-aminoimidazo[1,2a][1,3,5]triazin-4(1H)-one:6-amino-5-nitropyridin-2-one (P:Z)] to see how the intermolecular interaction energies and their energetic components (electrostatics, exchange-repulsion, induction/polarization, and London dispersion interactions) vary among the base pairs. We examined the Hoogsteen (HG) geometries in addition to the traditional Watson-Crick (WC) geometries. Coupled-cluster theory through perturbative triples [CCSD(T)] extrapolated to the complete basis set (CBS) limit and high-order symmetry-adapted perturbation theory (SAPT) at the SAPT2+(3)(CCD)δMP2/aug-cc-pVTZ level are used to estimate highly accurate noncovalent interaction energies. Electrostatic interactions are the most attractive component of the interaction energies, but the sum of induction/polarization and London dispersion is nearly as large, for all base pairs and geometries considered. Interestingly, the non-natural Hachimoji base pairs interact more strongly than the corresponding natural base pairs, by -21.8 (B:S) and -0.3 (P:Z) kcal mol^-1 in the WC geometries, according to CCSD(T)/CBS. This is consistent with the H-bond distances being generally shorter in the non-natural base pairs. The natural base pairs are energetically more stabilized in their Hoogsteen geometries than in their WC geometries. The Hoogsteen geometry makes the A:T base pair slightly more stable, by -0.8 kcal mol^-1, and it greatly stabilizes the G:C⁺ base pair, by -15.3 kcal mol^-1. The G:C⁺ stabilization is mainly due to the fact that C has typically added a proton when found in Hoogsteen geometries. By contrast, Hoogsteen geometries are substantially less favorable than WC geometries for non-natural Hachimoji base pairs, by 17.3 (B:S) and 13.8 (P:Z) kcal mol^-1.

NENCI-2021. I. A large benchmark database of non-equilibrium non-covalent interactions emphasizing close intermolecular contacts

In this work, we present NENCI-2021, a benchmark database of ∼8000 Non-Equilibirum Non-Covalent Interaction energies for a large and diverse selection of intermolecular complexes of biological and chemical relevance. To meet the growing demand for large and high-quality quantum mechanical data in the chemical sciences, NENCI-2021 starts with the 101 molecular dimers in the widely used S66 and S101 databases and extends the scope of these works by (i) including 40 cation-π and anion-π complexes, a fundamentally important class of non-covalent interactions that are found throughout nature and pose a substantial challenge to theory, and (ii) systematically sampling all 141 intermolecular potential energy surfaces (PESs) by simultaneously varying the intermolecular distance and intermolecular angle in each dimer. Designed with an emphasis on close contacts, the complexes in NENCI-2021 were generated by sampling seven intermolecular distances along each PES (ranging from 0.7× to 1.1× the equilibrium separation) and nine intermolecular angles per distance (five for each ion-π complex), yielding an extensive database of 7763 benchmark intermolecular interaction energies (E_int) obtained at the coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory. The E_int values in NENCI-2021 span a total of 225.3 kcal/mol, ranging from -38.5 to +186.8 kcal/mol, with a mean (median) E_int value of -1.06 kcal/mol (-2.39 kcal/mol). In addition, a wide range of intermolecular atom-pair distances are also present in NENCI-2021, where close intermolecular contacts involving atoms that are located within the so-called van der Waals envelope are prevalent-these interactions, in particular, pose an enormous challenge for molecular modeling and are observed in many important chemical and biological systems. A detailed symmetry-adapted perturbation theory (SAPT)-based energy decomposition analysis also confirms the diverse and comprehensive nature of the intermolecular binding motifs present in NENCI-2021, which now includes a significant number of primarily induction-bound dimers (e.g., cation-π complexes). NENCI-2021 thus spans all regions of the SAPT ternary diagram, thereby warranting a new four-category classification scheme that includes complexes primarily bound by electrostatics (3499), induction (700), dispersion (1372), or mixtures thereof (2192). A critical error analysis performed on a representative set of intermolecular complexes in NENCI-2021 demonstrates that the E_int values provided herein have an average error of ±0.1 kcal/mol, even for complexes with strongly repulsive E_int values, and maximum errors of ±0.2-0.3 kcal/mol (i.e., ∼±1.0 kJ/mol) for the most challenging cases. For these reasons, we expect that NENCI-2021 will play an important role in the testing, training, and development of next-generation classical and polarizable force fields, density functional theory approximations, wavefunction theory methods, and machine learning based intra- and inter-molecular potentials.

Do Secondary Electrostatic Interactions Influence Multiple Dihydrogen Bonds? AA-DD Array on an Amine-Borane Aza-Coronand: Theoretical Studies and Synthesis

Hydrogen bond plays a key role in a wide range of inorganic, organic, as well as biological systems. The understanding on how the chemical environment can affect this kind of interaction is crucial to predict its binding strength and consequently the robustness and the dynamic properties of many supramolecular systems. In this paper a new donor-acceptor complex was synthesized and characterized by SCXRD, showing for the first time in an organic system an AA-DD pattern of a particular hydrogen interaction, called dihydrogen bond. Over 250 functionals were computationally evaluated to select the best method to reproduce the binding interaction geometry of this new pattern. Moreover, a new vector force model was used to split the contribution of primary and secondary electrostatic interactions (SEIs), in order to evaluate how the latter one can modify the binding strength of this unusual hydrogen-hydrogen interaction.

Potential energy surface and bound states of the H2O-HF complex

We present the first global five-dimensional potential energy surface for the H₂O-HF dimer, a prototypical hydrogen bonded complex. Large scale ab initio calculations were carried out using the explicitly correlated coupled cluster approach with single- and double-excitations together with non-iterative perturbative treatment of triple excitations with the augmented correlation-consistent triple zeta basis sets, in which the water and hydrogen fluoride monomers were frozen at their vibrationally averaged geometries. The ab initio data points were fitted to obtain a global potential energy surface for the complex. The equilibrium geometry of the complex corresponds to the formation of a hydrogen bond with water acting as a proton acceptor and a binding energy of D_e = 3059 cm^-1 (8.75 kcal/mol). The energies and wavefunctions of the lowest bound states of the complex were computed using a variational approach, and the dissociation energies of both ortho-H₂O-HF (D₀ = 2089.4 cm^-1 or 5.97 kcal/mol) and para-H₂O-HF (D₀ = 2079.6 cm^-1 or 5.95 kcal/mol) were obtained. The rotational constant of the complex was found to be in good agreement with the available experimental data.

Ring-Stacking Water Clusters: Morphology and Stabilities

The structures and interaction energies of water clusters with ring stacking motifs are studied by using ab initio calculations. The structures of the water clusters are constructed by stacking either single rings or multi-rings of tetramer, pentamer, and hexamer. We found that, in the single-ring-stacking motif, the most stable isomers exhibit an alternative clockwise-anticlockwise stacking pattern. We also show that four-layer single-ring-stacking isomers are not energetically favorable in comparison with those of two-layer multi-ring-stacking isomers. The relative stability of the isomers is also analyzed in terms of H-bond strength and elastic distortions of the water molecules.

Assessment of Density Functional Methods for Geometry Optimization of Bimolecular van der Waals Complexes

We explore the suitability of three popular density functionals (B97-D3, B3LYP-D3, M05-2X) for producing accurate equilibrium geometries of van der Waals (vdW) complexes with diverse binding motifs. For these functionals, optimizations using Dunning's aug-cc-pVDZ basis set best combine accuracy and a reasonable computational expense. Each DFT/aug-cc-pVDZ combination produces optimized equilibrium geometries for 21 small vdW complexes of organic molecules (up to four non-hydrogen atoms total) that agree with high-level CCSD(T)/CBS reference geometries to within ±0.1 Å for the averages of the center-of-mass displacement and the mean least root-mean-squared displacement. The DFT/aug-cc-pVDZ combinations are also able to reproduce the optimal center-of-mass displacements interpolated from CCSD(T)/CBS radial potential energy surfaces in both NBC7x and HBC6 test sets to within ±0.1 Å. We therefore conclude that each of these denisty functional methods, together with the aug-cc-pVDZ basis set, is suitable for producing equilibrium geometries of generic nonbonded complexes.

Dissociation Energy of the H2O···HF Dimer

Even though (H₂O)₂ and (HF)₂ are arguably the most thoroughly characterized prototypes for hydrogen bonding, their heterogeneous analogue H₂O···HF has received relatively little attention. Here we report that the experimental dissociation energy ( D₀) of this important paradigm for heterogeneous hydrogen bonding is too large by 2 kcal mol^-1 or 30% relative to our computed value of 6.3 kcal mol^-1. For reference, computational procedures similar to those employed here to compute D₀ (large basis set CCSD(T) computations with anharmonic corrections from second-order vibrational perturbation theory) provide results within 0.1 kcal mol^-1 of the experimental values for (H₂O)₂ and (HF)₂. Near the CCSD(T) complete basis set limit, the electronic dissociation energy for H₂O···HF is ∼4 kcal mol^-1 larger than those for (H₂O)₂ and (HF)₂ (∼9 kcal mol^-1 for the heterogeneous dimer vs ∼5 kcal mol^-1 for the homogeneous dimers). Results reported here from symmetry-adapted perturbation theory computations suggest that this large difference is primarily due to the induction contribution to the interaction energy.

Sigma-hole carbon-bonding interactions in carbon-carbon double bonds: an unnoticed contact

In this manuscript, we combine high-level ab initio calculations on some small complexes and a CSD survey to analyze the existence of unprecedented noncovalent carbon bonds in X₂C[double bond, length as m-dash]CH₂Y systems (Y = electron-rich atom or group). The methylene group is usually seen as a weak hydrogen bond donor when interacting with an electron-rich atom. However, we demonstrate that when the electron-rich atom is located equidistant from the two H atoms and along the C[double bond, length as m-dash]C bond a σ-hole noncovalent carbon-bonding interaction is established, instead of a bifurcated hydrogen bond, as derived from Atoms-in-Molecules (AIM) and Natural bond orbital (NBO) analyses. The physical nature of the interaction has been analyzed using the Symmetry Adapted Perturbation Theory (SAPT) method. The results indicate that electrostatics is very important followed by either the induction or dispersion terms in anionic and neutral complexes, respectively. In addition the CSD analysis reveals the existence of such interactions, giving reliability to our calculations, which are much more numerous for neutral than for anionic Y systems.

Some Brief Notes on Theoretical and Experimental Investigations of Intramolecular Hydrogen Bonding

A review of selected literature data related to intramolecular hydrogen bonding in ortho-hydroxyaryl Schiff bases, ortho-hydroxyaryl ketones, ortho-hydroxyaryl amides, proton sponges and ortho-hydroxyaryl Mannich bases is presented. The paper reports on the application of experimental spectroscopic measurements (IR and NMR) and quantum-mechanical calculations for investigations of the proton transfer processes, the potential energy curves, tautomeric equilibrium, aromaticity etc. Finally, the equilibrium between the intra- and inter-molecular hydrogen bonds in amides is discussed.

Accurate Ionization Potentials and Electron Affinities of Acceptor Molecules I. Reference Data at the CCSD(T) Complete Basis Set Limit

In designing organic materials for electronics applications, particularly for organic photovoltaics (OPV), the ionization potential (IP) of the donor and the electron affinity (EA) of the acceptor play key roles. This makes OPV design an appealing application for computational chemistry since IPs and EAs are readily calculable from most electronic structure methods. Unfortunately reliable, high-accuracy wave function methods, such as coupled cluster theory with single, double, and perturbative triples [CCSD(T)] in the complete basis set (CBS) limit are too expensive for routine applications to this problem for any but the smallest of systems. One solution is to calibrate approximate, less computationally expensive methods against a database of high-accuracy IP/EA values; however, to our knowledge, no such database exists for systems related to OPV design. The present work is the first of a multipart study whose overarching goal is to determine which computational methods can be used to reliably compute IPs and EAs of electron acceptors. This part introduces a database of 24 known organic electron acceptors and provides high-accuracy vertical IP and EA values expected to be within ±0.03 eV of the true non-relativistic, vertical CCSD(T)/CBS limit. Convergence of IP and EA values toward the CBS limit is studied systematically for the Hartree-Fock, MP2 correlation, and beyond-MP2 coupled cluster contributions to the focal point estimates.

From unsuccessful H2-activation with FLPs containing B(Ohfip)3 to a systematic evaluation of the Lewis acidity of 33 Lewis acids based on fluoride, chloride, hydride and methyl ion affinities

From hard to soft: The ion affinities of a large set of 33 Lewis acids towards hard and soft bases were examined with a unified isodesmic approach.

Reactivity of BrCl, Br₂, BrOCl, Br₂O, and HOBr toward dimethenamid in solutions of bromide + aqueous free chlorine

HOBr, formed via oxidation of bromide by free available chlorine (FAC), is frequently assumed to be the sole species responsible for generating brominated disinfection byproducts (DBPs). Our studies reveal that BrCl, Br(2), BrOCl, and Br(2)O can also serve as brominating agents of the herbicide dimethenamid in solutions of bromide to which FAC was added. Conditions affecting bromine speciation (pH, total free bromine concentration ([HOBr](T)), [Cl(-)], and [FAC](o)) were systematically varied, and rates of dimethenamid bromination were measured. Reaction orders in [HOBr](T) ranged from 1.09 (±0.17) to 1.67 (±0.16), reaching a maximum near the pK(a) of HOBr. This complex dependence on [HOBr](T) implicates Br(2)O as an active brominating agent. That bromination rates increased with increasing [Cl(-)], [FAC](o) (at constant [HOBr](T)), and excess bromide (where [Br(-)](o)>[FAC](o)) implicate BrCl, BrOCl, and Br(2), respectively, as brominating agents. As equilibrium constants for the formation of Br(2)O and BrOCl (aq) have not been previously reported, we have calculated these values (and their gas-phase analogues) using benchmark-quality quantum chemical methods [CCSD(T) up to CCSDTQ calculations plus solvation effects]. The results allow us to compute bromine speciation and hence second-order rate constants. Intrinsic brominating reactivity increased in the order: HOBr ≪ Br(2)O < BrOCl ≈ Br(2) < BrCl. Our results indicate that species other than HOBr can influence bromination rates under conditions typical of drinking water and wastewater chlorination.

Cross-talk between amino acid residues and flavonoid derivatives: insights into their chemical recognition

Currently, there is a general consensus that flavonoids exert their antioxidant activity through their ability to interact with a broad range of proteins, enzymes and transcription factors rather than acting as conventional hydrogen-donating antioxidants. For this, the effect of different chemical groups of the conjugated flavonoid metabolites is apparently playing a pivotal role. Yet, many questions concerning the relevant molecular mechanisms still remain open. It is therefore crucial to gain a deeper insight into the amino acid residue-flavonoid interaction. Here we show extensive theoretical thermodynamic data and structural characteristics of the interaction of chalcone, genistein, epigallocatechin gallate, and quercetin and some of its metabolites with amino acid residues. By correlating (a) the binding energies of flavonoids-amino acid residues, (b) the hydrophobicity of amino acids, and (c) the abundance of amino acid residues in the binding sites of proteins, we can conclude that flavonoids appear to be strongly bonded to only few charged hydrophilic amino acids in the protein pockets, and rather weakly bonded to the majority of amino acid residues in the binding sites. This finding strongly impacts the understanding of the chemical recognition of flavonoids and their metabolites in their interaction with proteins and would contribute to a better design of further experimental studies. Particularly, the amino acids Phe, Leu, Ile and Trp seem to play a crucial role in the dynamics of flavonoid ligands in the binding sites of proteins.

How Accurate Can a Local Coupled Cluster Approach Be in Computing the Activation Energies of Late-Transition-Metal-Catalyzed Reactions with Au, Pt, and Ir?

To improve the accuracy of local coupled cluster (LCC) methods in computing activation energies, we propose herein a new computational scheme. Its applications to various types of late-transition-metal-catalyzed reactions involving Au, Pt, and Ir indicate that the new corrective approach for LCC methods can downsize the mean unsigned deviation and maximum deviation, from the CCSD(T)/CBS reference, to about 0.3 and 0.9 kcal/mol. Using this method, we also calibrated the performance of popular density functionals, with respect to the same test set of reactions. It is concluded that the best functional is the general-purpose double hybrid functional B2GP-PLYP. Other well-performing functionals include the "kinetic" functionals M06-2X and BMK, which have a large percentage of HF exchange, and general-purpose functionals like PBE0 and wB97X. Comparatively, general-purpose functionals like PBE0 and TPSSh perform much better than the tested "kinetic" functionals for Pt-/Ir-catalyzed reactions, while the opposite is true for Au-catalyzed reactions. In contrast, wB97X performs more uniformly in these two classes of reactions. These findings hint that even within the scope of late transition metals, different types of reactions may require different types of optimal DFT methods. Empirical dispersion correction of DFT was found to have a small or no effect on the studied reactions barriers.

Thermodynamics of chemical reactions with COSMO-RS: the extreme case of charge separation or recombination

Many technically relevant chemical processes in the condensed phase involve as elementary reactive steps the formation of ions from neutral species or, as the opposite, recombination of ions. Such reactions that generate or annihilate charge defy the standard gas phase quantum chemical treatment, and also continuum solvation models are only partially able to account for the right amount of stabilization in solution. In this work, for such types of reaction, a solvation treatment involving the COSMO-RS method is assessed, which leads to improved results, i.e., errors of only around 10 kJ/mol for both protic and aprotic solvents. The examples discussed here comprise protolysis reactions and organo halide heterolysis, for both of which a comparison with reliable experimental data is possible. It is observed that for protolysis, the quality of results does not strongly depend on the quantum chemical method used for energy calculation. In contrast, in the case of heterolytic carbon-chlorine bond cleavage, clearly better results are obtained for higher correlated (coupled cluster) methods or the density functional M06-2X, which is well known for its accuracy if applied to organic chemistry. This hints at least that the right answer is obtained for the right reason and not due to a compensation of errors from gas phase thermodynamics with those from the solvation treatment. Problems encountered with certain critical solvents or upon decomposing Gibbs free energies into heats or entropies of reaction are found to relate mostly to the parameterization of the H-bonding term within COSMO-RS.

Investigating the quartz (1010)/water interface using classical and ab initio molecular dynamics

Two different terminations of the (1010) surface of quartz (α and β) interacting with water are simulated by classical (CMD) (using two different force fields) and ab initio molecular dynamics (AIMD) and compared with previously published X-ray reflectivity (XR) experiments. Radial distribution functions between hydroxyl and water show good agreement between AIMD and CMD using the ClayFF force field for both terminations. The Lopes et al. (Lopes, P. E. M.; Murashov, V.; Tazi, M.; Demchuk, E.; MacKerell, A. D. J. Phys. Chem. B2006, 110, 2782-2792) force field (LFF), however, underestimates the extent of hydroxyl-water hydrogen bonding. The β termination is found to contain hydroxyl-hydroxyl hydrogen bonds; the quartz surface hydroxyl hydrogens and oxygens that hydrogen bond with each other exhibit greatly reduced hydrogen bonding to water. Conversely, the hydroxyl hydrogen and oxygens that are not hydrogen bonded to other surface hydroxyls but are connected to those that are show a considerable amount of hydrogen bonding to water. The electron density distribution of an annealed surface of quartz (1010) obtained by XR is in qualitative agreement with electron densities calculated by CMD and AIMD. In all simulation methods, the interfacial water peak appears farther from the surface than observed by XR. Agreement among AIMD, LFF, and XR is observed for the relaxation of the near-surface atoms; however, ClayFF shows a larger discrepancy. Overall, results show that for both terminations of (1010), LFF treats the near-surface structure more accurately whereas ClayFF treats the interfacial water structure more accurately. It is shown that the number of hydroxyl and water hydrogen bonds to the bridging Si-O-Si oxygens connecting the surface silica groups to the rest of the crystal is much greater for the α than the β termination. It is suggested that this may play a role in the greater resistance to dissolution of the β termination than that of the α termination.

The effect of polarization on multiple hydrogen-bond formation in models of self-assembling materials

We report density functional theory calculations at the B3LYP/D95(d,p) level on several different cyclic H-bonding dimers, where the monomers of each are connected by a pair of N-H···O=C H-bonding interactions, and the H-bonding donors and acceptors on each monomer are separated by polarizable spacers. Depending on the structures, the individual H-bonds vary in strength (enthalpy) by over a factor of four, from 2.41 to 10.99 kcal/mol. We attribute most of the variation in interaction energies to differences in the extent of polarization due to each of the H-bonds, which can either combine constructively or destructively. The dipole-dipole interactions between the pair of H-bonds also contribute somewhat to the relative stabilities. The relevance of these results to the design of self-assembling materials is discussed.

A direct comparison of azide and nitrile vibrational probes

The synthesis of 2'-azido-5-cyano-2'-deoxyuridine, N(3)CNdU (1), from trityl-protected 2'-amino-2'-deoxyuridine was accomplished in four steps with a 12.5% overall yield. The IR absorption positions and profiles of the azide and nitrile group of N(3)CNdU were investigated in 14 different solvents and water/DMSO solvent mixtures. The azide probe was superior to the nitrile probe in terms of its extinction coefficient, which is 2-4 times larger. However, the nitrile IR absorbance profile is generally less complicated by accidental Fermi resonance. The IR frequencies of both probes undergo a substantial red shift upon going from water to aprotic solvents such as THF or DMSO. DFT calculations supported the hypothesis that the molecular origin of the higher observed frequency in water is primarily due to hydrogen bonds between the probes and water molecules.