Clusters of hydrated methane sulfonic acid CH3SO3H.(H2O)n (n = 1-5): a theoretical study

Ab initio and density functional methods have been used to examine the structures and energetics of the hydrated clusters of methane sulfonic acid (MSA), CH3SO3H.(H2O)n (n = 1-5). For small clusters with one or two water molecules, the most stable clusters have strong cyclic hydrogen bonds between the proton of OH group in MSA and the water molecules. With three or more water molecules, the proton transfer from MSA to water becomes possible, forming ion-pair structures between CH3SO3- and H3O+ moieties. For MSA.(H2O)3, the energy difference between the most stable ion pair and neutral structures are less than 1 kJ/mol, thus coexistence of neutral and ion-pair isomers are expected. For larger clusters with four and five water molecules, the ion-pair isomers are more stable (>10 kJ/mol) than the neutral ones; thus, proton transfer takes place. The ion-pair clusters can have direct hydrogen bond between CH3SO3- and H3O+ or indirect one through water molecule. For MSA.(H2O)5, the energy difference between ion pairs with direct and indirect hydrogen bonds are less than 1 kJ/mol; namely, the charge separation and acid ionization is energetically possible. The calculated IR spectra of stable isomers of MSA.(H2O)n clusters clearly demonstrate the significant red shift of OH stretching of MSA and hydrogen-bonded OH stretching of water molecules as the size of cluster increases.     

Microsolvation of the sodium and iodide ions and their ion pair in acetonitrile clusters: a theoretical study

The structural and thermodynamic properties of Na+(CH3CN)n, I-(CH3CN)n, and NaI(CH3CN)n clusters have been investigated by means of room-temperature Monte Carlo simulations with model potentials developed to reproduce the properties of small clusters predicted by quantum chemistry. Ions are found to adopt an interior solvation shell structure, with a first solvation shell containing approximately 6 and approximately 8 acetonitrile molecules for large Na+(CH3CN)n and I-(CH3CN)n clusters, respectively. Structural features of Na+(CH3CN)n are found to be similar to those of Na+(H2O)n clusters, but those of I-(CH3CN)n contrast with those of I-(H2O)n, for which "surface" solvation structures were observed. The potential of mean force calculations demonstrates that the NaI ion pair is thermodynamically stable with respect to ground-state ionic dissociation in acetonitrile clusters. The properties of NaI(CH3CN)n clusters exhibit some similarities with NaI(H2O)n clusters, with the existence of contact ion pair and solvent-separated ion pair structures, but, in contrast to water clusters, both types of ion pairs adopt a well-defined interior ionic solvation shell structure in acetonitrile clusters. Whereas contact ion pair species are thermodynamically favored in small clusters, solvent-separated ion pairs tend to become thermodynamically more stable above a cluster size of approximately 26. Hence, ground-state charge separation appears to occur at larger cluster sizes for acetonitrile clusters than for water clusters. We propose that the lack of a large Na+(CH3CN)n product signal in NaI(CH3CN)n multiphoton ionization experiments could arise from extensive stabilization of the ground ionic state by the solvent and possible inhibition of the photoexcitation mechanism, which may be less pronounced for NaI(H2O)n clusters because of surface solvation structures. Alternatively, increased solvent evaporation resulting from larger excess energies upon photoexcitation or major solvent reorganization on the ionized state could account for the observed solvent-selectivity in NaI cluster multiphoton ionization.     

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.     

用密度泛函方法研究质子化苯基丙酮-水团簇

The structure and growth trend of the protonated acetophenone-water clusters have been investigated using the DFT-B3LYP method combined with the standard 6-31+G(d,p) basis set. In order to obtain more accurate single-point energy the B3LYP/6-311++G(3df,2p) method was adapted. The results show that the formation of H+C8H8O-H2O is a barrierless reaction process and the equilibrium distance between the proton and the O atom in C8H8O molecule is 1.015 A. For H+C8H8O-(H2O)n(n=1,2,3) clusters, the proton lies between the acetophenone molecule C8H8O and the water molecule H2O. The distance between the proton and the O atom of the C8H8O molecule increased from n=1 to n=3; C8H8O-H+-H2O can be regarded as an solvation shell. For H+C8H8O (H2O)n (n=4,5,6,7,8) clusters, the proton lies between the two H2O molecules forming a H5O2+ structure, C8H8O-H5O2+ is an important structure, which the other H2O molecules will attack from different sides.     

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.     

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.     

