Imaging H2O Photofragments in the Predissociation of the HCl−H2O Hydrogen-Bonded Dimer

Room temperature detection of the (H2)2 dimer

The hydrogen dimer, (H2)2, is among the most weakly bound van der Waals complexes and a prototype species for first principles ab initio studies. The detection of the (H2)2 infrared absorption spectrum was reported more than sixty years ago at a temperature of 20 K. Due to the sharp decrease of the (H2)2 abundance with temperature, detection at room temperature was generally considered hardly achievable. Here we report the first room temperature detection of partly resolved rotational structures of (H2)2 by cavity ring down spectroscopy at sub-atmospheric pressures, in the region of the first overtone band of H2 near 1.2 μm. The quantitative analysis of the absorption features observed around ten allowed or forbidden transition frequencies of the monomer provides insight on the structure of this elusive species and a benchmark for future theoretical studies.

Mechanism of ionic dissociation of HCl in the smallest water clusters

The dissociation of strong acids into water is a fundamental process in chemistry and biology. Determining the minimum number of water molecules that can result in an ionic dissociation of hydrochloric acid (HCl → H+ + Cl-) remains a challenging subject. In this study, the reactions of H2O with HCl(H2O)n-1 (HCl-H2O cluster), i.e., HCl(H2O)n-1 + H2O (n = 3-7), were investigated by using the direct ab initio molecular dynamics (AIMD) method. Direct AIMD calculations were performed to set the collision energy of H2O to zero for all trajectories. For n = 3, no reaction occurred. In contrast, HCl dissociated to H+ + Cl- at n = 4, forming a contact ion pair (cIP) and solvent-separated ion pair (ssIP) as products. The reactions were expressed as HCl(H2O)3 + H2O → H3O+(H2O)2Cl- (ssIP), and HCl(H2O)3 + H2O → H3O+(Cl-)(H2O)2 (cIP). The ion pair (IP) products were dependent on the collision site of H2O relative to HCl(H2O)3. For n = 5-7, both IPs were formed through the reaction between H2O and HCl(H2O)n-1 (n = 5-7). The reaction between HCl and (H2O)4 (HCl + (H2O)4 → HCl(H2O)4) was non-reactive in IP formation. The reaction mechanism was discussed based on the theoretical results.

Femtochemistry of bimolecular reactions from weakly bound complexes: computational study of the H + H'OD → H'OH + D or HOD + H' exchange reactions

A full-dimensional wavepacket propagation describing the bimolecular exchange reactions H + H'OD → H'OH + D or HOD + H' initiated by photolysis of HCl in the hydrogen-bound complex (HCl)⋯(HOD) is reported. The dynamics of this reaction is carried out with the MCTDH method on an ab initio potential energy surface (PES) of H₃O and the initial state is derived from the ground state wavefunction of the complex obtained by relaxation on its own electronic ground state ab initio PES. The description of the system makes use of polyspherical coordinates parametrizing a set of Radau and Jacobi vectors. The calculated energy- and time-resolved reaction probabilities show, owing to the large collision energies at play stemming from the (almost full) photolysis of HCl, that the repulsion between oxygen in the H'OD molecule and the incoming hydrogen atom is the main feature of the collision and leads to non-reactive scattering. No abstraction reaction products are observed. However, both exchange processes are still observable, with a preference in O-H' bond dissociation over that of O-D. The selectivity is reversed upon vibrational pre-excitation of the O-D stretching mode in the H'OD molecule. It is shown that, after the collision, the hydrogen atom of HCl does most likely not encounter the almost stationary chlorine atom again but we also consider the limit case where the H atom is forced to collide multiple times against H'OD as a result of being pushed back by the Cl atom.

Vibrational predissociation of the phenol-water dimer: a view from the water

