Hydrogen Bonding and OH-Stretch Spectroscopy in Water: Hexamer (Cage), Liquid Surface, Liquid, and Ice

Towards a unified description of the hydrogen bond network of liquid water: a dynamics based approach

The definition of a hydrogen bond (H-bond) is intimately related to the topological and dynamic properties of the hydrogen bond network within liquid water. The development of a universal H-bond definition for water is an active area of research as it would remove many ambiguities in the network properties that derive from the fixed definition employed to assign whether a water dimer is hydrogen bonded. This work investigates the impact that an electronic-structure based definition, an energetic, and a geometric definition of the H-bond has upon both topological and dynamic network behavior of simulated water. In each definition, the use of a cutoff (either geometric or energetic) to assign the presence of a H-bond leads to the formation of transiently bonded or broken dimers, which have been quantified within the simulation data. The relative concentration of transient species, and their duration, results in two of the three definitions sharing similarities in either topological or dynamic features (H-bond distribution, H-bond lifetime, etc.), however no two definitions exhibit similar behavior for both classes of network properties. In fact, two networks with similar local network topology (as indicated by similar average H-bonds) can have dramatically different global network topology (as indicated by the defect state distributions) and altered H-bond lifetimes. A dynamics based correction scheme is then used to remove artificially transient H-bonds and to repair artificially broken bonds within the network such that the corrected network exhibits the same structural and dynamic properties for two H-bond definitions (the properties of the third definition being significantly improved). The algorithm described represents a significant step forward in the development of a unified hydrogen bond network whose properties are independent of the original hydrogen bond definition that is employed.

Insights into hydrogen bonding via ice interfaces and isolated water

Water in a confined environment has a combination of fewer available configurations and restricted mobility. Both affect the spectroscopic signature. In this work, the spectroscopic signature of water in confined environments is discussed in the context of competing models for condensed water: (1) as a system of intramolecular coupled molecules or (2) as a network with intermolecular dipole-dipole coupled O-H stretches. Two distinct environments are used: the confined asymmetric environment at the ice surface and the near-isolated environment of water in an infrared transparent matrix. Both the spectroscopy and the environment are described followed by a perspective discussion of implications for the two competing models. Despite being a small molecule, water is relatively complex; perhaps not surprisingly the results support a model that blends inter- and intramolecular coupling. The frequency, and therefore the hydrogen-bond strength, appears to be a function of donor-acceptor interaction and of longer-range dipole-dipole alignment in the hydrogen-bonded network. The O-H dipole direction depends on the local environment and reflects intramolecular O-H stretch coupling.

Isomer-Specific Spectroscopy of Benzene-(H2O)n, n = 6,7: Benzene's Role in Reshaping Water's Three-Dimensional Networks

The water hexamer and heptamer are the smallest sized water clusters that support three-dimensional hydrogen-bonded networks, with several competing structures that could be altered by interactions with a solute. Using infrared-ultraviolet double resonance spectroscopy, we record isomer-specific OH stretch infrared spectra of gas-phase benzene-(H2O)(6,7) clusters that demonstrate benzene's surprising role in reshaping (H2O)(6,7). The single observed isomer of benzene-(H2O)6 incorporates an inverted book structure rather than the cage or prism. The main conformer of benzene-(H2O)7 is an inserted-cubic structure in which benzene replaces one water molecule in the S4-symmetry cube of the water octamer, inserting itself into the water cluster by engaging as a π H-bond acceptor with one water and via C-H···O donor interactions with two others. The corresponding D(2d)-symmetry inserted-cube structure is not observed, consistent with the calculated energetic preference for the S4 over the D(2d) inserted cube. A reduced-dimension model that incorporates stretch-bend Fermi resonance accounts for the spectra in detail and sheds light on the hydrogen-bonding networks themselves and on the perturbations imposed on them by benzene.

Efficient anharmonic vibrational spectroscopy for large molecules using local-mode coordinates

This article presents a general computational approach for efficient simulations of anharmonic vibrational spectra in chemical systems. An automated local-mode vibrational approach is presented, which borrows techniques from localized molecular orbitals in electronic structure theory. This approach generates spatially localized vibrational modes, in contrast to the delocalization exhibited by canonical normal modes. The method is rigorously tested across a series of chemical systems, ranging from small molecules to large water clusters and a protonated dipeptide. It is interfaced with exact, grid-based approaches, as well as vibrational self-consistent field methods. Most significantly, this new set of reference coordinates exhibits a well-behaved spatial decay of mode couplings, which allows for a systematic, a priori truncation of mode couplings and increased computational efficiency. Convergence can typically be reached by including modes within only about 4 Å. The local nature of this truncation suggests particular promise for the ab initio simulation of anharmonic vibrational motion in large systems, where connection to experimental spectra is currently most challenging.