Reactions of large water cluster anions with hydrogen chloride: formation of atomic hydrogen and phase separation in the gas phase

The reactions of water cluster anions (H2O)n-, n = 30-70, with hydrogen chloride have been studied by Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry. The first HCl taken up by the clusters is presumably ionically dissolved. The solvated electron recombines with the proton, which is thereby reduced to atomic hydrogen and evaporates from the cluster. This process is accompanied by blackbody radiation and collision induced loss of water molecules. Subsequent collisions lead to uptake of HCl and loss of H2O, yielding mixed clusters Cl-(HCl)m(H2O)n until they are saturated with HCl. Those saturated clusters lose H2O and HCl in a characteristic sequence. The final stage of the reaction, involving clusters with m = 0-4 and n = 0-6, is studied in detail with density functional theory calculations. The Cl-(HCl)4(H2O)6 cluster represents an example for supramolecular self-organization in the gas phase: it consists of a tetrahedral Cl-(HCl)4, connected on one side of the tetrahedron to a compact water hexamer.     

Micro-hydration of the MgNO3+ cation in the gas phase.

Coordination complexes of the magnesium nitrate cation with water [MgNO(3)(H(2)O)(n)](+) up to n=7 are investigated by experiment and theory. The fragmentation patterns of [MgNO(3)(H(2)O)(n)](+) clusters generated via electrospray ionization indicate a considerable change in stability between n=3 and 4. Further, ion-molecule reactions of mass-selected [MgNO(3)(H(2)O)(n)](+) cations with D(2)O reveal the occurrence of consecutive replacement of water ligands by heavy water, and in this respect the complexes with n=4 and 5 are somewhat more reactive than their smaller homologs with n=1-3 as well as the larger clusters with n=6 and 7. For the latter two ions, the theory suggests the existence of isomers, such as complexes with monodentate nitrato ligands as well as solvent-separated ion pairs with a common solvation shell. The reactions observed and the ion thermochemistry are discussed in the context of ab initio calculations, which also reveal the structures of the various hydrated cation complexes.     

Further insight into the relaxation dynamics of photoexcited I-(H2O)n clusters

First-principles molecular dynamics simulations of the excited-state dynamics of I-(H2O)3 have been performed to gain some insight into the general features of the relaxation process of photoexcited I-(H2O)n clusters. The relaxation of excited I-(H2O)3 is characterized by rapid motion of water molecules and slow recoil motion of the iodine atom. Both solvent reorganization and iodine atom motion appear to be important for interpreting the existing femtosecond photoelectron spectroscopy experimental results.     

Strong-field ionization and coulomb explosion of chlorine weakly bound to small water clusters

Clusters exhibit an enhancement in ionization rates under intense, ultrafast laser pulses compared to their molecular/atomic counterparts. Studies of ionization enhancement of weakly bound molecules to clusters have not been previously characterized and quantified. We demonstrate that weakly bound ClO to (H(2)O)(n) (n = 1-12) clusters and weakly bound HCl to (H(2)O)(n) (n = 1-12) clusters produce high atomic charge states of chlorine via Coulomb explosion. Density functional theory (DFT) was used to qualitatively compare the interaction energy of ClO with respect to the number of water molecules as well as HCl with respect to the number of water molecules. The chlorine ion signal intensity for each atomic charge state was observed to be dependent on the molecule-cluster bond strength. The observed ionization enhancement was quantified using semiclassical tunneling theory, and it was found that the Cl(3+-5+) and O(2+) charge states are enhanced in ionization. Possible mechanisms of ionization enhancement are explored for weakly bound chlorine species.     

