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.

Structure, Anharmonic Vibrational Frequencies, and Intensities of NNHNN(+).

A semiglobal potential energy surface (PES) and quartic force field (QFF) based on fitting high-level electronic structure energies are presented to describe the structures and spectroscopic properties of NNHNN(+). The equilibrium structure of NNHNN(+) is linear with the proton equidistant between the two nitrogen groups and thus of D(∞h) symmetry. Vibrational second-order perturbation theory (VPT2) calculations based on the QFF fails to describe the proton "rattle" motion, i.e., the antisymmetric proton stretch, due to the very flat nature of PES around the global minimum but performs properly for other modes with sharper potential wells. Vibrational self-consistent field/virtual state configuration interaction (VSCF/VCI) calculations using a version of MULTIMODE without angular momentum terms successfully describe this motion and predict the fundamental to be at 759 cm(-1). This is in good agreement with the value of 746 cm(-1) from a fixed-node diffusion Monte Carlo calculation and the experimental Ar-tagged result of 743 cm(-1). Other VSCF/VCI energies are in good agreement with other experimentally reported ones. Both double-harmonic intensity and rigorous MULTIMODE intensity calculations show the proton-transfer fundamental has strong intensity.

Ab Initio Quantum Approaches to the IR Spectroscopy of Water and Hydrates.

This Perspective highlights progress in ab initio quantum approaches to IR spectroscopy of water and hydrates. Here, "ab initio" refers to many-body potentials and dipole moment surfaces for flexible water and hydrates. Specifically, these are mathematical representations of two-body and three-body interactions based on permutationally invariant fitting of tens of thousands of ab initio electronic energies, a spectroscopically accurate one-body monomer potential, and four- and higher-body interactions described by the long-range interactions incorporated into, for example, the TTM3-F family of potentials. There are currently two such potentials of this type, denoted WHBB and MB-pol, which are being used in expanding applications. Here, the focus is on infrared spectroscopy, using the WHBB potential and dipole moment surface, with an embedded, local monomer quantum method to obtain vibrational energies and dipole transition moments. Comparisons are also made with the popular q-TIP4P/F potential. Brief mention is made of an application to small HCl-H2O clusters.

Communication: Spectroscopic consequences of proton delocalization in OCHCO⁺.

Even though quartic force fields (QFFs) and highly accurate coupled cluster computations describe the OCHCO(+) cation at equilibrium as a complex between carbon monoxide and the formyl cation, two notable and typical interstellar and atmospheric molecules, the prediction from the present study is that the equilibrium C(∞v) structure is less relevant to observables than the saddle-point D(∞h) structure. This is the conclusion from diffusion Monte Carlo and vibrational self-consistent field/virtual state configuration interaction calculations utilizing a semi-global potential energy surface. These calculations demonstrate that the proton "rattle" motion (ν6) has centrosymmetric delocalization of the proton over the D(∞h) barrier lying only 393.6 cm(-1) above the double-well OCHCO(+) C(∞v) minima. As a result, this molecule will likely appear D∞h, and the rotational spectrum will be significantly dimmer than the computed equilibrium 2.975 D center-of-mass dipole moment indicates. However, the proton transfer fundamental, determined to be at roughly 300 cm(-1), has a very strong intensity. This prediction as well as those of other fundamentals should provide useful guides for laboratory detection of this cation. Finally, it is shown that the two highest energy QFF-determined modes are actually in good agreement with their vibrational configuration interaction counterparts. These high-level quantum chemical methods provide novel insights into this fascinating and potentially common interstellar molecule.

Isolating the spectral signature of H3O(+) in the smallest droplet of dissociated HCl acid.

The centrally important role of acids in aqueous chemistry has stimulated the search for the smallest droplet of hydrochloric acid. Based on several independent quantum calculations, this appears to be the HCl(H2O)4 cluster, which dissociates into the so-called solvent ion pair (SIP), H3O(+)(H2O)3Cl(-). Experimental verification of this prediction via infra-red spectroscopy is a major challenge and despite several recent reports of this SIP, there remains uncertainty about these observations. In this report, we present a calculation of the IR spectrum of the SIP in a fashion that isolates the contribution from the signature hydronium ion, H3O(+). The computed spectrum indicates that the vibrational states of H3O(+) are highly mixed, resulting in dispersed spectral features between 1300 and 3000 cm(-1), with the region between 2100 and 2900 cm(-1) being especially rich. These predictions point out the complexity of the SIP spectrum and offer guidelines for experiment. The energies of the HCl stretch fundamentals for three minima of the undissociated HCl(H2O)4 cluster are also reported.

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.

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 trimer's 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 tetramer's anharmonic stretch frequencies. The optimized geometries are found to be in agreement with those of previous ab initio studies and the tetramer's anharmonic frequencies are computed within 11 cm(-1) of experiment.

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.