Resonant vibrational energy transfer in ice Ih

Fascinating anisotropy decay experiments have recently been performed on H2O ice Ih by Timmer and Bakker [R. L. A. Timmer, and H. J. Bakker, J. Phys. Chem. A 114, 4148 (2010)]. The very fast decay (on the order of 100 fs) is indicative of resonant energy transfer between OH stretches on different molecules. Isotope dilution experiments with deuterium show a dramatic dependence on the hydrogen mole fraction, which confirms the energy transfer picture. Timmer and Bakker have interpreted the experiments with a Förster incoherent hopping model, finding that energy transfer within the first solvation shell dominates the relaxation process. We have developed a microscopic theory of vibrational spectroscopy of water and ice, and herein we use this theory to calculate the anisotropy decay in ice as a function of hydrogen mole fraction. We obtain very good agreement with experiment. Interpretation of our results shows that four nearest-neighbor acceptors dominate the energy transfer, and that while the incoherent hopping picture is qualitatively correct, vibrational energy transport is partially coherent on the relevant timescale.

Determination of spectroscopic band shapes using second derivatives, part I: theory

The molecular spectra of water, aqueous solutions, hydrogen-bonded systems, and others have massive bands that contain many overlapping components. To decipher the spectra for molecular interpretation, it is necessary to separate these. Several attempts to do this have been made without clear success. To surmount some of the difficulties, we present a novel method, which consists of quantitatively evaluating the spectral band second-derivative profiles. This aids in the determination of the original band profiles: Gaussian, Lorentzian (Cauchy), and Gauss-Lorentz products. Then the number of components in a massive absorption, their shapes, and their positions can be determined. We tested the usefulness of the method in the visible region using calibration standards: a light emitting diode emission spectrum and a holmium chloride (HoCl2) solution. To verify its utility in the infrared region, we used liquid propanol, liquid acetonitrile, and aqueous acetone.

Infrared and Raman Spectroscopy of Liquid Water through "First-Principles" Many-Body Molecular Dynamics

Vibrational spectroscopy is a powerful technique to probe the structure and dynamics of water. However, deriving an unambiguous molecular-level interpretation of the experimental spectral features remains a challenge due to the complexity of the underlying hydrogen-bonding network. In this contribution, we present an integrated theoretical and computational framework (named many-body molecular dynamics or MB-MD) that, by systematically removing uncertainties associated with existing approaches, enables a rigorous modeling of vibrational spectra of water from quantum dynamical simulations. Specifically, we extend approaches used to model the many-body expansion of interaction energies to develop many-body representations of the dipole moment and polarizability of water. The combination of these "first-principles" representations with centroid molecular dynamics simulations enables the simulation of infrared and Raman spectra of liquid water under ambient conditions that, without relying on any ad hoc parameters, are in good agreement with the corresponding experimental results. Importantly, since the many-body energy, dipole, and polarizability surfaces employed in the simulations are derived independently from accurate fits to correlated electronic structure data, MB-MD allows for a systematic analysis of the calculated spectra in terms of both electronic and dynamical contributions. The present analysis suggests that, while MB-MD correctly reproduces both the shifts and the shapes of the main spectroscopic features, an improved description of quantum dynamical effects possibly combined with a dissociable water potential may be necessary for a quantitative representation of the OH stretch band.

Ions increase strength of hydrogen bond in water

Knowledge of water-water potential is important for an accurate description of water. Potential between two molecules depends upon the distance, relative orientation of each molecule and local environment. In simulation, water-water hydrogen bonds are handled by point-charge water potentials and by polarizable models. These models produce good results for bulk water being parameterized for such environment. Water around surfaces and in channels, however is different from bulk water. Using quantum-mechanical methods, hydrogen bond strength was calculated in the vicinity of different monoions. A simple empirical relationship was discovered between the maximum hydrogen bond and the electric field produced by ion.

Local-monomer calculations of the intramolecular IR spectra of the cage and prism isomers of HOD(D2O)5 and HOD and D2O ice Ih

Dilute mixtures of HOD in pure H2O and D2O ices and liquid have been used by experimentalists to focus on the spectrum and vibrational dynamics of the local OH and OD stretches and bend of HOD in these complex and highly heterogeneous environments. The hexamer version of the mixture is HOD(D2O)5. The cage isomer of this cluster was recently studied and analyzed theoretically using local-mode calculations of the IR spectrum by Skinner and co-workers. This and the further possibility of experimental investigation of this cluster have stimulated us to study HOD(D2O)5 using the three-mode, local-monomer model, with the ab initio WHBB dipole moment and potential energy surfaces. Both the cage and prism isomers of this cluster are considered. In addition to providing additional insight into the HOD portion of the spectrum, the spectral signatures of each D2O are also presented in the range of 1000-4000 cm(-1). The OH stretch bands of both the prism and cage isotopomers exhibit rich structures in the range of 3100-3700 cm(-1) that are indicative of the position of the HOD in these hexamers. A preliminary investigation of the site preference of the HOD is also reported for both cage and prism HOD(D2O)5 using an harmonic zero-point energy analysis of the entire cluster. This indicates that the energies of free-OH sites are lower than the ones of H-bonded OH sites. Finally, following our earlier work on the IR spectra of H2O ice models, we present IR spectra of pure D2O and HOD.

