The solvation of two electrons in the gaseous clusters of Na- (NH3)n and Li- (NH3)n

Alkali metal ammonia clusters, in their cationic, neutral, and anionic form, are molecular models for the alkali-ammonia solutions, which have rich variation of phases with the solvated electrons playing an important role. With two s electrons, the Na(-)(NH(3))(n) and Li(-)(NH(3))(n) clusters are unique in that they capture the important aspect of the coupling between two solvated electrons. By first principles calculations, we demonstrate that the two electrons are detached from the metal by n = 10, which produces a cluster with a solvated electron pair in the vicinity of a solvated alkali cation. The coupling of the two electrons leads to either the singlet or triplet state, both of which are stable. They are also quite distinct from the hydrated anionic clusters Na(-)(H(2)O)(n) and Li(-)(H(2)O)(n), in that the solvated electrons are delocalized and widely distributed among the solvent ammonia molecules. The Na(-)(NH(3))(n) and Li(-)(NH(3))(n) series, therefore, provide another interesting type of molecular model for the investigation of solvated electron pairs.     

Comparison of hydrated hydroperoxide anion (HOO-)(H2O)n clusters with alkaline hydrogen peroxide (HOOH)(OH-)(H2O)(n-1) clusters, n = 1-8, 20: an ab initio study

Hydroperoxide anion (HOO(-)), the conjugate base of hydrogen peroxide (HOOH), has been relatively little studied despite the importance of HOOH in commercial processes, atmospheric science, and biology. The anion has been shown to exist as a stable species in alkaline water. This project explored the structure of gas phase (HOO(-))(H(2)O)(n) clusters and identified the lowest energy configurations for n ≤ 8 at the B3LYP/6-311++G** level of theory and for n ≤ 6 at the MP2/aug-cc-pVTZ level of theory. As a start toward understanding equilibration between HOO(-) and HOOH in an alkaline environment, (HOOH)(OH(-))(H(2)O)(n-1) clusters were likewise examined, and the lowest energy configurations were determined for n ≤ 8 (B3LYP/6-311++G**) and n ≤ 6 (MP2/aug-cc-pVTZ). Some studies were also done for n = 20. The two species have very different solvation behaviors. In low energy (HOOH)(OH(-))(H(2)O)(n-1) clusters, HOOH sits on the surface of the cluster, is 4-coordinated (each O is donor once and acceptor once), and donates to the hydroxide ion. In contrast, in low energy (HOO(-))(H(2)O)(n) clusters, (HOO(-)) takes a position in the cluster center surrounded on all sides by water molecules, and its optimum coordination number appears to be 7 (one O is donor-acceptor-acceptor while the other is a 4-fold acceptor). For n ≤ 6 the lowest (HOOH)(OH(-))(H(2)O)(n-1) cluster lies 1.0-2.1 kcal/mol below the lowest (HOO(-))(H(2)O)(n) cluster, but the lowest clusters found for n = 20 favor (HOO(-))(H(2)O)(20). The results suggest that ambient water could act as a substantial kinetic brake that slows equilibration between (HOOH)(OH(-)) and (HOO(-))(H(2)O) because extensive rearrangement of solvation shells is necessary to restabilize either species after proton transfer.     

First principles study on the solvation and structure of C2O4(2-)(H2O)n, n = 6-12

The structures and energies of hydrated oxalate clusters, C2O4(2-)(H2O)n, n = 6-12, are obtained by density functional theory (DFT) calculations and compared to SO4(2-)(H2O)n. Although the evolution of the cluster structure with size is similar to that of SO4(2-)(H2O)n, there are a number of important and distinctive futures in C2O4(2-)(H2O)n, including the separation of the two charges due to the C-C bond in C2O4(2-), the lower symmetry around C2O4(2-), and the torsion along the C-C bond, that affect both the structure and the solvation energy. The solvation dynamics for the isomers of C2O4(2-)(H2O)12 are also examined by DFT based ab initio molecular dynamics.     

The elimination of a hydrogen atom in Na(H2O)n

