Elucidating Solvent Effects on Strong Intramolecular Hydrogen Bond: DFT-MD Study of Dibenzoylmethane in Methanol Solution

The dynamic aspect of solvation plays a crucial role in determining properties of strong intramolecular hydrogen bonds since solvent fluctuations modify instantaneous hydrogen-bonded proton transfer barriers. Previous studies pointed out that solvent-solute interactions in the first solvation shell govern the position of the proton but the ability of the electric field due to other solvent molecules to localize the proton remains an important issue. In this work, we examine the structure of the O−H⋅⋅⋅O intramolecular hydrogen bond of dibenzoylmethane in methanol solution by employing density functional theory-based molecular dynamics and quantum chemical calculations. Our computations showed that homogeneous electric fields with intensities corresponding to those found in polar solvents are able to considerably alter the proton transfer barrier height in the gas phase. In methanol solution, the proton position is correlated with the difference in electrostatic potentials on the oxygen atoms of dibenzoylmethane even when dibenzoylmethane-methanol hydrogen bonding is lacking. On a timescale of our simulation, the hydrogen bonding and solvent electrostatics tend to localize the proton on different oxygen atoms. These findings provide an insight into the importance of the solvent electric field on the structure of a strong intramolecular hydrogen bond.     

Prereactive Complexes of Dihalogens XY with Lewis Bases B in the Gas Phase: A Systematic Case for the Halogen Analogue B⋅⋅⋅XY of the Hydrogen Bond B⋅⋅⋅HX

Similar and yet also notably different are the B⋅⋅⋅XY and B⋅⋅⋅HX complexes in the gas phase, where B is a simple Lewis base, XY is a homo- or heterodihalogen molecule, and HX is a hydrogen halide. This is demonstrated, for example, by the structures of oxirane⋅⋅⋅ClF and oxirane⋅⋅⋅HCl (see picture). Both bonds are dominated by simple electrostatic interactions, but differ in terms of their propensity for nonlinearity.     

Nature and Strength of Weak O⋅⋅⋅O Interactions in Nitryl Halide Dimers

The use of real space functions and molecular graphs has pushed some chemists to wonder: Are interactions between negatively charged oxygen atoms possible? In this contribution we analyze whether there is a real interaction between oxygen atoms in nitryl halide dimers (XNO₂)₂ (X=F, Cl, Br and I) and in tetranitromethane and derivatives. Based on ab-initio and density functional theories (DFT) methods, we show these complexes are weakly stabilized. Energy decomposition analyses based on local molecular orbitals (LMOEDA) and interacting quantum atoms (IQA) reveal both dispersion and exchange play a crucial role in the stabilization of these complexes. Electron charge density and IQA analyses indicate that the oxygen atoms are connected by privileged exchange channels. In addition, electrostatic interactions between O and N atoms are also vital for the stabilization of the complexes. Finally, a reasonable explanation is given for the dynamic behavior of nitryl groups in tetranitromethane and derivatives.     

Probing the C?H⋅⋅⋅π Weak Hydrogen Bond in Anesthetic Binding: The Sevoflurane–Benzene Cluster

Cooperativity between weak hydrogen bonds can be revealed in molecular clusters isolated in the gas phase. Here we examine the structure, internal dynamics, and origin of the weak intermolecular forces between sevoflurane and a benzene molecule, using multi‐isotopic broadband rotational spectra. This heterodimer is held together by a primary C H⋅⋅⋅π hydrogen bond, assisted by multiple weak C H⋅⋅⋅F interactions. The multiple nonbonding forces hinder the internal rotation of benzene around the isopropyl C H bond in sevoflurane, producing detectable quantum tunneling effects in the rotational spectrum.     

Crystal and molecular structures of atropisomeric N‐aryl‐1,2,3,4‐tetrahydro‐3,3‐dimethyl‐2,4‐quinolinediones

The crystal structures of N‐aryl‐1,2,3,4‐tetrahydro‐3,3‐dimethyl‐2,4‐quinolinediones bearing methoxy‐ ( 1 ), methyl‐ ( 2 ), and chloro‐ ( 3 ) substituents in 2′‐position of the phenyl ring have been determined by X‐ray crystal structure analysis. The heterocyclic ring in 1–3 adopts an envelope conformation, with the smallest ring puckering in the ortho‐chloro derivative 3 . The N‐aryl ring is almost perpendicular with respect to the quinoline‐2,4‐dione ring. The corresponding dihedral angle values are 83.2(1)°, 80.0(9)°, and 83.4(2)° in 1, 2 and 3 , respectively. The hydrogen bond of C H⋅⋅⋅O type joins the molecules of the ortho‐methoxy derivative 1 into dimers. The supramolecular structure also contains two C H⋅⋅⋅π interactions that link the hydrogen‐bonded dimers into sheets. In ortho‐methyl derivative 2 , one C H⋅⋅⋅π interaction generates infinite chains, whereas two C H⋅⋅⋅O hydrogen bonds and three C H⋅⋅⋅π interactions in the ortho‐chloro derivative 3 form three‐dimensional framework. © 2008 Wiley Periodicals, Inc. Heteroatom Chem 19:325–331, 2008; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/hc.20436     

