Theoretical investigation of gas phase ethanol–(water)n (n = 1–5) clusters and comparison with gas phase pure water clusters (water)n (n = 2–6)

Various properties (such as optimal structures, structural parameters, hydrogen bonds, natural bond orbital charge distributions, binding energies, electron densities at hydrogen bond critical points, cooperative effects, and so on) of gas phase ethanol–(water)_n (n = 1–5) clusters with the change in the number of water molecules have been systematically explored at the MP2/aug‐cc‐pVTZ//MP2/6‐311++G(d,p) computational level. The study of optimal structures shows that the most stable ethanol‐water heterodimer is the one where exists one primary hydrogen bond (O? H…O) and one secondary hydrogen bond (C? H …O) simultaneously. The cyclic geometric pattern formed by the primary hydrogen bonds, where all the molecules are proton acceptor and proton donor simultaneously, is the most stable configuration for ethanol–(water)_n (n = 2–4) clusters, and a transition from two‐dimensional cyclic to three‐dimensional structures occurs at n = 5. At the same time, the cluster stability seems to correlate with the number of primary hydrogen bonds, because the secondary hydrogen bond was extremely weaker than the primary hydrogen bond. Furthermore, the comparison of cooperative effects between ethanol–water clusters and gas phase pure water clusters has been analyzed from two aspects. First of all, for the cyclic structure, the cooperative effect in the former is slightly stronger than that of the latter with the increasing of water molecules. Second, for the ethanol–(water)₅ and (water)₆ structure, the cooperative effect in the former is also correspondingly stronger than that of the latter except for the ethanol–(water)₅ book structure. © 2012 Wiley Periodicals, Inc.     

Geometries,vibrational frequencies,and electron affinities of X2Cl (X=C,Si,Ge) clusters

Ab initio quantum chemical calculations have been performed on X₂Cl^? and X₂Cl (X = C, Si, Ge) clusters. The geometrical structures, vibrational frequencies, electronic properties and dissociation energies are investigated at the Hartree–Fock (HF), Møller–Plesset second‐ and fourth‐order (MP2, MP4), CCSD(T) level with the 6‐311+G(d) basis set. The X₂Cl (X = C, Si, Ge) and X₂Cl^? (X = Si, Ge) take a bent shape obtained at the ground state, while C₂Cl^? has a linear structure. The impact on internal electron transfer between the X₂Cl and the corresponding anional clusters is studied. The three different types of electron affinities (EAs) at the CCSD(T) are reported. The most reliable adiabatic electronic affinities, obtained at the CCSD(T)/cc‐pvqz level of theory, are predicted to be 3.30, 2.62, and 1.98 eV for C₂Cl, Si₂Cl, and Ge₂Cl, respectively. The calculated EAs of C₂Cl and Ge₂Cl are in good agreement with theoretical results reported. The correlation effects and basis sets effects on the geometrical structures and dissociation energies are discussed. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

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     

The Evolution of Bonding and Thermodynamic Properties of Boron‐Doped Small Carbon Clusters: An Ab Initio Study

Theoretical studies on BC_n (n=1–6) clusters are carried out using density functional theory, Møller–Plesset second‐order perturbation theory (MP2), coupled‐cluster calculations including up to triple excitations (CCSD(T)), and higher‐level approaches. All possible isomers depending on the positions of the boron atom are generated and the lowest‐energy isomers are determined for doublet and quartet electronic states. The three potential evolution paths of the clusters are determined as a function of their size. The energetic and electronic consequences for the increased size of structures differ significantly, which leads to representatives of the ground electronic state from different structural groups. The ab initio calculated thermal functions allow enhancements to the available atomization energies and improve the agreement between the calculated and experimental heat content.     

Stability,Infrared Spectra and Electronic Structures of (ZrO2)n(n=3-6)Clusters:DFT Study

The stability, infrared spectra and electronic structures of (ZrO₂)_n (n=3–6) clusters have been investigated by using density‐functional theory (DFT) at B3LYP/6‐31G* level. The lowest‐energy structures have been recognized by considering a number of structural isomers for each cluster size. It is found that the lowest‐energy (ZrO₂)₅ cluster is the most stable among the (ZrO₂)_n (n=3–6) clusters. The vibration spectra of Zr? O stretching motion from terminal oxygen atom locate between 900 and 1000 cm^?1, and the vibrational band of Zr? O? Zr? O four member ring is obtained at 600–700 cm^?1, which are in good agreement with the experimental results. Mulliken populations and NBO charges of (ZrO₂)_n clusters indicate that the charge transfers occur between 4d orbital of Zr atoms and 2p orbital of O atoms. HOMO‐LUMO gaps illustrate that chemical stabilities of the lowest‐energy (ZrO₂)_n (n=3–6) clusters display an even‐odd alternating pattern with increasing cluster size.     

