Size-dependent charge-separation reaction for hydrated sulfate dianion cluster, SO(2-)(H2O)n, with n=3-7

The decrease in the reaction rate for the charge separation in SO(4) (2-)(H(2)O)(n) with increasing cluster size is examined by first-principles calculations of the energetics, activation barriers, and thermal stability for n=3-7. The key factor governing the charge separation is the difference in the strength of solvation interaction: while interaction with water is strong for the reactant SO(4) (2-) and the product OH(-), it is relatively weak for HSO(4) (-). It gives rise to a barrier for charge separation as SO(4) (2-) is transformed into HSO(4) (-) and OH(-), although the overall reaction energy is exothermic. The barrier is high when more than two H(2)O are left to solvate HSO(4) (-), as in the case of symmetric solvation structure and in the case of large clusters. The entropy is another important factor since the potential surface is floppy and the thermal motion facilitates the symmetric distribution of H(2)O around SO(4) (2-), which leads to the gradual reduction in reaction rate and the eventual switch-off of charge separation as cluster size increases. The experimentally observed products for n=3-5 are explained by the thermally most favorable isomer at each size, obtained by ab initio molecular-dynamics simulations rather than by the isomer with the lowest energy.     

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相似文献   

On the stability of IrCl6(3-) and other triply charged anions: solvent stabilization versus ionic fragmentation and electron detachment for the IrCl6(3-).(H2O)n n = 0-10 microsolvated clusters

The intrinsic gas-phase stability of the IrCl(6)(3-) trianion and its microsolvated clusters, IrCl(6)(3-).(H(2)O)(n) n = 1-10, have been investigated using density functional theory (DFT) calculations. Although IrCl(6)(3-) is known to exist as a stable complex ion in bulk solutions, our calculations indicate that the bare trianion is metastable with respect to decay via both electron detachment and ionic fragmentation. To estimate the lifetime of IrCl(6)(3-), we have computed the electron tunneling probability using an adaption of the Wentzel-Kramer-Brillouin theory and predict that the trianion will decay spontaneously via electron tunneling on a time scale of 2.4 x 10(-13) s. The global minimum structure for IrCl(6)(3-).H(2)O was found to contain a bifurcated hydrogen bond, whereas for IrCl(6)(3-).(H(2)O)(2), two low energy minima were identified; one involving two bifurcated water-ion hydrogen bonds and a second combining a bifurcated hydrogen bond with a water-water hydrogen bond. Clusters based on each of these structural motifs were obtained for all of the n = 3-10 systems, and the effect of solvation on the possible decay pathways was explored. The calculations reveal that solvation stabilizes IrCl(6)(3-) with respect to both electron detachment decay and ionic fragmentation, with the magnitude of the repulsive Coulomb barrier for ionic fragmentation increasing smoothly with sequential solvation. This study is the first to compare the propensity for electron detachment versus ionic fragmentation decay for a sequentially solvated triply charged anion.     

Sequential hydration energies of the sulfate ion, from determinations of the equilibrium constants for the gas-phase reactions: SO4(H2O)(n)2- = SO4(H2O)(n-1)2- + H2O

Sequential hydration energies of SO4(H2O)(n)2- were obtained from determinations of the equilibrium constants of the following reactions: SO4(H2O)(n)2- = SO4(H2O)(n-1)2- + H2O. The SO4(2-) ions were produced by electrospray and the equilibrium constants Kn,n-1 were determined with a reaction chamber attached to a mass spectrometer. Determinations of Kn,n-1 at different temperatures were used to obtain DeltaG0n,n-1, DeltaH0 n,n-1, and DeltaS0n,n-1 for n = 7 to 19. Interference of the charge separation reaction SO4(H2O)(n)2- = HSO4(H2O)(n-k)- + OH(H2O)(k-1)- at higher temperatures prevented determinations for n < 7. The DeltaS0n,n-1 values obtained are unusually low and this indicates very loose, disordered structures for the n > or = 7 hydrates. The DeltaH0n,n-1 values are compared with theoretical values DeltaEn,n-1, obtained by Wang, Nicholas, and Wang. Rate constant determinations of the dissociation reactions n,n - 1, obtained with the BIRD method by Wong and Williams, showed relatively lower rates for n = 6 and 12, which indicate that these hydrates are more stable. No discontinuities of the DeltaG0n,n-1 values indicating an unusually stable n = 12 hydrate were observed in the present work. Rate constants evaluated from the DeltaG0n,n-1 results also fail to indicate a lower rate for n = 12. An analysis of the conditions used in the two types of experiments indicates that the different results reflect the different energy distributions expected at the dissociation threshold. Higher internal energies prevail in the equilibrium measurements and allow the participation of more disordered transition states in the reaction.     

Theoretical study of the solvation of HgCl2, HgClOH, Hg(OH)2 and HgCl3(-): a density functional theory cluster approach

