Cooperativity and electron correlation effects on hydrogen bonding in infinite systems

Structural and electronic properties of hydrogen-bonded infinite chains of hydrogen cyanide and formamide molecules have been investigated by the ab initio crystal orbital method using several, partly highly polarized, atomic basis sets of increasing size at the Hartree–Fock (HF ) level and by including electron correlation effects in the second order of Møller–Plesset perturbation theory. The results obtained show that hydrogen bonding in molecular crystals of the type investigated is a highly cooperative phenomenon, both from the structural and energetic points of view. Comparison with clusters of up to four monomers demonstrate how various structural parameters converge toward their limiting values in the infinite system. The results obtained for infinite HCN chains show an excellent agreement with those observed for solid HCN, whereas the infinite formamide chain proves to be a reasonable model for the corresponding liquid phase. © 1994 John Wiley & Sons, Inc.     

SCF LCGO MO studies on the fluoronium ion FH 2 + and its hydrogen bonding interaction with hydrogen fluoride FH

The stability and geometrical structure of the fluoronium ion is investigated using the onedeterminant SCF LCAO MO method. The equilibrium geometry is characterized by a bond length of d(FH)=0.95 Å and a bond angle of 114.75°. The proton binding energy is determined to be 120.1 kcal/mole. The molecules FH ₃ ²⁺ and FH₃ are found to be unstable. A binding energy of 30.7 kcal/mole is obtained for the hydrogen bond formation between the systems FH ₂ ⁺ and FH. In the minimum energy structure the central proton is situated midway between the two F atoms in a symmetrical single minimum potential. The general behavior of the potential curves of the di-solvated proton involving NH₃, OH₂, and FH as solvent molecules is discussed. In all these cases double minimum potentials are found, if the equilibrium separation between the heavy atoms is larger than approximately 2.4 Å, and single minimum potential for separations smaller than this value.     

1H NMR studies of solvent effects on hydrogen bonding in some pyridine trifluoroacetates

The chemical shifts of hydrogen bonded protons in complexes of 11 substituted pyridines with trifluoroacetic acid were examined, in five dry solvents of different activity, with respect to proton transfer and aggregation effects. The results were correlated with ΔpK_a, the Kirkwood function and E _T parameters. The solvent effect on the intersection point obtained from the plot of the chemical shift of the hydrogen bonded protons against ΔpK_a can be used, similar to an isotopic effect, to differentiate strong hydrogen bonds. The aggregation of acid–base complexes can lead to downfield or upfield shifts; the variation of chemical shift with aggregation depends on the position of the proton in the hydrogen bridge.     

Analysis of nuclear quantum effects on hydrogen bonding

The impact of nuclear quantum effects on hydrogen bonding is investigated for a series of hydrogen fluoride (HF)n clusters and a partially solvated fluoride anion, F-(H2O). The nuclear quantum effects are included using the path integral formalism in conjunction with the Car-Parrinello molecular dynamics (PICPMD) method and using the second-order vibrational perturbation theory (VPT2) approach. For the HF clusters, a directional change in the impact of nuclear quantum effects on the hydrogen-bonding strength is observed as the clusters evolve toward the condensed phase. Specifically, the inclusion of nuclear quantum effects increases the F-F distances for the (HF)n=2-4 clusters and decreases the F-F distances for the (HF)n>4 clusters. This directional change occurs because the enhanced electrostatic interactions between the HF monomers become more dominant than the zero point energy effects of librational modes as the size of the HF clusters increases. For the F-(H2O) system, the inclusion of nuclear quantum effects decreases the F-O distance and strengthens the hydrogen bonding interaction between the fluoride anion and the water molecule because of enhanced electrostatic interactions. The vibrationally averaged 19F shielding constant for F-(H2O) is significantly lower than the value for the equilibrium geometry, indicating that the electronic density on the fluorine decreases as a result of the quantum delocalization of the shared hydrogen. Deuteration of this system leads to an increase in the vibrationally averaged F-O distance and nuclear magnetic shielding constant because of the smaller degree of quantum delocalization for deuterium.     

SCF-MO-LCGO studies on hydrogen bonding. The water dimer

