The Vibrational Spectrum of the Neutral (H2O)6 Precursor to the “Magic” (H2O)6- Cluster Anion by Argon-Mediated, Population-Modulated Electron Attachment Spectroscopy

Correlations between the Structures and Spectra of Protonated Water Clusters

Badger's rule-like correlations between OH stretching frequencies and intensities and the OH bond length are used to develop a spectral mapping procedure for studies of pure and protonated water clusters. This approach utilizes the vibrationally averaged OH bond lengths, which were obtained from diffusion Monte Carlo simulations that were performed using the general potential developed by Yu and Bowman. Good agreement is achieved between the spectra obtained using this approach and previously reported spectra for H+(H2O)n clusters, with n = 3, 4, and 5, as well as their perdeuterated analogues. The analysis of the spectra obtained by this spectral mapping approach supports previous work that assigned the spectrum of H+(H2O)6 to a mixture of Eigen and Zundel-like structures. Analysis of the calculated spectra also suggests a reassignment of the frequency of one of the transitions that involves the OH stretching vibration of the OH bonds in the hydronium core in the Eigen-like structure of H+(H2O)6 from 1917 cm-1 to roughly 2100 cm-1. For D+(D2O)6, comparison of the measured spectrum to those obtained by using the spectral mapping approach suggests that the carrier of the measured spectrum is one or more of the isomers of D+(D2O)6 that contain a four-membered ring and two flanking water molecules. While there are several candidate structures, the two flanking water molecules most likely form a chain that is bound to the hydronium core.

Microsolvation of Hot Ions: Spectroscopy and Statistical Mechanics of Phenide-Water Interactions

Thermal excitation alters the spectroscopic signatures of solvated ions and affects their interactions with neighboring molecules. By analyzing the photoelectron spectra of microhydrated phenide (Ph-), the temperatures of the Ph-·H2O and Ph-·(H2O)2 clusters from a hot ion source were determined to be 560 and 520 K, respectively, vs 700 K for unsolvated Ph-. Compared to theory predictions for cold clusters, the high temperature of the environment significantly reduces the average hydration stabilization of the ions and the corresponding band shifts in their spectra. The results are discussed in terms of a statistical model that describes the energy content of the intermolecular (IM) degrees of freedom of the cluster, ⟨EIM⟩. We show that over the entire solvation energy range, the density of states associated with the IM modes of Ph-·H2O, of which there are only 6, is more than an order of magnitude greater than that associated with the 27 internal vibrations of the core anion. The results suggest that the observed cluster temperatures are not determined by the ion source but represent the intrinsic properties of the clusters. The energetics and statistical mechanics of microsolvation limit the excitation that the IM degrees of freedom can sustain without significant solvent evaporation on the timescale of the experiment. The limit is expressed as a characteristic solvation temperature (CST), which is the maximum canonical temperature of a stable cluster ensemble. Driven by evaporative cooling, the terminal cluster temperature from a hot ion source will always be close to the cluster's CST. Only if the source temperature is lower than CST will the observed cluster temperature be determined by the source conditions. An approximate rule is proposed for estimating the characteristic temperature of any cluster using the inflection point on the ⟨EIM⟩ vs T curve.

Anharmonic Assignment of the Water Octamer Spectrum in the OH Stretch Region

We interface the quasi-classical trajectory approach with an ab initio potential energy surface for water to assign the vibrational spectroscopical features of the OH stretch region of the water octamer cluster, which is considered to be a precursor of ice. An attempt by Li et al. to assign their recent reference experiment involved lower-level calculations based on an ad hoc scaled harmonic approach. Differently from the conclusions of this previous assignment, which invoked the contribution of 5 conformers and a solvated form of the water heptamer in the spectrum, we find out that the spectroscopic features can be related to the 4 conformers of the octamer lying lower in energy.

Δ-Machine Learned Potential Energy Surfaces and Force Fields

There has been great progress in developing machine-learned potential energy surfaces (PESs) for molecules and clusters with more than 10 atoms. Unfortunately, this number of atoms generally limits the level of electronic structure theory to less than the "gold standard" CCSD(T) level. Indeed, for the well-known MD17 dataset for molecules with 9-20 atoms, all of the energies and forces were obtained with DFT calculations (PBE). This Perspective is focused on a Δ-machine learning method that we recently proposed and applied to bring DFT-based PESs to close to CCSD(T) accuracy. This is demonstrated for hydronium, N-methylacetamide, acetyl acetone, and ethanol. For 15-atom tropolone, it appears that special approaches (e.g., molecular tailoring, local CCSD(T)) are needed to obtain the CCSD(T) energies. A new aspect of this approach is the extension of Δ-machine learning to force fields. The approach is based on many-body corrections to polarizable force field potentials. This is examined in detail using the TTM2.1 water potential. The corrections make use of our recent CCSD(T) datasets for 2-b, 3-b, and 4-b interactions for water. These datasets were used to develop a new fully ab initio potential for water, termed q-AQUA.