Structures and dissociation channels of protonated mixed clusters around a small magic number: infrared spectroscopy of ((CH3)3N)n-H(+)-H2O (n = 1-3)

The magic number behavior of ((CH(3))(3)N)(n)-H(+)-H(2)O clusters at n = 3 is investigated by applying infrared spectroscopy to the clusters of n = 1-3. Structures of these clusters are determined in conjunction with density functional theory calculations. Dissociation channels upon infrared excitation are also measured, and their correlation with the cluster structures is examined. It is demonstrated that the magic number cluster has a closed-shell structure, in which the water moiety is surrounded by three (CH(3))(3)N molecules. The ion core (protonated site) of the clusters is found to be (CH(3))(3)NH(+) for n = 1-3, but coexistence of an isomer of the H(3)O(+) ion core cannot be ruled out for n = 3. Large rearrangement of the cluster structures of n = 2 and 3 before dissociation, which has been suggested in the mass spectrometric studies, is confirmed on the basis of the structure determination by infrared spectroscopy.     

Stepwise hydration of ionized aromatics. Energies, structures of the hydrated benzene cation, and the mechanism of deprotonation reactions

The stepwise binding energies (DeltaHdegree(n-1,n)) of 1-8 water molecules to benzene(.+) [Bz(.+)(H2O)n] were determined by equilibrium measurements using an ion mobility cell. The stepwise hydration energies, DeltaHdegree(n-1,n), are nearly constant at 8.5 +/- 1 kcal mol-1 from n = 1-6. Calculations show that in the n = 1-4 clusters, the benzene(.+) ion retains over 90% of the charge, and it is extremely solvated, that is, hydrogen bonded to an (H2O)n cluster. The binding energies and entropies are larger in the n = 7 and 8 clusters, suggesting cyclic or cage-like water structures. The concentration of the n = 3 cluster is always small, suggesting that deprotonation depletes this ion, consistent with the thermochemistry since associative deprotonation Bz(.+)(H2O)(n-1) + H2O-->C6H5. + (H2O)nH+ is thermoneutral or exothermic for n > or = 4. Associative intracluster proton transfer Bz(.+)(H2O)(n+1) + H2O-->C6H5.(H2O)nH+ would also be exothermic for n > or = 4, but lack of H/D exchange with D2O shows that the proton remains on C6H6(.+) in the observed Bz(.+)(H2O)n clusters. This suggests a barrier to intracluster proton transfer, and as a result, the [Bz(.+)(H2O)n]* activated complexes either undergo dissociative proton transfer, resulting in deprotonation and generation of (H2O)nH+, or become stabilized. The rate constant for the deprotonation reaction shows a uniquely large negative temperature coefficient of K = cT(-67+/-4) (or activation energy of -34+/- 1 kcal mol-1), caused by a multibody mechanism in which five or more components need to be assembled for 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.     

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.     

Mid-infrared characterization of the NH4 +(H2O)n clusters in the neighborhood of the n=20 "magic" number

Vibrational predissociation spectra are reported for size-selected NH4+ (H2O)n clusters (n=5-22) in the 2500-3900 cm(-1) region. We concentrate on the sharp free OH stretching bands to deduce the local H-bonding configurations of water molecules on the cluster surface. As in the spectra of the protonated water clusters, the free OH bands in NH4+ (H2O)n evolve from a quartet at small sizes (n<7), to a doublet around n=9, and then to a single peak at the n=20 magic number cluster, before the doublet re-emerges at larger sizes. This spectral simplification at the magic number cluster mirrors that found earlier in the H+(H2O)n clusters. We characterize the likely structures at play for the n=19 and 20 clusters with electronic structure calculations. The most stable form of the n=20 cluster is predicted to have a surface-solvated NH4+ ion that lies considerably lower in energy than isomers with the NH4+ in the interior.     

Photodissociation of hydrogen halide molecules on free ice nanoparticles