Theoretical investigation of the structures and spectroscopic properties of (H2O4)n (n = 1–4) clusters

Density functional theory and ab initio calculations were performed to elucidate the hydrogen interactions in (H₂O₄)_n (n = 1–4) clusters. The optimized geometries, binding energies, and harmonic vibrational frequencies were predicted at various levels of theory. The trans conformer of the H₂O₄ monomer was predicted to be the most stable structure at the CCSD(T)/aug‐cc‐pVTZ level of theory. The binding energies per H₂O₄ monomer increased in absolute value by 9.0, 10.1, and 11.8 kcal/mol from n = 2 to n = 4 at the MP2/cc‐pVTZ level of theory (after the zero‐point vibrational energy and basis set superposition error corrections). This result implies that the intermolecular hydrogen bonds were stronger in the long‐chain clusters, that is, the formation of the longer chain in the (H₂O₄)_n clusters was more energetically favorable.     

Ab initio study of Rg2I− (Rg = Ar,Kr, Xe)

The equilibrium geometries, vibrational frequencies, and dissociation energies of rare gas iodine clusters Rg₂I^?(Rg = Ar, Kr, Xe) were calculated at the Hartree–Fock (HF), second‐order Møller–Plesset (MP2), the coupled cluster method with single and double excitation and a noniterative correction for triple excitations method [CCSD(T)] levels. The title species have bent C_2v structure of about 60° angle. The electron correlation effects and relativistic effects on the geometry and stability were investigated at CCSD(T) level. Both effects stabilize title species. The calculated electron affinities are in good agreement with the experimental results available. The effect of high angular momentum functions (g and h) was studied. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Dynamics and structural changes of small water clusters on ionization

Despite utmost importance in understanding water ionization process, reliable theoretical results of structural changes and molecular dynamics (MD) of water clusters on ionization have hardly been reported yet. Here, we investigate the water cations [(H₂O)_{n = 2–6}⁺] with density functional theory (DFT), Möller–Plesset second‐order perturbation theory (MP2), and coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. The complete basis set limits of interaction energies at the CCSD(T) level are reported, and the geometrical structures, electronic properties, and infrared spectra are investigated. The characteristics of structures and spectra of the water cluster cations reflect the formation of the hydronium cation moiety (H₃O⁺) and the hydroxyl radical. Although most density functionals fail to predict reasonable energetics of the water cations, some functionals are found to be reliable, in reasonable agreement with high‐level ab initio results. To understand the ionization process of water clusters, DFT‐ and MP2‐based Born‐Oppenheimer MD (BOMD) simulations are performed on ionization. On ionization, the water clusters tend to have an Eigen‐like form with the hydronium cation instead of a Zundel‐like form, based on reliable BOMD simulations. For the vertically ionized water hexamer, the relatively stable (H₂O)₅⁺ (5sL4A) cluster tends to form with a detached water molecule (H₂O). © 2013 Wiley Periodicals, Inc.     

Probing the Properties of Polynuclear Superhalogens without Halogen Ligand via ab Initio Calculations: A Case Study on Double‐Bridged [Mg2(CN)5]−1 Anions

An ab initio study of the superhalogen properties of eighteen binuclear double‐bridged [Mg₂(CN)₅]^?1 clusters is reported herein by using various theoretical methods. High‐level CCSD(T) results indicate that all the clusters possess strong superhalogen properties owing to their high vertical electron detachment energies (VDEs), which exceed 6.8 eV (highest: 8.15 eV). The outer valence Green's function method provides inaccurate relative VDE values; hence, this method is not suitable for this kind of polynuclear superhalogens. Both the HF and MP2 results are generally consistent with the CCSD(T) level regarding the relative VDE values and—especially interesting—the average values of the HF and MP2 VDEs are extremely close to the CCSD(T) results. The distributions of the extra electrons of the anions are mainly aggregated into the terminal CN units. These distributions are apparently different from those of previously reported triple‐bridged isomers and may be the reason for the decreased VDE values of the clusters. In addition, comparisons of the VDEs of binuclear and mononuclear superhalogens as well as studies of the thermodynamic stabilities with respect to the detachment of various CN^?1 ligands are also performed. These results confirm that polynuclear structures with pseudohalogen ligands can be considered as probable new superhalogens with enhanced properties.     

Theoretical determination of the differential polarizability and anisotropy of alkaline earth oxide nanoclusters (BeO)n [n = 2–9]: The basis set and electron correlation effects

