The water hexamer: three-body interactions, structures, energetics, and OH-stretch spectroscopy at finite temperature

Water-Hydrocarbon Interactions in Anionic Pyrene Monohydrate

Interactions between water and polycyclic aromatic hydrocarbons are essential in many aspects of chemistry, from interstellar and atmospheric processes to interfacial hydrophobicity and wetting phenomena. Despite their growing importance, the intermolecular potentials of the water-hydrocarbon interactions are underdeveloped compared to the water-water potentials, and there are similarly few experimental probes that are sensitive to the details of the water-hydrocarbon potential. We present a combined experimental and computational study of anionic pyrene monohydrate, one of the simplest water/hydrocarbon clusters. The action spectrum in the OH region of the mass-selected cluster ion provides a rigorous benchmark for intermolecular potentials and computational methodologies. We identify missing intermolecular interactions and shortcomings in conventional dynamics calculations by comparing experimental data to density functional theory and classical molecular dynamics calculations. Kinetic trapping is prevalent, even for one water molecule and one pyrene molecule, leading to slow equilibration in conventional molecular dynamics calculations, even on nanosecond time scales and at low temperatures (50 K). At constant energy, temperature fluctuations for the pair of molecules are substantial. Immersing the system in a bath of soft spheres and employing parallel tempering alleviates kinetic trapping and dampens temperature fluctuations, bringing the system closer to the thermodynamic limit. With such augmented sampling, a simple, flexible water model reproduces the line width and the asymmetric broadening of the symmetric OH stretching mode, which we assign to spectral diffusion. In the OH stretching region, dynamics calculations predict a more intense antisymmetric peak than experiments observe but do not predict the bimodal split symmetric peak that the experiments show. Our work suggests that electronic polarization, missing in the empirical force field, is responsible for the first discrepancy and that quantum nuclear effects, captured neither in density functional theory nor in classical dynamics, may be responsible for the second.

Enhanced Hydrogen Bonds of the (H2O)n (n = 4-8) Clusters Confined in Uranyl Peroxide Cluster Na20(UO2)20(O2)30

Water is a basic resource and an essential component of living organisms. It often exhibits some novel properties under confinement. The water clusters (H₂O)_n (n = 4-8) confined in the cavity of uranyl peroxide cluster Na₂₀(UO₂)₂₀(O₂)₃₀ (U₂₀) have been computationally investigated by using ab initio molecular dynamics (AIMD) simulations and density functional theory (DFT) calculations in this study. The results show that the confined water clusters can form hydrogen bonds with the internal oxygen atoms (O_uranyl) of U₂₀, and their conformations changed significantly. The average lengths (2.553-2.645 Å) of hydrogen bonds in confined (H₂O)_n are shorter than those (2.731-2.841 Å) in the corresponding free water clusters. Moreover, these confined hydrogen bonds show better hydrogen bond patterns according to the quantified indices. The natural bond orbital (NBO) calculations determine that there is electron transferring from the U₂₀ to its interior (H₂O)_n. It is the main reason for enhancing hydrogen bond interactions among the confined water molecules because their oxygen atoms are more negatively charged and their hydrogen atoms are more positively charged. The quantum theory of atoms in molecules (QTAIM) and interacting quantum atoms (IQA) analyses indicate that the confined hydrogen bonds are more covalent, based on the significant electron density ρ(r) and local energy density H(r) at the bond critical points (BCPs), and the stronger energies of interatomic exchange interactions (V_xc). These findings may help to promote the communication of confined water clusters and enrich the understating of confined hydrogen bonds.

Implementation of the self-consistent phonons method with ab initio potentials (AI-SCP)

The self-consistent phonon (SCP) method allows one to include anharmonic effects when treating a many-body quantum system at thermal equilibrium. The system is then described by an effective temperature-dependent harmonic Hamiltonian, which can be used to estimate its various dynamic and static properties. In this paper, we combine SCP with ab initio (AI) potential energy evaluation in which case the numerical bottleneck of AI-SCP is the evaluation of Gaussian averages of the AI potential energy and its derivatives. These averages are computed efficiently by the quasi-Monte Carlo method utilizing low-discrepancy sequences leading to a fast convergence with respect to the number, S, of the AI energy evaluations. Moreover, a further substantial (an-order-of-magnitude) improvement in efficiency is achieved once a numerically cheap approximation of the AI potential is available. This is based on using a perturbation theory-like (the two-grid) approach in which it is the average of the difference between the AI and the approximate potential that is computed. The corresponding codes and scripts are provided.

