On differences between hydrogen bonding and improper blue-shifting hydrogen bonding.

Twenty two hydrogen-bonded and improper blue-shifting hydrogen-bonded complexes were studied by means of the HF, MP2 and B3LYP methods using the 6-31G(d,p) and 6--311 ++G(d,p) basis sets. In contrast to the standard H bonding, the origin of the improper blue-shifting H bonding is still not fully understood. Contrary to a frequently presented idea, the electric field of the proton acceptor cannot solely explain the different behavior of the H-bonded and improper blue-shifting H-bonded complexes. Compression of the hydrogen bond due to different attractive forces-dispersion or electrostatics--makes an important contribution as well. The symmetry-adapted perturbation theory (SAPT) has been utilized to decompose the total interaction energy into physically meaningful contributions. In the red-shifting complexes, the induction energy is mostly larger than the dispersion energy while, in the case of blue-shifting complexes, the situation is opposite. Dispersion as an attractive force increases the blue shift in the blue-shifting complexes as it compresses the H bond and, therefore, it increases the Pauli repulsion. On the other hand, dispersion in the red-shifting complexes increases their red shift.     

On the origin of red and blue shifts of X-H and C-H stretching vibrations in formic acid (formate ion) and proton donor complexes.

Complexes between formic acid or formate anion and various proton donors (HF, H(2)O, NH(3), and CH(4)) are studied by the MP2 and B3LYP methods with the 6-311++G(3df,3pd) basis set. Formation of a complex is characterized by electron-density transfer from electron donor to ligands. This transfer is much larger with the formate anion, for which it exceeds 0.1 e. Electron-density transfer from electron lone pairs of the electron donor is directed into sigma* antibonding orbitals of X--H bonds of the electron acceptor and leads to elongation of the bond and a red shift of the X--H stretching frequency (standard H-bonding). However, pronounced electron-density transfer from electron lone pairs of the electron donor also leads to reorganization of the electron density in the electron donor, which results in changes in geometry and vibrational frequency. These changes are largest for the C--H bonds of formic acid and formate anion, which do not participate in H-bonding. The resulting blue shift of this stretching frequency is substantial and amounts to almost 35 and 170 cm(-1), respectively.     

Contra-thermodynamic behavior in intermolecular hydrogen transfer of alkylperoxy radicals.

Quantum chemical investigation of bimolecular hydrogen transfer involving alkylperoxy radicals, a key reaction family in the free-radical oxidation of hydrocarbons, was performed to establish structure-reactivity relationships. Eight different reactions were investigated featuring four different alkane substrates (methane, ethane, propane and isobutane) and two different alkylperoxy radicals (methylperoxy and iso-propylperoxy). Including forward and reverse pairs, sixteen different activation energies and enthalpies of reaction were used to formulate structure-reactivity relationships to describe this chemistry. We observed that the enthalpy of formation of loosely bound intermediate states has a strong inverse correlation with the overall heat of reaction and that this results in unique contra-thermodynamic behavior such that more exothermic reactions have higher activation barriers. A new structure-reactivity relationship was proposed that fits the calculated data extremely well: E(A)=E(o)+alphaDeltaH(rxn) where alpha=-0.10 for DeltaH(rxn)<0, and alpha=1.10 for DeltaH(rxn)>0 and E(o)=3.05 kcal mol(-1).     

Understanding sterol-membrane interactions part I: Hartree-Fock versus DFT calculations of 13C and 1H NMR isotropic chemical shifts of sterols in solution and analysis of hydrogen-bonding effects

1H and 13C NMR chemical shifts are exquisitely sensitive probes of the local environment of the corresponding nuclei. Ultimately, direct determination of the chemical shifts of sterols in their membrane environment has the potential to reveal their molecular interactions and dynamics, in particular concerning the hydrogen-bonding partners of their OH groups. However, this strategy requires an accurate and efficient means to quantify the influence of the various interactions on chemical shielding. Herein the validity of Hartree-Fock and DFT calculations of the 13C and 1H NMR chemical shifts of cholesterol and ergosterol are compared with one another and with experimental chemical shifts measured in solution at 500 MHz. A computational strategy (definition of basis set, simpler molecular models for the sterols themselves and their molecular complexes) is proposed and compared with experimental data in solution. It is shown in particular that the effects of hydrogen bonding with various functional groups (water as a hydrogen-bond donor and acceptor, acetone) on NMR chemical shifts in CDCl3 solution can be accurately reproduced with this computational approach.     

