Ab initio studies of negative ion-molecule(s) clusters present in the atmosphere. II. OH−(H2O)n for n = 0, 2, 3, 4

Ab initio calculations are performed with 6–31G basis set to study the geometry and binding of the H₃O, H₅O, H₇O, and H₉O complexes. The H₃O complex is also investigated with the 6–31 G* basis set and MP 2 (Moller–Plesset perturbation theory of second order).     

Ab initio studies on the reaction of O2 with Ba(n) (n=2,5) clusters

Ab initio theoretical calculations have been performed to study the reaction of O(2) with Ba(n) (n=2,5) clusters. Our results show that O(2) can easily chemisorb and dissociate on small Ba(n) clusters and there is no obvious energy barrier in the process of the dissociation. The local magnetic moment contributed by oxygen must vanish during the intermediate states before the O(2) dissociation. Correspondingly, local magnetic moment only decreases from 2 mu(B) to about 1 mu(B) if O(2) molecularly adsorbs onto Ba(5) cluster. The electronic structure analysis indicates that the charge transfer from Ba(n) cluster to O(2) as well as the orbital hybridization between the cluster and the oxygen molecule may play a key role in O(2) dissociation.     

Ab initio studies of negative ion-molecule(s) clusters present in the atmosphere. III. OH−(CO2)n for n = 1, 2

Clusters formed by the OH^? ion and carbon dioxide are investigated using ab initio Hartree–Fock calculations, with 6–31 G and 6–31 G* basis sets. Geometries and binding energies are determined.     

Theoretical studies of O2−:(H2O)n clusters

The interaction of superoxide ion O₂^? with up to four water molecules [O₂^?: (H₂O)_n, n = 1, 2, 4] has been investigated using ab initio molecular orbital theory. The binding energy of O₂^?: H₂O is calculated to be ?20.6 kcal/mol in good agreement with gas phase experimental data. At the MP3/6-31G^* level the O₂^?:H₂O complex has a C_2v structure with a double (cyclic) hydrogen bond between O₂^? and H₂O. A C_s structure with a single hydrogen bond is only 0.7 kcal/mol less stable. Interaction of H₂O with the doubly occupied π^* orbital of O₂^? is preferred slightly over interaction with the singly occupied π^* orbital. Natural bond orbital analysis suggests that both electrostatic and charge transfer interactions are important in anionic complexes. The charge transfer occurs predominantly in the O₂^? → H₂O direction and is important in determining the relative stabilities of the different structures and states. Singly and doubly hydrogen-bonded structures for the O₂^?: (H₂O)₂ and O₂^?: (H₂O)₄ clusters were found to be similar in stability and the increase in binding of the cluster becomes smaller as each additional water molecule is added to the cluster.     

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.     

Structures, energetics, and vibrational spectra of H2O2...(H2O)n, n = 1-6 clusters: Ab initio quantum chemical investigations

Hydrogen-bonded heteroclusters of H(2)O(2)...(H(2)O)(n)(), with n varying from 1 through 6, have been investigated herein employing ab initio quantum chemical methods. For a given n, several energetically comparable conformers emerge as local minima on the potential energy surface (PES). All of the conformers obtained at restricted Hartree-Fock (RHF) and M?ller-Plesset second-order perturbation (MP2) levels of theory exhibit parallel trends in energy hierarchy. The effect of clustering by water on the modification in the vibrational frequencies has also been investigated and further, a many-body interaction-energy analysis is carried out providing insights into cooperativity in H(2)O(2)...(H(2)O)(n)() clusters.     

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.     

Ab initio studies of H2PXYH molecules (X,Y = O,S)

Ab initio molecular orbital theory is applied to the study of P? O and P? S bonding in the hypervalent phosphinic (H₂POOH), phosphinothioic (H₂POSH), and phosphinodithioic (H₂PSSH) acid molecules. Intramolecular proton exchange reactions are followed using the intrinsic reaction coordinate and Self-Consistent-Field energy localized orbitals. The P? O and PS bonds are characterized via force constants, phosphorus d orbital populations, and localized orbitals and are best described as either normal single bonds or dative bonds augmented by π back donation from the oxygen or sulfur lone pairs. The anions of these acids are also investigated, and they are found to contain only dative bonds to sulfur and oxygen.     

Ab initio and empirical energy landscapes of (MgF2)n clusters (n = 3, 4)

We explore the energy landscape of (MgF(2))(3) on both the empirical and ab initio level using the threshold algorithm. In order to determine the energy landscape and the dynamics of the trimer we investigate not only the stable isomers but also the barriers separating these isomers. Furthermore, we study the probability flows in order to estimate the stability of all the isomers found. We find that there is reasonable qualitative agreement between the ab initio and empirical potential, and important features such as sub-basins and energetic barriers follow similar trends. However, we observe that the energies are systematically different for the less compact clusters, when comparing empirical and ab initio energies. Since the underlying motivation of this work is to identify the possible clusters present in the gas phase during a low-temperature atom beam deposition synthesis of MgF(2), we employ the same procedure to additionally investigate the energy landscape of the tetramer. For this case, however, we use only the empirical potential.     

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.     

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.     

Ab initio study of electronic structures of Ptn clusters (n = 2–12)

Electronic structures of a series of Pt_n clusters (n = 2–12) have been calculated by ab initio method (SCF/MP2). The result shows that compared to that at the saturated coordination site, Pt at the unsaturated coordination site has lower Pt₆₅ electron occupancies, and the relationship between the electronic structure and the cluster size has been discussed. It is found that when the number of Pt in cluster reaches 7, electronic structures become more like the metallic ones. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 62 : 427–436, 1997     

