Ab initio investigation of water clusters (H2O)n (n = 2–34)

Various properties of typical structures of water clusters in the n = 2–34 size regime with the change of cluster size have been systematically explored. Full optimizations are carried out for the structures presented in this article at the Hartree–Fock (HF) level using the 6‐31G(d) basis set by taking into account the positions of all atoms within the cluster. The influence of the HF level on the results has been reflected by the comparison between the binding energies of (H₂O)_n (n = 2–6, 8, 11, 13, 20) calculated at the HF level and those obtained from high‐level ab initio calculations at the second‐order Møller–Plesset (MP2) perturbation theory and the coupled cluster method including singles and doubles with perturbative triples (CCSD(T)) levels. HF is inaccurate when compared with MP2 and CCSD(T), but it is more practical and allows us to study larger systems. The computed properties characterizing water clusters (H₂O)_n (n = 2–34) include optimal structures, structural parameters, binding energies, hydrogen bonds, charge distributions, dipole moments, and so on. When the cluster size increases, trends of the above various properties have been presented to provide important reference for understanding and describing the nature of the hydrogen bond. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Theoretical investigation of geometry and stability of small lithium‐iodide LinI (n = 2–6) clusters

We present theoretical investigation of the structural characteristics and stabilities of neutral and positively charged Li_nI (n = 2‐6) species. The structural isomers were found by using a randomized algorithm to search for minima structures, followed by B3LYP optimizations; the single‐point RCCSD(T)/cc‐pwCVTZ(‐PP) calculations were performed in order to compute relative energies, binding energies per atom, adiabatic and vertical ionization energies, and dissociation energies. Stability was compared to the pure lithium clusters; there is a typical odd‐even alternation; iodine doped clusters are more stable than pure lithium clusters. Lithium “cage” transfers its valence electron to the iodine atom to form neutral and cationic clusters. An electron departures the lithium cage upon ionization. An important reason for the larger stability of closed‐shell species is the existence of the HOMO 3c/2e natural bond orbitals. © 2013 Wiley Periodicals, Inc.     

Theoretical study of silicon–sulfur clusters (SiS2)n− (n=1–6)

The possible geometrical structures and relative stability of silicon–sulfur clusters (SiS₂) (n=1–6) are explored by means of density functional theory (DFT) quantum chemical calculations. We also compare DFT with second‐order Møller–Plesset (MP2) and Hartree–Fock (HF) methods. The effects of polarization functions, diffuse functions, and electron correlation are included in MP2 and B3LYP quantum chemical calculations, and B3LYP is effective in larger cluster structure optimization, so we can conclude that the DFT approach is useful in establishing trends. The electronic structures and vibrational spectra of the most stable geometrical structures of (SiS₂)_n⁻ are analyzed by B3LYP. As a result, the regularity of the (SiS₂)_n⁻ cluster growing is obtained, and the calculation may predict the formation mechanism of the (SiS₂)_n⁻ cluster. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 280–290, 2001     

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     

Theoretical prediction for the structures of gas phase lithium oxide clusters: (Li2O)n (n = 1–8)

We apply genetic algorithm combining directly with density functional method to search the potential energy surface of lithium‐oxide clusters (Li₂O)_n up to n = 8. In (Li₂O)_n (n = 1–8) clusters, the planar structures are found to be global minimum up to n = 2, and the global minimum structures are all three‐dimensional at n ≥ 3. At n ≥ 4, the tetrahedral unit (TU) is found in most of the stable structures. In the TU, the central Li is bonded with four O atoms in sp³ interactions, which leads to unusual charge transformation, and the probability of the central Li participating in the bonding is higher by adaptive natural density partitioning analysis, so the central Li is in particularly low positive charge. At large cluster size, distortion of structures is viewed, which breaks the symmetry and may make energy higher. The global minimum structures of (Li₂O)₂, (Li₂O)₆, and (Li₂O)₇ clusters are the most stable magic numbers, where the first one is planar and the later both have stable structural units of tetrahedral and C_4v. © 2012 Wiley Periodicals, Inc.     

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.     