The hydrogen bond in the water dimer is studied within the SCF-MO-LCAO framework, using a large Gaussian basis set to approximate the wavefunction. A geometry search restricted to structures with linear and bifurcated hydrogen bonds is performed and the associated potential energy curves are displayed. The minimum energy geometry of the water dimer is found to form a linear hydrogen bond with a hydrogen bond distance of 2.04 Å and a binding energy of 4.84 kcal/mole relative to the monomer (exp. 5.0 kcal/mole). No semistable structures are found. The charge density and charge density difference maps are discussed for the structure with a linear hydrogen bond for different subsystem (water) separations, including the minimum energy geometry. The dipole moment of the dimer is computed to be 1.69 a.u. The shift of the IR bands on hydrogen bond formation is explained qualitatively by comparing the potential energy curves of the hydrogen in the OH-bonds of the monomer and the dimer, and the intensity increase of the fundamental OH-stretching band is computed. The shift of the proton magnetic resonance signal is discussed qualitatively by inspecting the charge density change on hydrogen bond formation, and the average diamagnetic shielding is calculated.

---

Zusammenfassung Die Wasserstoffbrückenbindung im dimeren Wasser ist im Rahmen des SCF-MO-LCAO-Verfahrens untersucht worden, wobei die Wellenfunktion durch Gaußfunktionen angenähert wurde. Die Untersuchungen beschränkten sich auf Strukturen mit linearen und gegabelten Wasserstoffbrücken. Die zugehörigen Potentialkurven wurden berechnet und graphisch dargestellt. Danach besitzt das dimere Wasser in der energetisch stabilsten Form eine lineare Wasserstoffbrückenbindung mit einem Brückenbindungsabstand von 2,04 Å. Die Bildungsenergie aus zwei monomeren Wassermolekülen beträgt 4,84 Kcal/Mol. Es wurden keine semistabilen Strukturen anderer Geometrie gefunden. Für das dimere Wasser mit linearer Wasserstoffbrückenbindung sind für verschiedene Brückenbindungsabstände die Elektronendichten und die Elektronendichtedifferenzen, bezogen auf zwei ungestörte Wassermoleküle als Vergleichssystem, graphisch dargestellt und diskutiert worden. Das Dipolmoment des dimeren Wassers wurde zu 1,69 a.u. berechnet. Die Verschiebung der IR-Banden, die bei der Bildung der Wasserstoffbrückenbindung experimentell beobachtet wird, kann qualitativ aufgrund der entsprechenden Potentialkurven erklärt werden. Die dabei gleichfalls beobachtete Intensitätszunahme wurde berechnet. Die Verschiebung des Proton-Kernresonanzsignals konnte qualitativ an Hand der berechneten Ladungsdichtedifferenzen diskutiert und die diamagnetische Abschirmkonstante bestimmt werden.

---

Résumé Etude de la liaison hydrogène du dimère de l'eau dans le cadre SCF-MO-LCAO en utilisant une grande base de fonctions gaussiennes. La géométrie est recherchée parmi les structures à liaisons hydrogène linéaires et fourchues avec production des courbes d'énergie potentielle associées. La géométrie du dimère d'énergie minimum est du type liaison hydrogène linéaire de longueur 2,04 Å et d'énergie 4,84 Kcal/mole (exp. 5,0 Kcal/mole). On ne trouve pas de structures semi-stables. Les cartes de densité de charge et de densité différentielle sont discutées pour différentes distances de séparation du système à liaison hydrogène linéaire. Le moment dipolaire du dimère est évalué à 1,69 u.a. Le déplacement des bandes I.R par formation de liaison hydrogène est expliqué qualitativement en comparant les courbes d'énergie potentielle de la liaison OH dans le monomère et le dimère, et l'on calcule l'augmentation d'intensité de la bande fondamentale de vibration OH. Le déplacement du signal de résonance magnétique du proton est discuté qualitativement par inspection des variations de densité de charge par formation de la liaison hydrogène, et l'on a calculé l'écran diamagnétique moyen.

---

---

These studies were started while the author was visiting the IBM Research Laboratories, San Jose, California 95114, USA.     

Spectroscopic studies on proton transfer and hydrogen bonding in p-nitrophenol-donor systems

