Dissociative recombination of H+(H2O)3 and D+(D2O)3 water cluster ions with electrons: cross sections and branching ratios

Dissociative recombination (DR) of the water cluster ions H(+)(H(2)O)(3) and D(+)(D(2)O)(3) with electrons has been studied at the heavy-ion storage ring CRYRING (Manne Siegbahn Laboratory, Stockholm University). For the first time, absolute DR cross sections have been measured for H(+)(H(2)O)(3) in the energy range of 0.001-0.8 eV, and relative cross sections have been measured for D(+)(D(2)O)(3) in the energy range of 0.001-1.0 eV. The DR cross sections for H(+)(H(2)O)(3) are larger than previously observed for H(+)(H(2)O)(n) (n=1,2), which is in agreement with the previously observed trend indicating that the DR rate coefficient increases with size of the water cluster ion. Branching ratios have been determined for the dominating product channels. Dissociative recombination of H(+)(H(2)O)(3) mainly results in the formation of 3H(2)O+H (probability of 0.95+/-0.05) and with a possible minor channel resulting in 2H(2)O+OH+H(2) (0.05+/-0.05). The dominating channels for DR of D(+)(D(2)O)(3) are 3D(2)O+D (0.88+/-0.03) and 2D(2)O+OD+D(2) (0.09+/-0.02). The branching ratios are comparable to earlier DR results for H(+)(H(2)O)(2) and D(+)(D(2)O)(2), which gave 2X(2)O+X (X=H,D) with a probability of over 0.9.     

Structural rearrangements and magic numbers in reactions between pyridine-containing water clusters and ammonia

Molecular cluster ions H(+)(H(2)O)(n), H(+)(pyridine)(H(2)O)(n), H(+)(pyridine)(2)(H(2)O)(n), and H(+)(NH(3))(pyridine)(H(2)O)(n) (n = 16-27) and their reactions with ammonia have been studied experimentally using a quadrupole-time-of-flight mass spectrometer. Abundance spectra, evaporation spectra, and reaction branching ratios display magic numbers for H(+)(NH(3))(pyridine)(H(2)O)(n) and H(+)(NH(3))(pyridine)(2)(H(2)O)(n) at n = 18, 20, and 27. The reactions between H(+)(pyridine)(m)(H(2)O)(n) and ammonia all seem to involve intracluster proton transfer to ammonia, thus giving clusters of high stability as evident from the loss of several water molecules from the reacting cluster. The pattern of the observed magic numbers suggest that H(+)(NH(3))(pyridine)(H(2)O)(n) have structures consisting of a NH(4)(+)(H(2)O)(n) core with the pyridine molecule hydrogen-bonded to the surface of the core. This is consistent with the results of high-level ab initio calculations of small protonated pyridine/ammonia/water 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.     

Electrospray droplet impact/secondary ion mass spectrometry: cluster ion formation

The electrospray droplets that are sampled through an orifice into the vacuum chamber are accelerated by 10 kV and impact on the stainless steel substrate. The mass and the kinetic energy of electrospray droplets are roughly estimated to be a few 10(6) u and approximately 10(6) eV, respectively. The molecular ion M(+.) and the protonated molecule [M+H](+) are observed as secondary ions for chrysene and coronene deposited on the metal substrate (no matrix used). The ionization may take place in the shock wave generated by the high-momentum coherent collision between the droplet projectile and the solid sample. Cluster ions of H(+)(H(2)O)(n) and CF(3)COO(-)(H(2)O)(n), with n up to approximately 150, were observed as secondary ions formed by the electrospray droplet impact ionization (EDI) for 10(-2) M trifluoroacetic acid (TFA) aqueous solution. This indicates that the charged droplets that collide with the metal substrate with the kinetic energy of approximately 10(6) eV do not vaporize completely but are disintegrated into many tiny microdroplets. The ion signal intensity anomalies (i.e. magic numbers) were observed for the cluster ions of H(3)O(+)(H(2)O)(n) and CF(3)COO(-)(H(2)O)(n) for 10(-2) M TFA aqueous solution and of Cs(+)(H(2)O)(n), I(-)(H(2)O)(n), Cs(+)(CsI)(n), and I(-)(CsI)(n) for 10(-2) M CsI aqueous solution.     

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.     

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.     

