SCF-CI studies of correlation effects on hydrogen bonding and ion hydration

Previous single-determinant Hartree-Fock studies on the equilibrium structures and stabilities of H₂ O, H₃ O⁺ as well as of the monohydrated ionic systems Li⁺ · H₂O, F^? · H₂O and the hydrogen bonded water dimer, H₂ O · HOH, are extended by large scale configuration interaction calculations including all the possible single and double excitations arising from the canonical set of Hartree-Fock molecular orbitals. The correlation energy effects on the equilibrium geometrical parameters of the systems under consideration are found to be quite small. The contributions of the correlation energy to the total binding energies of the weakly interacting composed systems are obtained to be of the order of 1 kcal/mole, leading to a considerable increase of the hydrogen bond strength in F^? · H₂O and H₂O · HOH and to a small decrease of the binding energy in Li⁺ · H₂ O. The observed strengthening of the hydrogen bonding interaction due to correlation is shown to be partly compensated by the change in the vibrational zero-point energy of the composed systems compared to the non-interacting subsystems. Approximate force constants corresponding to the intersystem vibrations in Li⁺ · H₂O, F^? · H₂ O, and H₂O · HOH are deduced from the calculated potential curve data on the SCF and the CI level of accuracy.     

The role of hydrogen bonding in the crystal structures of zinc phosphate hydrates

The compounds alpha- and beta-hopeite have been synthesised by hydrothermal crystallisation from aqueous solution at 90 degrees C and 20 degrees C, respectively. The crystal structures of these polymorphic forms of zinc phosphate tetrahydrate (ZPT), Zn(3)(PO(4))(2).4 H(2)O, have been resolved. Single-crystal analysis proves that the main difference between the alpha and beta forms of ZPT is caused by the difference in orientation of one of the water molecules in the ZnO(6) octahedral network, indicating two different hydrogen-bonding patterns. A previously unknown hopeite, Zn(3)(HPO(4))(3).3 H(2)O (ZHPT), has been isolated and analysed. This helps to achieve a better understanding of the mechanism of formation of zinc phosphate compounds. Unambiguous identification of each phase is established by analysis of their unique thermal behaviour and thermodynamic interrelationship.     

Synergistic effect between metal coordination and hydrogen bonding in phosphate and halide recognition

The equilibrium constant for binding of dimethyl phosphate to a Co(III) complex in water increases from 6.2 to 210 M-1 upon addition of a single hydrogen bond between the bound phosphate and the metal complex. Crystal structure reveals that the hydrogen bond distance is 1.96 A. The synergistic effect between metal coordination and hydrogen bonding can also be observed for fluoride binding but not for bromide binding.     

The molecular mobility of 2-propenylphenols capable of hydrogen bonding

Hydrogen bonding and molecular mobility of 2-propenylphenols containing an aminomethyl fragment in the ortho position were studied by NMR. The results confirmed the importance of using relaxation times as a sensitive test for molecular dynamics studies.     

Hydrogen bonding mediated by key orbital interactions determines hydration enthalpy differences of phosphate water clusters

Electronic structure calculations have been carried out to provide a molecular interpretation for dihydrogen phosphate stability in water relative to that of metaphosphate. Specifically, hydration enthalpies of biologically important metaphosphate and dihydrogen phosphate with one to three waters have been computed with second-order M?ller-Plesset perturbation and density functional theory (B3LYP) with up to the aug-cc-pvtz basis set and compared to experiment. The inclusion of basis set superposition error corrections and supplemental diffuse functions are necessary to predict hydration enthalpies within experimental uncertainty. Natural bond orbital analysis is used to rationalize underlying hydrogen bond configurations and key orbital interactions responsible for the experimentally reported difference in hydration enthalpies between metaphosphate and dihydrogen phosphate. In general, dihydrogen phosphate forms stronger hydrogen bonds compared to metaphosphate due to a greater charge transfer or enhanced orbital overlap between the phosphoryl oxygen lone pairs, n(O), and the antibonding O-H bond of water. Intramolecular distal lone pair repulsion with the donor n(O) orbital of dihydrogen phosphate distorts symmetric conformations, which improves n(O) and sigma*(O-H) overlap and ultimately the hydrogen bond strength. Unlike metaphosphate, water complexed to dihydrogen phosphate can serve as both a hydrogen bond donor and a hydrogen bond acceptor, which results in cooperative charge transfer and a reduction of the energy gap between n(O) and sigma*(O-H), leading to stronger hydrogen bonds. This study offers insight into how orbital interactions mediate hydrogen bond strengths with potential implications on the understanding of the kinetics and mechanism in enzymatic phosphoryl transfer reactions.     