The interaction of p-nitrophenol with several electron-donors has been studied in aprotic and protic solvents by electronic absorption spectroscopy. The equilibrium data for different kinds of equilibria in 1,2-dichloroethane, n-butanol and acetonitrile have been obtained and discussed. The strong hydrogen bonding interaction between p-nitrophenol and a variety of electron donors has been investigated in heptane by electronic absorption spectroscopy. The interesting correlations obtained in terms of Mulliken's charge transfer model have been examined and extended in many cases.     

Nuclear quantum effects and hydrogen bonding in liquids

We have employed ab initio path integral molecular dynamics simulations to investigate the role of nuclear quantum effects on the strength of hydrogen bonds in liquid hydrogen fluoride. Nuclear quantum effects are shown to be responsible for a stronger hydrogen bond and an enhanced dipole-dipole interaction, which lead, in turn, to a shortening of the H...F intrachain distance. The simulation results are analyzed in terms of the electronic density shifts with respect to a purely classical treatment of the nuclei. The observed enhanced hydrogen-bond interaction, which arises from a coupling of intra- and intermolecular effects, should be a general phenomenon occurring in all hydrogen-bonded systems.     

The hydrogen bonding and hydration of 2'-OH in adenosine and adenosine 3'-ethyl phosphate

The 2'-OH group has major structural implications in the recognition, processing, and catalytic properties of RNA. We report here intra- and intermolecular H-bonding of 2'-OH in adenosine 3'-ethyl phosphate (1), 3'-deoxyadenosine (2), and adenosine (3) by both temperature- and concentration-dependent NMR studies, as well as by detailed endo ((3)J(H,H)) and exocyclic ((3)J(H,OH)) coupling constant analyses. We have also examined the nature of hydration and exchange processes of 2'-OH with water by a combination of NOESY and ROESY experiments in DMSO-d(6) containing 2 mol % HOD. The NMR-constrained molecular modeling (by molecular mechanics as well as by ab initio methods both in the gas and solution phase) has been used to characterize the energy minima among the four alternative dihedrals possible from the solution of the Karplus equation for (3)J(H2',OH) and (3)J(H3',OH) to delineate the preferred orientation of 2'-O-H proton in 1 and 2 as well as for 2'/3'-O-H protons in 3. The NMR line shape analysis of 2'-OH gave the DeltaG(H-bond)(298K) of 7.5 kJ mol(-1) for 1 and 8.4 kJ mol(-1) for 3; similar analyses of the methylene protons of 3'-ethyl phosphate moiety in 1 also gave comparable DeltaG(H-bond)(298K) of 7.3 kJ mol(-1). The donor nature of the 2'-OH in the intramolecular H-bonding in 3 is evident from its relatively reduced flexibility [-TDeltaS++](2'-OH) = -17.9(+/-0.5) kJ mol(-1)] because of the loss of conformational freedom owing to the intramolecular 2'O-H...O3' H-bonding, compared to the acceptor 3'-OH in 3 [-TDeltaS++](3'-OH) = -19.8 (+/- 0.6) kJ mol(-1)] at 298 K. The presence of intramolecular 2'-OH...O3' H-bonding in 3 is also corroborated by the existence of weak long-range (4)J(H2',OH3') in 3 (i.e., W conformation of H2'-C2'-C3'-O3'-H) as well as by (3)J(H,OH) dependent orientation of the 2'- and 3'-OH groups. The ROESY spectra for 1 and 3 at 308 K, in DMSO-d(6), show a clear positive ROE contact of both 2'- and 3'-OH with water. The presence of a hydrophilic 3'-phosphate group in 1 causes a much higher water activity in the vicinity of its 2'-OH, which in turn causes the 2'-OH to exchange faster, culminating in a shorter exchange lifetime (tau) for 2'-OH proton with HOD in 1 (tau2'-OH: 489 ms) compared to that in 3 (tau2'-OH: 6897 ms). The activation energy (E(a)) of the exchange with the bound-water for 2'- and 3'-OH in 3 (48.3 and 45.0 kJ mol(-1), respectively) is higher compared to that of 2'-OH in 1 (31.9 kJ mol(-1)), thereby showing that the kinetic availability of hydrated 2'-OH in 1 for any inter- and intramolecular interactions, in general, is owing to the vicinal 3'-phosphate residue. It also suggests that 2'-OH in native RNA can mediate other inter- or intramolecular interactions only in competition with the bound-water, depending upon the specific chemical nature and spatial orientation of other functions with potential for hydrogen bonding in the neighborhood. This availability of the bound water around 2'-OH in RNA would, however, be dictated by whether the vicinal phosphate is exposed to the bulk water or not. This implies that relatively poor hydration around a specific 2'-OH across a polyribonucleotide chain, owing to some hydrophobic microenvironmental pocket around that hydroxyl, may make it more accessible to interact with other donor or acceptor functions for H-bonding interactions, which might then cause the RNA to fold in a specific manner generating a new motif leading to specific recognition and function. Alternatively, a differential hydration of a specific 2'-OH may modulate its nucleophilicity to undergo stereospecific transesterification reaction as encountered in ubiquitous splicing of pre-mRNA to processed RNA or RNA catalysis, in general.     