Theoretical study of [Li(H2O)n]+ and [K(H2O)n]+ (n = 1−4) complexes

The geometries, successive binding energies, vibrational frequencies, and infrared intensities are calculated for the [Li(H₂O)_n]⁺ and [K(H₂O)_n]⁺ (n = 1?4) complexes. The basis sets used are 6-31G* and LANL 1DZ (Los Alamos ECP +DZ ) at the SCF and MP 2 levels. There is an agreement for calculated structures and frequencies between the MP 2/6-31G* and MP 2/LANL 1DZ basis sets, which indicates that the latter can be used for calculations of water complexes with heavier ions. Our results are in a reasonable agreement with available experimental data and facilitate experimental study of these complexes. © 1995 John Wiley & Sons, Inc.     

Ab initio study of neutral and charged SinNap(+) (n <or= 6, p <or= 2) clusters

Ab initio calculations in the framework of the density functional theory, with B3LYP functional, are performed to study the lowest-energy isomers of silicon sodium clusters Si(n)Na(p)(+) (n 相似文献   

Ab initio investigation of vibrational spectra of water-(CO2)n complexes (n = 1, 2)

In this paper, we have calculated using the ab initio method the IR vibrational spectra of complexes of CO2 formed with water (sp3 O-donating atom). Binding energies and structures of the CO2-H2O and water-(CO2)2 complexes have been determined at the second-order level of the Moller-Plesset perturbation theory (MP2) using Dunning's basis sets. The results are presented and critically discussed in terms of the nature of the water-CO2 interactions, electron donor acceptor (EDA) and weak O...H-O interactions. For water-(CO2)2 trimer, it is also shown that the contribution to the interaction energy of the irreducible three-bodies remains relatively negligible. We have analyzed the frequency shifts and the IR and Raman intensity variations under the complex formation. We have particularly emphasized the splitting of the 2 bending mode of CO2 and stretching modes of water, which have been revealed as the most pertinent probes to assess the nature of the forces involved in the different complexes. Finally, because water can play the role of Lewis base and acid as well, we found that weak O...H-O interactions can cooperate with EDA interactions in trimer, leading to very specific spectral signatures that are further discussed.     

Ab initio study of cyclic siloxanes (H2SiO)n: n = 3, 4, 5

The geometry optimizations for several conformations of tri-, tetra-, and pentacyclosiloxane (H₂SiO)_n (n = 3, 4, and 5) were carried out, and the relative stabilities were compared at the Hartree-Fock (HF) and second order perturbation theory (MP2) levels of theory using the 6–31G* and 6–311G(d, p) basis sets. At the highest levels of theory, the only minimum for n = 4 (D₄) occurs at the highly symmetric D_4h structure. In contrast, several, nearly isoenergetic, minima are found on the D₅ surface. These have C₁, C₂, C_s, and D_5h symmetries. While the C₁ structure has the lowest MP2/6–311G(d, p) energy, this species is predicted to be highly fluxional, and the distribution of isomers is dependent on temperature. © 1996 by John Wiley & Sons, Inc.     

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.     

An ab initio study of water clusters in gas phase and bulk aqueous media: (H2O)n,n=1–12

Geometries of several clusters of water molecules including single minimum energy structures of n‐mers (n=1–5), several hexamers and two structures of each of heptamer to decamer derived from hexamer cage and hexamer prism were optimized. One structural form of each of 11‐mer and 12‐mer were also studied. The geometry optimization calculations were performed at the RHF/6‐311G* level for all the cases and at the MP2/6‐311++G** level for some selected cases. The optimized cluster geometries were used to calculate total energies of the clusters in gas phase employing the B3LYP density functional method and the 6‐311G* basis set. Frequency analysis was carried out in all the cases to ensure that the optimized geometries corresponded to total energy minima. Zero‐point and thermal free energy corrections were applied for comparison of energies of certain hexamers. The optimized cluster geometries were used to solvate the clusters in bulk water using the polarized continuum model (PCM) of the self‐consistent reaction field (SCRF) theory, the 6‐311G* basis set, and the B3LYP density functional method. For the cases for which MP2/6‐311++G** geometry optimization was performed, solvation calculations in water were also carried out using the B3LYP density functional method, the 6‐311++G** basis set, and the PCM model of SCRF theory, besides the corresponding gas‐phase calculations. It is found that the cage form of water hexamer cluster is most stable in gas phase among the different hexamers, which is in agreement with the earlier theoretical and experimental results. Further, use of a newly defined relative population index (RPI) in terms of successive total energy differences per water molecule for different cluster sizes suggests that stabilities of trimers, hexamers, and nonamers in gas phase and those of hexamers and nonamers in bulk water would be favored while those of pentamer and decamer in both the phases would be relatively disfavored. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 90–104, 2001     

Ab initio studies on (H2O)14 - clusters: existence of surface and interior-bound extra electrons

A recent paper by Turi et al. [Science 309, 914 (2005)] suggests that the anionic water clusters smaller than (H2O)(45) (-) (at a low temperature) will only have surface-bound extra electrons and no internally bound electrons. Accordingly, (H2O)(14) (-) cluster isomers should only have surface-bound extra electrons. The ab initio results presented here, however, suggest that the (H2O)(14) (-) cluster isomers can have two distinct types of isomers with almost the same energy. The one type of isomer (type 1) has all the non-H-bonding H atoms (NHB H) directed outward and surface-bound extra electron while the other type (type 2) has a number of NHB H atoms directed toward cavity and has an interior-bound electron, and thus, contradicts the earlier quantum simulation results of Turi et al.