Strong Short‐Range Cooperativity in Hydrogen‐Bond Chains

Chains of hydrogen bonds such as those found in water and proteins are often presumed to be more stable than the sum of the individual H bonds. However, the energetics of cooperativity are complicated by solvent effects and the dynamics of intermolecular interactions, meaning that information on cooperativity typically is derived from theory or indirect structural data. Herein, we present direct measurements of energetic cooperativity in an experimental system in which the geometry and the number of H bonds in a chain were systematically controlled. Strikingly, we found that adding a second H‐bond donor to form a chain can almost double the strength of the terminal H bond, while further extensions have little effect. The experimental observations add weight to computations which have suggested that strong, but short‐range cooperative effects may occur in H‐bond chains.     

Intramolecular hydrogen bonding and cooperative interactions in carbohydrates via the molecular tailoring approach

In spite of many theoretical and experimental attempts for understanding intramolecular hydrogen bonding (H-bonding) in carbohydrates, a direct quantification of individual intramolecular H-bond energies and the cooperativity among the H-bonded networks has not been reported in the literature. The present work attempts, for the first time, a direct estimation of individual intramolecular O-H...O interaction energies in sugar molecules using the recently developed molecular tailoring approach (MTA). The estimated H-bond energies are in the range of 1.2-4.1 kcal mol(-1). It is seen that the OH...O equatorial-equatorial interaction energies lie between 1.8 and 2.5 kcal mol(-1), with axial-equatorial ones being stronger (2.0-3.5 kcal mol(-1)). The strongest bonds are nonvicinal axial-axial H-bonds (3.0-4.1 kcal mol(-1)). This trend in H-bond energies is in agreement with the earlier reports based on the water-water H-bond angle, solvent-accessible surface area (SASA), and (1)H NMR analysis. The contribution to the H-bond energy from the cooperativity is also estimated using MTA. This contribution is seen to be typically between 0.1 and 0.6 kcal mol(-1) when H-bonds are a part of a relatively weak equatorial-equatorial H-bond network and is much higher (0.5-1.1 kcal mol(-1)) when H-bonds participate in an axial-axial H-bond network.     

Characterizing the cooperativity in H-bonded amino structures

Density functional theory calculations were used to examine the effect of H-bond cooperativity on the magnitude of the NMR chemical shifts and spin-spin coupling constants in a C4h-symmetric G-quartet and in structures consisting of six cyanamide monomers. These included two ring structures (a planar C6h-symmetric structure and a nonplanar S6-symmetric structure) and two linear chain structures (a fully optimized planar Cs-symmetric chain and a planar chain structure where all intra- and intermolecular parameters were constrained to be identical). The NMR parameters were computed for the G-quartet and cyanamide structures, as well as for shorter fragments derived from these assemblies without reoptimization. In the ring structures and the chain with identical monomers, the intra- and intermolecular geometries of the cyanamides were identical, thereby allowing the study of cooperative effects in the absence of geometry changes. The magnitude of the |1JNH| coupling, 1H and 15N chemical shifts of the H-bonding amino N-H group, and the |h2JNN| H-bond coupling increased, whereas the size of the |1JNH| coupling of the non-H-bonded amino N-H bonds of the first amino group in the chain, which are roughly perpendicular to the H-bonding network, decreased in magnitude when H-bonding monomers were progressively added to extending ring or chain structures. These effects are attributed to electron redistribution induced by the presence of the nearby H-bonding guanine or cyanamide molecules.     

Theoretical and experimental study of the adsorption of neutral glycine on silica from the gas phase.