Ultrafast Water H-Bond Rearrangement in a Metal-Organic Framework Probed by Femtosecond Time-Resolved Infrared Spectroscopy

We investigated the water H-bond network and its dynamics in Ni₂Cl₂BTDD, a prototypical MOF for atmospheric water harvesting, using linear and ultrafast IR spectroscopy. Utilizing isotopic labeling and infrared spectroscopy, we found that water forms an extensive H-bonding network in Ni₂Cl₂BTDD. Further investigation with ultrafast spectroscopy revealed that water can reorient in a confined cone up to ∼50° within 1.3 ps. This large angle reorientation indicates H-bond rearrangement, similar to bulk water. Thus, although the water H-bond network is confined in Ni₂Cl₂BTDD, different from other confined systems, H-bond rearrangement is not hindered. The picosecond H-bond rearrangement in Ni₂Cl₂BTDD corroborates its reversibility with minimal hysteresis in water sorption.

Combinational Vibration Modes in H2O/HDO/D2O Mixtures Detected Thanks to the Superior Sensitivity of Femtosecond Stimulated Raman Scattering

Overtones and combinational modes frequently play essential roles in ultrafast vibrational energy relaxation in liquid water. However, these modes are very weak and often overlap with fundamental modes, particularly in isotopologues mixtures. We measured VV and HV Raman spectra of H₂O and D₂O mixtures with femtosecond stimulated Raman scattering (FSRS) and compared the results with calculated spectra. Specifically, we observed the mode at around 1850 cm^-1 and assigned it to H-O-D bend + rocking libration. Second, we found that the H-O-D bend overtone band and the OD stretch + rocking libration combination band contribute to the band located between 2850 and 3050 cm^-1. Furthermore, we assigned the broad band located between 4000 and 4200 cm^-1 to be composed of combinational modes of high-frequency OH stretching modes with predominantly twisting and rocking librations. These results should help in a proper interpretation of Raman spectra of aqueous systems as well as in the identification of vibrational relaxation pathways in isotopically diluted water.

Many-Body Contributions in Water Nanoclusters

Many-body interactions in water are known to be important but difficult to treat in atomistic models and often are included only as a correction. Polarizable models treat them explicitly, with long-range many-body potentials, within their classical approximation. However, their calculation is computationally expensive. Here, we evaluate how relevant the contributions to the many-body interaction associated with different coordination shells are. We calculate the global energy minimum, and the corresponding configuration, for nanoclusters of up to 20 water molecules. We find that including the first coordination shell, i.e., the five-body term of the central molecule, is enough to approximate within 5% the global energy minimum and its structure. We show that this result is valid for three different polarizable models, the Dang-Chang, the MB-pol, and the Kozack-Jordan potentials. This result suggests a strategy to develop many-body potentials for water that are reliable and, at the same time, computationally efficient.

Vibrational Signatures of HNO3 Acidity When Complexed with Microhydrated Alkali Metal Ions, M+·(HNO3)(H2O)n=5 (M = Li, K, Na, Rb, Cs), at 20 K

The speciation of strong acids like HNO₃ under conditions of restricted hydration is an important factor in the rates of chemical reactions at the air-water interface. Here, we explore the trade-offs at play when HNO₃ is attached to alkali ions (Li⁺-Cs⁺) with four water molecules in their primary hydration shells. This is achieved by analyzing the vibrational spectra of the M⁺·(HNO₃)(H₂O)₅ clusters cooled to about 20 K in a cryogenic photofragmentation mass spectrometer. The local acidity of the acidic OH group is estimated by the extent of the red shift in its stretching frequency when attached to a single water molecule. The persistence of this local structural motif (HNO₃-H₂O) in all of these alkali metal clusters enables us to determine the competition between the effect of the direct complexation of the acid with the cation, which acts to enhance acidity, and the role of the water network in the first hydration shell around the ions, which acts to counter (screen) the intrinsic effect of the ion. Analysis of the vibrational features associated with the acid molecule, as well as those of the water network, reveals how cooperative interactions in the microhydration regime conspire to effectively offset the intrinsic enhancement of HNO₃ acidity afforded by attachment to the smaller cations.