By a systematic examination on Na(H2O)n, with n = 4-7, 9, 10, and 15, we demonstrate that a hydrogen loss reaction can be initiated by a single sodium atom with water molecules. This reaction is similar to the well-known size-dependent intracluster hydrogen loss in Mg+(H2O)n, which is isoelectronic to Na(H2O)n. However, with one less charge on Na(H2O)n than that on Mg+(H2O)n, the hydrogen loss for Na(H2O)n is characterized by a higher barrier and a more flexible solvation shell around the metal ion, although the reaction should be accessible, as the lowest barrier is around 8 kcal/mol. Interestingly, the hydroxide ion OH- produced in the process is stabilized by the solvation of H2O molecules and the formation of an ion pair Na+(H2O)4(H2O)n-l-4[OH-(H2O)l]. The activation barrier is reduced as the unpaired electron in Na(H2O)n moves to higher solvation shells with increasing cluster size, and the reaction is not switched off for larger clusters. This is in sharp contrast to the reaction for Mg+(H2O)n, in which the OH- ion is stabilized by direct coordination with Mg2+ and the reaction is switched off for n > 17, as the unpaired electron moved to higher solvation shells. Such a contrast illustrates the important link between microsolvation environment and chemical reactivity in solvation clusters.     

The structures and electronic states of zinc-water clusters Zn(n)(H2O)(m) (n = 1-32 and m = 1-3)

Ab initio and Density Functional Theory (DFT) calculations have been carried out for zinc-water clusters Zn(n)-(H2O)(m) (n = 1-32 and m = 1-3, where n and m are the numbers of zinc atoms and water molecules, respectively) to elucidate the structure and electronic states of the clusters and the interaction of zinc cluster with water molecules. The binding energies of H2O to zinc clusters were small at n = 2-3 (2.3-4.2 kcal mol(-1)), whereas the energy increased significantly in n = 4 (9.0 kcal mol(-1)). Also, the binding nature of H2O was changed at n = 4. The cluster size dependency of the binding energy of H2O accorded well with that of the natural population of electrons in the 4p orbital of the zinc atom. In the larger clusters (n > 20), it was found that the zinc atoms in surface regions of the zinc cluster have a positive charge, whereas those in the interior region have a negative charge with the large electron population in the 4p orbital. The interaction of H2O with the zinc clusters were discussed on the basis of the theoretical results.     

In search of CS2(H2O)(n=1-4) clusters

Gaussian-3 and MP2/aug-cc-pVnZ methods have been used to calculate geometries and thermochemistry of CS(2)(H2O)n, where n=1-4. An extensive molecular dynamics search followed by optimization using these two methods located two dimers, six trimers, six tetramers, and two pentamers. The MP2/aug-cc-pVDZ structure matched best with the experimental result for the CS(2)(H2O) dimer, showing that diffuse functions are necessary to model the interactions found in this complex. For larger CS(2)(H2O)n clusters, the MP2/aug-cc-pVDZ minima are significantly different from the MP2(full)6-31G* structures, revealing that the G3 model chemistry is not suitable for investigation of sulfur containing van der Waals complexes. Based on the MP2/aug-cc-pVTZ free energies, the concentration of saturated water in the atmosphere and the average amount of CS(2) in the atmosphere, the concentrations of these clusters are predicted to be on the order of 10(5) CS(2)(H2O) clusters.cm(-3) and 10(2) CS(2)(H2O)(2) clusters.cm(-3) at 298.15 K. The MP2/aug-cc-pVDZ scaled harmonic and anharmonic frequencies of the most abundant dimer cluster at 298 K are presented, along with the MP2/aug-cc-pVDZ scaled harmonic frequencies for the CS(2)(H(2)O)(n) structures predicted to be present in a low-temperature molecular beam experiment.     

Ionization induced relaxation in solvation structure: a comparison between Na(H2O)n and Na(NH3)n