The vibrational predissociation (VP) dynamics of the phenol-water (PhOH-H₂O) dimer were studied by detecting H₂O fragments and using velocity map imaging (VMI) to infer the internal energy distributions of PhOH cofragments, pair-correlated with selected rotational levels of the H₂O fragments. Following infrared (IR) laser excitation of the hydrogen-bonded OH stretch fundamental of PhOH (Pathway 1) or the asymmetric OH stretch localized on H₂O (Pathway 2), dissociation to H₂O + PhOH was observed. H₂O fragments were monitored state-selectively by using 2+1 Resonance-Enhanced Multiphoton Ionization (REMPI) combined with time-of-flight mass spectrometry (TOF-MS). VMI of H₂O in selected rotational levels was used to derive center-of-mass (c.m.) translational energy (E_T) distributions. The pair-correlated internal energy distributions of the PhOH cofragments derived via Pathway 1 were well described by a statistical prior distribution. On the other hand, the corresponding distributions obtained via Pathway 2 show a propensity to populate higher-energy rovibrational levels of PhOH than expected from a statistical distribution and agree better with an energy-gap model. The REMPI spectra of the H₂O fragments from both pathways could be fit by Boltzmann plots truncated at the maximum allowed energy, with a higher temperature for Pathway 2 than that for Pathway 1. We conclude that the VP dynamics depends on the OH stretch level initially excited.

Relaxation of the H2O Overtone Bending Vibration in the Water Dimer···Hydroxyl Radical Complex

The relaxation mechanism of the overtone bending vibration in the collision of the water dimer with the vibrationally excited hydroxyl radical is studied by use of trajectory procedures. The transfer of the OH(v = 1) energy to the dimer stretches is followed by a near-resonant first overtone transition to the donor monomer. Nearly a quarter of the trajectories undergo a complex-mode collision, forming the (H₂O)₂···OH complex bound by a hydrogen bond with the lifetime ranging from a subpicosecond scale to >100 ps. The overtone vibration relaxes to the ground state, transferring approximately half of its energy to the dimer hydrogen-bonding (H₂O···H₂O) and the remaining half to the complex hydrogen-bonding (H₂O)₂···OH, via near-resonant pathways, each consisting of a series of intermolecular low-frequency vibrations.

Predissociation dynamics of the HCl-(H2O)3 tetramer: An experimental and theoretical investigation

The cyclic HCl-(H₂O)₃ tetramer is the largest observed neutral HCl-(H₂O)_n cluster. The vibrational predissociation of HCl-(H₂O)₃ is investigated by theory, quasiclassical trajectory (QCT) calculations, and experiment, following the infrared excitation of the hydrogen-bonded OH-stretch fundamental. The energetically possible dissociation pathways are HCl + (H₂O)₃ (Pathway 1) and H₂O + HCl-(H₂O)₂ (Pathway 2). The HCl and H₂O monomer fragments are observed by 2 + 1 resonance enhanced multiphoton ionization combined with time-of-flight mass spectrometry, and their rotational energy distributions are inferred and compared to the theoretical results. Velocity map images of the monomer fragments in selected rotational levels are used for each pathway to obtain pair-correlated speed distributions. The fragment speed distributions obtained by experiment and QCT calculations are broad and structureless, encompassing the entire range of allowed speeds for each pathway. Bond dissociation energies, D₀, are estimated experimentally from the endpoints of the speed distributions: 2100 ± 300 cm^-1 and 2400 ± 100 cm^-1 for Pathway 1 and Pathway 2, respectively. These values are lower but in the same order as the corresponding calculated D₀: 2426 ± 23 cm^-1 and 2826 ± 19 cm^-1. The differences are attributed to contributions from vibrational hot bands of the clusters that appear in the high-speed tails of the experimental pair-correlated distributions. Satisfactory agreement between theory and experiment is achieved when comparing the monomer fragments' rotational energies, the shapes of the speed distributions, and the average fragment speeds and center-of-mass translational energies. Insights into the dissociation mechanism and lifetime are gained from QCT calculations performed on a previously reported many-body potential energy surface. It is concluded that the dissociation lifetime is on the order of 10 ps and that the final trimer products are in their lowest energy cyclic forms.

Energy Transfer to the Hydrogen Bond in the (H2O)2 + H2O Collision