DFT studies on TinC2n (n=1–6) clusters

The geometric configurations and electronic structures of the Ti_nC_2n (n=1–6) clusters were studied by using the quantum chemical ab initio density functional theory (DFT) method. Our studies showed that these Ti_nC_2n (n=1–6) could grow gradually to form cyclic clusters through the subunits TiC₂ bonding to each other by C C or Ti C bond. The result could explain the existing experimental fact. The studies might also be helpful to the knowledge of the formation mechanism of the Met‐Cars. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 71: 313–318, 1999     

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     

Cooperativity between the Dihydrogen Bond and the N⋅⋅⋅HC Hydrogen Bond in LiH–(HCN)n Complexes

The cooperativity between the dihydrogen bond and the N???HC hydrogen bond in LiH–(HCN)_n (n=2 and 3) complexes is investigated at the MP2 level of theory. The bond lengths, dipole moments, and energies are analyzed. It is demonstrated that synergetic effects are present in the complexes. The cooperativity contribution of the dihydrogen bond is smaller than that of the N???HC hydrogen bond. The three‐body energy in systems involving different types of hydrogen bonds is larger than that in the same hydrogen‐bonded systems. NBO analyses indicate that orbital interaction, charge transfer, and bond polarization are mainly responsible for the cooperativity between the two types of hydrogen bonds.     

Theoretical study on the interaction of neutral and charged Tin (n = 1–7) clusters with one oxygen molecule

The interactions between the neutral and charged (?2, ?1, +1, and +2) Ti_n (n = 1–7) clusters and one O₂ molecule were investigated by density functional theory. The calculated results show that the oxygen molecule is dissociative on the neutral Ti_n clusters. Geometrically, the two O atoms are distributed at the two sides across the neutral Ti_n cluster for n = 1–4 and the oxygen atom favors the three‐fold hollow site for n = 5, 6, and 7. The binding energy per atom (E_b) and energy gap (E_gap) show higher stability and lower chemical activity of the neutral Ti_nO₂ (n = 1–7) systems compared with the corresponding Ti_n clusters. The adsorption energies (E_ad) exhibit a continuously ascending tendency except for n = 4. The results of the addition of different charges (?2, ?1, +1, and +2) on the most stable neutral Ti_nO₂ (n = 1–7) systems indicate that their geometries are usually perturbed. The stabilities of the neutral Ti_nO₂ systems are enhanced by adding one negative charge. The strongest interaction of the charged Ti_n clusters (?2, ?1, +1, and +2) with O₂ molecule is found at charge +2. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Pseudopolyrotaxanes of Cucurbit[6]uril: A Three‐Dimensional Network Self‐assembled by ClO4−(H2O)2 Water Clusters

A pseudorotaxane of cucurbit[6]uril (CB[6]) with guest molecule N,N′‐hexamethylenebis (pyrazinyl perchlorate) (BPHP) was synthesized and characterized by ¹H NMR spectra, IR, single crystal X‐ray diffraction analysis and thermogravimetric analysis. The structure of the pseudorotaxane (CB[6]·BPHP) is stabilized by host‐guest hydrogen bonds. Self‐assembly of the pseudorotaxane produces infinite one‐dimensional and two‐dimensional networks with intermolecular hydrogen bonds. In the molecular packing of the CB[6]·BPHP, ClO₄^?(H₂O)₂ water clusters serve as bridges to associate these pseudorotaxanes and form three‐dimensional networked pseudopolyrotaxane.     

Darstellung gemischter Kaliumhexabromochloro-osmate(IV), K2OsBrnCl6–n (n = 1…5) und Kaliumhexabromochloro-rhenate(IV), K2ReBrnCl6–n (n = 1…5)

