Structures and infrared spectra of fluoride-hydrogen sulfide clusters from ab initio calculations: F(-)-(H2S)n, n = 1-5

Clusters formed between a fluoride anion and several hydrogen sulfide molecules have been investigated via ab initio calculations at the MP2 level of theory, using Dunning's augmented correlation consistent basis sets. Optimised geometries, vibrational frequencies, and enthalpy changes for the ligand association reactions are presented for clusters with up to five H2S ligands interacting with a F- anion. The minimum energy structure for the 1:1 F(-)-H2S complex features proton transfer from the H2S to the F- anion, forming a planar C(s) symmetry FH...SH- structure. For the F(-)-(H2S)2 cluster, the FH...SH- core remains and is solvated by a perturbed H2S ligand. For the larger F(-)-(H2S)(3-5) clusters, in addition to the FH...SH(-)-(H2S)n cluster forms, other minima featuring a 'solvated F-' anion are predicted. Calculated infrared spectra for the minima of each cluster size are presented to aid in assigning spectra from future experimental studies.     

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.     

Computational study on the characteristics of the interaction in naphthalene...(H2X)n=1,2 (X = O,S) clusters

The characteristics of the interaction between the pi cloud of naphthalene and up to two H2O or H2S molecules were studied. Calculations show that clusters formed by naphthalene and one H2O or H2S molecule have similar geometric features, and also present similar interaction energies. Our best estimates for the interaction energy amount to -2.95 and -2.92 kcal/mol for H2O and H2S, respectively, as obtained with the CCSD(T) method. Calculations at the MP2 level employing large basis sets should be avoided because they produce highly overestimated interaction energies, especially for hydrogen sulfide complexes. The MPWB1K functional, however, provides values pretty similar to those obtained with the CCSD(T) method. Although the magnitude of the interaction is similar with both H2X molecules, its nature is somewhat different: the H2O complex presents electrostatic and dispersion contributions of similar magnitude, whereas for H2S the interaction is dominated by dispersion. In clusters containing two H2X molecules several minima were characterized. In water clusters, the total interaction energy is dominated by the presence of a O-H...O hydrogen bond and, as a consequence, structures where this contact is present are the most stable. However, clusters containing H2S show structures with no interaction between H2S moieties which are as stable as the hydrogen bonded ones, because they allow an optimal H2S...naphthalene interaction, which is stronger than the S-H...S contact. Therefore it is possible that in larger polycycles hydrogen sulfide molecules will be spread onto the surface maximizing S-H...pi interactions rather than aggregated, forming H2S clusters.     

Low-temperature adsorption of H2S on Ag(111)

H(2)S forms a rich variety of structures on Ag(111) at low temperature and submonolayer coverage. The molecules decorate step edges, exist as isolated entities on terraces, and aggregate into clusters and islands, under various conditions. One type of island exhibits a (√37×√37)R25.3° unit cell. Typically, molecules in the clusters and islands are separated by about 0.4 nm, the same as the S-S separation in crystalline H(2)S. Density functional theory indicates that hydrogen-bonded clusters contain two types of molecules. One is very similar to an isolated adsorbed H(2)S molecule, with both S-H bonds nearly parallel to the surface. The other has a S-H bond pointed toward the surface. The potential energy surface for adsorption and diffusion is very smooth.     

Theoretical studies on photoelectron and IR spectral properties of Br2.-(H2O)n clusters

We report vertical detachment energy (VDE) and IR spectra of Br2.-.(H2O)n clusters (n=1-8) based on first principles electronic structure calculations. Cluster structures and IR spectra are calculated at Becke's half-and-half hybrid exchange-correlation functional (BHHLYP) with a triple split valence basis function, 6-311++G(d,p). VDE for the hydrated clusters is calculated based on second order Moller-Plesset perturbation (MP2) theory with the same set of basis function. On full geometry optimization, it is observed that conformers having interwater hydrogen bonding among solvent water molecules are more stable than the structures having double or single hydrogen bonded structures between the anionic solute, Br2.-, and solvent water molecules. Moreover, a conformer having cyclic interwater hydrogen bonded network is predicted to be more stable for each size hydrated cluster. It is also noticed that up to four solvent H2O units can reside around the solute in a cyclic interwater hydrogen bonded network. The excess electron in these hydrated clusters is localized over the solute atoms. Weighted average VDE is calculated for each size (n) cluster based on statistical population of the conformers at 150 K. A linear relationship is obtained for VDE versus (n+3)(-1/3) and bulk VDE of Br2.- aqueous solution is calculated as 10.01 eV at MP2 level of theory. BHHLYP density functional is seen to make a systematic overestimation in VDE values by approximately 0.5 eV compared to MP2 data in all the hydrated clusters. It is observed that hydration increases VDE of bromine dimer anion system by approximately 6.4 eV. Calculated IR spectra show that the formation of Br2.--water clusters induces large shifts from the normal O-H stretching bands of isolated water keeping bending modes rather insensitive. Hydrated clusters, Br2.-.(H2O)n, show characteristic sharp features of O-H stretching bands of water in the small size clusters.     