The constant ionization potential for hydrated sodium clusters Na(H2O)n just beyond n=4, as observed in photoionization experiments, has long been a puzzle in violation of the well-known (n+1)(-1/3) rule that governs the gradual transition in properties from clusters to the bulk. Based on first principles calculations, a link is identified between this puzzle and an important process in solution: the reorganization of the solvation structure after the removal of a charged particle. Na(H2O)n is a prototypical system with a solvated electron coexisting with a solvated sodium ion, and the cluster structure is determined by a balance among three factors: solute-solvent (Na+-H2O), solvent-solvent (H2O-H2O), and electron-solvent (OH{e}HO) interactions. Upon the removal of an electron by photoionization, extensive structural reorganization is induced to reorient OH{e}HO features in the neutral Na(H2O)n for better Na+-H2O and H2O-H2O interactions in the cationic Na+(H2O)n. The large amount of energy released, often reaching 1 eV or more, indicates that experimentally measured ion signals actually come from autoionization via vertical excitation to high Rydberg states below the vertical ionization potential, which induces extensive structural reorganization and the loss of a few solvent molecules. It provides a coherent explanation for all the peculiar features in the ionization experiments, not only for Na(H2O)n but also for Li(H2O)n and Cs(H2O)n. In addition, the contrast between Na(H2O)n and Na(NH3)n experiments is accounted for by the much smaller relaxation energy for Na(NH3)n, for which the structures and energetics are also elucidated.     

Ab initio molecular dynamics studies of ionic dissolution and precipitation of sodium chloride and silver chloride in water clusters, NaCl(H2O)n and AgCl(H2O)n, n = 6, 10, and 14

An ab initio molecular dynamics method was used to compare the ionic dissolution of soluble sodium chloride (NaCl) in water clusters with the highly insoluble silver chloride (AgCl). The investigations focused on the solvation structures, dynamics, and energetics of the contact ion pair (CIP) and of the solvent-separated ion pair (SSIP) in NaCl(H(2)O)(n) and AgCl(H(2)O)(n) with cluster sizes of n = 6, 10 and 14. We found that the minimum cluster size required to stabilize the SSIP configuration in NaCl(H(2)O)(n) is temperature-dependent. For n = 6, both configurations are present as two distinct local minima on the free-energy profile at 100 K, whereas SSIP is unstable at 300 K. Both configurations, separated by a low barrier (<10 kJ mol(-1)), are identifiable on the free energy profiles of NaCl(H(2)O)(n) for n = 10 and 14 at 300 K, with the Na(+)/Cl(-) pairs being internally solvated in the water cluster and the SSIP configuration being slightly higher in energy (<5 kJ mol(-1)). In agreement with the low bulk solubility of AgCl, no SSIP minimum is observed on the free-energy profiles of finite AgCl(H(2)O)(n) clusters. The AgCl interaction is more covalent in nature, and is less affected by the water solvent. Unlike NaCl, AgCl is mainly solvated on the surface in finite water clusters, and ionic dissolution requires a significant reorganization of the solvent structure.     

Infrared predissociation spectroscopy of M+ (C6H6)(1-4)(H2O)(1-2)Ar(0-1) cluster ions, M = Li, Na

Infrared predissociation (IRPD) spectra of Li(+)(C(6)H(6))(1-4)(H(2)O)(1-2)Ar(0-1) and Na(+)(C(6)H(6))(2-4)(H(2)O)(1-2)Ar(1) are presented along with ab initio calculations. The results indicate that the global minimum energy structure for Li(+)(C(6)H(6))(2)(H(2)O)(2) has each water forming a π-hydrogen bond with the same benzene molecule. This bonding motif is preserved in Li(+)(C(6)H(6))(3-4)(H(2)O)(2)Ar(0-1) with the additional benzene ligands binding to the available free OH groups. Argon tagging allows high-energy Li(+)(C(6)H(6))(2-4)(H(2)O)(2)Ar isomers containing water-water hydrogen bonds to be trapped and detected. The monohydrated, Li(+) containing clusters contain benzene-water interactions with varying strength as indicated by shifts in OH stretching frequencies. The IRPD spectra of M(+)(C(6)H(6))(1-4)(H(2)O)(1-2)Ar are very different for lithium-bearing versus sodium-bearing cluster ions emphasizing the important role of ion size in determining the most favorable balance of competing noncovalent interactions.     

Size-dependent charge-separation reaction for hydrated sulfate dianion cluster, SO(2-)(H2O)n, with n=3-7

