Vibrational spectroscopy of small water complexes embedded in large liquid helium clusters

Revisiting the Spectroscopy of Water Dimer in Jets

Laser spectroscopy in jets is one of the main sources of structural data from molecular aggregates. Consequently, numerous and sophisticated experimental systems have been developed to extract precise information, which is usually interpreted in the light of quantum mechanical calculations. However, even with the most sophisticated experiments, it is sometimes difficult to interpret the experimental results. We present here the example of water dimer and how after almost 70 years, the assignment of its mass-resolved IR spectrum still generates controversy that extends toward the mechanism of ionization of water aggregates.

Water Radical Cations in the Gas Phase: Methods and Mechanisms of Formation, Structure and Chemical Properties

Water radical cations, (H₂O)_n^+•, are of great research interest in both fundamental and applied sciences. Fundamental studies of water radical reactions are important to better understand the mechanisms of natural processes, such as proton transfer in aqueous solutions, the formation of hydrogen bonds and DNA damage, as well as for the discovery of new gas-phase reactions and products. In applied science, the interest in water radicals is prompted by their potential in radiobiology and as a source of primary ions for selective and sensitive chemical ionization. However, in contrast to protonated water clusters, (H₂O)_nH⁺, which are relatively easy to generate and isolate in experiments, the generation and isolation of radical water clusters, (H₂O)_n^+•, is tremendously difficult due to their ultra-high reactivity. This review focuses on the current knowledge and unknowns regarding (H₂O)_n^+• species, including the methods and mechanisms of their formation, structure and chemical properties.

Spectra of the D2O dimer in the O-D fundamental stretch region: Vibrational dependence of tunneling splittings and lifetimes

The fundamental O-D stretch region (2600-2800 cm^-1) of the fully deuterated water dimer (D₂O)₂ is studied using a pulsed supersonic slit jet source and a tunable optical parametric oscillator source. Relatively high spectral resolution (0.002 cm^-1) enables all six dimer tunneling components to be observed, in most cases, for the acceptor asymmetric O-D stretch, the donor free O-D stretch, and the donor bound O-D stretch vibrations. The dominant acceptor switching tunneling splittings are observed to decrease moderately in the excited O-D stretch states, to roughly 75% of their ground state values, whereas the smaller donor-acceptor interchange splittings show more dramatic and irregular decreases. Excited state predissociation lifetimes, as determined from the observed line broadening, show large variations (0.2 ≤ τ ≤ 5 ns) depending on the vibrational state, K-value, and tunneling symmetry. Another very weak band is tentatively assigned to a combination mode involving an intramolecular O-D stretch plus an intermolecular twist overtone. Asymmetric O-D stretch bands of the mixed isotopologue dimers D₂O-DOH and D₂O-HOD are also observed and analyzed.

Infrared spectroscopy of the ν2 band of the water monomer and small water clusters (H2O)n=2,3,4 in helium droplets

Ro-vibrational transitions of water and water clusters (H₂O)_n=1,2,3,4 for the ν₂ bending vibration were observed and assigned to distinct structures.

Communication: Electron impact ionization of binary H2O∕X clusters in helium nanodroplets: an ab initio perspective

In a recent experiment (H(2)O)(n)∕X(m) binary clusters (where X = Ar, N(2), CO, CO(2), and several other molecules) were formed in superfluid helium nanodroplets and investigated by electron impact mass spectrometry [Liu et al., Phys. Chem. Chem. Phys. 13, 13920 (2011)]. The addition of dopant X was found to affect the branching ratio between H(3)O(+)(H(2)O)(n) and (H(2)O)(+)(n+2) formation. Specifically, the addition of CO increased the proportion of protonated water cluster ions, whereas dopants such as Ar, N(2), and CO(2), had the opposite effect. In this work ab initio calculations have been performed on [X(H(2)O)(2)](+) ions, where X = Ar, N(2), CO, and CO(2), to try and explain this distinct behavior. CO is found to be unique in that it forms a HOCO-H(3)O(+) unit in the most stable cationic complexes where the binding between HO and CO is stronger than that between H(3)O(+) and OH. Thus, on purely energetic grounds, loss of HOCO rather than CO should be the preferred fragmentation process. No comparable chemistry occurs when X = Ar, N(2), or CO(2) and so the co-dopant requires less energy to depart than OH. The calculations therefore account for the experimental observations and provide evidence that HOCO formation is induced in helium droplets containing (H(2)O)(n) clusters and co-doped with CO when subject to electron impact ionization.

A threshold-based approach to calorimetry in helium droplets: measurement of binding energies of water clusters

