Ab initio studies on the mechanism of the size-dependent hydrogen-loss reaction in Mg+(H2O)n

The mechanism of size-dependent intracluster hydrogen loss in the cluster ions Mg(+)(H(2)O)(n), which is switched on around n=6, and off around n=14, was studied by ab initio calculations at the MP2/6-31G* and MP2/6-31G** levels for n=1-6. The reaction proceeds by Mg(+)-assisted breaking of an H-O bond in one of the H(2)O molecules. The reaction barrier is dependent on both the cluster size and the solvation structure. As n increases from 1 to 6, there is a dramatic drop in the reaction barrier, from greater than 70 kcal mol(-1) for n=1 to less than 10 kcal mol(-1) for n=6. In the transition structures, the Mg atom is close to the oxidation state of +2, and H(2)O molecules in the first solvation shell are much more effective in stabilizing the transition structures and lowering the reaction barriers than H(2)O molecules in the other solvation shells. While the reaction barrier for trimer core structures with only three H(2)O molecules in the first shell is greater than 24 kcal mol(-1), even for Mg(+)(H(2)O)(6), it drops considerably for clusters with four-six H(2)O molecules in the first shell. The more highly coordinated complexes have comparable or slightly higher energy than the trimer core structures, and the presence of such high coordination number complexes is the underlying kinetic factor for the switching on of the hydrogen-loss reaction around n=6. For clusters with trimer core structures, the hydrogen loss reaction is much easier when it is preceded by an isomerization step that increases the coordination number around Mg(+). Delocalization of the electron on the singly occupied molecular orbital (SOMO) away from the Mg(+) ion is observed for the hexamer core structure, while at the same time this isomer is the most reactive for the hydrogen-loss reaction, with an energy barrier of only 2.7 kcal mol(-1) at the MP2/6-31G** level.     

Ab initio study of the structures and stabilities of the dimer of ethyl cation, (C(2)H(+)(5))(2) and related C(4)H(10)(2+) isomers.

Ab initio calculations at the MP4(SDTQ)/6-311G//MP2/6-31G level were performed to study the structures and stabilities of the dimer of ethyl cation, (C(2)H(+)(5))(2), and related C(4)H(10)(2+) isomers. Two doubly hydrogen bridged diborane type trans 1 and cis 2 isomers were located as minima. The trans isomer was found to be more favorable than cis isomer by only 0.6 kcal/mol. Several other minima for C(4)H(10)(2+) were also located. However, the global energy minimum corresponds to C-H (C(4) position) protonated 2-butyl cation 10. Structure 10 was computed to be substantially more stable than 1 by 31.7 kcal/mol. The structure 10 was found to be lower in energy than 2-butyl cation 13 by 34.4 kcal/mol.     

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.     

The structure of gas-phase [Al·nH2O]+: hydrated monovalent aluminium Al+ (H2O)n or hydride-hydroxide HAlOH+ (H2O)(n-1)?

Theoretical studies predict that [Al·nH(2)O](+) clusters are present as hydride-hydroxide species HAlOH(+)(H(2)O)(n-1) in gas-phase experiments, energetically favoured by 200 kJ mol(-1) over Al(+)(H(2)O)(n). After collisions with D(2)O, however, no H/D scrambling occurs between H(2)O and D(2)O in clusters with n > 38, indicating that large clusters are present as the higher-energy isomers Al(+)(H(2)O)(n).     

Isotope exchange and structural rearrangements in reactions between size-selected ionic water clusters, H3O(+)(H2O)(n) and NH4(+)(H2O)(n), and D2O

Hydrogen/deuterium exchange in reactions of H3O(+)(H2O)n and NH4(+)(H2O)n (1 < or = n < or = 30) with D2O has been studied experimentally at center-of-mass collisions energies of < or = 0.2 eV. For a given cluster size, the cross-sections for H3O(+)(H2O)n and NH4(+)(H2O)n are similar, indicating a structural resemblance and energetics of binding. For protonated pure water clusters, H3O(+)(H2O)n, reacting with D2O the main H/D exchange mechanism is found to be proton catalyzed. In addition the H/D scrambling becomes close to statistically randomized for the larger clusters. For NH4(+)(H2O)n clusters reacting with D2O, the main mechanism is a D2O/H2O swap reaction. The lifetimes of H3O(+)(H2O)n clusters have been estimated using RRKM theory and a plateau in lifetime vs. cluster size is found already at n = 10.     

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.     

