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.     

Undissociated versus dissociated structures for water clusters and ammonia-water clusters: (H2O)n and NH3(H2O)n-1 (n = 5, 8, 9, 21). Theoretical study

To understand the autoionization of pure water and the solvation of ammonia in water, we investigated the undissociated and dissociated (ion-pair) structures of (H2O) n and NH3(H2O)n-1 (n = 5, 8, 9, 21) using density functional theory (DFT) and second order Moller-Plesset perturbation theory (MP2). The stability, thermodynamic properties, and infrared spectra were also studied. The dissociated (ion-pair) form of the clusters tends to favor the solvent-separated ion-pair of H3O+/NH4+ and OH-. As for the NH3(H2O)20 cluster, the undissociated structure has the internal conformation, in contrast to the surface conformation for the (H2O)21 cluster, whereas the dissociated structure of NH3(H2O)20 has the surface conformation. As the cluster size of (H2O)n/NH3(H2O)n-1 increases, the difference in standard free energy between undissociated and dissociated (ion-pair) clusters is asymptotically well corroborated with the experimental free energy change at infinite dilution of H3O+/NH4+ and OH-. The predicted NH and OH stretching frequencies of the undissociated and dissociated (ion-pair) clusters are discussed.     

The identification of a solvated electron pair in the gaseous clusters of Na(-)(H2O)n and Li(-)(H2O)n

By first principles calculations, we explore the possibility that Na(-)(H(2)O)(n) and Li(-)(H(2)O)(n) clusters, which have been measured previously by photoelectron experiments, could serve as gas-phase molecular models for the solvation of two electrons. Such models would capture the electron-electron interaction in a solution environment, which is missed in the well-known anionic water clusters (H(2)O)(n) (-). Our results show that by n = 10, the two loosely bound s electrons in Li(-)(H(2)O)(n) are indeed detached from lithium, and they could exist in either the singlet (spin-paring) or the triplet (spin-coupling) state. In contrast, the two electrons would prefer to stay on the sodium atom in Na(-)(H(2)O)(n) and on the surface of the cluster. The formation of a solvated electron pair and the variation in solvation structures make these two cluster series interesting subjects for further experimental investigation.     

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.     

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.     

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.     

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.     

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.     

Density functional theory studies of inorganic metallocene multidecker V(n)(P(6))(n+1) (n=1-4) sandwich clusters

Motivated by the synthesis of the first entirely inorganic metallocene sandwich ion [eta(5)-Ti-(P(5))(2)](2-) [E. Urnezius et al. Science 295, 832 (2002)], we have designed a new inorganic metallocene sandwich [eta(6)-V-(P(6))(2)] and multidecker sandwich clusters up to V(4)(P(6))(5) by employing an all electron gradient-corrected density functional theory. The binding energies of the V(n)(P(6))(n+1) complexes increase rapidly from half sandwich to the smallest full sandwich and become gradually afterwards. The highest occupied and lowest unoccupied molecular orbital gap and the vertical ionization energy decrease with increasing cluster size. The V(n)(P(6))(n+1) clusters are nonferromagnetic and prefer the lowest available spin states. The smallest sandwich cluster, V(P(6))(2), has the high stability and might serve as a building block for one-dimensional inorganic polymers with high stabilities.     

Ab initio studies on Al(+)(H(2)O)(n), HAlOH(+)(H(2)O)(n-1), and the size-dependent H(2) elimination reaction

We report computational studies on Al(+)(H(2)O)(n), and HAlOH(+)(H(2)O)(n-1), n = 6-14, by the density functional theory based ab initio molecular dynamics method, employing a planewave basis set with pseudopotentials, and also by conventional methods with Gaussian basis sets. The mechanism for the intracluster H(2) elimination reaction is explored. First, a new size-dependent insertion reaction for the transformation of Al(+)(H(2)O)(n), into HAlOH(+)(H(2)O)(n-1) is discovered for n > or = 8. This is because of the presence of a fairly stable six-water-ring structure in Al(+)(H(2)O)(n) with 12 members, including the Al(+). This structure promotes acidic dissociation and, for n > or = 8, leads to the insertion reaction. Gaussian based BPW91 and MP2 calculations with 6-31G* and 6-31G** basis sets confirmed the existence of such structures and located the transition structures for the insertion reaction. The calculated transition barrier is 10.0 kcal/mol for n = 9 and 7.1 kcal/mol for n = 8 at the MP2/6-31G** level, with zero-point energy corrections. Second, the experimentally observed size-dependent H(2) elimination reaction is related to the conformation of HAlOH(+)(H(2)O)(n-1), instead of Al(+)(H(2)O)(n). As n increases from 6 to 14, the structure of the HAlOH(+)(H(2)O)(n-1) cluster changes into a caged structure, with the Al-H bond buried inside, and protons produced in acidic dissociation could then travel through the H(2)O network to the vicinity of the Al-H bond and react with the hydride H to produce H(2). The structural transformation is completed at n = 13, coincident approximately with the onset of the H(2) elimination reaction. From constrained ab initio MD simulations, we estimated the free energy barrier for the H(2) elimination reaction to be 0.7 eV (16 kcal/mol) at n = 13, 1.5 eV (35 kcal/mol) at n = 12, and 4.5 eV (100 kcal/mol) at n = 8. The existence of transition structures for the H(2) elimination has also been verified by ab initio calculations at the MP2/6-31G** level. Finally, the switch-off of the H(2) elimination for n > 24 is explored and attributed to the diffusion of protons through enlarged hydrogen bonded H(2)O networks, which reduces the probability of finding a proton near the Al-H bond.     