Energies, structures, and harmonic frequencies of small water clusters from the direct random phase approximation

The binding energies, structures, and vibrational frequencies of water clusters up to 20 molecules are computed at the direct random phase approximation (RPA) level of theory and compared to theoretical benchmarks. Binding energies of the WATER27 set, which includes neutral and positively and negatively charged clusters, are predicted to be too low in the complete basis set limit by an average of 7 kcal/mol (9%) and are worse than the results from the best density functional theory methods or from the Møller-Plesset theory. The RPA shows significant basis set size dependence for binding energies. The order of the relative energies of the water hexamer and dodecamer isomers is predicted correctly by the RPA. The mean absolute deviation for angles and distances for neutral clusters up to the water hexamer are 0.2° and 0.6 pm, respectively, using quintuple-ζ basis sets. The relative energetic order of the hexamer isomers is preserved upon optimization. Vibrational frequencies for these systems are underestimated by several tens of wavenumbers for large basis sets, and deviations increase with the basis set size. Overall, the direct RPA method yields accurate structural parameters but systematically underestimates binding energies and shows strong basis set size dependence.

Photoelectron Spectroscopy of the Water Dimer Reveals Unpredicted Vibrational Structure

Hydrogen bonds and proton transfer reactions can be considered as being at the very heart of aqueous chemistry and of utmost importance for many processes of biological relevance. Nevertheless, these processes are not yet well understood, even in seemingly simple model systems like small water clusters. We present a study of the photoelectron spectrum of the water dimer, revealing previously unresolved vibrational structure with 10-30 meV (80-242 cm^-1) typical splitting, in disagreement with a previous theoretical photoionization study predicting an apparent main vibrational progression with an ∼130 meV spacing [Kamarchik et al.; J. Chem. Phys. 2010, 132, 194311]. The observed vibrational structure and its deviation from the theoretical prediction is discussed in terms of known difficulties with calculations of strongly coupled anharmonic systems involving large amplitude motions. Potential contributions of the nonzero vibrational energy of the neutral water dimer at a finite experimental internal temperature are addressed. The internal temperature is estimated from the breakdown diagram associated with the dissociative ionization of the water dimer to be around to 130 K. This analysis also provides two additional, independently measured values for the 0 K appearance energy of the hydronium ion (H₃O⁺) from dissociative ionization of the water dimer.

Demystifying the Diffuse Vibrational Spectrum of Aqueous Protons Through Cold Cluster Spectroscopy

The ease with which the pH is routinely determined for aqueous solutions masks the fact that the cationic product of Arrhenius acid dissolution, the hydrated proton, or H⁺(aq), is a remarkably complex species. Here, we review how results obtained over the past 30 years in the study of H⁺⋅(H₂O)_n cluster ions isolated in the gas phase shed light on the chemical nature of H⁺(aq). This effort has also revealed molecular-level aspects of the Grotthuss relay mechanism for positive-charge translocation in water. Recently developed methods involving cryogenic cooling in radiofrequency ion traps and the application of two-color, infrared-infrared (IR-IR) double-resonance spectroscopy have established a clear picture of how local hydrogen-bond topology drives the diverse spectral signatures of the excess proton. This information now enables a new generation of cluster studies designed to unravel the microscopic mechanics underlying the ultrafast relaxation dynamics displayed by H⁺(aq).

Temperature evolution in IR action spectroscopy experiments with sodium doped water clusters

The combination of supersonic expansions with IR action spectroscopy techniques is the basis of many successful approaches to study cluster structure and dynamics. The effects of temperature and temperature evolution are important with regard to both the cluster synthesis and the cluster dynamics upon IR excitation. In the past the combination of the sodium doping technique with IR excitation enhanced near threshold photoionization has been successfully applied to study neutral, especially water clusters. In this work we follow an overall examination approach for inspecting the interplay of cluster temperature and cluster structure in the initial cooling process and in the IR excitation induced heating of the clusters. In molecular simulations, we study the temperature dependent photoionization spectra of the sodium doped clusters and the evaporative cooling process by water molecule ejection at the cluster surface. We present a comprehensive analysis that provides constraints for the temperature evolution from the nozzle to cluster detection in the mass spectrometer. We attribute the IR action effect to the strong temperature dependence of sodium solvation in the IR excited clusters and we discuss the effects of geometry changes during the IR multi-photon absorption process with regard to application prospects of the method.