Density functional theory and topological analysis on the hydrogen bonding interactions in cysteine‐thymine complexes

Hydrogen bonding interactions between amino acids and nucleic acid bases constitute the most important interactions responsible for the specificity of protein binding. In this study, complexes formed by hydrogen bonding interactions between cysteine and thymine have been studied by density functional theory. The relevant geometries, energies, and IR characteristics of hydrogen bonds (H‐bonds) have been systematically investigated. The quantum theory of atoms in molecule and natural bond orbital analysis have also been applied to understand the nature of the hydrogen bonding interactions in complexes. More than 10 kinds of H‐bonds including intra‐ and intermolecular H‐bonds have been found in complexes. Most of intermolecular H‐bonds involve O (or N) atom as H‐acceptor, whereas the H‐bonds involving C or S atom usually are weaker than other ones. Both the strength of H‐bonds and the structural deformation are responsible for the stability of complexes. Because of the serious deformation, the complex involving the strongest H‐bond is not the most stable structures. Relationships between H‐bond length (ΔR_X‐H), frequency shifts (Δv), and the electron density (ρ_b) and its Laplace (?²ρ_b) at bond critical points have also been investigated. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Arrangement and nature of intermolecular hydrogen bonding in complex biomolecular systems: modeling the vitamin C---L-alanine interaction

Vitamin C is known as an essential dietary supplement and implicated in diverse biological processes. We present here a theoretical study on the nature of hydrogen bonding of vitamin C in biological systems. For this reason, the complexes of vitamin C (VC) with neutral and zwitterionic L-alanine (as the simplest chiral amino acid) were studied at the MP2/6-311++G(d,p) level of theory. In the gas phase, neutral L-alanine leads to more stable complexes than the zwitterionic forms while the reverse is true in the aqueous phase. The complexes are formed via two hydrogen bond interactions, which result in a ring-like hydrogen-bonded networks. The nature of H-bonds was characterized in terms of natural bond orbital and quantum theory of atoms in molecule analyses (QTAIM). The H-bonds in the studied complexes were electrostatic in nature; however, in the case of shorter and directional H-bonds and ionic interactions, contributions of covalent character were also non-negligible. Natural energy decomposition analysis of hydrogen-bonded complexes reveals that the charge transfer and electrical components are the largest contributors for the interaction energies of complexes. Natural resonance theory analysis suggests higher resonance weight for charge-assisted interactions of vitamin C---alanine (zwitterionic) complexes, where the total interaction energy is considerably higher than that of neutral alanine.     

Conformational Equilibria of 2-Methoxypyridine⋅⋅⋅CO2: Cooperative and Competitive Tetrel and Weak Hydrogen Bonds

The rotational spectrum of 2-methoxypyridine⋅⋅⋅CO₂ was recorded and analysed employing a cavity-based Fourier transform microwave spectrometer, complemented with quantum chemical calculations which predicted three possible isomers within energies less than 1000 cm⁻¹. The two most stable isomers were observed in the pulsed jet, which are stabilized by a network of C⋅⋅⋅N/O tetrel and C−H⋅⋅⋅O weak hydrogen bonds. The relative population ratio of the two detected isomers was estimated to be N_I/N_II≈2.5. The competition and cooperation of the present non-covalent interactions in both isomers are discussed within the framework of Bader's quantum theory of atoms in molecules and Johnson's non-covalent interaction analyses. The study shows, that when looking for CO₂ adsorbents, one might prefer candidates with multiple interactions in one site over candidates with few but strong interactions.     

Charged versus Neutral Hydrogen‐Bonded Complexes: Is There a Difference in the Nature of the Hydrogen Bonds?

A theoretical study on some carboxylic acid dimers formed by positively or negatively charged molecules has been carried out by using DFT methods. The resulting dimers possess either a charge of +2 or ?2. In addition, the corresponding neutral complexes have also been considered. The electron density distribution described by the atoms in molecules and the natural bond orbital methods, as well as the electric field maps of the systems, have been analyzed and compared without finding significant differences between the neutral and ionic complexes. The interaction energy along the dissociation path of the charged dimers shows both a local minimum and a local maximum, defining a stability region between them. When this energetic profile is recalculated by removing the repulsion between the charged groups, it resembles to those of the neutral molecules. Hence, the characteristics of the charged dimers are similar to those of the neutral ones: the addition of a repulsion term for the charged groups permits to retrieve the energetic profiles dependence with the distance in the charged system. The interacting quantum atom (IQA) method has been used to calculate the interaction energy terms, including the classic Coulombic term between the whole molecules and the corresponding of the carboxylic acid groups. The IQA results show repulsive electrostatic interactions when the whole molecules are considered in the ionic complexes, but attractive ones between the carboxylic groups in both neutral and ionic complexes.     