Dielectron attachment and hydrogen evolution reaction in water clusters

Binding of excess electrons to nanosize water droplets, with a focus on the hitherto largely unexplored properties of doubly-charged clusters, were investigated experimentally using mass spectrometry and theoretically with large-scale first-principles simulations based on spin-density-functional theory, with all the valence electrons (that is, 8e per water molecule) and excess electrons treated quantum mechanically. Singly-charged clusters (H(2)O)(n)(-1) were detected for n = 6-250, and our calculated vertical detachment energies agree with previously measured values in the entire range 15 ≤ n ≤ 105, giving a consistent interpretation in terms of internal, surface and diffuse states of the excess electron. Doubly-charged clusters were measured in the range of 83 ≤ n ≤ 123, with (H(2)O)(n)(-2) clusters found for 83 ≤ n < 105, and mass-shifted peaks corresponding to (H(2)O)(n-2)(OH(-))(2) detected for n ≥ 105. The simulations revealed surface and internal dielectron, e(-)(2), localization modes and elucidated the mechanism of the reaction (H(2)O)(n)(-2) → (H(2)O)(n-2) (OH(-))(2) + H(2) (for n ≥ 105), which was found to occur via concerted approach of a pair of protons belonging to two water molecules located in the first shell of the dielectron internal hydration cavity, culminating in formation of a hydrogen molecule 2H(+) + e(-)(2) → H(2). Instability of the dielectron internal localization impedes the reaction for smaller (n < 105) doubly-charged clusters.     

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.     

Sequential hydration energies of the sulfate ion, from determinations of the equilibrium constants for the gas-phase reactions: SO4(H2O)(n)2- = SO4(H2O)(n-1)2- + H2O

Sequential hydration energies of SO4(H2O)(n)2- were obtained from determinations of the equilibrium constants of the following reactions: SO4(H2O)(n)2- = SO4(H2O)(n-1)2- + H2O. The SO4(2-) ions were produced by electrospray and the equilibrium constants Kn,n-1 were determined with a reaction chamber attached to a mass spectrometer. Determinations of Kn,n-1 at different temperatures were used to obtain DeltaG0n,n-1, DeltaH0 n,n-1, and DeltaS0n,n-1 for n = 7 to 19. Interference of the charge separation reaction SO4(H2O)(n)2- = HSO4(H2O)(n-k)- + OH(H2O)(k-1)- at higher temperatures prevented determinations for n < 7. The DeltaS0n,n-1 values obtained are unusually low and this indicates very loose, disordered structures for the n > or = 7 hydrates. The DeltaH0n,n-1 values are compared with theoretical values DeltaEn,n-1, obtained by Wang, Nicholas, and Wang. Rate constant determinations of the dissociation reactions n,n - 1, obtained with the BIRD method by Wong and Williams, showed relatively lower rates for n = 6 and 12, which indicate that these hydrates are more stable. No discontinuities of the DeltaG0n,n-1 values indicating an unusually stable n = 12 hydrate were observed in the present work. Rate constants evaluated from the DeltaG0n,n-1 results also fail to indicate a lower rate for n = 12. An analysis of the conditions used in the two types of experiments indicates that the different results reflect the different energy distributions expected at the dissociation threshold. Higher internal energies prevail in the equilibrium measurements and allow the participation of more disordered transition states in the reaction.     

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.     

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.     

Thermal effects on energetics and dynamics in water cluster anions (H2O)n(-)

The electron binding energies and relaxation dynamics of water cluster anions (H(2)O)(n)(-) (11 ≤ n ≤ 80) formed in co-expansions with neon were investigated using one-photon and time-resolved photoelectron imaging. Unlike previous experiments with argon, water cluster anions exhibit only one isomer class, the tightly bound isomer I with approximately the same binding energy as clusters formed in argon. This result, along with a decrease in the internal conversion lifetime of excited (H(2)O)(n)(-) (25 ≤ n ≤ 40), indicates that clusters are vibrationally warmer when formed in neon. Over the ranges studied, the vertical detachment energies and lifetimes appear to converge to previously reported values.     

Single photon ionization of van der Waals clusters with a soft x-ray laser: (SO2)n and (SO2)n(H2O)m