Microsolvation of hydrogen sulfide: exploration of H2S.(H2O)n and SH-.H3O+.(H2O)n-1 (n = 5-7) cluster structures on ab initio potential energy surfaces by the scaled hypersphere search method

The scaled hypersphere search method was applied to ab initio potential energy surfaces of the H2S.(H2O)n/SH-.H3O+.(H2O)n-1 system with n = 5-7. Local minima databases including 121, 326, and 553 structures for n = 5-7, respectively, were obtained based on calculations at the MP2/6-311++G(3df,2p)//B3LYP/6-31+G** level. In these small cluster sizes, the SH-.H3O+.(H2O)n-1 type is still unstable relative to the H2S.(H2O)n type, and the global minima for H2S.(H2O)n are very similar to those of pure water clusters of (H2O)n+1. Thermodynamic simulations based on the present databases showed a structure transition from the well-mixed (H2O)n+1-like global minimum at low temperatures to unmixed complexes between H2S and (H2O)n at high temperatures.     

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.     

Infrared photodissociation spectroscopy of Al(+)(CH(3)OH)(n) (n = 1-4)

Infrared photodissociation spectra of Al(+)(CH(3)OH)(n) (n = 1-4) and Al(+)(CH(3)OH)(n)-Ar (n = 1-3) were measured in the OH stretching region, 3000-3800 cm(-1). For n = 1 and 2, sharp absorption bands were observed in the free OH stretching region, all of which were well reproduced by the spectra calculated for the solvated-type geometry with no hydrogen bond. For n = 3 and 4, there were broad vibrational bands in the energy region of hydrogen-bonded OH stretching vibrations, 3000-3500 cm(-1). Energies of possible isomers for the Al(+)(CH(3)OH)(3),4 ions with hydrogen bonds were calculated in order to assign these bands. It was found that the third and fourth methanol molecules form hydrogen bonds with methanol molecules in the first solvation shell, rather than a direct bonding with the Al(+) ion. For the Al(+)(CH(3)OH)(n) clusters with n = 1-4, we obtained no evidence of the insertion reaction, which occurs in Al(+)(H(2)O)(n). One possible explanation of the difference between these two systems is that the potential energy barriers between the solvated and inserted isomers in the Al(+)(CH(3)OH)(n) system is too high to form the inserted-type isomers.     

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.     

Ultraviolet and infrared photodissociation of Si(+)(C6H6)n and Si(+)(C6H6)(n)Ar clusters

Ion-molecule complexes of the form Si(+)(C6H6)n and Si(+)(C6H6)(n)Ar are produced by laser vaporization in a pulsed nozzle cluster source. These clusters are mass-selected and studied with ultraviolet (355 nm) photodissociation and resonance-enhanced infrared photodissociation spectroscopy in the C-H stretch region of benzene. In the UV, Si(+)(C6H6)n clusters (n = 1-5) fragment to produce the Si(+)(C6H6)n mono-ligand species, suggesting that this ion has enhanced relative stability. IR photodissociation of Si(+)(C6H6)n complexes occurs by the elimination of benzene, while Si(+)(C6H6)(n)Ar complexes lose Ar. Resonances reveal C-H vibrational bands in the 2900-3300 cm(-1) region characteristic of the benzene ligand with shifts caused by the silicon cation bonding. The IR spectra confirm that the major component of the Si(+)(C6H6)n ions studied have the pi-complex structure rather than the isomeric insertion products suggested previously.     

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.     

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.     

Automated exploration of stable isomers of H(+)(H(2)O)(n) (n = 5-7) via ab initio calculations: An application of the anharmonic downward distortion following algorithm

Isomers of protonated water clusters H(+)(H(2)O)(n) (n = 5-7) have been explored on ab initio potential energy surfaces by means of the anharmonic downward distortion following algorithm. Totally, 9, 24, and 131 isomers for n = 5, 6, and 7, respectively, were located by the automatic exploration, and all of known important isomers previously reported by conventional geometry optimization approaches have been included in the present results. Moreover, structure transitions depending on n and temperature, which were observed by experimental studies, could be reproduced via thermodynamic simulation on the basis of the superposition approach and the present isomer sets.     

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.     

Periodic density functional theory study of methane activation over La(2)O(3): activity of O(2-), O(-), O(2)(2-), oxygen point defect,and Sr(2+)-doped surface sites