Photodissociation of water clusters doped with HX(X=Br,Cl), molecules has been studied in a molecular beam experiment. The HX(H2O)n clusters are dissociated with 193 nm laser pulses, and the H fragments are ionized at 243.07 nm and their time-of-flight distributions are measured. Experiments with deuterated species DBr(H2O)n and HBr(D2O)n suggest that the photodissociation signal originates from the presence of the HX molecule on the water cluster, but does not come directly from a photolysis of the HX molecule. The H fragment is proposed to originate from the hydronium molecule H3O. Possible mechanisms of the H3O production are discussed. Experimental evidence suggests that acidic dissociation takes place in the cluster, but the H3O+ ion remains rather immobile.     

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.     

Photochemistry of hydrogen halides on water clusters: simulations of electronic spectra and photodynamics, and comparison with photodissociation experiments

The photochemistry of small HX·(H(2)O)(n), n = 4 and 5 and X = F, Cl, and Br, clusters has been modeled by means of ab initio-based molecular simulations. The theoretical results were utilized to support our interpretation of photodissociation experiments with hydrogen halides on ice nanoparticles HX·(H(2)O)(n), n ≈ 10(2)-10(3). We have investigated the HX·(H(2)O)(n) photochemistry for three structural types: covalently bound structures (CBS) and acidically dissociated structures in a form of contact ion pair (CIP) and solvent separated pair (SSP). For all structures, we have modeled the electronic absorption spectra using the reflection principle combined with a path integral molecular dynamics (PIMD) estimate of the ground state density. In addition, we have investigated the solvent effect of water on the absorption spectra within the nonequilibrium polarizable continuum model (PCM) scheme. The major conclusion from these calculations is that the spectra for ionic structures CIP and SSP are significantly red-shifted with respect to the spectra of CBS structures. We have also studied the photodynamics of HX·(H(2)O)(n) clusters using the Full Multiple Spawning method. In the CBS structures, the excitation led to almost immediate release of the hydrogen atom with high kinetic energy. The light absorption in ionically dissociated species generates the hydronium radical (H(3)O) and halogen radical (X) within a charge-transfer-to-solvent (CTTS) excitation process. The hydronium radical ultimately decays into a water molecule and hydrogen atom with a characteristic kinetic energy irrespective of the hydrogen halide. We have also investigated the dynamics of an isolated and water-solvated H(3)O radical that we view as a central species in water radiation chemistry. The theoretical findings support the following picture of the HX photochemistry on ice nanoparticles investigated in our molecular beam experiments: HX is acidically dissociated in the ground state on ice nanoparticles, generating the CIP structure, which is then excited by the UV laser light into the CTTS states, followed by the H(3)O radical formation.     

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

Pure neutral (CO2)n clusters and mixed (CO2)n(H2O)m clusters are investigated employing time of flight mass spectroscopy and single photon ionization at 26.5 eV. The distribution of pure (CO2)n clusters decreases roughly exponentially with increasing cluster size. During the ionization process, neutral clusters suffer little fragmentation because almost all excess cluster energy above the vertical ionization energy is taken away by the photoelectron and only a small part of the photon energy is deposited into the (CO2)n cluster. Metastable dissociation rate constants of (CO2)n+ are measured in the range of (0.2-1.5) x 10(4) s(-1) for cluster sizes of 5< or =n< or =16. Mixed CO2-H2O clusters are studied under different generation conditions (5% and 20% CO2 partial pressures and high and low expansion pressures). At high CO2 concentration, predominant signals in the mass spectrum are the (CO2)n+ cluster ions. The unprotonated cluster ion series (CO2)nH2O+ and (CO2)n(H2O)2+ are also observed under these conditions. At low CO2 concentration, protonated cluster ions (H2O)nH+ are the dominant signals, and the protonated CO2(H2O)nH+ and unprotonated (H2O)n+ and (CO2)(H2O)n+ cluster ion series are also observed. The mechanisms and dynamics of the formation of these neutral and ionic clusters are discussed.