Coordination of methanol clusters to benzene: a computational study

Benzene-methanol cluster structures were investigated with theoretical chemistry methods to describe the microsolvation of benzene and the benzene-methanol azeotrope. Benzene-methanol (MeOH) clusters containing up to six methanol molecules have been calculated by ab initio [MP2/6-311++G(d,p)//MP2/6-31+G(d,p) + BSSE correction] method. The BSSE was found quite large with this basis set, hence, different extrapolation schemes in combination with the aug-cc-pVxZ basis sets have been used to estimate the complete basis set limit of the MP2 interaction energy [ΔE(MP2/CBS)]. For smaller clusters, n ≤ 3, DFT procedures (DFTB+, MPWB1K, M06-2X) have also been applied. Geometries obtained for these clusters by M06-2X and MP2 calculations are quite similar. Based on the MP2/CBS results, the most stable C(6)H(6)(MeOH)(3) cluster is characterized by a hydrogen bonded MeOH trimer chain interacting with benzene via π···H-O and O···H-C(benzene) hydrogen bonds. Larger benzene-MeOH clusters with n ≥ 4 consist of cyclic (MeOH)(n) subclusters interacting with benzene by dispersive forces, to be denoted by C(6)H(6) + (MeOH)(n). Interaction energies and cooperativity effects are discussed in comparison with methanol clusters. Besides MP2/CBS calculations, for selected larger clusters the M06-2X/6-311++G(d,p)//M06-2X/6-31+G(d,p) procedure including the BSSE correction was also used. Interaction energies obtained thereby are usually close to the MP2/CBS limit. To model the benzene-MeOH azeotrope, several structures for (C(6)H(6))(2)(MeOH)(3) clusters have been calculated. The most stable structures contain a tilted T-shaped benzene dimer interacting by π···H-O and O···H-C (benzene) hydrogen bonds with a (MeOH)(3) chain. A slightly less negative interaction energy results for a parallel displaced benzene sandwich dimer with a (MeOH)(3) chain atop of one of the benzene molecules.     

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.     

Infrared spectra of the Li+-(H2)n (n=1-3) cation complexes

The Li+-(H2)n n=1-3 complexes are investigated through infrared spectra recorded in the H-H stretch region (3980-4120 cm-1) and through ab initio calculations at the MP2/aug-cc-pVQZ level. The rotationally resolved H-H stretch band of Li+-H2 is centered at 4053.4 cm-1 [a -108 cm-1 shift from the Q1(0) transition of H2]. The spectrum exhibits rotational substructure consistent with the complex possessing a T-shaped equilibrium geometry, with the Li+ ion attached to a slightly perturbed H2 molecule. Around 100 rovibrational transitions belonging to parallel Ka=0-0, 1-1, 2-2, and 3-3 subbands are observed. The Ka=0-0 and 1-1 transitions are fitted by a Watson A-reduced Hamiltonian yielding effective molecular parameters. The vibrationally averaged intermolecular separation in the ground vibrational state is estimated as 2.056 A increasing by 0.004 A when the H2 subunit is vibrationally excited. The spectroscopic data are compared to results from rovibrational calculations using recent three dimensional Li+-H2 potential energy surfaces [Martinazzo et al., J. Chem. Phys. 119, 11241 (2003); Kraemer and Spirko, Chem. Phys. 330, 190 (2006)]. The H-H stretch band of Li+-(H2)2, which is centered at 4055.5 cm-1 also exhibits resolved rovibrational structure. The spectroscopic data along with ab initio calculations support a H2-Li+-H2 geometry, in which the two H2 molecules are disposed on opposite sides of the central Li+ ion. The two equivalent Li+...H2 bonds have approximately the same length as the intermolecular bond in Li+-H2. The Li+-(H2)3 cluster is predicted to possess a trigonal structure in which a central Li+ ion is surrounded by three equivalent H2 molecules. Its infrared spectrum features a broad unresolved band centered at 4060 cm-1.     

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.     

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.     

An ab initio investigation on (CO2)n and CO2(Ar)m clusters: geometries and IR spectra

An ab initio investigation on CO(2) homoclusters is done at MPWB1K6-31++G(2d) level of theory. Electrostatic guidelines are found to be useful for generating initial structures of (CO(2))(n) clusters. The ab initio minimum energy geometries of (CO(2))(n) with n=2-8 are T shaped, cyclic, trigonal pyramidal, tetragonal pyramidal, tetragonal bipyramidal, pentagonal bipyramidal, and pentagonal bipyramid with one CO(2) molecule attached to it. A test calculation on (CO(2))(20) cluster is also reported. The geometric parameters of the energetically most favored (CO(2))(n) clusters match quite well their experimental counterparts (wherever available) as well as those derived from molecular dynamics studies. The effect of clustering is quantified through the asymmetric C-O stretching frequency shift relative to the single CO(2) molecule. (CO(2))(n) clusters show an increasing blueshift from 1.8 to 9.6 cm(-1) on increasing number of CO(2) molecules from n=2 to 8. The energetics and geometries of CO(2)(Ar)(m) clusters have also been explored at the same level of theory. The geometries for m=1-6 show a predominant T type of the argon-CO(2) molecule interaction. Higher clusters with m=7-12 show that the argon atoms cluster around the oxygen atom after the saturation of the central carbon atom. The CO(2)(Ar)(m) clusters exhibit an increasing redshift in the C-O asymmetric stretch relative to CO(2) molecule of 0.7-5.6 cm(-1) with increasing number of argon atoms through m=1-8.     

Development of a 3-body:many-body integrated fragmentation method for weakly bound clusters and application to water clusters (H2O)(n = 3-10, 16, 17)

A 3-body:many-body integrated quantum mechanical (QM) fragmentation method for non-covalent clusters is introduced within the ONIOM formalism. The technique captures all 1-, 2-, and 3-body interactions with a high-level electronic structure method, while a less demanding low-level method is employed to recover 4-body and higher-order interactions. When systematically applied to 40 low-lying (H(2)O)(n) isomers ranging in size from n = 3 to 10, the CCSD(T):MP2 3-body:many-body fragmentation scheme deviates from the full CCSD(T) interaction energy by no more than 0.07 kcal mol(-1) (or <0.01 kcal mol(-1) per water). The errors for this QM:QM method increase only slightly for various low-lying isomers of (H(2)O)(16) and (H(2)O)(17) (always within 0.13 kcal mol(-1) of the recently reported canonical CCSD(T)/aug-cc-pVTZ energies). The 3-body:many-body CCSD(T):MP2 procedure is also very efficient because the CCSD(T) computations only need to be performed on subsets of the cluster containing 1, 2, or 3 monomers, which in the current context means the largest CCSD(T) calculations are for 3 water molecules, regardless of the cluster size.     

Treatment of dilute clusters of methanol and water by ab initio quantum mechanical calculations

Large molecular clusters can be considered as intermediate states between gas and condensed phases, and information about them can help us understand condensed phases. In this paper, ab initio quantum mechanical methods have been used to examine clusters formed of methanol and water molecules. The main goal was to obtain information about the intermolecular interactions and the structure of methanol/water clusters at the molecular level. The large clusters (CH(4)O...(H(2)O)(12) and H(2)O...(CH(4)O)(10)) containing one molecule of one component (methanol or water) and many (12, 10) molecules of the other component were considered. M?ller-Plesset perturbation theory (MP2) was used in the calculations. Several representative cluster geometries were optimized, and nearest-neighbor interaction energies were calculated for the geometries obtained in the first step. The results of the calculations were compared to the available experimental information regarding the liquid methanol/water mixtures and to the molecular dynamics and Monte Carlo simulations, and good agreement was found. For the CH(4)O...(H(2)O)(12) cluster, it was shown that the molecules of water can be subdivided into two classes: (i) H bonded to the central methanol molecule and (ii) not H bonded to the central methanol molecule. As expected, these two classes exhibited striking energy differences. Although they are located almost the same distance from the carbon atom of the central methanol molecule, they possess very different intermolecular interaction energies with the central molecule. The H bonding constitutes a dominant factor in the hydration of methanol in dilute aqueous solutions. For the H(2)O...(CH(4)O)(10) cluster, it was shown that the central molecule of water has almost three H bonds with the methanol molecules; this result differs from those in the literature that concluded that the average number of H bonds between a central water molecule and methanol molecules in dilute solutions of water in methanol is about two, with the water molecules being incorporated into the chains of methanol. In contrast, the present predictions revealed that the central water molecule is not incorporated into a chain of methanol molecules, but it can be the center of several (2-3) chains of methanol molecules. The molecules of methanol, which are not H bonded to the central water molecule, have characteristics similar to those of the methane molecules around a central water molecule in the H(2)O...(CH(4))(10) cluster. The ab initio quantum mechanical methods employed in this paper have provided detailed information about the H bonds in the clusters investigated. In particular, they provided full information about two types of H bonds between water and methanol molecules (in which the water or the methanol molecule is the proton donor), including information about their energies and lengths. The average numbers of the two types of H bonds in the CH(4)O...(H(2)O)(12) and H(2)O...(CH(4)O)(10) clusters have been calculated. Such information could hardly be obtained with the simulation methods.     

Protonated clathrate cages enclosing neutral water molecules: (H+)(H2O)21 and (H+)(H2O)28

This paper describes a systematic study on the clathrate structure of (H+)(H2O)21 using tandem mass spectrometry, vibrational predissociation spectroscopy, Monte Carlo simulations, and density functional theory calculations. We produced (H+)(H2O)n from a continuous corona-discharged supersonic expansion and observed three anomalies simultaneously at the cluster temperature near 150 K, including (1) the peak at n=21 is more intense than its neighboring ions in the mass spectrum, (2) the size-dependent dissociation fractions show a distinct drop for the 21-mer, and (3) the infrared spectrum of (H+)(H2O)21 exhibits only a single feature at 3699 cm(-1), corresponding to the free-OH stretching of three-coordinated water molecules. Interestingly, the anomalies appear or disappear together with cluster temperature, indicating close correlation of these three observations. The observations, together with Monte Carlo simulations and density functional theory calculations, corroborate the notion for the formation of a distorted pentagonal dodecahedral (5(12)) cage with a H2O molecule in the cage and a H3O+ ion on the surface for this "magic number" water cluster ion. The dodecahedral cage melts at higher temperatures, as evidenced by the emergence of a free-OH stretching feature at 3717 cm(-1) for the two-coordinated water in (H+)(H2O)21 produced in a warmer molecular beam. Extension of this study to larger clusters strongly suggests that the experimentally observed isomer of (H+)(H2O)28 is most likely to consist of a distorted protonated pentakaidecahedral (5(12)6(3)) cage enclosing two neutral water molecules.     

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.     

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

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

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.     

Prediction of accurate anharmonic experimental vibrational frequencies for water clusters, (H2O)n, n=2-5

Accurate anharmonic experimental vibrational frequencies for water clusters consisting of 2-5 water molecules have been predicted on the basis of comparing different methods with MP2/aug-cc-pVTZ calculated and experimental anharmonic frequencies. The combination of using HF/6-31G* scaled frequencies for intramolecular modes and anharmonic frequencies for intermolecular modes gives excellent agreement with experiment for the water dimer and trimer and are as good as the expensive anharmonic MP2 calculations. The water trimer, the cyclic Ci and S4 tetramers, and the cyclic pentamer all have unique peaks in the infrared spectrum between 500 and 800 cm-1 and between 3400 and 3700 cm-1. Under the right experimental conditions these different clusters can be uniquely identified using high-resolution IR spectroscopy.     

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.     

Blue-shifted A-H stretching modes and cooperative hydrogen bonding. 1. Complexes of substituted formaldehyde with cyclic hydrogen fluoride and water clusters

The structures and vibrational spectra of the intermolecular complexes formed by insertion of substituted formaldehyde molecules HRCO (R = H, Li, F, Cl) into cyclic hydrogen fluoride and water clusters are studied at the MP2/aug-cc-pVTZ computational level. Depending on the nature of the substituent R, the cluster type, and its size, the C-H stretching modes of HRCO undergo large blue and partly red shifts, whereas all the F-H and O-H stretching modes of the conventional hydrogen bonds are strongly red-shifted. It is shown that (i) the mechanism of blue shifting can be explained within the concept of the negative intramolecular coupling between C-H and C=O bonds that is inherent to the HRCO monomers, (ii) the blue shifts also occur even if no hydrogen bond is formed, and (iii) variation of the acceptor X or the strength of the C-H...X hydrogen bond may either amplify the blue shift or cause a transition from blue shift to red shift. These findings are illustrated by means of intra- and intermolecular scans of the potential energy surfaces. The performance of the negative intramolecular coupling between C-H and C=O bonds of H(2)CO is interpreted in terms of the NBO analysis of the isolated H(2)CO molecule and H(2)CO interacting with (H2O)n and (HF)n clusters.