van der Waals cluster (SO2)n is investigated by using single photon ionization of a 26.5 eV soft x-ray laser. During the ionization process, neutral clusters suffer a small fragmentation because almost all energy is taken away by the photoelectron and a small part of the photon energy is deposited into the (SO2)n cluster. The distribution of (SO2)n clusters decreases roughly exponentially with increasing cluster size. The photoionization dissociation fraction of I[(SO2)(n-1)SO+] / I[(SO2)n+] decreases with increasing cluster size due to the formation of cluster. The metastable dissociation rate constants of (SO2)n+ are measured in the range of (0.6-1.5) x 10(4) s(-1) for cluster sizes 5< or =n< or =16. Mixed SO2-H2O clusters are studied at different experimental conditions. At the condition of high SO2 concentration (20% SO2 partial pressure), (SO2)n+ cluster ions dominate the mass spectrum, and the unprotonated mixed cluster ions (SO2)nH2O+ (1< or =n< or =5) are observed. At the condition of low SO2 concentration (5% SO2 partial pressure) (H2O)nH+ cluster ions are the dominant signals, and protonated cluster ions (SO2)(H2O)nH+ are observed. The mixed clusters, containing only one SO2 or H2O molecule, SO2(H2O)nH+ and (SO2)nH2O+ are observed, respectively.     

Polymer imprinting with iron-oxo-hydroxo clusters: [Fe6O2(OH)2(O2CC(Cl)=CH2)12(H2O)2], [Fe6O2(OH)2(O2C-Ph-(CH)=CH2)12(H2O)2] and [{Fe(O2CC(Cl)=CH2)(OMe)2}10

We report the syntheses of imprinted polymers using iron-oxo-hydroxo clusters as templates. Three new iron clusters, [Fe(6)O(2)(OH)(2)(O(2)CC(Cl)=CH(2))(12)(H(2)O)(2)] (1), [{Fe(O(2)CC(Cl)=CH(2))(OMe)(2)}(10)] (2) and [Fe(6)O(2)(OH)(2)(O(2)C-Ph-(CH)=CH(2))(12)(H(2)O)(2)] (3) have been prepared from commercially-available carboxylic acids. Cluster-imprinted-polymers (CIPs) of 1, 2 and 3 were prepared with ethylene glycol dimethacrylate monomer, and of 1 with methyl methacrylate monomer. The imprinted sites within the CIPs were examined using EXAFS and diffuse reflectance UV/vis spectroscopy, demonstrating that the clusters 1, 2 and 3 were incorporated intact within the polymers. Extraction of the clusters from the CIPs imprinted with 1 and 3 gave new polymers that showed evidence of an imprinting effect.     

Exploring the rich energy landscape of sulfate-water clusters SO4(2-) (H2O)(n=3-7): an electronic structure approach

We present a reinvestigation of sulfate-water clusters SO4(2-) (H2O)(n=3-7), which involves several new aspects. Using a joint molecular mechanics/first principles approach, we perform exhaustive searches for stable cluster geometries, showing that the sulfate-water landscape is much richer than anticipated previously. We check the compatibility of the new structures with experiment by comparing vertical detachment energies (VDEs) calculated at the B3LYP/6-311++G** level of theory and determine the energetic ordering of the isomers at the RI-MP2/aug-cc-pVTZ level. Our results are bench-marked carefully against reference energies of estimated CCSD(T)/aug-cc-VTZ quality and VDEs of CCSD(T)/aug-cc-pVDZ quality. Furthermore, we calculate anharmonic vibrational corrections for up to the n = 6 clusters, which are shown to be significant for isomer energy ordering. We use energy decomposition analysis (EDA) based on the absolutely localized fragment (ALMO) expansion to gain chemical insight into the binding motifs.     

Tuning of the internal energy and isomer distribution in small protonated water clusters H(+)(H2O)(4-8): an application of the inert gas messenger technique

Infrared spectroscopy of gas-phase hydrated clusters provides us much information on structures and dynamics of water networks. However, interpretation of spectra is often difficult because of high internal energy (vibrational temperature) of clusters and coexistence of many isomers. Here we report an approach to vary these factors by using the inert gas (so-called "messenger")-mediated cooling technique. Protonated water clusters with a messenger (M), H(+)(H(2)O)(4-8)·M (M = Ne, Ar, (H(2))(2)), are formed in a molecular beam and probed with infrared photodissociation spectroscopy in the OH stretch region. Observed spectra are compared with each other and with bare H(+)(H(2)O)(n). They show clear messenger dependence in their bandwidths and relative band intensities, reflecting different internal energy and isomer distribution, respectively. It is shown that the internal energy follows the order H(+)(H(2)O)(n) > H(+)(H(2)O)(n)·(H(2))(2) > H(+)(H(2)O)(n)·Ar > H(+)(H(2)O)(n)·Ne, while the isomer-selectivity, which changes the isomer distribution in the bare system, follows the order H(+)(H(2)O)(n)·Ar > H(+)(H(2)O)(n)·(H(2))(2) > H(+)(H(2)O)(n)·Ne ~ (H(+)(H(2)O)(n)). Although the origin of the isomer-selectivity is unclear, comparison among spectra measured with different messengers is very powerful in spectral analyses and makes it possible to easily assign spectral features of each isomer.