Determining Partial Atomic Charges for Liquid Water: Assessing Electronic Structure and Charge Models

Partial atomic charges provide an intuitive and efficient way to describe the charge distribution and the resulting intermolecular electrostatic interactions in liquid water. Many charge models exist and it is unclear which model provides the best assignment of partial atomic charges in response to the local molecular environment. In this work, we systematically scrutinize various electronic structure methods and charge models (Mulliken, natural population analysis, CHelpG, RESP, Hirshfeld, Iterative Hirshfeld, and Bader) by evaluating their performance in predicting the dipole moments of isolated water, water clusters, and liquid water as well as charge transfer in the water dimer and liquid water. Although none of the seven charge models is capable of fully capturing the dipole moment increase from isolated water (1.85 D) to liquid water (about 2.9 D), the Iterative Hirshfeld method performs best for liquid water, reproducing its experimental average molecular dipole moment, yielding a reasonable amount of intermolecular charge transfer, and showing modest sensitivity to the local water environment. The performance of the charge model is dependent on the choice of the density functional and the quantum treatment of the environment. The computed molecular dipole moment of water generally increases with the percentage of the exact Hartree-Fock exchange in the functional, whereas the amount of charge transfer between molecules decreases. For liquid water, including two full solvation shells of surrounding water molecules (within about 5.5 Å of the central water) in the quantum chemical calculation converges the charges of the central water molecule. Our final pragmatic quantum chemical charge-assigning protocol for liquid water is the Iterative Hirshfeld method with M06-HF/aug-cc-pVDZ and a quantum region cutoff radius of 5.5 Å.

Isolating the Contributions of Specific Network Sites to the Diffuse Vibrational Spectrum of Interfacial Water with Isotopomer-Selective Spectroscopy of Cold Clusters

Decoding the structural information contained in the interfacial vibrational spectrum of water requires understanding how the spectral signatures of individual water molecules respond to their local hydrogen bonding environments. In this study, we isolated the contributions for the five classes of sites that differ according to the number of donor (D) and acceptor (A) hydrogen bonds that characterize each site. These patterns were measured by exploiting the unique properties of the water cluster cage structures formed in the gas phase upon hydration of a series of cations M⁺·(H₂O)_n (M = Li, Na, Cs, NH₄, CH₃NH₃, H₃O, and n = 5, 20-22). This selection of ions was chosen to systematically express the A, AD, AAD, ADD, and AADD hydrogen bonding motifs. The spectral signatures of each site were measured using two-color, IR-IR isotopomer-selective photofragmentation vibrational spectroscopy of the cryogenically cooled, mass selected cluster ions in which a single intact H₂O is introduced without isotopic scrambling, an important advantage afforded by the cluster regime. The resulting patterns provide an unprecedented picture of the intrinsic line shapes and spectral complexities associated with excitation of the individual OH groups, as well as the correlation between the frequencies of the two OH groups on the same water molecule, as a function of network site. The properties of the surrounding water network that govern this frequency map are evaluated by dissecting electronic structure calculations that explore how changes in the nearby network structures, both within and beyond the first hydration shell, affect the local frequency of an OH oscillator. The qualitative trends are recovered with a simple model that correlates the OH frequency with the network-modulated local electron density in the center of the OH bond.

Temperature and Nuclear Quantum Effects on the Stretching Modes of the Water Hexamer

The water hexamer has many low-lying isomers, e.g., ring, book, cage, and prism, shifting from two- to three-dimensional structures. We show that this dimensionality change is accompanied by a drop in the quantum nature of the cluster, as manifested in the red shift of the quantal OH stretching modes as compared with their classical counterparts. We obtain this "nuclear quantum effect" (NQE) as the mean deviation between the OH stretch frequencies from velocity autocorrelation Fourier transforms from classical trajectories on a high-level water potential (MB-pol) as compared with scaled harmonic frequencies from high-level quantum chemistry calculations. With a universal scaling factor, the predicted OH frequencies agree with experiment to a mean absolute deviation ≤10 cm-1, which allows unequivocal isomer assignments. By assuming temperature-independent NQEs, we produce the temperature dependence of the cage isomer OH stretch spectrum below 70 K, where it is the dominant structure. All bands widen and blue-shift with increasing temperature, most conspicuously the reddest mode, which thus constitutes a "vibrational thermometer".