Non-Classical Synthons: Supramolecular Recognition by S⋅⋅⋅O Chalcogen Bonding in Molecular Complexes of Riluzole

Classical examples of supramolecular recognition units or synthons are the ones formed by hydrogen bonds. Here, we report the ubiquity of a S⋅⋅⋅O chalcogen bonded synthon observed in a series of supramolecular complexes of the amyotrophic lateral sclerosis drug riluzole. Although the potential of higher chalcogens such as Se and Te to form robust and directional chalcogen bonded motifs is known, intermolecular sulfur chalcogen bonding is considered to be weak owing to the lower polarizability of S atoms. Here, the robustness and electronic nature of a S⋅⋅⋅O chalcogen bonding non-classical synthon, and the origin of its exceptional directionality have been explored. Bond orders of the drug–coformer chalcogen bonding are found to be as high as one third of a single bond, and they are largely ionic in nature. The contribution of the S⋅⋅⋅O chalcogen bonded motifs to the lattice energies of a series of crystals from the Cambridge Structural Database has been analyzed, showing they can be indeed significant, especially in molecules devoid of strong hydrogen bond donor groups.     

pi-H...O hydrogen bonds: multicenter covalent pi-H interaction acts as the proton-donating system

The MP2 method and the Pople-style basis sets 6-311++G(d,p), 6-311++G(2df,2pd), and 6-311++G(3df,3pd) were used to perform calculations on H3O+...C2H2 and C2H3+...C2H2 complexes and related species. Hydrogen bonds existing for the analyzed complexes were investigated as well as related pi-H...O --> pi...H-O and pi-H...pi --> pi...H-pi proton-transfer processes. For some of the complexes analyzed the multicenter pi-H interaction possessing the properties of a covalent bond acts as a proton donor; more generally it is classified as the Lewis acid. The quantum theory of "atoms in molecules" (QTAIM) was also applied to deepen the nature of these interactions in terms of characteristics of bond critical points. The pi-H...O, O-H...pi, and pi-H...pi interactions analyzed here may be classified as hydrogen bonds since their characteristics are the same as or at least similar to those of typical hydrogen bonds. H...pi interactions are common in crystal structures of organic and organometallic compounds. The analyses performed here show a continuum of such interactions since there are H...pi contacts possessing the characteristics of weak intermolecular interactions on the one hand and pi-H multicenter covalent bonds on the other. Ab initio and QTAIM results support the latter statements.     

Ab initio investigation of the complexes between bromobenzene and several electron donors: some insights into the magnitude and nature of halogen bonding interactions

Halogen bonding, a specific intermolecular noncovalent interaction, plays crucial roles in fields as diverse as molecular recognition, crystal engineering, and biological systems. This paper presents an ab initio investigation of a series of dimeric complexes formed between bromobenzene and several electron donors. Such small model systems are selected to mimic halogen bonding interactions found within crystal structures as well as within biological molecules. In all cases, the intermolecular distances are shown to be equal to or below sums of van der Waals radii of the atoms involved. Halogen bonding energies, calculated at the MP2/aug-cc-pVDZ level, span over a wide range, from -1.52 to -15.53 kcal/mol. The interactions become comparable to, or even prevail over, classical hydrogen bonding. For charge-assisted halogen bonds, calculations have shown that the strength decreases in the order OH- > F- > HCO2- > Cl- > Br-, while for neutral systems, their relative strengths attenuate in the order H2CS > H2CO > NH3 > H2S > H2O. These results agree with those of the quantum theory of atoms in molecules (QTAIM) since bond critical points (BCPs) are identified for these halogen bonds. The QTAIM analysis also suggests that strong halogen bonds are more covalent in nature, while weak ones are mostly electrostatic interactions. The electron densities at the BCPs are recommended as a good measure of the halogen bond strength. Finally, natural bond orbital (NBO) analysis has been applied to gain more insights into the origin of halogen bonding interactions.     

Theoretical Study on the Hydrogen Bonding Interactions in Complexes of 5‐Hydroxytryptamine with Water