The adsorption of neutral glycine onto amorphous silica was investigated both theoretically and experimentally. DFT calculations were performed at the BLYP-631++G** level using a cluster approach. Several possible configurations involving the formation of H bonds between glycine and one, two, or three silanol groups (SiOH) were considered. The most favorable bonding of glycine with one silanol group (45 kJ mol(-1)) occurs through the COOH moiety, thus forming a cycle in which the CO group is an H-bond acceptor whereas the acidic OH group is an H-bond donor. With two or three silanol groups, additional H bonds are formed between the amine moiety and the silanol groups, which leads to an increased adsorption energy (70 and 80 kJ mol(-1) for two and three silanol groups, respectively). Calculated nu(CO), delta(HNH), and delta(HCH) values are sensitive to the adsorption mode. A bathochromic shift of nu(CO) as compared to the nu(CO) of free glycine (calculated in the 1755-1790 cm(-1) range) is found for glycine in interaction with silanol(s). The more H bonds are formed between the COOH moiety and silanol groups, the higher the bathochromic shift. For delta(HNH), no shift is found for glycine adsorbed on one and two silanol groups (where the amine is either not bound or an H-bond donor), whereas a bathochromic shift is calculated with three silanols when the amine moiety is an H-bond acceptor. Experimental FTIR spectra performed at room temperature for glycine adsorbed at 160 degrees C on Aerosil amorphous silica exhibit bands at 1371, 1423, 1630, and 1699 cm(-1). The experimental/calculated frequencies have their best correspondence for glycine adsorbed on two silanol groups. It is important to note that the forms giving the best correspondence to experimental frequencies are the most stable ones.     

Spectroscopic and theoretical evidence for the cooperativity between red-shift hydrogen bond and blue-shift hydrogen bond in DMSO aqueous solutions

The cooperativity between red-shifted hydrogen bond and blue-shifting hydrogen bond in dimethyl sulfoxide aqueous solutions was studied by methods of quantum chemical calculations and infrared spectroscopy. The water molecule plays a different role in two types of hydrogen bonds: proton-donor in red-shifted hydrogen bond and proton-acceptor in blue-shifting hydrogen bond. The cooperativity is not prominent if the ring structure is formed through the OHcdots, three dots, centeredOS H-bond and CHcdots, three dots, centeredO(w) H-bond. However, if the methyl groups in the above ring structure participate in second CHcdots, three dots, centeredO(w) H-bond, the cooperativity is increased. The second CHcdots, three dots, centeredO(w) H-bond enhances OHcdots, three dots, centeredOS H-bond and weakens the first CHcdots, three dots, centeredO(w) H-bond.     

Cooperativity between OH...O and CH...O hydrogen bonds involving dimethyl sulfoxide-H2O-H2O complex

The cooperativity between the O-H...O and C-H...O hydrogen bonds has been studied by quantum chemical calculations at the MP2/6-311++G(d,p) level in gaseous phase and at the B3LYP/6-311++G(d,p) level in solution. The interaction energies of the O-H...O and C-H...O H-bonds are increased by 53 and 58%, respectively, demonstrating that there is a large cooperativity. Analysis of hydrogen-bonding lengths, OH bond lengths, and OH stretching frequencies also supports such a conclusion. By NBO analysis, it is found that orbital interaction plays a great role in enhancing their cooperativity. The strength increase of the C-H...O H-bond is larger than that of the O-H...O H-bond due to the cooperativity. The solvent has a weakening effect on the cooperativity.     

Cooperativity of conventional and unconventional hydrogen bonds involving imidazole

Correlated ab initio calculations are used to investigate the cooperativity of H‐bonds between imidazole and a pair of water molecules. H‐bonds comprise not only the conventional NH … O and OH … N types, but also CH … O and OH … φ (wherein a proton is donated to the delocalized π cloud lying above the aromatic ring). Conventional and OH … φ H‐bonds obey the normal principles of cooperativity, wherein these bonds are strengthened when a central molecule serves simultaneously as both proton donor and acceptor. In contrast, CH … O bonds do not appear to be amenable to such positive cooperativity. When placed in a polarizable medium, all H‐bonds weaken as the dielectric constant of the solvent grows. The qualitative aspects of the cooperativity are not affected by the medium, although some weakening is observed. Calculations also consider the effects of cooperativity on other aspects of the complexes, including the intermolecular distance, the effect on the covalent X‐H bond length, and IR and NMR spectral data. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Nature of one-dimensional short hydrogen bonding: bond distances, bond energies, and solvent effects