DFT and FT-IR analyses of hydrogen bonding in 3-substistuted-3-oxo-arylhydrazonopropanenitriles

Ab initio calculations, FT-IR and X-ray crystal analysis, indicated that the most stable configuration of 3-oxo-2-(phenylhydrazono)-3-(thien-2-yl)-propionitrile is the anti phenylhydrazone structure 1. Stability of such a conformation, over the possible E-form, 2, that would be stabilized by intramolecular hydrogen bonding, is due to interaction between electron-pair domains of the N, S and O atoms. However, the simulated and experimental IR frequency data indicated intermolecular hydrogen bonding between NH and CN, the latter being lowered to 2214 cm(-1). Studies on 3-oxo-3-phenyl-2-(phenylhydrazono)-propionitrile showed the same result, as well as, another intramolecular hydrogen association of the type N-H...O. This was clearly indicated by the absorbance of the carbonyl stretch at 1605 cm(-1). These data indicated the existence of a bifurcated hydrogen bond in 1a and a single intermolecular association in 1b.     

Elucidation of structural restraints for phosphate residues with different hydrogen bonding and ionization states

Solid state NMR spectroscopy and gauge including atomic orbital (GIAO) theoretical calculations were employed to establish structural restraints (ionization, hydrogen bonding, spatial arrangement) for O-phosphorylated l-threonine derivatives in different ionization states and hydrogen bonding strengths. These structural restraints are invaluable in molecular modeling and docking procedures for biological species containing phosphoryl groups. Both the experimental and the GIAO approach show that 31P delta ii chemical shift tensor parameters are very sensitive to the ionization state. The negative values found for the skew kappa are typical for -2 phosphates. The distinct span Omega values reflect the change of strength of hydrogen bonding. For species in the -1 ionization state, engaged in very strong hydrogen bonds, Omega is smaller than for a phosphate group involved in weak hydrogen bonding. For phosphates in the -2 ionization state, Omega is significantly smaller compared to -1 species, although the kappa for -1 samples never reaches negative values. For -1 phosphate residues, in the case when 1H one pulse and/or combined rotation and multiple pulse spectroscopy (CRAMPS) sequences fail and assignment of proton chemical shift is ambiguous, a combination of 1H-(13)C and 1H-(31)P 2D heteronuclear correlation (HETCOR) correlations is found to be an excellent tool for the elucidation of 1H isotropic chemical shifts. In addition, a 2D strategy using 1H-(1)H double quantum filter (DQF) correlations [a back-to-back (BABA) sequence in this work] is useful for analyzing the topology of hydrogen bonding. In the case of a multicenter phosphorus domain, 2D 31P-(31)P proton driven spin diffusion experiments give information about the spatial arrangement of the phosphate residues.     

Raman study of hydrogen bonding in the 2-haloethanols

The high frequency (3000–3700 cm⁻¹) Raman spectra of the 2-haloethanols (XCH₂CH₂OH₂ X  F, Cl, Br and I) were studied as a function of temperature in the neat liquid phase. The two bands in this region, at 3300 and 3600 cm⁻¹, were assigned to the valence mode of inter- and intramolecularly hydrogen bonded OH groups respectively. The intensity ratio, I_intra/Iinter, determined from the resolved band parameters, increases with temperature; calculated values of ΔH exhibit the trend, IBr>Cl>F. A simple picture is introduced which shows these results to be consistent with the order of intramolecular hydrogen bond strengths deduced from an earlier investigation of these alcohols in dilute solution.     

Intra- and intermolecular hydrogen bonding in 3-hydroxy- and 5-hydroxychromone

The crystal structures for 3-hydroxychromone and 5-hydroxychromone have been obtained. Both molecules exhibit intramolecular hydrogen bonding between the hydroxy group and the ketone oxygen atom. However, only the 3-hydroxy derivative contains hydrogen bonds between molecules. By comparing the current results with those obtained for the corresponding flavone derivatives, the effect of the B phenyl group on hydrogen bonding is inferred.     

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.     

SCF MO LCGO Studies on hydrogen bonding: The system NH3 · H2O

