John S. Mancini

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.

Effects of Zero-Point Delocalization on the Vibrational Frequencies of Mixed HCl and Water Clusters.

We demonstrate the significant effect that large-amplitude zero-point vibrational motion can have on the high-frequency fundamental vibrations of molecular clusters, specifically small (HCl)n-(H2O)m clusters. Calculations were conducted on a many-body potential, constructed from a mix of new and previously reported semiempirical and high-level ab initio potentials. Diffusion Monte Carlo simulations were performed to determine ground-state wave functions. Visualization of these wave functions indicates that the clusters exhibit delocalized ground states spanning multiple stationary point geometries. The ground states are best characterized by planar ring configurations, despite the clusters taking nonplanar configurations at their global minima. Vibrational calculations were performed at the global minima and the Diffusion Monte Carlo predicted configurations and also using an approach that spans multiple stationary points along a rectilinear normal-mode reaction path. Significantly better agreement was observed between the calculated vibrational frequencies and experimental peak positions when the delocalized ground state was accounted for.

Energetics and Predissociation Dynamics of Small Water, HCl, and Mixed HCl–Water Clusters

This Review summarizes recent research on vibrational predissociation (VP) of hydrogen-bonded clusters. Specifically, the focus is on breaking of hydrogen bonds following excitation of an intramolecular vibration of the cluster. VP of the water dimer and trimer, HCl clusters, and mixed HCl-water clusters are the major topics, but related work on hydrogen halide dimers and trimers, ammonia clusters, and mixed dimers with polyatomic units are reviewed for completion and comparison. The theoretical focus is on generating accurate potential energy surfaces (PESs) that can be used in detailed dynamical calculations, mainly using the quasiclassical trajectory approach. These PESs have to extend from the region describing large amplitude motion around the minimum to regions where fragments are formed. The experimental methodology exploits velocity map imaging to generate pair-correlated product translational energy distributions from which accurate bond dissociation energies of dimers and trimers and energy disposal in fragments are obtained. The excellent agreement between theory and experiment on bond dissociation energies, energy disposal in fragments, and the contributions of cooperativity demonstrates that it is now possible, with state-of-the-art experimental and theoretical methods, to make accurate predictions about dynamical and energetic properties of dissociating clusters.

Experiment and Theory Elucidate the Multichannel Predissociation Dynamics of the HCl Trimer: Breaking Up Is Hard To Do

The breaking of hydrogen bonds in molecular systems has profound effects on liquids, e.g., water, biomolecules, e.g., DNA, etc., and so it is no exaggeration to assert the importance of these bonds to living systems. However, despite years of extensive research on hydrogen bonds, many of the details of how these bonds break and the corresponding energy redistribution processes remain poorly understood. Here we report extensive experimental and theoretical insights into the breakup of two or three hydrogen bonds in the dissociation of a paradigm system of a hydrogen-bonded network, the ring HCl trimer. Experimental state-to-state vibrational predissociation dynamics of the trimer following vibrational excitation were studied by using velocity map imaging and resonance-enhanced multiphoton ionization, providing dissociation energies and product state distributions for the trimers breakup into three separate monomers or into dimer + monomer. Accompanying the experiments are high-level calculations using diffusion Monte Carlo and quasiclassical simulations, whose results validate the experimental ones and further elucidate energy distributions in the products. The calculations make use of a new, highly accurate potential energy surface. Simulations indicate that the dissociation mechanism requires the excitation to first relax into low-frequency motions of the trimer, resulting in the breaking of a single hydrogen bond. This allows the system to explore a critical van der Waals minimum region from which dissociation occurs readily to monomer + dimer.

A new many-body potential energy surface for HCl clusters and its application to anharmonic spectroscopy and vibration-vibration energy transfer in the HCl trimer.

The hydrogen bond has been studied by chemists for nearly a century. Interest in this ubiquitous bond has led to several prototypical systems emerging to studying its behavior. Hydrogen chloride clusters stand as one such example. We present here a new many-body potential energy surface for (HCl)n constructed from one-, two-, and three-body interactions. The surface is constructed from previous highly accurate, semiempirical monomer and dimer surfaces, and a new high-level ab initio permutationally invariant full-dimensional three-body potential. The new three-body potential is based on fitting roughly 52,000 three-body energies computed using coupled cluster with single, doubles, perturbative triples, and explicit correlation and the augmented correlation consistent double-ζ basis set. The first application, described here, is to the ring HCl trimer, for which the many-body representation is exact. The new potential describes all known stationary points of the trimer as well its dissociation to either three monomers or a monomer and a dimer. The anharmonic vibrational energies are computed for the three H-Cl stretches, using explicit three-mode coupling calculations and local-monomer calculations with Hückel-type coupling. Both methods produce frequencies within 5 cm(-1) of experiment. A wavepacket calculation based on the Hückel model and full-dimensional classical calculation are performed to study the monomer H-Cl stretch vibration-vibration transfer process in the ring HCl trimer. Somewhat surprisingly, the results of the quantum and classical calculations are virtually identical, both exhibiting coherent beating of the excitation between the three monomers. Finally, this representation of the potential is used to study properties of larger clusters, namely to compute optimized geometries of the tetramer, pentamer, and hexamer and to perform explicit four-mode coupling calculations of the tetramers anharmonic stretch frequencies. The optimized geometries are found to be in agreement with those of previous ab initio studies and the tetramers anharmonic frequencies are computed within 11 cm(-1) of experiment.