On the basis of recently synthesized calix[4]hydroquinone (CHQ) nanotubes which were self-assembled with infinitely long one-dimensional (1-D) short hydrogen bonds (SHB), we have investigated the nature of 1-D SHB using first-principles calculations for all the systems including the solvent water. The H-bonds relay (i.e., contiguous H-bonds) effect in CHQs shortens the H...O bond distances significantly (by more than 0.2 A) and increases the bond dissociation energy to a large extent (by more than approximately 4 kcal/mol) due to the highly enhanced polarization effect along the H-bond relay chain. The H-bonds relay effect shows a large increase in the chemical shift associated with the SHB. The average binding energies for the infinite 1-D H-bond arrays of dioles and dions increase by approximately 4 and approximately 9 kcal/mol per H-bond, respectively. The solvent effect (due to nonbridging water molecules) has been studied by explicitly adding water molecules in the CHQ tube crystals. This effect is found to be small with slight weakening of the SHB strength; the H...O bond distance increases only by 0.02 A, and the average binding energy decreases by approximately 1 kcal/mol per H-bond. All these results based on the first-principles calculations are the first detailed analysis of energy gain by SHB and energy loss by solvent effect, based on a partitioning scheme of the interaction energy components. These reliable results elucidate not only the self-assembly phenomena based on the H-bond relay but also the solvent effect on the SHB strength.     

Structure of liquid formic acid investigated by first principle and classical molecular dynamics simulations

The structure of liquid formic acid has been investigated by Car-Parrinello and classical molecular dynamics simulations, focusing on the characterization of the H-bond network and on the mutual arrangement of pairs of bonded molecules. In agreement with previous computational studies, two levels of H-bonded structures have been found. Small clusters, characterized by O-H...O bonds, are held together by weak C-H...O bonds to form large branched structures. From the ab initio simulation we infer the importance of cyclic H-bond dimer configurations, typical of the gas phase. Most of these dimer structures are however found to be embedded into H-bonded chains. When only O-H...O bonds are taken into account, linear H-bond chains are detected as basic structures of the liquid. More branched structures occur when C-H...O bonds are also considered. Regarding the arrangement of molecular pairs, we observed that O-H...O bonds favor the occurrence of configurations with parallel molecular planes, whereas no preferential orientation is observed for molecules forming C-H...O bonds.     

Properties of the H-bond network for two-dimensional lattice water model

A microscopic Hamiltonian of the hydrogen-bond network in two-dimensional lattice water is proposed, which describes the formation and disruption of the H bonds, their bending, and which satisfies the Bernal-Fowler rules [J. D. Bernal and R. H. Fowler, J. Chem. Phys. 1, 515 (1933)]. The thermodynamic properties of the H-bond network are studied using the method of many-particle irreducible distribution functions, which is a generalization of the Kikuchi cluster approach [R. Kikuchi, Phys. Rev. 81, 988 (1951)] and the Bethe-Peierls quasiactivities method [H. A. Bethe, Prog. R. Soc. A 150, 552 (1935)]. The temperature dependencies of the average number of H bonds per molecules, the contribution of the H bonds into the heat capacity of the system, and the parameters describing the correlations between the states of molecules on the neighboring sites are investigated. It is shown that depending on the magnitude of the interaction between the H bonds in the H-bond subsystem either smooth or sharp first-order phase transition can occur. The role of different factors in the formation of the properties of the H-bond network is discussed.     

Interaction of glycine with isolated hydroxyl groups at the silica surface: first principles B3LYP periodic simulation

The adsorption of a glycine molecule on a model silica surface terminated by an isolated hydroxyl group has been studied ab initio using a double-zeta polarized Gaussian basis set, the hybrid B3LYP functional, and a full periodic treatment of the silica surface/glycine system. The hydroxylated silica surface has been simulated using either a 2D slab or a single polymer strand cut out from the (001) surface of an all-silica edingtonite. A number of B3LYP-optimized structures have been found by docking glycine on the silica surface exploiting all possible hydrogen bond patterns. Whereas glycine is generally adsorbed in its neutral form, two structures show glycine adsorbed as a zwitterion, the surface playing the role of a "solid solvent" whereas intrastrand hydrogen bond cooperativity stabilizes the zwitterions. The adsorbed zwitterionic structures are no longer formed at a lower glycine coverage as simulated by enlarging the unit cell so as to break intrastrand hydrogen bonds, showing the importance of H-bond cooperativity in stabilizing the zwitterionic forms. Each structure has been characterized by computing its harmonic vibrational spectrum at the Gamma point, which also allowed us to calculate the free energy of adsorption. The experimental infrared features of chemical-vapor-deposited glycine on a silica surface are in agreement with those computed for glycine adsorbed in its neutral form and engaging three hydrogen bonds with the surface silanols, two of them involving the C=O bond and one originating from the glycine OH group. The NH(2) group plays only a minor role as a weak hydrogen bond donor.     