The determination of the solvation shell of Hg(II)-containing molecules and especially the interaction between Hg(II) and water molecules is the first requirement to understand the transmembrane passage of Hg into the cell. We report a systematic DFT study by stepwise solvation of HgCl(2) including up to 24 water molecules. In order to include pH and salinity effects, the solvation patterns of HgClOH, Hg(OH)(2) and HgCl(3)(-) were also studied using 24 water molecules. In all cases the hydrogen bond network is crucial to allow orbital-driven interactions between Hg(II) and the water molecules. DFT Born-Oppenheimer molecular dynamics simulations starting from the stable HgCl(2)-(H(2)O)(24) structure revealed that an HgCl(2)-(H(2)O)(3) trigonal bipyramid effective solute appears and then the remaining 21 water molecules build a complete first solvation shell, in the form of a water-clathrate. In the HgCl(2), HgClOH, Hg(OH)(2)-(H(2)O)(24) optimized structures Hg also directly interacts with 3 water molecules from an orbital point of view (three Hg-O donor-acceptor type bonds). All the other interactions are through hydrogen bonding. The cluster-derived solvation energies of HgCl(2), HgClOH and Hg(OH)(2) are estimated to be -34.4, -40.1 and -47.2 kcal mol(-1), respectively.     

Collision induced complex formation following electron capture of SO2-H2O complex interacting with argon atoms

Electron capture dynamics of SO(2)-H(2)O(Ar)(n) complexes (n = 0-2) have been investigated by means of direct ab initio molecular dynamics (MD) method in order to elucidate the effects of solvent argon on the reaction dynamics of SO(2)-H(2)O. The neutral complex of SO(2)-H(2)O has a C(s) symmetry, and the sulfur of SO(2) interacts with the oxygen of H(2)O with an eclipsed form. In the SO(2)-H(2)O(Ar)(n) complexes, the dipole of H(2)O interacts with the argon atoms in the most stable structure. Following the electron capture of the complex SO(2)-H(2)O, the complex anion SO(2)(-)(H(2)O) is dissociated directly into SO(2)(-) + H(2)O. On the other hand, the electron capture of SO(2)(H(2)O)(Ar)(n) argon complex (n = 1-2) leads to the anion-water complex SO(2)(-)(H(2)O) because the collision of H(2)O with the Ar atom causes a rebound of H(2)O from Ar atom to the SO(2)(-) anion. The argon solvent enhanced the SO(2)(-)(H(2)O) complex formation. The reaction mechanism of SO(2)(H(2)O) in the participation of argon atoms was discussed on the basis of the present results.     

The identification of a solvated electron pair in the gaseous clusters of Na(-)(H2O)n and Li(-)(H2O)n

By first principles calculations, we explore the possibility that Na(-)(H(2)O)(n) and Li(-)(H(2)O)(n) clusters, which have been measured previously by photoelectron experiments, could serve as gas-phase molecular models for the solvation of two electrons. Such models would capture the electron-electron interaction in a solution environment, which is missed in the well-known anionic water clusters (H(2)O)(n) (-). Our results show that by n = 10, the two loosely bound s electrons in Li(-)(H(2)O)(n) are indeed detached from lithium, and they could exist in either the singlet (spin-paring) or the triplet (spin-coupling) state. In contrast, the two electrons would prefer to stay on the sodium atom in Na(-)(H(2)O)(n) and on the surface of the cluster. The formation of a solvated electron pair and the variation in solvation structures make these two cluster series interesting subjects for further experimental investigation.     

Organic-inorganic hybrids based on novel bimolecular [Si2W22Cu2O78(H2O)]12- polyoxometalates and the polynuclear complex cations [Cu(ac)(phen)(H2O)]n n+ (n=2, 3)

The reaction of a monosubstituted Keggin polyoxometalate (POM) generated in situ with copper-phenanthroline complexes in excess ammonium or rubidium acetate led to the formation of the hybrid metal organic-inorganic compounds A7[Cu2(ac)2(phen)2(H2O)2][Cu3(ac)3(phen)3(H2O)3][Si2W22Cu2O78(H2O)].approximately 18 H2O (A=NH4+ (1), Rb+ (2); ac=acetate; phen=1,10-phenanthroline). These compounds are constructed from inorganic and metalorganic interpenetrated sublattices containing the novel bimolecular Keggin POM, [Si2W22Cu2O78(H2O)]12-, and Cu-ac-phen complexes, [Cu(ac)(phen)(H2O)]n n+ (n=2, 3). The packing of compound 1 can be viewed as a stacking of open-framework layers parallel to the xy plane built of hydrogen-bonded POMs, and zigzag columns of pi-stacked Cu-ac-phen complex cations running along the [111] direction. Magnetic and EPR results are discussed with respect to the crystal structure of the compounds. DFT calculations on [Cu(ac)(phen)(H2O)]n n+ cationic complexes have been performed, to check the influence of packing in the complex geometry and determine the magnetic exchange pathways.     