Trajectory procedures are used to study the collision between the vibrationally excited H₂O and the ground-state (H₂O)₂ with particular reference to energy transfer to the hydrogen bond through the inter- and intramolecular pathways. In nearly 98% of the trajectories, energy transfer processes occur on a subpicosecond scale (≤0.7 ps). The H₂O transfers approximately three-quarters of its excitation energy to the OH stretches of the dimer. The first step of the intramolecular pathway in the dimer involves a near-resonant first overtone transition from the OH stretch to the bending mode. The energy transfer probability in the presence of the 1:2 resonance is 0.61 at 300 K. The bending mode then redistributes its energy to low-frequency intermolecular vibrations in a series of small excitation steps, with the pathway which results in the hydrogen-bonding modes gaining most of the available energy. The hydrogen bonding in ∼50% of the trajectories ruptures on vibrational excitation, leaving one quantum in the bend of the monomer fragment. In a small fraction of trajectories, the duration of collision is longer than 1 ps, during which the dimer and H₂O form a short-lived complex through a secondary hydrogen bond, which undergoes large amplitude oscillations.

α-Cyclodextrin: How Effectively Can Its Hydrophobic Cavity Be Hydrated?

Cyclodextrins (CDs) are among the most widely used native host systems with ability to form inclusion complexes with various molecular objects. This ability is so strong that the "hydrophobic" CD cavity never remains empty, even in the guest-free state it is filled with water molecules. However, no consensus has been reached concerning both the total number of hydrating water molecules and their preferred binding location in the CDs. Several outstanding questions regarding the CD hydration still wait to be answered: (1) Which spots of the CD cavity ("hot spots") have the highest affinity for the guest water molecules? (2) How stable are water clusters inside the cavity? (3) Which mode of water binding, sequential or bulk, is thermodynamically more favored? (4) What is the upper limit of the number of water molecules bound inside the host cavity? (5) What factors do control the CD hydration process? Here, using αCD as a typical representative of the cyclodextrin family, we endeavor to answer these questions by combining experimental measurements (differential scanning calorimetry and thermogravimetry) with theoretical (DFT) calculations. Enthalpies of the αCD hydrate formation are evaluated and the role of different factors, such as the number and mode of binding (sequential vs bulk) of water molecules, type of hydrogen bonds established (water-water vs water-αCD), and the dielectric properties of the medium, on the complexation process is assessed. The results obtained shed light on the intimate mechanism of water binding to αCD and disclose the key factors governing the process.

Vibrational Predissociation of the HCl-(H2O)3 Tetramer

The vibrational predissociation of the HCl-(H₂O)₃ tetramer, the largest HCl-(H₂O)_n cluster for which HCl is not predicted to be ionized, is reported. This work focuses on the predissociation pathway giving rise to H₂O + HCl-(H₂O)₂ following IR laser excitation of the H-bonded OH stretch fundamental. H₂O fragments are monitored state selectively by 2 + 1 resonance-enhanced multiphoton ionization (REMPI) combined with time-of-flight mass spectrometry (TOF-MS). Velocity map images of H₂O in selected rotational levels are used to determine translational energy distributions from which the internal energy distributions in the pair-correlated cofragments are derived. From the maximum translational energy release, the bond dissociation energy, D₀ = 2400 ± 100 cm^-1, is determined for the investigated channel. The energy distributions in the fragments are broad, encompassing the entire range of allowed states. The importance of cooperative (nonpairwise) interactions is discussed.

Experimental and theoretical investigations of energy transfer and hydrogen-bond breaking in small water and HCl clusters