Disruption of hydrogen bond structure of water near charged electrode surfaces

The understanding of the hydrogen (H) bonded structure of water near charged surfaces is highly relevant in the context of several important areas of research, including electrochemistry, biochemistry, and geology. Past simulation studies have not yielded conclusive answers; while some suggest breakage of H bonds near a charged surface, others argue that H-bonding interactions can stabilize the structure of surface water even in the presence of high electric (E) fields. Recent experiments, on the other hand, suggest a partial breakdown of H-bond structure near a charged electrode. In all these studies, however, the conclusions regarding H bonding were drawn based on the density profile of hydrogen/oxygen atoms near the interface. In the present paper, we investigate this problem using a new theory that explicitly accounts for the influence of E field on the H-bond network of water near the solid-liquid interface. We find that the average number of H bonds per molecule in bulk increases from approximately 3.8 at E<10(5) V/m to approximately 3.95 at E=2x10(9) V/m (suggesting enhancement in H-bond network), while that near the electrode surface decreases from approximately 2.8 to a saturation value of approximately 2.0 (suggesting weakening of H-bond network).     

Hydrogen bonding in phosphonate cavitands: Investigation of host-guest complexes with ammonium salts

The H-bonding in alkylammonium complexes of phosphonate cavitands were studied by mass spectrometric methods and theoretical calculations. The alkylammonium ions included primary, secondary, and tertiary methyl- and ethylammonium ions. Their complexation with mono-, tetra-, and two di-phosphonate cavitands, which differ according to the number and position of H-bond acceptor P = O groups, was evaluated by using different competition experiments, energy-resolved CID, gas-phase H/D-exchange, and ligand-exchange reactions, together with ab initio theoretical optimization of the complexes. The phosphonate cavitands with two or more adjacent P = O groups were found to be selective towards secondary alkylammonium ions, due to simultaneous formation of two stable hydrogen bonds. In the ion-molecule reactions (both H/D- and ligand-exchange), the formation of two stable hydrogen bonds was observed either to slow down the reaction or to completely prevent it. This was, however, limited to situations where two hydrogen bonds are formed between the H-bond donor sites of the alkyl ammonium ion and the vicinal H-bond acceptor sites of the cavitand.     

Water versus acetonitrile coordination to uranyl. Density functional study of cooperative polarization effects in solution

Optimizations at the BLYP and B3LYP levels are reported for mixed uranyl-water/acetonitrile complexes [UO(2)(H(2)O)(5-n)(MeCN)(n)](2+) (n = 0-5), in both the gas phase and a polarizable continuum modeling acetonitrile. Car-Parrinello molecular dynamics (CPMD) simulations have been performed for these complexes in the gas phase, and for selected species (n = 0, 1, 3, 5) in a periodic box of liquid acetonitrile. According to structural and energetic data, uranyl has a higher affinity for acetonitrile than for water in the gas phase, in keeping with the higher dipole moment and polarizability of acetonitrile. In acetonitrile solution, however, water is the better ligand because of specific solvation effects. Analysis of the dipole moment of the coordinated water molecule in [UO(2)(H(2)O)(MeCN)(4)](2+) reveals that the interaction with the second-shell solvent molecules (through fairly strong and persistent O-H···N hydrogen bonds) causes a significant increase of this dipole moment (by more than 1 D). This cooperative polarization of water reinforces the uranyl-water bond as well as the water solvation via strengthened (UO(2))OH(2)···NCMe hydrogen bonds. Such cooperativity is essentially absent in the acetonitrile ligands that make much weaker (UO(2))NCMe···NCMe hydrogen bonds. Beyond the uranyl case, this study points to the importance of cooperative polarization effects to enhance the M(n+) ion affinity for water in condensed phases involving M(n+)-OH(2)···A fragments, where A is a H-bond proton acceptor and M(n+) is a hard cation.     