Preparation of Mixed Potassium Hexabromochloro Osmates(IV), K₂OsBr_nCl_6–n (n = 1…5), and Potassium Hexabromochloro Rhenates(IV), K₂ReBr_nCl_6–n (n = 1…5) Solutions of mixed hexabromochloro osmates and rhenates which were separated by ion column chromatography are extracted by tridodecylamine in ether. The potassium salts are precipitated from the extractions with potassium acetate and following further purification are characterized spectrophotometrically and by powder X-ray methods.     

Theoretical study of structure and vibrational properties of MgnOn (n = 3–10) clusters

The structure and harmonic vibrations of Mg_nO_n (n = 3–10) clusters have been investigated using density functional theory. All structures are found to be cumulenic D_nh rings (equal bonds, alternating angles), with one intense out‐of‐plane mode and three infrared (IR)‐active degenerate modes, of which the highest one is extremely intense and increases asymptotically to 1000 cm^?1 for n = 10 at the B3LYP/6‐311++G(2d,2p) level. Comparisons with C_2n clusters show that B_nN_n and Be_nO_n clusters, the structure and bonding type for the Mg_nO_n clusters are consistent with those of the C_2n (n = 3, 5, 7,…) clusters B_nN_n(n = 3–10) and Be_nO_n(n = 3–10) clusters. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Structure and energetic of Bn (n = 2–12) clusters: Electronic structure calculations

The electronic and geometric structures, total and binding energies, first and second energy differences, harmonic frequencies, point symmetries, and highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO–LUMO) gaps of small and neutral B_n (n = 2–12) clusters have been investigated using density functional theory (DFT), B3LYP with 6‐311++G(d,p) basis set. Linear, planar, convex, quasi‐planar, three‐dimensional (3D) cage, and open‐cage structures have been found. None of the lowest energy structures and their isomers has an inner atom; i.e., all the atoms are positioned at the surface. Within this size range, the planar and quasi‐planar (convex) structures have the lowest energies. The first and the second energy differences are used to obtain the most stable sizes. A simple growth path is also discussed with the studied sizes and isomers. The results have been compared with previously available theoretical and experimental works. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Theoretical study on AlnO2 (n = 1–10) clusters and O2 adsorption on the Al(111) surface

The structural properties of neutral and ionic Al_nO₂ (n = 1–10) clusters have been systematically investigated using the density functional method B3LYP with a standard 6‐311+G(d) basis set. The calculated results show that in the Al_nO, Al_nO₂, and Al_nO (n ≥ 3) clusters, O atoms tend to penetrate into the aluminum clusters with some Al atoms moving outward. The binding energies and natural charges populations indicate that the oxygen‐etching is generally stronger in the order Al < Al_n < Al for n < 3, and Al > Al_n > Al for n ≥ 3. To further understand the mechanism of interaction between Al and O₂, the adsorption of O₂ on the Al(111) surface was studied using the density functional theory with plane wave pseudopotential method. The calculated results are consistent with the experimental observation that the O₂ molecule would dissociate on the Al(111) surface and be adsorbed in adjacent hollow sites, forming a local structure of Al₃O–Al₃O. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Structure and stability of coinage metal fluoride and chloride clusters (MnFn and MnCln,M = Cu,Ag, or Au; n = 1–6)

Calculations in the framework of the density functional theory are performed to study the lowest‐energy isomers of coinage metal fluoride and chloride clusters (M_nF_n, M_nCl_n, M = Cu, Ag, or Au, n = 1–6). For all calculated species starting from the trimers the most stable structures are found to be cyclic arrangements. However, planar rings are favored in the case of metal fluorides whereas metal chlorides prefer nonplanar cycles. Calculated bond lengths and infrared frequencies are compared with the available experimental data. The nature of the bonding, involving both covalent and ionic contributions, is characterized. The stability and the fragmentation are also investigated. Trimers are found to be particularly stable when considering the Gibbs free energies. © 2012 Wiley Periodicals, Inc.