Helium droplet beam methods have emerged as a versatile technique that can be used to assemble a wide variety of atomic and molecular clusters. We have developed a method to measure the binding energies of clusters assembled in helium droplets by determining the minimum droplet sizes required to assemble and detect selected clusters in the spectrum of the doped droplet beam. The differences in the droplet sizes required between the various multimers are then used to estimate the incremental binding energies. We have applied this method to measure the binding energies of cyclic water clusters from the dimer to the tetramer. We obtain measured values of D(0) that are in agreement with theoretical estimates to within ∼20%. Our results suggest that this threshold-based approach should be generally applicable using either mass spectrometry or optical spectroscopy techniques for detection, provided that the clusters selected for study are at least as strongly bound as those of water, and that a peak in the overall spectrum of the beam corresponding only to the cluster chosen (at least in the vicinity of the threshold) can be located.

Ionization of doped helium nanodroplets: complexes of C60 with water clusters

Water clusters are known to undergo an autoprotonation reaction upon ionization by photons or electron impact, resulting in the formation of (H(2)O)(n)H(3)O(+). Ejection of OH cannot be quenched by near-threshold ionization; it is only partly quenched when clusters are complexed with inert gas atoms. Mass spectra recorded by electron ionization of water-doped helium droplets show that the helium matrix also fails to quench OH loss. The situation changes drastically when helium droplets are codoped with C(60). Charged C(60)-water complexes are predominantly unprotonated; C(60)(H(2)O)(4)(+) and (C(60))(2)(H(2)O)(4)(+) appear with enhanced abundance. Another intense ion series is due to C(60)(H(2)O)(n)OH(+); dehydrogenation is proposed to be initiated by charge transfer between the primary He(+) ion and C(60). The resulting electronically excited C(60)(+*) leads to the formation of a doubly charged C(60)-water complex either via emission of an Auger electron from C(60)(+*), or internal Penning ionization of the attached water complex, followed by charge separation within {C(60)(H(2)O)(n)}(2+). This mechanism would also explain previous observations of dehydrogenation reactions in doped helium droplets. Mass-analyzed ion kinetic energy scans reveal spontaneous (unimolecular) dissociation of C(60)(H(2)O)(n)(+). In addition to the loss of single water molecules, a prominent reaction channel yields bare C(60)(+) for sizes n=3, 4, or 6. Ab initio Hartree-Fock calculations for C(60)-water complexes reveal negligible charge transfer within neutral complexes. Cationic complexes are well described as water clusters weakly bound to C(60)(+). For n=3, 4, or 6, fissionlike desorption of the entire water complex from C(60)(H(2)O)(n)(+) energetically competes with the evaporation of a single water molecule.

Formation of even-numbered hydrogen cluster cations in ultracold helium droplets

Neutral hydrogen clusters are grown in ultracold helium nanodroplets by successive pickup of hydrogen molecules. Even-numbered hydrogen cluster cations are observed upon electron-impact ionization with and without attached helium atoms and in addition to the familiar odd-numbered H(n)(+). The helium matrix affects the fragmentation dynamics that usually lead to the formation of overwhelmingly odd-numbered H(n)(+). The use of high-resolution mass spectrometry allows the unambiguous identification of even-numbered H(n)(+) up to n approximately = 120 by their mass excess that distinguishes them from He(n)(+), mixed He(m)H(n)(+), and background ions. The large range in size of these hydrogen cluster ions is unprecedented, as is the accuracy of their definition. Apart from the previously observed magic number n=6, pronounced drops in the abundance of even-numbered cluster ions are seen at n=30 and 114, which suggest icosahedral shell closures at H(6)(+)(H(2))(12) and H(6)(+)(H(2))(54). Possible isomers of H(6)(+) are identified at the quadratic configuration interaction with inclusion of single and double excitations (QCISD)/aug-cc-pVTZ level of theory.

Blueshift and intramolecular tunneling of NH3 umbrella mode in 4He n clusters

We present diffusion Monte Carlo calculations of the ground and first excited vibrational states of NH(3) (4)He(n) for n< or =40. We use the potential energy surface developed by one of us [M. P. Hodges and R. J. Wheatley, J. Chem. Phys. 114, 8836 (2001)], which includes the umbrella mode coordinate of NH(3). Using quantum Monte Carlo calculations of excited states, we show that this potential is able to reproduce qualitatively the experimentally observed effects of the helium environment, namely, a blueshift of the umbrella mode frequency and a reduction of the tunneling splittings in ground and first excited vibrational states of the molecule. These basic features are found to result regardless of whether dynamical approximations or exact calculations are employed.

Close-coupling study of rotational energy transfer and differential scattering in H2O collisions with He atoms

Quantum close-coupling scattering calculations of rotational energy transfer (RET) of rotationally excited H(2)O due to collisions with He are presented for collision energies between 10(-6) and 1000 cm(-1) with para-H(2)O initially in levels 1(1,1), 2(0,2), 2(1,1), and 2(2,0) and ortho-H(2)O in levels 1(1,0), 2(1,2), and 2(2,1). Quenching cross sections and rate coefficients for state-to-state RET were computed. Both elastic and inelastic differential cross sections are also calculated and compared with relative experimental results giving generally good agreement in all cases, but less so for inelastic results. Significant differences in the computed collisional parameters, obtained on three different potential energy surfaces (PESs), were found particularly in the ultracold regime. In the thermal regime, the rate coefficients calculated on each of the surfaces are generally in better agreement and comparable, but typically larger, than those obtained in a previous calculation. Unfortunately, a lack of absolute differential or integral inelastic experimental data prevents firm determination of a preferred PES.