Unusually strong hydrogen bond cooperativity in particular (H2O)20 clusters

Drawing upon an intuitive charge-transfer-based picture of hydrogen bonding, we demonstrate that cooperativity effects acting in concert can lead to unusually strong hydrogen bonds in neutral water clusters. The structure, vibrational, and NMR properties of a (H2O)20 pentagonal dodecahedron cluster containing such a strong hydrogen bond were studied using second-order perturbation theory and density functional theory. The hydrogen bond length was found to be shorter than 2.50 Å. A large redshift of over 2000 cm-1 with respect to the isolated water molecule was predicted for the OH stretching frequency of the donor water molecule. A large downfield shift to 13.5 ppm of the isotropic part of the 1H magnetic shielding tensor together with an unusually large shielding anisotropy of 49.9 ppm was obtained. The hydrogen bond energy was calculated using symmetry-adapted perturbation theory and was found to be more than three times stronger than a typical hydrogen bond in liquid water.

Vibrational Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction

Vibrational spectroscopy is an essential tool in chemical analyses, biological assays, and studies of functional materials. Over the past decade, various coherent nonlinear vibrational spectroscopic techniques have been developed and enabled researchers to study time-correlations of the fluctuating frequencies that are directly related to solute-solvent dynamics, dynamical changes in molecular conformations and local electrostatic environments, chemical and biochemical reactions, protein structural dynamics and functions, characteristic processes of functional materials, and so on. In order to gain incisive and quantitative information on the local electrostatic environment, molecular conformation, protein structure and interprotein contacts, ligand binding kinetics, and electric and optical properties of functional materials, a variety of vibrational probes have been developed and site-specifically incorporated into molecular, biological, and material systems for time-resolved vibrational spectroscopic investigation. However, still, an all-encompassing theory that describes the vibrational solvatochromism, electrochromism, and dynamic fluctuation of vibrational frequencies has not been completely established mainly due to the intrinsic complexity of intermolecular interactions in condensed phases. In particular, the amount of data obtained from the linear and nonlinear vibrational spectroscopic experiments has been rapidly increasing, but the lack of a quantitative method to interpret these measurements has been one major obstacle in broadening the applications of these methods. Among various theoretical models, one of the most successful approaches is a semiempirical model generally referred to as the vibrational spectroscopic map that is based on a rigorous theory of intermolecular interactions. Recently, genetic algorithm, neural network, and machine learning approaches have been applied to the development of vibrational solvatochromism theory. In this review, we provide comprehensive descriptions of the theoretical foundation and various examples showing its extraordinary successes in the interpretations of experimental observations. In addition, a brief introduction to a newly created repository Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art multidimensional vibrational spectroscopy will be one of the most fruitful ways to study the structure and dynamics of chemical, biological, and functional molecular systems in the future.

Infrared spectroscopy of neutral water clusters at finite temperature: Evidence for a noncyclic pentamer

Infrared spectroscopic study of neutral water clusters is crucial to understanding of the hydrogen-bonding networks in liquid water and ice. Here we report infrared spectra of size-selected neutral water clusters, (H₂O) _n (n = 3-6), in the OH stretching vibration region, based on threshold photoionization using a tunable vacuum ultraviolet free-electron laser. Distinct OH stretch vibrational fundamentals observed in the 3,500-3,600-cm^-1 region of (H₂O)₅ provide unique spectral signatures for the formation of a noncyclic pentamer, which coexists with the global-minimum cyclic structure previously identified in the gas phase. The main features of infrared spectra of the pentamer and hexamer, (H₂O) _n (n = 5 and 6), span the entire OH stretching band of liquid water, suggesting that they start to exhibit the richness and diversity of hydrogen-bonding networks in bulk water.

Terahertz and mid-infrared spectroscopy of matrix-isolated clusters and matrix-sublimation ice of D2O

