Hydrogen bonding without borders: an atoms-in-molecules perspective

It is shown that the electron density at the hydrogen bond critical point increases approximately linearly with increasing stabilization energy in going from weak hydrogen bonds to moderate and strong hydrogen bonds, thus serving as an indicator of the nature and gradual change of strength of the hydrogen bond for a large number of test intermolecular complexes.     

Hydrogen bonding in phenol, water, and phenol-water clusters

Structure, stability, and hydrogen-bonding interaction in phenol, water, and phenol-water clusters have been investigated using ab initio and density functional theoretical (DFT) methods and using various topological features of electron density. Calculated interaction energies at MP2/6-31G level for clusters with similar hydrogen-bonding pattern reveal that intermolecular interaction in phenol clusters is slightly stronger than in water clusters. However, fusion of phenol and water clusters leads to stability that is akin to that of H(2)O clusters. The presence of hydrogen bond critical points (HBCP) and the values of rho(r(c)) and nabla(2)rho(r(c)) at the HBCPs provide an insight into the nature of closed shell interaction in hydrogen-bonded clusters. It is shown that the calculated values of total rho(r(c)) and nabla(2)rho(r(c)) of all the clusters vary linearly with the interaction energy.     

Hydrogen peroxide and ammonia on protonated ice clusters

A temperature controlled source for protonated water clusters has been combined with high-resolution mass spectroscopy to study the stability pattern of ice clusters and compounds with ammonia and hydrogen peroxide depending on temperature. The stability pattern of pure protonated ice shows the two well known peaks at 21 and 28 molecules and also less pronounced structure up to n=55. Ammonia and hydrogen peroxide do not destroy this pattern but shift it by a number of water molecules. The additives are therefore integrated in the persisting crystalline structure of the pure protonated ice. Based on this structural information, density functional theory calculations reveal that hydrogen peroxide and ammonia occupy surface positions on a dodecahedral 21-molecule cluster and are not caged in the center.     

Gas-phase calorimetry of protonated water clusters

Protonated water clusters with 60 to 79 molecules have been studied by nanocalorimetry. The technique is based on multi-collision excitations of the accelerated clusters with helium. The caloric curves indicate transitions that resemble those of water clusters charged by an excess electron, but the transition temperatures of the protonated clusters are higher.     

Hydrogen bonding in clusters and molecular crystals

Ab initio calculations on the Hartree–Fock level are used to study the effects of cooperativity in ternary complexes and in infinite chains on the example of hydrogen fluoride. SCF energy partitioning demonstrates that polarization forces dominate three-body energies in these aggregates. The crystal orbital method is applied to investigate the structural differences between dimers and molecular crystals.     

Fragmentation cross sections of protonated water clusters

We have measured fragmentation cross sections of protonated water cluster cations (H(2)O)(n=30-50)H(+) by collision with water molecules. The clusters have well-defined sizes and internal energies. The collision energy has been varied from 0.5 to 300 eV. We also performed the same measurements on deuterated water clusters (D(2)O)(n=5-45)D(+) colliding with deuterated water molecules. The main fragmentation channel is shown to be a sequential thermal evaporation of single molecules following an initial transfer of relative kinetic energy into internal energy of the cluster. Unexpectedly, that initial transfer is very low on average, of the order of 1% of collision energy. We evaluate that for direct collisions (i.e., within the hard sphere radius), the probability for observing no fragmentation at all is more than 35%, independently of cluster size and collision energy, over our range of study. Such an effect is well known at higher energies, where it is attributed to electronic effects, but has been reported only in a theoretical study of the collision of helium atoms with sodium clusters in that energy range, where only vibrational excitation occurs.     

Attachment cross sections of protonated water clusters

The attachment of water molecules onto size selected protonated water clusters has been experimentally investigated. Absolute attachment cross sections are measured as a function of cluster size, collision energy, and initial cluster temperature. Although thermal evaporation is ruled out in our experiment, attachment cross sections become significantly smaller than hard sphere cross sections as the collision energy increases. This feature is attributed to a transition from adiabatic to nonadiabatic regime. It is shown to be due to a dynamical effect: as the collision duration becomes shorter than the typical time required for collision energy redistribution into clusters internal energy, the attachment probability is reduced. We relate this typical time to the period of the main surface vibrational mode excited by the collisions. This hypothesis is further supported by results obtained with deuterated water clusters.     