The decrease in the reaction rate for the charge separation in SO(4) (2-)(H(2)O)(n) with increasing cluster size is examined by first-principles calculations of the energetics, activation barriers, and thermal stability for n=3-7. The key factor governing the charge separation is the difference in the strength of solvation interaction: while interaction with water is strong for the reactant SO(4) (2-) and the product OH(-), it is relatively weak for HSO(4) (-). It gives rise to a barrier for charge separation as SO(4) (2-) is transformed into HSO(4) (-) and OH(-), although the overall reaction energy is exothermic. The barrier is high when more than two H(2)O are left to solvate HSO(4) (-), as in the case of symmetric solvation structure and in the case of large clusters. The entropy is another important factor since the potential surface is floppy and the thermal motion facilitates the symmetric distribution of H(2)O around SO(4) (2-), which leads to the gradual reduction in reaction rate and the eventual switch-off of charge separation as cluster size increases. The experimentally observed products for n=3-5 are explained by the thermally most favorable isomer at each size, obtained by ab initio molecular-dynamics simulations rather than by the isomer with the lowest energy.     

Microsolvation of the dicyanamide anion: [N(CN)(2)(-)](H(2)O)n (n = 0-12)

Photoelectron spectroscopy is combined with ab initio calculations to study the microsolvation of the dicyanamide anion, N(CN)(2)(-). Photoelectron spectra of [N(CN)(2)(-)](H2O)n (n = 0-12) have been measured at room temperature and also at low temperature for n = 0-4. Vibrationally resolved photoelectron spectra are obtained for N(CN)(2)(-), allowing the electron affinity of the N(CN)2 radical to be determined accurately as 4.135 +/- 0.010 eV. The electron binding energies and the spectral width of the hydrated clusters are observed to increase with the number of water molecules. The first five waters are observed to provide significant stabilization to the solute, whereas the stabilization becomes weaker for n > 5. The spectral width, which carries information about the solvent reorganization upon electron detachment in [N(CN)(2)(-)](H2O)n, levels off for n > 6. Theoretical calculations reveal several close-lying isomers for n = 1 and 2 due to the fact that the N(CN)(2)(-) anion possesses three almost equivalent hydration sites. In all the hydrated clusters, the most stable structures consist of a water cluster solvating one end of the N(CN)(2)(-) anion.     

Reactions of hydrated electrons (H2O)n- with carbon dioxide and molecular oxygen: hydration of the CO2- and O2- ions

The gas-phase reactions of hydrated electrons with carbon dioxide and molecular oxygen were studied by Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry. Both CO2 and O2 react efficiently with (H2O)n- because they possess low-lying empty pi* orbitals. The molecular CO2- and O2- anions are concurrently solvated and stabilized by the water ligands to form CO2(-)(H2O)n and O2(-)(H2O)n. Core exchange reactions are also observed, in which CO2(-)(H2O)n is transformed into O2(-)(H2O)n upon collision with O2. This is in agreement with the prediction based on density functional theory calculations that O2(-)(H2O)n clusters are thermodynamically favored with respect to CO2(-)(H2O)n. Electron detachment from the product species is only observed for CO2(-)(H2O)2, in agreement with the calculated electron affinities and solvation energies.     

Isotope exchange in reactions between D2O and size-selected ionic water clusters containing pyridine, H+ (pyridine)m(H2O)n

Pyridine containing water clusters, H(+)(pyridine)(m)(H(2)O)(n), have been studied both experimentally by a quadrupole time-of-flight mass spectrometer and by quantum chemical calculations. In the experiments, H(+)(pyridine)(m)(H(2)O)(n) with m = 1-4 and n = 0-80 are observed. For the cluster distributions observed, there are no magic numbers, neither in the abundance spectra, nor in the evaporation spectra from size selected clusters. Experiments with size-selected clusters H(+)(pyridine)(m)(H(2)O)(n), with m = 0-3, reacting with D(2)O at a center-of-mass energy of 0.1 eV were also performed. The cross-sections for H/D isotope exchange depend mainly on the number of water molecules in the cluster and not on the number of pyridine molecules. Clusters having only one pyridine molecule undergo D(2)O/H(2)O ligand exchange, while H(+)(pyridine)(m)(H(2)O)(n), with m = 2, 3, exhibit significant H/D scrambling. These results are rationalized by quantum chemical calculations (B3LYP and MP2) for H(+)(pyridine)(1)(H(2)O)(n) and H(+)(pyridine)(2)(H(2)O)(n), with n = 1-6. In clusters containing one pyridine, the water molecules form an interconnected network of hydrogen bonds associated with the pyridinium ion via a single hydrogen bond. For clusters containing two pyridines, the two pyridine molecules are completely separated by the water molecules, with each pyridine being positioned diametrically opposite within the cluster. In agreement with experimental observations, these calculations suggest a "see-saw mechanism" for pendular proton transfer between the two pyridines in H(+)(pyridine)(2)(H(2)O)(n) clusters.     