Infrared intensity in small ammonia and water clusters

Helium droplet technique has been used in order to measure the strength of the infrared absorption in small ammonia and water clusters as a function of size. Hydrogen bonding in ammonia and water dimers causes an enhancement of the intensity of the hydrogen stretching bands by a factor of four and three, respectively. Two types of the hydrogen bonded clusters show different size dependence of the infrared intensity per hydrogen bond. In ammonia (NH3)2 and (NH3)3 it is close to the crystal value. In water clusters, it increases monotonically with cluster size being in tetramers, a factor of two smaller than in the ice. The measured infrared intensity in water clusters is found to be a factor of two to three smaller as compared to the results of numerical calculations.

Prediction of accurate anharmonic experimental vibrational frequencies for water clusters, (H2O)n, n=2-5

Accurate anharmonic experimental vibrational frequencies for water clusters consisting of 2-5 water molecules have been predicted on the basis of comparing different methods with MP2/aug-cc-pVTZ calculated and experimental anharmonic frequencies. The combination of using HF/6-31G* scaled frequencies for intramolecular modes and anharmonic frequencies for intermolecular modes gives excellent agreement with experiment for the water dimer and trimer and are as good as the expensive anharmonic MP2 calculations. The water trimer, the cyclic Ci and S4 tetramers, and the cyclic pentamer all have unique peaks in the infrared spectrum between 500 and 800 cm-1 and between 3400 and 3700 cm-1. Under the right experimental conditions these different clusters can be uniquely identified using high-resolution IR spectroscopy.

(HCl)2 and (HF)2 in small helium clusters: Quantum solvation of hydrogen-bonded dimers

We present a rigorous theoretical study of the solvation of (HCl)(2) and (HF)(2) by small ((4)He)(n) clusters, with n=1-14 and 30. Pairwise-additive potential-energy surfaces of He(n)(HX)(2) (X=Cl and F) clusters are constructed from highly accurate four-dimensional (rigid monomer) HX-HX and two-dimensional (rigid monomer) He-HX potentials and a one-dimensional He-He potential. The minimum-energy geometries of these clusters, for n=1-6 in the case of (HCl)(2) and n=1-5 for (HF)(2), correspond to the He atoms in a ring perpendicular to and bisecting the HX-HX axis. The quantum-mechanical ground-state energies and vibrationally averaged structures of He(n)(HCl)(2) (n=1-14 and 30) and He(n)(HF)(2) (n=1-10) clusters are calculated exactly using the diffusion Monte Carlo (DMC) method. In addition, the interchange-tunneling splittings of He(n)(HCl)(2) clusters with n=1-14 are determined using the fixed-node DMC approach, which was employed by us previously to calculate the tunneling splittings for He(n)(HF)(2) clusters, n=1-10 [A. Sarsa et al., Phys. Rev. Lett. 88, 123401 (2002)]. The vibrationally averaged structures of He(n)(HX)(2) clusters with n=1-6 for (HCl)(2) and n=1-5 for (HF)(2) have the helium density localized in an effectively one-dimensional ring, or doughnut, perpendicular to and at the midpoint of the HX-HX axis. The rigidity of the solvent ring varies with n and reaches its maximum for the cluster size at which the ring is filled, n=6 and n=5 for (HCl)(2) and (HF)(2), respectively. Once the equatorial ring is full, the helium density spreads along the HX-HX axis, eventually solvating the entire HX dimer. The interchange-tunneling splitting of He(n)(HCl)(2) clusters hardly varies at all over the cluster size range considered, n=1-14, and is virtually identical to that of the free HCl dimer. This absence of the solvent effect is in sharp contrast with our earlier results for He(n)(HF)(2) clusters, which show a approximately 30% reduction of the tunneling splitting for n=4. A tentative explanation for this difference is proposed. The implications of our results for the interchange-tunneling dynamics of (HCl)(2) in helium nanodroplets are discussed.

Infrared spectra of water clusters in krypton and xenon matrices

The infrared absorption spectra of the water molecules and small water clusters, (H(2)O)(n) with n = 2-6, trapped in solid argon, krypton, and xenon matrices have been investigated. The infrared bands of the water clusters with n = 5 and 6 in krypton and n = 3, 4, 5, and 6 in xenon matrices have been identified for the first time in the bonded OH stretching region. The frequency shifts in the bonded OH stretching band of the water dimer and trimer in xenon matrices show fairly large deviations to the red from the empirical correlation between the matrix shifts and the square root of the critical temperatures of the matrix material. The observed anomalous shifts suggest that the water dimer and trimer in solid xenon are trapped in multiple sites, and that the structures of the preferential trapping sites are different from those in argon and krypton matrices.