Infrared spectroscopy of small protonated water clusters at room temperature: an effective modes analysis

We perform infrared vibrational analysis on small protonated water clusters H(+)(H(2)O)(n), with n = 2, [ellipsis (horizontal)], 6, at room temperature. The absorption spectra are calculated based on classical trajectories obtained by the multistate empirical valence bond method. The analysis is carried out based on the effective modes analysis, which has been recently developed [Martinez et al., J. Chem. Phys. 125, 144106 (2006)] as generalization of the normal modes analysis. This technique enables us to decompose the full spectrum in maximally localized bands which are obtained by accounting for temperature and anharmonic effects. These effects are especially considered in the determination of the modes coupling. The spectra of the small clusters are interpreted by identifying the behavior of the excess charge, by understanding the role of hydrogen bonds, and by considering the effect of (micro-)solvation. Our results are presented by showing comparisons with other numerical methods and experimental measurements which are available in the literature.     

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.     

Structure of protonated water clusters: low-energy structures and finite temperature behavior

The structure of protonated water clusters H+(H2O)n (n=5-22) are examined by two Monte Carlo methods in conjunction with the OSS2 potential [L. Ojamae, I. Shavitt, and S. J. Singer J. Chem. Phys. 109, 5547 (1998)]. The basin-hopping method is employed to explore the OSS2 potential energy surface and to locate low-energy structures. The topology of the "global minimum," the most stable low-energy structure, changes from single ring to multiple ring to polyhedral cage as the cluster size grows. The temperature dependence of the cluster geometry is examined by carrying out parallel tempering Monte Carlo simulations. Over the temperature range we studied (25-330 K), all water clusters undergo significant structural changes. The trends are treelike structures dominating at high temperature and single-ring structures appearing in slightly lower temperatures. For n> or =7, an additional transition from single ring to multiple rings appears as the temperature decreases. Only for n> or =16 do polyhedral structures dominate the lowest temperature range. Our results indicate very dynamic structural changes at temperature range relevant to atmospheric chemistry and current experiments. The structures and properties of medium-sized protonated clusters in this temperature range are far from their global minimum cousins. The relevance of these findings to recent experiments and theoretical simulations is also discussed.     

Theoretical study of hydrogen bonded clusters of water and cyanic acid: Hydrogen bonding in terms of the molecular structure

Ab initio and density functional calculations were used to analyze the interaction between a molecule of cyanic acid (HOCN) and up to 4 molecules of water at the B3LYP/6-311++G(d,p) and MP2/6-311++G(d,p) computational levels. The cooperative effect (CE) is increased with the increasing size of the studied clusters. Red shifts of the H–O stretching frequency for complexes involving HOCN as an H-donor were predicted. The strength of the hydrogen bonds in terms of molecular structures could be deduced from a comparison of HOCN–H₂O with HCNO–H₂O, HONC–H₂O and HNCO–H₂O HB clusters. The atom in molecules (AIM) method was used to analyze the cooperative effects on topological parameters.     

Hydrogen bonding definitions and dynamics in liquid water

X-ray and neutron diffractions, vibrational spectroscopy, and x-ray Raman scattering and absorption experiments on water are often interpreted in terms of hydrogen bonding. To this end a number of geometric definitions of hydrogen bonding in water have been developed. While all definitions of hydrogen bonding are to some extent arbitrary, those involving one distance and one angle for a given water dimer are unnecessarily so. In this paper the authors develop a systematic procedure based on two-dimensional potentials of mean force for defining cutoffs for a given pair of distance and angular coordinates. They also develop an electronic structure-based definition of hydrogen bonding in liquid water, related to the electronic occupancy of the antibonding OH orbitals. This definition turns out to be reasonably compatible with one of the distance-angle geometric definitions. These two definitions lead to an estimate of the number of hydrogen bonds per molecule in liquid simple point charge/extended (SPC/E) water of between 3.2 and 3.4. They also used these and other hydrogen-bond definitions to examine the dynamics of local hydrogen-bond number fluctuations, finding an approximate long-time decay constant for SPC/E water of between 0.8 and 0.9 ps, which corresponds to the time scale for local structural relaxation.     

