This work considers how the properties of hydrogen bonded complexes, X-H· · · Y, are modified by the quantum motion of the shared proton. Using a simple two-diabatic state model Hamiltonian, the analysis of the symmetric case, where the donor (X) and acceptor (Y) have the same proton affinity, is carried out. For quantitative comparisons, a parametrization specific to the O-H· · · O complexes is used. The vibrational energy levels of the one-dimensional ground state adiabatic potential of the model are used to make quantitative comparisons with a vast body of condensed phase data, spanning a donor-acceptor separation (R) range of about 2.4 − 3.0Å, i.e., from strong to weak hydrogen bonds. The position of the proton (which determines the X-H bond length) and its longitudinal vibrational frequency, along with the isotope effects in both are described quantitatively. An analysis of the secondary geometric isotope effect, using a simple extension of the two-state model, yields an improved agreement of the predicted variation with R of frequency isotope effects. The role of bending modes is also considered: their quantum effects compete with those of the stretching mode for weak to moderate Hbond strengths. In spite of the economy in the parametrization of the model used, it offers key insights into the defining features of H-bonds, and semi-quantitatively captures several trends.

We consider how the infrared intensity of an O-H stretch in a hydrogen bonded complex varies as the strength of the H-bond varies from weak to strong. We obtain trends for the fundamental and overtone transitions as a function of donor-acceptor distance R, which is a common measure of H-bond strength. Our calculations use a simple two-diabatic state model that permits symmetric and asymmetric bonds, i.e. where the proton affinity of the donor and acceptor are equal and unequal, respectively. The dipole moment function uses a Mecke form for the free OH dipole moment, associated with the diabatic states. The transition dipole moment is calculated using one-dimensional vibrational eigenstates associated with the H-atom transfer coordinate on the ground state adiabatic surface of our model. Over 20-fold intensity enhancements for the fundamental are found for strong H-bonds, where there are significant non-Condon effects. The isotope effect on the intensity yields a nonmonotonic H/D intensity ratio as a function of R, and is enhanced by the secondary geometric isotope effect. The first overtone intensity is found to vary non-monotonically with H-bond strength; strong enhancements are possible for strong H-bonds. Modifying the dipole moment through the Mecke parameters is found to have a stronger effect on the overtone than the fundamental. We compare our findings with those for specific molecular systems analysed through experiments and theory in earlier works. Our model results compare favourably for strong and medium strength symmetric H-bonds. However, for weak asymmetric bonds we find much smaller effects than in earlier work.

If a deuterated molecule containing strong intramolecular hydrogen bonds is placed in a hydrogenated solvent it may preferentially exchange deuterium for hydrogen. This preference is due to the difference between the vibrational zero-point energy for hydrogen and deuterium. It is found that the associated fractionation factor Φ is correlated with the strength of the intramolecular hydrogen bonds. This correlation has been used to determine the length of the H-bonds (donor-acceptor separation) in a diverse range of enzymes and has been argued to support the existence of short low-barrier H-bonds. Starting with a potential energy surface based on a simple diabatic state model for H-bonds we calculate Φ as a function of the proton donor-acceptor distance R. For numerical results, we use a parameterization of the model for symmetric O-H· · · O bonds [Chem. Phys. Lett. 535, 196 (2012)]. We consider the relative contributions of the O-H stretch vibration, O-H bend vibrations (both in plane and out of plane), tunneling splitting effects at finite temperature, and the secondary geometric isotope effect. We compare our total Φ as a function of R with NMR experimental results for enzymes, and in particular with an empirical parametrization Φ(R), used previously to determine bond lengths. *

We present a study of the abstraction of alkyl hydrogen atoms from the β and α positions of ethanol by the CN radical in solution using the Empirical Valence Bond (EVB) method. We have built separate 2 × 2 EVB models for the H and H reactions, where the atom transfer is parameterized using ab initio calculations. The intra- and intermolecular potentials of the reactant and product molecules were modelled with the General AMBER Force Field, with some modifications. We have carried out the dynamics in water and chloroform, which are solvents of contrasting polarity. We have computed the potential of mean force for both abstractions in each of the solvents. They are found to have a small and early barrier along the reaction coordinate with a large energy release. Analyzing the solvent structure around the reaction system, we have found two solvents to have little effect on either reaction. Simulating the dynamics from the transition state, we also study the fate of the energies in the HCN vibrational modes. The HCN molecule is born vibrationally hot in the CH stretch in both reactions and additionally in the HCN bends for the H abstraction reaction. In the early stage of the dynamics, we find that the CN stretch mode gains energy at the expense of the energy in CH stretch mode.