Dissociative recombination of the deuterated acetaldehyde ion, CD3CDO(+): product branching fractions, absolute cross sections and thermal rate coefficient

Dissociative recombination of the deuterated acetaldehyde ion CD3CDO(+) has been studied at the heavy-ion storage ring CRYRING, located at the Manne Siegbahn Laboratory, Stockholm, Sweden. Product branching fractions together with absolute DR cross-sections were measured. The branching fractions were determined at a relative collision energy between the ions and the electrons of approximately 0 eV. With a probability of 34% the DR events resulted in no ruptures of bonds between heavy atoms (i.e. no breakage of the C-C bond or the C[double bond, length as m-dash]O bond). In the remaining 66% of the events one of the bonds between the heavy atoms was broken. The energy-dependent cross-section for the DR reaction was measured between approximately 0 and 1 eV relative kinetic energy. In the energy region between 1 meV and 0.2 eV the absolute cross section could be fitted by the expression sigma(E) = 6.8 x 10(-16)E(-1.28) cm(2), whereas in the energy interval between 0.2 and 1 eV the data were best fitted by sigma(E) = 4.1 x 10(-16)E(-1.60) cm(2). From these cross section data the thermal rate coefficient (as a function of the electron temperature), alpha(T) = 9.2 x 10(-7) (T/300)(-0.72) cm(3) s(-1) was obtained.     

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.     

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.     

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.     

Effect of formic acid addition on water cluster stability and structure

Computational chemistry simulations were performed to determine the effect that the addition of a single formic acid molecule has on the structure and stability of protonated water clusters. Previous experimental studies showed that addition of formic acid to protonated pure water results in higher intensities of large-sized clusters when compared to pure water and methanol-water mixed clusters. For larger, protonated clusters, molecular dynamics simulations were performed on H(+)(H(2)O)(n), H(+)(H(2)O)(n)CH(3)OH, and H(+)(H(2)O)(n)CHOOH clusters, 19-28 molecules in size, using a reactive force field (ReaxFF). Based on these computations, formic acid-water clusters were found to have significantly higher binding energies per molecule. Addition of formic acid to a water cluster was found to alter the structure of the hydrogen-bonding network, creating selective sites within the cluster, enabling the formation of new hydrogen bonds, and increasing both the stability of the cluster and its rate of growth.     

Alternative low-energy mechanisms for isotopic exchange in gas-phase D2O-H+(H2O)n reactions.

Molecular-dynamics (MD) trajectories and high-level ab initio methods have been used to study the low-energy mechanism for D(2)O-H(+)(H(2)O)(n) reactions. At low collisional energies, MD simulations show that the collisional complexes are long-lived and undergo fast monomolecular isomerization, converting between different isomers within 50-500 ps. Such processes, primarily involving water-molecule shifts along a water chain, require the surmounting of very-low-energy barriers and present sizable non- Rice-Ramsperger-Kassel-Marcus (RRKM) effects, which are interpreted as a lack of randomization of the internal kinetic energy. Interestingly, the rate of water shifts was found to increase upon increasing the size of the cluster. Based on these findings, we propose to incorporate the following steps into the mechanism for low-energy isotopic scrambling these D(2)O-H(+)(H(2)O)(n) reactions: a) formation of the collisional complex [H(+)(H(2)O)(n)D(2)O]* in a vibro-rotational excited state; b) incorporation of the heavy-water molecule in the cluster core as HD(2)O(+) by means of isomerization involving molecular shifts; c) displacement of solvation molecules from the first shell of HD(2)O(+) inducing de-deuteration (shift of a D(+) to a neighbor water molecule); d) reorganization of the clusters and/or expulsion of one of the isotopic variants of water (H(2)O, HDO or D(2)O) from the periphery of the complex.     

Unexpected size distribution of Ba(H2O)n clusters: why is the intensity of the Ba(H2O)1 cluster anomalously low?