Water is one of the most pervasive molecules on earth and other planetary bodies; it is the molecule that is searched for as the presumptive precursor to extraterrestrial life. It is also the paradigm substance illustrating ubiquitous hydrogen bonding (H-bonding) in the gas phase, liquids, crystals, and amorphous solids. Moreover, H-bonding with other molecules and between different molecules is of the utmost importance in chemistry and biology. It is no wonder, then, that for nearly a century theoreticians and experimentalists have tried to understand all aspects of H-bonding and its influence on reactivity. It is somewhat surprising, therefore, that several fundamental aspects of H-bonding that are particularly important for benchmarking theoretical models have remained unexplored experimentally. For example, even the binding strength between two gas-phase water molecules has never been determined with sufficient accuracy for comparison with high-level electronic structure calculations. Likewise, the effect of cooperativity (nonadditivity) in small H-bonded networks is not known with sufficient accuracy. An even greater challenge for both theory and experiment is the description of the dissociation dynamics of H-bonded small clusters upon acquiring vibrational excitation. This is because of the long lifetimes of many clusters, which requires running classical trajectories for many nanoseconds to achieve dissociation. In this Account, we describe recent progress and ongoing research that demonstrates how the combined and complementary efforts of theory and experiment are enlisted to determine bond dissociation energies (D0) of small dimers and cyclic trimers of water and HCl with unprecedented accuracy, describe dissociation dynamics, and assess the effects of cooperativity. The experimental techniques rely on IR excitation of H-bonded X-H stretch vibrations, measuring velocity distributions of fragments in specific rovibrational states, and determining product state distributions at the pair-correlation level. The theoretical methods are based on high-level ab initio potential energy surfaces used in quantum and classical dynamical calculations. We achieve excellent agreement on D0 between theory and experiments for all of the clusters that we have compared, as well as for cooperativity in ring trimers of water and HCl. We also show that both the long-range and the repulsive parts of the potential must be involved in bond breaking. We explain why H-bonds are so resilient and hard to break, and we propose that a common motif in the breaking of cyclic trimers is the opening of the ring following transfer of one quantum of stretch excitation to form open-chain structures that are weakly bound. However, it still takes many vibrational periods to release one monomer fragment from the open-chain structures. Our success with water and HCl dimers and trimers led us to embark on a more ambitious project: studies of mixed water and HCl small clusters. These clusters eventually lead to ionization of HCl and serve as prototypes of acid dissociation in water. Measurements and calculations of such ionizations are yet to be achieved, and we are now characterizing these systems by adding monomers one at a time. We describe our completed work on the HCl-H2O dimer and mention our recent theoretical results on larger mixed clusters.

Experimental and theoretical investigations of the dissociation energy (D0) and dynamics of the water trimer, (H2O)3

We report a joint experimental-theoretical study of the predissociation dynamics of the water trimer following excitation of the hydrogen bonded OH-stretch fundamental. The bond dissociation energy (D0) for the (H2O)3 → H2O + (H2O)2 dissociation channel is determined from fitting the speed distributions of selected rovibrational states of the water monomer fragment using velocity map imaging. The experimental value, D0 = 2650 ± 150 cm(-1), is in good agreement with the previously determined theoretical value, 2726 ± 30 cm(-1), obtained using an ab initio full-dimensional potential energy surface (PES) together with Diffusion Monte Carlo calculations [ Wang ; Bowman . J. Chem. Phys. 2011 , 135 , 131101 ]. Comparing this value to D0 of the dimer places the contribution of nonpairwise additivity to the hydrogen bonding at 450-500 cm(-1). Quasiclassical trajectory (QCT) calculations using this PES help elucidate the reaction mechanism. The trajectories show that most often one hydrogen bond breaks first, followed by breaking and re-forming of hydrogen bonds (often with different hydrogen bonds breaking) until, after many picoseconds, a water monomer is finally released. The translational energy distributions calculated by QCT for selected rotational levels of the monomer fragment agree with the experimental observations. The product translational and rotational energy distributions calculated by QCT also agree with statistical predictions. The availability of low-lying intermolecular vibrational levels in the dimer fragment is likely to facilitate energy transfer before dissociation occurs, leading to statistical-like product state distributions.

Communication: determination of the bond dissociation energy (D0) of the water dimer, (H2O)2, by velocity map imaging

The bond dissociation energy (D(0)) of the water dimer is determined by using state-to-state vibrational predissociation measurements following excitation of the bound OH stretch fundamental of the donor unit of the dimer. Velocity map imaging and resonance-enhanced multiphoton ionization (REMPI) are used to determine pair-correlated product velocity and translational energy distributions. H(2)O fragments are detected in the ground vibrational (000) and the first excited bending (010) states by 2 + 1 REMPI via the C̃ (1)B(1) (000) ← X̃ (1)A(1) (000 and 010) transitions. The fragments' velocity and center-of-mass translational energy distributions are determined from images of selected rovibrational levels of H(2)O. An accurate value for D(0) is obtained by fitting both the structure in the images and the maximum velocity of the fragments. This value, D(0) = 1105 ± 10 cm(-1) (13.2 ± 0.12 kJ/mol), is in excellent agreement with the recent theoretical value of D(0) = 1103 ± 4 cm(-1) (13.2 ± 0.05 kJ∕mol) suggested as a benchmark by Shank et al. [J. Chem. Phys. 130, 144314 (2009)].