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

Water clusters with internally solvated water molecules are wide-spread models that mimic the local environment of the condensed phase. The appearance of stable (H 2 O) 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 (H 2 O) 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.

Water clusters with internally solvated water molecules are widespread models that mimic the local environment of the condensed phase.The appearance of stable (H 2 O) 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 (H 2 O) 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.
Water clusters are of relevance to many fields of science like atmospheric and astrophysical chemistry 1,2 and have been reported to be trapped inside inorganic crystals, 3 polymers 4 and in crystal structures of hydrophobic proteins. 5In an aqueous solution the chemical reactivity is coupled to the subtle balance between the solvent-solute and solvent-solvent interactions that ultimately control the fleeting hydrogen bonding network.Water clusters provide a medium for probing the underlying solvent-solvent interactions in those aqueous environments and quantifying the properties of the hydrogen bonding interactions at the molecular level.Contrary to the fully connected network in liquid water, the structures of the first few water clusters exhibit arrangements in which all atoms are on the surface of the cluster, [6][7][8][9][10][11][12] a consequence of the need to maximize the number of hydrogen bonds in combination with the large surface-to-volume ratio that is characteristic of a small water cluster.4][15][16] However, all water molecules in these clusters are under-coordinated, having less than four neighbors, the number found on the average in bulk water.Previous theoretical studies 17,18 have identified a cluster size regime starting at around 17 molecules for which the transition between the ''all-surface'' to ''interior'' structures, viz. the appearance of stable cluster isomers with an interior, fully coordinated water molecule, may occur.The fact that experimental, size-resolved probes of either the infrared (IR) ''fingerprint'' region (3000 to 4000 cm À1 ) 14,15,[19][20][21] or the broadband rotational spectrum 11,22 have so far been limited to smaller sized clusters, has prevented the direct experimental probing of this important size regime.Therefore, these previous theoretical studies have been unchallenged by experiments, in contrast to the case of aqueous ionic clusters. 23,24Here we report size and isomer selective infrared (IR) spectra of (H 2 O) 20 clusters tagged with a sodium atom in the OH stretching region.The ensuing analysis will demonstrate that the isomer composition of this water cluster size is consistent with predictions of high level electronic structure calculations.
Hydrogen bonding interactions of fully coordinated water molecules that are characteristic of condensed environments consist of arrangements in which each molecule, on the average, both donates and accepts two hydrogen bonds (double donordouble acceptor configuration: DDAA). 25This connectivity leads to a complex, three-dimensional network of four-coordinated water molecules that is transient in the liquid 26 and static in solid water. 27Therefore, a special interest arises in identifying the emergence of water clusters which posses fully coordinated water molecules.In this context, the neutral (H 2 O) 20 clusters have previously served as a reference system for benchmarking theoretical methods, see ref. 28-31 and references cited therein.Additional reasons for studying clusters of this size is their potential connection to the ''magic number'' protonated water clusters, 23,29,32 the dodecahedral cage network of one of its families of isomers that is the building block of the structure I (sI) clathrate lattices 33,34 as well as their use as a model for cloud droplets in atmospheric reactions. 2he frequencies of the OH stretch oscillators can be used as a sensitive probe of the strength of the underlying hydrogen bonds: the stronger the intermolecular bond, the weaker the corresponding O-H covalent bond is, which leads to sizeable red shifts of such oscillators in the IR spectrum. 35To this end, vibrational spectroscopy of size-selected, neutral clusters is a powerful tool that probes the characteristic absorption patterns of the hydrogen bonded OH stretching vibrations 16 which, in conjunction with theoretical calculations, can be used to assign the underlying hydrogen bonded network via the ''structural-spectral correspondence'' principle. 36y applying such an experimental approach, the emergence of four-coordinated water molecules 37 for n 4 18 and the onset of crystallinity in larger clusters 38 have been characterized in previous experimental studies.We have recently reported an experimental approach, which allows to record precisely size-resolved infrared spectra in the OH stretching region of sodium tagged, neutral water clusters. 38,39The viability of the new method for arriving at structure assignments was indicated for sizes above n = 20, 39 but not rigorously proven.In the present work, the approach is further refined and coupled with first principles electronic structure calculations to assign the measured spectra for n = 20.
In the experiment, a molecular beam containing the clusters is formed in a skimmed supersonic expansion.The clusters are then doped with a single sodium atom in a pick-up cell and subsequently irradiated with IR photons inducing a signal enhancement to the UV photoionization, which is delayed by 80 ns.A single scan of the OH stretch region (2850 to 3800 cm À1 ) with a tuneable IR laser system provides simultaneously measured, size-resolved IR spectra for the complete cluster size distribution.The problem of UV photon induced fragmentation can be suppressed already in small clusters 40,41 due to the special ionization mechanism of sodium doped, hydrogen bonded clusters. 42For larger clusters (n 4 18) a spectral contamination by the sodium atom was not observed. 38In this study we used neon (Ne) and argon (Ar) instead of helium as seeding gas.The reason is the better suitability of Ne and Ar expansions to cool water clusters 41 and to relax them in their thermodynamical minima. 43Further experimental details are found in 38,40,42 and in the ESI, † in which also the reproducibility of the main features in the presented spectra is shown.
The electronic structure calculations were performed at the second-order Møller-Plesset perturbation theory (MP2) level of theory with Dunning's augmented, correlation-consistent basis sets (aug-cc-pVnZ) 44 of double through quadruple (n = D, T, Q) zeta quality allowing for the extrapolation of the energetics at the Complete Basis Set (CBS) limit (see ref. 28 for details).Geometries were optimized with the double and triple zeta basis sets whereas single point energies were obtained with the larger quadruple zeta set at the triple zeta geometries.Harmonic vibrational spectra and subsequent zero-point energy (ZPE) corrections are obtained with the aug-cc-pVDZ basis set.These calculation have been conducted on isomer B in the present study, the results for isomers A and C were taken from previous work. 28All electronic structure calculations were carried out with the NWChem suite of codes. 459][30][31] A ''family'' of isomers denotes clusters with the same oxygen network differing only in the position of the hydrogen atoms according to the Bernal-Fowler ice rules. 46he stability of the edge-sharing pentagonal prisms, first suggested using classical force fields, 47 was confirmed as the most stable family of minima from high level electronic structure calculations.The lowest isomer of this family 28 is shown in Fig. 1 (isomer A) together with the predicted infrared spectrum in the OH stretching region.It is a highly symmetric structure (as regards the oxygen atom network) and features in total eight water molecules in DDAA configurations, which are all placed at the surface of the cluster.Additionally, there are six water molecules having a single donor-double acceptor motif (DAA), which show a characteristic single peak in the most red-shifted part of the OH stretching spectrum.
An energetically competitive, (cf.Table 1) drop-like 31 isomer of (H 2 O) 20 is also considered (isomer B in Fig. 1).This isomer features a total of six fully coordinated and interconnected water molecules, five of which are on the surface and one DDAA-like fragment at the center of the cluster.Four out of the former five surface fragments are hydrogen bonded to the latter interior one.The surface of isomer B consists of four 4-, six 5-and two 6-member rings.Interconnected 4-and 6-member rings cover roughly one half of the surface area, whereas interconnected 5-member rings cover the other half.The most red-shifted part in the ''fingerprint'' OH stretching region of the IR spectra provides a definitive feature used for the structural identification of isomer B. There exists a characteristic pattern of seven DAA oscillators which arrange in three peaks of different intensities caused by the different environment.The dominant band is predicted to be more red-shifted than the corresponding single peak found in the edge-sharing pentagonal prism, isomer A. Other families of (H 2 O) 20 isomers that have been discussed in the past, like those of the facesharing pentagonal prisms (isomer C in Fig. 1), the fused cubes, and the dodecahedron are predicted to be significantly less stable 28 (cf.Table 1).Note that the ordering of the isomers in the various families does not change upon either increasing the basis set to CBS or including harmonic ZPE corrections, as seen from Table 1.The small energy difference between isomers A and B, obtained at the highest level of electronic structure theory considered in this study including ZPE corrections, suggests that we expect only two different cluster networks to contribute to the IR spectrum in the OH stretching region of (H 2 O) 20 , each showing one characteristic feature in the region of DAA oscillators.3 the experimental IR spectrum generated in the Ne expansions is compared with a simulated one, in which the ratio of isomer A (35%) and B (65%) is chosen such that similar intensities for the characteristic DAA peaks are obtained.Note that the frequencies in all simulated spectra are scaled by 0.96 to account for the harmonic approximation.The summed-up spectrum agrees well with the experimentally observed absorption pattern over the complete OH stretching region.In the experiment the strong double-donor single-acceptor (DDA) modes appear at slightly lower frequencies (3400 cm À1 ) compared to the predicted harmonic spectra.Such deviations are expected due to the uniform scaling of the harmonic frequencies. 18The assignment   Fig. 1 The most stable isomers of (H 2 O) 20 , as identified in this and previous work 28,31 and their predicted IR spectra in the OH stretch region.(A) Edgesharing pentagonal prism, (B) drop-like, (C) face-sharing pentagonal prism.
of the DAA double peak to isomers A and B would further become more definitive upon depopulation of one of the isomers in the experiment.This was indeed achieved by changing the seed gas and stagnation pressure.The middle panel of Fig. 3 shows the spectrum obtained in Ar-seeded expansions, compared to the predicted IR spectrum of isomer B alone.Again, the experimentally observed OH stretch absorption pattern is consistent with the theoretical prediction thus providing further credence to the structural assignment.The DAA peak characteristic to isomer A is completely suppressed in this experiment and the structure sensitive DAA region agrees well with the predicted absorption pattern of isomer B. Both experiment and calculation consistently show the separation of DAA oscillators appearing below 3180 cm À1 and DDAA oscillators appearing above 3200 cm À1 .The high intensity of the low-lying DDAA oscillator is probably caused by a Fermi resonance with the overtone of the OH bending motion appearing around 3200 cm À1 , as recently shown by Suhm and co-workers for small water clusters in spontaneous Raman scattering experiments. 48Similar to the expansion with Ne as seed gas, the strong DDA features appear at 3400 cm À1 indicating that the calculated DDA part of the spectrum is a bit too compressed.The absence of isomer A in the Ar-seeded expansions could indicate a slightly higher stability of isomer B. However, also a kinetic effect seems plausible, which would indicate a better relaxation of clusters to the thermodynamic minima in the Ne experiment.The cluster temperatures are estimated to be around 100 K for the Ne-seeded expansions; the Ar-seeded experiment is expected to produce colder clusters (see discussion in ref. 39 for temperature estimates).The isolation of isomer B achieved in the Ar-seeded expansion allows for an additional consistency test of our spectral assignments, facilitated by subtracting the spectrum recorded using Ar as seed gas from the one taken in the Ne-seeded expansion.Ideally, this difference spectrum, shown in the lower panel of Fig. 3, should reveal the characteristic absorption pattern of isomer A: the single DAA peak around 3140 cm À1 and a broad intensity gap between this peak and the absorption maximum of the DDAA oscillators, predicted to be around 3300 cm À1 .This highly characteristic feature is indeed found in the difference spectrum, with the gap being larger than predicted from the harmonic frequencies.Again, the (anharmonic) effect of the Fermi resonance with the overtone of the OH bending motion might explain this observation.Here, the close lying DDAA oscillators seem to be blue-shifted, away from the resonance around 3200 cm À1 , as found for the in-phase, hydrogen-bonded OH stretching motion in the cyclic water pentamer. 48The consistency of the difference spectrum with the calculated one of isomer A alone suggests that indeed the isomer families A and B dominate the isomer composition of Na(H 2 O) 20 in the experiments.Since isomer A has a highly symmetric oxygen network, the experimental IR spectrum should be sensitive to the influence of our sodium atom tag.We analyzed this effect and found only a weak influence on the general IR absorption pattern for sodium doped versions of isomer A (see ESI †).However, satellite peaks around the single DAA peak are predicted.Such spiky features in the DAA region are found in the Ne seed gas experiment (and in the difference spectrum), in which isomer A is present (see Fig. 3).The simulations further indicate that DDAA oscillators are slightly blue-shifted in the Na tagged cluster.This effect would be an alternative explanation for the relatively large absorption gap between the DAA and DDAA features of isomer A and the missing gap between DDAA and DDA oscillators in all spectra.In the discussion above we have presented theoretical and experimental evidence that isomers A and B, the edge-sharing pentagonal prism and a drop-like structure, are the most stable configurations of (H 2 O) 20 .The assignment is anchored on the energy ranking based on high-level first principles electronic structure calculations (cf.Table 1) and the consistency between isomer selective, experimental IR spectra and theoretical predictions (Fig. 3).With regard to the ongoing tedious experimental and theoretical exploration of the (H 2 O) 6 structures, initiated over almost two decades ago, 8,11,14,15,49 the current study represents a significant step towards expanding the size regime of neutral water clusters that can be probed experimentally.The drop like cluster is the first experimentally detected cluster with a fully coordinated interior water molecule featuring the local environment in liquid water.This feature was theoretically expected to occur first in water clusters with an odd number of water molecules 17 starting with n = 17.The similarity of the spectral envelope for n = 19, 20, 21 (see Fig. 2) may indicate that the neighboring cluster sizes of n = 20 also feature large fractions of isomers with interior water molecules.The DAA region and the corresponding free OH stretch features are similar for n = 21 in Fig. 2 and the spectrum assigned to the drop like isomer of n = 20 in the middle panel of Fig. 3.We note that a centred dodecahedron was discussed as global minimum for n = 21 17 and that one half of the surface of isomer B shows this topology of interconnected five member rings.In future the approach presented in this work may allow for further structural assignments of water clusters.However, at the moment it is not clear whether one or two energetically distinguished isomer families as found for n = 20 exist for water clusters containing more than 20 water molecules. 39 Infrared detection of (H 2 O) 20 isomers of exceptional stability: a drop-like and a face-sharing pentagonal prism cluster † Christoph C. Pradzynski, a Christoph W. Dierking, a Florian Zurheide, a Richard M. Forck, a Udo Buck, b Thomas Zeuch* a and Sotiris S. Xantheas* c This journal is © the Owner Societies 2014 Phys.Chem.Chem.Phys., 2014, 16, 26691--26696 | 26693 In Fig. 2 the IR spectrum of the sodium tagged (H 2 O) 20 cluster is shown together with simultaneously measured ones of cluster sizes n = 19, 21.Here, we have to concentrate on the marked region of the most red shifted part of the spectrum of the DAA stretch motion.Only for the Na(H 2 O) 20 cluster we observe a characteristic double peak, which corresponds to the specific DAA features of isomers A and B. In the upper panel of Fig.

Fig. 3
Fig. 3 Comparison of experiment with theory.Maximum intensity normalized OH stretch spectra of Na(H 2 O) 20 recorded in Ne (upper panel, offset 0.4) and Ar (middle panel, offset 0.8) seeded expansions are compared with simulated spectra of a mixture of isomer A (35%) and isomer B (65%).Expansion conditions: water at 323 to 353 K (0.12 to 0.47 bar) in Ar at 1.0 to 1.7 bar stagnation pressure, see Fig. 2 for Ne-seeded expansions and supplement for more details.Photoionization was performed at 360 or 388 nm.