Confocal Raman microscopic evidence for cyclic water Pentamer, at high temperatures in a supramolecular host of [Cu(cyclam)(N3)2]·4H2O

Here, we report the existence of a pentameric water cluster in the host framework of [Cu(cyclam)(N₃)2]·4H₂O, that is stable upto 167 °C, well above the boiling point of water. The pentameric cluster structure embedded in the host framework is evident from the single crystal studies. The high thermal stability is confirmed by TGA and temperature dependent confocal Raman microscopic studies, where loss of water bands is well captured between 167 and 170 °C, besides its existence through SCXRD studies. To the best of our knowledge, this is the first report where temperature dependent confocal Raman microscopic investigation is used to study the stability of water in crystal environment. The study promises that temperature dependent confocal Raman microscopy can be an efficient tool to investigate the existence and stability of small water clusters, precisely in restricted environments.

Infrared spectroscopic study of hydrogen bonding topologies in the smallest ice cube

The water octamer with its cubic structure consisting of six four-membered rings presents an excellent cluster system for unraveling the cooperative interactions driven by subtle changes in the hydrogen-bonding topology. Despite prediction of many distinct structures, it has not been possible to extract the structural information encoded in their vibrational spectra because this requires size-selectivity of the neutral clusters with sufficient resolution to identify the contributions of the different isomeric forms. Here we report the size-specific infrared spectra of the isolated cold, neutral water octamer using a scheme based on threshold photoionization using a tunable vacuum ultraviolet free electron laser. A plethora of sharp vibrational bands features are observed. Theoretical analysis of these patterns reveals the coexistence of five cubic isomers, including two with chirality. The relative energies of these structures are found to reflect topology-dependent, delocalized multi-center hydrogen-bonding interactions. These results demonstrate that even with a common structural motif, the degree of cooperativity among the hydrogen-bonding network creates a hierarchy of distinct species. The implications of these results on possible metastable forms of ice are speculated.

Spectroscopic studies of water clusters provide insight into the hydrogen bond structure of water and ice. The authors measure infrared spectra of neutral water octamers using a threshold photoionization technique based on a tunable vacuum-UV free electron laser, identifying two cubic isomers in addition to those previously observed.

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".

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.

Infrared Spectroscopy of Neutral Water Dimer Based on a Tunable Vacuum Ultraviolet Free Electron Laser

Infrared (IR) spectroscopy provides detailed structural and dynamical information on clusters at the fingerprint level. Herein, we demonstrate the capability of a tunable vacuum ultraviolet free electron laser (VUV-FEL) for selective detection of a wide variety of neutral water clusters and for recording the size-dependent IR spectra. The present technique does not require the presence of an ultraviolet chromophore or a dipole moment and is generally applicable for IR spectroscopy of neutral clusters free from confinement. To show the features of our technique, we report here the IR spectra of neutral water dimer in the OH stretch region, providing benchmarks for theoretical study of the accurate description of hydrogen bonding structures involved in liquid water and ice. Quantum mechanical calculations on a 12-dimensional ab initio potential energy surface are utilized to simulate the anharmonic vibrational spectra of water dimer. These results help to resolve the controversy of the exact vibrational assignment of each band feature of the water dimer.

Importance of self-interaction-error removal in density functional calculations on water cluster anions

Removing self-interaction errors in density functional approximations results in significantly improved vertical detachment energies of water anions and is essential for obtaining orbital energies consistent with electron binding energies.

Predicting stability limits for pure and doped dicationic noble gas clusters undergoing coulomb explosion: A parallel tempering based study

We have used a replica exchange Monte-Carlo procedure, popularly known as Parallel Tempering, to study the problem of Coulomb explosion in homogeneous Ar and Xe dicationic clusters as well as mixed Ar-Xe dicationic clusters of varying sizes with different degrees of relative composition. All the clusters studied have two units of positive charges. The simulations reveal that in all the cases there is a cutoff size below which the clusters fragment. It is seen that for the case of pure Ar, the value is around 95 while that for Xe it is 55. For the mixed clusters with increasing Xe content, the cutoff limit for suppression of Coulomb explosion gradually decreases from 95 for a pure Ar to 55 for a pure Xe cluster. The hallmark of this study is this smooth progression. All the clusters are simulated using the reliable potential energy surface developed by Gay and Berne (Gay and Berne, Phys. Rev. Lett. 1982, 49, 194). For the hetero clusters, we have also discussed two different ways of charge distribution, that is one in which both positive charges are on two Xe atoms and the other where the two charges are at a Xe atom and at an Ar atom. The fragmentation patterns observed by us are such that single ionic ejections are the favored dissociating pattern. © 2017 Wiley Periodicals, Inc.

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.

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.

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.

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.

A parallel tempering based study of Coulombic explosion and identification of dissociating fragments in charged noble gas clusters