An experimental and theoretical study on the reactivity of neutral Ba atoms with water clusters has been conducted to unravel the origin of the irregular intensity pattern observed in one-photon ionization mass spectra of a Ba(H(2)O)(n)/BaOH(H(2)O)(n-1) (n = 1-4) cluster distribution, which was generated in a laser vaporization-supersonic expansion source. The most remarkable irregular feature is the finding for n = 1 of a lower intensity for the Ba(+)(H(2)O)(n) peak with respect to that of BaOH(+)(H(2)O)(n-1), which is opposite to the trend for n = 2-4. Rationalization of the data required consideration of a distinct behavior of ground-state and electronically excited state Ba atoms in inelastic and reactive Ba + (H(2)O)(n) encounters that can occur in the cluster source. Within this picture, the generation of Ba(H(2)O)(n) (n > 1) association products results from stabilizing collisions with atoms of the carrier gas, which are favored by intramolecular vibrational redistribution that operates on the corresponding collision intermediates prior to stabilization; the latter is unlikely to occur for Ba + (H(2)O) encounters. Overall, this interpretation is consistent with additional in-source laser excitation and quenching experiments, which aimed to explore qualitatively the effect of perturbing the Ba atom electronic state population distribution on the observed intensity pattern, as well as with the energetics of various possible reactions for the Ba + H(2)O system that derive from high level ab initio calculations.     

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).     

Thermochemical stabilities and structures of the cluster ions OCS+, S2+, H+(OCS), and C2H5+ with OCS molecules in the gas phase

The gas-phase clustering reactions of OCS+, S2+, H+(OCS), and C2H5+ ions with carbonyl sulfide (OCS) molecules were studied using a pulsed electron-beam high-pressure mass spectrometer and applying density functional theory (DFT) calculations. In the cluster ions OCS+(OCS)(n) and H+(OCS)(OCS)(n), a moderately strong, here referred to as "semi-covalent", bond was formed with n = 1. However, the nature of bonding changed from semi-covalent to electrostatic with n = 1 --> 2. The bond energy of S2(+)(OCS) was determined experimentally to be 12.9 +/- 1 kcal/mol, which is significantly smaller than that of the isovalent S2(+)(CS2) complex (30.9 +/- 1.5 kcal/mol). DFT based calculations predicted the presence of several isomeric structures for H+(OCS)(OCS)(n) complexes. The bond energies in the C2H5+(OCS)(n) clusters showed an irregular decrease for n = 1 --> 2 and 7 --> 8. The nonclassical bridge structure for the free C2H5+ isomerized to form a semi-covalent bond with one OCS ligand, [H3CCH2...SCO]+, i.e., reverted to classical structure. However, the nonclassical bridge structure of C2H5+ was preserved in the cluster ions C2H5+(OCS)(n) below 140 K attributable to the lack of thermal energy for the isomerization. DFT calculations revealed that stability orders of the geometric isomers of H+(OCS)(OCS)(n) and C2H5+(OCS)(n) changed with increasing n values.     

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.     

Sequential bond energies and barrier heights for the water loss and charge separation dissociation pathways of Cd(2+)(H2O)n, n = 3-11

The bond dissociation energies for losing one water from Cd(2+)(H(2)O)(n) complexes, n = 3-11, are measured using threshold collision-induced dissociation in a guided ion beam tandem mass spectrometer coupled with a thermal electrospray ionization source. Kinetic energy dependent cross sections are obtained for n = 4-11 complexes and analyzed to yield 0 K threshold measurements for loss of one, two, and three water ligands after accounting for multiple collisions, kinetic shifts, and energy distributions. The threshold measurements are converted from 0 to 298 K values to give the hydration enthalpies and free energies for sequentially losing one water from each complex. Theoretical geometry optimizations and single point energy calculations are performed on reactant and product complexes using several levels of theory and basis sets to obtain thermochemistry for comparison to experiment. The charge separation process, Cd(2+)(H(2)O)(n) → CdOH(+)(H(2)O)(m) + H(+)(H(2)O)(n-m-1), is also observed for n = 4 and 5 and the competition between this process and water loss is analyzed. Rate-limiting transition states for the charge separation process at n = 3-6 are calculated and compared to experimental threshold measurements resulting in the conclusion that the critical size for this dissociation pathway of hydrated cadmium is n(crit) = 4.     

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.     