An ab-initio study of the C3H6-HX, C2H4-HX and C2H2-HX hydrogen-bonded complexes with X=F or Cl.

MP2/6-31G** ab-initio molecular orbital calculations have been performed to obtain geometries, H-bond energies and vibrational properties of the C3H6-HX, C2H4-HX and C2H2-HX H-bonded complexes with X=F or Cl. The more pronounced effects on the structural parameters of the isolated molecules due to complexation are verified to the CC and HX bond lengths, which are directly involved in the H-bond formation. They are increased after complexation. The calculated H-bond lengths for the hydrogen complexes for X=F are shorter than those for x-Cl by about 0.55 A, whereas the corresponding experimental value is 0.58 A. The H-bond energies are essentially determined by the nature of the proton donor molecule. For X=F, the AE mean value is 20 kJ/mol, whereas it is approximately 14.5 kJ/mol for X-Cl. The H-bond energies including zero-point corrections show a good correlation with the H-bond lengths. The more pronounced effect on the normal modes of the isolated molecules after complexation occurs to the H-X stretching mode. The H-X stretching frequency is shifted downward, whereas its IR intensity is much enhanced upon H-bond formation. The new vibrational modes arising from complexation show several interesting features.     

Cooperative interactions of hydrogen bonds in proton-transfer processes involving water molecules. Simulation of biochemical systems

Quantum-chemical calculations of molecular complexes simulating the proton channel of influenza A virus and the proton-transfer system of the active site of carboanhydrase enzyme were performed. These complexes comprise a proton-donor and a proton-acceptor groups bridged by a chain of water molecules. Calculations of the methylimidazole (H⁺)-H₂O-CH₃COO^? complex as a model of influenza M₂ virus revealed free translation motion of the water molecule between the donor and acceptor, as well as concerted proton transfer in both H bonds. The barrier for proton transfer is independent of the position of the bridging water molecule and varies linearly with the difference in the electrostatic potentials between the donor and acceptor. With elongation of the H-bond bridge between the donor and acceptor groups, the H-bond lengths and proton shifts in the chain links vary periodically. This process can be defined as an H-bond deformation wave (proton wave). It was shown that motion of one proton along the H bond is associated with vibrational motion of protons in other links, which results in wave propagation along the chain. The calculation results allowed the rate of the proton wave and the time of proton transfer from the donor to acceptor to be estimated.     

Direct observation and control of hydrogen-bond dynamics using low-temperature scanning tunneling microscopy

Hydrogen(H)-bond dynamics are involved in many elementary processes in chemistry and biology. Because of its fundamental importance, a variety of experimental and theoretical approaches have been employed to study the dynamics in gas, liquid, solid phases, and their interfaces. This review describes the recent progress of direct observation and control of H-bond dynamics in several model systems on a metal surface by using low-temperature scanning tunneling microscopy (STM). General aspects of H-bond dynamics and the experimental methods are briefly described in chapter 1 and 2. In the subsequent four chapters, I present direct observation of an H-bond exchange reaction within a single water dimer (chapter 3), a symmetric H bond (chapter 4) and H-atom relay reactions (chapter 5) within water–hydroxyl complexes, and an intramolecular H-atom transfer reaction (tautomerization) within a single porphycene molecule (chapter 6). These results provide novel microscopic insights into H-bond dynamics at the single-molecule level, and highlight significant impact on the process from quantum effects, namely tunneling and zero-point vibration, resulting from the small mass of H atom. Additionally, local environmental effect on H-bond dynamics is also examined by using atom/molecule manipulation with the STM.     

Two-way effects between hydrogen bond and intramolecular resonance effect: an ab initio study on complexes of formamide and its derivatives with water