We have established an apparatus for terahertz and mid-infrared spectroscopy in an ultrahigh vacuum and have measured absorption spectra of D₂O clusters trapped in solid Ar. To assign terahertz absorption peaks due to the D₂O dimer, trimer, and tetramer, the dependence of the spectrum on the annealing temperature and D₂O dilution was analyzed. The assignment was also examined by ab initio calculations with the use of the "our own N-layered integrated molecular orbital and molecular mechanics" method, where the flexibility of surrounding Ar atoms was systematically incorporated. We identified all the intermolecular fundamentals of the dimer and those with significant intensities of the trimer and tetramer, whose structural symmetries were revealed to be broken down. After isolating the D₂O clusters in solid Ar, we sublimated only Ar atoms to leave behind matrix-sublimation ice, which was found to be amorphous- or crystal-like depending on the formation conditions: the dilution and sublimation temperature. The crystallinity of matrix-sublimation ice was determined by decomposing its terahertz spectrum into the spectra of amorphous and crystalline ices. Since the crystallinity got higher by raising the dilution and sublimation temperature, the diffusion of the D₂O monomer on the surface of sublimating solid Ar was found to be crucial to the crystallization of the sublimation ice.

Terahertz VRT Spectroscopy of the Water Hexamer-h12 Cage: Dramatic Libration-Induced Enhancement of Hydrogen Bond Tunneling Dynamics

We report the assignment and analysis of 176 transitions belonging to a librational band of the (H₂O)₆ cage isomer near 525 cm^-1(15 THz). From a fit of the transitions to an asymmetric top model, we observe both dramatic changes in the rotational constants relative to the ground state, indicating significant nonrigidity, and striking enhancement in the tunneling motions that break and reform the hydrogen bonds in the cluster. This is the fifth water cluster system to display such an enhancement in the 15 THz librational region, the details of which may help to elucidate the hydrogen bond dynamics occurring in bulk liquid water.

Hydrogen bond network structures of protonated short-chain alcohol clusters

Protonated alcohol clusters enable extraction of the physical essence of the nature of hydrogen bond networks.

Capturing the Elusive Water Trimer from the Stepwise Growth of Water on the Surface of the Polycyclic Aromatic Hydrocarbon Acenaphthene

Polycyclic aromatic hydrocarbons (PAHs) are key players in reaction chemistry. While it is postulated that they serve as a basis for ice grains, there has been no direct detection of PAHs in astronomical environments. We aim to investigate the hydration of PAHs to set a foundation for the future exploration of potential ice formation pathways. We report results from chirped pulse Fourier transform microwave spectroscopy and quantum-chemical calculations for the PAH acenaphthene and acenaphthene complexed with up to four water molecules. The acenaphthene-(H₂O)₃ complex is of particular interest as the elusive cyclic water trimer was observed. It appears in a slightly distorted configuration when compared with the pure water trimer. This is explained by hydrogen-bond net cooperativity effects. Binding energies for the complexes are presented. Our results provide insight into the onset of complex aggregation that could be occurring in extraterrestrial environments as part of ice grain formation.

Exploring water radiolysis in proton cancer therapy: Time-dependent, non-adiabatic simulations of H+ + (H2O)1-6

To elucidate microscopic details of proton cancer therapy (PCT), we apply the simplest-level electron nuclear dynamics (SLEND) method to H⁺ + (H₂O)_1-6 at E_Lab = 100 keV. These systems are computationally tractable prototypes to simulate water radiolysis reactions—i.e. the PCT processes that generate the DNA-damaging species against cancerous cells. To capture incipient bulk-water effects, ten (H₂O)_1-6 isomers are considered, ranging from quasi-planar/multiplanar (H₂O)_1-6 to “smallest-drop” prism and cage (H₂O)₆ structures. SLEND is a time-dependent, variational, non-adiabatic and direct method that adopts a nuclear classical-mechanics description and an electronic single-determinantal wavefunction in the Thouless representation. Short-time SLEND/6-31G* (n = 1–6) and /6-31G** (n = 1–5) simulations render cluster-to-projectile 1-electron-transfer (1-ET) total integral cross sections (ICSs) and 1-ET probabilities. In absolute quantitative terms, SLEND/6-31G* 1-ET ICS compares satisfactorily with alternative experimental and theoretical results only available for n = 1 and exhibits almost the same accuracy of the best alternative theoretical result. SLEND/6-31G** overestimates 1-ET ICS for n = 1, but a comparable overestimation is also observed with another theoretical method. An investigation on H⁺ + H indicates that electron direct ionization (DI) becomes significant with the large virtual-space quasi-continuum in large basis sets; thus, SLEND/6-31G** 1-ET ICS is overestimated by DI contributions. The solution to this problem is discussed. In relative quantitative terms, both SLEND/6-31* and /6-31G** 1-ET ICSs precisely fit into physically justified scaling formulae as a function of the cluster size; this indicates SLEND’s suitability for predicting properties of water clusters with varying size. Long-time SLEND/6-31G* (n = 1–4) simulations predict the formation of the DNA-damaging radicals H, OH, O and H₃O. While “smallest-drop” isomers are included, no early manifestations of bulk water PCT properties are observed and simulations with larger water clusters will be needed to capture those effects. This study is the largest SLEND investigation on water radiolysis to date.