Results of gradient-corrected periodic density functional theory calculations are reported for hydrogen abstraction from methane at O(s)(2-), O(s)(-), O(2)(s)(2-) point defect, and Sr(2+)-doped surface sites on La(2)O(3)(001). The results show that the anionic O(s)(-) species is the most active surface oxygen site. The overall reaction energy to activate methane at an O(s)(-) site to form a surface hydroxyl group and gas-phase (*)CH(3) radical is 8.2 kcal/mol, with an activation barrier of 10.1 kcal/mol. The binding energy of hydrogen at an site O(s)(-) is -102 kcal/mol. An oxygen site with similar activity can be generated by doping strontium into the oxide by a direct Sr(2+)/La(3+) exchange at the surface. The O(-)-like nature of the surface site is reflected in a calculated hydrogen binding energy of -109.7 kcal/mol. Calculations indicate that surface peroxide (O(2(s))(2-)) sites can be generated by adsorption of O(2) at surface oxygen vacancies, as well as by dissociative adsorption of O(2) across the closed-shell oxide surface of La(2)O(3)(001). The overall reaction energy and apparent activation barrier for the latter pathway are calculated to be only 12.1 and 33.0 kcal/mol, respectively. Irrespective of the route to peroxide formation, the O(2)(s)(2-) intermediate is characterized by a bent orientation with respect to the surface and an O-O bond length of 1.47 A; both attributes are consistent with structural features characteristic of classical peroxides. We found surface peroxide sites to be slightly less favorable for H-abstraction from methane than the O(s)(-) species, with DeltaE(rxn)(CH(4)) = 39.3 kcal/mol, E(act) = 47.3 kcal/mol, and DeltaE(ads)(H) = -71.5 kcal/mol. A possible mechanism for oxidative coupling of methane over La(2)O(3)(001) involving surface peroxides as the active oxygen source is suggested.     

The importance of NO+ (H2O)4 in the conversion of NO+ (H2O)n to H3O+ (H2O)n: I. Kinetics measurements and statistical rate modeling

The kinetics for conversion of NO(+)(H(2)O)(n) to H(3)O(+)(H(2)O)(n) has been investigated as a function of temperature from 150 to 400 K. In contrast to previous studies, which show that the conversion goes completely through a reaction of NO(+)(H(2)O)(3), the present results show that NO(+)(H(2)O)(4) plays an increasing role in the conversion as the temperature is lowered. Rate constants are derived for the clustering of H(2)O to NO(+)(H(2)O)(1-3) and the reactions of NO(+)(H(2)O)(3,4) with H(2)O to form H(3)O(+)(H(2)O)(2,3), respectively. In addition, thermal dissociation of NO(+)(H(2)O)(4) to lose HNO(2) was also found to be important. The rate constants for the clustering increase substantially with the lowering of the temperature. Flux calculations show that NO(+)(H(2)O)(4) accounts for over 99% of the conversion at 150 K and even 20% at 300 K, although it is too small to be detectable. The experimental data are complimented by modeling of the falloff curves for the clustering reactions. The modeling shows that, for many of the conditions, the data correspond to the falloff regime of third body association.     

Femtosecond dynamics of Cu(H(2)O)(2)

The ultrafast relaxation dynamics of Cu(H(2)O)(2) is investigated using femtosecond photodetachment-photoionization spectroscopy. In addition, stationary points on the Cu(H(2)O)(2) anion, neutral, and cation potential energy surfaces are characterized by ab initio electronic structure calculations. Electron photodetachment from Cu(-)(H(2)O)(2) initiates the dynamics on the ground-state potential energy surface of neutral Cu(H(2)O)(2). The resulting Cu(H(2)O)(2) complexes experience large-amplitude H(2)O reorientation and dissociation. The time evolution of the Cu(H(2)O)(2) fragmentation products is monitored by time-resolved resonant multiphoton ionization. The parent ion, Cu(+)(H(2)O)(2), is not detected above background levels. The rise to a maximum of the Cu(+) signal from Cu(-)(H(2)O)(2), and the decay of the Cu(+)(H(2)O) signal from Cu(-)(H(2)O)(2) have similar tau approximately 10 ps time dependences to the corresponding signals from Cu(-)(H(2)O), but display clear differences at very short and long times. The experimental observations can be understood in terms of the following picture. Prompt dissociation of H(2)O from nascent Cu(H(2)O)(2) gives rise to a vibrationally excited Cu(H(2)O) complex, which dissociates to Cu+H(2)O due to coupling of H(2)O internal rotation to the dissociation coordinate. This prompt dissociation removes all intra-H(2)O vibrational excitation from the intermediate Cu(H(2)O) fragment, which quenches the long time vibrational predissociation to Cu+H(2)O previously observed in analogous experiments on Cu(-)(H(2)O).     

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.