Ab initio calculations up to MP2/aug-cc-pVTZ//MP2/cc-pVTZ level, including natural charge population and natural resonance theory analyses, have been carried out to study the two-way effects between hydrogen bond (H-bond) and the intramolecular resonance effect by using the H-bonded complexes of formamide ( FAO) and its derivatives ( FAXs, X represents the heavy atoms in the substituent groups, CH 2, NH, SiH 2, PH, and S) with water as models. Unlike NH 3 and NH 2CH 3 which prefer being H-bond acceptors ( HA) to form H-bond with water, the amino groups in the six monomers, because of the resonance effect, prefer being H-bond donors ( HD) rather HA. Six monomers can all form HD complexes with water, and only two ( FAC and FASi) with the weakest resonance effect are able to form HA complexes with water. The HD H-bond and resonance effect enhance each other (positive two-way effects) whereas the HA H-bond and resonance effect weaken each other (negative two-way effects). The H-bond energies in the six HD complexes are nearly linearly correlated with the weights of the dipolar resonance in Pauling's model and the N-C bond lengths; the correlation coefficients are 0.91 and 0.93, respectively. The positive two-way effects also happens in FAO-water complex, in which the FAO CO group serves as HA ( HA co ). Interestingly, when the HD and HA co H-bonds are present in FAO H-bond complex simultaneously, the enhancements are much more significant, and the energies of the two types of H-bonds are much larger than those when only one type of H-bond is present, reflecting the cooperative effects. By using the knowledge to the two-way effects, we computationally designed a molecule ( FAO- BH 3 ) to increase H-bond energy. Because of the oxygen lone pair donation to the empty pi orbital of BH 3, FAO- BH 3 has a much stronger resonance effect than FAO. As a result, the H-bond energy (-5.55 kcal/mol) in HD H 2O ... FAO- BH 3 complex is much greater than the -3.30 kcal/mol in the HD H 2O...FAO complex. The two-way effects can be rationalized as follows: the resonance effect leads to intramolecular charge shifts in the monomers which facilitate or prevent the charge donation or acceptation of their H-bond partners. Therefore, the H-bonds are strengthened or weakened. In reverse, the charge donations or acceptations of their H-bond partners facilitate or prevent the intramolecular charge shifts in the monomer moieties, which enhance or weaken the resonance effect. The understanding to the two-way effects may be helpful in drug design and refinement by modulating the H-bond strength and in building empirical H-bond models to study large biological molecules. The study supports Pauling's resonance model.     

Comparison between Hydrogen and Halogen Bonds in Complexes of 6-OX-Fulvene with Pnicogen and Chalcogen Electron Donors

Quantum chemical calculations are applied to complexes of 6-OX-fulvene (X=H, Cl, Br, I) with ZH₃/H₂Y (Z=N, P, As, Sb; Y=O, S, Se, Te) to study the competition between the hydrogen bond and the halogen bond. The H-bond weakens as the base atom grows in size and the associated negative electrostatic potential on the Lewis base atom diminishes. The pattern for the halogen bonds is more complicated. In most cases, the halogen bond is stronger for the heavier halogen atom, and pnicogen electron donors are more strongly bound than chalcogen. Halogen bonds to chalcogen atoms strengthen in the order O<S<Se<Te, whereas the pattern is murkier for the pnicogen donors. In terms of competition, most halogen bonds to pnicogen donors are stronger than their H-bond analogues, but there is no clear pattern with respect to chalcogen donors. O prefers a H-bond, while halogen bonds are favored by Te. For S and Se, I-bonds are strongest, followed Br, H, and Cl-bonds in that order.     

Supercritical Water is not Hydrogen Bonded

Thinking about water is inextricably linked to hydrogen bonds, which are highly directional in character and determine the unique structure of water, in particular its tetrahedral H-bond network. Here, we assess if this common connotation also holds for supercritical water. We employ extensive ab initio molecular dynamics simulations to systematically monitor the evolution of the H-bond network mode of water from room temperature, where it is the hallmark of its fluctuating three-dimensional network structure, to supercritical conditions. Our simulations reveal that the oscillation period required for H-bond vibrations to occur exceeds the lifetime of H-bonds in supercritical water by far. Instead, the corresponding low-frequency intermolecular vibrations of water pairs as seen in supercritical water are found to be well represented by isotropic van-der-Waals interactions only. Based on these findings, we conclude that water in its supercritical phase is not a H-bonded fluid.