Computational Prediction of 1H and 13C NMR Chemical Shifts for Protonated Alkylpyrroles: Electron Correlation and Not Solvation is the Salvation

Prediction of chemical shifts in organic cations is known to be a challenge. In this article we meet this challenge for α-protonated alkylpyrroles, a class of compounds not yet studied in this context, and present a combined experimental and theoretical study of the ¹³C and ¹H chemical shifts in three selected pyrroles. We have investigated the importance of the solvation model, basis set, and quantum chemical method with the goal of developing a simple computational protocol, which allows prediction of ¹³C and ¹H chemical shifts with sufficient accuracy for identifying such compounds in mixtures. We find that density functional theory with the B3LYP functional is not sufficient for reproducing all ¹³C chemical shifts, whereas already the simplest correlated wave function model, Møller–Plesset perturbation theory (MP2), leads to almost perfect agreement with the experimental data. Treatment of solvent effects generally improves the agreement with experiment to some extent and can in most cases be accomplished by a simple polarizable continuum model. The only exception is the NH proton, which requires inclusion of explicit solvent molecules in the calculation.     

Van der Waals interactions in aromatic systems: structure and energetics of dimers and trimers of pyridine.

Full geometry optimizations at the dispersion-corrected DFT-BLYP level of theory were carried out for dimers and trimers of pyridine. The DFT-D interaction energies were checked against results from single-point SCS-MP2/aug-cc-pVTZ calculations. Three stacked structures and a planar H-bonded dimer were found to be very close in energy (interaction energies in the range from -3.4 to -4.0 kcal mol(-1)). Two T-shaped geometries are higher lying, by about 1 kcal mol(-1), which is explained by the more favorable electrostatic interactions in the stacked and H-bonded arrangements. The DFT-D approach has proved to be a reliable and efficient tool to explore the conformational space of aromatic van der Waals complexes and furthermore provides interaction energies with errors of less than 10-20 % of DeltaE. Comparisons with previous results obtained by using only partially optimized model geometries strongly indicate that unconstrained optimizations are mandatory in such weakly bonded low-symmetry systems.     

Homodromic hydrogen bonds in low-energy conformations of single molecule cyclodextrins

Low-energy conformations of β-cyclodextrin under anhydrous conditions in the gas phase were investigated by DFT calculations. In these conformations, two homodromic hydrogen bond rings are formed with very short hydrogen bonds at the narrow side of the cyclodextrin ring and a second one at the wider side. These hypothetical conformations are not comparable to those conformations, which have been studied experimentally, forming inclusion complexes with small and medium-sized guest molecules, but their energy is significantly lower than the open conformations (ΔE = 10 kcal/mol).     

On the nature of the surprisingly small (red) shift in the halothane...acetone complex.

The halothane???acetone and fluoroform???acetone complexes are studied using the second‐order Møller–Plesset (MP2) method with a cc‐pVTZ basis set and the density functional theory (DFT) method with a TZVP basis set. Whereas halothane exhibits a small red shift upon complexation, fluoroform shows a pronounced blue shift. To explain this difference in behavior, we perform symmetry‐adapted perturbation theory (SAPT) and natural bond orbital (NBO) analyses. Although the composition of the total stabilization energy of each complex is different, that alone does not provide a satisfactory explanation for the difference in the spectral shifts. This difference is interpreted as a result of the interplay of the hyperconjugation and rehybridization mechanisms. The small and surprising red shift of the C? H stretching frequency of halothane, which resulted from the complexation of this species with acetone,is explained by the compensation of the two above‐mentioned mechanisms. On the other hand, the fluoroform???acetone complex exhibits a blue shift of the C? H stretching frequency upon complexation, the most likely reason for this shift being a concerted occurrence of the hyperconjugation and rehybridization mechanisms. The calculated shift of the C? H stretching vibration frequencies of halothane (+27 cm^?1) agree with the experimental value of +5 cm^?1.     