In this communication, we would like to test the feasibility of a parallel tempering based study of dissociation in dicationic noble gas clusters, namely, Ar(n)(2+), Kr(n)(2+), and Xe(n)(2+), where "n" is the size of the cluster units. We would like to find out the correct limit for sizes of each of these systems, above which the clusters stay intact as a single unit and does not dissociate into fragments by the process of Coulomb explosion. Moreover, we would also like to, for a specific case, i.e., Ar(n)(2+), study in detail the fragmentation patterns and point out the switchover from the non-fission way to the fission mechanism of dissociation. In all these calculations, we would like to analyse, how close we are in our predictions with that of experimental results. As a further check on the dissociating patterns found out by parallel tempering, we also conduct basin hopping based study on representative sizes of the clusters and find that parallel tempering, as used for this present work as an optimizer, is able to predict correct features when compared with other celebrated methods like the basin hopping algorithm.

IR Spectra of the Water Hexamer: Theory, with Inclusion of the Monomer Bend Overtone, and Experiment Are in Agreement

Signature IR spectra of isomers of the water hexamer in the spectral range 3000-3800 cm(-1) have been reported by experimentalists, but crucial theoretical interpretation has still not been definitive. Using ab initio potential and dipole moment surfaces and a fully coupled quantum treatment of the intramolecular modes, the ring and book are assigned to spectra obtained in the He nanodroplet and Ar tagging experiments, respectively. The overtone of the intramolecular bend at ca. 3200 cm(-1) is a new calculated feature that completes an important missing piece in previous experimental and theoretical comparisons and leads to a consistent assignment of these two experimental spectra. Calculated IR spectra for the lowest energy forms of the water heptamer and octomer are also presented and compared to experiment. In all the calculated spectra, the bend overtone is demonstrated to be a noticeable feature, and this is one important conclusion from the work. Also, the danger in using scaled double-harmonic spectra to assign spectra is demonstrated.

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

Using a newly developed and recently parameterized classical empirical simulation model for water that involves explicit three-body interactions, we determine the eleven most stable isomers of the water hexamer. We find that the lowest energy isomer is one of the cage structures, in agreement with far-IR and microwave experiments. The energy ordering for the binding energies is cage > glove > book > bag > chair > boat > chaise, and energies relative to the cage are in good agreement with CCSD(T) calculations. The three-body contributions to the cage, book, and chair are also in reasonable agreement with CCSD(T) results. The energy of each isomer results from a delicate balance involving the number of hydrogen bonds, the strain of these hydrogen bonds, and cooperative and anti-cooperative three-body interactions, whose contribution we can understand simply from the form of the three-body interactions in the simulation model. Oxygen-oxygen distances in the cage and book isomers are in good agreement with microwave experiments. Hydrogen-bond distances depend on both donor and acceptor, which can again be understood from the three-body model. Fully anharmonic OH-stretch spectra are calculated for these low-energy structures, and compared with shifted harmonic results from ab initio and density functional theory calculations. Replica-exchange molecular dynamics simulations were performed from 40 to 194 K, which show that the cage isomer has the lowest free energy from 0 to 70 K, and the book isomer has the lowest free energy from 70 to 194 K. OH-stretch spectra were calculated between 40 and 194 K, and results at 40, 63, and 79 K were compared to recent experiments, leading to re-assignment of the peaks in the experimental spectra. We calculate local OH-stretch cumulative spectral densities for different donor-acceptor types and compare to analogous results for liquid water.

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

We present a unified picture of how OH-stretch spectroscopy in water can be understood in terms of hydrogen bonding for the four systems listed in the title. To understand the strength, and hence OH-stretch frequency, of a hydrogen bond, it is crucial to consider the number of additional acceptor hydrogen bonds made by both the donor and acceptor molecules. This necessity for focusing on the hydrogen-bond environment of both donor and acceptor molecules follows from quantum chemical considerations and is related to the three-body interactions in water. Armed with this understanding we can make a detailed interpretation of the OH-stretch IR absorption spectrum of the cage conformer for HOD(D2O)5 and the imaginary part of the ssp OH-stretch sum-frequency spectrum of the surface of liquid D2O with dilute HOD.

Low-lying isomers and finite temperature behavior of (H2O)6−

(H2O)(6) (-) appears as a "magic" number water cluster in (H2O)(n) (-) mass spectra. The structure of the (H2O)(6) (-) isomer dominating the experimental population has been established only recently [N. I. Hammer et al., J. Phys. Chem. A 109, 7896 (2005)], and the most noteworthy characteristic of this isomer is the localization of the excess electron in the vicinity of a double-acceptor monomer. In the present work, we use a quantum Drude model to characterize the low-energy isomers and the finite temperature properties of (H2O)(6) (-). Comparison with ab initio calculations shows that the use of a water model employing distributed polarizabilities and distributed repulsive sites is necessary to correctly reproduce the energy ordering of the low-lying isomers. Both the simulations and the ab initio calculations predict that there are several isomers of (H2O)(6) (-) significantly lower in energy than the experimentally observed species, suggesting that the experimental distribution is far from equilibrium.