As an important lineage of metal oxides, alkaline earth oxides have attracted significant interest due to their unique properties and potential applications in material science and industry. In this article, we present the first ab initio (HF, MP2, and CCSD(T)) and density functional theory (TPSSh functional) investigation on the optical properties such as polarizabilities per atom (PPA), differential polarizability per unit (DPPU), and anisotropies of (BeO)_n [n = 2–9] nanoclusters as an illustrative example of alkaline earth oxides nanostructures. Basis set augmentation effects on the studied properties of BeO nanoclusters have been explored by using basis sets of triple‐zeta quality starting from 6‐311G with increasing completeness of the diffuse and polarization functions to the 6‐311+G(3df) basis. Checking carefully the basis set effects, it is shown that the 6‐311+G(3d) basis set provides the best compromise between the accuracy and computational cost. We found a decreasing trend for PPA values of BeO nanoclusters using all considered methods, indicating the strong electron delocalization with increasing cluster size. Moreover, in accordance with the energetic analysis of stability of BeO nanoclusters, the values of PPA show that the (BeO)₄ and (BeO)₆ clusters are the most stable magic numbers compared to the neighbors, satisfying the minimum polarizability principle. The computed values of DPPU demonstrate a strong binding effect in BeO nanoclusters. Taking into account the electron correlation correction (ECC), it is observed that the variations of ECC on dipole polarizabilities are almost smooth for clusters under study. © 2013 Wiley Periodicals, Inc.     

Search for the global minimum structures of AlB3H2n (n = 0 − 6) clusters

The global minimum structures of AlB₃H_2n (n = 0–6) clusters are determined using the stochastic search method at the B3LYP/6–31G level of theory. These initially specified geometries are recalculated using B3LYP and CCSD(T) methods using the 6–311++G^** basis set. The structural and electronic properties of the two lowest‐lying isomers are presented. The structural parameters obtained for aluminum borohydride are compared with the experimental and theoretical results. The H₂ fragmentation energies of the most stable isomers are investigated. Chemical bonding analyses for the global minimum of AlB₃H_2n (n = 0–6) clusters are performed using the adaptive natural density partitioning method. © 2014 Wiley Periodicals, Inc.     

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     

Accurate adsorption energies for small molecules on oxide surfaces: CH4/MgO(001) and C2H6/MgO(001)

A hybrid method is applied that combines second order Møller–Plesset perturbation theory (MP2) for cluster models with density functional theory for periodic (slab) models to obtain structures and energies for methane and ethane molecules adsorbed on the MgO(001) surface. Single point calculations are performed to estimate the effect of increasing the cluster size on the MP2 energies and to evaluate the difference between coupled cluster (CCSD(T)) and MP2 energies. The final estimates of the adsorption energies are 12.9 ± 1.3 and 18.9 ± 1.8 kJ/mol for CH₄ and C₂H₆, respectively. © 2016 Wiley Periodicals, Inc.     

Characterization of the potential energy surfaces of two small but challenging noncovalent dimers: (P2)2 and (PCCP)2

This work characterizes eight stationary points of the P₂ dimer and six stationary points of the PCCP dimer, including a newly identified minimum on both potential energy surfaces. Full geometry optimizations and corresponding harmonic vibrational frequencies were computed with the second‐order Møller–Plesset (MP2) electronic structure method and six different basis sets: aug‐cc‐pVXZ, aug‐cc‐pV(X+d)Z, and aug‐cc‐pCVXZ where X = T, Q. A new L‐shaped structure with C₂ symmetry is the only minimum for the P₂ dimer at the MP2 level of theory with these basis sets. The previously reported parallel‐slipped structure with C_2h symmetry and a newly identified cross configuration with D₂ symmetry are the only minima for the PCCP dimer. Single point energies were also computed using the canonical MP2 and CCSD(T) methods as well as the explicitly correlated MP2‐F12 and CCSD(T)‐F12 methods and the aug‐cc‐pVXZ (X = D, T, Q, 5) basis sets. The energetics obtained with the explicitly correlated methods were very similar to the canonical results for the larger basis sets. Extrapolations were performed to estimate the complete basis set (CBS) limit MP2 and CCSD(T) binding energies. MP2 and MP2‐F12 significantly overbind the P₂ and PCCP dimers relative to the CCSD(T) and CCSD(T)‐F12 binding energies by as much as 1.5 kcal mol^?1 for the former and 5.0 kcal mol^?1 for the latter at the CBS limit. The dominant attractive component of the interaction energy for each dimer configuration was dispersion according to several symmetry‐adapted perturbation theory analyses. © 2014 Wiley Periodicals, Inc.     

HF(H2O)n clusters with an excess electron: ab initio study

The structures of electron-bound and neutral clusters of HF(H2O)n (n=1-3) were optimized at the level of second-order Moller-Plesset perturbation theory (MP2). Then, the energies were studied using the coupled cluster singles, doubles, and perturbative triples correction [CCSD(T)] method. The vertical detachment energies of the electron-bound clusters for n=1-3 are 60, 180, and approximately 300 meV, respectively. In the case of the n=3, two structures are competing energetically. The electron-bound clusters for n=1 and 2 are 1.5 and 1.8 kcal/mol more stable than the neutral, while that for n=3 is 0.6-0.9 kcal/mol less stable. The excess electron is stabilized in the surface-bound state of the dipole oriented structures of the hydrated acid clusters. Vibrational spectra of the electron-bound clusters are discussed.     