Predicting C?H/π Interactions with Nonlocal Density Functional Theory

We examine the performance of a recently developed nonlocal density functional in predicting a model noncovalent interaction, namely the weak bond between an aromatic π system and an aliphatic C? H group. The new functional is a significant improvement over traditional density functionals, providing results which compare favorably to high‐level quantum‐chemistry techniques, but at considerably lower computational cost. Interaction energies in several model C? H/π systems are in good general agreement with coupled‐cluster calculations, though equilibrium distances are consistently overpredicted when using the revPBE functional for exchange. The new functional predicts changes in energy upon addition of halogen substituents correctly.     

Double Pancake Versus Long Chalcogen–Chalcogen Bonds in Six‐Membered C,N,S‐Heterocycles

The double “pancake” bonding in the dimers of the six‐membered heterocycles 1,3‐dithia‐2,4,6‐triazine ( 4 ) and 1,3‐dithia‐2,4‐diazine ( 16 ) were investigated by means of high‐level quantum chemical calculations (B3LYP and CCSD(T)). It was found that the S–S dimers, 20 a and 27 , are not the most stable isomers, but the dimers showing short S?N ( 21 a ) and S?C ( 25 , 28 ) bonds. An investigation of the 5‐phenyl‐1,3‐dithia‐2,4,6‐triazine ( 4 b ) yields that the syn dimer with two S?S bonds (2.57 Å) is the most stable one. In this dimer, the phenyl groups are placed on top of each other. The additional dispersion energy of the phenyl rings causes a stabilization of the syn‐S–S (C_2v‐like) isomer. As a result, two weak albeit relevant single S?S bonds (2.57 Å) are predicted. These findings contradict the recently published concept of double “pancake” bonding in the dimer 4 b ₂.     

On the Mechanism of Hydrogen Activation by Frustrated Lewis Pairs

We report herein a comprehensive theoretical study of the thermodynamics and kinetics of molecular hydrogen activation by frustrated Lewis pairs (FLPs). A series of intermolecularly combined boranes (Lewis acids) and phosphines (Lewis bases), with experimentally established different reactivities towards H₂, have been subjected to DFT and (SCS‐)MP2 calculations, and analyzed in terms of their structural properties, the energetics of association of the FLPs, and the kinetics of their interactions with H₂ and hydrogenation to the ion‐pair products. The analysis included the following steps: 1) assessment of the ability/inability of the Lewis species to preorganize into FLPs with an optimum arrangement of the acid and base sites for preconditioning the reaction with H₂, 2) comprehension of the different thermodynamics of hydrogenation of the selected FLPs by comparing the Gibbs energies of the overall reactions, and 3) estimation of the mechanism of the activation of H₂ by identifying the reaction steps and the associated kinetic barriers. The results of our studies correlate well with experimental findings and have clarified the reasons for the observed different reactivities of the investigated systems, ranging from reversible or nonreversible activation to no reaction with H₂. The derived predictions could assist the future design of Lewis acid–base systems with desired properties and applicability as metal‐free hydrogenation catalysts.     

Theoretical insights into the mechanism of acetylcholinesterase-catalyzed acylation of acetylcholine

Acylation of acetylcholine (ACh) catalyzed by acetylcholinesterase (AChE) has been studied using high-level theoretical calculations on a model system that mimics the reaction center of the enzyme, and compared with uncatalyzed acylation reaction. The geometries of all the intermediates and transition states, activation energies, and solvent effects have been calculated. The calculations predict simultaneous formation of two short-strong hydrogen bonds (SSHB) in the rate-determining transition state structures [the first SSHB involves the hydrogen atom of Ser-200 (H(s)) and another involves the hydrogen atom of His-440 (H(h))]. In the intermediate states, the H-bond corresponding to H(h) involves SSHB, whereas the one corresponding to H(s) does not.     