FTIR studies of intermolecular hydrogen bonding in halogenated ethanols

The effect of halogen substitution on intermolecular hydrogen-bonding in ethanol is studied. Specifically, Fourier-transform infrared (FTIR) spectra of ethanol, 2,2,2-trifluoroethanol (TFE), and 2,2,2-trichloroethanol dissolved in carbon tetrachloride are reported as a function of temperature and concentration. The spectral intensities corresponding to monomer, dimer, and multimer formation are used to determine the effect of halogen substitution on intermolecular hydrogen-bonding. The enthalpy for dimerization was found to evolve from -4.2+/-0.3 kcal/mol in ethanol to -6.8+/-1.0 kcal/mol in TFE. An opposite trend was observed for multimer formation with enthalpies of -3.7+/-0.5 in ethanol and -2.1+/-1.4 kcal/mol in TFE. The majority of this evolution is assigned to the ability of ethanols to form intramolecular hydrogen bonds involving the hydoxyl proton and the halogen substituents.     

Ion-water clusters, bulk medium effects, and ion hydration

Thermochemistry of gas-phase ion-water clusters together with estimates of the hydration free energy of the clusters and the water ligands are used to calculate the hydration free energy of the ion. Often the hydration calculations use a continuum model of the solvent. The primitive quasichemical approximation to the quasichemical theory provides a transparent framework to anchor such efforts. Here we evaluate the approximations inherent in the primitive quasichemical approach and elucidate the different roles of the bulk medium. We find that the bulk medium can stabilize configurations of the cluster that are usually not observed in the gas phase, while also simultaneously lowering the excess chemical potential of the ion. This effect is more pronounced for soft ions. Since the coordination number that minimizes the excess chemical potential of the ion is identified as the optimal or most probable coordination number, for such soft ions the optimum cluster size and the hydration thermodynamics obtained with and without account of the bulk medium on the ion-water clustering reaction can be different. The ideas presented in this work are expected to be relevant to experimental studies that translate thermochemistry of ion-water clusters to the thermodynamics of the hydrated ion and to evolving theoretical approaches that combine high-level calculations on clusters with coarse-grained models of the medium.     

Ab initio SCF-CI studies of the intermolecular interaction between two hydrogen molecules near the Van der Waals minimum

SCF-CI calculations have been used to study the intermolecular energy between two hydrogen molecules in four different geometrical configurations. The CI matrix was diagonalized using perturbation techniques. The importance of the perturbation expansion order upon the intermolecular energy could therefore be studied. The wave function includes all singly and doubly excited configurations. The natural orbitals have been determined and their relative importance on the intermolecular energy is considered.     

基于双嘧啶和双巴比妥酸的氢键组装体的理论研究

用AM1方法对双嘧啶和双巴比妥酸以及二者以1:1的摩尔比通过分子间多重氢键自组织形成的二体进行了几何构型优化,得到氢键键能。并在此基础上用INDO-CI方法讨论了三者的UV-Vis光谱,结果与实验值符合得很好,证实了多重氢键导致新的超分子聚集体的生成,并讨论了氢键在超分子聚集体形成过程中的作用本质。     

Modulation of hydrogen bonding upon ion binding: Insights into cooperativity