Electron impact ionization of water-doped superfluid helium nanodroplets: observation of He(H(2)O)(n)(+) clusters

Electron impact mass spectra have been recorded for helium nanodroplets containing water clusters. In addition to identification of both H(+)(H(2)O)(n) and (H(2)O)(n)(+) ions in the gas phase, additional peaks are observed which are assigned to He(H(2)O)(n)(+) clusters for up to n=27. No clusters are detected with more than one helium atom attached. The interpretation of these findings is that quenching of (H(2)O)(n)(+) by the surrounding helium can cool the cluster to the point where not only is fragmentation to H(+)(H(2)O)(m) (where m < or = n-1) avoided, but also, in some cases, a helium atom can remain attached to the cluster ion as it escapes into the gas phase. Ab initio calculations suggest that the first step after ionization is the rapid formation of distinct H(3)O(+) and OH units within the (H(2)O)(n)(+) cluster. To explain the formation and survival of He(H(2)O)(n)(+) clusters through to detection, the H(3)O(+) is assumed to be located at the surface of the cluster with a dangling O-H bond to which a single helium atom can attach via a charge-induced dipole interaction. This study suggests that, like H(+)(H(2)O)(n) ions, the preferential location for the positive charge in large (H(2)O)(n)(+) clusters is on the surface rather than as a solvated ion in the interior of the cluster.     

Hydration process of Na2- in small water clusters: photoelectron spectroscopy and theoretical study of Na2- (H2O)n

Photoelectron spectroscopy (PES) of Na2- (H2O)n (n < or = 6) was investigated to examine the solvation of sodium aggregates in small water clusters. The PES bands for the transitions from the anion to the neutral ground and first excited states derived from Na2 (1(1)Sigmag+) and Na2 (1(3)Sigmau+) shifted gradually to the blue, and those to the higher-excited states correlated to the 3(2)S + 3(2)P asymptote dropped down rapidly to the red and almost degenerated on the 1(3)Sigmau+-type band at n = 4. Quantum chemical calculations for n up to 3 showed that the spectra can be ascribed to structures where one of the Na atoms is selectively hydrated. From the electron distributions, it is found that the Na- -Na+(H2O)n- -type electronic state grows with increasing cluster size, which can be regarded as a sign of the solvation of Na2- with ionization of the hydrated Na.     

Infrared spectroscopy of Cu+(H2O)(n) and Ag+(H2O)(n): coordination and solvation of noble-metal ions

M(+)(H(2)O)(n) and M(+)(H(2)O)(n)Ar ions (M=Cu and Ag) are studied for exploring coordination and solvation structures of noble-metal ions. These species are produced in a laser-vaporization cluster source and probed with infrared (IR) photodissociation spectroscopy in the OH-stretch region using a triple quadrupole mass spectrometer. Density functional theory calculations are also carried out for analyzing the experimental IR spectra. Partially resolved rotational structure observed in the spectrum of Ag(+)(H(2)O)(1) x Ar indicates that the complex is quasilinear in an Ar-Ag(+)-O configuration with the H atoms symmetrically displaced off axis. The spectra of the Ar-tagged M(+)(H(2)O)(2) are consistent with twofold coordination with a linear O-M(+)-O arrangement for these ions, which is stabilized by the s-d hybridization in M(+). Hydrogen bonding between H(2)O molecules is absent in Ag(+)(H(2)O)(3) x Ar but detected in Cu(+)(H(2)O)(3) x Ar through characteristic changes in the position and intensity of the OH-stretch transitions. The third H(2)O attaches directly to Ag(+) in a tricoordinated form, while it occupies a hydrogen-bonding site in the second shell of the dicoordinated Cu(+). The preference of the tricoordination is attributable to the inefficient 5s-4d hybridization in Ag(+), in contrast to the extensive 4s-3d hybridization in Cu(+) which retains the dicoordination. This is most likely because the s-d energy gap of Ag(+) is much larger than that of Cu(+). The fourth H(2)O occupies the second shells of the tricoordinated Ag(+) and the dicoordinated Cu(+), as extensive hydrogen bonding is observed in M(+)(H(2)O)(4) x Ar. Interestingly, the Ag(+)(H(2)O)(4) x Ar ions adopt not only the tricoordinated form but also the dicoordinated forms, which are absent in Ag(+)(H(2)O)(3) x Ar but revived at n=4. Size dependent variations in the spectra of Cu(+)(H(2)O)(n) for n=5-7 provide evidence for the completion of the second shell at n=6, where the dicoordinated Cu(+)(H(2)O)(2) subunit is surrounded by four H(2)O molecules. The gas-phase coordination number of Cu(+) is 2 and the resulting linearly coordinated structure acts as the core of further solvation processes.     