On the Origin of the Enhanced Acidity of Chalcocyclopentadienes (Cyclopentadiene Chalcogenols) in the Gas Phase

The intrinsic acidity of chalcocyclopentadienes (CpXH; X=O, S, Se, Te) is investigated by high‐level G3B3 and G2 ab initio as well as B3LYP DFT calculations, which show that, independent of the nature of the heteroatom, all chalcocyclopentadienes are stronger acids in the gas phase than cyclopentadiene. However the acidity does not increase regularly down the group, and the acidity enhancement for Te derivatives is five times larger than for O derivatives, but only twice that of S‐containing compounds. The most favorable deprotonation process corresponds to loss of the proton attached to the heteroatom, with the sole exception of the 5‐substituted 1,3‐cyclopentadienes, for which the O and S derivatives are predicted to behave as carbon acids. No matter the nature of the heteroatom, the 1‐substituted 1,3‐cyclopentadienes are the strongest acids. The intrinsic acidity of all isomers, namely, 1‐substituted, 2‐substituted, and 5‐substituted 1,3‐cyclopentadienes, increases with increasing aromaticity of the anion formed on deprotonation, and therefore the Te compound is the strongest acid for the three series. However, the intrinsic acidity of chalcocyclopentadienes is not dictated by aromaticity, so that, in general, the most stable deprotonated species do not coincide with the most aromatic ones.     

Preparing Gold(I) for Interactions with Proton Donors: The Elusive [Au]⋅⋅⋅HO Hydrogen Bond

MP2 and DFT calculations with correlation consistent basis sets indicate that isolated linear anionic dialkylgold(I) complexes form moderately strong (ca. 10 kcal mol^?1) Au???H hydrogen bonds with single H₂O molecules as donors in the absence of sterically demanding substituents. Relativistic effects are critically important in the attraction. Such bonds are significantly weaker in neutral, strong σ‐donor N‐heterocyclic carbene (NHC) complexes (ca. 5 kcal mol^?1). The overall association (>11 kcal mol^?1), however, is strengthened by co‐operative, synergistic classical hydrogen bonding when the NHC ligands bear NH units. Further manipulation of the interaction by ligands positioned trans to the carbene, is possible.     

Highly Accurate CCSD(T) and DFT–SAPT Stabilization Energies of H‐Bonded and Stacked Structures of the Uracil Dimer

The CCSD(T) interaction energies for the H‐bonded and stacked structures of the uracil dimer are determined at the aug‐cc‐pVDZ and aug‐cc‐pVTZ levels. On the basis of these calculations we can construct the CCSD(T) interaction energies at the complete basis set (CBS) limit. The most accurate energies, based either on direct extrapolation of the CCSD(T) correlation energies obtained with the aug‐cc‐pVDZ and aug‐cc‐pVTZ basis sets or on the sum of extrapolated MP2 interaction energies (from aug‐cc‐pVTZ and aug‐cc‐pVQZ basis sets) and extrapolated ΔCCSD(T) correction terms [difference between CCSD(T) and MP2 interaction energies] differ only slightly, which demonstrates the reliability and robustness of both techniques. The latter values, which represent new standards for the H‐bonding and stacking structures of the uracil dimer, differ from the previously published data for the S22 set by a small amount. This suggests that interaction energies of the S22 set are generated with chemical accuracy. The most accurate CCSD(T)/CBS interaction energies are compared with interaction energies obtained from various computational procedures, namely the SCS–MP2 (SCS: spin‐component‐scaled), SCS(MI)–MP2 (MI: molecular interaction), MP3, dispersion‐augmented DFT (DFT–D), M06–2X, and DFT–SAPT (SAPT: symmetry‐adapted perturbation theory) methods. Among these techniques, the best results are obtained with the SCS(MI)–MP2 method. Remarkably good binding energies are also obtained with the DFT–SAPT method. Both DFT techniques tested yield similarly good interaction energies. The large magnitude of the stacking energy for the uracil dimer, compared to that of the benzene dimer, is explained by attractive electrostatic interactions present in the stacked uracil dimer. These interactions force both subsystems to approach each other and the dispersion energy benefits from a shorter intersystem separation.