The energy hypersurface of the system NH₃ · H₂O is investigated for a number of different internuclear geometries. In the minimum energy structure involving a linear hydrogen bond, NH₃ acts as proton acceptor. The binding energy of the system is calculated to be 6.28 kcal/mole and the bond distance d(NO) to be 3.07 Å. The potential energy curve of the inversion of the hydrogenbonded NH₃ is computed and discussed.

---

Zusammenfassung Die Energiehyperflächen des NH₃ · H₂O-Systems wurden für eine Anzahl von verschiedenen geometrischen Anordnungen untersucht. Im Falle der Struktur minimaler Energie wird eine lineare Wasserstoffbindung gebildet, NH₃ wirkt als Protonakzeptor. Die Berechnungen ergeben eine Bindungsenergie des Systems von 6,28 kcal/Mol und einen NO-Abstand von 3,07 Å. Außerdem wurde die Potentialkurve für die Inversion des über eine Wasserstoffbrücke gebundenen NH₃ berechnet und diskutiert.

---

---

It is a pleasure to thank our technical staff for the careful preparation of the input for the programs and for its skillful assistance in running the computer.     

The role of hydrogen bonding in water-metal interactions

The hydrogen bond interaction between water molecules adsorbed on a Pd <111> surface, a nucleator of two dimensional ordered water arrays at low temperatures, is studied using density functional theory calculations. The role of the exchange and correlation density functional in the characterization of both the hydrogen bond and the water-metal interaction is analyzed in detail. The effect of non local correlations using the van der Waals density functional proposed by Dion et al. [M. Dion, H. Rydberg, E. Schr?der, D. C. Langreth and B. I. Lundqvist, Phys. Rev. Lett., 2004, 92, 246401] is also studied. We conclude that the choice of this potential is critical in determining the cohesive energy of water-metal complexes. We show that the interaction between water molecules and the metal surface is as sensitive to the density functional choice as hydrogen bonds between water molecules are. The reason for this is that the two interactions are very similar in nature. We make a detailed analogy between the water-water bond in the water dimer and the water-Pd bond at the Pd <111> surface. Our results show a strong similarity between these two interactions and based on this we describe the water-Pd bond as a hydrogen bond type interaction. These results demonstrate the need to obtain an accurate and reliable representation of the hydrogen bond interaction in density functional theory.     

The role of hydrogen bonding in supercooled methanol

The role of hydrogen bonding on the microscopic properties of supercooled methanol has been analyzed by means of molecular dynamics simulations. Thermodynamic, structural, and dynamical properties have been investigated in supercooled methanol. The results have been compared with those of an ideal methanol-like system whose molecules have the same dipole moment as the methanol but lack sites for hydrogen bonding. Upon cooling the methanol samples, translational relaxation times increase more rapidly than reorientational ones. This effect is much more important when hydrogen bonds are suppressed. Suppression of hydrogen bonds also results in lower critical temperatures for diffusion and for several characteristic relaxation time constants. The anisotropy of individual dynamics and the existence of dynamical heterogeneities have also been investigated.     

The reaction of phosphate with aluminum hydroxide in relation with phosphate bonding in soils

The retention of phosphate in soils is to a large extent determined by insoluble metal-oxides. In this paper, model studies are described on the uptake and release of phosphate by gibbsite (γ-Al(OH)₃). At low phosphate concentration c_p, the uptake is a pure monolayer exchange adsorption on the edge surfaces of the gibbsite platelets. At high c_p and long reaction times, a precipitate is formed on top of this adsorbate. Kinetic equations are given to discriminate between these processes.Our analysis leads to arguments in favour of applying slow-release phosphate fertilizers. In addition we analyse the effect on the precipitation of some ionic admixtures such as K⁺, Na⁺ and NH⁺₄ that may be inadvertedly or intentionally present in the soil.     

Vibrations and hydrogen bonding in porphycene

Combined use of IR, Raman, neutron scattering and fluorescence measurements for porphycene isolated in helium nanodroplets, supersonic jet and cryogenic matrices, as well as for solid and liquid solutions, resulted in the assignments of almost all of 108 fundamental vibrations. The puzzling feature of porphycene is the apparent lack of the N-H stretching band in the IR spectrum, predicted to be the strongest of all bands by standard harmonic calculations. Theoretical modeling of the IR spectra, based on ab initio molecular dynamics simulations, reveals that the N-H stretching mode should appear as an extremely broad band in the 2250-3000 cm(-1) region. Coupling of the N-H stretching vibration to other modes is discussed in the context of multidimensional character of intramolecular double hydrogen transfer in porphycene. The analysis can be generalized to other strongly hydrogen-bonded systems.