Structures, energetics, vibrational spectra of NH4+ (H2O)(n=4,6) clusters: Ab initio calculations and first principles molecular dynamics simulations

Important structural isomers of NH(4) (+)(H(2)O)(n=4,6) have been studied by using density functional theory, Moller-Plesset second order perturbation theory, and coupled-cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. The zero-point energy (ZPE) correction to the complete basis set limit of the CCSD(T) binding energies and free energies is necessary to identify the low energy structures for NH(4) (+)(H(2)O)(n=4,6) because otherwise wrong structures could be assigned for the most probable structures. For NH(4) (+)(H(2)O)(6), the cage-type structure, which is more stable than the previously reported open structure before the ZPE correction, turns out to be less stable after the ZPE correction. In first principles Car-Parrinello molecular dynamics simulations around 100 K, the combined power spectrum of three lowest energy isomers of NH(4) (+)(H(2)O)(4) and two lowest energy isomers of NH(4) (+)(H(2)O)(6) explains each experimental IR spectrum.     

Structure and hydration of the C4H4●+ ion formed by electron impact ionization of acetylene clusters

Here we report ion mobility experiments and theoretical studies aimed at elucidating the identity of the acetylene dimer cation and its hydrated structures. The mobility measurement indicates the presence of more than one isomer for the C(4)H(4)(●+) ion in the cluster beam. The measured average collision cross section of the C(4)H(4)(●+) isomers in helium (38.9 ± 1 A?(2)) is consistent with the calculated cross sections of the four most stable covalent structures calculated for the C(4)H(4)(●+) ion [methylenecyclopropene (39.9 A?(2)), 1,2,3-butatriene (41.1 A?(2)), cyclobutadiene (38.6 A?(2)), and vinyl acetylene (41.1 A?(2))]. However, none of the single isomers is able to reproduce the experimental arrival time distribution of the C(4)H(4)(●+) ion. Combinations of cyclobutadiene and vinyl acetylene isomers show excellent agreement with the experimental mobility profile and the measured collision cross section. The fragment ions obtained by the dissociation of the C(4)H(4)(●+) ion are consistent with the cyclobutadiene structure in agreement with the vibrational predissociation spectrum of the acetylene dimer cation (C(2)H(2))(2)(●+) [R. A. Relph, J. C. Bopp, J. R. Roscioli, and M. A. Johnson, J. Chem. Phys. 131, 114305 (2009)]. The stepwise hydration experiments show that dissociative proton transfer reactions occur within the C(4)H(4)(●+)(H(2)O)(n) clusters with n ≥ 3 resulting in the formation of protonated water clusters. The measured binding energy of the C(4)H(4)(●+)H(2)O cluster, 38.7 ± 4 kJ/mol, is in excellent agreement with the G3(MP2) calculated binding energy of cyclobutadiene(●+)·H(2)O cluster (41 kJ/mol). The binding energies of the C(4)H(4)(●+)(H(2)O)(n) clusters change little from n = 1 to 5 (39-48 kJ/mol) suggesting the presence of multiple binding sites with comparable energies for the water-C(4)H(4)(●+) and water-water interactions. A significant entropy loss is measured for the addition of the fifth water molecule suggesting a structure with restrained water molecules, probably a cyclic water pentamer within the C(4)H(4)(●+)(H(2)O)(5) cluster. Consequently, a drop in the binding energy of the sixth water molecule is observed suggesting a structure in which the sixth water molecule interacts weakly with the C(4)H(4)(●+)(H(2)O)(5) cluster presumably consisting of a cyclobutadiene(●+) cation hydrogen bonded to a cyclic water pentamer. The combination of ion mobility, dissociation, and hydration experiments in conjunction with the theoretical calculations provides strong evidence that the (C(2)H(2))(2)(●+) ions are predominantly present as the cyclobutadiene cation with some contribution from the vinyl acetylene cation.