The impact due to the of presence of ions, such as Mg²⁺, Na⁺, H⁺, Cl^?, and OH^? on hydrogen bonded clusters of increasing size (water, formamide, and acetamide [n = 1–10]) in the context of associated cooperativity has been explored using density functinal theory (DFT) calculations. Sequential binding energies (SBE) rise on addition of monomer in case of parent clusters. SBE for ionic clusters are several times higher than that of parent clusters initially. This behavior is more dramatic on addition of either Mg²⁺ or H⁺ compared to other ions. Interestingly, SBE of both parent and ionic clusters approach nearly uniform values beyond n = 6 irrespective of kind of ion present in the cluster with the exception of magnesium. © 2013 Wiley Periodicals, Inc.     

A theoretical investigation of hydrogen bonding effects on oxygen and hydrogen chemical shielding tensors of aspirin

A computational investigation was carried out to characterize the ¹⁷O and ¹H chemical shielding (CS) tensors in crystalline aspirin. It was found that O–H⋯O and C–H⋯O hydrogen bonds around the aspirin molecule in the crystal lattice have a different influence on the calculated ¹⁷O and ¹H CS eigenvalues and their orientations in the molecular frame of axes. The calculations were performed with the BLYP, B3LYP, and M06 functionals employing 6-311++G(d,p) standard basis set. Calculated CS tensors were used to evaluate the ¹⁷O and ¹H chemical shift isotropy (δ_iso) and anisotropy (Δσ) in crystalline aspirin, which are in reasonable agreement with available experimental data. The difference between the calculated NMR parameters of the monomer and molecular clusters shows how much hydrogen-bonding interactions affect the CS tensors of each nucleus.     

Miscible blends through hydrogen bonding: effects on polymer properties

Miscible blends through hydrogen bonding have been intensively studied. The effects of a variety of miscible hydrogen bonded polymer blends on properties such as thermal and thermal oxidative stability, moisture sensitivity, modulus and glass transition temperature are discussed. In addition, the preparation of semi-interpenetrating polymer networks (IPNs) and studies of the effect of crosslinking on the miscibility in hydrogen bonded polymer blends are reviewed.     

Density functional studies on the effects of hydrogen bonding on the formation of a charge-transfer complex between p-benzoquinone and 2,6-dimethoxyphenol

The ability to form a ground-state charge-transfer (CT) complex between an electron acceptor, p-benzoquinone (BQ) and an electron donor, 2,6-dimethoxyphenol (DMOPh) was found to be enhanced by H-bonding of BQ to a hydrogen-bond donor, trifluoroacetic acid (TFA) and H-bonding DMOPh to a hydrogen-bond acceptor, 4-(N,N-dimethylamino)pyridine (DMAPy) [Chem. Phys. Lett. 2005, 401, 200]. Here is reported density functional theory (DFT) calculations to study the effect of H-bonding to electron donor and electron acceptor moieties on the ground-state CT complex formation ability between the aforementioned electron donor/acceptor pair. DFT calculations using B3LYP with the 6-311G(d,p) basis set show that the HOMO and LUMO energies of BQ drop on H-bonding to TFA through its C=O groups and the HOMO and LUMO energies of DMOPh increase on H-bonding to DMAPy via its O-H group. BQ molecules hydrogen-bonded as 1:1 and 1:2 complexes to TFA act as stronger acceptors than the bare molecule, while 1:1 complexes of DMOPh and DMAPy act as better donors. Vertical excitation energies for electronic transitions from the ground state to the first few excited states of BQ, DMOPh, DMAPy, and their different complexes have been investigated in the framework of time-dependent density functional theory (TD-DFT) to simulate and interpret experimental ultraviolet absorption spectra. Good agreement between experimental and calculated spectra is established. The enhancement of the CT complex formation ability between the BQ and DMOPh pair is favored by the strong H-bonding interaction of BQ with TFA as well as by the H-bonding interaction of DMOPh with DMAPy.     

Electrical effects on the vibrational transitions of hydrogen fluoride due to hydrogen bonding and applied fields

Derivative Hartree-Fock (DHF) theory has been used to study the transition dipole moment of hydrogen fluoride in several hydrogen-bonded complexes and in the presence of applied fields. DHF is an open-ended, analytical means for finding energy derivatives with respect to any number of parameters, and large basis set calculations going through the seventh dipole hyperpolarizability (ninth derivative) are reported. Using multipole moments, multipole polarizabilities and hyperpolarizabilities, the intermolecular electrical influence on vibrational transitions is analyzed.