Atoms-in-molecules analysis of extended hypervalent five-center, six-electron (5c-6e) C(2)Z(2)O interactions at the 1,8,9-positions of anthraquinone and 9-methoxyanthracene systems

To clarify the nature of five-center, six-electron (5c-6e) C(2)Z(2)O interactions, atoms-in-molecules (AIM) analysis has been applied to an anthraquinone, 1,8-(MeZ)(2)ATQ (1 (Z=Se), 2 (Z=S), and 3 (Z=O)), and a 9-methoxyanthracene system, 9-MeO-1,8-(MeZ)(2)ATC (4 (Z=Se), 5 (Z=S), and 6 (Z=O)), as well as 1-(MeZ)ATQ (7 (Z=Se), 8 (Z=S), and 9 (Z=O)) and 9-MeO-1-(MeZ)ATC (10 (Z=Se), 11 (Z=S), and 12 (Z=O)). The total electronic energy density (H(b)(r(c))) at the bond critical points (BCPs), an appropriate index for weak interactions, has been examined for 5c-6e C(2)Z(2)O and 3c-4e CZO interactions of the n(p)(O)sigma*(Z--C) type in 1-12. Some hydrogen-bonded adducts were also re-examined for convenience of comparison. The total electronic energy densities varied in the following order: OO (3: H(b)(r(c))=0.0028 au)=OO (6: 0.0028 au)>OO (9: 0.0025 au)> or =NNHF (0.0024 au)> or =OO (12: 0.0023 au)>H(2)OHOH (0.0015 au)>SO (8: 0.0013 au)=SO (2: 0.0013 au)> or =SO (11: 0.0012 au)=SO (5: 0.0012 au)>HFHF (0.0008 au)=SeO (10: 0.0008 au)=SeO (4: 0.0008 au)> or =SeO (1: 0.0007 au)> or =SeO (7: 0.0006 au)>HCNHF (-0.0013 au). H(b)(r(c)) values for SO were predicted to be smaller than the hydrogen bond of H(2)OHOH and H(b)(r(c)) values for SeO are very close to or slightly smaller than that for HFHF in both the ATQ and 9-MeOATC systems. In the case of Z=Se and S, H(b)(r(c)) values for 5c-6e C(2)Z(2)O interactions are essentially equal to those for 3c-4e CZO if Z is the same. The results demonstrate that two n(p)(O)sigma*(Z--C) 3c-4e interactions effectively connect through the central n(p)(O) orbital to form the extended hypervalent 5c-6e system of the sigma*(C--Z)n(p)(O)sigma*(Z--C) type for Z=Se and S in both systems. Natural bond orbital (NBO) analysis revealed that n(s)(O) also contributes to some extent. The electron charge densities at the BCPs, NBO analysis, and the total energies calculated for 1-12, together with the structural changes in the PhSe derivatives, support the above discussion.     

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.     

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.     

乙醇-水分子团簇C2H5OH(H2O)n (n=1-9)稳定结构的量子化学研究

采用密度泛函理论B3LYP方法, 在B3LYP/6-311++G(2d,2p)//B3LYP/6-311++G(d,p)基组水平上对乙醇-水分子团簇(C₂H₅OH(H₂O)_n (n=1-9))的各种性质进行研究, 如: 优化的几何构型、结构参数、氢键、结合能、平均氢键强度、自然键轨道(NBO)电荷分布、团簇的生长规律等. 结果表明, 从二维(2-D)环状结构到三维(3-D)笼状结构的过渡出现在n=5的乙醇-水分子团簇中. 此外, 利用团簇结合能的二阶差分、形成能、能隙等性质, 发现在n=6时乙醇-水分子团簇的最低能量结构稳定性较好, 可能为幻数结构. 最后, 为了进一步探讨氢键本质, 将C₂H₅OH(H₂O)_n (n=2-9)最低能量结构的各种性质与纯水分子团簇(H₂O)_n (n=3-10)比较, 结果表明前者与后者中的水分子之间氢键相似.     

Structures of [(CO2)n(H2O)m]- (n=1-4, m=1,2) cluster anions. I. Infrared photodissociation spectroscopy