Hydrogen isotope effects on covalent and noncovalent interactions: The case of protonated rare gas clusters

We investigate hydrogen isotope and nuclear quantum effects on geometries and binding energies of small protonated rare gas clusters (Rg $_n$ X $^ +$ , Rg = He,Ne,Ar, X = H,D,T, and $n$ = 1–3) with the any particle molecular orbital (APMO) MP2 level of theory (APMO/MP2). To gain insight on the impact of nuclear quantum effects on the different interactions present in the Rg $_n$ X $^ +$ systems, we propose an APMO/MP2 energy decomposition analysis scheme. For RgH $^ +$ ions, isotopic substitution leads to an increase in the stability of the complex, because polarization and charge transfer contributions increase with the mass of the hydrogen. In the case of Rg $_2$ H $^ +$ complexes, isotopic substitution results in a shortening and weakening of the rare gas‐hydrogen ion bond. For Rg $_3$ X $^ +$ complexes, the isotope effects on the rare gas binding energy are almost negligible. Nevertheless, our results reveal that subtle changes in the charge distribution of the Rg $_2$ X $^ +$ core induced by an isotopic substitution have an impact on the geometry of the Rg $_3$ X $^ +$ complex. © 2012 Wiley Periodicals, Inc.     

Eigen and Zundel forms of small protonated water clusters: structures and infrared spectra

The spectral properties of protonated water clusters, especially the difference between Eigen (H3O+) and Zundel (H5O2+) conformers and the difference between their unhydrated and dominant hydrated forms are investigated with the first principles molecular dynamics simulations as well as with the high level ab initio calculations. The vibrational modes of the excess proton in H3O+ are sensitive to the hydration, while those in H5O2+ are sensitive to the messenger atom such as Ar (which was assumed to be weakly bound to the water cluster during acquisitions of experimental spectra). The spectral feature around approximately 2700 cm-1 (experimental value: 2665 cm-1) for the Eigen moiety appears when H3O+ is hydrated. This feature corresponds to the hydrating water interacting with H3O+, so it cannot appear in the Eigen core. Thus, H3O+ alone would be somewhat different from the Eigen forms in water. For the Zundel form (in particular, H5O2+), there have been some differences in spectral features among different experiments as well as between experiments and theory. When an Ar messenger atom is introduced at a specific temperature corresponding to the experimental condition, the calculated vibrational spectra for H5O2+.Ar are in good agreement with the experimental infrared spectra showing the characteristic Zundel frequency at approximately 1770 cm-1. Thus, the effect of hydration, messenger atom Ar, and temperature are crucial to elucidating the nature of vibrational spectra of Eigen and Zundel forms and to assigning the vibrational modes of small protonated water clusters.     

Hydrogen-hydrogen bonding in planar biphenyl, predicted by atoms-in-molecules theory, does not exist

Based on an Atoms-in-Molecules (AIM) analysis, Matta et al. recently claimed evidence for the existence of hydrogen-hydrogen bonding between ortho-hydrogen atoms, pointing towards each other from adjacent phenyl groups in planar biphenyl. This AIM result is opposed to the classical view that nonbonded steric repulsion between the ortho-hydrogen atoms is responsible for the higher energy of the planar as compared to the twisted geometry of biphenyl. In the present work, we address the question if hydrogen-hydrogen bonding in biphenyl exists, as suggested by AIM, or not. To this end, we have analyzed the potential energy surface for internal rotation of biphenyl in terms of two interacting phenyl radicals using density functional theory (DFT) at BP86/TZ2P. A detailed analysis of the bonding mechanism and a quantitative bond energy decomposition in the framework of Kohn-Sham DFT show that Pauli (or overlap) repulsion, mainly between C(ortho)--H(ortho) phenyl MOs, prevents biphenyl from being planar and forces it to adopt a twisted equilibrium geometry. Furthermore, a derivative of biphenyl in which all four ortho-hydrogen atoms have been removed does adopt a planar equilibrium geometry. Thus, our results confirm the classical view of steric repulsion between ortho-hydrogen atoms in biphenyl and they falsify the hypothesis of hydrogen-hydrogen bonding.