Vibrational spectroscopy of water in hydrated lipid multi-bilayers. III. Water clustering and vibrational energy transfer

Water clustering and connectivity around lipid bilayers strongly influences the properties of membranes and is important for functions such as proton and ion transport. Vibrational anisotropic pump-probe spectroscopy is a powerful tool for understanding such clustering, as the measured anisotropy depends upon the time-scale and degree of intra- and intermolecular vibrational energy transfer. In this article, we use molecular dynamics simulations and theoretical vibrational spectroscopy to help interpret recent experimental measurements of the anisotropy of water in lipid multi-bilayers as a function of both lipid hydration level and isotopic substitution. Our calculations are in satisfactory agreement with the experiments of Piatkowski, Heij, and Bakker, and from our simulations we can directly probe water clustering and connectivity. We find that at low hydration levels, many water molecules are in fact isolated, although up to 70% of hydration water forms small water clusters or chains. At intermediate hydration levels, water forms a wide range of cluster sizes, while at higher hydration levels, the majority of water molecules are part of a large, percolating water cluster. Therefore, the size, number, and nature of water clusters are strongly dependent on lipid hydration level, and the measured anisotropy reflects this through its dependence on intermolecular energy transfer.

Effects of excluded volume and correlated molecular orientations on Förster resonance energy transfer in liquid water

Förster theory for the survival probability of excited chromophores is generalized to include the effects of excluded volume and orientation correlation in the molecular distribution. An analytical expression for survival probability was derived and written in terms of a few simple elementary functions. Because of the excluded volume, the survival probability exhibits exponential decay at early times and stretched exponential decay at later times. Experimental schemes to determine the size of the molecular excluded volume are suggested. With the present generalization of theory, we analyzed vibrational resonance energy transfer kinetics in neat water. Excluded volume effects prove to be important and slow down the kinetics at early times. The majority of intermolecular resonance energy transfer was found to occur with exponential kinetics, as opposed to the stretched exponential behavior predicted by Förster theory. Quantum yields of intra-molecular vibrational relaxation, intra-, and intermolecular energy transfer were calculated to be 0.413, 0.167, and 0.420, respectively.

Instantaneous normal mode analysis for intermolecular and intramolecular vibrations of water from atomic point of view

By exploiting the instantaneous normal mode (INM) analysis for models of flexible molecules, we investigate intermolecular and intramolecular vibrations of water from the atomic point of view. With two flexible SPC/E models, our investigations include three aspects about their INM spectra, which are separated into the unstable, intermolecular, bending, and stretching bands. First, the O- and H-atom contributions in the four INM bands are calculated and their stable INM spectra are compared with the power spectra of the atomic velocity autocorrelation functions. The unstable and intermolecular bands of the flexible models are also compared with those of the SPC/E model of rigid molecules. Second, we formulate the inverse participation ratio (IPR) of the INMs, respectively, for the O- and H-atom and molecule. With the IPRs, the numbers of the three species participated in the INMs are estimated so that the localization characters of the INMs in each band are studied. Further, by the ratio of the IPR of the H atom to that of the O atom, we explore the number of involved OH bond per molecule participated in the INMs. Third, by classifying simulated molecules into subensembles according to the geometry of their local environments or their H-bond configurations, we examine the local-structure effects on the bending and stretching INM bands. All of our results are verified to be insensible to the definition of H-bond. Our conclusions about the intermolecular and intramolecular vibrations in water are given.

Slow hydrogen-bond switching dynamics at the water surface revealed by theoretical two-dimensional sum-frequency spectroscopy

Using our newly developed explicit three-body (E3B) water model, we simulate the surface of liquid water. We find that the timescale for hydrogen-bond switching dynamics at the surface is about three times slower than that in the bulk. In contrast, with this model rotational dynamics are slightly faster at the surface than in the bulk. We consider vibrational two-dimensional (2D) sum-frequency generation (2DSFG) spectroscopy as a technique for observing hydrogen-bond rearrangement dynamics at the water surface. We calculate the nonlinear susceptibility for this spectroscopy for two different polarization conditions, and in each case we see the appearance of cross-peaks on the timescale of a few picoseconds, signaling hydrogen-bond rearrangement on this timescale. We thus conclude that this 2D spectroscopy will be an excellent experimental technique for observing slow hydrogen-bond switching dynamics at the water surface.