硼碳团簇BnC2 (n＝1～6)的理论研究

应用密度泛函理论在B3LYP/6-311＋G(d)水平上研究了硼碳团簇B_nC₂ (n＝1～6)的几何结构、生长机制和相对稳定性. 计算结果表明, 对于n＝2～6的簇, 平面多环状构型为最稳定的结构, 其中C原子分布于环的顶点、有尽可能多的三配位硼原子和尽可能多的B—C键. 碳原子作为杂原子倾向掺杂于团簇的顶点位置, 它的掺杂不改变硼团簇的主体结构. 与平面多环状结构相比, 随着簇尺寸的增大, 三维结构和线性链结构更不稳定. 在低能线性结构中, C原子位于链两侧的第二个位置. 计算的碎片分裂能、递增键能以及HOMO-LUMO能隙表明, B₄C₂为幻数簇.     

Electronic structure and stability of BP clusters: theoretical calculations for (BP)n (n = 2–4)

The geometry, electronic configurations, harmonic vibrational frequencies, and stability of the structural isomers of boron phosphide clusters have been investigated using density functional theory (DFT). CCSD(T) calculations show that the lowest‐energy structures are cyclic (IIt, IVs) with D_nh symmetry for dimers and trimers. The caged structure for B₄P₄ lie higher in energy than the monocyclic structure with D_2d symmetry (VIs). The B–P bond dominates the structures for many isomers, so that one preferred dissociation channel is loss of the BP monomer. The hybridization and chemical bonding in the different structures are also discussed. Comparisons with boron nitride clusters, the ground state structures of B_nP_n (n = 2, 3) clusters are analogous to those of their corresponding B_nN_n (n = 2, 3) counterparts. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Density Functional Theory Study on (Cl2GaN3)n(n=1–4) Clusters

The molecular geometries, vibrational properties, and thermodynamic properties of the clusters (Cl₂GaN₃)_n(n=1–4) have been predicted at the B3LYP/6‐311+G* level. The optimized clusters (Cl₂GaN₃)_n (n=2–4) all possess cyclic structures containing Ga N_α Ga linkages. The relationships between geometrical parameters and oligomerization degree n are discussed. The gas‐phase structures of the trimers prefer to exist in boat‐twisting conformation. As for the tetramer, the S₄ symmetry structure is the most stable. The infrared spectra are obtained and assigned by vibrational analysis. Thermodynamic properties derived from the infrared spectra on the basis of statistical thermodynamic principles are linearly correlated with the oligomerization degree n as well as the temperature. Meanwhile, thermodynamic analysis of the gas‐phase reaction suggests that the oligomerization is exothermic and favorable under high temperature.     

Cooperative effects in polymolecular nitrogen clusters

The structures and stabilities of the van der Waals clusters (N₂) _n (n = 2–8) have been evaluated using ab initio calculations at the MP2(full)/6–311+G* and CCSD(full)/6–311+G* level of theory. At n = 2–4, the formation of planar and three-dimensional structures is possible. At n > 4, only “globular” structures can be formed, whose formation energies increase with the cluster size. As the number of interacting molecules increases, the examined systems exhibit cooperative effects associated with non-additive enhancement of nonbonded intermolecular interactions.     

Theoretical study in [C2H4–Tl]+ and [C2H2–Tl]+ complexes

We studied the attraction between [C₂H_n] and Tl(I) in the hypothetical [C₂H_n–Tl]⁺ complexes (n = 2,4) using ab initio methodology. We found that the changes around the equilibrium distance C–Tl and in the interaction energies are sensitive to the electron correlation potential. We evaluated these effects using several levels of theory, including Hartree–Fock (HF), second‐order Møller–Plesset (MP2), MP4, coupled cluster singles and doubles CCSD(T), and local density approximation augmented by nonlocal corrections for exchange and correlation due to Becke and Perdew (LDA/BP). The obtained interaction energies differences at the equilibrium distance R_e (C–Tl) range from 33 and 46 kJ/mol at the different levels used. These results indicate that the interaction between olefinic systems and Tl(I) are a real minimum on the potential energy surfaces (PES). We can predict that these new complexes are viable for synthesizing. At long distances, the behavior of the [C₂H_n]–Tl⁺ interaction may be related mainly to charge‐induced dipole and dispersion terms, both involving the individual properties of the olefinic π‐system and thallium ion. However, the charge‐induced dipole term (R^?4) is found as the principal contribution in the stability at long and short distances. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007