On-the-fly ab intito calculations of anharmonic vibrational frequencies: local-monomer theory and application to HCl clusters.

We present an on-the-fly quantum mechanical method to obtain anharmonic vibrational frequencies for molecular clusters. The basis for the method is the local-monomer model, a "divide and conquer" approach to theoretical spectroscopy, previously applied using full-dimensional surfaces [Y. Wang and J. M. Bowman, J. Chem. Phys. 134, 154510 (2011)]. The model consists of performing a local normal-mode analysis for each monomer in a cluster in the field of the surrounding monomers. Anharmonic vibrational frequencies are then determined for each monomer by numerically solving the Schrödinger equation in terms of the local coordinates using ab initio energies obtained directly. Residual monomer-monomer coupling is accounted for using the Hückel-coupling extension [Y. Wang and J. M. Bowman, J. Chem. Phys. 136, 144113 (2012)]. In addition to the direct local-monomer approach, we propose and demonstrate a composite ab initio technique to reduce computational costs for calculating the anharmonic frequencies of large clusters. This technique utilizes two ab initio methods, a lower level of theory to compute geometries and perform harmonic analyses and a subsequent higher level of theory to compute the energies used in the anharmonic frequency calculations. We demonstrate the on-the-fly approach on hydrogen chloride clusters ranging in size from the dimer to the hexamer. Comparisons of the theoretical frequencies are made to previous experiments. We find the method to be an effective and computationally efficient approach to compute anharmonic frequencies.

Communication: A new ab initio potential energy surface for HCl-H2O, diffusion Monte Carlo calculations of D0 and a delocalized zero-point wavefunction.

We report a global, full-dimensional, ab initio potential energy surface describing the HCl-H2O dimer. The potential is constructed from a permutationally invariant fit, using Morse-like variables, to over 44,000 CCSD(T)-F12b∕aug-cc-pVTZ energies. The surface describes the complex and dissociated monomers with a total RMS fitting error of 24 cm(-1). The normal modes of the minima, low-energy saddle point and separated monomers, the double minimum isomerization pathway and electronic dissociation energy are accurately described by the surface. Rigorous quantum mechanical diffusion Monte Carlo (DMC) calculations are performed to determine the zero-point energy and wavefunction of the complex and the separated fragments. The calculated zero-point energies together with a De value calculated from CCSD(T) with a complete basis set extrapolation gives a D0 value of 1348 ± 3 cm(-1), in good agreement with the recent experimentally reported value of 1334 ± 10 cm(-1) [B. E. Casterline, A. K. Mollner, L. C. Ch'ng, and H. Reisler, J. Phys. Chem. A 114, 9774 (2010)]. Examination of the DMC wavefunction allows for confident characterization of the zero-point geometry to be dominant at the C(2v) double-well saddle point and not the C(s) global minimum. Additional support for the delocalized zero-point geometry is given by numerical solutions to the 1D Schrödinger equation along the imaginary-frequency out-of-plane bending mode, where the zero-point energy is calculated to be 52 cm(-1) above the isomerization barrier. The D0 of the fully deuterated isotopologue is calculated to be 1476 ± 3 cm(-1), which we hope will stand as a benchmark for future experimental work.

(rac)-1,1'-binaphthyl-based simple receptors designed for fluorometric discrimination of maleic and fumaric acids.

(rac)-1,1'-Binaphthyl-based simple receptors 1 and 2 have been designed, synthesized and studied theoretically. The receptors utilize naphthyridine as the binding motifs for complexation of dicarboxylic acids in CHCl(3). The emission of the BINOL moiety was monitored experimentally to ascertain the selectivity and sensitivity of the receptors. Receptor 1 distinguishes maleic acid from isomeric fumaric acid by exhibiting different fluorescence behavior and demonstrates stronger binding in the excited state. Modulation of the binding sites of 1 leads to a new receptor structure 2, which was found to be less efficient in distinguishing maleic from fumaric acid, fluorometrically. Both 1 and 2 also recognize other hydroxy di- and tricarboxylic acids. The binding interactions were monitored by (1)H NMR, fluorescence and UV-vis spectroscopic methods. Structures of apo-hosts, guests and host-guest complexes were determined using force-field based conformational searching. Low energy ensembles were grouped into geometrically similar families, and low energy structures from each family were verified using B3LYP/6-31G*/PB-SCRF(CHCl(3)) calculations. The atomistic calculations provide insight into the differential dicarboxylic acid binding behavior of receptors 1 and 2.