The infrared photodissociation spectra of [(CO(2))(n)(H(2)O)(m)](-) (n=1-4, m=1, 2) are measured in the 3000-3800 cm(-1) range. The [(CO(2))(n)(H(2)O)(1)](-) spectra are characterized by a sharp band around 3570 cm(-1) except for n=1; [(CO(2))(1)(H(2)O)(1)](-) does not photodissociate in the spectral range studied. The [(CO(2))(n)(H(2)O)(2)](-) (n=1, 2) species have similar spectral features with a broadband at approximately 3340 cm(-1). A drastic change in the spectral features is observed for [(CO(2))(3)(H(2)O)(2)](-), where sharp bands appear at 3224, 3321, 3364, 3438, and 3572 cm(-1). Ab initio calculations are performed at the MP2/6-311++G(**) level to provide structural information such as optimized structures, stabilization energies, and vibrational frequencies of the [(CO(2))(n)(H(2)O)(m)](-) species. Comparison between the experimental and theoretical results reveals rather size- and composition-specific hydration manner in [(CO(2))(n)(H(2)O)(m)](-): (1) the incorporated H(2)O is bonded to either CO(2) (-) or C(2)O(4) (-) through two equivalent OH...O hydrogen bonds to form a ring structure in [(CO(2))(n)(H(2)O)(1)](-); (2) two H(2)O molecules are independently bound to the O atoms of CO(2) (-) in [(CO(2))(n)(H(2)O)(2)](-) (n=1, 2); (3) a cyclic structure composed of CO(2) (-) and two H(2)O molecules is formed in [(CO(2))(3)(H(2)O)(2)](-).     

Dissolution thermochemistry of alkali metal dianion salts (M2X1, M = Li+, Na+, and K+ with X = CO3(2-), SO4(2-), C8H8(2-), and B12H12(2-))

The dissolution Gibbs free energies (ΔG°(diss)) of salts (M(2)X(1)) have been calculated by density functional theory (DFT) with Conductor-like Polarizable Continuum Model (CPCM) solvation modeling. The absolute solvation free energies of the alkali metal cations (ΔG(solv)(M(+))) come from the literature, which coincide well with half reduction potential versus SHE data. For solvation free energies of dianions (ΔG(solv)(X(2-))), four different DFT functionals (B3LYP, PBE, BVP86, and M05-2X) were applied with three different sets of atomic radii (UFF, UAKS, and Pauling). Lattice free energies (ΔG(latt)) of salts were determined by three different approaches: (1) volumetric, (2) a cohesive Gibbs free energy (ΔG(coh)) plus gaseous dissociation free energy (ΔG(gas)), and (3) the Born-Haber cycle. The G4 level of theory, electron propagator theory, and stabilization by dielectric medium were used to calculate the second electron affinity to form the dianions CO(3)(2-) and SO(4)(2-). Only the M05-2X/Pauling combination with the three different methods for estimating ΔG(latt) yields the expected negative dissolution free energies (ΔG°(diss)) of M(2)SO(4). Salts with large dianions like M(2)C(8)H(8) and M(2)B(12)H(12) reveal the limitation of using static radii in the volumetric estimation of lattice energies. The value of ΔE(coh) was very dependent on the DFT functional used.     

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.     

First steps towards dissolution of NaSO4- by water

NaSO(4)(-)(H(2)O)(n) (n = 0-4) clusters have been generated in the gas phase as model systems to simulate the first dissolution steps of sulfate salts in water; photoelectron spectroscopy and theoretical calculations indicate that the first three water molecules strongly interact with both Na(+) and SO(4)(2-), forming a three-water solvation ring to start to pry apart the Na(+)SO(4)(2-) contact ion pair.     

Stability and structure in alpha- and beta-Keggin Heteropolytungstates, [X(n)(+)W12O40)](8-n)(-), X = p-block cation

beta-[SiW(12)O(40)](4)(-) (C(3)(v) symmetry) is sufficiently higher in energy than its alpha-isomer analogue that effectively complete conversion to alpha-[SiW(12)O(40)](4)(-) (T(d)) is observed. By contrast, beta- and alpha-[AlW(12)O(40)](5)(-) (beta- and alpha-1; C(3)(v) and T(d), respectively) are sufficiently close in energy that both isomers are readily seen in (27)Al NMR spectra of equilibrated (alpha-beta) mixtures. Recently published DFT calculations ascribe the stability of beta-1 to an electronic effect of the large, electron-donating [AlO(4)](5)(-) (T(d)) moiety encapsulated within the polarizable, fixed-diameter beta-W(12)O(36) (C(3)(v)) shell. Hence, no unique structural distortion of beta-1 is needed or invoked to explain its unprecedented stability. The results of these DFT calculations are confirmed by detailed comparison of the X-ray crystal structure of beta-1 (beta-Cs(4.5)K(0.5)[Al(III)W(12)O(40)].7.5H(2)O; orthorhombic, space group Pmc2(1), a = 16.0441(10) A, b = 13.2270(8) A, c = 20.5919(13) A, Z = 4 (T = 100(2) K)) with previously reported structures of alpha-1, alpha- and beta-[SiW(12)O(40)](4)(-), and beta(1)-[SiMoW(11)O(40)](4)(-).     

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.     

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.