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.     

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.     

Global potential energy minima of C60(H2O)n clusters

Likely candidates for the global potential energy minima of C60(H2O)n clusters with n < or = 21 are found using basin-hopping global optimization. The potential energy surfaces are constructed using the TIP4P intermolecular potential for the water molecules, a Lennard-Jones water-fullerene potential, and a water-fullerene polarization potential, which depends on the first few nonvanishing C60 multipole polarizabilities. This combination produces a rather hydrophobic water-fullerene interaction. As a consequence, the water component of the lowest C60(H2O)n minima is quite closely related to low-lying minima of the corresponding TIP4P (H2O)n clusters. In most cases, the geometrical substructure of the water molecules in the C60(H2O)n global minimum coincides with that of the corresponding free water cluster. Exceptions occur when the interaction with C60 induces a change in geometry. This qualitative picture does not change significantly if we use the TIP3P model for the water-water interaction. Structures such as C60@(H2O)60, in which the water molecules surround the C60 fullerene, correspond to local minima with much higher potential energies. For such a structure to become the global minimum, the magnitude of the water-fullerene interaction must be increased to an unphysical value.     

Collision-induced dissociation studies of protonated ether--(H2O)n (n = 1-3) clusters

We have studied the protonated ether-(H2O)n (n = 1-3) complexes containing tetrahydrofuran, dimethyl, diethyl, dibutyl, and butylmethyl ethers using a flowing afterglow triple-quadrupole mass spectrometer. Collision-induced dissociation, CID, of all clusters with n = 1, 2 shows sequential water loss. The n = 3 cluster of dimethyl ether shows sequential water loss, while all other ether clusters display selective product formation. The CID spectra are interpreted based on known energetics, and theoretical studies of the dimethyl and diethyl ether systems.     

Infrared spectra of HC2-.(C2H2)n and O2-.(C2H2)n clusters (n=2-5)

We study the solvation of HC2- and O2- with acetylene ligands by means of midinfrared photodissociation spectroscopy in the CH stretching region, monitoring C2H2 evaporation upon infrared photon absorption by the parent cluster ions. Our findings are interpreted with the help of density functional theory. The infrared spectra indicate that while the binding generally occurs through ionic H bonds, there are two different classes of ligands which differ in their binding strength. This holds true for both core ions, even though their electronic structures and charge distributions are very different.     

Electron scattering on OH-(H2O)n clusters (n = 0-4)

The cross sections for electron scattering on OH-(H2O)n for n = 0-4 were measured from threshold to approximately 50 eV. All detachment cross sections were found to follow the classical prediction given earlier [Phys. Rev. Lett. 74, 892 (1995)] with a threshold energy for electron-impact detachment that increased upon sequential hydration, yielding values in the range from 4.5 eV +/- 0.2 eV for OH- to 12.10 eV +/- 0.5 eV for OH-(H2O)4. For n > or = 1, we found that approximately 80% of the total reaction events lead to electron detachment plus total dissociation of the clusters into the constituent molecules of OH and H2O. Finally, we observed resonances in the cross sections for OH-(H2O)3 and for OH-(H2O)4. The resonances were located at approximately 15 eV and were ascribed to the formation of dianions in excited states.     

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

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

Ab initio studies of F−(H2O)n and Cl−(H2O)n clusters for n = 1, 2

Ab initio Hartree–Fock calculations are performed on hydrates of the F^? and Cl^? ions using 6-31G, 6-31G**, and 6-21G basis sets. Geometries and binding energies are obtained. An estimate of the correlation energy is provided by an MP2/6-31G (Møller-Plesset second-order perturbation) calculation. Comparisons are made between the Cl^?(SO₂) and the Cl^?(H₂O) complexes.     

Vibrational spectroscopy of (SO4(2-)).(H2O)n clusters, n=1-5: harmonic and anharmonic calculations and experiment

The vibrational spectroscopy of (SO4(2-)).(H2O)n is studied by theoretical calculations for n=1-5, and the results are compared with experiments for n=3-5. The calculations use both ab initio MP2 and DFT/B3LYP potential energy surfaces. Both harmonic and anharmonic calculations are reported, the latter with the CC-VSCF method. The main findings are the following: (1) With one exception (H2O bending mode), the anharmonicity of the observed transitions, all in the experimental window of 540-1850 cm(-1), is negligible. The computed anharmonic coupling suggests that intramolecular vibrational redistribution does not play any role for the observed linewidths. (2) Comparison with experiment at the harmonic level of computed fundamental frequencies indicates that MP2 is significantly more accurate than DFT/B3LYP for these systems. (3) Strong anharmonic effects are, however, calculated for numerous transitions of these systems, which are outside the present observation window. These include fundamentals as well as combination modes. (4) Combination modes for the n=1 and n=2 clusters are computed. Several relatively strong combination transitions are predicted. These show strong anharmonic effects. (5) An interesting effect of the zero point energy (ZPE) on structure is found for (SO4(2-)).(H2O)(5): The global minimum of the potential energy corresponds to a C(s) structure, but with incorporation of ZPE the lowest energy structure is C2v, in accordance with experiment. (6) No stable structures were found for (OH-).(HSO4-).(H2O)n, for n相似文献   