Effects of solvation and core switching on the photoelectron angular distributions from (CO2)n(-) and (CO2)(n)(-).H2O

Photoelectron images are recorded in the photodetachment of two series of cluster anions, (CO(2))(n)(-), n=4-9 and (CO(2))(n)(-).H(2)O, n=2-7, with linearly polarized 400 nm light. The energetics of the observed photodetachment bands compare well with previous studies, showing evidence for switching between two anionic core structures: The CO(2)(-) monomer and covalent (CO(2))(2)(-) dimer anions. The systematic study of photoelectron angular distributions (PADs) sheds light on the electronic structure of the different core anions and indicates that solvation by several CO(2) molecules and/or one water molecule has only moderate effect on the excess-electron orbitals. The observed PAD character is reconciled with the symmetry properties of the parent molecular orbitals. The most intriguing result concerns the PADs showing remarkable similarities between the monomer and dimer anion cluster-core types. This observation is explained by treating the highest-occupied molecular orbital of the covalent dimer anion as a linear combination of two spatially separated monomeric orbitals.     

Proton mobility and stability of water clusters containing the bisulfate anion, HSO4(-)(H2O)n

Bisulfate water clusters, HSO(4)(-)(H(2)O)(n), have been studied both experimentally by a quadrupole time-of-flight mass spectrometer and by quantum chemical calculations. For the cluster distributions studied, there are some possible "magic number" peaks, although the increase in abundance compared to their neighbours is small. Experiments with size-selected clusters with n = 0-25, reacting with D(2)O at a center-of-mass energy of 0.1 eV, were performed, and it was observed that the rate of hydrogen/deuterium exchange is lower for the smallest clusters (n < 8) than for the larger (n > 11), with a transition taking place in the range n = 8-11. We propose that the protonic defect of the bisulfate ion remains rather stationary unless the degree of hydration reaches a given level. In addition, it was observed that H/D scrambling becomes close to statistically randomized for the larger clusters. Insight into this size dependency was obtained by B3LYP/6-311++G(2d,2p) calculations for HSO(4)(-)(H(2)O)(n) with n = 0-10. In agreement with experimental observations, these calculations suggest pronounced effectiveness of a 'see-saw mechanism' for pendular proton transfer with increasing HSO(4)(-)(H(2)O)(n) cluster size.     

Photoelectron spectroscopic study of the hydrated nucleoside anions: Uridine(-)(H(2)O)(n=0-2), cytidine(-)(H(2)O)(n=0-2), and thymidine(-)(H(2)O)(n=0,1)

The hydrated nucleoside anions, uridine(-)(H(2)O)(n=0-2), cytidine(-)(H(2)O)(n=0-2), and thymidine(-)(H(2)O)(n=0,1), have been prepared in beams and studied by anion photoelectron spectroscopy in order to investigate the effects of a microhydrated environment on parent nucleoside anions. Vertical detachment energies (VDEs) were measured for all eight anions, and from these, estimates were made for five sequential anion hydration energies. Excellent agreement was found between our measured VDE value for thymidine(-)(H(2)O)(1) and its calculated value in the companion article by S. Kim and H. F. Schaefer III.