Size-Resolved Photoelectron Anisotropy of Gas Phase Water Clusters and Predictions for Liquid Water

We report the first measurements of size-resolved photoelectron angular distributions for the valence orbitals of neutral water clusters with up to 20 molecules. A systematic decrease of the photoelectron anisotropy is found for clusters with up to 5-6 molecules, and most remarkably, convergence of the anisotropy for larger clusters. We suggest the latter to be the result of a local short-range scattering potential that is fully described by a unit of 5-6 molecules. The cluster data and a detailed electron scattering model are used to predict the anisotropy of slow photoelectrons in liquid water. Reasonable agreement with experimental liquid jet data is found.

Modeling Molecular Interactions in Water: From Pairwise to Many-Body Potential Energy Functions

Almost 50 years have passed from the first computer simulations of water, and a large number of molecular models have been proposed since then to elucidate the unique behavior of water across different phases. In this article, we review the recent progress in the development of analytical potential energy functions that aim at correctly representing many-body effects. Starting from the many-body expansion of the interaction energy, specific focus is on different classes of potential energy functions built upon a hierarchy of approximations and on their ability to accurately reproduce reference data obtained from state-of-the-art electronic structure calculations and experimental measurements. We show that most recent potential energy functions, which include explicit short-range representations of two-body and three-body effects along with a physically correct description of many-body effects at all distances, predict the properties of water from the gas to the condensed phase with unprecedented accuracy, thus opening the door to the long-sought "universal model" capable of describing the behavior of water under different conditions and in different environments.

Infrared detection of (H2O)20 isomers of exceptional stability: a drop-like and a face-sharing pentagonal prism cluster

Water clusters with internally solvated water molecules are widespread models that mimic the local environment of the condensed phase. The appearance of stable (H2O)n cluster isomers having a fully coordinated interior molecule has been theoretically predicted to occur around the n = 20 size range. However, our current knowledge about the size regime in which those structures become energetically more stable has remained hypothetical from simulations in lieu of the absence of precisely size-resolved experimental measurements. Here we report size and isomer selective infrared (IR) spectra of (H2O)20 clusters tagged with a sodium atom by employing IR excitation modulated photoionization spectroscopy. The observed absorption patterns in the OH stretching region are consistent with the theoretically predicted spectra of two structurally distinct isomers of exceptional stability: a drop-like cluster with a fully coordinated (interior) water molecule and an edge-sharing pentagonal prism cluster in which all atoms are on the surface. The drop-like structure is the first experimentally detected water cluster exhibiting the local connectivity found in liquid water.

Dielectric constant and low-frequency infrared spectra for liquid water and ice Ih within the E3B model

Two intrinsic difficulties in modeling condensed-phase water with conventional rigid non-polarizable water models are: reproducing the static dielectric constants for liquid water and ice Ih, and generating the peak at about 200 cm(-1) in the low-frequency infrared spectrum for liquid water. The primary physical reason for these failures is believed to be the missing polarization effect in these models, and consequently various sophisticated polarizable water models have been developed. However, in this work we pursue a different strategy and propose a simple empirical scheme to include the polarization effect only on the dipole surface (without modifying a model's intermolecular interaction potential). We implement this strategy for our explicit three-body (E3B) model. Our calculated static dielectric constants and low-frequency infrared spectra are in good agreement with experiment for both liquid water and ice Ih over wide temperature ranges, albeit with one fitting parameter for each phase. The success of our modeling also suggests that thermal fluctuations about local minima and the energy differences between different proton-disordered configurations play minor roles in the static dielectric constant of ice Ih. Our analysis shows that the polarization effect is important in resolving the two difficulties mentioned above and sheds some light on the origin of several features in the low-frequency infrared spectra for liquid water and ice Ih.

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.

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.