Photoionization and ab initio study of Ba(H2O)n (n = 1-4) clusters

An experimental and theoretical study of the photoionization energies (IE's) of Ba(H(2)O)(n) clusters containing up to n = 4 water molecules has been performed. The clusters were generated by a pick-up source combining laser vaporization with pulsed supersonic expansion, and then photoionized by radiation of 272.5-340 nm. The experimentally determined IE(e)'s for n = 1 to 4 are 4.56 ± 0.05, 4.26 ± 0.05, 3.90 ± 0.05 and 3.71 ± 0.05 eV. This cluster size dependence of IE is reproduced within ±0.06 eV employing the mPW1PW91 density-functional and CCSD(T, Full) quantum-chemical methods combined with the 6-311++G(d,p) basis set for the H and O atoms and three different relativistic effective core potentials for Ba atoms. The calculations indicate that the lowest energy hydration structures represent the most relevant contributions to both the vertical and adiabatic ionization energies. Experimental and theoretical evidence correlates with the progressive surface-delocalization of the electron from the hydration cavity around the Ba atom and suggests that the intra-cluster electron transfer is possible even for small aggregates.     

Density functional study of HO(H2O)n (n = 1–3) clusters

The hydrogen bonding complexes HO(H₂O)_n (n = 1–3) were completely investigated in the present study using DFT and MP2 methods at varied basis set levels from 6‐31++G(d,p) to 6‐311++G(2d,2p). For n = 1 two, for n = 2 two, and for n = 3 five reasonable geometries are considered. The optimized geometric parameters and interaction energies for various complexes at different levels are estimated. The infrared spectrum frequencies and IR intensities of the most stable structures are reported. Finally, thermochemistry studies are also carried out. The results indicate that the formation and the number of hydrogen bonding have played an important role in the structures and relative stabilities of different complexes. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Ab initio determination of the ionization potentials of water clusters (H2O)n (n = 2-6)

High-level quantum-chemical ab initio coupled-cluster and multiconfigurational perturbation methods have been used to compute the vertical and adiabatic ionization potentials of several water clusters: dimer, trimer, tetramer, pentamer, hexamer book, hexamer ring, hexamer cage, and hexamer prism. The present results establish reference values at a level not reported before for these systems, calibrating different computational strategies and helping to discard less reliable theoretical and experimental data. The systematic study with the increasing size of the water cluster allows obtaining some clues on the structure and reductive properties of liquid water.     

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

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

Global optimization of ionic Mg(n)F(2n) (n=1-30) clusters

The global optimization basin-hopping (BH) method has been used to locate the global minima (GM) of Mg(n)F(2n) (n=1-30) clusters using a Born-Mayer-type potential. Some of the GM were particularly difficult to find, requiring more than 1.5 x 10(4) BH steps. We have found that both the binding energy per MgF2 unit and the effective volume of the GM isomers increase almost linearly with n, and that cluster symmetry decreases with cluster size. The data derived from the BH runs reveal a growing density of local minima just above the GM as n increases. Despite this, the attraction basin around each GM is relatively large, since after all their atomic coordinates are randomly displaced by values as high as 2.0 bohrs, the perturbed structures, upon reoptimization, relax back to the GM in more than 50% of the cases (except for n=10 and 11). The relative stabilities derived from energy second differences suggest that n=8,10,13,15, and 20 are probably the magic numbers for these systems. Mass spectrum experiments would be very useful to clarify this issue.     

Stoichiometric and oxygen-rich M2O(n)- and M2O(n) (M = Nb, Ta; n = 5-7) clusters: molecular models for oxygen radicals, diradicals, and superoxides

We investigated the structures and bonding of two series of early transition-metal oxide clusters, M(2)O(n)(-) and M(2)O(n) (M = Nb, Ta; n = 5-7) using photoelectron spectroscopy (PES) and density-functional theory (DFT). The stoichiometric M(2)O(5) clusters are found to be closed shell with large HOMO-LUMO gaps, and their electron affinities (EAs) are measured to be 3.33 and 3.71 eV for M = Nb and Ta, respectively; whereas EAs for the oxygen-rich clusters are found to be much higher: 5.35, 5.25, 5.28, and 5.15 eV for Nb(2)O(6), Nb(2)O(7), Ta(2)O(6), and Ta(2)O(7), respectively. Structural searches at the B3LYP level yield triplet and doublet ground states for the oxygen-rich neutral and anionic clusters, respectively. Spin density analyses reveal oxygen radical, diradical, and superoxide characters in the oxygen-rich clusters. The M(2)O(7)(-) and M(2)O(7) clusters, which can be viewed to be formed by M(2)O(5)(-/0) + O(2), are utilized as molecular models to understand dioxygen activation on M(2)O(5)(-) and M(2)O(5) clusters. The O(2) adsorption energies on the stoichiometric M(2)O(5) neutrals are shown to be surprisingly high (1.3-1.9 eV), suggesting strong capabilities to activate O(2) by structural defects in Nb and Ta oxides. The PES data also provides valuable benchmarks for various density functionals (B3LYP, BP86, and PW91) for the Nb and Ta oxides.     

Structure and stability of (TiO2)n, (SiO2)n, and mixed Ti(m)Si(n-m)O(2n) [n = 2-5, m = 1 to (n - 1)] clusters

Structures, energetics, and vibrational spectra are investigated for small pure (TiO(2))(n), (SiO(2))(n), and mixed Ti(m)Si(n-m)O(2n) [n = 2-5, m = 1 to (n - 1)] oxide clusters by density functional theory (DFT). The BP86/ATZP level of theory is employed to obtain constitutional isomers of the oxide clusters. In accordance with previous studies, our calculations show three-dimensional compact structures are preferred for pure (TiO(2))(n) with oxo-stabilized higher hexavalent states, and linear chain structures are favored for pure (SiO(2))(n) with tetravalent states. However, the herein theoretically first reported mixed Ti(m)Si(n-m)O(2n) oxide clusters prefer either three-dimensional compact or linear chain structures depending upon the stoichiometry of the compound. Vibrational analysis of the important modes of some highly stable structures is provided. Coupled-cluster single and double excitation (with triples) [CCSD(T)] computed energy gaps for the TiO(2) dimers compare well with results from previous study. Excitation energies are computed by use of time-dependent (TD) DFT and equation-of-motion coupled-cluster calculations with singles and doubles (EOM-CCSD) for the most stable isomers.     

Theoretical study of Al(n) and Al(n)O (n = 2-10) clusters

The stable structures, energies, and electronic properties of neutral, cationic, and anionic clusters of Al(n) (n = 2-10) are studied systematically at the B3LYP/6-311G(2d) level. We find that our optimized structures of Al5(+), Al9(+), Al9(-), Al10, Al10(+), and Al10(-) clusters are more stable than the corresponding ones proposed in previous literature reports. For the studied neutral aluminum clusters, our results show that the stability has an odd/even alternation phenomenon. We also find that the Al3, Al7, Al7(+), and Al7(-) structures are more stable than their neighbors according to their binding energies. For Al7(+) with a special stability, the nucleus-independent chemical shifts and resonance energies are calculated to evaluate its aromaticity. In addition, we present results on hardness, ionization potential, and electron detachment energy. On the basis of the stable structures of the neutral Al(n) (n = 2-10) clusters, the Al(n)O (n = 2-10) clusters are further investigated at the B3LYP/6-311G(2d), and the lowest-energy structures are searched. The structures show that oxygen tends to either be absorbed at the surface of the aluminum clusters or be inserted between Al atoms to form an Al(n-1)OAl motif, of which the Al(n-1) part retains the stable structure of pure aluminum clusters.     

Association patterns in (HF)(m)(H2O)(n) (m + n = 2-8) clusters

In an attempt to understand the phase behavior of aqueous hydrogen fluoride, the clustering in the mixture is investigated at the molecular level. The study is performed at the mPW1B95/6-31+G(d,p) level of theory. Several previous studies attempted to describe the dissociation of HF in water, but in this investigation, the focus is only on the association patterns that are present in this binary mixture. A total of 214 optimized geometries of (HF)n(H2O)m clusters, with m + n as high as 8, were investigated. For each cluster combination, several different conformations are investigated, and the preferred conformations are presented. Using multiple linear regressions, the average strengths of the four possible H-bonding interactions are obtained. The strongest H-bond interaction is reported to be the H2O...H-F interaction. The most probable distributions of mixed clusters as a function of composition are also deduced. It is found that the larger (HF)n(H2O)m clusters are favored both energetically and entropically compared to the ones that are of size m + n < or = 3. Also, the clusters with equimolar contributions of HF and H2O are found to have the strongest interactions.     

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.     

Computational study on the negative electron affinities of NO2 -.(H2O)n clusters (n=0-30)

Here we report negative electron affinities of NO(2)(-).(H2O)n clusters (n=0-30) obtained from density functional theory calculations and a simple correction to Koopmans' theorem. The method relies on the calculation of the detachment energy of the monoanion and its highest occupied molecular orbital and lowest unoccupied molecular orbital energies, and explicit calculations on the dianion itself are avoided. A good agreement with resonances in the cross section for neutral production in electron scattering experiments is found for n=0, 1, and 2. We find several isomeric structures of NO(2)(-).(H2O)2 of similar energy that elucidate the interplay between water-water and ion-water interactions. The topology is predicted to influence the electron affinity by 0.5 and 0.4 eV for NO(2)(-).(H2O) and NO(2)(-).(H2O)2, respectively. The electron affinity of larger clusters is shown to follow a (n+delta)-1/3 dependence, where delta=3 represents the number of water molecules that in volume, could replace NO(2) (-).