The energies, geometries and harmonic vibrational frequencies of 1:1 5‐hydroxytryptamine‐water (5‐HT‐H₂O) complexes are studied at the MP2/6‐311++G(d,p) level. Natural bond orbital (NBO), quantum theory of atoms in molecules (QTAIM) analyses and the localized molecular orbital energy decomposition analysis (LMO‐EDA) were performed to explore the nature of the hydrogen‐bonding interactions in these complexes. Various types of hydrogen bonds (H‐bonds) are formed in these 5‐HT‐H₂O complexes. The intermolecular C4H5^5‐HT···O^w H‐bond in HTW3 is strengthened due to the cooperativity, whereas no such cooperativity is found in the other 5‐HT‐H₂O complexes. H‐bond in which nitrogen atom of amino in 5‐HT acted as proton donors was stronger than other H‐bonds. Our researches show that the hydrogen bonding interaction plays a vital role on the relative stabilities of 5‐HT‐H₂O complexes.     

Structures and Chemical Bonding in Antimony(III) Bromide Complexes with Pyridine

Weakly or “partially” bonded molecules are an important link between the chemical and van der Waals interactions. Molecular structures of six new SbBr₃-Py complexes in the solid state have been determined by single-crystal X-ray diffraction analysis. In all complexes all Sb atoms adopt a pseudo-octahedral coordination geometry which is completed by additional Sb⋅⋅⋅Br contacts shorter than the sum of the van der Waals radii, with Br−Sb⋅⋅⋅Br angles close to 180°. To reveal the nature of Sb–Br and Sb–N interactions, the DFT calculations were performed followed by the analysis of the electrostatic potentials, the orbital interactions and the topological analysis. Based on Natural Bond Orbital (NBO) analysis, the Sb–Br interactions range from the covalent bonds to the pnictogen bonds. A simple structural parameter, non-covalence criterion (NCC) is defined as a ratio of the atom-atom distance to the linear combination of sums of covalent and van der Waals radii. NCC correlates with E(2) values for Sb−N, Sb−Cl and Sb−Br bonds, and appears to be useful criterion for a preliminary evaluation of the bonding situation.     

An Experimental Exploration of C−H⋅⋅⋅X Hydrogen Bond in [CHCl3−X(CH3)2] Complexes (X=O,S, and Se)

Among the conglomeration of hydrogen bond donors, the C−H group is prevalent in chemistry and biology. In the present work, CHCl₃ has been selected as the hydrogen bond donor and are X(CH₃)₂ are the hydrogen bond acceptors. Formation of C−H⋅⋅⋅X hydrogen bond under the matrix isolation condition is confirmed by the observation of red-shift in the C−H stretching frequency of CHCl₃ and comparison with the simulated spectra. Stabilisation energy of all the three complexes is almost equal although the observed red-shift for the C−H⋅⋅⋅O complex is less compared to the C−H⋅⋅⋅S/Se complexes. The nature and origin of the hydrogen bond have been delineated using Natural Bond Orbital, Atoms in Molecules, Non-Covalent Interaction analyses, and Energy Decomposition Analysis. Charge transfer is found to be proportional to the observed red-shift. This work provides the first impression of C−H⋅⋅⋅Se hydrogen bond and its comparison with C−H⋅⋅⋅O/S hydrogen bond interaction under experimental condition.     

Spectroscopic and Crystallographic Characterization of the R3N+−C−H⋅⋅⋅X Interaction

As appreciation for nonclassical hydrogen bonds has progressively increased, so have efforts to characterize these interesting interactions. Whereas several kinds of C−H hydrogen bonds have been well-studied, much less is known about the R₃N⁺−C−H⋅⋅⋅X variety. Herein, we present crystallographic and spectroscopic evidence for the existence of these interactions, with special relevance to Selectfluor chemistry. Of particular note is the propensity for Lewis bases to engage in nonclassical hydrogen bonding over halogen bonding with the electrophilic F atom of Selectfluor. Further, the first examples of ¹H NMR experiments detailing R₃N⁺−C−H⋅⋅⋅X (X=O, N) hydrogen bonds are described.     

Hydrogen Bonding between Metal‐Ion Complexes and Noncoordinated Water: Electrostatic Potentials and Interaction Energies

The hydrogen bonding of noncoordinated water molecules to each other and to water molecules that are coordinated to metal‐ion complexes has been investigated by means of a search of the Cambridge Structural Database (CSD) and through quantum chemical calculations. Tetrahedral and octahedral complexes that were both charged and neutral were studied. A general conclusion is that hydrogen bonds between noncoordinated water and coordinated water are much stronger than those between noncoordinated waters, whereas hydrogen bonds of water molecule in tetrahedral complexes are stronger than in octahedral complexes. We examined the possibility of correlating the computed interaction energies with the most positive electrostatic potentials on the interacting hydrogen atoms prior to interaction and obtained very good correlation. This study illustrates the fact that electrostatic potentials computed for ground‐state molecules, prior to interaction, can provide considerable insight into the interactions.